The analysis of large-scale turbulence characteristics in the Indonesian seas derived from a regional model based on the Princeton Ocean Model

Turbulence characteristics in the Indonesian seas on the horizontal scale of order of 100 km were calculated with a regional model of the Indonesian seas circulation in the area based on the Princeton Ocean Model (POM). As is well known, the POM incorporates the Mellor–Yamada turbulence closure scheme. The calculated characteristics are: twice the turbulence kinetic energy per unit mass, <i>q</i>²; the turbulence master scale, &ell;; mixing coefficients of momentum, <i>K</i>_M; and temperature and salinity, <i>K</i>_H; etc. The analyzed turbulence has been generated essentially by the shear of large-scale ocean currents and by the large-scale wind turbulence. We focused on the analysis of turbulence around important topographic features, such as the Lifamatola Sill, the North Sangihe Ridge, the Dewakang Sill, and the North and South Halmahera Sea Sills. In general, the structure of turbulence characteristics in these regions turned out to be similar. For this reason, we have carried out a detailed analysis of the Lifamatola Sill region because dynamically this region is very important and some estimates of mixing coefficients in this area are available. <br><br> Briefly, the main results are as follows. The distribution of <i>q</i>² is quite adequately reproduced by the model. To the north of the Lifamatola Sill (in the Maluku Sea) and to the south of the Sill (in the Seram Sea), large values of <i>q</i>² occur in the deep layer extending several hundred meters above the bottom. The observed increase of <i>q</i>² near the very bottom is probably due to the increase of velocity shear and the corresponding shear production of <i>q</i>² very close to the bottom. The turbulence master scale, &ell;, was found to be constant in the main depth of the ocean, while &ell; rapidly decreases close to the bottom, as one would expect. However, in deep profiles away from the sill, the effect of topography results in the &ell; structure being unreasonably complicated as one moves towards the bottom. Values of 15 to 20 × 10^&minus;4 m² s^-1 were obtained for <i>K</i>_M and <i>K</i>_H in deep water in the vicinity of the Lifamatola Sill. These estimates agree well with basin-scale averaged values of 13.3 × 10^&minus;4 m² s^-1 found diagnostically for <i>K</i>_H in the deep Banda and Seram Seas (Gordon et al., 2003) and a value of 9.0 × 10^&minus;4 m² s^-1 found diagnostically for <i>K</i>_H for the deep Banda Sea system (van Aken et al., 1988). The somewhat higher simulated values can be explained by the presence of steep topography around the sill.

Validation of a regional Indonesian seas model based on a comparison between model and INSTANT transports

The International Nusantara Stratification and Transport (INSTANT) program measured currents through multiple Indonesian Seas passages simultaneously over a three-year period (from January 2004 to December 2006). The Indonesian Seas region has presented numerous challenges for numerical modelers - the Indonesian Throughflow (ITF) must pass over shallow sills, into deep basins, and through narrow constrictions on its way from the Pacific to the Indian Ocean. As an important region in the global climate puzzle, a number of models have been used to try and best simulate this throughflow. In an attempt to validate our model, we present a comparison between the transports calculated from our model and those calculated from the INSTANT in situ measurements at five passages within the Indonesian Seas (Labani Channel, Lifamatola Passage, Lombok Strait, Ornbai Strait, and Timor Passage). Our Princeton Ocean Model (POM) based regional Indonesian Seas model was originally developed to analyze the influence of bottom topography on the temperature and salinity distributions in the Indonesian seas region, to disclose the path of the South Pacific Water from the continuation of the New Guinea Coastal Current entering the region of interest up to the Lifamatola Passage, and to assess the role of the pressure head in driving the ITF and in determining its total transport. Previous studies found that this model reasonably represents the general long-term flow (seasons) through this region. The INSTANT transports were compared to the results of this regional model over multiple timescales. Overall trends are somewhat represented but changes on timescales shorter than seasonal (three months) and longer than annual were not considered in our model. Normal velocities through each passage during every season are plotted. Daily volume transports and transport-weighted temperature and salinity are plotted and seasonal averages are tabulated.