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

    一種浸潤凍結機制冰核測量裝置(FINDA)的搭建與應用*

    2021-11-12 09:17:04丁德平陳景川黃夢宇
    氣象學報 2021年5期
    關鍵詞:超純水液滴測量

    畢 凱 丁德平 楊 帥 陳景川 黃夢宇

    1. 中國氣象局華北云降水野外科學試驗基地,北京,101200

    2. 云降水物理研究和云水資源開發(fā)北京市重點實驗室,北京,100089

    3. 北京市人工影響天氣中心,北京,100089

    4. 北京大學環(huán)境科學與工程學院,北京,100871

    1 引 言

    大氣冰核是指懸浮在大氣中,能使過飽和水汽凝華或者過冷水滴凍結形成初始冰晶的微粒。它的存在大幅度提高了云中成冰溫度(Kanji,et al,2017),從而影響冷云生命周期、相態(tài)、光學特征和微物理特征(Beall,et al,2017),是全球氣候效應準確評估最大的不確定性來源之一(IPCC,2013)。對冰核的準確測量和結冰機理的認識直接影響對云降水的科學認識和模擬研究,進而影響天氣預報和氣候模式的準確度,制約著人工影響天氣作業(yè)、氣象災害防治(冰雹、強降水等)和氣候變化評估。

    20世紀40 年代開始,國際上有學者針對冰核不同成冰機制設計了多種裝置開展室內試驗和外場觀測。冰核的異質核化主要有4種機制,分別是凝華核化(水汽在冰核表面凝華成冰)、凝結凍結(水汽在冰核表面凝結而后凍結成冰)、浸潤凍結(冰核嵌入過冷水滴內使其轉換成一個冰粒子)和接觸凍結(過冷水滴同冰核碰撞并形成冰粒子)(盛裴軒等,2003)。研究表明,浸潤凍結是混合相態(tài)云中最主要的核化機制(Ansmann,et al,2008;Hoose,et al,2010;Murray,et al,2012),因此目前國際上對冰核的研究大部分是基于浸潤凍結機制開展的,測量方式分為在線云室法和離線凍滴法。在線云室法以在線連續(xù)流擴散云室(Continuous-Flow Diffusion Chamber,簡稱:CFDC)為主(Rogers,et al,2001;DeMott,et al,2003;Bi,et al,2019),移動的液滴經(jīng)過溫、濕度可控的云室發(fā)生凍結及增長,根據(jù)尺度差利用光電計數(shù)器檢測。另一種是離線凍滴法,采用液體樣本或用濾膜把大氣氣溶膠采集并配置成懸濁液,取大量等體積微小液滴放置在冷臺介質上進行降溫凍結,凍結比率和冰核積分濃度譜根據(jù)Vali(1971)提出的理論進行計算。基于離線凍滴法,學者們搭建了不同降溫率、不同液滴尺度和數(shù)量、不同凍結識別方式、不同冷臺尺寸的離線測量裝置,例如CSU-IS (Garcia,et al,2012)、INC (Wang,et al,2012)、LINDA (Stopelli,et al,2014)、μL-NIPI (Whale,et al,2015)、KIT-INSEKT(Schiebel,2017)、SIO-AIS (Beall,et al,2017)、IR-NIPI (Harrison,et al,2018)和DRINCZ (David,et al,2019)等設備,用于測量大氣、降水或特定物質中不同溫度區(qū)間的冰核濃度和活化溫度。

    中國從20世紀60年代開始冰核的觀測研究,普遍采用仿比格混合云室進行測量(游來光等,1964,2002;吳明林等,1986;牛生杰等,2000;周德平等,2012a;楊磊等,2013;Bi,et al,2018;Che,et al,2019,2021),但由于比格混合云室測量機制難確定,無法精確控溫控濕,限制了其應用。80年代后中國有學者制作了小型靜力擴散云室,通過濾膜在外場采集冰核后進行離線分析凝華和凝結凍結冰核(章澄昌,1987;楊紹忠等,1995;周德平等,2012b,2018;Li,et al,2017)。酆大雄等(2002)在靜力擴散云室冷臺的濾膜上滴加水滴測量了降水樣品中凍結核的含量,并提出可以通過凍結實驗測量不同區(qū)域降水中冰核含量的方法檢驗人工影響天氣作業(yè)效果。楊紹忠(2005)改進了水滴凍結實驗設備并進行了降水樣品中的冰核檢測,搭建了一種基于均勻水滴的平面冷臺裝置(楊紹忠等,2007),并被用于降水中生物冰核的檢測(Lu,et al,2016)。蘇航等(2014)在德國FRIDGE裝置(Klein,et al,2010)的基礎上改良搭建了新型擴散云室用于測量凝華或凝結凍結機制的冰核特征。Wang等(2012)搭建了基于單顆粒的冰核凝華和浸潤檢測裝置。Chen等(2018)基于顯微冷臺搭建了基于1 μl液滴的冷臺凍滴裝置測量浸潤機制冰核濃度。2017年,北京市人工影響天氣中心引進了新型在線連續(xù)流擴散云室(Continuous-Flow Diffusion Chamber Ice Activation Spectrometer,簡稱BJ-CFDC),并基于浸潤核化模式開展了北京區(qū)域冰核的連續(xù)流冰核外場觀測(Bi,et al,2019;畢凱等,2020a)。

    在線云室法和離線凍滴法各有優(yōu)缺點,目前還沒有任何設備測量能覆蓋混合相態(tài)云中溫度范圍和冰核濃度范圍(Harrison,et al,2018)。CFDC類在線云室法設備可以精確控制溫、濕度,并且實時測量,但受溫、濕度控制的限制,通常無法測量高活化溫度(高于-15℃)和低濃度(小于0.1 l-1)(DeMott,et al,2010)的冰核。離線凍滴法由于可長時間累積采樣(通常為幾小時至幾天),能夠捕獲并檢測高活化溫度和低濃度的冰核,但時間分辨率較低。液滴凍結受尺度的影響,使用大體積液滴(幾十微升)可以有效測量-15℃以上的冰核(Vali,1971),已經(jīng)被 用 于CSU-IS (Garcia,et al,2012)、DRINCZ(David,et al,2019)和IR-NIPI (Harrison,et al,2018)等多種基于聚丙烯酶標板(PCR板)的冷臺測量裝置上。目前已有的大滴凍滴設備通常測量冷臺底部溫度或者冷臺邊緣溫度代表液滴溫度(Hill,et al,2014;Beall,et al,2017),但由于冷臺面板的水平和垂直溫差,以及PCR板熱傳導差異,難以代表液滴實際溫度。LINDA設備通過把Pt1000傳感器改裝到單孔PCR中測溫以消除PCR板的溫度傳導造成的測量偏差,但使用了透光式凍結識別,需要在樣品中添加氯化鈉以提升識別效率(Stopelli,et al,2014)。Beall等(2017)通過模擬研究了冷臺垂直溫度偏差,Harrison等(2018)和Kunert等(2018)使用單一物質的凍結對冷臺某溫度下的水平溫度偏差進行了修正,但均未針對不同溫度下的溫差進行動態(tài)評估及修正。

    基于以上考慮,文中介紹了一種新的浸潤凍結核化機制冰核測量裝置(Freezing Ice Nucleation Detector Array,簡稱FINDA),用于檢測浸潤凍結核化機制的冰核數(shù)濃度譜,定量評估降水(雨、雪、雹等)和大氣中的冰核濃度。本研究對多通道溫度傳感器進行了改裝校準,設計了基于96孔聚丙烯酶標板(PCR板)的密封式實驗冷臺,并利用紅外熱成像儀對冷臺樣品水平溫差進行了測量及動態(tài)溫度校準,對降溫控制、凍結識別及數(shù)據(jù)解析進行了自動化軟件集成。作為對該裝置的性能檢驗,進行了超純水的凍結溫度檢驗,測量并分析討論了北京地區(qū)兩次降雪樣品的冰核譜。另外,通過濾膜采樣分析檢測了北京海坨山區(qū)大氣冰核的數(shù)濃度譜并與在線連續(xù)流擴散云室(BJ-CFDC)進行了對比。

    2 測量原理

    FINDA用于測量液體樣品中0—-25℃間隔0.25℃的冰核數(shù)濃度譜。測量原理是基于冷臺凍滴實驗定量反演冰核數(shù)濃度(Vali,1971,1994),即通過測量大量等體積小液滴的凍結比率隨溫度的變化,利用統(tǒng)計分析方法反演得到原始液體樣本中的冰核積分數(shù)濃度譜,已應用于多種大氣冰核測量裝置(楊紹忠等,2007;Hill,et al,2014;Chen,et al,2018;David,et al,2019)。這里的冰核積分數(shù)濃度是指高于該溫度下液體中含有的冰核數(shù)濃度之和(Vali,1971)。該方法基于以下假設:(1)所有取樣液滴體積相等,(2)取樣液滴和原始樣本中顆粒物數(shù)濃度相等,(3)每個顆粒物都有一個對應的特征凍結溫度,該溫度下能夠觸發(fā)液滴凍結,即忽略降溫率的影響,(4)顆粒物的特征凍結溫度在統(tǒng)計學上是隨機且獨立的,不受液滴中其他顆粒物的影響,(5)液滴凍結由內部包含的特征凍結溫度最高的顆粒物觸發(fā)(Levine,1950;Sear,2013)。

    基于上述理論假設,溫度為T時單個液滴的平均冰核個數(shù)( μ(T)),可表示為

    式中,Vliquid是 單個液滴的體積,T是溫度,Cinp(T)是懸濁液樣本在溫度為T時的冰核數(shù)濃度。

    溫度為T時液滴凍結的概率可用泊松分布(Hoel,1962)表示

    式中,P[x,μ(T)]是群滴中x個液滴發(fā)生凍結的概率。因此,群滴中未凍結的比率P[0,μ(T)],可簡化為

    未凍結的比率(funfrozen)可用測量的量表示為

    式中,N(T)是在過冷溫度T時凍結的液滴個數(shù),N0為總液滴數(shù)。由式(1)、(3)和(4)可計算得到冰核數(shù)濃度為

    式中,Cinp(T)為懸濁液樣本在溫度為T時的冰核積分數(shù)濃度。

    通過式(5)可得出,液體中冰核數(shù)濃度能夠通過統(tǒng)計大量等體積小液滴的凍結隨溫度的變化來計算得到。大氣中冰核數(shù)濃度也可通過采集一定體積的大氣氣溶膠到濾膜,然后把濾膜放置在一定體積的超純水中涮洗混勻配置成氣溶膠懸濁液,通過冷臺凍結實驗得到大氣中冰核數(shù)濃度計算公式如下

    式中,Vwater是用來涮洗濾膜的超純水體積,Vair是經(jīng)過濾膜的大氣采樣體積。根據(jù)式(5)和(6),增加液滴總數(shù)或者提高液滴體積能夠提高測量下限,通過減小液滴體積能夠提高測量上限,也可用超純水稀釋的方法拓寬測量的溫度區(qū)間。

    FINDA測量是基于統(tǒng)計分析法得到,因此需要評估結果的可靠性。對于群數(shù)滴小于40的測量樣本,可通過威爾遜區(qū)間計算95%置信區(qū)間(Brown,et al,2001)。根據(jù)Agresti等(1998)提出的二項式分布法得到溫度為T時的95%置信區(qū)間CI95%(T)(Levin,et al,2019),表示為

    式中,N(T)是 溫度T時液滴凍結的個數(shù),N0是液滴總數(shù)。

    3 裝置搭建與校準

    3.1 FINDA硬件搭建

    FINDA基于下沉托盤式冰核測量冷臺(畢凱等,2020b)升級改進而來。圖1為FINDA裝置的結構,主要包括溫度控制系統(tǒng)、一體式96孔冷臺、圖像采集系統(tǒng)、溫度測量系統(tǒng)、氮氣供氣系統(tǒng)、封閉式設備箱與高亮LED光源、控制與數(shù)據(jù)分析系統(tǒng)等。其中,溫度控制系統(tǒng)主體為高性能水浴循環(huán)制冷機(FP-50 HL,Julabo,德國),可通過編程在200—-50℃精確控溫(±0.01℃)。圖像采集系統(tǒng)為500萬像素的USB面陣掃描相機(GO-5100MUSB,JAI,日本)搭配變焦鏡頭(TV ZOOM LENS J6X11,佳能,日本),通過萬向支架固定在冷臺中央上方。黑色設備保護箱和雙高亮LED燈的設計提供了穩(wěn)定的光環(huán)境供相機識別凍結。一體式96孔冷臺根據(jù)實驗使用的聚丙烯酶標板(Polypropylene polymerase chain reaction Tray,簡稱PCR板)尺寸使用太空鋁材定制加工而成(圖2),可緊密貼合PCR板以減小熱傳導造成的溫差。本冷臺與Stopelli等(2014)、David等(2019)和畢凱等(2020b)設計的PCR板浸沒制冷劑中的方式相比,有效隔離了制冷液蒸發(fā)與冷凝對樣品識別造成的干擾。冷臺浸入制冷液部分和上部臺面之間使用聚四氟乙烯材料進行了隔離密封設計,避免了制冷傳導過快造成的空氣在面板上部邊緣冷凝對相機識別凍結造成的干擾。實驗之前,冷臺以5 l/min通氮氣20 min,最大程度上排除環(huán)境干擾。實驗過程中以1 l/min的流量供氮氣,去除密閉冷臺腔中空氣可能含有的雜質對樣品的污染。

    圖1 FINDA結構Fig.1 Freezing Ice Nucleation Detector Array (FINDA)

    FINDA溫度測量系統(tǒng)由1臺多通道無紙記錄儀(6000C,杭州美控自動化技術有限公司)和4個經(jīng)過改裝的高精度Pt100溫度傳感器(不確定性δ1為±0.1℃)組成(圖2),Pt100鉑電阻片直接改裝到了單孔PCR管中并固定在冷臺4個角測溫,溫度傳感器改裝后經(jīng)過了校準(3.2節(jié))。目前國際上大多使用測量PCR板底部冷臺溫度作為液滴溫度的方法(Beall,et al,2017),而David等(2019)把PCR板浸入制冷劑中對比了板孔中液滴溫度和制冷劑溫度,發(fā)現(xiàn)在0—-25℃有1—3℃的溫差。因此,由于PCR板材的熱傳導差異,冷臺溫度難以代表PCR樣品中液滴的真實溫度,而FINDA經(jīng)過改裝后的Pt100由于同樣位于PCR管中,更能真實代表樣品中液滴的溫度。

    圖2 一體式96孔冷臺結構 (a. 上斜視圖,b. 下斜視圖)Fig.2 Integrated 96 hole cold stage (a. top view,b. bottom view)

    為了盡可能排除環(huán)境對樣品的污染,樣品制備過程均在層流式超凈臺(Heraguard ECO 1.2,Thermo Fisher Scientific,美國)上進行(超凈臺空氣中顆粒物含量為0.67±1.72 cm-3),樣品制備所需的超純水由Milli-Q (IQ 7000,Milli-Q,德國)超純水系統(tǒng)供應(水中顆粒物含量<1 cm-3),超純水在使用前通過0.02 μm孔徑的針式過濾器(Puradisc 25 syringe filters,Whatman,英國)過濾。懸濁液經(jīng)過超聲震蕩清洗5 min,然后在漩渦混合器中以2800 轉/min的速度混合5 min,充分混合均勻。FINDA裝置使用了4個角的孔作為溫度測量位置,因此PCR板需要裁去4個角的孔。包括前期的樣品處理在內,F(xiàn)INDA每次試驗過程大概需要1 h。FINDA使用單通道高精度電動移液器(50—1000 μl,Eppendorf Xplorer plus,Eppendorf,德國)取樣50 μl (等效體積直徑2300 μm)液滴(誤差±6%,精確度±1%),每次測量時把2組不同稀釋倍數(shù)的多個等體積(通常32個50 μl 液滴)液滴放入一次性無菌96孔PCR板中,同時放入28個50 μl 超純水滴作為背景對照,以剔除操作過程和超純水中可能帶來的污染。PCR板放置到冷臺后,整個設備以1℃/min的恒定速率降溫,直到所有滴都凍結。整個實驗過程通過Labview軟件程序進行自動控制并保存數(shù)據(jù)。

    3.2 測量控制與數(shù)據(jù)分析系統(tǒng)

    測量控制與數(shù)據(jù)分析系統(tǒng)由基于Labview編寫的實時分析系統(tǒng)和基于Python編寫的數(shù)據(jù)解析軟件(畢凱等,2020c)組成。其中Labview實時分析系統(tǒng)根據(jù)液滴凍結前后光學特性的變化對各個液滴凍結進行監(jiān)測與顯示,同時以1 s的分辨率把4個溫度傳感器的數(shù)據(jù)和96孔的亮度值保存在“.csv”格式的數(shù)據(jù)文件中。軟件系統(tǒng)主要界面和算法設計流程如圖3。

    樣品放置到冷臺后啟動自動降溫過程,制冷機首先以最大功率從室溫降到0℃,維持10 min達到熱平衡。然后以1℃/min的降溫率進行降溫,直到-30℃。軟件基于用戶設置的背景參考圖片和液滴識別區(qū)參數(shù)提取每個液滴的位置(8×12位置陣列),液滴識別區(qū)通常設置為75×75像素,每個液滴的平均亮度灰度值被記錄下來。閾值設置完畢后啟動識別,以1 Hz的速度提取視頻中所有液滴位置的光亮度信息,并保存下來,與DRINCZ (15 s提取1張圖片;David,et al,2019)和IR-NIPI (20 s提取1張圖片;Harrison,et al,2018)等測量設備相比,采樣率大幅度提高,凍結溫度識別更精確。在監(jiān)控頁面,用戶可以實時查看亮度信息、溫度信息以及凍結的位置(圖3)。

    圖3 軟件界面及凍結識別算法流程 (a. 參數(shù)設置界面;b. 實時監(jiān)控界面,紅色為凍結的液滴位置;c. 軟件設計流程)Fig.3 Screenshot of FINDA software and Flow chart showing the detection of frozen droplet(a. parameters specification;b. real-time show,the red color stands for the frozen droplet;c. flow chart of software)

    圖4給出了一次降溫實驗中液滴凍結識別效果,由于液滴凍結前、后反射光亮度的差異被相機捕獲識別(圖4b),所有液滴位置光亮度在實驗過程中都經(jīng)過了亮度突變衰減(圖4a),亮度像素遞減率最大的位置對應著冰核核化事件(圖4c)。液滴凍結的過程對應著光亮度灰度值先增大后緩慢減小直到完全凍結后趨于穩(wěn)定(圖4d)。從圖4a也可以看出,不同液滴的亮度灰度由于位置差異存在系統(tǒng)性差異,但對每個液滴來說,凍結發(fā)生前和完全凍結后光亮度灰度值非常穩(wěn)定,這反映了測量系統(tǒng)的光環(huán)境穩(wěn)定性。液滴凍結不是瞬時完成的,每個液滴從冰核核化到完全凍結需要的時間差異較大,凍結的速度隨溫度降低而加快(圖4a),這可能與低溫下凍結釋放的潛熱熱量更容易被周圍環(huán)境吸收有關。

    圖4 基于反射光亮度差異的液滴凍結識別效果 (a. 92個液滴平均亮度灰度值隨溫度的變化;b. 凍結監(jiān)控界面截圖,圓圈中白色為未凍結液滴,圓圈中灰黑色為已凍結液滴;c. 92個液滴像素遞減率隨溫度的變化;d. 圖b中G8孔 (紅色圓圈) 的亮度像素隨溫度的變化)Fig.4 Frozen detection result based on differences in pixels of different droplets (a. variations of pixel values of 92 droplets with temperature;b. capture of software that shows unfrozen droplets by white circles and frozen droplets by dark circles;c. variations of decreasing rate of pixels for 92 droplets with temperature;d. pixel values of the G8 hole (red circle) with temperature)

    軟件系統(tǒng)中Labview程序僅用于記錄液滴的亮度信息和初步識別。冰核數(shù)濃度的反演和質量控制使用Python代碼實現(xiàn)。Python程序讀取Labview保存的數(shù)據(jù),進行識別分析、計算、質量控制以及繪圖等。由于同時保存了原始的亮度信號和液滴監(jiān)控視頻,因此可通過后處理程序對實驗過程進行全流程跟蹤。

    3.3 溫度測量與校準

    FINDA冷臺溫度為4個角孔的溫度測量平均值。為了提高測量精度,使用經(jīng)過中國計量科學研究院校準的標準溫度傳感器(不確定性為±0.008℃)在恒溫酒精浴槽(FP-50 HL,Julabo,德國)環(huán)境中每隔5℃穩(wěn)定10 min對4個改裝的Pt100傳感器進行標定(0—-30℃),并進行線性擬合校準(表1)。校準后測量值與平均值偏差小于±0.05℃(圖5),在Pt100(不確定性 δ1為±0.1℃)測量誤差范圍內。

    圖5 Pt100 溫度傳感器校準后溫度與4個傳感器平均溫度的偏差Fig.5 Deviations of calibrated Pt100 sensors from mean temperatures of 4 sensors

    表1 溫度擬合參數(shù)Tcal=aTraw+bTable 1 Parameters for temperature fitting with the functional form Tcal=aTraw+b

    由于冷臺及PCR板的熱傳導差異,冷臺不同位置孔中存在溫差。為了準確評估FINDA測量不確定性,使用高精度(640×480像素,0.04℃分辨率)視頻式( 30Hz 幀率)紅外熱像儀(T660,F(xiàn)LIR,瑞典)對冷臺面板溫度差進行測量與標定。在PCR板中放入50 μl 無水乙醇替代液體樣品,冷臺降溫至0℃穩(wěn)定10 min,然后以1℃/min降溫速率降至-30℃。提取96個孔中溫度值,對每個孔測量數(shù)據(jù)進行30 s平均(900個溫度值,對應約0.5℃)。

    圖6為四個角孔平均溫度為-10、-15和-20℃時冷臺面板的溫度分布,實驗樣品存在系統(tǒng)性水平溫差,冷臺中部溫度低于邊緣位置,溫度最高區(qū)域在右上角(圖6)。分布趨勢與Beall等(2017)利用有限元分析法模擬的AIS(Automated ice Spectrometer)冷臺一致。

    圖6 1℃/min降溫率下-10℃(a)、-15℃ (b) 和-20℃ (c) 時FINDA冷臺96孔中樣本溫度分布Fig.6 Temperature of 96 holes of FINDA cold stage at -10℃ (a),-15℃ (b) and -20℃ (c) with cooling rate of 1℃/min

    圖7為1℃/min降溫率下紅外熱像儀測量的96個孔的溫度與其中4個角孔位置的平均溫度之差的統(tǒng)計。圖7中96個孔的溫度與4個角位置溫度差(圖7a)的平均值和標準差(圖7b)均隨溫度降低而呈明顯的線性偏差,面板實際溫度比測量溫度低,在-25℃時偏差為-0.6℃(標準差為±0.3℃)。綜合溫度傳感器的誤差,F(xiàn)INDA溫度測量不確定性在0—-25℃測量范圍內約為±0.35—±0.75℃,國際上其他同類型設備不確定性通常為±0.9℃(David,et al,2019)、±0.9℃(Harrison,et al,2018)和±2.2℃(Beall,et al,2017)。針對冷臺水平溫差,采用兩種校準方法進行了對比分析。第一種采用整體校準法,利用冷臺的線性偏差對溫度測量結果進行線性擬合(表2)校正(圖9中橙、藍、綠、紅色圓點)。第二種采用單孔校準法,根據(jù)面板溫差對每個溫度下單孔的實際溫度進行了擬合校準(圖9黑色三角)。圖9表明,這兩種校準方法的測量結果無明顯差異。FINDA測量結果通常使用面板整體溫差動態(tài)校準,并添加溫度誤差(包括傳感器誤差和冷臺溫差誤差)。

    圖7 1℃/min降溫率下冷臺96孔溫度與4個角孔平均溫度偏差 (a) 及標準差 (b) 統(tǒng)計及擬合 (黑色直線)Fig.7 Statistics and fitting of temperature bias (a) and deviation (b) of the cold stage 96 holes and average temperature of four corners with cooling rate of 1℃/min(black line represents the fitting line)

    表2 樣品水平溫差及標準差擬合參數(shù)y=ax+bTable 2 Parameters for the fitting of temperature bias and deviation with the functional form y=ax+b

    4 設備測試與評估

    4.1 超純水凍結實驗

    液滴凍結法測量冰核通常使用超純水稀釋來拓寬溫度檢測區(qū)間(Vali,1971),大氣冰核檢測時采樣濾膜也需用超純水配置懸濁液,因此需要評估超純水本身含有的以及實驗過程引入的冰核對液滴凍結的影響。對10組(各28個50 μl 液滴)不同的超純水樣品進行了凍結實驗。Milli-Q超純水使用前通過0.02 μm過濾。樣品分別放置到無菌潔凈PCR板孔中在冷臺上以1 ℃/min的降溫率進行凍結實驗。

    超純水凍結溫度取決于液滴體積和降溫速率,同一超純水樣品1 μl 液滴凍結溫度(圖8藍色點)整體上比50 μl 液滴(圖8藍色實線)溫度低。對于大體積液滴(幾十微升),難以去除超純水和液滴接觸介質中所有雜質冰核,因此凍結溫度通常高于同質核化溫度(Koop,et al,2016)。經(jīng)典核化理論計算的50 μl 液滴的同質核化凍結中值溫度(T50,定義為50%的液滴個數(shù)凍結的溫度)約為-33℃(Harrison,et al,2018)。FINDA超純水T50平均為-25.6±0.64℃,初始凍結溫度平均為-21.27±1.45℃,全部凍結溫度平均為-26.93±0.63℃。凍結比率(ff)隨溫度的變化擬合如下

    式中,T單位為℃。

    盡管FINDA超純水凍結溫度與理論計算仍存在一定的差異,但整體上凍結溫度低于其他同樣使用50 μl 液滴的DRICZ (David,et al,2019)和IR_NIPI (Harrison,et al,2018)等設備(圖8綠色和橙色實線),也低于中國已有的10 μl 水滴凍結裝置(表3)。超純水凍結實驗結果與使用的水源、制備樣品的操作環(huán)境、冷臺所處的環(huán)境等因素有關。FINDA有效測量范圍0—-25℃內,尤其是-20℃以上的溫度下,測量結果基本不受超純水本身及實驗環(huán)境引入冰核的影響。為了進一步剔除實驗背景中的冰核影響,測量時放置28個50 μl 液滴進行平行對照實驗,并在數(shù)據(jù)分析時剔除背景冰核對樣品的影響。

    表3 FINDA與其他凍結設備 (使用幾十微升液滴) 超純水凍結試驗對比Table 3 Comparison of freezing temperature of the ultra-pure water experiment of FINDA with results using other instruments with tens μl droplet

    圖8 FINDA超純水凍結比率曲線及與其他設備對比(黑色實線為10組不同的超純水凍結結果,紅色實線為10組實驗的擬合曲線,藍色實線和點線分別為同一樣品50和1 μl液滴的凍結曲線,綠實線為DRINCZ設備 (David,et al,2019) 的擬合曲線,橙色實線為IR_NIPI設備 (Harrison,et al,2018)超純水凍結實驗的擬合曲線)Fig.8 Fraction of droplets frozen as a function of temperature of FINDA and comparison with results of other instruments (black solid line:10 groups of different ultra-pure water,red solid line:fitting line of 10 groups droplets freezing experiment of FINDA,blue solid line and dot line are 50 and 1 μl droplet experiments for the same sample,green solid line:fitting line of DRINCA instrument (David,et al,2019),orange solid line:fitting line of data from water experiment of IR_NIPI instrument (Harrison et al,2018))

    4.2 降雪中冰核檢測

    為了評估FINDA裝置對液體樣品中冰核的檢測性能,采集了2020年11月21日海坨山云物理綜合觀測站(畢凱等,2020a)和2021年1月19日北京市氣象局大院的降雪樣品進行冰核檢測。聚丙烯材質的雪采樣盒使用前經(jīng)過超純水清洗,雪樣采集完畢后放置到50 ml滅菌離心管(康寧,美國)內,雪樣的存儲和運輸均在-20℃環(huán)境下,以盡可能減少由于細菌生長或液體存儲造成的測量誤差(Stopelli,et al,2014)。樣品分析前在室溫下解凍,然后在漩渦混合器中以2800轉/min的速度混合5 min,為了拓寬測量溫度區(qū)間,同時進行了15、225和3375倍濃度稀釋。FINDA每次檢測2組32個樣品滴,同時使用28個超純水滴作為測量背景,每隔0.25℃冰核積分數(shù)濃度根據(jù)式(5)進行計算。其中11月21日樣品使用北京大學的顯微平板冷臺PKU-INA(Chen,et al,2018)進行檢測對比,PKU-INA冷臺試驗使用90個1 μl的液滴。

    圖9為11月21日10時樣品的檢測結果,原始濃度降雪樣品的凍結溫度在-8—-15℃(圖9紅色圓圈),通過稀釋方法檢測溫度可達到-25℃。由于樣品中冰核數(shù)濃度高出背景超純水3個數(shù)量級以上(圖9黑色點),可以忽略超純水和實驗過程引入的冰核對測量結果的影響。使用整體面板溫度偏差校準和單孔溫度偏差校準(圖9黑色圓圈)的測量結果無明顯差異。與PKU-INA測量結果(圖9粉色圓圈和淺藍色圓圈)相比,冰核譜型基本一致,偏差小于3℃。這個偏差小于Hiranuma等(2015)給出的使用6種不同的冷臺凍滴設備對同種物質測量得到的5℃偏差結果。因此,F(xiàn)INDA測量結果較為合理,且由于試驗液滴體積(50 μl )大于PKU-INA(1 μl ),能夠有效檢測降水中高核化溫度(高于-15℃)的冰核。

    圖9 2020年11月21日北京山區(qū)降雪中的冰核濃度與對比測量 (黑色實心點為實驗背景中冰核,紅色為原始濃度,綠色為15倍稀釋,藍色為225倍稀釋,橙色為3375倍稀釋,黑色圓圈為單孔溫度校準的結果;粉色為PKU-INA原始濃度測量結果,淺藍色為PKU-INA 225倍稀釋結果)Fig.9 Measurement and comparison of INPs concentration in snow water during the snow case on 21 November 2020 (black dots are background INPs,red circles are INPs in original snow water,green circles are INPs in 15 fold dilution sample,blue circles are INPs in 225 fold dilution sample,orange circles are INPs in 3375 fold dilution sample,pink circles are INPs of original snow water measured by PKU-INA,light blue circles are INPs in 225 fold dilution sample by PKU-INA)

    通過對比分析兩次降雪過程的雪水冰核譜(圖10a),可以看出,整體上原始濃度雪水的冰核譜均在Petters等(2015)總結的其他地區(qū)的冰核譜范圍內,與酆大雄等(2002)測量的北京地區(qū)4次降雪中冰核濃度范圍有一定的差異。FINDA測量的兩次降雪過程冰核譜差異較大,1月19日降雪過程雪水中含有較多高活化溫度的冰核。與11月21日小尺寸柱狀雪晶不同,1月19日雪花以枝狀大尺寸雪晶為主(圖10b和c),下落過程中有更多機會捕獲高活化溫度的冰核,另外也與不同的降雪氣團來源不同有一定的關系(David,et al,2019)。FINDA通過稀釋法可以拓寬檢測溫度達到-25℃,David等(2019)同樣通過雪水中冰核檢測對比驗證了稀釋法的可靠性。

    圖10 兩次降雪過程中冰核譜對比 (a. 紅色圓圈、橙色三角、粉色圓圈和紅色圓點為2021年1月19日樣品,藍色圓圈、綠色三角、淺藍色圓圈和藍色圓點為2020年11月21日樣品,灰色陰影為Petters等(2015)中匯總結果,黃色陰影為酆大雄等(2002)測量的北京原始濃度雪水中冰核范圍;b. 2020年11月21日降雪過程的柱狀冰晶;c. 2021年1月19日降雪過程的枝狀雪花)Fig.10 INPs concentration spectra for two snow cases (a. red circles,orange triangles,pink circles and red dots are for original,15 fold,225 fold and 3375 fold dilution snow samples for the case on 19 January 2021, the blue circles,green triangles,light blue circles and blue dots are for original,15 fold,225 fold and 3375 fold dilution snow samples for the case on 21 November 2020, the area shaded in gray shows results from Petters,et al (2015), the area shaded in yellow shows the fitting result of Feng (2002);b. column shape snowflakes on 21 November 2020;c. branch shape snowflakes on 19 January 2021)

    4.3 大氣中冰核檢測

    為了評估FINDA設備對大氣中冰核樣品的檢測性能,2021年2月16日在海坨山云物理綜合觀測站(畢凱等,2020a)使用0.2 μm孔膜(Whatman,英國)進行了大氣氣溶膠采樣。采樣器及膜前期處理方法與Schiebel (2017)一致。本次采樣共使用2張膜,其中1張膜以8 l/min的流量采集了2 h共960 l大氣氣溶膠樣本,另1張膜作為背景不進行采樣,以評估膜的處理、運輸、存儲等環(huán)節(jié)可能引入的冰核污染,同時使用BJ-CFDC穩(wěn)定在-20℃時浸潤模式(畢凱等,2020a)同步測量的冰核濃度進行對比。冰核檢測時把濾膜放置到10 ml超純水中,經(jīng)過超聲震蕩清洗5 min,然后在漩渦混合器中以2800轉/min的速度混合5 min,配置懸濁液,同時進行了15倍和225倍稀釋,每種樣品使用32個滴。冰核數(shù)濃度根據(jù)式(6)進行計算。

    通過不同濃度稀釋法,F(xiàn)INDA能夠檢測的冰核數(shù)濃度跨度4個數(shù)量級,核化溫度區(qū)間跨度20℃(圖11),這覆蓋了大部分混合云的溫度區(qū)間。與空白膜相比(圖11灰色圓圈),采樣膜同溫度下冰核數(shù)濃度高1個數(shù)量級,說明樣品制備、運輸、存儲和分析等環(huán)節(jié)引入冰核對檢測結果的影響較小。與BJ-CFDC在-20℃時測量結果相比,F(xiàn)INDA測量結果略偏高,這可能是由于BJ-CFDC前端使用了2.5 μm的粒子篩選器的原因,這與DeMott等(2017)在外場觀測中對比CFDC和多種離線凍滴設備的結果基本一致。盡管存在一定的誤差,但仍然低于大氣冰核本身測量的不確定性(DeMott,et al,2017)。FINDA通過大采樣體積累積可有效檢測低濃度(小于0.1 l-1)和生物冰核等高核化溫度(高于-15℃)冰核數(shù)濃度。

    圖11 2021年2月16日海坨山大氣中冰核濃度譜(黑色圓圈為空白濾膜背景冰核,紅色為原始濃度,綠色為15倍稀釋,藍色為225倍稀釋,紫色為BJ-CFDC穩(wěn)定在-20℃時同步測量結果)Fig.11 INPs concentration spectrum in the atmosphere on 16 February 2021 (black circles indicate background INPs concentration of the blank filter;red,green and blue circles are for original,15 fold dilution and 225 fold dilution samples;purple circles represent the BJ-CFDC result at -20℃at the same time of the filter sample)

    5 結論與討論

    介紹了一種新型用于檢測浸潤凍結機制冰核裝置(FINDA)的搭建、質量控制及初步應用,F(xiàn)INDA通過Labview軟件控制冷臺以1℃/min速度降溫,同時利用CCD相機實時監(jiān)測冷臺上96孔PCR板孔中50 μl液滴的反射光亮度以自動識別凍結事件,并計算得到0—-25℃間隔0.25℃的冰核濃度。通過進行超純水的凍結溫度檢驗,以及對降雪樣品和大氣采樣中冰核的檢測對比,評估了該裝置的性能,得到以下結論:

    (1)FINDA能 夠 有 效 測 量0—-25℃間隔0.25℃的冰核濃度。

    (2)FINDA溫度不確定性隨溫度降低而呈線性增大,在-25℃時最大為±0.75℃,與其他同類型的凍滴測量裝置相似。

    (3)FINDA的超純水初始凍結溫度平均為-21.27±1.45℃,中值凍結溫度平均為-25.6±0.64℃,表明實驗中由于超純水和樣品制備過程引入的冰核對測量結果的影響不大。

    (4)不同降雪過程雪水的冰核譜有所不同,同一樣品FINDA測量結果與PKU-INA(1 μl液滴)凍滴設備趨勢一致,偏差小于3℃,這與不同設備的測量方式有關,且偏差小于其他學者研究的多設備對比結果。

    (5)通過不同倍數(shù)稀釋法可獲得寬溫度區(qū)間大氣冰核濃度譜,與同步測量的在線連續(xù)流擴散云室(BJ-CFDC)相比,F(xiàn)INDA測量結果略偏高,這與BJ-CFDC使用了粒子篩選器有關,與國際上其他學者的研究結果一致。說明在大氣冰核測量中也有很大的優(yōu)勢。

    綜上所述,文中介紹的裝置可以有效測量大氣、降水或其他樣品溶液中的冰核數(shù)濃度譜和凍結溫度比率,測量結果較可靠。對于生物氣溶膠中的高活化溫度冰核檢測更有效。受液滴體積及超純水背景的影響,F(xiàn)INDA無法有效測量低于-25℃的冰核。但這部分可通過在線連續(xù)流擴散云室進行測量,F(xiàn)INDA可彌補在線連續(xù)流擴散云室(BJCFDC)在-15℃以上的測量缺陷,通過兩者協(xié)同觀測,獲得大氣中寬溫譜的冰核濃度譜。由于FINDA采用離線檢測,因此未來可以通過在中國不同地點采樣,在同一標準下獲得冰核的空間分布特征,幫助深入了解冰核分布特征,為預報和云模式的本地化改進及對云降水理論研究提供實驗數(shù)據(jù)和理論基礎。

    致 謝:感謝中國氣象局華北云降水野外科學試驗基地提供所有觀測資料。

    猜你喜歡
    超純水液滴測量
    Synthesis of new non-fluorous 2,2'-bipyridine-4,4'-dicarboxylic acid esters and their applications for metal ions extraction in supercritical carbon dioxide
    小水量超純水制備系統(tǒng)的最佳工藝選擇
    液滴間相互碰撞融合與破碎的實驗研究
    噴淋液滴在空氣環(huán)境下的運動特性
    把握四個“三” 測量變簡單
    滑動摩擦力的測量和計算
    滑動摩擦力的測量與計算
    測量
    雜銅冶煉廠超純水管路系統(tǒng)的設計
    氣井多液滴攜液理論模型研究
    河南科技(2013年10期)2013-08-12 05:55:10
    一本综合久久免费| 亚洲第一欧美日韩一区二区三区| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久人妻精品电影| 国产精品香港三级国产av潘金莲| 在线观看www视频免费| 亚洲av熟女| 免费女性裸体啪啪无遮挡网站| 成人午夜高清在线视频 | 可以在线观看毛片的网站| 中文字幕高清在线视频| 脱女人内裤的视频| 亚洲九九香蕉| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品色激情综合| 丁香六月欧美| 亚洲国产精品成人综合色| 麻豆成人午夜福利视频| 天天一区二区日本电影三级| 国产欧美日韩一区二区精品| 国产精品免费一区二区三区在线| 欧美性长视频在线观看| 国产97色在线日韩免费| 亚洲av日韩精品久久久久久密| 丰满的人妻完整版| 黄色女人牲交| 国产精品免费一区二区三区在线| 亚洲av日韩精品久久久久久密| 亚洲av日韩精品久久久久久密| 俄罗斯特黄特色一大片| 亚洲av电影在线进入| 黄色视频,在线免费观看| 久久精品国产亚洲av高清一级| 狠狠狠狠99中文字幕| 欧美日韩乱码在线| 91成年电影在线观看| 久久香蕉国产精品| 国产三级黄色录像| 久久亚洲真实| 久久午夜亚洲精品久久| 精品一区二区三区四区五区乱码| 成人国语在线视频| 亚洲天堂国产精品一区在线| 波多野结衣巨乳人妻| 女性被躁到高潮视频| 麻豆成人av在线观看| 亚洲成人久久性| 久久久久亚洲av毛片大全| 国产99白浆流出| 免费看a级黄色片| 黄色女人牲交| 久久精品影院6| 亚洲精品中文字幕一二三四区| 听说在线观看完整版免费高清| av超薄肉色丝袜交足视频| 精品一区二区三区四区五区乱码| 亚洲精品国产一区二区精华液| 正在播放国产对白刺激| 国产亚洲精品av在线| 日韩中文字幕欧美一区二区| 亚洲真实伦在线观看| 精品国产乱子伦一区二区三区| 成人18禁高潮啪啪吃奶动态图| 国产亚洲精品一区二区www| 免费在线观看影片大全网站| 少妇被粗大的猛进出69影院| 亚洲欧美一区二区三区黑人| 国产97色在线日韩免费| 久久久久久久久免费视频了| 一二三四在线观看免费中文在| 亚洲国产毛片av蜜桃av| 日韩视频一区二区在线观看| 欧美乱妇无乱码| 午夜老司机福利片| 1024手机看黄色片| 中文字幕人妻熟女乱码| 18禁黄网站禁片午夜丰满| 欧美丝袜亚洲另类 | svipshipincom国产片| 淫妇啪啪啪对白视频| 黑人操中国人逼视频| 欧美精品亚洲一区二区| 日本成人三级电影网站| www日本在线高清视频| 手机成人av网站| 99久久久亚洲精品蜜臀av| ponron亚洲| avwww免费| 国产精品,欧美在线| 哪里可以看免费的av片| 亚洲九九香蕉| 亚洲五月天丁香| 午夜精品在线福利| 最近在线观看免费完整版| 欧美黄色淫秽网站| 成人亚洲精品av一区二区| 欧美中文综合在线视频| 亚洲精品av麻豆狂野| 欧美乱色亚洲激情| 又紧又爽又黄一区二区| 国产成人欧美| 亚洲自拍偷在线| 国产午夜精品久久久久久| 久久中文字幕一级| 欧美成人性av电影在线观看| 亚洲成人久久爱视频| 国产成人av教育| 免费在线观看成人毛片| 国产亚洲欧美精品永久| 亚洲男人的天堂狠狠| 啦啦啦观看免费观看视频高清| netflix在线观看网站| 在线十欧美十亚洲十日本专区| 少妇熟女aⅴ在线视频| 国产精品98久久久久久宅男小说| 自线自在国产av| 国产视频一区二区在线看| 亚洲熟女毛片儿| 999久久久精品免费观看国产| 国产一区二区三区视频了| 国产精品99久久99久久久不卡| www日本在线高清视频| 黄色毛片三级朝国网站| 日本成人三级电影网站| 亚洲国产精品sss在线观看| www.999成人在线观看| av有码第一页| 高清毛片免费观看视频网站| 亚洲人成伊人成综合网2020| 精品一区二区三区av网在线观看| 曰老女人黄片| 大型黄色视频在线免费观看| 1024香蕉在线观看| 一二三四社区在线视频社区8| 大型av网站在线播放| 久久久国产成人精品二区| 黑人欧美特级aaaaaa片| 国产亚洲精品av在线| 夜夜看夜夜爽夜夜摸| 国产亚洲精品第一综合不卡| 亚洲第一青青草原| 亚洲精品粉嫩美女一区| 在线观看午夜福利视频| 可以免费在线观看a视频的电影网站| 亚洲 欧美一区二区三区| 午夜两性在线视频| 久久精品国产亚洲av香蕉五月| 1024香蕉在线观看| 亚洲av成人一区二区三| 欧美不卡视频在线免费观看 | 成人永久免费在线观看视频| 黄色 视频免费看| 久久精品国产亚洲av香蕉五月| 香蕉丝袜av| 午夜老司机福利片| 天堂影院成人在线观看| 日韩精品青青久久久久久| 免费女性裸体啪啪无遮挡网站| 欧美中文综合在线视频| 久久久久亚洲av毛片大全| 成人午夜高清在线视频 | 琪琪午夜伦伦电影理论片6080| 国产精品九九99| 在线观看66精品国产| 亚洲午夜精品一区,二区,三区| 久久精品91蜜桃| 成人三级黄色视频| 波多野结衣高清作品| 18美女黄网站色大片免费观看| 欧美人与性动交α欧美精品济南到| 在线观看一区二区三区| 90打野战视频偷拍视频| 男人操女人黄网站| 美女 人体艺术 gogo| 亚洲人成电影免费在线| 特大巨黑吊av在线直播 | 久久婷婷成人综合色麻豆| 两个人免费观看高清视频| 亚洲av第一区精品v没综合| 日本撒尿小便嘘嘘汇集6| 啦啦啦韩国在线观看视频| 999久久久国产精品视频| 中文字幕另类日韩欧美亚洲嫩草| 精品久久久久久久久久免费视频| 精品国产一区二区三区四区第35| 国产精品美女特级片免费视频播放器 | 亚洲真实伦在线观看| 一级片免费观看大全| 国内精品久久久久精免费| 男女午夜视频在线观看| 国产亚洲av嫩草精品影院| 免费在线观看日本一区| 两个人免费观看高清视频| 国产极品粉嫩免费观看在线| 国产乱人伦免费视频| 哪里可以看免费的av片| 丰满人妻熟妇乱又伦精品不卡| 国内少妇人妻偷人精品xxx网站 | 日本免费a在线| 日韩欧美国产在线观看| 国产亚洲精品第一综合不卡| 1024手机看黄色片| 国产乱人伦免费视频| 国产一区二区激情短视频| 99久久综合精品五月天人人| 国产高清视频在线播放一区| 欧美成人一区二区免费高清观看 | 欧美久久黑人一区二区| 精华霜和精华液先用哪个| 一本一本综合久久| 免费无遮挡裸体视频| 99在线人妻在线中文字幕| 精品午夜福利视频在线观看一区| 亚洲熟女毛片儿| 极品教师在线免费播放| 怎么达到女性高潮| 久热这里只有精品99| 午夜免费激情av| 午夜久久久久精精品| 国产v大片淫在线免费观看| 1024视频免费在线观看| 免费在线观看亚洲国产| 中文字幕av电影在线播放| 每晚都被弄得嗷嗷叫到高潮| 99热只有精品国产| 99riav亚洲国产免费| 91大片在线观看| 亚洲一区中文字幕在线| 嫩草影院精品99| 高清在线国产一区| 欧美日本亚洲视频在线播放| 欧美成狂野欧美在线观看| 国产免费男女视频| av有码第一页| 久久久久久久久免费视频了| 草草在线视频免费看| 一进一出抽搐动态| 曰老女人黄片| 亚洲第一青青草原| 女性生殖器流出的白浆| 男女做爰动态图高潮gif福利片| 久久久精品欧美日韩精品| 麻豆成人午夜福利视频| 国产激情久久老熟女| 国产熟女xx| 久久人人精品亚洲av| av在线播放免费不卡| 麻豆成人av在线观看| 亚洲专区国产一区二区| 亚洲美女黄片视频| a级毛片a级免费在线| 日韩大码丰满熟妇| 免费在线观看日本一区| 一边摸一边抽搐一进一小说| 很黄的视频免费| 国产精品电影一区二区三区| 亚洲国产欧美日韩在线播放| 色综合站精品国产| 黄色片一级片一级黄色片| 黑人操中国人逼视频| 国产成年人精品一区二区| 成人亚洲精品av一区二区| 此物有八面人人有两片| 午夜福利视频1000在线观看| 成人国产综合亚洲| 性色av乱码一区二区三区2| 日韩大尺度精品在线看网址| 宅男免费午夜| 欧美国产日韩亚洲一区| 一级a爱片免费观看的视频| 久99久视频精品免费| 精品少妇一区二区三区视频日本电影| 午夜亚洲福利在线播放| 91九色精品人成在线观看| 欧美黑人欧美精品刺激| 高清在线国产一区| 国产主播在线观看一区二区| 亚洲av电影不卡..在线观看| 亚洲国产精品999在线| 18禁黄网站禁片午夜丰满| 亚洲成人免费电影在线观看| 欧美成人一区二区免费高清观看 | 久久精品国产99精品国产亚洲性色| 悠悠久久av| 一级片免费观看大全| 国产伦一二天堂av在线观看| 淫秽高清视频在线观看| 国产激情欧美一区二区| 亚洲男人天堂网一区| 亚洲精品国产一区二区精华液| 亚洲真实伦在线观看| 1024手机看黄色片| 亚洲精品国产区一区二| 99久久无色码亚洲精品果冻| 欧美精品啪啪一区二区三区| 日韩欧美免费精品| 一本综合久久免费| 欧美日韩亚洲综合一区二区三区_| 我的亚洲天堂| 岛国视频午夜一区免费看| 成人亚洲精品av一区二区| 日韩免费av在线播放| 色精品久久人妻99蜜桃| 日韩精品中文字幕看吧| 51午夜福利影视在线观看| 热re99久久国产66热| 99热只有精品国产| 欧美久久黑人一区二区| 久久久久久久久久黄片| 久久午夜综合久久蜜桃| 少妇被粗大的猛进出69影院| 亚洲avbb在线观看| 中文字幕人成人乱码亚洲影| 国产国语露脸激情在线看| 中文资源天堂在线| 久久久久久国产a免费观看| 精品久久久久久久久久久久久 | 黑丝袜美女国产一区| 看免费av毛片| www.999成人在线观看| 一本精品99久久精品77| 午夜a级毛片| 欧美中文日本在线观看视频| 国产成人系列免费观看| 久久精品国产99精品国产亚洲性色| 波多野结衣高清无吗| 国产欧美日韩一区二区精品| 亚洲一区二区三区色噜噜| 亚洲精品久久国产高清桃花| 亚洲精品一卡2卡三卡4卡5卡| 欧美激情高清一区二区三区| 国产精品亚洲美女久久久| 一个人免费在线观看的高清视频| 亚洲av中文字字幕乱码综合 | 搡老熟女国产l中国老女人| 国产一区二区激情短视频| 高潮久久久久久久久久久不卡| 香蕉av资源在线| 日韩欧美在线二视频| 一a级毛片在线观看| 欧美中文日本在线观看视频| 波多野结衣av一区二区av| 精品一区二区三区四区五区乱码| 神马国产精品三级电影在线观看 | 一区二区三区国产精品乱码| 女性生殖器流出的白浆| 国产99久久九九免费精品| 日本精品一区二区三区蜜桃| xxxwww97欧美| 中文字幕av电影在线播放| 校园春色视频在线观看| 色综合婷婷激情| 男女下面进入的视频免费午夜 | 欧美人与性动交α欧美精品济南到| 国产精品电影一区二区三区| 国产真人三级小视频在线观看| 久久午夜亚洲精品久久| 欧美 亚洲 国产 日韩一| 久久精品影院6| 国产亚洲欧美精品永久| 好男人电影高清在线观看| 免费看十八禁软件| 国产成+人综合+亚洲专区| 亚洲成av片中文字幕在线观看| 国产亚洲欧美在线一区二区| 免费在线观看亚洲国产| 不卡一级毛片| 十八禁人妻一区二区| 真人做人爱边吃奶动态| 精品久久久久久,| 亚洲成av人片免费观看| 国产亚洲精品一区二区www| 1024香蕉在线观看| 成人亚洲精品av一区二区| 国产野战对白在线观看| 国产熟女午夜一区二区三区| 成年版毛片免费区| 久久中文字幕一级| 久久国产精品人妻蜜桃| 一区福利在线观看| 在线观看免费午夜福利视频| 熟女电影av网| 日韩视频一区二区在线观看| 亚洲欧美日韩高清在线视频| 国内揄拍国产精品人妻在线 | 一本综合久久免费| 国产三级黄色录像| 国产亚洲欧美精品永久| 国内少妇人妻偷人精品xxx网站 | 国产免费av片在线观看野外av| 亚洲免费av在线视频| 免费无遮挡裸体视频| 视频区欧美日本亚洲| 久久午夜综合久久蜜桃| 两个人免费观看高清视频| x7x7x7水蜜桃| 亚洲一区二区三区色噜噜| 成人免费观看视频高清| 欧美一级a爱片免费观看看 | 午夜免费观看网址| 亚洲欧美精品综合久久99| 成人国产一区最新在线观看| 黄片大片在线免费观看| 91在线观看av| 波多野结衣高清作品| 中出人妻视频一区二区| 欧美性长视频在线观看| 婷婷六月久久综合丁香| 性欧美人与动物交配| 日韩欧美 国产精品| 亚洲中文av在线| 很黄的视频免费| 国产精品1区2区在线观看.| 亚洲中文av在线| 人人妻人人澡人人看| 久久精品亚洲精品国产色婷小说| 日韩欧美三级三区| 少妇的丰满在线观看| 日本免费一区二区三区高清不卡| 最近最新中文字幕大全免费视频| 黄色a级毛片大全视频| 久久久久亚洲av毛片大全| 亚洲五月婷婷丁香| 欧美性长视频在线观看| 国产av不卡久久| 国产人伦9x9x在线观看| 午夜影院日韩av| 欧美日韩一级在线毛片| 亚洲五月色婷婷综合| 精品一区二区三区视频在线观看免费| 精品久久久久久,| 1024视频免费在线观看| 一级毛片女人18水好多| 老熟妇仑乱视频hdxx| 亚洲性夜色夜夜综合| 伦理电影免费视频| 成年版毛片免费区| 亚洲一区中文字幕在线| 久久久久久久久久黄片| 欧美 亚洲 国产 日韩一| 国内毛片毛片毛片毛片毛片| 久久人妻av系列| 欧美黑人精品巨大| 亚洲欧美激情综合另类| 丝袜在线中文字幕| 欧美色视频一区免费| 午夜视频精品福利| 非洲黑人性xxxx精品又粗又长| 国产亚洲精品第一综合不卡| 亚洲男人的天堂狠狠| 一进一出好大好爽视频| 国产激情欧美一区二区| 91av网站免费观看| 午夜激情av网站| 桃红色精品国产亚洲av| 亚洲全国av大片| xxxwww97欧美| 欧美黑人精品巨大| 我的亚洲天堂| 成在线人永久免费视频| 白带黄色成豆腐渣| 亚洲成av人片免费观看| 欧美又色又爽又黄视频| 一进一出抽搐动态| 大型黄色视频在线免费观看| 狂野欧美激情性xxxx| 在线观看一区二区三区| 精品久久久久久久久久免费视频| 91成年电影在线观看| 99国产精品一区二区三区| 亚洲精华国产精华精| 妹子高潮喷水视频| 国内精品久久久久久久电影| 日韩 欧美 亚洲 中文字幕| 天堂影院成人在线观看| 琪琪午夜伦伦电影理论片6080| 国产精品亚洲美女久久久| 欧美黑人巨大hd| 黄网站色视频无遮挡免费观看| 欧美日韩福利视频一区二区| 亚洲成人精品中文字幕电影| 伦理电影免费视频| av有码第一页| 天天一区二区日本电影三级| 大香蕉久久成人网| 人人妻,人人澡人人爽秒播| 人妻久久中文字幕网| 欧美激情 高清一区二区三区| 久久国产亚洲av麻豆专区| 国产黄色小视频在线观看| 一二三四在线观看免费中文在| 好男人在线观看高清免费视频 | 欧美zozozo另类| 成人特级黄色片久久久久久久| 免费在线观看黄色视频的| 欧美又色又爽又黄视频| www.精华液| 99精品欧美一区二区三区四区| 亚洲国产欧美网| 日韩av在线大香蕉| 久久婷婷成人综合色麻豆| 亚洲免费av在线视频| 午夜福利一区二区在线看| 免费观看精品视频网站| 精品国产乱码久久久久久男人| 国产精品98久久久久久宅男小说| 久9热在线精品视频| 91成人精品电影| 日韩欧美在线二视频| 国产av不卡久久| 搞女人的毛片| 成人手机av| 2021天堂中文幕一二区在线观 | 在线免费观看的www视频| 久久久久久久久中文| 欧美性猛交╳xxx乱大交人| netflix在线观看网站| 黑丝袜美女国产一区| 日本免费a在线| 人人澡人人妻人| avwww免费| 桃色一区二区三区在线观看| 好男人电影高清在线观看| 亚洲人成电影免费在线| 免费在线观看影片大全网站| 后天国语完整版免费观看| 国产成人av教育| 老熟妇乱子伦视频在线观看| 91av网站免费观看| 国产亚洲精品综合一区在线观看 | 婷婷六月久久综合丁香| 国产v大片淫在线免费观看| 国产在线观看jvid| 免费看日本二区| 亚洲美女黄片视频| 成年免费大片在线观看| 午夜福利在线在线| 成人欧美大片| 这个男人来自地球电影免费观看| 国产精品国产高清国产av| 免费看美女性在线毛片视频| 免费在线观看视频国产中文字幕亚洲| 黄频高清免费视频| 欧美精品啪啪一区二区三区| 人成视频在线观看免费观看| 色av中文字幕| 日韩国内少妇激情av| 成在线人永久免费视频| 欧美乱码精品一区二区三区| 久久久久久国产a免费观看| 怎么达到女性高潮| 成熟少妇高潮喷水视频| 久久久久久久午夜电影| 亚洲国产中文字幕在线视频| 亚洲中文字幕一区二区三区有码在线看 | 国产精品永久免费网站| 搡老岳熟女国产| 精品日产1卡2卡| 午夜福利一区二区在线看| 观看免费一级毛片| 岛国视频午夜一区免费看| 最近最新中文字幕大全电影3 | 欧美成人一区二区免费高清观看 | 老司机深夜福利视频在线观看| 在线观看免费日韩欧美大片| 亚洲成人精品中文字幕电影| 91老司机精品| 国产亚洲精品久久久久5区| 999久久久精品免费观看国产| 精品久久蜜臀av无| 变态另类丝袜制服| 男女床上黄色一级片免费看| 国产亚洲精品综合一区在线观看 | 天堂影院成人在线观看| 久久中文字幕人妻熟女| 国产亚洲精品久久久久久毛片| 男人操女人黄网站| 亚洲av成人不卡在线观看播放网| 性欧美人与动物交配| 亚洲精品久久成人aⅴ小说| 国产又黄又爽又无遮挡在线| 日韩三级视频一区二区三区| 亚洲av成人一区二区三| 久久精品国产亚洲av高清一级| 国产精品1区2区在线观看.| 国产精品 国内视频| 成人特级黄色片久久久久久久| 欧美日韩福利视频一区二区| 男女午夜视频在线观看| 熟女少妇亚洲综合色aaa.| 亚洲,欧美精品.| 看片在线看免费视频| 黄频高清免费视频| 精品不卡国产一区二区三区| 久9热在线精品视频| 女人高潮潮喷娇喘18禁视频| 色综合欧美亚洲国产小说| 久久婷婷人人爽人人干人人爱| 亚洲在线自拍视频| 久久久精品欧美日韩精品| 美女国产高潮福利片在线看| 这个男人来自地球电影免费观看| 日韩欧美 国产精品| 成人国语在线视频| 他把我摸到了高潮在线观看| 国产激情欧美一区二区| 国产黄色小视频在线观看| 又紧又爽又黄一区二区| 99精品久久久久人妻精品| 精品久久久久久,| 国产高清videossex| 在线天堂中文资源库| 国产精品亚洲美女久久久| 一边摸一边抽搐一进一小说| 男人舔女人的私密视频|