王海静, 张金流, 刘再华, 台永东

1)中国科学院地球化学研究所环境地球化学国家重点实验室, 贵州贵阳 550002;2)黄龙国家级风景名胜区管理局, 四川松潘 623300

基于水化学和同位素特征的四川黄龙沟泉群分类研究

王海静1), 张金流1), 刘再华1), 台永东2)

1)中国科学院地球化学研究所环境地球化学国家重点实验室, 贵州贵阳 550002;2)黄龙国家级风景名胜区管理局, 四川松潘 623300

本文采用水化学和同位素方法对四川黄龙沟沿途出露的7个泉点进行了分析。结果表明, 泉水水化学和同位素的时空变化反映了CO2逸出、钙华沉积和蒸发效应等诸多因素的共同影响, 是由这些泉水处在四个水循环转化段决定的。根据水化学和同位素特征可将这些泉水划分为三种不同的类型: 深部泉、表生泉和转化泉。这些认识将为四川黄龙沟景观水资源的管理和钙华景观的保护提供科学依据。

泉水; 水化学; 同位素; 景观水资源; 钙华景观保护; 四川黄龙

泉水是地下水的天然露头, 其地球化学特征反映了泉域的地质和气候背景条件以及地表地下水文过程。作为世界自然遗产的四川黄龙沟, 随着资源开发和旅游环境建设的发展, 人类活动影响下的水质已发生某些变化, 如钙华沉积减慢、砂化和变黑等问题相继出现(刘再华等, 2009; Wang et al., 2010),因此, 自然和人为影响背景的研究需要进一步深入:一方面是区域地质与地理背景, 包括岩石、水系和气候等方面的特征; 另一方面是监测和探明人类活动对自然环境的影响。目前, 区域内水化学和同位素的系列资料较少, 黄龙沟的水文地球化学研究工作还需进一步系统深入。本文主要研究了不同地貌背景条件下的黄龙沟泉水水化学和同位素的差异,以期揭示黄龙沟地表、地下水水化学主要是受何种因素控制和影响的, 以便为维系钙华景观的水资源管理提供科学依据。

1 研究区背景

黄龙沟补给区岩性总体以碳酸盐岩为主, 约占44%, 其次为板岩, 约占 27%, 砂岩及千枚岩各占18%和 10%, 上游的石炭-二叠碳酸盐岩层为钙华沉积的补给源区, 形成的钙华主要沉积在中下游三叠系和志留系的板岩、砂岩之上(卢国平等, 1992)。整个山体受垂直气候带的影响, 区内植被类型复杂,基本上包括了阔叶林、针叶林、灌丛、高山草甸及高寒荒漠等各种类型。土壤分布也具有明显的垂直分带性。该区属于典型的高原温带-亚寒带季风气候,冬季漫长, 年平均气温只有 1.1℃, 昼夜温差大, 降雨量 758.9 mm/a, 且集中在 5～9月份, 约占全年降水量的 70%。地表融雪水和泉水在沿途流动过程中进行着多次的地表地下转换, 形成了具有不同特征的复杂的地表地下水循环系统(李前银等, 2009)。

2 研究方法

沿黄龙沟自源头泉口至近涪江, 选取 7个主要实验观测泉点(图1)。主要采用仪器监测、现场滴定和样品室内测试相结合的方法。仪器使用德国WTW公司生产的Multiline-P3多参数水质监测仪, 现场测量泉水的pH、水温和电导率, 分辨率分别达到0.01 pH单位、0.1℃和1 μs/cm。仪器在测量之前进行校正, 其pH值用4和7两种标准缓冲溶液进行校正,而仪器的电导率用1412 μs/cm标准校正液校正。现场滴定使用德国Merck公司生产的碱度计和硬度计测定水的 HCO3-和 Ca2+浓度, 其精度分别为 6 mg/l和2 mg/l(王海静等, 2009; Wang et al., 2010)。

用于水化学分析和同位素测定的样品拿回中国科学院地球化学研究所环境地球化学国家重点实验室测定其 K+、Na+、Mg2+、 SO42-、Cl-含量, 检出下限是0.001 ppm; 碳、氢稳定同位素的测定误差分别小于0.15‰和0.5‰。

方解石和白云石饱和指数(SIC和SID), 以及CO2分压(pCO2)的计算利用WATSPEC软件进行(Wigley,1977)。

3 测试结果及其分析讨论

3.1 泉水水化学的特征及其季节变化规律分析

黄龙沟7个泉点的水化学特征总结于表1, 据此可将这些泉水分为3大类: 深部泉(1号泉)、表生泉(7 号泉)和转化泉(2～6 号泉)。

三大类泉水中表生泉的水温和 pCO2最低, pH值最高, 最低的电导率和 Ca2+、HCO3-浓度, 因其CO2主要来源于土壤和大气(表 1, 图 2), 深部泉因其CO2主要来源于地球深部(刘再华等, 2000, 2009),故具有最高的pCO2、最低的pH值、最高的电导率和 Ca2+、HCO3-浓度; 而转化泉主要是表生泉(地表融雪水)和深部泉水的混合, 故其水化学特征介于表生泉和深部泉水之间, 如表1。

图1 黄龙沟钙华景观形成的平/剖面示意图Fig. 1 Geological plan and section of Huanglong Ravine

表1 黄龙沟泉水水化学特征Table 1 Hydrochemical features of the Huanglong Springs

图2 黄龙沟泉水的水化学季节变化Fig. 2 Seasonal hydrochemical variations of the springs at Huanglong

从图2还可以看出深部泉的水温常年较为稳定,表生泉由于受温度较低的融雪水补给的影响使其水温比其他泉水偏低, 它与转化泉的水温均呈现出夏季升高冬季降低的季节变化特点。表生泉的 pCO2夏季高, 表现出表生岩溶系统相似的特征, 与土壤中CO2含量的变化有关(Liu et al., 2007)。深部泉则呈现冬季升高, 夏季降低的趋势, 可能是冬季来自于深部的 CO2较多, 受雨水稀释作用影响小的缘故, 还与冬季气温低, 水中的 CO2气体逸出速度较夏季缓慢有关。二次转化泉受其影响为主, 也表现出相似的季节变化特点。表生泉在夏季随着pCO2的升高, pH降低, 而深部泉的pH变化相对稳定, 在误差范围内。转化泉主要受温度的影响, 夏季 CO2逸出增加, 导致pH出现升高趋势。受其影响转化泉的 SIc也出现夏季升高的特点。深部泉由于较少受外界因素影响, SIc受水温和pH值的制约也相对稳定。表生泉的方解石饱和指数则受其夏季 pCO2升高、pH降低的影响而出现降低趋势。水中的主要阳离子 Ca2+的含量则在深部泉中为最高, 夏季受到雨水的少量稀释, 而它在表生泉中的含量在春季时偏高应该是与土壤中 CO2含量的增加有关, 夏季降低则主要是稀释效应影响为主(Liu et al., 2004, 2007)。由于沿途钙华的沉积使转化泉中 Ca2+的含量远低于深部泉, 季节上的变化则主要受稀释作用和较强的沉积作用而呈现出夏季降低的特点。

3.2 泉水碳和氢同位素的特征及其季节变化规律分析

由于不同碳库的δ13C差异较大, 所以利用碳同位素能很好的示踪泉水中碳的来源及演化(Fritz et al., 1989)。表生泉的δ13C最低, 反映出其碳来源的生物成因特征, 而其δ13C在雨季升高, 可能是因为有深部的重碳混入导致。深部泉水因CO2主要来源于深部非生物成因(刘再华等, 2000, 2009), 故其碳同位素值偏高, 且随着深部泉水出露后 CO2逸出,向下游流动过程中钙华大量沉积, 产生碳同位素的分馏作用, 导致下游转化泉水的碳同位素值进一步增加(图 3)。

雨水的δD值在春季偏高, 反映了该研究区春季降雨偏少而蒸发强烈的特点。主要受雨水补给的表生泉的δD值在夏季降低, 但比雨水滞后1～2个月,反映了表生泉水的循环速率还是较为缓慢的。转化泉的δD和深部泉的变化规律较为相似, 但深部泉的δD较低, 反映了其补给高程可能是最大的, 转化泉的δD较高应该主要是受蒸发作用的影响所致。

3.3 黄龙沟泉群的进一步分析

李前银等(2009)的研究表明, 黄龙景区水循环系统由钙华源泉岩溶地下水系统、地表水系统和地表、地下水转化系统三部分组成。若无沿途众多泉水的补给, 黄龙沟内溪水水化学应随源头泉水的距离增加单调增加或降低, 但每一循环段出露的转化泉打破了各项地球化学指标空间上的规律性(图 4、图 5), 不仅构成下一循环段水量上的“源泉”, 也对沟内溪水中各种离子再次进行补充; 每一循环段钙华的生长都随景观水流路径的延长、碳酸钙的不断析出而减弱, 景区最为艳丽壮观的五彩池、争艳彩池、金沙铺地和迎宾池都位于各循环段邻近“源泉”的区域。随着沉积作用的不断进行, 水中碳酸钙含量逐渐减少, 到黄龙沟口已与地表水相差无几, 相应地各循环转化段也显示出由泉口往下钙华沉积作用由弱至强再减弱的趋势, 如到每个循环阶段的最后一段, 钙华沉积已非常缓慢, 且大部分出现钙华砂化变黑等退化趋势。

图3 黄龙泉水碳、氢同位素的季节变化Fig. 3 Seasonal variations in carbon and hydrogen isotopic compositions of the springs at Huanglong

图4 黄龙沟泉水的水化学空间变化(图中1～7为泉水编号)Fig. 4 Spatial hydrochemical variations of the springs at Huanglong(1 to 7 of the fig is the number of spring)

图5 黄龙沟泉水同位素的空间变化(图中1～7为泉水编号)Fig. 5 Spatial variations in carbon and hydrogen isotopic compositions of the springs at Huanglong(1 to 7 of the fig is the number of spring)

下面, 根据泉水的地貌位置和泉水的地球化学特征, 探讨其空间分布意义, 如表2。

黄龙景区源头岩溶地下水(泉群I)出露地面后与地表融雪水混合形成沟内主要的景观水源, 在流动过程中不断沉积CaCO3形成钙华景观并沿途大量漏失转化为地下水,又在一定条件下以泉的方式出露形成转化泉, 周而复始, 在核心景区形成四个转化段,从而构成一个特殊的水循环及景观演化系统(图1)。

(1)五彩池—接仙桥转化段

1号泉群和7号泉出露后与地表融雪水混合, 一路漏失, 至马蹄海附近完全转入地下。漏失的水流在钙华与下伏冰碛层接触带汇集, 最终形成的泉水在接仙桥一带的沟谷两侧边缘出露, 从而完成了景观水的第一次转化。代表景观: 五彩池。随着 1号深部泉的出露, 水流向下的流动过程中, 水中 pCO2远大于大气中的pCO2, 导致水中CO2大量快速逸出,CO2分压急剧降低, 再加上地表融雪水的混入, 导致沿途水的 pH值显著升高, 方解石饱和指数也有较大的增加, 钙华大量沉积, 钙和重碳酸根的含量也有显著下降(表1, 1号泉和2号泉间的变化, 图4)。

(2)接仙桥—争艳池转化段

第一阶段漏失后形成的地下水在接仙桥2号泉点以下汇入 3号泉形成沟状水流, 流至宿云桥附近分流。东侧形成彩池群及钙华滩流并沿途漏失, 部分漏失水流在争艳彩池附近的三岔路口形成4号泉;西侧水流则形成艳丽无比的争艳彩池, 然后很快漏失, 并在彩池末端形成 5号泉, 从而完成了景观水的第二次转化。代表景观: 争艳池。第一阶段漏失的水量在此又重新得到补充, 由于泉水的出现, 水中的各项地球化学指标与第一段末段相比均有明显变化(图 4、图 5), 在经过地下一段距离的流动过程后, CO2得到了一定程度的补充, 在2、3号泉水出露处泉水水温较为稳定, 与第一循环段末端的地表水流相比, 2、3号转化泉的pCO2升高, pH值降低, 方解石饱和指数降低。

(3)争艳池—莲台飞瀑转化段

4号泉出露后与第二阶段的东侧剩余水流汇合,加上 5号泉, 共同沿平坦而倾斜的沟面倾泻而下,与此同时景观水流也在不断漏失, 到莲台飞瀑时,下渗水流在莲台飞瀑下形成 6号泉, 从而完成了景观水的第三次转化。代表景观: 金沙铺地。此阶段钙华沉积形态以钙华滩为主, 形成的钙华滩流是目前世界上发现的同类地质构造中, 状态最好、面积最大、距离最长的地表钙华滩流。4和 5号泉是这一循环段的主要补给水源, 在继承了两个主要补给泉水的各项特征的基础上, 在长约1.3 km的距离内水的地球化学特征在空间上体现出相对较为有规律性的变化: 随着水流的向下流动, 沿途pH值呈升高趋势, pCO2和电导率降低(图4)。

表2 黄龙沟泉群类型的划分及其基本特征Table 2 Classification of spring groups at Huanglong Ravine

(4)莲台飞瀑—黄龙沟口转化段

6号泉与该段地表水混合后向下游流动, 在形成和养护下游钙华景观的同时,一路漏失。漏失的水流多数以潜流的方式流入涪江, 从而完成了景观水的最后一次转化。代表景观: 迎宾池。在第三循环段的末段, 水量已大大减少, 水中的离子含量也大为降低, 钙华的沉积能力变得微弱, 此时 6号泉水的出露及时地补充了下游的水量, 并一定程度上增加了水中的pCO2及钙和重碳酸根的含量, 对最下游的迎宾彩池的形成和养护起到了非常重要的作用。

4 结论

黄龙沟地区泉水在地表和地下进行了多次转换,其过程主要受地形地貌的控制。通过对沟内沿途出露的7个泉点水化学和同位素时空变化特征(反映了CO2逸出、钙华沉积、蒸发效应等诸多因素的影响)的分析, 将泉水划分为三种类型: 深部泉、表生泉和转化泉, 它们受控于沟内四个水循环转化段: 五彩池—接仙桥转化段、接仙桥—争艳池转化段、争艳池—莲台飞瀑转化段和莲台飞瀑—黄龙沟口转化段。

总之, 本研究通过水化学和同位素时空变化规律的分析, 对于景区地表地下水的转换过程有了更清晰的认识, 这对于下一步可能进行的各景观段地表水的调配和防渗工程提供了重要参考, 同时也为有效防治景观水的污染提供了科学依据, 因为只有保持黄龙水资源量和质的相对稳定, 才能更好地保护黄龙景区的旅游资源。

李前银, 范崇荣. 2009. 黄龙景区水循环系统与景观演化研究[J].水文地质工程地质, (1): 108-112.

刘再华, 田友萍, 安德军, 王海静, 唐淑, 张金流, 孙海龙, 刘彦,张清明. 2009. 世界自然遗产——四川黄龙钙华景观的形成与演化[J]. 地球学报, 30(6): 841-847.

刘再华, 袁道先, 何师意, 张美良, 张加桂. 2000. 地热 CO2-水-碳酸盐岩系统的地球化学特征及其 CO2来源[J]. 中国科学(D辑), 30(2): 209-214.

卢国平, 李行健. 1992. 四川黄龙自然风景区冷水型钙华地表堆积地貌研究[J]. 成都地质学院学报, 19(4): 55-64.

王海静, 刘再华, 曾成, 刘香玲, 孙海龙, 安德军, 唐淑, 张清明.2009. 四川黄龙沟源头黄龙泉泉水及其下游溪水的水化学变化研究[J]. 地球化学, 38(3): 307-314.

FRITZ P, FONTES J C, FRAPE S K, LOUVAT D, MICHELOT J L,BALDERER W. 1989. The isotope geochemistry of carbon in groundwater at Stripa[J]. Geochim Cosmochim Acta, 53(8):1765-1775.

LI Qian-yin, FAN Chong-rong. 2009. A study of water circulation system and landscape evolution of Huanglong scenic spot[J].Hydrogeology and Engineering Geology, (1): 108-112(in Chinese with English abstract).

LIU Zai-hua, GROVES C, YUAN Dao-xian, MEIMAN J. 2004.South China Karst Aquifer Storm-Scale Hydrogeochemistry[J].Ground Water, 42(4): 491-499.

LIU Zai-hua, LI Qiang, SUN Hai-long, WANG Jin-liang. 2007.Seasonal, diurnal and storm-scale hydrochemical variations of typical epikarst springs in subtropical karst areas of SW China:soil CO2and dilution effects[J]. Journal of Hydrology,337(1-2): 207-223.

LIU Zai-hua, TIAN You-ping, AN De-jun, WANG Hai-jing, TANG Shu, ZHANG Jin-liu, SUN Hai-long, LIU Yan, ZHANG Qing-ming. 2009. Formation and evolution of the travertine landscape at Huanglong, Sichuan, one of the world natural heritages[J]. Acta Geoscientica Sinica, 30(6): 841-847(in Chinese with English abstract).

LIU Zai-hua, YUAN Dao-xian, HE Shi-yi, ZHANG Mei-liang,ZHANG Jia-gui. 2000. Geochemical features of the geothermal CO2-watercarbonate rock system and analysis on its CO2sources[J]. Science in China, 30(2): 209-214(in Chinese with English abstract).

LU Guo-ping, LI Xin-jian. 1992. A study on cold-water travertine surface depositional landforms in Huanglong scemic spot, Sichuan Province[J]. Journal of Chengdu College of geology,19(4): 55-64(in Chinese with English abstract).

WANG Hai-jing, LIU Zai-hua, ZENG Cheng, LIU Xiang-ling, SUN Hai-long, AN De-jun, TANG Shu, ZHANG Qing-ming. 2009.Hydrochemical variations of Huanglong Spring and the stream in Huanglong Ravine, Sichuan Province[J]. Geochimica, 38(3):307-314(in Chinese with English abstract).

WANG Hai-jing, LIU Zai-hua, ZHANG Jin-liu, SUN Hai-long, AN De-jun, FU Ru-xian, WANG Xiao-ping. 2010. Spatial and temporal hydrochemical variations of the spring-fed travertine-depositing stream in the Huanglong ravine, Sichuan,SW China[J]. Acta Carsologica, 39(2): 247-259.

WIGLEY T M L. 1977. WATSPEC: a computer program for determining equilibrium speciation of aqueous solutions[J]. British Geomorphological Research Group Technical Bulletin, 20: 1-46.

Classification of Spring Groups Based on Hydrochemical and Isotopic Features of Huanglong Ravine

WANG Hai-jing1), ZHANG Jin-liu1), LIU Zai-hua1), TAI Yong-dong2)
1)State Key Laboratory of Environmental Geochemistry, Institute of Geochemistry, Chinese Academy of Sciences,Guiyang, Guizhou550002;2)Huanglong National Scenic Spot Administration, Songpan, Sichuan623300

Hydrochemical and isotopic methods were used to analyze 7 springs issued along the Huanglong Ravine.It is found that spatial-temporal variations in hydrochemistry and isotopes of the 7 springs are related to the combined influence of CO2degassing, travertine deposition and evaporation, which are determined by four stages of water cycle. Three different spring groups were then identified based on these hydrochemical and isotopic features. These findings may provide scientific basis for the management of the landscape water resources and the protection of the travertine landscape.

spring; hydrochemistry; isotope; landscape water resource; travertine landscape protection;Huanglong Ravine in Sichuan

K928.9; P641.1; P597

A

10.3975/cagsb.2011.03.12

本文由中国科学院百人计划项目和国家自然科学基金项目(编号: 40872168)联合资助。

2011-02-27; 改回日期: 2011-04-05。责任编辑: 闫立娟。

王海静, 女, 1980年生。博士研究生, 主要从事岩溶环境地球化学研究。E-mail: whj0707@sina.com。

刘再华, 男, 1963年生。研究员, 博士生导师。长期从事岩溶作用动力学与全球变化研究。通讯地址: 550002, 贵州省贵阳市观水路46号。电话: 0851-5892338。E-mail: liuzaihua@vip.gyig.ac.cn。