70 万年以来西太平洋暖池区硅质生产力记录及其气候效应

唐正，熊志方,2，贾奇，秦秉斌，李铁刚,2

1. 自然资源部海洋地质与成矿作用重点实验室，青岛 266061

2. 崂山实验室海洋地质过程与功能实验室，青岛 266237

西太平洋暖池（WPWP）是全球海洋和气候系统的重要组分，也是大气圈热量和水汽的主要来源，因此可以驱动诸如厄尔尼诺-南方涛动（ENSO）的全球尺度气候事件。而西太平洋暖池的生产力水平被认为在调节第四纪全球大气CO2分压的变化上发挥着重要作用[1-3]。之前的研究认为，风尘和陆源有机质输入是WPWP 冰期-间冰期旋回尺度生物地球化学循环的首要控制因素[4-5]。也有研究发现，南大洋中层水（AAIW）可以将深层上涌的高营养水体通过“大洋隧道”过程水平传输至热带太平洋上层水体[6-7]，从而影响其冰期旋回中的硅质生产力变化。然而冰期旋回中WPWP 生产力的控制因素及其主次，目前依然存在争议。

初级生产者（硅藻）对溶解态硅酸（H4SiO4）的消耗是有效碳封存的首要先决条件[8]。硅藻生活在透光层，并利用硅酸合成生源蛋白石。生源蛋白石指示的生产力不仅受限于溶解态硅酸的供应[9-10]，铁的供应也可以限制单独硅藻细胞的物理状态和生产速率[11]。WPWP 沉积物中生源蛋白石的含量，不仅可以记录该海区的生产力水平和碳封存状态，也能记录陆源有机质、南源水体和风尘铁输入对WPWP的影响。本文对取自西菲律宾海的MD06-3047 站柱状沉积物开展了生源蛋白石含量分析，探讨了700 ka以来WPWP 的硅质生产力变化及其控制因素。

1 区域地质背景

西菲律宾海是西太平洋最大的边缘洋盆。沿着西太平洋的菲律宾沿岸，北赤道流（NEC）在15°N分流为向北的黑潮和向南的棉兰老流（MC）[12-14]。NEC 分叉点南北移动的年际变化，受ENSO 过程的强烈影响，在厄尔尼诺发生的年份分叉点位置位于最北端[14]。来自南大洋的AAIW 可以通过南赤道流（SEC）和新几内亚沿岸潜流进入到赤道中层水（EqPIW），并成为该水体的主要来源[7,15]。东亚夏季风和冬季风控制了研究区的气候，导致了吕宋火山岩和亚洲风尘向该地区的陆源输入。另外，研究区的表层海水温度（SST）几乎是全球最高，其季节变化最小，且年均SST 基本大于28°C[16]。研究区因较低的表层海水营养盐浓度，整体上呈现出较低的初级生产力（叶绿素a 浓度 0.04～0.07 mg/m3），而较厚的混合层也限制了富含营养物质的深层水体上涌[17-18]。

2 材料与方法

2.1 样品来源

重力柱状岩芯 MD06-3047，总长度890 cm，由中-法合作航次MARCO POLO 2/IMAGES XIV 于2006 年在西菲律宾洋盆的本哈姆海台（17°00.44'N、124°47.93'E，水深约2 510 m，图1）获取。该岩芯样品主要由黄褐色粉砂质黏土组成，未见明显生物扰动。岩芯位置位于吕宋岛以东240 km，该海区主要受亚热带东亚季风气候影响[19]。

图1 MD06-3047 岩芯站位信息图其中，NEC：北赤道流，KC：黑潮，MC：棉兰老流，WPS：西菲律宾海，EAWM：东亚冬季风，EASM：东亚夏季风。Fig.1 Location of core MD06-3047NEC: North Equatorial Current, KC: Kuroshio Current, MC: Mindanao Current, WPS: western Philippine Sea, EAWM: East Asian winter monsoon, EASM:East Asian summer monsoon.

2.2 分析方法

将岩芯按照2 cm 间隔分为445 个样品用于分析测试。样品在50°C 下被烘干，接着用5%的H2O2进行处理，然后用63 µm 孔径的筛子在去离子水下进行冲洗，以获取有孔虫壳体样品。每个样品挑选粒径大于150 µm 的浮游和底栖有孔虫壳体共计300 个起进行属种鉴定和计数，并进行了底栖有孔虫表生种Cibicidoides wuellerstorfi有孔虫稳定氧同位素测试和温跃层深度（DOT）的函数转换分析[3]。

生源蛋白石含量采用湿碱消解法测定[20]，该方法的平均误差为0.2%～1.0%，在自然资源部海洋地质与成矿作用重点实验室由钼酸蓝分光光度计进行测试。具体分析测试方法为：冻干样品研磨至200 目以下，准确称取约130～140 mg 样品粉末，置于 50 mL 离心管中。将 5 mL 10% 过氧化氢溶液加入离心管中，震荡并静置 30 min 以除去有机质；再在离心管中加入 5 mL 1∶9 盐酸，震荡反应 30 min以去除碳酸盐。随后，加入 40 mL 去离子水并以3 000 rpm 离心 10 min，滤出上清液，并重复该过程 3次后将样品烘干。在烘干后的样品中加入 40 mL 2 mol/L 的 Na2CO3溶液，充分混合后放入 85 ℃ 恒温水浴箱中提取生物硅。每小时从水浴箱中取出样品离心 15 min，取 125 μL 上清液加入钼酸铵溶液，配置还原剂使其显色。每次取样后用力摇晃试管使固体重新悬浮后，继续水浴加热提取，共提取 6次；该过程快速完成，以减少可溶硅在固体表面的不可逆损失。通过钼酸盐蓝光分光光度法分步测量提取物中的溶解硅，并根据硅的质量百分比随时间变化得出线性回归方程，其截距为沉积物中生源硅的含量。沉积物中的生源蛋白石含量按公式（Opal%= 2.4×BSi%）计算。

3 结果与讨论

3.1 年龄框架重建

MD06-3047 岩芯的年龄框架通过将该岩芯底栖有孔虫表生种C. wuellerstorfi的δ18O 曲线[3]和LR04 δ18O 标准曲线[21]进行图形比对来获得，并且考虑了浮游有孔虫Globigerinoides ruber(pink)的末现面。G. ruber(pink)的末现面出现在岩芯160 cm 处，对应年龄为120 kaBP[22]。该岩芯地层覆盖了海洋氧同位素期次MIS 1-17 期（图2），记录了WPWP 近700 ka的古海洋变化历史。

图2 MD06-3047 岩芯生源蛋白石含量特征A：LR04 氧同位素标准曲线[21]，B：MD06-3047 岩芯底栖有孔虫C. wuellerstorfi 氧同位素[3]，C：MD06-3047 岩芯生源蛋白石含量，D：MD06-3047 岩芯总有机碳含量[23]。图中灰色条带为间冰期。Fig.2 Characteristics of biogenic opal content in core MD06-3047A: LR04 oxygen isotope standard curve[21], B: benthic foraminifera[3], C: oxygen isotope of benthic foraminifera C. wuellerstorfi, D: total organic carbon content[23]. The gray bars in this figure are interglacial periods.

3.2 冰期旋回中的西太平洋暖池生源蛋白石含量特征

MD06-3047 岩芯的生源蛋白石含量在700 ka以来平均值为3.10%，最大值8.46%，最小值1.40%。在冰期-间冰期旋回中，其阶段性低值往往出现在间冰期；而在大多数间冰期的中后期到紧接着的冰期，对应着有孔虫氧同位素记录的全球冰量增加[3,21]，生源蛋白石含量往往呈现出阶段性上升的趋势（图2C）。其最大值出现在MIS 5.4 期，基本对应着MIS 5.5 期到MIS 5.4 期全球冰体积的迅速增加。在MIS 11 期以前，其冰期旋回尺度的变化幅度相对较小，而在MIS 11 期之后显著增大。与同为生产力指标的总有机碳含量（TOC）比较（图2D）[23]，在冰期旋回中整体变化趋势一致，即大致呈现出了冰期高、冰消期和间冰期低的特点。

3.3 西太平洋暖池硅质生产力控制因素

可能影响WPWP 硅质生产力的因素包括：风尘带来的Si 和Fe[2,23]、陆架物理剥蚀和硅酸盐风化导致的陆源有机质输入[24-25]、与类ENSO 过程相关的温跃层动力机制[26-27]和南大洋“硅溢漏”[28]。我们认为，表层海水硅藻生产力的勃发需要3 个条件：① 充足的硅酸盐供应；② 可供生物利用的铁的供给；③ 营养物质在表层水体较小的纵向空间和充分的滞留时间。根据上述研究提出的可能影响因素，结合硅藻生长所需条件，我们假设：亚洲陆源输入和南大洋为研究区供应硅酸；亚洲风尘Fe 的输入为水体提供Fe 供给；DOT 的深浅改变了混合层的深度，从而影响了硅酸和Fe 在表层水体的纵向空间和滞留时间。

通过与本岩芯风尘通量和溶解Fe 含量的比较（图3A，B）[23]，我们发现生源蛋白石和这两项风尘指标基本都在冰期表现为高值，而冰消期时蛋白石的低值也往往对应着风尘的低值。这说明风尘通过Fe 的输入对WPWP 生产力的影响符合我们的假设。而通过与本岩芯DOT 记录的比对，发现生源蛋白石在冰期时的高值往往对应DOT 较浅的时期；其低值期，DOT 也往往较深（图3D）。这和我们的假设一致，即较浅的温跃层可以减少营养物质和Fe 的垂向空间，增强其在透光层的滞留，为硅藻的生成创造时间窗口。而DOT 较深时，往往是北赤道流输送至WPWP 核心区的寡营养水体逐渐聚集、堆积，不利于硅藻的生长。为了验证南大洋“硅溢漏”对WPWP 硅质生产力的影响，我们比对了南大洋南极带上升流区的生源蛋白石记录和南半球西风带（SWW）记录（图3E，F）[29-30]，发现其最显著特征是在每个冰消期呈现出峰值，意味着冰消期SWW的南移和增强，导致了南极带风驱上升流的增强、富硅深层水体进入南极中层水和亚南极模态水并进一步向低纬的上层水体输送[29]。然而这一冰消期“硅溢漏”特征在WPWP 的生源蛋白石记录中并未得到响应。这说明冰消期南大洋“硅溢漏”未能影响到WPWP 的硅质生产力。另外，来自东菲律宾海的硅藻席硅同位素的研究也表明，末次盛冰期南大洋“硅溢漏”并未显著影响到研究区[18]。本岩芯生源蛋白石含量并未在末次盛冰期出现峰值，可能与年龄框架的不确定性、生源蛋白石的保存状况有关（图3C）。而跟海平面记录的比对（图4B）[31]，则进一步支持了我们的假设，即低海平面时裸露的东吕宋陆架的硅酸盐物理剥蚀和化学风化可能为研究区表层海水提供了硅酸和有机质，导致了生源蛋白石的增长，这与前人在这一海区的研究一致[25,32]。

图3 MD06-3047 岩芯生源蛋白石含量与潜在控制因素比对MD06-3047 相关参数：A：风尘通量[23]，B：溶解铁含量[23]，C：生源蛋白石含量，D：温跃层深度[3]；E：南大洋南极带ANT29/P1-03 岩芯生源蛋白石含量[29]；F：南极EDC 冰芯积累率[30]。图中灰色条带为间冰期，绿色半透明条带为冰消期。Fig.3 Comparison of bornite opal concentration and potential controls in core MD06-3047Proxies in MD06-3047: A: dust flux[23], B: dissolved iron content[23], C: biogenic opal concentration, D: thermocline depth[3]. E: biogenic opal content of core ANT29/P1-03 in the Antarctic Zone of the Southern Ocean[29], and F: accumulation rate in Antarctic Epica Dome C ice core[30]. The gray bars are interglacial periods, and the green translucent bars are the deglacial periods.

图4 MD06-3047 岩芯相关参数与全球气候指标比对MD06-3047 相关参数：A：风尘通量[23]，B：全球海平面变化[31]，C：生源蛋白石含量，D：总有机碳含量[23]，E：温跃层深度[3]；F：南极冰芯Epica Dome C 大气CO2 浓度[1]。图中灰色条带为间冰期，绿色半透明条带为冰消期。Fig.4 Comparison of MD06-3047 proxies with global climate proxiesMD06-3047 proxies: A: dust flux[23], B: global sea-level changes[31], C: biogenic opal concentration, D: total organic carbon concentration[23], E: depth of thermocline[3]; F: atmospheric CO2 concentration in the Epica Dome C ice core[1]. The gray bars are interglacial periods, the green translucent bars are the deglacial periods.

3.4 冰期WPWP 硅质生产力的增强与pCO2atm 的降低

冰期旋回中WPWP 的生源蛋白石含量大多在冰期中处于高值，这与TOC 的含量高值基本对应（图4C，D）。而且WPWP 硅质生产力的高值，往往对应大气CO2低值；其生产力增加的过程，则往往对应pCO2atm迅速降低的过程。WPWP 的硅质生产力，可能在冰期旋回中影响了大气CO2的浓度[4-5]。在第四纪的冰期，同时也是低海平面期，广泛存在的热带火山弧，导致了此时期可侵蚀沉积物的面积显著增大，并使暴露的未固结碎屑遭受侵蚀与风化（图4B）[25,32]。而从裸露的东吕宋陆架入海的巨量淡水，为研究区输入了大量陆源硅，并使生产力显著提高[25,32]。低海平面时期从热带火山弧向深海的硅质有效搬运，可能调节了WPWP 的生产力、碳埋藏和大气CO2浓度[24,32-33]。而亚洲风尘，则是WPWP的铁的首要来源，刺激了这些海区的生产力水平（图4A）[2,4,18,23]。冰期时热带火山弧附近裸露的陆架沉积物的物理剥蚀、硅酸盐风化以及风尘输入的增加，可使WPWP 西部边缘海的有机碳埋藏效率提高，有可能降低了冰期时的大气CO2分压（图4）[1-2,24-25,32]。而在这一过程中，较浅的DOT 水平使表层海水中营养物质的垂向空间变小、滞留时间增强，则进一步提高了这一海区的生产力水平和碳埋藏效率（图4E）。

3.5 冰消期南大洋“硅溢漏”在WPWP 的失效与硅质生产力的降低

紧接着冰期的冰消期，南大洋西风带和风驱上升流增强，南大洋深部的大量硅酸进入表层，并伴随南大洋中层水生成的增强向北水平输送，进入低纬大洋的上层水体[28-30]。然而南大洋“硅溢漏”并未导致WPWP 硅质生产力的提高（图4C），我们认为这与冰消期这一海区Fe 的输入水平有关。在冰消期，海平面迅速上升，导致东吕宋陆架重新被海水覆盖，显著降低了可侵蚀沉积物的面积，向研究区的硅质输入随之减少；而亚洲风尘输入也明显降低，从而降低了对这一海区的Fe 的输入（图3B，图4）。即使南大洋溢漏的硅酸到达WPWP 上层水体，此时亚洲风尘带来的Fe 含量的骤降，将使硅藻生长速率无法有效提高，从而降低了这一海区的硅质生产力，使南大洋“硅溢漏”失效。而冰消期的温跃层深度均处于较深水平（图4E），加厚的混合层使表层水体中的营养物质纵向空间更大、滞留时间更短，而且是自东向西的寡营养上层水体堆积，这使冰消期WPWP 的硅质生产力进一步降低。冰消期降低的WPWP 硅质生产力，将降低WPWP 的碳埋藏效率，并影响大气CO2浓度（图4F）。

4 结论

（1）700 ka 以来西太平洋暖池的硅质生产力变化呈现显著的冰期-间冰期旋回，基本在冰期较高，而间冰期较低。其主要控制因素可能是海平面变化引起的东吕宋陆架沉积物风化以及亚洲风尘输入，改变了该海区的陆源硅酸和铁的输入。温跃层深度的改变则可能强化了这一过程。南大洋中层水的“硅溢漏”可能无法显著影响WPWP 的生产力。

（2）冰期时的低海平面，导致热带火山弧附近裸露的陆架沉积物的物理剥蚀和硅酸盐风化，通过淡水输入为WPWP 提供了更多的硅酸；增强的风尘输入则为WPWP 提供了更多的Fe；冰期时较浅的DOT 使表层海水的营养物质垂向空间变小、滞留时间增多。这些因素使冰期的WPWP 生产力增高，有可能降低大气CO2分压。

（3）冰消期时迅速升高的海平面、亚洲风尘输入的减少以及温跃层变深，都使此时可能进入到WPWP 的南大洋“硅溢漏”失效，并使其硅质生产力显著降低。