• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    1961-2010年白龍江流域氣溫和降水量變化特征研究

    2014-09-21 02:23:49張曉曉徐浩杰
    水土保持研究 2014年4期
    關(guān)鍵詞:白龍江年際平均氣溫

    張曉曉, 張 鈺, 徐浩杰

    (蘭州大學(xué) 資源環(huán)境學(xué)院, 蘭州 730000)

    1961-2010年白龍江流域氣溫和降水量變化特征研究

    張曉曉, 張 鈺, 徐浩杰

    (蘭州大學(xué) 資源環(huán)境學(xué)院, 蘭州 730000)

    以1961—2010年白龍江流域氣溫和降水量資料為基礎(chǔ),綜合運(yùn)用線性趨勢(shì)分析、多項(xiàng)式回歸、Mann-Kendall趨勢(shì)性檢測(cè)和突變分析、復(fù)Morlet小波分析和R/S分析等數(shù)理統(tǒng)計(jì)方法,在季節(jié)和年際尺度上研究了氣溫和降水量的趨勢(shì)性、波動(dòng)性、突變性、周期性和持續(xù)性。結(jié)果表明:在研究時(shí)段內(nèi)四季和年平均氣溫均呈顯著上升趨勢(shì),20世紀(jì)80年代中期之前升溫緩慢,90年代以來(lái)升溫速率明顯加快并發(fā)生暖突變。年降水量呈“多—少—多—少”的波動(dòng)變化,整體微弱減少。春季和冬季降水量分別呈“多—少—多”和“少—多—少”變化,整體微弱增加。夏季和秋季分別呈“多—少—多—少”和“多—少”變化,整體微弱減少。氣溫和降水量存在7,11,22 a左右的準(zhǔn)周期變化。未來(lái)一段時(shí)間內(nèi),四季和年平均氣溫將可能繼續(xù)升高,夏季和冬季降水量將可能增加,而春季、秋季和年降水量將可能減少。

    白龍江; 氣候變化; 趨勢(shì); 突變; 周期

    氣候是自然環(huán)境的重要組成部分,其變化對(duì)人類生存、經(jīng)濟(jì)社會(huì)發(fā)展等將產(chǎn)生極其深遠(yuǎn)的影響。近百年來(lái)全球氣候正經(jīng)歷以變暖為主要特征的顯著變化,自1860年以來(lái)全球地表溫度升高了0.74℃,變暖幅度自20世紀(jì)90年代以來(lái)明顯加速,未來(lái)100 a全球氣溫可能升高1.1~6.4℃[1]。伴隨氣溫的顯著變化,全球降水量的時(shí)空變化格局亦發(fā)生明顯改變,熱帶地區(qū)和高緯度地區(qū)降水量趨于增多,副熱帶地區(qū)降水量趨于減少[2]。在全球氣候變化背景下,我國(guó)氣候也發(fā)生了顯著而深刻的變化。近百年來(lái)我國(guó)地表溫度氣溫升高約0.5~0.8℃,略高于同期全球升溫幅度平均值,且近50 a來(lái)增暖尤為明顯[3]。于此同時(shí),我國(guó)年降水量呈微弱減少趨勢(shì),其中西部地區(qū)(特別是西北地區(qū))降水量增加趨勢(shì)明顯,而西南地區(qū)降水量呈一定的減少趨勢(shì)[4]。

    在大氣環(huán)流控制下,受地形、自然條件、人類活動(dòng)的影響,氣候變化存在時(shí)空異質(zhì)性,加強(qiáng)區(qū)域氣候變化研究,對(duì)于評(píng)估區(qū)域氣候變化現(xiàn)狀、原因及其影響,指導(dǎo)區(qū)域經(jīng)濟(jì)社會(huì)和生態(tài)環(huán)境的可持續(xù)發(fā)展具有重要意義[5]。白龍江流域地處青藏高原與川西北高原交錯(cuò)帶,是中國(guó)西部亞熱帶和暖溫帶的分界線。受地形和氣候變化影響,區(qū)域生態(tài)環(huán)境脆弱,自然災(zāi)害頻發(fā),氣候系統(tǒng)的穩(wěn)定性嚴(yán)重影響著該區(qū)經(jīng)濟(jì)社會(huì)發(fā)展和生態(tài)環(huán)境保護(hù)[6]。白龍江流域氣溫變化20世紀(jì)80年代之前趨勢(shì)不明顯,80年代至90年代之間波動(dòng)明顯,呈下降趨勢(shì),90年代中后期氣溫年際變化呈顯著上升趨勢(shì),增溫率以冬季最大[7-8]。降水量年際變化90年代中期之前上下波動(dòng)趨勢(shì)不顯著,90年代中期以后呈下降趨勢(shì),并以夏季降水量遞減為主[7-8]。白龍江流域徑流受氣候變化和人類活動(dòng)影響,年際變化趨勢(shì)逐步減少[8]。目前關(guān)于白龍江流域氣候變化的研究較少,且多集中于氣象站點(diǎn)尺度的分析。由于研究?jī)?nèi)氣象站點(diǎn)分布稀少且空間分布不均勻,若干站點(diǎn)平均結(jié)果難以準(zhǔn)確反映研究區(qū)氣溫和降水量的整體變化情況。因此本文基于1961—2010年白龍江流域氣象站點(diǎn)資料,對(duì)氣溫和降水量進(jìn)行空間插值并計(jì)算流域氣溫和降水量的空間平均值,對(duì)其季節(jié)和年際變化特征進(jìn)行線性趨勢(shì)分析、平滑預(yù)測(cè)分析、周期性分析、突變分析和持續(xù)性分析,以揭示氣溫和降水量的變化特征,對(duì)理解氣候過(guò)渡帶氣候變化過(guò)程具有重要意義,可服務(wù)于該區(qū)經(jīng)濟(jì)社會(huì)發(fā)展和自然災(zāi)害預(yù)警等諸多方面。

    1 數(shù)據(jù)與方法

    1.1 研究區(qū)概況

    白龍江屬長(zhǎng)江二級(jí)支流,嘉陵江一級(jí)支流,河流全長(zhǎng)535 km,流域面積32 810 km2,發(fā)源于甘肅省碌曲縣與四川省若爾蓋縣交界的郎木寺鄉(xiāng),向東流經(jīng)碌曲、若爾蓋、迭部、舟曲、宕昌、武都、文縣、青川及廣元9個(gè)縣(市),于昭化縣舊城北部匯入嘉陵江,主要支流包括岷江、白水江、達(dá)拉河、大團(tuán)峪河、拱壩河和清江河。白龍江流域位于秦巴山地西部,青藏高原東緣,橫跨甘肅和四川兩省,受構(gòu)造運(yùn)動(dòng)和流水侵蝕作用影響,該區(qū)多高山峽谷地貌,地貌類型包括山地、丘陵和盆地,地理位置介于102°31′—105°45′E,32°10′—34°23′N,海拔介于450~4 900 m,相對(duì)高差1 000~2 000 m,地勢(shì)由西北向東南遞減。白龍江流域氣候類型復(fù)雜多樣,具有明顯的季風(fēng)氣候特征。武都水文站(104°55′E,33°24′N)以上流域?yàn)榘埥嫌?,屬溫帶濕?rùn)氣候區(qū)。武都至碧口水文站(105°15′E,32°45′N)流域?yàn)榘埥杏?,屬暖溫帶濕?rùn)氣候區(qū)。碧口水文站以下流域?yàn)榘埥掠?,屬于北亞熱帶濕?rùn)氣候區(qū)。流域年平均溫度約8~10℃,年降水量約500~900 mm,由東南向西北遞減,4—9月降水量約占全年降水量的85%,多為強(qiáng)降水形式。

    1.2 數(shù)據(jù)來(lái)源與研究方法

    選取1961—2010年白龍江流域及其周邊區(qū)域20個(gè)氣象站(連續(xù)觀測(cè)且未遷站)的月平均氣溫和月降水量數(shù)據(jù),具體包括瑪曲、碌曲、臨潭、卓尼、岷縣、郎木寺、迭部、宕昌、若爾蓋、禮縣、西和、成縣、武都、南坪、文縣、康縣、平武、青川、廣元、劍閣。由于研究?jī)?nèi)氣象站點(diǎn)分布稀少且空間分布不均勻,站點(diǎn)平均結(jié)果難以準(zhǔn)確反映研究區(qū)氣溫和降水量的整體變化情況。本文根據(jù)各氣象站點(diǎn)的經(jīng)度、緯度、海拔信息,逐年逐月對(duì)各氣象站點(diǎn)氣溫(降水量)和經(jīng)度、緯度、海拔進(jìn)行多元線性回歸,并對(duì)殘差值(真實(shí)值與擬合值之差)進(jìn)行樣條插值,得到時(shí)空分辨率為月與1 km的氣象柵格數(shù)據(jù)[9]。運(yùn)用ArcGIS空間統(tǒng)計(jì)模塊,逐年逐月統(tǒng)計(jì)白龍江流域氣溫和降水量的空間平均值(研究區(qū)內(nèi)若干柵格數(shù)據(jù)平均值)。采用交叉檢驗(yàn)法,以誤差平均值和誤差標(biāo)準(zhǔn)差為指標(biāo),對(duì)插值前后的數(shù)據(jù)進(jìn)行t檢驗(yàn),結(jié)果表明在99%的置信度下絕大部分?jǐn)?shù)據(jù)的顯著性概率均大于0.01,插值結(jié)果可靠。

    按照3—5月為春季、6—8月為夏季、9—11月為秋季、12月—翌年2月為冬季進(jìn)行季節(jié)劃分,研究白龍江流域氣溫和降水量的季節(jié)變化特征。運(yùn)用線性回歸和6次多項(xiàng)式擬合法,研究氣溫和降水量的多年變化趨勢(shì)及年際波動(dòng)規(guī)律。運(yùn)用Mann-Kendall(MK)法[10]和滑動(dòng)t檢驗(yàn)法[11]研究氣溫和降水量變化趨勢(shì)的顯著性及其突變性。運(yùn)用復(fù)Morlet小波分析[12]研究氣溫和降水量的周期變化。為了反映各準(zhǔn)周期的顯著性,本文采用0.05信度紅/白噪聲標(biāo)準(zhǔn)譜進(jìn)行檢驗(yàn)[13]。運(yùn)用R/S分析法[14]預(yù)測(cè)氣溫和降水量在未來(lái)一段時(shí)間內(nèi)可能的變化趨勢(shì)。

    1.2.1 Mann-Kendall法 當(dāng)時(shí)間序列x1,x2,…,xn相互獨(dú)立并具有相同連續(xù)分布時(shí),定義統(tǒng)計(jì)量:

    (1)

    式中:k=2,…,n;n——樣本數(shù);UFk——MK統(tǒng)計(jì)量,且UF1=0;秩序列Sk——第i時(shí)刻數(shù)值大于j時(shí)刻數(shù)值個(gè)數(shù)的累計(jì)數(shù);E(Sk)和var(Sk)——累計(jì)數(shù)Sk的均值和方差。

    利用MK統(tǒng)計(jì)值進(jìn)行趨勢(shì)顯著性檢驗(yàn)時(shí),通過(guò)兩尾檢驗(yàn)于正態(tài)分布表中查出一定顯著性水平下的臨界檢驗(yàn)值u(u0.05=±1.96,u0.01=±2.57 )。當(dāng)MK統(tǒng)計(jì)值大于臨界檢驗(yàn)值時(shí),認(rèn)為時(shí)間序列變化趨勢(shì)顯著,同時(shí)MK統(tǒng)計(jì)值為正時(shí)表示序列呈上升趨勢(shì),否則呈下降趨勢(shì)。MK統(tǒng)計(jì)值的絕對(duì)值越大,表明序列變化趨勢(shì)越顯著。將時(shí)間序列逆序,再重復(fù)MK計(jì)算過(guò)程,得到逆序列MK統(tǒng)計(jì)值UBk。在UFk曲線超過(guò)臨界值的前提下,若2組統(tǒng)計(jì)值僅有1個(gè)交叉點(diǎn)且位于臨界值之間,這點(diǎn)便是突變點(diǎn)的開(kāi)始。若交叉點(diǎn)位于臨界值之外,或者是存在多個(gè)明顯的交叉點(diǎn),則不確定是否為突變點(diǎn),需利用其它方法檢測(cè)突變點(diǎn)顯著性。

    1.2.2 Morlet小波分析法 對(duì)于時(shí)間序列f(t),Morlet小波分析的母函數(shù)為:

    ψ(t)=π-1/4eiω0te-t2/2

    (2)

    式中:ω0是角頻率,取ω0=6。

    離散小波變換公式:

    (3)

    式(3)中:“*”——復(fù)共軛;a——尺度因子(與周期和頻率有關(guān));b——平移因子(時(shí)間位置);Wf(a,b)——小波系數(shù);δt為資料序列時(shí)間間隔。

    小波方差為:

    Wp(a)=Wf(a,b)2

    (4)

    在一定的時(shí)間尺度下,小波方差表示時(shí)間序列在該尺度中周期波動(dòng)的強(qiáng)弱,小波方差隨尺度的變化過(guò)程能反映時(shí)間序列中所包含的各種時(shí)間尺度及其強(qiáng)弱隨尺度的變化特征,對(duì)應(yīng)峰值處的尺度即為該序列的準(zhǔn)周期。

    1.2.3R/S分析法 給定一時(shí)間序列Xi(i=1,2,…,n),對(duì)于任意正整數(shù)n,定義均值序列:

    (5)

    累積離差序列:

    (6)

    極差序列:

    Rξ=maxYξ-minYξ(ξ=1,2,…,n)

    (7)

    標(biāo)準(zhǔn)差序列:

    (8)

    (9)

    則表明所分析的時(shí)間序列存在Hurst現(xiàn)象,H稱為Hurst指數(shù)。若0.5

    2 結(jié)果與分析

    2.1 氣溫和降水量的年際變化特征

    白龍江流域多年平均氣溫為8.85℃。從年平均氣溫的線性變化趨勢(shì)來(lái)看(圖1),1961—2010年年平均氣溫呈極顯著上升趨勢(shì),年際變化率為+0.197℃/10a,且在未來(lái)一段時(shí)間內(nèi)增溫還將持續(xù)(表1)。多項(xiàng)式函數(shù)擬合的年平均氣溫曲線顯示,1961—1985年氣溫多年平均值為8.59℃,與整個(gè)研究時(shí)段相比處于低溫期,且呈微弱下降趨勢(shì),年際變化率為-0.046℃/10a,其年際波動(dòng)呈“降低(1961—1965)—升高(1966—1972)—降低(1973—1985年)”變化。1985年以后,氣溫急劇上升,年際變化率為+0.374℃/10a,1986—2010年的年平均氣溫較1961—1985年偏高0.52℃。從氣溫的年代際變化來(lái)看(圖1),20世紀(jì)60年代平均氣溫最低,70年代和80年代基本相當(dāng)且增溫緩慢,90年代以后氣溫快速升高,21世紀(jì)以來(lái)升溫速率進(jìn)一步加快,其平均氣溫較60年代偏高0.84℃。

    圖1 年平均氣溫的年際變化

    從氣溫的季節(jié)變化來(lái)看(圖2),春季平均氣溫在1961—1988年呈先升高(1961—1967年)再降低(1968—1988年)的波動(dòng)變化,1988年以后氣溫快速上升,2007年以后出現(xiàn)一定回落。夏季平均氣溫在1961—1986年呈“降低(1961—1965年)—升高(1966—1972年)—降低(1973—1985年)”的波動(dòng)變化,1986年以后氣溫快速升高,2001年以后趨于平穩(wěn)。秋季平均氣溫除在上世紀(jì)60年代中期以前出現(xiàn)短暫降溫以外,其余時(shí)間段均呈上升趨勢(shì)。冬季平均氣溫在整個(gè)研究時(shí)段內(nèi)均呈上升趨勢(shì)。從氣溫的線性變化趨勢(shì)及顯著性來(lái)看(圖2和表1),冬季升溫速率最快,其次是秋季、夏季和春季,其年際變化率分別為+0.335,+0.165,+0.158,+0.127℃/10a。冬季氣溫變化趨勢(shì)最為顯著,其次是夏季、秋季和春季。從四季氣溫的Hurst指數(shù)來(lái)看,未來(lái)一段時(shí)間內(nèi)四季氣溫仍將可能顯著上升(表1)。四季氣溫的年代際變化結(jié)果表明,四季平均氣溫自1990s以來(lái)快速增加,春季、夏季、秋季和冬季平均氣溫分別在1980s,1980s,1970s和1960s最低,春季、夏季、秋季和冬季平均氣溫在2000s最高(表2)。

    圖2 四季平均氣溫的年際變化表1 氣溫變化趨勢(shì)定量檢驗(yàn)與預(yù)測(cè)

    氣溫Mann-Kendall秩次相關(guān)系數(shù)M值趨勢(shì)Mα=0.05Mα=0.01顯著性Hurst指數(shù)未來(lái)趨勢(shì)春季2.1665升高±1.96±2.57顯著0.7271升高夏季2.9026升高±1.96±2.57極顯著0.6904升高秋季2.5346升高±1.96±2.57顯著0.6402升高冬季4.2577升高±1.96±2.57極顯著0.7383升高全年4.3916升高±1.96±2.57極顯著0.8252升高

    白龍江流域降水量的多年平均值為684.5 mm。從年降水量的線性變化趨勢(shì)來(lái)看(圖3),1961—2010年白龍江流域降水量呈微弱減少趨勢(shì),其年際變化率為-9.7 mm/10 a,且減少趨勢(shì)在未來(lái)一段時(shí)間內(nèi)還將持續(xù)(表3)。降水量在上世紀(jì)60年代偏多,70年代偏少,80年代有所回升,90年代和2000年以來(lái)再次偏少(表2)。多項(xiàng)式函數(shù)擬合的年降水量曲線顯示,白龍江流域年降水量經(jīng)歷了“多—少—多—少”的波動(dòng)變化。1961—1964年降水量整體偏多;1965—1973年降水量整體偏少,在1967年出現(xiàn)降水高值,在1965年出現(xiàn)降水低值;1974—1989年降水量整體偏多,在1981年出現(xiàn)降水量高值,在1986年、1987年出現(xiàn)降水低值;1990—2010降水量整體偏少,在1998年、2003年和2005年出現(xiàn)降水量高值,在1997年、2002年、2006年出現(xiàn)降水量低值。

    圖3 年降水量的年際變化

    從降水量的季節(jié)變化來(lái)看(圖4),春季、夏季、秋季、冬季降水量分別約占全年降水量的21.46%,52.06%,24.35%和2.13%。春季降水量在1961—1978年偏多,1979—1995年偏少,1996—2010年再次偏多,在1967年、1973年出現(xiàn)降水高值,在1962年、1979年出現(xiàn)降水低值。夏季降水量的年際波動(dòng)形式與年降水量十分一致,表明夏季降水量對(duì)年降水量影響顯著。秋季降水量在1961—1983年偏多,1984—2010年偏少,在1967,1975,1978年出現(xiàn)降水高值,在1972,1976,1987年出現(xiàn)降水低值。冬季降水量在1961—1982年偏少,1983—2005年偏多,2006—2010年再次偏少,在1990,1991,1994年出現(xiàn)降水高值,在1987,1999,2010年出現(xiàn)降水低值。從降水量的線性變化趨勢(shì)及顯著性來(lái)看,春季和冬季降水量呈微弱上升趨勢(shì),其年際變化率分別為+1.31,+0.8 mm/10 a,夏季和秋季降水量呈微弱下降趨勢(shì),其年際變化率為-5.1,-6.8 mm/10 a,春季降水量增加幅度最大,秋季降水量減少幅度最大。從四季降水量的Hurst指數(shù)來(lái)看,未來(lái)一段時(shí)間內(nèi)春季、秋季和年降水量將可能呈減少趨勢(shì),夏季和冬季降水量將可能呈增加趨勢(shì)(表2)。

    表3表示氣溫和降水量的年代際變化,從表3可以看出,四季降水量的年代際變化結(jié)果表明,春季降水量在1990s最低,2000s最高;夏季降水量在2000s最低,1980s最高;秋季降水量在1990s最低,1960s最高;冬季降水量在1970s最低,1990s最高。

    圖4 四季降水量的年際變化表2 降水量變化趨勢(shì)定量檢驗(yàn)與預(yù)測(cè)

    降水量Mann-Kendall秩次相關(guān)系數(shù)M值趨勢(shì)Mα=0.05Mα=0.01顯著性Hurst指數(shù)未來(lái)趨勢(shì)春季0.5939增加±1.96±2.57不顯著0.2069減少夏季-0.6943減少±1.96±2.57不顯著0.4587增加秋季-1.7315減少±1.96±2.57不顯著0.5960減少冬季1.4973增加±1.96±2.57不顯著0.5907增加全年-0.9787減少±1.96±2.57不顯著0.5173減少

    表3 氣溫和降水量的年代際變化

    2.2 氣溫和降水量的突變特征

    從氣溫和降水量Mann-Kendall突變檢測(cè)結(jié)果來(lái)看(圖5),四季和年平均氣溫的正序列曲線分別在2008年、2003年、2006年、1990年和2000年左右超過(guò)顯著性為0.05的置信度曲線,同時(shí)在置信度曲線區(qū)間內(nèi)僅存在一個(gè)交叉點(diǎn),表明春季、夏季、秋季、冬季、年平均氣溫分別在2002年、1993年、1994年、1990年、1996年發(fā)生暖突變。四季和年降水量在置信度曲線區(qū)間內(nèi)存在多個(gè)交叉點(diǎn),運(yùn)用5年滑動(dòng)t檢驗(yàn)對(duì)交叉點(diǎn)進(jìn)行顯著性為0.05的置信度檢驗(yàn),結(jié)果表明在研究時(shí)段內(nèi)除冬季降水量在1972年發(fā)生由少到多的顯著性突變外,其余時(shí)段降水量不存在顯著突變特征。

    圖5 氣溫和降水量的Mann-Kendall突變分析

    2.3氣溫和降水量的周期變化特征

    氣溫和降水量小波方差顯著性檢驗(yàn)結(jié)果表明氣溫和降水量不存在顯著變化周期,僅存在準(zhǔn)周期(圖6),其中春季平均氣溫準(zhǔn)周期為11 a,夏季平均氣溫準(zhǔn)周期為13 a,此外6~8 a的周期也比較明顯,秋季平均氣溫準(zhǔn)周期為8 a,冬季平均氣溫準(zhǔn)周期為11 a,此外6 a和22 a周期也較為明顯,年平均氣溫準(zhǔn)周期為12 a,其次還存在一個(gè)6 a左右的次級(jí)周期。春季降水量準(zhǔn)周期為7 a,夏季降水量準(zhǔn)周期為25 a,此外還存在6 a和15 a左右的次級(jí)周期,秋季降水量準(zhǔn)周期為13 a,次級(jí)周期為6 a,冬季和年降水量準(zhǔn)周期分別為18 a和25 a,同時(shí)也存在6~8 a的小周期。

    圖6 氣溫和降水量的Morlet小波方差

    3 結(jié)論與討論

    通過(guò)對(duì)白龍江流域1961—2010年氣溫和降水量的變化特征分析,氣溫呈顯著上升趨勢(shì),降水量總體呈波動(dòng)變化,但有微弱減少趨勢(shì),總體氣候趨于暖干,對(duì)全球變化響應(yīng)明顯。

    (1) 白龍江流域年平均氣溫上升趨勢(shì)顯著,年際變化率為0.197℃/10 a,其中1961—1985年氣溫呈微弱下降趨勢(shì),1985年以后氣溫急劇上升,并在1996年左右發(fā)生暖突變,突變后升溫速率進(jìn)一步加快。從年平均氣溫的年代際變化來(lái)看,20世紀(jì)60年代平均氣溫最低,70年代和80年代基本相當(dāng)且增溫緩慢,90年代以后氣溫快速升高,21世紀(jì)以來(lái)升溫速率進(jìn)一步加快。四季平均氣溫均呈顯著上升趨勢(shì),冬季升溫速率最快,其次是秋季、夏季和春季,其年際變化率分別為0.335,0.165,0.158,0.127℃/10 a。春季、夏季、秋季、冬季平均氣溫分別在2002年、1993年、1994年、1990年發(fā)生暖突變。四季平均氣溫分別在1980s、1980s、1970s和1960s最低,在2000s最高。

    (2) 白龍江流域年降水量線性變化趨勢(shì)不明顯,但有微弱減少趨勢(shì),年際變化率為-9.7 mm/10 a。年降水量經(jīng)歷了“多—少—多—少”的波動(dòng)變化。1961—1964年降水量整體偏多;1965—1973年降水量整體偏少,在1967年出現(xiàn)降水高值,在1965年出現(xiàn)降水低值;1974—1989年降水量整體偏多,在1981年出現(xiàn)降水量高值,在1986年、1987年出現(xiàn)降水低值;1990—2010降水量整體偏少,在1998年、2003年和2005年出現(xiàn)降水量高值,在1997年、2002年、2006年出現(xiàn)降水量低值。春季和冬季降水量呈微弱增加趨勢(shì),夏季和秋季降水量呈微弱減少趨勢(shì)。春季降水量在1990s最低,2000s最高;夏季降水量在2000s最低,1980s最高;秋季降水量在1990s最低,1960s最高;冬季降水量在1970s最低,1990s最高,且在1972年發(fā)生由少到多的顯著突變。

    (3) 白龍江流域氣溫主要存在8 a、11~13 a和22 a左右的準(zhǔn)周期,降水量主要存在7 a、13~18 a和25 a左右的準(zhǔn)周期。

    (4) 未來(lái)一段時(shí)間內(nèi),白龍江流域四季和年平均氣溫仍將可能繼續(xù)升高,夏季和冬季降水量將可能微弱增加,而春季、秋季和年降水量將可能微弱減少。

    近50 a來(lái)白龍江流域氣溫顯著升高,降水量微弱減少,這在一定程度上加劇了白龍江流域的暖干化發(fā)展態(tài)勢(shì),將影響到區(qū)域的經(jīng)濟(jì)發(fā)展、生態(tài)與環(huán)境等諸多方面,首先是水資源,干旱化趨勢(shì)引起地表徑流減少,將進(jìn)一步影響區(qū)域的水資源配置與工農(nóng)業(yè)生產(chǎn)[15]。其次,氣候變暖使該區(qū)生態(tài)環(huán)境發(fā)生著變化[16],氣候變暖使蒸騰蒸發(fā)量增大,降水變率增大,極端降水事件(旱澇災(zāi)害)的頻率和強(qiáng)度增加,自然災(zāi)害發(fā)生的可能性進(jìn)一步增加[17]。

    [1] IPCC. Summary for Policymakers of the Synthesis Report of the IPCC Fourth Assessment Report[R]. Cambridge, UK:Cambridge University Press,2007:5-6.

    [2] 施能,黃先香,楊揚(yáng).1948—2000年全球陸地年降水量場(chǎng)趨勢(shì)變化的時(shí)、空特征[J].大氣科學(xué),2003,27(6):971-982.

    [3] 唐國(guó)利,任國(guó)玉.近百年中國(guó)地表氣溫變化趨勢(shì)的再分析[J].氣候與環(huán)境研究,2005,10(4):791-798.

    [4] 王英,曹明奎,陶波,等.全球氣候變化背景下中國(guó)降水量空間格局的變化特征[J].地理研究,2006,25(6):1031-1040.

    [5] 秦大河,丁一匯,蘇紀(jì)蘭,等.中國(guó)氣候與環(huán)境演變?cè)u(píng)估(Ⅰ):中國(guó)氣候與環(huán)境變化及未來(lái)趨勢(shì)[J].氣候變化研究進(jìn)展,2005,1(1):4-9.

    [6] 柳金峰,黃江成,歐國(guó)強(qiáng),等.白龍江中游泥石流沖積扇坡耕地開(kāi)發(fā)潛力分析[J].水土保持研究,2011,18(1):92-96.

    [7] 趙艷萍,寧娜,馬金珠.白龍江流域近40 a氣候變化及徑流的響應(yīng)[J].節(jié)水灌溉,2012(5):6-10.

    [8] 張秀云,姚玉璧,王潤(rùn)元.白龍江流域氣候變化及其對(duì)水資源的影響[J].資源科學(xué),2009,31(8):1315-1320.

    [9] 徐浩杰,楊太保.1981—2010年柴達(dá)木盆地氣候要素變化特征及湖泊和植被響應(yīng)[J].地理科學(xué)進(jìn)展,2013,32(6):868-879.

    [10] Mann H B. Nonparametric tests against trend[J]. Econometrica,1945,13(3):245-259.

    [11] 魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M].北京:氣象出版社,2007.

    [12] Goupillaud P, Grossmann A, Morlet J. Cycle-octave and related transforms in seismic signal analysis[J]. Geoexploration,1984,23(1):85-102.

    [13] 尹繼堯,朱元清,宋治平,等.Morlet小波顯著性檢驗(yàn)和精度分析在地磁場(chǎng)和地震活動(dòng)性周期分析中的應(yīng)用[J].地震學(xué)報(bào),2011,33(5):663-671.

    [14] Hurst H E. Long-term storage capacity of reservoirs[J]. Trans. Amer. Soc. Civil Eng.,1951,116:770-808.

    [15] 馬正耀,胡興林,藍(lán)永超,等.1965—2010年白龍江上游徑流變化特征研究[J].冰川凍土,2011,33(3):612-618.

    [16] 李小亞,張勃,靳自寶.基于MODIS-NDVI的甘肅河?xùn)|地區(qū)植被覆蓋度動(dòng)態(tài)監(jiān)測(cè)[J].水土保持研究,2013,20(1):112-115.

    [17] 蘇布達(dá),姜彤,任國(guó)玉,等.長(zhǎng)江流域1960—2004年極端強(qiáng)降水時(shí)空變化趨勢(shì)[J].氣候變化研究進(jìn)展,2006,2(1):9-14.

    VariationCharacteristicsofTemperatureandPrecipitationintheBailongRiverBasinduring1961-2010

    ZHANG Xiao-xiao, ZHANG Yu, XU Hao-jie

    (CollegeofEarthandEnvironmentalSciences,LanzhouUniversity,Lanzhou730000,China)

    Based on temperature and precipitation datasets in the Bailong River basin during 1961—2010, the seasonal and inter-annual variation characteristics including tendency, fluctuation, mutation, period and persistence of temperature and precipitation were analyzed by using linear tendency and polynomial regression model, Mann-Kendall test, complex Morlet wavelet and R/S analysis. The results showed that the seasonal and annual mean temperature had a significant increasing trend during 1961—2010; the weak warming trend was detected before mid-1980s; the increasing rate of temperature was accelerated after a warming mutation in 1990s; the annual precipitation mainly experienced ‘high—low—high—low’ fluctuations. The seasonal precipitation had similar fluctuations except for the winter. The precipitation increased weakly in the spring and winter, while decreased insignificantly in the summer, autumn and year. The temperature and precipitation had periodic variation of 7, 11, 22 a.The potential trend of seasonal and annual mean temperature would persistent. The precipitation in the summer and winter will increase, but decrease in the spring, autumn and year in the near future.

    Bailong River; climate change; tendency; mutation; period

    2013-09-12

    :2013-10-03

    國(guó)家自然科學(xué)基金項(xiàng)目“黑河上游土壤水文異質(zhì)性觀測(cè)試驗(yàn)及其對(duì)山區(qū)水文過(guò)程的影響”(91125010)

    張曉曉(1988—),女,黑龍江齊齊哈爾人,碩士研究生,主要研究方向?yàn)楣こ趟膶W(xué)。E-mail:xxzhang2011@lzu.edu.cn

    張鈺(1963—),男,甘肅平?jīng)鋈?博士,副教授,主要研究方向?yàn)樗Y源規(guī)劃論證評(píng)價(jià)和水土保持。E-mail:gszhangyu@126.com

    P467

    :A

    :1005-3409(2014)04-0238-08

    猜你喜歡
    白龍江年際平均氣溫
    An Annotated List of Lizards (Sauria:Squamata) Recorded from the People’s Republic of China
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    烏蘭縣近38年氣溫變化特征分析
    白龍江畔菜花香
    從全球氣候變暖大背景看萊州市30a氣溫變化
    1981—2010年拐子湖地區(qū)氣溫變化特征及趨勢(shì)分析
    近50年來(lái)全球背景下青藏高原氣候變化特征分析
    白龍江畔
    白龍江流域生物多樣性保護(hù)面臨的主要威脅及工作對(duì)策
    亞洲夏季風(fēng)的年際和年代際變化及其未來(lái)預(yù)測(cè)
    18禁观看日本| 国产一区二区三区在线臀色熟女 | 飞空精品影院首页| 韩国精品一区二区三区| 精品人妻在线不人妻| 一本一本久久a久久精品综合妖精| 宅男免费午夜| 久久精品国产亚洲av高清一级| 夜夜爽天天搞| 99久久精品国产亚洲精品| 高清av免费在线| 精品久久久久久久久久免费视频 | 亚洲精品国产精品久久久不卡| 日韩视频一区二区在线观看| 国产精品偷伦视频观看了| aaaaa片日本免费| 欧美日韩亚洲国产一区二区在线观看 | 在线免费观看的www视频| 成年人免费黄色播放视频| 精品久久久久久久毛片微露脸| 久久精品熟女亚洲av麻豆精品| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利影视在线免费观看| 成人av一区二区三区在线看| av在线播放免费不卡| 欧美黑人欧美精品刺激| 欧美精品人与动牲交sv欧美| bbb黄色大片| 成人三级做爰电影| 中文字幕高清在线视频| 热re99久久精品国产66热6| 12—13女人毛片做爰片一| 午夜免费成人在线视频| 亚洲精品久久成人aⅴ小说| 丰满的人妻完整版| 久久天堂一区二区三区四区| 亚洲av电影在线进入| 51午夜福利影视在线观看| 亚洲精华国产精华精| 欧美日韩视频精品一区| √禁漫天堂资源中文www| 99国产精品99久久久久| 国产精品1区2区在线观看. | 成人av一区二区三区在线看| 精品视频人人做人人爽| 成在线人永久免费视频| 国产av又大| 国产成人欧美在线观看 | 91老司机精品| 欧美精品av麻豆av| 免费日韩欧美在线观看| 免费在线观看影片大全网站| 亚洲精品一二三| 国内毛片毛片毛片毛片毛片| 免费在线观看黄色视频的| 国产精品永久免费网站| 啦啦啦视频在线资源免费观看| 又黄又爽又免费观看的视频| 久久精品国产a三级三级三级| 欧美黑人精品巨大| 亚洲av电影在线进入| 这个男人来自地球电影免费观看| av网站在线播放免费| 欧美最黄视频在线播放免费 | 国产精品亚洲av一区麻豆| 久久影院123| 国产精品一区二区在线不卡| 亚洲午夜精品一区,二区,三区| 三上悠亚av全集在线观看| 精品国产一区二区三区四区第35| 日本黄色日本黄色录像| 免费女性裸体啪啪无遮挡网站| 国产又爽黄色视频| 亚洲av日韩精品久久久久久密| 大型av网站在线播放| 国产色视频综合| av片东京热男人的天堂| av天堂久久9| 女人爽到高潮嗷嗷叫在线视频| 日韩中文字幕欧美一区二区| 久久国产乱子伦精品免费另类| 淫妇啪啪啪对白视频| 亚洲人成伊人成综合网2020| 亚洲午夜理论影院| 国产精品亚洲一级av第二区| 捣出白浆h1v1| 国产欧美日韩精品亚洲av| 国产成人影院久久av| 狂野欧美激情性xxxx| 王馨瑶露胸无遮挡在线观看| av中文乱码字幕在线| 国产av精品麻豆| 国产麻豆69| 人妻丰满熟妇av一区二区三区 | 99久久国产精品久久久| 久久青草综合色| 建设人人有责人人尽责人人享有的| 国产不卡av网站在线观看| 在线观看午夜福利视频| 50天的宝宝边吃奶边哭怎么回事| 好男人电影高清在线观看| 美女高潮到喷水免费观看| 久久香蕉激情| 国产亚洲精品一区二区www | 久久青草综合色| 三上悠亚av全集在线观看| 99re6热这里在线精品视频| a在线观看视频网站| 亚洲情色 制服丝袜| 精品国产超薄肉色丝袜足j| 久久国产精品男人的天堂亚洲| 男女之事视频高清在线观看| 在线播放国产精品三级| 婷婷成人精品国产| 黄色视频,在线免费观看| 国产精品亚洲av一区麻豆| 美国免费a级毛片| 国产又爽黄色视频| 99riav亚洲国产免费| 久久人人爽av亚洲精品天堂| 国产精品影院久久| 国产精品久久久久久人妻精品电影| 午夜久久久在线观看| 亚洲专区国产一区二区| 中文字幕人妻丝袜一区二区| 国产亚洲欧美98| 中文字幕精品免费在线观看视频| 成年动漫av网址| 亚洲aⅴ乱码一区二区在线播放 | 精品一区二区三区视频在线观看免费 | 国产精品一区二区免费欧美| 好男人电影高清在线观看| 叶爱在线成人免费视频播放| 国产片内射在线| 国产99白浆流出| 亚洲精品粉嫩美女一区| 亚洲综合色网址| 国产精品免费视频内射| 视频区欧美日本亚洲| 国产精品亚洲av一区麻豆| 老司机在亚洲福利影院| 中文字幕人妻丝袜一区二区| 精品久久蜜臀av无| 超碰97精品在线观看| 视频区欧美日本亚洲| 国产成人欧美在线观看 | 国产av又大| 啦啦啦视频在线资源免费观看| 人人妻,人人澡人人爽秒播| 男女之事视频高清在线观看| 久久久久久久久免费视频了| 国产男靠女视频免费网站| 男女免费视频国产| 国产99久久九九免费精品| 亚洲av成人不卡在线观看播放网| 久久久久久久久久久久大奶| 精品高清国产在线一区| 国产蜜桃级精品一区二区三区 | 十八禁高潮呻吟视频| 亚洲国产精品一区二区三区在线| 免费观看精品视频网站| 天天影视国产精品| 后天国语完整版免费观看| 亚洲av片天天在线观看| 久久精品国产亚洲av香蕉五月 | 久久国产乱子伦精品免费另类| 母亲3免费完整高清在线观看| 脱女人内裤的视频| 巨乳人妻的诱惑在线观看| 一级,二级,三级黄色视频| 免费在线观看日本一区| 久久精品熟女亚洲av麻豆精品| 十八禁人妻一区二区| 啦啦啦在线免费观看视频4| 极品少妇高潮喷水抽搐| 18在线观看网站| 啦啦啦视频在线资源免费观看| 女人高潮潮喷娇喘18禁视频| 中文欧美无线码| 在线国产一区二区在线| 巨乳人妻的诱惑在线观看| 亚洲人成电影观看| 免费在线观看视频国产中文字幕亚洲| www.999成人在线观看| 日本vs欧美在线观看视频| 欧美人与性动交α欧美软件| 亚洲精品一二三| 国产极品粉嫩免费观看在线| 97人妻天天添夜夜摸| 深夜精品福利| 狂野欧美激情性xxxx| 亚洲中文字幕日韩| 亚洲精品中文字幕一二三四区| 最新美女视频免费是黄的| 亚洲全国av大片| 久久人妻熟女aⅴ| 亚洲精品国产色婷婷电影| 每晚都被弄得嗷嗷叫到高潮| 男女免费视频国产| 日日摸夜夜添夜夜添小说| 91成人精品电影| 久久天躁狠狠躁夜夜2o2o| 免费在线观看日本一区| 99re在线观看精品视频| 在线观看日韩欧美| 国产乱人伦免费视频| 深夜精品福利| 亚洲avbb在线观看| 亚洲 国产 在线| 久久精品国产亚洲av高清一级| 欧美大码av| 婷婷精品国产亚洲av在线 | 999久久久精品免费观看国产| 国产一区二区三区在线臀色熟女 | 少妇 在线观看| 欧美激情极品国产一区二区三区| 欧美国产精品一级二级三级| av线在线观看网站| 欧美亚洲 丝袜 人妻 在线| 成人精品一区二区免费| 亚洲成人国产一区在线观看| 在线av久久热| 12—13女人毛片做爰片一| 最近最新中文字幕大全免费视频| 日日摸夜夜添夜夜添小说| 欧美在线黄色| 一进一出抽搐动态| 亚洲欧美激情在线| 久久精品国产亚洲av高清一级| 国产精品一区二区免费欧美| 国产真人三级小视频在线观看| 中文欧美无线码| 欧美乱色亚洲激情| 亚洲精品粉嫩美女一区| 久久国产精品人妻蜜桃| 免费少妇av软件| 久久九九热精品免费| 黑人操中国人逼视频| 少妇粗大呻吟视频| 亚洲精品中文字幕一二三四区| 国产精品久久电影中文字幕 | 久久久久久久午夜电影 | 国产精品免费视频内射| 免费女性裸体啪啪无遮挡网站| 搡老乐熟女国产| 亚洲精品在线观看二区| 中文字幕人妻熟女乱码| 麻豆乱淫一区二区| 久久亚洲真实| 一本综合久久免费| 午夜福利免费观看在线| 亚洲欧美日韩高清在线视频| 精品国产一区二区三区四区第35| 99久久精品国产亚洲精品| 在线观看免费视频日本深夜| 亚洲av第一区精品v没综合| 国产激情欧美一区二区| 宅男免费午夜| 久久久国产一区二区| 亚洲专区中文字幕在线| 久久性视频一级片| 777久久人妻少妇嫩草av网站| 国产免费现黄频在线看| 亚洲精品久久午夜乱码| 他把我摸到了高潮在线观看| 亚洲精品国产一区二区精华液| 欧美+亚洲+日韩+国产| 国产精品 国内视频| 精品高清国产在线一区| 欧美日本中文国产一区发布| 国产aⅴ精品一区二区三区波| 色在线成人网| x7x7x7水蜜桃| 欧美 日韩 精品 国产| 一级,二级,三级黄色视频| 一本综合久久免费| 国产精品一区二区在线观看99| 伦理电影免费视频| 午夜福利一区二区在线看| 亚洲人成伊人成综合网2020| 一进一出抽搐gif免费好疼 | 免费在线观看日本一区| 亚洲av美国av| 亚洲欧洲精品一区二区精品久久久| videos熟女内射| 久久午夜亚洲精品久久| 麻豆成人av在线观看| 少妇裸体淫交视频免费看高清 | 不卡av一区二区三区| 99在线人妻在线中文字幕 | 久久久国产精品麻豆| 国产高清国产精品国产三级| 国产主播在线观看一区二区| 国产亚洲欧美精品永久| 亚洲成人免费av在线播放| 在线观看日韩欧美| 国产精品综合久久久久久久免费 | 国产91精品成人一区二区三区| 最新在线观看一区二区三区| 嫩草影视91久久| 国产色视频综合| 热99国产精品久久久久久7| 欧美激情极品国产一区二区三区| 夜夜躁狠狠躁天天躁| 人妻 亚洲 视频| 国产精品亚洲av一区麻豆| 嫩草影视91久久| 久久午夜亚洲精品久久| 激情视频va一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 80岁老熟妇乱子伦牲交| 亚洲精品乱久久久久久| 涩涩av久久男人的天堂| 国产国语露脸激情在线看| 又紧又爽又黄一区二区| 久久99一区二区三区| 色综合婷婷激情| 欧美黄色片欧美黄色片| 欧美中文综合在线视频| 国产精品久久久av美女十八| 男女之事视频高清在线观看| 国产成+人综合+亚洲专区| 一二三四社区在线视频社区8| 日韩欧美在线二视频 | 国产片内射在线| 国产精品久久久av美女十八| 国产区一区二久久| 好男人电影高清在线观看| 久久久国产成人免费| 在线观看免费视频日本深夜| 中亚洲国语对白在线视频| 欧美日韩瑟瑟在线播放| 大型av网站在线播放| 国产成人精品久久二区二区免费| 大香蕉久久网| 亚洲 欧美一区二区三区| 午夜福利,免费看| 欧美中文综合在线视频| 日本五十路高清| 麻豆国产av国片精品| 精品少妇一区二区三区视频日本电影| 中国美女看黄片| av免费在线观看网站| 国产亚洲精品久久久久5区| 男男h啪啪无遮挡| 亚洲av美国av| 脱女人内裤的视频| 一二三四社区在线视频社区8| 国产高清videossex| 国产淫语在线视频| 美女扒开内裤让男人捅视频| 精品视频人人做人人爽| 可以免费在线观看a视频的电影网站| 婷婷精品国产亚洲av在线 | 日韩精品免费视频一区二区三区| 久久久国产一区二区| 国产真人三级小视频在线观看| 两人在一起打扑克的视频| 每晚都被弄得嗷嗷叫到高潮| 免费观看a级毛片全部| 精品国产亚洲在线| 青草久久国产| 成人特级黄色片久久久久久久| 日韩人妻精品一区2区三区| 女性被躁到高潮视频| 黄频高清免费视频| 一本综合久久免费| 久久人妻福利社区极品人妻图片| 一本综合久久免费| 99国产精品一区二区蜜桃av | 久久久久国内视频| 麻豆乱淫一区二区| 成年人黄色毛片网站| 在线观看66精品国产| 日韩人妻精品一区2区三区| 亚洲欧美精品综合一区二区三区| 国产精品 国内视频| 欧美激情久久久久久爽电影 | 亚洲av日韩精品久久久久久密| 精品久久久久久久久久免费视频 | 亚洲熟妇熟女久久| 日日爽夜夜爽网站| 日韩三级视频一区二区三区| 亚洲 欧美一区二区三区| 亚洲精品美女久久久久99蜜臀| 亚洲欧美日韩高清在线视频| 日本黄色日本黄色录像| 亚洲精品av麻豆狂野| av有码第一页| 一级毛片高清免费大全| 久久久久久久午夜电影 | av中文乱码字幕在线| 亚洲欧美激情综合另类| 新久久久久国产一级毛片| 色婷婷av一区二区三区视频| av线在线观看网站| 成人国语在线视频| 精品国产一区二区三区四区第35| 欧美国产精品一级二级三级| 一进一出好大好爽视频| 亚洲在线自拍视频| av有码第一页| 看黄色毛片网站| 成人免费观看视频高清| 美国免费a级毛片| 精品国产一区二区三区久久久樱花| 欧美乱色亚洲激情| 高清黄色对白视频在线免费看| 无限看片的www在线观看| 狠狠婷婷综合久久久久久88av| 欧美 亚洲 国产 日韩一| 午夜亚洲福利在线播放| 日本五十路高清| 国产精品美女特级片免费视频播放器 | 91大片在线观看| 在线观看免费视频网站a站| 人妻久久中文字幕网| 久久久精品免费免费高清| 法律面前人人平等表现在哪些方面| 国产xxxxx性猛交| 亚洲精品乱久久久久久| 一区福利在线观看| 黑人猛操日本美女一级片| 国产蜜桃级精品一区二区三区 | 19禁男女啪啪无遮挡网站| 日日夜夜操网爽| 国产精品九九99| 亚洲 欧美一区二区三区| 女性生殖器流出的白浆| 亚洲专区字幕在线| 免费人成视频x8x8入口观看| 曰老女人黄片| 国产精品 国内视频| 国产视频一区二区在线看| 久久九九热精品免费| 久久婷婷成人综合色麻豆| 别揉我奶头~嗯~啊~动态视频| 老汉色∧v一级毛片| 国产极品粉嫩免费观看在线| 国产精品久久久av美女十八| 免费久久久久久久精品成人欧美视频| 久久精品国产亚洲av香蕉五月 | 国内毛片毛片毛片毛片毛片| 99国产精品一区二区三区| 欧美精品人与动牲交sv欧美| 久久香蕉精品热| 一级作爱视频免费观看| 免费看a级黄色片| 欧美一级毛片孕妇| 国产精品香港三级国产av潘金莲| 最新美女视频免费是黄的| 欧美激情极品国产一区二区三区| 美国免费a级毛片| 人成视频在线观看免费观看| 国产亚洲一区二区精品| 黄色片一级片一级黄色片| 国产一区有黄有色的免费视频| 国产激情久久老熟女| 亚洲第一青青草原| 国产激情欧美一区二区| 91av网站免费观看| 亚洲精品一卡2卡三卡4卡5卡| 大香蕉久久成人网| 日韩欧美国产一区二区入口| 久久久久视频综合| 亚洲,欧美精品.| 一级毛片精品| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品一卡2卡三卡4卡5卡| 又大又爽又粗| 伊人久久大香线蕉亚洲五| 精品无人区乱码1区二区| 99国产极品粉嫩在线观看| 老司机午夜福利在线观看视频| 久久精品aⅴ一区二区三区四区| 成人18禁高潮啪啪吃奶动态图| 精品人妻熟女毛片av久久网站| 黄色毛片三级朝国网站| 高清黄色对白视频在线免费看| 一边摸一边做爽爽视频免费| 成人黄色视频免费在线看| 韩国精品一区二区三区| 婷婷丁香在线五月| 午夜福利欧美成人| www.熟女人妻精品国产| 久久香蕉国产精品| 在线十欧美十亚洲十日本专区| 交换朋友夫妻互换小说| 成人特级黄色片久久久久久久| 国产成人一区二区三区免费视频网站| 在线视频色国产色| www.精华液| 人人妻人人爽人人添夜夜欢视频| 国产亚洲精品一区二区www | 成人特级黄色片久久久久久久| 在线av久久热| 久热这里只有精品99| 精品国内亚洲2022精品成人 | 高潮久久久久久久久久久不卡| videosex国产| 又黄又爽又免费观看的视频| 国产色视频综合| www.自偷自拍.com| 婷婷丁香在线五月| 日本精品一区二区三区蜜桃| 国产在视频线精品| 国产精品久久视频播放| 国产欧美日韩一区二区三区在线| 高清av免费在线| 一进一出好大好爽视频| 黄色毛片三级朝国网站| 国产精品久久久久久人妻精品电影| 黄网站色视频无遮挡免费观看| 99精品久久久久人妻精品| 精品久久久久久,| 亚洲国产中文字幕在线视频| 99热只有精品国产| 在线观看舔阴道视频| 久久天躁狠狠躁夜夜2o2o| 中文字幕制服av| 两人在一起打扑克的视频| 成年动漫av网址| 午夜免费鲁丝| 性色av乱码一区二区三区2| 精品福利观看| 久9热在线精品视频| videos熟女内射| 免费av中文字幕在线| 一二三四社区在线视频社区8| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品美女特级片免费视频播放器 | 欧美乱色亚洲激情| 亚洲国产欧美一区二区综合| 午夜视频精品福利| 免费久久久久久久精品成人欧美视频| 欧美黑人欧美精品刺激| 高清欧美精品videossex| 一区二区三区精品91| 99久久人妻综合| 男女高潮啪啪啪动态图| 亚洲性夜色夜夜综合| 精品福利观看| bbb黄色大片| av视频免费观看在线观看| 高潮久久久久久久久久久不卡| 国产亚洲精品久久久久5区| 国产无遮挡羞羞视频在线观看| 精品久久久久久电影网| 久久精品国产亚洲av香蕉五月 | 人成视频在线观看免费观看| 老汉色∧v一级毛片| 亚洲第一欧美日韩一区二区三区| 日本黄色日本黄色录像| 成年人午夜在线观看视频| 国产精华一区二区三区| 日韩成人在线观看一区二区三区| 人人妻人人添人人爽欧美一区卜| 国内毛片毛片毛片毛片毛片| 天堂俺去俺来也www色官网| 亚洲免费av在线视频| 99热只有精品国产| 亚洲自偷自拍图片 自拍| 一区在线观看完整版| 亚洲美女黄片视频| 淫妇啪啪啪对白视频| 最近最新中文字幕大全免费视频| 妹子高潮喷水视频| 18禁裸乳无遮挡动漫免费视频| 欧美黄色片欧美黄色片| 午夜福利影视在线免费观看| 中国美女看黄片| 人人妻,人人澡人人爽秒播| 精品第一国产精品| 天天添夜夜摸| av免费在线观看网站| 欧美丝袜亚洲另类 | 大型黄色视频在线免费观看| 91麻豆av在线| 哪里可以看免费的av片| 免费av不卡在线播放| 三级毛片av免费| 男人舔女人下体高潮全视频| 色综合婷婷激情| 成人无遮挡网站| 伊人久久大香线蕉亚洲五| 最新美女视频免费是黄的| 观看美女的网站| 国产亚洲精品久久久久久毛片| 成年免费大片在线观看| 黄片小视频在线播放| 国产熟女xx| www.www免费av| 国产野战对白在线观看| 国产在线精品亚洲第一网站| 日日干狠狠操夜夜爽| 97碰自拍视频| 精品久久久久久久久久久久久| 国产一区二区三区在线臀色熟女| 国产99白浆流出| 99久久综合精品五月天人人| 亚洲精品粉嫩美女一区| 九九在线视频观看精品| 亚洲欧美激情综合另类| 国产成年人精品一区二区| 国产精品国产高清国产av| 熟女人妻精品中文字幕| 99久久精品热视频| 国产免费男女视频| 欧美色视频一区免费| 国产伦精品一区二区三区视频9 | 99在线人妻在线中文字幕| 国内精品一区二区在线观看| 黄色片一级片一级黄色片| 岛国视频午夜一区免费看| 两性午夜刺激爽爽歪歪视频在线观看|