Evaluation of various vertical mixing parameterizations in a tropical Pacific Ocean GCM

Coupled general circulation models (GCMs) have had weak El Niٌo/Southern Oscillation variability that has been attributed to a diffuse thermocline in the modeled equatorial Pacific Ocean. Consequently, there have been many attempts to improve the thermocline by developing new or improved ocean vertical mixing schemes. This paper investigates the influence of gradient Richardson Number-based vertical mixing scheme profiles in a tropical Pacific Ocean GCM. It has been common for vertical mixing schemes to be assessed in tropical Pacific Ocean models that have a limited latitudinal domain bounded by zonal walls with sponge layers. However, recent work has shown that warm surface water can accumulate in these models and stop them from achieving the observed sharp equatorial thermocline. The present model employs a parameterized wall heat transport scheme that prevents warm surface water from accumulating. Thus we are able assess the influence of vertical mixing profiles in an ocean model that does not allow warm surface water to accumulate and influence the thermocline. s paper we evaluate the equatorial performance of three different Richardson number (Ri)-based vertical mixing profiles: an integer power (IP) profile based on the observations of Peters, Gregg and Toole; a form of the Pacanowski and Philander profile modified to have low background mixing; and the Max Planck Institute profile. With the accumulation of warm surface water prevented, each of these profiles is able to achieve a sharp thermocline. When compared with observations, the IP profile achieves a better upwelling velocity distribution. We also examine the influence on equatorial performance of very high mixing coefficients at low Richardson number, and of low background mixing coefficients.