南极湖泊沉积物中有机硫组成及其与铁硫化物的联系*

沈丽丽，孙婷婷，郭晓宇，黄 涛

(安徽大学资源与环境工程学院，湿地生态保护与修复安徽省重点实验室，合肥 230601)

水域生态系统中硫的转化与有机质成岩埋藏、营养元素循环、水体酸化、重金属生物有效性等紧密耦合，产生重要的环境效应[1-4]. 有机硫(organic sulfur，Sorg)包括由硫酸盐同化还原形成的生物有机硫(Sbio)和由有机质硫化形成的成岩有机硫(Sdiag)，是湖泊沉积物中硫的主要赋存形态之一，占比可高达90%[5-6]. 生物有机硫和成岩有机硫中的还原态硫功能团均可络合Hg+、Cu2+、Pb2+、Ag+及As3+等金属离子，对有毒有害元素的形态和生物有效性具有重要的抑制作用[7]. 成岩有机硫的形成可提高活性有机质的稳定性及其在沉积物中的埋藏保存，而且硫化过程记录了氧化还原条件和元素生物地球化学循环等重要的环境信息[8]. Sinninghe Damsté等[9]较早提出有机质成岩埋藏的自然硫化机制，即S2-优先与铁结合形成铁硫化物，当S2-含量超过铁硫化物计量配比时，便与有机质结合硫化形成有机硫，这主要发生在有机质早期成岩阶段. 湖泊沉积物中的活性有机质在好氧和厌氧条件下可降解形成温室气体二氧化碳和甲烷[10]，因此，开展硫转化及其与有机质结合、调控的相关研究对认识区域碳、硫循环及潜在的气候环境效应具有重要意义.

活性铁是指能被1 mol/L盐酸提取的铁，主要包括无定形或弱晶型Fe(Ⅲ)氧化物和Fe(Ⅱ)(黄铁矿除外)以及Fe3S4(含量很少，可忽略不计)[11]. 活性铁与硫化物可发生铁的硫化作用，生成酸可挥发性硫(acid volatile sulfur, AVS)和黄铁矿硫(pyrite sulfur, CRS). AVS和CRS是沉积物中还原性无机硫(RIS)的两个主要终产物[12-14]，它们的形成会抑制有机质硫化，而高含量活性铁对CRS的形成和有机质硫化则具有促进作用，使得有机质硫化受制于C-S-Fe耦合体系[15]. 黄铁矿化度(DOP)和硫化度(DOS)是评估沉积物中黄铁矿和铁硫化物的形成是受有机质还是活性铁控制的重要指标[16-17]. DOP值是表征沉积物氧化还原条件并识别黄铁矿形成的限制因素[18-19]，其判断标准为DOP>0.75为缺氧环境，0.45

硫同位素组成被广泛应用于水域生态系统中硫的来源及其转化循环过程研究[8,21-23]，在解析有机硫来源组成及其转化方面具有重要作用[15,24]. 沉积物中有机硫相对于同时存在的硫化物富集34S约为10‰，但在某些条件下可达30‰以上[25-27]，这主要与有机硫中生物有机硫和成岩有机硫的形成途径有关. 硫酸盐同化还原过程导致的硫同位素分馏效应仅为1‰～3‰[28]，形成的生物有机硫与其前体硫酸盐具有相似的硫同位素组成；而硫酸盐异化还原生成的硫化氢及其氧化/歧化产物中间态硫的34S亏损显著[8,23,29-30]，成岩有机硫由硫化氢/中间态硫与有机质化合形成，其硫同位素组成继承前体硫化氢/中间态硫偏负的信号. 总有机硫硫同位素组成取决于其生物有机硫和成岩有机硫的组成比例，介于硫化物和硫酸盐之间[31]. 因此，硫同位素是研究有机硫来源组成、转化以及成岩有机硫与铁硫化物联系的重要手段.

南极菲尔德斯半岛和阿德雷岛位于全球变暖最快的南极半岛西北部，随着冰川后退、陆地出露和雪冰消融，沿海岛屿形成了众多湖泊和集水区；其中，阿德雷岛一些湖泊接受了企鹅活动带来的大量有机质和营养元素[32]，成为温室气体重要的排放源. 我们前期对阿德雷岛受企鹅活动影响的Y2湖和未受影响的菲尔德斯半岛燕鸥(YO)湖沉积物中硫和重金属形态组成的分析研究[33-34]表明，有机硫是Y2湖沉积物中硫的主要赋存形态，企鹅活动带来的高含量有机质驱动了强烈的硫酸盐还原，进而限制了微量重金属的生物有效性. 而有关Y2、YO湖沉积物有机硫来源的具体组成及其与有机质、铁硫化物的联系等还不清楚. 因此，本研究在前期工作的基础上，利用总有机硫硫同位素组成、活性铁的DOP和DOS解析并揭示南极湖泊有机硫的来源组成及其与铁硫化物间的联系，对认识湖泊硫生物地球化学转化及其与有机质成岩埋藏、重金属生物有效性的关联提供了视角，并为进一步开展南极气候变化敏感区湖泊硫转化的气候环境效应研究奠定了基础.

1 材料与方法

1.1 研究区域与样品采集

南极菲尔德斯半岛和阿德雷岛位于乔治王岛西南部，面积约33 km2，是乔治王岛最大的无冰区，也是中国南极科学考察站长城站所在地. 当地气候寒冷，年平均气温为-2.2℃. 菲尔德斯半岛和阿德雷岛上的大部分湖泊是由冰川活动形成的，其中阿德雷岛Y2湖位于该岛西部，为一小型月牙状淡水湖，面积约3500 m2，水深在0.5 m左右，海拔高度约12 m，南侧邻海，地势东高西低. Y2湖周边历史时期有企鹅活动，湖泊沉积物中含有大量企鹅转移而来的有机质、营养盐和重金属污染物[32]；而菲尔德斯半岛的YO湖位于中国南极长城站东南方向，为一现代小型湖泊，面积约9000 m2，距离长城湾最近的海岸线约200 m；YO湖未受企鹅活动影响，营养元素含量较低.

本研究分析的Y2、YO湖沉积物为Chen等[33]研究的平行样. 样品于中国第29次南极科学考察期间(2012-2013年)分别采集于南极阿德雷岛的Y2湖和菲尔德斯半岛的YO湖(图1). 企鹅粪土沉积剖面Y2按1 cm间隔分成60个样，原生湖泊沉积剖面YO按1 cm间隔分成30个样. 在进行化学分析之前，一部分样品-20℃冷冻保存，而另一部分样品则在去除岩石碎片和生物残体后冷冻干燥并研磨过筛(孔径为0.125 mm). 本文选择Chen等[33]硫形态研究中的16个Y2和11个YO沉积物序列进行分析.

图1 南极半岛乔治王岛阿德雷岛和菲尔德斯半岛的研究区和采样点(参照Chen等[34])

1.2 样品分析与处理

1.2.1 活性铁分析 在氮气手套箱内称取解冻沉积物0.5 g，加入25 mL 1 mol/L HCl，在室温下振荡24 h，然后离心(6000 r/min, 15 min). 吸取25 mL上清液置于100 mL锥形瓶中，加入1 mL盐酸(1∶3)和1 mL质量分数为10%的盐酸羟胺，煮沸至15 mL，冷却至室温后转移至100 mL容量瓶并定容，混匀后移取1～5 mL于50 mL比色管中定容，使用邻菲啰啉分光光度法测定提取液中总活性铁Fe(HCl)的含量[19]，实验室内样品重复测定误差小于5%.

1.2.2 DOP和DOS的计算 黄铁矿化度(DOP)和硫化度(DOS)按下列公式[35]计算：

DOP=Fe(CRS)/[Fe(CRS)+Fe(HCl)]

(1)

DOS=[Fe(AVS)+Fe(CRS)]/[Fe(CRS)+Fe(HCl)]

(2)

式中，Fe (CRS)为黄铁矿结合态铁的含量，按FeS2的1∶2化学计量比CRS/2计算；Fe (HCl)为总活性铁含量；Fe(AVS)为假定AVS全部以FeS形式存在的硫化物结合态Fe的含量；AVS和CRS含量数据引自Chen等[33].

1.2.4 有机硫源占比计算 沉积物中Sorg由Sbio和Sdiag按一定比例组成，根据同位素质量平衡按公式(3)计算两种硫源的组成比例.

f·δ34Sdiag+ (1-f)·δ34Sbio=δ34Sorg

(3)

式中，f代表成岩有机硫占总有机硫的比例(Sdiag/Sorg)，1-f代表生物有机硫占总有机硫的比例(Sbio/Sorg)；δ34Sdiag和δ34Sbio分别是成岩有机硫和生物有机硫的硫同位素组成，δ34Sorg是沉积物总有机硫硫同位素组成. Y2沉积物中生物有机硫来源于海洋源的企鹅粪，与海水硫酸盐有相似的硫同位素组成，YO硫酸盐主要来自海盐飞沫[33]，其生物有机硫也具有海水硫酸盐硫同位素组成信号，即δ34Sbio≈20‰[37]. 而有机质硫化形成的成岩有机硫硫同位素继承了硫酸盐异化还原产物硫化氢的硫同位素组成信号，即δ34Sdiag≈δ34SAVS≈δ34SCRS. 其中，沉积物δ34SAVS和δ34SCRS数据引自Shen等[36].

2 结果

Y2和YO沉积物δ34Sorg、活性铁含量、DOP和DOS的垂直分布如图2所示. Y2剖面沉积物有机硫硫同位素组成呈3段式垂直分布特征，在1～15 cm较低(1.5‰～3.7‰，平均值为2.7‰)，在19～45 cm中等且波动较小(4.1‰～5.7‰，平均值为4.7‰)，在48～56 cm显著富集34S(8.5‰～11.4‰，平均值为9.6‰). YO剖面沉积物δ34Sorg除28 cm层位明显偏负外，从底部到表层呈现波动变负的趋势.

图2 Y2和YO沉积物中δ34Sorg、活性铁含量、DOP、DOS和Sdiag占比的剖面变化

Y2沉积物中总活性铁含量范围为15.34～42.41 mg/g，平均值为25.4 mg/g，其中在8～19 cm含量较高，其他层位含量平均值为21.67 mg/g且波动较小. YO沉积物活性铁含量为4.26～5.69 mg/g，平均值为4.86 mg/g，远小于Y2. Y2沉积物DOP (0.45±0.21)远低于DOS (5.10±4.13)，其DOS在8 cm层位达到较高值(11.82)，在底部48～56 cm 沉积物的DOP和DOS均变大且在48 cm处达到最高值(12.98). YO沉积物的DOP (0.49±0.16)和DOS (0.83±0.34)相差不大，随深度呈波动增大趋势. Y2剖面1～8 cm沉积物总有机硫中Sdiag占比均大于50%，平均值为58.9%；12～45 cm沉积物Sdiag的占比变化不大(42.7%～56.2%)，平均值为49.9%，底部48～56 cm沉积物Sdiag占比均小于50%，均值为41.0%. YO剖面除28 cm 沉积物Sdiag占比高达95.4%外，其他层位Sdiag占比范围为39.4%～59.2%，平均值为46.2%.

3 讨论

3.1 有机硫来源组成

YO剖面0～8 cm沉积物δ34Sorg值升高，与先前Chen等[33]硫形态指示的明显有机硫矿化相一致，研究表明表层沉积物的微量溶解氧能促进微生物分泌一些物质并有利于有机硫的矿化[41]；8～26 cm和30 cm层位沉积物δ34Sorg值变化不大，指示这些层位生物有机硫和成岩有机硫的组成比例稳定. YO剖面28 cm处沉积物δ34Sorg值显著降低，表明成岩有机硫组分显著增加；前期沉积物硫形态研究指出该层位的硫酸盐还原最强[33]，产生的硫化氢易与有机质结合形成大量的成岩有机硫.

3.2 生物有机硫与成岩有机硫的占比

有机质硫化成岩需要有足量的活性有机质，有机质不仅参与硫酸盐异化还原产生还原态硫，还直接与还原态硫结合形成有机硫[42-43]. Y2和YO沉积物有机硫中成岩有机硫所占的比例如图2所示. Y2沉积物0～8 cm成岩有机硫所占的比例大于50%，很可能是由于高活性有机质促进了成岩有机质硫化. 48～56 cm沉积物较低的成岩有机硫占比指示了大量的企鹅粪输入和沉积扰动导致生物有机硫占比的增加. YO沉积物0～26 cm和30 cm成岩有机硫的比例低于50%，前期研究表明这些层位硫酸盐还原较弱[33]，较低的硫化氢供给导致有机质硫化程度较低. YO沉积物28 cm成岩有机硫的比例高达95.4%，这是因为28 cm处硫酸盐还原强度大[33]，充足的硫化氢与沉积物中的有机物结合产生大量的成岩有机硫，其占比大幅增加.

3.3 有机质硫化与铁硫化的联系

Y2剖面沉积物低DOP值和高DOS值(图2)指示Y2容易发生铁的硫化过程；而YO剖面沉积物DOP>0.45，指示在还原环境中，大部分活性三价铁氧化物被还原为活性二价铁并有效转化为稳定的CRS. 目前，一些研究表明铁硫化物与成岩有机硫能同时形成[44-46]. 铁氧化物快速输入能提高中间态硫的产生，而这有可能同时促进有机质硫化和铁硫化物的形成[46-48]. 本研究中Y2剖面沉积物8～15 cm的低δ34Sorg值和高DOS指示该层位有机质硫化和铁的硫化同时发生，表明该层位高含量的活性铁有利于铁硫化物的形成且并促进了有机质硫化过程；Y2剖面48～56 cm沉积物高DOS值是企鹅活动带来的高有机质和强硫酸盐还原促进了铁的硫化作用，铁硫化物(包括黄铁矿)的大量形成对有机质硫化具有竞争性抑制作用[25,42]，从而降低了成岩有机硫的比例. YO剖面沉积物DOP和DOS随深度呈波动增加的趋势，即深部层位沉积物还原性无机硫更多以CRS为主；同时，28 cm处强硫酸盐还原提供了充足的硫化氢/中间态硫与有机质结合为成岩有机硫并降低δ34Sorg，成岩有机硫占比高达95.4%，表明该层位也存在有机质硫化与CRS的形成同时发生.

硫在沉积物有机质的成岩和埋藏中发挥着重要作用[49]，不同形态硫的转化会影响环境变化乃至生态安全[50-51]. 湖泊硫的转化循环涉及一系列氧化还原反应，并影响内源磷和铁的迁移释放、有机质的成岩埋藏、重金属的毒性和生物可利用性，产生重要的生态环境效应[17,34,52]. 环南极海岛、无冰区分布有众多接收海鸟、海兽传输的高有机质、营养盐和重金属污染的湖泊和集水区[53]，这类气候敏感区水生系统水文和氧化还原条件季节性变化显著，营养盐生物地球化学转化活跃. 对南极洲沿海地区海兽排泄物的研究中发现，海兽排泄物中的有机质具有很高的生物有效性，在冻融过程中，有机质会因微生物活性升高发生矿化过程，使二氧化碳或甲烷的排放显著增加[54]. 有机质硫化过程能够提高有机质的稳定性，有效抑制有机质的矿化分解，对碳转化循环产生实质影响. 本研究结果表明，有机硫组分的来源转化与有机质、活性铁紧密耦合，对区域湖泊有机质成岩埋藏、铁硫化物形成及其对金属生物有效性限制有直接指示作用，对进一步开展营养元素生物地球化学转化的气候环境效应评估具有重要的意义.

4 结论

1) Y2剖面底部生物有机硫占比的提高源于大量企鹅粪的输入，使沉积物中生物有机硫的增加超过了成岩有机硫；而YO剖面28 cm处的成岩有机硫占比的绝对优势源于强硫酸盐还原提供了充足的硫化氢/中间态硫作为有机质硫化的反应前体.

2) Y2和YO剖面部分层位沉积物有机质硫化与铁的硫化同时发生，主要与高活性铁含量和强硫酸盐还原过程相关.

致谢：本文研究样品来源于中国极地沉积物库(合肥).