基于土地利用松花江流域生態(tài)環(huán)境空間分異特征的研究——以哈爾濱段為例

龔文峰

(黑龍江大學(xué) 水利電力學(xué)院,哈爾濱 150086)

土地利用/土地覆被變化(LUCC)是全球環(huán)境變化的重要組成部分和主要原因之一[1],LUCC改變土地覆被狀況并影響許多生態(tài)過程[2],尤其對(duì)區(qū)域氣候、土壤、水量和水質(zhì)的影響是極其深刻的[3]。流域作為完整的自然地理單元[4],LUCC直接制約著流域內(nèi)的水文過程、水化學(xué)過程、生物過程等自然過程的發(fā)生發(fā)展[5]。然而,長(zhǎng)期以來,由于受傳統(tǒng)觀念及人口快速增長(zhǎng)等因素的影響[6],以流域開發(fā)為主的持續(xù)不斷的水利建設(shè)活動(dòng)[7],雖然促進(jìn)農(nóng)林牧業(yè)的發(fā)展,但對(duì)生態(tài)環(huán)境也帶來不利影響,特別是毀林開荒、陡坡耕種現(xiàn)象嚴(yán)重[6],結(jié)果造成土地利用狀況變化頻繁,水土流失逐年加劇、生態(tài)環(huán)境不斷惡化,嚴(yán)重制約著研究區(qū)域社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展。

松花江流域作為黑龍江省的重要糧食產(chǎn)區(qū),是研究資源環(huán)境動(dòng)態(tài)變化比較典型的區(qū)域。因此,在RS、GIS技術(shù)的支持下,以松花江流域哈爾濱段作為研究對(duì)象,在分析其土地利用/土地覆蓋變化的基礎(chǔ)上,從自然環(huán)境因子和社會(huì)因子定量完成對(duì)區(qū)域生態(tài)環(huán)境的綜合評(píng)價(jià),并探索區(qū)域生態(tài)環(huán)境演變的過程和空間分異特征,為生態(tài)環(huán)境恢復(fù)與重建、土地可持續(xù)利用決策提供科學(xué)依據(jù)[8-10],同時(shí)對(duì)流域水土流失治理和生態(tài)安全保護(hù)以及東北黑土區(qū)水土流失的防治規(guī)劃,實(shí)施退耕還林還草和農(nóng)林復(fù)合等措施都具有重要的意義。

1 研究區(qū)概況及研究方法

1.1 研究區(qū)域

松花江流域 (哈爾濱段)位于黑龍江省中南部,地處E125°40′17″～E129°33′2″,N45°12′6″～N49°40′18″,研究區(qū)域東部為中低山丘陵而西部為平原,地勢(shì)總趨勢(shì)是東高西低,年平均氣溫2.5～4.0℃,年平均無霜期 110～150 d,年降水量400～600 mm,年積溫在2 500～3 100℃,土地肥沃、地域遼闊、人口相對(duì)較少,適宜發(fā)展農(nóng)業(yè)生產(chǎn),是重要的糧食生產(chǎn)基地。

1.2 研究方法

1.2.1 圖件配準(zhǔn)及影像解譯

以1∶50 000地形圖為基礎(chǔ),在 ERDAS IMAGINE 9.2遙感圖像處理軟件中,采用二次多項(xiàng)式糾正方法對(duì)1996(6月15日)、2005年(8月27日)Landsat TM/ETM+圖像進(jìn)行幾何校正(平均誤差＜0.5個(gè)像元),并對(duì)圖像進(jìn)行區(qū)域裁剪。以 《中國(guó)土地分類系統(tǒng)》(2001)為標(biāo)準(zhǔn),在ArcGIS9.2環(huán)境下結(jié)合土地利用現(xiàn)狀和野外調(diào)查資料對(duì)影像進(jìn)行解譯并數(shù)字化,建立空間拓?fù)潢P(guān)系,生成土地利用的數(shù)據(jù)庫(kù)。土地利用方式分為耕地(水田和旱田)、林地 (有林地、灌木林和疏林地)、草地、水域(河流、湖泊和水庫(kù)池塘)、建設(shè)用地 (城鄉(xiāng)工礦居民用地和其它建設(shè)用地)和未利用地 (沙地、沼澤等)6個(gè)一級(jí)類型。

1.2.2 生態(tài)空間分異

以研究區(qū)的DEM為基礎(chǔ),運(yùn)用GIS空間分析和地統(tǒng)計(jì)分析功能,按5 km×5 km進(jìn)行小區(qū)劃分,采樣方式采用等間距[11-12],對(duì)松花江流域進(jìn)行分割,共劃分為22個(gè)小區(qū),建立LUCC變化的生態(tài)環(huán)境綜合評(píng)價(jià)單元——小流域,進(jìn)而綜合評(píng)價(jià)并統(tǒng)計(jì)每個(gè)小流域的生態(tài)環(huán)境質(zhì)量。

1.2.3 評(píng)價(jià)指標(biāo)體系的構(gòu)建

根據(jù)研究區(qū)域生態(tài)環(huán)境的特點(diǎn)和所獲環(huán)境信息的容量,按照綜合性、代表性、實(shí)用性[13-14]、主導(dǎo)因素和可操作性等原則[15],選擇反映土地生態(tài)環(huán)境的生物豐度、植被覆蓋、水系分布、土地質(zhì)量和人為活動(dòng)等方面的綜合指標(biāo),構(gòu)建區(qū)域生態(tài)環(huán)境狀況評(píng)價(jià)指標(biāo)體系,該體系涵蓋自然、農(nóng)業(yè)和人為生態(tài)環(huán)境等方面,并確定各指標(biāo)的土地利用類型的權(quán)重。

1)生物豐度指數(shù):反映區(qū)域內(nèi)的生物多樣性狀況,根據(jù)不同土地利用/覆被類型的生物豐度狀況,確定相應(yīng)的權(quán)重,計(jì)算公式為:生物豐度指數(shù)=(0.35×林地面積+0.21×草地面積+0.28×水域濕地面積+0.11×耕地面積+0.04×建設(shè)用地面積+0.01×未利用地面積)/區(qū)域面積。

2)植被覆蓋指數(shù):反映區(qū)域植被覆蓋的程度,計(jì)算公式為:植被覆蓋指數(shù)=(0.38×林地面積+ 0.34×草地面積+0.19×耕地面積+0.07×建設(shè)用地面積+0.02×未利用地面積)/區(qū)域面積。

3)水網(wǎng)密度指數(shù):反映區(qū)域內(nèi)河流、水域的豐富程度,計(jì)算公式為:水網(wǎng)密度指數(shù)=水體面積/區(qū)域面積。

4)土地退化指數(shù):水土流失是影響研究區(qū)域土地質(zhì)量的關(guān)鍵問題,采用土壤侵蝕程度來表征區(qū)域土地退化的程度。土地退化指數(shù)=(0.01×微度侵蝕面積+0.04×輕度侵蝕面積+0.25×中度侵蝕面積+0.7×重度侵蝕面積)/區(qū)域面積。

5)人為影響指數(shù)[16-17]:描述一定區(qū)域內(nèi)景觀受人為活動(dòng)的影響強(qiáng)度。

式中N為土地利用類型數(shù)量;Ai是第i種地類的面積;Pi為第i種地類的人為影響強(qiáng)度參數(shù);A為區(qū)域總面積。

HAI值在0～1之間變動(dòng),HAI的數(shù)值越大,表明人為活動(dòng)占優(yōu)勢(shì)的地類構(gòu)成越大,人類活動(dòng)影響強(qiáng)度越大;反之則小,見表1。

表1 各種人為影響強(qiáng)度參數(shù)Table 1 Human activity density index

1.2.4 評(píng)價(jià)指標(biāo)處理方法

為保持各指標(biāo)數(shù)據(jù)間的可比性[18],利用歸一化方法對(duì)各因素進(jìn)行處理[19-20],得出最終所需的各個(gè)指標(biāo)值,即:式中Xi為某一參評(píng)因子第i級(jí)值;Xmax為某一參評(píng)因子最高閾值;Xmin為某一參評(píng)因子最低閾值。

1.2.5 評(píng)價(jià)指標(biāo)權(quán)重的確定

采用Delphi和AHP兩種方法加權(quán)求平均值[17],確定各評(píng)價(jià)指標(biāo)的權(quán)重,結(jié)果見表2。

表2 研究區(qū)域生態(tài)環(huán)境評(píng)價(jià)指標(biāo)權(quán)重T able 2 Index weights of eco-environmental statusevaluation

1.2.6 評(píng)價(jià)方法

根據(jù)指標(biāo)體系、指標(biāo)權(quán)重以及實(shí)際研究的需要,松花江流域生態(tài)環(huán)境質(zhì)量綜合評(píng)價(jià)采用生態(tài)環(huán)境質(zhì)量指數(shù)法(EI),其計(jì)算方法為:EI=0.3×生物豐度指數(shù)+0.25×植被覆蓋指數(shù)+0.25×水網(wǎng)密度指數(shù)+0.2×(1-土地退化指數(shù))+0.1×人為影響指數(shù)。

2 結(jié)果與分析

2.1 土地利用變化分析

對(duì)研究區(qū)域兩期土地利用數(shù)據(jù)進(jìn)行分析得知(表3),土地利用類型以林地、耕地和草地為主,三者面積之和在各時(shí)期均占總面積的87%以上,這種特征決定了該區(qū)域以農(nóng)林牧業(yè)生產(chǎn)為主的特點(diǎn)。林地和水域面積增加,年變化率為1.314%和3.068%。水域增加主要?dú)w因于2005年受季節(jié)降水影響,林地的增加主要?dú)w功于 “天保工程”、松花江流域 “生態(tài)環(huán)境綜合治理”等多項(xiàng)生態(tài)修復(fù)工程的推行。耕地、草地、未利用地和建設(shè)用地的年變化率分別為-1.356%、-5.308%、-1.568%和-0.085%,耕地面積減少的最大達(dá)到1 230.2 hm2,草地資源下降最快,未利用地次之,建設(shè)用地面積減少最少。耕地減少很大程度受到退耕還林政策的影響,此外,部分耕地被水域淹沒也是其減少的重要原因。

表3 研究區(qū)域各土地利用類型面積及其變化Table 3 Change in area of different land use types and their dynamic degrees

利用馬爾科夫模型研究和分析土地利用空間變化過程(表4):林地、水域和建設(shè)用地具有較高的保留率;草地、未利用地和耕地具有較高的轉(zhuǎn)換率。林地增加主要源于草地、未利用地和耕地的轉(zhuǎn)化,其中草地轉(zhuǎn)化為林地的比例最大,為56.091%,其次是耕地,所占比例為7.999%,反映人類活動(dòng)正深刻影響著土地利用類型,其驅(qū)動(dòng)機(jī)制歸功于農(nóng)業(yè)結(jié)構(gòu)的調(diào)整和和區(qū)域生態(tài)環(huán)境的建設(shè);耕地的補(bǔ)充主要源于草地和未利用地 (沼澤地)轉(zhuǎn)化,表明在經(jīng)濟(jì)利益的驅(qū)使下,人們加大對(duì)該地類的開發(fā)利用強(qiáng)度,致使土地利用結(jié)構(gòu)趨向不合理,人為因素致使土地利用結(jié)構(gòu)和生態(tài)過程中造成生態(tài)環(huán)境退化;草地主要由未利用地(沼澤)轉(zhuǎn)化而來;水域變化不大,其補(bǔ)充主要源于耕地;建設(shè)用地主要來源于林地和草地的轉(zhuǎn)化;未利用地主要轉(zhuǎn)化為林地、耕地和草地。

表4 土地利用狀態(tài)轉(zhuǎn)移概率矩陣Table 4 Land use change matrix of source region of study area /%

2.2 生態(tài)環(huán)境的空間分異

研究結(jié)果表明(表5),各評(píng)價(jià)單元高低懸殊,差異較大,但是相對(duì)差異在逐步縮小。1996年評(píng)價(jià)單元的生態(tài)環(huán)境評(píng)價(jià)的最大值為7.5,最小值為1.846 8,最大值是最小值的4.06倍,2005年減為3.97倍,且最大和最小值分別為7.5和1.889 4。運(yùn)用統(tǒng)計(jì)分析軟件Static 6.0 for windows對(duì)評(píng)價(jià)單元進(jìn)行聚類分析,參照 《生態(tài)環(huán)境狀況評(píng)價(jià)技術(shù)規(guī)范(試行)》并結(jié)合實(shí)際情況,將該區(qū)域生態(tài)環(huán)境分為5級(jí),即優(yōu)、良、一般、較差和差。由圖1和表6可見:研究區(qū)域生態(tài)環(huán)境為優(yōu)的主要分布在通河縣、方正縣的南部和北部地區(qū)以及方正縣的西南部,該區(qū)域多為山地,林地占據(jù)較大的比例,植被覆蓋度高,人為的干擾度較小,土地退化程度較小,生態(tài)多樣性豐富,生態(tài)系統(tǒng)穩(wěn)定,其中所占比例分別為0.33%和6.38%,2005年的面積要遠(yuǎn)大于1996年的,其中增加的區(qū)域主要集中在方正縣的西南地區(qū);生態(tài)環(huán)境質(zhì)量為良的區(qū)域面積較大,該區(qū)域主要分布在木蘭縣、方正縣和通河縣,所占的比例分別為30.25%和27.67%,其中1996年比2005年多出745.62 km2,減少的區(qū)域分布在木蘭縣的西北部,該區(qū)域處于低山丘陵地帶,多為林地、未利用地、水域和部分耕地,植被覆蓋度較高,受到一定程度的人為干擾,生物多樣性也較豐富,林地和耕地之間的轉(zhuǎn)化是土地利用轉(zhuǎn)化的主要軌跡,是今后要繼續(xù)推行 “退耕還林”政策的重點(diǎn)區(qū)域,加強(qiáng)對(duì)土地合理利用格局的優(yōu)化配置;生態(tài)環(huán)境質(zhì)量一般的主要分布在巴彥縣、賓縣和呼蘭區(qū)的西部區(qū)域,該區(qū)域處于丘陵和平原區(qū)域的過渡地帶,且占據(jù)一定的比例,分別為 25.85%和31.04%,2005年的比例要稍高于1996年,多為耕地和部分林地、草地,植被蓋度一般,有一定程度的人為干擾,生態(tài)系統(tǒng)一般,土地利用主要集中在林地、草地和耕地之間的轉(zhuǎn)換,要加大對(duì)草地和棄耕地的治理和改造,促使生態(tài)環(huán)境的好轉(zhuǎn);土地生態(tài)環(huán)境較差的區(qū)域主要集中分布在雙城市、阿城區(qū)、呼蘭區(qū)和哈爾濱市等,該區(qū)域占據(jù)較大的比例,分別為34.28%和30.15%,多為耕地、水域和未利用地(濕地),有一定的植被蓋度,處于平原區(qū)域,農(nóng)村居民點(diǎn)分布較多,且人口增長(zhǎng)、城鎮(zhèn)化進(jìn)程較快,對(duì)耕地需求較大,濕地、草地和耕地之間的轉(zhuǎn)換較為明顯,土地壓力過重,土地利用趨向不合理,其中人口急劇增長(zhǎng)和城市化的迅速發(fā)展以及農(nóng)業(yè)結(jié)構(gòu)的調(diào)整是該區(qū)域生態(tài)環(huán)境一般的直接原因;生態(tài)環(huán)境差的區(qū)域主要集中在呼蘭區(qū)的東北部和雙城市的西南部,2005年比1996年減少的區(qū)域?yàn)殡p城市的西南部,雖然占據(jù)較小的比例,分別為9.3%和4.76%,多為耕地,植被蓋度較低,城鎮(zhèn)居民點(diǎn)分布較多,耕地和建設(shè)用地的強(qiáng)烈變化尤為明顯,部分區(qū)域依舊存在著濫墾草地、毀林開荒等現(xiàn)象,人為干擾程度強(qiáng),生態(tài)環(huán)境較為惡劣,更加速人類活動(dòng)與自然環(huán)境之間的矛盾。

表5 區(qū)域生態(tài)環(huán)境質(zhì)量分級(jí)Table 5 Eco-environmental status classification

表6 區(qū)域生態(tài)環(huán)境質(zhì)量分布及其變化Table 6 Change in area of eco-environmental status and their dynamic degrees

圖1 區(qū)域生態(tài)環(huán)境質(zhì)量分布圖(1996,2005年)Fig.1 Distribution map in area of eco-environmental status(1996,2005)

2.3 生態(tài)環(huán)境動(dòng)態(tài)分析

在ArcGIS的支持下,對(duì)兩期評(píng)價(jià)圖進(jìn)行疊加分析,在此基礎(chǔ)上檢測(cè)EI的變化情況,并根據(jù)其變化幅度進(jìn)行區(qū)域生態(tài)環(huán)境狀況的動(dòng)態(tài)分級(jí),分為3個(gè)等級(jí),即:略微變差、基本不變、略微變好(表7),同時(shí)生成土地生態(tài)環(huán)境狀況變化分級(jí)圖[17](圖2)。

表7 區(qū)域土地生態(tài)環(huán)境狀況變化分級(jí)Table 7 Classification of land eco-environmental status change degree

表8 土地生態(tài)環(huán)境分級(jí)百分比Table 8 Percentage of land eco-environmental grades

圖2 土地生態(tài)環(huán)境狀況變化分級(jí)Fig.2 Land ecological and environmental status change classification image

結(jié)合表8和圖2可以看出研究區(qū)域1996～2005年土地生態(tài)環(huán)境總體呈變好的趨勢(shì)。其中基本變好和略微變好的區(qū)域面積占總面積的21.39%,多分布在方正縣、雙城市的東南部、巴彥縣的東北部以及賓縣和木蘭縣的部分區(qū)域,該地區(qū)為山地和低山丘陵,林地和耕地占據(jù)一定比例,隨著 “天保工程”、松花江流域 “生態(tài)環(huán)境綜合治理”等多項(xiàng)生態(tài)修復(fù)工程的推行,耕地——林地轉(zhuǎn)變成為區(qū)域土地利用的主要軌跡,為生態(tài)環(huán)境的改善提供必要條件。因此,要繼續(xù)實(shí)施積極的農(nóng)業(yè)調(diào)整政策,加大林業(yè)監(jiān)督管理,制止不合理的開發(fā)利用活動(dòng),正確協(xié)調(diào)人地矛盾,實(shí)現(xiàn)土地資源的合理利用,進(jìn)而形成相對(duì)穩(wěn)定的生態(tài)結(jié)構(gòu)和生態(tài)功能;基本不變的面積占78.17%,處于絕對(duì)的優(yōu)勢(shì)地位,該地區(qū)為平原區(qū)和部分山地區(qū)域,是人們生活和從事農(nóng)業(yè)生產(chǎn)活動(dòng)的主要地段,耕地和林地占較大的比例,因此,要加強(qiáng)土地利用格局的合理優(yōu)化配置,逐步完善并形成良好的黑土地生態(tài)系統(tǒng);略微變差的地區(qū)為低山丘陵和平原區(qū)的交界處,生態(tài)環(huán)境比較脆弱,人類的干擾活動(dòng)比較頻繁,容易發(fā)生土壤侵蝕,是今后水土流失治理的重點(diǎn)。因此,在發(fā)展經(jīng)濟(jì)的同時(shí),要繼續(xù)推行 “退耕還林”政策,因地制宜營(yíng)造水土保持林、用材林、經(jīng)濟(jì)林等[21],實(shí)現(xiàn)土地合理利用,促進(jìn)區(qū)域的可持續(xù)發(fā)展。

3 結(jié) 論

1)該區(qū)域以農(nóng)林牧業(yè)生產(chǎn)為主,林地、耕地和草地三者面積之和在各時(shí)期均占總面積的87%;林地和水域面積增加,其它地類的面積減少,林地增加的最多,耕地減少的最多;土地利用轉(zhuǎn)化過程以耕地-林地、草地-林地和未利用地-林地為主,林地、水域和建設(shè)用地具有較高的保留率,草地、未利用地和耕地具有較高的轉(zhuǎn)換率。

2)通過對(duì)小流域統(tǒng)計(jì)結(jié)果進(jìn)行空間化可知,各評(píng)價(jià)單元高低懸殊,差異較大,1996年生態(tài)環(huán)境評(píng)價(jià)的最大和最小值分別為7.5和1.846 8,而2005年的為7.5和1.889 4;土地生態(tài)環(huán)境好多為山地區(qū)域,林地占據(jù)較大比例,生態(tài)系統(tǒng)穩(wěn)定;生態(tài)環(huán)境較好多為丘陵地帶,生態(tài)環(huán)境較差和差的區(qū)域多為平原地帶,人為干擾程度強(qiáng),生態(tài)環(huán)境較為惡劣。

3)從空間分布來看,1996～2005年土地生態(tài)環(huán)境總體呈變好的趨勢(shì)。其中變好和略微變好的區(qū)域面積占總面積的21.39%,基本不變的面積占78.17%,略微變差的地區(qū)占據(jù)較小的比例,為0.44%。

[1]陳百明,劉新衛(wèi),楊 紅.LUCC研究的最新進(jìn)展評(píng)述[J].地理科學(xué)進(jìn)展,2003,22(1):22-29.

[2]Turner M G.Landscape Ecology:the Effect of Pattern on Process[J].Annual Review of Ecology and Systematics,1989,20:171-197.

[3]郭旭東,陳利頂,傅伯杰.土地利用/土地覆被變化對(duì)區(qū)域生態(tài)環(huán)境的影響 [J].環(huán)境科學(xué)進(jìn)展,1999,7 (6):66-75.

[4]竇 燕,陳 曦,包 安.近40年和田河流域土地利用動(dòng)態(tài)變化及其生態(tài)環(huán)境效應(yīng) [J].干旱區(qū)地理, 2008,31(3):449-455.

[5]黃 方,劉湘南,劉 權(quán),等.遼河中下游流域土地利用變化及其生態(tài)環(huán)境效應(yīng) [J].水土保持通報(bào), 2004,24(6):18-21.

[6]王思遠(yuǎn).王光謙.陳志樣.黃河流域生態(tài)環(huán)境綜合評(píng)價(jià)及其演變[J].山地學(xué)報(bào),2004,22(2):133-139.

[7]吳貽名,李元紅,朱 強(qiáng),等.干旱區(qū)流域生態(tài)環(huán)境質(zhì)量現(xiàn)狀評(píng)價(jià)研究 [J].中國(guó)農(nóng)村水利水電,2001 (1):25-30.

[8]陳國(guó)強(qiáng),陳 鵬.海灣型城市土地利用動(dòng)態(tài)演變及其生態(tài)環(huán)境效應(yīng)研究 [J].水保保持研究,2008,15 (3):26-30.

[9]王瑞燕,趙庚星,周 偉,等.土地利用對(duì)生態(tài)環(huán)境脆弱性的影響評(píng)價(jià) [J].農(nóng)業(yè)工程學(xué)報(bào),2008,24 (12):215-129.

[10]梁吉順,張艷萍,馮井泉.花江流域植被覆蓋動(dòng)態(tài)監(jiān)測(cè)遙感分析—哈爾濱市區(qū)段為例[J].黑龍江水專學(xué)報(bào),2009,36(3):74-76.

[11]岳書平,張樹文,閆業(yè)超,等.公主嶺市土地利用變化的生態(tài)效應(yīng)研究 [J].資源科學(xué),2006,28(6): 161-166.

[12]岳書平,張樹文,閆業(yè)超,等.東北典型農(nóng)區(qū)土地利用變化的生態(tài)效應(yīng)研究—以公主嶺市為例 [J].干旱區(qū)資源與環(huán)境,2007,21(7):64-68.

[13]周華榮.新疆生態(tài)環(huán)境質(zhì)量評(píng)價(jià)指標(biāo)體系研究 [J].中國(guó)環(huán)境科學(xué),2000,20(2):150-153.

[14]李新琪,劉建軍,朱海涌,等.新疆環(huán)境質(zhì)量綜合評(píng)價(jià)研究[J].干旱環(huán)境監(jiān)測(cè),2003,17(2):82-85.

[15]馬乃喜.區(qū)域生態(tài)環(huán)境評(píng)價(jià)中的幾個(gè)理論問題 [J].西北大學(xué)學(xué)報(bào),1998,28(4):330-334.

[16]陳 浮,葛小平,陳 剛,等.城市邊緣區(qū)景觀變化與人為影響的空間分異研究 [J].地理學(xué)報(bào),2001, 21(3):210-216.

[l7]王靜愛,何春陽(yáng),董艷春,等.北京城鄉(xiāng)過渡區(qū)土地利用變換驅(qū)動(dòng)力分析[J].地球科學(xué)進(jìn)展,2002,17 (2):201-208.

[18]程晉南,趙庚星,李 紅,等.基于RS和GIS的土地生態(tài)環(huán)境狀況評(píng)價(jià)及其動(dòng)態(tài)變化[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(11):83-88.

[19]鄭新奇,王愛萍.基于RS與GIS的區(qū)域生態(tài)環(huán)境質(zhì)量綜合評(píng)價(jià)研究——以山東省為例[J].環(huán)境科學(xué)學(xué)報(bào),2000,20(4):489-493.

[20]仲嘉亮,謝 勇,朱海涌.塔里木河流域的生態(tài)環(huán)境質(zhì)量綜合評(píng)價(jià)研究 [J].干旱環(huán)境監(jiān)測(cè),2004,18 (4):203-207.

[21]劉寶元,閻百興,沈 波,等.東北黑土區(qū)農(nóng)地水土流失現(xiàn)狀與綜合治理對(duì)策 [J].中國(guó)水土保持科學(xué), 2008,6(1):1-8.