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

    基于CMAQ模式的中國大氣BaP遷移轉(zhuǎn)化模擬研究

    2016-10-13 19:45:10閆文君劉世杰何友江
    中國環(huán)境科學 2016年6期
    關(guān)鍵詞:氣相大氣顆粒

    閆文君,劉 敏,劉世杰,何友江,唐 偉,汪 青,李 曄,孟 凡*

    ?

    基于CMAQ模式的中國大氣BaP遷移轉(zhuǎn)化模擬研究

    閆文君1,2,劉 敏1,劉世杰2,何友江2,唐 偉2,汪 青1,李 曄1,孟 凡2*

    (1.華東師范大學地理科學學院,地理信息科學教育部重點實驗室,上海 200241;2.中國環(huán)境科學研究院大氣環(huán)境研究所,北京 100012)

    在空氣質(zhì)量模式CMAQ中增加了對大氣中苯并[a]芘(BaP)傳輸、轉(zhuǎn)化和沉降過程的模擬,新增了BaP氣相/顆粒相間分配模式(OA吸收模式)以及BaP化學轉(zhuǎn)換,并采用該模型模擬了2014年1、4、7、10月四個季節(jié)代表月中國大氣中BaP的濃度以及干濕沉降通量的空間分布特征.模擬結(jié)果與部分地區(qū)的實測數(shù)據(jù)進行了初步對比驗證,結(jié)果顯示BaP濃度的模擬值和測量值的量級基本一致,其比值基本在1~2個量級范圍內(nèi),表明模擬結(jié)果處于可接受的范圍.基于模型結(jié)果,還探討了BaP濃度與沉降季節(jié)變化特征與人為源排放的緊密聯(lián)系,說明濕沉降是BaP最主要的沉降途徑.BaP沉降的空間分布與其排放源和濃度具有密切的聯(lián)系,并且具有長距離輸送的特征等.

    CMAQ模式;BaP;季節(jié)變化;干沉降;濕沉降

    苯并[a]芘(BaP)是苯與芘稠合而成的一類五環(huán)的多環(huán)芳烴(PAHs),一般存在于機動車尾氣排放、住宅木材燃燒的濃煙、燃煤電廠的飛灰以及其他燃燒相關(guān)過程產(chǎn)生的可吸入粒子中[1].BaP作為第一個被發(fā)現(xiàn)的環(huán)境化學致癌物,具有很強的致癌特性,一般將它作為PAHs的代表.由于其可以在大氣中進行長距離輸送,導致在大陸范圍內(nèi)廣泛分布,同時其毒性特征導致存在嚴重的生態(tài)風險和人體健康風險,因此國際社會開始呼吁減少和消除PAHs排放的行動,并開展了大量相關(guān)研究.目前,EMEP(歐洲監(jiān)測和評估項目)下的國際合作項目MSC-E(氣象綜合中心—東部)已經(jīng)開始對PAHs等各種有毒物質(zhì)在歐洲的遠距離輸送以及沉降和轉(zhuǎn)化展開了相關(guān)模擬研究[2].而我國還沒有建立關(guān)于PAHs的監(jiān)測網(wǎng)絡(luò),數(shù)據(jù)獲取比較困難,不足以支撐PAHs的大區(qū)域管理控制,因此利用模式模擬PAHs等污染物在大氣中的傳輸、轉(zhuǎn)化和沉降刻不容緩.

    目前美國EPA開發(fā)并逐漸完善的MODELS-3/CMAQ模式(由中尺度氣象模式WRF、排放處理模式SMOKE和多尺度空氣質(zhì)量模式CMAQ組成)已全面考慮了常規(guī)氣態(tài)污染物和部分有毒空氣污染物,涉及氣溶膠、光化學氧化和酸沉降等多種污染問題,幾乎能夠模擬大氣中所有的物理化學過程[3],同時也考慮了實際大氣中不同物種之間的相互轉(zhuǎn)化和互相影響,不僅用于各種大氣污染物濃度的預(yù)測,而且可用于環(huán)境評估和環(huán)境控制決策[4-5]-[].

    CMAQ模式系統(tǒng)對PAHs的模擬在歐美等地已經(jīng)被廣泛應(yīng)用并進行了驗證,國外的研究重點主要集中在PAHs及其衍生物的來源解析與模擬[6]、大氣濃度與沉降量[7]、氣相—顆粒相間分配[8]以及煙霧箱模擬等.2008年Matthias等[9]第一次展示出了歐洲等地大氣中BaP的沉積域,表明CMAQ已經(jīng)擴展了模擬POPs的能力.而Meng等[10]和Zhang等[11]在CMAQ模式中增加了二類污染物(CMAQ-PCDD/Fs)氣相—顆粒相間分配機制以及二類污染物的化學過程,并用來模擬北美地區(qū)PCBs和PCDD/Fs的濃度和干濕沉降量.另外Bewersdorff等[12]分析了船舶排放對BaP的濃度和沉降分布的影響,并針對主要的船舶排放物NO,SO2和PM10對BaP濃度和沉降的影響進行了敏感性分析.

    國內(nèi)對大氣中PAHs的研究主要集中于來源、濃度、沉降以及氣相—顆粒相間分配[13]的實際觀測,對PAHs的模擬研究則多集中于基于逸度模型的多介質(zhì)歸趨模擬研究[14-15]-[]等.例如Tao等[16]利用改進的多介質(zhì)模型在自己建立的1km分辨率的排放清單基礎(chǔ)上模擬了天津地區(qū)多環(huán)芳烴的擴散遷移.Liu等[17]運用網(wǎng)格化的區(qū)域尺度上的多介質(zhì)逸度模型模擬了渤海灣地區(qū)BaP的遷移轉(zhuǎn)化,并對輸入?yún)?shù)做了敏感性測試.在CMAQ模式的應(yīng)用方面,國內(nèi)研究者目前主要著眼于對常規(guī)污染物CO、SO2、O3、NO、PM2.5等在大氣中傳輸與擴散以及能見度問題的研究,對POPs(持久性有機污染物,PAHs為其中一種)的模擬研究較少.目前僅有張玨等[18]利用改進的CMAQ模式模擬了2006年長三角地區(qū)PCDD/Fs在大氣中的輸送、轉(zhuǎn)化和沉降等演變過程.

    因此,本研究在CMAQ模式中增加了對大氣中BaP傳輸、轉(zhuǎn)化和沉降過程的模擬,新增了BaP氣相/顆粒相間分配模式(KOA吸收模式)以及BaP化學轉(zhuǎn)化,并采用該模式模擬了2014年1、4、7、10月四個季節(jié)代表月中國大氣中BaP的濃度以及干濕沉降通量的空間分布特征,以此來研究中國大氣中BaP的遷移轉(zhuǎn)化規(guī)律等.

    1 研究方法

    本研究采用CMAQ模式4.7版本,其中氣相化學反應(yīng)機制主要選擇SAPRC99化學機制[19].模擬區(qū)域采用Lambert地圖投影,投影中心坐標為35°N,110°E,標準緯線為24°N和46°N,覆蓋整個東亞地區(qū).CMAQ模式使用基于地形追隨的質(zhì)量垂直坐標系以及一重單相嵌套網(wǎng)格方案.網(wǎng)格的水平分辨率為36km′36km,網(wǎng)格數(shù)為200′160.

    1.1 氣象場及初始邊界條件

    在WRF模式中,氣象場的初始和邊界條件采用NCAR和NCEP的時間間隔為6h、空間分辨率為1°′1°的FNL再分析氣象數(shù)據(jù),地形和下墊面資料來自USGS 30s全球地形和MODIS下墊面資料.WRF的模擬結(jié)果可以通過CMAQ模式的氣象前處理模塊MCIP(氣象—化學處理端口)轉(zhuǎn)換為CMAQ可識別接收的格式.

    CMAQ模式中設(shè)置的初始和邊界條件在模式結(jié)果中有著重要的作用[20].在本研究中, CCTM(CMAQ化學—傳輸模塊)按月依次運行30d或31d,并且將每個月第1d的前1d作為初始場.敏感性測試表明,24h的自旋時間能夠使得初始條件的影響保持到最低[9].在整個區(qū)域里,將與CMAQ源代碼一起提供CO, O3,NO2,和SO2等氣相物種的典型配置文件作為邊界條件.

    1.2 排放數(shù)據(jù)

    本研究中,CMAQ模式所需的常規(guī)污染物排放清單國內(nèi)部分采用MEIC(中國多尺度排放清單模型)(2010年),國外部分采用REAS(亞洲區(qū)域排放清單)(2010年);BaP排放清單采用Shen等[21]于2013年公布的1960~2008年期間16種PAHs排放清單(2008年).利用SMOKE前處理模式對排放清單數(shù)據(jù)進行處理,將常規(guī)污染物排放清單和BaP排放清單進行整合,以模擬2014年1、4、7、10月四個季節(jié)代表月BaP的遷移轉(zhuǎn)化情況.SMOKE模式能夠?qū)⑴欧旁辞鍐螖?shù)據(jù)按照CMAQ模式所要求的時空精度和格式進行處理.通過SMOKE模式處理得到的BaP排放源強分布如圖1所示.由圖1可知其中華北平原、京津冀、長三角、珠三角等地區(qū)源排放較為集中.

    2 模式發(fā)展

    本研究基于CMAQ模式發(fā)展了BaP的相關(guān)算法,根據(jù)PAHs具有半揮發(fā)性的特征,在模式中添加了2種BaP物種形態(tài)[22]:氣相(BaP)和顆粒相(PBaP).在CMAQ模式中,氣溶膠或大氣中懸浮的顆粒物,根據(jù)粒徑大小可以分為3種模態(tài):愛根核膜(粒徑<0.1μm),積聚膜(粒徑為0.1~2.5μm),粗膜(粒徑為2.5~10μm).顆粒相BaP往往吸附在粒徑為0.1~2mm的粒子上[23],99%的BaP以積聚膜的形式存在,1%以愛根核膜的形態(tài)存在[5],因此假設(shè)模式中顆粒相BaP只以積聚膜的形態(tài)存在,并且默認該情況符合所有的有機氣溶膠.

    2.1 氣相/顆粒相間分配模式(KOA吸收模式)

    大氣中氣相和顆粒相之間的分配是POPs環(huán)境行為的一項重要特性.BaP在大氣中主要以顆粒相的形式存在[24].因此本研究分析了顆粒相BaP的濃度及干濕沉降通量.研究表明,通過吸附或吸收產(chǎn)生的平衡過程[25]或者動態(tài)傳質(zhì)過程[26]可以描述氣相、顆粒相在大氣中的分配.下面簡要地描述了本文所應(yīng)用的分配模式:OA分配模式[27].

    首先,物質(zhì)的氣相和顆粒相在大氣中的分配系數(shù)P可以估計為:

    式中:p是物質(zhì)顆粒相的濃度,ng/m3;g是物質(zhì)氣相的濃度,ng/m3;TSP是大氣中總懸浮顆粒物的濃度,μg/m3.

    根據(jù)Junge-Pankow吸附模式[28-29]可知,大氣中POPs在氣相和顆粒相之間的分配與辛醇—空氣分配系數(shù)OA直接相關(guān).Finizio等[30]對此提出了下列表達式:

    式中:om是粒子吸附的有機物所占的質(zhì)量分數(shù);o和om分別是辛醇和有機物的分子量,g/mol;o和om分別是辛醇和有機物中PAHs的活度系數(shù).Harner和Bidleman運用此模式提出了大氣中PAHs的P對KOA的回歸相依性,用方程表示如下[31-32]-[]:

    式中:和是經(jīng)驗常數(shù).采用Odabasi等[33]通過實驗得出的值:=?6.50,= 5382.

    2.2 BaP的去除過程

    模式考慮了BaP在大氣中的兩種去除過程:物理去除過程和化學轉(zhuǎn)化去除過程.BaP所依附的粒子發(fā)生干沉降時,它的化學反應(yīng)半衰期大約為10d[34],此時只考慮物理去除過程.夏季降水較多,濕沉降量較大,此時受到污染的城市大氣中BaP的半衰期僅為1.93h[35],說明大氣中的BaP被濕沉降過程有效地去除,這可能也是造成夏季環(huán)境大氣中測量的BaP具有較低濃度的原因.因此物理去除過程考慮了顆粒相BaP的干沉降過程以及氣相與顆粒相BaP的濕沉降過程.化學轉(zhuǎn)化去除過程是指BaP在大氣中發(fā)生氣象化學反應(yīng)而轉(zhuǎn)化成其他物質(zhì)進行清除.氣相BaP在大氣中會與OH和NO3自由基發(fā)生反應(yīng),根據(jù)反應(yīng)的重要性,本研究只考慮了氣相BaP與OH自由基的反應(yīng)[36].而顆粒相BaP一般發(fā)生的多相反應(yīng)此處不予考慮.

    3 模擬結(jié)果與模式驗證

    3.1 BaP濃度結(jié)果分析

    由圖2可知,濃度較高的區(qū)域主要分布在華北平原、長三角、珠三角等地區(qū),1月四川盆地一帶濃度也較高.因此,從空間分布來說,BaP的濃度分布與源排放的空間分布基本相同,源的高排放成為其濃度集中分布的主要原因,同時有比較明顯的擴散和傳輸效應(yīng).但是由于模式中假設(shè)源排放量無時間變化和季節(jié)變化,因此模擬的BaP濃度的時間變化基本取決于氣象條件的變化.還可看到,BaP的濃度分布具有明顯的季節(jié)變化,1月濃度最高;4月和10月次之;7月濃度最低.例如北京1月BaP濃度模擬值最高可達17.506ng/m3,7月濃度最低到0.092ng/m3.這是因為7月夏季降水較多,大量的顆粒相BaP隨著雨水沉降到地面,對BaP濃度進行有效地清除;而1月為冬季供暖季,人為源的排放量較多,同時邊界層的混合層高度降低,天氣系統(tǒng)穩(wěn)定,使大氣中的BaP較長時間停留在近地層中,另外還需要考慮的是冬季白天日照時間短,低溫能夠減少BaP的光氧化產(chǎn)物和相關(guān)化學反應(yīng)[23],從而造成BaP的濃度相對較高.Bewersdorff等[37]研究顯示,全年住宅供暖所產(chǎn)生的BaP排放量與大氣溫度有一定的線性關(guān)系, BaP的冬季總排放量約為夏季的10倍,而冬季增加的排放量也將造成年平均排放量的增加.

    3.2 模式驗證

    由于目前國內(nèi)BaP研究相對較少,且BaP不屬于大氣常規(guī)監(jiān)測污染物,因此很難獲得大量的觀測數(shù)據(jù).為了較好地評估模式模擬大氣中BaP的性能,表1中對比了2014年1、4、7、10四個季節(jié)代表月上海地區(qū)以及2015年1月1~10日北京、濟南、上海三地BaP濃度的模擬值相對于測量值的模擬效果,以此初步對模擬結(jié)果的可靠性進行驗證.本研究中模式模擬結(jié)果與實際測量值的吻合程度釆用標準化平均偏差(NMB)、標準化平均誤差(NME)、均方根誤差(RMSE)、相關(guān)系數(shù)(CORR)來評估.NMB能反映模擬值與測量值的平均偏離程度,NME反映模擬值與測量值的平均絕對誤差,RMSE反映模擬值與測量值的偏離程度.

    表1 2014~2015年上海地區(qū)BaP濃度的模擬月均值與測量值及相關(guān)統(tǒng)計指標對比Table 1 Statistical analysis of monthly averaged of modeled BaP concentrations with measurements in 2014~2015

    注:a目前采樣得到的測量值有2014年4、7、10月上海地區(qū)以及2015年1月1~10日北京、濟南、上海三地BaP濃度值,采樣時分別在3地各設(shè)置一個具有代表性的采樣點,每24h得到一個濃度值作為日均值.測量值為當月的月均值.

    從表1中可以看出2015年1月北京、濟南、上海的模擬值與測量值相關(guān)性較高,相關(guān)系數(shù)基本處于0.3~0.6之間,同時偏離度較小,模擬值約為測量值的2倍,說明模擬效果較好.2014年10月的模擬值和測量值之間的相關(guān)性最好,4月次之,7月較差,而7月和10月的模擬值偏離測量值的程度較大,偏離度甚至達到43.509,這表明模擬值與測量值仍然有很大的不確定性.BaP模式模擬的最大不確定性來源是排放源,也可能是由于模式中初始場、邊界場、氣象場以及化學機制的不確定性造成的.同時,觀測數(shù)據(jù)的時間和空間分辨率也遠小于模式,觀測數(shù)據(jù)的誤差和采樣點代表性也可能是模式與觀測比較差異的來源.通過2015年1月1~10日北京、濟南、上海三地的BaP濃度模擬日均值與同期測量值的對比(圖3)可知,模擬值與測量值的量級基本一致,模擬效果較好,但模擬值整體略高于實際測量值.研究表明,模擬結(jié)果的不確定性主要來自于排放清單的不確定性.由于排放清單的排放因子和活動數(shù)據(jù)有很大的不確定性,可能造成BaP排放清單的不確定性.一般來講排放清單的不確定性在2~5倍之間[38].但同時要考慮實驗條件的限制、采樣點不豐富以及經(jīng)緯度稀疏的因素,用于對比的實驗測量值本身也具有一定的誤差.

    3.3 BaP沉降結(jié)果分析

    由圖5可知,干沉降的分布區(qū)域基本集中在華中、華北、長三角等人口密集地區(qū),與濃度的空間分布基本一致,這表明BaP的干沉降空間分布主要取決于它的濃度分布,并且濃度高的地區(qū)干沉降量大,濃度低的地區(qū)干沉降量則小.同時,新疆、青海、西藏等偏遠地區(qū)的干沉降量普遍接近0,說明BaP在大氣中進行了長距離的遷移輸送,使得我國中部東部等沿海地區(qū)干沉降分布較為集中,干沉降量也較大.BaP的干沉降空間分布還具有明顯的季節(jié)變化,1月的干沉降范圍最大,干沉降量也大;4月和10月次之;7月的沉降范圍最小,干沉降量也最少.這主要是因為冬季燃煤較多,人為排放量大,大氣中的顆粒相BaP也較多,且低溫時BaP的濃度較高,因此干沉降量隨著溫度的降低而增加.

    通過與濕沉降的空間分布(圖6)對比可以看出,只有降水才會出現(xiàn)BaP的濕沉降,并且BaP的沉降分布基本符合其濃度分布.模擬期間全國的濕沉降分布范圍明顯大于干沉降分布,7月為夏季,降水最多最頻繁,分布范圍最大,1月和10月濕沉降分布最少,4月次之.另外,1月也有部分地區(qū)濕沉降量較大,是因為1月大氣中顆粒物較多,BaP濃度較大,一旦發(fā)生降水,則當日濕沉降量會較大,有時甚至會遠大于干沉降量.因此BaP濕沉降的空間分布形式主要取決于它的濃度和降水分布.

    從表2可見,三地的干沉降量均表現(xiàn)為1月冬季最高,7月夏季最低,2014年上海的沉降總量高于北京和濟南,這主要是因為上海地區(qū)大氣中水汽含量較高,降水量充足,因此濕沉降量較大,約占其沉降總量的84%.說明濕沉降是BaP最主要的沉降途徑,氣象條件對沉降量具有顯著影響.而三地的干沉降量則相差不大,這主要是因為BaP在大氣中能夠進行長距離遷移輸送的緣故.

    表2 2014年1、4、7、10月北京、濟南、上海三地BaP干濕沉降通量模擬值對比(μg/m2)Table 2 Comparison between modeled BaP dry deposition rates and wet deposition rates for January, April, July, and October 2014 in Beijing, Jinan and Shanghai (μg/m2)

    BaP是一種揮發(fā)性較小的化合物,在環(huán)境大氣中適當溫度下BaP質(zhì)量的主要部分與粒子結(jié)合在一起,因此沉降主要集中在在排放源的附近.研究表明超過80%的BaP排放量與私人燃燒以及公共發(fā)電廠的熱能產(chǎn)生有關(guān)[7].因此最大沉降量及最高沉降率一般會出現(xiàn)在人口密集分布的地區(qū),人口稀少的偏遠地區(qū)沉降較少.另外,沿海地區(qū)也有明顯的沉降,且沿海和外海之間存在很大的梯度,這主要是因為CMAQ模式模擬時只包含了陸地基準排放,而忽略了海洋上的船舶排放.沿海海岸附近的顯著源區(qū)在京津冀、山東半島、上海等人口分布密集的地區(qū).BaP沉降的空間分布與BaP的排放源、濃度具有密切的聯(lián)系.

    4 結(jié)論

    4.1 改進的CMAQ模式對大氣中BaP的模擬研究結(jié)果顯示,BaP的濃度分布與源排放的空間分布基本相同,源的高排放是其濃度集中分布的主要原因.BaP的濃度分布還具有明顯的季節(jié)變化,冬季最高,夏季最低.

    4.2 通過對模擬結(jié)果的一系列對比驗證可知,2015年1月1~10日北京、濟南、上海三地的BaP濃度模擬值與測量值量級基本一致,模擬值約為測量值的2倍,模擬效果較好.2014年10月的模擬值和測量值之間的相關(guān)性最好,而7月模擬值與測量值偏離度較高,但比值基本在1~2個量級范圍內(nèi),與國內(nèi)同類研究偏離水平相當,表明模擬結(jié)果處于可接受的范圍.

    4.3 BaP的干沉降空間分布主要取決于濃度分布,并具有明顯的季節(jié)變化,1月干沉降最大,7月最小,其沉降分布在模擬區(qū)域具有長距離輸送的特征.而濕沉降作為BaP最主要的沉降途徑,其空間分布主要取決于它的濃度和降水分布.BaP沉降的空間分布多集中在中部東部沿海等人口密集地區(qū),與BaP的排放源、濃度具有密切聯(lián)系.

    4.4 由于本次研究所應(yīng)用的排放清單尚未考慮各個物種的時間序列,因此BaP模擬值的驗證需要今后更大的努力進行時間序列的建模,以此來進一步探討B(tài)aP濃度與干濕沉降通量的時間變化特征.未來的研究工作中希望加強排放清單的制作,盡量保證排放源、排放因子與活動的數(shù)據(jù)的準確性,以更高的分辨率模擬區(qū)域內(nèi)的有機物.

    唐孝炎,張遠航,邵 敏.大氣環(huán)境化學 [M]. 北京:高等教育出版社, 2006:91-101.

    Wania F, Mackay D. Tracking the distribution of persistent organic pollutants [J]. Environ. Sci. Technol., 1996,30(9):390- 396.

    Jang C J C, Chang N B. Development and application of U.S. EPA's Models-3: An integrated "one-atmosphere" third- generation air quality modeling system [J]. Journal of the Chinese Institute of Environmental Engineering, 2000,10(1):19-34.

    Byun D W, Ching J K S. Science Algorithms of the EPA Models-3Community Multiscale Air Quality (CMAQ) Modeling System [Z]. EPA-600/R-99/030, US Environmental Protection Agency. US Government Printing Office, Washington, DC.See also /http://www.epa.gov/asmdnerl/models3/doc/science/science. html.

    San José R, Pérez J L, Callén M S, et al. BaP (PAH) air quality modelling exercise over Zaragoza (Spain) using an adapted version of WRF-CMAQ model [J]. Environmental Pollution, 2013,183:151-158.

    Zhang H, Ying Q. Secondary organic aerosol from polycyclic aromatic hydrocarbons in Southeast Texas [J]. Atmospheric Environment, 2012,55:279-287.

    Matthias V, Aulinger A, Quante M. CMAQ simulations of the Benzo(a)pyrene distribution over Europe for 2000and 2001 [J]. Atmospheric Environment, 2009,43(26):4078-4086.

    Simcik M F, Franz T P, Zhang H, et al. Gas-particle partitioning of PCBs and PAHs in the Chicago urban and adjacent coastal atmosphere: states of equilibrium [J]. Environmental Science & Technology, 1998,32(2):251-257.

    Matthias V, Aulinger A, Quante M. Adapting CMAQ to investigate air pollution in North Sea coastal regions [J]. Environmental Modelling & Software, 2008,23(3):356-368.

    Meng F, Zhang B, Gbor P, et al. Models for gas/particle partitioning, transformation and air/water surface exchange of PCBs and PCDD/Fs in CMAQ [J]. Atmospheric Environment, 2007,41(39):9111-9127.

    Zhang B, Meng F, Shi C, et al. Modeling the atmospheric transport and deposition of polychlorinated dibenzo-p-dioxins and dibenzofurans in North America [J]. Atmospheric Environment, 2009,43(13):2204-2212.

    Bewersdorff I, Aulinger A, Volker M, et al. Modeling of benzo(a)pyrene concentrations in north sea coastal areas: contribution of ship emissions [J]. Hrvatski Meteorolo?ki ?asopis, 2008,43(43/2):586-590.

    Wang W, Simonich S L M, Wang W, et al. Atmospheric polycyclic aromatic hydrocarbon concentrations and gas/particle partitioning at background, rural village and urban sites in the North China Plain [J]. Atmospheric Research, 2011,99(2):197- 206.

    汪祖丞,劉 敏,楊 毅,等.上海城區(qū)多環(huán)芳烴的多介質(zhì)歸趨模擬研究 [J]. 中國環(huán)境科學, 2011,31(6):984-990.

    劉 敏,王蘭輝.城市地表系統(tǒng)POPs多介質(zhì)循環(huán)過程與模擬研究進展 [J]. 中國環(huán)境科學, 2013,33(11):2018-2026.

    Tao S, Li X, Yang Y, et al. Dispersion modeling of polycyclic aromatic hydrocarbons from combustion of biomass and fossil fuels and production of coke in Tianjin, China [J]. Environmental Science & Technology, 2006,40(15):4586-4591.

    Liu S, Lu Y, Wang T, et al. Using gridded multimedia model to simulate spatial fate of Benzo[α]pyrene on regional scale [J]. Environment international, 2014,63:53-63.

    張 玨,孟 凡,何友江,等.長江三角洲地區(qū)大氣二噁英類污染物輸送—沉降模擬研究 [J]. 環(huán)境科學研究, 2011,24(12): 1393-1402.

    Carter W P L. Documentation of the SAPRC-99chemical mechanism for VOC reactivity assessment [J]. Contract, 2000, 92(329):95-308.

    Jiménez P, Parra R, Baldasano J M. Influence of initial and boundary conditions for ozone modeling in very complex terrains: A case study in the northeastern Iberian Peninsula [J]. Environmental Modelling & Software, 2007,22(9):1294-1306.

    Shen H, Huang Y, Wang R, et al. Global atmospheric emissions of polycyclic aromatic hydrocarbons from 1960 to 2008 and future predictions [J]. Environmental Science & Technology, 2013, 47(12):6415-6424.

    Halsall C J. Investigating the occurrence of persistent organic pollutants (POPs) in the arctic: Their atmospheric behaviour and interaction with the seasonal snow pack [J]. Environmental Pollution, 2004,128(1):163-175.

    Baek S O, Goldstone M E, Kirk P W W, et al. Phase distribution and particle size dependency of polycyclic aromatic hydrocarbons in the urban atmosphere [J]. Chemosphere, 1991,22(5):503-520.

    畢新慧,盛國英,譚吉華,等.多環(huán)芳烴(PAHs)在大氣中的相分布 [J]. 環(huán)境科學學報, 2004,24(1):101-106.

    Yamasaki H, Kuwata K, Miyamoto H. Effects of ambient temperature on aspects of airborne polycyclic aromatic hydrocarbons [J]. Environmental Science & Technology, 1982, 16(4):189-194.

    Odum H T. Ecological and general systems: an introduction to systems ecology [M]. University Press of Colorado, 1994.

    Pankow J F. An absorption model of gas/particle partitioning of organic compounds in the atmosphere [J]. Atmospheric Environment, 1994a,28(2):185-188.

    Junge C E. Basic considerations about trace constituents in the atmosphere as related to the fate of global pollutants [J]. Advances in Environmental Science and Technology 8 (Pt. 1), 1977:7-25.

    Pankow J F. Review and comparative analysis of the theories on partitioning between the gas and aerosol particulate phases in the atmosphere [J]. Atmospheric Environment (1967), 1987,21(11): 2275-2283.

    Finizio A, Mackay D, Bidleman T, et al. Octanol-air partition coefficient as a predictor of partitioning of semi-volatile organic chemicals to aerosols [J]. Atmospheric Environment, 1997,31(15): 2289-2296.

    Harner T, Bidleman T F. Measurements of octanol-air partition coefficients for polychlorinated biphenyls [J]. Journal of Chemical & Engineering Data, 1996,41(4):895-899.

    Harner T, Bidleman T F. Octanol-air partition coefficient for describing particle/gas partitioning of aromatic compounds in urban air [J]. Environmental Science & Technology, 1998,32(10): 1494-1502.

    Odabasi M, Cetin E, Sofuoglu A. Determination of octanol-air partition coefficients and supercooled liquid vapor pressures of PAHs as a function of temperature: application to gas-particle partitioning in an urban atmosphere [J]. Atmospheric Environment, 2006,40(34):6615-6625.

    The Staffs of the California Air Resources Board and the Office of Environmental Health Hazard-Assessment. Benzo[a]pyrene as a toxic air contaminant [R]. The California Air Resources Board and the Environmental Health Hazard-Assessment, 1994,(4): 39650-39662.

    Elena Mantseva, Alexander Malanichev, Nadezhda Vulykh. Poly aromatic hydrocarbons in the environment [R]. MSC-E Technical Note 9/2002.

    Keyte I J, Harrison R M, Lammel G. Chemical reactivity and long-range transport potential of polycyclic aromatic hydrocarbons-a review [J]. Chemical Society Reviews, 2013, 42(24):9333-9391.

    Bewersdorff I, Aulinger A, Matthias V, et al. The effect of temporal resolution of PAH emission data on transport and deposition patterns simulated with the Community Multiscale Air Quality modelling system (CMAQ) [J]. Meteorologische Zeitschrift, 2009,18(1):41-53.

    Berdowski J J M, Baas J, Bloos J P J, et al. The European emission inventory of heavy metals and persistent organic pollutants. Umweltforschungsplan des Bundesministers fur Umwelt, naturschutz und Reaktorsicherheit. Forschungsbericht 104 02 672/03, 1997.

    Aulinger A, Matthias V, Quante M. Introducing a partitioning mechanism for PAHs into the community multiscale air quality modeling system and its application to simulating the transport of benzo (a) pyrene over Europe [J]. Journal of Applied Meteorology and Climatology, 2007,46(11):1718-1730.

    * 責任作者, 研究員, mengfan@craes.org.cn

    Simulations of the atmospheric benzo[a]pyrene in China using CMAQ model

    YAN Wen-jun1,2, LIU Min1, LIU Shi-jie2, HE You-jiang2, TANG Wei2, WANG Qing1, LI Ye1, MENG Fan2*

    (1.Key Laboratory of Geographic Information Science, Ministry of Education, College of Geographic Sciences, East China Normal University, Shanghai 200241, China;2.Atmospheric Environment Institute, Chinese Research Academy of Environmental Sciences, Beijing 100012, China)., 2016,36(6):1681~1689

    The modules for simulating transport, transformation and deposition processes, as well as gas/particle partitioning (OAabsorption mechanism) processes of atmospheric BaP was implemented into the CMAQ, a regional air quality model developed by U.S. E.P.A. The spatial distributions of BaP concentrations, dry and wet deposition rates are simulated for January, April, July and October 2014 in China. The comparison between modeled and observed BaP concentrations shows the predicted values were in an acceptable range. The predicted deposition results show that the wet deposition is the dominant deposition pathway. The spatial distributions of the BaP depositions are closely related to that of the anthropogenic emissions and present obvious long-distance transportation characteristics.

    CMAQ;BaP;seasonal variation;dry deposition;wet deposition

    X512

    A

    1000-6923(2016)06-1681-09

    閆文君(1991-),女,甘肅敦煌人,華東師范大學地理科學學院碩士研究生,主要從事大氣環(huán)境模擬研究工作

    2015-11-06

    國家自然科學基金重點項目(41130525);國家自然科學基金面上項目(21177118,41371451);國家自然科學基金項目(41501530)

    猜你喜歡
    氣相大氣顆粒
    大氣的呵護
    軍事文摘(2023年10期)2023-06-09 09:15:06
    Efficacy and safety of Mianyi granules (免疫Ⅱ顆粒) for reversal of immune nonresponse following antiretroviral therapy of human immunodeficiency virus-1:a randomized,double-blind,multi-center,placebo-controlled trial
    氣相過渡金屬鈦-碳鏈團簇的研究
    要讓顆粒都歸倉
    心聲歌刊(2019年1期)2019-05-09 03:21:32
    疏風定喘顆粒輔料的篩選
    中成藥(2017年4期)2017-05-17 06:09:29
    大氣古樸揮灑自如
    大氣、水之后,土十條來了
    新型釩基催化劑催化降解氣相二噁英
    預(yù)縮聚反應(yīng)器氣相管“鼓泡”的成因探討
    氣相防銹技術(shù)在電器設(shè)備防腐中的應(yīng)用
    人人澡人人妻人| 久久香蕉激情| 午夜免费激情av| 日韩高清综合在线| 黄色视频,在线免费观看| 国产精品99久久99久久久不卡| 久久精品国产综合久久久| 亚洲欧美一区二区三区久久| 成人永久免费在线观看视频| 久久青草综合色| av国产精品久久久久影院| 在线观看舔阴道视频| 99国产精品一区二区三区| 少妇粗大呻吟视频| 久久久久久久午夜电影 | 亚洲国产欧美网| 999久久久精品免费观看国产| 久久精品aⅴ一区二区三区四区| 一个人免费在线观看的高清视频| 在线观看免费午夜福利视频| netflix在线观看网站| 欧美一区二区精品小视频在线| 免费在线观看完整版高清| 老鸭窝网址在线观看| 香蕉国产在线看| 国产aⅴ精品一区二区三区波| 成年版毛片免费区| 宅男免费午夜| 久久午夜亚洲精品久久| 伦理电影免费视频| 手机成人av网站| 一级a爱片免费观看的视频| 99精国产麻豆久久婷婷| 中出人妻视频一区二区| 9热在线视频观看99| 免费少妇av软件| 国产精品av久久久久免费| 国产成人av教育| 亚洲精品av麻豆狂野| 老汉色∧v一级毛片| 日韩欧美三级三区| 69精品国产乱码久久久| 脱女人内裤的视频| 身体一侧抽搐| 国产午夜精品久久久久久| 亚洲色图 男人天堂 中文字幕| 最近最新中文字幕大全电影3 | 欧美中文综合在线视频| 国产亚洲av高清不卡| 国产主播在线观看一区二区| 丝袜人妻中文字幕| 黑人欧美特级aaaaaa片| 久热爱精品视频在线9| avwww免费| а√天堂www在线а√下载| 欧美黄色淫秽网站| 9热在线视频观看99| 一区福利在线观看| 亚洲黑人精品在线| 老司机福利观看| 久久精品亚洲精品国产色婷小说| 日韩精品免费视频一区二区三区| 交换朋友夫妻互换小说| 亚洲国产中文字幕在线视频| 亚洲五月色婷婷综合| 亚洲精品久久成人aⅴ小说| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕色久视频| 老司机亚洲免费影院| 成人三级黄色视频| 老鸭窝网址在线观看| 丰满人妻熟妇乱又伦精品不卡| 精品少妇一区二区三区视频日本电影| 日韩三级视频一区二区三区| av片东京热男人的天堂| av有码第一页| 亚洲aⅴ乱码一区二区在线播放 | 国产一区二区在线av高清观看| 波多野结衣一区麻豆| 9色porny在线观看| 日本免费一区二区三区高清不卡 | 嫁个100分男人电影在线观看| 久久久久精品国产欧美久久久| 一二三四在线观看免费中文在| 两个人免费观看高清视频| 18禁裸乳无遮挡免费网站照片 | 久久人妻福利社区极品人妻图片| 日本a在线网址| 90打野战视频偷拍视频| 免费在线观看完整版高清| 黄频高清免费视频| 又大又爽又粗| 中文字幕人妻熟女乱码| 首页视频小说图片口味搜索| 亚洲熟妇熟女久久| 90打野战视频偷拍视频| 一级片免费观看大全| 精品熟女少妇八av免费久了| 久久精品国产99精品国产亚洲性色 | 老司机亚洲免费影院| 悠悠久久av| 精品一区二区三卡| 欧美中文综合在线视频| 在线观看一区二区三区激情| 亚洲片人在线观看| 深夜精品福利| 久久久久久久精品吃奶| 三级毛片av免费| 人妻丰满熟妇av一区二区三区| 激情在线观看视频在线高清| 男女下面进入的视频免费午夜 | 黄色女人牲交| 黄色视频不卡| 男人的好看免费观看在线视频 | 日韩人妻精品一区2区三区| 国产av精品麻豆| 日韩高清综合在线| 在线观看免费日韩欧美大片| 90打野战视频偷拍视频| 少妇裸体淫交视频免费看高清 | 亚洲精华国产精华精| 长腿黑丝高跟| 国产精品爽爽va在线观看网站 | 老司机午夜福利在线观看视频| 视频在线观看一区二区三区| 国产单亲对白刺激| 女人精品久久久久毛片| 成熟少妇高潮喷水视频| 99精品在免费线老司机午夜| 成年女人毛片免费观看观看9| 国产99白浆流出| 免费高清在线观看日韩| 夜夜爽天天搞| 欧美日韩福利视频一区二区| x7x7x7水蜜桃| 91麻豆av在线| 自拍欧美九色日韩亚洲蝌蚪91| 九色亚洲精品在线播放| 免费日韩欧美在线观看| 在线观看66精品国产| 亚洲全国av大片| ponron亚洲| 精品国产超薄肉色丝袜足j| 男女做爰动态图高潮gif福利片 | 女人被躁到高潮嗷嗷叫费观| 岛国视频午夜一区免费看| 免费人成视频x8x8入口观看| 嫁个100分男人电影在线观看| 女人被狂操c到高潮| av电影中文网址| 一区二区三区激情视频| 久久狼人影院| а√天堂www在线а√下载| 一区二区三区精品91| 久久这里只有精品19| 天天躁夜夜躁狠狠躁躁| 国产精品乱码一区二三区的特点 | 亚洲av熟女| 亚洲aⅴ乱码一区二区在线播放 | 国产黄色免费在线视频| 激情在线观看视频在线高清| 国产精品国产高清国产av| 九色亚洲精品在线播放| 欧美成人午夜精品| 好看av亚洲va欧美ⅴa在| 国产亚洲精品一区二区www| 天堂俺去俺来也www色官网| 国产欧美日韩一区二区三区在线| 99re在线观看精品视频| www.www免费av| 97超级碰碰碰精品色视频在线观看| 操出白浆在线播放| 精品午夜福利视频在线观看一区| 最近最新中文字幕大全电影3 | 国产成人欧美在线观看| 最好的美女福利视频网| 亚洲专区中文字幕在线| 神马国产精品三级电影在线观看 | 国产免费现黄频在线看| 国产av一区二区精品久久| 国产精品九九99| 免费看十八禁软件| 97超级碰碰碰精品色视频在线观看| 亚洲 欧美一区二区三区| 国产无遮挡羞羞视频在线观看| 午夜福利一区二区在线看| 母亲3免费完整高清在线观看| 麻豆av在线久日| 精品一区二区三区四区五区乱码| 久久国产精品人妻蜜桃| 亚洲国产精品合色在线| 人人妻,人人澡人人爽秒播| 日本 av在线| 亚洲av日韩精品久久久久久密| 黄片大片在线免费观看| 色综合欧美亚洲国产小说| 欧美乱色亚洲激情| 女人被躁到高潮嗷嗷叫费观| 欧美老熟妇乱子伦牲交| 男女下面插进去视频免费观看| 狂野欧美激情性xxxx| 亚洲午夜理论影院| 一边摸一边做爽爽视频免费| 久久国产精品男人的天堂亚洲| 国产成年人精品一区二区 | 成人精品一区二区免费| 动漫黄色视频在线观看| 少妇被粗大的猛进出69影院| 欧美日韩乱码在线| 日本一区二区免费在线视频| av国产精品久久久久影院| 天堂俺去俺来也www色官网| 精品一品国产午夜福利视频| 天天添夜夜摸| bbb黄色大片| 国产成人精品久久二区二区免费| 一区二区三区激情视频| 国产成年人精品一区二区 | 丝袜在线中文字幕| 国产亚洲精品第一综合不卡| 亚洲精品国产区一区二| 99国产精品免费福利视频| 亚洲精华国产精华精| 久久精品国产99精品国产亚洲性色 | 老司机靠b影院| 中文字幕人妻丝袜一区二区| 青草久久国产| 91大片在线观看| 999精品在线视频| 亚洲精品美女久久av网站| 久9热在线精品视频| 窝窝影院91人妻| 97碰自拍视频| 大陆偷拍与自拍| 操美女的视频在线观看| av电影中文网址| 免费少妇av软件| 一个人观看的视频www高清免费观看 | 男女下面插进去视频免费观看| xxxhd国产人妻xxx| 国产又色又爽无遮挡免费看| 三级毛片av免费| 夜夜看夜夜爽夜夜摸 | 久久国产亚洲av麻豆专区| 午夜福利一区二区在线看| 老司机午夜十八禁免费视频| 在线观看www视频免费| 一夜夜www| 看片在线看免费视频| 久久天堂一区二区三区四区| 午夜久久久在线观看| 伊人久久大香线蕉亚洲五| 又大又爽又粗| 水蜜桃什么品种好| 国产成人精品无人区| 国产1区2区3区精品| 国产免费av片在线观看野外av| 18禁黄网站禁片午夜丰满| 久久香蕉精品热| 狠狠狠狠99中文字幕| 中文字幕av电影在线播放| 欧美午夜高清在线| 色尼玛亚洲综合影院| 午夜影院日韩av| 久久午夜亚洲精品久久| 亚洲三区欧美一区| 国产黄色免费在线视频| 欧美国产精品va在线观看不卡| 欧美日本中文国产一区发布| 最近最新中文字幕大全电影3 | 一级a爱视频在线免费观看| 久久久精品欧美日韩精品| 精品久久久久久久久久免费视频 | 很黄的视频免费| 亚洲一区二区三区色噜噜 | 不卡av一区二区三区| 在线天堂中文资源库| 成人影院久久| 他把我摸到了高潮在线观看| 久久伊人香网站| 91成人精品电影| 久9热在线精品视频| 国产亚洲av高清不卡| 满18在线观看网站| 国产成+人综合+亚洲专区| 美女高潮到喷水免费观看| 国内久久婷婷六月综合欲色啪| 新久久久久国产一级毛片| 国产欧美日韩一区二区精品| 日本黄色视频三级网站网址| 麻豆av在线久日| 欧美日韩瑟瑟在线播放| 国产有黄有色有爽视频| 亚洲一区二区三区不卡视频| 国产成人精品在线电影| 中文字幕高清在线视频| 亚洲精品美女久久av网站| 久久影院123| 午夜福利欧美成人| av网站在线播放免费| 午夜影院日韩av| 女警被强在线播放| 欧美日韩亚洲高清精品| 日韩大码丰满熟妇| 老司机在亚洲福利影院| 大型av网站在线播放| 一本综合久久免费| 成在线人永久免费视频| 亚洲色图综合在线观看| 十八禁网站免费在线| 国产精品永久免费网站| 美女午夜性视频免费| 亚洲免费av在线视频| 好男人电影高清在线观看| 欧美色视频一区免费| 一级黄色大片毛片| 十八禁人妻一区二区| 九色亚洲精品在线播放| 国产精品久久视频播放| 女同久久另类99精品国产91| 亚洲国产欧美一区二区综合| 亚洲色图av天堂| 欧美成人性av电影在线观看| 久久精品亚洲熟妇少妇任你| 亚洲人成77777在线视频| 日韩大尺度精品在线看网址 | 一区二区三区精品91| 一级a爱视频在线免费观看| 国产成人精品久久二区二区91| 久久精品国产综合久久久| 19禁男女啪啪无遮挡网站| 国产精品一区二区免费欧美| 在线观看免费视频网站a站| 大码成人一级视频| 一级a爱片免费观看的视频| 久久久久久亚洲精品国产蜜桃av| 亚洲成人精品中文字幕电影 | 欧美一区二区精品小视频在线| 香蕉国产在线看| 夜夜看夜夜爽夜夜摸 | 亚洲国产精品合色在线| 交换朋友夫妻互换小说| 一二三四在线观看免费中文在| 欧美在线一区亚洲| 亚洲午夜精品一区,二区,三区| 久久精品国产亚洲av高清一级| 中文字幕精品免费在线观看视频| 中文字幕另类日韩欧美亚洲嫩草| 久久精品亚洲熟妇少妇任你| 亚洲aⅴ乱码一区二区在线播放 | bbb黄色大片| 国产精品综合久久久久久久免费 | 成年女人毛片免费观看观看9| 色在线成人网| 欧美成人免费av一区二区三区| 搡老岳熟女国产| 欧美丝袜亚洲另类 | 黄色视频,在线免费观看| 色播在线永久视频| 丝袜美腿诱惑在线| 免费日韩欧美在线观看| 老鸭窝网址在线观看| 欧美日韩av久久| 一个人观看的视频www高清免费观看 | 亚洲专区字幕在线| 日本免费a在线| 精品久久久精品久久久| 日本a在线网址| 久久精品亚洲熟妇少妇任你| 免费搜索国产男女视频| 国产精品亚洲一级av第二区| 亚洲五月婷婷丁香| 成人av一区二区三区在线看| 久99久视频精品免费| 少妇粗大呻吟视频| 国产麻豆69| 日本五十路高清| 亚洲成a人片在线一区二区| 国产三级在线视频| 狂野欧美激情性xxxx| 亚洲国产精品999在线| 欧美丝袜亚洲另类 | 人成视频在线观看免费观看| 又紧又爽又黄一区二区| 大香蕉久久成人网| 999精品在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 97碰自拍视频| 国产单亲对白刺激| 久久亚洲精品不卡| 久久人妻av系列| 日本精品一区二区三区蜜桃| 看片在线看免费视频| 国产成人精品在线电影| 超碰成人久久| 亚洲黑人精品在线| 夜夜看夜夜爽夜夜摸 | 在线观看66精品国产| 久久久久国内视频| 精品国产超薄肉色丝袜足j| 真人一进一出gif抽搐免费| 国产成人欧美| 性少妇av在线| 国产成人系列免费观看| 免费久久久久久久精品成人欧美视频| 国产蜜桃级精品一区二区三区| 国内久久婷婷六月综合欲色啪| 国产91精品成人一区二区三区| 黄片播放在线免费| 国产午夜精品久久久久久| 色综合婷婷激情| 黑丝袜美女国产一区| 高清欧美精品videossex| 男女之事视频高清在线观看| 手机成人av网站| 丝袜人妻中文字幕| 久久久久久人人人人人| 老司机午夜十八禁免费视频| 亚洲国产看品久久| 国产日韩一区二区三区精品不卡| 国产一区二区激情短视频| 一区二区三区国产精品乱码| 亚洲人成77777在线视频| 亚洲自偷自拍图片 自拍| 一夜夜www| 成年版毛片免费区| 久久精品亚洲av国产电影网| 91字幕亚洲| av片东京热男人的天堂| 啦啦啦免费观看视频1| 亚洲第一欧美日韩一区二区三区| 国产高清激情床上av| 精品第一国产精品| 久久精品影院6| 美女福利国产在线| 午夜亚洲福利在线播放| 在线观看舔阴道视频| 天堂中文最新版在线下载| 美女国产高潮福利片在线看| 在线观看午夜福利视频| 纯流量卡能插随身wifi吗| 国产精品秋霞免费鲁丝片| 最好的美女福利视频网| 日本三级黄在线观看| a级片在线免费高清观看视频| 50天的宝宝边吃奶边哭怎么回事| 视频区欧美日本亚洲| 高清欧美精品videossex| av超薄肉色丝袜交足视频| 国产99白浆流出| 日本wwww免费看| 亚洲成人精品中文字幕电影 | 村上凉子中文字幕在线| 国内久久婷婷六月综合欲色啪| 国产精品一区二区三区四区久久 | 真人一进一出gif抽搐免费| 三上悠亚av全集在线观看| 老司机午夜十八禁免费视频| 日本 av在线| 人人妻人人添人人爽欧美一区卜| 黄色视频,在线免费观看| a级毛片黄视频| 日韩大码丰满熟妇| 成人国产一区最新在线观看| 一级片'在线观看视频| 在线免费观看的www视频| 美女福利国产在线| 亚洲一码二码三码区别大吗| 久久香蕉精品热| 视频区图区小说| 巨乳人妻的诱惑在线观看| 欧美av亚洲av综合av国产av| 国产aⅴ精品一区二区三区波| 国产一区二区在线av高清观看| 欧美乱色亚洲激情| a在线观看视频网站| 男人舔女人的私密视频| 亚洲第一欧美日韩一区二区三区| 在线免费观看的www视频| 看黄色毛片网站| 99国产精品一区二区蜜桃av| 亚洲精品中文字幕在线视频| 欧美黄色淫秽网站| 丰满的人妻完整版| 久久久久久久久中文| а√天堂www在线а√下载| 亚洲三区欧美一区| av片东京热男人的天堂| 他把我摸到了高潮在线观看| 看黄色毛片网站| 交换朋友夫妻互换小说| 12—13女人毛片做爰片一| 久久久国产成人精品二区 | 黄色毛片三级朝国网站| 男女之事视频高清在线观看| 亚洲成人免费电影在线观看| 99久久国产精品久久久| 久久精品aⅴ一区二区三区四区| 久久青草综合色| 久久人妻熟女aⅴ| 黑人巨大精品欧美一区二区mp4| 视频区图区小说| av欧美777| av国产精品久久久久影院| 一夜夜www| 久久久久久免费高清国产稀缺| 国产成人一区二区三区免费视频网站| 99re在线观看精品视频| 久久青草综合色| 精品免费久久久久久久清纯| 不卡一级毛片| 成年女人毛片免费观看观看9| 在线观看一区二区三区激情| 老司机福利观看| 国产免费av片在线观看野外av| 88av欧美| 久久午夜综合久久蜜桃| 午夜两性在线视频| 久久久久久久久久久久大奶| 亚洲人成网站在线播放欧美日韩| 香蕉国产在线看| 老熟妇仑乱视频hdxx| 两人在一起打扑克的视频| 亚洲欧美日韩无卡精品| 久久久久国产一级毛片高清牌| avwww免费| √禁漫天堂资源中文www| 精品国产美女av久久久久小说| 丰满迷人的少妇在线观看| 黄片小视频在线播放| 丁香六月欧美| 大陆偷拍与自拍| 曰老女人黄片| 亚洲人成电影观看| 99久久久亚洲精品蜜臀av| 久久久精品国产亚洲av高清涩受| 黄色 视频免费看| x7x7x7水蜜桃| 久久狼人影院| 久久热在线av| 久久久久久大精品| 热99国产精品久久久久久7| 色精品久久人妻99蜜桃| 不卡一级毛片| 99久久精品国产亚洲精品| 巨乳人妻的诱惑在线观看| 成人亚洲精品av一区二区 | 在线观看免费日韩欧美大片| 国产欧美日韩一区二区三区在线| 欧美日韩视频精品一区| 国产欧美日韩一区二区三区在线| 日韩高清综合在线| 超碰97精品在线观看| 视频区图区小说| 国产在线观看jvid| 伦理电影免费视频| 成年人免费黄色播放视频| 看免费av毛片| 国产亚洲欧美在线一区二区| 国产熟女xx| 国产亚洲精品久久久久5区| 欧美激情久久久久久爽电影 | 国产精品久久久久久人妻精品电影| 欧美精品亚洲一区二区| 午夜免费激情av| 每晚都被弄得嗷嗷叫到高潮| 极品人妻少妇av视频| 亚洲精品国产色婷婷电影| 国产不卡一卡二| 人人妻人人添人人爽欧美一区卜| 国产一区二区三区在线臀色熟女 | 国产一卡二卡三卡精品| 国产成人av激情在线播放| 久久久久久大精品| 美女国产高潮福利片在线看| 男人的好看免费观看在线视频 | av有码第一页| 日韩国内少妇激情av| 日韩大码丰满熟妇| 女人被躁到高潮嗷嗷叫费观| 性欧美人与动物交配| 久久精品国产清高在天天线| 最新美女视频免费是黄的| 久久精品影院6| 亚洲av片天天在线观看| 中文字幕另类日韩欧美亚洲嫩草| 又黄又粗又硬又大视频| 大型av网站在线播放| 亚洲精品成人av观看孕妇| 超色免费av| 久久婷婷成人综合色麻豆| 热99re8久久精品国产| 亚洲欧美日韩无卡精品| 法律面前人人平等表现在哪些方面| 91麻豆精品激情在线观看国产 | 日本黄色日本黄色录像| 久久久久久久久中文| 一夜夜www| 麻豆成人av在线观看| 久久香蕉激情| 新久久久久国产一级毛片| 欧美成人免费av一区二区三区| 在线永久观看黄色视频| 一级毛片女人18水好多| 午夜福利免费观看在线| 日韩国内少妇激情av| 这个男人来自地球电影免费观看| 黑人巨大精品欧美一区二区mp4| 亚洲欧美激情综合另类| 99精品欧美一区二区三区四区| a级毛片黄视频| 免费在线观看亚洲国产| 99精品欧美一区二区三区四区| 欧美激情久久久久久爽电影 | 欧美乱码精品一区二区三区|