王彤,胡艺豪,贾奇,郭景腾,唐正,熊志方,2,李铁刚,2
1. 自然资源部第一海洋研究所,海洋地质与成矿作用重点实验室,青岛 266061
2. 崂山实验室海洋地质过程与环境功能实验室,青岛 266237
工业革命以来,化石燃料的燃烧导致大气二氧化碳分压(pCO2)的急剧升高,引起了人类社会对全球气候变暖的担忧。海洋是地球表生系统最为重要的碳储库之一,也是最大的活跃碳库,从而调节着大气pCO2的变化[1]。南大洋吸收了工业革命以来人为排放CO2的40%,是全球重要的碳汇海区[2]。古气候学研究进一步发现,轨道和千年尺度上大气pCO2的演化与南极温度、南大洋深部流通状况和表层生产力记录密切相关。这些结果表明,南大洋在调节海洋碳循环和全球气候演化中扮演着至关重要的角色[3]。
海洋浮游植物通过光合作用吸收CO2,将其转化为颗粒有机碳,其向深层水输出的过程中,一部分经有机质矿化生成“呼吸CO2(Respired CO2)”,封存在海洋深部,另一部分直接沉降到海底埋藏,从而降低了大气pCO2。南大洋深部水体富含大量营养物质与CO2,其上涌到表层时,会刺激浮游植物生长,并向大气中释放CO2。当通风作用直接排放的CO2的量大于生物泵吸收的CO2的量时,南大洋就会成为碳源,增加大气pCO2,反之,它会成为碳汇,使大气pCO2降低[4-7],这表明提高生产力可能会增加南大洋封存CO2的潜力。Jaccard 等发现南大洋存在两种不同的生产力模式(即冰期南大洋南极带生产力降低而亚南极带生产力增加,间冰期则相反),从而提出一个假说解释冰期-间冰期旋回中大气pCO2演化[7]。该假说认为,冰期广泛增加的风尘通量带来的铁施肥使得亚南极海域生产力增长,从而促进了冰期大气pCO2的降低;南极带内生产力的降低则可能反映了风驱上升流与冬季垂直混合的减弱,从而也导致了冰期大气pCO2的降低。间冰期,则正好相反。该假说成立的关键是确证南大洋南极带与亚南极带“镜像”的生产力演化模式,并揭示其控制机理。
威德尔海作为全球大洋深层水的主要源区之一,在链接大气、海洋表层与海洋深部之间的碳循环中扮演着重要的角色。表层海水对CO2的吸收主要发生在春夏季浮游植物勃发时期[8],深秋和冬季由于混合层加深与高CO2含量的绕极深层水上涌,导致表层海洋向大气释放CO2;随后,冬季海冰的覆盖导致表层海洋无法进一步向大气中释放CO2[9],这使得威德尔海存在碳收支的季节性不平衡。工业革命以前,威德尔海因季节性的碳收支不平衡,成为大气CO2的源[10]。工业革命以后,大气CO2浓度的迅速升高导致这种不平衡方向被改变,该海区成为碳汇[11];且这些被吸收的碳一部分随着南极底层水(AABW)被输送至全球深海封存,另一部分跟随着中层水向北传导至低纬海洋,影响着低纬海洋的气候环境[12]。总之,长时间尺度上,威德尔海能有效缓冲大气pCO2的上升。因此,重建地史时期威德尔海生产力的演化、探讨环境要素变化对生产力的影响,有助于更好地理解南大洋在过去大气pCO2演化中的作用,并预测南大洋在全球气候变暖背景下作出的响应。
研究区位于南大洋最大边缘海威德尔海西北部(图1a),该海区沉积物以等深流沉积与半远洋沉积为主[13-14]。威德尔海是AABW 的重要发源地之一[15],影响全球深海大部分区域的水文状况,从而调控全球气候[16]。威德尔环流是由风驱动的顺时针方向流动的沿岸流,其在自东向西流动的过程中因南极半岛的阻碍而转向北上,最终汇入德雷克海峡东部的南极绕极流(ACC)(图1b)。绕极深层水(CDW)流经南大西洋时在极锋附近发生上涌,一部分在艾克曼层北移成为南极中层水(AAIW)和亚南极模态水(SAMW)流向低纬海区[17],另一部分在威德尔海东部边界区侵入威德尔环流,并在流向威德尔海西南部海域时抬升,作为改良的暖深层水;其在随威德尔环流运移过程中侵入陆架,与陆架水混合,相互作用形成AABW,并沿陆坡流动,随深部环流向北输出(图1b)[18]。
图1 南极威德尔海站位分布、洋流与水团与硅酸盐含量a. 生产力模式对比中所涉及站位, b. 岩芯D5-12 位置及洋流分布,c. 硅酸盐浓度垂向剖面(39°W 附近)。图a 和b 中红色五角星标示岩芯采样位置,图b 中红色线示意图c 的经向截面位置。图1a 采用Mercator 投影方式。Fig.1 Map of the Antarctic Weddell Sea showing core location, ocean circulation and water masses, and dissolved-silicate concentrations a: the sites involved in the comparison of productivity models; b: location of Site D5-12 and ocean circulation; c: silicate concentration section along 39° W.The red stars indicate the sampling location, red dotted line indicates the location of the section of silicate concentrations shown in Fig.1c. Figure 1a is based on the Mercator projection.
磷酸盐和硝酸盐浓度在威德尔海的垂向分布十分相似,这主要归因于它们受控于相似的环流与生物地球化学过程[19]。磷酸盐和硝酸盐在表层水体中由于初级生产者的利用,其浓度呈现极小值;在表层水体下方由于有机质再矿化和富营养物的CDW 侵入,其浓度表现为极大值。与磷酸盐和硝酸盐相同,硅酸盐浓度也在表层水体中出现了极小值;不同的是,由于硅质浮游植物生物硅的溶解比有机质的矿化更难,硅酸盐浓度极大值出现的深度更深(图1c)。溶解氧含量的水深分布模式与上述营养物质相反,表层海洋中氧气的海-气交换以及植物的光合作用使得表层水体中溶解氧含量最大,随着深部有机物的再矿化,溶解氧含量出现极小值[19]。
沉积物岩芯ANT34/D5-12(以下简称D5-12)是2018 年中国第34 次南极科考航次“向阳红01”科考船采集的重力柱状岩芯。D5-12 孔位于威德尔海西北部,水深2 556 m,岩芯全长456 cm。该孔沉积物层理清晰,未见明显沉积间断,自上而下颜色变深,局部出现黑色纹层,整体可分为两段。岩芯上部0~265 cm 沉积物主要为含有孔虫粉砂质黏土,所含有孔虫优势种为浮游有孔虫Neogloboquadrina pachyderma(sinistral)。根据沉积物岩性差异可将其分为6 层,0~55 cm 为灰色含有孔虫粉砂质黏土,颜色分布均匀;55~86 cm 为浅灰色含有孔虫粉砂质黏土,颜色分布不均匀;86~120 cm 为灰黑色含有孔虫黏土质粉砂,自上而下颜色逐渐加深;120~258 cm 为灰黄色含粉砂质黏土,该层有孔虫含量相对较少,且随着深度增加逐渐减少,并出现丰富的冰筏碎屑;258~265 cm 为灰黑色夹黑色条带的粉砂质黏土,不见有孔虫但有大量冰筏碎屑出现。岩芯下部265~456 cm 沉积物主要为灰黄色黏土,颜色均匀,含较多放射虫及少量海绵骨针(图2a)。
图2 南极威德尔海D5-12 孔年龄模式图a. D5-12 孔岩芯岩性柱,b. D5-12 孔MS(磁化率)-深度变化曲线,c. EDML 冰芯nssCa2+(非海盐钙)通量曲线[28], d. D5-12 孔MS 与EDML 冰芯nssCa2+通量曲线对比,e. D5-12 孔线性沉积速率。图b 和c 中的红点代表年龄控制点,图e 中的水平虚线示意平均线性沉积速率。Fig.2 The age model of Core D5-12 in the Antarctic Weddell Sea a: Lithologic column of Core D5-12, b: MS-depth curve of Core D5-12, c: EDML nssCa2+ flux[28], d: comparison of nssCa2+ flux of EDML ice cores and MS of Core D5-12, e: linear sedimentation rate of Core D5-12. The red dots in Figures b and c represent age control points, the horizontal dashed line in Figure e indicates the mean linear sedimentation rate.
D5-12 孔在自然资源部第一海洋研究所对半剖分后,半管随即按 1 cm 步长进行物理参数的扫描,获得磁化率等数据。另一半管以1 cm 间隔分样,共获得456 个样品。进一步,以2 cm 间距取样进行生物硅(BSi)含量测试,共取样品228 个。BSi 的分析采用湿碱消解法[20],在中国海洋大学海底科学与探测技术教育部重点实验室完成。BSi 具体分析测试方法为:冻干样品研磨至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 次;该过程快速完成,以减少可溶硅在固体表面的不可逆损失。通过钼酸盐蓝光分光光度法测量分步提取物中的溶解硅,并根据硅的质量百分比随时间变化得出线性回归方程,其截距为沉积物中BSi 的含量。沉积物中生物硅(即生源蛋白石)含量按公式(Opal%= 2.4×BSi%)计算[21]。实验室该分析的长期精度(RSD)优于±3%。
沉积物准确定年是重建古环境的基础。目前,在南大洋高纬海域,对于年代小于40~50 ka 的地层,广泛采用生物壳体或有机质的AMS14C 实现绝对定年,以获取沉积物的年代框架。对于年代大于40~50 ka 的地层进行定年时,常使用有孔虫氧同位素与全球大洋底栖有孔虫氧同位素拟合记录LR04 曲线对比[22],并辅以AMS14C 测年来联合确定深海沉积物的年代框架,这种情况需要确保有孔虫壳体保存情况良好并且有孔虫沉积记录连续。在无法满足上述前提的情况下,也可使用其他特征记录指标(如磁化率和Ba/Ti 等)与典型气候曲线对比,并结合AMS14C 数据实现深海沉积物的年代框架确定[23]。岩芯D5-12 取样位置虽位于碳酸盐补偿深度以上,有孔虫壳体保存完整,但我们未发表的其他数据表明, 该孔中浮游有孔虫壳体N.pachyderma(sinistral)的δ18O 信号受局部海水温度与冰雪融水信号影响强烈,因而“擦拭”了δ18O 记录中的全球冰量演化信息,导致其无法与LR04 曲线进行对比定年。沉积物中有机质的AMS14C 绝对定年方法使用时,要求样品有机质含量高且来源单一,低含量的有机碳可能反映矿化程度严重,产生虚假年龄。D5-12 孔未发表的有机质AMS14C 测年结果显示,年代多次出现倒转现象,因而无法应用于该孔的年代模式确定。因此,本研究采用该孔沉积物中磁化率与南极冰芯风尘通量参数记录对比的方式,获得该站位年代模式。
南极冰芯的风尘通量记录保存了南半球过去的大气环流信息[24]。因此,南极冰芯的风尘颗粒浓度以及nssCa2+通量常被用于示踪大气环流[25]。冰期时,赤道和两极之间经向温度梯度更高,导致大气环流加强;间冰期则相反[4]。因此,风尘通量参数可以用于进行远距离年代地层对比,且在冰期-间冰期的时间尺度上不存在超前或者滞后,从而作为一种用于年代地层对比的类似于LR04 的有效标准曲线[26]。前人研究表明,南大洋深海沉积物中的磁化率主要来自于风尘搬运的磁性物质,风尘输入加强对应于磁化率升高,减弱对应于磁化率降低[27]。冰芯中nssCa2+通量是大气风尘输送强度的可靠指标[28-29],因此,南大洋沉积物中磁化率与冰芯中nssCa2+对比的方法已广泛应用于岩芯年龄模式的确定[23]。研究区威德尔海与南极EDML 冰芯站位的风尘沉积都主要来源于南美的巴塔哥尼亚地区[30-31];因此,本文将D5-12 孔沉积物磁化率(图2b)和南极EDML冰芯nssCa2+通量(图2c)进行对比,选取曲线中的峰值或曲线形态的拐点作为年龄控制点,建立该孔的年代框架;在此基础上,采用线性内插的方法计算线性沉积速率。结果表明,D5-12 孔MS 记录与EDML冰芯nssCa2+通量曲线对比良好,该孔底界年龄约为110 ka(图2d),平均沉积速率为5.7 cm/ka(图2e)。D5-12 孔年代模式表明该孔顶部缺失14 ka 以来的沉积物,我们推测可能的原因有二。其一, 这段时期(主要为全新世)研究区沉积速率较低,导致沉积层较薄,其在重力取样过程中由于较强的洋流作用流失了。其二,该孔附近的ODP113 航次的695、696 和967 站研究表明研究区在全新世存在沉积间断,普遍缺失沉积物[32]。
长期趋势上,D5-12 孔生源蛋白石含量在110~71 ka 较高,且表现为降低趋势,其值从9.11%减少到4.27%,平均值为3.21%;在71~14.3 ka 较低,且变化较为平稳,其值在0.28%~3.08%之间波动,平均值为1.35%。冰期-间冰期尺度上,生源蛋白石含量记录与气候旋回表现出明显的一致性,呈现出冰期(冷期)较低,间冰期(暖期)较高的特征,即在MIS 3 和5 期间较高,在MIS 2 和4 期较低(图3a)。
图3 南极威德尔海D5-12 孔生源蛋白石(Opal)、Ti、Opal/Ti 演化曲线a. 蛋白石含量,b. Ti 含量,c. 蛋白石/Ti 比值。Fig.3 The biogenic opal, titanium, and Opal/Ti records of Core D5-12 a: Biogenic opal contents, b: titanium contents, c: opal/Ti ratios.
利用岩芯中生源蛋白石含量指示生产力的前提是剔除其他主要组分对生物硅的稀释作用,具体措施是使用岩屑元素含量对生源蛋白石含量进行标准化[33]。Al、Ti 等元素在陆地风化过程中地球化学行为保守,且在海水中含量极低,因此海洋沉积物中Al 和Ti 通常被认为主要来源于陆源碎屑物质,可用作扣除稀释作用的标准化措施[31,34]。近年来研究发现,生源蛋白石可以通过“捕获”的方式将Al 富集在沉积物中,因此这种“过剩铝”反映的是生产力而非岩屑通量的变化[35-36]。鉴于此,选用陆源来源更稳定且受生物地球化学行为扰动较小的Ti 来对生源蛋白石含量进行标准化(图3b 和c)。如图3c 所示,经过Ti 标准化的蛋白石含量与初始的蛋白石含量变化趋势完全相同,因此排除了稀释效应的影响。另外,深部氧化还原环境可能会对生产力替代指标的保存产生影响,但生源蛋白石的保存与沉积物的氧化还原状态无关,其主要受控于深部水体硅酸的饱和程度[37]。尽管目前没有研究区MIS 5 期以来深部硅酸饱和程度(或浓度)的古海洋记录,但南大洋高的生源蛋白石埋藏效率对应于表层海水中较高的生源蛋白石雨率(rain rate),表明了沉积物中保存的生源蛋白石含量反映了其生产而非溶解信息[38]。综上所述,认为该孔生源蛋白石含量可靠地记录了研究区的生产力信息。
最近的研究显示,南大洋南极带通过成层化(表现为低生产力)封存深部“呼吸CO2”,可以解释0.8 Ma 以来冰期大气pCO2降幅的前半部分(约40~50 mg/L);而亚南极带通过风尘铁刺激的高生产力吸收大气CO2,则可解释后半部分[7,39]。这个大气pCO2演化驱动机制假说表明,不依赖于其他大洋,南大洋自身就能引起完全的大气pCO2冰期旋回。测试该假说合理性的前提就是要验证南大洋冰期南极带生产力低而亚南极带生产力高的演化模式。位于南大洋南极带内的D5-12 孔110 ka 以来生产力呈现出冰期低、间冰期高的演化模式(图4a),因此支持该假说。南极带这种生产力冰期旋回模式同样出现在其他岩芯或海域。例如,斯科舍海U1537 站的蛋白石含量记录[40]和威德尔海PS1768-8 孔的过剩Ba 通量记录[41]表明,暖期(MIS 5 和3 期)时生产力处于较高水平,但在冷期(MIS 4 和2 期),生产力处于较低水平(图4b 和c)。罗斯海RS15-GC40 和RS15-GC41 孔的蛋白石含量[42]也表现出这种冷期低、暖期高的生产力演化模式(图4d和e)。这些结果初步确证,南大洋南极带生产力呈现冰期低间冰期高的旋回模式。南极带这种生产力的演化模式很可能由南大洋西风带、海冰演化引起的营养物可利用性(Availability)变化所致(详述见第4.2 节)。
图4 南大洋南极带与亚南极带生产力冰期-间冰期旋回模式对比a. D5-12 孔生源蛋白石含量,b. U1537 站蛋白石含量[40],c. PS1768-8 孔过剩Ba 通量[41],d. RS15-GC41 孔蛋白石含量[42],e. RS15-GC40 孔蛋白石含量[42],f. PS2082-1 孔过剩Ba 通量[41],g. MD02-2588 孔蛋白石含量[45],h. PS97/093-2 孔蛋白石含量[44],i. DCR-1PC 孔蛋白石通量[46]。Fig.4 Comparison in productivity between the Antarctic zones and subantarctic zones in the Southern Ocean in glacial-interglacial cycles a: Opal contents of Core D5-12, b: opal contents of Core U1537[40], c: Baexcess flux of Core PS1768-8[41], d: opal contents of Core RS15-GC41[42], e: opal contents of Core RS15-GC40[42], f: Baexcess flux of Core PS2082-1[41], g: opal contents of Core MD02-2588[45], h: opal contents of Core PS97/093-2[44], i: opal flux of Core DCR-1PC[46] .
前人研究指出,南大洋不同海区对气候变化的反应可能是不同的[43],存在明显的生物地球化学分隔(Biogeochemical divide),使得南极带与亚南极带对气候产生不同的响应方式[17],这也体现在生产力演化模式上。东南太平洋德雷克海峡西入口处PS97/093-2 孔[44]和东南大西洋MD02-2588 孔[45]的蛋白石含量记录表明,亚南极带冷期(MIS 4 和2 期)生产力明显高于暖期(MIS 5 和1 期)(图4g 和h)。南大洋宇航员海DCR-1PC 孔的蛋白石通量记录[46]和威德尔海PS2082-1 孔的过剩Ba 通量记录[38]也表现出亚南极带内冷期高、暖期低的生产力演化特征,这种现象在PS2082-1 孔尤为明显(图4f 和i)。亚南极带这种生产力冰期高间冰期低的旋回模式主要受控于西风带强弱引起的风尘铁施肥效应变化所致[7,45,47]。综上所述,现存生产力的古海洋学记录表明,南大洋南极带和亚南极带生产力演化模式在冰期-间冰期尺度上呈“镜像”关系。
总体来看,海洋生物生产力最直接的内部影响因素就是硅酸盐、硝酸盐、磷酸盐等营养物以及铁等营养物限制因子的可利用性。对于南大洋而言,西风带控制的深部流通状况变化、海冰导致的光照时间和水体层化强度变化、风尘输入引起的铁施肥效应变化都能影响营养物的来源与利用[48-52],从而是生产力演化的外部控制因子。因此,本节通过MIS 5 期以来风尘、西风带和海冰演化记录与威德尔海生产力记录的对比,试图查明南大洋南极带生产力演化的控制机理。
4.2.1 风尘(铁)的影响
Gran 最先提出,铁是限制海洋浮游植物初级生产力的重要因子之一[53]。Martin 等发现冰芯中CO2浓度演化和风尘铁记录呈负相关,并由此提出了铁假说。一些研究认为冰期风尘通量的增加带来了额外的铁,从而解除了铁限制,提高了生物泵效率[49,54]。因此,假若南大洋的生产力演化模式受控于风尘铁,其应与风尘铁的供应模式相同,表现为在风尘铁输入较高的冰期,生产力较高;在风尘铁输入较低的间冰期,生产力较低。
威德尔海与东南极冰盖的风尘沉积主要来源于南美洲巴塔哥尼亚地区[30-31],且来自该源区的风尘在输运的过程中会途经斯科舍海域。因此,我们选取位于威德尔海另一站位ODP 1090[55]和斯科舍海U1537 站[40]的铁质量堆积速率(MAR)记录来代表影响威德尔海的风尘输入强度(图5b 和c)。D5-12 孔蛋白石含量记录总体上与ODP 1090 和U1537站的铁质量堆积速率记录呈反相位的关系:即冷期(MIS 4 和2 期)风尘铁输入增强对应生产力低,暖期(MIS 5 和3 期)则相反。这种对应关系与前人推断的南极带冰期风尘通量高但生产力低,间冰期风尘通量低但生产力高的认识相一致[27]。从长期趋势看,MIS 5 期以来D5-12 孔生源蛋白石含量呈减小趋势,而铁MAR 呈增加趋势;因此,两者的长期演化趋势也不满足南大洋生产力和风尘铁的同相位演化规律。这些结果表明,风尘及其铁的输入不是MIS 5 期以来威德尔海生产力演化的主控因素。另外,我们注意到,在考虑定年误差的情况下,生产力在MIS 5 期内86 ka、97 ka 和102 ka 表现的峰值与风尘铁输入通量可能有较好的对应关系。这种在千年尺度上的对应关系,已经超出了本文聚焦讨论轨道尺度上生产力演化及其机制的范围,因而将其另文讨论。
图5 D5-12 孔生源蛋白石含量记录与其他相关海洋-环境记录综合对比a. D5-12 孔生源蛋白石含量,b. ODP 1090 站Fe 质量堆积速率[55],c. U1537 站Fe 质量堆积速率[40],d. 德雷克海峡可分选粉砂和细砂组分平均粒径[57],e. 南极EDC 冰芯氘过剩[63],f. 底栖有孔虫Δδ13C[64],g. 南大洋大西洋扇面海冰强度演化曲线[52],h. 南极冰芯大气pCO2 演化曲线[65]。图f 中底栖有孔虫Δδ13C 为南大洋大西洋扇区ODP 1089 站和ODP 1088 站底栖有孔虫δ13C 之间的差值[64]。图g 中紫色曲线表示均值,紫色阴影表示四分位误差范围。Fig.5 Overall comparison in opal content between Core D5-12 and other relevant ocean-environment records a: Opal contents of Core D5-12, b: Fe mass accumulation rates (MAR) of ODP Site 1090[55], c: Fe MAR of Site U1537[40], d: the mean grain size of sortable silt plus fine sand fractions in the central Drake Passage as the proxy of ACC flow speed[57], e: Antarctic EDC ice core dln anomaly[63], f: benthic Δδ13C records[64].The benthic Δδ13C is the difference between the benthic δ13C from Site ODP 1089 and ODP 1088 in the Atlantic sector of the Southern Ocean[64]; g: stack of sea ice evolution in the Atlantic Southern Ocean [52], the purple curve marks the mean values and purple shading marks the interquartile ranges; h: compilation of the atmospheric pCO2 records in Antarctic ice cores[65]. In Figure f, the benthic Δδ13C is the difference between the benthic δ13C from Site ODP 1089 and ODP 1088 in the Atlantic sector of the Southern Ocean[64];In Figure g, the purple curve marks the mean values and purple shading marks the interquartile ranges.
4.2.2 西风带的影响
Trull 等认为南极带的生产力主要受控于营养物质而不是营养物限制因子铁的浓度,而西风带可以通过两种方式来影响南大洋深部流通状况[48,56],进而影响深部营养物向表层的输送。首先,西风带可以通过控制南极带海冰的幅度来影响南大洋深部流通状况[56]。冰期时,西风带向赤道方向移动,导致南极带海冰扩张,从而减弱了深部流通状况;间冰期时则相反。其次,西风带可以通过强迫ACC 来影响南大洋深部流通状况[57]。冰期时,西方带向北移动,一方面导致ACC 北移远离南极带,从而使南极带深部流通状况减弱,另一方面,北移的西风带也导致ACC 流速减弱,水体流动性变弱,从而导致深部流通性减弱;间冰期则相反[58]。总之,西风带不论是以何种方式影响南大洋深部流通状况,都表现为:冰期西风带北移,南极带深部流通状况减弱,不利于深部营养物上涌,导致生产力降低;间冰期,则导致生产力提高。
冰期旋回中,西风带南北移动、ACC 强度、深部流通状况和威德尔海生产力确实表现为上述的协变规律(图5a 和d-f)。具体表现为:MIS 4 和2 等冷期西风带的向北移动、ACC 强度的减弱和深部流通状况的减弱对应威德尔海生产力的降低;MIS 5 和3 期等暖期则相反。长期趋势上,MIS 5 期以来西风带总体北移,ACC 逐渐减弱,深部流通状况逐渐减弱;同时,威德尔海生产力总体降低。西风带移动和威德尔海生产力的长期演化趋势也符合西风带对生产力的影响机理。因此,MIS 5 期以来西风带的经向移动对威德尔海的生产力演化产生重要影响。
4.2.3 海冰的影响
威德尔海是南极海冰最为发育的边缘海,其东部海域海冰受季节主导,秋冬季迅速生长并受风力强迫随表层水运动,西部海域海冰受前一年海冰影响,是南大洋最厚海冰区之一[59-60]。海冰系统对海洋生产力的影响主要体现在两个方面。其一,海冰通过对调控光照条件影响着浮游植物生长条件[61]:温暖时期,海冰强度减弱、覆盖范围减小、持续时间缩短,使得光照增强,浮游植物的生产季节变长,促进生产力的提高;寒冷时期,则降低生产力。其二,海冰通过调节水体层化影响着浮游植物生长所需的营养物浓度[62]:寒冷时期,增加的海冰使近表层海水密度梯度增大,水体分层增加,阻碍了富含营养物的深层水上涌,导致表层水体营养盐匮乏而生产力降低;温暖时期,则生产力增加。总之,海冰强度增大导致生产力降低;海冰强度减少导致生产力升高。
尽管MIS 5 期以来,南大洋不同海域海冰的量级、持续时间与变化幅度存在差异,但在冰期-间冰期尺度上表现为MIS 5 和3 期海冰强度低,MIS 4 和2 期海冰强度高;在长期趋势上表现为MIS 5 期以来海冰强度总体升高[52]。很显然,MIS 5 期以来海冰强度的演化过程与趋势明显与威德尔海生产力记录一致,满足上述海冰影响生产力的机制。例如,MIS 5 和3 期低的海冰强度导致威德尔海高的生产力;MIS 5 期以来逐渐增加的海冰强度导致威德尔海总体降低的生产力。这些在不同时间尺度上海冰与生产力演化记录的耦合关系表明,海冰强度也对MIS 5 期以来威德尔海生产力有重要影响。
深部流通状况和表层生产力是控制海洋封存与释放CO2的两大因素,其最终链接到大气pCO2的变化[50,62]。深部流通状况减弱时,有利于“呼吸CO2”在海洋深部封存,导致大气pCO2的降低;深部流通状况加强时,海洋深部“呼吸CO2”最终释放到大气中,导致pCO2的升高。生产力提高时,一方面有利于更多的有机碳埋藏到海底,另一方面能有更多的有机碳矿化成更多的“呼吸CO2”封存在海洋深部,最终都导致大气pCO2的降低;生产力降低时,则导致大气pCO2的升高[50,62,66]。对于冰期旋回中的南大 洋,Jaccard 等和Martínez-García 等提出假说,认为冰期南大洋南极带流通状况减弱(即成层化加强)对“呼吸CO2”的封存导致了大气pCO2降幅的一半(约40~50 mg/L);冰期南大洋亚南极生产力加强对大气CO2的吸收能导致了大气pCO2降幅的另一半(约40~50 mg/L)[7,39,65]。在4.1 和4.2 节中,我们从威德尔海生产力的演化模式及影响因素视角,对该假说进行了验证。在本节中,从生产力与大气pCO2演化的关联视角,进一步对该假说进行验证。
依据上述讨论,MIS 5 和3 期等暖期,威德尔海生产力升高对应于大气pCO2升高;MIS 4 和2 期等冷期,威德尔海生产力降低对应大气pCO2降低(图5a 和h)。这种对应关系表明,威德尔海不可能通过生产力来控制大气pCO2的冰期旋回变化。我们注意到,威德尔海生产力演化实际上反映了南极带深部流通状况信息(详见第4.2 节)。冷期(MIS 4 和2 期)时,降低的生产力指示减弱的深部流通状况,对应降低的pCO2;暖期(MIS 5 和3 期)时,提高的生产力指示加强的深部流通状况,对应升高的pCO2。很显然,威德尔海处于南极带的研究区很可能通过深部流通状况的演化对大气pCO2冰期旋回产生贡献。进一步,我们注意到MIS 5 期以来威德尔海总体降低的生产力指示总体减弱的深部流通状况,正好与逐渐降低的大气pCO2对应(图5a 和h),这表明深部流通状况的减弱加强了“呼吸CO2”在深层水体中的封存,从而导致大气pCO2降低。这也从长期趋势上验证了南极带通过深部流通状况影响大气pCO2演化(图6)。综上所述,威德尔海生产力记录及其与大气pCO2等相关古环境参数的对比,从冰期旋回和长期趋势两种时间尺度上验证了 上 述Jaccard 等 和Martínez-García 等 有 关 大 气pCO2演化机制[7,39,65]的合理性。
图6 威德尔海生产力-深部流通状况-大气pCO2 耦合关系模式图图中箭头的粗细指示相应指标的强弱。Fig.6 Schematic presentation of the coupling relationship between productivity, deep convection and pCO2-atm of the Weddell Sea The thickness of arrows in the figure indicates proportionally the strength of the corresponding indicators.
(1)冰期-间冰期尺度上,MIS 5 期以来威德尔海生产力呈明显的冰期旋回,表现为暖期(MIS 5 和3 期)生产力高,冷期(MIS 4 和2 期)生产力低。长期趋势上,MIS 5 期以来威德尔海生产力总体上呈降低趋势。
(2)西风带的经向移动和海冰张缩主要通过影响深部流通状况,进而控制深部营养物进入表层的可利用性,最终驱动MIS 5 期以来威德尔海生产力的冰期旋回和长期趋势。风尘(铁)似乎对威德尔海生产力演化没有明显影响。
(3)MIS 5 期以来大气pCO2的冰期旋回和长期趋势变化与威德尔海生产力的耦合关系很可能受到南极带深部流通状况的调控,这与南大洋南极带驱动大气pCO2演化的机制假说相吻合。