高鳳杰 ， 馬泉來 ， 單培明 ， 劉 洋 ， 郭欣欣 ， 楊 雙 ， 張志民

(東北農(nóng)業(yè)大學(xué) 資源與環(huán)境學(xué)院，哈爾濱150030)

黑龍江省穆棱市土地利用/覆被變化及熱點分析

高鳳杰 ， 馬泉來 ， 單培明 ， 劉 洋 ， 郭欣欣 ， 楊 雙 ， 張志民

(東北農(nóng)業(yè)大學(xué) 資源與環(huán)境學(xué)院，哈爾濱150030)

土地利用/覆被變化及其過程特征是全球變化領(lǐng)域研究的重點和熱點。運用地統(tǒng)計學(xué)和3S技術(shù)，對黑龍江省穆棱市1991—2001年和2001—2013年土地利用/覆被變化及其空間演變特征進(jìn)行研究。結(jié)果表明：1991—2001年間，耕地是主要轉(zhuǎn)入地類，增加232.25 km2，林地作為最大轉(zhuǎn)出地類，減少290.35 km2，林地與耕地的轉(zhuǎn)化在土地利用變化中占最大比例，其他用地均發(fā)生變化，但面積小且散亂，整體上該時期土地利用呈現(xiàn)出開發(fā)混亂、效率低下的特點；2001—2013年間，受還林還草的影響，耕地減少179.01 km2，林地增加171.45 km2，其中坡耕地的退出占較大比例，且分布面積廣，該時段內(nèi)土地利用結(jié)構(gòu)得到優(yōu)化；土地利用變化的熱點區(qū)由2001年前的“兩大一小”模式，演變?yōu)橐允袇^(qū)和中心鄉(xiāng)鎮(zhèn)圍成的重點區(qū)域；1991—2013年間研究區(qū)自然因子變化較小，土地利用/覆被變化主要受退耕還林等國家政策的主導(dǎo)。

土地利用/覆被變化；地統(tǒng)計學(xué)；3S技術(shù)；黑龍江省穆棱市

0 引言

土地利用/覆被變化(LUCC)作為人類活動和自然過程相互作用最為密切的問題，一直是全球變化領(lǐng)域研究的重點和熱點，引起國內(nèi)外學(xué)者的廣泛關(guān)注[1-3]。“國際地圈生物圈計劃(IGBP)”和“全球變化人文因素計劃(IHDP)”共同制定了“土地利用/覆被變化(LUCC)科學(xué)研究計劃”和“全球土地計劃(GLP)”，并作為全球變化及人地關(guān)系研究領(lǐng)域的前沿問題。從現(xiàn)有成果看，有影響力的研究主要從全球和其他較大尺度上展開[4-5]，涉及土地利用變化的驅(qū)動力因子、區(qū)域土地利用的現(xiàn)狀分析及土地利用/覆被變化引起的景觀和生態(tài)環(huán)境效應(yīng)等響應(yīng)機(jī)制[6-8]，并且有在全球研究的框架之下逐步向中小尺度過渡的趨勢，其中，以典型流域和西北黃土丘陵等生態(tài)脆弱區(qū)為研究單元較多[9-10]，且多以定性分析和指標(biāo)計算方法為主，農(nóng)林牧交錯區(qū)是多種產(chǎn)業(yè)相互發(fā)展又深具干擾性的交錯地帶，尺度適中，也因其重要及特殊的生態(tài)區(qū)位而引起關(guān)注，并且取得了相關(guān)成果[11-12]。而結(jié)合地統(tǒng)計學(xué)的理論方法，在小單元內(nèi)研究土地利用空間格局演變可以更好地解釋其數(shù)量和空間變化特征，通過對其長時間序列土地利用變化的可視化演變分析，逐步推演其格局的演變過程和成因機(jī)制，并進(jìn)行未來情景演變的模擬預(yù)測，是未來研究的重點和發(fā)展方向。

黑龍江省穆棱市是東北地區(qū)典型的農(nóng)林交錯帶，在21世紀(jì)初開始實施退耕還林工程，土地利用/覆被變化劇烈，農(nóng)林用地呈現(xiàn)新的空間演變特征，選取該區(qū)域為研究范圍，采用遙感影像及相關(guān)統(tǒng)計數(shù)據(jù)，借助于3S技術(shù)，建立1991—2013年間土地利用數(shù)據(jù)庫，以1 km × 1 km柵格為單元，建立土地利用總體變化指數(shù)和土地利用類型面積變化指數(shù)，將不同時間點的土地利用數(shù)據(jù)及相關(guān)指數(shù)進(jìn)行空間展布，定量精確地獲取不同單元內(nèi)人類活動和社會發(fā)展對土地利用變化的影響強(qiáng)度及空間差異性，分析熱點區(qū)域變化特征，同時，結(jié)合研究區(qū)氣象及社會等要素對土地利用變化的驅(qū)動力進(jìn)行探討，為穆棱市區(qū)域發(fā)展的決策與管理提供理論依據(jù)。

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

1.1 研究區(qū)概況

穆棱市地處黑龍江省東南部，位于129°45′19″～130°58′07″E 、43°49′55″～45°07′16″N，南與吉林省接壤，東西部分別有老爺嶺和太平嶺南北走向分布，全境海拔232～1 111 m，平均海拔600 m，多為低山丘陵，屬溫帶大陸性季風(fēng)氣候。全市下轄六鎮(zhèn)兩鄉(xiāng)，市區(qū)在八面通鎮(zhèn)，總面積6 673 km2，總?cè)丝?0萬，生產(chǎn)總值超過100億元，是黑龍江省近些年發(fā)展速度較快的縣市之一。境內(nèi)河流眾多，穆棱河從南向北貫穿全境，林地在土地利用中占比較大，森林覆蓋率達(dá)60%以上，年均降水量530 mm，主要集中于6—8月份，占全年降水量的60%以上，耕地土層較薄，大部分為暗棕壤，是黑龍江省典型的農(nóng)林交錯區(qū)和生態(tài)敏感區(qū)。20世紀(jì)末，穆棱市人口快速增長，生活需求擴(kuò)大，以農(nóng)為綱思想影響廣泛，從而造成毀林開荒、超破耕作、水土流失、植被破壞等現(xiàn)象嚴(yán)重，林地、草地大面積減少，耕地驟增，土地資源遭到不合理開發(fā)和利用，水土流失情況嚴(yán)重，生態(tài)環(huán)境脆弱[13-14]。

1.2 數(shù)據(jù)來源

采用1991，2001年的Landsat5 TM影像和2013年的Landsat8 OLI影像，為了便于影像特征識別與處理，分別選取研究區(qū)9，10月份影像，空間分辨率為30 m，在ENVI 4.5中進(jìn)行波段融合、影像的配準(zhǔn)校正，在ArcGIS 10.0平臺以穆棱市所在的129°E經(jīng)線為中央經(jīng)線進(jìn)行投影轉(zhuǎn)換，建立解譯標(biāo)志，進(jìn)行人工目視解譯得到穆棱市不同時間點的土地利用類型分布圖，之后進(jìn)行研究區(qū)樣點選取和精度驗證，通過專家咨詢和野外調(diào)研選取驗證樣點，在野外通過手持GPS實現(xiàn)研究區(qū)土地利用類型的定點驗證，結(jié)果表明，土地利用類型判別準(zhǔn)確度在90%以上，滿足研究要求。同時，參照國家土地利用劃分標(biāo)準(zhǔn)并結(jié)合研究區(qū)實際狀況和本研究的目的，將研究區(qū)土地類型劃分為旱地、水田、林地、草地、水域和建設(shè)用地6種類型，以穆棱市地形圖、水系圖及行政區(qū)劃圖為輔助圖件數(shù)據(jù)，建立土地利用數(shù)據(jù)庫。相關(guān)社會經(jīng)濟(jì)統(tǒng)計數(shù)據(jù)主要來源于牡丹江地區(qū)統(tǒng)計年鑒和相關(guān)部門網(wǎng)站，氣象數(shù)據(jù)來源于中國氣象科學(xué)數(shù)據(jù)共享平臺。

2 研究方法

在時間尺度上，以退耕還林工程的實施為理論分割依據(jù)，以2001年為時間尺度分割點，劃分前后兩個時間段，即對退耕還林前后兩個時期展開研究。

2.1 土地利用轉(zhuǎn)移矩陣

運用1991—2001年和2001—2013年兩個時間段、3個時間點的穆棱市土地利用數(shù)據(jù)，通過ArcGIS 10.0中的 Intersect功能和在Excel中添加透視表對解譯得到的3期土地利用數(shù)據(jù)進(jìn)行空間疊加分析和統(tǒng)計，獲得穆棱市1991—2001年、2001—2013年前后兩個時間段的土地利用類型轉(zhuǎn)移矩陣。對比在研究區(qū)受到外部強(qiáng)烈干擾情況下，即退耕還林前后的土地利用/覆被變化的數(shù)量和結(jié)構(gòu)特征，并做進(jìn)一步分析。

2.2 土地利用變化指數(shù)

2.2.1 土地利用總體變化指數(shù)。土地利用總體變化指數(shù)反映區(qū)域土地利用變化的程度，本研究以公里網(wǎng)格為獨立單元，精確反映研究區(qū)不同地理單元內(nèi)部土地利用變化部分所占比例及土地利用變化的速率，整體上也可以反映人類活動對土地利用的干擾強(qiáng)度。公式為：

(1)

式中：Iluc為土地利用總體變化指數(shù)；uai,ubi分別為研究時段初年和末年第i類土地利用類型面積；D為柵格單元面積，除研究區(qū)邊緣外，其他均為1 km2。

2.2.2 土地利用類型面積變化指數(shù)。土地利用類型面積變化指數(shù)可以作為進(jìn)一步深入分析土地利用變化的狀況指標(biāo)，具體是指在單元格內(nèi)土地利用總體變化指數(shù)的基礎(chǔ)上進(jìn)一步得出每一單元內(nèi)發(fā)生變化的土地中各土地類型之間的轉(zhuǎn)化情況及其占單元格內(nèi)所有土地利用變化面積的比重。土地利用變化指數(shù)是局部中的整體反映，而土地利用類型面積變化指數(shù)是對前者的進(jìn)一步量化與解釋，使得區(qū)域土地利用變化更加明細(xì)化，公式如下[15-16]：

(2)

式中：Ilua為土地利用類型面積變化指數(shù)；uai,ubi分別為研究時段初年和末年第i類土地利用類型面積。

2.3 分級熱點探測

分級熱點探測是全局聚集性檢驗方法之一，用于確定在整個研究區(qū)域內(nèi)是否存在聚集性，通過定義一個“聚集單元”的“極限距離或閾值”，將其與每一個空間點對的距離進(jìn)行比較，距離小于極限距離時，該點被計入聚集單元[17]，在本研究中，采用定義極限距離進(jìn)行分析，同時也可以通過指定聚集單元點的數(shù)目來強(qiáng)化地理要素的聚集規(guī)則，以此類推得到研究區(qū)土地利用/覆被變化的熱點區(qū)域在空間上的轉(zhuǎn)移過程，是一種從整體上把握土地利用格局變化的直觀方法[18]。

3 結(jié)果與分析

3.1 土地利用現(xiàn)狀及轉(zhuǎn)移特征分析

結(jié)合表1、表2及土地利用現(xiàn)狀可以得出：在1991，2001，2013年中，林地一直是穆棱市主要的土地利用類型，分別占到72.22%，68.04%和71.27%，且主要分布于穆棱市南部和東部海拔較高的山區(qū)、半山區(qū)和丘陵地帶，集中連片分布，破碎度小。耕地面積在不同時期均占研究區(qū)總面積20%以上，集中分布在中部地勢平坦、水源豐富、居民點聚集的地區(qū)。

1991—2001年和2001—2013年土地利用轉(zhuǎn)移矩陣表明：1991—2001年，林地的面積不斷減少，共減少290.35 km2，其中主要轉(zhuǎn)化為耕地、草地和建設(shè)用地，轉(zhuǎn)化為耕地的面積占總比例的90%以上，同時也有少量的草地轉(zhuǎn)化為林業(yè)用地，建設(shè)用地面積小幅度增加，水域基本不變，這與同時期人類對研究區(qū)土地資源的不合理開發(fā)利用密切相關(guān)，在20世紀(jì)末，為提高糧食產(chǎn)量，研究區(qū)出現(xiàn)毀林開荒，亂砍濫伐現(xiàn)象，森林覆蓋面積急劇下降，也導(dǎo)致水土流失、土壤污染等災(zāi)害頻發(fā)，生態(tài)環(huán)境遭到破壞。2001年開始實施國家退耕還林工程，土地利用的格局和數(shù)量發(fā)生劇烈變化，呈現(xiàn)出新特點。2001年之后，由于國家政策干預(yù)，大量坡度在15°以上的、耕作條件差的耕地被還林還草，其中，旱地還林還草面積為171.51 km2，在研究區(qū)占主要地位。在此背景下，林地面積持續(xù)增加，截至2013年共增加171.45 km2，同時一些未利用地逐步演變?yōu)椴莸?，草地面積略有增加，水域面積在該時期增加約2 km2，由于農(nóng)村居民點的整治，建設(shè)用地的規(guī)范使用，使其利用率得到提高，騰出的空閑建設(shè)用地通過土地整治工程轉(zhuǎn)為其他用地，建設(shè)用地總面積有所下降，其開發(fā)使用更加科學(xué)合理。

表1 1991—2001年土地利用類型轉(zhuǎn)移矩陣 km2

表2 2001—2013年土地利用類型轉(zhuǎn)移矩陣 km2

3.2 土地利用變化特征分析

通過公式(1)和公式(2)，得出研究區(qū)1991—2013年間的土地利用總體變化指數(shù)和土地利用類型面積變化指數(shù)，并以區(qū)域內(nèi)交通、城鎮(zhèn)分布等信息為基礎(chǔ)，全面而準(zhǔn)確地反映了穆棱市土地利用/覆被變化特征。

根據(jù)土地利用總體變化指數(shù)計算結(jié)果，結(jié)合區(qū)域基礎(chǔ)信息可以看出，1991—2001年間土地利用的高變化區(qū)域主要集中在以穆棱市為中心、其西北部臨近雞西市區(qū)和沿鐵路公路交通干道向西南方向延伸的地帶，該地區(qū)地勢平坦，交通便利，人口聚集，耕地大面積連片分布，同時也是農(nóng)業(yè)與林業(yè)交錯地帶；2001—2013年間，土地利用變化中心明顯發(fā)生偏移，相比2001年前較分散的土地利用變化，2001年之后的土地利用變化情況呈現(xiàn)出中心聚集特征，城市已不是土地利用變化的最大區(qū)域，重心開始集中于距離城市比較近的中心鄉(xiāng)鎮(zhèn)，下城子鎮(zhèn)、馬橋河鎮(zhèn)和市區(qū)三者所圍成的區(qū)域內(nèi)，該區(qū)域是社會、經(jīng)濟(jì)和區(qū)域行政中心，多種優(yōu)勢要素聚集，未來發(fā)展空間大，同時也是鐵路、公路和水系交錯發(fā)達(dá)的地區(qū)，說明該區(qū)域的土地利用仍受交通的影響。

從土地利用類型面積變化指數(shù)計算結(jié)果可以得到，1991—2001年間，所有地類之間均有不同程度的相互轉(zhuǎn)化，其中林地轉(zhuǎn)耕地占轉(zhuǎn)化面積的50%以上，同時也有部分低覆蓋度的草地被開墾為耕地；其次在太平嶺和老爺嶺邊緣等難以開發(fā)利用地帶存在大比例的耕地轉(zhuǎn)為林草用地現(xiàn)象，但數(shù)量分布和所占比重小；這一時期在河流流域、水資源豐富地區(qū)旱地改水田的現(xiàn)象頻繁發(fā)生；建設(shè)用地變化在土地利用整體變化中十分微小，只略見于城市和鄉(xiāng)鎮(zhèn)邊緣區(qū)域。總之，該時期土地利用變化類型多種多樣。2001—2013年間，主要是耕地轉(zhuǎn)化為林草用地，且所占比重普遍很高；旱地改水田現(xiàn)象在數(shù)量和比重上均呈現(xiàn)增加趨勢；在局部地區(qū)出現(xiàn)草地轉(zhuǎn)為耕地和林地、建設(shè)用地增加等現(xiàn)象，但數(shù)量很少，且分布集中，主要分布于研究區(qū)中部地勢平坦區(qū)域及其邊緣的交錯地帶，這與該區(qū)域傳統(tǒng)的土地利用和社會發(fā)展分不開?？傊?，這一時期土地利用變化現(xiàn)象呈現(xiàn)出耕地轉(zhuǎn)為林草用地較多、其他地類間微小變動的特征。

3.3 熱點區(qū)域分析

借助于CrimeStat空間分析軟件，對穆棱市1991—2013年間兩個時間段的土地利用變化指數(shù)在同一參數(shù)標(biāo)準(zhǔn)下進(jìn)行分級熱點探測，得出1991—2001年和2001—2013年兩個時期土地利用變化的熱點區(qū)域空間分布圖(圖1)。1991—2001年間，穆棱市土地利用/覆被變化的熱點區(qū)域主要有3個，即“兩大一小”的區(qū)域，其中北部以市區(qū)為中心沿鐵路交通和穆棱河南北走向分布，同時也是研究區(qū)最大的熱點區(qū)域，其次為中部受國道交通影響較大，形成第二大熱點區(qū)域；穆棱市地處兩省交界處，在南部穆棱河源頭和吉林省交界地帶土地利用變化也較為突出。整體來看，該階段的土地利用變化并不受政策導(dǎo)向的影響，仍然是以肥沃的土地、豐富的水域和便利的交通等自然條件為先決要素，即自然要素決定土地資源的開發(fā)利用過程。2001—2013年間，穆棱市北部為土地利用/覆被變化的熱點區(qū)域，包括八面通鎮(zhèn)、馬橋河鎮(zhèn)和下城子鎮(zhèn)所形成的橢圓形區(qū)域，該區(qū)域包括了穆棱市郊和研究區(qū)的中心城鎮(zhèn)，同時是鐵路、公路、水系和地形等資源條件非常豐富的區(qū)域，構(gòu)成了研究區(qū)社會經(jīng)濟(jì)、交通發(fā)展的核心區(qū)域?？梢钥闯?，2001年之后的土地利用變化已經(jīng)受到社會發(fā)展、城市化等政策的干預(yù)，特別是近十年的新農(nóng)村建設(shè)、城鎮(zhèn)化發(fā)展更是加速了城市邊緣區(qū)、區(qū)域中心城鎮(zhèn)土地資源利用的劇烈變化，自然要素在土地利用中的主導(dǎo)地位已逐漸弱化?？傊?，隨著社會經(jīng)濟(jì)的發(fā)展，穆棱市土地資源的開發(fā)利用從亂開濫用轉(zhuǎn)向受政策導(dǎo)向的重點開發(fā)利用，土地資源得到高效合理配置。

圖1 1991—2013年穆棱市土地利用變化熱點圖

3.4 驅(qū)動力分析

3.4.1 氣候因子。氣候變化一般被認(rèn)為是土地利用覆被變化的主要驅(qū)動因素之一[16,19]。選取穆棱市1991—2012年間的年降水量、≥0.1 mm降水日數(shù)、年平均氣溫和平均相對濕度作為研究區(qū)自然驅(qū)動力因子進(jìn)行統(tǒng)計分析。從而得到穆棱市近22年間的降水量以600 mm線為中軸上下波動，在1999年和2011年出現(xiàn)較大波動點；≥0.1 mm降水日數(shù)和年平均氣溫在22年間均無較大的變化，相對穩(wěn)定；22年間穆棱市年平均相對濕度發(fā)生十分微小的波動，其值始終在65%左右。綜合幾項氣候因子在時間序列上的變化情況，穆棱市近22年間自然要素基本上沒有發(fā)生變化，也沒有氣溫驟降、降雨量突變等極端氣候的出現(xiàn)?？梢钥闯?，相對于穆棱市近22年的土地利用/覆被變化特征，氣溫、降雨等自然驅(qū)動力對其的影響是十分微弱的。

3.4.2 社會因子。人類活動和社會發(fā)展作為影響自然過程的主要要素受到廣泛的關(guān)注[20-22]。人類活動對研究區(qū)土地利用/覆被變化的影響具有時間段、變化大的特點。穆棱市近些年人口總數(shù)突破30萬，全市生產(chǎn)總值在2010年已超過100億元，人口、社會經(jīng)濟(jì)增長迅速；在土地利用方面，2001年前，相關(guān)政策規(guī)范尚未形成，人類活動受政策影響較少，對土地資源的開發(fā)利用依自然條件的優(yōu)劣而行。研究區(qū)經(jīng)歷多年隨意性的混亂開發(fā)，2001年初，穆棱市作為試點市開始實施退耕還林工程，坡度>15°耕地堅決退耕還林還草，中低農(nóng)業(yè)用地得到改善，不適宜農(nóng)業(yè)用地退出農(nóng)業(yè)用地，撂荒地進(jìn)行生態(tài)林營造，嚴(yán)格限制耕地面積和建設(shè)用地占用，農(nóng)民在這個過程中也得到國家經(jīng)濟(jì)補(bǔ)助，退耕還林，保護(hù)生態(tài)環(huán)境的熱情空前高漲。政策的實施改變了以往的土地利用情形，土地利用結(jié)構(gòu)在短期內(nèi)得到優(yōu)化調(diào)整，林業(yè)用地增加，耕地資源合理配置，建設(shè)用地限制性發(fā)展。由人類活動對土地利用產(chǎn)生的干擾隨著社會經(jīng)濟(jì)的發(fā)展更加趨于規(guī)范性和方向性，經(jīng)濟(jì)的發(fā)展促使越來越多的人口從第一產(chǎn)業(yè)向二、三產(chǎn)業(yè)轉(zhuǎn)移，土地利用的零碎狀態(tài)與盲目性逐漸消失，促進(jìn)了土地利用的規(guī)模化與規(guī)范性。同時，隨著科學(xué)發(fā)展理念的廣泛宣傳，人們環(huán)保意識不斷提高，河流湖泊等具有重要生態(tài)意義的土地類型和地表過程得到保護(hù)，經(jīng)過這一時期的調(diào)整發(fā)展，穆棱市的土地資源利用在結(jié)構(gòu)和數(shù)量上達(dá)到最佳狀態(tài)?？梢钥闯觯祟惿a(chǎn)活動和國家政策的實施是穆棱市過去幾十年土地利用/覆被變化的主要驅(qū)動因子。

4 結(jié)論

1)1991—2013年，穆棱市土地利用/覆被變化呈現(xiàn)出不同的空間特征和演變過程。以2001年為研究分割點，在實施退耕還林前，由于缺乏相關(guān)政策和可持續(xù)發(fā)展的意識，穆棱市土地利用呈現(xiàn)攤子大、土地資源開發(fā)利用無規(guī)律的特點，耕地突增，林草用地驟減，生態(tài)環(huán)境遭到破壞并惡化。退耕還林后，通過坡耕地還林還草，嚴(yán)格建設(shè)用地占用耕地等政策，耕地與林地之間轉(zhuǎn)化比例大，研究區(qū)森林覆蓋率得到提高，土地資源利用的適宜性進(jìn)一步提高，整體土地利用發(fā)生變化，但劇烈變化區(qū)域相對集中，多發(fā)生在城市和中心鄉(xiāng)鎮(zhèn)的周邊區(qū)域，區(qū)域內(nèi)中心鄉(xiāng)鎮(zhèn)的發(fā)展是新形勢下城鎮(zhèn)化發(fā)展的重點。

2)針對大時間尺度上的氣候變化和人類活動對土地利用/覆被變化的影響，在全球范圍、大流域尺度內(nèi)對氣候變化與土地利用/覆被變化的研究具有重大意義。從小空間尺度出發(fā)進(jìn)行同類研究，結(jié)果表明氣候等自然因素對小區(qū)域土地利用/覆被變化的影響十分有限，特別是氣候變化相對平穩(wěn)的區(qū)域范圍。本研究從1991年到2013年近23年的時間尺度分析表明，人類活動是造成區(qū)域土地利用發(fā)生劇烈變化和重心轉(zhuǎn)移的根本原因，特別是國家政策、人口增長和經(jīng)濟(jì)發(fā)展等人文因素。

3)研究表明，以公里網(wǎng)格為單元進(jìn)行土地類型面積變化指數(shù)的計算分析并進(jìn)行空間展布意義重大。在未來區(qū)域土地利用變化研究中，通過結(jié)合區(qū)域發(fā)展情況對小單元內(nèi)土地類型面積變化指數(shù)擬合分析，由點及面對區(qū)域內(nèi)土地利用變化做出模擬分析，以此為基礎(chǔ)，能夠進(jìn)一步準(zhǔn)確有效地補(bǔ)充和完善中小尺度的不同發(fā)展情景下土地利用/覆被變化的模擬預(yù)測及相關(guān)研究，也是土地科學(xué)領(lǐng)域未來研究的方向和重點。

[1] 后立勝，蔡運龍.土地利用/覆被變化研究的實質(zhì)分析與進(jìn)展評述[J].地理科學(xué)進(jìn)展，2004，23(6)：96-104.

[2] Su C H,Fu B J,Lu Y H,etal.Land Use Change and Anthropogenic Driving Force:A Case Study in Yanhe River Basin[J].Chinese Geographical Science,2011,21(5):587-599.

[3] Zhang K,Yu Z,Li X,etal.Land Use Change and Land Degradation in China from 1991 to 2001[J].Land Degradation and Development,2007，18(2):209-219.

[4] 劉紀(jì)遠(yuǎn)，匡文慧，張增祥，等.20世紀(jì)80年代末以來中國土地利用變化的基本特征與空間格局[J].地理學(xué)報，2014，69(1)：3-14.

[5] Liu J Y,Zhang Z X,Xu X L，etal.Spatial Patterns and Driving Forces of Land Use Change in China during the Early 21st Century[J].Journal of Geographical Sciences,2010,20(4):483-494.

[6] 李偲，韓桂紅，海米提·依米提，等.土地利用變化對喀納斯自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價值的影響[J].地域研究與開發(fā)，2011,30(3)：123-127.

[7] 李傳哲，于福亮，劉佳，等.近20年來黑河干流中游地區(qū)土地利用/覆被變化及驅(qū)動力定量研究[J].自然資源學(xué)報，2011，26(3)：353-363.

[8] 高鳳杰，雷國平,楊鳳海，等.牡丹江市退耕還林前后生態(tài)安全評價[J].農(nóng)業(yè)工程學(xué)報，2012，28(2)：206-214.

[9] 李小雁，許何也，馬育軍，等.青海湖流域土地利用/覆被變化研究[J].自然資源學(xué)報，2008，23(2)：285-296.

[10] 胡硯霞，黃進(jìn)良，王立輝.丹江口庫區(qū)1990—2010年土地利用時空動態(tài)變化研究[J].地域研究與開發(fā)，2013,32(3):133-137.

[11] 張偉科，封志明，楊艷昭，等.北方農(nóng)牧交錯帶土地利用/覆被變化分析——以西遼河流域為例[J].資源科學(xué)，2010，32(3)：573-579.

[12] 張繼平，常學(xué)禮，李健英，等.內(nèi)蒙古奈曼旗農(nóng)牧交錯區(qū)土地利用/覆被變化的區(qū)域分析[J].應(yīng)用生態(tài)學(xué)報，2008,19(3)：613-620.

[13] 戚伯華.調(diào)整土地利用結(jié)構(gòu)，構(gòu)建“大牡丹江”生態(tài)框架[J].國土與自然資源研究，2009(3)：40-42.

[14] 張連生.重整山河蔭后人——穆棱市退耕還林巡禮[J].中華兒女，2001(9)：62.

[15] 王介勇，劉彥隨.三亞市土地利用/覆被變化及其驅(qū)動機(jī)制研究[J].自然資源學(xué)報，2009，24(8)：1458-1464.

[16] 田錫文，王新軍，卡迪羅夫，等.近40 a凱拉庫姆庫區(qū)土地利用/覆蓋變化及景觀格局分析[J].農(nóng)業(yè)工程學(xué)報，2014，30(6)：232-241.

[17] 王勁峰，廖一蘭，劉鑫.空間數(shù)據(jù)分析教程[M].北京：科學(xué)出版社，2010.

[18] 王勁峰，葛詠，李連發(fā)，等.地理學(xué)時空數(shù)據(jù)分析方法[J].地理學(xué)報，2014，69(9)：1326-1345.

[19] 馬晴，李丁，廖杰，等.疏勒河中下游綠地土地利用變化及其驅(qū)動力分析[J].經(jīng)濟(jì)地理，2014，34(1)：148-155.

[20] 迪力夏提·阿卜杜薩拉木，阿不都克依木·阿布力孜，周玄德，等.干旱區(qū)典型綠洲耕地擴(kuò)張及其人為驅(qū)動研究——以新疆且末綠洲為例[J].地域研究與開發(fā)，2015，34(2)：131-136.

[21] 王云霞，南靈.西安市城鎮(zhèn)化進(jìn)程中人口與土地協(xié)調(diào)發(fā)展研究[J].地域研究與開發(fā)，2015，34(2)：85-90.

[22] 謝花林.典型農(nóng)牧交錯區(qū)土地利用變化驅(qū)動力分析[J].農(nóng)業(yè)工程學(xué)報，2008，24(10)：56-62.

Land Use/Cover Change and Hot Spots Analysis in Muling City of Heilongjiang Province

Gao Fengjie , Ma Quanlai , Shan Peiming ,Liu Yang , Guo Xinxin , Yang Shuang , Zhang Zhimin

(SchoolofResourceandEnvironment,NortheastAgriculturalUniversity,Harbin150030，China)

Land Use/Cover Change(LUCC) and its process characteristics is an important and hot issue in global change research filed. This study focused on spatial-temporal evolution of land use and cover change from 1991 to 2013 in Muling City of Heilongjiang Province, using geo-statistics and 3S technology. From 1991 to 2001, the area of cropland increased 232.25 km2, the forest decreased 290.35 km2and the changes in cropland and forest accounted were the big parts in the land use and cover change, others also had a change in area, but small and scattered, the exploit chaos and low efficiency were the characteristic of land use from 1991 to 2001 in Muling. From 2001 to 2013, the area of cropland decreased 179.01 km2, the forest increased 171.45 km2, and the exit of slop land had a great contribution for the change and a wide distribution, the structure of land use optimized in the period. The hot spots of land use was a key region consisted of city and key towns during 2001—2013 from the pattern of two big and a small regions during 1991—2001. From 2001 to 2013, nature factor changed slightly in Muling, land use and cover change was drived by the policy of Grain for Green and other land use policies.

land use/land cover change; geo-statistics; 3S technology; Muling City of Heilongjiang Province

2015-03-31；

2016-05-27

黑龍江省教育廳科學(xué)技術(shù)面上研究項目(12541037)

高鳳杰(1981-)，女，河北遷西縣人，副教授，博士， 主要從事資源環(huán)境遙感研究，(E-mail)gaojiecumt@126.com。

郭欣欣(1982-)，女，黑龍江齊齊哈爾市人，講師，博士，主要從事土地利用規(guī)劃研究，(E-mail)xinxin8504021@163.com。

F301.24

A

1003-2363(2016)04-0126-05