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

    冰晶異質核化對雷暴云電過程影響的數(shù)值模擬

    2021-09-22 01:39:40李璐瀅師正譚涌波汪海潮林曉彤
    熱帶氣象學報 2021年3期
    關鍵詞:核化冰晶雷暴

    李璐瀅,師正,譚涌波,汪海潮,林曉彤

    (南京信息工程大學氣象災害教育部重點實驗室/氣候與環(huán)境變化國際合作聯(lián)合實驗室/氣象災害預報預警與評估協(xié)同創(chuàng)新中心/中國氣象局氣溶膠與云降水重點開放實驗室,江蘇 南京210044)

    1 引言

    大量的觀測與試驗研究證明,雷暴云起電的主要機制是非感應起電機制,即當過冷水存在,大小冰相粒子碰撞而發(fā)生的電荷轉移機制[1-6],因此雷暴云的冰晶在起、放電過程中有著重要的影響。Lal等[7]認為冰晶直徑的增加有利于閃電的發(fā)生,但粒徑超過23~25μm的冰晶反而抑制了放電。Zheng等[8]通過日本冬季雷暴的觀測發(fā)現(xiàn),伴隨著弱雷達回波的反極性電荷結構主要是由零度層以上雪和冰晶之間的起電過程引起的。Ouyang等[9]通過WRF模式探討了霰粒密度對電荷結構的影響,結果表明隨著霰密度的增加,主正電荷區(qū)減弱,負的非感應起電率降低。冰相粒子的垂直分布也是通過直接影響非感應起電率,從而導致了雷暴云內電荷結構、放電特征及閃電頻次出現(xiàn)顯著變化[10-15]。實際上,冰相粒子的分布特征與核化過程息息相關,準確地描述冰晶核化方式更有利于提高大氣放電與云內動力微物理過程聯(lián)系的認知。

    在冰云的相關研究中發(fā)現(xiàn),冰晶粒子的形成途徑主要是同質核化與異質核化[16-20]。目前,對同質核化過程的研究認為,在環(huán)境溫度低于-37°C且相對濕度較高的條件下,可溶性氣溶膠作為云凝結核,參與過冷水的凍結形成冰晶[21]。而異質核化機制相對較復雜,如凝結凍結核化、接觸核化、浸潤核化、沉積核化等,多種氣溶膠可活化為冰核,如黑碳、礦物粉塵、細菌等[16,22-25],所以仍然存在較多的爭議。例如,凝結凍結和浸潤核化之間的界限似乎并不明確[26-27];沉積核化可能是由于粒子表面孔隙中水的凝結而形成冰晶,而這應屬于凝結凍結核化機制[28];黑碳難以通過浸潤核化對冰晶的生成造成影響[29]。

    近年來,隨著對冰核的觀測試驗和數(shù)值模式工作的展開,冰核參數(shù)化方案也趨于全面和復雜。例如早期研究認為冰核主要受到溫濕條件的限制,往往采用活化溫度或飽和度對冰核濃度進行計算[30],而后的研究發(fā)現(xiàn)不同類型的氣溶膠作為冰核時具有不同的活性[31],例如Phillips等[32]認為氣溶膠組分及氣溶膠表面積應考慮進核化方案中。Ervens等[33]的方案同樣肯定了冰核粒子表面特性的影響。DeMott等[34]發(fā)現(xiàn)大于0.5μm的氣溶膠更具有轉化為冰核的活性,并提出了與氣溶膠有關的計算冰核參數(shù)化公式。這些方案主要為應用于大尺度天氣和氣候數(shù)值模式開發(fā),適用于云模式的方案較少,且大部分采用參數(shù)化的方法,冰晶的產生或是只與溫度和過飽和度有關,或是未鏈接氣溶膠濃度、具體化學組分,模擬結果仍然與真實冰晶分布存在一定差距。此外,大部分研究僅限于核化過程對降水、云物理屬性等方面的影響,少數(shù)涉及雷暴云電過程的研究也尚未從冰晶具體核化途徑出發(fā),詳細探討氣溶膠與雷暴云間具體云物理過程的聯(lián)系。因此,為了更全面地了解氣溶膠具體作用于雷暴云的內在機理,本研究對原有冰晶核化方案進行改良,更換為Diehl等[35]異質核化方案,用于研究雷暴云的起電過程。該方案基于較新實驗結果和云氣塊模型,包含了典型的三類異質核化過程,同時考慮了冰核組分、各類冰核活化率等因素,能更真實地反映冰晶分布特征。利用該核化方案,本研究在已有的二維雷暴云起、放電模式中,通過對比不同異質核化途徑下云內冰相粒子分布的差異,探討不同異質核化過程影響雷暴云電過程的具體物理機制。

    2 數(shù)值模式介紹

    本研究采用譚涌波等[13,36-38]改進的二維積云模式,其中包含了雷暴云起、放電參數(shù)化方案。積云模式以云滴、雨滴、冰晶、霰、雹的混合比與數(shù)濃度為預報量,計算了云物理特征量、凝結蒸發(fā)過程、凍結融化過程等。耦合的雷暴云起電參數(shù)化方案中,主要考慮了非感應起電與感應起電過程。其中非感應起電采用的是Gardiner-Pereyra方案,描述了冰晶與霰粒子的碰撞導致電荷轉移[39];感應起電則采用Ziegler等[40]參數(shù)化方案,描述了霰、雹分別與云滴碰撞分離導致的電荷轉移。

    為深入分析冰晶異質核化過程,研究中選取Diehl等[35]的異質核化方案,涉及接觸核化,浸潤核化和沉積核化,用于取代原Fletcher[30]實驗室測得的冰核活化經驗公式,同時耦合氣溶膠作為冰核模塊,討論雷暴云微物理過程與電過程對不同冰晶異質核化過程的響應。下面對所采用的新型冰晶核化方案進行簡要描述。

    2.1 異質核化方案的改進

    2.1.1 接觸核化方案

    接觸凍結核化通常是指干燥的氣溶膠顆粒與液滴碰并凍結的過程。因此方案中重點考慮以下三點:(1)環(huán)境溫度低于凍結溫度;(2)冰核與過冷液滴相互碰并;(3)作為冰核的氣溶膠粒子為干燥狀態(tài)。不同類型的冰核對應不同的凍結臨界溫度,詳見文獻[35],此處不再贅述。當環(huán)境溫度低于閾值時,假設過冷液滴與冰核碰撞,發(fā)生接觸凍結核化:

    其中,Nc代表凍結生成冰晶的數(shù)濃度,F(xiàn)INc是氣溶膠活化為可接觸凍結的冰核的百分數(shù),Nliq為可與冰核碰撞的液滴數(shù)濃度,T為云內溫度,單位為℃,ac和bc為無量綱參數(shù),具體數(shù)值參考表1。需要注意的是,該方案僅在1≤(ac T+bc)≤0條件下成立。冰核組分設置為:32%的污染物顆粒、24%的長石、14%的伊利石、12%的高嶺石。

    表1 接觸核化參數(shù)公式

    2.1.2 凝結凍結和浸潤核化方案

    研究表明,浸潤凍結溫度與水成物粒子體積呈反比[41]。根據(jù)Diehl等[35]的描述,當達到啟動溫度時,浸潤凍結核化率為:

    其中,Ni代表在溫度T下生成的冰晶數(shù)濃度,mice是浸入液滴中冰核的質量,F(xiàn)INi是氣溶膠活化為可浸潤凍結的冰核的百分數(shù),K(T)為冰核溫度譜:

    這里,nact是在溫度T下單位質量活化冰核的數(shù)濃度,ai和bi為無量綱參數(shù),具體數(shù)值參考表2。當環(huán)境溫度降至最低溫度以下,核化率限制為最低溫度下的固定值。該方案中冰核組分與接觸核化方案一致。

    表2 浸潤核化參數(shù)公式

    凝結凍結過程與浸潤凍結類似,都是通過冰核浸入過冷液滴觸發(fā)凍結[42],所以實際上方程(2)和(3)已經隱含了凝結凍結過程。

    2.1.3 沉積核化方案

    Diehl等[23]的研究表明,冰面過飽和度的增加可有效提高冰核活化率。在方案中,沉積核化發(fā)生的條件為:(1)存在不參與接觸凍結過程的活化冰核;(2)環(huán)境溫度與冰面過飽和度達到閾值。沉積核化的計算公式:

    其中,Nd代表沉積核化生成的冰晶數(shù)濃度,Naerd是氣溶膠活化為可沉積凍結的冰核數(shù)濃度,sice為冰過飽和度,ad和b d為無量綱參數(shù),具體數(shù)值參考表3。由于沉積凍結方案中沒有高嶺石和污染物顆粒的參數(shù)化方程,所以將冰核組分設置為:32%的沙塵、24%的長石、14%的伊利石、12%的生物氣溶膠。

    表3 沉積核化參數(shù)公式

    2.2 初始背景場簡介

    研究選擇山地弱雷暴作為背景,其初始探測場的溫濕層結曲線和垂直風廓線見圖1(見下頁)。濕熱泡位于模擬域中心高度1 km的格點上,溫度擾動為3.5 K,濕度擾動為60%。模擬域范圍為76 km×20 km,分辨率為250 m×250 m。模擬時間為80 min,時間步長為2 s,地面氣溶膠背景選取污染的大氣環(huán)境,初始數(shù)濃度設置為1000 cm-3。研究中共設置三組模擬實驗,接觸方案僅包含接觸核化過程,浸潤方案包含了凝結凍結和浸潤核化過程,沉積方案僅包含沉積核化過程。此外,模式中還考慮了同質核化過程,并設置同一初始云凝結核濃度,因此三種方案的同質凍結過程一致,下文不再討論。

    圖1 試驗個例的環(huán)境溫濕層結(a)和垂直風廓線(b)

    3 模式結果與分析

    3.1 冰相粒子空間分布的變化

    不同核化過程控制初始冰晶的生成,從而在冰晶的分布高度與溫度區(qū)間顯示出較大的差異。圖2給出了冰晶數(shù)濃度和比含水量隨時間的變化。通過對比可發(fā)現(xiàn),接觸核化生成的冰晶數(shù)濃度和比含水量是最小的,主要集中在相對較高的溫度范圍(-20.0~-13.8°C),且持續(xù)生成的時間較短,只在24~37 min期間數(shù)濃度能夠達到1×105kg-1,40 min后冰晶主要通過同質核化產生。這主要是由于初期液態(tài)水主要分布在高溫區(qū),只有小部分氣溶膠能活化并轉化為冰晶。之后的降水、淞附等過程不斷消耗云中液滴,由于冰核粒徑較小,即使借助上升氣流,也難以與其他大的水成物粒子競爭云水,這些因素進一步降低了發(fā)生接觸核化的概率。浸潤核化方案中,相對較大的冰核活化溫度區(qū)間有效地促進了冰晶的發(fā)展,冰晶顯現(xiàn)出較廣泛的垂直分布。隨著冰晶大量的產生,貝吉隆過程不斷加強,促進了液滴的蒸發(fā),使冰晶的質量得到進一步的增長。所以在最初的發(fā)展中,冰晶的比含水量就超過3 g/kg。當雷暴云進入成熟階段,冰核活化釋放大量的潛熱,旺盛的上升氣流將冰晶向云頂-40°C以上的溫度層輸送,數(shù)濃度最高達到了6.87×106kg-1。對于沉積核化方案,由于受到冰面過飽和度和活化溫度的雙重限制,個例中直到25 min才有冰晶產生。盡管冰晶的最大比含水量不到3 g/kg,其數(shù)濃度仍然能達到106kg-1,這證明沉積核化主要作用在低溫區(qū)(-20℃以上),是中高層小冰晶的重要生成途徑之一。三個方案下冰晶分布的對比足以表明,浸潤核化是最重要的異質核化方式。這與之前的研究結果相吻合[42-45]。

    圖2 冰晶比含水量與數(shù)濃度的時空分布

    為了進一步探討核化過程分別對冰相粒子的貢獻,分別統(tǒng)計了三個方案下冰晶異質核化率(圖3),接觸凍結核化效率最低,該核化過程產生的冰晶主要分布在雷暴云的中下層3~5 km,4.8 km處核化率最高,同時也是浸潤核化最活躍的高度。不過,浸潤核的活化率遠高于接觸核,浸潤核化生成的冰晶質量比接觸核化高出兩個量級,且凍結高度從4 km延伸至9 km。而沉積核化的凍結高度更高,在7.4 km處達到最大。也就是說,-20~-40°C是異質核化效率最高的溫度區(qū)間,過低的溫度環(huán)境反而不利于異質核化的發(fā)展。

    圖3 冰晶異質核化率的垂直分布

    霰粒子的生成依賴于冰晶的轉化,而云水則是霰粒子增長的重要來源之一,所以不同核化方案中霰的分布取決于冰晶和云水條件,模擬結果如圖4所示。所有個例的霰粒最初都在-20°C附近產生,其中沉積核化個例生成的霰粒比含水量和數(shù)濃度最大,發(fā)展最旺盛。在溫度相對較高的過冷區(qū),霰粒子由于比含水量較大(最大比含水量為26.2 g/kg)而迅速降落。當下沉至2 km高度仍有部分質量較小的霰粒(約為5×104kg-1)流向云頂,與同質核化共同貢獻了云頂?shù)啮钡纳L。接觸核化方案中的冰晶轉化量較小(表4,為5.65×10-3g/kg),僅為沉積方案(1.71×10-2g/kg)的0.33倍,不過接觸核化過程主要發(fā)生在液態(tài)水含量相對較高的過冷區(qū),霰粒與過冷液滴的碰撞效率反而提高了,所以生成的霰略少于沉積方案,最大比含水量也能達到7.97 g/kg。而浸潤核化中霰粒主要集中在0~-20°C,最大比含水量僅為7.58 g/kg,是三個方案里分布范圍最小的。這主要是由于在浸潤核化過程消耗了大量的云水,同時削弱了水汽的向上輸入,抑制了高層霰的增長。所以,盡管接觸核化能通過提前釋放潛熱,增加對流云過冷區(qū)的不穩(wěn)定能量,使上升氣流得到加強,也難以形成利于霰粒增長的環(huán)境。

    表4 霰粒源匯項

    圖4 霰比含水量與數(shù)濃度的時空分布

    綜上所述,浸潤核的活化溫度范圍較廣,有利于冰晶的大量產生,該過程大量消耗云水,抑制了隨后霰的增長。沉積核化主要貢獻了低溫區(qū)冰晶,而云中層仍存在豐富水汽與液水,顯著提高了霰收集云滴的效率,表現(xiàn)為極高的霰比含水量。接觸核活化率為三者最低,盡管該個例下冰晶的生成量最少,但云水環(huán)境足以滿足霰粒增長的條件。也就是說,接觸核化以最少的冰晶生成量顯著地改變了云中水成物粒子的分布。這些結論表明,不同異質核化方式在雷暴云的電過程中可能產生相當大的影響。

    3.2 雷暴云起電特征與電荷結構

    冰晶與霰作為主要的荷電的粒子,由不同異質核化過程導致的冰相粒子微物理特征變化可能顯著影響起電過程。圖5為冰晶最大非感應起電率隨時間的變化。接觸與浸潤核化方案的非感應起電于22 min左右開始發(fā)生,而沉積方案由于冰晶與霰的生成稍晚,直至30 min才開始出現(xiàn)電荷分離。此外,由于所有個例中的冰晶粒子主要分布在-20°C以上區(qū)域,所以最大非感應起電率位于-20~-40°C,與冰晶分布的中心區(qū)域重合。其中沉積方案低溫區(qū)的非感應起電率最高,達到了1401.80 pc/(m3·s)。而浸潤方案中,可與冰晶碰撞的霰主要分布在0~-20°C溫區(qū),所以在高于反轉溫度(-15°C)區(qū)域的冰晶所攜帶負電荷量為三者最高。接觸核化生成的冰晶較少,非感應起電率主要依賴同質核化形成的冰晶,所以該方案下的起電并不強。

    圖5 最大非感應起電率隨時間的變化

    雷暴云的電荷結構主要是冰相粒子在空間里相互碰撞的結果,所以三種方案的電荷結構均呈現(xiàn)出不同的特點(圖6)??傮w上來看,電荷結構均由較復雜的三極性發(fā)展為偶極性,但隨著冰晶分布高度的增加,次正電荷區(qū)出現(xiàn)不同程度的減弱,主正電荷區(qū)明顯抬高。在云內電荷分離的初始階段,接觸方案與浸潤方案的電荷結構呈反偶極性,主電荷區(qū)主要位于3~6 km高度,隨后,在次正電荷區(qū)下出現(xiàn)了小范圍的負電荷團(約-0.5 nC/(m2·s))。對比該時刻的冰晶分布(圖2a、2b)發(fā)現(xiàn),這與異核化冰晶的形成高度相吻合,該高度層的環(huán)境溫度高于反轉溫度,使異質核化生成的冰晶顯負電。較于接觸方案,浸潤方案形成的冰晶數(shù)濃度更高,次正電荷區(qū)更弱,小范圍的負電荷團更強。隨著雷暴進入旺盛期,冰晶大量生成釋放潛熱,云體溫度進一步降低,在上升氣流的幫助下,三個方案下的主正電荷區(qū)均突破了-40°C溫度層,其中浸潤方案與沉積方案的主正電荷區(qū)分別在38 min和42 min達到9.5 km高度。對于浸潤方案,其冰晶分布范圍較廣,數(shù)濃度相對較高,而霰粒集中在5 km以下,導致云體中部本是荷負電的霰粒應聚集的區(qū)域,出現(xiàn)了超過0.5 nC/(m2·s)的正電荷團(冰晶)。而接觸方案與沉積方案中,底部的次正電荷區(qū)在該階段就已經消散,電荷結構提前演變?yōu)榕紭O性。雖然兩個方案的電荷結構發(fā)展類似,但成因完全相反。接觸方案的對流相對較弱,霰粒更早地降落至反轉溫度層高度以下,形成次正電荷區(qū),隨著霰粒降落至地面,次正電荷區(qū)消失。反觀沉積方案,冰晶主要在-20°C形成,旺盛的上升氣流不斷地將冰晶向云頂輸送,所以只有位于-20°C附近的霰粒能與冰晶碰撞起電(圖5c),直接導致了次正電荷區(qū)范圍的縮小,電荷量的減弱。此外,霰的下落速度本就比冰晶更快,中層充沛的云水促使霰??焖僭鲩L,進一步增加了霰的下落速度,大幅度降低了霰與冰晶碰撞的概率,從而造成了次正電荷區(qū)的消失。也就是說,不同異質核化生成的冰晶對空間電荷結構的影響主要體現(xiàn)在次正電荷區(qū)上。接觸方案生成的次正電荷區(qū)持續(xù)時間更長,分布更廣,浸潤方案次之,沉積方案最小。這與王夢旖等[15]對大氣冰核譜在雷暴云電過程的研究結果類似。

    圖6 三種方案下空間電荷結構分布

    4 結論

    研究通過改變異質核化方案進行敏感性試驗,重點分析了三種異質核化過程對冰相粒子的發(fā)展在雷暴云起電過程中的影響。

    (1)不同異質核化過程在雷暴云內有效高度存在差異。接觸核化生成的冰晶量最少,僅對雷暴云中下層3~5 km處的冰晶有貢獻,所以霰粒子數(shù)濃度并不高,但優(yōu)良的云水條件足以使霰粒充分增長;浸潤核化是最重要的異質核化方式,冰晶分布更廣泛,4.8 km處核化率最高,但大量的冰晶消耗了云水,抑制了霰的生長,導致較小的比含水量;沉積核化主要影響了云砧處的冰晶,而云中層仍存在豐富水汽與液水,顯著提高了霰收集云滴的效率,表現(xiàn)為極高的霰比含水量。

    (2)冰相粒子的分布差異直接影響了雷暴云非感應起電過程,從而表現(xiàn)出不同的電荷結構。總體上來看,電荷結構均由較復雜的三極性發(fā)展為偶極性。隨著冰相粒子的分布在接觸方案、浸潤方案和沉積方案里依次升高,最大非感應起電率也依次增加,主正電荷區(qū)依次升高,但較高與較低的冰晶分布高度都不利于維持次正電荷區(qū),反而使其更快地消散。

    本研究通過加入具體的異質核化過程,提高了模式的模擬能力,進一步揭示了氣溶膠作為不同冰核影響雷暴云微物理、起電過程的物理機制。由于文章篇幅的限制,本文只給出一類污染氣溶膠作為具有代表性的模擬結果,所以在未來的工作中,我們將進一步探討其它環(huán)境背景下不同類型冰核對雷暴云放電行為的影響。

    致謝:感謝中國氣象科學研究院胡志晉研究員提供了本文所采用的積雨云微物理框架。

    猜你喜歡
    核化冰晶雷暴
    冰晶世界,凝凍微觀的美麗
    手持式核化探測儀器發(fā)展現(xiàn)狀與應用展望
    新德里雷暴
    雪為什么只在冬天下
    人工影響天氣碘化銀催化劑研究進展
    阜新地區(qū)雷暴活動特點研究
    綠色科技(2018年22期)2019-01-15 05:36:08
    小泥人冰晶畫 蘊藏大商機
    廣西富川縣雷暴氣候特征分析
    冰晶奇域
    琴童(2016年5期)2016-05-14 15:43:58
    核化生醫(yī)學救援程序與訓練思考
    91字幕亚洲| 久久久久久国产a免费观看| 亚洲黑人精品在线| 成人亚洲精品av一区二区| 窝窝影院91人妻| 精品不卡国产一区二区三区| 国内揄拍国产精品人妻在线| 欧美zozozo另类| 国产视频内射| 久久精品亚洲精品国产色婷小说| 天天躁狠狠躁夜夜躁狠狠躁| 黄频高清免费视频| 性欧美人与动物交配| 欧美中文综合在线视频| 高清在线国产一区| 国产成人欧美在线观看| 无人区码免费观看不卡| 亚洲av中文字字幕乱码综合| 日本免费一区二区三区高清不卡| 国产午夜精品久久久久久| 美女cb高潮喷水在线观看 | 一进一出好大好爽视频| 丝袜人妻中文字幕| xxx96com| 999精品在线视频| 欧美日韩黄片免| 亚洲专区中文字幕在线| 美女午夜性视频免费| 色综合亚洲欧美另类图片| 国产精品,欧美在线| 我要搜黄色片| 国产精品 国内视频| 久99久视频精品免费| 亚洲一区高清亚洲精品| 校园春色视频在线观看| 亚洲第一电影网av| 欧美高清成人免费视频www| 免费无遮挡裸体视频| 欧美三级亚洲精品| 亚洲电影在线观看av| 每晚都被弄得嗷嗷叫到高潮| 男人舔女人的私密视频| 久久久久亚洲av毛片大全| 久久人人精品亚洲av| 久久亚洲精品不卡| 午夜免费观看网址| 久久香蕉国产精品| a级毛片a级免费在线| 91av网站免费观看| 日韩av在线大香蕉| 日韩高清综合在线| 午夜精品久久久久久毛片777| 免费在线观看成人毛片| 人妻久久中文字幕网| 亚洲 欧美一区二区三区| 亚洲电影在线观看av| 日韩欧美三级三区| 日韩有码中文字幕| 国产av麻豆久久久久久久| 日韩人妻高清精品专区| 91av网站免费观看| 蜜桃久久精品国产亚洲av| 757午夜福利合集在线观看| 亚洲美女黄片视频| 亚洲精品中文字幕一二三四区| 婷婷丁香在线五月| 国产视频一区二区在线看| 美女黄网站色视频| 亚洲人与动物交配视频| 亚洲中文日韩欧美视频| 91av网一区二区| 搡老妇女老女人老熟妇| 亚洲国产中文字幕在线视频| 欧美不卡视频在线免费观看| 中文字幕人妻丝袜一区二区| 很黄的视频免费| 亚洲av片天天在线观看| 国产真人三级小视频在线观看| 琪琪午夜伦伦电影理论片6080| 免费人成视频x8x8入口观看| 久久九九热精品免费| 女警被强在线播放| 色尼玛亚洲综合影院| 怎么达到女性高潮| av天堂中文字幕网| 日本五十路高清| 久久精品国产亚洲av香蕉五月| 深夜精品福利| 日韩三级视频一区二区三区| 国产一区二区三区在线臀色熟女| 成人特级av手机在线观看| 制服丝袜大香蕉在线| 成人午夜高清在线视频| 国产亚洲av嫩草精品影院| 99久久国产精品久久久| 成人亚洲精品av一区二区| 夜夜躁狠狠躁天天躁| 在线播放国产精品三级| 久久久久久久久中文| 一区二区三区激情视频| 国产精品免费一区二区三区在线| 一进一出抽搐动态| 亚洲国产精品合色在线| 99热只有精品国产| 很黄的视频免费| 老司机深夜福利视频在线观看| 精华霜和精华液先用哪个| 国产成人精品无人区| 久久伊人香网站| 日韩中文字幕欧美一区二区| 国产精品av久久久久免费| 亚洲第一欧美日韩一区二区三区| 一区二区三区激情视频| 亚洲国产欧美人成| www国产在线视频色| 国产精品一区二区三区四区免费观看 | 老司机深夜福利视频在线观看| 巨乳人妻的诱惑在线观看| 中亚洲国语对白在线视频| 亚洲人成网站高清观看| 日韩免费av在线播放| av在线天堂中文字幕| 国产v大片淫在线免费观看| 老司机在亚洲福利影院| 香蕉av资源在线| 两性夫妻黄色片| 窝窝影院91人妻| 国产三级在线视频| 国产成人影院久久av| 怎么达到女性高潮| 97人妻精品一区二区三区麻豆| 天天躁日日操中文字幕| 国产1区2区3区精品| av中文乱码字幕在线| 日日干狠狠操夜夜爽| 免费看日本二区| 一级毛片精品| 小说图片视频综合网站| 超碰成人久久| 成年女人看的毛片在线观看| 国产精品久久久久久亚洲av鲁大| 他把我摸到了高潮在线观看| 麻豆国产97在线/欧美| 久久久久久久精品吃奶| 久99久视频精品免费| 国内久久婷婷六月综合欲色啪| 中文字幕高清在线视频| 亚洲专区字幕在线| 亚洲成人免费电影在线观看| 19禁男女啪啪无遮挡网站| 99久久国产精品久久久| 亚洲欧美精品综合一区二区三区| 久久精品综合一区二区三区| 亚洲精品乱码久久久v下载方式 | 亚洲av成人av| 免费看日本二区| 日本成人三级电影网站| 久久久久久大精品| 国产高清videossex| 国产高清视频在线播放一区| 少妇裸体淫交视频免费看高清| 精品久久久久久久久久久久久| 999精品在线视频| 舔av片在线| 91麻豆av在线| 精品国产亚洲在线| 国产黄色小视频在线观看| 中文字幕最新亚洲高清| 我的老师免费观看完整版| 国内毛片毛片毛片毛片毛片| 久久午夜综合久久蜜桃| 色在线成人网| 麻豆国产av国片精品| 搡老熟女国产l中国老女人| 少妇丰满av| 国产亚洲av嫩草精品影院| 91在线观看av| 日韩欧美三级三区| e午夜精品久久久久久久| 亚洲专区国产一区二区| 国产精品女同一区二区软件 | 琪琪午夜伦伦电影理论片6080| 精品久久蜜臀av无| 国内少妇人妻偷人精品xxx网站 | 大型黄色视频在线免费观看| 欧美在线一区亚洲| 亚洲av片天天在线观看| 99精品欧美一区二区三区四区| aaaaa片日本免费| 国产黄片美女视频| 午夜福利免费观看在线| 色尼玛亚洲综合影院| 国内少妇人妻偷人精品xxx网站 | 99国产极品粉嫩在线观看| 久久天躁狠狠躁夜夜2o2o| 日日夜夜操网爽| 国产又色又爽无遮挡免费看| 国产成人aa在线观看| 白带黄色成豆腐渣| av在线天堂中文字幕| 中文字幕最新亚洲高清| www日本黄色视频网| 久久性视频一级片| 欧美日韩精品网址| 噜噜噜噜噜久久久久久91| 国产av一区在线观看免费| 午夜福利高清视频| 天天躁日日操中文字幕| 老司机深夜福利视频在线观看| 欧美午夜高清在线| 成年女人毛片免费观看观看9| 亚洲精品在线美女| 欧美激情在线99| 99久久精品一区二区三区| 一区福利在线观看| 国产成人欧美在线观看| 国产爱豆传媒在线观看| 国产亚洲精品久久久com| 亚洲av成人一区二区三| 国产精品,欧美在线| 美女高潮的动态| 动漫黄色视频在线观看| 午夜福利在线观看吧| 中文字幕av在线有码专区| 亚洲激情在线av| 亚洲欧美激情综合另类| 午夜视频精品福利| 亚洲五月天丁香| 91在线观看av| 久久久久久国产a免费观看| 亚洲国产精品久久男人天堂| 国产精品99久久99久久久不卡| 99久久综合精品五月天人人| 亚洲国产中文字幕在线视频| 久久久精品大字幕| 亚洲色图av天堂| 一a级毛片在线观看| 亚洲aⅴ乱码一区二区在线播放| 精品99又大又爽又粗少妇毛片 | 12—13女人毛片做爰片一| 亚洲av成人av| 禁无遮挡网站| 国产高清三级在线| 亚洲电影在线观看av| 国产成年人精品一区二区| 97超级碰碰碰精品色视频在线观看| 女同久久另类99精品国产91| 日韩欧美免费精品| 少妇的逼水好多| 1024香蕉在线观看| 99久国产av精品| 国内少妇人妻偷人精品xxx网站 | 欧美av亚洲av综合av国产av| 麻豆国产av国片精品| 网址你懂的国产日韩在线| 久久久精品欧美日韩精品| 亚洲在线观看片| 美女午夜性视频免费| 三级毛片av免费| 国产亚洲精品久久久久久毛片| 制服丝袜大香蕉在线| 国产av麻豆久久久久久久| 香蕉av资源在线| 婷婷丁香在线五月| 在线视频色国产色| 可以在线观看的亚洲视频| 黑人巨大精品欧美一区二区mp4| 久久中文字幕人妻熟女| 国产亚洲精品综合一区在线观看| 国产精品九九99| 色尼玛亚洲综合影院| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成a人片在线一区二区| 国产一级毛片七仙女欲春2| 国产成人精品无人区| 国产亚洲精品综合一区在线观看| 色综合站精品国产| 日韩大尺度精品在线看网址| 欧美黄色片欧美黄色片| 最近最新中文字幕大全电影3| 色在线成人网| 欧美日韩福利视频一区二区| 在线观看舔阴道视频| 又粗又爽又猛毛片免费看| 亚洲成av人片在线播放无| 一级黄色大片毛片| 美女黄网站色视频| 1024手机看黄色片| 精品久久久久久久人妻蜜臀av| 不卡av一区二区三区| 老熟妇仑乱视频hdxx| 90打野战视频偷拍视频| 婷婷丁香在线五月| 婷婷亚洲欧美| 热99在线观看视频| 性色avwww在线观看| 久久国产精品人妻蜜桃| 国产精品亚洲一级av第二区| 国产精品日韩av在线免费观看| 婷婷丁香在线五月| 麻豆久久精品国产亚洲av| 久久久久久久午夜电影| 99热这里只有精品一区 | 亚洲激情在线av| 18禁黄网站禁片午夜丰满| 老熟妇乱子伦视频在线观看| 十八禁网站免费在线| 桃色一区二区三区在线观看| 国产视频一区二区在线看| 99国产综合亚洲精品| av黄色大香蕉| 午夜福利在线观看免费完整高清在 | 午夜福利在线观看吧| 91老司机精品| 欧美一区二区精品小视频在线| 亚洲无线在线观看| av黄色大香蕉| 在线视频色国产色| 久久久久国内视频| 欧美日本视频| 美女黄网站色视频| 女警被强在线播放| 岛国在线观看网站| 在线免费观看的www视频| 国产精品美女特级片免费视频播放器 | 亚洲一区二区三区色噜噜| 女生性感内裤真人,穿戴方法视频| 国产成人一区二区三区免费视频网站| 亚洲国产精品合色在线| 长腿黑丝高跟| 91麻豆av在线| 熟女电影av网| 国产午夜精品久久久久久| 在线免费观看不下载黄p国产 | 动漫黄色视频在线观看| 亚洲av电影在线进入| 桃色一区二区三区在线观看| 999久久久精品免费观看国产| 他把我摸到了高潮在线观看| 可以在线观看的亚洲视频| 国产乱人伦免费视频| 久久香蕉国产精品| 婷婷丁香在线五月| 欧美激情久久久久久爽电影| 日韩欧美三级三区| 久久久国产成人免费| 美女被艹到高潮喷水动态| 婷婷精品国产亚洲av在线| 噜噜噜噜噜久久久久久91| 日本三级黄在线观看| 亚洲国产高清在线一区二区三| 51午夜福利影视在线观看| 国产成人欧美在线观看| 一个人看的www免费观看视频| 久久精品aⅴ一区二区三区四区| 亚洲欧美日韩高清专用| 欧美高清成人免费视频www| 国产真人三级小视频在线观看| 一级作爱视频免费观看| 精品久久久久久久人妻蜜臀av| 啦啦啦免费观看视频1| 国产精华一区二区三区| 99久久99久久久精品蜜桃| 成人三级黄色视频| 男女视频在线观看网站免费| 日本一本二区三区精品| 不卡一级毛片| 欧美zozozo另类| 亚洲人成网站高清观看| 精品国产三级普通话版| 无人区码免费观看不卡| 曰老女人黄片| 变态另类丝袜制服| 久久久久国产一级毛片高清牌| 欧美日韩综合久久久久久 | 亚洲18禁久久av| 女人被狂操c到高潮| 欧美日韩瑟瑟在线播放| 18禁国产床啪视频网站| 国产成人一区二区三区免费视频网站| 小说图片视频综合网站| 亚洲欧美一区二区三区黑人| 窝窝影院91人妻| 国产主播在线观看一区二区| 狂野欧美白嫩少妇大欣赏| 超碰成人久久| 免费电影在线观看免费观看| 国产麻豆成人av免费视频| 国产一区二区三区视频了| 国产成人精品无人区| 免费观看精品视频网站| 村上凉子中文字幕在线| 1024手机看黄色片| av天堂中文字幕网| 国产精品精品国产色婷婷| www.熟女人妻精品国产| 欧美日本亚洲视频在线播放| 亚洲精品久久国产高清桃花| 91在线观看av| 老汉色av国产亚洲站长工具| 97超级碰碰碰精品色视频在线观看| 老司机在亚洲福利影院| 十八禁网站免费在线| 丰满人妻熟妇乱又伦精品不卡| 好男人电影高清在线观看| 老鸭窝网址在线观看| 搡老熟女国产l中国老女人| 九九久久精品国产亚洲av麻豆 | 国内精品久久久久精免费| 人人妻人人澡欧美一区二区| 丝袜人妻中文字幕| 久久精品国产综合久久久| 国产蜜桃级精品一区二区三区| 美女高潮的动态| 99久久久亚洲精品蜜臀av| 国产精品精品国产色婷婷| 亚洲av电影不卡..在线观看| 12—13女人毛片做爰片一| 最近最新免费中文字幕在线| 国产探花在线观看一区二区| 久久久成人免费电影| 又爽又黄无遮挡网站| 熟妇人妻久久中文字幕3abv| 亚洲av成人av| 精品99又大又爽又粗少妇毛片 | 色在线成人网| 国产成人影院久久av| 日日夜夜操网爽| 搡老岳熟女国产| 亚洲av五月六月丁香网| 成年女人永久免费观看视频| 国产成人欧美在线观看| 亚洲国产高清在线一区二区三| 国产高清视频在线观看网站| 国产视频内射| 一级黄色大片毛片| 久久中文字幕人妻熟女| 国产 一区 欧美 日韩| 丁香欧美五月| 在线看三级毛片| 日本黄色视频三级网站网址| 免费大片18禁| 两个人看的免费小视频| 国产野战对白在线观看| 后天国语完整版免费观看| 精品福利观看| 国产精品精品国产色婷婷| 国产1区2区3区精品| 亚洲美女视频黄频| 最新在线观看一区二区三区| 国产精品99久久99久久久不卡| 偷拍熟女少妇极品色| 窝窝影院91人妻| 宅男免费午夜| 成人亚洲精品av一区二区| 极品教师在线免费播放| 亚洲欧洲精品一区二区精品久久久| 国产精品香港三级国产av潘金莲| 老司机福利观看| 最好的美女福利视频网| 热99re8久久精品国产| 亚洲 欧美 日韩 在线 免费| 观看美女的网站| 一区二区三区国产精品乱码| 美女高潮喷水抽搐中文字幕| 亚洲第一欧美日韩一区二区三区| 国产av在哪里看| 国产精品美女特级片免费视频播放器 | 美女大奶头视频| 99久久综合精品五月天人人| 午夜福利18| 热99re8久久精品国产| 国产爱豆传媒在线观看| 亚洲欧美日韩高清专用| 欧美日韩综合久久久久久 | 99热只有精品国产| 久久久国产成人精品二区| 18禁观看日本| 可以在线观看的亚洲视频| 亚洲专区中文字幕在线| 久99久视频精品免费| 黑人巨大精品欧美一区二区mp4| 最近最新中文字幕大全免费视频| 两性午夜刺激爽爽歪歪视频在线观看| 欧美性猛交╳xxx乱大交人| 国产精品野战在线观看| 岛国视频午夜一区免费看| 最近最新中文字幕大全电影3| 国产又色又爽无遮挡免费看| 欧美丝袜亚洲另类 | 午夜两性在线视频| 亚洲国产欧洲综合997久久,| 黄色丝袜av网址大全| 国产一区二区在线观看日韩 | 久久久久九九精品影院| 国产在线精品亚洲第一网站| 一区福利在线观看| 成熟少妇高潮喷水视频| 国产成人系列免费观看| 又爽又黄无遮挡网站| 久久亚洲精品不卡| 精品久久久久久久久久久久久| 午夜影院日韩av| 又爽又黄无遮挡网站| 精品久久久久久,| 很黄的视频免费| 精品久久久久久,| 少妇的丰满在线观看| 久久人人精品亚洲av| 国内精品久久久久久久电影| 国产高清视频在线播放一区| 亚洲成av人片免费观看| 国产午夜福利久久久久久| 老鸭窝网址在线观看| 网址你懂的国产日韩在线| 51午夜福利影视在线观看| 国产精品一区二区三区四区免费观看 | 91老司机精品| 深夜精品福利| 观看免费一级毛片| 国产精品久久久人人做人人爽| 伦理电影免费视频| 午夜免费观看网址| 级片在线观看| 久久天堂一区二区三区四区| 亚洲av电影不卡..在线观看| 少妇人妻一区二区三区视频| x7x7x7水蜜桃| av中文乱码字幕在线| 又爽又黄无遮挡网站| 久久久国产欧美日韩av| av片东京热男人的天堂| 女人高潮潮喷娇喘18禁视频| 桃色一区二区三区在线观看| 久久国产乱子伦精品免费另类| 免费看日本二区| 久久精品夜夜夜夜夜久久蜜豆| 国产精品 欧美亚洲| 久久精品aⅴ一区二区三区四区| 亚洲第一电影网av| 91麻豆精品激情在线观看国产| 久久人妻av系列| 免费高清视频大片| 麻豆一二三区av精品| 亚洲狠狠婷婷综合久久图片| 欧美乱色亚洲激情| 51午夜福利影视在线观看| 免费无遮挡裸体视频| 久久久久国内视频| 精品国产超薄肉色丝袜足j| 国产亚洲精品久久久com| 国产高清视频在线播放一区| 精品一区二区三区四区五区乱码| 男女午夜视频在线观看| 99久国产av精品| 婷婷六月久久综合丁香| 中文字幕高清在线视频| 国产成人影院久久av| 亚洲成人精品中文字幕电影| 国内精品久久久久精免费| 欧美zozozo另类| 欧美+亚洲+日韩+国产| 99久久成人亚洲精品观看| 俄罗斯特黄特色一大片| 一区二区三区激情视频| 亚洲第一欧美日韩一区二区三区| 床上黄色一级片| 亚洲欧美精品综合一区二区三区| xxxwww97欧美| 99热只有精品国产| 亚洲国产精品合色在线| 狂野欧美白嫩少妇大欣赏| 欧美午夜高清在线| 欧美绝顶高潮抽搐喷水| 欧美日韩一级在线毛片| 精品国产美女av久久久久小说| 99久久99久久久精品蜜桃| 美女cb高潮喷水在线观看 | 色尼玛亚洲综合影院| 99国产综合亚洲精品| 全区人妻精品视频| 成人av在线播放网站| 久久国产乱子伦精品免费另类| 日本三级黄在线观看| 天堂√8在线中文| 亚洲无线观看免费| 亚洲精华国产精华精| 两个人的视频大全免费| 国产69精品久久久久777片 | 99久久成人亚洲精品观看| 亚洲精品美女久久av网站| 欧美日本视频| 久久亚洲真实| 国产精品综合久久久久久久免费| 国产视频内射| 精品国产三级普通话版| 高潮久久久久久久久久久不卡| 国产精品爽爽va在线观看网站| 日本熟妇午夜| 日本一二三区视频观看| 久久久久久人人人人人| 国产亚洲欧美98| 在线永久观看黄色视频| 丁香六月欧美| 日日干狠狠操夜夜爽| 亚洲 欧美 日韩 在线 免费| 非洲黑人性xxxx精品又粗又长| 日韩 欧美 亚洲 中文字幕| 真实男女啪啪啪动态图| 亚洲狠狠婷婷综合久久图片| 白带黄色成豆腐渣| 麻豆成人午夜福利视频| 国产成人系列免费观看| 九色国产91popny在线| 亚洲av日韩精品久久久久久密|