李志祥，張藝凡，張和生

(1.太原理工大學(xué)礦業(yè)工程學(xué)院，山西 太原 030024；2.山西工程科技職業(yè)大學(xué)，山西 晉中 030619)

煤炭開(kāi)采在促進(jìn)我國(guó)社會(huì)經(jīng)濟(jì)發(fā)展的同時(shí)，也給區(qū)域生態(tài)環(huán)境造成了嚴(yán)重的破壞。據(jù)統(tǒng)計(jì)，我國(guó)井工煤礦每開(kāi)采1萬(wàn)t煤造成的土地塌陷面積為2 400 m2[1]；露天煤礦每開(kāi)采1萬(wàn)t煤損毀土地面積2 200 m2，露天煤礦開(kāi)采年均損毀和占用的土地面積多達(dá)100 km2[2]。煤炭開(kāi)采對(duì)生態(tài)環(huán)境造成的負(fù)效應(yīng)，引起了社會(huì)各界的廣泛關(guān)注，我國(guó)政府先后出臺(tái)了一系列政策法規(guī)，不斷加強(qiáng)煤礦區(qū)土地復(fù)墾工作。

隨著我國(guó)煤礦區(qū)生態(tài)修復(fù)工作的開(kāi)展，如何監(jiān)測(cè)和評(píng)價(jià)礦區(qū)生態(tài)系統(tǒng)成為學(xué)術(shù)界關(guān)注的焦點(diǎn)，學(xué)者們分別從土壤界面[3-4]和植被界面[5-6]研究了礦區(qū)生態(tài)擾動(dòng)特征。近年來(lái)，植被凈初級(jí)生產(chǎn)力(net primary productivity，NPP)研究為探究礦區(qū)生態(tài)擾動(dòng)過(guò)程提供了重要參考。NPP作為生態(tài)系統(tǒng)中的關(guān)鍵特征量，代表植物在單位時(shí)間和單位面積上的有機(jī)物凈積累量，是陸地生態(tài)系統(tǒng)中碳循環(huán)與能量流動(dòng)的重要組成部分[7]。NPP可以直接表示綠色植被在不同生長(zhǎng)條件下的生產(chǎn)能力，同時(shí)能夠反映區(qū)域自然生態(tài)環(huán)境質(zhì)量，是理想的生態(tài)指標(biāo)。侯湖平等[8]基于植被NPP對(duì)徐州礦區(qū)生態(tài)環(huán)境質(zhì)量進(jìn)行了評(píng)價(jià)，明確了植被NPP對(duì)礦區(qū)生態(tài)環(huán)境測(cè)度的可行性。陳廣洲等[9]以淮南煤礦區(qū)為例研究植被NPP的時(shí)空變化，發(fā)現(xiàn)年均NPP呈波動(dòng)緩慢增長(zhǎng)趨勢(shì)，且NPP分布格局與土地利用方式密切相關(guān)。

山西省是我國(guó)重要的煤炭能源基地，長(zhǎng)期的煤炭開(kāi)采給山西省生態(tài)環(huán)境造成強(qiáng)烈擾動(dòng)。據(jù)統(tǒng)計(jì)，截至2014年，山西省因采煤造成的沉陷區(qū)面積達(dá)到3 000 km2[10]。山西省進(jìn)行了廣泛的礦區(qū)生態(tài)修復(fù)活動(dòng)，學(xué)者們也研究了山西省平朔露天煤礦[11]、大同礦區(qū)[12]等礦區(qū)NPP變化特征及其影響因素。上述研究從礦區(qū)尺度揭示了山西省不同采煤方式下的生態(tài)環(huán)境效應(yīng)。然而，對(duì)于如何從省域尺度探究礦區(qū)土地覆被變化下的生態(tài)環(huán)境效應(yīng)的研究目前較少。如何定量監(jiān)測(cè)和評(píng)價(jià)省域尺度煤礦區(qū)土地覆被變化下的生態(tài)環(huán)境響應(yīng)過(guò)程，具有重要的理論和實(shí)踐價(jià)值。因此，本文以山西省六大煤田為研究對(duì)象，利用MOD17A3H數(shù)據(jù)和MOD12Q1數(shù)據(jù)，分析研究區(qū)植被NPP和土地覆被時(shí)空變化，引入NPP損益流，從定向和定量的角度分析省域尺度下礦區(qū)土地覆被變化對(duì)植被NPP的影響，以期為山西省省域礦區(qū)生態(tài)評(píng)價(jià)提供依據(jù)，為省域國(guó)土空間生態(tài)修復(fù)中的礦區(qū)生態(tài)修復(fù)提供參考。

1 研究區(qū)概況

山西省位于我國(guó)黃土高原東翼，是典型的資源型經(jīng)濟(jì)地區(qū)，總面積156 700 km2，緯度34°36′～40°44′N，經(jīng)度110°15′～114°32′E。境內(nèi)地貌類型以山地和丘陵為主，約占全省總面積的80%以上;氣候類型屬溫帶大陸性季風(fēng)氣候，年平均氣溫4～14 ℃，年降水量400～650 mm。 根據(jù)含煤地層的發(fā)育特征和及其構(gòu)造組合，把山西省分為大同煤田、寧武煤田、西山煤田、河?xùn)|煤田、霍西煤田、沁水煤田六大主要煤田[13](圖1)，其面積約為69 655 km2，占全省國(guó)土總面積的44.45%。截至2015年底，山西省共有煤礦1 078座，生產(chǎn)能力14.6×108t/a，累計(jì)生產(chǎn)原煤166×108t，同時(shí)山西省實(shí)施煤炭資源整合和企業(yè)重組，不斷淘汰落后產(chǎn)能，推行“藍(lán)天碧水”及造林綠化等生態(tài)環(huán)境修復(fù)治理工程，2005—2015年累計(jì)完成營(yíng)造林面積約30 000 km2，治理水土流失面積約17 000 km2。

圖1 研究區(qū)位置圖Fig.1 Location of study area

2 數(shù)據(jù)與方法

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

本文所用的NPP數(shù)據(jù)來(lái)自美國(guó)地質(zhì)調(diào)查局提供的MOD17A3H數(shù)據(jù)集，時(shí)間范圍為2005—2015年，空間分辨率為500 m。土地覆被數(shù)據(jù)采用MCD12Q1數(shù)據(jù)集，該數(shù)據(jù)集使用決策樹(shù)分類的方法對(duì)Modis Terra和Aqua反射數(shù)據(jù)進(jìn)行處理，得到分辨率為500 m的全球土地覆被分類結(jié)果[14]。參考通用于土地覆被類型的國(guó)際地圈-生物圈計(jì)劃(IGBP)分類方案，將研究區(qū)土地覆被類型分為耕地、林地、草地、灌叢、建設(shè)用地和其他土地六類。

為保證數(shù)據(jù)的可用性和分析結(jié)果的科學(xué)性，將覆蓋山西省六大煤田的3景MOD17A3H影像和3景MCD12Q1影像分別進(jìn)行鑲嵌，重投影在WGS_1984_UTM_zone_49N坐標(biāo)系，最后輸出空間分辨率為500 m的TIFF圖像。利用六大煤田矢量數(shù)據(jù)對(duì)得到的圖像進(jìn)行裁剪，去除無(wú)效值，得到山西省六大煤田植被NPP數(shù)據(jù)和土地覆被數(shù)據(jù)。

2.2 研究方法

2.2.1 趨勢(shì)分析

利用一元線性回歸分析對(duì)研究區(qū)2005—2015年植被NPP數(shù)據(jù)進(jìn)行逐象元線趨勢(shì)分析，用得到的斜率來(lái)反映其變化趨勢(shì)，計(jì)算公式見(jiàn)式(1)。

(1)

式中：θNPP為研究時(shí)段內(nèi)各象元NPP變化的一元線性回歸方程斜率；n為研究時(shí)段；NPP(i,j)為第i年j象元對(duì)應(yīng)的NPP值。θNPP>0表明植被NPP值呈上升趨勢(shì)；θNPP<0表明植被NPP值呈下降趨勢(shì)。參照李曉榮等[15]的方法對(duì)得到的斜率進(jìn)行顯著性檢驗(yàn)。

2.2.2 NPP損益流

引入動(dòng)態(tài)物質(zhì)變化中“流”的概念[16]，可以更加直觀地了解土地覆被變化所導(dǎo)致的NPP損益。當(dāng)某種土地覆被類型轉(zhuǎn)為另一種類型使得NPP總量增加時(shí)稱為“增益流”，反之稱為“損失流”，該土地覆被類型轉(zhuǎn)為其他土地類型所產(chǎn)生的NPP值變化之和為“NPP損益流”，當(dāng)值為正時(shí)，表示凈增益；值為負(fù)時(shí)，表示凈損失。其表達(dá)式見(jiàn)式(2)。

(2)

式中：NPPi為第i種地類的NPP損益流；NPP(i,j)為地類i轉(zhuǎn)為地類j時(shí)所產(chǎn)生的NPP變化。

3 結(jié)果與分析

3.1 山西省六大煤田2005—2015年NPP空間分布

山西省六大煤田2005—2015年NPP均值空間分布如圖2所示。整體來(lái)看，研究區(qū)NPP分布差異明顯，西北部較低，東南部較高。六大煤田多年NPP均值為335.03 gC/(m2·a)，其中大部分地區(qū)NPP均值分布在200～400 gC/(m2·a)，這些區(qū)域多為耕地和草地的覆蓋區(qū)，人類活動(dòng)相對(duì)頻繁，面積為52 131 km2，占研究區(qū)總面積的74.84%。NPP值高于400 gC/(m2·a)的地區(qū)主要位于沁水煤田東南部、河?xùn)|煤田南部以及寧武煤田中部，面積為14 759 km2，占研究區(qū)總面積的21.19%。 沁水煤田東南部、河?xùn)|煤田南部水熱條件充足，耕地和林地分布廣泛，寧武煤田中部分布有管芩山林區(qū)，屬于黃土高原土石山區(qū)的水源涵養(yǎng)和防護(hù)林區(qū)，植被覆蓋率高，因此NPP值較高。 NPP值低于200 gC/(m2·a)的地區(qū)主要分布于大同煤田北部、河?xùn)|煤田中西部，以及朔州、太原、孝義等礦業(yè)城市周邊，面積為2 765 km2，僅占研究區(qū)總面積的3.97%。大同煤田由于長(zhǎng)期高強(qiáng)度開(kāi)采，侏羅紀(jì)煤層已接近枯竭，加上該煤田所屬區(qū)為半干旱區(qū)，水熱條件差，受礦區(qū)開(kāi)采擾動(dòng)，植被恢復(fù)慢，而河?xùn)|煤田位于呂梁山脈西側(cè)，主要為黃土丘陵區(qū)，該區(qū)水土流失嚴(yán)重，氣候干旱，生態(tài)環(huán)境脆弱，煤炭產(chǎn)量卻逐年升高，因此NPP值偏低。

圖2 2005—2015年NPP均值空間分布圖Fig.2 Spatial distribution map of NPP mean value from 2005 to 2015

對(duì)六大煤田11年的植被NPP值進(jìn)行一元線性回歸，得到其線性回歸變化趨勢(shì)圖(圖3(a))。 研究區(qū)植被NPP值呈現(xiàn)增長(zhǎng)趨勢(shì)的面積為58 280.25 km2，占總面積的83.67%；呈減少趨勢(shì)的面積為11 374.75 km2，占總面積的16.33%。結(jié)合顯著性檢驗(yàn)結(jié)果(圖3(b))可以看出，植被NPP值呈增長(zhǎng)趨勢(shì)較快的地區(qū)主要分布在寧武煤田南部、河?xùn)|煤田中部、霍西煤田東部以及西山煤田東北部，多表現(xiàn)為極顯著增加，表明這些區(qū)域植被狀況明顯好轉(zhuǎn)，生態(tài)環(huán)境得到明顯提高。呈減少趨勢(shì)的地區(qū)主要位于沁水煤田和河?xùn)|煤田南部，其余零星分布于大同、朔州、太原、陽(yáng)泉等各城鎮(zhèn)周邊，顯著性檢驗(yàn)表明，這些地區(qū)變化不顯著，說(shuō)明生態(tài)環(huán)境沒(méi)有發(fā)生明顯惡化，應(yīng)繼續(xù)對(duì)其加強(qiáng)監(jiān)測(cè)，做好礦山生態(tài)修復(fù)工作。

圖3 六大煤田NPP線性回歸趨勢(shì)Fig.3 Linear regression trend of NPP in six coalfields

3.2 山西省六大煤田2005—2015年NPP時(shí)間變化

研究區(qū)NPP年際變化見(jiàn)圖4，植被NPP值總體呈上升趨勢(shì)，其中最大值為2012年的381.88 gC/(m2·a)，最小值為2007年的297.11 gC/(m2·a)。2005—2011年NPP值呈現(xiàn)波動(dòng)增長(zhǎng)的趨勢(shì)，這是這一時(shí)段內(nèi)六大煤田礦山開(kāi)采活動(dòng)和生態(tài)修復(fù)工程共同作用的結(jié)果。礦山開(kāi)采會(huì)使地表下沉，植被遭到破壞，進(jìn)而使植被NPP值降低，而生態(tài)修復(fù)工程可以有效地改善區(qū)域生態(tài)環(huán)境，對(duì)植被NPP產(chǎn)生促進(jìn)作用。2011—2012年NPP均值從324.80 gC/(m2·a)增長(zhǎng)到381.88 gC/(m2·a)，這主要得益于山西省煤炭限產(chǎn)以及生態(tài)修復(fù)工程和退耕還林政策的實(shí)施；2012—2015年NPP均值呈下降的趨勢(shì)，這是由于在山西省正處于轉(zhuǎn)型發(fā)展初期，人類對(duì)土地資源的不合理利用造成的。從NPP比例結(jié)構(gòu)上看，山西省六大煤田2005—2015年NPP值低于100 gC/(m2·a)的區(qū)域均不到1%，同時(shí)NPP值在100～200 gC/(m2·a))的區(qū)域在逐漸減少，而大于400 gC/(m2·a)的區(qū)域呈波動(dòng)增長(zhǎng)趨勢(shì)，這說(shuō)明研究區(qū)植被凈初級(jí)生產(chǎn)力高值區(qū)的比例逐漸增加，生態(tài)環(huán)境逐漸向好。11年來(lái)，六大煤田NPP總量從20 229.26 GgC增長(zhǎng)到23 786.61 GgC，年均增長(zhǎng)1.60%，但同期山西省煤炭總產(chǎn)量從5.5×108t增長(zhǎng)到9.75×108t，年均增長(zhǎng)達(dá)到7.02%。六大煤田的NPP總量隨煤炭產(chǎn)量的快速增長(zhǎng)而增長(zhǎng)，究其原因發(fā)現(xiàn)，研究時(shí)段內(nèi)山西省通過(guò)國(guó)土綠化、退耕還林等生態(tài)修復(fù)工程，森林覆蓋率從14%增長(zhǎng)到20.5%；同時(shí)《山西省礦產(chǎn)資源總體規(guī)劃(2016—2020年)》中指出，2008年以來(lái)山西省累計(jì)投入地質(zhì)環(huán)境修復(fù)治理資金21.53億元，修復(fù)治理礦山362座，治理面積達(dá)到600 km2左右，土地復(fù)墾率達(dá)到50%以上。這表明各項(xiàng)生態(tài)修復(fù)工程對(duì)植被NPP產(chǎn)生的促進(jìn)作用彌補(bǔ)了煤礦開(kāi)采對(duì)NPP的負(fù)面效應(yīng)，山西省六大煤田生態(tài)修復(fù)工作產(chǎn)生了顯著成效。

3.3 山西省六大煤田土地覆被時(shí)空變化

2005—2015年間，研究區(qū)土地覆被類型發(fā)生轉(zhuǎn)移的面積共9 006 km2，占研究區(qū)總面積的12.93%(表1)。耕地、林地、草地之間的相互轉(zhuǎn)化為該地區(qū)土地覆被轉(zhuǎn)移的主要特征。結(jié)合三期土地覆被分類圖(圖5)可以看出，耕地主要集中于六大煤田境內(nèi)的長(zhǎng)治盆地、太原盆地、臨汾盆地和大同盆地，少部分位于河?xùn)|煤田南部的鄉(xiāng)寧縣、吉縣和隰縣等地區(qū)。研究時(shí)段內(nèi)耕地轉(zhuǎn)為其他類型土地的面積共912.75 km2，其中大部分轉(zhuǎn)為草地，轉(zhuǎn)出面積為654.25 km2，占耕地總轉(zhuǎn)出面積的71.68%，這是在加快城市化進(jìn)程的背景下，農(nóng)村進(jìn)城人口增多，大量中低產(chǎn)田被棄耕拋荒，逐漸變?yōu)椴莸氐慕Y(jié)果；轉(zhuǎn)為林地的面積次之，占耕地總轉(zhuǎn)出面積的21.70%，這主要是由于六大煤田持續(xù)開(kāi)展退耕還林生態(tài)修復(fù)工程，大量耕地退耕為林地；在城鎮(zhèn)周邊，隨著城市擴(kuò)張以及工業(yè)發(fā)展的需要，部分耕地被逐步占用， 有57.75 km2的耕地轉(zhuǎn)為建設(shè)用地，占耕地總轉(zhuǎn)出面積的6.30%；由其他類型土地轉(zhuǎn)為耕地的面積共有4 494.5 km2，其中草地為主要的轉(zhuǎn)入源，占耕地總轉(zhuǎn)入面積的71.89%。耕地的轉(zhuǎn)入與轉(zhuǎn)出使其面積凈增長(zhǎng)3 581.75 km2，漲幅達(dá)到17.65%。林地主要分布于沁水煤田的太行山脈和河?xùn)|煤田東部的呂梁山脈一帶。2005—2015年林地轉(zhuǎn)為其他類型土地的面積共2 403 km2，由于管護(hù)缺失、成活率低等原因，大量林地退化轉(zhuǎn)為草地，占林地總轉(zhuǎn)出面積的43%；建設(shè)用地在擴(kuò)張的過(guò)程中也占用了12 km2的林地。 其他類型土地轉(zhuǎn)為林地的面積共2 551.5 km2，草地和耕地為主要轉(zhuǎn)入源，轉(zhuǎn)入面積分別為2 280 km2和198 km2。林地的轉(zhuǎn)入與轉(zhuǎn)出使其面積凈增長(zhǎng)148.65 km2，增長(zhǎng)了1.03%。 草地散布于六大煤田，且一直是研究區(qū)面積占比最大的土地類型。2005—2015年間草地面積凈減少3 882.25 km2，占比從2005年的49.62%降低到2015年的44.05%。建設(shè)用地面積的增加主要來(lái)源于對(duì)城鎮(zhèn)周邊耕地、林地和草地的占用，11年間共增長(zhǎng)84.25 km2，占地增長(zhǎng)了2.67%。灌叢和其他土地的面積由于基數(shù)小，研究時(shí)段內(nèi)增幅分別達(dá)到25.07%和40.13%?？傮w來(lái)說(shuō)，六大煤田土地覆被的主要變化表現(xiàn)為草地面積的減少和其他類型土地面積的增加，說(shuō)明研究區(qū)礦山生態(tài)修復(fù)工作雖然取得一定的成效，但未來(lái)仍應(yīng)繼續(xù)實(shí)施礦山生態(tài)修復(fù)工程，應(yīng)繼續(xù)加強(qiáng)對(duì)礦山生態(tài)修復(fù)工作的重視。

圖5 六大煤田不同時(shí)期土地覆被分類圖Fig.5 Land cover map of six coalfields in different periods

表1 2005—2015年土地利用轉(zhuǎn)移矩陣Table 1 Land use transfer matrix from 2005 to 2015

3.4 不同土地覆被類型對(duì)NPP的影響

將凈初級(jí)生產(chǎn)力數(shù)據(jù)與土地覆被數(shù)據(jù)進(jìn)行疊加分析可知，不同土地覆被類型NPP差異明顯。圖6展示了不同土地覆被類型在研究期內(nèi)的NPP變化趨勢(shì)，所有土地類型的NPP均出現(xiàn)波動(dòng)，但總體都呈現(xiàn)上升趨勢(shì)。 具體來(lái)看，2005—2015年山西省六大煤田不同土地覆被類型NPP均值大小依次為林地(404.87 gC/(m2·a))>灌叢(386.02 gC/(m2·a))>耕地(353.34 gC/(m2·a))>草地(303.64 gC/(m2·a))>其他土地(302.28 gC/(m2·a))>建設(shè)用地(238.69 gC/(m2·a))。這種差異一方面由于林地和灌叢自身植被覆蓋度大，光合作用效率高，另一方面在于研究區(qū)林地和灌叢多分布于山區(qū)，采礦活動(dòng)對(duì)這些植被造成的破壞較小，使之積累的生物量多，同時(shí)土地復(fù)墾和退耕還林等生態(tài)修復(fù)工程使林地面積11年內(nèi)增加148.5 km2，提高了林地的生產(chǎn)力。耕地在農(nóng)業(yè)技術(shù)發(fā)展的背景下，糧食產(chǎn)量明顯提高，一定程度上也促進(jìn)了有機(jī)物的積累。而研究區(qū)草地雖然分布廣泛，但受礦山開(kāi)采的影響，植被覆蓋度較低，加之對(duì)氣候環(huán)境變化敏感，總體退化嚴(yán)重，生產(chǎn)力相對(duì)偏低。

圖6 不同土地覆被類型NPP均值變化趨勢(shì)Fig.6 The change trend of NPP mean value of different land cover types

3.5 土地覆被變化對(duì)植被NPP的影響

將三個(gè)時(shí)段的NPP差值與土地覆被類型在對(duì)應(yīng)年份內(nèi)的轉(zhuǎn)移面積相乘，得到不同土地覆被變化下的NPP損益流(圖7)。研究時(shí)段內(nèi)，沒(méi)有發(fā)生變化的土地類型NPP總量均出現(xiàn)了增加，這是由于不同的地域格局和水熱條件的變化造成了NPP的變化。黃露等[17]在研究鄂爾多斯草地NPP時(shí)也得到了相似的結(jié)果。

圖7 六大煤田NPP損益流Fig.7 Profit and loss flow in six coalfields

2005—2010年，耕地轉(zhuǎn)為其他土地類型時(shí)，NPP總量?jī)魮p失2.94 GgC。 其中耕地轉(zhuǎn)為林地主要發(fā)生在沁水煤田的晉城市周邊，但2009年和2010年晉城市出現(xiàn)的極端氣候較多，抵消了耕地轉(zhuǎn)林地對(duì)NPP的促進(jìn)作用，導(dǎo)致NPP損失流為1.15 GgC。耕地向建設(shè)用地轉(zhuǎn)移時(shí)，轉(zhuǎn)入面積較小，NPP損失流為0.74 GgC，但單位面積損失量達(dá)到34.41 gC/(m2·a)，是單位面積損失量最大的轉(zhuǎn)移方式，這表明耕地向建設(shè)用地轉(zhuǎn)移會(huì)造成NPP的大量損失；林地轉(zhuǎn)為其他土地類型時(shí)，對(duì)NPP損失較大，NPP凈損失71.39 GgC；草地轉(zhuǎn)為耕地和林地時(shí)，NPP增益流為123.88 GgC，其中轉(zhuǎn)為耕地時(shí)增益流最大，占總增益流的92.44%；灌叢在轉(zhuǎn)為林地時(shí)，NPP增益流為1.09 GgC；其他土地轉(zhuǎn)為林地、草地時(shí)，由于植被覆蓋有所改善，NPP出現(xiàn)增益流，但增益較小。

2010—2015年，耕地轉(zhuǎn)為林地和草地時(shí)，NPP增益流為16.31 GgC；耕地轉(zhuǎn)為建設(shè)用地時(shí)，NPP持續(xù)損失，損失流(1.12 GgC)是2005—2010年的1.64倍，說(shuō)明研究區(qū)在這一時(shí)段內(nèi)受城市化進(jìn)程加快的影響，部分耕地被棄耕，逐漸變?yōu)槠渌恋睾凸鄥矔r(shí)，由于缺乏人工管理，植被整體的生產(chǎn)力水平降低，使NPP損失流為0.29 GgC，但損失較前一時(shí)段有所降低。林地轉(zhuǎn)為其他類型土地時(shí)，NPP總量?jī)魮p失33.17 GgC。草地在復(fù)墾為耕地和林地時(shí)，一方面由于復(fù)墾面積大，另一方面是人為促進(jìn)作用加強(qiáng)，有效地保證了植被的生長(zhǎng)，使NPP增益流達(dá)到111.37 GgC；草地轉(zhuǎn)為灌叢時(shí)，NPP損失流為0.35 GgC；灌叢和其他土地轉(zhuǎn)為林地和草地時(shí)，NPP總增益流為5.6 GgC。

2005—2015年間，不同土地類型轉(zhuǎn)為耕地和林地，使NPP總增益流(266.50 GgC)最大；林地轉(zhuǎn)為其他土地類型，使NPP總損失流(185.57 GgC)最大?？傮w來(lái)說(shuō)，不同土地覆被之間的相互轉(zhuǎn)換，使研究區(qū)2005—2015年間NPP凈增益216.35 GgC，說(shuō)明2005—2015年間山西省六大煤田生態(tài)環(huán)境逐漸好轉(zhuǎn)。

4 結(jié) 論

1) 山西省六大煤田2005—2015年植被NPP變化呈現(xiàn)波動(dòng)式漸進(jìn)增長(zhǎng)，NPP總量增加3 492.20 GgC，年均增速達(dá)到1.6%。空間分布表現(xiàn)出東高西低，南高北低的特征，NPP呈增加趨勢(shì)的區(qū)域占整個(gè)研究區(qū)的83.67%，高值區(qū)域面積在逐漸增加，而低值區(qū)域面積逐漸減少，說(shuō)明研究區(qū)的生態(tài)環(huán)境整體持續(xù)好轉(zhuǎn)。

2) 研究時(shí)段內(nèi)，土地覆被變化主要表現(xiàn)為耕地和林地面積持續(xù)增加，耕地增加了3 581.75 km2，林地增加了148.65 km2，草地面積持續(xù)減少，減少了3 882.25 km2，建設(shè)用地面積增加了84.25 km2，其他類型土地變化相對(duì)較小。

3) 2005年以來(lái)，不同土地覆被類型的NPP均值大小為林地>灌叢>耕地>草地>其他土地>建設(shè)用地。土地覆被類型未發(fā)生變化的區(qū)域植被凈初級(jí)生產(chǎn)力穩(wěn)步增長(zhǎng)，表明生態(tài)環(huán)境逐漸向好。增加林地和耕地的面積會(huì)使NPP總量產(chǎn)生明顯增益，但林地轉(zhuǎn)為耕地和草地以及耕地轉(zhuǎn)為建設(shè)用地會(huì)造成NPP的大量損失。

山西省生態(tài)修復(fù)工作在六大煤田取得了一定的成效，未來(lái)仍應(yīng)繼續(xù)實(shí)施礦區(qū)生態(tài)修復(fù)工程，具體措施如下所述。

1) 根據(jù)各礦區(qū)水土要素差異，因地制宜開(kāi)展生態(tài)修復(fù)。晉東南地區(qū)要在加強(qiáng)生態(tài)環(huán)境監(jiān)測(cè)工作的基礎(chǔ)上，繼續(xù)擴(kuò)大營(yíng)林造林面積；西北部礦區(qū)應(yīng)增加林地的生態(tài)用地比重，提高區(qū)域生態(tài)環(huán)境質(zhì)量。

2) 在改進(jìn)復(fù)墾技術(shù)的基礎(chǔ)上，充分利用生態(tài)系統(tǒng)的自然恢復(fù)能力，將人工修復(fù)和自然恢復(fù)相結(jié)合，提高土地復(fù)墾的效率。

3) 加強(qiáng)礦區(qū)生態(tài)用地管護(hù)，防止工業(yè)用地對(duì)生態(tài)用地的過(guò)度占用。