氣候變暖背景下青南牧區(qū)牧業(yè)生產(chǎn)關(guān)鍵期變化特征及預(yù)估研究

劉彩紅，李紅梅，張調(diào)風(fēng)

(1．青海省氣候中心，青海 西寧810001;2． 南京信息工程大學(xué) 大氣科學(xué)氣象災(zāi)害省部共建教育部重點(diǎn)實(shí)驗(yàn)室，江蘇 南京210044)

青南牧區(qū)是青海省的主要牧區(qū)，屬三江源自然保護(hù)區(qū)生態(tài)環(huán)境脆弱的地區(qū)，海拔均在3 500 m 以上，現(xiàn)有天然草地面積2 141．75 萬hm2，草地類型以高寒草甸為主體，約占草地面積的76%，草地既是發(fā)展草地畜牧業(yè)的基礎(chǔ)性資源，又是江河上游的主要生態(tài)屏障［1］。畜牧業(yè)是青南牧區(qū)經(jīng)濟(jì)發(fā)展的支柱產(chǎn)業(yè)，而且也是當(dāng)?shù)厝嗣袢罕娰囈陨婧桶l(fā)展的基礎(chǔ)。在全球變暖的大背景下，1961 ―2012 年青南牧區(qū)年平均氣溫和各季節(jié)氣溫顯著升高，降水量總體呈增多趨勢(shì)，但具有明顯的區(qū)域性差異［2］。由于氣候生態(tài)環(huán)境蠕變，使得畜牧業(yè)結(jié)構(gòu)格局發(fā)生了顯著變化［3］。目前已觀測(cè)到的草地植被退化、草地載畜能力降低［4］、牧草生育期延長等［2］都被認(rèn)為與氣候變化有關(guān)。因此，在已認(rèn)識(shí)到的氣候條件變化對(duì)青南牧業(yè)生產(chǎn)影響的基礎(chǔ)上，探討未來氣候變化亦或氣候持續(xù)變暖對(duì)牧業(yè)造成的可能影響，已成為學(xué)者、公眾和決策者共同關(guān)心的問題。

近年來，國內(nèi)學(xué)者從不同角度對(duì)青藏高原氣候變化及其對(duì)農(nóng)牧業(yè)的影響作了大量研究［5-9］。但目前針對(duì)青南牧區(qū)的研究工作較少，尤其是氣候變化對(duì)牧業(yè)影響的預(yù)評(píng)估更是鮮有報(bào)道，因此，本研究擬采用滑動(dòng)平均、雙線性插值、突變檢驗(yàn)、小波分析等多種統(tǒng)計(jì)方法，在分析青南牧區(qū)牧業(yè)生產(chǎn)關(guān)鍵期(牧草生長季、青草期、牲畜抓膘期及掉膘期)變化特征的基礎(chǔ)上，選擇政府間氣候變化專門委員會(huì)(Intergovernmental Panel on Climate Change，IPCC)第5 次評(píng)估報(bào)告公布的CMIP5(Coupled Model Intercomparison Project phase 5)計(jì)劃21 個(gè)全球模式數(shù)值模擬結(jié)果，通過建立統(tǒng)計(jì)模型研究未來典型濃度路徑(Representative Concentration Pathways，RCPs)情景下青南高寒牧區(qū)生產(chǎn)關(guān)鍵期的變化趨勢(shì)，以期為當(dāng)?shù)乜茖W(xué)安排牧業(yè)生產(chǎn)及適應(yīng)氣候變化工作提供決策參考。

1 資料與方法

1．1 資料

考慮到臺(tái)站遷徙及資料的穩(wěn)定性，本研究選取1961 ―2013 年青南牧區(qū)沱沱河、澤庫、河南、玉樹、雜多、囊謙、瑪沁、達(dá)日、久治9 個(gè)地面氣象站(圖1)逐日氣溫觀測(cè)資料。氣候標(biāo)準(zhǔn)期為1971 ― 2000年。未來情景數(shù)據(jù)采用國際耦合模式比較計(jì)劃第5階段(CMIP5)［10］中21 個(gè)全球大氣與海洋環(huán)流耦合模式數(shù)值模擬集合平均結(jié)果。溫室氣體排放情景RCPs 是目前國際采用的最新氣候情景(表1)，能夠完全反應(yīng)氣候公約中穩(wěn)定大氣溫室氣體濃度的目標(biāo)［11］。

1．2 研究方法

1．2．1 指標(biāo)確定

牧草生長季:以春季日平均氣溫≥0 ℃的初日為牧草生長季的始期，秋季日平均氣溫≥0 ℃的終日為牧草生長季的終止期，≥0 ℃的持續(xù)日數(shù)為牧草生長季。

牧草青草期:牧草青草期的開始期與日平均氣溫≥5 ℃的初日比較接近，青草期終日與日平均氣溫≥5 ℃的終日基本一致，因此把日平均氣溫≥5℃的持續(xù)日數(shù)作為牧草青草期。

牲畜抓膘期:春季，當(dāng)日平均氣溫≥5 ℃時(shí)，各類家畜擺脫“春乏期”，逐漸恢復(fù)膘情。秋季，當(dāng)日平均氣溫下降到－5 ～5 ℃時(shí)，牲畜的膘情大致維護(hù)平衡，也就是不再長膘［9］。因而，本研究中把日平均氣溫≤5 ℃初日到≥0 ℃的終日之間的天數(shù)定義為抓膘期。

圖1 青南牧區(qū)氣象站點(diǎn)的分布圖Fig．1 The distribution of meteorological stations in Qingnan Pastoral areas

表1 典型濃度路徑特征Table 1 Characteristics of the representative concentration pathways

牲畜掉膘期:當(dāng)日平均氣溫穩(wěn)定≤－5 ℃時(shí)，牧草干枯，營養(yǎng)成分下降，放牧牲畜以消耗體內(nèi)脂肪來抗御外界低溫而開始掉膘，故把日平均氣溫≤－5℃的初日稱為牲畜掉膘的開始日，≤－5 ℃的終日稱為牲畜掉膘的結(jié)束日，≤－5 ℃之間的天數(shù)為牲畜掉膘期。

1．2．2 統(tǒng)計(jì)方法

(1)滑動(dòng)平均法

界限溫度的初、終日采用5 日滑動(dòng)平均方法確定［12］。對(duì)樣本量為n 的序列x，其滑動(dòng)平均序列表示為:

經(jīng)過滑動(dòng)平均后，可以消除序列中由于偶然引述所引起的不規(guī)則變動(dòng)，序列中短于滑動(dòng)長度的周期大大削弱，從而呈現(xiàn)出現(xiàn)象在較長時(shí)間的基本發(fā)展趨勢(shì)。

(2)雙線性差值

本研究采用雙線性差值法［13］將未來氣候場預(yù)估格點(diǎn)數(shù)據(jù)插值到站點(diǎn)上，與插值方法名稱不同的是，這種插值方法是非線性的。

假設(shè)已知函數(shù)f 在Q11= (x1，y1)、Q12= (x1，y2)、Q21=(x2，y1)以及Q22=(x2，y2)4 個(gè)點(diǎn)的值，得到未知函數(shù)f 在點(diǎn)P(x，y)的值:

首先在x 方向進(jìn)行線性插值，得到

然后在y 方向進(jìn)行線性插值，得到

則:

(3)突變分析

Mann-Kendall 方法是一種較為客觀的突變檢測(cè)方法，特點(diǎn)是能夠較為準(zhǔn)確地確定符合統(tǒng)計(jì)檢驗(yàn)的突變時(shí)間，避免人為因素的影響［14］。

Mann-kendall 方法首先對(duì)樣本時(shí)間序列Ti(i=1，2，3，…，n)定義統(tǒng)計(jì)量:

在時(shí)間序列Sk隨機(jī)獨(dú)立的假定下，定義統(tǒng)計(jì)量:

式中，E(Sk)= n(n －1)/4，Var(Sk)= n(n － 1)(2n+5)/72．

對(duì)統(tǒng)計(jì)量Uk按順序記為UF，同時(shí)按照逆順序i=n，n－1，n－2，…，1 計(jì)算一次并取負(fù)值，記為UB。取顯著性水平α=0．05。

(4)小波分析

小波分析在時(shí)域和頻域上同時(shí)具有良好的局部性質(zhì)，是一種窗口大小固定不變但其形狀可改變，時(shí)間窗和頻率都可以改變的時(shí)頻局部化分析方法，被譽(yù)為分析信號(hào)的顯微鏡［15］，在氣候分析中被廣泛應(yīng)用，此處選取的小波變換為Morlet 小波，其形式為:

2 結(jié)果與分析

2．1 牧業(yè)生產(chǎn)關(guān)鍵期變化特征

1961 ―2013 年青南牧區(qū)9 個(gè)站牧業(yè)生產(chǎn)關(guān)鍵期多年平均值統(tǒng)計(jì)顯示，牧草生長季平均開始期為4 月14 日，結(jié)束期10 月19 日，各地生長季持續(xù)日數(shù)為171．3 ～203．8 d。近53 年青南牧區(qū)牧草生長季呈現(xiàn)出開始期提早、結(jié)束期推遲趨勢(shì)(圖2a)，生長季明顯延長，速率為1．1 ～5．6 d·10a－1，其中囊謙延長的最多(圖3)。牧草青草期的始期也呈提早趨勢(shì)，其中在久治最為明顯，速率達(dá)到7 d·10a－1，結(jié)束期除囊謙略提早外，其余地區(qū)趨于推遲，但變化幅度均不如開始期明顯，幅度在0．6 ～3．1 d·10a－1之間。近53 年來，青南牧區(qū)青草期持續(xù)日數(shù)在64．7 ～133．2 d(圖2b)，平均以4．8 d·10a－1的速度延長，其中久治、澤庫為高值區(qū)。

由于氣溫顯著上升，1961 －2013 年青南牧區(qū)牲畜抓膘始期提前、結(jié)束期推遲，致使持續(xù)日數(shù)顯著延長(P ＜0． 05)，幅度介于2．2 ～9． 2 d·10a－1(圖2c)，其中以久治表現(xiàn)的最為突出。從牲畜掉膘期的變化可見，掉膘期呈現(xiàn)出開始期推遲、結(jié)束期提早、持續(xù)日數(shù)縮短的趨勢(shì)(圖2d)，開始期平均推遲了0．7 ～4．3 d·10a－1，結(jié)束期提早了0．4 ～3．2 d·10a－1，持續(xù)日數(shù)縮短了2．2 ～7．4 d·10a－1，玉樹地區(qū)掉膘期縮短幅度最大。

從6 次多項(xiàng)式擬合的趨勢(shì)變化整體來看(圖2)，青南牧區(qū)牧草生長季與青草期、牲畜抓膘期以及掉膘期均具有明顯的年代際波動(dòng)，20 世紀(jì)60 年代初至70 年代末前3 個(gè)指標(biāo)經(jīng)過短暫的低值后開始回升，80 年代初到90 年代末表現(xiàn)為持續(xù)時(shí)間的偏短期，自90 年代末期后呈明顯的高值階段，而掉膘期呈相反的變化趨勢(shì)。

綜上所述，1961 ―2013 年青南牧區(qū)牧草生長季與青草期，以及牲畜抓膘期均表現(xiàn)為開始期提早、結(jié)束期推遲、持續(xù)日數(shù)延長的變化特征，而牲畜掉膘期呈現(xiàn)出開始日推遲、結(jié)束日提早、持續(xù)日數(shù)縮短的變化特點(diǎn)。尤其是近15 年這種變化趨勢(shì)更為明顯，有利于青南高寒牧區(qū)的牧業(yè)生產(chǎn)。

圖2 1961 －2013 年青南牧區(qū)牧草生長季、青草期、牲畜抓膘期及掉膘期的變化趨勢(shì)Fig．2 The tendency of the indicators of animal husbandry critical period in Qingnan Pastoral areas during 1961 －2013

圖3 1961 －2013 年青南牧區(qū)代表站牧業(yè)生產(chǎn)關(guān)鍵期變化特征Fig．3 Linear trend of the pasture growing duration，green grass duration，fatten duration and fat-loss duration of livestock in Qingnan Pastoral areas during 1961 －2013

2．2 突變分析

目前對(duì)氣候突變的檢測(cè)比較客觀、準(zhǔn)確的方法是Mann-Kendall 方法(簡稱M-K 法)，但也只能對(duì)均值突變的檢測(cè)有把握，本研究利用M-K 法對(duì)青南牧業(yè)關(guān)鍵期53 年來是否存在氣候突變進(jìn)行分析。

1961 －2013 年青南牧區(qū)關(guān)鍵期的突變檢驗(yàn)結(jié)果(圖4)顯示，1961 ―2013 年青南牧區(qū)牧草生長季、青草期、牲畜抓膘期和掉膘期的持續(xù)日數(shù)都發(fā)生了顯著的氣候突變(P ＜0．05)，并且突變均發(fā)生在1997/1998 年，牧草生長季、青草期及牲畜抓膘期(1997/1998 年)由相對(duì)偏短躍變?yōu)槠L，而牲畜抓膘期則由相對(duì)偏長躍變?yōu)槠?。這與李林等［16］對(duì)青海地區(qū)氣溫在1998 年由暖向顯著變暖突變分析結(jié)果是一致的，同時(shí)與青藏高原氣溫在1998 年有一次突變發(fā)生的檢測(cè)結(jié)果相吻合［17］?？梢?，青南牧區(qū)由于氣溫升高而導(dǎo)致的牧業(yè)關(guān)鍵期變化是對(duì)全球變暖的一種響應(yīng)。

2．3 周期分析

采用Morlet 小波為基小波的分析中，依據(jù)小波系數(shù)的實(shí)部圖結(jié)合模數(shù)分析，不僅能反映各個(gè)周期成分在局部時(shí)段的特征，而且模數(shù)代表不同參數(shù)的小波對(duì)總能量的貢獻(xiàn)，能清楚的反映出試驗(yàn)序列中各個(gè)周期的成分的強(qiáng)度隨時(shí)間的變化［18］。根據(jù)牧草生長季、青草期、牲畜抓膘期和掉膘期時(shí)間序列進(jìn)行Morlet 小波變換，作出青南牧區(qū)牧業(yè)關(guān)鍵期Morlet 小波實(shí)部系數(shù)時(shí)－頻分析圖(圖5)，圖5 反映出在不同階段的同一周期振蕩以及同一階段的不同周期振蕩所表現(xiàn)出來的強(qiáng)弱程度是不一樣的。近53年牧草生長季沒有表現(xiàn)出明顯的振蕩周期，對(duì)于10年以下相對(duì)較小的時(shí)間尺度，牧草青草期、牲畜抓膘期8 ～12 年的周期突出(圖5)，是一能量高值區(qū)。在20 世紀(jì)60 年代中期到90 年代，牲畜抓膘期和掉膘期4 ～6 年周期也十分明顯(圖5)，其他周期信號(hào)強(qiáng)度都較弱。而對(duì)于10 年及以上相對(duì)較大的時(shí)間尺度，牧草青草期、掉膘期16 ～18 年左右的振蕩周期較為明顯(圖5)，具有較強(qiáng)的信號(hào)，這與突變分析及年代際變化的研究結(jié)果存在很好的對(duì)應(yīng)關(guān)系。

3 未來青南牧區(qū)牧業(yè)生產(chǎn)關(guān)鍵期趨勢(shì)預(yù)估

3．1 未來牧區(qū)氣候變化趨勢(shì)

CMIP5 試驗(yàn)全球模式對(duì)青海高原氣候模擬能力較為穩(wěn)定［19］，給出21 個(gè)模式集合對(duì)未來20 年青南牧區(qū)氣溫、降水變化的預(yù)估結(jié)果。根據(jù)RCP2．6、RCP4．5 和RCP8．5 排放情景下青南牧區(qū)年平均氣溫距平、降水距平百分率時(shí)間變化曲線(圖6)可以看出，2015 －2035 年青南牧區(qū)在3 種排放情景下年平均氣溫均呈升高趨勢(shì)(圖4)，由于大氣中溫室氣體的逐步增加，氣溫增幅與排放強(qiáng)度成正比，與氣候基準(zhǔn)年(1971 －2000年)相比，低、中、高3種排放情景下氣溫分別升高1．78、1．80、1．92 ℃。對(duì)于降水的未來變化情景，未來20 年青海牧區(qū)年降水量呈微弱的增多趨勢(shì)(圖4)，與氣候基準(zhǔn)年相比平均偏多2．7% ～3．9%。可見，未來較為暖濕的氣候條件將可以在一定程度上緩解因熱量、水分條件不足對(duì)牧草和牲畜的脅迫作用。

圖4 青南牧區(qū)牧業(yè)關(guān)鍵期時(shí)間序列的Mann-Kendall 突變檢驗(yàn)Fig．4 Mann-Kendall mutation testing of animal husbandry critical period in Qingnan Pastoral areas during 1961 －2013

圖5 青南牧業(yè)生產(chǎn)關(guān)鍵期Morlet 小波實(shí)部系數(shù)時(shí)－頻分析圖Fig．5 Time-frequency distribution of Morlet wavelet transform for animal husbandry critical periodover Qingnan Pastoral areas

3．2 未來生產(chǎn)關(guān)鍵期可能變化

根據(jù)青南牧區(qū)牧草生長季、青草期、牲畜抓膘期和掉膘期持續(xù)日數(shù)與年平均氣溫的線性相關(guān)模型，分別計(jì)算了3 種溫室氣體排放情景下(RCP2． 6，RCP4．5，RCP8．5)未來20 年青南牧區(qū)牧業(yè)生產(chǎn)關(guān)鍵期持續(xù)日數(shù)的變化情況(表2，3)。

與1971 ―2000 年基準(zhǔn)氣候時(shí)段比較，RCP2．6情景下，未來20 年青海牧區(qū)牧草生長季、青草期、牲畜抓膘期平均分別延長11．5、18．2 和18．5 d，牲畜掉膘期縮短15．7 d;RCP4．5 情景下牧草生長季、青草期、抓膘期分別延長12．1、19．2 和19．7 d，掉膘期縮短16．6 d;RCP8．5 情景下，牧草生長季、青草期、牲畜抓膘期分別延長13．2、20．9 和21．5 d，掉膘期縮短18．1 d(表3)?？梢?，在未來氣候變暖的背景下，青南牧區(qū)牧草生長季、青草期、牲畜抓膘期延長，而掉膘期縮短的趨勢(shì)較為明顯，適宜的氣候條件對(duì)當(dāng)?shù)氐哪翗I(yè)生產(chǎn)較為有利。但值得注意的是，近年來青南也出現(xiàn)了草畜不平衡，草地退化等問題，因此，抓住有利時(shí)機(jī)，合理調(diào)整畜種結(jié)構(gòu)，有效控制牲畜數(shù)量，對(duì)青南發(fā)展牧業(yè)經(jīng)濟(jì)十分重要。

4 討論與結(jié)論

圖6 2015 －2035 年青南牧區(qū)年平均氣溫距平及年降水量距平百分率變化趨勢(shì)Fig．6 The annual average temperature and annual precipitation anomaly time curve in Qingnan Pastoral areas during 2015 －2035

表2 青南牧區(qū)牧業(yè)生產(chǎn)關(guān)鍵期的持續(xù)日數(shù)與年平均氣溫的線性相關(guān)模型Table 2 The linear correlation model between duration of animal husbandry critical periodand annual mean temperature in Qingnan Pastoral areas

此前許多學(xué)者對(duì)氣象條件與牧草生長和牧事活動(dòng)關(guān)系進(jìn)行了研究。黃德昌等［20］提出了川西北高原以草食為主的畜牧業(yè)生產(chǎn)應(yīng)遵循的氣候規(guī)律。楊志華［21］認(rèn)為水、熱條件是影響綿羊體重增長量多寡的主要因素。李英年［22］通過分析與氣象條件相適應(yīng)的各類牧事活動(dòng)的最佳配套時(shí)間指出，日平均氣溫穩(wěn)定通過0 ℃期間，家畜能吃上草籽和適口性很好的優(yōu)質(zhì)牧草，家畜膘情持續(xù)穩(wěn)定，直到日均氣溫從－5 ℃開始時(shí)，家畜膘情才略有下降。張葆成［23］研究表明，藏東牦牛、藏綿羊、藏山羊和黃牛生長發(fā)育適應(yīng)溫度范圍為3 ～22 ℃，在這種溫度條件下，牧畜的生理機(jī)能較好，利于育肥抓膘。黃愛纖等［24］研究顯示，牧草產(chǎn)量隨3 ―5 月氣溫升高有增加的趨勢(shì)，3 ―5 月月平均氣溫每增加1 ℃，干草產(chǎn)量增加585．1 kg·hm－2，但降水影響不明顯。黃德青等［25］研究指出祁連山北坡草地生物量與降水量及相對(duì)濕度具有顯著的正相關(guān)關(guān)系。張核真等［9］分析得到，近40 年藏西北地區(qū)由于積溫和降水量顯著增多，牧草生長季、青草期提前，結(jié)束期推遲，持續(xù)日數(shù)增多，大部分地區(qū)牲畜抓膘期日數(shù)顯著增多，而掉膘期日數(shù)顯著減少，這一結(jié)論與本研究結(jié)果相一致?？梢姡跉夂蜃兣谋尘跋?，高原大部分地區(qū)氣溫升高、降水增加，水、熱條件的改善有利于牧草及牲畜的生長［26-28］。

表3 未來20 年青南牧區(qū)牧業(yè)生產(chǎn)關(guān)鍵期變化Table 3 Variation of the duration of animal husbandry critical period in Qingnan Pastoral areas during 2015 －2035 d

本研究通過氣候統(tǒng)計(jì)方法，分析了青南牧業(yè)生產(chǎn)關(guān)鍵期的時(shí)空變化特征，并結(jié)合模式模擬數(shù)據(jù)預(yù)估了未來20 年在RCPs 情景下牧業(yè)生產(chǎn)關(guān)鍵期變化趨勢(shì)，得出:

1)1961 －2013 年，青南牧區(qū)牧草生長季與青草期，以及牲畜抓膘期均表現(xiàn)為開始期提早、結(jié)束期推遲、持續(xù)日數(shù)延長的變化特征，延長速率分別為2．5，4．8，4．6 d·10a－1，而掉膘期呈現(xiàn)出開始日推遲、結(jié)束日提早、持續(xù)日數(shù)縮短的變化特點(diǎn)，縮短幅度為4．0 d·10a－1。

2)1961 －2013 年，青南牧區(qū)牧草生長季、青草期、牲畜抓膘期和掉膘期的持續(xù)日數(shù)都發(fā)生了氣候突變，突變時(shí)間均發(fā)生在1997/1998 年，牧草生長季、青草期及牲畜抓膘期從突變開始(1997/1998年)由相對(duì)偏短躍變?yōu)槠L，而牲畜掉膘期則由相對(duì)偏長躍變?yōu)槠?。?3 年牧草生長季沒有表現(xiàn)出明顯的振蕩周期;牧草青草期8 ～12 年、16 ～18年的周期突出，是一能量高值區(qū);牲畜抓膘期4 ～6年、8 ～12 年周期信號(hào)強(qiáng)度較強(qiáng);掉膘期具有4 ～6年、16 ～18 年的振蕩周期。

3)2015 ―2035 年，預(yù)估RCP2．6 情景下，青海牧區(qū)平均牧草生長季、青草期、牲畜抓膘期分別延長11．5、18．2 和18．5 d，掉膘期縮短15．7 d;RCP4．5情景下，牧草生長季、青草期、牲畜抓膘期分別延長12．1、19．2 和19．7 d，掉膘期縮短16．6 d;RCP8．5情景下，牧草生長季、青草期、牲畜抓膘期分別延長13．2、20．9和21．5 d，掉膘期縮短18．1 d。

本研究總體上對(duì)青南牧區(qū)氣候變化對(duì)牧業(yè)關(guān)鍵期的影響進(jìn)行了分析，為當(dāng)?shù)啬翗I(yè)發(fā)展提供了一定的科學(xué)參考，但研究中沒有對(duì)牧草產(chǎn)量及牧業(yè)關(guān)鍵期變化的原因進(jìn)行詳盡的分析，有待于進(jìn)一步探討。此外在未來預(yù)估中，氣候系統(tǒng)本身的復(fù)雜性和模式本身在模擬能力及表述氣候系統(tǒng)內(nèi)部各種物理過程的不完善，相應(yīng)的預(yù)估結(jié)果存在著一定程度的不確定性。

［1］ 王江山．青海省天氣氣候［M］．北京:氣象出版社，2004．

［2］ 青海省氣候變化監(jiān)測(cè)評(píng)估中心．青海省氣候變化評(píng)估報(bào)告［M］．北京:氣象出版社，2012．

［3］ 趙小娟．青南牧區(qū)生態(tài)畜牧業(yè)發(fā)展現(xiàn)狀分析及對(duì)策［J］．青海畜牧獸醫(yī)雜志，2008，38(5):34-37．

［4］ 韓國軍，王玉蘭，房世波．近50 年青藏高原氣候變化及其對(duì)農(nóng)牧業(yè)的影響［J］．資源科學(xué)，2011，33(10):1969-1975．

［5］ 卓嘎，楊秀海，羅文紅．西藏那曲地區(qū)氣候變化與牧業(yè)生產(chǎn)的關(guān)系［J］．資源科學(xué)，2009，31(3):485-492．

［6］ 宋春橋，游松財(cái)，劉高煥，柯靈紅，鐘新科．那曲地區(qū)草地植被時(shí)空格局與變化及其人文因素影響研究［J］． 草業(yè)學(xué)報(bào)，2012，21(3):1-10．

［7］ 楊秀海，卓嘎，羅布．藏北高原氣候變化與植被生長狀況［J］．草業(yè)科學(xué)，2011，28(4):626-630．

［8］ 格桑，卓瑪．淺析西藏高原氣候變化及其對(duì)農(nóng)牧業(yè)生產(chǎn)的影響［J］．西藏科技，2007(4):53-55．

［9］ 張核真，路紅亞，洪健昌，馬鵬飛． 藏西北地區(qū)氣候變化及其對(duì)草地畜牧業(yè)的影響［J］． 干旱區(qū)研究，2013，30(2):308-314．

［10］ Xu C H，Xu Y．The projection of temperature and precipitation over China under RCP scenarios using a CMIP5［J］．Atmospheric and Oceanic Science Letters，2012，5(6):527-533．

［11］ 馮婧．多全球模式對(duì)中國區(qū)域氣候的模擬評(píng)估和預(yù)估［D］．南京:信息工程大學(xué)碩士論文，2012:3-11．

［12］ 姜會(huì)飛．農(nóng)業(yè)氣象觀測(cè)與數(shù)據(jù)分析［M］．北京:科學(xué)出版社，2009．

［13］ 任永建，洪國平，肖鶯．長江流域上游氣候變化的模擬評(píng)估及其未來50 年情景預(yù)估［J］． 長江流域資源與環(huán)境，2013，22(7):894-900．

［14］ 符淙斌，王強(qiáng)．氣候突變的定義和檢測(cè)方法［J］．大氣科學(xué)，1992，16(4):482-493．

［15］ Torrence C，Compo G P．A practical guide to wavelet analysis［J］．Bulletin of the American Meteorological Society，1998，79:61-78．

［16］ 李林，陳曉光，王振宇，徐維新，唐紅玉．青藏高原氣候變化的最新事實(shí)及其區(qū)域差異性研究［J］．氣候變化研究進(jìn)展，2010，6(3):117-126．

［17］ 丁一匯，張莉．青藏高原與中國其他地區(qū)氣候突變時(shí)間的比較［J］．大氣科學(xué)，2008，3(4):794-805．

［18］ 吳洪寶，吳蕾．氣候變率診斷和預(yù)測(cè)方法［M］．北京:氣象出版社，2005:230-232．

［19］ 趙天保，陳亮，馬柱國．CMIP5 多模式對(duì)全球典型干旱半干旱區(qū)氣候變化的模擬與預(yù)估［J］．科學(xué)通報(bào)，2014，59:1148-1163．

［20］ 黃德昌，岳安榮，石現(xiàn)文．川西北高原主要的牧事活動(dòng)與氣象條件［J］．中國農(nóng)業(yè)氣象，1995，16(1):37-40．

［21］ 楊志華．放牧綿羊膘情氣象條件的初步探討［J］．內(nèi)蒙古氣象，1980(20):38-42．

［22］ 李英年．高寒地區(qū)氣象條件與季節(jié)草場及牧事活動(dòng)［J］．青海草業(yè)，1995，4(3):26-31．

［23］ 張葆成．藏東牧畜氣候適應(yīng)性探討［A］．第27 屆中國氣象學(xué)會(huì)年會(huì)現(xiàn)代農(nóng)業(yè)氣象防災(zāi)減災(zāi)與糧食安全分會(huì)場論文集［C］．北京:中國氣象學(xué)會(huì)，2010．

［24］ 黃愛纖，張新躍，唐川江，張緒校，魯巖． 川西北牧區(qū)水熱條件與牧草產(chǎn)量的相關(guān)性［J］． 草業(yè)科學(xué)，2015，32(5):754-759．

［25］ 黃德青，于蘭，張耀生，趙新全．祁連山北坡草地生物量及其與氣象因子的關(guān)系［J］．草業(yè)科學(xué)，2011，28(8):1495-1501．

［26］ 沈永平，王國亞．IPCC 第一工作組第五次評(píng)估報(bào)告對(duì)全球氣候變化認(rèn)知的最新科學(xué)要點(diǎn)［J］．冰川凍土，2013，35(5):1068-1076．

［27］ 梁云．影響天山北坡中山帶天然牧草產(chǎn)量的主要?dú)夂蛞蜃臃治觯跩］．草業(yè)科學(xué)，1999，16(3):1-3．

［28］ 董鎖成，周長進(jìn)，王海英．三江源地區(qū)主要生態(tài)環(huán)境問題與對(duì)策［J］．自然資源學(xué)報(bào)，2002，17(6):713-720．