• <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)用
    久久6这里有精品| 日日啪夜夜爽| 欧美人与善性xxx| 亚洲激情五月婷婷啪啪| 国产成人福利小说| 亚洲精品日本国产第一区| 亚洲精品一区蜜桃| 大香蕉久久网| 国产一区二区在线观看日韩| 午夜亚洲福利在线播放| 久久精品久久久久久久性| 80岁老熟妇乱子伦牲交| 91久久精品电影网| av又黄又爽大尺度在线免费看| 18禁在线无遮挡免费观看视频| 男人狂女人下面高潮的视频| 伊人久久精品亚洲午夜| 亚洲成人中文字幕在线播放| 视频中文字幕在线观看| 99久国产av精品| 日韩欧美一区视频在线观看 | av天堂中文字幕网| 亚洲欧美精品自产自拍| 我的老师免费观看完整版| 欧美日韩综合久久久久久| 欧美不卡视频在线免费观看| 精品国产露脸久久av麻豆 | 国内揄拍国产精品人妻在线| 亚洲国产精品sss在线观看| 精品久久久久久电影网| 一区二区三区乱码不卡18| 亚洲精品乱码久久久v下载方式| 国产一区二区三区综合在线观看 | 国产精品一区二区三区四区免费观看| 看免费成人av毛片| 三级经典国产精品| 最近中文字幕2019免费版| 在线天堂最新版资源| 国内精品美女久久久久久| 久久久亚洲精品成人影院| 国产精品福利在线免费观看| 日韩伦理黄色片| 在线免费观看的www视频| 日韩强制内射视频| 一区二区三区四区激情视频| 美女内射精品一级片tv| 简卡轻食公司| 亚洲成人一二三区av| 国产 一区 欧美 日韩| 国产亚洲av片在线观看秒播厂 | 天堂中文最新版在线下载 | 免费看不卡的av| 亚洲成色77777| 最近中文字幕2019免费版| 国内少妇人妻偷人精品xxx网站| 两个人的视频大全免费| 国产 一区精品| av免费观看日本| 男女下面进入的视频免费午夜| 国产精品人妻久久久影院| 国产国拍精品亚洲av在线观看| 亚洲av成人av| 听说在线观看完整版免费高清| 六月丁香七月| 亚洲欧美清纯卡通| 男人爽女人下面视频在线观看| 午夜精品国产一区二区电影 | h日本视频在线播放| 纵有疾风起免费观看全集完整版 | 最近视频中文字幕2019在线8| 99久久中文字幕三级久久日本| 亚洲国产欧美人成| 女人十人毛片免费观看3o分钟| 国产成人aa在线观看| 如何舔出高潮| 国产色婷婷99| 一个人看视频在线观看www免费| 国产白丝娇喘喷水9色精品| 一区二区三区高清视频在线| 亚洲精品日韩av片在线观看| 草草在线视频免费看| 午夜激情久久久久久久| 成人漫画全彩无遮挡| 国语对白做爰xxxⅹ性视频网站| 三级国产精品片| 亚洲国产成人一精品久久久| 日韩欧美国产在线观看| 欧美+日韩+精品| 亚洲色图av天堂| 国产av在哪里看| av又黄又爽大尺度在线免费看| 一级a做视频免费观看| 免费观看性生交大片5| av国产久精品久网站免费入址| 好男人在线观看高清免费视频| 一级黄片播放器| a级一级毛片免费在线观看| 精品不卡国产一区二区三区| 人妻夜夜爽99麻豆av| 亚洲精品日韩在线中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 欧美精品国产亚洲| 国产av国产精品国产| 免费黄网站久久成人精品| 三级男女做爰猛烈吃奶摸视频| 免费av毛片视频| videossex国产| 中文在线观看免费www的网站| 肉色欧美久久久久久久蜜桃 | 女人十人毛片免费观看3o分钟| 亚洲久久久久久中文字幕| 波野结衣二区三区在线| 国产亚洲午夜精品一区二区久久 | 中文字幕久久专区| 日本与韩国留学比较| 少妇裸体淫交视频免费看高清| 女的被弄到高潮叫床怎么办| 精品国产露脸久久av麻豆 | av黄色大香蕉| 亚洲av日韩在线播放| 日韩不卡一区二区三区视频在线| 国产成人福利小说| 欧美性猛交╳xxx乱大交人| 免费观看性生交大片5| 国产成人a∨麻豆精品| av女优亚洲男人天堂| 啦啦啦啦在线视频资源| 成年女人在线观看亚洲视频 | 亚洲av男天堂| 97精品久久久久久久久久精品| 亚洲,欧美,日韩| 国内揄拍国产精品人妻在线| 国产精品一区www在线观看| 97超碰精品成人国产| 日本色播在线视频| 毛片一级片免费看久久久久| 美女被艹到高潮喷水动态| 久久99蜜桃精品久久| 91久久精品国产一区二区成人| 国产大屁股一区二区在线视频| 国产精品一区二区性色av| 淫秽高清视频在线观看| 亚洲av在线观看美女高潮| 国产免费福利视频在线观看| 搡老妇女老女人老熟妇| 日韩视频在线欧美| 街头女战士在线观看网站| 欧美zozozo另类| 国产 一区 欧美 日韩| 日韩人妻高清精品专区| 午夜福利视频精品| 天堂av国产一区二区熟女人妻| 国产三级在线视频| 日本一二三区视频观看| 欧美激情国产日韩精品一区| 国产亚洲91精品色在线| 国产精品蜜桃在线观看| 午夜福利高清视频| 午夜福利视频精品| 免费av观看视频| 最新中文字幕久久久久| 久久久久精品性色| h日本视频在线播放| 午夜激情久久久久久久| 内地一区二区视频在线| freevideosex欧美| 亚洲国产最新在线播放| 尾随美女入室| 亚洲av.av天堂| 亚洲成人久久爱视频| 伊人久久精品亚洲午夜| 麻豆国产97在线/欧美| 欧美3d第一页| 色网站视频免费| 有码 亚洲区| 国产精品蜜桃在线观看| 日韩av不卡免费在线播放| 搡老妇女老女人老熟妇| 国产激情偷乱视频一区二区| 七月丁香在线播放| 成人二区视频| 欧美精品一区二区大全| 一个人观看的视频www高清免费观看| 99热这里只有是精品50| 精华霜和精华液先用哪个| 午夜免费激情av| 亚洲欧美成人精品一区二区| 国产欧美另类精品又又久久亚洲欧美| 特大巨黑吊av在线直播| 男人舔奶头视频| 91精品伊人久久大香线蕉| 中文字幕久久专区| 精品久久久久久电影网| 国产一级毛片在线| 国产一级毛片七仙女欲春2| 美女xxoo啪啪120秒动态图| 久久精品国产鲁丝片午夜精品| 亚洲精品视频女| 免费观看性生交大片5| 91精品国产九色| 人妻制服诱惑在线中文字幕| 99热网站在线观看| 超碰97精品在线观看| 欧美三级亚洲精品| 少妇猛男粗大的猛烈进出视频 | 久久精品人妻少妇| 一个人观看的视频www高清免费观看| 国产淫语在线视频| 精品久久久久久久久亚洲| 好男人在线观看高清免费视频| 国产成人freesex在线| 久久久久久久久久黄片| 少妇人妻精品综合一区二区| 中文资源天堂在线| 七月丁香在线播放| 免费黄色在线免费观看| 欧美潮喷喷水| 亚州av有码| 男人狂女人下面高潮的视频| 久久精品人妻少妇| 在线 av 中文字幕| 高清午夜精品一区二区三区| 91在线精品国自产拍蜜月| 精品一区二区三卡| 免费av毛片视频| 少妇熟女欧美另类| 精品国产露脸久久av麻豆 | 真实男女啪啪啪动态图| 夜夜看夜夜爽夜夜摸| 精品久久国产蜜桃| 天堂俺去俺来也www色官网 | 日韩大片免费观看网站| 最后的刺客免费高清国语| videos熟女内射| 日韩一本色道免费dvd| 激情 狠狠 欧美| 天堂俺去俺来也www色官网 | 日韩强制内射视频| 人妻少妇偷人精品九色| 欧美97在线视频| 精品一区二区三区人妻视频| 国产成人免费观看mmmm| 欧美xxxx黑人xx丫x性爽| 免费黄网站久久成人精品| 女人久久www免费人成看片| 国产熟女欧美一区二区| 欧美97在线视频| 日韩精品青青久久久久久| 欧美成人一区二区免费高清观看| 精品少妇黑人巨大在线播放| 中国国产av一级| 日本一本二区三区精品| 3wmmmm亚洲av在线观看| 欧美3d第一页| 国产亚洲5aaaaa淫片| 免费观看的影片在线观看| 天天一区二区日本电影三级| 色哟哟·www| 高清毛片免费看| 久久99蜜桃精品久久| 国产淫片久久久久久久久| .国产精品久久| 欧美日韩视频高清一区二区三区二| 中文字幕久久专区| 久久久久久久久久黄片| 天堂影院成人在线观看| 亚洲怡红院男人天堂| 99热这里只有是精品在线观看| 肉色欧美久久久久久久蜜桃 | 午夜福利视频精品| 久久久久免费精品人妻一区二区| 国产亚洲91精品色在线| av线在线观看网站| 少妇被粗大猛烈的视频| 干丝袜人妻中文字幕| 边亲边吃奶的免费视频| 老司机影院成人| 国产白丝娇喘喷水9色精品| 深爱激情五月婷婷| 久久这里有精品视频免费| 国产 亚洲一区二区三区 | 干丝袜人妻中文字幕| 22中文网久久字幕| 男人和女人高潮做爰伦理| 欧美不卡视频在线免费观看| 亚洲欧美日韩卡通动漫| 成人高潮视频无遮挡免费网站| 国产精品日韩av在线免费观看| 十八禁国产超污无遮挡网站| 国产成人91sexporn| 欧美日韩亚洲高清精品| 国产一区二区三区av在线| 亚洲人与动物交配视频| 久久热精品热| 久久久久网色| 国产精品一区二区性色av| 欧美日韩国产mv在线观看视频 | 97超碰精品成人国产| 一级片'在线观看视频| 久久久久久久久久久丰满| 精品人妻视频免费看| 日韩 亚洲 欧美在线| 亚洲欧美精品自产自拍| 国产亚洲精品久久久com| 国产白丝娇喘喷水9色精品| 视频中文字幕在线观看| 国产精品99久久久久久久久| 久久综合国产亚洲精品| 伦精品一区二区三区| 91久久精品国产一区二区三区| 欧美一级a爱片免费观看看| 99热网站在线观看| 老师上课跳d突然被开到最大视频| 在线观看一区二区三区| 亚洲国产成人一精品久久久| 国产不卡一卡二| 干丝袜人妻中文字幕| 神马国产精品三级电影在线观看| 国产黄色视频一区二区在线观看| 国产女主播在线喷水免费视频网站 | 99热全是精品| 日韩av不卡免费在线播放| 久久久久网色| av网站免费在线观看视频 | 国产精品一区二区在线观看99 | kizo精华| 淫秽高清视频在线观看| 亚洲欧美成人精品一区二区| 99九九线精品视频在线观看视频| 人妻一区二区av| 激情 狠狠 欧美| 一级毛片aaaaaa免费看小| av又黄又爽大尺度在线免费看| 天堂影院成人在线观看| 日韩大片免费观看网站| 国产成人a区在线观看| 久久久久久伊人网av| 99久久人妻综合| 黄片wwwwww| 国产精品久久久久久久久免| 黄片wwwwww| 晚上一个人看的免费电影| 欧美日韩在线观看h| 99热这里只有是精品在线观看| 午夜日本视频在线| 久久精品国产亚洲网站| 国产av在哪里看| 久久精品国产亚洲网站| 中文欧美无线码| 麻豆乱淫一区二区| 国产三级在线视频| 国产精品久久久久久久久免| 男女下面进入的视频免费午夜| 特大巨黑吊av在线直播| 亚洲美女视频黄频| av黄色大香蕉| 亚洲精品色激情综合| 乱码一卡2卡4卡精品| 日日撸夜夜添| 一级爰片在线观看| 又粗又硬又长又爽又黄的视频| 国产免费又黄又爽又色| 国产成人福利小说| 亚洲精华国产精华液的使用体验| 成人漫画全彩无遮挡| 国产精品不卡视频一区二区| 校园人妻丝袜中文字幕| 国产成人freesex在线| 69人妻影院| 欧美bdsm另类| 中文字幕亚洲精品专区| 国内揄拍国产精品人妻在线| av黄色大香蕉| 亚洲av福利一区| 联通29元200g的流量卡| 高清av免费在线| 中文字幕久久专区| 精品一区在线观看国产| 国产伦理片在线播放av一区| 综合色丁香网| 伦理电影大哥的女人| 欧美一区二区亚洲| 亚洲最大成人中文| 美女内射精品一级片tv| 天堂中文最新版在线下载 | 美女高潮的动态| 亚洲欧美日韩无卡精品| 内地一区二区视频在线| 亚洲一级一片aⅴ在线观看| 亚洲av成人精品一区久久| 亚洲av福利一区| 成人一区二区视频在线观看| 日韩强制内射视频| 日本熟妇午夜| videossex国产| 夜夜看夜夜爽夜夜摸| 亚洲av免费高清在线观看| 两个人的视频大全免费| 韩国av在线不卡| 黄片无遮挡物在线观看| 97精品久久久久久久久久精品| 2022亚洲国产成人精品| 晚上一个人看的免费电影| 久久久久精品久久久久真实原创| 老师上课跳d突然被开到最大视频| 丰满人妻一区二区三区视频av| 久久久a久久爽久久v久久| 国产精品一区二区三区四区久久| 天堂av国产一区二区熟女人妻| 日韩强制内射视频| 人妻制服诱惑在线中文字幕| 丰满少妇做爰视频| 寂寞人妻少妇视频99o| 国内精品宾馆在线| 午夜日本视频在线| 晚上一个人看的免费电影| 国产精品av视频在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美三级三区| 赤兔流量卡办理| 免费看美女性在线毛片视频| 黄色一级大片看看| 国产在视频线精品| 日韩人妻高清精品专区| 韩国高清视频一区二区三区| 亚洲欧洲国产日韩| 一级爰片在线观看| 久久久精品欧美日韩精品| 99re6热这里在线精品视频| 亚洲国产精品成人综合色| 人体艺术视频欧美日本| 久久久久网色| av播播在线观看一区| 欧美极品一区二区三区四区| 在线免费观看的www视频| 美女被艹到高潮喷水动态| 青春草亚洲视频在线观看| 熟妇人妻久久中文字幕3abv| 欧美成人a在线观看| 成年女人看的毛片在线观看| 一级黄片播放器| 亚洲成人中文字幕在线播放| 精品人妻一区二区三区麻豆| 国产精品三级大全| 亚洲在久久综合| 亚洲婷婷狠狠爱综合网| 国产高潮美女av| 国产视频内射| 精品国产三级普通话版| 免费少妇av软件| 日韩强制内射视频| 亚洲色图av天堂| 毛片女人毛片| 免费看a级黄色片| 成年av动漫网址| 亚洲自偷自拍三级| 看免费成人av毛片| 午夜福利视频精品| 成人毛片60女人毛片免费| 99热全是精品| 国产精品av视频在线免费观看| 久久久久久久亚洲中文字幕| 男女国产视频网站| 人妻制服诱惑在线中文字幕| 国产亚洲av嫩草精品影院| 男女边摸边吃奶| 国产一级毛片在线| 在线免费观看的www视频| 精品久久久久久久人妻蜜臀av| 日韩成人av中文字幕在线观看| 免费av不卡在线播放| 成年人午夜在线观看视频 | 国产精品一区www在线观看| 黄片无遮挡物在线观看| 午夜久久久久精精品| 亚洲国产欧美人成| 亚洲欧美日韩东京热| 国语对白做爰xxxⅹ性视频网站| 午夜福利高清视频| 午夜福利成人在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 尾随美女入室| av卡一久久| 寂寞人妻少妇视频99o| 三级国产精品片| 久久精品久久久久久噜噜老黄| av在线老鸭窝| 欧美一区二区亚洲| 国产有黄有色有爽视频| 麻豆精品久久久久久蜜桃| 国产不卡一卡二| 午夜精品一区二区三区免费看| 嫩草影院精品99| 伊人久久精品亚洲午夜| 国产在视频线精品| 免费看a级黄色片| 亚洲国产精品国产精品| 久久久久久久午夜电影| 亚洲色图av天堂| 男人爽女人下面视频在线观看| 只有这里有精品99| 久久久久久久久久久丰满| 亚洲国产欧美人成| 看十八女毛片水多多多| 水蜜桃什么品种好| 少妇被粗大猛烈的视频| 欧美成人精品欧美一级黄| 免费电影在线观看免费观看| 97超碰精品成人国产| 精品久久久久久久末码| 男女啪啪激烈高潮av片| 亚洲图色成人| 亚洲精品国产成人久久av| 大香蕉久久网| 熟女人妻精品中文字幕| 热99在线观看视频| 亚洲人成网站在线观看播放| 久久韩国三级中文字幕| 男的添女的下面高潮视频| 成人亚洲精品av一区二区| 久久久国产一区二区| 伊人久久精品亚洲午夜| 久久久久网色| 26uuu在线亚洲综合色| 大又大粗又爽又黄少妇毛片口| 嫩草影院新地址| 免费观看av网站的网址| 亚洲一区高清亚洲精品| 黑人高潮一二区| 国产伦一二天堂av在线观看| 99久久精品国产国产毛片| 亚洲自拍偷在线| 久久久久久久亚洲中文字幕| 视频中文字幕在线观看| 国产精品一区二区性色av| 欧美成人精品欧美一级黄| 国产亚洲91精品色在线| 国产毛片a区久久久久| 国产精品久久久久久精品电影小说 | 精品少妇黑人巨大在线播放| 老司机影院成人| 国产一区二区在线观看日韩| 国产欧美日韩精品一区二区| 天天躁日日操中文字幕| 久久精品人妻少妇| av黄色大香蕉| 国产精品女同一区二区软件| 永久网站在线| 蜜桃亚洲精品一区二区三区| 床上黄色一级片| 国产成人精品福利久久| 51国产日韩欧美| 草草在线视频免费看| 80岁老熟妇乱子伦牲交| 欧美成人一区二区免费高清观看| 亚洲真实伦在线观看| 老女人水多毛片| 一个人免费在线观看电影| 我的老师免费观看完整版| 国产91av在线免费观看| 国产一区二区亚洲精品在线观看| 一个人看的www免费观看视频| 大香蕉97超碰在线| 日韩欧美 国产精品| 天堂俺去俺来也www色官网 | 免费观看精品视频网站| 黑人高潮一二区| 日韩一区二区视频免费看| 亚洲国产欧美人成| 尾随美女入室| 日日摸夜夜添夜夜爱| h日本视频在线播放| 大香蕉97超碰在线| 男的添女的下面高潮视频| 亚洲欧美日韩卡通动漫| 波野结衣二区三区在线| 国产老妇伦熟女老妇高清| 成人毛片60女人毛片免费| 午夜爱爱视频在线播放| 只有这里有精品99| 91精品国产九色| 亚洲欧洲国产日韩| 99热这里只有是精品50| 激情五月婷婷亚洲| 国产精品99久久久久久久久| 国产成人a∨麻豆精品| 欧美xxⅹ黑人| 麻豆乱淫一区二区| 最新中文字幕久久久久| 99热这里只有是精品50| a级毛片免费高清观看在线播放| 国产成人一区二区在线| 毛片一级片免费看久久久久| 国产午夜福利久久久久久| 午夜激情福利司机影院| 日本午夜av视频| 乱码一卡2卡4卡精品| 精品久久国产蜜桃| 精品少妇黑人巨大在线播放| 一级毛片我不卡| ponron亚洲| 99久国产av精品| 亚洲精品乱久久久久久| 18+在线观看网站| 久久6这里有精品| 免费观看av网站的网址| 久久精品国产亚洲网站| 久久这里只有精品中国| 欧美潮喷喷水| 久久99热这里只频精品6学生| av天堂中文字幕网| 秋霞在线观看毛片| 777米奇影视久久| 人人妻人人澡人人爽人人夜夜 | 亚洲欧美一区二区三区黑人 |