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

    地殼元素鐵的吸光貢獻(xiàn)對(duì)黑碳吸光增強(qiáng)估算的影響——以武漢為例

    2022-07-19 01:10:56孔少飛成春雷陶麗萍吳柏禧吳良斌
    中國(guó)環(huán)境科學(xué) 2022年7期
    關(guān)鍵詞:光吸收吸光依賴(lài)性

    譚 健,夏 瑞,吳 兌,孔少飛,陳 楠,成春雷,鄧 濤,陶麗萍,張 雪,吳柏禧,吳良斌,王 慶,吳 晟**

    地殼元素鐵的吸光貢獻(xiàn)對(duì)黑碳吸光增強(qiáng)估算的影響——以武漢為例

    譚 健1,2,3,夏 瑞1,2,3,吳 兌1,2,3,孔少飛4,5,陳 楠5,6*,成春雷1,2,3,鄧 濤7,陶麗萍1,2,3,張 雪1,2,3,吳柏禧1,2,3,吳良斌1,2,3,王 慶1,2,3,吳 晟1,2,3**

    (1.暨南大學(xué),質(zhì)譜儀器與大氣環(huán)境研究所,廣東 廣州 510632;2.粵港澳環(huán)境質(zhì)量協(xié)同創(chuàng)新聯(lián)合實(shí)驗(yàn)室,廣東 廣州 510632;3.廣東省大氣污染在線(xiàn)源解析系統(tǒng)工程技術(shù)研究中心,廣東 廣州 510632;4.中國(guó)地質(zhì)大學(xué)(武漢)環(huán)境學(xué)院,湖北 武漢 430074;5.湖北省大氣復(fù)合污染研究中心,湖北 武漢 430074;6.湖北省生態(tài)環(huán)境監(jiān)測(cè)中心站,湖北 武漢 430074;7.中國(guó)氣象局廣州熱帶海洋氣象研究所,廣東 廣州 510640)

    傳統(tǒng)方法使用固定的波長(zhǎng)吸收指數(shù)來(lái)估算地殼元素或棕碳吸光,但該方法只能處理僅存在兩種組分的情形(黑碳/棕碳,或黑碳/地殼元素)的吸光貢獻(xiàn),存在較大不確定性.為此,本文提出一種新的方法,利用武漢在線(xiàn)觀(guān)測(cè)數(shù)據(jù)(2021年2、3、8、9月),采用鐵作為地殼元素的示蹤物,利用最小相關(guān)系數(shù)法(MRS)得出地殼元素的吸光貢獻(xiàn),在扣除地殼元素吸光之后,再得到黑碳吸光增強(qiáng)系數(shù)(abs).結(jié)果顯示觀(guān)測(cè)期間370nm波段地殼元素吸光貢獻(xiàn)均值為12.3%,月均值范圍5.7%~15.5%,且波長(zhǎng)吸收指數(shù)與鐵的濃度正相關(guān),表明地殼元素吸光貢獻(xiàn)不可忽視.地殼元素吸光貢獻(xiàn)呈現(xiàn)出顯著的季節(jié)特征,呈現(xiàn)出冬季低,春季高的特點(diǎn).地殼元素吸光的分離前后計(jì)算的abs存在一定的差異,受到了黑碳和鐵的相關(guān)性的影響.觀(guān)測(cè)期間扣除地殼元素吸光后,abs均值為1.43±0.53,在季節(jié)性上呈現(xiàn)春夏高,秋冬低的特性,春季較高abs值與春季黑碳較高的老化程度有關(guān).abs對(duì)無(wú)機(jī)鹽和有機(jī)物的含量存在正的依賴(lài)關(guān)系,證明了這些包裹物質(zhì)對(duì)吸光增強(qiáng)的影響.

    礦塵吸光貢獻(xiàn);黑碳吸光增強(qiáng);最小相關(guān)系數(shù)法(MRS方法);波長(zhǎng)吸收指數(shù)()

    大氣氣溶膠可以吸收和散射太陽(yáng)輻射,對(duì)地球的輻射收支平衡產(chǎn)生重要影響[1].黑碳?xì)馊苣z(BC)是大氣氣溶膠中的主要吸光成分,主要來(lái)源于化石燃料和生物質(zhì)的不完全燃燒,可以直接在中低空吸收從紫外到近紅外較寬波段的太陽(yáng)光,造成正的輻射強(qiáng)迫[2].模式模擬和觀(guān)測(cè)研究表明,邊界層內(nèi)BC吸收太陽(yáng)輻射,導(dǎo)致抵達(dá)地表熱通量減少,形成了逆溫,進(jìn)而抑制邊界層的發(fā)展,加重大氣灰霾污染[3].

    排放到大氣中的BC經(jīng)歷老化過(guò)程,與有機(jī)碳或者二次無(wú)機(jī)鹽等污染物混合[4-5],其粒徑大小、形態(tài)結(jié)構(gòu)、吸濕性和光學(xué)特性等發(fā)生改變[6-7].BC與其他污染物混合會(huì)增強(qiáng)BC的光吸收能力,光吸收增強(qiáng)能力的大小與其混合狀態(tài)和外包裹物的多少有關(guān)[8-9].由于大氣中BC的性質(zhì)和狀態(tài)存在時(shí)空差異性,當(dāng)前對(duì)于BC輻射強(qiáng)迫的估計(jì)仍存在較大的不確定性[2].黑碳的光吸收增強(qiáng)通常使用吸光增強(qiáng)系數(shù)(abs)進(jìn)行表征.目前世界各地已有很多的學(xué)者通過(guò)現(xiàn)場(chǎng)觀(guān)測(cè)對(duì)BC的光吸收增強(qiáng)系數(shù)(abs)進(jìn)行研究報(bào)道,例如美國(guó)加利福尼亞州[10]、英國(guó)[8]、日本[11]和捷克[12],但不同研究之間abs存在一定的差異.最近的研究表明,不同地區(qū)abs的差異可能與BC粒子的混合狀態(tài)多樣性有關(guān)[13-16].我國(guó)的abs研究主要集中在華北平原[17-19]、長(zhǎng)三角[20-22]和珠三角地區(qū)[9,23-24].而中國(guó)中部城市武漢關(guān)于黑碳光吸收增強(qiáng)特性方面的研究較為薄弱.

    除了BC,礦物粉塵(MD)中的地殼元素對(duì)光吸收的也有一定的貢獻(xiàn),尤其是在農(nóng)村地區(qū)和受沙塵氣溶膠影響的時(shí)段[25-26].目前對(duì)MD光吸收的研究比較少,主要原因之一是缺乏直接定量的儀器和分析方法.我國(guó)是世界上受沙塵氣溶膠影響嚴(yán)重的區(qū)域之一,沙源地主要分布在西北和北部地區(qū)[27].前人研究發(fā)現(xiàn),在臨近沙源地的張掖市一次沙塵暴事件中,沙塵對(duì)光的吸收在總吸光的占比從26.7%增加到71.6%[26].在珠三角地區(qū)也觀(guān)測(cè)到了源于北方沙塵氣溶膠的粗粒子輸送導(dǎo)致的嚴(yán)重空氣污染事件[28],表明沙塵氣溶膠的影響不僅僅局限于其源地和周邊,對(duì)下風(fēng)向區(qū)域的影響也不容忽視.同時(shí)本地排放的道路揚(yáng)塵和建筑揚(yáng)塵對(duì)于MD也有一定的貢獻(xiàn)[29].當(dāng)前城市地區(qū)MD的光吸收貢獻(xiàn)因缺少合適的觀(guān)測(cè)手段而被忽略.使用單顆粒氣溶膠質(zhì)譜儀(SPAMS)研究發(fā)現(xiàn)北京BC與MD內(nèi)混合的粒子占比不到3%,表明MD與BC主要以外混形式存在[30].在計(jì)算abs的時(shí)候應(yīng)先剔除MD的光吸收貢獻(xiàn),否則可能會(huì)造成abs的高估.

    氣溶膠光學(xué)特性是影響大氣輻射平衡的重要因素.黑碳和棕碳是當(dāng)前研究較多的吸光物質(zhì),關(guān)于地殼元素的吸光貢獻(xiàn)卻鮮有報(bào)道.含鐵物質(zhì)是MD中的主要吸光組分,為此,本研究以武漢為例,提出一種新的估算MD吸光貢獻(xiàn)的方法,采用Fe作為示蹤物,利用最小相關(guān)系數(shù)(MRS)法計(jì)算MD的光吸收貢獻(xiàn),剔除掉MD的光吸收貢獻(xiàn)后再用MRS得出abs,從而獲得更為準(zhǔn)確的abs值,并分析其不同時(shí)間尺度變化和影響因素.

    1 外場(chǎng)觀(guān)測(cè)數(shù)據(jù)來(lái)源

    1.1 觀(guān)測(cè)站點(diǎn)和觀(guān)測(cè)時(shí)段

    圖1 觀(guān)測(cè)站點(diǎn)位置及主要觀(guān)測(cè)結(jié)果

    湖北省武漢市位于長(zhǎng)江中游,是中國(guó)中部地區(qū)的重要的經(jīng)濟(jì)中心和交通樞紐之一,本研究中使用的觀(guān)測(cè)數(shù)據(jù)來(lái)自湖北省大氣超級(jí)站,站點(diǎn)位于中國(guó)地質(zhì)大學(xué)(武漢)未來(lái)科技城校區(qū)“科教樓六”的5樓樓頂(114.61°E,30.46°N),如圖1所示.站點(diǎn)位置高于地面20m以上,周?chē)鸁o(wú)建筑物遮擋,站點(diǎn)周?chē)歉咝缕髽I(yè)區(qū)和居民區(qū),交通流量大,可以代表華中典型城市區(qū)域污染特征.本研究選取的觀(guān)測(cè)時(shí)間段為2021年2、3、8、9共四個(gè)月份,分別代表冬、春、夏、秋四個(gè)季節(jié).

    1.2 采樣儀器

    使用半連續(xù)在線(xiàn)ECOC分析儀(Model RT-4, Sunset Laboratory Inc, Tigard, Oregon, USA)來(lái)獲取PM2.5中有機(jī)碳(OC)和元素碳(EC)的小時(shí)濃度,進(jìn)樣切割頭為PM2.5,流量為8L/min.使用的方法為熱光透射法(TOT),是目前用來(lái)分析ECOC最為常用的方法之一,詳細(xì)的工作原理可參考文獻(xiàn)[31].

    使用AE-31型黑碳儀(AE-31, Magee Scientific Company, Berkeley, CA, USA)測(cè)得在7個(gè)光譜波段(370、470、520、590、660、880和950nm)下的氣溶膠光吸收系數(shù).AE-31型黑碳儀測(cè)得的是時(shí)間周期內(nèi)收集在膜上的氣溶膠光衰減系數(shù),與實(shí)際大氣中的氣溶膠吸收系數(shù)有差異,且存在負(fù)載效應(yīng),因此通常需要進(jìn)行數(shù)據(jù)校正,儀器詳細(xì)工作原理和數(shù)據(jù)校正可參考本課題組在珠三角的觀(guān)測(cè)研究[32].本文使用基于Weingartner等[33]的算法開(kāi)發(fā)的黑碳儀數(shù)據(jù)校正工具包[34]對(duì)AE-31黑碳儀數(shù)據(jù)進(jìn)行校正.基于儀器對(duì)比得到的經(jīng)驗(yàn)方程[1],本研究使用的校正因子(Cref)為3.23.波長(zhǎng)吸收指數(shù)()是表征氣溶膠光學(xué)特性的重要參數(shù),通常用于表征不同波長(zhǎng)下吸光的依賴(lài)性,可用以下公式計(jì)算:

    式中:1表示不同的波長(zhǎng),abs_λ表示對(duì)波長(zhǎng)的光吸收系數(shù).因?yàn)樽靥急憩F(xiàn)出在短波長(zhǎng)更強(qiáng)的光吸收導(dǎo)致其要大于BC,基于此特性常用于棕碳和黑碳吸光的分離[35].另外,值還可用于分析污染物來(lái)源和BC混合狀態(tài).

    PM2.5中水溶性離子成分(NH4+、NO3-、SO42-、Cl-、Na+、K+、Mg2+、Ca2+)和無(wú)機(jī)元素分別使用在線(xiàn)離子色譜儀(MARGA, Metrohm, Switzerland)和X射線(xiàn)熒光光譜儀(EHH-X200, Skyray Instruments Inc, China)測(cè)得,時(shí)間分辨率為小時(shí),采樣流量為16.67L/min.PM2.5質(zhì)量用微震蕩天平(1405DF, Thermo Fisher Scientific Inc, USA)測(cè)得,常規(guī)氣體(NO, NO2, SO2, O3)使用賽默飛世爾科技公司的儀器測(cè)得.

    2 數(shù)據(jù)分析方法

    本研究使用最小相關(guān)系數(shù)法(MRS)進(jìn)行數(shù)據(jù)分析,MRS方法基于EC示蹤法[36],最初被用來(lái)估算二次VOC[37],吳晟[38]對(duì)MRS方法的可靠性做了系統(tǒng)的數(shù)值模型驗(yàn)證并開(kāi)發(fā)了基于Igor Pro的MRS軟件.MRS法近年來(lái)被廣泛的用于二次有機(jī)碳(SOC)[31],黑碳光吸收增強(qiáng)系數(shù)(abs)[34]和二次棕碳光吸收(abs_BrCsec)[35]的計(jì)算.其基本原理主要是基于一次示蹤物和二次特征值之間的最小相關(guān)關(guān)系將一次組分和二次組分進(jìn)行分離.本文提出一種新的計(jì)算的abs優(yōu)化方法,即先使用元素鐵作為示蹤物用MRS解析出MD的吸光系數(shù),然后剩余吸光用于abs的計(jì)算.下面介紹詳細(xì)的計(jì)算過(guò)程.

    2.1 用MRS方法計(jì)算地殼元素Fe的光吸收系數(shù)

    在以往的研究中,BC和BrC被認(rèn)為是造成城市大氣環(huán)境中光吸收的兩個(gè)主要?dú)馊苣z組分,往往忽略了MD的吸光貢獻(xiàn)[35].鐵的氧化物是MD粒子中的主要吸光物質(zhì),并且MD的吸光截面系數(shù)與沙塵粒子中鐵氧化物含量成正比[39],因此本研究中選取元素鐵作為示蹤物來(lái)表征MD的吸光貢獻(xiàn).將總的光吸收系數(shù)abs_λ_total分為兩個(gè)部分,一部分是MD吸光abs_λ_MD,剩下的部分是非MD吸光abs_λ_NMD:

    式中:代表不同的波段,abs_λ_MD可以用以下公式計(jì)算:

    式中:MAEFe_是含鐵物質(zhì)在波長(zhǎng)處的質(zhì)量吸收系數(shù),[Fe]代表儀器測(cè)得的元素鐵的質(zhì)量濃度,計(jì)算abs_λ_MD的關(guān)鍵在于得到一個(gè)合適的MAEFe_,本研究使用MRS法來(lái)獲取MAEFe_,這與(OC/EC)pri的計(jì)算原理類(lèi)似.首先,讓MAEFe_在一個(gè)合理的范圍內(nèi)連續(xù)取值,然后每一個(gè)假設(shè)的MAEFe__h值(后綴h表示hypothetical,下同)可以根據(jù)公式4計(jì)算出對(duì)應(yīng)的非MD吸光abs_λ_NMD_h:

    然后通過(guò)計(jì)算得出的abs_λ_NMD_h和[Fe]相關(guān)系數(shù)2(abs_λ_NMD_h, Fe),基于abs_λ_NMD_h和Fe之間的相互獨(dú)立性可以認(rèn)為當(dāng)2最小時(shí)得到最合理的MAEFe_,然后代入公式3中計(jì)算出相應(yīng)的abs_λ_MD,詳細(xì)的計(jì)算過(guò)程和相應(yīng)的MRS圖形示例,以及MRS的各種用法之間的對(duì)比總結(jié)在表1中.

    表1 MRS計(jì)算SOC、Eabs和babs_MD

    2.2 黑碳光吸收增強(qiáng)的計(jì)算

    分離出abs_λ_MD之后,剩余的吸光部分認(rèn)為是碳質(zhì)氣溶膠造成的吸光,包括一次吸光(abs_λ_pri)和老化后增強(qiáng)的吸光(abs_λ_aging).其中一次吸光為老化前的BC核的光吸收系數(shù),abs_λ_aging為經(jīng)歷大氣老化過(guò)程后其他污染物與BC混合并包裹在BC粒子外導(dǎo)致光吸收系數(shù)增大的部分,其具體計(jì)算過(guò)程可參考已有的研究[23],本研究中的計(jì)算公式和MRS圖形示例參見(jiàn)表1.黑碳的光吸收增強(qiáng)系數(shù)abs__NMD可以定義如下:

    式中:MAENMD_是特定波長(zhǎng)下的老化后黑碳的質(zhì)量吸收系數(shù),可由觀(guān)測(cè)得到.MAEp_是通過(guò)MRS方法計(jì)算得到的未老化的黑碳的質(zhì)量吸收系數(shù),BC老化后會(huì)導(dǎo)致abs__NMD>1.

    2.3 MRS方法估算SOC

    EC示蹤法被廣泛用于計(jì)算SOC,基于一次排放的有機(jī)碳(POC)與EC之間的相關(guān)關(guān)系將其從OC中分離,然后用OC減去POC得到SOC,具體公式如下:

    式中:[OC]和[EC]都是觀(guān)測(cè)得到的濃度值,而(OC/ EC)pri由MRS方法計(jì)算得到,MRS方法計(jì)算SOC的可靠性經(jīng)過(guò)了數(shù)值模型[38]以及外場(chǎng)觀(guān)測(cè)數(shù)據(jù)[40]驗(yàn)證.數(shù)值模型的研究[38]中考慮了諸多因素及三種排放源情形.在考慮了這些影響因素和排放源情形下,MRS給出的SOC均顯著優(yōu)于OC/EC最小比值法和OC/EC百分位數(shù)法.前人的外場(chǎng)觀(guān)測(cè)研究中[40], MRS給出的SOC與WSOC,及AMS給出OOA,都具有很高的相關(guān)性,從側(cè)面驗(yàn)證了MRS方法的可靠性,因此本文選取MRS法計(jì)算SOC.計(jì)算(OC/EC)pri的MRS圖形示例見(jiàn)表1.為了減少不同季節(jié)排放特征變化帶來(lái)的影響,本文使用每個(gè)月計(jì)算出相應(yīng)的(OC/EC)pri.

    MRS計(jì)算所用的軟件是本課題組開(kāi)發(fā)的基于Igor Pro的MRS工具包[38].本研究中散點(diǎn)圖、時(shí)間序列圖、相關(guān)性圖、箱式圖由本課題組開(kāi)發(fā)的基于Igor Pro的ScatterPlot工具包[41]和HistBox工具包繪制[34].

    3 結(jié)果與討論

    3.1 觀(guān)測(cè)期間污染物特征總體概述

    如圖1(a)所示,觀(guān)測(cè)期間PM2.5平均濃度為(34.26±19.26)μg/m3,整體呈對(duì)數(shù)正態(tài)分布,數(shù)據(jù)質(zhì)量良好.PM2.5濃度低于北京(80μg/m3)、上海(53μg/ m3)[42]和2013年在武漢的觀(guān)測(cè)(65μg/m3)[43],較低的PM2.5濃度表明近年來(lái)武漢市環(huán)境空氣質(zhì)量得到有效提升.本研究觀(guān)測(cè)期間平均光吸收系數(shù)abs520_total為(13.34±6.53)Mm-1,遠(yuǎn)小于2012~2014年的觀(guān)測(cè)值(119.06Mm-1),表明武漢市近年來(lái)一次污染物排放顯著減少,這與在武漢的長(zhǎng)期觀(guān)測(cè)結(jié)果一致[44].由圖1可見(jiàn),主要二次無(wú)機(jī)鹽離子(NO3-, SO42-和NH4+)占比超過(guò)PM2.5中主要污染物組分總和的一半,表現(xiàn)出武漢以二次污染物為主的環(huán)境空氣污染特征.為進(jìn)一步探究PM2.5中主要污染物組分隨PM2.5濃度變化趨勢(shì),繪制出PM2.5中主要污染物組分占比和PM2.5濃度的依賴(lài)關(guān)系圖如圖2(a)所示.從圖中可以看出隨著污染程度的加重,NO3-的貢獻(xiàn)顯著增大,而有機(jī)物(OM)的貢獻(xiàn)逐漸減少,證明武漢地區(qū)高污染天氣可能與二次硝酸鹽的形成有關(guān).近期多個(gè)研究指出在我國(guó)大氣污染事件中NO3-對(duì)PM2.5的貢獻(xiàn)不斷增加[45-46],強(qiáng)調(diào)了在我國(guó)的空氣質(zhì)量管理中進(jìn)一步控制其前體物NO排放的必要性.

    圖2 觀(guān)測(cè)期間PM2.5組分占比對(duì)PM2.5濃度(a)和Ox濃度(b)的依賴(lài)性關(guān)系,其中OM=1.5*OC

    大氣氧化能力通常用O(O=O3+NO2)表征[47-48].如圖2(b)所示,隨著O的增加,有機(jī)氣溶膠和硫酸鹽占比呈現(xiàn)上升趨勢(shì),而硝酸鹽占比顯著下降,表明大氣中光化學(xué)氧化性的增強(qiáng)有利于有機(jī)物和硫酸鹽的形成并抑制硝酸鹽的形成.

    如表2所示,觀(guān)測(cè)期間OC和EC的平均質(zhì)量濃度分別為(5.70±2.24)μgC/m3和(1.47±0.77)μgC/m3, OC和EC濃度總體呈現(xiàn)對(duì)數(shù)正態(tài)分布(圖1).相比于中國(guó)其他站點(diǎn),如北京(4.1μgC/m3)[49]、香港(1.88μgC/m3)[50]、上海(2.46μgC/m3)[51]和廣州(2.67μgC/m3)[31],武漢的EC濃度處于較低的水平,但是要高于巴黎(1.4μgC/m3)[52]和紐約(0.95μgC/m3)[53].而OC濃度要低于香港(6.77μgC/m3)[50],北京(14μgC/ m3)[49]和石家莊(22.8 μgC/m3)[54].

    質(zhì)量吸收系數(shù)(MAE)定義為單位質(zhì)量濃度的氣溶膠光吸收系數(shù),是用來(lái)表示氣溶膠光吸收能力的重要光學(xué)參數(shù),MAE值越大,說(shuō)明單位質(zhì)量的光吸收能力越強(qiáng).如表2所示,本研究中觀(guān)測(cè)期間平均MAE520_NMD為(8.90±3.38) m2/g,要高于文獻(xiàn)中提出的在550nm處一次排放的BC的MAE(7.5± 1.2m2/g)[55],表明在武漢地區(qū)BC粒子存在明顯的老化現(xiàn)象.

    表3中總結(jié)了近年來(lái)在中國(guó)各地區(qū)進(jìn)行的有關(guān)于光吸收增強(qiáng)的研究.從表中可以看出中國(guó)大部分地區(qū)黑碳老化程度較高,abs值較高.本研究中觀(guān)測(cè)期間的平均abs520_NMD為1.43±0.53.黑碳光吸收增強(qiáng)系數(shù)高值主要存在中國(guó)北方地區(qū),例如北京(1.24~4.0)[18,56-58]、濟(jì)南(2.07和1.9)[19,59]、固城(1.32)[58]和運(yùn)城(2.25)[17].其他研究站點(diǎn)包括南京(1.6)[20]、西安(1.8)[60]、壽縣(2.3)[61]、泰山(2.0)[19]和青島(1.64和1.96)[62]也存在較高的abs值.只有兩個(gè)研究報(bào)道了較低的abs值,分別是廣州(1.07)[24]和泰州(1.13)[63].總的來(lái)說(shuō),我國(guó)各地區(qū)總體的abs值比較高,因此有必要繼續(xù)擴(kuò)大開(kāi)展有關(guān)黑碳老化特征的相關(guān)研究.表3中值得注意的是,傳統(tǒng)的方法獲取abs值不能兼顧觀(guān)測(cè)時(shí)長(zhǎng)和觀(guān)測(cè)時(shí)間分辨率.熱溶蝕器(TD)方法觀(guān)測(cè)時(shí)間分辨率高(分鐘級(jí)),但是觀(guān)測(cè)時(shí)長(zhǎng)較為有限(通常為數(shù)周).AFD方法可以用于歷史采樣膜從而獲得跨度較長(zhǎng)的abs時(shí)間序列,但是時(shí)間分辨率較低(通常為24hr).本研究使用的MRS方法可以兼顧觀(guān)測(cè)時(shí)長(zhǎng)和觀(guān)測(cè)時(shí)間分辨率,具有較好的推廣應(yīng)用潛力.例如,最近捷克的學(xué)者利用MRS方法實(shí)現(xiàn)了2013~2017年間歐洲中部背景站點(diǎn)的abs長(zhǎng)期變化趨勢(shì)的研究[12].

    表2 PM2.5組分和光學(xué)參數(shù)統(tǒng)計(jì)

    表3 我國(guó)不同地區(qū)Eabs研究對(duì)比

    3.2 污染物特征量季節(jié)性變化和日變化分析

    表2總結(jié)了不同季節(jié)污染物濃度和光學(xué)參數(shù)的均值.PM2.5表現(xiàn)出明顯的季節(jié)變化趨勢(shì),冬季濃度最高為(50.42±24.01)μg/m3,而夏季濃度(23.18± 11.03μg/m3)不足冬季濃度的一半,說(shuō)明夏季大氣污染程度較輕,與之前在武漢的研究一致[64].各季節(jié)PM2.5濃度日變化如圖3(a)所示,春季和冬季存在比較明顯的日變化趨勢(shì),上午和夜間濃度較高,10:00左右濃度開(kāi)始降低,直到14:00~16:00點(diǎn)出現(xiàn)谷值,而夏季和秋季日變化幅度較小.

    如表2所示,EC濃度在冬季最高,為(2.21±0.95) μgC/m3,夏季最低,為(1.15±0.51)μgC/m3,春季和秋季濃度差別不大,分別為(1.26±0.56)μgC/m3和(1.31± 0.48)μgC/m3.EC的日變化圖3(c)所示,冬季的日變化最為明顯,在04:00左右達(dá)到峰值,隨后開(kāi)始降低直到17:00達(dá)到谷值,春季、夏季和秋季的日變化趨勢(shì)類(lèi)似,變化相對(duì)平緩,峰值出現(xiàn)在08:00和凌晨,谷值出現(xiàn)在下午.EC屬于一次污染物,主要來(lái)源于不完全燃燒,因此EC的日變化主要與機(jī)動(dòng)車(chē)排放和邊界層變化有關(guān).OC濃度在秋季最高,春季最低,分別為(7.72±1.98)μgC/m3和(3.87±1.32)μgC/m3.OC的來(lái)源包括一次排放和二次生成,夏季SOC最高,平均濃度為(4.35±1.56)μgC/m3,而SOC/OC比值也高達(dá)0.73± 0.12,說(shuō)明夏季SOC生成較為活躍,推測(cè)與夏季的高輻射強(qiáng)度和較強(qiáng)的光化學(xué)反應(yīng)活性有關(guān).從圖3(e)中SOC/OC的日變化可以看出在春季、夏季和秋季在午后達(dá)到峰值,說(shuō)明午后SOC生成最為活躍,而冬季的峰值推遲到了20:00.如圖3(b)所示,OC在四個(gè)季節(jié)的日變化趨勢(shì)一致,均在17:00達(dá)到谷值,隨后開(kāi)始上升.

    二次無(wú)機(jī)鹽組分中,硝酸鹽表現(xiàn)出春冬季高夏秋季低的季節(jié)性特征,冬季和春季硝酸鹽平均濃度分別為(12.45±7.26)μg/m3和(10.46±5.90)μg/m3,夏季和秋季平均濃度低至為(2.42±2.58)μg/m3和(2.31± 1.87)μg/m3,推測(cè)可能是由于夏秋季節(jié)環(huán)境溫度較高,導(dǎo)致硝酸鹽的揮發(fā),影響了硝酸鹽在氣-粒之間的分配.如圖3(f)所示,硝酸鹽的四個(gè)季節(jié)日變化趨勢(shì)類(lèi)似,均表現(xiàn)為單峰型,冬季峰值出現(xiàn)在02:00,隨后濃度開(kāi)始降低直到19:00,而春季、夏季和秋季的日變化較為平緩,峰值出現(xiàn)在09:00,谷值出現(xiàn)在18:00,硝酸鹽的日變化也證實(shí)隨著下午溫度的升高,顆粒態(tài)硝酸鹽揮發(fā)導(dǎo)致濃度下降.相較硝酸鹽而言,硫酸鹽濃度的季節(jié)性差異不那么顯著(表2).硫酸鹽的日變化如圖3(g)所示,在冬季,硫酸鹽濃度變化呈現(xiàn)日間低夜間高的特點(diǎn).分析發(fā)現(xiàn)硫酸鹽在日變化尺度上與NO2只有在冬季存在正相關(guān)關(guān)系(=0.72),而其他季節(jié)均為負(fù)相關(guān)關(guān)系(春季=-0.44,夏季=-0.79,秋季=-0.85),這意味著冬季硫酸鹽的夜間高濃度現(xiàn)象,存在著與其他季節(jié)不同的生成路徑,可能與夜間高濕高污染物濃度條件下NO2氧化SO2的生成路徑有關(guān)[65].而在夏季和秋季,硫酸鹽峰值出現(xiàn)在中午,推測(cè)夏秋的日間硝酸鹽生成主要由光化學(xué)反應(yīng)主導(dǎo).

    圖4 不同季節(jié)EC與Fe的相關(guān)性

    黑碳光吸收增強(qiáng)系數(shù)abs520_NMD春季最高(1.78±0.72),其次是夏季(1.41±0.29)和冬季(1.33± 0.36),秋季(1.10±0.24)最低,盡管春季EC濃度并不是最高的,但是黑碳的老化程度最高,同時(shí)較高的MAE520_NMD(10.99±4.55m2/g)也進(jìn)一步證實(shí)春季BC的老化,如表2所示.此外,為了研究MD吸光解析對(duì)黑碳吸光增強(qiáng)計(jì)算的影響,本文還計(jì)算了未去除MD吸光時(shí)的abs520并列于表2中.對(duì)比分析發(fā)現(xiàn)abs520_NMD與abs520雖然觀(guān)測(cè)期間均值相近,但是在季節(jié)上存在明顯的差異.如表2所示,春季abs520(1.88±0.76)高于abs520_NMD(1.78±0.72),而非春季abs520略低于abs520_NMD.通過(guò)圖4中EC與Fe的相關(guān)性分析發(fā)現(xiàn)春季EC與Fe沒(méi)有相關(guān)性(2=0.05),推測(cè)是由于春季沙塵主導(dǎo)了Fe的來(lái)源,這種情況下不剔除MD吸光得出的abs520的會(huì)將MD吸光誤判為黑碳老化吸光貢獻(xiàn),導(dǎo)致abs520高估,然而,非春季EC與Fe存在弱相關(guān)性(2范圍0.13~0.26),推測(cè)部分Fe和EC來(lái)源相近,比如道路揚(yáng)塵,剎車(chē)盤(pán)磨損,以及部分燃燒源直接排放的Fe,在這種情況下在MRS計(jì)算abs520的過(guò)程中,EC同源的這部分Fe被歸為黑碳的一次吸光,導(dǎo)致MAEp的高估,進(jìn)而引起abs520低估.

    abs520_NMD日變化如圖3(h)所示,其中春季日變化最為顯著,表現(xiàn)出早上低、下午和晚上較高的特點(diǎn),早上較低的abs520_NMD可能與早高峰時(shí)期一次排放的BC增多有關(guān),而下午和夜間較高的abs520_NMD與二次污染物的形成有關(guān).氣溶膠波長(zhǎng)吸收指數(shù)可以用來(lái)表征氣溶膠在不同波長(zhǎng)下的吸光能力,通常認(rèn)為黑碳的AAE值為1左右,但是根據(jù)不同地點(diǎn)氣溶膠來(lái)源組成和粒徑分布的不同,往往會(huì)存在一定的差距.相比于BC,BrC在紫外波段表現(xiàn)出更強(qiáng)的光吸收,因此BrC的值大于BC.如表2所示,本文中觀(guān)測(cè)期間總體的370/880_total為0.79±0.17,分季節(jié)來(lái)看,冬季最高為0.87±0.18,夏季最低,為0.71± 0.15.

    3.3 地殼元素Fe的吸光貢獻(xiàn)

    本研究觀(guān)測(cè)期間Fe的平均質(zhì)量濃度為(0.31± 0.35)μg/m3(表2).Fe的濃度表現(xiàn)出明顯的季節(jié)性差異,其中春季(0.44±0.39μg/m3)濃度最高,其次是秋季(0.32±0.46μg/m3)和夏季(0.25±0.25μg/m3),冬季(0.23±0.15μg/m3)濃度最低,推測(cè)春季較高的Fe濃度與春季頻繁發(fā)生的北方沙塵活動(dòng)有關(guān)[66].觀(guān)測(cè)期間abs520_MD均值為(1.37±1.42) Mm-1,如表2所示.不同季節(jié)的abs520_MD如表2和圖5a所示,春季的abs520_MD最高,為(1.95±1.71) Mm-1,其次是秋季(1.40±1.84Mm-1)和夏季(1.07±0.82Mm-1),冬季最低(1.03±0.67Mm-1).盡管本研究中abs520_MD數(shù)值比較小,但是武漢總的光吸收系數(shù)本身不高,所以MD的光吸收還是占有一定的貢獻(xiàn).在520nm波段(表2),春季MD的光吸收占比能達(dá)到總的光吸收的14%,而夏季和秋季的貢獻(xiàn)占比分別為11%和12%,冬季占比較低,為5%.如圖5a所示,紫外波段(370nm)的結(jié)果稍高一點(diǎn),春夏秋冬四季占比分別為15%,12%, 14%和6%.總體而言,武漢市MD對(duì)總的光吸收貢獻(xiàn)不能忽略,尤其是在春季.

    與棕碳類(lèi)似,MD的光吸收貢獻(xiàn)偏向于太陽(yáng)短波輻射,因此其AAE要大于BC[39].作為MD中的主要光吸收物質(zhì),Fe與AAE370/880_total的依賴(lài)關(guān)系圖如圖5(b)所示,在觀(guān)測(cè)期間,AAE370/880_total表現(xiàn)出對(duì)Fe的濃度的正依賴(lài)性關(guān)系.從季節(jié)上來(lái)看,在四個(gè)季節(jié)也都表現(xiàn)出類(lèi)似的關(guān)系,表明本研究中MD的存在影響了氣溶膠的光學(xué)特性.之前的研究中計(jì)算沙塵的光吸收貢獻(xiàn)是基于A(yíng)AE模型,認(rèn)為AAEBC=1,基于此計(jì)算出BC的光吸收并將剩余的吸收部分歸為沙塵的光吸收[67].然而,這種方法存在很大的局限性,首先實(shí)際大氣中的AAEBC因其來(lái)源、形態(tài)特征和粒徑分布等的不同往往并不等于1[68],此外,BrC的存在也能導(dǎo)致AAE的增大,并不能完全歸因于沙塵.相比之下,本文中使用的示蹤法并沒(méi)有預(yù)先假設(shè)AAEBC和AAEMD,不同波段下的MAEFe是基于實(shí)測(cè)的數(shù)據(jù)利用MRS方法計(jì)算得出.本研究中計(jì)算出的AAE370/880_MD為1.13,要小于在在北京的研究報(bào)道的沙塵氣溶膠(1.89)[69],但是要大于表2中觀(guān)測(cè)期間總的光吸收計(jì)算得出的AAE370/880_total(0.79±0.17)和剔除MD的吸光部分后的AAE370/880_NMD(0.77±0.18),從側(cè)面證明了元素鐵作為示蹤物用于MRS計(jì)算MD吸光貢獻(xiàn)的合理性.在遠(yuǎn)距離傳輸過(guò)程中,不同粒徑段的沙塵氣溶膠的干濕沉降的效果不同,導(dǎo)致氣溶膠粒徑分布的改變[70].有研究表明,沙塵氣溶膠在經(jīng)過(guò)長(zhǎng)距離傳輸過(guò)程中,不同粒徑段中的Fe的濃度變化存在很大的差距[71].由于值與粒徑分布密切相關(guān),推測(cè)在長(zhǎng)距離傳輸?shù)倪^(guò)程中粒徑分布的變化是導(dǎo)致本研究中AAE370/880_MD值低于文獻(xiàn)中北方城市AAE370/880_MD值的原因之一.

    3.4 影響Eabs520_NMD的因素

    BC粒子排放到大氣環(huán)境中會(huì)不斷老化,老化后的BC形態(tài)特征、光學(xué)特性和組分特征等都會(huì)發(fā)生變化,最為重要的是會(huì)導(dǎo)致其光吸收能力的增強(qiáng),導(dǎo)致abs大于1.研究表明,BC的老化程度與其污染源排放特征、混合狀態(tài)、形態(tài)特征和BC表面包裹層的組分特征有關(guān)[13-14,72].本文中探究了abs與EC、SOC/OC、NO3-/EC、SO42-/EC和非難熔組分與EC的質(zhì)量比值(BC)之間的關(guān)系.如圖6(a)所示,abs520_NMD與EC之間表現(xiàn)出負(fù)的依賴(lài)性關(guān)系,即EC質(zhì)量濃度較低時(shí)abs520_NMD值較高,而隨著EC濃度的增加,abs520_NMD值降低,這與我們課題組在廣州的光吸收增強(qiáng)研究中的結(jié)果類(lèi)似[23],推測(cè)可能是當(dāng)EC濃度較高時(shí)主要由本地一次排放主導(dǎo),黑碳老化程度低,所以abs520_NMD值較低,而EC濃度較低時(shí),長(zhǎng)距離輸送的EC占比相對(duì)增加,這些長(zhǎng)距離傳輸?shù)腅C顆粒物老化程度較高,導(dǎo)致整體abs520_NMD值抬升.不同季節(jié)下abs520_NMD與EC質(zhì)量濃度之間的關(guān)系如圖7(a)所示,abs520_NMD與EC之間的負(fù)依賴(lài)性關(guān)系主要體現(xiàn)在春季和冬季.

    SOC也是影響abs的重要因素之一,本研究用SOC/OC來(lái)表示SOC生成的活躍程度.如圖6(b)所示,abs520_NMD與SOC/OC之間表現(xiàn)出明顯的正依賴(lài)性關(guān)系.分季節(jié)來(lái)看,abs520_NMD與SOC/OC之間的依賴(lài)性關(guān)系在春季和冬季更為顯著,而秋季沒(méi)有表現(xiàn)出依賴(lài)性關(guān)系,說(shuō)明秋季SOC對(duì)abs520_NMD的貢獻(xiàn)較低.

    最近在南京的研究發(fā)現(xiàn)abs與NO3-/EC之間存在正的依賴(lài)性關(guān)系[22],實(shí)驗(yàn)室研究也表明BC與NO3-的混合能夠有效增強(qiáng)BC的光吸收能力[4].本研究中abs520_NMD與NO3-/EC的依賴(lài)性關(guān)系圖如圖6(c)所示,呈現(xiàn)出正的依賴(lài)性關(guān)系.雖然夏季和秋季NO3-濃度較低,但圖7(c)顯示NO3-與abs520_NMD之間還是存在一定的依賴(lài)性關(guān)系,表明NO3-對(duì)BC光吸收增強(qiáng)有顯著的貢獻(xiàn).本研究觀(guān)測(cè)期間abs520_NMD與SO42-/EC的依賴(lài)關(guān)系如圖6(d)和7(d)所示,abs520_NMD隨著SO42-/EC的增大而增大,并且在春季最為顯著,而秋季表現(xiàn)出較弱的正依賴(lài)性關(guān)系,夏季和冬季abs520_NMD隨著SO42-/EC增長(zhǎng)先增加隨后趨于穩(wěn)定.為了探究二次無(wú)機(jī)鹽離子與abs之間的關(guān)系,abs520_NMD對(duì)SIA/EC(SIA=NO3-+SO42-+NH4+)的依賴(lài)關(guān)系如圖6(e)和7(e)所示,除了夏季依賴(lài)性關(guān)系較弱外,其他三個(gè)季節(jié)均發(fā)現(xiàn)隨著SIA/EC的增加abs520_NMD也增加的現(xiàn)象,這些證據(jù)表明二次無(wú)機(jī)鹽的生成能有效增強(qiáng)BC的光吸收能力.

    圖6 觀(guān)測(cè)期間Eabs520_NMD對(duì)各種影響因素的依賴(lài)性關(guān)系

    圖7 不同季節(jié)Eabs520_NMD與其影響因素之間的依賴(lài)性關(guān)系

    BC通常定義為與BC粒子內(nèi)混合的非BC物質(zhì)(也稱(chēng)非難熔性組分)與BC的質(zhì)量之比,黑碳吸光增強(qiáng)系數(shù)abs與BC之間的曲線(xiàn)斜率因BC的老化程度、形態(tài)特征、混合狀態(tài)和粒徑分布的差異存在不同程度的依賴(lài)性關(guān)系[10,15,73-74].本研究使用無(wú)機(jī)鹽和有機(jī)物的總和與EC的比值得到BC,來(lái)研究abs520_NMD與BC之間的依賴(lài)性關(guān)系.需要注意的是,本研究中因?yàn)椴⒉灰欢ㄋ械腟IA和有機(jī)物都與BC以?xún)?nèi)混合的形式存在,所以BC可能被高估.盡管如此,觀(guān)測(cè)期間abs520_NMD與BC還是存在明顯的正相關(guān)依賴(lài)趨勢(shì),abs520_NMD隨著B(niǎo)C增大而增大,如圖6(f)所示.從季節(jié)上來(lái)看,不同季節(jié)之間abs520_NMD與BC的依賴(lài)性趨勢(shì)相近,但是曲線(xiàn)的斜率存在季節(jié)性差異,推測(cè)可能與不同季節(jié)BC老化后的不同的粒徑分布和混合態(tài)有關(guān)[16,73-74].此外,從圖6(c~f)可見(jiàn),春季的NO3-/EC,SO42-/EC,SIA/EC和BC的取值范圍大于其他季節(jié),這也是造成春季abs520_NMD高于其它季節(jié)的原因之一.

    4 結(jié)論

    4.1 觀(guān)測(cè)期間PM2.5平均濃度為(34.26±19.26) μg/m3, 與以往文獻(xiàn)報(bào)道的武漢觀(guān)測(cè)結(jié)果相比有著顯著下降,平均光吸收系數(shù)abs520_total為(13.34±6.53) Mm-1,遠(yuǎn)小于2012~2014年的觀(guān)測(cè)值(119.06Mm-1),表明武漢市近年來(lái)空氣質(zhì)量得到了顯著改善.

    4.2 當(dāng)PM2.5升高時(shí),硝酸鹽的占比顯著上升,而有機(jī)物的占比顯著下降.當(dāng)O升高時(shí)出現(xiàn)了相反的現(xiàn)象,硝酸鹽的占比隨著O升高顯著下降,而有機(jī)物的占比隨著O升高顯著上升.

    4.3 武漢市區(qū)的MD吸光不可忽略.觀(guān)測(cè)期間的abs520_MD均值為(1.37±1.42) Mm-1,并存在顯著季節(jié)性變化,春季的abs520_MD最高,為(1.95±1.71) Mm-1,其次是秋季(1.40±1.84Mm-1)和夏季(1.07± 0.82Mm-1),冬季最低(1.03±0.67Mm-1).在520nm波段,春季MD吸光在總吸光的占比能達(dá)到14%,而夏季和秋季的貢獻(xiàn)占比分別為11%和12%,冬季占比較低,為5%.紫外波段(370nm)的結(jié)果稍高一點(diǎn),春夏秋冬四季占比分別為16%,12%,14%和6%.此外,分析發(fā)現(xiàn)370/880_total對(duì)Fe的濃度存在正依賴(lài)關(guān)系,這一系列結(jié)果證明了MD的吸光貢獻(xiàn)即使在遠(yuǎn)離沙源地的城市區(qū)域也是不可忽略的.

    4.4 本研究中觀(guān)測(cè)期間平均MAE520_NMD為(8.90± 3.38) m2/g.在剔除了MD吸光貢獻(xiàn)之后,得到的黑碳光吸收增強(qiáng)系數(shù)abs520_NMD為1.43±0.53, 這與近期我國(guó)其他城市所報(bào)道的abs值相近.MD吸光的分離對(duì)abs計(jì)算的影響與EC和Fe的相關(guān)性有關(guān),非春季EC和Fe存在的弱相關(guān)性會(huì)導(dǎo)致abs低估,而春季EC與Fe不存在相關(guān)性,如不剔除MD吸光,會(huì)導(dǎo)致abs高估.abs520_NMD存在著顯著的季節(jié)性變化,春季最高(1.78±0.72),其次是夏季(1.41±0.29)和冬季(1.33±0.36)和秋季(1.10±0.24).abs520_NMD對(duì)無(wú)機(jī)鹽和有機(jī)物的含量存在正的依賴(lài)關(guān)系,證明了這些包裹物質(zhì)對(duì)吸光增強(qiáng)的影響.

    [1] Wu D, Mao J T, Deng X J, et al. Black carbon aerosols and their radiative properties in the Pearl River Delta region [J]. Sci. China Ser. D., 2009,52(8):1152-1163.

    [2] Bond T C, Doherty S J, Fahey D W, et al. Bounding the role of black carbon in the climate system: A scientific assessment [J]. J. Geophys. Res., 2013,118(11):5380-5552.

    [3] Ding A J, Huang X, Nie W, et al. Enhanced haze pollution by black carbon in megacities in China [J]. Geophys. Res. Lett., 2016,43(6): 2873-2879.

    [4] Yuan C, Zheng J, Ma Y, et al. Significant restructuring and light absorption enhancement of black carbon particles by ammonium nitrate coating [J]. Environ. Pollut., 2020,262:114172.

    [5] Schnaiter M, Linke C, Mohler O, et al. Absorption amplification of black carbon internally mixed with secondary organic aerosol [J]. J. Geophys. Res., 2005,110(D19):D19204.

    [6] Zhang R Y, Khalizov A F, Pagels J, et al. Variability in morphology, hygroscopicity, and optical properties of soot aerosols during atmospheric processing [J]. Proc. Natl. Acad. Sci. USA, 2008,105(30): 10291-10296.

    [7] Pei X, Hallquist M, Eriksson A C, et al. Morphological transformation of soot: investigation of microphysical processes during the condensation of sulfuric acid and limonene ozonolysis product vapors [J]. Atmos. Chem. Phys., 2018,18(13):9845-9860.

    [8] Liu S, Aiken A C, Gorkowski K, et al. Enhanced light absorption by mixed source black and brown carbon particles in UK winter [J]. Nat. Commun., 2015,6(1):8435.

    [9] Sun J Y, Wu C, Wu D, et al. Amplification of black carbon light absorption induced by atmospheric aging: temporal variation at seasonal and diel scales in urban Guangzhou [J]. Atmos. Chem. Phys., 2020,20(4):2445-2470.

    [10] Cappa C D, Onasch T B, Massoli P, et al. Radiative absorption enhancements due to the mixing state of atmospheric black carbon [J]. Science, 2012,337(6098):1078-1081.

    [11] Ueda S, Nakayama T, Taketani F, et al. Light absorption and morphological properties of soot-containing aerosols observed at an East Asian outflow site, Noto Peninsula, Japan [J]. Atmos. Chem. Phys., 2016,16(4):2525-2541.

    [12] Mbengue S, Zikova N, Schwarz J, et al. Mass absorption cross-section and absorption enhancement from long term black and elemental carbon measurements: A rural background station in Central Europe [J]. Sci. Total Environ., 2021,794:148365.

    [13] Fierce L, Bond T C, Bauer S E, et al. Black carbon absorption at the global scale is affected by particle-scale diversity in composition [J]. Nat. Commun., 2016,7:12361.

    [14] Matsui H, Hamilton D S, Mahowald N M. Black carbon radiative effects highly sensitive to emitted particle size when resolving mixing-state diversity [J]. Nat. Commun., 2018,9(1):3446.

    [15] Cappa C D, Zhang X, Russell L M, et al. Light absorption by ambient black and brown carbon and its dependence on black carbon coating state for two California, USA cities in winter and summer [J]. J. Geophys. Res., 2019,124(3):1550-1577.

    [16] Wu Y, Cheng T, Zheng L. Light absorption of black carbon aerosols strongly influenced by particle morphology distribution [J]. Environ. Res. Lett., 2020,15(9):094051.

    [17] Cui X, Wang X, Yang L, et al. Radiative absorption enhancement from coatings on black carbon aerosols [J]. Sci. Total Environ., 2016,551:51-56.

    [18] Zhang Y, Zhang Q, Cheng Y, et al. Amplification of light absorption of black carbon associated with air pollution [J]. Atmos. Chem. Phys., 2018,18(13):9879-9896.

    [19] Bai Z, Cui X, Wang X, et al. Light absorption of black carbon is doubled at Mt. Tai and typical urban area in North China [J]. Sci. Total Environ., 2018,635:1144-1151.

    [20] Cui F, Chen M, Ma Y, et al. An intensive study on aerosol optical properties and affecting factors in Nanjing, China [J]. Journal of Environmental Sciences, 2016,40:35-43.

    [21] Chen D, Zhao Y, Lyu R, et al. Seasonal and spatial variations of optical properties of light absorbing carbon and its influencing factors in a typical polluted city in Yangtze River Delta, China [J]. Atmos. Environ., 2019,199:45-54.

    [22] Ma Y, Huang C, Jabbour H, et al. Mixing state and light absorption enhancement of black carbon aerosols in summertime Nanjing, China [J]. Atmos. Environ., 2020,222:117141.

    [23] 孫嘉胤,吳 晟,吳 兌,等.廣州城區(qū)黑碳?xì)馊苣z吸光增強(qiáng)特性研究 [J]. 中國(guó)環(huán)境科學(xué), 2020,40(10):4177-4189.

    Sun J Y, Wu C, Wu D, et al. The light absorbtion enhancement characteristics of black carbon aerosols in urban Guangzhou [J]. China Environmental Science, 2020,40(10):4177-4189.

    [24] Fu Y, Peng X, Guo Z, et al. Filter-based absorption enhancement measurement for internally mixed black carbon particles over southern China [J]. Sci. Total Environ., 2021,762:144194.

    [25] Mian C, Diehl T, Dubovik O, et al. Light absorption by pollution, dust, and biomass burning aerosols: a global model study and evaluation with AERONET measurements [J]. Ann. Geophys., 2009,27(9):3439-3464.

    [26] Wu X, Liu J, Wu Y, et al. Aerosol optical absorption coefficients at a rural site in Northwest China: The great contribution of dust particles [J]. Atmos. Environ., 2018,189:145-152.

    [27] Sun J, Zhang M, Liu T. Spatial and temporal characteristics of dust storms in China and its surrounding regions, 1960–1999: Relations to source area and climate [J]. J. Geophys. Res., 2001,106(D10):10325- 10333.

    [28] 吳 兌,吳 晟,李 菲,等.粗粒子氣溶膠遠(yuǎn)距離輸送造成華南嚴(yán)重空氣污染的分析 [J]. 中國(guó)環(huán)境科學(xué), 2011,31(4):540-545.

    Wu D, Wu C, Li F, et al. Air pollution episode in southern China due to the long range transport of coarse particle aerosol [J]. China Environmental Science, 2011,31(4):540-545.

    [29] Casotti Rienda I, Alves C A. Road dust resuspension: A review [J]. Atmos. Res., 2021,261:105740.

    [30] Xie C, He Y, Lei L, et al. Contrasting mixing state of black carbon- containing particles in summer and winter in Beijing [J]. Environ. Pollut., 2020,263:114455.

    [31] Wu C, Wu D, Yu J Z. Estimation and uncertainty analysis of secondary organic carbon using one-year of hourly organic and elemental carbon data [J]. J. Geophys. Res., 2019,124(5):2774-2795.

    [32] 程 丁,吳 晟,吳 兌,等.深圳市城區(qū)和郊區(qū)黑碳?xì)馊苣z對(duì)比研究 [J]. 中國(guó)環(huán)境科學(xué), 2018,38(5):1653-1662.

    Cheng D, Wu C, Wu D, et al. Comparative study on the characteristics of black carbon aerosol in urban and suburban areas of Shenzhen [J]. China Environmental Science, 2018,38(5):1653-1662.

    [33] Weingartner E, Saathoff H, Schnaiter M, et al. Absorption of light by soot particles: determination of the absorption coefficient by means of aethalometers [J]. J. Aerosol Sci., 2003,34(10):1445-1463.

    [34] Wu C, Wu D, Yu J Z. Quantifying black carbon light absorption enhancement with a novel statistical approach [J]. Atmos. Chem. Phys., 2018,18(1):289-309.

    [35] Wang Q, Han Y, Ye J, et al. High Contribution of secondary brown carbon to aerosol light absorption in the southeastern margin of Tibetan Plateau [J]. Geophys. Res. Lett., 2019,46(9):4962-4970.

    [36] Turpin B J, Huntzicker J J. Secondary formation of organic aerosol in the Los-Angeles Basin - a descriptive analysis of organic and elemental carbon concentrations [J]. Atmos. Environ., 1991,25(2):207-215.

    [37] Millet D B, Donahue N M, Pandis S N, et al. Atmospheric volatile organic compound measurements during the Pittsburgh Air Quality Study: Results, interpretation, and quantification of primary and secondary contributions [J]. J. Geophys. Res., 2005,110(D7):D07S.

    [38] Wu C, Yu J Z. Determination of primary combustion source organic carbon-to-elemental carbon (OC/EC) ratio using ambient OC and EC measurements: secondary OC-EC correlation minimization method [J]. Atmos. Chem. Phys., 2016,16(8):5453-5465.

    [39] Alfaro S C, Lafon S, Rajot J L, et al. Iron oxides and light absorption by pure desert dust: An experimental study [J]. J. Geophys. Res., 2004,109(D8):D08208.

    [40] 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]. Atmos. Chem. Phys., 2012,12(4):1811-1822.

    [41] Wu C, Yu J Z. Evaluation of linear regression techniques for atmospheric applications: the importance of appropriate weighting [J]. Atmos. Meas. Tech., 2018,11(2):1233-1250.

    [42] Zhang Y L, Cao F. Fine particulate matter (PM2.5) in China at a city level [J]. Sci. Rep., 2015,5(1):1-12.

    [43] Huang T, Chen J, Zhao W, et al. Seasonal variations and correlation analysis of water-soluble inorganic ions in PM2.5in Wuhan, 2013 [J]. Atmosphere, 2016,7(4):49.

    [44] Zheng H, Kong S, Zheng M, et al. A 5.5-year observations of black carbon aerosol at a megacity in Central China: Levels, sources, and variation trends [J]. Atmos. Environ., 2020,232:117581.

    [45] Xu Q, Wang S, Jiang J, et al. Nitrate dominates the chemical composition of PM2.5during haze event in Beijing, China [J]. Sci Total Environ, 2019,689:1293-1303.

    [46] Zhang Y, Vu T V, Sun J, et al. Significant changes in chemistry of fine particles in wintertime Beijing from 2007 to 2017: Impact of Clean Air Actions [J]. Environ. Sci. Technol., 2020,54(3):1344-1352.

    [47] Liu S C. Possible effects on tropospheric O3and OH due to No emissions [J]. Geophys. Res. Lett., 1977,4(8):325-328.

    [48] Canonaco F, Slowik J G, Baltensperger U, et al. Seasonal differences in oxygenated organic aerosol composition: implications for emissions sources and factor analysis [J]. Atmos. Chem. Phys., 2015,15(12): 6993-7002.

    [49] Ji D, Zhang J, He J, et al. Characteristics of atmospheric organic and elemental carbon aerosols in urban Beijing, China [J]. Atmos. Environ., 2016,125:293-306.

    [50] Wu C, Huang X H H, Ng W M, et al. Inter-comparison of NIOSH and IMPROVE protocols for OC and EC determination: implications for inter-protocol data conversion [J]. Atmos. Meas. Tech., 2016,9(9): 4547-4560.

    [51] Yao L, Huo J, Wang D, et al. Online measurement of carbonaceous aerosols in suburban Shanghai during winter over a three-year period: Temporal variations, meteorological effects, and sources [J]. Atmos. Environ., 2020,226:117408.

    [52] Bressi M, Sciare J, Ghersi V, et al. A one-year comprehensive chemical characterisation of fine aerosol (PM2.5) at urban, suburban and rural background sites in the region of Paris (France) [J]. Atmos. Chem. Phys., 2013,13(15):7825-7844.

    [53] Blanchard C L, Shaw S L, Edgerton E S, et al. Ambient PM2.5 organic and elemental carbon in New York City: changing source contributions during a decade of large emission reductions [J]. J. Air Waste Manage Assoc., 2021,71(8):995-1012.

    [54] Ji D, Gao M, Maenhaut W, et al. The carbonaceous aerosol levels still remain a challenge in the Beijing-Tianjin-Hebei region of China: Insights from continuous high temporal resolution measurements in multiple cities [J]. Environ. Int., 2019,126:171-183.

    [55] Bond T C, Bergstrom R W. Light absorption by carbonaceous particles: An investigative review [J]. Aerosol Sci. Technol., 2006,40(1):27-67.

    [56] Xu X, Zhao W, Zhang Q, et al. Optical properties of atmospheric fine particles near Beijing during the HOPE-J3A campaign [J]. Atmos. Chem. Phys., 2016,16(10):6421-6439.

    [57] Xie C, Xu W, Wang J, et al. Light absorption enhancement of black carbon in urban Beijing in summer [J]. Atmos. Environ., 2019,213: 499-504.

    [58] Sun J, Xie C, Xu W, et al. Light absorption of black carbon and brown carbon in winter in North China Plain: comparisons between urban and rural sites [J]. Sci. Total Environ., 2021,770:144821.

    [59] Chen B, Bai Z, Cui X, et al. Light absorption enhancement of black carbon from urban haze in Northern China winter [J]. Environ. Pollut., 2017,221:418-426.

    [60] Wang Q Y, Huang R J, Cao J J, et al. Mixing state of black carbon aerosol in a heavily polluted urban area of China: Implications for light absorption enhancement [J]. Aerosol Sci. Technol., 2014,48(7): 689-697.

    [61] Xu X, Zhao W, Qian X, et al. Influence of photochemical aging on light absorption of atmospheric black carbon and aerosol single scattering albedo [J]. Atmos. Chem. Phys., 2018,18:16829-16844.

    [62] Cao F, Zhang X, Hao C, et al. Light absorption enhancement of particulate matters and their source apportionment over the Asian continental outflow site and South Yellow Sea [J]. Environ. Sci. Pollut. R., 2021,28:8022–8035.

    [63] Wang T, Zhao G, Tan T, et al. Effects of biomass burning and photochemical oxidation on the black carbon mixing state and light absorption in summer season [J]. Atmos. Environ., 2021,248:118230.

    [64] 黃 凡,陳 楠,周家斌,等.2016~2017年武漢市城區(qū)大氣PM2.5污染特征及來(lái)源解析 [J]. 中國(guó)環(huán)境監(jiān)測(cè), 2019,35(1):17-25.

    Huang F, Chen N, Zhou J B,et al. Charateristics and source apportionment of PM2.5in urban areas of Wuhan during 2016~2017 [J]. Environmental Monitoring in China, 2019,35(1):17-25.

    [65] Wang G, Zhang R, Gomez M E, et al. Persistent sulfate formation from London Fog to Chinese haze [J]. Proc. Natl. Acad. Sci. USA, 2016, 113(48):13630-13635.

    [66] Li G, Lu D, Yang X, et al. Resurgence of sandstorms complicates China’s air pollution situation [J]. Environ. Sci. Technol., 2021, 55(17):11467-11469.

    [67] Fialho P, Hansen A D A, Honrath R E. Absorption coefficients by aerosols in remote areas: a new approach to decouple dust and black carbon absorption coefficients using seven-wavelength Aethalometer data [J]. J. Aerosol Sci., 2005,36(2):267-282.

    [68] Yuan J F, Huang X F, Cao L M, et al. Light absorption of brown carbon aerosol in the PRD region of China [J]. Atmos. Chem. Phys., 2016,16(3):1433-1443.

    [69] Yang M, Howell S G, Zhuang J, et al. Attribution of aerosol light absorption to black carbon, brown carbon, and dust in China - interpretations of atmospheric measurements during EAST-AIRE [J]. Atmos. Chem. Phys., 2009,9(6):2035-2050.

    [70] Mahowald N, Albani S, Kok J F, et al. The size distribution of desert dust aerosols and its impact on the Earth system [J]. Aeolian Research, 2014,15:53-71.

    [71] Mori I, Nishikawa M, Tanimura T, et al. Change in size distribution and chemical composition of kosa (Asian dust) aerosol during long-range transport [J]. Atmos. Environ., 2003,37(30):4253-63.

    [72] Lack D A, Cappa C D. Impact of brown and clear carbon on light absorption enhancement, single scatter albedo and absorption wavelength dependence of black carbon [J]. Atmos. Chem. Phys., 2010,10(9):4207-4220.

    [73] Fierce L, Onasch T B, Cappa C D, et al. Radiative absorption enhancements by black carbon controlled by particle-to-particle heterogeneity in composition [J]. Proc. Natl. Acad. Sci. USA, 2020, 117(10):5196-5203.

    [74] Zhai J, Yang X, Li L, et al. Absorption Enhancement of Black Carbon Aerosols Constrained by Mixing-State Heterogeneity [J]. Environ. Sci. Technol., 2022,56(3):1586-1593.

    Determination of crustal element iron light absorption contribution and effects on black carbon light absorption enhancement estimation: a case study in Wuhan.

    TAN Jian1,2,3, XIA Rui1,2,3, WU Dui1,2,3, KONG Shao-fei4,5, CHEN Nan5,6*, CHENG Chun-lei1,2,3, DENG Tao7, TAO Li-ping1,2,3, ZHANG Xue1,2,3, WU Bo-xi1,2,3, WU Liang-bin1,2,3, WANG Qing1,2,3, WU Cheng1,2,3**

    (1.Institute of Mass Spectrometry and Atmospheric Environment, Jinan University, Guangzhou 510632, China;2.Guangdong-Hongkong-Macau Joint Laboratory of Collaborative Innovation for Environmental Quality,Guangzhou 510632, China;3.Guangdong Provincial Engineering Research Center for On-line Source Apportionment System of Air Pollution, Guangzhou 510632, China;4.School of Environmental Studies, China University of Geosciences (Wuhan), Wuhan 430074, China;5.Research Centre for Complex Air Pollution of Hubei Province, Wuhan 430074, China;6.Eco-Environmental Monitoring Centre of Hubei Province, Wuhan 430074, China;7.Institute of Tropical and Marine Meteorology, China Meteorological Administration, Guangzhou 510640, China)., 2022,42(7):3033~3045

    The traditional method uses empirical absorption ?ngstr?m exponent to estimate the light absorption of crustal elements or brown carbon, and this method can only handle the two components scenario (blackcarbon/brown carbon or blackcarbon/crustal elements), leading to substantial uncertainties. This study introduces a new method that uses the Fe as the crustal element tracer. Using measurement data in Wuhan (February, March, August, and September 2021) by applying the minimum R-squared method (MRS), the crustal elements light absorption contribution can be obtained. After subtracting the crustal elements light absorption, absorption enhancement factor (abs) can be obtained. The results showed that the average crustal elements absorption contribution at 370nm during the observation period was 12.3%, with a monthly average ranging from 5.7% to 15.5%. Meanwhile, AAE exhibited a positive dependence on the concentration of Fe. These collective clues indicate that the crustal elements absorption contribution cannot be ignored. Crustal elements light absorption demonstrated a considerable seasonality, which was low in winter and high in spring. There are deviations in the calculatedabsbefore and after the separation of crustal elements light absorption, which is influenced by the correlation between blackcarbon and iron. During the observation period, after deducting the light absorption of crustal elements, the averageabswas 1.43±0.53, andabswas high in spring and summer, and low in autumn and winter. Highabsin spring was associated with a higher degree of BC aging. The analysis found thatabsdemonstrated a positive dependence on the content of inorganic salts and organic matters, which proves the influence of these coating materials on the enhancement of light absorption.

    mineral dust light absorption contribution;black carbon light absorption enhancement;minimum R-squared method (MRS);absorption ?ngstr?m exponent ()

    X513

    A

    1000-6923(2022)07-3033-13

    譚 健(1997-),男,湖南衡陽(yáng)人,暨南大學(xué)碩士研究生,主要從事黑碳?xì)馊苣z吸光增強(qiáng)方面研究.發(fā)表論文4篇.

    2021-11-22

    國(guó)家自然科學(xué)基金資助項(xiàng)目(41605002,41475004);廣東省科技創(chuàng)新戰(zhàn)略專(zhuān)項(xiàng)資金項(xiàng)目(2019B121205004);國(guó)家重點(diǎn)研發(fā)計(jì)劃(2019YFB2102904)

    * 責(zé)任作者, 正高級(jí)工程師, nanchen2020@126.com; ** 副研究員, wucheng.vip@foxmail.com

    猜你喜歡
    光吸收吸光依賴(lài)性
    金色的吸管
    非等熵 Chaplygin氣體極限黎曼解關(guān)于擾動(dòng)的依賴(lài)性
    關(guān)于N—敏感依賴(lài)性的迭代特性
    商情(2017年38期)2017-11-28 14:08:59
    N-月桂?;劝彼猁}性能的pH依賴(lài)性
    多元稀土硼化物Ce1-xNdxB6的制備及性能研究*
    功能材料(2016年1期)2016-05-17 03:38:24
    你把空氣全部吸光
    初中生(2016年1期)2016-04-14 20:13:15
    鑲嵌納米晶硅的氧化硅薄膜微觀(guān)結(jié)構(gòu)調(diào)整及其光吸收特性
    舒適護(hù)理在面部激素依賴(lài)性皮炎中的應(yīng)用
    半菁染料作為染料敏化太陽(yáng)能電池吸光材料的理論研究
    既能吸光又能發(fā)光的太陽(yáng)能電池
    在线观看免费视频日本深夜| 亚洲熟妇熟女久久| 亚洲av免费高清在线观看| 亚洲不卡免费看| 人人妻人人看人人澡| 一级黄色大片毛片| 日韩欧美一区二区三区在线观看| 欧美成人一区二区免费高清观看| 国产真实乱freesex| 久久精品国产亚洲av天美| 国产精品亚洲美女久久久| av在线蜜桃| 一进一出好大好爽视频| 国产成人影院久久av| 久久草成人影院| 亚洲精华国产精华液的使用体验 | 成人欧美大片| 亚洲久久久久久中文字幕| 色综合亚洲欧美另类图片| 啦啦啦韩国在线观看视频| 99久久无色码亚洲精品果冻| 亚洲内射少妇av| 美女高潮喷水抽搐中文字幕| 91麻豆精品激情在线观看国产| 久久久久久大精品| 韩国av在线不卡| 国产精品野战在线观看| 一区二区三区激情视频| 国产私拍福利视频在线观看| 99久久精品国产国产毛片| 直男gayav资源| 两性午夜刺激爽爽歪歪视频在线观看| 日韩精品有码人妻一区| 免费搜索国产男女视频| 99久国产av精品| 亚洲无线观看免费| 99国产极品粉嫩在线观看| 久久久久久国产a免费观看| 真人一进一出gif抽搐免费| 十八禁国产超污无遮挡网站| 国产精品久久久久久久电影| 床上黄色一级片| 日韩亚洲欧美综合| 亚洲中文日韩欧美视频| 日韩欧美精品v在线| 俄罗斯特黄特色一大片| 久久久国产成人精品二区| 韩国av一区二区三区四区| 日韩欧美国产一区二区入口| 国产伦人伦偷精品视频| 又爽又黄无遮挡网站| 久久久久免费精品人妻一区二区| avwww免费| 欧美zozozo另类| 极品教师在线视频| av福利片在线观看| 美女高潮喷水抽搐中文字幕| 亚洲av第一区精品v没综合| 男女视频在线观看网站免费| 99久久无色码亚洲精品果冻| 在线免费十八禁| 日韩av在线大香蕉| 亚洲精品日韩av片在线观看| 日本与韩国留学比较| 美女高潮喷水抽搐中文字幕| 欧美日韩乱码在线| 简卡轻食公司| 午夜久久久久精精品| 中文字幕高清在线视频| h日本视频在线播放| 国产高潮美女av| www.色视频.com| or卡值多少钱| 免费观看的影片在线观看| 极品教师在线视频| 搡老熟女国产l中国老女人| 成年女人看的毛片在线观看| 午夜福利18| 国产精品爽爽va在线观看网站| 嫁个100分男人电影在线观看| 日本撒尿小便嘘嘘汇集6| 久久精品国产亚洲av天美| 亚洲第一电影网av| 亚洲国产精品久久男人天堂| 观看免费一级毛片| 国产综合懂色| 欧美三级亚洲精品| 十八禁网站免费在线| 别揉我奶头~嗯~啊~动态视频| 免费观看精品视频网站| 亚洲国产精品sss在线观看| 国产高清有码在线观看视频| 老司机福利观看| 国产乱人伦免费视频| 国产精品一及| 非洲黑人性xxxx精品又粗又长| 91在线观看av| 国产主播在线观看一区二区| 午夜福利视频1000在线观看| 赤兔流量卡办理| 国产精品久久久久久亚洲av鲁大| 观看美女的网站| 成熟少妇高潮喷水视频| 精品不卡国产一区二区三区| АⅤ资源中文在线天堂| bbb黄色大片| 国产亚洲精品综合一区在线观看| 日韩精品有码人妻一区| 成人性生交大片免费视频hd| 搞女人的毛片| 丰满的人妻完整版| 亚洲成人久久爱视频| 尾随美女入室| 在线观看一区二区三区| 麻豆成人午夜福利视频| 久久久久久久久久成人| 免费av毛片视频| 给我免费播放毛片高清在线观看| 国产伦精品一区二区三区四那| 亚洲精品亚洲一区二区| 99国产精品一区二区蜜桃av| 国产 一区 欧美 日韩| 两个人的视频大全免费| 日本免费a在线| 一卡2卡三卡四卡精品乱码亚洲| 波野结衣二区三区在线| 亚洲熟妇熟女久久| 色视频www国产| 亚洲va日本ⅴa欧美va伊人久久| 日韩精品有码人妻一区| 18禁黄网站禁片午夜丰满| 日日摸夜夜添夜夜添av毛片 | 乱系列少妇在线播放| 变态另类成人亚洲欧美熟女| 久久精品国产亚洲av香蕉五月| 国模一区二区三区四区视频| 日日摸夜夜添夜夜添小说| 久久欧美精品欧美久久欧美| 久久久久久久精品吃奶| 看片在线看免费视频| 亚洲成人精品中文字幕电影| a在线观看视频网站| 国产探花在线观看一区二区| 五月伊人婷婷丁香| 一本一本综合久久| 老女人水多毛片| 欧美绝顶高潮抽搐喷水| 成人国产麻豆网| 搡老妇女老女人老熟妇| 日本黄色片子视频| 成人特级黄色片久久久久久久| 男女边吃奶边做爰视频| 99热这里只有精品一区| 亚洲 国产 在线| 男女之事视频高清在线观看| 精品人妻一区二区三区麻豆 | 精品人妻视频免费看| 亚洲五月天丁香| 国产探花在线观看一区二区| 国产午夜精品久久久久久一区二区三区 | 搡女人真爽免费视频火全软件 | 97超级碰碰碰精品色视频在线观看| 亚洲最大成人中文| 女生性感内裤真人,穿戴方法视频| 又爽又黄无遮挡网站| 欧美日韩综合久久久久久 | 一本精品99久久精品77| 午夜免费成人在线视频| 一本久久中文字幕| 国产探花极品一区二区| 亚洲内射少妇av| 蜜桃亚洲精品一区二区三区| 人妻夜夜爽99麻豆av| 最近视频中文字幕2019在线8| 国产久久久一区二区三区| 午夜福利在线观看免费完整高清在 | 日韩人妻高清精品专区| 午夜老司机福利剧场| 人人妻人人看人人澡| 精品国产三级普通话版| 国产中年淑女户外野战色| 欧美不卡视频在线免费观看| 欧美精品国产亚洲| 精品久久久久久成人av| 国产男靠女视频免费网站| 国产精品久久视频播放| 国产免费男女视频| 国产精品一区二区三区四区久久| 九九爱精品视频在线观看| 免费在线观看影片大全网站| 啪啪无遮挡十八禁网站| 日日夜夜操网爽| 美女大奶头视频| 午夜福利视频1000在线观看| 黄色一级大片看看| 日韩强制内射视频| 亚洲第一电影网av| 日本色播在线视频| 热99re8久久精品国产| 欧美激情久久久久久爽电影| 久久久久久大精品| 午夜精品一区二区三区免费看| 色在线成人网| 久久国内精品自在自线图片| 国产精品一区二区三区四区免费观看 | 久久久精品欧美日韩精品| 久久精品国产鲁丝片午夜精品 | 欧美极品一区二区三区四区| 国产一级毛片七仙女欲春2| 俄罗斯特黄特色一大片| 九色国产91popny在线| 日日啪夜夜撸| 亚洲av免费高清在线观看| 搞女人的毛片| 久久久久久大精品| 久久精品国产鲁丝片午夜精品 | 欧美国产日韩亚洲一区| 亚洲欧美日韩无卡精品| 午夜激情欧美在线| 色哟哟哟哟哟哟| 最近中文字幕高清免费大全6 | 欧美bdsm另类| 精品无人区乱码1区二区| 亚洲专区国产一区二区| 在线观看美女被高潮喷水网站| 国产又黄又爽又无遮挡在线| 午夜a级毛片| 欧美xxxx性猛交bbbb| .国产精品久久| 九九热线精品视视频播放| 亚洲第一电影网av| 少妇被粗大猛烈的视频| x7x7x7水蜜桃| 欧美日韩综合久久久久久 | 日日摸夜夜添夜夜添小说| 亚洲美女视频黄频| 性欧美人与动物交配| 亚洲中文字幕一区二区三区有码在线看| 嫩草影院新地址| 免费观看精品视频网站| 国产精品,欧美在线| 久久99热6这里只有精品| 亚洲第一区二区三区不卡| 欧美日韩乱码在线| 免费人成在线观看视频色| 男女做爰动态图高潮gif福利片| 国内精品久久久久久久电影| 日本黄大片高清| 岛国在线免费视频观看| 少妇的逼水好多| 一进一出好大好爽视频| 国产精品久久电影中文字幕| 少妇人妻一区二区三区视频| 很黄的视频免费| 国产精品久久久久久久电影| а√天堂www在线а√下载| 国产探花在线观看一区二区| 97超视频在线观看视频| 精品久久久久久久久久久久久| 亚洲乱码一区二区免费版| 精品欧美国产一区二区三| 欧美激情久久久久久爽电影| 校园人妻丝袜中文字幕| 亚洲三级黄色毛片| 亚洲乱码一区二区免费版| 变态另类丝袜制服| 一区二区三区高清视频在线| 精品久久久久久久久亚洲 | 亚洲av中文字字幕乱码综合| 三级国产精品欧美在线观看| 他把我摸到了高潮在线观看| 国产精华一区二区三区| 亚洲午夜理论影院| 一区二区三区高清视频在线| 久久久精品欧美日韩精品| 成人亚洲精品av一区二区| 国产精品98久久久久久宅男小说| 久久久久久大精品| 亚洲av免费高清在线观看| 乱码一卡2卡4卡精品| 99国产极品粉嫩在线观看| 午夜免费激情av| 天堂网av新在线| 国产三级中文精品| 我要搜黄色片| 亚洲国产精品久久男人天堂| 亚洲av五月六月丁香网| www.色视频.com| 国产成人aa在线观看| 欧美最黄视频在线播放免费| 亚洲av二区三区四区| 欧美黑人欧美精品刺激| 99热6这里只有精品| av福利片在线观看| 97超视频在线观看视频| 成人性生交大片免费视频hd| 精品久久久久久久人妻蜜臀av| 国产高潮美女av| 久久午夜福利片| 亚洲男人的天堂狠狠| 一区二区三区高清视频在线| 99九九线精品视频在线观看视频| 91精品国产九色| 国产精品国产高清国产av| 精品人妻视频免费看| 99在线视频只有这里精品首页| 99精品在免费线老司机午夜| 亚洲欧美日韩无卡精品| 日本 欧美在线| a级一级毛片免费在线观看| 成人特级黄色片久久久久久久| 淫秽高清视频在线观看| 免费黄网站久久成人精品| 身体一侧抽搐| 少妇猛男粗大的猛烈进出视频 | 真实男女啪啪啪动态图| 国产亚洲精品av在线| 九九久久精品国产亚洲av麻豆| 中文字幕精品亚洲无线码一区| 亚洲va日本ⅴa欧美va伊人久久| 俺也久久电影网| 精品午夜福利视频在线观看一区| 国产真实伦视频高清在线观看 | 99国产极品粉嫩在线观看| 成人av在线播放网站| 亚洲国产高清在线一区二区三| 毛片一级片免费看久久久久 | 九色成人免费人妻av| 国产视频内射| 亚洲av中文av极速乱 | 春色校园在线视频观看| 日韩,欧美,国产一区二区三区 | 亚洲经典国产精华液单| 久久亚洲真实| 免费av不卡在线播放| 丝袜美腿在线中文| 久久精品国产亚洲av涩爱 | 国产av不卡久久| 九色国产91popny在线| 亚洲综合色惰| 乱码一卡2卡4卡精品| 国产免费av片在线观看野外av| 成人亚洲精品av一区二区| 美女大奶头视频| 欧美3d第一页| 亚洲乱码一区二区免费版| 久久精品国产亚洲网站| 日本五十路高清| 三级毛片av免费| 欧美日韩黄片免| 免费观看的影片在线观看| 亚洲专区国产一区二区| 嫩草影院入口| 亚洲人成网站在线播放欧美日韩| 成人国产一区最新在线观看| 亚洲成人久久爱视频| 一区二区三区高清视频在线| 日韩 亚洲 欧美在线| 成人av一区二区三区在线看| 成人av在线播放网站| 成人av一区二区三区在线看| 欧美xxxx性猛交bbbb| 在线观看av片永久免费下载| 午夜激情欧美在线| 午夜影院日韩av| 夜夜夜夜夜久久久久| av天堂在线播放| 一级黄片播放器| 中文字幕av成人在线电影| 午夜爱爱视频在线播放| 男女啪啪激烈高潮av片| 中文字幕精品亚洲无线码一区| 欧美+亚洲+日韩+国产| 18禁黄网站禁片午夜丰满| 日韩人妻高清精品专区| 欧美色欧美亚洲另类二区| 久久欧美精品欧美久久欧美| 免费在线观看成人毛片| 深夜a级毛片| 91精品国产九色| 日韩欧美精品v在线| 亚洲精品456在线播放app | 国产高清三级在线| 变态另类丝袜制服| 免费av观看视频| 国产一区二区三区在线臀色熟女| 特级一级黄色大片| 99国产精品一区二区蜜桃av| 亚洲av二区三区四区| 精品久久久久久久末码| 国产精品久久久久久久电影| av女优亚洲男人天堂| 精品一区二区三区人妻视频| 可以在线观看的亚洲视频| 日韩人妻高清精品专区| 国产69精品久久久久777片| 深夜a级毛片| 精品日产1卡2卡| 在线观看一区二区三区| 老师上课跳d突然被开到最大视频| 精品乱码久久久久久99久播| 国产色婷婷99| 日韩强制内射视频| 国产精品亚洲美女久久久| 一夜夜www| 亚洲,欧美,日韩| 午夜亚洲福利在线播放| 网址你懂的国产日韩在线| 欧美黑人巨大hd| 白带黄色成豆腐渣| 一区二区三区免费毛片| 国产精品98久久久久久宅男小说| 久久久久九九精品影院| 亚洲专区中文字幕在线| АⅤ资源中文在线天堂| 色视频www国产| 亚洲色图av天堂| 亚洲18禁久久av| 非洲黑人性xxxx精品又粗又长| 看十八女毛片水多多多| 我的女老师完整版在线观看| 国产日本99.免费观看| 动漫黄色视频在线观看| 精品免费久久久久久久清纯| 99久久精品一区二区三区| 黄色女人牲交| 夜夜夜夜夜久久久久| 国产高清视频在线播放一区| 成人无遮挡网站| 欧美人与善性xxx| 亚洲内射少妇av| 国产综合懂色| 欧美bdsm另类| xxxwww97欧美| 国产乱人视频| 久久九九热精品免费| 亚洲经典国产精华液单| 超碰av人人做人人爽久久| 99久久精品一区二区三区| 日韩欧美三级三区| 黄色丝袜av网址大全| 亚洲精品日韩av片在线观看| 久久精品国产亚洲网站| 日本黄大片高清| 一本久久中文字幕| 色综合色国产| 超碰av人人做人人爽久久| 免费观看人在逋| 亚洲欧美日韩高清在线视频| 老司机福利观看| 婷婷色综合大香蕉| 中文字幕免费在线视频6| 国产一区二区在线观看日韩| 男人狂女人下面高潮的视频| 亚洲真实伦在线观看| 午夜爱爱视频在线播放| 精品一区二区三区人妻视频| 久久天躁狠狠躁夜夜2o2o| 国产乱人伦免费视频| 亚洲av第一区精品v没综合| 国产精品99久久久久久久久| 九色成人免费人妻av| 久久久精品大字幕| 欧美黑人欧美精品刺激| 国产 一区 欧美 日韩| 我的女老师完整版在线观看| 真人一进一出gif抽搐免费| 国产精品一区二区三区四区免费观看 | 欧美日韩亚洲国产一区二区在线观看| 欧美+亚洲+日韩+国产| 男插女下体视频免费在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 变态另类丝袜制服| 观看免费一级毛片| 欧美日韩精品成人综合77777| 久久午夜亚洲精品久久| 亚洲欧美日韩东京热| 国产一区二区三区在线臀色熟女| 十八禁国产超污无遮挡网站| 国产一区二区在线观看日韩| 人妻少妇偷人精品九色| 欧美不卡视频在线免费观看| 国产欧美日韩精品亚洲av| 九九爱精品视频在线观看| 国产精品久久久久久亚洲av鲁大| 性插视频无遮挡在线免费观看| 在线播放国产精品三级| 91av网一区二区| 97热精品久久久久久| 国产激情偷乱视频一区二区| 精品人妻视频免费看| 麻豆成人午夜福利视频| 黄色一级大片看看| 亚洲精品国产成人久久av| 精品久久久久久久久久久久久| 国内少妇人妻偷人精品xxx网站| 欧美在线一区亚洲| 看黄色毛片网站| АⅤ资源中文在线天堂| 日韩欧美在线二视频| 亚洲自拍偷在线| 噜噜噜噜噜久久久久久91| 又爽又黄a免费视频| 久久欧美精品欧美久久欧美| 一个人观看的视频www高清免费观看| 在线国产一区二区在线| 色综合亚洲欧美另类图片| 欧美丝袜亚洲另类 | 国产黄色小视频在线观看| 国产在线精品亚洲第一网站| av中文乱码字幕在线| 国产精品久久视频播放| 欧美人与善性xxx| www.www免费av| 日韩欧美国产一区二区入口| 亚洲人成网站在线播| 久久人人爽人人爽人人片va| 国产午夜精品论理片| 久久午夜福利片| 男女啪啪激烈高潮av片| 欧洲精品卡2卡3卡4卡5卡区| 十八禁网站免费在线| 欧美绝顶高潮抽搐喷水| 午夜免费激情av| av在线观看视频网站免费| 亚洲美女搞黄在线观看 | 精品人妻熟女av久视频| 成人三级黄色视频| 淫妇啪啪啪对白视频| 免费大片18禁| 一区二区三区激情视频| 国产69精品久久久久777片| 我的女老师完整版在线观看| 在线观看舔阴道视频| 一区二区三区免费毛片| 嫩草影视91久久| 国产久久久一区二区三区| 色哟哟·www| 色综合色国产| 18禁黄网站禁片午夜丰满| 一区福利在线观看| 99热精品在线国产| 两个人视频免费观看高清| 精品一区二区三区人妻视频| 91麻豆av在线| 亚洲精品粉嫩美女一区| 国产白丝娇喘喷水9色精品| 国产精品美女特级片免费视频播放器| 久久精品夜夜夜夜夜久久蜜豆| 国产探花在线观看一区二区| 99国产精品一区二区蜜桃av| 久久久久久九九精品二区国产| 亚洲人成网站在线播| 国产一区二区亚洲精品在线观看| 黄色日韩在线| 美女 人体艺术 gogo| 男女做爰动态图高潮gif福利片| 成人精品一区二区免费| 性插视频无遮挡在线免费观看| 亚洲精品亚洲一区二区| 别揉我奶头 嗯啊视频| 有码 亚洲区| 亚洲av二区三区四区| 很黄的视频免费| 色吧在线观看| 欧美绝顶高潮抽搐喷水| 麻豆成人午夜福利视频| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕av在线有码专区| 国模一区二区三区四区视频| 村上凉子中文字幕在线| 亚洲欧美日韩高清在线视频| 91久久精品国产一区二区成人| 亚洲 国产 在线| 又爽又黄a免费视频| 免费观看人在逋| 日韩欧美 国产精品| 亚洲国产精品久久男人天堂| 欧美性猛交╳xxx乱大交人| 国产熟女欧美一区二区| 91麻豆av在线| a级毛片免费高清观看在线播放| 欧美最新免费一区二区三区| 韩国av在线不卡| 国产精品人妻久久久影院| 国产精品国产高清国产av| 中文字幕精品亚洲无线码一区| 亚洲一级一片aⅴ在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 简卡轻食公司| 高清日韩中文字幕在线| 欧美激情国产日韩精品一区| 日韩亚洲欧美综合| 一区福利在线观看| 亚洲国产日韩欧美精品在线观看| 亚洲四区av| 黄片wwwwww| 国产亚洲精品综合一区在线观看| 午夜福利视频1000在线观看| 一级a爱片免费观看的视频| 真人做人爱边吃奶动态| 一本久久中文字幕| 3wmmmm亚洲av在线观看| 日韩在线高清观看一区二区三区 | 欧美最新免费一区二区三区| 成人鲁丝片一二三区免费| 国产av不卡久久| 麻豆成人av在线观看| 99久国产av精品| 男女边吃奶边做爰视频| 一级黄色大片毛片| 91精品国产九色| or卡值多少钱| 人妻夜夜爽99麻豆av| 国产久久久一区二区三区|