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

    兩種溫室氣體排放情景下中國汛期江淮暴雨低渦特征研究

    2016-04-16 05:05:31苗春生徐方姝王堅紅余鐘奇高義梅張旭南京信息工程大學氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心氣象災(zāi)害教育部重點實驗室南京00上海浦東新區(qū)氣象局上海00353遼寧省氣象服務(wù)中心沈陽066華風氣象傳媒集團有限公司北京0008
    大氣科學 2016年2期
    關(guān)鍵詞:預(yù)估時段暴雨

    苗春生 徐方姝, 3 王堅紅 余鐘奇 高義梅張旭南京信息工程大學氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心/氣象災(zāi)害教育部重點實驗室,南京00上海浦東新區(qū)氣象局,上海00353遼寧省氣象服務(wù)中心,沈陽066華風氣象傳媒集團有限公司,北京0008

    兩種溫室氣體排放情景下中國汛期江淮暴雨低渦特征研究

    苗春生1徐方姝1, 3王堅紅1余鐘奇2高義梅4張旭4
    1南京信息工程大學氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心/氣象災(zāi)害教育部重點實驗室,南京210044
    2上海浦東新區(qū)氣象局,上海200135
    3遼寧省氣象服務(wù)中心,沈陽110166
    4華風氣象傳媒集團有限公司,北京100081

    苗春生,徐方姝,王堅紅,等. 2016. 兩種溫室氣體排放情景下中國汛期江淮暴雨低渦特征研究 [J]. 大氣科學, 40 (2): 257?270.Miao Chunsheng, Xu Fangshu, Wang Jianhong, et al. 2016. Characteristics of heavy rainfall vortexes during the flood season in Jiang-Huai valley under two greenhouse gas emissions scenarios [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 40 (2): 257?270, doi: 10.3878/j.issn.1006-9895.1505.15109.

    本文基于一個水平分辨率為 50 km 的區(qū)域氣候模式RegCM4(Regional Climate Model,version 4.0)的模擬與預(yù)估結(jié)果,對我國汛期江淮暴雨低渦在氣候變化背景下的統(tǒng)計特征與合成結(jié)構(gòu)進行分析,進一步對兩種溫室排放情景下未來中國汛期的江淮暴雨低渦特征進行預(yù)估。結(jié)果表明:RegCM4模式對環(huán)境要素及低渦都具有一定的模擬能力,低渦的伸展高度、生命期及暴雨位置模擬結(jié)果與觀測較為接近,但模擬的低渦個數(shù)、最大暖區(qū)高度以及溫、濕要素分布均比實際略偏低,而風速和低渦的強度模擬則偏強;在未來兩種溫室排放情景預(yù)估方面,RCP4.5(Representative Concentration Pathways,簡稱RCP)典型濃度排放情景下,暴雨低渦數(shù)量比例減少,強度減弱,但低渦發(fā)展高度仍以850 hPa為主,生命期多為2 d以內(nèi),低渦雨區(qū)分布及最大暖區(qū)高度均與歷史時段相近;RCP8.5情景下,暴雨低渦比例明顯大于RCP4.5情景,低渦發(fā)展高度以700 hPa為主,生命期達3 d的增多,強度增強,最大暖區(qū)厚度范圍顯著伸展。兩種情景下均有低渦中溫度鋒區(qū)減弱,而濕度鋒區(qū)增強,但RCP8.5情景減弱與增強更顯著,顯示更高的溫室氣體排放將導致未來出現(xiàn)更強的暴雨低渦,造成伴隨暴雨的低渦災(zāi)害性天氣的增加,因此應(yīng)進一步深化對低渦暴雨災(zāi)害性天氣發(fā)展趨勢的研究。

    Founded by National Natural Science Foundation of China (Grant 41276033), National Science and Technology Support Program (Grant 2012BAH05B01), Project Supported by the Special Scientific Research Fund of Meteorological Public Welfare Profession of China (Grant GYHY 201206068), Climate Change Special Fund of China Meteorology Administration (Grant CCSF201318), Jiangsu Science and Technology Support Programs (Grants BE2012774, BE2014729), Superiority Discipline Construction Project of Jiangsu Universities and Colleges (PAPD)

    1 引言

    氣候變化是全球變化研究中的重要內(nèi)容和核心問題,它給全球和中國的社會、經(jīng)濟以及環(huán)境帶來重大影響,因此受到各國政府及民眾越來越多的關(guān)注與重視,而氣候模式對于氣候的模擬和預(yù)估,是氣候變化研究中的重要工具之一。由于氣候具有明顯的區(qū)域性特征,而全球環(huán)流模式在模擬區(qū)域氣候方面雖具有一定的模擬能力,但因其分辨率較低尚存在較大的不確定性,因此20世紀90年代以來區(qū)域氣候模式迅速發(fā)展起來且得到了較為廣泛的應(yīng)用,現(xiàn)在已經(jīng)成為氣候研究和業(yè)務(wù)預(yù)報的重要工具(Gao et al., 2001; 鞠麗霞和王會軍,2006; Yu et al., 2010)。相比全球氣候模式,區(qū)域氣候模式對區(qū)域氣候模擬表現(xiàn)出了明顯的優(yōu)勢(Liang et al., 2001;Patricola and Cook, 2007;張冬峰和石英,2012)。李巧萍和丁一匯(2004)的研究結(jié)果表明RegCM(Regional Climate Model)對我國北方長期的干旱氣候態(tài)有一定的模擬能力;很多學者又利用了RegCM2(version 2.0)、RegCM3和RegCM4等區(qū)域氣候模式對我國氣候進行了數(shù)值實驗,研究結(jié)果表明模式對溫度和降水均有一定的模擬能力,能較好地模擬它們的空間分布特征及時間變化趨勢(Gao et al., 2002;高學杰等,2003,2006;石英等,2010;廉麗姝和束炯,2007;翟穎佳和李耀輝,2013);Gao et al.(2013)利用RegCM4嵌套BCC_ CSM1.1模式對21世紀中國地區(qū)氣候變化進行了預(yù)估,結(jié)果表明RCP8.5(Representative Concentration Pathways,簡稱RCP)情景下增溫較RCP4.5更加明顯,而在降水方面兩個模式的結(jié)果則不同。

    可見,已有研究工作大部分都是對區(qū)域氣候模式產(chǎn)品的模擬效果進行評估以及對未來氣候變化進行預(yù)估,主要集中在溫度和降水兩個要素,而在全球氣候變暖的大背景下,中小尺度天氣系統(tǒng)對氣候變化也會有明顯響應(yīng),例如在維持時間、伸展高度、空間范圍、平均強度、發(fā)生頻率、發(fā)生區(qū)域和雨帶配置等方面都有所響應(yīng);而梅汛期暴雨是夏季影響中國長江中下游地區(qū)的主要氣象災(zāi)害之一,由于中尺度渦旋與暴雨密切相關(guān)(Sun et al., 2010),因此中尺度低渦是梅雨暴雨的主要中尺度天氣系統(tǒng)之一(苗春生等,2014a,2014b)。對長江流域梅汛期暴雨及暴雨低渦的研究一直是我國氣象工作者關(guān)心和研究的重點(高坤和徐亞梅,2001;孫建華等,2004;郭蕊等,2013)。如楊引明等(2010)和沈杭鋒等(2013)對長江下游的中尺度低渦活動進行了統(tǒng)計與普查,分析了低渦背景大尺度環(huán)境場及低渦與環(huán)境系統(tǒng)的關(guān)系;張敬萍等(2015)對夏季長江流域兩類中尺度渦旋進行了統(tǒng)計與合成研究,結(jié)果表明,與長江上游西南低渦相比,江淮大別山低渦多發(fā)生在6月,垂直伸展高度更淺薄,引發(fā)降水更強;更多的研究是對江淮暴雨低渦的結(jié)構(gòu)分析與動力機制的診斷(董佩明和趙思雄,2004;李博和趙思雄,2009;周玉淑和李博,2010;傅慎明等,2012;張元春等,2013)。因此結(jié)合區(qū)域氣候模式考慮不同溫室排放情景的影響,分析未來中國汛期江淮暴雨低渦的特征變化是有待開展的研究內(nèi)容。

    因此,本文是基于一個區(qū)域氣候模式RegCM4的氣候變化模擬和預(yù)估結(jié)果(Gao et al., 2013),進行了RCP8.5與RCP4.5兩種溫室排放情景對江淮暴雨低渦特征影響效果的分析研究。在檢驗?zāi)J綄Φ蜏u模擬能力的基礎(chǔ)上,探討未來兩種溫室排放情景下2020~2030年汛期江淮暴雨低渦的特征變化。

    2 資料與方法

    溫室氣體排放情景,是對未來氣候變化預(yù)估的基礎(chǔ)。本文采用新一代情景“典型濃度路徑”中RCP8.5與RCP4.5兩種情景(van Vuuren et al., 2011a,2011b;王紹武等,2012):RCP8.5 情景是2100年輻射強迫上升至8.5 W m?2,這是最高的溫室氣體排放情景。這個情景假定人口多、技術(shù)革新率不高、能源改善緩慢,導致長時間高能源需求及高溫室氣體排放,而缺少應(yīng)對氣候變化的政策;RCP4.5 情景是 2100 年輻射強迫穩(wěn)定在4.5 W m?2??紤]了與全球經(jīng)濟框架相適應(yīng)的,長期存在的全球溫室氣體和生存期短的物質(zhì)的排放,采用低端排放基準和中等減緩措施。

    文中江淮汛期指的是江淮地區(qū)(28°~35°N,111°~125°E)的5~7月,所用模式資料為區(qū)域氣候模式RegCM4產(chǎn)品,水平分辨率為50 km×50 km,垂直方向11層。該模式模擬范圍覆蓋整個中國及周邊地區(qū),模式邊界條件采用單項嵌套BCC_CSM1.1全球氣候系統(tǒng)模式獲得。本研究中所使用的時段分別為1995~2005年(當代),其溫室排放狀況為同時期觀測的溫室氣體狀態(tài);以及2020~2030年(未來),溫室排放狀況為RCP4.5 和RCP8.5兩種情景。用于檢驗?zāi)J綄Ξ敶鷼夂蚰M能力的資料為歐洲中心ERA-interim 1995~2005年的再分析資料(當代觀測)。

    低渦的識別是在研究區(qū)域內(nèi),對歐洲中心ERA-interim再分析資料以及RegCM4模擬產(chǎn)品進行地面至500 hPa各高度層流場分析,對繪制的各層流場圖,依據(jù)是否有閉合的氣旋式環(huán)流,判斷和識別低渦,在識別出的低渦范圍內(nèi)若出現(xiàn)暴雨[>50 mm (24 h)?1],則視為暴雨低渦。低渦中心通過對925~500 hPa各層中閉合流線最大外圈所在層的中心來確定;低渦垂直尺度以閉合環(huán)流發(fā)生的最高高度層來確定,即若500 hPa層上可見閉合環(huán)流則判斷低渦的垂直尺度為500 hPa(深厚);起始識別出渦旋的時間即為低渦出現(xiàn)的時間。

    3 汛期江淮地區(qū)背景要素特征及情景趨勢

    本研究主要分為兩個時段,歷史檢驗時段為1995~2005年,未來情景時段為2020~2030年。為了給出汛期江淮地區(qū)背景要素的定性變化趨勢,本文將溫度、濕度、風速和降水分別對兩個研究時段進行時間平均并在研究區(qū)域內(nèi)進行區(qū)域平均得到以下結(jié)果。

    從圖1a中的溫度圖中可以看出,在歷史時段內(nèi),垂直方向上在150 hPa以下各層,模式模擬的溫度都比觀測溫度偏小,且隨著高度的增加,模擬的溫度與觀測溫度的差距越來越大。而在未來時段內(nèi),兩種情景預(yù)估的平均溫度比模式模擬的歷史時段的平均溫度均有所增加,并且RCP8.5情景下的增溫大于RCP4.5情景下的增溫(圖1b)。在對流層低層水汽含量較大的區(qū)域,模式對歷史時段濕度的模擬與觀測資料的結(jié)果很相近,僅在700 hPa以上模擬結(jié)果較實況偏?。▓D1c)。而從圖1d的偏差圖中可以看出,未來兩種情景下濕度都相對歷史時段有所增加,且RCP8.5情景相較RCP4.5情景增濕更明顯,尤其是在對流層低層。

    從表1中可以看出RegCM4模式對江淮地區(qū)汛期近地面平均風速及總降水量的總體模擬效果及預(yù)估未來的總體趨勢。從平均風速來看,模式在歷史時段對江淮地區(qū)近地面平均風速的模擬結(jié)果是偏大的,而在未來溫室排放增加、溫度升高的情況下,近地面風速也有著相應(yīng)的增大,且RCP8.5情景下的增幅略大于RCP4.5。從降水方面來看,模式在歷史時段對江淮地區(qū)降水的模擬結(jié)果是偏小的,而在未來兩種情景下,降水量較歷史時段的模擬結(jié)果都有所增加,但RCP4.5情景下的模式預(yù)估結(jié)果降水量增大顯著,這與前面幾個要素的變化并不完全一致,說明并不是隨著溫度、濕度以及風速的增加程度越大,降水量的變化量也越大,從而也說明了降水量變化的復(fù)雜性。

    圖1 1995~2005年當代觀測與模式模擬(a)溫度(單位:K)和(c)濕度(單位:g kg?1)的垂直分布,2020~2030年未來兩種情景與歷史時段模擬(b)溫度差值(單位:K)和(d)濕度差值(單位:g kg?1)的垂直分布Fig. 1 Vertical section of observed (ERA-Interim) and simulated (Historical simulation) (a) temperature (units: K) and (c) humidity (units: g kg?1) during 1990–2005; vertical sections of differences between forecasted (2020–2030, under two emissions scenarios) and historical (1995–2005) (b) temperature (units: K) and (d) humidity (units: g kg?1)

    由上可見,歷史時段模式模擬的結(jié)果與觀測較為相近,說明模式對溫度、濕度、風速及降水量具有一定的模擬能力。而對未來兩種情景來說,溫度、濕度及近地面風速均有所增加,且RCP8.5情景相對RCP4.5情景下的增幅更加明顯。而降水在未來兩種情景下均有增加,但增加幅度不同。關(guān)于低渦暴雨增加幅度,以及江淮低渦其它特征變化將在后文第5節(jié)中分析討論。

    表1 不同時段汛期江淮地區(qū)平均降水量與風速Table 1 Average precipitation and average wind velocity in Jiang-Huai valley in different periods

    4 RegCM4模式對歷史汛期江淮暴雨低渦的模擬檢驗

    本節(jié)將對RegCM4模擬結(jié)果和觀測資料(ERA-interim)中1995~2005年汛期江淮地區(qū)低渦的一些基本統(tǒng)計特征以及結(jié)構(gòu)特征分別進行對比分析,從而檢驗RegCM4模式對汛期江淮暴雨低渦的模擬能力。

    4.1統(tǒng)計特征對比

    統(tǒng)計結(jié)果表明,觀測資料在11年中出現(xiàn)了51個低渦,而模式模擬結(jié)果則出現(xiàn)了39個,可見模式模擬出的低渦較觀測偏少。從低渦的發(fā)展高度來看(圖2a),觀測資料和模式模擬結(jié)果中低渦發(fā)展到700 hPa及以下的低渦比例都為84%,而發(fā)展到500 hPa的只有16%。說明雖然模式模擬的低渦個數(shù)比觀測少,但是在低渦發(fā)展高度方面,各高度層低渦個數(shù)占低渦總數(shù)的百分比相差不大,說明江淮地區(qū)汛期中尺度低渦大多數(shù)為淺薄低渦。同時從低渦的生命期來看(圖2b),觀測資料和模式模擬的低渦生命期在2 d以內(nèi)的分別占86% 和87%,說明汛期江淮地區(qū)的低渦系統(tǒng)主要表現(xiàn)為生命期較短的低渦活動。

    圖2 1995~2005 年汛期江淮地區(qū)低渦(a)發(fā)展高度(單位:hPa)、(b)生命期(單位:d),(a、b)中百分比為各項占低渦總數(shù)的比例;低渦暴雨區(qū)中心相對低渦中心的位置和距離示意圖:(c)當代觀測、(d)模式模擬Fig. 2 The (a) vertical development height (units: hPa) and (b) life span (units: d) of the low vortex in Jiang-Huai valley during the historical period (the percentage refers to the proportion relative to the total number of vortexes), and the rainstorm position relative to the center of the low vortex in the (c) observation and (d) simulated results

    1995~2005年汛期,觀測出現(xiàn)暴雨的比例達到49%,而模式模擬結(jié)果中此比例為54%。雖然模式模擬的總低渦數(shù)和暴雨數(shù)都比觀測資料的結(jié)果要小,但是出現(xiàn)暴雨的比例卻略高于觀測資料,而兩種資料的結(jié)果得出的暴雨比例都在50%左右,說明低渦的發(fā)生發(fā)展跟暴雨的關(guān)系是非常密切的。那么就低渦暴雨發(fā)生的位置(圖2c、d)而言,雖然每個低渦暴雨發(fā)生的位置、距離低渦中心的位置不盡相同,但是暴雨區(qū)主要發(fā)生在低渦的東側(cè),尤其是東南側(cè)。統(tǒng)計結(jié)果顯示觀測資料統(tǒng)計出的暴雨區(qū)中心更加密集,貼近于低渦中心,大部分在距離低渦中心400 km的范圍內(nèi),而模擬結(jié)果顯示的暴雨區(qū)中心則相對離散。

    從上面的分析可以看出,低渦的發(fā)生發(fā)展跟暴雨的關(guān)系非常密切,正確分析和認識汛期暴雨低渦結(jié)構(gòu)特征,對于了解低渦暴雨變化趨勢具有重要意義。

    4.2結(jié)構(gòu)特征對比

    為了進一步研究汛期江淮暴雨低渦的平均結(jié)構(gòu)特征,本文將出現(xiàn)暴雨的低渦進行了進一步挑選和動態(tài)合成,以每個渦旋發(fā)展最強盛的日期選為該渦旋合成時刻,并以渦旋中心為合成中心,東西各6個經(jīng)度,南北各6個緯度作為合成范圍,以合成后的渦旋作為11年平均的結(jié)果來討論其平均結(jié)構(gòu)特征。下面對觀測與模擬的兩種資料的低渦合成結(jié)果進行對比分析。

    4.2.1流場和渦度場

    對比兩種資料合成后各層的流場(圖3)可以看出,700 hPa及以下各層上均有明顯的氣旋性環(huán)流中心,對應(yīng)合成得到的中尺度暴雨低渦,且低渦自下而上逐漸向北傾斜,到500 hPa渦旋已消失,對應(yīng)對流層中層為一致的平直西風,由此可見這種汛期暴雨低渦是對流層中低層的系統(tǒng),較淺薄。雖然兩者在850 hPa和925 hPa的結(jié)果較為相似,但在700 hPa上實際合成的低渦緯向范圍較大,而模式模擬的低渦緯向范圍較觀測的偏小。

    圖3 觀測(左列)和模式模擬(右列)的多年平均的(a、b)500 hPa、(c、d)700 hPa及(e、f)850 hPa合成流場,“C”為850hPa層合成低渦中心,下同。橫坐標為距低渦中心的緯向的經(jīng)度(單位:°),縱坐標為距低渦中心的經(jīng)向的緯度(單位:°)Fig. 3 Dynamic synthesis of the stream fields at (a, b) 500 hPa, (c, d) 700 hPa, and (e, f) 850 hPa from the (a, c, e) observation and (b, d, f) simulation results, on average [the ‘C’ indicates the center of the vortex at 850 hPa, the same below; x-axis: zonal longitudes (units: longitudes); y-axis: meridional latitudes (units: latitudes)

    渦度場與流場上氣旋性環(huán)流相配合的有正渦度中心(圖4a、b),兩者在700 hPa及以下低渦中心附近均有正渦度中心,而觀測資料的正渦度中心相對低渦中心更偏東,且模式模擬的低渦強度比觀測的低渦偏強。從垂直分布圖(圖5a)中可以看出合成低渦在對流層低層對應(yīng)的是正渦度區(qū)和輻合區(qū),觀測低渦的低層正渦度最大值出現(xiàn)在700 hPa,而模式模擬的結(jié)果則出現(xiàn)在850 hPa,比觀測偏低。兩者均在400 hPa以上轉(zhuǎn)為負渦度,而觀測低渦的負渦度最大值強于模式模擬低渦的負渦度最大值。

    4.2.2散度和垂直速度場

    對于合成散度和垂直速度場(圖4c、d),在對流層低層700 hPa及以下,輻合區(qū)位于低渦中心東側(cè)及東南側(cè),500 hPa垂直上升運動區(qū)也與輻合中心相對應(yīng);而500 hPa及以上,低渦中心上空及東側(cè)則表現(xiàn)為輻散,200 hPa更強(圖略),這種低層輻合高層輻散的配置以及強上升運動有利于降水的產(chǎn)生,與降水產(chǎn)生在低渦東側(cè)相對應(yīng)。而模式模擬的輻合中心與垂直上升運動大值中心相對觀測結(jié)果都較為離散,這也解釋了模式模擬的暴雨中心相對低渦中心的位置較為離散,距離更遠。

    圖4?。╝、b)合成850 hPa渦度場(單位:10?5s?1),(c、d)850 hPa輻合場(陰影,單位:10?5s?1)和500 hPa垂直速度場(等值線,單位:10?3hPa s?1):(a、c)當代觀測;(b、d)模式模擬Fig. 4 Dynamic synthesis of (a, b) vorticity fields at 850 hPa (units: 10?5s?1) and (c, d) convergence fields at 850 hPa (shading; units: 10?5s?1) and vertical velocity at 500 hPa (contours; units: 10?3hPa s?1): (a, c) Observation; (b, d) simulation

    從散度場的垂直分布(圖5b)來看,觀測資料的無輻散層大約位于400 hPa左右,高層則轉(zhuǎn)為較強的輻散,輻散最大值出現(xiàn)在200 hPa以上;而模式模擬結(jié)果顯示低渦的無輻散層大約位于500 hPa,高層轉(zhuǎn)為強輻散,最大值位于300 hPa左右,且圖中顯示觀測資料的低層輻合略小于模式模擬結(jié)果,而高層輻散情況則相反。

    4.2.3溫度與濕度場

    對溫度場的結(jié)構(gòu)采用對空間平均值求離差的方法來討論(圖6a、b),具體做法是在合成范圍內(nèi),求出通過低渦中心附近的緯向剖面上相應(yīng)經(jīng)向上的平均值,再求出緯向剖面上各點值對此值的偏差。

    圖6顯示,兩種資料結(jié)果都表現(xiàn)為低渦附近700 hPa以下是冷空氣堆,之上到200 hPa都為暖性結(jié)構(gòu),這說明合成的暴雨低渦具有上暖下冷的結(jié)構(gòu),但模擬的低渦暖中心所在高度比觀測偏低。從850 hPa溫濕場的分布來看(圖略),暴雨低渦的東南側(cè)是一東北—西南走向的暖濕氣流帶,而北側(cè)地區(qū)有一干冷空氣自東北向西南;溫度場上,梅雨鋒鋒區(qū)和低壓區(qū)有一條等溫線較密集的弱鋒區(qū),濕度場上則為高濕區(qū)并且南北濕度的水平梯度較大,兩種資料的結(jié)果整體較為相似。

    4.2.4高低空急流

    從高低空急流分布表(表2)中可以看到對流層低層低渦的東南側(cè)有一條東北—西南走向的西南風大風速帶,合成風速在8 m s?1以上,該低空西南急流左側(cè)氣旋式切變對低渦及暴雨的發(fā)生發(fā)展具有重要的增強作用。同時觀測的最大風速達10 m s?1,模擬的最大風速則達到了14 m s?1,比觀測結(jié)果略偏大;在對流層高層存在近乎東西向的極大風速帶,實況低渦中心北側(cè)則存在大于32 m s?1的急流,且大風速范圍比模式模擬結(jié)果向南擴展更大。高空急流右側(cè)的反氣旋輻散流場,有利于低渦中的上升運動以及低渦的維持。

    通過上述分析來看,RegCM4區(qū)域氣候模式在低渦統(tǒng)計和結(jié)構(gòu)特征方面都具有一定的模擬能力;雖然有一些偏差,但在合理范圍內(nèi),因此考慮模式偏差隨時間平移,模式對未來不同情景的預(yù)估數(shù)據(jù)也具有相同的可信度。所以,下文將對預(yù)估數(shù)據(jù)進行分析,了解兩種情景下2020~2030年汛期江淮暴雨低渦的特征。

    圖5 低渦附近(a)渦度和(b)散度區(qū)域平均的垂直分布(單位:10?5s?1)Fig. 5 Vertical profiles of the (a) vorticity and (b) divergence averaged in the domain near the vortex (units: 10?5s?1)

    圖6 溫度偏差緯向垂直剖面圖(單位:K)(陰影區(qū)為正值區(qū)):(a)當代觀測;(b)模式模擬Fig. 6 The zonal vertical cross section of the temperature deviation (units: K)(shaded area is positive): (a) observation; (b) simulation

    表2 低渦附近高低空急流分布情況Table 2 The distribution of the upper-level and low-level jets over the vortexes

    5 未來兩種溫室排放情景下汛期江淮暴雨低渦特征預(yù)估

    5.1統(tǒng)計特征預(yù)估

    利用RegCM4模式未來2020~2030年兩種情景(RCP4.5和RCP8.5)的預(yù)估資料對江淮地區(qū)汛期(5~7月)出現(xiàn)的低渦進行統(tǒng)計,結(jié)果表明:RCP4.5情景下低渦出現(xiàn)46個,而RCP8.5情景下低渦出現(xiàn)51個,相對于模式對歷史時段模擬的低渦個數(shù)都有所增加,且RCP8.5情景則增加的更多,根據(jù)上一章的討論,我們認為模式對未來11年內(nèi)江淮汛期低渦個數(shù)的預(yù)估仍是偏少的。

    從未來低渦的發(fā)展高度來看(圖7a),低渦在700 hPa及以下的個數(shù)仍然達到90%左右,說明未來兩種情景下的低渦仍然是以淺薄系統(tǒng)為主。同歷史時段的模擬結(jié)果對比可知,RCP4.5情景下,925 hPa及850 hPa的低渦所占比例有所增加,尤其是在850 hPa層次上,增長了11%,而在700 hPa和500 hPa層次上低渦所占比例有所下降;但在RCP8.5情景下,只有700 hPa層次低渦比例有所增加,且增幅達到14%,而在850 hPa和925 hPa低渦比例明顯減少,500 hPa層次低渦比例保持不變。說明在RCP4.5情景下,低渦發(fā)展高度以850 hPa為主,而在RCP8.5情景下低渦發(fā)展高度則相對較高,以700 hPa為主。而從圖7b來看,RCP4.5情景下低渦生命期在3 d及以上的比例同歷史時段模式模擬的低渦生命期相同,說明在這種排放情景下,低渦仍以2 d內(nèi)的短生命期為主;而在RCP8.5排放情景下,低渦生命期達到3 d的比例明顯增加,從歷史時段的5%增加到了27.5%,從而說明RCP8.5排放情景下低渦的生命期有一定程度的延長。

    圖7 2020~2030年汛期江淮地區(qū)低渦(a)發(fā)展高度(單位:hPa)、(b)生命期(單位:d),(a、b)中百分比為各項占低渦總數(shù)的比例。低渦暴雨區(qū)中心相對低渦中心的位置和距離示意圖:(c)RCP4.5;(d)RCP8.5Fig. 7 (a) The vertical development height (units: hPa) and (b) the life span (units: d) of the low vortexes in the Jiang-Huai valley during the years 2020–2030 (the percentage is the proportion relative to the total number of low vortexes), and the rainstorm position relative to the center of the low vortexes under emissions scenarios (c) RCP4.5 and (d) RCP8.5

    從低渦暴雨方面來看,雖然未來情景下低渦的總數(shù)有所增加,但是出現(xiàn)暴雨的低渦比例卻有所減少,尤其是RCP4.5情景下,低渦暴雨比例較歷史模擬結(jié)果減少了一半以上,說明隨著溫室氣體排放量的增加,低渦暴雨出現(xiàn)的比例并沒有增加。再從暴雨位置的分布(圖7c、d)來看,RCP4.5情景下,暴雨低渦出現(xiàn)的個數(shù)雖然很少,但暴雨落區(qū)仍然在低渦中心的東部尤其是東南部,而其暴雨中心出現(xiàn)的位置非常離散,且距低渦中心400 km以外的低渦達50%,最遠可達800 km;在RCP8.5情景下,除個別暴雨位于低渦西南部,其它暴雨落區(qū)都在低渦中心的東部,而暴雨中心的分布相對于歷史時段模擬結(jié)果和RCP4.5結(jié)果都更加緊湊集中,且大部分都位于距低渦中心400 km以內(nèi)的位置,暴雨低渦數(shù)量較RCP4.5增多。

    5.2暴雨低渦的結(jié)構(gòu)特征預(yù)估

    為了進一步預(yù)估未來情景下汛期江淮暴雨低渦的結(jié)構(gòu)特征,將兩種情景下的暴雨低渦同樣進行動態(tài)合成。

    5.2.1流場和渦度場

    兩種情景下的流場特征如圖8所示,兩種情景的流場圖差別較大,RCP4.5情景下,只有在850 hPa才有明顯的閉合渦旋,而RCP8.5情景下在850 hPa 和700 hPa上都有明顯的閉合渦旋,且低渦的范圍較RCP4.5排放情景更大。到達500 hPa層次渦旋已經(jīng)消散,在RCP4.5情景下表現(xiàn)為近乎平直的西風,而RCP8.5情景下則表現(xiàn)為弱槽結(jié)構(gòu),因此未來情景下的低渦仍為對流層中低層的結(jié)構(gòu)。垂直方向上,低渦自下而上向北傾斜,且RCP4.5情景下的傾斜程度比RCP8.5情景下的傾斜程度更強。與歷史時段的模擬結(jié)果相比,RCP4.5情景下的低渦水平尺度略小,而RCP8.5情景下則略有增加。

    圖8 同圖3,但為RCP4.5(左列)和RCP8.5(右列)Fig. 8 Same as Fig. 3, but for emissions scenarios RCP4.5 (left column) and RCP8.5 (right column)

    從合成渦度場(圖9a、b)上來看,未來兩種情景下與低渦中心相對應(yīng)的正渦度中心位置相近,但強度差別較大,RCP4.5情景下的渦度場強度相對歷史模擬結(jié)果有所減弱,而高排放RCP8.5情景下則顯著增強,說明RCP8.5情景下的低渦強度比RCP4.5情景下的低渦強度更強。而在垂直方向上(圖5a中方形連線),RCP4.5情景下(空心方形線)400 hPa以下均為正渦度區(qū),而RCP8.5情景(實心方形線)正渦度區(qū)則達到了300 hPa,兩種情景下的正渦度最大值都出現(xiàn)在850~700 hPa,這與歷史時段的模擬結(jié)果相同,但RCP4.5情景下各層渦度值均小于模式對歷史時段的模擬結(jié)果,而RCP8.5情景在200 hPa以下的層次渦度值均大于模式對歷史時段的模擬結(jié)果,說明在垂直方向上,RCP4.5情景下低層正渦度區(qū)的強度有所減弱,但是高層負渦度區(qū)的強度有所增強,而RCP8.5情景的情況則剛好相反。

    5.2.2散度和垂直速度場

    如圖9c、d所示,兩種情景下對流層低層低渦中心及其東部和南部為強輻合區(qū)域,強輻合中心位于低渦中心的東部,而在高層200 hPa低渦中心附近則為強輻散區(qū),同時有垂直上升運動在此區(qū)與其配合,這種配置同歷史模擬的情況相似,同樣在低渦中心的東部有利于降水的產(chǎn)生。而RCP4.5情景下的低層輻合區(qū)域范圍明顯小于RCP8.5情景下低層輻合區(qū),且輻合強度也略小于RCP8.5情景(圖略)。兩者的垂直上升運動最大值較為接近,達到-8×10?3hPa s?1,比模擬的歷史結(jié)果有所增大。而且RCP8.5情景下散度大值中心的分布與歷史模擬的結(jié)果相比相對集中,因此降水中心分布也較歷史時段稍集中些。在垂直方向上(圖5b中方形連線),RCP8.5情景下(實心方形線)的低渦附近散度區(qū)域平均后的垂直分布與歷史模擬的結(jié)果極為相近,而RCP4.5情景下(空心方形線)則略有不同,表現(xiàn)為500 hPa以下輻合偏弱,400 hPa輻散偏強,之上近乎相同。

    圖9 同圖4,但為(a、c)RCP4.5和(b、d)RCP8.5Fig. 9 Same as Fig. 4, but for emissions scenarios (a, c) RCP4.5 and (b, d) RCP8.5

    5.2.3溫、濕場

    溫度離差與歷史模擬的結(jié)果相似(圖10a、b),兩種情景下700 hPa以下是冷空氣堆,之上到250 hPa都為暖性結(jié)構(gòu),這說明未來情景下合成的低渦仍然具有上暖下冷的結(jié)構(gòu)。相比歷史的模擬結(jié)果,RCP4.5情景下最大暖區(qū)高度與歷史相近,但強度增強;而RCP8.5情景下的最大暖區(qū)則在700~400 hPa之間,最大暖區(qū)厚度范圍顯著伸展。從溫濕場來看(圖略),兩種情景下低渦東南側(cè)為暖濕氣流帶,而北側(cè)則為干冷空氣。與歷史模擬結(jié)果對比,未來情景下的對流層低渦附近溫度梯度減小,濕度梯度增大,且濕度比歷史時段的濕度也有所增大,RCP8.5情景比RCP4.5情景增大的更多。

    圖10 未來兩種情景下溫度偏差緯向垂直剖面圖(單位:K)(陰影區(qū)為正值區(qū)):(a)RCP4.5;(b)RCP8.5Fig. 10 The zonal vertical cross sections of the temperature deviation (units: K) (shaded area is positive); under emissions scenarios (a) RCP4.5 and (b) RCP8.5

    5.2.4高低空急流

    從高低空急流分布表(表2下面兩行)來看,兩種情景下低渦東南側(cè)仍為低空西南風急流,由于模式對低空急流模擬偏大,未來低空急流的增幅應(yīng)弱于預(yù)估值。而從200 hPa高空風速場上來看,高空急流仍然在低渦北側(cè),但是模式對高空急流模擬偏小,因此未來高空急流的增幅應(yīng)強于預(yù)估值??紤]高低空急流的總效果,RCP8.5情景下低空西南風急流強于RCP4.5情景,更有利于淺薄渦旋的增強。對于高空急流的預(yù)估,無論RCP4.5情景還是RCP8.5情景,高空西風急流增幅有限,對低渦增強的作用基本維持類似歷史狀況。

    顯然,高溫室氣體排放下有低渦發(fā)生發(fā)展數(shù)量的增加,但是在RCP4.5情景下伴隨暴雨的低渦比例有明顯的減少,而RCP8.5情景下伴隨暴雨的低渦比例減幅較小。值得注意的是,更高的溫室氣體排放將促使伴隨暴雨的低渦發(fā)展更劇烈,因此對低渦暴雨災(zāi)害性天氣發(fā)展趨勢的研究有必要進一步深入,對未來情景下災(zāi)害性低渦暴雨預(yù)報提供一定的參考。

    6 結(jié)論與討論

    本文基于區(qū)域氣候模式RegCM4的模擬與預(yù)估結(jié)果,對1995~2005年汛期中國江淮地區(qū)的低渦進行氣候變化的模擬檢驗與未來兩種情景下2020~2030年的預(yù)估,旨在給出兩種溫室氣體排放情景下汛期江淮暴雨低渦的特征變化及其對比。結(jié)果表明:

    (1)RegCM4模式對汛期江淮地區(qū)平均溫度、濕度及風速等要素的模擬特點為:模擬的溫、濕要素均比觀測略微偏低一些,且隨高度增加偏差有所增大,而對風速的模擬則比觀測偏大;預(yù)估結(jié)果顯示未來兩個情景下各要素值均增大,且溫濕風的增大程度在RCP8.5情景下比RCP4.5情景更大,汛期平均總降水量有增加,但RCP4.5情景較RCP8.5更大些。

    (2)模式對汛期江淮低渦的模擬特點為:無論從低渦統(tǒng)計方面還是低渦結(jié)構(gòu)方面與觀測均較為相似。雖然模擬的低渦總數(shù)偏少,但是其統(tǒng)計特征的各項比例和暴雨低渦的結(jié)構(gòu)及診斷量都相似與接近,因此其模式模擬的暴雨低渦具有一定的可信度,可作為未來情景分析基礎(chǔ)。

    (3)未來兩種情景下低渦統(tǒng)計特征為:個數(shù)都有所增加,但是產(chǎn)生暴雨的低渦比例卻有所減少。RCP4.5情景下,暴雨落區(qū)仍然在低渦中心的東部尤其是東南部,但暴雨中心出現(xiàn)的位置非常離散,且距低渦中心400 km以外的低渦達50%。RCP8.5情景下,除個別低渦位于低渦西南部,其它低渦暴雨落區(qū)也都在低渦中心的東部,而暴雨中心的分布相對于歷史時段模擬結(jié)果和RCP4.5結(jié)果都更加緊湊集中,且大部分都位于距低渦中心400 km以內(nèi)的位置。

    (4)未來兩種情景下低渦結(jié)構(gòu)特征為:RCP4.5情景下,合成低渦發(fā)展高度以850 hPa為主,生命期多為2 d以內(nèi)短生命期,強度減弱。而RCP8.5情景下,暴雨低渦比例明顯大于RCP4.5情景,低渦發(fā)展高度以700 hPa為主,較RCP4.5情景升高;生命期達3 d的增多,更長;低渦渦度較RCP4.5強度明顯增強,其最大暖區(qū)厚度范圍顯著伸展。溫度鋒區(qū)減小而濕度鋒區(qū)進一步增強,近地面低渦低空西南風急流強于RCP4.5情景,更有利于淺薄渦旋的增強。高空西風急流增幅有限,對低渦增強的作用基本維持類似歷史狀況。因此RCP8.5較RCP4.5對暴雨低渦有更強的促發(fā)和增強效應(yīng)。

    綜上所述,在考慮模式不確定性基礎(chǔ)上,RegCM4模式的結(jié)果,具有一定的合理性及預(yù)測性,能夠定性的給出未來情景下低渦的變化特征與變化趨勢,高排放將促使伴隨暴雨的淺薄低渦發(fā)展更劇烈,因此應(yīng)進一步深化對低渦暴雨災(zāi)害性天氣發(fā)展趨勢的研究,對未來情景下災(zāi)害性低渦暴雨預(yù)報提供一定的參考。此外,模式的近期預(yù)估結(jié)果受自然變率影響較大,自然變率信號可能會部分掩蓋溫室氣體強迫,因此后續(xù)的工作將針對更長氣候時段進行溫室氣體排放對江淮暴雨低渦特征演變的影響研究。

    致謝感謝國家氣候中心提供的利用RegCM4區(qū)域氣候模式所進行的中國區(qū)域未來氣候變化模擬結(jié)果。

    (References)

    董佩明, 趙思雄. 2004. 引發(fā)梅雨鋒暴雨的頻發(fā)型中尺度低壓(擾動)的診斷研究 [J]. 大氣科學, 28 (6): 876–891.Dong Peiming, Zhao Sixiong. 2004. A diagnostic study of mesoscale lows (disturbances) on Meiyu front and associated heavy rainfall [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 28 (6): 876–891, doi:10.3878/j.issn. 1006-9895.2004.06.07.

    傅慎明, 于翡, 王東海, 等. 2012. 2010年梅雨期兩類東移中尺度渦旋的對比研究 [J]. 中國科學: 地球科學, 42 (8): 1282–1300.Fu Shenming, Yu Fei, Wang Donghai, et al. 2013. A comparison of two kinds of eastward-moving mesoscale vortices during the Mei-yu period of 2010 [J]. Science China Earth Sciences, 56 (2): 282–300, doi:10.1007/s11430-012-4420-5.

    高坤, 徐亞梅. 2001. 1999年6月下旬長江中下游梅雨鋒低渦擾動的結(jié)構(gòu)研究 [J]. 大氣科學, 25 (6): 740–756.Gao Kun, Xu Yamei. 2001. A simulation study of structure of mesovortexes along Meiyu front during 22–30 June 1999 [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 25 (6): 740–756, doi:10.3878/j.issn.1006-9895.2001.06.02.

    高學杰, 趙宗慈, 丁一匯, 等. 2003. 溫室效應(yīng)引起的中國區(qū)域氣候變化的數(shù)值模擬 II: 中國區(qū)域氣候的可能變化 [J]. 氣象學報, 61 (1): 29–38.Gao Xuejie, Zhao Zongci, Ding Yihui, et al. 2003. Climate change due to greenhouse effects in China as simulated by a regional climate model. Part II: Climate change [J]. Acta Meteorologica Sinica (in Chinese), 61 (1): 29–38, doi:10.11676/qxxb2003.004.

    高學杰, 徐影, 趙宗慈, 等. 2006. 數(shù)值模式不同分辨率和地形對東亞降水模擬影響的試驗 [J]. 大氣科學, 30 (2): 185–192.Gao Xuejie, Xu Ying, Zhao Zongci, et al. 2006. Impacts of horizontal resolution and topography on the numerical simulation of East Asian precipitation [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 30 (2): 185–192, doi:10.3878/j.issn.1006-9895.2006.02.01.

    Gao X J, Zhao Z C, Ding Y H, et al. 2001. Climate Change due to greenhouse effects in China as simulated by a regional climate model [J]. Adv. Atmos. Sci., 18 (6): 1224–1230, doi:10.1007/s00376-001-0036-y.

    Gao X J, Zhao Z C, Giorgi F. 2002. Changes of extreme events in regional climate simulations over East Asia [J]. Adv. Atmos. Sci., 19 (5): 927–942, doi:10.1007/s00376-002-0056-2.

    Gao X J, Wang M L, GIORGI F. 2013. Climate change over China in the 21st century as simulated by BCC_CSM1.1-RegCM4.0 [J]. Atmos. Oceanic Sci. Lett., 6 (5): 381–386, doi:10.3878/j.issn.1674-2834.13.0029.

    郭蕊, 苗春生, 張楠. 2013. 一次淮河流域梅雨鋒暴雨的大別山地形敏感性試驗 [J]. 大氣科學學報, 36 (5): 626–634.Guo Rui, Miao Chunsheng, Zhang Nan. 2013. Sensitivity experiments of effects of Dabie mountains terrain on Meiyu front rainstorm over Huaihe River basin [J]. Trans. Atmos. Sci. (in Chinese), 36 (5): 626–634, doi:10.3969/j.issn. 1674-7097.2013.05.014.

    鞠麗霞, 王會軍. 2006. 用全球大氣環(huán)流模式嵌套區(qū)域氣候模式模擬東亞現(xiàn)代氣候 [J]. 地球物理學報, 49: 52–60.Ju Lixia, Wang Huijun. 2006. Modern climate over East Asia simulated by a regional climate model nested in a global gridpoint general circulation model [J]. Chinese J. Geophys. (in Chinese), 2006, 49 (1): 52–60, doi:10.3321/j.issn:0001-5733.2006.01.008.

    李博, 趙思雄. 2009. 2007年入梅期由橫槽與低渦切變引發(fā)淮河流域強降水的診斷研究 [J]. 大氣科學, 33 (6): 1148–1164.Li Bo, Zhao Sixiong. 2009. A study of heavy rainfall resulting from transversal trough and shear line-low vortex in the Huaihe River during the Meiyu onset of 2007 [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 33 (6): 1148–1164, doi:10.3878/j.issn.1006-9895.2009.06.03.

    李巧萍, 丁一匯. 2004. 區(qū)域氣候模式對東亞季風和中國降水的多年模擬與性能檢驗 [J]. 氣象學報, 62 (2): 140–153.Li Qiaoping, Ding Yihui. 2004. Multi-year simulation of the East Asian monsoon and precipitation in China using a regional climate model and evaluation [J]. Acta Meteorologica Sinica (in Chinese), 62 (2): 140–153, doi:10.11676/ qxxb2004.016.

    廉麗姝, 束炯. 2007. 區(qū)域氣候模式對我國中、東部夏季氣候的數(shù)值模擬[J]. 熱帶氣象報, 23 (2): 162–170.Lian Lishu, Shu Jiong. 2007. Numerical simulation of summer climate over center and East China using a regional climate model [J]. J. Trop. Meteor. (in Chinese), 23 (2): 162–170.

    Liang X Z, Kunkel K E, Samel A N. 2001. Development of a regional climate model for U.S. Midwest applications. Part I: Sensitivity to buffer zone treatment [J]. J. Climate, 14 (23): 4363–4378, doi:10.1175/1520-0442(2001)014<4363:DOARCM>2.0.CO;2.

    苗春生, 劉維鑫, 王堅紅, 等 2014a. 梅雨期經(jīng)大別山兩側(cè)暴雨中尺度低渦對比分析 [J]. 高原氣象, 33 (2): 394–406.Miao Chunsheng, Liu Weixin, Wang Jianhong. 2014b. Comparison between two rainstorm meso-scale vortices circumambulated respectively from two sides ofmountain Dabieshan during Meiyu season [J]. Plateau Meteor. (in Chinese), 33 (2): 394–406, doi:10.7522/j.issn.1000-0534.2013.00195.

    苗春生, 吳旻, 王堅紅, 等. 2014b. 一次淺薄低渦暴雨過程數(shù)值模擬及發(fā)展機制分析 [J]. 氣象, 40 (1): 28–37.Miao Chunsheng, Wu Min, Wang Jianhong, et al. 2014a. Modeling of a shallow vortex heavy rainfall and analysis on its developing mechanism [J]. Meteor. Mon. (in Chinese), 40 (1): 28–37, doi:10.7519/j.issn.1000-0526.2014.01.004.

    Patricola C M, Cook K H. 2007. Dynamics of the West African monsoon under mid-Holocene precessional forcing: Regional climate model simulations [J]. J. Climate, 20(4): 694–716, doi:10.1175/JCLI4013.1.

    沈杭鋒, 翟國慶, 尹金方, 等. 2013. 長江下游梅汛期中尺度渦旋特征分析 [J]. 大氣科學, 37 (4): 923–932.Shen Hangfeng, Zhai Guoqing, Yin Jinfang, et al. 2013. Feature analysis of mesoscale vortex over lower reaches of Yangtze River during Meiyu period [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 37 (4): 923–932, doi:10.3878/j.issn. 1006-9895.2012.12075.

    石英, 高學杰, 吳佳, 等. 2010. 華北地區(qū)未來氣候變化的高分辨率數(shù)值模擬 [J]. 應(yīng)用氣象學報, 21 (5): 580–589.Shi Ying, Gao Xuejie, Wu Jia, et al. 2010. Simulating future climate changes over North China with a high resolution regional climate model [J]. Journal of Applied Meteorological Science (in Chinese), 21 (5): 580–589, doi:10.11898/ 1001-7313.20100507.

    孫建華, 張小玲, 齊琳琳, 等. 2004. 2002年6月20~24日梅雨鋒中尺度對流系統(tǒng)發(fā)生發(fā)展分析 [J]. 氣象學報, 62 (4): 423–438.Sun Jianhua, Zhang Xiaoling, Qi Linlin, et al. 2004. An analysis on MCSs in Meiyu front during 20–24 June 2002 [J]. Acta Meteor. Sinica (in Chinese), 62 (4): 423–438, doi:10.11676/qxxb2004.043.

    Sun J H, Zhao S X, Xu G K, et al. 2010. Study on a mesoscale convective vortex causing heavy rainy rainfall during the Mei-yu season in 2003 [J]. Adv. Atmos. Sci., 27 (5): 1195–1209, doi:10.1007/s00376-009-9156-6.

    van Vuuren D P, Edmonds J A, Kainuma M, et al. 2011a. A special issue on the RCPs [J]. Climatic Change, 109: 1–4, doi:10.1007/s10584-011-0157-y.

    van Vuuren D P, Edmonds J A, Kainuma M, et al. 2011b. The representative concentration pathways: An overview [J]. Climatic Change, 109: 5–31, doi:10.1007/s10584-011-0148-z.

    王紹武, 羅勇, 趙宗慈, 等. 2012. 新一代溫室氣體排放情景 [J]. 氣候變化研究進展, 8 (4): 305–307.Wang Shaowu, Luo Yong, Zhao Zongci, et al. 2012. New generation of scenarios of greenhouse gas emission [J]. Progressus Inquisitiones de Mutatione Climatis (in Chinese), 8 (4): 305–307, doi:10.3969/j.issn.1673-1719.2012.04.011.

    楊引明, 谷文龍, 趙銳磊, 等. 2010. 長江下游梅雨期低渦統(tǒng)計分析 [J].應(yīng)用氣象學報, 21 (1): 11–18.Yang Yinming, Gu Wenlong, Zhao Ruilei, et al. 2010. The statistical analysis of low vortex during Meiyu season in the lower reaches of the Yangtze [J]. J. Appl. Meteor. Sci. (in Chinese), 21 (1): 11–18, doi:10.3969/j.issn.1001-7313.2010.01.002.

    Yu E T, Wang H J, Sun J Q. 2010. A quick report on a dynamical downscaling simulation over china using the nested model [J]. Atmos. Ocean Sci. Lett., 3: 325–329, doi:10.1080/16742834.2010.11446886.

    翟穎佳, 李耀輝. 2013. 區(qū)域氣候模式RegCM4對我國西北東部到華北地區(qū)21世紀中期氣候變化的情景預(yù)估 [C]. 南京: 創(chuàng)新驅(qū)動發(fā)展提高氣象災(zāi)害防御能力——S18第四屆研究生年會. Zhai Yingjia, Li Yaohui. 2013. Scenario prediction on the climate changes from the eastern part of Northwest China to North China during the mid 21st century by regional climate model RegCM4 [C]. Nanjing: The 4th Annual Meeting of Postgraduates, S8: Innovate, Drive, Develop and Improve the Ability of Defending Meteorological Disasters.

    張冬峰, 石英. 2012. 區(qū)域氣候模式RegCM3對華北地區(qū)未來氣候變化的數(shù)值模擬 [J]. 地球物理學報, 55 (9): 2854–2866.Zhang Dongfeng, Shi Ying. 2012. Numerical simulation of climate changes over North China by the RegCM3 model [J]. Chinese J. Geophys. (in Chinese), 55 (9): 2854–2866, doi:10.6038/j.issn.0001-5733.2012.09.005.

    張敬萍, 傅慎明, 孫建華, 等. 2015. 夏季長江流域兩類中尺度渦旋的統(tǒng)計與合成研究 [J]. 氣候與環(huán)境研究, 23 (3): 319–336.Zhang Jingping, Fu Shenming, Sun Jianhua, et al. 2015. A statistical and compositional study on the two types of mesoscale vortices over the Yangtze River basin [J]. Climatic and Environmental Research (in Chinese), 23 (3): 319–336, doi:10.3878/j.issn.1006-9585.2015.14164.

    張元春, 孫建華, 徐廣闊, 等. 2013. 江淮流域兩次中尺度對流渦旋的結(jié)構(gòu)特征研究 [J]. 氣候與環(huán)境研究, 18 (3): 271–287.Zhang Yuanchun, Sun Jianhua, Xu Guangkuo, et al. 2013. Analysis on the structure of two mesoscale convective vortices over Yangtze–Huaihe River basin [J]. Climatic and Environmental Research (in Chinese), 18 (3): 271–287, doi:10.3878/j.issn.1006-9585.2012.11162.

    周玉淑, 李柏. 2010. 2003年7月8~9日江淮流域暴雨過程中渦旋的結(jié)構(gòu)特征分析 [J]. 大氣科學, 34 (3): 629–639.Zhou Yushu, Li Bai. 2010. Structural analyses of vortex causing torrential rain over the Changjiang–Huaihe River basin during 8 and 9 July 2003 [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 34 (3): 629–639, doi:10.3878/j.issn.1006-9895.2010.03.14.

    Characteristics of Heavy Rainfall Vortexes during the Flood Season in Jiang-Huai Valley under Two Greenhouse Gas Emission Scenarios

    MIAO Chunsheng1, XU Fangshu1, 3, WANG Jianhong1, YU Zhongqi2, GAO Yimei4, ZHANG Xu4
    1 Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters/Key Laboratory of Meteorological Disaster of Ministry of Education, Nanjing University of Information Science & Technology, Nanjing 210044

    2 Pudong New Area Meteorological Bureau of Shanghai Meteorological Service, Shanghai 200135
    3 Liaoning Province Meteorological Service Center, Shenyang 110166
    4 Huafeng Meteorological Media Group, Beijing 100081

    The aim of this study was to analyze the statistical and structure evolutive characteristics of heavy rainfall vortexes during the flood season in Jiang-Huai valley, China, against the background of climate change, and to provide projections under two emission scenarios for the period 2020–2030, by using the regional climate model RegCM4 (Regional Climate Model, version 4.0) with a horizontal resolution of 50 km. The results showed that RegCM4 can capture the environmental aspects and characteristics of this vortex type reasonably well. The simulated outputs were similar to reality in several ways, including vertical development height, lifetime, and the relative rainstorm position of the vortex. However, the simulated results were lower with respects to the number and the warmest section height of the vortex, and the background temperature and humidity, while the intensity of the vortex and velocity were higher. In terms of projection for the period 2020–2030, under the RCP4.5 (Representative Concentration Pathways) scenario, the proportion of vortexes with rainstorms was shown to decrease; the vertical development height of the vortex was mainly at 850 hPa; and the life span was mostly within 2 days. The distribution of heavy rainfall and the warmest zone height were close to those in the historical period, but the intensity was weaker. Meanwhile, under the RCP8.5 scenario, the percentage of vortexes with rainstorms was projected to be greater than that under the RCP4.5 scenario; the vertical development height of the vortex was mainly at 700 hPa; the vortexes with the life span of 3 days increased; the vortex intensity enhanced; and the thickness of the warmest zone extended significantly. The moisture fronts in the vortex strengthened, but temperature fronts weakened. The study reveals that more severe vortex precipitation may occur in the future under higher greenhouse gas emissions. Further research is therefore needed on the development trends of low vortexes with severe rainfall.

    Model RegCM4 (Regional Climate Model, version 4.0), Greenhouse gas, Emission scenarios, Flood season, Jiang-Huai valley, Heavy rainfall vortex

    國家自然科學基金項目41276033,國家科技支撐項目2012BAH05B01,公益性行業(yè)(氣象)專項項目 GYHY201206068,中國氣象局氣候變化專項江蘇氣候變化評估CCSF201318,江蘇科技支撐項目BE2012774、BE2014729,江蘇高校優(yōu)勢學科建設(shè)工程資助項目(PAPD)

    RegCM4(Regional Climate Model, version 4.0)模式典型濃度溫室排放情景江淮汛期暴雨低渦

    1006-9895(2016)02-0257-14

    P447 P466

    A

    10.3878/j.issn.1006-9895.1505.15109

    2015-01-13;網(wǎng)絡(luò)預(yù)出版日期 2015-05-15

    苗春生,男,1954年出生,教授,博士生導師,主要從事大氣中小尺度動力學。E-mail: 1597706505@qq.com, csmiao@nuist.edu.cn

    猜你喜歡
    預(yù)估時段暴雨
    美國銀行下調(diào)今明兩年基本金屬價格預(yù)估
    “80年未遇暴雨”襲首爾
    暴雨
    當暴雨突臨
    四個養(yǎng)生黃金時段,你抓住了嗎
    暴雨襲擊
    支點(2017年8期)2017-08-22 17:18:27
    史密斯預(yù)估控制在排焦控制中的應(yīng)用
    傍晚是交通事故高發(fā)時段
    分時段預(yù)約在PICC門診維護中的應(yīng)用與探討
    分時段預(yù)約掛號的實現(xiàn)與應(yīng)用
    国产一区二区三区在线臀色熟女| 日本在线视频免费播放| 两人在一起打扑克的视频| 国产97色在线日韩免费| 色av中文字幕| 久久伊人香网站| 国产aⅴ精品一区二区三区波| 国产精品 欧美亚洲| 国产99久久九九免费精品| 窝窝影院91人妻| 丁香六月欧美| 好看av亚洲va欧美ⅴa在| 激情在线观看视频在线高清| 又紧又爽又黄一区二区| 悠悠久久av| 人人妻,人人澡人人爽秒播| 日韩 欧美 亚洲 中文字幕| 他把我摸到了高潮在线观看| 国产私拍福利视频在线观看| 哪里可以看免费的av片| or卡值多少钱| 久久精品91无色码中文字幕| 99在线视频只有这里精品首页| 亚洲成国产人片在线观看| 丰满人妻熟妇乱又伦精品不卡| 两个人视频免费观看高清| 黑人巨大精品欧美一区二区mp4| 91成人精品电影| 国产精品综合久久久久久久免费| 巨乳人妻的诱惑在线观看| 男女床上黄色一级片免费看| 日韩欧美免费精品| 又大又爽又粗| 国产成人一区二区三区免费视频网站| 婷婷丁香在线五月| 精品国内亚洲2022精品成人| 又黄又粗又硬又大视频| 女同久久另类99精品国产91| 亚洲美女黄片视频| 日韩三级视频一区二区三区| 亚洲五月天丁香| 两人在一起打扑克的视频| 精品国产乱子伦一区二区三区| 亚洲专区字幕在线| 色哟哟哟哟哟哟| 草草在线视频免费看| 国内毛片毛片毛片毛片毛片| 国产av又大| 欧美中文综合在线视频| 看免费av毛片| 亚洲专区字幕在线| 久久国产亚洲av麻豆专区| 欧美国产精品va在线观看不卡| 久久精品国产综合久久久| 亚洲成人免费电影在线观看| 国产在线观看jvid| 女生性感内裤真人,穿戴方法视频| 中国美女看黄片| 日韩一卡2卡3卡4卡2021年| 精品第一国产精品| 又黄又爽又免费观看的视频| 韩国精品一区二区三区| 午夜免费观看网址| 亚洲欧美一区二区三区黑人| 精品熟女少妇八av免费久了| 国产av一区在线观看免费| 久久精品亚洲精品国产色婷小说| 在线观看www视频免费| 最好的美女福利视频网| 久久久久久久久久黄片| 成人午夜高清在线视频 | 色老头精品视频在线观看| 国产精品影院久久| 亚洲男人天堂网一区| 免费人成视频x8x8入口观看| 国产成人av教育| xxx96com| 在线十欧美十亚洲十日本专区| 日韩免费av在线播放| 国产av一区二区精品久久| 国产av不卡久久| 亚洲五月天丁香| 首页视频小说图片口味搜索| 一级毛片女人18水好多| 90打野战视频偷拍视频| 黄色视频不卡| 欧美日韩亚洲综合一区二区三区_| 男女视频在线观看网站免费 | 国产成人一区二区三区免费视频网站| 色尼玛亚洲综合影院| 欧美乱码精品一区二区三区| 精品一区二区三区视频在线观看免费| 99re在线观看精品视频| 国产一区二区在线av高清观看| 国产精品久久视频播放| 国产97色在线日韩免费| 免费观看精品视频网站| 欧美激情久久久久久爽电影| 波多野结衣高清无吗| 男女做爰动态图高潮gif福利片| 嫩草影院精品99| 亚洲精品中文字幕一二三四区| 母亲3免费完整高清在线观看| 夜夜爽天天搞| 满18在线观看网站| 国产麻豆成人av免费视频| 午夜久久久在线观看| 波多野结衣巨乳人妻| 亚洲熟女毛片儿| 啦啦啦免费观看视频1| 成人亚洲精品一区在线观看| 久久久久九九精品影院| www国产在线视频色| 一区二区三区国产精品乱码| 亚洲一区二区三区色噜噜| 国产高清激情床上av| 90打野战视频偷拍视频| 成年人黄色毛片网站| 欧美成人免费av一区二区三区| 日韩大尺度精品在线看网址| 日韩精品中文字幕看吧| 日韩欧美国产一区二区入口| 国产成年人精品一区二区| 欧美激情高清一区二区三区| 给我免费播放毛片高清在线观看| 亚洲免费av在线视频| 观看免费一级毛片| 99国产精品一区二区三区| 夜夜看夜夜爽夜夜摸| 精品国产一区二区三区四区第35| 日韩视频一区二区在线观看| 国产精华一区二区三区| 桃色一区二区三区在线观看| 久久国产精品人妻蜜桃| 又黄又爽又免费观看的视频| 最好的美女福利视频网| 天堂√8在线中文| 亚洲国产欧美一区二区综合| 国内精品久久久久久久电影| 伦理电影免费视频| 好看av亚洲va欧美ⅴa在| 国产一区二区三区视频了| 亚洲国产看品久久| 国产三级黄色录像| www.www免费av| 狂野欧美激情性xxxx| 男男h啪啪无遮挡| 亚洲中文字幕日韩| 亚洲自拍偷在线| 99在线人妻在线中文字幕| 俺也久久电影网| 欧美一级毛片孕妇| 真人做人爱边吃奶动态| 国产精品亚洲一级av第二区| 18禁黄网站禁片免费观看直播| 久久久久久国产a免费观看| 亚洲国产欧美一区二区综合| 日本一区二区免费在线视频| 午夜免费鲁丝| 亚洲精品一卡2卡三卡4卡5卡| 久久久国产欧美日韩av| 18禁国产床啪视频网站| 搡老妇女老女人老熟妇| 国产成人av教育| 91老司机精品| 成年版毛片免费区| 久久欧美精品欧美久久欧美| 亚洲成人国产一区在线观看| 母亲3免费完整高清在线观看| 十分钟在线观看高清视频www| 一级毛片女人18水好多| 欧美+亚洲+日韩+国产| 韩国精品一区二区三区| 美女大奶头视频| 国产一区二区激情短视频| 国产色视频综合| 黄色丝袜av网址大全| 日韩精品中文字幕看吧| 午夜福利18| 久久中文字幕人妻熟女| 亚洲成a人片在线一区二区| 精品久久久久久久人妻蜜臀av| 一区二区三区高清视频在线| 亚洲精品色激情综合| 久久亚洲真实| 日韩欧美在线二视频| 岛国在线观看网站| 亚洲成人精品中文字幕电影| 极品教师在线免费播放| 亚洲欧美激情综合另类| 精品久久久久久久人妻蜜臀av| 精品国产超薄肉色丝袜足j| 久久天堂一区二区三区四区| 久久久久久久精品吃奶| 国产区一区二久久| 亚洲午夜理论影院| 国产高清视频在线播放一区| 亚洲av成人av| 桃红色精品国产亚洲av| 久久伊人香网站| 亚洲午夜精品一区,二区,三区| 亚洲成人精品中文字幕电影| 亚洲性夜色夜夜综合| 免费在线观看视频国产中文字幕亚洲| 91九色精品人成在线观看| 人人妻,人人澡人人爽秒播| 久久久精品欧美日韩精品| 18禁黄网站禁片午夜丰满| 国产在线精品亚洲第一网站| 侵犯人妻中文字幕一二三四区| 久久精品夜夜夜夜夜久久蜜豆 | 巨乳人妻的诱惑在线观看| 天天添夜夜摸| 国产真实乱freesex| 天堂影院成人在线观看| 成人国产综合亚洲| 中文字幕人成人乱码亚洲影| 国产一区二区在线av高清观看| 久久久久九九精品影院| 一本精品99久久精品77| 熟女电影av网| 免费av毛片视频| 老鸭窝网址在线观看| 国产v大片淫在线免费观看| 成人永久免费在线观看视频| av中文乱码字幕在线| 91成年电影在线观看| 久久这里只有精品19| 精品高清国产在线一区| 国产成人欧美| 免费电影在线观看免费观看| 国产精品久久视频播放| 老司机午夜福利在线观看视频| 免费在线观看视频国产中文字幕亚洲| 久久人人精品亚洲av| 99国产综合亚洲精品| 久久久久久亚洲精品国产蜜桃av| 在线国产一区二区在线| 此物有八面人人有两片| 欧美日韩中文字幕国产精品一区二区三区| 麻豆成人午夜福利视频| 怎么达到女性高潮| 欧美三级亚洲精品| 中文字幕人妻丝袜一区二区| 免费人成视频x8x8入口观看| 可以免费在线观看a视频的电影网站| 成人永久免费在线观看视频| 2021天堂中文幕一二区在线观 | 国产人伦9x9x在线观看| 91九色精品人成在线观看| 人人妻人人澡欧美一区二区| 免费高清在线观看日韩| bbb黄色大片| 国产成人精品久久二区二区免费| 成人av一区二区三区在线看| 久久精品影院6| 欧美色视频一区免费| 日日爽夜夜爽网站| 国产高清视频在线播放一区| 99精品欧美一区二区三区四区| 成人手机av| 黄色a级毛片大全视频| 麻豆久久精品国产亚洲av| 亚洲熟妇中文字幕五十中出| 亚洲人成电影免费在线| 成人欧美大片| 免费看a级黄色片| 51午夜福利影视在线观看| 欧美一级毛片孕妇| 国产精品一区二区精品视频观看| 在线观看一区二区三区| 麻豆av在线久日| 国产av一区在线观看免费| 亚洲专区国产一区二区| 18禁国产床啪视频网站| tocl精华| 免费观看精品视频网站| 亚洲 国产 在线| videosex国产| 国产成人欧美在线观看| 国产亚洲精品久久久久久毛片| 他把我摸到了高潮在线观看| 欧美人与性动交α欧美精品济南到| 日韩大尺度精品在线看网址| 天天躁夜夜躁狠狠躁躁| 日韩视频一区二区在线观看| 欧美丝袜亚洲另类 | 在线观看www视频免费| 亚洲国产欧洲综合997久久, | 久久精品91蜜桃| 十八禁网站免费在线| 午夜老司机福利片| 真人做人爱边吃奶动态| 精品人妻1区二区| 两个人视频免费观看高清| 正在播放国产对白刺激| 熟女电影av网| 日韩欧美三级三区| 午夜精品久久久久久毛片777| 免费一级毛片在线播放高清视频| 亚洲 国产 在线| 国内毛片毛片毛片毛片毛片| 国产精品亚洲av一区麻豆| www日本在线高清视频| 欧美av亚洲av综合av国产av| 搡老妇女老女人老熟妇| 成人18禁在线播放| 成人特级黄色片久久久久久久| 男女床上黄色一级片免费看| 真人一进一出gif抽搐免费| 9191精品国产免费久久| 免费看a级黄色片| 两个人视频免费观看高清| 国内久久婷婷六月综合欲色啪| 久久香蕉精品热| 琪琪午夜伦伦电影理论片6080| 国产区一区二久久| 国语自产精品视频在线第100页| 精品乱码久久久久久99久播| 国产熟女午夜一区二区三区| 亚洲av成人一区二区三| 亚洲国产精品sss在线观看| 日韩av在线大香蕉| 色尼玛亚洲综合影院| 欧美三级亚洲精品| 国产一区二区激情短视频| 欧美黄色淫秽网站| 国产精品1区2区在线观看.| 国产精品一区二区三区四区久久 | 中文字幕av电影在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人av激情在线播放| 欧美日韩精品网址| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲成av片中文字幕在线观看| 一进一出抽搐动态| 美国免费a级毛片| 草草在线视频免费看| 制服诱惑二区| 91在线观看av| 搡老岳熟女国产| 91国产中文字幕| 精品国产亚洲在线| 国产又色又爽无遮挡免费看| 精品久久久久久久末码| 高潮久久久久久久久久久不卡| 少妇粗大呻吟视频| 国产蜜桃级精品一区二区三区| www.www免费av| 91老司机精品| 国产精品99久久99久久久不卡| 国产午夜精品久久久久久| 国产精品亚洲av一区麻豆| 变态另类成人亚洲欧美熟女| 波多野结衣高清作品| 亚洲人成伊人成综合网2020| 国产伦在线观看视频一区| 国产欧美日韩一区二区三| 精品熟女少妇八av免费久了| 久久精品国产综合久久久| 欧美av亚洲av综合av国产av| a级毛片在线看网站| 欧美最黄视频在线播放免费| 无遮挡黄片免费观看| 中亚洲国语对白在线视频| 国产精品野战在线观看| 丁香欧美五月| 高清毛片免费观看视频网站| 国产精品野战在线观看| 欧美在线黄色| tocl精华| 99国产精品一区二区三区| 黄片小视频在线播放| 老司机午夜十八禁免费视频| 777久久人妻少妇嫩草av网站| 中文字幕人妻丝袜一区二区| 免费在线观看成人毛片| 天天一区二区日本电影三级| 午夜免费观看网址| 国产成人av教育| 精品国内亚洲2022精品成人| 国产精品久久久av美女十八| 大香蕉久久成人网| av电影中文网址| 黄片小视频在线播放| 亚洲欧美激情综合另类| 国产精品一区二区精品视频观看| 啦啦啦 在线观看视频| 国产单亲对白刺激| 熟女电影av网| 俄罗斯特黄特色一大片| 级片在线观看| 日本五十路高清| 午夜免费激情av| 9191精品国产免费久久| 丰满的人妻完整版| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲va日本ⅴa欧美va伊人久久| 在线免费观看的www视频| 两人在一起打扑克的视频| 日韩中文字幕欧美一区二区| 亚洲性夜色夜夜综合| 夜夜躁狠狠躁天天躁| 99re在线观看精品视频| 久9热在线精品视频| 国产黄片美女视频| 色综合站精品国产| www.熟女人妻精品国产| 久久欧美精品欧美久久欧美| 国产成人影院久久av| 日本熟妇午夜| 黄色视频不卡| 99热这里只有精品一区 | 久久精品国产亚洲av高清一级| 亚洲在线自拍视频| 国产又色又爽无遮挡免费看| 日韩有码中文字幕| 麻豆成人av在线观看| 俄罗斯特黄特色一大片| 久久久久久久久免费视频了| 国产成人系列免费观看| 日韩有码中文字幕| 欧美成人午夜精品| 搡老妇女老女人老熟妇| 日日摸夜夜添夜夜添小说| 在线观看一区二区三区| 可以在线观看的亚洲视频| 在线观看66精品国产| 18禁美女被吸乳视频| 亚洲午夜理论影院| 久久精品影院6| 黄色丝袜av网址大全| 精品久久久久久,| 国产真实乱freesex| 女同久久另类99精品国产91| 日本a在线网址| 成人精品一区二区免费| 色综合亚洲欧美另类图片| 首页视频小说图片口味搜索| 99久久久亚洲精品蜜臀av| 国产精品久久久久久人妻精品电影| 丁香六月欧美| 1024香蕉在线观看| 中文字幕人妻丝袜一区二区| 18禁黄网站禁片午夜丰满| 欧美成人午夜精品| 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 久久久久久久久久黄片| 精品午夜福利视频在线观看一区| 男人舔奶头视频| 最近在线观看免费完整版| 一级a爱视频在线免费观看| 欧美激情高清一区二区三区| 99精品在免费线老司机午夜| 亚洲激情在线av| 亚洲男人的天堂狠狠| 2021天堂中文幕一二区在线观 | 欧美三级亚洲精品| 老鸭窝网址在线观看| 成人永久免费在线观看视频| 国产精品九九99| 久久国产精品人妻蜜桃| 一区福利在线观看| 国产区一区二久久| 久久欧美精品欧美久久欧美| 1024手机看黄色片| 亚洲专区国产一区二区| 最近最新中文字幕大全免费视频| 久久久久久久久中文| 人人妻人人看人人澡| 亚洲男人天堂网一区| 成人三级做爰电影| 国产真人三级小视频在线观看| 久久久久久亚洲精品国产蜜桃av| 动漫黄色视频在线观看| 男女视频在线观看网站免费 | 国产精品乱码一区二三区的特点| 99国产综合亚洲精品| 日日摸夜夜添夜夜添小说| 国产精品,欧美在线| 亚洲自偷自拍图片 自拍| 黄色视频不卡| 免费观看精品视频网站| 日本精品一区二区三区蜜桃| 欧美国产日韩亚洲一区| 村上凉子中文字幕在线| 久久久精品国产亚洲av高清涩受| 色哟哟哟哟哟哟| 成人手机av| 欧美成人午夜精品| 最近最新免费中文字幕在线| 亚洲色图av天堂| 夜夜爽天天搞| 国产乱人伦免费视频| 国产高清有码在线观看视频 | 久久久久久久久免费视频了| 日韩欧美一区二区三区在线观看| 成人18禁在线播放| 亚洲精品久久成人aⅴ小说| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲精品第一综合不卡| 老司机午夜十八禁免费视频| 日韩精品青青久久久久久| 一进一出抽搐动态| 好男人电影高清在线观看| 最近在线观看免费完整版| 久久精品91蜜桃| 91国产中文字幕| 亚洲熟女毛片儿| 大香蕉久久成人网| 亚洲中文字幕日韩| 一边摸一边抽搐一进一小说| 亚洲精品一卡2卡三卡4卡5卡| 日本免费a在线| 国产免费男女视频| or卡值多少钱| 看黄色毛片网站| 女性生殖器流出的白浆| 热re99久久国产66热| 精品久久久久久久毛片微露脸| 日日爽夜夜爽网站| 长腿黑丝高跟| 脱女人内裤的视频| 黄色视频不卡| 免费一级毛片在线播放高清视频| 国产精品日韩av在线免费观看| 一区二区三区国产精品乱码| 欧美激情 高清一区二区三区| 精品第一国产精品| 国产成年人精品一区二区| 精品久久久久久久久久免费视频| 亚洲成av片中文字幕在线观看| 欧美一级a爱片免费观看看 | 琪琪午夜伦伦电影理论片6080| 真人做人爱边吃奶动态| 91成年电影在线观看| 久久精品国产综合久久久| 欧美性猛交╳xxx乱大交人| 9191精品国产免费久久| 亚洲性夜色夜夜综合| 91九色精品人成在线观看| 国内毛片毛片毛片毛片毛片| 中文字幕精品免费在线观看视频| 黄色视频,在线免费观看| 韩国av一区二区三区四区| 99国产综合亚洲精品| 又黄又粗又硬又大视频| 久久久久久久精品吃奶| 精品福利观看| 亚洲五月天丁香| 黄网站色视频无遮挡免费观看| 国产欧美日韩一区二区精品| ponron亚洲| 国产精品美女特级片免费视频播放器 | 精华霜和精华液先用哪个| 成人三级做爰电影| 两个人免费观看高清视频| 国产成人啪精品午夜网站| 久久天堂一区二区三区四区| 男女那种视频在线观看| 十八禁网站免费在线| 免费一级毛片在线播放高清视频| 黑丝袜美女国产一区| 黑人操中国人逼视频| 一夜夜www| 黑人巨大精品欧美一区二区mp4| 久久九九热精品免费| 国产99白浆流出| 精品熟女少妇八av免费久了| 人人妻,人人澡人人爽秒播| 国产成人欧美在线观看| 国产高清激情床上av| xxx96com| 国产亚洲欧美在线一区二区| 亚洲午夜理论影院| 国产成人精品无人区| 国产视频一区二区在线看| 久久午夜亚洲精品久久| 国产亚洲精品久久久久久毛片| 久久性视频一级片| 日本免费一区二区三区高清不卡| 久久这里只有精品19| 一级a爱片免费观看的视频| 欧美国产日韩亚洲一区| 国产精品二区激情视频| 午夜日韩欧美国产| 日日干狠狠操夜夜爽| 久久香蕉精品热| 日本一区二区免费在线视频| 人人妻人人澡欧美一区二区| 特大巨黑吊av在线直播 | 在线av久久热| avwww免费| 成人18禁高潮啪啪吃奶动态图| 长腿黑丝高跟| 搡老熟女国产l中国老女人| 99国产精品一区二区蜜桃av| 午夜两性在线视频| av福利片在线| 国产成人系列免费观看| 天堂影院成人在线观看| 在线天堂中文资源库| 亚洲国产日韩欧美精品在线观看 | 性欧美人与动物交配| 久久99热这里只有精品18| 国内精品久久久久精免费| 国产精品爽爽va在线观看网站 | 男女下面进入的视频免费午夜 | 天天躁夜夜躁狠狠躁躁| 狂野欧美激情性xxxx| 国产欧美日韩一区二区三| 国产精品野战在线观看| 一级a爱片免费观看的视频| 日韩欧美国产在线观看|