珠江三角洲西缘晚第四纪沉积演化和最大海侵古岸线的重建

周良，王洋，杜学斌，卜建军，吴俊，吕万军，刘秀娟

1.中国地质调查局海口海洋地质调查中心，海口 570100

2.海洋地质资源湖北省重点实验室，中国地质大学，武汉 430074

3.中国地质调查局武汉地质调查中心，武汉 430200

0 引言

在21 世纪全球气候变暖的背景下，气温升高导致两极冰盖消融引起全球性海平面上升，对沿海地区造成一系列的影响，如沿海低地被淹没、海岸受侵蚀、海水入侵等地质环境问题，越来越引起人类社会的广泛关注[1]。晚第四纪全球性气候波动导致海平面升降变化，沿海地区形成了多套海相与陆相交替沉积的地层[2]。对沿海地区最大海侵古岸线的研究，不仅有助于我们了解沿海地区晚更新世以来的古气候、古环境演变，而且可以对未来发生海侵的影响提供对策和帮助。

前人对广东沿岸海侵的研究，主要集中在珠江三角洲和粤东韩江三角洲等区域[3-5]，对珠江三角洲以西到雷州半岛之间广阔的粤西地区研究较少。20世纪70年代前，由于缺少系统的年代学研究，认为珠江三角洲地区第四系只发育全新统[6]。70年代开始，随着14C、热释光等年代测试技术的发展，通过野外钻孔取心，从岩性、古地磁、孢粉、第四纪测年、古生物学等特征进行分析，确立了珠江三角洲开始的沉积年代为MIS3 期，晚更新世以来主要有2 套较大的沉积旋回，相对应2 次规模较大的海侵（MIS3 和全新世）和1 次大的海退[7-10]。珠江三角洲晚更新世最大海侵范围达到了太平—中堂—顺德一线以南[11]，全新世最大海侵范围更大，大致在黄埔—石楼—大良一线附近[12]。前人在韩江三角洲的研究证实了该地区第四系沉积始于晚更新世中期，韩江三角洲在晚更新世中期整体沉降，这期海侵在20 ka B.P.左右达到最大范围，最西到达贾里附近[13]，仅影响到三角洲外侧平原；全新世后沉积环境趋于稳定，冰后期以来海平面大幅上升，在7 800 a B.P.海侵达到最大范围，最大海侵古岸线东北段达到潮州附近[13-14]。

第四纪海侵造成中国沿海大部分低洼平原地区被淹没，最大海侵线甚至向陆延伸达300 km[15]。早期学者在开展调查时，主要以14C测年等测年方法确定地层年代，并根据生物组合和沉积特征来推测古海岸线位置[16]。或者利用贝壳堤、古海蚀崖等典型的古岸线地貌和沉积标志，来恢复古海平面等[17-18]。由于这些古岸线标志物不易完整的保存下来，导致恢复的古海岸线位置存在不确定性。因此前者仍是恢复最大海侵古岸线较为有效的方法，这种方法通过海相地层的分布来确定最大海侵古岸线，第四纪沿海地区发生周期性的海进海退，发育了海、陆交互的沉积地层，海相地层向陆方向延伸并减薄至尖灭，尖灭的位置即最大海侵古岸线的大致位置（不考虑沉积物压实和后期改造），而探测海相地层的分布主要有探地雷达[19]和钻孔两种手段[20]。

总得来说，前人的研究主要集中在晚第四纪海侵地层的时代归属、海侵的最大波及范围等问题，对最大海侵范围与构造升降、海平面变化之间的关系等问题的认识仍然不够深入，然而这些都是今后评价海侵对沿海地区影响的重要理论依据，值得我们进一步研究。本文结合前人研究，借助钻井资料和样品测试资料，对珠江三角洲西缘台山地区晚更新世以来的沉积记录进行综合分析，研究台山沿海地区晚第四纪沉积演化特征，并尝试重建最大海侵古岸线的位置，探讨其与新构造运动和海平面变化的关系，为研究本区新构造运动速率提供新证据、新思路。

1 区域地质概况

研究区位于广东省中南部、南海北岸、珠江三角洲西缘（图1a），东北侧百公里之外有珠江三角洲的崖门、鸡啼门、虎跳门及磨刀门四个河流入海口，西北侧有镇江河的入海口，广海湾的北部还有一个小型河流三夹海河入海。台山沿海一带的地貌主要是丘陵、低山，海岸弯曲，岛屿众多，具港湾海岸地貌特征。研究区中南部地形较平坦，发育开阔的三角洲冲积平原和海积平原；东部和南部多丘陵，地形起伏较大，地表为残坡积土体或不同程度的风化花岗岩。区内低山、丘陵、平原交错，海岸发育沙洲，地貌单元丰富（图1b）。

研究区划属华南地层区东江分区和沿海分区，地层发育，分布范围颇广。主要出露侏罗纪和第四纪地层[22]。第四纪地层以河口三角洲沉积为主，主要由晚更新世的礼乐组和全新世的桂洲组组成。礼乐组不整合于基岩风化壳之上，自下而上划分为石排组、西南镇组和三角组，时代属晚更新世；桂洲组指平行不整合覆于礼乐组或超覆于基岩风化壳之上的一套地层，自下而上可划分为杏坛组/三角组、横栏组、东升组、万顷沙组和灯笼沙组，时代属全新世。

研究区构造上属于华南褶皱系粤中凹陷南部的增城—台山隆断束（图1a），区内断裂构造较发育，按断裂构造产出形态主要发育NNE 向断裂，控制着河道延伸方向、古海岸线和第四系沉积物的展布。喜马拉雅运动期研究区第三纪一直处于隆升状态，一直到第四纪晚更新世才接受沉积。在地壳的继承性和差异性的新构造升降运动中，大多数学者认为研究区整体上主要以继承性构造抬升为主[21]。

2 材料与方法

项目于2017年1月到2018年7月在台山市陆域以及广海湾、川山群岛海域共采集30 个钻孔岩心。岩心经现场剖开，进行详细的岩性描述、分层、拍照。

选取ZK06、ZK11、ZK22、ZK29、ZK45 孔的10 个层位的样品和DZK02孔的3个层位样品分别送往北京大学第四纪年代测试实验室和美国Beta分析实验室进行14C测年，其中贝壳样品4件，含碳屑泥质沉积物样品4件，腐木和植物碎屑样品5件。测试结果通过CALIB 程序对14C 的放射性年龄数据校正为日历年龄，在校正曲线的选择上，贝壳等海洋无机碳样品选用的是Marine13 校正曲线，对于淤泥样品中的有机质，由于样品来自海陆交互的海岸带地区，校正曲线选择海洋和北半球大气混合曲线。根据Southonet al.[23]对中国南海地区区域海洋碳储库效应的研究，我们在校正过程中选用ΔR=-25±20 a 进行校正。采集ZK22、ZK25、ZK34 孔的8 份样品送往南京师范大学光释光实验室进行石英颗粒（4～11 μm）的光释光测年。采用XL3t950 型手持式XRF 矿石元素分析仪对DZK02 钻孔的岩心样品进行元素测试，采样间隔约10 cm，共测试样品340 个。使用SM-30 型便携式磁化率仪，以20 cm 为间隔对ZK06、ZK20、ZK33、DZK02 孔进行磁化率测试，共测试样品450 个。对DZK02 孔23 个沉积物样品进行微体化石（有孔虫）的处理、化石挑选后送至中国地质大学（武汉）生物地质与环境地质国家重点实验室进行鉴定。

3 结果

3.1 测年结果及岩性地层划分

14C和光释光测年结果显示年代随地层埋深的增大而逐渐变老（表1，2）。顶部人工填土沉积物年代基本限定在1 ka B.P.之内，其下部淤泥层沉积物年代在11～1 ka B.P.，大约在16～11 ka B.P.（末次冰期）期间形成了花斑黏土层，花斑黏土层下部的粉砂/黏土年代大约在38～16 ka B.P.，底部最老地层年代大致在50 ka B.P.（图2）。

图2 研究区钻孔N-S 和W-E 方向地层对比横剖面图Fig.2 Comparison of strata in N-S and W-E directions determined from nearby cores

表1 14C测年结果Table 1 AMS 14C dates

研究区钻孔整体上以海陆交互作用沉积为主，发育两期海相层，根据钻井取心岩性综合分析，研究区从基底往上，共划分为4 个沉积地层单元（表3），每个沉积单元岩性组合和识别特征明显。

表3 研究区岩性地层划分Table 3 Lithological stratigraphic division of study area

3.2 地球化学特征

蓝先洪等[24]分析珠江三角洲第四系和现代河流与浅海的表层样品地球化学数据，使用Sr和Ba的含量和比值判别海相与陆相等沉积特征，将Sr>100、Ba>400（单位为10-6，下同）用来指示珠江三角洲海相沉积物，Sr<50、Ba<300指示珠江三角洲陆相沉积物；（Sr/Ba）<0.5 时代表微咸水沉积环境，Sr/Ba 值在0.5～1.0时代表半咸水沉积环境，（Sr/Ba）>1.0时代表海水沉积环境。研究区利用元素含量和比值判断海陆相沉积时参考珠江三角洲的判别标准。Ba/Al 和Zn/Al代表的古生产力指标主要与东亚冬季季风和夏季季风的强化有关。

基岩和风化层以上至深度30 m，Ba<400，Sr<50，Sr/Ba 值均小于0.5，Ba/Al 值基本稳定在0.04～0.05，（Zn/Al）值基本稳定在0.005 左右；深度在21～30 m，Ba>400，Sr>50，仍处于微咸水环境，Ba/Al 值和Zn/Al值整体上出现了较为明显的升高，在26 m 处达到了一个峰值；在18～21 m 的花斑黏土层，各项指标在这个阶段都出现了明显的降低，Ba<300，Sr<50，（Sr/Ba）<0.5，（Ba/Al）<0.05，（Zn/Al）<0.01；之后18 m以上，整体上Ba、Sr 含量分别大于400、100，Sr/Ba 值波动上升，在10 m 处达到峰值，处于半咸水环境，之后逐渐减小，恢复到微咸水环境，Ba/Al 和Zn/Al 值开始升高，然后转为降低并逐步趋于稳定（图3）。

图3 DZK02 孔指相元素浓度分布图Fig.3 Concentrations of indicator elements in core DZK02

表2 光释光测年结果Table 2 OSL dates

3.3 有孔虫

经镜下鉴定统计，本批次样品中均见有数量不等的有孔虫化石，除个别样品中数量较丰富外，大多数样品中的有孔虫丰度及分异度均较低或极低。

钻孔下部深度30 m以下18～21 m的样品产出有孔虫数量较少，丰度和分异度均较低，近乎为0；深度28 m 以上有孔虫丰度和分异度升高，几乎全为底栖有孔虫，主要为个体较小的Ammonia aomoriensis,Ammonia beccarii vars.等弱海相性代表属种，偶见有Rotalidum annectens, Elphidium advenum, Elphidium hispidium等破碎壳体，应为搬运再沉积壳体。至花斑黏土层有孔虫丰度和分异度又逐渐降低，接近于0。总体来说，上部0～18 m有孔虫丰度和分异度较高，深度18～9 m丰度和分异度出现大幅度增大，后又迅速下降；深度6 m以上有孔虫分异度仍处于增大趋势。该段有孔虫以底栖类型为主，仅个别层位见有少量极细小的浮游类型壳体。底栖类型数量最多的为钙质透明壳的Elphidium、Rotalidum、Ammonia等属，瓷质壳类型的Quinqueloculina、Massilina、Spiroloculina等属亦有少量属种出现。优势属种组合为Elphidium hispidulum、Elphidium advenum、Rotalidium annectens、Ammonia beccarii等，优势属种在深度6～18 m 和21～30 m两部分较富集，该类型组合多为广盐性分子，主要出现于近岸的温暖浅水环境中，本化石群中偶见有一些海陆交互相环境中的典型分子（如Ammonia beccariivars.、Cribrononion subincertum等），同时又极少或未见浮游类型（图4）。

图4 DZK02 孔有孔虫优势种属分布及丰度、分异度变化特征Fig.4 Relative abundance of main species and total abundance and diversity of planktonic foraminiferal fauna in core DZK02

3.4 磁化率

气候温暖阶段，大洋环流增强、淡水注入增多会带入较多的外源碎屑矿物，在较强氧化条件下，磁性矿物浓度增加，磁化率较高一般指示较高的沉积水位[25]。钻孔下部磁化率普遍比较稳定，仅有部分轻微的波动，往上粉砂层磁化率有了一定幅度的波动上升，到花斑黏土层底部磁化率较低且稳定，斑黏土层上部至上部淤泥层磁化率开始显著升高，磁化率较离散，波动幅度较大（图5）。

图5 DZK02、ZK20、ZK33、ZK06 孔磁化率变化曲线Fig.5 Changes of magnetic susceptibility in cores DZK02, ZK20, ZK3 and ZK06

3.5 最大海侵古岸线的划定

若只考虑是否发育同一期海侵地层，那么晚更新世最大海侵时，最大海侵线达到了ZK46、ZK45、ZK06、ZK42、ZK13、ZK33、ZK11 和ZK01、ZK44、ZK43、ZK31、ZK32、ZK10 以及ZK27、ZK35、ZK37 和ZK26、ZK36、ZK38 之间；全新世最大海侵时，最大海侵线大致在ZK01、ZK44、ZK43、ZK13、ZK33、ZK10和ZK02、ZK23、ZK28、ZK31、ZK32、ZK09 以及ZK27、ZK36、ZK38 和ZK26、ZK40 之间（图6）。然而，由于研究区丘陵分布广泛，我们不难发现ZK34、ZK35、ZK37三个钻孔南北两侧均有丘陵地貌分布（图1b），海侵时海水会绕过海拔较高的山地丘陵地区而首先淹没海拔较低的平原区，结合地形特征我们推断ZK34、ZK35、ZK37 三个钻孔两侧应该存在2 条晚更新世最大海侵古岸线才是合理的，除上述外在ZK34、ZK35、ZK37以南还应该存在一条（图6）。同样，该区域全新世最大海侵古岸线理论上也应该存在两条，但是由于ZK34、ZK35、ZK37 三个钻孔以南缺少钻孔资料，我们只能确定ZK34、ZK35、ZK37 三个钻孔以南存在另一条全新世最大海侵古岸线，但无法确定其准确位置。因此，我们根据海侵地层尖灭点圈定最大海侵古岸线时，一定不能忽视地形地貌对最大海侵古岸线分布的影响。

晚更新世海侵最大海侵界线向陆最远位置在ZK46和ZK01之间，距现今岸线约15.4 km，全新世海侵最大海侵界线向陆最远位置在ZK01 和ZK02 之间，距现今岸线约16.3 km，两条最大海侵古岸线走势大致相同（图6），全新世最大海侵范围略大于晚更新世最大海侵范围，总得来说极为接近。

图6 研究区晚第四纪两次最大海侵古岸线位置及3 m 等高线分布Fig.6 Location of maximum transgression coastline in Late Quaternary and 3 m contour

4 讨论

4.1 研究区沉积演化特征

根据钻孔岩性特征，并结合测年数据、微体古生物分布、地球化学和磁化率等数据，研究区晚更新世以来的沉积演化主要分为5个阶段。

第1阶段（石排组）：50～38 ka B.P.

在此之前，研究区与珠江三角洲地区一样，一直处于风化、剥蚀、侵蚀的环境，基岩裸露，直到这个阶

段才开始接受第四纪沉积[6]。该时段沉积了较厚的冲积砂砾层，沉积地层中未发现各类微体古生物记录、贝壳碎片和腐木碎屑等，Ba、Sr、Sr/Ba值均处于低值，磁化率变化幅度小，代表了古河流的河流冲积作用。

第2阶段（西南镇组）：38～16 ka B.P.

晚更新世玉木亚间冰期，随着气候转暖，海平面波动上升（图7）,研究区发生了晚更新世以来第一次海侵，发育了以粉砂质黏土至淤泥质细砂为主的沉积地层，代表海相性种属的底栖有孔虫Elphidium hispidulum、Elphidium advenum较富集，Ba、Sr 元素值指示微咸水环境，磁化率开始出现突变，代表了河流与海水潮流动力相互作用的浅海相沉积环境，约26～24 ka B.P.到达本次海侵高海平面期，在此阶段中国东海、珠江三角洲、红海等地区基本都发育了此次海侵地层，且最高海平面期的年代也基本一致[26-28]，说明此次海侵具有一定的普遍性。

虽然众多学者关于广东沿海晚更新世海侵层的年代现在仍无定论，比如Zonget al.[29]认为此次海侵层经历了之后末次冰期时的暴露环境，地层遭受风化，导致测年数据比地层真实年龄普遍偏年轻，但仍存在较大争议，但研究区和珠江三角洲众多钻孔的测年数据[7-10]仍显示晚更新世这次海侵年代上属于MIS3期。

第3阶段（三角组）：16～11 ka B.P.

约22 ka B.P.以后，全球进入了末次冰期，全球海平面发生大幅下降（图7）[30]，研究区甚至南海北部都处于陆相暴露侵蚀的环境，该时期广泛发育了由风化作用形成的杂色花斑黏土，可见少量零星底栖有孔虫，壳体呈黄褐色，遭受过风化作用，Ba、Sr元素含量重回低值，磁化率降低至稳定，指示较强的氧化环境。前人研究认为，该时期全球最低海平面大致在16～18 ka B.P.[31]，与研究区末次冰期最低海平面时间基本相同。

第4阶段（横栏组）：11～1 ka B.P.

进入全新世，全球气候变暖，海平面开始迅速回升，研究区开始发生晚更新世以来第二次海侵，发育深色淤泥和粉细砂，钻孔中出现大量微体古动物群记录，优势属种组合为Elphidium hispidulum、

Elphidium advenum、Rotalidium annectens、Ammonia beccarii等海相属性底栖有孔虫，Ba、Sr、Sr/Ba 值指示半咸水环境，磁化率再次突变至高值，代表浅海潮坪沉积体系。至6～5 ka B.P.，海平面上升至最高值（图7），研究区除丘陵山地外，大部分平原地区均被海水淹没，形成了古河口湾。这次海侵发生于全新世大暖期，持续时间较长，范围也比晚更新世海侵更广，与珠江三角洲地区研究结果相一致[9,13-14]。

第5阶段（现代沉积层）：1 ka B.P.至今

4 000 年以来，全球海平面基本与现今海平面一致（图7），并伴随着小规模的波动变化，研究区现今海陆格局基本形成，开始接受现代沉积，该阶段的沉积环境会受到一定程度人类活动的影响。

图7 50 ka B.P.以来南海及其邻近海区海平面变化曲线与研究区沉积序列对应关系（修改自张虎南等[30]）Fig.7 Sea-level change in South China Sea and adjacent seas since 50 ka B.P., and relationship to sedimentary sequence(modified from Zhang et al.[30])

4.2 研究区晚第四纪最大海侵古岸线分布与构造运动、海平面变化的关系

第四纪以来，以冰期—间冰期交替为主要特征的全球气候变化造成海平面频繁波动[32]。根据古气候研究结果，晚更新世中晚期最高海平面大约低于现今海平面数20～60 m，全新世最高海平面高于现今海平面约2～5 m，两个海侵时期的最高海平面差距数十米[33-34]，两次最大海侵范围的极度接近与海侵时期实际全球海平面高度之间的显著差异不符，造成这种现象的主要原因可能是在全球海平面变化的大背景下，研究区构造抬升作用引起的可容空间逐渐减小造成的。研究区晚更新世以来的海进—海退，主要受两极冰盖体积控制的全球绝对(水动型)海面变化和该区构造垂直升降共同控制的，“水动型”和“地动型”海面变化相互叠加[35]，还构造运动的本来面目，是深入了解研究区晚更新世以来沉积演化特征的基础和关键。

目前众多学者对华南沿海晚更新世以来的海平面变化仍存在较多争议，其中主要分为两种观点。一种观点认为晚更新世到末次冰期华南海平面从-2 m或-4 m 降至-25 m 以下，进入早全新世海平面开始大幅度回升，约在距今4 ka 时接近现今海平面高度[18,36]。而另一种观点则认为40～10 ka B.P.以来海平面在-2 m至-12 m波动，并不存在低于-20 m的低海面时期[37]，全新世开始大幅度上升，距今6 ka 左右接近现代海面，距今4 ka 期间曾出现高于现今海平面2～3 m的高海面[38]。张虎南等[30]根据华南236个样品的年代、埋深和构造升降等数据，绘制出符合华南地区的晚第四纪海平面变化曲线，本文在计算和分析时主要参考此古海面资料。

构造升降速率的计算方法较多，前人在研究中通过海南、台湾海域的珊瑚礁等古标志物的分布高度，得到了全新世以来的构造抬升速率[18]，为研究新构造运动提供了一种新思路。本文根据前人的研究思路，引入了构造升降速率的计算公式：

其中：V代表垂直构造运动速率（+代表构造抬升，-代表构造沉降），A代表海相沉积物的海拔高度，E代表海相沉积物沉积时的古海面高度（高于现今海平面为负值）[37-38]。

根据研究区晚第四纪两次海侵层位的现代埋深，暂不考虑构造沉降和沉积物压实等引起的地层升降，从最大海侵古岸线向海一侧的钻井资料可以看出，研究区晚更新世海侵地层埋深的现代标高在-3.9～-12.7 m（图8a），这与众多学者研究所认为的晚更新世中晚期存在-20～-60 m 的最高海平面的结论不符合[18,36,38]，可以得到的是，晚更新世海侵以来，研究区域构造运动整体上应该发生了明显的地层抬升，ZK22和ZK35孔内深度23.5 m和11.9 m处的样品测年结果分别为29 704 a B.P.和28 148 a B.P.，据此推算的抬升速率分别为+0.32 mm/a和+0.70 mm/a，平均抬升速率为+0.51 mm/a（表4）。全新世海侵层的现今海拔为+1.1～-15.4 m（图8b），一般认为全新世最高海面高于现今海平面5 m 左右[30]，黄玉昆等[39]根据调查资料，认为高出现代海面的沉积物大部分位于上升地区，而在下降区（韩江三角洲、珠江三角洲、漠阳江河口等）海侵层都分布在现代平均潮面以下，ZK06和ZK11 孔内深度4.6 m 和5.5 m 处的样品测年结果分别为4 859 a B.P.和4 922 a B.P.，据此推算的沉降速率分别为-0.51 mm/a 和-0.71 mm/a，平均值为-0.61 mm/a（表4）。

表4 两期海侵沉积物样品构造升降速率计算结果Table 4 Rate of vertical tectonic motion calculated from sediment transgression

图8 研究区晚第四纪两期海侵地层的现代标高（a）全新世；（b）晚更新世Fig.8 Modern elevation of marine facies stratum(a)Holocene;(b)Late Pleistocene

本文计算结果没有考虑沉积物压实、静力均衡沉降以及后期侵蚀等非构造因素的影响，加上目前学术界比较认同不存在适应于全球性的海平面变化曲线[38]，在进行不同研究时选用的海平面曲线有所差异，而且在对晚更新世海相沉积物的计算中没有按照不同的年代数据细分为晚更新世以来和全新世以来两大时段，其计算所得抬升速率表示晚更新世以来升降相抵后的整体平均速率[36]，导致不同研究计算所得构造升降速率可能会存在一定的偏差，和前人研究结果相比，虽有一定偏差，倒也合理（表5）。

表5 台山地区新构造运动研究成果Table 5 Neotectonic movement in Taishan area, as determined in this study

掌握了研究区晚第四纪构造升降速率的基本信息，结合两次海侵时期最高海平面高度，消除构造升降和后期沉积的影响，可以反推出两次最大海侵古岸线的大致位置。全新世最大海侵时代约6 ka B.P.，当时海平面高于现今海平面约2 m（图7），估算全新世最大海侵以来研究区整体沉降了约4.2 m，根据钻孔估算全新世最大海侵以来研究区地层平均沉积厚度约为5.3 m，理论计算得到全新世最大海侵古岸线大致与研究区现今3 m等高线大致相符（图6）；晚更新世最大海侵时代约25 ka B.P.，当时海平面高度低于现今海平面约20 m（图7），晚更新世最大海侵以来研究区整体抬升了约12.8 m，根据钻孔估算晚更新世最大海侵以来研究区地层平均沉积厚度约为10 m，

因此，可以得到晚更新世最大海侵古岸线大致与研究区现今2.8 m等高线大致相符（图6）。这个结果进一步论证了研究区晚第四纪两次最大海侵古岸线走势相似、距离相近的主要原因，由于在全球海平面变化的大背景下，不同时期构造升降情况不一致，但整体上仍处于构造抬升的状况，从而导致的晚更新世海侵时期的古地势相较于全新世海侵时偏低，提供了发生海侵所具备的可容空间，即使最高海平面较低的情况下仍发育了范围与全新世海侵相近的海相沉积地层。

4.3 中国沿海晚第四纪海侵的新构造背景

晚第四纪以来，受冰期—间冰期控制，在中国沿海尤其是渤海、黄海和东海及华南沿海平原均发生了数次与全球气候和区域海平面变化相吻合的海侵事件[43]。前人关于渤海湾西侧和长江三角洲地区晚第四纪海侵研究已取得了大量的研究成果，晚第四纪以来渤海湾西侧华北平原发现了MIS5 期的沧州海侵、MIS3 期的献县海侵和MIS1 期的黄骅海侵，最大海侵范围献县海侵>沧州海侵>黄骅海侵[44-47]（图9a）。长江三角洲地区晚第四纪以来发生过3 次海侵，分别被命名为太湖海侵（晚更新世早期），滆湖海侵（晚更新世中晚期）和镇江海侵（全新世），其中镇江海侵范围最广，滆湖海侵次之，太湖海侵最小[48-50]（图9b）。渤海湾西岸和长江三角洲地区晚第四纪以来主要以持续的构造沉降为主[46,50]，前人研究发现，晚更新世早期本区古地势较高，物源供应较丰富，海水主要沿古河谷入侵，海侵规模较小[44-45]。但由于剥蚀作用和长期构造沉降，晚更新世中晚期和全新世期间该区极易受海平面波动影响而发生海侵[46]。即使MIS3 期海侵时最高海面低于全新世和晚更新世早期最高海面数十米，该时期最大海侵范围仍大于晚更新世早期时的海侵，这主要是由于构造沉降的结果。但在渤海湾西岸，虽然持续的构造沉降作用导致全新世古地势相对较低，但在8 500 a 左右黄河改道，从渤海入海[51]，为渤海的西海岸提供了大量的沉积物供给，从而补偿了构造沉降所增加的可容空间，导致黄骅海侵范围较小。

图9 晚第四纪最大海侵古岸线分布（a）渤海湾西岸；（b）长江三角洲（修改自Lin et al.[48]）Fig.9 Location of maximum coastline transgression in the Late Quaternary(a)west coast of Bohai Bay(b)Yangtze River Delta（modified from Lin et al.[48]）

福建沿海平原区晚更新世以来发生了两次海侵，为晚更新世的福州海侵和全新世的长乐海侵，福州海侵相对较弱，海水仅侵入到了福州市区一带，长乐海侵范围更广，其范围达到了福州西侧的闽侯一带[48,52]（图10a）。福建沿海整体处于新构造运动的上升区，但不同地区差异显著，表现为以闽江为界，南升北降[53]。闽江南侧地区海侵强度较弱，两次海侵最大范围相近，表现出与台山沿海相似的特征。在对珠江三角洲地区钻孔样品和微体古生物研究表明，晚更新世以来该区发育两套海侵地层[7-10]，全新世海侵主要发生在10～4 ka B.P.，最大海侵古岸线从现在的海岸向内陆延伸了约50 km[48]（图10b），晚更新世海侵最大范围小于全新世海侵[11,35]，此次海侵地层发育不普遍。珠江三角洲晚更新世海侵的时代归属有3～5 万年的晚更新世中晚期（相当于MIS3）[7,10]和10万年左右的晚更新世初期（相当于MIS5）[54]两种观点，同样也是影响三角洲地区新构造升降速率估算的重要因素。

图10 晚第四纪最大海侵古岸线分布（a）福建沿海；（b）珠江三角洲（修改自Lin et al.[48]）Fig.10 Location of maximum coastline transgression in the Late Quaternary(a)coastal zone of Fujian province；(b)Pearl River Delta(modified from Lin et al.[48])

我国海岸带跨越多个亚板块块体[55]，由于新构造差异活动，我国海岸带可分为山地丘陵海岸（福建沿海、台山）和平原海岸（渤海湾西岸、长江三角洲），前者基本上属于上升区，后者则属于下降区[48]。从以上粗略对比可以看出：1）中国沿海不同地区由于新构造背景不同，晚第四纪海侵发生次数和最大范围在空间上整体具有“北强南弱”的特点（表6）；2）中国沿海晚第四纪海侵发生的次数和最大范围差异显著，除不同构造单元存在着构造差异外，海平面高度和沉积物供给变化同样对海侵强度有着一定的影响，导致同一地区不同时期之间海侵强度差异并未完全遵循传统意义上“早弱晚强”[55]的规律（表6）；3）研究模拟表明，未来一段时间，全球海平面理论上以2.0 mm/a 作为预估的海平面上升速率较为合适[56]，加上人类对流域治理水平日趋成熟，流域沉积物供给会更加稳定，由此可见，今后海岸地区构造运动将成为影响海侵发生的重要因素。

表6 晚第四纪中国沿海部分地区海侵强度对比Table 6 Comparison of transgression intensity in Chinese coastal areas in the Late Quaternary

5 结论

（1）珠江三角洲西缘台山地区晚第四纪地层自下而上可划分为4个沉积单元（石排组、西南镇组、三角组、横栏组），经历了晚更新世中晚期的河流—海相交替沉积环境、末次冰期的海退暴露环境以及全新世滨海、浅海沉积环境。

（2）在MIS3和MIS1的高海面时期发生了2次主要的海侵事件，最大海侵古岸线最远分别向陆地延伸15.4 km 和16.3 km，两条最大海侵线距离接近，走势相同。两次海侵时期实际最高海平面高程之差达数十米，而两条最大海侵古岸线却极为接近，造成这种现象的原因主要是研究区晚更新世以来整体处于构造抬升状态，使得晚更新世海侵时的古地势较低，从而导致MIS3 期最高海平面较低的情况下仍发育了范围与全新世海侵相近的海相沉积地层。

（3）两期海侵沉积物的高程和年代资料揭示了研究区晚更新世以来整体抬升速率为+0.51 mm/a，全新世中晚期转变为沉降状态，沉降速率为-0.61 mm/a，MIS3 和MIS1 最大海侵古岸线分别与研究区当前2.8 m和3 m等高线相符。

（4）中国沿海晚第四纪海侵发生次数和最大范围在空间上整体具有“北强南弱”的特点，同一地区不同海侵之间不完全遵循“早弱晚强”的特点。构造运动、海平面变化和沉积物供给共同影响海侵的发生和最大海侵古岸线的分布，随着全球海平面上升速率和流域沉积物供给量趋于稳定，构造升降将会在未来成为应对海平面上升对海岸带地区的影响中不可忽视的重要影响因素。

致谢 样品测试过程中受到贺洋、赵珂、贾冀新、魏子谦、董鑫等同学的协助，在此一并感谢。