• <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ī)學救援程序與訓練思考
    高清不卡的av网站| 欧美黑人精品巨大| 51午夜福利影视在线观看| 国产熟女午夜一区二区三区| 制服人妻中文乱码| 手机成人av网站| 精品人妻1区二区| 高清视频免费观看一区二区| 国产色视频综合| 丰满人妻熟妇乱又伦精品不卡| 亚洲av片天天在线观看| 女人被躁到高潮嗷嗷叫费观| 五月开心婷婷网| 免费观看a级毛片全部| 搡老乐熟女国产| a 毛片基地| 国产精品久久久av美女十八| 国产福利在线免费观看视频| 中文精品一卡2卡3卡4更新| 中文字幕高清在线视频| 亚洲专区中文字幕在线| 欧美激情 高清一区二区三区| 久久99一区二区三区| 两人在一起打扑克的视频| 黄色毛片三级朝国网站| 乱人伦中国视频| 久久鲁丝午夜福利片| 久久99精品国语久久久| 午夜免费鲁丝| 国产高清视频在线播放一区 | 美女脱内裤让男人舔精品视频| 久9热在线精品视频| 国产日韩一区二区三区精品不卡| 亚洲精品在线美女| 啦啦啦 在线观看视频| 只有这里有精品99| 国产视频首页在线观看| 国产无遮挡羞羞视频在线观看| 国产精品国产av在线观看| 国产成人精品久久二区二区91| 女人精品久久久久毛片| 欧美性长视频在线观看| 国产真人三级小视频在线观看| 人体艺术视频欧美日本| 啦啦啦在线观看免费高清www| 亚洲国产最新在线播放| 精品亚洲成国产av| 夫妻性生交免费视频一级片| 各种免费的搞黄视频| 精品久久久精品久久久| 飞空精品影院首页| 一区二区日韩欧美中文字幕| 亚洲av国产av综合av卡| 黑人猛操日本美女一级片| 欧美日韩av久久| 日本欧美视频一区| 伊人亚洲综合成人网| 熟女少妇亚洲综合色aaa.| 一区二区三区乱码不卡18| 自线自在国产av| 在线天堂中文资源库| 69精品国产乱码久久久| 大片免费播放器 马上看| 一边摸一边做爽爽视频免费| 亚洲精品美女久久av网站| 纵有疾风起免费观看全集完整版| 日本午夜av视频| 在线精品无人区一区二区三| 免费在线观看影片大全网站 | 另类精品久久| www.自偷自拍.com| 男人添女人高潮全过程视频| 欧美黑人欧美精品刺激| 香蕉丝袜av| 亚洲精品日韩在线中文字幕| 免费在线观看黄色视频的| av片东京热男人的天堂| 久久久精品94久久精品| 国产精品欧美亚洲77777| 亚洲欧美日韩高清在线视频 | 99re6热这里在线精品视频| 亚洲一码二码三码区别大吗| 欧美日韩福利视频一区二区| 国产精品亚洲av一区麻豆| 国产片内射在线| 欧美精品亚洲一区二区| 热99久久久久精品小说推荐| 久久99精品国语久久久| 女人精品久久久久毛片| 国产精品免费视频内射| 男女床上黄色一级片免费看| 国产精品偷伦视频观看了| 一区二区av电影网| 老司机深夜福利视频在线观看 | 国产又色又爽无遮挡免| 男的添女的下面高潮视频| 日本猛色少妇xxxxx猛交久久| 大片免费播放器 马上看| 99久久人妻综合| 久久女婷五月综合色啪小说| 一边亲一边摸免费视频| 久久精品久久久久久久性| 人人妻人人爽人人添夜夜欢视频| 日本wwww免费看| 天堂8中文在线网| 啦啦啦啦在线视频资源| 日韩熟女老妇一区二区性免费视频| 国产麻豆69| 99久久99久久久精品蜜桃| 亚洲av在线观看美女高潮| 精品久久久久久久毛片微露脸 | 久久午夜综合久久蜜桃| 一区二区三区激情视频| 免费高清在线观看日韩| 少妇的丰满在线观看| 亚洲精品日本国产第一区| 美女中出高潮动态图| 男女下面插进去视频免费观看| 一级片'在线观看视频| 少妇 在线观看| 成年动漫av网址| 亚洲国产欧美一区二区综合| 人人澡人人妻人| 老司机靠b影院| 男女床上黄色一级片免费看| 99热网站在线观看| 老司机影院成人| 久久久精品免费免费高清| 欧美亚洲 丝袜 人妻 在线| 两个人看的免费小视频| 亚洲av在线观看美女高潮| 51午夜福利影视在线观看| 日本猛色少妇xxxxx猛交久久| 夜夜骑夜夜射夜夜干| 欧美97在线视频| 麻豆av在线久日| av天堂久久9| 性色av一级| 美女高潮到喷水免费观看| 99久久精品国产亚洲精品| 亚洲国产精品成人久久小说| 久久精品国产亚洲av涩爱| 免费在线观看黄色视频的| 亚洲欧美一区二区三区久久| 免费看不卡的av| 亚洲精品久久久久久婷婷小说| svipshipincom国产片| 韩国高清视频一区二区三区| 叶爱在线成人免费视频播放| www.熟女人妻精品国产| 丝袜在线中文字幕| 日本a在线网址| 久久久久精品国产欧美久久久 | 久久精品国产a三级三级三级| 久久人妻熟女aⅴ| 国产男人的电影天堂91| 热re99久久精品国产66热6| 日韩,欧美,国产一区二区三区| 一区二区av电影网| 国产片特级美女逼逼视频| 亚洲国产精品999| 中文字幕人妻熟女乱码| 欧美精品啪啪一区二区三区 | 夫妻性生交免费视频一级片| 国产亚洲欧美精品永久| 免费看av在线观看网站| 色94色欧美一区二区| 18在线观看网站| 国产高清国产精品国产三级| 欧美日韩黄片免| 国产精品九九99| 亚洲av成人精品一二三区| 男人爽女人下面视频在线观看| 国产伦理片在线播放av一区| 亚洲av男天堂| 天堂俺去俺来也www色官网| 亚洲久久久国产精品| 少妇人妻 视频| 老司机亚洲免费影院| 亚洲人成网站在线观看播放| 十分钟在线观看高清视频www| 下体分泌物呈黄色| 国产三级黄色录像| 午夜久久久在线观看| 日韩一区二区三区影片| av有码第一页| 久久久国产一区二区| 超碰97精品在线观看| 免费高清在线观看日韩| 中文字幕色久视频| 国产免费视频播放在线视频| 手机成人av网站| 欧美日韩亚洲高清精品| 免费少妇av软件| videosex国产| 亚洲天堂av无毛| 日韩一卡2卡3卡4卡2021年| 欧美国产精品一级二级三级| 亚洲欧美中文字幕日韩二区| 巨乳人妻的诱惑在线观看| 久久女婷五月综合色啪小说| 免费观看a级毛片全部| 久久国产精品人妻蜜桃| 国产精品99久久99久久久不卡| tube8黄色片| 国产有黄有色有爽视频| 色婷婷av一区二区三区视频| 亚洲情色 制服丝袜| 久久精品久久久久久久性| 岛国毛片在线播放| 欧美国产精品一级二级三级| 在线观看人妻少妇| 亚洲一卡2卡3卡4卡5卡精品中文| 精品少妇一区二区三区视频日本电影| 99国产综合亚洲精品| 黑人猛操日本美女一级片| 在线亚洲精品国产二区图片欧美| 18在线观看网站| 91精品伊人久久大香线蕉| 在线观看免费视频网站a站| 久久人人爽人人片av| 1024香蕉在线观看| 一本—道久久a久久精品蜜桃钙片| 日韩 亚洲 欧美在线| 中文字幕人妻丝袜一区二区| 少妇人妻 视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品美女久久久久99蜜臀 | 精品人妻在线不人妻| 婷婷色综合大香蕉| 久久久久国产精品人妻一区二区| 久久人妻熟女aⅴ| 久久国产精品人妻蜜桃| 在线观看人妻少妇| 亚洲av成人精品一二三区| 777久久人妻少妇嫩草av网站| 高清视频免费观看一区二区| 18在线观看网站| 色播在线永久视频| 欧美黑人欧美精品刺激| 日本猛色少妇xxxxx猛交久久| 嫩草影视91久久| 香蕉丝袜av| 黑人欧美特级aaaaaa片| 成人午夜精彩视频在线观看| 两人在一起打扑克的视频| 亚洲欧美日韩另类电影网站| 亚洲国产精品一区二区三区在线| 天堂8中文在线网| 日本欧美视频一区| 一级毛片 在线播放| 日韩中文字幕欧美一区二区 | 欧美性长视频在线观看| 老汉色av国产亚洲站长工具| 大码成人一级视频| av线在线观看网站| 免费观看a级毛片全部| 久久久久网色| 色网站视频免费| 日本av手机在线免费观看| 欧美97在线视频| 国产欧美日韩综合在线一区二区| 性色av一级| 久久精品国产亚洲av高清一级| 国产成人av激情在线播放| 青春草视频在线免费观看| 日本一区二区免费在线视频| 韩国高清视频一区二区三区| 51午夜福利影视在线观看| 真人做人爱边吃奶动态| 国产精品二区激情视频| 丰满人妻熟妇乱又伦精品不卡| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品秋霞免费鲁丝片| 一级黄色大片毛片| 久久中文字幕一级| 亚洲成人手机| 一本一本久久a久久精品综合妖精| 777米奇影视久久| 天天躁夜夜躁狠狠躁躁| 久久性视频一级片| 欧美成人午夜精品| a级毛片在线看网站| 精品卡一卡二卡四卡免费| 国产精品二区激情视频| 亚洲欧洲精品一区二区精品久久久| 女人高潮潮喷娇喘18禁视频| 中文字幕人妻熟女乱码| 免费看av在线观看网站| 亚洲人成电影观看| 一级a爱视频在线免费观看| 大陆偷拍与自拍| 黄片播放在线免费| 蜜桃在线观看..| 亚洲成国产人片在线观看| 婷婷成人精品国产| 亚洲五月婷婷丁香| 一级黄片播放器| 狠狠婷婷综合久久久久久88av| 亚洲,一卡二卡三卡| 亚洲人成网站在线观看播放| 日本猛色少妇xxxxx猛交久久| 精品亚洲乱码少妇综合久久| 国产成人av教育| 精品国产超薄肉色丝袜足j| 99国产精品一区二区蜜桃av | 久久精品成人免费网站| 亚洲欧美中文字幕日韩二区| 亚洲国产精品一区二区三区在线| 欧美国产精品va在线观看不卡| 性高湖久久久久久久久免费观看| 亚洲欧美日韩高清在线视频 | av欧美777| 性少妇av在线| 自线自在国产av| 成年人黄色毛片网站| 一级黄片播放器| 大型av网站在线播放| 咕卡用的链子| 色94色欧美一区二区| 久久精品久久精品一区二区三区| av欧美777| 亚洲精品乱久久久久久| 精品亚洲成国产av| 一级毛片 在线播放| 三上悠亚av全集在线观看| 久久九九热精品免费| 中文乱码字字幕精品一区二区三区| 国产欧美日韩一区二区三 | 日本av手机在线免费观看| 少妇被粗大的猛进出69影院| 国产视频首页在线观看| av视频免费观看在线观看| 日韩制服骚丝袜av| 成年av动漫网址| 久久99精品国语久久久| 亚洲三区欧美一区| 一区二区三区精品91| 狂野欧美激情性bbbbbb| 男女下面插进去视频免费观看| 我要看黄色一级片免费的| 久久女婷五月综合色啪小说| 久久青草综合色| 亚洲久久久国产精品| 搡老乐熟女国产| 欧美日韩亚洲高清精品| 狠狠婷婷综合久久久久久88av| 亚洲综合色网址| 777久久人妻少妇嫩草av网站| 欧美av亚洲av综合av国产av| 亚洲激情五月婷婷啪啪| 男女国产视频网站| 三上悠亚av全集在线观看| 免费观看av网站的网址| 精品一区在线观看国产| 美女视频免费永久观看网站| 精品国产乱码久久久久久男人| 午夜久久久在线观看| 好男人电影高清在线观看| 老鸭窝网址在线观看| 丁香六月欧美| 黑人巨大精品欧美一区二区蜜桃| 亚洲人成77777在线视频| 日本猛色少妇xxxxx猛交久久| 中文字幕亚洲精品专区| 国产免费福利视频在线观看| 久久ye,这里只有精品| av天堂在线播放| 婷婷色麻豆天堂久久| 午夜日韩欧美国产| 日韩欧美一区视频在线观看| av网站在线播放免费| 国产在线观看jvid| 成人三级做爰电影| 亚洲第一av免费看| 婷婷丁香在线五月| 亚洲精品久久成人aⅴ小说| 十分钟在线观看高清视频www| 亚洲av成人不卡在线观看播放网 | 亚洲,一卡二卡三卡| 欧美精品亚洲一区二区| 国产福利在线免费观看视频| 下体分泌物呈黄色| 交换朋友夫妻互换小说| a 毛片基地| av网站免费在线观看视频| a 毛片基地| 纵有疾风起免费观看全集完整版| 777久久人妻少妇嫩草av网站| 一级片免费观看大全| 亚洲国产精品成人久久小说| 一级片免费观看大全| 一二三四在线观看免费中文在| 母亲3免费完整高清在线观看| 一边摸一边抽搐一进一出视频| 伊人亚洲综合成人网| 亚洲国产精品成人久久小说| 亚洲欧美一区二区三区黑人| 国产精品99久久99久久久不卡| 亚洲av欧美aⅴ国产| 男人爽女人下面视频在线观看| av网站免费在线观看视频| 亚洲国产精品成人久久小说| 亚洲精品自拍成人| 另类精品久久| 日本一区二区免费在线视频| 激情视频va一区二区三区| 国产精品成人在线| 亚洲精品中文字幕在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 精品亚洲乱码少妇综合久久| 日韩中文字幕视频在线看片| 两个人免费观看高清视频| √禁漫天堂资源中文www| 久久亚洲精品不卡| 久久综合国产亚洲精品| 国产在线视频一区二区| 国产成人a∨麻豆精品| 黑丝袜美女国产一区| 9191精品国产免费久久| 黄频高清免费视频| 亚洲欧美一区二区三区国产| 亚洲av国产av综合av卡| 国产精品久久久久久精品古装| 欧美人与性动交α欧美精品济南到| 欧美 日韩 精品 国产| 免费在线观看黄色视频的| 亚洲国产精品国产精品| 成人国产av品久久久| 一级a爱视频在线免费观看| 亚洲一区二区三区欧美精品| 啦啦啦啦在线视频资源| 亚洲av成人精品一二三区| 精品人妻1区二区| 久久久久国产一级毛片高清牌| 日韩一本色道免费dvd| 久久久久网色| 汤姆久久久久久久影院中文字幕| 免费看av在线观看网站| 女警被强在线播放| 精品亚洲成国产av| 免费av中文字幕在线| 国产91精品成人一区二区三区 | 夫妻午夜视频| 亚洲国产精品999| 青青草视频在线视频观看| 欧美激情极品国产一区二区三区| 91九色精品人成在线观看| 满18在线观看网站| 一二三四在线观看免费中文在| 青草久久国产| 在线亚洲精品国产二区图片欧美| 99久久精品国产亚洲精品| 精品一区二区三卡| 久久影院123| 国产一区二区三区综合在线观看| 国语对白做爰xxxⅹ性视频网站| 久久国产精品大桥未久av| www日本在线高清视频| 亚洲黑人精品在线| 精品国产超薄肉色丝袜足j| 国产国语露脸激情在线看| 97在线人人人人妻| 啦啦啦在线观看免费高清www| 国产又色又爽无遮挡免| 日日夜夜操网爽| 人体艺术视频欧美日本| 色婷婷久久久亚洲欧美| 国产黄色免费在线视频| 亚洲国产成人一精品久久久| 欧美黄色淫秽网站| 午夜福利乱码中文字幕| 久久青草综合色| 亚洲人成77777在线视频| 中文乱码字字幕精品一区二区三区| 1024视频免费在线观看| 女人被躁到高潮嗷嗷叫费观| 可以免费在线观看a视频的电影网站| 国产主播在线观看一区二区 | 最近中文字幕2019免费版| 日本欧美国产在线视频| 久久99精品国语久久久| www.自偷自拍.com| av网站在线播放免费| 亚洲av美国av| 9191精品国产免费久久| 天天躁夜夜躁狠狠久久av| 建设人人有责人人尽责人人享有的| 最黄视频免费看| 亚洲中文字幕日韩| 日本wwww免费看| 精品国产一区二区三区四区第35| 50天的宝宝边吃奶边哭怎么回事| 久久人妻福利社区极品人妻图片 | 视频在线观看一区二区三区| 飞空精品影院首页| 亚洲国产欧美在线一区| 国产又爽黄色视频| 黑人欧美特级aaaaaa片| 久久中文字幕一级| 亚洲成人免费电影在线观看 | 国产淫语在线视频| 我的亚洲天堂| 中文字幕另类日韩欧美亚洲嫩草| 欧美人与性动交α欧美软件| 亚洲五月婷婷丁香| 涩涩av久久男人的天堂| av一本久久久久| 女性生殖器流出的白浆| 久久久久网色| 中文字幕人妻熟女乱码| 午夜影院在线不卡| 黑人欧美特级aaaaaa片| 成在线人永久免费视频| 秋霞在线观看毛片| 一区二区av电影网| 少妇被粗大的猛进出69影院| 久久国产亚洲av麻豆专区| cao死你这个sao货| 免费在线观看影片大全网站 | 国产爽快片一区二区三区| bbb黄色大片| 乱人伦中国视频| 亚洲精品国产一区二区精华液| 欧美人与善性xxx| 又黄又粗又硬又大视频| 亚洲激情五月婷婷啪啪| 成年av动漫网址| 狂野欧美激情性xxxx| 国产精品久久久人人做人人爽| 美女福利国产在线| 男女免费视频国产| 老汉色∧v一级毛片| 在线av久久热| 亚洲,欧美精品.| 久久青草综合色| 国产精品av久久久久免费| 男人爽女人下面视频在线观看| 99香蕉大伊视频| 天天影视国产精品| 9色porny在线观看| 成在线人永久免费视频| 国产精品熟女久久久久浪| 久久性视频一级片| 国产亚洲精品第一综合不卡| 久久人人爽av亚洲精品天堂| 叶爱在线成人免费视频播放| 麻豆乱淫一区二区| 少妇人妻久久综合中文| 少妇被粗大的猛进出69影院| 亚洲中文日韩欧美视频| 男女免费视频国产| 午夜精品国产一区二区电影| 免费久久久久久久精品成人欧美视频| 久久鲁丝午夜福利片| 国产精品熟女久久久久浪| 69精品国产乱码久久久| 日韩视频在线欧美| 中文字幕av电影在线播放| cao死你这个sao货| 欧美+亚洲+日韩+国产| 一边亲一边摸免费视频| 亚洲中文av在线| 久久久久国产精品人妻一区二区| 国产精品秋霞免费鲁丝片| kizo精华| av电影中文网址| 国产精品偷伦视频观看了| 国产亚洲av片在线观看秒播厂| 丰满迷人的少妇在线观看| 大陆偷拍与自拍| 国产日韩欧美视频二区| 成人国语在线视频| 久久亚洲国产成人精品v| 欧美 亚洲 国产 日韩一| 如日韩欧美国产精品一区二区三区| 国产真人三级小视频在线观看| 久久精品久久精品一区二区三区| 大型av网站在线播放| www.自偷自拍.com| 国产成人啪精品午夜网站| 亚洲av日韩精品久久久久久密 | 欧美xxⅹ黑人| 丰满饥渴人妻一区二区三| 亚洲成国产人片在线观看| 国产一区二区三区综合在线观看| www日本在线高清视频| 又大又黄又爽视频免费| 国产不卡av网站在线观看| 又黄又粗又硬又大视频| 中文字幕人妻熟女乱码| 考比视频在线观看| 日韩免费高清中文字幕av| 男人操女人黄网站| 国产淫语在线视频| 91精品伊人久久大香线蕉| 嫁个100分男人电影在线观看 | 9色porny在线观看| 日本a在线网址| 久久久久久久大尺度免费视频| 一级毛片电影观看| 美女脱内裤让男人舔精品视频| 精品国产超薄肉色丝袜足j| 91麻豆av在线| 国产精品一区二区在线观看99| 亚洲成人免费电影在线观看 | 亚洲久久久国产精品| 国产片特级美女逼逼视频| 国产国语露脸激情在线看| 国产伦理片在线播放av一区| 日本欧美国产在线视频| 亚洲av成人精品一二三区| 久久亚洲国产成人精品v| 日本欧美国产在线视频| 国产精品秋霞免费鲁丝片|