Warning: mkdir(): Permission denied in /home/virtual/lib/view_data.php on line 81

Warning: fopen(upload/ip_log/ip_log_2024-03.txt): failed to open stream: No such file or directory in /home/virtual/lib/view_data.php on line 83

Warning: fwrite() expects parameter 1 to be resource, boolean given in /home/virtual/lib/view_data.php on line 84
Characteristics of Tidal Current and Tidal Residual Current in the Chunsu Bay, Yellow Sea, Korea based on Numerical Modeling Experiments

Characteristics of Tidal Current and Tidal Residual Current in the Chunsu Bay, Yellow Sea, Korea based on Numerical Modeling Experiments

수치모델링 실험을 통한 서해 천수만의 조류와 조석잔차류 특성

Article information

J Korean Soc Coast Ocean Eng. 2013;25(4):207-218
*Young Jae Ro, Oceanography and Ocean Environmental Sciences, College of Natural Sciences, Chungnam National University, Daejeon 305-764, Korea, Tel: +82-42-821-6437, Fax: +82-42-822-8173, royoungj@cnu.ac.kr
Received 2013 June 13; Revised 2013 August 06; Accepted 2013 August 19.

Abstract

This study is based on a series of numerical modeling experiments to understand the circulation and its change in the Chunsu Bay (CSB), Yellow Sea of Korea. A skill analysis was performed for the tidal height and tidal current of the observation data using the amplitude and phase of the 4 major tidal constituents respectively for verification of modeling experimental results. As a result, most of the skill score was seen to be over 90%, so numerical model experiment results can be said to be in good agreement with the observed tidal height and tidal current. Tidal wave proceeded from the entrance of the CSB towards inside, and the tidal range gradually increased to the north. It took about 10 to 30 minutes for the tidal wave to reach to northern end. The tidal wave showed a characteristic to rotate counter-clockwise in the southern part. The tidal current flowed to the north-south direction along the bottom topography; the angle of the major axis appeared alongside the isobath. It showed the characteristics of reversing tidal current with the minor axis less than 10% of the major axis. The strength of the tidal residual current that is influenced by geographical factors including bathymetry and coastline showed the range of 1~30 cm/sec, greater in the south channel and smaller in northern Bay. Two pairs of cyclonic/anti-cyclonic eddies around Jukdo and 3~4 pairs of strong eddies at the southern part of CSB in hundreds of m to a few km size by relative vorticity derived from the tidal residual current.

Trans Abstract

수치 모델링 실험을 활용하여 서해 천수만의 해수 유동과 그 변화를 이해하기 위한 연구를 수행했다. 모델링 실험 결과에 대한 검증을 위해 관측 자료의 조위와 조류 각각 4대 분조의 진폭과 위상을 이용하여 스킬 분석을 실시했다. 그 결과 스킬 점수는 대부분 90%가 넘는 것으로 보아 수치 모델링 실험 결과는 관측된 조위와 조류가 양호하게 일치하는 것으로 나타났다. 천수만의 조석파는 만 입구에서 안쪽으로 진행되며 북부로 갈수록 조차는 점차 증가했다. 조석파가 북부까지 도달하는데 약 10~30분의 시간이 소요되었다. 남부에서 조석파는 반시계 방향으로 회전하는 특성을 보였다. 조류는 해저 지형을 따라 남-북 방향으로 흘렀으며, 장축의 각도는 등수심선과 나란히 나타났다. 조류타원의 단축이 장축의 10% 이하로 왕복성 조류의 특성을 보였다. 수심과 해안선 등 지형적 요인에 의해 좌우되는 조석잔차류의 크기는 1~30cm/sec의 범위를 보였고, 남쪽 수로에서 컸으며 만의 북부에서는 작았다. 조석잔차류로부터 유도된 상대와도를 통해 수 백 m에서 수 km 크기로 시계/반시계 방향으로 회전하는 와류를 확인했고, 죽도 주변에서 2쌍, 남부에서 형성된 3~4쌍의 강한 와류 특성을 파악했다.

1. Introduction

The Chunsu Bay (hereafter, CSB) is a shallow back Bay with its water depth within 25 m located at 36o23'~36o37'N,126o20'~126o30'E, on the western central coast of Korea, surrounded by Anmyeondo of Taean-gun, Kanwoldo of Seosan, and the West Sea Line of Boryeong (Fig. 1). As the reclamation projects proceeded in Seosan A·B zone, the tidal embankment construction clapboard was completed in 1983~1985 and two major channels through the northern part of CSB were blocked; as a result the CSB area was reduced from 380 km2 to 180 km2 and the flow velocity near the tidal embankment has plunged down (So et al., 1998; Lee et al., 2011). This change of velocity field affected various marine environment variables; the production of fish, crustacean and mollusks was reduced by about 1/4 from the end of 1970s to the late 1990s, from the change of the marine environment and sediment due to the geographic variation according to the construction of the tidal embankment and the influx of eutrophic freshwater. In addition, it was reported recently that the summer stratification in CSB generated low-rise anoxic layer (hypoxia) and threatened the ecosystem (Lee et al., 2012).

Fig. 1.

Map for the study area with locations for the various stations in the Chunsu Bay (CSB), Yellow Sea of Korea.

There are few researches until currently, however, on the physical process and characteristics in CSB and their impact to the ecosystem. Therefore, the importance of this study are as follows: first, the physical understanding of circulation (tidal current, density-driven current, and wind-driven current) with complex coastline and bottom topography considering the freshwater inflow from Kanwolho (hereafter, KW) and Bunamho (hereafter, BN), and the impact of the wind; second, understanding of changes in the ecosystem including eutrophication and the hypoxia, imminent regional issues from these geophysical processes.

The author has performing and reported a variety of researches recently for two years as follows by utilizing a numerical model in order to understand the physical characteristics of CSB with many of these issues. Hydrodynamic and hydrographic conditions in the CSB was observed and analyzed in detail using current meter records and water quality instrument by the present author (Jung et al., 2012b). Numerical modeling studies for the tide and tidal current in the CSB was also made by Jung et al. (2011a) and impact of the fresh water release on the circulation system in the CSB was investigated by Jung et al. (2011b, 2011c). Jung et al. (2012a, 2012b) also analyzed the salinity variation and stratification processes caused by freshwater release from the KW/BN based on numerical modeling study. Add to these results, Jung et al. (2012c) studied tracking patterns of freshwater using particle trajectory modeling experiments.

In the rapidly changing coastal marine environment, the inflow and behavior of freshwater, contaminants and other conservative substances are influenced by residual currents including the tidal residual current, wind-driven current and density-driven current (Yanagi, 1983). Therefore, in-depth researches are needed on the residual current influencing the net drift of substances in CSB that faces many problems in the marine environment and ecosystems in the summer due to the influx of a large amount of eutrophic freshwater from KW/BN. In the past, a variety of studies were performed in the world associated with the circulation in the coast and estuary on the characteristics of the tide and tidal current, and residual currents including wind-driven current and density-driven current, through observations and numerical models (Imasato, 1983; Goodrich et al., 1987; Dube et al., 1995; Guo and Valle-Levinson, 2008; Zhai et al., 2008); a variety of researches are also being conducted on the mechanism for the occurrence of residual currents (Kashiwai, 1984; Guo and Yanagi, 1996; Cai et al., 2003; Foreman et al., 2006).

In Korea, the physical researches on CSB were not much in the past 10 years compared to other Bays of Korea, except for a handful of papers including the study on the numerical modeling of circulation and diffusion of CSB and inshore of Yoo (1992), tidal phenomena change due to the tidal embankment construction of So et al. (1998), and CSB water quality prediction model development of Choi (2004). In this study, model configurations such as horizontal grid size, bottom topography, tidal constituents for model forcing and locations of open boundary, etc have improved compared to previous researches (Choi, 2004; Park and Oh, 1998; Seo et al., 1998; Yoo, 1992).

In this study, series of numerical modeling experiments were carried out and the results were analyzed to understand the tidal dynamics, the density-driven current caused by the freshwater discharge from the KW/BN tidal embankment, and the wind-driven current caused by the local wind. We intend to describe the results of intensive numerical modeling experiments in three separate papers in that in Paper 1, characteristics of tide and tidal current are analyzed in terms of tidal constituents; in Paper 2, the impact of summertime freshwater discharge on the modification of tidal current and its ellipse characteristics and change of 3D structure of salinity field; in Paper 3, the wind-driven current and its influence on the destratification of water column will be described. The objective of this paper is to analyze the characteristics of the tide, tidal current, and its residual current with associated vorticity field.

2. Numerical Model for the Chunsu Bay

2.1 Physical Setting and Bottom Topography

The CSB is located in the western coast of the Korean peninsula facing the Yellow Sea. It is a semi-closed north-south elongated bay blocked by the Seosan AB tidal embankment in the north. It is 25 km long by 8.5 km wide. The CSB is shallow with average depth of 10 m and maximum water depth is approximately 30 m at the entrance of the Bay shown in Fig. 1. A large amount of freshwater in two inland lakes (i.e., KW and BN) is discharged into the CSB during the summer season. Discharge rates for 3~4 hours were 400, 200 m3 /sec respectively. The annual 780 ×106m3 and 620 ×106m3 of freshwater entered the CSB from two tidal embankments (KW/BN) in 2010 and 2011 respectively. CSB is a typical shallow back bay with the bathymetry distribution ratio of 89.2% less than 20 m in water depth and 49.3% less than 10 m. As shown in Fig. 2 CSB was separated into northern, central, and southern regions.

Fig. 2.

Grid system for the numerical model for the CSB, Yellow Sea of Korea. And location of datasets for model calibration and validation.

Jung et al. (2013a), observed the multilayer current speeds and directions at the center CSB in summer 2010 by the ADCP. The current speed was observed up to 40 cm/sec at surface, and the direction of the main flow was the NNESSW direction. The residual current flowed northward over the entire layers, with longer flow distance in the mid-layer than the surface, which seems to be due to the influence of fresh water (density-driven current) flowing south from the north of CSB. It is interpreted that the freshwater reached to the center of CSB and affected the circulation in 2~3 days from KW/BN lakes. CSB is the region with summer season monsoon climate, having prevailing southerly winds during the summer, and northwesterly winds prevailed in the winter. In addition, a large amount of fresh water flows in from KW/BN lakes in the summer monsoon season. The dynamics of the water circulation in CSB of is typically dominated by tides and tidal currents, but its impact varies depending on the effect of wind-driven current and density-driven current due to freshwater inflow. Influence of fresh water from KW and BN tide embankment in to the CSB and local wind-driven current will be of major subjects in the next papers.

2.2 Model Specification and Numerical Schemes

The numerical model for CSB estuarine used in this study is ECOM3D, a three-dimensional, sigma coordinate, hydrostatic, primitive equation model derived from Princeton Ocean Model (Blumberg and Mellor, 1987). No salt and heat flux conditions are applied to the surface boundary. At the bottom, the quadratic bottom stress with a bottom drag coefficient is applied with no salt and heat fluxes and zero vertical velocity. The costal wall boundary is impenetrable, impermeable and no-slip. A wetting and drying schemes was applied to the model, which defines dry-cells as regions with a thin water column of fluid 0 cm.

Characteristics of the model for this study are as follows. Freshwater flows in from KW/BN lakes in CSB north as in Fig. 1; seawater circulates with the offshore through the south channel. Wide-area grid was formed by extending the south and west boundaries, in order to minimize the problems that can occur when CSB entrance and the open boundary is close enough. 79 ×136 variable grids were constructed with horizontal grid size 250~4000 m for modeling; and CSB are configured with finer equal grids of 250 m for detailed understanding of the internal circulation in the Bay (Fig. 2). We used the mean sea level (MSL) data for model depth. In order to get a better understanding of the seawater circulation in the surface layer and the bottom layer, it was divided into 9 vertical sigma layers; most surface layer grids affected by the freshwater inflow were constructed with the vertical grid size of 0.5 m or less.

The initial condition of the modeling experiments was specified for the summer case with water temperature, 24.5oC, salinity, 31.0 (psu) obtained from on-site observations. Southern and western open boundary conditions were imposed by 4 major tidal constituents (M2, S2, K1 and O1) to force the tidal oscillation. Also the radiation boundary condition (Orlanski, 1976) was imposed at ebbing phase of model simulations. As shown in Table 1, the amplitude and phase of the 4 major tidal constituents were specified by using the tidal chart of NAO99 (Mastumoto et al., 2000). The amplitude of the M2 tidal constituent was given 1.98~2.16 m along the southern boundary, and 1.98~2.14 m along thewestern boundary. Model performances were evaluated based on the skill analyses in terms of tidal harmonics and tidal current ellipse parameters (not shown in this paper).

Amplitude (m) and phases (deg.) of tidal constituents used for model forcings at open boundary points

2.3 Field Measurements and Datasets for Model Calibration and Validation

Long time series of ADCP measurements at 14 levels with 1 meter vertical interval, for the period of 32 days from July 29 to Aug 30, 2010, (Jung, et al., 2012b, 2013a). Time series of measured surface current and its residual current at Jukdo station are shown in Fig. 3. In addition, a total of eight times field water quality surveys were conducted by using YSI6600, in terms of six parameters such as temperature, salinity, dissolved oxygen, chlorophyll, turbidity, and pH in 2010 and 2011, strong salinity stratification occurred in interlayer in summer; hypoxia water masses less than 3 ppm emerged in bottom layer; and dissolved oxygen oversaturation (120~160%) and chlorophyll a surge (up to 57 μg/L) phenomenon were observed in surface layer. The current and physical properties data were used as initial and boundary conditions, and validation data for the numerical modeling experiment of the study area (Fig. 1).

Fig. 3.

Time series of raw current meter records (a) and the residual current (c) at surface during Jul. 28~Aug. 30, 2010 at Jukdo Station.

Harmonic analyses were carried out to estimate tidal constituents and ellipse parameters by using the program package (Pawlowicz et al., 2002). Twenty tidal constituents were obtained by using 18-day record of one hourly sampled elevation and/or current components. The shortest and longest constituents resolved were M10 (2.48 hr) and MSF (354.37 hr) in our harmonic analyses.

3. Simulation of Tide and Tidal Current

3.1 Model Assessment and Validation

Three cases of numerical experiments were designed to understand 1) the basic tidal and tidal current characteristics in CS10TN case, 2) the impacts of the KW/BN freshwater release on the circulation system in CS10TD case and 3) added on the wind force and its effect in CS10TDW. In this paper, Part 1, only the results of Run CS10T are presented.

The verification of the model is indispensable before the modeling experiment. Among modeling experiment results, the verification of amplitude and phase of the tidal height and tidal current was performed using past observed tidal record (Jung et al., 2013a) and paper (Choi, 2004). The conditions of the open boundary were given using NAO99 data, but tidal height conditions of the open boundary was partly adjusted until the skill score (hereafter, SS) of M2, S2 tidal constituent becomes 85% or more at the stations of each validation.

To assess the model performance quantitatively, we adopt the skill score suggested by Martin and McCutcheon (1999) which is define as the deviation of the relative error from unity in the Eqn. (1).

Here X can be any of characteristic parameters such as tidal amplitude, phase lag, or tidal ellipse parameters for the tidal elevation and velocity component.

Time series of sea level at the Boryung tidal station (Korea Hydrographic and Oceanographic Administration, abbreviated as KHOA) is shown in Fig. 4. To show the model performance of elevation to local tidal records (Boryung, KHOA), Kojeong and Kanwol (Choi, 2004)), SS of amplitudes and phase lags for the 4 major tidal constituents and are shown in Table 2. They show the average relative errors of 12.9% and 3.0% for the amplitude and phase lags, respectively. The skill scores demonstrate that the model performance is satisfactory. The tidal characteristics agree with previous works by Choi (2004).

Fig. 4.

Time series of sea level elevation at Boryung Station in the CB. Elevation compared to model result (CS10T).

Comparison between model output and observations of sea level at three points in terms of 4 tidal constituents (unit: amplitude (m), phase (deg.) and skill score (%))

Time series of measured and model-computed surface current speeds at Jukdo station are shown in Fig. 5. Velocity components exhibit semi-diurnal oscillation with the spring-neap modulation in Fig. 5. The results of skill analyses of current in terms of the 4 major tidal constituents for the Jukdo station are shown in Table 3. The SS of amplitude and phase of M2 and S2 were greater than 90%. The amplitude and phase of K1 and O1 were around 70~78% and 60~90%, respectively. The diurnal constituents show higher correlation to density-driven current caused by the fresh water discharge, which will be simulated and analyzed in next paper.

Fig. 5.

Time series of observed (ADCP) and model predict at Jukdo Station in the CB (a: U-component; b: V-component).

Comparison between model output and ADCP record for tidal current at Jukdo Station in terms of 4 tidal constituents (unit: amplitude (cm/s), phase (deg.) and skill score (%))

3.2 Simulation of Tide and Tidal Current

In Fig. 4, the time-series of the sea level shows the semidiurnal oscillation with fortnight modulation. The results of modeling experiments showed the spring range of 622 cm in Boryeong located at CSB entrance, 608 cm in adjacent Kojeong, and high 650 cm in KW north of CSB about 25 km away from Boryeong. The range at neap tide showed 210, 232 and 242 cm in Boryeong, Kojeong and KW respectively. Approximate Highest High Water (AHHW) was 728 cm in Boryeong, 718 cm in Kojeong and 760 cm in KW; High Water Ordinary Spring Tide (HWOST) was 675, 663 and 705 cm, respectively. In addition, High Water Ordinary Neap Tide (HWONT) was 469, 475 and 501 cm, respectively. KW, north of CSB showed larger spring and neap range than Kojeong, at the entrance of CSB by 42 and 10 cm, respectively. The tidal range gradually increases from the entrance to the northern limit of the CSB.

Fig. 6 show the tidal chart for the M2 and S2 constituents in the CSB. The amplitude of M2 and S2 tidal constituent showed arange of 202~223 cm and 93~112 cm, respectively. They show the increasing tendency towards the north from the bay entrance. On the other hand, the phase distributions show very complicated patterns near the entrance of the CSB probably due to the small islands. The phases of the M2 and S2 tidal constituent showed arange of 99~102o and 147~157o, in the southern part of Jukdo, respectively. The time lag from the entrance to the northern end is about 10~30 minutes. Simple calculation of the propagation velocity using 10 meter average depth will take around 33 minutes which agrees well with the phase difference in Figs. 6.

Fig. 6.

Co-tidal chart for M2 (a) and S2 (b) constituents in the CSB.

Tidal current in the CSB is complex due to a various geographical and physical factors such as the complex coastline, irregular bathymetry and open boundary. The northward velocity is approximately three times bigger than the east-west component at the central CSB, which is simply due to the fact that the CSB lying along the north-south. The current speed is fast at the south channel of the Bay entrance, but shows the gradual decrease towards the inner side of the Bay.

Fig. 7 show the horizontal distribution of simulated tidal current during flood/ebb of the spring/neap tide. The maximum ebb/flood current at spring tide shows the southward flow in Fig. 7a /northward in Fig. 7b. Both flood and ebb showed weak current velocity of 0.1~0.5 m/sec or less in the northern part, 0.3~1.5 m/sec at the central, 2~3 m/sec in the southern area; and a fast velocity of maximum 3~4 m/sec around the island at the Bay entrance. The horizontal distribution of maximum flood/ebb at neap tide is shown in Fig. 7c and 7d. The current speed showed approximately 3~8 times lower than the spring tide as 0.1 m/sec or less in the north part of the Bay, below 0.3 m/sec in the middle, and under 1.5 m/sec at the entrance of the Bay.

Fig. 7.

Horizontal distribution of the tidal current velocity field at spring ebb (a), spring flood phase (b) and neap ebb (c), neap flood phase (d).

Fig. 8 shows the horizontal distribution of the tidal ellipse for the M2 (a) and S2 (b) tidal constituents. Tidal ellipse characteristics of CSB can be summarized as follows. 1) Tidal current is recti-linear with ratio of minor/major axis very small less than 0.2. 2) The orientation of the major axis was in north-south direction along the isobath. 3) Around the Jukdo island, the current pattern shows more elliptical shape. 4) The major axis of the M2 tidal constituent was approximately 2 times bigger than that of S2. To represent the tidal ellipse parameters in the CSB, ten grid points are selected (Fig. 1) and Table 4 lists ellipse parameters in terms of semi-major, minor and inclination angle and phase lag for 4 major tidal constituents. Semi major axes for M2 tidal current range from 8.3 to 130.9 cm/sec, while semi minor axes for M2 from 0.6 to 9.6 cm/sec. Semi major axes for S2 tidal current range from 2.8 to 42.5 cm/sec, while semi minor axes for S2 from 0.5 to 9.4 cm/sec. At 10 stations, the angle of the major axis of the 4 tidal constituents shows the range of 74.4~129.4o in a counterclockwise direction around the x-axis.

Fig. 8.

Horizontal distribution of tidal ellipses for M2 (a) and S2 (b) tidal constituents for the Case CS10T.

Tidal ellipse parameters for 4 major tidal constituents at selected ten grid points

3.3 Tidal Residual Current Field

Tidal residual current appears in the process of the tidal interaction with the terrain, and is known to generate and develop due to Coriolis Effect, friction, and advection terms (Robinson, 1983). Fig. 9 shows the horizontal distribution the 17-day average of tidal residual currents including both spring tide and ebb tide, of which the distribution patterns are very complicated. The pattern of tidal residual current can be summarized as follows. 1) The velocity range of the residual current in northern CSB is 1~7 cm/sec, and a large 10~30 cm / sec in the south. 2) Paired in clockwise/counterclockwise direction, several tidal residual currents exist near Jukdo at the central CSB at a velocity of 5~15 cm/sec, and at the southern entrance to CSB at 20~30 cm/sec. Paired eddy will be discussed next in vorticity characteristics. 3) Similarly to the previous study results of So et al. (1998) and Choi (2004), northwestern flows are present along the coast in the waters near KW/BN seawall north of CSB. 4) Overall, mostly northward flows were dominant in eastern CSB, and southward flows were dominant in the west. 5) Northward residual current is distributed mostly in the middle of the Bay, and southward tidal residual current is distributed along the coast. 6) The tidal residual currents appear largely around islands. This pattern is deeply related with tidal residual currents generated by flows, past specific structures such as islands or breakwaters, studied by Imasato (1983) and Robinson (1981).

Fig. 9.

Horizontal distribution of tidal residual current field in the CSB at surface.

Table 5 is a summary of the basic statistics of 10 stations among tidal residual currents of the model run CS10T. In average residual currents, u-component range showed −6.0 ~12.6 cm/sec, and v-component range, 0.7~18.3 cm/sec respectively. Except for the residual currents with smaller size, the standard deviation of the u- and v-components were noticeable 1.1~3.8 and 1.8~13.3 respectively, which indicate that the tidal residual currents have greater the variations. Especially in stations 1 and 2 at CSB south, the v-comp has reached the maximum of 51.1 and 47.7 cm/sec respectively. Through the results of the average tidal residual current, the flood/ebb imbalance can be assumed to be dominated by geographic factors.

Basic statistics of tidal residual current of model run CS10T at the selected points

Vertical component of relative vorticity (hereafter, RV) is expressed in Eqn. (2)

In computation of RV, we estimated it by using center-finite difference scheme.

The horizontal distribution of the tidal residual current shown in Fig. 9 and the RV of Fig. 10 are highly relevant. In Fig. 10a, the red means cyclonic and the blue means anti-cyclonic eddy. RV range is −40×10−5~40×10−5 sec−1; (-) indicates anti-cyclonic and (+), cyclonic direction. Fig. 10b shows some patterns of eddy paired cyclonic and anti-cyclonic. Around Jukdo at central CSB, a total of 4 major eddies, 2 cyclonic and 2 anti-cyclonic, are present; and 2 small eddies are surrounding them. 3~4 pairs of main eddy exist as very strong flows in CSB south. From the flow direction of the tidal current in the CSB, cyclonic (+) eddy was formed when bathymetry deepens, anti-cyclonic (-) eddy was formed when shallow; anti-cyclonic eddy was formed in the southeast and northwest of the island, and cyclonic eddy was formed in southwest and northeast. This eddy formation in the CSB is consistent with the eddy formation from bathymetry changes and seabed friction proposed by Robinson (1981) and Maze et al. (1998), and with the eddy formation around the headland suggested in the study of Zimmerman (1979, 1981), Signell and Harris (2000).

Fig. 10.

Horizontal distribution of relative vorticity from residual current (a) and schematic patterns of tidal residual eddy (b) in the CB at surface.

4. Discussion and Summary

In this study, numerical modeling experiments were carried out in order to understand the characteristics of tide and tidal current, and its residual current. The model results were validated with the measured data for the validation of modeling experimental results, which were verified through skills analysis using the amplitude and phase of 4 major tidal constituent (M2, S2, K1 and O1) of current velocity and sealevel. The results of skill analyses for the 4 major tidal constituents of sea level showed that mostly 90%, and the phase, 95%. The results of skill analyses for tidal current for M2, S2 tidal constituents shows higher than 90%, while those for K1 and O1 constituents were around 70~78% and 60~90%.

The characteristics of the tidal height of CSB were able to be understood through tidal height time-series of Fig. 4, and co-tidal chart in Fig. 6. Tidal height fluctuations of CSB showed the characteristics of spring-neap with two week-periods, and characteristics of semi diurnal tide. The tidal wave propagates from the entrance of the Bay towards the inside, and rotates counter-clockwise in the southern part of CSB. Tidal range was gradually increased towards KW/BN. For each tidal constituent, time lag of tidal wave showed between 10 to 30 minutes from the Bay entrance to the northern end.

As for CSB tidal current, south-north direction of the flow was dominant along the seafloor topography, with fast flow of 2~3 m/sec at the narrow waterways of the CSB entrance. The farther north of CSB, the flow rate significantly reduced compared to southern part. According to the analysis result of the tidal ellipse of the 4 major tidal constituents, the minor axis was under the size of 10% of the major axis, showing the characteristics of reversing tide, with the angle of the major axis side-by-side with the isobath. In CSB, the size of the M2 tidal constituent was about twice larger than S2 tidal constituent.

In CSB, the size of the tidal residual current was 1~7 cm/sec in the northern part and 10~30 cm/sec in the South, with the distribution very complex, and largely dependent on the geographic factors including bathymetry and coastline. The tidal residual current tends to flow towards the northwest along the KW/BN tidal embankment in northern part, but much smaller than at the entrance. Several tidal residual currents were paired cyclonic/anti-cyclonic. Characteristics of the eddy rotation were examined through relative vorticity derived from the tidal residual current. Among several eddies, with the size of hundreds of m to a few km, that were formed by bathymetry changes, seabed friction, and certain structures, 2 pairs of cyclonic/anti-cyclonic major eddies were shown nearby Jukdo in the CSB center, and 3-4 pairs of strong major eddies were formed in the south. Small eddy that occurs locally has a very high correlation with the formation of tidal residual current.

CSB is the Bay with a variety of geographic variation for over 30 years such as shoreline changes, bathymetry and volume changes due to reclamation projects. These problems are closely related to the changes in tidal phenomena due to KW/BN tidal embankment construction that So et al. (1998) studied, and also have a close relationship with changes and reduction of biota due to eutrophic inflow of fresh water that Lee (1996) and Lee et al. (1997) suggested, and with the generation of bottom layer hypoxia stratification in summer that Lee et al. (2012) reported. Based on the numerical modeling experiments, we understand tide, tidal current and tidal residual current in the CSB; and is such as density-driven current formation from the influx of freshwater from KW/BN, salinity field changes and stratification, and wind-driven current due to the influence of local winds. Understanding these physical processes will lead to help solve more complicated issues such as changes in water quality, and formation of hypoxia, degradation of the ecosystems, etc. that the CSB currently suffer from.

Acknowledgements

This research was a part of the project titled ‘2011 ChungCheong Sea Grand Program’ funded by the Ministry of Oceans and Fisheries, Korea. The authors thank anonymous reviewers for their helpful comments.

References

1. Blumberg, A.F., Mellor, G.L. A Descriptive of a Three Dimensional Coastal Ocean Circulation Model. p1-16. In: Three-dimensional Coastal Ocean Models. Coastal Estuarine Sci 1987;4ed by N.S. Heaps, AGU; Washington, D.C:
2. Cai, S., Huang, Q., Long, X. Three-dimensional Numerical Model Study of the Residual Current in the South China Sea. Oceanol Acta 2003;26(5):597–607. 10.1016/S0399-1784(03)00053-7.
3. Choi, Y.H. Development of Water Quality Prediction Model in Chunsu Bay. PhD thesis. Chungnam Natl Univ 2004;:132.
4. Dube, S.K., Rao, A.D., Shinha P.C., Jain I. Implications of Climatic Variations in the Fresh Water Outflow in the Wind-induced Circulation of the Bay of Bengal. Atmospheric Env 1995;29(16):2133–2138. 10.1016/1352-2310(94)00238-G.
5. Foreman, M.G.G., Stucchi, D.J., Zhang, Y., Baptista, A.M. Estuarine and Tidal Currents in the Broughton Archipelago. Atmos Ocean 2006;44(1):47–63. 10.3137/ao.440104.
6. Goodrich, D., Boicourt, W., Hamilton, P., Pritchard, D. Wind-induced Destratification in Chesapeake Bay. J Phys Oceanogr 1987;17:2232–2240. 10.1175/1520-0485(1987)017<2232:WIDICB>2.0.CO;2.
7. Guo, X., Yanagi T. Seasonal Variation of Residual Current in Tokyo Bay, Japan- diagnostic Numerical Experiments. J Oceanogr 1996;52:597–616. 10.1007/BF02238323.
8. Guo, X., Valle-Levinson A. Wind Effects on the Lateral Structure of Density-driven Circulation in Chesapeake Bay. Cont Shelf Res 2008;28:2450–2471. 10.1016/j.csr.2008.06.008.
9. Imasato, N. What is Tide-induced Residual current? J Phy Oceanogr 1983;13:1307–1317. 10.1175/1520-0485(1983)013<1307:WITIRC>2.0.CO;2.
10. Jung, K.Y., Ro, Y.J., Kim, B.J. Numerical Modeling Experiments of Current Circulation in the Chunsu Bay, Yellow Sea, Korea during Summer Season. Proc of Spring Meeting, 2011 of the Korean Association of Ocean Sci and Tech Soc, 146 2001a.;
11. Jung, K.Y., Ro, Y.J., Kim, B.J. Impact of the Freshwater Release on the Tidal Circulation in the Chunsu Bay, Yellow Sea, Korea based on Numerical Model. Proc of PICES 2011 Annual Meeting Prog, 192 2011b.;
12. Jung, K.Y., Ro, Y.J., Kim, B.J. Influence of Freshwater Release on the Current System in the Chunsu Bay, Yellow Sea, Korea during summer season. Proc of Autumn Meeting, 2011 of the Korean Soc of Oceanogr 2011c;:105–106.
13. Jung, K.Y., Ro, Y.J., Kim, B.J. Salinity Variation and Stratification caused by Freshwater Input in the Chunsu Bay, Yellow Sea, Korea during Summer Season. Proc of Spring Meeting, 2012 of the Korean Association of Ocean Sci and Tech Soc, 188 2012a.;
14. Jung, K.Y., Ro, Y.J., Kim, B.J. Observation and Analysis of Hydrodynamic and Hydrography in the Chunsu Bay, Yellow Sea, Korea, 2010-2011. Proc of Spring Meeting, 2012 of the Korean Association of Ocean Sci and Tech Soc, 198 2012b.;
15. Jung, K.Y., Ro, Y.J., Kim, B.J. Tracking Patterns of Freshwater from Kanwol/Bunamho based on Particle Trajectory Modeling Experiments in the Chunsu Bay, Yellow Sea, Korea. Proc of Autumn Meeting, 2011 of the Korean Soc of Oceanogr 2012c;:75–76.
16. Jung, K.Y., Ro, Y.J., Kim, B.J. Tidal and Sub-tidal Current Characteristics in the Central Part of Chunsu Bay, Yellow Sea, Korea during the Summer Season. ‘The Sea’ J Korean Soc of Oceanogr 2013a;18(2):53–64. 10.7850/jkso.2013.18.2.53.
17. Kashiwai, M. Tidal Residual Circulation Produced by a Tidal Vortex. Part1. Life-history of a Tidal Vortex. J Oceanogr Soc Japan 1984;40(6):279–294. 10.1007/BF02302521.
18. Lee, J.S., Kim, K.H., Sim, J.H., Han, J.H., Choi, Y.H., Khang B.J. Massive Sedimentation of Fine Sediment with Organic Matter and Enhanced Benthic-pelagic Coupling by an Artificial Dyke in Semi-enclosed Chonsu Bay, Korea. Mar Pollut Bull 2012;64:153–163. 10.1016/j.marpolbul.2011.09.033.
19. Lee, T.W., Choi, M.S., Yang, S.Y., Ma, C.W., Ro, Y.J., Park, S.C. A Study on the Environment Investigation and Fishery Utilization in the Chunsu Bay. Final Report, The Province of Chungcheongnam-do 2011;:534.
20. Lee, T.W. Change in Species Composition of Fish in Chonsu Bay 1.Demersal Fish. Korean J. of Fisheries and Aquatic Sci 1996;29(1):71–83.
21. Lee, T.W., Moon, H.T., Choi, S.S. Change in Species Composition of Fish in Chonsu Bay 2. Surf Zone Fish. J Ichthyological Soc of Korea 1997;9(1):79–90.
22. Martin, J., McCutcheon, S.C. Hydrodynamics and Transport for Water Quality Modeling. Lewis Publishers 1999;:794.
23. Mastumoto, K., Takanezawa, T., Ooe, M. Ocean Tide Models Developed by Assimilating TOPEX/POSEIDON Altimeter Data into Hydrodynamical Model: A Global Model and Regional Model Around Japan. J Oceanogr 2000;56:567–581. 10.1023/A:1011157212596.
24. Maze, R., Langlois, G., Grosjean, F. Tidal Eulerian Residual Currents over a Slope, Analytical and Numerical Frictionless Models. J Phys Oceanogr 1998;28:1321–1332. 10.1175/1520-0485(1998)028<1321:TERCOA>2.0.CO;2.
25. Orlanski, I. A Simple boundary Condition for Unbounded Hyperbolic Flows. J Comut Phys 1976;21:251–269. 10.1016/0021-9991(76)90023-1.
26. Park K., Oh J.H. Calibration and Verification of a Hydrodynamic Model in Chunsu Bay and Adjacent Coastal Water. J of Korean Soc of Coastal and Ocean Eng 1998;10(3):109–119.
27. Pawlowicz, R., Beardsley, B., Lentz, S. Classical Tidal Harmonic Analysis including Error Estimates in MATLAB using T_TIDE. Comput. Geosci 2002;28:929–937. 10.1016/S0098-3004(02)00013-4.
28. Robinson, IS. Tidal Vorticity and Residual Circulation. Deep Sea Res 1981;28A(3):195–212. 10.1016/0198-0149(81)90062-5.
29. Robinson, IS. Tidally Induced residual Flows, In: Physical Oceanography of Coastal and Shelf Seas edited by B, Johns. 1983;:321–356. Elsevier; New York:
30. Signell, RP., Harris CK. Modeling Sand Bank Formation around Tidal Headlands. In: 6th International Conference of ASCE, New Orleans, LA, 3-5 Nov 1999 2000;:209–222.
31. So, J.K., Jung, K.T., Jang, W.C. Numerical Modeling of Tides and Tidal Currents Cuased by Embankment at Chunsu Bay. J of Korean Soc of Coastal and Ocean Eng 1998;10(4):151–164.
32. Yanagi, T. General Mechanism of the Tidal Residual Circulation. J Oceanogr Soc Japan 1983;35(6):241–252. 10.1007/BF02070259.
33. Yoo, I.H. Numerical Modeling of current and diffusion in Chunsu Bay. MS thesis, Chungnam Natl Univ 1992;:64.
34. Zhai, L., Sheng, J., Greatbatch, R.J. Baroclinic Dynamics of Wind-driven Circulation in a Stratified Bay: A Numerical Study using Models of Varying Complexity. Cont Shelf Res 2008;28:2357–2370. 10.1016/j.csr.2008.05.005.
35. Zimmerman, J. On the Euler-Lagrange Transformation and the Stoke’s Drift in the Presence of Oscillatory and Residual Currents. Deep-Sea Res 1979;26A:505–520. 10.1016/0198-0149(79)90093-1.
36. Zimmerman, J. Dynamics, Diffusion and Geomorphological Significance of Tidal Residual Eddies. Nature 1981;290:549–555. 10.1038/290549a0.

Article information Continued

Fig. 1.

Map for the study area with locations for the various stations in the Chunsu Bay (CSB), Yellow Sea of Korea.

Fig. 2.

Grid system for the numerical model for the CSB, Yellow Sea of Korea. And location of datasets for model calibration and validation.

Table 1.

Amplitude (m) and phases (deg.) of tidal constituents used for model forcings at open boundary points

Boundary Index M2 S2 K1 O1
Amp. Phase Amp. Phase Amp. Phase Amp. Phase
South Start 2 1.98 89.8 0.73 143.0 0.34 270.9 0.22 256.4
End 78 2.16 90.1 0.81 142.7 0.34 268.8 0.14 263.5
West Start 2 1.98 89.8 0.73 143.0 0.34 270.9 0.22 256.4
End 135 2.14 101.8 0.82 154.6 0.36 276.5 0.26 257.6

Fig. 3.

Time series of raw current meter records (a) and the residual current (c) at surface during Jul. 28~Aug. 30, 2010 at Jukdo Station.

Fig. 4.

Time series of sea level elevation at Boryung Station in the CB. Elevation compared to model result (CS10T).

Table 2.

Comparison between model output and observations of sea level at three points in terms of 4 tidal constituents (unit: amplitude (m), phase (deg.) and skill score (%))

Location M2 S2 K1 O1
Amp. Phase Amp. Phase Amp. Phase Amp. Phase
Boryung
(KHOA)
Obs. 2.10 94.4 1.07 145.1 0.35 281.0 0.26 250.5
Model 2.08 99.1 1.03 152.9 0.32 261.7 0.21 258.5
Skill Score 99.0 98.7 96.3 97.8 91.4 94.6 80.8 97.8
Obs. 2.37 88.6 0.83 140.8 0.38 271.4 0.27 246.0
Kojeong
(Choi,2004)
Model
(present)
2.10 99.7 0.94 149.3 0.33 261.9 0.22 258.6
Skill Score
vs Obs.
88.6 96.9 86.7 97.6 86.8 97.4 81.5 96.5
Kanwol
(Choi,2004)
Obs. 2.56 107.1 1.01 170.0 0.48 287.7 0.28 254.2
Model 2.23 101.6 1.02 156.7 0.33 263.5 0.22 258.8
Skill Score
vs Obs.
87.1 98.5 99.0 96.3 68.8 93.3 78.6 98.7
Average Skill Score 91.6 98.0 94.0 97.3 82.3 95.1 80.3 97.7

Table 3.

Comparison between model output and ADCP record for tidal current at Jukdo Station in terms of 4 tidal constituents (unit: amplitude (cm/s), phase (deg.) and skill score (%))

Location Component M2 S2 K1 O1
Amp. Phase Amp. Phase Amp. Phase Amp. Phase
Jukdo
(Jung,2012)
U-comp. Obs. 4.87 213.2 3.89 302.1 0.42 328.8 0.24 227.2
Model 5.34 206.8 3.51 280.8 0.53 49.2 0.31 11.3
Skill Score 90.3 98.2 90.2 94.1 73.8 77.7 70.8 60.0
V-comp. Obs. 13.63 26.7 6.23 103.1 1.93 304.9 0.49 95.3
Model 14.33 21.4 5.91 89.0 1.39 215.4 0.61 134.8
Skill Score 94.9 98.5 94.9 96.1 72.0 75.1 75.5 89.0

Fig. 5.

Time series of observed (ADCP) and model predict at Jukdo Station in the CB (a: U-component; b: V-component).

Fig. 6.

Co-tidal chart for M2 (a) and S2 (b) constituents in the CSB.

Fig. 7.

Horizontal distribution of the tidal current velocity field at spring ebb (a), spring flood phase (b) and neap ebb (c), neap flood phase (d).

Fig. 8.

Horizontal distribution of tidal ellipses for M2 (a) and S2 (b) tidal constituents for the Case CS10T.

Table 4.

Tidal ellipse parameters for 4 major tidal constituents at selected ten grid points

Station
(#)
M2 S2
Sema
(cm/s)
Semi
(cm/s)
inc.
(deg.)
pha
(deg.)
Sema
(cm/s)
Semi
(cm/s)
inc.
(deg.)
pha
(deg.)
1 116.3 8.0 105.9 344.2 37.9 7.7 109.9 91.0
2 130.9 9.6 77.8 339.7 42.5 9.4 81.5 88.2
3 51.9 3.4 95.9 344.1 17.3 3.2 97.8 98.8
4 47.4 3.1 98.1 350.8 15.7 3.6 102.0 99.7
5 52.8 4.3 88.3 356.0 18.4 4.0 93.4 107.8
6 45.4 3.5 96.0 352.5 15.6 3.3 100.0 100.5
7 32.5 2.0 115.6 3.0 11.3 2.1 118.9 113.3
8 19.6 1.5 99.1 359.7 6.7 1.5 99.6 108.0
9 8.3 0.6 110.7 357.6 2.8 0.5 112.3 101.0
10 18.0 1.1 129.4 352.1 6.3 1.1 133.2 92.4
K1 O1
1 11.5 0.9 101.7 238.1 5.1 0.9 93.6 114.4
2 12.2 0.9 74.0 236.8 6.0 0.9 74.4 111.1
3 5.6 0.5 90.5 254.7 2.4 0.4 90.8 133.0
4 4.8 0.4 96.4 248.5 1.6 0.4 97.5 130.8
5 6.2 0.6 90.7 251.9 3.2 0.6 94.9 130.9
6 4.9 0.4 90.3 248.0 2.4 0.4 86.9 125.0
7 4.2 0.3 115.6 258.6 2.3 0.3 118.9 140.3
8 2.3 0.2 96.1 252.8 1.3 0.2 99.7 131.6
9 0.9 0.1 99.3 241.6 0.5 0.1 121.1 124.5
10 11.5 0.9 101.7 238.1 5.1 0.9 93.6 114.4

Fig. 9.

Horizontal distribution of tidal residual current field in the CSB at surface.

Table 5.

Basic statistics of tidal residual current of model run CS10T at the selected points

Station
(#)
U-comp. (cm/s) V-comp. (cm/s)
Mean St. dev Mini Maxi Mean St. dev Mini Maxi
1 −4.6 3.0 −15.7 3.6 18.3 11.4 −3.4 51.1
2 12.6 3.5 4.8 28.3 6.9 13.3 −29.0 47.7
3 1.0 3.8 −21.7 18.4 0.7 6.6 −11.1 16.9
4 −6.0 2.7 −21.0 −0.9 7.1 6.5 −4.9 23.6
5 −1.2 1.1 −6.4 0.8 2.3 6.1 −9.6 17.6
6 3.0 1.5 −2.7 6.8 1.6 4.9 −7.3 15.7
7 −4.0 1.7 −8.4 0.5 4.6 3.9 −3.2 15.6
8 0.2 0.4 −1.0 1.1 2.4 2.3 −2.1 8.7
9 0.3 0.5 −0.9 1.5 1.3 0.8 −1.0 3.3
10 0.3 1.7 −4.8 3.9 2.8 1.8 −1.4 8.6

Fig. 10.

Horizontal distribution of relative vorticity from residual current (a) and schematic patterns of tidal residual eddy (b) in the CB at surface.