孔祥淮，劉 健，張 勇，陳 斌，徐 剛

（1.國土資源部 海洋油氣資源和環(huán)境地質(zhì)重點實驗室，山東 青島266071；2.青島海洋地質(zhì)研究所，山東 青島266071）

海岸帶地區(qū)海陸相互作用的研究是全球變化研究的重要領(lǐng)域。末次冰期以來的海岸帶地區(qū)河道演變，是在全球大幅度海退與海侵及全新世以來海面波動的驅(qū)動下發(fā)育起來的，是海陸相互作用的產(chǎn)物。我們有可能根據(jù)古河道的沉積物特征和幾何形態(tài)，分析古河道的河勢特征和水文特性，恢復(fù)古地理環(huán)境。許多學(xué)者對我國海域大陸架埋藏古河道進行過研究，我國近岸海域海底殘留的古河道，按古河系可劃分為古黃河河系、古長江河系、古臺閩河系［1］和古珠江河系［2-3］。南海北部大陸架淺部地層內(nèi)，廣泛發(fā)育有第四紀的埋藏古河道［2-3］；東部陸架區(qū)低海面時期存在大量古河道系統(tǒng)［4-11］，例如，長江口外有多條大型古河道系統(tǒng)，是末次冰期長江在東海陸架平原上的主要流路，古河道分布與現(xiàn)在海底帶狀高地形有對應(yīng)關(guān)系［12］。南黃海海底存在許多埋藏古河道、古湖泊、古三角洲，北部的古河道為古黃河水系，南部的為古長江水系，二者有可能一度匯合于南黃海中部［13-14］。埋藏古河道的分布、水系屬性、形成時代以及古河道的斷面類型、充填物特征、沉積相序等包含了大量的沉積環(huán)境特征和演化過程的信息［2-3，12-15］；古河道的存在，對海上工程施工、鉆井船就位、海上平臺安裝及海上管線的鋪設(shè)等造成許多影響［16］。因此，識別古河谷，建立古河谷的沉積模式，解釋其充填過程，對于了解中國陸架海沉積作用的發(fā)育歷史和古地理環(huán)境變遷具有重要意義，可為研究全球變化提供重要信息。同時，海底埋藏古河道的研究可以為海洋工程的建設(shè)提供基礎(chǔ)地質(zhì)信息，因此具有重要的現(xiàn)實意義。在本研究中，我們將淺地層地震剖面上出露的古河道形態(tài)特征進行定量分析，類比現(xiàn)代河流環(huán)境，采用不同河道形態(tài)的水文參數(shù)變化的經(jīng)驗方程，估算古河道的坡度、河曲波長、曲率和流量，從而追溯古河道形成時期的古環(huán)境。

1 地質(zhì)環(huán)境

研究區(qū)位于南黃海西部濱淺海區(qū)，主要入海河流有淮河、灌河、射陽河、沂河以及其它較小河流，主要海灣為海州灣。海岸線在蘇北廢黃河口處向海凸出，并不斷向陸蝕退。南黃海西部近岸海域水深較淺，向東北方向逐漸加深，在研究區(qū)東部水深達45～55m。海州灣中部出現(xiàn)由二級階地構(gòu)成的剝蝕平原，一級為水深20～30m，階地面坡度約為21″，陡坎坡度約為52″；另一級水深35～45m，階地面坡度約為26″，陡坎坡度約為1°14″。老黃河水下三角洲水深10～20m，地形平坦，但有沖刷淺洼地和沙脊地形。研究區(qū)以南緊鄰蘇北大型輻射狀沙脊群，由70多條沙脊和潮流通道組成，水深0～30m，為強潮流作用地帶。

研究區(qū)有多種沿岸流（由于河流稀釋作用而鹽度較低），包括山東南部沿岸流、江蘇北部沿岸流和黃海沿岸流。研究區(qū)東部鄰近黃海暖流（黑潮的一個分支）（圖1）［17］。圖1中為2006年淺剖測線分布和4個鉆孔位置。黃海暖流具有較高的水溫，北上途中與黃海沿岸流和朝鮮半島沿岸流相互作用，在暖流西側(cè)常形成氣旋型渦旋（冷水團或冷渦），其環(huán)流體系龐大，特別是南黃海中部的冷水團環(huán)流體系，為一較弱的水動力環(huán)境，是中部泥質(zhì)沉積區(qū)形成的主導(dǎo)因素［18］。黃河、長江和沿岸短源河流帶來的泥沙是研究區(qū)沉積物輸入的主要來源。

圖1 黃海及其鄰近海域水深和冬季洋流模式Fig.1 Bathymetry and circulation patterns in winter in the Yellow Sea

2 方 法

2.1 數(shù) 據(jù)

本研究主要是以2006年由“蘇太漁01155”船獲得的4 100km 淺地層剖面的解譯成果為基礎(chǔ)的。所采用的淺剖采集系統(tǒng)適合調(diào)查區(qū)淺海環(huán)境復(fù)雜的第四紀地層結(jié)構(gòu)的研究。導(dǎo)航系統(tǒng)采用差分GPS定位系統(tǒng)（精度為15m）。地震反射數(shù)據(jù)采用的激發(fā)能量為300J，頻帶為500～1 500 Hz。脈沖信號的接收采用20單元的水聽器，震源和接收系統(tǒng)拖于船后，航行速度小于9km／h。淺剖數(shù)據(jù)垂直分辨率為10cm，水平分辨率為2m，穿透深度為海底以下70m。時間-深度是根據(jù)聲波穿過水體和沉積物分別為1 500和1 550m／s進行轉(zhuǎn)換的。水深數(shù)據(jù)也進行了潮汐改正。調(diào)查區(qū)淺地層剖面測量主測線間距約為8km，聯(lián)絡(luò)測線間距為20km（圖1）。

為了校準地震剖面，在2007年和2008年采用重力取心技術(shù)由K407（勘查船）采集了4個取心鉆孔（圖1）。這些鉆孔資料可為淺地層剖面的解譯提供地層環(huán)境信息。

2.2 淺地層剖面解譯方法

觀察地震地層的反射終止形態(tài)（削截、上超、下超等）和反射結(jié)構(gòu)（即地震相），然后辨識地震相單元及其邊界。描述深切谷的幾何特征，包括地震相在縱向和橫向上的變化、相互連接、趨勢和尺度等。

在打印的紙介質(zhì)地震剖面上描繪出埋藏古河道，截切下伏地層的傾斜面識別為古河道的側(cè)緣，地層的疊置及其深度變化揭示其地質(zhì)歷史，下部的雜亂反射代表古河道殘留沉積，上部透明的近似平行的反射層組代表細粒沉積物充填，地層之間的交切關(guān)系揭示了地質(zhì)事件發(fā)生的相對時間，所繪出的圖像的幾何特征為揭示河道的形成過程提供豐富信息。

3 分 析

3.1 古河系的重建

研究區(qū)淺剖主測線間距大多數(shù)為8km，聯(lián)絡(luò)測線間距為20km。我們推測性地把出露的古河道斷面連接起來，以期重建古河道體系。首先，我們識別出古河道斷面的2 個側(cè)緣，然后，記錄其在剖面上的位置（mark線標號）。其次，我們繪制測線航跡圖，根據(jù)記錄的mark線標號在航跡圖上用粗線條標示出古河道斷面位置，粗線條的長度就代表河道斷面的寬度。最后，把這些代表古河道的粗線條推測性地連接起來勾勒出古河道體系（圖2）。然后，我們根據(jù)從紙介質(zhì)剖面上解譯得到的河道斷面統(tǒng)計數(shù)據(jù)，采用地質(zhì)作圖軟件Suefer 8.0在計算機上自動插值成圖。對比這2種用不同方法繪制的圖件，然后進行修改以期獲得盡可能接近真實的古河系。圖3就是對比分析后繪制的古河道體系。

圖2 推測連接的埋藏古河系Fig.2 The deduced buried-channel system

3.2 古河道形態(tài)特征

圖3 解譯的埋藏古河道軸線（紅線代表圖4中剖面的位置）Fig.3 Interpreted axes of shallowly buried incised channels（Red-highlighted lines indicate the location of seismic profiles shown in Fig.4）

研究區(qū)東西方向和南北方向測線淺地層剖面的解譯揭示了復(fù)雜的樹枝狀古河道網(wǎng)，它們埋藏于南黃海西部（江蘇省北部濱淺海區(qū)）海底以下數(shù)米到數(shù)十米。在30～50m 水深范圍內(nèi)，埋深一般＜5m；而在10～30m 水深范圍內(nèi)，埋深通常達20m，因其上被古三角洲覆蓋。在水深10～50m 的近50 000km2范圍內(nèi)，我們繪出了兩大河道網(wǎng)（圖2和圖3）。北部河道網(wǎng)的主干河道呈西南-東北走向，蜿蜒而行與現(xiàn)代海底等深線斜交；南部河道網(wǎng)的主干河道大致呈東西走向。2個河系均表現(xiàn)為北側(cè)支流較密（圖3）。這些河道在152 km 距離內(nèi)落差為31m，以0.012°的坡度向海傾斜。北部主河道橫截面向海逐漸變寬，達19km，南部主河道向海變寬為10km。

從解譯的古河系可以明顯看出，2個河系的支流分布是不對稱的，表現(xiàn)為梳狀格架（圖3）。大多數(shù)支流從北側(cè)交匯于主河道，且交結(jié)角接近90°。支流之間交匯角為銳角，一般＜60°。在兩個河系之間存在一個分水嶺，2個河系的分布范圍相差很大，北部河系流域面積較大，南部河系流域面積較小。

河道寬度從數(shù)百米到數(shù)千米，深度從幾米到十幾米不等。根據(jù)解譯的古河道形狀、大小和充填沉積物反射特征，古河道截面可以分成不同類型。一般說來，箱形截面主要存在于主干河道（圖4a和圖4b），透鏡狀截面主要見于二級河道（圖4c），而V 形截面通常在三級河道或更小的河道中出現(xiàn)（圖4d）。

河道充填單元的內(nèi)部反射結(jié)構(gòu)也不相同。在主干河道內(nèi)，強振幅雜亂反射常常出現(xiàn)在河道底部，其上為一套復(fù)雜的反射結(jié)構(gòu)，表現(xiàn)為加積、側(cè)積層組、侵蝕削截或不同傾斜角度的多期次層狀反射（圖4a和圖4b）。在某些情況下，常常會出現(xiàn)下部較老的層組被新的反射界面所削截。再往上部，可以清楚地看到一薄層近似平行和水平的沉積覆蓋層，其上下界面均為明顯的不整合接觸面。底界面削截下伏地層，而其頂界面又被一個S形進積反射楔狀體所覆蓋。另外，在河谷充填物的某些層位還可能有大片淺層氣屏蔽了下伏的地層反射。在中等河道或較小的支流河道里，往往出現(xiàn)反射層組由一側(cè)向另一側(cè)充填或者反射層組中部下凹現(xiàn)象（圖4c和圖4d）。

河道側(cè)緣常常以低角度（＜15°）傾斜面截切下伏近似水平的地層為特征，河道2個側(cè)緣有的對稱，有的不對稱，如河道的一邊側(cè)緣短而陡，另一側(cè)長而緩（圖4c）。主干河道河床底界面高低起伏說明下切深度不同（一系列較大規(guī)模的起伏可能意味著其為辮狀河）（圖4a和圖4b）。

圖4 淺地層剖面中的埋藏古河道斷面Fig.4 Cross-sections of the buried channels in shallow profiles

3.3 古河道形態(tài)定量分析

要想尋找古河道形態(tài)與古水文條件之間的關(guān)系，就要進行古河道形態(tài)測量分析。我們按照聲速1 550 m／s將雙程反射時間轉(zhuǎn)換成沉積地層的厚度，這個聲速是聲波穿過海底黏土質(zhì)砂等沉積物時的典型聲速［19］。根據(jù)對保存下來的古河道幾何形態(tài)的測量或推測的參數(shù)，我們可以利用現(xiàn)代河流和河口的水力學(xué)經(jīng)驗方程估算古流量。

根據(jù)淺地層剖面上觀察到的或者解譯的古河道形態(tài)特征，我們可以得出古水力學(xué)參數(shù)。在垂直于每條河道的方向上測量橫截面參數(shù)，測量參數(shù)包括寬度、平均深度和最大深度、河曲的直線長度和河曲路徑長度（圖5和表1）。相關(guān)計算公式如下：

式中，di為古河道實測深度，i=1～n。

式中，A 為古河道橫截面積，w 為河寬。

埋藏古河道截切下伏的平行地層，在河道底部是以雜亂反射為特征的充填單元，這可能代表了陸地河流過程的礫質(zhì)殘留沉積。另外，樹枝狀的河道網(wǎng)和測得的河道交結(jié)角度與陸源河流一致。而其上部主要以大型槽狀交錯層理和斜層理為主，見平行層理、波狀層理和極少量的水平層理，可能代表潮汐作用過程。研究區(qū)曾經(jīng)歷了從河流到河口、再到被海水淹沒的環(huán)境變遷。因此，我們認為，古河道充填是河流的侵蝕和沉積作用以及隨后受到潮汐改造作用的混合的結(jié)果。在我們的分析中，我們把河流和河口環(huán)境作為對蘇北陸架古河道水力學(xué)幾何特征的可能的解釋。因此，我們將其與現(xiàn)代邊緣海系統(tǒng)［21］和中等尺度的中緯度河流［22-23］的數(shù)學(xué)模型進行了對比。

圖5 典型河道斷面形態(tài)幾何參數(shù)的測量Fig.5 Measurements of the geometric parameters of typical channel sections

表1 南黃海西部濱淺海區(qū)古河道水力學(xué)參數(shù)Table 1 Hydraulic parameters of the buried paleo-channels in the western area of the South Yellow Sea

古環(huán)境水力學(xué)限制條件可作如下界定：河流古流量是給定河道斷面的上游集水區(qū)匯集到該斷面的平灘流量［24］，而在河流入?？诨蚝涌跒?，河流流量就是進潮量［25］。

潮汐河口流量可利用下面的冪函數(shù)進行估算：

式中，Q 為排水量，A 為橫截面積。我們在古流量估算中采用平均值α=0.96［21，26］。

古河系流量可利用Dury［23］的經(jīng)驗方程進行估算：

根據(jù)上面估算的流量，可以從下面的水流連續(xù)性方程求得河流的平均流速：

式中，v為平均流速。

然后，根據(jù)顆粒運動狀態(tài)、粒徑大小和水流流速三者之間的關(guān)系［27-28］估算古河道攜帶沉積物的能力。

表2 埋藏古河道古水流估算Table 2 Estimation of paleo-flow in the buried paleo-channels in the western area of the South Yellow Sea

3.4 河道古水流估算

埋藏古河道分析計算的水力學(xué)參數(shù)見表1。所繪制的埋藏古河道向海變寬變深（圖2和圖3）。較低的寬深比一般與V 形的較小支流相對應(yīng)，而較高的寬深比一般與較大的箱形的主干河道相對應(yīng)。

根據(jù)古河道橫截面積和前述的邊緣海系統(tǒng)水力學(xué)方程進行估算，古潮汐通道潛在的最大進潮量可能為67 000～630 000m3／s，相應(yīng)的平均流速為150～170cm／s（表2）；古河道潛在的最大流量可能為94 000～970 000m3／s，相應(yīng)的平均流速為200～260cm／s（表2）。實際上，以上計算出的水力參數(shù)應(yīng)大大地打個折扣，因為古河谷斷面中包含了一系列復(fù)雜的深切水道和充填單元，這些深切谷或許形成于不同時期，或許是某個深切谷隨著時間的推移發(fā)生了遷移，形成新河道，舊的河道被淤積充填。在我們繪出的主干河道中，根據(jù)內(nèi)部聲學(xué)反射特征已經(jīng)辨識出至少十幾處古水道遺跡（圖4a和b）。因此，我們采用測量參數(shù)（主要指河道寬度）的10%來估算古流量。那么，最大古進潮量約為3 820～14 700m3／s，相應(yīng)的流速約為140～150 cm／s；古河流流量約為4 600～19 000m3／s，相應(yīng)的平均流速約為170～190cm／s。

在邊緣海系統(tǒng)里，以上估算的水流速度足可以起動并搬運直徑為2～20mm 的顆粒（即從極粗砂到礫）作為床砂載荷；在河流系統(tǒng)里，估算的流速可以起動直徑達40mm 的顆粒（礫石）沿河床滾動。

3.5 河流分類

解譯的古河系河流類型的劃分是推斷古環(huán)境的一個重要參數(shù)，因此我們對照Rosgen［29］分類方案進行河流分類。Rosgen是根據(jù)以下參數(shù)進行分類的，如寬深比、曲折度和河床比降［29］。在我們的研究中，對解譯的埋藏古河道體系進行測量或計算得出的水力學(xué)變量可以直接與Rosgen分類方案中的形態(tài)特征進行對照。研究區(qū)古河道以低曲折度和比降、高寬深比為特征，與D 型（辮狀河）河流特征一致，因此應(yīng)當屬于“辮狀河”。

4 討 論

埋藏古河道一般由3個地震相單元組成（圖4a和b）。最低的地震相單元（介于界面T2-T3）是首先由河流下切作用形成的，因為它下切到更老的更新世地層里，是一組位于古河道底部的強振幅的沉積物充填反射和雜亂反射，很可能代表河流作用殘留的礫質(zhì)沉積物。另外，樹枝狀的河系結(jié)構(gòu)特征也意味著存在地表徑流的陸地過程。根據(jù)在50m 水深處采集的鉆孔SYS-0803巖心數(shù)據(jù)，這個沉積單元顯示大致向上變細的粒序，其中出現(xiàn)在底部的礫石碎屑為園狀，底部貝殼礫石層被解釋為侵蝕滯留沉積；而且在這一單元中還可觀察到植物根和鈣質(zhì)結(jié)核［30］。由該鉆孔巖心樣品的14C測年結(jié)果可知，古河道形成于距今約4.4萬a（圖4a），當時的海面位于現(xiàn)代海平面以下大約60～70m［26］，整個研究區(qū)出露地表。從充填反射結(jié)構(gòu)可以明顯看出（圖4a），在距今4.4萬a至大約1.3萬a期間，箱形古河槽中有許多水道遺跡，表現(xiàn)為或并列或疊置的反射結(jié)構(gòu)。這種反射結(jié)構(gòu)表明，基準面較低時，河流首先以下切侵蝕為主，然后，側(cè)向侵蝕發(fā)育，使河谷變寬，成為曲流河或辮狀河。這期間海面變化幅度約60～70 m（下降、再上升），但始終低于現(xiàn)代海平面50 m 以下［26，31-32］。以上估算的末次冰期河流流量變化范圍在4 600～19 000m3／s，這與長江當時的流量12 000～16 000m3／s［33］相當，但是低于現(xiàn)代長江流量（大通水文站統(tǒng)計數(shù)據(jù)表明，在1950～2000年期間平均流量為28 588m3／s［34］），遠遠高于現(xiàn)代黃河流量［34］。因此，我們推測，在末次冰期，研究區(qū)古環(huán)境與當時的長江流域相似。

中部地震相單元以界面T1 和T2 為其上下邊界（圖4a和b）。根據(jù)鉆孔SYS-0803 樣品14C 測年結(jié)果［30］，T1界面的年齡為距今10 975a，T2界面為12 600a，那時海面在現(xiàn)代海平面以下約40m［31-32］。這個地震相單元表現(xiàn)為近似平行和水平層組的加積疊置的內(nèi)部反射結(jié)構(gòu)（圖4a和b）。其下部邊界切割下伏地層，與下伏地層不整合接觸，上部邊界與上覆地層為角度不整合接觸。在距今1.16～1.13萬a期間，由于發(fā)生融水脈沖事件，海平面突然上升約15m［35］，研究區(qū)古河系逐漸被淹沒、填埋，在隨后的海進過程中，這些河道還受到海洋潮汐的改造作用。因此，該地區(qū)經(jīng)歷了由河流到河口，再到淺海的環(huán)境變遷。與SYS-0803巖心對照，這一單元由灰色細砂和深灰色粘土質(zhì)粉砂組成，整體上向上變細，生物擾動弱-中等。常見潮汐層理，含少量貝殼碎片。其頂界面T1對應(yīng)于中全新世最大海泛面。因此，這一沉積單元反映了從濱海向陸架過渡、水深逐漸加大的沉積環(huán)境變遷，形成時間約為距今1.3萬a至全新世的中期，對應(yīng)于冰后期海平面的上升［30］。

上部地震相單元（T0-T1）代表距今10 975a以來高海面時期的沉積層反射（圖4）。這個沉積單元是下超的S形進積反射結(jié)構(gòu)，推測其為強侵蝕作用使附近海區(qū)細粒沉積物再懸浮并部分地輸運到殘留河谷區(qū)沉積形成的。

5 結(jié) 論

根據(jù)淺地層剖面解譯，在末次冰期時，南黃海西部陸架區(qū)曾發(fā)育了2大河系網(wǎng)，均表現(xiàn)為類似梳子形狀的河系網(wǎng)絡(luò)結(jié)構(gòu)，北側(cè)支流多于南側(cè)支流。北部河系流域面積較大，主干河道由西南向東北蜿蜒而行，其西端位于1128-1855年老黃河（現(xiàn)已廢棄）入?？诟浇?。南部河系較小，主干河道大致東西走向，其西端位于現(xiàn)在的射陽河口以南。

在末次冰期，研究區(qū)2大古河系流域內(nèi)氣候濕潤，河流流速大，攜帶泥沙能力強。河流流量低于現(xiàn)代長江下游大通站實測流量，遠大于現(xiàn)代黃河流量，但與當時的長江流量相當。

在全新世海進時期，古河系被淹沒、填埋，并遭受海洋潮汐改造作用。隨著海面進一步升高，古河道被一薄層海洋砂質(zhì)沉積物覆蓋，但在地形上仍低于河谷兩岸。

在高海面時期，近岸的古河道被后來形成的水下三角洲所覆蓋，埋藏深度達20m。在離岸較遠的海域，古河道鄰近區(qū)域海底沉積物再懸浮、搬運，充填了古河道凹地，使其在地形上已完全看不出河道痕跡，古河道埋藏深度一般小于5m。

總之，自末次冰期以來，研究區(qū)經(jīng)歷了從裸露地表、發(fā)育河流、潮汐河口到淺海的環(huán)境變遷。

：

［1］ GENG X S.Paleo-river systems on the eastern shelf of China［J］.Marine Sciences，1981，（2）：22-27.耿秀山.中國東部陸架的海底古河系［J］.海洋科學(xué)，1981，（2）：22-27.

［2］ BAO C W.Buried ancient channels and deltas in the Zhujiang River mouth shelf area［J］.Marine Geology &Quaternary Geology，1995，15（2）：25-34.鮑才旺.珠江口陸架區(qū)埋藏古河道與古三角洲［J］.海洋地質(zhì)與第四紀地質(zhì)，1995，15（2）：15-34.

［3］ KOU Y Q，Du D L.Sedimentary features of shallow ancient river channels on the northern shelf of the South China Sea.［J］.Acta Geologica Sinica，1994，68（3）：268-277.寇養(yǎng)琦，杜德莉.南海北部陸架第四紀古河道的沉積特征［J］.地質(zhì)學(xué)報，1994，68（3）：268-277.

［4］ BUTENKO J，YE Y C，MILLIMAN J D.Morphology，sediments and Late Quaternary history of the East China Sea.In：Qingming J，Milliman J D，eds.Sedimentation on the Continental Shelf with Special Reference to the East China Sea［C］.Beijing：China Ocean Press，1983，653-677.

［5］ MILLIMAN J D，QIN Y S，PARK Y A.Sediments and sedimentary processes in the Yellow and East China Seas.In：Taira A，Masuda F，eds.Sedimentary Facies in the Active Plate Margin［C］.Tokyo：Terra Scientific Publishing Company，1989，233-249.

［6］ LIU Z X.Yangtze Shoal a modern tidal sand sheet in the northwestern part of the East China Sea［J］.Marine Geology，1997，137：321-330.

［7］ SAITO Y，KATAYAMA H，IKEHARA K，et al.Transgressive and highstand systems tracts and post-glacial transgression，the East China Sea［J］.Sedimentary Geology，1998，122（1-4）：217-232.

［8］ LIU Z X，BERN S，the L'ATALANTE Scientific Party.Paleochannels and paleodeltas in the continental shelf of the East China Sea［J］.Marine Geology &Quaternary Geology，2000，20（1）：9-14.劉振夏，BernéS，L'ATALANTE 科學(xué)考察組.東海陸架的古河道和古三角洲.海洋地質(zhì)與第四紀地質(zhì)，2000，20（1）：9-14.

［9］ LIU Z X，BernéS，Saito Y，et al.Quaternary seismic stratigraphy and paleoenvironments on the continental shelf of the East China Sea［J］.Journal of Asian Earth Sciences，2000，18（4）：441-452.

［10］ YOO D G，LEE C W，KIM S P，et al.Late Quaternary transgressive and highstand systems tracts in the northern East China Sea midshelf［J］.Marine Geology，2002，187：313-328.

［11］ BERN S，VAGNER P，GUICHARD F，et al.Pleistocene forced regressions and tidal sand ridges in the East China Sea［J］.Marine Geology，2002，188：293-315.

［12］ LI G X，LIU Y，YANG Z G，et al.The Changjiang paleochannels on the shelf plain of the East China Sea during the last glaciation［J］.SCIENCE IN CHINA Ser.D Earth Sciences，2005，35（3）：284-289.李廣雪，劉勇，楊子賡，等.末次冰期東海陸架平原上的長江古河道［J］.中國科學(xué)：D輯（地球科學(xué)），2005，35（3）：284-289.

［13］ QIN Y S，LI F，TANG B J，et al.Buried paleo-channel systems in the west of the South Yellow Sea［J］.Chinese Science Bulletin，1986，（24）：1887-1890.秦蘊珊，李凡，唐寶玨，等.南黃海西部埋藏古河系［J］.科學(xué)通報，1985，（24）：1887-1890.

［14］ LI F，YU J J，JIANG X H，et al.Study on Buried paleo-channel system in the South Yellow Sea［J］.Oceanologia et Limnologia sinica，1991，22（6）：501-508.李凡，于建軍，姜秀珩，等.南黃海埋藏古河系研究［J］.海洋與湖沼，1991，22（6）：501-508.

［15］ LI C X，ZHANG G J.Progress on high resolution sequence stratigraphy in incised paleovalley［J］.Advance in Earth Sciences，1996，11（2）：216-220.李從先，張桂甲.下切古河谷高分辨率層序地層學(xué)研究的進展［J］.地球科學(xué)進展，1996，11（2）：216-220.

［16］ WANG M T，ZHUANG Z Y，GE S L，et al.Sediment Characteristics of the Shallow Buried paleo-Channels in the North-central Liaodong Bay and Their Negative Effects on Nautical Engineering［J］.Journal of Oceanography of Huanghai &Bohai Seas，2000，18（2）：18-24.王明田，莊振業(yè)，葛淑蘭，等.遼東灣中北部淺層埋藏古河道沉積特征及對海上工程的影響［J］.黃渤海海洋，2000，18（2）：18-24.

［17］ LIU J，SAITO Y，WANG H，et al.Stratigraphic development during the Late Pleistocene and Holocene offshore of the Yellow River delta，Bohai Sea［J］.Journal of Asian Earth Sciences，2009，36（4-5）：318-331.

［18］ SHEN S X，LI A C，YUAN W.Low Energy Environment of the Central South Yellow Sea［J］.Oceanologia et Limnologia Sinica，1996，27（5）：518-523.申順喜，李安春，袁巍.南黃海中部的低能沉積環(huán)境［J］.海洋與湖沼，1996，27（5）：518-523.

［19］ ORSI T H，DUNN D A.Correlations between sound velocity and related properties of glacio-marine sediments：Barents Sea［J］.Geo-Marine Letters，1991，11（2）：79-83.

［20］ LEOPOLD L B，WOLMAN M G，MILLER J P.Fluvial Processes in Geomorphology［M］.San Francisco：Dover Pablications，1964：522.

［21］ FRIEDRICHS C T.Stability，shear stress and equilibrium cross-sectional geometry of sheltered tidal channels［J］.Coast.Res.，1995，11（4）：1062-1074.

［22］DURY G H.Attainable standards of accuracy in the retrodiction of palaeodischarge from channel dimensions［J］.Earth Surf.Process.Landf.，1985，10（3）：205-213.

［23］ DURY G H.Discharge prediction，present and former，from channel dimensions［J］.J.Hydrol.，1976，30（3）：219-245.

［24］ KNIGHTON D.Fluvial Forms and Processes：A New Perspective.Arnold Publishers，London.1998：389.

［25］ JARRETT J T.Tidal Prism-Inlet Area Relationships［R］.U.S.Army Coastal Engineering Research Center，Vicksburg，MS.1976，vol.3.

［26］ NORDFJORD S，GOFF J A，AUSTIN J A J，et al.Seismic Geomorphology of Buried Channel Systems on the New Jersey Outer Shelf：Assessing Past Environmental Conditions［J］.Marine Geology，2005，214（4）：339-364.

［27］ HJULSTROME F.Transportation of detritus by moving water.In：P D Trask ed：Recent Marine Sediments：A Symposium［C］，Spec.Pub.Soc.econ.Paleont.Miner，4，Tulsa，1939：5-31.

［28］ SUNDBORG A.The River Klaralven：A study of fluvial processes［J］.Geogr.Annr，1956，38（3）：238-316.

［29］ ROSGEN D L.A classification of natural rivers.Catena，1994，22：169-199.

［30］ LIU J，SAITO Y，KONG X H，et al.Delta development and channel incision during marine isotope stages 3and 2in the western South Yellow Sea.Marine Geology，2010，278：54-76.

［31］ LIU J P，MILLIMAN J D，GAO S，et al.Holocene development of the Yellow River's subaqueous delta，North Yellow Sea［J］.Mar.Geol.2004，209，45-67.

［32］ SAITO Y.Sea levels of the Last Glacial in the East China Sea Continental Shelf（In Japanese with English abstract）［J］.Quat.Res.1998，37（3）：235-242.

［33］ CAO G J，WANG J，ZHANG X Q，et al.Characteristics and runoff volume of the Yangtze River's paleovalley at Nanjing Reach in the Last Glacial Maximum［J］.Acta Geographica Sinica.2009，64（3）：331-338.曹光杰，王建，張學(xué)勤，等.末次盛冰期長江南京段河槽特征及古流量［J］.地理學(xué)報，2009，64（3）：331-338.

［34］ YOU X T，ZHU H.The discharge variation in the two main rivers of China during the recent half century and the relationship with the south-to-north water transfer project［J］.Advance in Earth Sciences，2002，17（6）：811-817.游性恬，朱禾.兩大江河流量的半世紀變化與“南水北調(diào)”［J］.地球科學(xué)進展，2002，17（6）：811-817.

［35］ LIU J，SAITO Y，WANG H，et al.Sedimentary evolution of the Holocene subaqueous clinoform off the Shandong Peninsula in the Yellow Sea［J］.Marine Geology，2007，236（3-4）：165-187.