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

    對(duì)流允許尺度區(qū)域氣候模擬的研究進(jìn)展

    2022-01-12 08:45:54陸云郭子悅湯劍平
    氣象科學(xué) 2021年6期
    關(guān)鍵詞:嵌套對(duì)流分辨率

    陸云 郭子悅 湯劍平

    (南京大學(xué) 大氣科學(xué)學(xué)院,南京 210023)

    引 言

    氣候模擬是理解區(qū)域氣候變化,進(jìn)行氣候預(yù)測(cè)及預(yù)估不可或缺的工具。在全球變暖背景下,區(qū)域氣候變化的預(yù)估及適應(yīng)是可持續(xù)發(fā)展所面臨的巨大挑戰(zhàn)[1]。全球氣候模式(Global Climate Model,GCM)和區(qū)域氣候模式(Regional Climate Model,RCM)在過(guò)去幾十年中經(jīng)歷了相當(dāng)大的發(fā)展[2]。雖然全球和區(qū)域氣候模式的分辨率在不斷提高,但網(wǎng)格距依然相對(duì)較粗(如區(qū)域氣候模式當(dāng)前分辨率都在20 km左右),無(wú)法顯式解析更小尺度的對(duì)流過(guò)程,因此需要采用積云對(duì)流參數(shù)化方案來(lái)表示(圖1)。盡管這種方案可以描述每個(gè)網(wǎng)格點(diǎn)上對(duì)流的平均特性,但對(duì)局地風(fēng)暴和短時(shí)強(qiáng)降水事件上的代表性較差,因此使用這種方案被認(rèn)為是區(qū)域氣候模擬不確定性的主要來(lái)源[3]。以往研究表明,對(duì)流參數(shù)化導(dǎo)致了模式出現(xiàn)一些常見(jiàn)誤差,如對(duì)流降水日循環(huán)的起始時(shí)間以及峰值不能很好地再現(xiàn)[3-5],對(duì)低降水強(qiáng)度事件的頻率有所高估,以及低估了每小時(shí)降水強(qiáng)度[3-4]等。隨著大氣數(shù)值模擬技術(shù)和高性能計(jì)算機(jī)的不斷發(fā)展,模式的水平分辨率也在不斷提高,當(dāng)分辨率提高至4 km及更高時(shí),大氣數(shù)值模式中可以避免使用對(duì)流參數(shù)化方案而直接顯式表示深對(duì)流過(guò)程,該尺度也被稱為“對(duì)流允許尺度(Convection Permitting Scale,CPS)”。這種高分辨率的區(qū)域氣候模式也能更好地獲取區(qū)域精細(xì)尺度的強(qiáng)迫信息,如在地形復(fù)雜地區(qū),對(duì)捕捉短時(shí)極端事件等方面具有明顯的增值能力[2,6]。目前對(duì)流允許尺度區(qū)域氣候模擬研究大多是通過(guò)對(duì)比對(duì)流參數(shù)化的區(qū)域氣候模擬來(lái)驗(yàn)證優(yōu)勢(shì),且大多研究區(qū)域都選擇歐洲和北美,以及非洲和亞洲的一些研究。在我國(guó),區(qū)域氣候模式已經(jīng)得到廣泛應(yīng)用,在對(duì)流允許尺度上的氣候模擬研究也在積極開(kāi)展。本文主要對(duì)以往關(guān)于對(duì)流允許區(qū)域氣候模式(Convection Permitting Regional Climate Model,CP-RCM)涉及的物理過(guò)程,外部驅(qū)動(dòng)以及模式評(píng)估結(jié)果進(jìn)行回顧,并著重關(guān)注最新CP-RCM的研究成果,以期對(duì)CP-RCM應(yīng)用于我國(guó)區(qū)域氣候模擬研究提供有益參考。

    圖1 從對(duì)流云尺度到全球尺度氣候建模Fig.1 Climate modeling from the convection-cloud scale to the global scale

    1 對(duì)流允許尺度區(qū)域氣候模擬的發(fā)展

    區(qū)域氣候模式的水平分辨率達(dá)到“對(duì)流允許尺度”時(shí)可以顯式表示深對(duì)流而不使用對(duì)流參數(shù)化方案,在這種分辨率下開(kāi)展的區(qū)域氣候模擬研究被稱為“對(duì)流允許區(qū)域氣候模擬”,研究表明這種對(duì)流允許尺度所需的水平格距一般小于5 km[7]。與“對(duì)流允許尺度”相似的術(shù)語(yǔ)有云可分辨尺度(Cloud Resolving Scale),對(duì)流分辨尺度(Convection Resolving Scale),等。云可分辨尺度模式一般指模擬云系統(tǒng)時(shí)可以分辨云尺度運(yùn)動(dòng)的二維或三維模式,典型的云可分辨尺度水平格距大約2 km,這個(gè)尺度可以解析深對(duì)流云系統(tǒng)中做大部分傳輸運(yùn)動(dòng)的大渦[8]。而達(dá)到“對(duì)流分辨尺度”,則需在云可分辨尺度基礎(chǔ)上再降低一個(gè)量級(jí),即達(dá)到大渦模擬的水平尺度才能更為真實(shí)的模擬較小的湍流夾帶以解析對(duì)流系統(tǒng)內(nèi)的降水過(guò)程[9-10]。淺積云高度達(dá)1~2 km,云的側(cè)邊緣發(fā)生的夾帶過(guò)程由大渦模擬能更好的表現(xiàn)出來(lái),CP-RCM對(duì)淺對(duì)流可能不足以解析[8,11]。因此目前大渦模擬也被廣泛用于淺對(duì)流參數(shù)化的評(píng)估和開(kāi)發(fā),但受計(jì)算資源的限制大渦模擬應(yīng)用于長(zhǎng)期區(qū)域氣候模擬還不太現(xiàn)實(shí)[11]。當(dāng)氣候模式達(dá)到CPS時(shí)可以顯示表示深對(duì)流,但不足解決湍流以及邊界層小尺度過(guò)程,4 km的氣候模式在尺度上并沒(méi)有完全解析對(duì)流過(guò)程,故稱為“對(duì)流允許”。

    自1990s開(kāi)始,CPS模式在數(shù)值天氣預(yù)報(bào)方面的應(yīng)用陸續(xù)開(kāi)展,有研究發(fā)現(xiàn)水平格距為2 km的精細(xì)網(wǎng)格足以顯式捕捉中尺度對(duì)流系統(tǒng)(Mesoscale Convective Systems,MCSs)和超級(jí)單體雷暴[8,12]。然而在區(qū)域氣候模擬方面,CPS模式的應(yīng)用稍晚,在近十幾年才有所發(fā)展,如2000s初期Herrmann, et al[13]對(duì)1986—1987年冬季地中海西北部開(kāi)展數(shù)值模擬(10 km和3 km)來(lái)研究深對(duì)流活動(dòng)中的中尺度過(guò)程,并發(fā)現(xiàn)3 km的CP-RCM更好地再現(xiàn)了中尺度結(jié)構(gòu)特征。近幾年CP-RCM的研究遍及歐洲、美洲、非洲、亞洲。利用CP-RCM進(jìn)行氣候模擬可避免與積云參數(shù)化方案相關(guān)的不確定性[14],其模擬深對(duì)流降水方面的增值能力已得到廣泛證實(shí)。除了顯式地分辨深對(duì)流外,CP-RCM還具有改善細(xì)尺度地形和變化的表面場(chǎng)表示的優(yōu)點(diǎn),這在山區(qū),沿海和城市地區(qū)等尤其有益[2-4,14]。

    對(duì)流允許尺度區(qū)域氣候模擬主要是通過(guò)GCM或全球再分析數(shù)據(jù)提供初始和側(cè)邊界條件,驅(qū)動(dòng)高分辨率的RCM進(jìn)行季節(jié)、年和年代際的長(zhǎng)期區(qū)域氣候模擬。近十幾年全球各大區(qū)域開(kāi)展的一些CP-RCM研究中所運(yùn)用的降尺度策略(表1)及其研究成果表明,CP-RCM最初在歐美地區(qū)開(kāi)展的比較多,如歐洲的一些國(guó)家利用COSMO-CLM(COSMO model in Climate Model)模式在阿爾卑斯山脈、德國(guó)、比利時(shí)等地區(qū)開(kāi)展了很多CP-RCM的研究。眾多試驗(yàn)表明CP-RCM明顯改善了對(duì)日降水周期特征的模擬,且與對(duì)流降水有關(guān)的環(huán)流場(chǎng)具有較好的一致性。高空間分辨率在一定程度上改善了對(duì)復(fù)雜地形降水場(chǎng)以及溫度空間變異性的表征[3,15-17]。北美大陸地區(qū)開(kāi)展基于WRF(Weather Research and Forecasting Model)的CP-RCM模擬,發(fā)現(xiàn)CP-RCM不僅能模擬出對(duì)流降水量大小,日周期以及相關(guān)的大氣環(huán)流,而且在未來(lái)全球變暖情景下,對(duì)北美地區(qū)水循環(huán)的變化,包括美國(guó)中部MCSs路徑相關(guān)的暖季降水的模擬表現(xiàn)良好[4,18-20]。此外,在非洲地區(qū)利用CP-RCM模擬結(jié)果也很好解釋出西非地區(qū)白天對(duì)流降水發(fā)生的有利天氣環(huán)流背景[21-22],以及非洲薩赫勒地區(qū)土壤水分與對(duì)流觸發(fā)的空間關(guān)系[23]。近年來(lái),東亞地區(qū)也開(kāi)展了基于對(duì)流允許尺度區(qū)域氣候模擬的研究工作,如對(duì)日本西南海域年和月平均降水量以及日和小時(shí)降水量頻率分布的評(píng)估[14]。我國(guó)CP-RCM的模擬研究工作也在陸續(xù)展開(kāi),如使用WRF模式動(dòng)力降尺度得到的高分辨率數(shù)據(jù)集很好地再現(xiàn)了我國(guó)長(zhǎng)三角地區(qū)(Yangtze River Delta,YRD)地區(qū)以及我國(guó)西北內(nèi)陸地區(qū)降水特點(diǎn)[24-25],同時(shí),有一些研究利用CP-RCM探討了城市熱島效應(yīng)對(duì)我國(guó)YRD地區(qū)氣候的影響[26]。CP-RCM模擬降水的優(yōu)越性在季風(fēng)區(qū)和干旱區(qū)得到證實(shí),同時(shí)對(duì)具有多樣下墊面的復(fù)雜地理環(huán)境的降水場(chǎng)特征也能很好捕獲[14,21-27]。

    表1 對(duì)流允許尺度區(qū)域氣候模擬概述Table 1 Overview of Convection Permitting Regional Climate simulation

    本文將主要從CP-RCM需要考慮的重要物理過(guò)程,外部驅(qū)動(dòng)以及其在模擬降水的增值能力方面的研究進(jìn)展和目前面臨的問(wèn)題展開(kāi)討論。

    2 CP-RCM模擬相關(guān)的物理過(guò)程

    隨著分辨率的提高,對(duì)流允許尺度區(qū)域氣候模擬中直接顯式表示深對(duì)流過(guò)程而不采用對(duì)流參數(shù)化方案,避免了積云對(duì)流參數(shù)化方案帶來(lái)的不確定性。而模式模擬中其他次網(wǎng)格尺度物理過(guò)程參數(shù)化方案的選擇有可能對(duì)CP-RCM氣候模擬效果產(chǎn)生明顯的影響,總結(jié)關(guān)于CP-RCM氣候模擬中一些至關(guān)重要的物理過(guò)程,主要包括:云微物理、土壤—大氣相互作用、邊界層過(guò)程和湍流[19,28],這些物理過(guò)程成為CP-RCM模擬不確定性的主要來(lái)源。

    2.1 與云微物理過(guò)程相關(guān)的不確定性

    Krueger[8]回顧了使用云分辨模式來(lái)模擬深對(duì)流和淺對(duì)流云過(guò)程的研究,發(fā)現(xiàn)由于模式顯式表示深對(duì)流過(guò)程,對(duì)流的開(kāi)始以及單個(gè)云狀態(tài)上的演化非常依賴于云微物理過(guò)程的參數(shù)化。云微物理過(guò)程是模擬氣溶膠—云相互作用不確定性的主要來(lái)源之一,特別是對(duì)混合相和冰相云的研究[28]。一些研究認(rèn)為在云微物理參數(shù)化方案中引入額外的水凝物類型,例如霰和冰雹,或使用雙參微物理方案模擬云粒子的數(shù)濃度,在改善氣溶膠與云之間相互作用的模擬上有很大潛力,有利于提高對(duì)冬季地形降水和高云覆蓋的模擬能力[29-31]。但也有研究認(rèn)為單參方案能夠比雙參方案更好地捕捉降水特征和反射率分布等方面,因?yàn)樵陔p參方案中雨滴的大小分布可能變得對(duì)許多未知過(guò)程非常敏感[32]。FENG, et al[19]和GUO, et al[33]分別利用WRF模式對(duì)美國(guó)暖季和中國(guó)夏季進(jìn)行對(duì)流允許尺度(4 km網(wǎng)格間距)區(qū)域氣候模擬,分析了不同微物理方案對(duì)降水模擬的影響。其結(jié)果表明,應(yīng)用不同云微物理方案的CP-RCM均能很好地模擬出降水日變化及其傳播特征,但是不同方案對(duì)于水凝物不同處理導(dǎo)致降水量的模擬存在差異。這些研究說(shuō)明CP-RCM對(duì)與云微物理相關(guān)過(guò)程的敏感性很大,云微物理參數(shù)化效果可能在不同案例間也存在差別。此外云微物理過(guò)程也會(huì)影響CP-RCM中的輻射過(guò)程,云輻射反饋在很大程度上依賴于云量和云輻射特性的準(zhǔn)確表示,而當(dāng)前云—?dú)馊苣z—輻射相互作用的參數(shù)化仍然知之甚少,并且在CP-RCM中存在較多的不確定性[28]。CP-RCM對(duì)云微物理過(guò)程的表示具有高度敏感性,而氣溶膠對(duì)氣候系統(tǒng)會(huì)產(chǎn)生廣泛影響,其對(duì)云層的相互作用的參數(shù)方案需要更詳細(xì)的研究,因此在區(qū)域氣候模式中改進(jìn)氣溶膠與云相互作用的表示非常重要。

    2.2 對(duì)邊界層和湍流的處理方法

    除了云微物理過(guò)程外,邊界層和湍流過(guò)程也是影響CP-RCM的重要物理過(guò)程。由于CPS遠(yuǎn)不足以解析湍流尺度,地面和大氣之間在次網(wǎng)格尺度上熱量、水分和動(dòng)量的垂直傳輸,需要使用行星邊界層(Planetary Boundary Layer,PBL)方案來(lái)參數(shù)化。區(qū)域氣候模式對(duì)不同PBL參數(shù)化方案敏感,PBL方案本身結(jié)構(gòu)性偏差以及對(duì)參數(shù)的敏感性為CP-RCM帶來(lái)不確定性[34]。Cohen, et al[35]總結(jié)了WRF模式中 PBL方案的主要特點(diǎn),并通過(guò)對(duì)美國(guó)冷季龍卷風(fēng)事件采用的PBL方案進(jìn)行測(cè)試,發(fā)現(xiàn)使用PBL的非局地參數(shù)化的WRF模擬產(chǎn)生了比只使用局地方案更精確的結(jié)果。Kouadio, et al[36]對(duì)西非的研究結(jié)果也表明非局地PBL在模式降水模擬中有增值能力。Larson, et al[37]在云分辨模型中加入一種云和湍流的參數(shù)化(Clubb參數(shù)化),發(fā)現(xiàn)對(duì)毛毛雨混合比的預(yù)測(cè)取得了很好的模擬效果。這些對(duì)參數(shù)化評(píng)估的研究有助于我們?cè)谶M(jìn)行CP-RCM模擬時(shí)選取較為合適的PBL方案。為擺脫常規(guī)的邊界層參數(shù)化對(duì)氣候模式的發(fā)展限制,一些超參數(shù)化模式也得到了發(fā)展,通過(guò)在每個(gè)GCM網(wǎng)格列中嵌入小的(250 m×20 m)云解析模式(Cloud Resolving Model,CRM),可以顯式捕獲邊界層湍流和夾帶過(guò)程[38]。這種在湍渦網(wǎng)格分辨率下對(duì)邊界層和深對(duì)流云進(jìn)行的全局模擬取得很好的效果,對(duì)未來(lái)區(qū)域氣候模擬邊界層和湍流的處理提供了思路。

    2.3 由土壤—大氣相互作用產(chǎn)生的影響

    在CP-RCM氣候模擬中陸面過(guò)程也是一個(gè)關(guān)鍵部分。眾所周知,土壤—大氣相互作用對(duì)許多大氣過(guò)程很重要,比如土壤濕度動(dòng)態(tài)變化會(huì)影響著感熱和潛熱的分配,甚至影響邊界層狀態(tài),從而對(duì)中尺度環(huán)流產(chǎn)生影響。CP-RCM通過(guò)精細(xì)分辨率來(lái)模擬地表與大氣的能量交換和相互作用,會(huì)改進(jìn)區(qū)域氣候模擬能力[39-40]。Taylor, et al[23]研究認(rèn)為,4 km的CP-RCM模擬可以再現(xiàn)土壤濕度和對(duì)流觸發(fā)之間的空間關(guān)系,表明陸面過(guò)程對(duì)中尺度過(guò)程發(fā)生發(fā)展的重要性。許多相關(guān)研究發(fā)現(xiàn)CP-RCM可以更真實(shí)地模擬這種反饋,但模式對(duì)土壤水分變化非常敏感,因?yàn)橥寥浪值淖兓ㄟ^(guò)影響地氣能量交換影響大氣邊界層,會(huì)導(dǎo)致模擬結(jié)果產(chǎn)生較大的不確定性[41-42],因此也成為CP-RCM模擬不確定性的主要來(lái)源之一[41]。CP-RCM對(duì)土壤水分—降水反饋模擬的良好表現(xiàn)可能來(lái)源于其精細(xì)網(wǎng)格尺度對(duì)局部效應(yīng)的捕獲,而土壤或地表中的許多諸如蒸發(fā)過(guò)程之類的高度非線性過(guò)程仍然沒(méi)有得到解決,陸—?dú)怦詈蠈?duì)土壤水分動(dòng)態(tài)的敏感性表明,對(duì)于CP-RCM來(lái)說(shuō),精確識(shí)別大氣所涉及的復(fù)雜反饋過(guò)程仍然是個(gè)挑戰(zhàn)。

    3 城市氣候與未來(lái)區(qū)域氣候模擬

    城市化引起的土地覆蓋變化和溫室氣體排放的增加對(duì)區(qū)域氣候產(chǎn)生很大影響。城市化通過(guò)改變地表特征影響陸地表面與大氣之間的能量和水交換,同時(shí)城市建筑材料通過(guò)改變反照率來(lái)影響輻射平衡,進(jìn)而影響區(qū)域氣候[43-44]。污染物排放通過(guò)氣溶膠—輻射相互作用來(lái)影響地表和大氣輻射通量和輻射加熱過(guò)程[26,45]。眾多研究表明更高的分辨率更真實(shí)地代表了在非均勻地表上與大氣水分和能量交換以及對(duì)變化的陸—?dú)怦詈袭a(chǎn)生的大氣過(guò)程反饋的響應(yīng),能夠再現(xiàn)土地利用/覆蓋引起的局地氣候的變化[20,39]。有研究表明當(dāng)RCM空間分辨率達(dá)到3 km,大城市區(qū)域可以得到足夠詳細(xì)的分辨,引入多個(gè)城市土地利用類別會(huì)有效地改善對(duì)城市熱島的模擬[15]。許多研究也發(fā)現(xiàn)CP-RCM的精細(xì)網(wǎng)格尺度對(duì)于模擬城市氣候變化具有明顯的優(yōu)勢(shì)[26,46],為改進(jìn)城市化對(duì)未來(lái)氣候變化影響的評(píng)估開(kāi)辟了道路。我國(guó)長(zhǎng)三角地區(qū)處于典型的季風(fēng)氣候區(qū),且城市化水平較高,在該地區(qū)開(kāi)展了很多CP-RCM模擬。例如用WRF模式在對(duì)流允許尺度模擬長(zhǎng)三角區(qū)域極端氣候發(fā)現(xiàn),CP-RCM能夠模擬出夏季降水的空間分布特征以及強(qiáng)降水事件的時(shí)間演變特征[26]。利用CP-RCM開(kāi)展的城市化區(qū)域氣候效應(yīng)的研究,也再現(xiàn)了長(zhǎng)三角地區(qū)氣溫和降水的變化特征以及由于城市化引發(fā)的城市熱島效應(yīng)以及影響機(jī)制[26,46]。諸多研究表明在對(duì)流允許尺度上,土地利用和覆蓋變化產(chǎn)生的局部效應(yīng)比以往所認(rèn)為的更為重要,CP-RCM精細(xì)的網(wǎng)格分布通過(guò)捕獲地形和土地覆蓋相關(guān)的區(qū)域強(qiáng)迫和過(guò)程,使模式能夠更好解決相關(guān)的物理過(guò)程及其與大氣的相互作用,從而增強(qiáng)了區(qū)域氣候模擬能力[20,27]。

    隨著氣候變暖,大氣水汽含量也顯著上升,未來(lái)的降水事件將變得更加強(qiáng)烈[47],對(duì)未來(lái)極端事件的預(yù)估很有必要。CP-RCM模擬所采用的分辨率可以把中尺度系統(tǒng)更好地表現(xiàn)出來(lái),其對(duì)當(dāng)下極端天氣氣候事件更真實(shí)的模擬能力為未來(lái)氣候預(yù)測(cè)增加了信心,因此CP-RCM也廣泛用于區(qū)域氣候變化的預(yù)估[18,48-49]?;贑P-RCM的模擬預(yù)測(cè)發(fā)現(xiàn)在RCP8.5排放情景下日本極端溫度表現(xiàn)出廣泛增加,而7月極端降水的增加則出現(xiàn)區(qū)域性[48]。Broucke, et al[50]基于CP-RCM開(kāi)展的未來(lái)極端降水預(yù)測(cè)研究也發(fā)現(xiàn)夏季極端降水的變化信號(hào)可能與地形有關(guān),與RCM相比,CP-RCM模擬出更大的夏季極端降水增幅,而兩者對(duì)平坦地區(qū)小時(shí)降水強(qiáng)度的模擬出現(xiàn)相反的信號(hào)。在對(duì)北美地區(qū)當(dāng)前和未來(lái)的CP-RCM模擬也發(fā)現(xiàn),對(duì)流允許尺度區(qū)域氣候模擬并不一定能提供更多有益的未來(lái)極端事件氣候變化信息[18]。在未來(lái)極端降水的模擬方面,CP-RCM模擬并沒(méi)有出現(xiàn)普適性的結(jié)果,需要開(kāi)展多個(gè)CP-RCM的集合預(yù)測(cè)來(lái)減少氣候變化背景下的極端事件預(yù)估的不確定性[51]。

    4 CP-RCM外部驅(qū)動(dòng)的影響

    在CP-RCM模擬中,大尺度信息是通過(guò)初始和側(cè)邊界條件來(lái)驅(qū)動(dòng)模式的。影響CP-RCM模擬的外部驅(qū)動(dòng)因子包括嵌套方案,初始和側(cè)邊界條件以及模擬區(qū)域大小和模式運(yùn)行的spin-up時(shí)間等。值得注意的是,這些因素的影響都是交織在一起的,有時(shí)很難對(duì)這些誤差分量進(jìn)行量化和區(qū)分。

    為了達(dá)到對(duì)流允許尺度的分辨率,CP-RCM往往通過(guò)嵌套網(wǎng)格的方式來(lái)實(shí)現(xiàn),嵌套可以是單向的也可以是雙向的。單向嵌套是指高分辨率的較小區(qū)域所需的側(cè)邊界條件由粗分辨率的較大區(qū)域提供,而計(jì)算結(jié)果不反饋到較大區(qū)域。而雙向嵌套是指粗、細(xì)分辨率模擬同時(shí)運(yùn)行,并且細(xì)網(wǎng)格將其計(jì)算值反饋到粗網(wǎng)格的區(qū)域[52]。Prein, et al[16]表明,在東阿爾卑斯山進(jìn)行的CP-RCM氣候模擬在單向和雙向嵌套方案之間沒(méi)有顯著差異。另外一些研究發(fā)現(xiàn),對(duì)于沿海的海洋模擬,雙向嵌套模擬海岸大尺度環(huán)流比單向嵌套的效果要好[53-54]。值得注意的是,在動(dòng)力降尺度過(guò)程中,對(duì)于有限的區(qū)域內(nèi)將GCM輸出縮減到區(qū)域尺度所采用的“跳躍分辨率”(粗分辨率與細(xì)分辨率網(wǎng)格距之比)并沒(méi)有一致的認(rèn)定。有研究發(fā)現(xiàn),小尺度特征的發(fā)展需要較長(zhǎng)的spin-up時(shí)間,需要考慮降尺度過(guò)程中跳躍分辨率對(duì)小尺度的發(fā)展的影響[55]。而通過(guò)使用多重嵌套降低了空間spin-up對(duì)小尺度發(fā)展的制約,給小尺度的發(fā)展產(chǎn)生了積極的影響[56]。跳躍分辨率的選擇又與嵌套次數(shù)有關(guān),利用不同嵌套步驟進(jìn)行比利時(shí)臭氧濃度的動(dòng)力降尺度(25 km到3 km)研究中發(fā)現(xiàn),增加中間嵌套步驟與從RCM到CP-RCM的直接嵌套模擬結(jié)果相似,而后者更加節(jié)省計(jì)算資源[57]。Brisson, et al[30]還認(rèn)為在ERA-Interim降尺度至CPS時(shí)增加7 km網(wǎng)格嵌套,模擬沒(méi)有產(chǎn)生明顯的影響,因此可以省略這一嵌套步驟,但是他們認(rèn)為25 km中間嵌套對(duì)于模擬降水正確表示至關(guān)重要。目前尚不清楚不同嵌套方案是如何影響CP-RCM模擬的,對(duì)于嵌套方法的選擇,既希望帶來(lái)好的模擬效果,又能更好的節(jié)約計(jì)算資源,這可能需要對(duì)特定事件進(jìn)行廣泛評(píng)估。

    由于對(duì)流過(guò)程是高度非線性的,并且一般持續(xù)時(shí)間較短,因此它們的預(yù)測(cè)受到不確定因素的強(qiáng)烈影響,例如模式的初始條件、參數(shù)化設(shè)置以及側(cè)邊界的驅(qū)動(dòng)場(chǎng)等[16,25,55-61]。一些研究評(píng)估了初始和側(cè)邊界條件不確定性對(duì)CP-RCM模擬的影響,發(fā)現(xiàn)對(duì)流尺度初始條件的不確定性強(qiáng)烈依賴于天氣尺度強(qiáng)迫[62-63]。在12 h以內(nèi)的天氣預(yù)報(bào)中,初始條件產(chǎn)生較大影響,而在長(zhǎng)期區(qū)域氣候模擬中側(cè)邊界條件的不確定性則產(chǎn)生越來(lái)越大的影響[62]。Tennant, et al[63]使用不同的方法來(lái)初始化Met Office有限區(qū)域?qū)α髟试S集合預(yù)測(cè)系統(tǒng),相比于目前流行的動(dòng)力降尺度方法和只更新側(cè)邊界驅(qū)動(dòng)的運(yùn)行方法有明顯優(yōu)勢(shì)。近期全球可變網(wǎng)格模式的發(fā)展有利于解決在有限區(qū)域模式中出現(xiàn)的側(cè)邊界條件問(wèn)題。可變網(wǎng)格模式能夠在區(qū)域上平滑地改變分辨率,有助于協(xié)調(diào)區(qū)域模擬和驅(qū)動(dòng)場(chǎng)之間在側(cè)邊界附近的差異,減少模擬誤差并且大大降低成本[59-64]。

    由于受到粗分辨率的初始和側(cè)邊界條件的驅(qū)動(dòng)強(qiáng)迫,CP-RCM需要一些spin-up時(shí)間以產(chǎn)生細(xì)尺度的特征,同時(shí)區(qū)域大小設(shè)置對(duì)CP-RCM也有較大影響,因此在運(yùn)行CP-RCM時(shí)需要考慮模擬域的大小以及模擬開(kāi)始的時(shí)間的選擇[30]。Brisson, et al[30]研究了CP-RCM模擬區(qū)域大小和嵌套方法對(duì)降水模擬的影響,發(fā)現(xiàn)對(duì)流允許尺度的嵌套區(qū)域并沒(méi)有比其嵌套的粗網(wǎng)格區(qū)域有更好的模擬效果,這主要是因?yàn)楦叻直媛实哪M需要有足夠的spin-up時(shí)間。和RCM一樣,CP-RCM也需要仔細(xì)選擇模擬區(qū)域大小以實(shí)現(xiàn)適當(dāng)?shù)慕党叨?如果域尺寸很小則CP-RCM可能沒(méi)有足夠的時(shí)間和空間來(lái)生成小尺度細(xì)節(jié),而如果模擬域很大則CP-RCM和側(cè)邊界條件之間可能存在明顯的不匹配,特別是在流出邊界處[52]。諸多研究表明區(qū)域氣候模擬對(duì)選擇用于計(jì)算的域的大小和位置敏感。因此許多人嘗試使用譜逼近來(lái)消除模擬結(jié)果對(duì)區(qū)域選擇的依賴性。對(duì)于區(qū)域模式大小達(dá)到幾千公里的模擬試驗(yàn),長(zhǎng)波的譜逼近有助于精確模擬大尺度環(huán)流以及消除邊界的虛假影響,同時(shí)對(duì)極端氣候和地面氣候時(shí)空變化的模擬有很大改善[65-66]。但是在模擬區(qū)域相對(duì)較小的情況下,譜逼近的應(yīng)用對(duì)模擬效果影響不大[67]。

    5 CP-RCM的增值能力

    CP-RCM的優(yōu)點(diǎn)是能夠顯式分辨深對(duì)流過(guò)程而避免使用對(duì)流參數(shù)化方案,同時(shí)高分辨率也改善了對(duì)地形和其他陸面強(qiáng)迫的表示,這些都是CP-RCM增值能力的驅(qū)動(dòng)因子。與GCM輸出相比,CP-RCM并不會(huì)在所有地區(qū),季節(jié)和物理量模擬中都有明顯的增值能力[2]。春季和秋季由于局地混合和大尺度強(qiáng)迫條件,CP-RCM的增值能力并不明顯,而在夏季對(duì)流降水占主導(dǎo)地位[3,50],研究表明CP-RCM對(duì)夏季降水的空間分布形式、降水日變化和極端降水的強(qiáng)度等方面的模擬同大尺度氣候模式相比改進(jìn)明顯[2-4,16,19]。

    5.1 夏季降水的日變化

    在CP-RCM的增值方面,對(duì)夏季降水日循環(huán)起始以及峰值的捕獲能力在許多研究中得到證實(shí)。Fosser, et al[3]研究了不同分辨率的區(qū)域氣候模式對(duì)德國(guó)降水日變化的模擬能力,發(fā)現(xiàn)關(guān)閉了對(duì)流參數(shù)化方案的CP-RCM模擬與實(shí)際觀測(cè)更為接近,且與對(duì)流降水有關(guān)的大尺度環(huán)流之間存在更好的一致性。同樣地,基于WRF模式對(duì)美國(guó)大陸的CP-RCM模擬試驗(yàn)?zāi)芨诱鎸?shí)地模擬出與夏季降水相關(guān)的低空急流和大氣環(huán)流,并且更精確地再現(xiàn)了對(duì)流降水出現(xiàn)的時(shí)間以及從洛基山脈到大平原的雨帶日循環(huán)的東傳信號(hào)[4]。MCSs是造成對(duì)流強(qiáng)降水重要的天氣系統(tǒng),在美國(guó)中部[19]以及美國(guó)西南部盆地[20]的CP-RCM模擬再次說(shuō)明CP-RCM對(duì)MCSs模擬的可信度。在非洲地區(qū)的模擬試驗(yàn)也表明CP-RCM能較好地再現(xiàn)西非地區(qū)降雨的日變化特征,而且能夠模擬出主導(dǎo)降雨的MCSs的生命周期[21-22]。

    在中國(guó)開(kāi)展的CP-RCM模擬研究也表明其能夠改善模式對(duì)降水日變化特征的模擬能力。ZHU, et al[24]運(yùn)用WRF模式在4 km對(duì)流允許分辨率下準(zhǔn)確地再現(xiàn)了強(qiáng)降水帶的空間范圍和方向以及每小時(shí)降雨量的晝夜時(shí)間序列。WRF能預(yù)測(cè)出2個(gè)降水峰值,但是對(duì)峰值略微過(guò)高預(yù)測(cè)。FU, et al[68]在中國(guó)東部大別山地區(qū)的模擬研究發(fā)現(xiàn),CP-RCM對(duì)降水的雙峰特征進(jìn)行了較好的模擬,與觀測(cè)結(jié)果在量級(jí)上有很好的吻合,但對(duì)降水峰值出現(xiàn)時(shí)間的捕獲稍有延后。在喜馬拉雅中部地區(qū)CP-RCM模擬真實(shí)地再現(xiàn)了降水日周期及其峰值出現(xiàn)的時(shí)間,相比于粗分辨率模式,對(duì)地形氣候模擬效果更好[27]。眾多研究表明CP-RCM更好地反映地形強(qiáng)迫以及解析大氣中尺度過(guò)程,在對(duì)流降水的日變化上良好的捕獲能力得到廣泛證實(shí)。

    5.2 平均降水和極端降水

    一些研究表明CP-RCM與較粗網(wǎng)格的RCM模擬相比在日平均降水量上并沒(méi)有得到普遍改善[3-5],如SUN, et al[4]對(duì)美國(guó)大平原夏季平均降水分布模擬的評(píng)估表明4 km和25 km的模擬在某些地區(qū)存在相似的降水偏差,這可能是模式系統(tǒng)性偏差所致,而CP-RCM在對(duì)降水分布的表現(xiàn)總體上是優(yōu)于RCM的,特別是在局部變率的細(xì)節(jié)表示方面。在次季節(jié)/季節(jié)/年降水量和溫度的空間和時(shí)間分布上,CP-RCM模擬性能優(yōu)于較粗網(wǎng)格的模式[14,18,24-25]。Hart, et al[69]對(duì)南非亞熱帶地區(qū)長(zhǎng)達(dá)10 a的模擬表明,CP-RCM能準(zhǔn)確捕捉熱帶—溫帶區(qū)云帶的年循環(huán),并且對(duì)當(dāng)?shù)叵募窘涤昴M的濕偏差顯著減小。雖然深對(duì)流參數(shù)化在CP-RCM中已經(jīng)關(guān)閉,淺對(duì)流方案的影響依然很重要,淺對(duì)流方案起到穩(wěn)定對(duì)流層的作用,一定程度上阻礙了對(duì)流的發(fā)展和降水的產(chǎn)生,從而有利于減少濕偏差[49]。通過(guò)對(duì)關(guān)閉和不關(guān)閉對(duì)流參數(shù)化的RCM模擬對(duì)比研究,發(fā)現(xiàn)在采用對(duì)流參數(shù)方案的粗分辨率RCM中往往會(huì)高估低強(qiáng)度降水頻率,這可能是RCM在對(duì)流層低層的模擬偏濕冷,有利于低云形成從而高估低強(qiáng)度降水事件,而相反在CP-RCM中,由于直接顯式表示了深對(duì)流過(guò)程這種情況得到改善[3]。

    CP-RCM對(duì)極端降水和強(qiáng)降水的模擬也有一定的增值能力。Murata, et al[14]通過(guò)2 km網(wǎng)格間距的CP-RCM很好地再現(xiàn)了降水強(qiáng)度低于20 mm·h-1的降水頻率,這與Fosser, et al[3]在德國(guó)西南部以2.8 km分辨率所做試驗(yàn)的結(jié)果一致。在極端事件方面,有研究表明高分辨率(1.3 km)模擬的德國(guó)當(dāng)前極端降水與觀測(cè)較為接近,其對(duì)極端降水的模擬比7 km有更好的表現(xiàn)[49]。Gensini, et al[70]的研究也說(shuō)明對(duì)美國(guó)大陸部分地區(qū)的極端天氣(龍卷風(fēng)、陣風(fēng)和冰雹)的變化以4 km分辨率的CP-RCM進(jìn)行模擬效果良好。在全球變暖背景下,極端降水變得更短時(shí)更強(qiáng)烈,可能造成嚴(yán)重的自然災(zāi)害,CP-RCM對(duì)小時(shí)極端降水的表示方面能力對(duì)未來(lái)極端氣候的預(yù)估增強(qiáng)了信心[71]。

    盡管CP-RCM的水平分辨率很高,可以合理地模擬出MCSs造成的極端強(qiáng)降水的相關(guān)特征,但如前所述,MCSs屬性在不同微物理參數(shù)化方案中差異很大,云微物理方案對(duì)MCSs的不同處理不僅改變了MCSs對(duì)流尺度動(dòng)力學(xué),而且對(duì)其生命周期有很大影響。CP-RCM中MCSs的相關(guān)特性對(duì)不同微物理參數(shù)化方案很敏感,CP-RCM對(duì)MCSs模擬的增值能力強(qiáng)烈依賴于云微物理方案的選擇[19]。

    5.3 地面氣溫

    相對(duì)于夏季對(duì)流降水的模擬,CP-RCM對(duì)模擬地面氣溫時(shí)空分布的增值能力并不顯著。Leutwyler, et al[5]對(duì)歐洲大陸進(jìn)行了為期10 a的CP-RCM模擬,發(fā)現(xiàn)模式對(duì)夏季地面氣溫度的模擬有較大的暖偏差,這可能是由于CP-RCM中顯式地表示深對(duì)流云從而增加了地表太陽(yáng)輻射。LIU, et al[18]的研究也發(fā)現(xiàn)WRF模式對(duì)美國(guó)大陸大部分時(shí)間地面氣溫模擬存在暖偏差,其中美國(guó)中部地區(qū)夏季的溫度偏差最顯著。同樣Brisson, et al[17]也認(rèn)為CP-RCM對(duì)地面氣溫模擬方面的改進(jìn)不太明顯。雖然Prein, et al[16]的研究發(fā)現(xiàn)CP-RCM對(duì)夏季2 m高度的平均溫度模擬有所改進(jìn),然而他們認(rèn)為這種改進(jìn)很難歸因于CP-RCM的優(yōu)越性,因?yàn)檫@也可以通過(guò)非常簡(jiǎn)單的高度校正來(lái)實(shí)現(xiàn)。

    6 挑戰(zhàn)和展望

    本文基于對(duì)CP-RCM相關(guān)文獻(xiàn)的系統(tǒng)性回顧,討論了CP-RCM主要涉及的物理過(guò)程,外部驅(qū)動(dòng)條件以及氣候(降水和溫度)的模擬能力等方面。綜合來(lái)看,對(duì)流允許尺度區(qū)域氣候模擬顯著提升了我們對(duì)區(qū)域氣候變化的理解,尤其是在地形復(fù)雜多變的地區(qū),CP-RCM在對(duì)流降水的再現(xiàn)能力方面改進(jìn)明顯。但CP-RCM的模擬能力還有待進(jìn)一步提高,比如陣雨或者小雨的模擬不準(zhǔn)確,并且對(duì)極端降水容易高估。通過(guò)對(duì)CP-RCM涉及物理過(guò)程和外部驅(qū)動(dòng)條件等方面研究的探討,發(fā)現(xiàn)CP-RCM在建模過(guò)程中仍然存在很多挑戰(zhàn)和不確定性。

    首先,通過(guò)對(duì)其物理過(guò)程的討論總結(jié),發(fā)現(xiàn)CP-RCM對(duì)次網(wǎng)格尺度過(guò)程的表示相當(dāng)敏感,與許多未知物理過(guò)程相關(guān)聯(lián),如邊界層,微物理過(guò)程和陸面過(guò)程等參數(shù)化方案,與物理參數(shù)化相關(guān)的參數(shù)敏感性以及系統(tǒng)性偏差成為CP-RCM主要不確定性來(lái)源,未來(lái)對(duì)其中物理成因進(jìn)行更全面的探索以及模式參數(shù)方案的改進(jìn)很有必要,特別是對(duì)云—?dú)馊苣z—輻射相互作用的理解。此外,還需重新評(píng)估它們?cè)贑P-RCM上的適用性,比如有研究指出CPS模式對(duì)次網(wǎng)格湍流方案中混合長(zhǎng)度值的具有敏感性,僅僅增加分辨率并不一定有助于對(duì)對(duì)流模擬的提升,可能需要重新考慮次網(wǎng)格湍流混合公式才能更好地描述在100~1 000 m的網(wǎng)格長(zhǎng)度范圍內(nèi)對(duì)流云的性質(zhì)[72]。同時(shí),除了對(duì)降水模擬能力的評(píng)估外,應(yīng)進(jìn)一步探索CP-RCM在其他相關(guān)領(lǐng)域,包括局地風(fēng)、水文、陸—?dú)庀嗷プ饔?、蒸發(fā)蒸騰、湍流通量以及云和輻射等過(guò)程的應(yīng)用能力[6]。

    其次,CP-RCM在降尺度過(guò)程中初始和側(cè)邊界條件,嵌套方法以及模擬域大小的選擇對(duì)其模擬氣候的能力也有很大影響。較長(zhǎng)時(shí)間尺度的區(qū)域氣候模擬主要受驅(qū)動(dòng)場(chǎng)提供的側(cè)邊界影響,初始條件影響不大,高時(shí)間分辨率側(cè)邊界條件能帶來(lái)較好的模擬效果,而區(qū)域模擬和驅(qū)動(dòng)場(chǎng)之間在側(cè)邊界附近的差異可通過(guò)變網(wǎng)格技術(shù)得到改善。相較于單向嵌套,雙向嵌套對(duì)海岸附近大尺度環(huán)流模擬較好,并且適當(dāng)增加嵌套步驟有利于對(duì)小尺度特征的描述。此外在進(jìn)行CP-RCM時(shí),實(shí)驗(yàn)區(qū)域不宜太小,譜逼近的應(yīng)用會(huì)對(duì)模擬效果帶來(lái)積極影響,同時(shí)模式還需要一些spin-up時(shí)間來(lái)達(dá)到平衡狀態(tài)。

    此外,由于沒(méi)有足夠的地面測(cè)量數(shù)據(jù)來(lái)描繪降水的高分辨率空間分布特征,而基于衛(wèi)星反演的降水?dāng)?shù)據(jù)集的空間分辨率比CP-RCM分辨率又粗糙很多,因此開(kāi)發(fā)高時(shí)空分辨率的觀測(cè)數(shù)據(jù)集來(lái)評(píng)估模擬也很有必要。

    最近,有研究利用全球可變網(wǎng)格模式開(kāi)展了熱帶區(qū)域?qū)α髟试S尺度模擬研究,發(fā)現(xiàn)模式對(duì)熱帶對(duì)流和遙相關(guān)模擬有很好的效果[73]。將變網(wǎng)格技術(shù)應(yīng)用于CP-RCM的模擬研究也是未來(lái)發(fā)展的趨勢(shì)[64]。計(jì)算資源的大量需求將CP-RCM模擬限制在短期和小區(qū)域中,CP-RCM氣候模擬能夠在更大的區(qū)域和更長(zhǎng)的時(shí)間段上運(yùn)行與高性能科學(xué)計(jì)算的發(fā)展密切相關(guān)。

    猜你喜歡
    嵌套對(duì)流分辨率
    例析“立幾”與“解幾”的嵌套問(wèn)題
    齊口裂腹魚集群行為對(duì)流態(tài)的響應(yīng)
    基于嵌套Logit模型的競(jìng)爭(zhēng)性選址問(wèn)題研究
    EM算法的參數(shù)分辨率
    原生VS最大那些混淆視聽(tīng)的“分辨率”概念
    基于深度特征學(xué)習(xí)的圖像超分辨率重建
    一種改進(jìn)的基于邊緣加強(qiáng)超分辨率算法
    基于ANSYS的自然對(duì)流換熱系數(shù)計(jì)算方法研究
    二元驅(qū)油水界面Marangoni對(duì)流啟動(dòng)殘余油機(jī)理
    一種基于區(qū)分服務(wù)的嵌套隊(duì)列調(diào)度算法
    老师上课跳d突然被开到最大视频| 别揉我奶头 嗯啊视频| 国产精品麻豆人妻色哟哟久久 | 国产熟女欧美一区二区| 色综合亚洲欧美另类图片| 嫩草影院新地址| 性色avwww在线观看| 最近中文字幕高清免费大全6| 欧美日本视频| 真实男女啪啪啪动态图| 少妇的逼水好多| 人妻制服诱惑在线中文字幕| 午夜福利在线观看免费完整高清在| 国产午夜福利久久久久久| 日韩精品青青久久久久久| 久久99蜜桃精品久久| 777米奇影视久久| 亚洲av电影在线观看一区二区三区 | 水蜜桃什么品种好| 久久精品夜色国产| 亚洲精品456在线播放app| 亚洲天堂国产精品一区在线| 国产精品99久久久久久久久| 亚洲av日韩在线播放| 26uuu在线亚洲综合色| 国产成人福利小说| 亚洲av电影在线观看一区二区三区 | 最新中文字幕久久久久| 亚洲熟女精品中文字幕| 五月玫瑰六月丁香| 天美传媒精品一区二区| 18禁裸乳无遮挡免费网站照片| 国产老妇伦熟女老妇高清| 精品一区二区三区视频在线| 色视频www国产| 成年av动漫网址| 亚洲图色成人| 亚洲精品乱码久久久v下载方式| 国产一级毛片在线| 高清在线视频一区二区三区| 亚洲人成网站高清观看| 欧美日本视频| 欧美区成人在线视频| 2022亚洲国产成人精品| 精品久久久久久久久亚洲| 免费电影在线观看免费观看| 男女国产视频网站| 91久久精品国产一区二区成人| 午夜福利视频1000在线观看| 久久精品国产亚洲网站| 国产一级毛片在线| 日本一本二区三区精品| 免费少妇av软件| 免费观看精品视频网站| 日本-黄色视频高清免费观看| 我的女老师完整版在线观看| 在线 av 中文字幕| 国产色婷婷99| 汤姆久久久久久久影院中文字幕 | 国产一级毛片在线| av天堂中文字幕网| 五月玫瑰六月丁香| 日韩 亚洲 欧美在线| 国产欧美日韩精品一区二区| 国产三级在线视频| 超碰av人人做人人爽久久| 天堂俺去俺来也www色官网 | 欧美性猛交╳xxx乱大交人| 99热6这里只有精品| 免费av观看视频| 亚洲精品久久久久久婷婷小说| 国产精品爽爽va在线观看网站| 97精品久久久久久久久久精品| 久久这里有精品视频免费| 午夜福利高清视频| 亚洲av一区综合| 又爽又黄无遮挡网站| av在线亚洲专区| 91aial.com中文字幕在线观看| 久久综合国产亚洲精品| 久久久久久久久久久丰满| 三级男女做爰猛烈吃奶摸视频| 菩萨蛮人人尽说江南好唐韦庄| 午夜福利在线观看吧| 极品少妇高潮喷水抽搐| 日本黄大片高清| 国产av不卡久久| 免费观看的影片在线观看| 亚洲欧美清纯卡通| 国产av在哪里看| 久久精品国产自在天天线| 九九在线视频观看精品| 亚洲第一区二区三区不卡| 久久99蜜桃精品久久| 99热网站在线观看| 国产高清国产精品国产三级 | 网址你懂的国产日韩在线| 欧美日韩视频高清一区二区三区二| 亚洲婷婷狠狠爱综合网| 美女高潮的动态| 国产一区有黄有色的免费视频 | 3wmmmm亚洲av在线观看| 亚洲欧洲日产国产| 午夜免费激情av| 18禁动态无遮挡网站| 午夜免费男女啪啪视频观看| 黄色配什么色好看| 22中文网久久字幕| 噜噜噜噜噜久久久久久91| 国产精品国产三级国产专区5o| 91狼人影院| 精华霜和精华液先用哪个| 69av精品久久久久久| av黄色大香蕉| 日韩av在线免费看完整版不卡| 国产精品精品国产色婷婷| 一级av片app| 白带黄色成豆腐渣| 国产精品人妻久久久影院| 成人av在线播放网站| 国产成人午夜福利电影在线观看| 亚洲国产av新网站| 国产片特级美女逼逼视频| 小蜜桃在线观看免费完整版高清| 男人舔女人下体高潮全视频| 久久久久网色| 成年女人在线观看亚洲视频 | 91狼人影院| 亚洲精品乱码久久久v下载方式| 精品少妇黑人巨大在线播放| 亚洲精品乱码久久久久久按摩| 亚洲一级一片aⅴ在线观看| 色播亚洲综合网| 国产单亲对白刺激| 啦啦啦韩国在线观看视频| 91狼人影院| 在线观看av片永久免费下载| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品久久午夜乱码| 成人性生交大片免费视频hd| av线在线观看网站| 国产一区二区三区综合在线观看 | 国语对白做爰xxxⅹ性视频网站| 久久久久久久大尺度免费视频| 国产午夜精品论理片| 国产av在哪里看| 久久久精品欧美日韩精品| 在线免费十八禁| 亚洲人成网站在线观看播放| 日韩精品有码人妻一区| 水蜜桃什么品种好| .国产精品久久| 日本wwww免费看| 欧美日韩视频高清一区二区三区二| 91久久精品电影网| 亚洲精品成人久久久久久| 国产 亚洲一区二区三区 | 中文精品一卡2卡3卡4更新| 一二三四中文在线观看免费高清| a级毛色黄片| 人人妻人人看人人澡| 亚洲电影在线观看av| 欧美bdsm另类| 天美传媒精品一区二区| 国产精品一二三区在线看| 波多野结衣巨乳人妻| 91久久精品电影网| 亚州av有码| 一级毛片我不卡| 中文字幕人妻熟人妻熟丝袜美| 三级男女做爰猛烈吃奶摸视频| 热99在线观看视频| 久久国内精品自在自线图片| 日日摸夜夜添夜夜添av毛片| 精品一区二区免费观看| 最近的中文字幕免费完整| 亚洲欧美日韩东京热| 久久久久久久久大av| 亚洲国产色片| 国产白丝娇喘喷水9色精品| 国产亚洲91精品色在线| 蜜桃亚洲精品一区二区三区| 日本-黄色视频高清免费观看| 乱人视频在线观看| 午夜精品一区二区三区免费看| 国产三级在线视频| 2021少妇久久久久久久久久久| 久久久欧美国产精品| eeuss影院久久| 在现免费观看毛片| 成人高潮视频无遮挡免费网站| 五月玫瑰六月丁香| 欧美成人午夜免费资源| 久久99蜜桃精品久久| 亚洲va在线va天堂va国产| 免费看不卡的av| 成人国产麻豆网| 亚洲国产精品成人久久小说| 九草在线视频观看| 国产精品日韩av在线免费观看| 亚洲精品国产av蜜桃| 1000部很黄的大片| 好男人视频免费观看在线| 国产激情偷乱视频一区二区| www.色视频.com| 最近的中文字幕免费完整| 久久久久久久久久久免费av| 一级毛片电影观看| 久久精品国产亚洲av天美| ponron亚洲| 久久韩国三级中文字幕| 国产午夜精品论理片| 亚洲精品乱久久久久久| 九九爱精品视频在线观看| 中文天堂在线官网| 18禁裸乳无遮挡免费网站照片| 日日摸夜夜添夜夜添av毛片| 一边亲一边摸免费视频| 国产在线一区二区三区精| 国产黄色免费在线视频| 性色avwww在线观看| 91狼人影院| 爱豆传媒免费全集在线观看| 国产久久久一区二区三区| av网站免费在线观看视频 | 青春草视频在线免费观看| 亚洲精品成人久久久久久| 人妻一区二区av| 免费电影在线观看免费观看| 成年av动漫网址| 亚洲激情五月婷婷啪啪| 99久国产av精品国产电影| 一个人看视频在线观看www免费| 亚洲av福利一区| 晚上一个人看的免费电影| 中文天堂在线官网| 毛片女人毛片| 精品酒店卫生间| av一本久久久久| 国产综合懂色| 国产伦精品一区二区三区视频9| 国产黄色视频一区二区在线观看| 肉色欧美久久久久久久蜜桃 | 91精品国产九色| 男的添女的下面高潮视频| 大片免费播放器 马上看| 国产午夜精品一二区理论片| 久久人人爽人人爽人人片va| videos熟女内射| 国产真实伦视频高清在线观看| 亚洲av免费高清在线观看| 国产一区二区三区综合在线观看 | 亚洲欧美日韩无卡精品| 国产美女午夜福利| av在线蜜桃| 性插视频无遮挡在线免费观看| 亚洲欧美成人精品一区二区| 男女边吃奶边做爰视频| 天堂网av新在线| 久久精品夜夜夜夜夜久久蜜豆| 国产成人aa在线观看| 国产成人freesex在线| 婷婷色麻豆天堂久久| 欧美xxⅹ黑人| 精华霜和精华液先用哪个| 97热精品久久久久久| 国产在线男女| 久久久久久九九精品二区国产| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av天堂中文字幕网| 高清av免费在线| 国产精品美女特级片免费视频播放器| 中文在线观看免费www的网站| 日本黄大片高清| 国产精品女同一区二区软件| 美女cb高潮喷水在线观看| 成人欧美大片| 免费av不卡在线播放| 日韩中字成人| av在线播放精品| 乱系列少妇在线播放| 内地一区二区视频在线| 色综合色国产| 日韩精品青青久久久久久| 国产精品爽爽va在线观看网站| 如何舔出高潮| 精品国产一区二区三区久久久樱花 | 中国国产av一级| 成人二区视频| 在线a可以看的网站| 少妇的逼好多水| 国产精品一二三区在线看| 欧美最新免费一区二区三区| 国产黄a三级三级三级人| 免费看av在线观看网站| 91狼人影院| 免费少妇av软件| 国产午夜精品论理片| 欧美激情久久久久久爽电影| 亚洲伊人久久精品综合| 99久久精品国产国产毛片| 国产在线男女| 精华霜和精华液先用哪个| 国语对白做爰xxxⅹ性视频网站| 可以在线观看毛片的网站| 亚洲成人精品中文字幕电影| 成人午夜高清在线视频| 国产大屁股一区二区在线视频| 成人毛片a级毛片在线播放| 欧美一区二区亚洲| 亚洲不卡免费看| 91久久精品电影网| 一个人免费在线观看电影| 国产色爽女视频免费观看| 久久久午夜欧美精品| 一个人观看的视频www高清免费观看| 国产免费又黄又爽又色| 亚洲欧美精品专区久久| 久久精品久久久久久久性| 日韩av在线免费看完整版不卡| 好男人在线观看高清免费视频| 欧美日韩综合久久久久久| 在线观看免费高清a一片| 国产探花极品一区二区| 综合色丁香网| 搞女人的毛片| 国产黄片视频在线免费观看| 不卡视频在线观看欧美| 深爱激情五月婷婷| 色5月婷婷丁香| 狂野欧美激情性xxxx在线观看| 久久久久久久亚洲中文字幕| 亚洲一级一片aⅴ在线观看| 久久久久久久久久久丰满| 男人舔女人下体高潮全视频| 嘟嘟电影网在线观看| 国产av在哪里看| 午夜激情福利司机影院| 日日啪夜夜撸| 欧美丝袜亚洲另类| 亚洲精品亚洲一区二区| 人妻系列 视频| 久久热精品热| 国产中年淑女户外野战色| 亚洲精品aⅴ在线观看| 亚洲四区av| 精品熟女少妇av免费看| 日韩亚洲欧美综合| 成人毛片60女人毛片免费| 久久国内精品自在自线图片| 三级国产精品欧美在线观看| 午夜爱爱视频在线播放| 亚洲精品视频女| 三级毛片av免费| 午夜福利高清视频| 国产美女午夜福利| 日产精品乱码卡一卡2卡三| 亚洲精品第二区| 国内少妇人妻偷人精品xxx网站| 一本一本综合久久| 久久久久久久久久人人人人人人| 亚洲内射少妇av| 免费电影在线观看免费观看| 亚洲国产欧美人成| 亚洲成人一二三区av| 中文资源天堂在线| 免费看a级黄色片| 啦啦啦中文免费视频观看日本| 国产成人精品一,二区| 91av网一区二区| 美女内射精品一级片tv| 久久久欧美国产精品| 七月丁香在线播放| 亚洲一区高清亚洲精品| 亚洲欧美精品自产自拍| 啦啦啦中文免费视频观看日本| 啦啦啦韩国在线观看视频| av.在线天堂| 天美传媒精品一区二区| 精品国内亚洲2022精品成人| 日本wwww免费看| 国产精品一二三区在线看| 精品不卡国产一区二区三区| 听说在线观看完整版免费高清| 国产成人a区在线观看| 日韩欧美三级三区| 久久精品熟女亚洲av麻豆精品 | 中文乱码字字幕精品一区二区三区 | 青青草视频在线视频观看| 国产午夜精品论理片| 天美传媒精品一区二区| 国产白丝娇喘喷水9色精品| 嫩草影院入口| 天堂中文最新版在线下载 | 一级黄片播放器| 2022亚洲国产成人精品| 91精品国产九色| 久久精品久久久久久噜噜老黄| 精品99又大又爽又粗少妇毛片| 在线观看一区二区三区| 在线观看人妻少妇| 精品久久久久久久久av| 国产伦精品一区二区三区四那| av一本久久久久| 尾随美女入室| 亚洲国产精品成人综合色| 成年版毛片免费区| 最近手机中文字幕大全| 国产永久视频网站| a级毛色黄片| 亚洲国产精品sss在线观看| 亚洲国产成人一精品久久久| 亚洲国产精品专区欧美| 亚洲三级黄色毛片| 听说在线观看完整版免费高清| 精品少妇黑人巨大在线播放| 亚洲婷婷狠狠爱综合网| 亚洲av日韩在线播放| 国产成人aa在线观看| 亚洲精品久久久久久婷婷小说| 美女高潮的动态| 美女主播在线视频| videos熟女内射| 亚洲av福利一区| 亚洲av中文字字幕乱码综合| 观看美女的网站| 白带黄色成豆腐渣| 99久国产av精品国产电影| 男的添女的下面高潮视频| 少妇熟女欧美另类| 欧美日韩精品成人综合77777| 丝瓜视频免费看黄片| 日韩欧美三级三区| 欧美变态另类bdsm刘玥| 一个人免费在线观看电影| 99久久九九国产精品国产免费| 成人无遮挡网站| 国内精品美女久久久久久| 最近的中文字幕免费完整| 亚洲成人精品中文字幕电影| 日本三级黄在线观看| 亚洲av免费在线观看| 日韩一区二区视频免费看| 一级二级三级毛片免费看| 国产成人精品福利久久| 春色校园在线视频观看| 久久精品人妻少妇| 极品教师在线视频| 一级毛片 在线播放| 免费人成在线观看视频色| 舔av片在线| www.av在线官网国产| 一区二区三区免费毛片| 国产高潮美女av| 色视频www国产| 免费播放大片免费观看视频在线观看| 观看免费一级毛片| 97在线视频观看| 亚洲性久久影院| 在线观看免费高清a一片| 97超碰精品成人国产| 日本一本二区三区精品| 国产一区二区在线观看日韩| 99热这里只有是精品在线观看| 午夜福利高清视频| 久久久色成人| 你懂的网址亚洲精品在线观看| 一区二区三区免费毛片| 国产精品一区二区三区四区免费观看| 18禁在线播放成人免费| 国产精品人妻久久久久久| 日本猛色少妇xxxxx猛交久久| 亚洲av成人精品一区久久| 一级毛片 在线播放| 久久精品国产鲁丝片午夜精品| 少妇丰满av| 晚上一个人看的免费电影| 亚洲欧美日韩无卡精品| 国产精品一二三区在线看| 国产黄色免费在线视频| 亚洲人与动物交配视频| 22中文网久久字幕| 亚洲人与动物交配视频| av网站免费在线观看视频 | 精品久久久噜噜| 97超视频在线观看视频| 日日啪夜夜爽| 国产精品人妻久久久影院| 国产成人午夜福利电影在线观看| 国产av国产精品国产| 欧美三级亚洲精品| 国产亚洲一区二区精品| 亚洲av男天堂| 午夜免费男女啪啪视频观看| 毛片一级片免费看久久久久| ponron亚洲| 欧美 日韩 精品 国产| 久久久久久久久久黄片| 国产午夜精品论理片| 亚洲国产最新在线播放| freevideosex欧美| 国产在线一区二区三区精| 国产精品三级大全| 伦精品一区二区三区| videossex国产| 日韩av不卡免费在线播放| 中文字幕亚洲精品专区| 又粗又硬又长又爽又黄的视频| 欧美人与善性xxx| 最后的刺客免费高清国语| 男人舔奶头视频| 亚洲精品乱码久久久v下载方式| 床上黄色一级片| 不卡视频在线观看欧美| 男人和女人高潮做爰伦理| 国产精品无大码| 中文在线观看免费www的网站| 淫秽高清视频在线观看| 波野结衣二区三区在线| 国产成人aa在线观看| 少妇丰满av| 成人欧美大片| 午夜福利视频精品| 久久精品人妻少妇| 色尼玛亚洲综合影院| 国产亚洲午夜精品一区二区久久 | 看非洲黑人一级黄片| 免费少妇av软件| 日本av手机在线免费观看| 国产午夜福利久久久久久| 超碰97精品在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产伦一二天堂av在线观看| 精品久久久精品久久久| 亚洲真实伦在线观看| 老师上课跳d突然被开到最大视频| 久久久精品欧美日韩精品| www.av在线官网国产| 久久99热这里只频精品6学生| 联通29元200g的流量卡| 亚洲国产日韩欧美精品在线观看| 在线免费观看的www视频| 中文字幕人妻熟人妻熟丝袜美| 又大又黄又爽视频免费| 免费看日本二区| 日韩伦理黄色片| 久久精品国产鲁丝片午夜精品| 亚洲欧洲国产日韩| 中文在线观看免费www的网站| 国产 亚洲一区二区三区 | 免费黄色在线免费观看| 亚洲精品自拍成人| 久久99热6这里只有精品| 亚洲自拍偷在线| 亚洲人成网站在线播| 草草在线视频免费看| 国产在视频线在精品| 777米奇影视久久| 精品久久久久久久人妻蜜臀av| 我要看日韩黄色一级片| 精品久久久久久久久久久久久| 一级a做视频免费观看| 免费无遮挡裸体视频| 国产精品美女特级片免费视频播放器| 亚洲怡红院男人天堂| 看黄色毛片网站| 亚洲va在线va天堂va国产| 亚洲成人精品中文字幕电影| 夜夜爽夜夜爽视频| 色播亚洲综合网| 午夜日本视频在线| 观看免费一级毛片| 视频中文字幕在线观看| 1000部很黄的大片| 亚洲精品乱码久久久v下载方式| or卡值多少钱| av.在线天堂| 女的被弄到高潮叫床怎么办| 亚洲欧洲日产国产| 国产老妇女一区| 一级av片app| 国产精品三级大全| 精品不卡国产一区二区三区| 免费看av在线观看网站| 一边亲一边摸免费视频| 国产精品一区二区在线观看99 | 免费观看av网站的网址| 欧美 日韩 精品 国产| 国产精品日韩av在线免费观看| 精品不卡国产一区二区三区| 亚洲精品中文字幕在线视频 | 乱码一卡2卡4卡精品| 天美传媒精品一区二区| 国内少妇人妻偷人精品xxx网站| 精品久久久久久久久久久久久| 观看美女的网站| 午夜福利高清视频| 乱码一卡2卡4卡精品| 91在线精品国自产拍蜜月| 久久精品人妻少妇| 十八禁国产超污无遮挡网站| 女人十人毛片免费观看3o分钟| 亚洲av中文av极速乱| 国产一区有黄有色的免费视频 | 久久久久免费精品人妻一区二区| 在现免费观看毛片| 26uuu在线亚洲综合色| 亚洲欧美成人精品一区二区| 男人和女人高潮做爰伦理| 亚洲欧美成人综合另类久久久| 日日摸夜夜添夜夜添av毛片| 精品欧美国产一区二区三| 亚洲怡红院男人天堂| 国产精品99久久久久久久久| 国产精品美女特级片免费视频播放器|