王一凱，張蒙蒙，段衛(wèi)虎，楊凱悅，周汝良，葉江霞

(西南林業(yè)大學(xué)，云南 昆明 650224)

基于DOMAIN模型的干熱河谷稀樹灌木草叢空間分布

王一凱，張蒙蒙，段衛(wèi)虎，楊凱悅，周汝良，葉江霞

(西南林業(yè)大學(xué)，云南 昆明 650224)

以溫度、DEM(數(shù)字高程模型)、NDVI(歸一化植被指數(shù))3個因子為參數(shù)，運用DOMAIN模型對元江縣干熱河谷稀樹灌木草叢空間分布進行模擬，實現(xiàn)了其空間分布的可視化表達。模擬結(jié)果表明：稀樹灌木草叢群落主要分布于除元江縣城外的河谷兩側(cè)，其中北部分布較為集中，南部較為分散。在自然條件及人為因素雙重干擾下，元江縣干熱區(qū)域的面積呈現(xiàn)出擴大趨勢，尤其是元江北部龍?zhí)多l(xiāng)以北、青龍廠鎮(zhèn)以南的地區(qū)存在更高的演化風(fēng)險，建議將上述區(qū)域及元江國家級自然保護區(qū)章巴片區(qū)列為重點森林保護區(qū)域。

干熱河谷；稀樹灌木草叢；DOMAIN模型；空間分布；元江縣

干熱河谷在干旱區(qū)系劃分中屬于局部干旱生境，稀樹灌木草叢群落為干熱河谷的主要植被類型[1]。群落主要以草叢為主，灌木和喬木散生。主要分布在云南怒江、金沙江、瀾滄江、元江等大河及其支流：綠汁江、普渡河、南盤江、勐河、龍川江、漁泡江等深度切割、地形比較封閉的地段[2]，組成我國特殊的干熱河谷生態(tài)環(huán)境類型區(qū)。干熱河谷的出現(xiàn)和分布具有明顯的自然背景，標示著自然地帶性環(huán)境的特殊性。

干熱河谷是典型的生態(tài)脆弱帶，生態(tài)系統(tǒng)退化嚴重，植被稀疏，森林大面積萎縮，系統(tǒng)結(jié)構(gòu)不穩(wěn)定，生態(tài)功能低下，土壤侵蝕劇烈，造林極端困難地段和生態(tài)綜合整治的難點區(qū)域。在過往的研究中，大多集中于土壤水分[3-4]、土壤組成[5-7]、植被組成[8-9]等方面及與水熱條件(降水、蒸發(fā)量)的關(guān)系；同時在其成因分析方面，何永彬等[10-11]把干熱河谷的成因總結(jié)為氣候、植被、土地結(jié)構(gòu)、人為干擾4個原因；高文學(xué)等[12]、柴宗新等[13]對該區(qū)域生態(tài)恢復(fù)等進行研究，提出一系列植被恢復(fù)的方法，例如：建立水土保持經(jīng)濟林、封山育林、加強溝谷植被覆蓋度等措施，但沒有對該區(qū)域主要的生態(tài)類型進行空間分布的模擬研究。

從以前學(xué)者研究中得知，稀樹灌木草叢的形成與環(huán)境條件有著密切的關(guān)系，特別是降水和蒸發(fā)量。本研究從環(huán)境梯度因子方面對其進行空間分布的模擬，以元江縣為研究區(qū)域，使用海拔、溫度、歸一化植被指數(shù)(NDVI)來模擬其空間分布。以期能對干熱河谷區(qū)域的空間動態(tài)變化提供可靠依據(jù)，以便更準確地預(yù)測稀樹灌草叢的發(fā)展趨勢，從而對生物多樣性的保護、生態(tài)治理提供重要的實踐指導(dǎo)。

1 材料與方法

1.1 研究數(shù)據(jù)

1)溫度場數(shù)據(jù)。采用MODIS11級產(chǎn)品MOD11數(shù)據(jù)，主要為地表溫度和輻射率，空間分辨率1 km。應(yīng)用MODIS Reprojection Tools (MRT) 軟件對下載的數(shù)據(jù)進行數(shù)據(jù)格式轉(zhuǎn)換、數(shù)據(jù)拼接與裁剪。同時使用云南省134個氣象站點15 a的年均數(shù)據(jù)進行訂正(訂正模型通過置信度為95%的統(tǒng)計檢驗，其中回歸方程F值為108.587)，得到云南省溫度基線數(shù)據(jù)，綜合考慮后續(xù)數(shù)據(jù)處理與現(xiàn)有數(shù)據(jù)的條件，整個研究數(shù)據(jù)空間分辨率選為100 m，同時為后續(xù)的空間分布模型預(yù)測，統(tǒng)一采用地理坐標系Beijing 1954，投影系統(tǒng)為：Beijing 1954 3 degree Degree_GK_CM_102E。最后通過ArcGIS的掩膜提取功能裁剪，得到研究區(qū)的溫度基線數(shù)據(jù)。

2)NDVI數(shù)據(jù)。NDVI數(shù)據(jù)來自寒旱區(qū)科學(xué)數(shù)據(jù)中心(http：//westdc.westgis.ac.cn)，通過ArcGIS軟件的柵格計算功能，得到云南省15 a的年均NDVI數(shù)據(jù)(1985—2010年)，并通過重采樣和裁剪之后得到研究區(qū)的15 a平均NDVI數(shù)據(jù)。

3)空間分辨率為100 m的數(shù)字高程模型(DEM)數(shù)據(jù)。數(shù)據(jù)類型為GRID格式，利用地形圖數(shù)字化得到等高線，借助ArcGIS中的模塊生成不規(guī)則三角網(wǎng)(TIN)，裁剪后最終得到研究區(qū)的DEM數(shù)據(jù)。

4)云南1∶25萬的基礎(chǔ)地理信息數(shù)據(jù)，包含縣級以上行政區(qū)劃及村莊、公路以及河流數(shù)據(jù)。由項目組多年數(shù)據(jù)積累，行政區(qū)劃以及影響人工解譯得到。

1.2 軟件平臺

本研究使用的DOMAIN模型在ModEco v3.02軟件中實現(xiàn)[14]。該軟件由美國加州大學(xué)的Guo Qinghua教授與北京大學(xué)的劉瑜教授共同開發(fā)完成，集成了本研究所使用的多種模型，是研究物種分布與模擬的組合軟件。該模型通過計算環(huán)境空間中待確定類別點到已知類別點的最大相似度來判定其所屬的類別[15]。此外，還使用了ArcGIS 10.1作為環(huán)境梯度因子的計算和制圖軟件。

1.3 稀樹灌木草叢適生條件以及關(guān)鍵因子分析

1.3.1 環(huán)境梯度因子 本研究采用的環(huán)境梯度因子主要是數(shù)字高程模型(DEM)和溫度數(shù)據(jù)。元江縣平均海拔在1800～2000 m，四周山地與元江河谷相對形成低洼環(huán)境，稀樹灌木草叢一般生長在海拔800～1400 m的河谷兩側(cè)，同時處于東南暖濕氣流與西南暖濕氣流的背風(fēng)雨影區(qū)，對形成干熱河谷氣候創(chuàng)造了條件，形成該區(qū)域復(fù)雜的植被類型。因為兩側(cè)高大山體的阻擋，氣流越過高山后，從山頂下降過程中釋放大量的熱量，造成山體從山頂?shù)缴降诇囟仍絹碓礁?，濕度越來越小。氣溫隨著海拔的變化每下降1000 m，溫度上升約6.5 ℃，也就是說，當(dāng)從2500 m的哀牢山山系頂部下降到河谷底部時，溫度變化約13～15 ℃左右，這使得河谷溫度高于周圍其他區(qū)域的溫度，氣候干燥，元江河谷地帶年均溫在20 ℃以上，最冷月均溫16.8 ℃，最熱月均溫28.6 ℃，≥10 ℃的年積溫達8700 ℃，所以本研究在進行空間模擬時使用的DEM和對其生境影響程度很大的溫度作為參數(shù)變量。

1.3.2 NDVI 歸一化植被指數(shù)(NDVI)是反映土地覆蓋植被狀況的一種遙感指標，而且具有便于大尺度研究，可視化效果好等優(yōu)點，在數(shù)據(jù)處理過程中有較強的抗干擾能力，是植被變化監(jiān)測的重要方法[16]。因為干熱河谷區(qū)域受氣候影響，植被覆蓋度較低，尤其是在干季，主要以扭黃茅、黃花稔、疏序黃荊、霸王鞭等旱生植被為主。散生有厚皮樹、蝦子花等低矮的喬木和灌木。大部分區(qū)域0.3

圖1 樣本點分布

1.4 樣本點分布

圖2 模擬結(jié)果分布

本研究共使用295個樣本點，樣本采集時間為2013年6月10—17日，采集地點主要集中于元江國家級自然保護區(qū)、農(nóng)場七大隊、深溝、扒垤龍、施底、二臺坡、坡垤、昆土、熱水塘、觀音灘、西拉河和干壩五隊(圖1)。所采用的儀器工具主要有手持GPS，通過樣線與樣地相結(jié)合的方法記錄樣點的經(jīng)緯度與海拔高度。由于干熱河谷區(qū)域沿元江及其支流的兩岸分布，樣點的采集主要沿河公路以及縣級、鄉(xiāng)級公路采集，因為手持GPS定位儀的經(jīng)緯度信息與實際的經(jīng)緯度信息有差異，以及樣點可達性和分布區(qū)域的限制，本研究通過元江縣的二調(diào)數(shù)據(jù)進行了經(jīng)緯度的訂正，得到樣本點集合，樣點數(shù)據(jù)有代表性。同時模型還采用了環(huán)境梯度與植被指數(shù)作為關(guān)鍵因子，因此樣點分布對預(yù)測結(jié)果影響較小。

2 結(jié)果與分析

2.1 預(yù)測結(jié)果

通過DOMAIN模型利用295個樣本點與DEM、NDVI以及溫度基線數(shù)據(jù)進行計算，模擬得到元江縣稀樹灌木草叢的分布狀況(圖2)。干熱區(qū)域稀樹灌木草叢面積柵格數(shù)量為80062個，面積約為800.62 km2，占總面積的29.448%。在進行數(shù)據(jù)處理過程中，根據(jù)軟件的數(shù)據(jù)統(tǒng)計結(jié)果，統(tǒng)計出295個樣本點所對應(yīng)的坐標、海拔、溫度和NDVI值，便于直觀的從數(shù)值上體現(xiàn)它們之間的相關(guān)性，部分數(shù)據(jù)見表1。

模型檢驗真陽性率(TPR)，也稱為靈敏度檢驗，值越高模型精度越好[17]，DOMAIN模型的TPR=0.9833，同時該模型還通過受試者工作特征曲線(ROC)的檢驗，該方法計算不同閾值正確模擬存在的百分率和正確模擬不存在的百分率，然后將其分別在Y軸和X軸上表示，通過比較曲線和45°線(表示物種處于隨機分布狀態(tài))之間的面積(Area Under the ROC Curve，AUC)確定模型的精度，該方法的評價標準為：0.50

由圖2可知，稀樹灌木草叢主要分布于元江河谷的兩側(cè)，元江縣中部有一個類似于“][”的凹形區(qū)域，這一部分凹形區(qū)域是元江縣城。主要呈現(xiàn)北多南少的分布狀態(tài)，北部分布相對集中，但在南部地區(qū)分布相對分散。在南部相對分散的原因是因為受西南暖濕氣流的影響，在迎風(fēng)坡水熱條件充足，大部分分布著中山濕性闊葉林、半濕潤常綠闊葉林、季風(fēng)常綠闊葉林以及少部分的山地季雨林，是元江國家級自然保護區(qū)的章巴望鄉(xiāng)臺片區(qū)。而在背風(fēng)坡則有少面積的干熱區(qū)域存在，面積不大，但是該區(qū)域人類活動頻繁，植被破壞嚴重，植被覆蓋度低，區(qū)域破碎化程度較高，雖然在進行采樣的時候并未有樣點分布，但是根據(jù)模型預(yù)測結(jié)果以及結(jié)合遙感影像來看，該區(qū)域確實有變?yōu)楦蔁釁^(qū)域的潛在變化，尤其是在與墨江縣交界的區(qū)域。同時模擬結(jié)果顯示在元江縣城南部的山區(qū)有一條帶狀分布的區(qū)域，這一區(qū)域土地破碎，植被破壞后難以恢復(fù)，逐漸形成次生的稀樹灌木草叢群落。

在元江縣北部，干熱區(qū)域相對集中，在元江流經(jīng)該縣的下游部分面積較大，通過對比遙感影像，與模擬的結(jié)果基本相同。從遙感影像上看，北部地區(qū)還有一片覆蓋度較高的闊葉林，不過該區(qū)域南部山地也出現(xiàn)小面積的干熱區(qū)域，對比模擬結(jié)果該區(qū)域南部也有一點零星的分布，但不集中，也從另一面說明模擬結(jié)果的可靠性。

表1 樣本點所對應(yīng)各因子的值

2.2 空間變化分析

從模擬結(jié)果分析，造成這種后果的原因除了與當(dāng)?shù)乇旧淼淖匀画h(huán)境有關(guān)外，還與人類活動的影響有密切關(guān)系，主要表現(xiàn)在以下幾個方面：首先，土地利用方式對干熱區(qū)域的影響。在元江河谷附近的低矮山坡上分布著大面積的香蕉林、芒果林以及火龍果園等經(jīng)濟作物，使得本來就極易受侵蝕的地表抗蝕能力更加弱化。同時農(nóng)地數(shù)量的增大使得在農(nóng)業(yè)作業(yè)過程中，本來土質(zhì)不高的土壤性質(zhì)再發(fā)生變化，這會加劇土壤侵蝕，而土壤侵蝕反過來又進一步改變著土壤的性質(zhì)。其次，是人們的生態(tài)保護意識淡薄，對干熱河谷氣候的成因了解不夠，不計后果的從事農(nóng)業(yè)經(jīng)濟活動，對自己所生活的環(huán)境造成的影響程度。最后，是生態(tài)恢復(fù)工程的滯后性，往往都是在發(fā)現(xiàn)該區(qū)域生態(tài)系統(tǒng)失衡時才會采取措施，如封山育林、加強溝谷地區(qū)的植被覆蓋度等，但是這樣做往往沒有辦法彌補對生態(tài)環(huán)境的破壞，次生植被的保水固土能力大不如前，使得生態(tài)恢復(fù)難上加難。

3 結(jié)論與討論

3.1 結(jié)論

環(huán)境是社會進步的基礎(chǔ)，經(jīng)濟發(fā)展的前提，根據(jù)模擬預(yù)測的結(jié)果對元江縣的干熱區(qū)域進行分析，元江縣干熱區(qū)域稀樹灌木草叢群落的面積呈現(xiàn)出擴大趨勢，尤其是元江北部龍?zhí)多l(xiāng)以北、青龍廠鎮(zhèn)以南的地區(qū)存在更高的演化風(fēng)險，建議將上述區(qū)域及元江國家級自然保護區(qū)章巴片區(qū)列為重點森林保護區(qū)域。

同時模擬結(jié)果也從另一個方面說明該區(qū)域土壤侵蝕、生態(tài)破壞嚴重，森林稀疏，植被覆蓋度進一步下降，對于該區(qū)域植被恢復(fù)的難度加大。以往的生態(tài)恢復(fù)工程的進度遠遠不及經(jīng)濟發(fā)展對環(huán)境的破壞速度，造成的后果極其嚴重，應(yīng)花費更多的財力、人力、物力來進行再修復(fù)，避免造成“破壞—修復(fù)—再破壞—再修復(fù)”的惡性循環(huán)。

3.2 討論

在全球變暖的今天，世界范圍內(nèi)的溫度升高，可能會對本來就形式嚴峻的干熱河谷區(qū)域造成更大影響。受整體的氣候影響，河谷區(qū)域的溫度會升高，意味著降雨減少，蒸發(fā)量增大，那么河谷干熱面積會增大，對于稀樹灌木草叢群落中那些散生的喬木和灌木，它們的生境將更加困難，對于其他植物群落影響也是必然的。

目前，對干熱河谷稀樹灌木草叢的研究，主要從群落生態(tài)學(xué)的角度進行，對其空間分布與其生境的研究比較少。本研究中預(yù)測得到元江縣干熱河谷稀樹灌木草叢分布區(qū)，通過加強對稀樹灌木草叢的保育管理可以達到對干熱河谷區(qū)域生態(tài)保護的效果，模擬得到的稀樹灌木草叢可作為一項指標，結(jié)合地形、氣候及人類干擾等因素來綜合分析其生境變化，從而為該區(qū)域生物多樣性保護和生態(tài)恢復(fù)提供理論依據(jù)和技術(shù)指導(dǎo)。但是也存在一些不足，比如模型的選取上，不同的模型會有不同的結(jié)果，為了使研究結(jié)果更加精確，對于模型的優(yōu)選和其空間的變化趨勢進一步的研究也是今后的方向。

4 建議

針對該區(qū)域特殊的氣候條件與植被類型，不僅要考慮當(dāng)下的經(jīng)濟利益，更要考慮長遠發(fā)展的利益。應(yīng)多角度、全方位、深層次認識到保護干熱河谷生態(tài)環(huán)境的重要性和緊迫性，良好的生態(tài)系統(tǒng)對干熱河谷生態(tài)退化及植被類型災(zāi)變有明顯的減緩作用。因此，本文建議從生態(tài)恢復(fù)的角度防治干熱河谷災(zāi)變。

4.1 遵循可持續(xù)發(fā)展原則，選擇適生物種

在地力低下的河谷山坡上，要將植被恢復(fù)與森林資源的培育、開發(fā)利用相結(jié)合，走生態(tài)林業(yè)道路，選擇適生的喬木物種，借鑒在過去生態(tài)恢復(fù)工程的成功經(jīng)驗，營建高價值的印度黃檀和紫膠寄主林，培育優(yōu)良的寄主久樹(Schleicheraoleosa)，用以生產(chǎn)優(yōu)質(zhì)的紫膠；同時栽植新銀合歡與桉樹混交林，以期對坡地起到水土保持的作用，以及農(nóng)用小徑材、薪炭材水土保持林等，以達到水土保持林與防護林相結(jié)合的多目標人工生態(tài)系統(tǒng)。

4.2 建立生態(tài)恢復(fù)評價體系，進行河谷坡地改良技術(shù)的研究

對于河谷坡地，水土保持是生態(tài)恢復(fù)的重要評價指標，而選擇耐旱力強、速生、多功能的多年生植物是關(guān)鍵。因此，干熱河谷兩側(cè)的坡地資源的改良技術(shù)尤為重要，應(yīng)按照其所受自然限制因素及其強度加以區(qū)別劃分。主要有溫度、土壤成分、土層厚度、裸巖所占面積、地表保水程度等作為評價指標，從綜合自然因素原則為指標定量化分析，在此基礎(chǔ)上進行生態(tài)恢復(fù)評價體系。但是本研究沒有進行專門的指標體系研究，將作為后續(xù)的研究工作。

4.3 改變土地利用方式，提高土地利用率

在經(jīng)濟發(fā)展的同時要兼顧生態(tài)環(huán)境的平衡發(fā)展，多開展普法宣傳工作，從人們的思想意識上改變其固有的思維方式，使其意識到環(huán)境保護的重要性，同時政府要做好土地利用規(guī)劃工作，提高土地利用率，發(fā)展三面溝+護埂林草+逐年往下翻土的農(nóng)林牧逐年梯帶化技術(shù)。對于地力低下的區(qū)域，要加強植被恢復(fù)，政府要確?；謴?fù)工作先行者的地位，帶動群眾一起來保護河谷流域的植被。

[1]趙琳，郎南軍，鄭科，等.云南干熱河谷生態(tài)環(huán)境特性研究[J].林業(yè)調(diào)查規(guī)劃，2006，31(3)：114-117.

[2]盧培澤，朱彤，何永彬，等.云南干熱河谷坡地資源改良利用途徑與技術(shù)的研究[J].熱帶地理，2001，21(2)：118-120.

[3]黃成敏，何毓蓉.云南省元謀干熱河谷土壤水分的動態(tài)變化[J].山地研究，1997，15(4)：234-238.

[4]李昆，陳玉德.元謀干熱河谷人工林地的水分輸入與土壤水分研究[J].林業(yè)科學(xué)研究，1995，8(6)：651-657.

[5]張信寶，陳玉德.云南元謀干熱河谷區(qū)不同巖土類型荒山植被恢復(fù)研究[J].應(yīng)用與環(huán)境生物學(xué)報，1997，3(1)：13-18.

[6]聶小軍，張建輝，劉剛才，等.金沙江干熱河谷侵蝕陡坡植被恢復(fù)對土壤質(zhì)量的影響[J].生態(tài)環(huán)境，2008，17(4)：1636-1640.

[7]田園，王靜，胡燕.我國干熱河谷地區(qū)土壤侵蝕研究進展[J].中國水土保持，2013(6)：51-54.

[8]金振洲.滇川干熱河谷種子植物區(qū)系研究[J].廣西植物，1999，19(1)：1-14.

[9]南嶺，郭芬芬，王小丹，等.云南元謀干熱河谷區(qū)典型植被恢復(fù)模式的水土保持效應(yīng)[J].安徽農(nóng)業(yè)科學(xué)，2011，39(9)：5168-5171.

[10]何永彬，盧培澤，朱彤.橫斷山——云南高原干熱河谷形成原因[J].資源科學(xué)，2000，22(5)：69-72.

[11]盧培澤，朱彤，何永彬.云南紅河流域的水土流失與重點防治[J].云南環(huán)境科學(xué)，2000，19(1)：21-23.

[12]高文學(xué)，王志和.金沙江干熱河谷稀樹灌草叢植被恢復(fù)方式研究[J].林業(yè)調(diào)查規(guī)劃，2005，30(3)：78-91.

[13]柴宗新，范建容.金沙江干熱河谷植被恢復(fù)的思考[J].山地學(xué)報，2001，19(4)：381-384.

[14]馮永康.湖北省神農(nóng)架溫性及寒溫性針葉林空間分布的多模型預(yù)測及比較[D].南京：南京大學(xué)，2010.

[15]Carpenter，Q.Gillison，A.N..Winter，Domain-a flexible modeling procedure for mappingpotential distributions of plants and animals[J].Biodiversity and Conservation，1993(2)：667-680.

[16]田慶久，閔祥軍.植被指數(shù)研究進展[J].地球科學(xué)進展，1998，13(4)：327-333.

[17]Engler，R.Guisan，A.，Rechsteiner，L.An improved approach for predicting the distribution of rare and endangered species from occurrence and pseudo-absence data[J].AppL EcoL，2004(41)：263-274.

[18]Swets，K.A.Measuring the accuracy ofdiagnostic systems[J].Science，1988(240)：1285-1293.

Spatial Distribution of Savanna Shrub & Grass in Dry-hot Valley on the Base of the DOMAIN Model

WANG Yi-kai，ZHANG Meng-meng，DUAN Wei-hu，YANG Kai-yue，ZHOU Ru-liang，YE Jiang-xia

(SouthwestForestryUniversity，Kunming650224，Yunnan，China)

Taking temperature，DEM and NDVI as the parameters，using the DOMAIN model，we simulate the Spatial Distribution of Savanna shrub& grass in Dry-hot Valley of YuanJiang and realize visualization of the spatial data.The simulation results show：Savanna shrub& grass is mainly distributed in the both sides of the valley，which is concentrated in the northern and dispersed in the northern.Dry-hot areas present the tendency of expansion under the general influence of natural and artificial factors，espesially for those areas with higher evolution risks，such as the north of LongTan township，the south of QingLongchangtownship.We suggest that deserve high priority in these regions and ZhangBa District of YuanJiang National Nature Reserve.

dry-hot valleys (DHV)；savanna shrub& grass；DOMAIN model；spatial distribution；Yuanjiang county

2015-01-05；

2015-02-11

云南省科技創(chuàng)新人才計劃項目(2014HC01)；國家林業(yè)局林業(yè)公益性行業(yè)科研專項項目(天地觀測數(shù)據(jù)同化的森林火險預(yù)報新技術(shù)， 201404402-2)

王一凱(1986—)，女，山西忻州人，西南林業(yè)大學(xué)碩士研究生，從事3S技術(shù)應(yīng)用與研究。E-mail：wangyikail@126.com。

葉江霞，四川成都人，西南林業(yè)大學(xué)副教授，從事3S技術(shù)在生態(tài)學(xué)、森林保護學(xué)中的應(yīng)用。E-mail：yjx125@163.com。

10.13428/j.cnki.fjlk.2015.04.025

S717.1

A

1002-7351(2015)04-0116-05