閩江河口潮灘季節性沖淤變化格局及其控制機制

王愛軍，李海琪，葉翔，梁灝深，張望澤,5，吳水蘭,6，冉暢，陶舒琴，劉子同,7，于謙

1.自然資源部第三海洋研究所海洋與海岸地質研究室，廈門 361005

2.福建省海洋物理與地質過程重點實驗室，廈門 361005

3.自然資源部海峽西岸海島海岸帶生態系統野外科學觀測研究站，廈門 361005

4.上海海洋大學海洋科學學院，上海 200090

5.福州大學先進制造學院，泉州 362251

6.南京信息工程大學海洋科學學院，南京 210044

7.汕頭大學海洋科學研究院，汕頭 515063

8.南京大學海岸與海島開發教育部重點實驗室，南京 210023

潮灘是在潮汐作用顯著、細顆粒沉積物供應豐富條件下形成的一種地貌類型，位于陸地、海洋、大氣及人類活動相互作用的敏感地帶。潮灘廣泛分布在除基巖海岸和高緯度（70°～73°以上）海岸外的沿海地區，并且隨著潮差、波浪、物質供應、植被類型等不同而表現出顯著的空間差異[1-4]。目前全球 潮 灘 面 積 約 為127 921 km2，主 要 分 布 在 亞 太 地區，如印尼、中國、印度等國家[5]。我國的潮灘規模大且分布廣泛，其中長江口-杭州灣及以北沿海地區主要為平原型潮灘，而長江口-杭州灣以南的浙、閩、粵、桂等沿海地區主要為港灣型潮灘[6]。從20 世紀初開始，各國學者對潮灘沉積和地貌特征、潮灘沉積物輸運和堆積的過程、鹽沼植被對潮灘沉積過程的影響、極端事件對潮灘地貌演化影響等方面進行了系統研究和總結[7-16]。進入21 世紀，隨著現場觀測技術和數值模擬的不斷改進，我國科研人員開始關注潮灘沉積過程及其控制機制、植被-泥沙-地貌耦合的生物地貌過程、極端事件影響下的潮灘沉積物輸運與地貌演化過程、潮灘沉積層序成因機制及潮灘演化趨勢模擬等方面，取得了豐碩的成果[17-33]。

潮灘作為重要的海岸帶生態系統，不僅在藍色碳匯方面發揮著重要作用[34-35]，而且在海岸防護方面起到至關重要的作用[36]。然而，隨著海平面不斷上升、地面沉降、流域及海岸帶人類活動的不斷加劇，潮灘面積逐漸減少[5]，其中人類活動是主要驅動因素[37]。潮灘面積的減少不僅導致藍色碳匯能力減弱[5,35]，而且引起諸多地區海岸遭受不同程度侵蝕[38]。高強度的人類活動顯著影響了潮灘局部沉積物輸運過程及地貌演化過程[20,24]，而流域輸沙量的銳減則引起了潮灘沖淤格局的改變[39]。隨著流域及河口地區人類活動強度的不斷增大，河口系統狀態正在發生轉換[40]，直接影響著河口潮灘未來演化趨勢及其功能發揮。

本文選擇我國東南沿海典型的山溪性中小型河流河口潮灘作為研究區域，通過沉積物采樣、定點沖淤觀測和沉積動力學觀測等手段，初步分析了河口潮灘的季節性沖淤變化過程及其控制機制，以期為科學評估潮灘的海岸防護和潮灘生態系統保護與修復能力提供科學依據。

1 研究區概況

閩江發源于福建與江西省交界的武夷山東麓及仙霞嶺南麓，總體呈NW-SE 向橫貫福建中北部，經福州匯入東海南部，全長2 872 km，其中干流全長近577 km，流域面積為60 992 km2，多年平均徑流量為605.5×108m3，多年平均輸沙量為750×104t，并且表現出明顯的洪、枯季變化，即洪季（3—8 月）徑流量約占全年的76%，輸沙量約占全年的92%[41]。近年來隨著流域水庫建設的不斷增加，平均輸沙量已經下降為原來的三分之一[42]。

受地轉偏向力及臺灣海峽地形效應作用，閩江河口區潮汐作用強，實測最大潮差達到7.04 m，平均潮差4.46 m（梅花站），屬正規半日潮，河口內潮流基本呈往復流，其流向與河槽線一致，河口外潮流具有一定的旋轉流特征[41]。閩江口外海面開闊，夏季以西南風為主，冬季以東北風為主，波浪作用強，波浪的主要形式為風浪及涌浪同時存在的混合浪，多年平均波高為1.1 m，歷史最大波高為6.5 m[41]。此外，閩江口海域是我國東南沿海地區遭受臺風影響最為嚴重的區域之一，平均每年影響的臺風次數為5.7 個，最大臺風增水值可達2.52 m[43]。

閩江河口地區的潮灘主要分布在川石島南部、瑯岐島東部以及梅花水道南側海岸（圖1a），其中瑯岐島潮灘寬度為600～1 100 m（圖1b），灘面平緩?，F場調查結果顯示（圖1c），瑯岐島潮灘由岸向海主要分布有互花米草、海三棱藨草、光灘等，其中互花米草鹽沼灘面平均坡度約為0.36°，海三棱藨草鹽沼灘面平均坡度約為0.14°，潮間帶中下部的光灘平均坡度約為0.08°，河道邊緣地區的灘面坡度約為1.84°。

圖1 研究區位置（a）、沉積物采樣站位（b）、剖面形態及觀測站位圖（c）Fig.1 Sketch map of study area (a), sediment sampling sites (b), morphological character of profile and observation sites of tidal flat of Langqi Island (c)

2 研究方法

2.1 野外采樣與觀測

在瑯岐島東部潮灘選擇一條剖面利用RTK 對剖面進行了灘面高程測量，并在剖面上由岸向海設置了10 個采樣站位（圖1b），分別于2022 年8 月和2023 年2 月采集了表層1 cm 厚度的沉積物。為了掌握潮間帶不同空間位置的沖淤變化特征，在互花米草鹽沼邊緣（L03 站）、海三棱藨草鹽沼中部（L05 站）及光灘外緣（L09 站）設置了3 個站位，安裝了沖淤觀測架，在觀測架上設置12 個點，定期測量每個點距離灘面的高度，然后對各站位每次的測量數據進行平均，獲得各站位的沖淤變化信息（圖1c）。

為了解研究區水動力狀況，在海三棱藨草內部靠近互花米草邊緣設置一個沉積動力學觀測站位（圖1c，L04 站），利用Nortek 公司生產的聲學多普勒海流計（ADV，VECTOR 6MHz）分別于2022 年1 月25 日至2 月1 日、7 月21—31 日開展近底部潮流觀測，儀器探頭距離灘面30 cm，觀測點位于灘面以上15 cm 處；觀測設置為脈沖模式進行采樣，采樣間隔為10 min，其中冬季采樣頻率為4 Hz，每次采集1 024 組數據，夏季采樣頻率為16 Hz，每次采集4 096組數據；冬季利用北京海洲賽維科技有限公司生產的波潮儀（TWaves）進行淹沒水深觀測，以脈沖模式進行采樣，以10 min 為采樣間隔，采樣頻率設置為8 Hz，每次采集時長為5 min，取平均值。

2.2 樣品分析與數據處理

2.2.1 沉積物粒度分析

在實驗室內將采集回來的常溫保存樣品充分混合，取約2 g 重的沉積物加入濃度為5%的H2O2溶液以去除有機質，再加入濃度為1 mol/L 的HCl溶液以去除碳酸鹽，用蒸餾水清洗至溶液為中性，加入濃度為5‰的六偏磷酸鈉溶液（(NaPO3)6）浸泡24 h，并利用超聲波震蕩1 min，使樣品充分混和、分散后，用英國馬爾文公司生產的 Mastersizer 2000 型激光粒度儀進行粒度分析，并用儀器自帶軟件導出中值粒徑和各粒級組分數據，采用謝帕德分類體系對沉積物進行分類[44]。

2.2.2 懸沙濃度校對

利用在L04 站采集的表層沉積物在室內配置不同懸沙濃度（SSC）的水體與ADV 采集的聲學信號進行相關分析，結果發現兩者之間具有很好的相關性（圖2）。因此，根據該回歸方程可以將ADV 觀測獲得的信噪比（Signal-Noise Ratio，SNR）轉化成時間序列的SSC。

圖2 懸沙濃度與ADV 觀測獲得的信噪比的相關關系Fig.2 Relationship between SSC and SNR recorded by ADV

2.2.3 沉積動力學參數計算

利用儀器自帶軟件將觀測數據導出，然后再利用該軟件的數據轉換功能，將原始數據轉換成數據文件以便計算水平方向上的平均流速和波浪參數。

潮間帶地區作為極淺水環境，整個水層都屬于邊界層范圍內，平均水平流速在垂向上的分布遵循對數分布[45]：

其中，U(z)是平均水平流速度，u*是摩阻速度；κ是卡門常數（κ=0.4)；z是距底床高度，z0是底部粗糙長度，τb是潮流引起的底部切應力，ρw是海水的密度。一般在邊界層范圍內，可以利用對數剖面方法計算潮流引起的底部切應力。

有效波高（Hs）采用如下公式計算[46]：

式中，Sη是水位的功率譜。將ADV 測得的壓力轉換到水位時，需要考慮衰減，根據線性波理論，衰減系數可以表示為[47]：

其中，k是波數，h是平均水深，z是壓力傳感器的深度（負值）。波浪引起的底部切應力（τw）可以用如下公式計算[45]：

其中，fw為波浪摩擦因子，Uw為波浪軌跡速度，可以表達為波高（Hs）、波周期（T）和水深（h）的函數[45]：

其中k為波數，可以表達為：k= 2π/L，L為波長，可以用圖解法計算得到[48]。波浪摩擦因子可以表達為：

式中A為半軌跡距離，可以表達為：A=UwT/2π。

由流速和懸沙濃度數據計算近底層懸沙瞬時水平輸運率f(t)及潮周期內懸沙凈水平輸運通量F(t)[45]：

其中U(z,t)和C(z,t) 分別表示距離海底高度z(m)處t時刻的瞬時水平流速（ua，m/s）和懸沙濃度(SSC, mg/L)，Δt表示實測時刻代表時間長度。

研究區表層沉積物冬季以粗顆粒沉積物為主，夏季以細顆粒沉積物為主，因此，在計算沉積物臨界侵蝕切應力時分別選擇適合非粘性沉積物和粘性沉積物的計算公式。非粘性沉積物臨界侵蝕切應力時選擇如下計算公式[45]：

其 中，θ是Shields 參 數，θcr是 沉 積 物 臨 界 起 動Shields 參 數，s=ρs/ρw，ρs是 沉 積 物 密 度（取 值 為2 650 kg/m3），ρw是海水密度（取值為1 025 kg/m3），d是 沉 積 物 粒 徑， g 是 重 力 加 速 度（ 取 值 為9.81 m/s2）,D*是一個與沉積物粒徑有關的無量綱參數，ν是水體的運動黏滯系數，取值為1.36×10-6m2/s。對于黏性沉積物，其臨界侵蝕切應力計算用如下公式[49]：

其中，β是與顆粒粒徑有關的系數，根據試驗資料[49]，在中值粒徑小且表層沉積物含水量較大時，取β=0.3，W為含水量，定義為樣品中水的質量與沉積物質量之比，根據采集的現場表層底質樣品的分析獲得。

3 研究結果

3.1 沉積物粒度的季節分布特征

粒度分析結果顯示（圖3、圖4），瑯岐島潮灘表層沉積物粒度組分以粉砂和砂為主，黏土含量相對較少，呈現出顯著的時空差異。從季節方面看，冬季表層沉積物砂含量顯著增大，黏土含量顯著減小，沉積物中值粒徑Φ 值較??；夏季則相反，砂含量顯著減小，黏土含量和粉砂明顯增大，沉積物中值粒徑Φ 值較大（表1）。從空間上看，潮間帶上部的互花米草鹽沼內表層沉積物組分以粉砂為主，平均含量為60.6%，砂含量平均為21.0%，黏土含量平均為18.4%，中值粒徑平均為5.60 Φ；在潮間帶中部的海三棱藨草鹽沼，沉積物組分以砂為主，平均含量為57.1%，粉砂含量平均為32.4%，黏土含量平均為10.5%，中值粒徑平均為4.31 Φ；在潮間帶下部的光灘地區，表層沉積物組分以粉砂為主，平均含量為58.7%，砂含量平均為23.0%，黏土含量平均為18.3%，中值粒徑平均為5.56 Φ（圖3、圖4）。

表1 瑯岐島潮間帶不同季節表層沉積物粒度組分及中值粒徑Table 1 Seasonal variations of composition and median grain-size of surficial sediment

圖3 瑯岐島潮灘夏季（左）和冬季（右）表層沉積物粒度組成Fig.3 Surficial sediment composition along the profile of tidal flat in summer (left) and winter (right)

圖4 瑯岐島潮灘不同季節表層沉積物中值粒徑Fig.4 Seasonal variations of median grain-size of surficial sediment along the profile of tidal flat

根據Shepard 三角分類圖，瑯岐島潮間帶表層沉積物類型以黏土質粉砂（YT）為主，并且呈現出顯著的季節變化，其中夏季表層沉積物黏土質粉砂占80%，僅在潮間帶中部出現砂質粉砂（ST）和砂-粉砂-黏土（STY）；冬季以砂（S）和砂質粉砂為主，其次為粉砂質砂（TS）（圖5）。

圖5 瑯岐島潮灘不同季節表層沉積物類型Fig.5 Seasonal variations of surficial sediment type along the profile of tidal flat based on Shepard classification

3.2 潮灘侵蝕-淤積的時空變化特征

潮間帶不同空間位置的沖淤觀測結果顯示，瑯岐島潮間帶灘面沖淤表現出顯著的時空變化特征?；セ撞蓰}沼邊緣灘面表現出總體持續淤積的變化特征，僅在冬末至初春季節出現輕微侵蝕（0.55 cm），觀測期間年淤積厚度為1.22 cm（圖6a）。海三棱藨草中部灘面在初夏至初秋季節表現為持續淤積，淤積厚度為14.13 cm；初秋至冬末季節表現為持續侵蝕，侵蝕深度為20.13 cm；冬末至春末期間呈現出緩慢淤積，淤積厚度為2.54 cm；整個初夏至翌年春末表現為凈侵蝕，侵蝕深度為3.42 cm（圖6b）。光灘外緣灘面在初夏至秋末季節表現為持續侵蝕，侵蝕深度為5.90 cm；秋末至翌年春末季節表現為持續淤積，侵蝕深度為5.38 cm，整個初夏至翌年春末表現為凈侵蝕，侵蝕深度為0.52 cm（圖6c）。

圖6 潮間帶不同站位沖淤變化過程正值為淤積，負值為沖刷。Fig.6 Variations of erosion and accretion along the section in different stations Positive value means accretion, negative value means erosion.

3.3 潮間帶水動力特征

瑯岐島潮間帶中上部L04 站位淹沒水深觀測結果顯示（圖7a、8a），冬季和夏季觀測期間大潮最大淹沒水深分別為2.20 和1.86 m，小潮最大淹沒水深分別為1.37 和0.90 m。潮流觀測結果顯示（圖7b、8b），冬季小潮觀測期間近底部潮流流速為0.008 ～0.127 m/s，大潮期間為0.004 ～0.187 m/s；夏季小潮觀測期間近底部潮流流速為0.001 ～0.046 m/s，大潮期間為0.001 ～0.060 m/s；無論冬季還是夏季，潮周期內潮流流速最大值均出現在落潮后期。波浪計算結果顯示（圖7c、8c），冬季和夏季觀測期間有效波高最大值分別為0.63 和0.20 m，平均值分別為0.26 和0.07 m，冬季波浪明顯較夏季強。底部切應力計算結果顯示（圖7d、8d），冬季和夏季觀測期間潮流引起的底部切應力最大值分別為0.551 和0.057 N/m2，平均值分別為0.067 和0.003 N/m2；波浪引起的底部切應力最大值分別為1.885 和0.293 N/m2，平均值分別為0.652 和0.097 N/m2，冬季潮流和波浪引起的底部切應力明顯大于夏季。

圖7 潮灘中上部光灘冬季近底部水動力時間序列變化過程a：淹沒水深，b：近底部流速，c：有效波高，d：底部切應力，e：近底部懸沙濃度。Fig.7 Time-series of near-bottom hydrodynamic processes in winter on the mid-upper tidal flat a: Inundation height, b: near-bottom mean current speed, c: significant wave height, d: bottom shear stress induced by tidal current and wave, e: near-bottom suspended sediment concentration.

3.4 潮間帶近底部懸沙濃度變化

懸沙濃度計算結果顯示，瑯岐島潮間帶中上部L04 站位近底部懸沙濃度呈現出顯著的時間變化特征（圖7e、8e）。在季節尺度上，冬季觀測期間近底部懸沙濃度為17.4 ～844.5 mg/L，平均值為306.8 mg/L，而夏季觀測期間近底部懸沙濃度為3.4 ～607.8 mg/L，平均值為138.5 mg/L，冬季懸沙濃度顯著高于夏季。在大-小潮周期尺度上，小潮期間近底部懸沙濃度相對較低，大潮期間相對較高，但潮周期平均值相差不大，如冬季大潮和小潮觀測期間近底部懸沙濃度平均值分別為296.5 和268.6 mg/L，而夏季大潮和小潮觀測期間近底部懸沙濃度平均值分別為158.1 和123.7 mg/L。在潮周期尺度上，漲潮初期和落潮后期相對較高，高平潮期間懸沙濃度相對較低。

3.5 潮間帶近底部懸沙輸運

近底部懸沙東向水平輸運率計算結果呈現出顯著的時間變化特征（圖9），冬季觀測期間近底部懸沙在東向上水平輸運率為-0.118 ～0.051 kg·m-2·s-1，平均值為0.011 kg·m-2·s-1；夏季觀測期間近底部懸沙水平輸運率為-0.006 ～0.009 kg·m-2·s-1，平均值為0.001 kg·m-2·s-1。

每個潮周期內的東-西方向上的凈懸沙輸運量計算結果顯示（圖9），瑯岐島潮間帶中部潮周期內近底部懸沙水平凈輸運通量呈現出明顯的時間變化特征。在大-小潮時間尺度上，小潮期間潮周期內懸沙凈水平輸運通量顯著小于大潮，而凈輸運方向則呈現出季節性差異：小潮期間冬季呈現出向岸凈輸運和向海凈輸運交替出現，而夏季則整體呈現向海凈輸運的變化特征；中潮期間冬季總體表現為向海凈輸運，而夏季則整體呈現向陸凈輸運的變化特征；大潮期間冬季總體呈現向岸凈輸運，而夏季則呈現向岸凈輸運和向海凈輸運交替出現的變化特征。在季節尺度上，冬季潮周期內懸沙凈水平輸運通量為-124.30～243.24 kg/m，觀測期間凈輸運通量為289.49 kg/m；夏季周期內懸沙凈水平輸運通量僅為-38.07 ～4.13 kg/m，觀測期間凈輸運通量僅為-58.73 kg/m，即冬季潮周期內凈水平輸運通量顯著大于夏季，并且冬季懸沙凈向海輸運，而夏季凈向陸輸運。

4 討論

河口潮灘是由地質過程以及水動力過程、沉積物供應和生態過程的相互作用形成的，其地貌形態及沉積物組成受動力、物源及生物活動控制，具有高度的時空變化特征[4,50]。大量的河流入海泥沙進入河口后，在水動力的作用下發育了廣泛的潮灘[6,51]，但由于我國河流泥沙入海主要集中在夏季，冬季入海泥沙通量很小[41]，因此，潮灘表層沉積物也表現出顯著的季節差異[52-54]。閩江河流泥沙入海通量也主要集中在夏季[55]，本文斷面表層沉積物分析結果（圖3、圖4）和沖淤觀測結果（圖6）也顯示，夏季期間河口潮灘捕獲了大量的河流入海沉積物，潮灘為淤積狀態，潮灘表層沉積物以粉砂和黏土為主，平均粒徑較??；冬季期間河流入海泥沙很少，河口潮灘很少能捕獲到來自河流的泥沙，潮灘表層沉積物在水動力的作用下不斷發生改造，最終細顆粒沉積物被搬運帶走，潮灘發生侵蝕，表層沉積物以砂和粉砂為主。關于沉積物供應變化對閩江口潮灘沉積物組成的影響綜合分析詳見李海琪等[56]，本文不再贅述，但需要指出的是，隨著閩江入海泥沙通量的顯著減少[55]，河流輸沙季節變化對瑯岐島潮灘季節性沖淤變化的影響將會減小。本文將主要從水動力方面闡述閩江河口潮灘沖淤的季節性變化機制。

影響潮灘沉積和地貌演化的動力過程主要包括潮汐、波浪、跨岸/沿岸流、風致環流、潮不對稱等[50,57]，尤其是在極淺水環境下，潮流、波浪是影響潮灘沉積物分布的主要動力因素[27-28]。沉積物是否發生侵蝕和淤積，主要是看水動力引起的底部切應力是否會超過臨界切應力，當水動力引起的底部切應力大于沉積物的臨界侵蝕切應力，則沉積物發生再懸浮，灘面呈現侵蝕狀態；當水動力引起的底部切應力小于泥沙臨界沉降切應力，則泥沙發生沉降，灘面呈現淤積狀態。根據本文觀測期間采集的表層沉積物粒度分析結果，冬季表層沉積物砂含量達到了93.9%，中值粒徑為2.85 Φ，屬于典型的非粘性沉積物，而夏季表層沉積物砂、粉砂和黏土含量分別為24.4%、54.2%和21.4%，中值粒徑為6.17 Φ，屬于典型的粘性沉積物。因此，根據相應的計算公式可知，觀測站位冬季和夏季表層沉積物臨界侵蝕切應力分別為0.158 和0.107 N/m2。對比結果發現（圖7d、8d），瑯岐島潮灘中上部地區近底層流速相對較小，由此產生的底部切應力也相對較小，對沉積物的活動性影響不明顯，但冬季觀測期間波浪引起的底部切應力顯著大于表層沉積物的臨界侵蝕切應力，表明冬季觀測期間觀測站位附近灘面發生強烈的侵蝕；夏季觀測期間波浪引起的底部切應力的變化大部分時刻均小于表層沉積物臨界侵蝕切應力，表明夏季觀測期間觀測站位附近灘面表層沉積物不易發生侵蝕。在長江口及江蘇沿海潮間帶觀測結果也顯示，在這種極淺水環境下波浪作用顯著影響潮灘沖淤過程[28,58]。研究表明[59]，懸沙臨界沉 降 切 應 力（τcrd）一 般 為0.06～0.1 N/m2，根 據Whitehouse 等的推薦，本文取τcrd=0.08 N/m2（其值適用于細顆粒泥沙），則對比結果顯示（圖8d），夏季觀測期間大部分時刻由潮流和波浪引起的底部切應力小于臨界沉降切應力，表明觀測期間觀測站位附近水體懸沙易發生沉降，導致灘面出現淤積。近底部懸沙凈輸運通量計算結果顯示，冬季觀測期間泥沙凈向海輸運，而夏季則凈向岸輸運。該結果表明，波浪是引起潮間帶表層沉積物發生侵蝕的主要動力，而潮流則是搬運泥沙的主要動力。漲落潮不對稱現象是河口海岸地區重要的水動力特征，是影響區域泥沙不對稱輸運的重要因素[60]，在河口海岸地區泥沙輸運和地貌演變過程中扮演著十分重要的角色[57,61]。由圖8 和圖9 可知，冬季漲潮期間懸沙濃度總體大于落潮，但由于波浪引起的底部切應力顯著大于懸沙臨界沉降切應力，導致懸沙很少發生沉降，漲潮流帶來的懸沙除部分被潮間帶上部的互花米草鹽沼植被捕獲外，其他懸沙又在落潮流的作用下向海輸運，而研究區無論冬季還是夏季，潮流均以落潮流占優，因此，落潮期間的底層懸沙輸運通量總體大于漲潮，潮周期內表現為凈向海輸運；夏季，雖然是落潮流占優，但漲潮期間的底層懸沙濃度總體大于落潮，僅在大潮期間波浪作用較強時期出現落潮懸沙濃度大于漲潮的情況，因此漲潮期間的懸沙輸運通量總體大于落潮，并且在多數時刻底部切應力小于臨界沉降切應力，懸沙沿途發生沉降。該結論得到了本文現場沖淤觀測結果的支持（圖6），同時也與其他河口地區潮灘沖淤過程一致[52,62]。

圖8 潮灘中上部光灘夏季近底部水動力時間序列變化過程a：淹沒水深，b：近底部流速，c：有效波高，d：底部切應力，e：近底部懸沙濃度。Fig.8 Time-series of near-bottom hydrodynamic processes in summer on the mid-upper tidal flat a: Iinundation height, b: near-bottom mean current speed, c: significant wave height, d: bottom shear stress induced by tidal current and wave, e: near-bottom suspended sediment concentration.

圖9 近底部瞬時懸沙水平輸運率及潮周期內懸沙水平輸運凈通量a：冬季，b：夏季。正值為向海輸運，負值為向岸輸運；T1—T19 為觀測期間的潮周期數。Fig.9 Seaward transport rate and net flux of near-bottom instantaneous suspended sediment in a tidal cycle in winter (a) and summer (b)Positive values indicate seaward transport, negative ones are landward transport; T1—T19: the number of tidal cycle during the observation.

本文觀測結果顯示，瑯岐島潮間帶中部地區，冬季觀測期間的潮流和波浪作用（圖7b、c）均明顯較夏季強（圖8b、c），這主要一方面是由于研究區冬季主要受東北季風控制，而夏季主要受西南季風控制，冬季期間風速大于夏季，由此引起研究區冬季期間潮流和波浪作用較夏季強[43,63]；另一方面，則與潮間帶地貌的季節性變化有關?，F場調查發現，瑯岐島潮間帶中部夏季被海三棱藨草覆蓋，而冬季則變為光灘?，F場觀測和數值模擬結果均表明，鹽沼植被具有顯著的弱流、消波作用[64-66]?，樶獚u潮灘夏季由于生長了大面積的海三棱藨草，起到了一定的緩流作用，因此，夏季觀測期間最大潮流流速和平均潮流流速均小于沒有鹽沼植被覆蓋的冬季，波浪作用也是夏季顯著小于冬季，表明瑯岐島潮間帶海三棱藨草的生長具有一定弱流和顯著的消波作用。根據不同站位沖淤變化對比分析（圖6），在觀測期間互花米草鹽沼邊緣灘面沖淤變化幅度很小，并且與其他兩個站位的沖淤變化趨勢沒有明顯的一致性，說明該地區灘面沖淤主要受水動力和互花米草的影響。位于海三棱藨草覆蓋區域和光灘區域灘面沖淤變化總體表現為相反的變化趨勢，這可能與海三棱藨草生長和消亡會引起潮間帶中部地區沖淤變化有關，即在有海三棱藨草覆蓋時，在動力-植被-泥沙的耦合作用下，漲潮帶來的泥沙被植被捕獲并沉積下來[30]，退潮水體懸沙濃度降低，導致在光灘地區退潮水流攜沙能力有所增強，在一定程度上可能會引起潮間帶下部光灘的侵蝕；當沒有植被覆蓋時，潮間帶中部遭受顯著侵蝕，侵蝕下來的泥沙在被落潮流攜帶至河口的過程中沿途發生沉降，引起潮間帶下部灘面淤積。上述分析表明，鹽沼植被在潮灘沖淤變化過程和潮灘發育演變中發揮著重要作用。

綜上所述，瑯岐島潮間帶中上部潮灘沖淤變化的動力學機制為：冬季，由于潮間帶中部的海三棱藨草消失，潮間帶被水淹沒期間的波浪作用強，由此引起的底部切應力顯著大于表層沉積物的臨界侵蝕切應力，表層沉積物發生再懸浮，再懸浮起來的泥沙在潮流的作用下向海輸運，導致灘面發生侵蝕；夏季潮間帶中部被海三棱藨草覆蓋，潮流和波浪均減弱，由此引起的底部切應力在大部分時刻小于表層沉積物的臨界侵蝕切應力和臨界沉降切應力，水體中的懸沙易發生沉降，導致灘面發生淤積。

5 結論

（1）瑯岐島潮灘表層沉積物粒度組分以粉砂和砂為主，黏土含量相對較少，呈現出顯著的時空差異，夏季主要以粉砂和黏土組成的細顆粒沉積物為主，冬季則主要以砂和粉砂組成的粗顆粒沉積物為主。

（2）瑯岐島潮間帶上部的互花米草鹽沼總體以持續淤積為主，潮間帶中部和下部表現出周期性的沖淤特征，總體以夏季淤積、冬季侵蝕為總特征，并且潮間帶中部的沖淤變化幅度較潮間帶下部大。

（3）沉積動力觀測與計算結果表明，瑯岐島潮間帶中上部近底部潮流流速相對較小，但波浪作用強，近底部懸沙濃度相對較高，其中冬季潮流流速、波浪作用及懸沙濃度均大于夏季，冬季潮周期內近底部懸沙以凈向海輸運為主，而夏季則以凈向岸輸運為主。

（4）瑯岐島潮灘沖淤變化主要受動力過程控制，冬季因潮間帶海三棱藨草消失，波浪作用強，導致表層沉積物易發生再懸浮，再懸浮的泥沙在較強的潮流作用下凈向海輸運；夏季因潮間帶有海三棱藨草生長，波浪作用較弱，水體中的懸沙易發生沉降，較弱的潮流將泥沙凈向岸輸運，導致潮間帶發生淤積。

致謝：黃書仁、黃思添、余永澤、陳?；?、劉三善參與了野外采樣工作，黃書仁參與了實驗室樣品分析，謹致謝忱！