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

    1960-2014年祁連山南坡及其附近地區(qū)降水時(shí)空變化特征

    2018-06-15 06:44:58付建新曹廣超李玲琴曹生奎唐仲霞楊曉敏刁二龍
    水土保持研究 2018年4期
    關(guān)鍵詞:南坡降水強(qiáng)度祁連山

    付建新, 曹廣超, 李玲琴, 曹生奎, 唐仲霞,楊曉敏, 蔣 剛, 虞 敏, 袁 杰, 刁二龍

    (1.青海師范大學(xué) 青海省自然地理與環(huán)境過(guò)程重點(diǎn)實(shí)驗(yàn)室, 西寧 810008; 2.青海師范大學(xué) 地理科學(xué)學(xué)院, 西寧 810008)

    IPCC報(bào)告表明全球變暖已經(jīng)是不容爭(zhēng)辯的事實(shí),全球變暖對(duì)降水有可能產(chǎn)生影響[1]。秦大河等[2]認(rèn)為全球變暖可能造成中緯度地區(qū)降水增加,前人對(duì)于祁連山區(qū)的降水研究主要從時(shí)間尺度與空間尺度進(jìn)行了分析。時(shí)間尺度包括小時(shí)、季節(jié)、年等角度。小時(shí)尺度:殷雪蓮等[3]對(duì)祁連山區(qū)2003年7月25日與2005年6月30日兩次降水過(guò)程從環(huán)流背景、動(dòng)力機(jī)制以及地形特征等方面做了研究。季節(jié)尺度:劉雪梅等[4]對(duì)祁連山區(qū)夏季降水的日變化特征進(jìn)行了分析;陳乾等[5]對(duì)祁連山區(qū)2006年夏季降水過(guò)程從3類(lèi)氣流型進(jìn)行了分析。年尺度:強(qiáng)芳等[6]基于格點(diǎn)數(shù)據(jù)對(duì)祁連山區(qū)1961—2012年面雨量的特征進(jìn)行了分析;魏鋒等[7]使用小波變換方法分析了祁連山地區(qū)17個(gè)測(cè)站1960—2004年5—9月逐日降水量的特征;張小明等[8]對(duì)祁連山近45 a異常降水的時(shí)空特征進(jìn)行了研究;張存杰等[9]采用EOF和REOF等方法,研究了近40 a來(lái)祁連山附近氣溫和降水的時(shí)空分布特征。空間尺度包括流域、分段、整個(gè)區(qū)域。流域尺度:藍(lán)永超等[10]研究了黑河山區(qū)匯流區(qū)降水對(duì)全球變暖的響應(yīng);牛赟等[11]對(duì)祁連山大野口流域降水特征進(jìn)行了分析。分段尺度:王寧練等[12]對(duì)祁連山中段北坡最大降水高度帶進(jìn)行了分析。整個(gè)區(qū)域:賈文雄[13]對(duì)近50 a來(lái)祁連山及河西走廊降水的時(shí)空變化進(jìn)行了研究;湯懋蒼[14]對(duì)祁連山區(qū)降水的地理分布特征進(jìn)行了闡述。對(duì)降水成因的研究主要從經(jīng)緯度、地形、大氣環(huán)流等角度。陳少勇等[15]認(rèn)為祁連山主體的降水多于周?chē)貐^(qū),祁連山東部降水比西部穩(wěn)定;李巖瑛等[16]分析了不同降水強(qiáng)度的時(shí)空分布特征及其與海拔的關(guān)系;賈文雄[17]認(rèn)為祁連山的氣溫與地理位置的關(guān)系存在相關(guān)性關(guān)系,降水與地理位置的關(guān)系較為復(fù)雜;程鵬等[18]分析了高空西風(fēng)急流對(duì)祁連山區(qū)降水的影響。也有學(xué)者從降水與祁連山區(qū)水源涵養(yǎng)的關(guān)系角度進(jìn)行了闡述。常學(xué)向等[19]研究了祁連山林區(qū)森林對(duì)降水的截留作用;張學(xué)龍等[20]分析了祁連山青海云杉林截留對(duì)降水的分配效應(yīng)。前人對(duì)于祁連山區(qū)降水主要集中在時(shí)空視角、成因分析與對(duì)水源涵養(yǎng)的影響等方面,在方法上使用傳統(tǒng)與信息技術(shù)結(jié)合的較少,本文主要對(duì)祁連山南坡附近的降水特征進(jìn)行分析,使用傳統(tǒng)氣象分析方法與GIS與SPPSS軟件制圖分析方法并用進(jìn)行研究。

    祁連山區(qū)的生態(tài)是國(guó)家當(dāng)前與今后重點(diǎn)保護(hù)與建設(shè)的對(duì)象,該區(qū)域是“絲綢之路經(jīng)濟(jì)帶”重要的水源涵養(yǎng)功能區(qū)[21],所以在新形勢(shì)下,研究祁連山南坡及其附近地區(qū)降水的變化具有重要的理論與現(xiàn)實(shí)意義。

    1 研究區(qū)概況

    本文研究區(qū)概況見(jiàn)圖1,核心區(qū)位于青海省境內(nèi),但是核心區(qū)氣象站點(diǎn)較少,所以增加了外圍區(qū)域的氣象站點(diǎn),外圍區(qū)為核心區(qū)的南北兩側(cè),分別位于青、甘兩省。祁連山南坡,海拔2 257~5 235 m,是我國(guó)重要的西北—東南走向的一系列山脈,是我國(guó)重要的地理分界線之一[22],具有高原大陸性氣候特征[23]。祁連山南坡是青海省重要的水源涵養(yǎng)保護(hù)區(qū),氣候的變化直接會(huì)對(duì)當(dāng)?shù)刂脖坏壬鷳B(tài)環(huán)境造成影響,通過(guò)研究祁連山南坡及其附近地區(qū)降水的變化為進(jìn)一步研究祁連山南坡水源涵養(yǎng)區(qū)的生態(tài)系統(tǒng)等研究做鋪墊。

    圖1祁連山南坡及其附近地區(qū)氣象站點(diǎn)分布

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

    本文數(shù)據(jù)來(lái)源于中國(guó)氣象數(shù)據(jù)網(wǎng)(http:∥data.cma.cn/),選取1960—2014年祁連山南坡及其附近地區(qū)19個(gè)氣象站點(diǎn)的降水日數(shù)、降水量、降水強(qiáng)度數(shù)據(jù),為了保證研究結(jié)果的準(zhǔn)確性,所以對(duì)數(shù)據(jù)進(jìn)行一致性檢驗(yàn)的驗(yàn)證,對(duì)區(qū)域內(nèi)平均降水序列與各站點(diǎn)降水的原始序列的相關(guān)性進(jìn)行了檢測(cè),相關(guān)系數(shù)均在0.793以上,說(shuō)明降水?dāng)?shù)據(jù)表現(xiàn)出較高的相關(guān)性,所以數(shù)據(jù)不僅能反映各站點(diǎn)的氣候變化,同時(shí)也能反映出區(qū)域內(nèi)的氣候變化;降水日數(shù)為日降水量≥0.1 mm日數(shù)的總和,降水強(qiáng)度為各等級(jí)降水總量與其對(duì)應(yīng)降水日數(shù)之比。研究方法包括線性趨勢(shì)法、相關(guān)分析法、多項(xiàng)式趨勢(shì)法、5 a滑動(dòng)平均,Mann-Kendall與滑動(dòng)t檢驗(yàn)[24-25]法,利用ArcGIS 10.0對(duì)降水日數(shù)與降水強(qiáng)度的空間變化做制圖分析,對(duì)于文中顯著性檢驗(yàn),使用SPSS統(tǒng)計(jì)學(xué)軟件中的相關(guān)性分析模塊,對(duì)時(shí)間序列與原序列變量之間的相關(guān)系數(shù)進(jìn)行不同程度的顯著性檢驗(yàn)(p<0.05,p<0.01)。

    3 降水的時(shí)間序列變化

    3.1 多年降水日數(shù)與降水強(qiáng)度年際變化趨勢(shì)

    祁連山區(qū)的降水受到3種氣流的影響:其一來(lái)自太平洋的東南季風(fēng)帶來(lái)的濕潤(rùn)水汽,其二受到來(lái)自大西洋冷濕氣流的影響,夏季還要受到翻越高大青藏高原之后的印度洋暖濕氣團(tuán)的影響[26]。

    祁連山南坡及其附近地區(qū)1960—2014年多年降水日數(shù)(圖2A)為緩慢增長(zhǎng)趨勢(shì)(0.030 d/a),其平均值為77.67 d,降水日數(shù)大于平均值的年份占49.09%,小于平均值的年份占50.91%;降水日數(shù)最多的為93.42 d(1988年),最少的為64.89 d(1997年),二者相差28.53 d;90 d以上占5.45%,80 d以上占32.73%;從5 a滑動(dòng)平均來(lái)看,多年降水日數(shù)20世紀(jì)90年代以后較以前變化大,大致分為2個(gè)大周期(20世紀(jì)80年代中后期—21世紀(jì)初期、21世紀(jì)初期以來(lái))與3個(gè)小周期(20世紀(jì)60年代中期—70年代中期、70年代中期—80年代初期、80年代初期—后期)。降水強(qiáng)度(圖2B)整體上為緩慢波動(dòng)增長(zhǎng)趨勢(shì)(0.009 mm/a),降水強(qiáng)度的平均值為3.247 mm/a,52.73%的年份降水強(qiáng)度超過(guò)了均值,47.27%的年份降水強(qiáng)度低于均值,最大值為3.79mm/a(2014年),最小值為2.63 mm/a(1962年)二者相差1.16 mm/a;從5 a滑動(dòng)平均來(lái)看,其波動(dòng)較降水日數(shù)小,期間存在2個(gè)大的周期(20世紀(jì)60年代中期—70年代中期、90年代初期—21世紀(jì)初期),同時(shí)還并存若干小的周期。

    圖2祁連山南坡及其附近地區(qū)多年平均降水日數(shù)與降水強(qiáng)度時(shí)間變化

    3.2 降水日數(shù)與降水強(qiáng)度的季節(jié)變化趨勢(shì)

    從四季降水的線性回歸圖(圖3)得出,春季(3月—5月)、夏季(6月—8月)、秋季(9月—11月)、冬季(12月—翌年2月)的降水日數(shù)傾斜率分別為:-0.000 2,0.002 3,-0.000 04,0.008 2 d/a,除了冬季降水日數(shù)為增加趨勢(shì),其余季節(jié)變化趨勢(shì)不大,四季的降水日數(shù)的增加趨勢(shì)均小于多年平均降水日數(shù)(0.034 9 d/a);降水強(qiáng)度傾斜率分別為:0.005 5 mm/a(p<0.01),0.013 6 mm/a(p<0.01),0.007 2 mm/a(p<0.05),0.005 mm/a(p<0.01),均通過(guò)了顯著性檢驗(yàn),四季降水強(qiáng)度均為增加趨勢(shì),上升的速率由大到小依次為:夏、秋、春、冬,夏季降水強(qiáng)度的增長(zhǎng)率高于多年平均降水強(qiáng)度的增長(zhǎng)率(0.008 8 mm/a),說(shuō)明夏季降水強(qiáng)度對(duì)全年降水貢獻(xiàn)最大。

    春季降水日數(shù)(圖3A)中最大值為8.404 d(1967年),最小值為3.737 d(2013年),二者相差4.667 d,其中大于平均值(5.949 d)的占到47.27%,小于平均值的占到52.73%;5 a滑動(dòng)平均顯示3個(gè)升降周期:20世紀(jì)60年代中期—70年代末、70年代末—21世紀(jì)初期、21世紀(jì)初期—研究期末,其中第2個(gè)周期最長(zhǎng),包括了1個(gè)顯著上升區(qū)(20世紀(jì)70年代末—90年代初)。夏季降水日數(shù)(圖3C)整體變化與春季相似,變化不明顯,不同的是夏季略微上升,春季略微下降;最大值14.754 d(1983年)與最小值9.702 d(1962年)相差5.052 d,其中大于平均值(11.908 d)的占到45.45%,;20世紀(jì)80年代中期為1個(gè)降水日數(shù)高峰期。秋季降水日數(shù)(圖3E)變化趨勢(shì)基本不變,平均值為6.383 d,大于與小于平均值所占比重大致相當(dāng),基本上在平均值上下波動(dòng),20世紀(jì)70年代中期—90年代中期為1個(gè)較長(zhǎng)下降區(qū)間,20世紀(jì)90年代中期—21世紀(jì)10年代初期為1個(gè)較長(zhǎng)的上升區(qū)間。冬季降水日數(shù)(圖3G)的變化趨勢(shì)較其他季節(jié)的變化大,呈現(xiàn)增長(zhǎng)趨勢(shì)(0.008 2 d/a),2008年降水日數(shù)最多(5.035 d),比最少降水日數(shù)(1.456 d)多3.579 d,小于平均值的比重比大于平均值的比重多16.36%,5 a滑動(dòng)平均表現(xiàn)4個(gè)升降區(qū)間。

    春季降水強(qiáng)度(圖3B)為上升趨勢(shì)(0.005 5 mm/a),春季降水強(qiáng)度的平均值為2.521 mm/a,41.82%的年份變化趨勢(shì)大于均值,58.18%的年份變化趨勢(shì)小于均值,最大值3.799 mm/a(1985年)與最小值差1.440 mm/a(1995年)差值較大,為2.359 mm/a;5 a平均滑動(dòng)顯示兩個(gè)大的周期,即20世紀(jì)60年代中期—80年代中期與20世紀(jì)80年代中期—21世紀(jì)初。夏季降水強(qiáng)度(圖3D)的傾斜率(0.013 6 mm/a)比春季、秋季、冬季降水強(qiáng)度的傾斜率分別高0.008,0.006,0.009 mm/a;夏季降水強(qiáng)度范圍為3.343~5.014 mm/a,平均值為4.298 mm/a,大于均值的比重為52.73%,小于均值的比重為47.27%,大致從20世紀(jì)90年代降水強(qiáng)度開(kāi)始大于均值;5 a滑動(dòng)平均表現(xiàn)出5個(gè)周期,其中包括了多個(gè)小的波動(dòng)。秋季降水強(qiáng)度(圖3F)的增長(zhǎng)率為0.072 mm/a,平均值為3.196 mm/a,期間包括1個(gè)大的周期與若干個(gè)小的周期,大的周期為20世紀(jì)80年代—90年代中期。冬季降水強(qiáng)度(圖3H)的增長(zhǎng)率(0.005 mm/a)是四季中最低的,其均值為0.899 mm,大于均值的年份與小于均值的年份之比為5∶6,20世紀(jì)60年代中期—70年代中期為低于均值的1個(gè)較長(zhǎng)周期,21世紀(jì)開(kāi)始了新的1個(gè)較長(zhǎng)周期。

    圖3祁連山南坡及其附近地區(qū)降水日數(shù)與降水強(qiáng)度的季節(jié)變化

    3.3 季節(jié)降水量的突變檢驗(yàn)

    春季降水量(圖4A):正序列UF與反序列UB的交點(diǎn)是1982年,交點(diǎn)且位于±1.96臨界線之間,所以1982年為突變開(kāi)始的年份,這與賈文雄等[27]對(duì)祁連山東部季節(jié)降水突變檢驗(yàn)的結(jié)論一致,1986年降水量增加趨勢(shì)明顯,但是之后UF與UB基本沒(méi)有超過(guò)臨界線,UF只在1993年超過(guò)了臨界線,所以突變并不明顯;夏季降水量(圖4B):正序列曲線UF的值于1973年大于0,之后降水量波動(dòng)上升,期間UF與UB有有效的交點(diǎn)為3個(gè),1972年、2008年、2010年,結(jié)合滑動(dòng)t檢驗(yàn),當(dāng)n1=n2=5,10時(shí),在a=0.01的顯著性檢驗(yàn),可以確定1972年為突變年份開(kāi)始年份,1987—1991年、1993—2001年、2006—2009年、2010—2014年為突變時(shí)間段,一直到2010年突變明顯加快;秋季降水量(圖4C):突變開(kāi)始的年份為2003年,但是正序列曲線UF一直沒(méi)有超過(guò)臨界線,突變不明顯;冬季降水量(圖4D):正序列UF與反序列UB的交點(diǎn)共有4個(gè),為1989年、1993年、1996年、1998年,結(jié)合滑動(dòng)t檢驗(yàn),當(dāng)n1=n2=5,9時(shí),在a=0.05的顯著性檢驗(yàn),確定為1998年為突變開(kāi)始的年份,從1990年降水增多,尤其在2005年之后,增加趨勢(shì)迅速。

    圖4祁連山南坡及其附近地區(qū)四季降水量突變檢驗(yàn)

    3.4 降水日數(shù)與降水強(qiáng)度的年代際變化趨勢(shì)

    年代際年平均降水日數(shù)(圖5)由大到小依次為:81.03 d(20世紀(jì)80年代),79.97 d(2000—2014年),77.01 d(20世紀(jì)70年代),76.89 d(20世紀(jì)60年代),74.70 d(20世紀(jì)90年代),20世紀(jì)80年代與2000—2014年高于年降水日數(shù)平均值(77.67 d);20世紀(jì)80年代增長(zhǎng)最快,為1.071 d/a,20世紀(jì)90年代為負(fù)增長(zhǎng)(-0.326 d/a),二者相差1.397 d/a,年代際增長(zhǎng)率波動(dòng)較大。20世紀(jì)60,70,80,90年代、2000-2014年的平均降水強(qiáng)度分別3.055 mm,3.219 mm,3.144 mm,3.308 mm,3.448 mm,可以看出60—80年代降水強(qiáng)度較低,90年代之后降水強(qiáng)度增加明顯,尤其進(jìn)入21世紀(jì)降水強(qiáng)度增加愈加突出,這與中國(guó)60—80年代中期降水偏少,80年代中后期增多的趨勢(shì)相似[28];降水強(qiáng)度增長(zhǎng)率最高的是0.381 mm/10 a(20世紀(jì)90年代),增長(zhǎng)率最低的是0.074 mm/10 a(20世紀(jì)70年代),大致在0~0.4 mm/10 a。

    4 降水的空間變化

    4.1 多年平均降水日數(shù)與降水量空間分布

    從圖6中可以看出,多年平均降水量與平均降水日數(shù)均呈現(xiàn)帶狀分布,大致為西北—東南走向,由東南向西北遞減,主要因?yàn)橄募撅L(fēng)東南季風(fēng)對(duì)此區(qū)域影響由東南向西北減弱;研究區(qū)多年平均降水日數(shù)與降水量的均值分別為80.76 d與267.79 mm,降水日數(shù)排在前4位即降水日數(shù)大于100 d的站點(diǎn)為烏鞘嶺、門(mén)源、野牛溝、祁連、剛察,這些站點(diǎn)均位于東南部,降水量大于400 mm的站點(diǎn)為門(mén)源、野牛溝、祁連、烏鞘嶺,這4個(gè)站點(diǎn)同樣均位于東南部;位于西北部的高臺(tái)、酒泉降水日數(shù)大約分別為45,40 d,降水量大約分別為109,88 mm;祁連山南坡與北坡的降水量之比為4∶1,南坡多年降水量最多與最少的站點(diǎn)分別為門(mén)源(520.86 mm)、德令哈(180.39 mm),北坡多年降水量最多與最少的站點(diǎn)分別為永昌(204.29 mm)、酒泉(87.51 mm);南坡多年降水日數(shù)最多與最少的站點(diǎn)分別為烏鞘嶺(136.53 d/a)、德令哈(52.53 d/a),北坡多年降水日數(shù)最多與最少的站點(diǎn)分別為永昌(7.82 d/a)、民勤(39.46 d/a),南坡降水明顯多于北坡,因?yàn)槠钸B山南坡為東南季風(fēng)的迎風(fēng)坡,來(lái)自西北太平洋的暖濕氣流順著湟水谷地沿著祁連山南坡不斷爬升,受到地形的抬升作用,在一定海拔高度降水量增加,南坡形成了多雨區(qū),北坡處于東南季風(fēng)的背風(fēng)坡,水汽受到山地阻擋,成為雨影區(qū),氣候類(lèi)型為溫帶大陸性氣候,海洋水汽影響弱,降水量少。

    圖5祁連山南坡及其附近地區(qū)降水日數(shù)與降水強(qiáng)度的年代際變化

    圖6祁連山南坡及其附近地區(qū)多年平均降水量與平均降水日數(shù)空間分布

    4.2 降水日數(shù)與降水強(qiáng)度年際變化趨勢(shì)空間分布

    從(圖7A)可知,1960—2014年各站點(diǎn)的年降水日數(shù)為正增長(zhǎng)的站點(diǎn)數(shù)量為10個(gè),占到全部站點(diǎn)52.63%,為負(fù)增長(zhǎng)的站點(diǎn)數(shù)量為9個(gè),占到全部站點(diǎn)47.37%,正負(fù)增長(zhǎng)之比大體相當(dāng);增長(zhǎng)最快與最慢的站點(diǎn)分別為德令哈(0.365 d/a)與托勒(0.004 d/a),下降最快與最慢的站點(diǎn)分別為門(mén)源(-0.218 d/a)與西寧(-0.012 d/a);核心研究區(qū)傾斜率的平均值為-0.09 d/a,呈下降趨勢(shì),外圍研究區(qū)傾斜率的平均值為0.051 d/a,呈上升趨勢(shì);祁連山南坡站點(diǎn)傾斜率平均值為0.026 d/a,北坡站點(diǎn)傾斜率平均值為0.020 d/a,南北坡降水日數(shù)基本相當(dāng);只有3個(gè)站點(diǎn)通過(guò)了顯著性檢驗(yàn),德令哈為p<0.01,門(mén)源、烏鞘嶺為p<0.05。

    由圖7B可知,1960—2014年各站點(diǎn)的年降水強(qiáng)度除了民和(-0.023 mm/10 a)與恰卜恰(-0.002 mm/10 a)為負(fù)增長(zhǎng)外,其余站點(diǎn)均為正增長(zhǎng),占全部站點(diǎn)89.47%;增長(zhǎng)最快與最慢的站點(diǎn)分別為都蘭(0.207 mm/10 a)與張掖(0.037 mm/10 a),二者相差0.170 mm/10 a;傾斜率>0.1 mm/10 a的站點(diǎn)數(shù)量為8個(gè),占到全部站點(diǎn)42.11%,增長(zhǎng)率<0.1 mm/10 a的站點(diǎn)數(shù)量為11個(gè),占到全部站點(diǎn)57.89%;核心研究區(qū)3個(gè)站點(diǎn)傾斜率>0.1 mm/10 a,野牛溝(0.155 mm/10 a)、托勒(0.154 mm/10 a)、祁連(0.117 mm/10 a),傾斜率平均值為0.120 mm/10 a,外圍研究區(qū)傾斜率平均值為0.079 mm/10 a,核心研究區(qū)的增長(zhǎng)速度是外圍研究區(qū)的1.524倍,增長(zhǎng)較快;祁連山南坡站點(diǎn)傾斜率平均值為0.100 1 mm/10 a,北坡站點(diǎn)傾斜率平均值為0.066 mm/10 a,所以南坡增長(zhǎng)速度大約是北坡的1.513倍;63.16%的站點(diǎn)通過(guò)顯著性檢驗(yàn),其中p<0.01的站點(diǎn)為5個(gè),占全部站點(diǎn)26.32%,占通過(guò)檢驗(yàn)站點(diǎn)41.67%;p<0.05的站點(diǎn)為7個(gè),占全部站點(diǎn)36.84%,占通過(guò)檢驗(yàn)站點(diǎn)58.33%。

    圖7祁連山南坡及其附近地區(qū)降水日數(shù)與降水強(qiáng)度年際變化空間分布

    4.3 降水日數(shù)與降水強(qiáng)度的季節(jié)變化空間分布

    各站點(diǎn)的春季降水日數(shù)(圖8A)正負(fù)增長(zhǎng)之比為11∶8,分別占到全部站點(diǎn)57.89%與42.11%;正增長(zhǎng)的范圍為0.002~0.062 d/a,負(fù)增長(zhǎng)的范圍為-0.105 8~-0.000 4 d/a,增長(zhǎng)最快與最慢的站點(diǎn)分別為德令哈與高臺(tái),下降最快與最慢的站點(diǎn)分別為門(mén)源與張掖,增長(zhǎng)最快與下降最快的站點(diǎn)與年際變化的站點(diǎn)一致;核心研究區(qū)傾斜率的平均值為-0.049 d/a,呈下降趨勢(shì),外圍研究區(qū)傾斜率的平均值為0.008 d/a,呈略微上升趨勢(shì);祁連山南坡站點(diǎn)傾斜率平均值為-0.008 d/a,北坡站點(diǎn)傾斜率平均值為0.003 d/a,南北坡降水日數(shù)變化趨勢(shì)差別不大;門(mén)源、都蘭、剛察3個(gè)站點(diǎn)通過(guò)了顯著性檢驗(yàn)(p<0.05)。夏季降水日數(shù)(圖8C)正負(fù)增長(zhǎng)的站點(diǎn)數(shù)量基本相當(dāng),增長(zhǎng)最快的站點(diǎn)0.203 d/a(德令哈)與下降最快的站點(diǎn)-0.112 d/a(門(mén)源)相差0.315 d/a;核心研究區(qū)野牛溝、祁連、門(mén)源呈下降趨勢(shì),托勒呈增長(zhǎng)趨勢(shì),整體上呈下降趨勢(shì)(-0.025 d/a);外圍研究區(qū)整體上呈上升趨勢(shì)(0.023 d/a);祁連山南坡站點(diǎn)整體上呈上升趨勢(shì)(0.029 d/a),北坡整體上呈下降趨勢(shì)(-0.016 d/a);通過(guò)顯著性檢驗(yàn)的站點(diǎn)占到全部站點(diǎn)21.05%,其中德令哈通過(guò)(p<0.01)的檢驗(yàn)。秋季降水日數(shù)(圖8E)為正增長(zhǎng)的站點(diǎn)數(shù)量為10個(gè),負(fù)增長(zhǎng)的為9個(gè),整體上為略微下降趨勢(shì)(-0.005 d/a);變化最大的兩個(gè)站點(diǎn)為烏鞘嶺(0.079 d/a)與野牛溝(-0.072 d/a),二者相差0.151 d/a;只有烏鞘嶺通過(guò)了(p<0.05)相關(guān)性檢驗(yàn)。冬季降水日數(shù)(圖8G)的正增長(zhǎng)、零增長(zhǎng)、負(fù)增長(zhǎng)三者之比為15∶1∶3,整體上表現(xiàn)出正增長(zhǎng)(0.024 d/a),增長(zhǎng)最快的是烏鞘嶺(0.108 d/a),降低最快的是-0.034 d/a,二者相差0.074 d/a;核心研究區(qū)整體上為增長(zhǎng)趨勢(shì)(0.014 d/a),托勒是唯一零增長(zhǎng)的站點(diǎn);3個(gè)站點(diǎn)通過(guò)了顯著性檢驗(yàn)。

    各站點(diǎn)的春季(圖8B)、夏季(圖8D)、秋季(圖8F)、冬季(圖8H)(以下皆為此順序)降水強(qiáng)度增長(zhǎng)率之比分別為15∶4,18∶1,13∶6,18∶1,增長(zhǎng)率最高的站點(diǎn)分別為都蘭(0.251 mm/10 a)、永昌(0.252 mm/10 a)、烏鞘嶺(0.629 mm/10 a)、高臺(tái)(0.155 mm/10 a),下降最快的站點(diǎn)分別為茶卡(-0.074 mm/10 a)、民和(-0.102 mm/10 a)、酒泉(-0.095 mm/10 a)、托勒(-0.003 mm/10 a);核心研究區(qū)的平均增長(zhǎng)率分別為0.069,0.171,0.096,0.027 mm/10 a;通過(guò)顯著性檢驗(yàn)的站點(diǎn)數(shù)量依次為0,3,2,10;其中通過(guò)p<0.01檢驗(yàn)的站點(diǎn)數(shù)量為1(秋季),6(冬季),春、夏均未通過(guò)此檢驗(yàn),通過(guò)p<0.05檢驗(yàn)的站點(diǎn)數(shù)量為3(夏季),1(秋季),4(冬季),春季未通過(guò)此檢驗(yàn)。

    4.4 降水日數(shù)與降水強(qiáng)度的年代際變化空間分布

    20世紀(jì)60年代年均降水日數(shù)(圖9A)84.21%的站點(diǎn)為正增長(zhǎng)(0.188~1.546 d/a),15.79%的站點(diǎn)為負(fù)增長(zhǎng)(-0.418~-0.273 d/a),整體上表現(xiàn)出正增長(zhǎng)(均值為0.708 d/a);武威與門(mén)源成為增長(zhǎng)最快與下降最快的的站點(diǎn),二者相差1.127 d/a;核心研究區(qū)整體上為下降趨勢(shì)(均值為-0.197 d/a);其中武威與高臺(tái)通過(guò)了p<0.05的相關(guān)性檢驗(yàn)。20世紀(jì)60年代年均降水強(qiáng)度(圖9B)13個(gè)站點(diǎn)為正增長(zhǎng)(0.127~1.238 mm/10 a),6個(gè)站點(diǎn)為負(fù)增長(zhǎng)(-0.288~-0.234 mm/10 a),整體呈現(xiàn)正增長(zhǎng)(0.300 mm/10 a);永昌與民和分別為增長(zhǎng)最快與最慢的站點(diǎn),門(mén)源與德令哈為下降最快與最慢的站點(diǎn);核心研究區(qū)整體上呈上升趨勢(shì)(均值為0.171 mm/10 a);通過(guò)顯著性檢驗(yàn)的站點(diǎn)為祁連(p<0.01)、永昌(p<0.01)、都蘭(p<0.05)。20世紀(jì)70年代年均降水日數(shù)(圖9C)正負(fù)增長(zhǎng)之比為14∶5,整體上呈正增長(zhǎng)(均值為0.578 d/a);核心研究區(qū)整體上呈略微下降趨勢(shì)(均值為-0.001 5 d/a),外圍研究區(qū)呈增長(zhǎng)趨勢(shì)(均值為0.733 d/a)。20世紀(jì)70年代年均降水強(qiáng)度(圖9D)整體上表現(xiàn)為上升趨勢(shì)(均值為0.066 mm/10 a),核心研究區(qū)的傾斜率為0.326 mm/10 a,比全部站點(diǎn)的傾斜率大0.260 mm/10 a,外圍區(qū)的傾斜率為-0.003 1 mm/10 a,呈現(xiàn)略微下降的趨勢(shì);通過(guò)相關(guān)性檢驗(yàn)的比例為21.05%。

    圖8祁連山南坡及其附近地區(qū)降水日數(shù)與降水強(qiáng)度季節(jié)變化空間分布

    20世紀(jì)80年代年均降水日數(shù)(圖9E)變化范圍為-0.709~3.794 d/a,二者相差3.085 d/a,變化差異較大,整體上呈現(xiàn)增長(zhǎng)趨勢(shì)(均值為1.071 d/a);核心研究區(qū)的傾斜率為0.641 d/a,外圍研究區(qū)的傾斜率為1.186 d/a;正增長(zhǎng)的站點(diǎn)數(shù)量為13,負(fù)增長(zhǎng)的站點(diǎn)數(shù)量為6;剛察與托勒為增長(zhǎng)最快與最慢的站點(diǎn);剛察與野牛溝、茶卡、民和分別通過(guò)了p<0.01與p<0.05相關(guān)性檢驗(yàn)。20世紀(jì)80年代年均降水強(qiáng)度(圖9F)整體上呈增長(zhǎng)趨勢(shì)(0.249 mm/10 a),南坡傾斜率為0.304 mm/10 a,北坡傾斜率為0.130 mm/10 a,南坡較北坡增加的快;正增長(zhǎng)的比例為63.16%,負(fù)增長(zhǎng)的比例為36.84%;只有門(mén)源通過(guò)了相關(guān)性檢驗(yàn)(p<0.05)。20世紀(jì)90年代年均降水日數(shù)(圖9G)68.42%的站點(diǎn)為負(fù)增長(zhǎng),下降最快的為民和(-1.606 d/a),核心研究區(qū)整體上為負(fù)增長(zhǎng)(均值為-0.406 d/a),外圍研究區(qū)也為負(fù)增長(zhǎng)(均值為-0.253 d/a)。20世紀(jì)90年代年均降水強(qiáng)度(圖9H)正負(fù)增長(zhǎng)之比為12∶7,正增長(zhǎng)的范圍為0.015~1.383 mm/10 a,增長(zhǎng)最快與最慢的站點(diǎn)分別為剛察與門(mén)源,負(fù)增長(zhǎng)的范圍為-0.910~-0.081 mm/10 a,下降最快與最慢的站點(diǎn)分別為酒泉與都蘭;其中6個(gè)站點(diǎn)通過(guò)了相關(guān)性檢驗(yàn)(p<0.05)。2000—2014年年均降水日數(shù)(圖9I)的平均值為0.24 d/a,正增長(zhǎng)范圍為0.246~1.021 d/a,增長(zhǎng)最快的為茶卡,民勤為零增長(zhǎng),正增長(zhǎng)范圍為-0.496~-0.068 d/a,門(mén)源下降最快;茶卡與恰卜恰通過(guò)了p<0.05的相關(guān)性檢驗(yàn)。2000—2014年年均降水強(qiáng)度(圖9J)正負(fù)增長(zhǎng)之比為13∶6,平均值為0.133 mm/10 a,正增長(zhǎng)的范圍為0.037~0.812 mm/10 a,負(fù)增長(zhǎng)的范圍為-0.564~-0.050 mm/10 a;核心研究區(qū)的傾斜率為0.310 mm/10 a,外圍研究區(qū)的傾斜率為0.086 mm/10 a;南坡的傾斜率為0.202 mm/10 a,北坡的傾斜率為0.016 mm/10 a;其中野牛溝通過(guò)了p<0.01的相關(guān)性檢驗(yàn)。

    圖9祁連山南坡及其附近地區(qū)降水日數(shù)與降水強(qiáng)度年代際變化空間分布

    5 討 論

    (1) 時(shí)間序列變化。多年降水日數(shù)與降水強(qiáng)度整體上表現(xiàn)為波動(dòng)緩慢增長(zhǎng)趨勢(shì);祁連山區(qū)的降水主要受到太平洋暖濕氣流的影響,同時(shí)也受到來(lái)自大西洋冷濕氣流的影響,降水年際變化大,夏季降水強(qiáng)度對(duì)全年貢獻(xiàn)最大;20世紀(jì)80年代中后期之后較60—80年代中期的降水強(qiáng)度有明顯增加,這種趨勢(shì)與全國(guó)降水特征一致;季節(jié)突變減壓發(fā)現(xiàn)除了秋季突變不明顯,春、夏、冬發(fā)生突變的年份分別為1982年、1972年、1998年。

    (2) 空間分布規(guī)律。多年平均降水量與平均降水日數(shù)受到東南季風(fēng)的影響,均呈現(xiàn)出西北—東南走向的帶狀分布,各個(gè)站點(diǎn)降水分布呈現(xiàn)的走向與山脈走向一致,南坡為東南季風(fēng)的迎風(fēng)坡降水明顯多于處在背風(fēng)坡的北坡;影響降水的因素諸多,海拔、地形、迎風(fēng)坡與背風(fēng)坡、大氣環(huán)流等;1960—2014年各站點(diǎn)的年降水日數(shù)為正負(fù)增長(zhǎng)的站點(diǎn)數(shù)量之比大體相當(dāng),為10∶9,89.47%的站點(diǎn)年降水強(qiáng)度為正增長(zhǎng),這與張存杰等[29]研究結(jié)論基本一致;夏季降水強(qiáng)度增長(zhǎng)最明顯。

    (3) 20世紀(jì)80年代中后期降水強(qiáng)度有所增加,尤其進(jìn)入21世紀(jì)增加更為明顯,這會(huì)導(dǎo)致河流徑流量增加,比如黑河干流的徑流量從20世紀(jì)偏枯狀態(tài)到21世紀(jì)年徑流量出現(xiàn)回升的現(xiàn)象[30]。

    參考文獻(xiàn):

    [1] Wentz F J, Ricciardulli L, Hilburn K, et al. How much more rain will global warming bring[J]. Science, 2007,317(5835):233-235.

    [2] 秦大河, Thomas Stocker. IPCC第五次評(píng)估報(bào)告第一工作組報(bào)告的亮點(diǎn)結(jié)論[J].氣候變化研究進(jìn)展,2014,10(1):1-6.

    [3] 殷雪蓮,郭建華,董安祥,等.沿祁連山兩次典型強(qiáng)降水天氣個(gè)例對(duì)比分析[J].高原氣象,2008,27(1):184-192.

    [4] 劉雪梅,張明軍,王圣杰,等.2008—2014年祁連山區(qū)夏季降水的日變化特征及其影響因素[J].地理學(xué)報(bào),2016,71(5):754-767.

    [5] 陳乾,張久林,陳添宇.祁連山區(qū)2006年夏季降水過(guò)程的天氣分析[J].干旱氣象,2008,26(4):22-28.

    [6] 強(qiáng)芳,張明軍,王圣杰,等.基于格點(diǎn)數(shù)據(jù)的1961—2012年祁連山面雨量特征分析[J].地理學(xué)報(bào),2015,70(7):1125-1136.

    [7] 魏鋒,王勁松,李寶梓,等.祁連山近45 a5—9月日降水氣候特征[J].干旱氣象,2010,28(3):285-290.

    [8] 張小明,魏鋒,陸燕.祁連山近45 a年降水異常的氣候特征[J].干旱氣象,2006,24(3):35-41.

    [9] 張存杰,郭妮.祁連山區(qū)近40年氣候變化特征[J].氣象,2002,28(12):33-39.

    [10] 藍(lán)永超,丁永建,康爾泗.近50年來(lái)黑河山區(qū)匯流區(qū)溫度及降水變化趨勢(shì)[J].高原氣象,2004,23(5):723-727.

    [11] 牛赟,劉賢德,敬文茂,等.祁連山大野口流域氣溫、降水、河川徑流特征分析[J].干旱區(qū)地理:漢文版,2014,37(5):931-938.

    [12] 王寧練,賀建橋,蔣熹,等.祁連山中段北坡最大降水高度帶觀測(cè)與研究[J].冰川凍土,2009,31(3):395-403.

    [13] 賈文雄.近50年來(lái)祁連山及河西走廊降水的時(shí)空變化[J].地理學(xué)報(bào),2012,67(5):631-644.

    [14] 湯懋蒼.祁連山區(qū)降水的地理分布特征[J].地理學(xué)報(bào),1985(4):323-332.

    [15] 陳少勇,董安祥,韓通.祁連山東、西部夏季降水量時(shí)空分布的差異及其成因研究[J].南京氣象學(xué)院學(xué)報(bào),2007,30(5):715-719.

    [16] 李巖瑛,張強(qiáng),許霞,等.祁連山及周邊地區(qū)降水與地形的關(guān)系[J].冰川凍土,2010,32(1):52-61.

    [17] 賈文雄.祁連山氣候的空間差異與地理位置和地形的關(guān)系[J].干旱區(qū)研究,2010,27(4):607-615.

    [18] 程鵬,陶健紅,張新榮,等.高空西風(fēng)急流對(duì)祁連山區(qū)降水的影響[J].氣象科技,2007,35(4):489-494.

    [19] 常學(xué)向,趙愛(ài)芬,王金葉,等.祁連山林區(qū)大氣降水特征與森林對(duì)降水的截留作用[J].高原氣象,2002,21(3):274-280.

    [20] 張學(xué)龍,羅龍發(fā),敬文茂,等.祁連山青海云杉林截留對(duì)降水的分配效應(yīng)[J].山地學(xué)報(bào),2007,25(6):678-683.

    [21] 張強(qiáng),杜志成.絲綢之路經(jīng)濟(jì)帶上區(qū)域生態(tài)安全評(píng)價(jià)研究:以祁連山冰川與水涵養(yǎng)生態(tài)功能區(qū)為例[J].生態(tài)經(jīng)濟(jì):中文版,2016,32(10):169-173.

    [22] 王海軍,張勃,靳曉華,等.基于GIS的祁連山區(qū)氣溫和降水的時(shí)空變化分析[J].中國(guó)沙漠,2009,29(6):1196-1202.

    [23] 戴聲佩,張勃.基于GIS的祁連山植被NDVI對(duì)氣溫降水的旬響應(yīng)分析[J].生態(tài)環(huán)境學(xué)報(bào),2010,19(1):140-145.

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

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

    [26] 丁永建,葉柏生,周文娟.40年來(lái)西北干旱區(qū)黑河流域降水時(shí)空分布特征[J].冰川凍土,1999,21(1):42-48.

    [27] 賈文雄,何元慶,李宗省,等.祁連山區(qū)氣候變化的區(qū)域差異特征及突變分析[J].地理學(xué)報(bào),2008,63(3):257-269.

    [28] 秦大河,丁一匯,蘇紀(jì)蘭,等.中國(guó)氣候與環(huán)境演變(上卷)[M].北京:科學(xué)出版社,2005.

    [29] 張存杰,郭妮.祁連山區(qū)近40年氣候變化特征[J].氣象,2002, 28(12):33-39.

    [30] 藍(lán)永超,康爾泗,金會(huì)軍,等.黑河出山徑流量年際變化特征和趨勢(shì)研究[J].冰川凍土,1999,21(1):49-53.

    猜你喜歡
    南坡降水強(qiáng)度祁連山
    圖志
    發(fā)展(2024年1期)2024-04-25 03:45:30
    祁連山下
    南坡
    近59年江漢平原降水氣候變化特征分析
    不同降水強(qiáng)度下風(fēng)廓線雷達(dá)譜矩特征與測(cè)風(fēng)準(zhǔn)確性分析
    珠峰南坡傳統(tǒng)攀登路線
    2007—2017年大連市降水及降水日數(shù)變化特征分析
    祁連山南坡不同耕地類(lèi)型土壤理化特征性分析
    祁連山草原:如夢(mèng)如幻近高天
    祁連山
    黃河之聲(2016年20期)2016-02-21 11:55:33
    av免费在线看不卡| 亚洲丝袜综合中文字幕| 91aial.com中文字幕在线观看| 午夜免费男女啪啪视频观看| 国产精品国产三级国产专区5o | 日韩 亚洲 欧美在线| 亚洲av.av天堂| 午夜视频国产福利| 日本三级黄在线观看| 久久久久性生活片| 天美传媒精品一区二区| 日日啪夜夜撸| 网址你懂的国产日韩在线| 欧美激情国产日韩精品一区| 亚洲国产精品成人久久小说| 三级国产精品欧美在线观看| 91狼人影院| 亚洲在线自拍视频| 夜夜爽夜夜爽视频| 久久久a久久爽久久v久久| 国产色婷婷99| 99热精品在线国产| 中文字幕亚洲精品专区| 午夜激情欧美在线| 日本av手机在线免费观看| 国产精品国产三级国产av玫瑰| 久久精品91蜜桃| 久久久久久久久久久丰满| 午夜久久久久精精品| 久久久久性生活片| 1024手机看黄色片| 嫩草影院入口| 天天躁日日操中文字幕| 日本黄色视频三级网站网址| 欧美97在线视频| 91精品国产九色| 国产成人aa在线观看| 91aial.com中文字幕在线观看| 看非洲黑人一级黄片| 永久网站在线| 熟女电影av网| 特级一级黄色大片| 国产精华一区二区三区| 成人国产麻豆网| 国产免费一级a男人的天堂| 最近中文字幕2019免费版| 18+在线观看网站| a级毛色黄片| 麻豆av噜噜一区二区三区| 久久久久久久久久久丰满| 国产高清国产精品国产三级 | 欧美+日韩+精品| 欧美成人免费av一区二区三区| 免费观看a级毛片全部| 欧美成人免费av一区二区三区| 精品99又大又爽又粗少妇毛片| 精品人妻熟女av久视频| 亚洲人与动物交配视频| 99久久精品国产国产毛片| 搡老妇女老女人老熟妇| 黄片无遮挡物在线观看| 97人妻精品一区二区三区麻豆| 欧美一区二区精品小视频在线| 麻豆久久精品国产亚洲av| 九九在线视频观看精品| 国产女主播在线喷水免费视频网站 | 边亲边吃奶的免费视频| 波野结衣二区三区在线| 国产亚洲精品av在线| 女的被弄到高潮叫床怎么办| 国产精品蜜桃在线观看| 五月玫瑰六月丁香| 精品久久久久久久久av| 国产av不卡久久| 99久久中文字幕三级久久日本| 免费不卡的大黄色大毛片视频在线观看 | 亚洲不卡免费看| 波多野结衣高清无吗| 国产精品国产高清国产av| 国产精品人妻久久久久久| 久久久久久久久久久丰满| 日韩亚洲欧美综合| 一区二区三区免费毛片| 久久精品国产自在天天线| 亚洲国产精品sss在线观看| 三级男女做爰猛烈吃奶摸视频| 国产精品99久久久久久久久| 国产成人一区二区在线| 如何舔出高潮| 国产人妻一区二区三区在| 一二三四中文在线观看免费高清| 国产精品一及| 亚洲av成人av| 欧美三级亚洲精品| 网址你懂的国产日韩在线| 国产一区有黄有色的免费视频 | 69人妻影院| 噜噜噜噜噜久久久久久91| 夫妻性生交免费视频一级片| 免费观看在线日韩| 人妻夜夜爽99麻豆av| 大又大粗又爽又黄少妇毛片口| 亚洲欧美一区二区三区国产| 美女xxoo啪啪120秒动态图| 亚洲av一区综合| 国产精华一区二区三区| 成人鲁丝片一二三区免费| 天堂网av新在线| 亚洲国产成人一精品久久久| 亚洲不卡免费看| 亚洲图色成人| 免费在线观看成人毛片| 亚洲国产精品成人久久小说| 亚洲欧美中文字幕日韩二区| 欧美3d第一页| 免费黄网站久久成人精品| 亚洲国产精品国产精品| 青春草视频在线免费观看| 国产高清不卡午夜福利| 国产精品无大码| 国产一区二区在线av高清观看| 欧美成人a在线观看| 久久精品久久久久久久性| 大又大粗又爽又黄少妇毛片口| 欧美最新免费一区二区三区| 国产乱人视频| 国产精品,欧美在线| 丝袜美腿在线中文| 狂野欧美白嫩少妇大欣赏| 久久久a久久爽久久v久久| 热99re8久久精品国产| 级片在线观看| 色5月婷婷丁香| 国产精品综合久久久久久久免费| 精品久久国产蜜桃| 亚洲欧美一区二区三区国产| 亚洲国产日韩欧美精品在线观看| 国产精品综合久久久久久久免费| 九九热线精品视视频播放| 午夜福利视频1000在线观看| 亚洲在线自拍视频| 日本-黄色视频高清免费观看| av免费在线看不卡| 亚洲国产精品成人久久小说| 亚洲美女搞黄在线观看| 欧美日韩国产亚洲二区| 成人毛片60女人毛片免费| av播播在线观看一区| 草草在线视频免费看| 亚洲国产欧洲综合997久久,| 性色avwww在线观看| 99热这里只有是精品在线观看| 免费av观看视频| 国内少妇人妻偷人精品xxx网站| 亚洲一级一片aⅴ在线观看| 亚洲国产最新在线播放| 嘟嘟电影网在线观看| 波多野结衣巨乳人妻| 日本免费一区二区三区高清不卡| 精品久久国产蜜桃| 干丝袜人妻中文字幕| 熟妇人妻久久中文字幕3abv| 男的添女的下面高潮视频| 国产亚洲午夜精品一区二区久久 | 日韩成人av中文字幕在线观看| 久久精品国产亚洲av天美| 日本-黄色视频高清免费观看| 欧美性猛交黑人性爽| 18禁动态无遮挡网站| 色综合色国产| 九色成人免费人妻av| 桃色一区二区三区在线观看| 伊人久久精品亚洲午夜| 精品久久久久久电影网 | 伦理电影大哥的女人| 日本色播在线视频| 亚洲高清免费不卡视频| videos熟女内射| 亚洲中文字幕日韩| 亚洲18禁久久av| 国内精品宾馆在线| 国产伦精品一区二区三区四那| 久久精品人妻少妇| 男女国产视频网站| 亚洲av不卡在线观看| 国产老妇女一区| 国产精品三级大全| 乱码一卡2卡4卡精品| 天天躁夜夜躁狠狠久久av| 一区二区三区四区激情视频| 午夜精品一区二区三区免费看| 国产老妇女一区| 成年女人永久免费观看视频| 免费看av在线观看网站| 免费无遮挡裸体视频| 搡女人真爽免费视频火全软件| 91精品伊人久久大香线蕉| 免费黄网站久久成人精品| 白带黄色成豆腐渣| 最近最新中文字幕免费大全7| 国产在视频线精品| 真实男女啪啪啪动态图| 日本黄大片高清| 欧美精品国产亚洲| 建设人人有责人人尽责人人享有的 | 乱系列少妇在线播放| 亚洲精华国产精华液的使用体验| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 久久久国产成人精品二区| 日韩国内少妇激情av| 欧美bdsm另类| 卡戴珊不雅视频在线播放| 国产91av在线免费观看| 亚洲电影在线观看av| 国产亚洲91精品色在线| 男人狂女人下面高潮的视频| av天堂中文字幕网| 五月玫瑰六月丁香| 免费看av在线观看网站| 精品国产露脸久久av麻豆 | 天堂网av新在线| 精品人妻偷拍中文字幕| 综合色av麻豆| 欧美成人午夜免费资源| 两个人的视频大全免费| 久久久久久久久久黄片| 国产在视频线精品| 小蜜桃在线观看免费完整版高清| 黄片wwwwww| 亚洲五月天丁香| videos熟女内射| 中文字幕av成人在线电影| 免费观看人在逋| 亚洲精品456在线播放app| h日本视频在线播放| 色吧在线观看| 韩国av在线不卡| 秋霞伦理黄片| 久久亚洲精品不卡| 如何舔出高潮| 男人舔奶头视频| 午夜福利在线观看吧| 十八禁国产超污无遮挡网站| 久久久午夜欧美精品| 女人十人毛片免费观看3o分钟| 直男gayav资源| 精品熟女少妇av免费看| 深爱激情五月婷婷| 婷婷色av中文字幕| 国产免费福利视频在线观看| 激情 狠狠 欧美| 亚洲久久久久久中文字幕| 成年女人看的毛片在线观看| 最近视频中文字幕2019在线8| 久久久久性生活片| 在线观看66精品国产| 欧美精品国产亚洲| 麻豆乱淫一区二区| 国产午夜精品论理片| 亚洲一级一片aⅴ在线观看| 一级二级三级毛片免费看| 亚洲国产成人一精品久久久| 老司机影院毛片| 热99在线观看视频| 国产片特级美女逼逼视频| 少妇的逼水好多| 26uuu在线亚洲综合色| 大香蕉97超碰在线| 在线播放无遮挡| 免费播放大片免费观看视频在线观看 | 青青草视频在线视频观看| 国产精品久久电影中文字幕| 欧美又色又爽又黄视频| 女人十人毛片免费观看3o分钟| 真实男女啪啪啪动态图| 大话2 男鬼变身卡| 美女xxoo啪啪120秒动态图| 国产在线一区二区三区精 | 精品久久久久久久末码| 少妇熟女aⅴ在线视频| АⅤ资源中文在线天堂| 男女那种视频在线观看| 国内揄拍国产精品人妻在线| 一级黄片播放器| av在线蜜桃| 欧美成人a在线观看| 国产伦一二天堂av在线观看| 综合色av麻豆| 亚洲人成网站在线观看播放| 国产一区二区亚洲精品在线观看| 亚洲av成人精品一区久久| 国产一区二区在线观看日韩| 桃色一区二区三区在线观看| 欧美日韩综合久久久久久| 欧美一区二区国产精品久久精品| 久热久热在线精品观看| 国产精品.久久久| 亚洲av免费在线观看| 国产精品国产三级国产专区5o | a级毛色黄片| 亚洲成人av在线免费| 亚洲内射少妇av| 久久99热6这里只有精品| 熟女电影av网| 亚洲av电影在线观看一区二区三区 | 午夜福利成人在线免费观看| 国产亚洲5aaaaa淫片| 欧美色视频一区免费| 一级黄色大片毛片| 天天躁日日操中文字幕| 国产精品国产高清国产av| 我要搜黄色片| 一级毛片aaaaaa免费看小| 国产午夜精品论理片| 国产av不卡久久| 老师上课跳d突然被开到最大视频| 深夜a级毛片| 一级黄片播放器| 国产女主播在线喷水免费视频网站 | 亚洲欧美日韩卡通动漫| 少妇的逼水好多| 亚洲av成人av| 联通29元200g的流量卡| 久久人人爽人人爽人人片va| 国产成人精品一,二区| 久久亚洲精品不卡| 高清毛片免费看| 少妇人妻一区二区三区视频| 成年免费大片在线观看| 亚洲乱码一区二区免费版| 简卡轻食公司| 中文天堂在线官网| 欧美色视频一区免费| 中文亚洲av片在线观看爽| 一个人观看的视频www高清免费观看| 国产免费男女视频| 国产精品女同一区二区软件| 91精品一卡2卡3卡4卡| 国产成人福利小说| 一级毛片aaaaaa免费看小| 亚洲精品色激情综合| 天堂影院成人在线观看| 久久久成人免费电影| 久久久久久久国产电影| 国内揄拍国产精品人妻在线| 高清毛片免费看| 少妇丰满av| 国产精品久久久久久久久免| 欧美日韩精品成人综合77777| 秋霞在线观看毛片| 日韩一本色道免费dvd| 中文乱码字字幕精品一区二区三区 | 日韩欧美 国产精品| 禁无遮挡网站| 永久网站在线| 成人性生交大片免费视频hd| 成人毛片60女人毛片免费| 精品欧美国产一区二区三| 欧美日韩一区二区视频在线观看视频在线 | 成年女人看的毛片在线观看| 好男人在线观看高清免费视频| 亚洲国产高清在线一区二区三| 午夜老司机福利剧场| 三级经典国产精品| av国产久精品久网站免费入址| 免费av不卡在线播放| 亚洲国产精品成人综合色| 日韩国内少妇激情av| 日本av手机在线免费观看| 内射极品少妇av片p| 欧美一区二区亚洲| 久久精品国产亚洲av涩爱| 国产黄色小视频在线观看| 国产精品人妻久久久久久| 久久久久久大精品| 日日摸夜夜添夜夜添av毛片| 国产黄片视频在线免费观看| h日本视频在线播放| 亚洲高清免费不卡视频| 你懂的网址亚洲精品在线观看 | 汤姆久久久久久久影院中文字幕 | 别揉我奶头 嗯啊视频| 99久久精品热视频| 久久久久久久亚洲中文字幕| 一边亲一边摸免费视频| 国产精品电影一区二区三区| 中文字幕熟女人妻在线| 不卡视频在线观看欧美| 两个人视频免费观看高清| 看免费成人av毛片| 亚洲高清免费不卡视频| 天堂中文最新版在线下载 | 九色成人免费人妻av| 可以在线观看毛片的网站| 亚洲av成人av| 欧美高清成人免费视频www| 综合色av麻豆| 伊人久久精品亚洲午夜| 18禁在线无遮挡免费观看视频| 小说图片视频综合网站| 插阴视频在线观看视频| 国产精品一区二区性色av| 能在线免费观看的黄片| 国产三级在线视频| 国产91av在线免费观看| 最近最新中文字幕免费大全7| 亚洲在线自拍视频| 久久久久久久午夜电影| 国产伦理片在线播放av一区| 岛国在线免费视频观看| 欧美区成人在线视频| 春色校园在线视频观看| 成人毛片60女人毛片免费| 国产黄色视频一区二区在线观看 | 99久国产av精品| 日本色播在线视频| 国产av一区在线观看免费| 日韩三级伦理在线观看| 小蜜桃在线观看免费完整版高清| 在线观看美女被高潮喷水网站| 人妻制服诱惑在线中文字幕| 一级毛片我不卡| 寂寞人妻少妇视频99o| 精品人妻熟女av久视频| 国产精品无大码| 久久久久久久国产电影| 丰满人妻一区二区三区视频av| 男女那种视频在线观看| a级毛片免费高清观看在线播放| 日韩在线高清观看一区二区三区| 嫩草影院精品99| 国产精品乱码一区二三区的特点| 97热精品久久久久久| 国产av码专区亚洲av| 久久国内精品自在自线图片| 91久久精品电影网| 人妻系列 视频| 欧美日韩精品成人综合77777| 别揉我奶头 嗯啊视频| 爱豆传媒免费全集在线观看| 中文字幕av在线有码专区| av在线天堂中文字幕| 国产免费男女视频| 成人无遮挡网站| 国产免费福利视频在线观看| 日韩一区二区三区影片| 最新中文字幕久久久久| 十八禁国产超污无遮挡网站| 国产亚洲最大av| 色5月婷婷丁香| 国产免费视频播放在线视频 | 国产私拍福利视频在线观看| 真实男女啪啪啪动态图| 日韩强制内射视频| av福利片在线观看| 99久久精品国产国产毛片| 精品人妻偷拍中文字幕| 久热久热在线精品观看| 久久久久久伊人网av| 十八禁国产超污无遮挡网站| a级毛片免费高清观看在线播放| 欧美人与善性xxx| 亚洲精品影视一区二区三区av| 色网站视频免费| 国产精品嫩草影院av在线观看| 国产精品一区二区性色av| 麻豆乱淫一区二区| 搡女人真爽免费视频火全软件| 亚洲在线自拍视频| av在线天堂中文字幕| 看十八女毛片水多多多| 2021少妇久久久久久久久久久| 男人舔奶头视频| 中文字幕制服av| 在线免费观看的www视频| 国产欧美日韩精品一区二区| 青青草视频在线视频观看| 国产精品一区二区三区四区免费观看| 一区二区三区乱码不卡18| 永久网站在线| 美女高潮的动态| 精品久久久久久电影网 | 国产精品爽爽va在线观看网站| av在线亚洲专区| 亚洲精品影视一区二区三区av| 亚洲电影在线观看av| 少妇熟女欧美另类| 亚洲欧美成人综合另类久久久 | 亚洲av电影在线观看一区二区三区 | 如何舔出高潮| 亚洲18禁久久av| 久久久久免费精品人妻一区二区| 欧美xxxx黑人xx丫x性爽| 国产爱豆传媒在线观看| 欧美丝袜亚洲另类| 菩萨蛮人人尽说江南好唐韦庄 | 久久久久网色| 日本色播在线视频| 亚洲av成人精品一区久久| 国产高清国产精品国产三级 | 成人av在线播放网站| 亚洲精品日韩av片在线观看| 99热这里只有是精品50| 99九九线精品视频在线观看视频| 精品不卡国产一区二区三区| 国产淫片久久久久久久久| 日本与韩国留学比较| 级片在线观看| 舔av片在线| 在线观看66精品国产| 狂野欧美白嫩少妇大欣赏| 国产不卡一卡二| 免费观看精品视频网站| 久久精品久久精品一区二区三区| 国内精品美女久久久久久| 日本黄色片子视频| 丰满人妻一区二区三区视频av| 日韩在线高清观看一区二区三区| 成人无遮挡网站| 日韩精品青青久久久久久| 亚洲人成网站在线播| 国产成人精品婷婷| 少妇熟女aⅴ在线视频| 日韩亚洲欧美综合| 日本午夜av视频| 十八禁国产超污无遮挡网站| 一级毛片久久久久久久久女| 亚洲电影在线观看av| 欧美丝袜亚洲另类| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品嫩草影院av在线观看| 国产一区二区三区av在线| 亚洲人成网站高清观看| 我要搜黄色片| 亚洲中文字幕一区二区三区有码在线看| 欧美3d第一页| 五月玫瑰六月丁香| 日韩,欧美,国产一区二区三区 | 国产三级在线视频| 十八禁国产超污无遮挡网站| 在线观看美女被高潮喷水网站| 日本-黄色视频高清免费观看| 超碰97精品在线观看| av又黄又爽大尺度在线免费看 | 综合色av麻豆| 精品国产一区二区三区久久久樱花 | 国产又色又爽无遮挡免| 国模一区二区三区四区视频| 只有这里有精品99| videos熟女内射| 免费看美女性在线毛片视频| 日本-黄色视频高清免费观看| 免费看日本二区| 长腿黑丝高跟| 白带黄色成豆腐渣| 热99re8久久精品国产| 久久99热这里只有精品18| 日韩在线高清观看一区二区三区| 最近手机中文字幕大全| 2021天堂中文幕一二区在线观| 国产色爽女视频免费观看| 校园人妻丝袜中文字幕| 亚洲精品aⅴ在线观看| 亚洲国产成人一精品久久久| 欧美一区二区国产精品久久精品| 亚洲成人av在线免费| 色视频www国产| 美女黄网站色视频| 免费看光身美女| 成人性生交大片免费视频hd| 九九久久精品国产亚洲av麻豆| 日韩在线高清观看一区二区三区| 美女cb高潮喷水在线观看| 免费观看a级毛片全部| av专区在线播放| 一个人观看的视频www高清免费观看| 日韩精品青青久久久久久| videos熟女内射| 日本av手机在线免费观看| 91久久精品电影网| 亚洲精品一区蜜桃| 青青草视频在线视频观看| 我的老师免费观看完整版| 亚洲美女搞黄在线观看| 有码 亚洲区| 成人亚洲欧美一区二区av| 久久久久性生活片| 国产一区亚洲一区在线观看| 国产成人一区二区在线| 麻豆成人午夜福利视频| 精品久久久久久电影网 | 国语对白做爰xxxⅹ性视频网站| 国产真实伦视频高清在线观看| 国产精品麻豆人妻色哟哟久久 | 亚洲不卡免费看| av在线天堂中文字幕| 1000部很黄的大片| 色尼玛亚洲综合影院| 一级毛片久久久久久久久女| 欧美精品一区二区大全| 尤物成人国产欧美一区二区三区| 国产黄色小视频在线观看| 天天躁夜夜躁狠狠久久av| 中文在线观看免费www的网站| 69人妻影院| 亚洲av中文av极速乱| 国产成人a区在线观看| 少妇猛男粗大的猛烈进出视频 | 国产一级毛片七仙女欲春2| 蜜桃久久精品国产亚洲av| 美女xxoo啪啪120秒动态图| 日韩欧美 国产精品| 中文字幕人妻熟人妻熟丝袜美| 99九九线精品视频在线观看视频|