• <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)研究
    人工降雪
    91大片在线观看| 亚洲自拍偷在线| 国产一区二区三区综合在线观看| 亚洲色图 男人天堂 中文字幕| 高清av免费在线| 中文字幕精品免费在线观看视频| 久久久国产欧美日韩av| 51午夜福利影视在线观看| 两性夫妻黄色片| 国产aⅴ精品一区二区三区波| 美女国产高潮福利片在线看| 亚洲美女黄片视频| 极品教师在线免费播放| 性欧美人与动物交配| 亚洲精品av麻豆狂野| 一二三四在线观看免费中文在| 久久国产精品男人的天堂亚洲| 亚洲性夜色夜夜综合| 精品熟女少妇八av免费久了| 男男h啪啪无遮挡| 亚洲一区二区三区色噜噜 | 满18在线观看网站| 亚洲国产中文字幕在线视频| 国产乱人伦免费视频| 天天影视国产精品| 精品福利永久在线观看| 色尼玛亚洲综合影院| 色老头精品视频在线观看| 久久国产乱子伦精品免费另类| 91成年电影在线观看| 国产欧美日韩综合在线一区二区| 热re99久久国产66热| ponron亚洲| 亚洲成av片中文字幕在线观看| 啪啪无遮挡十八禁网站| 天堂影院成人在线观看| 国产亚洲精品久久久久久毛片| 午夜a级毛片| 狠狠狠狠99中文字幕| 成人国语在线视频| 手机成人av网站| 国产又爽黄色视频| 日韩av在线大香蕉| 欧美中文日本在线观看视频| 美女国产高潮福利片在线看| 18禁国产床啪视频网站| 国产精品成人在线| 熟女少妇亚洲综合色aaa.| 久久青草综合色| 丰满饥渴人妻一区二区三| 久久精品国产清高在天天线| 男人操女人黄网站| 久久久久久免费高清国产稀缺| 国产精品香港三级国产av潘金莲| 亚洲精华国产精华精| 首页视频小说图片口味搜索| 国产成+人综合+亚洲专区| 欧美日本亚洲视频在线播放| 免费搜索国产男女视频| 成人18禁高潮啪啪吃奶动态图| 变态另类成人亚洲欧美熟女 | 日本a在线网址| 首页视频小说图片口味搜索| 水蜜桃什么品种好| 中文字幕人妻丝袜一区二区| 精品国产一区二区久久| 国产精品久久久久成人av| 精品久久久久久,| 精品国产乱码久久久久久男人| 操美女的视频在线观看| 91在线观看av| avwww免费| 一个人观看的视频www高清免费观看 | 亚洲成国产人片在线观看| 咕卡用的链子| 在线观看一区二区三区激情| 18禁美女被吸乳视频| 欧美中文日本在线观看视频| 国产精品国产高清国产av| 国产不卡一卡二| 亚洲精华国产精华精| 伦理电影免费视频| 91字幕亚洲| 亚洲av日韩精品久久久久久密| 国产成人av教育| 99久久99久久久精品蜜桃| 涩涩av久久男人的天堂| 最近最新免费中文字幕在线| xxx96com| 神马国产精品三级电影在线观看 | 精品少妇一区二区三区视频日本电影| 天堂动漫精品| e午夜精品久久久久久久| 午夜福利影视在线免费观看| 国产蜜桃级精品一区二区三区| 又大又爽又粗| 国产激情欧美一区二区| 99久久人妻综合| 叶爱在线成人免费视频播放| 成人18禁在线播放| 老司机靠b影院| 国产一区二区激情短视频| 免费在线观看亚洲国产| 欧美日韩av久久| 国产蜜桃级精品一区二区三区| 在线观看免费日韩欧美大片| 一级毛片女人18水好多| 丁香欧美五月| 青草久久国产| 在线观看舔阴道视频| 国产成人系列免费观看| av网站免费在线观看视频| 高清在线国产一区| 亚洲 国产 在线| 中文字幕人妻丝袜制服| 中文字幕最新亚洲高清| 日韩精品青青久久久久久| 99精国产麻豆久久婷婷| 欧美一级毛片孕妇| 亚洲一区二区三区色噜噜 | 午夜免费鲁丝| 99久久人妻综合| 狠狠狠狠99中文字幕| 成人18禁在线播放| 日韩欧美国产一区二区入口| 欧美日韩乱码在线| 午夜免费鲁丝| 男女午夜视频在线观看| 久久久久久久久久久久大奶| 久久国产亚洲av麻豆专区| 韩国精品一区二区三区| 久久亚洲真实| 亚洲精品久久午夜乱码| 国产xxxxx性猛交| 日韩人妻精品一区2区三区| 国产日韩一区二区三区精品不卡| 欧美日本亚洲视频在线播放| 大香蕉久久成人网| 亚洲av成人一区二区三| 国产精品国产av在线观看| 成人手机av| 真人做人爱边吃奶动态| 久久亚洲精品不卡| 18禁黄网站禁片午夜丰满| 国产免费男女视频| 99精品欧美一区二区三区四区| 国产国语露脸激情在线看| 美国免费a级毛片| 国产主播在线观看一区二区| 日本一区二区免费在线视频| 国产精品免费视频内射| 90打野战视频偷拍视频| 美女扒开内裤让男人捅视频| 99久久久亚洲精品蜜臀av| 黄色视频,在线免费观看| 亚洲国产中文字幕在线视频| 999精品在线视频| 国产av精品麻豆| 亚洲avbb在线观看| 久久久久久免费高清国产稀缺| 久久久久久大精品| 色尼玛亚洲综合影院| 每晚都被弄得嗷嗷叫到高潮| 国产成人精品久久二区二区91| 欧美久久黑人一区二区| 免费少妇av软件| 久久精品亚洲av国产电影网| 757午夜福利合集在线观看| 91麻豆精品激情在线观看国产 | 久久狼人影院| 69精品国产乱码久久久| 九色亚洲精品在线播放| 精品久久久久久久毛片微露脸| 色综合欧美亚洲国产小说| 久久香蕉精品热| 国产亚洲欧美在线一区二区| svipshipincom国产片| 每晚都被弄得嗷嗷叫到高潮| 一级片免费观看大全| 高潮久久久久久久久久久不卡| 人人妻,人人澡人人爽秒播| 一进一出好大好爽视频| 精品一区二区三区四区五区乱码| 精品少妇一区二区三区视频日本电影| 少妇裸体淫交视频免费看高清 | 国产真人三级小视频在线观看| 这个男人来自地球电影免费观看| av电影中文网址| 在线永久观看黄色视频| 欧美日韩福利视频一区二区| 国产精品国产高清国产av| 亚洲熟妇中文字幕五十中出 | 黄色女人牲交| 视频区图区小说| 欧美日韩瑟瑟在线播放| 琪琪午夜伦伦电影理论片6080| 老司机午夜福利在线观看视频| 成人特级黄色片久久久久久久| 久久久久久久午夜电影 | 老汉色av国产亚洲站长工具| 免费久久久久久久精品成人欧美视频| 波多野结衣高清无吗| av超薄肉色丝袜交足视频| 脱女人内裤的视频| 亚洲激情在线av| 亚洲人成伊人成综合网2020| 亚洲 欧美 日韩 在线 免费| 亚洲一码二码三码区别大吗| 在线观看免费高清a一片| 一a级毛片在线观看| 后天国语完整版免费观看| 亚洲专区字幕在线| 亚洲精品粉嫩美女一区| 国产人伦9x9x在线观看| 欧美在线黄色| 日韩有码中文字幕| 国产精品亚洲一级av第二区| 51午夜福利影视在线观看| xxxhd国产人妻xxx| 18禁国产床啪视频网站| 亚洲国产中文字幕在线视频| 波多野结衣av一区二区av| 十八禁人妻一区二区| 女人爽到高潮嗷嗷叫在线视频| 久久狼人影院| 首页视频小说图片口味搜索| 亚洲国产欧美日韩在线播放| www.www免费av| 欧美日韩国产mv在线观看视频| 1024香蕉在线观看| 精品日产1卡2卡| 99久久久亚洲精品蜜臀av| 日韩av在线大香蕉| 国产免费男女视频| 久久久精品欧美日韩精品| 国产高清国产精品国产三级| 在线观看一区二区三区| 首页视频小说图片口味搜索| 三上悠亚av全集在线观看| 免费搜索国产男女视频| 欧美性长视频在线观看| 国产成人精品在线电影| xxxhd国产人妻xxx| 日韩欧美一区二区三区在线观看| 久久久久久亚洲精品国产蜜桃av| 三上悠亚av全集在线观看| 久久精品91无色码中文字幕| 国产精品国产高清国产av| 久久精品人人爽人人爽视色| 夜夜躁狠狠躁天天躁| 美女大奶头视频| 一级作爱视频免费观看| 激情视频va一区二区三区| 久久人妻熟女aⅴ| 脱女人内裤的视频| 亚洲男人天堂网一区| 婷婷丁香在线五月| 韩国精品一区二区三区| 在线观看免费午夜福利视频| 动漫黄色视频在线观看| 国产精品久久视频播放| 又黄又爽又免费观看的视频| 中文字幕人妻丝袜制服| 91九色精品人成在线观看| 亚洲久久久国产精品| 法律面前人人平等表现在哪些方面| 少妇粗大呻吟视频| 男人的好看免费观看在线视频 | 成人黄色视频免费在线看| 国产乱人伦免费视频| 嫩草影院精品99| 精品一区二区三区av网在线观看| 免费av毛片视频| 性欧美人与动物交配| 三上悠亚av全集在线观看| 超色免费av| 色老头精品视频在线观看| 午夜福利在线免费观看网站| 精品高清国产在线一区| 亚洲熟妇熟女久久| 淫妇啪啪啪对白视频| 亚洲国产看品久久| 欧美成人免费av一区二区三区| 999久久久国产精品视频| 嫁个100分男人电影在线观看| 99国产精品一区二区三区| 国产亚洲欧美在线一区二区| 啦啦啦在线免费观看视频4| 国产精品久久久人人做人人爽| 国产精品亚洲av一区麻豆| 女人被躁到高潮嗷嗷叫费观| 欧美成狂野欧美在线观看| 九色亚洲精品在线播放| 久久亚洲真实| 欧美不卡视频在线免费观看 | 一级a爱视频在线免费观看| 国产有黄有色有爽视频| 在线观看舔阴道视频| 亚洲精品国产精品久久久不卡| 成人影院久久| 精品人妻在线不人妻| 首页视频小说图片口味搜索| 久久人人精品亚洲av| 99久久国产精品久久久| 伊人久久大香线蕉亚洲五| 成人手机av| 又紧又爽又黄一区二区| 国产精品 国内视频| 一本大道久久a久久精品| 精品久久久久久久毛片微露脸| 脱女人内裤的视频| 老汉色av国产亚洲站长工具| 脱女人内裤的视频| 欧美激情极品国产一区二区三区| 亚洲精品久久午夜乱码| 在线天堂中文资源库| 成人18禁高潮啪啪吃奶动态图| 日本 av在线| 欧美成人免费av一区二区三区| 国产三级在线视频| 十分钟在线观看高清视频www| 午夜精品国产一区二区电影| 久久精品亚洲av国产电影网| 精品国产超薄肉色丝袜足j| 日韩高清综合在线| 成年版毛片免费区| 一级片免费观看大全| 色婷婷久久久亚洲欧美| 男女下面插进去视频免费观看| 高清欧美精品videossex| 日韩视频一区二区在线观看| 91九色精品人成在线观看| 久久久久久大精品| 欧美乱码精品一区二区三区| 国产有黄有色有爽视频| 欧美 亚洲 国产 日韩一| 嫩草影院精品99| 国产成人精品久久二区二区免费| 久久久水蜜桃国产精品网| 人人妻人人添人人爽欧美一区卜| 亚洲第一青青草原| 人人澡人人妻人| 少妇被粗大的猛进出69影院| 国产片内射在线| 久热这里只有精品99| 高清av免费在线| 国产xxxxx性猛交| 精品久久蜜臀av无| 在线观看一区二区三区| 黑人猛操日本美女一级片| 欧美黑人精品巨大| 天堂影院成人在线观看| 国产伦人伦偷精品视频| 亚洲专区字幕在线| 满18在线观看网站| 水蜜桃什么品种好| a在线观看视频网站| 麻豆成人av在线观看| 母亲3免费完整高清在线观看| 老熟妇乱子伦视频在线观看| 最近最新免费中文字幕在线| 亚洲成人免费电影在线观看| 91国产中文字幕| 夜夜看夜夜爽夜夜摸 | 欧美精品啪啪一区二区三区| 免费搜索国产男女视频| 亚洲精品成人av观看孕妇| 亚洲国产欧美一区二区综合| 亚洲精品一卡2卡三卡4卡5卡| 人成视频在线观看免费观看| 精品少妇一区二区三区视频日本电影| www.www免费av| 亚洲欧美激情在线| 欧美成狂野欧美在线观看| 999久久久精品免费观看国产| 日韩大码丰满熟妇| 在线观看www视频免费| 如日韩欧美国产精品一区二区三区| 一区二区三区激情视频| 久9热在线精品视频| 悠悠久久av| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩精品青青久久久久久| 国产麻豆69| 免费看a级黄色片| 欧美黄色淫秽网站| 美女午夜性视频免费| 亚洲成国产人片在线观看| 丁香欧美五月| 国产精品永久免费网站| 日韩 欧美 亚洲 中文字幕| 亚洲精品国产一区二区精华液| 一区福利在线观看| 国产高清videossex| 另类亚洲欧美激情| 精品国产乱子伦一区二区三区| 国产精品一区二区三区四区久久 | 18禁黄网站禁片午夜丰满| 97碰自拍视频| 十分钟在线观看高清视频www| 亚洲国产精品999在线| 成人三级做爰电影| 国产蜜桃级精品一区二区三区| 国产高清videossex| 午夜福利,免费看| 午夜免费鲁丝| 午夜福利在线免费观看网站| 久久精品亚洲熟妇少妇任你| 一级片免费观看大全| 女同久久另类99精品国产91| 亚洲色图 男人天堂 中文字幕| 久久久久久人人人人人| e午夜精品久久久久久久| 国产成人系列免费观看| 国产精品永久免费网站| 一二三四社区在线视频社区8| 日韩大尺度精品在线看网址 | 很黄的视频免费| 免费看a级黄色片| 99久久综合精品五月天人人| 免费av毛片视频| 人人澡人人妻人| 久久香蕉国产精品| av电影中文网址| 亚洲av熟女| 又紧又爽又黄一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品一区av在线观看| 久久精品国产99精品国产亚洲性色 | 亚洲午夜精品一区,二区,三区| 啪啪无遮挡十八禁网站| 国产99白浆流出| 国产极品粉嫩免费观看在线| 久久久国产欧美日韩av| 国内久久婷婷六月综合欲色啪| 亚洲精品一二三| 午夜免费鲁丝| 久久中文字幕人妻熟女| 色婷婷久久久亚洲欧美| a在线观看视频网站| 免费在线观看影片大全网站| 亚洲免费av在线视频| 欧美在线黄色| 国产高清videossex| 国产精品美女特级片免费视频播放器 | 老熟妇仑乱视频hdxx| 日韩一卡2卡3卡4卡2021年| 国产在线观看jvid| 久久久久精品国产欧美久久久| 亚洲av美国av| 精品国产国语对白av| 最近最新中文字幕大全免费视频| 亚洲中文字幕日韩| 免费在线观看黄色视频的| 国产男靠女视频免费网站| 老司机午夜福利在线观看视频| av视频免费观看在线观看| 黄色怎么调成土黄色| 嫩草影院精品99| 美女国产高潮福利片在线看| 香蕉国产在线看| 久久久久亚洲av毛片大全| 精品日产1卡2卡| a级片在线免费高清观看视频| 欧美最黄视频在线播放免费 | 看片在线看免费视频| 黄色视频,在线免费观看| 老熟妇乱子伦视频在线观看| 欧美日韩视频精品一区| 色在线成人网| 最新美女视频免费是黄的| 久久久久久久午夜电影 | 欧美乱妇无乱码| 日韩欧美国产一区二区入口| 欧美另类亚洲清纯唯美| 久久精品影院6| 精品久久久精品久久久| 久久久久国内视频| 日韩精品青青久久久久久| 国产亚洲欧美精品永久| 天堂动漫精品| 亚洲自偷自拍图片 自拍| 中国美女看黄片| 日本免费一区二区三区高清不卡 | 欧美国产精品va在线观看不卡| 一级,二级,三级黄色视频| 成人av一区二区三区在线看| 国产单亲对白刺激| 88av欧美| 无限看片的www在线观看| 久久久久国内视频| 亚洲精品在线观看二区| www.熟女人妻精品国产| videosex国产| 亚洲一区中文字幕在线| 如日韩欧美国产精品一区二区三区| 精品人妻1区二区| 免费av中文字幕在线| 成年人免费黄色播放视频| 超碰成人久久| 午夜精品在线福利| 亚洲国产精品sss在线观看 | 在线观看免费视频网站a站| 怎么达到女性高潮| 亚洲精品中文字幕一二三四区| 在线观看www视频免费| 亚洲成a人片在线一区二区| 在线永久观看黄色视频| 国产成人精品久久二区二区免费| 啪啪无遮挡十八禁网站| 18禁观看日本| 午夜成年电影在线免费观看| 19禁男女啪啪无遮挡网站| 欧美 亚洲 国产 日韩一| 久久性视频一级片| 日韩成人在线观看一区二区三区| 91大片在线观看| 黄色视频,在线免费观看| 侵犯人妻中文字幕一二三四区| 99久久国产精品久久久| av视频免费观看在线观看| 美女福利国产在线| 一级,二级,三级黄色视频| 好男人电影高清在线观看| 老司机福利观看| 国产精品二区激情视频| 欧美性长视频在线观看| 国产欧美日韩一区二区精品| 精品国产一区二区三区四区第35| 一级毛片女人18水好多| 国产成人精品久久二区二区91| 亚洲国产欧美日韩在线播放| 国产亚洲欧美精品永久| 国产三级黄色录像| 欧美不卡视频在线免费观看 | 男女下面进入的视频免费午夜 | 午夜老司机福利片| 午夜福利一区二区在线看| 免费一级毛片在线播放高清视频 | 国产精品永久免费网站| 亚洲av电影在线进入| 久热爱精品视频在线9| 精品国产一区二区久久| 午夜免费观看网址| 成人永久免费在线观看视频| 国产色视频综合| 在线观看一区二区三区激情| 免费观看精品视频网站| 狠狠狠狠99中文字幕| 变态另类成人亚洲欧美熟女 | 国产熟女xx| 桃红色精品国产亚洲av| 亚洲国产毛片av蜜桃av| 淫妇啪啪啪对白视频| а√天堂www在线а√下载| 12—13女人毛片做爰片一| av网站在线播放免费| 五月开心婷婷网| 一二三四在线观看免费中文在| 免费在线观看亚洲国产| 免费在线观看完整版高清| 两性午夜刺激爽爽歪歪视频在线观看 | 久久国产亚洲av麻豆专区| 夜夜爽天天搞| 日韩精品青青久久久久久| av天堂久久9| 精品国产乱码久久久久久男人| 黑人巨大精品欧美一区二区mp4| e午夜精品久久久久久久| 99精品欧美一区二区三区四区| 欧美国产精品va在线观看不卡| 久久香蕉激情| 午夜a级毛片| 精品卡一卡二卡四卡免费| 国产精品偷伦视频观看了| 国产又爽黄色视频| 一级毛片高清免费大全| 高清毛片免费观看视频网站 | 三级毛片av免费| 欧美老熟妇乱子伦牲交| 国产精品乱码一区二三区的特点 | 色播在线永久视频| 在线观看免费日韩欧美大片| 啦啦啦免费观看视频1| 啦啦啦在线免费观看视频4| 午夜福利在线观看吧| 国产男靠女视频免费网站| 50天的宝宝边吃奶边哭怎么回事| 亚洲七黄色美女视频| 纯流量卡能插随身wifi吗| 国产精品 欧美亚洲| 亚洲九九香蕉| 色哟哟哟哟哟哟| 国产精品 欧美亚洲| 久久中文字幕人妻熟女| 亚洲成人久久性| 亚洲精品国产色婷婷电影| 香蕉丝袜av| 老司机靠b影院| 夫妻午夜视频| av在线天堂中文字幕 | 中亚洲国语对白在线视频| 村上凉子中文字幕在线| xxxhd国产人妻xxx| 1024视频免费在线观看| 久久久国产欧美日韩av| 黄色毛片三级朝国网站| 国产成人影院久久av| 亚洲国产精品sss在线观看 | 十八禁网站免费在线| 国产精品免费视频内射| 国产单亲对白刺激| 亚洲av日韩精品久久久久久密| 亚洲精品在线美女| www.精华液|