農(nóng)墾區(qū)土壤侵蝕及生態(tài)效益監(jiān)測項目技術(shù)方案

農(nóng)墾區(qū)土壤侵蝕及生態(tài)效益監(jiān)測項目技術(shù)方案總體實施方案查閱國內(nèi)外土壤侵蝕模型相關(guān)文獻，研究土壤侵蝕模型，并根據(jù)模型運行要求，收集新疆兵團五期（80年代至今）基礎數(shù)據(jù)，遙感數(shù)據(jù)，成果數(shù)據(jù)和相關(guān)統(tǒng)計資料等，結(jié)合兩次野外調(diào)查監(jiān)測數(shù)據(jù)，以一期（2010年）數(shù)據(jù)為例，提取模型所需相關(guān)因子進行模型試運行，調(diào)整相關(guān)參數(shù)，對結(jié)果進行驗證分析；基于成功運行的土壤侵蝕模型，開發(fā)模型軟件，應用于其他四期的土壤侵蝕監(jiān)測中；以五期多源遙感影像為數(shù)據(jù)源，提取農(nóng)墾區(qū)界線及土地利用類型，分析其變化規(guī)律，并與提取的相應時期的土壤侵蝕結(jié)果對比，分析農(nóng)墾區(qū)對土壤侵蝕的影響狀況。綜合各種相關(guān)資料，選取研究示范區(qū)水土保持生態(tài)效益評價指標，建立指標體系，并提取五期水土保持生態(tài)效益指標數(shù)據(jù)結(jié)果，結(jié)合相應時期的水土保持措施情況，分析水土保持生態(tài)效益的變化規(guī)律，從而完成項目工作報告。技術(shù)流程土壤風蝕定量監(jiān)測土壤風蝕模型選取與修正通過對國內(nèi)外文獻的大量查閱，發(fā)現(xiàn)國外模型中沒有一種適合研究區(qū)特點的模型，繼而轉(zhuǎn)向國內(nèi)研究，近年來，隨著國內(nèi)學者在野外風蝕觀測、風洞實驗模擬和數(shù)學模型理論研究方面的不斷積累和努力，國內(nèi)學者站在我國自身角度也提出了不少模型，最終基于中科院寒區(qū)旱區(qū)環(huán)境與工程研究所提出的風洞實驗模型基礎上，選擇了該模型的大田推廣模型。該模型是通過各類地表在不同植被和風力條件下計算土壤風蝕模數(shù)的估算方法。具體說，是將非風蝕地表類型（如水體、道路、城鎮(zhèn)用地等）以外的易風蝕地表劃分為耕地、林（草）地（包含退耕地）、沙地三大類，分別針對不同地表類型估算土壤風蝕量。原模型共分為三個子模型，耕地模型，林草地模型和沙地模型，分別為：（1）耕地模型（2）林草地模型（3）沙地模型結(jié)合研究區(qū)概況和風蝕影響因素的綜合分析，對三個子模型分別進行了修正。由于研究區(qū)極端干旱農(nóng)墾區(qū)的氣候特點和耕作制度（耕地采用引水灌溉方法），加入土壤含水量因子（根據(jù)土壤含水率等級，對土壤含水量因子賦值），對風蝕起到阻礙作用；不同的土地利用類型中，土壤質(zhì)地（粒徑大于0.84mm，則對風蝕起到阻礙作用，根據(jù)粒徑含量的百分比進行S因子的賦值（0-1之間））不同，對風蝕的阻礙作用也不同，故加入；干旱農(nóng)墾區(qū)為保護性農(nóng)田，在農(nóng)田田埂周圍種植防護林，對風蝕同樣起到阻礙作用，在公路兩邊，設計草方格等沙障措施，對風蝕起到阻礙作用，故加入P因子（根據(jù)防護林和沙障對風蝕的影響強度和影響范圍，給P因子賦值）。修正后的土壤風蝕定量模型仍然由三個模型公式組成，分別為耕地模型、林草地模型和沙地模型。具體模型公式如下：（1）耕地模型修正后為：式中：為耕地模型區(qū)土壤風蝕模數(shù)[t/hm2·a]；是大于臨界侵蝕風速的第j級風速（指氣象站統(tǒng)計的整點風速）；為風蝕活動發(fā)生月份內(nèi)風速為Ui的累積時間（min）；第一個風速（大于臨界侵蝕風速）等級為5.0-6.0m/s，取平均值5.5m/s，第二個風速等級為6.0-7.0m/s，取平均值6.5m/s，……。是最大風速（指氣象站統(tǒng)計的整點風速）；a1、b1、c1為與土壤類型有關(guān)的常數(shù)項，暫取-9.208、0.018、1.955；為尺度修訂系數(shù)，約為0.0018；A為風速修訂系數(shù)，計算耕地模型區(qū)風蝕模數(shù)時需要利用當?shù)貧庀笳居涗浀娘L速資料，風洞軸線風速與氣象站10米觀測高度測定的風速有一定的差距，因此應將氣象站記錄的風速換算為風洞軸線風速，在假定風洞實驗動力相似得到保證的前提下，A是與下墊面有關(guān)的風速修訂系數(shù)。（2）林草地模型修正后為：式中：為林草地模型區(qū)林草地土壤風蝕模數(shù)[t/hm2·a]；是大于臨界侵蝕風速的第j級風速（指氣象站統(tǒng)計的整點風速）；為風蝕活動發(fā)生月份內(nèi)風速為的累積時間（min）；為為尺度修訂系數(shù)，約為0.0018；A為風速修訂系數(shù)；a2、b2、c2為與土壤類型有關(guān)的常數(shù)項，暫取2.4869、-0.0014、-54.9472。（3）沙地模型修正后為：式中：為沙地模型區(qū)土壤風蝕模數(shù)[t/hm2·a]；是大于臨界侵蝕風速的第j級風速（指氣象站統(tǒng)計的整點風速）；沙地=5.5m/s；為風蝕活動發(fā)生月份內(nèi)風速為的累積時間（min）。為尺度修訂系數(shù)，約為0.0018，A為風速修訂系數(shù)；VC為植被覆蓋度（%）；a3、b3、c3為與土壤類型有關(guān)的常數(shù)項，暫取6.1689、-0.0743、-27.9613。關(guān)于修正模型的說明：第一，所采用的風蝕模型沒有考慮地形的作用。在風蝕過程中，地形的作用遠不如在水蝕過程中那樣顯著。第二，所采用的風蝕模型沒有考慮林地和草地的區(qū)別，即將林、草地風蝕模數(shù)按相同的風蝕模數(shù)計算。第三，模型將除沙地以外的未利用地作為相當植被蓋度下的林（草）地處理。修正后的模型更能反映干旱農(nóng)墾區(qū)風蝕的特點，能更準確真實地反映其風蝕狀況，提高模型在極端干旱農(nóng)墾區(qū)的計算精度。根據(jù)模型運行特點，將模型區(qū)總共劃分為四大類，分別為耕地模型區(qū)，林草地模型區(qū)（包含未利用地，不包含沙地），沙地模型區(qū)和非風蝕地表模型區(qū)。但通常在河流水面、灘涂、水庫水面、建筑用地、坑塘水面、農(nóng)田利用地等非風蝕地表類型，往往是風沙堆積和沙塵沉降的區(qū)域，幾乎不會發(fā)生風蝕，其土壤風蝕量可視為零，故可采用GIS的空間分析功能，將以上地表類型摳除，使其不參與風蝕量計算中，一方面提高了模型運算的效率，另一方面也提高了模型計算精度，因此本報告中將只討論耕地模型區(qū)、林草地模型區(qū)和沙地模型區(qū)三種類型。野外調(diào)查野外調(diào)查共分三次第一次（野外調(diào)查為主）：2010年10月28日至2010年11月17日，按照隨機抽樣的方式，產(chǎn)生13個抽樣單元，對新疆兵團農(nóng)一師墾區(qū)相關(guān)指標進行調(diào)查，包括土壤特性、土壤含水量、樣區(qū)內(nèi)的植被蓋度及高度情況、作物高度及作物殘茬高度、連續(xù)觀測記錄每月測釬高度變化情況。第二次（野外調(diào)查為主）:2011年4月9日至2011年6月17日，按照隨機抽樣的方式，產(chǎn)生48個抽樣單元，對新疆兵團農(nóng)一師墾區(qū)相關(guān)指標進行調(diào)查，包括周圍地表狀況、土地利用、植被蓋度、土壤含水量、土壤質(zhì)地、地表粗糙度、風蝕措施等。第三次（資料收集為主）：2011年7月-2012年6月包括各團年鑒資料、各團氣象數(shù)據(jù)、各團團志資料、土地利用變化。數(shù)據(jù)預處理1、數(shù)據(jù)源獲取1）遙感影像數(shù)據(jù)（1）1976年左右的LandsatMSS數(shù)據(jù)（8、9、10月）（2）1990年左右數(shù)據(jù)（LandsatTM）（7、8、9、10月）（3）2000年左右數(shù)據(jù)（LandsatETM）（6、8、9、10月）（4）2006年左右數(shù)據(jù)（LandsatETM）（5、6、9、10月）（5）2006年SPOT數(shù)據(jù)（7，8月份）（6）2010年HJ數(shù)據(jù)（3-6月份和9-11月份）2）矢量數(shù)據(jù)（1）2006年土地利用數(shù)據(jù)（人工數(shù)字化結(jié)果）（2）2006年植被蓋度數(shù)據(jù)（人工數(shù)字化結(jié)果）（3）2006年土壤侵蝕數(shù)據(jù)（人工數(shù)字化結(jié)果）（4）農(nóng)一師團界數(shù)據(jù)3）野外調(diào)查監(jiān)測數(shù)據(jù)（1）1976、1990、2000、2006和2010年一團和五團氣象數(shù)據(jù)（2）2010和2011年風蝕模型因子野外監(jiān)測數(shù)據(jù)（地表粗糙度、土壤質(zhì)地、植被蓋度、土壤含水量、風蝕措施）4）新疆兵團農(nóng)一師年鑒與統(tǒng)計數(shù)據(jù)2、數(shù)據(jù)預處理數(shù)據(jù)收集包括所有應用于項目的基礎數(shù)據(jù)，遙感數(shù)據(jù)，監(jiān)測數(shù)據(jù)和相關(guān)統(tǒng)計資料等。數(shù)據(jù)預處理包括文檔資料整理，信息統(tǒng)計，空間數(shù)據(jù)信息的屬性信息檢查修正，投影信息統(tǒng)一，遙感數(shù)據(jù)的預處理等。五期遙感數(shù)據(jù)（1976、1990、2000、2006和2010年）經(jīng)過幾何校正或正射校正，波段疊加，輻射定標，反射率計算，植被指數(shù)計算，拼接，裁剪進行預處理；氣象數(shù)據(jù)和農(nóng)一師年鑒通過手工錄入計算機，得到電子版數(shù)據(jù)；2010和2011年風蝕模型因子野外監(jiān)測數(shù)據(jù)（地表粗糙度、土壤質(zhì)地、植被蓋度、土壤含水量、風蝕措施）經(jīng)過處理后錄入計算機得到電子版，再在ArcGIS中建立相應字段并為整個研究區(qū)所有土地利用類型賦值。土壤風蝕定量模型因子提取需要提取的土壤風蝕模型因子有風力因子，土壤含水量因子，土壤質(zhì)地因子，植被蓋度因子，地表粗糙度因子，風蝕措施因子等。土地利用類型圖人工目視解譯五期遙感影像圖獲得。風力因子計算為了對氣象站逐日4次（2時、8時、14時、20時）觀測的風速風向數(shù)據(jù)進行空間內(nèi)插和外推計算。首先采用線性插值的方法，得到24小時整點風速估計值，然后按1m/s間隔統(tǒng)計全年≥5m/s的各等級風速的累積時間。根據(jù)公式（Ⅰ）計算相鄰兩個整點時刻的風速：（Ⅰ）式中：t1、t2—逐日4次風速風向統(tǒng)計表中相鄰兩個風速記錄的時刻（2時、8時、14時、20時）；t—t1～t2之間的整點時間；Ui—與ti相對應的整點風速（m/s）；Ut1—與t1相對應的風速（m/s）；Ut2—與t2相對應的風速（m/s）；按1m/s間隔統(tǒng)計全年≥5m/s（5.0～5.9m/s、6.0～6.9m/s、7.0～7.9m/s、……最高風速統(tǒng)計到該氣象站點記錄到的最大風速）的各等級風速的累積時間。如統(tǒng)計全年5.0～5.9m/s風速發(fā)生的頻數(shù)，每次代表1小時，按下式換算為該等級風速的累積時間。累積時間=頻數(shù)×1hr/次×60min/hr為便于計算，介于5.0～5.9m/s的風速等級取平均值U＝5.5m/s，6.0～6.9m/s的風速等級取平均值U＝6.5m/s，依此類推?？紤]到風力因子在空間上的變化是漸進的，即風速和風向在空間上是逐漸變化的，相鄰氣象站之間的風速和風向變化呈逐漸過渡關(guān)系。基于以上原理，運用距離倒數(shù)算法自主開發(fā)程序進行空間插值，并將每個氣象站點作為空間上的一個點，生成1976、1990、2000、2006和2010年3～6月和9～11月的每月不同風力等級風速累計時間的30m分辨率柵格圖。土壤含水量因子計算對于2010年土壤含水量因子的計算，可根據(jù)野外實測根據(jù)公式（Ⅱ）計算土壤含水量(Ⅱ)式中：m0—烘干空鋁盒質(zhì)量（g）；m1—烘干前鋁盒及土樣質(zhì)量（g）；m2—烘干后鋁盒及土樣質(zhì)量（g）；對于1976年、1990年、2000年和2006年的土壤含水量推算，采用如下估算方法：由于塔里木盆地平原農(nóng)業(yè)區(qū)為干旱區(qū)，原始土壤多為荒漠土，普遍含鹽，這里年降水量通常小于100mm，年蒸發(fā)量2200-2900mm，農(nóng)業(yè)全部依靠灌概，是典型的干早灌溉農(nóng)業(yè)區(qū)。土壤含水量影響因素較多，主要有氣象（特別是日照數(shù)和降雨量）、作物、土壤質(zhì)地、灌溉田間管理，但在本研究中，只考慮灌溉、蒸散發(fā)和降雨量三個主要因素。依據(jù)是否灌溉，首先將土地利用類型圖分成兩大類，一種為非灌溉土壤，一種為灌溉土壤。因為研究區(qū)域范圍較大，不可能對每年每月每天每個點的土壤含水量進行實地測定，因此考慮采用氣候資料，通過降雨、蒸發(fā)等數(shù)據(jù)建立經(jīng)驗模型，關(guān)鍵在于蒸散發(fā)的計算，它是水量平衡中的重要參數(shù)，對模型參數(shù)優(yōu)選結(jié)果和模擬精度的影響很大。對于非灌溉土壤類型（主要有有林地、水土保持林、天然草地、荒草地、鹽堿地、裸土地、裸巖石礫地、固定沙地、半固定沙地和流動沙地十類），只考慮降水和蒸發(fā)對其土壤含水量的影響。對于灌溉土壤類型（主要有灌溉水田、水澆地、果園和人工草地四類），只考慮降水、蒸發(fā)和灌溉三個主要因素對土壤含水量的影響。少量降雨對土壤含水量影響不大，但降雨量較大時，在短時間（一般持續(xù)1-2天）對土壤水分影響比較明顯。而灌溉通常情況下則可以持續(xù)7天左右。據(jù)以上分析，土壤含水量因子提取流程圖如下：圖2-1土壤含水量因子提取流程圖（1）對于非灌溉性土壤，首先根據(jù)平均氣溫，找出氣溫與蒸發(fā)量之間的關(guān)系，通過擬合回歸，關(guān)系式如下：ET=49.158*EXP(0.0835*T)（式中：ET為蒸發(fā)量，T為平均氣溫，單位為℃，ET單位為mm）。由此得出蒸發(fā)量，然后根據(jù)2010年的降水量、蒸發(fā)量和實測土壤含水量值，通過線性擬合回歸，得到關(guān)系式如下：y=0.0183x1-0.0009x2+1.2771，（式中：y表示土壤含水量，x1表示降水量，x2表示蒸發(fā)量，所有值的單位均為mm）。（2）對于灌溉性土壤，以2010年降雨量、蒸發(fā)量和灌溉量作為自變量，土壤含水量作為因變量，擬合回歸公式如下：y=0.1216x1-0.0248x2+0.0100x3+6.7640，（式中：y表示土壤含水量，x1表示降雨量，x2表示蒸發(fā)量，x3表示灌溉量）。土壤質(zhì)地因子計算土壤中粒徑小于0.84mm的顆粒最容易發(fā)生風蝕（Chepil，1955）。在衡量不同土地利用類型的抗風蝕能力時，值越大表明該土地利用類型的土壤抗風蝕能力越強。通過對野外調(diào)查監(jiān)測過程中選取的幾種具有代表性的土地利用類型的土樣，通過環(huán)刀法沿S型取五處表層（0-5cm深度內(nèi)）土壤，進行取樣分析，用20目手動土壤篩篩選出小于0.84mm土壤顆粒，然后用天平稱重，從而得到其土壤粒徑組成情況以及不同土地利用類型的不可蝕性顆粒（>0.84mm）所占質(zhì)量百分比。式中：m1—不可蝕性顆粒的質(zhì)量（g）；m—土樣質(zhì)量（g）；地表粗糙度因子計算地表粗糙度Z0是耕地模型運行必不可少的因子。隨著作物的生長，高度和蓋度會發(fā)生變化，其地表粗糙度值也會相應發(fā)生變化，在水澆地（棉花）、小麥和新墾地三種土地利用類型下于2010年3-6月和9-11月進行連續(xù)監(jiān)測。對于2010年地表粗糙度的計算，采用如下方法：先對風速求平均，再根據(jù)公式（Ⅲ）得到下墊面地表粗糙度。（Ⅲ）其中，Z1,Z2為兩個觀測高度，V1,V2為兩個觀測高度上的風速。對于1976年、1990年、2000年和2006年的地表粗糙度，采用作物高度估算的方法。首先通過野外調(diào)查實測，結(jié)合研究區(qū)農(nóng)作物種植結(jié)構(gòu)狀況，記錄主要作物（棉花、水稻、玉米、小麥）每月生長高度。然后對每月作物的高度求平均，估測每月作物平均高度。查閱《空氣動力學粗糙度》文獻，得知植被覆蓋地表的空氣動力學粗糙度公式，其中玉米地表粗糙度Z0=0.025*H^(1.1)，其他作物地表粗糙度Z0=0.04*H^(1.417)。參考第一次全國水利普查不同下墊面粗糙度值，耙平無壟農(nóng)田地表粗糙度為0.06，留茬無壟農(nóng)田地表粗糙度為0.18，通過對研究區(qū)調(diào)查和分析，水稻（10月至次年3月）、棉花（11月至次年3月）、玉米（10月至次年4月）、小麥（9月至10月）收割后均留茬，水稻和棉花在4月上旬至中旬留茬，4月中旬至下旬翻耕，故地表粗糙度取兩者平均，確定為0.12。植被蓋度因子計算包括五期（1976、1990、2000、2006和2010年）遙感數(shù)據(jù)的幾何校正或正射校正，波段疊加，輻射定標，反射率計算，植被指數(shù)計算，拼接，裁剪，植被蓋度計算，地表覆蓋因子計算等。植被蓋度因子（C）描述了植被對土壤風蝕的影響。C因子代表了農(nóng)作物、土地覆蓋、生物量、干擾土特性等因子，主要反映植被、作物對土壤侵蝕量的影響，C因子都小于或等于1，其相對變化范圍較其它因子大。在植被蓋度提取預處理工作完成以后，我們選擇歸一化差分植被指數(shù)（NDVI）反演植被蓋度的方法。首先需要計算歸一化植被指數(shù)，該指數(shù)是近紅外波段與紅光波段反射率比值的一種變換形式，即：（IV）具體的處理過程中會出現(xiàn)NDVI大于1和小于-1的情況，故在ERDAS中建模，CONDITIONAL{(ndvi<-1)-1,(ndvi>1)1,(-1=<ndvi<=1)ndvi)}，將大于1的數(shù)據(jù)改為1，小于-1的數(shù)據(jù)改為-1。然后，通過基于植被指數(shù)的方法來反演植被蓋度，其常用方法為——像元二分法，公式(V)計算植被蓋度。（V）式中，NDVImax和NDVImin為純植被與純土壤的植被指數(shù)；NDVI為待求象元的植被指數(shù)。在此，通過NDVI直方圖（可由ERDAS影像窗口中的attribute導出統(tǒng)計）確定植被整個生長季NDVI的最大值和最小值，分別為NDVImax和NDVImin。最大植被指數(shù)和最小植被指數(shù)的確定，非常關(guān)鍵，其精度將影響到研究區(qū)植被蓋度提取的準確性，在實際操作過程中，由于3~6和9~11月植被長勢變化較大，因而反映到植被蓋度反演圖像中的亮度最亮和亮度最低的位置以及范圍就會有所不同，為保證驗證精度，我們通過野外實地調(diào)查并拍照的方法，找到最亮像元和最暗像元所在位置，根據(jù)其當?shù)刂脖桓采w度的實際情況確定出NDVImax和NDVImin。模型運行要求每個月的植被蓋度數(shù)據(jù)，對于遙感數(shù)據(jù)資料不全的1976年、1990年、2000年和2006年植被蓋度因子的計算，采用2010年各月植被蓋度數(shù)據(jù)為參考，為不同地物類型計算蓋度平均值，并根據(jù)土地利用數(shù)據(jù)，為不同地物類型賦平均值作為該月份的研究區(qū)域的植被蓋度數(shù)據(jù)。通過前四期相應月份的植被蓋度數(shù)據(jù)與2010年相應月份植被蓋度數(shù)據(jù)的映射，得到前四期各月的植被蓋度數(shù)據(jù)。風蝕措施因子計算風蝕措施因子就是指為有效防止土壤風蝕，所采取的措施，主要分為兩大類人工措施和自然措施。風蝕措施包括工程措施、林草措施和保護性耕作措施；對于風沙地區(qū)工程措施主要采用機械沙障和化學固沙等方法；林草措施即生物措施，主要是通過造林、種草以及封山育林、育草等措施，提高植被覆蓋率；保護性耕作措施較常見的有帶狀間作、溝壟耕作、深耕免耕、留茬覆蓋、豆科牧草種植，如苜蓿等，從而改良土壤結(jié)構(gòu)，增強土壤抗風蝕能力。新疆兵團農(nóng)一師墾區(qū)地處塔克拉瑪干沙漠北緣，根據(jù)野外調(diào)查和搜集資料總結(jié)出新疆兵團農(nóng)一師墾區(qū)采取的防治土壤風蝕措施包括經(jīng)濟林、農(nóng)田防護林網(wǎng)、護路林、防風固沙林、沙漠公路兩側(cè)機械沙障、溝壟耕作措施、地膜覆蓋措施、帶狀間作措施、作物留茬、深耕和種植豆科牧草11種風蝕措施。參照國內(nèi)外相關(guān)文獻和研究，為每種風蝕措施類型確定其估算值，并給每種土地利用類型賦相應的風蝕措施值，從而將風蝕措施因子面域化。土壤風蝕模型運行圖2-2土壤風蝕模型運行方式針對修正后的土壤風蝕定量模型中不同模型參數(shù)特點，進行大量野外調(diào)查、取樣分析以及實地風蝕量監(jiān)測，并通過新疆兵團農(nóng)一師大風和沙塵暴天氣歷史資料分析，確立風蝕發(fā)生時段為3-6月和9-11月，從而為下一步模型修正結(jié)果的驗證、模型尺度修正系數(shù)的調(diào)整，以及各風蝕影響因子的提取奠定基礎。由于遙感數(shù)據(jù)資料的限制，模型運行步長設置為以月為時間尺度，三個子模型分別計算3-6月和9-11月的月土壤風蝕模數(shù)，累加得到年土壤風蝕模數(shù)，最后乘以各自模型區(qū)的面積，從而計算出整個農(nóng)墾區(qū)的年土壤風蝕量。土壤風蝕模型驗證本研究選取的土壤風蝕模型主要是對各類地表土壤風蝕模數(shù)的估算，修正后能夠預報新疆地區(qū)不同風力、不同植被狀況下的月土壤風蝕量，進而預報不同年份的土壤風蝕量。具體說，是將易風蝕地表劃分為耕地、林草地、沙地三大類，分別針對不同地表類型估算土壤風蝕量。為了對修正模型的適用性進行檢驗，通過野外調(diào)查通過數(shù)據(jù)監(jiān)測，研究區(qū)氣象、水文、土壤數(shù)據(jù)，用土壤風蝕模型計算出每月的土壤風蝕量，將研究區(qū)2010年土壤風蝕模型計算結(jié)果和2010年野外調(diào)查實際測量的月土壤風蝕量進行對比分析。通過模擬結(jié)果與2010年實際測量的月土壤風蝕量的對比分析，來驗證模型的準確度。土壤風蝕動態(tài)變化分析以五期多源遙感影像為數(shù)據(jù)源，提取研究區(qū)界線及種植類型，分析其變化規(guī)律，與提取的相應時期土壤侵蝕結(jié)果對比分析，分析其對土壤侵蝕的影響狀況，并結(jié)合野外實測數(shù)據(jù)與土壤風蝕模型運行結(jié)果數(shù)據(jù)，分析研究區(qū)因子動態(tài)變化、土壤風蝕模數(shù)和風蝕強度動態(tài)變化情況。墾區(qū)建設對土壤侵蝕影響研究結(jié)合野外實測數(shù)據(jù)與土壤風蝕模型運行結(jié)果數(shù)據(jù)，分析研究區(qū)5期墾區(qū)的面積變化與土壤風蝕的變化，總結(jié)墾區(qū)建設對土壤侵蝕影響變化規(guī)律。水土保持生態(tài)效益評價研究區(qū)水土保持及荒漠化防治概況兵團墾區(qū)水土保持的具體治理工作主要為防風治沙、防洪治水、改良鹽堿地三大方面。表現(xiàn)在以下幾方面：(1)防風治沙體系趨于完善，防護范圍不斷擴大；(2)防洪護岸，河道整治力度不斷加大；(3)改良鹽堿地技術(shù)提高，防止鹽堿化措施加強（資料來源于《新疆生產(chǎn)建設兵團水土保持規(guī)劃-報批稿》）。兵團自成立以來，在原有的荒漠戈壁上營造450萬hm2的綠洲生態(tài)系統(tǒng)，極大地改善了兵團墾區(qū)范圍內(nèi)水土流失狀況，綠洲生態(tài)系統(tǒng)的建成有效地阻止了風沙對人類的侵害，阻止洪水對人類的威脅，為兵團的建設發(fā)展營造了良好的生態(tài)環(huán)境。2000-2010年水土保持治理措施規(guī)劃中提出：新增治理措施面積87.04萬hm2，其中新建改造農(nóng)田防護林5.4萬hm2，新建改造防風固沙林15.6萬hm2，新建經(jīng)濟林2.64萬hm2，人工種草0.6萬hm2，改良草場3萬hm2，鹽堿地改良17.9萬hm2，退耕還林(草)5.2萬hm2，封禁36.02萬hm2，其它措施0.68萬hm2。工程措施有護岸工程200km，小型水保工程360處（資料來源于《新疆生產(chǎn)建設兵團水土保持規(guī)劃-報批稿》）。水土保持生態(tài)效益及評價的意義生態(tài)效益是水土流失治理綜合效益的重要組成部分，是生物種群能力、物質(zhì)轉(zhuǎn)化效率及維持生態(tài)環(huán)境穩(wěn)定能力的反映。水土保持生態(tài)效益是指水土保持對人類和生態(tài)環(huán)境在有序結(jié)構(gòu)維持和動態(tài)平衡保持方面輸出效益之和，包括控制土壤侵蝕、改良土壤、調(diào)節(jié)氣候，涵養(yǎng)水源，減少災害、保存物種、改善水資源環(huán)境條件等。水土流失生態(tài)效益綜合表現(xiàn)為：保持水土，改善生態(tài)環(huán)境，增強生態(tài)系統(tǒng)的穩(wěn)定性、承載力，最終實現(xiàn)系統(tǒng)內(nèi)經(jīng)濟、社會效益的增加。進行生態(tài)效益的綜合評價，可以考察水土流失治理及荒漠化防治綜合目標的實現(xiàn)程度。生態(tài)效益評價指標體系構(gòu)建1、生態(tài)效益評價指標選取原則（1）地理區(qū)域性差異原則；（2）主導因素原則；（3）具體指標的可操作性原則；（4）具體指標的不可替代性原則。根據(jù)研究區(qū)水土流失及荒漠化的發(fā)生條件、發(fā)展程度及治理情況，選取最能反映研究區(qū)水土保持及荒漠化防治的效益作用的指標。2、水土保持生態(tài)效益評價指標體系（1）耕地面積占區(qū)域總面積百分比耕地面積占區(qū)域總面積百分比=耕地（水澆地和灌溉水田）面積/研究區(qū)域總面積（2）林草覆蓋率林草覆蓋率=林草（林地、草地、園地）面積/研究區(qū)域總面積（3）鹽堿化面積百分比鹽堿地面積百分比=鹽堿地（重度鹽堿、中度鹽堿、輕度鹽堿）面積/研究區(qū)域總面積（4）沙化面積百分比沙化面積百分比=荒漠化土地（沙地、鹽堿地和裸巖石礫地）/研究區(qū)域總面積（5）沙塵暴日數(shù)風沙日包括浮塵日、揚沙日和沙塵暴日。是荒漠化治理后減少就地起沙、荒漠化程度減輕的顯著指標。和災害風降低率、土壤風蝕量減少率等相關(guān)因子相比，后者適用的尺度偏小，數(shù)據(jù)獲取困難。本報告采用沙塵暴日數(shù)作為評價指標。（數(shù)據(jù)來源于：18:81/scbweb/index.htm）（6）土壤風蝕模數(shù)土壤風蝕模數(shù)的大小直觀地反映了土壤的風蝕程度。土壤風蝕模數(shù)的大小可由耕地、林草地、沙地經(jīng)驗模型計算得到。（7）風蝕面積占區(qū)域總面積百分比風蝕面積占區(qū)域總面積百分比=受風蝕土地面積（處水庫、坑塘水面外的土地面積）/區(qū)域土地總面積表2-1水土保持生態(tài)效益評價指標體系目標層（A）制約層（B）指標層（C）水土保持生態(tài)效益（A）維持生態(tài)系統(tǒng)穩(wěn)定性（B1）耕地面積占區(qū)域總面積百分比（C11）林草覆蓋率（C12）改良土壤功能（B2）鹽堿化面積百分比（C21）沙化面積百分比（C22）防沙治沙功能（B3）沙塵暴日數(shù)（C31）土壤風蝕模數(shù)（C32）土壤風蝕面積占區(qū)域總面積百分比（C33）評價指標權(quán)重確定由于各指標因子在指標體系中的作用不同，對于生態(tài)環(huán)境狀況影響程度的有差異性，為了區(qū)分其對系統(tǒng)影響的差異性，常采用加權(quán)評價法。迄今為止，對權(quán)重的確定問題已經(jīng)進行了大量的研究，有以研究人員的實踐經(jīng)驗和主觀判斷為主來確定權(quán)重的，也有用各種數(shù)學方法來確定權(quán)重，本研究報告主要是采用層次分析方法（AHP）來確定縣級生態(tài)狀況的各評價參數(shù)的權(quán)重，咨詢有關(guān)專家意見，構(gòu)成判斷矩陣，獲得各層指標權(quán)數(shù)及隨機一致性率值。1、層次分析法概述層次分析法（AHP（AnalyticHierarchyProcess））：是由美國匹茲堡大學教授T.LSaaty最早提出的一種定性與定量分析結(jié)合的多目標評價決策方法。它將決策者對復雜系統(tǒng)的評價思維過程數(shù)學化，其基本思路是決策者通過把復雜問題分解為若干層次和要素，在各要素之間簡單地進行比較、判斷和計算，然后再轉(zhuǎn)為對這些元素的整體權(quán)重進行排序判斷，最后確立各元素的權(quán)重。該方法的原理和應用步驟如下：（1）構(gòu)建層次結(jié)構(gòu)模型構(gòu)建層次結(jié)構(gòu)模型之前，首先要對研究的對象進行系統(tǒng)分析，根據(jù)既定目標把所考慮的因素按其屬性層次化、條理化，構(gòu)造系統(tǒng)的層次模型。一般的決策層次分析模型可以分為三層，即最高層、中間層和最底層。最高層是目標層，代表要解決問題的預期目標或理想結(jié)果；中間層是準則層（也稱為條件層），是對上層目標的具體化和度量，對下層而言，則是用來評價各種備選方案的優(yōu)劣；最底層是方案層，是實現(xiàn)目標的方案、方法和途徑。（2）構(gòu)造判斷矩陣各個評價指標權(quán)重的確定在綜合評價中占有非常重要的位置，權(quán)重的大小對評估結(jié)果十分重要，它反映了各指標的相對重要性。其含義為:在縱向上，評價指標權(quán)重反映了該指標變化對生態(tài)綜合效益變化所起作用的大?。辉跈M向上，評價指標權(quán)重表示了該指標在同一評價指標層次中所處的地位，權(quán)重確定的合理與否將直接影響評判結(jié)果。依據(jù)層次分析法，請教有關(guān)從事該領(lǐng)域的專家，結(jié)合新疆兵團中現(xiàn)有的生態(tài)環(huán)境，構(gòu)造兩兩比較判斷矩陣。判斷矩陣中的標度值依據(jù)Saaty提出的1-9及其倒數(shù)作為衡量尺度的標度方法給出。表2-2層次分析法評價指標重要程度標度值及其含義標度含義1表示第i個因素與第j個因素的影響相同3表示第i個因素比第j個因素的影響稍強5表示第i個因素比第j個因素的影響強7表示第i個因素比第j個因素的影響明顯強9表示第i個因素比第j個因素的影響絕對地強2、4、6、8表示第i個因素相對于第j個因素的影響介于上述兩個相鄰等級之間倒數(shù)設因素i與因素j比較得判斷，則因素j與因素i的重要性之比為ɑji=1/ɑij（3）層次單排序?qū)哟螁闻判蛑饕嬎闩袛嗑仃嚨淖畲筇卣鞲捌鋵奶卣飨蛄亢吞卣鞲?，即計算滿足：(1)由矩陣理論可知道，式中為判斷矩陣唯一非零的最大特征根，為其所對應的特征向量，則式3-1可改寫為：(2)式（3-2）的解經(jīng)過歸一化處理后，即為同一層次各因素對于上一層次某一準則相對重要性的排序權(quán)重。a、權(quán)重的計算1）設一與判斷矩陣同階的正規(guī)化初始向量，例如：；2）對于3）令4）對于預先給定的精度，當成立時則停止計算，取，即為所求的特征向量。否則繼續(xù)進行第二步計算。最后取(3)b、計算一致性指標，檢查判斷思維的一致性當判斷矩陣A的階數(shù)為1或2時，判斷矩陣總是一致的。當判斷矩陣A的階數(shù)大于2時，計算一致性指標=(4)是一致性或判斷可靠性的度量，為了度量不同階的判斷矩陣是否具有滿意的一致性，還要引入判斷矩陣的平均隨機一致性指標值，對于1～12階判斷矩陣的值如下表所示。表2-3多階判斷矩陣平均隨機一致性指標RI值矩陣階數(shù)123456789101112000.580.9021.411.451.491.511.54計算一致性比例：當時，認為判斷矩陣的一致性是可以接受的。（4）層次總排序為了得到遞階層次結(jié)果中每層次中所有元素相對于總目標的相對權(quán)重，需把前一步計算的結(jié)果進行適當?shù)慕M合，已計算出總排序的相對權(quán)重，并進行層次和結(jié)構(gòu)一致性檢驗。因此，要由上至下逐層進行，最后得出最低層次元素，即決策方案優(yōu)先順序的相對權(quán)重和整個遞階層次模型的判斷一致性。2、研究區(qū)評價指標權(quán)重計算（1）構(gòu)建層次模型如上表2-1所示。（2）構(gòu)建判斷矩陣表2-4A-B判斷矩陣計算表后者前者維持生態(tài)系統(tǒng)穩(wěn)定性改良土壤功能防沙治沙功能一致性比維持生態(tài)系統(tǒng)穩(wěn)定性1.00000.33330.25000.0707改良土壤功能3.00001.00000.3333防沙治沙功能4.00003.00001.0000表2-5B1-C判斷矩陣計算表后者前者農(nóng)墾區(qū)面積占區(qū)域總面積百分比林草覆蓋面積一致性比耕地面積占區(qū)域總面積百分比1.00000.25000.0000林草覆蓋面積4.00001.0000表2-6B2-C判斷矩陣計算表后者前者風蝕期平均土壤濕度鹽堿地面積占區(qū)域總面積百分比荒漠化土地面積占區(qū)域總面積百分比一致性比鹽堿化面積占區(qū)域總面積百分比3.00001.00000.50000.0000沙化土地面積占區(qū)域總面積百分比3.00002.00001.0000表2-7B3-C判斷矩陣計算表后者前者沙塵暴日數(shù)土壤風蝕模數(shù)土壤風蝕面積占區(qū)域總面積百分比一致性比沙塵暴日數(shù)1.00000.20002.00000.0279土壤風蝕模數(shù)5.00001.00006.0000土壤風蝕面積占區(qū)域總面積百分比0.50000.16671.0000（3）權(quán)重計算表2-8水土保持生態(tài)效益評價指標權(quán)重匯總表目標層（A）制約層（B）指標層（C）水土保持生態(tài)效益（A）維持生態(tài)系統(tǒng)穩(wěn)定性（B1）(0.1172)耕地面積占區(qū)域總面積百分比（C11）(0.0234)林草覆蓋率（C12）(0.0938)改良土壤功能（B2）(0.2684)鹽堿地面積百分比（C21）(0.1342)沙化面積百分比（C22）(0.1342)防沙治沙功能（B3）(0.6144)沙塵暴日數(shù)（C31）(0.1058)土壤風蝕模數(shù)（C32）(0.446)土壤風蝕面積占區(qū)域總面積百分比（C33）(0.0627)評價指標數(shù)據(jù)標準化由于本文選取的評價指標原始數(shù)據(jù)的量綱不同，即具體的數(shù)據(jù)單位不同，無法進行比較，必須對原始數(shù)據(jù)進行標準化處理。數(shù)據(jù)標準化的方法有多種，在本文中采用極差化處理。具體計算公式如下：式中：、分別為觀測值的最大值和最小值。評價指標中鹽堿地面積百分比、荒漠化面積百分比和沙塵暴日數(shù)、土壤風蝕模數(shù)、和土壤風蝕面積占區(qū)域總面積百分比都為負向指標，因此在進行同量綱化處理之前需對其進行同向化處理，本報告采用的是倒數(shù)法。綜合評價模型選擇由單項評價值計算綜合評價值的方法及綜合評價的方法有很多，常用的有計分法、指數(shù)法、功效系數(shù)法等，本報告采用加權(quán)指數(shù)法計算：其中，為單項指標權(quán)重；為標準化后的單項指標值。評價等級劃分1、等級劃分結(jié)果將標準化后的指標作用分值乘以對應的權(quán)重，逐級分層次計算，最后可得到水土保持生態(tài)效益的綜合指數(shù)得分。綜合指數(shù)得分多少與水土保持及荒漠化治理效益水平呈正相關(guān)，得分越高，治理效益越好。根據(jù)研究，可以把效益分析結(jié)果可以劃分為5個不同的等級及其相對應的綜合指數(shù)見表。表2-9水土保持生態(tài)效益等級劃分表水土保持生態(tài)效益等級劃分綜合指數(shù)得分優(yōu)0.8-1良0.6-0.8中0.4-0.6較差0.2-0.4差0-0.22、等級說明對于已經(jīng)劃分的評價等級結(jié)果，其不同的等級反映了不同的荒漠化治理效益。表2-10水土保持生態(tài)效益分級級別優(yōu)良中較差差0.8-10.6-0.80.4-0.60.2-0.40-0.2狀態(tài)和治理前相比，生態(tài)系統(tǒng)已經(jīng)穩(wěn)定，鹽堿化和荒漠化土地基本消失。風蝕面積小、強度弱。和治理前相比，生態(tài)系統(tǒng)基本穩(wěn)定，鹽堿化和荒漠化土地大部分消失，風蝕面積相對較小、強度較弱。和治理前相比，生態(tài)系統(tǒng)極趨于穩(wěn)定，荒漠化和鹽堿化土地得到有效治理，風蝕面積減小、強度減弱。和治理前相比，生態(tài)系統(tǒng)尚不穩(wěn)定，土地荒漠化和鹽堿化較為嚴重，風蝕面積較大、強度較強。