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

    1960-2011年河西走廊地表干濕狀況的時(shí)空變化及影響因素

    2015-03-14 01:00:46楊永春孫彥猛
    水土保持通報(bào) 2015年1期
    關(guān)鍵詞:河西走廊

    劉 揚(yáng), 楊永春,2, 張 軻, 孫彥猛

    (1.蘭州大學(xué) 資源環(huán)境學(xué)院, 甘肅 蘭州 730000; 2.蘭州大學(xué) 西部環(huán)境教育部

    重點(diǎn)實(shí)驗(yàn)室, 甘肅 蘭州 730000; 3.蘭州大學(xué) 信息科學(xué)與工程學(xué)院, 甘肅 蘭州 730000)

    1960-2011年河西走廊地表干濕狀況的時(shí)空變化及影響因素

    劉 揚(yáng)1, 楊永春1,2, 張 軻3, 孫彥猛3

    (1.蘭州大學(xué) 資源環(huán)境學(xué)院, 甘肅 蘭州 730000; 2.蘭州大學(xué) 西部環(huán)境教育部

    重點(diǎn)實(shí)驗(yàn)室, 甘肅 蘭州 730000; 3.蘭州大學(xué) 信息科學(xué)與工程學(xué)院, 甘肅 蘭州 730000)

    摘要:[目的] 為了探討河西走廊地區(qū)地表干濕狀況對(duì)氣候變化的響應(yīng)。[方法] 利用河西走廊地區(qū)14個(gè)氣象站1960—2011年逐日氣象數(shù)據(jù),基于FAO推薦的Penman—Monteith模型計(jì)算河西走廊各氣象站的參考蒸散量和濕潤(rùn)指數(shù),采用Sen’s斜率和Mann—Kendall趨勢(shì)檢驗(yàn)方法,分析河西地區(qū)地表干濕的時(shí)空變化特征,探討影響濕潤(rùn)指數(shù)變化的主導(dǎo)因子。[結(jié)果] 河西走廊地表多年平均干濕指數(shù)為0.00~0.56,且52 a來河西走廊濕潤(rùn)指數(shù)整體呈增加趨勢(shì);濕潤(rùn)指數(shù)的變化具有明顯地域性差異。各季節(jié)濕潤(rùn)指數(shù)亦整體呈增加趨勢(shì),以冬季增幅最大,秋、春次之,而夏季變化最小。[結(jié)論] 河西走廊地區(qū)氣候?qū)儆诟珊担?2 a來濕潤(rùn)指數(shù)略有增大,其主要原因是降水的增加和風(fēng)速的減小。

    關(guān)鍵詞:參考蒸散發(fā); 濕潤(rùn)指數(shù); Mann—Kendall趨勢(shì)檢驗(yàn); Sen’s斜率; 河西走廊

    在全球變暖背景下,環(huán)境狀況與區(qū)域差異是研究全球變化的區(qū)域響應(yīng)、土地沙漠化、土地覆被變化等環(huán)境變動(dòng)的基礎(chǔ)[1],氣候因子間的相互作用與制約產(chǎn)生了陸地表層的干濕狀況。濕潤(rùn)指數(shù)可客觀反映地表水熱收支狀況,是衡量區(qū)域干濕程度的重要指標(biāo),一般而言,區(qū)域的干濕狀況可由蒸發(fā)與降水之關(guān)系決定[2-3]。近年來,在中國(guó)西北地區(qū)、青藏高原等一些區(qū)域針對(duì)地表蒸散及干濕狀況及變化規(guī)律開展了廣泛研究[4-9],西北干旱半干旱區(qū)近50 a實(shí)測(cè)記錄顯示氣候環(huán)境出現(xiàn)了0.2 ℃/10 a的趨勢(shì)升溫[4,8,10-11],降水顯著增加,冰川消融加速,河川徑流增加,植被有所改善等現(xiàn)象[12]。河西走廊作為內(nèi)陸干旱化研究的氣候變化極度敏感區(qū)和生態(tài)脆弱區(qū),氣溫和降水增幅表現(xiàn)出明顯的空間差異[8,13],水資源受氣候變化影響而面臨諸多問題。本文利用河西走廊地區(qū)主要?dú)庀笳举Y料,結(jié)合彭曼模型計(jì)算區(qū)域尺度的參考蒸散量和濕潤(rùn)指數(shù),探討河西走廊地表干濕變化的程度及時(shí)空分布特征,為河川徑流變化研究、綠洲農(nóng)業(yè)結(jié)構(gòu)調(diào)整等提供科學(xué)依據(jù),尤其是為水資源的生態(tài)脆弱性評(píng)價(jià)提供科學(xué)的氣候變化背景,同時(shí)對(duì)揭示西北地區(qū)由暖干向暖濕轉(zhuǎn)型的氣候科學(xué)集成研究有重要的現(xiàn)實(shí)價(jià)值。

    1數(shù)據(jù)資料與模型

    1.1 資料選取

    選取河西走廊的敦煌、安西、玉門鎮(zhèn)、梧桐溝、酒泉、金塔、鼎新、高臺(tái)、張掖、山丹、永昌、民勤、武威和烏鞘嶺14個(gè)氣象站的1960—2011年(其中梧桐溝站的數(shù)據(jù)為1966—1988年;金塔站的數(shù)據(jù)為1989—2011年)逐日平均氣溫、最高氣溫、最低氣溫、平均相對(duì)濕度、平均風(fēng)速和日照時(shí)數(shù)等氣象觀測(cè)數(shù)據(jù),部分缺測(cè)數(shù)據(jù)由鄰近站點(diǎn)數(shù)據(jù)按差值訂正法進(jìn)行插補(bǔ),按照(FAO)Penman—Monteith模型在Matlab平臺(tái)編程計(jì)算各氣象站的逐日參考作物蒸散量,并基于降水?dāng)?shù)據(jù)計(jì)算月濕潤(rùn)指數(shù),進(jìn)而統(tǒng)計(jì)逐季(3—5月為春季,6—8月為夏季,9—11月為秋季,12—次年2月為冬季)和逐年地表濕潤(rùn)指數(shù)的時(shí)間序列。

    1.2 參考蒸散與濕潤(rùn)指數(shù)的計(jì)算

    本文首先采用Penman—Monteith模型[14],計(jì)算參考作物蒸散發(fā)ETo:

    (1)

    式中:ETo——參考蒸散量(mm/d);Δ——飽和水汽壓溫度曲線斜率(kPa/℃);G——地表以下的土壤熱通量〔MJ/(m2·d)〕;γ——干濕常數(shù)(kPa/℃);U2——2 m高處的風(fēng)速(m/s);es,ea——為飽和水汽壓和實(shí)際水汽壓(kPa);Rn——凈輻射〔MJ/(m2·d)〕,由凈短波輻射Rns與凈長(zhǎng)波輻射Rnl之差得出,即:

    Rn=Rns-Rnl

    (2)

    Yin等[15]利用中國(guó)81個(gè)氣象站的逐月輻射觀測(cè)數(shù)據(jù)建立了適用于中國(guó)凈短波輻射經(jīng)驗(yàn)公式,即:

    (3)

    (4)

    式中:α——地表反照率; FAO推薦取0.23[14],Tmax,Tmin——日最高氣溫和最低氣溫(K);σ——Stefan—Boltzman常數(shù)〔4.903×10-9MJ/(K4m2·d)〕;Rso——晴空太陽(yáng)輻射(MJ/d);n——日照時(shí)數(shù)(h);N——可照時(shí)數(shù)(h)。

    本文中計(jì)算ETo以日尺度計(jì)算,然后累加獲得月尺度的ETo。就濕潤(rùn)指數(shù)而言,定義月濕潤(rùn)指數(shù)W為月降水量P與月參考蒸散ETo之比,即:

    (5)

    式中:W——月地表濕潤(rùn)指數(shù);P——月降水量(mm); ETo——月參考蒸散量(mm)。

    1.3 Sen’s斜率及M-K趨勢(shì)檢驗(yàn)

    采用Sen’s斜率、10 a滑動(dòng)平均,并結(jié)合非參數(shù)Mann—Kendall統(tǒng)計(jì)(簡(jiǎn)稱M—K)趨勢(shì)檢驗(yàn)法[16]量化濕潤(rùn)指數(shù)的年、季變化率及變化程度,并通過ArcGIS的IDW插值得到濕潤(rùn)指數(shù)變化趨勢(shì)的空間分布,從而進(jìn)行地表干濕狀況時(shí)間序列的空間差異分析。

    Sen’s斜率是一種非參數(shù)趨勢(shì)斜率計(jì)算方法,不受序列奇異值的影響,能很好的反應(yīng)序列趨勢(shì)的變化程度,而滑動(dòng)平均可在一定程度上消除序列頻繁隨機(jī)起伏的影響,使水文變化的趨勢(shì)性或階段性更直觀、明顯,Sen’s斜率K的公式如下[17]:

    K=median(Qi)

    (6)

    (7)

    式中:i——常數(shù)取值為(1,2…,N), N——時(shí)間序列長(zhǎng)度; xj和xk則分別是j和k的序列之,且滿足j>k。

    M—K法[17]是用來評(píng)估水文氣候要素時(shí)間序列趨勢(shì)的非參數(shù)檢驗(yàn)方法,可以檢驗(yàn)時(shí)間序列變化趨勢(shì)的顯著性,以適用范圍廣,人為性少,定量化程度高而著稱,其檢驗(yàn)統(tǒng)計(jì)量公式為:

    (8)

    式中:sign()——符號(hào)函數(shù),Xi-Xj小于、等于或者大于0時(shí),sign(Xi-Xj)分別為-1,0和1。就S而言:

    (9)

    Z=0(S=0)

    (10)

    (11)

    Mann—Kendall趨勢(shì)檢驗(yàn)的統(tǒng)計(jì)量Z為正值表示增加趨勢(shì),負(fù)值表示減少趨勢(shì)。Z的絕對(duì)值在≥1.28,1.96,2.32時(shí),表明其顯著性分別通過了信度90%,95%和99%的顯著性檢驗(yàn)。

    2結(jié)果分析

    2.1 地表多年平均濕潤(rùn)指數(shù)的空間分布

    河西走廊地區(qū)多年平均濕潤(rùn)指數(shù)在0.00~0.56 (圖1),地表相對(duì)干旱,其中玉門、安西及敦煌的濕潤(rùn)指數(shù)在0.00~0.05,為河西地區(qū)最為干旱的區(qū)域。整體而言,地表干濕狀況總體表現(xiàn)出明顯的西北—東南差異,即濕潤(rùn)指數(shù)以敦煌為最低值,由此向東向南大致呈增加趨勢(shì),烏鞘嶺的濕潤(rùn)指數(shù)最高,西北部疏勒河流域(玉門、安西、敦煌)、中部黑河流域(山丹、張掖、高臺(tái)、酒泉、金塔、鼎新、梧桐溝)、東南部石羊河流域(武威、永昌、民勤、烏鞘嶺)的平均濕潤(rùn)指數(shù)分別為0.05,0.11及0.28。

    圖1 1960-2011年河西走廊地表平均濕潤(rùn)指數(shù)的空間分布

    2.2 地表干濕狀況的時(shí)空變化特征

    2.2.1年均濕潤(rùn)指數(shù)的空間變化整體平均而言,河西走廊年均濕潤(rùn)指數(shù)(14個(gè)氣象站均值)呈波狀起伏的增大趨勢(shì)(圖2),52a均值為0.136,最高為0.210 (2007年),最低為0.066 (1962年),變化幅度較大,20世紀(jì)60,70年代濕潤(rùn)指數(shù)距平為負(fù)距平,80年代后為正距平(表1),且自70年代以來有逐年增大趨勢(shì),尤其是1987年以來,大部分以正距平為主,表明在1987年之前地表相對(duì)干旱,之后有變濕趨勢(shì)。濕潤(rùn)指數(shù)年際變化的Sen’s斜率K為0.001 2,變化傾向率為0.009/10a,且氣候傾向率為增加趨勢(shì),表明在氣溫升高和降水量增加的背景下,河西走廊地表干旱狀況有明顯緩解。

    圖2 地表年均濕潤(rùn)指數(shù)的年際變化趨勢(shì)(1960-2011年)

    河西走廊地表濕潤(rùn)指數(shù)年際變化趨勢(shì)(圖3)及其顯著性(Mann—Kendall統(tǒng)計(jì)量Z,圖4)呈現(xiàn)一定的區(qū)域差異,濕潤(rùn)指數(shù)距平在不同年代的變化特征亦存在空間差異(圖5)。河西走廊整體表現(xiàn)為變濕趨勢(shì),年均濕潤(rùn)指數(shù)的Sen’s斜率K在0.000 0~0.005 5,由西北到東南,其變濕趨勢(shì)增大。具體而言,西北部疏勒河流域呈不明顯的變濕趨勢(shì),Sen’s斜率K在0.000 0~0.000 8;中部黑河流域表現(xiàn)為南北差異,自北向南變濕加強(qiáng),除酒泉以北地區(qū)變化不明顯外,其他大部分地區(qū)呈明顯變濕趨勢(shì),以張掖綠洲變濕最為顯著;東南部石羊河流域Sen’s斜率K在0.000 5~0.005 5,除民勤以北地區(qū)不明顯變濕,其他大部分地區(qū)顯著變濕,其中烏鞘嶺以南變濕最明顯。

    圖3 河西走廊地表年均濕潤(rùn)指數(shù)年際變化趨勢(shì)的空間分布

    圖4  河西走廊地表年均濕潤(rùn)指數(shù)年際變化Mann-Kendall趨勢(shì)檢驗(yàn)的空間分布(1960-2011年)

    表1 河西走廊地表年均濕潤(rùn)指數(shù)的年代際距平(1960-2011年)

    圖5地表年均濕潤(rùn)指數(shù)距平的空間分布

    2.2.2季均濕潤(rùn)指數(shù)的空間變化就季節(jié)而言(圖6),河西走廊春、夏、秋、冬4季多年平均濕潤(rùn)指數(shù)分別為0.087,0.224,0.182和0.066,表明河西走廊整體以夏季最濕,秋季次之,春冬最小。季均濕潤(rùn)指數(shù)年際變化的Sen’s斜率K分別為0.001 5,0.000 4,0.002 5和0.002 8,暗示4個(gè)季節(jié)濕潤(rùn)指數(shù)均呈上升趨勢(shì),其中冬季增大幅度最大,秋季和春季增大趨勢(shì)明顯,但夏季增大趨勢(shì)不太明顯。春、冬季均濕潤(rùn)指數(shù)年代距平變化與年均濕潤(rùn)指數(shù)的變化趨勢(shì)一致(表2),即濕潤(rùn)指數(shù)都在60,70年代為負(fù)距平,80年代以來為正距平,特別是2000年以后地表較大程度地變濕。夏季,濕潤(rùn)指數(shù)在80年代初之前呈波狀緩慢增大趨勢(shì),之后出現(xiàn)1983,1993年2個(gè)極大值點(diǎn)及1991,2001和2010年3個(gè)極小值點(diǎn),且正負(fù)距平交替出現(xiàn),表明這一時(shí)段內(nèi)濕潤(rùn)指數(shù)波動(dòng)起伏但整體平緩增加。秋季,濕潤(rùn)指數(shù)呈波狀起伏變化,在80年代中期之前呈波狀緩慢下降趨勢(shì),之后波狀上升趨勢(shì)明顯,正負(fù)距平波動(dòng)變化較大。

    表2 河西走廊地表季均濕潤(rùn)指數(shù)的年代際距平

    河西走廊各季節(jié)濕潤(rùn)指數(shù)的變化趨勢(shì)(Sen’s斜率K,圖7)及其顯著性存在空間差異。春季,地表干濕狀況由西北到東南變濕趨勢(shì)增加,西北部疏勒河流域的Sen’s斜率K在0.000 4~0.001 0,中部黑河流域除鼎新以北地區(qū)、張掖氣象站點(diǎn)周圍地區(qū)呈不明顯的變濕趨勢(shì)外,其他大部分地區(qū)也顯著變濕,Sen’s斜率K在0.000 1~0.001 6,東南部石羊河流域季均濕潤(rùn)指數(shù)年際變化的Sen’s斜率K在0.000 5~0.007 0,除民勤以北地區(qū)不明顯變濕,其他大部分地區(qū)顯著變濕,其中烏鞘嶺以南變濕最明顯;夏季,平原地區(qū)除永昌氣象站點(diǎn)周圍地區(qū)呈明顯的變濕趨勢(shì)外,絕大部分地區(qū)均呈不顯著的變濕趨勢(shì),Sen’s斜率K在0.000 0~0.000 6;秋季的變化趨勢(shì)空間分布與春季類似,西北部疏勒河流域得Sen’s斜率K在0.000 3~0.000 8,中部黑河流域除山丹氣象站點(diǎn)周圍地區(qū)呈明顯的變濕趨勢(shì)外,其他絕大部分地區(qū)均呈不明顯變濕趨勢(shì),Sen’s斜率K在0.000 3~0.003 0,東南部石羊河流域季均濕潤(rùn)指數(shù)年際變化的Sen’s斜率K在0.000 9~0.010 0,大部分地區(qū)呈不明顯變濕,但烏鞘嶺以南變濕最明顯;冬季,中部的黑河中游的變濕趨勢(shì)最明顯,Sen’s斜率K在0.000 5~0.004 5;除此之外,在西部的疏勒河流域也有明顯的變濕趨勢(shì)。

    圖6 1960-2011年地表季均濕潤(rùn)指數(shù)的年際變化趨勢(shì)

    圖7 1960-2011年河西走廊地表季均濕潤(rùn)指數(shù)年際變化趨勢(shì)(Sen’s斜率K)的空間分布

    2.3 濕潤(rùn)指數(shù)變化的影響因子及相關(guān)分析

    從濕潤(rùn)指數(shù)的定義公式來看,濕潤(rùn)指數(shù)的變化在很大程度上取決于降水量和潛在蒸散量的變化速率,而潛在蒸散量又受到各種氣象要素綜合作用的影響,如氣溫、風(fēng)速、日照時(shí)數(shù)、相對(duì)濕度等,為了探討河西走廊濕潤(rùn)指數(shù)增加的原因,除了選取直接相關(guān)的降水量以外,還選取了影響蒸發(fā)的6個(gè)氣象要素(風(fēng)速、平均氣溫、最高氣溫、最低氣溫、相對(duì)濕度和日照時(shí)數(shù)),根據(jù)公式(7)計(jì)算主要?dú)庀笳緷駶?rùn)指數(shù)與各氣象要素的Pearson相關(guān)系數(shù)(表3)。

    就年平均相關(guān)系數(shù)而言,近52 a濕潤(rùn)指數(shù)與平均風(fēng)速、平均氣溫、平均最高氣溫、平均最低氣溫、降水量、平均相對(duì)濕度、日照時(shí)數(shù)的相關(guān)系數(shù)分別為:-0.24,0.09,-0.07,0.22,0.98,0.46以及-0.21。各站與降水量、相對(duì)濕度呈現(xiàn)顯著的正相關(guān),西北部疏勒河流域、中部黑河流域和東南部石羊河流域氣象站年W均與相對(duì)濕度相關(guān)性很強(qiáng)(分別為0.44,0.45,0.47),表明對(duì)于干旱少雨的地區(qū)來說,濕潤(rùn)指數(shù)增減與相對(duì)濕度的變化息息相關(guān);濕潤(rùn)指數(shù)與降水量的空間分布規(guī)律非常相似,與降水量的相關(guān)系數(shù)接近于1(通過了0.01的置信度檢驗(yàn)),說明雖然河西走廊降水量少,蒸發(fā)強(qiáng)烈,但近52 a降水量在研究時(shí)段內(nèi)呈顯著增加趨勢(shì)(降水增率為3.94 mm/10 a),降水仍對(duì)地表濕度起主導(dǎo)作用。濕潤(rùn)指數(shù)與日照時(shí)數(shù)、平均風(fēng)速之間存在明顯的負(fù)相關(guān),表明日照時(shí)數(shù)和平均風(fēng)速的減少對(duì)濕潤(rùn)指數(shù)增加有重要作用,而與氣溫的相關(guān)程度較弱,說明了氣溫的增高或降低對(duì)濕潤(rùn)指數(shù)的影響較小。

    表3 河西走廊主要?dú)庀笳緷駶?rùn)指數(shù)的影響因子及相關(guān)系數(shù)

    注:*,**分別表示通過0.05,0.01的置信度檢驗(yàn);U為年平均風(fēng)速;Tmean為年平均氣溫;Tmax為年均最高氣溫;Tmin為年均最低氣溫;P為年降水量;Rh為年平均相對(duì)濕度;n為年日照時(shí)數(shù)。

    3討 論

    地表干濕的估算對(duì)于區(qū)域氣候變化研究、水資源評(píng)價(jià)與合理利用、農(nóng)業(yè)灌溉和水利工程建設(shè)均有重要意義。因此,全球及區(qū)域尺度的地表干濕特征、機(jī)制與模擬一直是生態(tài)、水文、氣象和農(nóng)業(yè)等領(lǐng)域共同關(guān)注的熱點(diǎn)問題,貫穿于氣候變化、水量平衡、旱澇極端事件、農(nóng)業(yè)灌溉與生態(tài)修復(fù)等各方面研究。然而,地表干濕又受到多種環(huán)境因素和生物因素的控制[2,18],過程甚是復(fù)雜多變,特別在較大空間尺度上,更是難以測(cè)量或預(yù)測(cè)。一般而言,區(qū)域尺度的干濕度多利用降水與潛在蒸散發(fā)之比進(jìn)行量度,而潛在蒸散發(fā)作為地球系統(tǒng)水循環(huán)的重要環(huán)節(jié),是熱量與水量平衡的關(guān)鍵參數(shù),在生物圈—水圈—大氣圈中發(fā)揮著不可或缺的作用,與降水共同決定著區(qū)域的干濕狀況,并直接影響著全球與區(qū)域尺度的水量平衡和氣候變化[1]。事實(shí)上,潛在蒸發(fā)的估算方法眾多,例如常用的Priestley—Taylor法[19]等,而本文利用的P—M方法系聯(lián)合國(guó)糧農(nóng)組織推薦的方法[14],具有較好的可信度[18]。在全球變暖的背景下,中國(guó)變濕地區(qū)面積大于變干地區(qū)面積,但是干旱化進(jìn)程大于濕潤(rùn)化進(jìn)程。近十幾年來西北大部分地區(qū)氣候環(huán)境悄然發(fā)生了重大變化,出現(xiàn)了氣溫持續(xù)上升、降水顯著增加、冰川消融加速、河川徑流增加、湖泊水位上升、面積擴(kuò)大、洪澇災(zāi)害頻發(fā)、植被有所改善等現(xiàn)象[12],而河西地區(qū)又是甘肅省農(nóng)業(yè)發(fā)展的核心區(qū)域,其氣候變化對(duì)地表干濕的影響研究更為重要。2003年施雅風(fēng)等[12]提出了西北地區(qū)從“暖干”即將轉(zhuǎn)型為“暖濕”的觀點(diǎn),引起了學(xué)術(shù)界的廣泛關(guān)注?;诖?,作者利用降水與蒸散發(fā)之比,推出河西地區(qū)地表整體上有變濕趨勢(shì),這與前人的研究結(jié)果較為一致[6,9]。而馬柱國(guó)等[20]利用Palmer干旱指數(shù)推出西北西部地表為變濕趨勢(shì)、西北地區(qū)東部表現(xiàn)為變干趨勢(shì)的結(jié)論;這可能是由于其計(jì)算方法與本文的方法略有差異所導(dǎo)致。與此同時(shí),本文基于相關(guān)分析表明,河西地區(qū)的濕潤(rùn)指數(shù)增大主要是由于降水的增大和風(fēng)速的減小所致。曹雯等[21]通過敏感性分析表明,風(fēng)速對(duì)西北地區(qū)蒸散發(fā)的變化貢獻(xiàn)最大,河西走廊[9]、天山[10]等地的研究結(jié)果亦支持這一結(jié)論,即風(fēng)速的減小是導(dǎo)致蒸散發(fā)減小的主要原因,本文研究發(fā)現(xiàn),研究區(qū)風(fēng)速主要呈現(xiàn)減小趨勢(shì),故導(dǎo)致了河西地區(qū)的蒸散發(fā)減小,這就使得區(qū)域內(nèi)的水分支出部分減小,同時(shí)降水的略微增大,綜合使得地表濕潤(rùn)指數(shù)增大。

    4結(jié) 論

    (1) 河西走廊地表多年平均干濕指數(shù)為0.00~0.56,皆較為干旱,區(qū)內(nèi)呈現(xiàn)從西北—東南逐漸增大的趨勢(shì),即疏勒河流域濕潤(rùn)指數(shù)最小,石羊河流域最大。就季節(jié)而言,河西走廊濕潤(rùn)指數(shù)以夏季最大,秋季次之,春冬最小。

    (2) 整體而言,1960—2011年,河西走廊年均濕潤(rùn)指數(shù)以0.009/10 a的速率呈現(xiàn)出增大趨勢(shì)。具體而言,濕潤(rùn)指數(shù)的增大主要發(fā)生在走廊東部的石羊河流域和中部的黑河流域的南部,而走廊西部的疏勒河流域之濕潤(rùn)指數(shù)沒有顯著變化。

    (3) 1960—2011年,各季節(jié)濕潤(rùn)指數(shù)亦皆表現(xiàn)出增大趨勢(shì),其中以冬季增大幅度最大,而夏季的增大幅度最小。具體而言,春季、秋季濕潤(rùn)指數(shù)變化速率的空間分布與年濕潤(rùn)指數(shù)變化的空間分布相似。而夏季除永昌外,其余地區(qū)變濕趨勢(shì)不甚明顯;冬季則以黑河中游的變濕趨勢(shì)最大,并且疏勒河流域亦有明顯變濕。

    (4) 盡管研究時(shí)段內(nèi)氣溫明顯升高,然降水的增多和風(fēng)速的減小是導(dǎo)致研究區(qū)內(nèi)濕潤(rùn)指數(shù)增大的主要原因。

    [參考文獻(xiàn)]

    [1]吳紹洪,尹云鶴,鄭度,等.近30 a中國(guó)陸地表層干濕狀況研究[J].中國(guó)科學(xué):D輯,2005,35(3):276-283.

    [2]Ma Ning, Wang Naiang, Zhao Liqiang, et al. Observation of mega-dune evaporation after various rain events in the hinterland of Badain Jaran desert [J]. Chinese Science Bulletin, 2014, 59(2):162-170.

    [3]郭燕紅,張寅生,馬穎釗,等. 藏北羌塘高原雙湖地表熱源強(qiáng)度及地表水熱平衡[J].地理學(xué)報(bào),2014,69(7):983-992.

    [4]馬寧,王乃昂,朱金峰,等.巴丹吉林沙漠周邊地區(qū)近50 a來氣候變化[J].中國(guó)沙漠,2011,31(6):1541-1547.

    [5]馬寧,王乃昂,王鵬龍,等.黑河流域參考蒸散的時(shí)空變化特征及影響因素的定量分析[J].自然資源學(xué)報(bào),2012,27(6):975-989.

    [6]黃小燕,張明軍.中國(guó)西北地區(qū)地表干濕變化及影響因素[J].水科學(xué)進(jìn)展,2011,22(2):151-159.

    [7]王乃昂,馬寧,陳紅寶,等.巴丹吉林沙漠腹地降水特征初步分析[J].水科學(xué)進(jìn)展,2013,24(2):153-160.

    [8]馬寧,王乃昂,李卓侖,等. 1960—2009年巴丹吉林沙漠南北緣氣候變化分析[J]. 干旱區(qū)研究,2011,28(2):242-250.

    [9]賈文雄,何元慶,李宗省,等.祁連山及河西走廊地表干濕變化的時(shí)空分布[J].地球科學(xué):中國(guó)地質(zhì)大學(xué)學(xué)報(bào),2010,35(2):268-276.

    [10]張明軍,李瑞雪,賈文雄,等.中國(guó)天山山區(qū)潛在蒸發(fā)量的時(shí)空變化[J].地理學(xué)報(bào),2009,64(7):798-806.

    [11]秦大河,丁一匯,王紹武,等.中國(guó)西部生態(tài)環(huán)境變化與對(duì)策建議[J].地球科學(xué)進(jìn)展,2002,17(3):314-319.

    [12]施雅風(fēng),沈永平,李棟梁,等.中國(guó)西北氣候由暖干向暖濕轉(zhuǎn)型的特征和趨勢(shì)探討[J].第四紀(jì)研究,2003,23(2):152-164.

    [13]賈文雄,何元慶,李宗省,等.祁連山及河西走廊氣候變化的時(shí)空分布特征[J].中國(guó)沙漠,2008,28(6):1151-1155.

    [14]Allen R G, Pereira L S, Raes D, et al. Crop evapotranspiration: Guidelines for computing crop water requirements [M]. Rome: United Nations Food and Agriculture Organization, 1998.

    [15]Yin Yunhe, Wu Shaohong, Zheng Du, et al, Radiation calibration of FAO56 Penman—Monteith model to estimate reference crop evapotranspiration in China [J]. Agricultural Water Management, 2008,95(1):77-84.

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

    [17]易湘生,尹衍雨,李國(guó)勝,等.青海三江源地區(qū)近52 a來的氣溫變化[J].地理學(xué)報(bào),2011,66(11):1451-1465.

    [18]劉紹民,孫中平,李小文,等.蒸散量測(cè)定與估算方法的對(duì)比研究[J].自然資源學(xué)報(bào),2003,18(2):161-167.

    [19]Priestley C H B, Taylor R J. On the assessment of the surface heat flux and evaporation using large-scale parameters[J]. Monthly Weather Review, 1979,100(2),81-92.

    [20]馬柱國(guó),符淙斌.1951—2004年中國(guó)北方干旱化的基本事實(shí)[J].科學(xué)通報(bào),2006,51(20):2429-2439.

    [21]曹雯,申雙和,段春峰.西北地區(qū)生長(zhǎng)季參考作物蒸散變化成因的定量分析[J].地理學(xué)報(bào),2011,66(3):407-415.

    Temporal and Spatial Variation of Humidity and Its Influential Factors in Hexi Corridor During 1960-2011

    LIU Yang1, YANG Yongchun1,2, ZHANG Ke3, SUN Yanmeng3

    (1.CollegeofEarthandEnvironmentalScience,LanzhouUniversity,Lanzhou,Gansu730000,China;

    2.KeyLaboratoryofWesternChina’sEnvironmentalSystems,LanzhouUniversity,Lanzhou,Gansu

    730000,China; 3.CollegeofInformationScience&Engineering,LanzhouUniversity,Lanzhou,Gansu730000,China)

    Abstract:[Objective] To explore the response of land surface humidity to the regional climate change in Hexi Corridor.[Methods] The daily meteorological data from 14 stations during the period from 1960 to 2011 were used to calculate the referential evapotranspiration and humid index for each site based on the Penman—Monteith model. The Sen’s slope and Mann—Kendall trend test were used to analyze the temporal and spatial variation of humidity and its influential factors in the study area.[Results] The humid index of Hexi Corridor was around 0.00~0.56, and the humid index increased during the past 52 years. The humidity conditions showed a spatial variation due to the differences in climate and terrain. The seasonal humid index showed an increasing trend as a whole, with the greatest increases in winter, while smallest changes in summer. [Conclusion] Hexi Corridor is subject to arid climate, while the humid index increased slightly in the past 52 years due to the increase in precipitation and decrease in wind speed.

    Keywords:reference evapotranspiration; humid index; Mann-Kendall trend test; Sen’s slope; Hexi Corridor

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1000-288X(2015)01-0054-07

    中圖分類號(hào):P426.2

    通信作者:楊永春(1969—),男(漢族),陜西省白水縣人,博士,教授,博士生導(dǎo)師,主要從事城市與區(qū)域發(fā)展研究。E-mail:yangych@lzu.edu.cn。

    收稿日期:2013-12-06修回日期:2014-01-10

    資助項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目“制度非均衡供給與后進(jìn)轉(zhuǎn)型:漸進(jìn)制度改革背景下的中國(guó)西部城市生長(zhǎng)轉(zhuǎn)型的格局、模式與動(dòng)力機(jī)制”(41171143)

    第一作者:劉揚(yáng)(1987—),女(漢族),河北省衡水市人,博士研究生,研究方向?yàn)槌鞘协h(huán)境對(duì)氣候變化的響應(yīng)。E-mail:tifiyang@gmail.com。

    猜你喜歡
    河西走廊
    《河西走廊通史》編纂工作推進(jìn)會(huì)在蘭州召開
    絲綢之路(2024年1期)2024-04-26 09:38:17
    河西走廊—塔克拉瑪干沙漠邊緣阻擊戰(zhàn)正式啟動(dòng)
    河西走廊的風(fēng)
    在河西走廊聆聽
    文苑(2020年8期)2020-09-09 09:30:44
    河西走廊寶卷文化傳統(tǒng)的經(jīng)驗(yàn)價(jià)值與現(xiàn)實(shí)意義
    民國(guó)時(shí)期哈薩克族在河西走廊的活動(dòng)述論
    河西走廊沙漠開發(fā)區(qū)飼用甜高粱高效栽培技術(shù)
    河西走廊沙塵暴日數(shù)變化及其與氣候因素的關(guān)系
    河西走廊戈壁主要類型及其分布
    河西走廊5種沙棘果主要活性成分及對(duì)ROS清除作用比較分析
    国产熟女欧美一区二区| 精品人妻1区二区| 操美女的视频在线观看| 99久久精品国产亚洲精品| 免费久久久久久久精品成人欧美视频| 欧美另类一区| 丁香六月欧美| 午夜91福利影院| 久久久精品区二区三区| 1024香蕉在线观看| 一级,二级,三级黄色视频| 最新在线观看一区二区三区 | 少妇 在线观看| 丁香六月天网| 少妇裸体淫交视频免费看高清 | 天天躁狠狠躁夜夜躁狠狠躁| 国产精品一二三区在线看| 狠狠精品人妻久久久久久综合| 欧美精品高潮呻吟av久久| 亚洲中文av在线| 亚洲av欧美aⅴ国产| 成人手机av| 又大又黄又爽视频免费| 国产精品人妻久久久影院| 伦理电影免费视频| 亚洲国产最新在线播放| 亚洲国产欧美日韩在线播放| 欧美变态另类bdsm刘玥| 在线看a的网站| 日本午夜av视频| 亚洲综合色网址| 啦啦啦在线免费观看视频4| 亚洲一区中文字幕在线| tube8黄色片| 一二三四社区在线视频社区8| 精品国产一区二区三区久久久樱花| 久久久久久久久久久久大奶| 亚洲人成电影免费在线| 亚洲精品国产色婷婷电影| 啦啦啦中文免费视频观看日本| a级毛片黄视频| 老司机靠b影院| 欧美成人精品欧美一级黄| 波多野结衣av一区二区av| 国产伦理片在线播放av一区| 精品国产乱码久久久久久男人| 日本色播在线视频| 久久久久精品国产欧美久久久 | 午夜av观看不卡| 一本色道久久久久久精品综合| 中国国产av一级| 午夜福利影视在线免费观看| 国产亚洲av高清不卡| 狠狠精品人妻久久久久久综合| 一区二区日韩欧美中文字幕| 女人久久www免费人成看片| 久久99一区二区三区| 欧美日韩视频精品一区| 大片电影免费在线观看免费| 亚洲精品av麻豆狂野| 曰老女人黄片| videos熟女内射| www.999成人在线观看| 亚洲国产精品一区二区三区在线| 99香蕉大伊视频| 高清av免费在线| 啦啦啦视频在线资源免费观看| 超碰成人久久| 国产高清视频在线播放一区 | 深夜精品福利| 中文字幕人妻丝袜制服| 免费观看a级毛片全部| 欧美日韩亚洲国产一区二区在线观看 | 国产一区二区 视频在线| 美国免费a级毛片| 亚洲精品国产区一区二| 欧美黑人精品巨大| 久久精品国产亚洲av高清一级| 亚洲七黄色美女视频| 日韩伦理黄色片| 亚洲精品国产色婷婷电影| 2021少妇久久久久久久久久久| 免费看av在线观看网站| 亚洲精品一二三| 黄频高清免费视频| 国产麻豆69| 黄色a级毛片大全视频| 亚洲午夜精品一区,二区,三区| 99国产精品99久久久久| av线在线观看网站| 国产精品九九99| 亚洲一码二码三码区别大吗| 国产视频首页在线观看| 老司机在亚洲福利影院| 男女无遮挡免费网站观看| 亚洲精品美女久久av网站| 天天躁夜夜躁狠狠久久av| 国产一区二区在线观看av| 亚洲国产最新在线播放| 又紧又爽又黄一区二区| 午夜激情av网站| 真人做人爱边吃奶动态| av天堂久久9| 国产成人欧美在线观看 | 赤兔流量卡办理| 欧美在线黄色| 精品人妻一区二区三区麻豆| 黄色 视频免费看| 老司机靠b影院| 天天躁夜夜躁狠狠躁躁| 丝瓜视频免费看黄片| 国产深夜福利视频在线观看| 一二三四社区在线视频社区8| 纵有疾风起免费观看全集完整版| 国产成人91sexporn| 高潮久久久久久久久久久不卡| 亚洲国产精品成人久久小说| 青春草视频在线免费观看| 成人亚洲精品一区在线观看| 一区在线观看完整版| 丰满迷人的少妇在线观看| 欧美精品亚洲一区二区| 人妻 亚洲 视频| 最近最新中文字幕大全免费视频 | 国产亚洲一区二区精品| 久久毛片免费看一区二区三区| 国产亚洲精品久久久久5区| 国产极品粉嫩免费观看在线| 国产精品一区二区在线不卡| 一级a爱视频在线免费观看| 亚洲欧美色中文字幕在线| 90打野战视频偷拍视频| 亚洲国产av影院在线观看| 精品少妇久久久久久888优播| 久热爱精品视频在线9| tube8黄色片| 黑丝袜美女国产一区| 99热国产这里只有精品6| 在线观看国产h片| 国产老妇伦熟女老妇高清| 自线自在国产av| 国产不卡av网站在线观看| 久久久久久久久久久久大奶| av天堂久久9| 丝袜脚勾引网站| 国产免费又黄又爽又色| 久久人妻福利社区极品人妻图片 | 亚洲国产精品国产精品| 国产三级黄色录像| 91国产中文字幕| 久久久精品94久久精品| 狂野欧美激情性bbbbbb| 一区福利在线观看| 国产亚洲精品久久久久5区| 日本vs欧美在线观看视频| 黄色毛片三级朝国网站| 999精品在线视频| 久久av网站| 日韩av在线免费看完整版不卡| 多毛熟女@视频| 手机成人av网站| 久久久国产欧美日韩av| 你懂的网址亚洲精品在线观看| 国产成人91sexporn| 久久国产亚洲av麻豆专区| 另类精品久久| 黄色一级大片看看| 性少妇av在线| 一本大道久久a久久精品| 男人舔女人的私密视频| 男女床上黄色一级片免费看| 国产精品成人在线| 国产成人av激情在线播放| 欧美日韩黄片免| 国产又色又爽无遮挡免| 又大又爽又粗| 国产成人av激情在线播放| 久热爱精品视频在线9| 精品一区二区三卡| 亚洲精品自拍成人| 国产成人啪精品午夜网站| 久久久久网色| 搡老岳熟女国产| 亚洲人成77777在线视频| 最黄视频免费看| 久久人人爽av亚洲精品天堂| 在线亚洲精品国产二区图片欧美| 久久久国产精品麻豆| 男人操女人黄网站| 国产人伦9x9x在线观看| 老汉色∧v一级毛片| 捣出白浆h1v1| 18禁黄网站禁片午夜丰满| av国产精品久久久久影院| 久久久久久人人人人人| 90打野战视频偷拍视频| 亚洲免费av在线视频| 亚洲色图综合在线观看| 国产精品香港三级国产av潘金莲 | 18禁裸乳无遮挡动漫免费视频| 欧美日韩国产mv在线观看视频| av线在线观看网站| 天天操日日干夜夜撸| 国产亚洲午夜精品一区二区久久| 精品国产超薄肉色丝袜足j| 成人18禁高潮啪啪吃奶动态图| 我的亚洲天堂| 国产在线一区二区三区精| 波野结衣二区三区在线| 日韩中文字幕视频在线看片| 777久久人妻少妇嫩草av网站| 久久精品国产亚洲av涩爱| 欧美人与性动交α欧美精品济南到| 啦啦啦中文免费视频观看日本| 嫩草影视91久久| 夫妻性生交免费视频一级片| 18禁裸乳无遮挡动漫免费视频| 国产1区2区3区精品| 亚洲国产欧美网| 国产片特级美女逼逼视频| 一二三四在线观看免费中文在| 91字幕亚洲| 久久毛片免费看一区二区三区| 成人亚洲欧美一区二区av| 在线观看免费日韩欧美大片| 女警被强在线播放| 高清av免费在线| 9色porny在线观看| 久久精品国产综合久久久| 日韩 亚洲 欧美在线| 999久久久国产精品视频| 国产又爽黄色视频| 日本黄色日本黄色录像| 久久99一区二区三区| 精品国产乱码久久久久久小说| 啦啦啦中文免费视频观看日本| 丝瓜视频免费看黄片| 中文字幕另类日韩欧美亚洲嫩草| 美女扒开内裤让男人捅视频| 999久久久国产精品视频| 美女脱内裤让男人舔精品视频| 爱豆传媒免费全集在线观看| 少妇 在线观看| 欧美激情 高清一区二区三区| 69精品国产乱码久久久| 午夜日韩欧美国产| 一个人免费看片子| 人成视频在线观看免费观看| 伦理电影免费视频| 老汉色∧v一级毛片| 午夜福利在线免费观看网站| 亚洲精品乱久久久久久| 精品福利观看| 国产又色又爽无遮挡免| 后天国语完整版免费观看| 日日夜夜操网爽| www.自偷自拍.com| 免费观看人在逋| svipshipincom国产片| 亚洲熟女毛片儿| 19禁男女啪啪无遮挡网站| 日本wwww免费看| 国产精品久久久久成人av| 在线观看免费视频网站a站| 大香蕉久久成人网| 热re99久久国产66热| 亚洲第一av免费看| 成人黄色视频免费在线看| 久久国产亚洲av麻豆专区| videos熟女内射| 天天操日日干夜夜撸| 国产一区二区 视频在线| 黑人巨大精品欧美一区二区蜜桃| 成年美女黄网站色视频大全免费| 欧美日韩视频高清一区二区三区二| 在线亚洲精品国产二区图片欧美| 人成视频在线观看免费观看| 国产成人av激情在线播放| 日韩电影二区| 制服诱惑二区| 热99国产精品久久久久久7| 久久亚洲国产成人精品v| 黄色片一级片一级黄色片| 老司机深夜福利视频在线观看 | 亚洲,欧美,日韩| 久久精品久久久久久噜噜老黄| 免费日韩欧美在线观看| 国产不卡av网站在线观看| 亚洲欧洲国产日韩| 美女国产高潮福利片在线看| 黄片小视频在线播放| 亚洲人成电影免费在线| 老司机亚洲免费影院| 国产91精品成人一区二区三区 | 咕卡用的链子| 搡老岳熟女国产| 久久精品亚洲av国产电影网| 男男h啪啪无遮挡| 在线观看一区二区三区激情| 超色免费av| 免费一级毛片在线播放高清视频 | 亚洲av电影在线进入| 各种免费的搞黄视频| 97人妻天天添夜夜摸| 欧美黄色片欧美黄色片| 日韩制服骚丝袜av| 女性被躁到高潮视频| 亚洲成人国产一区在线观看 | 一边亲一边摸免费视频| 亚洲熟女毛片儿| 天天添夜夜摸| 亚洲国产精品一区三区| av在线老鸭窝| 婷婷色综合大香蕉| 天天躁日日躁夜夜躁夜夜| 99久久人妻综合| 欧美日本中文国产一区发布| 久久久国产一区二区| 久久人人爽人人片av| 亚洲av日韩精品久久久久久密 | 看免费av毛片| 久久久久久久国产电影| 国产有黄有色有爽视频| 国产成人啪精品午夜网站| 国产亚洲一区二区精品| 久久中文字幕一级| 亚洲情色 制服丝袜| 丝瓜视频免费看黄片| 色精品久久人妻99蜜桃| 精品久久久精品久久久| 男女床上黄色一级片免费看| 777米奇影视久久| 国产精品免费视频内射| 久久中文字幕一级| 人人妻人人添人人爽欧美一区卜| tube8黄色片| 麻豆av在线久日| 9191精品国产免费久久| 国产欧美日韩精品亚洲av| 午夜激情av网站| 久久女婷五月综合色啪小说| 久久久久视频综合| 悠悠久久av| 国产有黄有色有爽视频| 午夜福利影视在线免费观看| 在线天堂中文资源库| 最新在线观看一区二区三区 | 老汉色av国产亚洲站长工具| 18禁裸乳无遮挡动漫免费视频| 激情视频va一区二区三区| 黑丝袜美女国产一区| 免费不卡黄色视频| 精品国产国语对白av| 自线自在国产av| 亚洲专区中文字幕在线| 少妇人妻 视频| 黑人欧美特级aaaaaa片| 亚洲成人免费电影在线观看 | 性色av乱码一区二区三区2| 少妇人妻 视频| 精品卡一卡二卡四卡免费| 性色av一级| 精品少妇黑人巨大在线播放| 观看av在线不卡| 国产视频首页在线观看| 一区二区三区激情视频| 另类亚洲欧美激情| 日韩一本色道免费dvd| 99热网站在线观看| 久久国产亚洲av麻豆专区| 国产爽快片一区二区三区| 精品免费久久久久久久清纯 | av片东京热男人的天堂| 97精品久久久久久久久久精品| 一级毛片黄色毛片免费观看视频| 天堂8中文在线网| 国产亚洲av片在线观看秒播厂| 最新的欧美精品一区二区| 99国产精品99久久久久| 国产又色又爽无遮挡免| 免费观看a级毛片全部| 亚洲精品日本国产第一区| 成人免费观看视频高清| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久人人人人人| 久久久国产精品麻豆| 国产亚洲欧美精品永久| 亚洲伊人色综图| 免费高清在线观看视频在线观看| 曰老女人黄片| 欧美日韩综合久久久久久| 亚洲国产精品一区三区| 人人妻,人人澡人人爽秒播 | 9色porny在线观看| 午夜福利,免费看| 精品欧美一区二区三区在线| 在线观看免费视频网站a站| 在线观看免费午夜福利视频| 国产日韩欧美视频二区| 成人手机av| 亚洲av五月六月丁香网| 中文字幕另类日韩欧美亚洲嫩草| 啦啦啦韩国在线观看视频| 黄色视频不卡| 50天的宝宝边吃奶边哭怎么回事| 欧美+亚洲+日韩+国产| 午夜激情福利司机影院| 精品一区二区三区视频在线观看免费| 色播亚洲综合网| 亚洲色图av天堂| 久久中文字幕人妻熟女| 久久久久久久久中文| 国产免费av片在线观看野外av| 色尼玛亚洲综合影院| 午夜福利在线在线| 香蕉国产在线看| 婷婷精品国产亚洲av| 欧美日韩亚洲综合一区二区三区_| 精品国产乱子伦一区二区三区| 一区福利在线观看| 午夜福利成人在线免费观看| 亚洲国产毛片av蜜桃av| 老司机福利观看| 国产精品综合久久久久久久免费| 青草久久国产| 777久久人妻少妇嫩草av网站| 丝袜美腿诱惑在线| 在线国产一区二区在线| 日韩三级视频一区二区三区| 一本精品99久久精品77| 精华霜和精华液先用哪个| 99热6这里只有精品| 久久久久久免费高清国产稀缺| 久久中文字幕人妻熟女| 丁香六月欧美| 免费在线观看影片大全网站| 性色av乱码一区二区三区2| 99国产精品99久久久久| 黄色视频不卡| 亚洲在线自拍视频| 久久精品夜夜夜夜夜久久蜜豆 | 午夜亚洲福利在线播放| 日本一本二区三区精品| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲精华国产精华精| 久久精品国产综合久久久| 中文字幕人妻丝袜一区二区| 国产又黄又爽又无遮挡在线| 久久精品aⅴ一区二区三区四区| 一本精品99久久精品77| 国产成年人精品一区二区| 亚洲国产精品久久男人天堂| 中文字幕人成人乱码亚洲影| 女人被狂操c到高潮| 午夜福利在线在线| 麻豆成人午夜福利视频| 精品少妇一区二区三区视频日本电影| 性色av乱码一区二区三区2| 欧美乱码精品一区二区三区| 久久久久久国产a免费观看| 黑人欧美特级aaaaaa片| 岛国视频午夜一区免费看| 黄色丝袜av网址大全| avwww免费| 后天国语完整版免费观看| 婷婷精品国产亚洲av| 日韩欧美免费精品| 亚洲成人免费电影在线观看| 大型av网站在线播放| 精品不卡国产一区二区三区| 99在线视频只有这里精品首页| 97碰自拍视频| 精品欧美一区二区三区在线| 色综合亚洲欧美另类图片| 亚洲 欧美一区二区三区| 精品久久久久久久毛片微露脸| 搡老熟女国产l中国老女人| 国产成人精品久久二区二区91| 又黄又粗又硬又大视频| 日本免费a在线| 中文字幕精品免费在线观看视频| 国产精品爽爽va在线观看网站 | 精品免费久久久久久久清纯| 国产精品一区二区精品视频观看| 国产熟女午夜一区二区三区| 91成人精品电影| 国产区一区二久久| 人人妻人人澡欧美一区二区| 午夜福利成人在线免费观看| 久久狼人影院| 特大巨黑吊av在线直播 | 国产精品美女特级片免费视频播放器 | 国产高清视频在线播放一区| 国产成人影院久久av| 精品少妇一区二区三区视频日本电影| 亚洲成av片中文字幕在线观看| 97超级碰碰碰精品色视频在线观看| 黄色成人免费大全| 一个人免费在线观看的高清视频| 国产精品永久免费网站| 亚洲五月色婷婷综合| 亚洲精品美女久久av网站| 亚洲成a人片在线一区二区| 国产黄片美女视频| av欧美777| 满18在线观看网站| 国产又黄又爽又无遮挡在线| 啦啦啦观看免费观看视频高清| 嫩草影视91久久| 99国产精品一区二区蜜桃av| 亚洲人成网站高清观看| 欧美日本亚洲视频在线播放| 久久久久国内视频| 国产1区2区3区精品| 国产一区在线观看成人免费| 色播亚洲综合网| 国产又黄又爽又无遮挡在线| 嫩草影视91久久| 18禁黄网站禁片午夜丰满| 久久久久久九九精品二区国产 | 91成人精品电影| 国产麻豆成人av免费视频| 日本免费一区二区三区高清不卡| 欧美色视频一区免费| 久热爱精品视频在线9| 1024视频免费在线观看| 日本成人三级电影网站| 桃红色精品国产亚洲av| 国产精品98久久久久久宅男小说| cao死你这个sao货| 老司机深夜福利视频在线观看| 51午夜福利影视在线观看| 搡老熟女国产l中国老女人| 免费看美女性在线毛片视频| 又黄又粗又硬又大视频| 中文字幕另类日韩欧美亚洲嫩草| 岛国在线观看网站| 日本a在线网址| 老司机在亚洲福利影院| 欧美在线一区亚洲| 中文资源天堂在线| 精品国产乱码久久久久久男人| 国产免费男女视频| 性色av乱码一区二区三区2| 两性午夜刺激爽爽歪歪视频在线观看 | 国产亚洲精品av在线| 中文字幕另类日韩欧美亚洲嫩草| 欧美精品啪啪一区二区三区| 国产成人精品无人区| 黄色成人免费大全| 久久九九热精品免费| 成人特级黄色片久久久久久久| 禁无遮挡网站| 免费在线观看影片大全网站| 最近在线观看免费完整版| 每晚都被弄得嗷嗷叫到高潮| 男人操女人黄网站| 一进一出抽搐gif免费好疼| 亚洲成a人片在线一区二区| 午夜影院日韩av| 国产精华一区二区三区| 999久久久精品免费观看国产| 亚洲成av人片免费观看| 精品国产超薄肉色丝袜足j| 最近在线观看免费完整版| 国产成人欧美| 国产亚洲av嫩草精品影院| 国产午夜福利久久久久久| 一进一出抽搐动态| 日韩精品青青久久久久久| xxxwww97欧美| 一区二区三区激情视频| 黄片大片在线免费观看| 精品国产亚洲在线| 欧美成狂野欧美在线观看| 欧美日韩黄片免| 别揉我奶头~嗯~啊~动态视频| 欧美一区二区精品小视频在线| 免费在线观看黄色视频的| av天堂在线播放| 久99久视频精品免费| 欧美三级亚洲精品| 在线观看日韩欧美| 亚洲美女黄片视频| 欧美另类亚洲清纯唯美| 91麻豆精品激情在线观看国产| 人妻丰满熟妇av一区二区三区| 国产一区二区激情短视频| 黑人操中国人逼视频| 中文字幕另类日韩欧美亚洲嫩草| 国产精品国产高清国产av| 欧美成人一区二区免费高清观看 | 欧美成人一区二区免费高清观看 | 欧美丝袜亚洲另类 | 一级毛片精品| 久久久久久大精品| 搡老熟女国产l中国老女人| 国内毛片毛片毛片毛片毛片| 久久精品亚洲精品国产色婷小说| 亚洲欧美一区二区三区黑人| 欧美绝顶高潮抽搐喷水| 日日干狠狠操夜夜爽| 国语自产精品视频在线第100页| www.www免费av| 老熟妇仑乱视频hdxx| 在线av久久热| 国产av一区二区精品久久| svipshipincom国产片| 无遮挡黄片免费观看| 99在线视频只有这里精品首页| 99国产极品粉嫩在线观看| 女人被狂操c到高潮| 久热这里只有精品99| 精品久久久久久久久久免费视频| 国产精品日韩av在线免费观看|