漢瑞英，趙志平，肖能文，史娜娜，張風(fēng)春，高曉奇，劉高慧

(中國(guó)環(huán)境科學(xué)研究院，北京 100012)

0 引 言

被喻為“地球癌癥”的荒漠化是當(dāng)今世界關(guān)注的十大環(huán)境問(wèn)題之一[1-2]，中國(guó)的荒漠化形勢(shì)十分嚴(yán)峻[3]。根據(jù)第五次全國(guó)荒漠化和沙化土地監(jiān)測(cè)結(jié)果，全國(guó)荒漠化面積為261.2×104km2，占國(guó)土面積的27.2%。1992年聯(lián)合國(guó)環(huán)境與發(fā)展大會(huì)后，根據(jù)聯(lián)合國(guó)決議，部分國(guó)家政府形成《聯(lián)合國(guó)關(guān)于在發(fā)生嚴(yán)重干旱和/或荒漠化的國(guó)家特別是在非洲防治荒漠化的公約》(簡(jiǎn)稱(chēng)《聯(lián)合國(guó)防治荒漠化公約》)。中國(guó)作為第一批簽署該公約的國(guó)家之一，非常重視荒漠化的防治。黨的十九大報(bào)告亦明確提出“開(kāi)展國(guó)土綠化行動(dòng)，推進(jìn)荒漠化、石漠化、水土流失綜合治理”[4]。土地沙化是中國(guó)西北地區(qū)荒漠化的主要類(lèi)型之一[5]。

目前國(guó)內(nèi)外學(xué)者針對(duì)沙化動(dòng)態(tài)變化、沙化風(fēng)險(xiǎn)評(píng)估進(jìn)行了大量研究工作。早期的沙化動(dòng)態(tài)變化研究?jī)H僅是進(jìn)行局地沙化監(jiān)測(cè)[6]。近年來(lái)隨著遙感技術(shù)的發(fā)展，遙感數(shù)據(jù)[7]和沙丘活動(dòng)指數(shù)[8-9]被廣泛應(yīng)用于土地沙化動(dòng)態(tài)變化研究中。沙化風(fēng)險(xiǎn)研究多采用選取的沙漠化敏感性指標(biāo)進(jìn)行分級(jí)評(píng)價(jià)[10-11]，并基于評(píng)價(jià)結(jié)果來(lái)指示區(qū)域環(huán)境質(zhì)量，因此,識(shí)別沙源區(qū)域?qū)τ谏郴芯渴种匾?。Rayegan等利用陸地衛(wèi)星Landsat8遙感數(shù)據(jù)得到植被覆蓋、土壤水分和土地利用風(fēng)侵敏感性圖，進(jìn)而獲取沙源地潛在空間分布[12]。伊塔娜通過(guò)對(duì)比分析3期TM遙感影像，提取內(nèi)蒙古烏梁素海地區(qū)土地沙化信息[13]。侯閣通過(guò)對(duì)1984～2013年NDVI遙感影像分析，并結(jié)合景觀指數(shù)分析方法，探究了羌塘地區(qū)高寒草地沙化的時(shí)空動(dòng)態(tài)變化[14]。但上述研究都沒(méi)有分析沙源擴(kuò)散路徑，無(wú)法揭示沙化動(dòng)態(tài)變化趨勢(shì)，目前國(guó)內(nèi)外仍缺乏大尺度沙源擴(kuò)散路徑的相關(guān)研究。

為了更好地揭示宏觀尺度沙源擴(kuò)散變化，本文應(yīng)用“源-匯”景觀生態(tài)學(xué)理論，將沙源擴(kuò)散動(dòng)態(tài)過(guò)程融入到靜態(tài)格局中。2006年，陳利頂?shù)仁状翁岢觥霸?匯”景觀生態(tài)學(xué)理論[15]?！霸础笔鞘挛锵蛲鈹U(kuò)散的起點(diǎn)與物質(zhì)基礎(chǔ)，是能促進(jìn)生態(tài)過(guò)程發(fā)展的景觀類(lèi)型[16]；“匯”是阻止延緩生態(tài)過(guò)程發(fā)展的景觀類(lèi)型[17-18]。目前，“源-匯”理論已廣泛應(yīng)用于生態(tài)網(wǎng)絡(luò)構(gòu)建、環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)、景觀生態(tài)安全格局構(gòu)建等領(lǐng)域[19-21]。本文以新疆和田地區(qū)為研究對(duì)象，綜合考慮土壤風(fēng)蝕、土壤類(lèi)型、濕潤(rùn)指數(shù)、土壤侵蝕等因子，通過(guò)生態(tài)阻力面(Urban Expansion Ecological Resistance，UEER)模型提取沙源地，利用ArcGIS10.2軟件Cost Distance工具提取研究區(qū)潛在沙源擴(kuò)散路徑，同時(shí)采用重力模型篩選出沙源擴(kuò)散的重要路徑，分析沙源變化擴(kuò)散過(guò)程，厘清和田地區(qū)沙源擴(kuò)散規(guī)律，識(shí)別沙源擴(kuò)散節(jié)點(diǎn)，并提出在Ⅰ類(lèi)沙源擴(kuò)散節(jié)點(diǎn)建設(shè)綠洲防護(hù)區(qū)來(lái)阻礙沙源地的流動(dòng)擴(kuò)散，以期為沙源動(dòng)態(tài)變化研究提供一種新的思路，并為區(qū)域沙化防治提供參考。

1 研究區(qū)概況與數(shù)據(jù)來(lái)源

1.1 研究區(qū)概況

和田地區(qū)位于新疆塔里木盆地西南部，地理坐標(biāo)范圍為34°22′N(xiāo)～39°38′N(xiāo)和77°31′E～84°55′E，包含皮山縣、墨玉縣、和田縣、洛浦縣、策勒縣、于田縣、民豐縣及和田市(圖1)。研究區(qū)山勢(shì)起伏大，南依昆侖山和喀喇昆侖山，北臨塔克拉瑪干沙漠，大部分市縣處于被沙漠包圍和半包圍的狀態(tài)，是典型的內(nèi)陸干旱區(qū)。研究區(qū)位于中亞干旱區(qū)腹地，常年干旱少雨[22]，地表土質(zhì)疏松，流動(dòng)沙丘不斷侵占綠洲、草地、耕地[23]。年降水量為28.9～47.1 mm，年潛在蒸發(fā)量為2 198～2 790 mm[24]，日照時(shí)數(shù)為2 470～3 000 h，光照充足，無(wú)霜凍期。在氣候變化和人類(lèi)活動(dòng)的雙重影響下，該地區(qū)土地沙化嚴(yán)重，生態(tài)安全受到威脅。

1.2 數(shù)據(jù)來(lái)源

圖1 新疆和田地區(qū)土地利用現(xiàn)狀Fig.1 Present Situation of Land Use in Hotan Area of Xinjiang

2016年土地利用數(shù)據(jù)(空間分辨率為30 m×30 m)和濕潤(rùn)指數(shù)數(shù)據(jù)(空間分辨率為1 km×1 km)來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心；NDVI遙感數(shù)據(jù)采用中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)鏡像網(wǎng)站(http:∥www.gscloud.cn)提供的2016年逐月數(shù)據(jù)，利用最大值合成法獲取年NDVI遙感數(shù)據(jù)(空間分辨率為500 m×500 m)；土壤濕度因子(2016年)、日均風(fēng)速、降水、溫度和日照時(shí)數(shù)等來(lái)源于中國(guó)氣象科學(xué)數(shù)據(jù)共享網(wǎng)(http:∥cdc.cma.gov.cn)，用于計(jì)算風(fēng)蝕因子；土壤類(lèi)型和土壤理化性狀、雪覆蓋因子數(shù)據(jù)來(lái)源于中國(guó)西部環(huán)境與生態(tài)科學(xué)數(shù)據(jù)中心網(wǎng)站(http:∥westdc.westgis.ac.cn)，分辨率為1 km。

2 分析方法

2.1 基于生態(tài)阻力面模型的沙源潛在擴(kuò)散途徑模擬

最小阻力模型(Minimal Cumulative Resistance，MCR)被廣泛應(yīng)用于城市建設(shè)、保護(hù)地規(guī)劃等景觀格局分析以及生態(tài)網(wǎng)絡(luò)構(gòu)建等領(lǐng)域[25-29]。本文采用經(jīng)最小阻力模型改進(jìn)后的生態(tài)阻力面模型，即引入不同等級(jí)源的相對(duì)阻力因子(K)[30]，以沙源地為“源”，以阻礙沙源擴(kuò)散為“匯”，構(gòu)建適合于新疆和田地區(qū)的沙源擴(kuò)散模擬模型。其表達(dá)式為

(1)

式中：FUEER表示沙源擴(kuò)散最小生態(tài)累積阻力面值；f函數(shù)表示最小累積阻力與生態(tài)適宜性的負(fù)相關(guān)關(guān)系；min表示某景觀單元對(duì)不同的沙源取累積阻力最小值；Dij表示從源j到柵格單元i的空間距離；Ri表示柵格單元i對(duì)運(yùn)動(dòng)過(guò)程的阻力系數(shù)；Kj表示源j所屬等級(jí)的阻力因子，沙源風(fēng)險(xiǎn)等級(jí)越高，擴(kuò)展能力越強(qiáng)，相對(duì)阻力因子就越小。

2.1.1 沙源地提取方法

研究區(qū)起沙風(fēng)險(xiǎn)評(píng)價(jià)綜合考慮區(qū)域土壤風(fēng)蝕、土壤侵蝕、土壤類(lèi)型、濕潤(rùn)指數(shù)4個(gè)因子，基于經(jīng)典層次分析(AHP)法確定評(píng)價(jià)因子權(quán)重，通過(guò)Arc-GIS10.2軟件空間疊加運(yùn)算生成起沙風(fēng)險(xiǎn)評(píng)價(jià)圖層，并采用自然斷點(diǎn)法將其分為低度、較低、中度、較高和高度5個(gè)擴(kuò)散風(fēng)險(xiǎn)等級(jí)。根據(jù)風(fēng)險(xiǎn)等級(jí)劃分結(jié)果及區(qū)域生態(tài)狀況，提取高度風(fēng)險(xiǎn)等級(jí)區(qū)域作為沙源地。對(duì)確定權(quán)重的單因子進(jìn)行一致性檢驗(yàn)。當(dāng)一致性比例(CR)小于0.1時(shí)，表示判斷矩陣一致性較好，反之則需要調(diào)整。其具體表達(dá)式為

(2)

(3)

式中：n為判斷矩陣階數(shù)；λmax為判斷矩陣特征值；RI為平均隨機(jī)一致性指標(biāo)，當(dāng)階數(shù)為4時(shí)，RI=0.94；CI為一致性指標(biāo)。

(1)利用風(fēng)蝕模型(RWEQ)開(kāi)展土壤風(fēng)蝕評(píng)估。風(fēng)蝕模型充分考慮了氣候、土壤可蝕性、土壤結(jié)皮、地表粗糙度及植被覆蓋狀況等，其因子覆蓋全面。通過(guò)修正的風(fēng)蝕方程計(jì)算多因素影響下的區(qū)域土壤轉(zhuǎn)運(yùn)總量，可作為該區(qū)域的土壤風(fēng)蝕量[31]。在理想數(shù)據(jù)情況下，風(fēng)蝕模型能很好地完成對(duì)土壤侵蝕情況的預(yù)測(cè)，計(jì)算公式[32]為

Qmax=109.8WFEFSCFK′C

(4)

S=150.71(WFEFSCFK′C)-0.371

(5)

(6)

式中：Qmax表示風(fēng)力的最大輸沙能力；WF表示氣象因子；EF表示土壤可蝕性因子；SCF表示土壤結(jié)皮因子；K′表示土壤糙度因子；C表示植被因子；S為關(guān)鍵地塊長(zhǎng)度；SL為土壤損失量；X表示距上風(fēng)向距離，本次計(jì)算取50 m。

土壤風(fēng)蝕評(píng)估結(jié)果如圖2(a)所示。

氣象因子計(jì)算公式為

(7)

式中：WS2為2 m處風(fēng)速；WSt為2 m處臨界風(fēng)速(假定為5 m·s-1)；N為風(fēng)速的觀測(cè)次數(shù)(一般為500次)；Nd為試驗(yàn)的天數(shù)；ρ為空氣密度；g為重力加速度；SW為土壤濕度因子,無(wú)量綱；SD為雪覆蓋因子。

土壤可蝕性因子計(jì)算公式為

式中：Sa為土壤砂粒含量(質(zhì)量分?jǐn)?shù)，下同)；Si為土壤粉砂含量；Cl為黏土含量；OM為有機(jī)質(zhì)含量；w(CaCO3)為碳酸鈣含量。

土壤結(jié)皮因子計(jì)算公式為

(9)

式中：Cl為5.0%～39.3%;OM為0.32%～4.74%。

土壤糙度因子計(jì)算公式為

(10)

(11)

式中：Crr為隨機(jī)粗糙度，取值為0；Kr為地形粗糙長(zhǎng)度；ΔH為距離L范圍內(nèi)的海拔高程差。

(2)土壤侵蝕評(píng)估。區(qū)域土地特征信息評(píng)價(jià)主要參考全國(guó)生態(tài)系統(tǒng)格局和全國(guó)1∶1 000 000數(shù)字化土壤圖中的數(shù)據(jù)，基于水利部發(fā)布的《土壤侵蝕分類(lèi)分級(jí)標(biāo)準(zhǔn)》(SL 190—2007)[33]，將土壤侵蝕等級(jí)分類(lèi)合并為輕度、中度、較強(qiáng)、強(qiáng)度與極強(qiáng)5個(gè)等級(jí)[圖2(b)]。

(3)土壤類(lèi)型評(píng)估?；哪苯幼饔玫膶?duì)象是土壤，土壤類(lèi)型不僅是土壤抗風(fēng)蝕能力的綜合體現(xiàn)，而且是影響土地風(fēng)蝕的重要因素[13]。不同類(lèi)型的土壤及其母質(zhì)對(duì)于各種侵蝕營(yíng)力的抵抗作用不同。根據(jù)研究區(qū)土壤類(lèi)型，利用ArcGIS10.2軟件工具箱Reclassify功能，將不同類(lèi)型的土壤劃分為5類(lèi)[圖2(c)]起沙等級(jí)：林地、冰川和建設(shè)用地劃為不起風(fēng)沙區(qū)；草甸土和灌淤土劃為不易起風(fēng)沙區(qū)；石質(zhì)土和棕漠土劃為較易起風(fēng)沙區(qū)；新積土和粗骨土劃為易起風(fēng)沙區(qū)；風(fēng)沙土劃為極易起風(fēng)沙區(qū)。

(4)濕潤(rùn)指數(shù)。濕潤(rùn)指數(shù)適用于旬以上尺度干濕研究[34]，能很好地反映和田地區(qū)氣候干濕狀況變化特征[35][圖2(d)]。

區(qū)域環(huán)境因子和區(qū)域生態(tài)背景共同控制著沙源擴(kuò)散，因此，確定各因子對(duì)沙源擴(kuò)散的貢獻(xiàn)在沙源風(fēng)沙擴(kuò)散等級(jí)評(píng)價(jià)中尤其重要。應(yīng)用經(jīng)典層次分析法確定各指標(biāo)權(quán)重(表1)，計(jì)算判斷矩陣一次性比例，得出CR=0.018，該值小于0.1，符合判斷矩陣一致性檢驗(yàn)結(jié)果。

表1 判斷矩陣Tab.1 Judgement Matrix

2.1.2 生態(tài)阻力面構(gòu)建方法

考慮不同土地利用景觀類(lèi)型的景觀阻力、植被覆蓋情況，經(jīng)專(zhuān)家打分確定不同生境斑塊的阻力大小(表2)。

（5）砼工程中砼標(biāo)號(hào)分多種強(qiáng)度等級(jí)和不同的配合比，一般情況應(yīng)該把砼配合比單獨(dú)歸納到一章內(nèi)，套定額時(shí)，再根據(jù)不同標(biāo)號(hào)在本章節(jié)中查找換算。但考慮到系統(tǒng)內(nèi)填寫(xiě)檢修計(jì)劃委托書(shū)的是設(shè)備點(diǎn)檢員，定額編制組在編制砼澆筑子目時(shí)，直接將廠內(nèi)常用的強(qiáng)度等級(jí)混凝土按照不同標(biāo)號(hào)編制在砼澆筑子目?jī)?nèi)。這樣，使用者就不用選定子目后再進(jìn)行查找、換算操作了。

圖2 土壤風(fēng)蝕、土壤侵蝕、土壤類(lèi)型和濕潤(rùn)指數(shù)分級(jí)評(píng)價(jià)等級(jí)分布Fig.2 Grade Distributions for Grading Evaluation of Soil Drifting, Soil Erosion, Soil Type and Moisture Index

沙源擴(kuò)散必須克服阻力。不同等級(jí)的沙源擴(kuò)張能力不同，如沙漠腹地的沙源擴(kuò)張速度明顯快于覆有灌叢、草地的沙源，風(fēng)蝕模型考慮了沙源地等級(jí)。基于ArcGIS10.2軟件柵格計(jì)算工具，分別生成一級(jí)源地、二級(jí)源地、三級(jí)源地的生態(tài)阻力面，通過(guò)對(duì)3幅生態(tài)阻力面圖層進(jìn)行疊加分析，計(jì)算各柵格單元的最小值，最終生成沙源擴(kuò)散生態(tài)阻力面(圖3)。

2.2 基于重力模型的沙源重要擴(kuò)散路徑提取

重力模型普遍應(yīng)用于生態(tài)網(wǎng)絡(luò)中重要生態(tài)廊道的提取[36-37]，生態(tài)廊道的重要性可用生態(tài)斑塊之間的相互作用力大小進(jìn)行度量[38]。在已構(gòu)建的沙源擴(kuò)散網(wǎng)絡(luò)基礎(chǔ)上，基于重力模型計(jì)算出沙源斑塊之間相互作用力大小，結(jié)合沙源地等級(jí)，定量評(píng)價(jià)每條擴(kuò)散路徑的相對(duì)重要程度。重力模型計(jì)算公式為

(12)

式中：Gab為斑塊a和斑塊b的相互作用力；Na、Nb分別表示斑塊a和斑塊b的權(quán)重值；Dab為斑塊a和斑塊b之間擴(kuò)散路徑阻力的標(biāo)準(zhǔn)化值；Pa為斑塊a的阻力值；Pb為斑塊b的阻力值；Sa為斑塊a的面積；Sb為斑塊b的面積；Lab為斑塊a與斑塊b之間的生態(tài)廊道累積阻力值；Lmax為研究區(qū)所有擴(kuò)散路徑阻力最大值。

基于ArcGIS10.2軟件空間提取工具計(jì)算沙源斑塊之間累積阻力值，采用重力模型計(jì)算沙源斑塊之間的相互作用力大小，選取相互作用力較大的沙源擴(kuò)散路徑。

3 結(jié)果分析

3.1 沙源地劃定

表2 生態(tài)阻力面賦值及權(quán)重Tab.2 Values and Weights of Ecological Resistance Surface

圖3 沙源擴(kuò)散生態(tài)阻力面Fig.3 Ecological Resistance Surface of Sand Source Diffusion

圖4 沙源擴(kuò)散等級(jí)評(píng)價(jià)及沙源地分級(jí)圖Fig.4 Maps of Grading Evaluation of Sand Source Diffusion and Classification of Sand Source Area

3.2 提取沙源擴(kuò)散路徑

基于沙源擴(kuò)散生態(tài)阻力面和提取的沙源地，采用最小費(fèi)用路徑法，生成沙源斑塊之間的沙源擴(kuò)散路徑，得到和田地區(qū)最小生態(tài)功能耗費(fèi)梯度數(shù)據(jù)，形成由多條沙源潛在擴(kuò)散路徑組成的沙源擴(kuò)散網(wǎng)絡(luò)。研究區(qū)共有128條沙源潛在擴(kuò)散路徑[圖5(a)]，路徑密度較高的區(qū)域集中分布于和田地區(qū)西北部的墨玉縣、皮山縣、于田縣中部、民豐縣南部。于田縣到民豐縣沙源擴(kuò)散阻力較大，擴(kuò)散路徑僅兩條，主要是因?yàn)樵搮^(qū)域的林地、草地阻礙了沙源擴(kuò)散。對(duì)生成的沙源擴(kuò)散路徑設(shè)置閾值，提取最小生態(tài)功能耗費(fèi)路徑，剔除其中經(jīng)過(guò)同一斑塊的冗余路徑，最后生成斑塊間相互作用較強(qiáng)的最優(yōu)路徑，即沙源擴(kuò)散重要路徑，共71條[圖5(b)]。

圖5 沙源擴(kuò)散路徑空間分布Fig.5 Spatial Distributions of Diffusion Paths of Sand Sources

研究區(qū)西北部沙源擴(kuò)散作用最強(qiáng)。該區(qū)域位于塔克拉瑪干沙漠腹地，沙丘植被稀少，沙源主要以流動(dòng)沙丘為主。中部沙源等級(jí)較低，沙源擴(kuò)散路徑較少，沙化主要是流沙侵蝕綠洲和耕地而形成，沙源擴(kuò)散作用較弱。東南部沙源擴(kuò)散作用較集中，主要沙化形式是流沙侵占裸土地、裸巖石地和稀疏草地。

3.3 沙源擴(kuò)散節(jié)點(diǎn)識(shí)別

在構(gòu)建沙源擴(kuò)散網(wǎng)絡(luò)時(shí)，除了面積較大的沙源斑塊，沙源擴(kuò)散路徑之間的交點(diǎn)、路徑和最小阻力路徑的交點(diǎn)也是沙源擴(kuò)散匯集的地方，即沙源擴(kuò)散節(jié)點(diǎn)。基于和田地區(qū)沙源擴(kuò)散的重要路徑識(shí)別，研究區(qū)沙源擴(kuò)散節(jié)點(diǎn)共有48個(gè)(圖6)。為有效阻礙沙源在節(jié)點(diǎn)區(qū)域匯集，將沙源擴(kuò)散節(jié)點(diǎn)劃分為Ⅰ類(lèi)和Ⅱ類(lèi)沙源擴(kuò)散節(jié)點(diǎn)。Ⅰ類(lèi)沙源擴(kuò)散節(jié)點(diǎn)31個(gè)，主要分布在有植被覆蓋的綠洲、荒漠交錯(cuò)區(qū)域，這些區(qū)域也是防沙治沙的關(guān)鍵區(qū)域；Ⅱ類(lèi)沙源擴(kuò)散節(jié)點(diǎn)17個(gè)，主要分布在沒(méi)有植被覆蓋的沙地或裸土地，沙源流動(dòng)活躍。因此，建議將識(shí)別的31個(gè)Ⅰ類(lèi)沙源擴(kuò)散節(jié)點(diǎn)建設(shè)為綠洲防護(hù)區(qū)，可以阻礙沙源地之間的交流，防止沙漠對(duì)綠洲和農(nóng)田的侵蝕。

圖6 沙源擴(kuò)散節(jié)點(diǎn)分布Fig.6 Distribution of Diffusion Nodes of Sand Sources

4 結(jié) 語(yǔ)

在綠色“一帶一路”建設(shè)背景下，荒漠化防治是保障絲綢之路生態(tài)安全和可持續(xù)發(fā)展的重要措施[39]，沙源擴(kuò)散路徑的識(shí)別可以為荒漠化防治提供有效方法。本文嘗試將“源-匯”景觀生態(tài)學(xué)理論應(yīng)用于沙源擴(kuò)散過(guò)程，通過(guò)生態(tài)阻力面模型構(gòu)建新疆和田地區(qū)沙源擴(kuò)散路徑，同時(shí)運(yùn)用重力模型獲取沙源擴(kuò)散的重要路徑，并識(shí)別沙源擴(kuò)散節(jié)點(diǎn)，厘清了和田地區(qū)沙源擴(kuò)散規(guī)律，將沙化動(dòng)態(tài)過(guò)程融入到荒漠靜態(tài)景觀。不同于生態(tài)流的擴(kuò)散，沙源擴(kuò)散阻力較大的地區(qū)反而有利于生態(tài)流擴(kuò)散，本文將提取的沙源地分級(jí)，在構(gòu)建阻力面時(shí)引入不同等級(jí)源地的相對(duì)阻力因子，綜合考慮土地利用、植被覆蓋和生態(tài)等重要性因素，最終生成沙源擴(kuò)散生態(tài)阻力面。根據(jù)生態(tài)阻力面模型結(jié)果，一級(jí)源地?cái)U(kuò)散作用最強(qiáng)，二級(jí)源地次之，三級(jí)源地?cái)U(kuò)散作用最弱，但從荒漠化防治角度上看，三級(jí)源地主要分布在綠洲、荒漠交錯(cuò)區(qū)域，是防沙治沙最關(guān)鍵的區(qū)域，而分布在沙漠腹地的沙源雖然流動(dòng)多活躍，但不是防治沙化優(yōu)先區(qū)。

(1)和田地區(qū)一級(jí)源地占提取沙源總面積的73.57%，面積占比最大，集中分布在研究區(qū)西北部塔克拉瑪干沙漠腹地和民豐縣南部，是沙源擴(kuò)散過(guò)程的主要“源”物質(zhì)；二級(jí)源地占沙源總面積的17.35%，主要分布在于田縣中部、皮山縣中部；三級(jí)源地占沙源總面積的9.07%，零星分散在綠洲、荒漠交錯(cuò)區(qū)域，是沙化防治重點(diǎn)關(guān)注區(qū)域。

(2)和田地區(qū)西北部沙源擴(kuò)散作用最強(qiáng)，東南部沙源擴(kuò)散作用較集中，中部沙源等級(jí)較低，沙源擴(kuò)散作用最弱。

(3)和田地區(qū)共有71條沙源擴(kuò)散重要路徑，31個(gè)Ⅰ類(lèi)沙源擴(kuò)散節(jié)點(diǎn)，建議將其建設(shè)為綠洲防護(hù)區(qū)，可以阻礙沙源地之間的交流，防止沙漠對(duì)綠洲和農(nóng)田的侵蝕。

(4)沙源地阻力面的生成受諸多因素影響，很難準(zhǔn)確量化。由于本研究數(shù)據(jù)的可獲得性，沙源擴(kuò)散生態(tài)阻力面構(gòu)建未能考慮荒漠植被防風(fēng)固沙能力、地表徑流量等因素，劃定的沙源地也沒(méi)有進(jìn)行現(xiàn)場(chǎng)調(diào)研修正。下一步將在本研究基礎(chǔ)上進(jìn)行沙源擴(kuò)散生態(tài)阻力面的不可替代性評(píng)估，同時(shí)識(shí)別沙源類(lèi)型，并基于不同類(lèi)型的沙源地和沙源擴(kuò)散路徑實(shí)施相應(yīng)的荒漠區(qū)域沙化防治措施。