劉佳駒,趙雨順,黃 香,郭懷成*

?

雅魯藏布江流域水化學(xué)時(shí)空變化及其控制因素

劉佳駒1,趙雨順2,黃 香3,郭懷成1*

(1.北京大學(xué)環(huán)境科學(xué)與工程學(xué)院,北京 100871；2.天津大學(xué)環(huán)境科學(xué)與工程學(xué)院,天津 300072；3.西藏大學(xué)化學(xué)與環(huán)境學(xué)院,西藏 拉薩 850000)

為研究雅魯藏布江全流域水化學(xué)特征及其控制因素,2016年先后采集雅魯藏布江干流及其支流豐、平、枯三個(gè)水期水樣212個(gè),綜合運(yùn)用數(shù)理統(tǒng)計(jì)、Piper三線圖、Gibbs 模型和離子比等方法,分析了全流域河水的水文地球化學(xué)特征,并探討了流域的水化學(xué)演化規(guī)律.結(jié)果表明,河水中陽離子均以 Ca2+、Na+和Mg2+為主;陰離子以HCO3-和 SO42-為主;TDS平均值為204.51mg/L,礦化度較低;水化學(xué)類型以 HCO3·SO4(SO4·HCO3)-Ca·Mg(Mg·Ca)型水為主;雅魯藏布江主干河流內(nèi)的主要離子含量變化存在波動(dòng)且季節(jié)性變化顯著.空間上,河水中絕大多數(shù)離子遵循著先增大后減少的趨勢;水化學(xué)樣品均分布在 Gibbs 模型左中部,說明該流域水化學(xué)離子組成受巖石風(fēng)化作用控制;主成分分析及相關(guān)分析表明雅魯藏布江流域水化學(xué)組分受人為影響作用控制,方解石、白云巖的風(fēng)化及硫酸溶解也起到非常重要的作用.流域三個(gè)水期的重金屬元素的含量絕大多數(shù)能滿足地表水I類水體的要求.

水化學(xué)；雅魯藏布江；時(shí)空變化；水質(zhì)

青藏高原被稱為亞洲水塔.其為世界近1/3人口的生活、生產(chǎn)提供淡水資源[1].由于地形、地質(zhì)和氣候等自然條件以及人類生產(chǎn)、生活等的影響,水質(zhì)受多因素、多層次和多尺度相互綜合影響,水質(zhì)的變化會(huì)對(duì)水生生態(tài)系統(tǒng)產(chǎn)生一系列的生態(tài)效應(yīng),能在很大程度上反映流域的變化特征[2].作為世界的“第三極”,青藏高原在通過陸地和海洋能量交換發(fā)生的元素地球化學(xué)循環(huán)中占有非常重要的比例[1].研究河流和湖泊水循環(huán)和水化學(xué)特征以及水資源和環(huán)境的保護(hù)具有重要的科學(xué)意義.此外,在水化學(xué)研究中揭示最主要離子的來源以及其控制因素非常重要[3].

近年來,眾多國內(nèi)外學(xué)者開始對(duì)雅魯藏布江流域展開了較為深入的研究并取得了豐碩的成果,主要分自然影響與人為影響兩方面:隨著全球氣候變暖, 青藏高原的冰川面積減少,導(dǎo)致的水文水生態(tài)環(huán)境也在發(fā)生著變化,進(jìn)而影響了高原河流的水化學(xué)性質(zhì)[4-5].Adachi等[6]對(duì)高原東部雅魯藏布江流域主要化學(xué)組成的研究表明,巖石風(fēng)化風(fēng)化和侵蝕過程對(duì)流域水化學(xué)性質(zhì)具有強(qiáng)烈的影響;沈大軍等[7]對(duì)雅魯藏布江流域水資源進(jìn)行了評(píng)價(jià),認(rèn)為雅魯藏布江外流水系河水泥沙含量較小,河水礦化度較小,水質(zhì)良好;李紅敬等[8]通過對(duì)雅魯藏布江的水質(zhì)分析發(fā)現(xiàn),在6~10月,溶氧值低、水溫較高,12~2月,溶氧值高、水溫較低、8月懸浮物含量高,下游較中游的pH值和懸浮物含量低;Huang等[9]對(duì)雅魯藏布江全流域的水質(zhì)進(jìn)行了評(píng)估,并基于統(tǒng)計(jì)學(xué)對(duì)水化學(xué)的控制因素進(jìn)行了分析.一些學(xué)者的研究表明,雅魯藏布江的水化學(xué)特征變化一定程度上受到人類活動(dòng)的影響[10]:由于青藏高原是西風(fēng)、印度洋季風(fēng)和東亞季風(fēng)的交匯點(diǎn),大氣環(huán)流也導(dǎo)致了印度和尼泊爾等南亞國家的工農(nóng)業(yè)污染物跨界傳輸.高原上的人為活動(dòng)和氣候變化都將影響雅魯藏布江的水化學(xué)過程,并對(duì)青藏高原南部和印度和孟加拉國下游地區(qū)的人類生產(chǎn)生活構(gòu)成風(fēng)險(xiǎn)[8].Huang等[11]對(duì)雅魯藏布江流域水質(zhì)研究表明, Cd,Hg,Ni和Cu等重金屬含量相對(duì)較高,污染被認(rèn)為是由當(dāng)?shù)夭傻V業(yè)造成的.在鹽湖中,還發(fā)現(xiàn)了高濃度的重金屬元素,如鉛和砷.這可能是由于地?zé)岬茸匀贿^程引發(fā)[10].總的來說,雅魯藏布江水化學(xué)特性還沒有經(jīng)過系統(tǒng)深入研究,對(duì)全流域、季節(jié)性變化的研究較為缺乏.相對(duì)于我國其它流域,有關(guān)雅魯藏布江流域水化學(xué)特征、化學(xué)風(fēng)化過程,尤其是季節(jié)性變化及影響因素的研究尚不夠全面、細(xì)致、深入.因此,有必要繼續(xù)開展和深入研究雅魯藏布江全流域化學(xué)風(fēng)化和元素地球化學(xué)循環(huán)過程,為將來探索“第三極”以及全球氣候環(huán)境變化提供基礎(chǔ).

本文基于2016年對(duì)我國西藏境內(nèi)雅魯藏布江全流域干支流三個(gè)水文期的河水采樣,結(jié)合已有的流域環(huán)境背景資料, 分析報(bào)道了水體樣品中主要陰陽離子的含量及其時(shí)空分布特征,探討了水化學(xué)的主控因素,進(jìn)一步研究雅江水質(zhì)和“亞洲水塔”的水質(zhì)變遷提供基礎(chǔ)數(shù)據(jù),為了解流域水環(huán)境和水資源提供科學(xué)依據(jù).

1 研究區(qū)域概況及樣品采集與分析

1.1 研究區(qū)概況

圖1 青藏高原雅魯藏布江流域和采樣斷面分布

研究區(qū)地處西藏的南部,岡底斯山脈、念青唐古拉山脈以南,喜馬拉雅山脈以北的地區(qū).雅魯藏布江發(fā)源于西藏自治區(qū)西南部喜馬拉雅山北麓的杰馬央宗冰川,自西向東流經(jīng)仲巴、日喀則、澤當(dāng)?shù)鹊?在米林縣派區(qū)附近,形成一大拐彎,折向南流,經(jīng)中國巴昔卡出境,流入印度洋的孟加拉灣.流域面積為24萬km2,沿雅魯藏布江干流谷地,海拔高度為2800~4500m[12].

雅魯藏布江流域下游地區(qū)高溫多雨.降水的年際變化較小,年內(nèi)分配很不均勻.年均降水量從河流上游200mm左右至最下游的5000mm左右.全年平均降水主要集中在7~9月,達(dá)到全年降水量的50%~80%.年均最高溫度出現(xiàn)在6月,下游地區(qū)多出現(xiàn)在7月,年均溫度最低出現(xiàn)在1月.降水最多的月份,其冰雪融水補(bǔ)給河流的水量也大.此外,該流域還具有枯水期水量較大、懸移質(zhì)泥沙含量少、河水溫度低、河水礦化度小、總硬度低等特點(diǎn)[13].

雅魯藏布江有5條流域面積超過10000km2的支流,如,帕隆藏布、年楚河、尼洋河、拉薩河、多雄藏布.尼洋河是雅魯藏布江的一級(jí)支流,尼洋河所處地區(qū)的氣候溫和濕潤.年均氣溫8℃左右,年降水量 600~900mm,年內(nèi)徑流量相對(duì)集中,主要集中在每年的6~9月.拉薩河位于西藏中南部全長568km,流域面積31760km2.拉薩河流域面積僅占西藏自治區(qū)總面積的2.7%,而流域內(nèi)的人口、耕地卻約占全自治區(qū)的15%,是西藏工、農(nóng)、牧業(yè)集中的地區(qū).

1.2 水樣采集與分析

依據(jù)《水質(zhì)采樣技術(shù)指導(dǎo)》(HJ 494-2009)相關(guān)內(nèi)容以及雅魯藏布江流域地貌特征,并兼顧各支流匯入位置,在干流及其支流布設(shè)67個(gè)采樣點(diǎn),保證每個(gè)樣品點(diǎn)具有一定的代表性[11].采樣分別在2016年9月15日~10月8日(豐水期),2016年11月8~11月26日(平水期)和2017年4月22日~5月7日(枯水期)分3次采集樣品.如圖1所示.每個(gè)時(shí)間段樣本在半個(gè)月內(nèi)從下游到上游連續(xù)采集.為保持采樣條件的一致性,選擇24h內(nèi)沒有降雨的天數(shù)作為采集樣本日[14].

樣品收集之前,每個(gè)樣品瓶用硝酸洗滌,之后用蒸餾水沖洗.將水樣品收集在2L聚丙烯瓶中,然后立即通過0.45mm醋酸纖維素膜過濾器.將樣品置于干凈的自密封聚乙烯(PE)塑料袋中,并保存在4℃的冰箱中直至測試.為了確保數(shù)據(jù)的準(zhǔn)確性,所有樣品瓶在收集水樣之前用收集的濾液洗滌2~3次.為避免污染和破壞,將樣品置于冰箱中冷凍運(yùn)輸.將由水,過濾器和沉淀物組成的樣品運(yùn)送到實(shí)驗(yàn)室.將水樣在4℃的冰箱中儲(chǔ)存,并在自然風(fēng)化條件下將沉淀樣品保存在無塵實(shí)驗(yàn)室中[11].

雅魯藏布江河水基礎(chǔ)理化性質(zhì),陰陽離子及重金屬的測量儀器及方法如表1所示.在實(shí)驗(yàn)過程中,儀器由國家黃土標(biāo)準(zhǔn)物質(zhì)(GBW07408)標(biāo)定,并由標(biāo)準(zhǔn)溶液控制.陽離子和陰離子的分析精度優(yōu)于2%.總陽離子和陰離子的總和表明與使用探針測量的TDS值(Tz++Tz-=1.02TDS,2=0.98)匹配度良好,這表明數(shù)據(jù)質(zhì)量很高[11].

所有統(tǒng)計(jì)計(jì)算均使用Excel 2007(Microsoft Office)和spss18.0進(jìn)行.實(shí)驗(yàn)結(jié)果圖使用Origin8.0制作[11].首先,通過統(tǒng)計(jì)檢驗(yàn)確定離子元素的含量是否遵循正態(tài)分布.然后,使用Pearson相關(guān)分析評(píng)估離子元素濃度的相關(guān)關(guān)系.

表1 基礎(chǔ)理化性質(zhì)、陰陽離子及重金屬測量儀器及方法

2 結(jié)果和討論

2.1 河水水化學(xué)總體特征

2.1.1 基礎(chǔ)理化性質(zhì) 通過對(duì)常規(guī)水化學(xué)參數(shù)進(jìn)行分析,雅魯藏布江流域水樣品pH值、T、EC、DO、TDS變化范圍及均值如表2所示.雅魯藏布江流域三期水體的水溫為0.0~20.5℃,平水期平均溫度最低(4.2℃);雅魯藏布江及其支流的pH值變化幅度較小,介于7.97~9.03之間,平均值為8.46,整體呈弱堿性,碳酸鹽系統(tǒng)是天然水中優(yōu)良的緩沖系統(tǒng),它對(duì)避免天然水 pH 值急劇變化起緩沖作用.豐水季節(jié)河水 pH 值(8.6)略低于枯水期(8.7)及平水期(8.8),這可能的原因是豐水期降雨導(dǎo)致水體的pH值減小;3個(gè)水期TDS變化很大,13.30~451.00mg/L,平均值為183.95mg/L,高于世界河流的平均值(120mg/L)[15](表2);電導(dǎo)率(EC)可以反映水體的離子強(qiáng)度,雅魯藏布江流域三期水體的電導(dǎo)率變化范圍均值為18.5~632.0μS/cm,均值為252.3μS/cm;在溫度和鹽度的影響下,雅魯藏布江及其支流的DO在枯水期最高(均值為10.85mg/L),在豐水期最低(均值為8.1mg/L)(表2),這可能是水體中DO的含量與空氣中的氧的分壓,水溫有相關(guān),枯水期的溫度水體較低,因此導(dǎo)致DO含量高.雅魯藏布江流域主要河道的陽離子濃度依次為Ca2+>Na+>Mg2+> K+.Ca2+為是優(yōu)勢陽離子,占陽離子總量的31.3%.陰離子濃度順序?yàn)?HCO3->SO42->Cl-.其中HCO3-和SO42-為主要陰離子,分別占陰離子總量的72.35%和14.32%.Ca2+和SO42-的平均濃度分別是全球河流平均濃度的3~4以上[15](表3).

表2 基礎(chǔ)理化性質(zhì)

表3 雅魯藏布江及世界大河主要離子含量統(tǒng)計(jì)(mg/L)

2.1.2 水化學(xué)類型 河流水化學(xué)特征主要受天然補(bǔ)給來源、化學(xué)風(fēng)化以及人類活動(dòng)的影響.決定河流水化學(xué)特征的主要過程是降水對(duì)地表組成物質(zhì)的淋溶過程、主支流交匯處的混合過程、懸浮物的吸附和解析過程以及蒸發(fā)濃縮而導(dǎo)致的變質(zhì)等過程.Piper三線圖可以用來分析河水水化學(xué)成分的演化規(guī)律,結(jié)果如圖2所示.雅魯藏布江及其支流的樣品離子含量分布接近CO32-+HCO3-軸分布(圖2),陰離子主要由HCO3-和SO42-組成.在陽離子三角圖中,采樣點(diǎn)位于Ca2++Mg2+軸附近,表明該區(qū)域的陽離子主要為Ca2+和Mg2+,這可能反映了流域水化學(xué)性質(zhì)主要受碳酸鹽風(fēng)化.流域以SO4·HCO3(HCO3·SO4)- Ca·Mg(Mg·Ca)型為主.在雅魯藏布江及其支流中,80%以上的Ca2+, Mg2+和HCO3-產(chǎn)生于碳酸鹽風(fēng)化過程中.在河流的中上部,尤其是年楚河附近,Na+的比例高達(dá)29%,Cl-高達(dá)31%.說明蒸發(fā)巖的溶解是其主要控制因素.

圖2 雅魯藏布江干流及其支流河水陰陽離子三角圖

圖3 雅魯藏布江流域豐、平、枯3季水體重金屬濃度箱線圖

2.1.3 重金屬分析 青藏高原近30a的金屬采礦業(yè)影響了該地區(qū)的地表水質(zhì)量.重金屬含量的增加會(huì)對(duì)漁業(yè)和飲用水安全產(chǎn)生影響.如圖3所示雅魯藏布江流域在豐、平、枯3個(gè)水期的重金屬濃度的箱線圖.研究表明,Al(均值為62.65mg/L)和Fe(均值為60.8mg/L)濃度較高.Cu,Zn,Ag,Cd和Cr的濃度普遍較低,符合中華人民共和國地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB 3838-2002)的I類水標(biāo)準(zhǔn). 雅魯藏布江流域三個(gè)水期的As元素的含量大部分能滿足地表水I類水體的要求.但是,羊巴井地?zé)崽锵掠螛悠吩谟昙竞涂菟诘腁s濃度分別為170.5ug/L和144.0ug/L,超過了地表水環(huán)境類水V標(biāo)準(zhǔn).青藏高原沉積物中As的濃度是中國地殼豐度的11.6倍.地?zé)峄顒?dòng)釋放的氣液導(dǎo)致As含量高.其含量的超標(biāo)可能會(huì)影響該地區(qū)居民的飲水安全[18].90.2%水樣中汞的濃度達(dá)到地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)的I類水標(biāo)準(zhǔn).上游瑪玖藏布采樣點(diǎn)As濃度極高(枯水期為181.2ug/L,豐水期為217.58ug/L,平水期為182.1ug/L),這可能與采礦業(yè)以及As金屬的熔煉有關(guān)[19].但雅魯藏布江干流及其支流重金屬濃度季節(jié)變化不明顯.

2.2 河水主要離子時(shí)空變化特征

圖4 雅魯藏布江干流豐,平,枯三季主要離子濃度的變化

徑流補(bǔ)給來源及氣候的顯著差異是影響季風(fēng)區(qū)河流水文、水化學(xué)季節(jié)性特征的重要因素.通過分析雅魯藏布江干流30個(gè)采樣點(diǎn)的離子濃度,結(jié)果表明,該河流的水化學(xué)特征具有顯著的季節(jié)差異,如圖4所示.雅魯藏布江干流陽、陰離子濃度變化趨勢總體表現(xiàn)為豐水期最低.豐水季節(jié)流域河水中大部分陰陽離子平均濃度(mg/L)均小于枯水期及平水期,離子濃度值下降明顯,較枯水期離子濃度減少約 20.97mg/L,較平水期減少約為50.37mg/L.(表3)陽離子變化以Ca2+為主,枯水期較平水期減小了約9.35mg/L,其余陽離子變化不明顯.流域陰離子濃度總量變化幅度遠(yuǎn)遠(yuǎn)超過了陽離子,HCO3-變化最為明顯,豐水期較平水期減少了約37.42mg/L.豐水期降水量明顯增加,而離子來源沒有明顯變化的情況下,形成的稀釋效應(yīng).相反,平水期和枯水期的流速相對(duì)較慢,水體與巖石接觸時(shí)間較長,增加了巖石的溶解.在枯水期和豐水期,雅魯藏布江干流的陽離子濃度依次為,Ca2+>Na+>Mg2 +>K+,但在平水期,表現(xiàn)為Ca2+>Mg2+>Na+>K+,這種變化在中游和下游最為明顯(圖4).在3個(gè)水期,水樣中陰離子濃度變化趨勢均為HCO3->SO42->Cl->NO3-.

從空間上看,雅魯藏布江干流離子濃度從上游至下游含量變化存在波動(dòng),但總體呈先增后減的趨勢(圖4).最上游采樣點(diǎn)(Y17、Y16)測得的TDS濃度與整個(gè)干流相比明顯較低.這一結(jié)果可能是由于冰川和融雪稀釋了離子濃度.雅江上游地區(qū)表現(xiàn)為較為干燥的草原氣候,其蒸發(fā)量明顯超過降雨量,使得水體在一定程度上得到蒸發(fā)濃縮,大量裸地、灘涂的存在使得地表大面積暴露,導(dǎo)致當(dāng)?shù)仫L(fēng)化侵蝕作用強(qiáng)烈,大量鹽離子得以釋放進(jìn)入水體.因此,Y14和Y15采樣點(diǎn)離子濃度相比前段劇烈增加,增速最快的離子為Ca2+、Na+、HCO3-、SO42-.從源頭向下游地區(qū)過渡的過程中隨著降水增加,徑流量得以補(bǔ)充,同時(shí),灌木、針葉林、闊葉林的組合植被類型減輕了水土流失,最終導(dǎo)致了水體中的離子濃度呈現(xiàn)逐漸下降的趨勢.對(duì)比枯水季節(jié)與豐水季節(jié)流域離子空間特征分布發(fā)現(xiàn)(圖4),豐水季節(jié)陽離子總濃度值均低于枯水季節(jié).其中,上游降幅最大,約為 29.65mg/L,中游和下游降幅分別約為16.95mg/L和14.69mg/L,表明上游陽離子含量季節(jié)性差異最大,中游和下游依次降低.雅江流域從上游至下游離子濃度季節(jié)性差異呈逐步降低的趨勢.

2.3 離子來源分析

2.3.1 相關(guān)性分析 相關(guān)分析通常用于揭示主要離子的來源.表4顯示了雅魯藏布江流域各化學(xué)成分之間的相關(guān)矩陣.TDS與Ca2+,Mg2+,K+,Na+,HCO3-, Cl-和SO42-呈顯著正相關(guān)(<0.01),表明這些物質(zhì)都是TDS的來源.其中,TDS與Ca2+和HCO3-的相關(guān)性最顯著,相關(guān)系數(shù)分別為0.874和0.924.變化趨勢相同,統(tǒng)計(jì)學(xué)上表明Ca2+和HCO3-是TDS的主要來源.HCO3-與Ca2+和Mg2+有顯著相關(guān)性(<0.01).由于碳酸鹽和硫酸鹽礦物的風(fēng)化溶解作用,SO42-與 Ca2+、Mg2+的相關(guān)性也很強(qiáng)說明它們之間有可能來自同一來源,可能來源于碳酸鹽和硫酸鹽礦物的風(fēng)化溶解.Na+、K+都與Cl-有強(qiáng)相關(guān)性,說明它們有可能來自蒸發(fā)巖的溶解.

表4 離子相關(guān)性分析

2.3.2 主成分分析 為了進(jìn)一步探討雅魯藏布江流域水體主要離子與地質(zhì)背景之間的關(guān)系,了解主要離子的控制因素,應(yīng)用主成分分析法(PCA)對(duì)雅魯藏布江流域3個(gè)水期主要離子間的關(guān)系進(jìn)行分析.基于3個(gè)水期的212份水樣,通過因子分析得到3個(gè)主要公共因素.如圖5所示,3個(gè)水期的第一因子與K+,Na+,Cl-,Li,SiO2之間存在很強(qiáng)的相關(guān)性,這與高原各類蒸發(fā)巖的溶解有關(guān)[19].第二因子與Mg2+, SO42-,Ca2+和HCO3-存在較高相關(guān)性.說明方解石、白云巖的風(fēng)化及硫酸溶解的相關(guān)貢獻(xiàn)較大.其中碳酸鹽風(fēng)化作用是Ca2+的主要來源,80%~90%的溶解Ca2+和Mg2+是由雅魯藏布江流域碳酸鹽風(fēng)化產(chǎn)生的.第三個(gè)因素與Mn和Zn有較強(qiáng)相關(guān)性,表明與采礦業(yè)等人類活動(dòng)相關(guān),但人類活動(dòng)的影響并不明顯.其中流域豐水期及平水期的三個(gè)主要因素累計(jì)方差貢獻(xiàn)率為75.14%.第1,2,3因子貢獻(xiàn)率分別達(dá)到34.5%,26.2%和14.4%.枯水期的累計(jì)方差貢獻(xiàn)率為66.67%.對(duì)比平水期及豐水期可以發(fā)現(xiàn),枯水期累積方差百分比(66.67%)比平水期及豐水期(75.14%)略有降低,但離子負(fù)荷值增高,說明豐水期及平水期的主要離子來源的變化更加復(fù)雜(表5).

圖5 雅魯藏布江流域主要離子主成分分析

表5 雅魯藏布江流域3個(gè)水期離子和元素的主成分分析(PCA)

圖6 雅魯藏布江流域主要離子吉布斯圖

Y-雅魯藏布江干流,L-拉薩河,N-尼洋河,NY-年楚河,P-帕隆藏布江,D-多雄藏布江

2.3.3 吉布斯分析 Gibbs 設(shè)計(jì)了半對(duì)數(shù)坐標(biāo)圖,縱坐標(biāo)為水樣TDS的對(duì)數(shù)坐標(biāo),橫坐標(biāo)分別為陽離子質(zhì)量濃度的比值(Na+)/(Na++Ca2+)或者陰離子質(zhì)量濃度的比值(Cl-)/(Cl-+HCO3-)[21]. Gibbs圖顯示Na+/(Na++Ca2+)(<0.38)和Cl-/(Cl-+HCO3-)(<0.21)比值較低, TDS值中等,說明雅魯藏布江干流及其支流的離子化學(xué)性質(zhì)主要是由巖石風(fēng)化控制.所采集水樣遠(yuǎn)離“降水主導(dǎo)”和“蒸發(fā)結(jié)晶”作用控制區(qū)域,且河水的Cl-/Na+平均值為0.65,遠(yuǎn)低于世界平均海水比例(Cl-/Na+=1.15),這些均表明研究區(qū)內(nèi)大氣降水和蒸發(fā)結(jié)晶對(duì)水體中離子質(zhì)量濃度的影響不是主要控制因素.

3 結(jié)論

雅魯藏布江流域的地表水質(zhì)呈弱堿性.河流的平均TDS為204.51mg/L,高于全球平均水平.主要陽離子是Ca2+,Na+和Mg2+,而主要陰離子是HCO3-和SO42-.流域大部分地區(qū)的水化學(xué)特征為Ca2+-HCO3-型.大部分重金屬元素均能達(dá)到中華人民共和國地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB 3838-2002)的I類水標(biāo)準(zhǔn).然而,受到上游采礦及金屬冶煉的影響,有少數(shù)采樣點(diǎn)的As含量超標(biāo),中游羊八井地?zé)崽锔浇鼧悠稟s濃度較高,則與地?zé)峄顒?dòng)有關(guān).季節(jié)變化上,雅魯藏布江流域的水化學(xué)特征具有顯著的季節(jié)性差異,其中,豐水期的離子濃度最低,枯水期和平水期的離子濃度相對(duì)較高.在空間上,雅魯藏布江干流的離子濃度從上游到下游呈逐步降低的趨勢,支流中的離子濃度波動(dòng)顯著.主成分分析和吉布斯模型表明,雅魯藏布江流域的離子組成主要受巖石風(fēng)化作用控制.總體來說,雅魯藏布江絕大部分采樣斷面的水質(zhì)良好,極個(gè)別采樣斷面重金屬超標(biāo).在全球氣候變化的背景及人類活動(dòng)增多的雙重壓力下,雅魯藏布江的水化學(xué)特征會(huì)受到影響從而間接會(huì)對(duì)青藏高原的人類健康構(gòu)成風(fēng)險(xiǎn).本文的研究結(jié)論對(duì)進(jìn)一步研究雅江水質(zhì)和“亞洲水塔”的水質(zhì)變遷提供了基礎(chǔ),為了解流域水環(huán)境和水資源的科學(xué)管理提供基礎(chǔ)性科學(xué)依據(jù).

[1] Li F, Zhang Y, Xu Z, et al. The impact of climate change on runoff in the southeastern Tibetan Plateau [J]. Journal of Hydrology, 2013, 505(505):188-201.

[2] Cullaj A, Hasko A, Miho A, et al. The quality of Albanian natural waters and the human impact [J]. Environment International, 2005, 31(1):133-152.

[3] 樂嘉祥,王德春.中國河流水化學(xué)特征[J]. 地理學(xué)報(bào), 1963,(1):3-15.

[4] Wang S L, Jin H J, Li S X, et al. Permafrost degradation on the Qinghai-Tibet Plateau and its environmental impacts. [J]. Permafrost & Periglacial Processes, 2000,11(1):43–53.

[5] Wu W, Xu S, Yang J, et al. Silicate weathering and CO2, consumption deduced from the seven Chinese rivers originating in the Qinghai- Tibet Plateau [J]. Chemical Geology, 2008,249(3):307-320.

[6] Adachi K, Tainosho Y. Characterization of heavy metal particles embedded in tire dust. [J]. Environment International, 2004,30(8):1009.

[7] 沈大軍,陳傳友.青藏高原水資源及其開發(fā)利用[J]. 自然資源學(xué)報(bào), 1996,(1):8-14.

[8] 李紅敬,張 娜,林小濤.西藏雅魯藏布江水質(zhì)時(shí)空特征分析[J]. 河南師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010,38(2):126-130.

[9] X. Huang ?, Sillanp?? M, Gjessing E T, et al. Water quality in the southern Tibetan Plateau: chemical evaluation of the Yarlung Tsangpo (Brahmaputra) [J]. River Research & Applications, 2011,27(1):113–121.

[10] Yao T, Thompson L, Yang W, et al. Different glacier status with atmospheric circulations in Tibetan Plateau and surroundings [J]. Nature Climate Change, 2012,2(9):663-667.

[11] Huang X, Sillanp?? M, Gjessing E T, et al. Water quality in the Tibetan Plateau: major ions and trace elements in the headwaters of four major Asian rivers [J]. Science of the Total Environment, 2009,407(24):6242

[12] Jiang L, Yao Z, Wang R, et al. Hydrochemistry of the middle and upper reaches of the Yarlung Tsangpo River system: weathering processes and CO2consumption [J]. Environmental Earth Sciences, 2015,74(3):2369-2379.

[13] Huang Xiang, Sillanp??, M, Gjessing, E. T, et al. Environmental impact of mining activities on the surface water quality in Tibet: Gyama valley [J]. Science of the Total Environment, 2010,408(19):4177-4184.

[14] Pant R R, Zhang F, Rehman F U, et al. Spatiotemporal variations of hydrogeochemistry and its controlling factors in the Gandaki River Basin, Central Himalaya Nepal [J]. Science of the Total Environment, 2017,622-623:770.

[15] Gaillardet J, Dupré B, Louvat P, et al. Global silicate weathering and CO2, consumption rates deduced from the chemistry of large rivers [J]. Chemical Geology, 1999,159(1–4):3-30.

[16] Sarin M M, Krishnaswami S, Dilli K, et al. Major ion chemistry of the Ganga- Brahmaputra river system: Weathering processes and fluxes to the Bay of Bengal [J]. Geochimica Et Cosmochimica Acta, 1989,53(5):997-1009.

[17] Ahmad T, Khanna P P, Chakrapani G J, et al. Geochemical characteristics of water and sediment of the Indus river, Trans- Himalaya, India: constraints on weathering and erosion [J]. Journal of Asian Earth Sciences, 1998,16(2/3):333-346.

[18] Thomas J, Joseph S, Thrivikramji K P, et al. Seasonal variation in major ion chemistry of a tropical mountain river, the southern Western Ghats, Kerala, India [J]. Environmental Earth Sciences, 2014,71(5): 2333-2351.

[19] Dekov V M, Komy Z, Araújo F, et al. Chemical composition of sediments, suspended matter, river water and ground water of the Nile (Aswan-Sohag traverse). [J]. Science of the Total Environment, 1997,201(3):195.

[20] Jiang L, Yao Z, Wang R, et al. Hydrochemistry of the middle and upper reaches of the Yarlung Tsangpo River system: weathering processes and CO2consumption [J]. Environmental Earth Sciences, 2015,74(3):2369-2379.

[21] Gibbs R J. Mechanisms Controlling World Water Chemistry [J]. Science, 1971,172(3985):870-2.

[22] 劉 劍,姚治君,陳傳友.雅魯藏布江徑流變化趨勢及原因分析[J]. 自然資源學(xué)報(bào), 2007,(3):471-477.

[23] 楊志剛,卓 瑪,路紅亞,等.1961~2010年西藏雅魯藏布江流域降水量變化特征及其對(duì)徑流的影響分析[J]. 冰川凍土, 2014,36(1):166-172.

[24] 陳西平.雅魯藏布江水環(huán)境的探討[J]. 環(huán)境科學(xué), 1989,(5):85-89.

[25] Diamantini E, Lutz S R, Mallucci S, et al. Driver detection of water quality trends in three large European river basins [J]. Science of the Total Environment, 2018,612:49-62.

[26] Zayed J, Guessous A, Lambert J, et al. Estimation of annual Mn emissions from MMT source in the Canadian environment and the Mn pollution index in each province [J]. Science of the Total Environment, 2003,312(1-3):147.

[27] Jiang L, Yao Z, Wang R, et al. Hydrochemistry of the middle and upper reaches of the Yarlung Tsangpo River system: weathering processes and CO2consumption [J]. Environmental Earth Sciences, 2015,74(3):2369-2379.

致謝：本研究由西藏大學(xué)及西藏自治區(qū)各級(jí)政府支持下完成,在此表示感謝.

Spatiotemporal variations of hydrochemistry and its controlling factors in the Yarlung Tsangpo River.

LIU Jia-ju1, ZHAO Yu-shun2, HUANG Xiang3, GUO Huai-cheng1*

(1.College of Environmental Sciences and Engineering, Peking University, Beijing 100871, China；2.College of Environmental Sciences and Engineering, Tianjin University, Tianjin 300072, China；3.Department of Chemistry and Environmental Science, Tibet University, Lhasa 850000, China)., 2018,38(11)：4289~4297

In order to study the characteristics of the water chemistry of the Yarlung Tsangpo River and its control factors, 212 water samples from the mainstream of the Yarlung Tsangpo River and its tributaries were collected in the three water periods in 2016. The mathematical statistics, Piper three-line map, Gibbs model, and ion ratio were analyzed. The hydrochemical characteristics of rivers in the whole basin, and the evolution of water chemistry in the basin. The results showed that the cations in the river water were dominated by Ca2+, Na+, and Mg2+; the anions were dominated by HCO3-and SO42-; the average value of TDS was 204.51mg/L, the degree of salinity is relatively low; the type of water chemistry is dominated by HCO3·SO4(SO4·HCO3)-Ca·Mg (Mg·Ca) type water; The main ion concentration changes in the main stream of Yarlung Tsangpo fluctuated, and the seasonal changes were significant. Spatially, most ions in the river follow the trend of increasing first and then decreasing. The chemical samples of water are distributed in the middle left part of the Gibbs model, indicating that the chemical composition of the watershed in the basin is controlled by rock weathering. Principal component analysis and related analysis show that the Yarlung Tsangpo River Basin The chemical composition of water is controlled by human influence, and weathering of calcite, dolomite and dissolution of sulfuric acid also play a very important role. The vast majority of heavy metals in the three watersheds of the basin meet the requirements of Class I water bodies in surface waters.

hydrochemical；Yarlung Tsangpo River；spatio-temporal variations；river water quality

X826

A

1000-6923(2018)11-4289-09

劉佳駒(1991-),男,河北保定人,北京大學(xué)博士研究生,研究方向?yàn)榄h(huán)境規(guī)劃.發(fā)表論文3篇.

2018-04-02

西藏環(huán)境水化學(xué)調(diào)查及飲用水安全評(píng)價(jià)技術(shù)研究項(xiàng)目(22015FY111000)

* 責(zé)任作者, 教授, ghc1111@pku.edu.cn