王立楊,付玉成,楊麗中,石洪源*,尤再進(jìn) ,馮興如
(1.魯東大學(xué) 土木工程學(xué)院,山東 煙臺(tái) 264025;2.大連海事大學(xué) 交通運(yùn)輸工程學(xué)院,遼寧 大連 116026;3.中國(guó)科學(xué)院海洋研究所 海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室,山東 青島 266071)
海灘是海陸交互的重要地帶,臺(tái)風(fēng)在海灘剖面的演變中起著重要作用[1],其引發(fā)風(fēng)暴的增水和大浪,在短時(shí)間內(nèi)會(huì)造成海灘的嚴(yán)重侵蝕和形態(tài)變化[2],海灘對(duì)臺(tái)風(fēng)的響應(yīng)與臺(tái)風(fēng)路徑、海岸形態(tài)、相對(duì)于海岸的風(fēng)向和波浪等許多因素有關(guān)[3]。臺(tái)風(fēng)強(qiáng)烈作用后的海灘需要長(zhǎng)時(shí)間的恢復(fù),短則數(shù)月,長(zhǎng)則數(shù)年[4],在臺(tái)風(fēng)后的1~2月內(nèi)是海灘的快速恢復(fù)期[5],之后要經(jīng)過(guò)長(zhǎng)達(dá)數(shù)年的調(diào)整,最終會(huì)恢復(fù)到其歷史平均形態(tài)[6]。
在全球氣候變暖大背景下,臺(tái)風(fēng)頻率有逐年增加的趨勢(shì)[7],海灘對(duì)臺(tái)風(fēng)的響應(yīng)愈加引起研究者的重視,Larson和Kraus[8]提出了風(fēng)暴侵蝕海岸模型并給出了經(jīng)驗(yàn)剖面數(shù)據(jù),Basco[9]通過(guò)對(duì)37個(gè)海灘的觀測(cè)分析發(fā)現(xiàn)臺(tái)風(fēng)對(duì)不同位置的海灘造成的侵蝕以及臺(tái)風(fēng)后的恢復(fù)也均不相同。國(guó)內(nèi)對(duì)海灘的風(fēng)暴響應(yīng)研究也很重視,Qi等[10]發(fā)現(xiàn)不同地貌類型的海灘對(duì)臺(tái)風(fēng)的響應(yīng)也不同,反射性越強(qiáng)的海灘,在臺(tái)風(fēng)中的變化越明顯;陳子燊[11]探究了弧形海岸對(duì)臺(tái)風(fēng)的響應(yīng)過(guò)程,提出海灘剖面的變化主要由灘肩進(jìn)退、后濱侵蝕、水下沙壩運(yùn)移等空間過(guò)程組成;蔡鋒等[12]探究了熱帶氣旋前進(jìn)方向兩側(cè)海灘風(fēng)暴響應(yīng)特征的差異,認(rèn)為位于臺(tái)風(fēng)前進(jìn)路徑右側(cè)的海灘對(duì)臺(tái)風(fēng)的響應(yīng)速度更快,遭受的侵蝕破壞程度更大;童宵嶺等[13]通過(guò)探究臺(tái)風(fēng)“??钡卿洉r(shí)浙江皇城海灘剖面的變化過(guò)程,認(rèn)為岬角突出程度、相對(duì)于臺(tái)風(fēng)的走向是海灘對(duì)臺(tái)風(fēng)響應(yīng)強(qiáng)度的重要影響因素。
但是目前針對(duì)同一海灘的不同區(qū)域?qū)ε_(tái)風(fēng)的響應(yīng)研究較少,而且關(guān)于海灘對(duì)臺(tái)風(fēng)的響應(yīng)研究,國(guó)內(nèi)主要聚焦于東南沿海地區(qū)的海灘,鮮有關(guān)注山東海灘。山東擁有3 000 km 長(zhǎng)的海岸線,而且山東岸灘侵蝕災(zāi)害頻發(fā),風(fēng)暴潮、巨浪等極端天氣造成了巨大的經(jīng)濟(jì)損失,僅“利奇馬”過(guò)境對(duì)山東省就造成了多達(dá)26.31億的經(jīng)濟(jì)損失[14]。因此臺(tái)風(fēng)對(duì)楮島海灘影響的研究,可為山東海灘保護(hù)研究提供重要參考。
本研究利用臺(tái)風(fēng)前后楮島海灘固定剖面的連續(xù)測(cè)量和沉積物數(shù)據(jù),通過(guò)研究同一沙灘不同剖面在臺(tái)風(fēng)下的變化特征,探究其對(duì)“利奇馬”的響應(yīng)特征,為楮島海灘的開發(fā)與養(yǎng)護(hù)提供科學(xué)依據(jù)。
研究區(qū)楮島(122°31′~122°34′E,37°00′~37°03′N)位于山東半島魯東膠南隆起區(qū)域,東臨黃海、南臨石島灣、北接桑溝灣(圖1)。楮島面積約0.75 km2,沙灘長(zhǎng)約2.7 km,前濱帶寬43.2 m,后濱寬33.4 m,沉積物主要以細(xì)砂為主,為優(yōu)質(zhì)的天然砂質(zhì)海岸。研究區(qū)屬于溫帶季風(fēng)氣候,常風(fēng)向?yàn)镹-NW,平均風(fēng)速6.4~6.6 m/s,年降水量分配不均,夏季降水量約占全年總降水量的80%[15]。楮島所處海域主要以風(fēng)浪為主,基本不受涌浪的影響[16],主要波向是SW 和SEE 向,年平均有效波高0.41 m,年最大有效波高1.47 m,年最大波高3.6 m,年平均周期5.02 s。
圖1 研究區(qū)位置及“利奇馬”臺(tái)風(fēng)路徑Fig.1 The location of the study area and the route of Typhoon Lekima
2019年第9 號(hào)臺(tái)風(fēng)“利奇馬”是當(dāng)年登錄我國(guó)的最強(qiáng)臺(tái)風(fēng),于2019-08-04T14:00在太平洋上生成,中心坐標(biāo)(131°54′E,17°24′N),中心最大風(fēng)力為18 m/s,強(qiáng)度為熱帶風(fēng)暴級(jí),8月7日23:00加強(qiáng)為超強(qiáng)臺(tái)風(fēng),8月10日01:45在臺(tái)州溫嶺城南鎮(zhèn)沿海登錄,登錄時(shí)風(fēng)速達(dá)到了16 級(jí)(52 m/s),登陸后的“利奇馬”穿過(guò)浙江和江蘇,8月11日12:00從連云港附近出海,沿著偏北方向移動(dòng),于當(dāng)天20時(shí)在青島市黃島區(qū)再次登陸,中心風(fēng)力為23 m/s,強(qiáng)度為熱帶風(fēng)暴級(jí),風(fēng)圈的影響半徑最大約330 km。自此楮島進(jìn)入其影響范圍內(nèi),直至被停止編號(hào)[14]。
通過(guò)楮島外中國(guó)科學(xué)院7 號(hào)浮標(biāo)(122°34′48″E,37°03′36″N)獲得“利奇馬”登陸前后研究區(qū)的波浪、風(fēng)數(shù)據(jù)(圖2),浮標(biāo)可收集波高、波向、風(fēng)速、風(fēng)向等波浪氣象要素。波動(dòng)時(shí)分辨率可達(dá)1 cm,海洋起伏的交叉靈敏度小于3%,數(shù)據(jù)以1.28 Hz的頻率記錄,每0.5 h處理一次。8月正常狀態(tài)下,研究區(qū)海域風(fēng)向以NW 向和WEW 向?yàn)橹?臺(tái)風(fēng)影響期間研究區(qū)海域風(fēng)向主要為S向和SSE 向(圖2),于8月11日16:00研究區(qū)最大風(fēng)速達(dá)到了18.6 m/s,是常風(fēng)速的3倍。臺(tái)風(fēng)影響期間研究區(qū)最大波高達(dá)到了5.7 m(圖2),主要浪向?yàn)镾E向(圖3)。距離研究區(qū)最近的成山頭潮位站(122°00′E,37°18′N)在風(fēng)暴過(guò)程中最大的風(fēng)暴增水為60 cm[14],達(dá)到了當(dāng)?shù)氐乃{(lán)色警戒潮位。
圖2 臺(tái)風(fēng)登陸前后最大波高逐時(shí)變化過(guò)程Fig.2 Rose chart of wave and wind directions in the study area under typhoon and normal conditions
圖3 臺(tái)風(fēng)和正常狀態(tài)下研究區(qū)波向及風(fēng)向玫瑰圖Fig.3 Variations of maximum wave height with time before and after the typhoon
為了獲取臺(tái)風(fēng)前后楮島海灘的地形變化情況,于2019-08-10—13使用測(cè)高儀(北斗海達(dá)TS7 RTK)對(duì)海灘進(jìn)行觀測(cè)。在楮島長(zhǎng)約2.7 km 的沙灘設(shè)置10個(gè)固定監(jiān)測(cè)剖面(圖4),岸邊硬化路面打固定樁,每次測(cè)量前進(jìn)行點(diǎn)放樣操作,確保測(cè)量的是固定剖面。自岸邊固定樁開始手持儀器,選擇在最低潮位時(shí)進(jìn)行測(cè)量,測(cè)量點(diǎn)之間距離設(shè)置為2 m,沖刷區(qū)則縮短測(cè)量距離,設(shè)置為50 cm。
圖4 剖面位置Fig.4 Section location
在海灘剖面的沖刷區(qū)、灘肩以及沙丘底腳位置3個(gè)固定點(diǎn)采集沙灘表面沉積物,取沙剖面為剖面N01、剖面N06和剖面N10號(hào)。取樣時(shí)先除去沙灘表面的垃圾等附著物,采集約200 g的沉積物,在實(shí)驗(yàn)室中對(duì)砂樣烘干,使用電動(dòng)振篩分儀對(duì)所取樣本篩分,篩分范圍為-2.25Ф~4Ф(Ф=-log2D,D為中值粒徑),并采用Folk&Ward公式計(jì)算沉積物粒度參數(shù):平均粒徑(Md)和分選系數(shù)(σi)[17]。
將各個(gè)剖面的歷次測(cè)量圖疊加起來(lái)(基于1985國(guó)家高程基準(zhǔn)),對(duì)比分析海灘剖面的地形變化特征,同時(shí),通過(guò)計(jì)算剖面的單寬蝕積量(UED,Unit Erosional Deformation)[10]來(lái)表示海灘侵蝕堆積情況(圖5):
圖5 臺(tái)風(fēng)作用下海灘剖面變化[10]Fig.5 Schematic diagram of beach profile changes under a typhoon[10]
式中:UED 表示剖面的侵蝕、淤積量;X為水平坐標(biāo);Z為剖面高程。
水動(dòng)力與地形的相互作用導(dǎo)致了海灘狀態(tài)的變化,海灘狀態(tài)所發(fā)生的改變可以清晰地反映出海灘在臺(tái)風(fēng)中的變化。海灘狀態(tài)參數(shù)的變化與波浪及泥沙粒徑有密切關(guān)聯(lián),Wright和Short[18]提出使用無(wú)量綱沉降速率(Ω)來(lái)劃分海灘的類型:
式中:Hb為破波波高(m);Ws為泥沙沉降速率(m/s);T為波浪周期(s)。
由于臺(tái)風(fēng)影響期間研究區(qū)域的破波波高數(shù)據(jù)難以獲取,可通過(guò)海灘坡度、沉積物粒徑等易得數(shù)據(jù)間接推導(dǎo)。Sunamura[19]通過(guò)因次分析得到海灘坡度計(jì)算公式:
式中:tana為海灘坡度;g為重力加速度;D為沉積物中值粒徑(mm)。
依照《海洋傾倒區(qū)選劃技術(shù)導(dǎo)則》[20]規(guī)定的泥沙沉速公式和研究區(qū)沉積物平均粒徑大小,選擇D為0.15~1.50 mm,利用岡恰諾夫過(guò)渡區(qū)沉降經(jīng)驗(yàn)公式計(jì)算沉降速率w(cm/s):
式中:ρs=2.65 g/cm3為沙的密度;ρ=1.02 g/cm3為海水密度;t為溫度(℃)。
由式(2)~(4)聯(lián)立可得[21]:
當(dāng)Ω<1時(shí),海灘坡度較陡趨于反射型;當(dāng)1≤Ω<6時(shí),海灘型態(tài)介于兩者之間為過(guò)渡型;當(dāng)Ω≥6時(shí),海灘坡度平緩趨于消散型。
3.1.1 海灘剖面變化
臺(tái)風(fēng)“利奇馬”明顯影響了楮島海灘,海灘剖面地形地貌都發(fā)生了改變。剖面N01~N04類型相似(圖6),4個(gè)剖面灘肩明顯,后濱寬度較大。臺(tái)風(fēng)后,剖面類型基本保持原狀或者略顯沙壩式斷面雛形,后濱部分略微侵蝕,灘肩蝕退量很小,侵蝕位置主要集中在高、中潮帶,水下部分略微下蝕,海灘后濱有堆積沙丘形成,單寬侵淤量相差不大,自西向東分別為-9.3,-16.4,-9.2和-10.6 m3/m。
圖6 剖面N01~剖面N04地形變化Fig.6 Topographic changes for N01-N04 profiles
剖面N05和N06(圖7)后濱發(fā)育有窄而尖的風(fēng)暴灘肩。臺(tái)風(fēng)過(guò)后,剖面N05后濱有堆積沙丘形成,灘肩蝕退;剖面N06后濱灘面略微堆積,灘肩被沖刷至基本消失。2個(gè)剖面前濱和潮間帶侵蝕嚴(yán)重,低潮帶形成了淤積與侵蝕相間的起伏地形,單寬淤積量自西向東分別為:-12.9和-17.7 m3/m。由圖7可知,海灘形態(tài)由灘肩剖面向沙壩剖面轉(zhuǎn)換在剖面N05和N06之間開始,并且在N05剖面以西,沙灘剖面形態(tài)相對(duì)穩(wěn)定。
圖7 剖面N05和剖面N06地形變化Fig.7 Topographic changes for N05-N06 profiles
海灘東側(cè)4個(gè)剖面(剖面N07~剖面N10)(圖8)對(duì)“利奇馬”的響應(yīng)非常強(qiáng)烈,所有剖面均發(fā)生了嚴(yán)重的侵蝕,灘面有大量的侵蝕陡坎發(fā)育,最大落差為0.4 m。高潮帶整體下蝕嚴(yán)重,前濱泥沙被海浪沖刷帶離致使灘面凹陷,向海移動(dòng)的泥沙在水下堆積成沙壩,剖面形態(tài)轉(zhuǎn)變?yōu)榈湫偷纳硥问狡拭妗?個(gè)剖面的單寬淤積量分別為:-19.8、-17.0,5.2和-7.0 m3/m,侵蝕量自西南至東北呈現(xiàn)逐漸減少的趨勢(shì)。
圖8 剖面N07~剖面N10地形變化Fig.8 Topographic changes for N07-N10 profiles
3.1.2 海灘沉積物變化
研究區(qū)正常狀態(tài)下和臺(tái)風(fēng)前后海灘表層沉積物特征見表1。由表1可知,在無(wú)臺(tái)風(fēng)月份,海灘沉積物粒徑波動(dòng)很小,臺(tái)風(fēng)期間,研究區(qū)海灘的沉積物發(fā)生了顯著的改變,沉積物粗化,沖刷區(qū)非常顯著,沙灘的平均粒徑范圍由臺(tái)風(fēng)前的2.02Ф~3.37Ф變?yōu)?.36Ф~2.43Ф。在臺(tái)風(fēng)前,各個(gè)剖面的沉積物分選系數(shù)集中在0.64~1.08,以分選系數(shù)較好(0.50<σi<0.71)到中等(0.71<σi<1.00)為主。經(jīng)臺(tái)風(fēng)作用分選系數(shù)范圍變?yōu)?.62~2.09,沉積物分選性變差。
表1 楮島海灘表層沉積物特征Table 1 Characteristics of surface sediments from Chudao beach
綜上所述,3個(gè)剖面沖刷區(qū)區(qū)域的粒徑變化均變粗,而且幅度最大,灘肩區(qū)域次之,沙丘底腳區(qū)域變化最小,分選系數(shù)(σi)的變化同平均粒徑Ф值表現(xiàn)較為一致。臺(tái)風(fēng)期間強(qiáng)大的動(dòng)力條件對(duì)灘面產(chǎn)生的侵蝕作用將表層細(xì)顆粒沉積物帶走,使下層粗顆粒沉積物裸露(圖9),是導(dǎo)致臺(tái)風(fēng)后粒度整體粗化的重要原因。
圖9 2019-08-13臺(tái)風(fēng)過(guò)后灘面情況Fig.9 The situation of the beach after the typhoon in 2019-08-13
3.1.3 海灘地形動(dòng)力狀態(tài)變化
在極端波況下隨著海灘對(duì)臺(tái)風(fēng)的響應(yīng),其地形動(dòng)力狀態(tài)也會(huì)發(fā)生變化。臺(tái)風(fēng)前后,研究區(qū)海灘剖面N01、剖面N06和剖面N10的海灘地形狀態(tài)特征值見表2。
表2 利奇馬前后楮島海灘剖面類型Table 2 Beach profile types of the Chudao before and after the Typhoon Lekima
臺(tái)風(fēng)過(guò)后,剖面N01、剖面N06和剖面N10的Ω 值都變大,地形狀態(tài)相比臺(tái)風(fēng)前更加消散。剖面N01的Ω值變化最小,剖面N10變大幅度最大,自西南向東北呈現(xiàn)逐漸增大的態(tài)勢(shì)。臺(tái)風(fēng)期間研究區(qū)的最大波高達(dá)到了5.7 m,研究區(qū)海洋動(dòng)力環(huán)境增強(qiáng),海灘迅速響應(yīng),原本屬于消散型的剖面N01、剖面N06和剖面N10呈現(xiàn)更加消散的狀態(tài),灘肩消失,坡度減小,以緩沖大浪帶來(lái)的能量。剖面N01左側(cè)有小型船塢的遮蔽削弱了波浪的強(qiáng)度以及有寬大的灘肩致使其受臺(tái)風(fēng)的影響相對(duì)較小。
在臺(tái)風(fēng)作用下,楮島海灘各個(gè)剖面均遭侵蝕,并且侵蝕主要集中于灘肩和潮間帶。前濱表層細(xì)砂被大浪帶走,致使灘面下的大量礫石露出,由于風(fēng)暴影響期間主要浪向?yàn)镾SE 向和S向,而沙灘為NNE-SSW 走向,臺(tái)風(fēng)造成的大浪幾乎直接沖擊所有剖面,致使臺(tái)風(fēng)對(duì)海灘的所有剖面均造成了一定程度的破壞。
根據(jù)尤再進(jìn)等[22]給出的沙丘-海灘穩(wěn)定性因素分析,海灘西南側(cè)的剖面N01~剖面N05由于灘肩較寬,沙丘底腳高程較高,坡度較緩,穩(wěn)定性相對(duì)較好,故在臺(tái)風(fēng)期間雖有侵蝕但蝕量較小,剖面形態(tài)變化相對(duì)不明顯。海灘東北側(cè)的剖面N06~剖面N10沒(méi)有灘肩、坡度陡、穩(wěn)定性差,故在臺(tái)風(fēng)期間侵蝕情況比較嚴(yán)重,剖面形態(tài)變化明顯,成為典型的風(fēng)暴型剖面。
在臺(tái)風(fēng)期間,研究區(qū)風(fēng)浪作用以南向?yàn)橹?在海灘延伸方向上存在著由西南向東北的沿岸流。依照沿岸流計(jì)算公式[16]:
式中:αb為波峰線與岸線的夾角。
隨著海灘方位的改變,斜向波入射角αb變小,在αb變小的位置也存在沿岸流流速的減小梯度。臺(tái)風(fēng)過(guò)境中,高速的沿岸流將會(huì)在海岸拐點(diǎn)以南海區(qū)激浪帶和濺浪帶內(nèi)的泥沙攜帶,并且自西南向東北搬運(yùn)。然而隨著沿岸流在岸線方向改變處速度變緩,導(dǎo)致流體切應(yīng)力變小,因此粒徑較大的泥沙在岸線曲率最大點(diǎn)處沉降下來(lái),故剖面N10沖刷區(qū)粒徑變化最大,由2.10Ф變?yōu)?.85Ф,且分選由中等變?yōu)椴睢?/p>
海灘對(duì)風(fēng)暴的響應(yīng)與風(fēng)暴路徑的走向有密切的聯(lián)系,楮島海灘呈SW—NE 走向,灣口朝向偏東南,SE向浪對(duì)其地形變化影響最大,楮島海灘的剖面走向與臺(tái)風(fēng)“利奇馬”走向近乎垂直,SE 向的大浪近乎正面沖擊海灘,使得其剖面產(chǎn)生了較大的沖淤變化。熱帶氣旋的低層風(fēng)場(chǎng)是不對(duì)稱的,在前進(jìn)方向的右側(cè),由于風(fēng)向和前進(jìn)方向一致,移速和風(fēng)速相互疊加,這就造成熱帶氣旋前進(jìn)方向右側(cè)的風(fēng)力比較大,臺(tái)風(fēng)所誘發(fā)的臺(tái)風(fēng)浪能量較高。楮島海灘位于“利奇馬”前進(jìn)方向的右側(cè),并且其沿海為向岸浪,這就使得波浪的破壞作用加劇,導(dǎo)致研究區(qū)海灘形態(tài)變化十分顯著,海灘東北側(cè)的剖面N06~剖面N10變?yōu)樯硥问狡拭妗?/p>
本文通過(guò)對(duì)“利奇馬”登陸前后楮島海灘剖面的觀測(cè)和表層沉積物樣品研究,探討該處海灘對(duì)臺(tái)風(fēng)的響應(yīng)特征,得出以下主要結(jié)論:
①在臺(tái)風(fēng)“利奇馬”的影響下,楮島海灘遭受侵蝕,不同岸段對(duì)臺(tái)風(fēng)的響應(yīng)表現(xiàn)出差異性:西南側(cè)剖面N01~剖面N05形態(tài)變化不大,東北側(cè)剖面N06~剖面N10形態(tài)由灘肩式向沙壩式轉(zhuǎn)變。這種差異性與剖面形態(tài)有著密切關(guān)聯(lián)。由于在剖面N10斜向波入射角αb變小,導(dǎo)致沿岸流在此處速度變緩,粒徑較大的沉積物在此處沉降下來(lái),致使剖面N10沖刷區(qū)粒徑變化最大。
②在臺(tái)風(fēng)作用下,楮島海灘的表層沉積物變粗,沖刷區(qū)尤為明顯,臺(tái)風(fēng)后平均粒徑由2.02Ф~3.37Ф變?yōu)?.36Ф~2.43Ф。沉積物的分選性變差,分選系數(shù)由0.64~1.08變?yōu)?.62~2.09。海灘剖面的無(wú)量綱沉降速率(Ω)增大,表明海灘狀態(tài)朝著更加消散狀態(tài)轉(zhuǎn)化,以緩沖和適應(yīng)高強(qiáng)度的海洋動(dòng)力條件。