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

    基于衛(wèi)星和地面觀測的中國典型城市群對流層內(nèi)臭氧時空變化特征

    2022-03-24 08:36:48郭云飛包云軒劉端陽
    環(huán)境科學研究 2022年3期
    關鍵詞:對流層珠三角臭氧

    郭云飛,包云軒,劉端陽

    1. 南京信息工程大學氣象災害預報和評估協(xié)同創(chuàng)新中心,江蘇 南京 210044

    2. 江蘇省無錫學院物聯(lián)網(wǎng)設備超融合應用與安全工程研究中心,江蘇 無錫 214105

    3. 江蘇省氣象科學研究所,中國氣象局交通氣象重點實驗室,江蘇 南京 210008

    4. 南京氣象科技創(chuàng)新研究院,江蘇 南京 210008

    臭氧(O3)是大氣中重要的微量氣體,90%左右位于平流層,其對大氣化學和環(huán)流過程起重要作用[1].對流層O3占10%左右,高濃度O3不僅對空氣質量造成影響,還對人體健康和陸地生態(tài)系統(tǒng)產(chǎn)生危害[2-3].自2013年我國《大氣污染防治行動計劃》頒布以來,各地區(qū)PM2.5排放得到有效控制[4]. 然而,隨著很多城市O3前體物-NOx和VOCs排放的增多,O3污染日益凸顯,O3已逐漸取代顆粒物成為大氣中的首要污染物[5]. 因此,及時準確地掌握大范圍尺度下對流層和近地面O3的時空分布及其演變特征,對于解決O3污染問題、提高空氣質量具有重要意義.

    目前,許多學者對我國重大城市群(京津冀、長三角、珠三角城市群等)地面O3濃度的時空分布特征開展研究[6-9],結果表明,這些城市群均存在嚴重的O3污染問題并呈逐年加重趨勢. 但研究多僅局限于單點或地面基站的監(jiān)測,且地域代表性常受到制約.隨著空間遙感技術的發(fā)展,衛(wèi)星遙感為獲取全球或區(qū)域尺度的O3探測數(shù)據(jù)提供了可能,國內(nèi)學者常使用美國EOS-Aura衛(wèi)星搭載的OMI臭氧監(jiān)測儀數(shù)據(jù)和中國“風云三號”氣象衛(wèi)星搭載的紫外臭氧柱總量探測儀(TOU)數(shù)據(jù)對大氣O3柱含量變化展開討論[10-11].但多數(shù)研究偏重整個大氣柱O3柱總量的分析,因平流層O3含量占比較高且受緯度制約性較大,大氣O3柱總量也表現(xiàn)出隨緯度升高而增加的特征[12],很難體現(xiàn)出地形和人類排放因子對O3柱總量的貢獻情況,因此需進一步分析對流層內(nèi)O3柱總量的時空分布狀況;同時,已有研究多局限于某一地區(qū)或城市,針對對流層O3柱總量演變趨勢的研究較少.

    Ziemke等[13]通過殘差法反演得到的OMI-MLS(臭氧監(jiān)測儀-微波臨邊探測器)對流層O3柱總量數(shù)據(jù)在中國地區(qū)的適用性較高,可以揭示近地面的O3污染情況[14-16]. Ziemke等[17]研究發(fā)現(xiàn),北半球對流層O3柱總量在北美大陸和亞洲大陸東部的中緯度地區(qū)較高,而京津冀、長三角和珠三角城市群作為中國東部中低緯度地區(qū)較大的經(jīng)濟中心,O3污染問題受到學者的廣泛關注. 因此,該研究擬利用OMI-MLS衛(wèi)星數(shù)據(jù)來研究京津冀、長三角和珠三角城市群對流層O3柱總量的時空分布特征和演變趨勢,同時結合地面監(jiān)測的O3數(shù)據(jù),分析各城市近地面O3濃度的分布及污染狀況,以期為中國城市群環(huán)境治理和O3減排措施的制定提供科學依據(jù).

    1 資料與方法

    1.1 城市群劃分

    選取我國東部三大典型城市群開展研究,包括京津冀城市群13個城市、長三角城市群26個城市和珠三角城市群9個城市,具體各城市群城市名稱及空間分布如圖1所示.

    圖 1 中國典型城市群空間位置概況Fig.1 Overview of the spatial locations of typical urban agglomerations in China

    1.2 數(shù)據(jù)來源

    1.2.1對流層O3柱總量數(shù)據(jù)

    采用OMI-MLS對流層O3柱總量數(shù)據(jù)對城市群對流層O3柱總量進行對比分析,該數(shù)據(jù)集是由Ziemke等[13]通過殘差法反演得到,即以同時搭載在Aura衛(wèi)星上的臭氧監(jiān)測儀(OMI)探測到的大氣O3柱總量和微波臨邊探測器(MLS)探測的平流層以上O3廓線數(shù)據(jù)為基礎,通過計算二者之差可得出對流層各測點O3柱總量. 數(shù)據(jù)采用多步森單位DU (dobson unit,1 DU=2.69×1016cm-1),空間分辨率為1.00°×1.25°,由美國國家航空航天局(NASA,National Aeronautics and Space Administration in USA)戈達德太空飛行中心(Goddard Space Flight Center,https://acd-ext.gsfc.nasa.gov/Data_services/cloud_slice/new_data.html)提供,目前更新了2004年10月-2020年12月對流層O3柱濃度逐月數(shù)據(jù). 經(jīng)過過濾處理,剔除了多云、陰天等惡劣天氣下質量較差的數(shù)據(jù),但對于部分地區(qū)惡劣天氣出現(xiàn)頻率較大的月份,因數(shù)據(jù)過濾日數(shù)多,其代表性會受到一定限制. 整體上數(shù)據(jù)集經(jīng)過全球各國的驗證和采用,尤其在東亞地區(qū)準確性較高、可靠性較好[18-20]. 該研究選用的數(shù)據(jù)集時間范圍為2005年1月-2020年12月,空間覆蓋范圍為上述我國三大城市群.

    1.2.2近地面大氣層O3監(jiān)測數(shù)據(jù)

    京津冀、長三角和珠三角城市群地面O3監(jiān)測數(shù)據(jù)來自中國環(huán)境監(jiān)測總站的全國城市空氣質量實時發(fā)布平臺(http://106.37.208.233:20035),共包括京津冀13個城市(80個監(jiān)測站)、長三角26個城市(137個監(jiān)測站)、珠三角9個城市(56個監(jiān)測站)的O3逐小時濃度數(shù)據(jù),時間范圍為2015年1月-2020年12月,通過取城市所有監(jiān)測站點的O3濃度平均值來代表該城市的總體空氣質量狀況.

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

    1.3.1時空分布計算方法

    基于對流層O3柱總量數(shù)據(jù)可以計算三大城市群年均值、季均值和月均值,其中季節(jié)劃分標準為春季(3-5月)、夏季(6-8月)、秋季(9-11月)和冬季(12月-翌年2月).

    對于地面O3監(jiān)測數(shù)據(jù),根據(jù)《環(huán)境空氣質量標準》(GB 3095-2012)二級標準限值和《環(huán)境空氣質量評價技術規(guī)范(試行)》(HJ 663-2013)規(guī)定,計算出O3日最大8 h滑動平均濃度(簡稱“O3-8 h濃度”),以該年第90百分位數(shù)O3-8 h濃度限值(160 μg/m3)來判定O3是否年超標. 依據(jù)《環(huán)境空氣質量指數(shù)(AQI)技術規(guī)定(試行)》(HJ 633-2012)將O3-8 h濃度劃分為5個等級,分別為優(yōu)(1~100 μg/m3)、良(101~160 μg/m3)、輕度污染(161~215 μg/m3)、中度污染(216~265 μg/m3)、重度及以上污染(>265 μg/m3),對5個等級分別賦予不同的權重系數(shù)(0、1、2、3和5),以創(chuàng)建反映地區(qū)O3污染輕重的指數(shù)OPI,計算公式:

    式中:OPI為O3污染指數(shù);Di(1)、Di(2)、Di(3)、Di(4)、Di(5)為統(tǒng)計出的O3-8 h濃度5個等級的天數(shù),d;pi(1)、pi(2)、pi(3)、pi(4)、pi(5)分別為5個等級對應的權重系數(shù);n代表年份. OPI污染指數(shù)越高,表明O3污染越嚴重. 利用Jenks自然斷點法對結果分級,其原理是運用了聚類的思想,使組中的差異最小化,外部組和組之間的差異最大化,從而對相似值進行最有效的區(qū)分.

    1.3.2演變趨勢分析方法

    為了探究三大城市群對流層O3柱總量的長期變化趨勢,該文采用一元線性回歸模型[21]加以研究,計算公式:

    式中:t為年份序列;y為t年對流層O3柱總量的估計值,DU;α為對流層O3柱總量的年均變化量,DU;β為回歸模型在Y軸上的截距,DU. 若α>0,則對流層O3柱總量呈上升趨勢;若α<0,對流層O3柱總量為下降趨勢.α的絕對值可以反映年均變化的快慢,通過相關系數(shù)檢驗法對趨勢結果做顯著性驗證.

    2 結果與討論

    2.1 OMI-MLS對流層O3柱總量分析

    2.1.1多年整體變化狀況

    2005-2020年京津冀、長三角和珠三角城市群對流層O3柱總量的月變化情況如圖2所示. 由圖2可見:三大城市群之間對流層O3柱總量月變化有較大差異,其中在京津冀和長三角城市群表現(xiàn)為單峰特征,而在珠三角城市群呈雙峰的變化特征. 京津冀城市群對流層O3柱總量在1-6月穩(wěn)定增長,并于7月達到峰值(52.5 DU),之后迅速降低,10月后下降速率減緩并趨于穩(wěn)定;長三角城市群對流層O3柱總量上升速率顯著大于下降速率,峰值出現(xiàn)在5月(47.8 DU),隨后一直呈降低趨勢,7-9月下降速率得到一定減緩;珠三角城市群第一次峰值出現(xiàn)在4月(44.2 DU),5-7月迅速下降,并于8月達到谷值,雖然在9月出現(xiàn)了次峰,但峰值僅為第一次峰值的80%. 三大城市群之間對流層O3柱總量峰值表現(xiàn)為京津冀城市群>長三角城市群>珠三角城市群的特征,并且峰值出現(xiàn)時間隨著緯度的增加逐漸后延.

    圖 2 2005—2020年京津冀、長三角和珠三角城市群對流層O3柱總量月變化情況Fig.2 Monthly variation of tropospheric O3 column amount in the urban agglomerations of Beijing-Tianjin-Hebei, Yangtze River Delta and Pearl River Delta from 2005 to 2020

    三大城市群對流層O3柱總量在空間上同樣表現(xiàn)出顯著差異. 由圖3可見:京津冀城市群對流層O3柱總量最高,長三角城市群次之,珠三角城市群最低. 京津冀城市群對流層O3柱總量變化范圍為37.0~41.6 DU,呈由東南向西北減小的環(huán)帶狀分布特征,其中東南部的滄州市對流層O3柱總量最高;長三角城市群對流層O3柱總量變化范圍為36.9~40.8 DU,表現(xiàn)出隨緯度升高而增加的特征;珠三角城市群對流層O3柱總量變化范圍為35.2~36.1 DU,呈西北高于東南的分布特征. Li[22]研究發(fā)現(xiàn),對流層O3柱總量的空間分布受地形影響較大,低海拔地區(qū)往往較高. 京津冀西北部和長三角南部多為山地、丘陵地貌,海拔相對較高,對流層O3柱總量低于京津冀東南部和長三角北部的平原地區(qū). 地形對對流層O3柱總量空間分布的影響可歸結為大氣環(huán)流作用,低海拔地區(qū)易出現(xiàn)逆溫現(xiàn)象,造成O3、NO2等污染物的堆積,且日照時間較長,加快了光化學反應速率[14]. 長三角城市群對流層O3柱總量略低于京津冀城市群,這可能由于長三角城市群東面臨海,內(nèi)陸產(chǎn)生的O3等污染物會被海上大量的清潔空氣所稀釋[23]. 而處在華南的珠三角城市群雖多為平原,但平原面積較小,三面受到山地丘陵包圍,且受季風輸送與稀釋的影響較大,對流層O3柱總量較低,海拔對該地區(qū)對流層O3柱總量空間分布的影響較弱.

    圖 3 2005?2020年京津冀、長三角和珠三角城市群對流層O3柱總量空間分布及海拔Fig.3 Spatial distribution of tropospheric O3 column amount and altitude in the urban agglomerations of Beijing-Tianjin-Hebei, Yangtze River Delta and Pearl River Delta from 2005 to 2020

    2.1.2季節(jié)性變化特征

    進一步分析三大城市群不同季節(jié)對流層O3柱總量的空間分布發(fā)現(xiàn),城市群間、同一城市群不同季節(jié)的空間分布有顯著差異,空間異質性較大(見圖4).京津冀城市群四季對流層O3柱總量平均值分別為39.6、50.0、36.8和32.1 DU,均表現(xiàn)為從東南向西北逐漸減小的趨勢,變化梯度較為均勻. 夏季對流層O3柱總量最高,其中東南部最高,達52.4 DU;春季對流層O3柱總量明顯低于夏季,秋季低于春季,冬季最低. 長三角城市群四季對流層O3柱總量平均值分別為42.2、44.4、37.6和32.3 DU,空間分布季節(jié)性差異較大. 夏季對流層O3柱總量最高,冬季最低,而春季則略低于夏季,秋季又明顯低于春季;夏、秋兩季對流層O3柱總量均表現(xiàn)為南低北高的特征,其中夏季地區(qū)差異較明顯(南北相差10.6 DU). 珠三角城市群四季對流層O3柱總量平均值分別為42.2、34.4、34.4和32.0 DU,空間上呈現(xiàn)出較明顯的緯向梯度分布特征. 春季對流層O3柱總量明顯高于其余3季,且地區(qū)差異較?。幌?、秋兩季對流層O3柱總量均呈南低北高,且夏季略低于秋季的特征;冬季表現(xiàn)出與夏、秋兩季相反的空間分布格局.

    圖 4 2005—2020年京津冀、長三角和珠三角城市群對流層O3柱總量的季節(jié)性空間分布情況Fig.4 Seasonal spatial distribution of tropospheric O3 column amount in the urban agglomerations of Beijing-Tianjin-Hebei, Yangtze River Delta and Pearl River Delta from 2005 to 2020

    結合已有研究來分析三大城市群對流層O3柱總量季節(jié)性差異的原因. 京津冀和長三角城市群對流層O3柱總量的季節(jié)性分布受溫度影響較大,夏季溫度高,太陽輻射強,有利于光化學反應的進行;而春、秋兩季溫度相差較小,但對流層O3柱總量均在春季高于秋季,這可能與春季中低緯度地區(qū)較為頻繁的平流層-對流層O3交換有關[24]. 研究[25]表明,對流層頂折疊引發(fā)的平流層-對流層O3交換是春季對流層O3柱總量增加的重要來源. 與Wang等[26]發(fā)現(xiàn)廣西壯族自治區(qū)對流層O3柱總量在夏季出現(xiàn)極低值的現(xiàn)象類似,珠三角城市群對流層O3柱總量在夏季明顯降低,與亞夏季風影響下O3被稀釋并往內(nèi)陸遷移有關[27];另外,當相對濕度較高時,水汽所含的·OH、HO2·等自由基會迅速將O3分解成O2,從而降低O3濃度[28]. 珠三角城市群夏季溫度雖然較高,但頻繁出現(xiàn)的降水天氣會抑制O3的形成并促進其分解,并且降水天氣也對衛(wèi)星數(shù)據(jù)的反演造成一定影響. 在全球范圍內(nèi),對流層O3柱總量在春季、夏季較高[29],其中,中高緯度地區(qū)多以夏季最高[30],而中低緯度地區(qū)多以春季最高[31-32],與筆者研究的三大城市群對流層O3柱總量的季節(jié)性變化特征一致.

    2.1.3年際變化特征

    2005-2020年三大城市群對流層O3柱總量的年際變化情況見圖5. 由圖5可見,2005-2020年三大城市群對流層O3柱總量均呈逐年顯著增長趨勢(P<0.01),年均增加值分別為0.25、0.28和0.27 DU,長三角城市群增長最快,京津冀城市群最慢. 其中,京津冀城市群對流層O3柱總量的波動范圍為37.4~41.4 DU,2009-2015年波動較大,2016年以后變化平穩(wěn)但維持在一個較高的范圍,最高峰出現(xiàn)在2013年;長三角城市群波動范圍為36.8~41.8 DU,雖然在2014-2016年有所下降,但2016年后迅速上升,并在2018年和2019年出現(xiàn)峰值,2020年對流層O3柱總量出現(xiàn)突降,可能與新冠肺炎疫情期間O3前體物排放減少有關[33];珠三角城市群對流層O3柱總量波動范圍為33.0~39.0 DU,高低振蕩較為劇烈,2019年出現(xiàn)峰值,2020年同樣受新冠肺炎疫情影響出現(xiàn)了明顯的下降趨勢.

    圖 5 2005—2020年京津冀、長三角和珠三角城市群對流層O3柱總量的年際和四季變化情況Fig.5 Interannual and seasonal variations of tropospheric O3 column amount in the urban agglomerations of Beijing-Tianjin-Hebei, Yangtze River Delta and Pearl River Delta from 2005 to 2020

    進一步分析不同季節(jié)對流層O3柱總量變化趨勢,以確定各季節(jié)對對流層O3柱總量的影響. 京津冀城市群四季對流層O3柱總量年均增加值分別為0.22、0.19、0.29、0.29 DU,夏季對流層O3柱總量波動最小且增長最慢,秋、冬兩季波動較大且增長最快. 長三角城市群四季對流層O3柱總量年均增加值分別為0.24、0.39、0.32和0.18 DU,在對流層O3柱總量最高的夏季增長最快,秋季次之. 珠三角城市群四季對流層O3柱總量年均增加值分別為0.39、0.30、0.29和0.16 DU,春季對流層O3柱總量最高且增長最快,夏、秋兩季的變化趨勢相似,冬季增長最慢. 綜上,京津冀城市群在對流層O3柱總量較低的秋冬季年均增長最快,而長三角和珠三角城市群分別在對流層O3柱總量最高的夏季和春季增長最快. Zhang等[34]研究發(fā)現(xiàn),全球范圍內(nèi)對流層O3柱總量增長最快的地區(qū)從地中海東部沿岸延伸到印度和東南亞,再向東越過太平洋地區(qū). Ziemke等[35]發(fā)現(xiàn),近38年來東亞地區(qū)對流層O3柱總量年均增長由2005年之前的0.12~0.14 DU增至0.3 DU,與筆者研究結論基本一致.

    2.2 近地面O3濃度的分析

    2.2.1O3-8 h濃度與對流層O3柱總量的對比

    2015-2020年三大城市群對流層O3柱總量與O3-8 h濃度對比見圖6. 由圖6可見,京津冀和長三角城市群對流層O3柱總量與O3-8 h濃度相關性較好(P<0.01). 2015-2020年京津冀城市群對流層O3柱總量均呈倒“V”型,O3-8 h濃度在多數(shù)年份也呈現(xiàn)此特征. 長三角城市群對流層O3柱總量與O3-8 h濃度在多數(shù)年份呈“M”型,個別年份表現(xiàn)為單峰型. 而珠三角城市群對流層O3柱總量與O3-8 h濃度相關性不明顯,對流層O3柱總量在每年4月達到峰值,該時段O3-8 h濃度較低;而秋季因受來自大陸上風向的區(qū)域傳輸及局地污染導致O3聚集[20],珠三角城市群O3-8 h濃度出現(xiàn)峰值,對流層O3柱總量則維持在較低水平. 珠三角城市群對流層O3柱總量與O3-8 h濃度相關性較差,一方面與對流層O3柱總量受平流層-對流層O3交換過程及季風環(huán)流擾動關系密切;另一方面,不可控的環(huán)流形勢和氣象因素會造成不同年份O3-8 h濃度的月變化出現(xiàn)差異,并且會對反演的衛(wèi)星數(shù)據(jù)準確性產(chǎn)生限制,從而影響二者相關性.

    圖 6 2015?2020年京津冀、長三角和珠三角城市群對流層O3柱總量與O3-8 h濃度的對比Fig.6 Comparison of tropospheric O3 column amount and O3-8 h concentration in the urban agglomerations of Beijing-Tianjin-Hebei, Yangtze River Delta and Pearl River Delta from 2015 to 2020

    2.2.2O3-8 h濃度和第90百分位O3-8 h濃度的空間分布情況

    利用Jenks自然斷點法對2015-2020年三大城市群所有城市O3-8 h濃度及第90百分位O3-8 h濃度進行空間分類(見圖7). 由圖7可見,各城市群不同年份O3-8 h濃度的空間分布差異明顯,第90百分位O3-8 h濃度的空間分布與O3-8 h濃度空間分布成正比.

    圖 7 2015—2020年京津冀、長三角和珠三角城市群O3-8 h濃度和第90百分位O3-8 h濃度的空間分布情況Fig.7 Spatial distribution of O3-8 h concentration and the 90th percentile O3-8 h concentration in the urban agglomerations of Beijing-Tianjin-Hebei, Yangtze River Delta and Pearl River Delta from 2015 to 2020

    京津冀城市群O3-8 h濃度范圍為64.9~119.6 μg/m3,第90百分位O3-8 h濃度范圍為108~216 μg/m3,空間分布上由2016年之前的北高南低轉變?yōu)槟细弑钡? 2015年和2016年O3-8 h濃度較低,這兩年O3污染程度最輕. 2017-2019年O3-8 h濃度較高,其中2018年濃度最高(區(qū)域平均值為110.9 μg/m3),2017年次之(區(qū)域平均值為107.8 μg/m3),且2017-2019年所有城市O3均未達標,污染較為嚴重,與王玫等[36]研究結論一致. 2017年保定市、滄州市、廊坊市和邢臺市O3-8 h濃度均超過110 μg/m3,這4個城市與唐山市的第90百分位O3-8 h濃度突破了200 μg/m3.2018年共8個城市O3-8 h濃度超過110 μg/m3,其中保定市、石家莊市、滄州市、天津市和邢臺市第90百分位O3-8 h濃度均達到200 μg/m3. 研究[37]表明,顆粒物的非均相反應會降低O3濃度,減排背景下顆粒物濃度的下降對O3的靜作用可能導致了O3-8 h濃度的升高. 2019年O3污染有所改善,O3-8 h濃度較高的區(qū)域主要集中在京津冀城市群的南部城市. 2020年O3-8 h濃度明顯下降,該年張家口市和承德市O3年達標,污染得到緩解.

    長三角城市群O3-8 h濃度范圍為39.0~118.8 μg/m3,第90百分位O3-8 h濃度范圍為61~201 μg/m3,其中2015年和2016年呈由沿海向內(nèi)陸減小的布局,2017-2020年則表現(xiàn)為北高南低. 2015年和2016年O3-8 h濃度相對較低,其中2015年O3-8 h濃度以塊狀分布為主,年超標城市均位于東部沿海地區(qū);2016年東北部的城市O3-8 h濃度有所降低,而西部安徽省部分城市O3-8 h濃度增加迅速;2017年O3-8 h濃度為近6年最高(區(qū)域平均值達106.7 μg/m3),O3-8 h高濃度中心開始西移,西北部的滁州市O3-8 h濃度最高(118.8 μg/m3),此外中部以北的湖州市、南通市、上海市、嘉興市、揚州市、泰州市和鎮(zhèn)江市的O3-8 h濃度也均超過110 μg/m3,而第90百分位O3-8 h濃度揚州市最高,2017年僅不到1/3的城市O3年達標. 常爐予等[38]模擬發(fā)現(xiàn):2017年受西太平洋副高控制的環(huán)流形勢更有利于O3污染的形成,這一年的副高較往年偏強且西伸脊點更西,這可能導致了長三角城市群O3-8 h濃度較高;2018年和2019年O3-8 h高濃度中心進一步西移,O3-8 h濃度呈降低趨勢,東部沿海城市O3污染得到明顯改善,而西部安慶市、池州市等城市則由O3達標轉變?yōu)槌瑯耍?020年因新冠肺炎疫情導致O3-8 h濃度進一步降低,達標城市明顯增多,O3污染得到有效改善. 長三角城市群O3-8 h高濃度中心的西移,這是由于沿海城市O3前體物的排放減少,而西部安徽省部分城市經(jīng)濟的迅速發(fā)展造成O3-8 h濃度升高[39].

    珠三角城市群O3-8 h濃度范圍為72.4~112.2 μg/m3,第90百分位O3-8 h濃度范圍為127~197 μg/m3,其中2015-2019年O3-8 h濃度呈逐年增加的趨勢,2020年O3污染有所改善. 2015年O3-8 h和第90百分位O3-8 h濃度均為近6年最低,僅東莞市O3年超標.2016-2019年O3-8 h濃度同比均有所上升,增速分別為3.8%、10.9%、0.6%和3.8%,其中位于中部的江門市、中山市和佛山市增長迅猛,至2019年,珠三角O3-8 h濃度區(qū)域平均值達100.4 μg/m3;并且O3年超標城市數(shù)量也逐年增加,至2019年,僅有惠州市和深圳市O3年達標. 珠三角城市群O3-8 h濃度在中部地區(qū)上升速率明顯高于外圍地區(qū),原因為中部地區(qū)NO2濃度的顯著下降使得O3滴定效應減弱[40]. 2020年所有城市O3-8 h濃度有所下降,其中深圳市和肇慶市的O3-8 h濃度最低,該年僅江門市未達標.

    2.2.3污染等級空間分布

    按式(1)計算出各城市O3污染指數(shù)OPI,其變化范圍為0.34~0.77,并利用Jenks自然斷點法將污染分為很輕、較輕、中等、較重和嚴重5個等級(見圖8).由圖8可見:京津冀城市群O3污染指數(shù)最高,且呈現(xiàn)出南高北低的分布狀況,其中,保定市、廊坊市、滄州市、衡水市、邢臺市和邯鄲市O3污染較嚴重,秦皇島市污染相對較輕;長三角城市群的O3污染指數(shù)空間分布表現(xiàn)為中心城市向外圍擴散,城市間差異較大的特征,其中,揚州市、湖州市和嘉興市O3污染較重,其次為中北部的鹽城市、泰州市、南京市、鎮(zhèn)江市、常州市和無錫市,而中西部的銅陵市和宣城市污染指數(shù)較低;珠三角城市群整體上O3污染程度較輕,其中惠州市和深圳市的O3污染指數(shù)較低.

    圖 8 2015?2020年京津冀、長三角和珠三角城市群O3污染指數(shù)的空間分布情況Fig.8 Spatial distribution of O3 pollution index in the urban agglomerations of Beijing-Tianjin-Hebei, Yangtze River Delta and Pearl River Delta from 2015 to 2020

    3 結論

    a) 基于OMI-MLS的對流層O3柱總量月峰值和年均值均呈京津冀城市群>長三角城市群>珠三角城市群的特征,空間分布受海拔影響較大,低海拔地區(qū)往往對流層O3柱總量較高. 季節(jié)分布上,京津冀城市群對流層O3柱總量在夏季(50.0 DU)最高,四季均表現(xiàn)為從西北向東南增大的趨勢;長三角城市群對流層O3柱總量同樣在夏季(44.4 DU)最高,表現(xiàn)出由低緯向高緯度增大的緯向帶狀分布特征;珠三角城市群對流層O3柱總量在春季(42.2 DU)明顯高于其他三季,空間上地區(qū)差異較小.

    b) 三大城市群對流層O3柱總量均呈現(xiàn)出逐年顯著增大的趨勢,年均增加值分別為0.25、0.28和0.27 DU. 所有季節(jié)對流層O3柱總量均表現(xiàn)出增長趨勢,京津冀城市群在對流層O3柱總量較低的秋冬季年均增長(0.29 DU)最快,而長三角和珠三角城市群分別在對流層O3柱總量最高的夏季和春季增長最快.

    c) 京津冀和長三角城市群對流層O3柱總量與O3-8 h濃度月均值均有良好的相關性,而珠三角城市群二者相關性較差.

    d) O3-8 h濃度表現(xiàn)為京津冀城市群>長三角城市群>珠三角城市群的特征,其中京津冀城市群在2018年O3-8 h濃度(110.9 μg/m3)最高,空間分布上由2016年之前的北高南低轉變?yōu)槟细弑钡?,多?shù)城市O3污染較重且年達標率較低;長三角城市群O3-8 h濃度在2017年(106.7 μg/m3)最高,2016年起O3-8 h高濃度中心由東北逐漸往西南部內(nèi)陸轉移,東南沿海城市O3年達標率顯著提高;珠三角城市群O3污染程度最輕,達標城市較多,但O3-8 h濃度逐年上升并于2019年達最高值(100.4 μg/m3),中心城市上升速率遠大于外圍城市.

    參考文獻(References):

    [1]李剛,譚言科,李崇銀,等.近30年北半球冬季臭氧總量分布特征及其與平流層溫度的關系[J].地球物理學報,2015,58(5):1475-1491.

    LI G,TAN Y K,LI C Y,et al.The distribution characteristics of total ozone and its relationship with stratospheric temperature during boreal winter in the recent 30 years[J].Chinese Journal of Geophysics,2015,58(5):1475-1491.

    [2]EBI K L,MCGREGOR G.Climate change,tropospheric ozone and particulate matter,and health impacts[J].Environmental Health Perspectives,2008,116(11):1449-1455.

    [3]REN W,TIAN H Q,TAO B,et al.Impacts of tropospheric ozone and climate change on net primary productivity and net carbon exchange of China′s forest ecosystems[J].Global Ecology and Biogeography,2011,20(3):391-406.

    [4]李紅,彭良,畢方,等.我國PM2.5與臭氧污染協(xié)同控制策略研究[J].環(huán)境科學研究,2019,32(10):1763-1778.

    LI H,PENG L,BI F,et al.Strategy of coordinated control of PM2.5and ozone in China[J].Research of Environmental Sciences,2019,32(10):1763-1778.

    [5]WANG T,XUE L K,BRIMBLECOMBE P,et al.Ozone pollution in China:a review of concentrations,meteorological influences,chemical precursors,and effects[J].Science of the Total Environment,2017,575:1582-1596.

    [6]李娜,周濤,劉小雪,等.廊坊市臭氧污染特征及其與氣象因素的關系[J].環(huán)境工程技術學報,2021,11(2):217-225.

    LI N,ZHOU T,LIU X X,et al.Characteristics of ozone pollution and its relationship with meteorological factors in Langfang City[J].Journal of Environmental Engineering Technology,2021,11(2):217-225.

    [7]栗澤苑,楊雷峰,華道柱,等.2013—2018年中國近地面臭氧濃度空間分布特征及其與氣象因子的關系[J].環(huán)境科學研究,2021,34(9):2094-2104.

    LI Z Y,YANG L F,HUA D Z,et al.Spatial pattern of surface ozone and its relationship with meteorological variables in China during 2013-2018[J].Research of Environmental Sciences,2021,34(9):2094-2104.

    [8]ZHANG A,LIN J H,CHEN W H,et al.Spatial-temporal distribution variation of ground-level ozone in China′s Pearl River Delta metropolitan region[J].International Journal of Environmental Research and Public Health,2021,18(3):872.

    [9]TANG G Q,ZHU X W,XIN J Y,et al.Modelling study of boundarylayer ozone over Northern China.Part I:ozone budget in summer[J].Atmospheric Research,2017,187:128-137.

    [10]王維和,張興贏,安興琴,等.風云三號氣象衛(wèi)星全球臭氧總量反演和真實性檢驗結果分析[J].科學通報,2010,55(17):1726-1733.

    [11]單源源,李莉,劉瓊,等.基于OMI數(shù)據(jù)的中國中東部臭氧及前體物的時空分布[J].環(huán)境科學研究,2016,29(8):1128-1136.

    SHAN Y Y,LI L,LIU Q,et al.Spatial-temporal distribution of ozone and its precursors over central and eastern China based on OMI data[J].Research of Environmental Sciences,2016,29(8):1128-1136.

    [12]于瑞新,劉旻霞,李亮,等.長三角地區(qū)近15年大氣臭氧柱濃度時空變化及影響因素[J].環(huán)境科學學報,2021,41(3):770-784.

    YU R X,LIU M X,LI L,et al.Spatial and temporal variation of atmospheric ozone column concentration and influencing factors in the Yangtze River Delta Region in recent 15 years[J].Acta Scientiae Circumstantiae,2021,41(3):770-784.

    [13]ZIEMKE J R,CHANDRA S,DUNCAN B N,et al.Tropospheric ozone determined from Aura OMI and MLS:evaluation of measurements and comparison with the Global Modeling Initiative′s Chemical Transport Model[J].Journal of Geophysical Research Atmospheres,2006,111(D19):D19303.

    [14]曹曉云,祁棟林,肖建設,等.基于衛(wèi)星觀測的青海高原對流層臭氧時空分布特征研究[J].環(huán)境科學學報,2021,41(5):1640-1648.

    CAO X Y,QI D L,XIAO J S,et al.Temporal and spatial distribution characteristics of tropospheric ozone in Qinghai Plateau based on satellite observations[J].Acta Scientiae Circumstantiae,2021,41(5):1640-1648.

    [15]BEIG G,SINGH V.Trends in tropical tropospheric column ozone from satellite data and MOZART model[J].Geophysical Research Letters,2007,34(17):L17801.

    [16]NOREEN A,KHOKHAR M F,ZEB N,et al.Spatio-temporal assessment and seasonal variation of tropospheric ozone in Pakistan during the last decade[J].Environmental Science and Pollution Research,2018,25(9):8441-8454.

    [17]ZIEMKE J R,CHANDRA S,LABOW G J,et al.A global climatology of tropospheric and stratospheric ozone derived from Aura OMI and MLS measurements[J].Atmospheric Chemistry and Physics,2011,11(17):9237-9251.

    [18]SCHOEBERL M R,ZIEMKE J R,BOJKOV B,et al.A trajectorybased estimate of the tropospheric ozone column using the residual method[J].Journal of Geophysical Research Atmospheres,2007,112(D24):D24S49.

    [19]NASSAR R,LOGAN J A,WORDEN H M,et al.Validation of Tropospheric Emission Spectrometer (TES) nadir ozone profiles using ozonesonde measurements[J].Journal of Geophysical Research Atmospheres,2008,113(D15):D15S17.

    [20]張倩倩,張興贏.基于衛(wèi)星和地面觀測的2013年以來我國臭氧時空分布及變化特征[J].環(huán)境科學,2019,40(3):1132-1142.

    ZHANG Q Q,ZHANG X Y.Ozone spatial-temporal distribution and trend over China since 2013:insight from satellite and surface observation[J].Environmental Science,2019,40(3):1132-1142.

    [21]魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預測技術[M].北京:氣象出版社,1999.

    [22]LI J.Pollution trends in China from 2000 to 2017:a multi-sensor view from space[J].Remote Sensing,2020,12(2):208.

    [23]徐曉斌,林偉立.衛(wèi)星觀測的中國地區(qū)1979—2005年對流層臭氧變化趨勢[J].氣候變化研究進展,2010,6(2):100-105.

    XU X B,LIN W L.1979-2005 trends of tropospheric ozone over China based on the satellite data[J].Advances in Climate Change Research,2010,6(2):100-105.

    [24]王飛,劉毅,蔡兆男,等.利用GOME-2衛(wèi)星數(shù)據(jù)反演對流層臭氧[J].遙感技術與應用,2016,31(2):316-323.

    WANG F,LIU Y,CAI Z N,et al.Tropospheric ozone products retrieved from GOME-2[J].Remote Sensing Technology and Application,2016,31(2):316-323.

    [25]COOPER O R,MOODY J L.Meteorological controls on ozone at an elevated eastern United States regional background monitoring site[J].Journal of Geophysical Research:Atmospheres,2000,105(D5):6855-6869.

    [26]WANG Y P,YU C,TAO J H,et al.Spatio-temporal characteristics of tropospheric ozone and its precursors in Guangxi,South China[J].Atmosphere,2018,9(9):355.

    [27]LIU Y,WANG Y,LIU X,et al.Tibetan middle tropospheric ozone minimum in June discovered from GOME observations[J].Geophysical Research Letters,2009,36(5):L05814.

    [28]齊冰,牛彧文,杜榮光,等.杭州市近地面大氣臭氧濃度變化特征分析[J].中國環(huán)境科學,2017,37(2):443-451.

    QI B,NIU Y W,DU R G,et al.Characteristics of surface ozone concentration in urban site of Hangzhou[J].China Environmental Science,2017,37(2):443-451.

    [29]COOPER O R,PARRISH D D,ZIEMKE J,et al.Global distribution and trends of tropospheric ozone:an observation-based review[J].Elementa:Science of the Anthropocene,2014.doi:10.12952/journal.elementa.000029.

    [30]OSTERMAN G B,KULAWIK S S,WORDEN H M,et al.Validation of Tropospheric Emission Spectrometer (TES)measurements of the total,stratospheric,and tropospheric column abundance of ozone[J].Journal of Geophysical Research Atmospheres,2008,113(D15):D15S16.

    [31]CHAN C Y,CHAN L Y,CUI H,et al.Origin of the springtime tropospheric ozone maximum over East China at Linan in 2001[J].Tellus B:Chemical and Physical Meteorology,2011,55(5):982-992.

    [32]JIA J,LADST?TTER-WEIENMAYER A,HOU X W,et al.Tropospheric ozone maxima observed over the Arabian Sea during the pre-monsoon[J].Atmospheric Chemistry and Physics,2017,17(8):4915-4930.

    [33]曠雅瓊,鄒忠,張秀英,等.長三角地區(qū)大氣污染物對新冠肺炎封城的時空響應特征[J].環(huán)境科學學報,2021,41(4):1165-1172.

    KUANG Y Q,ZOU Z,ZHANG X Y,et al.Spatial-temporal responses of atmospheric pollutants to the COVID-19 lockdown across the Yangtze River Delta Region[J].Acta Scientiae Circumstantiae,2021,41(4):1165-1172.

    [34]ZHANG Y Q,COOPER O R,GAUDEL A,et al.Tropospheric ozone change from 1980 to 2010 dominated by equatorward redistribution of emissions[J].Nature Geoscience,2016,9(12):875-879.

    [35]ZIEMKE J R,OMAN L D,STRODE S A,et al.Trends in global tropospheric ozone inferred from a composite record of TOMS/OMI/MLS/OMPS satellite measurements and the MERRA-2 GMI simulation[J].Atmospheric Chemistry and Physics,2019,19(5):3257-3269.

    [36]王玫,鄭有飛,柳艷菊,等.京津冀臭氧變化特征及與氣象要素的關系[J].中國環(huán)境科學,2019,39(7):2689-2698.

    WANG M,ZHENG Y F,LIU Y J,et al.Characteristics of ozone and its relationship with meteorological factors in Beijing-Tianjin-Hebei Region[J].China Environmental Science,2019,39(7):2689-2698.

    [37]余益軍,孟曉艷,王振,等.京津冀地區(qū)城市臭氧污染趨勢及原因探討[J].環(huán)境科學,2020,41(1):106-114.

    YU Y J,MENG X Y,WANG Z,et al.Driving factors of the significant increase in surface ozone in the Beijing-Tianjin-Hebei Region,China,during 2013-2018[J].Environmental Science,2020,41(1):106-114.

    [38]常爐予,許建明,瞿元昊,等.上海市臭氧污染的大氣環(huán)流客觀分型研究[J].環(huán)境科學學報,2019,39(1):169-179.

    CHANG L Y,XU J M,QU Y H,et al.Study on objective synoptic classification on ozone pollution in Shanghai[J].Acta Scientiae Circumstantiae,2019,39(1):169-179.

    [39]黃小剛,邵天杰,趙景波,等.長三角城市群臭氧濃度的時空分異及驅動因素[J].長江流域資源與環(huán)境,2019,28(6):1434-1445.

    HUANG X G,SHAO T J,ZHAO J B,et al.Spatio-temporal differentiation of ozone concentration and its driving factors in Yangtze River Delta urban agglomeration[J].Resources and Environment in the Yangtze Basin,2019,28(6):1434-1445.

    [40]趙偉,高博,盧清,等.2006—2019年珠三角地區(qū)臭氧污染趨勢[J].環(huán)境科學,2021,42(1):97-105.

    ZHAO W,GAO B,LU Q,et al.Ozone pollution trend in the Pearl River Delta region during 2006-2019[J].Environmental Science,2021,42(1):97-105.

    猜你喜歡
    對流層珠三角臭氧
    養(yǎng)殖面積逐年減少!珠三角這條魚曾風靡一時,為何現(xiàn)在漁民都不愿養(yǎng)殖它?
    文印室內(nèi)臭氧散發(fā)實測分析及模擬
    郴州地區(qū)對流層頂氣候概況
    制造業(yè)“逃離”珠三角?
    商周刊(2017年17期)2017-09-08 13:08:58
    珠三角世界級城市群建設進入施工期
    看不見的污染源——臭氧
    利用臭氧水防治韭菜遲眼蕈蚊
    中國蔬菜(2016年8期)2017-01-15 14:23:43
    實時干涉測量中對流層延遲與鐘差精修正建模
    載人航天(2016年4期)2016-12-01 06:56:24
    成都地區(qū)2005~2015年對流層NO2柱濃度趨勢與時空分布
    臭氧分子如是說
    丁香六月欧美| 黄频高清免费视频| 国产精品自产拍在线观看55亚洲 | 成人亚洲精品一区在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 老司机靠b影院| 亚洲欧美激情在线| 丝袜脚勾引网站| 中文字幕精品免费在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲九九香蕉| 精品人妻熟女毛片av久久网站| 老司机影院毛片| videosex国产| 十八禁网站免费在线| 日韩熟女老妇一区二区性免费视频| 午夜福利在线免费观看网站| 69av精品久久久久久 | 国产亚洲精品久久久久5区| 国产成人精品无人区| 亚洲精品自拍成人| 如日韩欧美国产精品一区二区三区| 热99re8久久精品国产| 欧美黄色淫秽网站| 国产精品免费大片| 国产av国产精品国产| netflix在线观看网站| 老熟妇乱子伦视频在线观看 | 1024香蕉在线观看| 两个人免费观看高清视频| 国产精品.久久久| 国产免费现黄频在线看| 天天添夜夜摸| 视频在线观看一区二区三区| 欧美精品一区二区免费开放| 超色免费av| 精品熟女少妇八av免费久了| 日韩一卡2卡3卡4卡2021年| 人妻久久中文字幕网| 91老司机精品| 中文字幕精品免费在线观看视频| 欧美精品亚洲一区二区| 老司机靠b影院| 亚洲精品一区蜜桃| 亚洲专区中文字幕在线| 97在线人人人人妻| 国产xxxxx性猛交| av视频免费观看在线观看| 久久精品国产亚洲av香蕉五月 | 高清av免费在线| 国产片内射在线| 高清在线国产一区| 精品国产乱子伦一区二区三区 | 肉色欧美久久久久久久蜜桃| 亚洲五月婷婷丁香| 黄片播放在线免费| 久久久久久人人人人人| 日韩视频一区二区在线观看| 狂野欧美激情性xxxx| 亚洲专区字幕在线| 国产日韩欧美亚洲二区| 一级片免费观看大全| 熟女少妇亚洲综合色aaa.| 免费日韩欧美在线观看| 久久天堂一区二区三区四区| 国产色视频综合| 久久热在线av| 18禁观看日本| 狠狠婷婷综合久久久久久88av| 操美女的视频在线观看| 免费日韩欧美在线观看| 老熟妇仑乱视频hdxx| 国产精品香港三级国产av潘金莲| 69av精品久久久久久 | 欧美日韩av久久| 国产真人三级小视频在线观看| 欧美在线一区亚洲| 中国国产av一级| 亚洲熟女毛片儿| a级片在线免费高清观看视频| 亚洲激情五月婷婷啪啪| 国产精品一区二区在线不卡| 久久久久久人人人人人| 免费观看av网站的网址| 成年女人毛片免费观看观看9 | 2018国产大陆天天弄谢| 黑人巨大精品欧美一区二区mp4| 人妻 亚洲 视频| 中亚洲国语对白在线视频| www.熟女人妻精品国产| 黑人巨大精品欧美一区二区mp4| 黄片小视频在线播放| 在线观看免费日韩欧美大片| 黄片播放在线免费| av在线app专区| 香蕉丝袜av| 久久久国产一区二区| 下体分泌物呈黄色| 在线天堂中文资源库| 欧美日韩福利视频一区二区| 99久久综合免费| 日本五十路高清| 亚洲九九香蕉| 久久午夜综合久久蜜桃| 丝袜在线中文字幕| 国内毛片毛片毛片毛片毛片| 欧美另类一区| 黄色 视频免费看| 亚洲精品粉嫩美女一区| 亚洲av男天堂| 精品福利观看| 黄色视频不卡| 老司机午夜十八禁免费视频| 一本大道久久a久久精品| av在线老鸭窝| 纵有疾风起免费观看全集完整版| 亚洲精品国产色婷婷电影| 一级a爱视频在线免费观看| 国产免费av片在线观看野外av| 精品福利永久在线观看| 丰满迷人的少妇在线观看| 午夜视频精品福利| 永久免费av网站大全| 国产片内射在线| 国产一级毛片在线| 亚洲七黄色美女视频| 国产亚洲午夜精品一区二区久久| 韩国高清视频一区二区三区| 国产亚洲精品久久久久5区| 夫妻午夜视频| 婷婷成人精品国产| 国产精品一区二区在线观看99| 亚洲免费av在线视频| www日本在线高清视频| 国产av一区二区精品久久| tocl精华| 一本大道久久a久久精品| 日韩欧美一区视频在线观看| 欧美日韩福利视频一区二区| 水蜜桃什么品种好| 人人妻人人澡人人爽人人夜夜| 成人国产av品久久久| 国产亚洲av高清不卡| 欧美黑人欧美精品刺激| 在线av久久热| 国产亚洲一区二区精品| 中文字幕另类日韩欧美亚洲嫩草| 久久久久精品国产欧美久久久 | 男女边摸边吃奶| 久久久久久久久免费视频了| 国产野战对白在线观看| 美女大奶头黄色视频| 亚洲av电影在线观看一区二区三区| 婷婷色av中文字幕| 亚洲精品中文字幕一二三四区 | 丰满饥渴人妻一区二区三| 日韩熟女老妇一区二区性免费视频| 一本—道久久a久久精品蜜桃钙片| 亚洲专区字幕在线| 十八禁网站网址无遮挡| 美女扒开内裤让男人捅视频| 亚洲精品国产av蜜桃| 天堂中文最新版在线下载| 一级毛片电影观看| 99国产极品粉嫩在线观看| 在线看a的网站| 搡老乐熟女国产| 午夜视频精品福利| 男女边摸边吃奶| 丁香六月欧美| 美女中出高潮动态图| 黄片大片在线免费观看| 日韩免费高清中文字幕av| 午夜免费观看性视频| 麻豆国产av国片精品| 精品国产国语对白av| 亚洲国产看品久久| 亚洲色图 男人天堂 中文字幕| 亚洲av片天天在线观看| 美女高潮喷水抽搐中文字幕| 2018国产大陆天天弄谢| 12—13女人毛片做爰片一| 国产有黄有色有爽视频| 午夜免费鲁丝| 久久精品熟女亚洲av麻豆精品| 国产黄频视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩制服骚丝袜av| 亚洲熟女毛片儿| 老司机影院毛片| 丰满人妻熟妇乱又伦精品不卡| 婷婷成人精品国产| 中亚洲国语对白在线视频| 少妇猛男粗大的猛烈进出视频| 午夜影院在线不卡| 成年人黄色毛片网站| 9191精品国产免费久久| 一二三四社区在线视频社区8| 狂野欧美激情性xxxx| 不卡av一区二区三区| 97在线人人人人妻| 成人亚洲精品一区在线观看| 69精品国产乱码久久久| 国产精品 欧美亚洲| 搡老熟女国产l中国老女人| 午夜影院在线不卡| 99精品久久久久人妻精品| www日本在线高清视频| 9191精品国产免费久久| 80岁老熟妇乱子伦牲交| 国产97色在线日韩免费| 秋霞在线观看毛片| 亚洲欧美精品自产自拍| 在线观看免费高清a一片| 无遮挡黄片免费观看| 十八禁人妻一区二区| 国产成人精品久久二区二区免费| 亚洲第一青青草原| av网站在线播放免费| 美女视频免费永久观看网站| 国产精品 国内视频| 国产成人精品久久二区二区免费| 99re6热这里在线精品视频| 在线av久久热| 国产深夜福利视频在线观看| 久久性视频一级片| 亚洲国产精品成人久久小说| 欧美日韩亚洲综合一区二区三区_| 男女国产视频网站| 99国产极品粉嫩在线观看| 最新的欧美精品一区二区| 亚洲国产欧美网| 久久精品国产亚洲av香蕉五月 | 香蕉国产在线看| 老司机深夜福利视频在线观看 | 一本久久精品| 国产欧美亚洲国产| 久久影院123| 女人高潮潮喷娇喘18禁视频| 成人黄色视频免费在线看| 国产老妇伦熟女老妇高清| 国产亚洲欧美在线一区二区| 亚洲欧美日韩高清在线视频 | 最近最新免费中文字幕在线| 高清欧美精品videossex| 亚洲精品美女久久av网站| 人妻久久中文字幕网| 久久人人97超碰香蕉20202| 日韩人妻精品一区2区三区| 美女福利国产在线| 国产精品二区激情视频| 男女高潮啪啪啪动态图| 国产精品香港三级国产av潘金莲| 亚洲国产av新网站| 久久久久网色| 悠悠久久av| 午夜91福利影院| 欧美 日韩 精品 国产| 亚洲av日韩精品久久久久久密| 国产亚洲午夜精品一区二区久久| 天天影视国产精品| 久久精品久久久久久噜噜老黄| 五月天丁香电影| 欧美变态另类bdsm刘玥| 一级,二级,三级黄色视频| 免费av中文字幕在线| 韩国高清视频一区二区三区| 亚洲男人天堂网一区| 国产欧美日韩一区二区精品| 日韩一区二区三区影片| 精品福利永久在线观看| www日本在线高清视频| 国产免费现黄频在线看| 日日摸夜夜添夜夜添小说| 80岁老熟妇乱子伦牲交| 久久久国产一区二区| 午夜免费成人在线视频| 免费在线观看黄色视频的| 少妇 在线观看| 国内毛片毛片毛片毛片毛片| 国产成人a∨麻豆精品| 亚洲国产精品999| 亚洲综合色网址| 精品亚洲乱码少妇综合久久| 国产成人一区二区三区免费视频网站| 亚洲熟女毛片儿| 交换朋友夫妻互换小说| 欧美 亚洲 国产 日韩一| 欧美激情极品国产一区二区三区| 天堂8中文在线网| 国产淫语在线视频| av网站在线播放免费| 新久久久久国产一级毛片| av福利片在线| 久久久久久人人人人人| 国产一区二区三区av在线| 丝袜脚勾引网站| 免费女性裸体啪啪无遮挡网站| 91精品国产国语对白视频| 男女床上黄色一级片免费看| 久久影院123| 一区二区日韩欧美中文字幕| 亚洲精品乱久久久久久| tube8黄色片| 亚洲伊人久久精品综合| 午夜免费成人在线视频| 日本vs欧美在线观看视频| 久久久久久人人人人人| 亚洲人成电影免费在线| 纯流量卡能插随身wifi吗| 最新的欧美精品一区二区| 老熟妇乱子伦视频在线观看 | 久久久久视频综合| 精品亚洲乱码少妇综合久久| 性少妇av在线| 国产在线一区二区三区精| 亚洲人成电影观看| 亚洲欧美激情在线| 国产精品偷伦视频观看了| 手机成人av网站| 少妇裸体淫交视频免费看高清 | 捣出白浆h1v1| 日日夜夜操网爽| 久久中文看片网| 男女高潮啪啪啪动态图| 国产精品国产av在线观看| 国产男人的电影天堂91| 国产在线观看jvid| 少妇的丰满在线观看| 久久精品国产a三级三级三级| 日韩有码中文字幕| 午夜视频精品福利| 18禁黄网站禁片午夜丰满| 12—13女人毛片做爰片一| 制服诱惑二区| 久久久久久免费高清国产稀缺| 日韩大码丰满熟妇| 亚洲欧美激情在线| 成年人黄色毛片网站| 韩国高清视频一区二区三区| 制服诱惑二区| 中文字幕制服av| 波多野结衣一区麻豆| 另类精品久久| 欧美日韩成人在线一区二区| 999久久久精品免费观看国产| 黄片大片在线免费观看| 婷婷色av中文字幕| 亚洲国产欧美一区二区综合| 精品少妇内射三级| 青春草视频在线免费观看| 又大又爽又粗| 久久国产精品人妻蜜桃| 亚洲伊人色综图| 国产成人精品久久二区二区91| 9色porny在线观看| 91成人精品电影| 18在线观看网站| 午夜91福利影院| 国产av一区二区精品久久| 国产精品秋霞免费鲁丝片| 欧美日韩精品网址| 亚洲色图综合在线观看| 国产在线一区二区三区精| 啪啪无遮挡十八禁网站| 日日爽夜夜爽网站| 免费少妇av软件| 亚洲美女黄色视频免费看| 在线十欧美十亚洲十日本专区| 欧美97在线视频| 天天躁夜夜躁狠狠躁躁| 中文字幕人妻丝袜制服| 欧美一级毛片孕妇| 日韩一区二区三区影片| 欧美日韩中文字幕国产精品一区二区三区 | 欧美黄色淫秽网站| 精品国产超薄肉色丝袜足j| 国产精品av久久久久免费| av天堂久久9| 新久久久久国产一级毛片| 制服诱惑二区| 中文字幕精品免费在线观看视频| 国产淫语在线视频| 日韩有码中文字幕| 亚洲专区国产一区二区| 亚洲av电影在线进入| 老司机福利观看| 操出白浆在线播放| 女性被躁到高潮视频| 男女床上黄色一级片免费看| 国产精品亚洲av一区麻豆| 日韩一区二区三区影片| 欧美激情极品国产一区二区三区| 搡老乐熟女国产| 成年动漫av网址| 午夜老司机福利片| 国产视频一区二区在线看| 天堂中文最新版在线下载| 一级毛片女人18水好多| h视频一区二区三区| 日本一区二区免费在线视频| 女性被躁到高潮视频| 国产欧美亚洲国产| 爱豆传媒免费全集在线观看| 99热全是精品| 美女脱内裤让男人舔精品视频| 欧美另类一区| 丝袜喷水一区| 国产黄频视频在线观看| 超碰成人久久| 午夜福利免费观看在线| 青青草视频在线视频观看| 国产主播在线观看一区二区| 性少妇av在线| 丝袜在线中文字幕| 亚洲精华国产精华精| 丝袜美足系列| 免费在线观看完整版高清| 两人在一起打扑克的视频| 国产无遮挡羞羞视频在线观看| 男女免费视频国产| 窝窝影院91人妻| 亚洲伊人色综图| 欧美日韩成人在线一区二区| 夜夜骑夜夜射夜夜干| 国产一卡二卡三卡精品| 成年女人毛片免费观看观看9 | 女人久久www免费人成看片| 欧美久久黑人一区二区| 日韩电影二区| 久久久久久亚洲精品国产蜜桃av| 欧美日韩黄片免| av超薄肉色丝袜交足视频| 亚洲自偷自拍图片 自拍| 成人亚洲精品一区在线观看| 国产精品一区二区免费欧美 | 午夜成年电影在线免费观看| 十八禁高潮呻吟视频| 成人黄色视频免费在线看| 亚洲专区国产一区二区| 久久久久国内视频| 国产成人欧美| 性少妇av在线| 久久精品aⅴ一区二区三区四区| 久久国产精品人妻蜜桃| 90打野战视频偷拍视频| 亚洲一区中文字幕在线| 欧美亚洲 丝袜 人妻 在线| 婷婷色av中文字幕| 亚洲av男天堂| 考比视频在线观看| 日日爽夜夜爽网站| 日本a在线网址| 九色亚洲精品在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻,人人澡人人爽秒播| 老司机午夜福利在线观看视频 | 久久精品人人爽人人爽视色| 在线 av 中文字幕| 亚洲美女黄色视频免费看| 日韩制服丝袜自拍偷拍| 精品欧美一区二区三区在线| 日日摸夜夜添夜夜添小说| 在线精品无人区一区二区三| 搡老熟女国产l中国老女人| 欧美精品一区二区大全| 一区二区三区激情视频| 三级毛片av免费| 少妇精品久久久久久久| 欧美大码av| 女人久久www免费人成看片| 中文字幕精品免费在线观看视频| 丝袜人妻中文字幕| 久久久久久久国产电影| 国产成人欧美在线观看 | 日韩视频在线欧美| 国产成人av激情在线播放| 久久国产精品大桥未久av| 9热在线视频观看99| 国产精品偷伦视频观看了| 老熟女久久久| 2018国产大陆天天弄谢| 国产一区二区三区综合在线观看| 精品高清国产在线一区| 欧美精品一区二区大全| 午夜视频精品福利| 老司机福利观看| 人人妻人人爽人人添夜夜欢视频| av国产精品久久久久影院| 99国产综合亚洲精品| www日本在线高清视频| 久久精品亚洲av国产电影网| 涩涩av久久男人的天堂| av天堂久久9| 一个人免费在线观看的高清视频 | 99精国产麻豆久久婷婷| 黑人巨大精品欧美一区二区mp4| 国产极品粉嫩免费观看在线| 久久狼人影院| 国产精品免费大片| 久热这里只有精品99| 色综合欧美亚洲国产小说| 无遮挡黄片免费观看| 一本—道久久a久久精品蜜桃钙片| 亚洲国产成人一精品久久久| 大香蕉久久成人网| 欧美成人午夜精品| 美女主播在线视频| 男人添女人高潮全过程视频| 免费在线观看日本一区| 满18在线观看网站| 妹子高潮喷水视频| 最新在线观看一区二区三区| 黄片小视频在线播放| 在线观看舔阴道视频| 国产精品 国内视频| 亚洲专区中文字幕在线| 法律面前人人平等表现在哪些方面 | 性少妇av在线| 国产91精品成人一区二区三区 | 一个人免费看片子| 青青草视频在线视频观看| 精品国产乱码久久久久久男人| 国产成人欧美在线观看 | 97在线人人人人妻| 久久影院123| 99国产综合亚洲精品| 99热全是精品| 女人精品久久久久毛片| 亚洲男人天堂网一区| 国产成+人综合+亚洲专区| 这个男人来自地球电影免费观看| 97在线人人人人妻| 亚洲国产精品一区三区| 中文字幕另类日韩欧美亚洲嫩草| 另类精品久久| 91麻豆av在线| 日韩三级视频一区二区三区| 欧美日韩av久久| 视频区图区小说| 成年动漫av网址| 国产免费视频播放在线视频| 欧美精品一区二区大全| 自线自在国产av| 国产片内射在线| 99精国产麻豆久久婷婷| 动漫黄色视频在线观看| 777米奇影视久久| 十八禁网站免费在线| www日本在线高清视频| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩精品免费视频一区二区三区| 真人做人爱边吃奶动态| 中国美女看黄片| 丝瓜视频免费看黄片| 亚洲精品国产色婷婷电影| 最近中文字幕2019免费版| 秋霞在线观看毛片| 成人国语在线视频| 深夜精品福利| a级毛片黄视频| 99久久国产精品久久久| kizo精华| av电影中文网址| 欧美日韩视频精品一区| 美女午夜性视频免费| 啦啦啦中文免费视频观看日本| 精品福利永久在线观看| 乱人伦中国视频| 国产高清国产精品国产三级| 国产精品av久久久久免费| 国产无遮挡羞羞视频在线观看| 777米奇影视久久| 五月开心婷婷网| 黑人巨大精品欧美一区二区mp4| 欧美黄色淫秽网站| 日本黄色日本黄色录像| av免费在线观看网站| 亚洲国产av新网站| 亚洲av美国av| 亚洲国产精品成人久久小说| 亚洲国产av新网站| 高清视频免费观看一区二区| 欧美性长视频在线观看| 国产国语露脸激情在线看| 国产免费av片在线观看野外av| 一区二区三区乱码不卡18| 国产黄频视频在线观看| 美女中出高潮动态图| 精品卡一卡二卡四卡免费| 国产成人影院久久av| 国产亚洲精品第一综合不卡| 免费在线观看日本一区| 1024视频免费在线观看| 精品亚洲成a人片在线观看| 久久久久久久精品精品| 国产精品 国内视频| 91成人精品电影| 一区二区三区四区激情视频| 久久久久久久久久久久大奶| 国产在线一区二区三区精| 国产日韩一区二区三区精品不卡| 在线看a的网站| 人人妻,人人澡人人爽秒播| 50天的宝宝边吃奶边哭怎么回事| 国产成人精品无人区| 精品亚洲乱码少妇综合久久| 久久久久久久精品精品| 男女边摸边吃奶| 国产精品久久久av美女十八| 91精品伊人久久大香线蕉| 大香蕉久久成人网| 一区福利在线观看| 精品国产国语对白av| 69精品国产乱码久久久|