Regional primitive equation modeling and analysis of the polymode data set

Abstract

A regional, hybrid coordinate, primitive equation (PE) model is applied to a 60-day period of the POLYMODE data set. The initialization techniques and open boundary conditions introduced by Spall and Robinson are shown to produce stable, realistic, and reasonably accurate hindcasts for the 2-month data set. Comparisons with quasi-geostrophic (QG) modeling studies indicate that the PE model reproduced the jet formation that dominates the region more accurately than did the QG model. When the PE model used boundary conditions that were partially adjusted by the QG model, the resulting fields were very similar to the QG fields, indicating a rapid degradation of small-scale features near the boundaries in the QG calculation. A local term-by-term primitive equation energy and vorticity analysis package is also introduced. The full vorticity, horizontal divergence, kinetic energy, and available gravitational energy equations are solved diagnostically from the output of the regional PE model. Through the analysis of a time series of horizontal maps, the dominant processes in the flow are illustrated. The individual terms are also integrated over the region of jet formation to highlight the net balances as a function of time. The formation of the deep thermocline jet is shown to be due to horizontal advection through the boundary, baroclinic conversion in the deep thermocline and vertical pressure work, which exports the deep energy to the upper thermocline levels. It is concluded here that the PE model reproduces the observed jet formation better than the QG model because of the increased horizontal advection and stronger vertical pressure work. Although the PE model is shown to be superior to the QG model in this application, it is believed that both PE and QG models can play an important role in the regional study of mid-ocean mesoscale eddies.