Observations of Inflow of Philippine Sea Surface Water into the South China Sea through the Luzon Strait

Transcription

1 113 Observations of Inflow of Philippine Sea Surface Water into the South China Sea through the Luzon Strait LUCA R. CENTURIONI AND PEARN P. NIILER Scripps Institution of Oceanography, La Jolla, California DONG-KYU LEE Busan National University, Busan, South Korea (Manuscript received 19 February 2003, in final form 12 June 2003) ABSTRACT Velocity observations near the surface made with Argos satellite-tracked drifters between 1989 and 2002 provide evidence of seasonal currents entering the South China Sea from the Philippine Sea through the Luzon Strait. The drifters cross the strait and reach the interior of the South China Sea only between October and January, with ensemble mean speeds of ms 1 and daily mean westward speeds that can exceed 1.65 ms 1. The majority of the drifters that continued to reside in the South China Sea made the entry within a westward current system located at 20 N that crossed the prevailing northward Kuroshio path. In other seasons, the drifters looped across the strait within the Kuroshio and exited along the south coast of Taiwan. During one intrusion event, satellite altimeters indicated that, directly west of the strait, anticyclonic and cyclonic eddies resided, respectively, north and south of the entering drifter track. The surface currents measured by the crossing drifters were much larger than the Ekman currents that would be produced by an 8 10 m s 1 northeast monsoon, suggesting that a deeper westward current system, as seen in historical watermass analyses, was present. 1. Introduction This work describes observations of velocity made at a nominal depth of 15 m with satellite-tracked drifters and provides new direct evidence of seasonal near-surface flow from the Philippine Sea into the South China Sea through the Luzon Strait (Fig. 1). Such flow is frequently described as a westward branch of the Kuroshio. Mainly hydrographic methods have provided evidence of the seasonal penetration of the Philippine Sea water into the South China Sea (e.g., Fang et al. 1998). The surface circulation south and east of the Luzon Strait is dominated by strong and persistent subtropical current systems. At the surface, the yearly mean Pacific North Equatorial Current bifurcates at 13 N near the east coast of Luzon to form the northward-flowing Kuroshio and the southward-flowing Mindanao Current (Nitani 1972; Toole et al. 1990; Qu and Lukas 2003). The near-surface bifurcation latitude moves between 11 N in May and 14.5 N in November, and at depth it is even farther north of its surface expression (Nitani 1972; Qu and Lukas 2003). Qiu and Lukas (1996) have noted that the interannual variations of the bifurcation Corresponding author address: Luca R. Centurioni, Scripps Institution of Oceanography, 9500 Gilman Dr., La Jolla, CA lcenturioni@ucsd.edu latitude can be successfully computed from the winddriven vorticity dynamics of linear and nonlinear reduced-gravity circulation models. At 18 N the Kuroshio is a well-formed, northwardflowing western boundary current concentrated entirely west of 124 E; its high-speed core is positioned at 123 E, and its baroclinic structure is evident in the upper 600 m (Toole et al. 1990; Qu et al. 1998). Before reaching Taiwan, the Kuroshio encounters the Luzon Strait, which is the deepest passage from the Pacific Ocean to the South China Sea. At the southern portion of the strait, the Kuroshio takes a westward set and makes a detour into the South China Sea through the deepest channels of the Luzon Strait: the Balintany Channel and south of Babuyan Island. West and north of Batan Island, the Kuroshio flows within the Bashi Channel until it reaches the southeast coast of Taiwan (Gilson and Roemmich 2002). Between the north coast of Luzon and the southeast coast of Taiwan, the Kuroshio is occasionally referred to as a loop current because it makes a loop, or excursion, into the South China Sea. The largest loop occurs between October and January when the winds are dominated by the northeast monsoon. A description of the loop current can be found, for example, in Nitani s (1972) maps of the geomagnetic electrokinetograph data and dynamic heights compiled from several hydrographic sur American Meteorological Society

2 114 JOURNAL OF PHYSICAL OCEANOGRAPHY VOLUME 34 FIG. 1. Bathymetry of the study region (Smith and Sandwell 1997, bathymetry version 8.2). veys. Li et al. (1998) used hydrographic data to show the looping of the Kuroshio in the Luzon Strait in summer months in the upper 200 m. Between October and March the warm surface water of Philippine Sea origin contrasts with the surrounding colder South China Sea water. This condition was exploited by Farris and Wimbush (1996), who observed the occurrence of Kuroshio surface loops and intrusions in autumn and winter with sea surface temperature (SST) maps made from Advanced Very High Resolution Radiometer data. As defined by the SST, the loop intrudes most severely into the South China Sea during the October January period, when significant amounts of Philippine Sea water are also found below the surface. For example, Shaw (1989) used conductiv- TABLE 1. Summary of the Lagrangian statistics in the region bounded, in the zonal direction, by 120 and 135 E and, in the meridional direction, by 10 and 25 N. Here T is the Lagrangian timescale, L is the Lagrangian length scale, and denotes the standard deviation. The first figure in the columns labeled with N T is the number of 6-h-interval time series longer than 20 days for which the velocity autocovariance converges, and the second figure is the total number of time series examined. Here E and N in the subscripts refer to the zonal and meridional direction, respectively, and N P is the total number of velocity observations in the region. JFM is for Jan Mar, AMJ is for Apr Jun, JAS is for Jul Sep, OND is for Oct Dec, and All is for all months. Months T E (days) (T E ) (days) L E (10 4 m) (L E ) T N (10 4 m) N TE (days) (T N ) (days) L N (10 4 m) (L N ) (10 4 m) N TN N P JFM AMJ JAS OND All /95 93/96 98/102 90/90 259/ /95 93/96 96/102 88/90 255/

3 115 FIG. 2. Six-hour-interval positions of the drifters, color-coded in accordance with the local instantaneous speed. The data shown here were collected between Sep 1987 and May The black lines represent the 500-m depth contours. ity temperature depth (CTD) profiles to demonstrate that Philippine Sea water extends down to 500 m on the continental slope west of the southern tip of Taiwan in the March August period. Shaw (1991), by analysis of a larger hydrographic dataset, tracked water of Kuroshio origin as far as west of 115 E in the upper 250 m and along the continental slope of the northern South China Sea. Qu (2002) reached similar conclusions about intrusions upon inspection of the distribution of oxygen concentration. Although the inflow is apparent from watermass properties, diapycnal mixing in the upper South China Sea and the outflow needed to achieve a mass balance of the basin are not known well. Multilayer ocean general circulation model simulations of the transport of upper-layer water through the Luzon Strait are discussed in Metzger and Hurlburt (1996, 2001). In these simulations, the inflow is a result of the seasonally varying basin-scale circulation of the entire western subtropical and tropical Pacific. These works associate the seasonal inflow with the lowering of sea level caused by the increased cyclonic circulation of the South China Sea that results from the positive wind stress curl of the northeast monsoon. Sheremet (2001) used a barotropic model to address the local parameter dependence of a western boundary current that encounters a gap along the coast. The parameters are the breadth of the gap relative to the width of the Munk viscous boundary current and the Reynolds number of the flow, defined as the ratio between the transport rate of mass per unit depth and the lateral eddy diffusion coefficient. The general theoretical results are that a strong current leaps the gap and a weak current loop propagates westward into the gap, forming eddies in the intrusion. Even in the absence of direct evidence of the circulation patterns (but see Farris and Wimbush 1996), observations and theory are supporting the view that water of Philippine Sea origin can reach the interior of the South China Sea. In numerical studies of the circulation of the South China Sea, the processes that cause the inflow appear to be associated with, or produce, mesoscale eddies. Observational evidence for this mesoscale variability has recently become available from hydrography (Li et al. 1998). In this work, we present and discuss observations of the surface circulation in the western tropical Pacific and in the South China Sea and we report on the observed flow from the Philippine Sea to the South China Sea. 2. The data Direct velocity measurements in the surface mixed layer were obtained with Argos satellite-tracked drifters

4 116 JOURNAL OF PHYSICAL OCEANOGRAPHY VOLUME 34 FIG. 3. Seasonal displacement diagrams of the drifters: (a) Jan Mar, (b) Apr Jun, (c) Jul Sep, and (d) Oct Dec. Contours of bottom depth are 200 m. drogued at a nominal depth of 15 m. Between November of 1986 and May of 2002, 380 instruments were deployed or drifted into the region delimited by E and N. Only 14 deployments were made in the South China Sea. For a description of the drifters, wind slip correction, and data-processing methods see Niiler (2001). In this analysis, the positions of each drifter, interpolated to 6-h intervals, were used to construct the velocity vectors, which are called instantaneous velocities. The estimated accuracy of the velocity measurements in a 10 m s 1 wind is 10 2 ms 1 (Niiler et al. 1995). The Lagrangian statistics and the ensemble mean velocity field were calculated in the region delimited by 120 and 135 E longitude and by 10 and 25 N latitude. For each drifter trajectory that was longer than 20 days, we computed the Lagrangian decorrelation timescale and length scale in the zonal and in the meridional directions (e.g., Freeland et al. 1975; Poulain and Niiler 1989). Nonconverging autocovariance functions (the ones without zero crossings) were discarded. Seasonal Lagrangian timescale and length scale were also computed. The results are summarized in Table 1. The dimensions of each averaging cell were chosen to be the eastward and northward Lagrangian space scales, respectively (and up to 2 times these values for the sparser seasonal datasets). Within each cell, the independent estimates of the zonal and meridional velocities were computed by time averaging the 6-h velocity time series in bins with amplitude given by the upper bound of the largest Lagrangian timescale [T (T) in Table 1]. The ensemble mean velocity and its standard error ellipse referred to the principal axes (Freeland et al. 1975) were calculated for each cell. The North Equatorial Current bifurcation latitudes in each season were chosen where the zonal average (between 127 and 130 E) of the seasonal ensemble mean northward velocity components was zero. Unless stated otherwise, all uncertainties are assumed to be 1 standard deviation. 3. Flow of the Philippine Sea water through the Luzon Strait The ensemble of the individual drifter tracks (Fig. 2) shows how the North Equatorial Current intensifies to

5 117 FIG. 4. Displacement diagram of the drifters that crossed the line E (the dashed line labeled D ) at least once; 29 drifters satisfied this condition. The trajectories of the drifters that moved back to the Pacific Ocean are plotted in red. The others reached the interior of the South China Sea and are plotted in blue. See main text for the meaning of A, B, and C. The 200-m contours of bottom depth are shown with a solid black line. form the Kuroshio, made visible by instantaneous speeds in excess of 0.7 m s 1, east of Luzon between 16.5 and 18.5 N. To the north, the drifters in this flow move northwest and enter the Luzon Strait principally through the channels north (Balitany Channel) and south of Babuyan Island. Two drifters moved very rapidly westward on the continental shelf north of Luzon. Some drifters within the South China Sea continued to exhibit instantaneous speeds in excess of 0.6 m s 1 within the current system that appears to surround its entire perimeter. The seasonal displacement diagram of the drifters (Fig. 3) shows that a westward flow through the Luzon Strait occurs only between October and March. Of the 29 drifters that crossed the meridional line at E between Luzon and Taiwan at least once, 15, as shown in blue in Fig. 4, (5 of which had lost the drogue) moved to the interior of the South China Sea and 6 reached the broad continental-shelf area south of Vietnam. Daily average speeds in excess of 0.8 m s 1 were observed when the drifters traveled southward over the narrow shelf area located along the southeast Vietnamese coast (marked A in Fig. 4). The drifters whose positions are plotted in red on Fig. 4 recrossed the strait eastward and moved north, some of them following a small loop pattern southwest of the south coast of Taiwan (marked B in Fig. 4). Five of these drifters had lost the drogue, implying that the position error accumulates at a rate of at most 10 km day 1 downwind because of the slip of the float relative to a drogued drifter (Niiler 2001), an error that tends to push the tracks southwestward during the northeast monsoon. Some of the red drifters were trapped in a low-velocity region (marked C in Fig. 4) near the southeast tip of Taiwan before rejoining the northward Kuroshio flow. In the proximity of the Luzon Strait there were data in 13 of the 17 yr analyzed. In 10 of the 13 yr, the drifters crossed the strait between October and January and remained in the South China

6 118 JOURNAL OF PHYSICAL OCEANOGRAPHY VOLUME 34 FIG. 5. Daily average velocity for drifters crossing the line E. The thick vectors correspond to the drifters that reached the interior of the South China Sea (with a mean speed of m s 1 ). The thinner vectors correspond to the drifters that crossed the line back toward the Philippine Sea. The thin dashed lines are the bottom depth (200 m) contours. Sea (the blue drifters in Fig. 4). In 8 of the 13 yr, the drifters made the small loop across the Luzon Strait between October and March (the red drifters in Fig. 4). Five drifters (all with their drogues attached) reached the strait by traveling westward and crossing directly over the mean Kuroshio path. Most of the trajectories north of 15 N (Fig. 4) indicate that besides the seasonal westward current into the Luzon Strait there is, to the east of the strait, a vigorous eddy field. The velocities at the crossing line are plotted in Fig. 5. The average speed computed only with data from the drifters that reached the interior of the South China Sea is ms 1. The maximum, almost westward, speed of m s 1 (the instrumental error is used here) was observed at 20.7 N, E in December of The ensemble mean velocity field, computed with a resolution of approximately 0.5 from all of the data east and southeast of the Luzon Strait (Fig. 6), shows that the Kuroshio system is relatively well mapped. The northward boundary current forms near the Philippines coast where the North Equatorial Current separates in two branches. The latitude of this bifurcation point, computed from the annual mean surface (15 m) velocity field, is located at N (dashed line in Fig. 6), and it varies seasonally by , reaching its northernmost position of N in the October December period (not shown), in general agreement with the results of Qu and Lukas (2003). The Kuroshio accelerates as mass is added from the east, and, by 14.8 N, the average speed increases above m s 1 (computed in box 1, Fig. 6). A second concentrated stream appears to join the Kuroshio from the east at about 20 N. The mean westward speed of this current is m s 1 (computed in box 2, Fig. 6) and increases to m s 1 if only autumn data are considered (average computed in a slightly larger box, not shown). This stream is mainly sampled by the drifters that crossed the Luzon Strait and moved into the South China Sea in the October January period, as depicted and discussed in Figs. 3 and 4. This relatively strong current that crosses the Kuroshio path appears in both the global and autumn ensemble means. However, the interannual variability of this westward current structure is not well determined by the drifter data, implying that the estimate of the standard error is based on potentially time-aliased data and should be considered as a lower bound. The mesoscale activity between 15 December 1997 and 18 January 1998 can be seen in the sea level anomaly of Centre National d études Spatiales Archivage, Validation et Interprétation des Données des Satellites Océanographiques (CNES/AVISO) Ocean Topography Experiment (TOPEX)/Poseidon European Remote Sensing Satellite (ERS) merged altimeter data of sea surface height anomaly (SLA) and the Tropical Rainfall Measuring Mission Microwave Imager (TMI) 3-dayaveraged SST data (Fig. 7). During this period, one drifter crossed the Luzon Strait. The water warmer than 26 C from the Philippine Sea intruded north of 20 N and is visible in Fig. 7 because the surrounding South China Sea surface water was colder than 24 C. The most obvious intrusion of warm water along the northern boundary of the South China Sea extended westward to E, with a weaker continuation east of 116 E. An anticyclonic eddy (A) was centered at approximately (21.6 N, E) near the 200-m isobath. In the Luzon Strait and in the South China Sea, south of approximately 20.5 N, the sea level was lower than the average and a cyclonic eddy was located directly north of Luzon (B). The drifter track in Fig. 7b crosses the Luzon Strait westward along the border between the cold and the warm water that lies between the anticyclonic and cyclonic eddies at approximately 20.5 N. The trajectory in Fig. 7b is shown 5 days before and after the date of the SLA and SST observations. The average

7 119 FIG. 6. Mean velocity and standard-error ellipses plotted on the same scale: Nov 1986 May Averaging grid (lon lat) is , approximately corresponding to the Lagrangian length scales. Here out of measurements are independent (based on a decorrelation timescale of 5.3 days). Gray cells contain less than 15 (at least 3) independent measurements. The 500-m contours of bottom depth are shown with a thin dashed black line. The latitude of the North Equatorial Current bifurcation point is ( ) N and is marked with the thick dashed black line. The mean velocities were computed in boxes 1 and 2 (see main text). speed of the drifter for the first 5 days, west of 118 E, was m s 1. For the last 5 days, east of 118 E, when the drifter began to leave the warm water, the average speed was m s 1. Figures 7c and 7d show that the warm-water intrusion separated from the warm water along the southern coast of Taiwan by being entrapped in the westward-moving anticyclonic eddy. The drifter moved to the south and west in a cyclonic circulation around a cyclonic feature (C) that propagated westward in the northern part of the South China Sea basin. 4. Discussion and summary Velocity observations made with drifters drogued to 15-m depth describe a seasonally recurring inflow of surface water from the Philippine Sea into the South China Sea through the Luzon Strait between October and January. The permanent looplike intrusion of the Kuroshio described, for example, by Nitani (1972) and, with numerical models, Metzger and Hurlburt (2001) does not appear in the drifter data between April and September. The loop is confined to the north of the Luzon Strait in the remaining months, when its westward extension reaches to E. A net westward Ekman transport through the Strait of Luzon, computed from Special Sensor Microwave Imager (SSM/I)-derived winds (Fig. 8) exists only during the northeast monsoon, reaching a maximum of 0.6 Sv (1 Sv 10 6 m 3 s 1 ). Figure 8 also shows that the drifters enter the South China Sea preferentially during the early part of the onset of the northeast monsoon. The drifters that crossed the strait were entrained by the Kuroshio east of Luzon or crossed the Kuroshio path inside the Luzon Strait within a strong westward jet. The ensemble mean velocity field computed between October and December (not shown) indicates that a nearly zonal westward flow of m s 1 occurred east of the Luzon Strait at 20 N. The sampling of the drifters was not uniform in time; however, this jet appears between 18.5 and 20.5 N in all of the years (1989, 1992, 1994, and 1995) during which there were data in the region. The

8 120 JOURNAL OF PHYSICAL OCEANOGRAPHY VOLUME 34 FIG. 7. Maps of SLA [10-day average, contoured with thick dashed (continuous) lines for negative (positive) values, in meters] and SST (3-day average): (a) 15 Dec 1997, (b) 25 Dec 1997, (c) 4 Jan 1998, and (d) 14 Jan The purple marks show the trajectory of a drifter (5 days before and after the nominal date). Contours of bottom depth (200 m) are shown with a thin black dashed line. near-surface Ekman currents produced by 10 m s 1 northeast monsoon winds would be at most 0.1 m s 1 (Ralph and Niiler 1999). The much stronger surface currents observed with the drifters indicate that a deeper current system to the west must be present during the FIG. 8. Ekman transport across the Luzon Strait into the South China Sea computed between 1992 and 2001 from monthly average SSM/I satellite winds (solid line, right scale), and number of drifters per month entering the South China Sea from Kuroshio (histogram, left scale). Two yearly cycles are repeated. northeast monsoon [as hypothesized, e.g., by Shaw (1989, 1991)] and that, in the Luzon Strait region, entrains Kuroshio water. The study of Sheremet (2001) predicts that the Kuroshio can penetrate the South China Sea when its strength is reduced. This weakening occurs in the autumn months (see, e.g., Qu et al. 1998), when the North Equatorial Current system shifts northward. The depth of the penetration and the flow regime must also depend on the topography of the Luzon Strait, a problem not dealt with well by Sheremet s (2001) theory or the low-vertical-resolution models of Metzger and Hurlburt (1996, 2001). Acknowledgments. The help of Sharon Lukas, Mayra Pazos, and Jessica Redman for processing and recovering the most recent drifter data is gratefully acknowledged. TMI data and images are produced by Remote Sensing Systems and are sponsored by NASA s Earth Science Information Partnerships (ESIP; a federation of information sites for earth science) and by NASA s TRMM Science Team. SSM/I data and images are produced by Remote Sensing Systems and are sponsored by the NASA Pathfinder Program for early Earth Observing System (EOS) products. This work was supported at the Scripps Institution of Oceanography by NOAA Grant NOAA-NA-17R1231.

9 121 REFERENCES Fang, G., W. Fang, Y. Fang, and K. Wang, 1998: A survey of studies on the South China Sea upper ocean circulation. Acta Oceanogr. Taiwan, 37, Farris, A., and M. Wimbush, 1996: Wind-induced intrusion into the South China Sea. J. Oceanogr., 52, Freeland, H. J., P. B. Rhines, and T. Rossby, 1975: Statistical observations of the trajectory of neutrally buoyant floats in the North Atlantic. J. Mar. Res., 33, Gilson, J., and D. Roemmich, 2002: Mean and temporal variability in Kuroshio geostrophic transport south of Taiwan ( ). J. Oceanogr., 58, Li, L., W. D. Nowlin, and S. Jilan, 1998: Anticyclonic rings from the Kuroshio in the South China Sea. Deep-Sea Res., 45A, Metzger, E. J., and H. E. Hurlburt, 1996: Coupled dynamics of the South China Sea, the Sulu Sea, and the Pacific Ocean. J. Geophys. Res., 101, , and, 2001: The nondeterministic nature of Kuroshio penetration and eddy shedding in the South China Sea. J. Phys. Oceanogr., 31, Niiler, P. P., 2001: The World Ocean surface circulation. Ocean Circulation and Climate Observing and Modeling the Global Ocean, J. Church, G. Siedler, and J. Gould, Eds., Academic Press, , A. S. Sybrandy, K. Bi, P. M. Poulain, and D. Bitterman, 1995: Measurements of the water following capability of holey-sock and TRISTAR drifters. Deep-Sea Res., 42A, Nitani, H., 1972: Beginning of the Kuroshio. Kuroshio: Its Physical Aspects, H. Stommel and K. Yashida, Eds., University of Tokyo Press, Poulain, P. M., and P. P. Niiler, 1989: Statistical analysis of the surface circulation in the California current system using satellitetracked drifters. J. Phys. Oceanogr., 19, Qiu, B., and R. Lukas, 1996: Seasonal and interannual variability of the North Equatorial Current, the Mindanao Current, and the Kuroshio along the Pacific western boundary. J. Geophys. Res., 101, Qu, T., 2002: Evidence for water exchange between the South China Sea and the Pacific Ocean through the Luzon Strait. Acta Oceanol. Sinica, 21, , and R. Lukas, 2003: The bifurcation of the North Equatorial Current in the Pacific. J. Phys. Oceanogr., 33, 5 18., H. Mitsudera, and T. Yamagata, 1998: On the western boundary currents in the Philippine Sea. J. Geophys. Res., 103, Ralph, E. A., and P. P. Niiler, 1999: Wind-driven currents in the tropical Pacific. J. Phys. Oceanogr., 29, Shaw, P., 1989: The intrusion of water masses into the sea southwest of Taiwan. J. Geophys. Res., 94, , 1991: The seasonal variation of the intrusion of the Philippine Sea water into the South China Sea. J. Geophys. Res., 96, Sheremet, V., 2001: Hysteresis of a western boundary current leaping across a gap. J. Phys. Oceanogr., 31, Smith, W. H. F., and D. T. Sandwell, 1997: Global seafloor topography from satellite altimetry and ship depth soundings. Science, 277, Toole, J. M., R. C. Millard, Z. Wang, and S. Pu, 1990: Observations of the Pacific North Equatorial Current bifurcation at the Philippine coast. J. Phys. Oceanogr., 20,