浮游有孔虫表层水种不同形态类型壳体的Sr/Ca 记录

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

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

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

浮游有孔虫壳体元素组成作为重要的海洋环境参数替代性指标，在古海洋学研究中发挥着重要的作用。随着浮游有孔虫壳体元素比值测试技术的快速发展，通过对浮游有孔虫壳体元素比值的分析，取得了一系列重要成果[1-3]。其中，Globigerinoides ruber（白色）和Trilobatus sacculifer壳体作为定量重建地质历史时期表层海水温度和盐度的重要信息载体，得到了广泛的研究和应用[1,4-5]。然而，在利用G. ruber和T. sacculifer壳体进行测试时，通常忽略G. ruber和T. sacculifer不同形态类型的壳体对重建结果的可能影响。

大部分海洋沉积物样品中G. ruber壳体存在两种形态类型，分别为狭义种（sensu stricto, s.s.）和广义种（sensu lato, s.l.）。根据分类学标准，G. rubers.s.壳体的主要形态为：一个近球形房室对称地生长于原有结构上，而形成具有高角度拱形的较大口孔；而G. rubers.l. 壳体具有相对紧凑的结构特征：一个扁平的房室不对称地生长于原有结构之上，从而形成具有中等角度的拱形和相对较小的口孔[6]。研究发现浮游有孔虫G.ruber壳体的两种形态类型具有不同的稳定同位素组成，且G. rubers.s.的生活水深浅于G. rubers.l.，因此两种形态类型的G. ruber壳体记录的海洋环境信息可能存在差异[4,7-8]。此外，也有研究发现热带和亚热带海区G. rubers.s.和G.rubers.l.壳体的元素比值也存在差异[4,9]。

T. sacculifer作为另外一个重要的浮游有孔虫混合层水种，其壳体在古海洋学研究中也得到了广泛的应用[10-11]。尽管从分类学角度来看，T. sacculifer的分类更加复杂，但其在形态方面的主要区分依据为是否具有最后一个似袋状房室[12]。因此，T.sacculifer从形态学上可以分为T. sacculifer（with sac）和T. sacculifer（without sac）。Elderfield 等[11]和Anand 等[13]通过对取自大西洋的沉积岩芯以及沉积物捕获器样品的分析，发现T. sacculifer（with sac）和T. sacculifer（without sac）壳体的元素比值同样存在一定的差异。

已有研究表明有孔虫壳体的Sr/Ca 比值在第四纪以来存在着明显的冰期-间冰期变化特征，可能是指示海水Sr/Ca 水平的潜在替代性指标[14-15]。此外，有孔虫壳体Sr/Ca 比值的变化可能可以指示第四纪冰期旋回中的海平面变化[16]。也有研究发现浮游有孔虫壳体的Sr/Ca 受海水温度和盐度等因素的影响[17]。因此，有孔虫壳体Sr/Ca 是古海洋学研究的潜在指标之一。本文通过西菲律宾海MD06-3047B岩芯中G. rubers.s.和G. rubers.l.以及T. sacculifer（with sac）和T. sacculifer（without sac）壳体Sr/Ca 比值的测试分析，探讨它们之间是否存在着显著性差异，并分析不同形态类型壳体Sr/Ca 比的影响因素，为未来利用两个浮游有孔虫表层水种在该区域开展古海洋学工作提供借鉴。

1 材料与方法

1.1 研究材料

MD06-3047B 岩芯（17°00.44′N、124°47.93′E）位于吕宋岛以东约240 km 的西菲律宾海本哈姆高原（图1a），水深2 510 m。该沉积岩芯沉积连续, 没有发现明显的沉积间断以及浊流沉积层，沉积柱状样主要由黄色粉砂质泥组成。根据前人研究，西菲律宾海现代碳酸盐溶跃面深度约为3 400 m[18]，MD06-3047B 孔位于海区溶跃面深度之上，因此该沉积物岩芯中有孔虫保存程度较好[19]。在本次研究中，我们选取钻孔岩芯上部60 cm，按4 cm 间隔取样，取得15 个层位的样品。每个层位样品分别挑选30～50 枚粒径范围为250～300 μm 的G. rubers.s.、G. rubers.l.、T. sacculifer（with sac）和T. sacculifer（without sac）壳体（图2）。并对这59 个有孔虫样品（G. rubers.s.有一层位缺失）进行Sr/Ca 比值测试。尽管G. ruber和T. sacculifer的生活水深存在差异，但两者的平均钙化深度均位于混合层内[7,13]，该层内海水温度和盐度随深度的变化较小（图1b）。

图1 MD06-3047B 孔站位信息a：年均表层海水温度信息，b：MD06-3047B 孔附近的垂向温度与盐度信息。数据引自WOA13[20]。Fig.1 The location of Site MD06-3047Ba: The annual mean sea surface temperature, b: profiles of water temperature and salinity around Site MD06-3047B.The data are from World Ocean Atlas 2013[20].

图2 G. ruber 和T. sacculifer 不同形态类型壳体的显微照片Fig.2 Microphotographs of G. ruber and T. sacculifer morphotypes

1.2 壳体元素测试

有孔虫壳体的清洗主要依据Barker 等[21]的方法。首先，将有孔虫壳体在显微镜下压碎，保证每个房室均被打开；分别对每一个有孔虫样品用去离子水进行超声清洗5 次、乙醇（优级纯）超声清洗2次、去离子水清洗2 次，用以去除黏土等；用加热的H2O2缓冲溶液进行氧化处理，并用去离子水进行清洗，以去除有机质等；在镜下剔出非有孔虫壳体碎片的杂物（暗色矿物、絮尘等）；转移并进行淋洗保存等待上机测试。元素测试在中国科学院海洋研究所电感耦合等离子发射光谱仪（ICP-OES）上进行。通过对标准溶液（Sr/Ca=1.20 mmol/mol）进行45 次重复测试分析，得到Sr/Ca 测试的标准偏差为1.1% (1σ)。

1.3 统计分析

为了从统计学角度分析G.ruber和T. sacculifer不同形态类型壳体的Sr/Ca 比值差异，我们依照Antonarakou 等[8]的方法，对MD06-3047B 岩芯中上述4 类有孔虫壳体的Sr/Ca 测试结果进行韦尔奇检验。首先，假设相比较的两组数据均值结果相当，如果G.rubers.s.和G. rubers.l.或T. sacculifer（with sac）和T. sacculifer（without sac）壳体的Sr/Ca 结果相当，即接受虚假设（H=H0），说明G.ruber或T.sacculifer不同形态类型的Sr/Ca 比值的差异不大；相反，如果对比结果存在显著差异，即拒绝虚假设（H=Ha），说明不同形态类型壳体的Sr/Ca 存在显著差异。

2 壳体Sr/Ca 结果

MD06-3047B 孔的年龄框架由Jia 等建立[19]，主要依据全球大洋底栖有孔虫氧同位素堆叠曲线[22]，并辅以粉红色G. ruber末现面（～120 ka）作为参考点[23]而确立。本次研究的样品时间跨度约48 ka，覆盖了MIS 3-1。图3 所示为MD06-3047B 孔MIS 3 期以来的G.rubers.s. 和G. rubers.l. 以及T.sacculifer（with sac）和T. sacculifer（without sac）壳体的Sr/Ca 比值。48 ka 以来G.rubers.s.和G. rubers.l.壳体的Sr/Ca 整体上具有相同的变化趋势，其差值变化范围为-0.006～0.022 mmol/mol，平均差值约0.006 mmol/mol。G. ruber（白色）的两种形态类型壳体的Sr/Ca 比值并没有表现出明显的阶段性高低变化规律，但G.rubers.l.壳体Sr/Ca 波动幅度相对较大。T. sacculifer（with sac）和T. sacculifer（without sac）壳体Sr/Ca 比值存在差异，整体上T. sacculifer（with sac）壳体的Sr/Ca 比值相对较高，两者的差值变化范围为-0.008～0.034 mmol/mol，平均约0.017mmol/mol。

3 讨论

3.1 G.ruber 和T. sacculifer 不同形态类型壳体Sr/Ca差异

如图3 所示，G.rubers.s.和G. rubers.l.壳体Sr/Ca变化趋势较为一致，平均差值仅约0.006 mmol/mol，小于Sr/ Ca 比值的测试误差（± 0.011 mmol/mol）。T.sacculifer（with sac）和T. sacculifer（without sac）壳Sr/Ca 平均差值为约0.017 mmol/mol，大于测试误差（± 0.011 mmol/mol）。同时，统计学韦尔奇检验结果也显示G.rubers.s.和G. rubers.l. 壳体Sr/Ca 平均值差异结果不具有显著差异；而T. sacculifer（with sac）和T. sacculifer（without sac）壳Sr/Ca 平均值差异显著（p＜0.05， 表1） 。 综上， 我们判断在MD06-3047B 站位的附近海区G. ruber不同形态类型壳体Sr/ Ca 比值的差异较小；而T. sacculifer不同形态类型壳体Sr/Ca 的差异较大。因此，在西菲律宾海区对浮游有孔虫表层水种Sr/Ca 比值进行测试分析时，如在样品量不足的情况下，可以选择G. ruber的不同形态类型壳体，但需要尽量选择T. sacculifer（with sac）或T. sacculifer（without sac）的单一形态类型壳体，以免造成结果偏差。

表1 MD06-3047B 孔G. ruber s.s.和G. ruber s.l.以及T. sacculifer（with sac）和T. sacculifer（without sac）壳体Sr/Ca 平均值（mmol/mol）以及韦尔奇检验结果（p＜0.05）Table 1 Mean shell Sr/Ca (mmol/mol) of G. ruber s.s., G. ruber s.l., T. sacculifer (with sac), and T. sacculifer (without sac) of Core MD06-3047B with the results of the Welch’s t-test at p＜0.05 level

3.2 不同形态类型壳体Sr/Ca 差异的影响因素

浮游有孔虫壳体Sr/Ca 主要受到海水Sr/Ca、海水温度、盐度和溶解作用等因素的影响[17,24-27]。根据前人关于海水Sr/Ca 的冰期-间冰期变化特征的研究可知，海水Sr/Ca 在冰期时高，而在间冰期时低[28-29]。这一特征与我们的结果并不完全一致（图4），特别是MIS 3 期的Sr/Ca 整体低于MIS 1 期。此外，不同形态类型壳体的Sr/Ca 变化并没有表现出完全一致的变化，也说明其他因素在其中发挥作用。因此，海水Sr/Ca 可能不是影响研究区浮游有孔虫壳体Sr/Ca 变化特征的唯一因素。有研究表明，G.ruber壳体Sr/Ca 可能受到海水温度和盐度的影响，而受pH 的影响较小[26]，随着海水温度和盐度的升高，G.ruber壳体Sr/Ca 呈增大趋势。T. sacculifer壳体Sr/Ca 同样受到温度和盐度的影响[27]，海水温度升高，T. sacculifer壳体Sr/Ca 越大，而盐度越高，Sr/Ca越小。

图4 MD06-3047B 孔G. ruber 和T. sacculifer 不同形态类型壳体Sr/Ca 与表层海水温度和盐度记录[19]的对比MIS 为氧同位素期次，T-I 为末次冰消期。Fig.4 Comparison of shell Sr/Ca of G. ruber and T. sacculifer morphotypes with sea surface temperature and salinity[19] from Core MD06-3047BMIS: marine isotope stage, T-I: the last deglacial stage.

如图3 和图4 所示，同一粒径范围下，MIS 3 期以来G. ruber两种形态类型的壳体Sr/ Ca 十分相似，说明两者的变化可能受到相同影响因素的控制。而T. sacculifer（with sac）和T. sacculifer（without sac）壳体Sr/Ca 存在显著差异（本文3.1 节），可能指示二者受不同因素的影响。进一步将G.rubers.s.、G. rubers.l.、T. sacculifer（with sac）和T. sacculifer（without sac）壳体Sr/Ca 进行对比可发现（图4），G.rubers.s.、G. rubers.l.和T. sacculifer（without sac）壳体Sr/Ca 呈现较为一致的变化趋势，并两两进行线性相关分析，发现三组记录之间具有较好的相关性，因此这两个种的3 种形态类型壳体Sr/Ca 的记录可能受到相同因素的影响。为方便分析，将这三组记录进行堆叠平均（Sr/Castack），并与同站位48 ka以来的表层海水温度和盐度等古海洋学记录[19]进行对比和线性相关分析。结果显示Sr/Castack与表层海水温度呈现线性正相关（图5a），从整体趋势上，G.rubers.s.、G. rubers.l.和T. sacculifer（without sac）的Sr/Ca 的增大对应表层海水温度的升高（图4）。而Sr/Castack与表层海水盐度替代性指标（δ18Osw-ice）之间不具有明显的相关性（图5b），并且G.rubers.s.、G. rubers.l.和T. sacculifer（without sac）壳体Sr/Ca 与表层海水盐度的变化趋势也存在较大差异（图4）。因此，研究区G.rubers.s.、G. rubers.l.和T.sacculifer（without sac）壳体的Sr/Ca 变化可能主要受到表层海水温度的影响，而受到盐度的影响较小。其中，Sr/Ca 的高值并未完全出现在MIS 1 期，而整体出现在末次冰消期，即MIS 2 期向MIS 1 期的过渡阶段（图4）。这可能是由于太平洋区域存在显著的末次冰消期表层海水温度显著增暖的特征[30]。此外，MD06-3047B 孔年龄框架由底栖有孔虫氧同位素建立，而根据前人工作，研究区SST 在变化特征上超前于底栖有孔虫氧同位素的变化[31]。

图5 MD06-3047B 孔浮游有孔虫表层水种Sr/Ca 与表层海水温度和盐度记录[19]的线性相关分析SST 为表层海水温度，δ18Osw-ice 为表层海水盐度替代性指标（高值指示高盐），Sr/Castack 为G.ruber s.s.、G. ruber s.l.和T. sacculifer（without sac）壳体Sr/Ca 记录的堆叠结果。Fig.5 Linear correlation of shell Sr/Ca ratio of planktonic surface-water-dwelling foraminifer species and sea surface temperature and salinity[19] from Core MD06-3047BSST: sea surface temperature, δ18Osw-ice: the proxy of sea surface salinity (high δ18Osw-ice means higher salinity), Sr/Castack: the stack of Sr/Ca records of G.ruber s.s., G. ruber s.l., and T. sacculifer (without sac).

由于T. sacculifer（with sac）壳体Sr/Ca 与G.rubers.s.、G. rubers.l.和T. sacculifer（without sac）的变化趋势存在明显差异（图4），故将其单独进行分析。如图5c 和图5d 所示，T. sacculifer（with sac）壳体Sr/Ca 与表层海水温度记录无显著相关性，而与表层海水盐度呈反相关。因此，T. sacculifer（with sac）壳体的Sr/Ca 可能主要受表层海水盐度的影响，这一关系与Dissard 等的研究结果一致[27]。

4 结论

通过对西太平洋暖池北部边缘海区MD06-3047B 孔中浮游有孔虫表层水种G.ruber（G.rubers.s.与G. rubers.l.）和T. sacculifer（T. sacculifer（with sac）与T. sacculifer（without sac））壳体的Sr/ Ca 进行分析，发现MIS 3 期以来，G. ruber不同形态类型壳体的Sr/ Ca 差异较小；而T. sacculifer不同形态类型壳体的Sr/ Ca 相差较大。因此，在利用G. ruber和T. sacculifer壳体的Sr/Ca 结果重建古海洋信息的过程中，如在样品量有限的条件下，可以选择G.ruber壳体不同形态类型进行测试，但应尽量选择T.sacculifer单一形态类型壳体。不同形态类型壳体Sr/Ca 与海水温度和盐度古海洋学记录对比显示，研究区G.rubers.s.、G. rubers.l.和T. sacculifer（without sac）壳体Sr/Ca 可能主要受海水温度的影响；T. sacculifer（with sac）壳体Sr/Ca 主要受到盐度的影响。

致谢：感谢中法合作MARCO POLO 2 航次的全体工作人员在取样过程中提供的帮助。