康志強(qiáng) , 陳 駿, 袁道先, 何師意, 鄧 艷, 陳 旸,劉媛媛, 姜光輝, 張勤軍

1)表生地球化學(xué)教育部重點(diǎn)實(shí)驗(yàn)室/南京大學(xué)地球科學(xué)與工程學(xué)院, 江蘇南京 210093;

2)中國(guó)-東盟地學(xué)合作中心(南寧), 廣西南寧 530023;

3)中國(guó)地質(zhì)科學(xué)院巖溶地質(zhì)研究所, 廣西桂林 541004;

4)廣西壯族自治區(qū)地質(zhì)調(diào)查院, 廣西南寧 530023

近些年來(lái), 我國(guó)在石漠化治理方面做出了諸多的探索研究(蔣忠誠(chéng)等, 2009; 曹建華等, 2016; 李強(qiáng)等, 2017; 何霄嘉等, 2019; 肖林穎等, 2021), 形成了一套基本成熟的石漠化治理模式(曹建華等, 2008;Jiang et al., 2014), 取得了明顯的效果(蔣忠誠(chéng)等,2008; 杜文鵬等, 2019; 王克林等, 2019), 石漠化蔓延的勢(shì)頭也得到了有效的遏制(國(guó)家林業(yè)和草原局,2018), 且持續(xù)向好發(fā)展。2019年 2月 12日, 美國(guó)航天局(NASA)發(fā)布消息稱(chēng), 由于中國(guó)和印度的人類(lèi)活動(dòng)的影響, 在過(guò)去20年中, 地球變得越來(lái)越綠了, 尤其在中國(guó)最為顯著, 占到全球新增綠葉面積的25%(Chen et al., 2019)。巖溶作用過(guò)程是由地下水驅(qū)動(dòng)的地質(zhì)循環(huán)過(guò)程, 是大氣圈、巖石圈、生物圈、水圈相互作用最密切的地帶(袁道先, 1993, 1995;Yuan, 1997b), 因而也是巖溶區(qū)地球關(guān)鍵帶的最主要構(gòu)成部分。地球關(guān)鍵帶空間界限范圍為上到植被冠層, 下到地下水含水層底部, 包含近地表的生物圈、大氣圈、整個(gè)土壤圈以及水圈和巖石圈近地表部分(National Research Council, 2001; Lin et al.,2005; Sparks, 2005; Lin, 2010)。尤其在巖溶區(qū), 關(guān)鍵帶中水循環(huán)過(guò)程顯得格外重要, 其不僅能促進(jìn)巖溶作用過(guò)程本身, 還進(jìn)一步促進(jìn)了大氣、植被和碳酸鹽巖之間的碳循環(huán)過(guò)程。研究表明, 巖溶區(qū)森林植被可以增強(qiáng)巖溶小流域的降水-徑流過(guò)程(Mcculloch and Robinson, 1993), 同時(shí), 植被自身生長(zhǎng)過(guò)程還會(huì)影響巖溶作用過(guò)程(曹建華和袁道先, 1999; 何師意等, 2001; 楊安學(xué)和彭云, 2007)。

本項(xiàng)研究在我國(guó)西南巖溶區(qū)選擇了2個(gè)不同植被覆蓋條件的巖溶關(guān)鍵帶中的地下水系統(tǒng)為研究對(duì)象, 通過(guò)對(duì)巖溶水循環(huán)特征及水化學(xué)特征的對(duì)比,研究森林植被對(duì)地下水循環(huán)過(guò)程和水化學(xué)組分的影響。

1 研究區(qū)概況及取樣點(diǎn)布置

選擇的兩個(gè)研究點(diǎn)地域相近, 均處于亞熱帶季風(fēng)氣候區(qū), 具有高溫多雨、降雨集中、濕度大等特點(diǎn)。但其所處的地貌單元和植被覆蓋條件不同, 其巖溶關(guān)鍵帶的結(jié)構(gòu)和功能也具有明顯的不同。其中貴州茂蘭國(guó)家級(jí)自然保護(hù)區(qū)位于貴州省東南部, 地表標(biāo)高在 450～1100 m之間, 處于云貴高原向廣西盆地過(guò)渡的斜坡地帶。森林覆蓋率在90%以上。出露地層為石炭系中統(tǒng)灰?guī)r、白云巖。地貌以峰叢洼地、谷地為主。該區(qū)人口密度較小, 受人類(lèi)活動(dòng)的干預(yù)程度甚微。板寨地下河發(fā)育于保護(hù)區(qū)內(nèi), 為獨(dú)立的巖溶地下水系統(tǒng)(圖 1-R1), 僅受大氣降水的補(bǔ)給。在石炭系中統(tǒng)碳酸鹽巖和下統(tǒng)砂頁(yè)巖的地層界面處發(fā)育巖溶管道, 系統(tǒng)內(nèi)地下水在地下河出口處全部排向地表。板寨自動(dòng)化監(jiān)測(cè)站建于地下河出口下游, 其匯水面積共計(jì)27.26 km2。

柳州市農(nóng)耕區(qū)位于廣西壯族自治區(qū)中北部柳州市區(qū)周?chē)? 處于廣西平原桂中盆地內(nèi), 標(biāo)高在90～150 m之間。為峰叢谷地地貌, 地勢(shì)較平坦, 出露地層為石炭系至二疊系碳酸鹽巖。谷底內(nèi)為主要的農(nóng)耕區(qū), 人口分布較為稠密, 農(nóng)業(yè)活動(dòng)較為發(fā)達(dá)。山上植被稀疏, 以灌木為主。官村地下河位于峰叢洼地向峰叢谷地的過(guò)渡區(qū), 亦為獨(dú)立的巖溶地下水系統(tǒng)(圖 1-R2), 僅受大氣降水的補(bǔ)給。其流域總面積約28.90 km2, 全長(zhǎng)約7 km, 規(guī)模與板寨地下河相似, 平面展布形態(tài)呈簡(jiǎn)單樹(shù)枝狀, 構(gòu)造呈南北向展布。地下河水力坡度較緩, 約為 13.3‰, 主體流向自北向南, 于大良—潭頭谷地邊緣排泄, 形成地下河出口, 在地下河出口處修建了自動(dòng)化監(jiān)測(cè)站,獲取了自2015年8月至2018年12月共計(jì)3年多的連續(xù)數(shù)據(jù)。柳州農(nóng)耕區(qū)其他平原地區(qū)巖溶發(fā)育整體較弱, 地下水一般以分散排泄為主。

圖1 研究區(qū)地理位置圖Fig. 1 Geographical location of the study areas

2 數(shù)據(jù)采集

2010年, 在茂蘭國(guó)家級(jí)自然保護(hù)區(qū)布置了7個(gè)水化學(xué)取樣點(diǎn), 取樣頻次每次/月。于 2015年 9月至2016年4月, 在柳州市周邊SS1區(qū)、SS2和SS3等3個(gè)區(qū)域, 在豐水季節(jié)和枯水季節(jié)2批次共采集水樣39組。在板寨地下河出口和官村地下河出口處均建有監(jiān)測(cè)站, 安裝澳大利亞 Greenspan公司生產(chǎn)的 CDTP300多參數(shù)水質(zhì)監(jiān)測(cè)儀, 配備雨量筒和矩形薄壁堰。同步監(jiān)測(cè)指標(biāo)為水體水溫、水位、電導(dǎo)率、pH值以及降水量, 其精度分別為 0.1 ℃、0.01 m、0.01 μs/cm、0.1 pH單位和0.5 mm。碳同位素標(biāo)準(zhǔn)樣品為V-PDB, 測(cè)試儀器為MAT 235氣體穩(wěn)定同位素質(zhì)譜儀, 測(cè)定精度分別為δ13C≤0.01‰。

3 結(jié)果及討論

3.1 森林植被的氣候效應(yīng)

1984年茂蘭喀斯特森林科學(xué)考察過(guò)程中, 發(fā)現(xiàn)森林小氣候特征較明顯, 主要表現(xiàn)為云霧較多, 陰雨天氣出現(xiàn)頻繁, 日照時(shí)數(shù)相對(duì)較少, 且年均氣溫低于荔波縣氣象站 3℃左右, 降水量明顯較大(周政賢, 1987)。2009年, 中國(guó)地質(zhì)科學(xué)院巖溶地質(zhì)研究所在茂蘭保護(hù)區(qū)建立了板寨水文監(jiān)測(cè)站(海拔530 m), 對(duì)林區(qū)降水量進(jìn)行連續(xù)監(jiān)測(cè), 發(fā)現(xiàn)茂蘭保護(hù)區(qū)降水量值通常明顯大于周?chē)鷧^(qū)域。如2010年,茂蘭自然保護(hù)區(qū)降水量達(dá)1 710.0 mm, 而荔波縣氣象站(海拔432.9 m)降水量?jī)H為1 058.9 mm, 說(shuō)明茂蘭自然保護(hù)區(qū)森林的氣象水文效應(yīng)明顯, 森林激發(fā)降水量值占全年降水量的 38.08%, 顯著高于預(yù)測(cè)值(表 1)。

表1 茂蘭保護(hù)區(qū)與荔波氣象站降水情況對(duì)比表(單位: mm)Table 1 Comparison of precipitation data in Maolan karst forest and Libo station (unit: mm)

3.2 植被對(duì)流域水循環(huán)的調(diào)控作用

通過(guò)對(duì)多場(chǎng)暴雨過(guò)程中降水-徑流曲線分析(圖2), 發(fā)現(xiàn)板寨地下河和官村地下河的在雨季暴漲暴落現(xiàn)象非常明顯, 但去噪后退水過(guò)程流量曲線衰減部分?jǐn)M合效果基本一致, 說(shuō)明每場(chǎng)暴雨事件的流量衰減方式相同。從降水-徑流曲線上可以看出, 雖然板寨地下河流量的峰值較大, 但官村地下河對(duì)降水的響應(yīng)比板寨地下河迅速, 說(shuō)明森林植被能起到有效的截流消洪作用。

圖2 板寨(A)和官村(B)地下河降水-流量過(guò)程曲線Fig. 2 Rainfall-flow dynamic of Banzhai (A) and Guancun(B) subterranean river system

從地下河流量的衰減曲線來(lái)看, 初期地下河水流量衰減速率較快, 水文地質(zhì)固有條件控制洪水的消減過(guò)程。但在退水后期, 板寨地下河洪水過(guò)后退水曲線的衰減過(guò)程非常平滑, 說(shuō)明從管道流向基巖裂隙流的排水過(guò)渡過(guò)程中, 森林植被發(fā)揮了重要的調(diào)節(jié)作用; 而官村地下河流域缺乏森林植被覆蓋,其退水曲線的拐點(diǎn)十分明顯, 表現(xiàn)了由小型管道流直接向基巖裂隙流的過(guò)渡。

3.3 植被對(duì)巖溶作用過(guò)程的調(diào)控作用

研究表明, 海相碳酸鹽巖礦物δ13C值介于0‰至2.0‰之間(Drever, 1997); C3植物分布區(qū)土壤CO2的δ13C 值介于–23‰ ～ –40‰之間, 平均值為–27‰;C4植物區(qū)土壤CO2的δ13C值介于–9‰ ～ –19‰之間,平均值為–12‰(Smith and Epstein, 1971)。根據(jù)巖溶動(dòng)力學(xué)方程式, 巖溶水中碳有兩個(gè)來(lái)源, 其一為碳酸鹽巖礦物, 其二來(lái)自大氣(包括自由空氣和土壤空氣), 且大氣碳源和碳酸鹽礦物的碳源為1:1(B?gli, 1980; Dreybrodt, 1980; Yuan, 1997a)。因此有研究者認(rèn)為在良好的植被覆蓋條件下, 海相碳酸鹽巖分布區(qū)巖溶水中δ13CDIC的預(yù)測(cè)值應(yīng)該為–11‰(Grootes et al., 1969; Mook, 1986)。在柳州市農(nóng)耕區(qū), 受酸雨和施肥等影響, 大部分巖溶水點(diǎn)的值δ13CDIC介于–8‰ ～ –10‰之間(表 2), 統(tǒng)計(jì)值為(–9.08±4.70)‰(圖 3), 值域明顯較寬。而在森林覆蓋條件良好的茂蘭森林保護(hù)區(qū), 巖溶水中δ13CDIC值均小于–11‰(康志強(qiáng)和何師意, 2011), 統(tǒng)計(jì)值為(–14.18±1.64)‰, 其值整體較茂蘭自然保護(hù)區(qū)偏負(fù)。通過(guò)對(duì)比發(fā)現(xiàn), 森林區(qū)參與巖溶作用的大氣CO2的碳同位素受森林植被的控制左右較為明顯。

表2 柳州農(nóng)耕區(qū)分散流系統(tǒng)地下水中HCO濃度和δ13CDIC值Table 2 HCO concentrations and δ13CDIC values of groundwater in the Liuzhou decentralized flow system

表2 柳州農(nóng)耕區(qū)分散流系統(tǒng)地下水中HCO濃度和δ13CDIC值Table 2 HCO concentrations and δ13CDIC values of groundwater in the Liuzhou decentralized flow system

樣品序號(hào)HCO– 3/(mg/L)δ13CDIC /‰樣品序號(hào)HCO– 3/(mg/L)δ13CDIC/‰WS01 246 –10.50 WS21 227 –8.60 WS02 248 –8.96 WS22 133 –10.45 WS03 270 –4.86 WS23 232 –9.50 WS04 161 –7.22 WS24 132 –12.40 WS05 141 –10.56 WS25 234 –8.49 WS06 166 –9.10 WS26 235 –12.62 WS07 171 –6.04 WS27 240 –9.82 WS08 172 –7.74 WS28 246 –11.70 WS09 224 –10.80 WS29 191 –10.83 WS10 236 –9.18 WS30 248 –6.78 WS11 191 –7.71 WS31 283 –9.51 WS12 284 –10.02 WS32 287 –9.24 WS13 191 –6.12 WS33 183 –11.66 WS14 168 –10.46 WS34 301 –8.58 WS15 219 –9.03 WS35 271 –10.67 WS16 229 –8.44 WS36 176 –10.86 WS17 222 –8.28 WS37 210 –10.55 WS18 225 –9.66 WS38 261 –9.95 WS19 245 –9.70 WS39 235 –9.24 WS20 149 –9.54

圖3 研究區(qū)δ13CDIC值與水體HCO含量關(guān)系圖Fig. 3 Relationship between the δ13CDIC and HCO in the karst water

4 結(jié)論

不同類(lèi)型巖溶關(guān)鍵帶中地下水循環(huán)模式不同,其中的森林植被對(duì)地下水系統(tǒng)內(nèi)部水循環(huán)和水化學(xué)組分均有明顯的影響, 具體表現(xiàn)在以下方面:

(1)植被通過(guò)蒸騰作用而在林區(qū)激發(fā)二次降水,使林區(qū)降水量明顯增大;

(2)在降水-徑流曲線上, 森林覆蓋區(qū)降水洪峰產(chǎn)生比較滯后, 而地下水排泄過(guò)程中退水曲線較為平滑, 說(shuō)明森林植被具有明顯的截流消洪作用;

(3)相比柳州農(nóng)耕區(qū), 茂蘭自然保護(hù)區(qū)地下水中HCO–3濃度整體上較大, 碳穩(wěn)定同位素也明顯偏負(fù),其兩者的值域也明顯較小。說(shuō)明森林植被對(duì)巖溶關(guān)鍵帶地下水中無(wú)機(jī)碳的貢獻(xiàn)比較大, 并能有效增強(qiáng)巖溶作用。

Acknowledgements:

This study was supported by Scientific Research and Technology Development Program of Guangxi(Nos. GuikeAB110004 and GuikeAB18126062), National Key Research and Development Program of China (Nos. 2016YFC0502302 and 2016YFC0502500), Science and Technology Base and Talent Special Project of Guangxi (No. AD17129052),and National Natural Science Foundation of China(Nos. 41502260 and 41472239).