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

    50年來中國臭氧總量時空變化特征

    2022-07-18 13:56:12謝靜晗肖鐘湧
    中國環(huán)境科學(xué) 2022年7期
    關(guān)鍵詞:臭氧總量平均值

    謝靜晗,李 颯,肖鐘湧

    50年來中國臭氧總量時空變化特征

    謝靜晗,李 颯,肖鐘湧*

    (集美大學(xué)港口與海岸工程學(xué)院,福建 廈門 361021)

    利用來自世界臭氧與紫外輻射數(shù)據(jù)中心的中國區(qū)域6個地基觀測站點(diǎn)數(shù)據(jù),對多傳感器再分析遙感數(shù)據(jù)進(jìn)行驗(yàn)證,并基于驗(yàn)證后的遙感數(shù)據(jù)分析了1971~2020年中國區(qū)域臭氧總量不同尺度的時空變化特征.結(jié)果表明,50a來中國區(qū)域臭氧總量呈現(xiàn)輕微的下降趨勢.年平均臭氧總量在1978年和1993年分別出現(xiàn)最大值(347.5±53.8)DU和最小值(291.9±29.5)DU,在1971~1978年、1978~1993年、1993~2020年,這3個時段年平均臭氧總量在整個中國區(qū)域分別是增長、減少、增長.月平均臭氧總量隨季節(jié)變化呈現(xiàn)出正弦曲線形態(tài),在3月和10月分別出現(xiàn)峰值(約338DU)和谷值(約285DU).中國區(qū)域臭氧總量在空間上呈現(xiàn)由東北向西南遞減的緯向條帶狀分布.在40°N以北的東北部地區(qū),該值可達(dá)360DU以上.中國區(qū)域50a月平均臭氧總量同樣呈現(xiàn)緯向條帶狀分布.此外,時間變異系數(shù)和空間變異系數(shù)隨季節(jié)的變化規(guī)律相似,夏季最小,接著依次是秋季和春季,冬季最大.即臭氧總量的變化和空間差異在夏季都最小.50a期間,不同時段、不同區(qū)域臭氧總量的變化趨勢各不相同.在1971~1978年,臭氧總量的增長量和增長率都呈現(xiàn)由北向南遞減的緯向條帶狀分布.在40°N以北的相對高值地區(qū)最大增加了56DU,約為16%;而在30°N以南的相對低值地區(qū),最小增加了12DU,約為5%.在1978~1993年,減少量和減少率也呈現(xiàn)由北向南遞減的緯度地帶性.在40°N以北的相對高值地區(qū)最大減少了93DU,約為22%;而在30°N以南的相對低值地區(qū),最小減少了11DU,約為4%.在1993~2020年,西北地區(qū)出現(xiàn)最大增長,增長量為18DU,約為6%;東南地區(qū)出現(xiàn)最小增長,增長量為4DU,約為1%.

    臭氧總量;時序變化;季節(jié)變化;空間特征;中國區(qū)域

    臭氧是大氣中的痕量氣體,主要分布在平流層和對流層.其中平流層臭氧約占大氣臭氧總量的90%,余下的分布在對流層[1-3].雖然臭氧在大氣中的體積占比不到0.001%,但卻對地球影響深遠(yuǎn)[4].平流層臭氧能夠大量吸收太陽紫外輻射,從而充當(dāng)?shù)厍蛏谋Wo(hù)傘,被稱為“好臭氧”[5-8].此外,平流層臭氧對紫外輻射的大量吸收使得平流層大氣增溫,進(jìn)而影響大氣環(huán)流[9-11].對流層臭氧被稱為“壞臭氧”,因其強(qiáng)氧化性,超過一定濃度會危害人類健康以及降低農(nóng)作物產(chǎn)量,是主要的大氣污染物之一[12-16].

    目前,大氣臭氧的研究主要基于地面和衛(wèi)星兩種平臺[17-19].由于地面平臺的觀測儀器易于維護(hù)和標(biāo)定,能夠?qū)Τ粞踹M(jìn)行長時間序列的定點(diǎn)觀測,其臭氧數(shù)據(jù)的穩(wěn)定性和時序性較好[17,20].但地面平臺的建設(shè)受到地理環(huán)境的制約,在高山、海洋等地區(qū)無法建立觀測站,因此其臭氧數(shù)據(jù)全球覆蓋率較差.衛(wèi)星平臺很好的彌補(bǔ)了這一缺陷[21],其臭氧數(shù)據(jù)具備全球覆蓋率高、數(shù)據(jù)充足等特點(diǎn)[21-23].自1970年BUV衛(wèi)星傳感器(搭載于Nimbus-4衛(wèi)星)發(fā)射以來, TOMS和SBUV(Nimbus-7)、GOME(ERS-2)、SCIAMACHY(Envisat)等衛(wèi)星傳感器隨后相繼發(fā)射,形成一系列觀測平臺[24].鑒于地面和衛(wèi)星這兩種平臺的觀測各具特點(diǎn)并存在互補(bǔ)性,有些研究基于地面觀測數(shù)據(jù)檢驗(yàn)衛(wèi)星反演數(shù)據(jù)的可靠性[25-26],有些研究更進(jìn)一步,運(yùn)用驗(yàn)證或校正后的衛(wèi)星數(shù)據(jù)分析臭氧總量的時空變化[27-28].

    臭氧總量的時空變化深受人類活動影響.20世紀(jì)以來,氯氟烴類(CFCs)被大量應(yīng)用于制冷系統(tǒng)[29],且在70年代被發(fā)現(xiàn)能長期存留于大氣并破壞平流層臭氧[30-31].為保護(hù)臭氧層,聯(lián)合國環(huán)境署制定《維也納公約》且呼吁會員國簽署《蒙特利爾議定書》并對其進(jìn)行多次修訂[32-33].在經(jīng)歷了這一系列與臭氧息息相關(guān)的事件至今,研究長時間序列的臭氧總量變化顯得尤為重要,而中國這一區(qū)域受到廣泛關(guān)注.已有研究調(diào)查中國區(qū)域一定時間序列的臭氧總量時空變化[27,34-35],也有研究針對中國局部地區(qū)[1,9,11,28,36-37],但罕有研究探討中國區(qū)域迄今為止長時間序列的臭氧總量變化特征.此外,已有研究較少使用MSR-2這一高精準(zhǔn)度的多傳感器再分析遙感數(shù)據(jù)[9].鑒于此,本研究主要基于MSR-2分析中國區(qū)域1971~2020年這50a來臭氧總量在不同時空尺度的變化特征,為充分了解中國區(qū)域臭氧總量情況提供參考,也為臭氧相關(guān)政策的制定提供支撐.

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

    1.1 數(shù)據(jù)

    本研究的數(shù)據(jù)源為多傳感器再分析(MSR-2)資料(https://www.temis.nl/protocols/o3field/o3field_ msr2.php)以及世界臭氧與紫外輻射數(shù)據(jù)中心(WOUDC)的全球共享數(shù)據(jù).MSR-2是BUV、TOMS、SBUV、GOME、SCIAMACHY、OMI、GOME-2等衛(wèi)星傳感器的臭氧總量觀測數(shù)據(jù)融合產(chǎn)生的再分析數(shù)據(jù)[24].已有研究結(jié)果表明BUV、TOMS、SBUV、GOME、SCIAMACHY、OMI、GOME-2等衛(wèi)星反演的臭氧總量產(chǎn)品,可以構(gòu)建一個長時間序列的數(shù)據(jù)集進(jìn)行區(qū)域臭氧總量的研究[4, 27-28].文件格式為HDF,空間分辨率為1.5°×1°,時間分辨率為天.地面觀測數(shù)據(jù)來自WOUDC,本研究使用的中國區(qū)域6個地面站點(diǎn)的空間分布如圖1所示.結(jié)合表1可知,這6個站點(diǎn)分別是河北香河(XIANGHE)、云南昆明(KUNMING)、青海瓦里關(guān)(MT. WALIGUAN)、浙江臨安(LINAN)、黑龍江龍鳳山(LONGFENG SHAN)以及西藏拉薩(Lhasa).其中,前2個地面站點(diǎn)使用的是Dobson分光光度計,后4個地面站點(diǎn)使用的是Brewer分光光度計.采用臭氧總量數(shù)據(jù)的時間跨度為50a(1971~2020年),單位是DU(Dobson Unit).

    圖1 地基觀測站點(diǎn)分布

    底圖源自國家測繪地理信息局網(wǎng)站(http://zwfw.nasg.gov.cn/approval)審圖號為GS(2016)1554號的標(biāo)準(zhǔn)地圖制作,底圖無修改,下同

    表1 6個地基觀測站點(diǎn)信息

    1.2 方法

    本研究基于中國區(qū)域6個地面站點(diǎn)的觀測數(shù)據(jù)對衛(wèi)星反演數(shù)據(jù)的可靠性進(jìn)行檢驗(yàn).驗(yàn)證數(shù)據(jù)取樣的方法是將地面站點(diǎn)觀測的日平均臭氧總量與一個像元值相對應(yīng),然后進(jìn)行線性回歸分析.這些像元值是由衛(wèi)星過境的瞬時觀測值反演得到的.首先,從HDF文件中讀取衛(wèi)星臭氧總量,將其轉(zhuǎn)化為ArcGIS的矢量數(shù)據(jù)存儲文件Shapefile,然后把日尺度的數(shù)據(jù)(共有18263d)合成月尺度和年尺度的數(shù)據(jù),最后利用最近鄰法插值得到分辨率為1°×1°的柵格數(shù)據(jù).基于這些柵格數(shù)據(jù),分析了1971~2020年這50a來中國區(qū)域臭氧總量的時空變化特征.在時間特征的分析中,計算了整個中國區(qū)域的臭氧總量在日尺度、月尺度和年尺度的平均值.在空間特征的分析中,計算了中國區(qū)域的各個像元在月尺度與年尺度的平均值以及年際間的增長率.

    此外,臭氧總量在時間上的波動情況和在空間上的異質(zhì)性分別用時間變異系數(shù)和空間變異系數(shù)來表征.變異系數(shù)(CV)為標(biāo)準(zhǔn)差()與平均數(shù)()的比值,計算公式:

    CV=/

    式中:CV為時間變異系數(shù)時和分別根據(jù)時序變化的數(shù)據(jù)來計算;而CV為空間變異系數(shù)時和分別根據(jù)區(qū)域空間分布的數(shù)據(jù)來計算.

    2 結(jié)果與分析

    2.1 驗(yàn)證

    本研究利用從WOUDC獲取的中國區(qū)域6個地面站點(diǎn)觀測的臭氧總量,對MSR-2資料衛(wèi)星反演的臭氧總量的可靠性進(jìn)行檢驗(yàn).圖2展示了二者的線性回歸分析,回歸參數(shù)如表2所示.6個站點(diǎn)線性回歸的相關(guān)系數(shù)均較高,為0.87~0.97;斜率和截距分別為0.73~0.99,2.57~75.31.其中相關(guān)性最強(qiáng)的是河北香河(站點(diǎn)208),相關(guān)系數(shù)為0.97,斜率和截距分別為0.99和2.57,標(biāo)準(zhǔn)誤差為9.21;最低的是浙江臨安(站點(diǎn)325),相關(guān)系數(shù)為0.87,斜率和截距分別為0.73和75.31,標(biāo)準(zhǔn)誤差為12.00DU.除了衛(wèi)星反演的算法存在誤差外,進(jìn)行比較驗(yàn)證的地面觀測值和衛(wèi)星反演值在時間和空間上也未完全匹配,大氣中臭氧總量時間變化和空間分布不均勻,因此地面觀測和衛(wèi)星反演存在一定的誤差.另外,從圖中可以得出,在較低值時,地面觀測臭氧總量普遍低于衛(wèi)星反演臭氧總量;在較高值時,地面觀測臭氧總量普遍高于衛(wèi)星反演臭氧總量.這可能是由于衛(wèi)星反演算法在計算臭氧總量時,先計算斜氣柱臭氧總量,再通過斜氣柱臭氧總量計算垂直氣柱臭氧總量.在這個過程中,較高臭氧總量值被抑制,而較低臭氧總量值被放大.而地面觀測的儀器直接計算出垂直氣柱臭氧總量.驗(yàn)證分析結(jié)果表明MSR-2資料適合應(yīng)用于中國區(qū)域臭氧總量時空變化特征的研究.

    表2 地面觀測與相應(yīng)衛(wèi)星觀測臭氧總量的線性回歸參數(shù)

    2.2 時序變化

    圖3展示了中國區(qū)域日平均和年平均臭氧總量的變化特征.圖3(b)中的空間變異系數(shù)為中國區(qū)域年平均臭氧總量像元標(biāo)準(zhǔn)差與平均值的比值,描述空間異質(zhì)性;時間變異系數(shù)為中國區(qū)域1年12個月平均臭氧總量標(biāo)準(zhǔn)差與平均值的比值,描述不同年份臭氧總量的季節(jié)變化情況.由圖3可知,50a來中國區(qū)域臭氧總量呈現(xiàn)輕微下降的趨勢.從1971年的(310.4±41.1)DU,空間變異系數(shù)為 0.123,時間變異系數(shù)為0.066,到2020年的(303.8±30.8)DU,空間變異系數(shù)為0.122,時間變異系數(shù)為0.075;年平均臭氧總量減少了6.6DU,空間異質(zhì)性無明顯變化,但年內(nèi)波動出現(xiàn)增大.在1978年出現(xiàn)最大值(347.5±53.8)DU,空間變異系數(shù)為0.155,時間變異系數(shù)為0.098,空間異質(zhì)性和年內(nèi)波動達(dá)到最大.在1971~2020年這50a間,中國區(qū)域年平均臭氧總量最大值出現(xiàn)在1978年的可能原因是70年代初期氯氟烴類被發(fā)現(xiàn)會在大氣中長期存留并破壞平流層臭氧[30-31],國際上開始廣泛關(guān)注臭氧問題并逐步采取行動保護(hù)臭氧層.此外,我國環(huán)境保護(hù)意識從1972年開始萌芽,并在 1973 年召開第一次全國環(huán)境保護(hù)會議,從此開始了環(huán)境保護(hù)工作[38],這也為臭氧的積累創(chuàng)造有利條件.而臭氧的積累會使得臭氧的空間異質(zhì)性和年內(nèi)波動變大.1978年12月十一屆三中全會后,中國開始實(shí)行改革開放,工業(yè)化和城市化進(jìn)程加快[16].在發(fā)展過程中,人為活動排放氯氟烴類等臭氧消耗物質(zhì)[39],導(dǎo)致臭氧總量開始減少.

    中國區(qū)域年平均臭氧總量在1993年出現(xiàn)最小值(291.9±29.5)DU,空間變異系數(shù)為0.096,時間變異系數(shù)為0.041,臭氧總量空間異質(zhì)性和年內(nèi)波動達(dá)到最小.最小值出現(xiàn)的年份與韋惠紅等[34]的研究結(jié)果一致.韋惠紅等將中國劃分為7個區(qū)域并分析其上空臭氧總量的變化,發(fā)現(xiàn)這些區(qū)域都在1993年達(dá)到最小值.這可能是由于1991年6月菲律賓Pinatubo火山爆發(fā)產(chǎn)生大量氣溶膠[40-41],高濃度的氣溶膠使得自由基ClO的濃度升高,并長期存留在大氣中,劇烈消耗了臭氧[42].世界氣象組織在評估臭氧消耗時指出全球臭氧總量在1992~1993年間達(dá)到了最低值[43].臭氧大尺度的大量消耗可能是導(dǎo)致1993年中國區(qū)域臭氧總量空間異質(zhì)性和年內(nèi)波動小的主要原因.從1993年往后,臭氧總量有所回升,這主要是因?yàn)镻inatubo火山爆發(fā)所產(chǎn)生的氣溶膠在隨后幾年被逐漸清除[33].此外,《蒙特利爾議定書》及其修正案中有關(guān)限制臭氧消耗物質(zhì)生產(chǎn)和使用的規(guī)定也起了成效.研究表明,等效平流層氯化物從90年代中后期以來,呈現(xiàn)下降趨勢[33].最后,工業(yè)生產(chǎn)等活動排放的臭氧進(jìn)入大氣也是致使臭氧總量增加的原因之一.

    圖3 中國區(qū)域日均和年均臭氧總量變化特征

    2.3 季節(jié)變化

    圖4和圖5分別為中國區(qū)域多年月平均臭氧總量變化及其空間分布.圖4中的空間變異系數(shù)為中國區(qū)域50a月平均臭氧總量像元標(biāo)準(zhǔn)差與平均值的比值,描述空間異質(zhì)性;時間變異系數(shù)為中國區(qū)域50a月平均臭氧總量標(biāo)準(zhǔn)差與平均值的比值,描述不同月份臭氧總量的變化情況.由圖可知,中國區(qū)域多年月平均臭氧總量存在明顯的時空變化特征.月平均臭氧總量隨季節(jié)變化呈現(xiàn)出正弦曲線形態(tài).峰值和谷值分別出現(xiàn)在3月和10月,約為338和285DU,相差53DU.這一變化模式與張瑩等[27]對中國區(qū)域30年(1981~2010年)月平均臭氧總量變化研究的結(jié)果相符.臭氧總量的季節(jié)變化主要受光化學(xué)反應(yīng)和動力輸送的影響[44-47].太陽輻射是光化學(xué)反應(yīng)的成分之一,直接影響臭氧的生成和消耗[44].臭氧主要生成于熱帶地區(qū),通過大氣環(huán)流,輸送到中高緯度地區(qū)[44-45].臭氧總量的變化模式反映出從3~10月,臭氧處于消耗狀態(tài),這是因?yàn)樘栞椛漭^強(qiáng)時,光化學(xué)損失起主導(dǎo)作用;其余時間臭氧處于積累狀態(tài),這主要是受到動力輸送的影響[47].

    圖4 中國區(qū)域多年月平均臭氧總量變化

    圖5 多年月平均臭氧總量空間分布

    圖6 時間變異系數(shù)的空間分布

    由時間變異系數(shù)可知,臭氧總量在夏季(6~8月)的變化最小,接著依次是秋季(9~11月)和春季(3~5月),冬季(12~2月)最大.這與Chen等[44]對1988~1997年瑞典西南部的臭氧時空變化研究的結(jié)果一致.空間變異系數(shù)呈現(xiàn)“V”字型分布,夏季,特別是8月,臭氧總量的空間差異最小,接著依次是秋季和春季,冬季空間差異最大.這可能是由于北方地區(qū)受到強(qiáng)大的西風(fēng)帶的影響,且其強(qiáng)度在冬季大于夏季,使得臭氧總量在季節(jié)上的波動以及空間異質(zhì)性都是冬季大于夏季.

    整體而言,中國區(qū)域50a月平均臭氧總量東北高西南低,且在春季和冬季呈現(xiàn)出明顯的由東北向西南遞減的緯向條帶狀分布.在40°N以北的東北部地區(qū),比如黑龍江,吉林以及內(nèi)蒙古東北部,多年月平均臭氧總量全國最高,且在冬季和春季明顯高于其他地區(qū),在 1~4月能夠達(dá)到400DU以上.青藏高原地區(qū)在夏季和秋季呈現(xiàn)出一個明顯的低值中心,且秋季的值低于夏季.這與肖鐘勇等[28]熊思章等[36]的研究結(jié)果相符.

    圖6描述了12個月時間變異系數(shù)的空間分布,時間變異系數(shù)是相同空間位置50個月的標(biāo)準(zhǔn)差與平均值的比值,用于描述不同空間位置在時間波動的情況.整體上,中國區(qū)域50a月平均臭氧總量在冬季波動最大,夏季最小,北方地區(qū)的波動大于南方地區(qū).冬季北方地區(qū)的臭氧總量年際波動明顯大于夏季南方地區(qū)的主要原因是北方地區(qū)受到強(qiáng)大西風(fēng)帶的影響,且其強(qiáng)度在冬季大于夏季.

    2.4 空間特征

    2.4.1 空間分布特征 圖7為中國區(qū)域臭氧總量的多年平均值、月尺度變異系數(shù)、年尺度變異系數(shù)的空間分布,多年平均值是利用50a的年平均進(jìn)行計算.從圖7可以看出,中國區(qū)域臭氧總量的50a平均值由東北向西南遞減,呈現(xiàn)很明顯的緯向條帶狀分布.這與沈凡卉等[35]對將近30年(1979~2005年)中國區(qū)域臭氧總量的空間分布的研究結(jié)果相符.在40°N以北的東北部地區(qū),臭氧總量可達(dá)360DU以上;臭氧總量在青藏高原出現(xiàn)明顯的低值區(qū),臭氧總量約為280DU,比同緯度的其他區(qū)域偏低.形成臭氧總量這一空間分布特征的主要原因是高緯度地區(qū)的臭氧層較厚,臭氧保存條件較好;且其對流層較薄,對臭氧層的影響較小;再加上臭氧由處于高勢能位的低緯度地區(qū)向處于低勢能位的高緯度地區(qū)輸送,使得高緯度地區(qū)的臭氧總量高于低緯度地區(qū).青藏高原臭氧低谷的形成,除了各種臭氧前體物排放到大氣破壞了臭氧層這一因素外,地形是主要的影響因素.首先,青藏高原區(qū)的平均高程在4000m以上,這一海拔高度縮短了臭氧氣柱,進(jìn)而減少了臭氧總量[28].再者,高海拔的高原地形所引起的動力過程也促成了其上空臭氧低谷的產(chǎn)生[48-49].最后,青藏高原高空受南亞高壓控制,有著很強(qiáng)的反氣旋環(huán)流,使得臭氧由高原向四周輻散[50],導(dǎo)致青藏高原高空形成臭氧低谷.

    圖7 50a臭氧總量的平均值、月和年尺度變異系數(shù)的空間分布

    月尺度變異系數(shù)是基于600個月的月平均,描述了在600個月中像元的波動情況,變異系數(shù)大波動大.由圖7(b)可知,月尺度變異系數(shù)在空間分布上南北梯度顯著,呈現(xiàn)條帶狀分布.在40°N以北的東北部地區(qū),月尺度變異系數(shù)最大,可達(dá)0.1以上,說明該地區(qū)臭氧總量的月際波動最大.值得注意的是,月尺度變異系數(shù)并非單一的由東北向西南遞減.東南沿海地區(qū),比如海南和臺灣兩省以及福建廣東的沿海地區(qū),月尺度變異系數(shù)略高于與其相鄰的內(nèi)陸地區(qū).可能原因是東南沿海地區(qū)受海陸風(fēng)影響較大,使得其臭氧總量的月際波動略大于相鄰的內(nèi)陸地區(qū).年尺度變異系數(shù)是基于50a的年平均,如圖7(c)所示,年尺度變異系數(shù)在空間上呈現(xiàn)由北向南遞減的條帶狀分布.在40°N以北的地區(qū),比如黑龍江、吉林、內(nèi)蒙古北部和新疆北部,年尺度變異系數(shù)最大,約達(dá)0.04以上,說明該地區(qū)臭氧總量年際波動最大.造成南北年際波動差異的主要原因可能是北方地區(qū)臭氧總量高于南方地區(qū),且受到強(qiáng)大西風(fēng)帶的影響,使得其波動更大.對比月尺度變異系數(shù)和年尺度變異系數(shù)可知,中國區(qū)域臭氧總量的月際波動大于年際波動.

    2.4.2 臭氧增長變化的空間特征 圖8(a)~圖8(d)依次展示了中國區(qū)域在1971年、1978年、1993年、2020年臭氧總量平均值的空間分布.其中,1971年和2020年是本文分析的50a時間序列的起止年份,1978年和1993年分別是這50a中臭氧總量平均值最大值和最小值出現(xiàn)的年份.圖8(e)~圖8(h)依次是中國區(qū)域在1978年與1971年、1993年與1978年、2020年與1993年、2020年與1971年臭氧總量平均值差值的空間分布,圖8(i~l)為相應(yīng)的增長率的分布情況.

    如圖8(e)和(i),1971~1978年中國區(qū)域臭氧總量平均值呈增長趨勢,且增長量和增長率都呈現(xiàn)由北向南遞減的緯向條帶狀分布.在40°N以北的相對高值地區(qū)最大增加了56DU,約為16%;而在30°N以南的相對低值地區(qū),最小增加了12DU,約為5%.增長量和增長率在相對高值地區(qū)更大的可能原因是大氣環(huán)流將大部分低緯度地區(qū)增長的臭氧輸送到更高緯度地區(qū),使得臭氧總量在更高的緯度增長得更多.

    如圖8(f)和(j),1978~1993年整個區(qū)域的臭氧總量呈減少趨勢,且減少量和減少率呈現(xiàn)由北向南遞減的緯度地帶性.在40°N以北的相對高值地區(qū)最大減少了93DU,約為22%;而在30°N以南的相對低值地區(qū),最小減少了11DU,約為4%.減少量和減少率在相對高值地區(qū)更大的可能原因是大氣環(huán)流將1991年6月菲律賓Pinatubo火山噴發(fā)產(chǎn)生的大量氣溶膠輸送到更高緯度地區(qū),造成該地區(qū)的臭氧被劇烈消耗.

    圖8 (a)1971、(b)1978、(c)1993、(d)2020臭氧總量平均值及(e)1978~1971、(f)1993~1978、(g)2020~1993、(h)2020~1971臭氧總量平均值差值及(i)1978~1971增長率、(j)1993~1978增長率、(k)2020~1993增長率、(l)2020~1971增長率

    如圖8(g)和(k),1993~2020年的變化特征與上述兩個時段不同,無明顯的緯度地帶性.臭氧總量在整個區(qū)域都是增長的,且在西北地區(qū)出現(xiàn)最大增長,增長量為18DU,約為6%;東南地區(qū)出現(xiàn)最小增長,增長量為4DU,約為1%.增長量和增長率在西北地區(qū)更大的原因可能是我國于2000年開始加強(qiáng)西部開發(fā),大力開展的經(jīng)濟(jì)建設(shè)、工業(yè)生產(chǎn)等活動使得該地區(qū)人為排放的臭氧增多.而我國東南地區(qū)經(jīng)濟(jì)較為發(fā)達(dá),雖然工業(yè)生產(chǎn)排放的臭氧增加,但臭氧前體物(如汽車排放的氮氧化合物)也增加,二者相互消耗抵消,使得東南地區(qū)的臭氧增長量和增長率均無西北地區(qū)大.

    1971~2020年,中國區(qū)域臭氧總量平均值整體呈現(xiàn)輕微下降趨勢(圖3),但在空間上的變化卻不同.如圖8(h)和(l)所示,臭氧總量的變化率和變化率均呈現(xiàn)緯度地帶性.但在40°N以北的相對高值地區(qū),臭氧總量平均值在這一時段呈減少趨勢;最大減少了27DU,約為8%.而在30°N以南的相對低值地區(qū),臭氧總量平均值呈增長趨勢;最大增長了14DU,約為5%.在這50a的時間序列中,臭氧總量平均值在相對高值地區(qū)減少,而在相對低值地區(qū)增長的主要原因是大氣環(huán)流將1991年6月菲律賓Pinatubo火山噴發(fā)產(chǎn)生的大量氣溶膠輸送到相對高值地區(qū),造成該地區(qū)的臭氧被劇烈消耗.雖然《蒙特利爾議定書》及其修正案中有關(guān)限制臭氧消耗物質(zhì)生產(chǎn)和使用的規(guī)定起了成效,但是還不足以彌補(bǔ)相對高值地區(qū)之前的劇烈消耗;卻補(bǔ)足了相對低值地區(qū)被消耗的臭氧,還使其獲得增長.另外,我國在2007年禁止生產(chǎn)氯氟烴類[51],但是仍有庫存,同年在東北地區(qū)仍檢測到一定程度的源自庫存物質(zhì)的氯氟烴類的排放[52].這進(jìn)一步導(dǎo)致了臭氧總量平均值在40°N以北的相對高值地區(qū)出現(xiàn)最大減少.

    3 結(jié)論

    3.1 MSR-2資料衛(wèi)星反演的臭氧總量與地面觀測臭氧總量的相關(guān)性強(qiáng),適合應(yīng)用于中國區(qū)域臭氧總量時空變化特征的研究.

    3.2 50a來中國區(qū)域的臭氧總量呈現(xiàn)輕微的下降趨勢.年平均臭氧總量在1978年和1993年分別出現(xiàn)最大值(347.5±53.8)DU和最小值(291.9±29.5)DU.在1971~1978年、1978~1993年、1993~2020年這3個時段臭氧總量平均值在整個中國區(qū)域的變化呈單向,分別是增長、減少、增長.月平均臭氧總量隨季節(jié)變化呈現(xiàn)正弦曲線形態(tài).峰值和谷值分別出現(xiàn)在3月和10月,約為338和285DU,相差53DU.

    3.3 中國區(qū)域臭氧總量的50a平均值由東北向西南遞減,呈現(xiàn)緯向條帶狀分布.在40°N以北的東北部地區(qū),該值可達(dá)360DU以上.中國區(qū)域50a月平均臭氧總量同樣呈現(xiàn)緯向條帶狀分布.此外,時間變異系數(shù)和空間變異系數(shù)隨季節(jié)的變化規(guī)律相似,夏季最小,接著依次是秋季和春季,冬季最大,即臭氧總量的變化和空間差異在夏季均最小.

    3.4 50a期間,不同時段、不同區(qū)域臭氧總量的變化趨勢各不相同.在1971~1978年,臭氧總量的增長量和增長率都呈現(xiàn)由北向南遞減的緯向條帶狀分布.在40°N以北的相對高值地區(qū)最大增加了56DU,約為16%;而在30°N以南的相對低值地區(qū),最小增加了12DU,約為5%.在1978~1993年,減少量和減少率也呈現(xiàn)由北向南遞減的緯度地帶性.在40°N以北的相對高值地區(qū)最大減少了93DU,約為22%;而在30°N以南的相對低值地區(qū),最小減少了11DU,約為4%.在1993~2020年,西北地區(qū)出現(xiàn)最大增長,增長量為18DU,約為6%;東南地區(qū)出現(xiàn)最小增長,增長量為4DU,約為1%.

    [1] 曹曉云,祁棟林,肖建設(shè),等.基于衛(wèi)星觀測的青海高原對流層臭氧時空分布特征研究 [J]. 環(huán)境科學(xué)學(xué)報, 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.

    [2] Chen L, Yu B, Chen Z, et al. Investigating the temporal and spatial variability of total ozone column in the Yangtze River Delta using satellite data: 1978-2013 [J]. Remote Sensing, 2014,6(12):12527- 12543.

    [3] Bencherif H, Toihir A M, Mbatha N, et al. Ozone variability and trend estimates from 20-years of ground-based and satellite observations at Irene Station, South Africa [J]. Atmosphere, 2020,11(11):1216.

    [4] Al-Salihi A, Hassan Z. Temporal and spatial variability and trend investigation of total ozone column over Iraq employing remote sensing data: 1979-2012 [J]. International Letters of Chemistry, Physics and Astronomy, 2015,53:1-18.

    [5] Li X, Li Y, Ma Q, et al. Spatial and temporal variation of total and tropospheric ozone columns over China [J]. Int. Arch. Photogramm. Remote Sens. Spatial Inf. Sci., 2019,XLII-3/W9:113-118.

    [6] Andersen S O, Halberstadt M L, Borgford-Parnell N. Stratospheric ozone, global warming, and the principle of unintended consequences- An ongoing science and policy success story [J]. Journal of the Air & Waste Management Association, 2013,63(6):607-647.

    [7] 馬明亮.對流層臭氧時空分析影響因素研究及近地面臭氧估算 [J]. 測繪學(xué)報, 2020,49(11):1507.

    Ma M L. Temporal and spatial analysis of tropospheric ozone, simulation of influencing factors and surface ozone estimation [J]. Acta Geodaetica et Cartographica Sinica, 2020,49(11):1507.

    [8] 閆建軍,金蓮姬,王革麗.平流層臭氧非線性響應(yīng)的機(jī)理研究 [J]. 氣候與環(huán)境研究, 2012,17(5):639-645.

    Yan J J, Jin L J, Wang G L. A Study of the mechanism of nonlinear responses in stratospheric ozone [J]. Climatic and Environmental Research (in Chinese), 2012,17(5):639-645.

    [9] 趙川鴻,趙鵬國,周筠珺.西南地區(qū)臭氧空間分布及變化趨勢 [J]. 氣象科學(xué), 2018,38(2):149-156.

    Zhao C H, Zhao P G, Zhou Y J. The temporal and spatial distributions of the total ozone in Southwest China [J]. Journal of the Meteorological Sciences, 2018,38(2):149-156.

    [10] 呂達(dá)仁,卞建春,陳洪濱,等.平流層大氣過程研究的前沿與重要性 [J]. 地球科學(xué)進(jìn)展, 2009,24(3):221-228.

    LV D L, Bian J C, Chen H B, et al. Frontiers and significance of research on stratospheric processes [J]. Advances in Earth Science, 2009,24(3):221-228.

    [11] 張秀年,常有禮,陳 輝,等.中國低緯度地區(qū)大氣臭氧總量的變化特征 [J]. 熱帶氣象學(xué)報, 2004,(3):321-327.

    Zhang X N, Chang Y L, Chen H, et al. Variation of total ozone in lower latitude area of China [J]. Journal of Tropical Meteorology, 2004, (3):321-327.

    [12] Ainsworth E A. Understanding and improving global crop response to ozone pollution [J]. The Plant Journal, 2017,90(5):886-897.

    [13] Nouchi I, Ito O, Harazono Y, et al. Effects of chronic ozone exposure on growth, root respiration and nutrient uptake of rice plants [J]. Environmental Pollution, 1991,74(2):149-164.

    [14] Cooley D RManning W J. The impact of ozone on assimilate partitioning in plants: A review [J]. Environmental Pollution, 1987, 47(2):95-113.

    [15] Nuvolone D, Petri DVoller F.The effects of ozone on human health [J]. Environmental Science and Pollution Research, 2018,25(9):8074- 8088.

    [16] 王 萍,劉 濤,楊國林,等.中國主要城市臭氧濃度的時空變化特征 [J]. 遙感信息, 2019,34(4):121-127.

    Wang P, Liu T, Yang G L, et al. Spatial and temporal variation of ozone concentration in major cities of China [J]. Remote Sensing Information, 2019,34(4):121-127.

    [17] 鄭向東,韋小麗.中國4個地點(diǎn)地基與衛(wèi)星臭氧總量長期觀測比較 [J]. 應(yīng)用氣象學(xué)報, 2010,21(1):1-10.

    Zheng X D, Wei X L. Long-term total ozone comparisons between space-based and ground-based observations at 4sites in China [J]. Journal of Applied Meteorological Science, 2010,21(1):1-10.

    [18] Kuttippurath J, Kumar P, Nair P J, et al. Accuracy of satellite total column ozone measurements in polar vortex conditions: Comparison with ground-based observations in 1979~2013 [J]. Remote Sensing of Environment, 2018,209:648-659.

    [19] Nair P J, Froidevaux L, Kuttippurath J, et al. Subtropical and midlatitude ozone trends in the stratosphere: Implications for recovery [J]. Journal of Geophysical Research: Atmospheres, 2015,120(14): 7247-7257.

    [20] Peng X, Shen H, Zhang L, et al. Spatially continuous mapping of daily global ozone distribution (2004~2014) with the Aura OMI sensor [J]. Journal of Geophysical Research: Atmospheres, 2016,121(21):12,702- 12,722.

    [21] Loyola D G, Coldewey-Egbers R M, Dameris M, et al. Global long-term monitoring of the ozone layer – a prerequisite for predictions [J]. International Journal of Remote Sensing, 2009,30 (15/16):4295-4318.

    [22] Orphal J, Staehelin J, Tamminen J, et al. Absorption cross-sections of ozone in the ultraviolet and visible spectral regions: Status report 2015 [J]. Journal of Molecular Spectroscopy, 2016,327:105-121.

    [23] 劉姝媛,朱 彬,蘇繼峰,等.近30a北極平流層臭氧的季節(jié)和年際變化特征 [J]. 大氣科學(xué)學(xué)報, 2011,34(5):538-546.

    Liu S Y, Zhu B, Su J F, et al. Seasonal and interannual characteristics of ozone over Arctic stratosphere in recent 30years [J]. Trans Atmos Sci, 2011,34(5):538-546.

    [24] Van Der A R J, Allaart M A F, Eskes H J. Extended and refined multi sensor reanalysis of total ozone for the period 1970~2012 [J]. Atmos. Meas. Tech., 2015,8(7):3021-3035.

    [25] Bai K, Liu C, Shi R, et al. Global validation of FY-3A total ozone unit (TOU) total ozone columns using ground-based Brewer and Dobson measurements [J]. International Journal of Remote Sensing, 2013, 34(14):5228-5242.

    [26] Ma M, Shi R, Gao W.Validation of the OMI-TOMS and OMI-DOAS total ozone column data using ground-based observations over China [C]//Proceedings of SPIE Volume 9610, Remote Sensing and Modeling of Ecosystems for Sustainability XII. San Diego, California, United States, 2015.

    [27] 張 瑩,高 玚,祝善友,等.近30a中國上空臭氧總量時空變化遙感監(jiān)測分析 [J]. 地球信息科學(xué)學(xué)報, 2014,16(6):971-978.

    Zhang Y, Gao Y, Zhu S Y, et al. Variation of total ozone over China for 30years analyzed by multi-source satellite remote sensing data [J]. Journal of Geo-information Science, 2014,16(6):971-978.

    [28] 肖鐘湧,江 洪.利用遙感監(jiān)測青藏高原上空臭氧總量30a的變化 [J]. 環(huán)境科學(xué), 2010,31(11):2569-2574.

    Xiao Z Y, Jiang H. Variation of total ozone over Tibetain Plateau during 30years monitored by remotely sensed data [J]. Environmental Science, 2010,31(11):2569-2574.

    [29] Manzer L E.The CFC-Ozone Issue: Progress on the development of alternatives to CFCs [J]. Science, 1990,249(4964):31-35.

    [30] Lovelock J E. Atmospheric fluorine compounds as indicators of air movements [J]. Nature, 1971,230(5293):379-379.

    [31] Molina M J, Rowland F S. Stratospheric sink for chlorofluoromethanes: Chlorine atom-catalysed destruction of ozone [J]. Nature, 1974,249(5460):810-812.

    [32] WMO. Scientific assessment of ozone depletion: 2006 [R]. Global Ozone Research and Monitoring Project-Report No. 50,2006.

    [33] WMO. Scientific assessment of ozone depletion: 2018 [R]. Global Ozone Research and Monitoring Project-Report No. 58,2018.

    [34] 韋惠紅,鄭有飛.我國臭氧總量的時空分布特征 [J]. 南京氣象學(xué)院學(xué)報, 2006,(3):390-395.

    Wei H H, Zheng Y F. Analysis of the temporal and spatial distributions of the total ozone over China [J]. Journal of Nanjing Institute of Meteorology, 2006,(3):390-395.

    [35] 沈凡卉,王體健,沈 毅,等.中國近30a臭氧氣候場特征 [J]. 大氣科學(xué)學(xué)報, 2011,34(3):288-296.

    Shen F H, Wang T J, Shen Y, et al. Characteristics of ozone in China in recent thirty years [J]. Trans. Atmos. Sci., 2011,34(3):288-296.

    [36] 熊思章,陳權(quán)亮.青藏高原上空臭氧的時空演變特征 [J]. 成都信息工程大學(xué)學(xué)報, 2020,35(6):671-677.

    Xiong S Z, Chen Q L. Temporal and spatial evolution of ozone over the Tibetan Plateau [J]. Journal of Chengdu University of Information Technology, 2020,35(6):671-677.

    [37] 姚 青,馬志強(qiáng),郝天依,等.京津冀區(qū)域臭氧時空分布特征及其背景濃度估算 [J]. 中國環(huán)境科學(xué), 2021,41(11):4999-5008.

    Yao Q, Ma Z Q, Hao T Y, et al. Temporal and spatial distribution characteristics and background concentration estimation of ozone in Beijing-Tianjin-Hebei region [J]. China Environmental Science, 2021,41(11):4999-5008.

    [38] 王金南,董戰(zhàn)峰,蔣洪強(qiáng),等.中國環(huán)境保護(hù)戰(zhàn)略政策70年歷史變遷與改革方向 [J]. 環(huán)境科學(xué)研究, 2019,32(10):1636-1644.

    Wang J N, Dong Z F, Jiang H Q, et al. Historical evolution and reform of China’s environmental strategy and policy during the past seventy years (1949~2019) [J]. Research of Environmental Sciences, 2019,32(10):1636-1644.

    [39] 王 鳳,張劍波,馮金敏,等.中國歷史上和未來CFC-11和CFC-12的排放量計算 [J]. 環(huán)境科學(xué)學(xué)報, 2010,30(9):1758-1765.

    Wang F, Zhang J B, Feng J M, et al. Estimated historical and future emissions of CFC-11and CFC-12 in China [J]. Acta Scientiae Circumstantiae, 2010,30(9):1758-1765.

    [40] Hofmann D J, Oltmans S J, Komhyr W D, et al. Ozone loss in the lower stratosphere over the United States in 1992~1993: Evidence for heterogeneous chemistry on the Pinatubo aerosol [J]. Geophysical Research Letters, 1994,21(1):65-68.

    [41] Solomon S, Portmann R W, Garcia R R, et al. Ozone depletion at mid-latitudes: Coupling of volcanic aerosols and temperature variability to anthropogenic chlorine [J]. Geophysical Research Letters, 1998,25(11):1871-1874.

    [42] Portmann R W, Solomon S, Garcia R R, et al. Role of aerosol variations in anthropogenic ozone depletion in the polar regions [J]. Journal of Geophysical Research: Atmospheres, 1996,101(D17): 22991-23006.

    [43] WMO. Scientific assessment of ozone depletion: 2002 [R]. Global Ozone Research and Monitoring Project-Report No. 47,2002.

    [44] Chen D, Nunez M. Temporal and spatial variability of total ozone in Southwest Sweden revealed by two ground-based instruments [J]. International Journal of Climatology, 1998,18(11):1237-1246.

    [45] Tung K K, Yang H. Dynamic variability of column ozone [J]. Journal of Geophysical Research: Atmospheres, 1988,93(D9):11123-11128.

    [46] Aesawy A M, Mayhoub A B, Sharobim W M. Seasonal variation of photochemical and dynamical components of ozone in subtropical regions [J]. Theoretical and Applied Climatology, 1994,49(4):241- 247.

    [47] Antón M, Bortoli D, Costa M J, et al. Temporal and spatial variabilities of total ozone column over Portugal [J]. Remote Sensing of Environment, 2011,115(3):855-863.

    [48] 劉 煜,李維亮.青藏高原臭氧低谷的加深及其可能的影響 [J]. 氣象學(xué)報, 2001,(1):97-106.

    Liu Y, Li W L. Deepening of ozone valley over Tibetan Plateau and its possible influences [J]. Acta Meteorologica Sinica, 2001,(1):97-106.

    [49] 劉 煜,李維亮,周秀驥.青藏高原臭氧變化趨勢的預(yù)測 [J]. 中國科學(xué)(D輯:地球科學(xué)), 2001,(S1):308-311.

    Liu Y, Li W L, Zhou X J. Prediction of ozone variation trend over Tibetain Plateau [J]. Science in China (Series D), 2001,(S1):308-311.

    [50] 王喜紅,李興生.中國地區(qū)對流層臭氧變化和分布的數(shù)值研究 [J]. 氣象學(xué)報, 1998,(3):78-93.

    Wang X H, Li X S. A numerical study of the variations and distributions of tropospheric ozone and its precursors over China [J]. Acta Meteorologica Sinica, 1998,(3):78-93.

    [51] Wan D, Xu J, Zhang J, et al. Historical and projected emissions of major halocarbons in China [J]. Atmospheric Environment, 2009, 43(36):5822-5829.

    [52] Xue L, Wang T, Simpson I J, et al. Vertical distributions of non-methane hydrocarbons and halocarbons in the lower troposphere over northeast China [J]. Atmospheric Environment, 2011,45(36): 6501-6509.

    致謝:本研究數(shù)據(jù)來自美國戈達(dá)地球科學(xué)數(shù)據(jù)和信息服務(wù)中心(Goddard Earth Sciences Data and Information Services Center),在此表示感謝!

    Investigating the temporal and spatial variabilities of total ozone over China for the past 50 years.

    XIE Jing-han, LI Sa, XIAO Zhong-yong*

    (College of Harbour and Coastal Engineering, Jimei University, Xiamen 361021, China)., 2022,42(7):2977~2987

    Based on the Multi-Sensor Reanalysis version 2data, which was validated by the data from 6ground-based observation stations belonging to the World Ozone and Ultraviolet Radiation Data Centre, this study analyzed the temporal and spatial variabilities of total ozone over China for the period 1971~2020. It was found that the total ozone over China showed a slight downward trend in the past 50years. The annual mean total ozone reached a maximum value of (347.5 ± 53.8) DU in 1978 and a minimum one of (291.9 ± 29.5) DU in 1993. In the three sub-periods, 1971~1978, 1978~1993, and 1993~2020, the annual mean total ozone over China completely increased, decreased, and increased, respectively. The seasonal variability of the monthly mean total ozone presented a sinusoidal wave with a maximum value of about 338DU in March and a minimum value of about 285DU in October. The spatial distribution of the annual mean total ozone over China in the past 50years showed a latitudinal pattern by decreasing gradually from northeast areas to southwest areas. Moreover, the values could be greater than 360DU in the northeast China with latitudes over 40°N. The spatial distribution of monthly mean total ozone also showed a latitudinal pattern. In addition, the temporal and spatial variation coefficients showed similar seasonal variability, with the minimum value in summer, followed by the values in autumn, spring, and winter. Namely, the variability of total ozone and its spatial difference was minimum in summer. Finally, the trend of total ozone varied with periods and areas in the past 50years. During the period 1971~1978, the increase value and rate of annual mean total ozone showed a latitudinal pattern by decreasing gradually from north areas to south areas. In the areas with relatively high annual mean total ozone, such as the north China with latitude over 40°N, the maximum increase value was 56DU, around 16%. Yet, in the areas with relatively low annual mean total ozone, such as the south China with latitude below 30°N, the minimum increase value was 12DU, around 5%. During the period 1978~1993, the decreased value and rate of annual mean total ozone also showed a latitudinal pattern. In the areas with relatively high annual mean total ozone, such as the north China with latitudes over 40°N, the maximum decrease value was 93DU, around 22%. Yet, in the areas with relatively low annual mean total ozone, such as the south China with latitude below 30°N, the minimum decrease value was 11DU, around 4%. During the period 1993~2020, the maximum increase value (18DU, about 6%) was found in the northwest areas. In contrast, the minimum increase value (4DU, about 1%) was found in the southeast areas.

    total ozone;temporal variability;seasonal variability;spatial characteristics;China

    X511

    A

    1000-6923(2022)07-2977-11

    謝靜晗(1988-),女,福建廈門人,講師,博士,主要從事大氣環(huán)境遙感研究.發(fā)表論文3篇.

    2021-12-01

    集美大學(xué)國家基金培育計劃項(xiàng)目資助(ZP2021019);福建省自然科學(xué)基金計劃資助項(xiàng)目(2022J01817);福建省自然科學(xué)基金計劃資助項(xiàng)目(2021J01839);福建省中青年教師教育科研項(xiàng)目(JAT200282);集美大學(xué)科研基金(ZQ2019025)

    * 責(zé)任作者, 副教授, zyxiao@jmu.edu.cn

    猜你喜歡
    臭氧總量平均值
    “平均值代換”法在數(shù)學(xué)解題中的應(yīng)用
    “十三五”期間山西省與10省簽約糧食總量2230萬噸
    文印室內(nèi)臭氧散發(fā)實(shí)測分析及模擬
    2020年全國農(nóng)民工總量比上年減少517萬人
    為何化肥淡儲總量再度增加
    看不見的污染源——臭氧
    利用臭氧水防治韭菜遲眼蕈蚊
    中國蔬菜(2016年8期)2017-01-15 14:23:43
    臭氧分子如是說
    平面圖形中構(gòu)造調(diào)和平均值幾例
    基于電流平均值的改進(jìn)無功檢測法
    電測與儀表(2014年6期)2014-04-04 11:59:46
    高清视频免费观看一区二区| 中文字幕制服av| 欧美三级亚洲精品| 国产在线一区二区三区精| 一二三四中文在线观看免费高清| 亚洲综合色惰| 亚洲欧美一区二区三区黑人 | 美女脱内裤让男人舔精品视频| 免费人成在线观看视频色| 能在线免费看毛片的网站| 各种免费的搞黄视频| 九色成人免费人妻av| 国产成人av激情在线播放 | 毛片一级片免费看久久久久| 午夜激情久久久久久久| 久久久精品94久久精品| xxxhd国产人妻xxx| 免费观看性生交大片5| 99九九在线精品视频| 国产精品久久久久久精品电影小说| 天堂俺去俺来也www色官网| 精品人妻熟女av久视频| 亚洲av综合色区一区| 国产日韩欧美亚洲二区| 久热这里只有精品99| 亚洲国产成人一精品久久久| 日日摸夜夜添夜夜添av毛片| 国产精品偷伦视频观看了| 国产亚洲午夜精品一区二区久久| 免费黄色在线免费观看| 91精品一卡2卡3卡4卡| 母亲3免费完整高清在线观看 | 日韩熟女老妇一区二区性免费视频| 一区在线观看完整版| 亚洲av福利一区| 婷婷色综合www| 99久久人妻综合| 日韩熟女老妇一区二区性免费视频| 久久精品国产自在天天线| 人人妻人人添人人爽欧美一区卜| 国产欧美日韩一区二区三区在线 | 99九九线精品视频在线观看视频| 国产乱来视频区| 少妇被粗大的猛进出69影院 | 中文欧美无线码| 久久午夜综合久久蜜桃| 亚洲av不卡在线观看| 99re6热这里在线精品视频| 汤姆久久久久久久影院中文字幕| 在线观看美女被高潮喷水网站| 日韩一区二区三区影片| 十八禁网站网址无遮挡| 美女主播在线视频| 性色avwww在线观看| av线在线观看网站| 天堂中文最新版在线下载| 一级毛片 在线播放| 日韩免费高清中文字幕av| 亚洲少妇的诱惑av| 久久99一区二区三区| 精品久久蜜臀av无| 欧美xxxx性猛交bbbb| 久热这里只有精品99| 亚洲精品乱码久久久久久按摩| 777米奇影视久久| 精品久久久精品久久久| 99热这里只有精品一区| 熟女电影av网| 亚洲国产欧美在线一区| 在线免费观看不下载黄p国产| 日韩一区二区视频免费看| 交换朋友夫妻互换小说| 亚洲av电影在线观看一区二区三区| 99热6这里只有精品| 国产精品国产av在线观看| 男女边摸边吃奶| 久久久精品94久久精品| 久久99热6这里只有精品| 尾随美女入室| 成人国产麻豆网| 18禁在线无遮挡免费观看视频| 久久精品国产亚洲av天美| 老司机影院毛片| www.色视频.com| 丰满少妇做爰视频| 成年女人在线观看亚洲视频| av播播在线观看一区| 黄片播放在线免费| 高清在线视频一区二区三区| 不卡视频在线观看欧美| 91国产中文字幕| 亚洲情色 制服丝袜| 秋霞在线观看毛片| 亚洲精品美女久久av网站| 一区在线观看完整版| 久久97久久精品| 在线观看www视频免费| 熟女av电影| 汤姆久久久久久久影院中文字幕| 青青草视频在线视频观看| 蜜桃国产av成人99| 大陆偷拍与自拍| av免费在线看不卡| 亚洲欧美一区二区三区国产| 亚洲国产日韩一区二区| 桃花免费在线播放| 国产精品99久久99久久久不卡 | 国产日韩一区二区三区精品不卡 | 日本-黄色视频高清免费观看| av又黄又爽大尺度在线免费看| 最后的刺客免费高清国语| 国产乱人偷精品视频| 人人妻人人澡人人爽人人夜夜| 久久久久久伊人网av| 高清欧美精品videossex| 欧美97在线视频| 国产一区二区三区综合在线观看 | 亚洲国产精品成人久久小说| 尾随美女入室| 久久女婷五月综合色啪小说| 亚洲美女视频黄频| 中文字幕av电影在线播放| 国产乱来视频区| 日韩视频在线欧美| 国产成人精品福利久久| 亚洲一区二区三区欧美精品| 9色porny在线观看| 久久久久精品性色| 高清毛片免费看| 女性生殖器流出的白浆| 美女主播在线视频| 多毛熟女@视频| 中文乱码字字幕精品一区二区三区| 亚洲av在线观看美女高潮| 免费av不卡在线播放| 亚洲国产精品专区欧美| 免费观看无遮挡的男女| 久久精品久久久久久久性| 久久精品国产a三级三级三级| av专区在线播放| 狠狠婷婷综合久久久久久88av| 观看av在线不卡| 久久青草综合色| 免费看av在线观看网站| 在线观看人妻少妇| 成人亚洲欧美一区二区av| 日韩熟女老妇一区二区性免费视频| 亚洲精品视频女| 中文精品一卡2卡3卡4更新| 久久人人爽av亚洲精品天堂| 午夜激情av网站| 青春草视频在线免费观看| 亚洲欧美一区二区三区黑人 | 伊人久久国产一区二区| 国产毛片在线视频| 国产男女超爽视频在线观看| 欧美人与性动交α欧美精品济南到 | 国产成人精品在线电影| 日产精品乱码卡一卡2卡三| 亚洲综合色惰| av国产精品久久久久影院| 韩国高清视频一区二区三区| 欧美变态另类bdsm刘玥| 亚洲精品一区蜜桃| 女的被弄到高潮叫床怎么办| 日日啪夜夜爽| videos熟女内射| 亚洲经典国产精华液单| 国产伦精品一区二区三区视频9| 色视频在线一区二区三区| 欧美xxxx性猛交bbbb| 成人午夜精彩视频在线观看| 亚洲欧美成人精品一区二区| 国产黄色免费在线视频| 在线观看国产h片| 国产视频内射| 国产探花极品一区二区| 69精品国产乱码久久久| 国产精品久久久久久av不卡| 午夜精品国产一区二区电影| 伊人久久精品亚洲午夜| 婷婷色综合大香蕉| 国产黄片视频在线免费观看| 婷婷成人精品国产| 欧美人与善性xxx| 男人操女人黄网站| 欧美日韩视频高清一区二区三区二| 免费看av在线观看网站| 久久久久久久久久成人| 免费人妻精品一区二区三区视频| 国产国拍精品亚洲av在线观看| 水蜜桃什么品种好| 少妇的逼水好多| av在线观看视频网站免费| 国产视频首页在线观看| 国产日韩欧美视频二区| 精品人妻在线不人妻| 日韩av在线免费看完整版不卡| 午夜久久久在线观看| 黄色视频在线播放观看不卡| 久久这里有精品视频免费| 男女边摸边吃奶| 亚洲av电影在线观看一区二区三区| 久久99热6这里只有精品| 麻豆成人av视频| 哪个播放器可以免费观看大片| 久久久久久伊人网av| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩国产mv在线观看视频| videos熟女内射| 成人毛片a级毛片在线播放| xxx大片免费视频| av国产精品久久久久影院| 黑人巨大精品欧美一区二区蜜桃 | 另类亚洲欧美激情| 亚洲av中文av极速乱| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品美女久久av网站| 最黄视频免费看| 久久97久久精品| 在线亚洲精品国产二区图片欧美 | 亚洲色图综合在线观看| 亚洲欧美中文字幕日韩二区| 国产日韩欧美在线精品| 亚洲成人av在线免费| 久久久久精品久久久久真实原创| 91久久精品国产一区二区三区| 老司机影院成人| 大又大粗又爽又黄少妇毛片口| 日本黄大片高清| 久久久久久伊人网av| 亚洲色图 男人天堂 中文字幕 | 久久婷婷青草| av.在线天堂| 日日撸夜夜添| 性色avwww在线观看| 国产日韩欧美亚洲二区| 国产免费现黄频在线看| 欧美人与性动交α欧美精品济南到 | 亚洲精品日本国产第一区| 午夜影院在线不卡| 中文天堂在线官网| 午夜久久久在线观看| 欧美日韩av久久| h视频一区二区三区| 一级,二级,三级黄色视频| 欧美精品国产亚洲| 国产老妇伦熟女老妇高清| 亚洲精品久久成人aⅴ小说 | 精品久久久久久久久av| av女优亚洲男人天堂| 日韩精品免费视频一区二区三区 | 国产成人一区二区在线| 国产精品一国产av| 99九九线精品视频在线观看视频| 蜜桃国产av成人99| 日韩,欧美,国产一区二区三区| 国产午夜精品久久久久久一区二区三区| 久久97久久精品| 黑人巨大精品欧美一区二区蜜桃 | 亚洲一区二区三区欧美精品| 国产精品人妻久久久影院| 亚洲情色 制服丝袜| 妹子高潮喷水视频| 中文字幕免费在线视频6| 97超碰精品成人国产| 成人国产麻豆网| 色婷婷久久久亚洲欧美| 热99久久久久精品小说推荐| 亚洲av中文av极速乱| 热re99久久国产66热| 成人综合一区亚洲| 纵有疾风起免费观看全集完整版| 精品久久久久久久久av| 国产视频内射| 久久久久久久久久人人人人人人| 国产精品国产三级国产av玫瑰| 久久久久网色| 热99久久久久精品小说推荐| 观看av在线不卡| 美女福利国产在线| 亚洲国产精品国产精品| 国产高清不卡午夜福利| 桃花免费在线播放| av电影中文网址| 亚洲精品一二三| 91国产中文字幕| 大片免费播放器 马上看| 精品一区二区三卡| 亚洲一区二区三区欧美精品| 最近中文字幕高清免费大全6| 新久久久久国产一级毛片| 狠狠婷婷综合久久久久久88av| 国产亚洲欧美精品永久| 麻豆乱淫一区二区| kizo精华| 日本与韩国留学比较| av国产精品久久久久影院| 亚洲少妇的诱惑av| 久久久久久久久久久久大奶| 黄片播放在线免费| 女的被弄到高潮叫床怎么办| 久久久久久久久久人人人人人人| xxx大片免费视频| 亚洲欧美日韩另类电影网站| 热99久久久久精品小说推荐| 最黄视频免费看| 国产免费一区二区三区四区乱码| 国产日韩欧美在线精品| 日产精品乱码卡一卡2卡三| 精品少妇久久久久久888优播| 亚洲av在线观看美女高潮| 寂寞人妻少妇视频99o| √禁漫天堂资源中文www| 欧美亚洲日本最大视频资源| 免费高清在线观看日韩| 18禁裸乳无遮挡动漫免费视频| 久久国内精品自在自线图片| 精品一区二区三卡| 亚洲精品自拍成人| 最近中文字幕高清免费大全6| 久久久久久久久久久丰满| 免费观看在线日韩| 国产亚洲午夜精品一区二区久久| 免费观看的影片在线观看| 成年女人在线观看亚洲视频| 大香蕉久久成人网| 久久青草综合色| 亚洲国产精品国产精品| 两个人的视频大全免费| 成人18禁高潮啪啪吃奶动态图 | 国产黄片视频在线免费观看| 91成人精品电影| 一级,二级,三级黄色视频| 黄色毛片三级朝国网站| 精品视频人人做人人爽| 精品国产乱码久久久久久小说| 亚洲精品久久午夜乱码| 国产乱人偷精品视频| 亚洲一区二区三区欧美精品| 国产成人免费无遮挡视频| 日日啪夜夜爽| 国产永久视频网站| 久久国产亚洲av麻豆专区| 成人午夜精彩视频在线观看| 国产国语露脸激情在线看| av电影中文网址| 人人妻人人爽人人添夜夜欢视频| 精品久久久久久电影网| 色婷婷久久久亚洲欧美| 亚洲无线观看免费| 久久久久久久久大av| 人妻系列 视频| 免费av不卡在线播放| 黑人猛操日本美女一级片| 国产精品蜜桃在线观看| 国产免费现黄频在线看| 精品国产一区二区久久| 免费大片18禁| 国产日韩欧美亚洲二区| 美女大奶头黄色视频| 亚洲精品久久午夜乱码| 丝袜喷水一区| 婷婷色麻豆天堂久久| 色视频在线一区二区三区| 午夜福利视频在线观看免费| 午夜福利网站1000一区二区三区| 少妇的逼好多水| 国产精品久久久久久精品古装| 亚洲精品美女久久av网站| 极品人妻少妇av视频| 国产视频内射| 一级毛片黄色毛片免费观看视频| 亚洲一级一片aⅴ在线观看| 亚洲精品日本国产第一区| 久久久久国产网址| 日韩中文字幕视频在线看片| 天堂俺去俺来也www色官网| 色94色欧美一区二区| 亚洲内射少妇av| 免费黄网站久久成人精品| 男女啪啪激烈高潮av片| 少妇猛男粗大的猛烈进出视频| 国产色婷婷99| 亚洲一区二区三区欧美精品| 91在线精品国自产拍蜜月| 伊人久久国产一区二区| 永久免费av网站大全| 色视频在线一区二区三区| 日日啪夜夜爽| 久久99一区二区三区| 久久久久人妻精品一区果冻| 中文字幕制服av| 欧美日韩在线观看h| 久久精品久久精品一区二区三区| 免费少妇av软件| 免费观看a级毛片全部| a 毛片基地| 美女脱内裤让男人舔精品视频| 人人妻人人添人人爽欧美一区卜| 久久国产精品男人的天堂亚洲 | av免费在线看不卡| 欧美 日韩 精品 国产| 亚洲激情五月婷婷啪啪| 中国国产av一级| 亚洲精品第二区| .国产精品久久| 91久久精品国产一区二区三区| 国产免费又黄又爽又色| 亚洲国产av影院在线观看| 韩国高清视频一区二区三区| 久久99一区二区三区| 亚洲综合精品二区| 黄色配什么色好看| 激情五月婷婷亚洲| 午夜影院在线不卡| 视频区图区小说| 免费看av在线观看网站| 国产精品 国内视频| 国产黄片视频在线免费观看| 国产爽快片一区二区三区| 免费看光身美女| 国国产精品蜜臀av免费| 亚洲av在线观看美女高潮| 国产精品一国产av| 国产精品一区www在线观看| 国产高清不卡午夜福利| 亚洲国产精品一区三区| av专区在线播放| 又粗又硬又长又爽又黄的视频| 精品99又大又爽又粗少妇毛片| 免费高清在线观看视频在线观看| 国产 精品1| 亚洲欧美一区二区三区黑人 | 久久久午夜欧美精品| 国产精品熟女久久久久浪| 亚洲人成77777在线视频| 日韩成人av中文字幕在线观看| 18+在线观看网站| 欧美国产精品一级二级三级| 青春草视频在线免费观看| 日韩大片免费观看网站| 美女视频免费永久观看网站| 亚洲激情五月婷婷啪啪| 在现免费观看毛片| 肉色欧美久久久久久久蜜桃| 能在线免费看毛片的网站| 精品酒店卫生间| 国精品久久久久久国模美| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久精品区二区三区| 国产视频首页在线观看| 亚洲国产毛片av蜜桃av| 亚洲综合色惰| 午夜精品国产一区二区电影| 最新的欧美精品一区二区| 色哟哟·www| 91成人精品电影| 免费播放大片免费观看视频在线观看| 麻豆乱淫一区二区| 人成视频在线观看免费观看| 午夜91福利影院| 日韩熟女老妇一区二区性免费视频| 国产综合精华液| 久久毛片免费看一区二区三区| 国产男女内射视频| 丁香六月天网| 国产极品粉嫩免费观看在线 | 亚洲精华国产精华液的使用体验| 高清av免费在线| 亚洲欧美色中文字幕在线| 菩萨蛮人人尽说江南好唐韦庄| 女的被弄到高潮叫床怎么办| 欧美日韩在线观看h| 国产精品一二三区在线看| 午夜日本视频在线| 曰老女人黄片| 天堂中文最新版在线下载| 夜夜骑夜夜射夜夜干| 久久精品人人爽人人爽视色| 一区二区三区四区激情视频| 在线精品无人区一区二区三| 国国产精品蜜臀av免费| 男女无遮挡免费网站观看| 黑人高潮一二区| 久久久久久久亚洲中文字幕| 久久久久视频综合| 91成人精品电影| 亚洲综合色惰| 亚洲少妇的诱惑av| 九九久久精品国产亚洲av麻豆| 蜜桃国产av成人99| 亚洲精品日韩av片在线观看| 性色av一级| 国产黄频视频在线观看| 少妇人妻精品综合一区二区| 交换朋友夫妻互换小说| 一二三四中文在线观看免费高清| 美女中出高潮动态图| 高清av免费在线| 99热这里只有精品一区| 亚洲欧美一区二区三区黑人 | 制服诱惑二区| 国产日韩欧美视频二区| 国产精品 国内视频| 简卡轻食公司| 天天躁夜夜躁狠狠久久av| 下体分泌物呈黄色| 亚洲精品aⅴ在线观看| 国产老妇伦熟女老妇高清| 三级国产精品欧美在线观看| 中文字幕精品免费在线观看视频 | a级毛片在线看网站| 亚洲精品日本国产第一区| 人妻夜夜爽99麻豆av| 大香蕉97超碰在线| 国产有黄有色有爽视频| 夜夜骑夜夜射夜夜干| 亚洲欧美精品自产自拍| 亚洲av免费高清在线观看| 亚洲欧美精品自产自拍| 青青草视频在线视频观看| 在线精品无人区一区二区三| 人人澡人人妻人| 91精品三级在线观看| 久久 成人 亚洲| 国产白丝娇喘喷水9色精品| 精品卡一卡二卡四卡免费| 国产男女超爽视频在线观看| 99九九在线精品视频| 天堂中文最新版在线下载| av女优亚洲男人天堂| 日韩av不卡免费在线播放| 久久精品国产亚洲网站| 超色免费av| 国产精品一二三区在线看| 哪个播放器可以免费观看大片| 亚洲精品乱码久久久v下载方式| 91aial.com中文字幕在线观看| 青春草国产在线视频| 天天躁夜夜躁狠狠久久av| 欧美xxxx性猛交bbbb| 亚洲av中文av极速乱| 在线观看www视频免费| 午夜av观看不卡| 内地一区二区视频在线| 九九久久精品国产亚洲av麻豆| 国产一区有黄有色的免费视频| 啦啦啦在线观看免费高清www| 大又大粗又爽又黄少妇毛片口| 一级毛片电影观看| 色5月婷婷丁香| 天堂8中文在线网| 久久精品人人爽人人爽视色| 国语对白做爰xxxⅹ性视频网站| 满18在线观看网站| 熟女人妻精品中文字幕| 久久免费观看电影| 热re99久久精品国产66热6| 少妇被粗大猛烈的视频| 久久久久久久久久久丰满| 日韩大片免费观看网站| 日韩一区二区视频免费看| 中文字幕最新亚洲高清| 亚洲在久久综合| 国产免费一区二区三区四区乱码| 日本黄大片高清| 欧美日韩av久久| 国产精品偷伦视频观看了| 亚洲av不卡在线观看| 我的女老师完整版在线观看| 伊人久久国产一区二区| 午夜激情久久久久久久| 免费看不卡的av| 国产黄片视频在线免费观看| 国产高清不卡午夜福利| 女性被躁到高潮视频| 一级,二级,三级黄色视频| 免费看av在线观看网站| 国产精品秋霞免费鲁丝片| 国产日韩欧美在线精品| av天堂久久9| 成年人免费黄色播放视频| 王馨瑶露胸无遮挡在线观看| 欧美xxⅹ黑人| 99热6这里只有精品| 最黄视频免费看| 丝袜喷水一区| av一本久久久久| 免费大片黄手机在线观看| 日韩免费高清中文字幕av| 国产成人精品久久久久久| 亚洲欧美日韩另类电影网站| 欧美人与性动交α欧美精品济南到 | 欧美精品高潮呻吟av久久| a级毛片免费高清观看在线播放| 久久久久久久久久成人| 精品卡一卡二卡四卡免费| 边亲边吃奶的免费视频| 亚洲欧美日韩另类电影网站| 精品国产乱码久久久久久小说| 亚洲av成人精品一区久久| 亚洲av电影在线观看一区二区三区| 女人久久www免费人成看片| 老司机亚洲免费影院| 制服诱惑二区| 免费黄频网站在线观看国产| av一本久久久久| 看十八女毛片水多多多| 国产精品嫩草影院av在线观看| 久久久午夜欧美精品| 自线自在国产av| 午夜久久久在线观看| 国产在线视频一区二区| 久久国内精品自在自线图片|