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

    怒江上游流域降雪識別及其演變趨勢和原因分析

    2018-03-27 21:29:56劉少華嚴(yán)登華秦天玲翁白莎盧亞靜
    水利學(xué)報 2018年2期
    關(guān)鍵詞:降雪量相態(tài)降雪

    劉少華 ,嚴(yán)登華 ,王 浩 ,秦天玲 ,翁白莎 ,盧亞靜

    (1.流域水循環(huán)模擬與調(diào)控國家重點實驗室,北京 100038;2.中國水利水電科學(xué)研究院 水資源研究所,北京 100038;3.長江勘測規(guī)劃設(shè)計研究院,湖北 武漢 430010)

    1 研究背景

    降水是水循環(huán)過程中最活躍的要素,其相態(tài)變化直接影響地表能量和水循環(huán)過程[1-2]。降雪在地表的累積不僅增加了地表反射率,改變了地表的能量過程;而且延遲了產(chǎn)匯流時間,改變水循環(huán)過程[3-4]。青藏高原是氣候變化的敏感區(qū),典型的高寒氣候使得其降雪量在降水中占據(jù)了較大的比例,在以增溫為背景的氣候變化作用下其降雪時空分布均發(fā)生了顯著變化,進而改變了融雪過程的時空特征,對其水循環(huán)過程產(chǎn)生了深刻的影響[5]。

    當(dāng)前基于降水形成機理的相態(tài)識別研究以氣象學(xué)研究為主,即以大氣層溫度廓線為基礎(chǔ),根據(jù)大氣層溫度及厚度的梯度變化識別降水過程中的相態(tài)變化[6-8]。如許愛華等[9]對2005年我國的一次寒潮事件研究表明,925 hPa氣溫<-2℃則可作為降雪的預(yù)報判據(jù)。李江波等[10]對2007年河北地區(qū)一次春季強寒潮研究指出:0℃層高度低于950 hPa、1000 hPa溫度<2℃且925 hPa溫度<-2℃時,降雨轉(zhuǎn)變?yōu)榻笛┗蛘哂陫A雪;高洋等[11]對2008年我國南方地區(qū)的凍雨過程研究指出,700 hPa附近溫度最高達到4℃以上,有利于形成降雨;在850 hPa以下氣溫明顯<0℃,有利形成降雪。然而,由于氣象系統(tǒng)的復(fù)雜性和大氣層相關(guān)資料的稀缺性,以及氣象學(xué)與水文學(xué)研究時空尺度的差異性,基于氣象學(xué)研究的降雪識別尚難以應(yīng)用于水循環(huán)研究。水文學(xué)研究中通常采用統(tǒng)計學(xué)方法對降水相態(tài)進行識別,常用方法可分為以下兩種:(1)單溫度閾值法,當(dāng)氣溫高于臨界溫度時認(rèn)為降水以液態(tài)(降雨)發(fā)生,當(dāng)氣溫低于臨界溫度是認(rèn)為降水以固態(tài)(降雪)發(fā)生[12];(2)雙溫度閾值法,當(dāng)氣溫高于臨界高溫時認(rèn)為降雨概率是100%,當(dāng)氣溫低于臨界低溫時則認(rèn)為降雪概率是100%,當(dāng)氣溫介于臨界溫度之間時認(rèn)為降雪概率隨溫度呈線性變化[13-14],或根據(jù)臨界溫度區(qū)間降雪事件與氣溫的統(tǒng)計關(guān)系確定降雪事件發(fā)生概率[15-2]。由于不同區(qū)域氣候條件和下墊面系統(tǒng)的差異性,單溫度閾值難以適應(yīng)不同區(qū)域的氣候特征,加之降水過程受局部環(huán)境影響顯著,單溫度閾值降雪識別必然會導(dǎo)致較大誤差[16]。雙溫度閾值法在當(dāng)前降雪識別研究中應(yīng)用較廣泛,相關(guān)學(xué)者根據(jù)其研究區(qū)域的氣候特征提出了多種雙溫度閾值法,其中指數(shù)方程形式的雙溫度閾值法在青藏高原地區(qū)的降雪識別中精度較高[17]。

    本文選取青藏高原東部的怒江上游流域為研究區(qū)域,以流域氣象站點觀測氣溫和降水相態(tài)資料為基礎(chǔ),采用指數(shù)方程形式的雙溫度閾值方法進行流域降雪的識別、驗證及趨勢分析,并通過設(shè)置氣溫和降水的控制情景方案,量化識別氣溫和降水的變化對流域降雪量變化的貢獻。

    2 研究區(qū)概況及數(shù)據(jù)

    2.1 研究區(qū)概況怒江上游流域位于青藏高原東南部,流域界于東經(jīng)91°10′—96°20′和北緯29°55′—32°50′之間、面積約76 000 km2、平均海拔在4 000 m以上。流域大部分區(qū)域?qū)俑咴瓉喓畮Ъ撅L(fēng)濕潤與半濕潤氣候區(qū),多年平均日氣溫為0.26℃,月平均氣溫從11月至次年3月均在0℃以下。流域多年平均降水量556 mm,其中80%的降水量集中在6—9月,多年平均降雪量不足降水總量的20%,但由于流域內(nèi)積雪分布廣泛,積雪融水是流域春季徑流的重要組成部分,且在溫升作用下發(fā)生了顯著的變化。識別怒江上游流域降雪過程,量化氣候變化對降雪量的影響,對于提高積雪融水模擬精度、明晰徑流的年內(nèi)分配特征、支撐變化環(huán)境下流域水資源綜合管理具有重要意義。

    2.2 數(shù)據(jù)資料我國公開降水相態(tài)資料相對匱乏,由國家氣象信息中心提供的日降水資料集(V3.0)中僅1979年以前降水資料中標(biāo)記了降水相態(tài)(降雨/降雪/雨夾雪)。怒江上游流域氣象站點稀疏,是典型的缺資料地區(qū)[18],流域內(nèi)部及周邊僅有6個氣象站具有降水相態(tài)記載(圖1),統(tǒng)計各氣象站點年均氣溫、降水量以及降水相變的臨界氣溫(根據(jù)每0.5℃區(qū)間氣溫與降雪比例關(guān)系得到)見表1。由于單個站點系列較短、樣本代表性差,研究綜合流域6個氣象站的1979年以前逐日降水相態(tài)和平均氣溫資料,統(tǒng)計臨界氣溫內(nèi)(-4~8℃)氣溫與降雪比例關(guān)系,并擬合得到流域降雪識別指數(shù)方程。

    圖1 怒江上游流域位置及氣象站點分布

    3 研究方法

    3.1 降雪量識別方法氣象站點的降水相態(tài)標(biāo)記資料中無法得知雨夾雪中降雪的含量,統(tǒng)計研究區(qū)氣象站點不同氣溫下降雪量占降雪和降水總量的比值(不考慮雨夾雪情況下的氣溫與降雪比例關(guān)系),繪制其降雪比例與日平均氣溫的關(guān)系(圖2(a))。當(dāng)溫度小于臨界低溫(Tmin)時,降水完全以降雪形態(tài)發(fā)生;當(dāng)溫度大于臨界高溫(Tmax)時,降水完全以降雨形態(tài)發(fā)生;臨界溫度內(nèi)(Tmin~Tmax),降雪比例與日平均氣溫成指數(shù)函數(shù)關(guān)系[2-17],可概化為如下方程:

    表1 氣象站點相關(guān)信息

    式中:PS為降雪比例;T為日平均氣溫;a和b均為指數(shù)方程的經(jīng)驗參數(shù),可以根據(jù)研究區(qū)域氣溫與降雪比例的統(tǒng)計關(guān)系采用最小二乘法估算得到。

    在此基礎(chǔ)上,假定相同溫度下雨夾雪中降雪比例也滿足該指數(shù)方程,并統(tǒng)計研究區(qū)域不同氣溫下降雪量、雨夾雪和降雨量占總降水量的比值,繪制降水各相態(tài)與日平均氣溫的關(guān)系(圖2(b))。進而,把各氣象站點雨夾雪中估計降雪量和觀測降雪量之和作為其總降雪量,并以此為據(jù)驗證采用指數(shù)方程識別總降水量得到的降雪量。

    圖2 降水各相態(tài)比例與氣溫關(guān)系

    3.2 降雪量變化趨勢檢驗根據(jù)各氣象站點逐日氣溫和降水量,采用指數(shù)方程識別其逐日降雪量,并累計得到其年降雪量過程。通過線性回歸分析各氣象站點年平均氣溫、年降水量和年降雪量變化趨勢,并采用Mann-Kendall檢驗法對其變化趨勢進行顯著性檢驗。Mann-Kendall檢驗法是氣象水文分析中最常用的一種非參數(shù)檢驗方法[19],其計算公式如下:

    式中:x為進行趨勢檢驗的時間序列;n為時間序列的長度。統(tǒng)計量Z符合正態(tài)分布,在給定的α置信水平上,當(dāng),則拒絕零假設(shè),表明該時間序列的變化趨勢滿足顯著性檢驗;Z>0,表明序列具有上升趨勢,反之為下降趨勢。

    3.3 氣溫和降水對降雪量影響的分解降雪是降水在一定氣候條件下的特殊形態(tài),降雪量的變化是降水量變化的直接體現(xiàn),但同時也受到以增溫為主導(dǎo)的氣候變化的影響。為了量化氣候變化背景下氣溫和降水的變化對降雪量變化的貢獻,根據(jù)研究區(qū)域氣候變化的突變點把其氣候要素時間序列劃分為兩個階段:氣候變化前和氣候變化后,并以氣候變化前多年平均降雪量作為基準(zhǔn)值(S0),以氣候變化后多年平均降雪量作為氣候變化后降雪量(SC)。同時,設(shè)置不同氣溫和降水量控制方案,分別研究氣溫和降水量變化對降雪量變化的貢獻。

    (1)氣溫控制方案。對氣候變化后逐日平均氣溫去趨勢化修正得到無趨勢逐日平均氣溫系列。具體步驟:①計算得到氣候變化后年平均氣溫的增長速率;②根據(jù)年逐日平均氣溫減去該年平均氣溫增量(年平均氣溫增長率乘以年序號),得到相應(yīng)年份去趨勢化逐日氣溫系列。根據(jù)去趨勢化逐日平均氣溫和降水量,重新應(yīng)用指數(shù)方程對氣候變化后逐日降雪量進行估計,得到控制氣溫條件下多年平均降雪量(ST)。

    (2)降水控制方案。由于指數(shù)方程可知,當(dāng)溫度不變時降水中降雪的比例是定值,即統(tǒng)計意義上溫度不變降雪率(降雪量/降水量)相同。根據(jù)氣候變化后多年平均降雪率乘以氣候變化前多年平均降水量得到控制降水條件下多年平均降雪量(SP)。根據(jù)以上控制條件的降雪量,可分別得到氣溫、降水量以及兩者共同影響對降雪量變化的貢獻率:

    其中,CT、CP和CP&T分別是氣溫、降水以及兩者共同影響對降雪量變化的貢獻率,需要注意的是此處貢獻率是指由于氣溫和降水的改變導(dǎo)致的降雪變化率(相對于氣候變化前多年平均降雪量S0),所以其理論取值范圍為[-1,+∞)(負(fù)值表示該要素的增加/升高對降雪量的增加具有負(fù)貢獻),一般情況下三者之和不為1。

    4 案例應(yīng)用

    圖3 各氣象站點總降雪量與估計年降雪量

    圖4 各氣象站點多年月總降雪量與估計月降雪量

    4.1 降雪量的識別及驗證綜合怒江上游流域6個氣象站點1979年以前氣溫和降水相態(tài)資料,統(tǒng)計其氣溫與降雪比例關(guān)系,采用指數(shù)方程對其進行擬合,以流域臨界氣溫范圍(-4~8℃)作為擬合區(qū)間,擬合計算步長為0.5℃。根據(jù)最小二乘法估算擬合指數(shù)方程參數(shù)a和b的值分別為-6.11和1.81,殘差平方和為0.0038。通過指數(shù)方程估計流域各氣象站點1979年以前逐日降雪量,并統(tǒng)計其歷年估計降雪量與總降雪量關(guān)系(圖3)可知:除那曲站外各氣象站點估計年降雪量相關(guān)系數(shù)均在0.95以上,那曲站估計年降雪量相關(guān)系數(shù)略低,但仍表現(xiàn)出與總降雪量較好的一致性。各氣象站點估計年降雪量的累計相對誤差范圍均在5%以內(nèi),僅那曲站相對誤差為正值。其原因是流域臨界氣溫與氣象站點實際臨界氣溫存在一定偏差:流域臨界氣溫較那曲站臨界氣溫偏高(指數(shù)方程向右便移),導(dǎo)致那曲站估計降雪量偏高;相反,流域臨界氣溫較其他站點臨界氣溫偏低(指數(shù)方程向左偏移),導(dǎo)致其估計降雪量偏低??傮w而言,該降雪識別方法可能低估怒江上游流域降雪量,但整體能夠較好的估計各氣象站點的年降雪量。其次,統(tǒng)計氣象站點各月多年平均估計降雪量(圖4)顯示,各氣象站點降雪量主要集中在春夏(4—6月)和夏秋交替月份(9—10月),其原因是怒江上游流域春季和冬季降水量較少、夏季氣溫較高,不利于形成降雪。而春夏和夏秋交替月份具有一定規(guī)模的降水量,且其平均氣溫均在臨界溫度內(nèi),容易形成降雪。春夏交替季節(jié)各氣象站點降雪量出現(xiàn)的峰值從流域東南向西北具有一定滯后現(xiàn)象,其中丁青站降雪峰值出現(xiàn)最早在4月份,安多站降雪峰值出現(xiàn)最晚在6月份,其他站點降雪峰值均出現(xiàn)在5月份,且降雪峰值月份的估計降雪量明顯偏低。其可能原因是怒江上游流域春夏交替季節(jié),伴隨氣溫升高季風(fēng)攜帶大量水汽由東南向西北快速推進,使得流域西北部雨季較東南部滯后;同時,高海拔地形作用下該月份冷暖氣流交匯頻繁易形成區(qū)域短歷時降雪,該降雪過程受氣溫影響較?。?0],導(dǎo)致基于氣溫和降雪比例統(tǒng)計關(guān)系的指數(shù)方程普遍低估流域降雪量。

    4.2 降雪量演變趨勢分析降雪是大氣在一定的濕度、溫度以及區(qū)域環(huán)流場等綜合條件下形成的固態(tài)降水,其形成條件雖復(fù)雜多樣,但并不隨氣象要素演變而變化。根據(jù)氣溫和降雪比例統(tǒng)計關(guān)系擬合得到的降雪識別指數(shù)方程則是對降雪形成條件的定量概化。因而,可根據(jù)流域氣象站點1979年以前資料獲得的流域降雪識別指數(shù)方程對各氣象站點1980—2016年逐日降雪量進行識別,并采用線性回歸分析1979年前后各氣象站點年平均氣溫、年降水量、年降雪量變化趨勢(圖5)。結(jié)果可知:流域整體上氣溫呈升高趨勢,且1979年以后氣溫升高趨勢更為明顯;1979年以前流域整體降水量變化趨勢不明顯,各氣象站點降水趨勢不完全一致,1979年以后安多和雜多站降水量呈減少趨勢,其余氣象站點的降水量均呈現(xiàn)出增加趨勢。與此同時,1979年以前流域降雪量呈增加趨勢,僅那曲站降雪量表現(xiàn)出微弱的減少趨勢,1979年以后所有氣象站點降雪量均表現(xiàn)出減少趨勢。根據(jù)Mann-Kendall趨勢檢驗結(jié)果(表2)可知,各氣象站點1980—2016年平均氣溫、年降水量和年降雪量的變化趨勢與線性回歸趨勢基本一致,其中平均氣溫均呈顯著增加趨勢,但降水量和降雪量變化趨勢均不顯著。

    表2 各氣象站點年1980—2016年均氣溫、年降水量和年降雪量變化趨勢的顯著性檢驗結(jié)果(Z值)

    圖5 各氣象站點估計年降雪量、年平均氣溫和年降水量變化趨勢

    4.3 氣溫和降水對降雪量影響分析研究采用現(xiàn)有研究普遍認(rèn)可的1980年作為怒江上游流域氣候變化的轉(zhuǎn)折點[21-24],采用已獲得的流域降雪指數(shù)識別方程,對氣象站點1980年后不同控制方案下的多年平均降雪量進行估計,并量化識別氣候變化后流域氣象站點氣溫和降水變化對降雪量變化的貢獻率(圖6),結(jié)果表明:各氣象站點1979年以后多年平均降雪量較之前多年平均降雪量均有所增加。各氣象站點氣溫升高對于降雪量表現(xiàn)出明顯的減少作用。除丁青站外其他氣象站點1979年以后多年平均降水量均大于之前,并對其降雪量具有明顯的增加作用。丁青站1979年以后多年平均降水量略低于之前,因而其對降雪量呈現(xiàn)出了減少作用。怒江上游流域氣溫和降水對降雪量的影響呈東南向西北方向逐漸增強分布規(guī)律,其可能原因是流域的平均氣溫和年降水量均呈西南向東北遞減分布特征,氣溫越低、降水量越少降雪量對氣溫和降水變化的響應(yīng)越敏感。氣溫和降水量共同影響對各氣象站點降雪量均表現(xiàn)出增加作用。原因是氣候變化背景下流域低溫條件下降水概率增加(圖7),對于各氣象站點降雪量增加具有顯著促進作用。

    圖6 各氣象站點1979年以后多年平均降雪量增加率及其分解

    圖7 1979年前后降水隨氣溫變化的概率分布

    5 結(jié)語

    (1)基于降雪比例和日平均氣溫的降雪識別指數(shù)方程,能夠在年和月尺度上較好估計怒江上游流域各氣象站點降雪量。各氣象站點估計降雪量與總降雪量相關(guān)系數(shù)均在0.86以上,相對誤差均在5%以內(nèi)。

    (2)1979年以后,怒江上游流域各氣象站點氣溫均呈顯著增加趨勢,各氣象站點降水量和降雪量變化趨勢并不一致,流域整體上年降水量和年降雪量分別呈增加和減少趨勢,但趨勢并不顯著。

    (3)怒江上游流域氣象站點氣溫升高具有減少降雪量的作用,且作用強度從流域東南向西北逐漸增強。氣候變化背景下,低溫條件下降水概率增加對各氣象站點的降雪量的增加具有顯著促進作用,但1979年以后氣溫的顯著增加對降雪量的影響占據(jù)主導(dǎo)作用,使得流域降雪量整體呈減少趨勢。

    [1]LOTH B,GRAF H F,OBERHUBER J M.Snow cover model for global climate simulations[J].Journal of Geo?physical Research:Atmospheres,1993,98(D6):10451-10464 .

    [2]DING B,YANG K,QIN J,et al.The dependence of precipitation types on surface elevation and meteorological conditions and its parameterization[J].Journal of Hydrology,2014,513(11):154-163 .

    [3]CLARK M P,SLATER A G,BARRETT A P,et al.Assimilation of snow covered area information into hydrologic and land-surface models[J].Advances in Water Resources,2006,29(8):1209-1221 .

    [4]BOX J,F(xiàn)ETTWEIS X,STROEVE J,et al.Greenland ice sheet albedo feedback:thermodynamics and atmospher?ic drivers[J].The Cryosphere,2012,6(4):821-839 .

    [5]姚檀棟,秦大河,沈永平,等.青藏高原冰凍圈變化及其對區(qū)域水循環(huán)和生態(tài)條件的影響[J].自然雜志,2013,35(3):179-186 .

    [6]漆梁波,張瑛.中國東部地區(qū)冬季降水相態(tài)的識別判據(jù)研究[J].氣象,2012,38(1):96-102.

    [7]張琳娜,郭銳,曾劍,等 .北京地區(qū)冬季降水相態(tài)的識別判據(jù)研究[J].高原氣象,2013,32(6):1780-1786.

    [8]徐輝,宗志平 .一次降水相態(tài)轉(zhuǎn)換過程中溫度垂直結(jié)構(gòu)特征分析[J].高原氣象,2014,33(5):1272-1280.

    [9]許愛華,喬林,詹豐興,等.2005年3月一次寒潮天氣過程的診斷分析[J].氣象,2006,32(3):49-55.

    [10]李江波,李根娥,裴雨杰,等 .一次春季強寒潮的降水相態(tài)變化分析[J].氣象,2009,35(7):87-94.

    [11]高洋,吳統(tǒng)文,陳葆德.2008年1月我國南方凍雨過程的熱力異常及其形成原因[J].高原氣象,2011,30(6):1526-1533.

    [12]YANG Z L,DICKINSON R E,ROBOCK A,et al.Validation of the snow submodel of the biosphere-atmosphere transfer scheme with Russian snow cover and meteorological observational data[J].Journal of Climate,1997,10(2):353-373.

    [13]WIGMOSTA M S,VAIL L W,LETTENMAIER D Pv.A distributed hydrology-vegetation model for complex ter?rain[J].Water Resources Research,1994,30(6):1665-1679 .

    [14]CHEN R,LU S,KANG E,et al.A distributed water-heat coupled model for mountainous watershed of an inland river basin of Northwest China(I)model structure and equations[J].Environmental Geology,2008,53(6):1299-1309.

    [15]DAI A.Temperature and pressure dependence of the rain-snow phase transition over land and ocean[J].Geo?physical Research Letters,2008,35(12):62-77 .

    [16]韓春壇,陳仁升,劉俊峰,等.固液態(tài)降水分離方法探討[J].冰川凍土,2010,32(2):249-256.

    [17]LIU S,YAN D,QIN T,et al.Precipitation phase separation schemes in the Naqu River basin,eastern Tibetan plateau[J].Theoretical&Applied Climatology,2016.doi:10.1007/s00704-016-1967-7 .

    [18]劉昌明,白鵬,王中根,等 .稀缺資料流域水文計算若干研究:以青藏高原為例[J].水利學(xué)報,2016,47(3):272-282.

    [19]劉少華,嚴(yán)登華,翁白莎,等 .近50a中國≥10℃有效積溫時空演變[J].干旱區(qū)研究,2013,30(4):689-696.

    [20]徐祥德,陶詩言,王繼志,等.青藏高原-季風(fēng)水汽輸送“大三角扇型”影響域特征與中國區(qū)域旱澇異常的關(guān)系[J].氣象學(xué)報,2002,60(3):257-266.

    [21]馮新靈,羅隆誠,馮自立.中國近50年降水變化趨勢及突變的Hurst指數(shù)試驗[J].干旱區(qū)地理,2009,32(6):859-866.

    [22]尹云鶴,吳紹洪,陳剛.1961—2006年我國氣候變化趨勢與突變的區(qū)域差異[J].自然資源學(xué)報,2009,24(12):2147-2157.

    [23]劉春蓁,占車生,夏軍,等.關(guān)于氣候變化與人類活動對徑流影響研究的評述[J].水利學(xué)報,2014,45(4):379-393.

    [24]朱麗華,范廣洲,華維.全球變暖背景下青藏高原夏季氣溫在對流層上下反相變化及其與降水和環(huán)流的關(guān)系[J].大氣科學(xué),2015,39(6):1250-1262.

    猜你喜歡
    降雪量相態(tài)降雪
    1970—2019年內(nèi)蒙古大興安嶺林區(qū)降雪特征分析
    1961—2016年中國天山不同級別降雪事件變化特征分析
    降雪量與積雪深度
    少兒科技(2018年7期)2018-08-04 22:35:11
    川西高原東部兩次連續(xù)強降雪過程對比分析
    SBS改性瀝青相態(tài)結(jié)構(gòu)的參數(shù)化表征方法
    2016年1月19~21日山南地區(qū)南部一次降雪過程分析
    西藏科技(2016年10期)2016-09-26 09:02:08
    四川省降水相態(tài)識別判據(jù)研究
    PS/PLA共混物的相態(tài)結(jié)構(gòu)及其發(fā)泡行為研究
    中國塑料(2015年4期)2015-10-14 01:09:21
    泡沫油WINPROP相態(tài)研究
    人工降雪
    中国美女看黄片| 咕卡用的链子| av又黄又爽大尺度在线免费看| 十八禁人妻一区二区| 久久久久精品人妻al黑| 久久精品亚洲熟妇少妇任你| 久久午夜综合久久蜜桃| 国产有黄有色有爽视频| 久久久久精品人妻al黑| 午夜福利一区二区在线看| 久久久久久久国产电影| 久久久精品区二区三区| 女人被躁到高潮嗷嗷叫费观| 丰满饥渴人妻一区二区三| 亚洲专区国产一区二区| 12—13女人毛片做爰片一| 女人高潮潮喷娇喘18禁视频| 一边摸一边做爽爽视频免费| 桃花免费在线播放| 国产精品自产拍在线观看55亚洲 | 热re99久久精品国产66热6| 亚洲精品国产色婷婷电影| 国产xxxxx性猛交| 精品少妇久久久久久888优播| 黑人欧美特级aaaaaa片| 各种免费的搞黄视频| 日本一区二区免费在线视频| 日本vs欧美在线观看视频| 中文字幕色久视频| 欧美另类亚洲清纯唯美| av天堂久久9| 精品国产乱子伦一区二区三区 | 免费观看av网站的网址| 亚洲中文字幕日韩| 国产老妇伦熟女老妇高清| 久久国产亚洲av麻豆专区| 国产成人精品久久二区二区免费| 人妻人人澡人人爽人人| 一本大道久久a久久精品| 色婷婷久久久亚洲欧美| 五月天丁香电影| 日本猛色少妇xxxxx猛交久久| 久久女婷五月综合色啪小说| 99精品欧美一区二区三区四区| 巨乳人妻的诱惑在线观看| 丁香六月天网| 美女中出高潮动态图| 亚洲视频免费观看视频| 国产精品.久久久| 亚洲av欧美aⅴ国产| 亚洲天堂av无毛| 777米奇影视久久| 两个人看的免费小视频| 黄色a级毛片大全视频| 久久香蕉激情| 精品亚洲成a人片在线观看| 91精品国产国语对白视频| 亚洲第一欧美日韩一区二区三区 | 免费日韩欧美在线观看| 亚洲欧美精品综合一区二区三区| 婷婷色av中文字幕| 日本欧美视频一区| 日本黄色日本黄色录像| 性色av乱码一区二区三区2| 最近最新免费中文字幕在线| 亚洲男人天堂网一区| 欧美 日韩 精品 国产| 成年人午夜在线观看视频| 日韩 亚洲 欧美在线| 丁香六月欧美| 国产精品久久久久成人av| 少妇 在线观看| 新久久久久国产一级毛片| 国产精品香港三级国产av潘金莲| 久久久久久人人人人人| 俄罗斯特黄特色一大片| 国产欧美日韩一区二区三区在线| 久久精品国产a三级三级三级| 巨乳人妻的诱惑在线观看| 欧美黑人欧美精品刺激| 肉色欧美久久久久久久蜜桃| 男人操女人黄网站| 欧美黄色淫秽网站| 国产麻豆69| 日韩制服骚丝袜av| 91成年电影在线观看| 脱女人内裤的视频| 他把我摸到了高潮在线观看 | 久久久精品区二区三区| av国产精品久久久久影院| 欧美少妇被猛烈插入视频| 夜夜夜夜夜久久久久| 亚洲天堂av无毛| 男女之事视频高清在线观看| 日韩中文字幕欧美一区二区| 夜夜夜夜夜久久久久| 最近最新免费中文字幕在线| 成年人黄色毛片网站| 久久久精品国产亚洲av高清涩受| 可以免费在线观看a视频的电影网站| 欧美黄色淫秽网站| 大陆偷拍与自拍| 电影成人av| 菩萨蛮人人尽说江南好唐韦庄| 在线看a的网站| 在线观看免费午夜福利视频| 99久久精品国产亚洲精品| 国产成人av激情在线播放| 国产亚洲午夜精品一区二区久久| 亚洲情色 制服丝袜| 又黄又粗又硬又大视频| 在线观看www视频免费| 国产精品久久久久久精品古装| 亚洲少妇的诱惑av| 男女无遮挡免费网站观看| 欧美在线黄色| 久久av网站| 成人亚洲精品一区在线观看| 这个男人来自地球电影免费观看| 国产精品亚洲av一区麻豆| 久久久久国产一级毛片高清牌| 国产区一区二久久| 国产又爽黄色视频| 激情视频va一区二区三区| 12—13女人毛片做爰片一| 午夜福利视频在线观看免费| 97在线人人人人妻| 黄网站色视频无遮挡免费观看| 老司机靠b影院| 久久毛片免费看一区二区三区| 国产亚洲欧美在线一区二区| 性色av乱码一区二区三区2| 在线观看免费日韩欧美大片| 国产成人免费观看mmmm| 久久国产精品人妻蜜桃| 蜜桃国产av成人99| 一个人免费在线观看的高清视频 | 亚洲自偷自拍图片 自拍| 亚洲熟女毛片儿| 免费日韩欧美在线观看| 婷婷丁香在线五月| 久久国产精品大桥未久av| 国产深夜福利视频在线观看| 亚洲av成人一区二区三| 少妇人妻久久综合中文| bbb黄色大片| 日韩有码中文字幕| 最黄视频免费看| 精品熟女少妇八av免费久了| 日韩大码丰满熟妇| 肉色欧美久久久久久久蜜桃| 日韩视频一区二区在线观看| 五月开心婷婷网| 人妻 亚洲 视频| 国产精品一二三区在线看| 50天的宝宝边吃奶边哭怎么回事| 欧美日韩亚洲综合一区二区三区_| 一边摸一边抽搐一进一出视频| 国产精品免费视频内射| 久久久精品区二区三区| 巨乳人妻的诱惑在线观看| 亚洲男人天堂网一区| 90打野战视频偷拍视频| 国产成人精品久久二区二区91| 国产成人精品在线电影| 一二三四社区在线视频社区8| 黄片小视频在线播放| 国产三级黄色录像| 最新的欧美精品一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 精品久久久久久电影网| 老汉色av国产亚洲站长工具| 欧美xxⅹ黑人| 爱豆传媒免费全集在线观看| 国产亚洲午夜精品一区二区久久| 亚洲色图综合在线观看| 日韩免费高清中文字幕av| 水蜜桃什么品种好| 色婷婷久久久亚洲欧美| 自拍欧美九色日韩亚洲蝌蚪91| 黑人欧美特级aaaaaa片| 久久精品亚洲熟妇少妇任你| 在线观看www视频免费| 在线观看一区二区三区激情| 亚洲人成77777在线视频| 18禁黄网站禁片午夜丰满| 91老司机精品| 一级毛片精品| www.熟女人妻精品国产| 夫妻午夜视频| 国产成人av激情在线播放| 五月天丁香电影| 水蜜桃什么品种好| 三上悠亚av全集在线观看| 精品人妻在线不人妻| 汤姆久久久久久久影院中文字幕| 一本综合久久免费| 亚洲男人天堂网一区| 女人精品久久久久毛片| 国产精品熟女久久久久浪| 国产亚洲av片在线观看秒播厂| 一区福利在线观看| 中文字幕制服av| 国产亚洲av片在线观看秒播厂| 久久天堂一区二区三区四区| 人人妻人人爽人人添夜夜欢视频| 亚洲视频免费观看视频| 亚洲,欧美精品.| 国产精品av久久久久免费| 69精品国产乱码久久久| 视频区图区小说| 欧美xxⅹ黑人| av网站在线播放免费| 我要看黄色一级片免费的| 麻豆乱淫一区二区| 国产亚洲av片在线观看秒播厂| 宅男免费午夜| 成人亚洲精品一区在线观看| 国内毛片毛片毛片毛片毛片| 亚洲av国产av综合av卡| 亚洲av电影在线观看一区二区三区| 天天躁日日躁夜夜躁夜夜| 亚洲 欧美一区二区三区| 亚洲国产欧美一区二区综合| 99久久精品国产亚洲精品| 视频区图区小说| 免费看十八禁软件| 制服人妻中文乱码| 曰老女人黄片| 美女中出高潮动态图| 男女床上黄色一级片免费看| 老鸭窝网址在线观看| 中文字幕人妻丝袜一区二区| 动漫黄色视频在线观看| 久久国产精品男人的天堂亚洲| 久久精品国产亚洲av香蕉五月 | 欧美黑人精品巨大| 男女之事视频高清在线观看| 日韩中文字幕视频在线看片| 老司机亚洲免费影院| 日韩电影二区| 日韩欧美免费精品| 18禁裸乳无遮挡动漫免费视频| 欧美 亚洲 国产 日韩一| 在线观看人妻少妇| 女人久久www免费人成看片| 日本wwww免费看| 在线观看舔阴道视频| 欧美在线黄色| 大香蕉久久网| 下体分泌物呈黄色| 国产在线观看jvid| 五月开心婷婷网| 精品国产乱码久久久久久小说| 精品国内亚洲2022精品成人 | 12—13女人毛片做爰片一| 午夜精品国产一区二区电影| 成人影院久久| 国产精品熟女久久久久浪| 麻豆乱淫一区二区| 精品人妻在线不人妻| 日本a在线网址| bbb黄色大片| 亚洲精品av麻豆狂野| 婷婷色av中文字幕| av免费在线观看网站| e午夜精品久久久久久久| 精品一区二区三区av网在线观看 | 菩萨蛮人人尽说江南好唐韦庄| 桃花免费在线播放| 国产精品欧美亚洲77777| 精品少妇一区二区三区视频日本电影| 美女高潮喷水抽搐中文字幕| 99香蕉大伊视频| 美女扒开内裤让男人捅视频| 久久国产精品影院| 深夜精品福利| 欧美成狂野欧美在线观看| xxxhd国产人妻xxx| 99国产极品粉嫩在线观看| 日韩欧美一区视频在线观看| 十八禁高潮呻吟视频| 在线精品无人区一区二区三| 久久久久久久久久久久大奶| 一本—道久久a久久精品蜜桃钙片| 色播在线永久视频| 日韩精品免费视频一区二区三区| 制服诱惑二区| www日本在线高清视频| 另类亚洲欧美激情| tocl精华| av国产精品久久久久影院| 亚洲全国av大片| 亚洲欧美一区二区三区黑人| 黑人猛操日本美女一级片| 又紧又爽又黄一区二区| 少妇被粗大的猛进出69影院| 欧美xxⅹ黑人| 成人三级做爰电影| 国产精品熟女久久久久浪| 亚洲第一av免费看| 人妻 亚洲 视频| 亚洲全国av大片| a级片在线免费高清观看视频| 最近最新免费中文字幕在线| svipshipincom国产片| 激情视频va一区二区三区| 中文字幕人妻丝袜一区二区| 一区二区av电影网| 国产在视频线精品| 国产在线一区二区三区精| 大香蕉久久网| av超薄肉色丝袜交足视频| 成人影院久久| 亚洲全国av大片| 男女之事视频高清在线观看| √禁漫天堂资源中文www| av不卡在线播放| 午夜视频精品福利| 捣出白浆h1v1| 日本av免费视频播放| 人人澡人人妻人| 日韩制服丝袜自拍偷拍| 国产精品.久久久| 18禁黄网站禁片午夜丰满| 99久久综合免费| 老汉色∧v一级毛片| 亚洲av电影在线进入| 国产97色在线日韩免费| 女人爽到高潮嗷嗷叫在线视频| 五月天丁香电影| 男女无遮挡免费网站观看| 欧美精品一区二区大全| 国产欧美日韩综合在线一区二区| 国产99久久九九免费精品| 精品视频人人做人人爽| 一进一出抽搐动态| 我的亚洲天堂| 国产一区有黄有色的免费视频| 波多野结衣一区麻豆| 国产在视频线精品| 老汉色av国产亚洲站长工具| 日韩制服骚丝袜av| 亚洲午夜精品一区,二区,三区| 国产一卡二卡三卡精品| 69精品国产乱码久久久| 91精品伊人久久大香线蕉| 老熟妇仑乱视频hdxx| 在线亚洲精品国产二区图片欧美| 免费观看av网站的网址| 天堂8中文在线网| 久久精品国产a三级三级三级| 电影成人av| 国产成人a∨麻豆精品| a级片在线免费高清观看视频| 亚洲av男天堂| 欧美 日韩 精品 国产| 男女边摸边吃奶| 精品少妇内射三级| 另类亚洲欧美激情| 日韩一卡2卡3卡4卡2021年| 男人爽女人下面视频在线观看| 亚洲精品美女久久av网站| 亚洲av欧美aⅴ国产| 91精品伊人久久大香线蕉| 美女高潮喷水抽搐中文字幕| 亚洲精品国产一区二区精华液| 精品亚洲乱码少妇综合久久| 大香蕉久久成人网| 精品卡一卡二卡四卡免费| 国产欧美日韩一区二区三 | 亚洲精品国产一区二区精华液| 黄色视频在线播放观看不卡| 他把我摸到了高潮在线观看 | 最新在线观看一区二区三区| 一二三四在线观看免费中文在| 免费高清在线观看日韩| 天堂8中文在线网| 午夜免费观看性视频| 国产精品秋霞免费鲁丝片| 不卡一级毛片| 老鸭窝网址在线观看| 天天躁日日躁夜夜躁夜夜| 亚洲第一青青草原| 日韩有码中文字幕| 在线观看舔阴道视频| 国产又色又爽无遮挡免| 久久精品久久久久久噜噜老黄| 手机成人av网站| 大香蕉久久网| 91老司机精品| 啦啦啦免费观看视频1| 99香蕉大伊视频| 水蜜桃什么品种好| 捣出白浆h1v1| 中文精品一卡2卡3卡4更新| 成人国产一区最新在线观看| 欧美日韩成人在线一区二区| 岛国在线观看网站| 久久人妻福利社区极品人妻图片| 日日夜夜操网爽| 深夜精品福利| 欧美精品av麻豆av| 精品一区二区三卡| 精品国产乱码久久久久久小说| 亚洲成人免费电影在线观看| 中亚洲国语对白在线视频| 久久精品国产综合久久久| 国产免费现黄频在线看| 91大片在线观看| 飞空精品影院首页| 香蕉国产在线看| 精品少妇一区二区三区视频日本电影| 亚洲av欧美aⅴ国产| 不卡av一区二区三区| 99精国产麻豆久久婷婷| 精品少妇久久久久久888优播| 国产又色又爽无遮挡免| 精品久久蜜臀av无| 丝瓜视频免费看黄片| 国产av又大| 一二三四在线观看免费中文在| 人成视频在线观看免费观看| 高清av免费在线| 91字幕亚洲| 亚洲精品日韩在线中文字幕| 一本大道久久a久久精品| 黄色视频不卡| 午夜成年电影在线免费观看| 国产野战对白在线观看| 欧美精品一区二区免费开放| 国产精品99久久99久久久不卡| 亚洲性夜色夜夜综合| 欧美另类亚洲清纯唯美| 久久国产精品男人的天堂亚洲| 亚洲伊人久久精品综合| 女人被躁到高潮嗷嗷叫费观| 欧美成狂野欧美在线观看| 啦啦啦啦在线视频资源| 久久精品成人免费网站| 一本大道久久a久久精品| 丝袜美腿诱惑在线| 咕卡用的链子| 日韩熟女老妇一区二区性免费视频| 中文字幕最新亚洲高清| 男人添女人高潮全过程视频| 国产黄频视频在线观看| 一区二区三区乱码不卡18| 精品国产一区二区三区久久久樱花| 久久 成人 亚洲| 欧美日本中文国产一区发布| 亚洲欧美日韩高清在线视频 | 巨乳人妻的诱惑在线观看| 欧美黑人精品巨大| 亚洲免费av在线视频| 妹子高潮喷水视频| 高清黄色对白视频在线免费看| 青春草亚洲视频在线观看| 老熟妇仑乱视频hdxx| 国产精品偷伦视频观看了| 亚洲五月色婷婷综合| 啪啪无遮挡十八禁网站| 免费日韩欧美在线观看| 精品国产乱码久久久久久男人| 亚洲成人手机| 国产在线免费精品| 国产成人啪精品午夜网站| 免费少妇av软件| 狂野欧美激情性xxxx| www.自偷自拍.com| av福利片在线| 老司机靠b影院| 超色免费av| 一边摸一边做爽爽视频免费| 老司机影院毛片| 美女福利国产在线| 欧美老熟妇乱子伦牲交| 大片电影免费在线观看免费| 精品人妻1区二区| 国产一区二区三区av在线| 精品少妇一区二区三区视频日本电影| 老鸭窝网址在线观看| 亚洲专区国产一区二区| 久久久久国内视频| 啦啦啦啦在线视频资源| 69精品国产乱码久久久| 男人舔女人的私密视频| 亚洲精品自拍成人| 国产视频一区二区在线看| 中文字幕另类日韩欧美亚洲嫩草| 国产成+人综合+亚洲专区| 成年美女黄网站色视频大全免费| 国产成人免费无遮挡视频| h视频一区二区三区| 2018国产大陆天天弄谢| 在线观看一区二区三区激情| 亚洲激情五月婷婷啪啪| 激情视频va一区二区三区| 51午夜福利影视在线观看| 欧美日韩国产mv在线观看视频| 午夜福利在线免费观看网站| 高清欧美精品videossex| 51午夜福利影视在线观看| 97人妻天天添夜夜摸| 99国产精品一区二区蜜桃av | 亚洲熟女精品中文字幕| 国产日韩欧美视频二区| 欧美在线一区亚洲| 一级黄色大片毛片| 亚洲美女黄色视频免费看| 欧美成人午夜精品| 制服诱惑二区| 久久精品国产a三级三级三级| 十分钟在线观看高清视频www| 巨乳人妻的诱惑在线观看| 各种免费的搞黄视频| 国产亚洲一区二区精品| 一个人免费在线观看的高清视频 | 老熟女久久久| 老司机午夜十八禁免费视频| 人人妻,人人澡人人爽秒播| 欧美黄色淫秽网站| 亚洲专区中文字幕在线| h视频一区二区三区| 日韩中文字幕视频在线看片| 女人被躁到高潮嗷嗷叫费观| 考比视频在线观看| 国产精品自产拍在线观看55亚洲 | 精品免费久久久久久久清纯 | 国产成人av教育| 午夜福利视频在线观看免费| 国产亚洲欧美精品永久| 国产不卡av网站在线观看| 少妇人妻久久综合中文| 国产成人影院久久av| 一区二区三区四区激情视频| 女人久久www免费人成看片| 91精品国产国语对白视频| 久久久水蜜桃国产精品网| 亚洲国产欧美网| 亚洲精品美女久久av网站| 丰满少妇做爰视频| 一级片免费观看大全| 999精品在线视频| 日本vs欧美在线观看视频| 美女午夜性视频免费| 成人18禁高潮啪啪吃奶动态图| 中文字幕人妻丝袜一区二区| 亚洲av片天天在线观看| 亚洲精品中文字幕一二三四区 | 午夜免费鲁丝| 女人久久www免费人成看片| av电影中文网址| 最近最新免费中文字幕在线| 黄色 视频免费看| 亚洲国产av新网站| 男人添女人高潮全过程视频| 18禁裸乳无遮挡动漫免费视频| 老熟女久久久| 视频在线观看一区二区三区| 亚洲精品第二区| 18在线观看网站| 动漫黄色视频在线观看| 精品国产乱码久久久久久小说| 少妇 在线观看| 久久九九热精品免费| 亚洲三区欧美一区| 亚洲天堂av无毛| 免费日韩欧美在线观看| 亚洲国产成人一精品久久久| 色婷婷av一区二区三区视频| 欧美日韩av久久| 国产成人欧美| 成年人黄色毛片网站| 亚洲少妇的诱惑av| 黑人猛操日本美女一级片| 人妻一区二区av| 丝袜美足系列| 亚洲七黄色美女视频| 中文字幕人妻熟女乱码| 动漫黄色视频在线观看| 亚洲天堂av无毛| 自线自在国产av| 国产男女内射视频| 免费黄频网站在线观看国产| 黑人猛操日本美女一级片| 国产欧美日韩精品亚洲av| 夜夜骑夜夜射夜夜干| 亚洲人成电影观看| 999精品在线视频| 日韩制服骚丝袜av| 在线观看人妻少妇| 成人三级做爰电影| 亚洲男人天堂网一区| 免费一级毛片在线播放高清视频 | 少妇粗大呻吟视频| 亚洲色图 男人天堂 中文字幕| 午夜免费成人在线视频| 精品国产一区二区久久| 精品一区二区三卡| 成人18禁高潮啪啪吃奶动态图| 亚洲国产日韩一区二区| 一本久久精品| 亚洲人成电影观看| 脱女人内裤的视频| 岛国在线观看网站| 久久久久国内视频| 三上悠亚av全集在线观看| 久久这里只有精品19| 99国产极品粉嫩在线观看| 国产精品欧美亚洲77777| 黄色视频不卡| 免费日韩欧美在线观看| 如日韩欧美国产精品一区二区三区| 亚洲精品中文字幕在线视频|