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

    二次有機(jī)氣溶膠估算方法研究進(jìn)展

    2014-05-02 11:03:14閆才青李小瀅王雪松張遠(yuǎn)航
    中國(guó)環(huán)境科學(xué) 2014年3期
    關(guān)鍵詞:氣溶膠貢獻(xiàn)生物質(zhì)

    鄭 玫,閆才青,李小瀅,王雪松,張遠(yuǎn)航

    (北京大學(xué)環(huán)境科學(xué)與工程學(xué)院,環(huán)境模擬與污染控制國(guó)家重點(diǎn)聯(lián)合實(shí)驗(yàn)室,北京 100871)

    二次有機(jī)氣溶膠估算方法研究進(jìn)展

    鄭 玫,閆才青,李小瀅,王雪松,張遠(yuǎn)航*

    (北京大學(xué)環(huán)境科學(xué)與工程學(xué)院,環(huán)境模擬與污染控制國(guó)家重點(diǎn)聯(lián)合實(shí)驗(yàn)室,北京 100871)

    我國(guó)區(qū)域復(fù)合污染日趨嚴(yán)重,二次有機(jī)氣溶膠(SOA)是細(xì)顆粒物PM2.5的重要組成之一.本研究總結(jié)了當(dāng)前國(guó)內(nèi)外定量估算SOA的主要方法,包括EC示蹤法、WSOC法、受體模型法、示蹤物產(chǎn)率法、數(shù)值模擬法等.通過(guò)對(duì)各種方法的原理介紹與比較,指出1)EC示蹤法、WSOC法、受體模型法簡(jiǎn)單易行,與在線連續(xù)觀測(cè)數(shù)據(jù)結(jié)合可估算高時(shí)間分辨率的 SOA濃度,但受限于對(duì)當(dāng)?shù)卦醋V和特定物種的了解;2)示蹤物產(chǎn)率法雖分析技術(shù)難度高,但可針對(duì)不同來(lái)源的SOA估算;3)數(shù)值模擬可獲得大尺度SOA的空間分布.針對(duì)國(guó)內(nèi)外的最新研究成果進(jìn)行簡(jiǎn)單概述,不同方法的研究表明中國(guó)二次有機(jī)氣溶膠是總有機(jī)氣溶膠的重要部分,人為源VOCs對(duì)二次有機(jī)氣溶膠發(fā)揮著重要貢獻(xiàn).本文旨在為未來(lái)我國(guó)二次有機(jī)氣溶膠的研究提供基礎(chǔ)信息和研究思路.

    二次有機(jī)氣溶膠;細(xì)顆粒物;定量估算方法;來(lái)源;中國(guó)

    有機(jī)氣溶膠(OA)包括污染源直接排放的一次有機(jī)氣溶膠(POA),以及二次轉(zhuǎn)化形成的二次有機(jī)氣溶膠(SOA).SOA是由天然源VOCs和人為源 VOCs氣態(tài)前體物與大氣中的氧化劑(如O3、OH、NO3)發(fā)生光化學(xué)反應(yīng)后,生成的半揮發(fā)和難揮發(fā)性有機(jī)物,經(jīng)過(guò)氣-粒轉(zhuǎn)化或凝結(jié)到已有顆粒物上形成的有機(jī)氣溶膠[1-2].

    研究表明,SOA是大氣氣溶膠的重要組成部分.在美國(guó)東南部和加拿大,大約 10%~80%的有機(jī)氣溶膠為二次來(lái)源[3-4];在中國(guó),顆粒物中二次來(lái)源的有機(jī)物約占總有機(jī)物的 30%~95%[5-7].由于含氧、含氮等極性官能團(tuán)的引入,二次有機(jī)氣溶膠具有更強(qiáng)的極性和吸濕性,對(duì)能見(jiàn)度降低、灰霾形成、氣候變化具有重要影響.

    然而,由于對(duì)SOA復(fù)雜的大氣形成過(guò)程、凝結(jié)/分配機(jī)制以及化學(xué)組成尚缺乏全面的認(rèn)識(shí),目前還不具備統(tǒng)一的對(duì) SOA直接測(cè)量的分析手段,SOA的識(shí)別及其來(lái)源貢獻(xiàn)的研究仍是國(guó)際氣溶膠研究領(lǐng)域的熱點(diǎn)和難點(diǎn)[8-9].除利用煙霧箱在特定的條件下直接模擬SOA的生成外,大氣環(huán)境中SOA濃度的估算一般均采用間接方法.目前常用的SOA估算方法大致可概括為兩大類,即基于觀測(cè)數(shù)據(jù)的估算與基于源排放清單的模擬.前者主要包括EC示蹤法、WSOC法、受體模型法、示蹤物-產(chǎn)率法;而后者主要指數(shù)值模擬法.本文綜述了國(guó)內(nèi)外SOA定量估算的主要方法,總結(jié)國(guó)內(nèi)SOA的最新研究結(jié)果,以期為評(píng)估我國(guó)SOA對(duì)細(xì)顆粒物的貢獻(xiàn)以及相關(guān)控制對(duì)策的制定提供依據(jù).

    1 基于觀測(cè)數(shù)據(jù)的估算法

    1.1 EC示蹤法

    元素碳EC來(lái)自化石燃料和生物質(zhì)等含碳燃料不完全燃燒的直接排放;一次燃燒排放 EC的同時(shí)也伴隨產(chǎn)生有機(jī)碳OC.因此,EC可作為指示燃燒排放OC的示蹤物種[3,10],根據(jù)EC可估算一次有機(jī)碳(POC),二次有機(jī)碳(SOC)則可經(jīng)由式(1)、(2)計(jì)算得到[11]∶

    其中,OCpri, OCnc, (OC/EC)pri, OCsec, OCtot分別代表一次來(lái)源的有機(jī)碳、一次非燃燒源有機(jī)碳、一次燃燒源排放的OC與EC比值、二次來(lái)源的有機(jī)碳和測(cè)量的總有機(jī)碳.EC與OCtot可通過(guò)野外在線觀測(cè)或?qū)嶒?yàn)室離線測(cè)量獲取,常用儀器包括美國(guó)Sunset實(shí)驗(yàn)室的半連續(xù)或離線OC/EC分析儀(Sunset Laboratory Inc., Model 3型)、以及美國(guó)沙漠研究所研制的 OC/EC分析儀(DRI Model 2001型).

    EC示蹤法亦稱為排放清單OC/EC比值法或最小OC/EC比值法[12],該法的關(guān)鍵在于OCnc與(OC/EC)pri的估算. (OC/EC)pri的獲取主要通過(guò)以下幾種方式∶(1)各種一次源排放清單中OC/EC比值的平均值[13];(2)將光化學(xué)反應(yīng)較弱(光照強(qiáng)度弱、臭氧濃度低)、一次源排放占主導(dǎo)的相對(duì)穩(wěn)定的大氣環(huán)境下,環(huán)境樣品的OC、EC數(shù)據(jù)做線性相關(guān)分析[11],回歸曲線的截距即為OCnc,斜率則被視為(OC/EC)pri[14],OCnc通常被忽略不計(jì)[11];(3)觀測(cè)期間環(huán)境OC與EC最小觀測(cè)比值[15],即利用統(tǒng)計(jì)學(xué)方法,對(duì)觀測(cè)所得的OC/EC比值做頻率分布,選擇一定置信區(qū)間(一般為5%或10%)的最小比值作為一次源排放的OC/EC比值[16].

    EC示蹤法由于僅需OC與EC數(shù)據(jù),方法簡(jiǎn)單而被廣泛應(yīng)用[16-24].Cao等[18]在全國(guó)14個(gè)城市站點(diǎn)進(jìn)行冬、夏季同步觀測(cè),利用EC示蹤法估算表明中國(guó)城市大氣環(huán)境中 SOC約占總碳的23.7%~41.5%; Zhang等[19]對(duì)中國(guó)18個(gè)城市、鄉(xiāng)村、遠(yuǎn)郊站點(diǎn)長(zhǎng)達(dá)一年的 OC、EC數(shù)據(jù)分析表明,城市與鄉(xiāng)村大氣中分別有 48%~62%和53%~83%的OC來(lái)自于二次轉(zhuǎn)化.

    基于EC示蹤法開(kāi)展的中國(guó)SOC研究(見(jiàn)圖1)表明∶在空間分布上,北京[5]、天津[20]、濟(jì)南[21]、廈門[22]等地 SOC濃度明顯高于珠三角[23]、福建[24]地區(qū),內(nèi)陸城市高于沿海城市,北方高于南方;從季節(jié)變化上,各大城市秋、冬季SOC濃度明顯高于春、夏季,且冬季最高,這可能是由于冬季供暖燃煤等人為源排放增加,半揮發(fā)性及揮發(fā)性SOA前體物排放增加[5],且擴(kuò)散條件差,降水量少,濕去除不明顯,較低的溫度也有利于氣-粒轉(zhuǎn)換的發(fā)生[21].此外,污染物的區(qū)域傳輸[23]和氣象條件(尤其是風(fēng)向[22-23])均會(huì)影響二次有機(jī)污染物的分布.

    EC示蹤法適用于污染源種類和源排放量相對(duì)穩(wěn)定的條件.然而,實(shí)際環(huán)境大氣常處于動(dòng)態(tài)變化中,即使取最小的OC/EC比值仍難以避免二次有機(jī)物的干擾[12];同時(shí),氣象條件、長(zhǎng)距離傳輸?shù)榷紩?huì)影響OC/EC比值[22]. 借助高時(shí)間分辨率的CO數(shù)據(jù)對(duì)OC進(jìn)行標(biāo)準(zhǔn)化,可減少因新鮮排放、稀釋、大氣混合等導(dǎo)致的OC變化[25-26]. Hu等[27]利用統(tǒng)計(jì)方法,先假定OC/EC比值范圍,以0.001為步長(zhǎng)選擇不同 OC/EC比值作為(OC/EC)pri,分別用EC示蹤法求解,獲得的各組SOC與實(shí)測(cè)EC相關(guān)系數(shù)最小時(shí)對(duì)應(yīng)的 OC/EC比值被視為(OC/EC)pri;Yu等[28]利用排放/傳輸模型(如多尺度綜合空氣質(zhì)量模型,CMAQ),結(jié)合高時(shí)間分辨率排放清單,追蹤模擬(OC/EC)pri,并與測(cè)量獲取的環(huán)境OC、EC數(shù)據(jù)結(jié)合,為區(qū)分POC和SOC提供了新思路.此外,不同OC、EC分析方法(熱光反射法 TOR和熱光透射法 TOT)會(huì)導(dǎo)致測(cè)量的OC/EC比值不同[29],增加估算結(jié)果的不確定性.利用環(huán)境OC與EC濃度估算SOC時(shí)需評(píng)估分析方法的差異,盡量獲取接近一次燃燒源的OC/EC比值.

    圖1 中國(guó)各大城市基于EC示蹤法的SOC濃度分布[5,20-24]Fig.1 SOC mass concentration in different cities of China based on the EC tracer method[5,20-24]

    在 EC示蹤法的基礎(chǔ)上,一些研究者利用其他物種(如CO)代替EC作為一次燃燒源示蹤物,或同時(shí)采用多個(gè)物種加以限制(如多元回歸法)[25-26,30],引申了 EC示蹤法.如 Blanchard等[26]利用CO和EC作為POC指示物,SO42-、NO3-和臭氧作為 SOC指示物,通過(guò)多元回歸法,估算POC與SOC. Pachon等[30]在Blanchard基礎(chǔ)上增加K+指示生物質(zhì)燃燒排放的POC,與CO、EC共同作為總POC的指示物.

    1.2 WSOC法

    水溶性有機(jī)物(WSOC)是總有機(jī)碳的重要組分部分.大氣顆粒物中WSOC主要來(lái)自生物質(zhì)燃燒排放[31]和二次轉(zhuǎn)化[32].總 WSOC扣除生物質(zhì)燃燒貢獻(xiàn)即可估算二次轉(zhuǎn)化來(lái)源的WSOC;而在生物質(zhì)燃燒可忽略的地區(qū),WSOC被用來(lái)指示二次有機(jī)碳[33].

    生物質(zhì)燃燒產(chǎn)生的WSOC可通過(guò)測(cè)量該排放源的示蹤物左旋葡聚糖(levoglucosan)來(lái)推算

    [34].較 K+而言,左旋葡聚糖對(duì)生物質(zhì)燃燒的指示更為獨(dú)特,且在有機(jī)溶劑可提取的有機(jī)物種中豐度高,因而是生物質(zhì)燃燒源較為理想的示蹤物選擇.此方法的原理如式(3)、(4)所示∶

    式中∶WSOCM表示測(cè)量的總 WSOC濃度; WSOCBB與(WSOC/levoglucosan)BB分別代表生物質(zhì)燃燒源排放的WSOC濃度及WSOC與左旋葡聚糖的濃度比值,而 Clevoglucosan則為環(huán)境中左旋葡聚糖的濃度.

    WSOC法所需檢測(cè)的物種相對(duì)較少,在實(shí)際應(yīng)用中簡(jiǎn)單易行.利用在線水萃取裝置(Particle into Liquid Sampler, PILS)[33]或在線氣體-氣溶膠收集儀(Online Gas and Aerosol Collector, GAC)連接TOC分析儀,可獲取連續(xù)在線的WSOC濃度,從而提供高時(shí)間分辨率SOC估算的基礎(chǔ)數(shù)據(jù).

    WSOC法假設(shè)WSOC僅來(lái)自于二次轉(zhuǎn)化和生物質(zhì)燃燒,忽略了 WSOC的其他一次來(lái)源,如機(jī)動(dòng)車排放[35].方法中(WSOC/levoglucosan)BB需借助當(dāng)?shù)厣镔|(zhì)燃燒源的源譜獲取,該比值具有明顯的地區(qū)差異性,且隨燃燒生物質(zhì)的種類和燃燒條件等不同.例如,美國(guó)亞特蘭大的研究[36-37]顯示,計(jì)劃性生物質(zhì)燃燒產(chǎn)生的(WSOC/ levoglucosan)BB約為9.1±2.9,而野火燃燒排放產(chǎn)生的比值約為 4.5~6.8.此外,WSOC法忽略了部分水溶性較低的 SOC貢獻(xiàn),可能會(huì)導(dǎo)致對(duì) SOC的低估.

    1.3 受體模型法

    化學(xué)質(zhì)量平衡模型(Chemical Mass Balance, CMB)、正定矩陣分解模型(Positive Matrix Factorization, PMF)是目前大氣顆粒物來(lái)源研究中廣泛應(yīng)用的兩種受體模型.二者主要用于解析一次源對(duì)受體點(diǎn)顆粒物的濃度貢獻(xiàn),但近幾年也逐漸與其它手段結(jié)合同時(shí)解析一次與二次源對(duì)顆粒物或有機(jī)碳的濃度貢獻(xiàn).

    1.3.1 CMB-MM 法 化學(xué)質(zhì)量平衡模型結(jié)合分子示蹤物(CMB-MM)法是在CMB中利用有機(jī)分子作為源的示蹤物種定量各種源對(duì)顆粒物或OC貢獻(xiàn)的方法[4].在解析OC來(lái)源的研究中,解析結(jié)果中由一次源貢獻(xiàn)的OC之和視為POC,總OC與POC之差可間接作為SOC[38],估算如下∶

    其中,Sj代表一次源j對(duì)OC濃度的貢獻(xiàn),OCtot代表測(cè)量的總OC.

    CMB-MM模型法解析OC來(lái)源時(shí)輸入數(shù)據(jù)包括受體點(diǎn)有機(jī)示蹤物種的濃度及其不確定性,以及污染源成分譜;通常采用有效方差最小二乘法確定各污染源的貢獻(xiàn)值.污染源成分譜則盡可能應(yīng)用當(dāng)?shù)氐脑闯煞肿V,美國(guó)已建立了可共享的揮發(fā)性有機(jī)物和顆粒物的污染源成分譜數(shù)據(jù)庫(kù) SPECIATE,然而,我國(guó)源譜工作的開(kāi)展地區(qū)及研究的源譜種類均相對(duì)稀少,且主要集中在北方及東部沿海發(fā)達(dá)地區(qū),如北京、天津、濟(jì)南、南京等城市已針對(duì)燃煤源、機(jī)動(dòng)車源、塵源(土壤塵、道路塵)等源類建立了當(dāng)?shù)氐脑闯煞肿V,而冶煉源、生物質(zhì)燃燒、餐飲源等源類的研究較為有限[39].

    除一次污染源成分譜外,部分研究嘗試在CMB模型里引入虛擬的“純”O(jiān)C源成分譜,作為SOC源譜,解析獲得SOC和其他一次源類的貢獻(xiàn)[40].Shi等[41]基于 CMB模型結(jié)果,通過(guò)嵌套的迭代估算方法扣除SOC,并將扣除全部SOC后的受體成分譜與源成分譜納入 CMB模型,估算POC與SOC濃度,及各個(gè)源的貢獻(xiàn)值.CMB-MM法建立在所有一次污染源類均被鑒別的基礎(chǔ)上,若CMB模型的輸入源譜缺少某些POC貢獻(xiàn)源類,則會(huì)導(dǎo)致對(duì)SOC的高估.

    14C同位素示蹤技術(shù)可以區(qū)分現(xiàn)代生物源和化石燃料源對(duì)大氣含碳顆粒物的貢獻(xiàn)[42].主要原理基于大氣14CO2中的14C通過(guò)光合作用可儲(chǔ)存在新形成的含碳物質(zhì)中;而化石燃料的形成時(shí)間遠(yuǎn)遠(yuǎn)大于14C的半衰期(約5730年),化石燃料中不存在14C.通過(guò)測(cè)量14C,可定量估算顆粒物中來(lái)自化石燃料碳和現(xiàn)代碳的貢獻(xiàn).Ding等[42]利用14C技術(shù)解析顆粒物中總碳的現(xiàn)代源與化石燃料源的貢獻(xiàn)量,并將CMB-MM模型解析結(jié)果中木柴燃燒、餐飲、植物碎屑等源類統(tǒng)歸為一次現(xiàn)代源,柴油車、汽油車、焦化廠等源類劃分為一次化石燃料源;總碳與一次源碳解析結(jié)果中兩類源的貢獻(xiàn)量分別做差值,即對(duì)應(yīng)SOC中現(xiàn)代源和化石燃料源的貢獻(xiàn).

    1.3.2 AMS-PMF法 與CMB法不同,PMF法是基于統(tǒng)計(jì)分析的方法,需輸入大量顆粒物化學(xué)物種濃度數(shù)據(jù),以及各化學(xué)物種的測(cè)量偏差;但不需源成分譜的詳細(xì)信息,避免了源譜測(cè)試本身的不確定性,以及借用外地源譜而導(dǎo)致的估算結(jié)果偏差.利用PMF模型可解析出顆粒物的一次源類因子貢獻(xiàn),若將分配到具有二次來(lái)源特征因子(如硫酸鹽、硝酸鹽等)中的有機(jī)碳判定為 SOC,則可簡(jiǎn)單估算SOC的濃度.Yuan等[43]利用該法估算香港SOC平均濃度為4.25μgC/m3,約占香港年均總OC的46%.

    近年來(lái)氣溶膠質(zhì)譜技術(shù)在大氣環(huán)境領(lǐng)域得到廣泛應(yīng)用.AMS(Aerosol Mass Spectrometer, Aerodyne Inc., USA)可提供高時(shí)間分辨率的超細(xì)顆粒物PM1的化學(xué)成分信息,滿足PMF需大量數(shù)據(jù)的要求.通過(guò)篩選 AMS提供的高質(zhì)量分辨率有機(jī)氣溶膠離子碎片(m/z約1~300),輸入PMF受體模型進(jìn)行來(lái)源解析;由 AMS有機(jī)氣溶膠質(zhì)譜圖的時(shí)間序列,即單位時(shí)間分辨率內(nèi)大氣顆粒物平均有機(jī)氣溶膠質(zhì)譜圖,可分解出多個(gè)有機(jī)氣溶膠因子/源類質(zhì)譜圖、各因子/源類貢獻(xiàn)濃度,可區(qū)分顆粒物中代表一次排放的類似烴類有機(jī)氣溶膠(Hydrocarbon-like Organic Aerosol, HOA)和二次生成的含氧有機(jī)氣溶膠(Oxygenated Organic Aerosol, OOA)的貢獻(xiàn)[44-45].

    AMS質(zhì)譜圖中較為重要的有機(jī)碎片包括m/z 43、m/z 44、m/z 55和m/z 57,通過(guò)外場(chǎng)觀測(cè)和煙霧箱實(shí)驗(yàn)已證明這些碎片及其相對(duì)豐度可用于區(qū)分顆粒物源類.其中m/z 57 (C4H9+)通常來(lái)源于飽和直鏈烴類分子的斷裂,可代表一次源排放的顆粒有機(jī)物;m/z 44 (CO2+)一般形成于羧基在電子轟擊電離下的斷裂,可認(rèn)為來(lái)自于含氧有機(jī)物,而m/z 43 (C3H7+和C2H3O+)可來(lái)自于一次來(lái)源和二次/氧化來(lái)源的貢獻(xiàn)[46].

    黃曉峰等[9]利用 AMS-PMF法解析出深圳冬季大氣細(xì)粒子的 5類因子,將其中對(duì)一次有機(jī)氣溶膠特征離子碎片峰貢獻(xiàn)大的因子中有機(jī)物加和作為POA,而將分配到具有二次來(lái)源特征因子(如SO42-、NO3-等)中的有機(jī)物加和作為SOA,發(fā)現(xiàn)深圳冬季 SOA 濃度水平約為(9.41± 6.33)μg/m3,占總有機(jī)物質(zhì)量的40%,與POA相比空間上濃度變化小,體現(xiàn)了區(qū)域性二次污染物特征.Xiao等[47]利用AMS-PMF法估算表明PRD地區(qū)夏季SOA約占總有機(jī)物的53%~66%.

    此外,傅立葉變換紅外 FTIR技術(shù)可通過(guò)對(duì)有機(jī)官能團(tuán)進(jìn)行定性與定量來(lái)估算 SOA[48]. FTIR提供的官能團(tuán)信息與AMS有機(jī)離子碎片的PMF因子解析,對(duì)比表明兩種方法在SOA源類型的識(shí)別、各源/因子的質(zhì)量分?jǐn)?shù)和時(shí)間序列等方面都具有良好的一致性[49].

    1.4 示蹤物-產(chǎn)率法

    除上述針對(duì)SOA總量的估算方法外,近年來(lái)SOA來(lái)源研究成為控制二次有機(jī)污染的關(guān)鍵.異戊二烯和單萜烯是自然界中排放量較高的兩種天然源VOCs,揮發(fā)性芳香族化合物則是重要的人為源VOCs.近年來(lái),煙霧箱模擬實(shí)驗(yàn)中異戊二烯、單萜烯、倍半萜烯和芳香烴等VOCs的氧化產(chǎn)物在外場(chǎng)觀測(cè)中被檢測(cè)到[50-51].美國(guó)環(huán)保署(USEPA)建立了一種利用大氣顆粒物中可檢測(cè)的VOCs氧化產(chǎn)物作為二次有機(jī)示蹤物來(lái)估算人為源及天然源SOA的方法,即示蹤物-產(chǎn)率法[52].

    該方法基于煙霧箱實(shí)驗(yàn),首先通過(guò)GC/MS測(cè)定每種VOC前體物在煙霧箱中反應(yīng)生成的一系列有機(jī)示蹤物的濃度,將所有SOA有機(jī)示蹤物的濃度之和,除以煙霧箱反應(yīng)后產(chǎn)生的總氣溶膠(即總 SOA)的濃度,即為SOA的質(zhì)量分?jǐn)?shù)(FSOA),見(jiàn)式(7).

    利用煙霧箱獲得的OM/OC比值,經(jīng)過(guò)SOA與SOC的轉(zhuǎn)換,轉(zhuǎn)換為SOC質(zhì)量分?jǐn)?shù)(FSOC)[52],并假定煙霧箱實(shí)驗(yàn)中生成的二次有機(jī)示蹤物與對(duì)應(yīng)前體物產(chǎn)生SOA的質(zhì)量關(guān)系與實(shí)際大氣中形成的相同[52-53].從而,通過(guò)測(cè)定實(shí)際大氣中各種前體物產(chǎn)生的SOA示蹤物濃度借助FSOA(FSOC)可估算某一特定 VOC前體物產(chǎn)生的總 SOA (SOC)濃度,并推算不同來(lái)源 SOA(SOC)對(duì) OA (OC)的貢獻(xiàn).

    目前煙霧箱實(shí)驗(yàn)獲得的異戊二烯、α-蒎烯、β-丁香烯、甲苯等前體物對(duì)應(yīng)的 FSOC值[54]已被應(yīng)用于廣州[6]、北京[7]、香港[54]3個(gè)超大城市SOC的估算.研究表明,人為源VOCs前體物對(duì)北京和廣州等地(不論城、郊)SOC均有重要貢獻(xiàn);廣州人為源芳香烴前體物對(duì) SOC的貢獻(xiàn)高達(dá)79%,較北京(45%~54%)更為顯著,均遠(yuǎn)高于以天然源VOCs貢獻(xiàn)為主的香港地區(qū).

    國(guó)內(nèi)對(duì)SOA示蹤物-產(chǎn)率法的應(yīng)用較少,除了上述研究外,Fu等[55]測(cè)定了泰山上以α/β-蒎烯和 β-丁香烯為前體物的 SOA濃度,分別可占到總有機(jī)物的 10%和 9.8%.目前,該法主要考慮了以異戊二烯、單萜烯、倍半萜烯等天然源VOCs為前體物,以及甲苯等人為源芳香烴為前體物產(chǎn)生的 SOA.然而,僅僅考慮這幾種主要的前體物可能導(dǎo)致 SOA含量的低估.此外,該法忽略了云中過(guò)程形成的氧化性產(chǎn)物(如乙二醛等),這些產(chǎn)物分配到云滴、霧、濕氣溶膠中,經(jīng)過(guò)水相反應(yīng),轉(zhuǎn)化為羧酸類物質(zhì),具有較強(qiáng)的 SOA生成潛力

    [56].未來(lái)的研究應(yīng)側(cè)重于增加和完善對(duì)人為源前體物的二次有機(jī)示蹤物的定性與定量研究,確定更符合實(shí)際大氣環(huán)境的SOA質(zhì)量分?jǐn)?shù),減少示蹤物-產(chǎn)率法的不確定性.

    2 基于排放清單的數(shù)值模擬法

    空氣質(zhì)量模型(CMAQ、GEOS-CHEM、CAMx等)也逐漸被用于SOA的模擬,成為研究大尺度 SOA空間分布和源貢獻(xiàn)的一個(gè)重要手段和熱點(diǎn)方向[57].與受體模型估算 SOA的方法不同,空氣質(zhì)量模型利用煙霧箱實(shí)驗(yàn)獲得的SOA產(chǎn)率、氣-粒分配系數(shù)等信息建立 SOA模塊[25],在輸入準(zhǔn)確的源排放清單和氣象資料條件下,通過(guò)對(duì)影響SOA濃度水平的各種大氣過(guò)程(如VOCs前體物的排放和氧化過(guò)程、半揮發(fā)性有機(jī)物的氣-粒分配過(guò)程,以及傳輸、擴(kuò)散、沉降等)的模擬,提供SOA濃度的時(shí)空變化及其與排放源之間的關(guān)系.

    空氣質(zhì)量模式對(duì)SOA濃度的模擬值往往低于觀測(cè)值,其中一個(gè)重要原因是模式中對(duì)SOA反應(yīng)機(jī)制及前體物的了解尚有待完善.早期的CMAQ版本僅將單萜烯、長(zhǎng)鏈烷烴和芳香烴視為SOA來(lái)源[58],且未考慮在酸性和高NOx條件下 SOA的產(chǎn)率和生成速率的變化等[1];CMAQ 4.7版本考慮了芳香烴類VOCs前體物在氮氧化物存在下的SOA產(chǎn)率,增加了苯、異戊二烯、倍半萜烯等SOA前體物,以及乙二醛和丙酮醛的云中氧化反應(yīng)、顆粒態(tài)的低聚反應(yīng)、酸對(duì)異戊二烯生成SOA的催化作用等機(jī)制[59].GEOS-Chem模型可用于模擬由生物源萜烯、異戊二烯,以及人為活動(dòng)和生物質(zhì)燃燒排放的芳香烴類等產(chǎn)生的半揮發(fā)性含氧有機(jī)物通過(guò)平衡分配產(chǎn)生的SOA[60],同時(shí),還用于模擬乙二醛、丙酮醛的水相反應(yīng)生成的SOA[61];近年來(lái),GEOS-Chem模型中已增加了二羰基化合物、半揮發(fā)性和中間體VOCs生成SOA的機(jī)制[62-63].

    此外,減少SOA前體物源排放清單的不確定性[63],填補(bǔ)模式中有機(jī)氣溶膠缺失的重要源、細(xì)化源清單的空間、時(shí)間分辨率,都會(huì)使得模擬結(jié)果更能夠接近實(shí)測(cè)值.

    程艷麗等[64]利用二維空氣質(zhì)量模式,對(duì)珠三角區(qū)域尺度SOA進(jìn)行了模擬研究,表明SOA的生成具有明顯的光化學(xué)反應(yīng)特征,濃度高值出現(xiàn)在14∶00左右.Fu等[62]利用GEOS-Chem模型模擬中國(guó)大陸SOA濃度,表明中國(guó)地面21%的OC來(lái)自于二次源,而東部地區(qū)夏季高達(dá) 62%的 OC為二次來(lái)源;并指出中國(guó)最重要的SOC前體物是異戊二烯、單萜烯和芳香烴.Jiang等[65]利用WRF/Chem結(jié)合SOA模型(SORGAM)研究表明中國(guó)人為源對(duì)總SOA的貢獻(xiàn)約為35%,春、夏、秋、冬各季分別為41%、26%、39%、59%;模擬值對(duì)SOA濃度低估約為0~75%.

    3 各種估算方法的比較

    SOA的估算方法多樣,已從SOA總量的估算發(fā)展到針對(duì)各種特定前體物產(chǎn)生SOA濃度的細(xì)化估算.但現(xiàn)有方法各有利弊,仍缺乏普遍適用的“統(tǒng)一”方法(表1).

    總的來(lái)說(shuō),基于觀測(cè)數(shù)據(jù)的估算方法,如 EC示蹤法、WSOC法、AMS-PMF法等,簡(jiǎn)單易行,可結(jié)合高時(shí)間分辨率的在線數(shù)據(jù),快速估算環(huán)境大氣中 SOC(SOA)的總體濃度;但由于分析物種有限,未能建立 SOC(SOA)濃度與源之間的關(guān)系.相比之下,14C技術(shù)結(jié)合受體模型CMB-MM法可區(qū)分化石燃料源和現(xiàn)代生物源 SOA,示蹤物-產(chǎn)率法則可針對(duì)性地評(píng)估某一類特定 VOCs前體物對(duì) SOA的濃度貢獻(xiàn),而受限于煙霧箱實(shí)驗(yàn)對(duì)VOCs光化學(xué)產(chǎn)物的準(zhǔn)確認(rèn)識(shí),及SOA示蹤物的測(cè)量分析困難,示蹤物-產(chǎn)率法仍未得到廣泛的應(yīng)用.

    基于觀測(cè)數(shù)據(jù)的估算方法依賴于外場(chǎng)觀測(cè)的時(shí)間跨度和空間尺度,往往只能反映某一或某些觀測(cè)點(diǎn)位的大氣環(huán)境;而基于源清單的數(shù)值模擬估算可覆蓋較大的空間范圍,據(jù)某地區(qū)或某區(qū)域的源排放清單和氣象場(chǎng),可進(jìn)行大尺度范圍、高時(shí)間分辨率SOA濃度模擬及來(lái)源解析;空氣質(zhì)量模型模擬的準(zhǔn)確與否,與源清單、氣象場(chǎng)和SOA生成模塊完善程度緊密相關(guān).

    目前的研究多采用單一的估算方法,針對(duì)多種方法的對(duì)比研究相對(duì)較少.Pachon等[30]對(duì) EC示蹤法、回歸法、CMB模型法、PMF模型法等對(duì)SOA的估算結(jié)果進(jìn)行了較為綜合的比較,發(fā)現(xiàn)CMB模型評(píng)估結(jié)果最高(SOC占 OC的 46%), PMF模型評(píng)估結(jié)果最低(占OC的26%);而根據(jù)誤差傳遞、季節(jié)評(píng)估、日變化研究等對(duì)各種估算方法的不確定性和估算結(jié)果進(jìn)行比較發(fā)現(xiàn),回歸法具有最小的不確定性;在進(jìn)行長(zhǎng)時(shí)間序列SOC的估算時(shí),回歸法具有優(yōu)勢(shì).

    國(guó)內(nèi)針對(duì)二次有機(jī)氣溶膠的濃度水平、來(lái)源貢獻(xiàn)的研究尚處于起步階段.目前應(yīng)用最為廣泛的SOA估算方法是EC示蹤法,其他方法的使用相對(duì)較少.Xiao等[47]初步嘗試?yán)猛煌鈭?chǎng)觀測(cè)數(shù)據(jù)比較了EC示蹤法、WSOC法和AMS-PMF法的估算結(jié)果,發(fā)現(xiàn)三種方法具有較高的一致性.Yuan等[43]在香港的研究中,對(duì)比表明EC示蹤法所得夏、冬兩季SOC較PMF模型法分別高估約 70%~212%和 4%~43%,主要是由于缺乏隨時(shí)空變化的一次混合源OC/EC比值.Ding等[6]在珠三角的研究發(fā)現(xiàn),夏季示蹤物產(chǎn)率法與 EC-示蹤法的估算結(jié)果雖不一致但相關(guān)性很好,而在生物質(zhì)燃燒明顯的秋冬季節(jié),由于生物質(zhì)燃燒源中OC/EC比值高,EC示蹤法高估,明顯高于示蹤物-產(chǎn)率法.然而,目前國(guó)內(nèi)外針對(duì)同一觀測(cè)期間的多種方法(尤其是基于觀測(cè)和污染源排放清單方法)的比較,以及對(duì) OC質(zhì)量濃度閉合的研究均十分缺乏.

    表1 各種SOC估算方法比較Table 1 Comparison of different methods for estimating SOC levels

    4 結(jié)語(yǔ)

    在選擇 SOA估算方法時(shí),應(yīng)根據(jù)研究地區(qū)所具備的各估算方法輸入條件及所需達(dá)到的預(yù)期目標(biāo)來(lái)綜合考慮.對(duì)源清單較為完善的地區(qū),可結(jié)合空氣質(zhì)量模型或源清單 OC/EC比值法進(jìn)行估算;在源成分譜較完善的地區(qū),可利用CMB-MM受體模型法或WSOC法開(kāi)展;而在有機(jī)物種數(shù)據(jù)相對(duì)缺乏的地區(qū),則可選擇EC示蹤法進(jìn)行;在排放源以有機(jī)物貢獻(xiàn)為主的地區(qū),示蹤物-產(chǎn)率法、AMS-PMF法可提供源類別多樣化和高時(shí)間分辨率的 SOA信息.如果關(guān)注某區(qū)域大尺度內(nèi) SOA的濃度水平,則可采用數(shù)值模擬的手段來(lái)研究.未來(lái)的研究有待將多種 SOA的估算方法對(duì)比和綜合利用,并對(duì)各種方法的不確定性進(jìn)行評(píng)估.

    [1] Jang M, Czoschke N M, Lee S, et al. Heterogeneous atmospheric aerosol production by acid-catalyzed particle-phase reactions [J]. Science, 2002,298:814-817.

    [2] 謝紹東,于 淼,姜 明.有機(jī)氣溶膠的來(lái)源與形成研究現(xiàn)狀 [J].環(huán)境科學(xué)學(xué)報(bào), 2006,26(12):1933—1939.

    [3] Turpin B J, Huntzicker J J. Identification of secondary organic aerosol episodes and quantitation of primary and secondary organic aerosol concentrations during SCAQS [J]. Atmospheric Environment, 1995,29(23):3527-3544.

    [4] Zheng M, Cass G R, Ke L, et al. Source apportionment of daily fine particulate matter at Jefferson Street, Atlanta, GA, during summer and winter [J]. Journal of the Air and Waste Management Association, 2007,57(2):228-242.

    [5] Dan M, Zhuang G S, Li X X, et al. The characteristics of carbonaceous species and their sources in PM2.5in Beijing [J]. Atmospheric Environment, 2004,38(21):3443-3452.

    [6] Ding X, Wang X M, Gao B, et al. Tracer-based estimation of secondary organic carbon in the Pearl River Delta, south China [J]. Journal of Geophysical Research-Atmospheres, 2012,117, D05313,doi:10.1029/2011JD016596.

    [7] Guo S, Hu M, Guo Q F, et al. Primary sources and secondary formation of organic aerosols in Beijing, China [J]. Environmental Science and Technology, 2012,46(18):9846-9853.

    [8] Rodrigo S A, Raul G E M S, Manuel A L G. Estimations of primary and secondary organic carbon formation in PM2.5aerosols of Santiago City, Chile [J]. Atmospheric Environment, 2009,43(13):2125-2131.

    [9] 黃曉峰,趙倩彪,何凌燕,等.基于氣溶膠質(zhì)譜的二次有機(jī)氣溶膠識(shí)別 [J]. 中國(guó)科學(xué):化學(xué), 2010,40(10):1550-1557.

    [10] Cabada J C, Pandis S N, Subramanian R, et al. Estimating the secondary organic aerosol contribution to PM2.5using the EC tracer method [J]. Aerosol Science and Technology, 2004,38(51):140-155.

    [11] Lim H J, Turpin B J. Origins of primary and secondary organic aerosol in Atlanta: results of time-resolved measurements during the Atlanta supersite experiment [J]. Environmental Science and Technology, 2002,36(21):4489-4496.

    [12] 葉文媛,吳 琳,馮銀廠,等.大氣中二次有機(jī)氣溶膠估算方法研究進(jìn)展 [J]. 安全與環(huán)境學(xué)報(bào), 2011,11(1):127-130.

    [13] Gray H A, Cass G R, Huntzicker J J, et al. Characteristics of atmospheric organic and elemental carbon particle concentrations in Los Angeles [J]. Environmental Science and Technology, 1986,20(6):580-589.

    [14] Saylor R D, Edgerton E S, Hartsell B E. Linear regression techniques for use in the EC tracer method of secondary organic aerosol estimation [J]. Atmospheric Environment, 2006,40(39):7546-7556.

    [15] Na K, Sawant A A, Song C, et al. Primary and secondary carbonaceous species in the atmosphere of Western Riverside County, California [J]. Atmospheric Environment, 2004,38(9):1345-1355.

    [16] Lin P, Hu M, Deng Z, et al. Seasonal and diurnal variations of organic carbon in PM2.5in Beijing and the estimation of secondary organic carbon [J]. Journal of Geophysical Research-Atmospheres, 2009,114,D00G11,doi:10.1029/2008JD010902.

    [17] Turpin B J, Huntzicker J J, Larson S M, et al. Los Angeles summer midday particulate carbon: primary and secondary aerosol [J]. Environmental Science & Technology, 1991,25(10):1788-1793.

    [18] Cao J J, Lee S C, Chow J C, et al. Spatial and seasonal distributions of carbonaceous aerosols over China [J]. Journal of Geophysical Research-Atmospheres, 2007,112,D22S11,doi:10. 1029/2006JD008205.

    [19] Zhang X Y, Wang Y Q, Zhang X C, et al. Carbonaceous aerosol composition over various regions of China during 2006 [J]. Journal of Geophysical Research-Atmospheres, 2008,113, D14111,doi:10.1029/2007JD009525.

    [20] Gu J X, Bai Z P, Liu A X, et al. Characterization of atmospheric organic carbon and element carbon of PM2.5and PM10at Tianjin, China [J]. Aerosol and Air Quality Research, 2010,10:167-176.

    [21] Yang L X, Zhou X H, Wang Z, et al. Airborne fine particulate pollution in Jinan, China: Concentrations, chemical compositions and influence on visibility impairment [J]. Atmospheric Environment, 2012,55:506-514.

    [22] Zhang F W, Zhao J P, Chen J S, et al. Pollution characteristics of organic and elemental carbon in PM2.5in Xiamen, China [J]. Journal of Environmental Sciences, 2011,23(8):1342-1349.

    [23] Huang H, Ho K F, Lee S C, et al. Characteristics of carbonaceous aerosol in PM2.5: Pearl Delta River Region, China [J]. Atmospheric Research, 2012,104-105:227-236.

    [24] Zhang F W, Xu L L, Chen J S, et al. Chemical characteristics of PM2.5during haze episodes in the urban of Fuzhou, China [J]. Particuology, 2012,11(3):264-272.

    [25] Lee S, Wang Y H, Russell A G. Assessment of secondary organic carbon in the southeastern United States: A review [J]. Journal of the Air and Waste Management Association, 2010,60(11):1282-1292.

    [26] Blanchard C L, Hidy G M, Tanenbaum S, et al. Carbon in southeastern U.S. aerosol particles: Empirical estimates of secondary organic aerosol formation [J]. Atmospheric Environment, 2008,42(27):6710-6720.

    [27] Hu W W, Hu M, Deng Z Q, et al. The characteristics and origins of carbonaceous aerosol at a rural site of PRD in summer of 2006 [J]. Atmospheric Chemistry and Physics, 2012,12:1811-1822.

    [28] Yu S C, Dennis R L, Bhave P V, et al. Primary and secondary organic aerosols over the United States: estimates on the basis of observed organic carbon (OC) and elemental carbon (EC), and air quality modeled primary OC/EC ratios [J]. AtmosphereEnvironment, 2004,38(31):5257-5268.

    [29] Cheng Y, Zheng M, He K B, et al. Comparison of two thermal-optical methods for the determination of organic carbon and elemental carbon: Results from the southeastern United States [J]. Atmosphere Environment, 2011,45(11):1913-1918.

    [30] Pachon J E, Balachandran S, Hu Y T, et al. Comparison of SOC estimates and uncertainties from aerosol chemical composition and gas phase data in Atlanta [J]. Atmospheric Environment, 2010,44(32):3907-3914.

    [31] Mukai H, Ambe Y. Characterization of a humic acid-like brown substance in airborne particulate matter and tentative identification of its origin [J]. Atmospheric Environment, 1986, 20(5):813-819.

    [32] Saxena P, Hildemann L M. Water-Soluble organics in atmospheric particles: A critical review of the literature and application of thermodynamics to identify candidate compounds [J]. Journal of Atmospheric Chemistry, 1996,24(1):57-109.

    [33] Weber R J, Sullivan A P, Peltier R E, et al. A study of secondary organic aerosol formation in the anthropogenic-influenced southeastern United States [J]. Journal of Geophysical Research-Atmosphere, 2007,112,D13302,doi:10.1029/2007JD008408.

    [34] Fraser M P, and Lakshmanan K. Using levoglucosan as a molecular marker for the long-range transport of biomass combustion aerosols [J]. Environmental Science and Technology, 2000,34(21):4560-4564.

    [35] Fisseha R, Dommen J D, Gaeggeler K, et al. Online gas and aerosol measurement of water soluble carboxylic acids in Zurich [J]. Journal of Geophysical Research-Atmospheres, 2006,111, D12316,doi:10.1029/2005JD006782.

    [36] Yan B, Zheng M, Hu Y T, et al. Organic composition of carbonaceous aerosols in an aged prescribed fire plume [J]. Atmospheric Chemistry and Physics, 2008,8:6381-6394.

    [37] Yan B. Characterization and source apportionment of ambient PM2.5in Atlanta, Georgia: On-road emission, biomass burning and SOA impact [D]. United States: Georgia Institute of Technology, 2009.

    [38] Zheng M, Ke L, Edgerton E S, et al. Spatial distribution of carbonaceous aerosol in the southeastern United States using molecular markers and carbon isotope data [J]. Journal of Geophysical Research-Atmospheres, 2006,111,D10S06,doi:10.1029/2005JD006777.

    [39] 鄭 玫,張延君,閆才青,等.上海 PM2.5工業(yè)源譜的建立 [J]. 中國(guó)環(huán)境科學(xué), 2013,33(8):1354-1359.

    [40] Lee S, Liu W, Wang Y H, et al. Source apportionment of PM2.5:comparing PMF and CMB results for four ambient monitoring sites in the southeastern United States [J]. Atmospheric Environment, 2008,42(18):4126-4137.

    [41] Shi G L, Tian Y Z, Zhang Y F, et al. Estimation of the concentrations of primary and secondary organic carbon in ambient particulate matter: application of the CMB-iteration method. Atmospheric Environment, 2011,45:5692-5698.

    [42] Ding X, Zheng M, Edgerton E S, et al. Contemporary or fossil origin: Split of estimated secondary organic carbon in the southeastern United States [J]. Environmental Science and Technology, 2008,42(24):9122-9218.

    [43] Yuan Z B, Yu J Z, Lau A K H, et al. Application of positive matrix factorization in estimating aerosol secondary organic carbon in Hong Kong and its relationship with secondary sulfate [J]. Atmospheric Chemistry and Physics, 2006,6:25-34.

    [44] Huang X F, He L Y, Hu M, et al. Highly time-resolved chemical characterization of atmospheric submicron particles during 2008Beijing Olympic Games using an Aerodyne High-Resolution Aerosol Mass Spectrometer [J]. Atmospheric Chemistry and Physics, 2010,10(18):8933-8945.

    [45] Sun Y L, Zhang Q, Anastasio C, et al. Insights into secondary organic aerosol formed via aqueous-phase reactions of phenolic compounds based on high resolution mass spectrometry [J]. Atmospheric Chemistry and Physics, 2010,10(10):4809-4822.

    [46] Zhang Q, Worsnop D R, Canagaratna M R, et al. Hydrocarbonlike and oxygenated organic aerosols in Pittsburgh: insights into sources and processes of organic aerosols [J]. Atmosphere Chemistry and Physics, 2005,5(12):3289-3311,doi:10.5194/acp-5-3289-2005.

    [47] Xiao R, Takegawa N, Zheng M, et al. Characterization and source apportionment of submicron aerosol with aerosol mass spectrometer during the PRIDE-PRD 2006campaign [J]. Atmospheric Chemistry and Physics, 2011,11(14):6911-6929.

    [48] Maria S F, Russell L M, Turpin B J, et al. FTIR measurements of functional groups and organic mass in aerosol samples over the Caribbean [J]. Atmospheric Environment, 2002,36(33):5185-5196.

    [49] Liu S, Ahlm L, Day D A, et al. Secondary organic aerosol formation from fossil fuel sources contribute majority of summertime organic mass at Bakersfield [J]. Journal of Geophysical Research, 2012, 117,D00V26,doi:10.1029/2012JD018170.

    [50] Jaoui M, Kleindienst T E, Lewandowski M, et al. Identification and quantification of aerosol polar oxygenated compounds bearing carboxylic or hydroxyl groups. 2. Organic tracer compounds from monoterpenes[J]. Environmental Science and Technology, 2005,39(15):5661-5673.

    [51] Szmigielski R, Surratt J D, Gómez-González Y, et al. 3-methyl-1,2,3-butanetricarbocylic acid: An atmospheric tracer for terpene secondary organic aerosol [J]. Geophysical Research Letters, 2007,34(24),L24811,doi:10.1029/2007GL031338.

    [52] Kleindienst T E, Jaoui M, Lewandowski M, et al. Estimates of the contributions of biogenic and anthropogenic hydrocarbons to secondary organic aerosol at a southeastern US location [J]. Atmospheric Environment, 2007,41(37):8288-8300.

    [53] Kleindienst T E, Jaoui M, Lewandowski M, et al. The formation of SOA and chemical tracer compounds from the photoxidation of naphthalene and its methyl analogs in the presence and absence of nitrogen oxides [J]. Atmospheric Chemistry and Physics, 2012,12,8711-8726.

    [54] Hu D, Bian Q J, Li T W Y, et al. Contributions of isoprene, monoterpenes, β-caryophyllene, and toluene to secondary organic aerosols in Hong Kong during the summer of 2006 [J]. Journal of Geophysical Research, 2008,113,D22206,doi:10.1029/2008JD010437.

    [55] Fu P Q, Kawamura K, Kanaya Y, et al. Contributions of biogenic volatile organic compounds to the formation of secondary organic aerosols over Mt. Tai, Central East China [J]. Atmospheric Environment, 2010,44(38):4817-4826.

    [56] Lim Y B, Tan Y, Turpin B J. Chemical insights, explicit chemistry and yields of secondary organic aerosol from methylglyoxal and glyoxal [J]. Atmospheric Chemistry and Physics Discussion, 2013,13,4687-4725.

    [57] Lane T E, Donahue N M, Pandis S N. Simulating secondary organic aerosol formation using the volatility basis-set approach in a chemical transport model [J]. Atmospheric Environment, 2008,42(32):7439-7451.

    [58] Binkowski F S, Roselle S J. Model-3community multiscale air quality (CMAQ) model aerosol component 1. Model description [J]. Journal of Geophysical Research-Atmospheres, 2003,108,D6,4183,doi:10.1029/2001JD001409.

    [59] Carlton A G, Bhave P V, Napelenok S L, et al. Model representation of secondary organic aerosol in CAMQv4.7 [J]. Environmental Science & Technology, 2010,44(22):8553-8560.

    [60] Henze D K, Seinfeld J H, Ng N L, et al. Global modeling of secondary organic aerosol formation from aromatic hydrocarbons:high-vs. low-yield pathways [J]. Atmospheric Chemistry and Physics, 2008,8:2405—2421.

    [61] Fu T M, Jacob D J, Heald C L. Aqueous-phase reactive uptake of dicarbonyls as a source of organic aerosol over eastern North America [J]. Atmospheric Environment, 2009,43(10):1814—1822.

    [62] Fu T M, Cao J J, Zhang X Y, et al. Carbonaceous aerosols in China: top-down constraints on primary sources and estimation of secondary contribution [J]. Atmospheric Chemistry and Physics, 2012,12(5):2725-2746.

    [63] Pye H O T and Seinfeld J H. A global perspective on aerosol from low-volatility organic compounds [J]. Atmospheric Chemistry and Physics, 2010,10:4377—4401.

    [64] 程艷麗,李湉湉,白郁華,等.珠江三角洲區(qū)域大氣二次有機(jī)氣溶膠的數(shù)值模擬 [J]. 環(huán)境科學(xué), 2009,30(12):3441-3447.

    [65] Jiang F, Liu Q, Huang X X, et al. Regional modeling of secondary organic aerosol over China using WRF/Chem [J]. Journal of Aerosol Science, 2012,43(1):57-73.

    [66] Pun B K, Wu S Y, Seigneur C, et al. Uncertainties in modeling secondary organic aerosols: Three-dimensional modeling studies in Nashville/Western Tennessee [J]. Environmental Science and Technology, 2003,37(16):3647-3661.

    《中國(guó)環(huán)境科學(xué)》被Ei收錄

    根據(jù) Ei總部 2013年頒布的期刊收錄情況,《中國(guó)環(huán)境科學(xué)》已被 Ei數(shù)據(jù)庫(kù)作為源期刊收錄,詳見(jiàn)http://www.chinaeidata.com/periodical.htm

    《中國(guó)環(huán)境科學(xué)》編輯部

    2013-03-14

    A review of methods for quantifying secondary organic aerosol.


    ZHENG Mei, YAN Cai-qing, LI Xiao-ying, WANG Xue-song, ZHANG Yuan-hang*
    (State Key Joint Laboratory of Environmental Simulation and Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing 100871, China). China Environmental Science, 2014,34(3):555~564

    Regional air pollution is complex and becomes increasingly important in China. Among many others, secondary organic aerosol (SOA) is one of the most important components of PM2.5. This paper discusses various methods for quantifying SOA in the atmosphere (including methods based on the EC tracer, WSOC, receptor model, the SOA tracers, and air quality model), presents the basic principle of each method and points out that 1) the EC-tracer method, the WSOC method and the receptor model method are relatively simple and convenient, but limited by the availability of local source profiles and some specific tracers; 2) the SOA-tracer method is analytically challenging but can supply source-specific SOA information; and 3) the air quality model method can provide large scale spatial distribution of SOA. This paper also summarizes the most recent results of SOA research in China and abroad and indicates that SOA is important in organic aerosol, and anthropogenic VOCs play a significant role in SOA formation in China. The primary purpose of this review is to provide basic and integrated information and suggestion for future directions of SOA study in China.

    secondary organic aerosol;fine particulate matter;quantification method;source;China

    X513

    :A

    :1000-6923(2014)03-0555-10

    鄭 玫(1969-),女,貴州貴陽(yáng)人,教授,從事一次和二次來(lái)源的有機(jī)氣溶膠物化特性及大氣顆粒物來(lái)源解析的研究工作.發(fā)表論文65篇.

    2013-05-03

    國(guó)家自然科學(xué)基金資助項(xiàng)目(21190050,41175102);上海市環(huán)保科研項(xiàng)目(45210128-0432)

    * 責(zé)任作者, 教授, yhzhang@pku.edu.cn

    猜你喜歡
    氣溶膠貢獻(xiàn)生物質(zhì)
    中國(guó)共產(chǎn)黨百年偉大貢獻(xiàn)
    生物質(zhì)揮發(fā)分燃燒NO生成規(guī)律研究
    能源工程(2021年5期)2021-11-20 05:50:44
    《生物質(zhì)化學(xué)工程》第九屆編委會(huì)名單
    《造紙與生物質(zhì)材料》(英文)2020年第3期摘要
    為加快“三個(gè)努力建成”作出人大新貢獻(xiàn)
    氣溶膠傳播之謎
    氣溶膠中210Po測(cè)定的不確定度評(píng)定
    貢獻(xiàn)榜
    海洋貢獻(xiàn)2500億
    商周刊(2017年6期)2017-08-22 03:42:37
    四川盆地秋季氣溶膠與云的相關(guān)分析
    免费高清在线观看视频在线观看| 国产色婷婷99| or卡值多少钱| 日韩伦理黄色片| 国产成人免费观看mmmm| 亚洲真实伦在线观看| 日本与韩国留学比较| 深夜a级毛片| 蜜臀久久99精品久久宅男| 日本午夜av视频| 2022亚洲国产成人精品| 高清午夜精品一区二区三区| 男女那种视频在线观看| 欧美xxxx性猛交bbbb| 51国产日韩欧美| 国产大屁股一区二区在线视频| 一个人看视频在线观看www免费| 亚洲精品久久久久久婷婷小说| 一边亲一边摸免费视频| 久久精品夜色国产| 欧美日本视频| 免费av不卡在线播放| 真实男女啪啪啪动态图| a级毛色黄片| 如何舔出高潮| 尤物成人国产欧美一区二区三区| 亚洲国产精品国产精品| 欧美一级a爱片免费观看看| 成人综合一区亚洲| 成人综合一区亚洲| 欧美激情国产日韩精品一区| h日本视频在线播放| 成年女人看的毛片在线观看| 婷婷色麻豆天堂久久| 久久精品国产自在天天线| 免费av毛片视频| 夫妻午夜视频| 国产av国产精品国产| 国产精品熟女久久久久浪| 韩国高清视频一区二区三区| 久久精品国产鲁丝片午夜精品| 午夜视频国产福利| 国产白丝娇喘喷水9色精品| 纵有疾风起免费观看全集完整版 | 国产一区二区在线观看日韩| 亚洲国产欧美在线一区| 男女国产视频网站| 久久午夜福利片| 国产乱来视频区| 精品99又大又爽又粗少妇毛片| 国产乱人偷精品视频| 精品人妻一区二区三区麻豆| 联通29元200g的流量卡| 国产一区二区三区av在线| 色综合亚洲欧美另类图片| 麻豆国产97在线/欧美| 69人妻影院| 国语对白做爰xxxⅹ性视频网站| 亚洲内射少妇av| 亚洲国产精品专区欧美| 国国产精品蜜臀av免费| 女人十人毛片免费观看3o分钟| 欧美xxxx黑人xx丫x性爽| 成人无遮挡网站| 亚洲高清免费不卡视频| 久久久久精品性色| 欧美 日韩 精品 国产| 真实男女啪啪啪动态图| 日韩成人伦理影院| 天堂中文最新版在线下载 | 亚洲av不卡在线观看| 婷婷色av中文字幕| 建设人人有责人人尽责人人享有的 | 九九爱精品视频在线观看| 日日啪夜夜爽| 天堂av国产一区二区熟女人妻| 久久久欧美国产精品| h日本视频在线播放| 人体艺术视频欧美日本| 国产av在哪里看| 久久精品夜夜夜夜夜久久蜜豆| 国产黄片视频在线免费观看| 久久久久久久久久久丰满| 最近中文字幕2019免费版| 免费av毛片视频| 午夜福利视频1000在线观看| 午夜精品国产一区二区电影 | 亚洲国产精品专区欧美| 爱豆传媒免费全集在线观看| 欧美另类一区| 亚洲最大成人av| 精品久久久精品久久久| 国产69精品久久久久777片| 国产大屁股一区二区在线视频| 国产精品99久久久久久久久| 日本免费在线观看一区| 日本与韩国留学比较| 亚洲av免费在线观看| 少妇裸体淫交视频免费看高清| 六月丁香七月| 秋霞在线观看毛片| 亚洲av成人精品一区久久| 亚洲精品一二三| 九九在线视频观看精品| 亚洲欧美精品专区久久| 日韩大片免费观看网站| 亚洲精品久久午夜乱码| 国产片特级美女逼逼视频| 亚洲av电影不卡..在线观看| 黄色欧美视频在线观看| 久久久久国产网址| 日韩一本色道免费dvd| 国产不卡一卡二| 狂野欧美激情性xxxx在线观看| 18禁在线无遮挡免费观看视频| a级毛色黄片| 日韩视频在线欧美| 成人亚洲精品一区在线观看 | 一区二区三区免费毛片| 日本黄色片子视频| 又黄又爽又刺激的免费视频.| 精品少妇黑人巨大在线播放| 永久免费av网站大全| 草草在线视频免费看| 青青草视频在线视频观看| 天天躁日日操中文字幕| 国产精品国产三级国产专区5o| 真实男女啪啪啪动态图| 亚洲欧美成人精品一区二区| 天堂中文最新版在线下载 | 1000部很黄的大片| 久久久成人免费电影| or卡值多少钱| 国产在视频线精品| 国产一级毛片在线| 亚洲丝袜综合中文字幕| 欧美最新免费一区二区三区| 久久久午夜欧美精品| 99久久中文字幕三级久久日本| 观看美女的网站| 午夜爱爱视频在线播放| 中文在线观看免费www的网站| 91精品一卡2卡3卡4卡| 永久免费av网站大全| 日日撸夜夜添| 亚洲国产精品成人综合色| a级毛色黄片| 秋霞伦理黄片| 午夜福利在线观看免费完整高清在| 在线观看一区二区三区| 久久久精品欧美日韩精品| 成年人午夜在线观看视频 | 欧美日韩综合久久久久久| 免费观看a级毛片全部| 久久精品综合一区二区三区| 好男人在线观看高清免费视频| 精品久久久久久久久久久久久| 日韩三级伦理在线观看| 男女视频在线观看网站免费| 免费在线观看成人毛片| 国产一区有黄有色的免费视频 | 成人国产麻豆网| 欧美日韩视频高清一区二区三区二| 亚洲成人久久爱视频| 观看免费一级毛片| 大香蕉久久网| ponron亚洲| 欧美日韩精品成人综合77777| 九色成人免费人妻av| 国精品久久久久久国模美| 51国产日韩欧美| 大又大粗又爽又黄少妇毛片口| 99视频精品全部免费 在线| eeuss影院久久| 欧美一区二区亚洲| 狂野欧美白嫩少妇大欣赏| 午夜精品一区二区三区免费看| 三级男女做爰猛烈吃奶摸视频| 亚洲精品日韩av片在线观看| 亚洲av电影在线观看一区二区三区 | 日本猛色少妇xxxxx猛交久久| 日韩欧美精品免费久久| 精品99又大又爽又粗少妇毛片| 免费av观看视频| 成人毛片60女人毛片免费| 亚洲av.av天堂| 午夜亚洲福利在线播放| 亚洲精品456在线播放app| 国产又色又爽无遮挡免| 国产乱人视频| 欧美另类一区| 亚洲人与动物交配视频| 精品久久久噜噜| 国产有黄有色有爽视频| 日韩制服骚丝袜av| 国产高清有码在线观看视频| 国产中年淑女户外野战色| 一级毛片aaaaaa免费看小| 肉色欧美久久久久久久蜜桃 | 国产国拍精品亚洲av在线观看| 亚洲综合精品二区| 又爽又黄a免费视频| 国产成人精品福利久久| 亚洲欧美日韩无卡精品| 舔av片在线| 非洲黑人性xxxx精品又粗又长| 午夜免费男女啪啪视频观看| 男人舔奶头视频| 国国产精品蜜臀av免费| 国产精品美女特级片免费视频播放器| 欧美一级a爱片免费观看看| 国产午夜福利久久久久久| 我的女老师完整版在线观看| 久久国内精品自在自线图片| 高清毛片免费看| 亚洲,欧美,日韩| 日本黄大片高清| 免费av毛片视频| 2021少妇久久久久久久久久久| 亚洲国产精品国产精品| 最近中文字幕2019免费版| 国产一区亚洲一区在线观看| 狠狠精品人妻久久久久久综合| 国产精品综合久久久久久久免费| 国产精品久久视频播放| 久久久久久久久久成人| 成人特级av手机在线观看| 直男gayav资源| 亚洲激情五月婷婷啪啪| 亚洲国产精品成人综合色| 久久精品国产鲁丝片午夜精品| 国产成人精品福利久久| 中文天堂在线官网| 亚州av有码| 国模一区二区三区四区视频| 亚洲精品乱码久久久v下载方式| 日韩,欧美,国产一区二区三区| 中文欧美无线码| 99久久精品国产国产毛片| 日日摸夜夜添夜夜添av毛片| 天堂√8在线中文| 亚洲av不卡在线观看| 22中文网久久字幕| 精品久久久久久久末码| 黄色配什么色好看| 国产精品精品国产色婷婷| 亚洲人成网站高清观看| 日韩制服骚丝袜av| 欧美极品一区二区三区四区| 国产三级在线视频| 三级毛片av免费| 三级男女做爰猛烈吃奶摸视频| 深爱激情五月婷婷| 成人午夜精彩视频在线观看| 免费电影在线观看免费观看| 熟妇人妻不卡中文字幕| 黄片无遮挡物在线观看| 97精品久久久久久久久久精品| 三级男女做爰猛烈吃奶摸视频| 国产精品女同一区二区软件| 小蜜桃在线观看免费完整版高清| 国产av不卡久久| 精品久久久久久成人av| 黄片无遮挡物在线观看| 免费在线观看成人毛片| 成人毛片60女人毛片免费| 亚洲精品影视一区二区三区av| 91精品一卡2卡3卡4卡| 亚洲国产欧美在线一区| 最近最新中文字幕大全电影3| 美女大奶头视频| 亚洲欧美中文字幕日韩二区| 国产精品久久久久久久电影| 六月丁香七月| 成年版毛片免费区| 日本免费a在线| 成人一区二区视频在线观看| 爱豆传媒免费全集在线观看| 99热这里只有精品一区| 美女大奶头视频| 久久精品夜色国产| 国产永久视频网站| 日本爱情动作片www.在线观看| 午夜精品在线福利| 亚洲国产精品专区欧美| 色哟哟·www| 在线免费十八禁| 看免费成人av毛片| 热99在线观看视频| 一区二区三区乱码不卡18| 亚洲丝袜综合中文字幕| 亚洲精品亚洲一区二区| 日本免费a在线| 日日摸夜夜添夜夜爱| 亚洲av不卡在线观看| 高清视频免费观看一区二区 | 久久6这里有精品| 欧美 日韩 精品 国产| 亚洲av二区三区四区| 2018国产大陆天天弄谢| 国产日韩欧美在线精品| 日韩三级伦理在线观看| 国产色爽女视频免费观看| 精品国产三级普通话版| 久久久久久久大尺度免费视频| 婷婷色麻豆天堂久久| 国产淫片久久久久久久久| 最近中文字幕2019免费版| 69av精品久久久久久| 日韩成人伦理影院| 国产精品精品国产色婷婷| 美女被艹到高潮喷水动态| 亚洲av免费在线观看| 久久久久久久久久久丰满| av国产免费在线观看| 韩国av在线不卡| 国产免费一级a男人的天堂| 亚洲va在线va天堂va国产| 亚洲av二区三区四区| 69人妻影院| 人妻一区二区av| 淫秽高清视频在线观看| 亚洲人成网站在线观看播放| 精品一区在线观看国产| 免费看不卡的av| 欧美精品国产亚洲| 国产探花在线观看一区二区| 成年女人看的毛片在线观看| 男插女下体视频免费在线播放| 伦理电影大哥的女人| 久久99热这里只频精品6学生| 小蜜桃在线观看免费完整版高清| av一本久久久久| 国产淫语在线视频| 色哟哟·www| 欧美97在线视频| 日韩三级伦理在线观看| 天堂网av新在线| 爱豆传媒免费全集在线观看| 精品久久久噜噜| 日韩欧美三级三区| 国产免费又黄又爽又色| 日韩大片免费观看网站| 日韩欧美国产在线观看| 99久久中文字幕三级久久日本| 久久精品国产亚洲av天美| 在线观看一区二区三区| av线在线观看网站| 一级毛片 在线播放| 精品久久久久久久久久久久久| 国产精品av视频在线免费观看| 午夜精品国产一区二区电影 | 国产 亚洲一区二区三区 | 久久草成人影院| 免费高清在线观看视频在线观看| 能在线免费观看的黄片| 日韩精品青青久久久久久| av女优亚洲男人天堂| 欧美日韩一区二区视频在线观看视频在线 | 国产亚洲午夜精品一区二区久久 | 91久久精品电影网| 国产成人aa在线观看| 1000部很黄的大片| 一区二区三区免费毛片| 一级毛片黄色毛片免费观看视频| 2021少妇久久久久久久久久久| 日本免费在线观看一区| 成人漫画全彩无遮挡| 午夜免费男女啪啪视频观看| 激情五月婷婷亚洲| 看黄色毛片网站| 久久久久久久午夜电影| 中文资源天堂在线| 天美传媒精品一区二区| 国产精品.久久久| 丝袜美腿在线中文| 国产白丝娇喘喷水9色精品| 卡戴珊不雅视频在线播放| 女人久久www免费人成看片| 丝袜美腿在线中文| 欧美日韩综合久久久久久| 嫩草影院入口| 三级男女做爰猛烈吃奶摸视频| 毛片一级片免费看久久久久| 精品国产一区二区三区久久久樱花 | 高清午夜精品一区二区三区| 男人和女人高潮做爰伦理| 寂寞人妻少妇视频99o| 男人和女人高潮做爰伦理| 亚洲av日韩在线播放| 国产精品熟女久久久久浪| 两个人的视频大全免费| 乱系列少妇在线播放| 亚洲精品乱码久久久久久按摩| 国产黄色视频一区二区在线观看| av在线观看视频网站免费| 国产亚洲一区二区精品| eeuss影院久久| 蜜桃久久精品国产亚洲av| 亚洲怡红院男人天堂| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久久大av| 尤物成人国产欧美一区二区三区| 五月玫瑰六月丁香| 波多野结衣巨乳人妻| av一本久久久久| 久久亚洲国产成人精品v| 少妇裸体淫交视频免费看高清| 久久精品国产亚洲av涩爱| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩东京热| 亚洲国产精品成人综合色| av福利片在线观看| 精品欧美国产一区二区三| 纵有疾风起免费观看全集完整版 | 亚洲第一区二区三区不卡| 在线 av 中文字幕| 亚洲激情五月婷婷啪啪| 男人狂女人下面高潮的视频| 国产伦精品一区二区三区视频9| 极品少妇高潮喷水抽搐| av又黄又爽大尺度在线免费看| 午夜精品一区二区三区免费看| 精品久久国产蜜桃| 亚洲,欧美,日韩| 亚洲aⅴ乱码一区二区在线播放| 夫妻性生交免费视频一级片| 国产亚洲精品久久久com| 日韩欧美 国产精品| 免费黄色在线免费观看| 亚洲国产精品sss在线观看| 精品99又大又爽又粗少妇毛片| 亚洲国产精品专区欧美| 国产探花在线观看一区二区| 中国国产av一级| 亚洲欧美日韩东京热| 大片免费播放器 马上看| 夜夜看夜夜爽夜夜摸| 熟妇人妻不卡中文字幕| 日日啪夜夜撸| 极品教师在线视频| 99久久中文字幕三级久久日本| 夜夜爽夜夜爽视频| 搡老乐熟女国产| 国产精品99久久久久久久久| 午夜福利视频精品| 免费观看在线日韩| 午夜激情欧美在线| 3wmmmm亚洲av在线观看| 天堂影院成人在线观看| 亚洲av男天堂| 亚洲av成人精品一区久久| 久久久久久国产a免费观看| 中国美白少妇内射xxxbb| 搡老妇女老女人老熟妇| 国产精品一区二区在线观看99 | 成人高潮视频无遮挡免费网站| 卡戴珊不雅视频在线播放| 亚洲国产高清在线一区二区三| 午夜精品在线福利| 亚洲精品日韩在线中文字幕| 国产一级毛片在线| 免费电影在线观看免费观看| videos熟女内射| 亚洲av成人av| 亚洲,欧美,日韩| 精品久久久久久久末码| 中文天堂在线官网| 蜜桃亚洲精品一区二区三区| 国产成人91sexporn| 啦啦啦韩国在线观看视频| 18禁在线播放成人免费| 三级男女做爰猛烈吃奶摸视频| 国产片特级美女逼逼视频| 人人妻人人澡欧美一区二区| 亚洲人成网站在线播| 99久久中文字幕三级久久日本| 青春草视频在线免费观看| 国产视频首页在线观看| 免费无遮挡裸体视频| 亚洲精品国产av成人精品| 久久精品人妻少妇| 狠狠精品人妻久久久久久综合| 爱豆传媒免费全集在线观看| 中文资源天堂在线| 97热精品久久久久久| 一区二区三区高清视频在线| 蜜桃亚洲精品一区二区三区| 亚洲av免费在线观看| 日韩在线高清观看一区二区三区| 天天躁日日操中文字幕| 三级国产精品欧美在线观看| 国产亚洲一区二区精品| 国产精品久久久久久精品电影| 美女被艹到高潮喷水动态| 校园人妻丝袜中文字幕| 九九在线视频观看精品| 99热全是精品| 国产一区二区在线观看日韩| 亚洲三级黄色毛片| 国产黄片视频在线免费观看| 国产有黄有色有爽视频| 国产精品国产三级国产专区5o| 亚洲美女视频黄频| 欧美成人午夜免费资源| 我的女老师完整版在线观看| 亚洲欧美成人精品一区二区| 97热精品久久久久久| 中文字幕av成人在线电影| 成人二区视频| 美女黄网站色视频| 人妻一区二区av| 男女视频在线观看网站免费| 天堂中文最新版在线下载 | 国产v大片淫在线免费观看| 免费不卡的大黄色大毛片视频在线观看 | 精品久久久久久久久亚洲| 99热这里只有是精品50| 视频中文字幕在线观看| 秋霞伦理黄片| 久久6这里有精品| 免费观看在线日韩| a级毛色黄片| 九九久久精品国产亚洲av麻豆| 国产色爽女视频免费观看| 99九九线精品视频在线观看视频| 国产av码专区亚洲av| 国产精品国产三级国产专区5o| 亚洲国产精品专区欧美| 亚洲精华国产精华液的使用体验| 成人无遮挡网站| 日韩欧美 国产精品| 18禁动态无遮挡网站| 亚洲av免费在线观看| 久久久久久久久中文| av在线播放精品| 亚洲欧美精品自产自拍| 亚洲国产精品成人久久小说| 少妇丰满av| 日韩欧美一区视频在线观看 | 亚洲精品中文字幕在线视频 | 国产精品一区二区三区四区久久| 伦精品一区二区三区| 亚洲四区av| 黄片无遮挡物在线观看| 亚洲在线观看片| 少妇裸体淫交视频免费看高清| 亚洲电影在线观看av| 欧美3d第一页| 国产男人的电影天堂91| 国产一区二区三区综合在线观看 | 又粗又硬又长又爽又黄的视频| 美女被艹到高潮喷水动态| or卡值多少钱| 尾随美女入室| 搡老妇女老女人老熟妇| 看免费成人av毛片| 一夜夜www| www.av在线官网国产| 一个人观看的视频www高清免费观看| 免费观看av网站的网址| 国产人妻一区二区三区在| 亚洲欧美精品专区久久| 欧美高清性xxxxhd video| 国产精品人妻久久久影院| kizo精华| av卡一久久| 国产午夜福利久久久久久| 国产精品久久久久久精品电影小说 | 亚洲精品乱久久久久久| 成年女人在线观看亚洲视频 | 91久久精品国产一区二区三区| 久久久久免费精品人妻一区二区| 久久6这里有精品| 夫妻性生交免费视频一级片| 亚洲精品国产av成人精品| 纵有疾风起免费观看全集完整版 | 成人美女网站在线观看视频| 欧美精品国产亚洲| 国产成人a区在线观看| 国产亚洲午夜精品一区二区久久 | 亚洲av中文字字幕乱码综合| 一级毛片电影观看| 久久亚洲国产成人精品v| 免费大片18禁| 不卡视频在线观看欧美| 在线观看免费高清a一片| 成年版毛片免费区| 一区二区三区乱码不卡18| 天堂中文最新版在线下载 | 久久6这里有精品| 欧美日韩综合久久久久久| 乱系列少妇在线播放| 国产三级在线视频| 99久久人妻综合| 国产黄a三级三级三级人| 国产伦一二天堂av在线观看| 男女视频在线观看网站免费| 国产精品久久久久久久电影| 狂野欧美白嫩少妇大欣赏| 2022亚洲国产成人精品| 又爽又黄无遮挡网站| 日韩欧美精品v在线| 在线免费十八禁| 亚洲电影在线观看av| 亚洲欧美清纯卡通| av网站免费在线观看视频 | 中文字幕人妻熟人妻熟丝袜美| av又黄又爽大尺度在线免费看| 成年女人在线观看亚洲视频 | 一二三四中文在线观看免费高清| .国产精品久久| 亚洲四区av| 神马国产精品三级电影在线观看| 干丝袜人妻中文字幕|