晚上新世以来菲律宾海盆XT4孔黏土矿物特征及其古环境意义

丁雪，胡邦琦,2,3，徐方建，郭建卫，崔汝勇，易亮

1. 中国地质调查局青岛海洋地质研究所，青岛 266071

2. 青岛海洋科学与技术试点国家实验室海洋矿产资源评价与探测技术功能实验室，青岛 266237

3. 中国科学院地球环境研究所黄土与第四纪国家重点实验室，西安 710061

4. 中国石油大学（华东）地球科学与技术学院，青岛 266580

5. 同济大学海洋与地球科学学院海洋地质国家重点实验室，上海 200092

风尘“源-汇”过程是指源区风尘释放、区域至全球尺度的输送和沉降，是地球不同圈层间（陆地-大气-海洋）相互作用的重要纽带，也是影响全球变化及其环境效应的一个关键过程[1-3]。亚洲是第二大风尘释放区，每年排放约1～8亿t，约占全球风尘排放量的25%。东亚地区释放的风尘约30%重新沉降在沙漠区或其周边，20%的风尘会在中国内陆沉降（如黄土高原等），剩下的50%则被输送到北太平洋甚至更远的地区[2]。风尘“源-汇”过程又敏感地响应于区域和全球气候-环境变化，在不同气候边界条件下，受多个物理、化学和生物过程的动态影响[1-3]。

西太平洋低纬度地区的菲律宾海是一个典型风尘汇聚区，受到国内外众多学者的高度关注[4]。前人已对菲律宾海第四纪沉积物中风尘物质贡献定量化、源区识别及其搬运动力开展了大量研究[5-9]，发现菲律宾海沉积物主要是火山和风尘的双端元混合产物，但对风尘源区究竟是亚洲东部沙漠地带还是亚洲中部沙漠地带以及其具体携带动力究竟是东亚冬季风还是西风带仍存在一些争议[10-11]。最近利用沉积捕获器获得的西菲律宾海季节性悬浮颗粒物Sr-Nd同位素数据，证实了风尘组分主要来自鄂尔多斯沙漠（约80%）且主要发生在冬季风爆发期（即冬—春季节）[12]。现代冬季西菲律宾海上空风尘物质的矿物和地球化学特征也表明鄂尔多斯沙漠是现代风尘的主要源区，而东亚冬季风是风尘搬运的重要动力[13]。此外，前人研究采用的站位大多位于菲律宾海吕宋岛附近海域（如班哈姆隆起），离吕宋岛距离较近（约250 km），接受了大量吕宋岛火山物质（50%以上，最高达90%）[6,9]，可能较大程度上反映了热带季风降雨、赤道太平洋ENSO相位或海平面升降的影响[14-18]。相反，菲律宾海帕里西维拉海盆的PV090102孔（距吕宋岛约1 200 km，火山灰成分只占10%～30%）的粒度和Sr-Nd同位素记录则更清晰记录了近2 Ma以来亚洲风尘沉积的演化历史，更好地反映了晚第四纪以来亚洲大陆干旱化历史和大气环流异常[19-20]。因此，菲律宾海远离吕宋岛的深海沉积物是重建构造尺度亚洲风尘演化历史的良好载体（图1）。

通过对比风尘颗粒及其母岩中较为特征的信息有助于深入理解亚洲风尘的源-汇过程。理想的物源示踪指标应具备以下特征，传输过程中不受风力分选的影响，沉积之后不受早期成岩改造，能代表风尘的主要成分信号[21-23]。黏土矿物作为深海沉积物的主要组成部分，其形态、含量、特征参数及组合类型不仅能够有效指示物源，而且可以记录源区气候和环境信息。因此，明确深海沉积物中黏土矿物的来源和搬运模式可为研究源区古环境演化和推测古大气环流模式提供关键证据，已在黄海[24]、南海[25]、菲律宾海[15]以及日本海[26]的相关研究中得到了广泛应用，取得了一系列重要成果。本文选取菲律宾海盆XT4孔沉积物开展了黏土矿物研究工作，并结合前人研究成果，在明确黏土矿物来源的基础上，探讨晚上新世以来亚洲风尘的演化历史及其可能的影响机制。

1 材料和方法

2018年在西菲律宾海盆南部使用重力采样器取得 XT4 孔岩芯（16°04.3′N、133°29.1′E）（图 1），长4.28 m，水深5 372 m，黄褐色致密状黏土，岩性均一，顶部含有少量铁锰结核，按2 cm间隔进行取样，并开展下一步测试分析工作。XT4孔位于太平洋碳酸盐补偿深度以下，样品中钙质生物基本都被溶解，难以用有孔虫壳体氧同位素来确定年代，故采用古地磁年代学进行了年代标定[27]。通过XT4孔的磁极性序列与国际地磁年代标尺进行精细对比，初步建立了XT4孔年代框架，该孔54 cm处对应于中更新世（0.781 Ma）的开始，第四纪底界位于322 cm（2.581 Ma），外插得到底部年龄大致为3.723 Ma。

沉积物样品前处理方法如下：称取约5g沉积物全样，置于离心管中，加入15%的H2O2水浴60 ℃加热约1 h去除有机质，再加入30 mL HCl（1M）水浴60 ℃去除碳酸盐，最终获得碎屑组分，后续开展粒度分析（214个样品，2 cm间隔）和黏土矿物分析（107个样品，4 cm间隔）。粒度分析在青岛海洋地质研究所海洋地质实验检测中心进行，所用设备为MS2000型激光粒度分析仪，测量范围为0.5～2 000 μm，重复测量的相对误差小于2%。黏土矿物分析测试在中国石油大学（华东）完成，所用分析仪器为荷兰帕纳科公司X’Pert PRO MPD，角度重现性：± 0.000 1°，角度精度：0.002 5°；测试精度：误差小于10%。利用离心法从碎屑组分中提取＜2 μm的黏土组分，制成自然风干定向片（N片）和乙二醇饱和片（EG 片，60 ℃，12 h），部分样品制成 490 ℃ 加热片（4 h）。依据不同黏土矿物在3种测试片（自然定向片、乙二醇饱和片和高温加热片）中衍射峰发生的不同变化对黏土矿物进行识别。蒙脱石（001）衍射峰在自然定向片中的范围是12～15Å，而在乙二醇饱和片中其（001）衍射峰可移至17Å左右，这是判断是否存在蒙脱石并进一步区分蒙脱石和绿泥石的主要依据。伊利石的10Å和5Å衍射峰在自然定向片中比较明显，经乙二醇饱和后10Å和5Å衍射峰的位置和强度变化很小。高岭石和绿泥石在7Å和3.5Å附近都存在特征峰，若3.5Å附近的特征峰表现为双峰，证实高岭石和绿泥石同时存在。经高温加热后，高岭石矿物晶格被破坏，变为非晶质物质，因而其衍射峰消失，绿泥石衍射峰也会变弱甚至消失。

根据Biscaye方法，利用Jade 6.5软件在乙二醇曲线上蒙脱石（17Å）、伊利石（10Å）、绿泥石（7Å）衍射峰峰面积各自乘以其强度因子1、4和2，开展黏土矿物相对含量计算。绿泥石和高岭石两者的含量比例则依据绿泥石（3.54Å）和高岭石（3.58Å）的衍射峰面积比值求得，最后将4种黏土矿物的特征峰面积校正为含量100%。此外，伊利石结晶度可由乙二醇饱和曲线上10Å衍射峰的半峰宽表示，其值越低表示伊利石矿物晶格有序度越高，指示源区水解作用较弱，气候寒冷、干燥。伊利石化学指数的计算是在乙二醇饱和曲线上根据5Å/10Å峰面积比确定的，比值小于0.5为富Fe-Mg伊利石，为物理风化结果；比值大于0.5为富Al伊利石，代表强烈的水解作用。

图1 研究区地理背景及XT4孔站位示意图中国北部沙漠（NCDs）包括：GD-古尔班通古特沙漠、ODSL-浑善达克沙地、HBSL-浑伦布尔沙地、HSL-科尔沁沙地；中国中西部沙漠（WCCDs）包括：TKD-塔克拉玛干沙漠、QD-柴达木沙漠、TD-腾格里沙漠、BJD-巴丹吉林沙漠；中国东部沙漠（ECDs）包括：HD-库布齐沙漠和MUD-毛乌素沙漠。Fig.1 Geographical background of the study area and location of Core XT4The northern Chinese deserts (NCDs, include GD, Gurbantunggut Desert; ODSL, Onqin Daga Sandy Land; HBSL, Hunlun Buir Sandy Land; and HSL,Horqin Sandy Land); western and central China deserts (WCCDs, e.g., TKD, Taklimakan Desert; QD, Qaidam Desert; TD, Tenger Desert; and BJD,Badain Jaran Desert); and East China deserts (ECDs, e.g., HD, Hobq Desert; and MUD, Mu Us Desert).

2 结果

XT4孔碎屑组分的粒度参数如图2所示，主要为粉砂（50.17%～81.49%，平均为66.47%），其次为黏土（17.65%～49.83%，平均为33.10%），大部分层位不含砂，少数层位含砂（0.1%～3.52%）。XT4孔碎屑组分的中值粒径变化范围为3.92 ～14.53 μm，平均值为7.31 μm。3.7 Ma以来XT4孔粒度参数随深度增加呈现出规律性变化，3.7～2.9 Ma期间中值粒径较粗，黏土含量较少，粉砂含量较多，砂含量主要出现在这一时期；2.9 Ma以来则相反，中值粒径逐步减小，黏土含量逐渐增加，粉砂含量逐渐减少，砂含量基本消失。

XT4孔的黏土矿物组成如图3所示，主要由蒙脱石、伊利石、高岭石和绿泥石组成，其中伊利石含量最高，为10%～58%，平均值为39%；蒙脱石次之，含量为11%～70%，平均值为28%；绿泥石和高岭石含量相当，含量分别为7%～55%和0～31%，平均值分别为18%和15%。伊利石结晶度为0.183°～0.482°Δ2θ，平均值为 0.272°Δ2θ，绝大部分都小于 0.4°Δ2θ，说明 XT4 孔伊利石结晶度很好；伊利石化学指数变化范围为0.15～0.68，平均值为0.35，为富Fe-Mg伊利石。蒙脱石含量从3.7 Ma开始逐步降低，直到1.5 Ma左右蒙脱石含量基本保持稳定；与之相反，伊利石含量在3.7～2.7 Ma期间逐步增加，在2.7～1.5 Ma期间伊利石含量快速增加，1.5 Ma之后基本稳定，仅伴随小幅波动变化。高岭石和绿泥石含量除少数几个层位波动剧烈外，整体变化幅度较小，且变化趋势基本一致。XT4孔底部的伊利石结晶度和化学指数波动较大，随后逐步稳定，且表现出长期减小的趋势。

图2 XT4孔的岩性和中值粒径垂向变化Fig.2 Vertical variations of lithology and median grain size of Core XT4

图3 XT4孔年代模型（a）和黏土矿物组成（b-e）及伊利石结晶度（f）、化学指数（g）垂向变化粗实线为局部加权回归拟合。Fig.3 Age model (a), vertical variations of clay minerals assemblages (b-e), illite crystallinity (f) and chemical index (g) of Core XT4Solid line is local weighted regression fitting.

3 讨论

3.1 菲律宾海盆XT4孔黏土矿物的物源判别

明确黏土矿物的来源及其搬运动力机制是进行古气候、古环境研究的重要前提。前人利用黏土矿物、元素地球化学和Sr-Nd同位素等多种示踪方法对菲律宾海盆邻近海域（奄美三角盆地、四国海盆和帕里西维拉海盆等）沉积物来源开展了综合研究工作，表明菲律宾海沉积物是亚洲大陆风尘和周围火山岛弧的双端元混合产物[11,19,28-31]。黏土矿物组成的三角图解显示XT4孔黏土矿物来源与邻近海域（如帕里西维拉海盆、四国海盆、西菲律宾海盆）类似，基本上符合亚洲大陆风尘与周边火山岛弧的二端元混合模型特征（图4），与前人结果一致。

3.1.1 蒙脱石

海洋沉积物中的蒙脱石来源主要有两种，一种是陆源碎屑，源岩在大陆温暖至半干旱气候条件下发生水解作用而成；另一种是海底基性火山物质经长时间海解作用而成。XT4孔蒙脱石含量平均为28%，最高可达70%，远高于长江、黄河、黄海及东海的蒙脱石含量[24]，且台湾河流中几乎不含蒙脱石[32]，因此，大陆河流输入物质（如长江、黄河、台湾等）不可能是XT4孔蒙脱石的主要来源。吕宋岛的蒙脱石来源于岛弧基性火山物质的蚀变，含量高达60%以上[33]，但是考虑吕宋岛河流物质受到自东向西的北赤道流阻隔，很难向XT4孔提供大量的蒙脱石。已有研究表明菲律宾海海底基性火山物质可风化形成自生蒙脱石[34]，前人研究利用扫描电镜观察也发现附近海域的蒙脱石形状不规则且边缘模糊，符合自生蒙脱石的特征[30]。综上所述，我们认为XT4孔中蒙脱石是周边火山岛弧物质的风化和蚀变产物。此外，XT4孔中蒙脱石与陆源成因的伊利石存在显著负相关（图5a），这也说明蒙脱石并非来源于亚洲大陆，而是主要来源于周边火山岛弧物质。

3.1.2 伊利石

伊利石是陆源火成岩、变质岩或沉积岩风化蚀变产物，形成于典型的干冷气候条件下。由于吕宋岛河流沉积物中伊利石含量很少[33]，因此，XT4孔中高含量的伊利石（平均为39%，最高为58%）不可能来自于包括吕宋岛在内的菲律宾群岛。XT4孔中伊利石化学指数变化范围为0.15～0.68，平均值为0.35，指示了干冷气候条件下的物理风化作用；伊利石结晶度变化范围为 0.183°～0.482°Δ2θ，平均为0.272°Δ2θ，与中国内陆黄土和古土壤伊利石结晶度接近（分别为 0.22°～0.33°Δ2θ和 0.22°～0.42°Δ2θ）[35]，说明研究区伊利石源区与中国内陆的黄土可能有关。另外，虽然伊利石含量在台湾河流沉积物中较高（＞65%）[32]，但是其伊利石结晶度值一般都小于0.23°Δ2θ，而伊利石的化学指数平均值为 0.49[32]，与XT4孔的结果存在显著差异。因此，台湾河流虽然每年向海输入大量富伊利石的沉积物，却不是XT4孔的伊利石主要源区。前人研究利用扫描电镜观察发现邻近海域的伊利石大多为不规则片状且磨圆度较好，进一步说明伊利石并不是自生成因[30]。综上所述，我们认为XT4孔沉积物的伊利石主要来自于亚洲大陆，风力吹扬为其主要的输运方式。

图4 XT4孔与潜在物源区的蒙脱石-伊利石-绿泥石三角图解[29]Fig.4 Ternary diagram of the major clay minerals smectite-illite-chlorite of Core XT4 and the potential source areas[29]

3.1.3 绿泥石

绿泥石多形成于高、中纬度区干冷的气候环境中[36]，深海沉积物中的绿泥石大多为陆源碎屑成因，但局部地区也有少量绿泥石是火山物质蚀变而成[37]。晚上新世以来XT4孔中绿泥石与蒙脱石变化相反（图3），可见绿泥石与蒙脱石来源不同。前已述及，伊利石主要源于亚洲大陆，因此我们认为绿泥石也主要来源于亚洲大陆。

3.1.4 高岭石

高岭石是酸性岩石在潮湿温暖气候下被强烈淋滤形成的，如吕宋岛、珠江、红河等的高岭石含量相对较丰富[25]。海洋的碱性、弱碱性条件不适合高岭石生成，因此，大洋中的高岭石一般都为陆源成因[38]。XT4孔沉积物中高岭石平均含量为15%，与亚洲大陆黄土中高岭石含量（6%～16%）相当，表明高岭石主要来自亚洲大陆风尘。

图5 XT4孔沉积物不同黏土矿物相关图解Fig.5 Correlation diagrams between clay minerals in Core XT4

3.2 菲律宾海盆XT4孔黏土矿物指示的古气候意义

不同类型的黏土矿物的形成条件具有明显差异，如中、基性火山物质在温暖且干湿交替的环境下容易形成蒙脱石；而低纬度地区炎热潮湿、化学风化较强的环境中易形成高岭石；而伊利石和绿泥石往往是干旱寒冷条件下物理风化的产物[37,39]。由于黏土矿物的含量是一个相对值，单个黏土矿物的含量变化很难直接用来指示古环境变化，而不同黏土组合及其比值可以有效消除黏土矿物之间的稀释效应，被广泛用于指示古气候变化[5,15]。伊利石结晶度和化学指数也是气候变化的典型指标，当气温升高、降雨增加时，化学风化增强，伊利石的结晶度变差，化学指数增加，反之当气候较为干燥寒冷时，伊利石结晶度较好，往往生成Fe-Mg伊利石[37]。

前文讨论可知XT4孔中黏土矿物主要由伊利石和蒙脱石组成，其中伊利石主要来自亚洲大陆风尘，而蒙脱石主要由周围海底火山灰蚀变后形成。因此，XT4孔的伊利石/蒙脱石比值代表了亚洲大陆风尘物质对于火山物质的相对变化，可用于指示菲律宾海盆亚洲风尘贡献的演化历史。晚上新世以来XT4孔的伊利石/蒙脱石比值表现出明显的阶段变化（图 6a）：3.7～2.7 Ma期间，伊利石/蒙脱石比值较低，且表现为逐步降低的趋势；2.7～ 0.9 Ma期间，伊利石/蒙脱石比值先快速变大，随后在高值大幅波动；0.9 Ma以来，伊利石/蒙脱石比值表现为逐步降低的趋势。XT4孔伊利石/蒙脱石比值所反映的亚洲风尘贡献变化历史与北太平洋ODP 885/886孔风尘通量[40]（图6d）（受年代框架误差的限制，存在一定的偏差）、日本海IODP 1422C孔沉积物钾含量[41]（图6e）的变化趋势大体一致，均表现为在约2.7 Ma以来风尘通量或贡献逐步升高，与该时期北半球高纬度地区首次出现大规模冰盖的时间节点一致（图6g）。这些沉积记录共同反映了北半球高纬度地区冰川扩张，在加剧亚洲中纬度地区干旱化的同时，也通过冰川剥蚀、寒冻风化、河谷深切等高山剥蚀过程促进了细颗粒物质的产生，随后通过风尘释放将这些信号进一步传送到下风带地区（如日本海、菲律宾海、北太平洋等）。

黄土高原的黄土-红黏土风尘沉积序列是反映东亚季风系统形成演变的良好地质记录，黄土粒度和磁化率被广泛用于指示东亚冬季风和夏季风的演化特征[42-43]。甘肃灵台剖面的黄土粒度记录[44]（图6b）和黄土堆积通量记录[45]均指示东亚冬季风强度在北半球永久性冰盖出现后（～2.7 Ma）明显增强，随后保持高频波动趋势，与XT4孔伊利石/蒙脱石比值的变化趋势存在一定的差异。有意思的是，XT4孔伊利石/蒙脱石比值在各阶段内部（即3.7～2.7、2.7～0.9、0.9 Ma以来）的变化趋势，与灵台黄土磁化率[44]所指示的东亚夏季风变化表现为一定的反相关系（图6c）。也就是说，灵台黄土磁化率所记录的东亚夏季风越强，XT4孔的伊利石/蒙脱石比值（代表亚洲风尘贡献程度）越低，反之则相反。事实上，粉尘释放是一个复杂的风沙物理过程，其是否发生及发生的强度由气流状况和下垫面状况共同决定[12,46-47]。风场为颗粒的侵蚀、搬运提供了动力来源，又称动力因子，而下垫面状况决定了地表抵抗风力侵蚀的能力大小，称为抗蚀因子。黄土高原及周边沙漠处于干旱-半干旱地区，年均降雨量200～400 mm，夏季降雨（7—9月）是主要水汽来源，可占全年降水的60%～80%。因此，东亚冬季风所携带风尘通量可能也受控于下垫面状况，也就是说细粒级物质的供应量可能也是控制东亚冬季风影响区风尘释放的重要环境因子之一[21]。据此，我们推断晚上新世以来，除了全球变冷的大背景外，东亚夏季风可能通过降雨变化深刻影响着中国干旱-半干旱地区下垫面状况，进而决定了沙尘源丰富度。

图6 XT4孔伊利石/蒙脱石比值（a）与灵台黄土东亚冬季风指标[44]（b）、灵台黄土东亚夏季风指标[44]（c）、北太平洋ODP 885/886孔风尘通量[40]（d）、日本海IODP 1422C孔沉积物钾含量[41]（e）、热带太平洋经向温度梯度变化[48]（f）和深海底栖有孔虫氧同位素[49]（g）对比粗实线为局部加权回归拟合。Fig.6 Comparison of (a) illite/smectite ratio from Core XT4 in the Western Philippine Sea (this study), (b) the mean grain size of quartz particles (MGSQ) from Lingtai loess as a proxy for East Asia Winter Monsoon[44], (c) East Asia Summer Monsoon index based on magnetic susceptibility from Lingtai loess[44], (d) mass accumulation rate of eolian dust at Core ODP 885/886[40], (e) the K (wt.%) record at site IODP U1422 in Japan Sea[41], (f) Equatorial Pacific Ocean zonal SST gradient[48], and (g) the stacked global benthic δ18O record of LR04[49]solid line is local weighted regression fitting.

现代观测资料和数值模拟结果均表明在拉尼娜和厄尔尼诺现象之间交替发生的厄尔尼诺-南方涛动循环（ENSO）主导着东亚降水的年际变化[50-52]。在拉尼娜（厄尔尼诺）期间，赤道太平洋Walker环流增强（减弱），西北太平洋副热带高压向北（南）迁移，夏季风降雨锋面北侵（南撤），中国干旱-半干旱地区的降雨量增加（减少）[53-54]。进一步研究则表明在更长时间尺度上（全新世、轨道尺度、构造尺度），类ENSO系统可能起到类似于热带气候变化“放大器”的作用，将热带气候系统较小的扰动放大，通过海洋和大气环流的传输进而影响东亚地区降雨强弱和空间变化[15,55]。晚上新世以来太平洋经向温度梯度（图6f）逐步变大，表明太平洋类ENSO系统相位逐步趋向于拉尼娜相[48,56]，通过控制热带辐合带（ITCZ）的南北迁移，可能影响着中国干旱-半干旱地区的降雨强弱和空间分布。

4 结论

菲律宾海盆XT4孔的黏土矿物以伊利石为主，平均含量为39%，蒙脱石含量次之，平均为28%，绿泥石平均含量为18%，高岭石平均含量为15%。XT4孔的伊利石主要来自亚洲大陆风尘，而蒙脱石主要由海底火山灰蚀变后形成，因此该孔伊利石/蒙脱石比值可用于指示菲律宾海盆亚洲风尘的演化历史。晚上新世以来XT4孔的伊利石/蒙脱石比值表现出的三阶段变化，反映了北半球高纬度地区冰川的扩张导致全球持续变冷，亚洲中纬度地区干旱化加剧，通过风尘释放将这些信号进一步传送到下风带地区。除了全球变冷的大背景外，细粒级物质供应量可能也是控制东亚冬季风影响区风尘释放的重要环境因子之一。进一步分析表明，晚上新世以来赤道太平洋经向温度梯度的变化，通过海洋和大气环流的传输控制着东亚地区降雨的空间变化，进而影响着该地区的风尘释放强弱。