晉建霞 張勝利 陸 斌 馬國棟

(西北農(nóng)林科技大學，楊凌，712100)

西安是世界著名的歷史文化名城，位于陜西關中平原的中部，是我國中西部地區(qū)重要的科研、高等教育、國防科技工業(yè)和高新技術產(chǎn)業(yè)基地，但人均占有地表水資源量遠遠低于國際公認的維持一個地區(qū)可持續(xù)發(fā)展人均占有量1 000 m3的臨界值［1］。隨著西安都市圈的加快建設，西安將建成為國際大都市［2］，其對水的需求將與日俱增，因此，優(yōu)質(zhì)的水源對人類健康、城市發(fā)展至關重要［3］。西安市城市用水主要從發(fā)源于太白山國家級自然保護區(qū)的黑河和石頭河取水，特別是近年來實施的石頭河、黑河引水工程。秦嶺北麓位于秦嶺分水嶺至關中平原南緣之間，具有重要的生態(tài)功能，是關中地區(qū)的生態(tài)屏障和水源涵養(yǎng)地，對西安甚至關中地區(qū)社會經(jīng)濟可持續(xù)發(fā)展影響深遠。但多年來重開發(fā)輕保護，尤其是多年來礦產(chǎn)開發(fā)、旅游開發(fā)、污染型“五小企業(yè)”生產(chǎn)［4］，使排污加劇，從而對水源地水環(huán)境安全構成極大威脅。西安市是陜西的政治、經(jīng)濟、文化中心，其供水水源地水質(zhì)安全是陜西省政府十分關注的重大問題。

森林生態(tài)系統(tǒng)水質(zhì)空間變化及其機理是森林水文學研究的重要內(nèi)容，目前在這方面的研究還較少。而太白山蒿坪保護區(qū)位于秦嶺北麓黑河和石頭河之間林區(qū)的中部，地處陜西省眉縣境內(nèi)，其氣候、土壤、植被、樹種等在水源地林區(qū)具有較好的代表性，因此，文中以該區(qū)天然林為對象，采用比較分析方法，系統(tǒng)探討雨水進去森林生態(tài)系統(tǒng)后水質(zhì)的空間變化及其機理，理清其影響水質(zhì)的關鍵過程，以便為西安市供水水源地保護與生態(tài)建設提供科學依據(jù)。

1 試驗地概況

研究區(qū)域位于陜西省寶雞市眉縣營頭鎮(zhèn)的太白山國家級自然保護區(qū)蒿坪林區(qū)(圖1)。地理坐標為東經(jīng) 107°41'～107°43'、北緯 33°04'～34°05'，海拔1 300～2 800 m，地處我國西北部暖溫帶，屬典型的暖溫帶半濕潤氣候，有著明顯的氣候垂直帶，海拔1 000～2 000 m區(qū)域的年均溫為11.4℃，海拔2 000 m以上年均溫低于8℃。東、北部的黑河和石頭河基本上流入渭河(圖2)。

區(qū)內(nèi)森林資源豐富，森林面積45 725 hm2，森林覆蓋率81.2%，總蓄積705萬m3。植被垂直帶譜發(fā)育十分典型，自下而上依次為落葉闊葉林帶、針葉林帶和高山灌叢草甸帶。主要樹種有華山松(PinusarmandiiFranch.)、栓 皮 櫟 (QuercusvariabilisBlume)、銳齒櫟(Quercusalienavar.acuteserrata)、油松(Pinustabuliformis)、青榨槭(Acerdavidiifranch)等。土壤類型眾多，有高山草甸土、山地暗棕壤、山地沼澤土等6個土類，16個亞類。

圖1 蒿坪林區(qū)位置圖

圖2 蒿坪林區(qū)水系布局圖(圖中黑線代表河流)

2 材料與方法

2.1 試驗設計與水樣采集

大氣降水進入流域森林生態(tài)系統(tǒng)后，依次通過林冠層、枯枝落葉層到達土壤層表面，并分別形成林內(nèi)雨和枯透水。在形成林內(nèi)雨過程中，還伴有樹干莖流的產(chǎn)生，但其對水質(zhì)總體影響較小［5］，本研究不將其單獨列出。到達土壤層表面的降水，幾乎全部下滲補充土壤水分或形成壤中流，這是因為森林土壤表層富含有機質(zhì)，質(zhì)地疏松，空隙大，利于入滲［6］。

本試驗采用對比大氣降水，林內(nèi)雨，枯透水和土壤 A(0≤H(土壤深)＜20 cm)、B(20 cm≤H＜50 cm)層滲透水，溪流水以及流域出口徑流水的化學成分質(zhì)量濃度變化，來探討水質(zhì)的空間變化。因此，需要采集的水樣類型為大氣降水，林內(nèi)雨，枯透水，土壤A、B層滲透水，溪流水和流域出口徑流水。降水的收集點設在蒿坪保護站的空曠地帶，共布置2個點，附近無樹木、房屋等。分別在針闊混交林、闊葉混交林和銳齒櫟林內(nèi)布設林內(nèi)雨收集器，每個林子布設2個點，各點代表的面積基本相同，分別用口徑20 cm、高15 cm的塑料桶接取，其頂部用帶孔的聚氯乙烯薄膜罩住，以防雜物進入;枯透水分別用口徑25 cm、高10 cm的塑料盆接取，其頂部用帶孔的聚乙烯薄膜罩住，薄膜上放置與林內(nèi)厚度相似的枯落物，盡量不要破壞枯落物的結構，采樣時，將平行的水樣混合，裝入聚乙烯塑料瓶中，大概500 mL;土壤A、B層的水樣分別在集水槽的引水管管口接取，同樣裝入聚乙烯塑料瓶中，大概500 mL;另外，在典型支溝口修建復合式量水槽，并且在其三角出口接取溪流水500 mL;典型的流域出口修建巴歇爾量水槽，并且在其下游出口接取流域出口徑流水500 mL。溪流水和流域出口徑流水樣一般在雨后8～24 h采集。

2.2 水樣測試項目與方法

pH值采用電位法測定;NO-3和NH+4采用AA3型連續(xù)流動化學分析儀測定;采用原子吸收分光光度計火焰法測定 K+、Na+、Ca2+和 Mg2+，石墨爐法測定Cu2+和總Cr;采用等離子體發(fā)射光譜儀測定Cd2+、Zn2+、Pb2+和 Al3+。

2.3 數(shù)據(jù)處理

研究期間，共采集水樣10次，63個水樣，其中大氣降水6個、林內(nèi)雨10個、枯透水10個、土壤A層滲透水10個、土壤B層滲透水7個、溪流水10個、流域出口徑流水10個。大氣降水、林內(nèi)雨、枯透水、土壤A層滲透水、土壤B層滲透水、溪流水和流域出口徑流水中化學成分質(zhì)量濃度分別為各次降水同類水樣測定結果的平均值。

3 結果與分析

3.1 不同層次水樣pH值變化

在研究期間，大氣降水pH的平均值為7.29(表1)，水質(zhì)偏堿性。降水通過林冠層和枯枝落葉層后，pH值依次降低，其中，林內(nèi)雨的pH平均值降低幅度較大;經(jīng)過土壤層后，pH的平均值升高到7.45，并且土壤B層相對土壤A層的升高幅度最大，為0.54，這是由于H+與黏土礦物或有機質(zhì)吸附的陽離子發(fā)生交換反應，減少了H+的質(zhì)量濃度［7］。流域出口徑流水pH的平均值達到7.96，這是由降水與土壤基質(zhì)和溝道巖石接觸發(fā)生分解交換作用導致的［8］。

表1 不同層次水樣pH值及各化學成分質(zhì)量濃度

3.2 不同層次水樣中N、N質(zhì)量濃度變化

3.3 不同層次水樣中 K +、Na+、Ca2+、Mg2+質(zhì)量濃度變化

大氣降水K+的平均質(zhì)量濃度為0.525 mg·L-1(表1)。降水通過林冠層后，K+的質(zhì)量濃度增加了1.016 mg·L-1，且林內(nèi)雨中K+的質(zhì)量濃度為大氣降水中的 3 倍，Cappellato et al.［11］的研究也說明 K+在穿透雨增加很多，因為K+是移動性極強的元素，在水中極易溶解，極易從林冠的葉子和枝上溶脫［12］;經(jīng)過枯枝落葉層后，K+的質(zhì)量濃度增加了2.027 mg·L-1，薛立等［13］也發(fā)現(xiàn)枯落物為K+源，因此林冠層和枯枝落葉層都對K+表現(xiàn)為淋溶;通過土壤層后，K+質(zhì)量濃度升高到11.307 mg·L-1，說明土壤中含有豐富的K+;但是在溪流水和流域出口徑流水中的K+質(zhì)量濃度相對土壤層是非常低，流域出口徑流水K+質(zhì)量濃度為 0.996 mg·L-1，說明 K+在森林中循環(huán)時，森林生態(tài)系統(tǒng)起重要作用。袁東海等［14］研究表明植被覆蓋度低，無水保措施的試驗區(qū)K+流失量最大，而秦嶺北麓植被覆蓋度高(81.2%)，有效地防止K+流失，維持了K+的平衡。

大氣降水 Na+的平均質(zhì)量濃度為0.478 mg·L-1(表1)。當降水通過林冠層后，Na+質(zhì)量濃度降為最低值，為0.089 mg·L-1，說明林冠層對 Na+具有吸附作用;通過枯枝落葉層后，Na+質(zhì)量濃度為林內(nèi)雨的2倍多，說明枯落物對Na+表現(xiàn)為淋溶;經(jīng)土壤層后，Na+升高到最大值，5.589 mg·L-1，尤其是土壤B層(相對土壤A層增幅達4.316 mg·L-1)，說明土壤層對Na+有淋溶作用;但流域出口徑流水Na+的質(zhì)量濃度降為2.680 mg·L-1，這是由于溝道附近樹木少，降水直接進入流域出口徑流水中，對Na+離子產(chǎn)生了稀釋作用。

降水中Ca2+的平均質(zhì)量濃度為4.257 mg·L-1(表1)。Ca2+質(zhì)量濃度變化與Na+相似，也是林冠層吸附，與陳步峰等［12］得出的林冠層對降水中Ca2+具有吸貯效應相一致;枯枝落物層對Ca2+具有淋溶作用;從土壤層到溪流水后，Ca2+的質(zhì)量濃度升高到14.104 mg·L-1，這主要是由于土壤中的 CaCO3轉化為可溶于水的Ca(HCO3)2，從而使溪流水中Ca2+質(zhì)量濃度升高［15］。流域出口徑流水Ca2+的質(zhì)量濃度為10.585 mg·L-1，相對溪流水有所降低，這可能是因為徑流水的Ca2+多以CaCO3存在于巖石中。

大氣降水Mg2+的平均質(zhì)量濃度很低，為0.187 mg·L-1(表1)。降水通過林冠層后，Mg2+的質(zhì)量濃度增加了0.104 mg·L-1，經(jīng)枯枝落葉層后，Mg2+的質(zhì)量濃度為林內(nèi)雨的6倍多，說明林冠層和枯枝落葉層對Mg2+表現(xiàn)為淋溶，尤其是枯枝落葉層;通過土壤層后，Mg2+的變化不大，說明土壤層對Mg2+的影響不大，這與張勝利等［5］的研究一致;但是當流入支溝溪流后，Mg2+的質(zhì)量濃度升高到4.222 mg·L-1，這與 Mg2+主要來源于巖石［5］有關。相對大氣降水來說，Mg2+的質(zhì)量濃度是上升的，各層次對Mg2+均表現(xiàn)為淋溶。

3.4 不同層次水樣中 C u2+、總 C r、Al3+、Cd2+、Pb2+、Zn2+質(zhì)量濃度變化

大氣降水Cu2+的平均質(zhì)量濃度為3.709μg·L-1(表2)。當降水通過林冠層后，Cu2+的質(zhì)量濃度降低了18%，這與梁翠萍等［16］的研究一致;通過枯枝落葉層后，Cu2+的質(zhì)量濃度相對于林內(nèi)雨增加了1倍多，說明枯枝落葉能夠促進Cu2+向水中遷移;通過土壤層后，Cu2+的質(zhì)量濃度相對枯透水降低了25%，說明土壤對 Cu2+具有固定作用。溪流水中Cu2+的質(zhì)量濃度達到3.095μg·L-1;流域出口徑流水Cu2+的質(zhì)量濃度又增到5.572μg·L-1，可能是由于枯枝落葉釋放的Cu2+進入到流域出口徑流中。從大氣降水到流域出口徑流，Cu2+的質(zhì)量濃度是增加的，但各層次Cu2+質(zhì)量濃度均在我國生活飲用水水質(zhì)標準規(guī)定的1 000μg·L-1范圍內(nèi)。

大氣降水中總Cr的平均質(zhì)量濃度為0.360μg·L-1(表2)?？侰r質(zhì)量濃度在各層次的變化與Cu2+相似。降水通過林冠層后，總Cr的質(zhì)量濃度降低了42%;經(jīng)過枯枝落葉層后，總Cr的質(zhì)量濃度為林內(nèi)雨的4倍;通過土壤層后，相對枯透水下降了64%;溪流水中總Cr的質(zhì)量濃度達到0.257μg·L-1;流域出口徑流水總Cr的質(zhì)量濃度相對溪流水增加1倍，原因與Cu2+相似。在各層次中總Cr的質(zhì)量濃度都很低，符合我國生活飲用水水質(zhì)標準規(guī)定的50 μg·L-1。

大氣降水中Al3+的平均質(zhì)量濃度為92.181μg·L-1(表2)。大氣降水通過林冠層后，Al3+的質(zhì)量濃度增加了25.047μg·L-1，說明林冠層對 Al3+具有釋放作用;通過枯枝落葉層后，Al3+的質(zhì)量濃度相對林內(nèi)雨降低了46%，說明枯枝落葉層對Al3+有明顯的固定作用;經(jīng)過土壤層后，Al3+又升高到160.300 μg·L-1，并且土壤A層和B層的質(zhì)量濃度相近，說明土壤層對Al3+具有釋放作用;流域出口徑流水Al3+的質(zhì)量濃度相對土壤層降低了12%。歐陽學軍等［7］對鼎湖山地表水的研究發(fā)現(xiàn)地表水和30 cm深土壤溶液中的Al3+質(zhì)量濃度分別是大氣降水的5倍和8倍，得出地表水中的Al3+主要來源于酸雨對土壤的淋溶，而本研究中流域出口徑流水Al3+的質(zhì)量濃度是大氣降水的1.53倍，間接說明了太白山?jīng)]有受到酸雨的影響。

表2 不同層次水樣中各金屬元素質(zhì)量濃度 μg·L-1

大氣降水中Cd2+的平均質(zhì)量濃度很低，為1.864 μg·L-1(表 2)。降水通過林冠層后，Cd2+的質(zhì)量濃度降低1.4%，為1.838μg·L-1;經(jīng)枯枝落葉層后，Cd2+質(zhì)量濃度相對林內(nèi)雨降低12.2%，成為1.613μg·L-1，說明林冠層和枯枝落葉層對 Cd2+的作用不明顯，這與楊小梅［17］的研究一致;通過土壤層后，Cd2+質(zhì)量濃度增加到 2.281 μg·L-1，這與森林土壤枯枝落葉層和有機質(zhì)層的淋溶對土壤重金屬元素有很大的解吸作用有關［6］;在流域出口徑流水中Cd2+的質(zhì)量濃度卻降到 1.283 μg·L-1(pH=7.96)，這是因為當pH＞7.50時，Cd2+主要以氧化物與黏土礦物結合態(tài)及殘留態(tài)存在［18］。整體來說，各層次水樣的Cd2+都在我國生活飲用水水質(zhì)標準5μg·L-1范圍內(nèi)。

大氣降水中Pb2+的平均質(zhì)量濃度為26.612μg·L-1(表2)，超過了我國生活飲用水水質(zhì)標準規(guī)定的10μg·L-1。降水經(jīng)過林冠層后，Pb2+的質(zhì)量濃度降到6.476 μg·L-1，為大氣降水的 24.3%，因為植物可以通過枝葉吸收并累積Pb2+［19］;通過枯枝落葉層后，Pb2+又升高到 17.875 μg·L-1，說明枯落物對Pb2+具有淋溶作用;經(jīng)過土壤層后，Pb2+又降低到5.130 μg·L-1，為枯透水的 28.7%，歐陽學軍等［7］也得出森林的存在及其土壤對Pb2+有凈化作用，土壤對Pb2+的吸附強度與有機質(zhì)質(zhì)量濃度呈正相關［20］，而秦嶺山區(qū)森林植被下發(fā)育的土壤有機質(zhì)積累強烈［21］，所以，土壤層對重金屬的吸附能力較強;使得流域出口徑流水的Pb2+質(zhì)量濃度為5.259μg·L-1，符合我國生活飲用水標準。

Zn2+在各層次的質(zhì)量濃度變化趨勢與Pb2+的相同。通過林冠層后，Zn2+降低到0.803μg·L-1(表2)，為大氣降水的7.7%，幾乎全被林冠層吸收，與田大倫［22］、方晰［23］等研究發(fā)現(xiàn)杉木林林冠對 Zn2+有吸收作用相一致;通過枯枝落葉層后，Zn2+的質(zhì)量濃度增加到11.120μg·L-1，說明枯落物對 Zn2+具有釋放作用;經(jīng)過土壤層后，Zn2+的質(zhì)量濃度降到6.238 μg·L-1，說明土壤層對 Zn2+有固定作用;流域出口徑流水的Zn2+質(zhì)量濃度為5.657μg·L-1，各層次水樣的Zn2+質(zhì)量濃度都遠遠小于我國生活飲用水水質(zhì)標準1 000μg·L-1。

4 結束語

通過分析太白山蒿坪保護區(qū)森林生態(tài)系統(tǒng)對水質(zhì)空間變化的影響，得出西安市供水水源地森林生態(tài)系統(tǒng)水質(zhì)空間變化特征為：從大氣降水到流域出口徑流，pH值升高，貢獻較大者為土壤層。但是降水通過林冠層和枯枝落葉層后，pH值降低，N的質(zhì)量濃度也降低，枯枝落葉層是吸收N的關鍵層次，土壤層對N具有淋溶作用。林冠層和土壤層是吸收或固定N的關鍵層次，枯枝落葉層對N具有淋溶作用。相對大氣降水來說，流域出口徑流水的N質(zhì)量濃度上升，N質(zhì)量濃度降低。森林生態(tài)系統(tǒng)各層次對K+離子都表現(xiàn)為淋溶，尤其是林冠層和土壤層;林冠層是吸收Na+和Ca2+的關鍵層次，枯枝落葉層和土壤層則表現(xiàn)為淋溶;林冠層和枯枝落葉層對Mg2+具有淋溶作用，尤其是枯枝落葉層，土壤層對其影響不大。相對大氣降水來說，流域出口徑流 K+、Na+、Ca2+、Mg2+質(zhì)量濃度均表現(xiàn)為上升。大氣降水通過林冠層后，Cu2+、總 C r、Cd2+、Pb2+和Zn2+質(zhì)量濃度均降低，Al3+質(zhì)量濃度上升;經(jīng)過枯枝落葉層后，Cu2+、總 C r、Pb2+和 Z n2+質(zhì)量濃度均上升，Cd2+和Al3+質(zhì)量濃度下降;經(jīng)過土壤層后，Cu2+、總 C r、Pb2+和 Z n2+質(zhì)量濃度均下降，Cd2+和 A l3+質(zhì)量濃度上升，整體來說林冠層和土壤層對重金屬的吸附作用明顯，相對大氣降水來說，流域出口徑流水的Cu2+、總 Cr和 Al3+質(zhì)量濃度上升，Cd2+、Pb2+和 Zn2+質(zhì)量濃度下降，符合我國生活飲用水水質(zhì)標準。

［1］陳文郁.西安市城市飲用水源地安全狀況淺析［J］.水資源與水工程學報，2009，20(3)：152-154.

［2］王永勝，張定青.西安市秦嶺北麓村鎮(zhèn)生態(tài)化建設規(guī)劃初探：以周至縣為例［J］.華中建筑，2010(12)：126-130.

［3］仵彥卿.西安市渭濱傍河水源地地下水質(zhì)變化預測［J］.西安理工大學學報，2000，16(4)：348-352.

［4］肖玲，王書轉，張健，等.秦嶺北麓主要河流的水質(zhì)現(xiàn)狀調(diào)查與評價［J］.干旱區(qū)資源與環(huán)境，2008，22(1)：74-78.

［5］張勝利，李光錄.秦嶺火地塘森林生態(tài)系統(tǒng)不同層次的水質(zhì)效應［J］.生態(tài)學報，2007，27(5)：1838-1844.

［6］張勝利.秦嶺火地塘林區(qū)森林生態(tài)系統(tǒng)對降水中重金屬的作用［J］.林業(yè)科學，2009，45(11)：55-42.

［7］歐陽學軍，周國逸，黃忠良，等.鼎湖山森林地表水水質(zhì)狀況分析［J］.生態(tài)學報，2002，22(9)：1373-1379.

［8］張勝利.秦嶺火地塘森林水質(zhì)的季節(jié)性變化特征［J］.環(huán)境科學，2008，29(2)：316-321.

［9］王青山，何利平.土壤有機質(zhì)與氮素供應的相關關系［J］.山西林業(yè)科技，2003，9(增刊)：25-27.

［10］張勝利.南水北調(diào)中線工程水源林生態(tài)系統(tǒng)對水質(zhì)的影響：以秦嶺南坡中山地帶火地塘林區(qū)為例［J］.林業(yè)科學，2008，44(2)：15-22.

［11］Cappellato R，Peters N E.Dry deposition and canopy leaching rates in deciduous and coniferous forests of Georgia Piedmont：an assessment of regression model［J］.J Hydro，1995，169(1/4)：131-150.

［12］陳步峰，陳勇，尹光天，等.珠江三角洲城市森林植被生態(tài)系統(tǒng)水質(zhì)效應研究［J］.林業(yè)科學研究，2004，17(4)：453-460.

［13］薛立，傅靜丹，鄭衛(wèi)國，等.有無枯落物覆蓋對加勒比松林地表徑流及其K流失的影響［J］.林業(yè)科學研究，2010，23(4)：510-514.

［14］袁東海，王兆騫，郭新波，等.紅壤小流域不同利用方式土壤鉀素流失特征研究［J］.水土保持通報，2003，23(3)：16-20.

［15］石福臣，李鳳英，蔡體久，等.不同森林群落類型溪流水化學特征的季節(jié)動態(tài)［J］.應用生態(tài)學報，2008，19(4)：717-722.

［16］梁翠萍，張勝利.秦嶺銳齒櫟林和油松林林冠層對大氣降雨水質(zhì)的影響［J］.東北林業(yè)大學學報，2011，39(2)：53-56.

［17］楊小梅.縉云山森林生態(tài)系統(tǒng)水質(zhì)效應及評價研究［D］.北京：北京林業(yè)大學，2010.

［18］廖敏，黃常用，謝正苗，等.pH對鎘在土水系統(tǒng)中的遷移和形態(tài)的影響［J］.環(huán)境科學學報，1999，19(1)：81-86.

［19］薛皎亮，劉紅霞，謝映平.城市空氣中鉛在國槐樹體內(nèi)的積累［J］.中國環(huán)境科學，2000，20(6)：536-539.

［20］廖自基.環(huán)境中微量重金屬元素的污染危害與遷移轉化［M］.北京：科學出版社，1989：85.

［21］常慶瑞，雷梅，馮立孝，等.秦嶺北坡土壤發(fā)生特性與系統(tǒng)分類［J］.土壤學報，2002，39(2)：227-235.

［22］田大倫.會同廣坪林區(qū)降水和杉木林內(nèi)雨的養(yǎng)分含量［J］.中南林學院學報，2002，22(3)：9-13.

［23］方晰，田大倫，項文化，等.第2代杉木人工林微量元素的積累、分配及其生物循環(huán)特征［J］.生態(tài)學報，2003，23(7)：1313-1320.