随着大气中CO₂浓度从1750年的280×10^-6逐渐上升到2005年的380×10^-6[1]，使得的环境问题中碳循环和温室气体(CO₂、O₃、N₂O、CH₄)问题更为突出。全球碳循环主要是指碳在岩石圈、水圈、气圈和生物圈之间不同形态之间的相互转换和运移。在全球碳循环的研究中，随着岩溶地区碳“汇”和“源”的研究逐步深入^[2-3]，发现在岩溶作用下产生溶解性无机碳(Dissolved Inorganic Carbon，以下简称DIC)并存储于水体中，具有碳“汇”效应^[4-5]。DIC是整个碳循环的一个重要组成部分。在整个过程中水、大气和基岩相互作用主要方程式可以简化为：

(1)气与水相互作用：

CO₂(g)+H₂O(aq)→ H₂CO₃(aq)

(2)水与岩相互作用：

H₂CO₃(aq)+CaCO₃(s)→2HCO₃^-(aq)+Ca²⁺(aq)

大气CO₂溶于水体后在水体中形成DIC的存在形式分别为：CO₂、H₂CO₃、HCO₃^-、CO₃^2-等主要形态平衡混合物。由于pH值控制着水体中DIC不同形态所占的比例，在pH值为7～9时，水体中的DIC主要是以HCO₃^-的形式存在^[4]。随着碳稳定同位素分析精度的提高，研究水循环中不同形态碳及测定同位素能够揭示碳“源”、“汇”以及碳通量等问题^[6-9]。最初测试水体中DIC的前处理方法主要是BaCl₂沉淀法^[10]，通过将液体样品中DIC转化为固体样品来测试沉淀物(碳酸盐)的碳同位素值，但已经有学者指出BaCl₂沉淀法在快速沉淀过程中存在着同位素分馏^[11]。Matthews等^[12]开发研究了以惰性气体为载气携带测试气体进入质谱仪的高真空室离子源室。随着这种测试技术的发展和完善，美国Thermo Fisher公司研发了水体中DIC的碳同位素分析的前处理装置GasBench。该装置能够直接测定水中DIC的碳同位素，具有准确度高、分析速度较快等特点，已经在国内实验室得到广泛的应用^[13-16]。为了准确测定水体DIC中极易逸出的游离CO₂，我国研究人员针对其特点提出了准确的测试方法^[17]，并且对室内测试DIC含量和碳同位素的方法都进行了详细研究^[18]，但是并没有指出如何解决野外样品前处理过程对水体中所有DIC引起的碳同位素分馏。

水体中DIC是全球碳循环中的重要研究对象，因此能够精确地分析测定水中的DIC碳同位素，对于探究碳在水圈中的运移、循环机制等都具有十分重要的科学意义。本文对现有水体中DIC碳同位素测试的三种前处理方法(BaCl₂沉淀法、医用无菌高密度聚乙烯瓶装样、GasBenchⅡ顶空样品瓶野外直接生成CO₂气体)进行定量对比分析，研究不同前处理方法对水体中DIC的测试结果的影响和引起的碳同位素分馏大小，以寻求建立一套操作简单、准确测定水体中DIC碳同位素的前处理方法。

1.   实验部分

1.1   水样的采集和水样理化性质测定

为了能够对比前处理方法对DIC的影响，本次研究选择不同浓度的DIC样品进行对比，采集桂林地区盘龙洞洞穴滴水和地下河天窗样品。桂林市地处于低纬度地区，湘桂走廊南端，东经109°36'至111°29'，北纬24°15'至26°23'，平均海拔150 m。该地区属于亚热带季风气候，四季分明，雨量充沛，夏季的季风降水主要受到夏季风的影响，年均气温19.5℃，年均降水量1868 mm^[19]。

盘龙洞处于桂林市南部38 km处的报安村(如图 1所示)，所处为典型的亚热带温润地区岩溶峰丛洼地。洞穴围岩为上泥盆统融县组(D_3r)灰岩，灰岩质纯，局部有页岩夹层。洞穴内年均气温19.5℃，洞穴滴水的水温为19.3～21.2℃，洞穴空气的相对湿度为90%～98%^[20]。在野外使用Merck碱度试剂盒滴定地下河天窗(样品编号P_地)和洞穴滴水(样品编号P₆)的HCO₃^-浓度；使用WTW Multi 3420多参数水质分析仪现场测定采集水样的pH值和温度；使用testo 435-2多功能测量仪检测样品，采集周围环境的大气压强值。在实验室测试样品时，使用testo 435-2多功能测量仪现场测量环境温度和大气压强，结果列于表 1。

表  1  野外水样HCO₃^-、pH、温度和实验室温度、压强测量结果

Table  1.  The measurement results of HCO₃^-,pH,temperature in field samples and the temperature and pressure in laboratory

样品编号	HCO3-浓度 c/(mmol·L-1)	pH值	温度 θ/℃	大气压 p/Pa
P地	4.1	7.28	19.5	995.9
P6	7.1	7.85	17.8	995.9
实验室	-	-	18.5	998.2

下载: 导出CSV 

| 显示表格

图  1  研究点位置示意图及地形图

Figure  1.  The location of the study site and the topographic map

下载: 全尺寸图片 幻灯片

1.2   水样的前处理

1.2.1   BaCl₂沉淀法

用于采集样品的600 mL聚乙烯瓶用稀盐酸浸泡24 h后，用超纯水(Mill-Q advantage A10超纯水机纯化，电阻率18.2 MΩ·cm，下同)洗净烘干，备用。在野外采集样品时先用水样冲洗样品瓶三次，加入200 mL左右过滤后的水样，再加入6 mL的2 mol/L NaOH溶液(用市售纯度 > 96%的分析纯NaOH粉末配制)至pH=12后继续加入过量的BaCl₂粉末(市售分析纯，纯度 > 99.5%)，最后再加满水样用parafilm封口膜密封。野外采集P_地和P₆样品(含平行样品)总共4份避光保存带回实验室处理。带回实验室的δ¹³C_DIC野外水样在低温中静置24 h后迅速过滤、烘干，得到不纯的BaCO₃样品。

1.2.2   医用无菌高密度聚乙烯瓶装样

市售医用无菌25 mL聚乙烯瓶，用稀盐酸浸泡24 h后，用超纯水洗净烘干，备用。在野外采集样品时先用水样润洗三次，采集的同时加入HgCl₂来淬灭水中微生物。装样时尽量避免样品瓶中有气泡，同时用parafilm封口膜密封瓶口。野外采集P_地和P₆样品(含平行样品)总共4份，避光保存，带回实验室低温保存以备测试。

1.2.3   GasBenchⅡ顶空样品瓶装样

GasBenchⅡ顶空样品瓶和1 mL医用注射器均用稀盐酸浸泡24 h后，用超纯水洗净烘干，备用。顶空样品瓶放入GasBenchⅡ样品盘中利用自动进样器加入脱水100%的磷酸(德国Merck公司)。为了避免排空时间较偏短导致的顶空瓶内残余空气中的CO₂和排空过程中外部空气少量回流对测试结果的影响^[21]，因此在加入脱水100%磷酸后用高纯He气(纯度 > 99.999%)吹扫540 s。在野外样品采集时先用水样润洗注射器三次后，取600 μL样品缓慢注射至顶空样品瓶内，注射样品时避免样品受瓶内气压影响喷至顶空样品瓶盖口导致测样时将水汽带入测试系统。野外采集P_地和P₆样品(含平行样品)总共4份带回实验室测试。顶空样品瓶带回实验室后在平衡18 h后开始测量。

1.3   水样溶解性无机碳同位素检测

碳同位素测试仪器为GasBenchⅡ-IRMS系统(美国Thermo Fisher公司)。GasBenchⅡ前处理装置包括：GC-PAL自动进样器；PoraPlot Q色谱柱(25 m×0.32 mm)；恒温样品盘(控制温度±0.1℃)。IRMS为MAT-253同位素质谱仪(美国Thermo Finnigan公司)。

δ¹³C分析计算公式为：

水样和粉末沉淀物样品的碳同位素测试结果均为相对于V-PDB标准，水样品和粉末沉淀物的碳同位素分析精度小于0.2‰。为了避免人为因素所产生的误差对实验结果造成影响，每个样品均有平行对比样品。平行对比样品测试数据均有较好的重现性。所有样品均由中国地质科学院岩溶地质研究所测试中心测定。

2.   结果与讨论

2.1   前处理方法测试结果对比

所有水样的DIC的碳同位素测试结果列于表 2。通过BaCl₂沉淀法测定不纯BaCO₃粉末的碳同位素，P_地为-14.30‰，P₆为-16.06‰，平行样品之间的测试结果偏差小于0.06‰。医用聚乙烯瓶采集水样测定的DIC的碳同位素P_地为-14.54‰，P₆为-16.39‰，平行样品之间的测试结果偏差小于0.03‰。顶空样品瓶采集水样的DIC的碳同位素P_地为-14.61‰，P₆为-16.41‰，平行样品之间的测试结果偏差小于0.06‰。

表  2  不同前处理方法的碳同位素分析结果

Table  2.  Analytical results of carbon isotope determined with different sample pretreatment methods

样品编号	δ13CV-PDB/‰
	BaCl2沉淀法	医用聚乙烯瓶采样	顶空样品瓶采样
P地	-14.30	-14.54	-14.67
P地-平行样品	-14.36	-14.52	-14.61
P6	-16.06	-16.39	-16.41
P6-平行样品	-16.05	-16.36	-16.39

下载: 导出CSV 

| 显示表格

本次研究所采集的两种岩溶区域富含DIC的水样的测试结果分布图如图 2所示。图 2数据分布指示了医用聚乙烯瓶和顶空样品瓶采集样品测试的DIC碳同位素值均有良好的重现性，并且两种前处理方法的重现性偏差均小于0.1‰；而BaCl₂沉淀法测定的碳同位素值明显偏正于后两种处理方法：地下河样品(P_地)碳同位素值偏正0.24‰，平行样品偏正0.26‰；洞穴滴水样品(P₆)的碳同位素值偏正0.33‰，平行样品偏正0.29‰，主要由于沉淀法前处理过程中只是将水体中碳同位素值相对偏重的HCO₃^-和CO₃^2-沉淀，水体中碳同位素值相对偏轻的游离CO₂脱气逸出。已有学者指出DIC中游离CO₂的碳同位素值明显偏负于水中剩余成分DIC的碳同位素值^[17]，为本次研究的BaCl₂沉淀偏正于其他的前处理方法提供了实验数据支持。

图  2  不同前处理方法测试的碳同位素分布图

Figure  2.  Carbon isotope distribution with different sample pretreatment methods

下载: 全尺寸图片 幻灯片

BaCl₂沉淀法中加入NaOH溶液使得HCO₃^-转化为CO₃^2-便于沉淀，这种方法使得水中DIC平衡方程式：CO₂+H₂O↔H₂CO₃↔H⁺+HCO₃^-↔2H⁺+CO₃^2-的动态平衡破坏，CO₂、HCO₃^-的溶解度发生改变，导致相对富集¹²C的溶解性游离CO₂逸出，或大气中相对富集¹³C的CO₂溶解进入水体，引起碳同位素分馏。

2.2   水体中DIC碳同位素最佳前处理方法的确定

通过对比三种前处理方法的碳同位素结果显示，BaCl₂沉淀法明显偏正于医用聚乙烯瓶和顶空样品瓶采集水样方法。由此可见，BaCl₂沉淀法引起了岩溶水体中相对富集¹²C的溶解性游离CO₂逸出，使得DIC的碳同位素测试结果偏正。由于医用聚乙烯瓶采集的水样带回实验室及时测试，并且野外水样采集环境的温度与实验室的温度仅仅相差1℃，大气压强的差异也仅为2.3 Pa，在短时间存于细微改变的外界环境中并未引起CO₂、HCO₃^-的溶解度发生变化。但是，如果实验室环境温度(例如夏季和冬季采集样的样品，温差超过10℃)和大气压强相对于野外样品采集环境发生较大改变时，肯定会导致医用聚乙烯瓶开盖的瞬间水中CO₂、HCO₃^-的溶解度改变,发生脱气作用，引起碳同位素分馏。由此可见，针对不同季节的样品采集，医用聚乙烯瓶采集水样并不是最好的前处理方法。

上述研究表明，医用聚乙烯瓶和顶空样品瓶采集样品都能满足DIC碳同位素的测试要求。顶空样品瓶在野外将样品在高纯氦气的背景下酸化，使DIC全部组分转化为CO₂气体，以免采集、运输、保存过程中由于人为操作等多种外界因素影响各组分的动态平衡导致的同位素分馏；并且，顶空样品瓶为硼硅玻璃制成，能防止游离和酸化生成的CO₂气体在采集、运输、保存过程中通过样品瓶壁“渗出”所引起的同位素分馏。可见，用顶空样品瓶采集水样的方法能够避免由于外界环境条件变化引起CO₂、HCO₃^-的溶解度发生改变，避免发生CO₂逸出或大气中CO₂溶解入水体引起碳同位素分馏，因此最佳的水样前处理方法是利用GasBenchⅡ顶空样品瓶直接产生CO₂气体，可为高分辨率监测水体中DIC的碳同位素变化提供精确的前处理方法。

3.   结语

通过对比传统BaCl₂沉淀法、医用无菌高密度聚乙烯瓶、GasBenchⅡ顶空样品瓶装样的DIC碳同位素前处理方法的实验数据，综合分析总结出以下结论。

(1) 对于pH值在7～8，HCO₃^-浓度在4～7 mmol/L的岩溶水体，传统BaCl₂沉淀法的测试结果明显偏正于医用无菌高密度聚乙烯瓶装样品、GasBenchⅡ顶空样品瓶直接产气的测试结果。地下河(P_地)和洞穴滴水(P₆)的碳同位素值分别偏正0.26‰、0.33‰，主要是由于水中相对富集¹²C，溶解性游离CO₂逸出，使得BaCl₂沉淀法测试结果相对偏于其他前处理方法。

(2) 应用GasBenchⅡ顶空样品瓶采集样品，能够避免外界环境条件变化导致CO₂、HCO₃^-的溶解度发生变化所引起的碳同位素分馏，确保测试结果的准确性。GasBenchⅡ顶空样品瓶直接产生CO₂气体，能够高分辨率监测水中DIC的碳同位素变化，精确地反映其时间演化规律。