金輝，李晉陵，韓偉宏，邢虹娟，郭軍玲

（1.山西省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與經(jīng)濟(jì)研究所，山西太原030006；2.山西省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與資源研究所，山西太原030031）

基于GIS的土地利用變化對(duì)礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響

金輝1，李晉陵1，韓偉宏1，邢虹娟1，郭軍玲2

（1.山西省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與經(jīng)濟(jì)研究所，山西太原030006；2.山西省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與資源研究所，山西太原030031）

運(yùn)用GIS技術(shù)對(duì)山西繁峙縣金方圓礦區(qū)土地復(fù)墾方案（2013—2028年）的土地利用數(shù)據(jù)進(jìn)行處理，對(duì)研究區(qū)進(jìn)行生態(tài)系統(tǒng)服務(wù)價(jià)值估算，分析了礦區(qū)土地利用類型變化及其所引起的生態(tài)系統(tǒng)服務(wù)價(jià)值的變化。結(jié)果表明，2013—2028年間土地利用類型變化較大，其他草地、采礦用地、其他林地和裸地呈減少趨勢(shì)，有林地、灌木林地、人工牧草地呈增加趨勢(shì)；土地利用類型變化主要是其他草地和采礦用地向人工牧草地、灌木林地和有林地轉(zhuǎn)變，其中，以其他草地和采礦用地減少為主，人工牧草地呈平穩(wěn)增加趨勢(shì)；期間生態(tài)系統(tǒng)服務(wù)價(jià)值呈現(xiàn)先下降后升高的趨勢(shì)，研究區(qū)內(nèi)生態(tài)系統(tǒng)服務(wù)價(jià)值占比以林地和草地為主。礦區(qū)土地復(fù)墾方案的實(shí)施，驅(qū)使土地利用類型向著耕地、有林地和草地等生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)較高的地類轉(zhuǎn)變，不斷提高礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值總量，對(duì)保障區(qū)域生態(tài)安全和生態(tài)系統(tǒng)可持續(xù)發(fā)展具有重要意義。

GIS；土地利用；礦區(qū)；生態(tài)系統(tǒng)服務(wù)；價(jià)值

礦產(chǎn)資源是社會(huì)發(fā)展和經(jīng)濟(jì)建設(shè)的重要物質(zhì)基礎(chǔ)，礦產(chǎn)資源的合理、高效開發(fā)利用極大地推動(dòng)了社會(huì)經(jīng)濟(jì)發(fā)展，但伴隨而來(lái)的礦區(qū)土地退化和環(huán)境惡化等給生態(tài)環(huán)境帶來(lái)了巨大的壓力[1-2]。土地復(fù)墾能夠有效、快速且明顯地恢復(fù)和提高土地可持續(xù)利用能力，提高土地生產(chǎn)率和利用率，保護(hù)和增值自然資源，進(jìn)而增強(qiáng)經(jīng)濟(jì)、生態(tài)、社會(huì)綜合效益[3]。土地復(fù)墾規(guī)劃是以最大程度恢復(fù)礦山生態(tài)系統(tǒng)功能為根本目標(biāo)，對(duì)礦山生態(tài)恢復(fù)起到指導(dǎo)和約束的作用，同時(shí)規(guī)劃的實(shí)施將有效提升礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值，極大地緩解礦區(qū)土地破壞給環(huán)境帶來(lái)的壓力[4]。Costanza等較早對(duì)全球生態(tài)系統(tǒng)服務(wù)與自然資本價(jià)值估算進(jìn)行了研究，推動(dòng)了關(guān)于生態(tài)系統(tǒng)服務(wù)在全球的廣泛研究，我國(guó)眾多學(xué)者也分別對(duì)青藏高原、陸地、平原等不同生態(tài)系統(tǒng)進(jìn)行了評(píng)價(jià)[5-10]。生態(tài)系統(tǒng)服務(wù)功能是指自然生態(tài)系統(tǒng)所提供的可以滿足和維持人類生活所需的條件和過(guò)程[11-13]，土地利用/覆蓋變化能改變生態(tài)系統(tǒng)的結(jié)構(gòu)、過(guò)程和功能，進(jìn)而影響生態(tài)系統(tǒng)服務(wù)[14-16]，因此，針對(duì)山西工礦區(qū)待復(fù)墾面積大、復(fù)墾率低的現(xiàn)狀[17-18]，可對(duì)復(fù)墾前后不同時(shí)期生態(tài)系統(tǒng)服務(wù)進(jìn)行定量測(cè)算以反映礦區(qū)土地利用結(jié)構(gòu)的變化以及礦山生態(tài)環(huán)境治理所取得的效果，引導(dǎo)生態(tài)環(huán)境朝著生態(tài)系統(tǒng)服務(wù)價(jià)值提升的土地利用方向進(jìn)行恢復(fù)治理，對(duì)保障區(qū)域生態(tài)安全和可持續(xù)發(fā)展具有至關(guān)重要意義，促進(jìn)縣域經(jīng)濟(jì)與生態(tài)環(huán)境健康可持續(xù)發(fā)展[19-20]。

本研究以忻州市繁峙縣金方圓礦區(qū)為研究對(duì)象，應(yīng)用土地利用現(xiàn)狀及復(fù)墾規(guī)劃資料，定量計(jì)算和分析礦區(qū)復(fù)墾過(guò)程中土地利用變化所導(dǎo)致的生態(tài)系統(tǒng)服務(wù)價(jià)值的變化，從生態(tài)經(jīng)濟(jì)價(jià)值的角度研究礦區(qū)生態(tài)復(fù)墾為環(huán)境帶來(lái)的貢獻(xiàn)，為礦區(qū)土地生態(tài)恢復(fù)治理提供科學(xué)決策依據(jù)。

1 研究區(qū)域概況

繁峙縣金方圓礦業(yè)有限公司鐵礦位于繁峙縣城214°方向，直距約15 km處，官地村西南3.5 km處，行政區(qū)劃隸屬于繁峙縣巖頭鄉(xiāng)管轄。礦區(qū)中心地理坐標(biāo)為東經(jīng)113°15′13″，北緯39°04′55″，面積為2.49 km2。屬暖溫帶半干旱季風(fēng)氣候，四季分明，夏秋濕潤(rùn)多雨，冬季干燥寒冷。多年平均氣溫6.9℃，最高氣溫35.5℃，最低氣溫-32.5℃；多年平均降水量565.8 mm，最大年降水量為914.4 mm，最小年降水量為265.7 mm，降水主要集中在7—9月，約占全年降水量的70%；年蒸發(fā)量為1 600～1 800 mm，年日照2700～2900h，≥10℃積溫2600～3500℃，無(wú)霜期120 d。地勢(shì)高峻而起伏不平，地形切割強(qiáng)烈，土壤分布主要類型為中性粗骨土和石灰性褐土。項(xiàng)目區(qū)所在巖頭鄉(xiāng)總面積315 km2，2014年巖頭鄉(xiāng)財(cái)政收入4 654萬(wàn)元，人均收入4 351元，人均耕地0.16 hm2。

2 數(shù)據(jù)來(lái)源與研究方法

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

本研究以《繁峙縣金方圓礦業(yè)有限公司鐵礦項(xiàng)目土地復(fù)墾方案》（2013—2028年）、2013年度繁峙縣土地利用現(xiàn)狀圖、遙感影像、野外實(shí)地調(diào)查為主要數(shù)據(jù)來(lái)源，同時(shí)根據(jù)需要搜集的繁峙縣2010—2014年糧食總產(chǎn)量、糧食作物播種面積、農(nóng)業(yè)總產(chǎn)值等社會(huì)經(jīng)濟(jì)數(shù)據(jù)來(lái)源于《山西省統(tǒng)計(jì)年鑒》（2011—2015年）。

2.2 研究方法

2.2.1 土地利用動(dòng)態(tài)變化土地利用動(dòng)態(tài)模型反映研究時(shí)間域內(nèi)土地利用類型面積的變化速度。即在某一段時(shí)期內(nèi)，某種土地利用類型面積變化量與初始值的比值。其表達(dá)式如下。

其中，K為研究時(shí)間區(qū)段內(nèi)土地利用類型時(shí)間區(qū)間變化率；Ua，Ub分別為研究時(shí)間區(qū)段初期和末期的土地利用類型面積。本研究中以時(shí)間區(qū)間為單位，所以，K值為該研究區(qū)域某種土地利用類型面積的區(qū)間變化強(qiáng)度。

為反映某土地利用類型的變化劇烈程度，采用公式（2）表示某一土地類型變化的動(dòng)態(tài)度（Rss）。

其中，ΔUin為在一定時(shí)期內(nèi)其他地類轉(zhuǎn)變?yōu)樵擃愋屯恋氐拿娣e之和；ΔUout為在一定時(shí)期內(nèi)某一地類轉(zhuǎn)變?yōu)槠渌仡惖拿娣e之和；Rss≥0，當(dāng)Rss趨于0時(shí)，說(shuō)明在研究時(shí)段內(nèi)，該類型土地面積變化極小，Rss數(shù)值越大，表明該類型土地面積的變化程度愈劇烈。

2.2.2 土地生態(tài)系統(tǒng)服務(wù)價(jià)值參考Costanza的生態(tài)系統(tǒng)服務(wù)價(jià)值計(jì)算公式，結(jié)合繁峙縣社會(huì)經(jīng)濟(jì)數(shù)據(jù)和土地利用等相關(guān)資料，對(duì)生態(tài)服務(wù)價(jià)值進(jìn)行參數(shù)修正，計(jì)算得出礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值。

式中，ESV代表生態(tài)服務(wù)價(jià)值；AR代表R類土地利用類型的面積（hm2）；VCR為土地利用類型的生態(tài)服務(wù)系數(shù)（元/（hm2·a））。

2.3 技術(shù)路線

本研究參考Costanza對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的估算模型，結(jié)合社會(huì)經(jīng)濟(jì)數(shù)據(jù)、礦區(qū)土地利用數(shù)據(jù)和生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估體系，引用謝高地等[21]對(duì)中國(guó)生態(tài)系統(tǒng)服務(wù)價(jià)值估算當(dāng)量表，將重新分類的土地利用類型與Costanza等研究全球生態(tài)系統(tǒng)分類進(jìn)行比對(duì)，對(duì)單位面積各地類生態(tài)系統(tǒng)服務(wù)功能進(jìn)行估價(jià)，即確定項(xiàng)目區(qū)各類土地利用類型的生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)，得出區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值多元一次函數(shù)關(guān)系式；利用遙感影像、空間數(shù)據(jù)和ArcGIS軟件對(duì)各項(xiàng)空間數(shù)據(jù)和屬性數(shù)據(jù)進(jìn)行處理，建立礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)價(jià)數(shù)據(jù)庫(kù)，得出礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值匯總數(shù)據(jù)及相關(guān)圖件，以便進(jìn)一步分析研究區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值的變化情況，評(píng)估復(fù)墾規(guī)劃目標(biāo)對(duì)生態(tài)系統(tǒng)的貢獻(xiàn)價(jià)值量，技術(shù)路線如圖1所示

。

3 結(jié)果與分析

3.1 土地利用狀況分析

由表1和圖2可知，項(xiàng)目區(qū)土地現(xiàn)狀以其他草地和采礦用地為主，占項(xiàng)目區(qū)總面積的67.08%，其次是林地、裸地和農(nóng)村道路，同時(shí)有少量耕地、村莊分布。根據(jù)復(fù)墾方案的規(guī)劃，2013—2028年間土地利用方式轉(zhuǎn)變最主要的是：（1）其他草地、采礦用地減少，其中，其他草地面積減少31.07 hm2，采礦用地面積減少12.58hm2，其他林地和裸地減少面積較小，僅為3.48，0.31 hm2；（2）人工牧草地、灌木林地和有林地增加，其中，人工牧草地面積增加35.80 hm2，有林地和灌木林地面積分別增加了4.79，6.83 hm2。通過(guò)規(guī)劃實(shí)施，對(duì)礦區(qū)損毀土地復(fù)墾，在保證項(xiàng)目實(shí)施后耕地面積不減少的基本原則下，將已損毀土地根據(jù)復(fù)墾適宜性恢復(fù)為林地和人工牧草地，有效增加項(xiàng)目區(qū)生態(tài)綠地面積。

表12013 —2028年礦區(qū)土地利用的年際變化hm2

在礦區(qū)生產(chǎn)過(guò)程中實(shí)施“邊生產(chǎn)、邊治理”的策略，2013—2019年以土地占用為主，復(fù)墾工作逐步開展，采礦用地面積增加27.18 hm2用于礦區(qū)開采作業(yè)，因土地占用，其他草地面積減少31.07 hm2，其他林地面積減少3.48 hm2，旱地和裸地面積分別減少0.22，0.31 hm2，期間通過(guò)恢復(fù)治理，人工牧草地增加7.70hm2，灌木林地面積增加0.19hm2；2019—2025年土地利用方式轉(zhuǎn)變因工程復(fù)墾而由建設(shè)用地向生態(tài)用地轉(zhuǎn)變，期間采礦用地面積減少了20.39 hm2，復(fù)墾為人工牧草地面積14.69 hm2，有林地和灌木林地面積分別增加2.29，3.42 hm2；2025—2028年為項(xiàng)目區(qū)復(fù)墾收官階段，所有損毀土地都將恢復(fù)為農(nóng)用地，本階段采礦用地面積再減少19.37 hm2，恢復(fù)旱地面積0.22 hm2，有林地面積2.50 hm2，灌木林地3.22 hm2，人工牧草地13.41 hm2，至此因礦區(qū)開采作業(yè)損毀土地均已完成生態(tài)恢復(fù)，主要以人工牧草地為主要恢復(fù)地類，有效增加了礦區(qū)植被覆蓋。

3.2 土地利用變化率和動(dòng)態(tài)度分析

從圖3可以看出，2013—2019年，礦區(qū)采礦用地面積變化率最大，為45.3%，其次是其他草地，為-29.02%，旱地、灌木林地、其他林地和裸地變化率較??；2019—2025年，礦區(qū)地類面積變化率最大的是人工牧草地，為190%，灌木林地和采礦用地面積變化率分別為15.72%和-23.39%；2025—2028年，礦區(qū)地類面積變化率最大的是有林地，為104.60%，其次是人工牧草地，為59.89%，采礦用地面積變化率為-29%，灌木林地和旱地變化幅度較小，分別為12.79 %和4.99%。

2013—2028年礦區(qū)土地利用類型面積變化速度和地類相關(guān)性較大，復(fù)墾初地類變化速度較快的是采礦用地，動(dòng)態(tài)度達(dá)到7.55；隨著規(guī)劃的實(shí)施，2019—2025年間地類變化速度較快的是人工牧草地，動(dòng)態(tài)度為31.80；規(guī)劃實(shí)施末期地類向生態(tài)用地轉(zhuǎn)變，變化速度較快的是有林地，動(dòng)態(tài)度達(dá)到34.87，其次是人工牧草地，為19.96，其他地類變化動(dòng)態(tài)度較小。整個(gè)復(fù)墾期間建設(shè)用地的變化除其他草地和裸地因礦區(qū)開采作業(yè)階段損毀外，而后地類再?zèng)]有發(fā)生變化。

3.3 生態(tài)系統(tǒng)服務(wù)價(jià)值分析

為了反映礦區(qū)土地利用類型對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響，利用2010—2012年繁峙縣糧食作物播種面積3.65萬(wàn)hm2，年均糧食產(chǎn)量7.12萬(wàn)t，年均農(nóng)業(yè)總產(chǎn)值5 651.69元/hm2，生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子則參考謝高地等[21]的中國(guó)陸地生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值當(dāng)量表，1個(gè)生態(tài)服務(wù)價(jià)值當(dāng)量因子的經(jīng)濟(jì)價(jià)值量等于當(dāng)年全國(guó)平均糧食單產(chǎn)市場(chǎng)價(jià)值的1/7，同時(shí)根據(jù)土地利用現(xiàn)狀分類定義，對(duì)林地、草地的生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行系數(shù)修訂，計(jì)算出繁峙縣1個(gè)生態(tài)服務(wù)價(jià)值當(dāng)量因子的經(jīng)濟(jì)價(jià)值為807.38元，基于此建立多元一次函數(shù)式來(lái)推算得出繁峙縣生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值表（表2），利用地理信息系統(tǒng)構(gòu)建數(shù)據(jù)庫(kù)，得出礦區(qū)復(fù)墾前、后生態(tài)系統(tǒng)服務(wù)價(jià)值單位面積分布（圖4）。

表2 繁峙縣生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值元/（hm2·a）

根據(jù)繁峙縣生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值表和土地利用面積計(jì)算出了復(fù)墾方案實(shí)施階段（2013，2019，2025，2028年）礦區(qū)的生態(tài)系統(tǒng)服務(wù)價(jià)值。由表3可知，2013—2028年間生態(tài)系統(tǒng)服務(wù)價(jià)值呈先降低后增加的趨勢(shì)，即2013年因礦區(qū)開采作業(yè)階段損毀部分土地，導(dǎo)致礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值總量降低，而后隨著復(fù)墾規(guī)劃實(shí)施，高生態(tài)系統(tǒng)服務(wù)價(jià)值土地利用類型土地面積逐漸增加，2013年礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值為110.73萬(wàn)元，2028年礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值為131.07萬(wàn)元，年均增加1.36萬(wàn)元。其中，農(nóng)田的生態(tài)系統(tǒng)服務(wù)價(jià)值保持不變，其占比由2013年的2.33%減少到2028年的1.97%；林地作為礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值唯一增加的土地利用類型，增加了13.18萬(wàn)元，占比由原來(lái)的57.54%提高到2028年的58.66%，林地生態(tài)系統(tǒng)服務(wù)價(jià)值增加的主要原因是：一方面源自于林地面積的增加，另一方面林地是礦區(qū)內(nèi)對(duì)生態(tài)系統(tǒng)服務(wù)功能貢獻(xiàn)最大的土地利用類型；荒漠生態(tài)系統(tǒng)服務(wù)價(jià)值因面積減少0.31 hm2而減少0.01萬(wàn)元；2013年草地生態(tài)系統(tǒng)服務(wù)價(jià)值占比是39.57%，隨著規(guī)劃的實(shí)施，草地生態(tài)系統(tǒng)服務(wù)價(jià)值占比呈現(xiàn)先降低后增加的趨勢(shì)，2028年草地生態(tài)系統(tǒng)服務(wù)價(jià)值占比為38.90%；2013—2028年農(nóng)田生態(tài)系統(tǒng)服務(wù)價(jià)值因?yàn)槠渌仡惐徽加脤?dǎo)致占比提高，隨著復(fù)墾的實(shí)施其他高生態(tài)價(jià)值地類得到恢復(fù)補(bǔ)充，使得農(nóng)田占礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值比例下降至1.97%。

表32013 ，2019，2025，2028年礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值

4 結(jié)論

本研究對(duì)《繁峙縣金方圓礦業(yè)有限公司鐵礦項(xiàng)目土地復(fù)墾方案》中土地利用類型的變化情況進(jìn)行評(píng)估，結(jié)果顯示，2013—2028年礦區(qū)土地利用類型總體變化是由損毀土地向宜耕宜農(nóng)的方向轉(zhuǎn)變，其中變化最大的是其他草地和采礦用地?fù)p毀后恢復(fù)為人工牧草地；從變化過(guò)程來(lái)看，由初期項(xiàng)目占用各個(gè)地類的損毀導(dǎo)致土地利用類型向建設(shè)用地的轉(zhuǎn)移，復(fù)墾方案實(shí)施后進(jìn)入生態(tài)系統(tǒng)恢復(fù)階段，土地利用類型根據(jù)工程安排，逐步開始向著旱地、林地、草地等生態(tài)系統(tǒng)服務(wù)價(jià)值較高的土地利用類型轉(zhuǎn)變，這也符合國(guó)家對(duì)礦業(yè)用地恢復(fù)治理的基本原則，因此，提高土地利用水平是恢復(fù)和提高區(qū)域生態(tài)系統(tǒng)服務(wù)功能的最有效途徑。

礦區(qū)土地復(fù)墾方案的實(shí)施將大幅提升礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值，由2013年的110.73萬(wàn)元增加到2028年的131.07萬(wàn)元，增幅達(dá)18.36%，年均增加1.36萬(wàn)元，其中為礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值貢獻(xiàn)最大的是林地，面積增加8.14 hm2，其生態(tài)系統(tǒng)服務(wù)價(jià)值增加13.18萬(wàn)元，占比增加1.12百分點(diǎn)。草地生態(tài)系統(tǒng)服務(wù)價(jià)值增加7.17萬(wàn)元，為礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值貢獻(xiàn)占比減少0.66百分點(diǎn)。礦區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值的提升得益于建設(shè)用地和荒漠的損毀后的生態(tài)恢復(fù)，通過(guò)工程、技術(shù)措施和資金投入，使單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值較低的建設(shè)用地和未利用地轉(zhuǎn)變?yōu)檗r(nóng)用地。

基于時(shí)空動(dòng)態(tài)監(jiān)測(cè)評(píng)估的區(qū)域生態(tài)系統(tǒng)服務(wù)功能和價(jià)值的評(píng)價(jià)，是評(píng)估區(qū)域生態(tài)系統(tǒng)服務(wù)真實(shí)價(jià)值的基礎(chǔ)和依據(jù)，有助于從科學(xué)角度對(duì)礦區(qū)土地復(fù)墾方案進(jìn)行評(píng)估，可作為復(fù)墾方案的可行性和合理性的重要指標(biāo)。

［1］邱文瑋.礦區(qū)生態(tài)服務(wù)功能價(jià)值的評(píng)估模型及其應(yīng)用研究：以徐州九里礦區(qū)為例[D].北京：中國(guó)礦業(yè)大學(xué)，2014.

［2］張耿杰，白中科，喬麗，等.平朔礦區(qū)生態(tài)系統(tǒng)服務(wù)功能價(jià)值變化研究[J].資源與產(chǎn)業(yè)，2008（6）：8-12.

［3］劉慧芳，洪堅(jiān)平.礦區(qū)損毀土地復(fù)墾適宜性評(píng)價(jià)：以靈石縣夏門煤礦為例[J].山西農(nóng)業(yè)科學(xué)，2014，42（3）：272-275．

［4］邢育剛.煤礦區(qū)地面沉陷引發(fā)的生態(tài)服務(wù)價(jià)值變化與生態(tài)修復(fù)對(duì)策研究[D].太原：山西大學(xué)，2013.

［5］余嘉琦，李鋼，趙華，等.江蘇省沛縣土地利用變化及其生態(tài)服務(wù)價(jià)值研究[J].江蘇農(nóng)業(yè)科學(xué)，2015，43（6）：371-376.

［6］鄭江坤，余新曉，賈國(guó)棟，等.密云水庫(kù)集水區(qū)基于LUCC的生態(tài)服務(wù)價(jià)值動(dòng)態(tài)演變[J].農(nóng)業(yè)工程學(xué)報(bào)，2010，26（9）：315-320．

［7］肖強(qiáng)，肖洋，歐陽(yáng)志云，等.重慶市森林生態(tài)系統(tǒng)服務(wù)功能價(jià)值評(píng)估[J].生態(tài)學(xué)報(bào)，2014，34（1）：216-223.

［8］程琳，李鋒，鄧華鋒.中國(guó)超大城市土地利用狀況及其生態(tài)系統(tǒng)服務(wù)動(dòng)態(tài)演變[J].生態(tài)學(xué)報(bào)，2011，31（20）：6194-6203.

［9］陳春陽(yáng)，戴君虎，王煥炯，等.基于土地利用數(shù)據(jù)集的三江源地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值變化[J].地理科學(xué)進(jìn)展，2012（7）：970-977.

［10］石龍宇，崔勝輝，尹鍇，等.廈門市土地利用/覆被變化對(duì)生態(tài)系統(tǒng)服務(wù)的影響[J].地理學(xué)報(bào)，2010（6）：708-714.

［11］王偉，陸健健.生態(tài)系統(tǒng)服務(wù)功能分類與價(jià)值評(píng)估探討[J].生態(tài)學(xué)雜志，2005，24（11）：1314-1316.

［12］Costanza R.The value of the world's ecosystem services and natural capital[J].Nature，1997，387：253-260.

［13］Ouyang Z Y，Wang R S，Zhao J Z.Ecosystem services and their economic valuation[J].Chinese Journal of Applied Ecology，1999，10（5）：635-640.

［14］Priess J A，Mimler M，Klein A M，et al.Linking deforestation scenarios to pollination services and economic returns in coffee agroforestry systems[J].Ecological Applications，2007，17：407-417.

［15］Ricketts TH，RegetzJ，Steffan-Dewenter I，et al.Landscape effects on crop pollination services:are there general patterns?[J].Ecology Letters，2008，11：499-515.

［16］吳海珍，阿如旱，郭田保，等.基于RS和GIS的內(nèi)蒙古多倫縣土地利用變化對(duì)生態(tài)服務(wù)價(jià)值的影響[J].地理科學(xué)，2011，31（1）：110-116.

［17］暢功民，李建華，張強(qiáng).山西省工礦區(qū)土地復(fù)墾現(xiàn)狀及發(fā)展對(duì)策研究[J].山西農(nóng)業(yè)科學(xué)，2011，39（9）：972-974，980.

［18］盧朝東，李建華，靳東升，等.山西省礦區(qū)土地復(fù)墾與生態(tài)重建現(xiàn)狀及發(fā)展對(duì)策[J].山西農(nóng)業(yè)科學(xué)，2013，41（2）：152-154，182.

［19］胡振琪.山西省煤礦區(qū)土地復(fù)墾與生態(tài)重建的機(jī)遇和挑戰(zhàn)[J].山西農(nóng)業(yè)科學(xué)，2010，38（1）：42-45，64.

［20］劉小翠，白中科，馬蕭，等.露天開采及土地復(fù)墾對(duì)生態(tài)系統(tǒng)服務(wù)功能價(jià)值影響變化的對(duì)比研究：以內(nèi)蒙古勝利東二號(hào)露天煤礦為例[J].資源與產(chǎn)業(yè)，2010（4）：74-78.

［21］謝高地，魯春霞，冷允法，等.青藏高原生態(tài)資產(chǎn)的價(jià)值評(píng)估[J].自然資源學(xué)報(bào)，2003，18（2）：189-196.

Effect of Land Use Change on Ecosystem Services Value in Mining Area Based on GIS

JINHui1，LI Jinling1，HANWeihong1，XINGHongjuan1，GUOJunling2
（1.Institute of Agricultural Resources&Economy，Shanxi Academy of Agricultural Sciences，Taiyuan 030006，China；2.Institute of Agricultural Environment&Resources，Shanxi Academy of Agricultural Sciences，Taiyuan 030031，China）

This paper studied the land use pattern change of Jinfangyuan mining area in Fanshi county of Shanxi province in 2013-2028.It evaluated ecosystem service value of the study area and analyzed the changes of land use pattern and ecosystem service value based on the land use changes.The results showed that land use patterns had remarkable changes between 2013 to 2028.Other grassland,mining land,other woodland and bare land appeared decreasing trend,while forest land,shrub land and artificial grassland showed increasing trend.The land use patterns were changed mainly by other grassland and mining land to artificial grassland,shrub land and woodland.Other grassland and mining land substantially reduced,artificial grassland steadily increased year by year. Meanwhile,ecosystem service value showed decreased trend at first and then increased.Ecosystem service values were mainly controlled by woodland and grassland in the study area.The implementation of mining land reclamation scheme promoted transformation of land use patterns to farmland,woodland and grassland,which ecosystem service value coefficients were higher.Finally,the total ecosystem service value in the mining area improved greatly,which was especially important for the protection of regional ecological security and the sustainable development of ecosystem.

GIS;land use;mining area;ecosystem service;values

F301.24

A

1002-2481（2016）11-1669-06

10.3969/j.issn.1002-2481.2016.11.22

2016-07-08

山西省農(nóng)業(yè)科學(xué)院科技攻關(guān)項(xiàng)目（YGG1430）

金輝（1985-），男，河南濮陽(yáng)人，助理研究員，主要從事基于GIS的農(nóng)業(yè)資源高效利用研究工作。郭軍玲為通信作者。