侯志雄， 井長(zhǎng)青， 王公鑫， 郭文章， 趙葦康

（1.新疆農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，新疆草地修復(fù)與環(huán)境信息重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830052；2.新疆農(nóng)業(yè)大學(xué)草業(yè)學(xué)院，烏魯木齊 830052）

全球氣候變化是全球變化與陸地生態(tài)系統(tǒng)（global change and terrestrial ecosystem，GCTE）研究的重要內(nèi)容，得到國(guó)內(nèi)外學(xué)者的高度重視[1]。全球氣候變化主要體現(xiàn)在對(duì)生態(tài)環(huán)境系統(tǒng)的變化、資源條件和生態(tài)脆弱敏感區(qū)的植被覆蓋度等方面的影響[2]。研究表明，伴隨著全球氣候變暖，中高緯度地區(qū)植被活動(dòng)顯著增強(qiáng)[3]。植被是地表植物群落的總稱，作為生態(tài)系統(tǒng)的重要組成部分，是土壤、大氣和水分的連接紐帶[4]，在促進(jìn)土壤形成與發(fā)育等方面起重要作用[5]。植被覆蓋度（fractional vegetation cover, FVC）是植被在地面的垂直投影面積占總面積的百分比[6]，可反映植被稀疏程度及地表狀況。植被作為表征全球氣候變化的綜合指示器，直接影響人類賴以生存的環(huán)境，植被結(jié)構(gòu)和功能對(duì)氣候變化的響應(yīng)能夠在一定程度上反映全球氣候變化的趨勢(shì)。植被與氣候的關(guān)系是生態(tài)系統(tǒng)研究過程中的重要內(nèi)容，因此，了解北疆天然草地植被與氣候的相互作用機(jī)制具有重要的科學(xué)意義[7]。

20 世紀(jì)70 年代，遙感技術(shù)快速發(fā)展，已成為人類大范圍、實(shí)時(shí)監(jiān)測(cè)地表植被強(qiáng)有力的工具[8]。Rouse 等[9]1974 年首次提出了歸一化植被指數(shù)（normalized difference vegetation index，NDVI）的概念，NDVI是近紅外反射（near infrared reflection，NIR）與紅光（red，R）之差除以NIR與R之和[10]。由于NDVI 計(jì)算簡(jiǎn)單，且可以盡可能的減少儀器定標(biāo)、地形、云等影響造成的誤差，對(duì)低植被覆蓋度識(shí)別敏感，是目前運(yùn)用最廣泛的植被指數(shù)[11]。Lamchin 等[12]研究了蒙古國(guó)植被覆蓋度的變化特征，分析了研究區(qū)植被覆蓋度的變化趨勢(shì)。李雙雙等[13]研究了陜甘寧地區(qū)植被的時(shí)空變化特征，指出近10 年植被覆蓋度呈波動(dòng)增長(zhǎng)。Wang 等[14]通過分析南非植被覆蓋度與同期月、季、年降雨量的關(guān)系，指出降雨量與不同植被類型間存在一定相關(guān)性，不同植被類型對(duì)月降雨量存在滯后性，總植被覆蓋度對(duì)年降雨量具有一定的敏感度。Weiss 等[15]分析了1990—2000 年NDVI 數(shù)據(jù)與同期氣象數(shù)據(jù)之間的關(guān)系，指出美國(guó)新墨西哥州的氣候與干旱、半干旱地區(qū)不同類型植被的關(guān)系。Milich 等[16]在2000 年研究了氣候因子與西非牧場(chǎng)之間的相關(guān)性，表明生長(zhǎng)季牧場(chǎng)植被與降水具有顯著相關(guān)性。張戈麗等[17]采用遙感反演呼倫貝爾近30 年的草地覆蓋度變化，分析表明該區(qū)草地覆蓋度與降雨的響應(yīng)程度較好。李曉光等[18]反演了鄂爾多斯地區(qū)植被覆蓋度與氣候之間的關(guān)系，表明在干旱區(qū)草地覆蓋度對(duì)降水的響應(yīng)呈顯著正相關(guān)，在冬季草地覆蓋度對(duì)氣溫的響應(yīng)呈顯著正相關(guān)。丁明軍等[19]分析了青藏高原區(qū)植被覆蓋度與氣候因子之間的相關(guān)性，表明植被與氣溫、降水的響應(yīng)度較好。目前對(duì)植被覆蓋度的研究主要集中在年際和生長(zhǎng)季，對(duì)生長(zhǎng)季不同月份間草地植被覆蓋度變化的研究相對(duì)較少。因此，研究北疆6、7、8和9 月天然草地植被覆蓋度變化及其與氣象因子的關(guān)系，對(duì)保護(hù)該地區(qū)草原生態(tài)、發(fā)展畜牧業(yè)和維持草地生態(tài)系統(tǒng)的穩(wěn)定性具有重要意義。

北疆位于我國(guó)西北邊沿，地處亞歐大陸腹地，包括天山以北、阿爾泰山以南的區(qū)域，是新疆經(jīng)濟(jì)繁榮的中心區(qū)域，屬于溫帶大陸性干旱/半干旱氣候，區(qū)域內(nèi)干旱缺水，生態(tài)環(huán)境極其脆弱[20-21]。本文選取1998—2018 年MODIS-NDVI 數(shù)據(jù)，結(jié)合同期北疆氣象數(shù)據(jù)和土地利用數(shù)據(jù)，揭示北疆天然草地植被覆蓋度時(shí)空分布特征及其與降水、氣溫、日照時(shí)數(shù)、相對(duì)濕度的關(guān)系，為北疆天然草地的可持續(xù)利用以及為干旱區(qū)生態(tài)環(huán)境治理提供參考依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

北疆地處亞歐大陸腹地，草地面積廣闊，資源豐富，是我國(guó)重要的天然草地牧區(qū)之一，地理位置介于43°23′—49°10′N，79°57′—91°32′E 之間，在行政區(qū)劃上包括烏魯木齊市、克拉瑪依市、石河子市、北屯市、昌吉回族自治州、伊犁哈薩克自治州、阿勒泰地區(qū)、塔城地區(qū)、博爾塔拉蒙古自治州等地區(qū)[22]。“三山夾兩盆”構(gòu)成了新疆獨(dú)特的地形地貌，從南到北依次為昆侖山脈、塔里木盆地（塔克拉瑪干沙漠）、天山山脈、準(zhǔn)噶爾盆地（古爾班通古特沙漠）、阿爾泰山。其中，北疆包括了阿爾泰山、準(zhǔn)噶爾盆地（古爾班通古特沙漠）、天山山脈（圖1）。北疆地形復(fù)雜，以準(zhǔn)噶爾盆地為中心，三面環(huán)山，形成明顯的地形單元，地勢(shì)西高東低，南北高中間低。中部的準(zhǔn)噶爾盆地，南邊廣闊，北邊窄小，海拔500 m 以上[23]；北部的阿爾泰山平均海拔在3 000 m 以上，植被垂直分異比較顯著；南部的天山平均海拔約4 000 m，由一系列平行山脈組成，占新疆全區(qū)總面積的1/3，山脈內(nèi)常年積雪，發(fā)育眾多冰川，南北兩側(cè)河流分布密度大[24]。

圖1 北疆高程圖Fig. 1 Altitude distribution in Northern Xinjiang

1.2 數(shù)據(jù)來源與處理

1.2.1 歸一化植被指數(shù)（NDVI）數(shù)據(jù) 從中國(guó)科學(xué)院資源環(huán)境與數(shù)據(jù)中心（www.resdc.cn）下載的NDVI數(shù)據(jù)為月度植被指數(shù)（NDVI）空間分布數(shù)據(jù)集[25]，時(shí)間為1998—2018 年，空間分辨率為1 km，是基于連續(xù)時(shí)間序列的SPOT/VEGETATION NDVI 衛(wèi)星遙感數(shù)據(jù)，采用最大值合成法生成，并在ArcGIS中對(duì)其進(jìn)行投影轉(zhuǎn)換、掩膜提取。

1.2.2 土地利用數(shù)據(jù) 土地利用數(shù)據(jù)為來源于中國(guó)科學(xué)院資源環(huán)境與數(shù)據(jù)中心（www.resdc.cn）中國(guó)土地利用遙感監(jiān)測(cè)數(shù)據(jù)，該數(shù)據(jù)集基于Landsat TM/ETM、Landsat 8 遙感影像，通過人工目視解譯生成，每5年為1個(gè)監(jiān)測(cè)期，空間分辨率為1 km。

我國(guó)的土地利用一級(jí)類型主要分為6個(gè)，該分類整體準(zhǔn)確性為92.9%[26-27]，分別是耕地、林地、草地、水域、建設(shè)用地和未利用土地。本研究主要使用一級(jí)分類中的草地作為研究區(qū)草地范圍確定的依據(jù)，選取包含研究時(shí)間段內(nèi)1995、2000、2005、2010、2015和2020年的土地利用數(shù)據(jù)，使用ArcGIS軟件進(jìn)行重分類，篩選出每年度的草地范圍，之后對(duì)各年度草地范圍進(jìn)行疊加，在研究期間，只要某處有草地出現(xiàn)，即作為草地范圍進(jìn)行研究。

1.2.3 氣象數(shù)據(jù) 氣象數(shù)據(jù)主要來源于國(guó)家科技基礎(chǔ)條件平臺(tái)—國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)中心（http://www.geodata.cn），包括逐月降水量、逐月日照時(shí)數(shù)、逐月相對(duì)濕度、逐月平均氣溫空間插值數(shù)據(jù)集，分辨率為1 km，并使用ArcGIS 進(jìn)行投影變換、掩膜、運(yùn)算及分析。

1.2.4 DEM 數(shù)據(jù) 從中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http://www.resdc.cn）下載DEM 數(shù)據(jù)，分辨率為90 m。利用ArcGIS 軟件進(jìn)行拼接、掩膜等預(yù)處理后獲得研究所需數(shù)據(jù)。

1.3 研究方法

1.3.1 像元二分模型 像元二分模型法是間接計(jì)算草地植被覆蓋度的最常見方法，通常使用NDVI直接提取植被覆蓋信息。該模型假設(shè)影像中的像元只由草地植被和非草地植被組成。本研究選擇NDVI 作為參數(shù)，在ENVI 中分別將置信區(qū)間95%和5%的NDVI 值作為NDVIveg和NDVIsoil。其計(jì)算公式如下。

式中，F(xiàn)VCi為第i年的草地植被覆蓋度；i為年份；NDVIsoil為裸地土壤或非草地植被區(qū)域的NDVI 值；NDVIveg為有草地植被完全覆蓋區(qū)域的NDVI 值；NDVI 為像元草地植被覆蓋部分和非草地植被覆蓋部分的加權(quán)平均。

1.3.2 變化率分析 常使用變化率對(duì)影像進(jìn)行逐像元變化率分析，本研究NDVI 的變化率計(jì)算公式如下。

式中，i為年序號(hào)；NDVIi為第i年的NDVI 值；n為樣本量；θslope是NDVI在研究時(shí)間內(nèi)的變化率，其為正值說明草地狀況改善，為負(fù)值則說明草地狀況趨于退化。

1.3.3 相關(guān)性分析 運(yùn)用相關(guān)分析法研究天然草地植被覆蓋度與各氣象因子之間的關(guān)系，用相關(guān)系數(shù)R來具體描述，計(jì)算公式如下。

式中，Rxy為變量x、y的相關(guān)系數(shù)；xi為第i年的草地植被覆蓋度；yi為第i年氣象因子；xˉ為1998—2015 年月平均草地植被覆蓋度；yˉ為1998—2015年各氣象因子的平均值。

2 結(jié)果與分析

2.1 北疆1998—2018 年月平均天然草地植被覆蓋度時(shí)空變化

2.1.1 北疆天然草地植被覆蓋度的空間分布 由像元二分模型計(jì)算出1998—2018 年北疆各月天然草地平均植被覆蓋度之后，求得1998—2018 年北疆天然草地6、7、8、9 月平均植被覆蓋度，結(jié)果（圖2）表明，6 月平均天然草地植被覆蓋度為64%，7 月平均天然草地植被覆蓋度為62.2%，8 月平均天然草地植被覆蓋度為61.5%，9月平均天然草地植被覆蓋度為61.7%?？傮w來看，北疆天然草地植被覆蓋度呈現(xiàn)為中部和東部較低，北部、西部和南部較高的空間分布趨勢(shì)，可能是北部地區(qū)受北冰洋水汽的影響[28]，降水量較充足；西部和南部地帶山頂常年積雪覆蓋，積雪融化使水文條件較好；而中部干旱地帶降水少且蒸發(fā)強(qiáng)烈，存在大范圍的沙漠。

圖2 1998—2018年北疆天然草地月平均植被覆蓋度Fig. 2 Monthly average FVC of natural grassland in Northern Xinjiang from 1998 to 2018

2.1.2 北疆天然草地植被覆蓋度的時(shí)間變化 將1998—2018年長(zhǎng)時(shí)間序列天然草地植被覆蓋度分為1998—2000、2000—2005、2005—2010、2010—2015和2015—2018年5個(gè)時(shí)段，對(duì)研究區(qū)天然草地植被的像元數(shù)值逐年取平均值，計(jì)算得出1998—2018年和5個(gè)時(shí)段北疆天然草地植被覆蓋度的月平均時(shí)間變化折線圖，結(jié)果（圖3）表明，北疆天然草地6 月平均植被覆蓋度呈現(xiàn)波動(dòng)下降的生長(zhǎng)趨勢(shì)，在5 個(gè)時(shí)間段中，1998—2000 年天然草地6 月平均植被覆蓋度較高，達(dá)65.5%；2015—2018年天然草地6月平均植被覆蓋度較低，為62.1%。其中1998—2000、2000—2005、2005—2010 年的6 月平均天然草地植被覆蓋度高于1998—2018 年6 月；2010—2015、2015—2018 年的6 月平均天然草地植被覆蓋度低于1998—2018年6月。

圖3 1998—2018年北疆天然草地分時(shí)段月平均植被覆蓋度Fig. 3 Monthly average FVC of natural grassland in Northern Xinjiang from 1998 to 2018

北疆天然草地7 月平均植被覆蓋度呈現(xiàn)波動(dòng)下降的的生長(zhǎng)趨勢(shì)，在5 個(gè)時(shí)間段中，1998—2000 年天然草地7 月平均植被覆蓋度較高，達(dá)64.3%；2010—2015 年天然草地7 月平均植被覆蓋度較低，為61%。其中1998—2000、2000—2005 年的7 月平均天然草地植被覆蓋度高于1998—2018 年7 月；2005—2010、2010—2015、2015—2018年的7月平均天然草地植被覆蓋度低于1998—2018年7月。

北疆天然草地8 月平均植被覆蓋度呈現(xiàn)波動(dòng)下降的的生長(zhǎng)趨勢(shì)，在5 個(gè)時(shí)間段中，2000—2005 年天然草地8 月平均植被覆蓋度較高，達(dá)63.4%；2015—2018年天然草地8月平均植被覆蓋度較低，為58.4%。其中1998—2000、2000—2005、2005—2010年的8 月平均天然草地植被覆蓋度高于1998—2018年8月；2010—2015、2015—2018年的8月平均天然草地植被覆蓋度低于1998—2018年8月。

北疆天然草地9 月平均植被覆蓋度呈現(xiàn)波動(dòng)下降的的生長(zhǎng)趨勢(shì)，在5 個(gè)時(shí)間段中，1998—2000 年天然草地9 月平均植被覆蓋度較高，達(dá)64.6%；2015—2018年天然草地9月平均植被覆蓋度 較 低，為57.5%。其 中1998—2000、2000—2005、2005—2010 年的9 月平均天然草地植被覆蓋度高于1998—2018年9月；2010—2015、2015—2018年的9月平均天然草地植被覆蓋度低于1998—2018年9月。

2.2 北疆天然草地NDVI變化率

圖4 是北疆1998—2018 年月平均天然草地NDVI 變化率空間分布。從空間上看，6 月研究區(qū)總面積45.7%天然草地的NDVI 呈上升趨勢(shì)，54.3%呈下降趨勢(shì)；7 月，55.9%天然草地NDVI 呈上升趨勢(shì)，44.1%呈下降趨勢(shì)；8 月，33.4%天然草地NDVI 呈上升趨勢(shì)，66.6%呈下降趨勢(shì)；9 月，29.2%天然草地NDVI 呈上升趨勢(shì)，70.8%呈下降趨勢(shì)?？傮w來看，1998—2018 年的6、8、9 月北疆天然草地NDVI 呈增加趨勢(shì)的區(qū)域均小于呈減少趨勢(shì)的區(qū)域，且減少趨勢(shì)的區(qū)域變化率多處于-0.01～0.00 之間；而1998—2018 年7 月北疆天然草地NDVI 呈增加趨勢(shì)的區(qū)域大于呈減少趨勢(shì)的區(qū)域，且呈增加趨勢(shì)區(qū)域的變化率多處于0.00～0.03之間。

圖4 1998—2018年北疆天然草地NDVI月平均變化率Fig. 4 Monthly average change rate of NDVI of natural grassland in Northern Xinjiang from 1998 to 2018

2.3 北疆氣象因子的空間變化特征

北疆多年月均氣溫有明顯的空間分異特征，區(qū)域內(nèi)氣溫空間差異較大。如圖5 所示，1998—2015年的6月月均氣溫為10～25 ℃，平均19.77 ℃；7月為12～26 ℃，平均21.28 ℃；8月為12～25 ℃，平均19.87 ℃；9月為7～18 ℃，平均13.83 ℃。氣溫受海拔影響明顯，準(zhǔn)噶爾盆地北部的阿爾泰山和南部的天山山脈、西部的塔爾巴哈臺(tái)山和阿拉套山年均氣溫較低，而盆地西南部的艾比湖為氣溫高值中心。

圖5 1998—2015年北疆氣象因子空間變化Fig. 5 Spatial variation of meteorological factors in Northern Xinjiang from 1998 to 2015

北疆多年月均降水有明顯的空間分異特征，整體上呈現(xiàn)南北高中間低，西部邊緣高、中部和東部低，自北西南三面向中部遞減的趨勢(shì)。1998—2015 年6、7、8 和9 月北疆月均降水量分別為9～107、15～125、10～116和5～69 mm。北部主要受大西洋和北冰洋的水汽影響；南部的天山山脈和博羅科努山、西部的塔爾巴哈臺(tái)山和阿拉套山整體海拔高，受垂直地帶性分布的影響，加之山頂多年冰雪覆蓋，多年平均月降水量較高；中部和東部以古爾班通古特沙漠為中心三面環(huán)山且遠(yuǎn)離海洋，水汽難以到達(dá)，因此多年平均月降水量較低。

北疆多年月均日照時(shí)數(shù)具有明顯的空間分異特征，1998—2015年6、7、8和9月北疆多年月均日照時(shí)數(shù)分別為250～335、263～340、259～327 和238～286 h?？傮w來看，1998—2015年6、7、8月的月均日照時(shí)數(shù)呈現(xiàn)北高南低，自北向南遞減趨勢(shì)；9月呈現(xiàn)東南高西南低，自東南向西南遞減趨勢(shì)。

北疆多年月均相對(duì)濕度具有明顯的空間分異特征，1998—2015 年6、7、8 和9 月月均相對(duì)濕度分別為35%～65%、36%～67%、35%～66%和35%～67%。總體來看，北疆多年月均相對(duì)濕度呈現(xiàn)出中部和東南低北部和西南高，自中部和東南向四周遞增趨勢(shì)。

2.4 北疆天然草地植被覆蓋度與氣象因子的關(guān)系

2.4.1 天然草地植被覆蓋度與氣溫的關(guān)系 如圖6 所示，1998—2015 年6 月占研究區(qū)總面積30.1%的天然草地植被覆蓋度與平均氣溫呈正相關(guān)，69.9%的天然草地植被覆蓋度與平均氣溫呈負(fù)相關(guān)；7月，18.4%的天然草地植被覆蓋度與平均氣溫呈正相關(guān)，81.6%的天然草地植被覆蓋度與平均氣溫呈負(fù)相關(guān)；8月，33.1%的天然草地植被覆蓋度與平均氣溫呈正相關(guān)，66.9%的天然草地植被覆蓋度與平均氣溫呈負(fù)相關(guān)；9月，77.1%的天然草地植被覆蓋度與平均氣溫呈正相關(guān)，22.9%的天然草地植被覆蓋度與平均氣溫呈負(fù)相關(guān)。整體來看，1998—2015 年6、7 和8 月平均氣溫與天然草地植被覆蓋度呈正相關(guān)的區(qū)域小于呈負(fù)相關(guān)的區(qū)域；9月平均氣溫與天然草地植被覆蓋度呈正相關(guān)的區(qū)域大于呈負(fù)相關(guān)的區(qū)域。

圖6 1998—2015年北疆天然草地植被覆蓋度與氣溫的相關(guān)分析Fig. 6 Correlation analysis of natural grassland vegetation coverage and temperature in Northern Xinjiang from 1998 to 2015

2.4.2 天然草地植被覆蓋度與降水的關(guān)系 如圖7 所示，1998—2015 年6 月占研究區(qū)總面積65.7%的天然草地植被覆蓋度與降水呈正相關(guān)，34.3%的天然草地植被覆蓋度與降水呈負(fù)相關(guān)；7 月，74.6%的天然草地植被覆蓋度與降水呈正相關(guān)，25.4%的天然草地植被覆蓋度與降水呈負(fù)相關(guān)；8月，75.4%的天然草地植被覆蓋度與降水呈正相關(guān)，24.6%的天然草地植被覆蓋度與降水呈負(fù)相關(guān)；9月，30.3%的天然草地植被覆蓋度與降水呈正相關(guān)，69.7%的天然草地植被覆蓋度與降水呈負(fù)相關(guān)。整體上來看，1998—2015年6、7和8月的月降水與天然草地植被覆蓋度呈正相關(guān)的區(qū)域面積大于呈負(fù)相關(guān)的區(qū)域；9月降水與天然草地植被覆蓋度呈正相關(guān)的區(qū)域小于呈負(fù)相關(guān)的區(qū)域。

圖7 1998—2015年北疆天然草地植被覆蓋度與降水的相關(guān)分析Fig. 7 Correlation analysis of natural grassland vegetation coverage and precipitation in Northern Xinjiang from 1998 to 2015

2.4.3 天然草地植被覆蓋度與日照時(shí)數(shù)的關(guān)系如圖8 所示，1998—2015 年6 月占研究區(qū)總面積59.7%的天然草地植被覆蓋度與日照時(shí)數(shù)呈正相關(guān)，40.3%的天然草地植被覆蓋度與日照時(shí)數(shù)呈負(fù)相關(guān)；7 月，24.7%的天然草地植被覆蓋度與日照時(shí)數(shù)呈正相關(guān)，75.3%的天然草地植被覆蓋度與日照時(shí)數(shù)呈負(fù)相關(guān)；8 月，34.4%的天然草地植被覆蓋度與日照時(shí)數(shù)呈正相關(guān)，65.6%的天然草地植被覆蓋度與日照時(shí)數(shù)呈負(fù)相關(guān)；9 月，80.8%的天然草地植被覆蓋度與日照時(shí)數(shù)呈正相關(guān)，19.2%的天然草地植被覆蓋度與日照時(shí)數(shù)呈負(fù)相關(guān)。整體上來看，1998—2015 年6 和9 月月均日照時(shí)數(shù)與天然草地植被覆蓋度呈正相關(guān)的區(qū)域大于呈負(fù)相關(guān)的區(qū)域，7和8月月均日照時(shí)數(shù)與天然草地植被覆蓋度呈正相關(guān)的區(qū)域小于呈負(fù)相關(guān)的區(qū)域。

圖8 1998—2015年北疆天然草地植被覆蓋度與月日照時(shí)數(shù)的相關(guān)分析Fig. 8 Correlation analysis of natural grassland vegetation coverage and monthly sunshine hours in Northern Xinjiang from 1998 to 2015

2.4.4 天然草地植被覆蓋度與相對(duì)濕度的關(guān)系如圖9 所示，1998—2015 年6 月占研究區(qū)總面積70.8%的天然草地植被覆蓋度與相對(duì)濕度呈正相關(guān)，29.2%的天然草地植被覆蓋度與相對(duì)濕度呈負(fù)相關(guān)；7 月，80.6%的天然草地植被覆蓋度與相對(duì)濕度呈正相關(guān)，19.4%的天然草地植被覆蓋度與相對(duì)濕度呈負(fù)相關(guān)；8 月，81.0%天然草地植被覆蓋度與相對(duì)濕度呈正相關(guān)，19.0%的天然草地植被覆蓋度與相對(duì)濕度呈負(fù)相關(guān)；9 月，49.9%天然草地植被覆蓋度與相對(duì)濕度呈正相關(guān)，50.1%的天然草地植被覆蓋度與相對(duì)濕度呈負(fù)相關(guān)。整體上來看，1998—2015 年6、7 和8 月相對(duì)濕度與天然草地植被覆蓋度呈正相關(guān)的區(qū)域遠(yuǎn)大于呈負(fù)相關(guān)的區(qū)域；9 月相對(duì)濕度與天然草地植被覆蓋度呈正相關(guān)的區(qū)域略微小于呈負(fù)相關(guān)的區(qū)域。

圖9 1998—2015年北疆天然草地植被覆蓋度與月相對(duì)濕度的相關(guān)分析Fig. 9 Correlation analysis of natural grassland vegetation coverage and monthly relative humidity in Northern Xinjiang from 1998 to 2015

3 討論

研究表明，氣候變暖是導(dǎo)致北半球中高緯度地區(qū)植被覆蓋度變化的主要因素[29-31]。本研究表明，北疆天然草地植被覆蓋度的變化主要受降水和相對(duì)濕度變化的影響，這可能與北疆地處干旱、半干旱地區(qū)有關(guān)[32]；氣溫、降水、日照時(shí)數(shù)、相對(duì)濕度均為北疆天然草地植被覆蓋度的重要影響因子，且多表現(xiàn)為促進(jìn)作用，但相比較來說，降水和相對(duì)濕度對(duì)北疆天然草地植被覆蓋度的促進(jìn)作用更加明確。從時(shí)間上來看，北疆天然草地植被覆蓋度與降水量、相對(duì)濕度等氣象因子關(guān)系密切，說明在干旱、半干旱區(qū)水分條件是影響天然草地植被覆蓋度變化的重要因素，水分增多有利于草地植被生長(zhǎng)，與周丹等[33]研究結(jié)果一致，表明降水量對(duì)草地植被覆蓋度有重要影響。另外，氣溫和日照時(shí)數(shù)也是影響北疆天然草地植被覆蓋度重要?dú)庀笠蜃?，這在9 月表現(xiàn)得最明顯，可能與9 月氣溫升高減緩或推遲了草地植被的枯黃期有關(guān)[34]。

總體來看，北疆天然草地在6、7、8、9 月的平均植被覆蓋度受降水、氣溫、日照時(shí)數(shù)、相對(duì)濕度和海拔的影響，空間上呈現(xiàn)降水較少、氣溫較高的中部和東部天然草地植被覆蓋度較低，而海拔較高、降水較多且相對(duì)濕度較高的北部、西部和南部天然草地植被覆蓋度較高，天然草地植被覆蓋度總體上由中部和東部向北部、西部和南部三面波動(dòng)增加。北部地區(qū)受北冰洋水汽的影響[29]，其降水量較充足，西部和南部地帶山頂常年積雪覆蓋，積雪融化水文條件較好；中部干旱地帶降水少且蒸發(fā)強(qiáng)烈，存在大范圍的沙漠。21 年間北疆6、7、8、9 月的平均天然草地植被覆蓋度呈不顯著波動(dòng)降低的趨勢(shì)，其中月均天然草地覆蓋度表現(xiàn)為6 月＞7 月＞9 月＞8 月。本研究尚欠缺從多角度分析影響天然草地植被覆蓋度的因子，對(duì)于利用多種植被指數(shù)相結(jié)合進(jìn)行反演天然草地植被覆蓋度和從多種方面分析不同地形和人文因子等對(duì)天然草地植被覆蓋度的影響仍有待進(jìn)一步研究。