周曉曦，阿布都熱合曼·哈力克

(新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830046)

0 引 言

氣候變化是全球環(huán)境變化的重要體現(xiàn)，也是改變河川徑流量的主要因素[1]。氣候變化對(duì)徑流的影響呈現(xiàn)出多樣性和復(fù)雜性的特征，主要是通過氣溫和降水的變化影響河流徑流的補(bǔ)給、蒸發(fā)和流域的匯水等[2]。在劇烈變化的全球大環(huán)境下，氣溫上升顯著，水文過程的下墊面條件隨之改變，導(dǎo)致一系列水資源問題，尤其是對(duì)河流徑流的影響日益趨增，所以近年來關(guān)于氣候變化對(duì)徑流的影響已經(jīng)成為國內(nèi)外研究學(xué)者的熱點(diǎn)問題。在國外，Jiang C[3]等根據(jù)Budyko型方程來分析氣候變化和人類活動(dòng)對(duì)徑流的貢獻(xiàn)，結(jié)果表明氣候變化是影響徑流的主要驅(qū)動(dòng)因素；Zhang L[4]等分析了引發(fā)了黃河及其支流中游徑流量的顯著下降的原因主要是大規(guī)模的植被恢復(fù)和氣候變化。而在國內(nèi)，楊霞[5]利用SWAT模型研究了烏倫古河流域氣候變化對(duì)徑流影響，得出降水不變，流域氣溫升高1℃，徑流減少1.60%，徑流與氣溫呈負(fù)相關(guān)關(guān)系；胡浩云等[6]利用模型對(duì)9種氣候情景和土地利用組合進(jìn)行徑流模擬，揭示了氣候變化對(duì)流域的徑流量的影響程度不同；劉貴花等[7]對(duì)贛江徑流進(jìn)行氣候變化和人類活動(dòng)定量分析研究，表明降水是影響贛江流域徑流年際變化的主要因素。

新疆位于中國西北干旱區(qū)，而且也是對(duì)氣候變化響應(yīng)比較敏感的地區(qū)之一。根據(jù)相關(guān)文獻(xiàn)表明，西北干旱地區(qū)的水資源系統(tǒng)主要以冰雪融水為基礎(chǔ)，在全球氣候變化背景下，表現(xiàn)出的脆弱性非常明顯[8]。氣候變化對(duì)西北干旱區(qū)水資源影響的研究成為社會(huì)各界關(guān)注的焦點(diǎn)。新疆南部地區(qū)是我國冰雪對(duì)河流徑流補(bǔ)給較豐富的地區(qū)，在徑流對(duì)氣候變化的響應(yīng)中，冰雪補(bǔ)給徑流對(duì)氣候變化最為敏感，因此是影響河流徑流變化的最主要自然因素。受全球氣候變化的影響，新疆近年來氣溫有明顯的上升趨勢(shì)，河流徑流量呈現(xiàn)增加的勢(shì)頭?；谶@種情況，不少學(xué)者在新疆地區(qū)徑流變化趨勢(shì)研究方面開展了比較全面的研究工作，尤其針對(duì)氣溫、降水和冰川融化與徑流變化之間的關(guān)系作了系統(tǒng)分析[8-10]。而對(duì)于地處沙漠南緣的和田地區(qū)氣候變化對(duì)于和田河的相關(guān)研究較少，所以本文以和田河為研究區(qū)域，基于DPS、excel及標(biāo)準(zhǔn)值法采用數(shù)理方法來分析和田河的徑流變化特征，認(rèn)識(shí)其演化的周期性和趨勢(shì)性，不僅對(duì)認(rèn)識(shí)資源環(huán)境系統(tǒng)具有重要意義，而且對(duì)區(qū)域可持續(xù)性發(fā)展提供一定的科學(xué)依據(jù)。

1 研究區(qū)概況

和田河是塔里木河的第三大支流、是塔里木河的主要源流之一。和田河流域位于塔里木盆地南緣、地理位置介于東經(jīng) 77°25′～81°43′、北緯34°52′～40°28′、流域總面積48 870 km2[11]。流域包括的行政區(qū)為：和田行政區(qū)的和田縣、市、墨玉縣、洛浦縣、皮山縣和策勒縣的一小部分；阿克蘇行政區(qū)的阿瓦提縣一部分。其中東支玉龍喀什河發(fā)源與昆侖山北坡、河長504 km、多年平均徑流量22.5 億m3、徑流量較大、習(xí)慣上作為和田河主流。西支喀拉喀什河發(fā)源于昆侖山和喀喇昆侖山、河長808 km、多年平均徑流量21.7 億m3。兩河在闊什拉克(北緯38°05′、東經(jīng)80°33′)附近匯合后稱為和田河、從兩支流匯合口至和田河口末端(北緯40°29′10″、東經(jīng)80°56′40″)河長319 km。

2 資料和方法

研究區(qū)1953-2014年氣候資料降水(0.1 mm)、氣溫、濕度(1%)均在中國氣象數(shù)據(jù)共享網(wǎng)下載，和田河徑流量來自野外考察資料(同古孜洛克站和烏魯瓦提站)相加求年均值。

本文主要采取以下幾種方法：

(1) 線性趨勢(shì)方法。本文以時(shí)間序列為自變量，氣候條件為因變量，建立一元回歸方程。設(shè)y為某一氣象變量，t為時(shí)間(本文為年或月)，建立y與t之間的一元線性回歸方程，

y=a+bx

(1)

把b*10稱為變化傾向率，單位為℃/10 a或mm/10 a等，b值是反應(yīng)上升或者下降的變化趨勢(shì)，b>0表示在計(jì)算時(shí)段內(nèi)呈上升趨勢(shì)，b<0表示呈下降的趨勢(shì)，b值的絕對(duì)值的大小可以度量其演變趨勢(shì)上升或者下降的程度[12]。

(2)Mann-Kendall檢驗(yàn)法。Mann-Kendall是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法，在此利用統(tǒng)計(jì)軟件 DPS 對(duì)和田水文站1953-2014年的觀測(cè)徑流數(shù)據(jù)進(jìn)行M-K檢驗(yàn)，檢驗(yàn)的置信水平設(shè)為0.01。假定時(shí)間序列變量為X1，X2,…,Xn，n為時(shí)間序列的長度，構(gòu)造秩序列Sk:

(2)

其中：

式中：Xi，Xj分別為i，j年相應(yīng)的測(cè)量值，且i﹥j；在時(shí)間序列隨即獨(dú)立的假定下,定義統(tǒng)計(jì)量UFk：

(3)

式中UF1=0，E(Sk)，Var(Sk)是累計(jì)數(shù)Sk的均值和方差，可由下式求出：

E(Sk)=K(K-1)/4

Var(Sk)=k(k-1)(2k+5)/72

(4)

式中：k=2,3,…,n。

UFi為標(biāo)準(zhǔn)正態(tài)分布，是按照時(shí)間序列計(jì)算出的統(tǒng)計(jì)量序列,給定顯著性水平a，若|UFi|﹥Ua，則表明趨勢(shì)存在明顯的趨勢(shì)變化。按時(shí)間序列逆序Xn，Xn-1，…,X1，重復(fù)上述步驟，且使UBK=UFK(K=n，n-1，1)，UB1=0。

若UF或UB大于0，則表明序列呈上升趨勢(shì)，反之呈下降趨勢(shì)，當(dāng)他們超過臨界值時(shí)，表明上升或下降的趨勢(shì)顯著。且交點(diǎn)在臨界區(qū)域內(nèi)，那么交點(diǎn)對(duì)應(yīng)的時(shí)刻便是突變開始的時(shí)間[13]。

(3)功率譜分析。譜分析是非常重要的時(shí)間序列和空間序列分析方法。在時(shí)間序列分析中，這種方法叫作功率譜分析；在空間序列分析中，該方法叫作波譜分析。本文是采用功率譜分析來計(jì)算氣候因子的周期性，是根據(jù)陳彥光《基于Excel的地理數(shù)據(jù)分析》方法：

②快速傅里葉變換(FFT)：計(jì)算離散傅里葉變換的快速方法，有按時(shí)間抽取的FFT算法和按頻率抽取的FFT算法。前者是將時(shí)域信號(hào)序列按偶奇分排，后者是將頻域信號(hào)序列按偶奇分排。它們都借助于的兩個(gè)特點(diǎn)：一是周期性；二是對(duì)稱性。

③計(jì)算譜密度P(f)：

P(f)=1/T|S(f)|2=1/T(a2+b2)

(5)

式中：a為復(fù)數(shù)的實(shí)部(real part)；b為虛部(imaginary part)；T為樣本路徑長度，在本文實(shí)則補(bǔ)充0后的數(shù)據(jù)序列長度，對(duì)于本例，T=64。但是該方法比較繁瑣，所以本文使用的是Excel工程函數(shù)“imAbs”即，返回復(fù)數(shù)的絕對(duì)值(模)。

④繪制頻譜圖，首先要計(jì)算頻率，因?yàn)槎呦嗖畛?shù)倍(2π)，則f=時(shí)序/T，采用的頻率變化范圍 0～1，則繪制的頻譜圖是對(duì)稱的，實(shí)際上，另一半是多余的。因此，我們可以以對(duì)稱點(diǎn)f=0.5 為界，截取前面一半的數(shù)據(jù)，如果存在周期，則在尖峰突出的最大點(diǎn)可以找到，最后根據(jù)T1=1/f1(6)，算出周期[14]。

(4)變差系數(shù)Cv。Cv值反映年徑流量總體離散程度，Cv值大，年徑流的年際變化劇烈，這對(duì)水利資源的利用不利，而且易發(fā)生洪澇災(zāi)害；Cv值小，則年徑流量的年際變化小，有利于徑流資源用。

(7)

式中：n是觀測(cè)年數(shù)；ki是第i年的年徑流變化率。

為了消除年徑流量中Cv值的影響，本文運(yùn)用公式(K-1)/Cv(8)對(duì)徑流量進(jìn)行歸一化處理，其中為呈現(xiàn)流量的階段性變化情況，計(jì)算了其差積曲線，縱坐標(biāo)表示為∑(K-1)/Cv，K為模比系數(shù)。模比系數(shù)差積曲線主要反映徑流或降水序列是否經(jīng)歷豐水和枯水過程。差積曲線的上升段反映豐水期，下降段反映枯水期[15]。

3 和田河徑流量變化特征

由圖1可見，和田河年徑流量變化是呈波動(dòng)上升的趨勢(shì)，傾向率為0.3123，其中，1953-1977年徑流變化波動(dòng)幅度相對(duì)較小，而1978-2014年徑流變化幅度是相對(duì)較大，總體來看和田河徑流量是不斷地波動(dòng)上升的，這也和鄒朝望等的相關(guān)研究相符合[16]。從圖2和田河年徑流量模比系數(shù)差積曲線圖可見，1953-1975和1988-2001年是枯水期、1976-1987和2002-2014年是豐水期。雖然和田河有豐水期和枯水期，但是他的變化在零刻度線以下所以年際變化不明顯。

圖1 和田河徑流量年際變化曲線圖 Fig. 1 Hotan River runoff interannual variation graph

圖2 和田河年徑流量模比系數(shù)差積曲線圖Fig. 2 Hotan River annual runoff coefficient differential mode than the plot graph

4 和田河流域氣候變化特征

4.1 流域降水變化特征

如圖3可知，在0.05的顯著性水平下，曲線UF基本是在0以上，說明和田河流域歷年降水量是呈現(xiàn)增加的趨勢(shì)，但是波動(dòng)也比較大。UB和UF均未超過臨界直線，二者在臨界范圍之內(nèi)有交點(diǎn)，交點(diǎn)年份是2001年，說明2001是突變年份，但是和田站1953-2014年降水增加趨勢(shì)不顯著。從圖4可以看出，和田河流域頻率突出點(diǎn)較多，參差不齊，但是只有最高點(diǎn)是f=0.19，通過了在0.05顯著水平下的檢驗(yàn)，所以根據(jù)公式(6)得出T=5，說明和田河降水變化的周期是5 a。

總體來說，通過突變檢驗(yàn)可以得知：和田河流域降水呈增加的趨勢(shì)但是，波動(dòng)性比較大，增加不明顯，顯著性水平不高。通過功率譜分析，流域降水變化周期是5 a但突出點(diǎn)較多，周期變化幅度不大。

圖3 降水突變分析圖Fig. 3 Precipitation mutation analysis chart

圖4 和田河流域降水歷年功譜率圖Fig.4 Hotan River Basin precipitation rate over the years the power spectrum of FIG

4.2 流域氣溫變化特征

如圖5可知，在0.05的顯著性水平下，曲線UF從1978年以后都在0值以上而且超過了臨界線，說明和田河流域歷年溫度是呈現(xiàn)增加的趨勢(shì)及變化比較顯著。據(jù)曲線UF和曲線UF相交于1997年，說明1997年是突變年份，從1997年以后UF值逐漸增大，且超過了0.05顯著性水平。從圖6可以看出，頻率最高點(diǎn)是f=0.08，而且只有這個(gè)點(diǎn)通過了顯著性檢驗(yàn)，所以根據(jù)公式(6)得出T=12.5，說明和田河降水變化的周期大約是13 a。

總體來說，通過突變檢驗(yàn)得出：和田河流域溫度呈增加的趨勢(shì)而且顯著性水平高，從突變年份到2014溫度在持續(xù)遞增，這也符合近來全球氣候變化一致。通過功譜率分析得出：流域溫度變化周期是13 a，周期變化比較大。

圖5 氣溫突變分析圖Fig.5 Temperature mutation analysis chart

圖6 和田河流域溫度歷年功譜率圖Fig.6 Hotan River Basin over the years the power spectrum of temperature rate diagram

5 氣候變化對(duì)和田河徑流量的影響

5.1 影響徑流量變化的主要?dú)夂蛞蜃臃治?/h3>

為了對(duì)所得到的不同數(shù)據(jù)資料進(jìn)行比較及進(jìn)行相關(guān)性檢驗(yàn)分析和建立逐步回歸方程：即對(duì)數(shù)據(jù)做標(biāo)準(zhǔn)化處理是第一步，因?yàn)樵跉庀蟮拿總€(gè)影響因素當(dāng)中，每個(gè)要素的單位不相同，標(biāo)準(zhǔn)差和平均值也就會(huì)不一樣，計(jì)算出的結(jié)果也會(huì)有差異，所以為了使它們比較能在同一水平下進(jìn)行，常使用標(biāo)準(zhǔn)化的數(shù)據(jù)處理方法，讓他們成為無單位的變量在同一水平上，即消除量綱，這種變量稱為標(biāo)準(zhǔn)化變量。常用的數(shù)據(jù)標(biāo)準(zhǔn)化方法有最大值法、最小值法、標(biāo)準(zhǔn)差法等。本文采用如下標(biāo)準(zhǔn)化公式：

運(yùn)用公式對(duì)序列x1、x2、x3,…,xn進(jìn)行變換：

(9)

將和田河流域1953-2014年的氣溫、降水、濕度、徑流量進(jìn)行變換得到標(biāo)準(zhǔn)化值Ti(氣溫)、Pi(降水)、Ei(濕度)、Ri(徑流量)，然后繪出如圖7～圖9所示。

圖7 和田河流域年徑流量與氣溫標(biāo)準(zhǔn)化數(shù)據(jù)序列變化趨勢(shì)圖Fig.7 Annual runoff and temperature trends standardized data sequence Hotan River Basin

圖9 和田河流域年徑流量與濕度標(biāo)準(zhǔn)化數(shù)據(jù)序列趨勢(shì)圖Fig.9 Annual runoff and humidity standardized data sequence trend the Hotan River Basin

由圖6可知，和田河流域的氣溫和年徑流量呈現(xiàn)出總體增加的趨勢(shì)，而且變化趨勢(shì)非常接近，從20世紀(jì)80年代以后豐水階段和溫度的增長是一致的，以前枯水階段和溫度的增長也是一致的，因?yàn)楹吞锖訌搅髁康脑鲩L或減少主要是依靠冰雪融水補(bǔ)給，這也和全球溫度的升高有相關(guān)聯(lián)系。從圖7可知，年徑流量和降水的變化趨勢(shì)也基本一致但是降水相比年徑流量變化幅度較大，由上圖可看出和田地區(qū)年降水量偏少，主要是該地區(qū)深居內(nèi)陸，屬于荒漠氣候，降水稀少。從圖8可知，雖然濕度相比徑流量變化幅度較大，但是20世紀(jì)80年代前后就是一個(gè)分界線，80年代以前兩者的變化趨勢(shì)一致，以后變化幅度就有明顯的差異，主要是因?yàn)樵摰貐^(qū)受氣候變化和人類活動(dòng)的影響較小，80年代以后隨著全球氣候的變化，溫度升高，年降水量減少，蒸發(fā)量增加，濕度相對(duì)的減小所致。

5.2 氣候變化對(duì)徑流的影響

本文利用和田水文站1953-2014年的年降水序列、年平均溫度序列以及年濕度序列，與和田河1953-2014年的年徑流量序列在DPS中進(jìn)行逐步回歸分析并進(jìn)行顯著性檢驗(yàn)，結(jié)果得出和田河年徑流量與年平均溫度的相關(guān)系數(shù)為0.437，和年均降水量的相關(guān)系數(shù)為0.088，和年濕度的相關(guān)系數(shù)為-0.315，溫度和降水量通過了顯著水平a= 0.05的檢驗(yàn)，但是濕度沒能通過顯著水平的檢驗(yàn)，又因?yàn)闇囟葘?duì)徑流量的相關(guān)系數(shù)大于降水量對(duì)徑流量的相關(guān)系數(shù)，所以只選年均溫度被選入回歸方程?；貧w方程為：

Y徑流量=5.692 8X溫度-31.394 7

(10)

運(yùn)用以上回歸方程來進(jìn)行溫度對(duì)徑流量影響的量化研究，分析如表1所示。

表1 和田河流域不同溫度對(duì)徑流量的影響Tab.1 Effects of different temperatures on theHotan River Basin runoff

注：X為年均溫度，ΔX代表和田河氣溫距多年均值增加或減少的百分率；Y表示由式算出的對(duì)應(yīng)溫度的和田河年徑流量；ΔY表示和田河年徑流量距多年均值增加或減少的百分率。

由上表1可以看出，當(dāng)和田河流域溫度偏多或偏少2%時(shí)，和田河年徑流量會(huì)偏多3.2%或偏少0.8%，溫度偏多或偏少4%時(shí)，和田河徑流量會(huì)偏多6.8%或偏少6.7%；當(dāng)溫度偏多或偏少6%時(shí)，徑流量會(huì)偏多10.3%或偏少10.7%的響應(yīng)。

6 結(jié)果與結(jié)論

西北干旱區(qū)社會(huì)經(jīng)濟(jì)發(fā)展、生態(tài)環(huán)境的安全，以及未來經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展主要是受制約與水資源的變化。氣候變化是一個(gè)復(fù)雜的、多樣的系統(tǒng)問題，包含要素多，關(guān)系錯(cuò)綜復(fù)雜。隨著全球氣候變暖而引起的水資源的變化，將會(huì)使干旱區(qū)在開發(fā)和利用資源的過程中生態(tài)保護(hù)和經(jīng)濟(jì)發(fā)展的矛盾日益尖銳。影響徑流量的因素也很多，除了主要受氣候變化的影響，另外流域下墊面對(duì)徑流量也有影響，尤其是隨著社會(huì)經(jīng)濟(jì)的日益發(fā)展人類活動(dòng)對(duì)氣候變化和徑流量的變化也起著越來越重要的作用，而徑流量對(duì)氣候變化的響應(yīng)也是很復(fù)雜的。另一方面，由于受資料條件的限制，人類活動(dòng)因素在研究過程中未加以考慮，這需要在后繼工作中進(jìn)一步加以研究 本文通過和田河流域1953-2014年62 a的氣候資料(溫度、降水、濕度)以及水文站徑流資料進(jìn)行分析研究，得出以下結(jié)論：

(1)對(duì)和田河流域的降水、溫度做Mann-kendall檢驗(yàn)得知，降水變化雖然在臨界線之內(nèi)有突變點(diǎn)，但是降水降水變化顯著性不明顯，增加趨勢(shì)不大，溫度變化的顯著性比較大，1997年是突變年份，從1997年以后UF值逐漸增大，且超過了0.05顯著性水平，而且增加趨勢(shì)大。

(2)對(duì)和田河流域的降水、溫度做功率譜分析得出，和田河流域頻率突出點(diǎn)較多，但是只有最高點(diǎn)是f=0.19，通過了在0.05顯著水平下的檢驗(yàn)，所以根據(jù)公式(6)得出T=5，說明和田河降水變化的周期是5年。溫度頻率最高點(diǎn)是f=0.08，而且只有這個(gè)點(diǎn)通過了顯著性檢驗(yàn)，所以根據(jù)公式1.6得出T=12.5，說明和田河降水變化的周期大約是13 a。

(3)對(duì)和田河流域62 a的氣候因子(濕度、氣溫、降水量)序列和徑流進(jìn)行相關(guān)性分析，濕度同徑流間的相關(guān)系數(shù)為-0.315，降水量同徑流間的相關(guān)系數(shù)為0.088，溫度同徑流見的相關(guān)系數(shù)為0.437，根據(jù)相關(guān)系數(shù)顯著性檢驗(yàn)，降水和溫度的相關(guān)系數(shù)均通過了a=0.05的置信檢驗(yàn)，降水、溫度和徑流之間是正相關(guān)，說明他們是促使徑流量增加的要素；濕度與徑流之間是負(fù)相關(guān)，濕度的大小往往徑流的大小相反。

(4)和田河年徑流變差系數(shù)是0.28，反映出和田河徑流的多年變化的年變差系數(shù)相對(duì)較大，徑流年際變化不穩(wěn)定，總趨勢(shì)是在波動(dòng)上升，傾向率是0.312 3。其中1953-1975和1988-2001年是枯水期、1976-1987和2002-2014年是豐水期。雖然和田河有豐水期和枯水期但是他的變化在零刻度線以下所以年際變化不明顯。和田河流域溫度偏多或偏少2%時(shí)，和田河年徑流量會(huì)偏多或偏少-3.8%～2.9%，溫度偏多或偏少4%時(shí)，和田河徑流量會(huì)偏多或偏少-7%～6.5%；當(dāng)溫度偏多或偏少6%時(shí)，徑流量會(huì)偏多或偏少-11%～10%的響應(yīng)。

□

[1] David L C，Davisa J，Marta M J，et al．Space versus place incomplex human-natural systems: Spatial and multi-level models of tropical land use and cover change (LUCC) in Guatemala[J]．E-cological Modelling，2012，229:64-75．

[2] 張利平，于松延，段堯彬，等.氣候變化和人類活動(dòng)對(duì)永定河流域徑流變化影響定量研究[J].氣候變化研究進(jìn)展，2013， 9(6):391-397.

[3] Jiang C，Xiong L H，Wang D B, et al．Separating the impacts of climate change and human activities on runoff using the Budyko-type equations with time-varying parameters[J]．Journal of Hydrology， 2015，522:326-338．

[4] Zhang L，Podlasly C，F(xiàn)eger K et al．Different land management measures and climate change impacts on the runoff——A simple empirical method derived in a mesoscale catchment on the Loess Plateau[J]．Journal of Arid Environments，2015，120:42-50．

[5] 楊 霞. 基于SWAT模型的烏倫古河流域氣候變化對(duì)徑流影響研究[J]. 新疆環(huán)境保護(hù),2015,(1):45-50.

[6] 胡浩云,王 鋒,李宗華,等. 不同時(shí)段氣候變化對(duì)徑流的影響研究——以漳河上游濁漳干流流域?yàn)槔齕J]. 中國農(nóng)村水利水電,2015,(8):102-104,109.

[7] 劉貴花,齊述華,朱婧瑄,等. 氣候變化和人類活動(dòng)對(duì)鄱陽湖流域贛江徑流影響的定量分析[J]. 湖泊科學(xué),2016,(3):682-690.

[8] 陳亞寧，李 稚，范煜婷，等.西北干旱區(qū)氣候變化對(duì)水文水資源影響研究進(jìn)展[J]. 地理學(xué)報(bào)，2014，69(9):1 295-1 304.

[9] 梁立功.新疆塔河流域氣溫、降水和徑流變化特征[J]. 人民黃河，2014，36(8):24-27.

[10] 努爾比耶·艾合麥提托合提. 氣候變化與人類活動(dòng)對(duì)開都河徑流量的影響研究[D]. 烏魯木齊：新疆大學(xué),2015.

[11] 李衛(wèi)紅,黎 楓,陳忠升,等. 和田河流域平原耗水驅(qū)動(dòng)力與適宜綠洲規(guī)模分析[J]. 冰川凍土,2011,(5):1 161-1 168.

[12] 周 麗,謝舒蕾,范 君,等. 內(nèi)江市近43年來日照時(shí)數(shù)變化趨勢(shì)預(yù)測(cè)研究[J]. 內(nèi)江師范學(xué)院學(xué)報(bào),2015,(12):37-42,65.

[13] 唐啟義. DPS數(shù)據(jù)處理系統(tǒng)-實(shí)驗(yàn)設(shè)計(jì)、統(tǒng)計(jì)分析及數(shù)據(jù)挖掘[M]. 北京：科學(xué)出版社，2010:992,995.

[14] 陳彥光. 基于Excel的地理數(shù)據(jù)分析[M]. 北京：科學(xué)出版社, 2010.

[15] 白淑英,顧海敏,史建橋,等. 近50 a長江流域暴雨日數(shù)時(shí)空變化分析[J]. 長江流域資源與環(huán)境,2015,(7):1 255-1 262.

[16] 鄒朝望，尹耀鋒，孫媛媛．氣候變化下南疆地區(qū)徑流變化趨勢(shì)分析[J]．人民珠江，2016，37(1): 60-63．