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

    南京夏季氣溶膠吸濕增長因子和云凝結核的觀測研究

    2015-11-18 01:28:54顧雪松孔少飛
    中國環(huán)境科學 2015年2期
    關鍵詞:吸濕性氣溶膠飽和度

    李 琦,銀 燕*,顧雪松,2,袁 亮,孔少飛,江 琪,陳 魁,李 力

    (1.南京信息工程大學,中國氣象局氣溶膠與云降水重點開放實驗室,江蘇 南京 210044;2.無錫市環(huán)境監(jiān)測中心站,江蘇 無錫 214023)

    南京夏季氣溶膠吸濕增長因子和云凝結核的觀測研究

    李 琦1,銀 燕1*,顧雪松1,2,袁 亮1,孔少飛1,江 琪1,陳 魁1,李 力1

    (1.南京信息工程大學,中國氣象局氣溶膠與云降水重點開放實驗室,江蘇 南京 210044;2.無錫市環(huán)境監(jiān)測中心站,江蘇 無錫 214023)

    為了更加全面地研究長三角地區(qū)氣溶膠的理化特性,尤其是吸濕和活化特性,于2013年8月在南京市區(qū)對氣溶膠的吸濕增長因子(GF)和云凝結核(CCN)展開相關觀測研究.使用串聯(lián)電遷移率顆粒物吸濕粒徑分析儀(H-TDMA)觀測32~350nm氣溶膠在相對濕度為90%條件下的吸濕性參數(shù),使用云凝結核計數(shù)器(CCNC)觀測過飽和度在0.2%~0.8%的CCN數(shù)濃度.結果表明,不同氣溶膠粒子的吸濕增長行為均表現(xiàn)出較為明顯的雙峰分布,即一個強吸濕模態(tài)和一個弱吸濕模態(tài),且吸濕性在不同粒徑(愛根核模態(tài)和積聚模態(tài))上存在較為明顯的不同,不溶性物質(zhì)和二次氣溶膠所占比重較大,并且在穩(wěn)定的天氣條件下,氣溶膠的混合狀態(tài)表現(xiàn)為由外混向內(nèi)混發(fā)展的過程.觀測期間該區(qū)域CCN的平均數(shù)濃度為13776(0.6%)cm-3,比沿海區(qū)、山區(qū)、干旱地區(qū)及清潔城市地區(qū)要高很多.其日變化表現(xiàn)為中午時刻出現(xiàn)峰值,影響因素主要與光化學反應有關.同時25日出現(xiàn)的輕霧過程對CCN有較為明顯的清除作用.通過吸濕性參數(shù)計算得到的CCN數(shù)濃度和實際觀測得到的CCN數(shù)濃度進行了閉合實驗,結果顯示出較好的相關性,表明將未飽和條件下觀測得到的吸濕性參數(shù)帶入到K?hler方程中,即可預測過飽和條件下氣溶膠的活化能力.

    南京市;城市氣溶膠;吸濕增長;云凝結核

    目前,中國多數(shù)城市的氣溶膠污染現(xiàn)象日趨嚴重,氣溶膠顆粒物的增加直接導致了空氣質(zhì)量的下降以及霧霾現(xiàn)象的頻發(fā).氣溶膠顆粒物對人體健康、大氣環(huán)境甚至氣候的影響在很大程度上依賴于其譜分布、化學組成、混合方式以及吸濕特性.其中吸濕特性是聯(lián)系氣溶膠微物理、化學參數(shù)的橋梁和紐帶之一,更是氣溶膠光學性質(zhì)的決定性參數(shù)之一.環(huán)境氣溶膠的吸濕增長會影響粒子尺度分布及其光學特性[1],同樣也會影響其活化成為云凝結核的能力[2-3],因此,它也成為決定氣溶膠間接效應的關鍵因素之一.除此之外,與氣溶膠健康效應相關聯(lián)的顆粒物在呼吸道中的位置及其沉積速率也和氣溶膠的吸濕特性直接相關[4-5].

    組成環(huán)境大氣氣溶膠的化學組分可大致分為無機組分和有機組分,其中具有水溶性的組分對于氣溶膠的吸濕特性有著重要的影響.目前對于細粒子的化學分析顯示,不管是來自人為排放或是自然排放的純氣溶膠或混合氣溶膠,都含有大量的硫酸鹽、硝酸鹽與氯鹽,這些無機鹽類都具有吸濕的特性,但隨著觀測技術的發(fā)展,人們逐漸認識到氣溶膠中的有機組分,尤其是可溶性有機組分對氣溶膠吸濕性的重要影響.現(xiàn)有的研究顯示[6-7],無機組分氣溶膠吸濕生長較為明顯,有機組分對氣溶膠吸濕生長因子的影響因不同組分而異,并且當有機組分所占比例上升時氣溶膠的吸濕性反而下降.

    在不飽和條件下吸濕性二次氣溶膠(硫酸鹽、硝酸鹽、銨鹽等)可以作為CCN并進一步成為云滴從而影響云的理化特性,其激活過程可以用寇拉理論來解釋[8].其中,為了能夠?qū)⒘W拥奈鼭窈突罨癄顟B(tài)更好地聯(lián)系起來,Petters等

    [2]和Wex等[9]分別對寇拉公式中某些未知量變換歸納,將所有與化學有關的信息一并歸納入吸濕性參數(shù)k和ρion中,來表征粒子的吸濕活化能力.吸濕性參數(shù)k和ρion在K?hler方程中并非是一個關于相對濕度RH的常數(shù),隨著RH變化而變化的k和ρion對于準確預測CCN數(shù)濃度起著重要的作用.

    近年來,觀測氣溶膠吸濕和活化特性所使用的儀器種類繁多,其中串聯(lián)電遷移率顆粒物吸濕粒徑分析儀(H-TDMA)作為觀測氣溶膠吸濕增長因子和混合狀態(tài)的儀器,被廣泛應用于相關領域的研究[10-23];云凝結核計數(shù)器(CCNC)則是基于水汽和熱擴散速率差別的原理搭建的恒流熱梯度云室,也被大量用于觀測地面和高空的CCN數(shù)濃度[24-31].隨著觀測手段的進步,我國對氣溶膠吸濕增長和活化特性的觀測研究逐漸增多,諸如華北地區(qū)[12]、珠江三角洲等[20],但對長江三角洲區(qū)域內(nèi)的相關觀測和研究仍較為少見.長江三角洲地區(qū)作為我國四大霾嚴重地區(qū)之一[32],霾日數(shù)呈現(xiàn)出逐漸增多的年際變化特征[33].而其中南京的灰霾天數(shù)從2000年的112d增加到2010年的212d,能見度則從1965年的6.5km下降到2005年的4km[34],環(huán)境空氣質(zhì)量呈現(xiàn)出明顯的下降趨勢.

    因此本研究使用H-TDMA和CCNC于2013年夏季對南京市區(qū)亞微米級氣溶膠進行連續(xù)觀測,并分析了在較高相對濕度條件下氣溶膠粒子的吸濕增長因子和混合狀態(tài),不同過飽和度下CCN數(shù)濃度的日變化趨勢,以及通過吸濕性參數(shù)計算CCN數(shù)濃度的閉合實驗,進一步討論氣溶膠粒子在典型城市背景下的構成特點和演變規(guī)律.

    1 資料與方法

    1.1 觀測時間、地點和使用儀器

    本次觀測點設在南京高等職業(yè)技術學校內(nèi),地理坐標為32.0°N,118.7°E,海拔高度為6m,觀測時間為2013年8月1~30日.觀測點所在區(qū)域?qū)儆谀暇┦猩虅罩行膮^(qū)域,其特點為居民建筑居多,道路和車輛密集,并離工業(yè)排放源較遠.

    本次觀測使用的儀器為美國MSP公司生產(chǎn)的串聯(lián)電遷移率顆粒物吸濕粒徑分析儀(H-TDMA)和由DMT公司生產(chǎn)的云凝結核計數(shù)器(CCNC).H-TDMA可用來觀測氣溶膠譜分布、吸濕增長因子和氣溶膠的混合狀態(tài),CCNC則能夠給出不同過飽和度下的CCN數(shù)濃度.觀測期間,2臺儀器均被置于室溫為25℃左右的室內(nèi),其中H-TDMA進行譜掃描時的粒徑范圍是15~700nm,進行吸濕觀測時的粒徑設為32,50,80,110,162,238,350nm,CCNC則設定了4個過飽度(0.2%、0.4%、0.6%、0.8%),且每個過飽和度下的觀測時間為10min.

    H-TDMA的工作原理可以簡述為2個SMPS(Scanning Mobility Particle Sizer Spectrometer)中間加一個吸濕單元.先通過繼電中和器(Kr85,TSI Inc.)使粒子帶電,然后通過差分遷移率分析儀(DMA1-Differential Mobility Analyzer,Model 3081L,TSI Inc)篩選出不同干粒徑(D0)的氣溶膠粒子,再經(jīng)過膜滲透式加濕裝置將樣流加濕至設定濕度,最后通過DMA2來測量吸濕增長后粒子(Dp)的數(shù)濃度譜分布.通過設置的初始粒徑(D0)和吸濕增長后的粒徑(Dp)就可以計算得到吸濕增長因子GF:

    CCNC的工作原理則是基于圓柱形連續(xù)氣流縱向熱梯度云室,利用水汽擴散速率大于熱擴散速率,在云室的垂直中心達到不同的過飽和狀態(tài),而通過改變圓柱筒壁上的溫度梯度,就可以得到不同的過飽和度.

    1.2 數(shù)據(jù)處理方法

    對于H-TDMA的吸濕過程,由于氣溶膠樣流在DMA中會與鞘流相交匯,粒子在此過程中會擴散稀釋,以及管路中的損失等其它影響,觀測直接得到的GF分布,稱為MDF,并不能代表真實的GF概率密度分布(PDF),GF-MDF與GF-PDF相比,在GF中值位置、不同模態(tài)的粒子個數(shù)比(NF)以及模態(tài)譜型的半寬(σ)上均有差異.因此需要對GF-MDF作訂正以得到GF-PDF,常用的方法之一為TDMAFit算法——先利用流量、DMA參數(shù)和粒子干粒徑計算傳輸擴散增寬的影響,并使用對數(shù)正態(tài)分布擬合測得譜分布得到MDF,最后在MDF中修正傳輸擴散的影響.使用TDMAFit算法來處理每個H-TDMA的掃描數(shù)據(jù),并將得到的GF進行歸一化處理,即,其中g代表GF,c(g)則表示GF-PDF.為了方便處理TDMAFit的結果,Gysel等[35]給出了如下常用統(tǒng)計量:

    式中:gmean、σ和Na,b分別表示等效粒徑的GF平均值、GF標準差和a<GF<b下的粒子個數(shù).

    對于CCNC的原始數(shù)據(jù),使用Fortran編程處理,對不同過飽和度下的CCN數(shù)濃度進行小時平均,并考慮到儀器開機后需要一定時間才能達到穩(wěn)定狀態(tài)以及CCNC在不同過飽和度換檔時存在的不穩(wěn)定問題,因此對儀器開機后和每次過飽和度換檔時的前5min數(shù)據(jù)做了相應的剔除,并剔除了“Temps Stabilized”為0時的數(shù)據(jù).

    1.3 吸濕模態(tài)的劃分

    由于環(huán)境氣溶膠的化學成分和混合狀態(tài)比較復雜,因此在實際觀測中往往會出現(xiàn)不同的吸濕模態(tài).通常根據(jù)粒子吸濕性的不同,會將其劃分為幾乎不吸濕類(NH)、弱吸濕類(LH)和強吸濕類(MH).本文使用可溶性物質(zhì)體積比-ε來定義不同吸濕模態(tài)的邊界.假設在相同濕度條件下,將環(huán)境氣溶膠等效為由一定體積的可溶性無機鹽均勻包裹著一定體積的不溶性碳,而在城市背景下的可溶性無機鹽又多為硫酸銨,因此即假定硫酸銨為粒子的無機鹽部分.通過計算相同粒徑的純硫酸銨的吸濕增長因子,從而得到該粒徑下粒子的可溶性物質(zhì)等效體積比),進而可以根據(jù)不同的εAS值(等同GF值)來劃分不同的吸濕模態(tài):

    幾乎不吸濕模態(tài):

    弱吸濕模態(tài):

    式中為32nm粒子在RH=90%下的GF換算值,對于不同粒子以及不同相對濕度這些值均會有所變化.本文即依據(jù)上述分類來劃分其相關的其他吸濕性參數(shù)(σMH、σLH、NFMH、NFLH).

    2 結果與分析

    2.1 觀測期間氣象要素的變化特征

    圖1 觀測期間氣象要素和NCCN時間序列Fig.1 Temporal variation of meteorological parameters and number concentration of CCN during the observation period

    受到全球氣候變暖、城市熱島效應和西太平洋副熱帶高壓的持續(xù)影響,2013年南京夏季出現(xiàn)了連續(xù)的高溫現(xiàn)象,不僅高溫日數(shù)(≥35℃的天數(shù))持續(xù)時間長,而且一些天中的極端高溫值甚至超過了40℃.此外,由于中高緯度的西風氣流較為平直,東西向環(huán)流增強,南北向環(huán)流減弱,因此北方的冷空氣難以南下與南京及其下游的暖空氣交會,從而導致降水的減少.圖1為觀測期間氣象要素和不同過飽和度下NCCN隨時間的變化趨勢,可以看出,觀測期間多為偏南風,風速較低,最大值不超過4m/s,不利于污染物的擴散和傳輸.溫度和相對濕度均呈現(xiàn)出明顯的日變化,高溫日數(shù)有16d,并在8月10日出現(xiàn)最高溫度41℃.同時在較高的相對濕度下,粒子吸濕長大,增加了光的散射強度,從而導致能見度的降低.觀測期間的CCN數(shù)濃度在高過飽和度下(SS=0.6%和SS=0.8%)出現(xiàn)了2個明顯的峰值,分別為8月6日和9日的中午時刻,且該時刻的氣象條件為高溫、低濕和高能見度的清潔條件,有利于新粒子的生成[36].8月19~23日CCN數(shù)濃度較低,是因為在該時段內(nèi)有陣性降水過程,對氣溶膠粒子有清除作用.8月25日和26日有較為明顯的輕霧過程,表現(xiàn)為相對濕度高,能見度的日變化較大,CCN數(shù)濃度在上、下午也出現(xiàn)了明顯的差異.

    2.2 亞微米級氣溶膠吸濕增長因子的觀測

    Swietlicki等[37]認為,在城市背景下,亞微米級氣溶膠粒子會表現(xiàn)出不同的吸濕特性,其吸濕增長因子的概率分布常常出現(xiàn)2個或3個不同的模態(tài).根據(jù)吸濕性的不同,一般將這些模態(tài)定義為幾乎不吸濕模態(tài)、弱吸濕模態(tài)和強吸濕模態(tài)3類.其中,幾乎不吸濕模態(tài)通常含有大量燃燒過程排放的細粒子,來源包括機動車尾氣和人為活動排放,這部分粒子一般含有較高比例的不溶物,如黑碳及部分不可溶有機物.弱吸濕模態(tài)通常含有老化的黑碳、部分氧化的有機物和一些凝結在碳表面的無機鹽等,其中老化過程主要是水汽的凝結和顆粒物相互間的碰并作用.強吸濕模態(tài)則往往含有較高比例的可溶性無機鹽,如硫酸鹽、硝酸鹽等.

    由于幾乎不吸濕和弱吸濕模態(tài)的分界線并不明顯,且隨著粒徑增大,兩者在整體吸濕性中所占比例又逐漸降低,因此本文將這兩個模態(tài)一并歸為弱吸濕模態(tài)進行相關統(tǒng)計.表1即為對整個觀測期間分檔氣溶膠吸濕性參數(shù)的定量統(tǒng)計,對應不同粒徑(32,50,80,110,162,238,350nm),弱吸濕增長因子gLH的值分別為1.030、1.039、1.044、1.058、1.059、1.080、1.050,強吸濕增長因子gMH的值分別為1.270、1.315、1.343、1.355、1.372、1.405、1.438,發(fā)現(xiàn)gLH和gMH的值在小粒徑上同2004年北京夏季[15]的觀測值相當,但在大粒徑上卻較小,平均值均小于2012年南京夏季(南京信息工程大學)的觀測值[23],分析原因可能是因為南京市區(qū)遠離工業(yè)排放區(qū),并且觀測期間風速較小,限制了氣溶膠的遠距離輸送,從而使該區(qū)域的氣溶膠構成以局地排放為主,多為交通和人為活動排放,直接導致不可溶成分比例的增加,進而導致整體吸濕性的下降,特別是積聚模態(tài)的粒子,數(shù)值上表現(xiàn)為平均吸濕增長因子的平均值(GF_mean)較低.同時,由于夏季太陽輻射增強,隨之發(fā)生的光化學反應產(chǎn)生了大量的新粒子,這些由氣粒轉(zhuǎn)化產(chǎn)生的超細粒子對小粒徑段(核模態(tài)和愛根核模態(tài))氣溶膠的吸濕性有一定的貢獻.同時GF_mean的變化在愛根核和積聚模態(tài)間存在著一個較為明顯的躍遷,這表明了不同模態(tài)的化學組分和混合狀態(tài)往往是不同的.

    表1 觀測期間分檔氣溶膠粒子吸濕性統(tǒng)計結果(RH=90%±1%)Table 1 The statistical results of size-resolved aerosol hygroscopic properties during the observation.

    圖2為觀測期間50,110,238nm的平均GF-PDF,其分布特征表現(xiàn)為明顯的雙峰分布,這與Chen等[38]和Santarpia等[39]的研究結果相似.由圖2可知,50nm粒子弱吸濕模態(tài)的面積要大于238nm的,同時大粒子峰的寬度要大于小粒子的,說明大粒子的外混程度更高.圖中GF=1.20是強弱吸濕態(tài)的分界線,由于弱吸濕態(tài)粒子正處于粒子老化的短暫階段,并且其粒子組成具有明顯的日變化特征,因此隨著相對濕度的進一步增大,部分與強吸濕態(tài)重疊區(qū)域的粒子會繼續(xù)吸濕長大,從而進入強吸濕態(tài),可稱其為過渡態(tài)粒子.而在重疊區(qū)域的右側,強吸濕峰左緣會受弱吸濕峰延伸段的影響被抬高.

    圖2 觀測期間平均GF-PDF, RH=90%Fig.2 Growth factor probability density function (PDF)for 50nm、110nm、238nm particles exposed at 90% relative humidity

    圖3 分檔氣溶膠在強、弱吸濕態(tài)的吸濕增長因子分布Fig.3 The growth factor of the more-hygroscopic and the low-hygroscopic group at each dry diameter

    圖3為觀測期間分檔氣溶膠在強、弱吸濕態(tài)的吸濕增長因子分布,圖中實線為εAS=0.24,認為干粒徑越大的粒子中含有越多的具有較強吸濕性的可溶性無機鹽,相比于華北平原的觀測結果,εAS在愛根核模態(tài)和積聚模態(tài)的值均較小,這可能與大氣中污染物類型不同有關,南京市區(qū)的污染物多為汽車尾氣和人為活動排放,而在華北平原,因為有較多的熱電站,因此在大氣中會有較多的含硫化合物[40],如硫酸銨,并通過凝結效應使愛根核粒子的吸濕性增強.同時由于觀測期間風速較小,所以該區(qū)域的粒子構成受到遠距離輸送的影響較小,同樣,較小的風速也限制了污染物的擴散,使其更易于積累,因而吸濕性的表現(xiàn)也多由粒子間相互碰并導致化學成分和混合狀態(tài)發(fā)生變化所致.其混合方式相比于2012年夏季在南京信息工程大學的觀測結果更趨向于內(nèi)混狀態(tài).

    由圖4可見,隨著粒徑的增長,強吸濕態(tài)的粒子個數(shù)有逐漸增大的趨勢,這與Massling等[15]的研究結果相一致.由于大量汽車尾氣和人為活動排放多為小粒子,并且成分多為黑碳及不可溶有機物等弱吸濕類,因此在小于100nm的愛根核模態(tài)中,弱吸濕顆粒物占了更大的比例.同時,發(fā)現(xiàn)32nm粒子在強吸濕態(tài)個數(shù)的標準差最大,表明其日變化最大,可能原因是由于光化學反應活躍導致中午時刻的細粒子數(shù)目增加,這是除汽車尾氣和人為活動外產(chǎn)生細粒子的又一個主要源.同時,認為這些由氣粒轉(zhuǎn)化生成的二次氣溶膠(SOA)對弱吸濕部分也具有一定的貢獻,這也是弱吸濕態(tài)粒子個數(shù)有較高比例并具有較大日變化的部分原因,因為環(huán)境大氣中一次氣溶膠和二次氣溶膠都具有很明顯的日變化特征.

    圖4 分檔氣溶膠在強、弱吸濕態(tài)的數(shù)目比例Fig.4 The number fraction of the more-hygroscopic and the low-hygroscopic group at each dry diameter

    通過研究80nm粒子在中午時刻和傍晚下班時刻的吸濕增長因子,發(fā)現(xiàn)中午時刻的GF_mean要高于傍晚,這可能是由于中午受到溫度以及光照的影響,環(huán)境氣溶膠的組分發(fā)生改變,導致吸濕性組分增多,從而影響其吸濕增長因子,而傍晚時刻由于處于下班高峰期,受到機動車排放的影響,并根據(jù)先前對機動車尾氣吸濕性的實驗室研究,發(fā)現(xiàn)柴油車及汽油車吸濕性為近憎水性,生長因子在1.0左右[7],因此影響了氣溶膠的吸濕性.此外,王軒等[7]研究中也發(fā)現(xiàn)處于相對較強吸濕性時,硝酸根、硫酸根和鉀離子的相對比例要高于低吸濕性的,而鋁離子、銨鹽以及硅酸根的相對比例則較低.并且大氣中鋁離子以及硅酸鹽多來自揚塵[41],大氣中的硫酸鹽和硝酸鹽可能來源于及機動車的排放.鋁離子可以看作是氣溶膠一次排放組分,而硫酸鹽和硝酸鹽則多為二次粒子.因此影響環(huán)境氣溶膠吸濕性的主要成分中,硫酸鹽和硝酸鹽類是增強吸濕性的,而來自揚塵的鋁離子和硅酸鹽則是降低吸濕性的.

    2.3 CCN譜日變化特征

    通過對CCN數(shù)濃度在不同過飽和度下的觀測及分析,將進一步解釋夏季南京市區(qū)氣溶膠粒子的構成特點和變化特征.觀測期間較高的氣溶膠數(shù)濃度(NPM2.5=2466.05±967.65cm-3)導致了較高的CCN數(shù)濃度,圖5為觀測期間不同過飽和度下CCN數(shù)濃度的日變化趨勢,NCCN隨過飽和度的增加而增加,其中最高的CCN數(shù)濃度(SS= 0.8%)達到20849.84cm-3,最低為(SS=0.2%)2181.05cm-3.整個觀測期間,4個過飽和度下(0.2%、0.4%、0.6%和0.8%)的CCN數(shù)濃度平均值和標準差分別為(2524±181), (8367±628),(13776±1047), (18278± 1469)cm-3.表2總結了近年來我國使用DMT- CCNC在各地區(qū)的觀測結果,天氣條件均為晴天.其中本次南京觀測的CCN數(shù)濃度同石家莊的觀測值相當,比武清觀測的要低,但比其他沿海區(qū)、山區(qū)、干旱地區(qū)及清潔城市地區(qū)的數(shù)濃度要高很多.外場觀測表明,大部分的CCN 是自然或人為產(chǎn)生的硫酸鹽氣溶膠.在大陸,尤其是工業(yè)區(qū),工業(yè)排放的SO2造成大氣中硫酸鹽含量的增加,同時,局地污染源和遠距離污染源對S的貢獻也比較大,但沙塵氣溶膠由于其本身的不可溶成分不利于S在氣溶膠中的富集[28].因此南京市區(qū)的CCN數(shù)濃度要低于污染物較多的工業(yè)城市,但由于人類活動頻繁會使CCN數(shù)濃度高于其他清潔地區(qū).

    圖5 觀測期間CCN數(shù)濃度日變化趨勢Fig.5 The diurnal variation of the concentration of cloud condensation nuclei (CCN)

    表2 不同地區(qū)CCN觀測情況統(tǒng)計Table 2 Summary of CCN concentrations measured in different regions in China

    分析指出,導致CCN數(shù)濃度發(fā)生變化的原因主要有局地排放、夏季光化學反應增強導致的新粒子生成事件、邊界層高度日變化導致的氣溶膠粒子的擴散和聚集.其中CCN數(shù)濃度在高過飽和度下出現(xiàn)較為明顯的峰值是在中午時刻,表明在該時段內(nèi)由于太陽輻射增強,并在環(huán)境大氣中有大量氣相前體物存在的前提下,往往會發(fā)生新粒子生成事件,生成的新粒子多為粒徑較小的細粒子或超細粒子,這部分粒子和吸濕增長成為愛根核模態(tài)的部分粒子增加了過飽和度0.8%下的CCN數(shù)濃度.有相關研究同樣指出,凝結核25%來源于次生源[42],即在大氣光化學作用和其他化學過程的作用下,促使大氣中的氣體成分轉(zhuǎn)化成小的質(zhì)粒.同時,交通排放呈現(xiàn)出早晚兩個高峰,早晨時段的峰值主要是由于交通排放和陸面揚塵在風速小和高度低的邊界層內(nèi)不易擴散導致,夜間時段由于太陽輻射消失,溫度降低,逆溫層開始形成,再加之下班高峰期產(chǎn)生的交通排放,使污染物容易在邊界層內(nèi)停留聚集,因此峰值較大.而相比于第2個峰值,第1個峰值雖然也主要由上班高峰期中交通排放造成,但此時太陽輻射逐漸增強,逆溫層開始消失,邊界層中熱力對流加強,混合層快速增長,污染物不易在近地層停留積聚,因此峰值相對較低.

    不同過飽和度間的CCN數(shù)濃度也存在一個較大的差值,表明該地區(qū)粒子的構成主要為集中在愛根核模態(tài)的細粒子,并且日變化特征顯著,這在對粒子吸濕性的討論中也有所表述.氣象要素對CCN數(shù)濃度的影響主要表現(xiàn)為較大風速對污染物的擴散和降水的濕清除作用,從觀測期間對風速和降水量的統(tǒng)計結果可以看出,除8月19~23日間偶有短時降水外,其余時間均無明顯降水過程,并且觀測期間近地層的平均風速僅為(1.25 ±0.64)m/s,不利于污染物的擴散,這也是導致CCN數(shù)濃度偏高的原因之一.同時,觀測到25日清晨出現(xiàn)了一次輕霧過程,對CCN數(shù)濃度的影響表現(xiàn)為上、下午的NCCN出現(xiàn)了較為明顯的差異,這也驗證了在高濕條件下一部分CCN能活化成為云霧滴,即云霧過程對CCN存在清除作用.

    2.4 CCN與吸濕增長的閉合

    除了通過CCN擬合公式來預測不同過飽和度下的CCN數(shù)濃度,還可以通過對吸濕增長的觀測來計算CCN數(shù)濃度.其理論是通過引入一個單參數(shù)kappa(k)來表征粒子的吸濕和活化能力,將實際觀測得到的吸濕性參數(shù)帶入K?hler方程中,即可以計算出不同模態(tài)的k值,從而得到不同粒徑的k-PDF分布,假定觀測設定的各粒徑檔即為臨界干粒徑,即可通過公式(8)得到不同過飽度下各粒徑粒子活化所需達到的κsc值:

    結合k-PDF就可以計算出不同粒徑粒子在不同過飽和度下的活化率,用Sigmoidal公式對活化率曲線進行擬合,即可得到不同過飽和度下的中值粒徑,然后分別選取相同時段下氣壓空腔氣溶膠探頭地基型(PCASP-X2)和SMPS大于活化粒徑的數(shù)濃度和CCNC直接觀測得到的數(shù)濃度進行線性擬合,得到如圖6a、6b的關系.進一步分析,發(fā)現(xiàn)通過吸濕性參數(shù)計算的CCN數(shù)濃度和實際觀測得到的CCN數(shù)濃度存在著較好的相關性,即可通過未飽和條件下觀測得到的吸濕性參數(shù)來預測過飽和條件下的CCN數(shù)濃度,但本文中計算的相關系數(shù)并不是很理想,其主要原因可能是:(1)通過擬合得到的活化率曲線在判斷活化粒徑上并不是很準確;(2)測量氣溶膠譜的儀器——PCASP-X2和SMPS在不同粒徑檔上的劃分并不細致,因此不能準確找到相對應的活化粒徑,進而不能準確計算出大于該活化粒徑下的氣溶膠數(shù)濃度(NCN).

    圖6 由吸濕性參數(shù)結合氣溶膠譜計算的NCCN和實際觀測的NCCN之間的閉合Fig.6 NCCNcalculated from hygroscopic parameters combined with aerosol spectrum compared with measured NCCN

    3 結論

    3.1 夏季南京市區(qū)氣溶膠粒子的組成主要以細粒子為主,來源多為交通和人為活動排放以及新粒子的生成,其數(shù)濃度的日變化主要依賴于局地排放、光化學反應以及邊界層高度的變化.

    3.2 通過對亞微米級氣溶膠吸濕增長因子的觀測發(fā)現(xiàn),南京市區(qū)不溶性物質(zhì)所占比重較大,尤其是在小于100nm的愛根核模態(tài),這部分不溶物主要來源于局地排放.并且認為夏季太陽輻射的增強使得光化學反應活躍,通過氣粒轉(zhuǎn)化過程,生成大量的二次氣溶膠,這部分氣溶膠對粒子整體的吸濕性也有一定的貢獻.

    3.3 對CCN數(shù)濃度的日變化分析可以得到,在高過飽和度下表現(xiàn)出明顯的日變化趨勢,表明粒子的組成多以細粒子為主且易受到不同程度局地排放變化的影響.同時,不同過飽和度下CCN數(shù)濃度的較大差值及光化學反應會直接產(chǎn)生部分小的質(zhì)粒這兩點也解釋了夏季氣溶膠粒子在南京市區(qū),即典型城市背景下的構成特點.而氣象要素對CCN數(shù)濃度的影響,主要是降水的濕清除作用和風的擴散作用,由于觀測期間無明顯降水過程,且平均風速較低,致使污染物不易擴散,從而使CCN數(shù)濃度始終保持在較高的水平.

    3.4 通過使用H-TDMA觀測得到的氣溶膠吸濕性參數(shù)來計算CCN數(shù)濃度,并與CCNC觀測得到的CCN數(shù)濃度相比較,這種閉合方案通過寇拉理論試圖將氣溶膠吸濕和活化特性聯(lián)系到一起.雖然在計算過程中存在著一定的誤差,但也表現(xiàn)出了較好的相關性,即用未飽和條件下氣溶膠的吸濕性參數(shù)可以較好的預測過飽和條件下能夠活化成為云凝結核的氣溶膠數(shù)濃度.

    [1]Gassó S, Hegg D A, Covert D S, et al. Influence of humidity on the aerosol scattering coefficient and its effect on the upwelling radiance during ACE-2 [J]. Tellus B, 2000,52(2):546-567.

    [2]Petters M D, Kreidenweis S M. A single parameter representation of hygroscopic growth and cloud condensation nucleus activity[J]. Atmospheric Chemistry and Physics, 2007,7(8):1961-1971.

    [3]Svenningsson B, Rissler J, Swietlicki E, et al. Hygroscopic growth and critical supersaturations for mixed aerosol particles of inorganic and organic compounds of atmospheric relevance [J]. Atmospheric Chemistry and Physics, 2006,6(7):1937-1952.

    [4]Chan H K, Eberl S, Daviskas E, et al. Changes in lung deposition of aerosols due to hygroscopic growth: a fast SPECT study [J]. Journal of aerosol medicine, 2002,15(3):307-311.

    [5]L?ndahl J, Massling A, Pagels J, et al. Size-resolved respiratory-tract deposition of fine and ultrafine hydrophobic and hygroscopic aerosol particles during rest and exercise [J]. Inhalation toxicology, 2007,19(2):109-116.

    [6]Carrico C M, Kreidenweis S M, Malm W C, et al. Hygroscopic growth behavior of a carbon-dominated aerosol in Yosemite National Park [J]. Atmospheric Environment, 2005,39(8):1393-1404.

    [7]王 軒.氣溶膠吸濕特性研究 [D]. 北京:中國環(huán)境科學研究院,2010.

    [8]Seinfeld J H, Pandis S N. Atmospheric chemistry and physics:from air pollution to climate change [M]. John Wiley & Sons,2012.

    [9]Wex H, Hennig T, Salma I, et al. Hygroscopic growth and measured and modeled critical super-saturations of an atmospheric HULIS sample [J]. Geophysical research letters,2007,34:2818.

    [10]Gysel M, Crosier J, Topping D O, et al. Closure study between chemical composition and hygroscopic growth of aerosol particles during TORCH2 [J]. Atmospheric Chemistry and Physics, 2007,7(24):6131-6144.

    [11]Henning S, Ziese M, Kiselev A, et al. Hygroscopic growth and droplet activation of soot particles: uncoated, succinic or sulfuric acid coated [J]. Atmospheric Chemistry and Physics, 2012,12(10):4525-4537.

    [12]Liu P F, Zhao C S, G?bel T, et al. Hygroscopic properties of aerosol particles at high relative humidity and their diurnal variations in the North China Plain [J]. Atmospheric Chemistry and Physics, 2011,11(7):3479-3494.

    [13]Lance S, Raatikainen T, Onasch T B, et al. Aerosol mixing state,hygroscopic growth and cloud activation efficiency during MIRAGE 2006 [J]. Atmospheric Chemistry and Physics, 2013,13(9):5049-5062.

    [14]Massling A, Leinert S, Wiedensohler A, et al. Hygroscopic growth of sub-micrometer and one-micrometer aerosol particles measured during ACE-Asia [J]. Atmospheric Chemistry and Physics, 2007,7(12):3249-3259.

    [15]Massling A, Stock M, Wehner B, et al. Size segregated water uptake of the urban submicrometer aerosol in Beijing [J]. Atmospheric Environment, 2009,43(8):1578-1589.

    [16]Meier J, Wehner B, Massling A, et al. Hygroscopic growth ofurban aerosol particles in Beijing (China) during wintertime: a comparison of three experimental methods [J]. Atmospheric Chemistry and Physics, 2009,9(18):6865-6880.

    [17]Mircea M, Facchini M C, Decesari S, et al. Importance of the organic aerosol fraction for modeling aerosol hygroscopic growth and activation: a case study in the Amazon Basin [J]. Atmospheric Chemistry and Physics, 2005,5(11):3111-3126.

    [18]Su H, Rose D, Cheng Y F, et al. Hygroscopicity distribution concept for measurement data analysis and modeling of aerosol particle mixing state with regard to hygroscopic growth and CCN activation [J]. Atmospheric Chemistry and Physics, 2010,10(15):7489-7503.

    [19]Svenningsson B, Rissler J, Swietlicki E, et al. Hygroscopic growth and critical supersaturations for mixed aerosol particles of inorganic and organic compounds of atmospheric relevance [J]. Atmospheric Chemistry and Physics, 2006,6(7):1937-1952.

    [20]Tan H, Yin Y, Gu X, et al. An Observational Study of the Hygroscopic Properties of Aerosols over the Pearl River Delta Region [J]. Atmospheric Environment, 2013,77:817-826.

    [21]Ye X, Ma Z, Hu D, et al. Size-resolved hygroscopicity of submicrometer urban aerosols in Shanghai during wintertime [J]. Atmospheric Research, 2011,99(2):353-364.

    [22]Zamora I R, Jacobson M Z. Measuring and modeling the hygroscopic growth of two humic substances in mixed aerosol particles of atmospheric relevance [J]. Atmospheric Chemistry and Physics, 2013,13(17):8973-8989.

    [23]吳奕霄,銀 燕,顧雪松,等.南京北郊大氣氣溶膠的吸濕性觀測研究 [J]. 中國環(huán)境科學, 2014,34(8):1938-1949.

    [24]Baumgardner D, Raga G B, Muhlia A. Evidence for the formation of CCN by photochemical processes in Mexico City [J]. Atmospheric Environment, 2004,38(3):357-367.

    [25]Kuwata M, Kondo Y, Miyazaki Y, et al. Cloud condensation nuclei activity at Jeju Island, Korea in spring 2005 [J]. Atmospheric Chemistry and Physics, 2008,8(11):2933-2948.

    [26]Reade L, Jennings S G, McSweeney G. Cloud condensation nuclei measurements at Mace Head, Ireland, over the period 1994-2002[J]. Atmospheric research, 2006,82(3):610-621.

    [27]石立新,段 英.華北地區(qū)云凝結核的觀測研究 [J]. 氣象學報,2007,65(4):644-652.

    [28]岳巖裕,牛生杰,桑建人,等.干旱區(qū)云凝結核分布及其影響因子的觀測研究 [J]. 中國環(huán)境科學, 2010,30(5):593-598.

    [29]趙永欣,牛生杰,呂晶晶,等.2007年夏季我國西北地區(qū)云凝結核的觀測研究 [J]. 高原氣象, 2010,(4):1043-1049.

    [30]王婷婷.華北地區(qū)云凝結核特性研究 [D]. 北京:中國氣象科學研究院, 2011.

    [31]李 力,銀 燕,顧雪松,等.黃山地區(qū)不同高度云凝結核的觀測分析 [J]. 大氣科學, 2014,38(3):410-420.

    [32]陳 魁,銀 燕,魏玉香,等.南京大氣PM2.5中碳組分觀測分析[J]. 中國環(huán)境科學, 2010,30(8):1015-1020.

    [33]史 軍,崔林麗,賀千山,等.華東霧和霾日數(shù)的變化特征及成因分析 [J]. 地理學報, 2010,65(5):533-542.

    [34]王 玨,孟維宸,張奇漪,等.南京市灰霾期間顆粒物污染的主要影響源識別 [J]. 環(huán)境保護科學, 2012,38(4):6-11.

    [35]Gysel M, McFiggans G B, Coe H. Inversion of tandem differential mobility analyser (TDMA) measurements [J]. Journal of Aerosol Science, 2009,40(2):134-151.

    [36]Kulmala M, Vehkam?ki H, Pet?j? T, et al. Formation and growth rates of ultrafine atmospheric particles: a review of observations[J]. Journal of Aerosol Science, 2004,35(2):143-176.

    [37]Swietlicki E, HANSSON H C, H?meri K, et al. Hygroscopic properties of submicrometer atmospheric aerosol particles measured with H-TDMA instruments in various environments— A review [J]. Tellus B, 2008,60(3):432-469.

    [38]Chen L Y, Jeng F T, Chen C C, et al. Hygroscopic behavior of atmospheric aerosol in Taipei [J]. Atmospheric Environment, 2003,37(15):2069-2075.

    [39]Santarpia J L, Li R, Collins D R. Direct measurement of the hydration state of ambient aerosol populations [J]. Journal of Geophysical Research:Atmospheres(1984~2012), 2004,109(D18):18209.

    [40]Zhang X Y, Wang Y Q, Niu T, et al. Atmospheric aerosol compositions in China: spatial/temporal variability, chemical signature, regional haze distribution and comparisons with global aerosols [J]. Atmospheric Chemistry and Physics, 2012,12(2):779-799.

    [41]華 蕾,郭 婧,徐子優(yōu),等.北京市主要PM10排放源成分譜分析[J]. 中國環(huán)境監(jiān)測, 2006,(6):64-71.

    [42]Brock J R. Condensational growth of atmospheric atmospheric aerosols [J]. Journal of Colloid and Interface Science, 1972,39(1):32-36.

    An observational study of aerosol hygroscopic growth factor and cloud condensation nuclei in Nanjing in summer.

    LI Qi1, YIN Yan1*, GU Xue-song1,2, YUAN Liang1, KONG Shao-fei1, JIANG Qi1, CHEN Kui1, LI Li1
    (1.Key Laboratory for Aerosol-Cloud-Precipitation, China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing 210044, China;2.Wuxi Environmental Monitoring Center Station, Wuxi 214023, China).

    China Environmental Science, 2015,35(2):337~346

    In order to study the physical and chemical properties of aerosols in Yangtze River delta (YRD) more comprehensively, especially the hygroscopic growth and activation feature. In this study, aerosol hygroscopic growth factor (GF) and cloud condensation nuclei (CCN) in urban Nanjing were studied in August 2013. A Hygroscopic Tandem Differential Mobility Analyzer (H-TDMA) was utilized to observe the hygroscopic parameters at 90% relative humidity(RH) for particles with dry diameter between 32~350nm. And a Cloud Condensation Nuclei Counter (CCNC) was employed to measure the number concentration of CCN (NCCN) at supersaturation of 0.2%~0.8% simultaneously. The statistical results show that, the hygroscopic growth behavior of aerosol particles is shown as a distinct bimodal distribution, and can be classified into a dominant more-hygroscopic (MH) group and a smaller less-hygroscopic (LH)group. In addition, the hygroscopic growth factor of particles of different diameters, especially the Aitken mode and accumulation mode particles, exhibited clear differences. Meanwhile, insoluble substances and secondary aerosols accounted for a large proportion. It was also found that, in a stable weather condition, aerosol tended to be in an internally mixed state. During the period of observation, the average NCCNwas 13776cm-3(at 0.6% supersaturation), which was higher than other coastal, mountains, arid and clean urban areas. The peak diurnal variation of NCCNappeared at noon,which was mainly related to the photochemical reactions. The observations also show that the light fog process on August 25 had a clear wet scavenging effect on CCN. A comparison of the NCCNcalculated from hygroscopic parameter-k with actual measured NCCN, indicating that combination of the hygroscopic parameter-k under sub-saturation from actualmeasurement and the K?hler equation may well predict the activation ability of aerosols to CCN under supersaturation.

    Nanjing;urban aerosol;hygroscopic growth;cloud condensation nuclei

    X513

    A

    1000-6923(2015)02-0337-10

    李 琦(1990-),男,內(nèi)蒙古呼和浩特市人,南京信息工程

    2014-06-10

    國家自然科學基金重點項目(41030962);江蘇高校優(yōu)勢學科建設工程資助項目(PAPD)

    * 責任作者, 教授, yinyan@nuist.edu.cn

    大學碩士研究生,主要從事氣溶膠吸濕和活化特征研究.

    猜你喜歡
    吸濕性氣溶膠飽和度
    糖臬之吻
    氣溶膠傳播之謎
    氣溶膠中210Po測定的不確定度評定
    添加不同輔料牛膝浸膏粉體表面特征與吸濕性的相關性
    中成藥(2018年2期)2018-05-09 07:19:44
    中藥浸膏粉吸濕性的評價方法
    中成藥(2017年5期)2017-06-13 13:01:12
    黃芪多糖防潮樹脂復合物的制備及其吸濕性
    中成藥(2017年3期)2017-05-17 06:08:53
    四川盆地秋季氣溶膠與云的相關分析
    制作一個泥土飽和度測試儀
    巧用有機物的不飽和度
    柔情粉色
    優(yōu)雅(2014年4期)2014-04-18 10:18:25
    99久久精品热视频| 免费av不卡在线播放| 哪里可以看免费的av片| 人妻夜夜爽99麻豆av| 在线国产一区二区在线| 一卡2卡三卡四卡精品乱码亚洲| 久久午夜亚洲精品久久| 两个人视频免费观看高清| 18禁国产床啪视频网站| 亚洲自拍偷在线| 亚洲色图av天堂| 在线国产一区二区在线| 日韩大尺度精品在线看网址| 成人亚洲精品av一区二区| www.www免费av| 黄色视频,在线免费观看| 亚洲,欧美精品.| 免费av毛片视频| 床上黄色一级片| 国产亚洲欧美在线一区二区| 亚洲片人在线观看| 日本免费一区二区三区高清不卡| 又黄又爽又免费观看的视频| 两个人看的免费小视频| 国产69精品久久久久777片| а√天堂www在线а√下载| 99精品欧美一区二区三区四区| 人人妻人人澡欧美一区二区| 日本成人三级电影网站| 国产亚洲精品av在线| 国产综合懂色| 狂野欧美激情性xxxx| 99riav亚洲国产免费| 乱人视频在线观看| 五月玫瑰六月丁香| 欧美日本亚洲视频在线播放| 亚洲av五月六月丁香网| 麻豆一二三区av精品| 日韩欧美在线二视频| 三级毛片av免费| 亚洲精品久久国产高清桃花| 90打野战视频偷拍视频| 男女床上黄色一级片免费看| 啦啦啦免费观看视频1| 麻豆久久精品国产亚洲av| 成人鲁丝片一二三区免费| 女生性感内裤真人,穿戴方法视频| 亚洲精品亚洲一区二区| 国产精品亚洲美女久久久| 成人亚洲精品av一区二区| 中国美女看黄片| 国产精品一区二区三区四区久久| 国产精品亚洲一级av第二区| 欧美日韩瑟瑟在线播放| 高潮久久久久久久久久久不卡| 国产成人啪精品午夜网站| 亚洲第一欧美日韩一区二区三区| 日韩欧美在线二视频| 国产日本99.免费观看| 亚洲欧美日韩高清专用| 亚洲五月天丁香| 超碰av人人做人人爽久久 | 亚洲精品成人久久久久久| 国产一级毛片七仙女欲春2| 婷婷精品国产亚洲av| 亚洲国产欧美网| 日日夜夜操网爽| 亚洲在线观看片| 蜜桃久久精品国产亚洲av| 欧美3d第一页| 国产成人aa在线观看| a在线观看视频网站| 亚洲av第一区精品v没综合| 久久精品综合一区二区三区| 99在线人妻在线中文字幕| 91麻豆精品激情在线观看国产| 精品久久久久久,| 日本免费一区二区三区高清不卡| 成人高潮视频无遮挡免费网站| 又黄又粗又硬又大视频| 熟妇人妻久久中文字幕3abv| 欧美最新免费一区二区三区 | 欧美大码av| 丁香六月欧美| 午夜精品一区二区三区免费看| 中出人妻视频一区二区| 51国产日韩欧美| 国产成人影院久久av| 女人十人毛片免费观看3o分钟| 久久精品综合一区二区三区| 亚洲欧美日韩卡通动漫| 天堂√8在线中文| 搡女人真爽免费视频火全软件 | 高潮久久久久久久久久久不卡| 亚洲国产中文字幕在线视频| 成人鲁丝片一二三区免费| 在线免费观看不下载黄p国产 | 久久国产精品影院| 国产探花极品一区二区| netflix在线观看网站| 啦啦啦韩国在线观看视频| 黑人欧美特级aaaaaa片| 日韩欧美国产在线观看| 国内久久婷婷六月综合欲色啪| 99久国产av精品| 久久草成人影院| 首页视频小说图片口味搜索| 成人午夜高清在线视频| 精品国产三级普通话版| 亚洲成人中文字幕在线播放| 久久国产乱子伦精品免费另类| 日韩免费av在线播放| 精品一区二区三区人妻视频| 日韩中文字幕欧美一区二区| 波野结衣二区三区在线 | 99国产极品粉嫩在线观看| ponron亚洲| 国产乱人伦免费视频| 香蕉久久夜色| 天堂av国产一区二区熟女人妻| 99热精品在线国产| 中文字幕久久专区| 国产成人av激情在线播放| 黑人欧美特级aaaaaa片| 久久久久久久亚洲中文字幕 | 精品日产1卡2卡| 夜夜夜夜夜久久久久| 在线免费观看不下载黄p国产 | 国产视频一区二区在线看| www日本在线高清视频| 在线十欧美十亚洲十日本专区| 亚洲va日本ⅴa欧美va伊人久久| 特级一级黄色大片| 99热精品在线国产| 脱女人内裤的视频| 级片在线观看| 高清在线国产一区| e午夜精品久久久久久久| 亚洲五月婷婷丁香| 免费搜索国产男女视频| 淫秽高清视频在线观看| 免费av毛片视频| 啦啦啦免费观看视频1| 久久久久亚洲av毛片大全| 国产激情偷乱视频一区二区| 在线播放国产精品三级| 老司机午夜福利在线观看视频| x7x7x7水蜜桃| 日韩亚洲欧美综合| 嫁个100分男人电影在线观看| 母亲3免费完整高清在线观看| 午夜亚洲福利在线播放| 国产黄a三级三级三级人| 色噜噜av男人的天堂激情| 人人妻人人澡欧美一区二区| 国产老妇女一区| 亚洲国产高清在线一区二区三| 悠悠久久av| 国产不卡一卡二| 国产精品,欧美在线| 国产黄a三级三级三级人| 亚洲av一区综合| а√天堂www在线а√下载| 亚洲精华国产精华精| 亚洲黑人精品在线| 欧美乱妇无乱码| 中文字幕精品亚洲无线码一区| 一级黄片播放器| 欧美极品一区二区三区四区| 久久久久精品国产欧美久久久| 日日夜夜操网爽| 天天一区二区日本电影三级| 一区二区三区激情视频| 国产高清视频在线播放一区| 在线观看一区二区三区| 亚洲熟妇熟女久久| 深爱激情五月婷婷| a在线观看视频网站| 别揉我奶头~嗯~啊~动态视频| 性欧美人与动物交配| 观看免费一级毛片| 国产成+人综合+亚洲专区| 免费在线观看影片大全网站| 人人妻,人人澡人人爽秒播| 97人妻精品一区二区三区麻豆| 精品乱码久久久久久99久播| 精品久久久久久久毛片微露脸| 欧美黑人巨大hd| 免费在线观看成人毛片| 中文在线观看免费www的网站| 欧美+日韩+精品| 国产极品精品免费视频能看的| 麻豆久久精品国产亚洲av| 日韩人妻高清精品专区| 国产一区在线观看成人免费| 久久婷婷人人爽人人干人人爱| 中文字幕人妻丝袜一区二区| a在线观看视频网站| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美在线乱码| 每晚都被弄得嗷嗷叫到高潮| 久久精品夜夜夜夜夜久久蜜豆| 久久久国产精品麻豆| www.www免费av| 51国产日韩欧美| 久久久精品大字幕| 熟妇人妻久久中文字幕3abv| 日韩欧美 国产精品| 日韩成人在线观看一区二区三区| 日本 av在线| 有码 亚洲区| 亚洲aⅴ乱码一区二区在线播放| 日本一二三区视频观看| 在线免费观看不下载黄p国产 | 国产亚洲精品久久久久久毛片| 日韩欧美 国产精品| av天堂在线播放| 九九久久精品国产亚洲av麻豆| 国产一区二区在线av高清观看| av视频在线观看入口| 国产精品嫩草影院av在线观看 | 色综合婷婷激情| 亚洲精品一卡2卡三卡4卡5卡| ponron亚洲| 国产一区在线观看成人免费| 一二三四社区在线视频社区8| 免费大片18禁| 十八禁人妻一区二区| 欧美色视频一区免费| 性欧美人与动物交配| 国产在线精品亚洲第一网站| 国产午夜精品论理片| 亚洲国产高清在线一区二区三| 午夜免费成人在线视频| 精品国产美女av久久久久小说| 天天添夜夜摸| 黄片大片在线免费观看| 美女高潮的动态| 观看免费一级毛片| 内射极品少妇av片p| 久久精品91蜜桃| 97人妻精品一区二区三区麻豆| 免费电影在线观看免费观看| 香蕉丝袜av| 国产精品嫩草影院av在线观看 | 91九色精品人成在线观看| 在线a可以看的网站| 男人舔奶头视频| 丰满的人妻完整版| 久久久久久大精品| 久久久久国内视频| 两人在一起打扑克的视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 免费看美女性在线毛片视频| 18美女黄网站色大片免费观看| 非洲黑人性xxxx精品又粗又长| 级片在线观看| 日韩欧美在线乱码| 成年女人看的毛片在线观看| 国产av在哪里看| 中文字幕人妻丝袜一区二区| 久久久久久久精品吃奶| 久99久视频精品免费| av黄色大香蕉| 亚洲专区国产一区二区| 久久亚洲精品不卡| 在线观看一区二区三区| 婷婷精品国产亚洲av在线| 美女高潮喷水抽搐中文字幕| 一夜夜www| 日本撒尿小便嘘嘘汇集6| 亚洲在线自拍视频| 婷婷亚洲欧美| 午夜免费成人在线视频| 中文字幕人妻丝袜一区二区| 一个人观看的视频www高清免费观看| 午夜激情欧美在线| 国产成人av激情在线播放| 欧美日韩精品网址| 久久国产乱子伦精品免费另类| 亚洲国产欧洲综合997久久,| 夜夜爽天天搞| 亚洲成a人片在线一区二区| 高清在线国产一区| 欧美日韩瑟瑟在线播放| 午夜福利在线在线| 免费高清视频大片| 国产精品久久久久久人妻精品电影| 日本撒尿小便嘘嘘汇集6| 国产精品电影一区二区三区| 国产男靠女视频免费网站| 日本免费a在线| 久久久色成人| 亚洲不卡免费看| 伊人久久精品亚洲午夜| 老熟妇仑乱视频hdxx| 99精品欧美一区二区三区四区| 男女视频在线观看网站免费| 性欧美人与动物交配| 天堂网av新在线| 一个人看的www免费观看视频| 91av网一区二区| 噜噜噜噜噜久久久久久91| 欧美日韩一级在线毛片| 亚洲av成人精品一区久久| 在线观看66精品国产| 美女免费视频网站| 乱人视频在线观看| 在线观看免费午夜福利视频| 国产高清三级在线| 麻豆成人av在线观看| 少妇丰满av| 国产成人福利小说| 色精品久久人妻99蜜桃| 日本黄色片子视频| 成人18禁在线播放| 久久精品国产亚洲av涩爱 | 1024手机看黄色片| 亚洲国产精品久久男人天堂| 99久久精品一区二区三区| 成人av一区二区三区在线看| 在线观看美女被高潮喷水网站 | 91av网一区二区| 国产高清视频在线观看网站| 色播亚洲综合网| 国产中年淑女户外野战色| 最新美女视频免费是黄的| 国产免费男女视频| 日韩人妻高清精品专区| 91麻豆精品激情在线观看国产| 国产黄色小视频在线观看| 午夜精品久久久久久毛片777| www.熟女人妻精品国产| 久久精品亚洲精品国产色婷小说| 日韩欧美 国产精品| 久久99热这里只有精品18| 日本a在线网址| 在线免费观看不下载黄p国产 | 国产爱豆传媒在线观看| 啦啦啦免费观看视频1| 五月玫瑰六月丁香| 午夜福利在线观看吧| 亚洲人成伊人成综合网2020| 国产色爽女视频免费观看| 网址你懂的国产日韩在线| 亚洲欧美日韩无卡精品| 亚洲无线在线观看| www.www免费av| 内地一区二区视频在线| 搞女人的毛片| 欧美又色又爽又黄视频| 男女床上黄色一级片免费看| 中亚洲国语对白在线视频| 久久精品国产99精品国产亚洲性色| 日韩大尺度精品在线看网址| 看黄色毛片网站| 老熟妇仑乱视频hdxx| 亚洲,欧美精品.| 欧美色视频一区免费| 日日夜夜操网爽| 日韩人妻高清精品专区| 久久国产乱子伦精品免费另类| 岛国视频午夜一区免费看| 一本久久中文字幕| 啦啦啦韩国在线观看视频| 91久久精品电影网| 舔av片在线| 精品99又大又爽又粗少妇毛片 | 欧美性猛交╳xxx乱大交人| 18禁黄网站禁片午夜丰满| 亚洲国产日韩欧美精品在线观看 | 麻豆国产av国片精品| 日日夜夜操网爽| 人人妻人人澡欧美一区二区| 国内毛片毛片毛片毛片毛片| 色尼玛亚洲综合影院| 精品一区二区三区人妻视频| 国产精品自产拍在线观看55亚洲| 精品久久久久久成人av| 亚洲精品在线美女| 亚洲无线在线观看| 成人亚洲精品av一区二区| 亚洲电影在线观看av| 精品久久久久久久毛片微露脸| 欧美一级毛片孕妇| 国产乱人伦免费视频| 精品国产三级普通话版| a在线观看视频网站| 脱女人内裤的视频| 村上凉子中文字幕在线| 亚洲五月天丁香| 真实男女啪啪啪动态图| 亚洲精品在线美女| 国产精品一区二区三区四区久久| 午夜激情欧美在线| 天天一区二区日本电影三级| 他把我摸到了高潮在线观看| 高清毛片免费观看视频网站| 啦啦啦观看免费观看视频高清| 久久欧美精品欧美久久欧美| 18禁在线播放成人免费| 又爽又黄无遮挡网站| 老司机在亚洲福利影院| 欧美黄色片欧美黄色片| 熟女人妻精品中文字幕| 国产精品野战在线观看| 日韩 欧美 亚洲 中文字幕| 美女 人体艺术 gogo| 婷婷精品国产亚洲av| 一个人免费在线观看的高清视频| www.熟女人妻精品国产| 精品电影一区二区在线| 天堂动漫精品| 少妇人妻精品综合一区二区 | 国产真实伦视频高清在线观看 | 国产探花极品一区二区| www日本黄色视频网| 桃红色精品国产亚洲av| 啦啦啦观看免费观看视频高清| 99国产极品粉嫩在线观看| 日韩av在线大香蕉| 亚洲在线自拍视频| 国产av在哪里看| 在线观看美女被高潮喷水网站 | 久久精品国产清高在天天线| 色精品久久人妻99蜜桃| 99精品在免费线老司机午夜| 国产午夜精品久久久久久一区二区三区 | 国产成人啪精品午夜网站| 亚洲中文字幕一区二区三区有码在线看| 久久久久久大精品| 熟女少妇亚洲综合色aaa.| 亚洲一区二区三区色噜噜| 舔av片在线| 丰满人妻熟妇乱又伦精品不卡| 91久久精品电影网| av视频在线观看入口| 99热这里只有是精品50| 好男人电影高清在线观看| 欧美激情在线99| 免费观看精品视频网站| 婷婷精品国产亚洲av在线| 亚洲五月天丁香| a级一级毛片免费在线观看| 成人三级黄色视频| 岛国在线免费视频观看| 一级黄片播放器| av天堂中文字幕网| 国产国拍精品亚洲av在线观看 | 可以在线观看毛片的网站| 97人妻精品一区二区三区麻豆| 男人的好看免费观看在线视频| 91九色精品人成在线观看| 国产欧美日韩一区二区三| 女警被强在线播放| 制服丝袜大香蕉在线| 内地一区二区视频在线| 日韩免费av在线播放| 男女视频在线观看网站免费| 国产伦精品一区二区三区四那| 一a级毛片在线观看| 成人性生交大片免费视频hd| 欧美激情久久久久久爽电影| 变态另类丝袜制服| 欧美中文综合在线视频| 国产三级中文精品| 最近最新免费中文字幕在线| a级一级毛片免费在线观看| 久久精品夜夜夜夜夜久久蜜豆| 女警被强在线播放| 一级黄色大片毛片| 午夜亚洲福利在线播放| 亚洲欧美日韩高清在线视频| 搡老岳熟女国产| 欧美日韩一级在线毛片| 国产精品99久久99久久久不卡| 久久精品综合一区二区三区| 亚洲第一电影网av| 亚洲精品456在线播放app | 精品久久久久久久末码| 99久久99久久久精品蜜桃| 在线免费观看的www视频| 成年女人毛片免费观看观看9| 人妻夜夜爽99麻豆av| 嫩草影院入口| 9191精品国产免费久久| 99热只有精品国产| 欧美大码av| 在线视频色国产色| 给我免费播放毛片高清在线观看| 搡老熟女国产l中国老女人| 国产精品久久久人人做人人爽| bbb黄色大片| 一个人看的www免费观看视频| 国产三级黄色录像| 亚洲av不卡在线观看| 1000部很黄的大片| av视频在线观看入口| 脱女人内裤的视频| www日本黄色视频网| 麻豆成人午夜福利视频| 欧美日韩瑟瑟在线播放| 变态另类丝袜制服| 亚洲国产中文字幕在线视频| 日韩免费av在线播放| 亚洲精品影视一区二区三区av| 小蜜桃在线观看免费完整版高清| 热99在线观看视频| 非洲黑人性xxxx精品又粗又长| 三级男女做爰猛烈吃奶摸视频| 国产亚洲欧美在线一区二区| 久久伊人香网站| 99国产极品粉嫩在线观看| 99riav亚洲国产免费| 欧美+日韩+精品| 国产欧美日韩精品亚洲av| 欧美日韩亚洲国产一区二区在线观看| 一区二区三区高清视频在线| 99久久无色码亚洲精品果冻| 丁香六月欧美| 香蕉av资源在线| 成人鲁丝片一二三区免费| 我的老师免费观看完整版| 啪啪无遮挡十八禁网站| 国产午夜精品论理片| 久久久精品欧美日韩精品| 美女被艹到高潮喷水动态| 免费高清视频大片| 中文字幕久久专区| 宅男免费午夜| 精品午夜福利视频在线观看一区| 神马国产精品三级电影在线观看| 最近在线观看免费完整版| 在线观看免费视频日本深夜| 99国产精品一区二区蜜桃av| 亚洲国产日韩欧美精品在线观看 | 精品人妻一区二区三区麻豆 | av在线蜜桃| 午夜免费激情av| 在线观看66精品国产| 欧美在线一区亚洲| 91九色精品人成在线观看| 国产一区二区亚洲精品在线观看| 国产在线精品亚洲第一网站| 国产精品久久久人人做人人爽| 欧美成人免费av一区二区三区| 欧美日本亚洲视频在线播放| 美女 人体艺术 gogo| 91久久精品电影网| 成年人黄色毛片网站| 成人特级黄色片久久久久久久| 最近最新免费中文字幕在线| 757午夜福利合集在线观看| 在线观看美女被高潮喷水网站 | 国产美女午夜福利| 91久久精品电影网| 99精品久久久久人妻精品| 一区二区三区高清视频在线| 最近在线观看免费完整版| 麻豆久久精品国产亚洲av| 亚洲av免费高清在线观看| 欧美三级亚洲精品| 欧美乱妇无乱码| 国产极品精品免费视频能看的| 免费观看的影片在线观看| 熟女电影av网| 男插女下体视频免费在线播放| 中文亚洲av片在线观看爽| 一边摸一边抽搐一进一小说| 一级黄色大片毛片| 久久天躁狠狠躁夜夜2o2o| 一级黄色大片毛片| 一边摸一边抽搐一进一小说| 又紧又爽又黄一区二区| 一边摸一边抽搐一进一小说| 国产老妇女一区| 18禁裸乳无遮挡免费网站照片| 日本成人三级电影网站| 日韩高清综合在线| 亚洲成人免费电影在线观看| 美女cb高潮喷水在线观看| 欧美日韩黄片免| 美女cb高潮喷水在线观看| 国产极品精品免费视频能看的| 欧美成狂野欧美在线观看| 999久久久精品免费观看国产| 给我免费播放毛片高清在线观看| 亚洲成人久久性| 中文字幕精品亚洲无线码一区| 午夜影院日韩av| 国产精品 欧美亚洲| 岛国视频午夜一区免费看| 国产欧美日韩精品一区二区| 亚洲av电影在线进入| 免费看光身美女| 亚洲黑人精品在线| 熟女电影av网| 好男人在线观看高清免费视频| 欧美一区二区精品小视频在线| 少妇人妻精品综合一区二区 | 一个人免费在线观看的高清视频| 午夜免费男女啪啪视频观看 | 国产精品女同一区二区软件 | 91在线精品国自产拍蜜月 | 香蕉丝袜av| 国产毛片a区久久久久| 色精品久久人妻99蜜桃| 国产精品久久久久久久电影 | av福利片在线观看| 久久久久九九精品影院| 不卡一级毛片| 一进一出抽搐动态| 神马国产精品三级电影在线观看| 桃色一区二区三区在线观看| 18禁黄网站禁片免费观看直播|