田皓予，佟玲，余國安，Arika Bridhikitti

（1.中國農(nóng)業(yè)大學水利與土木工程學院，北京100083；2.中國科學院地理科學與資源研究所陸地水循環(huán)及地表過程重點實驗室，北京100101；3.瑪希隆大學北碧府校區(qū)跨學科研究學院環(huán)境工程與災害管理系，泰國71150）

土地利用是人類出于經(jīng)濟、居住、娛樂、保護和政府目的等對陸地空間的利用[1]。土地利用常改變流域水文過程和養(yǎng)分遷移，有時直接形成非點源污染，進而引起河流水質(zhì)退化[2-3]，影響流域生態(tài)環(huán)境[4]。因此，分析流域水質(zhì)與土地利用的響應關系對于認識水環(huán)境變化成因和趨勢，構(gòu)建適宜的土地利用格局，促進流域水資源保護和生態(tài)環(huán)境修復有重要的科學價值和實踐意義。

對土地利用的分析通?；谀骋豢臻g尺度開展。近年來，國內(nèi)外學者基于不同空間尺度分析河流水質(zhì)與土地利用的關系。如Rodrigues等[5]基于全流域尺度分析巴西Córregoágua Limpa流域土地利用對河流水質(zhì)的影響。Shukla等[6]從流域和行政區(qū)兩個空間尺度分析了印度Ganga河流域水質(zhì)和土地利用之間的關系。Xiao等[7]從100～2 000 m站點緩沖區(qū)選擇4個尺度分析太湖流域水質(zhì)與土地利用的相關關系。Chen等[8]開展了從100 m到2 000 m緩沖區(qū)共11種不同空間尺度下北京市水質(zhì)與土地利用關系的研究。由于土地利用的多尺度格局差異，基于土地利用方式的河流水質(zhì)研究存在不確定性[9]，其結(jié)果（相關性）受空間分析尺度大小的影響，但哪個尺度的土地利用格局能夠更好地解釋與水質(zhì)之間的關系存在不同的結(jié)論[10]，有必要考慮從不同空間尺度開展分析。Ding等[11]分析東江河流域不同尺度下水質(zhì)與土地利用的關系后發(fā)現(xiàn)，流域尺度比緩沖區(qū)尺度能更好地解釋土地利用對河流水質(zhì)的影響；而Xu等[12]對淮河流域7種不同空間尺度下水質(zhì)和土地利用的關系研究發(fā)現(xiàn)，半徑小于20 km的緩沖區(qū)能夠更好地解釋土地利用對氨氮和溶解氧濃度變化的影響。可以看出，在分析河流水質(zhì)與土地利用關系時，不同流域適宜的空間分析尺度存在差異。因此，有必要對不同空間尺度進行對比，以尋找水質(zhì)與土地利用二者相關性更顯著的分析尺度。

河流水質(zhì)退化是泰國面臨的重要生態(tài)環(huán)境問題，其三分之一的河流處于退化或被污染狀態(tài)[13-14]。農(nóng)業(yè)營養(yǎng)富集、生活污水和工業(yè)廢水排放導致河流水體出現(xiàn)藻華現(xiàn)象[15-16]。在泰國的25個流域中有9個處于水資源短缺狀態(tài)，位于東北部的蒙河流域就是其中之一。作為湄公河最大支流，蒙河每年向湄公河注入約250億m3的水量[17]，其水質(zhì)對湄公河中下游水環(huán)境具有重要影響。蒙河流域是泰國重要的農(nóng)業(yè)區(qū)，也是泰國人口第二多的流域，聚集大量農(nóng)業(yè)人口[15]。近幾十年來由于經(jīng)濟發(fā)展和高密度人類活動影響，土地開發(fā)利用加劇，由此帶來的非點源污染問題愈顯突出，河流水環(huán)境也因此受到擾動[18]。

基于以上背景，本文以蒙河流域為研究對象，在之前研究[19]的基礎上應用莫蘭指數(shù)（Moran′s I）等空間分析方法和Spearman相關分析，探究小流域和不同緩沖區(qū)尺度下土地利用對河流水質(zhì)的影響，以期促進流域可持續(xù)發(fā)展，為平衡土地利用和水環(huán)境保護提供參考。

1 材料與方法

1.1 研究區(qū)概況

蒙河（Mun River）位于泰國東北部，大致位于14°～16°N，101°～105°E之間。蒙河發(fā)源于呵叻府附近的考艾國家公園，流經(jīng)武里南、素林、四色菊、烏汶等十府，最終匯入湄公河。蒙河干流長約800 km，流域面積約8.2萬km2（圖1）。流域氣候類型屬于濕潤熱帶氣候，常年溫度不低于18℃，平均年降雨量1300～1500 mm，受亞熱帶季風影響，5月中旬至10月中旬降雨頻繁，一般稱之為雨季，其中強降雨一般出現(xiàn)在8月或9月。11月至次年4月稱之為旱季，旱季降雨稀少，與雨季降雨量差異極大[20]。流域內(nèi)地形復雜，以平原為主，并分布有高原和山地。主要土壤類型為沙質(zhì)壤土和壤質(zhì)砂土等。研究區(qū)主要種植作物為水稻、玉米、木薯、甘蔗和橡膠林，有少量的養(yǎng)殖業(yè)。

圖1泰國蒙河流域位置、水系及水質(zhì)監(jiān)測站點示意圖Figure 1 Location of Mun River basin，drainage system and water quality monitoring sites

1.2 數(shù)據(jù)來源

蒙河水質(zhì)數(shù)據(jù)來自泰國污染管制署（Pollution Control Department，Thailand，簡稱PCD），數(shù)據(jù)序列為1997—2017年。為反映河流水質(zhì)的季節(jié)性變化，PCD一般在每年2、5、8、11月進行監(jiān)測。由下游到上游分布有MU01～MU20（蒙河的水質(zhì)監(jiān)測站點簡稱MU）共21個水質(zhì)監(jiān)測站點（圖1，包括MU17.1站點），其中MU19和MU20站點的監(jiān)測數(shù)據(jù)較少，因此選擇MU01～MU18共19個站點數(shù)據(jù)進行分析。從物理、化學、生物學3大類指標中選擇8項典型水質(zhì)指標開展分析，包括溶解氧（DO）、生化需氧量（BOD）、總大腸菌群（TCB）、糞大腸菌群（FCB）、總磷（TP）、硝態(tài)氮（NO-3-N）、氨氮（NH3-N）和固體懸浮物（SS）。由于蒙河流域為農(nóng)業(yè)區(qū)，基本無采礦區(qū)且工業(yè)生產(chǎn)水平較低，有毒重金屬（如Hg、Pb、As、Cd、Cr）含量幾乎全部低于泰國地表水質(zhì)量標準的三類水質(zhì)限值[21]，因此本文暫不考慮有毒重金屬指標。

土地利用數(shù)據(jù)來源于泰國土壤資源及土地發(fā)展部門獲取的2000、2007年和2015年三期矢量數(shù)據(jù)。元數(shù)據(jù)中有城鎮(zhèn)建設用地、林地、農(nóng)業(yè)用地、水域、混合用地5個一級分類，共28個二級分類，根據(jù)元數(shù)據(jù)中的二級分類將土地利用類別進行重新合并與劃分，分為水田、旱地、其他農(nóng)業(yè)用地、林地、城鎮(zhèn)建設用地、水域、混合用地共7個類別（圖2）。

1.3 分析方法

1.3.1 空間自相關分析

空間自相關在自然界普遍存在，地理信息系統(tǒng)使得探索自然變量的空間自相關性十分便利，如采用空間自相關指標，基于ArcGIS 10.2平臺進行水質(zhì)序列空間自相關指標的計算?？臻g自相關是同一自然變量在不同空間位置上的相關關系，空間自相關指標能反映水環(huán)境中污染物的分布狀況和各地區(qū)所含污染物情況的相關關系[22]。全局指標采用單一值來反映區(qū)域環(huán)境變量的空間自相關程度，用于分析整個研究區(qū)域的空間相關性模式。本文采用全局莫蘭指數(shù)（Global Moran′s I）來診斷流域?qū)崪y序列的總體空間分布結(jié)構(gòu)，計算公式為：

式中：I為全局自相關度量指標；zˉ為序列的平均值；zi、zj分別為水質(zhì)序列在站點i、j的數(shù)值(i≠j)；Wij為zi與zj之間的距離權(quán)重；N為研究的空間區(qū)域數(shù)。

莫蘭指數(shù)是一個有理數(shù)，經(jīng)方差歸一化后，其值歸一化到-1～1之間：I＞0表示空間正相關性，即水質(zhì)指標濃度高（低）值在空間上趨于集聚；I＜0表示空間負相關性，即某站點與周邊站點水質(zhì)指標濃度存在空間差異；I=0，說明不存在空間相關性，即流域各站點的水質(zhì)檢測序列在空間上呈隨機分布[23]。此外，檢測的空間自相關性還需進行顯著性檢驗，標準化的顯著性統(tǒng)計指標計算如公式（2），一般情況下，取顯著性水平為0.01和0.05，相應的標準方差Z值為2.58和1.96。

式中：Var（I）為方差；E（I）為期望。

1.3.2 不同空間尺度的建立

基于ArcGIS 10.2平臺的水文模塊，根據(jù)流域DEM數(shù)據(jù)和19個水質(zhì)監(jiān)測站點將整個流域劃分為19個小流域（圖3a），將某一站點與相鄰上游站點之間的區(qū)域稱為小流域（如MU09與MU08之間的區(qū)域為MU08站對應的小流域），從總體上分析每個水質(zhì)站點對應其上游小流域的土地利用對出水口水質(zhì)的影響。

圖2蒙河流域土地利用類別分布Figure 2 Land use in the Mun River basin

此外，本文還選擇3個緩沖區(qū)尺度進行分析。不同地形條件下緩沖區(qū)劃分方法不同：在坡度較大的流域，多以監(jiān)測站點上溯并向河道兩岸拓展一定寬度建立帶狀緩沖區(qū)[24-26]，而對于河網(wǎng)較為密集的平原流域，應以監(jiān)測站點為圓心，不同距離為半徑建立圓形緩沖區(qū)[27-28]?？紤]到蒙河流域大多為平原地區(qū)，且流域內(nèi)水網(wǎng)交織，故采用圓形緩沖區(qū)法。因此，本文選取小流域和5、2、1 km圓形緩沖區(qū)共4個空間尺度（圖3b和圖3c），并與土地利用圖疊加得到不同空間尺度下不同土地利用類型的面積占比。

圖3不同空間尺度示意（以MU19站點為例）Figure 3 Picture of different spatial scales（Take the site MU19 as an example）

1.3.3 統(tǒng)計分析

基于SPSS23.0軟件利用Spearman相關性分析研究2015年不同季節(jié)、不同空間尺度下水質(zhì)指標與土地利用類型面積占比之間的相關關系，使用雙尾檢驗在α=0.01和α=0.05顯著性水平下計算Spearman相關系數(shù)。

2 結(jié)果與分析

2.1 水質(zhì)指標分析

從1997—2017年19個水質(zhì)監(jiān)測站點水質(zhì)指標的統(tǒng)計分析結(jié)果（表1）和各水質(zhì)指標多年平均濃度的空間分布（圖4）可以發(fā)現(xiàn)，BOD多年平均濃度呈現(xiàn)出上游高、下游低的特點，DO多年平均濃度自上游向下游呈增大趨勢，不同站點間差異相對較小。TCB、FCB濃度分布特征大體相同，不同站點間差異極大，變異系數(shù)分別為3.39和6.02。TP、SS濃度呈現(xiàn)出上游高、下游低的特征，上游TP濃度高于《地表水環(huán)境質(zhì)量標準》（GB 3838—2002）中Ⅲ類水TP限值0.20 mg·L-1，是上游水體的主要污染物。NH3-N濃度的空間分布特征與TP相反，NO-3-N濃度上游和下游較高，中游較低。

為探究水質(zhì)時空變化特征且為避免時段過長掩蓋水質(zhì)指標的空間自相關性，將研究時段劃分為A（1997—2007年）和B（2008—2017年）兩個時段進行水質(zhì)指標的空間自相關檢測。根據(jù)公式（1）計算各時段各指標多年平均全局Moran′s I值，如表2所示。從時段A到時段B，DO、SS的Moran′s I值分別由0.70、0.23變?yōu)?.81、0.54，在時段B都通過P＜0.01的顯著性水平檢驗，說明DO和SS的空間分布表現(xiàn)出顯著的聚集性特征且在時段B空間聚集程度有所上升；BOD序列的Moran′s I值從時段A的0.46變?yōu)闀r段B的0.32，其空間正相關由極顯著變?yōu)椴伙@著，說明流域BOD的空間分布趨于局部化，聚集性特征減弱。TPNH3-N序列由時段A的空間負相關變?yōu)闀r段B的空間正相關，說明流域水質(zhì)營養(yǎng)指標空間差異性越來越小，逐漸轉(zhuǎn)變?yōu)榫奂Ｊ?。TCB、FCB序列的Moran′s I值從時段A到時段B有所增大，說明其離散程度變小。

2.2 土地利用結(jié)構(gòu)分析

農(nóng)業(yè)用地（水田、旱地和其他農(nóng)業(yè)用地）是蒙河流域最主要的土地利用類型（圖2），占流域總面積的73%以上。林地次之，2000、2007、2015年農(nóng)業(yè)用地和林地面積之和分別占流域總面積的92.10%、87.24%、86.55%，第3位是城鎮(zhèn)建設用地，水域面積最少。結(jié)合2000、2007、2015年3期土地利用數(shù)據(jù)與GIS工具，提取并統(tǒng)計蒙河流域19個站點對應不同空間尺度的土地利用面積占比情況（圖5，其中每個站點相應的土地利用類型面積占比柱狀圖從左至右3列依次表征2000、2007、2015年）。大體上，與2000年相比，2015年水田、旱地和林地面積有所減少，其他農(nóng)業(yè)用地、水域和混合用地面積有所增加。在小流域尺度土地利用分類中，農(nóng)業(yè)用地所占比例最大，尤其是水田（MU01和MU02站點對應的小流域位于蒙河流域東側(cè)邊緣，該處林地較多；MU05與MU06距離很近，因此劃分出的MU05所對應的小流域面積很小，導致小流域內(nèi)只有城鎮(zhèn)建設用地、水域和混合用地）。

表1水質(zhì)指標濃度的統(tǒng)計特征Table 1 The statistical feather values of water quality indicators concentration

圖4 1997—2017年蒙河代表性水質(zhì)指標平均濃度空間變化Figure 4 Spatial variations in the average concentration of representative water quality indicators in the Mun River basin from 1997 to 2017

表2 1997—2017年水質(zhì)空間全局Moran′s I值Table 2 Moran′s Ivalue of water quality from 1997 to 2017

圖5基于不同空間尺度的土地利用類型面積占比時空變化Figure 5 Spatiotemporal variations in the proportion of land use types based on different spatial scales

統(tǒng)計2015年不同空間尺度下各土地利用類型平均面積占比情況如表3所示，可以看出各類土地利用類型面積占比在不同空間尺度下呈現(xiàn)明顯差異，且基于不同空間尺度提取得到的面積占比大小排序也不同。其中，水田和旱地：小流域＞5 km緩沖區(qū)＞2 km緩沖區(qū)＞1 km緩沖區(qū)；林地、水域和其他農(nóng)業(yè)用地：1 km緩沖區(qū)＞2 km緩沖區(qū)＞小流域＞5 km緩沖區(qū)；城鎮(zhèn)建設用地和混合用地：1 km緩沖區(qū)＞2 km緩沖區(qū)＞5 km緩沖區(qū)＞小流域。產(chǎn)生這種差異的主要原因在于：（1）小流域尺度較其他3個緩沖區(qū)尺度大得多，更接近整個流域?qū)嶋H的土地利用分布格局；（2）各水質(zhì)監(jiān)測站點主要位于交通便利、易于采樣的位置，一般距離城鎮(zhèn)、水域、道路和橋梁等較近，因此在緩沖區(qū)范圍內(nèi)的城鎮(zhèn)建設用地、水域、混合用地等占比更大，且緩沖區(qū)尺度越小，其占比越大。此外，在4種不同空間尺度下水田都為優(yōu)勢地類。

2.3 水質(zhì)指標與土地利用類型的相關分析

為進一步探究水質(zhì)與土地利用變化的響應關系，將2015年19個水質(zhì)站點監(jiān)測的水質(zhì)數(shù)據(jù)和對應的不同空間尺度下土地利用類型的面積占比進行關聯(lián)，利用SPSS 23.0軟件進行Spearman相關性分析，如表4所示。本研究只對5種主要土地利用類型進行分析，即水田、旱地、林地、城鎮(zhèn)建設用地和水域，其他農(nóng)業(yè)用地和混合用地中包含的植被類型種類較多，難以明確具體對水質(zhì)的影響，因此暫不進行分析。蒙河不同季節(jié)水質(zhì)參數(shù)濃度變化較大，因此分為汛期（6月至10月）和非汛期（11月至次年5月中旬）兩個階段分別對水質(zhì)指標平均濃度與不同空間尺度土地利用類型面積占比進行分析。

由表4可以看出，2015年蒙河流域水體DO含量與林地面積占比呈正相關，相關系數(shù)在5 km緩沖區(qū)時最高，為0.69（P＜0.01），與城鎮(zhèn)建設用地面積占比為負相關；BOD含量與農(nóng)業(yè)用地面積占比呈正相關，與林地面積占比呈負相關，且在5 km緩沖區(qū)尺度下相關性最優(yōu)，相關系數(shù)為-0.58；TCB和FCB含量與城鎮(zhèn)建設用地面積占比呈顯著正相關，且非汛期的相關關系更強；TP含量與農(nóng)業(yè)用地和小流域尺度下的城鎮(zhèn)建設用地面積占比呈正相關，而與林地、水域和緩沖區(qū)尺度（5、2、1 km）下的城鎮(zhèn)建設用地面積占比呈負相關，與林地的相關系數(shù)絕對值最大為0.72，出現(xiàn)在5 km緩沖區(qū)尺度內(nèi)；-N和NH3-N含量與農(nóng)業(yè)用地和小流域尺度下的城鎮(zhèn)建設用地面積占比呈負相關，與林地、水域和緩沖區(qū)尺度下的城鎮(zhèn)建設用地面積占比呈正相關。受降雨徑流影響，汛期-N和NH3-N含量與水田、林地面積占比的相關性優(yōu)于非汛期；SS含量與城鎮(zhèn)建設用地面積占比呈負相關，與農(nóng)業(yè)用地面積占比呈正相關，與林地和水域的關系表現(xiàn)出季節(jié)性差異，汛期與二者呈正相關，非汛期與二者呈負相關。

3 討論

3.1 河流水質(zhì)特征及土地利用類型組分的關系

分析顯示，在全流域尺度上，蒙河水體DO、BOD濃度具有顯著的空間正相關特征，呈聚集性分布，DO濃度流域下游高，BOD濃度流域上游高。這與其他受強人類活動（尤其是水利工程）影響的流域存在差異，如在我國淮河流域，水體DO濃度在空間上呈離散分布[29]，主要原因是淮河流域閘壩工程較多，天然河道被其分割，且很多排污口散亂地分布于河道兩側(cè)，而蒙河受水利工程影響較小。

許多研究發(fā)現(xiàn)，河流水體養(yǎng)分含量與流域農(nóng)業(yè)活動密切相關[8，30-31]。與其他區(qū)域的研究結(jié)果（如我國云南撫仙湖[32]）相似，蒙河流域水體營養(yǎng)指標（TP、NO-3-N、NH3-N）與農(nóng)業(yè)用地面積比例的相關性十分顯著，其中水體TP含量與農(nóng)業(yè)用地面積占比具有明顯的正相關關系，而-N和NH3-N含量與農(nóng)業(yè)用地占比存在明顯的負相關關系，這與蒙河流域農(nóng)業(yè)用地利用類型和主導作物種類對氮、磷的吸收能力有重要聯(lián)系。蒙河流域土地利用的主導類型為農(nóng)業(yè)用地，其中水田面積占流域總面積的50%以上[18]。水稻等農(nóng)作物生長過程中常大量施用氮、磷肥。相關研究表明水稻吸收氮的能力很強而吸收磷的能力相對較弱[33]，如超級雜交稻對氮素的吸收率為177.69～189.09 kg·hm-2，而對磷素的吸收率僅為36.94～39.80 kg·hm-2[34]。未被作物吸收的磷肥由于降雨濺蝕、水流入滲淋溶和徑流的沖刷進入河流，從而造成汛期水體磷含量顯著升高，而水稻根系有效吸收和土壤滯留等相關生物化學過程則顯著消減了地表徑流中的氮類營養(yǎng)物質(zhì)。水體-N和NH3-N含量較低還可能與氮素的遷移轉(zhuǎn)化有關，農(nóng)田中有大量的微生物，進而造成不同形態(tài)氮素之間的硝化、反硝化反應[5，35]。此外，流域上游的農(nóng)業(yè)用地面積占比較下游大，而NO-3-N和NH3-N的含量則在下游更高，所以農(nóng)業(yè)用地占比與河流氮含量之間呈負相關關系。

表3不同空間尺度下各土地利用類型面積平均占比（%）Table 3 Average proportion of land use types area at different spatial scales（%）

表4基于不同空間尺度的2015年水質(zhì)指標與土地利用類型的相關性Table 4 Correlation between water quality indicators and land use types in 2015 based on different spatial scales

不同空間尺度下FCB含量與城鎮(zhèn)建設用地面積占比呈正相關關系。緩沖區(qū)尺度下水體TCB含量也與城鎮(zhèn)建設用地呈正相關，與我國撫仙湖[32]和北運河流域[36]的研究結(jié)果類似。不同空間尺度上水體TCB和FCB含量與城鎮(zhèn)建設用地的相關系數(shù)均高于其他土地利用類型，說明作為人類活動密集區(qū)的城鎮(zhèn)對水體中細菌污染物的引入有重要貢獻。非汛期城鎮(zhèn)建設用地面積比例與細菌類指標的相關系數(shù)大于汛期。城鎮(zhèn)建設用地代表著高強度的人類活動，生活污水中含有大量的TCB、FCB等細菌類污染物以及氮類營養(yǎng)物質(zhì)。隨著城鎮(zhèn)和建設用地的擴張，居民生活垃圾和污水、工業(yè)廢水的排放相應增加，對河流水質(zhì)造成影響。

3.2 水質(zhì)與土地利用關系的尺度效應

空間尺度效應是水質(zhì)與土地利用關系分析中的一個重要問題，但對于最佳空間分析尺度范圍仍存在爭議[8]。有研究認為，小流域尺度下的河流水質(zhì)與土地利用的關系較緩沖區(qū)尺度下更為顯著[37-38]，也有研究發(fā)現(xiàn)緩沖區(qū)尺度的土地利用格局能更好地解釋河流水質(zhì)變化[12，39]?；诿珊恿饔虻姆治鲲@示，大部分水質(zhì)指標與土地利用在小流域和不同緩沖區(qū)尺度上表現(xiàn)出一致的相關性，但相較于小流域尺度，緩沖區(qū)尺度（尤其5 km半徑緩沖區(qū)）下更多的水質(zhì)指標與土地利用類型占比之間的相關關系通過P＜0.05的顯著性水平檢驗，小流域及5、2、1 km緩沖區(qū)尺度下具有顯著相關關系的指標數(shù)量在汛期和非汛期分別為7、7、7、5項和5、12、6、6項。非汛期水質(zhì)與土地利用之間顯著相關的指標數(shù)量多于汛期，其主要原因是由于非汛期水質(zhì)與城鎮(zhèn)建設用地顯著相關的指標數(shù)量較汛期多，因為非汛期城鎮(zhèn)居民生活污水等點源污染對水質(zhì)的影響比汛期更大。因此推測，蒙河流域5 km緩沖區(qū)尺度下的土地利用能夠更好地解釋與水質(zhì)之間的關系，這與Yadav等[15]對蒙河流域的研究結(jié)果相似。各流域水質(zhì)與土地利用的相關關系存在不同的最佳空間分析尺度，這可能主要是受流域獨特的自然屬性影響[40]。

4 結(jié)論

（1）蒙河流域河流水質(zhì)空間自相關檢測結(jié)果顯示：DO、BOD、TP、NO-3-N、NH3-N、SS序列在流域站點間呈一般或顯著的空間正相關，而TCB、FCB主要來自人類和畜禽的廢物排放造成的點源污染，其濃度高值分布較為分散，表現(xiàn)出空間差異性。

（2）水質(zhì)與土地利用之間的相關關系在不同空間尺度存在一致性，但緩沖區(qū)尺度下的土地利用格局能夠更好地解釋蒙河水質(zhì)變化，其中5 km緩沖區(qū)范圍內(nèi)水質(zhì)指標與土地利用的關系更為顯著。

（3）蒙河水質(zhì)與流域農(nóng)業(yè)用地（尤其是水田）之間的關系最為顯著，其中BOD和TP含量與農(nóng)業(yè)用地面積占比呈顯著正相關，說明農(nóng)業(yè)用地對河流水體BOD和磷的污染負荷起到了“源”的作用。水體DO含量與城鎮(zhèn)建設用地面積占比呈負相關，而BOD、TCB和FCB含量與其呈正相關，說明城鎮(zhèn)建設用地對這3類水質(zhì)指標的負荷起到“源”的作用。相反，DO含量與林地呈正相關，BOD和TP含量與其呈負相關，說明林地對水體BOD和TP污染負荷起到截留、吸附的作用。

（4）蒙河流域內(nèi)主要土地利用類型為農(nóng)業(yè)用地，其面積占流域總面積的73%以上，河流水質(zhì)主要受農(nóng)業(yè)非點源污染的影響，還受城鎮(zhèn)化發(fā)展和人口密度增大帶來的廢水排放的影響。