晚第四纪全球性气候波动导致海平面升降变化，环渤海沿岸历经沧桑，形成了多套海相与陆相交替沉积的地层。广饶地区最大海侵古岸线的研究，对于了解晚更新世以来环渤海内陆地区的沉积演化，古气候、古环境变迁及海平面变化具有重要意义。

前人对渤海沿岸海侵研究中，主要集中于渤海湾西岸、现代黄河三角洲和莱州湾滨岸地带。通过野外钻孔取芯，从岩性、古地磁、孢粉、第四纪测年、古生物学等特征进行分析，研究滨岸古环境变迁，海侵期次，海侵界线等问题^{[1, 2, 3, 4, 5, 6, 7]}。赵松岭等^[4]、苍树溪等^[2]通过对渤海湾西岸钻孔的研究，阐明了晚更新世至全新世三次海侵事件发生的期次、规模、界线等；王绍鸿^[1]最早对莱州湾西岸晚第四纪海相层进行了研究，论述了海侵发育的最大界线及其沉积环境；汪品先等^[3]对我国东部海侵进行了研究，认为我国第四纪海侵从更新世到全新世逐渐加强，并具有地区性差异；韩有松等^[5]对渤海莱州湾滨海平原第四纪进行了研究，阐述了海侵发育的层位、期次及大致界线；张祖陆等^{[19, 22]}对莱州湾南岸平原黄土阜进行了研究，受海侵作用影响黄土阜呈带状沿海岸线东西向分布。目前，前人研究中很少涉及内陆最大海侵古岸线地带海侵沉积演化，对海侵与风成黄土沉积演化的时空展布，及其与古气候、古环境变迁的对应关系等也未见相关报道。

结合前人研究，借助钻井资料和样品测试资料，对广饶地区晚更新世以来的海侵沉积演化进行沉积学综合分析，尝试建立以广饶研究区为代表的渤海西南内陆最大海侵古岸线地带海侵沉积演化模式，并对晚更新世以来所对应的古气候变迁进行分析。

广饶研究区位于现代黄河三角洲平原的南部，莱州湾南岸平原的西侧，处于南部山区丘陵冲积平原向滨岸海积平原过渡的地带，区内发育小清河、淄河两条水系(图 1)。广饶地区处于暖温带半湿润季风气候区，年平均气温12.6℃，年降水量约600 mm，夏季降水约占全年的76%^[12]。历次黄河改道都未波及到研究区，沉积物源上受黄河影响较小。由引黄济青工程揭露的剖面所确定的黄河冲积物与莱州湾南岸平原沉积物的界线是在广饶县以北，小清河北岸^[25]。构造位置上，研究区位于东营凹陷南斜坡南缘地带，东靠广饶凸起，南临鲁西隆起，西、北接博兴和牛庄洼陷。其所属的构造区域延续了新生代以来的构造趋势，具有整体沉降，局部隆升的特点，第四系地层厚度260 m到400 m不等^[12]。

图 1 广饶地区地理位置与最大海侵古岸线界线(据韩有松,1994,资料修改) Fig. 1 The geographical location of the study area and the maximum transgressive paleoshoreline boundaries (after Han, 1994, modified)

图选项

晚第四纪间冰期，研究区所在的渤海西南沿岸发生过多期海侵事件，沉积了多套海相地层，并在近地表留下了古海蚀崖、侵蚀痕、沿岸堤和贝壳堤等海侵地貌遗迹^{[4, 8]}。距今最近的三期是：8~4 ka B.P.的全新世黄骅海侵，发育第Ⅰ海侵层；40~28 ka B.P.的晚更新世中晚期献县海侵，发育了第Ⅱ海侵层；102~70 ka B.P.的晚更新世早期沧州海侵，发育第Ⅲ海侵层。由于研究区构造区域和地理位置的特殊性，再加上晚第四纪全球性气候波动所导致的海平面升降变化，第四系浅层经历了复杂的沉积演化。

研究区有钻孔取芯井4口(CQ1孔、CQ2孔、CQ3孔全取芯，CQ4孔未完全取芯)和静力触探井88口，采用旋转机械钻机取芯获得岩芯，每10 cm取样测试分析。根据沉积物的粒度测试和微体古生物的鉴定，结合宏观的沉积构造等特征，判别沉积环境。沉积物粒度特征是判别沉积时自然地理环境和水动力条件的良好指标，运用青岛中科院海洋研究所Master Sizer 2000型激光粒度仪测得粒度数据，经处理得到各种粒度参数，判别沉积类型、沉积环境等。微体古生物在中国石油大学(华东)古生物实验室鉴定，微体古生物对生活环境非常敏感，不同的沉积环境具有不同的种属特征，通过分析沉积物中介形虫、有孔虫、轮藻等微体古生物组合特征，恢复海相、陆相沉积环境。同时，借用前人研究的孔2、孔1、孔10资料^[1]，进行纵向上沉积相对比分析。研究区CQ1孔、CQ2孔、CQ3孔单井相分析如下：

(1) CQ1孔单井相特征

CQ1孔(37°06′27.5″ N，118°28′56.9″ E)位于研究区北部小清河附近(图 1)，钻孔高程+9 m，钻深30 m，取芯率91%，0~2 m为近地表沉积的耕作土，未取芯。该孔整体上以海陆交互作用沉积为主，发育两期海相层，根据钻井取芯岩性综合分析，共划分为9个沉积地层单元(图 2)：

图 2 CQ1井综合柱状图 Fig. 2 The integrated histogram of CQ1

图选项

第1沉积单元，为河流相沉积。埋深2.0~4.7 m，下部发育灰黑色泥，含螺化石，为河漫湖泊沉积；上部以泥质粉砂为主，为河漫滩沉积。

第2沉积单元，为古土壤层。埋深4.7~6.0 m，黄褐色泥，含钙质结核，直径0.5 cm。

第3沉积单元，沉积第一海相层。埋深6.0~6.8 m，下部发育青灰色粉砂，含大量贝壳碎片；上部青灰色泥，贝壳碎片稀少。

第4沉积单元，河流相沉积。埋深6.8~16 m，自下向上由棕黄色细砂、粉砂过渡为黄褐色泥质粉砂、泥和灰黑色泥。由河道滞留沉积向上过渡为边滩沉积、河漫沉积，构成了典型的河流相“二元结构”。

第5沉积单元，为黄土沉积层。埋深16.0~17.0 m，岩性以黄褐色粉砂质泥为主，底部含钙质结核，直径1.0 cm。

第6沉积单元，沉积第二海相层。埋深17.0~23.6 m，以黄褐色粉砂为主，发育贝壳层。

第7沉积单元，河流相沉积。埋深23.6~26.5 m，表现为正粒序，由棕黄色粉砂向上过渡为黄褐色泥质粉砂、泥和灰黑色泥，河流相“二元结构”明显。

第8沉积单元，为风成黄土层。埋深26.5~29.0 m，以棕黄色粉砂质泥、泥为主，含钙质结核，直径1.0 cm。

第9沉积单元，河流相沉积。埋深29.0~30.0 m，底部为棕黄色粉砂，上部发育灰黑色泥沼沉积。

(2) CQ2孔单井相特征

CQ2孔(37°03′51.9″ N，118°28′32.7″ E)位于研究区中北部(图 1)，高程+12 m，钻深30 m，取芯率89%。0~3 m为近地表沉积，未取芯。根据钻井取芯岩性分析，共划分为7个沉积地层单元，仅在中下部发育第二期海相层。具体分析如下：

第1沉积单元，古土壤层。埋深3.0~5.0 m，岩性为棕黄色泥。

第2沉积单元地层，河流相沉积。埋深5.0~7.5 m，上部发育河漫滩沉积，岩性以棕黄色含泥粉砂为主；下部是湖沼沉积，发育黑色泥沼层，可见螺化石。第一海相层海退期，河流进积沉积。

第3沉积单元，古土壤层。埋深7.5~15.0 m，以黄褐色泥、粉砂质泥为主，见少数螺。

第4沉积单元，古土壤—风成黄土层。埋深15.0~19.2 m，上部岩性以棕黄色泥质粉砂为主；下部黄褐色粉砂质泥，含钙质结核。

第5沉积单元，第二海相层。埋深19.2~20.1 m，岩性以泥质粉砂为主，少量贝壳碎片。

第6沉积单元，风成黄土—古土壤层。埋深20.1~25.7 m，上部发育黄褐色粉砂质泥，含钙质结核；下部发育棕黄色泥质粉砂。

第7沉积单元，风成黄土—古土壤层。埋深25.7~30.0 m，上部发育棕黄色粉砂质泥和黄褐色粉砂质泥，含钙质结核；下部发育棕黄色泥质粉砂。

(3) CQ3孔单井相特征

CQ3孔(37°04′13.9″ N，118°25′43.3″ E)位于研究区中南部(图 1)，高程+10 m，钻深30 m，取芯率93%，地表以下0~3.0 m未取芯。根据钻井取芯岩性分析，划分为12个地层单元，以棕黄色、褐黄色的黄土沉积为主，黄土—古土壤旋回沉积为特征，未见海相沉积地层。具体划分如下：

第1沉积单元，近地表黄土层。埋深3.0~3.8 m，以棕黄色、黄褐色泥为主，质地疏松，含红棕色斑块黏土。

第2沉积单元，湖沼沉积。发育灰黑色泥沼层，埋深3.8~4.6 m，岩性灰黑色、灰棕色泥，含少量螺化石。

第3沉积单元，黄土层。埋深4.6~6.5 m，岩性以棕黄色粉砂质泥为主，致密，零星分布黄褐色的黄铁矿质点，5.0 m处发育钙质结核，呈层状分布。

第4沉积单元，泥沼沉积。埋深6.5~7.8 m，灰黑色泥为主，含螺化石，下部灰黑色与棕黄色泥混合沉积。

第5沉积单元，黄土层。埋深7.8~11.0 m，岩性以棕黄色、黄褐色的粉砂质泥、泥为主，在9.0 m和10.8 m分别发育钙质结核，似圆状，直径1~3 cm不等。

第6沉积单元，黄土层。埋深11.0~15.8 m，整体泛黄色，11.0~12.8 m，岩性较粗，疏松，以泥质粉砂、粉砂为主；12.8~15.8 m，岩性细，致密，发育泥和粉砂质泥。

第7沉积单元，古土壤层。埋深15.8~16.2 m，红棕色泥、粉砂质泥为主，发育零散灰白色斑块(未完全固结的钙质结核)。

第8沉积单元，黄土层。埋深16.2~17.5 m，棕黄粉砂质泥，致密坚硬，17.2 m发育钙质结核层。

第9沉积单元，黄土—古土壤层。埋深17.5~22.4 m，整体棕红色，岩性上下分两段：下部，20.0~22.4 m以粉砂质泥、泥为主，质地疏松，含有微小的蜗牛化石和黄褐色的铁质；上部17.5~20.0 m以粉砂、泥质粉砂为主，颜色变深，发育河道沉积的粉砂层。

第10沉积单元，黄土—古土壤层，同上单元，埋深22.4~24.0 m，为黄土—古土壤旋回型沉积层，下部是致密的黄棕色粉砂质泥的黄土沉积，含钙质结核；上部红棕色的古土壤层，岩性疏松，顶部以粉砂，泥质粉砂为主，属于河道沉积。

第11沉积单元，黄土—古土壤层。同上单元，埋深24.0~27.0 m，底部是含钙质结核的黄土层沉积，向上颜色变深，属古土壤层，顶部发育河道沉积的黄棕色粉砂层。岩性整体向上变粗。

第12沉积单元，河流相沉积。埋深27.0~30.0 m，以黄棕色粉砂和含细沙的粉砂为主，为古河道沉积。

通过对研究区钻井取芯分析，研究区北部小清河附近，地表以下30 m内自上而下发育第一、第二两期海相地层(图 2)，并向研究区南部逐渐尖灭，到南部过渡为风成黄土沉积。研究区南部的黄土沉积区大致以埋深17.0 m为界分两部分：上部以黄土层沉积为主；下部以黄土—古土壤层旋回型沉积为主，并发育多层河道粉砂。

钻孔综合分析得出，研究区北部CQ1与CQ2孔之间存在两期最大海侵古岸线。以古海岸线为界，将研究区自北向南划分为三部分：研究区北部为海陆交互作用沉积区，自上而下发育第一、第二海相层，并在中部、底部沉积黄土层；研究区中部为海向陆过渡区，仅发育第二海相层；研究区南部为黄土沉积区，海侵层完全尖灭，以风成黄土沉积为主，并发育河流相粉砂。

山东黄土沉积范围广泛，从鲁中山麓边缘、山前平原、冲积平原到滨海大陆架都有分布，按自然地理区位可大致分为两大部分，边缘海陆架黄土(滨岸平原与岛屿区)和鲁中山区黄土。山东省第四纪黄土沉积地层，主要划分为大站组和羊栏河组^[18]。大站组，发育黄土层，古土壤层不发育或少发育，在临淄厚7.7 m，青州傅家庄厚3.2 m，由研究区向南部丘陵沉积地层逐渐变薄，在地质年代与黄土高原的马兰黄土相对应，热释光测年在68 000±5 000年~12 927±11 293年之间，为晚更新世玉木冰期沉积。羊栏河组，山东最老的黄土堆积层，发育多层古河道粉砂，以黄土层与古土壤层交替旋回沉积为特征。在青州西南傅家庄黄土热释光年龄为22.7±1.8万年~79.3±6.4万年，属中更新世里斯冰期沉积，沉积年代上可与黄土高原离石黄土进行地层对比^[18]。

研究区黄土沉积物颜色以黄色为基调，主要为灰黄、棕黄和棕红色，岩性由泥质粉砂和粉砂质泥组成，伴生钙质结核，常见蜗牛化石，垂向上与古土壤交替互层，野外剖面见垂直节理(图 3)。

图 3 研究区黄土层野外特征 A．黄土层野外M剖面(图1)；B．钙质结核；C．蜗牛化石 Fig. 3 Characteristics of loess layer in field

图选项

沉积物粒度特征作为研究较为成熟的指标被广泛应用于风成沉积的研究中^[21]，不同沉积成因具有不同的粒度特征。研究区黄土沉积物粒度频率曲线(图 4a)与河流相(图 4b)存在明显的差异，与前人研究成熟的青州黄土^[30]形态接近，粒径在0~100 μm之间，曲线呈双峰形，主峰粒度区间在40~60 μm，次峰在1 μm左右，曲峰分叉且尖锐，具典型的风成沉积特征。研究区黄土沉积物(图 4a)与本文采样测得的青州云门山剖面黄土粒度特征非常相似(图 4c)，平均粒径15~22 μm，两者都比西峰黄土(平均粒径10~16 μm)^[30]粗，分选差，频率曲线形态(图 4d)也有差异，说明与黄土高原黄土有不同的物源组成。

图 4 研究区黄土粒度频率曲线 A．CQ3孔不同深度样品；B．河流相；C．青州黄土；D．西峰黄土(据彭淑贞,2007) Fig. 4 The loess grain size frequency curves

图选项

据钻孔取芯分析，将广饶地区出现的黄土(以CQ3孔为例)与前人研究的周边(临淄、青州、临朐、蓬莱)风成黄土进行了对比(图 5)。研究区南部黄土沉积区在17 m以下发育多期古河道粉砂层，岩性以粉砂、泥质粉砂为主，并发育多层偏宗红色古土壤层，¹⁴C测年为中更新世，符合羊栏河组的特征。据此，以17 m为界将研究区南部黄土层划分为两部分：5~17 m为大站组，晚更新世沉积，以粉砂质泥为主，粉砂较少，以发育黄土层为特征；17~30m为羊栏河组，未见底，中更新世沉积，岩性以泥质粉砂为主，发育多层古土壤，伴有河道粉砂沉积。

图 5 研究区与周边地区黄土地层对比图(据王世进,1999资料修改) 1.坡积层及湖沼层；2.大站组(黄土层)；3.羊栏河组(古土壤及黄土层)；4.粉砂、砂、砾石层；5.玄武岩；6.基岩；7.热释光法取样点及年龄/a Fig. 5 The loess stratigraphic correlation and classification(after Wang,1999, modified)

图选项

地质历史时期黄土沉积范围比现今所认识的广大的多。不同学者通过对环渤海滨浅海的钻孔研究发现，不同深度也相继沉积含钙质结核的硬黏土质黄土层^{[7, 16]}。同时，通过微体古生物鉴定，在距海侵古岸线较远的CQ3孔，埋深20.5 m的黄土层发现了半咸水的海相微体古生物Zlyocypris sp.，表明冰期黄土沉积期有海相沉积物的混入。在晚第四纪冰期时期，气候寒冷干燥，末次冰期最盛时，中国东部海平面最大下降到-150~-160 m，海岸线到达朝鲜半岛的东侧，整个渤、黄海陆架和东海的部分陆架出露^[13]，渤海西、南内陆盛行西北风和偏北风，在冰期强劲的风力作用下，出露的渤海陆架松散沉积物，向陆一侧搬运堆积，并且强劲的北向风所携带的西北高原粉尘，在渤海西、南岸适宜地貌区沉降，经后期改造演化，形成现今的黄土堆积体^{[19, 30]}。

渤海西南沿岸海侵的研究，不同学者已在海侵层沉积年代，发育层位，最大古海岸线位置等方面取得了广泛性认识^{[3, 6, 10]}。王绍鸿等^[1]对莱州湾南岸平原西侧海侵的研究认为，晚更新世以来发生的两次海侵(献县海侵、黄骅海侵)在最大海侵面时曾分别到达广饶、博兴地区。祁雅丽等^[17]运用地面雷达测试技术，对渤海西南内陆全新世大暖期最大海侵线进行了精确重建，得出的全新世最大海侵古岸线经过广饶研究区北部，与研究区钻井取芯分析所得到的古岸线位置一致(图 1)。

研究区北部小清河附近，自上而下发育两期海相层，与黄土层存在鲜明的沉积特征差异。第一海相层，埋深6.0~6.8 m，岩性以青灰色粉砂、泥为主，含大量贝壳碎片；第二海相层，埋深18.0~23.6 m，岩性以青灰色、深褐色粉砂为主，贝壳完整。王绍鸿等^[1]借助1~10孔对莱州湾南岸西侧海相层进行研究，确定了海相层发育的期次、年代、埋深^[1]，之后不同的学者进行了证实^{[5, 17]}。其中孔2与研究区CQ1孔距离非常接近(图 1)，在统一高程下，两孔海相层埋深相同(图 6)，在沉积年代上可以认为相同时期的海相地层。以CQ1孔海相层作为研究区的对比标志层，与周边前人研究的黄金寨9孔、乔家庄10孔和牛庄、王家岗地区(图 1)的海相地层进行横向对比具有可行性。

图 6 研究区与周边地区海侵地层对比图(据韩有松,1994; 王绍鸿,1979资料修改) Fig. 6 The transgression stratigraphic correlation and classification (modified after Han, 1994; Wang, 1979)

图选项

从对比中可以看到，研究区第一海相层对应前人研究得出的10~4 ka B.P.的全新世(第Ⅰ)海侵层，为黄骅海侵沉积；第二海相层对应40~28 ka B.P.的晚更新世晚期(第Ⅱ)海侵层，为献县海侵沉积。两期海侵层由研究区向现代海岸逐渐增厚，第二海侵层比第一海侵层沉积厚，规模大，经历的演化时间久(图 6)。

晚更新世以来环渤海地区发生的不同期次海侵，在沉积环境、海平面变化、古气候变迁等方面存在明显的差异^{[1, 3, 13, 14]}。

(1) 第一海侵层特征

广饶地区发育的第一海侵层对应渤海沿岸的第Ⅰ海侵层，为全新世黄骅海侵沉积。在渤海湾西南沿岸第Ⅰ海相层分布比较广泛，根据研究区的钻孔分析，并结合前人研究^{[1, 5, 23]}，最大古海岸线大致分布于乔庄、博兴、草桥、牛头镇、南河、侯镇、永安、李家埠、夏店及新河一带(图 1)。其海侵范围比晚更新世早期的沧州海侵要大，比晚更新世中期的献县海侵要小，在莱州湾南岸平原东部地区，全新世的海侵范围比晚更新世中期的大。

根据最大海侵古岸线海侵层位埋深，暂不考虑构造沉降和沉积压实等引起的地层沉降，取芯资料显示，研究区第一海侵层最大海侵面的高程为2.0~2.5 m，这与众多学者研究所认为的全新世海侵存在2.0~3.0 m的高海平面的结论相吻合^[11]。可以看出，全新世海侵以来，区域构造运动趋于稳定，没有发生明显的地层沉降。

通过研究区钻孔取样微体古生物的鉴定，结合1~10钻孔^[1]分析(图 1)，全新世黄骅海侵沉积环境主要以水深小于15 m的潮间带—浅海为主。其中第一海侵层有孔虫化石群较多的五种是：嗜温转轮虫(Ammonia tepida Cushman)、连接转轮虫(Ammonia annectens Parker and Jones)、显脐转轮虫(Ammonia flevensis Hofker)、圆形短五块虫(Quinqueloculina akneriana rotunda Gerke)、易变筛九字虫(Cribrononion incertum Williamson)，这些有孔虫在现代渤海海洋环境中主要生活在水深小于20 m的水域。而且，第一海侵层中45%的软体动物化石群以潮间带种为主，浅海种最少，只占10%^{[1, 15, 26]}。

(2) 第二海侵层特征

研究区第二海侵层与环渤海沿岸献县(第Ⅱ)海侵层相对应，为距今4.2~2.8万年间晚更新世中期海侵沉积。该海相层分布范围广泛，经钻孔证实，在渤海西南沿岸古海岸线大致从乔庄西侧，向南经博兴，过研究区广饶县以北，向东一直延伸到寿光寒桥和昌邑，至新河、灰埠北和虎头崖附近^{[1, 5]}(图 1)。

根据海侵层在研究区最大海侵古岸线位置埋深，第二海侵层最大海侵面的现代标高是-14.6~-8.0 m。事实上，晚更新世以来的区域构造沉降对海侵层的埋深具有直接的影响。魏灵等利用渤海沿岸众多钻孔微体生物化石的组合特征来恢复海水深度，重建40~28 ka B.P.献县海侵最大海侵面的现代标高是-11.0~-5.0 m^[11]，与研究区第二海侵层的最大海平面现代标高相差3.6~3.0 m。可见，自晚更新世中期第二海侵层沉积以来，受构造沉降和沉积压实作用，地层沉降了约3.0 m。

在古生物学特征上，研究区第二海侵层有孔虫化石群较多的五种是：嗜温转轮虫、易变筛九字虫、缝裂企虫(Elphidium mangellanisum Heron-Allen and Earland)、整洁五块虫(Quinqueloculina bellatula Bandy)、异地企虫(Elphidium advenum Cushman)，这些有孔虫主要生活在水深15~25 m的浅海区。在软体动物化石群上，第二海侵层与第一层海侵层正相反，浅海物种最多，占总种数的34%，潮间带物种最少，占23%^{[1, 15, 26]}。可见，第二海侵层以水深10~20 m的滨岸浅海环境为主，比第一海侵层沉积水体要深，海侵幅度更大，演化持续时间更久。

研究区地表以下30 m内，存在两期最大海侵的古岸线，以古岸线为界，南部以风成黄土沉积为主，北部发育两期海相层(黄骅海侵、献县海侵)。晚第四纪冰期时期，偏北风所携带的西北内陆干旱地区和裸露的边缘海陆架的沙尘，在山地丘陵，山前平原及渤海西南沿岸适宜位置沉降，形成大面积黄土堆积带。间冰期气温回升，海平面上升，发生海侵事件，侵蚀改造前冰期黄土。

晚更新世发育的36 ka B.P.献县海侵，在最大海侵面时到达研究区，侵蚀、超覆于中更新世羊栏河组黄土层之上，沉积第二海侵层。之后，进入晚更新世末次冰期最盛时期(玉木Ⅲ冰期)，沉积大站组黄土。伴随最后一次冰期结束，转入全新世大暖期，发生了黄骅海侵，侵蚀大站组黄土，并沉积第一海侵层，再次经历了黄土—海侵沉积旋回。

在垂向上，两期海侵层与黄土—古土壤层呈叠置互层关系，海侵层穿插超覆于南部的黄土层中，与黄土层呈不等时沉积接触(图 7)。在莱州湾南岸平原黄土层被海侵侵蚀，形成海蚀黄土崖，称为“黄土阜”，呈带状沿海岸线分布^{[19, 22]}。

图 7 研究区海侵、黄土沉积演化模式 Fig. 7 The deposition evolution model of transgressions and loess

图选项

海侵的侵蚀改造作用对于山东省黄土分布格局和沉积演化具有重要影响。间冰期海侵侵蚀冰期黄土，在最大海侵位置沉积了海侵黄土的古岸线，使山东省近地表黄土分布的北缘与全新世最大海侵古岸线的界线相一致。

笔者把这种海侵沉积地貌称为“上超型海蚀黄土古岸线”。其含义是，间冰期气候回暖，海平面上升，渤海陆架及内陆平原黄土沉积区发生海侵剥蚀作用，造成海侵侵蚀上超于先前冰期沉积的较老黄土层之上，在最大海侵位置形成了海侵黄土古岸线的演化模式(图 7)。

研究区海侵、黄土的沉积演化与晚第四纪以来的古气候变化、古环境变迁息息相关。晚第四纪全球气候变化频繁，出现多次冰期与间冰期的气候冷暖旋回，这种气候的旋回性变化在我国晚第四纪古里雅冰芯氧同位素记录、渤海钻孔古地磁和海侵研究中都有相应体现^{[13, 14, 20, 28]}。苍树溪等^[2]对渤海湾西岸中更新世以来古生态、古气候、古地理的研究认为，布伦赫斯正极性世(距今大约69万年)以来，发生的五次短期地磁极倒转事件与大约以10万年为周期的古气候摆动，微体古生物组合更替等相对应。

古气候波动的周期性在末次冰期和全新世存在千年和百年尺度的变化周期^{[13, 14]}，并具有“旋回加旋回”的特点^{[24, 27, 29]}。海侵演化过程中的表现是，在长周期的海平面升降变化中，叠加了中短期的海平面波动，即长周期的海退存在短周期的海侵。学者对渤海湾海侵的研究已证实，在全新世大暖期海侵事件的海退期，有短时期的海侵发生^{[13, 14]}。

晚第四纪气候演变是个非常复杂的问题，是多种因素综合作用的结果。根据研究区海侵与黄土沉积演化关系(图 2)，对应不同的沉积演化特征，结合渤海沿岸古气候变化特征^{[1, 6, 14, 27]}，将研究区晚更新世以来的气候演变划分为6个时期，具体划分如下：

(1) 干冷期，对应玉木Ⅱ冰期，海平面下降，侵蚀基准面下降，发育深切河流，在干冷季风下沉积羊栏河组黄土，埋深-17 m以下。

(2) 温暖湿润期，为玉木冰期上亚间冰期，气候回暖，海平面上升，发育献县海侵，在最大海侵面时到达研究区，侵蚀改造前期黄土；同时，由于侵蚀基准面的上升，河流动力较弱，流水不畅，伴随湖沼相沉积，埋深-17~-9 m。

(3) 干冷偏凉期，对应玉木Ⅲ冰期，气候转冷，但平均气温比玉木Ⅱ冰期略高些，海平面最大下降到-100 m以下，渤海、黄海陆架裸露，在风力作用下沉积大站组黄土；河流相在研究区发育较弱，埋深-9~+1 m。

(4) 温暖期，进入全新世，气候转暖，发育黄骅海侵，侵蚀大站组黄土；由于流路不畅，河流动力较弱，湖沼相发育广泛，埋深+1~+2 m。

(5) 潮湿湿润期，海侵后退期，伴随陆上水位基准线下降，在先期海侵的低洼地带淤积成湖，发育大面积湖沼沉积，覆盖在研究区南部大站组黄土和北部第一海侵层之上，¹⁴C测年3 000 a B.P.左右，埋深+2~+3 m。

(6) 现代期，与近代气候相近，受人类活动影响较大，一般地表以下2 m或3 m。

(1) 研究区地表以下30 m内，南部沉积风成黄土，北部发育两期海相层(10~4 ka B.P.的黄骅海侵层和40~28 ka B.P.的献县海侵层)，侵蚀上超于南部的黄土层中，在渤海西南岸内陆平原形成“上超型海蚀黄土”的演化模式。

(2) 对研究区钻孔取芯研究证实，40~28 ka B.P.的献县海侵最大海侵界线并没有越过广饶县城，而是位于广饶县城以北，草桥、小清河以南，对此进行了确定(图 1)。

(3) 研究区第一海侵层存在海拔2~3 m的高海侵面，但第二海侵层以水深10~20 m的滨岸浅海环境为主，比以潮间带—浅海沉积环境为主的第一海侵层沉积水体要深，海水陆泛的幅度更大，持续时间更久。

(4) 根据研究区海侵与黄土的沉积演化特征，将研究区晚更新世以来的气候演化划分为6个阶段：干冷期—温暖湿润期—干冷偏凉期—温暖期—潮湿湿润期—现代期，分别对应海侵与黄土演化的不同气候类型。

但是，关于两次海侵的争论远没有停止。前人对环渤海沿岸海侵研究也发现^{[2, 11]}，40~28 ka B.P.的献县海侵比10~4 ka B.P.的黄骅海侵，海水陆泛的幅度更大，海水更深，与本研究得出的结论一致。但与全球气候处于末次冰期的间冰阶，冰川部分消融的观点不切合。研究区海侵与黄土的时空沉积演化关系仅是环渤海地区古气候演变的物质载体，而且影响两期海侵古气候变化差异性的原因是多因素的，今后还应从构造沉降差异，区域气候波动差异，古地貌，古地理和古环境等进行更加深入的多角度研究。

致谢 对参加野外工作的张光德、徐钰、李帅、晋同杰等同志，和样品测试中给予热忱帮助的青岛中科院海洋研究所，在此一并致谢。