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

    中國東部AOD等級變化及與東亞夏季風(fēng)的聯(lián)系

    2018-09-25 01:19:10馬奮華管兆勇
    中國環(huán)境科學(xué) 2018年9期
    關(guān)鍵詞:季風(fēng)距平短波

    馬奮華,管兆勇

    ?

    中國東部AOD等級變化及與東亞夏季風(fēng)的聯(lián)系

    馬奮華,管兆勇*

    (南京信息工程大學(xué)氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室,氣象災(zāi)害預(yù)報(bào)預(yù)警與評估協(xié)同創(chuàng)新中心,江蘇 南京 210044)

    利用2000年3月~2017年1月Terra衛(wèi)星反演的最新版本C06的MODIS氣溶膠光學(xué)厚度(AOD)、NCEP fnl全球業(yè)務(wù)分析數(shù)據(jù)、CMAP降水、CERES SYN1deg Ed4月平均資料、NCEP/NCAR再分析資料等,依據(jù)區(qū)域平均AOD距平的顯著年際變化特征,將歷年AOD劃分為5個(gè)不同的污染等級,并探討了不同污染等級的局地成因及其與季風(fēng)環(huán)流的聯(lián)系.結(jié)果表明,夏季中國東部地區(qū)AOD異常偏大與地面風(fēng)速小、中高層季風(fēng)環(huán)流系統(tǒng)不利于氣溶膠擴(kuò)散有關(guān).氣溶膠污染最重(最輕)時(shí),氣溶膠與到達(dá)地面的太陽短波輻射、地面氣溫、風(fēng)速及降水存在密切聯(lián)系;在其他污染等級下,除了與到達(dá)地面的短波輻射有密切關(guān)系外,與其他量的聯(lián)系并不明顯.氣溶膠污染嚴(yán)重時(shí),到達(dá)地面的太陽短波輻射相對減小,地面氣溫異常偏低,低層大氣冷卻,地面風(fēng)速減小,地面降水呈現(xiàn)南多北少的變化特征;反之,氣溶膠污染較輕時(shí),到達(dá)地面的太陽短波輻射相對增加,地面氣溫異常偏高,低層大氣加熱,地面風(fēng)速增大,地面降水呈現(xiàn)南少北多的變化特征.

    中國東部地區(qū);AOD;東亞夏季風(fēng)

    近幾十年來,中國東部氣溶膠污染日趨嚴(yán)重[1-4],城市地區(qū)氣溶膠排放濃度最高達(dá)偏遠(yuǎn)地區(qū)排放的2.5倍[5],東部站點(diǎn)AOD是西部站AOD的2倍[6]. PM2.5與PM10[7]、黑碳[8-9]、礦物沙塵[10]等污染物的濃度亦均呈現(xiàn)較高的態(tài)勢,嚴(yán)重影響了大氣環(huán)境.同時(shí),中國東部受季風(fēng)的影響,季風(fēng)活動異??稍斐蓢?yán)重的氣象災(zāi)害,影響人們的生產(chǎn)活動和日常生活[11].因此,氣溶膠與季風(fēng)被認(rèn)為是影響中國東部氣候環(huán)境的兩大重要影響因子[12-14].

    研究表明,氣溶膠可影響季風(fēng)變動.氣溶膠通過對輻射強(qiáng)迫產(chǎn)生影響[15-18],致使海陸熱力差異改變,從而影響季風(fēng).20世紀(jì)70年代以來,東亞夏季風(fēng)減弱,長江中下游地區(qū)降水偏多、華北地區(qū)降水偏少,出現(xiàn)南澇北旱的年代際變化現(xiàn)象[19-23].吸收性氣溶膠增加可能是導(dǎo)致中國中東部夏季氣候產(chǎn)生突變的原因[24-25],而散射性氣溶膠(硫酸鹽氣溶膠)的直接氣候效應(yīng)顯著地減少了到達(dá)地面的太陽輻射,導(dǎo)致海陸溫差縮小,使東亞夏季風(fēng)強(qiáng)度減弱[26],也可能是各種氣溶膠綜合作用的結(jié)果[27].數(shù)值模擬結(jié)果也表明,如果模式同時(shí)考慮散射性氣溶膠與吸收性氣溶膠,結(jié)果也顯示出氣溶膠增加可使東亞夏季風(fēng)減弱[28-29].

    氣溶膠時(shí)空變化也受季風(fēng)影響.模式研究結(jié)果表明,夏季風(fēng)能夠影響氣溶膠的濃度及空間分布[30-35].模式模擬結(jié)果表明,夏季風(fēng)是氣溶膠濃度變化的主要影響因子[30],而不是氣溶膠的局地排放[31]. PM2.5的濃度在弱夏季風(fēng)年要高于強(qiáng)夏季風(fēng)年[32],對南亞夏季風(fēng)而言,也有類似的結(jié)論[31].由季風(fēng)環(huán)流變化引起的風(fēng)場變化是氣溶膠遠(yuǎn)距離輸送的主要?dú)庀笠蜃覽35-36],而季風(fēng)帶來的濕空氣也是氣溶膠濃度下降的主要原因[34],夏季風(fēng)的季節(jié)進(jìn)程對氣溶膠的再分布也產(chǎn)生影響[37].

    盡管上述有關(guān)氣溶膠與季風(fēng)的關(guān)系有了大量研究,但不同條件下氣溶膠污染的程度是不同的.如果將污染分成若干等級,這些不同等級與季風(fēng)環(huán)流異常是什么關(guān)系,目前尚不是非常清楚.由于氣溶膠光學(xué)厚度表征了大氣柱中氣溶膠對光的消散作用,其量值在一定程度上反映出大氣中氣溶膠濃度的大小,且有研究表明,AOD的變化特征基本上與氣溶膠排放源的年(代)際變化相一致[38-41].因此,AOD被用來分析氣溶膠與季風(fēng)變化的聯(lián)系是一較好的選擇.本文將對AOD進(jìn)行分級,并分析其與季風(fēng)環(huán)流異常的聯(lián)系,旨在認(rèn)識氣溶膠與季風(fēng)相互作用.

    1 資料與方法

    為了研究中國東部夏季AOD的分布特征及其與季風(fēng)環(huán)流異常的聯(lián)系,本文選取2000年3月~2017年1月的數(shù)據(jù)資料.

    AOD資料:搭載于Terra衛(wèi)星上的中分辨率成像光譜儀(MODIS)是美國地球觀測系統(tǒng)(EOS)計(jì)劃中用于觀測全球生物和物理過程的重要儀器.它具有36個(gè)中等分辨率水平(0.25um~1um)的光譜波段,每1~2d對地球表面觀測一次.本文選取MODIS最新版本C06的大氣三級產(chǎn)品中的月平均資料,空間分辨率為1°×1°. C06資料較C05資料在地面反照率估計(jì)、氣溶膠模式選擇及云篩查方案等方面均有所改進(jìn)[42-43],尤其是植被與非植被地表混合的地表反照率估計(jì)得到改進(jìn),使氣溶膠數(shù)據(jù)集的覆蓋范圍從干旱半干旱地區(qū)擴(kuò)展到整個(gè)大陸地區(qū),且與地表觀測網(wǎng)數(shù)據(jù)集(Aerosol Robotic Network,簡稱AERONET)比較,氣溶膠的反演誤差減小.本文基于其中的550nm波長范圍的深藍(lán)與暗目標(biāo)算法相結(jié)合反演的AOD進(jìn)行分析.

    由于上午星MODIS/Terra與下午星MODIS/ Aqua過境時(shí)間不同,所測得的AOD可能不同.本文分別對MODIS/Terra與MODIS/Aqua夏季區(qū)域平均AOD的歷年變化分析時(shí)發(fā)現(xiàn),兩星的年際變化一致,且兩者差異不大.由于Terra衛(wèi)星觀測時(shí)間序列相對長了2a,因此本文分析時(shí)選用了Terra衛(wèi)星的數(shù)據(jù)資料.

    美國國家環(huán)境預(yù)報(bào)中心(NCEP)提供的業(yè)務(wù)全球分析資料:空間分辨率為1°×1°,為6h一次的逐時(shí)資料,垂直方向?yàn)?1層,分析時(shí)將逐時(shí)資料處理為月平均數(shù)據(jù)資料.本文主要利用地面10m處的風(fēng)速、地面2m處氣溫,等壓面風(fēng)場等物理量進(jìn)行分析.

    云地球輻射能量系統(tǒng)(Clouds and the Earth’s Radiant Energy System,簡稱CERES)CERES SYN1deg Ed4月平均資料,空間分辨率為1°×1°.

    CMAP降水資料[44],空間分辨率為2.5°×2.5°.

    NCEP/NCAR再分析月平均資料[45],空間分辨率為2.5°×2.5°,垂直方向?yàn)?7層,主要用于計(jì)算東亞夏季風(fēng)指數(shù).

    本文所用主要方法包括距平分析、相關(guān)分析及檢驗(yàn)等.夏季定義為6~8月,而夏季平均指6~8月平均.

    2 結(jié)果與討論

    2.1 夏季不同污染等級劃分與AOD分布特征

    2.1.1 AOD變化與風(fēng)速的關(guān)系 夏季中國東部處于污染物濃度高值區(qū).從近17a夏季平均AOD分布可看出,中國東部包括山東、京津冀以南、山西、江蘇、安徽、河南、湖北、浙江以及江西與湖南北部,AOD值均在0.6以上,AOD呈西南-東北向帶狀分布,且一直延伸到山東沿海.此時(shí),地面上10m處風(fēng)速在1.5m/s以下,比周圍及東南沿海風(fēng)速均偏小.對流層低層850hPa風(fēng)場由冬季盛行的偏北風(fēng)轉(zhuǎn)變?yōu)橄募臼⑿械哪巷L(fēng),越赤道西南氣流從孟加拉灣流經(jīng)南海時(shí)風(fēng)速減小,同時(shí)由西南氣流轉(zhuǎn)變?yōu)槟巷L(fēng)氣流,從華南經(jīng)江淮流域到達(dá)華北時(shí)風(fēng)速逐漸減小,尤其是在AOD大值區(qū),風(fēng)速達(dá)到最小,在2m/s以下.同時(shí)低層印度西南季風(fēng)攜帶水汽經(jīng)南海輸送至華南,水汽通量減弱,在黃河流域以北水汽通量減小到2g/s,達(dá)到最小.

    圖1 2000~2016年多年平均的夏季各變量分布

    方框表示中國東部地區(qū)

    較大的AOD可能使到達(dá)地表的太陽輻射減小,并減小海陸溫差,因而不利于夏季風(fēng)增強(qiáng)(如孫家仁等[26,28]).而弱的地面風(fēng)速亦不利于污染物的擴(kuò)散,導(dǎo)致AOD增大(如Yan等[32]).進(jìn)一步分析AOD與地面及850hPa風(fēng)速的相關(guān)系數(shù)分布(圖2),發(fā)現(xiàn)中國東部地區(qū)AOD與地面及850hPa風(fēng)速大小呈負(fù)相關(guān)關(guān)系.尤其是AOD與地面風(fēng)速的負(fù)相關(guān)關(guān)系(圖2a)在江蘇大部分地區(qū)及山東西南與河南接壤處很強(qiáng),相關(guān)系數(shù)達(dá)到-0.60(通過0.05置信水平檢驗(yàn)).說明AOD異常偏高,地面風(fēng)速小;AOD異常偏低,地面風(fēng)速大;反之亦然.這從觀測資料分析的角度進(jìn)一步說明風(fēng)速越小,越不利于氣溶膠等污染物擴(kuò)散.

    圖2 2000~2016年夏季平均AOD異常與地面及850風(fēng)速異常的相關(guān)系數(shù)分布

    方框表示中國東部地區(qū),黑點(diǎn)通過0.05置信水平檢驗(yàn)

    上述分析表明,AOD變化與地面和對流層低層風(fēng)速變化關(guān)系密切.但在不同AOD大小的變化情況下風(fēng)速及大氣環(huán)流則不能很好地刻畫.為此,首先依據(jù)AOD數(shù)值變化劃分污染程度等級.

    2.1.2 污染程度的等級劃分 對中國東部地區(qū)[110~121°E,27.5~40°N]求AOD歷年夏季區(qū)域平均值(記為[AOD]),并給出歷年[AOD]的距平變化(圖3).從歷年AOD變化特征看,[AOD]距平分布總體呈現(xiàn)先增長后下降的趨勢,但年際變化也很明顯.有些年份如2007、2008、2010、2011、2014年,[AOD]值比平均值高出0.1以上,甚至個(gè)別年份高出平均值0.2以上;有些年份如2000、2013、2016年,[AOD]值比平均值低0.1,尤其是2016年夏季比平均值低0.2以上.這說明近17a,中國東部區(qū)域平均AOD值的年際差異比較大.

    依據(jù)[AOD]距平值的不同大小范圍,將近17a氣溶膠的污染程度劃分為5級(表1),即污染嚴(yán)重、污染重、污染較重、污染輕、污染較輕,相應(yīng)地定義污染等級為Level-Ⅴ、Level-Ⅳ、Level-Ⅲ、Level-Ⅱ、Level-Ⅰ,并對應(yīng)相應(yīng)的年份.本文從這5個(gè)等級出發(fā),研究季風(fēng)環(huán)流對氣溶膠的影響.

    按照5個(gè)不同等級的年份,合成分析了不同污染等級下中國東部夏季AOD距平分布(圖4).由圖4可以明顯看出,隨著AOD從Levlel-Ⅰ變化到Level-Ⅴ,中國東部區(qū)域尤其是北京、天津、河北南部以及山東半島以西AOD距平從負(fù)值逐漸變?yōu)檎?且距平值越來越大,到Level-Ⅴ級時(shí),整個(gè)中國東部及周邊地區(qū)均為正距平區(qū),AOD北高南低,與Wang等[45]分析結(jié)果一致.Level-Ⅴ與Level-Ⅰ相比,AOD距平差達(dá)到0.6(均通過了0.10置信水平距平差值檢驗(yàn)),說明中國東部AOD在不同分級之間不僅區(qū)域平均差異大,而且區(qū)域分布特征也明顯不同.

    圖3 2000~2016年中國東部夏季區(qū)域平均值距平的逐年變化

    直線為線性趨勢,曲線為二次擬合分布,2稱為決定系數(shù),表示二次曲線的擬合程度高低

    表1 中國東部氣溶膠污染分級

    圖4 不同污染等級下中國東部夏季AOD距平的分布

    方框表示中國東部地區(qū),黑點(diǎn)表示通過0.10置信水平檢驗(yàn)

    2.2 不同等級污染的局地成因及可能影響

    中國東部污染不同等級與地面10m處風(fēng)速距平大小有關(guān).由圖5可見,污染較輕時(shí)(Level-Ⅰ),整個(gè)中國東部地面風(fēng)速距平值為正,最大正距平達(dá)到0.3m/s.隨著污染的加重(Level-Ⅱ到Level-Ⅳ),風(fēng)速距平分布變化明顯,但變化規(guī)律不明顯.到污染最重時(shí)(Level-Ⅴ),除津京冀一帶整個(gè)中國東部風(fēng)速距平變化為負(fù)值,尤其是東部沿海地區(qū)風(fēng)速距平變化最大,與歷年平均風(fēng)速相比減小最明顯.對比Level-Ⅴ與Level-Ⅰ兩級地面風(fēng)速距平分布可以說明,中國東部地區(qū)AOD異常偏高時(shí),對應(yīng)地面風(fēng)速小,AOD異常偏低時(shí),對應(yīng)地面風(fēng)速大.反之也說明,地面風(fēng)速大,AOD小;地面風(fēng)速小,AOD大.地面風(fēng)速大小,是AOD異常分布的影響因子之一,風(fēng)速大,氣溶膠容易擴(kuò)散;風(fēng)速小,氣溶膠容易積聚.但在津京冀一帶,污染較輕與污染較重的情況下,地面風(fēng)速距平大小與其他區(qū)域表現(xiàn)出相反的變化特征,即污染較輕時(shí),地面風(fēng)速距平值為負(fù),污染較重時(shí),地面風(fēng)速距平值為正,這說明尚有其他原因?qū)е翧OD的異常.

    圖5 不同污染等級下中國東部夏季地面10m處地面風(fēng)速距平的分布

    方框表示中國東部地區(qū),黑點(diǎn)表示通過0.10置信水平檢驗(yàn)

    圖6 不同污染等級下中國東部夏季晴空地面向下的短波輻射距平的分布

    方框表示中國東部地區(qū),黑點(diǎn)表示通過0.10置信水平檢驗(yàn)

    不考慮云的影響,氣溶膠污染明顯使得到達(dá)地面的太陽短波輻射通量減小.由圖6可見,污染較輕時(shí)(Level-Ⅰ),到達(dá)地面的太陽短波輻射通量距平值為正,正距平為8W/m2以上的區(qū)域幾乎覆蓋了整個(gè)中國東部地區(qū),與AOD負(fù)距平分布區(qū)域范圍幾乎一致(圖4,Level-Ⅰ).隨著污染的逐步加劇(Level-Ⅱ~Level-Ⅳ),到達(dá)地面的太陽短波輻射通量距平分布發(fā)生明顯變化,這種變化特征基本上與AOD距平分布呈一致的反相變化關(guān)系,即在AOD負(fù)距平區(qū),短波輻射通量為正距平,AOD正距平區(qū),短波輻射通量為負(fù)距平區(qū).達(dá)到污染最嚴(yán)重時(shí)(Level-Ⅴ),整個(gè)中國東部短波輻射通量為負(fù)距平區(qū),污染最重的北京、天津、河北南部及河南北部,短波輻射通量負(fù)距平值達(dá)最大,為-6W/m2.這種AOD與輻射異常的反相分布的相關(guān)關(guān)系很好地反映了氣溶膠對到達(dá)地面的太陽短波輻射通量的影響.但反過來,因短波輻射的減少而導(dǎo)致大氣層結(jié)穩(wěn)定度的改變進(jìn)而影響AOD變化的過程仍需深入探討.

    不同污染等級下的輻射變化可影響中國東部夏季地面2m氣溫的距平分布.由圖7可見,氣溶膠污染較輕時(shí)(Level-Ⅰ),地面氣溫距平值為正.污染逐步加劇時(shí)(Level-Ⅱ到Level-Ⅳ),地面氣溫距平分布發(fā)生變化,但與AOD距平分布(圖4,Level-Ⅱ到Level-Ⅳ)變化對應(yīng)關(guān)系較差.當(dāng)氣溶膠污染嚴(yán)重時(shí)(Level-Ⅴ),中國東部地區(qū)除京津冀南邊小范圍為正距平外,其他地區(qū)地面氣溫距平均為負(fù)值.對比Level-Ⅰ與Level-Ⅴ等級下AOD與地面氣溫的距平分布特征發(fā)現(xiàn),AOD異常偏高時(shí),地面氣溫相對下降;AOD異常偏低時(shí),地面氣溫相對增高.這說明因氣溶膠污染,使到達(dá)地面的太陽短波輻射減小,造成地面氣溫降低.反過來也說明,當(dāng)?shù)孛鏆鉁鼐嗥街递^大,地面氣溫較高時(shí),地面邊界層內(nèi)對流容易發(fā)生,氣溶膠較易向高層擴(kuò)散,可能使區(qū)域氣溶膠濃度降低;當(dāng)?shù)孛鏆鉁鼐嗥街递^小,地面氣溫較低時(shí),地面邊界層內(nèi)對流不易發(fā)生,氣溶膠不易向高層擴(kuò)散,可能使區(qū)域氣溶濃度增加.同時(shí)從上述分析發(fā)現(xiàn),氣溶膠對地面氣溫的影響,取決于氣溶膠污染的程度.

    圖7 不同污染等級下中國東部夏季地面2m處氣溫距平的分布

    方框表示中國東部地區(qū),黑點(diǎn)表示通過0.10置信水平檢驗(yàn)

    圖8 不同污染等級下中國東部夏季地面降水距平的分布

    方框表示中國東部地區(qū),黑點(diǎn)表示通過0.10置信水平檢驗(yàn)

    不同污染等級與地面降水距平分布亦存在一定聯(lián)系.由圖8可見,在污染較輕時(shí)(Level-Ⅰ),中國東部江淮流域?yàn)榻邓?fù)距平區(qū),江淮流域以北至京津冀陜一帶為降水正距平區(qū),降水呈現(xiàn)“南少北多”的分布特征.而污染嚴(yán)重時(shí)(Level-Ⅴ),中國東部長江流域?yàn)榻邓嗥絽^(qū),淮河流域及以北地區(qū)為降水負(fù)距平區(qū),降水呈現(xiàn)南多北少的分布特征,這種分布特征與Menon等[25]數(shù)值試驗(yàn)結(jié)果一致,說明異常氣溶膠變化與降水異常存在著可能的聯(lián)系.但在其他污染等級下(Level-Ⅱ~Level-Ⅳ),AOD距平分布與降水距平分布關(guān)系不是很明顯.

    綜合以上分析可見,AOD異常偏高時(shí)(北高南低),氣溶膠污染使得到達(dá)地面的太陽短波輻射減小,地面氣溫下降,大氣冷卻,地面風(fēng)速減小,地面降水呈現(xiàn)南多北少的變化特征;AOD異常偏低時(shí)(北低南高),到達(dá)地面的太陽短波輻射增加,地面氣溫相對增高,大氣加熱,地面風(fēng)速增加,地面降水呈現(xiàn)南少北多的變化特征.這些結(jié)論從年際變化和觀測證據(jù)角度,進(jìn)一步說明氣溶膠直接效應(yīng)造成中國東部氣候變化的事實(shí).

    2.3 氣溶膠不同污染等級與季風(fēng)環(huán)流的聯(lián)系

    2.3.1 夏季風(fēng)強(qiáng)度與氣溶膠的聯(lián)系 為探討東亞夏季風(fēng)強(qiáng)度與氣溶膠的聯(lián)系,計(jì)算了2000~2016年夏季反映東西海陸熱力差對照指數(shù)(施能等[47];趙平等[48]),緯向風(fēng)距平差指數(shù)(張慶云等[49];Wang等[50]),經(jīng)向風(fēng)距平指數(shù)(Wu等[51];Wang等[52]),既考慮海陸熱力對照又考慮緯向風(fēng)切變的指數(shù)(祝叢文等[53])以及風(fēng)速大小指數(shù)(Wang等[54]),并求歷年各指數(shù)標(biāo)準(zhǔn)化距平與歷年中國東部區(qū)域平均[AOD]標(biāo)準(zhǔn)化距平的相關(guān)系數(shù)(表2).由表2可見,這8種夏季風(fēng)指數(shù)與AOD的相關(guān),有4種為負(fù)相關(guān)關(guān)系,分別為SNZQYWB及WYF指數(shù),有4種為正相關(guān)關(guān)系,分別為ZPWAMZCW及WHJ指數(shù),且各指數(shù)與AOD的相關(guān)系性都不是很大,均沒有通過0.10置性水平檢驗(yàn).說明各指數(shù)在表征東亞夏季風(fēng)強(qiáng)度方面雖然都有各自的優(yōu)勢,且能分析判斷東亞夏季風(fēng)的強(qiáng)弱,但在判斷大氣環(huán)流場對AOD的影響方面卻不是最好的表征量.這也從另一方面說明,大氣環(huán)流對AOD的影響是多方面綜合作用的結(jié)果,影響比較復(fù)雜.

    表2 各類夏季風(fēng)指數(shù)與AOD的相關(guān)系數(shù)

    2.3.2 季風(fēng)環(huán)流與氣溶膠AOD變化的聯(lián)系 不同污染等級下850hPa與200hPa旋轉(zhuǎn)風(fēng)與輻散風(fēng)的距平分布呈現(xiàn)較為復(fù)雜的特征.由圖9可見,污染較輕時(shí)(Level-Ⅰ),中國東部低層850hPa受異常反氣旋式環(huán)流西側(cè)的異常東南氣流所控制,高層200hPa上受平直的異常西風(fēng)氣流控制,這樣低層利于暖濕氣流而高層利于干冷氣流存在,利于造成對流層中上層對流活動;污染嚴(yán)重時(shí)(Level-Ⅴ),對流層低層受氣旋式異常環(huán)流控制,高層受異常西南-南風(fēng)氣流控制,再加上地面風(fēng)速減小,氣溫降低,抑制了季風(fēng)區(qū)對流的發(fā)生,利于對流層低層污染的積聚.

    注意到在Level-Ⅴ等級和Level-Ⅰ等級上環(huán)流異常的分布型并未明顯相反,這說明氣溶膠污染與季風(fēng)環(huán)流的關(guān)系并不是一種線性的變化關(guān)系,這也是為什么表2中各季風(fēng)指數(shù)與AOD相關(guān)系數(shù)較小的原因.

    黑色(細(xì)灰)代表U,V風(fēng)場均沒通過(均通過)0.1置信水平檢驗(yàn),粗灰代表U,V其中一個(gè)量通過0.1置信水平檢驗(yàn)

    2.4 討論

    由于重污染與輕污染發(fā)生時(shí)季風(fēng)環(huán)流異常型并不完全相反,說明AOD變化與季風(fēng)環(huán)流變化并無顯著的線性關(guān)系,這也說明了雖然東亞夏季風(fēng)指數(shù)在表征東亞夏季風(fēng)強(qiáng)度方面都有各自的優(yōu)勢,但在判斷環(huán)流異常對AOD影響方面卻不是最好的表征量.如何提取AOD年際變化和季風(fēng)環(huán)流變化密切聯(lián)系的相關(guān)信息仍是一個(gè)有待解決的問題.

    還需要說明,由于MODIS氣溶膠光學(xué)厚度資料時(shí)間序列長度不足,且在AOD反演過程中可能存在一定誤差,文中結(jié)果尚需用模式或更長的時(shí)間序列進(jìn)一步驗(yàn)證.

    3 結(jié)論

    3.1 通過對AOD距平年際變化不同等級的劃分,將中國東部污染劃分為5個(gè)由輕至重的不同等級,并分析了AOD與環(huán)流異常的聯(lián)系.當(dāng)中國東部AOD值較高時(shí),地面風(fēng)速比周邊地區(qū)小,對流層低層流場比周邊地區(qū)弱,這種負(fù)相關(guān)關(guān)系的變化特征在AOD與地面風(fēng)速的相關(guān)分布上表現(xiàn)最明顯.

    3.2 氣溶膠污染最重與最輕時(shí)(Level-Ⅴ與Level-Ⅰ級),氣溶膠與到達(dá)地面的太陽短波輻射、地面氣溫、風(fēng)速及降水存在密切的聯(lián)系,其他污染等級下(Level-Ⅱ到Level-Ⅳ級),除了與到達(dá)地面的短波輻射有聯(lián)系外,與其他變量的關(guān)系并不十分密切.

    3.3 AOD異常偏高時(shí)(北高南低),到達(dá)地面的太陽短波輻射相對偏小,地面氣溫相對偏低,地面風(fēng)速減小,地面降水呈現(xiàn)南多北少的變化特征;相反,AOD異常偏低時(shí)(北低南高),到達(dá)地面的太陽短波輻射相對偏大,地面氣溫相對偏高,地面風(fēng)速加大,地面降水呈現(xiàn)“南少北多”的變化特征.

    [1] 羅云峰,呂達(dá)仁,李維亮.近30年來中國地區(qū)大氣氣溶膠光學(xué)厚度的變化特征[J]. 科學(xué)通報(bào), 2000,45(5):549–554.

    [2] 李成才,毛節(jié)泰,劉啟漢.利用MODIS研究中國東部地區(qū)氣溶膠光學(xué)厚度的分布和季節(jié)變化特征[J]. 科學(xué)通報(bào), 2003,48(19):2094– 2100.

    [3] 毛節(jié)泰,李成才.氣溶膠輻射特性的觀測研究[J]. 氣象學(xué)報(bào), 2005,63(5):622–635.

    [4] 段 婧,毛節(jié)泰.長江三角洲大氣氣溶膠光學(xué)厚度分布和變化趨勢研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 2007,27(4):537–543.

    [5] Zhang X Y, Wang Y Q, Niu T, et al. Atmospheric aerosol compositions in China: Spatial/temporal variability, chemical signature, regional haze distribution and comparisons with global aerosols [J]. Atmospheric Chemistry & Physics, 2012,11:26571–26615.

    [6] 張志薇,王宏斌,張 鐳,等.中國3個(gè)AERONET站點(diǎn)氣溶膠微物理特性分析及比較[J]. 中國環(huán)境科學(xué), 2014,34(8):1927–1937.

    [7] Oanh N T K, et al. Particulate air pollution in six Asian cities: Spatial and temporal distributions, and associated sources [J]. Atmospheric Environment, 2006,40:3367–3380.

    [8] Novakov T, Kirchstetter T W, Sinton J E, et al. A. Large historical changes of fossil– fuel black carbon aerosols [J]. Geophysical Research Letters, 30:1324, doi:10.1029/2002GL016345,2003.

    [9] 王志立,張 華,郭品文.南亞地區(qū)黑碳?xì)馊苣z對亞洲夏季風(fēng)的影響[J]. 高原氣象, 2009,28(2):419–424.

    [10] IPCC. Fifth Assessment Report, Climate Change 2013: The Physical Science Basis [R]. New York: Cambridge University Press, 2013.

    [11] Wang H J, Sun J Q, Chen H P, et al. Extreme climate in China: Facts, simulation and projection [J]. Meteorologische Zeitschrift, 2012, 21(3):279–304,doi:10.1127/0941–2948/2012/0330.

    [12] au K M W, Ramanathan V, Wu G X et al. The joint aerosol–monsoon experiment [J]. Bulletin of the American Meteorological society, 2008,89(3):369–383, DOI:10.1175/BAMS–89–3–369.

    [13] Lau K M W. The aerosol–monsoon climate system of Asia: A new paradigm [J]. Journal of Meteorological Research., 2016,30(1):1–11, doi:10.1007/s13351–015–5999–1.

    [14] Li Z Q, Lau K M W, Ramanathan V, et al. Aerosol and monsoon climate interactions over Asia [J]. Reviews of Geophysics, 2016,297: 2250–2253,doi:/10.1002/2015RG000500.

    [15] Ramanathan V, Coauthors. Indian Ocean Experiment: An integrated assessment of the climate forcing and effects of the great Indo_Asian haze [J]. Journal of Geophysical Research, 2001a,106,D22,28371– 28398.

    [16] Conant W C, Seinfeld J H, Wang J, et al. A model for the radiative forcing during ACE–Asia derived from CIRPAS Twin Otter and R/V Ronald H.Brown data and comparision with observation [J]. Journal of Geophysical Research, 108(D23),8661, DOI: 10.1029/2002JD003260, 2003.

    [17] Ramanathan V, Crutzen P J. Newdirection: Atmospheric brown clouds [J]. Atmospheric Environment, 2003,37:4033–4035.

    [18] Chung C E, Ramanthan V, Kim D, et al. Global anthropogenic aerosol direct forcing derived from satellite and ground–based observations [J]. Journal of Geophysical Research, 2005,110,D24207,doi:10.1029/ 2005JD006356.

    [19] Yang F, Lau K M. Trend and variability of China precipitation in spring and summer:Linkage to sea–surface temperature [J]. International Journal of Climatology, 2004,24:1625–1644.

    [20] Ding Y H, Wang Z, Sun Y. Interdecadal variation of the summer precipitation in East China and its association with decreasing Asian summer monsoon. Part I: observed evidences [J]. International Journal of Climatology, 2008,28:1139–1161.

    [21] Ding Y H, Liu Y, Sun Y, et al. Weakening of the Asina summer monsoon and its impact on the precipitation pattern in China [J]. Water Resources Research, 2010,26:423–439.

    [22] Jin D C, Guan Z Y. Summer Rainfall Seesaw between Hetao and the Middle and Lower Reaches of the Yangtze River and Its Relationship with the North Atlantic Oscillation [J]. Journal of Climate, 2017, 30:6629–6243,DOI:10.1175/JCLI–D–16–0760.1.

    [23] Jin D C, Hameed S N, Huo L W. Changes in ENSO Teleconnection over the Western Pacific Impacts the Eastern China Precipitation Dipole [J]. Journal of Climate, 2016,29:7587–7598,DOI:10.1175/ JCLI–D–16–0235.1.

    [24] Xu Q. Abrupt change of the middle–summer climate in central east China by the influence of atmospheric pollution [J]. Atmospheric Environment, 2001,35(30):5029–5040.

    [25] Menon S, Hansen J, Nazarenko L, et al. Climate effects of black carbon aerosols in China and India [J]. Science, 2002,297:2250–2253, doi:/10.1126/science.1075159.

    [26] 孫家仁,劉 煜:中國區(qū)域氣溶膠對東亞夏季風(fēng)的可能影響(I):硫酸鹽氣溶膠的影響[J]. 氣候變化研究進(jìn)展, 2008,2:111–116.

    [27] Gu Y, Liou K N, Xue Y, et al. Climatic effects of different aerosol types in China simulated by the University of California, Los Angeles atmospheric general circulation model [J]. Journal of Geophysical Research, 2006,111,D15201,doi:10.1029/2005JD006312.

    [28] 孫家仁,劉 煜.中國區(qū)域氣溶膠對東亞夏季風(fēng)的可能影響(Ⅱ):黑碳?xì)馊苣z及其與硫酸鹽氣溶膠的綜合影響[J]. 氣候變化研究進(jìn)展, 2008b,4(3),161-166.

    [29] Zhang H, Wang Z L, Wang Z Z. Simulation of direct radiative forcing of aerosols and their effects on East Asian climate using an interactive AGCM–aerosol coupled system [J]. Climate Dynamics, 2012,38: 1675–1693,DOI:10.1007/s00382–001–1131–0.

    [30] Zhang L, Liao H, Li J P. Impacts of Asian summer monsoon on seasonal and interannual variations of aerosols over eastern China [J]. Journal of Geophysical Research, 2010a,115,D00K05,doi:10.1029/ 2009JD012299.

    [31] Yan L, Liu X, Yang P, et al. Study of the impact of summer monsoon circulation on spatial distribution of aerosol in East Asia based on numerical simulations [J]. Journal of Applied Meteorology and Climatology, 2011,50(11):2270–2282,doi:10.1175/2011JAMC–D–11– 06.1.

    [32] Zhu J, Liao H, Li J. Increases in aerosol concentrations over eastern China due to the decadal–scale weakening of the East Asian summer monsoon [J]. Geophysical Research Letters, 2012,39,L09809,doi: 10.1029/2012GL051428.

    [33] Zhang L, Liao H, Li J P. Impact of the southeast Asian summer monsoon strength on the outflow of aerosols from South Asia [J]. Annals of Geophysics, 2010b,28:277–287.

    [34] Kim D, Chin M, Yu H B, et al. Sources, sinks, andtransatlantic transport of North Africandust aerosol: A multimodel analy sisand comparison with remotesensing data [J]. Journal of Geophysical Research: Atmosphere, 2014,119,6259–6277,doi:10.1002/2013JD02109 9.

    [35] Shahid M Z, Liao H, Li J P, et al. Seasonal variations of aerosols in Pakistan: Contributions of domestic anthropogenic emissions and transboundary transport [J]. Aerosol and Air Quality Research, 2015,15:1580–1600,doi:10.4209/aaqr.2014.12.0332.

    [36] 薛文博,付 飛,王金南,等.中國PM2.5跨區(qū)域傳輸特征數(shù)值模擬研究[J]. 中國環(huán)境科學(xué), 2014,34(6):1361–1368.

    [37] 孫 一,管兆勇,馬奮華,等.夏季東亞地區(qū)AOD與地面太陽輻射變化的聯(lián)系及季風(fēng)環(huán)流異常:季節(jié)趨勢影響[J]. 大氣科學(xué)學(xué)報(bào), 2015, 38(2):165–174.

    [38] Bond T C, Bhardwaj E, Dong R, et al. Historical emission of black and organic carbon aerosol from energy–related combusition, 1850~2000 [J]. Global Biogeochemical Cycles, 2017,21,GB2018,doi:10.1029/ 2006GB002840.

    [39] Smith S J, Aardenne J, Klimont. Anthropogenic sulfur dioxide emissions: 1850~2005 [J]. Atmospheric Chemistry & Physics, 2011, 11:1101–1116,doi:10.5194/acp–11–1101–2011.

    [40] Zhang H, Wang Z L, Guo P W, et al. A modeling study of the effects of direct radiative forcing due to carbonaceous aerosol on the climate in East Asia [J]. Advances in Atmospheric Sciences, 2009,26(1):57–66.

    [41] Diehl T, Heil A, Chin M, et al.Anthropogenic, biomass burning, and volcanic emissions of black carbon, organic carbon, and SO2from 1980 to 2010 for hindcast model experiments, Atmos [J]. Atmospheric Chemistry and Physics Discussions, 2012,12,24895–24954,doi: 10.5194/acpd–12–24895–2012.

    [42] Hsu N C, Jeong M–J, Bettenhausen C, et al. Enhanced Deep Blue aerosol retrieval algorithm: The second generation [J]. Journal of Geophysical Research: Atmosphere, 2013,118,9296–9315,doi: 10.1002/jgrd.50712.

    [43] Sayer A M, Munchak L A, Hsu N C, et al. MODIS Collection 6aerosol products: Comparison between Aqua’s e–Deep Blue, Dark Target, and “merged” data sets, and usage recommendations [J]. Journal of Geophysical Research, 2014,119,13:965–13,989,doi:10.1002/ 2014JD022453.

    [44] Xie P, Arkin P A. Global precipitation: A 17–year monthly analysis based on gauge observations, satellite estimates, and numerical model outputs [J]. Bulletin of the American Meteorological Society, 1997, 78:2539–2558.

    [45] Kalnay E, Kanamitsu M, Kistler R, et al. The NCEP/NCAR 40–year reanalysis project [J]. Bulletin of the American Meteorological Society, 1996,7:437–471.

    [46] Wang H L, Xie X N, Yan L B, et al. Impact of East Asian summer monsoon circulation on the regional aerosol distribution in observations and models [J]. Theoreticaland Applied Climatology, 2017,doi:10.1007/s00704-017-2199-1.

    [47] 施 能,朱乾根,吳彬貴.近40年東亞夏季風(fēng)及我國夏季大尺度天氣氣候異常[J]. 大氣科學(xué), 1996,20(5):575–583.

    [48] 趙 平,周自江.東亞副熱帶夏季風(fēng)指數(shù)及其與降水的關(guān)系[J]. 氣象學(xué)報(bào), 2005,63(6):933–941.

    [49] 張慶云,陶詩言,陳烈庭.東亞夏季風(fēng)指數(shù)的年際變化與東亞大氣環(huán)流 [J]. 氣象學(xué)報(bào), 2003,61(5):559–568.

    [50] Wang B, Fan Z. Choice of south Asian summer monsoon indices [J]. Bulletin of the American Meteorological Society, 1999,80:629–638.

    [51] Wu A M, Ni Y Q. The influence of Tibetan Plateau on the interannual variability of Asian monsoon [J]. Advances in Atmospheric Sciences, 1997,14(4):491–504.

    [52] Wang Y, Wang B, Oh J H. Impact of preceeding El Nino on the East Asian summer atmosphere circulation [J]. Journal of the Meteorological Society of Japan, 2001,79:575–589.

    [53] 祝從文,何金海,吳國雄.東亞季風(fēng)指數(shù)及其與大尺度熱力環(huán)流年際變化的關(guān)系[J]. 氣象學(xué)報(bào), 2000,58(4):391–402.

    [54] Wang H J. The instability of the East Asian summer monsoon–ENSO relations [J]. Advances in Atmospheric Sciences, 2002,19(1):1–11.

    致謝:所用MODIS氣溶膠光學(xué)厚度資料由MODIS大氣資料中心提供(https://modis-atmosphere.gsfc.nasa.gov/);NCEPfnl資料、CMAP降水資料及NCEP/NCAR資料由NOAA的地球系統(tǒng)研究實(shí)驗(yàn)室提供(https://www.esrl.noaa.gov/psd/data/gridded/tables/monthly.html);所用輻射資料為云地球輻射能量系統(tǒng)資料(https://ceres.larc.nasa.gov/ products- info.php?product=SYN1deg),利用Grads軟件軟件繪制了本文插圖.

    Features of graded AOD in East China in association with East Asian summer monsoon anomalies.

    MA Fen-hua, GUAN Zhao-yong*

    (Key Laboratory of Meteorological Disaster of Ministry of Education, Joint International Research Laboratory of Climate and Environment Change, Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science & Technology, Nanjing 210044, China)., 2018,38(9):3201~3210

    Using the data of the latest-released aerosol optical thickness (AOD) (Collection 6) as derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) on the Terra spacecraft, the Final Operational Global Analysis data (NCEP FNL) from National Centers for Environmental Prediction, Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP), CERES syn1deg Ed4monthly average data, NCEP/National Center for Atmospheric Research (NCAR) reanalysis data from March 2000 to February 2017, we have investigated the features of variations of graded AOD and their related circulation anomalies. As per the interannual variations of AOD averaged over the eastern part of China, the AOD is divided into 5 different levels, which corresponds to 5pollution grades. The results show that the anomalies of AOD in East China are larger when the surface wind speed is weaker together with the unfavorable aerosol diffusion condition of the anomalous monsoon circulation in the middle and high troposphere. When the aerosol pollution is strongest (weakest), variations of AOD is strongly related to the variations of the solar radiation at the earth surface, surface temperature, wind speed, and precipitation. At other AOD levels, the air pollution is significantly related to the surface solar radiation rather than other variables. When the air pollution is severe, the solar radiation at the earth surface is relatively lower, along with the lower surface temperature. The atmosphere in lower troposphere is cooling, and the wind speed is weakening. Meanwhile, the more precipitation occurs in the southern and less in the northern part of the eastern China. On the other hand, when the air pollution is very light, the solar radiation that reaches the earth surface will increase, and the surface air temperature gets higher. At this time, the surface wind will get stronger, and less precipitation will be received in the southern and more in the northern part of East China.

    East China;aerosol optical depth (AOD);East Asian summer monsoon

    X16

    A

    1000-6923(2018)09-3201-10

    馬奮華(1974-),女,山西忻州人,副編審,碩士,研究方向?yàn)閰^(qū)域氣候變化.發(fā)表論文3篇.

    2018-01-14

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)項(xiàng)目(2011CB403406);江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項(xiàng)目(PAPD)

    * 責(zé)任作者, 教授, guanzy@nuist.edu.cn

    猜你喜歡
    季風(fēng)距平短波
    颶風(fēng)Edouard(2014)暖心結(jié)構(gòu)的多資料對比分析
    綠水青山圖——海洋季風(fēng)的贊歌
    近40年阿里地區(qū)云量和氣溫的年際變化
    西藏科技(2018年9期)2018-10-17 05:51:30
    戶撒刀
    樂海短波
    人民音樂(2016年1期)2016-11-07 10:02:42
    甘肅省降水和冰雹天氣氣候分析
    工運(yùn)短波
    工運(yùn)短波
    綠野短波
    萬馬奔騰
    中華奇石(2015年9期)2015-07-09 18:32:58
    老熟妇乱子伦视频在线观看| 色精品久久人妻99蜜桃| 国产一区在线观看成人免费| 九九在线视频观看精品| 波多野结衣巨乳人妻| 成年女人毛片免费观看观看9| 免费看十八禁软件| 久久99热这里只有精品18| 在线国产一区二区在线| 一本久久中文字幕| 一a级毛片在线观看| 一个人看视频在线观看www免费 | 成人性生交大片免费视频hd| 日韩 欧美 亚洲 中文字幕| 伊人久久精品亚洲午夜| 亚洲国产精品久久男人天堂| 1024手机看黄色片| 亚洲熟妇熟女久久| 国产视频一区二区在线看| 色综合婷婷激情| 亚洲av不卡在线观看| 日韩欧美精品免费久久 | 亚洲精品在线美女| 国内精品一区二区在线观看| 亚洲国产中文字幕在线视频| 国产精品久久久久久人妻精品电影| 男女做爰动态图高潮gif福利片| 午夜两性在线视频| 亚洲av中文字字幕乱码综合| 国产美女午夜福利| 中文字幕av成人在线电影| 亚洲一区高清亚洲精品| 亚洲国产精品成人综合色| 中文字幕人妻丝袜一区二区| 两个人的视频大全免费| 国产一区二区三区在线臀色熟女| 亚洲自拍偷在线| 在线观看免费视频日本深夜| 免费一级毛片在线播放高清视频| 女警被强在线播放| 亚洲国产精品999在线| 中文亚洲av片在线观看爽| 精品国内亚洲2022精品成人| 97超级碰碰碰精品色视频在线观看| 久久久久久久精品吃奶| 人人妻人人澡欧美一区二区| 国产亚洲精品av在线| 丰满人妻熟妇乱又伦精品不卡| 一级作爱视频免费观看| 中文字幕熟女人妻在线| 亚洲成人精品中文字幕电影| 亚洲va日本ⅴa欧美va伊人久久| 成熟少妇高潮喷水视频| 亚洲国产精品成人综合色| 国产黄a三级三级三级人| 美女高潮喷水抽搐中文字幕| 精品国产超薄肉色丝袜足j| xxx96com| xxxwww97欧美| 一边摸一边抽搐一进一小说| 99在线人妻在线中文字幕| 18禁黄网站禁片免费观看直播| 一个人看视频在线观看www免费 | 亚洲国产精品合色在线| 一本精品99久久精品77| 国产中年淑女户外野战色| 一本久久中文字幕| 长腿黑丝高跟| 好男人电影高清在线观看| 亚洲国产日韩欧美精品在线观看 | 午夜免费成人在线视频| 精品不卡国产一区二区三区| 亚洲成av人片在线播放无| 亚洲人成网站在线播| 十八禁人妻一区二区| 叶爱在线成人免费视频播放| 国产精品一及| 午夜福利18| 日韩欧美在线乱码| 国产97色在线日韩免费| 欧美成人a在线观看| 国产一区二区亚洲精品在线观看| 国产精品香港三级国产av潘金莲| 一区二区三区激情视频| 日日摸夜夜添夜夜添小说| 欧美3d第一页| 欧美成狂野欧美在线观看| 天美传媒精品一区二区| 长腿黑丝高跟| 亚洲美女视频黄频| 久久久久久人人人人人| 亚洲国产日韩欧美精品在线观看 | 亚洲一区二区三区色噜噜| 亚洲av成人精品一区久久| 非洲黑人性xxxx精品又粗又长| 香蕉丝袜av| 国产精品久久久人人做人人爽| 亚洲精品在线观看二区| 网址你懂的国产日韩在线| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲精品一区二区www| 午夜精品一区二区三区免费看| 热99re8久久精品国产| 女人十人毛片免费观看3o分钟| 亚洲内射少妇av| 在线观看一区二区三区| 欧美又色又爽又黄视频| 国产极品精品免费视频能看的| 免费看光身美女| 亚洲人成网站在线播| 国产淫片久久久久久久久 | 女同久久另类99精品国产91| 国内精品一区二区在线观看| 中文字幕高清在线视频| 一个人免费在线观看电影| 女生性感内裤真人,穿戴方法视频| 欧美午夜高清在线| 男人舔女人下体高潮全视频| 成年女人看的毛片在线观看| 看黄色毛片网站| 法律面前人人平等表现在哪些方面| 热99在线观看视频| 手机成人av网站| 久久精品国产亚洲av涩爱 | 人妻夜夜爽99麻豆av| 亚洲国产精品合色在线| 国产伦在线观看视频一区| 精品久久久久久久久久免费视频| 蜜桃亚洲精品一区二区三区| 少妇裸体淫交视频免费看高清| 日韩欧美一区二区三区在线观看| 成人三级黄色视频| 国产三级在线视频| 我的老师免费观看完整版| 99国产精品一区二区三区| 两个人看的免费小视频| 国产日本99.免费观看| 麻豆成人午夜福利视频| 日本一二三区视频观看| 欧美在线一区亚洲| 丝袜美腿在线中文| 亚洲av日韩精品久久久久久密| 亚洲精品日韩av片在线观看 | 国产真人三级小视频在线观看| 精品久久久久久,| 午夜免费成人在线视频| 999久久久精品免费观看国产| 国产午夜精品论理片| 国产色婷婷99| 天堂影院成人在线观看| 蜜桃久久精品国产亚洲av| 又爽又黄无遮挡网站| bbb黄色大片| a在线观看视频网站| 欧美日韩乱码在线| 噜噜噜噜噜久久久久久91| 国产亚洲欧美在线一区二区| 国产黄a三级三级三级人| 九色成人免费人妻av| 少妇的逼好多水| 免费电影在线观看免费观看| 亚洲自拍偷在线| 女人高潮潮喷娇喘18禁视频| 波多野结衣高清无吗| 日本黄大片高清| 婷婷精品国产亚洲av在线| 欧美日韩黄片免| 少妇高潮的动态图| 国产激情欧美一区二区| 日本一本二区三区精品| 亚洲在线观看片| 啦啦啦观看免费观看视频高清| 尤物成人国产欧美一区二区三区| bbb黄色大片| 亚洲欧美日韩高清专用| 欧美成人一区二区免费高清观看| 欧美在线一区亚洲| av视频在线观看入口| 高清毛片免费观看视频网站| 日韩av在线大香蕉| 欧美性猛交╳xxx乱大交人| 少妇的逼好多水| 久久香蕉国产精品| 乱人视频在线观看| av视频在线观看入口| 午夜老司机福利剧场| 露出奶头的视频| 亚洲内射少妇av| 18禁国产床啪视频网站| 色在线成人网| 亚洲欧美日韩无卡精品| avwww免费| 男女做爰动态图高潮gif福利片| 久久精品国产亚洲av香蕉五月| 精品国产三级普通话版| 啦啦啦免费观看视频1| 波多野结衣高清无吗| 最好的美女福利视频网| 亚洲成人精品中文字幕电影| 久久精品夜夜夜夜夜久久蜜豆| 热99在线观看视频| 久久这里只有精品中国| 亚洲精品一区av在线观看| 悠悠久久av| 欧美+亚洲+日韩+国产| 久久人人精品亚洲av| 国产精品香港三级国产av潘金莲| 久久精品国产亚洲av香蕉五月| 免费电影在线观看免费观看| 99久久精品热视频| 最好的美女福利视频网| 中文字幕人成人乱码亚洲影| 国产伦精品一区二区三区视频9 | 久久久成人免费电影| 五月伊人婷婷丁香| 亚洲片人在线观看| 男女做爰动态图高潮gif福利片| 久久精品国产清高在天天线| 床上黄色一级片| 国产一区二区三区在线臀色熟女| 一个人免费在线观看电影| 欧美黑人巨大hd| 俺也久久电影网| 国产精品久久视频播放| 欧美成狂野欧美在线观看| 午夜激情欧美在线| 少妇的丰满在线观看| 国内久久婷婷六月综合欲色啪| 国产 一区 欧美 日韩| 亚洲人成网站高清观看| 久久久久久久亚洲中文字幕 | 成年女人毛片免费观看观看9| 色尼玛亚洲综合影院| 日韩欧美国产在线观看| av专区在线播放| 欧美中文综合在线视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲在线自拍视频| 精品久久久久久久人妻蜜臀av| 99国产精品一区二区三区| 一个人看的www免费观看视频| 日本黄大片高清| 精品电影一区二区在线| 国内精品美女久久久久久| 国产高潮美女av| 久久国产乱子伦精品免费另类| 男女下面进入的视频免费午夜| 久久精品国产清高在天天线| 天天添夜夜摸| 亚洲欧美一区二区三区黑人| 最近最新中文字幕大全免费视频| 午夜福利在线观看吧| 国产成人aa在线观看| 老司机午夜福利在线观看视频| 最近最新中文字幕大全电影3| 小蜜桃在线观看免费完整版高清| 国产伦一二天堂av在线观看| 国产伦在线观看视频一区| 亚洲va日本ⅴa欧美va伊人久久| 欧美区成人在线视频| 国产爱豆传媒在线观看| 欧美av亚洲av综合av国产av| 欧美午夜高清在线| 丰满乱子伦码专区| 国产精品日韩av在线免费观看| 欧美性猛交黑人性爽| 69人妻影院| 桃红色精品国产亚洲av| 免费在线观看亚洲国产| 中文字幕人妻熟人妻熟丝袜美 | av中文乱码字幕在线| 精品久久久久久,| 久久精品国产亚洲av香蕉五月| 九色成人免费人妻av| 在线观看美女被高潮喷水网站 | 午夜免费成人在线视频| 亚洲片人在线观看| 88av欧美| 欧美性猛交╳xxx乱大交人| 最后的刺客免费高清国语| 一级黄片播放器| av黄色大香蕉| 亚洲国产欧美网| 亚洲av二区三区四区| 亚洲电影在线观看av| 男人舔女人下体高潮全视频| 亚洲精品久久国产高清桃花| 亚洲第一欧美日韩一区二区三区| 在线播放国产精品三级| 在线播放无遮挡| 午夜两性在线视频| 日韩欧美精品v在线| 床上黄色一级片| 最近最新中文字幕大全免费视频| 国产成人av教育| 久久亚洲真实| 青草久久国产| 在线观看美女被高潮喷水网站 | 午夜福利高清视频| 日本 av在线| 亚洲无线在线观看| 女人十人毛片免费观看3o分钟| 日韩高清综合在线| 可以在线观看的亚洲视频| 亚洲不卡免费看| 一本综合久久免费| 欧美一区二区精品小视频在线| 色综合亚洲欧美另类图片| 欧美一区二区国产精品久久精品| 少妇裸体淫交视频免费看高清| 一本久久中文字幕| 国产69精品久久久久777片| 久久久精品大字幕| 日韩精品中文字幕看吧| 欧美又色又爽又黄视频| 午夜a级毛片| 国产亚洲欧美在线一区二区| 亚洲欧美日韩卡通动漫| 在线观看日韩欧美| 国产久久久一区二区三区| 无限看片的www在线观看| 好男人电影高清在线观看| 丰满人妻一区二区三区视频av | 国产97色在线日韩免费| 久久国产精品人妻蜜桃| 国产免费男女视频| 国内精品美女久久久久久| 国产精品av视频在线免费观看| 免费大片18禁| 可以在线观看的亚洲视频| 午夜久久久久精精品| 特级一级黄色大片| 亚洲精品亚洲一区二区| 热99在线观看视频| 国内毛片毛片毛片毛片毛片| 99在线视频只有这里精品首页| 国产在视频线在精品| 白带黄色成豆腐渣| 国产成人系列免费观看| a级毛片a级免费在线| 亚洲成人久久爱视频| 国产高清三级在线| 9191精品国产免费久久| 日本成人三级电影网站| 色噜噜av男人的天堂激情| 欧美性感艳星| 亚洲成人久久爱视频| 国内精品久久久久精免费| 精品不卡国产一区二区三区| 成人三级黄色视频| 中文字幕人妻丝袜一区二区| 丝袜美腿在线中文| 午夜激情福利司机影院| 成人欧美大片| 人妻夜夜爽99麻豆av| 在线免费观看的www视频| 国产精品久久久久久人妻精品电影| 久久久国产成人免费| 亚洲国产中文字幕在线视频| 国产午夜精品论理片| 国产一区二区三区视频了| 精品国产美女av久久久久小说| 91在线精品国自产拍蜜月 | 全区人妻精品视频| 高清日韩中文字幕在线| 男女视频在线观看网站免费| 麻豆国产av国片精品| 国产精品久久电影中文字幕| 国产av麻豆久久久久久久| 亚洲成人久久性| 亚洲18禁久久av| 欧美乱妇无乱码| 欧美大码av| 中文资源天堂在线| 一个人看视频在线观看www免费 | 69av精品久久久久久| www国产在线视频色| 又黄又粗又硬又大视频| 亚洲美女黄片视频| 国内毛片毛片毛片毛片毛片| tocl精华| 欧洲精品卡2卡3卡4卡5卡区| 国产高清三级在线| 老司机深夜福利视频在线观看| 精品人妻偷拍中文字幕| 国产91精品成人一区二区三区| 亚洲av熟女| 精品久久久久久成人av| 欧美日韩精品网址| 国产三级在线视频| www.999成人在线观看| 免费看美女性在线毛片视频| 久久九九热精品免费| 久久久久久久亚洲中文字幕 | 欧美乱码精品一区二区三区| 久久久久久九九精品二区国产| 国产高清视频在线观看网站| 偷拍熟女少妇极品色| 欧美色欧美亚洲另类二区| 色综合婷婷激情| 91麻豆精品激情在线观看国产| 桃红色精品国产亚洲av| 国产精品一区二区免费欧美| 亚洲成人久久性| a级毛片a级免费在线| 最新在线观看一区二区三区| 日韩国内少妇激情av| 最新中文字幕久久久久| 99国产精品一区二区蜜桃av| 亚洲一区高清亚洲精品| 亚洲天堂国产精品一区在线| 又黄又爽又免费观看的视频| 国内揄拍国产精品人妻在线| 亚洲aⅴ乱码一区二区在线播放| 婷婷丁香在线五月| 亚洲色图av天堂| 香蕉av资源在线| 青草久久国产| 国内精品久久久久精免费| 午夜福利在线观看免费完整高清在 | 精品人妻1区二区| 男女之事视频高清在线观看| 亚洲国产精品久久男人天堂| 变态另类成人亚洲欧美熟女| 啦啦啦免费观看视频1| 岛国在线免费视频观看| 日本在线视频免费播放| 精品一区二区三区视频在线观看免费| 亚洲熟妇熟女久久| 熟女人妻精品中文字幕| 亚洲一区二区三区不卡视频| 婷婷精品国产亚洲av在线| 女人十人毛片免费观看3o分钟| 欧美黄色淫秽网站| 手机成人av网站| 丰满的人妻完整版| 久久午夜亚洲精品久久| 啦啦啦观看免费观看视频高清| av天堂在线播放| 啪啪无遮挡十八禁网站| xxx96com| 国产av在哪里看| 一区二区三区免费毛片| 亚洲在线观看片| 中文字幕av成人在线电影| 亚洲无线在线观看| 久久精品国产99精品国产亚洲性色| 老司机在亚洲福利影院| 丰满人妻一区二区三区视频av | 欧美国产日韩亚洲一区| 色哟哟哟哟哟哟| 国产精品精品国产色婷婷| 校园春色视频在线观看| 国产老妇女一区| 桃红色精品国产亚洲av| 两人在一起打扑克的视频| 亚洲av五月六月丁香网| 欧美日韩中文字幕国产精品一区二区三区| 国产精品美女特级片免费视频播放器| 无限看片的www在线观看| 国产蜜桃级精品一区二区三区| 麻豆成人午夜福利视频| 一个人看视频在线观看www免费 | 夜夜夜夜夜久久久久| 九九在线视频观看精品| 九九热线精品视视频播放| 99热这里只有精品一区| 亚洲中文字幕一区二区三区有码在线看| 国产亚洲欧美在线一区二区| 一卡2卡三卡四卡精品乱码亚洲| av天堂中文字幕网| 国产三级黄色录像| 亚洲男人的天堂狠狠| av在线蜜桃| 国产高清有码在线观看视频| 国产成人a区在线观看| 神马国产精品三级电影在线观看| 免费搜索国产男女视频| 可以在线观看毛片的网站| 一级a爱片免费观看的视频| 国产亚洲av嫩草精品影院| 在线播放无遮挡| 最新中文字幕久久久久| 成人性生交大片免费视频hd| 国产精品乱码一区二三区的特点| 亚洲av免费在线观看| 首页视频小说图片口味搜索| 亚洲中文字幕日韩| 免费看美女性在线毛片视频| 精品一区二区三区av网在线观看| 久久人人精品亚洲av| 午夜免费激情av| 又黄又粗又硬又大视频| 99热只有精品国产| 亚洲avbb在线观看| 一区福利在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲中文字幕日韩| 免费观看人在逋| 一级毛片女人18水好多| 在线观看一区二区三区| 国产麻豆成人av免费视频| 国产黄片美女视频| 久久久久久久午夜电影| svipshipincom国产片| 又粗又爽又猛毛片免费看| 99国产精品一区二区蜜桃av| 免费观看精品视频网站| 一a级毛片在线观看| 在线免费观看不下载黄p国产 | 国产一区二区在线av高清观看| 亚洲精品456在线播放app | 天天添夜夜摸| 国产色爽女视频免费观看| 村上凉子中文字幕在线| 1024手机看黄色片| h日本视频在线播放| 亚洲av电影在线进入| 国产色爽女视频免费观看| 亚洲精品日韩av片在线观看 | 法律面前人人平等表现在哪些方面| 中文字幕人成人乱码亚洲影| 欧美日韩一级在线毛片| 免费在线观看影片大全网站| 制服丝袜大香蕉在线| 国产精品野战在线观看| 他把我摸到了高潮在线观看| 性欧美人与动物交配| 午夜福利高清视频| 一本久久中文字幕| 国产精品国产高清国产av| 国产一区二区三区在线臀色熟女| 国产在线精品亚洲第一网站| av欧美777| 成人鲁丝片一二三区免费| 国产高清视频在线播放一区| 欧美xxxx黑人xx丫x性爽| 中文资源天堂在线| 99在线人妻在线中文字幕| 国内久久婷婷六月综合欲色啪| 色在线成人网| 少妇人妻精品综合一区二区 | 在线观看一区二区三区| 国产精品美女特级片免费视频播放器| 国产精品99久久99久久久不卡| 超碰av人人做人人爽久久 | 国产成+人综合+亚洲专区| 亚洲内射少妇av| 国产色爽女视频免费观看| 婷婷六月久久综合丁香| 精品午夜福利视频在线观看一区| 又粗又爽又猛毛片免费看| 成年人黄色毛片网站| 少妇高潮的动态图| 亚洲av一区综合| 亚洲国产欧美网| 麻豆国产97在线/欧美| 久久人妻av系列| 久久精品国产亚洲av香蕉五月| 99久久综合精品五月天人人| 中文字幕av成人在线电影| 18+在线观看网站| 亚洲精品色激情综合| www.熟女人妻精品国产| 国产不卡一卡二| 国产69精品久久久久777片| 午夜福利高清视频| 给我免费播放毛片高清在线观看| 亚洲人成网站在线播放欧美日韩| 久久亚洲真实| 婷婷亚洲欧美| 少妇裸体淫交视频免费看高清| 欧美日本视频| 成人特级av手机在线观看| 一级毛片女人18水好多| 男女床上黄色一级片免费看| 日韩欧美国产在线观看| 草草在线视频免费看| 欧美成人免费av一区二区三区| 午夜免费成人在线视频| 国产精品美女特级片免费视频播放器| 精品一区二区三区人妻视频| 午夜福利视频1000在线观看| 欧美区成人在线视频| 一边摸一边抽搐一进一小说| 国语自产精品视频在线第100页| 欧美激情在线99| 日韩精品中文字幕看吧| 国产伦人伦偷精品视频| 亚洲一区二区三区不卡视频| 欧美性猛交╳xxx乱大交人| 中文字幕高清在线视频| 女生性感内裤真人,穿戴方法视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久久久久久午夜电影| 国产精品永久免费网站| 国产精品久久久人人做人人爽| 日本三级黄在线观看| 日本 av在线| 黄色丝袜av网址大全| 男女做爰动态图高潮gif福利片| 校园春色视频在线观看| 午夜福利18| 久久久国产成人精品二区| 2021天堂中文幕一二区在线观| 国内精品久久久久久久电影| ponron亚洲| 久久久国产成人免费| 日韩欧美 国产精品| 在线观看免费午夜福利视频| 国内精品美女久久久久久| 国产欧美日韩一区二区精品| 三级国产精品欧美在线观看| 最近在线观看免费完整版| 一级黄色大片毛片|