朱中華，王雄師，柳金權

(1.甘肅省水利水電學校, 蘭州 730021；2.甘肅省科學院地質自然災害防治研究所, 蘭州 730000； 3.華亭縣水務局，甘肅 華亭 744100)

0 引 言

在全球氣候變暖的背景下，氣候極端事件和各種災害顯著增加，尤其是干旱災害[1]；同時溫度逐漸升高使得地表蒸發(fā)量增加, 干旱將會有進一步發(fā)展和加重趨勢[2-4]。土壤濕度作為地表主要物理參量之一，其對地表能量和水分再分配有著重要的作用，且能夠充分反映地表水分的收入和支出情況。20世紀80年代以來，采用田間實測方法對區(qū)域土壤水分狀況進行研究發(fā)展很快，并且取得了一定研究成果[5-14]。陳少勇等[15,16]分析了黃土高原土壤濕度的地域和時間分布特征以及土壤濕度的變化規(guī)律和土壤濕度對生態(tài)的影響，指出影響土壤濕度的主要影響因素是降水和氣溫。蔣沖等[17]運用線性趨勢和相關分析等方法對土壤相對濕度資料進行分析，指出黃土高原農田代表站點農田土壤濕度和主要氣象要素的年、月變化以及土壤濕度對氣候變化的響應。嚴麗等[18]采用線性趨勢法分析隴東主要氣象要素及土壤濕度的變化特征，指出20世紀90年代以來隴東各層土壤濕度總體上均呈現(xiàn)出下降趨勢，春季是土壤濕度減少最明顯的季節(jié)，夏季土壤濕度變化趨勢與春季具有類似規(guī)律，但變化率明顯低于春季。相對國內而言，國外學者針對土壤水分提出了一系列的計算模型。Mehrotra[19]針對印度中部的三個流域探討了影響氣候變化對水文循環(huán)的各個部分，即地表徑流，土壤水分和蒸散印度中部的三個流域，并假設場景降水和溫度變化作為輸入一個概念性的降雨徑流模型。Rushton等[20]對各種氣候條件是可能使用每日土壤水分平衡進行了概念模型和計算模型研究。通過國內外學者對土壤濕度的研究取得了一系列成果，并在實際的農業(yè)生產應用中取得了較好的效果；然而針對土壤濕度的研究，大部分只是分析了其與氣溫和降水等之間的關系，或者提出一些土壤模型的計算，并沒有考慮到土壤濕度與當?shù)氐膹搅髁筷P系，一個地區(qū)的徑流量對不同深度的土壤濕度在特定地區(qū)特定時間段內會有決定性的影響；同時干旱區(qū)小區(qū)域的土壤濕度的研究較少，但在實際的農業(yè)生產中，土壤濕度的研究對于一個小區(qū)域是非常有意義和必要性的。因此，文章以甘肅省靖遠縣為研究對象，基于靖遠縣農業(yè)站的降水、氣溫和徑流量資料，采用Matlab軟件、DPS數(shù)據(jù)處理系統(tǒng)、最小二乘法、SPSS軟件和Mann-Kendall方法等進行分析和計算及檢驗，最終得到靖遠縣不同作物各層土壤濕度對氣候變化響應關系的結果。這可以進一步發(fā)展節(jié)水灌溉農業(yè)，提高水資源利用率，進而緩解日益突出的水資源矛盾有十分重要的現(xiàn)實指導意義，同時，也可以為解決干旱區(qū)農業(yè)發(fā)展所面臨的水資源缺乏地區(qū)問題提供一定的借鑒和參考。

1 研究區(qū)概況

靖遠縣位于黃河上游，甘肅省中東部，處于東經104°13′～105°15′，北緯36°～37°15′。屬黃土高原溝壑區(qū)，地勢西高東低，由西北向東南傾斜，南北長125 km，東西寬106 km，總面積5 809.4 km2，海拔1 300～3 017 m之間。靖遠縣屬溫帶干旱半干旱氣候，每年均氣溫8.9 ℃左右，每年極端最高氣溫35.1 ℃，每年均降水量240 mm，每年蒸發(fā)量1 634 mm，每年平均日照時數(shù)2 696 h，無霜期165 d[21]。靖遠縣耕地面積7.67 萬hm2，其中有效灌溉面積3.66 萬hm2，小麥播種面積1.11 萬hm2左右，玉米播種面積1.02 萬hm2左右，水稻播種面積0.27 萬hm2左右等。靖遠縣總人口47.65萬人，其中農業(yè)人口43萬人。其中包括回、藏、滿、蒙古、東鄉(xiāng)等少數(shù)民族。

2 實驗數(shù)據(jù)

作物地段土壤濕度(重量含水量)觀測資料來自靖遠縣農業(yè)氣象觀測站，觀測時間序列為1981-2000年，種植作物類型分別為春小麥和春玉米，觀測土壤層次依次為0～10、10～20、20～30、30～40和40～50 cm。根據(jù)干旱區(qū)主要作物的需水規(guī)律，春小麥主要生長期為3-7月，春玉米主要生長期為4-9月，因此，年土壤濕度資料為不同作物生長發(fā)育期各旬資料的平均值。氣溫和降水資料序列為1956-2002年，來自靖遠縣氣象站，徑流量為還原后的天然出山徑流量，徑流序列為1956-2000年。通過分析1956-2002年氣溫、降水及徑流的變化規(guī)律和特征，其中，重點分析了1981-2000年20年的氣候變化和出山徑流變化特征，然后詳細地探討了不同層次土壤濕度的變化和分布特征。最終對土壤濕度與氣溫、降水和徑流進行了相關性分析和討論。

3 研究方法

3.1 土壤濕度

甘肅省氣象局根據(jù)省內各地的土壤質地，考慮實際氣候情況，結合多年的基礎研究及實際使用經驗，參照國家指標，確定了土壤濕度被用來代表甘肅省土壤旱澇檢驗指標；土壤濕度是土壤的干濕程度，即土壤的實際含水量，即土壤含水量與田間持水量的比值(公式1)。

(1)

3.2 最小二乘法

基于最小二乘法的應用背景和原理及算法；由于Matlab是一種向量化變成語言，因此在Matlab語言中對線性模型參數(shù)進行估計，進而可得到線性參數(shù)估計值[22,23]。

3.3 Mann-Kendall趨勢檢驗法

Mann-Kendall方法是一種檢測序列變化趨勢的方法[24-26]，其優(yōu)點在于所研究序列不需要遵從一定的分布，也不受少數(shù)異常值的干擾，更適合于類型變量和順序變量，計算比較簡單。文章采用該方法檢驗土壤濕度、氣溫和降水量的變化。

3.4 數(shù)據(jù)處理

通過采用Matlab軟件和DPS數(shù)據(jù)處理系統(tǒng)及Excel軟件，對研究區(qū)1981-2000年生長季內平均土壤濕度、平均氣溫和降水量的變化趨勢進行分析，同時運用Mann-Kendall方法進行趨勢檢驗，并應用最小二乘法和SPSS軟件對研究區(qū)土壤濕度與氣溫和降水量做線性擬合和相關性分析。

4 結果與分析

4.1 生長季作物土壤濕度變化

干旱半干旱區(qū)0～50 cm土層中土壤有效水分含量無論是在自然還是灌溉條件下，其變異率總是最大[27]。針對研究區(qū)土壤特點，土層土壤水分年變化過程在50 cm以下比較穩(wěn)定；因此，文章選取了研究區(qū)土層兩種作物生長季0～50 cm土壤濕度的變化趨勢及特征。

土壤濕度與氣象、水文要素以及地表覆蓋、土壤類型等有關。一般來說，同一土壤類型種植同樣的作物其土壤濕度高，農作物可以利用的土壤水分就多。同一作物類型地表各層土壤濕度隨時間的變化規(guī)律基本一致，而不同作物類型地表土壤濕度變化則差別較大。1981-2000年間，春小麥生長季(3-7月)相對濕度在不同土層深度年際變化趨勢較為相似，但也具有差異性；同時土壤濕度的年際變化具有3 a左右的周期。其中土層深度0～10和10～20 cm的土壤相對濕度呈現(xiàn)波動下降的趨勢，尤其是1988年到1990年下降幅度最大。土層深度20～30和30～40及40～50 cm的土壤相對濕度是波動增加的趨勢，特別是1990年到1993區(qū)間增加幅度最大(圖1)。春小麥各不同深度土壤濕度40～50 cm最大，除土壤表層(0～10 cm)受自然或人為因素干擾較大外，20～30 cm土壤濕度最小。

圖1 研究區(qū)春小麥生長季各層土壤土壤濕度變化Fig.1 The each soil moisture changes of the spring wheat growing season in the study area

春玉米生長季(4-9月)不同土層深度的相對濕度總體呈現(xiàn)波動上升趨勢，特別是自20世紀90年代以來，各層土壤濕度變幅增大，增加趨勢明顯加快；也存在變化周期為3年的年際變化；另外春玉米地塊各層土壤濕度的年際變化規(guī)律并不完全一致。其中0～10 cm增加幅度最大，尤其是1981-1985年、1989-1994年、1995-1998年；而1994-1995年下降幅度也大；由于人為因素的干擾，在1998年土壤相對濕度值達到最大值。土層深度10～20和20～30 cm的相對濕度在1989年到1994年增加的幅度最大。土層深度30～40 cm的相對濕度在1996-1999年之間的上升趨勢明顯。土層深度40～50 cm的土壤相對濕度1988年到1990年增加幅度最大(圖2)。

圖2 研究區(qū)春玉米生長季各層土壤濕度變化Fig.2 The each soil moisture changes of the spring corn growing season in the study area

4.2 氣候變化

氣溫和降水是兩個最重要的氣候因素，其直接影響糧食安全生產。氣溫狀況在作物的生長季節(jié)直接影響糧食的產量及質量。在干旱半干旱區(qū)，大部分區(qū)域是依靠天吃飯的雨養(yǎng)農業(yè)[28-30]。因此，對于研究區(qū)的氣溫和降水變化研究對靖遠縣糧食生產是有非常重要的意義。

4.2.1 氣 溫

1981-2000年，靖遠縣春小麥生長季(3-7月)氣溫變化趨勢與年平均氣溫變化趨勢基本一致，呈現(xiàn)周期式波動上升趨勢，年際變化周期約為3 a左右，1981年、1985年、1990年、1994年、1997年、2000年為周期變化峰值，而1983年、1986年、1989年、1992年、1996年、1999年為周期變化變化的谷值(圖3)。1983-1996年之間的4個周期波動振幅較小，變化相對平穩(wěn)。靖遠縣春小麥生長季(3-7月)平均氣溫在14.04～16.46 ℃之間。

圖3 研究區(qū)春小麥氣溫變化Fig.3 The temperature of spring wheat changes in the study area

研究區(qū)春玉米生長季(4-9月)氣溫變化趨勢與年平均氣溫變化趨勢基本一致，呈現(xiàn)增加趨勢，年際變化周期約為3 a左右，1981年、1985年、1989年、1994年、1998年、2000年為周期變化最大值，而1984年、1986年、1988年、1992年、1997年、1999年為周期變化變化的最小值(圖4)。1984-1992年之間的3個周期內變化幅度相對較小。靖遠縣春玉米生長季(4-9月)平均氣溫在17.27～19.28 ℃之間，較春小麥生長季溫度高出月3 ℃左右，原因在于研究區(qū)4-9月之間溫度普遍較高，且溫度在4月份后靖遠縣溫度逐步升高。

圖4 研究區(qū)春玉米氣溫變化Fig.4 The temperature of spring corn changes in the study area

4.2.2 降 水

研究區(qū)1981-2000年春小麥作物生長季(3-7月)降水與年平均降水趨勢變化很相似，總體呈現(xiàn)周期式波動增加，且年際變化較大；同時存在3～4 a的波動周期，在波動周期內，降水量呈現(xiàn)降低～升高～降低～升高～降低。1981年、1985年、1988年、1991年、1996年、1990年是周期變化的峰值，1982年、1987年、1989年、1995年、1997年、2000年為周期變化谷值(圖5)。靖遠縣春小麥生長季降水量在1982-1985年增加幅度最大，其中1982年達到降水量最小值(51.60 mm)；1997年到1999年降水量增加了115.70 mm，波動幅度也相對較大，其中1999年降水量為211.40 mm，為最大值。春玉米生長季降水量也存在著三個波動減少區(qū)間，其中1986-1987年降水量減少了60.7 mm；1988-1989年間降水量減少了65.8 mm；1999-2000年降水量下降幅度最大，減少了122.90 mm(圖5)。

圖5 研究區(qū)春小麥降水的變化Fig.5 The precipitation of spring wheat changes in the study area

靖遠縣1981-2000年春玉米作物生長季(4-9月)降水與年平均降水趨勢變化一致，降水量波動增加，且年際變化較大；波動周期約為2～3 a，在波動周期內，降水量呈現(xiàn)降低～升高～降低～升高～降低(圖6)；另外春玉米生長季降水量占年降水量比例均超過了80%，由于研究區(qū)降水主要集中在5-9月。1981年、1985年、1988年、1990年、1992年、1995年、1999年為變化周期內的最大值，1982年、1987年、1989年、1991年、1994年、1997年、2000年是波動周期內的最小值。1982-1985年間春玉米生長季降水量波動幅度最大，降水量增加了139.80 mm，其中1985年降水量達到了最大值(361.90 mm)，占年降水量的86.83%；然1982年降水量達到了最小值(122.10 mm)，占年降水量的89.19%。1997-1999年降水量增加了77.10 mm，增加幅度也較大。春玉米生長季降水量減少幅度有三個區(qū)間較大，其中從1985年到1987年間下降波動幅度最大，降水量減少了191.50 mm；1992年到1944年降水量減少了43.10 mm；1995-1997年降水量減少了93.2 mm；其余降水量波動周期減少值較小(圖6)。

圖6 研究區(qū)春玉米降水的變化Fig.6 The precipitation of spring corn changes in the study area

4.3 趨勢檢驗

利用Mann-Kendall方法對靖遠縣春小麥生長季(3-7月)和春玉米生長季(4-9月)土壤濕度以及氣溫和降水量的變化趨勢進行檢驗(圖7，圖8)。

圖7 春小麥生長季不同土層濕度和氣溫及降水的M-K檢驗Fig.7 The each soil humidity and air temperature and precipitation were inspect by using Mann-Kendall method to the spring wheat growing season

圖8 春玉米生長季不同土層濕度和氣溫及降水的M-K檢驗Fig.8 The each soil humidity and air temperature and precipitation were inspect by using Mann-Kendall method to the spring corn growing season

1981-2000年間，春小麥生長季(3-7月)表層土壤深度0～10 cm土壤濕度增加趨勢比土層深度10～20、20～30、30～40、40～50 cm和氣溫及降水的增加趨勢明顯，尤其在1986年后增加趨勢顯著。土層深度30～40 cm的土壤濕度在1982-1992年間呈下降趨勢；作物生長季氣溫在20世紀80年代初呈現(xiàn)下降趨勢，一直到20世紀90年代氣溫開始回升，逐漸增加。

1981-2000年間，春玉米生長季(4-9月)的土壤濕度和氣溫及降水呈現(xiàn)增加趨勢。其中土層深度0～10 cm在1995年呈現(xiàn)在增加趨勢，1998年增加趨勢顯著；10～20 cm土壤濕度1995年增加趨勢明顯；其余土層深度的土壤濕度增加趨勢不是很顯著。作物生長季氣溫在20世紀80年代初出現(xiàn)下降趨勢，到20世紀80年代末氣溫開始上升。

4.4 土壤濕度與氣溫和降水的關系

各層土壤濕度年際變化均有不同程度的增加(除春小麥地20～30 cm為負值外，其他各層斜率均為正值)，垂直方向上，隨著土層深度的增加，土壤濕度逐漸降低，但不同作物稍有不同。對于春小麥地塊來說，20～30 cm深度土壤濕度呈現(xiàn)微弱下降趨勢非常值得注意，因為這可能與春小麥屬于高耗水性作物，且其根系多集中分布在20～30 cm深度范圍內，多年持續(xù)種植同一種作物有關。對于春玉米來說，20～30、30～40和40～50 cm這三層的土壤濕度基本保持一致，有微弱的增加趨勢(表1)。不同深度土壤濕度多年平均值則與種植作物類型有一定的關系，春小麥地塊土壤濕度從上到下依次為低～高～低變化即呈“凸型”變化趨勢，春玉米地塊剛好與之相反，土壤濕度從上到下依次為高～低～高變化即呈“凹型”變化趨勢(表2)。

表1 土壤濕度各層年變化系數(shù)Tab.1 Variation coefficients of soil moisture in different soil layers

表2 不同作物土壤濕度各層年平均值Tab.2 The average soil moisture of each layer in different crops

土壤濕度的變化不僅與當?shù)氐臍夂蚝退臈l件有關，而且還與種植作物類型和土壤類型有一定的關系。隨著降水和徑流量的增加，春小麥地土壤濕度總體上略有增加趨勢，但不同層次的土壤濕度其增加幅度不同。降水和徑流量均與春小麥地各層土壤濕度呈正相關關系；氣溫與土層20～30 cm的土壤濕度負相關關系，氣溫與其他土層的土壤濕度呈正相關關系。另外，不同的影響因子對土壤濕度影響的層次也各不相同，其中氣溫和降水與0～10 cm土壤濕度相關性最大；徑流對20～30 cm土壤濕度影響最大(表3)。

氣溫和降水與春玉米地各層土壤濕度呈正相關關系，而徑流量則與土壤濕度呈負相關關系(表3)。其中，氣溫和降水與土層10～20 cm的土壤濕度相關性最大，徑流對30～40 cm土壤濕度影響最大。此外，氣溫與土層10～20和20～30 cm土壤濕度相關關系達到極顯著水平，而與0～10 cm土壤濕度相關關系達到顯著水平。

氣溫與土壤濕度的關系隨種植作物類型和土層深度不同而存在一定的差異；降水與土壤濕度的關系與種植作物類型和土層深度均無關，但它們的關系隨作物類型和土層深度敏感程度有一定差別；而徑流量與土壤濕度的關系隨作物類型不同有著完全相反的表現(xiàn)，這并不能說明徑流與土壤濕度關系因作物類型而定，更主要的是因距離河流遠近及地下水水位有著密切關系。此外，土壤濕度對氣候變化響應的敏感性和滯后性隨著季節(jié)和資料序列間隔等因素的不同也有一定的變化。

表3 氣溫、降水和徑流量與春小麥(春玉米)地塊土壤濕度相關性Tab.3 Correlation coefficients for spring wheat (spring corn) field between temperature,precipitation,runoff and soil moisture in different layers

注：**表示極顯著水平，*表示顯著水平。

5 結 語

通過對靖遠縣不同作物生長季土壤濕度和氣溫及降水進行研究，采用Matlab軟件和DPS數(shù)據(jù)處理系統(tǒng)及Microsoft Excel數(shù)據(jù)處理軟件進行分析，并結合最小二乘法進行線性參數(shù)估算，利用Mann-Kendall方法檢驗對不同作物各層土壤濕度和氣溫及降水的變化趨勢進行檢驗；同時運用SPSS軟件對不同土層的土壤濕度與氣溫、降水、徑流做相關分性分析，得出的主要結論如下。

同一作物類型地表各層土壤濕度隨時間的變化規(guī)律基本一致，而不同作物類型各層土壤濕度變化則差別較大。不同深度土壤濕度多年平均值則與種植作物類型有一定的關系，春小麥地塊土壤濕度從上到下依次為低～高～低變化即呈“凸型”變化趨勢，春玉米地塊剛好與之相反，土壤濕度從上到下依次為高～低～高變化即呈“凹型”變化趨勢。

氣溫與土壤濕度的關系隨種植作物類型和土層深度不同而存在一定的差異；降水與土壤濕度的關系與種植作物類型和土層深度均無關，但它們的關系隨作物類型和土層深度敏感程度有一定差別；而徑流量與土壤濕度的關系隨作物類型不同有著完全相反的表現(xiàn)，這并不能說明徑流與土壤濕度關系因作物類型而定，更主要的是因距離河流遠近及地下水水位有著密切關系。此外，土壤濕度對氣候變化響應的敏感性和滯后性隨著季節(jié)和資料序列間隔等因素的不同也有一定的變化。

[1] 劉 瑜,趙爾旭,黃 瑋,等. 2005年初夏云南嚴重干旱的診斷分析[J].熱帶氣象學報,2007,23(8):50-56.

[2] 安芷生,符淙斌. 全球變化科學的進展[J].地球科學進展,2001,16(5):671-680.

[3] 蘇志珠,盧 琦,吳 波,等. 氣候變化和人類活動對我國荒漠化的可能影響[J].中國沙漠,2006,26(3):329-335.

[4] 袁文平,周廣勝. 標準化降水指標與Z指數(shù)在我國應用的對比分析[J].植物生態(tài)學報,2004,28(4):523-529.

[5] 李旭輝,馮秀藻.西北地區(qū)旱地農田土壤濕度資料開發(fā)問題初探[J].干旱區(qū)農業(yè)研究,1989,(1)1:34-39.

[6] 馬柱國,魏和林,符淙斌. 土壤濕度與氣候變化關系的研究進展與展望[J]. 地球科學進展,1999,14(3):299-305.

[7] 周 濤,史培軍,范一大.中國北方土壤濕度變化趨勢及人類活動對其影響研究[J].北京師范大學學報,2002,38(1):131-137.

[8] 王 磊,文 軍,韋志剛,等. 中國西北區(qū)西部土壤濕度及其氣候響應[J]. 高原氣象,2008,27(6):1 257-1 266.

[9] 張文君,周天軍,宇如聰.中國土壤濕度的分布與變化I.多種資料間的比較[J].大氣科學,2009,32(3):581- 597.

[10] 郭維棟,馬柱國,姚永紅. 近50 年中國北方土壤濕度的區(qū)域演變特征[J].地理學報,2003,58(增刊):83-90.

[11] 林 潔,陳效民,張 勇. 氣候變化與土壤濕度關系的研究進展[J]. 土壤通報,2012,43(5):1 271-1 276.

[12] Mileham L, TaylorR, Thompson J. Impact of rainfall distribution on the parameterisation of a soil moisture balance model of groundwater recharge in equatorial Africa [J].Journal of Hydrology, 2008,359(1):46-58.

[13] Bennie J, Huntley B, Wiltshire A, et al. Slope, aspect and climate: Spatially explicit and implicit models of topographic microclimate in chalk grassland [J]. Ecological Modelling, 2008,216(1):47-59.

[14] Moroizumi T, Hamada H, Sukchan S, et al. Soil water content and water balance in rainfed fields in northeast Thailand [J].Agricultural Water Management, 2009,96(1):160-166.

[15] 陳少勇,董安祥. 中國黃土高原土壤濕度的氣候響應[J]. 中國沙漠,2008,28(1):66-72.

[16] 陳少勇,郭凱忠,董安祥. 黃土高原土壤濕度變化規(guī)律研究[J]. 高原氣象,2008,27(3):530-537.

[17] 蔣 沖,王飛,穆興民,等. 黃土高原農田土壤濕度演變及其與氣候變化的響應關系[J]. 干旱地區(qū)農業(yè)研究,2012,30(3):234-243.

[18] 嚴 麗,王 飛,蔣 沖,等. 隴東黃土高原農田土壤濕度演變對氣候變化的響應[J].水土保持通報2012,32(3):11-16.

[19] Mehrotra R. Sensitivity of runoff, soil moisture and reservoir design to climate change in central indian river basins[J],Climatic Change,1999,42(4):725-757.

[20] Rushton K R， M Eilers V H， Carter R C. Improved soil moisture balance methodology for recharge estimation[J]. Journal of Hydrology, 2006, 318(1-4):379-399.

[21] 楊 義,張寶軍,丁貞玉,等. 50 年來靖遠縣氣候變化趨勢分析[J].干旱區(qū)研究,2007,24(4):532-536.

[22] 王 可,毛志伋. 基于Matlab實現(xiàn)最小二乘曲線擬合[J]. 北京廣播學院學報(自然科學版),2005,12(2):52-56.

[23] 許弟春. 關于最小二乘法的參數(shù)估計問題探討[J]. 長春師范學院學報(自然科學版),2009,28(2):21-23.

[24] 符淙斌,王 強. 氣候突變的定義和檢測方法[J].大氣科學,1992,16(4):482-493.

[25] 馬柱國,魏和林,符淙斌. 中國東部區(qū)域土壤濕度的變化及其與氣候變率的關系[J].氣象學報,2000,58(3):278-287.

[26] 劉東飛,舒和平,馬金珠. 白龍江流域極端降水年際變化及對氣候變化的響應[J].水電能源科學,2015,33(3):1-5.

[27] 劉長民,趙凡衍. 旱地農田土壤水分含量變化特征研究[J].西北林學院學報,1995,10(增刊):148-152.

[28] 郭海英,楊興國,黃 斌,等. 隴東黃土高原主要農作物生長狀況評定指標體系的建立[J]. 干旱地區(qū)農業(yè)研究,2005,23(6):12-16.

[29] 吳愛敏,郭江勇,王勁松. 中國西北地區(qū)伏期干旱指數(shù)及干旱分析[J].干旱區(qū)研究,2007,24(2):227-233.

[30] 黃 斌,郭江勇,張洪芬,等. 隴東玉米拔節(jié)至抽雄期降水與產量及生物量的相關性分析[J]. 干旱地區(qū)農業(yè)研究,2007,25(1):172-175.