田 元 范德江① 張喜林 陳 彬 王 亮 逄 悅

(1. 中國海洋大學(xué)海洋地球科學(xué)學(xué)院 青島 266100; 2. 國土資源部青島海洋地質(zhì)研究所 青島 266071;3. 國家海洋局第三海洋研究所 廈門 361005)

粒度是沉積物的基本屬性, 它是沉積物分類和命名的依據(jù), 如海洋沉積物Shepard分類方法和Folk的分類方法, 它們都是依據(jù)沉積物粒級(jí)組成進(jìn)行了分類(Shepard, 1954; Folk et al, 1970)。粒度也是進(jìn)行沉積環(huán)境恢復(fù)和重建的重要指標(biāo), 其中粒級(jí)參數(shù)、粒級(jí)分布和統(tǒng)計(jì)學(xué)圖表在沉積相和古環(huán)境研究中被廣泛應(yīng)用(Krumbein, 1934; Folk et al, 1957; Passega,1957; Visher, 1969; Clark, 1976; Syvitski, 1991; Prins et al, 2000; Davis et al, 2002)。近年來出現(xiàn)了敏感粒級(jí)的概念, 它特指在巖芯中變化強(qiáng)烈的粒級(jí)組分, 并借助粒度-標(biāo)準(zhǔn)偏差法、主成分因子分析法、擬合函數(shù)法和端元粒度模型提取敏感性粒級(jí)(孫東懷等, 2001;陳國成等, 2007; 張存勇等, 2009; Weltje, 1997;Boulay et al, 2003)。學(xué)者們?cè)噲D利用敏感性粒級(jí)組成作為特定地質(zhì)作用的指標(biāo)(肖尚斌等, 2005a; Liu et al,2010b; Fan et al, 2011; Huang et al, 2011; Hu et al,2012), 如肖尚斌等(2005a)通過敏感粒度組分并反演了 8ka來的東亞季風(fēng)變化, 向榮等(2006)依據(jù)巖芯中敏感性粒級(jí)組分重建了末次冰期以來東亞冬季風(fēng)的變化; 然而, 李云海等(2010)對(duì)閩浙沿岸泥質(zhì)區(qū)南部沉積中心的沉積物冬、夏季表層樣的分析, 認(rèn)為冬季和夏季敏感粒級(jí)相同, 且該區(qū)域的表層樣敏感粒級(jí)和柱狀樣敏感粒級(jí)相似。迄今為止, 不同學(xué)者對(duì)敏感性粒級(jí)的構(gòu)成以及替代性作用的認(rèn)識(shí)不一致, 缺乏對(duì)敏感性粒級(jí)組分成因機(jī)制的了解。

東海內(nèi)陸架泥質(zhì)區(qū)高海平面以來沉積環(huán)境穩(wěn)定,是在東亞季風(fēng)、長江輸入、海洋動(dòng)力環(huán)境等因素共同影響下沉積作用的結(jié)果; 該區(qū)沉積體厚度大、連續(xù)性好, 成為研究古環(huán)境信息的重要載體。前人利用沉積物巖芯沉積學(xué)記錄開展了較多的冰后期以來古環(huán)境、古氣候方面的研究, 取得了重要的成果(肖尚斌等,2005b; 向榮等, 2006; 劉升發(fā)等, 2010; Wang et al,2005; Liu et al, 2010a; Liu et al, 2010b), 其中個(gè)別學(xué)者涉及到敏感性粒級(jí)的古氣候恢復(fù)研究(向榮等,2006)。前人的研究多集中在單個(gè)巖芯的敏感粒級(jí)提取與解析, 對(duì)內(nèi)陸架不同區(qū)域的對(duì)比研究較少見, 對(duì)東海內(nèi)陸架泥質(zhì)區(qū)從北到南不同沉積位置的敏感粒級(jí)組成的差異不清楚, 對(duì)敏感性粒級(jí)的形成機(jī)制未進(jìn)行深入分析。為此, 本文通過對(duì)東海內(nèi)陸架泥質(zhì)區(qū)3根巖芯(圖 1)敏感性粒級(jí)的對(duì)比研究, 闡明該區(qū)的敏感性粒級(jí)構(gòu)成, 探討敏感性粒級(jí)的形成機(jī)制, 為更好地利用敏感性粒級(jí)進(jìn)行古環(huán)境和古氣候研究提供依據(jù)。

1 研究區(qū)概況

東海西靠中國大陸, 東臨琉球群島, 北接黃海,南以臺(tái)灣海峽與南海相連, 面積約 70×104km2, 平均水深為 130m。東海大陸架是世界上最為寬廣的陸架之一, 平均坡度僅為4′17″, 平均水深78m。50—60m等深線將其分為內(nèi)陸架和外陸架兩部分, 內(nèi)陸架海底地形底質(zhì)相對(duì)復(fù)雜, 外陸架相對(duì)平緩。

沿著東海內(nèi)陸架發(fā)育著名的泥質(zhì)沉積帶, 它發(fā)育在水深50m左右的近海海域, 緊靠浙江、福建海岸(秦蘊(yùn)珊等, 1982; 郭志剛等, 2003; Saito et al,1995)(圖1)。泥質(zhì)區(qū)向外則是相對(duì)粗顆粒的被改造殘留沉積物(Liu et al, 2007)。長江輸入東海的巨量沉積物以及該區(qū)的海洋動(dòng)力環(huán)境共同塑造了東海內(nèi)陸架現(xiàn)代沉積作用。由河流帶來的巨量泥沙輸送入海并沉積于河口以及內(nèi)陸架等海域, 形成厚度不一的陸源碎屑沉積層(秦蘊(yùn)珊, 1987; 李家彪, 2008)。內(nèi)陸架區(qū)域沉積動(dòng)力環(huán)境主要受控閩浙沿岸流和臺(tái)灣暖流(郭志剛等, 2000)。其中閩浙沿岸流隨季節(jié)變化, 夏季因東南季風(fēng)盛行流向北,冬季由于偏北季風(fēng)盛行而向南運(yùn)移(秦蘊(yùn)珊, 1987); 而臺(tái)灣暖流在東海中內(nèi)陸架終年存在,沿50—100m等深線向北流動(dòng)(Su et al, 1987)。冬季閩浙沿岸流將再懸浮的長江口門及附近的沉積物向南輸運(yùn), 又由于臺(tái)灣暖流的“頂托”作用阻擋了其攜帶的物質(zhì)向東海中外陸架輸送, 使得沉積物沉積在123°E以西的閩浙沿岸, 形成東海內(nèi)陸架泥質(zhì)沉積區(qū)(楊作升等, 1992; 郭志剛等, 1999), 是東海現(xiàn)代陸架的沉積中心, 是長江入海物質(zhì)的“匯”(Yang et al,1994; Saito et al, 1995)。

2 采樣與研究方法

2.1 研究站位和樣品采集

沉積物巖芯D1、D2、D3由“東方紅2號(hào)”科學(xué)考察船分別在2006年、2009年、2011年的三個(gè)航次中利用重力取樣器取得。其中, 巖芯 D1(122.6167°E,31.0167°N)位于長江口泥質(zhì)區(qū), 水深 20.4m, 巖芯長度為 226cm。巖芯 D2(122.4529°E, 29.2088°N)、D3(121.8690°E, 27.7997°N)位于閩浙沿岸泥質(zhì)沉積區(qū),水深分別為 40m、53.8m, 巖芯長度均為 180cm。巖芯站位見圖1。所采集的沉積物巖芯于室內(nèi)進(jìn)行分樣,從中間剖分, 按照 0.25cm 的間隔分樣, 并把樣品裝入塑料樣品袋中, 室溫保存。

圖1 研究區(qū)及采樣站位分布圖。底圖據(jù)Su等(1987)和Liu等(2006)重繪Fig.1 Map of sampling sites. The base map was modified from Su(1987) and Liu et al(2006)

2.2 粒度分析方法

采用激光粒度分析儀進(jìn)行沉積物的粒度分析。實(shí)驗(yàn)步驟如下: 取0.5—1g左右的樣品放入100mL的干凈燒杯中, 加入 5mL濃度為 30%的 H2O2, 靜置 24h以去除有機(jī)質(zhì); 加入3mL分散劑(0.5mol/L六偏磷酸鈉)超聲 30min確保足夠分散并待上機(jī)測(cè)試。實(shí)驗(yàn)儀器為英國Malvern公司的Mastersizer 2000 型激光粒度分析儀, 該儀器測(cè)量范圍為 0.02—2000μm, 粒級(jí)分辨率為0.01Φ, 重復(fù)測(cè)量的相對(duì)誤差＜3%。實(shí)驗(yàn)在中國海洋大學(xué)海底科學(xué)與探測(cè)教育部重點(diǎn)實(shí)驗(yàn)室完成。粒度分級(jí)采用伍登-溫德華-Φ 值標(biāo)準(zhǔn), 粒度數(shù)據(jù)以 0.25Φ的采集間隔導(dǎo)出, 粒度參數(shù)采用 McManus矩法公式(McManus,1988)計(jì)算。

3 結(jié)果與討論

3.1 巖芯的粒度特征

巖芯的粒度具體統(tǒng)計(jì)結(jié)果見表1。三根巖芯沉積物組成均以粘土和粉砂為主, 其中砂平均含量均小于 3%, 粉砂平均含量在 67%—71%之間, 粘土平均含量在27%—31%之間。各巖芯平均粒度參數(shù)中, 中值粒徑 Md 7.1—8.9μm, 平均粒徑 Mz 6.5—10.9μm;分選δ均較差; 偏態(tài)Sk僅D2為負(fù)偏態(tài), 其余為正偏態(tài); 峰態(tài)Ku均為寬峰態(tài)。

從平均粒徑垂向變化圖(圖2)上可以看出, 3根巖芯的平均粒徑主要分布在6—8μm之間, 值得注意的是, 各巖芯均有明顯峰值出現(xiàn)。

表1 巖芯粒度參數(shù)Tab.1 Grain size distribution in the sediments from the cores

圖2 平均粒徑垂向分布圖Fig.2 The vertical profiles of the mean grain-size distribution

3.2 敏感粒級(jí)提取

采用粒徑-標(biāo)準(zhǔn)偏差法(Boulayet al, 2003)對(duì)粒度分析結(jié)果進(jìn)行敏感粒級(jí)的提取, 各個(gè)巖芯的粒級(jí)-標(biāo)準(zhǔn)偏差曲線如圖3所示。根據(jù)圖3確定各個(gè)巖芯的敏感粒級(jí), 巖芯 D1、D2、D3粒級(jí)-標(biāo)準(zhǔn)偏差曲線基本一致, 它們可識(shí)別出3個(gè)敏感粒級(jí)區(qū)間, 分別稱之為敏感粒級(jí)1、敏感粒級(jí)2和敏感粒級(jí)3(表2), 其中敏感粒級(jí)2在所有3個(gè)巖芯中標(biāo)準(zhǔn)偏差值最大, 表明其在巖芯中的波動(dòng)強(qiáng)烈, 具有最好的環(huán)境變化指示作用。從長江口向南, 敏感粒級(jí)區(qū)間逐漸向細(xì)顆粒方向偏移的趨勢(shì), 以敏感粒級(jí) 2、3的界線的偏移更加明顯, 巖芯 D1、D2、D3 的該界線分別為 272.62、229.25、114.63μm。敏感粒級(jí) 1、2的界線亦有偏移, 但是偏移幅度很小。這種偏移現(xiàn)象與長江入海沉積物向南搬運(yùn)過程中發(fā)生的沉積分異作用有關(guān)。

圖3 粒徑-標(biāo)準(zhǔn)偏差曲線圖Fig.3 The standard deviation vs. grain-size of the cores

表2 巖芯敏感粒級(jí)(單位μm)Tab.2 The sensitive grain-size groups

3.3 敏感性粒級(jí)形成機(jī)制分析

東海內(nèi)陸架沉積物以陸源碎屑占據(jù)絕對(duì)優(yōu)勢(shì),沉積物粒級(jí)構(gòu)成受到牽引流搬運(yùn)和沉積作用的總體控制。為此, 可以利用沉積物粒度分布形式揭示沉積物的搬運(yùn)特征, 進(jìn)而探討敏感粒級(jí)的形成機(jī)制。

3.3.1 搬運(yùn)特征分析 沉積物粒度的概率累積曲線能直觀地展現(xiàn)沉積物搬運(yùn)特征(Visher, 1969;Glaisteret al, 1974; 鄭浚茂等, 1980; 袁靜等, 2011)。三根巖芯的概率累積曲線表現(xiàn)形式一致, 包括三種情況: 一段式、兩段式、三段式(圖 4)。一段式主要出現(xiàn)于巖芯中的平均粒徑較細(xì)的常態(tài)沉積物中, 以懸浮式搬運(yùn)為主; 兩段式主要出現(xiàn)于巖芯中較粗夾層的沉積物中, 分為兩個(gè)搬運(yùn)次總體, 即跳躍次總體和懸浮次總體, 截點(diǎn)在5—6Φ之間; 三段式見于巖芯的個(gè)別層位之中, 除了懸浮搬運(yùn)次總體外, 粗粒段出現(xiàn)典型的推移載荷, 推移次總體和懸浮次總體的截點(diǎn)在4Φ附近。

為進(jìn)一步揭示該區(qū)沉積物的機(jī)械搬運(yùn)形式, 對(duì)該處巖芯沉積物粒度的C-M圖式進(jìn)行分析(圖5)。在粒度的C-M圖中, C值與樣品中最粗顆粒的粒徑相當(dāng),代表了水動(dòng)力攪動(dòng)開始搬運(yùn)沉積物的最大能量; M值代表了水動(dòng)力的平均能量。典型的牽引流型 C-M圖包含N-O-P-Q-R-S等線段, 代表滾動(dòng)、遞變懸浮、均勻懸浮等沉積物搬運(yùn)形式(Passega, 1977)。因?yàn)楸緟^(qū)三根巖芯長度都在 250cm以下, 據(jù)該區(qū)的沉積速率估算它們的形成時(shí)間都在最近250年以內(nèi), 而這期間該處沉積環(huán)境基本穩(wěn)定(Yanget al, 1994; 郭志剛等,2000, 2003; 肖尚斌等, 2005a), 所以以巖芯不同層位樣品集合形成C-M圖解(圖4)。三個(gè)巖芯的C-M圖具有基本相同的形式, 各個(gè)巖芯樣品的散點(diǎn)落在3個(gè)區(qū)域, 即圖中的L-1、L-2、L-3。L-1的C、M值呈現(xiàn)同步變化, 形態(tài)大致平行于 C=M 線, 為遞變懸浮搬運(yùn)次總體。L-2樣品點(diǎn)較集中, C值變化部明顯, 因?yàn)樗o鄰L-1, 且M值比L-1小, 分析認(rèn)為屬于均勻懸浮搬運(yùn)次總體。相比于L-1、L-2, L-3樣品點(diǎn)較分散, 特點(diǎn)是M值變化不大, 并且M值位于均勻懸浮次總體的 M 值范圍之內(nèi), 分析認(rèn)為該次總體屬于細(xì)顆粒沉積物中所含的少量的粗顆粒, 它們屬于生物碎屑, 以懸浮或者推移方式進(jìn)行搬運(yùn)。需要注意的是D1、D2、D3三個(gè)巖芯遞變懸浮次總體、均勻懸浮次總體的分界分別在 14、12和 12μm, 與敏感粒級(jí) 1、2的分界一致; 三個(gè)巖芯的均勻懸浮次總體和生物碎屑的分界分別在200、200、和100μm左右, 也與敏感粒級(jí)2、3的分界基本一致。

3.3.2 敏感性粒級(jí)的成因 根據(jù)前文分析可知,研究區(qū)沉積物組成分為三個(gè)敏感粒級(jí),粒級(jí) 1(0—12μm 或 0—14μm)、粒級(jí) 2(12—200μm 或 14—100μm)和粒級(jí)3(大于200μm或大于100μm), 接下來對(duì)這三個(gè)敏感粒級(jí)的成因進(jìn)行討論。

敏感粒級(jí) 1: 它構(gòu)成了本區(qū)沉積物的主體部分,粒級(jí)區(qū)間與概率累積曲線上的懸浮載荷相對(duì)應(yīng)(圖6)。該粒級(jí)是本區(qū)沉積物中最細(xì)的組成部分, 概率累積曲線圖指示其為懸浮載荷, C-M圖進(jìn)一步證實(shí)該沉積物屬于均勻懸浮載荷, 它可以隨著水體進(jìn)行長距離搬運(yùn), 海洋動(dòng)力弱的海況下有利于該粒級(jí)沉積物的沉積。夏季段, 浙閩沿岸流衰退, 風(fēng)浪能量弱, 有利于該粒級(jí)沉積物的沉積和保存; 冬季時(shí)段, 雖然沿岸流和風(fēng)浪作用強(qiáng)烈, 但是由于絮凝作用、生物過程以及海洋鋒面過程等的存在, 也能促使該粒級(jí)的沉積物發(fā)生沉積。

敏感粒級(jí) 2: 其含量僅次于敏感粒級(jí) 1, 且在巖芯中波動(dòng)最明顯, 粒級(jí)區(qū)間與概率累積曲線上的推移載荷(兩段式)以及懸移載荷(一段式)相對(duì)應(yīng)(圖6)。雖然概率累積曲線指示該粒級(jí)既有推移載荷又有懸移載荷, 但是從 C-M 圖顯示該粒級(jí)屬于遞變懸浮搬運(yùn)次總體, 說明這部分沉積物在水柱中的含量是變化的, 并且存在與底質(zhì)沉積物交換的可能。該粒級(jí)沉積物概率累積曲線斜率大、分選好, 反映水動(dòng)力能量較大。該海域冬季時(shí)期受到來自西伯利亞冷空氣的影響, 風(fēng)浪大, 同時(shí)冬季浙閩沿岸流強(qiáng)勁, 可以認(rèn)為該部分沉積物主要為冬季沿岸流搬運(yùn)并沉積而成的。敏感粒級(jí)2既出現(xiàn)于平均粒徑較粗的層位中, 也出現(xiàn)于平均粒徑較細(xì)的常態(tài)層位之中, 可以認(rèn)為該粒級(jí)為該區(qū)正常海況環(huán)境下的沉積產(chǎn)物, 這與肖尚斌等(2005a)認(rèn)為大于45μm的粗粒組分為風(fēng)暴流攜帶的沉積物的認(rèn)識(shí)有所不同。

敏感粒級(jí) 3: 它是本區(qū)最粗的沉積物, 出現(xiàn)在平均粒徑較細(xì)的沉積物之中, 粒級(jí)區(qū)間與概率累積圖中的推移載荷相對(duì)應(yīng)(圖 6)。它在概率累積曲線上表現(xiàn)為典型的推移載荷, 但是在 C-M 圖上則表現(xiàn)為近平行于C軸的較寬的帶狀區(qū)域, 與牽引流的推移載荷段明顯不同。聯(lián)系到該海域具有高的初級(jí)產(chǎn)生力(王亮, 2014), 硅藻、甲藻以及底棲生物發(fā)育, 在分樣中也見到保存程度不一的生物碎片, 為此認(rèn)為該敏感粒級(jí)主要屬于生物成因。由于是生物來源, 它可以直接來自該處底棲、浮游生物的硬體, 或者來自周邊再懸浮生物碎片, 它們與該處的水動(dòng)力條件不一致。

圖4 三根巖芯沉積物粒度概率累積曲線圖Fig.4 The probability cumulative curves of sediment grain sizes of the cores

圖5 三根巖芯沉積物粒度C-M圖Fig.5 C-M diagrams of grain sizes of the cores

圖6 概率累積曲線與粒徑-標(biāo)準(zhǔn)偏差曲線對(duì)比圖Fig.6 Comparison between probability cumulative curve and standard deviation-grain size

4 結(jié)論

東海內(nèi)陸架現(xiàn)代沉積物發(fā)育大致相同的三個(gè)敏感性粒級(jí), 分別是敏感粒級(jí) 1(＜12μm)、敏感粒級(jí)2(12—225μm)和敏感粒級(jí) 3(＞225μm)。從北到南敏感粒級(jí)區(qū)間具有逐漸向細(xì)顆粒方向偏移的趨勢(shì), 以敏感性粒級(jí)2、3的界線的偏移較明顯, 敏感性粒級(jí)1、2的界線偏移幅度很小。這種偏移現(xiàn)象與長江入海沉積物向南搬運(yùn)過程中發(fā)生的沉積分異作用有關(guān)。

敏感粒級(jí)1和2分別由均勻懸浮次總體、遞變懸浮次總體構(gòu)成, 代表不同水動(dòng)力環(huán)境下的搬運(yùn)沉積的結(jié)果。而敏感粒級(jí)3則由生物過程產(chǎn)生, 難以反映水動(dòng)力條件。結(jié)合該區(qū)的海洋動(dòng)力特征, 敏感粒級(jí) 2可較好地指示冬季海洋動(dòng)力的強(qiáng)弱, 并進(jìn)而揭示東亞冬季風(fēng)強(qiáng)度的潛在意義。

王 亮, 2014. 東海典型泥質(zhì)區(qū)高分辨沉積記錄及其對(duì)氣候環(huán)境變化的響應(yīng). 青島: 中國海洋大學(xué)博士學(xué)位論文,101—105

向 榮, 楊作升, Saito Y等, 2006. 濟(jì)州島西南泥質(zhì)區(qū)近2300a來環(huán)境敏感粒度組分記錄的東亞冬季風(fēng)變化. 中國科學(xué)D輯 地球科學(xué), 36(7): 654—662

劉升發(fā), 石學(xué)法, 劉焱光等, 2010. 中全新世以來東亞冬季風(fēng)的東海內(nèi)陸架泥質(zhì)沉積記錄. 科學(xué)通報(bào), 55(14):1387—1396

孫東懷, 安芷生, 蘇瑞俠等, 2001. 古環(huán)境中沉積物粒度組分分離的數(shù)學(xué)方法及其應(yīng)用. 自然科學(xué)進(jìn)展, 11(3): 269—276

李云海, 陳 堅(jiān), 黃財(cái)賓等, 2010. 浙閩沿岸南部泥質(zhì)沉積中心表層沉積物粒度特征及其季節(jié)性差異. 沉積學(xué)報(bào), 28(1):150—157

李家彪, 2008. 東海區(qū)域地質(zhì). 北京: 海洋出版社, 304—390

楊作升, 郭志剛, 王兆祥等, 1992. 黃東海陸架懸浮體向其東部深海區(qū)輸送的宏觀格局. 海洋學(xué)報(bào), 14(2): 81—90

肖尚斌, 李安春, 陳木宏等, 2005a. 近8ka東亞冬季風(fēng)變化的東海內(nèi)陸架泥質(zhì)沉積記錄. 地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報(bào), 30(5): 573—581

肖尚斌, 李安春, 2005b. 東海內(nèi)陸架泥區(qū)沉積物的環(huán)境敏感粒度組分. 沉積學(xué)報(bào), 23(1): 122—129

張存勇, 馮秀麗, 2009. 連云港近岸海域沉積物粒度空間分布特征及其分析. 海洋學(xué)報(bào), 31(4): 120—127

陳國成, 鄭洪波, 李建如等, 2007. 南海西部陸源沉積粒度組成的控制動(dòng)力及其反映的東亞季風(fēng)演化. 科學(xué)通報(bào),52(23): 2768—2776

鄭浚茂, 王德發(fā), 孫永傳, 1980. 黃驊拗陷幾種砂體的粒度分布特征及其水動(dòng)力條件的初步分析. 石油實(shí)驗(yàn)地質(zhì), 2(2):9—20

秦蘊(yùn)珊, 鄭鐵民, 1982. 東海大陸架沉積物分布特征的初步探討. 見: 秦蘊(yùn)珊主編. 黃東海地質(zhì). 北京: 科學(xué)出版社,31—51

秦蘊(yùn)珊, 趙一陽, 陳麗蓉等, 1987. 東海地質(zhì). 北京: 科學(xué)出版社, 29—34

袁 靜, 楊學(xué)君, 路智勇等, 2011. 東營凹陷鹽22塊沙四上亞段砂礫巖粒度概率累積曲線特征. 沉積學(xué)報(bào), 29(5):815—824

郭志剛, 楊作升, 1999. 東海陸架北部泥質(zhì)區(qū)沉積動(dòng)力過程的季節(jié)性變化. 青島海洋大學(xué)學(xué)報(bào), 29(3): 507—513

郭志剛, 楊作升, 曲艷慧等, 2000. 東海陸架泥質(zhì)區(qū)沉積地球化學(xué)比較研究. 沉積學(xué)報(bào), 18(2): 284—289

郭志剛, 楊作升, 范德江等, 2003. 長江口泥質(zhì)區(qū)的季節(jié)性沉積效應(yīng). 地理學(xué)報(bào), 58(4): 591—597

Boulay S, Colin C, Trentesaux A et al, 2003. Mineralogy and sedimentology of Pleistocene sediment in the South China Sea (ODP site 1144). In: Prell W L, Wang P, Blum P et al eds. Proceedings of the Ocean Drilling Program. Scientific Results. 184: 1—21

Clark M W, 1976. Some methods for statistical analysis of multimodal distributions and their application to grain-size data. Journal of the International Association for Mathematical Geology, 8(3): 267—282

Davis J C, Sampson R J, 1975. Statistics and data analysis in geology. The Journal of Geology, 83(1): 141—142

Fan D J, Qi H Y, Sun X X et al, 2011. Annual lamination and its sedimentary implications in the Yangtze River delta inferred from high-resolution biogenic silica and sensitive grain-size records. Continental Shelf Research, 31(2): 129—137

Folk R L, Andrews P B, Lewis D W, 1970. Detrital sedimentary rock classification and nomenclature for use in New Zealand.New Zealand Journal of Geology and Geophysics, 13(4):937—968

Folk R L, Ward W C, 1957. Brazos River bar: A study in the significance of grain size parameters. Journal of Sedimentary Research, 27(1): 3—26

Glaister R P, Nelson H W, 1974. Grain-size distributions, an aid in facies identification. Bulletin of Canadian Petroleum Geology, 22(3): 203—240

Hu B Q, Yang Z S, Zhao M X et al, 2012. Grain size records reveal variability of the East Asian Winter Monsoon since the Middle Holocene in the Central Yellow Sea mud area,China. Science China Earth Sciences, 55(10): 1656—1668

Huang J, Li A C, Wan S M, 2011. Sensitive grain-size records of Holocene East Asian summer monsoon in sediments of northern South China Sea slope. Quaternary Research, 75(3):734—744

Krumbein W C, 1934. Size frequency distributions of sediments.Journal of Sedimentary Research, 4(2): 65—77

Liu J, Saito Y, Kong X H et al, 2010a. Sedimentary record of environmental evolution off the Yangtze River estuary, East China Sea, during the last ～13, 000 years, with special reference to the influence of the Yellow River on the Yangtze River delta during the last 600 years. Quaternary Science Reviews, 29(1718): 2424—2438

Liu J P, Li A C, Xu K H et al, 2006. Sedimentary features of the Yangtze River-derived along-shelf clinoform deposit in the East China Sea. Continental Shelf Research, 26(17—18):2141—2156

Liu J P, Xu K H, Li A C et al, 2007. Flux and fate of Yangtze River sediment delivered to the East China Sea.Geomorphology, 85(3—4): 208—224

Liu S F, Shi X F, Liu Y G et al, 2010b. Records of the East Asian winter monsoon from the mud area on the inner shelf of the East China Sea since the mid-Holocene. Chinese Science Bulletin, 55(21): 2306—2314

McManus J, 1988. Grain size determination and interpretation.Techniques in Sedimentology, 408: 63—85

Passega R, 1957. Texture as characteristic of clastic deposition.AAPG Bulletin, 41(9): 1952—1984

Passega R, 1977. Significance of CM diagrams of sediments deposited by suspensions. Sedimentology, 24(5): 723—733

Prins M A, Postma G, Weltje G J, 2000. Controls on terrigenous sediment supply to the Arabian Sea during the late Quaternary: The Makran continental slope. Marine Geology,169(3—4): 351—371

Saito Y, Yang Z S, 1995. Historical change of the Huanghe(Yellow River) and its impact on the sediment budget of the East China Sea. In: Tsunogai S, Iseki K, Koike I et al eds.Proceedings of International Symposium on Global Fluxes of Carbon and its Related Substances in the Coastal Sea-Ocean Atmosphere System. Sapporo, Japan: Hokkido University Press, 7—12

Shepard F P, 1954. Nomenclature based on sand-silt-clay ratios.Journal of Sedimentary Research, 24(3): 151—158

Su J L, Pan Y Q, 1987. On the shelf circulation north of Taiwan.Acta Oceanologica Sinica, 6(S1): 1—20

Syvitski J P M, 1991. Factor analysis of size frequency distributions: Significance of factor solutions based on simulation experiments. In: Syvitski J ed. Principles,Methods, and Application of Particle Size Analysis.Cambridge, UK: Cambridge University Press, 249—263

Visher G S, 1969. Grain size distributions and depositional processes.Journal of Sedimentary Research, 39(3): 1074—1106

Wang Y J, Cheng H, Edwards R L et al, 2005. The Holocene Asian monsoon: Links to solar changes and North Atlantic climate. Science, 308(5723): 854—857

Weltje G J, 1997. End-member modeling of compositional data:Numerical-statistical algorithms for solving the explicit mixing problem. Mathematical Geology, 29(4): 503—549

Yang Z S, Saito Y, Guo Z G et al, 1994. Distal mud area as a material sink in the East China Sea. In: Tsunogai S, Iseki K,Koike I et al eds. Proceedings of International Symposium on Global Fluxs of Carbon and Its Related Substances in the Coastal Sea-ocean Atmosphere System. Sapporo, Japan:Hokkido University Press, 1—6