Theory and Method of Determining Erosion Rate in Limestone With 36Cl
-
摘要: 从理论上研究了方解石中36Cl的5种起源: (1)近地表Ca散裂和散裂成因的中子捕获; (2)由Ca产生的负μ介子捕获; (3)负μ介子捕获以后的中子发射; (4)由快μ介子引起的光致蜕变反应产生的中子; (5)由U裂变和(α, n)反应产生的中子.论述了36Cl在灰岩深度剖面上的4个分布特征: (1)36Cl浓度随深度增加而减少; (2)深度能够自然地区分由不同生成反应形成的36Cl贡献; (3)为了把测量的36Cl浓度与Ca的μ介子捕获及在35Cl上μ介子产生的中子捕获区别开来, 样品必须具有一定N(Ca)/ N(Cl)比值的范围; (4)甚至在最干旱条件下, 灰岩的侵蚀也是不可避免的, 并且在稳定状态下起重要作用.详细描述了36Cl的采样方法和AMS分析技术.研究了测定灰岩侵蚀速率的数学模型, 并且对一些灰岩表面侵蚀速率作了评价.Abstract: Five origins of 36Cl in calcite were studied theoretically: (1)near-surface production by Ca spallation and spallogenic neutron capture, (2)negative muon capture by calcium, (3)neutron emission following negative muon capture, (4)neutrons from photodisintegration reactions initiated by fast muons, (5) neutron from uranium fission and(α, n)reactions. Four characteristics of 36Cl distribution for deep profile in limestone have been discussed: (1)36Cl concentration is reduced with depth, (2)depth naturally separates 36Cl contributions from the various production reactions, (3)in order to partition measured 36Cl concentrations into contributions from muon capture on Ca and capture of muon-produced neutrons on 35Cl, the sample set must possess a range of Ca/Cl ratios, (4)limestone erosion is inevitable under even the most arid conditions and plays a significant role in determining steady-state concentrations.Sampling method of 36Cl and AMS analytical technique are described in detail.Mathematical model of determining erosion rate in limestone is studied and some surface erosion rate of limestone are evaluated.
-
Key words:
- surface erosion rate of limestone /
- 36Cl /
- accelerator mass spectrometry
-
自从1946年Libby及其同事们发现由于大气中宇宙射线核素相互作用可以生成宇宙成因同位素14C以后, 世界各国科学家便开始了大规模的宇宙成因同位素的研究.特别是20世纪70年代末用超灵敏加速器质谱计(AMS) 测定宇宙成因同位素获得成功, 开辟了这一研究领域的新纪元.其中应用10Be和26Al来研究岩石表面暴露年龄和侵蚀速率开展得比较广泛[1~11], 少数学者应用3He和21Ne研究了这一问题.应用36Cl开展这一研究的学者也比较多.但是, 其研究对象多数是年轻的低钾火山岩, 应用35Cl (n, γ)36Cl这一核反应原理测量原地产生的36Cl浓度, 只有少数学者把灰岩(方解石) 作为研究对象, 利用40Ca (μ-, α)36Cl和35Cl (n, γ)36Cl这两种核反应原理测定原地产生的36Cl浓度.
1984年, Kubik[1]应用AMS测定了灰岩及古生物群样品中的36Cl浓度.尽管灰岩中36Cl浓度很低, 澳大利亚学者Stone等[3, 4, 6]一直坚持用36Cl来研究灰岩表面的暴露年龄和侵蚀速率, 而且取得了显著的效果.
在我国还没有开展过用宇宙成因核素36Cl研究岩石表面侵蚀速率的工作, 特别是没有开展应用36Cl研究灰岩岩石表面侵蚀速率的工作.
2001年, 我们申请到了“用36Cl研究清江流域灰岩侵蚀速率及成壤作用和碳循环意义”的国家自然科学基金项目.此处发表的是我们的部分研究成果.
1. 方解石中36Cl起源
方解石中36Cl的起源有以下5种: (1) 近地表Ca散裂和散裂成因中子捕获.由高能量次级宇宙射线产生的Ca散裂是地壳顶部数米内方解石中36Cl的主要来源.在接近地表处, 散裂是次级中子的主要来源, 一旦热化, 次级中子就被35Cl捕获而产生36Cl. (2) 由Ca产生的负μ介子捕获.方解石中最重要的μ介子引起36Cl产生机制是Ca产生的负μ介子捕获.对于Ca的捕获而言, 由下式确定:
式中: *符号表示激发生成物核素, 退激时生成36Cl而放出α粒子:
将两种反应综合起来, 可表示为40Ca (μ-, α)36Cl. (3) 负μ介子捕获以后的中子发射.方解石成分中的中子产量可以由Ca、C和O核素以及每次捕获的平均中子产量所形成的μ介子捕获的数量和分布给出. (4) 由快μ介子引起的光致蜕变反应产生的中子.快速μ介子的减速和电子地下碰撞可以产生韧致辐射γ辐射, 这个辐射部分有足够的能量引起(γ, n) 核蜕变. (5) 由铀裂变和(α, n) 反应产生的中子.由U和Th系列衰变产生的238U和(α, n) 反应的自发裂变在岩石中可引起非宇宙成因的地下中子源.在富含Ca的灰岩中比花岗岩类岩石中这类中子要低得多, 因为花岗岩中U和Th含量高.一般来说, 由U裂变和(α, n) 反应生成的中子而产生的36Cl浓度可以忽略不计.但在深部由于宇宙成因的36Cl大大减少, 其贡献就要增大, 必须予以考虑.
2. 36Cl在灰岩深度剖面上的分布特征
一般有如下特征: (1) 灰岩中36Cl浓度随着深度增加而减少. (2) 深度的大小能够自然地区分由不同生成反应形成的36Cl贡献.散裂和散裂成因次级中子捕获在深度剖面的顶部500~600 g/cm2 (即2 m) 处其36Cl产量占优势; 500~6 000 g/cm2 (即2~20 m) 处μ介子捕获反应是36Cl的主要来源(但取决于Cl的浓度和侵蚀速率ε的大小); 在20 m以下, 由快μ介子产生的中子捕获形成的36Cl占优势; 在更大的深度上, 地下放射成因的36Cl产量超过了宇宙成因的36Cl产量. (3) 灰岩中的Cl浓度可影响其中的36Cl产量, Cl浓度低的样品可制约方解石中每一个停留μ-的36Cl产量, 而Cl浓度高的样品却制约了由μ介子捕获产生的平均中子产量.因此, 所采集的样品必须具有一定Ca/Cl比值的范围. (4) 在最干旱条件下, 灰岩的侵蚀也是不可避免的, 并且在稳定状态浓度下起重要作用.在大多数情况下, 不可能独立地用地质证据来判定研究地点上的侵蚀速率.但是, 这个困难能够用测量剖面顶部样品来解决.由于在剖面顶部数米内散裂超过了μ介子的产量, 并且超过这个深度范围散裂就大大减少.在这里36Cl浓度对于侵蚀来说是灵敏的, 但仅仅很微弱地取决于μ介子的反应.由于近地表的反应速率已被许多学者研究和校正[4, 7, 9], 所以接近剖面顶部的36Cl浓度可以用来测量岩石表面的侵蚀速率.
3. 采样方法和36Cl测试技术
3.1 采样方法
(1) 选择采样地点的原则.一般考虑以下几点: ①在纯灰岩或大理岩分布区采样, 岩石表面侵蚀很慢, 而且已经达到了稳定状况, 即侵蚀速率(ε) 趋于常数, 在地下一定深度上36Cl浓度趋于不变.这一条是最重要的公认的假设.②采样点岩石裸露, 或者覆盖层很薄, 至少现代看来覆盖层很薄.而且, 地形坡度要 < 20°.③取样深度从地表到2 m深处(~500 g/cm2), 一般不超过5 m (~1 400 g/cm2).④N (Ca) /N (Cl) 比值有一定范围, 即Cl浓度要低, Ca浓度要大.⑤通过地质研究和地质年龄测定, 要知道连续侵蚀的年龄, 即要知道岩石表面的暴露年龄.⑥采样方便.采样地点通常选择在采石场, 道路切面或溶洞内.
(2) 采样类型和测定方法.①36Cl样品, 用AMS测定; ②Cl浓度测定样品, 用离子色谱测定; ③Ca, Mg, B, Sm, Gd, Fe, U, Th样品, 用ICP和ICP-MS测定; ④灰岩密度(ρ) 测定样品; ⑤14C年龄测定样品用液闪仪或AMS测定; ⑥不平衡U系年龄测定样品, 用α-谱仪或热离子质谱计测定.
(3) 采样技术.①用专门设计的水平钻取样.②在同一水平上可钻2~3个孔, 如果采样点不规则, 可在1.5 m×2.5 m范围内移动, 采样量1~2 kg.③采样点的位置和高程都用经纬仪和红外光学测定仪测量, 精度达±1 m.采样场剖面若不垂直则要测量其倾角.④采样间隔靠近地表时要小一点, 样品点要密一些, 随着深度加大, 增大采样间隔.为了满足侵蚀速率计算要求, 一般采样深度0~5.0 m即可.因为由μ介子反应产生的36Cl比26Al和10Be更丰富.⑤良好的暴露年龄和侵蚀速率的评价可在下述情况下取得.自从岩石形成以来, 表面的总侵蚀损失小于50~100 g/cm2 (即0.2~0.5 m).在这一深度上, 热中子流在地下可达到最大值.随着侵蚀加深, 年分辨降低, 然而侵蚀速率可仍旧被精确地测定.⑥方解石中36Cl积累可以测定100~100×106 a的灰岩表面暴露年龄, 同时可以确定1~1 000 μm/a的侵蚀速率.⑦本方法可以研究多阶段的侵蚀历史, 例如末次冰期以前侵蚀速率ε1和目前的侵蚀速率ε2.
3.2 36Cl测试技术[10]
目前, 用AMS测定36Cl的极限已达5×103原子/g或者N (36Cl) /N (Cl) 比值已达5×10-15.
测定对象有: 玄武岩中的钙长石、方解石、基岩中的斜长石和辉石类等富含Ca的矿物.
(1) AgCl样品制备.①压碎样品.压碎样品, 过筛, 使用150~250μm的样品.实验证明, 使用150~250 μm样品或250~500 μm样品对36Cl的测量结果没有明显的差别.②将粉碎样品混合均匀.用四分法取样备用.制备AgCl样品时通常需要5~15 g样品(一般不超过20 g).③除掉大气污染的36Cl.首先用去离子水冲洗样品, 再用2%的HNO3洗涤两次, 并且过滤.对方解石来说, 在室温条件下, 要加入足够的2%的HNO3, 以便溶解掉颗粒表面最外层的10%~15%.④溶解样品.用逐渐加入2 mol/L HNO3的方法选择性地溶解样品, 并保持样品溶液总酸度在0.5 mol/L以下.样品溶解后, 过滤, 把溶解的Ca (NO3) 2收集起来, 作氯化物和36Cl测定.⑤加入氯化物载体.为了增加最终的样品量和减少36Cl浓度计算中的氯化物浓度测量误差的扩散, 要加入1 mg左右的氯化物载体(NH4Cl).最佳的选择是加入如下所述的载体, 即N (Cl载体) /N (Cl岩石) 的比值要与来自Ca和35Cl的36Cl产量成比例.通常, 加入载体的总量要比AMS测量中保持高的N (36Cl) /N (Cl) 比值的最优条件要低2~3倍, 以便保持Cl浓度的灵敏度, 遵循由中子捕获而引起的36Cl产量, 以便进行校正.⑥加入AgNO3使AgCl沉淀出来.⑦除硫.⑧纯化AgCl沉淀.⑨AgCl样品保存.通常要求给出4~8 mg的AgCl供AMS测定使用, 实际上装入靶的样品只要1~2 mg.AgCl固体样品要用铝箔包好, 置于烘箱中, 在100 ℃条件下保持24 h以上.⑩按上述步骤制备空白样品.
(2) 加速器质谱计(AMS) 测定.为了取得准确可靠的测试结果, 通常要求延长30 min计数时间.要求最低检测极限为5×103原子/g, 或N (36Cl) /N (Cl) 比值< 5×10-15.
4. 数学模型研究及侵蚀速率计算
4.1 数学模型研究(校正策略)
数学模型研究[4]的目的是修复影响μ介子反应引起的36Cl产量中的36Cl浓度剖面.为此, 要求有一个36Cl产量(P) 与深度相关的模型, 它包括所有有影响的反应和剖面暴露历史的知识.
在灰岩侵蚀表面以下, 36Cl浓度可聚集到由侵蚀速率(ε) 支配的稳定状态, 在深度z处稳定状态的36Cl浓度由下式计算:
式中: ρ为灰岩的密度; λ为36Cl的衰变常数; t为岩石表面暴露年龄.
4.2 侵蚀速率计算和评价
在确定侵蚀速率时公认的假设[2, 4]是, 岩石长期暴露地表后, 由于宇宙射线的照射, 地下产生的同位素浓度(例如36Cl, 10Be, 26Al等) 达到稳定状态, 此时, 侵蚀速率也为常数.在这种情况下, 侵蚀速率可以用单个同位素测量结果来确定, 并且可提供已知的产生剖面.
对于仅由散裂反应产生的36Cl同位素来说, 使用指数产生剖面、稳定状态浓度和侵蚀速率之间的相互关系可以计算侵蚀速率.这个直截了当的模型是一个与10Be和26Al情况很相似的模型.在大多数侵蚀速率的研究中, 被众多的研究者使用过.但是, 由于忽略了由μ介子引起的36Cl产量, 常可导致较大的误差.
在由热中子捕获引起的36Cl同位素产量的情况下, 地表以下最大产生速率在50~100 g/cm2 (即0.2~0.5 m) 处, 在稳定状态下, 36Cl同位素浓度和侵蚀速率关系式由Liu等[7]给出.
Stone等[3]根据灰岩表面中36Cl浓度测量结果计算了澳大利亚一些地区的侵蚀速率.在干旱的Nullarbor Plain测定的灰岩侵蚀速率为4.5 μm/a, 这一结果与用26Al和10Be对含石英岩石(砂岩) 的测量结果很相似.在东部高原潮湿温和而适中的地区Strickland Range测定的侵蚀速率为184 μm/a.在Territory北部Tindal灰岩区测定的侵蚀速率为14 μm/a, 似乎与地区的季风气候有关.在澳大利亚东部高地用36Cl测定的侵蚀速率比现有的地貌评价要偏大, 但是与欧洲北部和北美的冷—温和地区所测定的侵蚀速率更加相似.
-
[1] Kubik P W. Accelerator mass spectrometry of 36Cl in limestone and some paleontological samples using completely stripped ions[J]. Nucl Instr Meth, 1984. B5: 326-330. [2] 李义连, 王焰新, 张江华, 等. 娘子关泉域岩溶水硫酸盐污染的地球化学模拟分析[J]. 地球科学——中国地质大学学报, 2000, 25(5): 467-471. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200005004.htmLI Y L, WANG Y X, ZHANG J H, et al. Geochemical modeling of sulfate pollution in karst water of Niangziguan region, Shanxi Province, China[J]. Earth Science — Journal of China University of Geosciences, 2000, 25(5): 467-471. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200005004.htm [3] Stone J, Fifield L K, Allan G L, et al. Limestone erosion measurement with cosmogenic chlorine-36 in calcite preliminary results from Australia[J]. Nucl Instr Meth, 1994, B92: 311-316. [4] Stone J, Allan G L, Fifield L K, et al. Cosmogenic chlorine-36 from calcium spallation[J]. Geochim Cosmochim Acta, 1996, 60: 679-692. doi: 10.1016/0016-7037(95)00429-7 [5] 黄俊华, 胡超涌, 周群峰. 湖北清江和尚洞石笋的高分辨率碳氧同位素及古气候研究[J]. 地球科学——中国地质大学学报, 2000, 25(5): 505-509. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200005013.htmHUANG J H, HU C Y, ZHOU Q F. High-resolution carbon and oxygen isotope records from stalagmite and palaeoclimate in Heshangdong cave, Qingjiang, Hubei Province[J]. Earth Science—Journal of China University of Geosciences, 2000, 25(5): 505-509. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200005013.htm [6] Stone J, Evans J M, Fifield L K, et al. Cosmogenic chlorine-36 production in calcite by muons[J]. Geochim Cosmochim Acta, 1998, 62(3): 433-454. doi: 10.1016/S0016-7037(97)00369-4 [7] Liu B, Phillips F M, Fabryka-Martin J T, et al. Cosmogenic 36Cl accumulation in unstable landforms, 1. Effects of the thermal neutron distribution[J]. Water Resour Res, 1991, 30: 3115-3125. [8] 祁士华, 盛国英, 傅家谟, 等. 鼎湖山自然保护区大气气溶胶中的PAHs[J]. 地球科学——中国地质大学学报, 2001, 26(1): 83-87. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200101015.htmQI S H, SHENG G Y, FU J M, et al. PAHs in aerosols at Dinghushan natural protection zone[J]. Earth Science — Journal of China University of Geosciences, 2001, 26 (1): 83-87. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200101015.htm [9] Phillips F M, Zreda M G, Flinsch M R, et al. A reevaluation of cosmogenic 36Cl production rate in terrestrial rocks[J]. Geophys Res Lett, 1996, 23: 949-951. doi: 10.1029/96GL00960 [10] 刘存富, 王佩仪, 周炼, 等. 河北平原第四系地下水36Cl年龄研究[J]. 水文地质工程地质, 1993, 6: 35-38. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG199306009.htmLIU C F, WANG P Y, ZHOU L, et al. Cl dating of Quaternary groundwater in Hebei plain[J]. Hydrogeology and Engineering Geology, 1993, 6: 35-38. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG199306009.htm [11] 杨银湖, 黄正发. 京珠线湖北省南段岩溶地质问题与勘察对策[J]. 地球科学——中国地质大学学报, 2001, 26 (4): 361-364. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200104007.htm -

计量
- 文章访问数: 3458
- HTML全文浏览量: 123
- PDF下载量: 3
- 被引次数: 0