李紅梅

(青海省氣候中心，青海 西寧 810001)

植被是陸地生態(tài)系統(tǒng)中對(duì)氣候變化響應(yīng)最敏感的組分，在一定程度上是氣候變化的指示器。在當(dāng)前氣候變化背景下，世界各地植被生長發(fā)育狀況和地域分布特征都發(fā)生了明顯改變，但隨著各地氣候變化特征的不同，植被所作出的響應(yīng)亦存在較大的差異[1-6]。有學(xué)者研究表明20 世紀(jì)80 年代初至90 年代末，氣候變暖導(dǎo)致了中國青藏高原和北半球高緯度地區(qū)植被活動(dòng)的顯著增強(qiáng)[7-10]。也有學(xué)者指出2000 年以后北半球多個(gè)地區(qū)植被活動(dòng)呈減弱趨勢，氣溫持續(xù)升高所導(dǎo)致的干旱化過程可能是驅(qū)動(dòng)這些地區(qū)植被活動(dòng)變化的主要原因[11]。綜合以上分析可以發(fā)現(xiàn)，氣候變化在不同時(shí)間尺度和空間尺度上分析結(jié)果存在較大的差異，或者由于利用的資料源不同，所分析的結(jié)果也存在很大的不同，由此引起的對(duì)植被的影響評(píng)估結(jié)果也存在很大的差異性，尤其是在地形條件復(fù)雜多樣的青藏高原表現(xiàn)更為突出。

柴達(dá)木盆地位于青海高原西北部，是我國海拔最高的高原型盆地，盆地四周高山環(huán)抱，地貌復(fù)雜多樣，垂直分異明顯，整個(gè)盆地大致以海拔3350 m為界，可分為盆地內(nèi)干旱荒漠區(qū)和盆地四周高寒區(qū)[12]。該區(qū)域降水稀少, 蒸發(fā)強(qiáng)烈，氣候干燥，生態(tài)環(huán)境脆弱，是整個(gè)青藏高原升溫最顯著的地方，也是全球氣候變化影響最敏感的地區(qū)[13-15]。近年來隨著全球氣候變暖，眾多學(xué)者對(duì)柴達(dá)木盆地氣候變化特征進(jìn)行了研究[16-22]，并指出，隨著柴達(dá)木盆地降水的增多，該區(qū)正經(jīng)歷著由暖干向暖濕的氣候轉(zhuǎn)型[23-26]。然而，在此氣候轉(zhuǎn)型背景下，對(duì)本地條件極其脆弱的草地植被產(chǎn)生了何種影響，加之，柴達(dá)木盆地各地氣候條件差異較大，氣候變化對(duì)不同區(qū)域、不同植被類型所產(chǎn)生的影響是否相同，目前對(duì)這些問題的研究尚顯不足。 本研究利用1961-2016年氣象觀測資料和1982-2016年遙感監(jiān)測資料，評(píng)估柴達(dá)木盆地氣候變化對(duì)不同區(qū)域、不同植被類型歸一化植被指數(shù)(normalized difference vegetation index, NDVI)的影響，以及氣候變化對(duì)植被演替的影響趨勢，通過研究氣候變化對(duì)植被的影響，在一定程度上可以預(yù)估未來植被可能變化的趨勢，為今后柴達(dá)木盆地實(shí)施生態(tài)環(huán)境保護(hù)和開展適應(yīng)氣候變化工作等提供參考。

1 材料與方法

1.1 資料來源

氣象數(shù)據(jù)利用柴達(dá)木盆地1961-2016年茫崖、冷湖、大柴旦、德令哈、天峻、格爾木、諾木洪、都蘭、茶卡和小灶火共10個(gè)氣象站逐日氣溫、降水等資料，所有資料均經(jīng)過嚴(yán)格的質(zhì)量控制，具有較高的可信度。

NDVI值采用GIMMS和MODIS兩個(gè)衛(wèi)星源的數(shù)據(jù)。GIMMS 數(shù)據(jù)是美國航空航天局(National Aeronautics and Space Administration, NASA)全球監(jiān)測與模型研究組(Global Inventor Modeling and Mapping Studies, GIMMS)發(fā)布的NDVI半月最大值合成值(maximum value composites, MVC), 空間分辨率為 8 km×8 km, 時(shí)段為1982年5月至2006年12月。MODIS數(shù)據(jù)采用青海省遙感中心 EOS/MODIS系統(tǒng)接收的數(shù)據(jù), 空間分辨率為 250 m×250 m，時(shí)段為2002年1月至2016年12月。兩種NDVI數(shù)據(jù)集都已經(jīng)過幾何精糾正、輻射校正、大氣校正等預(yù)處理, 并采用最大值合成法減少云、大氣、太陽高度角等的影響。

1.2 方法

1.2.1遙感資料的處理 由于GIMMS和MODIS數(shù)據(jù)采用了不同的傳感器，利用1982-2006年GIMMS資料和2001-2016年MODIS資料相重疊的2001-2006年NDVI值進(jìn)行相關(guān)性分析。分析發(fā)現(xiàn)在不同區(qū)域(盆地四周高寒區(qū)和盆地內(nèi)干旱荒漠區(qū))以及不同植被類型(溫性荒漠、低地草甸、高寒草原)的GIMMS 數(shù)據(jù)和NDVI數(shù)據(jù)均存在較高的相關(guān)性，都通過了置信度為0.01的極顯著檢驗(yàn)(圖1)。根據(jù)GIMMS和MODIS兩者關(guān)系，利用2007-2016年GIMMS的NDVI值插補(bǔ)MODIS NDVI值，建立1982-2016年柴達(dá)木盆地四周高寒區(qū)、盆地內(nèi)干旱荒漠區(qū)、低地草甸、溫性荒漠、高寒草原植被MODIS NDVI值序列。

1.2.2遙感資料的插補(bǔ)模型 根據(jù)不同區(qū)域(盆地四周高寒區(qū)和盆地內(nèi)干旱荒漠區(qū))和不同植被類型(低地草甸類、溫性荒漠類、高寒草原類)GIMMS NDVI和MODIS NDVI值的相關(guān)性，利用GIMMS NDVI值對(duì)MODIS NDVI值進(jìn)行插補(bǔ)，具體插補(bǔ)模型見表1。

圖1 2001-2006年盆地四周高寒區(qū)(a)、盆地內(nèi)干旱荒漠區(qū)(b)、低地草甸類(c)、溫性荒漠類(d)和高寒草原類(e)GIMMS NDVI和MODIS NDVI相關(guān)性Fig.1 The correlation between GIMMS NDVI and MODIS NDVI in high cold region around Qaidam Basin (a), arid desert region in the Qaidam Basin (b), lowland meadow region (c), warm desert region (d) and alpine meadow region (e)

類型Type模型ModelF值FvalueR2盆地四周高寒區(qū)ThecoldregionaroundbasinY=1．547X+0．184935．063??0．9524盆地內(nèi)干旱荒漠區(qū)TheariddesertregioninbasinY=2．717X-0．014330．426??0．9019低地草甸類LowlandmeadowY=1．181X+0．048643．673??0．9758溫性荒漠類WarmdesertY=0．493X+0．0204110．309??0．9154高寒草原類AlpinemeadowY=0．414X+0．021432．517??0．9697

注：Y為MODIS NDVI，X為GIMMS NDVI，**表示通過置信度為0.01的極顯著檢驗(yàn)。

Note:Yexpress MODIS NDVI,Xexpress GIMMS NDVI, **express through highly significant test of the reliability of 0.01.

1.2.3植被演替分析方法 采用改進(jìn)的綜合順序分類法分析植被演替特征，綜合順序分類法是以任繼周院士為代表，在其草原發(fā)生與發(fā)展理論的指導(dǎo)下，參考并吸收世界各國草原分類方法的優(yōu)點(diǎn)，提出的一種草原分類方法[27-28]。1995年胡自治等[29]對(duì)綜合順序分類法進(jìn)行了新的改進(jìn)，使劃分結(jié)果更趨完善，該方法以>0 ℃的年積溫和濕潤度指數(shù)K作為草地類型的劃分因子，比較適合于劃分青海高原的草地類型[30-32]。

2 結(jié)果與分析

2.1 植被NDVI變化特征分析

2.1.1不同區(qū)域植被變化趨勢 1982-2016年柴達(dá)木盆地四周高寒區(qū)植被NDVI值呈明顯增加趨勢，趨勢系數(shù)為0.008/10年(圖2a)，并通過顯著性水平0.05的檢驗(yàn)。從長期變化曲線來看以1999年為界大致分為兩個(gè)階段，其中1982-1998年NDVI值總體呈略微下降趨勢，而1999年以來植被NDVI值變化較為平穩(wěn)，且維持在一個(gè)較高水平。

1982-2016年柴達(dá)木盆地內(nèi)干旱荒漠區(qū)植被總體呈增加趨勢，趨勢系數(shù)為0.007/10年，通過顯著性水平0.01的檢驗(yàn)。干旱荒漠區(qū)NDVI值具有階段性變化，1996年以前年際間變化幅度較大，呈增加趨勢，1996年以后變化較為平穩(wěn)(圖2b)。

圖2 1982-2016年柴達(dá)木盆地四周高寒區(qū)(a)、盆地內(nèi)干旱荒漠區(qū)(b)植被NDVI變化趨勢Fig.2 The change trend of NDVI in the high cold region around Qaidam Basin (a) and arid desert region in the Qaidam Basin (b) from 1982 to 2016

2.1.2不同植被類型NDVI變化特征 低地草甸類、溫性荒漠類、高寒草原類3種植被類型是柴達(dá)木盆地主要植被類型，1982-2016年這3種植被類型NDVI值均呈增加趨勢，趨勢系數(shù)分別為0.001/10年、0.003/10年和0.005/10年，其中只有高寒草原類變化趨勢通過置信度0.01的極顯著水平檢驗(yàn)。主要是因?yàn)楦吆菰悘V泛分布在柴達(dá)木盆地四周高寒區(qū)，受氣溫升高、降水增多的有利影響，NDVI值增加較為明顯，而在盆地中分布較多的溫性荒漠類受氣候干旱的影響，NDVI值變化不明顯(圖3)。

2.2 氣候變化對(duì)植被NDVI的影響

2.2.1氣候條件對(duì)不同類型植被NDVI的影響 以柴達(dá)木盆地低地草甸類、溫性荒漠類、高寒草原類為代表性植被，分析氣候變化對(duì)這3種植被NDVI的影響。受柴達(dá)木盆地植被生長季高溫少雨的影響，低地草甸類、溫性荒漠類、高寒草原類3種主要植被類型生長均受降水和蒸發(fā)量影響較大，其中低地草甸類、高寒草原類與降水量呈線性相關(guān)，隨著降水量的增加，植被NDVI值呈增加趨勢(圖4a，c)。而降水對(duì)溫性荒漠類植物的影響較為復(fù)雜，在降水量開始增加時(shí)，隨著降水量的增加NDVI值呈增加趨勢，降水增加到一定程度后NDVI值不再增加(圖4b)。

蒸發(fā)量對(duì)植被NDVI值具有較明顯的負(fù)影響(圖4d～f)，通過相關(guān)分析表明，低地草甸類、溫性荒漠類與蒸發(fā)量的相關(guān)系數(shù)分別為-0.49和-0.56，通過置信度為0.01的顯著性檢驗(yàn)，高寒草原NDVI與蒸發(fā)量的相關(guān)系數(shù)為0.40，通過了置信度為0.05的顯著性檢驗(yàn)。

2.2.2植被NDVI的氣候模擬模型 根據(jù)1982-2016年低地草甸類、溫性荒漠類和高寒草原類NDVI值和4-8月降水量和蒸發(fā)量變化特征，建立不同植被類型氣候模擬模型，見表2。

利用1982-2016年4-8月降水量和蒸發(fā)量模擬低地草甸類、溫性荒漠類和高寒草原類NDVI值，從圖5中可以看出，3種類型草地NDVI值的氣候模擬模型都能很好地模擬出歷年變化趨勢。

圖3 1982-2016年柴達(dá)木盆地低地草甸類(a)、溫性荒漠類(b)和高寒草原類(c)NDVI變化趨勢Fig.3 The NDVI change trend of lowland meadow (a), warm desert (b) and alpine meadow (c) from 1982 to 2016

圖4 低地草甸類NDVI、溫性荒漠類NDVI、高寒草原類NDVI與4-8月降水量(a、b、c)以及年蒸發(fā)量(d、e、f)相關(guān)性Fig.4 The correlation between NDVI of lowland meadow (a), warm desert (b) and alpine meadow (c) and the sum of precipitation from April to August, and the correlation between NDVI of lowland meadow (d), warm desert (e) and alpine meadow (f) and the annual evaporation

2.3 植被類型演替特征

利用綜合順序分類法分析柴達(dá)木盆地1961-2016年植被演替特征，根據(jù)相關(guān)學(xué)者對(duì)柴達(dá)木盆地氣溫變化特征的研究發(fā)現(xiàn)，該地區(qū)在1987年發(fā)生由冷向暖的突變[22]，因此以氣溫發(fā)生突變的1987年為界和以10年為間隔，計(jì)算不同年代植被演替方向。>0 ℃年積溫和濕潤度K值是確定植被類型的氣象因子，本研究主要基于這兩個(gè)區(qū)劃因子進(jìn)行植被演替特征的分析。

表2 低地草甸類、溫性荒漠類和高寒草原類NDVI氣候模擬模型Table 2 The climate simulation model of NDVI of lowland meadow, warm desert and alpine meadow

注：Y為NDVI值，X1為4-8月降水量，X2為4-8月蒸發(fā)量，**表示通過置信度為0.01的極顯著檢驗(yàn)。

Note:Yexpress NDVI,X1express precipitation from April to August,X2expresse evaporation from April to August， ** express through highly significant test of the reliability of 0.01.

圖5 低地草甸類(a)、溫性荒漠類(b)和高寒草原類(c)實(shí)際NDVI值與模擬值對(duì)比Fig.5 Comparison of actual and simulated NDVI values of lowland meadow (a), warm desert (b) and alpine meadow (c)

2.3.1氣象因子變化特征 1961-2016年>0 ℃年積溫呈顯著上升趨勢，趨勢系數(shù)達(dá)每10年上升85.4 ℃，與柴達(dá)木盆地年平均氣溫變化趨勢一致，1998年以來>0 ℃年積溫呈明顯增加趨勢，1961-1997年平均>0 ℃年積溫為2157.6 ℃，而1998-2016年平均值上升為2452.4 ℃，前后兩個(gè)時(shí)段平均值相差294.8 ℃(圖6a)。同時(shí)，從>0 ℃年積溫變化圖也可以看出，在1994年左右柴達(dá)木盆地的熱量級(jí)由寒溫帶轉(zhuǎn)為微溫帶。

從空間變化趨勢來看，柴達(dá)木盆地的西部積溫增加率大于東部地區(qū)，其中茫崖、格爾木一帶>0 ℃年積溫增加幅度較大，平均每10年達(dá)到183.8 ℃，而柴達(dá)木盆地的東部天峻、烏蘭等地增加率相對(duì)較小，平均每10年增加34.1 ℃(圖6b)。

受柴達(dá)木盆地降水量增多的影響，1961-2015年濕潤度K值呈微弱增加趨勢，整體來看，90年代以前濕潤度K值變化幅度較大，位于歷史前4高的極值均出現(xiàn)在這一時(shí)段；1990年以后濕潤度K值變化平穩(wěn)，尤其是進(jìn)入21世紀(jì)以來增加明顯，且一直在較高水平(圖7a)。

柴達(dá)木盆地各地濕潤度指數(shù)變化趨勢不盡相同，呈明顯的東高西低特征，東部地區(qū)的德令哈、烏蘭、都蘭為一個(gè)大值區(qū)，濕潤度指數(shù)增加系數(shù)為每10年增加0.052，而柴達(dá)木盆地西部的格爾木一帶濕潤度指數(shù)呈降低趨勢，平均每10年降低0.028(圖7b)。

圖6 1961-2016年柴達(dá)木盆地>0 ℃年積溫變化(a)、空間變率分布(b)Fig.6 >0 ℃ accumulated temperature change (a), and its variability spatial distribution (b) from 1961 to 2016 in Qaidam Basin

圖7 1961-2015年柴達(dá)木盆地濕潤度K變化(a)、濕潤度K空間變率分布(b)Fig.7 K index change (a), and its variability spatial distribution (b) from 1961 to 2015 in Qaidam Basin

2.3.2以氣溫突變點(diǎn)為界分析植被演替 柴達(dá)木盆地氣溫在1987年前后發(fā)生了突變，因此以1987年為界分析氣候突變前后植被演替特征。

在氣象要素進(jìn)行空間插值時(shí)，發(fā)現(xiàn)氣候突變前后的>0 ℃年積溫和濕潤度K值與海拔存在較高的相關(guān)性，為了提高空間插值精度，采用表3所示的回歸方程進(jìn)行插值。

利用表3中>0 ℃年積溫和濕潤度K值的插值方程對(duì)柴達(dá)木盆地各格點(diǎn)進(jìn)行插值，并根據(jù)綜合順序分類法對(duì)氣候突變前和突變后植被進(jìn)行分類。從表4可以看出，氣候突變植被類型發(fā)生了變化，突變前植被類型主要有寒冷潮濕類、寒冷濕潤類、寒溫潮濕類、冷溫干旱類、冷溫微干類、冷溫微潤類、冷溫濕潤類、微溫極干類和微溫干旱類，而氣候突變后植被類型主要有寒冷潮濕類、寒溫潮濕類、冷溫干旱類、冷溫微干類、冷溫微潤類、冷溫濕潤類、冷溫潮濕類、微溫極干類和微溫干旱類。氣溫突變前后雖然植被類型數(shù)量都為9類，但突變前的寒冷濕潤植被類型消失，氣溫突變后新增了冷溫潮濕類。

由植被類型發(fā)生改變的區(qū)域可以看出，在盆地四周高寒區(qū)和盆地內(nèi)干旱荒漠區(qū)的交錯(cuò)區(qū)氣候變化對(duì)植被的影響最大，植被類型易發(fā)生改變。

從各植被類型面積變化(表5)來看，氣溫突變前后熱量級(jí)發(fā)生了顯著改變，1987年前植被類型以寒溫類為主，而1987年以后轉(zhuǎn)為以微溫類植被為主。從濕潤度面積變化結(jié)果來看，植被主要朝著濕潤的方向發(fā)展，尤以極干類植被變化較為顯著，由1987年以前的12.85×104m2減少為之后的9.67×104m2，減少了近25%。

表3 >0 ℃年積溫、濕潤度K回歸方程Table 3 Regression equation of >0 ℃ accumulated temperature and K index

注：表中HB表示海拔，**表示通過顯著性水平0.01的相關(guān)性檢驗(yàn)。

Note: HB express altitude, **express through highly significant test of the reliability of 0.01.

表4 氣溫突變前后主要植被類型Table 4 Main vegetation types before and after temperature change

表5 氣溫突變前后不同類型植被分布面積變化Table 5 The area of different vegetation types of before and after temperature change (×104 m2)

2.3.3以年代際變化為界分析植被演替 >0 ℃年積溫和濕潤度K是植被類型分布的決定性因素，受這兩要素的影響，50多年來，植被類型發(fā)生了變化。將1961-2010年以10年為間隔進(jìn)行植被類型演變分析(表6)，從表中可以看出柴達(dá)木盆地寒溫帶植被類型逐步減少，微溫帶植被類型逐步增多，同時(shí)濕潤度級(jí)別逐漸升高，尤其是進(jìn)入21世紀(jì)以來，出現(xiàn)了微溫微干和微溫微潤的植被類型。

表6 不同年代不同類型植被分布面積變化Table 6 The area of different vegetation types in different interdecadal (×104 m2)

3 討論與結(jié)論

溫性荒漠類植被對(duì)降水變化響應(yīng)的特征主要是因?yàn)闇匦曰哪愔饕猿瞪嗄竞桶牍嗄緸閮?yōu)勢種，這些植被長期生長于干旱的環(huán)境，比較適合于旱生生境，適度的干旱有利于其生長發(fā)育，當(dāng)降水量增大到一定程度后，植被類型就會(huì)演替到合適潮濕環(huán)境的植被生長，這和相關(guān)的研究結(jié)論比較一致[33]。

通過對(duì)柴達(dá)木盆地不同區(qū)域植被NDVI指數(shù)和不同類型植被NDVI指數(shù)變化特征、不同時(shí)期植被演替特征分析，主要得出以下結(jié)論：

1)按不同海拔、植被分布特征等將柴達(dá)木盆地分為盆地四周高寒區(qū)和盆地內(nèi)干旱荒漠區(qū)。1982-2016年兩區(qū)域植被NDVI指數(shù)均呈顯著上升趨勢，趨勢系數(shù)分別為0.009/10年、0.007/10年，其中四周高寒區(qū)植被NDVI變化趨勢通過顯著性水平0.05的檢驗(yàn)，而盆地內(nèi)干旱荒漠植被NDVI指數(shù)變化趨勢通過顯著性水平0.01的檢驗(yàn)。

2)從不同植被類型NDVI指數(shù)變化特征可以看出，1982-2016年溫性荒漠類和低地草甸類植被類型NDVI呈微弱的增加趨勢，但沒有通過顯著性檢驗(yàn)，而高寒草原類植被NDVI指數(shù)變化明顯，趨勢系數(shù)為0.0058/10年，通過顯著性水平0.01的檢驗(yàn)。

3)柴達(dá)木盆地氣候極其干旱，年降水量和蒸發(fā)量是影響植被生長發(fā)育的重要?dú)庀笠蜃?。因此年降水量與植被NDVI值呈明顯的正相關(guān)關(guān)系，而蒸發(fā)量與植被NDVI值呈顯著的負(fù)相關(guān)關(guān)系。但值得探討的是溫性荒漠類植被NDVI指數(shù)與降水量呈非線性相關(guān)，當(dāng)降水量增加到一定程度后，植被NDVI指數(shù)不再增加，這可能是該植被類型比較適合干旱環(huán)境，當(dāng)生存環(huán)境發(fā)生較大的改變時(shí)，則會(huì)抑制生長發(fā)育，被其他類型植被所代替。

4)柴達(dá)木盆地在1987年發(fā)生氣溫由冷向暖的突變，以氣溫發(fā)生突變的1987年為界，1987年前后植被類型雖然均為9類，但突變后植被朝著溫暖化及濕潤化的方向發(fā)展。從植被類型發(fā)生變化的區(qū)域分布來看，在盆地四周高寒區(qū)和盆地內(nèi)干旱荒漠區(qū)的交錯(cuò)區(qū)氣候變化對(duì)植被演替的影響最為明顯。

5)從年代際變化特征分析結(jié)果表明，寒溫帶植被類型逐漸減少，微溫帶植被類型逐漸增多，同時(shí)濕潤度級(jí)別有所升高，尤其是進(jìn)入21世紀(jì)以來，出現(xiàn)了微溫微干和微溫微潤的植被類型。

References：

[1] Zhang X Z, Zheng J Y, He F N,etal. Spatial differences of variations in autumn fractional vegetation coverage in eastern China during 1982-2006. Journal of Natural Resources, 2013, 28(1): 28-37.

張學(xué)珍, 鄭景云, 何凡能, 等. 1982-2006年中國東部秋季植被覆蓋變化過程的區(qū)域差異. 自然資源學(xué)報(bào), 2013, 28(1): 28-37.

[2] Mu S J, Li J L, Chen Y Z,etal. Spatial differences of variations of vegetation coverage in Inner Mongolia during 2001-2010. Acta Geographica Sinica, 2012, 67(9): 1255-1268.

穆少杰, 李建龍, 陳奕兆, 等. 2001-2010年內(nèi)蒙古植被覆蓋度時(shí)空變化特征. 地理學(xué)報(bào), 2012, 67(9): 1255-1268.

[3] Ren Z C, Zhu H Z, Zhang D G,etal. Temporal and spatial response of vegetation net primary productivity to climate change in Buryatiya Republic Russia. Journal of Natural Resources, 2011, 26(5): 790-801.

任正超, 朱華忠, 張德罡, 等. 俄羅斯布里亞特共和國植被NPP對(duì)氣候變化的時(shí)空響應(yīng). 自然資源學(xué)報(bào), 2011, 26(5): 790-801.

[4] Tucker C J, Newcomb W W, Los S O,etal. Mean and inter-year variation of growing seasons normalized difference vegetation index for the Sahel 1981-1989. International Journal of Remote Sensing, 1991, 12: 1133-1135.

[5] Yang L M, Bruce K W, Larry L T,etal. An analysis of relationships among climate forcing and time-integrated NDVI of grasslands over the U.S. Northern and Central Great Plains. Remote Sensing of Environment, 1998, 65(1): 25-37.

[6] Yang W, Yang Z, Merchant J W. An assessment of AVHRR NDVI eco-climatological relations in Nebraska, USA. International Journal of Remote Sensing, 1997, 18(10): 2161-2180.

[7] Myneni R B, Keeling C D, Tucker C J,etal. Increased plant growth in the northern high latitudes from 1981 to 1991. Nature, 1997, 386: 698-702.

[8] Zhou L M, Tucker C J, Kaufmann R K,etal. Variations in northern vegetation activity inferred from satellite data of vegetation index during 1981 to 1999. Journal of Geophysical Research, 2001, 106: 20069-20083.

[9] Xu X K, Chen H, Jason L K. Temporal and spatial variation of vegetation cover and its cause analysis in the Tibetan Plateau under climate warming. Science in China, 2008, 53(4): 456-462.

徐興奎, 陳紅, Jason L K. 氣候變暖背景下青藏高原植被覆蓋特征的時(shí)空變化及其成因分析. 科學(xué)通報(bào), 2008, 53(4): 456-462.

[10] Zhang G L, Ouyang H, Zhang X Z,etal. Vegetation change and its responses to climatic variation based on eco-geographical regions of Tibetan Plateau. Geographical Research, 2010, 29(11): 2004-2016.

張戈麗, 歐陽華, 張憲洲, 等. 基于生態(tài)地理分區(qū)的青藏高原植被覆被變化及其對(duì)氣候變化的響應(yīng). 地理研究, 2010, 29(11): 2004-2016.

[11] Hua W, Fan G Z, Zhou D W,etal. Preliminary analysis on the relationships between Tibetan Plateau NDVI change and its surface heat source and precipitation of China. Science in China Series D-Earth Sciences, 2008, 51(5): 677-685.

[12] The Editorial Committee of the Ecological Protection and Economic Cycle in Qaidam. The ecological protection and economic cycle in Qaidam. Xining: Qinghai People’s Press, 2013.

《柴達(dá)木生態(tài)保護(hù)與循環(huán)經(jīng)濟(jì)》編輯委員會(huì). 柴達(dá)木生態(tài)保護(hù)與循環(huán)經(jīng)濟(jì). 西寧: 青海人民出版社, 2013.

[13] Li H M, Li Z W, Wang Z Y,etal. Comparative analysis of the climate abrupt changing time in different ecological function regions of Qinghai Province. Journal of Glaciology and Geocryology, 2012, 34(6): 1388-1393.

李紅梅， 李作偉， 王振宇， 等. 青海高原不同生態(tài)功能區(qū)氣候突變時(shí)間的比較分析. 冰川凍土， 2012, 34(6): 1388-1393.

[14] Chen X G, Li L, Zhu X D,etal. Regional differences of climate change in Qinghai Province and its contributing factors. Advances in Climate Change Research, 2009, 5(5): 249-254.

陳曉光, 李林, 朱西德, 等. 青海省氣候變化的區(qū)域性差異及其成因研究. 氣候變化研究進(jìn)展, 2009, 5(5): 249-254.

[15] Li L, Chen X G, Wang Z Y,etal. Climate change and its regional differences over the Tibetan Plateau. Advances in Climate Change Research, 2010, 6(3): 181-186.

李林, 陳曉光, 王振宇, 等. 青藏高原區(qū)域氣候變化及其差異性研究. 氣候變化研究進(jìn)展, 2010, 6(3): 181-186.

[16] Shi X H, Zhao Y N, Dai S,etal. Research on climatic change of Qaidam Basin since 1961. Journal of Desert Research, 2005, 25(1): 123-128.

時(shí)興合, 趙燕寧, 戴升, 等. 柴達(dá)木盆地40多年來的氣候變化研究. 中國沙漠, 2005, 25(1): 123-128.

[17] Li Y F, Yang T B. The temperature ascends in the Qaidam Basin in the past 50 years and global warming. Journal of Shangrao Normal College, 2005, 25(3): 105-109.

李永飛, 楊太保. 近50年來柴達(dá)木盆地升溫與全球變暖. 上饒師范學(xué)院學(xué)報(bào), 2005, 25(3): 105-109.

[18] Chen B S, Pan A D, Yang M Z. Spatial and temporal characteristics and change trend of climatic elements of Qaidam Basin in recent 50 years. Journal of Arid Land Resources and Environment, 2010, 24(5): 117-123.

陳碧珊, 潘安定, 楊木壯. 近50年柴達(dá)木盆地氣候要素分布特征及變化趨勢分析. 干旱區(qū)資源與環(huán)境, 2010, 24(5): 117-123.

[19] Tang H Y, Li X F. Preliminary analysis on maximum and minimum temperature change in Qinghai plateau during recent 40 years. Plateau Meteorology, 1999, 18(2): 230-235.

唐紅玉, 李錫福. 青海高原近40年來最高和最低溫度變化趨勢的初步分析. 高原氣象, 1999, 18(2): 230-235.

[20] Wang J, Huang Q H, Bai C G,etal. Tendency of the quaternary climatic change in Qaidam Basin and its causal mechanism. Scientia Geographica Sinica, 2002, 22(1): 34-38.

王建, 黃巧華, 柏春廣, 等. 2.5Ma以來柴達(dá)木盆地的氣候干濕變化特征及其原因. 地理科學(xué), 2002, 22(1): 34-38.

[21] Lü S N, Li D L, Wen J,etal. Analysis on periodic variations and abrupt change of air temperature over Qinghai-Xizang Plateau under global warming. Plateau Meteorology, 2010, 29(6): 1378-1385.

呂少寧, 李棟梁, 文軍, 等. 全球變暖背景下青藏高原氣溫周期變化與突變分析. 高原氣象, 2010, 29(6): 1378-1385.

[22] Li L, Shen H Y, Li H M,etal. Regional differences of climate change in Qaidam Basin and its contributing factors. Journal of Natural Resources, 2015, 30(4): 641-650.

李林, 申紅艷, 李紅梅, 等. 柴達(dá)木盆地氣候變化的區(qū)域顯著性及其成因研究. 自然資源學(xué)報(bào), 2015, 30(4): 641-650.

[23] Dai S, Shen H Y, Li L,etal. Analysis on climate transition characteristic from warm-dry to warm-wet in Qaidam Basin. Plateau Meteorology, 2013, 32(1): 211-220.

戴升, 申紅艷, 李林, 等. 柴達(dá)木盆地氣候由暖干向暖濕轉(zhuǎn)型的變化特征分析. 高原氣象, 2013, 32(1): 211-220.

[24] Shi Y F, Shen Y P, Hu R J. Preliminary study on signal, impact and foreground of climatic shift from warm-dry to warm-humid in Northwest China. Journal of Glaciology and Geocryology, 2000, 24(3): 219-226.

施雅風(fēng), 沈永平, 胡汝驥. 西北氣候由暖干向暖濕轉(zhuǎn)型的信號(hào)、影響和前景的初步探討. 冰川凍土, 2000, 24(3): 219-226.

[25] Li L, Wang Z Y, Qin N S,etal. Characteristics and tendency forecasting of dry/wet climate evolution over Qaidamu Basin in recent 1100 years. Plateau Meteorology, 2005, 24(3): 326-330.

李林, 王振宇, 秦寧生, 等. 近1100 年來柴達(dá)木盆地干濕氣候演變特征及趨勢預(yù)測. 高原氣象, 2005, 24(3): 326-330.

[26] Li Y P, Yang T B. Study on abrupt and periodic changes of temperature and precipitation in the Qaidam Basin. Geography and Geo-Information Science, 2007, 23(3): 105-108.

李遠(yuǎn)平, 楊太保. 柴達(dá)木盆地氣溫、降水突變與周期特征分析. 地理與地理信息科學(xué), 2007, 23(3): 105-108.

[27] Ren J Z, Hu Z Z, Mu X D. Bioclimatic indicators of the first class classification of grassland types in China. Journal of Gansu Agricultural University, 1965, (1): 48-64.

任繼周, 胡自治, 牟新待. 我國草原類型第一級(jí)分類的生物氣候指標(biāo). 甘肅農(nóng)業(yè)大學(xué)學(xué)報(bào), 1965, (1): 48-64.

[28] Ren J Z, Hu Z Z, Mu X D,etal. Comprehensive sequential classification of grassland and its significance in grassland genesis. Grassland in China, 1980, 3(1): 12-24.

任繼周, 胡自治, 牟新待, 等. 草原綜合順序分類法及其草原發(fā)生學(xué)意義. 中國草原, 1980, 3(1): 12-24.

[29] Hu Z Z, Gao C X. Improvement of the comprehensive and sequential classification system of grasslands indices if grassland classes and index chart. Acta Prataculturae Sinica, 1995, 4(3): 1-7.

胡自治, 高彩霞. 草原綜合順序分類法的新改進(jìn)-類的劃分指標(biāo)及其分類的檢索圖. 草業(yè)學(xué)報(bào), 1995, 4(3): 1-7.

[30] Du T Y. Classifying native grasslands of Qinghai Province with the integrated classification system. Pratacultural Science, 1992, 5: 28-32.

杜鐵英. 用綜合順序分類法對(duì)青海草地分類的探討. 草業(yè)科學(xué), 1992, 5: 28-32.

[31] Li H M, Ma Y S. Application on classification of Qinghai grassland by advanced comprehensive and sequential classification. Acta Prataculturae Sinica, 2009, 18(2): 76-82.

李紅梅, 馬玉壽. 改進(jìn)的綜合順序分類法在青海草地分類中的應(yīng)用. 草業(yè)學(xué)報(bào), 2009, 18(2): 76-82.

[32] Zhang Y L, Wei S C. Research on classification of Inner Mongolia grassland by complex-ordinal system. Grassland of China, 1990, 18(2): 14-20.

張永亮, 魏紹誠. 用綜合順序分類法對(duì)內(nèi)蒙古草原分類的研究. 中國草地, 1990, 5: 14-20.

[33] Du Q, Sun S Z. Vegetation and its utilization in Chaidamu region. Beijing: Science Press, 1990: 7-11.

杜慶, 孫世洲. 柴達(dá)木地區(qū)植被及其利用. 北京: 科學(xué)出版社, 1990: 7-11.