Loading [MathJax]/jax/output/SVG/jax.js

石灰土对硫酸型酸雨缓冲过程模拟及碳汇效应研究

赵光帅, 黄奇波, 朱义年, 李腾芳, 普政功. 石灰土对硫酸型酸雨缓冲过程模拟及碳汇效应研究[J]. 中国岩溶, 2022, 41(5): 796-807. doi: 10.11932/karst20220504
引用本文: 赵光帅, 黄奇波, 朱义年, 李腾芳, 普政功. 石灰土对硫酸型酸雨缓冲过程模拟及碳汇效应研究[J]. 中国岩溶, 2022, 41(5): 796-807. doi: 10.11932/karst20220504
ZHAO Guangshuai, HUANG Qibo, ZHU Yinian, LI Tengfang, PU Zhenggong. Simulation of buffering process and carbon sink effect of lime soil on sulfuric acid rain[J]. Carsologica Sinica, 2022, 41(5): 796-807. doi: 10.11932/karst20220504
Citation: ZHAO Guangshuai, HUANG Qibo, ZHU Yinian, LI Tengfang, PU Zhenggong. Simulation of buffering process and carbon sink effect of lime soil on sulfuric acid rain[J]. Carsologica Sinica, 2022, 41(5): 796-807. doi: 10.11932/karst20220504

石灰土对硫酸型酸雨缓冲过程模拟及碳汇效应研究

  • 基金项目: 广西自然科学基金重点基金项目(2018GXNSFDA281036);中国地质调查项目(DD20221758)
详细信息
    作者简介: 赵光帅(1989-),男,博士研究生,助理研究员,主要从事岩溶碳循环与全球变化研究。E-mail: zhaoguangshuai@mail.cgs.gov.cn
    通讯作者: 黄奇波(1982-),男,博士,研究员,主要从事岩溶水文地质、岩溶碳循环与全球变化研究。E-mail: qbohuang0108@163.com
  • 中图分类号: X517;X142

Simulation of buffering process and carbon sink effect of lime soil on sulfuric acid rain

More Information
  • 硫酸型酸雨沉降至地表经石灰土缓冲后,参与碳酸盐岩溶蚀及对岩溶碳汇的影响尚不明确,严重制约了我国岩溶碳汇效应的准确评估。本研究通过设置不同土层厚度条件下pH 4.5的硫酸型酸雨淋滤实验,以明确石灰土对硫酸型酸雨的缓冲过程及关键控制因素。结果表明:石灰土对酸雨的缓冲作用主要发生在表层(10 cm),淋出液中Ca2+、Mg2+HCO3含量在淋溶初期均表现快速降低,当淋溶量(土壤水达饱和后)为1 020 mL时Ca2+、Mg2+HCO3淋失量趋于稳定,稳定淋失量分别为20 mg·L−1、6 mg·L−1、40 mg·L−1。淋出液中被酸雨H+交换出的Ca2+、Mg2+仅占很小一部分,土壤水溶性Ca2+、Mg2+是淋出液中Ca2+、Mg2+的主要部分,开放系统中,大气和土壤CO2溶于降雨形成H2CO3不仅增加碳汇,且H2CO3解离产生的H+交换土壤中交换态Ca2+、Mg2+,造成Ca2+、Mg2+的淋失量不容忽视。不同厚度石灰土中交换性钙镁可缓冲酸雨容量均大于土壤碳酸钙矿物可缓冲容量,前者是后者的1.17~1.59倍。相同酸度、同一降雨量(土壤水达饱和后)下土壤盐基离子参与酸雨缓冲产生的碳汇量约为碳酸盐矿物风化缓冲产生碳汇量的2.1倍,不同厚度(≥10 cm)石灰土产生的碳汇量大致相等。根据本次实验及桂林市降雨数据计算,桂林市质纯石灰岩风化残积土壤区(土层厚度≥10 cm),土壤盐基离子参与酸雨缓冲每年可产生0.59 ~0.93 mol·m−2的碳汇通量。

  • 长期以来,地质碳汇一直被公认为脱碳方案的一个重要部分,其中碳酸盐岩的碳酸溶蚀产生岩溶碳汇是地质碳汇的主要组成部分,然而近年来越来越多的学者发现硫酸广泛参与了流域岩石矿物的化学风化[1-4],硫酸参与碳酸盐岩溶蚀不消耗大气或土壤中的CO2,当环境条件改变时该过程产生的HCO3会因碳酸盐矿物沉淀而向大气净释放CO2,成为大气CO2[5-6]。我国西南岩溶区是岩溶碳汇的重要发生地,然而硫酸对碳酸盐岩侵蚀向大气释放的CO2通量相当于每年西南碳酸盐岩风化消耗CO2总通量的33%[5]。另一方面,岩溶区石灰土具有富钙偏碱的特征,pH接近中性或中性以上[7],阳离子交换量较高,对大气酸沉降具有较高的缓冲容量及阈值[8]。当土壤中交换性盐基离子不能充分交换H+时,H+才会溶蚀风化土壤或表层带碳酸盐岩[9]。因此,大气酸沉降的H+可能优先与土壤中交换性阳离子进行交换,而并没有参与碳酸盐岩的溶蚀[6],进而不能形成大量的因碳酸盐矿物沉淀而向大气释放CO2HCO3,产生减汇。

    土壤对外来离子有一定的中和缓冲能力[10-11],酸沉降进入土壤引起土壤内部一系列的化学反应,包括黏土矿物表面吸附的阳离子与酸雨H+交换、元素淋溶、黏土矿物风化、土壤酸化等,致使土壤可给性营养元素(钙、镁、钾、钠)流失,甚至释放或活化某些毒性元素(铝、镉等)[12-13]。耿增超等[14]依据土壤pH范围划分了碳酸盐缓冲体系、交换性盐基和可变电荷位点缓冲体系、铝(氢)氧化物缓冲体系3个土壤酸缓冲体系。廖柏寒等[15]认为酸雨进入土壤后按缓冲先后过程存在如下三个缓冲体系:阳离子交换缓冲、氢氧化铝水解缓冲和原生矿物风化缓冲。前人对酸雨的土壤缓冲性能已有较多的研究,且多运用模拟酸雨土柱淋溶实验探究土壤对酸雨的缓冲过程[15-18],而针对质纯石灰岩风化形成石灰土对酸雨缓冲性能的研究较少[19-20],由于石灰土有较强的酸雨缓冲能力,酸雨来源的H+经过石灰土层的缓冲后是否还具有参与碳酸盐岩溶蚀而产生减汇效应的研究并未开展,这不仅关系到全球岩溶碳汇效应计算结果的准确性和可信度,更是当前岩溶碳汇研究必须回答的科学问题。

    中国酸雨区是世界三大酸雨区之一[21-22],2003年以来中国酸雨区面积占国土面积比例为5.5%~15.6%,酸雨频率在50%以上的城市比例为7.9%~27.1%[23],其中强酸雨地区多出现在四川省东部、贵州省、广西壮族自治区北部[24]。本文以广西壮族自治区北部桂林市恭城县嘉会镇质纯厚层石灰岩风化残疾土壤为供试土壤,设置10~80 cm不同土层厚度条件下相同酸度(pH=4.5)硫酸型酸雨的淋滤实验,通过收集检测淋出液,以明确石灰土对硫酸型酸雨的缓冲效应及关键控制因素,揭示土壤对酸雨的缓冲机理和缓冲容量。研究结果有助于岩溶碳汇潜力的准确评估,确立我国岩溶碳汇的国际地位,对准确评价地质碳汇在碳中和中的作用意义重大。

    实验用土为桂林市恭城县嘉会镇绵羊脚村融县组(D3r)灰岩风化土壤,该处基岩为质纯微晶灰岩、生物粒屑灰岩夹少量白云质灰岩,在厚层灰岩风化形成土壤的地区具有很好的代表性。2021年7月分层采集供试土壤并用环刀法现场监测土壤容重和含水率,土壤剖面见照片1,采集土壤经风干研磨过100目样筛后送至自然资源部岩溶地质资源环境监督检测中心检测分析,土壤水溶性Ca2+、Mg2+和土壤交换性钙、镁采用电感耦合等离子光谱仪(赛默飞世尔iCAP 7600 ICP-OES)检测,土壤碳酸钙采用气量法检测,绝对偏差<2×10−3,检测结果见表1

    图 1.  供试土壤自然剖面
    Figure 1.  Natural profile of tested soil
    表 1.  供试土壤基本理化性质
    Table 1.  Physico-chemical properties of soil samples
    石灰土厚度/cm
    1020304050607080
    pH 6.62 6.66 6.51 6.63 6.62 6.52 6.54 6.92
    容重/g·cm−3 1.47 1.37 1.24 1.33 1.58 1.68 1.63 1.62
    含水率/% 23.2 26.6 24.7 27.9 29.7 27.0 25.6 27.4
    水溶性Ca2+/×10-6 5.75 11.70 8.97 10.90 8.89 3.10 11.80 8.02
    交换性钙/cmol·kg−1 7.4 4.7 4.2 4.6 6.4 4.8 5.0 5.3
    水溶性Mg2+/×10-9 3.48 3.89 4.77 6.42 5.90 4.69 5.90 6.02
    交换性镁/cmol·kg−1 2.4 1.9 1.7 2.1 2.4 2.5 2.7 2.9
    碳酸钙/×10-3 6.83 8.34 2.77 3.00 3.70 5.55 7.30 5.09
     | Show Table
    DownLoad: CSV

    以10 cm为一个土壤层,采集10~80 cm,8组层位土壤,将土壤晾晒风干研磨过20目筛子做实验用土。采样时现场分层监测土壤容重D(环刀法,未烘干状态)、含水率W(环刀法)。设置8组平行淋滤实验,每组实验石灰土厚度依次为10 cm、20 cm、···、80 cm,在实验室用超纯水和浓硫酸溶液配置pH为4.5的淋滤液,H+浓度3.16×10−5 mol·L−1。实验用淋滤柱为直径20 cm的圆柱形有机玻璃柱,底部设出水阀口,用于收集淋出液,出水孔上部放置3 cm厚石英砂,采用蠕动泵驱动淋滤实验,设置淋滤速率0.03 mL·s−1(蠕动泵转速单位为2),单位面积石灰土淋滤速率0.34 mL·h−1·cm−2,当淋滤液淹没土柱表层土壤约1 cm后关停蠕动泵,观察淋滤液液面位置,通过开关蠕动泵使淋滤液液面始终位于土柱表层土壤上方1 cm处。具体实验步骤如下:

    (1)由底层至顶层按层位将研磨过目后土壤按原状土顺序填入淋滤柱内,每层(10 cm)放置土壤的质量依如下公式计算:

    M=Dπ10210÷(1+W)

    式中:M为需添加的晒干过目后土壤质量;D为原始土壤容重;W为原始土壤含水率。

    (2)单层土壤样放置入淋滤柱后,添加M质量的纯水,M的计算公式:

    M=Dπ10210÷(1+W)W

    式中:M为每层土壤需添加的纯水质量;D为原始土壤容重;W为原始土壤含水率。

    (3)单层土壤添加完纯水后将土壤压平至该层土壤层位刻度(每层土壤层位厚10 cm),然后加入上层土壤,重复步骤(1)、(2)。使淋滤柱内土壤基本保持原状。

    (4)土壤添加完毕后,每组淋滤柱表层土壤上覆一层滤纸,并用塑料薄膜封住淋滤柱,静置24 h,然后开启蠕动泵,进行石灰土淋滤实验。实验装置见照片2

    图 2.  淋滤实验装置
    Figure 2.  Experimental equipment for leaching

    以检测水中Ca2+、Mg2+、K+、Na+HCO3SO24NO3、Cl 离子所需样品量(170 mL)为标准,在淋滤柱底部采集淋出液并记录采样时间,所有水样采样后放入4 ℃的保温箱中,3 d 内送至自然资源部岩溶地质资源环境监督检测中心检测。Ca2+、Mg2+、K+、Na+采用电感耦合等离子光谱仪(赛默飞世尔iCAP 7600 ICP-OES)检测;SO24NO3、Cl 采用离子色谱仪(万通883)检测;HCO3采用滴定法检测,相对标准偏差1.2%。

    实验共收集水样240组,10~80 cm厚度石灰土淋出液各30组(一组水样170 mL,淋溶总量5100 mL),淋出液主要水化学指标统计结果见表2

    表 2.  石灰土淋出液水化学分析结果 (mg·L−1)
    Table 2.  Hydrochemical analysis results of leachate from lime soil (mg·L−1)
    土层厚度/cmpHCa2+Mg2+K+Na+Cl
    10
    (n=30)
    最大值7.71127.7944.448.066.9060.6333.35547.5913.94
    最小值7.1515.685.722.970.5928.2919.833.991.48
    平均值7.5227.6810.003.691.3050.2624.3161.333.18
    标准差0.1727.079.871.331.329.523.11138.603.18
    变异系数2.3097.7998.6136.18101.7718.9412.79225.98100.30
    20
    (n=30)
    最大值7.4669.9716.6112.515.3644.4634.83249.8414.56
    最小值6.4217.084.406.210.7716.1717.144.861.90
    平均值7.0030.207.848.481.6131.6028.0181.154.36
    标准差0.2815.344.022.270.979.256.6186.152.77
    变异系数4.0050.7951.3226.8360.2029.2623.61106.1563.56
    30
    (n=30)
    最大值7.5355.7912.046.775.0544.4640.76175.9810.50
    最小值6.7319.044.644.100.9228.2920.444.771.90
    平均值7.1833.798.195.421.7939.2830.3477.234.62
    标准差0.2510.542.470.940.864.407.1754.161.86
    变异系数3.5131.1830.1617.3048.0911.2023.6370.1340.40
    40
    (n=30)
    最大值7.5748.209.304.825.2460.6329.72119.707.74
    最小值6.8821.555.293.361.2634.3616.5213.023.11
    平均值7.2329.607.103.992.0142.9821.7359.984.37
    标准差0.147.491.530.540.925.384.5141.551.04
    变异系数1.9525.3221.5413.4345.6812.5320.7769.2823.75
    50
    (n=30)
    最大值7.6050.549.544.376.0364.6730.32115.028.77
    最小值6.8219.554.742.811.1636.3816.775.723.10
    平均值7.2325.586.143.252.0844.5320.6941.804.14
    标准差0.226.911.220.411.086.243.2431.251.43
    变异系数3.1026.9919.8912.7051.9714.0215.6774.7634.54
    60
    (n=30)
    最大值7.5047.8410.433.727.2652.5533.96125.599.41
    最小值6.8519.245.162.281.3832.3416.603.523.60
    平均值7.1624.846.672.742.2342.8420.7839.904.63
    标准差0.186.281.270.361.206.593.8732.111.48
    变异系数2.4525.2919.1012.9753.7815.3918.6480.4832.06
    70
    (n=30)
    最大值7.4543.1811.032.625.8748.5032.28109.339.60
    最小值6.8917.635.561.631.3928.2918.603.303.60
    平均值7.2122.937.111.922.1938.7421.6640.084.92
    标准差0.155.431.270.220.957.603.1029.481.42
    变异系数2.0223.6817.8711.4643.2819.6214.3273.5628.88
    80
    (n=30)
    最大值7.3830.828.251.083.3844.4618.7862.978.53
    最小值6.7417.106.040.751.6028.2914.3316.364.88
    平均值7.0821.317.390.892.0235.3015.4145.095.81
    标准差0.192.890.660.070.384.491.0815.330.87
    变异系数2.6413.548.908.3318.8412.737.0133.9915.03
    注:n为样品数量。
     | Show Table
    DownLoad: CSV

    土壤中的盐基离子通常以可溶态、交换态或与矿物结合态存在,因此,盐基离子淋失的过程是与溶解、交换和水解作用相联系的[25]。由图1可知,10~80 cm土柱淋滤实验的淋溶初期,淋出液中Ca2+、Mg2+含量均表现快速降低,且均在淋溶量为1 020 mL后趋于稳定,Ca2+稳定流出量约为20 mg·L−1,Mg2+稳定流出量约为6 mg·L−1,不同厚度土柱淋出液同一离子表现相同的变化趋势说明酸雨的缓冲可能主要发生在10 cm(表层)土壤层,俞元春等[26]的研究结果也显示表层土对酸雨的缓冲性能远大于底层土。土壤每缓冲掉酸雨中一个当量的H+就会淋失一个当量的盐基离子[15],根据淋滤液H+浓度可知土壤中被交换淋失的[Ca2++Mg2+](当量)值为3.16×10−5 mol·L−1,该值仅为Ca2+、Mg2+稳定流出当量值(1.5×10−3 mol·L−1)的2.1%,说明淋出液中Ca2+、Mg2+可能主要为水溶性Ca2+、Mg2+淋失。以Ca2+、Mg2+稳定流出量和表1数据计算各土壤层位水溶性Ca2+、Mg2+含量及可参与淋溶量(表3),可知在本次实验设计的5 100 mL淋溶量内,10~30 cm土柱土壤水溶性Ca2+、Mg2+可能已基本淋溶完,并可能有碳酸盐矿物参与淋溶风化。

    图 1.  Ca2+、Mg2+淋溶动态变化
    Figure 1.  Leaching dynamic change of Ca2+ and Mg2+
    表 3.  不同厚度土壤Ca2+、Mg2+含量及可参与淋溶量
    Table 3.  Contents of Ca2+ and Mg2+ in soil with different thicknesses and the amount that can be involved in leaching
    石灰土厚度/cm
    10 20304050607080
    单层水溶性Ca2+/mg 21.48 39.78 28.03 35.61 34.04 12.84 48.12 32.04
    累积水溶性Ca2+/mg 21.48 61.26 89.29 124.90 158.94 171.78 219.90 251.94
    水溶性Ca2+可参与淋溶量/L 1.07 3.06 4.46 6.24 7.95 8.59 11.00 12.60
    单层水溶性Mg2+/mg 13.00 13.23 14.91 20.97 22.59 19.43 24.06 24.05
    累积水溶性Mg2+/mg 13.00 26.22 41.13 62.10 84.70 104.13 128.19 152.24
    水溶性Mg2+可参与淋溶量/L 2.17 4.37 6.86 10.35 14.12 17.36 21.37 25.37
     | Show Table
    DownLoad: CSV

    淋出液pH在6.42~7.71范围内变化,淋溶流失的无机碳主要以HCO3形式存在。由图2可知,10~30 cm土柱淋出液中HCO3含量随淋溶量的增加而增加,说明表层石灰土中无机碳受酸雨淋溶影响较为强烈。40~80 cm土柱淋出液中HCO3含量则表现先快速降低,在淋溶量为680 mL时趋于平稳,稳定流出量约在40 mg·L−1,以桂林市多年平均降雨量1 886 mm[27]计算可知研究区每年碳流失量约为14.8 g·m−2,该值约为陆地生态系统每年平均碳流失量[28]的2.2倍,可见岩溶区碳流失是陆地生态系统碳流失的主要发生场所,且以表层土壤碳淋失最为强烈。

    图 2.  HCO3SO24淋溶动态变化
    Figure 2.  Leaching dynamic change of HCO3 and SO24

    淋溶初始阶段不同厚度土柱淋出液中SO24含量随淋溶量增加均呈现快速下降,在淋溶量为1 020 mL时20~80 cm土柱淋出液SO24含量转为缓慢上升,该过程的临界点与Ca2+、Mg2+淋溶临界点相同,而10 cm土柱淋出液SO24含量微弱上升后继续下降,前文“Ca2+、Mg2+淋溶动态变化”分析可知10 cm土柱淋溶过程中可能有碳酸盐矿物参与淋溶,碳酸盐矿物淋溶风化后形成的Ca2+、Mg2+SO24结合形成溶解度低的硫酸盐淀积于土壤中,致使淋出液SO24含量下降。不同土柱淋溶初始阶段SO24含量快速下降可能与土壤中水溶性SO24随淋滤液流失有关,当淋溶量达某一阀值(约1 020 mL)时,土壤中水溶性SO24可能已基本淋失完,随着淋溶的继续进行,大量的SO24由于土壤胶体的专性吸附和静电吸附[14]滞留于土壤中,土壤胶体的吸附量逐渐达饱和,致使淋出液中SO24含量逐渐升高,且随土层厚度增加,淋出液中SO24含量增幅逐渐变缓。

    土壤中钙、镁的存在形态主要有矿物态、交换态、水溶态和少量的有机结合态,交换态钙约占土壤全钙量20%~30%,是土壤盐基总量的主要部分,镁在土壤中主要以矿物态存在,可占土壤全镁量的70%~90%[14]。由图3可知,淋出液中Ca2++Mg2+当量值远大于淋滤液H+交换出的盐基离子理论当量值(3.16×10−2 meq·L−1),这主要是因为淋出液中Ca2+、Mg2+除了被交换淋出的Ca2+、Mg2+,还含有大量的土壤水溶性Ca2+、Mg2+,当酸雨持续沉降时,土壤中原生或次生矿物风化消耗H+释放Ca2+、Mg2+,实现对酸沉降的缓冲[29]

    图 3.  淋出液[Ca2++Mg2+]与HCO3毫克当量关系
    Figure 3.  The equivalent relation between [Ca2++Mg2+] and HCO3 in the leaching solution

    水中HCO3来源主要为大气或土壤中CO2溶于水解离和碳酸盐矿物受质子侵蚀风化[30],反应方程如下:

    CO2+H2O=H2CO3=H++HCO3
    H++CaxMg(1x)CO3=Ca2+x+Mg2+(1x)+HCO3

    在石灰土的缓冲下,反应式1中生成的H+可能并不参与碳酸盐矿物的分解(反应式2),而优先参与交换土壤中的盐基离子[15],反应过程如下:

    \textemdash+H+\textemdashH++

    碳酸盐岩风化残积土中钙、镁含量可占其他化学成分的90%以上[7],土壤中盐基离子含量可估算为Ca2+、Mg2+含量。淋滤液为超纯水配制,无机碳含量可忽略不计,淋出液中HCO3主要为大气、土壤CO2溶于淋滤液解离生成(反应式1),该反应生成的H+交换土壤胶体中的盐基离子(反应式3),产生等当量的Ca2+、Mg2+,故淋出液中水溶性Ca2+、Mg2+当量值可用如下公式计算:

    [Ca2++Mg2+]水溶性=[Ca2++Mg2+]实测 [Ca2++Mg2+]H+[HCO3]实测 

    土壤中被H+交换出的Ca2+、Mg2+当量值为:

    [Ca2++Mg2+]=[Ca2++Mg2+]H+交换 +[HCO3]实测 

    由降水pH值计算可知淋出液中被H+交换淋失的Ca2+、Mg2+当量浓度为:

    [Ca2++Mg2+]H+=3.16×102meqL1

    根据公式(1)、(2)及表1表3数据计算土壤淋出液水溶性和交换性Ca2+、Mg2+当量值及淋失比,计算结果见表4

    表 4.  不同厚度土壤累积Ca2+、Mg2+含量及淋失比
    Table 4.  Accumulative content of Ca2+ and Mg2+ in soil with different thicknesses and its leaching-loss ratio
    石灰土厚度/cm
    1020304050607080
    累积土壤水溶性Ca2+/meq 1.07 3.06 4.46 6.24 7.95 8.59 11.00 12.60
    累积土壤水溶性Mg2+/meq 1.08 2.19 3.43 5.18 7.06 8.68 10.68 12.69
    累积土壤水溶性Ca2++Mg2+/meq 2.16 5.25 7.89 11.42 15.00 17.27 21.68 25.28
    淋出液水溶性Ca2++Mg2+/meq 6.95 8.23 8.65 6.81 5.25 5.43 5.47 5.46
    水溶性Ca2++Mg2+淋失比/% 321.76 156.76 109.63 59.63 35.00 31.44 25.23 21.60
    累积土壤交换性Ca2+/mg 11056 17448 22698 28709 38511 46466 54622 63091
    累积土壤交换性Ca2+/meq 553 872 1135 1435 1926 2323 2731 3155
    累积土壤交换性Mg2+/mg 2151 3702 4977 6623 8829 11315 13957 16738
    累积土壤交换性Mg2+/meq 179 308 415 552 736 943 1163 1395
    累积土壤交换性Ca2++Mg2+/meq 732 1180 1550 1987 2662 3266 3894 4550
    淋出液交换性Ca2++Mg2+/meq 4.36 2.80 3.44 3.75 3.88 3.74 3.40 3.11
    交换性Ca2++Mg2+淋失比/% 0.60 0.24 0.22 0.19 0.15 0.11 0.09 0.07
     | Show Table
    DownLoad: CSV

    表4可知,10~30 cm土柱土壤水溶性Ca2+、Mg2+淋失比均大于100%,说明淋出液中Ca2+、Mg2+有水溶性Ca2+、Mg2+、酸雨H+交换淋失Ca2+、Mg2+之外的其他来源,另外10~30 cm土柱淋出液中HCO3含量均随淋溶量升高而增加,说明淋滤液中H+或CO2溶解解离产生的H+可能与土壤中碳酸盐矿物反应生成HCO3,同时产生2倍当量的Ca2+、Mg2+(反应式2)。表层土壤特殊的物理结构可能会有优先流产生[31],淋滤液中H+不能很好的渗流至土壤中每个盐基离子交换点,虽然土壤中有足够可被H+交换的盐基离子,解离的H+仍有可能侵蚀土壤碳酸盐矿物产生一定的碳汇效应。淋出液中被酸雨H+交换出的Ca2+、Mg2+仅占很小一部分,土壤中水溶性Ca2+、Mg2+是淋出液中Ca2+、Mg2+的主要部分,在开放系统中,大气和土壤CO2溶于降雨形成H2CO3不仅增加碳汇,且H2CO3解离产生的H+交换土壤中交换态Ca2+、Mg2+,造成Ca2+、Mg2+的淋失量不容忽视。

    由反应式2、3可知,土壤受H+侵蚀产生的Ca2+、Mg2+HCO3存在一定的比例关系,若淋出液中[Ca2++Mg2+](当量值)显著偏离H+溶蚀碳酸盐曲线,则可能有大量的水溶性Ca2+、Mg2+参与淋溶。由图3可知,10 cm土柱有大量的水溶性Ca2+、Mg2+淋失,这可能因为表层土壤较易产生优先流,土壤对离子的吸附淀积作用较弱。当淋出液[HCO3]大于0.8 meq·L−1时,10 cm土柱淋出液离子当量关系接近H2CO3溶蚀碳酸盐曲线,说明在优先流作用下,H+的可交换点显著下降,造成以碳酸盐矿物溶蚀辅助缓冲淋滤液中H+

    交换性盐基阳离子、氢氧化铝和原生或次生矿物是土壤缓冲酸雨的主要反应物[15]。本次实验供试土壤缓冲酸雨反应物主要为土壤中交换性钙镁和碳酸盐矿物。土壤交换性钙镁的酸雨缓冲容量依据公式3及表5数据计算。若酸雨中k1mol H+k2mol CO2溶于降水产生的H+侵蚀碳酸盐矿物(反应式4、反应式5),依据(公式4)-(公式8)及表1数据可知碳酸盐矿物溶蚀产生的Ca2+当量值及可缓冲的酸雨容量,计算统计结果见表5

    表 5.  不同厚度石灰土可缓冲酸雨(pH=4.5)容量
    Table 5.  Capacity of buffer acid rain(pH=4.5) in different thicknesses of lime
    石灰土厚度/cm
    1020304050607080
    土壤交换性Ca2++Mg2+/meq 732 1 180 1 550 1 987 2 662 3 266 3 894 4 550
    淋出液交换性Ca2++Mg2+/meq 4.36 2.80 3.44 3.75 3.88 3.74 3.40 3.11
    H交换态/mm(pH=4.5) 27 269 68 449 73 184 86 061 111 434 141 836 186 019 237 625
    土壤碳酸钙以钙计/g 10.20 21.55 25.01 28.93 34.60 43.79 55.70 63.84
    溶蚀矿物态钙/g 0.09 0.06 0.07 0.08 0.08 0.07 0.07 0.06
    H矿物态/mm(pH=4.5) 18408 58 336 58 030 58 735 70 247 101 606 129 240 172 815
    H交换态/ H矿物态 1.48 1.17 1.26 1.47 1.59 1.40 1.44 1.38
     | Show Table
    DownLoad: CSV
    H交换态 =10×5100×[Ca2++Mg2+]土壤102× π ×[Ca2++Mg2+]淋出液

    式中:H交换态为土壤交换性钙镁可缓冲pH4.5的酸雨容量(mm);5 100为本次实验设计的总淋溶量(mL);[Ca2++Mg2+]土壤为土壤交换性钙镁当量值(meq);[Ca2++Mg2+]淋出液为淋出液交换性钙镁当量值(meq)。

    k2CO2+k2H2O=k2H2CO3==k2H++k2HCO3
    (k1+k2)H++(k1+k2)CaxMg(1x)CO3=(k1+k2)Ca2+x+(k1+k2)Mg2+(1x)+(k1+k2)HCO3

    上述反应产生的HCO3满足如下关系:

    [HCO3]实测 =k1+2k2
    k1=3.16×105 molL1(pH)
    [Ca2+]溶蚀 =2k1+2k2

    土壤原生或次生碳酸钙可缓冲酸雨容量为:

    H矿物态 =10×5100 ×MCa土壤102× π ×MCa溶蚀
    MCa 溶蚀 =40×1/2[Ca2+]溶蚀 ×1/1000

    式中:H矿物态为土壤碳酸钙矿物可缓冲pH=4.5的酸雨容量(mm);5100为本次实验设计总淋溶量(mL);MCa土壤为土壤碳酸钙(以钙计)含量(g);MCa溶蚀为土壤碳酸钙(以钙计)的溶蚀量(g)。

    表5可知,10~80 cm厚石灰土中交换性钙镁可缓冲27 269~23 7625 mm pH为4.5的酸性降水,土壤碳酸钙矿物可缓冲18 408~172 815 mm pH为4.5的酸性降水,随土层厚度增加土壤交换性钙镁和碳酸钙可缓冲酸性降水量逐渐增加。对比可知,土壤交换性钙镁对酸雨的缓冲容量是碳酸盐矿物(碳酸钙)缓冲容量的1.17~1.59倍,这一方面因为土壤中交换性钙镁含量远大于土壤碳酸钙含量,另一方面土壤交换性钙镁优先参与缓冲酸雨中H+,但由于原状土壤的特殊物理结构,淋滤液不能均匀的渗流到土壤中的离子交换点,致使交换性钙镁淋失不均匀,局部伴随有碳酸钙淋溶风化。实验结果显示土壤交换性盐基离子具有巨大的酸雨缓冲能力,一定强度的酸性降水并不能很快将交换性盐基离子全部消耗掉,甚至需要持续几十年的酸性降水才能完全消耗掉土壤中交换性盐基离子。

    然而天然条件下石灰土对酸性降水的缓冲能力往往远大于实验计算的结果[32],由于天然条件的复杂性,本研究对石灰土缓冲酸雨容量的计算设置了几个限制条件:(1)只考虑酸雨中H+输入带来的交换缓冲反应,未考虑天然降水中H+之外的其他成分的影响;(2)输入的H+与交换性盐基离子的交换是完全的;(3)未考虑土表植被凋落物、土壤养分循环、土壤生物作用、人为输入物质等的影响。

    酸雨沉降经石灰土缓冲后,H+可能全部参与交换土壤中盐基离子,也可能侵蚀部分碳酸盐矿物甚至全部与碳酸盐矿物反应。由反应式1、3可知若酸雨中a mmol H+b mmol CO2溶于降水解离产生的H+全部参与土壤盐基离子交换,则会产生by mmol的碳汇,计算公式如下:

    by=30i=1[(HCO3)实测 ×0.17]

    a mmol H+b mmol CO2溶解解离产生的H+全部与碳酸盐矿物反应,依据反应式4-5及公式4可知产生的碳汇量为:

    bt=30i=1{1/2[(HCO3)实测 a]×0.17}

    式中:bybt为两种酸雨缓冲机制产生的碳汇量(mmol);(HCO3)实测为淋出液HCO3浓度(mmol·L−1);a为3.16×10−2 mmol·L−1 (依降水pH 计算),每个土柱共收集30组淋出液,每组淋出液0.17 L。

    依据公式9和公式10计算土壤盐基离子、碳酸盐矿物分别参与酸雨缓冲作用产生的碳汇量(表6),可知土壤盐基离子参与酸雨缓冲产生的碳汇量约为碳酸盐矿物风化缓冲产生碳汇量的2.1倍。相同酸度、同一降雨量(土壤水达饱和后)下不同厚度石灰土产生的碳汇量大致相等,说明石灰土对酸雨的缓冲作用主要发生在表层,相同降雨条件下石灰土的存在能显著提高碳汇量。

    表 6.  酸雨沉降下不同厚度石灰土产生的碳汇通量
    Table 6.  The carbon sink flux of lime soil with different thicknesses caused by acid rain deposition
    石灰土厚度/cm
    1020304050607080
    碳汇量by/mmol 4.20 2.64 3.28 3.59 3.72 3.58 3.24 2.95
    碳汇量bt/mmol 2.02 1.24 1.56 1.72 1.78 1.71 1.54 1.40
    碳汇量by/碳汇量bt 2.08 2.13 2.10 2.09 2.09 2.09 2.10 2.11
    碳汇量CY/mol·m−2·y−1 0.93 0.59 0.73 0.80 0.83 0.79 0.72 0.65
    碳汇量CT/mol·m−2·y−1 0.45 0.28 0.35 0.38 0.39 0.38 0.34 0.31
     | Show Table
    DownLoad: CSV

    研究区(桂林市)多年平均降雨量1 886 mm,降雨pH年均值范围4.85~5.23,酸雨频率范围42.5%~74.9%[24],本文取降雨pH=4.5、酸雨频率60%、降雨量1 886 mm计算石灰土缓冲下,酸雨沉降引起的碳汇量。

    土壤盐基离子(Y)、碳酸盐矿物(T)分别参与酸雨缓冲产生的碳汇量计算公式如下:

    CY=h×δ×π×r2×by×1015100×104π×r2×103=h×δ×by5100
    CT=h×δ×π×r2×bt×1015100×104π×r2×103=h×δ×bt5100

    式中:CYCT分别为土壤盐基离子、碳酸盐矿物参与酸雨缓冲产生的碳汇量(mol·m−2·y−1);h为降雨量(mm);δ为酸雨频率;r为淋滤柱半径(cm);bybt分别为本次实验土壤盐基离子、碳酸盐矿物参与酸雨缓冲产生的碳汇量(mmol);5 100为本次实验设计的总淋溶量(mL)。

    根据公式11和公式12计算研究区石灰土缓冲下酸雨沉降产生的碳汇量(表6),可知桂林市质纯石灰岩风化残积土壤区(土层厚度≥10 cm),土壤盐基离子参与酸雨缓冲每年可产生0.59~0.93 mol·m−2的碳汇通量。当土壤盐基离子淋失殆尽或无土层覆盖,原生或次生碳酸盐矿物参与酸雨缓冲每年可产生0.28~0.45 mol·m−2的碳汇通量。

    (1)石灰土对酸雨的缓冲作用主要发生在表层(10 cm),表层土壤特殊的物理结构可能会引起优先流,并致使H+的可交换点显著下降,造成以碳酸盐矿物溶蚀辅助缓冲淋滤液中H+。淋出液中被酸雨H+交换出的Ca2+、Mg2+仅占很小一部分,土壤中水溶性Ca2+、Mg2+是淋出液中Ca2+、Mg2+的主要部分,在开放系统中,大气和土壤CO2溶于降雨形成H2CO3不仅增加碳汇,且造成土壤中交换态Ca2+、Mg2+一定程度的流失;

    (2)不同厚度石灰土中交换性钙镁可缓冲酸雨容量均大于土壤碳酸钙矿物可缓冲容量,前者是后者的1.17~1.59倍,随土层厚度增加土壤交换性钙镁和碳酸钙可缓冲酸性降水量逐渐增加。土壤交换性盐基离子具有巨大的酸雨缓冲能力,一定强度的酸性降水并不能很快将交换性盐基离子消耗掉,甚至需要持续几十年的酸性降水才能完全消耗掉土壤中交换性盐基离子;

    (3)相同酸度、同一降雨量(土壤水达饱和后)下土壤盐基离子参与酸雨缓冲产生的碳汇量约为碳酸盐矿物风化缓冲产生碳汇量的2.1倍,不同厚度(≥10 cm)石灰土产生的碳汇量差值较小。根据本次实验及桂林市降雨数据计算,桂林市质纯石灰岩风化残积土壤区(土层厚度≥10 cm),土壤盐基离子参与酸雨缓冲每年可产生0.59~0.93 mol·m−2的碳汇通量。当土壤盐基离子淋失殆尽或无土层覆盖,原生或次生碳酸盐矿物参与酸雨缓冲每年可产生0.28~0.45 mol·m−2的碳汇通量。

  • 图 1  供试土壤自然剖面

    Figure 1. 

    图 2  淋滤实验装置

    Figure 2. 

    图 1  Ca2+、Mg2+淋溶动态变化

    Figure 1. 

    图 2  HCO3SO24淋溶动态变化

    Figure 2. 

    图 3  淋出液[Ca2++Mg2+]与HCO3毫克当量关系

    Figure 3. 

    表 1  供试土壤基本理化性质

    Table 1.  Physico-chemical properties of soil samples

    石灰土厚度/cm
    1020304050607080
    pH 6.62 6.66 6.51 6.63 6.62 6.52 6.54 6.92
    容重/g·cm−3 1.47 1.37 1.24 1.33 1.58 1.68 1.63 1.62
    含水率/% 23.2 26.6 24.7 27.9 29.7 27.0 25.6 27.4
    水溶性Ca2+/×10-6 5.75 11.70 8.97 10.90 8.89 3.10 11.80 8.02
    交换性钙/cmol·kg−1 7.4 4.7 4.2 4.6 6.4 4.8 5.0 5.3
    水溶性Mg2+/×10-9 3.48 3.89 4.77 6.42 5.90 4.69 5.90 6.02
    交换性镁/cmol·kg−1 2.4 1.9 1.7 2.1 2.4 2.5 2.7 2.9
    碳酸钙/×10-3 6.83 8.34 2.77 3.00 3.70 5.55 7.30 5.09
    下载: 导出CSV

    表 2  石灰土淋出液水化学分析结果 (mg·L−1)

    Table 2.  Hydrochemical analysis results of leachate from lime soil (mg·L−1)

    土层厚度/cmpHCa2+Mg2+K+Na+Cl
    10
    (n=30)
    最大值7.71127.7944.448.066.9060.6333.35547.5913.94
    最小值7.1515.685.722.970.5928.2919.833.991.48
    平均值7.5227.6810.003.691.3050.2624.3161.333.18
    标准差0.1727.079.871.331.329.523.11138.603.18
    变异系数2.3097.7998.6136.18101.7718.9412.79225.98100.30
    20
    (n=30)
    最大值7.4669.9716.6112.515.3644.4634.83249.8414.56
    最小值6.4217.084.406.210.7716.1717.144.861.90
    平均值7.0030.207.848.481.6131.6028.0181.154.36
    标准差0.2815.344.022.270.979.256.6186.152.77
    变异系数4.0050.7951.3226.8360.2029.2623.61106.1563.56
    30
    (n=30)
    最大值7.5355.7912.046.775.0544.4640.76175.9810.50
    最小值6.7319.044.644.100.9228.2920.444.771.90
    平均值7.1833.798.195.421.7939.2830.3477.234.62
    标准差0.2510.542.470.940.864.407.1754.161.86
    变异系数3.5131.1830.1617.3048.0911.2023.6370.1340.40
    40
    (n=30)
    最大值7.5748.209.304.825.2460.6329.72119.707.74
    最小值6.8821.555.293.361.2634.3616.5213.023.11
    平均值7.2329.607.103.992.0142.9821.7359.984.37
    标准差0.147.491.530.540.925.384.5141.551.04
    变异系数1.9525.3221.5413.4345.6812.5320.7769.2823.75
    50
    (n=30)
    最大值7.6050.549.544.376.0364.6730.32115.028.77
    最小值6.8219.554.742.811.1636.3816.775.723.10
    平均值7.2325.586.143.252.0844.5320.6941.804.14
    标准差0.226.911.220.411.086.243.2431.251.43
    变异系数3.1026.9919.8912.7051.9714.0215.6774.7634.54
    60
    (n=30)
    最大值7.5047.8410.433.727.2652.5533.96125.599.41
    最小值6.8519.245.162.281.3832.3416.603.523.60
    平均值7.1624.846.672.742.2342.8420.7839.904.63
    标准差0.186.281.270.361.206.593.8732.111.48
    变异系数2.4525.2919.1012.9753.7815.3918.6480.4832.06
    70
    (n=30)
    最大值7.4543.1811.032.625.8748.5032.28109.339.60
    最小值6.8917.635.561.631.3928.2918.603.303.60
    平均值7.2122.937.111.922.1938.7421.6640.084.92
    标准差0.155.431.270.220.957.603.1029.481.42
    变异系数2.0223.6817.8711.4643.2819.6214.3273.5628.88
    80
    (n=30)
    最大值7.3830.828.251.083.3844.4618.7862.978.53
    最小值6.7417.106.040.751.6028.2914.3316.364.88
    平均值7.0821.317.390.892.0235.3015.4145.095.81
    标准差0.192.890.660.070.384.491.0815.330.87
    变异系数2.6413.548.908.3318.8412.737.0133.9915.03
    注:n为样品数量。
    下载: 导出CSV

    表 3  不同厚度土壤Ca2+、Mg2+含量及可参与淋溶量

    Table 3.  Contents of Ca2+ and Mg2+ in soil with different thicknesses and the amount that can be involved in leaching

    石灰土厚度/cm
    10 20304050607080
    单层水溶性Ca2+/mg 21.48 39.78 28.03 35.61 34.04 12.84 48.12 32.04
    累积水溶性Ca2+/mg 21.48 61.26 89.29 124.90 158.94 171.78 219.90 251.94
    水溶性Ca2+可参与淋溶量/L 1.07 3.06 4.46 6.24 7.95 8.59 11.00 12.60
    单层水溶性Mg2+/mg 13.00 13.23 14.91 20.97 22.59 19.43 24.06 24.05
    累积水溶性Mg2+/mg 13.00 26.22 41.13 62.10 84.70 104.13 128.19 152.24
    水溶性Mg2+可参与淋溶量/L 2.17 4.37 6.86 10.35 14.12 17.36 21.37 25.37
    下载: 导出CSV

    表 4  不同厚度土壤累积Ca2+、Mg2+含量及淋失比

    Table 4.  Accumulative content of Ca2+ and Mg2+ in soil with different thicknesses and its leaching-loss ratio

    石灰土厚度/cm
    1020304050607080
    累积土壤水溶性Ca2+/meq 1.07 3.06 4.46 6.24 7.95 8.59 11.00 12.60
    累积土壤水溶性Mg2+/meq 1.08 2.19 3.43 5.18 7.06 8.68 10.68 12.69
    累积土壤水溶性Ca2++Mg2+/meq 2.16 5.25 7.89 11.42 15.00 17.27 21.68 25.28
    淋出液水溶性Ca2++Mg2+/meq 6.95 8.23 8.65 6.81 5.25 5.43 5.47 5.46
    水溶性Ca2++Mg2+淋失比/% 321.76 156.76 109.63 59.63 35.00 31.44 25.23 21.60
    累积土壤交换性Ca2+/mg 11056 17448 22698 28709 38511 46466 54622 63091
    累积土壤交换性Ca2+/meq 553 872 1135 1435 1926 2323 2731 3155
    累积土壤交换性Mg2+/mg 2151 3702 4977 6623 8829 11315 13957 16738
    累积土壤交换性Mg2+/meq 179 308 415 552 736 943 1163 1395
    累积土壤交换性Ca2++Mg2+/meq 732 1180 1550 1987 2662 3266 3894 4550
    淋出液交换性Ca2++Mg2+/meq 4.36 2.80 3.44 3.75 3.88 3.74 3.40 3.11
    交换性Ca2++Mg2+淋失比/% 0.60 0.24 0.22 0.19 0.15 0.11 0.09 0.07
    下载: 导出CSV

    表 5  不同厚度石灰土可缓冲酸雨(pH=4.5)容量

    Table 5.  Capacity of buffer acid rain(pH=4.5) in different thicknesses of lime

    石灰土厚度/cm
    1020304050607080
    土壤交换性Ca2++Mg2+/meq 732 1 180 1 550 1 987 2 662 3 266 3 894 4 550
    淋出液交换性Ca2++Mg2+/meq 4.36 2.80 3.44 3.75 3.88 3.74 3.40 3.11
    H交换态/mm(pH=4.5) 27 269 68 449 73 184 86 061 111 434 141 836 186 019 237 625
    土壤碳酸钙以钙计/g 10.20 21.55 25.01 28.93 34.60 43.79 55.70 63.84
    溶蚀矿物态钙/g 0.09 0.06 0.07 0.08 0.08 0.07 0.07 0.06
    H矿物态/mm(pH=4.5) 18408 58 336 58 030 58 735 70 247 101 606 129 240 172 815
    H交换态/ H矿物态 1.48 1.17 1.26 1.47 1.59 1.40 1.44 1.38
    下载: 导出CSV

    表 6  酸雨沉降下不同厚度石灰土产生的碳汇通量

    Table 6.  The carbon sink flux of lime soil with different thicknesses caused by acid rain deposition

    石灰土厚度/cm
    1020304050607080
    碳汇量by/mmol 4.20 2.64 3.28 3.59 3.72 3.58 3.24 2.95
    碳汇量bt/mmol 2.02 1.24 1.56 1.72 1.78 1.71 1.54 1.40
    碳汇量by/碳汇量bt 2.08 2.13 2.10 2.09 2.09 2.09 2.10 2.11
    碳汇量CY/mol·m−2·y−1 0.93 0.59 0.73 0.80 0.83 0.79 0.72 0.65
    碳汇量CT/mol·m−2·y−1 0.45 0.28 0.35 0.38 0.39 0.38 0.34 0.31
    下载: 导出CSV
  • [1]

    Han G L, Liu C Q. Water geochemistry controlled by carbonate dissolution: a study of the river waters draining karst-dominated terrain, Guizhou Province, China[J]. Chemical Geology, 2004, 204(1-2):1-21. doi: 10.1016/j.chemgeo.2003.09.009

    [2]

    Ding H, Lang Y C, Liu C Q, et al. Chemical characteristics and δ34S-SO42− of acid rain: Anthropogenic sulfate deposition and its impacts on CO2 consumption in the rural karst area of southwest China[J]. Geochemical Journal, 2013, 47(6):625-638. doi: 10.2343/geochemj.2.0293

    [3]

    Huang Q B, Qin X Q, Liu P Y, et al. Impact of sulfuric and nitric acids on carbonate dissolution, and the sociated deficit of CO2 uptake in the upper-middle reaches of the Wujiang River, China[J]. Journal of Contaminant Hydrology, 2017, 203:18-27. doi: 10.1016/j.jconhyd.2017.05.006

    [4]

    于奭. 西江流域河流化学风化及其碳汇效应研究[D]. 武汉: 中国地质大学, 2019.

    YU Shi. Study on chemical weathering and carbon sink in Xijiang basin[D]. Wuhan: China University of Geosciences, 2019.

    [5]

    刘丛强, 蒋颖魁, 陶发祥, 郎赟超, 李思亮. 西南喀斯特流域碳酸盐岩的硫酸侵蚀与碳循环[J]. 地球化学, 2008, 37(4):404-414. doi: 10.3321/j.issn:0379-1726.2008.04.014

    LIU Congqiang, JIANG Yingkui, TAO Faxiang, LANG Yunchao, LI Siliang. Chemical weathering of carbonate rocks by sulfuric acid and the carbon cycling in Southwest China[J]. Geochimica, 2008, 37(4):404-414. doi: 10.3321/j.issn:0379-1726.2008.04.014

    [6]

    黄奇波, 覃小群, 程瑞瑞, 李腾芳, 刘朋雨. 硫酸型酸雨参与碳酸盐岩溶蚀的研究进展[J]. 中国岩溶, 2019, 38(2):149-156.

    HUANG Qibo, QIN Xiaoqun, CHENG Ruirui, LI Tengfang, LIU Pengyu. Research progress of sulfuric acid rain participating in the dissolution of carbonate rocks[J]. Carsologica Sinica, 2019, 38(2):149-156.

    [7]

    袁道先, 蔡桂鸿. 岩溶环境学[M]. 重庆: 重庆出版社, 1988: 24-27

    YUAN Daoxian, CAI Guihong. The science of karst environment[M]. Chongqing: Chongqing Publishing House, 1988: 24-27

    [8]

    Matschullat J, Andreae H, Lessmann D, et al. Catchment acidification-from the top down[J]. Environmental Pollution, 1992, 77(2-3):143-150. doi: 10.1016/0269-7491(92)90070-Q

    [9]

    黄芬, 肖琼, 尹伟璐, 胡刚. 岩溶系统中土壤氮肥施用对岩溶碳汇的影响[J]. 中国岩溶, 2014, 33(4):405-411. doi: 10.11932/karst20140403

    HUANG Fen, XIAO Qiong, YIN Weilu, HU Gang. The effects of using N-fertilizers in soil on karst carbon sink in karst system[J]. Carsologica Sinica, 2014, 33(4):405-411. doi: 10.11932/karst20140403

    [10]

    Hartikainen H. Soil response to acid percolation: Acid-base buffering and cation leaching[J]. Journal of Environmental Quality, 1996, 25(4):638-645.

    [11]

    廖柏寒, 蒋青. 酸沉降与我国南方森林土壤的酸化[J]. 农业环境保护, 2002, 21(2):110-114.

    LIAO Baihan, JIANG Qing. Acid deposition and acidification of forest soils in southern China[J]. Agro-environmental Protection, 2002, 21(2):110-114.

    [12]

    岑慧贤, 王树功, 仇荣亮, 马灵芳. 模拟酸雨对土壤盐基离子的淋溶释放影响[J]. 环境污染与防治, 2001, 23(1):13-16. doi: 10.3969/j.issn.1001-3865.2001.01.005

    CEN Huixian, WANG Shugong, QIU Rongliang, MA Lingfang. Effect of simulating acid rain on cation release of some soils[J]. Environmental Pollution & Control, 2001, 23(1):13-16. doi: 10.3969/j.issn.1001-3865.2001.01.005

    [13]

    杨忠芳, 余涛, 唐金荣, 朱翠娟, 宗思锋, 张娇, 张建新, 申志军. 湖南洞庭湖地区土壤酸化特征及机理研究[J]. 地学前缘, 2006, 13(1):105-112.

    YANG Zhongfang, YU Tao, TANG Jinrong, ZHU Cuijuan, ZONG Sifeng, ZHANG Jiao, ZHANG Jianxin, SHEN Zhijun. A study of the characteristics and mechanisms of soil acidification in the Dongting Lake region in Hunan Province, South China[J]. Earth Science Frontiers, 2006, 13(1):105-112.

    [14]

    耿增超, 贾宏涛. 土壤学[M]. 第二版. 北京: 科学出版社, 2020.

    GENG Zengchao, JIA Hongtao. Soil Science[M]. Second Edition. Beijing: Science Press, 2020.

    [15]

    廖柏寒, 李长生. 土壤对酸沉降缓冲机制探讨[J]. 环境科学, 1989, 10(1):30-34. doi: 10.3321/j.issn:0250-3301.1989.01.001

    LIAO Baihan, LI Changsheng. Study on buffering mechanism of acid precipitation in soil[J]. Environmental Science, 1989, 10(1):30-34. doi: 10.3321/j.issn:0250-3301.1989.01.001

    [16]

    Liu K H, Mansell R S, Rhue R D. Cation removal during application of acid solutions into air-dry soil columns[J]. Soil Science Society of America Journal, 1990, 54(6):1747-1753. doi: 10.2136/sssaj1990.03615995005400060040x

    [17]

    凌大炯, 章家恩, 黄倩春, 韩维栋, 欧阳颖. 模拟酸雨对砖红壤盐基离子迁移和释放的影响[J]. 土壤学报, 2007, 44(3):444-450. doi: 10.3321/j.issn:0564-3929.2007.03.010

    LING Dajiong, ZHANG Jiaen, HUANG Qianchun, HAN Weidong, OU Yangying. Influences of simulated acid rain on leaching and release of base ions in latosol[J]. Acta Pedologica Sinica, 2007, 44(3):444-450. doi: 10.3321/j.issn:0564-3929.2007.03.010

    [18]

    于俊霞, 焦燕, 杨文柱, 刘立家, 宋春妮, 于亚泽. 外源盐对盐碱土壤无机碳淋溶特征的影响[J]. 环境科学学报, 2021, 41(6):2358-2368. doi: 10.13671/j.hjkxxb.2020.0500

    YU Junxia, JIAO Yan, YANG Wenzhu, LIU Lijia, SONG Chunni, YU Yaze. Influence of exogenous salt on inorganic carbon leaching in saline-alkali soil[J]. Acta Scientiae Circumstantiae, 2021, 41(6):2358-2368. doi: 10.13671/j.hjkxxb.2020.0500

    [19]

    李春龙, 赵家梅, 龙偲, 陈中吉, 周运超, 张春来. 模拟酸雨条件下石灰土-碳酸盐岩体系的碳汇效应[J]. 中国岩溶, 2014, 33(1):51-56. doi: 10.3969/j.issn.1001-4810.2014.01.008

    LI Chunlong, ZHAO Jiamei, LONG Si, CHEN Zhongji, ZHOU Yunchao, ZHANG Chunlai. Carbon sink effect of simulated acid rain in lime soil and carbonate rock system[J]. Carsologica Sinica, 2014, 33(1):51-56. doi: 10.3969/j.issn.1001-4810.2014.01.008

    [20]

    刘炜, 周运超, 张春来. 石灰土盐基离子迁移对模拟酸雨的响应[J]. 中国岩溶, 2018, 37(3):336-342.

    LIU Wei, ZHOU Yunchao, ZHANG Chunlai. Response of base cations migration of lime soil to simulated acid rain[J]. Carsologica Sinica, 2018, 37(3):336-342.

    [21]

    吴春笃, 戴竞, 阿琼, 陈诗龙, 解清杰. 酸雨的研究现状及新兴趋式的可视化分析[J]. 安全与环境学报, 2019, 19(1):344-353. doi: 10.13637/j.issn.1009-6094.2019.01.051

    WU Chundu, DAI Jing, A Qiong, CHEN Shilong, XIE Qingjie. Visualization analysis of the research status and emerging trends of the acid rain application[J]. Journal of Safety and Environment, 2019, 19(1):344-353. doi: 10.13637/j.issn.1009-6094.2019.01.051

    [22]

    Pan Y, Birdsey R A, Fang J, et al. A large and persistent carbon sink in the world’s forests[J]. Science, 2011, 333(6045):988-993. doi: 10.1126/science.1201609

    [23]

    解淑艳, 王胜杰, 于洋, 刘丽, 张凤英. 2003-2008年全国酸雨状况变化趋势研究[J]. 中国环境监测, 2020, 36(4):80-88.

    XIE Shuyan, WANG Shengjie, YU Yang, LIU Li, ZHANG Fengying. Analysis on acid rain status and its change trends in China from 2003 to 2008[J]. Environmental Monitoring in China, 2020, 36(4):80-88.

    [24]

    李陵, 黎泳珊, 彭良, 周艺, 柴发合, 莫招育, 陈志明, 李红. 桂林市酸雨污染特征及防治对策研究[J]. 环境科学研究, 2020, 33(6):1393-1401. doi: 10.13198/j.issn.1001-6929.2020.03.11

    LI Ling, LI Yongshan, PENG Liang, ZHOU Yi, CHAI Fahe, MO Zhaoyu, CHEN Zhiming, LI Hong. Pollution characteristics and control countermeasures of acid rain in Guilin[J]. Research of Environmental Sciences, 2020, 33(6):1393-1401. doi: 10.13198/j.issn.1001-6929.2020.03.11

    [25]

    陈怀满. 环境土壤学[M]. 第三版. 北京: 科学出版社, 2018.

    CHEN Huaiman. Environmental soil science[M]. Third Edition. Beijing: Science Press, 2018.

    [26]

    俞元春, 丁爱芳, 胡笳, 孟磊. 模拟酸雨对土壤酸化和盐基迁移的影响[J]. 南京林业大学学报, 2001, 25(2):39-42.

    YU Yuanchun, DING Aifang, HU Jia, MENG Lei. Effects of simulated acid rain on soil acidification and base ions transplant[J]. Journal of Nanjing Forestry University, 2001, 25(2):39-42.

    [27]

    黄奇波, 吴华英, 程瑞瑞, 李腾芳, 罗飞, 赵光帅, 李小盼. 桂林岩溶区石灰土壤对酸雨缓冲作用的观测及其对岩溶碳汇的指示意义[J/OL]. 地球学报, https://kns.cnki.net/kcms/detail/11.3474.p.20220121.1039.002.html.

    HUANG Qibo, WU Huaying, CHENG Ruirui, LI Tengfang, LUO Fei, ZHAO Guangshuai, LI Xiaopan. Buffering effect of lime soil on acid rain and its influence on the evaluation of the karst carbon sink effect[J/OL]. Acta Geoscientica Sinica, https://kns.cnki.net/kcms/detail/11.3474.p.20220121.1039.002.html.

    [28]

    Schulze E D, Luyssaert S, Ciais P, et al. Importance of methane and nitrous oxide for Europe’s terrestrial greenhouse-gas balance[J]. Nature Geoscience, 2009, 2(12):842-850. doi: 10.1038/ngeo686

    [29]

    廖柏寒, 戴昭华. 土壤对酸沉降的缓冲能力与土壤矿物的风化特征[J]. 环境科学学报, 1991, 11(4):425-431. doi: 10.13671/j.hjkxxb.1991.04.007

    LIAO Baihan, DAI Zhaohua. Soil buffering capacity to acid precipitation and weathering characteristics of soil minerals[J]. Acta Scientiae Circumstantiae, 1991, 11(4):425-431. doi: 10.13671/j.hjkxxb.1991.04.007

    [30]

    章程, 汪进良, 肖琼, 郭永丽, 苗迎. 桂林潮田河水溶解无机碳昼夜变化与通量[J]. 地球学报, 2021, 42(4):555-564. doi: 10.3975/cagsb.2020.110301

    ZHANG Cheng, WANG Jinliang, XIAO Qiong, GUO Yongli, MIAO Ying. Day and night variations of dissolved inorganic carbon and flux in Chaotian River, Guilin, Guangxi[J]. Acta Geoscientica Sinica, 2021, 42(4):555-564. doi: 10.3975/cagsb.2020.110301

    [31]

    肖佩文, 肖保华. 土柱淋洗实验研究进展及其在土壤有机碳迁移研究中的应用展望[J]. 地球与环境, 2021, 49(1):106-114. doi: 10.14050/j.cnki.1672-9250.2020.48.096

    XIAO Peiwen, XIAO Baohua. Research progress of soil column leaching experiment and its application prospect in soil organic carbon migration[J]. Earth and Environment, 2021, 49(1):106-114. doi: 10.14050/j.cnki.1672-9250.2020.48.096

    [32]

    Mcfee W W, Kelly J M, Beck R H. Acid precipitation effects on soil pH and base saturation of exchange sites[J]. Water, Air and Soil Pollution, 1977, 7(3):401-408.

  • 加载中
    Created with Highcharts 5.0.7访问量Chart context menu近一年内文章摘要浏览量、PDF下载量统计信息摘要浏览量PDF下载量2024-052024-062024-072024-082024-092024-102024-112024-122025-012025-022025-032025-0400.250.50.7511.25Highcharts.com
    Created with Highcharts 5.0.7Chart context menu访问类别分布DOWNLOAD: 1.3 %DOWNLOAD: 1.3 %摘要: 98.7 %摘要: 98.7 %DOWNLOAD摘要Highcharts.com
    Created with Highcharts 5.0.7Chart context menu访问地区分布其他: 20.7 %其他: 20.7 %其他: 1.5 %其他: 1.5 %[]: 1.5 %[]: 1.5 %上海: 2.6 %上海: 2.6 %伊利诺伊州: 0.4 %伊利诺伊州: 0.4 %北京: 2.7 %北京: 2.7 %十堰: 0.2 %十堰: 0.2 %南京: 1.1 %南京: 1.1 %南昌: 0.2 %南昌: 0.2 %哥伦布: 0.2 %哥伦布: 0.2 %嘉兴: 0.2 %嘉兴: 0.2 %天津: 2.4 %天津: 2.4 %宁波: 0.2 %宁波: 0.2 %宣城: 0.4 %宣城: 0.4 %常州: 0.7 %常州: 0.7 %广州: 0.4 %广州: 0.4 %张家口: 1.1 %张家口: 1.1 %扬州: 1.6 %扬州: 1.6 %昆明: 0.7 %昆明: 0.7 %杭州: 0.2 %杭州: 0.2 %格兰特县: 0.4 %格兰特县: 0.4 %桂林: 0.4 %桂林: 0.4 %武汉: 1.3 %武汉: 1.3 %漯河: 4.0 %漯河: 4.0 %石家庄: 0.4 %石家庄: 0.4 %秦皇岛: 0.5 %秦皇岛: 0.5 %芒廷维尤: 29.4 %芒廷维尤: 29.4 %芝加哥: 5.5 %芝加哥: 5.5 %莫斯科: 2.6 %莫斯科: 2.6 %衡阳: 0.2 %衡阳: 0.2 %西宁: 11.0 %西宁: 11.0 %西雅图: 0.5 %西雅图: 0.5 %贵阳: 0.9 %贵阳: 0.9 %赣州: 0.4 %赣州: 0.4 %运城: 1.1 %运城: 1.1 %邯郸: 0.2 %邯郸: 0.2 %郑州: 0.7 %郑州: 0.7 %重庆: 0.5 %重庆: 0.5 %长沙: 0.4 %长沙: 0.4 %长治: 0.5 %长治: 0.5 %黔西南: 0.4 %黔西南: 0.4 %其他其他[]上海伊利诺伊州北京十堰南京南昌哥伦布嘉兴天津宁波宣城常州广州张家口扬州昆明杭州格兰特县桂林武汉漯河石家庄秦皇岛芒廷维尤芝加哥莫斯科衡阳西宁西雅图贵阳赣州运城邯郸郑州重庆长沙长治黔西南Highcharts.com

(5)

(6)

计量
  • 文章访问数:  1635
  • PDF下载数:  147
  • 施引文献:  0
出版历程
收稿日期:  2022-03-20
刊出日期:  2022-10-25

目录

  • 表 1.  供试土壤基本理化性质
    Table 1.  Physico-chemical properties of soil samples
    石灰土厚度/cm
    1020304050607080
    pH 6.62 6.66 6.51 6.63 6.62 6.52 6.54 6.92
    容重/g·cm−3 1.47 1.37 1.24 1.33 1.58 1.68 1.63 1.62
    含水率/% 23.2 26.6 24.7 27.9 29.7 27.0 25.6 27.4
    水溶性Ca2+/×10-6 5.75 11.70 8.97 10.90 8.89 3.10 11.80 8.02
    交换性钙/cmol·kg−1 7.4 4.7 4.2 4.6 6.4 4.8 5.0 5.3
    水溶性Mg2+/×10-9 3.48 3.89 4.77 6.42 5.90 4.69 5.90 6.02
    交换性镁/cmol·kg−1 2.4 1.9 1.7 2.1 2.4 2.5 2.7 2.9
    碳酸钙/×10-3 6.83 8.34 2.77 3.00 3.70 5.55 7.30 5.09
     | Show Table
    DownLoad: CSV
  • 表 2.  石灰土淋出液水化学分析结果 (mg·L−1)
    Table 2.  Hydrochemical analysis results of leachate from lime soil (mg·L−1)
    土层厚度/cmpHCa2+Mg2+K+Na+Cl
    10
    (n=30)
    最大值7.71127.7944.448.066.9060.6333.35547.5913.94
    最小值7.1515.685.722.970.5928.2919.833.991.48
    平均值7.5227.6810.003.691.3050.2624.3161.333.18
    标准差0.1727.079.871.331.329.523.11138.603.18
    变异系数2.3097.7998.6136.18101.7718.9412.79225.98100.30
    20
    (n=30)
    最大值7.4669.9716.6112.515.3644.4634.83249.8414.56
    最小值6.4217.084.406.210.7716.1717.144.861.90
    平均值7.0030.207.848.481.6131.6028.0181.154.36
    标准差0.2815.344.022.270.979.256.6186.152.77
    变异系数4.0050.7951.3226.8360.2029.2623.61106.1563.56
    30
    (n=30)
    最大值7.5355.7912.046.775.0544.4640.76175.9810.50
    最小值6.7319.044.644.100.9228.2920.444.771.90
    平均值7.1833.798.195.421.7939.2830.3477.234.62
    标准差0.2510.542.470.940.864.407.1754.161.86
    变异系数3.5131.1830.1617.3048.0911.2023.6370.1340.40
    40
    (n=30)
    最大值7.5748.209.304.825.2460.6329.72119.707.74
    最小值6.8821.555.293.361.2634.3616.5213.023.11
    平均值7.2329.607.103.992.0142.9821.7359.984.37
    标准差0.147.491.530.540.925.384.5141.551.04
    变异系数1.9525.3221.5413.4345.6812.5320.7769.2823.75
    50
    (n=30)
    最大值7.6050.549.544.376.0364.6730.32115.028.77
    最小值6.8219.554.742.811.1636.3816.775.723.10
    平均值7.2325.586.143.252.0844.5320.6941.804.14
    标准差0.226.911.220.411.086.243.2431.251.43
    变异系数3.1026.9919.8912.7051.9714.0215.6774.7634.54
    60
    (n=30)
    最大值7.5047.8410.433.727.2652.5533.96125.599.41
    最小值6.8519.245.162.281.3832.3416.603.523.60
    平均值7.1624.846.672.742.2342.8420.7839.904.63
    标准差0.186.281.270.361.206.593.8732.111.48
    变异系数2.4525.2919.1012.9753.7815.3918.6480.4832.06
    70
    (n=30)
    最大值7.4543.1811.032.625.8748.5032.28109.339.60
    最小值6.8917.635.561.631.3928.2918.603.303.60
    平均值7.2122.937.111.922.1938.7421.6640.084.92
    标准差0.155.431.270.220.957.603.1029.481.42
    变异系数2.0223.6817.8711.4643.2819.6214.3273.5628.88
    80
    (n=30)
    最大值7.3830.828.251.083.3844.4618.7862.978.53
    最小值6.7417.106.040.751.6028.2914.3316.364.88
    平均值7.0821.317.390.892.0235.3015.4145.095.81
    标准差0.192.890.660.070.384.491.0815.330.87
    变异系数2.6413.548.908.3318.8412.737.0133.9915.03
    注:n为样品数量。
     | Show Table
    DownLoad: CSV
  • 表 3.  不同厚度土壤Ca2+、Mg2+含量及可参与淋溶量
    Table 3.  Contents of Ca2+ and Mg2+ in soil with different thicknesses and the amount that can be involved in leaching
    石灰土厚度/cm
    10 20304050607080
    单层水溶性Ca2+/mg 21.48 39.78 28.03 35.61 34.04 12.84 48.12 32.04
    累积水溶性Ca2+/mg 21.48 61.26 89.29 124.90 158.94 171.78 219.90 251.94
    水溶性Ca2+可参与淋溶量/L 1.07 3.06 4.46 6.24 7.95 8.59 11.00 12.60
    单层水溶性Mg2+/mg 13.00 13.23 14.91 20.97 22.59 19.43 24.06 24.05
    累积水溶性Mg2+/mg 13.00 26.22 41.13 62.10 84.70 104.13 128.19 152.24
    水溶性Mg2+可参与淋溶量/L 2.17 4.37 6.86 10.35 14.12 17.36 21.37 25.37
     | Show Table
    DownLoad: CSV
  • 表 4.  不同厚度土壤累积Ca2+、Mg2+含量及淋失比
    Table 4.  Accumulative content of Ca2+ and Mg2+ in soil with different thicknesses and its leaching-loss ratio
    石灰土厚度/cm
    1020304050607080
    累积土壤水溶性Ca2+/meq 1.07 3.06 4.46 6.24 7.95 8.59 11.00 12.60
    累积土壤水溶性Mg2+/meq 1.08 2.19 3.43 5.18 7.06 8.68 10.68 12.69
    累积土壤水溶性Ca2++Mg2+/meq 2.16 5.25 7.89 11.42 15.00 17.27 21.68 25.28
    淋出液水溶性Ca2++Mg2+/meq 6.95 8.23 8.65 6.81 5.25 5.43 5.47 5.46
    水溶性Ca2++Mg2+淋失比/% 321.76 156.76 109.63 59.63 35.00 31.44 25.23 21.60
    累积土壤交换性Ca2+/mg 11056 17448 22698 28709 38511 46466 54622 63091
    累积土壤交换性Ca2+/meq 553 872 1135 1435 1926 2323 2731 3155
    累积土壤交换性Mg2+/mg 2151 3702 4977 6623 8829 11315 13957 16738
    累积土壤交换性Mg2+/meq 179 308 415 552 736 943 1163 1395
    累积土壤交换性Ca2++Mg2+/meq 732 1180 1550 1987 2662 3266 3894 4550
    淋出液交换性Ca2++Mg2+/meq 4.36 2.80 3.44 3.75 3.88 3.74 3.40 3.11
    交换性Ca2++Mg2+淋失比/% 0.60 0.24 0.22 0.19 0.15 0.11 0.09 0.07
     | Show Table
    DownLoad: CSV
  • 表 5.  不同厚度石灰土可缓冲酸雨(pH=4.5)容量
    Table 5.  Capacity of buffer acid rain(pH=4.5) in different thicknesses of lime
    石灰土厚度/cm
    1020304050607080
    土壤交换性Ca2++Mg2+/meq 732 1 180 1 550 1 987 2 662 3 266 3 894 4 550
    淋出液交换性Ca2++Mg2+/meq 4.36 2.80 3.44 3.75 3.88 3.74 3.40 3.11
    H交换态/mm(pH=4.5) 27 269 68 449 73 184 86 061 111 434 141 836 186 019 237 625
    土壤碳酸钙以钙计/g 10.20 21.55 25.01 28.93 34.60 43.79 55.70 63.84
    溶蚀矿物态钙/g 0.09 0.06 0.07 0.08 0.08 0.07 0.07 0.06
    H矿物态/mm(pH=4.5) 18408 58 336 58 030 58 735 70 247 101 606 129 240 172 815
    H交换态/ H矿物态 1.48 1.17 1.26 1.47 1.59 1.40 1.44 1.38
     | Show Table
    DownLoad: CSV
  • 表 6.  酸雨沉降下不同厚度石灰土产生的碳汇通量
    Table 6.  The carbon sink flux of lime soil with different thicknesses caused by acid rain deposition
    石灰土厚度/cm
    1020304050607080
    碳汇量by/mmol 4.20 2.64 3.28 3.59 3.72 3.58 3.24 2.95
    碳汇量bt/mmol 2.02 1.24 1.56 1.72 1.78 1.71 1.54 1.40
    碳汇量by/碳汇量bt 2.08 2.13 2.10 2.09 2.09 2.09 2.10 2.11
    碳汇量CY/mol·m−2·y−1 0.93 0.59 0.73 0.80 0.83 0.79 0.72 0.65
    碳汇量CT/mol·m−2·y−1 0.45 0.28 0.35 0.38 0.39 0.38 0.34 0.31
     | Show Table
    DownLoad: CSV