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

    大氣環(huán)流模式CAM4對(duì)亞洲氣候的高分辨率模擬

    2016-03-21 03:17:02沙瑩瑩石正國(guó)
    地球環(huán)境學(xué)報(bào) 2016年3期
    關(guān)鍵詞:中亞地區(qū)東亞尺度

    沙瑩瑩,石正國(guó)

    (1.中國(guó)科學(xué)院地球環(huán)境研究所 黃土與第四紀(jì)地質(zhì)國(guó)家重點(diǎn)實(shí)驗(yàn)室,西安 710061;2.中國(guó)科學(xué)院大學(xué),北京 100049;3.中國(guó)科學(xué)院青藏高原地球科學(xué)卓越中心,北京 100101)

    大氣環(huán)流模式CAM4對(duì)亞洲氣候的高分辨率模擬

    沙瑩瑩1,2,石正國(guó)1,3

    (1.中國(guó)科學(xué)院地球環(huán)境研究所 黃土與第四紀(jì)地質(zhì)國(guó)家重點(diǎn)實(shí)驗(yàn)室,西安 710061;2.中國(guó)科學(xué)院大學(xué),北京 100049;3.中國(guó)科學(xué)院青藏高原地球科學(xué)卓越中心,北京 100101)

    中小尺度地形不僅在區(qū)域尺度上影響降水和環(huán)流形勢(shì)分布,還會(huì)影響季風(fēng)、西風(fēng)急流等大尺度的天氣系統(tǒng)。新的全球模式分辨率得到明顯提高,本研究通過對(duì)比分析CAM4高分辨率配置下的模擬結(jié)果與高分辨率APHRO降水和ERA-Interim再分析資料,評(píng)估高分辨率CAM4對(duì)亞洲地區(qū)氣候系統(tǒng)的模擬性能。結(jié)果表明,該模式不僅對(duì)亞洲季風(fēng)-干旱氣候系統(tǒng)的大尺度分布特征有良好的模擬能力,還能很好地模擬出季風(fēng)區(qū)和內(nèi)陸干旱半干旱區(qū)年均和季節(jié)平均降水的高、低值中心以及水平風(fēng)場(chǎng)、地表溫度和海平面氣壓的區(qū)域差異性特征,尤其是山脈地形附近的降水極值和風(fēng)場(chǎng)變化。此外,CAM4與氣象資料間在不同空間尺度上的空間相關(guān)性均很好,其中地表溫度的相關(guān)系數(shù)在0.9以上,降水在亞洲地區(qū)超過0.7。

    CAM4;高分辨率;數(shù)值模擬;亞洲季風(fēng);內(nèi)陸干旱

    亞洲地區(qū)氣候類型復(fù)雜,不僅包括典型的印度季風(fēng)區(qū)和東亞季風(fēng)區(qū),還包括內(nèi)陸干旱和半干旱區(qū)等。其獨(dú)特的氣候系統(tǒng)形成與太陽(yáng)輻射、海陸分布和南北極冰蓋等外部強(qiáng)迫因子密切相關(guān)。青藏高原作為全球最高最大的地形,其對(duì)亞洲氣候的影響也早在20世紀(jì)50年代就引起關(guān)注(Bolin,1950;葉篤正,1952;Flohn,1957)。通過大量的氣候模擬研究,證明青藏高原對(duì)大尺度亞洲季風(fēng)系統(tǒng)的加強(qiáng)和內(nèi)陸干旱化的加劇等均起到不可忽視的作用(Hahn and Manabe,1975;Broccoli and Manabe,1992;Kutzbatch et al,1993;An et al,2001;Liu and Yin,2002;Molnar et al,2010;Shi et al,2011;Wu et al,2012;Jiang et al,2015a)。

    相比于對(duì)青藏高原大地形氣候效應(yīng)的研究,針對(duì)較小地形和山脈氣候效應(yīng)的研究卻很少。這其中一部分原因是由于大部分模式的水平分辨率較低,除青藏高原主體外其他山脈不能很好地在模式中得以體現(xiàn),模式往往忽略了小地形的作用,因此也無法進(jìn)行青藏高原主體外較小山脈地形的氣候效應(yīng)的模擬研究。此外,分辨率低的觀測(cè)資料也往往會(huì)使結(jié)果失真,導(dǎo)致對(duì)較小尺度的氣候特征認(rèn)識(shí)出現(xiàn)偏差(Xie et al,2006;Wu et al,2014)。亞洲地區(qū)的氣候特征受季風(fēng)和西風(fēng)環(huán)流以及復(fù)雜地形的綜合影響,開展亞洲區(qū)域氣候特征及其變化的研究也需要高分辨率的觀測(cè)資料。

    模式水平分辨率的提高能有效改進(jìn)模式對(duì)天氣尺度以及中尺度氣候系統(tǒng)的氣候態(tài)和各氣象要素及其變率的模擬能力(Christensen and Kuhry,2000;Gao et al,2006;Jiang et al,2015b)。通過比較政府間氣候變化專門委員會(huì)(IPCC)評(píng)估報(bào)告的77個(gè)耦合全球氣候模式發(fā)現(xiàn),模式水平分辨率的提高有效增強(qiáng)了其對(duì)中國(guó)地區(qū)年均以及季節(jié)平均溫度和降水氣候態(tài)的模擬能力,此外,對(duì)春季和夏季氣溫以及冬季降水年際變率的模擬能力得以加強(qiáng)(Jiang et al,2015b)。我國(guó)年均降水的雨帶分布也更接近觀測(cè)形態(tài),如低分辨下1000 mm年總降水量等值線位置偏北的狀況得到明顯改善,東亞大尺度降水分布以及東亞夏季風(fēng)期間的較小尺度對(duì)流過程也受到水平分辨率的顯著影響(Gao et al,2006)。

    隨著模式分辨率的提高,一些針對(duì)較小地形的模擬試驗(yàn)表明小地形對(duì)區(qū)域乃至大尺度亞洲季風(fēng)環(huán)流以及內(nèi)陸干旱區(qū)降水等也有重要影響(Xie et al,2006;Wang and Tan,2014;Liu et al,2015)。印度季風(fēng)區(qū)各季風(fēng)降水的對(duì)流中心的形成和所處位置與位于南亞地區(qū)的中尺度山脈(寬度不超過500 km)有關(guān)(Xie et al,2006)。Wu et al(2014) 通過模擬發(fā)現(xiàn),位于中印半島西北部的Arakan山脈對(duì)印度季風(fēng)的爆發(fā)有重要作用,甚至對(duì)東亞夏季風(fēng)的開始也有影響。蒙古高原較青藏高原面積小、高度低,位于偏北的位置,然而有山/無山的敏感性試驗(yàn)表明蒙古高原對(duì)大氣準(zhǔn)定常行星波的形成、東亞西風(fēng)急流的加強(qiáng)以及東亞冬季風(fēng)的加強(qiáng)作用與青藏高原可比(Shi et al,2014;Sha et al,2015)。

    氣候模式分辨率的提高能夠更好地在模式中體現(xiàn)出除青藏高原外的中尺度山脈和更小尺度的地形,有利于模擬出更接近實(shí)際的大氣環(huán)流型和降水特征,從而進(jìn)一步研究小地形的氣候效應(yīng)。目前,利用高分辨率氣候模式模擬亞洲地區(qū)的區(qū)域性氣候特征的研究工作還較少,模式分辨率的提高能否更好地反映出亞洲地區(qū)復(fù)雜的氣候分布特征亟需進(jìn)行評(píng)估。新的全球模式采用有限元差分方法,分辨率也得到了明顯提高。因此,本研究利于大氣環(huán)流模式CAM4高分辨率配置下建立工業(yè)革命前控制試驗(yàn),通過與高分辨率觀測(cè)及再分析資料對(duì)比,評(píng)估CAM4對(duì)亞洲主要?dú)夂蛳到y(tǒng)的模擬能力,并分析亞洲地區(qū)復(fù)雜的氣候特征及可能的影響因素。

    1 數(shù)據(jù)與模式介紹

    降水?dāng)?shù)據(jù)來自日本APHRODITE(Asian Precipitation-Highly Resolved Observational Data Integration Towards Evaluation of Water Resources)建立的高分辨率逐日亞洲陸地降水?dāng)?shù)據(jù)集(簡(jiǎn)稱APHRO),水平分辨率為0.25°×0.25°,時(shí)間從1951年到2007年。該降水?dāng)?shù)據(jù)集能準(zhǔn)確表征我國(guó)雨帶的季節(jié)性移動(dòng),且降水量與臺(tái)站資料一致(韓振宇和周天軍,2012)。因此可以用來研究亞洲地區(qū)降水的氣候特征。由于該降水?dāng)?shù)據(jù)集分為季風(fēng)區(qū)、中亞、俄羅斯和日本四個(gè)子數(shù)據(jù)集,因此在分析前將這四個(gè)子集數(shù)據(jù)進(jìn)行整合,使降水?dāng)?shù)據(jù)覆蓋整個(gè)亞洲地區(qū)。

    海平面氣壓、溫度和風(fēng)場(chǎng)等數(shù)據(jù)來自歐洲中期天氣預(yù)報(bào)中心(ECMWF)的ERA-Interim再分析月均資料,水平分辨率為0.75°×0.75°,時(shí)間從1979年至2014年。該再分析資料能很好地表現(xiàn)海表湍流通量(Brunke et al,2011),其環(huán)流和溫度場(chǎng)等在我國(guó)的適用性也較好(白磊等,2013;曾剛等,2013;高路和郝璐,2014)。為保持降水和再分析資料各變量在時(shí)間上的一致性,本文均選取1979年至2007年共29年進(jìn)行分析。

    本文所用模式為美國(guó)國(guó)家大氣研究中心(NCAR)的全球大氣環(huán)流模式CAM4。CAM4是NCAR的第六代全球大氣環(huán)流模式,相比于CAM3,CAM4缺省的動(dòng)力內(nèi)核由原來的譜內(nèi)核改為有限體積內(nèi)核,CAM4在深對(duì)流方案、北極云量模擬、輻射接口以及計(jì)算可擴(kuò)展性等方面也均做了顯著改進(jìn),有效提高了對(duì)ENSO的模擬能力(Neale et al,2008,2010)。此外,CAM4還包括了陸面模塊CLM4,為大氣模式提供陸面狀況和陸氣之間的能量、動(dòng)量和水汽交換等下邊界條件(Oleson et al,2010)。本研究的模擬試驗(yàn)水平分辨率為0.47°×0.63°(緯向×經(jīng)向),垂直方向上分為26層,模式層頂氣壓為20 hPa,CO2濃度為284.7 mg · kg-1。模擬使用現(xiàn)代氣候平均的海洋表面溫度和海冰范圍作為海洋邊界條件。模式共積分6年,取后5年平均進(jìn)行分析,研究區(qū)域取0° — 60°N,40° — 160°E。

    文章主要分析高分率大氣環(huán)流模式CAM4對(duì)亞洲主要?dú)夂蛳到y(tǒng)的模擬能力,因此需要與高分辨率的觀測(cè)資料進(jìn)行對(duì)比,由于APHRO陸地降水資料的水平分辨率與ERA-Interim再分析資料不同,因此采用雙線性插值方法將APHRO降水資料處理成0.75°×0.75°水平格點(diǎn)資料,插值結(jié)果與原始數(shù)據(jù)在氣候態(tài)上結(jié)果一致。此外,APHRO為逐日資料,而ERA-Interim再分析和CAM4輸出為逐月資料,故而在分析前將ARHRO逐日降水資料處理成月均資料。本文主要進(jìn)行氣候態(tài)分析,夏季平均為6月、7月和8月平均,冬季則為12月、1月和2月。

    2 結(jié)果分析

    2.1 地表溫度

    2.1.1 夏季平均

    熱帶太平洋和印度洋洋面溫度一般在28℃左右。30°N以北的太平洋區(qū)域海表溫度自南向北逐漸遞減(圖1a)。陸地上溫度分布比較復(fù)雜,地表溫度最高的地區(qū)為巴基斯坦中南部,而青藏高原上則為溫度極小值區(qū)。印度季風(fēng)區(qū)溫度梯度較小,地表溫度高值區(qū)位于印度北部和西北部,東亞東部有一暖舌,從20°N向北延伸至東亞東部40°N。我國(guó)新疆有一個(gè)暖中心,中心位于塔里木盆地。中亞地區(qū)地表溫度比明顯高于我國(guó)西北地區(qū),除了位于巴基斯坦的溫度極大值中心外,在里海東部還存在一個(gè)溫度高值區(qū)。由于受海拔高度的影響,青藏高原以及天山周邊的地表溫度梯度很大。

    CAM4模擬的地表溫度梯度最大值區(qū)域同樣位于青藏高原和天山周邊,南側(cè)經(jīng)向溫度梯度為正值而北側(cè)為負(fù)值(圖1b)。東亞東部大陸溫度梯度小,等溫線從南部向北延伸至40°N。我國(guó)新疆南部存在一個(gè)高值區(qū)。中亞地區(qū)的地表溫度較我國(guó)西北地區(qū)高,里海東部有一個(gè)溫度高值中心。整體上看,CAM4模擬的地表溫度在中亞南部以及印度北部地區(qū)偏低,而在東亞東部地區(qū)偏高,但分布形態(tài)與再分析資料一致。

    2.1.2 冬季平均

    冬季海表溫度比陸地高,青藏高原上較同緯度低海拔區(qū)域的溫度低,形成一個(gè)地表溫度低值中心,東亞大陸東部等溫線呈西南—東北向傾斜,地表溫度自南向北逐漸減小,溫度梯度比夏季大(圖2a)。青藏高原仍為地表溫度低值區(qū),但在蒙古國(guó)西北部也出現(xiàn)一個(gè)低值中心,該位置與西伯利亞高壓的中心位置保持一致。除此之外,新疆南部為另一個(gè)地表溫度的低值區(qū),溫度低于-10℃,可見,新疆南部冬季和夏季地表溫度變化劇烈。然而,40°N以南同緯度的中亞地區(qū)地表溫度卻高于0℃,里海為一地表溫度高值中心,溫度超過4℃。冬季地表溫度梯度最大的區(qū)域與夏季有所不同,除青藏高原南側(cè)的極大值區(qū)域外,東亞40°N左右,120° —140°E也為一個(gè)地表溫度梯度高值區(qū),而由于冬季南疆地表溫度的迅速減小,青藏高原北側(cè)的經(jīng)向溫度梯度比夏季小(絕對(duì)值)。

    圖1 夏季平均地表溫度/ ℃(a)ERA-Interim再分析資料夏季平均;(b)CAM4模式結(jié)果Fig.1 Surface temperature /℃ in summer (a) ERA-Interim reanalysis; (b) CAM4 results

    圖2b為CAM4模擬的冬季地表溫度,整體上,海表溫度高于陸地溫度,而中亞地區(qū)溫度高于東亞地區(qū)??梢钥吹剑J讲粌H很好地模擬出了大尺度地表溫度分布形態(tài),對(duì)地表溫度低值和高值中心,以及溫度梯度高值區(qū)也有很好的模擬能力。

    2.2 海平面氣壓

    2.2.1 夏季平均

    夏季時(shí)亞洲大陸上的氣壓較低,在大尺度上沒有明顯的氣壓中心,中印半島以及整個(gè)東亞地區(qū)海平面氣壓基本一樣,等壓線稀疏(圖3a)。太平洋上存在一個(gè)閉合高壓,中心位于30° — 40°N的東太平洋,該高壓中心一直向西延伸至西太平洋。印度洋上海平面氣壓呈東西向水平分布,從南向北氣壓減小,中亞南部有一個(gè)閉合的低壓中心,向北氣壓又逐步增大。CAM4模擬的太平洋副高位置與ERA-Interim再分析資料一致,但強(qiáng)度略大(圖3b)。此外,也模擬出印度洋及中亞地區(qū)海平面的東西向分布特征以及中亞南部的低壓帶。

    2.2.2 冬季平均

    冬季陸地氣壓整體上高于海洋,亞洲大陸和太平洋上均存在兩個(gè)明顯的閉合氣壓中心。西伯利亞高壓中心位于貝加爾湖和巴爾喀什湖之間,該大陸性冷高壓在冬季影響整個(gè)東亞大陸,從西北向東亞東南部方向氣壓逐漸減?。▓D4a)。阿留申低壓位于30°N以北的北太平洋區(qū)域。20°N以南的熱帶地區(qū)氣壓則呈東西向水平分布??梢钥吹?,CAM4模擬的西伯利亞高壓和阿留申低壓位置與ERA-Interim再分析資料一致(圖4b),但強(qiáng)度偏強(qiáng),這可能與再分析資料以及模式的準(zhǔn)確性有關(guān)。

    圖2 冬季平均地表溫度 /℃(a)ERA-Interim再分析資料冬季平均;(b)CAM4模式結(jié)果Fig.2 Surface temperature /℃ in winter (a) ERA-Interim reanalysis; (b) CAM4 results

    2.3 降水和850 hPa風(fēng)場(chǎng)

    2.3.1 年平均特征

    從APHRO降水觀測(cè)資料來看(圖5a),印度和東亞季風(fēng)區(qū)降水最多,印度半島西海岸、印度東北部靠近喜馬拉雅山脈地區(qū)、孟加拉灣北部以及中印半島西部海岸均是印度季風(fēng)降水的大值區(qū),而印度西北部和東南部、中亞半島中部降水相對(duì)較少,東亞季風(fēng)區(qū)降水的大值區(qū)位于我國(guó)東南部沿海,從東南向西北逐漸減少,日本和朝鮮半島降水也較多。對(duì)于內(nèi)陸干旱半干旱區(qū),降水的極小值中心位于我國(guó)新疆以及伊朗、巴基斯坦和阿富汗三國(guó)的交界處,年均降水量在80 mm以下。值得注意的是,位于我國(guó)西部邊界的哈薩克斯坦等國(guó)的東部地區(qū)以及伊朗西北部則為干旱區(qū)降水最多的區(qū)域,年降水量超過300 mm。

    CAM4模擬的印度和東亞季風(fēng)區(qū)降水與APHRO觀測(cè)結(jié)果高度一致(圖5b)。幾個(gè)降水高值區(qū)和低值區(qū)均有清晰的體現(xiàn),如印度半島西海岸和印度東北部靠近喜馬拉雅山脈的降水極大值中心。對(duì)于內(nèi)陸干旱半干旱區(qū),CAM4模擬能力也很好。尤其是新疆地區(qū)的降水極小值中心以及我國(guó)西部邊界各國(guó)東部相對(duì)較高的降水量區(qū)域。

    由此可見,高分辨率的大氣環(huán)流模式不僅能夠模擬出大尺度的季風(fēng)和干旱半干旱區(qū)的大尺度降水特征,也能夠模擬出氣候系統(tǒng)內(nèi)部的差異性特征。這對(duì)于深入認(rèn)識(shí)各氣候系統(tǒng)的特征及其發(fā)展和形成過程及機(jī)理起到至關(guān)重要的作用。

    2.3.2 夏季平均

    與年平均降水相比,印度和東亞季風(fēng)區(qū)降水增加明顯,印度半島西海岸和印度東北部靠近喜馬拉雅山脈的尼泊爾等地區(qū)、孟加拉國(guó)以及中印半島西海岸為降水高值區(qū)(圖6a)。我國(guó)夏季降水從東南沿海向西北內(nèi)陸逐漸減少。對(duì)于內(nèi)陸干旱半干旱區(qū)而言,包括新疆在內(nèi)的我國(guó)西北地區(qū)夏季降水均較年平均多,但除巴基斯坦以外的其他中亞國(guó)家夏季降水則較年平均少。也就是說,夏季對(duì)我國(guó)西北內(nèi)陸來說是雨季而對(duì)中亞地區(qū)來說是旱季(除巴基斯坦)。降水極小值區(qū)出現(xiàn)在里海以東、巴基斯坦和我國(guó)新疆以西的中亞大部分地區(qū),而新疆塔克拉瑪干沙漠地區(qū)降水增多。CAM4模擬的夏季降水分布與APHRO結(jié)果一致(圖6b)。模擬出印度西海岸和靠近喜馬拉雅山脈的尼泊爾等地區(qū)、中印半島西部等降水高值區(qū),以及夏季位于里海以東中亞地區(qū)的降水極小值區(qū)域,我國(guó)西北地區(qū)夏季降水量高于年平均,新疆塔克拉瑪干沙漠地區(qū)降水增多,年平均所在的極小值中心也明顯縮小。

    圖3 夏季海平面氣壓 / hPa(a)ERA-Interim再分析資料夏季平均;(b)CAM4模式結(jié)果Fig.3 Sea level pressure / hPa in summer (a) ERA-Interim reanalysis; (b) CAM4 results

    印度、東亞季風(fēng)區(qū)與干旱半干旱區(qū)的低層風(fēng)場(chǎng)差異明顯(圖6a)。來自印度洋的西南氣流與阿拉伯海的西風(fēng)氣流一起向東,巴基斯坦和印度西北部為西南風(fēng),而印度半島主要為偏西風(fēng),偏西氣流經(jīng)過印度半島分為兩支,一支北上后流向印度東北部、尼泊爾地區(qū),另一支與孟加拉灣的西南氣流一起繼續(xù)向東、向北流去,孟加拉國(guó)為強(qiáng)的南西南風(fēng),西南風(fēng)流經(jīng)中印半島和我國(guó)南海繼續(xù)向北到達(dá)東亞地區(qū),與來自西太平洋的東南風(fēng)共同影響包括我國(guó)東部、朝鮮半島以及日本在內(nèi)的東亞季風(fēng)區(qū)。對(duì)東亞季風(fēng)區(qū)而言,強(qiáng)的偏南風(fēng)從大陸東部進(jìn)入西部?jī)?nèi)陸地區(qū)并逐漸減弱,至100°E左右基本不受季風(fēng)影響,100°E以西主要受西風(fēng)氣流控制。西北氣流沿阿爾泰山脈進(jìn)入我國(guó)新疆,一部分氣流繼續(xù)向東、向南,另一部分氣流則轉(zhuǎn)為東、東北風(fēng)進(jìn)入塔里木盆地,這一風(fēng)場(chǎng)的轉(zhuǎn)向可能和天山以及昆侖山脈的阻擋有關(guān)。中亞地區(qū)的風(fēng)場(chǎng)結(jié)構(gòu)與我國(guó)西北地區(qū)不同,除巴基斯坦外均為很強(qiáng)的偏北風(fēng)。CAM4模擬的低層風(fēng)場(chǎng)與ERA-Interim再分析資料一致(圖6b),尤其是印度季風(fēng)區(qū)的西南風(fēng)和西風(fēng)以及東亞季風(fēng)區(qū)的偏南風(fēng)。此外,模式也很好地模擬出我國(guó)新疆地區(qū)風(fēng)場(chǎng)的轉(zhuǎn)向特征。中亞地區(qū)夏季的偏北風(fēng)也有較好的體現(xiàn),僅低估了35°N以南阿富汗、伊朗地區(qū)的偏北風(fēng)。

    圖4 冬季海平面氣壓/ hPa(a)ERA-Interim再分析資料冬季平均;(b)CAM4模式結(jié)果Fig.4 Sea level pressure / hPa in winter (a) ERA-Interim reanalysis; (b) CAM4 results

    2.3.3 冬季平均

    季風(fēng)區(qū)冬季降水明顯少于年平均,極小值區(qū)域包括印度半島西部、中印半島西部,而我國(guó)30°N以南的東南部地區(qū)和印度北部地區(qū)則為季風(fēng)區(qū)冬季降水高值區(qū)域(圖7a)。我國(guó)整個(gè)西北地區(qū)冬季降水也少于年平均值,而除巴基斯坦外的中亞國(guó)家冬季降水則高于年平均,降水高值區(qū)靠近我國(guó)西部邊界。即:冬季為中亞干旱區(qū)的雨季而為季風(fēng)區(qū)和我國(guó)干旱半干旱區(qū)的旱季。CAM4不僅模擬出季風(fēng)區(qū)和干旱半干旱區(qū)總體的降水分布形態(tài),也對(duì)降水極值中心有很好的反映(圖7b)。如位于印度半島北部和西部的降水極大、極小值區(qū),我國(guó)東南部降水高值區(qū),以及我國(guó)西北部降水極小值區(qū)和中亞東部的降水高值區(qū)等。

    圖5 年平均降水(陰影,mm · day-1)和850 hPa風(fēng)場(chǎng)(矢量箭頭,m · s-1)(a)APHRO降水資料與ERA-Interim再分析資料年平均;(b)CAM4模式結(jié)果Fig.5 Annual mean precipitation (shaded, mm · day-1) and 850 hPa wind (vectors, m · s-1) (a) APHRO precipitation data and ERA-Interim reanalysis; (b) CAM4 results

    冬季,30°N以北的整個(gè)東亞冬季風(fēng)區(qū)為典型的西北風(fēng),以我國(guó)東北、華北地區(qū)為最強(qiáng),沿西南方向逐漸減弱,弱的西、西北風(fēng)與西南風(fēng)在30°N以南的東亞東南部地區(qū)匯合(圖7a)。東亞100°E以西主要為西風(fēng)氣流,但我國(guó)新疆地區(qū)風(fēng)速很小。中亞地區(qū)緯向西風(fēng)氣流比我國(guó)西北地區(qū)強(qiáng),此外,35°N以北區(qū)域還受較強(qiáng)的偏南風(fēng)影響。來自中亞30°N附近的西風(fēng)進(jìn)入印度半島,在其北部和東北部為西北風(fēng)流入孟加拉灣,一支流向緬甸進(jìn)入我國(guó)云南,另一支在20°N以南轉(zhuǎn)向,與來自我國(guó)南海以及中印半島南部的東風(fēng)匯合,流向阿拉伯海以及印度洋。圖7b為CAM4模擬的冬季低層風(fēng)場(chǎng),模式能很好地模擬東亞冬季風(fēng)區(qū)的強(qiáng)西北風(fēng)、中亞35°N以北的西南風(fēng)、以及印度半島和阿拉伯海的偏東風(fēng)。對(duì)我國(guó)西北新疆地區(qū)的風(fēng)場(chǎng)形態(tài)也有良好的模擬能力。但整體上看,模擬的東亞地區(qū)西北風(fēng)和印度季風(fēng)區(qū)的東風(fēng)的強(qiáng)度均偏強(qiáng),此外,中亞南部為偏北風(fēng),這與ERA-Interim在該區(qū)域的風(fēng)場(chǎng)特征存在偏差。

    總體來說,CAM4對(duì)低層風(fēng)場(chǎng)的模擬能力很好,不僅在季風(fēng)-干旱大尺度上表現(xiàn)出很好的一致性,還能準(zhǔn)確模擬出季風(fēng)、干旱系統(tǒng)內(nèi)部的差異。

    2.4 空間相關(guān)系數(shù)

    空間相關(guān)系數(shù)可以進(jìn)一步定量化評(píng)估CAM4對(duì)溫度和降水等近地面氣象要素空間分布的模擬能力。為此,我們同樣利用雙線性插值方法將CAM4模擬結(jié)果插值到0.75°×0.75°的水平格點(diǎn)上,從而與ERA-Interim再分析資料以及插值后的APHRO降水資料在水平格點(diǎn)上保持一致。表1是CAM4與氣象資料在夏季和冬季時(shí)的地表溫度、海平面氣壓、850 hPa緯向風(fēng)和經(jīng)向風(fēng)以及降水在全球、北半球、亞洲、印度季風(fēng)區(qū)、東亞季風(fēng)區(qū)、我國(guó)西北以及中亞地區(qū)的空間相關(guān)系數(shù),其中,也計(jì)算了年均降水在各區(qū)域的空間相關(guān)系數(shù)。這里,亞洲地區(qū)指研究區(qū)域0° — 60°N,40° — 160°E,此外,還包括四個(gè)子區(qū)域:印度季風(fēng)區(qū)為10° — 30°N,70° — 90°E,東亞季風(fēng)區(qū)為20° — 45°N,110° — 140°E,我國(guó)西北地區(qū)為35° — 45°N,80° — 105°E,中亞地區(qū)為30° — 45°N,50° — 75°E。

    圖6 夏季平均降水(陰影,mm · day-1)和850 hPa風(fēng)場(chǎng)(矢量箭頭,m · s-1)(a)APHRO降水資料與ERA-Interim再分析資料夏季平均;(b)CAM4模式結(jié)果Fig.6 precipitation (shaded, mm · day-1) and 850 hPa wind (vectors, m · s-1) in summer (a) APHRO precipitation data and ERA-Interim reanalysis; (b) CAM4 results

    可以看到,CAM4對(duì)各變量空間分布形態(tài)均有很好的模擬能力,均通過99%的顯著性檢驗(yàn)。整體而言,各變量在全球、北半球和亞洲區(qū)域尺度上的空間相關(guān)系數(shù)均在0.7以上。CAM4對(duì)地表溫度的模擬能力最好,與觀測(cè)資料的相關(guān)系數(shù)在0.993 — 0.892,其次為海平面氣壓和850 hPa緯向風(fēng)(二者相當(dāng))、850 hPa經(jīng)向風(fēng),最后為降水場(chǎng)。對(duì)比夏季和冬季的空間相關(guān)系數(shù),地表溫度冬季時(shí)在除我國(guó)西北以外的各區(qū)域其相關(guān)系數(shù)均比夏季高,夏季時(shí)四個(gè)子區(qū)域相關(guān)系數(shù)由高到低分別為:我國(guó)西北、印度地區(qū)、東亞季風(fēng)區(qū)和中亞地區(qū),夏季則為:東亞季風(fēng)區(qū)、印度地區(qū)、中亞和我國(guó)西北。各區(qū)域海平面氣壓夏季的相關(guān)性均比冬季高,夏季時(shí)相關(guān)系數(shù)由高到低分別為:印度、我國(guó)西北、中亞和東亞季風(fēng)區(qū),冬季為印度、東亞季風(fēng)區(qū)、我國(guó)西北和中亞地區(qū)。冬季的850 hPa緯向風(fēng)在全球、北半球和包括東亞季風(fēng)區(qū)的亞洲地區(qū)相關(guān)性比夏季高,而印度季風(fēng)區(qū)、我國(guó)西北和中亞地區(qū)則夏季相關(guān)性高于冬季,其中四個(gè)子區(qū)域在夏季的相關(guān)系數(shù)排序?yàn)椋ǜ咧档降椭担河《?、我?guó)西北、中亞和東亞季風(fēng)區(qū),冬季為:東亞季風(fēng)區(qū)、印度、我國(guó)西北和中亞地區(qū)。850 hPa經(jīng)向風(fēng)則表現(xiàn)為除東亞季風(fēng)區(qū)和我國(guó)西北外夏季的相關(guān)系數(shù)均比冬季高,夏季四個(gè)區(qū)域的相關(guān)系數(shù)順序?yàn)椋ǜ咧档降椭担河《取⒅衼?、我?guó)西北和東亞地區(qū),其中東亞地區(qū)的相關(guān)系數(shù)為0.495,而其他地區(qū)則均超過0.78,冬季順序?yàn)椋褐衼?、我?guó)西北、東亞季風(fēng)區(qū)和印度地區(qū)。值得注意的是,降水空間相關(guān)系數(shù)的季節(jié)性特征與850 hPa緯向風(fēng)相同,即:印度季風(fēng)區(qū)、我國(guó)西北和中亞地區(qū)夏季高于冬季,其他區(qū)域則冬季高于夏季。夏季相關(guān)系數(shù)由高到低分別為:中亞、我國(guó)西北、印度地區(qū)和東亞季風(fēng)區(qū),冬季則為:東亞季風(fēng)區(qū)、印度、中亞和我國(guó)西北地區(qū)。就年均降水的區(qū)域分布來看,東亞地區(qū)的空間相關(guān)系數(shù)最高,達(dá)到0.902,其次分別為中亞地區(qū)、我國(guó)西北和印度地區(qū)。

    圖7 冬季平均降水(陰影,mm · day-1)和850 hPa風(fēng)場(chǎng)(矢量箭頭,m · s-1)(a)APHRO降水資料與ERA-Interim再分析資料冬季平均;(b)CAM4模式結(jié)果Fig.7 precipitation (shaded, mm · day-1) and 850 hPa wind (vectors, mm · day-1) in winter (a) APHRO precipitation data and ERA-Interim reanalysis; (b) CAM4 results

    表1 現(xiàn)代氣象資料與CAM4氣候態(tài)的空間相關(guān)系數(shù)Tab.1 Spatial correlation between meteorological and CAM4 data

    可以看到,CAM4對(duì)各變量空間分布形態(tài)均有很好的模擬能力,均通過99%的顯著性檢驗(yàn)。整體而言,各變量在全球、北半球和亞洲區(qū)域尺度上的空間相關(guān)系數(shù)均在0.7以上。CAM4對(duì)地表溫度的模擬能力最好,與觀測(cè)資料的相關(guān)系數(shù)在0.993 — 0.892,其次為海平面氣壓和850 hPa緯向風(fēng)(二者相當(dāng))、850 hPa經(jīng)向風(fēng),最后為降水場(chǎng)。對(duì)比夏季和冬季的空間相關(guān)系數(shù),地表溫度冬季時(shí)在除我國(guó)西北以外的各區(qū)域其相關(guān)系數(shù)均比夏季高,夏季時(shí)四個(gè)子區(qū)域相關(guān)系數(shù)由高到低分別為:我國(guó)西北、印度地區(qū)、東亞季風(fēng)區(qū)和中亞地區(qū),夏季則為:東亞季風(fēng)區(qū)、印度地區(qū)、中亞和我國(guó)西北。各區(qū)域海平面氣壓夏季的相關(guān)性均比冬季高,夏季時(shí)相關(guān)系數(shù)由高到低分別為:印度、我國(guó)西北、中亞和東亞季風(fēng)區(qū),冬季為印度、東亞季風(fēng)區(qū)、我國(guó)西北和中亞地區(qū)。冬季的850 hPa緯向風(fēng)在全球、北半球和包括東亞季風(fēng)區(qū)的亞洲地區(qū)相關(guān)性比夏季高,而印度季風(fēng)區(qū)、我國(guó)西北和中亞地區(qū)則夏季相關(guān)性高于冬季,其中四個(gè)子區(qū)域在夏季的相關(guān)系數(shù)排序?yàn)椋ǜ咧档降椭担河《?、我?guó)西北、中亞和東亞季風(fēng)區(qū),冬季為:東亞季風(fēng)區(qū)、印度、我國(guó)西北和中亞地區(qū)。850 hPa經(jīng)向風(fēng)則表現(xiàn)為,除東亞季風(fēng)區(qū)和我國(guó)西北外夏季的相關(guān)系數(shù)均比冬季高,夏季四個(gè)區(qū)域的相關(guān)系數(shù)順序?yàn)椋ǜ咧档降椭担河《取⒅衼?、我?guó)西北和東亞地區(qū),其中東亞地區(qū)的相關(guān)系數(shù)為0.495,而其他地區(qū)則均超過0.78,冬季順序?yàn)椋褐衼?、我?guó)西北、東亞季風(fēng)區(qū)和印度地區(qū)。值得注意的是,降水空間相關(guān)系數(shù)的季節(jié)性特征與850 hPa緯向風(fēng)相同,即:印度季風(fēng)區(qū)、我國(guó)西北和中亞地區(qū)夏季高于冬季,其他區(qū)域則冬季高于夏季,夏季相關(guān)系數(shù)由高到低分別為:中亞、我國(guó)西北、印度地區(qū)和東亞季風(fēng)區(qū),冬季則為:東亞季風(fēng)區(qū)、印度、中亞和我國(guó)西北地區(qū)。就年均降水的區(qū)域分布來看,東亞地區(qū)的空間相關(guān)系數(shù)最高,達(dá)到0.902,其次分別為中亞地區(qū)、我國(guó)西北和印度地區(qū)。

    3 結(jié)論

    印度和東亞季風(fēng),以及內(nèi)陸干旱和半干旱的氣候分布是亞洲地區(qū)最基本的大尺度氣候特征,這一大尺度特征在很多分辨率較低的大氣/氣候模式中也有良好的模擬(Broccoli and Manabe,1992;Kutzbatch et al,1993;Kitoh,2004),然而,低分辨率模式模擬的結(jié)果卻無法準(zhǔn)確反映降水和氣溫等物理量的高、低值中心以及環(huán)流和氣壓場(chǎng)的梯度分布等較小尺度特征,這不僅會(huì)導(dǎo)致對(duì)局地氣候的認(rèn)識(shí)出現(xiàn)偏差,也導(dǎo)致可能忽略小尺度變化對(duì)大尺度季風(fēng)-干旱環(huán)境的重要作用。高分辨率的大氣環(huán)流模式CAM4不僅能模擬出季風(fēng)-干旱分布的大尺度降水、流場(chǎng)和溫度場(chǎng)等特征,還能很好地模擬出較小尺度的氣候差異。主要結(jié)論如下:

    (1)從降水分布來看,印度半島西海岸、喜馬拉雅山南部、孟加拉國(guó)和中印半島西部等亞洲季風(fēng)降水的高值中心在CAM4中得以良好的體現(xiàn),此外,模式也能反映出夏季中亞降水低于我國(guó)內(nèi)陸而冬季多于我國(guó)內(nèi)陸的亞洲內(nèi)陸地區(qū)的差異性特征。

    (2)高分辨率CAM4對(duì)季節(jié)性流場(chǎng)的模擬更接近實(shí)際分布,如印度季風(fēng)區(qū)夏季低層風(fēng)場(chǎng)在印度洋、印度半島和孟加拉灣及其北部陸地以及中印半島的區(qū)域性差異,以及東亞大陸東部的西南風(fēng)逐漸向西北轉(zhuǎn)為東南風(fēng)的特征,尤其對(duì)夏季新疆地區(qū)風(fēng)場(chǎng)轉(zhuǎn)向的模擬,此外,也體現(xiàn)出中亞地區(qū)夏季為偏北風(fēng)而冬季主要為偏南風(fēng)的季節(jié)性差異。

    (3)從溫度場(chǎng)和氣壓場(chǎng)來看,除海陸間的大尺度差異以及冬季的西伯利亞高壓和夏季的副熱帶高壓等大尺度系統(tǒng)外,CAM4也能模擬出季節(jié)性地表溫度高、低值中心以及等溫線密集區(qū)(溫度梯度高值區(qū))。

    (4)高分辨率CAM4的近地表各氣象要素的空間分布特征與氣象觀測(cè)資料間相關(guān)關(guān)系很高,其中地表溫度的空間相關(guān)系數(shù)最高,在不同區(qū)域上均能達(dá)到0.9以上,海平面氣壓和850 hPa緯向風(fēng)其次,850 hPa經(jīng)向風(fēng)比緯向風(fēng)低,降水最低。

    從高分辨率APHRO降水資料和ERA-Interim再分析資料可以看到,無論是季風(fēng)降水還是干旱區(qū)降水,降水的高值區(qū)往往位于山脈迎風(fēng)面而低值區(qū)位于山脈背風(fēng)面,這一特征與印度夏季風(fēng)的對(duì)流活動(dòng)中心位置和對(duì)應(yīng)山脈地形間的關(guān)系一致(Xie et al,2006;Wu et al,2014)。Gao at al(2006)通過對(duì)比不同地形分辨率的模擬試驗(yàn)后也指出,隨著地形分辨率的提高,模擬的降水分布更接近實(shí)際。除此之外,水平風(fēng)場(chǎng)在山脈附近會(huì)發(fā)生明顯的變化(李斐等,2012;Li and Zhang,2012),從而通過溫度平流等作用直接影響到溫度分布,進(jìn)一步改變環(huán)流場(chǎng)的結(jié)構(gòu)和垂直運(yùn)動(dòng),最終對(duì)大尺度氣候系統(tǒng)產(chǎn)生影響(Shi et al,2014;Sha et al,2015)。因此,深入認(rèn)識(shí)亞洲季風(fēng)-干旱氣候系統(tǒng)須采用高分辨率氣候模式,分析區(qū)域差異產(chǎn)生的原因以及對(duì)大尺度氣候系統(tǒng)的作用,從而避免對(duì)模式結(jié)果的理解產(chǎn)生偏差。降水及環(huán)流特征的區(qū)域差異在高分辨率CAM4的模擬結(jié)果中得以體現(xiàn),初步來看,近地面風(fēng)在遇到山脈阻擋后往往會(huì)發(fā)生繞流和爬坡,因此有利于山前降水,此外,也可能改變水汽的輸送路徑和水汽的輻合輻散,但降水和環(huán)流場(chǎng)等氣象要素和氣候系統(tǒng)與山脈地形的關(guān)系需要進(jìn)一步進(jìn)行敏感性試驗(yàn)以及對(duì)觀測(cè)資料的深入分析。

    白 磊, 王維霞, 姚亞楠, 等. 2013. ERA-Interim和NCEP/ NCAR 再分析數(shù)據(jù)氣溫和氣壓值在天山山區(qū)適用性分析 [J].沙漠與綠洲氣象, 7(3): 51 – 56. [Bai L, Wang W X, Yao Y N, et al. 2013. Reliability of NCEP/NCAR and ERA-Interim reanalysis data on Tianshan mountainous area [J].Dersert and Oasis Meteorology, 7(3): 51 – 56.]

    曾 剛, 伯忠凱, 倪東鴻, 等. 2013. 多套大氣再分析資料的南亞高壓強(qiáng)度變化特征及其與海表溫度異常關(guān)系的比較分析 [J].大氣科學(xué)學(xué)報(bào), 36(5): 577 – 585. [Zeng G, Bo Z K, Ni D H, et al. 2013. Comparison analysis of South Asia high intensity variation and its relation to SSTA in atmospheric multi-reanalysis data [J].Transactions of Atmospheric Sciences, 36(5): 577 – 585.]

    高 路, 郝 璐. 2014. ERA-Interim 氣溫?cái)?shù)據(jù)在中國(guó)區(qū)域的適用性評(píng)估 [J].亞熱帶資源與環(huán)境學(xué)報(bào), 9(2): 75 – 81. [Gao L, Hao L. 2014. Verifi cation of ERA-Interim reanalysis data over China [J].Journal of Subtropical Resources and Environment, 9(2): 75 – 81.]

    韓振宇, 周天軍. 2012. APHRODITE 高分辨率逐日降水資料在中國(guó)大陸地區(qū)的適用性[J].大氣科學(xué), 36(2): 361 – 373. [Han Z Y, Zhou T J. 2012. Assessing the quality of APHRODITE high-resolution daily precipitation dataset over contiguous China [J].Chinese Journal of Atmospheric Sciences, 36(2): 361 – 373.]

    李 斐, 李建平, 李艷杰, 等. 2012. 青藏高原繞流和爬流的氣候?qū)W特征 [J].大氣科學(xué), 36(6): 1236 – 1252. [Li F, Li J P, Li Y J, et al. 2012. Climatological characteristics of fl ow around and fl ow over the Tibetan Plateau [J].Chinese Journal of Atmospheric Sciences, 36(6): 1236 – 1252.]

    葉篤正. 1952. 西藏高原對(duì)于大氣環(huán)流影響的季節(jié)變化 [J].氣象學(xué)報(bào), 23: 33 – 47. [Yeh T C. 1952. The seasonal variation of the infl uence of Tibetan Plateau on the general circulation [J].Acta Meteoro Sinica, 23: 33 – 47.]

    An Z S, Kutzbach J E, Prell W L, et al. 2001. Evolution of Asian monsoons and phased uplift of the Himalaya-Tibetan plateau since Late Miocene times [J].Nature, 411(6833): 62 – 66.

    Bolin B. 1950. On the influence of the Earth's orography on the general character of the westerlies [J].Tellus A, 2(3): 184 – 195.

    Broccoli A J, Manabe S. 1992. The effects of orography on midlatitude Northern Hemisphere dry climates [J].Journal of Climate, 5(11): 1181 – 1201.

    Brunke M A, Wang Z, Zeng X B, et al. 2011. An assessment of the uncertainties in ocean surface turbulent fl uxes in 11 reanalysis, satellite-derived, and combined global datasets [J].Journal of Climate, 24(21): 5469 – 5493.

    Christensen J H, Kuhry P. 2000. High-resolution regional climate model validation and permafrost simulation for the East European Russian Arctic [J].Journal of Geophysical Research, 105(D24): 29647 – 29658.

    Flohn H. 1957. Large-scale aspects of the “summer monsoon”in south and east Asia [J].Journal of the Meteorological Society of Japan, 75: 180 – 186.

    Gao X J, Xu Y, Zhao Z C, et al. 2006. On the role of resolution and topography in the simulation of East Asia precipitation [J].Theoretical and Applied Climatology, 86(1 / 2/ 3/ 4): 173 – 185.

    Hahn D G, Manabe S. 1975. The role of mountains in the south Asian monsoon circulation [J].Journal of the Atmospheric Sciences, 32(8): 1515 – 1541.

    Jiang D B, Tian Z P, Lang X M. 2015b. Reliability of climate models for China through the IPCC Third to Fifth Assessment Reports [J].International Journal of Climatology, 36(3): 1114 – 1133.

    Jiang D B, Yu G, Zhao P, et al. 2015a. Paleoclimate modeling in China: A review [J].Advances in Atmospheric Sciences, 32(2): 250 – 275.

    Kitoh A. 2004. Effects of mountain uplift on East Asian summer climate investigated by a coupled atmosphereocean GCM [J].Journal of Climate, 17(4): 783 – 802.

    Kutzbach J E, Prell W L, Ruddiman W F. 1993. Sensitivity of Eurasian climate to surface uplift of the Tibetan Plateau [J].The Journal of Geology, 101(2): 177 – 190.

    Li Q, Zhang R H. 2012. Seasonal variation of climatological bypassing fl ows around the Tibetan Plateau [J].Advances in Atmospheric Sciences, 29: 1100 – 1110.

    Liu X D, Sun H, Miao Y F, et al. 2015. Impacts of uplift of northern Tibetan Plateau and formation of Asian inland deserts on regional climate and environment [J].Quaternary Science Reviews, 116: 1 – 14.

    Liu X D, Yin Z Y. 2002. Sensitivity of East Asian monsoon climate to the uplift of the Tibetan Plateau [J].Palaeogeography, Palaeoclimatology, Palaeoecology, 183(3): 223 – 245.

    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 [J].Annual Review of Earth and Planetary Sciences, 38(1): 77 – 102.

    Neale R B, Richter J H, Jochum M, et al. 2010. Description of the NCAR Community Atmosphere Model (CAM 4.0) [R]. Technical Note (NCAR/TN-486-STR), National Center for Atmospheric Research.

    Neale R B, Richter J H, Jochum, M. 2008. The impact of convection on ENSO: From a delayed oscillator to a series of events [J].Journal of climate, 21(22): 5904 – 5924.

    Oleson K W, Lawrence D M, Gordon B, et al. 2010. Technical description of version 4.0 of the Community Land Model (CLM) [R]. Technical Note NCAR/TN-478+STR, National Center for Atmospheric Research.

    Sha Y Y, Shi Z G, Liu X D, et al. 2015. Distinct impacts of the Mongolian and Tibetan plateaus on the evolution of the East Asian monsoon [J].Journal of Geophysical Research, 120: 4764 – 4782.

    Shi Z G, Liu X D, An Z S, et al. 2011. Simulated variations of eolian dust from inner Asian deserts at the mid-Pliocene, last glacial maximum, and present day: contributions from the regional tectonic uplift and global climate change [J].Climate Dynamics, 37: 2289 – 2301.

    Shi Z G, Liu X D, Liu Y M, et al. 2014. Impact of Mongolian Plateau versus Tibetan Plateau on the westerly jet over North Pacifi c Ocean [J].Climate Dynamics, 44(11 / 12): 3067 – 3076.

    Wang Q W, Tan Z M. 2014. Multi-scale topographic control of southwest vortex formation in Tibetan Plateau region in an idealized simulation [J].Journal of Geophysical Research, 119(20): 11543 – 11561.

    Wu C H, Hsu H H, Chou M D. 2014. Effect of the Arakan Mountains in the northwestern Indochina Peninsula on the late May Asian monsoon transition [J].Journal of Geophysical Research, 119(18): 10769 – 10779.

    Wu G X, Liu Y, He B, et al. 2012. Thermal controls on the Asian summer monsoon [J].Scientific Reports, 2, doi: 10.1038/srep00404.

    Xie S P, Xu H, Saji N H, et al. 2006. Role of narrow mountains in large-scale organization of Asian monsoon convection [J].Journal of Climate, 19(14): 3420 – 3429.

    A high resolution simulation on Asian climate by Community Atmosphere Model 4

    SHA Yingying1,2, SHI Zhengguo1,3
    (1. State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment, Chinese Academy of Sciences, Xi'an 710061, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 3. CAS Center for Excellence in Tibetan Plateau Earth Sciences, Beijing 100101, China)

    Background, aim, and scopeAsian climate is characterized by the Indian monsoon, the East Asian monsoon and inland aridity, which is complex and sensitive to global change. The formation and variation of the climate system are correlated with lots of external forcings, such as solar radiation, landsea distribution and ice sheet. The Tibetan Plateau is also demonstrated to have a vital effect on the Asian monsoons and inland aridity. Besides the large scale plateau, the meso- and small-scale mountains also show great infl uence. These mountains not only affect the regional precipitation and circulation, but also infl uence large scale weather systems such as monsoons and westerly jet fl ow. The climate model CAM4 is with much higher resolution, it can refl ect a more realistic topography and resolute those meso- and small-scale mountains to evaluate the climate effect of them.Materials and methodsThe precipitation data is a daily database from Asian Precipitation-Highly Resolved Observational Data IntegrationTowards Evaluation of Water Resources of Japan. Its spatial resolution is 0.25°×0.25° from 1951 to 2007. Other meteorological variables such as sea level pressure, horizontal wind and temperature is from ERA-Interim reanalysis data of European Center for Medium-Range Weather Forecasts(ECMWF). The reanalysis data is from 1979 to present with a resolution of 0.75°×0.75°. Besides the meteorological data, the climate model we used to simulate is the CAM4. The CAM4 is the sixth generation of the NCAR atmospheric general climate model. The horizontal resolution is 0.47°×0.63° (latitude×longitude) in the simulation.ResultsWith the apparently increased resolution of global climate model, the modeling result conducted by CAM4 with high resolution is compared with APHRO precipitation and ERAInterim reanalysis data to evaluate whether and in what extent the simulating capability is improved. (1) The CAM4 with high resolution can well simulates the large scale Asian monsoon — aridity climate pattern. (2) More importantly, the Asian monsoon precipitation are well presented, including the Indian monsoon rain centers over the west coast of the Indian subcontinent, the southern Himalaya mountains, Bangladesh and the west coast of Indochina peninsula, and the East Asian monsoon rain gradient decreases from southeast to northwest. The CAM4 also shows more precipitation in winter and less precipitation in summer over arid inland of the Central Asia, which is different compared to inland China, where the most raining season is summer. (3) The atmosphere circulation simulated in high-resoluted CAM4 is more realistic, such as the 850 hPa wind difference over the Indian ocean, the Indian subcontinent, the Bay of Bengal and the Indochina peninsula of the Indian monsoon region, the wind direction turning from southwesterly to southeasterly over east China of the East Asian monsoon region, and the summer bypassing fl ows over the south Tianshan mountain in Xinjiang. (4) Besides, the surface temperature and sea level pressure patterns of the Asian monsoon and interior aridity are also well simulated, including the large scale climate systems such as the Siberian high, the Aleutian low and the western Pacifi c high, and the seasonal variations of the meso-scale pressure and temperature centers and gradients. (5) The spatial correlations show that the correlation between CAM4 and meteorological observation are signifi cant on global, North Hemisphere and continent scale. The surface temperature has a very high correlation between high-resoluted CAM4 model and meteorological data. The correlation coeffi cient of surface temperature is above 0.9. The correlation coeffi cient of sea level pressure, 850 hPa zonal wind and meridional wind is above 0.9, 0.9 and 0.8 on global scale, respectively. The correlation coeffi cient of precipitation over Asia is no less than 0.7.DiscussionThe results suggest that the precipitation center is always located over upwind regions no matter over monsoon or inland arid areas while the downwind regions always get minimum rainfall. The character is accordance with the relationship between the convective activity and mountains. Besides, the wind field also changes near mountainous regions and obviously different from other regions. The related temperature distribution can be directly infl uenced by temperature advection change caused by wind variation. As a result, the atmosphere circulation and vertical structure are changed, which can also infl uence large scale climate systems. Lots of studies show that more realistic precipitation and wind fi eld can be simulated with higher topography resolution in climate models. Therefore, high-resoluted climate models need to be used to study the detailed monsoon and arid climate and inner difference over Asia and study the relationship between precipitation, atmosphere circulation and mountains.ConclusionsThe results show that the high-resoluted CAM4 model can improve the simulation skill over Asian, where the climate system is characterized by complex monsoons and aridity. With the improved resolution, the differences between the Indian monsoon and the East Asian monsoon, and between the arid Central Asia and inland China are obviously presented. Besides, it suggests a great infl uence of meso- and small-scale mountains on convection centers and precipitation distribution.Recommendations and perspectivesThe high-resoluted CAM4 model is very helpful for detailed study of the complex Asian climate and to deep our understanding of the infl uence of the mountain topogaraphy. The sensitivity experiments related to meso-scale mountain need to be designed to further fi nger out its climate effect.

    CAM4; high resolution; simulation; Asian monsoon; interior aridity

    SHI Zhengguo, E-mail: shizg@ieecas.cn

    10.7515/JEE201603006

    2015-11-23;錄用日期:2015-12-24

    Received Date:2015-11-23;Accepted Date:2015-12-24

    中國(guó)科學(xué)院戰(zhàn)略重點(diǎn)研究項(xiàng)目(XDB03020601);國(guó)家自然科學(xué)基金項(xiàng)目(41290255,41572160);國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2013CB955904);黃土與第四紀(jì)地質(zhì)國(guó)家重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目(SKLLQG1416)

    Foundation Item:Strategic Priority Research Program of Chinese Academy of Sciences (XDB03020601); National Natural Science Foundation of China (41290255, 41572160); National Basic Research Program of China (2013CB955904); Program of State Key Laboratory of Loess and Quaternary Geology (SKLLQG1416)

    石正國(guó),E-mail: shizg@ieecas.cn

    猜你喜歡
    中亞地區(qū)東亞尺度
    民心相通語(yǔ)境下中華文化在中亞地區(qū)的傳播及其影響
    “東亞漢詩(shī)史(多卷本)”簡(jiǎn)介
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    我校東亞研究院一行應(yīng)邀訪問韓國(guó)東亞大學(xué)
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    未定的東亞多邊體制與能動(dòng)的中國(guó)東亞政策
    淺談中亞地區(qū)EPC項(xiàng)目當(dāng)?shù)毓蛦THSE管理
    當(dāng)前中亞地區(qū)安全局勢(shì)的特點(diǎn)及走勢(shì)分析
    9
    中國(guó)與中亞地區(qū)油氣資源合作開發(fā)模式與前景分析
    国产v大片淫在线免费观看| 此物有八面人人有两片| 免费av毛片视频| 亚洲av第一区精品v没综合| 国产精品一区www在线观看 | 日日撸夜夜添| 啪啪无遮挡十八禁网站| 亚洲人成网站在线播| 国产免费男女视频| 国产一区二区在线观看日韩| 欧美性感艳星| 国产精品伦人一区二区| 国产免费男女视频| 欧美日韩精品成人综合77777| 最近视频中文字幕2019在线8| 亚洲四区av| 久久久久精品国产欧美久久久| 伊人久久精品亚洲午夜| eeuss影院久久| 观看美女的网站| 久久久久性生活片| 波多野结衣巨乳人妻| 麻豆久久精品国产亚洲av| 久久这里只有精品中国| 中国美女看黄片| 中文字幕免费在线视频6| 内射极品少妇av片p| 亚洲真实伦在线观看| 久久久久久久午夜电影| 很黄的视频免费| 日日摸夜夜添夜夜添小说| 麻豆久久精品国产亚洲av| 美女高潮喷水抽搐中文字幕| 国产欧美日韩一区二区精品| 99热精品在线国产| 久久精品91蜜桃| 国产一级毛片七仙女欲春2| 黄色女人牲交| 国产精品免费一区二区三区在线| 日韩高清综合在线| 淫妇啪啪啪对白视频| 国产亚洲精品久久久久久毛片| eeuss影院久久| 欧美成人性av电影在线观看| 欧美高清性xxxxhd video| 亚洲最大成人中文| 婷婷六月久久综合丁香| 一本久久中文字幕| 免费人成在线观看视频色| 国产一区二区三区视频了| 国产精品免费一区二区三区在线| 成人av在线播放网站| 日韩大尺度精品在线看网址| 男人舔奶头视频| 女生性感内裤真人,穿戴方法视频| 国产欧美日韩一区二区精品| 亚洲午夜理论影院| av国产免费在线观看| 观看美女的网站| 男人狂女人下面高潮的视频| 白带黄色成豆腐渣| 一a级毛片在线观看| 亚洲七黄色美女视频| 人妻久久中文字幕网| 看免费成人av毛片| 国产精品免费一区二区三区在线| 国产精品一区www在线观看 | 色精品久久人妻99蜜桃| 午夜精品一区二区三区免费看| 久久精品91蜜桃| 国产精品一区二区三区四区久久| 午夜激情福利司机影院| 熟妇人妻久久中文字幕3abv| 91av网一区二区| 日本爱情动作片www.在线观看 | 国产成人一区二区在线| 久久亚洲精品不卡| 欧美成人a在线观看| 别揉我奶头 嗯啊视频| 一级黄色大片毛片| 99久久中文字幕三级久久日本| www.www免费av| 日韩欧美国产在线观看| 国产高清视频在线播放一区| 黄色视频,在线免费观看| 99热这里只有是精品在线观看| av在线亚洲专区| 午夜福利18| 制服丝袜大香蕉在线| 国产精品伦人一区二区| 九九在线视频观看精品| 亚洲内射少妇av| 无人区码免费观看不卡| 黄色配什么色好看| 在线免费观看的www视频| 男插女下体视频免费在线播放| 日本一本二区三区精品| 麻豆久久精品国产亚洲av| 国产精品三级大全| 国内精品宾馆在线| 日日摸夜夜添夜夜添av毛片 | 黄色丝袜av网址大全| 欧美日韩综合久久久久久 | 精品午夜福利视频在线观看一区| 又紧又爽又黄一区二区| 亚洲精品粉嫩美女一区| 免费看av在线观看网站| 在线a可以看的网站| 亚洲av美国av| 在线看三级毛片| 尾随美女入室| 黄色女人牲交| 日韩 亚洲 欧美在线| 在线国产一区二区在线| avwww免费| 亚洲精华国产精华精| 久久人妻av系列| 国产一区二区三区视频了| 男女边吃奶边做爰视频| 99久国产av精品| 亚洲av二区三区四区| 国产私拍福利视频在线观看| 久久精品国产亚洲网站| 白带黄色成豆腐渣| 久久久久久国产a免费观看| 很黄的视频免费| 亚洲欧美日韩东京热| 久久精品夜夜夜夜夜久久蜜豆| 97热精品久久久久久| 十八禁国产超污无遮挡网站| 午夜福利在线在线| 国模一区二区三区四区视频| 赤兔流量卡办理| 久久人人精品亚洲av| 日日干狠狠操夜夜爽| 亚洲图色成人| 久久精品人妻少妇| 夜夜爽天天搞| 国产精品电影一区二区三区| 少妇丰满av| 99riav亚洲国产免费| 黄色欧美视频在线观看| 国产男靠女视频免费网站| 狠狠狠狠99中文字幕| 亚洲精品乱码久久久v下载方式| 国产亚洲精品久久久久久毛片| 床上黄色一级片| 99热这里只有是精品在线观看| 成人国产综合亚洲| 免费高清视频大片| 99久国产av精品| 中文亚洲av片在线观看爽| 九色成人免费人妻av| 两人在一起打扑克的视频| 色综合色国产| 国产黄a三级三级三级人| xxxwww97欧美| 亚洲中文字幕日韩| 国内毛片毛片毛片毛片毛片| 观看免费一级毛片| 亚洲精品一区av在线观看| av视频在线观看入口| 午夜久久久久精精品| 欧美一区二区精品小视频在线| av在线天堂中文字幕| 欧美另类亚洲清纯唯美| 日本黄大片高清| 亚洲av日韩精品久久久久久密| 99久国产av精品| 久久精品91蜜桃| 我要看日韩黄色一级片| 老女人水多毛片| 男人舔奶头视频| 久久6这里有精品| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产| 极品教师在线免费播放| 国产色爽女视频免费观看| 日韩一本色道免费dvd| 一夜夜www| 国产综合懂色| 一级a爱片免费观看的视频| 18禁裸乳无遮挡免费网站照片| av专区在线播放| 日本色播在线视频| 久久久国产成人免费| 中文资源天堂在线| 国产精品一区二区三区四区久久| 99热这里只有精品一区| 国产精品一区二区三区四区免费观看 | 小说图片视频综合网站| 日本黄大片高清| 精华霜和精华液先用哪个| 国产精品女同一区二区软件 | 夜夜看夜夜爽夜夜摸| 亚洲精品日韩av片在线观看| 综合色av麻豆| 精品无人区乱码1区二区| 精品一区二区三区av网在线观看| 黄色一级大片看看| 人人妻人人澡欧美一区二区| 国产国拍精品亚洲av在线观看| 亚洲精品久久国产高清桃花| 一卡2卡三卡四卡精品乱码亚洲| 伊人久久精品亚洲午夜| 久久欧美精品欧美久久欧美| 精品久久国产蜜桃| 国产久久久一区二区三区| 午夜福利在线观看免费完整高清在 | 麻豆成人av在线观看| 哪里可以看免费的av片| 内地一区二区视频在线| 嫩草影视91久久| 中文字幕熟女人妻在线| 久久久久性生活片| 国产精品伦人一区二区| 别揉我奶头 嗯啊视频| 在线观看舔阴道视频| 日本成人三级电影网站| 搡女人真爽免费视频火全软件 | 我的女老师完整版在线观看| 一个人观看的视频www高清免费观看| 又爽又黄无遮挡网站| 韩国av一区二区三区四区| 日日啪夜夜撸| 亚洲av日韩精品久久久久久密| 免费在线观看成人毛片| 国产伦在线观看视频一区| 亚洲美女黄片视频| 给我免费播放毛片高清在线观看| 国产午夜精品论理片| 深夜精品福利| 久久久久久久久久久丰满 | 久久天躁狠狠躁夜夜2o2o| 日韩,欧美,国产一区二区三区 | 成年女人看的毛片在线观看| 婷婷六月久久综合丁香| 亚洲aⅴ乱码一区二区在线播放| 蜜桃亚洲精品一区二区三区| 欧美人与善性xxx| 精品人妻1区二区| 久久久久久久午夜电影| 亚洲人与动物交配视频| 国产精品野战在线观看| 久久久久久九九精品二区国产| 九色成人免费人妻av| 韩国av在线不卡| 欧美绝顶高潮抽搐喷水| 在线观看一区二区三区| 国产中年淑女户外野战色| av天堂中文字幕网| 色播亚洲综合网| 中文字幕av在线有码专区| 色综合亚洲欧美另类图片| 两人在一起打扑克的视频| 国产午夜福利久久久久久| 女人十人毛片免费观看3o分钟| a级一级毛片免费在线观看| 五月伊人婷婷丁香| 成人鲁丝片一二三区免费| 人妻夜夜爽99麻豆av| 99久久无色码亚洲精品果冻| 亚洲va日本ⅴa欧美va伊人久久| 色哟哟·www| 婷婷精品国产亚洲av在线| 亚洲va在线va天堂va国产| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久国产a免费观看| 中文字幕人妻熟人妻熟丝袜美| 久久人人精品亚洲av| 亚洲国产色片| 动漫黄色视频在线观看| 国产精品久久视频播放| 99久久九九国产精品国产免费| 琪琪午夜伦伦电影理论片6080| 老师上课跳d突然被开到最大视频| 亚洲国产精品sss在线观看| 国产真实乱freesex| 国产精品98久久久久久宅男小说| 69人妻影院| 波多野结衣巨乳人妻| 国产精品人妻久久久影院| 精品久久久久久成人av| 久久精品国产自在天天线| 韩国av在线不卡| 亚洲成人久久性| 欧美成人一区二区免费高清观看| 成人国产综合亚洲| av视频在线观看入口| 在线播放无遮挡| 少妇的逼好多水| 国产av一区在线观看免费| 美女免费视频网站| 老司机午夜福利在线观看视频| 久久精品人妻少妇| 国产精品美女特级片免费视频播放器| 亚洲成人久久性| 午夜免费激情av| h日本视频在线播放| 99九九线精品视频在线观看视频| 免费一级毛片在线播放高清视频| АⅤ资源中文在线天堂| 精品久久久久久久久久免费视频| 国产精品国产三级国产av玫瑰| 内地一区二区视频在线| 男女下面进入的视频免费午夜| 两人在一起打扑克的视频| 日本免费一区二区三区高清不卡| 好男人在线观看高清免费视频| 三级男女做爰猛烈吃奶摸视频| 最近中文字幕高清免费大全6 | 欧美黑人欧美精品刺激| 观看免费一级毛片| 99热只有精品国产| 久久久久久久亚洲中文字幕| 99久久精品热视频| 日韩欧美在线二视频| 国产精品无大码| 韩国av在线不卡| 国产爱豆传媒在线观看| 成人av在线播放网站| 国产极品精品免费视频能看的| 中亚洲国语对白在线视频| 亚洲欧美日韩东京热| 18禁在线播放成人免费| 五月伊人婷婷丁香| 欧美中文日本在线观看视频| av天堂在线播放| 亚洲国产色片| 国内少妇人妻偷人精品xxx网站| 春色校园在线视频观看| 一区二区三区激情视频| 毛片女人毛片| 精品一区二区三区人妻视频| 成年人黄色毛片网站| 在线a可以看的网站| 一个人观看的视频www高清免费观看| 亚洲精华国产精华液的使用体验 | 在线观看免费视频日本深夜| 高清日韩中文字幕在线| 麻豆国产av国片精品| 亚洲,欧美,日韩| 精品乱码久久久久久99久播| 春色校园在线视频观看| 老司机深夜福利视频在线观看| 国产一区二区三区av在线 | 色综合色国产| 国产成年人精品一区二区| 亚洲欧美激情综合另类| 波野结衣二区三区在线| 久久久久久久午夜电影| 国产精品久久电影中文字幕| 亚洲成人久久性| 日日干狠狠操夜夜爽| 啦啦啦观看免费观看视频高清| 99国产极品粉嫩在线观看| 久久天躁狠狠躁夜夜2o2o| 国产在线男女| 久久99热这里只有精品18| 欧美色欧美亚洲另类二区| 在线观看免费视频日本深夜| 亚洲性夜色夜夜综合| 日本三级黄在线观看| 国产免费一级a男人的天堂| a级毛片a级免费在线| 99热只有精品国产| 久久精品国产99精品国产亚洲性色| 身体一侧抽搐| 校园人妻丝袜中文字幕| 国产又黄又爽又无遮挡在线| 搡老岳熟女国产| 麻豆国产av国片精品| 日日摸夜夜添夜夜添小说| 色吧在线观看| 三级毛片av免费| 欧美最新免费一区二区三区| 男女边吃奶边做爰视频| 国产久久久一区二区三区| 午夜福利高清视频| 国产精品久久久久久av不卡| netflix在线观看网站| 热99在线观看视频| 一级黄片播放器| 国产视频一区二区在线看| 变态另类成人亚洲欧美熟女| 国产一区二区激情短视频| 日本在线视频免费播放| 亚洲 国产 在线| 18禁裸乳无遮挡免费网站照片| 中文字幕高清在线视频| 日日撸夜夜添| 一级av片app| 国内精品美女久久久久久| 日本免费a在线| 欧美潮喷喷水| 男人舔女人下体高潮全视频| 在线免费十八禁| 婷婷亚洲欧美| 国产男人的电影天堂91| 欧美丝袜亚洲另类 | 我的女老师完整版在线观看| 97超级碰碰碰精品色视频在线观看| 一区二区三区激情视频| 欧美性感艳星| 啦啦啦韩国在线观看视频| 丰满的人妻完整版| 午夜爱爱视频在线播放| 一级黄色大片毛片| 少妇裸体淫交视频免费看高清| 国产色婷婷99| 3wmmmm亚洲av在线观看| 欧美xxxx性猛交bbbb| 又紧又爽又黄一区二区| 日本爱情动作片www.在线观看 | 赤兔流量卡办理| 国产精品一及| 婷婷丁香在线五月| 22中文网久久字幕| 嫁个100分男人电影在线观看| 国产伦在线观看视频一区| 精品久久国产蜜桃| 国产av一区在线观看免费| 久久精品久久久久久噜噜老黄 | 亚洲欧美激情综合另类| 国产精品久久久久久久电影| 亚洲七黄色美女视频| 精品久久久久久久人妻蜜臀av| 亚洲av免费在线观看| 成人亚洲精品av一区二区| av中文乱码字幕在线| 亚洲四区av| 国产视频一区二区在线看| 亚洲av熟女| 久久精品国产亚洲av涩爱 | 色综合色国产| 我要搜黄色片| 99久久精品一区二区三区| 日本精品一区二区三区蜜桃| 久久精品国产鲁丝片午夜精品 | 国产单亲对白刺激| 成人国产一区最新在线观看| 免费人成视频x8x8入口观看| 少妇的逼水好多| 亚洲自偷自拍三级| 国产在视频线在精品| 无遮挡黄片免费观看| 乱系列少妇在线播放| 婷婷精品国产亚洲av| 岛国在线免费视频观看| 精品国内亚洲2022精品成人| 国产精品电影一区二区三区| 精华霜和精华液先用哪个| 日韩 亚洲 欧美在线| 永久网站在线| 性插视频无遮挡在线免费观看| 亚洲成av人片在线播放无| 十八禁网站免费在线| 在线播放国产精品三级| eeuss影院久久| 国产单亲对白刺激| 嫩草影视91久久| 好男人在线观看高清免费视频| 精品福利观看| 校园人妻丝袜中文字幕| 日韩欧美在线乱码| 少妇的逼好多水| 老女人水多毛片| 久久婷婷人人爽人人干人人爱| 麻豆av噜噜一区二区三区| 亚洲专区中文字幕在线| 亚洲欧美清纯卡通| 婷婷精品国产亚洲av| 精品免费久久久久久久清纯| 久久久久久久久中文| 欧美黑人欧美精品刺激| 国产精品永久免费网站| 精品人妻视频免费看| 午夜激情欧美在线| 999久久久精品免费观看国产| 亚洲七黄色美女视频| 国产午夜精品论理片| 亚洲精品国产成人久久av| 免费人成在线观看视频色| 亚洲av免费在线观看| 一本精品99久久精品77| 午夜福利欧美成人| 色视频www国产| 婷婷精品国产亚洲av在线| 很黄的视频免费| av中文乱码字幕在线| 精品无人区乱码1区二区| 国产高清激情床上av| 欧美3d第一页| 人人妻人人澡欧美一区二区| 亚洲人与动物交配视频| 欧美一区二区精品小视频在线| 久久午夜亚洲精品久久| 亚洲专区国产一区二区| 日日干狠狠操夜夜爽| 黄色配什么色好看| 国产黄色小视频在线观看| 欧美绝顶高潮抽搐喷水| 亚洲狠狠婷婷综合久久图片| 国产精品人妻久久久久久| 十八禁网站免费在线| 免费观看精品视频网站| 亚洲av免费在线观看| 国产成人a区在线观看| 国产91精品成人一区二区三区| 精品人妻一区二区三区麻豆 | 少妇的逼水好多| 欧美日韩中文字幕国产精品一区二区三区| 九色国产91popny在线| 美女黄网站色视频| 国产又黄又爽又无遮挡在线| 美女高潮喷水抽搐中文字幕| 久久国产精品人妻蜜桃| 亚洲自偷自拍三级| 国产高清三级在线| 国产精品国产三级国产av玫瑰| 久久人妻av系列| www.色视频.com| 久久精品国产鲁丝片午夜精品 | 欧美高清性xxxxhd video| 一a级毛片在线观看| 一卡2卡三卡四卡精品乱码亚洲| 精品人妻视频免费看| 亚洲精品456在线播放app | 国产日本99.免费观看| 国产乱人视频| 少妇人妻精品综合一区二区 | 亚洲成人精品中文字幕电影| 午夜激情福利司机影院| 午夜免费成人在线视频| 99热这里只有是精品在线观看| 午夜影院日韩av| 美女高潮的动态| 久久国产乱子免费精品| 国产精品伦人一区二区| 3wmmmm亚洲av在线观看| 极品教师在线视频| 热99在线观看视频| 日本在线视频免费播放| 97超视频在线观看视频| 联通29元200g的流量卡| 男女边吃奶边做爰视频| 91在线观看av| 免费看美女性在线毛片视频| 黄色一级大片看看| 99久久无色码亚洲精品果冻| 在线观看66精品国产| 午夜福利在线观看吧| 全区人妻精品视频| 国产精华一区二区三区| 两个人的视频大全免费| 成人性生交大片免费视频hd| 国内精品久久久久精免费| 欧美又色又爽又黄视频| 亚洲一区高清亚洲精品| 悠悠久久av| 欧美成人一区二区免费高清观看| 国产精品乱码一区二三区的特点| 免费人成视频x8x8入口观看| 亚洲av五月六月丁香网| 亚洲精品影视一区二区三区av| 精品福利观看| 18禁裸乳无遮挡免费网站照片| 午夜福利在线观看免费完整高清在 | 中文字幕av在线有码专区| 国产一区二区亚洲精品在线观看| 99久久精品一区二区三区| 丰满乱子伦码专区| 中文资源天堂在线| 两个人视频免费观看高清| 亚洲av二区三区四区| 成人无遮挡网站| 国产白丝娇喘喷水9色精品| 久久天躁狠狠躁夜夜2o2o| 不卡一级毛片| 能在线免费观看的黄片| 国产午夜精品久久久久久一区二区三区 | 亚洲自拍偷在线| 免费看日本二区| 黄色女人牲交| 淫妇啪啪啪对白视频| 女生性感内裤真人,穿戴方法视频| 亚洲国产日韩欧美精品在线观看| 国产精品女同一区二区软件 | 岛国在线免费视频观看| 在线观看免费视频日本深夜| 久久久色成人| 久久久午夜欧美精品| 免费看av在线观看网站| 婷婷六月久久综合丁香| 亚洲三级黄色毛片| 国国产精品蜜臀av免费| 成人永久免费在线观看视频| 亚洲国产精品成人综合色| 国产精品久久久久久久电影| 好男人在线观看高清免费视频| av天堂在线播放| 老司机午夜福利在线观看视频| 国产精品久久久久久精品电影| 黄色欧美视频在线观看| 中文字幕熟女人妻在线| 老熟妇乱子伦视频在线观看| 在线观看午夜福利视频| 亚洲专区中文字幕在线| 亚洲一区高清亚洲精品| 国产精品av视频在线免费观看| 欧美一级a爱片免费观看看| 久久久久久大精品| 看免费成人av毛片| 日韩高清综合在线| 国产女主播在线喷水免费视频网站 |