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

    基于WRF模式的青藏高原斜坡和平臺加熱影響亞洲夏季風(fēng)的模擬研究

    2016-11-16 00:56:27王子謙段安民李茂善何編
    地球物理學(xué)報 2016年9期
    關(guān)鍵詞:喜馬拉雅山脈季風(fēng)南亞

    王子謙, 段安民, 李茂善, 何編

    1 中山大學(xué)大氣科學(xué)學(xué)院, 廣州 510275 2 中國科學(xué)院大氣物理研究所大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國家重點(diǎn)實(shí)驗(yàn)室(LASG), 北京 100029 3 中國科學(xué)院寒旱區(qū)陸面過程與氣候變化重點(diǎn)實(shí)驗(yàn)室, 蘭州 730000

    ?

    基于WRF模式的青藏高原斜坡和平臺加熱影響亞洲夏季風(fēng)的模擬研究

    王子謙1, 2, 段安民2, 李茂善3, 何編2

    1 中山大學(xué)大氣科學(xué)學(xué)院, 廣州 510275 2 中國科學(xué)院大氣物理研究所大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國家重點(diǎn)實(shí)驗(yàn)室(LASG), 北京 100029 3 中國科學(xué)院寒旱區(qū)陸面過程與氣候變化重點(diǎn)實(shí)驗(yàn)室, 蘭州 730000

    青藏高原大地形的熱力強(qiáng)迫作用對亞洲夏季風(fēng)的形成和發(fā)展具有重要的影響.本文利用較高分辨率的WRF區(qū)域模式,探討了高原不同區(qū)域(斜坡和平臺)的地形加熱分別對南亞夏季風(fēng)和東亞夏季風(fēng)的影響.結(jié)果表明:高原南部喜馬拉雅山脈的斜坡地形加熱對其周圍局地的環(huán)流形勢和降水影響十分明顯,是南亞夏季風(fēng)北支分量形成和維持的主導(dǎo)因子,也是斜坡上氣流爬坡和降水發(fā)生的必要條件.斜坡加熱對東亞夏季風(fēng)也有明顯的增強(qiáng)作用,它不僅加強(qiáng)了中國東部低空西南季風(fēng)環(huán)流,還會造成北部南下的異常干冷空氣的響應(yīng).斜坡上的地形加熱作用也是對流層高層暖中心位置維持在斜坡上空的一個重要原因.而高原平臺加熱對季風(fēng)環(huán)流和降水的影響雖然沒有喜馬拉雅山脈斜坡加熱那么顯著,但是對南亞夏季風(fēng)的影響范圍更廣,對經(jīng)向哈得來環(huán)流影響更明顯,能夠調(diào)控高原以外更遠(yuǎn)處熱帶洋面上的西南季風(fēng)環(huán)流.通過比較高原不同區(qū)域地形加熱條件下的多種季風(fēng)指數(shù),進(jìn)一步表明了高原地形加熱對南亞和東亞夏季風(fēng)均有增強(qiáng)作用,但是高原不同區(qū)域的地形加熱對兩類夏季風(fēng)子系統(tǒng)又會產(chǎn)生不一樣的影響.

    青藏高原; 地形加熱; 亞洲夏季風(fēng); WRF模式

    1 引言

    亞洲夏季風(fēng)系統(tǒng)主要包括南亞夏季風(fēng)和東亞夏季風(fēng)兩個子系統(tǒng),是全球最為典型的季風(fēng)系統(tǒng),其異?;顒又苯佑绊懼蟻喖皷|亞地區(qū)各國的天氣氣候變化.在全球變暖的大背景下,亞洲夏季風(fēng)的演變規(guī)律變得更為復(fù)雜(Webster, 2006),季風(fēng)區(qū)的極端天氣氣候事件也逐漸增多(Herring et al., 2014).夏季風(fēng)異常造成的洪澇、干旱、高溫等災(zāi)害常給人類帶來重大經(jīng)濟(jì)損失,因此其形成、維持以及變異的物理機(jī)制一直以來備受科學(xué)家的關(guān)注.

    以往大量研究表明亞洲夏季風(fēng)所對應(yīng)的降水和環(huán)流變化與青藏高原(簡稱高原)的動力和熱力強(qiáng)迫作用有著密切的聯(lián)系(Flohn, 1957; 葉篤正等, 1979; Wu and Zhang, 1998; Zhao and Chen, 2001; 李躍清, 2003; 梁瀟云等, 2006; 于樂江和胡敦欣, 2008; Xu et al., 2010; Wang et al., 2014b).春夏季,由于高原的地表感熱加熱較強(qiáng),使得其熱力作用對亞洲夏季風(fēng)的影響尤為重要(吳國雄等, 1997; 毛江玉和吳國雄, 2006; Wu et al., 2007; 何金海等, 2011).然而,近年來也有部分研究工作(Boos and Kuang, 2010, 2013; Molnar et al., 2010)對高原的熱力強(qiáng)迫作用提出了質(zhì)疑,強(qiáng)調(diào)影響南亞夏季風(fēng)的主導(dǎo)因子并非是高原的熱力作用.他們利用美國的CAM(Community Atmosphere Model)及CESM(Community Earth System Model)全球模式將高原的平臺部分去除而只保留狹長的喜馬拉雅山脈,模擬得到的南亞夏季風(fēng)和完整地形條件下的結(jié)果彼此相似,由此認(rèn)為南亞夏季風(fēng)形成和維持的主導(dǎo)因子是喜馬拉雅山脈的機(jī)械隔斷作用而不是高原的熱力作用.Wu等(2012)利用中國科學(xué)院大氣物理研究所的SAMIL(Spectral Atmospheric Model of IAP/LASG)全球大氣模式重新設(shè)計試驗(yàn),所得結(jié)果仍然認(rèn)為高原熱力作用對南亞夏季風(fēng)的形成更為重要.該研究工作指出Boos和Kuang(2010)的試驗(yàn)設(shè)計中喜馬拉雅山斜坡上的感熱加熱依然存在,意味著保留了高原地形的感熱氣泵作用,因而不能由此否認(rèn)高原熱力作用對南亞夏季風(fēng)影響的主導(dǎo)地位.Wu等(2012)的模擬結(jié)果表明如果將高原的平臺部分去除而只保留其南側(cè)喜馬拉雅山脈,同時也去除喜馬拉雅山脈地表的感熱加熱,那么南亞夏季風(fēng)的北支分量就無法形成.

    從以上研究背景可以看出,在南亞夏季風(fēng)的形成和維持上,高原的熱力和動力作用哪個是主導(dǎo)因子仍然存在爭議.同時,前人研究結(jié)果也表明只保留喜馬拉雅山脈就能得到全地形下的南亞夏季風(fēng)形態(tài),這說明高原不同區(qū)域的熱力強(qiáng)迫會對南亞夏季風(fēng)或者其他亞洲夏季風(fēng)子系統(tǒng)造成不一樣的影響.高原南部斜坡上的地形加熱存在顯著影響的同時,喜馬拉雅山脈以北的平臺區(qū)域的地形加熱作用又會產(chǎn)生怎樣的影響?哪個區(qū)域的熱力強(qiáng)迫更重要一些?這些問題在以往研究中認(rèn)識的還十分有限.而且除了對南亞夏季風(fēng)的影響,不同區(qū)域的高原地形加熱對東亞夏季風(fēng)的影響也需進(jìn)一步探討.此外,值得提出的是,前人的數(shù)值模擬結(jié)果都是基于分辨率較粗(180~200 km)的全球氣候系統(tǒng)模式來分離喜馬拉雅山和高原平臺部分并研究其分別對南亞夏季風(fēng)的影響,這其中難免會出現(xiàn)由于模型過于粗糙而對模擬結(jié)果帶來不確定性,尤其是針對地形十分復(fù)雜的高原地區(qū).因此,本文將利用較高分辨率的WRF(Weather Research and Forecasting)區(qū)域模式來進(jìn)一步認(rèn)識高原大地形對亞洲夏季風(fēng)的影響,研究高原不同區(qū)域的地形加熱對其周邊季風(fēng)區(qū)大氣環(huán)流和降水的影響差異,并理解其中的物理過程.

    2 模式和試驗(yàn)設(shè)計

    相比于以往相關(guān)工作所用的全球模式,區(qū)域模式具有更高的分辨率以及更為合理的中小尺度動力和物理過程,是研究具有復(fù)雜地形的青藏高原以及復(fù)雜海陸分布的亞洲季風(fēng)區(qū)的有力工具.隨著WRF模式的不斷發(fā)展和完善,近年來越來越多的研究者將它作為區(qū)域氣候模式用于亞洲季風(fēng)的模擬(Kim and Hong, 2010; Yu et al., 2011; Yang et al., 2012; Wang et al., 2014a).本文選用WRF(Version 3.4.1, Skamarock et al., 2008)模式作為研究工具,所用到的具體物理參數(shù)化方案包括:WSM6微物理過程,Grell-Devenyi積云對流方案,Noah陸面模式,BouLac行星邊界層方案,Goddard短波輻射方案,以及RRTM長波輻射方案.模擬范圍覆蓋了大部分亞洲區(qū)域及鄰近海洋區(qū)域(圖1),試驗(yàn)采用蘭伯特地圖投影,中心點(diǎn)位于(30°N,95°E),模式水平分辨率為45 km,南北和東西方向格點(diǎn)數(shù)分別為133和195(模式邊界緩沖區(qū)為5個格點(diǎn)),垂直方向分35層,最頂層為50 hPa.模式初始場和邊界場均來源于NCEP提供的每日四次(世界時00時,06時,12時,18時)及水平分辨率為1°×1°的全球最終分析場(NCEP-FNL),積分過程中邊界條件每六小時更新一次.采用NOAA提供的高分辨率逐日最優(yōu)插值海表溫度(OISST)作為模式的海表強(qiáng)迫場,該資料水平分辨率為0.25°×0.25°.在對模擬結(jié)果驗(yàn)證時還用到了NASA提供的TRMM-3B42降水?dāng)?shù)據(jù),其水平分辨率為0.25°×0.25°.

    前人相關(guān)的研究在分離喜馬拉雅山脈試驗(yàn)設(shè)計中,將高原區(qū)域同一經(jīng)度的地形高度達(dá)到該經(jīng)度上最大值的2/3的北側(cè)全部設(shè)為0,剩下的地形則為喜馬拉雅山脈(Boos and Kuang, 2010; Wu et al., 2012);或者直接以4 km地形高度為界,4 km分界線以南部分為喜馬拉雅山斜坡,而4 km以北為高原平臺部分(何編, 2012).在分辨率較高的WRF模式地形條件下,無論是每個經(jīng)度上海拔最大值的2/3還是直接4 km高度都要遠(yuǎn)低于喜馬拉雅山脈的最大高度,所以在本文的試驗(yàn)設(shè)計中有必要重新定義高原地區(qū)喜馬拉雅山脈及其以北平臺部分的劃分,以免低估高原南部斜坡的加熱作用.如圖2a所示,將每個經(jīng)度上地形高度最大的點(diǎn)作為邊界(如75°E和85°E),如果最大高度所在的格點(diǎn)位于高原偏中部地區(qū)(即喜馬拉雅山脈以北),那么就取該格點(diǎn)以南的第二最大高度為邊界(如92°E).利用該方法得到如圖2b所示的喜馬拉雅山脈斜坡地區(qū)和高原平臺地區(qū),文章將以此來研究高原不同區(qū)域的地形加熱對亞洲夏季風(fēng)的影響.

    針對高原不同區(qū)域的地形加熱設(shè)計了四組模擬試驗(yàn)(表1),分別是參考試驗(yàn)(CTL)、高原南部喜馬拉雅山脈斜坡無地形加熱試驗(yàn)(HIM_NS)、高原平臺無地形加熱試驗(yàn)(PL_NS)以及整個高原無地形加熱試驗(yàn)(TP_NS).修改地形加熱的具體方案是保留地形本身不變的條件下,在模式的每個積分步上將所選區(qū)域格點(diǎn)上所釋放到大氣中的感熱加熱減為0,同時保持地表能量平衡不被破壞(也就是在大氣熱力學(xué)方程中地表溫度的垂直耗散項(xiàng)為0,模式積分過程中陸面模塊仍然計算感熱通量并允許其影響地表溫度).每組試驗(yàn)共包括從2003—2008年6個夏季,模式采用分段積分方式,從每年的5月1日00時(本文所用時間均為世界時)開始積分,到當(dāng)年8月31日18時模擬結(jié)束,主要分析后三個月(夏季JJA)的模擬結(jié)果.圖3a—3d給出了各組試驗(yàn)中高原地區(qū)夏季的地表感熱分布情況,由參考試驗(yàn)CTL可知,夏季高原的地形抬升加熱大值區(qū)主要位于高原的西北部地區(qū),超過100 W·m-2,而高原東南部由于季風(fēng)降水的影響導(dǎo)致地表感熱加熱較小,最小值小于40 W·m-2.

    表1 試驗(yàn)方案設(shè)計Table 1 Detail of experimental design

    圖1 WRF模式模擬區(qū)域和地形高度(陰影,單位:m)圖中虛線方框?yàn)楦咴舾行栽囼?yàn)區(qū).Fig.1 WRF model domain and terrain height (shading, unit: m)Dashed box is the area for sensitivity experiment over the Tibetan Plateau (TP).

    圖2 (a) WRF模式中不同經(jīng)度上地形高度隨緯度的變化(箭頭標(biāo)注的是喜馬拉雅山脈北邊界); (b) 高原地區(qū)(23°N—40°N, 70°E—105°E,對應(yīng)圖1中的虛線框)的地形高度(單位:m),藍(lán)色三角形標(biāo)注的是每個經(jīng)度上地形的最高點(diǎn),藍(lán)線為喜馬拉雅山脈和高原平臺的分界線, 分界線以南并且地形大于300 m的區(qū)域?yàn)橄柴R拉雅山脈斜坡區(qū)域(黑點(diǎn)區(qū)),分界線以北且地形大于1500 m的區(qū)域?yàn)楦咴脚_區(qū)域(白點(diǎn)區(qū))Fig.2 (a) Topography height at different longitude. Arrows denote the north boundary of the Himalayas (HIM). (b) Terrain height (unit: m) over the sensitivity experiment area (23°N—40°N, 70°E—105°E). Blue triangles denote the highest point along every longitude, and the blue curve presents the boundary between HIM and TP′s platform. The area over the south of blue curve and meanwhile topography is above 300 m denotes southern slope of HIM (black dotted region), and the area over the north of blue curve and meanwhile topography is above 1500 m denotes the TP′s platform (white dotted region)

    圖3 各組試驗(yàn)高原地區(qū)夏季的感熱加熱分布及其差異(單位:W·m-2),(a) CTL, (b) HIM_NS, (c) PL_NS, (d) TP_NS, (e) CTL減HIM_NS, (f) CTL減PL_NS, (g) CTL減TP_NS其中(e—g)中打點(diǎn)區(qū)表示差值通過95%信度檢驗(yàn).Fig.3 Summer mean surface sensible heating (unit: W·m-2) for experiments of (a) CTL, (b) HIM_NS, (c) PL_NS, and (d) TP_NS; (e) CTL minus HIM_NS; (f) CTL minus PL_NS; (g) CTL minus TP_NSDotted regions in (e—g) denote the statistical significance of differences above 95% level.

    圖4 夏季降水(陰影,單位:mm·day-1)和850 hPa風(fēng)場(矢量,單位:m·s-1),(a) CTL, (b) HIM_NS, (c) PL_NS, (d) TP_NS, (e) TRMM降水及NCEP-FNL風(fēng)場,(f) CTL減HIM_NS,(g) CTL減PL_NS,(h) CTL減TP_NS其中(f—h)中紅色矢量箭頭和打點(diǎn)區(qū)分別表示風(fēng)場和降水差異通過95%信度檢驗(yàn).Fig.4 Summer mean precipitation (shading, unit: mm·day-1) and wind field at 850 hPa (vectors, unit: m·s-1) for experiments of (a) CTL, (b) HIM_NS, (c) PL_NS, (d) TP_NS, and (e) OBS (wind is from NCEP-FNL, precipitation is from TRMM); (f) CTL minus HIM_NS; (g) CTL minus PL_NS; (h) CTL minus TP_NSDotted regions denote the statistical significance of precipitation anomalies above 95% level, and wind differences significant above 95% level are plotted in red color in (f—h).

    3 模擬結(jié)果分析

    3.1 高低層季風(fēng)環(huán)流、降水及溫度場的響應(yīng)

    圖4給出觀測及各組試驗(yàn)?zāi)M的亞洲夏季風(fēng)降水和850 hPa低空環(huán)流場,由參考試驗(yàn)(圖4a)和觀測數(shù)據(jù)(圖4e)的結(jié)果可知,相比于粗網(wǎng)格的全球氣候系統(tǒng)模式(Boos and Kuang, 2010; Wu et al., 2012),WRF模式能更好地模擬出阿拉伯海、孟加拉灣、南海以及東亞的夏季風(fēng)降水和環(huán)流特征.尤其在青藏高原區(qū)域,模式能再現(xiàn)高原南側(cè)的降水大值中心以及高原主體由東南向西北遞減的降水分布特征,這些結(jié)果在一定程度上可提高之后的敏感性試驗(yàn)的可靠性.從其他敏感性試驗(yàn)的模擬結(jié)果及其相對于參考試驗(yàn)CTL的差異場可以看出,高原不同區(qū)域的地形加熱作用對亞洲夏季風(fēng)具有不一樣的影響.參考試驗(yàn)CTL與無南坡地形加熱試驗(yàn)HIM_NS的差異場表明,喜馬拉雅山脈斜坡地形加熱對其周圍附近的環(huán)流形勢影響十分顯著,它主要加強(qiáng)了繞高原西側(cè)南下的干冷氣流以及高原以南地區(qū)暖濕氣流向北的輻合.從高緯南下的異常干冷空氣同時還增大了印度西北部及巴基斯坦地區(qū)的地表感熱加熱(圖3e),導(dǎo)致該區(qū)域繞高原西南側(cè)的氣旋式環(huán)流進(jìn)一步加強(qiáng),而異常干冷空氣的加強(qiáng)也造成印度西北部及巴基斯坦附近的降水明顯減少(圖4f).此外,圖4f同時也表明喜馬拉雅山脈斜坡加熱有利于高原南部暖濕氣流的爬坡輻合,進(jìn)而使得夏季高原南部的降水驟增,如果沒有喜馬拉雅山脈斜坡的加熱作用,高原南部的降水基本上消失(圖4b).至于東亞夏季風(fēng)的響應(yīng),喜馬拉雅山脈斜坡加熱同樣加強(qiáng)了中國東部低空西南季風(fēng)環(huán)流,以及北部南下的異常干冷空氣,導(dǎo)致異常的冷暖氣流交匯,進(jìn)而使得中國東部的夏季風(fēng)降水增多(圖4f).北部南下的異常干冷空氣主要是由高原加熱在對流層中高層激發(fā)Rossby波列所引起(Wang et al., 2014b),位于高原以東即中國北部上空的氣旋式異常環(huán)流具有相當(dāng)正壓結(jié)構(gòu)(圖6e).氣旋式異常環(huán)流西側(cè)盛行西北氣流,相應(yīng)的在低層華北地區(qū)出現(xiàn)南下的異常干冷空氣.

    全球氣候模式的結(jié)果顯示高原平臺的熱力作用對850 hPa風(fēng)場影響不明顯,對南亞夏季風(fēng)環(huán)流的調(diào)控較弱(何編, 2012).從本文PL_NS試驗(yàn)結(jié)果來看,我們在該試驗(yàn)中將高原平臺的感熱加熱去除,其造成的南亞季風(fēng)環(huán)流和降水響應(yīng)強(qiáng)度的確沒有喜馬拉雅山脈斜坡加熱被抑制(HIM_NS)的試驗(yàn)結(jié)果那么顯著.但是,平臺加熱對南亞夏季風(fēng)的影響范圍更大,圖4g可以看出高原以外的阿拉伯海及孟加拉灣海洋上的西南氣流有明顯的加強(qiáng),這樣造成更多的水汽輸送使得孟加拉灣北部、印度大陸北部以及印度西岸的降水都有所增加.同時,高原平臺加熱的去除直接導(dǎo)致平臺上局地的降水消失.與喜馬拉雅山斜坡加熱類似,平臺加熱同樣有利于東亞地區(qū)西南季風(fēng)環(huán)流的北上以及中國東部季風(fēng)降水的增強(qiáng);但是對北方南下的異常干冷氣流影響不顯著,因而造成中國東部的降水差異沒有HIM_NS試驗(yàn)結(jié)果明顯.此外,當(dāng)整個高原的地形加熱去除后(圖4h),其結(jié)果基本上是斜坡和平臺的綜合效應(yīng);大地形加熱作用顯著增強(qiáng)了高原西側(cè)繞流南下的干冷空氣、熱帶海洋及南亞西南季風(fēng)氣流向高原南部輻合、以及中國東部北上的東亞夏季風(fēng)環(huán)流.這樣的環(huán)流形勢促使印度西北部以及巴基斯坦地區(qū)附近感熱加熱明顯增強(qiáng)(圖3g)以及降水減少,而高原及中國東部降水顯著增多.

    基于以上分析可知高原的地形抬升加熱是造成高原降水的根本原因,一旦加熱被抑制,高原斜坡和平臺上的降水就基本消失.當(dāng)然,不難發(fā)現(xiàn)地形加熱抑制后高原斜坡上仍有少量的降水發(fā)生(圖4b和4d),我們已有的研究表明這主要是由微物理相變過程造成的(Wang et al., 2015).圖5進(jìn)一步給出了高原地區(qū)及周邊的近地表(σ=0.98)環(huán)流場.同時,基于錢永甫等(1988)及張耀存和錢永甫(1999)在進(jìn)行高原動力作用數(shù)值試驗(yàn)中對模式近地面風(fēng)場所采取的計算方法,將地表環(huán)流分解為繞流分量和爬流分量.由圖5b可知,氣候平均狀態(tài)下高原近地表為一個氣旋式環(huán)流,高原南部喜馬拉雅山斜坡以及高原東側(cè)斜坡均存在明顯的爬流.尤其是高原南坡,來自熱帶海洋的暖濕氣流通過爬坡抬升,并在高原上空輻合進(jìn)而產(chǎn)生降水.當(dāng)喜馬拉雅山脈的感熱加熱去除后,高原南坡的爬流強(qiáng)度很快減弱(圖5d),同時從孟加拉灣北上的暖濕氣流遇到高原地形后基本都轉(zhuǎn)變?yōu)槔@流,使得繞流強(qiáng)度明顯較CTL模擬結(jié)果更強(qiáng)(圖5a和5c).因此,沒有喜馬拉雅山脈的斜坡加熱,熱帶的暖濕氣流無法爬坡北上產(chǎn)生抬升運(yùn)動以及激發(fā)對流,進(jìn)而抑制高原南坡的降水發(fā)生.如果只去除平臺加熱而保留喜馬拉雅山脈的加熱,高原地表環(huán)流形勢轉(zhuǎn)變?yōu)橐粋€反氣旋式環(huán)流(圖5f),說明夏季高原近地層的氣旋式環(huán)流場主要受平臺地表加熱所控制.然而,高原南部的爬流依然存在,進(jìn)而使得斜坡上的降水得以維持,但是爬流的強(qiáng)度明顯減弱;同時,平臺無加熱直接造成高原東側(cè)的爬流特征消失,切斷了平臺上東側(cè)的水汽輸送.如果整個高原的加熱均去除,高原近地表完全被反氣旋環(huán)流控制,并在邊界處向外輻散下滑,整個高原周邊的爬流基本消失,水汽無法向高原輸送,而南部的繞流也類似于HIM_NS試驗(yàn)結(jié)果將得到加強(qiáng)(圖5g和5h).

    由于具有高海拔特征,高原的地形加熱效應(yīng)可以直接或間接影響到對流層中上層.所以,高原地形的熱力強(qiáng)迫作用除了影響低層的季風(fēng)環(huán)流和降水,同時也必然影響對流層高層的環(huán)流形勢和溫度場結(jié)構(gòu).根據(jù)參考試驗(yàn)CTL的模擬結(jié)果,夏季高原上空表現(xiàn)為強(qiáng)大的反氣旋式環(huán)流(南亞高壓),其中心位于高原南坡上空,溫度場上對應(yīng)一個顯著的暖中心(圖6a).當(dāng)喜馬拉雅山脈的感熱加熱去除后,暖中心明顯減弱并往南移至印度大陸北部上空(圖6b),可見斜坡上的地形加熱對高層的增暖起到了重要的貢獻(xiàn)作用,并且使得暖中心的位置維持在斜坡上空.從CTL和HIM_NS的差異場(圖6e)也可看出,喜馬拉雅山脈加熱主要位于高原南坡上空,并對應(yīng)一個反氣旋環(huán)流異常,即對南亞高壓有增強(qiáng)作用,同時該加熱中心在其下游東北方向激發(fā)出一個異常的氣旋式環(huán)流,這和以往研究結(jié)果類似,即高原的熱力強(qiáng)迫作用可激發(fā)大尺度波動并向下游傳播(Wang B et al., 2008a; Wang Z Q et al., 2014b).高原平臺的加熱作用在高層表現(xiàn)的范圍更廣,覆蓋了整個高原地區(qū),且強(qiáng)度也比喜馬拉雅山脈的加熱大,異常中心位于高原的西部(這主要由夏季高原地表感熱大值區(qū)位于西北部決定的(圖3a和3b)).由于平臺加熱的范圍和強(qiáng)度都較大,這可能使得高原平臺的地形加熱可以影響到更大尺度范圍的海陸熱力差異以及經(jīng)圈環(huán)流,進(jìn)而造成之前圖3g所示的高原平臺加熱可以影響更遠(yuǎn)處的熱帶海洋上的西南季風(fēng)環(huán)流.此外,當(dāng)去除整個高原的地形加熱時,對高層的影響更為顯著,高層的暖中心降低超過4 K,南亞高壓的強(qiáng)度也被大大削弱(圖6d和6g).

    圖6 夏季400~200 hPa質(zhì)量權(quán)重溫度場(陰影,單位:K)和風(fēng)場(矢量,單位:m·s-1)(a) CTL, (b) HIM_NS, (c) PL_NS, (d) TP_NS, (e) CTL減HIM_NS,(f) CTL減PL_NS,(g) CTL減TP_NS.其中(e—g)中深色矢量箭頭和打點(diǎn)區(qū)分別表示風(fēng)場和氣溫差異通過95%信度檢驗(yàn)區(qū).Fig.6 Mass-weighted vertical mean temperature (shading, unit: K) and wind (vectors, unit: m·s-1) for 400~200 hPa in summer time(a) CTL, (b) HIM_NS, (c) PL_NS, (d) TP_NS, (e) CTL minus HIM_NS, (f) CTL minus PL_NS, and (g) CTL minus TP_NS. Dotted regions in (e—g) denote the statistical significance of temperature anomalies above 95% level, and wind differences significant above 95% level are plotted in black color.

    圖7 緯向(85°E—90°E)平均的垂直速度(陰影,單位:0.02 Pa·s-1)和垂直環(huán)流v&50×(-ω)(矢量箭頭,單位:0.02 m·Pa·s-2)的剖面圖(a) CTL; (b) HIM_NS; (c) PL_NS; (d) TP_NS.Fig.7 Pressure-latitude cross-sections of the vertical circulation (vectors, v and -50×ω, unit: 0.02 m·Pa·s-2) and vertical velocity (shading, unit: 0.02 Pa·s-1) averaged from 85°E to 90°E for the experiments

    圖8 高原不同地形加熱條件下的各季風(fēng)指數(shù)(已標(biāo)準(zhǔn)化)(a) 南亞夏季風(fēng); (b) 東亞夏季風(fēng).Fig.8 Standardized monsoon indices in different experiments(a) South Asian summer monsoon (SASM); (b) East Asian summer monsoon (EASM).

    3.2 季風(fēng)垂直結(jié)構(gòu)分布

    高原不同區(qū)域的熱力強(qiáng)迫可引起對流層高層溫度場結(jié)構(gòu)的不同變化,隨之而變的溫度經(jīng)向梯度必然也會影響季風(fēng)的垂直環(huán)流結(jié)構(gòu).圖7是南亞季風(fēng)區(qū)緯向(85°E—90°E)平均的垂直環(huán)流場,由參考試驗(yàn)CTL的模擬結(jié)果可知,南亞夏季風(fēng)的垂直上升氣流分為南北兩支,南支位于20°N左右,而北支位于喜馬拉雅山脈的斜坡上.另外,高原北坡也有一上升支,但是其強(qiáng)度比南亞季風(fēng)的北支分量明顯偏弱,原因是北部的上升氣流主要來自高緯地區(qū)的干冷空氣,上升運(yùn)動伴隨著的降水及其釋放的凝結(jié)潛熱加熱較弱,這種降水加熱和上升運(yùn)動的正反饋?zhàn)饔脹]有高原南坡上空那樣顯著.圖7b是去除喜馬拉雅山脈斜坡地形加熱后的試驗(yàn)結(jié)果,和Wu等(2012)的結(jié)論一致,南亞季風(fēng)的北支分量基本消失,說明高原南部斜坡的熱力強(qiáng)迫的確對南亞夏季風(fēng)北支的形成和維持起主導(dǎo)作用.對于南支分量,海陸熱力差異對其起決定作用(Wu et al., 2012),但是不難發(fā)現(xiàn)去除斜坡加熱后,南支分量相比于參考試驗(yàn)CTL明顯向北拓寬.這一現(xiàn)象與繞高原西側(cè)南下的干冷空氣減弱有關(guān),喜馬拉雅山脈斜坡無加熱大大削弱了這支西北氣流(圖4f),從而使得南部海洋上的暖濕氣流北推,最終造成南亞夏季風(fēng)的南支分量向北拓展.由圖4f的降水差異可知喜馬拉雅山脈無加熱反而增加了印度中部(20°N—25°N)的降水量,但這并不能認(rèn)為是整個南亞夏季風(fēng)的增強(qiáng),其實(shí)這也是剛剛解釋的南支分量向北拓展造成的結(jié)果,而事實(shí)上這時候南亞季風(fēng)的北支分量已經(jīng)消失.如果只去除高原平臺的加熱(圖7c),可以發(fā)現(xiàn)南亞夏季風(fēng)的南北支分量都依然存在,跟CTL模擬的結(jié)構(gòu)基本保持一致,只是上升氣流的強(qiáng)度有所改變.值得注意的是,無論是單獨(dú)去除喜馬拉雅山脈斜坡的加熱還是單獨(dú)去除高原平臺的加熱,均造成各自局地垂直上升環(huán)流的削弱或消失,同時還分別使得對應(yīng)保留地形加熱的平臺或斜坡上的垂直上升氣流增強(qiáng),如圖7b所示高原平臺上的垂直上升氣流以及圖7c中高原南部斜坡上的上升氣流都比圖7a中CTL的結(jié)果要偏強(qiáng).不難發(fā)現(xiàn),高原南部斜坡上的感熱氣泵效應(yīng)使得氣流在斜坡上升而同時在平臺上下沉,而平臺上的感熱氣泵作用使得氣流在平臺上為上升運(yùn)動同時在斜坡下沉,所以兩部分加熱同時存在會相互抵消一定強(qiáng)度的氣流上升運(yùn)動(圖7a),王在志(2005)基于大氣環(huán)流模式的理想地形試驗(yàn)也得到了類似的結(jié)果.

    當(dāng)整個高原地區(qū)的感熱加熱去除后,南亞夏季風(fēng)的北支分量基本被抑制,此時南支分量同樣向北拓展變寬,與HIM_NS的結(jié)果非常類似;當(dāng)然,不同于HIM_NS的是高原北部加熱的去除使得北坡的上升氣流也隨之消失(圖7d).另一方面,關(guān)于東亞夏季風(fēng)垂直環(huán)流的響應(yīng)相對簡單一些,與圖4中國東部低層偏南氣流的響應(yīng)類似,無論是喜馬拉雅山脈還是高原平臺,它們的地形加熱作用都使得東亞夏季風(fēng)垂直上升區(qū)向北加強(qiáng),推進(jìn)季風(fēng)系統(tǒng)的向北發(fā)展(圖略).

    3.3 不同季風(fēng)指數(shù)的對比

    為了進(jìn)一步闡明高原不同區(qū)域的地形加熱對亞洲夏季風(fēng)強(qiáng)度的影響差異,我們還計算了幾種常用的夏季風(fēng)指數(shù)(圖8).SA1為夏季南亞區(qū)域(5°N—20°N, 60°E—110°E)平均的高低空(850 hPa減200 hPa)緯向風(fēng)差異(Webster et al., 1998);SA2為印度季風(fēng)區(qū)及其東側(cè)區(qū)域(10°N—30°N, 70°E—110°E)平均的高低空(850 hPa減200 hPa)經(jīng)向風(fēng)差異(Goswami et al., 1999);SA3為南印度(5°N—15°N, 60°E—80°E)和北印度(20°N—30°N, 70°E—90°E)之間低空850 hPa緯向風(fēng)差異(Wang et al., 2001).以上三種季風(fēng)指數(shù)代表了南亞夏季風(fēng)區(qū)不同的動力學(xué)特征,例如分別表征季風(fēng)緯向環(huán)流的反轉(zhuǎn)(SA1)、季風(fēng)經(jīng)向哈得來環(huán)流(SA2)以及印度季風(fēng)槽強(qiáng)度(SA3).另外,根據(jù)Wang等(2008b)的定義,東亞夏季風(fēng)指數(shù)可利用低空偏南風(fēng)來描述,EA1為夏季中國東部(20°N—40°N, 110°E—120°E)平均的850 hPa經(jīng)向風(fēng)速,這代表了向北的水汽輸送強(qiáng)度;EA2為中國東部南(20°N—30°N, 110°E—125°E)北(30°N—40°N, 110°E—125°E)之間850 hPa經(jīng)向風(fēng)差異,代表了東亞地區(qū)水汽輻合的強(qiáng)度及其南北位置.如果各組敏感性試驗(yàn)得出的季風(fēng)指數(shù)相比于參考試驗(yàn)變化的幅度越大,表明該區(qū)域的地形加熱對季風(fēng)的影響就越強(qiáng).由圖8可知,除了SA3印度季風(fēng)槽強(qiáng)度指數(shù)隨高原地形加熱的減弱而增大,其他季風(fēng)指數(shù)均變小(表明高原地形加熱對南亞和東亞夏季風(fēng)均有增強(qiáng)作用),其中整個高原地形加熱去除后(TP_NS)指數(shù)的變化幅度最大,而各個指數(shù)在HIM_NS和PL_NS試驗(yàn)中的變化特征卻也不盡相同.通過對比HIM_NS和PL_NS試驗(yàn)中的各季風(fēng)指數(shù)表明,高原南部喜馬拉雅山脈斜坡加熱對南亞夏季風(fēng)緯向風(fēng)的垂直切變影響更大(SA1HIM_NS< SA1PL_NS);而平臺加熱對經(jīng)向哈得來環(huán)流的調(diào)控更強(qiáng)(SA2HIM_NS> SA2PL_NS),經(jīng)向環(huán)流的范圍更大,所以這也與前文的分析結(jié)果吻合(平臺加熱及其對環(huán)流的影響范圍更廣).至于SA3為什么會在高原地形加熱的削弱時反而增大,由前文圖4及圖5可知,高原的地形加熱削弱后,孟加拉灣北上的氣流爬坡減弱,而沿著高原南坡向西的繞流卻增強(qiáng)了,所以造成印度北部的季風(fēng)槽加強(qiáng);且斜坡加熱造成的SA3變化幅度明顯強(qiáng)于高原平臺作用,這是因?yàn)槿コ逼录訜岷蠡疽种屏藲饬鞯呐榔?,進(jìn)而全部轉(zhuǎn)為繞流,大大加深了印度季風(fēng)槽.在東亞區(qū)域,因?yàn)樾逼录訜岵粌H加強(qiáng)了中國東部南方的偏南氣流,而且也激發(fā)了北方異常的偏北氣流南下(圖4f),南北異常氣流的交匯使得EA1指數(shù)相比CTL變化幅度較?。欢脚_加熱沒有異常偏北氣流南下的響應(yīng)(圖4g),所以平臺加熱PL_NS中EA1指數(shù)的變化幅度要大于HIM_NS.但是對于表示中國東部水汽輻合強(qiáng)度的EA2指數(shù),HIM_NS的變化幅度要比平臺加熱PL_NS大,造成的降水差異也更為顯著(圖4f).

    綜上所述,高原的大地形熱力強(qiáng)迫對亞洲夏季風(fēng)的形成和維持的確有很重要的影響,尤其是喜馬拉雅山脈斜坡的加熱,它對南亞夏季風(fēng)北支分量的形成起著決定性的作用.另外,多種季風(fēng)指數(shù)的比較也表明了高原不同區(qū)域的地形加熱分別對南亞夏季風(fēng)和東亞夏季風(fēng)有著不一樣的影響.

    4 結(jié)論與討論

    本文利用較高分辨率的WRF區(qū)域氣候模式,模擬研究了高原不同區(qū)域(高原南部喜馬拉雅山脈斜坡、高原平臺以及整個高原地區(qū))的地形加熱對南亞夏季風(fēng)和東亞夏季風(fēng)所對應(yīng)的大氣環(huán)流、降水以及溫度場結(jié)構(gòu)的不同影響.主要結(jié)論如下:

    (1) 高原南部喜馬拉雅山脈的斜坡地形加熱對其周圍局地的環(huán)流形勢和降水影響十分明顯,是南亞夏季風(fēng)北支分量形成和維持的主導(dǎo)因子,也是斜坡上氣流爬坡和降水發(fā)生的根本原因.斜坡加熱對東亞夏季風(fēng)也有明顯的增強(qiáng)作用,它不僅加強(qiáng)了中國東部低空西南季風(fēng)環(huán)流,還會引起北部南下的異常干冷空氣的響應(yīng).此外,斜坡上的地形加熱作用也是對流層高層暖中心增強(qiáng)及其位置維持在斜坡上空的一個重要原因.

    (2) 高原平臺加熱對季風(fēng)環(huán)流和降水的影響雖然沒有喜馬拉雅山脈斜坡加熱那么顯著,但是對南亞夏季風(fēng)的影響范圍更廣,對經(jīng)向哈得來環(huán)流影響更明顯,能夠調(diào)控高原以外更遠(yuǎn)處熱帶洋面上的西南季風(fēng)環(huán)流.同樣,平臺加熱對東亞夏季風(fēng)也有增強(qiáng)作用.

    (3) 當(dāng)單獨(dú)去除高原南坡喜馬拉雅山脈斜坡加熱時,將造成斜坡局地垂直上升流的消失,同時又使得高原平臺上的垂直上升氣流增強(qiáng);反之當(dāng)單獨(dú)去除高原平臺加熱時,將造成平臺局地的垂直上升流消失,同時使得高原南部斜坡上的上升支加強(qiáng).

    (4) 通過比較高原不同地形加熱條件下的多種季風(fēng)指數(shù),表明了高原地形加熱對南亞和東亞夏季風(fēng)均有增強(qiáng)作用(地表加熱的去除使得季風(fēng)指數(shù)均減小),但是高原不同區(qū)域的地形加熱分別對兩類夏季風(fēng)子系統(tǒng)又會產(chǎn)生不一樣的影響.比如,高原南部的斜坡加熱對南亞季風(fēng)緯向環(huán)流的垂直切變指數(shù)(SA1)影響更大,而高原平臺加熱對表征經(jīng)向哈得來環(huán)流的SA2指數(shù)影響更顯著;平臺加熱對東亞夏季風(fēng)EA1指數(shù)的影響幅度更大,而斜坡加熱對表示中國東部水汽輻合強(qiáng)度的EA2指數(shù)影響更顯著.

    文章結(jié)論指出高原地形加熱是局地降水發(fā)生的基本條件,一旦地表感熱加熱被抑制,降水就基本消失了.換句話說,高原感熱加熱的去除同時也削弱了與降水相關(guān)的凝結(jié)潛熱加熱,使得整個大氣熱源得到削弱.因?yàn)橄募靖咴值亟邓饕l(fā)生在斜坡上(圖4a),當(dāng)斜坡上的地表感熱加熱去除后造成的大氣熱源異常要明顯強(qiáng)于平臺加熱去除后的大氣熱源異常,所以這也可以進(jìn)一步解釋為什么高原斜坡加熱的氣候響應(yīng)要強(qiáng)于平臺加熱.在今后的研究工作中有必要繼續(xù)探討高原上空單獨(dú)的潛熱加熱對亞洲夏季風(fēng)產(chǎn)生的影響.此外,本文的研究結(jié)果是基于WRF區(qū)域大氣模式得到的,模式中海表溫度是人為給定的強(qiáng)迫場,沒有考慮大氣對海洋的影響.然而,高原地形加熱導(dǎo)致的海面風(fēng)場變化(如圖4h)必然會引起海表溫度的改變,同時海溫的改變又會進(jìn)一步反饋于大氣,并造成季風(fēng)環(huán)流和降水的異常.那么,海氣耦合在高原影響亞洲夏季風(fēng)的過程中會發(fā)揮怎樣的調(diào)控作用?這種作用顯著嗎?針對這些問題,下一步利用海氣耦合模式來進(jìn)一步模擬研究高原大地形加熱對亞洲夏季風(fēng)的影響.

    致謝 感謝中國科學(xué)院大氣物理研究所吳國雄院士對本研究提出的寶貴建議.

    Boos W R, Kuang Z M. 2010. Dominant control of the South Asian monsoon by orographic insulation versus plateau heating.Nature, 463(7278): 218-223.

    Boos W R, Kuang Z M. 2013. Sensitivity of the South Asian monsoon to elevated and non-elevated heating.ScientificReports, 3: 1192, doi: 10.1038/srep01192.

    Flohn H. 1957. Large-scale aspects of the “summer monsoon” in South and East Asia.JournaloftheMeteorologicalSocietyofJapan, 75: 180-186.

    Goswami B N, Krishnamurthy V, Annmalai H. 1999. A broad-scale circulation index for the interannual variability of the Indian summer monsoon.QuarterlyJournalofRoyalMeteorologicalSociety, 125(554): 611-633. He B. 2012. Numerical simulation and mechanism study on the impacts of Tibetan Plateau thermodynamic forcing over South Asian summer monsoon [Ph. D. thesis] (in Chinese). Nanjing: Nanjing University of Information Science & Technology. He J H, Xu H M, Zhong S S, et al. 2011. Characteristics of Atmospheric Heat Source over the Tibetan Plateau and Its Impacts and Possible Mechanisms (in Chinese). Beijing: China Meteorological Press.

    Herring S C, Hoerling M P, Peterson T C, et al. 2014. Explaining extreme events of 2013 from a climate perspective.BulletinoftheAmericanMeteorologicalSociety, 95(9): S1-S104.

    Kim E J, Hong S Y. 2010. Impact of air-sea interaction on East Asian summer monsoon climate in WRF.JournalofGeophysicalResearch, 115: D19118, doi: 10.1029/2009JD013253. Li Y Q. 2003. Surface heating in the Tibetan Plateau and general circulation over it and their relations with the prediction of drought-flood at its eastern side.ChineseJournalofAtmosphericSciences(in Chinese), 27(1): 107-114. Liang X Y, Liu Y M, Wu G X. 2006. Roles of tropical and subtropical land-sea distribution and the Qinghai-Xizang Plateau in the formation of the Asian summer monsoon.ChineseJournalofGeophysics(in Chinese), 49(4): 983-992. Mao J Y, Wu G X. 2006. Impacts of anomalies of thermal state over the Qinghai-Xizang Plateau and sea surface temperature on interannual variability of the Asian monsoon seasonal transition.ChineseJournalofGeophysics(in Chinese), 49(5): 1279-1287. Molnar P, Boos W R, Battisti D S. 2010. Orographic controls on climate and paleoclimate of Asia: Thermal and mechanical roles for the Tibetan Plateau.AnnualReviewofEarthandPlanetarySciences, 38(1): 77-102. Qian Y F, Yan H, Wang Q Q, et al. 1988. Numerical Study of the Effect of Topography in the Planetary Atmosphere (in Chinese). Beijing: Science Press, 105. Skamarock W C, Klemp J B, Dudhia J, et al. 2008. A description of the advanced research WRF version 3. NCAR Technical Note NCAR/TN-475+STR, DOI: 10.5065/D68S4MVH.

    Wang B, Wu R G, Lau K M. 2001. Interannual variability of the Asian summer monsoon: contrasts between the Indian and the western North Pacific-East Asian monsoons.JournalofClimate, 14(20): 4073-4090. Wang B, Bao Q, Hoskins B, et al. 2008a. Tibetan Plateau warming and precipitation changes in East Asia.GeophysicalResearchLetters, 35: L14702, doi: 10.1029/2008GL034330.

    Wang B, Wu Z W, Li J P, et al. 2008b. How to measure the strength of the East Asian summer monsoon.JournalofClimate, 21(17): 4449-4463.

    Wang Z Q, Duan A M, Wu G X. 2014a. Impacts of boundary layer parameterization schemes and air-sea coupling on WRF simulation of the East Asian summer monsoon.ScienceChina:EarthScience, 57(7): 1480-1493, doi: 10.1007/s11430-013-4801-4. Wang Z Q, Duan A M, Wu G X. 2014b. Time-lagged impact of spring sensible heat over the Tibetan Plateau on the summer rainfall anomaly in East China: case studies using the WRF model.ClimateDynamics, 42(11-12): 2885-2898, doi: 10.1007/s00382-013-1800-2.

    Wang Z Q, Duan A M, Wu G X, et al. 2015. Mechanism for occurrence of precipitation over the southern slope of the Tibetan Plateau without local surface heating.InternationalJournalofClimatology, doi: 10.1002/joc.4609.

    Webster P J. 2006. The coupled monsoon system.∥ The Asian Monsoon. Berlin Heidelberg: Springer, 3-66.

    Wu G X, Li W P, Guo H. 1997. Sensible heat driven air-pump over the Tibetan-Plateau and its impacts on the Asian Summer Monsoon.∥ Ye D Z ed. Collections on the Memory of Zhao Jiuzhang (in Chinese). Beijing: Science Press.

    Wu G X, Zhang Y S. 1998. Tibetan Plateau forcing and the timing of the monsoon onset over South Asia and the South China Sea.MonthlyWeatherReview, 126: 913-927.

    Wu G X, Liu Y M, Zhang Q, et al. 2007. The influence of mechanical and thermal forcing by the Tibetan Plateau on Asian climate.JournalofHydrometeorology, 8(4): 770-789, doi: 10.1175/JHM609.1.

    Wang Z Z. 2005. Numerical study of the Tibetan Plateau effects on weather and climate [Ph. D. thesis] (in Chinese). Beijing: Graduate University of Chinese Academy of Sciences.

    Wu G X, Liu Y M, He B, et al. 2012. Thermal controls on the Asian Summer Monsoon.ScientificReports, 2: 404, doi: 10.1038/srep00404.

    Xu X D, Lu C G, Shi X H, et al. 2010. Large-scale topography of China: A factor for the seasonal progression of the Meiyu rainband?JournalofGeophysicalResearch, 115: D02110, doi: 10.1029/2009JD01244.

    Yang H W, Wang B, Wang B. 2012. Reduction of systematic biases in regional climate downscaling through ensemble forcing.ClimateDynamics, 38(3-4): 655-665.

    Ye D Z, Gao Y X, Zhou M Y, et al. 1979. Meteorology of the Qinghai-Xizang (Tibet) Plateau (in Chinese). Beijing: Science Press. Yu E T, Wang H J, Gao Y Q, et al. 2011. Impacts of cumulus convective parameterization schemes on summer monsoon precipitation simulation over China.ActaMeteorologicaSinica, 25(5): 581-592. Yu L J, Hu D X. 2008. Role of snow depth in spring of Tibetan Plateau in onset of South China Sea summer monsoon.ChineseJournalofGeophysics(in Chinese), 51(6): 1682-1694.

    Zhang Y C, Qian Y F. 1999. Numerical studies on the effects of the critical height of Qinghai-Xizang Plateau uplift on the atmosphere.ActaMeteorologicaSinica(in Chinese), 57(2): 157-167. Zhao P, Chen L X. 2001. Climatic features of atmospheric heat source/sink over the Qinghai-Xizang Plateau in 35 years and its relation to rainfall in China.ScienceinChinaSeriesD:EarthSciences, 44(9): 858-864.

    附中文參考文獻(xiàn)

    何編. 2012. 青藏高原熱動力強(qiáng)迫對南亞夏季風(fēng)影響的數(shù)值模擬研究和機(jī)理分析[博士論文]. 南京: 南京信息工程大學(xué).

    何金海, 徐海明, 鐘珊珊等. 2011. 青藏高原大氣熱源特征及其影響和可能機(jī)制. 北京: 氣象出版社.

    李躍清. 2003. 青藏高原地面加熱及上空環(huán)流場與東側(cè)旱澇預(yù)測的關(guān)系. 大氣科學(xué), 27(1): 107-114.

    梁瀟云, 劉屹岷, 吳國雄. 2006. 熱帶、副熱帶海陸分布與青藏高原在亞洲夏季風(fēng)形成中的作用. 地球物理學(xué)報, 49(4): 983-992.

    毛江玉, 吳國雄. 2006. 青藏高原熱狀況和海溫異常對亞洲季風(fēng)季節(jié)轉(zhuǎn)換年際變化的影響. 地球物理學(xué)報, 49(5): 1279-1287.

    錢永甫, 顏宏, 王謙謙等. 1988. 行星大氣中地形效應(yīng)的數(shù)值研究. 北京: 科學(xué)出版社, 105.

    王在志. 2005. 青藏高原天氣氣候效應(yīng)的數(shù)值模擬研究[博士論文]. 北京: 中國科學(xué)院研究生院.

    吳國雄, 李偉平, 郭華. 1997. 青藏高原感熱氣泵和亞洲夏季風(fēng).∥趙九章紀(jì)念文集. 北京: 科學(xué)出版社.

    葉篤正, 高由禧, 周明煜等. 1979. 青藏高原氣象學(xué). 北京: 科學(xué)出版社.

    于樂江, 胡敦欣. 2008. 青藏高原春季積雪在南海夏季風(fēng)爆發(fā)過程中的作用. 地球物理學(xué)報, 51(6): 1682-1694.

    張耀存, 錢永甫. 1999. 青藏高原隆升作用于大氣臨界高度的數(shù)值研究. 氣象學(xué)報, 57(2): 157-167.

    (本文編輯 何燕)

    Influences of thermal forcing over the slope/platform of the Tibetan Plateau on Asian summer monsoon: Numerical studies with the WRF model

    WANG Zi-Qian1,2, DUAN An-Min2, LI Mao-Shan3, HE Bian2

    1SchoolofAtmosphericSciences,SunYat-SenUniversity,Guangzhou510275,China2StateKeyLaboratoryofNumericalModelingforAtmosphericSciencesandGeophysicalFluidDynamics,InstituteofAtmosphericPhysics,ChineseAcademyofSciences,Beijing100029,China3KeyLaboratoryforLandSurfaceProcessandClimateChangeinColdandAridRegions,ChineseAcademyofSciences,Lanzhou730000,China

    With global warming, the evolution of Asian summer monsoon (ASM) becomes more complicated, while the thermal effect of the Tibetan Plateau (TP) is an important forcing to the variability of weather and climate in ASM areas. The regional climate model (RCM) is a useful tool in research of the regional climate change and has higher resolutions that can represent topography and land surface processes more accurately compared with atmospheric general circulation models (AGCMs). In this study, the impact of thermal forcing over different terrains of TP on ASM (including South and East ASM (SASM and EASM)) is investigated using the Weather Research and Forecasting (WRF) models. Results indicate that the local circulation and precipitation around the Himalayas (HIM) are significantly influenced by the surface heating over the HIM′s southern slope, which is a dominant factor for the formation of the north branch of SASM. Meanwhile, the climbing moist airflow and precipitation over the southern slope of TP are mainly induced by HIM′s thermal forcing. Due to HIM′s slope heating, the upper-level troposphere warm center is steadily located over the HIM area and the EASM is also intensified obviously (characterized both by the enhanced low-level southwesterly over East China and the enhanced southward anomalous dry-cold northerly). As to the surface heating over the TP′s platform region (PL), although its influence on the summer monsoon circulation and precipitation is weaker than the HIM′s, it induces a wider response of SASM and a stronger influence on the meridional Hadley circulation. The PL′s heating is able to regulate low-level southwesterly over the remote tropical ocean. Furthermore, the comparison of multiple monsoon indices reveals that both HIM sloping heating and PL heating can intensify SASM or EASM, but the impacts of thermal forcing over different terrains on the two subsystems of ASM are distinct.

    Tibetan Plateau; Thermal forcing; Asian summer monsoon; WRF model

    10.6038/cjg20160904.

    廣東省自然科學(xué)基金(2015A030310224),國家重大科學(xué)研究計劃(2014CB953900),國家自然科學(xué)基金(91337216, 41605038)共同資助.

    王子謙,男,助理研究員,主要從事青藏高原氣候動力學(xué)及海-氣相互作用研究.E-mail:wangziq5@mail.sysu.edu.cn

    10.6038/cjg20160904

    P461

    2015-09-27,2016-07-08收修定稿

    王子謙, 段安民, 李茂善等. 2016. 基于WRF模式的青藏高原斜坡和平臺加熱影響亞洲夏季風(fēng)的模擬研究. 地球物理學(xué)報,59(9):3175-3187,

    Wang Z Q, Duan A M, Li M S, et al. 2016. Influences of thermal forcing over the slope/platform of the Tibetan Plateau on Asian summer monsoon: Numerical studies with the WRF model.ChineseJ.Geophys. (in Chinese),59(9):3175-3187,doi:10.6038/cjg20160904.

    猜你喜歡
    喜馬拉雅山脈季風(fēng)南亞
    喜瑪拉雅山脈
    江南詩(2023年1期)2023-04-15 12:40:40
    斑頭雁:看我如何翻越高山
    南亞磷肥市場需東繼續(xù)提升
    中國熱科院南亞所 研發(fā)成功菠蘿機(jī)械化種植機(jī)
    綠水青山圖——海洋季風(fēng)的贊歌
    The Yeti
    南亞煉油化工市場一覽
    戶撒刀
    南亞高壓對西藏夏季降水的影響
    西藏科技(2016年5期)2016-09-26 12:16:41
    萬馬奔騰
    中華奇石(2015年9期)2015-07-09 18:32:58
    国产 一区精品| 日韩不卡一区二区三区视频在线| 天堂网av新在线| 九九爱精品视频在线观看| 18禁动态无遮挡网站| 国产精品久久久久久精品古装| 精品99又大又爽又粗少妇毛片| 乱系列少妇在线播放| 色视频www国产| 特大巨黑吊av在线直播| av线在线观看网站| 国产精品99久久99久久久不卡 | 香蕉精品网在线| 18禁裸乳无遮挡动漫免费视频 | 一级爰片在线观看| 菩萨蛮人人尽说江南好唐韦庄| 尾随美女入室| 一级av片app| 91精品国产九色| 久久97久久精品| 亚洲国产精品专区欧美| 亚洲第一区二区三区不卡| 色综合色国产| 久久久亚洲精品成人影院| 插逼视频在线观看| 国产精品麻豆人妻色哟哟久久| 伊人久久国产一区二区| 美女脱内裤让男人舔精品视频| 免费播放大片免费观看视频在线观看| 日本爱情动作片www.在线观看| 日韩一区二区三区影片| 午夜免费观看性视频| 亚洲欧美精品自产自拍| 国国产精品蜜臀av免费| 色婷婷久久久亚洲欧美| 国产 一区 欧美 日韩| 成人免费观看视频高清| 美女国产视频在线观看| 婷婷色综合www| 国产一区二区亚洲精品在线观看| 熟女电影av网| 女人被狂操c到高潮| 美女内射精品一级片tv| 一级毛片我不卡| videossex国产| 色视频www国产| 国产伦精品一区二区三区四那| 欧美97在线视频| 欧美97在线视频| 国产在线男女| 人妻 亚洲 视频| 日日摸夜夜添夜夜添av毛片| 午夜免费观看性视频| 国产欧美亚洲国产| 高清视频免费观看一区二区| 嫩草影院精品99| 69av精品久久久久久| 最近2019中文字幕mv第一页| 日韩欧美精品免费久久| 丝袜脚勾引网站| 性色av一级| 免费黄色在线免费观看| 国产精品三级大全| 亚洲国产av新网站| 国产一区亚洲一区在线观看| 麻豆国产97在线/欧美| 18禁在线播放成人免费| 日本爱情动作片www.在线观看| 免费观看无遮挡的男女| 午夜老司机福利剧场| 一本一本综合久久| 国产人妻一区二区三区在| 国产精品伦人一区二区| 精品99又大又爽又粗少妇毛片| 欧美xxxx性猛交bbbb| 美女高潮的动态| 久久精品国产亚洲网站| 美女高潮的动态| 最近2019中文字幕mv第一页| 国产一级毛片在线| 国产精品蜜桃在线观看| 中文字幕制服av| 乱码一卡2卡4卡精品| 在线播放无遮挡| 黄片无遮挡物在线观看| 18禁在线无遮挡免费观看视频| av在线app专区| 国产精品一区二区三区四区免费观看| 日韩精品有码人妻一区| 亚洲精品国产av成人精品| 99热这里只有是精品50| 亚洲一区二区三区欧美精品 | 国产高清三级在线| 午夜日本视频在线| 一区二区三区四区激情视频| 最近最新中文字幕大全电影3| 亚洲色图av天堂| 亚洲国产av新网站| 人妻少妇偷人精品九色| 亚洲欧美日韩东京热| 少妇人妻 视频| 男男h啪啪无遮挡| 久久精品国产a三级三级三级| 高清午夜精品一区二区三区| 免费播放大片免费观看视频在线观看| av一本久久久久| 久久久久网色| 高清视频免费观看一区二区| 久久久久久伊人网av| 亚洲高清免费不卡视频| 欧美变态另类bdsm刘玥| 久久国产乱子免费精品| 国产成人免费无遮挡视频| 精品亚洲乱码少妇综合久久| 欧美精品一区二区大全| 亚洲精品国产色婷婷电影| 欧美丝袜亚洲另类| 成年版毛片免费区| 国产亚洲91精品色在线| 男人狂女人下面高潮的视频| 日本-黄色视频高清免费观看| 久久人人爽人人爽人人片va| 嘟嘟电影网在线观看| 国产毛片a区久久久久| 欧美bdsm另类| 毛片女人毛片| 亚洲国产精品999| 亚洲精品日本国产第一区| 精品人妻一区二区三区麻豆| 中文字幕av成人在线电影| 偷拍熟女少妇极品色| 狠狠精品人妻久久久久久综合| 亚洲成色77777| 日韩大片免费观看网站| 狂野欧美激情性bbbbbb| 日韩国内少妇激情av| 精品一区二区免费观看| 狂野欧美激情性bbbbbb| 久久女婷五月综合色啪小说 | 热re99久久精品国产66热6| 大香蕉久久网| 2018国产大陆天天弄谢| 欧美成人精品欧美一级黄| 国产白丝娇喘喷水9色精品| 亚洲av.av天堂| 国产老妇伦熟女老妇高清| 久久久欧美国产精品| 边亲边吃奶的免费视频| 爱豆传媒免费全集在线观看| 少妇人妻精品综合一区二区| 亚洲欧美成人精品一区二区| 国产真实伦视频高清在线观看| 久久久成人免费电影| 性插视频无遮挡在线免费观看| 下体分泌物呈黄色| kizo精华| 国产一区二区在线观看日韩| 久久综合国产亚洲精品| 亚洲欧美成人综合另类久久久| 久久影院123| 中国三级夫妇交换| av在线天堂中文字幕| 国产亚洲91精品色在线| 欧美潮喷喷水| 久久精品熟女亚洲av麻豆精品| 欧美成人午夜免费资源| 成人亚洲精品av一区二区| 男人狂女人下面高潮的视频| 亚洲人成网站在线播| 日韩一区二区视频免费看| 成人亚洲精品一区在线观看 | 欧美3d第一页| 日韩一区二区视频免费看| 成人亚洲欧美一区二区av| 免费看a级黄色片| 国产成人福利小说| 男插女下体视频免费在线播放| 99热这里只有是精品50| 自拍偷自拍亚洲精品老妇| 国产毛片在线视频| 热99国产精品久久久久久7| 18+在线观看网站| 亚洲自拍偷在线| 一区二区三区免费毛片| 国产精品国产av在线观看| 最新中文字幕久久久久| 国产黄色视频一区二区在线观看| 黄色日韩在线| 欧美少妇被猛烈插入视频| 久热久热在线精品观看| 亚洲精品aⅴ在线观看| av一本久久久久| 亚洲不卡免费看| 成年女人看的毛片在线观看| 我的女老师完整版在线观看| 中国三级夫妇交换| 国产成人a区在线观看| 97超碰精品成人国产| 在线a可以看的网站| 日韩欧美精品v在线| 九色成人免费人妻av| 卡戴珊不雅视频在线播放| 狂野欧美激情性bbbbbb| 在线亚洲精品国产二区图片欧美 | kizo精华| 91久久精品电影网| av又黄又爽大尺度在线免费看| 欧美日韩国产mv在线观看视频 | av在线app专区| 日韩一本色道免费dvd| 欧美日韩视频高清一区二区三区二| 久久久久久久大尺度免费视频| 大陆偷拍与自拍| 色哟哟·www| 成人二区视频| 成年女人看的毛片在线观看| 日本爱情动作片www.在线观看| 极品少妇高潮喷水抽搐| 午夜福利在线观看免费完整高清在| 久久久久久久午夜电影| 欧美日韩亚洲高清精品| 人妻夜夜爽99麻豆av| 久久综合国产亚洲精品| 亚洲av在线观看美女高潮| 一级毛片黄色毛片免费观看视频| 色视频在线一区二区三区| 99久久人妻综合| 人妻夜夜爽99麻豆av| 欧美日韩视频高清一区二区三区二| 久久久久久久久久久丰满| 国产高潮美女av| 国产精品女同一区二区软件| 三级国产精品片| 男人和女人高潮做爰伦理| 亚洲精品国产色婷婷电影| 精品久久久精品久久久| 国产高清三级在线| 全区人妻精品视频| 欧美精品国产亚洲| 大片免费播放器 马上看| 国产探花在线观看一区二区| 中文字幕av成人在线电影| 97在线视频观看| 久久久久久久大尺度免费视频| 看十八女毛片水多多多| 精品午夜福利在线看| 最近2019中文字幕mv第一页| 欧美变态另类bdsm刘玥| 国产极品天堂在线| 亚洲精品久久久久久婷婷小说| 久久ye,这里只有精品| 建设人人有责人人尽责人人享有的 | 全区人妻精品视频| 91久久精品国产一区二区三区| 91精品国产九色| 日韩,欧美,国产一区二区三区| 精华霜和精华液先用哪个| 欧美成人a在线观看| 亚洲精品一区蜜桃| 久久久久久久久久久免费av| 日日撸夜夜添| 国产人妻一区二区三区在| 99久久精品一区二区三区| 中文字幕av成人在线电影| 国产成年人精品一区二区| 日日摸夜夜添夜夜爱| 亚洲精品国产av蜜桃| 伊人久久精品亚洲午夜| 日日撸夜夜添| 久久精品国产亚洲网站| 成人亚洲精品一区在线观看 | 欧美日韩精品成人综合77777| 国产成人免费无遮挡视频| 国产 一区 欧美 日韩| 身体一侧抽搐| 亚洲国产欧美人成| 一级av片app| 亚洲国产精品999| 99热这里只有是精品在线观看| 日韩一区二区视频免费看| 一级毛片黄色毛片免费观看视频| 在线a可以看的网站| 亚洲欧美一区二区三区黑人 | 中文乱码字字幕精品一区二区三区| 人人妻人人看人人澡| 精品国产乱码久久久久久小说| 久久久久九九精品影院| 国产高清三级在线| 亚洲av欧美aⅴ国产| 激情 狠狠 欧美| 国产精品蜜桃在线观看| 色网站视频免费| 国产精品三级大全| 欧美3d第一页| 国产白丝娇喘喷水9色精品| 国产成人91sexporn| 亚洲国产最新在线播放| 尤物成人国产欧美一区二区三区| 中文天堂在线官网| 成人特级av手机在线观看| 亚洲精品中文字幕在线视频 | 看非洲黑人一级黄片| 王馨瑶露胸无遮挡在线观看| 久久久欧美国产精品| 亚洲伊人久久精品综合| 好男人在线观看高清免费视频| 欧美日本视频| 国产精品av视频在线免费观看| 热99国产精品久久久久久7| 日韩一区二区视频免费看| 婷婷色综合大香蕉| 国产精品嫩草影院av在线观看| 最近最新中文字幕大全电影3| 国产在线男女| 久久99热6这里只有精品| 人妻制服诱惑在线中文字幕| 亚洲av男天堂| 一本一本综合久久| 少妇的逼水好多| 欧美日韩在线观看h| 99热6这里只有精品| 成人午夜精彩视频在线观看| 国产精品福利在线免费观看| 91精品伊人久久大香线蕉| 国产大屁股一区二区在线视频| 国产爱豆传媒在线观看| 一区二区三区四区激情视频| 国产亚洲av片在线观看秒播厂| 成人毛片a级毛片在线播放| 国产精品一及| 中文欧美无线码| 嫩草影院新地址| 成人鲁丝片一二三区免费| 亚洲欧美日韩东京热| 亚洲精品色激情综合| 久久人人爽人人片av| 九草在线视频观看| 国产男人的电影天堂91| 亚洲国产精品成人久久小说| 尾随美女入室| 久久精品久久久久久噜噜老黄| 国内精品宾馆在线| 尾随美女入室| www.色视频.com| 日本一二三区视频观看| av专区在线播放| 一区二区三区精品91| 嫩草影院入口| 99久久精品一区二区三区| 国产精品久久久久久精品电影| 九九爱精品视频在线观看| 亚洲精品成人久久久久久| 国产精品一二三区在线看| 午夜日本视频在线| 一级黄片播放器| 久久热精品热| 三级国产精品片| 日日啪夜夜撸| 国产视频内射| 亚洲成人中文字幕在线播放| 亚洲av不卡在线观看| 少妇被粗大猛烈的视频| 国产欧美日韩一区二区三区在线 | 日韩 亚洲 欧美在线| 久久热精品热| 三级国产精品片| 国产精品偷伦视频观看了| 99热6这里只有精品| 亚洲av国产av综合av卡| 亚洲成人一二三区av| 日韩一本色道免费dvd| 51国产日韩欧美| 午夜福利高清视频| 国产乱来视频区| 亚洲av中文av极速乱| 久久人人爽人人片av| 蜜臀久久99精品久久宅男| 免费黄色在线免费观看| 国产精品国产三级专区第一集| 国产视频内射| 精品久久久久久久久亚洲| 国产91av在线免费观看| 少妇人妻一区二区三区视频| 晚上一个人看的免费电影| 2021天堂中文幕一二区在线观| 久久精品久久久久久久性| 男女边摸边吃奶| 狂野欧美激情性bbbbbb| 亚洲欧美精品专区久久| 男人狂女人下面高潮的视频| 亚洲欧美成人综合另类久久久| 久久精品国产a三级三级三级| 97人妻精品一区二区三区麻豆| 狂野欧美白嫩少妇大欣赏| 国产精品伦人一区二区| 丝瓜视频免费看黄片| 大陆偷拍与自拍| 大码成人一级视频| 成年免费大片在线观看| 亚洲av二区三区四区| 亚洲国产欧美人成| av黄色大香蕉| 久久国产乱子免费精品| 黄色视频在线播放观看不卡| 大香蕉久久网| 成人漫画全彩无遮挡| tube8黄色片| 久久精品久久久久久噜噜老黄| 在线免费观看不下载黄p国产| 亚洲第一区二区三区不卡| 一区二区三区四区激情视频| 久久久欧美国产精品| 在线观看三级黄色| 亚洲精品国产av蜜桃| 亚洲精品亚洲一区二区| 春色校园在线视频观看| 国产一区二区三区综合在线观看 | 国产免费又黄又爽又色| 久久久久久久大尺度免费视频| 看黄色毛片网站| 少妇高潮的动态图| 亚洲精品国产成人久久av| 精品人妻一区二区三区麻豆| 深夜a级毛片| 麻豆成人午夜福利视频| 国产熟女欧美一区二区| 免费观看无遮挡的男女| 亚洲精品一二三| 久久人人爽人人爽人人片va| 极品少妇高潮喷水抽搐| 日韩伦理黄色片| 国产中年淑女户外野战色| av在线蜜桃| 国产在视频线精品| h日本视频在线播放| 啦啦啦中文免费视频观看日本| 亚洲国产日韩一区二区| av国产久精品久网站免费入址| 制服丝袜香蕉在线| 少妇熟女欧美另类| 国产真实伦视频高清在线观看| 九九爱精品视频在线观看| 色吧在线观看| 少妇的逼水好多| 日本av手机在线免费观看| 黄色日韩在线| 久久久久精品性色| 国产真实伦视频高清在线观看| 亚洲不卡免费看| av又黄又爽大尺度在线免费看| 欧美变态另类bdsm刘玥| 亚洲,一卡二卡三卡| 最近中文字幕2019免费版| 亚洲四区av| 国内精品美女久久久久久| 亚洲成人久久爱视频| 国产一区二区三区av在线| 精品99又大又爽又粗少妇毛片| 99久久中文字幕三级久久日本| 97在线人人人人妻| 久久99热这里只有精品18| 午夜福利视频1000在线观看| 亚洲精品日本国产第一区| 国产在线男女| 精品午夜福利在线看| 国产乱人偷精品视频| 97人妻精品一区二区三区麻豆| 久久久久性生活片| 久久久亚洲精品成人影院| 国产精品一区www在线观看| 在线天堂最新版资源| 久久午夜福利片| 国产白丝娇喘喷水9色精品| 久久精品久久精品一区二区三区| 国产精品国产av在线观看| 国产成年人精品一区二区| 亚洲自偷自拍三级| 精品视频人人做人人爽| 久久影院123| 女人久久www免费人成看片| 一级毛片我不卡| 国内少妇人妻偷人精品xxx网站| 久久精品久久久久久久性| 国内揄拍国产精品人妻在线| 国产高清三级在线| 我要看日韩黄色一级片| 久久女婷五月综合色啪小说 | 青春草视频在线免费观看| 久久人人爽av亚洲精品天堂 | 欧美成人一区二区免费高清观看| 成人美女网站在线观看视频| 搡老乐熟女国产| 亚洲精品一二三| 欧美成人a在线观看| 1000部很黄的大片| 嘟嘟电影网在线观看| 丰满人妻一区二区三区视频av| 国产精品国产av在线观看| 色播亚洲综合网| 两个人的视频大全免费| 欧美成人精品欧美一级黄| 日韩三级伦理在线观看| 麻豆成人av视频| 亚洲va在线va天堂va国产| 亚洲av中文av极速乱| 蜜桃久久精品国产亚洲av| 国产精品国产三级国产av玫瑰| 国产伦精品一区二区三区四那| 一区二区三区四区激情视频| 国产精品一区二区性色av| 又爽又黄无遮挡网站| 日韩精品有码人妻一区| 午夜免费观看性视频| 日日摸夜夜添夜夜爱| 亚洲欧美一区二区三区黑人 | 最近的中文字幕免费完整| 99热国产这里只有精品6| 亚洲三级黄色毛片| av在线播放精品| 国产高清不卡午夜福利| 国产老妇伦熟女老妇高清| 成人黄色视频免费在线看| 国产精品麻豆人妻色哟哟久久| 噜噜噜噜噜久久久久久91| 久久久久久久久久人人人人人人| 久久人人爽人人片av| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕久久专区| 欧美zozozo另类| 91久久精品电影网| 亚洲久久久久久中文字幕| 亚洲欧洲日产国产| 亚洲色图av天堂| 亚洲人与动物交配视频| 精品久久久久久电影网| 亚洲欧美成人精品一区二区| 日韩免费高清中文字幕av| 国产精品福利在线免费观看| 国产精品麻豆人妻色哟哟久久| 欧美日韩视频高清一区二区三区二| 欧美三级亚洲精品| 国产一级毛片在线| 一级片'在线观看视频| 在线精品无人区一区二区三 | 国产高清三级在线| 久久久久久久精品精品| 国产成人精品久久久久久| 欧美丝袜亚洲另类| 成年av动漫网址| 欧美bdsm另类| 高清毛片免费看| 男女那种视频在线观看| 777米奇影视久久| 日韩一区二区视频免费看| 激情五月婷婷亚洲| 国产精品不卡视频一区二区| 少妇的逼好多水| 欧美变态另类bdsm刘玥| 日韩成人av中文字幕在线观看| 国产成人精品福利久久| av国产久精品久网站免费入址| 一级毛片黄色毛片免费观看视频| 国语对白做爰xxxⅹ性视频网站| 亚洲av二区三区四区| 能在线免费看毛片的网站| 日本黄大片高清| 国产中年淑女户外野战色| 国产一区二区三区av在线| 高清毛片免费看| 97超视频在线观看视频| 日韩视频在线欧美| 丰满乱子伦码专区| 在线精品无人区一区二区三 | 一级av片app| 一边亲一边摸免费视频| 欧美最新免费一区二区三区| 成年免费大片在线观看| 中文精品一卡2卡3卡4更新| 亚洲最大成人中文| 亚洲人与动物交配视频| av免费在线看不卡| 国产一区二区三区综合在线观看 | 91狼人影院| 好男人在线观看高清免费视频| 最新中文字幕久久久久| 国产精品不卡视频一区二区| 日本欧美国产在线视频| 日韩不卡一区二区三区视频在线| 精品一区二区三卡| 美女xxoo啪啪120秒动态图| 日韩,欧美,国产一区二区三区| 草草在线视频免费看| 在现免费观看毛片| 中文字幕av成人在线电影| 麻豆乱淫一区二区| 视频区图区小说| 乱码一卡2卡4卡精品| 国产一区二区亚洲精品在线观看| 伦精品一区二区三区| 亚洲国产精品成人综合色| 免费观看在线日韩| 久久久久网色| 亚洲真实伦在线观看| 九九爱精品视频在线观看| 777米奇影视久久| 另类亚洲欧美激情| 欧美一区二区亚洲| 免费大片黄手机在线观看| 一区二区三区免费毛片| 欧美日韩视频精品一区| 国产在线男女| 2021少妇久久久久久久久久久| 国内揄拍国产精品人妻在线| av在线亚洲专区| 久久午夜福利片|