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

    四川盆地大氣氨與氮氧化物排放對細顆粒物污染的影響及減排潛力

    2022-08-24 00:20:26楊顯玉呂雅瓊潘月鵬王式功張小玲易家俊喬玉紅
    中國環(huán)境科學 2022年8期
    關鍵詞:川南四川盆地氮氧化物

    楊顯玉,呂雅瓊,邵 平*,潘月鵬,樊 晉,王式功,張小玲,易家俊,喬玉紅,吳 鍇

    四川盆地大氣氨與氮氧化物排放對細顆粒物污染的影響及減排潛力

    楊顯玉1,呂雅瓊2,邵 平1*,潘月鵬3,樊 晉1,王式功1,張小玲1,易家俊1,喬玉紅4,5,吳 鍇1

    (1.成都信息工程大學大氣科學學院,高原大氣與環(huán)境四川省重點實驗室,四川 成都 610225;2.中國科學院成都山地災害與環(huán)境研究所,四川 成都 610041;3.中國科學院大氣物理研究所,大氣邊界層物理和大氣化學國家重點實驗室,北京 100029;4.重慶市生態(tài)環(huán)境科學研究院,重慶 401147;5.重慶市城市大氣環(huán)境綜合觀測與污染防控重點實驗室,重慶 401147)

    基于IASI、OMI與TROPOMI衛(wèi)星數(shù)據(jù)識別了2008~2019年四川盆地氨與氮氧化物柱濃度的變化趨勢,并進一步采用空氣質量模型CMAQ對2019年冬季四川盆地氨排放的大氣環(huán)境影響進行了研究,評估了氨與氮氧化物單獨減排及氨與氮氧化物協(xié)同減排情景下對四川盆地顆粒物污染的影響.結果表明:2008~2012年四川盆地氮氧化物排放逐年升高隨后在2013~2019年迅速下降,而氨柱濃度在2008~2013年期間較為穩(wěn)定,自2014年起迅速增長.四川盆地氨排放的高值區(qū)主要集中在人為活動強烈的成都及周邊地區(qū)和川南城市群以及農業(yè)源氨排放主導的川西北地區(qū).銨根離子在川南城市群的PM2.5當中占比高達11.4%,而對川西地區(qū)城市的PM2.5貢獻較低.敏感性實驗結果表明,氨與氮氧化物協(xié)同減排50%能有效降低大氣中硝酸銨與硫酸銨的濃度,從而減少細顆粒物污染,改善四川盆地區(qū)域環(huán)境空氣質量.

    氨;細顆粒物(PM2.5);四川盆地;WRF-CMAQ

    氮氧化物(NO)與氨(NH3)不僅是大氣中重要的痕量氣體,還是細顆粒物(PM2.5)形成的關鍵組分[1].其中,全球大氣NH3濃度的增加趨勢具有多重環(huán)境影響:作為全球氮循環(huán)的一部分[2],NH3的大氣沉降會導致水體富營養(yǎng)化,破壞生態(tài)系統(tǒng)平衡[3-4].此外,NH3在大氣中與酸性氣體快速反應,形成PM2.5中二次無機氣溶膠的主要組分—硫酸銨和硝酸銨[5-7].

    近幾十年來持續(xù)增長的NH3排放已成為我國乃至全球大氣環(huán)境治理的難題.自2013年大氣十條實施以來,我國NO和SO2排放總量已明顯下降,但對NH3的排放控制措施仍未得到充分實施.我國2007~2018年氨排放總量一直在10Mt上下浮動,NH3減排效果并不顯著[8].目前,國內外學者針對我國氨排放時空分布特征、年際排放變化、行業(yè)來源解析等方面已開展了一系列研究[9-11].然而,已有研究大多基于衛(wèi)星資料與地基觀測的方法針對我國京津冀、長江三角洲和珠江三角洲等地區(qū)開展NH3排放清單估算及其環(huán)境影響評估.董文煊等[12]采用排放因子法估算了1994~2006年中國各省份的人為源氨排放量,發(fā)現(xiàn)我國農業(yè)源氨排放在人為源氨排放中占到90%,其它排放源包括生物質燃燒、交通排放、化工生產等僅占10%.此外,人口密集的四川盆地是我國氨排放的熱點區(qū)域,僅2006年氨總排放量達到1115.4kt,僅次于河南(1783.6kt)、山東(1569.3kt).馮小瓊等[13]基于四川省農業(yè)氨活動水平數(shù)據(jù)采用排放因子法建立了2005~2014年農業(yè)源氨排放清單,指出四川省2005~2014年氨排放總量雖有下降,但2012年谷值排放也高達919kt.

    盡管已有研究基于自下而上與排放因子構建的排放清單對四川盆地氨排放量進行估算且分析氨排放來源,但對四川盆地的大氣氨濃度時空分布及變化特征的研究相對較少,尤其缺乏對氨減排政策所導致的大氣環(huán)境影響的評估[17-20].基于此,本文選取四川盆地為研究對象,結合衛(wèi)星資料、高分辨率排放清單與WRF-CMAQ空氣質量模型評估了氮氧化物與氨排放長期變化趨勢,通過比較不同程度的氮氧化物與氨排放削減比例下細顆粒物以及重要二次組分濃度的變化,分析四川盆地氨的減排潛力,以期為有效推進PM2.5污染治理和制定區(qū)域大氣污染協(xié)同管控策略提供有力科學依據(jù)和決策參考.

    1 數(shù)據(jù)與方法

    1.1 OMI與TROPOMI衛(wèi)星NO2柱濃度

    Ozone Monitoring Instrument (OMI)是美國國家航空航天局(NASA)于2004年發(fā)射的Aura衛(wèi)星所搭載的紫外可見光譜儀.OMI每天覆蓋全球1次,軌道掃描寬度為2600km,空間分辨率為13km×24km,過境時間約為當?shù)貢r間13:45.本研究采用由荷蘭皇家氣象研究院反演的QA4ECVversion 1.1OMI對流層NO2垂直柱濃度產品[21].

    TROPOspheric Monitoring Instrument (TROPOMI)是由歐洲航天局(ESA)為哥白尼哨兵5號發(fā)射,空間分辨率為7km×3.5km(2019年8月6日后分辨率為5.5km×3.5km)[24].

    1.2 IASI衛(wèi)星NH3柱濃度

    IASI-A儀器搭載在2006年發(fā)射的MetOp-A衛(wèi)星上,波譜覆蓋范圍為645~2769cm-1,每天覆蓋全球2次,空間分辨率約為12km×12km,并在當?shù)貢r間09:30和21:30穿過赤道,提供高精度大氣痕量氣體產品(NH3等).因日間數(shù)據(jù)誤差相比夜間數(shù)據(jù)較小,本文僅采用日間過境時間的IASI NH3測量值[25].

    1.3 模型設置

    本研究采用WRFv3.9.1-CMAQv5.3.2空氣質量模型對四川盆地城市群的大氣污染物進行模擬.模擬區(qū)域采用3層嵌套[16],空間分辨率分別為27,9和3km,其中最內層覆蓋四川盆地的主要區(qū)域,如圖1所示.WRF 模型的物理參數(shù)化方案如表1所示,其初始場與邊界場數(shù)據(jù)均采用美國國家環(huán)境預報中心(NCEP)提供的6h一次、1°分辨率的 FNL全球再分析資料.模擬起始時間為2018年12月29日00:00,結束時間為2019年3月1日00:00.為減少WRF模型數(shù)值模擬長期運行累積的系統(tǒng)誤差,本研究采用Spin-up時間設置為72h并每間隔2周初始化一次進行模擬.WRF共設置30層垂直氣壓層且頂層氣壓為50hPa,其中2km以內加密以更好地模擬邊界層內污染物傳輸.WRF與CMAQ模型均采用Lambert投影坐標系.CMAQ的氣相化學和氣溶膠機制分別采用CB06及AERO06機制.CMAQ最外層區(qū)域的初始場和邊界場均由全球大氣化學傳輸模型MOZART-4提供[26].人為源排放清單采用清華大學開發(fā)的2017年中國多尺度高分辨率排放清單MEIC (http://meicmodel. org/),并基于人口密度、交通路網(wǎng)及GDP分布數(shù)據(jù)利用ArcGIS對清單進行時空分配[27-28].天然源揮發(fā)性有機物排放通過MEGANv2.1模型模擬得到[29].

    表1 WRF物理參數(shù)化方案設置

    圖1 WRF-CMAQ模型最內層模擬區(qū)域

    審圖號:GS(2019)-1822

    1.4 敏感性實驗設置

    為評估不同的人為源減排策略對四川盆地細顆粒物污染的影響,本研究設定4個排放情景:基準情景,即不對排放清單進行調整;單一的NH3減排情景(I),即削減排放清單中50%的氨排放;單一的NO減排情景(II),即對排放清單中NO削減50%;NH3與NO的協(xié)同減排情景(III),即對排放清單中的氨與NO同步削減50%.通過CMAQ模型對不同削減策略下的四川盆地2019年冬季的PM2.5濃度及其化學組分進行模擬,對比分析不同排放情景下PM2.5及其化學組分的濃度變化,探究單一減排與協(xié)同減排策略對進一步減少PM2.5濃度的有效性.需要指出的是,NH3與NO協(xié)同減排50%是在碳中和政策愿景下結合清潔能源、發(fā)展電動汽車等關鍵技術發(fā)展的可達目標.

    1.5 模擬評估方法

    為評估WRF模型模擬氣象場的準確性,本研究采用中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/)的地面逐小時氣象觀測資料對WRF模擬的氣象要素進行評估. CMAQ模擬結果則利用四川省環(huán)境監(jiān)測總站的地面污染物監(jiān)測數(shù)據(jù)進行對比并計算統(tǒng)計指標.統(tǒng)計參量采用平均偏差(MB),標準化平均偏差(NMB),均方根誤差(RMSE),一致性指數(shù)(IOA).各個統(tǒng)計量的定義如下:

    式中:PO分別為各個要素的模擬值(μg/m3)與觀測值(μg/m3).

    2 結果與討論

    2.1 四川盆地2008~2019年NO2與NH3柱濃度長期變化趨勢

    歸因于快速的城市化與工業(yè)化進程以及“十一五”規(guī)劃期間對NO排放缺乏管控措施,2008~2012年期間四川盆地NO2呈顯著升高趨勢,且2012年年均NO2柱濃度高達3.2′1015molec/cm2(圖2).然而,該時段NH3柱濃度在6.9′1015~7.6′1015molec/cm2范圍內波動,變化較為平穩(wěn).自2013年“大氣十條”實施以來,四川盆地NO2柱濃度逐年顯著下降,2019年區(qū)域平均柱濃度為2.7′1015molec/cm2,較2012年峰值水平降低了15.6%.然而,NH3柱濃度呈現(xiàn)出與NO2柱濃度相反的變化趨勢,2013~2017年NH3柱濃度持續(xù)升高且2018~2019年均維持在較高水平.此外,2017年IASI NH3年均柱濃度高達12.5′1015molec/cm2,相比2012年柱濃度升高了78.0%.這一顯著的NH3柱濃度升高一方面可能歸因于2013~ 2019年對氨排放仍缺乏有效的管控策略,另一方面則可能是由于自2013年“大氣十條”頒布以來對SO2和NO2嚴格的排放控制導致NH3和酸性氣體之間的中和反應被削弱,導致增加大氣中的NH3濃度.此外,城市地區(qū)來源于機動車的NH3排放增加也可能對這一現(xiàn)象具有重要貢獻[30].然而,2013~2017年四川盆地NH3柱濃度迅速升高的成因仍有待進一步研究.

    2.2 四川盆地2019年NO2與NH3衛(wèi)星柱濃度及排放特征

    已有研究表明,四川盆地城市群NO排放主要由人為源排放主導,土壤、閃電等自然源排放的影響相對較小[31].由圖3可知,NO2柱濃度的高值區(qū)廣泛分布于四川盆地城市群且年均值大于5′1015molec/cm2的區(qū)域主要集中于各城市的主城區(qū).成都市及周邊的眉山市,德陽市,資陽市構成了位于四川盆地西部的污染高值區(qū),即成都平原污染帶.另一高值區(qū)則主要分布在重慶市主城區(qū)及周邊區(qū)域,表明四川盆地NO污染呈現(xiàn)明顯的區(qū)域性特征.川南城市群的內江市,自貢市,瀘州市和樂山市為盆地內的次高值區(qū)(約為4′1015molec/cm2).盆地北部的南充市和廣元市則較為清潔,僅城市中心區(qū)域NO2柱濃度較高.不同于NO2柱濃度的分布,IASI衛(wèi)星不僅在城市中心觀測到氨柱濃度的明顯高值(大于30′1015molec/ cm2),還在四川盆地各城市群的農村及郊區(qū)也呈現(xiàn)出廣泛的氨熱點排放分布,尤其是成都平原與川南城市群.這主要是因為農業(yè)源是四川盆地氨排放的主要貢獻源.徐晨曦等[32]利用人為源活動水平數(shù)據(jù)構建了四川地區(qū) 2017年人為源大氣污染物排放清單發(fā)現(xiàn)農業(yè)源的氨排放貢獻占比高達97.1%.值得注意的是,自下而上方法構建的氨排放清單顯示川西北地區(qū)因畜牧業(yè)發(fā)達具有較高的氨排放量,但IASI衛(wèi)星的NH3柱濃度并沒有很好地識別出川西北地區(qū)(阿壩藏族自治州)的氨排放熱點.這可能是由于衛(wèi)星的過境時間與柱濃度反演方法的不確定性所致.

    年份

    圖3 2019年四川盆地對流層NO2與NH3年均柱濃度空間分布

    進一步分析發(fā)現(xiàn)盡管四川盆地冬季NH3柱濃度明顯低于春季與夏季,但由于NO冬季生命周期較長,NO2柱濃度于冬季出現(xiàn)峰值,在不利氣象條件下對區(qū)域灰霾污染事件具有重要作用[33-34].此外,已有研究表明四川盆地區(qū)域PM2.5污染事件多發(fā)于冬季,其中成都平原與川南城市群較為嚴重[16].因此,本研究選取2019年冬季(1~2月)并采用CMAQ模型對四川盆地冬季氨氮排放減排的有效性進行進一步評估.

    2.3 模型性能評估

    由表2可知,WRF模型模擬的氣溫與觀測的氣溫呈現(xiàn)出高度一致性,平均偏差(MB)為-0.16℃且IOA高達0.96.WRF對相對濕度的模擬存在一定偏差,這主要是由于四川盆地冬季持續(xù)的高濕環(huán)境被WRF模型低估所致.此外,風速的平均偏差為-1.12m/s且RMSE為4.88m/s,表明WRF對10m風速的模擬有所高估,這可能是因為WRF模型的下墊面類型有所偏差以及WRF參數(shù)化方案誤差.

    表2 WRF模型模擬與氣象站點觀測的氣象要素對比評估

    圖4 2019年1~2月四川盆地各城市PM2.5觀測與CMAQ模擬結果對比

    CMAQ模擬的PM2.5濃度采用四川省環(huán)境監(jiān)測總站的地面污染物監(jiān)測數(shù)據(jù)進行評估,如圖4與圖5所示.各城市的模擬值與觀測值之間偏差相對較小,尤其是遂寧和內江的平均偏差(MB)僅為6.2與-1.9μg/m3.此外,CMAQ能較好地模擬2019年1~2月四川盆地各個城市PM2.5濃度變化趨勢,各城市IOA除成都外均高于0.6.需要指出的是,CMAQ對重污染事件期間的PM2.5濃度峰值存在一定程度的低估,這可能是由于排放清單的不確定性和WRF模擬的氣象場存在一定偏差所致[33].

    2.4 氮氧化物與氨排放對四川盆地冬季PM2.5濃度的影響

    由圖6可知,重慶市主城區(qū)、成都平原城市群及川南城市群顆粒物污染呈連片狀分布,冬季PM2.5濃度均高于70μg/m3,這主要是由于居高不下的工業(yè)源與交通源排放以及冬季頻繁的不利氣象條件(逆溫,靜穩(wěn)天氣形勢,靜小風頻繁等)所致[34].川東北城市群冬季PM2.5濃度則相對較低,僅達州市區(qū)出現(xiàn)濃度高值,郊區(qū)與農村的PM2.5濃度普遍低于35μg/m3.

    進一步分析PM2.5化學組分特征發(fā)現(xiàn),銨根離子占比較高的區(qū)域為川南城市群、成都平原以及川西北地區(qū),這一現(xiàn)象與氨排放的空間格局密切相關.具體而言,川西北地區(qū)農業(yè)發(fā)達,銨根離子主要源于畜牧業(yè)與養(yǎng)殖業(yè)的氨態(tài)肥料的廣泛使用,導致該區(qū)域氨排放量相對較高,PM2.5中銨根占比高達11.4%.相比于川西北地區(qū),成都及周邊地區(qū)的銨根離子不僅受到農業(yè)源的影響,還源于城市地區(qū)的交通源與工業(yè)排放.這一結果與Yang等[16]對成都市生態(tài)環(huán)境科學研究院超級站的PM2.5化學組分分析結果一致.周子航等[17]基于2015年四川省活動水平數(shù)據(jù)建立了2015年1km×1km人為源大氣排放清單,指出除農業(yè)源以外,由機動車尾氣排放占主導的道路移動源以及工藝過程源是成都及周邊地區(qū)氨的主要來源(占比約為14%).高濃度硝酸根離子(>20μg/m3)主要集中在城市中心,如成都市區(qū)、川南城市群、重慶市主城區(qū).這主要是由于城市中心的機動車尾氣以及工業(yè)源排放的氮氧化物水平較高,且加之氣象條件有利于硝酸根生成.需要注意的是,自2013年“大氣十條”實施以來,盡管MEIC清單與OMI衛(wèi)星NO2柱濃度數(shù)據(jù)均表明四川盆地NO總排放量在2013~2020年期間逐年下降,但城市地區(qū)機動車尾氣及工業(yè)源排放的NO排放量仍處于較高水平[32].王譽曉等[35]結合大氣污染物排放趨勢對2013~2017年四川省大氣污染防治成效進行了評估,指出2017年四川盆地NO的排放總量相較2013年減少了15.4×104t,但仍高達52.2×104t.此外,截止2019年底,國家統(tǒng)計局數(shù)據(jù)表明成都市機動車保有量已高達576萬輛,位居全國城市機動車保有量第2.因此交通源NO排放的進一步嚴格管控是四川盆地冬季灰霾污染治理的關鍵因素.

    2.5 氨氮減排情景對四川盆地冬季PM2.5濃度的影響

    由圖7a可知,單一的氨排放減少50%策略可對四川盆地冬季PM2.5濃度產生廣泛影響,導致川南城市群等氨排放熱點地區(qū)的月均PM2.5降低了4~ 6μg/m3.相比之下,削減50%的NO排放的作用則更為明顯,能有效減少四川盆地的PM2.5污染,尤其是川南城市群及重慶市主城區(qū)等氮氧化物排放高值區(qū)的PM2.5減少了10μg/m3(圖7b).然而,對于川東北城市群與成都及周邊地區(qū)的影響相對較小,平均變化為5μg/m3.該現(xiàn)象可能是由于這些地區(qū)氮氧化物排放量過高導致即使削減了50%的NO排放的情況下,大氣中仍有足夠的NO被氧化并在富氨環(huán)境下中和生成硝酸鹽.分析NH3和NO協(xié)同減排情景下PM2.5的濃度變化可知(圖7c),同步削減50%的氨和NO導致四川盆地2019年1~2月的平均PM2.5濃度降低了25%,其中川南城市群與重慶市主城區(qū)PM2.5濃度減幅高達14~16μg/m3,減排效果最為顯著,且模擬時段內PM2.5峰值較基準情景下降了20%,實現(xiàn)了PM2.5削峰的作用.然而,協(xié)同減排情景對成都市、德陽市、綿陽市及川東北城市群影響并不明顯,僅降低4~6μg/m3的PM2.5濃度.此外,即使在氨和NO協(xié)同減排50%的情景下,冬季四川盆地多個城市的PM2.5仍處于較高水平(月均濃度高于50μg/m3).進一步分析表明盡管氨和NO協(xié)同減排導致PM2.5中二次無機組分的濃度顯著下降,但有機氣溶膠組分因VOCs排放保持在較高水平而成為PM2.5的主要化學組分.因此,在協(xié)同減排氨和NO的同時應進一步加強對VOCs的排放控制.

    本研究評估了協(xié)同控制氨與氮氧化物排放情景下對PM2.5污染的有效性,強調了氨排放的減排潛力及必要性.然而,值得注意的是,氨排放的減少可能導致酸雨問題加劇[36],且城市地區(qū)的氨來源復雜,已有研究表明機動車尾氣的氨排放已逐漸成為城市大氣中氨的主導來源[30].因此,氨減排導致的環(huán)境影響及氨排放清單的精準估算仍有待進一步研究.

    3 結論

    3.1 四川盆地氨排放的高值區(qū)主要集中在人為活動強烈的成都及周邊地區(qū)和川南城市群以及農業(yè)源氨排放主導的川西北地區(qū).由于富氨環(huán)境的影響,銨根離子在川南城市群的PM2.5當中占比較高,高達11.4%.而對川西地區(qū)的PM2.5貢獻較低.

    3.2 在氨與氮氧化物同步減少50%的情形下,四川盆地2019年冬季PM2.5濃度減少了25%.說明該策略能有效降低PM2.5中的硝酸銨與硫酸銨等二次無機組分濃度,從而減少細顆粒物污染,改善四川盆地冬季灰霾污染.

    [1] Erisman J W,Sutton M A,Galloway J N,et al. How a century of ammonia synthesis changed the world [J]. Nature Geoscience,2008,1:636–639.

    [2] Fowler D,Coyle M,Skiba U,et al. The global nitrogen cycle in the twenty-first century: introduction [J]. Philosophical transactions of the Royal Society of London. Series B,Biological Sciences,2013,368 (1621).

    [3] Krupa S V. Effects of atmospheric ammonia (NH3) on terrestrial vegetation: A review [J]. Environmental Pollution,2003,124(2):179- 221.

    [4] Camargo J A,Alonso á. Ecological and toxicological effects of inorganic nitrogen pollution in aquatic ecosystems: A global assessment [J]. Environment International,2006,32(6):831-849.

    [5] Aneja V P,Schlesinger W H,Erisman J W. Effects of agriculture upon the air quality and climate: Research,policy,and regulations [J]. Environmental Science & Technology,2009,43(12):4234–4240.

    [6] Wu K,Zhu S P,Liu Y M,et al. Modeling ammonia and its uptake by secondary organic aerosol over China. [J]. Journal of Geophysical Research: Atmospheres,2021,126(7).

    [7] Nowak J B,Neuman J A,Bahreini R,et al. Ammonia sources in the California South Coast air basin and their impact on ammonium nitrate formation [J]. Geophysical Research Letters,2012,39(7).

    [8] 劉學軍,沙志鵬,宋 宇,等.我國大氣氨的排放特征、減排技術與政策建議 [J]. 環(huán)境科學研究,2021,34(1):149-157.

    Liu X J,Sha Z P,Song Y,et al. China's atmospheric ammonia emission characteristics,mitigation options and policy recommendations [J].Research of Environmental Sciences,2021,34(1):149-157.

    [9] Zhang L,Chen Y,Zhao Y,et al. Agricultural ammonia emissions in China: Reconciling bottom-up and top-down estimates [J]. Atmospheric Chemistry and Physics,2018,18(1):339-355.

    [10] Kong L,Tang X,Zhu J,et al. Improved Inversion of Monthly Ammonia Emissions in China Based on the Chinese Ammonia Monitoring Network and Ensemble Kalman Filter [J]. Environmental Science & Technology,2019,53(21):12529–12538.

    [11] Liu L,Zhang X Y,Anthony Y. Estimating global surface ammonia concentrations inferred from satellite retrievals [J]. Atmospheric Chemistry and Physics,2019,19:12051-12066.

    [12] 董文煊,邢 佳,王書肖.1994~2006年中國人為源大氣氨排放時空分布 [J]. 環(huán)境科學,2010,31(7):1457-1463.

    Dong W X,Xing J,Wang S X.Temporal and spatial distribution of anthropogenic ammonia emissions in China: 1994~2006 [J]. Environmental Science,2010,31(7):1457-1463.

    [13] 馮小瓊,陳軍輝,姜 濤,等.四川省2005~2014年農業(yè)源氨排放清單及分布特征 [J]. 環(huán)境科學,2017,38(7):2728-2737.

    Feng X Q,Chen J H,Jiang T,et al. Agricultural ammonia emission inventory and its spatial distribution in Sichuan Province from 2005 to 2014 [J]. Environmental Science,2017,38(7):2728-2737.

    [14] 曹庭偉,吳 鍇,康 平,等.成渝城市群臭氧污染特征及影響因素分析 [J]. 環(huán)境科學學報,2018,38(4):1275-1284.

    Cao T W,Wu K,Kang P,et al. Study on ozone pollution characteristics and meteorological cause of Chengdu-Chongqing urban agglomeration [J]. Acta Scientiae Circumstantiae,2018,38(4): 1275-1284.

    [15] 吳 鍇,康 平,于 雷,等.2015~2016年中國城市臭氧濃度時空變化規(guī)律研究 [J]. 環(huán)境科學學報,2018,38(6):2179-2190.

    Wu K,Kang P,Yu L,et al. Pollution status and spatio-temporal variations of ozone in China during 2015~2016 [J]. Acta Scientiae Circumstantiae,2018,38(6):2179-2190.

    [16] Yang X Y,Lu Y Q,Zhu X S,et al. Formation and evolution mechanisms of severe haze pollution in the Sichuan Basin,Southwest China [J]. Aerosol and Air Quality Research,2020,20(11):2557–2567.

    [17] 周子航,鄧 也,譚欽文,等.四川省人為源大氣污染物排放清單及特征 [J]. 環(huán)境科學,2018,39(12):5344-5358.

    Zhou Z H,Deng Y,Tan Q W,et al.Emission inventory and characteristics of anthropogenic air pollutant sources in the Sichuan Province [J]. Environmental Science,2018,39(12):5344-5358.

    [18] 馮小瓊,陳軍輝,姜 濤,等.四川盆地2018年畜禽養(yǎng)殖氨排放清單研究 [J]. 四川環(huán)境,2020,39(5):24-31.

    Feng X Q,Chen J H,Jiang T,et al.,Study on livestock ammonia emission inventory of Sichuan Basin in 2018 [J]. Sichuan Environment,2020,39(5):24-31.

    [19] 李富春,韓圣慧,楊 俊,等.川渝地區(qū)農業(yè)生態(tài)系統(tǒng)NH3排放 [J]. 環(huán)境科學,2009,30(10):2823-2831.

    Li F C,Han S H,Yang J,et al. Agro-ecosystem ammonia emission in Sichuan-Chongqing Region [J]. Environmental Science,2009,30(10): 2823-2831.

    [20] 馮小瓊,王幸銳,何 敏,等.四川省2012年人為源氨排放清單及分布特征 [J]. 環(huán)境科學學報,2015,35(2):394-401.

    Feng X Q,Wang X R,He M,et al. A 2012-based anthropogenic ammonia emission inventory and its spatial distribution in Sichuan Province [J]. Acta Scientiae Circumstantiae,2015,35(2):394-401.

    [21] Boersma K,Eskes H,Richter A,et al. Improving algorithms and uncertainty estimates for satellite NO2retrievals: results from the quality assurance for the essential climate variables (QA4ECV) project [J]. Atmospheric Measurement Techniques,2018,11:6651-6678.

    [22] van Geffen J,Boersma K F,Eskes H,et al. S5P TROPOMI NO2slant column retrieval: method,stability,uncertainties and comparisons with OMI [J]. Atmospheric Measurement Techniques,2020,13(3):1315– 1335.

    [23] Veefkind J P,Aben I,Mcmullan K,et al. TROPOMI on the ESA Sentinel-5Precursor: A GMES mission for global observations of the atmospheric composition for climate,air quality and ozone layer applications [J]. Remote Sensing of Environment,2012,120:70-83.

    [24] Goldberg D L,Lu Z F,David G S,et al. Enhanced capabilities of TROPOMI NO2: Estimating NOfrom North American Cities and Power Plants [J]. Environmental Science & Technology,2019,53(21): 12594-12601.

    [25] Van Damme M,Clarisse L,Whitburn S,et al. Industrial and agricultural ammonia point sources exposed [J]. Nature,2018,564: 99-103.

    [26] Emmons L K,Walters S,Hess P G,et al. Description and evaluation of the model for ozone and related chemical tracers,Version 4 (MOZART-4) [J]. Geoscience Model Development,2010,25.

    [27] 鄭君瑜,張禮俊,鐘流舉,等.珠江三角洲大氣面源排放清單及空間分布特征 [J]. 中國環(huán)境科學,2009,29(5):455-460.

    Zheng J,Zhang L,Zhong L,et al. Area source emission inventory of air pollutant and its spatial distribution characteristics in Pearl River Delta [J]. China Environmental Science,2009,29(5):455-460.

    [28] Zheng B,Tong D,Li M,et al. Trends in China’s anthropogenic emissions since 2010 as the consequence of clean air actions [J]. Atmospheric Chemistry and Physics,2018,18(19):14095-14111.

    [29] Guenther A,Jiang X,Heald C L,et al. The model of emissions of gases and aerosols from nature version 2.1 (MEGAN2.1): An extended and updated framework for modeling biogenic emissions. [J]. Geoscientific Model Development,2012,5(6):1471-1492.

    [30] Farren N,Davison J,Rose R,et al. Underestimated ammonia emissions from road vehicles [J]. Environmental Science & Technology,2020,54,15689-15697.

    [31] Wu K,Wang Y R,Qiao Y H,et al. Drivers of 2013-2020ozone trends in the Sichuan Basin,China: Impacts of meteorology and precursor emission changes [J]. Environmental Pollution,2022,300,118914.

    [32] 徐晨曦,陳軍輝,李 媛,等.四川省基于第二次污染源普查數(shù)據(jù)的人為源大氣污染源排放清單及特征 [J]. 環(huán)境科學,2020,41(10): 4482-4494.

    Xu C X,Chen J H,Li Y,et al.Emission inventory and characteristics of anthropogenic air pollution sources based on second pollution source census data in Sichuan Province [J]. Environmental Science,2020,41(10):4482-4494.

    [33] 楊顯玉,易家俊,呂雅瓊,等.成都市及周邊地區(qū)嚴重臭氧污染過程成因分析 [J]. 中國環(huán)境科學,2020,40(5):2000-2009.

    Yang X Y,Yi J J,Lu Y Q,et al. Characteristics and formation mechanism of a severe O3episode in Chengdu and surrounding areas [J]. China Environmental Science,2020,40(5):2000-2009.

    [34] 危詩敏,馮鑫媛,王式功,等.四川盆地多層逆溫特征及其對大氣污染的影響 [J]. 中國環(huán)境科學,2021,41(3):1005-1013.

    Wei S M,Feng X Y,Wang S G.,et al. Characteristics of multi-layer inversions in Sichuan Basin and their influences on air pollution [J]. China Environmental Science,2021,41(3):1005-1013.

    [35] 王譽曉,張云秋,王式功.2013~2017年四川省大氣污染防治成效評估 [J]. 蘭州大學學報(自然科學版),2020,56(3):388-395.

    Wang Y X,Zhang Y Q,Wang S G,et al. Results and analysis of air pollution prevention and control in Sichuan Province from 2013 to 2017 [J]. Journal of Lanzhou University,2020,56(3):388-395.

    [36] Liu M X,Huang X,Song Y,et al. Ammonia emission control in China would mitigate haze pollution and nitrogen deposition,but worsen acid rain [J]. Proceedings of the National Academy of Sciences,2019,116(16):7760-7765.

    致謝:感謝清華大學MEIC團隊提供的多尺度高分辨率排放清單.感謝四川省生態(tài)環(huán)境科學研究院對本研究的合作支持.

    Impacts of ammonia and nitrogen oxides emissions on fine particle matter pollution and mitigation potential over the Sichuan Basin.

    YANG Xian-yu1,LYU Ya-qiong2,SHAO Ping1*,PAN Yue-peng3,FAN Jin1,WANG Shi-gong1,ZHANG Xiao-ling1,YI Jia-jun1,QIAO Yu-hong4,5,WU Kai1

    (1.Plateau Atmosphere and Environment Key Laboratory of Sichuan Province,School of Atmospheric Sciences,Chengdu University of Information Technology,Chengdu 610225,China;2.Institute of Mountain Hazards and Environment,Chinese Academy of Sciences,Chengdu 610041,China;3.State Key Laboratory of Atmospheric Boundary Layer Physics and Atmospheric Chemistry (LAPC),Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029,China;4.Chongqing Institute of Eco-Environmental Science,Chongqing 401147,China;5.Chongqing Key Laboratory of Urban Atmospheric Environment Observation and Pollution Control,Chongqing 401147,China).,2022,42(8):3502~3511

    In this study,satellite data derived from IASI,OMI and TROPOMI were adopted to infer the trend of ammonia and NOemissions over 2008~2019 and the weather research forecast model (WRF) and the Community Multiscale Air Quality (CMAQ) were used to assess the environmental consequence of ammonia emissions over the Sichuan Basin in the winter of 2019. Furthermore,sensitivity experiments were performed to evaluate impacts of ammonia or NOemission reduction and ammonia and nitrogen oxides coordinated emission reduction scenarios on fine particles. The results show that nitrogen oxides emissions increased interannually from 2008 to 2012 and then decreased rapidly from 2013 to 2019,while ammonia columns were relatively stable over 2008~2013 and increased significantly from 2014 onwards. The hotspots of ammonia emissions in the Sichuan Basin were mainly concentrated in southern Sichuan city cluster and Chengdu and surrounding areas where anthropogenic activities were intense. The elevated ammonia emissions were also observed in the northwestern Sichuan city cluster dominated by agricultural sources. Ammonia ions accounted for a relatively high fraction of PM2.5in southern Sichuan city cluster,up to 11.4%,while its contribution to PM2.5in the cities of western Sichuan were relatively low. Sensitivity experiments indicate that the coordinated emission regulation in simultaneously reducing 50% NOand ammonia could effectively diminish the concentrations of ammonium nitrate and ammonium sulfate in PM2.5,therefore reduce fine particulate pollution and improve the regional ambient air quality across the Sichuan Basin.

    ammonia;fine particle matter(PM2.5);Sichuan Basin;WRF-CMAQ

    X513

    A

    1000-6923(2022)08-3502-10

    2022-01-24

    國家自然科學基金資助項目(42175174,42105167);國家重點研發(fā)計劃(2017YFC0212700);高原大氣與環(huán)境四川省重點實驗室開放課題(PAEKL-2020-C6)

    * 責任作者,講師,pingshao@cuit.edu.cn

    楊顯玉(1982-),男,副教授,主要從事大氣環(huán)境化學與陸氣相互作用研究.發(fā)表論文15篇.

    猜你喜歡
    川南四川盆地氮氧化物
    四川盆地極端降水演變特征及擬合
    大數(shù)據(jù)在油氣勘探開發(fā)中的應用——以川南頁巖氣田為例
    礦產勘查(2020年9期)2020-12-25 02:53:52
    川南地區(qū)五峰組-龍馬溪組3500米以淺探明萬億方頁巖氣田
    礦產勘查(2020年5期)2020-12-19 18:25:11
    川南頁巖氣水平井鉆井技術難點與對策
    川南古敘礦區(qū)煤層氣資源有利區(qū)優(yōu)選
    四川盆地秋季氣溶膠與云的相關分析
    四川盆地土壤濕度時空分布及影響因子分析
    四川盆地西南部短時強降水天氣特征分析
    通過高壓直接噴射實現(xiàn)高效率和低氮氧化物排放的氫燃燒方式
    氮氧化物吸附催化器與選擇性催化還原裝置組合排放控制系統(tǒng)用的先進催化劑
    精品人妻熟女av久视频| 青春草亚洲视频在线观看| 欧美日韩精品成人综合77777| 亚洲欧美精品专区久久| 亚洲精品日韩av片在线观看| 中国三级夫妇交换| 久久久国产一区二区| 夫妻午夜视频| 欧美成人a在线观看| 最近手机中文字幕大全| 国产人妻一区二区三区在| 亚洲国产精品成人久久小说| 国产视频首页在线观看| 一个人看视频在线观看www免费| 亚洲最大成人中文| 成人亚洲精品一区在线观看 | 久久久亚洲精品成人影院| 久久久久久久国产电影| 欧美日韩亚洲高清精品| 国产欧美日韩精品一区二区| 看免费成人av毛片| 天天躁夜夜躁狠狠久久av| 日本wwww免费看| 国产一区二区亚洲精品在线观看| 欧美成人a在线观看| 欧美xxxx性猛交bbbb| 麻豆精品久久久久久蜜桃| 国产乱人偷精品视频| 熟女av电影| 我的老师免费观看完整版| 精品人妻视频免费看| 日本wwww免费看| 国产av不卡久久| 一级毛片电影观看| 九九久久精品国产亚洲av麻豆| 国产综合懂色| 国产亚洲最大av| 日韩av免费高清视频| 国产精品.久久久| 麻豆久久精品国产亚洲av| 日韩,欧美,国产一区二区三区| 成人毛片a级毛片在线播放| 看非洲黑人一级黄片| 黑人高潮一二区| 卡戴珊不雅视频在线播放| 久久精品国产亚洲网站| 亚洲在线观看片| 99久久九九国产精品国产免费| 日韩人妻高清精品专区| 国产精品一及| 18禁在线播放成人免费| 日韩 亚洲 欧美在线| 国产亚洲av片在线观看秒播厂| 九草在线视频观看| 男女啪啪激烈高潮av片| 亚洲欧美成人综合另类久久久| 久久久成人免费电影| 亚洲精品乱久久久久久| 一级毛片黄色毛片免费观看视频| 久热久热在线精品观看| 久久久久性生活片| 啦啦啦啦在线视频资源| 亚洲久久久久久中文字幕| 精品一区二区三卡| 高清视频免费观看一区二区| 97在线视频观看| 男女啪啪激烈高潮av片| 菩萨蛮人人尽说江南好唐韦庄| 噜噜噜噜噜久久久久久91| 国产黄片视频在线免费观看| 色婷婷久久久亚洲欧美| 国产91av在线免费观看| 国产在线男女| 国产成人freesex在线| 永久免费av网站大全| 精品一区二区免费观看| 久久ye,这里只有精品| 卡戴珊不雅视频在线播放| 乱系列少妇在线播放| 青青草视频在线视频观看| 精品一区在线观看国产| 国产视频首页在线观看| 男女国产视频网站| 国产视频首页在线观看| 精品久久久久久久末码| 狂野欧美白嫩少妇大欣赏| 蜜桃久久精品国产亚洲av| 80岁老熟妇乱子伦牲交| 不卡视频在线观看欧美| 三级男女做爰猛烈吃奶摸视频| 99精国产麻豆久久婷婷| 高清毛片免费看| 久久久久久久久久久丰满| 国产精品熟女久久久久浪| 人妻系列 视频| 一级av片app| 91精品一卡2卡3卡4卡| 男女边摸边吃奶| 国产成人精品婷婷| 午夜福利视频精品| 麻豆精品久久久久久蜜桃| 99热这里只有精品一区| 精品人妻一区二区三区麻豆| 亚洲不卡免费看| 97在线人人人人妻| 午夜免费男女啪啪视频观看| 伊人久久精品亚洲午夜| 九草在线视频观看| 大话2 男鬼变身卡| 免费高清在线观看视频在线观看| 国产国拍精品亚洲av在线观看| 五月天丁香电影| 免费黄频网站在线观看国产| 国产黄片美女视频| 亚洲自拍偷在线| xxx大片免费视频| 又大又黄又爽视频免费| 久久久午夜欧美精品| 特大巨黑吊av在线直播| 色哟哟·www| 夫妻午夜视频| 国产成人a∨麻豆精品| 中文乱码字字幕精品一区二区三区| 成年免费大片在线观看| tube8黄色片| 国产精品99久久久久久久久| 欧美日韩综合久久久久久| 国产永久视频网站| 岛国毛片在线播放| h日本视频在线播放| 在线观看免费高清a一片| 亚洲最大成人av| 亚洲av免费高清在线观看| 尤物成人国产欧美一区二区三区| 不卡视频在线观看欧美| av.在线天堂| 又大又黄又爽视频免费| 少妇的逼好多水| 久久久久久久精品精品| 亚洲精品国产av成人精品| 亚洲av不卡在线观看| 亚洲最大成人中文| 美女主播在线视频| 国产在视频线精品| 久久久久国产精品人妻一区二区| 久久午夜福利片| 欧美zozozo另类| 人妻少妇偷人精品九色| 尤物成人国产欧美一区二区三区| 99久久九九国产精品国产免费| 免费观看av网站的网址| 国产 精品1| 久久99热6这里只有精品| 又爽又黄a免费视频| 免费看不卡的av| 国产白丝娇喘喷水9色精品| 亚洲真实伦在线观看| 爱豆传媒免费全集在线观看| 免费观看性生交大片5| 日本欧美国产在线视频| 性色avwww在线观看| 欧美bdsm另类| 搡老乐熟女国产| 国模一区二区三区四区视频| 99热这里只有是精品在线观看| 建设人人有责人人尽责人人享有的 | 久久精品久久精品一区二区三区| 国产成人aa在线观看| 一级爰片在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲成人中文字幕在线播放| 超碰av人人做人人爽久久| 国产乱来视频区| 欧美变态另类bdsm刘玥| 午夜福利网站1000一区二区三区| 亚洲av免费在线观看| 一本久久精品| 精品亚洲乱码少妇综合久久| 亚洲色图av天堂| 精品国产乱码久久久久久小说| 午夜免费鲁丝| 国产69精品久久久久777片| 国产精品人妻久久久久久| 神马国产精品三级电影在线观看| 精品久久久久久久久av| 亚洲最大成人手机在线| 男女那种视频在线观看| 成人亚洲欧美一区二区av| 又黄又爽又刺激的免费视频.| 国产综合懂色| 欧美激情久久久久久爽电影| 国产午夜精品久久久久久一区二区三区| 久久人人爽人人爽人人片va| 97在线视频观看| 美女cb高潮喷水在线观看| 女的被弄到高潮叫床怎么办| 伊人久久精品亚洲午夜| 免费观看a级毛片全部| 听说在线观看完整版免费高清| 中文乱码字字幕精品一区二区三区| 亚洲美女搞黄在线观看| 欧美激情国产日韩精品一区| 国产真实伦视频高清在线观看| 久久久久久久午夜电影| 国产老妇伦熟女老妇高清| 国产色婷婷99| 精品国产露脸久久av麻豆| 少妇人妻久久综合中文| 一区二区三区免费毛片| 嘟嘟电影网在线观看| 日韩成人av中文字幕在线观看| 久久精品久久久久久噜噜老黄| 成人特级av手机在线观看| 免费观看无遮挡的男女| 三级国产精品片| 亚洲国产日韩一区二区| 亚洲av.av天堂| 在线免费十八禁| 欧美人与善性xxx| 成人国产麻豆网| 春色校园在线视频观看| 国产欧美日韩精品一区二区| 国产毛片在线视频| 特级一级黄色大片| 国产精品一及| 国产片特级美女逼逼视频| 欧美老熟妇乱子伦牲交| 亚洲精品成人久久久久久| 午夜免费鲁丝| 亚洲精品国产av蜜桃| 中国美白少妇内射xxxbb| 国产爱豆传媒在线观看| av又黄又爽大尺度在线免费看| 欧美xxxx黑人xx丫x性爽| 高清毛片免费看| 国产老妇女一区| 少妇丰满av| 啦啦啦中文免费视频观看日本| 亚洲精品成人av观看孕妇| 人人妻人人爽人人添夜夜欢视频 | 97超视频在线观看视频| 成人无遮挡网站| 日韩 亚洲 欧美在线| 99九九线精品视频在线观看视频| 91久久精品国产一区二区三区| 亚洲成人精品中文字幕电影| 一边亲一边摸免费视频| 精品一区二区三区视频在线| 联通29元200g的流量卡| 永久免费av网站大全| av又黄又爽大尺度在线免费看| 国产av国产精品国产| 涩涩av久久男人的天堂| 成年女人看的毛片在线观看| av免费在线看不卡| 国产成人a∨麻豆精品| 日韩一区二区三区影片| 伊人久久精品亚洲午夜| 3wmmmm亚洲av在线观看| 亚洲精品国产成人久久av| 五月天丁香电影| 亚洲人成网站在线播| 亚洲最大成人中文| 欧美亚洲 丝袜 人妻 在线| 日日撸夜夜添| 免费不卡的大黄色大毛片视频在线观看| 成人黄色视频免费在线看| 最近中文字幕高清免费大全6| 亚洲久久久久久中文字幕| 国产欧美日韩精品一区二区| 欧美3d第一页| 成人亚洲欧美一区二区av| 欧美精品国产亚洲| 少妇人妻精品综合一区二区| 欧美97在线视频| 国产爱豆传媒在线观看| 久久人人爽av亚洲精品天堂 | 最近手机中文字幕大全| 国产精品人妻久久久久久| 亚洲欧美日韩无卡精品| 少妇裸体淫交视频免费看高清| videossex国产| 男人舔奶头视频| 自拍欧美九色日韩亚洲蝌蚪91 | 麻豆久久精品国产亚洲av| 精品一区在线观看国产| 日本黄色片子视频| 黄色一级大片看看| 久久99蜜桃精品久久| 男男h啪啪无遮挡| 国产欧美日韩精品一区二区| 少妇人妻 视频| 国产成人a∨麻豆精品| 国产91av在线免费观看| 国产熟女欧美一区二区| 国产成年人精品一区二区| 午夜激情福利司机影院| 五月开心婷婷网| 国产成人福利小说| 亚洲欧美日韩卡通动漫| 在线亚洲精品国产二区图片欧美 | 亚洲人与动物交配视频| 老师上课跳d突然被开到最大视频| 波野结衣二区三区在线| 久久久成人免费电影| 国产老妇伦熟女老妇高清| 中文资源天堂在线| 18禁裸乳无遮挡动漫免费视频 | 亚洲精品aⅴ在线观看| 久久久久久久午夜电影| av黄色大香蕉| 久久99蜜桃精品久久| 高清欧美精品videossex| 亚洲性久久影院| 国产成年人精品一区二区| 男人狂女人下面高潮的视频| 大话2 男鬼变身卡| 国产黄片视频在线免费观看| 日本三级黄在线观看| 男女啪啪激烈高潮av片| 亚洲aⅴ乱码一区二区在线播放| videos熟女内射| 美女被艹到高潮喷水动态| 青春草国产在线视频| 亚洲经典国产精华液单| 日韩强制内射视频| 涩涩av久久男人的天堂| 国产精品无大码| 成人无遮挡网站| 免费人成在线观看视频色| 午夜日本视频在线| 成人无遮挡网站| 熟女av电影| 七月丁香在线播放| 国产精品99久久99久久久不卡 | 国产伦在线观看视频一区| 寂寞人妻少妇视频99o| 久久精品综合一区二区三区| 国产精品蜜桃在线观看| 国产高清不卡午夜福利| 丰满乱子伦码专区| 老师上课跳d突然被开到最大视频| 男人舔奶头视频| 欧美xxxx黑人xx丫x性爽| 亚洲三级黄色毛片| 2021天堂中文幕一二区在线观| 亚洲,一卡二卡三卡| 亚洲国产成人一精品久久久| 波野结衣二区三区在线| 中国三级夫妇交换| 亚洲一区二区三区欧美精品 | 少妇人妻一区二区三区视频| 久久99蜜桃精品久久| 2018国产大陆天天弄谢| 国产精品一二三区在线看| 美女脱内裤让男人舔精品视频| 亚洲国产欧美人成| 国产69精品久久久久777片| 最近中文字幕2019免费版| 一级毛片电影观看| 欧美日韩亚洲高清精品| 国产女主播在线喷水免费视频网站| 午夜福利高清视频| 人体艺术视频欧美日本| 国产探花极品一区二区| 国产精品av视频在线免费观看| 高清欧美精品videossex| 夜夜看夜夜爽夜夜摸| 观看美女的网站| 一级毛片aaaaaa免费看小| 亚洲欧美成人综合另类久久久| 亚洲精品第二区| 成人亚洲欧美一区二区av| 黄色怎么调成土黄色| 日韩一区二区视频免费看| 亚洲欧洲国产日韩| 久久精品夜色国产| 插逼视频在线观看| 我要看日韩黄色一级片| 亚洲天堂国产精品一区在线| 国产老妇女一区| 小蜜桃在线观看免费完整版高清| 神马国产精品三级电影在线观看| 久久久久九九精品影院| www.av在线官网国产| 国产亚洲午夜精品一区二区久久 | 国产又色又爽无遮挡免| 国产一区二区三区av在线| 日韩欧美精品免费久久| 精品久久久久久久人妻蜜臀av| 亚洲熟女精品中文字幕| 亚洲最大成人av| 亚洲精品国产色婷婷电影| 国产视频内射| 99久久九九国产精品国产免费| 亚洲国产欧美人成| 亚洲图色成人| av.在线天堂| 国产精品精品国产色婷婷| 真实男女啪啪啪动态图| 哪个播放器可以免费观看大片| 亚洲欧美清纯卡通| 黄色配什么色好看| 久久久久久伊人网av| 麻豆乱淫一区二区| 波野结衣二区三区在线| 免费大片18禁| 亚洲av成人精品一区久久| 国内揄拍国产精品人妻在线| 欧美丝袜亚洲另类| 亚洲av日韩在线播放| 亚洲国产最新在线播放| 精品酒店卫生间| 观看美女的网站| 欧美日本视频| 亚洲国产欧美人成| 日韩电影二区| 亚洲av免费在线观看| 内地一区二区视频在线| 久久久亚洲精品成人影院| 国产精品爽爽va在线观看网站| 一边亲一边摸免费视频| 麻豆乱淫一区二区| 国产成人精品久久久久久| 日本av手机在线免费观看| 国产精品麻豆人妻色哟哟久久| 日产精品乱码卡一卡2卡三| 好男人在线观看高清免费视频| 伦理电影大哥的女人| 丰满乱子伦码专区| 国产亚洲av片在线观看秒播厂| 亚洲精品,欧美精品| 日本午夜av视频| 久久久久久久久大av| 大香蕉97超碰在线| www.av在线官网国产| 舔av片在线| 乱码一卡2卡4卡精品| 91aial.com中文字幕在线观看| 国产在视频线精品| 人妻制服诱惑在线中文字幕| 国产国拍精品亚洲av在线观看| 狠狠精品人妻久久久久久综合| 亚洲电影在线观看av| 天天躁日日操中文字幕| 国产午夜精品一二区理论片| 亚洲国产精品专区欧美| 18禁动态无遮挡网站| 国产av码专区亚洲av| h日本视频在线播放| 一边亲一边摸免费视频| 久久国内精品自在自线图片| 国产成人freesex在线| 午夜福利在线在线| 青春草亚洲视频在线观看| 日本一本二区三区精品| 久久热精品热| 国产真实伦视频高清在线观看| 综合色丁香网| 美女xxoo啪啪120秒动态图| 蜜桃久久精品国产亚洲av| 中文天堂在线官网| 亚洲精品影视一区二区三区av| 精品久久国产蜜桃| 亚洲精品日韩av片在线观看| 日韩亚洲欧美综合| h日本视频在线播放| 免费观看性生交大片5| 国产伦精品一区二区三区四那| 一级毛片 在线播放| 成人亚洲精品av一区二区| 亚洲欧美日韩无卡精品| 又爽又黄无遮挡网站| 97人妻精品一区二区三区麻豆| 成人鲁丝片一二三区免费| 中文字幕av成人在线电影| 日本欧美国产在线视频| 夫妻性生交免费视频一级片| 国产亚洲最大av| 精品久久久久久久人妻蜜臀av| 一本色道久久久久久精品综合| 国产91av在线免费观看| 国产午夜精品久久久久久一区二区三区| 丝袜脚勾引网站| 日本免费在线观看一区| 丝袜美腿在线中文| 99re6热这里在线精品视频| 九色成人免费人妻av| 国产成人精品婷婷| 下体分泌物呈黄色| av在线蜜桃| 亚洲欧美中文字幕日韩二区| 国产有黄有色有爽视频| 日日摸夜夜添夜夜爱| 亚洲欧美日韩东京热| 天美传媒精品一区二区| 在线播放无遮挡| 欧美国产精品一级二级三级 | 亚洲三级黄色毛片| 国产精品秋霞免费鲁丝片| 在线亚洲精品国产二区图片欧美 | 青青草视频在线视频观看| 免费看不卡的av| 亚洲av电影在线观看一区二区三区 | 国产男女超爽视频在线观看| av国产精品久久久久影院| 欧美+日韩+精品| 在线观看美女被高潮喷水网站| 成年女人在线观看亚洲视频 | 听说在线观看完整版免费高清| 国产高清三级在线| 国产精品国产三级国产av玫瑰| av在线老鸭窝| 2018国产大陆天天弄谢| 成人美女网站在线观看视频| 少妇的逼好多水| 午夜免费男女啪啪视频观看| 久久综合国产亚洲精品| www.av在线官网国产| 亚洲婷婷狠狠爱综合网| 亚洲经典国产精华液单| 欧美日韩国产mv在线观看视频 | 亚洲精品自拍成人| 国产乱来视频区| eeuss影院久久| 尤物成人国产欧美一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 天堂网av新在线| 五月玫瑰六月丁香| 性色av一级| 成人毛片a级毛片在线播放| 麻豆精品久久久久久蜜桃| av国产免费在线观看| 直男gayav资源| av又黄又爽大尺度在线免费看| 日韩欧美一区视频在线观看 | 精品久久久精品久久久| 97超碰精品成人国产| 青春草亚洲视频在线观看| 欧美 日韩 精品 国产| 白带黄色成豆腐渣| 欧美zozozo另类| 国产高清不卡午夜福利| 免费观看的影片在线观看| 亚洲av日韩在线播放| 国产片特级美女逼逼视频| 国产午夜福利久久久久久| 国产精品嫩草影院av在线观看| 在线看a的网站| 久久久色成人| 久久热精品热| 久久久久国产精品人妻一区二区| 99热6这里只有精品| 午夜福利网站1000一区二区三区| 丝袜脚勾引网站| 两个人的视频大全免费| 美女cb高潮喷水在线观看| 身体一侧抽搐| 一级黄片播放器| 亚洲天堂av无毛| 狂野欧美激情性bbbbbb| 精品亚洲乱码少妇综合久久| 99热国产这里只有精品6| 内地一区二区视频在线| 国产一区亚洲一区在线观看| 久久精品熟女亚洲av麻豆精品| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲三级黄色毛片| 99热网站在线观看| 男女国产视频网站| a级一级毛片免费在线观看| 亚洲国产欧美人成| 精品少妇黑人巨大在线播放| 毛片一级片免费看久久久久| 在线天堂最新版资源| 国产 精品1| 国产极品天堂在线| 一级av片app| 亚洲伊人久久精品综合| 老女人水多毛片| 又粗又硬又长又爽又黄的视频| 成人午夜精彩视频在线观看| 我要看日韩黄色一级片| 日本爱情动作片www.在线观看| 高清视频免费观看一区二区| 欧美激情在线99| av福利片在线观看| 亚洲最大成人手机在线| av在线观看视频网站免费| 久久精品国产鲁丝片午夜精品| 久久精品国产自在天天线| 精品亚洲乱码少妇综合久久| 亚洲美女视频黄频| 日本三级黄在线观看| 一本色道久久久久久精品综合| 大又大粗又爽又黄少妇毛片口| 91久久精品国产一区二区三区| 人人妻人人爽人人添夜夜欢视频 | 亚洲av男天堂| tube8黄色片| 又大又黄又爽视频免费| 在线观看国产h片| 欧美老熟妇乱子伦牲交| 中文天堂在线官网| 激情 狠狠 欧美| 国产伦在线观看视频一区| 91在线精品国自产拍蜜月| 特级一级黄色大片| 亚洲高清免费不卡视频| 1000部很黄的大片| 日本黄色片子视频| 黄片wwwwww| 亚洲人成网站高清观看| 久久久精品欧美日韩精品| 国产成人freesex在线|