任 巖，張 飛，王 娟，張 月，楊勝天, Abduwasit Ghulam, Kung Hsiangte

(1：新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院，烏魯木齊 830046)(2：新疆大學(xué)綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830046)(3：新疆智慧城市與環(huán)境建模普通高校重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830046)(4：北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院，北京 100875)(5：美國(guó)圣路易斯大學(xué)可持續(xù)發(fā)展中心，圣路易斯 63108)(6：美國(guó)孟菲斯大學(xué)地球科學(xué)系，孟菲斯 38152)

新疆艾比湖流域地表水豐水期和枯水期水質(zhì)分異特征及污染源解析

任 巖1,2,3，張 飛1,2,3，王 娟1,2,3，張 月1,2,3，楊勝天1,4, Abduwasit Ghulam5, Kung Hsiangte6

(1：新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院，烏魯木齊 830046)(2：新疆大學(xué)綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830046)(3：新疆智慧城市與環(huán)境建模普通高校重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830046)(4：北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院，北京 100875)(5：美國(guó)圣路易斯大學(xué)可持續(xù)發(fā)展中心，圣路易斯 63108)(6：美國(guó)孟菲斯大學(xué)地球科學(xué)系，孟菲斯 38152)

采用2015年艾比湖流域54個(gè)采樣點(diǎn)的10個(gè)地表水水質(zhì)指標(biāo)數(shù)據(jù)，首先利用水質(zhì)指數(shù)模型(WQI)和地統(tǒng)計(jì)學(xué)方法對(duì)流域水質(zhì)污染情況進(jìn)行全局評(píng)價(jià)，然后利用層次聚類法、判別分析法和因子分析法分析艾比湖流域地表水豐水期和枯水期水質(zhì)分異特征. 在水質(zhì)時(shí)空分異特征研究的基礎(chǔ)上，利用主成分回歸分析法對(duì)艾比湖流域水質(zhì)進(jìn)行污染源解析. 結(jié)果表明：艾比湖流域豐水期WQI值介于38～70之間，枯水期WQI值介于31～71之間，艾比湖流域豐水期的地表水水質(zhì)污染情況比枯水期嚴(yán)重，而艾比湖、博爾塔拉河和精河靠近艾比湖湖區(qū)的河道污染程度均比其他河道嚴(yán)重. 由聚類分析和判別分析得出艾比湖流域豐水期和枯水期的水質(zhì)采樣點(diǎn)在空間上均被分成A、B兩組，A組包括艾比湖湖區(qū)西部、奎屯河、古爾圖河和四棵樹河，B組包括艾比湖湖區(qū)東部、精河和博爾塔拉河. 艾比湖流域豐水期和枯水期的水體主要受到化學(xué)需氧量、溶解氧、氨氮和懸浮物濃度等指標(biāo)的影響，B組水質(zhì)污染指標(biāo)的值相比于A組的值偏高，B組區(qū)域內(nèi)存在高污染企業(yè)，艾比湖流域水環(huán)境治理工作需主要集中在B組所包括的艾比湖湖區(qū)、博爾塔拉河和精河. ④艾比湖湖區(qū)、精河和博爾塔拉河地表水體的污染主要來(lái)自于有機(jī)物污染和營(yíng)養(yǎng)物質(zhì)污染，其次為工礦業(yè)污染；而奎屯河、古爾圖河、四棵樹河地表水體的污染主要來(lái)自于有機(jī)物污染，其次為營(yíng)養(yǎng)物質(zhì)污染，生物污染的影響較為微弱. 該研究結(jié)果可為艾比湖流域地表水水環(huán)境改善和治理提供一定參考.

艾比湖流域；多元統(tǒng)計(jì)方法；水質(zhì)指數(shù)；時(shí)空分布；污染源解析

內(nèi)陸河流作為區(qū)域居民生活污水、工業(yè)廢水和地表徑流排放的主要載體最易遭受到污染和破壞[1-2]. 近幾十年來(lái)，政府已逐步開展了大量?jī)?nèi)陸河流域水質(zhì)監(jiān)測(cè)項(xiàng)目，并獲得大量的水質(zhì)監(jiān)測(cè)數(shù)據(jù). 由于各個(gè)監(jiān)測(cè)指標(biāo)及監(jiān)測(cè)點(diǎn)之間存在復(fù)雜的相互影響，導(dǎo)致大量監(jiān)測(cè)數(shù)據(jù)并不能充分利用與分析，這也給水質(zhì)專家和地方?jīng)Q策者如何采取有效措施管理和改善水環(huán)境提出了挑戰(zhàn)[3-5]. 因此，從大量的環(huán)境監(jiān)測(cè)數(shù)據(jù)中挖掘有用的信息，探索水質(zhì)的時(shí)空分布模式，識(shí)別潛在污染源和污染源的解析能夠提高人們對(duì)區(qū)域環(huán)境狀況的認(rèn)識(shí)，幫助決策者建立高效合理的水環(huán)境管理方案[6].

近年來(lái)，各類數(shù)學(xué)及統(tǒng)計(jì)評(píng)價(jià)方法被廣泛用于水質(zhì)評(píng)價(jià)、時(shí)空分異、潛在污染源識(shí)別和解析的研究. 水質(zhì)指數(shù)(water quality index，WQI)模型與其他水質(zhì)模型相比具有簡(jiǎn)單易行、靈活多變等特點(diǎn)，對(duì)地表水體及水體富營(yíng)養(yǎng)化污染尤其敏感，在水質(zhì)綜合評(píng)價(jià)中取得了令人滿意的效果[7]. 多元統(tǒng)計(jì)技術(shù)相比于近來(lái)被應(yīng)用較多的神經(jīng)網(wǎng)絡(luò)、模糊綜合評(píng)價(jià)和灰色評(píng)價(jià)法等方法更為簡(jiǎn)便、普適性廣，它不僅能夠?qū)?fù)雜的多源數(shù)據(jù)降維簡(jiǎn)化，又能保證主要信息不丟失，提高分析結(jié)果的可靠性，在實(shí)際應(yīng)用中可以獲得很好的效果[8]. 聚類分析、判別分析、因子分析和主成分回歸分析作為傳統(tǒng)的多元統(tǒng)計(jì)技術(shù)，在水質(zhì)時(shí)空分異特征及潛在污染源識(shí)別、解析中得到廣泛應(yīng)用[9]. 國(guó)內(nèi)外學(xué)者利用多元統(tǒng)計(jì)技術(shù)分別對(duì)太湖流域[10]、陜西灃河[11]、沁河流域[12]、洞庭湖[13]、土耳其中部近海[14]、溫瑞塘河流域[15]以及德國(guó)北部低洼地區(qū)[16]的污染物時(shí)空分異特征及潛在污染源識(shí)別進(jìn)行了研究，并取得了滿意的效果. 內(nèi)陸河流域的水質(zhì)時(shí)空分異特征和污染源解析研究較成熟[17-21]，也有學(xué)者利用多元統(tǒng)計(jì)方法評(píng)估了近海海域的地表水質(zhì)量[22]，分析了地下水的水質(zhì)狀況[23]，此外，Mohammadpour等[24]利用WQI和改進(jìn)模型分析了人工濕地方面水質(zhì)的時(shí)空特征. 但這些研究割裂了時(shí)間與空間相互作用的機(jī)制，只是單獨(dú)探討了時(shí)間或空間上分異特征及污染源識(shí)別，忽略了時(shí)間對(duì)空間分布規(guī)律及其污染源的影響.

艾比湖流域是古絲綢之路的要道，同時(shí)也是準(zhǔn)噶爾盆地、天山北坡經(jīng)濟(jì)帶，水環(huán)境的好壞直接影響到第二條亞歐大陸橋和國(guó)家級(jí)二類口岸阿拉山口的發(fā)展以及該流域的可持續(xù)發(fā)展[25]. 作為內(nèi)陸河流域，艾比湖流域?qū)^(qū)域工農(nóng)業(yè)發(fā)展及生態(tài)環(huán)境保護(hù)起著重要作用. 以往的研究多是對(duì)艾比湖流域進(jìn)行水質(zhì)評(píng)價(jià)，而沒(méi)有分析水質(zhì)的時(shí)空分異特征和污染源解析，因此本文選擇5月和10月的地表水采樣點(diǎn)，分析艾比湖流域豐水期和枯水期的地表水水污染時(shí)空分異特征，并聯(lián)合水質(zhì)時(shí)空因子進(jìn)行污染源解析. 故本文首先根據(jù)艾比湖流域河網(wǎng)水質(zhì)監(jiān)測(cè)數(shù)據(jù)，應(yīng)用水質(zhì)指數(shù)WQI模型對(duì)艾比湖流域河網(wǎng)水質(zhì)進(jìn)行綜合評(píng)價(jià)；然后采用聚類分析對(duì)艾比湖流域豐水期和枯水期的水質(zhì)采樣點(diǎn)進(jìn)行空間分組，同時(shí)利用判別分析驗(yàn)證聚類分析空間分組結(jié)果的可信度；其次基于聚類分析結(jié)果，采用因子分析法研究艾比湖流域水質(zhì)時(shí)空分異特征及識(shí)別顯著性污染指標(biāo)；最后結(jié)合時(shí)空聯(lián)合因子分析的結(jié)果采用主成分回歸分析中基于受體的源分配模型，識(shí)別不同時(shí)期不同區(qū)域水環(huán)境污染源以及主要污染源對(duì)河道水質(zhì)的影響，為艾比湖流域地表水水環(huán)境治理提供科學(xué)依據(jù).

1 研究區(qū)域和研究方法

1.1 研究區(qū)概況

艾比湖流域(43°38′～45°52′N，79°53′～85°02′E)包括奎屯河、四棵樹河、古爾圖河、精河、阿恰勒河、大河沿子河、博爾塔拉河7條河流，其中博河、精河和奎屯河是主要的3大支流，地跨博爾塔拉蒙古自治州的博樂(lè)市、溫泉縣和精河縣，塔城地區(qū)的烏蘇市和托里縣南部，伊犁直屬的奎屯市和克拉瑪依的獨(dú)山子區(qū). 流域面積50621 km2，其中山地面積24317 km2，平原區(qū)面積25762 km2，湖泊面積542 km2. 流域西、南、北三面環(huán)山，中間為谷地平原，東部有尾閭艾比湖，是一個(gè)具有典型干旱區(qū)山地-綠洲-荒漠生態(tài)環(huán)境特點(diǎn)的區(qū)域，湖區(qū)多年平均降水不足100 mm，而山區(qū)降水相對(duì)較多，是湖區(qū)的主要水源補(bǔ)給區(qū)[26]. 近年來(lái)，由于人類大規(guī)模的墾荒活動(dòng)，導(dǎo)致艾比湖不斷涸縮，湖水量減少，目前只有博爾塔拉河、精河、奎屯河，分別從西、南、東3個(gè)方向注入艾比湖，成為湖水的主要來(lái)源. 湖面呈橢圓狀，面積650 km2，平均水深為2～3 m，湖面海拔189 m[27]. 艾比湖流域包括博爾塔拉蒙古自治州、奎屯市、獨(dú)山子、烏蘇市、托里縣，是北疆地區(qū)重要的棉花、糧食、畜牧業(yè)基地和石油化工基地. 研究區(qū)水質(zhì)監(jiān)測(cè)采樣點(diǎn)和水系分布如圖1所示.

1.2 數(shù)據(jù)收集

2015年5月和10月分別采集艾比湖流域54個(gè)采樣點(diǎn)(圖1)的地表水樣品. 選取化學(xué)需氧量(COD)、溶解氧(DO)濃度、五日生化需氧量(BOD5)、總氮(TN)濃度、總磷(TP)濃度、氨氮(NH3-N)濃度、pH、懸浮物(SS)濃度、色度、濁度10個(gè)指標(biāo)進(jìn)行測(cè)定. 所有樣本的采集及檢測(cè)分析嚴(yán)格按照《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838-2002)中的方法進(jìn)行.

圖1 艾比湖流域地理位置Fig.1 Geographical location of Lake Ebinur Watershed

1.3 研究方法

1.3.1WQI模型 水質(zhì)指數(shù)是綜合了多個(gè)水質(zhì)參數(shù)的歸一化無(wú)量綱數(shù)，需要考慮的參數(shù)取決于水的用途，常考慮的水質(zhì)參數(shù)包括DO濃度、BOD5、COD、pH值、SS濃度、NH3-N濃度、TN濃度和TP濃度等[24]，該模型較為成熟，在國(guó)內(nèi)也得到廣泛應(yīng)用[28]. 本文主要利用該模型對(duì)整個(gè)流域的水質(zhì)污染情況進(jìn)行全局評(píng)價(jià)，參考Mohammadpour等[24]關(guān)于馬來(lái)西亞濱州濕地的綜合水質(zhì)模型WQI，選擇了DO濃度、BOD5、COD濃度、pH、SS濃度、NH3-N濃度6個(gè)水質(zhì)參數(shù)，將這些水質(zhì)參數(shù)變量根據(jù)分類指數(shù)函數(shù)(SI)轉(zhuǎn)變成非空間變量，不同參數(shù)的SI計(jì)算公式見表1.

表1 水質(zhì)指數(shù)的分類指數(shù)函數(shù)[24]

aX代表的變量單位是mg/L，其中pH無(wú)量綱.

WQI模型計(jì)算公式為[23]：

WQI=0.22SIDO+0.19SIBOD5+0.16SICOD+0.15SINH3-N+0.16SISS+0.12SIpH

(1)

式中，SI代表分類指數(shù). 模型根據(jù)輸入水質(zhì)監(jiān)測(cè)數(shù)據(jù)將采樣點(diǎn)水質(zhì)歸一化到0～100之間，并根據(jù)水體污染的程度，從0～100將水質(zhì)由差到極好劃分為6類，分別對(duì)應(yīng)嚴(yán)重污染(<31.0)、污染(31.0～51.9)、輕度污染(51.9～76.5)、較好(76.5～92.7)和好(>92.7)[24].

1.3.2 層次聚類分析 聚類分析(CA)是根據(jù)對(duì)象距離遠(yuǎn)近或相似性大小進(jìn)行分類的多元統(tǒng)計(jì)方法. 本文利用該方法對(duì)采樣點(diǎn)的水質(zhì)數(shù)據(jù)做空間分類，采用的層次聚類分析(HCA)是應(yīng)用最廣泛的聚類方法，該法以逐次聚合的方式將距離最近或者最相似的對(duì)象聚成一個(gè)類簇，直至最后聚成一類[29-30]. 本文采用的計(jì)算方法是歐氏距離平方和離差平方法.

聚類分析要求數(shù)據(jù)符合正態(tài)分布，在進(jìn)行分析之前，需要利用K-S非參數(shù)檢驗(yàn)[29]對(duì)采樣點(diǎn)水質(zhì)數(shù)據(jù)進(jìn)行正態(tài)檢驗(yàn). 結(jié)果顯示，采樣點(diǎn)各監(jiān)測(cè)指標(biāo)除BOD5外均以95%或更高的可信度服從正態(tài)分布；采樣點(diǎn)各指標(biāo)偏離正態(tài)分布，進(jìn)行自然對(duì)數(shù)轉(zhuǎn)換后[30]，各監(jiān)測(cè)指標(biāo)除pH和SS外均能夠以95%或更高的可信度服從正態(tài)分布. 為了消除變量單位量綱的影響，需要對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理(均值為0，方差為1).

1.3.3 判別分析 判別分析(DA) 是多變量統(tǒng)計(jì)分析中用于判別樣品所屬類別的一種統(tǒng)計(jì)分析方法. 可以用來(lái)判別CA分析結(jié)果和識(shí)別顯著性的污染指標(biāo). 此方法分為3類：標(biāo)準(zhǔn)式、前進(jìn)式和后退式. 相比之下，后退式DA方法具有更好的指標(biāo)降維能力和判別能力[31]. 所以，本文采用該方法對(duì)1.3.2節(jié)的聚類分析做驗(yàn)證分析，利用后退式DA方法對(duì)原始數(shù)據(jù)進(jìn)行時(shí)空差異性分析，并采用交叉驗(yàn)證法(Cross-validation)檢驗(yàn)此方法的判別能力. 相應(yīng)的判別函數(shù)表達(dá)式為：

(2)

式中，f(Gi)表示第i類的判別函數(shù)，Ci為第i類的固有常數(shù)，n表示參與判別分析指標(biāo)個(gè)數(shù)，pij表述第i類第j個(gè)指標(biāo)值，wij表示對(duì)應(yīng)的判別系數(shù).

1.3.4 因子分析法 因子分析(FA)是一種既可以降低變量維數(shù)，又可以對(duì)變量進(jìn)行分類的廣泛應(yīng)用的方法[32]. 其實(shí)質(zhì)是從多個(gè)實(shí)測(cè)的原變量中提取出較少的、互不相關(guān)的抽象綜合指標(biāo)，即因子. 每個(gè)原變量可用這些提取出的公共因子的線性組合表示，具體參見公式(3). 同時(shí)，根據(jù)各個(gè)因子對(duì)原變量的影響大小，也可將原變量劃分為等同于因子數(shù)目的類數(shù). 在水質(zhì)分析中，此方法主要用于提取污染因子和識(shí)別污染源[33]，本文在聚類分析和判別分析的基礎(chǔ)上，利用該方法分析各類水質(zhì)指標(biāo)的主要污染因子.

Zk=ak1·F1+ak2·F2+,…,+akm·Fm+Uk(k=1、2、3、…、n)

(3)

式中，Zk表示第k個(gè)變量的標(biāo)準(zhǔn)化分?jǐn)?shù)，akm表示第m個(gè)因子相對(duì)于第k個(gè)變量的因子載荷，F(xiàn)m表示第m個(gè)公共因子，Uk為Zk的唯一因素，m表示所有變量公共因子的數(shù)目，n表示變量的個(gè)數(shù).

因子分析要求變量間要有一定的相關(guān)性，所以在數(shù)據(jù)分析之前需要利用巴特萊檢驗(yàn)(KMO)[29]檢驗(yàn)數(shù)據(jù)應(yīng)用因子分析的可行性. 同時(shí)，考慮到水質(zhì)指標(biāo)數(shù)量級(jí)上的差異，需要對(duì)數(shù)據(jù)再進(jìn)行標(biāo)準(zhǔn)化(均值為0，方差為1). 一般來(lái)說(shuō)，分析結(jié)果只選取特征值大于1的因子[33]. 因子分析及判別分析均采用原始數(shù)據(jù).

1.3.5 主成分回歸分析 主成分回歸分析法(APCS-MLR)是一種基于因子得分評(píng)價(jià)各個(gè)因子對(duì)各個(gè)變量貢獻(xiàn)的統(tǒng)計(jì)方法. 其原理是將變量值與因子得分進(jìn)行多元線性回歸，根據(jù)回歸參數(shù)得到針對(duì)各個(gè)因子的估計(jì)值，從而確定因子對(duì)各個(gè)變量的貢獻(xiàn)，具體參見公式(4)[34]. 目前，這種方法在水質(zhì)評(píng)價(jià)中用于計(jì)算污染因子對(duì)各個(gè)評(píng)價(jià)指標(biāo)的貢獻(xiàn)率[35]，本文在因子分析的基礎(chǔ)上，對(duì)各個(gè)主要污染因子分析水質(zhì)污染的貢獻(xiàn)率.

(4)

式中，Ms表示第s個(gè)變量的實(shí)測(cè)值，aos表示第s個(gè)變量的多元回歸的常數(shù)項(xiàng)，Aps表示第p個(gè)因子對(duì)第s個(gè)變量的回歸系數(shù)，APCSp表示調(diào)整后的第p個(gè)因子的分?jǐn)?shù)，n表示因子個(gè)數(shù).Aps·APCSp表示第p個(gè)因子對(duì)Ms的質(zhì)量濃度貢獻(xiàn)，所有樣本的Aps·APCSp的平均值就是表示因子平均絕對(duì)貢獻(xiàn)率. 研究中多元統(tǒng)計(jì)分析及水體綜合水質(zhì)空間分析采用的是Microsoft Excel 2010、SPSS 19和ArcGIS 10.0軟件.

2 結(jié)果與分析

2.1 水質(zhì)指標(biāo)的統(tǒng)計(jì)特征

5月屬于豐水期，10月屬于枯水期[14,36-37]. 總體上來(lái)看，艾比湖流域的地表水水環(huán)境污染比較嚴(yán)重. 豐水期，COD、TP、TN、NH3-N、BOD5和DO濃度均值分別為46.93、0.07、1.97、0.63、1.85和4.46 mg/L，根據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》，DO屬于Ⅳ類水質(zhì)標(biāo)準(zhǔn)，TN屬于Ⅴ類水質(zhì)標(biāo)準(zhǔn)，NH3-N屬于Ⅲ類水質(zhì)標(biāo)準(zhǔn)，TP屬于Ⅱ類水質(zhì)標(biāo)準(zhǔn)，COD超出Ⅴ類水質(zhì)標(biāo)準(zhǔn)的1.17倍，BOD5則分別屬于Ⅰ類水質(zhì)標(biāo)準(zhǔn). 枯水期，COD、TP、TN、NH3-N、BOD5和DO濃度的均值分別為15.58、0.08、1.06、0.08、2.09和8.29 mg/L，與《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中相應(yīng)指標(biāo)Ⅴ類水標(biāo)準(zhǔn)進(jìn)行比較可知，DO屬于Ⅰ類水標(biāo)準(zhǔn)，TN屬于Ⅳ類水標(biāo)準(zhǔn)，NH3-N屬于Ⅰ類水質(zhì)標(biāo)準(zhǔn)，TP屬于Ⅱ類水質(zhì)標(biāo)準(zhǔn)，COD、 BOD5則分別屬于Ⅰ類和Ⅱ類水質(zhì)標(biāo)準(zhǔn)(表2). 綜上所述，豐水期比枯水期的水質(zhì)污染嚴(yán)重.

2.2WQI水質(zhì)綜合評(píng)價(jià)

通過(guò)計(jì)算得到，豐水期WQI值的范圍是38～70，枯水期WQI值的范圍是31～71. 按照采樣點(diǎn)的數(shù)量計(jì)算，豐水期有18.50%處于污染，83.30%處于輕度污染. 枯水期有5.20%處于污染，94.70%處于輕度污染. 總體來(lái)看，豐水期的艾比湖流域水質(zhì)污染情況較為嚴(yán)重，相比較來(lái)看枯水期水質(zhì)情況較好，普遍只有輕度污染.

從WQI采樣斷面豐水期和枯水期的空間分布(圖2)可知，豐水期的污染區(qū)域主要分布在艾比湖湖區(qū)，精河、博爾塔拉河和奎屯河靠近艾比湖的支流；枯水期的污染區(qū)域主要分布在艾比湖湖區(qū)；豐水期的污染區(qū)域大于枯水期；枯水期的輕度污染區(qū)域大于豐水期. 總體來(lái)看，豐水期和枯水期的艾比湖湖區(qū)和周邊的河道污染程度要比外圍河道嚴(yán)重，主要分布在艾比湖湖區(qū)、博爾塔拉河、精河和奎屯河靠近艾比湖的河道. 豐水期的污染情況要比枯水期嚴(yán)重，這可能因?yàn)榘群饔虻呢S水期是施肥和噴灑農(nóng)藥的主要季節(jié)，艾比湖流域主要支流流經(jīng)精河縣、博樂(lè)市、溫泉縣和烏蘇縣等地的農(nóng)田綠洲，這些地區(qū)的農(nóng)耕活動(dòng)會(huì)影響沿途支流的水質(zhì)[38-39].

表2 豐水期和枯水期水質(zhì)指標(biāo)的統(tǒng)計(jì)描述

圖2 WQI采樣斷面豐水期和枯水期空間分布Fig.2 The spatial distribution of WQI during wet and dry seasons

2.3 水質(zhì)污染特征的時(shí)空相似性與差異性

豐水期空間尺度的聚類分析和判別分析表明，空間上可以將采樣點(diǎn)分為2組(圖3)：A組包括P1、P2、P4、P5、P6、P7、P10、P11、P12、P14、P15、P17、P18、P19、P22、P30、P31、P32、P33、P34、P35、P45、P49和P50這24個(gè)采樣點(diǎn)，主要位于奎屯河、古爾圖河、四棵樹河和精河的部分支流；B組包括P3、P8、P9、P13、P16、P20、P21、P23、P24、P25、P26、P27、P28、P29、P36、P37、P38、P39、P40、P41、P42、P43、P44、P46、P47、P48、P51、P52、P53和P54這30個(gè)采樣點(diǎn)，主要位于博爾塔拉河、精河和艾比湖湖區(qū). 采樣點(diǎn)聚類分析結(jié)果的判別分析交叉驗(yàn)證正確率達(dá)到95.70%(表3)，誤判區(qū)的區(qū)域主要集中在艾比湖湖區(qū)和3大主要支流的連接河段，總體分類結(jié)果較好.

表3 水質(zhì)污染的空間判別分析*

*卡方檢驗(yàn)，Sig.小于0.05.

枯水期空間尺度的聚類分析和判別分析表明，空間上也是將采樣點(diǎn)分為2組(圖3)：A組包括Q1、Q4、Q6、Q7、Q10、Q11、Q12、Q19、Q21、Q23、Q27、Q28、Q29、Q32、Q34、Q35、Q37、Q41、Q48、Q49、Q50這21個(gè)采樣點(diǎn)，主要位于奎屯河、古爾圖河和四棵樹河；B組包括Q2、Q3、Q5、Q8、Q9、Q13、Q14、Q15、Q16、Q17、Q18、Q20、Q22、Q24、Q25、Q26、Q30、Q31、Q33、Q36、Q38、Q39、Q40、Q42、Q43、Q44、Q45、Q46、Q47、Q51、Q52、Q53、Q54這33個(gè)采樣點(diǎn)，主要位于博爾塔拉河、精河和艾比湖湖區(qū). 采樣點(diǎn)聚類分析結(jié)果的判別分析交叉驗(yàn)證正確率達(dá)到86.20%(表3)，誤判區(qū)的區(qū)域主要集中在艾比湖湖區(qū)和3大主要支流的連接河段，總體分類結(jié)果較好.

圖3 艾比湖流域不同時(shí)期水質(zhì)污染的空間尺度聚類分析Fig.3 Spatial cluster analysis of water pollution in Lake Ebinur Watershed

根據(jù)艾比湖流域豐水期和枯水期的水質(zhì)聚類結(jié)果，做了A、B兩組分區(qū)的空間分布圖(圖1). 從圖中可以看出，豐水期A組和枯水期A組主要分布在奎屯河、古爾圖河、四棵樹河和艾比湖湖區(qū)西部；豐水期B組和枯水期B組主要分布在博河、精河和艾比湖湖區(qū)東部.

從艾比湖流域水質(zhì)污染物同一時(shí)期不同分組的差異性(圖4)可以看出：豐水期A組的COD和DO濃度較高，污染主要來(lái)源于有機(jī)物，主要反映了與有機(jī)物相關(guān)的因素對(duì)地表水體的影響. 該區(qū)域地勢(shì)較為平坦，河道水流長(zhǎng)年滯流，水體復(fù)氧能力差[40]，這可能是COD和DO濃度較高的主要原因. 豐水期B組的SS和COD濃度較高，該區(qū)域范圍主要包括溫泉縣、博樂(lè)市和精河縣的兩大支流博爾塔拉河和精河，沿線有工礦企業(yè)等，污染主要來(lái)源于工礦企業(yè)的廢水污染超標(biāo)(數(shù)據(jù)來(lái)自博州水利局官網(wǎng)http:∥slj.xjboz.gov.cn/). 枯水期A組與豐水期A組的情況類似，COD和DO濃度相對(duì)較高. 枯水期B組也主要分布在博河、精河，該區(qū)域的COD、DO和NH3-N濃度均較高，首先因?yàn)樵搮^(qū)域有工礦企業(yè)和污水處理廠，其次由于是枯水期流域水量較小，水流較緩，該區(qū)域?qū)儆邴}漬化比較嚴(yán)重的地區(qū)，10月地表植被已經(jīng)枯萎，導(dǎo)致鹽漬化程度更加嚴(yán)重，這也是COD濃度較高的原因. NH3-N濃度較高主要是因?yàn)橹Я髁鹘?jīng)農(nóng)業(yè)區(qū)，10月是春夏作物和秋冬作物的交替時(shí)期，主要污染物來(lái)源于農(nóng)業(yè)營(yíng)養(yǎng)物質(zhì)、畜牧業(yè)污染、生活污水和污水處理廠等的影響[1].

總體來(lái)看，艾比湖流域的水質(zhì)污染情況在空間上具有一定的相似性，豐水期和枯水期的COD濃度均相較于其他水質(zhì)指標(biāo)偏高，B組水質(zhì)污染指標(biāo)的值相比于A組的值較高. 在空間上也有一定差異性，豐水期A組的COD濃度相較于B組偏高， B組的BOD5、NH3-N和SS濃度相較于A組偏高；枯水期B組的COD、BOD5和NH3-N濃度比A組偏高，A組的SS和DO濃度比B組偏高. 艾比湖流域地表水環(huán)境治理主要集中在艾比湖湖區(qū)，需要引清入渠和河道疏通等，同時(shí)需要對(duì)兩大支流流經(jīng)的工礦企業(yè)污水排放等加強(qiáng)監(jiān)督和管理.

圖4 艾比湖流域水質(zhì)指標(biāo)同一時(shí)期不同分組的差異性Fig.4 The water quality indexes of different groups in the same period of Lake Ebinur Watershed

2.4 時(shí)空聯(lián)合因子對(duì)污染源的識(shí)別和解析

空間聚類結(jié)果判別分析所建立的判別函數(shù)對(duì)不同聚類分組進(jìn)行判別驗(yàn)證，發(fā)現(xiàn)空間上A組和B組的判別正確率分別為91.40%、91.70%，總判別正確率為95.70%(表3). 說(shuō)明同一區(qū)域在不同時(shí)期或者同一時(shí)期內(nèi)不同區(qū)域的水質(zhì)指標(biāo)分布規(guī)律存在明顯差異. 同時(shí)，從艾比湖流域水質(zhì)指標(biāo)各空間分組在不同時(shí)期的空間差異(圖5)可以看出，空間A組和B組在不同時(shí)期內(nèi)的污染水平有著較大的差異. 因此，有必要將時(shí)間及空間聚類分析結(jié)果有機(jī)結(jié)合，對(duì)不同區(qū)域在不同時(shí)段下的污染源進(jìn)行解析，為水環(huán)境管理和改善提供更為詳盡的理論依據(jù).

圖5 艾比湖流域水質(zhì)指標(biāo)同一分組不同時(shí)期的差異性Fig.5 The water quality indexes of the same group in different periods of Lake Ebinur Watershed

2.4.1 空間分組各時(shí)段因子分析及污染源識(shí)別 不同空間分組各時(shí)段的水質(zhì)指標(biāo)相關(guān)系數(shù)矩陣特征值、貢獻(xiàn)率和累計(jì)貢獻(xiàn)率，如表4所示. 根據(jù)表4選出主因子，并對(duì)初始主因子載荷矩陣進(jìn)行方差極大化正交旋轉(zhuǎn)得到載荷矩陣因子(表5)，選擇載荷矩陣因子≥0.65[41]的歸為主因子，一般是由多項(xiàng)水質(zhì)指標(biāo)組合而成. 利用KMO檢驗(yàn)和Bartlett球形方法對(duì)地表水體標(biāo)準(zhǔn)化的監(jiān)測(cè)數(shù)據(jù)進(jìn)行相關(guān)矩陣檢驗(yàn). 豐水期、枯水期的空間A組和B組的KMO值均是0.52，Bartlett的顯著性也是0.00，滿足P<0.05置信度，檢驗(yàn)結(jié)果表明因子分析是有效的(表5). 根據(jù)特征因子>1的原則[33]對(duì)空間A組在豐水期和枯水期分別提取3個(gè)因子(最大方差旋轉(zhuǎn))，累積解釋方差分別為82.31%和72.53%(表5)；根據(jù)特征因子>1的原則[33]對(duì)空間B組在豐水期和枯水期分別提取3個(gè)和2個(gè)因子(最大方差旋轉(zhuǎn))，累積解釋方差分別為74.30%和76.85%(表5).

表4 各水質(zhì)指標(biāo)相關(guān)系數(shù)矩陣特征值、貢獻(xiàn)率和累計(jì)貢獻(xiàn)率

豐水期A組：F1的方差貢獻(xiàn)率為36.74%，表征因子為COD和BOD5，并且都呈正相關(guān)，同時(shí)由表6可知COD和BOD5之間有顯著相關(guān)性，Pearson相關(guān)系數(shù)都在0.7以上，表明F1主要代表生物、化學(xué)需氧量相關(guān)污染源，而此類污染主要來(lái)源于有機(jī)物，因此F1主要反映了與有機(jī)物相關(guān)的生活用水和第三產(chǎn)業(yè)廢水的污染對(duì)地表水體的影響；F2的方差貢獻(xiàn)率為26.95%，表征因子為NH3-N，呈正相關(guān)，此類主要來(lái)源于過(guò)量的營(yíng)養(yǎng)物質(zhì)污染，根據(jù)塔城地區(qū)環(huán)境保護(hù)局官網(wǎng)(http:∥www.tcdqhb.gov.cn)公布的2015年上半年廢水排放污染企業(yè)名單可知，本區(qū)有烏蘇市天然番茄制品有限公司，因此可以推斷F2主要因?yàn)楣I(yè)和生活污水排入河道，導(dǎo)致營(yíng)養(yǎng)物質(zhì)過(guò)剩，同時(shí)也跟這個(gè)時(shí)期大量的農(nóng)業(yè)活動(dòng)有關(guān)；F3的方差貢獻(xiàn)率為18.62%，表征因子為SS，呈正相關(guān)，推斷F3主要來(lái)源于傍河而立的高污染企業(yè)廢水廢料的排放.

枯水期A組：F1的方差貢獻(xiàn)率為31.30%，表征因子為COD和pH，其中COD呈正相關(guān)，pH呈負(fù)相關(guān)，由表6可知COD和pH之間也呈顯著負(fù)相關(guān)，說(shuō)明它們之間也是相互影響的，因此F1代表著有機(jī)物污染源，隨著有機(jī)物污染的加重，水質(zhì)pH值呈現(xiàn)負(fù)增加；F2的方差貢獻(xiàn)率為23.17%，表征因子為BOD5和DO，其中BOD5呈負(fù)相關(guān)，DO呈正相關(guān)，由表6可知兩者相關(guān)性一般，說(shuō)明兩者之間不影響，跟豐水期A組的F1基本一致，不同的是本時(shí)段DO在第2因子中占有較大的載荷，說(shuō)明這一時(shí)段水流量減少，河道中出現(xiàn)滯留，水體復(fù)氧能力差，導(dǎo)致水體中DO濃度減少，因此F2代表著生物、化學(xué)污染物；F3的方差貢獻(xiàn)率為18.06%，表征因子為SS，根據(jù)伊犁州環(huán)境保護(hù)局官網(wǎng)(http:∥www.ylepa.gov.cn/index.aspx)公布的2015年下半年廢水排放污染企業(yè)名單可知，本區(qū)有新疆天玉生物科技有限公司，推斷本區(qū)水質(zhì)污染主要源于生物科技公司的廢水、廢料排放.

表5 時(shí)空聯(lián)合因子分析因子載荷矩陣

表6 豐水期和枯水期水質(zhì)指標(biāo)間的Pearson相關(guān)系數(shù)

*表示在0.05水平(雙側(cè))上顯著相關(guān); **表示在0.01水平(雙側(cè))上顯著相關(guān).

豐水期B組：F1的方差貢獻(xiàn)率為33.64%，表征因子為COD、BOD5和DO，均呈正相關(guān)，由表7可知COD和BOD5呈顯著相關(guān)，兩者與DO都呈顯著正相關(guān)，跟枯水期A組的F2基本一致，不同的是本時(shí)段的COD在第2因子中占有較大的載荷，說(shuō)明隨著BOD5的增加，DO、COD也跟著增加，因此F1主要代表著營(yíng)養(yǎng)物質(zhì)、耗氧有機(jī)物以及生物、化學(xué)需氧量相關(guān)的污染物；F2的方差貢獻(xiàn)率為22.97%，表征因子為NH3-N，與其他各水質(zhì)指標(biāo)的相關(guān)性都一般(表6)，表明F2代表著季節(jié)性污水排放的營(yíng)養(yǎng)物質(zhì)過(guò)量污染；F3的方差貢獻(xiàn)率為17.69%，表征因子為SS，同豐水期A組的F2，本區(qū)域博爾塔拉蒙古自治州西部礦業(yè)有限公司和博爾塔拉蒙古自治州西域福樂(lè)礦業(yè)有限責(zé)任公司的排放物分別超標(biāo)2倍和4倍(數(shù)據(jù)來(lái)源于博州環(huán)境保護(hù)局官網(wǎng)http:∥hbj.xjboz.gov.cn/hjzl.htm)，因此推斷F3主要代表工礦業(yè)排放物的污染.

枯水期B組：F1的方差貢獻(xiàn)率為54.75%，表征因子為COD、BOD5、NH3-N和SS，由表6可知，4種指標(biāo)的相關(guān)性均較高，跟豐水期A組的F1基本一致，不同的是本時(shí)段NH3-N和SS在第1因子中占有較大的載荷，說(shuō)明COD和BOD5濃度的增加會(huì)影響NH3-N和SS濃度，但因?yàn)楸緟^(qū)域溫泉縣供排水有限責(zé)任公司排出的廢水中SS濃度超過(guò)城鎮(zhèn)污水處理廠污染物排放標(biāo)準(zhǔn)的1倍(數(shù)據(jù)來(lái)自于博州水利局官網(wǎng)公開監(jiān)測(cè)數(shù)據(jù))，因此F1主要代表營(yíng)養(yǎng)物質(zhì)、耗氧有機(jī)物和工礦企業(yè)廢水污染；F2的方差貢獻(xiàn)率為22.11%，表征因子為DO，因此F2主要代表生物化學(xué)污染源.

依據(jù)因子分析，在豐水期和枯水期艾比湖流域主要受到COD、DO、NH3-N和SS濃度等指標(biāo)的影響，空間B組整體比空間A組污染因子多，豐水期比枯水期的污染情況嚴(yán)重. 從水質(zhì)時(shí)空聯(lián)合因子分析來(lái)看，艾比湖流域水環(huán)境治理工作需主要集中在空間B組所包括的艾比湖湖區(qū)、精河和博爾塔拉河部分河道.

2.4.2 污染源解析和貢獻(xiàn)率估算 通過(guò)上述分析，確定了污染源的數(shù)量和特征，進(jìn)一步利用主成分回歸分析法計(jì)算污染源貢獻(xiàn)率. E/O是回歸分析中估計(jì)值與實(shí)測(cè)值的比值，主要集中在1.00左右，表明了回歸結(jié)果較為準(zhǔn)確[2]. A、B兩組中SS濃度和pH兩種水質(zhì)指標(biāo)對(duì)各項(xiàng)污染源的貢獻(xiàn)率的R2均小于0.5，這說(shuō)明這兩種指標(biāo)對(duì)各項(xiàng)污染源的貢獻(xiàn)率不大，那么主要對(duì)COD、BOD5、DO和NH3-N 4種水質(zhì)指標(biāo)的污染貢獻(xiàn)率進(jìn)行分析. 空間A組地表水體的污染主要來(lái)源于有機(jī)物污染和營(yíng)養(yǎng)物質(zhì)污染，其中有機(jī)物污染主要影響 BOD5和COD濃度，定量化的污染貢獻(xiàn)率分別為57.60%和55.30%，營(yíng)養(yǎng)物質(zhì)污染主要影響COD、DO和NH3-N濃度，貢獻(xiàn)率分別為41.60%、30.30%和50.80%；空間B組地表水體的污染主要來(lái)源于有機(jī)物污染，其次為營(yíng)養(yǎng)物質(zhì)污染，生物污染的影響較為微弱；有機(jī)物污染主要影響COD、BOD5和DO濃度，貢獻(xiàn)率分別為64.40%、47.60%和20.20%，營(yíng)養(yǎng)物質(zhì)的污染主要影響該區(qū)域NH3-N和COD濃度，其影響貢獻(xiàn)率為63.10%和21.40%(表7).

表7 艾比湖流域水質(zhì)指標(biāo)各污染源貢獻(xiàn)率

3 討論

本文選擇WQI模型，首次用于評(píng)價(jià)艾比湖流域的水質(zhì)污染情況，主要因?yàn)樵撃Ｐ团c其他水質(zhì)指數(shù)模型相比的優(yōu)勢(shì)是運(yùn)算簡(jiǎn)單、參數(shù)設(shè)定可以結(jié)合研究區(qū)水質(zhì)保護(hù)提升目標(biāo)而設(shè)定，而且也已經(jīng)被國(guó)內(nèi)外專家學(xué)者成熟應(yīng)用在地表水和干旱半干旱區(qū)的水質(zhì)研究. 其中，李義祿等[2]和楊柳等[28]利用WQI模型研究?jī)?nèi)陸湖泊地表水污染時(shí)空分異特征，Bassam等[42]利用該模型研究干旱半干旱亞地區(qū)伊拉克西北部Sinjar的水質(zhì)污染情況. 研究區(qū)艾比湖流域也屬于干旱半干旱地區(qū)，與伊拉克的生態(tài)環(huán)境相似，所以本文選擇WQI模型來(lái)研究艾比湖流域的地表水水質(zhì)污染狀況.

根據(jù)前人相關(guān)研究可知，遼河流域[43]、晉江流域[44]、三峽流域[45]和太湖流域[10]等內(nèi)陸河流域在水質(zhì)污染源解析的研究較多，但是處于干旱區(qū)的艾比湖流域在這方面的研究較少，因此本文嘗試探討該地區(qū)的水質(zhì)污染源解析以及水質(zhì)時(shí)空分異特征. 由彌艷等[27]和博樂(lè)市水文水資源勘測(cè)大隊(duì)的調(diào)查研究[46]可知，艾比湖流域兩大主要支流(精河和博河)的入湖量在豐水期和枯水期差異較大，因此本文選擇豐水期和枯水期兩個(gè)典型時(shí)期對(duì)艾比湖流域的水質(zhì)污染情況進(jìn)行評(píng)價(jià)，這樣能夠較全面地解析艾比湖湖區(qū)及周邊水質(zhì)的污染源. 本文針對(duì)艾比湖流域豐水期和枯水期的水質(zhì)分異特征及污染源解析結(jié)果得出，艾比湖流域豐水期的水質(zhì)污染比枯水期嚴(yán)重，這與彌艷等[36,39]在2008年的研究結(jié)果一致. 彌艷等研究還發(fā)現(xiàn)2008年艾比湖豐水期的NH3-N、TP、TN濃度均較高，而且TN屬于V類水質(zhì)標(biāo)準(zhǔn)，整個(gè)艾比湖TP、TN濃度達(dá)到了富營(yíng)養(yǎng)化水平；但本文的研究發(fā)現(xiàn)：雖然豐水期的TN仍然達(dá)到Ⅴ類水質(zhì)標(biāo)準(zhǔn)，但NH3-N屬于Ⅲ類水質(zhì)標(biāo)準(zhǔn)，TP屬于Ⅱ類水質(zhì)標(biāo)準(zhǔn). 與2008年相比，2015年艾比湖流域的水質(zhì)明顯得到改善，這與國(guó)家的水質(zhì)改善和環(huán)保政策密不可分，也說(shuō)明兩個(gè)保護(hù)區(qū)的設(shè)立以及國(guó)家和地區(qū)的相關(guān)保護(hù)措施效果顯著.

針對(duì)本文的研究?jī)?nèi)容和結(jié)果，從自然因素和人文因素兩方面對(duì)艾比湖流域的地表水水質(zhì)污染原因進(jìn)行討論.

1)自然因素：艾比湖是新疆最大的鹽水湖，根據(jù)張飛等[26]、吳超存[47]、董煜等[48]和王娟等[49]的研究可知：近20年來(lái)，溫度、降水量和蒸發(fā)量逐年增長(zhǎng)，冰川融水導(dǎo)致上游徑流量增加，需水量的增加導(dǎo)致下游徑流量減少，流域徑流量總體上變化不明顯，但水資源總量在減少. 流域內(nèi)耕地面積增加，水體、林草地面積減少，水資源需求量增加，致使流入艾比湖的主要河流水量減少甚至奎屯河出現(xiàn)斷流的現(xiàn)象，艾比湖湖面萎縮、湖底裸露，湖內(nèi)水流較緩，湖周鹽漬化現(xiàn)象嚴(yán)重，周邊沙漠化面積增加. 這些自然生態(tài)因素導(dǎo)致艾比湖湖區(qū)及周邊地表水體復(fù)氧能力較差，水體中耗氧有機(jī)物和營(yíng)養(yǎng)物質(zhì)過(guò)剩，使該區(qū)域污染較為嚴(yán)重.

2)人為因素：根據(jù)吉仁塔[50]的研究可知：由于社會(huì)經(jīng)濟(jì)的發(fā)展，博爾塔拉河和精河流域內(nèi)人口增加，耕地面積增加，作物施肥量和灌溉增加，流域內(nèi)大量興建水庫(kù)，城市供水量和排水量增加，博爾塔拉河和精河兩大支流流經(jīng)的博樂(lè)市、溫泉縣和精河縣3個(gè)地區(qū)陸續(xù)出現(xiàn)廢水和排污的高污染企業(yè)，使得該區(qū)域比流域的其他區(qū)域污染更嚴(yán)重. 根據(jù)博州、伊犁州和塔城地區(qū)環(huán)境保護(hù)局官網(wǎng)公布的廢水排污超標(biāo)企業(yè)可知，豐水期和枯水期的B組(主要是博爾塔拉河和精河)區(qū)域內(nèi)的博樂(lè)市自來(lái)水公司污水處理廠、中糧屯河博州糖業(yè)有限公司和新疆騰博熱力有限公司等企業(yè)的排污和廢水排放導(dǎo)致該區(qū)域地表水體中生物污染、有機(jī)物污染和營(yíng)養(yǎng)物質(zhì)過(guò)剩；本區(qū)域的博爾塔拉蒙古自治州西部礦業(yè)有限公司、博爾塔拉蒙古自治州西域福樂(lè)礦業(yè)有限責(zé)任公司和精河縣新石石灰有限公司等企業(yè)的廢水廢料排放導(dǎo)致了水體中金屬物質(zhì)和懸浮物超標(biāo)，引起污染. 豐水期和枯水期A組(主要是奎屯河、古爾圖河和四棵樹河)區(qū)域內(nèi)的烏蘇市天然番茄制品有限公司和新疆天玉生物科技有限公司等企業(yè)的廢水、廢料排放導(dǎo)致本區(qū)域的地表水體中生物、化學(xué)污染嚴(yán)重. 豐水期的5月正是農(nóng)耕時(shí)期，農(nóng)民大量的施肥噴藥、灌溉和排水等農(nóng)業(yè)活動(dòng)造成水體中N、P濃度增高，COD、BOD5和NH3-N濃度等水質(zhì)指標(biāo)較高，使得豐水期較枯水期污染嚴(yán)重.

總之，本文利用地統(tǒng)計(jì)學(xué)方法及多元統(tǒng)計(jì)技術(shù)對(duì)艾比湖流域地表水水質(zhì)污染進(jìn)行時(shí)空分異特征研究和污染源解析，該方法雖然可以最大程度地避免信息丟失，提高分析結(jié)果的可靠性，但它是一種經(jīng)驗(yàn)識(shí)別，與SWAT、GWLF和HSPF等水質(zhì)模型方法[51]相比，缺少對(duì)水質(zhì)污染的量化預(yù)測(cè)，在污染源解析方面有一定局限性，所以今后考慮使用專門模型(如SPARROW模型)對(duì)水質(zhì)污染進(jìn)行量化分析和預(yù)測(cè).

[1] Zhang Yitao, Liu Hongbin, Wang Hongyuanetal. A bibliometric analysis of status and trend of international research on field nitrogen application effects on nitrogen losses and water quality.ActaEcologicaSinica, 2016, 36(20): 1-13. [張亦濤, 劉宏斌, 王洪媛等. 農(nóng)田施氮對(duì)水質(zhì)和氮素流失的影響. 生態(tài)學(xué)報(bào), 2016, 36(20): 1-13.]

[2] Li Yilu, Zhang Yuhu, Jia Haifengetal. Spatio-temporal characteristics and source identification of water pollutants in ancient town of Suzhou.ActaScientiaeCircumstantiae, 2014, (4): 1032-1044. [李義祿, 張玉虎, 賈海峰等. 蘇州古城區(qū)水體污染時(shí)空分異特征及污染源解析. 環(huán)境科學(xué)學(xué)報(bào), 2014, (4): 1032-1044.]

[3] Zhang Yuezhen, Dong Pingguo, Li Yanrongetal. Spatial variation regulation of water quality for centralized drinking water source in Wuwei City.JournalofDrainageandIrrigationMachineryEngineering(JDIME), 2015, (12): 1078-1084, 1104. [張?jiān)抡? 董平國(guó), 李彥榮等. 武威市集中式飲用水源地水質(zhì)空間分異規(guī)律. 排灌機(jī)械工程學(xué)報(bào), 2015, (12): 1078-1084, 1104.]

[4] Su S, Zhi J, Lou Letal. Spatio-temporal patterns and source apportionment of pollution in Qiantang River (China) using neural-based modeling and multivariate statistical techniques.PhysicsandChemistryoftheEarth,PartsA/B/C, 2011, 36(9): 379-386.

[5] Xu HS, Xu ZX, Wu Wetal. Assessment and spatiotemporal variation analysis of water quality in the Zhangweinan River Basin, China.ProcediaEnvironmentalSciences, 2012, 13(3): 1641-1652.

[6] Wang X, Cai Q, Ye Letal. Evaluation of spatial and temporal variation in stream water quality by multivariate statistical techniques: A case study of the Xiangxi River basin, China.QuaternaryInternational, 2012, 282: 137-144.

[7] Terrado M, Barceló D, Tauler Retal. Surface-water-quality indices for the analysis of data generated by automated sampling networks.TrACTrendsinAnalyticalChemistry, 2010, 29(1): 40-52.

[8] Akkoyunlu A, Akiner ME. Pollution evaluation in streams using water quality indices: A case study from Turkey’s Sapanca Lake Basin.EcologicalIndicators, 2012, 18: 501-511.

[9] Espejo L, Kretschmer N, Oyarzún Jetal. Application of water quality indices and analysis of the surface water quality monitoring network in semiarid North-Central Chile.EnvironmentalMonitoringandAssessment, 2012, 184(9): 5571-5588.

[10] Hu Kaiming, Li Bing, Wang Shuietal. Spatial distribution characteristics of water quality pollution in the Lake Taihu basin, Jiangsu Province.JLakeSci, 2014, 26(2): 200-206. DOI: 10.18307/2014.0205. [胡開明, 李冰, 王水等. 太湖流域(江蘇省)水質(zhì)污染空間特征. 湖泊科學(xué), 2014, 26(2): 200-206.]

[11] Li Yayao, Guo Wenjing, Dong Wenetal. Water environmental characteristics and pollution constitution for the Fenghe River in Shaanxi Province.ChinaRuralWaterandHydropower, 2014, (6): 70-76. [李亞嬌, 郭雯婧, 董雯等. 陜西灃河水環(huán)境特征及污染構(gòu)成解析. 中國(guó)農(nóng)村水利水電, 2014, (6): 70-76.]

[12] Wang Cuiyu, Yang Yonghui, Zhou Fengetal. Spatio-temporal characteristics and source identification of water pollutants in River Qinhe Basin.ActaScientiaeCircumstantiae, 2012, (9): 2267-2278. [王翠榆, 楊永輝, 周豐等. 沁河流域水體污染物時(shí)空分異特征及潛在污染源識(shí)別. 環(huán)境科學(xué)學(xué)報(bào), 2012, (9): 2267-2278.]

[13] Wang Lijing, Wang Xing, Liu Lusanetal. Multivariate analysis of water factors in Dongting Lake.ResearchofEnvironmentalSciences, 2013, (1): 1-7. [王麗婧, 汪星, 劉錄三等. 洞庭湖水質(zhì)因子的多元分析. 環(huán)境科學(xué)研究, 2013, (1): 1-7.]

[14] Kowalkowski T, Zbytniewski R, Szpejna Jetal. Application of chemometrics in river water classification.WaterResearch, 2006, 40(4): 744-752.

[15] Ma Xiaoxue, Wang Lachun, Liao Lingling. Spatio-temporal characteristics and source identification of water pollutants in Wenruitang River Watersheld.EnvironmentalScience, 2015, (1): 64-71. [馬小雪, 王臘春, 廖玲玲. 溫瑞塘河流域水體污染時(shí)空分異特征及污染源識(shí)別. 環(huán)境科學(xué), 2015, (1): 64-71.]

[16] Lam QD, Schmalz B, Fohrer N. Assessing the spatial and temporal variations of water quality in lowland areas, Northern Germany.JournalofHydrology, 2012, 438: 137-147.

[17] Bu H, Tan X, Li Setal. Temporal and spatial variations of water quality in the Jinshui River of the South Qinling Mts., China.EcotoxicologyandEnvironmentalSafety, 2010, 73(5): 907-913.

[18] Yang YH, Zhou F, Guo HCetal. Analysis of spatial and temporal water pollution patterns in Lake Dianchi using multivariate statistical methods.EnvironmentalMonitoringandAssessment, 2010, 170(1/2/3/4): 407-416.

[19] Zhang Wangchou, Li Xiaoxiu, Wang Xiaoyanetal. Temporal and spatial variations of water pollution in Wuqing section of Beiyunhe River.ActaScientiaeCircumstantiae, 2012, (4): 836-846. [張汪壽, 李曉秀, 王曉燕等. 北運(yùn)河武清段水污染時(shí)空變異特征. 環(huán)境科學(xué)學(xué)報(bào), 2012, (4): 836-846.]

[20] Li Wenzan, Li Xuyong, Wang Huiliangetal. Spatial distribution of the main contaminations in aquatic environment in Fuyang River.ActaScientiaeCircumstantiae, 2012, 11: 2814-2819. [李文贊, 李敘勇, 王慧亮等. 滏陽(yáng)河主要水環(huán)境污染物空間分布特性研究. 環(huán)境科學(xué)學(xué)報(bào), 2012, 11: 2814-2819.]

[21] Xie Sen, He Liansheng, Tian Xuedaetal. Study on temporal and spatial distribution patterns of water quality in Chaohu Lake.ChineseJournalofEnvironmentalEngineering, 2010, (3): 531-539. [謝森, 何連生, 田學(xué)達(dá)等. 巢湖水質(zhì)時(shí)空分布模式研究. 環(huán)境工程學(xué)報(bào), 2010, (3): 531-539.]

[22] Akbal F, Gürel L, Bahadr Tetal. Multivariate statistical techniques for the assessment of surface water quality at the mid-Black Sea coast of Turkey.Water,Air&SoilPollution, 2011, 216(1/2/3/4): 21-37.

[23] Singh EJK, Gupta A, Singh NR. Groundwater quality in Imphal West district, Manipur, India, with multivariate statistical analysis of data.EnvironmentalScienceandPollutionResearch, 2013, 20(4): 2421-2434.

[24] Mohammadpour R, Shaharuddin S, Zakaria NAetal. Prediction of water quality index in free surface constructed wetlands.EnvironmentalEarthSciences, 2016, 75(2): 1-12.

[25] Zhou Mei, Zhang Fei, Tashpolat Tiyipetal. Spatial distribution of hydro-chemistry characteristics and water quality in Ebinur Lake Basin.EnvironmentalChemistry, 2015, (9): 1694-1702. [周梅, 張飛, 塔西甫拉提·特依拜等. 艾比湖流域水化學(xué)指標(biāo)平面分布特征. 環(huán)境化學(xué), 2015, (9): 1694-1702.]

[26] Zhang Fei, Wang Juan, Tashpolat Tiyipetal. The spatial and temporal dynamic changes and driving forces in the surface area of Ebinur Lake from 1998-2013.ActaEcologicaSinica, 2015, (9): 2848-2859. [張飛, 王娟, 塔西甫拉提·特依拜等. 1998-2013年新疆艾比湖湖面時(shí)空動(dòng)態(tài)變化及其驅(qū)動(dòng)機(jī)制. 生態(tài)學(xué)報(bào), 2015, (9): 2848-2859.]

[27] Mi Yan, Chang Shunli, Shi Qingdongetal. Study on the effect of agricultural non-point source pollution to water environment of the Ebinur Lake Basin during high flow period.ArdZoneResearch, 2010, (2): 278-283. [彌艷, 常順利, 師慶東等. 農(nóng)業(yè)面源污染對(duì)豐水期艾比湖流域水環(huán)境的影響. 干旱區(qū)研究, 2010, (2): 278-283.]

[28] Yang Liu, Song Jianfei, Song Boetal. Primary pollutant water quality identification index method and its application to comprehensive evaluation of river water quality.EnvironmentalScience&Technology, 2015, (11): 239-245. [楊柳, 宋健飛, 宋波等. 主要污染物水質(zhì)標(biāo)識(shí)指數(shù)法在河流水質(zhì)評(píng)價(jià)的應(yīng)用. 環(huán)境科學(xué)與技術(shù), 2015, (11): 239-245.]

[29] Lattin J, Carroll D, Green P eds. Analyzing multivariate date. New York: Duxbury, 2003: 1-200.

[30] Chen Yanguang ed. Geography mathematical methods: Fundamentals and applications. Beijing: Science Press, 2011: 1-484. [陳彥光. 地理數(shù)學(xué)方法: 基礎(chǔ)和應(yīng)用. 北京: 科學(xué)出版社, 2011: 1-484.]

[31] Shrestha S, Kazama F. Assessment of surface water quality using multivariate statistical techniques: A case study of the Fuji river basin, Japan.EnvironmentalModelling&Software, 2007, 22(4): 464-475.

[32] Davis JC ed. Statistics date analysis in geology. New York: Wiley, 1973: 1-656.

[33] Huang Jinliang, Huang Yaling, Li Qingshengetal. Preliminary analysis of spatiotemporal variation of water quality and its influencing factors in the Jiulong River Watershed.EnvironmentalScience, 2012, (4): 1098-1107. [黃金良, 黃亞玲, 李青生等. 流域水質(zhì)時(shí)空分布特征及其影響因素初析. 環(huán)境科學(xué), 2012, (4): 1098-1107.]

[34] Zhou F, Guo H, Liu L. Quantitative identification and source apportionment of anthropogenic heavy metals in marine sediment of Hong Kong.EnvironmentalGeology, 2007, 53(2): 295-305.

[35] Qu Mingkai, Li Weidong, Zhang Chuanrongetal. Source apportionment of soil heavy metal Cd based on the combination of receptor model and geostatistics.ChinaEnvironmentalScience, 2013, (5): 854-860. [瞿明凱, 李衛(wèi)東, 張傳榮等. 基于受體模型和地統(tǒng)計(jì)學(xué)相結(jié)合的土壤鎘污染源解析. 中國(guó)環(huán)境科學(xué), 2013, (5): 854-860.]

[36] Mi Yan, Chang Shunli, Shi Qingdongetal. Aquatic environmental quality assessment in Ebinur Lake Catchment during high flow period, 2008.JLakeSci, 2009, 21(6): 891-894. DOI: 10.18307/2009.0621. [彌艷, 常順利, 師慶東等. 艾比湖流域2008年豐水期水環(huán)境質(zhì)量現(xiàn)狀評(píng)價(jià). 湖泊科學(xué), 2009, 21(6): 891-894.]

[37] Dai Guofei, Zhang Wei, Peng Ningyanetal. Study on distribution of N and P pollutants and risk of cyanobacteria bloom in Poyang Lake and waters around the lake during drought periods.EcologyandEnvironmentalSciences, 2015, 24(5): 838-844. [戴國(guó)飛, 張偉, 彭寧彥等. 枯水期鄱陽(yáng)湖及其濱湖水體氮磷等污染物分布與藻華風(fēng)險(xiǎn)研究. 生態(tài)環(huán)境學(xué)報(bào), 2015, 24(5): 838-844.]

[38] Dong Yu, Haimiti Yimiti, Bayin Chahanetal. Analysis on concentration and correlation of runoff and precipitation in BoErtala River basin.JournalofWaterResources&WaterEngineering, 2014, (4): 28-34. [董煜, 海米提·依米提, 巴音查汗等. 博爾塔拉河徑流與降水的集中性及關(guān)聯(lián)性分析. 水資源與水工程學(xué)報(bào), 2014, (4): 28-34.]

[39] Mi Yan. Investigation and assessment of agricultural non-point source pollution in Ebinur Lake Basin based on 3S [Dissertation]. Urumqi: Xinjiang University, 2010. [彌艷. 基于3S技術(shù)的艾比湖流域農(nóng)業(yè)非點(diǎn)源污染對(duì)水環(huán)境的影響研究[學(xué)位論文]. 烏魯木齊: 新疆大學(xué), 2010.]

[40] Yan Yixin, Jiang Xiaoxin, Ruan Xiaohuaetal. Water pollution characteristics and control measures in cities of plain river network area.WaterResourcesProtection, 2008, (5): 1-3, 29. [嚴(yán)以新, 蔣小欣, 阮曉紅等. 平原河網(wǎng)區(qū)城市水污染特征及控制對(duì)策研究. 水資源保護(hù), 2008, (5): 1-3, 29.]

[41] Fu Qing, Zheng Yan, Wang Shanjun. Water quality and spatial and temporal distribution characteristics of surface drinking water sources in China.JournalofBasicScienceandEngineering, 2015, (5): 886-894. [付青, 鄭彥, 王山軍. 全國(guó)地表飲用水源水質(zhì)狀況及時(shí)空分布特征. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào), 2015, (5): 886-894.]

[42] Bassam BFA, Fadhil AM. An analysis of groundwater in Sinjar Plain (Northwest of Iraq) using WQI model.JournalofChinaUniversityofGeosciences, 2001, 26(4): 291-295.

[43] Zhao Jie, Xu Zongxue, Liu Xingcaietal. Source apportionment in the Liao River Basin.ChinaEnvironmentalScience, 2013, (5): 838-842. [趙潔, 徐宗學(xué), 劉星才等. 遼河河流水體污染源解析. 中國(guó)環(huán)境科學(xué), 2013, (5): 838-842.]

[44] Chen Haiyang, Teng Yanguo, Wang Jinshengetal. Pollution load and source apportionment of non-point source nitrogen and phosphorus in Jinjiang River watershed.TransactionsoftheChineseSocietyofAgriculturalEngineering, 2012, (5): 213-219. [陳海洋, 滕彥國(guó), 王金生等. 晉江流域非點(diǎn)源氮磷負(fù)荷及污染源解析. 農(nóng)業(yè)工程學(xué)報(bào), 2012, (5): 213-219.]

[45] Ran Guihua, Ge Jiwen, Miao Wenjieetal. Spatial-temporal differentiation of water quality in Gufu River of Three Gorges Reservoir.ActaEcologicaSinica, 2013, (17): 5385-5396. [冉桂花, 葛繼穩(wěn), 苗文杰等. 三峽庫(kù)區(qū)古夫河水質(zhì)時(shí)空分異特征. 生態(tài)學(xué)報(bào), 2013, (17): 5385-5396.]

[46] Zeng Qingjing. Recent water volume flowing into Ebinur Lake and its change.AridLandGeography, 1995, (2): 17-23. [曾慶江. 艾比湖近期入湖水量及其變化. 干旱區(qū)地理, 1995, (2): 17-23.]

[47] Wu Chaocun. On the change law of hydrological characteristics of Aibi Lake Basin.EnergyandEnergyConservation, 2015, (5): 101-102, 115. [吳超存. 艾比湖流域的水文特征變化規(guī)律. 能源與節(jié)能, 2015, (5): 101-102, 115.]

[48] Dong Yu, Haimiti Yimiti. Runoff characteristics and responses to precipitation in the Ebinur Lake Basin: A case study of Boertala River.ResearchofSoilandWaterConservation, 2014, (2): 94-99. [董煜, 海米提·依米提. 艾比湖流域徑流水文特征及其對(duì)降水變化響應(yīng)——以博爾塔拉河為例. 水土保持研究, 2014, (2): 94-99.]

[49] Wang Juan, Zhang Fei, Zhang Yueetal. Correlation between the spatial water quality and land use/cover in the Ebinur Lake area.ActaEcologicaSinica, 2016, 36(24): 7971-7980. DOI: 10.5846/stxb201505271059. [王娟, 張飛, 張?jiān)碌? 艾比湖區(qū)域水質(zhì)空間分布特征及其與土地利用/覆被類型的關(guān)系. 生態(tài)學(xué)報(bào), 2016, 36(24): 7971-7980.]

[50] Ji Renta. The Ebinur Lake basin ecological environment problems and Countermeasures.ForestryofXinjiang, 2016, (1): 22-24. [吉仁塔. 艾比湖流域生態(tài)環(huán)境問(wèn)題與治理措施. 新疆林業(yè), 2016, (1): 22-24.]

[51] Wu Zaixing, Wang Xiaoyan. SPARROW-Watershed model and research progress.EnvironmentalScience&Technology, 2010, (9): 87-90, 139. [吳在興, 王曉燕. 流域空間統(tǒng)計(jì)模型SPARROW及其研究進(jìn)展. 環(huán)境科學(xué)與技術(shù), 2010, (9): 87-90, 139.]

Spatio-temporal characteristics and source identification of surface water pollutants in Lake Ebinur Watershed, Xinjiang

REN Yan1,2,3, ZHANG Fei1,2,3**, WANG Juan1,2,3, ZHANG Yue1,2,3, YANG Shengtian1,4, ABDUWASIT Ghulam5& KUNG Hsiangte6

(1:CollegeofResourcesandEnvironmentScience,XinjiangUniversity,Urumqi830046,P.R.China)(2:KeyLaboratoryofOasisEcology,XinjiangUniversity,Urumqi830046,P.R.China)(3:KeyLaboratoryofXinjiangWisdomCityandEnvironmentModeling,Urumqi830046,P.R.China)(4:CollegeofGeographyandRemoteSensingScience,BeijingNormalUniversity,Beijing100875,P.R.China)(5:CenterforSustainability,SaintLouisUniversity,St.Louis63108,USA)(6:DepartmentofEarthSciences,UniversityofMemphis,Memphis38152,USA)

Spatio-temporal characteristics and source identification of water pollutants are important for the management and protection of water resources. In this study, water-quality index (WQI) modeling, hierarchical cluster analysis, discriminant analysis, factor analysis and principle analysis of multiply linear regression (APCS-MLR) were used to analyze the water quality data sets from river network during 2015, including 10 parameters monitored at 54 different sites in Lake Ebinur Watershed of Xinjiang. Using GIS to study the spatio-temporal characteristics and source identification of surface water pollutants were the major goals of this study. The results showed that: The value ofWQIvaries between 38 and 70 during wet season when the water levels were high, while the value was between 31 and 71 during dry season. The Lake Ebinur water quality deterioration was more significant during wet season, and the water pollution status at Lake Ebinur, Bortala River and Jinghe River near Lake Ebinur were more serious than other tributaries. By cluster analysis and discriminant analysis, the basin could be divided into two groups. The group A included the eastern area of Lake Ebinur, Kuitun River, Guertu River and Sikeshu River. The group B included the western area of Lake Ebinur, Jinghe River and Bortala River. Water quality at the B group of Lake Ebinur Watershed was mainly influenced by chemical oxygen demand, dissolved oxygen, ammonia nitrogen and suspended solids concentrations during the dry season, more serious than the group A because there are some high polluting enterprises at the group B region, while the water environment governance mainly cancentrated in the B region. In Lake Ebinur, Jinghe River and Bortala River, water pollutants mainly came from organic and nutrient pollution, followed by mining and industrial pollution. In Kuitun River, Guertu Rive and Sikeshu River, surface water pollutants mainly came from organic and biological pollution, followed by nutrient pollution, while biological pollution effects were relatively weak. Results of the water quality study would provide valuable and useful information for the future management of water resources in Lake Ebinur.

Lake Ebinur Watershed; multivariate statistical methods; water quality index; spatio-temporal distribution; source identification

國(guó)家自然科學(xué)基金項(xiàng)目(41361045)和國(guó)家自然科學(xué)基金-新疆聯(lián)合基金項(xiàng)目(U1503302)聯(lián)合資助. 2016-07-10收稿；2016-10-31收修改稿. 任巖(1989～)，女，碩士研究生；E-mail: arielrenyan@163.com.

； E-mail: zhangfei3s@163.com.

DOI 10.18307/2017.0512