張艷珍,王釗齊,楊 悅,李建龍,張 穎,章超斌

(南京大學(xué)生命科學(xué)學(xué)院生態(tài)學(xué)系，江蘇 南京 210023)

草地是陸地生態(tài)系統(tǒng)的組成成分,在全球碳循環(huán)中扮演重要角色,不僅有保持水土、防風(fēng)固沙等生態(tài)功能,更是草原區(qū)畜牧業(yè)生產(chǎn)發(fā)展、經(jīng)濟(jì)進(jìn)步的重要環(huán)節(jié)[1-3]。蒙古高原作為東北亞地區(qū)的一個(gè)特殊地貌單元,遠(yuǎn)離海洋,草地是其主要植被覆蓋類型,荒漠化是其重要的地貌特征,是全球變化研究中的必要環(huán)節(jié)[4-7]。近40年來，由氣候變化和人類活動(dòng)引起的蒙古高原草地退化形式加劇,裸地面積有擴(kuò)大趨勢(shì),荒漠化十分嚴(yán)重[8-11]。草地退化一方面直接影響內(nèi)外蒙古廣大牧民群眾的生活水平和社會(huì)經(jīng)濟(jì)的健康發(fā)展;另一方面,作為我國(guó)北方重要的生態(tài)屏障,草地退化會(huì)加重高原荒漠化、沙化和水土流失,造成嚴(yán)重的生態(tài)環(huán)境問題。

草地植被覆蓋變化可表征草地植被長(zhǎng)勢(shì)、生態(tài)系統(tǒng)健康程度以及植被退化狀況，獲取草地植被覆蓋度變化特征在區(qū)域生態(tài)環(huán)境評(píng)價(jià)中有重要意義[12]。遙感監(jiān)測(cè)手段隨著計(jì)算機(jī)和遙感等技術(shù)的發(fā)展已經(jīng)成為大尺度區(qū)域乃至全球范圍植被覆蓋變化監(jiān)測(cè)研究的主要手段[10]。歸一化植被指數(shù)(normalized difference vegetation index,NDVI)由于與植被覆蓋度存在極顯著線性相關(guān)關(guān)系[13],是可以反映植被生長(zhǎng)分布及其年際變化且可用來進(jìn)行植被變化監(jiān)測(cè)的一種遙感數(shù)據(jù)[14]。在現(xiàn)有的NDVI數(shù)據(jù)中,EOS/MODIS NDVI 的光譜性能穩(wěn)定,輻射定標(biāo)精度優(yōu)越,空間相關(guān)性高,能較好地估算出高原地區(qū)的草地生長(zhǎng)狀況[15-18]。

草地退化是一種人為因素和氣候因素等對(duì)草地生態(tài)系統(tǒng)進(jìn)行干擾超過自我調(diào)節(jié)閾值,使其出現(xiàn)逆向演替的現(xiàn)象[19]。目前，關(guān)于草地退化監(jiān)測(cè)和評(píng)價(jià)的方法主要有實(shí)測(cè)法和遙感監(jiān)測(cè)法等。實(shí)測(cè)法適用于小范圍內(nèi)的植被退化監(jiān)測(cè),而針對(duì)大尺度區(qū)域,遙感技術(shù)則成為了重要手段。近年來,有不少針對(duì)我國(guó)重要牧區(qū)的草地退化狀況研究。薛存芳等[20]利用實(shí)測(cè)法選擇了種群自然生殖枝高度、草地植被覆蓋度和草地生物量3個(gè)評(píng)價(jià)指標(biāo),加權(quán)得到草地退化指數(shù)GDI,結(jié)果顯示,與2002年相比，2006年中國(guó)內(nèi)蒙古地區(qū)草地植被退化狀況呈現(xiàn)出整體改善、局部惡化的狀況。曹旭娟等[21]利用遙感監(jiān)測(cè)法基于NDVI反演藏北地區(qū)草地植被覆蓋度,以覆蓋度作為判斷草地退化的標(biāo)準(zhǔn),綜合計(jì)算逐年草地退化指數(shù),得出該地區(qū)草地退化程度分布結(jié)果。由于蒙古高原覆蓋面積較大,涉及到兩個(gè)行政單元,本研究采用遙感監(jiān)測(cè)法,以草地植被覆蓋度為指標(biāo)監(jiān)測(cè)草地退化,建立蒙古高原草地退化評(píng)價(jià)指標(biāo),參照國(guó)家標(biāo)準(zhǔn)(GB19377-2003)[22],以1982-1985年植被狀況為基準(zhǔn),研究蒙古高原1986-1999年和2000-2013年(2000年以來，我國(guó)陸續(xù)實(shí)施了“退牧還草”、“京津源風(fēng)沙治理”等工程,蒙古國(guó)則實(shí)施草地管理及土地法等舉措，為判斷生態(tài)治理工程的實(shí)際意義選擇了2000年為前后兩段時(shí)期的節(jié)點(diǎn))的植被退化時(shí)空分布特征,通過精確定量監(jiān)測(cè)和分析蒙古高原草地退化程度及其時(shí)空分布,積累蒙古高原草地退化基本數(shù)據(jù),有利于高原草地資源的可持續(xù)發(fā)展和生態(tài)環(huán)境的改善。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)是蒙古高原主體部分(圖1)，包括蒙古國(guó)和中國(guó)內(nèi)蒙古自治區(qū)兩個(gè)行政單元33個(gè)行政區(qū)(蒙古國(guó)21個(gè)省市，中國(guó)內(nèi)蒙古12個(gè)盟市),位于87°43′-126°04′ E和37°22′-53°20′ N,面積約為2.74×106km2,其中蒙古國(guó)和中國(guó)內(nèi)蒙古自治區(qū)的面積分別為1.56×106和1.18×106km2。蒙古高原位于歐亞大陸中部,東抵大興安嶺西麓,西達(dá)阿爾泰山脈和薩彥嶺,南部以陰山山脈為界，北臨薩彥嶺、肯特山脈和亞布洛夫山脈[23]。平均海拔約1 580 m,地勢(shì)西高東低，年降水量約269 mm,各地分布不均,東北部地區(qū)最多,年總量300～400 mm,降水量隨著離海洋距離的增加由東北向西南逐漸減小,西南部地區(qū)降水量?jī)H約100 mm。蒙古高原年均溫是4.12 ℃,分布特征基本與降水量一致。高原自然植被覆蓋受氣候因子影響,由東向西南依次有森林、森林草原、草甸草原、典型草原、荒漠草原和戈壁荒漠,生態(tài)環(huán)境多樣且脆弱[24]。

圖1 研究區(qū)位置圖及國(guó)際生物圈計(jì)劃(IGBP)植被類型Fig. 1 Location of the study area and the land use map from International Geosphere-Biosphere Programme(IGBP)

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

1982-2006年GIMMS NDVI數(shù)據(jù)來自美國(guó)國(guó)家航天航空局(NASA),時(shí)間分辨率為15 d,空間分辨率為8 km,該數(shù)據(jù)集制作過程中已經(jīng)經(jīng)過輻射校正、圖像幾何精校正、除云和除壞線等處理,數(shù)據(jù)誤差小且精度較高,已被廣泛應(yīng)用于大尺度區(qū)域植被覆蓋變化研究中;2001-2013年EOS/MODIS NDVI數(shù)據(jù)來自NASA地球觀測(cè)系統(tǒng)EOS/MODIS NDVI數(shù)據(jù)中的MOD13A2級(jí)數(shù)據(jù)產(chǎn)品,時(shí)間分辨率為16 d,空間分辨率1 km。土地覆蓋分類數(shù)據(jù)采用MODIS數(shù)據(jù)MOD12C1產(chǎn)品中的IGBP 2001年數(shù)據(jù),通過投影變換、重采樣、裁剪和類型合并等預(yù)處理,生成空間分辨率1 km的,包括林地、灌叢、草地、農(nóng)耕地和稀疏植被5類,并提取草地覆蓋類型用于分析草地植被覆蓋度面積動(dòng)態(tài)和草地退化程度分析。

1.3 研究方法

1.3.1植被覆蓋度(vegetation fraction,Fv)的計(jì)算 植被覆蓋度是反映草地資源是否發(fā)生退化的直觀指標(biāo),由于植被蓋度與歸一化植被指數(shù)(NDVI)存在極顯著線性相關(guān)關(guān)系[18],本研究利用像元二分模型近似估算植被覆蓋度,具體計(jì)算公式如下:

(1)

式(1)中:Fv表示草地的覆蓋度,NDVImin代表研究范圍內(nèi)的NDVI最小值(或裸地NDVI值),NDVImax代表研究范圍內(nèi)的NDVI最大值(或純植被NDVI值)。

1.3.2草地覆蓋度的變化趨勢(shì)計(jì)算 本研究利用線性回歸方程在像元基礎(chǔ)上分析蒙古高原1986-1999年和2000-2013年最大草地覆蓋度的變化趨勢(shì),具體的計(jì)算公式如下:

(2)

式(2)中:S代表趨勢(shì)斜率,n為研究年限14,i=1為1986或2000年,i=2為1999或2013年,Ci為第i年的最大覆蓋度值,S為正值則草地覆蓋度呈上升趨勢(shì),若為負(fù)值則草地覆蓋度呈下降趨勢(shì)。關(guān)于變化趨勢(shì)顯著性分析則使用F檢驗(yàn),顯著性檢驗(yàn)結(jié)果僅代表變化趨勢(shì)可置信程度的高低,而與其變化的快慢沒有關(guān)系[25]。顯著性分析計(jì)算公式為:

(3)

(4)

(5)

1.3.3草地退化評(píng)價(jià)指標(biāo) 建立基于遙感監(jiān)測(cè)草地退化的評(píng)價(jià)指標(biāo)選擇原則既要考慮遙感原理和數(shù)據(jù)特點(diǎn),也要有科學(xué)性、可行性和針對(duì)性。由于尚未有針對(duì)蒙古高原整體或是蒙古國(guó)內(nèi)的草地退化評(píng)價(jià),而內(nèi)、外蒙古植被類型和環(huán)境條件具有區(qū)域延伸性,參照我國(guó)國(guó)家標(biāo)準(zhǔn)(GB 19377-2003)，將20世紀(jì)80年代初研究區(qū)域未退化草地植被特征定為未退化草地標(biāo)準(zhǔn)[20，26]。本研究以蒙古高原1982-1985年間的最大草地植被覆蓋度為未退化草地基準(zhǔn),將草地退化程度分為5個(gè)級(jí)別,覆蓋度高于未退化草地覆蓋度的90%的退化程度評(píng)分為1分,覆蓋度為未退化草地覆蓋度的75%～90%評(píng)分為2分,覆蓋度為未退化草地覆蓋度的60%～75%評(píng)分為3分,覆蓋度為未退化草地覆蓋度的30%～60%評(píng)分為4分,覆蓋度小于未退化草地覆蓋度的30%為評(píng)分為5分。在此基礎(chǔ)上,選用草地退化指數(shù)(grassland degradation index，GDI)判斷并監(jiān)測(cè)蒙古高原的草地退化情況[26]。GDI的具體計(jì)算公式如下:

(6)

式(6)中:GDI為草地退化指數(shù);Di為草地退化等級(jí)i的評(píng)分;Ai為草地退化等級(jí)i的分布面積;A為研究區(qū)的草地植被總面積。根據(jù)計(jì)算得出的GDI值,參照蒙古高原草地退化總和判別指標(biāo)(表1),得出草地退化等級(jí)結(jié)果。

1.4 數(shù)據(jù)分析

采用MODIS Reprojection Tools(MRT)軟件對(duì)下載的數(shù)據(jù)進(jìn)行格式轉(zhuǎn)換、數(shù)據(jù)拼接和裁剪,采用最大值合成法消除部分大氣、云和太陽高度角的干擾,為進(jìn)一步消除干擾,采用三點(diǎn)平滑方法對(duì)兩種NDVI數(shù)據(jù)進(jìn)行進(jìn)一步修正,得到研究區(qū)逐年NDVI數(shù)據(jù)。植被覆蓋度的反演以及草地退化程度級(jí)別等利用ENVI 5.1中的Band Math工具計(jì)算;遙感制圖在ArcGIS 10.2中進(jìn)行;統(tǒng)計(jì)數(shù)據(jù)的統(tǒng)計(jì)分析和圖表制作利用Excel 2010和Sigma Plot 12 .5實(shí)現(xiàn)。

表1 蒙古高原草地退化綜合判別指標(biāo)Table 1 Integrated evaluation indices of grassland degradation on the Mongolia Plateau

2 結(jié)果與分析

2.1 蒙古高原草地植被覆蓋度時(shí)空分布特征和變化趨勢(shì)分析

蒙古高原1986-1999年和2000-2013年草地植被覆蓋度空間分布整體呈現(xiàn)由東北向西南逐漸減小的趨勢(shì)(圖2)。1986-1999年和2000-2013年平均覆蓋度分別為15.60%和18.43%,總體呈增加趨勢(shì)。其中，1986-1999年草地覆蓋度顯著和極顯著增加的區(qū)域主要分布在蒙中國(guó)杭愛山南麓、中國(guó)陰山山脈地區(qū)、中國(guó)赤峰市中南部和通遼市南部少數(shù)區(qū)域;2000-2013年草地覆蓋度顯著(P<0.05)和極顯著(P<0.01)增加的區(qū)域主要分布在蒙古國(guó)杭愛山脈東部、中國(guó)肯特山脈區(qū)域、中國(guó)東方省東部地區(qū)、中國(guó)興安盟東部和中國(guó)鄂爾多斯南部和東部一帶。從變化率來看，2000年前變化率較大區(qū)域主要集中分布于中國(guó)呼和浩特及其鄰近區(qū)域和通遼市、赤峰市南部地帶，2000年后變化率較大區(qū)域則主要分布在蒙古國(guó)東方省和蘇赫巴托爾省地區(qū)。在年際變化過程中(表2),1986-1999年間,草地覆蓋度面積最大的是30%～60%,在中國(guó)內(nèi)蒙古、蒙古國(guó)和蒙古高原分別達(dá)到1 225.74、2 193.53和3 415.42 km2;面積最小的是90%～100%,中國(guó)內(nèi)蒙古、蒙古國(guó)和蒙古高原面積分別為14.24、1.36和15.60 km2。2000-2013年間,草地覆蓋度面積較大的是30%～60%和60%～75%,蒙古高原分別達(dá)到2 487.42和2 237.61 km2;其次依次是75%～90%、0～30%和90%～100%,面積分別為1 357.53、455.07和82.06 km2。整體來看,草地植被覆蓋度0～30%和30%～60%呈下降趨勢(shì),60%以上覆蓋度的面積均呈上升趨勢(shì)。從變化速率顯著性分析的面積統(tǒng)計(jì)可以看出(表3),兩個(gè)研究時(shí)段草地覆蓋度均呈現(xiàn)上升趨勢(shì),1986-1999年無顯著(P>0.05)增加、顯著(P<0.05)增加和極顯著(P<0.01)增加的草地面積分別為3 706.37、298.86和189.67 km2;2000-2013年無顯著增加、顯著增加和極顯著增加的面積分別為4 526.06、443.32和92.46 km2。

圖2 蒙古高原草地覆蓋度時(shí)空分布、變化趨勢(shì)和顯著性分析分析Fig. 2 The space distribution features and change trend analysis of grassland vegetation coverage in the Mongolia Plateau

ESD：極顯著減少(P<0.01)；SD：顯著減少(P<0.05)；NSD：減少不顯著(P>0.05)；NSI：增加不顯著(P>0.05)；SI：顯著增加(P<0.05)；ESI：極顯著增加(P<0.01)。

ESD: Extreme significantly decrease(P<0.01); SD: Significantly decrease(P<0.05); NSD: No significantly decrease(P>0.05); NSI: No significantly increase(P>0.05); SI: Significantly increase(P<0.05); ESI: Extreme significantly increase(P<0.01).

表2 蒙古高原不同草地植被覆蓋度的面積變化Table 2 Change of area statistics of grassland coverage in the Mongolia Plateau

表3 蒙古高原草地植被覆蓋度變化率顯著性面積統(tǒng)計(jì)Table 3 Area statistics of grassland coverage about different significance levels in the Mongolia Plateau

2.2 蒙古高原草地植被覆蓋度差值分析

研究時(shí)段為1986-2013年,為達(dá)到蒙古高原草地覆蓋在不同時(shí)段變化情況直觀判斷的目的,選擇1986年作為基準(zhǔn)比較年,將1995、2001、2007和2013年的草地覆蓋度分別與1986年的草地覆蓋度進(jìn)行差值計(jì)算(圖3)。然后，根據(jù)實(shí)際情況將草地覆蓋度變動(dòng)幅度進(jìn)行區(qū)間面積統(tǒng)計(jì)(表4)。相對(duì)于1986年,1995、2001和2007年草地覆蓋度依次變差,而2013年草地覆蓋度最好,且高于1995基準(zhǔn)比較年差值,表現(xiàn)為初期一般、中期變差、后期變好且好于初期的情況(圖3)?？臻g分布來看,變化最為劇烈的草地區(qū)域主要集古中在蒙國(guó)肯特山脈附近、東方省、扎布汗省和蘇赫巴托爾省東部,中國(guó)內(nèi)蒙古呼倫貝爾西部、興安盟東部扎賚特旗附近區(qū)域和錫林浩特中東部地區(qū);相對(duì)穩(wěn)定的草地區(qū)域集中于蒙古國(guó)高原北部蒙古國(guó)薩彥嶺、南部陰山山脈區(qū)域以及中國(guó)赤峰市和通遼市南部。

圖3 1986年與1999、2001、2007、2013年蒙古高原草地植被覆蓋度差值圖Fig. 3 The differential chart of grassland vegetation coverage in the Mongolia Plateau between 1986 and 1999, 2001, 2007 and 2013

從覆蓋度變動(dòng)幅度來看(表4),相對(duì)1986年,草地覆蓋度面積負(fù)變化區(qū)域主要集中于[-0.1,-0.001),少數(shù)集中在[-0.2,-0.1),特定時(shí)期少數(shù)區(qū)域集中表現(xiàn)為[-0.3,-0.2);草地覆蓋度面積正變化區(qū)域主要集中于[0.001,0.1)和[0.1,0.2)。其中,2001年較1986年草地覆蓋度負(fù)變化的面積最大,達(dá)到16 730.99 km2,其次依次是2007年的15 843.23 km2、1995年的13 684.29 km2和2013年的13 160.72 km2;2013年較1986年草地覆蓋度正變化面積最大,達(dá)到8 496.5 km2,其次依次為1995年、2007年和2001年,面積分別為7 972.93、5 813.99和4 926.23 km2。

2.3 蒙古高原草地退化狀況時(shí)空分布特征

通過計(jì)算草地退化指數(shù)GDI,得出1986-1999年和2000-2013年蒙古高原草地退化程度分布特征(圖4)。從空間分布來看,退化程度區(qū)域性明顯,兩階段草地植被退化均呈現(xiàn)東北向西南依次逐漸嚴(yán)重趨勢(shì)。未退化和輕度退化區(qū)域主要集中于高原北部蒙古國(guó)薩彥嶺區(qū)域、中國(guó)大興安嶺區(qū)域;高原中部區(qū)域主要表現(xiàn)出中度和重度退化,極重度退化出現(xiàn)在草地與非草地過渡區(qū)域。

1986-1999年重度退化草地面積所占百分比最大,為50.35%,其次依次為極重度退化、中度退化、輕度退化和未退化,面積比例分別是20.37%、18.88%、9.50%和0.90%(圖5);2000-2013年中度和重度退化草地面積百分比超過了30%,分別是39.42%和30.34%,未退化草地面積最小,僅占1.95%。結(jié)果表明，草地退化狀況由極重度、重度向中度和輕度變化,草地退化狀況轉(zhuǎn)好。

3 討論和結(jié)論

本研究以覆蓋度作為表征草地退化的主要指標(biāo),從2000年前后14年蒙古高原草地年際變化可以看出,草地有明顯的恢復(fù)趨勢(shì)，中國(guó)內(nèi)蒙古、蒙古國(guó)和蒙古高原覆蓋度上升斜率分別為0.23、0.28和0.26。中國(guó)內(nèi)蒙古大興安嶺地區(qū)、蒙古國(guó)肯特山脈及薩彥嶺東部山區(qū)主要植被類型是森林，植被覆蓋度最高,均高于60%;南部和西南部與荒漠區(qū)過渡帶草地覆蓋度最低,小于30%。近年來，針對(duì)蒙古高原植被覆蓋狀況的研究已有不少(表5),張宏斌等[27]利用NDVI分析中國(guó)內(nèi)蒙古植被的研究結(jié)果表明，2000-2008年草原植被總體趨于恢復(fù)。包剛等[28]基于GIMMS NDVI 分析蒙古高原植被覆蓋狀況認(rèn)為，1982-2006年高原植被年平均上升趨勢(shì)為0.000 4，蒙古國(guó)的年均上升趨勢(shì)(0.000 5)比中國(guó)內(nèi)蒙古的年均上升趨勢(shì)(0.000 3)更加顯著。周錫飲等[29]研究蒙古高原地區(qū)1981-2006年植被覆蓋的時(shí)空變化顯示,草原區(qū)呈現(xiàn)上升趨勢(shì),森林區(qū)和荒漠區(qū)植被有下降趨勢(shì)。本研究1986-1999年草地植被覆蓋度顯著和極顯著增加的高原面積分別為2 98.86和189.67 km2;2000-2013年分別為443.32和92.46 km2，增加面積也呈現(xiàn)出上升趨勢(shì)。

表4 蒙古高原草地植被覆蓋度正負(fù)變化面積統(tǒng)計(jì)Table 4 Positive and negative variation statistics of grassland vegetation coverage area in the Mongolia Plateau km2

圖4 1986-1999年和2000-2013年蒙古高原草地退化分布圖Fig. 4 Grassland degradation on the Mongolia Plateau during the periods 1986-1999 and 2000-2013

圖5 1986-1999和2000-2013年蒙古高原草地退化程度面積百分比統(tǒng)計(jì)圖Fig. 5 Statistics of grassland degradation degree areas on the Mongolia Plateau during the periods 1987-1999 and 2000-2013

研究區(qū) Studyarea 研究時(shí)段Timescale覆蓋度Fv歸一化植被指數(shù)NDVI凈初級(jí)生產(chǎn)力NPP植被狀況描述 Vegetationdescription 參考文獻(xiàn)Reference蒙古高原MongoliaPlateau2002-20120．32-0．52先退化后改善趨勢(shì)Atreadtoimproveafterdegradation[30]蒙古高原MongoliaPlateau1982-20060．26-0．30植被整體呈改善趨勢(shì)Atrendofimprovement[10]蒙古高原MongoliaPlateau1982-20060．2038高原植被覆蓋呈上升趨勢(shì)Atrendofimprovement[28]中國(guó)內(nèi)蒙古InnerMongolia，China2000-2010草原區(qū)Grassland0．4Fv整體上呈現(xiàn)上升趨勢(shì)Fvshowedanoverallascendanttrend[31]中國(guó)內(nèi)蒙古InnerMongolia，China2001-2010281．3g·(m2·a)-1NPP呈上升趨勢(shì)NPPshowedanoverallascendanttrend[32]中國(guó)鄂爾多斯Ordos，China2001-20130．14-0．20生長(zhǎng)季呈增加趨勢(shì)Anincreasetrendinthegrowingseason[33]中國(guó)錫林郭勒草原XilinGolgrassland，China2000-20100～0．5kg·m-2波動(dòng)變化，呈現(xiàn)增加趨勢(shì)Fluctuationchange，showedatrendofincrease[34]中國(guó)渾善達(dá)克沙地OtindagSandyLand，China2000-20140．2-0．4生長(zhǎng)季平均值呈上升趨勢(shì)Thelengthofthegrowingseasonisincreasing[35]

從草地退化程度的空間分布特征來看,草地退化區(qū)域性較為明顯,極重度退化區(qū)域主要集中在高原西南部草地植被和荒漠區(qū)的過渡地帶,調(diào)查發(fā)現(xiàn),該區(qū)有大量礦產(chǎn)資源,是主要礦區(qū);草地退化程度由西南向東北方向逐漸減弱,除人為因素外,另一方面是由于蒙古高原典型的干旱半干旱氣候類型,水熱分布不均,地帶性明顯。高原東部和北部分別受太平洋和北冰洋水汽影響,年降水量最多,隨著離兩個(gè)海洋距離的增加,降水量逐漸減少[36]。蒙古高原草地退化程度區(qū)域性明顯,1986-1999年和2000-2013年兩階段草地植被退化程度均呈現(xiàn)東北向西南依次逐漸嚴(yán)重趨勢(shì)。未退化和輕度退化區(qū)域主要集中于高原北部蒙古國(guó)薩彥嶺區(qū)域、中國(guó)大興安嶺區(qū)域;高原中部區(qū)域主要表現(xiàn)出中度和重度退化,極重度退化出現(xiàn)在草地與非草地過渡區(qū)域。草地植被退化程度明顯改善與21世紀(jì)以來我國(guó)實(shí)施退耕還林還草、退牧還草和京津風(fēng)沙源治理、蒙古國(guó)實(shí)施土地法等一系列環(huán)境治理工程的實(shí)施息息相關(guān)[31]。本研究結(jié)果表明，草地退化狀況有由極重度、重度向中度和輕度變化的趨勢(shì),退化狀況轉(zhuǎn)好,在蒙古高原草地地帶性植被條件下,利用草地退化指數(shù)方法可以評(píng)估草地退化基況和程度，不過，鑒于評(píng)價(jià)因子相對(duì)單一，今后研究中應(yīng)該考慮多因子綜合評(píng)價(jià)，從土壤、植被和生產(chǎn)力等多方面著手分析。

References:

[1] Ojima D S,Dirks B O,Glenn E P,Owensby C E,Scurlock J O.Assessment of budget for grasslands and drylands of the world.Water,Air,and Soil Pollution,1993,70(1/4):95-109.

[2] White R P,Murry S,Rohweder M.Grassland Ecosystems.Washington DC,USA:World Resources Institute,2000.

[3] 劉桂香.草原環(huán)境質(zhì)量檢測(cè)評(píng)價(jià)現(xiàn)狀、問題及對(duì)策.中國(guó)草地學(xué)報(bào),2009,31(3):8-12.

Liu G X.Current situation of monitoring,assessment,problems and countermeasures of grassland eco-environmental quality.Chinese Journal of Grassland,2009,31(3):8-12.(in Chinese)

[4] 王銀柱,戴良先,劉曉英,劉剛,武高林.我國(guó)草原現(xiàn)狀及其生態(tài)恢復(fù)途徑初探.草業(yè)與畜牧,2007(5):29-32.

Wang Y Z,Dai L X,Liu X Y,Liu G,Wu G L.Grassland status and ecological restoration approaches in China.Caoye Yu Xumu,2007(5):29-32.(in Chinese)

[5] 張凱,高會(huì)旺,張仁建,朱燕君,王躍思.我國(guó)沙塵的來源、移動(dòng)路徑及對(duì)東部海域的影響.地球科學(xué)進(jìn)展,2005,20(6):627-636.

Zhang K,Gao H W,Zhang R J,Zhu Y J,Wang Y S.Sources and movement routes of sand-dust aerosols and their impact probabilities on China seas in 2000-2002.Advances in Earth Science,2005,20(6):627-636.(in Chinese)

[6] 孫艷玲,郭鵬,延曉東,趙天寶.內(nèi)蒙古植被覆蓋變化及其與氣候、人類活動(dòng)的關(guān)系.自然資源學(xué)報(bào),2010,25(3):407-414.

Sun Y L,Guo P,Yan X D,Zhao T B.Dynamics of vegetation cover and its relationship with climate change and human activities in Inner Mongolia.Journal of Natural Resources,2010,25(3):407-414.(in Chinese)

[7] 師華定,周錫飲,孟凡浩,白鶴鳴.30年來蒙古國(guó)和內(nèi)蒙古的LUCC區(qū)域分異.地球信息科學(xué)學(xué)報(bào),2013,15(5):719-725.

Shi H D,Zhou X Y,Meng F H,Bai H M.Mongolia and Inner Mongolia LUCC regional differentiation over the past 30 years.Journal of Geo-information Science,2013,15(5):719-725.(in Chinese)

[8] Wang L,Zhen L,Liu X L,Batkhishig O,Wang Q X.Comparative studies on climate changes and influencing factors in central Mongolia Plateau Region.Geographical Research,2008,27(1):171-180.

[9] 馬梅,張圣微,魏寶成.錫林郭勒草原近30年草地退化的變化特征及其驅(qū)動(dòng)因素分析.中國(guó)草地學(xué)報(bào),2017,39(4):86-93.

Ma M,Zhang S H,Wei B C.Temporal and spatial pattern of grassland degradation and its determinants for recent 30 years in Xilingol.Chinese Journal of Grassland,2017,39(4):86-93.(in Chinese)

[10] 戴琳,張麗,王昆,王仁禮.蒙古高原植被變化趨勢(shì)及其影響因素.水土保持通報(bào),2014,34(5):218-225.

Dai L,Zhang L,Wang K,Wang R L.Vegetation changing trend and its affecting factors in Mongolian Plateau.Bulletin of Soil and Water Conservation,2014,34(5):218-225.(in Chinese)

[11] 馬崇勇，張卓然，杜桂林，單艷敏，季彥華.內(nèi)蒙古草原蝗害與綠色防控技術(shù)應(yīng)用現(xiàn)狀.草業(yè)科學(xué),2017,34(5):1113-1120.

Ma C Y,Zhang Z R,Du G L,Shan Y M,Ji Y H.Locust plague and its present situation after the application of green control technology in grasslands of Inner Mongolia.Pratacultural Science,2017,34(5):1113-1120.(in Chinese)

[12] 毛飛,盧志光,張佳華,侯英雨.近20年藏北地區(qū)AVHRR NDVI與氣候因子的關(guān)系.生態(tài)學(xué)報(bào),2007,27(8):3198-3205.

Mao G,Lu Z G,Zhang J H,Hou Y Y.Relations between AVHRR NDVI and climate factors in Northern Tibet in recent 20 years.Acta Ecologica Sinica,2007,27(8):3198-3205.(in Chinese)

[13] 王茜，陳瑩，阮璽睿，李小冬，王小利，史健宗，王志偉.1982-2012年中國(guó)NDVI變化及其與氣候因子的關(guān)系.草地學(xué)報(bào),2017,25(4):691-700.

Wang Q，Chen Y，Ruan X R，Li X D，Wang X L，Shi J Z，Wang Z W.The changes of NDVI in China from 1982 to 2012 and its relationship with climatic factors.Acta Agrestia Sinica,2017,25(4):691-700.(in Chinese)

[14] 邢文淵,李大平,王蕾,石玉,肖然,肖繼東.巴里坤草原面積動(dòng)態(tài)及其氣象因素分析.草業(yè)科學(xué),2014,31(3):408-414.

Xing W Y,Li D P,Wang L,Shi Y,Xiao R,Xiao J D.Dynamic variation of grassland area and meteorological influence analysis in Barkol.Pratacultural Science,2014,31(3):408-414.(in Chinese)

[15] Schulz C,Koch R,Cierjacks A,Kleinschmit B.Land change and loss of landscape diversity at the Caatingaphytogeographical domain:Analysis of pattern-process relationships with MODIS land cover products(2001-2012).Journal of Arid Environments,2017,136:54-74.

[16] Zhou W,Gang C,Zhou L,Chen Y Z,Li J L,Ju W M,Odeh I.Dynamic of grassland vegetation degradation and its quantitative assessment in the northwest China.Acta Oecologica,2014,55(2):86-96.

[17] Zhou W,Gang C,Chen Y,Mu S J,Sun Z G,Li J L.Grassland coverage inter-annual variation and its coupling relation with hydrothermal factors in China during 1982-2010.Journal of Geographical Sciences,2014,24(4):593-611.

[18] 杜子濤,占玉林,王長(zhǎng)耀.基于NDVI序列影像的植被覆蓋變化研究.遙感技術(shù)與應(yīng)用,2008,23(1):47-50.

Du Z T,Zhan Y L,Wang C Y.Study on vegetation-cover changes based on NDVI serial images.Remote Sensing Technology and Application,2008,23(1):47-50.(in Chinese)

[19] 任又成,張志斌.三江源生態(tài)逆向演變現(xiàn)狀調(diào)查及防治對(duì)策.安徽農(nóng)業(yè)科學(xué),2014,15(1):4720-4763.

Ren Y C,Zhang Z B.Sanjiangyuan Ecological Reverse Evolution Survey and Countermeasures.Journal of Anhui Agricultural Sciences,2014,15(1):4720-4763. (in Chinese)

[20] 曹旭娟,干珠扎布,梁艷,高清竹,張勇,李玉娥,萬運(yùn)帆,旦久羅布.基于NDVI的藏北地區(qū)草地退化時(shí)空特征分析.草業(yè)學(xué)報(bào),2016,25(3):1-8.

Cao X J,Ganzhuzhabu,Liang Y,Gao Q Z,Zhang Y,Li Y E,Wan Y F,Danjiuluobu.Temporal and spatial distribution of grassland degradation in northern Tibet based on NDVI.Acta Prataculturae Sinica,2016,25(3):1-8.(in Chinese)

[21] 胡云峰,徐芝英.蒙古高原地形與植被指數(shù)的特征尺度及多尺度耦合關(guān)系.地理科學(xué),2014,34(12):1511-1517.

Hu Y F,Xu Z Y.Characteristic scales and the multi-scale correlation between DEM and NDVI factors:A case study in the Mongolian Plateau.Scientia Geographica Sinica,2014,34(12):1511-1517.(in Chinese)

[22] 中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn)GB19377:天然草地退化、沙化、鹽漬化的分級(jí)標(biāo)準(zhǔn).中華人民共和國(guó)國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局,2003.

GB19377-2003:Parameters for degradation,sandification and salification of rangelands.General Administration of quality supervision,inspection and Quarantine of the People’s Republic of China,2003.(in Chinese)

[23] 胡云峰,徐芝英.蒙古高原地形與植被指數(shù)的特征尺度及多尺度耦合關(guān)系.地理科學(xué),2014,34(12):1511-1517.

Hu Y F,Xu Z Y.Characteristic scales and the multi-scale correlation between DEM and NDVI factors:A case study in the Mongolian Plateau.Scientia Geographica Sinica,2014,34(12):1511-1517.(in Chinese)

[24] Bao G,Bao Y H,Qin Z H,Xin X P,Bao Y L,Bayarsaikan S,Zhou Y,Chuntai B.Modeling net primary productivity of terrestrial ecosystems in the semi-arid climate of the Mongolian Plateau using LSWI-based CASA ecosystem model.International Journal of Applied Earth Observation and Geoinformation,2016,46:84-93.

[25] Chen X Q,Tan Z J,Schwartz M D,Xu C X.Determining the growing season of land vegetation on the basis of plant phenology and satellite data in northern China.International Journal of Biometeorology,2000,44(2):97-101.

[26] 高清竹,李玉娥,林而達(dá),江村旺扎,萬運(yùn)帆,熊偉,王寶山,李文福.藏北地區(qū)草地退化的時(shí)空分別特征.地理學(xué)報(bào),2005,60(6):965-973.

Gao Q Z,Li Y E,Lin E D,Jiang C W Z,Wan Y F,Xiong W,Wang B S,Li W F.Temporal and spatial distribution of grassland degradation in northern Tibet.Acta Geographic Sinia,2005,60(6):965-973.(in Chinese)

[27] 張宏斌,唐華俊,楊桂霞,李剛,陳寶瑞,辛?xí)云?2000-2008年內(nèi)蒙古草原MODIS NDVI時(shí)空特征變化.農(nóng)業(yè)工程學(xué)報(bào),2009,25(9):168-175.

Zhang H B,Tang H J,Yang G X,Li G,Chen B R,Xin X P.Change of spatial-temporal characteristics based on MONDIS NDVI data in Inner Mongolia grassland from 2000 to 2008.Society of Agricultural Engineering,2009,25(9):168-175.(in Chinese)

[28] 包剛,秦志豪,包玉海,周義.1982-2006年蒙古高原植被覆蓋時(shí)空變化分析.中國(guó)沙漠,2013,33(3):918-927.

Bao G,Qin Z H,Bao Y H,Zhou Y.Spatial-temporal changes of vegetation cover in Mongolia Plateau during 1982-2006. Journal of Desert Research,2013,33(3):918-927.(in Chinese)

[29] 周錫飲,師華定,王秀茹,孟凡浩.蒙古高原近30年來土地利用時(shí)空變化特征與動(dòng)因分析.浙江農(nóng)業(yè)學(xué)報(bào),2012,24(6):1102-1110.

Zhou X Y,Shi H D,Wang X R,Meng F H.Study on the temporal and spatial dynamic changes of land use and driving forces analysis of Mongolia Plateau in recent 30 years.Acta Agriculturae Zhejiangensis,2012,24(6):1102-1110.(in Chinese)

[30] 繆麗娟,蔣沖,何斌,劉強(qiáng),朱楓,崔雪鋒.近10年來蒙古高原植被覆蓋變化對(duì)氣候的響應(yīng). 生態(tài)學(xué)報(bào),2014,34(5):1295-1301.Miu L J,Jiang C,He B,Liu Q,Zhu F,Cui X F. Response of vegetation coverage to climate change in Mongolian Plateau during recent 10 years.Acta Ecologica Sinica,2014,34(5):1295-1301.(in Chinese)

[31] 穆少杰,李建龍,陳奕兆,剛成誠(chéng),周偉,居為民.2001-2010年內(nèi)蒙古植被覆蓋度時(shí)空變化特征,地理學(xué)報(bào),2012,67(9):1255-1267.

Mu S J,Li J L,Chen Y Z,Gang C C,Zhou W,Ju W M.Spatial Differences of variations of vegetation coverage in Inner Mongolia during 2001-2010.Acta Geographica Sinica,2012,67(9):1255-1267.(in Chinese)

[32] 穆少杰,李建龍,周偉,楊紅飛,章超斌,居為民.2001-2010年內(nèi)蒙古植被凈初級(jí)生產(chǎn)力的時(shí)空格局及其與氣候的關(guān)系. 生態(tài)學(xué)報(bào),2013,33(12):3752-3764.

Mu S J,Li J L,Zhou W,Yang H F,Zhang C B,Ju W M.Spatial-temporal distribution of net primary productivity and relationship with climate factors in Inner Mongolia from 2001 to 2010.Acta Ecologica Sinica,2013,33(12):3752-3764.(in Chinese)

[33] 額爾敦格日樂,包剛,包玉龍,杭玉玲，布仁吉日嘎拉.2001-2013年西鄂爾多斯國(guó)家級(jí)自然保護(hù)區(qū)植被覆蓋變化.水土保持研究,2016,23(1):110-116.

Erdengerel,Bao G,Bao Y L,Hang Y L,Burenjirigala.Vegetation Coverage Changes in West Ordos National Natural Reserves During the Period from 2001 to 2013.Research of Soil and Water Conservation,2016,23(1):110-116.(in Chinese)

[34] 劉海江,尹思陽,孫聰,等. 2000-2010年錫林郭勒草原NPP時(shí)空變化及其氣候響應(yīng).草業(yè)科學(xué),2015,32(11):1709-1720.

Liu H J,Yin S Y,Sun C,Peng F L,Zhou P.Temporal and spatial variation of net primary productivity (NPP) and its responses with climatic changes in the Xilingol grassland from 2000 to 2010.Pratacultural Science,2015,32(11):1709-1720.(in Chinese)

[35] 元志輝,包剛,銀山,雷軍，包玉海，薩楚拉. 2000-2014年渾善達(dá)克沙地植被覆蓋變化研究.草業(yè)學(xué)報(bào),2016,25(1):33-46.

Yuan Z H,Bao G,Yin S，Lei J,Bao Y H,Sa C L.Vegetation changes in Otindag sand Country during 2000-2014.Acta Parataculturae Sinica,2016,25(1):33-46.(in Chinese)

[36] 包剛,包玉海,覃志豪,周義,Shiirev-Adiya.近10年蒙古高原植被覆蓋變化及其對(duì)氣候的季節(jié)響應(yīng).地理科學(xué),2013,33(5):613-621.

Bao G,Bao Y H,Qin Z H,Zhou Y,Shiirev-Adiya.Vegetation cover changes in Mongolian Plateau and Its response to seasonal climate change in recent 10 years. Scientia Geographica Sinica,2013,33(5):613-621.(in Chinese)