章馨謠，戴志軍*，馮凌旋，黃祖明

(1.華東師范大學(xué) 河口海岸學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200241；2.上海市水利工程設(shè)計(jì)研究院有限公司,上海 200061)

潮流沖刷槽作為河口區(qū)域重要的地貌類型[1],是漲潮流或落潮流強(qiáng)烈沖刷作用下形成和維持的負(fù)地形,在世界各有潮河口普遍存在,如英國(guó)泰晤士河口[2],阿根廷布蘭加河口[3],新西蘭Manukau河口[4],我國(guó)長(zhǎng)江口[5]、珠江口[6]及杭州灣北岸[7-9]。相對(duì)穩(wěn)定的潮流沖刷槽給沿海城市或港口的發(fā)展提供了較好的水運(yùn)資源,但徑潮流復(fù)合作用以及人類活動(dòng)如沿岸圍墾、航道工程及建設(shè)堤壩等[7,10-11]亦在不同程度改變潮流沖刷槽的邊界形態(tài),進(jìn)而對(duì)河床及碼頭安全產(chǎn)生影響。顯然,研究沖刷槽地貌及沖淤過(guò)程對(duì)河口或海岸港口選址、航道維護(hù)和治理等都具有重大意義。

杭州灣是典型的喇叭型強(qiáng)潮河口灣,平均水深8～10 m。海灣北岸物質(zhì)組成普遍較細(xì),以粉砂質(zhì)黏土為主[12-14]。受控于強(qiáng)勁漲潮流塑造作用,杭州灣北岸疏松第四紀(jì)沉積層發(fā)育有金山、全公亭及乍浦等一系列深槽[8]。其中位于金山區(qū)上海石油化工廠附近海域的金山深槽是在漲潮流作用下形成的河口漲潮沖刷坑[5],為杭州灣北岸系列深槽的東部起點(diǎn)[9]。深槽呈NE—SW 向延伸,東部始于大、小金山之間,西至上海石化廠前沿(圖1)。深槽在潮流與長(zhǎng)江水沙變化等多種因素的影響下一直處在動(dòng)態(tài)變化過(guò)程。先前不少學(xué)者對(duì)金山深槽的形成及演變進(jìn)行研究,陳吉余等[5]率先提出潮流沖刷槽的概念,稱其為“河槽中的河槽”,并指出金山到澉浦一帶貼近北岸發(fā)育有一個(gè)規(guī)模巨大的沖刷槽和沖刷坑。之后,部分學(xué)者對(duì)金山深槽的地貌特征[15]、潮流特征[16]、沉積特征[17]、演變模式及其穩(wěn)定性[9,18-22]等作了探討。然而,隨著上海石油化工廠、港口碼頭、濱海旅游區(qū)等一系列重要沿海工程沿深槽所在位置的北部海岸構(gòu)建,特別在2017年環(huán)杭州灣大灣區(qū)建設(shè)推進(jìn),杭州灣北岸經(jīng)濟(jì)發(fā)展與建設(shè)出現(xiàn)新的發(fā)展,這將嚴(yán)重影響深槽地貌的變化狀態(tài)。深槽的快速?zèng)_淤亦將危及杭州灣北岸海堤穩(wěn)定,進(jìn)而對(duì)該地區(qū)港航及石化重工企業(yè)造成風(fēng)險(xiǎn)。開(kāi)展金山深槽岸段動(dòng)力地貌沖淤過(guò)程的研究,其價(jià)值不言而喻。據(jù)此,本文利用金山深槽岸段1989—2018年的高時(shí)間分辨率水深地形數(shù)據(jù),對(duì)水沙條件變化下的金山深槽動(dòng)態(tài)變遷過(guò)程進(jìn)行分析,進(jìn)而揭示金山深槽動(dòng)力地貌變化機(jī)制,探究深槽未來(lái)發(fā)展態(tài)勢(shì),以期為沿岸工程規(guī)劃和海灣經(jīng)濟(jì)發(fā)展提供科學(xué)依據(jù)。

圖1 金山咀及金山深槽位置Fig.1 The location of Jinshanzui and Jinshan Trough

1 數(shù)據(jù)來(lái)源與研究方法

本文的數(shù)據(jù)來(lái)自1989—2011年多幅杭州灣金山咀附近海圖地形資料:1989年(交通部上海海上安全監(jiān)督局,1∶100 000)、1997年(中華人民共和國(guó)交通部安全監(jiān)督局,1∶100 000)、2002年(中華人民共和國(guó)海事局,1∶25 000)、2007 年(中華人民共和國(guó)海事局,1∶25 000)、2011 年(中華人民共和國(guó)海事局,1∶25 000)以及金山咀附近2005年、2006年、2008年、2009年、2010年、2012年、2013年、2015年、2018年共9次實(shí)測(cè)斷面水深數(shù)據(jù)(圖1)。斷面水深數(shù)據(jù)采用單波束測(cè)深儀進(jìn)行采集,其主要特點(diǎn)是利用單點(diǎn)連續(xù)測(cè)量形成斷面水深數(shù)據(jù),結(jié)合多個(gè)斷面數(shù)據(jù)達(dá)到由點(diǎn)到面的區(qū)域地形數(shù)據(jù)測(cè)量[23]。同時(shí),本文還分析了長(zhǎng)江入?？刂普敬笸ㄕ?989—2018年連續(xù)的水沙變化數(shù)據(jù),數(shù)據(jù)來(lái)自水利部長(zhǎng)沙水利委員會(huì)發(fā)布的長(zhǎng)江泥沙公報(bào)(http:∥www.cjw.gov.cn/zwzc/bmgb/)。

為系統(tǒng)闡述金山深槽岸段的動(dòng)力地貌演化過(guò)程,本文將深槽的研究范圍定義為20 m 等深線包絡(luò)區(qū)域,定義深潭為30 m 等深線包絡(luò)區(qū)域。首先,利用ArcGIS平臺(tái)對(duì)海圖資料進(jìn)行數(shù)字化,將不同時(shí)期的海圖地形數(shù)據(jù)和實(shí)測(cè)斷面水深數(shù)據(jù)校正到統(tǒng)一的WGS 1984 坐標(biāo)系與理論深度基準(zhǔn)面。然后,采用Natural Neighbour插值方法構(gòu)建金山深槽及其附近岸段不同年份DEM 模型,同時(shí)利用Contour工具提取20和30 m 等深線。最后,利用Minus工具制作沖淤變化圖,Cut Fill模塊對(duì)相同的研究區(qū)域進(jìn)行處理,獲取DEM 間的體積變化,計(jì)算沖淤量分析金山深槽及其鄰近海床沖淤特征。進(jìn)一步利用MeanCenter工具以水深為權(quán)重計(jì)算深槽重心[24]分析深槽變遷狀態(tài)。此外,量算不同時(shí)間20 m 等深線的包絡(luò)范圍,分析金山深槽面積及體積變化特征。

2 結(jié)果分析

2.1 金山深槽及鄰近區(qū)海床沖淤特征

因金山深槽變遷頻繁,為更好地反映金山深槽的海床沖淤過(guò)程,對(duì)金山咀附近覆蓋深槽的區(qū)域及鄰近區(qū)海床進(jìn)行沖淤分析。1989—2018年,金山咀附近海床沖淤變化明顯,整個(gè)區(qū)域被沖刷近162.7×106m3,平均侵蝕厚度超過(guò)2.4 m(圖2)。其中,金山深槽內(nèi)大部分區(qū)域沖刷嚴(yán)重,局部最大沖刷厚度超過(guò)12 m。同時(shí),即使深潭(30 m 等深線包絡(luò)范圍)區(qū)域大部分淤積,但在深槽北部靠近海岸邊坡的20和30 m 等深線之間發(fā)生大面積沖刷,沖刷厚度為8～12 m。

圖2 1989—2018年金山咀附近沖淤變化Fig.2 Bathymetric changes of the Jinshanzui during 1989 to 2018

1989—1997年,金山咀附近整體淤漲,淤積厚度約0.92 m,年均淤積體積約為7.7×106m3,尤其是金山深槽西北部以及中部向東南延伸區(qū)域有兩處達(dá)到3 m 以上的淤漲區(qū),東南部淤漲區(qū)斜向貫穿金山深槽。此外,深槽內(nèi)部存在多處中心沖刷厚度在3 m 以上較為分散的沖刷體于深槽內(nèi)部呈串珠狀分布(圖3a)。

1997—2007年期間,金山咀附近海床總沖刷近2 m,年均侵蝕體積超過(guò)13.2×106m3(圖4),深槽內(nèi)部出現(xiàn)一條明顯沖刷帶(圖3)。其中1997—2002年深槽中下部初步形成兩處沖刷區(qū)(圖3b)。2002年后深槽中下部沖刷體連接形成沖刷帶。此后沖刷帶開(kāi)始向岸移動(dòng)(圖3c～圖3d)。2005年后海岸邊坡附近由淤積轉(zhuǎn)為沖刷,深槽東部、中部、西部存在3處沖刷深度超過(guò)3 m 的沖刷區(qū)(圖3d)。以2005年為界,在此年份之前深槽與岸線相連的北部為淤積區(qū)域,深槽以南為沖刷區(qū)域。此后深槽由淤轉(zhuǎn)沖,且深槽與岸堤之間的海床整體也由淤積變?yōu)榍治g,深槽外向海一側(cè)由侵蝕轉(zhuǎn)變?yōu)橛贊q。

2007—2018年,金山咀海床仍處于大沖大淤狀態(tài),不同年份地形資料表明年沖淤深度達(dá)到1 m 以上(圖4)。其中2007—2009年金山咀附近除城市沙灘沿岸出現(xiàn)條帶狀淤積外,海床大幅度沖刷,金山深槽以南全區(qū)域侵蝕,深槽內(nèi)部侵蝕嚴(yán)重,多處沖刷超過(guò)3 m 且僅殘留直徑小于1 km 大小不均一的淤積體(圖3e)。在2009—2011年,先前淤積的城市沙灘沿線目前變?yōu)榍治g,而先前侵蝕的區(qū)域都出現(xiàn)淤積,年均淤積體積達(dá)到68.9×106m3(圖3f)。2007—2011年雖逐年沖淤變動(dòng)幅度較大,但由于“沖刷-淤積”的交替變化,4年間總體沖淤變化幅度微弱,年均淤積厚度為0.035 m。2011年后金山咀附近經(jīng)歷了“沖刷-淤積-沖刷”演變,2011—2013年沖淤變化幅度較大(圖3g),其中2011—2013年沖刷深度＞1.5 m,2013—2015年期間發(fā)生補(bǔ)償性淤漲,淤積厚度超過(guò)1 m。但深槽內(nèi)部以沖刷為主,2015年后沖刷程度明顯減弱,年均沖刷約0.24 m(圖3h和圖3i)。

圖3 1989—2018年期間沖淤變化Fig.3 Bathymetric changes of the Jinshanzui during 1989－2018

圖4 1989—2018年金山咀附近年平均沖淤變化Fig.4 The yearly changes in accretion and erosion near the Jinshanzui during 1989－2018

為進(jìn)一步分析金山咀附近的水深變化,統(tǒng)計(jì)了20及30 m 以淺所測(cè)水深點(diǎn)出現(xiàn)頻率分布(圖5)可知,20 m 以淺水深點(diǎn)的變化趨勢(shì)與沖淤變化趨勢(shì)類似,1997年20 m 以上水深點(diǎn)頻率下降的同時(shí)金山咀附近海床整體淤積。1997—2006年20 m 以上水深點(diǎn)占比逐步上升,此后2007—2014年期間20 m 以上水深點(diǎn)的數(shù)量表現(xiàn)出兩次“增加-減小”的周期性變化,符合在其期間內(nèi)金山咀附近“沖深-淤淺”的旋回變化特征。然而2009、2012、2015和2018年20 m 以深的水深點(diǎn)數(shù)量明顯高于其他年份,和先前年份的地形沖淤比較,表現(xiàn)為強(qiáng)侵蝕特征。可表征深潭的30 m 以上水深點(diǎn)數(shù)量變化幅度更大,呈現(xiàn)出反復(fù)增減的變化狀態(tài),顯然深潭區(qū)域沖淤變化明顯,海床穩(wěn)定性差。2015年后超過(guò)30 m 的水深點(diǎn)數(shù)量趨于穩(wěn)定,同樣,在此期間金山咀附近發(fā)生持續(xù)性沖刷但沖刷強(qiáng)度較弱。1989—2018年期間,金山咀附近水深超過(guò)20 m的區(qū)域呈增長(zhǎng)趨勢(shì)。

圖5 金山深槽20 m 和30 m 以上水深點(diǎn)出現(xiàn)頻率Fig.5 The frequencies of water depth points above 20 m and 30m within the Jinshan Through

2.2 金山深槽變遷特征

分析金山深槽不同時(shí)期水深地形圖(圖6)發(fā)現(xiàn),深槽區(qū)域地貌形態(tài)發(fā)生了明顯變化。深槽在東西方向上變化較小,在南北方向上延伸明顯,向岸延伸0.5～1.3 km,向南延展0.4～1.2 km,尤其是深槽東部向北擴(kuò)張約1.1 km。同時(shí),1989年深槽內(nèi)部存在東部、中部及西部三處水深超過(guò)30 m 的深潭,中部深潭呈兩端寬、中間窄的“∞”形,深潭中心部分深度最深(圖6a),1997年中部深潭消失(圖6b)。2002年30 m 以上的深潭區(qū)域明顯增多且深度加深,西部與中部開(kāi)始連成一片且向岸移動(dòng)明顯(圖6c)。2002年后中西部深潭開(kāi)始經(jīng)歷“連接-分離-連接”的反復(fù)演變,東部深潭開(kāi)始淤積變淺不斷減小,2009年基本消失(圖6d和圖6e)。2011年由于強(qiáng)烈的沖刷作用,東部深潭沖刷回到2007年相似形態(tài)(圖6d和圖6f)。2012年金山附近岸段大范圍淤積,導(dǎo)致東部深潭再次淤漲,地貌形態(tài)與2009年形似(圖6e和圖6g)。此后金山深槽以輕微沖刷為主且沖刷速率有所下降,深潭面積略有增大(圖6h～圖6g)。

圖6 金山深槽不同時(shí)期水深地形Fig.6 Bathymetric maps in the Jinshan Trough from 1989 to 2018

由深槽內(nèi)部深泓線擺動(dòng)(圖7a)特點(diǎn)可知,1989—2018年期間深泓線中部無(wú)明顯變化,且基本保持穩(wěn)定狀態(tài)。深槽主泓線貼近北岸發(fā)育,深泓線擺動(dòng)主要發(fā)生在東西兩側(cè)。1997、2012和2018年深泓線有比較明顯的擺動(dòng)。1989—1997年上海石化工廠區(qū)前深泓線的擺動(dòng)幅度最大,擺幅超過(guò)1 km,2002年后該處深泓線回歸到原來(lái)位置。2012年深槽東側(cè)深泓線明顯向岸移動(dòng),隨后同樣擺動(dòng)回原位。2018年城市沙灘東側(cè)前沿深泓線向南移動(dòng)超過(guò)600 m。長(zhǎng)期以來(lái)深槽主軸擺幅小,較為穩(wěn)定。

以水深為權(quán)重,計(jì)算出近30 a深槽重心的經(jīng)緯度坐標(biāo),并利用ArcGIS 的標(biāo)準(zhǔn)差橢圓工具創(chuàng)建1989—2018年深槽重心的一級(jí)標(biāo)準(zhǔn)差橢圓(圖7b),分析1989—2018年金山深槽重心的移動(dòng)狀態(tài)。其標(biāo)準(zhǔn)差橢圓的方向?yàn)镹E 36°48′,說(shuō)明深槽重心整體呈NE—SW 方向分布。標(biāo)準(zhǔn)差橢圓長(zhǎng)軸長(zhǎng)約670 m,短軸為NE—SW 方向,長(zhǎng)約200 m。長(zhǎng)軸與短軸之比超過(guò)3∶1,意味著深槽重心在NE—SW 方向上的離散性大于NE—SW 方向,重心主要在NE—SW 方向上移動(dòng),所以重心的移動(dòng)方向與深泓線的擺動(dòng)方向具有一致性。

圖7 1989—2018年金山深槽深泓線及重心變化Fig.7 Changes in the thalweg migration and barycenter of the Jinshan Trough from 1989 to 2018

同時(shí),在1989—2002年,深槽重心經(jīng)度增大且有明顯東移趨勢(shì)。2002—2006年重心向西大幅度移動(dòng)到1989年重心以西位置,此后深槽重心一直在東西方向反復(fù)遷移。2010—2011年重心位置變化幅度最大,2011年重心移動(dòng)到最東側(cè),此后重心開(kāi)始向西移動(dòng),2013年基本回到2010年同一經(jīng)度,2013年后重心經(jīng)度有所增大?？傮w而言,1989—2018年重心經(jīng)度無(wú)明顯的變化規(guī)律,重心在東西方向上不斷擺動(dòng)(圖8a)。深槽重心緯度總體呈減小的趨勢(shì),深槽在南北方向上的移動(dòng)幅度大于東西方向(圖8b)。2005—2007 年、2010—2013年重心的經(jīng)緯度有類似的變化趨勢(shì),均經(jīng)歷了大幅的變化后回到原來(lái)的位置附近,其變化可能與此期間的沖淤波動(dòng)變化有關(guān)。整體上深槽重心的位置變化不大,一直在城市沙灘附近移動(dòng)。

圖8 1989—2017年金山深槽重心坐標(biāo)變化Fig.8 Variations in longitude coordinate and latitude coordinate of the Jinshan Trough barycenter from 1989 to 2018

2.3 金山深槽面積與體積變化過(guò)程

金山深槽沖淤動(dòng)態(tài)變化與深槽面積及體積的擴(kuò)大和縮小密切關(guān)聯(lián)。自1989年以來(lái),金山深槽內(nèi)部沖刷和淤積交迭發(fā)生,金山深槽的面積與體積也在不斷變化(圖9)。1989—2018年期間金山深槽20 m 等深線在不斷地變化移動(dòng),相應(yīng)的包絡(luò)面積有增有減,但整體上處于擴(kuò)大態(tài)勢(shì)(圖9a)；20 m 等深線以深的區(qū)域面積年均增長(zhǎng)率約為0.32 km2/a。具體而言,與1989年比較,1997年金山深槽面積銳減至不足10 km2,隨后保持波動(dòng)性擴(kuò)張趨勢(shì),在2009年深槽面積最大,之后面積有所減小,2014年后深槽面積再次增加,至2018年深槽面積達(dá)到20.8 km2,較1989年面積增長(zhǎng)了近1倍。

金山深槽面積變化表征其在水平空間上的變化狀態(tài),深槽體積變化則可以表現(xiàn)其在空間垂向上的沖淤變化(圖9b)。自1989—2018年深槽體積一直處于波動(dòng)狀,但總體處于擴(kuò)大狀態(tài),年均增長(zhǎng)速率超過(guò)1.9×106m3/a。其間,金山深槽體積由1997年87.9×106m3約縮小為1989年的一半,而在2002年體積增加至124.0×106m3,體積增加近3倍,2005年深槽體積減小至57.2×106m3。2005—2014年期間深槽體積反復(fù)增減變動(dòng)且變動(dòng)幅度較大,2008—2009年期間深槽體積增長(zhǎng)了67.8×106m3,而在2009—2010年深槽體積減少了64.8×106m3。類似地,2012年增長(zhǎng)了52.7×106m3,而又于2013年減小了51.7×106m3。2014年后深槽體積持續(xù)增長(zhǎng),年均增長(zhǎng)率達(dá)18.2×106m3/a。

圖9 1989—2018年金山深槽面積與體積變化狀態(tài)Fig.9 Yearly variations in area and volume of the Jinshan Trough from 1989 to 2018

3 討論

3.1 長(zhǎng)江來(lái)沙的影響

杭州灣緊靠長(zhǎng)江入海口的南面,是錢塘江河口的口外海濱,與錢塘江和長(zhǎng)江物質(zhì)交換頻繁。長(zhǎng)江的多年平均徑流量約為8 931億m3,而錢塘江的多年平均徑流量?jī)H為220.5億m3；長(zhǎng)江多年平均輸沙量為368 000億萬(wàn)t,而錢塘江的年均輸沙量?jī)H為長(zhǎng)江的0.7%。且錢塘江流域來(lái)沙主要堆積在浙江省海鹽縣澉浦鎮(zhèn)以上河段,部分由杭州灣南部下泄,僅有極少部分在潮流作用下沿杭州灣北岸東輸[25]。長(zhǎng)江入海泥沙在潮流作用下經(jīng)杭州灣北岸輸入杭州灣內(nèi),是杭州灣北岸主要的物質(zhì)來(lái)源。近些年,受堤防修筑及加固的影響,尤其在2003年三峽大壩開(kāi)始調(diào)蓄運(yùn)行后,長(zhǎng)江入海泥沙量急劇減少,致使長(zhǎng)江口進(jìn)入杭州灣北岸的泥沙減少,進(jìn)而影響北岸潮灘的發(fā)育過(guò)程[26-28],這從圖10a中長(zhǎng)江大通站年均輸沙量變化情況也可看出(圖10a)。

1989—2018年期間金山深槽的面積與體積整體呈現(xiàn)擴(kuò)大態(tài)勢(shì),由深槽面積及體積變化與大通站的年均輸沙量關(guān)系進(jìn)一步分析結(jié)果(圖10)可知,金山深槽的面積和體積變化與大通站年均輸沙量之間的關(guān)系密切。長(zhǎng)江來(lái)沙量減少深槽面積擴(kuò)張?bào)w積增大,反之,面積和體積都減小。相較于體積,金山深槽的面積對(duì)長(zhǎng)江來(lái)沙量變化的響應(yīng)更為強(qiáng)烈,金山深槽作為河口漲潮槽,其展寬作用大于刷深作用[29]。鄭璐等[22]利用2002和2012年的2期水深地形數(shù)據(jù)研究亦發(fā)現(xiàn),入海泥沙減少引起2002—2012年深槽整體處于沖刷狀態(tài)。而本文利用此期間內(nèi)的12期連續(xù)高分辨率水深數(shù)據(jù),還揭示了金山深槽及其附近海床在近十年內(nèi)為整體沖刷的狀態(tài),但各年間處于沖刷—淤積不斷交迭演化中(圖4)。這可能是因?yàn)榻鹕缴畈畚挥诤贾轂潮卑段鱾?cè),東側(cè)部分區(qū)域沖刷產(chǎn)生的泥沙進(jìn)入深槽附近岸段,在一定程度上緩解了由長(zhǎng)江入海泥沙減小對(duì)金山深槽體積變化的影響[21],此外,下文要討論的杭州灣北岸大規(guī)模周期性侵蝕—淤積波轉(zhuǎn)移對(duì)深槽地貌也會(huì)產(chǎn)生影響。

圖10 金山深槽體積變化與大通站年均輸沙量的關(guān)系Fig.10 Relationship between yearly volume changes of the Jinshan Trough and yearly sediment discharges at the Datong Station

3.2 杭州灣侵蝕-淤積波對(duì)金山深槽地貌的影響

杭州灣動(dòng)力地貌演變受到長(zhǎng)江入海水沙補(bǔ)給以及杭州灣沿岸水動(dòng)力條件的制約。在不同時(shí)段和不同岸段,沖淤動(dòng)態(tài)不同,這種沖刷和淤積作用在時(shí)間和空間上的組合,便形成了侵蝕－淤積波的波動(dòng)[30]。南匯咀向海一側(cè)存在一個(gè)沙尖外指的高位沙體(俗稱“搖頭沙”),其為侵蝕/淤積波的源頭?！皳u頭沙”水下部分受長(zhǎng)江口南槽及杭州灣兩股水流相互作用的影響而發(fā)生變動(dòng)。當(dāng)南槽水流作用強(qiáng)于杭州灣水流時(shí),其沙體向南突出誘發(fā)淤積波,南匯南灘在“搖頭沙”的掩護(hù)作用下淤漲,杭州灣北岸岸段在夾沙力增大的潮流作用下發(fā)生沖刷；當(dāng)杭州灣水流強(qiáng)時(shí),沙體向東或東南擴(kuò)展誘發(fā)侵蝕波,大量泥沙于“搖頭沙”堆積,南匯南灘沖刷后退,沖刷產(chǎn)生的泥沙在杭州北岸發(fā)生堆積發(fā)生淤積[30]。在侵蝕/淤積波的作用下杭州灣北岸岸灘此沖彼淤、相互消長(zhǎng),金山深槽及其附近海床沖淤變化頻繁。沿岸泥沙在強(qiáng)潮流作用下以波狀沙體向西移動(dòng)[13],將本文研究結(jié)果與前人對(duì)南匯動(dòng)力地貌的研究結(jié)果進(jìn)一步比較發(fā)現(xiàn)[31-32],當(dāng)南匯附近發(fā)生沖刷時(shí),金山深槽及其附近岸段相應(yīng)發(fā)生沖刷。進(jìn)一步比較不同年代沖淤圖中侵蝕體與淤積體的變化發(fā)現(xiàn),金山深槽內(nèi)部除2007—2009年全面沖刷(圖3e)以及2009—2011年全面淤積外(圖3f),侵蝕與淤積區(qū)域均成片交迭出現(xiàn),且侵蝕淤積幅度相對(duì)均勻。

在侵蝕/淤積波的影響下,金山深槽地貌形態(tài)也表現(xiàn)出波動(dòng)性的特點(diǎn),呈現(xiàn)出周期性的演化,并可以劃分為4個(gè)階段:階段1(1989—1997年),受深槽及其附近海床整體淤積的影響,深槽內(nèi)中部深潭淤積消失,深槽的各個(gè)方向均有所縮窄；階段2(1997—2002年),深槽內(nèi)沖刷加劇,于南北方向擴(kuò)展,中部深潭發(fā)育且與西側(cè)連接,東部深潭面積擴(kuò)大,深槽向岸運(yùn)移；階段3(2002—2013年),深槽再次發(fā)育增大,深槽內(nèi)部的地貌狀態(tài)經(jīng)歷了“沖深體積增大－淤積體積減小”反復(fù)循環(huán)變化,同時(shí)中西部深潭區(qū)域也在不斷經(jīng)歷局部填充分離以及再次沖刷連接的過(guò)程,東部深潭保持較為穩(wěn)定狀態(tài)；階段4(2013—2018年),深槽向西南方向移動(dòng),南北向展寬,中部深潭萎縮消失,西部深潭向中部延伸(圖11)。

圖11 金山深槽演變模式示意圖Fig.11 The geomorphologic cycle patterns of the Jinshan Trough

4 結(jié)論

基于1989年以來(lái)多年的海圖資料以及實(shí)測(cè)地形數(shù)據(jù),分析了金山深槽的沖淤變化與地貌演變特征,進(jìn)而探討了金山咀前沿河床演變及其穩(wěn)定性,得出以下結(jié)論:

①自1989年以來(lái)近30 a期間,金山深槽及其附近海床沖淤變化明顯,受侵蝕/淤積波影響呈沖刷、淤積周期性循環(huán)的地貌演變模式,最大沖淤深度超過(guò)10 m,整體以侵蝕為主。

②1989—2018年期間,金山深槽面積與體積整體上呈增長(zhǎng)趨勢(shì),其中面積和體積年均增長(zhǎng)率分別約為0.32 km2/a、1.94×106m3/a。長(zhǎng)江來(lái)沙急劇減少是導(dǎo)致金山深槽面積與體積擴(kuò)大的主要因素。

③金山深槽深泓線和重心在NW—SE方向上擺動(dòng),與深槽延伸方向一致,但擺幅較小。深槽在中部和西部深槽變動(dòng)較為活躍,中部和西部深潭多次經(jīng)歷“連接－分離－連接”的過(guò)程。深槽于南北方向上延伸明顯,從而引起深槽與堤岸之間的邊坡處于不穩(wěn)定狀態(tài),若金山深槽進(jìn)一步向岸移動(dòng),金山附近岸段的穩(wěn)定性則將面臨風(fēng)險(xiǎn)。