1982-2015 年蒙古國植被覆蓋變化及其與氣候變化和人類活動的關系

杜佳夢，包剛，2，3*，佟斯琴，2，黃曉君，2，溫都日娜，美麗，包玉海，2

（1.內蒙古師范大學地理科學學院，內蒙古呼和浩特010022；2.內蒙古師范大學內蒙古自治區遙感與地理信息系統重點實驗室，內蒙古 呼和浩特010022；3.寧夏大學西北土地退化與生態恢復省部共建國家重點實驗室培育基地，寧夏 銀川750021）

全球氣候變化與陸地生態系統間的關系（global change and terrestrial ecosystem，GCTE）是全球變化研究的重要內容之一［1］。植被是陸地生態系統的重要組成部分，連接著土壤、水圈和大氣圈，促進各個圈層物質遷移和能量交換，在氣候調節、陸地碳循環、水土保持等方面具有重要作用［2?3］。氣候變化是改變植被生長、結構和功能的主要因素之一［4］。植被對氣候變化的響應研究已成為全球變化研究的主要內容之一［5?6］。

蒙古國是蒙古高原的主要組成部分［7］，也是中國“一帶一路”合作戰略和中蒙俄經濟走廊的重要區域［8］。蒙古國位于干旱半干旱地區，加之其相對高的海拔和大陸性氣候，生態系統較為脆弱。植被類型主要以草原為主，占國土面積的79%，包括草甸草原、典型草原和荒漠草原，而森林約占國土面積的10%，主要分布在氣候相對濕潤的北部高緯度地區和山區［9］。此外，蒙古國是世界上人口稀少的國家之一［10］，蒙古國國家統計局數據顯示，2019年蒙古國人口為329.69 萬，平均人口密度為2.1 人·km?2，人類活動對生態系統的影響相對較少，加之其獨特的地理位置、廣闊的草原及其多種植被類型等為區域尺度植被生長狀況及其氣候變化響應研究創造了良好的研究條件。蒙古國南部與中國內蒙古接壤，其生態環境、自然資源和社會人文等與中國有著緊密關聯，屬于中國北部最重要的綠色生態屏障。近些年蒙古國陸地生態系統受氣候變化影響顯著，如湖泊面積萎縮［11］、干旱加?。?2］、森林草原火災［13］和森林病蟲害［14］等自然災害頻發，可能直接或間接影響中國北方生態環境狀況。因此，有必要在蒙古國國家尺度上開展近30 多年來植被覆蓋變化及其與氣候要素和人類活動因素關系的研究。

本研究利用1982?2015 年全球監測與模型研究組（global inventory modeling and mapping studies，GIMMS）第三代NDVI（normalized difference vegetation index）數據（NDVI3g），結合同期溫度、降水量及植被類型數據，分析蒙古國植被NDVI 的年際、季節變化特征，在此基礎上探討其對氣候變化和人類活動的響應。旨在揭示不同時空尺度下蒙古國植被生長變化趨勢及其主要驅動，為理解該地區生態環境的總體發展趨勢和對中國“三北”地區生態環境建設以及防災減災等方面提供參考依據。

1 材料與方法

1.1 研究區概況

蒙古國位于亞洲中部（41°35′?52°09′N，87°44′?119°56′E），屬內陸高原國家，國土面積為約156 萬km2［10］。地勢由西向東逐漸降低（圖1a），平均海拔在1580 m 左右［7］。氣候類型為典型的大陸性干旱半干旱氣候［9］，冬季漫長而寒冷，夏季炎熱干燥，年平均溫度在?8～6 ℃波動，由南向北降低（圖1b）。年降水量的分布與溫度相反（圖1c），從北部山區的>350 mm 降低到南部戈壁荒漠地區的<50 mm［15］。受氣候分布的影響，蒙古國植被類型由北向南依次為森林、草原和荒漠稀疏植被［16］（圖1d）。在氣候和植被類型的影響下，草地畜牧業是蒙古國最主要的經濟基礎，支撐著80%以上的農牧業產值［17］。

圖1 蒙古國高程（a）、溫度（b）、降水（c）及植被類型（d）的空間分布格局Fig.1 Spatial distribution pattern of elevation（a），temperature（b），precipitation（c）and vegetation type（d）in Mongolia

1.2 數據來源及預處理

本研究使用的遙感數據為1982?2015 年間的GIMMS NDVI3g V1 數據（https：//ecocast.arc.nasa.gov/da?ta/pub/gimms/3g.v1/），時間和空間分辨率分別為15 d 和0.0833°［18］。研究區大部分植被在冬季停止生長或受積雪覆蓋的影響等，每年選取4?10 月（即生長季）的NDVI 數據用于研究［16］。對每月兩景數據進行最大值合成處理［19?20］，獲得月尺度的NDVI 時間序列數據集，以進一步消除15 d NDVI 中存在的云、氣溶膠和其他非植被因素等的影響。在此基礎上，對4?5、6?8 和9?10 月NDVI 分別進行平均值計算，得到春、夏和秋季NDVI［21］，用于分析不同季節NDVI 變化趨勢及其氣候響應。對生長季3 個季節NDVI 進一步計算均值獲得生長季平均NDVI，分析年NDVI 的時空變化趨勢。由于在植被覆蓋度極低的地區NDVI 容易受到土壤背景的影響而不能準確反映植被覆蓋的實際信息，在此將多年平均NDVI 小于0.075 的地區定為無植被區，不做分析［10］。

氣象數據由蒙古國科學院地理與地質生態研究所提供，為蒙古國133 個氣象臺站的1982 年1 月1 日?2015年12 月31 日溫度和降水量數據。對日氣象數據進行處理得到月數據后，在ArcGIS 10.2 環境下對氣象數據進行克里金插值，獲得時間和空間分辨率與NDVI 數據一致的柵格化氣象數據集，用于像元尺度的NDVI 變化趨勢及其對氣候的響應研究。

植被類型數據來源于蒙古國國家地圖集。對地圖集中的植被類型數據進行掃描、幾何校正、數字化等處理，獲得森林、草甸草原、典型草原、荒漠草原、戈壁荒漠和高山草地6 種植被類型，用于分析不同植被類型NDVI 變化趨勢及其氣候響應。牲畜數量和人口數據均從蒙古國國家統計信息服務網（http：//www.1212.mn/）獲取。

1.3 研究方法

采用一元線性回歸法（公式1）在區域尺度和像元尺度上分別分析研究區1982?2015 年生長季和生長季3 個季節（春、夏和秋季）NDVI 時空變化趨勢［22］。

式中：b代表斜率，即變化趨勢；NDVIi表示第i年的生長季或各季節NDVI。根據斜率和顯著性水平值，將變化趨勢分為以下4 個等級，分別為顯著增加（b>0，P<0.05）；不顯著增加（b>0，P>0.05）；顯著減小（b<0，P<0.05）；不顯著減?。╞<0，P>0.05）。

為了分析研究區植被覆蓋變化（NDVI）對溫度和降水量的響應特征，分別對生長季和生長季三個季節平均NDVI 與溫度和降水量數據進行偏相關性分析［5］，先計算各要素的簡單相關關系，再計算偏相關關系。

式中：x為NDVI 值；y為溫度和降水；Rxy為x和y間的相關系數；xi為第i年或季節的平均NDVI 值；yi為第i年或季節的平均溫度或降水量為研究時段內年或季節平均NDVI 值為研究時段內年或季節平均溫度或降水量。

式中：Rab，c為將變量c 固定后變量a 與變量b 間的偏相關系數，rab、rac和rbc分別表示變量a 與b、變量a 與c、變量b 與c的相關系數。

在氣象因子和NDVI 進行偏相關分析的基礎上，再采用殘差分析法［23］，探討人類活動對植被生長的影響。在像元尺度上建立溫度和降水的回歸模型獲得逐年NDVI 的預測值（公式4）。然后，將從遙感影像提取的NDVI 值假設為實測值，減去NDVI 預測值，得到逐年NDVI 殘差（公式5），以表示人類活動引起的植被生長變化［20］。

式中：NDVI和NDVI′分別為遙感觀測的NDVI 實測值和模型計算的NDVI 預測值；T為生長季平均溫度；P為生長季累計降水量；ε為殘差。

2 結果與分析

2.1 蒙古國生長季多年平均NDVI 空間格局

1982?2015 年蒙古國植被生長季平均NDVI 整體上呈由北向南逐漸減小的空間分布特征（圖2）。蒙古國北部地區主要分布著森林植被［24］，NDVI 值大于0.6。由北至南植被類型依次跨越草甸草原、典型草原、荒漠草原和戈壁荒漠稀疏植被［16］，NDVI 值也隨植被類型的變化由0.8 下降至0.1 以下。為了更直觀體現研究區NDVI 南北地域差異，選取3 條南北剖面，分析了NDVI 值隨剖面的變化趨勢（圖3）。蒙古國植被NDVI 在3 個剖面上從南向北的上升速率分別為0.2°/10° N，0.7°/10° N 和0.7°/10° N，具有顯著的緯度地帶性分布特征，尤其在蒙古國中部和東部此特征更加顯著（R=0.91，P<0.01；R=0.93，P<0.01）。而西部緯度地帶性相對較弱（R=0.47，P<0.01）（圖3a，3b 和3c），可能與其較高的海拔有關。在蒙古國西部分布著杭愛山、戈壁阿爾泰山和蒙古阿爾泰山等山脈［25］，在一定程度上可能模糊了此區域植被NDVI 的緯度地帶性分布。從高程梯度上的變化看（圖3d），隨著海拔的升高，蒙古國ND?VI 以0.2/100 m 的速率降低，但存在3 個不同階段，即500～1300 m NDVI 隨海拔上升而減少；與此相反，在1300～2200 m 間NDVI 隨海拔的上升而增加；而在2200 m 之后，隨著海拔升高NDVI 值下降趨勢明顯。其原因可能是因為隨著海拔的升高溫度逐漸降低，導致植被生長所需的水熱條件下降，NDVI 降低。蒙古國34 年生長季平均NDVI 為0.26，低于內蒙古2001?2010 年平均NDVI 值（0.38）［5］，說明蒙古國植被平均生產力總體上低于內蒙古平均生產力。

圖2 1982-2015年蒙古國植被生長季平均NDVI空間分布格局Fig.2 Spatial distribution pattern of mean growing season NDVI in Mongolia from 1982 to 2015

圖3 1982-2015 年蒙古國植被生長季平均NDVI 緯度（a～c）和海拔（d）梯度上的變化特征Fig.3 Variations of growing season NDVI along latitude（a-c）and altitude（d）gradients in Mongolia from 1982 to 2015

2.2 蒙古國生長季平均NDVI 變化及其氣候響應

2.2.1 區域尺度上的NDVI 變化趨勢及其氣候響應 1982?2015 年，蒙古國植被生長季NDVI 介于0.24～0.29，每年以0.0005（R=0.43，P<0.05）的速率呈顯著增加趨勢，最低和最高值分別出現在2007 和2012 年（圖4）。然而，過去34 年間，NDVI分別經歷了3 個階段的變化趨勢：1）1982?1994 年的 顯 著 增 加 趨 勢，其增加速率為0.0017·a?1（R=0.64，P<0.05）；2）1994?2007 年的顯著下降趨勢，下降速率為0.0016·a?1（R=?0.73，P<0.01）；3）2007?2015 年的顯著增加趨勢，增加速率為0.0044·a?1（R=0.80，P<0.01），且其增加幅度明顯大于第一階段。蒙古國生長季NDVI 變化趨勢1994 年發生轉折（即從上升到下降趨勢的轉折），2007 年發生第二次轉折，由下降轉為上升趨勢，且上升速率高于90 年代之前的上升速率。從圖4 可看出，降水是改變蒙古國植被NDVI 變化及其轉折的主要驅動力，即蒙古國生長季降水量同樣在1982?1994、1994?2007 和2007?2015 年分別呈增加、減小和再增加趨勢，與NDVI 的波動和趨勢高度一致，兩者之間的偏相關性為0.74（P<0.01）。而蒙古國生長季溫度總體上呈持續升高趨勢，其中，1982?2000 年的增溫趨勢最明顯，2000 年以來在保持高溫狀態的同時，增溫態勢略有緩慢。其與NDVI 的偏相關系數為0.60（P<0.01）。需要指出的是，蒙古國NDVI 最低值出現的2007 年正好對應于降水較低和溫度較高的年份，說明兩者導致的干旱是可能導致2007 年NDVI 最低的直接因素（圖4）。

圖4 1982-2015年蒙古國生長季NDVI、降水量和溫度的年際變化Fig.4 Interannual variation of NDVI，precipitation and tem?perature in growing season of Mongolia from 1982 to 2015

圖5 蒙古國1982-2015（a）、1982-1994（b）、1994-2007（c）和2007-2015 年（d）生長季NDVI 變化趨勢的空間分布Fig.5 Spatial distribution of growing season NDVI change trend in Mongolia from 1982 to 2015（a），1982 to 1994（b），1994 to 2007（c）and 2007 to 2015（d）

2.2.2 蒙古國生長季平均NDVI 變化趨勢的空間格局及其氣候響應 從NDVI 變化趨勢的空間分布特征看，1982?2015 年，蒙古國生長季NDVI 呈增加和減小趨勢的面積分別占其國土面積的66.6%和33.4%（表1），其中顯著增加（P<0.05）的地區為32.8%，主要集中分布在東部開闊的草原區和零星分布在西部和北部山區；而顯著減小（P<0.05）的地區僅占研究區8.9%，零星分布在中部和西南部地區（圖5a）。與國家尺度的變化趨勢一致（圖4），1982?1994 年蒙古國生長季NDVI 主要以增加趨勢為主，約占研究區88.1%，其中顯著增加的面積為22.8%，主要分布在西部大湖盆地、杭愛山脈、西南部戈壁阿爾泰山脈和北部森林地區；呈減小趨勢的面積約為11.9%，但大多并不顯著（圖5b）。1994?2007 年，蒙古國生長季NDVI 主要以減小為主，約占80.4%，其中顯著減小的面積為22.9%，主要集中在中東部典型草原區及與稀疏植被的交界處；而增加趨勢的面積為19.6%，零星分布在西部、北部和南部地區（圖5c）。2007?2015 年，NDVI 呈增加趨勢的面積（78.2%）同樣遠高于呈減小趨勢的面積（21.8%），其中呈顯著增加和顯著減小趨勢的面積分別為26.8%和2.2%。呈增加趨勢的地區主要分布于蒙古國東部、中部和西部一線上，而減小趨勢的地區分布在其兩側，即蒙古國南部和北部地區，但大多像元的減小趨勢并不顯著（圖5d）。像元尺度的相關分析表明，約86.7%地區的NDVI 與降水量呈正相關（其中顯著正相關為55.1%），主要分布在除北部森林和西部高山地區外的蒙古國大多地區（圖6a）。NDVI 與溫度相關性的空間分布總體與降水的空間格局相反，即在與降水呈正相關的地區NDVI 與溫度呈負相關（40.3%），而在與降水呈負相關的地區NDVI 與溫度呈正相關（59.7%）（圖6b）。在北部與溫度正相關（與降水呈負相關）的地區可能與該地區降水充足但年平均溫度較低有關［10］（圖1），即溫度的升高有利于這些地區植被生長［2，26］，而降水量的增加可能會增加云量和降低溫度而抑制植被生長［10，27］。其余大多與溫度呈負相關（與降水正相關）的地區可能與溫度的升高將增加蒸散量而加劇這些原本就水資源匱乏地區的干旱，從而抑制有關植被生長［28］。為了對比降水和溫度對NDVI 產生影響的大小，在此對兩個偏相關系數求平方，再計算兩者差值（即降水偏相關系數平方減去溫度偏相關系數平方）（圖6c）。可以看出，在南部、中部和東部草原地區降水的正效應大于溫度的負作用（75.2%），而北部森林分布區及其鄰近地區溫度的正作用大于降水的負作用（24.8%）。這說明蒙古國大多地區主要處于干旱半干旱氣候區，降水對植被生長的影響范圍遠大于溫度的影響范圍。

圖6 1982-2015 年蒙古國生長季NDVI 與降水（a）和溫度（b）的偏相關性及兩個偏相關系數的平方之差（c）Fig.6 Partial correlation between growing season NDVI and precipitation（a），temperature（b）and the difference be?tween the squares of the two partial correlation coefficients（c）in Mongolia from 1982 to 2015

2.3 蒙古國NDVI 和對氣候變化響應的季節差異

2.3.1 區域尺度上的各季節NDVI 變化及其氣候響應 過去34 年，蒙古國生長季各季節平均NDVI 都呈現出增加趨勢（圖7）。其中，春季NDVI 的增加趨勢最為明顯，每年以0.0005（P<0.05）的速率增加，且其波動與春季溫度的波動高度一致（圖7a），兩者偏相關系數為0.64（P<0.01）（圖7d）。這說明溫度是影響蒙古國植被春季生長最主要的氣候因子，與北半球大多地 區植被 變化及其 主要驅動力高度 一致［3，26?29］。需要指出的是蒙古國植被春季NDVI 變化趨勢沒有明顯的轉折，而隨著春季溫度的持續升高NDVI 也呈持續增加趨勢。1982?2015 年蒙古國春季降水量也呈增加趨勢，其中，只有1988、2002、2003、2010 和2014 年為春季豐水年份。在34 年的觀測尺度中，夏季NDVI 的增加速率約為0.0003·a?1（P=0.27）（圖7b），且其波動與生長季NDVI 的波動（圖4）高度一致，也經歷了3 個階段的變化趨勢，即1）1982?1994 年呈顯著增加趨勢（速率為0.0025·a?1，R=0.57，P<0.05）；2）1994?2007 年的顯著下降趨勢（速率為0.0034·a?1，R=?0.73，P<0.01）；3）2007?2015 年的顯著增加趨勢（速率為0.0068·a?1，R=0.72，P<0.05）。這可能與夏季是包括蒙古國在內的北半球植被生長的旺季有關。Bao 等［30］的研究表明，蒙古高原植被生產力的70%以上都集中在夏季，其變化直接決定年生產力的變化趨勢，而不是由春季返青期提前或秋季枯黃期推遲導致的生長季延長引起［31?32］。同樣，與生長季NDVI 變化趨勢的主要驅動力一致，夏季降水量是改變夏季NDVI 變化的最主要氣候因子，兩者偏相關系數為0.65（P<0.01）（圖7d）。從氣候對植被生長影響的滯后效應看，盡管不顯著，春季降水的增加將促進夏季NDVI 的增加（R=0.17，P=0.36），而春季溫度的升高將抑制夏季NDVI 的增加（R=?0.10，P=0.67）（圖7d）。Bao 等［30］的研究還表明，蒙古高原秋季植被生產力約占年總生產力的20%，其變化對年總生產力的影響比春季生產力要大。過去34 年蒙古國秋季平均NDVI（0.30）略大于春季（0.22），且呈現出不顯著的增加趨勢（速率為0.0003·a?1，R=0.17，P=0.35）。與春季NDVI 一致，溫度是改變蒙古國秋季NDVI 的主要氣候驅動因子（R=0.42，P<0.05）（圖7d），但其強度低于春季NDVI 與春季溫度的相關性。盡管不顯著，秋季降水量對秋季NDVI 的影響也不可忽略。從圖7c 可以看出，2002?2010 年秋季NDVI明顯比其他年份較低，與其降水量負距平年份高度一致。值得注意的是，夏季降水與秋季NDVI 呈顯著的正相關關系（R=0.68，P<0.01）（圖7d），其原因是夏季降水的增加將增加夏季植被生長，然而夏季植被的生長對秋季生長具有重要的正“遺產效應”。進一步的統計分析也表明，夏季與秋季NDVI 之間存在顯著的正相關（R=0.76，P<0.01）關系，而春季NDVI 與夏季和秋季NDVI 之間并沒有顯著的相關關系（春季NDVI 與夏季NDVI：R=0.22，P=0.22；春季NDVI 與秋季NDVI：R=0.12，P=0.52）。

表1 蒙古國不同時段NDVI 變化趨勢的面積百分比Table 1 Area percentage of NDVI change trend in Mongolia during different periods（%）

圖7 1982-2015 年蒙古國春季（a）、夏季（b）和秋季（c）NDVI、降水量和溫度的年際變化及不同季節NDVI 與氣候因子的偏相關關系（d）Fig.7 The interannual changes of NDVI，precipitation and temperature in spring（a），summer（b），autumn（c）and the partial correlation between NDVI and climate factors in different seasons（d）from 1982 to 2015 in Mongolia

圖8 1982-2015 年蒙古國各季節NDVI 與溫度和降水偏相關性的空間分布Fig.8 Spatial distribution of partial correlation between NDVI and precipitation and temperature in various seasons in Mongolia from 1982 to 2015

2.3.2 蒙古國各季節NDVI 對氣候響應的空間格局 在空間上，春季NDVI 與春季溫度呈正相關的面積（87.0%）遠大于呈負相關的面積（13.0%），其中顯著正相關的地區主要分布在蒙古國北部杭愛山脈、薩彥嶺、肯特山和西部阿爾泰山脈等高海拔和高緯度地區（圖8a）。而春季NDVI 與降水呈正相關的面積（58.2%）和呈負相關的面積（41.8%）相差不大。呈正相關的地區主要分布在蒙古國東部開闊的草原地區，而呈負相關的地區分布在西部和西南地區（圖8b）。夏季NDVI 與溫度的相關性主要呈負相關（55.9%）（圖8c），較高的溫度可能通過增加蒸發造成水分不足，從而在一定程度上抑制植被生長［5］。與前面在區域尺度上的討論一致，由于夏季植被生產力在全年生產力中的比重高達70%以上［30］，夏季NDVI 與降水的空間分布（圖8d）與生長季觀測到的結果（圖6a）高度一致，說明夏季氣候因子是決定全年植被生長的最主要因子［33］，而春季和秋季氣候因子主要對植被年生長開始時期和結束時期的生長產生影響［34］，而對全年生長的直接影響相對要小［30］。秋季NDVI 與溫度的相關性總體上以負相關為主（54.6%），分布在蒙古國西部、中部和東部地區，而呈正相關的地區集中在北部寒冷地區（圖8e）。在空間分布上，秋季NDVI 和降水的關系大致與秋季NDVI和溫度的關系相反（圖8f），即北部寒冷地區主要以負相關為主，此區域年均降水量相對較多，降水的增加會減少日照時數，從而抑制植被生長。而西部、中部和東部草原區主要以正相關為主。

2.3.3 蒙古國不同植被類型NDVI 對氣候的響應從圖9 可以看出，在整個生長季內蒙古國除荒漠草原外，其他5 種植被類型與溫度均呈正相關，其中森林、典型草原和高山草地為顯著正相關（圖9a）。與降水的相關分析表明，森林植被NDVI 與降水呈負相關（?0.14），其余植被類型呈正相關（圖9b）。在生長季3 個季節中，所有植被類型春季NDVI 與溫度呈顯著正相關（除了戈壁荒漠稀疏植被不顯著外），說明溫度是蒙古國所有植被類型春季開始生長的最主要驅動力。其中春季溫度對草原植被NDVI 的影響程度從大到小依次為草甸草原、典型草原、荒漠草原，說明在熱量較低的地區，春季溫度對植被的影響更重要［16］。降水也是對草原植被的春季開始生長和后續生長具有重要促進作用（圖9b）。由于夏季是3 個季節中溫度最高的季節，總體上抑制植被的生長，而夏季降水對植被生長具有重要的正效應（除森林外）。秋季氣溫和降水對秋季所有植被類型的NDVI 影響都不是很明顯，而秋季NDVI 對夏季降水的滯后效應較顯著（圖9c），與Piao 等［21］的中國溫帶草原和荒漠草原秋季NDVI 與前一季的降水量均呈正相關關系的研究結果基本一致。

圖9 不同植被類型NDVI 與溫度和降水之間的偏相關關系Fig.9 Partial correlations between NDVI of different vegeta?tion types with temperature and precipitation

2.4 蒙古國生長季平均NDVI 對人類活動的響應

圖10 蒙古國NDVI 殘差變化趨勢（a）及其顯著性（b）的空間分布Fig.10 Spatial distribution of NDVI residual variation trend（a）and its significance（b）in Mongolia

氣候是導致植被覆蓋變化的最重要的因素，但人類活動也通過放牧、開采、植樹造林和自然保護區的建立等來改變植被覆蓋變化［9，25］。從殘差分析結果看（圖10），在1982?2015 年，蒙古國NDVI 殘差值呈上升趨勢（代表人類活動對植被的改善作用）的面積為74.2%，其中顯著增加的區域約為31.3%，主要分布于蒙古國東部與中國內蒙古自治區接壤的草原區，在其他研究中也發現相似結果，可能是該區域牲畜總數和人口密度呈下降趨勢的原因［36］。而殘差值呈下降趨勢（代表人類活動對植被的破壞作用）的地區零星分布于人口密度相對較大的蒙古國中部草原區和西部高山地區，這可能與這些地區放牧、開礦等人類活動密切相關［9，37］。從圖11可看出，1982?2015 年蒙古國牲畜數量從247.65 萬頭大幅增加到559.80 萬頭，2001?2015 年，人口從243.24 萬增加到305.78 萬（蒙古國國家統計局）。特別是2002?2009 年蒙古國牲畜數量急劇上升，在一定程度上與此時間段的NDVI 較低值吻合（圖4）。

圖11 蒙古國牲畜數量和人口的變化趨勢Fig.11 Trends of livestock quantity and population in Mongolia

3 討論

本研究表明，1982?2015 年蒙古國植被生長季平均NDVI 空間上呈由北向南逐漸減小的分布特征（圖2），可能是北部地區受北冰洋水汽的影響［38］，其降水量較充足，但隨著向南離北冰洋距離的增加，降水量也逐漸降低。34 年間蒙古國生長季NDVI 呈顯著增加趨勢（圖4），然而，在研究期間，這種增長趨勢并不是持續的，而是有兩次轉折。第一次趨勢轉折發生在1994 年，即從上升到下降趨勢的轉折，總體上與北半球中高緯度觀測到的結論高度一致。如Peng 等［39］、Piao 等［40］和張學珍［41］的研究表明，在中國、歐亞大陸溫帶地區和北半球中緯度地區的植被NDVI 分別以1990、1997 和1994 年為轉折點，呈現先增加后下降的變化趨勢。第二次趨勢轉折發生在2007 年，由下降轉為上升趨勢。這說明在氣候變暖背景下的區域NDVI 變化趨勢不是持續的增加或減小，而可能每隔10～15 年時間尺度上發生一次轉折性變化，其背后的驅動機制較為復雜，涉及多個環境因素和人類活動的共同作用［42］。

過去大多研究表明，氣候變暖是導致北半球中高緯度大多地區植被NDVI 變化及其轉折的主要因素［39?42］。而本研究中蒙古國植被生長季NDVI 變化主要受降水變化的影響（圖6）。這可能與蒙古國地處干旱半干旱地區有關［5，43?44］。另外，森林植被生長與降水呈負相關，說明在整個生長季內降水的增加可能通過云量、降低溫度和輻射等來抑制分布于相對濕潤和寒冷地區的森林植被生長［26］。而溫度是影響蒙古國植被春季和秋季生長最主要的氣候因子（圖7），這可能與春季溫度的升高導致植被返青期提前［37］，秋季溫度的升高可以減緩或推遲秋季植被枯黃期有關［34，45］，即返青期的提前或枯黃期推遲將延長植被生長天數而增加植被NDVI［30，32］。但春季溫度的升高抑制夏季植被的生長，這與Buermann 等［46］的研究結論較一致，其在北美地區的研究表明，春季溫度的升高使植被返青期提前，但返青期的提前將增加蒸散發而加劇生態系統干旱，最終導致降低夏季植被的生長。然而植被覆蓋變化在受到氣候影響的同時，人類活動也對其產生了不可忽視的影響。魏云潔等［17］的研究表明，蒙古國國民生產總值的26.2%來自農牧業，而農牧業產值中的80%來自畜牧業，說明支撐畜牧業生產力的草原放牧活動主要分布在人口和牲畜相對密集的蒙古國中部和西部地區，這將降低植被生產力，在空間上與殘差值下降趨勢基本吻合（圖10）。同時蒙古國東部植被明顯的綠化趨勢，與Zhou 等［18］的研究一致?？赡苁怯捎诿晒艊鴸|部擁有蒙古達烏爾保護地［47］和東方蒙古草原自然保護區［48］等國家自然保護區的原因。此外，色楞格省、布爾干省和中央省這3 個省植被綠化趨勢增加的原因可能與蒙古國政府發起“第三次開墾荒地運動”的計劃有關［49］。

4 結論

利用1982?2015 年的GIMMS NDVI 數據、同期月氣溫和降水量數據及植被類型數據，系統分析了蒙古國植被覆蓋年際、季節變化特征，并從氣候和人類活動的角度分析其變化原因。通過研究得出了如下結論。

1）從34 年生長季平均NDVI 的空間分布看，蒙古國NDVI 由南向北逐漸增加，在西、中、東3 條剖面緯度每升高10° NDVI 分別增加0.2、0.7 和0.7，具有顯著的緯度地帶性分布特征。從海拔梯度上，海拔每升高100 m，ND?VI 下降0.2。

2）從時間變化趨勢看，過去34 年，蒙古國生長季NDVI 呈顯著的增加趨勢（0.0005·a?1），但經歷了3 個階段性變化趨勢，分別為1982?1994 年的顯著增加（0.0017·a?1），1994?2007 年的顯著下降（0.0016·a?1）和2007?2015 年的顯著增加（0.0044·a?1）。NDVI 變化趨勢及其階段性變化與降水量波段高度一致，兩者偏相關系數為0.74（P<0.01），說明降水量變化是蒙古國植被覆蓋變化的主要原因。從區域差異看，除北部小部分森林分布區NDVI 受溫度影響外，其他大多地區的NDVI 主要受降水的控制。

3）從季節差異看，總體上3 個季節NDVI 均呈增加趨勢，其中由于夏季是植被生長的旺季，其NDVI 波動與生長季NDVI 波動高度一致，即夏季降水量變化引起的NDVI 是決定全年植被覆蓋變化的最直接因素。與夏季不同，春季和秋季溫度是春、秋季植被生長的主要因素。各季節NDVI 及其與氣象因素的關系具有明顯的空間異質性。

4）殘差分析表明，人類活動主要對蒙古國中西部地區的NDVI 變化產生負影響，與蒙古國人口密度的空間分布基本一致。