• 沒有找到結果。

Evaluation of Tidal Removal Method Using Phase Average Technique from ADCP Surveys along the Peng-Hu Channel in the Taiwan Strait

N/A
N/A
Protected

Academic year: 2021

Share "Evaluation of Tidal Removal Method Using Phase Average Technique from ADCP Surveys along the Peng-Hu Channel in the Taiwan Strait"

Copied!
11
0
0

加載中.... (立即查看全文)

全文

(1)

NOTES AND CORRESPONDENCE

Evaluation of Tidal Removal Method Using Phase Average Technique

from ADCP Surveys along the Peng-Hu Channel in the Taiwan Strait

Yu-Chia Chang

1, 2

, Ruo-Shan Tseng

1, 2, *

, and Cho-Teng Liu

3 1

Department of Marine Biotechnology and Resources, National Sun Yat-Sen University, Kaohsiung 804, Taiwan, ROC 2

Kuroshio Research Group, Asia-Pacific Ocean Research Center, National Sun Yat-Sen University, Kaohsiung 804, Taiwan, ROC 3Institute of Oceanography, National Taiwan University, Taipei 106, Taiwan, ROC

Received 25 September 2006, accepted 16 August 2007

ABSTRACT

Three cruises with shipboard Acoustic Doppler Current Profiler(ADCP) were performed along a transect across the Peng-hu Channel (PHC) in the Taiwan Strait during 2003 - 2004 in order to investigate the feasibility and accuracy of the phase-averaging method to eliminate tidal components from shipboard ADCP measurement of currents. In each cruise measurement was repeated a number of times along the transect with a specified time lag of either 5, 6.21, or 8 hr, and the repeated data at the same location were averaged to eliminate the tidal currents; this is the so-called “phase-averaging method”. We employed 5-phase-averaging, 4-phase-averaging, 3-phase-averaging, and 2-phase-averaging methods in this study. The residual currents and volume transport of the PHC derived from various phase-averaging methods were intercompared and were also compared with results of the least-square harmonic reduction method proposed by Simpson et al. (1990) and the least-square interpolation method using Gaussian function (Wang et al. 2004). The estimated uncertainty of the residual flow through the PHC derived from the 5-phase-averaging, 4-phase-averaging, 3-phase-averaging, and 2-phase-averaging methods is 0.3, 0.3, 1.3, and 4.6 cm s-1, respectively. Procedures for choosing a best phase average method to remove tidal currents in any particular region are also suggested.

Key words: Tidal current, Phase-averaging method, Shipboard ADCP, Residual flow

Citation: Chang, Y. C., R. S. Tseng, and C. T. Liu, 2008: Evaluation of tidal removal method using phase average technique from ADCP surveys along the Peng-Hu Channel in the Taiwan Strait. Terr. Atmos. Ocean. Sci., 19, 433-443, doi: 10.3319/TAO.2008.19.4.433(Oc)

1. INTRODUCTION

Shipboard Acoustic Doppler Current Profiler(ADCP) has been used extensively to measure three-dimensional flow structures and their vertical profiles in oceans for over twenty years. If one is interested only in the circulation of residual flow, tidal components have to be eliminated from total flow quantity as measured by ADCP. Generally speak-ing there are three basic ways of dealspeak-ing with this problem. The first approach is based on using repeated surveys along a given transect to construct time series of currents at se-lected sites, so that harmonic analysis may be performed on the data (Simpson et al. 1990). This kind of ADCP measure-ment is capable of providing a synoptic map of the tidal and

residual current field that is not obtainable by conventional means. However, the main limitation of this method is relatively small spatial coverage, which does usually not exceed 30 km in length (Geyer and Signell 1990; Vennell 1994). The second approach is to fit specific functions to the ADCP data to determine the spatially varying am-plitudes and phases of major tidal constituents and the spatially varying magnitude of the residual flow. Using polynomial or biharmonic functions, this method has been applied successfully in the Yellow Sea (Candela et al. 1992) and the southwest coast of Vancouver Island (Fore-man and Freeland 1991). A similar method using Gaussian interpolation functions with local support was proposed by Wang et al. (2004), and they were able to reconstruct the spatial and temporal flow field variability and estimate the

* Corresponding author

(2)

volume transport through the Peng-hu Channel (PHC), the gate of the Taiwan Strait. The third approach relies on designing the sampling survey in a way that permits mea-suring at the same location several times within a tidal cycle in order to safely average out the tide. This method is sometimes termed as the “phase-averaging method” be-cause the repeated measurements are designed to be at some specified phases of the semidiurnal or diurnal tidal cycle. Katoh et al. (1996, 2000) first applied this method to remove diurnal and semidiurnal tidal flows from observed flows by performing four round-trip ADCP surveys along each transect during a diurnal tidal period of 24 hr and 50 min. This method is referred to as “4R-phase-averaging method” in this study. Note that “4R” stands for four round-trip surveys, and residual flows at any position along the transect can be derived by averaging out of the eight re-peated ADCP data. A similar but simpler approach was used in the Taiwan Strait to eliminate only dominant semidiurnal tidal currents by making two repeated current measurements along the same transect with the time sepa-ration being 6 hr and 12 min (Liu et al. 2000; Jan and Chao 2003). It is termed the “2-phase-averaging method” in this study. In cases of limited availability of ship time, the phase-averaging method is fairly convenient and easy to use. In addition to the “4R-phase-averaging” and “2-averaging” methods, some other schemes of phase-averaging methods are also proposed and implemented in this study. Details of these phase-averaging methods will be described in a later section.

The objective of this study is to investigate the fea-sibility and accuracy of the phase-averaging method to eliminate tidal components from shipboard ADCP measure-ment of currents. Several schemes of phase averaging such as 2-phase, 3-phase, 4-phase, and 5-phase will be employed in this study and their results will be intercompared and also compared with two other techniques of tidal removal, i.e., the methods following Simpson et al. (1990) and Wang et al. (2004).

This paper is organized as follows. In the second sec-tion, the geography and general tidal characteristics of the region are described. In section 3, the various phase-averag-ing methods and their workphase-averag-ing principles are presented. The observations of shipboard ADCP and data reduction of the three cruises are briefly described in section 4. In section 5, uncertainties of phase-averaging methods are estimated, and various tidal removal techniques are compared using ADCP data collected in the PHC. Finally, a discussion and conclu-sion are presented.

2. GEOGRAPHY, TIDES, AND HYDROGRAPHY The Peng-hu Channel (PHC) is a funnel-shaped deep passage separating the Taiwan Island and the Peng-hu Ar-chipelago (Fig. 1). To the north, the PHC is about 40 km

wide and 100 m deep and is connected to the Chang-yuen Rise, which is a shallow seamount of roughly 40 m in depth extending from the western coast of central Taiwan and dividing the Taiwan Strait into its northern and southern parts. The southern portion of the PHC becomes wider and deeper further southward along the channel where it is con-nected to the continental shelves of the northern South China Sea (SCS). To the west of the southern PHC lies the Taiwan Banks with its depth ranging between 15 and 40 m.

Flow in the PHC is largely controlled by tides, migration of water masses (Kuroshio, South China Sea water, and China Coastal Current) and wind-driven currents. Early ob-servations from moored current meters at 20 m above the bottom from March to July (Chuang 1985, 1986) have indi-cated that a persistent northward residual flow exists in the PHC. Chuang found the northward flow in the PHC has a mean value of 0.27 cm s-1and a tidal amplitude of roughly 1.0 m s-1. Lagrangian observations of surface flow patterns from two satellite-tracked drifters also provided further evi-dence of the northward flow in the PHC during October (Tseng and Shen 2003). Analysis of shipboard ADCP data in the PHC revealed that the currents there are essentially barotropic and characterized by strong semidiurnal tides and mean flow (Wang et al. 2004). They found the along-chan-nel averaged tidal velocity amplitudes are 1.2 m s-1 and 0.32 m s-1 for semidiurnal and diurnal tides, respectively, and the mean current is 0.73 m s-1northward.

The origin of the northward flow in the PHC has been a research topic for many earlier studies. Wang and Chern (1988) suggested that the northward flow is the major inflow

Fig. 1. Map of the Taiwan Strait showing the ship tracks of velocity measurement for the three cruises in the PHC, the position of moored ADCP (marked as M), and wind station of Dongee Island (marked as DG).

(3)

entering the Taiwan Strait based on hydrographic data. Chen (2003) indicated that waters across the Taiwan Strait all flow northward in summer, and the waters originate from the South China Sea or the Kuroshio. In winter, on the other hand, strong and steady northeastern monsoon push the fresh, cold, nutrient-rich China coastal water to flow south-ward along the western part of the Taiwan Strait. On the eastern part of the Strait, including the PHC, some salty and warm, but nutrient-poor SCS or Kuroshio water flows north-ward. From analysis of eight cruises of shipboard current measurements along the PHC, Jan and Chao (2003) found that the through flow transports in the PHC vary seasonally, being minimal (around 0) during the peak winter monsoon and increase to around 1.5 Sv to the north as the southwest monsoon prevails.

3. PHASE-AVERAGING METHODS

In order to eliminate tidal currents from observed cur-rents, four types of phase-averaging methods are used in this study and their results are inter-compared. These methods are called “2-phase averaging”, “3-phase averaging”, “4 (or 4R)-phase averaging”, and “5-phase averaging”, and their schematic diagrams are shown in Fig. 2. The 2-phase aver-aging is designed to remove the semidiurnal (M2) tidal flows, and the two repeats (as the survey 1 and 3 in Fig. 2a) are separated by 6.21 hours so that they are out of phase by half a cycle of the semidiurnal tide (Jan and Chao 2003). The 3-phase averaging is designed to remove both the

semi-diurnal and semi-diurnal tidal flows. Current measurements along the same transect are repeated three times (as the survey 1, 3, and 5 in Fig. 2b), and the time spacing between survey 1 and 3 (or between survey 3 and 5) is roughly 8 hours. Theoreti-cally, the average of these three repeats will safely eliminate the semidiurnal and diurnal tidal flows. Based on a similar principle, current measurements along the same transect are repeated four times in the 4-phase averaging, i.e., the survey 1, 3, 5, and 7 as illustrated in Fig. 2c. The neighboring re-peats are separated by roughly 6.21 hours so that the average of these four surveys will eliminate both the semidiurnal and diurnal tidal flows. Note that the 4-phase averaging method is simply the extension of the 2-phase averaging method so that the number of surveys and the total surveying time of the former are exactly twice those of the latter. An adaptation of 4-phase averaging is the 4R-phase averaging in which four round-trip ADCP surveys are performed along the same transect with uniform time spacing during a diurnal tidal pe-riod of 24.84 hour (Katoh et al. 2000). That is, each survey along the transect should take exactly 3.105 hr. The average of all eight surveys (as the survey 1, 2, 3, 4, 5, 6, 7, and 8 in Fig. 2c) will eliminate both the semidiurnal (in period of 12.42 hour) and diurnal (in period of 24.84 hour) tidal flows. Finally, in 5-phase averaging the five repeated surveys (as the survey 1, 3, 5, 7, and 9 in Fig. 2d) are separated by a time spacing of approximately 5 hours. The average of these five repeats will eliminate both the semidiurnal and diurnal tidal flows. It is worth mentioning that neither the 3-phase aver-aging nor 5-phase averaver-aging methods has ever been proposed

Fig. 2. Various phase-averaging methods used in this study. They are: (a) 2-phase averaging, (b) 3-phase averaging, (c) 4 (and 4R)-phase averaging, and (d) 5-phase averaging. See text for more details.

(a)

(b)

(4)

or used before; we will test their feasibility in this study. In order to examine the accuracy of various phase aver-aging methods and to estimate errors induced by the imper-fectness of these methods, we have conducted numerical computations which applied the phase averaging methods to a composite of tidal currents. The estimation bias, denoted by VEB, which represents the error of residual flows from

phase averaging methods can be estimated from the follow-ing mathematical expression:

(1)

where Cq,wq, andFqare the amplitude, angular frequency

and phase of the qth constituent tidal current, respectively, the subscript q represents the qth constituent; M is the num-ber of tidal constituents used in the analysis; t is the time and

Dt is the time interval between neighboring surveys; N

repre-sents the total number of surveys to be averaged (N = 2, 3, 4, 5 for 2-phase, 3-phase, 4-phase, and 5-phase averaging, respectively); and the subscript i represents the ith repeated survey. Classification of tidal types is based on the conven-tional usage of “form ratio” F = (K1+ O1) / (M2+ S2) of the

sum of the amplitudes of the two main diurnal constituents to that of the two main semi-diurnal amplitudes (Pond and Pickard 1983). In our calculation we assumed that the ampli-tude of M2and S2(also for K1and O1) is equal, and the sum

of M2, S2, K1, and O1is 100 cm s-1. Four different tidal types

are considered here, that is, F = 0 (semi-diurnal tides, K1= O1

= 0, M2= S2= 50 cm s-1), F = 1, (mixed, mainly semi-diurnal

tides, K1= O1= M2= S2= 25 cm s-1), F = 2 (mixed, mainly

di-urnal tides, K1= O1= 33 cm s-1, M2= S2= 16.5 cm s-1), and F =

3 (diurnal tides, K1= O1= 37.5 cm s-1, M2= S2= 12.5 cm s-1).

Considering that the 2-phase averaging is applied to a semi-diurnal tide (F = 0), VEBcan be described as:

(2)

where t2= t1+Dt , w1,w2,F1, andF2are the angular

frequ-ency and phase of M2and S2, respectively. Since, cos A + cos B = 2

2 2

cos(A+B)cos(A B- ), Eq. (2) can be reorganized as:

(3)

For the case ofw1Dt = p, or Dt = T1/2 = 6.21 hr, Eq. (3) can be

further simplified as: VEB =C2 2 1t + t -2 2 cos (w w D F2) cos (w2 ) 2 Dt

. This means thatVEBis a function of C2,F2, and Dt. Generally, VEBof the 3-phase, 4-phase, and 5-phase

averaging are also functions of the harmonic constants of the four tidal constituents (C1, C2, C3, C4,F1,F2,F3,F4)

andDt. Our calculations of VEBwere performed for various Dt ranging between 0 and 24 hr with a time increment of 1

min over a total period of 30 days. The averaged VEBover

the 30-day period for anyDt is almost invariant with Fq,

and is plotted in Figs. 3a - d againstDt for 2-phase, 3-phase, 4-phase, and 5-phase averaging, respectively. The VEB

de-rived from various phase averaging methods has minimum values at two optimumDt for cases of different tidal form ratios (see Table 1). Several conclusions can be drawn from Fig. 3. The 2-phase averaging will yield a minimum

VEBof about 0.9 cm s-1at roughly 6.1 hr for the case of F =

0. For the other cases of F > 0, the 2-phase averaging appar-ently in not a good tide removal method because VEBis

al-ways over 10 cm s-1, which is a significant portion of the total tidal amplitude (Fig. 3a). The 3-phase, 4-phase, and 5-phase averaging methods all produce good detiding ef-fect for all tidal types of semi-diurnal, mixed and diurnal tides at aDt of about 8, 6.1, and 5 hr, respectively (the arrows in Figs. 3b - d). Their corresponding VEBat the first optimum Dt are about 1.8, 1.5, and 1.1 cm s-1

(see Table 1). Note that this optimumDt may shift slightly for different combina-tions of main tidal constituents. A deviation of±0.1 hr from the optimumDt for each phase averaging method and form ratio will increase slightly the VEB. In addition to the firstDt

which produces a minimum VEB, a second optimumDt also

exist for various phase averaging methods, e.g., 18.3, 16.3, 18.3, and 9.8 hr for 2-phase, 3-phase, 4-phase, and 5-phase averaging, respectively, and their corresponding

VEB are about 2.8, 3.2, 2.8, and 2.5 cm s-1 (Table 1).

Therefore, computation results as shown in Fig. 3 and Table 1 provide a theoretical and operational basis to de-termine theDt for various phase averaging methods.

4. OBSERVATIONS AND DATA ANALYSIS

Three field surveys were carried out along zonal tran-sects across the PHC during 2003 ~ 2004. Flow velocities were measured by using the shipboard ADCP manufactured by RD Instruments with 150 KHz for R/V OR1 or OR3. The bin depth was set to 4 m, and the averaging time interval was 1 min. The ship’s absolute movement was determined with bottom tracking, and the cruise track is shown in Fig. 1. Situated at the main section of the PHC, the transects of these three cruises are somewhat different but are parallel to each other and along approximately 23°30’N with the length of the transect about 30 km. Sampling spanned the period 15 - 18 September 2003, 29 - 31 March 2004, and 27 - 30

(5)

July 2004 for three cruises. Each cruise consisted of a num-ber of repeated transects with theirDt designed specifically to remove the tidal signals by using various phase-averaging methods.

Table 2 lists the observational period, detiding methods that were used, and other relevant information of three cruises in this study. In the first cruise of September 2003, ten round-trip surveys were made along the same PHC transect, and they consisted of 3-phase averaging, 5-phase averaging, and 2-phase averaging. In the second cruise of March 2004,

four round-trip surveys were made along the transect to test the 2-phase averaging, 4-phase averaging, and 4R-phase av-eraging methods. Note that the 4-phase and 4R-phase averag-ing is different in that the former is based on the data average of four repeated surveys while the latter is based on the data average of four round-trip (or eight) surveys. Altogether, six sets of 2-phase averaged results and one set of 4R-phase aver-aged results can be obtained from the westbound and east-bound surveys of this cruise. In the third cruise of July 2004, eleven round-trip surveys were made along the transect,

Table 2. Shipboard ADCP measurement periods, detiding methods, and transports for the three cruises in this study.

Measurement periods Methods round-trip surveys R/V Transports (Sv)

2003/09/15 ~ 2003/09/18 2, 3, and 5 phase, S90 10 OR-1 1.52 ~ 1.87 2004/03/29 ~ 2004/03/31 2 and 4 phase, S90 04 OR-3 1.19 ~ 1.70 2004/07/27 ~ 2004/07/30 3 and 5 phase, S90, W04 11 OR-3 1.22 ~ 1.58

The positive transport is toward the north. (1 Sv = 106m3s-1)

Fig. 3. The mean VEBversus theDt for: (a) 2-phase, (b) 3-phase, (c) 4-phase, and (d) 5-phase averaging methods in different cases of tidal form ratio (F = 0 ~ 3).

Table 1. The mean residual of tidal current (RT, cm s-1) from: 2-phase, 3-phase, 4-phase, and 5-phase averaging methods for different cases of tidal form ratio (F = 0 ~ 3). TwoDt in the parenthesis produce two minimum VEB.

Method \ Form ratio F = 0 F = 1 F = 2 F = 3

2 phase (Dt = 6.1, 18.3 hr) 0.9 2.8 * * * * * *

3 phase (Dt = 8.1, 16.3 hr) 2.0 3.2 1.7 2.3 1.8 2.5 1.8 2.7 4 phase (Dt = 6.1, 18.3 hr) 1.4 2.8 1.5 2.5 1.4 3.1 1.6 3.3 5 phase (Dt = 4.9, 9.8 hr) 1.0 2.5 1.1 1.8 1.1 1.7 1.3 1.7

(6)

which consisted of 5-phase averaging and 3-phase averag-ing surveys. Four sets of 5-phase averaged results and six sets of 3-phase averaged results can be obtained from the westbound and eastbound surveys of this cruise.

Although the observational strategy of the three cruises was originally developed to test the various phase-averaging methods, the repeated surveys along a transect permit the ap-plication of stationary analysis techniques (Simpson et al. 1990) or the least squares interpolation techniques using Gaussian interpolation functions with local support (Wang et al. 2004). The former method is abbreviated as S90, and the latter is W04 in this study. These two methods will be used where appropriate to obtain the residual currents at each segment along the transect, and the results of the mean flow profiles and calculated volume transport can be com-pared with those of the phase-averaging methods.

The first bin of the ADCP data is 10 and 10.6 m for R/V OR1 and OR3, respectively. Velocities of the surface layer from the first bin up to the sea surface are assumed to be con-stant and equal to the velocity of the first bin. The deepest bin of the ADCP data is determined by the following criteria of data screening: the acceptable data are characterized by percentage with good being greater than 85% and the bin depth being smaller than 85% of the mean water depth. Ve-locity of the bottom layer is obtained from linear interpola-tion between the last-bin velocity and the zero bottom veloc-ity. The velocity data along the transect are divided into 10 segments and averaged over each segment.

The averaged velocity is assigned at the midpoint of each segment. For the phase-averaging method, the velocity of each segment and bin for the corresponding tracks were averaged to get the mean velocity. The volume transport through the transect is calculated using the mean velocity perpendicular to the transect as:

(4)

where n is the segment number; k is the layer index; Kiis

the bottom layer index of the ith segment; v is the v-compo-nent residual velocity; d is the width of each segment; and Dz is the layer thickness.

Hourly wind data, which were collected at an islet weather station, Dongee, near the south of the PHC transect (DG in Fig. 1), were acquired from the Central Weather Bu-reau of Taiwan for further data analysis.

5. RESULTS

5.1 Estimate of Errors for Phase-Averaging Methods in the PHC

For the four phase-averaging methods used in this study,

the 2-phase averaging is the semidiurnal averaging, and the other three methods (3-phase, 4-phase, and 5-phase) are the diurnal averaging. Therefore, it is already clear to conclude that there are some portions of diurnal tidal currents which cannot be eliminated by the 2-phase averaging, and the con-sequent error of mean velocity and volume transport of the 2-phase averaging will be larger than the other diurnal aver-aging methods. We estimate the errors of various phase-averaging methods based on a dataset of currents collected near the PHC as follow.

Theoretically, Eq. (1) is only applicable for any fixed location, and the harmonic constants of Cq and Fq are

location-dependent. Practically, it is quite safe to assume that CqandFqare invariant along a transect line or studied

region and can be obtained from a nearby moored current re-cord. The current data was measured by a bottom-mounted ADCP which was deployed at a station about 90 km south of the PHC (M in Fig. 1). The water depth of the station is approximately 300 m, the bin depth of the ADCP was set to 10 m, and the period of observation is from 22 June to 19 November 1996. Harmonic analysis suggests that tidal cur-rent amplitude of M2, O1, K1, S2, and N2for the dominant v-component are 24, 8, 6, 6, and 4 cm s-1, respectively, and the corresponding phases for the five constituents are 243°, 87°, 137°, 282°, and 228°. The blue line in Fig. 4 shows the composite tidal currents over a period of 30 days, which were reconstructed from the amplitude and phase of the five major constituents. The VEB for various phase-averaging

methods can be estimated from Eq. (1), in which theDt is set to be 6.21, 8, 6.21, and 5 hr, respectively for 2-phase, 3-phase, 4-phase, and 5-phase averaging. The green, red, cyan and purple lines in Fig. 4 show VEBversus time for

2-phase, 3-phase, 4-phase, and 5-phase averaging,

respec-Fig. 4. Time series of v component reconstructed from five tidal con-stituents (O1+ K1+ M2+ S2+ N2) at station M in the PHC, and VEBafter applying various phase-averaging methods.

(7)

tively. Note that the VEBin the PHC from the 2-phase

averag-ing ranges from about±2.5 cm s-1at the neap tide to about ±10 cm s-1

at the spring tide. On the other hand, the VEBin the

PHC from the 4-phase and 5-phase averaging remain quite small (less than 2 cm s-1) irrespective of the lunar time. For the 3-phase averaging, the VEB ranges between -2.5 and

2.5 cm s-1. The mean of absolute VEBis 4.6, 1.3, 0.3, and

0.3 cm s-1, respectively for 2-phase, 3-phase, 4-phase, and 5-phase averaging (Fig. 4). The consequent uncertainty of volume transport through the PHC is 0.13, 0.04, 0.01, and 0.01 Sv. It can be concluded that for the PHC region the error of the VEBand the consequent error of volume transport

in-duced by the semi-diurnally 2-phase averaging method are significantly larger than those induced by the diurnally aver-aging methods of 3-phase, 4-phase, and 5-phase averaver-aging. Among the three diurnally averaging methods used in this study, 4-phase and 5-phase averaging yield slightly better re-sults than 3-phase averaging.

5.2 Velocity Profiles and Volume Transports

Now we turn our attention to the phase-averaged results of the shipboard ADCP data from three cruises. In the first experiment from September 15 to 18 of 2003, three phase-averaging methods were conducted in the PHC sequentially, i.e., in the order of 3-phase averaging, 5-phase averaging, and 2-phase averaging. That is, a total of ten repeated sur-veys were done along the transect. For each time span the corresponding phase-averaging analysis was applied to the

ADCP data to obtain the residual velocity at each segment and vertical bin. In addition, the stationary analysis tech-nique of S90 was also applied to the twenty repeated ob-servations at each location to obtain the mean velocity. A comparison of the 5-phase averaged, 3-phase averaged, 2-phase averaged, and S90-derived velocity contours of the v-component in the transect is provided in Fig. 5. The flow velocities are essentially northward in transport during this cruise. Zonal currents (u-component) are much weaker and therefore are not shown here. It can be seen from the velocity contour plots that a core structure centered in the middle of the transect was developed. The northward velocity in the upper 40 m layer can reach over 80 cm s-1in the 3-phase av-eraged, 5-phase avav-eraged, and S90-derived velocity con-tours. The reason for the apparent similarity between Figs. 5a, b, and d can be attributed to their common nature of diurnally averaging analysis. On the other hand, velocity contours from the 2-phase averaging, which eliminate only the semidiurnal tidal currents, are of somewhat different shape (Fig. 5c) such as having a higher core velocity of over 1 m s-1. Also note that the S90 analysis was made based on the least-square interpo-lation over the entire time span of this experiment, rather than the individual time span of each phase-averaging analysis. As a consequence, the velocity contours derived from the S90 analysis (Fig. 5d) are smoother and in a sense as a combina-tion of the other three contours.

In the second experiment from March 29 to 31 of 2004, four round-trip surveys along the PHC transect were con-ducted to test the 4R-phase averaging method. Note that the

Fig. 5. Velocity contours (m s-1) of PHC from: (a) 3-phase, (b) 5-phase, (c) 2-phase averaging, and (d) S90 analysis (Simpson et al. 1990) in Septem-ber, 2003.

(a) 3 phase (b) 5 phase

(8)

time span of each survey (either eastbound or westbound) along the transect was maintained constant (3.06 hr). There-fore, six sets of individual 2-phase averaged results can also be obtained from this experiment in addition to the 4R-phase averaged result. Figure 6 shows the contours of mean veloc-ity derived from the six 2-phase averaging analysis, the 4R-phase averaging analysis, and the stationary analysis technique of S90. All contours have similar pattern, with the maximum core velocity reaching around 0.6 m s-1near the central, upper layer of the PHC. The diurnally-averaged

velocity contours of 4R-phase averaging have the most resemblance to those of the S90 analysis, which was also fit-ted to both M2and K1functions.

In the third experiment from July 27 to 30 of 2004, eleven round-trip surveys were conducted along the PHC transect. The first six round-trips were designed to run in 5-phase averaging scheme, and the next five round-trips were designed to run in 3-phase averaging scheme. The extra surveys, more than that required for each phase-averaging method, allow us to obtain more results to make

Fig. 6. Velocity contours (m s-1) of PHC from: (a) 4R phase, (b) 21phase, (c) 22phase, (d) 23phase, (e) 24phase, (f) 25phase, (g) 26phase averaging,

and (h) S90 (Simpson et al. 1990) in March, 2004.

(a) 4R phase (b) 21phase

(c) 22phase (d) 23phase

(e) 24phase (f) 25phase

(9)

comparison. In addition to the phase-averaging method, the stationary analysis technique of S90 and the least-square interpolation technique using Gaussian functions of W04 were also applied to the twenty-two observations at each location to obtain the mean velocity. Figure 7 shows the

mean velocity contours of PHC for each phase-averaging run and from the technique of S90 and W04. All contours show a maximum speed of about 1 m s-1at the central core region of the PHC.

Because the transects of these three cruises are not

Fig. 7. Velocity contours (m s-1) of PHC from: (a) 51phase, (b) 52phase, (c) 53phase, (d) 54phase, (e) 31phase, (f) 32phase, (g) 33phase, (h) 34phase,

(i) 35phase, (j) 36phase averaging, (k) S90 (Simpson et al. 1990), and (l) W04 (Wang et al. 2004) in July, 2004.

(a) (b) (c)

(d) (e) (f)

(g) (h) (i)

(10)

completely identical, having a maximum separation of about 20 km apart in between the transects (Fig. 1), it will only be qualitative to compare the detided results of these three cruises. The maximum velocities of the central core in the PHC transect for these three cruises are in accord with the wind speed. For instance, the maximum core velocities are over 100, 80, and 60 cm s-1flowing toward the north re-spectively for the July 2004, September 2003, and March 2004 cruises during which the local winds are approxi-mately 6 m s-1south wind, 1 m s-1north wind, and 5 m s-1 north wind. This finding was confirmed previously by Jan and Chao (2003) and Wang et al. (2004).

6. DISCUSSION

The phase average method can be applied fairly conve-niently to either smaller or larger transects or spatial cover-age. For shorter channels such as PHC with a length of no more than about 50 km, 2-phase, 3-phase, 4-phase, or even 5-phase averaging is all feasible with a maximum ship speed of 12 knots. In this case the detiding methods of S90 and W04 are also applicable. For larger study regions and longer transects, S90 and W04 may not be useful because usually at least more than 8 observation at any fixed location are needed for least-square interpolation (Geyer and Signell 1990). The phase average method is, however, still good to use in this case. To choose an appropriate phase averaging scheme for any study region, we suggest the following

pro-cedures contained here in this example for the region of Taiwan Strait. The length of a transect across the Taiwan Strait is about 150 km, and the nominal ship speed is 10 knots. It will take approximately 8 hours for the ship to finish this transect. The dominant tidal current amplitude for the Taiwan Strait is 31 cm s-1(M2), 9 cm s-1(S2), 9 cm s-1(K1), and 7 cm s-1(O1) according to the moored current data re-ported by Lin et al. (2005). The mean VEBfor various phase

average methods as a function of Dt can then be obtained from Eq. (1) and is plotted in Figs. 8a - d for 2-phase, 3-phase, 4-phase, and 5-phase averaging, respectively. From Fig. 8, some optimumDt and the corresponding VEB

of some phase averaging schemes, which satisfy the condi-tion ofDt ³ 2T (16 hours), where T is the time to finish one transect across the Taiwan Strait and their total required ship time are listed below:

(i) 2-phase:Dt = 18.4 hr, VEB= 5 cm s-1, total required time

= 2Dt = 36.8 hr

(ii) 3-phase:Dt = 16.5 hr, VEB= 1.5 cm s-1, total required

time = 3Dt = 49.5 hr

(iii) 4-phase:Dt = 18.4 hr, VEB= 1.5 cm s-1, total required

time = 4Dt = 73.6 hr

(iv) 5-phase:Dt = 17.3 hr, VEB= 2 cm s-1, total required time

= 5Dt = 86.5 hr, or Dt = 19.9 hr, VEB= 1.5 cm s-1, total

re-quired time = 5Dt = 99.5 hr

A trade-off is usually to be made between the phase

Fig. 8. The mean VEBwith respect to theDt for: (a) 2-phase, (b) 3-phase, (c) 4-phase, and (d) 5-phase averaging methods in the Taiwan Strait based on the reconstructed tidal currents (Lin et al. 2005).

(11)

averaging scheme (and thus the total required ship time) and the VEB. For this example, in the Taiwan Strait, the

3-phase averaging appears to be a fairly reasonable and fea-sible detiding method because total required ship time is less while still obtaining good accuracy. The errors of the re-sidual velocity and the net volume transport induced by the 3-phase averaging are±1.5 cm s-1and±0.14 Sv in the Tai-wan Strait. One of our earlier cruises of August 2002 was conducted using the 3-phase averaging method in the Tai-wan Strait and the results of mean flows and transport will be reported elsewhere.

Acknowledgements This work was supported by the National Science Council of the Republic of China under contract NSC92-2611-M-110-006. Prof. Yu-Huai Wang kindly provided us the Matlab programs for least-square interpolation using Gaussian function. We also would like to thank Prof. David T. Y. Tang for providing us the bottom-mounted ADCP data. This paper represents Con-tribution number 12 of Kuroshio Research Group, Asia-Pacific Ocean Research Center.

REFERENCES

Candela, J., R. C. Beardsley, and R. Limeburner, 1992: Separa-tion of tidal and subtidal current in shipboard-mounted acoustic Doppler current profiler observations. J.

Geo-phys. Res., 97, 769-788.

Chen, C. T. A., 2003: Rare northward flow in the Taiwan Strait in winter: A note. Cont. Shelf Res. 23, 387-391.

Chuang, W. S., 1985: Dynamics of subtidal flow in the Taiwan Strait. J. Oceanogr. Soc. Jpn., 41, 65-72.

Chuang, W. S., 1986: A note on the driving mechanisms of cur-rent in the Taiwan Strait. J. Oceanogr. Soc. Jpn., 42, 355-361.

Foreman, M. G. G. and H. J. Freeland, 1991: A comparison of techniques for tide removal from ship-mounted acoustic

Doppler measurements along the southwest coast of Vancouver Island. J. Geophys. Res., 96, 17007-17021. Geyer, W. R. and R. Signell, 1990: Measurements of tidal flow

around a headland with a shipboard acoustic Doppler cur-rent profiler. J. Geophys. Res., 95, 3189-3197.

Jan, S. and S. Y. Chao, 2003: Seasonal variation of volume transport in the major inflow region of the Taiwan Strait: The Penghu Channel. Deep-Sea Res. II, 50, 1117-1126. Katoh, O., K. Teshima, O. Abe, H. Fujita, K. Miyaji, K.

Morinaga, and N. Nakagawa, 1996: Process of the Tsu-shima Current formation revealed by ADCP measure-ments in summer. J. Oceanogr., 52, 491-507.

Katoh, O., K. Mornaga, and N. Nakagawa, 2000: Current distri-butions in the southern East China Sea in summer. J.

Geophys. Res., 105, 8565-8573.

Lin, S. F., T. Y. Tang, S. Jan, and C. J. Chen, 2005: Taiwan strait current in winter. Cont. Shelf Res., 25, 1023-1042. Liu, K. K., T. Y. Tang, G. C. Gong, L. Y Chen, and F. K. Shiah,

2000: Cross-shelf and along-shelf nutrient fluxes derived from flow fields and chemical hydrography observed in the southern East China Sea off northern Taiwan. Cont.

Shelf Res., 20, 493-523.

Simpson, J. H., E. G. Mitchelson-Jacob, and A. E. Hill., 1990: Flow structure in a channel from acoustic Doppler current profiler. Cont. Shelf Res., 10, 589-603.

Tseng, R. S. and Y. T. Shen, 2003: Lagrangian observations of surface flow patterns in the vicinity of Taiwan. Deep-Sea

Res. II, 50, 1107-1115.

Vennell, R., 1994: Acoustic Doppler current profiler measure-ments of tidal phase and amplitude in Cook Strait, New Zealand. Cont. Shelf Res., 14, 353-364.

Wang, J. and C. S. Chern, 1988: On the Kurshio branch in the Taiwan strait during wintertime. Prog. Oceanogr. , 21, 469-491.

Wang, Y. H., L. Y. Chiao, K. M. M. Lwiza, and D. P. Wang, 2004: Analysis of flow at the gate of Taiwan Strait. J. Geophys.

數據

Fig. 1. Map of the Taiwan Strait showing the ship tracks of velocity measurement for the three cruises in the PHC, the position of moored ADCP (marked as M), and wind station of Dongee Island (marked as DG).
Fig. 2. Various phase-averaging methods used in this study. They are: (a) 2-phase averaging, (b) 3-phase averaging, (c) 4 (and 4R)-phase averaging, and (d) 5-phase averaging
Table 2. Shipboard ADCP measurement periods, detiding methods, and transports for the three cruises in this study.
Fig. 4. Time series of v component reconstructed from five tidal con- con-stituents (O 1 + K 1 + M 2 + S 2 + N 2 ) at station M in the PHC, and V EB after applying various phase-averaging methods.
+5

參考文獻

相關文件

A mixture-energy-consistent 6-equation two-phase numerical model for fluids with interfaces, cavitation and evaporation waves. Modelling phase transition in metastable

Attributes characterized in this campaign for the turboshaft engine emissions include, for particulate phase, engine smoke number (SN) (which currently is the only accepted method

The main tool in our reconstruction method is the complex geometri- cal optics (CGO) solutions with polynomial-type phase functions for the Helmholtz equation.. This type of

Recommendation 14: Subject to the availability of resources and the proposed parameters, we recommend that the Government should consider extending the Financial Assistance

• Non-vanishing Berry phase results from a non-analyticity in the electronic wave function as function of R.. • Non-vanishing Berry phase results from a non-analyticity in

Let T ⇤ be the temperature at which the GWs are produced from the cosmological phase transition. Without significant reheating, this temperature can be approximated by the

◆ Understand the time evolutions of the matrix model to reveal the time evolution of string/gravity. ◆ Study the GGE and consider the application to string and

For ASTROD-GW arm length of 260 Gm (1.73 AU) the weak-light phase locking requirement is for 100 fW laser light to lock with an onboard laser oscillator. • Weak-light phase