商志文，田立柱，李建芬，陳永勝，王宏

（中國地質(zhì)調(diào)查局 天津地質(zhì)調(diào)查中心，天津 300170）

作為深入內(nèi)陸的陸表海灣，渤海灣受到海陸交互作用的強(qiáng)烈影響。晚更新世以來，渤海灣及相鄰沿海低地經(jīng)歷了3次大的海陸交替過程。其中，全新世海侵地層保存相對完整，信息記錄最為豐富，是開展地質(zhì)環(huán)境演變研究的良好對象。但以往全新世海侵時(shí)的海陸作用研究多是從地貌演化和史料記載入手，僅近10余年來筆者所在研究小組對該地區(qū)地層與古環(huán)境演變進(jìn)行了新一輪研究（李鳳林等，2004；李建芬等，2004；閻玉忠等，2006；范昌福等，2005、2008；商志文等 ， 2010a、 b；Fan et al， 2011；陳 永 勝 等 ，2012a、b）。

硅藻是海洋中最重要的微體植物之一，因水體物理、化學(xué)及水動(dòng)力條件的差異，其種類、數(shù)量及組合特征存在明顯的多樣性。因此，可以通過沉積物中化石硅藻的特征推斷古海洋環(huán)境和古氣候條件。近年來，已有眾多海域及海陸過渡相化石硅藻用于第四紀(jì)古環(huán)境重建的研究實(shí)例（王開 發(fā) 等 ， 1995、 2002a、 b、 2003；Sato et al，2001、2003；支崇遠(yuǎn)等，2003；李超等，2004；Freund et al，2004；黃玥等，2006；冉利華等，2008；李冬玲 等，2009；Fukumoto，2011；Saegusa et al，2011）。但渤海灣地區(qū)，硅藻與古環(huán)境的研究仍是極少涉及的領(lǐng)域（商志文等，2010a、b）。

本文試圖通過CH114孔（圖1）的沉積學(xué)和微體古生物學(xué)研究，劃分沉積相，進(jìn)而結(jié)合年代學(xué)數(shù)據(jù)，揭示渤海灣西北部淺海區(qū)全新世古環(huán)境演變、海洋影響與沉積速率特征，并討論它們對環(huán)境演化的貢獻(xiàn)。

圖1 CH114位置圖

1 材料和方法

CH114孔位于渤海灣西北部、岸外約3 km的淺海區(qū)，是天津漢沽海域活動(dòng)斷層地球物理調(diào)查的驗(yàn)證孔之一。巖心反映了渤海灣西北部近岸淺海區(qū)淺地層物質(zhì)組成、層序與地質(zhì)結(jié)構(gòu)。鉆孔采用機(jī)械鉆探，全取心，巖心直徑108 mm，平均采取率＞90%。鉆探進(jìn)尺依靠測深儀進(jìn)行實(shí)時(shí)潮位水深校正，以保證進(jìn)尺的準(zhǔn)確性。由測深儀測量并經(jīng)潮位校正后計(jì)算的孔口高程是-3.9 m（國家85高程）。CH114孔全新統(tǒng)厚16 m，為本文的研究對象。

CH114孔9個(gè)測年樣品均采自巖心中相對完整的貝殼或泥炭，并在北京大學(xué)考古文博學(xué)院AMS（加速器質(zhì)譜）14C實(shí)驗(yàn)室進(jìn)行測試。測定值以5568半衰期計(jì)年，測定結(jié)果按照國際14C委員會(huì)通行方法進(jìn)行分餾效應(yīng)、海洋貯存庫效應(yīng)和大氣效應(yīng)校正后以日歷年齡表達(dá)（cal BP）。

依據(jù)地層變化，CH114孔全新世沉積物共采取硅藻分析樣品32個(gè)，取樣間距10～60 cm。樣品按照常規(guī)方法處理（Richard，1986），過程依次為：30%雙氧水去除有機(jī)質(zhì)、10%鹽酸去除鈣質(zhì)、比重為2.38的ZnBr2重液浮選兩次、加拿大樹膠封片。制成玻片后，一般在400倍（個(gè)體較小的屬種用1 000倍的油鏡）顯微鏡下鑒定屬種，每個(gè)樣品鑒定硅藻殼體約200個(gè)。

粒度分析方法為激光粒度分析。粒級標(biāo)準(zhǔn)為尤登-溫德華氏等比制φ值粒級標(biāo)準(zhǔn)；平均粒徑、分選系數(shù)、偏態(tài)和峰態(tài)等粒度參數(shù)運(yùn)用?？?沃德公式計(jì)算。

2 結(jié)果

2.1 年代學(xué)與沉積速率特征

9個(gè)AMS14C年齡由下向上表現(xiàn)出依次年輕的自恰性，且2組取自相同層位的測年結(jié)果基本一致（表1、圖2），與相鄰鉆孔CH19孔相同沉積層位的年齡具有較好的可比性（商志文等，2010a），結(jié)果可信，顯示該段沉積始于全新世（表1、圖2）。根據(jù)AMS14C年齡結(jié)果，表2給出了各時(shí)段的沉積速率。需要說明的是：（1）沉積速率計(jì)算時(shí)，10.9 m和11.2 m處的年齡，為該處取自同一位置的兩個(gè)樣品年齡數(shù)據(jù)的平均值。（2）12 m以下缺乏年齡數(shù)據(jù)，因此，全新世底界參考與該孔相距約3 km的CH19孔（圖1）年齡數(shù)據(jù)11 180 cal BP（商志文等，2010a）作為該孔全新世的開始（表2中用斜體表示參考年齡數(shù)據(jù)及相應(yīng)的沉積速率）。

表1 CH114孔AMS14C測年結(jié)果

圖2 CH114深度-年齡分布圖

2.2 硅藻分布與組成

對CH114孔巖心的系統(tǒng)分析研究表明，除個(gè)別層位硅藻貧乏外，大部分樣品硅藻比較豐富，32個(gè)樣品中共發(fā)現(xiàn)17屬28種。

表2 渤海灣西北岸CH114孔沉積速率

從每個(gè)樣品中選取含量最高的前3種硅藻作為該樣品的優(yōu)勢硅藻種，CH114孔的主要優(yōu)勢種有10種：Cyclotella striata/stylorum，Actinoptychus undulatus，Coscinodiscus argus， Coscinodiscus perforatus， Surirella armoricana， Coscinodiscus radiatus， Paralia sulcata， Coscinodiscus subconcavus，Grammatophora oceanica 和 Thalassionema nitzschioides。

Cyclotella striata/stylorum為潮間帶-沿岸種，潮間帶地區(qū)含量高達(dá)30%以上，最多達(dá)80%，自岸向海數(shù)量減少，水深大于30 m海區(qū)，數(shù)量一般小于15%。這兩個(gè)種在中國北黃海沿岸、黃海、東海、南海北部灣等地均有記錄（蔣輝，1987），為研究區(qū)第一優(yōu)勢硅藻種（商志文等，2012）。

Actinoptychus undulatus為廣布性的沿岸底棲種，該種的最高含量出現(xiàn)在南黃海和東海水深40 m左右的淺海區(qū)，含量大于30%，大于或小于這個(gè)水深，數(shù)量明顯減少（蔣輝，1987）。在渤海灣潮間帶-淺海區(qū)表層沉積物中，潮間帶地區(qū)含量低于淺海區(qū)（商志文，2012）。

Coscinodiscus argus、Coscinodiscus radiatus 為廣布性浮游種，在渤海灣潮間帶-淺海區(qū)表層沉積物中，含量隨水深的變化不明顯（商志文等，2012）。

Coscinodiscus perforatus、 Surirella armoricana，Coscinodiscus subconcavus均為沿岸底棲種。在渤海灣潮間帶-淺海區(qū)表層沉積物中，Coscinodiscus perforatus在海河北部潮間帶至2 m水深的淺海區(qū)含量最高。Surirella armoricana、Coscinodiscus subconcavus主要分布在淺海區(qū)，且隨著水深加大，含量逐漸增多（商志文等，2012）。

Paralia sulcata為典型的淺海種，是我國表層沉積物中分布最廣泛的種類之一，從岸線到淺海，數(shù)量逐漸增加，而從淺海到深海區(qū)或海槽，含量則逐漸減少，水深50～100 m最適合其生長（蔣輝，1987）。在渤海灣潮間帶-淺海區(qū)表層沉積物中，該種主要分布在淺海區(qū)，且隨著水深的增加，數(shù)量逐漸增多（商志文等，2012）。

Thalassionema nitzschioides為亞熱帶浮游種，分布區(qū)域鹽度較高，代表海水影響較強(qiáng)的沉積環(huán)境，是渤海灣10～20 m水深表層沉積硅藻優(yōu)勢種之一（商志文等，2012）。

2.3 海洋影響分析與古環(huán)境重建

利用Tilia軟件對CH114孔的10個(gè)優(yōu)勢硅藻種進(jìn)行聚類分析（圖3），共劃分出6個(gè)硅藻組合帶。同時(shí)，本文將亞熱帶浮游種Thalassionema nitzschioides與潮間帶沿岸種Cyclotella striata/stylorum含量的比值，作為判定海洋影響強(qiáng)弱的指標(biāo)（MI），比值越大，表明海洋影響越強(qiáng)（MI為0，表明Thalassionema nitzschioides的缺失，但并不代表完全沒有海洋影響，只是海洋作用相對較弱而已）。各帶的硅藻組合特征及古環(huán)境如下：

組合1帶，埋藏深度為16～11.2 m。其中，16～14.3 m，灰黃棕色粉砂，夾泥質(zhì)條帶和碳質(zhì)條帶，波狀復(fù)合層理，底部為淡水貝殼的富集層。本段未見硅藻，為陸相環(huán)境。14.3～11.2 m，暗橄欖棕與灰黃棕色粘土質(zhì)粉砂，夾薄層炭質(zhì)條帶和炭斑，發(fā)育生物潛穴。本段只見零星的圓篩藻屬硅藻碎片，代表了受海水影響的鹽沼低地環(huán)境。

組合2帶，埋藏深度為11.2～10.4 m。黑棕色貝殼碎屑與粉砂質(zhì)砂的混雜堆積，無層理，與下伏地層侵蝕接觸。本帶硅藻屬種相對較少，以Cyclotella striata/stylorum-Coscinodiscus perforatus-Coscinodiscus subconcavus為組合特征，3個(gè)優(yōu)勢種的含量＞80%，見Actinoptychus undulatus，Coscinodiscus argus，Coscinodiscus radiatus，Thalassionema nitzschioides和Grammatophora oceanica等。該帶沉積物顯示環(huán)境較為動(dòng)蕩，硅藻保存較差，屬種相對單一，硅藻組合并不能完全反映當(dāng)時(shí)的沉積環(huán)境。年代學(xué)數(shù)據(jù)顯示該段地層沉積于6 646～4 280 cal BP（表1、圖2），當(dāng)時(shí)海平面與現(xiàn)在相當(dāng)或略高（Chpell et al.，1996），海岸線位于現(xiàn)代岸線以西約40～50 km處（李世瑜，1962；王穎，1964；趙希濤，1980），而CH114孔位于現(xiàn)代岸線以東約3 km的淺海區(qū)，因此該孔在當(dāng)時(shí)應(yīng)位于遠(yuǎn)離岸線的淺海區(qū)，水深約15 m（孔口高程-3.9 m+埋深約10.8 m）。由于遠(yuǎn)離岸線，機(jī)械沉積速率緩慢，為軟體動(dòng)物繁盛提供了相對適宜的環(huán)境。另一方面，這一深度仍然位于壞天氣浪基面以上，該貝殼碎屑層分選差、不具層理、硅藻保存不良，充分顯示了動(dòng)蕩的水動(dòng)力特征。據(jù)此推斷該段為水深較大的淺海環(huán)境，且貝屑混雜堆積可能受到風(fēng)暴事件影響。

圖3 CH114孔柱狀地層剖面、硅藻組合特征、年齡與海洋影響變化

組合3帶，埋藏深度為10.4～9.6 m。黑棕色粘土質(zhì)粉砂，偶見貝殼及貝殼碎片。本帶硅藻屬種豐富，以Thalassionema nitzschioides-Cyclotella striata/stylorum-Coscinodiscus perforatus為組合特征，三者的平均含量分別為35%、25%和13%。另見Coscinodiscus argus、Coscinodiscus radiatus、Coscinodiscussubconcavus、 Paralia sulcata 和 Grammatophora oceanica等。本帶中Thalassionema nitzschioides為第一優(yōu)勢種，且潮間帶沿岸種含量在地層中所見硅藻的各組合帶中最低，指示該帶為受海水影響較強(qiáng)（MI達(dá)到最大峰值）、水深較大的淺海環(huán)境，推測水深＞10 m。

組合4帶，埋深為9.6～4.7 m。粒度向上逐漸變粗，下部為黑棕色粘土質(zhì)粉砂，含砂質(zhì)透鏡體。中部為黑棕色粉砂與粘土質(zhì)粉砂互層，具波狀復(fù)合層理。上部為黃棕色粉細(xì)砂，發(fā)育小型交錯(cuò)層理與生物潛穴。本帶硅藻屬種豐富，以Cyclotella striata/stylorum-Coscinodiscus perforatus-Thalassionema nitzschioides為組合特征，三者的平均含量分別為28%、23%和9%。另見Actinoptychus undulatus、 Coscinodiscus argus、 Coscinodiscus radiatus、Coscinodiscus subconcavus、Paralia sulcata，Surirella armoricana和Grammatophora oceanica等。本帶硅藻組合以淺海種為主，潮間帶沿岸種Cyclotella striata/stylorum和Coscinodiscus perforatus的含量較組合3帶略有增加，Thalassionema nitzschioides仍為本帶的優(yōu)勢種之一，指示了較強(qiáng)的海水影響（MI出現(xiàn)3次峰值），綜合判斷為淺海環(huán)境，推測水深7～10 m，鹽度較組合3帶降低。

組合5帶，埋藏深度為4.7～1.65 m。黑棕色粉砂，粒度向上變粗，含大量的貝殼碎屑、砂團(tuán)及砂質(zhì)條帶，生物擾動(dòng)明顯。本帶硅藻種屬豐富，以Coscinodiscus perforatus-Cyclotella striata/stylorum-Coscinodiscus argus為組合特征，三者的平均含量分別為34%、29.5%和13%；另見Coscinodiscus subconcavus（3.5%）、Coscinodiscus radiatus（11%）和Grammatophora oceanica（6%）等。與組合4帶相比，Cyclotella striata/stylorum和沿岸種Coscinodiscus perforatus含量增加，淺海種Paralia sulcata和亞熱帶浮游種Thalassionema nitzschioides消失。綜合判斷該帶為淺海環(huán)境，推測水深范圍為2～5 m，海洋影響減弱，鹽度進(jìn)一步降低。

組合6帶，埋藏深度為1.65～0 m。濁黃棕色粘土質(zhì)粉砂，發(fā)育生物潛穴，含砂質(zhì)透鏡體與貝殼碎屑。本帶硅藻屬種豐富，以Cyclotella striata/stylorum -Coscinodiscus perforatus-Coscinodiscus subconcavus為組合特征，3個(gè)優(yōu)勢種的含量＞80%，見 Actinoptychus undulatus、 Coscinodiscus argus、Coscinodiscus radiatus和Grammatophora oceanica等。硅藻組合以潮間帶沿岸種為主，含淺海種與少量浮游種，判斷為近岸淺海環(huán)境，根據(jù)優(yōu)勢硅藻種在渤海灣表層沉積物中的分布特征，推測該帶的水深范圍為0～2 m，鹽度進(jìn)一步降低。

3 討論

3.1 沉積環(huán)境與沉積速率變化

CH114孔揭示該孔全新世以來經(jīng)歷了陸相→受海水影響的鹽沼低地環(huán)境→受風(fēng)暴強(qiáng)力事件影響的淺海環(huán)境→淺海環(huán)境的演化過程。

全新世初期，海水尚未到達(dá)CH114孔所在地，仍為陸相環(huán)境。隨著海面快速上升，經(jīng)歷了受海水影響較弱的鹽沼低地環(huán)境后（埋深12 m處的年齡為9 455 cal BP），轉(zhuǎn)為淺海環(huán)境。之后海面快速上升，物源貧乏，9 455～3 634 cal BP期間機(jī)械沉積速率低，特別是 6 646～4 280 cal BP的約2 300年間，僅沉積了0.8 m厚的受風(fēng)暴事件影響的貝殼碎屑與泥砂混雜沉積物，平均速率為0.03 cm/a。3 634 cal BP之后，仍為淺海環(huán)境，但水深逐漸減小，直至現(xiàn)代的1.5 m。其中，3 634～2 641 cal BP時(shí)，沉積速率較大，平均為0.5 cm/a；2 641～424 cal BP，沉積速率減小，僅為 0.12 cm/a；424 cal BP以后，平均速率為0.38 cm/a（圖4）。

圖4 渤海灣CH114鉆孔沉積速率變化

沉積速率的階段性變化與河流供給有較好的對應(yīng)關(guān)系。全新世初期，海水尚未到達(dá)CH114孔所在地，河流供給物源相對充足，為沉積速率相對較高時(shí)期。海水約9 ka cal BP到達(dá)渤海灣西岸后，海面持續(xù)上升，岸線迅速后退，形成了巨大的可容空間，但9～4 ka cal BP期間的海面上升速率遠(yuǎn)高于這一階段的沉積速率。CH114孔在9.5～3.6 ka cal BP為沉積速率最低時(shí)期，說明在此期間盡管有足夠的可容空間，但沉積供給僅維持了低水平的相對穩(wěn)定，河流作用不明顯。距今約3 ka cal BP，中國東部相對海面已經(jīng)穩(wěn)定于現(xiàn)代海面高度（沈明潔，2002；Zong，2004）。3.6～2.6 ka cal BP期間，平均沉積速率上升至0.5 cm/a，距離CH114孔最近的黃河三角洲超級葉瓣5發(fā)育于 3.0～2.5 ka BP（Xue et al，1993），該時(shí)段古黃河三角洲的發(fā)育與CH114孔沉積速率的增加在時(shí)空上基本對應(yīng)，說明CH114孔在該時(shí)段沉積速率的增大可能與黃河輸入加強(qiáng)有關(guān)。2.6～0.4 ka cal BP期間，平均沉積速率下降至0.12 cm/a，揭示該時(shí)段河流輸入作用較弱。0.4 ka cal BP以來，平均沉積速率上升至0.38 cm/a。公元1194年，黃河南徙，海河得以在其現(xiàn)代河口處獨(dú)立入海（黃盛璋，1982），該時(shí)段沉積速率的增加與海河獨(dú)立入海在時(shí)空上基本對應(yīng)，揭示其可能與海河輸入增強(qiáng)有關(guān)。

3.2 海洋影響曲線與全新世5 000 cal BP以來氣候變化的對比

CH114孔海洋影響曲線特征是（圖3、5）：

（1）埋深10.4～11.2 m，為貝殼碎屑與粉砂質(zhì)砂的混雜堆積，MI＜0.07，表明海洋影響較弱，但是，此時(shí)沉積環(huán)境不穩(wěn)定，硅藻保存相對較差，屬種單一。因此，硅藻組合并不能完全反映當(dāng)時(shí)的沉積環(huán)境。埋深10.2 m粗貝殼碎屑層之上的泥質(zhì)沉積層底部MI達(dá)到最大峰值，指示海洋作用很強(qiáng)。之后向上至埋深7.2 m處MI整體上呈波動(dòng)減小的趨勢。埋深6 m處MI增大，表明一次海水影響增強(qiáng)的短暫過程。之后，又一次呈波動(dòng)減弱趨勢；埋深4.7 m以上，MI為0，指示海洋影響更弱。

（2）依據(jù)已有14C年齡采用線性內(nèi)插方法，賦予CH114孔巖心硅藻-海洋影響曲線以時(shí)間意義，便獲得了該孔硅藻所記錄的全新世5 000 cal BP以來海洋影響強(qiáng)弱的波動(dòng)曲線，曲線反映存在4次海洋影響增強(qiáng)的時(shí)期（圖3、5）。

前人多認(rèn)為全新世5000 cal BP以來共有4次氣候變冷事件（事件3～0），但不同研究者對于各次事件的起迄時(shí)間尚未完全統(tǒng)一。本文綜合了前人的研究成果（Hong et al，2009；王紹武，2009；Mayewski et al， 2004；Bond et al， 1997；Wang，1994；Wang et al，1997，2010；葛全勝 等，2006；竺可楨，1973；Ge et al，2010；鄭景云等，2005；王宏，1996），重點(diǎn)參考了與本研究區(qū)相鄰的中國東北部、東部地區(qū)各次寒冷事件的起迄時(shí)間，將這4次寒冷事件分別定為4 300～3 800 cal BP（事件 3）、3 550～3 350 cal BP 和 3 050～2 900 cal BP（事件2：細(xì)分為2次）、1 550～1 000 cal BP（事件1）及480～90 cal BP（事件0），寒冷事件之間為溫暖期（圖5）。

圖5 CH114孔海洋影響曲線與全新世5 000年以來氣候變化對比圖

海洋影響的波動(dòng)變化與全新世氣候變化對比顯示（圖5），由硅藻記錄的海洋影響增強(qiáng)的4次峰值期多與氣候溫暖期對應(yīng)（峰值1，與相關(guān)溫暖期不同相位，尚無法解釋）。表明研究區(qū)海洋影響的增強(qiáng)與氣候變暖具有一定的正相關(guān)性。但是，二者變化的起迄時(shí)間是否嚴(yán)格一致，仍需今后進(jìn)一步深入探討。

4 結(jié)論

（1）CH114孔巖心全新世硅藻經(jīng)系統(tǒng)研究，共發(fā)現(xiàn)17屬28種，與沉積巖石學(xué)、年代學(xué)（AMS14C）研究相結(jié)合，將該孔自下而上劃分為6個(gè)硅藻組合帶，分別對應(yīng)全新世初期陸相（組合1帶下部）、受海水影響的鹽沼低地環(huán)境（組合1帶上部）、6 646～4 280 cal BP年間受風(fēng)暴強(qiáng)力事件影響的淺海環(huán)境（組合2帶）和4 280年以來水深不斷變淺的淺海環(huán)境（組合3-6帶）。

（2）CH114孔沉積速率自下而上呈高-低-高-低-高的階段性變化特征，這與研究區(qū)的河流供給有較好的對應(yīng)關(guān)系：即3.6～2.6 ka cal BP和0.4 ka cal BP以來兩個(gè)時(shí)段相對較高的沉積速率，分別與黃河三角洲超級葉瓣5的形成和海河獨(dú)立入海相對應(yīng)。

（3）由硅藻記錄的海洋影響的波動(dòng)變化與全新世氣候變化對比顯示：研究區(qū)海洋影響的增強(qiáng)與氣候變暖具有一定的正相關(guān)性。

致謝：北京大學(xué)考古文博學(xué)院AMS（加速器質(zhì)譜）14C實(shí)驗(yàn)室提供AMS14C測年數(shù)據(jù)，野外鉆探由裴艷東、王福和田立柱完成，在此一并致謝。

Bond G，Showers W，Cheseby M，et al，1997.A prevasive millennial-scale cycle in north Atlantic Holocene and glacial climates.Science，278（14）：1257-1266.

Chapell J，Omura A，Esat T，et al，1996.Reconciliation of late Quaternary sea levels derived from coral terraces at Huon Peninsula with deep sea oxygen isotope records.Earth and Planetary Science Letters，141：227-236.

Fan C F，Koniger P，Wang H，et al，2011.Ligament increments of pacific oyster are reliable independent proxies for seasonality in the western Bohai Bea， China.Paleogeography， Paleoclimatology，Paleoecology，299：437-448.

Freund H，Gerdes G，Streif H，et al，2004.The indicative meaning of diatoms，pollen and botanical macrofossils for the reconstruction of palaeoenvironments and sea-level fluctuations along the coast of Lower Saxony；Germany.Quaternary International，112：71-87.

Fukumoto Y，2011.Mid-late Holocene paoeoenvironment in Karako lowland，western Japan，inferred from diatom analysis.Quaternary International，230：115-121.

Ge Q S，Zheng J Y，Hao Z X，et al，2010.Temperature variation through 2000 years in China： an uncertainly analysis of reconstruction and regional differences.Geophysical Research Letters，37：1-5.

Hong Y T， Hong B， Lin Q H， et al， 2009.Synchronous climate anomalies in the western North Pacific and North Atlantic regions during the last 14，000 years.Quaternary Science Reviews，28：840-849.

Mayewski P A，Rohling E E，Stager J C，et al，2004.Holocene climate variability.Quaternary Research，62：243-255.

Richard W，1986.Battarbee.Diatom analysis，Handbook of Holocene Palaeoecology and Palaeohydrology[M].John Wiley&Sons Ltd，527-570.

Saegusa Y， Sugai T， Ogami T， et al， 2011.Reconstruction of Holocene environmental changes in the Kiso-Ibi-Nagara compound river delta，Nobi Plain，central Japan，by diatom analyses.Quaternary International，230：67-77.

Sato H，Okuno J，Katoh S，2003.Holocene crustal movement along the coast of western Kobe and the 1995 Kobe Earthquake，Japan.Quaternary Science Reviews，22：891-897.

Sato H，Okuno J，Nakada M，2001.Holocene uplift derived from relative sea-level records along the coast of western Kobe，Japan.Quaternary Science Reviews，20：1459-1474.

Wang H，1994.Palaeoenvironment of Holocene Chenier and Oyster Reefs in the Bohai Bay（China）.PhD Dissertation，Vrije Universiteit Brussel，1-245.

Wang H，Shang Z W，Li J F，et al，2010.Holocene shoreline changes and marine impacts on the muddy coast，western Bohai Bay，China.Geological Bulletin of China，29（5）：627-640.

Wang H，Van Strydonck M，1997.Chronology of Holocene cheniers and oyster reefs on the coast of Bohai Bay，China.Quaternary Research，47：192-205.

Xue C T，1993.Historical Changes in the Yellow River Delta，China.Marine Geology，113：321-329.

Zong Y Q，2004.Mid-Holocene Sea-Level Highstand along the Southeast Coast of China.Quaternary International，117：55-67.

陳永勝，王宏，李建芬，等，2012b．渤海灣西岸BT113孔35ka以來的沉積環(huán)境演化與海陸作用．吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），42（增刊）：205-216.

陳永勝，王宏，裴艷東，等，2012a．渤海灣西岸晚第四紀(jì)海相地層劃分及地質(zhì)意義．吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），42（3）：652-663.

范昌福，李建芬，王宏，等，2005.渤海灣西北岸大吳莊牡蠣礁測年與古環(huán)境變化.地質(zhì)調(diào)查與研究，28（2）：124-129.

范昌福，王宏，裴艷東，等，2008.渤海灣西北岸濱海湖埋藏牡蠣礁古生態(tài)環(huán)境.海洋地質(zhì)與第四紀(jì)地質(zhì)，28（1）：33-41.

葛全勝，王順兵，鄭景云，2006.過去5000年中國氣溫變化序列重建.自然科學(xué)進(jìn)展，16（6）：689-696.

黃盛璋，1982．歷史時(shí)期的水系變遷，海河．中國科學(xué)院《中國自然地理》編輯委員會(huì)．中國自然地理，歷史自然地理．北京：科學(xué)出版社，38-86.

黃玥，冉莉華，蔣輝，2006.南海北部陸坡晚更新世末期硅藻及其古環(huán)境意義.海洋地質(zhì)與第四紀(jì)地質(zhì)，26（4）：7-13.

蔣輝，1987.我國某些常見化石硅藻的環(huán)境分析.植物學(xué)報(bào)，29（4）：440-448.

李超，藍(lán)東兆，石學(xué)法，等，2004.東海陸架與沖繩海漕晚第四紀(jì)沉積硅藻的對比研究.廈門大學(xué)學(xué)報(bào)（自然科學(xué)版），43（6）：836-841.

李冬玲，范昌福，黃玥，等，2009.渤海灣西北岸中全新世埋藏牡蠣礁的硅藻記錄及古環(huán)境意義.海洋通報(bào)，28（3）：22-27.

李鳳林，王宏，閻玉忠，等，2004.渤海灣西岸濱海平原晚第四紀(jì)以來的沉積間斷.地質(zhì)調(diào)查與研究，27（3）：177-183.

李建芬，王宏，李鳳林，等，2004.渤海灣牡蠣礁平原中部興坨剖面全新世地質(zhì)環(huán)境變遷.地質(zhì)通報(bào)，23（2）：169-176.

李世瑜，1962.古代渤海灣西部海岸遺跡及地下文物的初步調(diào)查研究.考古，12：652-657.

冉莉華，蔣輝，Knudsen K L，等，2008.北大西洋北部中全新世以來古海洋環(huán)境記錄.海洋通報(bào)，27（5）：39-46.

商志文，范昌福，李冬玲，等，2010b.硅藻組合指示的渤海灣西北岸兩個(gè)牡蠣礁體生長環(huán)境的差異.海洋地質(zhì)與第四紀(jì)地質(zhì)，30（5）：33-39.

商志文，田立柱，王宏，等，2010a.渤海灣西北部CH19孔全新統(tǒng)硅藻組合、年代學(xué)與古環(huán)境.地質(zhì)通報(bào)，29（5）：675-681.

商志文，田立柱，王宏，等，2012.渤海灣中北部表層沉積硅藻分布及環(huán)境指示意義.中國地質(zhì)，39（4）：1099-1107.

沈明潔，謝志仁，朱誠，2002．中國東部全新世以來海平面波動(dòng)特征探討．地球科學(xué)進(jìn)展，17（6）：886-894.

王宏，1996.渤海灣全新世貝殼堤和牡蠣礁的古環(huán)境.第四紀(jì)研究，1：71-79.

王穎，1964.渤海灣西北殼堤與古岸線問題.南京大學(xué)學(xué)報(bào)，8（3）：424-440（附圖版3幀）.

王開發(fā)，陸繼軍，鄭玉龍，1995.福建沿岸第四紀(jì)孢粉、硅藻組合及其古環(huán)境意義.微體古生物學(xué)報(bào)，12（4）：388-397.

王開發(fā)，鄭玉龍，支崇遠(yuǎn)，等，2002b.東海南北陸緣（莆、泉段）全新世沉積硅藻.古生物學(xué)報(bào)，41（2）：273-279.

王開發(fā)，支崇遠(yuǎn)，陶明華，2003.東海陸緣（浙南段）晚第四紀(jì)硅藻的發(fā)現(xiàn)及古環(huán)境分析.微體古生物學(xué)報(bào)，20（4）：350-357.

王開發(fā)，支崇遠(yuǎn)，鄭玉龍，等，2002a.東海陸緣（閩北段）晚第四紀(jì)沉積的硅藻學(xué)研究.沉積學(xué)報(bào)，20（1）：135-143.

王紹武，2009.全新世氣候.氣候變化研究進(jìn)展，5（4）：247-248.

閻玉忠，王宏，李鳳林，等，2006.渤海灣西岸BQ1孔揭示的沉積環(huán)境與海面波動(dòng).地質(zhì)通報(bào)，25（3）：357-382.

趙希濤，1980.渤海灣西岸全新世海岸線變遷，華北斷塊區(qū)的形成與發(fā)展.北京：科學(xué)出版社：302-309.

鄭景云，滿志敏，方修琦，等，2005.魏晉南北朝時(shí)期的中國東部溫度變化.第四紀(jì)研究，25（2）：129-140.

支崇遠(yuǎn)，王開發(fā)，藍(lán)東兆，等，2003.閩南第四紀(jì)晚期沉積硅藻組合與古環(huán)境研究.微體古生物學(xué)報(bào)，20（3）：244-252.

竺可楨，1973.中國近五千年來氣候變遷的初步研究.中國科學(xué)，16：168-189.