• <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柱濃度趨勢與時空分布
    臭氧分子如是說
    在线观看一区二区三区| 精品久久久久久,| 性欧美人与动物交配| 日韩一本色道免费dvd| 国产毛片a区久久久久| 成人国产综合亚洲| 日本黄色视频三级网站网址| 久久久久国产精品人妻aⅴ院| 亚洲无线观看免费| 国产精品99久久久久久久久| 熟女电影av网| 日韩欧美精品v在线| 久久中文看片网| АⅤ资源中文在线天堂| 国产精品久久久久久av不卡| 夜夜夜夜夜久久久久| 窝窝影院91人妻| 免费av不卡在线播放| 欧美日韩精品成人综合77777| 中文字幕熟女人妻在线| 乱系列少妇在线播放| 精品一区二区三区av网在线观看| 1000部很黄的大片| 琪琪午夜伦伦电影理论片6080| 成年人黄色毛片网站| 精品久久久久久久久久久久久| 欧美日韩瑟瑟在线播放| 国产伦精品一区二区三区视频9| 亚洲真实伦在线观看| 欧美日韩精品成人综合77777| 国产免费av片在线观看野外av| 亚洲精品粉嫩美女一区| 日韩欧美精品免费久久| 成人鲁丝片一二三区免费| 真实男女啪啪啪动态图| 免费无遮挡裸体视频| 国产视频一区二区在线看| 免费看光身美女| 久久久久国内视频| 国产精品综合久久久久久久免费| 亚洲欧美日韩东京热| 国产真实乱freesex| 1024手机看黄色片| 久久午夜亚洲精品久久| 亚洲最大成人手机在线| 亚洲va在线va天堂va国产| 毛片一级片免费看久久久久 | 国产精品久久久久久久电影| 国产高清不卡午夜福利| 欧美日韩瑟瑟在线播放| 欧美精品啪啪一区二区三区| 国产探花极品一区二区| 搞女人的毛片| 身体一侧抽搐| 久久久久久久久大av| 狂野欧美激情性xxxx在线观看| 一边摸一边抽搐一进一小说| 久久久国产成人精品二区| 国产高清三级在线| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩无卡精品| 嫩草影院精品99| 毛片女人毛片| 精品人妻一区二区三区麻豆 | 久久久久久久午夜电影| 男人和女人高潮做爰伦理| 极品教师在线视频| 亚洲精华国产精华液的使用体验 | 免费看日本二区| 免费人成在线观看视频色| 黄色日韩在线| 最新中文字幕久久久久| 午夜亚洲福利在线播放| x7x7x7水蜜桃| 亚洲欧美日韩高清在线视频| 日韩中文字幕欧美一区二区| 色综合亚洲欧美另类图片| 国国产精品蜜臀av免费| 午夜精品久久久久久毛片777| 免费在线观看日本一区| 成人特级黄色片久久久久久久| 色在线成人网| 床上黄色一级片| 一区二区三区四区激情视频 | 成人三级黄色视频| 国产精品爽爽va在线观看网站| 国产欧美日韩一区二区精品| 亚洲三级黄色毛片| 两个人的视频大全免费| 国产精品三级大全| 国产真实乱freesex| 国产 一区 欧美 日韩| 黄色配什么色好看| av国产免费在线观看| 亚洲最大成人av| 午夜精品在线福利| 中文字幕人妻熟人妻熟丝袜美| 国产探花极品一区二区| 国产不卡一卡二| 18禁在线播放成人免费| 99热6这里只有精品| 欧美最黄视频在线播放免费| 女同久久另类99精品国产91| 国产女主播在线喷水免费视频网站 | 国产亚洲精品av在线| 亚洲天堂国产精品一区在线| 久99久视频精品免费| 日本黄色片子视频| 久久久国产成人免费| 午夜精品在线福利| 久久久色成人| 可以在线观看的亚洲视频| 久久久国产成人免费| 午夜爱爱视频在线播放| 桃红色精品国产亚洲av| 午夜福利高清视频| 午夜影院日韩av| av天堂在线播放| 又粗又爽又猛毛片免费看| 99国产极品粉嫩在线观看| 波野结衣二区三区在线| 久久人妻av系列| 免费人成在线观看视频色| 小蜜桃在线观看免费完整版高清| 韩国av一区二区三区四区| 亚洲最大成人中文| 国产精品久久久久久久电影| 欧美日韩黄片免| 国产三级中文精品| 村上凉子中文字幕在线| 五月伊人婷婷丁香| 99热网站在线观看| 亚洲熟妇熟女久久| 99九九线精品视频在线观看视频| 亚洲性夜色夜夜综合| 最后的刺客免费高清国语| 国内精品宾馆在线| 中文亚洲av片在线观看爽| h日本视频在线播放| 国产午夜精品论理片| 亚洲成av人片在线播放无| 黄色一级大片看看| 波野结衣二区三区在线| 大型黄色视频在线免费观看| 色吧在线观看| 亚洲五月天丁香| 一边摸一边抽搐一进一小说| 在线看三级毛片| 天堂av国产一区二区熟女人妻| 不卡一级毛片| 在线免费观看的www视频| 国模一区二区三区四区视频| 久久草成人影院| 免费无遮挡裸体视频| 亚州av有码| bbb黄色大片| 一个人看视频在线观看www免费| 久久久久久伊人网av| 久久久国产成人免费| 国产精品亚洲美女久久久| 国产主播在线观看一区二区| 51国产日韩欧美| 亚洲va日本ⅴa欧美va伊人久久| 九九久久精品国产亚洲av麻豆| 亚洲av免费高清在线观看| 亚洲成av人片在线播放无| 久久这里只有精品中国| 美女 人体艺术 gogo| 不卡视频在线观看欧美| 天堂√8在线中文| 欧美色欧美亚洲另类二区| 国产乱人伦免费视频| 亚洲熟妇熟女久久| 国产精品永久免费网站| 人人妻,人人澡人人爽秒播| 丝袜美腿在线中文| 看片在线看免费视频| 久久国产乱子免费精品| 18禁黄网站禁片午夜丰满| 欧美日韩瑟瑟在线播放| 91久久精品国产一区二区三区| 成年免费大片在线观看| 永久网站在线| 免费人成视频x8x8入口观看| 小蜜桃在线观看免费完整版高清| 国产精品,欧美在线| 三级毛片av免费| 欧美日韩瑟瑟在线播放| 淫秽高清视频在线观看| 淫妇啪啪啪对白视频| 久久久午夜欧美精品| 最近最新免费中文字幕在线| 欧美激情国产日韩精品一区| 床上黄色一级片| 欧美另类亚洲清纯唯美| 国内少妇人妻偷人精品xxx网站| 两个人视频免费观看高清| av在线老鸭窝| 国产精品无大码| 日韩高清综合在线| 黄色配什么色好看| 国产精品精品国产色婷婷| 国产精品永久免费网站| 神马国产精品三级电影在线观看| 免费观看的影片在线观看| 在线播放无遮挡| 国产免费av片在线观看野外av| 午夜激情欧美在线| 日本-黄色视频高清免费观看| a级毛片免费高清观看在线播放| 蜜桃久久精品国产亚洲av| 国产女主播在线喷水免费视频网站 | 日本 av在线| 国产成人aa在线观看| 少妇高潮的动态图| 欧美性猛交黑人性爽| 亚洲国产精品久久男人天堂| 色在线成人网| 欧美xxxx性猛交bbbb| 精品福利观看| 日本免费a在线| 男女那种视频在线观看| 精品国产三级普通话版| 麻豆av噜噜一区二区三区| 熟女电影av网| 欧美日韩中文字幕国产精品一区二区三区| 麻豆国产av国片精品| 日韩,欧美,国产一区二区三区 | 国产精品伦人一区二区| 日韩人妻高清精品专区| 久久国内精品自在自线图片| 亚洲精华国产精华液的使用体验 | 国产亚洲精品av在线| 成人国产一区最新在线观看| 男人的好看免费观看在线视频| 中文资源天堂在线| 久久精品夜夜夜夜夜久久蜜豆| 色精品久久人妻99蜜桃| 亚洲第一电影网av| 老师上课跳d突然被开到最大视频| 男人和女人高潮做爰伦理| 国产视频内射| 精品国内亚洲2022精品成人| 老熟妇仑乱视频hdxx| 免费电影在线观看免费观看| 高清毛片免费观看视频网站| 国产精品久久久久久久久免| 欧美日韩国产亚洲二区| 色在线成人网| 搡老妇女老女人老熟妇| 欧美3d第一页| 亚洲va日本ⅴa欧美va伊人久久| 国产伦精品一区二区三区视频9| 男人的好看免费观看在线视频| 国产精品电影一区二区三区| 999久久久精品免费观看国产| www.www免费av| 99精品在免费线老司机午夜| 欧美日韩黄片免| 国内毛片毛片毛片毛片毛片| 亚洲欧美精品综合久久99| 国内精品一区二区在线观看| 91av网一区二区| 亚洲欧美日韩卡通动漫| 免费黄网站久久成人精品| 波野结衣二区三区在线| 亚洲精品色激情综合| 久久久国产成人免费| 国产精品不卡视频一区二区| 国产 一区精品| 日韩av在线大香蕉| 亚洲第一电影网av| 12—13女人毛片做爰片一| 亚洲精品国产成人久久av| 色视频www国产| 人妻制服诱惑在线中文字幕| 国产精品一及| 最近在线观看免费完整版| av在线亚洲专区| 99久久精品热视频| 欧美色视频一区免费| 日本熟妇午夜| 美女黄网站色视频| 色播亚洲综合网| 波多野结衣高清作品| 真人一进一出gif抽搐免费| 国产探花在线观看一区二区| 老熟妇仑乱视频hdxx| 舔av片在线| 热99re8久久精品国产| 亚洲在线自拍视频| 男人舔奶头视频| 一夜夜www| 两人在一起打扑克的视频| 成人永久免费在线观看视频| 搡老岳熟女国产| 一边摸一边抽搐一进一小说| 天堂影院成人在线观看| 一区二区三区高清视频在线| 麻豆国产97在线/欧美| 狂野欧美白嫩少妇大欣赏| 国语自产精品视频在线第100页| 欧美黑人欧美精品刺激| 一级黄片播放器| 两性午夜刺激爽爽歪歪视频在线观看| 此物有八面人人有两片| 欧美日韩乱码在线| 亚洲av成人精品一区久久| 日本与韩国留学比较| 日本-黄色视频高清免费观看| 欧美性感艳星| 乱码一卡2卡4卡精品| 色在线成人网| 在线观看午夜福利视频| 日日夜夜操网爽| 无遮挡黄片免费观看| 亚洲最大成人手机在线| 国产亚洲精品av在线| 久久中文看片网| 性欧美人与动物交配| 中文亚洲av片在线观看爽| 免费大片18禁| 国产又黄又爽又无遮挡在线| 国产精品久久久久久亚洲av鲁大| 噜噜噜噜噜久久久久久91| 有码 亚洲区| 男人狂女人下面高潮的视频| 亚洲专区国产一区二区| 少妇的逼好多水| 成年女人毛片免费观看观看9| 国产高清不卡午夜福利| 天堂影院成人在线观看| 韩国av在线不卡| 日本黄色视频三级网站网址| 亚洲av免费在线观看| 国产成人影院久久av| 亚洲av中文字字幕乱码综合| 2021天堂中文幕一二区在线观| 精品无人区乱码1区二区| 老司机深夜福利视频在线观看| 久久精品久久久久久噜噜老黄 | 日韩欧美免费精品| 最近最新中文字幕大全电影3| 国产午夜福利久久久久久| 欧美黑人欧美精品刺激| 国产精品一区www在线观看 | 国产欧美日韩精品亚洲av| 又爽又黄无遮挡网站| 老熟妇乱子伦视频在线观看| 久久久久久久久久黄片| 国产淫片久久久久久久久| 在线播放无遮挡| 51国产日韩欧美| 最近最新中文字幕大全电影3| 1000部很黄的大片| 日日撸夜夜添| 久久国产乱子免费精品| 久久国内精品自在自线图片| 赤兔流量卡办理| 精品久久久久久久久亚洲 | 亚洲专区国产一区二区| 久久欧美精品欧美久久欧美| av福利片在线观看| 亚洲av.av天堂| 成人av在线播放网站| 少妇高潮的动态图| 琪琪午夜伦伦电影理论片6080| 久久精品国产99精品国产亚洲性色| av在线蜜桃| 欧美色视频一区免费| 亚洲avbb在线观看| 日本成人三级电影网站| 亚洲熟妇中文字幕五十中出| 国产麻豆成人av免费视频| 亚洲最大成人中文| 琪琪午夜伦伦电影理论片6080| 久久久国产成人精品二区| 婷婷精品国产亚洲av| 麻豆成人午夜福利视频| 日韩 亚洲 欧美在线| 少妇人妻精品综合一区二区 | 婷婷精品国产亚洲av在线| 国产欧美日韩精品亚洲av| 国内精品久久久久久久电影| 亚洲av五月六月丁香网| 1024手机看黄色片| 午夜精品一区二区三区免费看| 欧美日韩综合久久久久久 | 丰满的人妻完整版| 亚洲成人久久爱视频| 精品福利观看| 精品一区二区免费观看| 永久网站在线| 色综合站精品国产| 午夜久久久久精精品| 国产高清视频在线播放一区| 免费搜索国产男女视频| av专区在线播放| 免费av毛片视频| 国产人妻一区二区三区在| av在线亚洲专区| 一区福利在线观看| 成年免费大片在线观看| 国产精品98久久久久久宅男小说| 久久久久久久午夜电影| 亚洲成人久久爱视频| 国内精品久久久久精免费| 三级国产精品欧美在线观看| 久久人人爽人人爽人人片va| 女人被狂操c到高潮| 久久久久国产精品人妻aⅴ院| 少妇被粗大猛烈的视频| 久9热在线精品视频| 日本熟妇午夜| 丰满乱子伦码专区| 精品久久久久久久久久久久久| 欧美+亚洲+日韩+国产| 日日干狠狠操夜夜爽| 亚洲三级黄色毛片| 久久精品久久久久久噜噜老黄 | 性插视频无遮挡在线免费观看| 成人一区二区视频在线观看| 一区福利在线观看| 麻豆av噜噜一区二区三区| 亚洲中文字幕日韩| 99国产精品一区二区蜜桃av| 中文字幕av在线有码专区| 国产高清视频在线观看网站| 久久精品国产亚洲av天美| 亚洲成人中文字幕在线播放| 国产美女午夜福利| 永久网站在线| 性色avwww在线观看| 99久久九九国产精品国产免费| 亚洲性夜色夜夜综合| 日日夜夜操网爽| 久久婷婷人人爽人人干人人爱| 亚洲va在线va天堂va国产| 特级一级黄色大片| 久久6这里有精品| 亚洲av二区三区四区| 一区二区三区免费毛片| 国产一区二区三区视频了| 无人区码免费观看不卡| eeuss影院久久| 在线播放国产精品三级| 亚洲18禁久久av| 波多野结衣巨乳人妻| 免费看美女性在线毛片视频| 黄色日韩在线| 禁无遮挡网站| 欧美又色又爽又黄视频| 精品午夜福利视频在线观看一区| 美女黄网站色视频| 欧美高清成人免费视频www| 狠狠狠狠99中文字幕| 免费av不卡在线播放| 高清毛片免费观看视频网站| 亚洲精品亚洲一区二区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲综合色惰| 日日摸夜夜添夜夜添小说| 国产精品自产拍在线观看55亚洲| 又黄又爽又免费观看的视频| 亚洲av免费高清在线观看| 神马国产精品三级电影在线观看| 性欧美人与动物交配| 国产伦在线观看视频一区| 极品教师在线免费播放| 少妇高潮的动态图| 欧美成人a在线观看| 欧美三级亚洲精品| 国产日本99.免费观看| 一区福利在线观看| 成人一区二区视频在线观看| 我要搜黄色片| 国产熟女欧美一区二区| 亚洲人成伊人成综合网2020| 久久精品国产自在天天线| 中文在线观看免费www的网站| 床上黄色一级片| 欧美+日韩+精品| 91在线精品国自产拍蜜月| 天天一区二区日本电影三级| 欧美日韩乱码在线| 啦啦啦啦在线视频资源| 欧美潮喷喷水| 五月玫瑰六月丁香| 乱码一卡2卡4卡精品| 男女下面进入的视频免费午夜| 两人在一起打扑克的视频| 日韩人妻高清精品专区| 亚洲成人久久爱视频| 12—13女人毛片做爰片一| 国产精品美女特级片免费视频播放器| 在线观看66精品国产| 久久午夜福利片| 看十八女毛片水多多多| 国产伦人伦偷精品视频| 欧美人与善性xxx| 国内毛片毛片毛片毛片毛片| av在线亚洲专区| 乱人视频在线观看| 免费av观看视频| 欧美性感艳星| 久久久久久久久久成人| 日本爱情动作片www.在线观看 | 国产精品人妻久久久久久| 久久久久久大精品| 真人做人爱边吃奶动态| 色哟哟·www| 日日夜夜操网爽| 有码 亚洲区| 午夜视频国产福利| 91久久精品国产一区二区三区| 窝窝影院91人妻| 亚洲性夜色夜夜综合| 日本熟妇午夜| 男女之事视频高清在线观看| 国产成人一区二区在线| 观看美女的网站| 日日摸夜夜添夜夜添小说| 露出奶头的视频| 久久精品国产亚洲av涩爱 | 黄色日韩在线| 久久天躁狠狠躁夜夜2o2o| 国产精品98久久久久久宅男小说| 亚洲五月天丁香| 变态另类成人亚洲欧美熟女| 18禁裸乳无遮挡免费网站照片| 中文亚洲av片在线观看爽| 日本欧美国产在线视频| 成人av一区二区三区在线看| 午夜精品久久久久久毛片777| 在线观看av片永久免费下载| 亚洲人成伊人成综合网2020| 午夜福利18| 色av中文字幕| 一区福利在线观看| 国产久久久一区二区三区| 日韩精品有码人妻一区| 高清在线国产一区| 亚洲av免费在线观看| 窝窝影院91人妻| 色综合站精品国产| 国产伦在线观看视频一区| 最近视频中文字幕2019在线8| 中文资源天堂在线| 亚洲性夜色夜夜综合| 可以在线观看毛片的网站| 无遮挡黄片免费观看| 少妇的逼好多水| 亚洲av成人精品一区久久| 色吧在线观看| 精品国内亚洲2022精品成人| 嫩草影院精品99| 伊人久久精品亚洲午夜| 国产精品自产拍在线观看55亚洲| 中文资源天堂在线| 我要搜黄色片| 老熟妇乱子伦视频在线观看| 精品免费久久久久久久清纯| 黄色配什么色好看| 中国美女看黄片| 日本黄大片高清| 成人av在线播放网站| 99国产极品粉嫩在线观看| 国产精品久久久久久久久免| 婷婷六月久久综合丁香| 精品人妻一区二区三区麻豆 | 精品不卡国产一区二区三区| 好男人在线观看高清免费视频| 日本黄大片高清| 午夜激情福利司机影院| 亚洲美女黄片视频| 欧美3d第一页| 少妇被粗大猛烈的视频| 精品久久久久久久久久久久久| 国产精品一区二区三区四区久久| 啦啦啦韩国在线观看视频| 一级毛片久久久久久久久女| 亚洲av五月六月丁香网| 色精品久久人妻99蜜桃| av中文乱码字幕在线| 九色国产91popny在线| 国产高清有码在线观看视频| ponron亚洲| 一本精品99久久精品77| 我要搜黄色片| 国产极品精品免费视频能看的| 别揉我奶头~嗯~啊~动态视频| av.在线天堂| 国产精品福利在线免费观看| 99热这里只有是精品在线观看| 美女大奶头视频| 中出人妻视频一区二区| 我的女老师完整版在线观看| 少妇的逼水好多| 久久久久久久精品吃奶| 身体一侧抽搐| 麻豆av噜噜一区二区三区| 国产伦精品一区二区三区视频9| 欧美一级a爱片免费观看看| 色综合站精品国产| 亚洲国产精品合色在线| 亚洲va在线va天堂va国产| 免费在线观看影片大全网站| xxxwww97欧美| 亚洲精品成人久久久久久| 如何舔出高潮| 美女免费视频网站| 久久久久久久久久久丰满 | 两个人视频免费观看高清| 在线天堂最新版资源| 少妇人妻精品综合一区二区 |