| 摘要: | 海岸帶作為陸地、海洋與大氣的交界區域，是陸源物質進入海洋的主要通道。進入近岸海域的陸源物質，一旦其濃度超過了海域涵容能力，將造成物質循環的平衡狀態改變，進而衝擊海域的生態環境。影響海域涵容能力的因素中，最重要的是水體混合與擴散效率。受近岸水動力影響，水體混合與擴散特性在時間上非定常、空間上非均勻，從理論、數值模式計算以了解其細部特徵都是非常困難的。若能藉由有效的觀測增進對於陸源物質在近岸海域的輸送與移運過程之了解，對陸源物質的擴散效率及其與水動力的相互作用特性的掌握將會更準確。
依據地貌分類，海岸帶包括海灘、濕地、河口與瀉湖等區域，從水動力特性角度分類，海岸帶可劃分為碎波帶、內陸棚、外陸棚區域。對於海氣通量研究，濕地及瀉湖扮演重要角色，對於海陸通量研究，陸源物質在河口及內陸棚的混合與擴散則扮演重要角色。水體的混合與擴散方面，傳統水體交換率、滯留時間等手段無法呈現河口混合行為的細部結構，也無法反映層化效應的影響。本研究研製漂流浮球陣列作為主要觀測工具，分別就內陸棚海域的混合與擴散特性、河口區域層化效應對混合的影響進行研究。研究成果分述如下：
（1）本研究採商用GPS定位晶片研製低成本且適用於近岸海域觀測的漂流浮球陣列。本研究改進定位精度、提高採樣頻率、減低製作成本、減少體積與重量、改進資料回傳及交互操作方式。針對低成本GPS晶片在低速下會發生定位跳躍的問題，本文基於ARIMA模式建立資料處理與品管方法，消除了該問題的影響，提高數據品質。透過在大型斷面水槽的漂流浮球觀測測試與校驗實驗，與ADV流速計之觀測結果比對之RMSE約為0.02 m/s。
（2）基於Stokes Boundary Layer理論，提出了一個新的內陸棚海域混合與擴散研究模型，作為理論框架。該理論中描述邊界層內速度剖面變化的參數κ1與描述流速反轉時間提前的參數κ2是兩個可用於分析海域混合效率的中間參數。進而本研究設計現場實驗檢驗並建立了一個更具可行性的內陸棚海域混合與擴散行為觀測方法。在桃園觀音近岸海域進行的不同海況下的漂流浮球陣列觀測實驗，證實理論所預測的潮流轉向時間早於水位轉向時間，隨著離岸距離越近時間提前量會增加；沿岸方向之流速會隨著離岸越近而快速下降。以漂流浮球陣列實測資料代入理論擬合得到不同時空下的κ1與κ2組合。這兩個參數顯著相關，κ2=6.02*κ1-0.0008（r = 0.97, p < 0.05）。進一步直接計算得實驗海域的υ，結果發現速度衰減項引致的υ約為相位差異項引致的3至12倍，二者均隨量距離的增加成指數形式降低。透過對比不同海況下的擬合結果，本文發現當波浪較小時，Stokes Boundary Layer理論的解釋度佳，說明此時潮流在近岸混合中扮演主要的角色；當波浪較大時，Stokes Boundary Layer理論的解釋度降低，是由於波浪引致的底邊界層增加了底床粗糙度，額外增強了混合強度所致。
（3）本文以淡水河口及核三廠溫排水口作為研究區域，探究層化效應在不同時空尺度上對河口混合的影響。淡水河口區域以海研二號研究船進行了多個航次的調查，基於三維溫鹽觀測資料，以其中2016年10月2202航次資料作為高層化效應之代表例，並以2017年3月2225航次資料作為垂向均勻混合之案例，比對二者差異。基於營養鹽濃度資料分析發現淡水河口在垂向均勻混合時的營養鹽濃度衰減速率是層化效應強時的約2倍。淡水河口DIN濃度隨鹽度衰減的速率比珠江河口快6%至46%，漂流浮球陣列觀測之淡水河口延散係數比同等規模河口的大1至2個量級。淡水河口的層化效應會阻礙營養鹽向上傳遞及溶解氧向下傳遞，改變河口氮循環進而可引發河口初級生產力下降等生態問題，本研究提供了觀測證據。
基於溫排水口航次調查之三維溫鹽資料分析顯示溫排水口漲潮期間層化效應強，退潮期間垂向均勻混合。基於溫度漂流浮球量測之水溫資料，在垂向均勻混合時溫排水口的溫度衰減速率e-folding scale是層化效應強時的約2倍。基於漂流浮球陣列觀測資料計算的溫排水口剪力流延散係數，發現漲、退潮期間的剪力流延散係數差異可達26.0%。本文透過溫排水口漂流浮球陣列觀測實驗量化了剪力流延散係數隨著Richardson數變化的趨勢。;As the boundary area between land, sea and atmosphere, the coastal zone is the main channel for terrestrial materials transferring into the ocean. The balance of the material cycle will be changed once the terrestrial materials exceed the capacity of the ocean self-purification, which will cause ecological pollution problems. Among the factors affecting the capacity of self- purification, the most important factor is the mixing and diffusion efficiency of water. The near-shore hydrodynamics as well as the mixing and diffusion characteristics are unsteady in time and spatially non-uniform. As a result, it is very difficult to grasp the detail information of mixing and diffusion behavior by theory analysis or by numerical model. If effective observations can be applied to improve the understanding of the mixing and transport processes of terrestrial materials in coastal waters, the evaluation of diffusion efficiency of terrestrial materials and its interaction with hydrodynamics will be more accurate.
According to the classification of landforms, the coastal zone includes beaches, wetlands, estuaries and lagoons. From the perspective of hydrodynamic characteristics, the coastal zones can be divided into surf-zone, inner shelf, middle shelf and outer shelf. For air-sea flux research, wetlands and lagoons have important impacts; for the study of the mixing and diffusion of terrestrial materials in the sea, estuaries and inner shelf play important roles. The study of mixing and diffusion of water has traditionally started with the exchange rate and residence time of water. However, these coefficients can neither represent the details of the estuary mixing behavior, nor represent the influence of the stratification effect. This paper develops the sea surface drifter array as the main observation tool to carry out researches on a) the mixing and diffusion characteristics of the inner shelf and, b) the influence of the stratification effect on the estuary mixing, specific as follows.
(1) This study developed the miniature sea surface drifters with the commercial low-cost GNSS positioning chip. On account of shallow water depth and small hydrodynamic temporal-spatial scale of inner shelf or estuaries, this study improves the positioning accuracy, increases the sampling frequency, lower the production cost, reduces the volume and weight, and improves the data post back method and interaction operation system of sea surface drifter. Aiming at the problem of positioning jumping under low speed, this study established a data processing and quality control method based on Auto Regressive Integrated Moving Average model. The influence of this problem is eliminated and the data quality of sea surface drifter is improved after applying this QC method. All the drifters are calibrated in the 200-m length water channel and the RMSE is ~0.02 m/s compared with the Acoustic Doppler Velocimeter which is installed in the middle of the water channel.
(2) Based on Stokes Boundary Layer theory, a new research framework of inner shelf mixing and diffusion is proposed. The parameter κ1 describing the change of velocity profile in the boundary layer and the parameter κ2 describing the current reversal time are two intermediate parameters that can be used to analyze the mixing efficiency of the sea area. And then this study designed the field experiments and establisded a more feasible observation method of mixing and diffusion behavior in inner shelf. This paper conducted several drifter observation experiments under different sea conditions in the coastal waters of Guanyin, Taoyuan City. It was found that the time of tide reversal is earlier than that of water level reversal, and the time difference increase as the distance to the coast decreases. The alongshore current velocity decreases rapidly as it gets closer to the coast line. Above-mentioned two parameters, κ1 and κ2, were fitted by the field drifter data. These two parameters were significantly correlated, κ2=6.02*κ1-0.0008 (r = 0.97, p < 0.05). Further, the eddy viscosity coefficient of the study area was calculated based on these two parameters. It is found that the eddy viscosity coefficient caused by the velocity attenuation term is about 3 to 12 times that of the phase difference term. Both of them decrease exponentially with increasing distance to the coast line. By comparing the fitting results under different sea conditions, this study found that the Stokes Boundary Layer theory has a good interpretation when the waves are small (Hs < 0.3 m), indicating that the tide also plays a role in coastal mixing. When the waves are large, Stokes Boundary Layer stays at the lower degree of interpretation.
(3) This study selects the Danshui estuary and thermal discharge outlet of No.3 nuclear power plant as the research areas to explore the influence of stratification effect on mixing at different temperial and spatial scales. A number of voyage surveys were conducted with the OR-2 research vessel around Danshui estuary. The 2202 voyage conducted in October 2016 was selected as a representative example of strong stratification condition and the 2225 vyage conducted in March 2017 was selected as an example of verticle well mixing condition based on the 3-D water temperature and salinity observation data. It is found that the decay rate of nutrient concentration under weak stratification condition is about 2 times of that under strong stratification condition in the Danshui estuary. The decay rate of DIN concentration of Danshui estuary is 6% to 46% larger than that of the Pearl River estuary, and the dispersion coefficient of Danshui estuary is one to two orders of magnitude larger than that of an equivalent estuary. The stratification effect of the Danshui estuary will hinder the upward transfer of nutrients and the downward transfer of dissolved oxygen, and then changes the nitrogen cycle of estuary, which will lead to ecological problems such as the decline of primary productivity in the estuary. This study provides observational evidence.
The analysis of the 3-D water temperature and salinity data shows that the stratification effect of the water around thermal discharge outlet is strong during flood tide and weak during ebb tide. The temperature decay rate under weak stratification condition is 2 times faster than that under strong stratification condition based on the e-folding scale. The shear dispersion coefficient is calculated based on the drifter data. And it is found that the difference of dispersion coefficients under different stratification strength can reach 26.0%. This study quantify the variation trend of dispersion coefficient with the increasing Richardson number. |
