山地縣域聚落生態位分布特征及其與地形因子的關系

劉延國

摘要為探究少數民族山區不同地形因子對聚落生態位分布空間格局的影響，以岷江上游少數民族山區為例，基于DEM和遙感影像及GIS技術，研究了山地縣域聚落生態位在不同地形因子上的分布特征，結合6個景觀指數的計算，定性和定量分析了聚落生態位格局與地形因子的相關性.結果表明：聚落生態位總面積較小，僅占區域面積的3.98%，單個面積也較小，平均僅為0.61 km2，分布存在明顯的垂直梯度特征，且集中分布于15～35°之間的緩斜坡及斜坡，45～135°與225～315°的半陽坡以及135～225°的陽坡面上；景觀格局指數也顯示出聚落生態位分布隨地形變化的梯度特征.由此可見，高程、坡度和坡向這三類地形因子對山地聚落生態位的大小、形狀和空間分布起著重要的影響作用.

關鍵詞生態位；山地聚落；地形；景觀格局；岷江上游

中圖分類號K90； P94文獻標識碼A文章編號10002537（2017）01001709

山區聚落作為山區群眾生存與生活的基本空間組織單元，是人類與周圍地質、地貌、大氣、水體、土壤和植被等環境要素之間交流的平臺，其生態位反映聚落在山地垂直方向上所處的空間位置及人類生活所能利用的資源空間[13]，也是長期以來山區群眾利用自然資源和適應自然環境的智慧結晶與外在表現形式，蘊含了豐富的生態和地理信息.山區（地）作為以地形為主導的復雜生態系統，其土地利用景觀垂直變化梯度大約是水平變異梯度的1000倍[46]，是影響局部地區土地利用變化最直接和最關鍵的因素，在一定程度上決定著土地的利用類型[714].地形因子作為山區土地利用格局分布的重要影響因素，其變化直接影響到地面的物質流動與能量轉化，對人類活動的限制性和對不同地類空間分布的自然選擇性，使不同地類及聚落在地形梯度上的分布呈現出規律變化特征[6，10].數字高程模型（DEM）作為研究空間變化的基礎數據，蘊含著大量的地形結構和特征信息[4，5，15].基于DEM和遙感影像及GIS技術，通過提取高質量的高程、坡度和坡向等地形因子信息，定量研究山地縣域聚落生態位在不同地形因子上的分布特征，可以進一步厘清山區聚落生態位的分布格局特征，為山區聚落優化布局提供依據.同時，鑒于山區多是貧困集中區的特點[16，17]，而聚落作為貧困的直接載體，其研究結果對理解區域貧困及開展減貧策略的制定也有重要意義.

1區域概況及數據來源

1.1區域概況

岷江上游地處四川盆地丘陵山地向川西北高原的過渡地帶，地質構造復雜，斷裂發育，新構造運動隆升強烈，是青藏高原東麓自然地理環境的典型代表.地理位置介于31°26′～33°16′N和102°59′～104°14′E之間，地勢由西北向東南降低，地表切割由北向南加劇，屬于中深切割地形，平均高差在1 000 m以上，中、高、低山所占比例相差懸殊全國罕見，且該區域具有水平地帶性和垂直地帶性的雙重氣候特征，并在局地焚風效應的影響下，干旱河谷

的隱域性分異現象明顯.岷江上游行政區劃隸屬于四川省阿壩藏族羌族自治州，其范圍大致與汶川、理縣、茂縣、黑水及松潘五縣的行政界線相一致，面積2.24萬km2，人口39.2萬人，農業人口占75%以上，少數民族占總人口的89%以上，是四川省內藏族、羌族和回族聚落的集中分布區，我國主要的羌族聚居區.在藏—羌—漢—回多元民族交融的發展進程中，形成了我國西南地區一條重要的民族廊道以及自然景觀帶.岷江上游地區在多源文化背景下，長期以來的人類生計活動對山地生態環境的持續干擾，加上活躍地震帶的疊加影響，使其成為長江上游典型的生態環境脆弱區，川西農牧民貧困連片集中區（4個省級貧困縣，1個國家級貧困縣黑水縣），該區一直是扶貧開發關注的重點區域，也是學術研究的熱點區[1820].

1.2數據來源

本研究采用的遙感數據包括3個方面：（1）SPOT5（2009～2010年，2.5 m全色和10 m多光譜波段）影像數據，成像傳感器均HRG2，數據格式均為GEOTIFF，產品級別均為1 A；（2）數字高程模型（DEM）；（3）Google Earth三維立體影像及在線免費數據.采用田野調查、3S技術相結合的方法，開展對岷江上游的聚落生態位的定量信息提取.

2數據處理與研究方法

2.1數據處理

2.1.1遙感影像預處理采用ERDAS IMAGINE 9.2對SPOT遙感影像進行幾何校正及分辨率融合.以實地采取的控制點為地理參考，通過Google Earth三維立體影像對多光譜和全色波段數據進行校正，實現多光譜和全色波段數據的融合.

2.1.2聚落生態位信息提取聚落生態位是以聚落為中心，占有一定地理空間的實體范圍，包含房屋建筑、農田、經濟林地、部分斜坡草地以及道路等要素的人工—自然綜合體，故以融合影像為數據源，輔助區域Google Earth影像，在ArcGIS 10.2中采用目視解譯法提取岷江上游山區聚落生態位空間信息.

2.1.3基于DEM的高程、坡度和坡向分級利用DEM 數據，采用ArcGIS中的Spatial Analyst模塊下的表面分析分別提取高程、坡度和坡向信息.高程是影響土地利用方式的重要環境因子，隨著高程的增加，大氣濕度、溫度都會發生明顯的變化，導致土地的利用方式隨高程的變化發生一定的變化.坡度對土地利用有著重要的影響：小于6°的坡地，可以作為農用旱地或牧草地；大于6°的坡地，易產生強烈的侵蝕，需修筑梯田或采用水保耕作法等水保措施；25°是退耕還林還牧界限；溝坡地的地面坡度大部分在35°以上，該類土地應以種草造林為主要利用方式，以保護邊坡的穩定性，防止崩塌、錯落等重力侵蝕發生；45°是植樹造林的上限[15].坡向表征了地面任何一點高程值改變量的最大變化方向.為了便于對聚落生態位在不同地形因子上進行特征分析，結合相關研究成果及研究區實際[2，3，15，21]，對聚落生態位高程、坡度、坡向分級如表1.

2.2研究方法

景觀指數是指能夠高度概括景觀格局信息，反映結構組成和空間配置等方面特征的簡單定量指標[2224].根據研究需要，分別從斑塊水平、類型水平和景觀水平3個層次，選取斑塊數量、斑塊密度、斑塊平均面積、面積加權形狀指數、面積加權分維數、斑塊分離度指數，研究聚落生態位景觀格局狀況.以提取到的聚落生態位信息為數據源，在ArcGIS 10.2中將提取到的聚落生態位矢量數據柵格化，用于在Fragstats 4.2中進行景觀指數的計算，輔助SPSS 13.0進行統計分析.

式中，n為聚落生態位數量，Pi為聚落生態位周長，ai是聚落生態位面積，A為縣域聚落生態位總面積.

面積加權分維數（AWMPFD）： AWMPFD運用分維理論來測量斑塊和景觀的空間形狀復雜性，是反映景觀格局總體特征的重要指標.AWMPFD=1代表斑塊形狀最為簡單的正方形或圓形，AWMPFD=2代表周長最復雜的斑塊類型，

采用下列公式[2224]：

3結果與分析

通過ArcGIS 10.2對岷江上游山區聚落生態位的提取結果（表2）顯示，聚落生態位總計1 667個，區域總體聚落生態位占區域總面積比例都較小，最小的理縣僅占209%，一定程度上反映山區適宜的聚落生態位面積區域較小；汶川平均海拔最低，松潘處在高原面上，平均海拔最高，區域垂直分異特征明顯，隨著區域平均海拔的升高，聚落生態位密度變小，平均面積呈增大趨勢，汶川最小，松潘平均最大超過113公頃，說明高海拔區域聚落需要更大的資源空間以維持其生計.

3.1聚落生態位在不同地形上的分布特征

3.1.1在不同海拔梯度上的分布特征將DEM和聚落生態位疊加，統計聚落生態位在不同海拔梯度上的分布，結果如表3所示.

隨著海拔的升高，聚落生態位面積逐步增大；汶川縣聚落生態位整體分布海拔最低，松潘縣聚落生態位整體分布海拔最高，且海拔3 200 m以上聚落生態位面積占25%以上，屬典型的高山聚落；汶川、理縣、茂縣及黑水縣，聚落大都分布在1 300～2 800 m之間，處在河谷暖溫帶半干旱氣候（2 000 m以下）、溫帶半干旱河谷氣候（2 000～2 500 m）以及山地寒溫帶氣候（2 500～3 000 m）帶上；松潘縣聚落主要分布在2 800～3 200 m之間，處于山地亞寒帶氣候（3 000～4 000 m）帶上.

3.1.2聚落生態位在不同坡度上的分布特征將坡度和聚落生態位疊加，統計聚落生態位在坡度上的分布，結果如表4.

（1）岷江上游聚落生態位集中分布于15°～35°之間的緩斜坡及斜坡山地上，分布數量最多，面積占比最大；聚落生態位平均面積隨坡度的升高逐漸降低，說明坡度高的山地很難找到合適的生態位，平均周長隨著坡度升高也降低，形狀趨于簡單；松潘縣處在高原臺面上，15°以下緩地分布聚落面積占40%以上，汶川、理縣、茂縣及黑水15°以下緩地分布都在10%左右.

（2）<6°的平緩地上，松潘縣聚落生態位面積最大，平均面積最大，高于茂縣，最低的為黑水縣在2%以下；6°～15°的緩地上，松潘縣聚落生態位占據明顯優勢，面積占比接近25%，其次為汶川，最低的為茂縣、黑水都在9%左右；15°～25°的緩斜坡上，汶川、茂縣及黑水的面積超過30%，其中茂縣數量占比超過30%，最少的理縣數量超過20%；25°～35°的斜坡上，汶川、理縣、茂縣及黑水的聚落生態位分布面積最大，其中黑水面積占比超過40%，數量超過37%，除松潘以外數量比例都接近30%；>35°以上的斜坡上，平均面積都較小，其中汶川聚落生態位分布數量最多，其次是黑水，松潘最少，理縣面積占比及數量占比都最高，分別為2719%及34.51%.

3.1.3聚落生態位在不向坡向上的分布特征將坡向和聚落生態位疊加，統計聚落生態位在坡向上的分布，結果如表5.

（1）岷江上游聚落生態位主要分布于45°～135°和225°～315°以及135°～225°之間半陽坡及陽坡山地上，半陽坡山地分布最為集中，且聚落生態位平均面積較小；陰坡聚落生態位也占據一定比例，汶川、理縣、茂縣總面積占比超過20%，最高的理縣達32%，平均面積也較大.

（2）0°～45°和315°～360°的陰坡上，汶川縣數量分布最多達116個，面積占比最高的為理縣，松潘縣數量及面積占比都最少；45°～135°和225°～315°的半陽坡上，聚落分布最多，最高的松潘超過58%，最少的理縣也超過44%，面積占比也都在40%以上；135°～225°的陽坡上，黑水與松潘占明顯優勢，最小的汶川面積占比也超過22%.

3.2不同地形因子聚落生態位景觀格局指數特征

3.2.1總體景觀格局特征岷江上游聚落生態位總體景觀格局計算結果如表6，岷江上游山區聚落生態位占國土面積比例都較小，平均僅為3.98%，單個聚落生態位的面積也較小，最大的僅為1.14 km2；隨著區域整體平均海拔的升高，聚落生態位密度變小，平均面積變大，高寒區域受作物生長季節及產出的限制[18，20，25]，要滿足生產生活需要的作物產出，目前只有通過擴大作物種植面積來取得，因此聚落生態位面積低海拔區域大.總體而言整個岷江上游地區聚落聚集度明顯，隨著區域平均海拔的升高，聚落區域集中分布，且形狀趨于不規則，說明高海拔區域聚落生態位處于自然狀態.

3.2.2不同海拔梯度上聚落生態位集中分布區景觀格局特征不同海拔梯度岷江上游聚落生態位集中分布景觀格局特征如表7.

（1）隨著海拔梯度的升高，聚落生態位密度逐漸降低，分離度指數逐漸變小，在一定程度上說明高海拔聚落集中趨勢明顯，符合便于勞動互助與交往等方面的社會性需要特征，聚落面積逐漸變大，且形狀趨于不規則.

（2）1 300～2 200 m海拔區間，理縣聚落生態位比汶川更聚集，且面積更小，形狀更復雜；2 200～2 800 m海拔區間，高原面上的松潘聚落生態更集中，面積更大，形狀也比汶川、理縣、茂縣更為復雜；2 800～3 200 m海拔區間，松潘縣聚落生態位更復雜也更集中，其中黑水縣在此區段生態位密度最大.

（3）低海拔區域的汶川而言，隨著海拔梯度升高，聚落生態位形狀趨于復雜，而高原面上的松潘聚落形狀指數逐漸減小，形狀趨于簡單.

3.2.3不同坡度聚落生態位集中分布區景觀格局特征不同坡度岷江上游聚落生態位集中分布景觀格局特征如表8所示.

（1）隨著坡度的升高，聚落生態位密度逐漸升高，集中趨勢明顯，且形狀趨于不規則.

（2）6°～15°區間段，低海拔區域的汶川比高原面上的松潘更集聚，密度更大，且形狀更規則，說明汶川聚落生態位在緩地上的人為干預程度高；15°～25°的緩斜坡上，茂縣聚落生態位形狀最復雜，集聚程度而言，汶川最為集中，松潘較松散，密度而言，高原面上的松潘依然最小；25°～35°的斜坡上，汶川、茂縣及黑水縣聚落生態位數量分布較多，且比理縣及松潘縣要松散；茂縣聚落生態位形狀最復雜；>35°的斜坡上，依然分布有數量較多聚落，除松潘高原面上，其余4縣域聚落生態位密度大，面積小，形狀比低坡度聚落生態位復雜.

3.2.4不同坡向聚落生態位集中分布區景觀格局特征不同坡向上岷江上游聚落生態位集中分布景觀格局特征如表9所示.

（1）聚落生態位在半陽坡分布最多，密度最大，聚落形狀從陰坡到半陽坡及陽坡逐漸趨于簡單，從平均低海拔的汶川到高原面上的松潘，集中趨勢明顯.

（2）0°～45°及315°～360°的陰坡上，汶川分布聚落最多，密度最大，松潘最少，且形狀最復雜，茂縣聚落生態位形狀也比較復雜；45°～135°及 225°～315°的半陽坡上，聚落生態位分布最多，密度最大，形狀較陰坡簡單；135°～225°的陽坡上，集聚度明顯，形狀較半陽坡簡單，從另外一個角度說明山區陽坡適宜生態位面積較小，半陽坡為山區聚落的適宜生態位.

3.3結果分析

從聚落生態位在不同地形上的分布特征來看，具有明顯的海拔梯度特征，隨著海拔梯度升高，聚落數量減小，面積逐漸增大，不論在哪一個海拔區間，聚落生態位多選擇在坡度較平緩的半陽坡及陽坡區域，這與岷江上游的山地特征及特殊的氣候特征相吻合，也是少數民族群眾在首先滿足基本的生存需要前提下，提高生產效率與及生活便捷水平的必然選擇；從聚落生態位景觀格局指數特征來看，基本與不同地形上的分布特征相一致，聚落生態位密度逐漸降低，分離度指數逐漸變小，符合便于勞動互助與交往等方面的社會性需要特征，聚落生態位形狀趨于不規則，主要是高海拔區域，受作物生長季節及產出的限制[18，20，25]，要足夠的滿足生產生活需要的作物產出，目前只有通過擴大作物種植面積來取得，因受地域空間限制程度較小，生態位更多處于自然狀態，形狀更不規則.

綜上，岷江上游山區聚落生態位分布及其空間格局特征是長期以來少數民族群眾利用自然資源和適應自然環境的產物，是其智慧的結晶與外在表現形式，但由于不合理的人類生計活動及活躍地震帶的影響，岷江上游也成為長江上游典型的生態環境脆弱區和川西農牧民貧困連片集中區，到2020年實現區域全部脫貧，是四川省及國家的現實要求與戰略目標，本文的研究結論也能從一個側面為區域精準扶貧策略的制定提供一定參考.

參考文獻：

[1]MACKENZIE A， BALL A S， VIRDEE S R. Ecology[M]. Beijing： Science Press， 2007.

[2]馬旭，王青，丁明濤，等.岷江上游山區聚落生態位及其模型[J].生態與農村環境學報， 2012，28（5）：574578.

[3]王青，石敏球，郭亞琳，等.岷江上游山區聚落生態位垂直分異研究[J].地理學報， 2013，68（11）：15601567.

[4]LAURE W. Human development and environment in the Andes： A geoecological overview[J]. Mount Res Devel， 1993，13（2）：157166.

[5]KRNER C. The use of “altitude” in ecological research[J]. Trends Ecol Evol， 2007，22（11）：569574.

[6]哈凱，丁慶龍，門明新，等.山地丘陵區土地利用分布及其與地形因子關系—以河北省懷來縣為例[J].地理研究， 2015，34（5）：909921.

[7]喬青，高吉喜，王維.川滇農牧交錯區地形特征對土地利用空間格局的影響[J].長江流域資源與環境， 2009，18（9）：812818.

[8]孫然好，陳利頂，張百平，等.山地景觀垂直分異研究進展[J].應用生態學報， 2009，20（7）：16171624.

[9]黃海洋，楊慶媛，王成.基于DEM的土地利用類型與地形因子關系研究：以重慶市酉陽縣麻旺鎮為例[J].西南大學學報（自然科學版）， 2009，31（4）：159164.

[10]梁發超，劉黎明.基于地形梯度的土地利用類型分布特征分析：以湖南省瀏陽市為例[J].資源科學， 2010， 32（11）：21382144.

[11]張錦明，游雄.地形起伏度最佳分析區域研究[J].測繪科學技術學報， 2011，28（5）：369373.

[12]郭洪峰， 許月卿， 吳艷芳. 基于地形梯度的土地利用格局與時空變化分析： 以北京市平谷區為例[J]. 經濟地理， 2013，33（1）：160166.

[13]李丹，劉丹丹，趙金祥.基于DEM的山區土地利用變化分析[J].水土保持研究， 2014，21（1）：6670.

[14]國巧真，寧曉平，王志恒，等.地形地貌對半山區土地利用動態變化影響分析——以天津市薊縣為例[J].國土資源遙感， 2015，27（1）：153159.

[15]趙玲，吳良林，莫建飛.基于DEM的桂西北土地利用與地形關系特征分析[J].地理空間信息， 2011，9（1）：103105.

[16]陳全功，程蹊.少數民族山區長期貧困與發展型減貧政策研究[M].北京：科學出版社， 2014.

[17]國家統計局住戶調查辦公室.中國農村貧困監測報告2011[M].北京：中國統計出版社， 2011.

[18]陳國階，涂建軍，樊宏，等. 岷江上游生態建設的理論與實踐[M]. 重慶：西南師范大學出版社， 2006.

[19]劉邵權. 農村聚落生態研究： 理論與實踐[M]. 北京： 中國環境科學出版社， 2006.

[20]印開蒲.百年追尋：見證中國西部環境變遷[M]. 北京：中國大百科全書出版社， 2010.

[21]范建容，張子瑜，李立華.四川省山地類型界定與山區類型劃分[J].地理研究， 2015，34（1）：6573.

[22]王宗明， 陳銘， 宋開山. 三江平原別拉洪河流域濕地農田化過程中濕地——農田景觀梯度時空特征分析[J]. 水土保持學報， 2008，22（1）：194198.

[23]齊偉，曲衍波，劉洪義，等.區域代表性景觀格局指數篩選與土地利用分區[J].中國土地科學， 2009，23（1）：3337.

[24]鄭新奇，付梅臣.景觀格局空間分析技術及其應用[M].北京：科學出版社， 2010.

[25]阿壩州人民政府. 阿壩州年鑒[M]. 成都： 巴蜀書社， 2014.