• <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
    搡老乐熟女国产| 欧美日韩在线观看h| 国产精品福利在线免费观看| 欧美日韩视频精品一区| 最新中文字幕久久久久| 女人精品久久久久毛片| 久久久久久伊人网av| 色吧在线观看| 久久人人爽人人爽人人片va| 亚洲av电影在线观看一区二区三区| 在线观看国产h片| 日日爽夜夜爽网站| 日本爱情动作片www.在线观看| 国产日韩一区二区三区精品不卡 | 91成人精品电影| 男人狂女人下面高潮的视频| 国产无遮挡羞羞视频在线观看| 欧美丝袜亚洲另类| 啦啦啦中文免费视频观看日本| 精品久久久久久久久av| 少妇的逼好多水| 久久影院123| 黄色毛片三级朝国网站 | www.av在线官网国产| 亚洲成人av在线免费| 亚洲欧美精品自产自拍| 边亲边吃奶的免费视频| 亚洲精品aⅴ在线观看| 亚洲一区二区三区欧美精品| 中文欧美无线码| 日本vs欧美在线观看视频 | 丝袜喷水一区| 两个人免费观看高清视频 | 在线观看美女被高潮喷水网站| 日韩成人av中文字幕在线观看| 嫩草影院新地址| 一级毛片我不卡| 简卡轻食公司| 亚洲人与动物交配视频| 色视频www国产| 晚上一个人看的免费电影| 国产精品蜜桃在线观看| 蜜桃久久精品国产亚洲av| 国产亚洲av片在线观看秒播厂| 成人综合一区亚洲| 国产欧美日韩精品一区二区| 国产欧美日韩精品一区二区| 少妇的逼水好多| 日韩欧美一区视频在线观看 | 黑人猛操日本美女一级片| av又黄又爽大尺度在线免费看| 高清av免费在线| 丰满迷人的少妇在线观看| 高清不卡的av网站| 最近手机中文字幕大全| 男人爽女人下面视频在线观看| a 毛片基地| 超碰97精品在线观看| 久久久久久久大尺度免费视频| 99热这里只有是精品在线观看| 麻豆乱淫一区二区| 国产精品一区www在线观看| 亚洲av综合色区一区| 国产精品一二三区在线看| 人妻系列 视频| 久久久久久久大尺度免费视频| 97在线人人人人妻| 青春草视频在线免费观看| 亚洲在久久综合| 国产成人精品久久久久久| 国产精品偷伦视频观看了| 极品人妻少妇av视频| 噜噜噜噜噜久久久久久91| 国产高清有码在线观看视频| 大码成人一级视频| 久久久国产一区二区| 久久ye,这里只有精品| 一本久久精品| 国产日韩欧美视频二区| 美女视频免费永久观看网站| 极品教师在线视频| 高清毛片免费看| 午夜影院在线不卡| 色5月婷婷丁香| 女性生殖器流出的白浆| 国产精品一区二区性色av| a级毛片在线看网站| 免费高清在线观看视频在线观看| 日韩伦理黄色片| 在线观看三级黄色| 天美传媒精品一区二区| 亚洲高清免费不卡视频| 中文乱码字字幕精品一区二区三区| 国产老妇伦熟女老妇高清| 久久99热6这里只有精品| av在线app专区| av视频免费观看在线观看| 成人免费观看视频高清| 亚洲成人av在线免费| 亚洲精品日韩在线中文字幕| 免费人妻精品一区二区三区视频| 成人国产麻豆网| 国产亚洲午夜精品一区二区久久| 久久人人爽av亚洲精品天堂| 最近手机中文字幕大全| 午夜激情福利司机影院| 蜜桃久久精品国产亚洲av| 久久国产精品男人的天堂亚洲 | 激情五月婷婷亚洲| 久久 成人 亚洲| 精品人妻熟女av久视频| 国产欧美日韩一区二区三区在线 | 丝袜喷水一区| 18+在线观看网站| 成人亚洲精品一区在线观看| 国产熟女欧美一区二区| 午夜av观看不卡| 2018国产大陆天天弄谢| 国产高清国产精品国产三级| 九色成人免费人妻av| 亚洲av.av天堂| 18禁在线播放成人免费| 亚洲婷婷狠狠爱综合网| 亚洲欧美一区二区三区国产| 婷婷色av中文字幕| 中文字幕人妻熟人妻熟丝袜美| 中国美白少妇内射xxxbb| 国产av精品麻豆| 在线看a的网站| 91aial.com中文字幕在线观看| 18禁裸乳无遮挡动漫免费视频| 一级毛片久久久久久久久女| 最后的刺客免费高清国语| 在线免费观看不下载黄p国产| 老女人水多毛片| 午夜福利,免费看| 精品熟女少妇av免费看| 欧美精品人与动牲交sv欧美| 国产精品秋霞免费鲁丝片| 成人特级av手机在线观看| 日日摸夜夜添夜夜爱| 成年人午夜在线观看视频| 久久久久久久国产电影| 欧美变态另类bdsm刘玥| 九草在线视频观看| 久久久亚洲精品成人影院| 91aial.com中文字幕在线观看| 亚洲精品乱久久久久久| 国产综合精华液| 在线 av 中文字幕| 国产片特级美女逼逼视频| 日韩一区二区三区影片| 国产在线一区二区三区精| 久久亚洲国产成人精品v| 精品国产露脸久久av麻豆| 男女啪啪激烈高潮av片| 欧美日韩视频精品一区| 国产精品免费大片| 天堂8中文在线网| 热99国产精品久久久久久7| 久久久久久人妻| 久久亚洲国产成人精品v| 免费观看av网站的网址| 亚洲天堂av无毛| 精品久久国产蜜桃| 久久6这里有精品| 国产精品一二三区在线看| 日韩视频在线欧美| 亚洲国产精品999| 午夜福利在线观看免费完整高清在| 欧美精品一区二区大全| av福利片在线| 成人影院久久| a 毛片基地| 在线观看一区二区三区激情| 久久鲁丝午夜福利片| 汤姆久久久久久久影院中文字幕| 亚洲伊人久久精品综合| 人人澡人人妻人| 亚洲欧美中文字幕日韩二区| 免费黄网站久久成人精品| 成人午夜精彩视频在线观看| 你懂的网址亚洲精品在线观看| 亚洲精品日韩av片在线观看| 久久国产精品男人的天堂亚洲 | 噜噜噜噜噜久久久久久91| 精品卡一卡二卡四卡免费| 国产视频内射| 亚洲人成网站在线观看播放| 欧美日本中文国产一区发布| 国内少妇人妻偷人精品xxx网站| 色婷婷av一区二区三区视频| 久久人人爽av亚洲精品天堂| 欧美xxxx性猛交bbbb| 男女无遮挡免费网站观看| 欧美最新免费一区二区三区| 精华霜和精华液先用哪个| 夫妻午夜视频| 人妻人人澡人人爽人人| 国产精品久久久久久久久免| 国产欧美日韩精品一区二区| 亚洲国产精品一区二区三区在线| 精品久久久久久电影网| 国产高清有码在线观看视频| 日韩一区二区三区影片| 一级毛片aaaaaa免费看小| 在线天堂最新版资源| 国产精品无大码| 午夜久久久在线观看| 日本av手机在线免费观看| 中文乱码字字幕精品一区二区三区| 国产精品伦人一区二区| 亚洲欧美成人精品一区二区| 夜夜看夜夜爽夜夜摸| 日韩一本色道免费dvd| 中文精品一卡2卡3卡4更新| 亚洲精品国产av成人精品| 久久精品国产亚洲网站| 免费少妇av软件| 黄色配什么色好看| 男人舔奶头视频| 一级毛片 在线播放| 中国国产av一级| 成人无遮挡网站| 九九爱精品视频在线观看| 国产毛片在线视频| 婷婷色麻豆天堂久久| 国产精品人妻久久久影院| 亚洲成人手机| 一级av片app| 中文资源天堂在线| 亚洲精品aⅴ在线观看| 国产男人的电影天堂91| 国产亚洲精品久久久com| 美女中出高潮动态图| 欧美日韩综合久久久久久| 街头女战士在线观看网站| 亚洲欧洲精品一区二区精品久久久 | 99热网站在线观看| 日本猛色少妇xxxxx猛交久久| 18禁在线无遮挡免费观看视频| 亚洲欧美日韩另类电影网站| 精品亚洲成a人片在线观看| 亚洲国产精品999| 伊人亚洲综合成人网| 亚州av有码| 国产成人a∨麻豆精品| 天堂俺去俺来也www色官网| av在线app专区| 高清视频免费观看一区二区| 日韩欧美 国产精品| 国产成人91sexporn| 久久综合国产亚洲精品| 国产中年淑女户外野战色| 国产精品人妻久久久久久| 男女边摸边吃奶| 一本色道久久久久久精品综合| 免费观看性生交大片5| 国产精品久久久久久精品古装| 麻豆成人av视频| 91午夜精品亚洲一区二区三区| 国产国拍精品亚洲av在线观看| 99九九在线精品视频 | 国产av精品麻豆| 成人国产麻豆网| 国产黄片视频在线免费观看| 一级黄片播放器| 在线观看免费视频网站a站| 日韩熟女老妇一区二区性免费视频| av福利片在线| 校园人妻丝袜中文字幕| 狂野欧美激情性bbbbbb| 精品一区二区三卡| 大陆偷拍与自拍| 观看av在线不卡| 午夜免费鲁丝| 国产日韩欧美视频二区| 日韩成人伦理影院| 国产极品天堂在线| 国产男人的电影天堂91| 婷婷色综合www| videos熟女内射| 久久鲁丝午夜福利片| 美女主播在线视频| 一级毛片久久久久久久久女| 中文字幕制服av| 99国产精品免费福利视频| 只有这里有精品99| 嫩草影院新地址| 亚洲精品第二区| 另类亚洲欧美激情| 在线观看av片永久免费下载| 国产亚洲91精品色在线| 国语对白做爰xxxⅹ性视频网站| 国产一区有黄有色的免费视频| 国产有黄有色有爽视频| 亚州av有码| 久久久久网色| 亚洲三级黄色毛片| 成人毛片60女人毛片免费| 欧美精品国产亚洲| 欧美xxxx性猛交bbbb| 在线看a的网站| 国产精品人妻久久久影院| 简卡轻食公司| 亚洲人成网站在线观看播放| 91久久精品电影网| 激情五月婷婷亚洲| 久久久久国产精品人妻一区二区| 99久久精品一区二区三区| 亚洲av中文av极速乱| 婷婷色综合大香蕉| 人人妻人人看人人澡| 青春草国产在线视频| 国精品久久久久久国模美| 国产永久视频网站| 一个人免费看片子| 高清欧美精品videossex| 成人二区视频| 内射极品少妇av片p| 国产 一区精品| 六月丁香七月| 精品一区二区三卡| 日产精品乱码卡一卡2卡三| 9色porny在线观看| 美女视频免费永久观看网站| 欧美日本中文国产一区发布| 国产成人91sexporn| 中文字幕精品免费在线观看视频 | 精品国产一区二区久久| 夜夜骑夜夜射夜夜干| 精品一品国产午夜福利视频| 男女边摸边吃奶| 久久人人爽av亚洲精品天堂| 中文字幕免费在线视频6| 99精国产麻豆久久婷婷| 各种免费的搞黄视频| 国产精品国产三级国产专区5o| 精品亚洲成a人片在线观看| 午夜福利视频精品| 中文字幕人妻熟人妻熟丝袜美| 高清午夜精品一区二区三区| 成人国产av品久久久| a级毛片免费高清观看在线播放| 久久人人爽人人片av| 亚洲欧美日韩东京热| 亚洲欧洲国产日韩| 亚洲精品国产av成人精品| 国产在线男女| 久久精品久久久久久噜噜老黄| 一本大道久久a久久精品| av专区在线播放| 黄片无遮挡物在线观看| 久久国产亚洲av麻豆专区| 一级爰片在线观看| 久久久久久久久久久久大奶| 亚洲精品国产成人久久av| 香蕉精品网在线| 你懂的网址亚洲精品在线观看| 国产亚洲午夜精品一区二区久久| 妹子高潮喷水视频| 岛国毛片在线播放| 亚州av有码| 欧美丝袜亚洲另类| 99精国产麻豆久久婷婷| 九色成人免费人妻av| 国产亚洲一区二区精品| 成人亚洲欧美一区二区av| 日韩av免费高清视频| 丝袜在线中文字幕| 久久人妻熟女aⅴ| 免费久久久久久久精品成人欧美视频 | 三级国产精品片| 一级黄片播放器| 午夜免费观看性视频| 国产精品福利在线免费观看| 色视频在线一区二区三区| 国产日韩欧美亚洲二区| 日韩一本色道免费dvd| 久久99热这里只频精品6学生| 久久精品国产a三级三级三级| 亚洲欧美一区二区三区黑人 | 亚洲,一卡二卡三卡| 国产免费又黄又爽又色| 日韩av在线免费看完整版不卡| 日日摸夜夜添夜夜爱| 精品亚洲乱码少妇综合久久| 成人毛片a级毛片在线播放| 久久ye,这里只有精品| 99久久人妻综合| 国产黄色免费在线视频| 亚洲色图综合在线观看| 精品亚洲成a人片在线观看| 男人舔奶头视频| 日日爽夜夜爽网站| 在线观看免费视频网站a站| 精华霜和精华液先用哪个| 国产精品嫩草影院av在线观看| 在线精品无人区一区二区三| 亚洲av.av天堂| 丰满饥渴人妻一区二区三| 亚洲成色77777| 欧美少妇被猛烈插入视频| 特大巨黑吊av在线直播| 天天躁夜夜躁狠狠久久av| 人人妻人人添人人爽欧美一区卜| 欧美精品高潮呻吟av久久| 汤姆久久久久久久影院中文字幕| 99久久精品国产国产毛片| 看免费成人av毛片| xxx大片免费视频| 色视频www国产| 午夜福利网站1000一区二区三区| 亚洲色图综合在线观看| 久久婷婷青草| 搡女人真爽免费视频火全软件| 国产美女午夜福利| 成人亚洲精品一区在线观看| 又粗又硬又长又爽又黄的视频| 91精品国产九色| 成人毛片a级毛片在线播放| 久久国产亚洲av麻豆专区| 国产午夜精品一二区理论片| 美女国产视频在线观看| 麻豆乱淫一区二区| 肉色欧美久久久久久久蜜桃| 国产精品一区二区在线不卡| av黄色大香蕉| 欧美日韩在线观看h| 国产老妇伦熟女老妇高清| 人妻制服诱惑在线中文字幕| 亚洲精品国产成人久久av| 成人美女网站在线观看视频| 最新中文字幕久久久久| 51国产日韩欧美| 全区人妻精品视频| 久久久久久久精品精品| 精品视频人人做人人爽| 亚州av有码| 久久人人爽人人片av| 欧美区成人在线视频| 99九九线精品视频在线观看视频| 免费少妇av软件| 又大又黄又爽视频免费| 国产 精品1| 三级经典国产精品| 欧美精品国产亚洲| 女的被弄到高潮叫床怎么办| 美女大奶头黄色视频| 美女内射精品一级片tv| 国产成人一区二区在线| 国产男女内射视频| 久久精品国产a三级三级三级| 亚洲真实伦在线观看| 99久国产av精品国产电影| 国产成人精品福利久久| 伦理电影免费视频| av在线app专区| 亚洲美女黄色视频免费看| 我要看黄色一级片免费的| 在线精品无人区一区二区三| 国产男人的电影天堂91| 国产在线一区二区三区精| 国产精品国产三级专区第一集| 日韩免费高清中文字幕av| 久久精品久久精品一区二区三区| 另类亚洲欧美激情| 中国三级夫妇交换| 亚洲av不卡在线观看| 街头女战士在线观看网站| 我要看黄色一级片免费的| 黄色配什么色好看| 国产精品人妻久久久影院| 亚洲国产成人一精品久久久| 你懂的网址亚洲精品在线观看| 国产黄片视频在线免费观看| 日韩在线高清观看一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 国产精品久久久久久久久免| 亚洲精品视频女| 我要看黄色一级片免费的| 性高湖久久久久久久久免费观看| 国产精品女同一区二区软件| 国产av一区二区精品久久| 亚洲精品乱久久久久久| 午夜精品国产一区二区电影| 日韩av不卡免费在线播放| 国产精品熟女久久久久浪| h日本视频在线播放| 免费观看性生交大片5| 精品久久久噜噜| 亚洲精品乱久久久久久| 精品99又大又爽又粗少妇毛片| 深夜a级毛片| 国产日韩欧美视频二区| 日韩一本色道免费dvd| 国产亚洲精品久久久com| 国产男女超爽视频在线观看| 久久 成人 亚洲| 欧美日本中文国产一区发布| 在线观看国产h片| 国产伦精品一区二区三区四那| xxx大片免费视频| 一级爰片在线观看| av免费在线看不卡| 久久精品国产鲁丝片午夜精品| 国产欧美日韩综合在线一区二区 | 97在线人人人人妻| 26uuu在线亚洲综合色| 美女视频免费永久观看网站| 国产精品偷伦视频观看了| 51国产日韩欧美| 纵有疾风起免费观看全集完整版| 我要看黄色一级片免费的| 久久国产精品大桥未久av | 不卡视频在线观看欧美| 激情五月婷婷亚洲| 精品视频人人做人人爽| av又黄又爽大尺度在线免费看| 99久久人妻综合| 精品国产乱码久久久久久小说| 秋霞在线观看毛片| 天天操日日干夜夜撸| 99久国产av精品国产电影| 国产熟女欧美一区二区| 王馨瑶露胸无遮挡在线观看| 欧美+日韩+精品| 欧美日韩视频高清一区二区三区二| 蜜桃在线观看..| 一级毛片电影观看| 黄色欧美视频在线观看| 精品午夜福利在线看| 乱码一卡2卡4卡精品| 亚洲怡红院男人天堂| 成年av动漫网址| 91aial.com中文字幕在线观看| 少妇熟女欧美另类| 国产永久视频网站| 丝袜脚勾引网站| 久久人妻熟女aⅴ| 国产亚洲精品久久久com| 人体艺术视频欧美日本| 久久韩国三级中文字幕| 亚洲国产日韩一区二区| 午夜免费鲁丝| 国产探花极品一区二区| 欧美日韩视频高清一区二区三区二| 国产精品不卡视频一区二区| 成人影院久久| 伦精品一区二区三区| 观看美女的网站| www.av在线官网国产| 中文字幕av电影在线播放| 久久av网站| 国产黄片视频在线免费观看| 亚洲在久久综合| 国产高清有码在线观看视频| 九草在线视频观看| 午夜福利视频精品| 久久久久久久久久人人人人人人| 99九九在线精品视频 | 六月丁香七月| 老司机影院毛片| 春色校园在线视频观看| 大话2 男鬼变身卡| 在线 av 中文字幕| 伊人亚洲综合成人网| 免费黄色在线免费观看| 在线观看免费日韩欧美大片 | 性色av一级| 美女主播在线视频| 精品99又大又爽又粗少妇毛片| 亚洲av成人精品一二三区| 精品亚洲成a人片在线观看| 亚洲精品乱久久久久久| 国产成人精品一,二区| 久久久国产一区二区| 午夜免费观看性视频| 国产视频首页在线观看| 亚洲性久久影院| 日韩成人av中文字幕在线观看| 国产一区二区三区综合在线观看 | 全区人妻精品视频| 在线 av 中文字幕| 日韩一本色道免费dvd| 天堂俺去俺来也www色官网| 日本与韩国留学比较| 三上悠亚av全集在线观看 | 晚上一个人看的免费电影| 久久精品久久久久久噜噜老黄| 亚洲精品一区蜜桃| 极品人妻少妇av视频| 日本欧美视频一区| 男女国产视频网站| 熟女av电影| 视频中文字幕在线观看| 黄色一级大片看看| 国语对白做爰xxxⅹ性视频网站| 在线看a的网站| 在线免费观看不下载黄p国产| 欧美精品一区二区大全| 少妇熟女欧美另类| 日韩伦理黄色片| 极品少妇高潮喷水抽搐| 精品国产一区二区久久| 又黄又爽又刺激的免费视频.| 国产精品三级大全| 午夜福利,免费看| 六月丁香七月| 夜夜看夜夜爽夜夜摸| 黄色视频在线播放观看不卡| 久久99精品国语久久久| 中文字幕制服av| 十分钟在线观看高清视频www | 欧美精品人与动牲交sv欧美|