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

    近22 a人為氣溶膠對(duì)東亞地區(qū)夏季風(fēng)環(huán)流影響的數(shù)值模擬研究

    2014-09-09 07:30:09陳明誠(chéng)孫照渤沈新勇王勇吉振明
    大氣科學(xué)學(xué)報(bào) 2014年4期

    陳明誠(chéng),孫照渤,沈新勇,王勇,吉振明

    (1.南京信息工程大學(xué) 氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210044;2.中國(guó)科學(xué)院 大氣物理研究所云降水物理與強(qiáng)風(fēng)暴重點(diǎn)實(shí)驗(yàn)室,北京 100029;3.中國(guó)科學(xué)院 青藏高原研究所 青藏高原環(huán)境變化與地表過(guò)程重點(diǎn)實(shí)驗(yàn)室,北京 100085)

    近22a人為氣溶膠對(duì)東亞地區(qū)夏季風(fēng)環(huán)流影響的數(shù)值模擬研究

    陳明誠(chéng)1,2,孫照渤1,沈新勇1,2,王勇1,吉振明3

    (1.南京信息工程大學(xué) 氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210044;2.中國(guó)科學(xué)院 大氣物理研究所云降水物理與強(qiáng)風(fēng)暴重點(diǎn)實(shí)驗(yàn)室,北京 100029;3.中國(guó)科學(xué)院 青藏高原研究所 青藏高原環(huán)境變化與地表過(guò)程重點(diǎn)實(shí)驗(yàn)室,北京 100085)

    運(yùn)用區(qū)域氣候模式RegCM3耦合入一個(gè)化學(xué)過(guò)程,對(duì)東亞地區(qū)三類(lèi)人為排放氣溶膠(硫酸鹽、黑碳和有機(jī)碳)的時(shí)空分布特征及其對(duì)夏季風(fēng)環(huán)流的影響進(jìn)行了數(shù)值模擬研究。模擬結(jié)果顯示,氣溶膠的引入會(huì)引起東亞地區(qū)夏季850 hPa風(fēng)場(chǎng)發(fā)生改變,我國(guó)江淮以東洋面上空出現(xiàn)了一個(gè)氣旋式距平環(huán)流中心,中心以西的偏北風(fēng)氣流將削弱東亞地區(qū)夏季西南季風(fēng)。通過(guò)討論春季中國(guó)地區(qū)氣溶膠濃度與夏季東亞地區(qū)850 hPa經(jīng)向風(fēng)的時(shí)滯關(guān)系,以及夏季中國(guó)地區(qū)氣溶膠濃度與同期東亞地區(qū)850 hPa經(jīng)向風(fēng)的關(guān)系,可以發(fā)現(xiàn),春、夏季中國(guó)地區(qū)氣溶膠濃度均與夏季東亞地區(qū)850 hPa經(jīng)向風(fēng)有很好的負(fù)相關(guān)關(guān)系,當(dāng)春季中國(guó)北方和夏季中國(guó)南方地區(qū)氣溶膠濃度增加時(shí),中國(guó)東部地區(qū)夏季偏南季風(fēng)減弱。這可能與氣溶膠改變了大氣層頂和地表的輻射強(qiáng)迫,進(jìn)而引起了海陸氣壓差異和位勢(shì)高度場(chǎng)的變化有關(guān)。

    人為氣溶膠;東亞夏季風(fēng);季風(fēng)環(huán)流;數(shù)值模擬;相關(guān)分析

    0 引言

    我國(guó)地處東亞季風(fēng)區(qū),由于東亞季風(fēng)的年際和年代際變化較大,因此夏季旱澇和冬季寒潮等重大氣候?yàn)?zāi)害在我國(guó)發(fā)生頻繁,所造成的經(jīng)濟(jì)損失和人員傷亡往往也十分嚴(yán)重,而這些嚴(yán)重的氣候?yàn)?zāi)害與東亞季風(fēng)系統(tǒng)的異常變化間有著密切的聯(lián)系。黃榮輝等(2003)在研究了我國(guó)主要的氣候?yàn)?zāi)害時(shí)空分布特征后指出,近20余年來(lái)華北地區(qū)春、夏季極易發(fā)生干旱,長(zhǎng)江流域及華南地區(qū)洪澇災(zāi)害頻發(fā)。韓晉平和王會(huì)軍(2007)利用ERA和NCEP/NCAR再分析資料以及部分觀測(cè)資料研究后同樣發(fā)現(xiàn),東亞夏季風(fēng)在20世紀(jì)80年代以來(lái)減弱這一現(xiàn)象普遍存在,主要表現(xiàn)為中國(guó)東部地區(qū)出現(xiàn)偏北風(fēng)異常。相關(guān)研究表明,近20 a來(lái)工業(yè)化、城市化等人類(lèi)活動(dòng)發(fā)展迅速的亞洲地區(qū)正逐漸成為氣溶膠排放的高值區(qū)(Hansen et al.,1998;沈志寶和魏麗,1999;王喜紅等,2002;Chung et al.,2002;陳麗和銀燕,2008;鐘昌琴等,2008)。由于氣溶膠在改變大氣及地表能量平衡的同時(shí)影響到云微物理過(guò)程,并進(jìn)而對(duì)東亞季風(fēng)氣候系統(tǒng)產(chǎn)生重要影響,而氣溶膠與東亞季風(fēng)系統(tǒng)間的相互作用依然具有很大的未知性,因此更好地理解兩者間的相互關(guān)系具有很大的科學(xué)和社會(huì)利益。

    大氣中的氣溶膠組成成分比較復(fù)雜,但大量研究表明硫酸鹽、黑碳和有機(jī)碳這三類(lèi)人為排放氣溶膠對(duì)亞洲季風(fēng)系統(tǒng)正起著越來(lái)越重要的影響作用(孫家仁和劉煜,2008a,2008b;王志立等,2009;吉振明等,2010)。關(guān)于亞洲季風(fēng)異常變化與氣溶膠的相互作用,國(guó)內(nèi)外學(xué)者通過(guò)外場(chǎng)觀察和數(shù)值模擬方法已進(jìn)行了一系列的研究并取得許多具有應(yīng)用價(jià)值的結(jié)論:Massimo et al.(2008)對(duì)南亞地區(qū)夏季風(fēng)演變與氣溶膠的大尺度不規(guī)則分布研究中揭示了南亞夏季風(fēng)的增強(qiáng)與季風(fēng)爆發(fā)前期氣溶膠的大規(guī)模分布間存在著時(shí)空響應(yīng)關(guān)系。Leon and Ulrike(2002)在研究了熱帶降水與氣溶膠的間接影響作用后指出氣溶膠分布的時(shí)空變化可對(duì)大氣低層緯向環(huán)流及降水產(chǎn)生實(shí)質(zhì)性的影響。Xu(2001)認(rèn)為從20世紀(jì)90年代末中國(guó)夏季連續(xù)出現(xiàn)的“南澇北旱”現(xiàn)象是由于人為排放硫化物氣溶膠的輻射影響下引起的雨帶南移造成的。在氣溶膠對(duì)東亞夏季風(fēng)環(huán)流的影響方面。王志立等(2009)通過(guò)模擬南亞地區(qū)黑碳?xì)馊苣z對(duì)亞洲夏季風(fēng)氣候效應(yīng)的影響時(shí)發(fā)現(xiàn)南亞地區(qū)黑碳?xì)馊苣z減弱了東亞夏季風(fēng)。孫家仁和劉煜(2008a,2008b)利用新一代GCM CAM3.0模式離線耦合一個(gè)氣溶膠同化系統(tǒng)分別模擬了中國(guó)區(qū)域內(nèi)硫酸鹽氣溶膠和黑碳?xì)馊苣z對(duì)東亞季風(fēng)的可能影響后進(jìn)一步指出,要綜合考慮印度半島及其周邊地區(qū)與中國(guó)區(qū)域氣溶膠對(duì)東亞夏季風(fēng)強(qiáng)度的影響作用。吉振明等(2010)使用區(qū)域氣候模式RegCM3,通過(guò)長(zhǎng)時(shí)間連續(xù)積分模擬了亞洲地區(qū)主要?dú)馊苣z的時(shí)空分布,結(jié)果表明模式能夠較好地模擬中國(guó)地區(qū)氣溫和降水的分布,對(duì)該地區(qū)氣溶膠空間分布也具有一定的模擬能力。同時(shí),研究還發(fā)現(xiàn)氣溶膠對(duì)地面氣溫和降水有著重要的影響作用。

    相關(guān)研究表明,在模擬東亞季風(fēng)氣候時(shí)應(yīng)選取具有高分辨率的區(qū)域氣候模式(高學(xué)杰等,2006;Zhou and Yu,2006;Gao et al.,2008;吉振明等,2010)。本次試驗(yàn)中,使用耦合入化學(xué)過(guò)程的區(qū)域氣候模式RegCM3,在引入三類(lèi)主要人為排放氣溶膠(硫酸鹽、黑碳和有機(jī)碳)后,對(duì)亞洲地區(qū)季風(fēng)環(huán)流系統(tǒng)進(jìn)行長(zhǎng)時(shí)間的模擬,模擬結(jié)果顯示了三類(lèi)氣溶膠各自的輸送特征及綜合影響下的輻射強(qiáng)迫效果。最后,在前人的基礎(chǔ)之上,運(yùn)用相關(guān)性分析和合成差值檢驗(yàn)等統(tǒng)計(jì)學(xué)方法,對(duì)中國(guó)地區(qū)不同區(qū)域人為氣溶膠排放量與東亞夏季風(fēng)環(huán)流間的相關(guān)關(guān)系進(jìn)行分析研究,以期望找出氣溶膠對(duì)東亞夏季風(fēng)的可能影響,并得出前期及同期的預(yù)報(bào)信號(hào)。

    1 模擬試驗(yàn)方案

    使用的模型為意大利國(guó)際理論物理中心區(qū)域氣候模式RegCM3(Pal et al.,2007),該模式耦合了Solmon et al.(2006)在Qian et al.(2001)方案基礎(chǔ)上發(fā)展出的氣溶膠化學(xué)過(guò)程,對(duì)氣溶膠在大氣中的傳輸過(guò)程進(jìn)行了更精確的模擬。計(jì)算氣溶膠輻射效應(yīng)時(shí),采用Kiehl et al.(1996)的計(jì)算方法,對(duì)模式劃分出的19個(gè)波段(0.2~5.0 μm)分別給予不同的光學(xué)特性,以期得到更為全面的結(jié)果。

    為方便于分析人為氣溶膠對(duì)東亞地區(qū)夏季風(fēng)環(huán)流的綜合影響效果,本次模擬運(yùn)算中共設(shè)計(jì)兩套試驗(yàn)方案:不包含氣溶膠的控制試驗(yàn)(CON)和包含了三類(lèi)人為排放氣溶膠(硫酸鹽、黑碳和有機(jī)碳)的敏感性試驗(yàn)(SEN)。在控制性試驗(yàn)中,模式的初始場(chǎng)以及側(cè)邊界場(chǎng)由NCAR/NCEP再分析數(shù)據(jù)驅(qū)動(dòng),海溫采用NOAA提供的OISST資料(Reynolds et al.,2002)。為考慮模擬區(qū)域包含中國(guó)、中南半島、印度半島及其周邊的亞洲大部分地區(qū),將模擬區(qū)域中心點(diǎn)取為(99.0°E,32.0°N),東西向和南北向格點(diǎn)數(shù)分別為192×160,水平分辨率取為50 km。在敏感性試驗(yàn)中,在前者的基礎(chǔ)之上引入了三類(lèi)人為排放氣溶膠,用于突出氣溶膠的氣候效應(yīng),氣溶膠模擬所采用的排放清單為日本全球變化前沿研究中心制作的REAS資料(Ohara et al.,2007),該排放清單分辨率為0.5°×0.5°。模擬過(guò)程中參數(shù)化方案的選擇可參閱吉振明等(2010)的研究成果。

    圖1 亞洲地區(qū)夏季平均200 hPa緯向風(fēng)場(chǎng)(a、b)和850 hPa水平風(fēng)場(chǎng)(c、d)的觀測(cè)(a、c)與模擬(b、d)結(jié)果(單位:m/s)Fig.1 (a,c)Observed and (b,d)simulated summer mean (a,b)200 hPa zonal wind field and (c,d)850 hPa horizontal wind field over Asia(units:m/s)

    兩套模擬試驗(yàn)積分時(shí)間均自1987年1月1日起至2009年12月31日止,其中1987年的模擬結(jié)果作為模擬初始化階段不予分析,選取其后的22 a作為分析時(shí)段,并將敏感性試驗(yàn)與控制性試驗(yàn)結(jié)果之差(SEN-CON)作為三類(lèi)人為氣溶膠的綜合氣候效應(yīng)進(jìn)行討論分析。

    2 模擬結(jié)果分析

    2.1 模擬性能檢驗(yàn)

    在對(duì)模擬結(jié)果進(jìn)行分析之前,首先對(duì)模式的模擬效果進(jìn)行檢驗(yàn),將控制性試驗(yàn)的結(jié)果與NCEP再分析結(jié)果做對(duì)比分析。圖1分別給出了亞洲地區(qū)夏季22 a平均200 hPa緯向風(fēng)場(chǎng)及850 hPa水平風(fēng)場(chǎng)的模擬和檢驗(yàn),通過(guò)對(duì)比可以看出,夏季亞洲200 hPa上中高緯地區(qū)以偏西風(fēng)氣流為主,低緯地區(qū)則轉(zhuǎn)為偏東風(fēng)氣流,緯向風(fēng)速零值線約位于25°N附近(圖1a、b)。在夏季對(duì)流層中低層850 hPa水平風(fēng)場(chǎng)中(圖1c、d),南亞地區(qū)盛行較強(qiáng)的偏西氣流,夏季風(fēng)槽區(qū)位于印度半島東側(cè),東亞地區(qū)則以偏西南氣流的經(jīng)向北上為主??傮w而言,模式模擬結(jié)果與再分析資料基本保持一致,能夠模擬出亞洲夏季風(fēng)區(qū)的基本氣候特征。

    2.2 氣溶膠的傳輸機(jī)制

    近22 a來(lái),亞洲地區(qū)作為氣溶膠排放量增長(zhǎng)速度最快的地區(qū)之一,其空間分布的差異對(duì)區(qū)域大氣環(huán)流的改變起著十分重要的作用。為了了解模式在東亞夏季風(fēng)爆發(fā)前期(春季)及爆發(fā)過(guò)程中(夏季)對(duì)氣溶膠分布狀況和水平風(fēng)場(chǎng)影響下其輸送特征的模擬效果,分別對(duì)春、夏季亞洲地區(qū)三類(lèi)氣溶膠的柱濃度含量分布特征進(jìn)行了模擬研究。

    圖2 春季(a、c、e)、夏季(b、d、f)氣溶膠平均柱濃度含量(陰影;單位:mg/m2)及850 hPa水平風(fēng)場(chǎng)(箭矢;單位:m/s) a,b.硫酸鹽;c,d.黑碳;e,f.有機(jī)碳Fig.2 The (a,c,e)spring and (b,d,f)summer aerosol mean column burden(shaded areas;units:mg/m2) and 850 hPa wind field (arrows;units:m/s) a,b.sulfate;c,d.black carbon;e,f.organic carbon

    春季,三類(lèi)人為排放氣溶膠(硫酸鹽、黑碳和有機(jī)碳)的分布情況與REAS資料(Ohara et al.,2007)給出的結(jié)果基本保持一致(圖略),印度半島地區(qū)的氣溶膠柱濃度含量中心區(qū)出現(xiàn)在東部沿海地區(qū),而在東亞地區(qū)氣溶膠的濃度高值區(qū)主要位于我國(guó)的南方和四川盆地地區(qū)(圖2a、c、e)。結(jié)合亞洲地區(qū)春季850 hPa風(fēng)場(chǎng)可以看出,東亞地區(qū)西太副高脊線呈東北—西南走向(圖略),高壓脊南側(cè)的低緯地區(qū)有較強(qiáng)勁的偏東氣流自海洋吹向大陸,這支偏東氣流在經(jīng)過(guò)南海上空時(shí)分為兩支:一支于副高西側(cè)轉(zhuǎn)為偏南氣流北上進(jìn)入我國(guó);另一支繼續(xù)向西前進(jìn),在孟加拉灣海域北上,同印度半島地區(qū)的偏西氣流合并為西南氣流,經(jīng)中南半島進(jìn)入我國(guó)。兩支偏南氣流于我國(guó)南方地區(qū)匯合北上,至長(zhǎng)江流域與來(lái)自高緯的偏北氣流相遇,轉(zhuǎn)為向東傳播。在這樣的大氣環(huán)流背景下,推測(cè)印度半島東部地區(qū)的濃度中心主要為半島中北部工業(yè)排放區(qū)大量排放人為氣溶膠和偏西氣流輸送共同導(dǎo)致的結(jié)果,而我國(guó)南方地區(qū)由于自身氣溶膠排放量較高,來(lái)自低緯的偏南氣流輸送則有可能進(jìn)一步促進(jìn)該地區(qū)氣溶膠的堆積。四川盆地地區(qū)由于高山、丘陵等地形較多,在地形的輻合作用影響下污染物不易擴(kuò)散,西南氣流輸送而來(lái)的氣溶膠易于在此沉積,從而形成了東亞地區(qū)另一個(gè)氣溶膠含量的高值中心。綜合來(lái)看,較高的氣溶膠排放量和區(qū)域的輻合作用(包括地形作用和風(fēng)場(chǎng)自身的輻合輻散效果)是形成春季亞洲地區(qū)人為氣溶膠濃度高值區(qū)的兩個(gè)主要原因。

    到了夏季,三類(lèi)人為氣溶膠的空間分布特征與春季有所不同。中國(guó)東部、中南半島以及印度半島的北部地區(qū)均有分布,高濃度中心區(qū)域出現(xiàn)在四川盆地、中國(guó)北方地區(qū)以及青藏高原南側(cè)、印度北部的狹長(zhǎng)地帶(圖2b、d、f)。對(duì)比春季的氣溶膠分布情況發(fā)現(xiàn),在東亞地區(qū)氣溶膠有明顯的向北輸送的特征,我國(guó)長(zhǎng)江以北的華北大部分地區(qū)出現(xiàn)了與四川盆地地區(qū)同量級(jí)的高值中心,在氣溶膠區(qū)域排放量季節(jié)性變化較小時(shí),這種空間分布的季節(jié)性差異有可能由夏季風(fēng)盛行風(fēng)向的變化所導(dǎo)致。在夏季850 hPa風(fēng)場(chǎng)中,東亞夏季風(fēng)建立,低層以西南氣流的經(jīng)向北上為主。自南向北的偏南風(fēng)氣流伴隨著大量的氣溶膠向北輸送,因此夏季中國(guó)北方地區(qū)氣溶膠柱濃度含量較春季有了顯著地升高。

    2.3 氣溶膠對(duì)東亞地區(qū)夏季風(fēng)的影響

    大氣中氣溶膠濃度變化所引起的地氣系統(tǒng)能量平衡間產(chǎn)生的擾動(dòng)稱(chēng)之為“輻射強(qiáng)迫”。不同種類(lèi)氣溶膠粒子由于自身光學(xué)特性的不同,對(duì)地氣系統(tǒng)產(chǎn)生的輻射強(qiáng)迫效應(yīng)也會(huì)有所不同,硫酸鹽氣溶膠以散射和反射大氣長(zhǎng)波輻射及太陽(yáng)短波輻射為主,到達(dá)地面的太陽(yáng)輻射減弱,因此對(duì)整層及地表均產(chǎn)生負(fù)輻射強(qiáng)迫,從而起到增加大氣靜力穩(wěn)定度,穩(wěn)定大氣層結(jié)的效果,而黑碳?xì)馊苣z在反射的同時(shí)也會(huì)吸收一部分太陽(yáng)輻射,對(duì)大氣起到加熱作用。實(shí)際情況中的氣溶膠輻射強(qiáng)迫作用往往十分復(fù)雜,綜合考慮其混合后所產(chǎn)生的總輻射性質(zhì)或許能夠得到更真實(shí)的結(jié)果(孫家仁和劉煜,2008a,2008b)。圖3給出了模式模擬的春、夏季綜合了三類(lèi)氣溶膠后在大氣層頂及地表的短波輻射強(qiáng)迫分布情況,可以看出,綜合考慮了三類(lèi)人為排放氣溶膠之后的輻射效應(yīng)均為負(fù)值,春季大氣層頂負(fù)輻射效應(yīng)的中心位于四川盆地地區(qū),整個(gè)負(fù)輻射效應(yīng)區(qū)覆蓋了我國(guó)東部沿海地區(qū),對(duì)應(yīng)的地表負(fù)輻射效應(yīng)區(qū)范圍更廣,向北可延伸至朝鮮半島、日本等地,向南則擴(kuò)展到整個(gè)亞洲的低緯地區(qū),在孟加拉灣地區(qū)出現(xiàn)了負(fù)輻射強(qiáng)迫中心,中心強(qiáng)度-20 W/m2以上。夏季的短波輻射強(qiáng)迫效果在大氣層頂及地表較春季更為明顯,我國(guó)的華北和渤海地區(qū)強(qiáng)度有明顯的增強(qiáng),江淮地區(qū)出現(xiàn)了第二個(gè)極值中心,中心最大強(qiáng)度與四川盆地地區(qū)相同,可達(dá)-30 W/m2以上。

    氣溶膠的間接氣候效應(yīng)相較于直接效應(yīng)更為復(fù)雜,由于氣溶膠粒子可以通過(guò)凝結(jié)核作用,在云形成過(guò)程中參與到云的生成、演化、消散過(guò)程中每一環(huán)節(jié),使得空氣中顆粒濃度增加,云滴半徑減小,改變?cè)频蔚脑莆⑽锢斫Y(jié)構(gòu)和光學(xué)特性,達(dá)到調(diào)整空氣中液態(tài)水含量及云層厚度的效果,并最終影響到云的生命期和降水率,因此氣溶膠的間接效應(yīng)充滿(mǎn)著不確定性,本次試驗(yàn)中只考慮氣溶膠的直接氣候效應(yīng)。

    氣溶膠的輻射強(qiáng)迫效應(yīng)將引起自地表至高層中大氣運(yùn)動(dòng)的變化,為此對(duì)氣溶膠引起的地表氣壓及大氣中位勢(shì)高度場(chǎng)變化進(jìn)行分析。圖4a、b分別給出了春、夏季東亞地區(qū)敏感性試驗(yàn)與控制實(shí)驗(yàn)結(jié)果中地表氣壓之差,即氣溶膠對(duì)地表氣壓的影響效果,結(jié)果顯示:中國(guó)東部沿海地區(qū)地表氣壓普遍降低,夏季降壓趨勢(shì)更加顯著,華北和華南地區(qū)為降壓中心區(qū),最大降壓可達(dá)0.6 hPa。因此,海陸氣壓差異將有所減小,氣溶膠對(duì)東亞夏季風(fēng)會(huì)起到一定的削弱作用。

    圖3 春季(a、c)和夏季(b、d)模擬的大氣層頂(a、b)及地面(c、d)平均短波輻射強(qiáng)迫(單位:W/m2)Fig.3 Simulated (a,c)spring and (b,d)summer mean short wave radiative forcing on (a,b)the top of atmosphere and (c,d)the surface(units:W/m2)

    高空中,位勢(shì)高度場(chǎng)的變化與地表氣壓變化間有著很好的對(duì)應(yīng)關(guān)系。圖4c、d、e、f分別給出了春、夏季850 hPa和500 hPa位勢(shì)高度場(chǎng)的變化情況,春季我國(guó)北方至日本上空為正位勢(shì)高度距平,中心位于朝鮮半島上空,500 hPa上最大正位勢(shì)高度距平為1.5 gpm,同時(shí)發(fā)現(xiàn),在(130°E,50°N)地區(qū)以北出現(xiàn)了負(fù)位勢(shì)高度距平中心,結(jié)合春季大氣環(huán)流形勢(shì)(圖2a),這種南北向位勢(shì)高度距平分布將加強(qiáng)高緯地區(qū)對(duì)流層中的西風(fēng)氣流,偏強(qiáng)的歐亞西風(fēng)將阻礙夏季中國(guó)南方西南氣流北上,從而造成東亞夏季風(fēng)偏弱。夏季,朝鮮半島上空的中心轉(zhuǎn)變?yōu)樨?fù)位勢(shì)高度距平,并且中心略有南移至長(zhǎng)江出海口外,中心最大值超過(guò)了-2.5 gpm,這個(gè)負(fù)距平區(qū)域覆蓋了我國(guó)東南沿海地區(qū),同時(shí)在(110°E,25°N)附近有一正位勢(shì)高度中心生成,中心最大值達(dá)到了1.5 gpm。

    在夏季850 hPa風(fēng)場(chǎng)(圖1d)上,亞洲低緯地區(qū)盛行西風(fēng),經(jīng)過(guò)印度半島,孟加拉灣和中南半島后,在南海北部地區(qū)與越赤道氣流匯合轉(zhuǎn)為由南向北運(yùn)動(dòng)。因此,夏季我國(guó)東部地區(qū)以經(jīng)向氣流為主,在平坦的下墊面條件下盛行西南季風(fēng)。將敏感性試驗(yàn)得到的結(jié)果減去控制性試驗(yàn)的結(jié)果作為三類(lèi)人為排放氣溶膠所引起的氣候效應(yīng)(SEN-CON)。注意到在氣溶膠影響下(圖5),在我國(guó)東部江淮以東洋面上空出現(xiàn)了一個(gè)氣旋式距平環(huán)流,該環(huán)流中心位于長(zhǎng)江出海口外,在環(huán)流中心西側(cè)偏北氣流的影響下,我國(guó)東部地區(qū)夏季風(fēng)被削弱,出現(xiàn)了較為明顯的偏北風(fēng)距平,氣旋式環(huán)流南部的偏西氣流使得我國(guó)南方地區(qū)的西南氣流中西風(fēng)增強(qiáng),因此西太副高位置在氣溶膠的影響下可能較為偏南。張慶云等(2003)在研究了東亞夏季風(fēng)指數(shù)的年際變化與東亞大氣環(huán)流的關(guān)系時(shí)曾經(jīng)指出,夏季西太副高位置偏南時(shí),東亞地區(qū)熱帶季風(fēng)槽低壓環(huán)流減弱,梅雨鋒低壓環(huán)流加強(qiáng),在850 hPa距平風(fēng)場(chǎng)上30°N以北地區(qū)會(huì)出現(xiàn)氣旋性距平環(huán)流(圖略)。該氣旋性距平環(huán)流與圖5中的結(jié)果較為相近,因此認(rèn)為氣溶膠的影響效果有可能導(dǎo)致夏季東亞地區(qū)西太副高位置的南移,并對(duì)夏季風(fēng)環(huán)流起到一定的削弱作用。

    大氣氣壓場(chǎng)的改變直接導(dǎo)致風(fēng)場(chǎng)的變化。東亞地區(qū),尤其是中緯度地區(qū)所呈現(xiàn)出的西部偏高(110°E以西地區(qū)為正距平),東部偏低(中國(guó)東部及沿海區(qū)域?yàn)樨?fù)距平)態(tài)勢(shì)會(huì)加強(qiáng)我國(guó)江淮以東洋面上氣旋式距平環(huán)流中心南側(cè)的西風(fēng)氣流(圖5),促進(jìn)了該氣旋式環(huán)流的發(fā)展,并使得西太副高位置南移,阻礙了夏季西南季風(fēng)的向北推進(jìn)。同時(shí),氣旋式距平環(huán)流西側(cè)的偏北氣流將進(jìn)一步削弱我國(guó)東部地區(qū)的夏季偏南季風(fēng)。

    圖4 春季(a、c、e)和夏季(b、d、f)氣溶膠對(duì)地表氣壓(a、b;單位:hPa)、850 hPa位勢(shì)高度(c、d;單位:gpm)和500 hPa位勢(shì)高度(e、f;單位:gpm)影響的模擬Fig.4 The simulated (a,c,e)spring and (b,d,f)summer aerosol effects on (a,b)surface pressure(units:hPa),(c,d)850 hPa potential height and (e,f)500 hPa potential height(units:gpm)

    圖5 氣溶膠對(duì)東亞地區(qū)夏季850 hPa風(fēng)場(chǎng)的影響(單位:m/s)Fig.5 Aerosol effect on summer 850 hPa wind field in East Asia(units:m/s)

    2.4 中國(guó)地區(qū)人為氣溶膠排放與東亞夏季風(fēng)的關(guān)系

    為便于分析氣溶膠對(duì)東亞地區(qū)夏季風(fēng)的可能影響,找出前期及同期的預(yù)報(bào)信號(hào),將中國(guó)境內(nèi)氣溶膠排放源較為集中的地區(qū)劃分為三個(gè)子區(qū)域(吉振明等,2010),分別為:中國(guó)北方地區(qū)(109.3~124.3°E,33.4~42.4°N)、中國(guó)南方地區(qū)(109.3~124.3°E,21.5~33.4°N)和四川盆地地區(qū)(105.5~109.3°E,27.7~32.6°N,)。在分別對(duì)1988—2009年春季三個(gè)子區(qū)域氣溶膠柱濃度含量做標(biāo)準(zhǔn)化時(shí)間序列(圖略)后,發(fā)現(xiàn)它們均有比較明顯的年代際變化,在20世紀(jì)80年代末至90年代中期,氣溶膠的柱濃度含量相對(duì)偏低,而從90年代末至今的10 a里有明顯的增加。

    圖6 春季中國(guó)北方地區(qū)氣溶膠柱濃度含量區(qū)域平均值與東亞地區(qū)夏季850 hPa經(jīng)向風(fēng)的相關(guān)分布(a)以及春季中國(guó)北方地區(qū)氣溶膠柱濃度含量高、低值年的東亞夏季850 hPa經(jīng)向風(fēng)合成差值分布(b;單位:m/s;陰影表示通過(guò)0.05信度的顯著性檢驗(yàn))Fig.6 (a)The correlation coefficients between the spring column burden averaged in North China and summer 850 hPa meridional wind in East Asia;(b)The composite difference of East Asia summer 850 hPa meridional wind in years with aerosol column burden at high and low values in North China in spring(units:m/s;the shaded areas denote the significance at 95% confidence level)

    分別計(jì)算了春季三個(gè)子區(qū)域氣溶膠柱濃度含量與敏感性試驗(yàn)中東亞地區(qū)夏季850 hPa經(jīng)向風(fēng)分量v的相關(guān)系數(shù)后看出,中國(guó)中東部地區(qū)均為相關(guān)系數(shù)負(fù)值區(qū),其中中國(guó)北方地區(qū)氣溶膠含量與偏南風(fēng)的相關(guān)性最為顯著,25°N以北的中國(guó)大陸地區(qū)為負(fù)相關(guān)區(qū),相關(guān)系數(shù)高值區(qū)位于25~35°N間的中國(guó)中東部地區(qū),中心值超過(guò)了-0.6,通過(guò)了0.05信度的顯著性檢驗(yàn)(圖6a)。因此,當(dāng)春季我國(guó)北方地區(qū)氣溶膠排放量偏高時(shí),當(dāng)年夏季我國(guó)中東部地區(qū)夏季風(fēng)被削弱;反之,當(dāng)氣溶膠的柱濃度含量偏低時(shí),夏季850 hPa偏南風(fēng)將偏強(qiáng)。在合成差值t檢驗(yàn)中,將中國(guó)北方地區(qū)氣溶膠柱濃度含量區(qū)域平均值大于0.5的年份定義為氣溶膠排放量高值年;將小于-0.5的年份定義為排放量低值年(圖7),據(jù)此得到6個(gè)高值年(1997、1998、2000、2002、2008、2009年)和6個(gè)低值年(1989、1992、1994、1995、2001、2004年),對(duì)相應(yīng)的夏季850 hPa經(jīng)向風(fēng)場(chǎng)進(jìn)行了正、負(fù)位相年的合成分析后,結(jié)果顯示:春季,合成差值t檢驗(yàn)的結(jié)果與相關(guān)系數(shù)的分布形勢(shì)較為一致(圖6b),我國(guó)中緯度地區(qū)夏季經(jīng)向風(fēng)明顯減小,表明在春季我國(guó)北方地區(qū)氣溶膠排放量不斷增大的年際變化背景下,東亞夏季風(fēng)環(huán)流在逐漸削弱。

    夏季,氣溶膠柱濃度含量的年代際變化有著與春季相似的特征,三個(gè)子區(qū)域氣溶膠排放量依然體現(xiàn)出逐漸增長(zhǎng)的趨勢(shì)(圖略)。相關(guān)性研究表明東亞大陸地區(qū)為負(fù)相關(guān)區(qū),中國(guó)南方地區(qū)氣溶膠柱濃度含量與其相關(guān)效果最明顯,如圖8a所示,相關(guān)性分析表明負(fù)相關(guān)系數(shù)高值區(qū)出現(xiàn)在長(zhǎng)江中下游地區(qū),中心值也在-0.6以上,且強(qiáng)度較春季更大,這表明當(dāng)夏季我國(guó)東部地區(qū)850 hPa經(jīng)向風(fēng)偏大時(shí),南方地區(qū)氣溶膠濃度偏低,而當(dāng)夏季中國(guó)南方地區(qū)氣溶膠排放量增加,氣溶膠柱濃度含量升高時(shí),東亞夏季風(fēng)就會(huì)被削弱。夏季中國(guó)南方地區(qū)氣溶膠柱濃度含量(圖9)高、低值年統(tǒng)計(jì)如下:高值年份為1996、1998、1999、2000、2001、2002、2005、2008年,共8 a;低值年份為1988、1989、1990、1991、1992、1994、1995、2006年,同樣為8 a,同期的差值t檢驗(yàn)的結(jié)果更進(jìn)一步證明了上述結(jié)論(圖8b)。

    圖7 1988—2009年春季中國(guó)北方地區(qū)氣溶膠柱濃度含量區(qū)域平均值的逐年變化Fig.7 Interannual variations of spring column burden averaged in North China during 1988—2009

    圖8 夏季中國(guó)南方地區(qū)氣溶膠柱濃度含量區(qū)域平均值與東亞地區(qū)夏季850 hPa經(jīng)向風(fēng)的相關(guān)分布(a)以及夏季中國(guó)南方地區(qū)氣溶膠柱濃度含量高、低值年的東亞地區(qū)夏季850 hPa經(jīng)向風(fēng)合成差值分布(b;單位:m/s;陰影表示通過(guò)0.05信度的顯著性檢驗(yàn))Fig.8 (a)The correlation coefficients between the summer column burden averaged in South China and summer 850 hPa meridional wind in East Asia;(b)The composite difference of East Asia summer 850 hPa meridional wind in years with aerosol column burden at high and low values in South China in summer(units:m/s;the shaded areas denote the significance at 95% confidence level)

    圖9 1988—2009年夏季中國(guó)南方地區(qū)氣溶膠柱濃度含量區(qū)域平均值的逐年變化Fig.9 Interannual variations of summer column burden averaged in South China during 1988—2009

    3 結(jié)論和討論

    1)東亞地區(qū)春季氣溶膠主要分布在中國(guó)南方和四川盆地地區(qū),夏季由于低層以西南氣流的經(jīng)向北上為主,氣溶膠向北輸送,因此在華北地區(qū)出現(xiàn)高值中心,而四川盆地地區(qū)由于地形作用也使得其成為氣溶膠的主要沉積區(qū)。

    2)春、夏季氣溶膠在大氣層頂及地表的短波輻射強(qiáng)迫分布情況與氣溶膠濃度間保持了較好的一致性,綜合考慮了三類(lèi)人為排放氣溶膠之后我國(guó)東部地區(qū)均以負(fù)輻射強(qiáng)迫效應(yīng)為主,夏季效果更為顯著,輻射強(qiáng)迫中心區(qū)域位于我國(guó)四川盆地及江淮流域地區(qū)。模式在引入氣溶膠的影響作用后,在我國(guó)江淮地區(qū)以東洋面上空出現(xiàn)了一個(gè)氣旋式距平環(huán)流中心,該中心位于長(zhǎng)江出??谕?其西側(cè)偏北氣流的影響將會(huì)有效地削弱我國(guó)東部地區(qū)夏季西南季風(fēng)。這可能與氣溶膠減小了海陸氣壓差異,并引起夏季大氣位勢(shì)高度場(chǎng)“西高東低”的變化態(tài)勢(shì)有關(guān),這種變化會(huì)進(jìn)一步加強(qiáng)我國(guó)江淮以東洋面上的氣旋式環(huán)流,迫使西太副高位置南移,削弱東亞夏季風(fēng)。

    3)相關(guān)性分析及合成差值檢驗(yàn)結(jié)果表明,中國(guó)東部地區(qū)夏季850 hPa經(jīng)向風(fēng)與春季中國(guó)北方及夏季中國(guó)南方地區(qū)氣溶膠濃度間相關(guān)顯著,當(dāng)氣溶膠濃度增加時(shí),中國(guó)東部地區(qū)夏季偏南季風(fēng)減弱。

    需要提出的是,由于氣溶膠及其氣候效應(yīng),尤其是氣溶膠的間接效應(yīng)十分復(fù)雜,充滿(mǎn)著不確定性,這加大了對(duì)氣溶膠氣候效應(yīng)模擬的難度。本次試驗(yàn)僅僅考慮了氣溶膠的直接氣候效應(yīng),這無(wú)疑是不全面的。同時(shí),不同氣候模式參數(shù)化方案間的差異以及不同氣溶膠排放源數(shù)據(jù)間的差別也是很大的,因此,在今后的研究過(guò)程中需要進(jìn)一步對(duì)不同參數(shù)化方案下,不同源強(qiáng)數(shù)據(jù)驅(qū)動(dòng)下的模擬結(jié)果進(jìn)行對(duì)比分析,以期對(duì)氣溶膠的氣候效應(yīng)獲得更深刻的認(rèn)識(shí)。

    陳麗,銀燕.2008.礦物氣溶膠遠(yuǎn)程傳輸過(guò)程中的吸收增溫效應(yīng)對(duì)云和降水的影響[J].高原氣象,27(3):628-636.

    高學(xué)杰,徐影,趙宗慈,等.2006.數(shù)值模式不同分辨和地形對(duì)東亞降水模擬影響的試驗(yàn)[J].大氣科學(xué),30(2):185-192.

    韓晉平,王會(huì)軍.2007.東亞夏季風(fēng)年代際變化的若干重要特征及兩份大氣再分析資料的異同[J].地球物理學(xué)報(bào),50(6):1666-1676.

    黃榮輝,陳際龍,周連童,等.2003.關(guān)于中國(guó)重大氣候?yàn)?zāi)害與東亞氣候系統(tǒng)之間關(guān)系的研究[J].大氣科學(xué),27(4):770-787.

    吉振明,高學(xué)杰,張冬峰,等.2010.亞洲地區(qū)氣溶膠及其對(duì)中國(guó)區(qū)域氣候影響的數(shù)值模擬[J].大氣科學(xué),34(2):262-274.

    沈志寶,魏麗.1999.中國(guó)西北大氣沙塵對(duì)地氣系統(tǒng)和大氣輻射加熱的影響[J].高原氣象,18(3):425-435.

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

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

    王喜紅,石廣玉,馬曉燕.2002.東亞地區(qū)對(duì)流層人為硫酸鹽輻射強(qiáng)迫及其溫度響應(yīng)[J].大氣科學(xué),26(6):751-760.

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

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

    鐘昌琴,張美根,朱凌云.2008.東亞地區(qū)春季硫氧化物濃度分布及其收支研究[J].高原氣象,27(1):122-127.

    Chung C E,Ramanathan V,Kiehl J T.2002.Effects of the South Asian absorbing haze on the northeast monsoon and surface-air heat exchange[J].J Climate,15:2462-2476.

    Gao X J,Shi Y,Song R Y,et al.2008.Reduction of future monsoon precipitation over China:Comparison between a high resolution RCM simulation and the driving GCM[J].Meteor Atmos Phys,100:73-86.

    Hansen J,Sato M,Ruedy R,et al.1998.Global warming in the twenty first century:An alternative scenario[C]//Proceeding of the National Academy of Science of the United States of America.97(18):9875-9880.doi:10.1073/pnas.170278997.

    Kiehl J T,Hack J J,Bonan G B,et al.1996.Description of the NCAR community climate model(CCM3)[R]//Tech.Rep.NCAR/TN-420+STR.Boulder:National Center for Atmospheric Research.

    Leon D R,Ulrike L.2002.Tropical rainfall trends and the indirect aerosol effect[J].J Climate,15:2103-2116.

    Massimo B,Sumant N,Lau K M.2008.Absorbing aerosol and summer monsoon evolution over South Asia:An observational portrayal[J].J Climate,21:3221-3239.

    Ohara T,Akimoto H,Kurokawa J,et al.2007.An Asian emission inventory of anthropogenic emission sources for the period 1980—2020[J].Atmos Chem Phys,7:6843-6902.

    Pal J S,Giorgi F,Bi X Q,et al.2007.Regional climate modeling for the developing world:The ICTP RegCM3 and RegCNET[J].Bull Amer Meteor Soc,88(9):1395-1409.

    Qian Y,Giorgi F,Huang Y,et al.2001.Regional simulation of anthropogenic sulfur over East Asia and its sensitivity to model parameters[J].Tellus B,53:171-191.

    Reynolds R W,Rayner N A,Smith T M,et al.2002.An improved in situ and satellite SST analysis for climate[J].J Climate,15:1609-1625.

    Solmon F,Giorgi F,Liousse C.2006.Aerosol modeling for regional climate studies:Application to anthropogenic particles and evaluation over a European/African domain[J].Tellus B,58:51-72.

    Xu Qun.2001.Abrupt change of the mid-summer climate in central east China by the influence of atmospheric pollution[J].Atmos Environ,35:5029-5040.

    Zhou T J,Yu R C.2006.Twentieth century surface air temperature over China and the globe simulation by coupled climate models[J].J Climate,19(22):5843-5858.

    (責(zé)任編輯:張福穎)

    Simulationofanthropogenicaerosols’effectonEastAsiansummermonsooninrecent22years

    CHEN Ming-cheng1,2,SUN Zhao-bo1,SHEN Xin-yong1,2,WANG Yong1,JI Zhen-ming3

    (1.Key Laboratory of Meteorological Disaster of Ministry of Education,Nanjing University of Information Science and Technology,Nanjing 210044,China;2.Key Laboratory of Cloud-Precipitation Physics and Sever Storms,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029,China;3.Institute of Tibetan Plateau Research,Key Laboratory of Tibetan Environment Changes and Land Surface Processes,Chinese Academy of Sciences,Beijing 100085,China)

    The anthropogenic aerosols(sulfate,black carbon and organic carbon) and their effects on the East Asian summer monsoon have been simulated by using a regional climate model RegCM3 coupled with a chemistry aerosol model.Results show that aerosol can change summer 850 hPa wind over East Asia.A cyclonic anomaly circulation center appears on the estuary of the Yangtze River and the northerly air flow on its west will weaken the summer southwestly monsoon in East Asian.Investigation has been conducted on the relationship between spring and summer aerosol column burden over China and summer 850 hPa meridional wind in East Asia.Evidence shows that the aerosol column burden in both spring and summer bears a close negative correlation with summer 850 hPa meridional wind in East Asia.The summer southerly monsoon in East China weakens when the spring aerosol column burden over North China or the summer aerosol column burden over South China increases.The reasons may be the short wave radiative forcing caused by aerosol effects on the top of the atmosphere and the surface,which leads to the change of air pressure and potential height between sea and land.

    anthropogenic aerosol;East Asian summer monsoon;monsoon circulation;numerical simulation;correlation analysis

    2012-01-05;改回日期2012-04-26

    國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)項(xiàng)目(2011CB403405;2013CB430103);國(guó)家自然科學(xué)基金資助項(xiàng)目(41375058;41175065);江蘇高等學(xué)校優(yōu)秀科技創(chuàng)新團(tuán)隊(duì)計(jì)劃項(xiàng)目(PIT2012)

    沈新勇,博士,教授,研究方向?yàn)橹谐叨葰庀髮W(xué)及氣溶膠的氣候效應(yīng),shenxy@nuist.edu.cn.

    10.13878/j.cnki.dqkxxb.20120105005.

    1674-7097(2014)04-0395-10

    P401

    A

    10.13878/j.cnki.dqkxxb.20120105005

    陳明誠(chéng),孫照渤,沈新勇,等.2014.近22 a人為氣溶膠對(duì)東亞地區(qū)夏季風(fēng)環(huán)流影響的數(shù)值模擬研究[J].大氣科學(xué)學(xué)報(bào),37(4):395-404.

    Chen Ming-cheng,Sun Zhao-bo,Shen Xin-yong,et al.2014.Simulation of anthropogenic aerosols’ effect on East Asian summer monsoon in recent 22 years[J].Trans Atmos Sci,37(4):395-404.(in Chinese)

    精品久久久久久久久久免费视频 | 国产精品自产拍在线观看55亚洲 | 国产欧美日韩综合在线一区二区| 日本精品一区二区三区蜜桃| 国产不卡一卡二| 亚洲精品在线观看二区| 亚洲 国产 在线| 午夜精品国产一区二区电影| 超碰97精品在线观看| 在线观看免费视频日本深夜| 国产精品电影一区二区三区 | 久久精品国产清高在天天线| 久久天堂一区二区三区四区| 两个人免费观看高清视频| 80岁老熟妇乱子伦牲交| 日韩欧美一区视频在线观看| 欧美av亚洲av综合av国产av| 免费看a级黄色片| 色在线成人网| 美女午夜性视频免费| 亚洲国产毛片av蜜桃av| 欧美日韩黄片免| 又黄又粗又硬又大视频| 午夜福利影视在线免费观看| 又黄又爽又免费观看的视频| 又黄又爽又免费观看的视频| 国产精品1区2区在线观看. | 人人妻人人爽人人添夜夜欢视频| 大香蕉久久网| 日韩欧美免费精品| 亚洲黑人精品在线| 一级片免费观看大全| 麻豆乱淫一区二区| 黄片播放在线免费| 成人精品一区二区免费| 久久国产亚洲av麻豆专区| 午夜免费鲁丝| 人妻丰满熟妇av一区二区三区 | 久久影院123| av国产精品久久久久影院| 少妇猛男粗大的猛烈进出视频| 免费一级毛片在线播放高清视频 | 大香蕉久久成人网| 亚洲人成电影免费在线| aaaaa片日本免费| www.精华液| 国产日韩一区二区三区精品不卡| 我的亚洲天堂| 免费一级毛片在线播放高清视频 | 这个男人来自地球电影免费观看| 18在线观看网站| 69精品国产乱码久久久| 最新在线观看一区二区三区| 国产精品欧美亚洲77777| 午夜免费鲁丝| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利欧美成人| 欧美激情 高清一区二区三区| 欧美日本中文国产一区发布| 叶爱在线成人免费视频播放| 最新的欧美精品一区二区| 日韩有码中文字幕| 性少妇av在线| www.熟女人妻精品国产| 黄色怎么调成土黄色| 精品高清国产在线一区| av电影中文网址| 欧美国产精品va在线观看不卡| 女同久久另类99精品国产91| 老汉色∧v一级毛片| 新久久久久国产一级毛片| 久久精品熟女亚洲av麻豆精品| 久久精品国产清高在天天线| 黑丝袜美女国产一区| av网站免费在线观看视频| 麻豆国产av国片精品| 国产亚洲欧美精品永久| 久久影院123| 男男h啪啪无遮挡| 色老头精品视频在线观看| 黄色 视频免费看| 午夜福利视频在线观看免费| 宅男免费午夜| 美女高潮到喷水免费观看| 狠狠狠狠99中文字幕| 午夜精品在线福利| 日本wwww免费看| 成人三级做爰电影| 久久影院123| 国产精品久久电影中文字幕 | 又紧又爽又黄一区二区| 国产精品影院久久| 人人妻人人爽人人添夜夜欢视频| 性色av乱码一区二区三区2| 亚洲色图 男人天堂 中文字幕| 女人精品久久久久毛片| 在线播放国产精品三级| 亚洲成人手机| 在线观看66精品国产| 老熟妇仑乱视频hdxx| 免费少妇av软件| 99在线人妻在线中文字幕 | 欧美国产精品va在线观看不卡| 黄色 视频免费看| 精品国产一区二区三区四区第35| 高清视频免费观看一区二区| 视频区欧美日本亚洲| 1024视频免费在线观看| 日本vs欧美在线观看视频| av一本久久久久| netflix在线观看网站| 国产亚洲av高清不卡| 中文字幕人妻丝袜一区二区| 精品一品国产午夜福利视频| 国产三级黄色录像| a在线观看视频网站| 十分钟在线观看高清视频www| 黑人猛操日本美女一级片| 一级毛片高清免费大全| 中文字幕精品免费在线观看视频| 精品一区二区三区视频在线观看免费 | 国产成人免费观看mmmm| 在线播放国产精品三级| 99久久99久久久精品蜜桃| 一区二区三区激情视频| 精品一区二区三卡| 天堂动漫精品| 亚洲av熟女| 免费不卡黄色视频| 免费在线观看日本一区| a级毛片黄视频| 纯流量卡能插随身wifi吗| 国产又色又爽无遮挡免费看| 老司机影院毛片| 99riav亚洲国产免费| 首页视频小说图片口味搜索| 正在播放国产对白刺激| 久久精品亚洲精品国产色婷小说| 欧美精品一区二区免费开放| 国产99久久九九免费精品| 在线免费观看的www视频| a级毛片在线看网站| 成人av一区二区三区在线看| 热99re8久久精品国产| 欧美精品啪啪一区二区三区| 国产精品98久久久久久宅男小说| 波多野结衣av一区二区av| 在线永久观看黄色视频| 久久久国产精品麻豆| 亚洲五月天丁香| 夜夜夜夜夜久久久久| 国产高清国产精品国产三级| 757午夜福利合集在线观看| videos熟女内射| 国产精品.久久久| 欧美亚洲 丝袜 人妻 在线| 国产精品亚洲av一区麻豆| 成人精品一区二区免费| 一级黄色大片毛片| 国产在线一区二区三区精| 国产一卡二卡三卡精品| 18禁美女被吸乳视频| 精品熟女少妇八av免费久了| 校园春色视频在线观看| 一级毛片高清免费大全| 丰满饥渴人妻一区二区三| 午夜福利欧美成人| 久久久久久久午夜电影 | 大型av网站在线播放| 午夜老司机福利片| 黄片播放在线免费| 亚洲色图av天堂| 热re99久久精品国产66热6| 午夜免费成人在线视频| 啦啦啦 在线观看视频| 亚洲五月天丁香| 777米奇影视久久| 真人做人爱边吃奶动态| 一边摸一边做爽爽视频免费| 国产又色又爽无遮挡免费看| 成人亚洲精品一区在线观看| 制服诱惑二区| 亚洲avbb在线观看| 午夜福利影视在线免费观看| 一进一出抽搐gif免费好疼 | 久久久国产欧美日韩av| 丁香欧美五月| 午夜福利在线免费观看网站| 成人黄色视频免费在线看| 免费av中文字幕在线| 好看av亚洲va欧美ⅴa在| 国产精品1区2区在线观看. | 99re在线观看精品视频| 精品久久久久久,| 日本黄色日本黄色录像| 视频在线观看一区二区三区| 中出人妻视频一区二区| 韩国av一区二区三区四区| 999久久久国产精品视频| 午夜福利,免费看| 一级,二级,三级黄色视频| 精品一区二区三区四区五区乱码| 国产高清videossex| 97人妻天天添夜夜摸| 搡老岳熟女国产| 精品国产一区二区三区久久久樱花| 久久香蕉激情| 久久午夜综合久久蜜桃| 又黄又爽又免费观看的视频| 亚洲伊人色综图| 国产一区二区激情短视频| 一区在线观看完整版| 黄色a级毛片大全视频| 国产亚洲av高清不卡| 国产蜜桃级精品一区二区三区 | 在线十欧美十亚洲十日本专区| 精品亚洲成国产av| 国产一区二区三区在线臀色熟女 | 国产精品二区激情视频| а√天堂www在线а√下载 | 国产野战对白在线观看| 女性生殖器流出的白浆| 国产欧美日韩综合在线一区二区| 最新美女视频免费是黄的| 老司机在亚洲福利影院| 搡老岳熟女国产| 欧美日韩亚洲综合一区二区三区_| 国产精品美女特级片免费视频播放器 | 久久久久久亚洲精品国产蜜桃av| 国产精品电影一区二区三区 | 桃红色精品国产亚洲av| 男女高潮啪啪啪动态图| 97人妻天天添夜夜摸| 亚洲片人在线观看| 久久99一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 午夜福利影视在线免费观看| 亚洲伊人色综图| 久久国产精品大桥未久av| 麻豆国产av国片精品| 国产精品电影一区二区三区 | 国产欧美日韩一区二区精品| xxx96com| 久久久久久免费高清国产稀缺| 国产不卡av网站在线观看| 99riav亚洲国产免费| 久久久久精品国产欧美久久久| 国产精品一区二区精品视频观看| 久久精品aⅴ一区二区三区四区| 高清黄色对白视频在线免费看| 十分钟在线观看高清视频www| 精品国产国语对白av| 色在线成人网| 亚洲av电影在线进入| 如日韩欧美国产精品一区二区三区| 精品第一国产精品| 亚洲专区国产一区二区| 国产免费av片在线观看野外av| av超薄肉色丝袜交足视频| 色尼玛亚洲综合影院| 午夜两性在线视频| 黄色女人牲交| 超色免费av| 黄色a级毛片大全视频| 黄色 视频免费看| 大香蕉久久网| 人人妻人人澡人人看| 热re99久久国产66热| 欧美 日韩 精品 国产| 亚洲熟妇中文字幕五十中出 | 18在线观看网站| 国产成人av激情在线播放| 国产一区二区三区视频了| 精品人妻熟女毛片av久久网站| 中文欧美无线码| 国产99久久九九免费精品| 999精品在线视频| 日本黄色日本黄色录像| 欧美日韩一级在线毛片| 黑人操中国人逼视频| 日韩欧美国产一区二区入口| 黄色视频,在线免费观看| 亚洲免费av在线视频| 高清欧美精品videossex| 免费少妇av软件| 最近最新中文字幕大全免费视频| 法律面前人人平等表现在哪些方面| 国产乱人伦免费视频| 久久国产精品男人的天堂亚洲| 男女下面插进去视频免费观看| 亚洲精品国产色婷婷电影| 亚洲精品中文字幕在线视频| 黄色视频不卡| 美女 人体艺术 gogo| 亚洲色图综合在线观看| 男人舔女人的私密视频| 日韩熟女老妇一区二区性免费视频| 很黄的视频免费| 9191精品国产免费久久| 亚洲第一av免费看| 啪啪无遮挡十八禁网站| 黑人欧美特级aaaaaa片| 国产一区二区三区视频了| 久久久水蜜桃国产精品网| 国产极品粉嫩免费观看在线| 国产亚洲精品一区二区www | av福利片在线| 91九色精品人成在线观看| 午夜福利欧美成人| 最新的欧美精品一区二区| 日日摸夜夜添夜夜添小说| 日本五十路高清| 国产伦人伦偷精品视频| 熟女少妇亚洲综合色aaa.| 中文亚洲av片在线观看爽 | 丰满迷人的少妇在线观看| ponron亚洲| 日韩欧美免费精品| 91麻豆av在线| 亚洲三区欧美一区| 亚洲精品粉嫩美女一区| 免费在线观看亚洲国产| 两性午夜刺激爽爽歪歪视频在线观看 | 777米奇影视久久| 一边摸一边做爽爽视频免费| 在线天堂中文资源库| 无人区码免费观看不卡| 亚洲午夜理论影院| 一a级毛片在线观看| 亚洲熟女精品中文字幕| 色播在线永久视频| 日韩 欧美 亚洲 中文字幕| 这个男人来自地球电影免费观看| xxxhd国产人妻xxx| 无遮挡黄片免费观看| 久久天躁狠狠躁夜夜2o2o| 最近最新中文字幕大全电影3 | 亚洲中文字幕日韩| 19禁男女啪啪无遮挡网站| 91精品国产国语对白视频| 熟女少妇亚洲综合色aaa.| 动漫黄色视频在线观看| 午夜免费成人在线视频| 啪啪无遮挡十八禁网站| 久久久久久久久免费视频了| www.自偷自拍.com| 久久99一区二区三区| 国产激情欧美一区二区| 久久久国产欧美日韩av| 黄片大片在线免费观看| 在线国产一区二区在线| 亚洲片人在线观看| 欧美色视频一区免费| 在线观看舔阴道视频| 另类亚洲欧美激情| bbb黄色大片| 老司机亚洲免费影院| 天堂√8在线中文| 国产97色在线日韩免费| 九色亚洲精品在线播放| 亚洲成av片中文字幕在线观看| 老司机午夜十八禁免费视频| 757午夜福利合集在线观看| 高清视频免费观看一区二区| 老司机午夜十八禁免费视频| a级毛片黄视频| 欧美日本中文国产一区发布| 亚洲第一欧美日韩一区二区三区| 动漫黄色视频在线观看| 国产av精品麻豆| 熟女少妇亚洲综合色aaa.| 老汉色∧v一级毛片| 香蕉丝袜av| 多毛熟女@视频| 999久久久国产精品视频| 91av网站免费观看| 日韩欧美在线二视频 | 一区在线观看完整版| av不卡在线播放| 天堂动漫精品| 777米奇影视久久| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品中文字幕一二三四区| av天堂久久9| 新久久久久国产一级毛片| 欧美精品人与动牲交sv欧美| 99精品在免费线老司机午夜| 女人高潮潮喷娇喘18禁视频| 国产精品 国内视频| 在线观看免费高清a一片| 日韩 欧美 亚洲 中文字幕| 男人操女人黄网站| 欧美人与性动交α欧美软件| a级片在线免费高清观看视频| 国产欧美日韩一区二区精品| 又黄又爽又免费观看的视频| 一二三四社区在线视频社区8| 在线av久久热| 久久国产精品影院| 一个人免费在线观看的高清视频| 校园春色视频在线观看| 12—13女人毛片做爰片一| 亚洲国产精品一区二区三区在线| 精品少妇久久久久久888优播| 久久国产乱子伦精品免费另类| 青草久久国产| 两性午夜刺激爽爽歪歪视频在线观看 | 视频在线观看一区二区三区| 人人澡人人妻人| 亚洲国产精品sss在线观看 | 亚洲av电影在线进入| 国产高清国产精品国产三级| 女人高潮潮喷娇喘18禁视频| av超薄肉色丝袜交足视频| 亚洲av片天天在线观看| 黄色女人牲交| 老司机午夜福利在线观看视频| 欧美日韩黄片免| 天天添夜夜摸| 丰满迷人的少妇在线观看| 波多野结衣av一区二区av| 欧美亚洲 丝袜 人妻 在线| 一本综合久久免费| 亚洲av美国av| 成人三级做爰电影| 一进一出抽搐gif免费好疼 | 18在线观看网站| 亚洲专区字幕在线| 亚洲精品美女久久久久99蜜臀| 人妻久久中文字幕网| 国产精品综合久久久久久久免费 | 最新美女视频免费是黄的| 新久久久久国产一级毛片| 亚洲熟妇中文字幕五十中出 | www.熟女人妻精品国产| 丝瓜视频免费看黄片| 十八禁人妻一区二区| 欧美日韩瑟瑟在线播放| 成人亚洲精品一区在线观看| 久久午夜综合久久蜜桃| 午夜影院日韩av| av福利片在线| 女同久久另类99精品国产91| 午夜亚洲福利在线播放| 精品高清国产在线一区| 国产精品国产av在线观看| 亚洲免费av在线视频| ponron亚洲| 这个男人来自地球电影免费观看| 欧美 亚洲 国产 日韩一| 老司机午夜十八禁免费视频| 麻豆乱淫一区二区| 亚洲黑人精品在线| 99久久99久久久精品蜜桃| 美女高潮喷水抽搐中文字幕| 国产成人啪精品午夜网站| 国产亚洲精品第一综合不卡| 午夜精品久久久久久毛片777| 亚洲午夜精品一区,二区,三区| 俄罗斯特黄特色一大片| 黄色毛片三级朝国网站| 丰满的人妻完整版| 婷婷丁香在线五月| 午夜视频精品福利| 超色免费av| 欧美不卡视频在线免费观看 | 国产精品综合久久久久久久免费 | 亚洲国产精品合色在线| 露出奶头的视频| 男男h啪啪无遮挡| 在线天堂中文资源库| 男女高潮啪啪啪动态图| 国产成人系列免费观看| 无人区码免费观看不卡| ponron亚洲| 麻豆国产av国片精品| 一二三四在线观看免费中文在| 精品国产一区二区三区久久久樱花| 日日爽夜夜爽网站| 黄片小视频在线播放| 高清av免费在线| 无遮挡黄片免费观看| 国产精品亚洲av一区麻豆| 老熟妇仑乱视频hdxx| 黄片大片在线免费观看| 水蜜桃什么品种好| 电影成人av| 岛国在线观看网站| 18禁美女被吸乳视频| 国产成人免费无遮挡视频| 高清视频免费观看一区二区| 国产精品 欧美亚洲| 国产又爽黄色视频| 久久这里只有精品19| 丝袜在线中文字幕| 视频在线观看一区二区三区| 欧美国产精品va在线观看不卡| 视频区欧美日本亚洲| 亚洲色图综合在线观看| 午夜福利,免费看| 中文字幕最新亚洲高清| 久久天堂一区二区三区四区| 欧美中文综合在线视频| 一a级毛片在线观看| 久久久国产成人精品二区 | 国产不卡一卡二| 91麻豆av在线| 欧美精品高潮呻吟av久久| 国产97色在线日韩免费| 人人妻,人人澡人人爽秒播| 99热国产这里只有精品6| 欧美日韩亚洲高清精品| 69精品国产乱码久久久| 久久香蕉精品热| 嫁个100分男人电影在线观看| 久久国产精品影院| 超色免费av| 国产视频一区二区在线看| 电影成人av| 亚洲国产精品合色在线| 男女之事视频高清在线观看| netflix在线观看网站| 国产成人精品在线电影| 久久青草综合色| 日本黄色日本黄色录像| 午夜免费鲁丝| 亚洲熟妇中文字幕五十中出 | 看片在线看免费视频| 一二三四社区在线视频社区8| 久久热在线av| 亚洲成av片中文字幕在线观看| 看免费av毛片| 99在线人妻在线中文字幕 | 岛国毛片在线播放| 一级毛片高清免费大全| 岛国毛片在线播放| 99国产精品99久久久久| 午夜福利免费观看在线| av国产精品久久久久影院| 精品人妻熟女毛片av久久网站| 丰满迷人的少妇在线观看| 亚洲成人国产一区在线观看| 国产日韩一区二区三区精品不卡| 亚洲男人天堂网一区| 怎么达到女性高潮| 成人免费观看视频高清| 99久久99久久久精品蜜桃| 精品国产一区二区三区久久久樱花| 亚洲专区字幕在线| 大香蕉久久网| 国产在线精品亚洲第一网站| 一二三四社区在线视频社区8| 黄片播放在线免费| 99re6热这里在线精品视频| 淫妇啪啪啪对白视频| 亚洲av成人av| 成人精品一区二区免费| 99国产综合亚洲精品| 热99国产精品久久久久久7| 亚洲国产精品一区二区三区在线| 国产精品久久久久成人av| 欧美精品人与动牲交sv欧美| 97人妻天天添夜夜摸| 国产精品成人在线| 久久久久久人人人人人| 免费在线观看亚洲国产| 欧美精品高潮呻吟av久久| 亚洲avbb在线观看| 久热这里只有精品99| 欧美大码av| 黄色视频不卡| 91九色精品人成在线观看| 51午夜福利影视在线观看| 国产精品一区二区精品视频观看| 男人的好看免费观看在线视频 | 脱女人内裤的视频| 91成年电影在线观看| 国产亚洲精品一区二区www | 欧美精品啪啪一区二区三区| 亚洲av欧美aⅴ国产| 亚洲精品成人av观看孕妇| 成人手机av| 精品电影一区二区在线| e午夜精品久久久久久久| 一级a爱片免费观看的视频| 丰满人妻熟妇乱又伦精品不卡| 一本一本久久a久久精品综合妖精| 国产成人免费无遮挡视频| 叶爱在线成人免费视频播放| www日本在线高清视频| 亚洲精华国产精华精| 久久国产精品影院| 一区在线观看完整版| 亚洲少妇的诱惑av| 亚洲精品国产一区二区精华液| 国产精品免费一区二区三区在线 | 好男人电影高清在线观看| 18禁裸乳无遮挡动漫免费视频| 久热这里只有精品99| 免费观看精品视频网站| 老司机影院毛片| a级毛片黄视频| 亚洲黑人精品在线| 老汉色∧v一级毛片| 两人在一起打扑克的视频| 国产精品久久久av美女十八| 欧美大码av| 久久久国产精品麻豆| 国产精品国产高清国产av | 亚洲成人免费电影在线观看| av福利片在线| 亚洲欧美色中文字幕在线| 亚洲欧洲精品一区二区精品久久久| 国产欧美日韩精品亚洲av| 精品久久久精品久久久|