江蘇中部海岸潮溝形態對灘涂圍墾的響應

時海東，沈永明*，康敏

(1.南京師范大學 地理科學學院，江蘇 南京 210023)

?

江蘇中部海岸潮溝形態對灘涂圍墾的響應

時海東1，沈永明1*，康敏1

(1.南京師范大學 地理科學學院，江蘇 南京 210023)

摘要：根據2001-2013年江蘇中部海岸的遙感影像和野外觀測數據，對灘涂圍墾影響下潮溝的形態演變過程進行分析。結果顯示，研究區潮溝在大規模灘涂圍墾影響下整體有明顯的退化和消亡趨勢，潮溝集水面積由2001年的293.94 km2逐漸降低到2013年的103.40 km2，潮溝長度隨著其集水面積的縮減呈冪函數減小(R2=0.98)，而潮溝寬度隨著潮溝長度和集水面積的變小呈指數方式變窄(R2分別為0.87、0.93)，總體上研究區每圍墾1 km2灘涂，潮溝長度和寬度分別減小約2.60 km和0.67 m，在整個灘涂圍墾過程中，落潮后期潮灘的主要排水方式逐漸由潮溝排水向灘面排水方式轉變。江蘇中部海岸潮溝形態在統計上具有多尺度分形結構，總體上研究區分緯值從大到小依次為潮溝分維值、各潮溝系統分維值、單支潮溝分維值。研究區潮溝的分維值(無標度區間25～210內)在灘涂圍墾影響下呈波動性降低，由2001年的1.26逐漸降低到2013年的1.13，各潮溝系統的分維值隨著其分汊個數減少呈線性降低(R2=0.89)，各單支潮溝的分維值隨著其蜿蜒性減小呈線性降低(R2=0.86)。根據潮溝的退化和消亡過程反推潮溝的起源和發育過程，本研究認為研究區最早出現的一類潮溝可能是陸源水流侵蝕型潮溝，而灘面水流沖刷型和潮流輻聚侵蝕型潮溝是隨著潮灘淤積變寬過程逐漸發育的，其最初的發育形態可能是潮灘下部的沖刷凹地。

關鍵詞：潮溝；灘涂圍墾；地貌演變；江蘇中部海岸

1引言

作為潮灘上典型的沉積地貌單元，潮溝伴隨潮灘的形成而出現，在潮灘沉積和水流作用下逐漸發育成熟[1—2]。潮溝雖然在世界各地河口三角洲、平原海岸、海灣、潟湖以及沙洲淺灘上廣泛發育，但多數已處于均衡或準均衡狀態[3—5]，獲得其早期發育過程和后期消亡過程的實測數據較困難[6—7]。較多學者根據衛片、航片和野外觀測數據對潮溝形態進行描述[8—9]，分析潮溝發育中期的形態演變過程[10—11]，對影響潮溝發育的環境因子(如潮灘地形、沉積物粒徑、水動力、鹽沼植被、底棲動物、風暴潮、海平面變化等[12—18])進行探討，或根據概念模型[19—21]、數學模型[22—24]、實驗室模擬數據[25—27]對潮溝起源和早期發育過程進行研究，關于潮溝退化或消亡過程的報道相對較少。近年來，隨著各地灘涂圍墾活動的快速進行，潮灘上潮溝形態發生較大變化，尤其在沖淤較強的海岸地區，其整體上呈現退化和消亡趨勢，并在局部產生新的潮溝系統[28—32]，這為研究潮溝早期發育過程或后期消亡過程提供了新信息源。陳才俊[28]、黃海軍和樊輝[29]、Hood[30]、Symonds和Collins[31]、吳德力等[32]分析了灘涂圍墾前后潮溝的形態變化和發育特征，而大規模圍墾期間潮溝的形態演變過程有待進一步研究。

江蘇中部海岸是我國典型的淤長型粉砂淤泥質中、強潮平原海岸，寬闊的潮灘上潮溝系統廣泛發育，是研究我國海岸潮溝地貌特征的重要區域之一[13]，近年來隨著該區灘涂圍墾活動的持續進行，潮溝形態呈現出一定的退化和消亡趨勢。本文選擇江蘇中部從斗龍港至梁垛河口之間面積約為470 km2的潮灘為研究區，擬通過統計2000年以來各圍墾階段潮溝的長度、寬度、密度、集水面積、分形維數等參數信息，定量描述該區潮溝的退化和消亡過程，分析灘涂圍墾活動對潮溝形態變化的影響，并根據新形成的潮溝系統和野外觀測數據，對潮溝發育和形態演變過程進行探討。

2數據與方法

2.1數據來源與處理

本文所用數據主要來源于美國國家航空航天局NASA的陸地衛星Landsat TM/ETM+遙感影像，行列號為119/37，影像的空間和時間分辨率分別為30 m×30 m和16 d，成像時間見表1。并分別于2012年12月和2013年4月對江蘇大豐、如東的潮灘潮溝地貌在匡圍前后的形態變化特征進行觀測，以輔助影像的解譯和分析。所用影像統一采用橫軸墨卡托51°N帶投影(UTM 51N)，WGS84基準面，使用ENVI 4.7軟件對影像進行幾何精校正(精度控制在一個像元以內)。由于研究區潮溝線狀地物特征明顯，選用TM/ETM+5、TM/ETM+4和TM/ETM+3波段進行RGB標準假彩色合成，結合野外觀測在ArcGIS 9.3中進行常規的目視解譯[33]，對提取要素(圍墾區、潮溝、潮灘等)進行矢量化，統計不同時期各矢量要素的變化特征。

表1　主要遙感影像的成像時間

2.2研究方法

在描述潮溝平面形態的眾多參數中，潮溝長度、寬度、集水面積、潮溝密度和分形維數等對灘涂圍墾的響應較大，本文主要對這些參數在圍墾期間的變化特征進行了統計。其中，潮溝長度由潮溝中軸線法[34]進行描述和統計；考慮到潮位變化對潮溝口寬度的較大影響，潮溝寬度根據潮溝口最低部(喇叭口以上)位置潮溝主干的寬度進行計算；潮溝集水區域[35—36]基于研究區潮溝系統的分汊性和匯水特征來描述(分汊性：每個潮溝系統由一條主潮溝和一系列不同級別的支潮溝組成，主支潮溝相互連通呈樹枝狀分布；匯水特征：落潮后期灘面薄層水流向附近的潮溝，細潮溝里的水流向支潮溝，繼而匯入主潮溝)，如圖1所示，將潮溝口、各支潮溝的尾梢頂點及堤壩連線所圍區域表示為潮溝系統的集水區域，集水區域的邊界在向海一側緊鄰潮溝，向岸一側到達圍墾堤壩前沿，左右兩側則以支潮溝尾梢頂點的連線為準，集水區域的面積即為潮溝的集水面積；潮溝密度[37]是描述潮溝發育密集程度的參數，本文分別用潮溝長度與其集水面積之比(潮溝密度Ⅰ)和潮溝長度與潮灘面積之比(潮溝密度Ⅱ)來定義；潮溝分維值[38—41]是度量潮溝形態復雜性的參數，用數盒子方法進行統計。

圖1　2002年研究區的潮溝系統及其集水區域Fig.1　Tidal creek systems and their catchments of the study area in 2002

3結果分析

3.1潮溝長度對灘涂圍墾的響應

如圖2所示，隨著灘涂圍墾面積的增加，研究區潮溝長度呈減小趨勢，總體上平均每圍墾1 km2灘涂，潮溝長度減小約2.60 km。根據潮溝長度在圍墾期間的減小速率將灘涂圍墾過程大致分為3個時期：高潮帶圍墾時期(2001—2005年)，中潮帶中上部圍墾時期(2005—2008年)和中潮帶下部至低潮帶圍墾時期(2008—2013年)。3個時期平均每圍墾1 km2灘涂，潮溝長度分別減小約1.18、5.25和1.71 km，其減小速率可能與潮溝密度的分布特征有關。圖3為3個時期研究區潮溝的分布及變化情況。

① y=-1.18x+847.69，R2=0.799 6;② y=-5.25x+1 379.10，R2=0.944 3;③ y=-1.71x+737.81，R2=0.882 8;④ y=-2.60x+952.24，R2=0.949 9圖2　潮溝形態對灘涂圍墾的響應Fig.2　The respond of tidal creeks’ pattern to the reclamation

根據圖2和圖3可知，灘涂圍墾活動對潮溝長度變化產生了較大影響。其具體影響方式包括兩個方面(圖4)，一是受灘涂圍墾活動直接影響，圍墾匡圍的潮溝長度不斷增大，研究區平均每圍墾1 km2灘涂，圍墾匡圍消失的潮溝長度增加約3.32 km；二是受灘涂圍墾活動間接影響，堤壩外新生了部分潮溝，使研究區堤壩外的潮溝長度有所增大，但每年新增的潮溝長度隨著圍墾進程逐漸變小。

3.2潮溝寬度對灘涂圍墾的響應

如圖5所示，隨著灘涂圍墾面積的增加，研究區無徑流潮溝寬度呈減小趨勢，含徑流潮溝寬度先是減小，后期有所增大，總體上平均每圍墾1 km2灘涂，低潮帶每條潮溝的寬度平均減小約0.67 m。含徑流潮溝寬度在中潮帶下部至低潮帶圍墾時期有增大的趨勢，這可能與灘涂圍墾使河流作用逐漸代替潮汐作用控制了含徑流潮溝寬度有關。

根據圖5可知，灘涂圍墾活動對潮溝寬度變化產生了較大影響。根據研究區潮溝寬度與潮溝長度和集水面積之間存在的線性關系(圖6)，可將其影響過程表述如下：灘涂圍墾使潮溝長度和集水面積不斷縮小，落潮后期潮溝所接受的水量隨之減少，水動力減弱，潮溝寬度變窄。

3.3潮溝密度對灘涂圍墾的響應

如圖7所示，隨著灘涂圍墾面積的增加，研究區潮溝密度總體上呈先增加后減小趨勢，這種變化趨勢可能與研究區潮溝密度分布分帶性特征有關。圍墾初期高潮帶、中潮帶和低潮帶潮溝密度分別為0.88、1.87和1.07 km/km2，即高潮帶潮溝密度較小，高潮帶圍墾使堤壩外潮溝密度呈增大趨勢；中潮帶潮溝密度較大，中潮帶圍墾使堤壩外潮溝密度呈減小趨勢；低潮帶多為潮溝主干，低潮帶圍墾使潮溝密度Ⅰ(單位集水區內潮溝長度)有所增大，使潮溝密度Ⅱ(單位潮灘上潮溝長度)繼續減小，其中的差異可能與低潮帶圍墾較大程度地減小了潮溝集水面積有關。

如圖8所示，研究潮溝長度∑L與潮灘面積或集水面積A之間存在著冪函數關系，∑L=aAb。在潮溝長度與集水面積的關系中，參數b的變化范圍在1±0.03內，這與Marani等[37]及吳曉東和高抒[15]研究結果較一致，在一定程度上說明本文統計潮溝集水面積所用的方法具有一定的可靠性。根據圖1可知，研究區潮溝的集水區域多集中在潮灘中上部，在高潮帶和中潮帶平均每圍墾1 km2灘涂，潮溝集水面積減少約0.96 km2，即在潮灘上部約96%的灘面薄層水由潮溝系統排入外海。若用潮溝集水面積與潮灘總面積之比來評估潮溝排水作用的大小，那么在高潮帶向中低潮帶圍墾過程中潮溝排水作用大小由最初的70%逐漸降低到40%左右，這在一定程度上說明大規模灘涂圍墾使潮溝排水作用逐漸減弱，落潮后期潮灘的主要排水方式由潮溝排水逐漸向灘面排水方式轉變，從側面反映了研究區潮溝系統的退化和消亡趨勢。

圖3　研究區潮溝的分布及變化情況Fig.3　The distribution and variation of tidal creek in the study area

① y=3.32x-111.08，R2=0.998 3；② y=89.04ln(x)-253.80，R2=0.732 4圖4　潮溝長度與圍墾面積之間的關系Fig.4　The relationship between the tidal creeks’ length and reclamation area

① y=-21.588x+348 2.3，R2=0.979 1；② y=-6.225 1x+2 043.8，R2=0.979 2；③ y=5.869 2x-250.3，R2=0.861 0；④ y=0.076 3x2-28.42x+3 582.3，R2=0.972 7；⑤ y=-7.754 4x+1 558.2，R2=0.985 0；⑥ y=-8.406 9x+2 048.8，R2=0.959 8；⑦ y=-4.139 7x+1 310.4，R2=0.796 0；⑧ y=-4.286x+1 356.8，R2=0.913 2圖5　潮溝寬度對灘涂圍墾的響應Fig.5　The respond of tidal creeks’ width to the reclamation

圖6　潮溝寬度與上游潮溝長度和集水面積之間的關系Fig.6　The relationship between the tidal creeks’ width and the upstream tidal creeks’ length or watershed area

圖7　潮溝密度對灘涂圍墾的響應Fig.7　The respond of tidal creeks’ density to the reclamation

圖8　潮溝長度與潮灘面積和集水面積之間的關系Fig.8　The relationship between the tidal creeks’ length and the area of the tidal flat or the watershed area of tidal creek

3.4潮溝分維值對灘涂圍墾的響應

如圖9所示，隨著灘涂圍墾面積的增加，研究區潮溝分維值在一條直線和一條拋物線相交區域內呈波動性下降，總體上平均每圍墾1 km2灘涂，潮溝分維值減小約0.006維。潮溝分維值的快速降低反映了研究區潮溝系統的退化和消亡趨勢[40]，而潮溝分形維數在降低過程中所表現的波動性變化可能與潮溝密度分布和灘涂圍墾有關。

① y=-0.006x+1.295 5，R2=0.962 8；② y=-5×10-6x2+0.01x+1.232 9，R2=0.979 9圖9　研究區潮溝分維值對灘涂圍墾的響應Fig.9　The respond of tidal creeks’ fractal dimension value to the reclamation in the study area

在計算潮溝分維值時，用正方形盒子將研究區潮溝覆蓋起來，發現所得非空盒子數N(r)與盒子邊長r之間存在一定的標度關系N(r)～r-D，其中1.129 6

圖10　潮溝分維值與潮溝分汊性和蜿蜒性之間的關系Fig.10　The relationship between the tidal creeks’ fractal dimension value and the tidal creeks’ bifurcation or meandering

4討論

4.1灘涂圍墾對潮溝形態的影響

潮溝各形態參數之間是一個相互關聯的體系，在該體系中，任何一個參數的改變對其他參數都會產生一定影響，其影響程度一般由各形態參數之間的相關性所決定，而其影響過程往往是由潮灘水流對沉積物的搬運和堆積來完成的。潮溝長度、寬度和集水面積之間的相關性較大(R2分別為0.98和0.93)，灘涂圍墾切斷了部分潮溝長度，侵占了其集水面積，使落潮后期潮溝的納潮水量隨之減小，低潮帶水動力強度減弱，潮溝淤積作用增強，寬度變窄。潮溝系統的分維值與潮溝分汊性和潮溝密度的相關性較大(R2=0.89)，灘涂圍墾切斷了部分支潮溝，使潮溝交匯點個數隨之減少，分維值總體呈現降低趨勢，灘涂圍墾過程也改變了堤壩外潮灘水動力環境，使局部潮灘區新生了部分潮溝，新生潮溝可能是引起研究區潮溝分維值呈波動性降低的原因，而研究區潮溝密度和分汊性的潮灘分帶性差異，即中潮帶潮溝的密度和分汊率明顯大于高潮帶和低潮帶，可能是導致圍墾中期潮溝分維值波動幅度變大的原因。

4.2潮溝的起源和發育過程

江蘇中部海岸潮溝系統根據成因分為灘面水流沖刷型、陸源水流侵蝕繼承型和潮流輻聚侵蝕型潮溝[12]。若根據潮溝退化和消亡過程反推該區潮溝的起源和發育過程，陸源水流侵蝕繼承型潮溝可能是研究區最早出現的一類潮溝，它們是從入海河口(如斗龍港、王港河、川東港和東臺河等)隨著海岸淤積過程逐漸過渡延伸而來，并在后期大規模灘涂圍墾影響下又逐漸向入海河口退化，是研究區最先出現和最后消失的一類潮溝。灘面水流沖刷型和潮流輻聚侵蝕型潮溝則是在潮灘淤長過程中逐漸發育的，其早期發育過程可描述如下(圖11)：隨著海岸的淤積過程潮灘寬度不斷增大，落潮后期灘面水流對低潮灘的沖刷強度逐漸增強，當潮灘寬度增大到一定程度后，低潮灘下部開始出現較多的沖刷凹地，一些凹地在潮灘淤長過程中因匯水性較差而淤積消失，一些匯水性較好的凹地則不斷向岸延展伸長，并在潮灘微地貌的影響下出現分汊型沖溝，各沖溝支汊不斷溯源侵蝕發育成灘面水流沖刷型或潮流輻聚侵蝕型潮溝系統。影像分析發現隨著灘涂圍墾面積的不斷增加研究區堤壩外潮灘寬度快速變小，當潮灘平均寬度減小到3～5 km時，低潮帶下部開始出現較多新發育的沖溝(圖3，2005年)，這些沖溝隨著堤壩外潮灘寬度繼續減小而逐漸消失，即研究區海岸潮灘寬度淤長到3～5 km的階段可能是灘面水流沖刷型或潮流輻聚侵蝕型潮溝開始快速發育的時期。

5結論

(1)根據潮溝的分汊性和匯水特征，將研究區潮溝口、各支潮溝尾梢頂點及堤壩連線所圍區域表示為潮溝系統的集水區域，使用這種方法所計算的研究區潮溝集水面積與潮溝長度和潮溝寬度之間存在較好的相關性(R2分別為0.98和0.93)。

(2)潮溝各形態參數之間是一個相互關聯的體系。潮溝長度隨著其集水面積的縮減呈冪函數減小(R2=0.98)，潮溝寬度隨著潮溝長度和集水面積的變小呈指數方式變窄(R2分別為 0.87和0.93)；潮溝系統分維值隨著其分汊個數減少而線性降低(R2=0.89)，單支潮溝分維值隨著其蜿蜒性減小而線性降低(R2=0.86)。

(3)研究區潮溝形態在大規模灘涂圍墾影響下呈現明顯的退化和消亡趨勢。總體上每圍墾1 km2灘涂，潮溝長度、寬度和分維值分別減小約2.60 km、0.67 m和0.006維；用潮溝集水面積與潮灘面積之比來評估潮溝排水作用的大小，得出落潮后期潮灘主要排水方式由潮溝排水逐漸向灘面排水方式轉變。

參考文獻：

[1]崔承琦，印萍. 黃河三角洲潮灘發育時空譜系[J]. 青島海洋大學學報，1994，12(S3): 51-61．

Cui Chengqi，Yin Ping. The morphologic phylogenetic temporal special system of the tidal flat of the Yellow River Delta[J]. Journal of Ocean University of Qingdao，1994，12(S3): 51-61．

[2]Wolanski E. Estuarine Ecohydrology[M]. Amsterdam: Elsevier Press，2007.

[3]張忍順. 淤泥質潮灘均衡態——以江蘇輻射沙洲內緣區為例[J]. 科學通報，1995，40(4):347-350．

Zhang Renshun. Muddy tidal flat——Taking the inner of radiation shoal Jiangsu area as an example[J]. Chinese Science Bulletin，1995，40(4):347-350．

[4]謝東風，高抒，汪亞平. 砂質底質潮汐水道均衡態模擬初探[J]. 海洋學報，2006，28(6):86-93．

Xie Dongfeng，Gao Shu，Wang Yaping. A preliminary modeling study on morphodynamic equilibrium of tidal channels floored with sandy sediments[J]. Haiyang Xuebao，2006，28(6):86-93．

[5]Toffolon M，Lanzoni S. Morphological equilibrium of short channels dissecting the tidal flats of coastal lagoons[J]. J Geophys Res，2010，115: F04036．

[6]Perillo G M，Iribarne O O. New mechanisms studied for creek formation in tidal flats: from crabs to tidal channels[J]. Eos，Transactions American Geophysical Union，2003，84(1): 1-5．

[7]Vandenbruwaene W，Meire P，Temmerman S. Formation and evolution of a tidal channel network within a constructed tidal marsh[J]. Geomorphology，2012(151/152): 114-125．

[8]Novakowski K I，Torres R，Gardner L R，et al. Geomorphic analysis of tidal creek networks[J].Water Resources Research，2004，40(5): W05401．

[9]Hughes Z J. Tidal channels on tidal flats and marshes[M]//Principles of Tidal Sedimentology. Netherlands:Springer，2012: 269-300．

[10]Shi Z，Lamb H F，Collin R L. Geomorphic change of salt marsh tidal creek networks in the Dyfi Estuary，Wales[J]. Mar Geol，1995，128: 73-83．

[11]Hood W G. Tidal channel meander formation by depositional rather than erosional processes: examples from the prograding Skagit River Delta (Washington，USA)[J]. Earth Surface Processes and Landforms，2010，35(3): 319-330．

[12]邵虛生. 潮溝成因類型及其影響因素的探討[J].地理學報，1988，43(1): 35-43．

Shao Xusheng. Genetic classification of tidal creek and factors affecting its development[J]. Acta Geographica Sinica，1988，43(1):35-43．

[13]張忍順，王雪瑜. 江蘇省淤泥質海岸潮溝系統[J]. 地理學報，1991，46(2)：195-206．

Zhang Renshun，Wang Xueyu. Tidal creek system on tidal mud flat of Jiangsu Province[J]. Acta Geographica Sinica，1991，46(2):195-206．

[14]沈永明，張忍順，王艷紅. 互花米草鹽沼潮溝地貌特征[J]. 地理研究，2003，22(4):520-527．

Shen Yongming，Zhang Renshun，Wang Yanhong. The tidal creek character in salt marsh of spartina alterniflora Loisel on strong tide coast[J]. Geographical Reseach，2003，22(4):520-527．

[15]吳曉東，高抒. 長江口潮間帶九段沙淺灘潮水溝形態分析[J]. 海洋學報，2012，34(6):126-132．

Wu Xiaodong，Gao Shu. A morphological analysis of tidal creek network patterns on the Jiuduansha Shoal in the Changjiang Estuary[J]. Haiyang Xuebao，2012，34(6):126-132．

[16]Chen Y，Thompson C E L，Collins M B. Saltmarsh creek bank stability: Biostabilisation and consolidation with depth[J]. Continental Shelf Research，2012，35: 64-74．

[17]Wilson C A，Hughes Z J，Fitzgerald D M. The effects of crab bioturbation on Mid-Atlantic saltmarsh tidal creek extension: Geotechnical and geochemical changes[J]. Estuarine，Coastal and Shelf Science，2012，106: 33-44．

[18]Rizzetto F，Tosi L. Rapid response of tidal channel networks to sea-level variations (Venice Lagoon，Italy)[J]. Global and Planetary Change，2012，92: 191-197．

[19]Perillo G M E，Iribarne O O. Processes of tidal channel development in salt and freshwater marshes[J]. Earth Surface Processes and Landforms，2003，28(13): 1473-1482．

[20]Hood W G. A conceptual model of depositional，rather than erosional，tidal channel development in the rapidly prograding Skagit River Delta (Washington，USA)[J]. Earth Surface Processes and Landforms，2006，31(14): 1824-1838．

[21]D’Alpaos A，Lanzoni S，Marani M，et al. Spontaneous tidal network formation within a constructed salt marsh: observations and morphodynamic modelling[J]. Geomorphology，2007，91(3): 186-197．

[22]D’Alpaos A，Lanzoni S，Marani M，et al. Tidal network ontogeny: Channel initiation and early development[J]. Journal of Geophysical Research: Earth Surface (2003-2012)，2005，110(F2)．

[23]Hood W G. Scaling tidal channel geometry with marsh island area: A tool for habitat restoration，linked to channel formation process[J]. Water Resources Research，2007，43(3):W03409.

[24]Maanen B V，Coco G，Bryan K R. Sensitivity of modelled channel network formation to environmental conditions and initial bathymetry[J]. EGU General Assembly，2010，12：457.

[25]Tambroni N，Bolla Pittaluga M，Seminara G. Laboratory observations of the morphodynamic evolution of tidal channels and tidal inlets[J]. Journal of Geophysical Research: Earth Surface (2003-2012)，2005，110(F4):doi:10.1029/2004JF000243．

[26]Stefanon L，Carniello L，D’Alpaos A，et al. Experimental analysis of tidal network growth and development[J]. Continental Shelf Research，2010，30(8): 950-962．

[27]Vlaswinkel B M，Cantelli A. Geometric characteristics and evolution of a tidal channel network inexperimental setting[J]. Earth Surface Processes and Landforms，2011，36(6): 739-752．

[28]陳才俊. 江蘇中部海堤大規模外遷后的潮水溝發育[J]. 海洋通報，2001，20(6):71-79．

Chen Caijun. Change in tide after mudflat being enclose in the middle coast in Jiangsu pricince[J]. Marine Science Bulletin，2001，20(6):71-79．

[29]黃海軍，樊輝. 黃河三角洲潮灘潮溝近期變化遙感監測[J]. 地理學報，2004，59(5):723-730．

Huang Haijun，Fan Hui. Change detection of tidal flats and tidal creeks in the Yellow River delta using Landsat TM/ETM+images[J]. Acta Geographica Sinica ，2004，59(5):723-730．

[30]Hood W G. Indirect environmental effects of dikes on estuarine tidal channels: thinking outside of the dike for habitat restoration and monitoring[J]. Estuaries，2004，27(2): 273-282．

[31]Symonds A M，Collins M B. The establishment and degeneration of a temporary creek system in response to managed coastal realignment: The Wash，UK[J]. Earth Surface Processes and Landforms，2007，32(12): 1783-1796．

[32]吳德力，沈永明，方仁建. 江蘇中部海岸潮溝的形態變化特征[J]. 地理學報，2013，68(7): 955-965．

Wu Deli，Shen Yongming，Fang Renjian. A morphological analysis of tidal creek network patterns on the central Jiangsu coast[J]. Acta Geographica Sinica，2013，68(7):955-965．

[33]唐文周，孫國清，陳丙咸，等. 利用陸地衛星圖像和數字資料對江蘇海岸帶進行分類的研究[J]. 南京大學學報(自然科學版)，1982(4): 939-950．

Tang Wenzhou，Sun Guoqing，Chen Bingxian，et al. Classification of coast-zone in Jiangsu Province using Landsat image and digital data[J]. Academic Journal Eletronic Publishing House，1982(4): 939-950．

[34]劉燕春，張鷹. 遙感中軸線法在江蘇輻射沙洲潮溝演變監測中的應用[J]. 海洋科學，2011，35(2):72-76．

Liu Yanchun，Zhang Ying. Application of remote sensing medial axis method in investigation of dynamic changes of tidal creek in radial sandbanks offshore Jiangsu Province[J]. Marine Sciences，2011，35(2):72-76．

[35]Rinaldo A，Fagherazzi S，Lanzoni S，et al. Tidal networks 2. Watershed delineation and comparative network morphology[J]. Water Resources Research，1999，35(12): 3905-3917．

[36]Novakowski K I，Torres R，Gardner T L，et al. Geomorphic analysis of tidal creek networks[J]. Water Resources Research，2004，40(2): 1-13．

[37]Marani M，Belluco E，D’Alpaos A，et al. On the drainage density of tidal networks[J]. Water Resources Research，2003，39(2): 1040．

[38]Fagherazzi S，Bortoluzzi A，Dietrich W E，et al. Tidal network 1. Automatic network extraction and preliminary scaling feature from digital terrain maps[J]. Water Resources Research，1999，35(12): 3891-3904．

[39]崔承琦，李師湯，孫曉霞，等. 黃河三角洲海岸岸線和潮水溝體系發育及其分維值研究——黃河三角洲潮灘海岸時空譜系研究Ⅲ[J]. 海洋通報，2001，20(6):60-70．

Cui Chengqi，Li Shitang，Sun Xiaoxia，et al. Resergh on the development of the Huanghe River Delta coastline and the tidal creek system as well as fractional dimension—researgh on space-time lineage of the tidal bank coast of the Huanghe River Delta Ⅲ[J]. Marine Science Bulletin，2001，20(6):60-70．

[40]孫效功，趙海虹，崔承琦. 黃河三角洲潮灘潮溝體系的分維特征[J]. 海洋與湖沼，2001，32(1):75-80．

Sun Xiaogong，Zhao Haihong，Cui Chengqi. The fractal characteristics of tidal flat and tidal creek system in the Huanghe River Delta[J]. Oceanologia et Limnologia Sinica，2001，32(1):75-80．

[41]Angeles G R，Perillo G M E，Piccolo M C，et al. Fractal analysis of tidal channels in the Bahia Blanca Estuary (Argentina)[J]. Geomorphology，2004，57：263-274．

收稿日期：2015-03-23；

修訂日期：2015-05-26。

基金項目：氣候變化重大基礎研究項目(2013CB956503);國家自然科學基金項目(40973053)。

作者簡介：時海東(1990—), 男, 江蘇省東臺市人, 主要研究方向為海岸地貌學。E-mail:15189832427@163.com *通信作者：沈永明, 教授。E-mail:yongmsh@163.com

中圖分類號:P748

文獻標志碼：A

文章編號：0253-4193(2016)01-0106-10

Rapid response of tidal creek network patterns to the reclamation on the central Jiangsu coast

Shi Haidong1，Shen Yongming1，Kang Min1

(1.CollegeofGeographicScience，NanjingNormalUniversity，Nanjing210023，China)

Abstract:Tidal creek system，one of the most particular geomorphology units in muddy flat，is widely developed on central Jiangsu coast. This paper focuses on their evolutionary process in response to reclamation activities using TM/ETM+data of 2001-2013 and field investigation. The watershed area of tidal creeks decreased from 293.94 km2 in 2001 to 103.40 km2 in 2013 following sequential reclamation，which caused the decreasing of tidal creeks’ length with power function and tidal creeks’ width with exponential function. The average change rate of tidal creeks’ length and width was 2.60 km and 0.67 m with 1 km2 reclamation. The drainage pattern of tidal flat was changed from tidal creeks to beach，which implied that the tidal creeks systems in study area have a retrograde and demise tendency. The fractal dimension value D of each single tidal creeks and tidal creek systems were calculated from box-counting method. Results of calculation and analysis indicated that the fractal dimension value D in the region is greater than that of each tidal creek systems，followed by each single tidal creek. The fractal dimension value D in the region is positively correlated with the bifurcation ratio of tidal creek systems，and decreased from 1.26 in 2001 to 1.13 in 2013 with influences of reclamation. The fractal dimension value D of single tidal creek is positively correlated with the meandering of tidal creek. According to the retrograde and demise process of tidal creeks studied the results suggested that the tidal creek inheriting the runoff from land was the earliest tidal creeks，followed by the tidal creek scoured by the currents on tidal flat and tidal creek scoured by the converged tidal currents，whose primary form could be eroded indentation.

Key words:tidal creeks; reclamation; morphologic evolution; Jiangsu coast

時海東，沈永明，康敏. 江蘇中部海岸潮溝形態對灘涂圍墾的響應[J]. 海洋學報，2016，38(1): 106-115，doi：10.3969/j.issn.0253-4193.2016.01.010

Shi Haidong，Shen Yongming，Kang Min. Rapid response of tidal creek network patterns to the reclamation on the central Jiangsu coast[J]. Haiyang Xuebao，2016，38(1):106-115，doi：10.3969/j.issn.0253-4193.2016.01.010