1. Introduction
A vertically quasi-uniform layer of temperature (isothermal layer), salinity (isohaline layer), and density (mixed layer) is often observed in the upper ocean. Underneath it, there exists a layer with strong vertical gradient such as the thermocline (in temperature), halocline (in salinity), and pycnocline (in density). The isothermal layer depth (HT) is not necessarily identical to the isohaline layer depth (HS) or the mixed layer depth (HD). There are areas of the World Ocean where HT is deeper than HS or HD (Lindstrom et al. 1987; Lukas and Lindstrom 1991; Sprintall and Tomczak 1992; You 1995; Kara et al. 2000). The layer difference between halocline and thermocline is defined as the barrier layer (BL), which has strong salinity stratification and weak (or neutral) temperature stratification (Fig. 1). The barrier layer thickness (BLT) is originally defined by the difference between HT and HS. However, in recent studies BLT is often referred to as the difference between HT and HD (e.g., Kara et al. 2000).
Less turbulence in the BL than in the isohaline layer due to strong salinity stratification isolates the isohaline layer from cool thermocline water. This process affects the ocean heat budget (Swenson and Hansen 1999) and the heat exchange with the atmosphere. Studies on BL and related physical mechanisms have been devoted to the open ocean such as in the equatorial Pacific with BLT from 20 to 40 m (e.g., Lukas and Lindstrom 1991; You 1998) and the North Pacific with BLT larger than 50 m (Kara et al. 2000). However, little attention has been given to BL that occurs in regional seas. We use the conductivity–temperature–depth (CTD) data from the U.S. Navy's Master Oceanographic Observational Data Set (MOODS) for the Sulu and Celebes Seas to demonstrate the existence and variability of the BL in regional seas. The outline of this paper is as follows. A description of the Sulu/Celebes Seas is given in section 2. The MOODS–CTD data are depicted in section 3. Criteria for determining isothermal and mixed layers are discussed in section 4. A description of BL in the Sulu and Celebes Seas is given in sections 5 and 6. Physical processes for the BL formation are depicted in section 7. In section 8 we present the conclusions.
2. Sulu and Celebes Seas
The Sulu and Celebes Seas are the two major marginal seas in the outer southeastern Asia region (Fig. 2), which consists of an island arc stretching across near 5150 km along the equator at about 15°N–11°S and 94°–141°E, and which has two contrasting zones (Arief 1998). The Sulu Sea belongs to the western zone, where approximately 15%–20% of the global freshwater discharge enters (Toole 1987) through high annual rainfall, between 2 and 4 m yr−1 (ASEAN Subcommittee on Climatology 1982). The Celebes Sea belongs to the eastern zone, which is composed of deep-basin chains of depth 1000–4000 m and filled mainly by western Pacific water masses. This zone has annual rainfall less than 2 m yr−1 (ASEAN Subcommittee on Climatology 1982).
3. MOODS-CTD data
The temperature and salinity data for the Sulu Sea (8°–12°N, 117°–123°E) and the Celebes Sea (0°–6°N, 117°–126°E) acquired from MOODS went through quality control procedures such as a min–max check (e.g., disregarding any temperature data less than −2°C and greater than 35°C), error anomaly check (e.g., rejecting temperature data deviating more than 7°C from climatology), ship-tracking algorithm (screening out data with obvious ship position errors), max number limit (limiting a maximum number of observations within a specified and rarely exceeded space–time window), and buddy check (deleting contradicting data). The climatological data set used for the quality control is the U.S. Navy's Generalized Digital Environmental Model (GDEM) climatological temperature and salinity dataset. After quality control, there are 221 (278) CTD stations in the Sulu (Celebes) Sea. The majority of the CTDs were nominally capable of reaching a maximum depth of 2000 m.
The temporal and spatial distribution of MOODS–CTD data in the Sulu/Celebes Seas is irregular. Certain periods and areas are very well sampled while others lack enough observations to gain any meaningful insights. For example, there is a minimum (maximum) number of CTD stations in the Sulu Sea in May (July). There are no observations in the whole Sulu (Celebes) Sea in certain years (e.g., 1942–46, 1952–63 for the Sulu Sea, and 1979–88 for the Celebes Sea) and there are many observations in other years (e.g., 75 profiles in the Sulu Sea in 1949, and 52 profiles in the Celebes Sea in 1941). Spatial and temporal irregularities along with the lack of data in certain regions must be carefully weighted in order to avoid statistically induced variability.
4. Criteria for determining HT and HD
a. Existing criteria
There are two types of criteria, difference and gradient, for determining HT and HD in the upper ocean. The difference criterion requires the deviation of temperature (density) from its surface value to be smaller than a certain fixed value. The gradient criterion requires the vertical derivative of temperature (density) to be smaller than a certain fixed value.
1) Difference criteria
2) Gradient criteria
b. Criteria used in this study
c. Threshold for determining BL
In the CTD measurements there are generally two thermometers, whose accuracy is about 0.02°C. Thus, the error in computing HT is around 1 m using the gradient criterion (4). The density σt is computed from temperature and salinity. The error in computing HD may be estimated by 2 m. Thus, the threshold for the BL is given as 3 m in this study.
5. Existence of BL in the Sulu Sea
a. General statistics
The two depths, HT and HD, are obtained from 221 CTD profiles using the gradient criteria (4). Monthly CTD stations (denoted by the symbol “○”) and the stations with the BL occurrence (HD < HT, represented by the symbol “+”) in the Sulu Sea (Fig. 3) show a rather frequent occurrence of BL in the Sulu Sea. For example, among 16 (6) CTD stations in September (November), there are 15 (6) stations where a BL occurs. The rate of occurrence reaches 94% in September and 100% in November (Table 1). When BL occurs, the BLT is computed as the difference between HT and HD. The climatological characteristics are outlined as follows: (i) a BL occurs most frequently in September–November (fall) with a frequency of 90%–100% and least frequently in April and May (spring) with a frequency of 50%; (ii) The mixed layer depth (MLD) has a minimum value of 3 m in May and a maximum value of 30.5 m in February; and (iii) the BLT has large values in September–October (39.7–47.6 m) and small values of 9–11.3 m in April and May. It is noticed that the sample size is not sufficiently large; thus the statistical features mentioned here may have bias, especially in May.
b. Seasonal variability
Four groups of CTD data (Fig. 4) with relatively high sampling rates are selected to represent winter (January–February 1949, February–March 1972) and summer (July–September 1949, October 1973). Scatter diagrams of HT and HD (Fig. 5) show that the BL occurs quite often in both summer and winter and is more evident in the summer. For example, the MLD is much less than the isothermal depth (ITD) at almost all stations in October 1973. The maximum thickness of the BL (80 m) is found at Station 6 (see Fig. 4) with a MLD of 18 m and ITD of 98 m. Comparison among the four periods shows the BL has a relatively thin BLT in July–September 1949.
6. Existence of BL in the Celebes Sea
a. General statistics
The two depths (HT and HD) for the Celebes Sea are obtained from processing 179 CTD profiles using the gradient criteria (5). Monthly CTD stations (denoted by the symbol “○”) and the stations with the BL occurrence (HD < HT, represented by the symbol “+”) in the Celebes Sea (Fig. 6) show a rather frequent occurrence of BL in the Celebes Sea. For example, among 14 CTD stations in December, there are 13 stations where BL occurs. The rate of occurrence reaches 93% (Table 2). The BLT has a maximum value of 62.0 m in April and a minimum value of 9.6 m in June. The climatological features of a BL are outlined as follows: (i) the BL occurs most frequently in December with a frequency of 93% and least frequently in April with a frequency of 36%; (ii) the MLD has a minimum value of 7.3 m in March and a maximum value of 28.1 m in September; and (iii) the BLT has a maximum value in March–April (49.7–62.0 m) and a minimum value of 9.6 m in June. These numbers might not have statistical significance because of the small sample size.
b. Seasonal variability
Four groups of CTD data (Fig. 7) with relatively high sampling rates are selected to represent winter (January 1941), spring (May–June 1972), summer (July 1941), and fall (October 1978). Scatter diagrams of HT and HD (Fig. 8) show that the BL occurs quite often in winter and fall and less often in spring. For example, HD is much less than HT at almost all stations in January 1941. At Station 16, HD and HT are 25 and 100 m, respectively, with a BLT about 75 m. Among the four periods, the BL has relatively thin BLT in July 1941 and May–June 1972.
7. Discussion of physical processes causing the BL formation
a. Two mechanisms for the occurrence of BL
Two major factors determining HT (HD) are 1) surface winds and net heat (buoyancy) flux and 2) thermal (density) stratification underneath the isothermal (mixed) layer. Thus, the condition (HT > HD) is caused by the surface freshwater flux (precipitation excess over evaporation) and strong salinity stratification underneath the mixed layer. Thus, the BL can be classified as (i) rain formed, and (ii) stratification formed (weaker thermal and strong salinity stratification).
1) Rain-formed mechanism
2) Stratification-formed mechanism
b. Sulu Sea
The Comprehensive Ocean–Atmosphere Data Set (COADS) provides monthly mean sea surface wind stress and net heat and freshwater fluxes (Fig. 9). The surface wind stress, net heat flux (F), and net freshwater flux (P − E) have seasonal variability in the Sulu Sea. The wind stress strengthens in December–February [6.57–6.67 (×10−2 N m−2)] and weakens in April–November [0.63–3.66 (×10−2 N m−2)]. The net heat flux is negative in November–January (−25.4 to 0 W m−2, surface cooling) and becomes large positive in March–May (60–110.3 W m−2, strong surface warming). The surface freshwater flux (P − E) is weak in January–May (<0.6 m yr−1) and strong in June–November (>1 m yr−1).
In April and May the Sulu Sea is under the influence of weak winds [1.3–3 (×10−2 N m−2)], strong net heat flux (60–110.3 W m−2), and weak surface freshwater flux (0.52–0.56 m yr−1). The strong heat flux and weak winds cause a shallow isothermal layer (small HT). The weak freshwater flux cannot strengthen the salinity stratification. Thus, the BL occurs least frequently in April and May.
In September and October the Sulu Sea is under influence of weak winds [0.63–2.19 (×10−2 N m−2)], weak net heat flux (24.6–29.2 W m−2), and strong surface freshwater flux (3.1 m yr−1 in September). These conditions lead to the detrainment regime for the upper ocean and favor the rain-formed mechanism (excess precipitation over evaporation).
In December the Sulu Sea is under influence of strong winds (6.57 × 10−2 N m−2), surface cooling (−25.4 W m−2), and relatively weak surface freshwater flux (0.8 m yr−1). These conditions favor the entrainment regime in the upper ocean and the stratification-formed mechanism for the occurrence of a BL.
c. Celebes Sea
The surface wind stress, net heat flux (F), and net freshwater flux (P − E) have less seasonal variability in the Celebes Sea than in the Sulu Sea (Fig. 9). The wind stress varies between 0.45 × 10−2 (Nov) and 2.6 × 10−2 N m−2 (Feb). The surface net heat and freshwater fluxes are positive year-round with F varying from 17.5 (Jan) to 87.4 W m−2 (Mar) and (P − E) varying from 0.17 (Apr) to 2.1 m yr−1 (Jan). Thus, the Celebes Sea is under the influence of weak winds, surface warming, and surface freshwater flux all year. These conditions favor the rain-formed mechanism.
In December–January the Celebes Sea experiences weak net heat flux (17.5–24.5 W m−2) and strong surface freshwater flux (1.14–2.06 m yr−1). These conditions may cause the detrainment regime for the upper ocean and favor the rain-formed mechanism (excess precipitation over evaporation).
In April, the Celebes Sea experiences strong net heat flux (80.4 W m−2) and weak surface freshwater flux (0.17 m yr−1). The strong heat flux and weak winds cause a shallow isothermal layer (small HT). The weak freshwater flux strengthens the salinity stratification slightly. Thus, the BL occurs least frequently in April.
8. Summary
Atmospheric forcing (rain-formed mechanism) and ocean physical processes (stratification-formed mechanism) may cause the occurrence of a barrier layer (deeper thermocline than pycnocline) in the Sulu/Celebes Seas. Using the U.S. Navy's Master Oceanographic Observational Data Set, we have shown barrier layer thickness up to 47.6 m (62.0 m) in the Sulu (Celebes) Sea.
The barrier layer appears quite often in both Sulu and Celebes Seas. The estimated occurrence rate varies from 50% (April and May) to 100% (November) in the Sulu Sea and from 36% (April) to 93% (December) in the Celebes Sea.
In the Sulu Sea, the mixed layer depth has a minimum value of 3 m in May and a maximum value of 30.5 m in February. The barrier layer thickness has large values in September–November (39.7–47.6 m) and small values of 9–11.3 m in April and May.
In the Celebes Sea, the mixed layer depth has a minimum value of 7.3 m in March and a maximum value of 28.1 m in September. The barrier layer thickness has large values in March–April (49.7–62.0 m) and small value of 9.6 m in June. Relatively higher precipitation and lower wind stress than for the Sulu Sea are found in winter (see Fig. 9).
In the Sulu Sea, both rain-formed and stratification-formed mechanisms are responsible for the barrier layer formation. However, in the Celebes Sea the rain-formed mechanism is responsible for the barrier layer formation.
Acknowledgments
This work was funded by the Naval Oceanographic Office, the Office of Naval Research, and the Naval Postgraduate School. Liu was also funded by the Chinese National Key Program for Developing Basic Science (G1999043807) and Chinese Special Program “2001DIA50041” of the National Ministry of Science and Technology.
REFERENCES
Arief, D., 1998: Outer Southeast Asia: A region of deep straits, including the Benda Sea coastal segment (13°S). The Sea, A. R. Robinson and K. H. Brink, Eds., Regional Studies and Syntheses, Vol. 11, John Wiley and Sons, 507–522.
ASEAN Subcommittee on Climatology, 1982: The ASEAN Climatic Atlas. Directorate of National Mapping, Kuala Lumpur, Malaysia, 104 pp.
Bathen, K. H., 1972: On the seasonal changes in the depth of the mixed layer in the North Pacific Ocean. J. Geophys. Res., 77 , 7138–7150.
Chu, P. C., 1993: Generation of low frequency unstable modes in a coupled equatorial troposphere and ocean mixed layer. J. Atmos. Sci., 50 , 731–749.
Chu, P. C., and R. W. Garwood Jr.,, 1991: On the two-phase thermodynamics of the coupled cloud-ocean mixed layer. J. Geophys. Res., 96 , 3425–3436.
Chu, P. C., R. W. Garwood Jr.,, and P. Muller, 1990: Unstable and damped modes in coupled ocean mixed layer and cloud models. J. Mar. Syst., 1 , 1–11.
Defant, A., 1961: Physical Oceanography. Vol. 1. Pergamon, 729 pp.
Kara, A. B., P. A. Rochford, and H. E. Hurlburt, 2000: Mixed layer depth variability and barrier layer formation over the North Pacific Ocean. J. Geophys. Res., 105 , 16783–16801.
Lindstrom, E., R. Lukas, R. Fine, E. Firing, S. Godfrey, G. Meyers, and M. Tsuchiya, 1987: The Western Equatorial Pacific Ocean Circulation Study. Nature, 330 , 533–537.
Lukas, R., and E. Lindstrom, 1991: The mixed layer of the western equatorial Pacific Ocean. J. Geophys. Res., 96 , 3343–3357.
Miller, J. R., 1976: The salinity effect in a mixed layer ocean model. J. Phys. Oceanogr., 6 , 29–35.
Monterey, G., and S. Levitus, 1997: Seasonal Variability of Mixed Layer Depth for the World Ocean. NOAA Atlas NESDIS 14, 100 pp.
Ohlmann, J. C., D. A. Siegel, and C. Gautier, 1996: Ocean mixed layer radiant heating and solar penetration: A global analysis. J. Climate, 9 , 2265–2280.
Spall, M. A., 1991: A diagnostic study of wind- and buoyancy-driven North Atlantic circulation. J. Geophys. Res., 96 , 18509–18518.
Sprintall, J., and M. Tomczak, 1992: Evidence of barrier layer in the surface layer of tropics. J. Geophys. Res., 97 , 7305–7316.
Swenson, M. S., and D. V. Hansen, 1999: Tropical Pacific Ocean mixed layer budget: The Pacific cold tongue. J. Phys. Oceanogr., 29 , 69–81.
Toole, J. M., 1987: WOCE, interbasin exchanges, and marginal sea overflows. Eos, Trans. Amer. Geophys. Union, 68 (1) 2–3. 11.
Wyrtki, K., 1964: The thermal structure of the eastern Pacific Ocean. Dstch. Hydrogr. Z., 8A , (Suppl.),. 6–84.
You, Y., 1995: Salinity variability and its role in the barrier-layer formation during TOGA–COARE. J. Phys. Oceanogr., 25 , 2778–2807.
You, Y., 1998: Rain-formed barrier layer of the western equatorial Pacific warm pool: A case study. J. Geophys. Res., 103 , 5361–5378.
Monthly occurrence frequency and climate features of the barrier layer in the Sulu Sea
Monthly occurrence frequency and climate features of the barrier layer in the Celebes Sea