0 引言
黑潮是北太平洋最强的西边界流,由北赤道暖流在菲律宾群岛以东分支汇合形成,沿台湾东侧北上,经琉球群岛进入冲绳海槽,最终将巨量热量和盐分输送至中高纬度地区(Jayne
et al.,
2009; Fenies
et al.,
2023)。其主轴位于西太平洋海气热交换最强烈的区域,对东亚气候、上层海洋热结构以及该地区的海洋沉积物分布有着强烈的影响( Jian
et al.,
2000; Qiu,
2001; Hu
et al.,
2015; Yang
et al.,
2024)。同时,黑潮有可能通过海—气相互作用显著影响东亚季风系统,进而影响中国东部的降水格局和极端气候事件(Sasaki
et al.,
2012; 殷明等,
2016;Gan
et al.,
2019)。
近年来,前人多基于沉积和地化指标重建黑潮的活动规律与演化路径。根据有孔虫碳同位素梯度(
Δδ13C)及表层海水温度分别从温跃层深度、水体分层与水团温盐性质等角度揭示上层水团结构对黑潮强度演化的响应(Ujiié and Ujiié,
1999; Jian
et al.,
2000; Ujiié
et al.,
2003;Li
et al.,
2020)。对沉积物粒度和Sr-Nd同位素的分析指示黑潮的流速大小和移动路径(McCave
et al.,
1995; Diekmann
et al.,
2008,Dou
et al.,
2010; Zheng
et al.,
2016)。同时也有利用沉积物的磁学特征讨论粗颗粒磁性矿物的输运能力对黑潮的响应,认为黑潮增强时粗颗粒组分输入增加(Wang
et al.,
2025b)。综合上述指标研究发现,黑潮强度在末次冰盛期整体较弱,自冰消期开始呈现阶段性增强趋势,至全新世达到强盛状态(Shi
et al.,
2014)。然而,已有研究多侧重单一维度,在黑潮活动时上层水团结构和流速与搬运能力的综合特征认识相对不足,不同钻孔不同指标重建的黑潮演化时空响应差异显著。如最新重建结果揭示,末次冰盛期黑潮主路径可能偏移或未进入冲绳海槽,冰后期随海平面回升重新入侵并持续增强(Wang
et al.,
2025a),且黑潮的强度变化与路径迁移存在显著的空间异质性(Huang
et al.,
2025)。
此外,现代观测与数值模拟证实,在年际—年代际尺度上,黑潮的变化可以通过调节向西太平洋暖池的热量积累及向东亚大陆的水汽输送,对东亚季风降水产生一定影响(Qiu and Chen,
2010;Taguchi
et al.,
2010; Gan
et al.,
2019)。近期定量研究进一步显示,黑潮对东亚降水的影响在中国东部梅雨带尤为显著: 东海黑潮增温可通过强化对流层低层风场辐合与水汽聚集增强梅雨降水(Qiao
et al.,
2024),过去40年间黑潮温度升高对梅雨期降水增强的贡献率可达37%(Qiao
et al.,
2025)。在长江流域的全新世气候演变研究中,也发现了降水与黑潮强度从早全新世显著波动向中晚全新世振荡性减少的阶段性转变模式(Chen
et al.,
2024)。但是在较长时间尺度,在轨道—千年尺度上黑潮演化与中国东部降水之间直接耦合关系及其主导机制的认识仍有待进一步厘清。
基于此,本研究针对冲绳海槽中南部获取的沉积钻孔CF5,通过对海水表层种(Globigerinoides ruber,G.ruber)和次表层种(Pulleniatina obliquiloculata,P.ob)有孔虫碳同位素分析,结合沉积物粒度及表征磁性矿物颗粒大小和含量变化的磁学参数,构建黑潮综合活动指标,系统重建末次盛冰期至早全新世期间冲绳海槽中南部黑潮的演化特征,揭示轨道—千年尺度上黑潮演变与中国东部降水的耦合关系,并结合热带太平洋东西海温梯度记录,探讨黑潮活动的驱动机制,以期为深入理解西边界流在东亚水文气候变化中的作用提供关键的地质证据。
1 区域地质背景
冲绳海槽位于琉球岛弧与亚洲大陆边缘之间,其形成与菲律宾海板块和欧亚板块的碰撞与俯冲密切相关(Sibuet
et al.,
1987; Diekmann
et al.,
2008; Liu
et al.,
2016)。海槽呈弯曲状自北向南延伸,最深处可达2500 m,水深一般超过1000 m(Diekmann
et al.,
2008)。海槽两侧具有陡峭的边缘,底部则呈现平坦形态,由琉球岛弧与琉球海沟构成完整的海沟岛弧盆地体系(Li
et al.,
2023)。起源于西太平洋赤道附近的高温、高盐表层暖流,沿菲律宾东部海域转向之后沿台湾东岸北上,形成黑潮,是北太平洋副热带环流的核心西边界流(Hu
et al.,
2015)。黑潮在进入冲绳海槽后沿大陆坡北上(
图 1),其主流轴位置和水动力强度对该区域的海洋环境具有主导性控制作用(Wang
et al.,
2016; Hsueh,
2000)。黑潮夏季表层水温可达28 ℃以上,盐度通常在34.5‰~35‰之间,流速快、流量大、具有较强的搬运能力,是将低纬度的高热量输送到中高纬度的重要输送带。现代观测数据表明,黑潮在冲绳海槽段的表层流速可达1.5~2.0 m/s,并受到季节性风应力场、海洋温盐结构等因素的调节(Sasaki
et al.,
2012; Zhao
et al.,
2015)。在冰期—间冰期气候转换过程中,黑潮的强度和空间上的南、北摆动,与El Niño等海—气耦合过程相关,深刻影响着中国东部沿海区域的气候和沉积过程(Hu
et al.,
2015)。
研究区沉积环境受物质来源、海洋动力条件和区域构造活动等多重因素控制(Ujiié
et al.,
2003; Bian
et al.,
2010; Dou
et al.,
2010; 姜韬等,
2015)。沉积物主要来自邻近的东海陆架、台湾造山带以及琉球火山岛弧,通过河流输入、重力搬运和洋流输送进入海槽(Bian
et al.,
2010; Zheng
et al.,
2016; 王越奇等,
2018; Wang
et al.,
2025a)。黑潮作为该区域的主控动力,其流速和路径变化直接影响沉积物的输运、分选和堆积过程: 流速增强时,水动力增强可将粗颗粒物质向更远区域搬运; 流速减弱时,沉积物沉降速率和堆积厚度相应变化(Ellwood and Ledbetter,
1977; Qiu,
2001; Qiu and Chen,
2010; Zhao
et al.,
2015)。这种水动力调控在沉积记录中表现为粒度组成、磁学参数及矿物特征的变化,为重建区域古海洋环境提供了重要载体(Thompson and Oldfield,
1986)。
2 材料与方法
2.1 研究材料
本研究选取的沉积柱状样(CF5钻孔)位于冲绳海槽中南部(125.33°E,26.57°N;水深: 1015 m),钻孔总长514 cm,是由青岛海洋地质研究所于2020年使用重力柱状采样器获取。钻孔岩性以青灰色粉砂质黏土为主,沉积连续,在岩性上未发现明显的沉积扰动和间断,处于现代黑潮主流轴影响范围内,是研究黑潮演化的理想站位。将钻孔沿中间切分之后,以约2 cm为间距,连续采集249个样品进行磁学测试; 以相同间距同步采集249个样品用于粒度分析,测试过程中个别样品损耗未获得有效数据,但不影响整体粒度曲线的趋势与解释; 以约6 cm为间距,采集76个样品进行有孔虫碳同位素分析。
2.2 研究方法
2.2.1 碳同位素分析
在CF5钻孔样品中挑选浮游有孔虫表层种
Globigerinoides ruber(
G. ruber)和温跃层种
P. obliquiloculata(
P. ob)壳体进行稳定碳同位素分析。测试样品个体要求壳体完整、表面洁净且无明显溶蚀痕迹。挑选出的壳体经玻璃板轻压破碎后,在10% H
2O
2中浸泡12 h去除有机质,随后在乙醇溶液中用超声波清洗器(40 Hz)震荡3~5 s去除残留颗粒,重复清洗3次,直至上清液不再浑浊,最后置于60 ℃烘箱中烘干。将预处理后的有孔虫个体转移至Gas bench Ⅱ样品瓶中,加入100% 磷酸过夜反应,生成的CO
2气体导入Thermo Fisher Delta V Plus稳定同位素质谱仪进行测定。测试工作在中山大学地球科学与工程学院古地磁实验室完成。测试结果按国际标准转换为Vienna Pee Dee Belemnite(VPDB)量表,根据NBS 18和04405标准进行精度校准,分析精度优于±0.06‰。由于黑潮活动强烈时,水体分层增强,不同水深有孔虫的同位素比值的差值能够反映水体结构的变化,我们采用$\mathit{\Delta }{\mathit{\delta }}^{13}{{{\mathit{C}}_{\mathit{G}.}}_{\mathit{r}\mathit{u}\mathit{b}\mathit{e}\mathit{r}-\mathit{P}.}}_{\mathit{o}\mathit{b}}$=${\mathit{\delta }}^{13}{{\mathit{C}}_{\mathit{G}.}}_{\mathit{r}\mathit{u}\mathit{b}\mathit{e}\mathit{r}}$-${\mathit{\delta }}^{13}{{\mathit{C}}_{\mathit{P}.}}_{\mathit{o}\mathit{b}}$,反映表层与温跃层之间的碳同位素差异。
Δδ13C增大指示温跃层加深、上层海洋分层加强,通常对应黑潮增强阶段; 反之则指示温跃层变浅、分层减弱(Ujiié
et al.,
2003)。
2.2.2 粒度分析
称取约0.2~0.3 g沉积物样品(经预实验确定该质量可满足仪器遮光度要求),依次加入10%双氧水和10%盐酸溶液去除有机质与碳酸盐,反应完全后用超纯水离心清洗至pH呈中性。上机测试前加入0.05 mol/L六偏磷酸钠溶液并超声震荡5 min使颗粒充分分散。粒度分析在中山大学地球科学与工程学院海洋地质实验室完成,使用Mastersizer 3000型激光粒度分析仪。测试过程中仪器搅拌器转速为2500转/min,并伴随15 s超声处理,样品之间执行超声清洗以防交叉污染。粒度测量范围为0~3500 μm,重复测量相对误差小于5%。
2.2.3 磁学参数测量
对CF5钻孔所有样品进行等温剩磁(IRM)测量。首先给样品施加1 T的脉冲场使样品饱和磁化,测得样品的饱和等温剩磁(SIRM),随后逐步施加梯度交变场退磁(2.5~80 mT)并测量剩磁。测量工作在中国地质大学(武汉)地球科学学院环境磁学实验室完成,使用2G755-4k U-Channel SQUID超导岩石磁力仪,灵敏度为1.0×10-12 Am2。
对所有样品在AGICO KLY-5 Kappabridge磁化率仪进行同相位(in-phase)和异相位(out-of-phase)体积磁化率的测试(工作频率1220 Hz,外加场强400 A/m)。异相位磁化率源于磁性颗粒在交变磁场中的磁弛豫滞后效应,其数值大小对处于超顺磁/单畴边界的细小磁性颗粒尤为敏感,功能上与频率磁化率(
χfd)等效,可指示超顺磁颗粒的相对丰度(Hrouda
et al.,
2013)。实验在中山大学地球科学与工程学院古地磁实验室完成。
3 年代框架
为获取精确的沉积物年龄信息,作者根据沉积物粒度的变化趋势,从CF5沉积物柱状样中自上而下选取岩性均一、粒度较细的层位,挑选18个浮游有孔虫壳体样品分别送往美国Beta实验室和中国科学院广州地球化学研究所进行放射性碳同位素加速质谱仪测年(AMS14C)。
测年数据显示,不同层位AMS
14C年龄存在明显的新、老倒置现象(
表 1中带
*的样品),这与冲绳海槽高沉积速率地区普遍存在的沉积物二次搬运过程密切相关。末次冰消期以来,海平面快速回升导致东海陆架被大面积淹没,陆架上原有的老沉积物被侵蚀再搬运至冲绳海槽再沉积(Dou
et al.,
2010; Zheng
et al.,
2016),造成了不同层位中浮游有孔虫壳体年龄的新老倒置。Broecker等(
1989)指出,在高沉积速率的海洋沉积物中,再搬运作用是造成有孔虫AMS
14C年龄异常的主要原因之一,剔除明显偏老的非原位年龄是建立可靠年代框架的必要步骤。
在年龄控制点的筛选过程中,本研究遵循以下标准: (1)优先选取自下而上年龄随深度递减的测年点,符合地层叠覆律; (2)剔除明显偏离整体年龄—深度趋势的异常偏老样品(
表 1中带
*标记),这些样品的偏老年龄由二次搬运的有孔虫加入所致。综合上述标准,从18个测年点中剔除9个受二次搬运影响的样品,最终选取9个测年样品建立年龄深度模型(见
表 1)。测年结果使用CALIB 8.1.0软件进行日历校正,校正曲线选用Marine 20(Heaton
et al.,
2020),应用当地海洋储库效应校正值
ΔR=-150±59
14C yr(Wang
et al.,
2025b)。年代—深度模型基于Bacon v2.5.6构建(默认分层厚度为5 cm;Blaauw and Christen,
2011),采用贝叶斯统计方法并考虑沉积速率的时变性。结果表明(
图 2),CF5钻孔记录了自约18.9—2.9 ka期间较连续的沉积过程,沉积速率介于20.8~49.6 cm/kyr之间,在约7.0 ka时达到峰值,平均沉积速率为33.2 cm/kyr。
鉴于本研究聚焦末次盛冰期至早全新世黑潮的阶段性演化及其气候效应,且钻孔上部(约6 ka以来)仅有2个年龄控制点(8 cm和48 cm,
表 1),沉积速率出现异常波动,可能存在由生物扰动、垂向混入导致的AMS
14C年龄偏差(Broecker
et al.,
1989; Zheng
et al.,
2016)。此外,浮游有孔虫碳同位素数据在约6 ka以来分辨率较低,难以支持精细的古海洋重建。因此,讨论分析主要围绕约18.9—6 ka时段展开,中晚全新世黑潮变率及其与东亚降水的耦合关系有待今后更高分辨率和更精确年代控制的工作进一步厘清。
4 结果
4.1 有孔虫碳同位素
浮游有孔虫
G.ruber主要生活在表层水体,而
P.ob是温跃层种属,生活在约100 m左右水深,两者壳体
Δδ13C值的变化反映了上层水团中碳的混合均一程度和水体分层情况(Jian
et al.,
2000; Jing
et al.,
2024)。由于黑潮以高温、高盐和低营养盐为特征,当黑潮侵入时,水体中营养盐组分相应会降低,因此表层浮游有孔虫壳体的碳同位素相应地会偏正,表层和温跃层之间碳同位素的梯度将加大(Ujiié
et al.,
2003; Chen
et al.,
2011; Xu
et al.,
2021)。CF5钻孔中
G. ruber碳同位素呈现3段波动变化: 阶段Ⅰ(18.9—12 ka)逐渐升高; 阶段Ⅱ(12—8.8 ka)波动降低; 阶段Ⅲ(8.8—6 ka)快速波动恢复到与之前相近的水平(
图 3-a)。
P. ob的总体规律与
G. ruber相反,伴随更为显著的千年尺度的波动(
图 3-b)。即当表层浮游有孔虫碳同位素偏正时,温跃层种属的碳同位素则偏负,二者之间的梯度
Δδ13C加大(
图 3-c)。这种变化反映了当黑潮入侵时,温跃层加厚,海水分层增强,使表层与次表层的碳同位素之间的混合作用变弱,梯度增加。
4.2 黑潮活动磁学指标构建
冲绳海槽粗颗粒碎屑主要源自东海陆架和台湾造山带(Bian
et al.,
2010; Dou
et al.,
2010),黑潮流速变化直接影响其搬运能力,其增强时搬运更多粗颗粒物质(Diekmann
et al.,
2008; Li
et al.,
2014; Zheng
et al.,
2016)。在同一类型的磁性矿物中,剩磁矫顽力的大小可以反映颗粒的粗细变化。粗颗粒磁铁矿矫顽力低,其剩磁在低场区间可以被快速清洗; 细颗粒磁铁矿矫顽力较高,携带的剩磁在较高的交变场才能逐步清除(Thompson and Oldfield,
1986)。分析CF5钻孔样品SIRM的交变退磁曲线谱(
图 4-a),可以发现在2.5~10 mT和30~60 mT之间,退磁曲线的衰减过程存在2个显著的斜率,前者快速的衰减代表了粗颗粒磁性矿物,而后者则反映了较细磁性颗粒携带剩磁的变化。基于这一现象,SIR
M@(30-60/(2.5-10)mT比值(以下简称S
soft比值)可以反映沉积物中磁性矿物颗粒粗、细组分含量的相对变化,其可能表示了搬运动力强弱变化,可以作为黑潮强度的替代指标,即S
soft比值越低,指示粗颗粒磁性矿物相对含量越高,反映黑潮流速增大、搬运能力增强。但是细颗粒磁性矿物可能会受早期成岩作用的影响而溶解,使原始的颗粒大小信息发生变化,本研究对早期成岩作用可能的效应进行分析评估。
早期成岩过程中的还原作用优先溶解细小的超顺磁(SP)颗粒(Karlin and Levi,
1983),若本钻孔沉积物经历了显著的还原改造,对SP颗粒丰度敏感的
χfd值则会显著减小。由于SIRM只表示大于SP颗粒磁性矿物携带的剩磁,如果还原作用显著,
χfd-SIRM之间就不可能呈现线性关系。分析CF5钻孔所有样品的
χfd-SIRM关系,二者在全孔范围内呈现显著的线性正相关关系(
图 4-b),表明SP颗粒与其他颗粒均主要受碎屑输入控制,SP颗粒保存良好,早期成岩溶解作用对磁性矿物粒级组成的改造有限。这一结果与本钻孔较高的沉积速率(平均33.2 cm/kyr)及区域已有研究表明冲绳海槽磁学参数变化主要受控于碎屑输入而非成岩改造(Liu
et al.,
2007;Li
et al.,
2019)的结果相一致。此外,区域物源研究表明,冲绳海槽粗颗粒磁铁矿主要源自台湾造山带的碎屑输入(Li
et al.,
2019),黑潮动力变化调控的是该物源端元的输入通量,而非物源区的根本性改变。因此,磁学信号的变化本质上反映的是水动力强度对物源混合比例的调制。
在上述可靠性检验的基础上,S
soft比值呈现4阶段变化特征(
图 5-a): 阶段Ⅰ(18.9—17 ka),S
soft比值升高,说明此时黑潮流速低,搬运能力弱; 阶段Ⅱ(17—13 ka)其值逐渐减低,表示粗颗粒磁性组分含量逐渐增加,黑潮流速加大、搬运能力逐渐增强; 阶段Ⅲ(13—10.5 ka)期间再次增加,阶段Ⅳ(10.5—6 ka)又开始整体上降低,反映了黑潮流速的再次减弱和增强的过程。
综合分析有孔虫碳同位素值(
图 5-b)和粗颗粒磁性矿物组分的变化,可以发现尽管两者在总体上指示了黑潮活动在末次冰消期之前(约17 ka)较弱,之后逐渐增强,全新世基本保持小幅波动的特征,但是二者反映的黑潮入侵和搬运动力大小方面,在千年尺度上存在一定差异。即当黑潮入侵时,其流速并不一定显著增强,这一特征在末次冰消期和早全新世体现的尤为明显。为了统一表示其活动特征,本研究尝试将磁学参数和有孔虫碳同位素梯度指标各自归一化,然后等间距插值。考虑到S
soft比值与黑潮强度呈负相关,对其进行反向处理:
SIRMrev=1-
SIRMnorm 后求取二者的平均值作为黑潮活动的一个综合指标(
IKC=0.5×
Δδ13Cnorm+0.5×
SIRMrev)(
SIRMnorm为归一化后的剩磁比值,
SIRMrev为反向处理后的
SIRMnorm;
Cnorm为归一化后的有孔虫碳同位素差值,
IKC为构建的黑潮综合指标)(
图 5-c)。尽管该指标由于碳同位素数据分辨率相对不足存在一定误差,但在总体上仍能反映黑潮在冲绳海槽活动的历史。
4.3 碎屑物质粒度变化
为进一步验证上述指标的可靠性,准确重建黑潮演化信息,作者通过AnalySize软件(Paterson and Heslop,
2015)中的非参数化端元建模算法(HALS-NMF)对钻孔沉积物粒度进行分解。经测试,当端元数为2时拟合效果最佳,所有样品的拟合决定系数(
R2)均大于0.98(平均约0.99),角距离(
θ)普遍较小,表明两端元模型已能充分表征粒度数据的主要变异特征分析。其中EM1组分众数粒径约6~12 μm,代表细颗粒组分,EM2组分众数粒径约40~100 μm,代表粗颗粒组分。发现粗粒组分(EM2
coarse)含量与黑潮综合指标
IKC在长时间尺度上呈现一致的波动趋势(
图 5-d),均记录到了末次盛冰期的弱水动力、冰消期的波动增强及早全新世的强盛状态。且在冷热事件期间均表现出一致的波动特征,表明粒度对千年尺度的黑潮减弱或增强事件同样具有清晰响应。尽管在个别时间如约16 ka,粒度中的粗粒组分含量突然大幅下降,而
IKC仅小幅回落,这可能反映了冰消期初期黑潮水团增强先于搬运能力建立的异步响应,
Δδ13C的上升部分抵消了
IKC的降幅。这进一步支持了
Δδ13C和S
soft作为黑潮演化代用指标的可靠性。
5 讨论
5.1 末次盛冰期至早全新世冲绳海槽中南部黑潮的阶段性演化
综合分析S
soft比值和有孔虫碳同位素特征,本研究可以在整体上刻画黑潮在冲绳海槽18.9~6 ka期间在入侵和搬运动力等方面的演化特征。尽管由于可能源自沉积二次搬运的影响,AMS
14C存在一定的不确定性,但是从年龄模型分析,仍然可以在总体上约束数千年尺度的黑潮演化的讨论,尤其冰消期以来几个关键的气候事件,如BA暖期(14.7—12.9 ka)和YD冷事件(12.9—11.7 ka),有较好的年龄控制。具体可以将黑潮在此时期内的活动规律分为以下3个阶段(
图 6)。
阶段Ⅰ(18.9—17 ka): 末次盛冰期弱黑潮阶段。此阶段
IKC逐渐降低,
Δδ13C较小、S
soft比值较高,指示温跃层较浅、上层海洋分层较弱,黑潮流速低、搬运能力弱。区域钻孔的粒度指标同样显示水动力减弱(Diekmann
et al.,
2008),本研究中粒度含量的同步变化更为明显。此时海平面较现代低约120 m(Lambeck
et al.,
2014; Li
et al.,
2014),东海陆架大面积暴露。
阶段Ⅱ(17—10.5 ka): 冰消期波动增强阶段。IKC呈波动上升趋势,指示黑潮活动逐渐增强,温跃层开始加深、上层海洋分层加强,流速增大。期间可进一步识别出H1事件后的缓慢增强、BA暖期的快速增强以及YD事件时的短暂减弱。
其中BA暖期
Δδ13C值相对减小,但粗颗粒组分增加,反映黑潮流速增大而温跃层实际变浅。同时,综合区域对比结果显示,在此期间黑潮响应存在空间差异: 南部记录显示黑潮搬运能力增强,中部和北部相对减弱(Diekmann
et al.,
2008; Shi
et al.,
2014; Zheng
et al.,
2014,
2016),可能反映了黑潮主流轴的摆动或不同纬度对驱动力的敏感性差异(Jin
et al.,
2010;武杰等,
2020)。磁学指标与不同区域钻孔粒度指标呈现一定的差异,说明了从南向北,黑潮搬运粗碎屑能力的减弱和不同碎屑颗粒对黑潮活动的响应特征不同,沉积物粒度的可排序组分,可能更多地混合了来自陆架的碎屑物质,在一定程度上弱化了黑潮信息。从表层海水温度梯度重建的黑潮南、北空间扩展范围看(Liu
et al.,
2022),在BA暖期黑潮向北扩展和之后的向南收缩(Yamamoto,
2009),与磁学参数指示的黑潮活动强度也一致,向北扩展时搬运能力增强,而向南收缩时搬运能力减弱。
YD事件期间,IKC明显回落,与区域各钻孔记录的黑潮减弱特征一致,但CF5钻孔Δδ13C在YD初期出现的短暂升高可能反映了黑潮在该时期水体温盐性质与强度的差异。即该时期冲绳海槽的水体物理性质体现的黑潮特征非常显著,但是流动速度偏弱,搬运能力下降。
此外,冰消期前期,随着
Δδ13C逐渐升高,南、中、北部区域基于沉积物粒度的指标尽管存在多个高频波动(Fenies
et al.,
2023),但没有体现出黑潮逐渐增强的特征。这种差异可能源于东海陆架碎屑物质和由黑潮搬运的碎屑物质的混合影响,弱化了黑潮搬运可排序粒径碎屑物质携带的信息。
阶段Ⅲ(10.5—6 ka): 早全新世强盛阶段。
IKC在YD事件后快速反弹,至约9.5 ka左右达到整个序列的最高值,指示黑潮输运能力达到最强,对应温跃层达到最深、上层海洋分层最为稳定。这一时期黑潮将最强盛的暖水输送至冲绳海槽,表层增温显著,稳定的浮力层抑制了垂直混合,导致温跃层持续加深(Ujiié
et al.,
2003; Kubota
et al.,
2010)。此时海平面已上升至接近现代水平,东海陆架被完全淹没,黑潮路径趋于稳定,对冲绳海槽的影响达到最强(Xu
et al.,
2021)。区域对比显示,南部和中部钻孔在8.4 ka事件期间同样记录了黑潮的短暂减弱,而北部响应较弱。
5.2 黑潮演化对中国东部降水的影响及驱动机制探讨
黑潮作为北太平洋副热带环流的核心西边界流,通过调控西太平洋暖池热量积累及向东亚大陆的水汽输送,对东亚季风降水产生重要影响(Liu
et al.,
2024)。现代观测与模拟研究已揭示,在年际—年代际尺度上,黑潮强度与ENSO背景态密切相关: 类La Niña态下黑潮增强,为东亚夏季风提供更充沛的水汽来源; 类El Niño态下黑潮减弱,东亚降水相应减少(Qiu and Chen,
2010; Taguchi
et al.,
2010)。且其对中国东部梅雨带的影响尤为显著,研究已证实东海黑潮增温通过强化水汽辐合增强梅雨降水(Qiao
et al.,
2024,
2025)。然而,这一“ENSO-黑潮—东亚降水”遥相关机制在轨道—千年尺度上是否依然稳定,尚需地质记录检验。因此,本研究选取梅雨带内大九湖(109.9~110.2°E,31.5~31.6°N,西部边缘; Xiao
et al.,
2015; Sun
et al.,
2019)、巢湖(117.4°E,31.6°N,核心区; Li
et al.,
2018)和新街遗址(119.7°E,31.4°N,核心区; Lu
et al.,
2019)3条孢粉降水重建记录与本研究对比分析,以检验长时间尺度上黑潮与区域降水的关联。
基于冲绳海槽中南部钻孔CF5重建的黑潮综合活动指数(
IKC,
图 7-a)显示,末次盛冰期至早全新世期间冲绳海槽中南部黑潮演化记录与热带太平洋东西海表温度梯度ΔSST(
图 7-b)高度同步。定量分析表明:
IKC与ΔSST呈显著正相关(Pearson
r=0.60,
p<0.0001;Spearman
ρ=0.44,
p=0.002;
n=47),证实现代观测所揭示的ENSO对黑潮的调制作用在更长尺度上同样有效,热带太平洋背景态是调控黑潮长期演变的关键因素(Bond
et al.,
2001)。结合3条降水记录的时间覆盖,将黑潮与降水的对比划分为3个时段。
阶段Ⅰ(18.9—17 ka): 此阶段仅大九湖记录有效。
IKC、ΔSST与大九湖降水均呈降低趋势(
图 7-a至
7-d),指示弱黑潮、类El Niño背景态与低降水的耦合。此时西太平洋暖池热量积累减少、副热带环流减弱(Sun
et al.,
2005; Liu
et al.,
2024),同时暴露的陆架改变了黑潮的流动路径,使其更偏向东侧(Ujiié
et al.,
2003; 赵子荟等,
2023),向东亚大陆的水汽输送通量显著降低。
阶段Ⅱ(17—12 ka): 该时期西太平洋暖池热量开始积累,驱动黑潮呈现波动增强趋势(Cane,
2005; Qiu and Chen,
2005; Xie
et al.,
2009; Taguchi
et al.,
2010),
IKC和ΔSST大幅波动上升,BA暖期出现增强,反映了变暖背景下副热带环流系统增强,随后在YD期间短暂回落。大九湖降水整体小幅上升,这表明随着冰盖消融和海平面上升,黑潮主流轴逐渐西移,向陆地输送暖湿水汽(Diekmann
et al.,
2008; Zheng
et al.,
2014,
2016),为东亚夏季风提供了充沛的水汽来源,从而驱动了区域降水的增加。然而,大九湖降水在H1事件后短时间内和YD期间与
IKC存在明显解耦,且新街遗址降水记录在YD期间恰与
IKC、ΔSST同步下降(
图 7-c),这一差异暗示了黑潮对降水影响的空间不均一性及其他气候驱动因子的补偿作用。
阶段Ⅲ(12—6 ka): 此阶段根据巢湖与新街记录重建的2条降水曲线与
IKC、ΔSST呈现近乎同步的演化——冰消期末至早全新世同步波动上升,在约10 ka和7.5 ka先后达到峰值,是西太平洋暖池充分发展、东西海温梯度达到峰值、海平面上升至接近现代水平的综合结果(Linsley
et al.,
2010; Xu
et al.,
2021)。此时,东海陆架被完全淹没,黑潮路径趋于稳定,西太平洋暖池为西边界流的增强提供了充足的热量积累,其对冲绳海槽的影响达到最强。梅雨带核心区2条独立记录与
IKC具有高度一致性,表明黑潮对梅雨带降水的调控在该时段内占主导地位。这种“强黑潮—多降水”的对应关系,清晰地揭示了在适宜的海平面和气候背景下,强盛的西边界流能够有效地将热带海洋的水分向北输送,为东亚夏季风降水提供关键的水汽支持。
综合上述分段对比,本研究发现黑潮对中国东部降水的影响存在显著的空间不均一性。在YD期间,靠近冲绳海槽的梅雨带核心区(新街,119.7°E)降水与IKC同步下降,对黑潮减弱响应敏感; 而西部边缘(大九湖,109.9~110.2°E)降水未显著下降,对黑潮信号的响应相对弱化。12—6 ka期间核心区降水与IKC高度同步、而17—12 ka期间边缘区降水与IKC多次解耦的对比,进一步印证了这一空间分异。究其机制,黑潮对中国东部降水的影响存在自东向西的梯度变化。
在靠近冲绳海槽的梅雨带核心区,降水变化主要受热带太平洋背景态主导下的黑潮水汽输送调控。该区域距黑潮水汽源区近,太平洋水汽在降水来源中占比较高,ENSO-黑潮信号能够直接、高效地传递至降水记录中,表现为
IKC与降水的同步耦合。在远离冲绳海槽的梅雨带西部边缘,黑潮水汽输送随距离增加而衰减,降水的水汽来源更为多元,南海水汽的贡献显著增强(江志红等,
2013)。因此,该区域降水对黑潮信号的响应相对弱化,而季风系统自身的强度变化、高纬过程的遥相关调制等其他因子的作用则更为凸显。例如,H1事件后短期内和11.5—10.5 ka期间黑潮与大九湖降水的解耦,可能正是由于高纬融冰事件导致AMOC减弱、通过大气遥相关激发东亚夏季风增强(Jian
et al.,
2000; Stott
et al.,
2002; Wassenburg
et al.,
2021; Zhang
et al.,
2021),季风增强带来的水汽输送增加部分抵消了黑潮减弱或滞后增强的影响,从而在降水记录中掩盖了黑潮信号。类似地,H1事件后短期内黑潮增强而降水未同步响应,亦可能因AMOC减弱背景下东亚夏季风系统尚未完全建立(Cheng
et al.,
2016),对黑潮水汽的“承接”能力有限。此外,东亚季风系统自身对黑潮也存在反馈式调制: 夏季风增强可促进黑潮延伸体北扩,冬季风增强则可能抑制黑潮系统(Sasaki
et al.,
2012; Yu
et al.,
2019),这种调制在全新世尤为显著(Zheng
et al.,
2014)。在梅雨带边缘区,季风对降水的直接调控可能进一步弱化了黑潮信号的辨识度。
综上所述,末次盛冰期至早全新世期间黑潮强度与梅雨带降水在轨道尺度上呈“强黑潮—多降水”的耦合关系,但黑潮对中国东部降水的影响具有显著的空间异质性: 自东向西随距离增加而逐渐减弱,梅雨带核心区由ENSO-黑潮主导,边缘区则受多因子协同调制。在千年尺度上,边缘区降水与黑潮的解耦反映了季风系统独立调制和高纬过程介入的叠加效应。这种“热带太平洋背景态主导长期趋势,季风调制高频波动,水汽输送空间衰减”的多因子协同机制,为理解黑潮—降水关系的复杂模式提供了完整的理论框架,深入认识这一模式,对预测未来区域水循环变化具有重要意义。
6 结论
基于冲绳海槽中南部CF5钻孔的表层浮游有孔虫(G. ruber)和温跃层种有孔虫(P. ob)壳体的碳同位素、磁学参数及粒度分析,本研究构建了黑潮综合活动指数(IKC),重建了末次盛冰期至早全新世(约18.9~6 ka)黑潮的阶段性演化历史,并结合梅雨带内多条古降水记录,探讨了黑潮对中国东部降水的影响及驱动机制。主要结论如下:
1)末次盛冰期至早全新世期间冲绳海槽中南部黑潮演化呈现显著的3阶段特征: 末次盛冰期(18.9—17 ka)黑潮整体较弱,温跃层浅、搬运能力低; 冰消期(17—12 ka)黑潮波动增强,BA暖期出现显著增强脉冲,YD期间短暂减弱; 冰消期末至早全新世(12—6 ka)黑潮持续增强,约9.5 ka前后达到鼎盛,温跃层最深、输运能力最强。
2)黑潮的长期演化受热带太平洋背景态主导。IKC与热带太平洋东西海表温度梯度(ΔSST)呈显著正相关,类La Niña背景态有利于黑潮增强,类El Niño背景态则对应黑潮减弱,表明ENSO对黑潮的调制作用在轨道—千年尺度上稳定存在。
3)黑潮对中国东部降水的影响具有显著的空间不均一性。在靠近冲绳海槽的梅雨带核心区(如新街、巢湖),降水受ENSO-黑潮水汽输送主导,与IKC高度同步; 而在远离海槽的西部边缘(如大九湖),黑潮信号随距离衰减,降水响应相对弱化,更易受南海水汽补给、季风系统独立变化及高纬融冰事件等因子的叠加调制。这种“热带太平洋背景态主导长期趋势、季风系统调制高频波动、水汽输送空间衰减”的多因子协同机制,为理解黑潮—降水关系的复杂模式提供了完整的理论框架。
*国家自然科学基金项目(42330504)
国家自然科学基金项目(42274090)
国家自然科学基金项目(42576072)
山东省自然科学基金项目(ZR2024MD018)
中国地质调查项目(DD20230647)
中国地质调查项目(DD20240301801)