Processing math: 100%
  • 全国中文核心期刊
  • 中国科技核心期刊
  • 中国科学引文数据库收录期刊
  • 世界期刊影响力指数(WJCI)报告来源期刊
  • Scopus, CA, DOAJ, EBSCO, JST等数据库收录期刊

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

岩溶试片的岩性差异对估算岩溶速率和碳通量的影响

张小琴 罗维均 王彦伟 蔡先立 吕伊娜 王世杰

张小琴,罗维均,王彦伟,等. 岩溶试片的岩性差异对估算岩溶速率和碳通量的影响[J]. 中国岩溶,2022,41(3):477-487 doi: 10.11932/karst20220314
引用本文: 张小琴,罗维均,王彦伟,等. 岩溶试片的岩性差异对估算岩溶速率和碳通量的影响[J]. 中国岩溶,2022,41(3):477-487 doi: 10.11932/karst20220314
ZHANG Xiaoqin, LUO Weijun, WANG Yanwei, CAI Xianli, LYU Yina, WANG Shijie. Effects of rock tablet lithology difference on estimation of rock dissolution rate and carbon flux[J]. CARSOLOGICA SINICA, 2022, 41(3): 477-487. doi: 10.11932/karst20220314
Citation: ZHANG Xiaoqin, LUO Weijun, WANG Yanwei, CAI Xianli, LYU Yina, WANG Shijie. Effects of rock tablet lithology difference on estimation of rock dissolution rate and carbon flux[J]. CARSOLOGICA SINICA, 2022, 41(3): 477-487. doi: 10.11932/karst20220314

岩溶试片的岩性差异对估算岩溶速率和碳通量的影响

doi: 10.11932/karst20220314
基金项目: 中国科学院战略性先导科技专项(B类)(XDB400202001);国家自然科学基金项目(41673121);贵州省高层次创新型人才培养计划“十”层次人才项目(黔科合平台人才[2016]5648);贵州省省级科技计划项目支助(黔科合基础-ZK[2022]- 564)
详细信息
    作者简介:

    张小琴(1997- ),女,硕士研究生,主要从事喀斯特碳循环研究。E-mail: zhangxiaoqin@mail.gyig.ac.cn

    通讯作者:

    罗维均(1978- ),男,博士,研究员,主要研究方向为喀斯特关键带物质循环。E-mail: luoweijun@vip.gyig.ac.cn

  • 中图分类号: P642.25

Effects of rock tablet lithology difference on estimation of rock dissolution rate and carbon flux

  • 摘要: 研究岩性差异对溶蚀速率的影响有助于提高溶蚀试片法估算岩溶碳汇强度的精确度。本文以贵州省普定县为研究区,将埋放地的主要基岩类型(石灰岩与白云岩)制成标准尺寸的试片,并将其埋设于不同土地利用类型和土壤深度下,经过4个水文年的监测后将估算的结果与前人在同一区域使用标准溶蚀试片的研究结果进行对比分析,结果表明:(1)在相同气候和土壤环境条件下,岩性对溶蚀试片的溶蚀速率有显著影响,且溶蚀速率与岩石中CaO含量呈正相关关系,与MgO含量呈负相关关系;(2)石灰岩与白云岩试片溶蚀速率的差异程度受土地利用及埋放深度的调控,整体上石灰岩溶蚀速率比白云岩溶蚀速率大14%;(3)不同岩性试片估算的岩溶碳汇强度相差较大,标准溶蚀试片估算的结果比埋放地基岩试片估算的结果高。故使用溶蚀试片法估算区域岩溶碳通量时应考虑埋放地基岩的岩性,或者对基于标准溶蚀试片的估算结果进行校正,才能准确反映区域尺度真实的岩溶碳通量大小。

     

  • 自工业革命以来,化石燃料燃烧及土地利用变化等人类活动使得大气中CO2浓度剧增,从1750年的277 ×10−6上升到2021年的416.87×10−6[1]。人类活动的碳排放40%直接增加了大气中的CO2浓度, 30%被海洋所吸收,还有10%左右的剩余部分被称为遗失碳汇,可能分布在陆地生态系统中[2-3]。而陆地生态系统组成十分复杂,确定具体的碳汇分布难度较大,特别是在喀斯特地区,地表的二维三元结构使得碳循环路径更多、过程更为复杂,碳循环的研究需考虑水-岩-土-气-生5个子系统[4]。近来的很多研究表明,碳酸盐岩的化学风化过程将在短期和长期产生碳汇[5-6]。20世纪90年代,科学家利用石灰岩溶蚀试片法、水化学法和扩散边界层理论估算出全球范围内岩溶作用每年的碳回收量为0.22~0.608 Pg C,约占未知碳汇的1/3[7]。随着研究的加深,发现碳酸盐岩溶解的无机碳可在水循环的参与下与水生生物光合作用产生联系[7-9],并且得出该过程的碳汇强度为0.5 Pg C·a−1,约占陆地遗漏碳汇的50%[10]。这些研究为全球碳循环研究提供了理论补充,表明对岩溶碳通量的研究不容忽视[11-12]

    在定量评估岩溶作用速率及其碳汇强度的研究中,标准溶蚀试片法是最简单且常用的。在中国,标准溶蚀试片为采自桂林泥盆系融县组纯的石灰岩,经加工制成标准尺寸的岩石试片(圆饼形,直径40 mm,厚3 mm),将其放置于空气中、土壤表面或地下不同深度(−5 cm、−20 cm、−50 cm),一定时间内取回称重计算其溶蚀量[13],并采用以点带面的方式,直接计算区域内的溶蚀速率,从而简化了因碳酸盐岩种类太多带来的复杂计算 [14]。自20世纪90年代,学者通过埋放标准溶蚀试片进行溶蚀实验,研究土壤水、CO2、有机质、微生物等环境因素对碳酸盐岩溶蚀速率的影响[15-19],并估算区域内岩溶碳汇强度[20-23]。但是,由于喀斯特地区地质背景具有高度的异质性,如果均使用标准溶蚀试片(取自桂林泥盆系的石灰岩)来研究岩溶碳汇效应,会产生较大的误差[24-25]。基于此,本文以中国贵州省普定县为研究区,选择当地基岩(石灰岩与白云岩)制成与标准试片形状和大小相同的试片,并参考标准试片的野外试验方法进行埋设,经过四个水文年的试片收集和称重,对该区不同土地利用类型的岩溶作用速率及岩溶碳汇强度进行估算,通过与标准溶蚀试片的相关研究进行对比分析,并探讨差异形成原因,为准确估算喀斯特地区的碳汇量提供科学依据。

    普定县位于贵州省中偏西部(105°27'49"E~105°58'51"E,26° 26' 36" N~26° 31' 42" N),海拔1 100~1 600 m。属于亚热带季风湿润气候。试验期间年均气温16.1 ℃,日最高温度29 ℃,最低温度−0.4 ℃,年平均降水量在1 079.3到1 457.7 mm。普定县为高原丘陵、山地地形,县域地势为中间低南北高,主要为峰丛、溶洞、暗河、洼地等喀斯特地貌。县域碳酸盐岩地层出露面积约占总面积的84.27%,主要为石炭系、二叠系、三叠系,其中以三叠系出露面积最大,约占60.7%。全县土壤种类较多,主要为石灰土、黄壤、水稻土、山地黄棕壤、山地灌丛草甸土等,其中,石灰土占比最大,约为63.7%。普定县境内代表性植被为常绿栎林、常绿落叶混交林和马尾松林,灌丛和草地主要为藤刺灌丛、火棘、白茅、青蒿等植物[26]

    研究区位于普定县马官镇北部的天龙山(26°14'40″N,105°45'45″E)和白岩镇东部的讲义村(26°15′57″N,105°50′26″E)(图1)。天龙山区域为喀斯特峰丛洼地地貌,基岩为三叠系关岭组中段(T2g2),岩性主要是薄层至中厚灰岩、泥灰岩[27],土壤类型为黑色石灰土,土层较薄,地面平均岩石裸露率为44.7%,土被不连续[28],山体陡峭,森林植被保存较好,主要以圆果化香树、云南鼠刺、窄叶石栎为优势种的常绿落叶阔叶混交林[29];农田以水田和旱地为主,主要种植水稻、玉米和油菜。讲义村区域的基岩为三叠系关岭组上段(T2g3),岩性主要为紫红色微晶白云岩[27],土壤类型为黑色石灰土,土被不连续,土层厚度在10~50 cm。次生林地以椤木石楠为优势种的常绿落叶混交林,山体草甸主要的植物为白茅,水田交替种植水稻和油菜,旱地种植玉米和蔬菜。

    图  1  研究区及试片埋放位置图[27]
    Figure  1.  Map of study area and buried sites of rock tablet

    将研究区的基岩加工制成标准尺寸的溶蚀试片,选取典型的土地利用类型(次生林地、灌丛、草地、旱地、水田),每个土地利用类型设置一个埋放点,埋放深度为地下5 cm、20 cm、50 cm,每层埋放4组试片,每组为一片白云岩(讲义村基岩)和一片石灰岩(天龙山基岩)。试片于2016年12月进行埋放,后期每年取出一组测定其溶蚀量。埋放点概况如表1,土壤剖面如图2图3

    表  1  不同土地利用类型埋放点概况
    Table  1.  Overview of buried sites of different land use types
    地点土地利用类型土壤剖面基本情况位置及其他
    天龙山 次生林 A层为25 cm黑色石灰土,碎石比30%,B层为黄土未见底 半山腰处
    灌丛 碎石比4%,草根发达, A层为18 cm黑色石灰土,B层为黄色土壤,不见底 山脚处,土壤疏松度低
    水田 A层38 cm黑色石灰土,B层为黄色土,未见底 洼地,南部为峰丛,水稻交替种植,夏季长时间处于淹水状态
    旱地 A层为40 cm为黑色石灰土,B层为黄色土未见底 洼地,四周为峰丛,种植玉米及蔬菜,人为干扰较多
    讲义村 次生林 A层35 cm,为黑色石灰土,根系发达,B层未见底 山腰处,植被丰富,枯枝落叶多
    灌丛 碎石比20%; A层21 cm,黑色石灰土,B层黄色土壤,未见底 山脚处,土壤疏松度低
    草地 碎石比50%,A层20~30 cm黑色石灰土,根系发达,B层厚度10~15 cm,以下为基岩 山坡草地上,土壤疏松度低
    旱地 A层5~20 cm,为黑色石灰土;B层为黄色土壤,未见底 山下洼地处,种植玉米及蔬菜,人为干扰较多
    水田 A层43 cm,为黑色石灰土,以下为黄色土壤 山下洼地处,油菜地、水稻交替种植,夏季长时间处于淹水状态
    下载: 导出CSV 
    | 显示表格
    图  2  天龙山不同土地利用类型土壤剖面
    Figure  2.  Soil profiles of different land use types in Tianlong mountain
    图  3  讲义村不同土地利用类型土壤剖面
    Figure  3.  Soil profiles of different land use types in Jiangyi

    每片试片埋放前用纯净水洗净并烘干,称重两次取平均值,记录下初始重量 W1;埋放后取出时再用纯净水洗净,放入105 ℃的烘箱烘干24 h,再次称重两次取平均值W2。称重精度为0.01 mg。计算试片日均溶蚀速率[13, 30]

    ER=(W1W2)×10000/(T×S)
    (1)

    式中:ER为日均单位面积溶蚀量,即溶蚀速率 (mg·m−2·d−1);W1为试片初始重量(mg);W2为试片埋放后重量(mg);T 为埋放天数(d);S为试片表面积(约28. 91 cm2)。采用以下公式计算两地不同土地利用下的岩溶碳汇强度:

    F=ER×365×103×R×M1/M2
    (2)
    F=ER×365×103×R×2M1/M2
    (3)

    公式(2)为天龙山(石灰岩区)岩溶碳汇强度计算公式,公式(3)为讲义(白云岩区)岩溶碳汇强度计算公式。式中,F为碳汇强度,单位是tCO2·km−2·a−1ER为岩石试片的日单位面积溶蚀量,单位是mg·m−2·d−1R为岩石试片的碳酸盐岩纯度;M1为CO2的分子量44;M2为石灰岩试片岩石成分CaCO3的分子量100或白云岩试片岩石成分CaMg (CO3)2的分子量184。

    将天龙山和讲义村两埋放地4年内回收的全部试片计算日均溶蚀速率,并取平均值得到表2

    表  2  不同土地利用下不同深度的试片日均溶蚀速率
    Table  2.  Daily average dissolution rate of rock tablets at different soil depths under different land uses
    地点试片岩性埋放深度/cm次生林灌丛草地旱地水田
    天龙山石灰岩日均溶蚀速率
    mg·m−2·d−1
    地下5119.94138.69106.84186.99
    地下2018.91171.64159.80
    地下508.73131.3269.75
    平均值49.19136.60138.85
    白云岩日均溶蚀速率
    mg·m−2·d−1
    地下555.09101.9179.12144.75
    地下2017.0556.46143.21125.42
    地下503.9542.12110.0537.92
    平均值25.3666.83110.79102.70
    讲义村石灰岩日均溶蚀速率
    mg·m−2·d−1
    地下535.2528.4469.2863.08110.13
    地下2035.9321.7946.11114.19121.32
    地下5044.0420.3130.6385.90136.47
    平均值38.4123.5148.6887.72122.64
    白云岩日均溶蚀速率
    mg·m−2·d−1
    地下515.8012.6537.0943.5478.91
    地下2014.6311.3718.4267.7187.43
    地下5019.656.957.8356.15110.54
    平均值16.6910.3321.1155.8092.29
    注:溶蚀速率为4年内试片日均溶蚀速率的平均值;−为试片缺失。
    下载: 导出CSV 
    | 显示表格

    表2可以得出不同土地利用类型不同层位试片的溶蚀速率特征如下:

    (1)不同土地利用类型地下0~50 cm范围内,天龙山石灰岩试片的平均溶蚀速率为水田(138.85 mg·m−2·d−1)>旱地(136.60 mg·m−2·d−1)>次生林(49.19 mg·m−2·d−1),白云岩试片的溶蚀速率为旱地(110.79 mg·m−2·d−1)>水田(102.70 mg·m−2·d−1)>灌丛(66.83 mg·m−2·d−1)>次生林(25.36 mg·m−2·d−1)。总体来看,天龙山两种岩性试片的溶蚀速率都表现为水田和旱地的溶蚀速率最大,次生林的溶蚀速率最小。讲义村两种岩性试片的溶蚀速率表现为水田>旱地>草地>次生林>灌丛,但不同土地利用类型的溶蚀速率差异较明显。

    (2)相同土地利用类型、不同层位二种试片的溶蚀速率差异较大,规律不一。具体表现为:在天龙山,不同埋放深度下,埋放在次生林和水田的石灰岩和白云岩试片都表现为地下5 cm处溶蚀速率最大,地下20 cm次之,地下50 cm处最小;埋放在旱地的石灰岩和白云岩试片则表现为地下20 cm处溶蚀速率最大,地下50 cm处次之,最小为地下5 cm处。其中,石灰岩和白云岩试片溶蚀速率最大值都出现在水田地下5 cm处,分别为186.99 mg·m−2·d−1和144.75 mg·m−2·d−1,最小值出现在次生林地下50 cm处,分别为8.73 mg·m−2·d−1和3.95 mg·m−2·d−1。在讲义村各埋试点,埋放在灌丛和草地的石灰岩和白云岩试片的溶蚀速率都呈现出随埋深增加逐渐减小的趋势;而埋放在水田中的两种试片的溶蚀速率都表现出随埋深增加逐渐增大;旱地里两种试片的溶蚀速率在地下20 cm处最大,地下50 cm处次之,地下5 cm处最小;埋放在次生林中石灰岩和白云岩试片的溶蚀速率与埋放深度的关系不大。在讲义村,两种试片的溶蚀速率最大出现在水田地下50 cm处,最小出现在灌丛的地下50 cm处。

    (3)天龙山和讲义村各土地利用类型及各埋放深度,都表现出石灰岩试片的溶蚀速率比白云岩试片的溶蚀速率高。天龙山次生林地下5 cm处石灰岩试片比白云岩试片溶蚀速率大65 mg·m−2·d−1左右,埋深在地下20 cm和50 cm处两种试片的溶蚀速率差异不大;在灌丛中,地下5 cm处石灰岩试片比白云岩试片溶蚀速率大36.78 mg·m−2·d−1;在水田和旱地中,两种试片的溶蚀速率差异较稳定,在25~45 mg·m−2·d−1之间。在讲义村,埋在旱地土壤下20 cm处石灰岩试片溶蚀速率比白云岩试片溶蚀速率大47 mg·m−2·d−1左右,除此之外,其他土地利用类型石灰岩和白云岩试片的溶蚀速率差异都稳定在10 ~30 mg·m−2·d−1之间。

    3.2.1   石灰岩与白云岩溶蚀速率关系

    对研究区4年内所有的白云岩与石灰岩试片的日均溶蚀速率进行相关性分析(图4)发现,绝大部分的点分布在y=x线上方,K>1(K=石灰岩溶蚀速率/白云岩溶蚀速率),即石灰岩试片的溶蚀速率整体大于白云岩试片的溶蚀速率,统计分析结果得出,石灰岩溶蚀速率比白云岩溶蚀速率大14%左右,相关性系数为0.85。

    图  4  石灰岩试片与白云岩试片溶蚀速率关系
    Figure  4.  Relationship between dissolution rates of limestone and dolomite rock tablets
    3.2.2   土地利用类型与溶蚀速率的关系

    将4年内的试片日均溶蚀速率的平均值(表2)按土地利用类型进行划分,对3个土壤深度下石灰岩与白云岩试片的溶蚀速率进行相关性分析(图5),在石灰岩溶蚀速率大于白云岩溶蚀速率的基础上(K>1),同一土地利用类型的3个埋放深度的K值不一致,即3个埋放深度的石灰岩与白云岩溶蚀速率差异程度不同,3个点趋势线的斜率越接近1,说明该土地利用类型下3个埋放深度对石灰岩与白云岩溶蚀速率差异的影响越小。天龙山的次生林、旱地和水田3个埋放点趋势线的斜率分别为2.28、1.01和1.08,讲义村的旱地、次生林、草地、灌丛、水田3个埋放点趋势线的斜率分别为2.11、1.78、1.31、1.14、0.80。

    图  5  不同土地利用类型下石灰岩与白云岩试片溶蚀速率关系
    Figure  5.  Relationship between dissolution rate of limestone and dolomite rock tablets under different land use types
    3.2.3   土壤埋放深度与溶蚀速率的关系

    同理,将4年内的试片日均溶蚀速率的平均值(表2)按埋放深度进行分类,对研究区全部土地利用类型的石灰岩与白云岩试片的溶蚀速率进行相关性分析(图6),在石灰岩溶蚀速率大于白云岩溶蚀速率的基础上(K>1),同一埋放深度下,不同土地利用类型的K值不一致,埋放点趋势线斜率越接近1,说明研究区内不同土地利用类型对石灰岩与白云岩溶蚀速率的影响越小。天龙山的土壤下5 cm、20 cm、50 cm处埋放点趋势线的斜率分别为0.84、1.24、1.11,讲义村的土壤下5 cm、20 cm、50 cm处埋放点趋势线的斜率分别为1.20、1.30、1.08。

    图  6  不同土壤深度条件下石灰岩与白云岩试片溶蚀速率关系
    Figure  6.  Relationship between dissolution rate of limestone and dolomite rock tablets at different soil depths
    3.3.1   埋放地的岩溶碳汇强度

    本研究在天龙山和讲义村均使用了石灰岩和白云岩两种岩性试片,天龙山基岩为石灰岩,讲义村基岩为白云岩,根据天龙山各土地利用类型下的石灰岩试片的平均溶蚀速率及其主要化学成分(表3),由公式(2)计算出天龙山不同土地利用类型的岩溶碳汇强度。同理,由公式(3)计算出讲义村不同土地利用类型的岩溶碳汇强度(表4)。

    表  3  不同岩石试片主要化学成分
    Table  3.  Main chemical components of different rock tablets
    岩性组分/%
    SiO2CaOMgOK2ONa2OCO2
    石灰岩6.5647.102.260.4570.03441.58
    白云岩1.2630.1421.40<0.010.02346.41
    标准溶蚀试片[31]55.610.16
    下载: 导出CSV 
    | 显示表格
    表  4  不同土地利用类型地下平均岩溶碳汇强度
    Table  4.  Average karst carbon sink intensity under different land use types
    地点项目次生林灌丛草地旱地水田
    天龙山 试片溶蚀速率/mg·m−2·d−1 49.19 136.60 138.85
    岩溶碳汇强度/tCO2·km−2·a−1 7.47 20.73 21.07
    讲义村 试片溶蚀速率/mg·m−2·d−1 16.69 10.33 21.11 55.80 92.29
    岩溶碳汇强度/tCO2·km−2·a−1 2.84 1.75 3.40 9.48 15.68
    下载: 导出CSV 
    | 显示表格
    3.3.2   不同岩性试片碳汇估算对比

    本研究中,同一埋放点埋放的两种岩性试片其溶蚀速率差异显著,根据公式(2)和(3)计算得出天龙山各土地利用类型地下50 cm处两种试片对应的岩溶碳汇强度(表5),次生林、旱地和水田下石灰岩试片的岩溶碳汇强度分别为1.33 tCO2·km−2·a−1、19.93 tCO2·km−2·a−1、10.59 tCO2·km−2·a−1,均大于白云岩试片估算的结果。曾成等[17]在普定县陈旗村使用标准溶蚀试片(采自广西桂林七星岩上泥盆统融县组)的研究结果表明:旱地和水田的地下50 cm处的岩溶碳汇强度分别为22.49 tCO2·km−2·d−1和13.02 tCO2·km−2·a−1,均高于本研究的石灰岩试片和白云岩试片的估算结果,灌丛的岩溶碳汇为2.60 tCO2·km−2·a−1,比本研究白云岩试片估算的结果低。可见在环境条件相同时,试片岩性的差异会造成岩溶碳汇强度估算的差异,总体上标准试片法估算的结果较大。

    表  5  地下50 cm处不同岩性试片估算的岩溶碳汇强度
    Table  5.  Karst carbon sink intensity at 50 cm underground estimated by different lithology rock tablets
    岩性不同土地利用类型岩溶碳汇强度/ tCO2·km−2·a−1
    次生林灌丛旱地水田
    石灰岩1.3319.9310.59
    白云岩0.687.1518.696.44
    标准溶蚀试片2.6022.4913.02
    下载: 导出CSV 
    | 显示表格

    岩溶作用是CO2溶于水形成弱酸溶液溶解可溶性岩石的化学溶蚀过程,因此试片的溶蚀作用强度直接受土壤环境中水份含量和CO2浓度影响。此外,有研究认为试片的溶蚀作用还受到土壤有机质含量、植被状况以及土壤微生物含量等因素的影响。土壤有机质含量越高时,土壤pH越低,越容易使试片溶蚀[17];植物根系分泌的有机酸可降低土壤pH以及提高土壤碳酸酐酶活性,从而促进岩石化学风化[33];土壤微生物的活性和多样性也可加速碳酸钙类岩石的溶蚀[32]。另外,土壤中的碳酸盐岩角砾会与降雨后下渗的土壤水发生溶蚀反应,从而削弱了继续下渗的土壤水的溶蚀能力,使深部的溶蚀试片的溶蚀量减少。本研究中,次生林、灌丛、水田、旱地和草地的土壤水分含量、土壤孔隙度、植被状况、凋落物含量以及各埋放点土壤的碎石比均存在较大差异(表1),故同一岩性的试片在不同土地利用类型下溶蚀速率差异显著。

    在相同环境条件下,不同岩性岩石的矿物组成和结构决定了自身的溶蚀速率[33]。研究区内的基岩主要由方解石(CaCO3)和白云石(CaMg(CO3)2)两种矿物组成,两种矿物晶体结构不同,因此其溶蚀作用也存在一定差异。石灰岩的主要矿物成分为方解石,白云岩则为白云石。在溶蚀过程中,石灰岩沿分散孤立的晶间空隙溶蚀,形成很多溶蚀坑,白云岩则沿节理裂隙进行溶蚀,使裂隙不断扩大而崩解[34]。一般情况下,碳酸盐岩类型的溶蚀比为:石灰岩>云灰岩>泥质云灰岩>大理岩>泥质灰岩>灰云岩>泥质灰云岩>白云岩>泥质白云岩[35]。本研究中,石灰岩和白云岩试片的主要化学成分如表3,结合曾成等[22]在普定县使用标准溶蚀试片的试验结果,将试片的化学成分与其溶蚀速率进行相关性分析,得出在水田和旱地两种土地利用类型下,岩石溶蚀速率与CaO含量呈正相关,与MgO含量呈现负相关,相关性系数均大于0.90(图7),且水田和旱地的地下试片溶蚀速率受MgO、CaO含量影响的程度不同,这可能与两种土地利用类型的土壤环境有关。

    图  7  溶蚀速率与岩石试片化学成分含量的关系
    Figure  7.  Relation between dissolution rate and chemical composition content of rock tablets

    不同土地利用类型及不同土壤深度下,石灰岩与白云岩溶蚀速率的比值不同,其原因除与试片岩性有关外,还与试片埋放点的土壤环境有关,是各因素的共同作用的最终结果。同一土地利用类型,在不同深度下,试片的溶蚀速率主要受不同深度土壤环境条件的影响,从图5分析得出,在天龙山,地下埋放深度对埋放在水田和旱地的石灰岩与白云岩试片溶蚀速率差异程度的影响较小,对埋放在次生林下的影响相对较大;而在讲义村,埋放深度对草地下石灰岩与白云岩试片溶蚀速率差异程度的影响较小,而对旱地和次生林下的试片影响相对较大。同理,在同一区域的相同土壤深度下,受不同土地利用类型的土壤环境条件的影响,石灰岩与白云岩的溶蚀速率差异程度相差较大,从图6分析得出,研究区内地下5 cm和20 cm处土地利用类型对石灰岩与白云岩溶蚀速率差异程度的影响较大,50 cm处的最小。总之,影响试片溶蚀速率的主要因素是水和CO2的含量及其与岩石的接触面积,当反应条件已达到充分时,溶蚀速率就主要由岩性来决定。

    使用溶蚀试片法估算的岩溶碳汇强度是在已知试片溶蚀速率的基础上进行的,岩溶碳汇强度与溶蚀速率两者有一定的正相关性,但相关程度受碳酸盐岩纯度的影响。根据表5得出在水田和旱地下,标准溶蚀试片法估算的岩溶碳汇强度比本研究采用的石灰岩与白云岩试片估算的岩溶碳汇强度都大,本研究石灰岩试片的岩溶碳汇强度又比白云岩的岩溶碳汇强度大,结合三种试片的化学成分(表3)认为,岩性差异是导致相同土地利用类型下岩溶碳汇强度差异的主要原因。这和黄奇波等[22]在山西吕梁山岩溶流域的研究结果一致,使用当地岩性的试片估算的岩溶碳汇强度(1.738 tCO2·km−2·a−1)明显低于标准试片的计算结果(1.821 tCO2·km−2·a−1)。而与水化学法相比,多项研究表明,采用溶蚀试片法估算得出的岩溶碳汇强度比采用水化学法计算的结果要低,仅为水化学法计算结果的1/5到1/6[22, 36],这除了与流域内土壤中存在的原生和次生碳酸盐有关外,可能还与当地基岩复杂的岩性有关,标准溶蚀试片不能完全代表流域内的基岩,从而造成标准试片法计算的结果偏低。

    综上,岩性及土壤环境条件对试片的溶蚀速率影响较大,当前的研究虽未能将每个因素的影响进行量化,但得出了岩石的溶蚀速率与自身CaO含量呈一定的正相关关系,虽然在不同环境条件下正相关程度有一定的差异,但总的看来,当基岩CaO 含量相差17%时,其溶蚀速率相差14%,结合基岩的碳酸盐岩纯度可计算岩溶碳汇强度。

    (1)使用埋放地的基岩试片估算了研究区不同土地利用类型的溶蚀速率及岩溶碳汇强度,并与前人使用标准溶蚀试片在该区的研究进行对比分析,发现不同土地利用类型、不同埋放深度等条件下,石灰岩试片比白云岩试片岩溶速率大14%,产生原因主要是受溶蚀试片CaO和MgO的相对含量控制;

    (2)在不同研究区,由于岩性差异,采用标准溶蚀试片估算岩溶碳汇强度会造成结果的高估或低估,而基于埋放地基岩岩性的试片更能准确地代表当地的岩溶速率;

    (3)溶蚀试片法与水化学法计算的岩溶碳汇结果往往存在差异,可能与采用的标准试片法岩性与当地基岩岩性不同有关,需要进一步研究确认。

  • 图  1  研究区及试片埋放位置图[27]

    Figure  1.  Map of study area and buried sites of rock tablet

    图  2  天龙山不同土地利用类型土壤剖面

    Figure  2.  Soil profiles of different land use types in Tianlong mountain

    图  3  讲义村不同土地利用类型土壤剖面

    Figure  3.  Soil profiles of different land use types in Jiangyi

    图  4  石灰岩试片与白云岩试片溶蚀速率关系

    Figure  4.  Relationship between dissolution rates of limestone and dolomite rock tablets

    图  5  不同土地利用类型下石灰岩与白云岩试片溶蚀速率关系

    Figure  5.  Relationship between dissolution rate of limestone and dolomite rock tablets under different land use types

    图  6  不同土壤深度条件下石灰岩与白云岩试片溶蚀速率关系

    Figure  6.  Relationship between dissolution rate of limestone and dolomite rock tablets at different soil depths

    图  7  溶蚀速率与岩石试片化学成分含量的关系

    Figure  7.  Relation between dissolution rate and chemical composition content of rock tablets

    表  1  不同土地利用类型埋放点概况

    Table  1.   Overview of buried sites of different land use types

    地点土地利用类型土壤剖面基本情况位置及其他
    天龙山 次生林 A层为25 cm黑色石灰土,碎石比30%,B层为黄土未见底 半山腰处
    灌丛 碎石比4%,草根发达, A层为18 cm黑色石灰土,B层为黄色土壤,不见底 山脚处,土壤疏松度低
    水田 A层38 cm黑色石灰土,B层为黄色土,未见底 洼地,南部为峰丛,水稻交替种植,夏季长时间处于淹水状态
    旱地 A层为40 cm为黑色石灰土,B层为黄色土未见底 洼地,四周为峰丛,种植玉米及蔬菜,人为干扰较多
    讲义村 次生林 A层35 cm,为黑色石灰土,根系发达,B层未见底 山腰处,植被丰富,枯枝落叶多
    灌丛 碎石比20%; A层21 cm,黑色石灰土,B层黄色土壤,未见底 山脚处,土壤疏松度低
    草地 碎石比50%,A层20~30 cm黑色石灰土,根系发达,B层厚度10~15 cm,以下为基岩 山坡草地上,土壤疏松度低
    旱地 A层5~20 cm,为黑色石灰土;B层为黄色土壤,未见底 山下洼地处,种植玉米及蔬菜,人为干扰较多
    水田 A层43 cm,为黑色石灰土,以下为黄色土壤 山下洼地处,油菜地、水稻交替种植,夏季长时间处于淹水状态
    下载: 导出CSV

    表  2  不同土地利用下不同深度的试片日均溶蚀速率

    Table  2.   Daily average dissolution rate of rock tablets at different soil depths under different land uses

    地点试片岩性埋放深度/cm次生林灌丛草地旱地水田
    天龙山石灰岩日均溶蚀速率
    mg·m−2·d−1
    地下5119.94138.69106.84186.99
    地下2018.91171.64159.80
    地下508.73131.3269.75
    平均值49.19136.60138.85
    白云岩日均溶蚀速率
    mg·m−2·d−1
    地下555.09101.9179.12144.75
    地下2017.0556.46143.21125.42
    地下503.9542.12110.0537.92
    平均值25.3666.83110.79102.70
    讲义村石灰岩日均溶蚀速率
    mg·m−2·d−1
    地下535.2528.4469.2863.08110.13
    地下2035.9321.7946.11114.19121.32
    地下5044.0420.3130.6385.90136.47
    平均值38.4123.5148.6887.72122.64
    白云岩日均溶蚀速率
    mg·m−2·d−1
    地下515.8012.6537.0943.5478.91
    地下2014.6311.3718.4267.7187.43
    地下5019.656.957.8356.15110.54
    平均值16.6910.3321.1155.8092.29
    注:溶蚀速率为4年内试片日均溶蚀速率的平均值;−为试片缺失。
    下载: 导出CSV

    表  3  不同岩石试片主要化学成分

    Table  3.   Main chemical components of different rock tablets

    岩性组分/%
    SiO2CaOMgOK2ONa2OCO2
    石灰岩6.5647.102.260.4570.03441.58
    白云岩1.2630.1421.40<0.010.02346.41
    标准溶蚀试片[31]55.610.16
    下载: 导出CSV

    表  4  不同土地利用类型地下平均岩溶碳汇强度

    Table  4.   Average karst carbon sink intensity under different land use types

    地点项目次生林灌丛草地旱地水田
    天龙山 试片溶蚀速率/mg·m−2·d−1 49.19 136.60 138.85
    岩溶碳汇强度/tCO2·km−2·a−1 7.47 20.73 21.07
    讲义村 试片溶蚀速率/mg·m−2·d−1 16.69 10.33 21.11 55.80 92.29
    岩溶碳汇强度/tCO2·km−2·a−1 2.84 1.75 3.40 9.48 15.68
    下载: 导出CSV

    表  5  地下50 cm处不同岩性试片估算的岩溶碳汇强度

    Table  5.   Karst carbon sink intensity at 50 cm underground estimated by different lithology rock tablets

    岩性不同土地利用类型岩溶碳汇强度/ tCO2·km−2·a−1
    次生林灌丛旱地水田
    石灰岩1.3319.9310.59
    白云岩0.687.1518.696.44
    标准溶蚀试片2.6022.4913.02
    下载: 导出CSV
  • [1] Fortunat Joos, Renato Spahni. Rates of change in natural and anthropogenic radiative forcing over the past 20000 years[J]. Proceedings of the National Academy of Sciences of the United States of America, 2008, 105(5):1425-1430.
    [2] Pierre Friedlingstein, Michael O'Sullivan, Matthew W. Jones, Robbie M. Andrew, Judith Hauck, Are Olsen, et al. Global Carbon Budget 2020[J]. Earth System Science Data, 2020, 12(4):3269-3340.
    [3] PIAO Shilong, HE Yue, WANG Xuhui, CHEN Fahu. Estimation of China's terrestrial ecosystem carbon sink: Methods, progress and prospects[J]. Science China Earth Sciences, 2022, 65(4):641-651.
    [4] 王世杰, 刘再华, 倪健, 闫俊华, 刘秀明. 中国南方喀斯特地区碳循环研究进展[J]. 地球与环境, 2017, 45(01): 2-9.

    WANG Shijie, LIU Zaihua, NI Jian, YAN Junhua, LIU Xiuming. A review of research progress and future prospective of carbon cycle in karst area of south China. Erath and environment. 2017, 45(1): 2-9.
    [5] LIU Zaihua, Wolfgang DREYBRODT, WANG Haijing. A new direction in effective accounting for the atmospheric CO2 budget: Considering the combined action of carbonate dissolution, the global water cycle and photosynthetic uptake of DIC by aquatic organisms[J]. Earth Science Reviews, 2010, 99(3-4):162-172.
    [6] ZHOU Guoqing, HUANG Jingjin, TAO Xiaodong, LUO Qingli, ZHANG Rongting, LIU Zaihua. Overview of 30 years of research on solubility trapping in Chinese karst[J]. Earth-Science Reviews, 2015, 146:183-194.
    [7] YUAN Daoxian. The carbon cycle in karst[J]. Zeitschrift fuer Geomorphologie, 1997, 108:91-102.
    [8] LIU Zaihua, ZHAO Jinbo. Contribution of carbonate rock weathering to the atmospheric CO2 sink[J]. Environmental Geology, 2000, 39(9):1053-1058.
    [9] LIU Zaihua, Wolfgang Dreybrodt. Significance of the carbon sink produced by H2O-carbonate-CO2-aquatic phototroph interaction on land[J]. Science Bulletin, 2015, 60(2):182-191.
    [10] LIU Zaihua, G. L. Macpherson, Chris Groves, Jonathan B. Martin, YUAN Daoxian, ZENG Sibo. Large and active CO2 uptake by coupled carbonate weathering[J]. Earth-Science Reviews, 2018, 182:42-49.
    [11] 袁道先. 碳循环与全球岩溶[J]. 第四纪研究, 1993, 13(1):1-6. doi: 10.3321/j.issn:1001-7410.1993.01.001

    YUAN Daoxian. Carbon cyclce and global karst[J]. Quaternary Sciences, 1993, 13(1):1-6. doi: 10.3321/j.issn:1001-7410.1993.01.001
    [12] CAO Jianhua, YUAN Daoxian, Chris GROVES, HUANG Fen, YANG Hui, LU Qian. Carbon Fluxes and Sinks: the Consumption of Atmospheric and Soil CO2 by Carbonate Rock Dissolution[J]. Acta Geologica Sinica, 2012, 86(4):963-972.
    [13] 袁道先, 蔡桂鸿. 岩溶环境学[M]. 重庆: 重庆出版社, 1988: 62-64.

    YUAN Daoxian, CAI Guihong. The science of karst environment[M]. Chongqing: Chongqing Publication House, 1988: 62-64.
    [14] Kristina Krklec, David Domínguez-Villar, Dražen Perica. Use of rock tablet method to measure rock weathering and landscape denudation[J]. Earth-Science Reviews, 2021, 212(4):103449.
    [15] ZHANG Cheng, PEI Jianguo, XIE Yunqiu, CAO Jianghua, WANG Lanling. Impact of land use covers upon karst processes in a typical Fengcong depression system of Nongla, Guangxi, China[J]. Environmental Geology, 2008, 55(8):1621-1626.
    [16] 谢芳, 傅瓦利, 王晓阳, 蒲 鹏, 张蕾, 谭波, 彭景涛, 甄晓君. 重庆中梁山碳酸盐岩溶蚀速率对季节的响应研究[J]. 中国岩溶, 2010, 29(4):410-413. doi: 10.3969/j.issn.1001-4810.2010.04.010

    XIE Fang, FU Wali, WANG Xiaoyang, PU Peng, ZHANG Lei, TAN Bo, PENG Jingtao, ZHEN Xiaojun. Study on response of carbonate rock dissolution rate to the season in Zhongliang Mountain, Chongqing[J]. Carsologica Sinica, 2010, 29(4):410-413. doi: 10.3969/j.issn.1001-4810.2010.04.010
    [17] 刘文, 张强, 贾亚男. 夏季不同土地利用方式下的溶蚀作用研究: 以重庆青木关岩溶槽谷区为例[J]. 中国岩溶, 2012, 31(1):1-6. doi: 10.3969/j.issn.1001-4810.2012.01.001

    LIU Wen, ZHANG Qiang, JIA Yanan. Karstification under different lan-use patterns in sumer: A case study in the Qingmuguan karst valley, Chongqing[J]. Carsologica Sinica, 2012, 31(1):1-6. doi: 10.3969/j.issn.1001-4810.2012.01.001
    [18] 蓝家程, 傅瓦利, 彭景涛, 周小萍, 肖时珍, 袁波. 不同土地利用方式土下岩溶溶蚀速率及影响因素[J]. 生态学报, 2013, 33(10):3205-3212.

    LAN Jiacheng, FU Wali, PENG Jingtao, ZHOU Xiaoping, XIAO Shizhen, YUAN Bo. Dissolution rate under soil in karst areas and the influencing factors of different land use patterns[J]. Acta Ecologica Sinica, 2013, 33(10):3205-3212.
    [19] SHAO Mingyu, ZHANG Liankai, LIU Pengyu, CAO Jianhua, QIN Xiaoqun, HUANG Qibo, LUO Mingming, ZHANG Chunlai. Influential factors and spatial suitability of the method of limestone tablets in karst carbon cycle study in China[J]. Carbonates and Evaporites, 2020, 35(3):85.
    [20] 徐胜友, 蒋忠诚. 我国岩溶作用与大气温室气体CO2源汇关系的初步估算[J]. 科学通报, 1997, 42(9):953-956. doi: 10.3321/j.issn:0023-074X.1997.09.019

    XU Shengyou, JIANG Zhongcheng. Preliminary estimation of the relationship between karstification and the source sink of atmospheric greenhouse gas CO2 in China[J]. Chinese Science Bulletin, 1997, 42(9):953-956. doi: 10.3321/j.issn:0023-074X.1997.09.019
    [21] 罗健, 蒋勇军, 胡毅军, 李林立, 刘文. 亚高山表层岩溶泉域土壤溶蚀速率季节变化及碳汇量估算: 以重庆金佛山水房泉流域为例[J]. 中国岩溶, 2011, 30(4):443-448. doi: 10.3969/j.issn.1001-4810.2011.04.015

    LUO Jian, JIANG Yongjun, HU Yijun, LI Linli, LIU Wen. Seasonal changes of soil dissolution rate and estimations on carbon sequestration in the subalpine epikarst spring zone: A case study in the Shuifang Spring in Jinfoshan, Chongqing[J]. Carsologica Sinica, 2011, 30(4):443-448. doi: 10.3969/j.issn.1001-4810.2011.04.015
    [22] 曾成, 赵敏, 杨睿, 刘再华. 岩溶作用碳汇强度计算的溶蚀试片法和水化学径流法比较: 以陈旗岩溶泉域为例[J]. 水文地质工程地质, 2014, 41(1):106-111.

    ZENG Cheng, ZHAO Min, YANG Rui, LIU Zaihua. Comparison of karst processes-related carbon sink intensity calculated by carbonate rock tablet test and solute load method: A case study in the Chenqi karst spring system[J]. Hydrogeology & Engineering Geology, 2014, 41(1):106-111.
    [23] 闫伟, 曾成, 肖时珍, 蓝家程, 何江湖, 何春, 代林玉, 狄永宁. 湿润亚热带典型白云岩流域不同土地利用下的试片溶蚀速率及岩溶碳汇[J]. 地球与环境, 2021, 49(5): 529-538.

    YAN Wei, ZENG Cheng, XIAO Shizhen, LAN Jiacheng, HE Jianghu, HE Chun, DAI Linyu, DI Yongning. Dissolution rate and karst carbon sink of different land use in typical dolomite watershed with humid subtropical weather[J]. Earth and Environment, 2021, 49(5): 529-538.
    [24] 王文娟, 蓝芙宁, 蒋忠诚, 覃小群, 劳文科. 湖南大龙洞流域不同岩性不同土地利用类型条件下碳酸盐岩试片的溶蚀速率[J]. 中国岩溶, 2013, 32(1):29-33. doi: 10.3969/j.issn.1001-4810.2013.01.005

    WANG Wenjuan, LAN Funing, JIANG Zhongcheng, QIN Xiaoqun, LAO Wenke. Corrosion rate of carbonate tablet under diverse land use and lithology in the Dalongdong basin, Hunan[J]. Carsologica Sinica, 2013, 32(1):29-33. doi: 10.3969/j.issn.1001-4810.2013.01.005
    [25] 黄奇波, 覃小群, 刘朋雨, 蓝芙宁, 张连凯. 不同岩性试片溶蚀速率差异及意义[J]. 地球与环境, 2015, 43(4):379-385. doi: 10.14050/j.cnki.1672-9250.2015.04.001

    HUANG Qibo, QIN Xiaoqun, LIU Pengyu, LAN Funing, ZHANG Liankai. Dissolution on rate and it’s significance of different lithological tablets[J]. Earth and Environment, 2015, 43(4):379-385. doi: 10.14050/j.cnki.1672-9250.2015.04.001
    [26] 杨成华, 陈景艳, 丁访军. 普定县天龙山喀斯特森林群落植物多样性对比研究[J]. 贵州林业科技, 2012, 40(2):1-7.

    YANG Chenghua, CHEN Jinyan, DING Fangjun. Study on the plant diversity of karst forest community in Tianlong mountain of Puding County[J]. Guizhou Forestry Science and Technology, 2012, 40(2):1-7.
    [27] 俞锦标, 潘瑞鸿. 贵州省普定县碳酸盐岩沉积相带成岩作用与岩溶发育的研究[J]. 南京大学学报(自然科学版), 1982(3):787-798,823.

    YU Jinbiao, PAN Ruihong. Sedimentary facies zones, diagenesis and karst development of carbonate rocks in Pudding, west Guizhou[J]. Journal of Nanjing University(Natural Science), 1982(3):787-798,823.
    [28] 胡刚, 张忠华, 程安云, 刘立斌, 吴洋洋, 倪健. 黔中天龙山喀斯特次生林林分空间结构的量化与分析[J]. 地球与环境, 2017, 45(1):25-31. doi: 10.14050/j.cnki.1672-9250.2017.01.004

    HU Gang, ZHANG Zhonghua, CHENG Anyun, LIU Libin, WU Yangyang, NI Jian. Characterizing and analyzing stand spatial structure of a northern subtropical karst secondary forest in Tianlong mountain of central Guizhou Province, China[J]. Earth and Environment, 2017, 45(1):25-31. doi: 10.14050/j.cnki.1672-9250.2017.01.004
    [29] 习新强, 赵玉杰, 刘玉国, 王欣, 高贤明. 黔中喀斯特山区植物功能性状的变异与关联[J]. 植物生态学报, 2011, 35(10):1000-1008.

    XI Xinqiang, ZHAO Yujie, LIU Yuguo, WANG Xin, GAO Xianming. Variation and correlation of plant functional traits in karst area of central Guizhou Province, China[J]. Chinese Journal of Plant Ecology, 2011, 35(10):1000-1008.
    [30] Luo Mingming, ZHOU Hong, Liang Yongping, Chen Zhihua, Chen Rubing, Li Xinlong, Hamza Jakada, et al. Horizontal and vertical zoning of carbonate dissolution in China[J]. Geomorphology, 2018, 322:66-75.
    [31] 龚自珍, 黄庆达. 碳酸盐岩岩块野外溶蚀速度试验[J]. 中国岩溶, 1984(2):17-26.

    GONG Zizhen, HUANG Qingda. Field corrosion rate tests on carbonate rocks[J]. Carsologica Sinica, 1984(2):17-26.
    [32] 李永双, 范周周, 国辉, 周金星, 彭霞薇. 菌剂添加对不同树种根际土壤微生物及碳酸钙溶蚀的影响[J]. 中国岩溶, 2020, 39(6):854-862.

    LI Yongshuang, FAN Zhouzhou, GUO Hui, ZHOU Jinxing, PENG Xiawei. Effects of microorganisms agent addition on soil microbes in different rhizosphere soils and calcium carbonate dissolution[J]. Carsologica Sinica, 2020, 39(6):854-862.
    [33] 李强, 何媛媛, 曹建华, 梁建宏, 朱敏洁. 植物碳酸酐酶对岩溶作用的影响及其生态效应[J]. 生态环境学报, 2011, 20(12):1867-1871. doi: 10.3969/j.issn.1674-5906.2011.12.015

    LI Qiang, HE Yuanyuan, CAO Jianhua, LIANG Jianhong, ZHU Minjie. The plant carbonic anhydrase at karst area and its ecological effects[J]. Ecology and Environmental Sciences, 2011, 20(12):1867-1871. doi: 10.3969/j.issn.1674-5906.2011.12.015
    [34] SHE Min, Shou Jianfeng, SHEN Anjiang, PAN Liyin, HU Anping, HU Yuanyuan. Experimental simulation of dissolution law and porosity evolution of carbonate rock[J]. Petroleum Exploration and Development, 2016, 43(4):616-625.
    [35] 何宇彬, 金玉璋, 李康. 碳酸盐岩溶蚀机理研究[J]. 中国岩溶, 1984, 3(2):12-16.

    HE Yubin, JIN Yuzhang, LIKang. An experimental study of carbonate rock corrosion mechanism[J]. Carsologica Sinica, 1984, 3(2):12-16.
    [36] 黄奇波, 覃小群, 刘朋雨, 康志强, 唐萍萍. 半干旱区岩溶碳汇原位监测方法适宜性研究[J]. 吉林大学学报(地球科学版), 2015, 45(1):240-246. doi: 10.13278/j.cnki.jjuese.201501206

    HUANG Qibo, QIN Xiaoqun, LIU Pengyu, KANG Zhiqiang TANG Pingping. Applicability of karst carbon sinks calculation methods in semi-arid climate environment[J]. Joural of Jilin University(Earth Science Edition), 2015, 45(1):240-246. doi: 10.13278/j.cnki.jjuese.201501206
  • 期刊类型引用(2)

    1. 彭郭智,蒲俊兵. 关于溶蚀试片法计算溶蚀速率公式的一点认识. 中国岩溶. 2025(01): 113-123 . 本站查看
    2. 李源永,周华,丁访军,吴鹏,叶雨艳,崔迎春. 贵州典型岩溶区不同植被类型下碳酸盐岩溶蚀速率及影响因素. 地球与环境. 2024(04): 489-499 . 百度学术

    其他类型引用(5)

  • 加载中
图(7) / 表(5)
计量
  • 文章访问数:  1088
  • HTML浏览量:  881
  • PDF下载量:  72
  • 被引次数: 7
出版历程
  • 收稿日期:  2022-02-20
  • 刊出日期:  2022-06-25

目录

/

返回文章
返回