張 淼，劉俊國，趙 旭，陳忞忞

（北京林業(yè)大學(xué) 自然保護(hù)區(qū)學(xué)院，北京100083）

濕地被譽(yù)為“地球之腎”，具備水質(zhì)調(diào)節(jié)、水源供給、洪水調(diào)蓄、生物資源、休閑游憩等多種生態(tài)系統(tǒng)服務(wù)功能。然而，濕地也是最易受人類活動(dòng)干擾的脆弱生態(tài)系統(tǒng)［1］。由于人類活動(dòng)導(dǎo)致的過度墾殖、不合理水資源利用、過量的污染物排放，使得全球濕地退化嚴(yán)重，濕地面積急劇減小，濕地生態(tài)環(huán)境日趨惡化［2］。近年來，我國政府通過濕地生態(tài)修復(fù)、退耕還湖、保護(hù)區(qū)建立等政策措施在一定程度上減緩了濕地退化的趨勢。但由于監(jiān)測條件有限，人類活動(dòng)對濕地的影響很難量化評價(jià)［3－4］。

為了對濕地進(jìn)行從簡單到復(fù)雜的全方位評價(jià)，20世紀(jì)90年代美國環(huán)保局（EPA）提出了濕地評價(jià)的“三級”評價(jià)體系：景觀尺度評價(jià)的Level 1方法、實(shí)地快速評價(jià)的Level 2方法、生物物理化學(xué)指標(biāo)評價(jià)的 Level 3方法［5］。2005年 Brown等［6］提出的景觀開發(fā)強(qiáng)度方法（Landscape Development Intensity，以下簡稱“LDI方法”），屬于EPA“三級”方法中的Level 1方法。LDI方法認(rèn)為，人類活動(dòng)對自然環(huán)境的干擾反映在其對自然生態(tài)系統(tǒng)人為的能量輸入上。對一種土地利用類型輸入的能量越大，表示人類對生態(tài)系統(tǒng)健康干擾越大，生態(tài)系統(tǒng)健康程度越低。該方法直接利用遙感數(shù)據(jù)，通過人類對土地利用的干擾情況評價(jià)濕地健康。由于簡單易行，需要數(shù)據(jù)量較小，并能為管理者快速提供生態(tài)系統(tǒng)健康狀況，已經(jīng)在美國佛羅里達(dá)州［7］、美國俄亥俄州［8］等地區(qū)得到廣泛應(yīng)用。但目前在中國的相關(guān)研究和應(yīng)用幾乎處于空白。

本文結(jié)合研究區(qū)土地開發(fā)現(xiàn)狀，將LDI方法應(yīng)用于我國干旱區(qū)典型湖泊濕地黃旗海的生態(tài)系統(tǒng)健康評價(jià)中，并通過考慮各類土地開發(fā)活動(dòng)引起的能源消耗隨時(shí)間變化，實(shí)現(xiàn)LDI系數(shù)的時(shí)間序列核算，從而克服原方法只適合于單一年份或者土地利用類型變化較小的濕地評價(jià)，無法對演變中的濕地進(jìn)行評價(jià)的問題。

1 研究區(qū)概況

黃旗海位于內(nèi)蒙古自治區(qū)烏蘭察布市察哈爾右翼前旗中南部，屬于內(nèi)陸閉塞湖泊，地理位置為東經(jīng)113°10′—113°23′，北緯 40°48′—40°55′，海拔高度1 266m。黃旗海處于典型的北方干旱地區(qū)，年平均氣溫4.7℃，年平均降水量369.7mm，而年平均潛在蒸發(fā)量卻高達(dá)1 957.4mm［9］。補(bǔ)給水源主要來自霸王河和泉玉林河等11條河溝［10］，河水只進(jìn)不出，礦物質(zhì)沉積，湖水呈堿性。黃旗海地表平坦，略向湖心傾斜，由全新世湖積淤泥質(zhì)黏土組成。洼地邊緣有沼澤、洼地中土壤鹽漬化較嚴(yán)重，地下水位較高，隱域性的低濕洼地上生長著芨芨草、寸草苔等雜類草為主的草甸植被，是野生濕地鳥類理想的棲息地。歷史上湖水最大面積近130km2，平均湖深2m左右，其中水深1m的湖面約占70%［11］。近年來，由于工農(nóng)業(yè)用水量不斷增加，流域內(nèi)興修水利，黃旗海湖面面積逐漸縮小，工業(yè)廢水不斷向湖內(nèi)排放，使湖水水質(zhì)惡化。2008年黃旗海徹底干涸，2012年少有恢復(fù)。

2 研究方法與數(shù)據(jù)

2.1 LDI系數(shù)核算

LDI方法利用遙感、GIS工具確定研究區(qū)土地利用空間分布，進(jìn)而計(jì)算空間上不同土地利用類型的LDI系數(shù)，并以各土地利用占總體土地面積百分比作為權(quán)重計(jì)算總體LDI值［6］：

式中：LDItotal——研究區(qū)總體LDI值；%LUi——第i種土地利用類型的面積占土地總面積的百分比；LDIi——第i種土地利用類型的LDI系數(shù)。

LDI系數(shù)的計(jì)算方法為：（1）調(diào)查某種土地利用類型因人類活動(dòng)而產(chǎn)生的能源消耗項(xiàng)目（除可更新能源之外）；（2）核算各能源消耗項(xiàng)目所產(chǎn)生的能值；（3）各土地利用類型消耗的能值取自然對數(shù)并標(biāo)準(zhǔn)化到1至10之間，即為LDI系數(shù)［6］。其中，1代表處于完全自然狀態(tài)下的土地，10代表被人類高度開發(fā)的土地［6］。研究區(qū)總體LDI值越大，表明人類干擾強(qiáng)度越強(qiáng)烈。

數(shù)據(jù)方面，以黃旗海湖區(qū)作為研究對象，對湖區(qū)周圍6個(gè)緩沖區(qū)寬度（100，200，300，400，500，600 m）進(jìn)行評價(jià)。根據(jù)Landsat MSS、TM、ETM影像解譯得到1976年、1987年、1993年、1999年、2006年、2010年的土地利用類型分類數(shù)據(jù)［12］。利用ArcGIS對遙感影像的結(jié)果進(jìn)行統(tǒng)計(jì)，得到黃旗海湖區(qū)空間柵格大小為30m×30m的主要土地利用類型（濕地、裸地、草地、以及耕地）。

2.2 各土地利用類型能耗

LDI方法的實(shí)質(zhì)是核算不同土地利用類型的能值輸入量（可更新能源除外），以LDI系數(shù)表示［13］。能值是生產(chǎn)產(chǎn)品或服務(wù)所需要的直接及間接能量，以太陽能焦耳作為統(tǒng)一的衡量單位［（sej／（hm2·a）］［14］。能值使用大自然的價(jià)值體系，基于有效能量的流動(dòng)，從環(huán)境、社會(huì)和經(jīng)濟(jì)等方面全面客觀地衡量商品及服務(wù)的價(jià)值，為全面有效客觀地評價(jià)生態(tài)系統(tǒng)健康提供統(tǒng)一的計(jì)量標(biāo)準(zhǔn)［15］。

各土地利用類型的能值輸入可歸為四種類型：源于自然系統(tǒng)的可更新能源，不可更新能源，來自于社會(huì)環(huán)境的物資材料和服務(wù)（圖1）［16］。以農(nóng)田為例，可更新能源包括太陽能、雨水化學(xué)能、雨水勢能、風(fēng)能等［17］；不可更新能源為凈表土損失能；物資材料為殺蟲劑、磷酸鹽、氮肥等，服務(wù)輸入例如勞工等［18］。LDI方法認(rèn)為除可更新能源之外，其它能值的輸入代表人類對生態(tài)系統(tǒng)健康的干擾［6］：能值越高，生態(tài)系統(tǒng)健康程度越低。

圖1 通用生產(chǎn)系統(tǒng)的能值系統(tǒng)關(guān)系

為了使LDI方法能夠適用于研究區(qū)，本研究對Brown等［6］的能耗核算方法進(jìn)行相應(yīng)調(diào)整，首先，本文加入了農(nóng)業(yè)非點(diǎn)源污染的能耗。農(nóng)業(yè)非點(diǎn)源污染是指農(nóng)田中的土粒、氮素、磷、農(nóng)藥及其它有機(jī)或無機(jī)污染物在降水或灌溉過程中，通過農(nóng)田地表徑流、農(nóng)田排水和地下滲漏，排入水體所造成的污染［19］。全國第一次污染源普查結(jié)果顯示，農(nóng)業(yè)源污染（大部分為非點(diǎn)源污染）對污染物總量（COD、氮、磷）的貢獻(xiàn)高達(dá)46%［20］?？紤]黃旗海農(nóng)業(yè)非點(diǎn)源污染的能耗是因?yàn)椋海?）黃旗海湖區(qū)污染較為嚴(yán)重，2011年的調(diào)查結(jié)果顯示，黃旗海水質(zhì)為劣Ⅴ類；（2）黃旗海周邊地區(qū)基本不存在工業(yè)，但農(nóng)田較多?？梢耘袛嗷室约稗r(nóng)藥的使用所產(chǎn)生的農(nóng)業(yè)非點(diǎn)源污染是黃旗海的主要污染源。本文采用的非點(diǎn)源污染計(jì)算公式如下［21］：

式中：NPS——非點(diǎn)源污染量；H——化肥（N、P、K）的施用量；10%——假設(shè)根區(qū)水分平衡后污染物流入水體的百分比［21］。本文假設(shè)黃旗海各項(xiàng)土地利用類型都存在非點(diǎn)源污染輸入。

其次，開發(fā)LDI方法的最初目的是為管理者提供一種快速評估地區(qū)生態(tài)系統(tǒng)健康狀況的工具［4］。Brown所提出的LDI方法，因?yàn)槠淠芎捻?xiàng)目固定不變，只能評價(jià)單一年份或者土地利用類型變化較小濕地的生態(tài)系統(tǒng)健康。本文考慮了能耗項(xiàng)目隨時(shí)間的變化，使結(jié)果能夠反映濕地生態(tài)系統(tǒng)健康的變化趨勢。能值核算需要的數(shù)據(jù)包括：能耗項(xiàng)目，輸入量以及能值轉(zhuǎn)換率，能耗項(xiàng)目根據(jù)現(xiàn)場調(diào)研結(jié)合文獻(xiàn)查閱獲得［22－23］，輸入量主要來自《內(nèi)蒙古統(tǒng)計(jì)年鑒》［24］，能值轉(zhuǎn)換率根據(jù)Brown結(jié)果［6］獲得。

3 結(jié)果與分析

3.1 各土地利用類型面積變化

各土地利用類型的面積及占總面積的百分比見表1?？梢钥闯?，1976年黃旗海湖區(qū)大部分面積以濕地為主（88%），伴有部分草地（12%），耕地和裸地面積為0。隨著時(shí)間的推移，濕地面積占比例逐漸縮小，到2006年達(dá)到最小值（19%），到2010年有所恢復(fù)（33%）。自1993年建立“黃旗海自然保護(hù)區(qū)管理站”以來，已經(jīng)禁止在湖區(qū)內(nèi)部放牧和耕作，但違禁放牧和違法耕作仍時(shí)有發(fā)生。表現(xiàn)在面積上，草地面積逐漸增大，到2010年其所占比例達(dá)到最大值45%。較大面積的草地也間接引發(fā)了湖區(qū)內(nèi)部違禁放牧的情況。湖區(qū)內(nèi)部耕作情況同樣沒有完全禁止，在2006年其所占比例達(dá)到最大值7%。值得注意的是，2006年和2010年湖區(qū)內(nèi)開始出現(xiàn)裸地。通過表1可以發(fā)現(xiàn)，由于草地與耕地的面積之和在這兩年幾乎沒有變化，可以判斷裸地是由濕地轉(zhuǎn)化而來。這一轉(zhuǎn)化在1999—2006年間較高（58%的濕地轉(zhuǎn)化為裸地），而2006—2010年又發(fā)生了逆向轉(zhuǎn)化（43%的裸地轉(zhuǎn)化為濕地）。將裸地轉(zhuǎn)化為濕地應(yīng)是今后濕地恢復(fù)的重點(diǎn)。

表1 黃旗海湖區(qū)土地利用類型及面積

3.2 各土地利用類型能耗

對各研究年各土地利用類型能耗及相應(yīng)的LDI系數(shù)隨時(shí)間變化進(jìn)行核算。表2顯示了2010年各土地利用類型能耗項(xiàng)目及其產(chǎn)生的單位面積能值。可以看出，各土地利用類型中，耕地單位面積能值最高，表明其受人類干擾程度最大。其次是草地與裸地，濕地最低。這是因?yàn)椋卦谒型恋乩妙愋椭心苤递斎腠?xiàng)目最多，且數(shù)值均較大。而農(nóng)藥和化肥的直接輸入是耕地總體能耗最大的主要原因。按照Brown的假設(shè)，濕地由于是自然系統(tǒng)并沒有能源輸入，其單位面積能值為0［6］。但Brown并沒有考慮農(nóng)業(yè)非點(diǎn)源等排污活動(dòng)對濕地的間接影響。通過計(jì)算農(nóng)業(yè)非點(diǎn)源污染，黃旗海湖區(qū)內(nèi)的濕地在2010年產(chǎn)生6.2×1014sej／（hm2·a）的能值消耗。對于草地與裸地來說，其單位面積能值的計(jì)算結(jié)果同樣由于考慮農(nóng)業(yè)非點(diǎn)源污染而增大。

表2 2010年各土地利用類型能值

根據(jù)各土地利用類型能值得到相應(yīng)的LDI系數(shù)（圖2）。LDI系數(shù)的大小排列與各土地利用類型能值大小一致。其中，濕地的LDI系數(shù)最低（1.00～1.77），耕地最高（4.65～5.51）。時(shí)間變化方面，1976—2010年間所有土地利用類型的LDI系數(shù)都呈現(xiàn)增加的趨勢。其中濕地增幅最大，LDI系數(shù)增加了77%。這表明濕地對人類干擾最為敏感，惡化程度最快。其次是裸地（31%）、耕地（18%）與草地（10%）。

圖2 各土地利用類型景觀開發(fā)強(qiáng)度系數(shù)值

3.3 LDI系數(shù)空間變化

選取黃旗海濕地及其周圍100～600m緩沖區(qū)，將每一個(gè)柵格當(dāng)作一個(gè)景觀單元，得到黃旗海湖區(qū)（包含600m緩沖帶）在1976—2010年間6a的LDI系數(shù)空間分布圖（圖3）。根據(jù)Brown的研究，LDI系數(shù)在1～2之間屬于自然類型濕地。而LDI系數(shù)為1屬于純自然濕地，表明濕地不受任何人類干擾，生態(tài)系統(tǒng)健康程度最佳［5］。從圖3可以看出，黃旗海濕地從內(nèi)向外LDI系數(shù)逐漸增大。邊緣LDI系數(shù)最大，且隨時(shí)間逐漸增加，2010年部分邊界LDI系數(shù)達(dá)到最大值5.51。而自然類型濕地主要分布在湖區(qū)內(nèi)部，其面積總體呈現(xiàn)減少的趨勢 。而純自然類型濕地的變化則顯示了人類干擾活動(dòng)從湖區(qū)邊緣向內(nèi)部的延伸：1976年黃旗海湖區(qū)近80%土地的LDI系數(shù)為1，整個(gè)湖區(qū)受人類干擾程度非常低。1976—1993年，純自然類型濕地的面積逐漸變少，從77%下降到42%。1999以后純自然類型濕地已經(jīng)消失。

3.4 黃旗海濕地及其緩沖區(qū)生態(tài)系統(tǒng)健康變化評價(jià)

圖4展示了湖區(qū)及6個(gè)緩沖帶的LDI值變化趨勢。黃旗海濕地的LDI值在包含不同距離緩沖區(qū)下變化趨勢基本一致，LDI值隨著緩沖區(qū)的擴(kuò)大而逐漸變大。緩沖區(qū)大小與LDI值呈顯著正相關(guān)（p＜0.001）。參考Brown對LDI值進(jìn)行的分類［6］，得出黃旗海濕地生態(tài)系統(tǒng)健康評價(jià)標(biāo)準(zhǔn)（表3）?？梢钥闯?，1976年和1987年6個(gè)緩沖區(qū)的LDI值均低于2.0，說明1976—1987年黃旗海湖區(qū)及其緩沖區(qū)處于生態(tài)系統(tǒng)健康的狀態(tài)。1987—1993年間LDI值呈現(xiàn)增長趨勢，1993年除湖區(qū)外，其它6個(gè)緩沖區(qū)LDI值均超過2.0，黃旗海濕地處于生態(tài)系統(tǒng)健康程度由健康轉(zhuǎn)為一般轉(zhuǎn)折點(diǎn)。1999年，湖區(qū)及其緩沖區(qū)均已超過2.0，但LDI值變化不大。從1999年開始，LDI值又呈急劇上升趨勢，接近于3.0；2006—2010年，LDI值有略下降趨勢。1999年以后黃旗海濕地健康程度為“一般”。黃旗海濕地生態(tài)系統(tǒng)健康程度總體上呈現(xiàn)逐漸惡化的趨勢。

圖3 黃旗海濕地（包含湖區(qū)以及周邊600m緩沖帶）的LDI空間分布

表3 LDI值評價(jià)濕地生態(tài)系統(tǒng)健康

圖4 六期不同緩沖區(qū)LDI值

4 結(jié)論

本文采用LDI方法，結(jié)合遙感與GIS工具，對黃旗海濕地及其緩沖區(qū)1976—2010年6a間的空間LDI系數(shù)以及總體LDI值進(jìn)行核算，并進(jìn)一步評價(jià)其生態(tài)系統(tǒng)健康變化。結(jié)果顯示，1993—1999年屬于生態(tài)系統(tǒng)健康轉(zhuǎn)變的過渡時(shí)期，在這一時(shí)期之前黃旗海濕地內(nèi)部大部分地區(qū)仍屬于純自然類型濕地，人類干擾活動(dòng)多發(fā)生在湖區(qū)邊緣，生態(tài)系統(tǒng)處于“健康”狀態(tài)。之后由于人類活動(dòng)的加劇，其影響從邊緣延伸到湖區(qū)內(nèi)部所有地區(qū)，使得1999—2010年間生態(tài)系統(tǒng)處于“一般”的狀態(tài)。由于黃旗海濕地位于季風(fēng)邊緣區(qū)，屬于氣候上的過渡地帶和生態(tài)系統(tǒng)上的脆弱地帶，對人類和環(huán)境變化的響應(yīng)十分敏感。因此本研究對研究氣候過渡帶典型濕地與人類活動(dòng)的響應(yīng)關(guān)系，濕地的管理與修復(fù)具有重要意義。

本文研究結(jié)果與實(shí)地考察調(diào)研結(jié)果相符合。黃旗海濕地LDI值1976—1993年呈上升趨勢，通過實(shí)地考察及走訪政府部門和周邊村民得知，1976—1993年黃旗海周邊村民放牧問題較嚴(yán)重，并且化肥廠和糖廠對黃旗海的水污染較嚴(yán)重；1993—1999年國家實(shí)施禁牧政策，并在1993年建立了黃旗海濕地自然保護(hù)區(qū)，黃旗海周邊工廠逐漸停止運(yùn)行，放牧問題和農(nóng)業(yè)污染問題都得到了一定的制止，這對人類干擾活動(dòng)產(chǎn)生了一定的抑制作用。1999—2006年，黃旗海周邊地區(qū)農(nóng)業(yè)發(fā)展迅速，人類活動(dòng)干擾加劇，2007年開始黃旗海湖泊逐漸干涸，成為季節(jié)性湖泊。此后，黃旗海自治區(qū)級自然保護(hù)區(qū)加大對黃旗海的保護(hù)力度，抑制濕地生態(tài)系統(tǒng)惡化，2010年，黃旗海湖泊面積也有少量增加，黃旗海濕地LDI值有所下降?？梢钥闯?，雖然黃旗海建立了自然保護(hù)區(qū)，但湖區(qū)不合理的農(nóng)田耕作、草地放牧以及非點(diǎn)源污染仍然是造成LDI值逐年上升的主要原因。雖然設(shè)立了自治區(qū)級自然保護(hù)區(qū)，但只是進(jìn)行了基本的圈地保護(hù)工作，保護(hù)區(qū)內(nèi)放牧以及耕作仍然無法全面禁止，這也解釋了為什么保護(hù)區(qū)內(nèi)有耕地這種土地利用類型。因此，首先應(yīng)采取措施有效禁止湖區(qū)內(nèi)部放牧及耕作現(xiàn)象，給予周邊居民一定的生態(tài)補(bǔ)償是一種可行的選擇。其次，緩沖區(qū)大小與LDI值呈顯著正相關(guān)，這說明黃旗海濕地不具備有效的湖岸帶緩沖區(qū)，可以阻止農(nóng)業(yè)非點(diǎn)源污染與泥沙向湖區(qū)內(nèi)部的運(yùn)移作用。因此，建議開展生態(tài)護(hù)岸建設(shè)，在湖區(qū)周邊設(shè)置有效的湖岸帶緩沖區(qū)。

［1］ 蔣衛(wèi)國，李京，李加洪，等.遼河三角洲濕地生態(tài)系統(tǒng)健康評價(jià)［J］.生態(tài)學(xué)報(bào)，2005，25（3）：408－414.

［2］ 楊永興.國際濕地科學(xué)研究的主要特點(diǎn)，進(jìn)展與展望［J］.地理科學(xué)進(jìn)展，2002，21（2）：111－120.

［3］ Martha A S，Eric D S，Joshua N C，et al.A practical guide for the development of wetland assessment method：the california experience［J］.Journal of the American Water Resources Association，2006，2（1）：157－175.

［4］ 胡小龍，李錦榮，薛博，等.農(nóng)牧交錯(cuò)區(qū)濕地景觀動(dòng)態(tài)變化分析：以內(nèi)蒙古多倫縣為例［J］.水土保持研究，2013，20（1）：116－122.

［5］ Stapanian M A，Waite T A，Krzys G，et al.Rapid assessment indicator of wetland integrity as an unintended predictor of avian diversity［J］.Hydrobiologia，2004，520（1／3）：119－126.

［6］ Brown M T，Vivas M B.Landscape development intensity index［J］.Environmental Monitoring and Assessment，2005，101（1／3）：289－309.

［7］ Lane C R，Brown M T.Diatoms as indicators of isolated herbaceous wetland condition in Florida，USA ［J］.Ecological Indicators，2007，7（3）：521－540.

［8］ Mack J J.Landscape as a predictor of wetland condition：an evaluation of the Landscape Development Index（LDI）with a large reference wetland dataset from Ohio［J］.Environmental Monitoring and Assessment，2006，120（1／3）：221－241.

［9］ 趙雙吉，高文海，李建平，等.黃旗海歷史演變的分析及發(fā)展方向［J］.內(nèi)蒙古農(nóng)業(yè)科技，2002（1）：27－29.

［10］ 李華章，劉清泗，汪家興.內(nèi)蒙古高原黃旗海、岱海全新世湖泊演變研究［J］.湖泊科學(xué)，1992（1）：31－39.

［11］ 李栓科，王濤.全新世內(nèi)蒙南緣黃旗海湖面的波動(dòng)［J］.湖泊科學(xué)，1993（4）：324－334.

［12］ 陳忞忞，劉俊國.黃旗海流域景觀格局動(dòng)態(tài)分析［J］.西北林學(xué)院學(xué)報(bào)，2014（待刊）.

［13］ Reiss K C，Brown M T，Lane C R.Characteristic community structure of Florida′s subtropical wetlands：the Florida wetland condition index for depressional marshes，depressional forested，and flowing water forested wetlands［J］.Wetlands Ecology and Management，2010，18（5）：543－556.

［14］ Odum H T.Environmental Accounting：Emergy and Environmental Decision Making［J］.New York：John Wiley and Sons，1996：20－50.

［15］ 李雙成，傅小鋒，鄭度.中國經(jīng)濟(jì)持續(xù)發(fā)展水平的能值分析［J］.自然資源學(xué)報(bào)，2001，16（4）：297－305.

［16］ 梁春玲，谷勝利.南四湖濕地生態(tài)系統(tǒng)能值分析與區(qū)域發(fā)展［J］.水土保持研究，2012，19（2）：185－189.

［17］ 韓瑛，馮文勇.寧夏移民區(qū)種植業(yè)生態(tài)系統(tǒng)的能值分析［J］.水土保持研究，2010，17（6）：263－267.

［18］ 藍(lán)盛芳，欽佩.生態(tài)系統(tǒng)的能值分析［J］.應(yīng)用生態(tài)學(xué)報(bào)，2001，12（1）：129－131.

［19］ 張水龍，莊季屏.農(nóng)業(yè)非點(diǎn)源污染研究現(xiàn)狀與發(fā)展趨勢［J］.生態(tài)學(xué)雜志，1998，17（6）：52－56.

［20］ 張蓉.農(nóng)用地利用方式與非點(diǎn)源污染控制政策研究［D］.杭州：浙江大學(xué)，2012.

［21］ Chapagain A K，Hoekstra A Y，Savenije H H G，et al.The water footprint of cotton consumption：An assessment of the impact of worldwide consumption of cotton products on the water resources in the cotton producing countries［J］.Ecological Economics，2006，60（1）：186－203.

［22］ 郭婧.典型草原放牧生態(tài)系統(tǒng)能值分析與評價(jià)［D］.呼和浩特：內(nèi)蒙古大學(xué)，2009.

［23］ 王瑗.流域尺度景觀格局時(shí)空演變與生態(tài)系統(tǒng)健康評價(jià)研究［D］.長春：東北師范大學(xué)，2012.

［24］ 內(nèi)蒙古統(tǒng)計(jì)年鑒［M］.呼和浩特：中國統(tǒng)計(jì)出版社，1993—2010.