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

    RegCM4.0不同對流參數(shù)化方案對山東省氣溫和降水的模擬

    2016-02-15 02:17:53孟祥新雒佳麗王麗娟
    干旱氣象 2016年6期
    關(guān)鍵詞:對流降水量降水

    商 林,孟祥新,雒佳麗,王麗娟

    (1.山東省氣候中心,山東 濟南 250031;2.蘭州大學大氣科學學院,半干旱氣候變化教育部重點實驗室,甘肅 蘭州 730000;3.中國氣象局蘭州干旱氣象研究所,甘肅 蘭州 730020)

    RegCM4.0不同對流參數(shù)化方案對山東省氣溫和降水的模擬

    商 林1,孟祥新1,雒佳麗2,王麗娟3

    (1.山東省氣候中心,山東 濟南 250031;2.蘭州大學大氣科學學院,半干旱氣候變化教育部重點實驗室,甘肅 蘭州 730000;3.中國氣象局蘭州干旱氣象研究所,甘肅 蘭州 730020)

    利用NCAR/NNRP2的每日4次再分析資料、NOAA的周平均OI_WK再分析資料和山東省122個觀測站逐日氣溫和降水資料,使用RegCM4.0區(qū)域氣候模式,選取MIT-Emanuel、Grell和Kuo3種對流參數(shù)化方案,對山東省1990—2009年氣溫和降水進行了數(shù)值模擬。結(jié)果表明:3種對流參數(shù)化方案均能模擬出山東省平均氣溫和降水量的年際變化,且對氣溫的模擬效果總體好于降水。模擬的年平均氣溫偏低,其中Grell試驗模擬結(jié)果最低,而Ema和Kuo試驗模擬的年平均氣溫在內(nèi)陸地區(qū)存在冷偏差(約-1℃),在沿海地區(qū)存在暖偏差(約+1.5℃)。模擬的年降水量偏少,其中Kuo試驗模擬值最少。該模式對降水的模擬效果冬季最好,夏季最差。其中,夏季降水量的模擬值Ema試驗較觀測值偏多13.2%,Grell試驗較觀測值偏少17.8%。Grell和Kuo試驗,對夏季和秋季降水量的模擬值與觀測值的相關(guān)系數(shù)最高。

    RegCM4.0;對流參數(shù)化方案;氣溫;降水

    商 林,孟祥新,雒佳麗,等.RegCM4.0不同對流參數(shù)化方案對山東省氣溫和降水的模擬[J].干旱氣象,2016,34(6):995-1002,[SHANG Lin,MENG Xiangxin,LUO Jiali,et al.Simulation of Temperature and Precipitation Based on RegCM4.0 with Different Convection Parameterization Schemes in Shandong Province During 1990-2009[J].Journal of Arid Meteorology,2016,34(6):995-1002],DOI:10.11755/j.issn.1006-7639 (2016)-06-0995

    引 言

    區(qū)域氣候模式可用于區(qū)域尺度上的氣候模擬和預測,是研究區(qū)域氣候變化重要且有效的手段之一。目前基于區(qū)域氣候模式發(fā)展和模擬應(yīng)用的研究已有很多[1-5]。積云對流是區(qū)域氣候模式中重要的非絕熱加熱過程之一,其時空尺度是一般模式無法分辨的次網(wǎng)格尺度。同時,由于沒有觀測資料,且受計算機技術(shù)方面的局限性,在數(shù)值模擬和預報中一般采用積云對流參數(shù)化方案。但由于未知量個數(shù)多于方程個數(shù),出現(xiàn)了閉合假設(shè)問題,因此選擇合適的閉合假設(shè)成為積云參數(shù)化的一個關(guān)鍵。

    國內(nèi)外對區(qū)域氣候模式RegCM的積云對流參數(shù)化方案進行了大量研究。Giorgi等[6]指出,區(qū)域氣候模式RegCM的參數(shù)化過程是針對特定的氣候狀況和分辨率設(shè)定的,對于相同的參數(shù)設(shè)置,模擬結(jié)果會因模擬區(qū)域不同而有很大差別,因此,針對模擬區(qū)域選取合適的參數(shù)化方案要進行大量反復的調(diào)試。眾多研究顯示,RegCM模式的模擬性能對積云對流參數(shù)化方案有較大的敏感性[7-11],不同的對流參數(shù)化方案對中國東部季風區(qū)和江淮流域降水在不同年份表現(xiàn)效果并不一致[12-17]。這可能是由于選擇個例不同造成對流參數(shù)化方案的表現(xiàn)效果不同。但同時也說明,對于哪種對流參數(shù)化方案更適合中國東部季風區(qū)還沒有一致的結(jié)論,仍然是一個值得研究的問題。中國東部季風區(qū)夏季降水集中,選擇合適的對流參數(shù)化方案尤為重要。

    以往模擬研究的陸地區(qū)域多為全國或者流域尺度,且水平分辨率較低。高學杰等[18]使用RegCM2對東亞地區(qū)降水的模擬結(jié)果顯示,東亞地區(qū)降水的模擬效果取決于模式的水平分辨率,水平分辨率越高,模擬效果越好。山東地處南北氣候過渡帶,旱澇既不與華北地區(qū)同步,也和江淮地區(qū)不一致[19],因此,有必要使用RegCM對山東省進行高分辨率的區(qū)域氣候模擬和應(yīng)用研究。胡軼佳等[20]的研究結(jié)果顯示,RegCM3.0中 Kuo和Grell對流參數(shù)化方案對1998年中國東部降水場形勢模擬效果的好壞隨季節(jié)交替出現(xiàn)。山東省屬季風氣候,降水集中,冬季以降雪為主,春秋季以非對流性降水為主,夏季對流性降水與非對流性降水的比重都很大,因此,不同對流參數(shù)化方案對山東省不同季節(jié)的氣候模擬效果不一致。針對以上問題,本文使用RegCM4.0區(qū)域氣候模式,模擬了山東省1990—2009年氣溫和降水,檢驗采用不同對流參數(shù)化方案的RegCM4.0模式對山東省的模擬能力,揭示不同對流參數(shù)化方案在不同季節(jié)的模擬效果,為該模式對山東省進行氣候預測和診斷提供基礎(chǔ)。

    1 資料與方法

    1.1 資 料

    所用資料包括:(1)模式的初始場和邊界場使用1989年11月—2009年12月每日4次的NCAR/NNRP2再分析資料,水平分辨率為2.5°×2.5°;(2)海表溫度使用同期美國國家海洋和大氣管理局(NOAA)的周平均OI_WK再分析資料,水平分辨率為1°×1°;(3)山東省氣象信息中心提供的1990年1月—2009年12月全省122個觀測站(圖1)逐日氣溫和降水資料,用于檢驗模式模擬結(jié)果。為便于分析,以12月至翌年2月為冬季,3—5月為春季,6—8月為夏季,9—11月為秋季生成序列。

    圖1 山東省122個觀測站分布Fig.1 Spatial distribution of122 meteorological observation stations in Shandong Province

    1.2 插值方法

    采用Cressman方法對山東省 122個觀測站資料進行插值,該方法是將離散點內(nèi)插到規(guī)則格點上而引起誤差較小的一種逐步訂正內(nèi)插方法[21],被廣泛應(yīng)用于氣象領(lǐng)域診斷分析和數(shù)值預報方案的客觀分析 中[22]。

    模式模擬結(jié)果與觀測數(shù)據(jù)進行比較和相關(guān)性分析時,對模擬數(shù)據(jù)采用雙線性內(nèi)插法(Bilinear Interpolation),將格點資料插值到站點上。該方法取(x,y)點周圍的4鄰點,在y方向內(nèi)插2次,再在x方向內(nèi)插1次,得到(x,y)點的值f(x,y)。

    2 區(qū)域氣候模式及試驗方案

    使用的模式為2010年6月發(fā)行的RegCM4.0,是RegCM相對較新的版本,其動力內(nèi)核基于中尺度模式MM5,模式垂直方向上使用地形坐標σ。與之前的版本相比,該版本更新包括了新的陸面過程、行星邊界層、海—氣通量方案和混合對流方案,修改了原有的輻射傳輸和邊界層方案,并全面升級了模式代碼,以提高靈活性和可移植性。模式的最高水平分辨率為10 km×10 km,山東省南北寬約420 km,東西長約800 km,RegCM4.0模式可以完全滿足對山東省的模擬需求。

    模擬的區(qū)域中心點設(shè)為(120°E,37°N),水平分辨率為10 km,南北方向格點數(shù)為120,東西方向格點數(shù)為180。模式垂直分層為18層,頂層氣壓為5.0 hPa。大氣模塊的積分步長為30 s,陸地模塊的積分步長為90 s。方案選擇依照以往研究結(jié)果[23-25],側(cè)邊界方案為指數(shù)松弛方案,側(cè)邊界每6 h輸入一次以驅(qū)動區(qū)域氣候模式。投影選擇蘭伯特投影,海洋通量參數(shù)化方案采用Zeng方案。需要指出的是,RegCM4.0的BATS陸面過程方案中新增了2種土地利用類型來描述城市和郊區(qū)的環(huán)境,而山東的城鎮(zhèn)化率已經(jīng)超過50%[26],因此本文選取BATS陸面過程方案來描述山東省的陸面過程。

    本文共設(shè)計了3個試驗:Ema試驗,選擇Massachusetts Institute of Technology(MIT)-Emanuel積云對流參數(shù)化方案,該方案的云模式是隨機混合模式,主要假設(shè)云尺度的上升氣流作為對流層輸送的主要因子[27];Grell試驗,選擇Grell積云對流參數(shù)化方案,該方案依賴于不穩(wěn)定化速率的閉合假設(shè),認為云內(nèi)有上升和下沉引起的2個穩(wěn)定氣流,除環(huán)流頂部和底部外,在云內(nèi)與環(huán)境空氣之間無直接的混合[28];Kuo試驗,選擇Kuo積云對流參數(shù)化方案,該方案采用質(zhì)量、水汽輻合型假設(shè),即水汽輻合總量的一部分用于凝結(jié)和降水,另一部分則用于濕潤環(huán)境空氣[29]。模式初始時間一般選擇天氣系統(tǒng)較為穩(wěn)定的秋冬季節(jié)[30]。積分時間為1989年11月1日—2009年12月31日;初始化階段為1989年11 月1日—12月31日,1990年1月1日—2009年12 月31日用于結(jié)果分析。山東省地勢較為平緩,海拔較高的地區(qū)位于山東半島內(nèi)陸以及中部的泰山山脈(圖2)。

    圖2 山東省地形高度分布(單位:m)Fig.2 Altitude of terrain in Shandong Province(Unit:m)

    3 山東省氣候概況

    泰山位于山東中部,東西長約200 km,南北寬約20 km,泰山站位于主峰玉皇頂東南方向的日觀峰上,海拔1 534 m,是我國建站最早的高山氣象站之一。圖3給出1990—2009年山東省站點觀測的年平均氣溫和年降水量。其中,圖3a和圖3b為包括泰山站在內(nèi)的122個站點數(shù)據(jù),可看出,在泰山站及其周邊存在年平均氣溫的異常低值中心和年降水量的高值區(qū)。這是因為當臺站資料分布極為不均時,倘若在一大片區(qū)域內(nèi)只有一個較小集中區(qū)域溫度異常低或者降水異常多時,采用Cressman方法插值后該區(qū)域周圍會出現(xiàn)虛假的圓形插值區(qū)域,且范圍過大,插值的邊界也變得不合理[31]。因此在分析山東省氣候特征時,一般不考慮泰山站(圖3c、圖3d)。圖3c顯示,山東省年平均氣溫受緯度、地形和海陸分布的影響,表現(xiàn)出由東北向西南逐漸升高的空間分布特征,年平均氣溫為12~15℃。其中,在內(nèi)陸地區(qū),低緯度地區(qū)年平均氣溫高于高緯度地區(qū),南北溫差達1.5℃;山東半島內(nèi)陸地區(qū)海拔較高,年平均氣溫低于同緯度其他地區(qū),而半島東部沿海地區(qū)受海洋影響較大,年平均氣溫較低,低于12.5℃。山東省降水分布表現(xiàn)為南多北少,平均年降水量在500~1 000 mm之間,半島南部地區(qū)受山區(qū)地形和海洋暖濕氣流的影響,降水較多(圖3d)。綜上所述,山東省氣溫和降水分布受緯度、地形和海洋的共同影響。

    圖3 1990—2009年山東122(a、b)和121(c、d)個觀測站的年平均氣溫(a、c,單位:℃)和年降水量(b、d,單位:mm)的空間插值(圖a、b中圓圈為泰山站所在位置)Fig.3 The spatial interpolation of the observed annualmean temperature(a,c,Unit:℃)and annual precipitation (b,d,Unit:mm)of 122(a,b)and 121(c,d)stations in Shandong Province during 1990-2009 (The circles in fig.a(chǎn) and fig.b represent the location of Taishan station)

    4 RegCM4.0對山東省氣溫和降水的模擬

    圖4給出1990—2009年RegCM4.0模擬的年平均氣溫和年降水量。與觀測相比(圖3c),RegCM4.0可模擬出山東內(nèi)陸低緯度地區(qū)年平均氣溫高于高緯度地區(qū),半島內(nèi)陸地區(qū)氣溫較低(圖4a、圖4c、圖4e),但3個試驗的模擬結(jié)果在泰山地區(qū)均存在一冷中心,且整體存在一定的冷偏差。這與已有的研究結(jié)果“RegCM3/4模式對東亞地區(qū)氣溫的模擬存在系統(tǒng)性的冷偏差”[32-34]一致。其中,Ema 和Kuo試驗在山東內(nèi)陸地區(qū)年平均氣溫偏低約1℃,半島內(nèi)陸地區(qū)相差不大,而Grell試驗結(jié)果偏冷更為明顯。同時,模式對近海岸地區(qū)的模擬能力較差,山東北部近海岸地區(qū)氣溫明顯偏低,而半島東部近海岸地區(qū)氣溫偏高約1.5℃。

    與年降水量觀測值(圖3d)相比,RegCM4.0模式基本上可模擬出山東省降水量南多北少的特征,但降水量值明顯偏少。其中,Grell試驗模擬的降水量整體約為觀測值的一半;Kuo試驗對山東西部的降水特征模擬能力不足,年降水量<200 mm,遠低于觀測的600 mm,但對其他地區(qū)的模擬結(jié)果與觀測值相差不大;Ema試驗的模擬結(jié)果較好,模擬的年降水量總體上比觀測值小200 mm。值得注意的是,Ema試驗的降水模擬結(jié)果在泰山山脈存在一降水偏多中心,而在山東南部偏少400 mm。

    圖4 Ema試驗(a、b)、Grell試驗(c、d)和 Kuo試驗(e、f)模擬的1990—2009年山東省年平均氣溫(a、c、e,單位:℃)和年降水量(b、d、f,單位:mm)Fig.4 The simulated annualmean temperature(a,c,e,Unit:℃)and precipitation(b,d,f,Unit:mm)of Ema(a,b),Grell(c,d)and Kuo(e,f)tests in Shandong Province during 1990-2009

    圖5為1990—2009年山東省觀測和模擬的年平均氣溫和年降水量變化。可以看出,與觀測時間序列相比,Ema試驗、Grell試驗和Kuo試驗均能夠模擬出年平均氣溫和年降水量的變化趨勢,但模擬值與觀測值存在一定的偏差。對年平均氣溫而言(圖5a),3個試驗的模擬值較觀測值偏低,其中Kuo試驗偏低0.7℃,Ema試驗偏低0.8℃,Grell試驗偏低最多為1.2℃。就年降水量而言(圖5b),Grell和Kuo試驗的模擬結(jié)果明顯低于觀測值,而Ema試驗的模擬值與觀測值較接近(個別年份除外)。另外,3個試驗模擬的年平均氣溫與觀測值的相關(guān)系數(shù)均>0.97,表明該模式對氣溫的模擬效果較好,因此,下文僅著重分析 RegCM4.0對降水的模擬情況。

    圖5 觀測和模擬的1990—2009年山東省年平均氣溫(a)和年降水量(b)的變化Fig.5 The annual variations of the simulated and observated average temperature(a)and precipitation(b)in Shandong Province during 1990-2009

    山東屬于暖溫帶季風氣候,降水年內(nèi)分布很不均勻,全年降水量60%~70%集中在夏季,春秋季降水較少,冬季降水最少。表1是觀測和模擬的山東省降水量的季節(jié)分布,可看出,各試驗均能模擬出降水量的季節(jié)分布特征。其中,Ema試驗模擬的夏季降水量較觀測值偏多13.2%,其余季節(jié)均偏少,且秋季偏少最多;Kuo試驗模擬的各季節(jié)降水量均較觀測值偏少;Grell試驗模擬的夏季降水量偏少17.8%,其他季節(jié)降水量偏少更為明顯。3個試驗對比發(fā)現(xiàn),Grell試驗模擬的季節(jié)降水量除夏季外,其他季節(jié)與觀測值相差甚遠。圖6為觀測和模擬的山東省降水量的逐月分布,可看出各試驗均能模擬出降水量的月分布,但不同試驗的模擬效果不盡相同??傮w上,觀測值大于模擬值,但Ema試驗模擬的夏季降水量高于觀測值,其中,6月和7月降水量分別偏多10 mm和40 mm。這主要是由于MITEmanuel對流參數(shù)化方案夏季模擬的西南氣流偏北分量較大,水汽更容易由西南氣流輸送到山東地區(qū),導致該地區(qū)夏季降水模擬偏多;Kuo和Grell對流參數(shù)化方案夏季在孟加拉灣存在較強的氣旋性環(huán)流,且西太平洋副熱帶高壓偏強,導致北方地區(qū)夏季降水偏弱[15]。春季,Ema試驗的模擬結(jié)果與觀測值(107 mm)最為接近,偏少約15 mm,其中在降水量相對較多的4月和5月模擬效果較好。秋季,Kuo試驗的模擬結(jié)果與觀測值(117 mm)最為接近,偏少約21 mm,且在秋季各月的模擬值都比較穩(wěn)定。

    表1 1990—2009年山東省各季節(jié)模擬的降水量與觀測值對比Tab.1 The comparison between the observed and simulated precipitation in different seasons in Shandong Province during 1990-2009

    圖6 觀測和3個試驗模擬的1990—2009年山東省月降水量Fig.6 The observed and simulated precipitation in differentmonths in Shandong Province during 1990-2009

    表2是3個試驗模擬的降水量與觀測值的相關(guān)系數(shù),且相關(guān)系數(shù)均通過0.01顯著性檢驗??芍?個試驗對夏季降水模擬的偏差較大,夏季降水量的模擬值與觀測值的相關(guān)系數(shù)最低,其中Grell試驗模擬的夏季降水量與觀測值的相關(guān)系數(shù)最高為0.517,而Ema試驗和Kuo試驗模擬值與觀測值的相關(guān)系數(shù)均<0.5;模式模擬的冬季降水量與觀測值的相關(guān)系數(shù)最高,Ema試驗的模擬值與觀測值的相關(guān)系數(shù)最高達0.779;春季和秋季,Kuo試驗的模擬值與觀測值的相關(guān)系數(shù)最高,分別為0.595和0.702。就全年而言,Grell試驗對年降水模擬的結(jié)果與觀測值的相關(guān)系數(shù)最高,達0.709。

    表2 1990—2009年山東省降水量的模擬值與觀測值的相關(guān)系數(shù)Tab.2 The correlation coefficients between the simulated and observed precipitation in Shandong Province during 1990-2009

    表1顯示,夏季和秋季是山東降水量最多和次多的季節(jié),因此對夏秋季降水量的模擬值與觀測值的相關(guān)性進行空間分析(圖7)。由圖7可見,降水模擬值與觀測值的相關(guān)系數(shù)秋季高于夏季,說明該模式對夏季降水的模擬能力較弱。夏季,Grell試驗的相關(guān)系數(shù)普遍高于Ema和Kuo試驗,相關(guān)系數(shù)在0.45~0.6之間,其中山東東南部沿海、西南部和北部的部分內(nèi)陸地區(qū)相關(guān)系數(shù)較低。相對于降水觀測值,Ema試驗夏季的偏差小于Grell試驗(表1),但Grell試驗夏季的相關(guān)系數(shù)較高。秋季,Kuo試驗的相關(guān)系數(shù)高于Ema和Grell試驗,大部分地區(qū)的相關(guān)系數(shù)在0.75以上,但對山東南部地區(qū)的模擬效果相對較差,部分地區(qū)相關(guān)系數(shù)在0.6以下。

    圖7 1990—2009年Ema(a、b)、Grell(c、d)和Kuo(e、f)試驗模擬的山東省夏季(a、c、e)和秋季(b、d、f)降水量與觀測值的相關(guān)系數(shù)Fig.7 The correlation coefficients between the observed and simulated precipitation by Ema(a,b),Grell(c,d)and Kuo(e,f)tests in summer(a,c,e)and autumn(b,d,f)in Shandong Province during 1990-2009

    5 結(jié) 論

    (1)RegCM4.0的3種對流參數(shù)化方案均能夠模擬出山東省年平均氣溫和年降水量隨時間的變化,總體上對氣溫的模擬效果好于降水。模式模擬的年平均氣溫較觀測偏低,其中Grell試驗的模擬值最低;Ema和Kuo試驗模擬的年平均氣溫在山東內(nèi)陸地區(qū)存在冷偏差(約-1℃),在沿海地區(qū)存在暖偏差(約+1.5℃)。模式模擬的年降水量偏少,其中Kuo試驗的模擬值最少。

    (2)RegCM4.0模式對山東冬季降水的模擬效果最好,夏季最差,但各對流參數(shù)化方案的模擬效果在不同季節(jié)表現(xiàn)并不一致。夏季,Ema試驗的模擬值較觀測值偏多13.2%,Grell試驗的模擬值較觀測值偏少17.8%。雖然Grell試驗模擬的夏季降水偏差更大,但其與觀測值的相關(guān)系數(shù)最高;春季和秋季,Kuo試驗的模擬值與觀測值的相關(guān)系數(shù)最高;冬季,Ema試驗的模擬值與觀測值的相關(guān)系數(shù)最高??梢?,選擇最優(yōu)的對流參數(shù)化方案,以提高山東省不同季節(jié)的氣候模擬效果。

    [1]趙宗慈,羅勇.二十世紀九十年代區(qū)域氣候模擬研究進展[J].氣象學報,1998,56(2):225-246.

    [2]劉鴻波,張大林,王斌.區(qū)域氣候模擬研究及其應(yīng)用進展[J].氣候與環(huán)境研究,2006,11(5):649-668.

    [3]石英.中國區(qū)域氣候變化的高分辨率數(shù)值模擬[D].北京:中國氣象科學研究院,2007.

    [4]Ding Y H,Shi X L,Liu Y M,et al.Multi-year simulations and experimental seasonal predictions for rainy seasons in China by using a nested regional climatemodel(RegCM_NCC).Part I:Sensitivity study[J].Advances in Atmospheric Sciences,2006,23(3):323-341.

    [5]Ding Y H,Liu Y M,Shi X L,et al.Multi-Year Simulations and Experimental Seasonal Predictions for Rainy Seasons in China by U-sing a Nested Regional Climate Model(RegCM_NCC)PartⅡ:The experimental seasonal prediction[J].Advances in Atmospheric Sciences,2006,23(4):487-503.

    [6]Giorgi F,Marinucci M R.A investigation of the sensitivity of simulated precipitation tomodel resolution and its implications for climate studies[J].Monthly Weather Review,1996,124(1):148-166.

    [7]劉一鳴,丁一匯.修正的質(zhì)量通量積云對流方案及其模擬試驗研究I:方案介紹及對1991年洪澇過程的模擬[J].氣象學報,2001,59(1):10-22.

    [8]Dash SK,Shekhar M S,Singh G P.Simulation of Indian summer monsoon circulation and rainfall using RegCM3[J].Theoretical and Applied Climatology,2006,86(1):161-172.

    [9]鮑艷,呂世華,左洪超,等.RegCM3模式在西北地區(qū)的應(yīng)用研究Ⅱ:區(qū)域選擇及參數(shù)化方案的敏感性[J].冰川凍土,2006,28 (2):175-182.

    [10]曹杰,張秀年,尤亞磊,等.區(qū)域氣候模式(RegCM3)中積云對流參數(shù)化方案在縱向嶺谷區(qū)的適用性研究[J].科學通報,2007,52(S2):101-106.

    [11]Zou L,Zhou T.Sensitivity of a regional ocean-atmosphere coupled model to convection parameterization overwestern North Pacific[J].Journal of Geophysical Research Atmospheres,2011,116 (D18):597-616.

    [12]劉永強,丁一匯,趙宗慈.1991年江淮特大異常降水的區(qū)域氣候模擬.趙宗慈.中國短期氣候預測的模式研究[M].北京:氣象出版社,1996:106-120.

    [13]潘勁松,翟國慶,高坤.區(qū)域氣候模擬中多種對流參數(shù)化方案的比較研究[J].大氣科學,2002,26(2):206-220.

    [14]劉曉東,江志紅,羅樹如,等.RegCM3模式對中國東部夏季降水的模擬試驗[J].大氣科學學報,2005,28(3):351-359.

    [15]王美麗.RegCM4.4區(qū)域氣候模式對中國當代氣候及極端事件的高分辨率數(shù)值模擬[D].北京:中國氣象科學研究院,2015.

    [16]李建云,王漢杰.RegCM3積云參數(shù)化方案對中國南方夏季強降水過程模擬的影響[J].氣候與環(huán)境研究,2008,13(2):149 -160.

    [17]李清華.區(qū)域氣候模擬中對流參數(shù)化方案的比較研究[D].蘭州:蘭州大學,2006.

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

    [19]苑文華,張玉潔,孫茂璞,等.山東省降水量與不同強度降水日數(shù)變化對干旱的影響[J].干旱氣象,2010,28(1):35-40.

    [20]胡軼佳,鐘中,閔錦忠.兩種積云對流參數(shù)化方案對1998年區(qū)域氣候季節(jié)變化模擬的影響研究[J].大氣科學,2008,32(1):90-100.

    [21]Cressman G P.An operational objective analysis system[J].Monthly Weather Review,1959,87(10):367-374.

    [22]王躍山.客觀分析和四維同化——站在新世紀的回望(Ⅱ)客觀分析的主要方法(2)[J].氣象科技,2001,29(3):1-11.

    [23]朱濤.RegCM3.0與RegCM4.0對中國區(qū)域氣候模擬的對比分析[D].南京:南京信息工程大學,2012.

    [24]張井貝.區(qū)域氣候模式RegCM4對中國地區(qū)的數(shù)值模擬[D].上海:華東師范大學,2013.

    [25]鄒靖,謝正輝.RegCM4中陸面過程參數(shù)化方案對東亞區(qū)域氣候模擬的影響[J].氣象學報,2012,70(6):1312-1326.

    [26]張東紅.《大眾日報》山東城鎮(zhèn)化焦點問題系列報道[D].濟南:山東大學,2013.

    [27]Emanuel K A.A Scheme for Representing Cumulus Convection in Large-Scale Models[J].Journal of the Atmospheric Sciences,1991,48(21):2313-2329.

    [28]Grell G A.Prognostic Evaluation of Assumptions Used by Cumulus Parameterizations[J].Monthly Weather Review,1993,121(3):764-787.

    [29]Anthes R A.A Cumulus Parameterization Scheme Utilizing a One -Dimensional Cloud Model[J].Monthly Weather Review,1977,105(3):270-286.

    [30]鄭益群,高峻嶺,高藝.東亞區(qū)域氣候模擬研究初探[M].北京:氣象出版社,2010.

    [31]張紅杰,馬清云,吳煥萍,等.氣象降水分布圖制作中的插值算法研究[J].氣象,2009,35(11):131-136.

    [32]李巧萍,丁一匯.區(qū)域氣候模式對東亞季風和中國降水的多年模擬與性能檢驗[J].氣象學報,2004,62(2):140-153.

    [33]張冬峰,高學杰,趙宗慈.RegCM3區(qū)域氣候模式對中國氣候的模擬[J].氣候變化研究進展,2005,1(3):119-121.

    [34]鞏崇水,段海霞,李耀輝,等.RegCM4模式對中國過去30 a氣溫和降水的模擬[J].干旱氣象,2015,33(3):379-385.

    Simulation of Temperature and Precipitation Based on RegCM 4.0 w ith Different Convection Parameterization Schemes in Shandong Province During 1990-2009

    SHANG Lin1,MENG Xiangxin1,LUO Jiali2,WANG Lijuan3

    (1.Shandong Provincial Climate Centre,Ji'nan 250031,China;2.College of Atmospheric Sciences,Lanzhou University,Key Laboratory of Semi-Arid Climate Change of the Ministry of Education,Lanzhou 730000,China;3.Institute of Arid Meteorology,CMA,Lanzhou 730020,China)

    Based on the NCAR/NNRP2 reanalysis data,theweekly average OI_WK reanalysis data from the American NationalOceanic and Atmospheric Administration(NOAA),the daily temperature and precipitation data at 122 meteorological stations in Shandong Province,the temperature and precipitation in Shandong Province during 1990-2009 were simulated by using the regional climate model(RegCM4.0)with different convection parameterization schemes,including MIT-Emanuel,Grell and Kuo parameterization schemes.The results showed that RegCM4.0 with three convection parameterization schemes could simulate reasonably the annual changes of temperature and precipitation in Shandong Province during1990-2009,and the effectof simulation on temperature was better than precipitation.The simulated values of annualmean temperaturewere lower than the observations,and the annualmean temperatures simulated by RegCM4.0 with the Grell scheme were lowest,while thatwith the Ema and Kuo schemeswere 1℃ lower than the observation in the inland areas and were 1.5℃higher than the observation in the coastal areas of Shandong Province.The annual precipitations simulated by RegCM4.0 were less than the observations,and thatwith the Kuo schemewas least.In additional,the simulation effect of precipitation was the best in winter and theworst in summer.The simulated precipitation with the Ema scheme in summer was 13.2%more than the observation,while thatwith the Grell scheme was 17.8%less than the observation.The correlation coefficients of summer and autumn precipitation between the simulation and observation were highestwith the Grell and Kuo schemes.

    RegCM4.0;convection parameterization scheme;temperature;precipitation

    1006-7639(2016)-06-0995-08

    10.11755/j.issn.1006-7639(2016)-06-0995

    P435

    A

    2016-02-24;改回日期:2016-05-31

    山東省氣象局青年基金(2015SDQN11)資助

    商林(1987-),女,山東聊城人,博士,工程師,主要從事數(shù)值模式發(fā)展和研究工作.E-mail:shangl2005@163.com

    猜你喜歡
    對流降水量降水
    繪制和閱讀降水量柱狀圖
    齊口裂腹魚集群行為對流態(tài)的響應(yīng)
    黑龍江省玉米生長季自然降水與有效降水對比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    降水量是怎么算出來的
    啟蒙(3-7歲)(2019年8期)2019-09-10 03:09:08
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    1988—2017年呼和浩特市降水演變特征分析
    降水現(xiàn)象儀模擬軟件設(shè)計與實現(xiàn)
    基于小波變換的三江平原旬降水量主周期識別
    基于ANSYS的自然對流換熱系數(shù)計算方法研究
    二元驅(qū)油水界面Marangoni對流啟動殘余油機理
    国产精品福利在线免费观看| 丝瓜视频免费看黄片| 亚洲av成人精品一区久久| 久久久久久久久久成人| 久久热精品热| 国产精品av视频在线免费观看| 我要看日韩黄色一级片| 人妻系列 视频| 男女无遮挡免费网站观看| 成人免费观看视频高清| 26uuu在线亚洲综合色| 久久久久人妻精品一区果冻| 少妇熟女欧美另类| h视频一区二区三区| 免费av不卡在线播放| 国产 一区精品| 夫妻午夜视频| 我要看黄色一级片免费的| 内地一区二区视频在线| 国产一区二区三区综合在线观看 | 99久久综合免费| 国产无遮挡羞羞视频在线观看| 国产av一区二区精品久久 | 天堂8中文在线网| 国产美女午夜福利| 精品一区二区三卡| 22中文网久久字幕| 国产高清国产精品国产三级 | 黄片无遮挡物在线观看| 日韩成人av中文字幕在线观看| 亚洲欧美一区二区三区国产| 午夜免费观看性视频| 国产免费福利视频在线观看| 我的女老师完整版在线观看| 免费观看的影片在线观看| 直男gayav资源| 免费大片黄手机在线观看| 久久99蜜桃精品久久| 夫妻午夜视频| 亚洲欧美日韩另类电影网站 | 18禁裸乳无遮挡动漫免费视频| 欧美成人a在线观看| 欧美 日韩 精品 国产| 色5月婷婷丁香| 亚洲激情五月婷婷啪啪| 亚洲精品国产成人久久av| 老司机影院成人| 老女人水多毛片| 波野结衣二区三区在线| 一级毛片 在线播放| 99久久综合免费| 蜜桃久久精品国产亚洲av| 日韩国内少妇激情av| 日本与韩国留学比较| 精品一区在线观看国产| 国产精品一区www在线观看| 国产色婷婷99| 永久免费av网站大全| 日韩中文字幕视频在线看片 | 国产黄色免费在线视频| 成人特级av手机在线观看| 水蜜桃什么品种好| 免费看不卡的av| 日韩av不卡免费在线播放| 国产男女内射视频| 日本免费在线观看一区| 亚洲av.av天堂| 六月丁香七月| 免费看av在线观看网站| 久久久久久久久久久免费av| www.色视频.com| 国产成人a∨麻豆精品| 中文字幕亚洲精品专区| .国产精品久久| 国产精品秋霞免费鲁丝片| 黄色怎么调成土黄色| 热re99久久精品国产66热6| 尾随美女入室| 国产一区二区在线观看日韩| 能在线免费看毛片的网站| 亚洲精品一区蜜桃| 国产精品一区二区在线观看99| 成人国产麻豆网| 男男h啪啪无遮挡| 国产91av在线免费观看| 国产高清不卡午夜福利| 黑人高潮一二区| 亚洲综合精品二区| 大又大粗又爽又黄少妇毛片口| 成人毛片60女人毛片免费| 777米奇影视久久| 在线免费十八禁| 青春草视频在线免费观看| 成人毛片60女人毛片免费| 麻豆国产97在线/欧美| 一本一本综合久久| 91精品一卡2卡3卡4卡| 久久99蜜桃精品久久| 欧美日韩视频高清一区二区三区二| 国产乱来视频区| 久久av网站| 国产精品无大码| 人人妻人人添人人爽欧美一区卜 | 欧美精品人与动牲交sv欧美| 舔av片在线| 婷婷色av中文字幕| 精品国产一区二区三区久久久樱花 | 免费av中文字幕在线| 国产黄色视频一区二区在线观看| 色5月婷婷丁香| 久久国产乱子免费精品| 日韩欧美精品免费久久| 成人一区二区视频在线观看| 亚洲精品亚洲一区二区| 免费观看性生交大片5| 欧美 日韩 精品 国产| 女人久久www免费人成看片| 草草在线视频免费看| 一区二区三区免费毛片| 免费人成在线观看视频色| 成人毛片a级毛片在线播放| 国产精品三级大全| 男人狂女人下面高潮的视频| 99视频精品全部免费 在线| 日韩中文字幕视频在线看片 | 欧美日韩视频精品一区| 亚洲婷婷狠狠爱综合网| 久久这里有精品视频免费| 午夜日本视频在线| 男男h啪啪无遮挡| 中文在线观看免费www的网站| 国产免费福利视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 天美传媒精品一区二区| 直男gayav资源| 午夜福利在线观看免费完整高清在| 一边亲一边摸免费视频| 国产一级毛片在线| 久久国产精品大桥未久av | 日本一二三区视频观看| 蜜桃在线观看..| 亚洲欧美精品专区久久| 嫩草影院入口| 日韩欧美一区视频在线观看 | 偷拍熟女少妇极品色| 一级av片app| 精品酒店卫生间| 亚洲伊人久久精品综合| 日韩 亚洲 欧美在线| 97超碰精品成人国产| 欧美日韩在线观看h| 91久久精品国产一区二区成人| 亚洲图色成人| 男女边摸边吃奶| 丰满迷人的少妇在线观看| 亚洲人成网站在线播| 成人亚洲欧美一区二区av| 精品午夜福利在线看| 久久久久精品性色| 九草在线视频观看| 中国三级夫妇交换| 免费黄色在线免费观看| 99国产精品免费福利视频| 国产精品av视频在线免费观看| 高清欧美精品videossex| 精品一品国产午夜福利视频| 久久99蜜桃精品久久| 久久韩国三级中文字幕| h视频一区二区三区| 国产精品福利在线免费观看| 天天躁日日操中文字幕| 在线观看一区二区三区激情| 日本av免费视频播放| 少妇猛男粗大的猛烈进出视频| 国产一区有黄有色的免费视频| 男女啪啪激烈高潮av片| 嫩草影院入口| 中文字幕亚洲精品专区| 人妻制服诱惑在线中文字幕| 国产免费一级a男人的天堂| 国产伦精品一区二区三区四那| 内射极品少妇av片p| 街头女战士在线观看网站| 国产精品免费大片| 国产精品久久久久久av不卡| 嘟嘟电影网在线观看| 亚洲久久久国产精品| 国产高潮美女av| 国产精品免费大片| 夜夜骑夜夜射夜夜干| 久久久久久久久久久丰满| 十八禁网站网址无遮挡 | 寂寞人妻少妇视频99o| 亚洲国产av新网站| 亚洲av电影在线观看一区二区三区| 在线免费十八禁| 久久精品熟女亚洲av麻豆精品| 狠狠精品人妻久久久久久综合| 中国美白少妇内射xxxbb| 日韩av在线免费看完整版不卡| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久色成人| 成人美女网站在线观看视频| 熟女电影av网| 亚洲自偷自拍三级| av卡一久久| 国产在线一区二区三区精| 永久网站在线| 久久6这里有精品| 美女中出高潮动态图| 欧美日韩视频精品一区| 在线 av 中文字幕| 国产白丝娇喘喷水9色精品| 一边亲一边摸免费视频| 高清在线视频一区二区三区| 国产黄频视频在线观看| 18禁在线无遮挡免费观看视频| 青春草视频在线免费观看| 哪个播放器可以免费观看大片| av国产久精品久网站免费入址| 欧美变态另类bdsm刘玥| 日韩大片免费观看网站| 国产69精品久久久久777片| 99视频精品全部免费 在线| 国产精品爽爽va在线观看网站| 直男gayav资源| 男人舔奶头视频| 欧美成人a在线观看| 国产av一区二区精品久久 | 色视频www国产| 插逼视频在线观看| 91久久精品国产一区二区成人| 一级毛片aaaaaa免费看小| 深夜a级毛片| 国产精品久久久久成人av| 91精品一卡2卡3卡4卡| 我的女老师完整版在线观看| 亚洲精品国产av蜜桃| 天天躁夜夜躁狠狠久久av| 男男h啪啪无遮挡| 亚洲人成网站在线观看播放| 精品午夜福利在线看| 亚洲精品成人av观看孕妇| 青春草视频在线免费观看| 91午夜精品亚洲一区二区三区| 99热国产这里只有精品6| 少妇人妻 视频| 日韩强制内射视频| 香蕉精品网在线| 欧美激情极品国产一区二区三区 | 五月玫瑰六月丁香| 亚洲精品日本国产第一区| 亚洲三级黄色毛片| 午夜福利网站1000一区二区三区| 97在线视频观看| 国产淫片久久久久久久久| h视频一区二区三区| 在线观看av片永久免费下载| 欧美少妇被猛烈插入视频| 建设人人有责人人尽责人人享有的 | 在线观看免费视频网站a站| 亚洲av中文av极速乱| 欧美日韩视频高清一区二区三区二| 91精品国产九色| 亚洲成人av在线免费| 亚洲自偷自拍三级| 97精品久久久久久久久久精品| 人人妻人人爽人人添夜夜欢视频 | 亚洲av日韩在线播放| 一级毛片电影观看| 日日摸夜夜添夜夜爱| 亚洲精品国产av成人精品| 精品人妻一区二区三区麻豆| 国产精品久久久久久久电影| 亚洲国产欧美人成| 亚洲精品456在线播放app| 国产欧美亚洲国产| 18禁裸乳无遮挡动漫免费视频| 99热6这里只有精品| 日韩,欧美,国产一区二区三区| 国产伦精品一区二区三区视频9| 亚洲精品日本国产第一区| 99久久精品一区二区三区| 2022亚洲国产成人精品| 国产精品女同一区二区软件| 中文字幕制服av| 国产成人一区二区在线| 自拍偷自拍亚洲精品老妇| 在线观看美女被高潮喷水网站| 久久久久久久国产电影| 欧美日韩综合久久久久久| 亚洲国产欧美人成| 国产久久久一区二区三区| 亚洲精品视频女| 91久久精品国产一区二区三区| 欧美极品一区二区三区四区| 最近中文字幕2019免费版| 97在线视频观看| 久久精品国产a三级三级三级| 国产永久视频网站| 视频中文字幕在线观看| 亚洲精品一区蜜桃| 亚洲真实伦在线观看| 最近中文字幕高清免费大全6| 99热这里只有精品一区| 夜夜看夜夜爽夜夜摸| 少妇的逼好多水| 18禁动态无遮挡网站| 久久久亚洲精品成人影院| 久久久久久久精品精品| 男人添女人高潮全过程视频| 最黄视频免费看| 男人和女人高潮做爰伦理| 久久久久国产网址| 又粗又硬又长又爽又黄的视频| 精品一区二区三区视频在线| 午夜免费男女啪啪视频观看| 国产真实伦视频高清在线观看| 高清不卡的av网站| 欧美97在线视频| 欧美日韩在线观看h| 99久国产av精品国产电影| 亚洲精品一二三| 视频区图区小说| 亚洲国产色片| 免费看光身美女| 亚洲欧洲日产国产| 亚洲不卡免费看| 在线看a的网站| 黄色怎么调成土黄色| 熟妇人妻不卡中文字幕| 久久久久久久久久人人人人人人| 麻豆国产97在线/欧美| 久久精品国产亚洲av天美| 亚洲av欧美aⅴ国产| 中文天堂在线官网| 丰满乱子伦码专区| 韩国av在线不卡| 国产精品一二三区在线看| 国产高清不卡午夜福利| 国产精品一区二区在线不卡| 高清在线视频一区二区三区| 国产精品福利在线免费观看| 大香蕉97超碰在线| 久热这里只有精品99| 亚洲精品日韩在线中文字幕| 欧美xxⅹ黑人| 人体艺术视频欧美日本| 国产精品国产三级国产av玫瑰| 中文字幕精品免费在线观看视频 | 国产成人精品一,二区| 精品少妇黑人巨大在线播放| 最黄视频免费看| videossex国产| 免费人成在线观看视频色| 国产爱豆传媒在线观看| 人人妻人人添人人爽欧美一区卜 | 十八禁网站网址无遮挡 | 91久久精品国产一区二区三区| 日韩一区二区三区影片| 亚洲三级黄色毛片| 人妻夜夜爽99麻豆av| av网站免费在线观看视频| 夫妻性生交免费视频一级片| 日产精品乱码卡一卡2卡三| 亚洲欧美精品自产自拍| 视频中文字幕在线观看| 午夜免费观看性视频| 汤姆久久久久久久影院中文字幕| 高清黄色对白视频在线免费看 | 亚洲一级一片aⅴ在线观看| 亚洲一级一片aⅴ在线观看| 国产精品一区www在线观看| 在线观看免费视频网站a站| 少妇丰满av| 国产在视频线精品| 亚洲,欧美,日韩| 国产一区二区在线观看日韩| 国产成人aa在线观看| 在线播放无遮挡| 日韩av在线免费看完整版不卡| 久久国产亚洲av麻豆专区| a级毛色黄片| 欧美变态另类bdsm刘玥| 精品一区二区免费观看| 五月天丁香电影| 久久久久久久久久久免费av| 深夜a级毛片| 国产免费一区二区三区四区乱码| 欧美日本视频| 在线观看一区二区三区| 亚洲国产色片| 亚洲经典国产精华液单| 一区二区av电影网| 又黄又爽又刺激的免费视频.| 日日啪夜夜撸| a 毛片基地| 精品国产露脸久久av麻豆| 啦啦啦视频在线资源免费观看| 女性生殖器流出的白浆| 免费观看性生交大片5| 亚洲av不卡在线观看| 日韩不卡一区二区三区视频在线| 欧美 日韩 精品 国产| 日韩伦理黄色片| a级毛色黄片| 亚洲精品日韩av片在线观看| 狂野欧美白嫩少妇大欣赏| 午夜日本视频在线| 国产成人aa在线观看| 少妇高潮的动态图| 啦啦啦视频在线资源免费观看| 久久精品熟女亚洲av麻豆精品| 丰满乱子伦码专区| 午夜免费观看性视频| 久久久久人妻精品一区果冻| 80岁老熟妇乱子伦牲交| 少妇人妻久久综合中文| 97在线人人人人妻| 边亲边吃奶的免费视频| 成人黄色视频免费在线看| 欧美日韩视频精品一区| 国产成人a∨麻豆精品| 久久久久久久大尺度免费视频| 欧美精品人与动牲交sv欧美| 自拍欧美九色日韩亚洲蝌蚪91 | 这个男人来自地球电影免费观看 | 日韩制服骚丝袜av| 简卡轻食公司| 内射极品少妇av片p| 黄片wwwwww| 欧美xxxx性猛交bbbb| 亚洲中文av在线| 亚洲精品色激情综合| 建设人人有责人人尽责人人享有的 | 波野结衣二区三区在线| 美女cb高潮喷水在线观看| 午夜福利网站1000一区二区三区| 国产精品嫩草影院av在线观看| 亚洲av欧美aⅴ国产| 国产乱来视频区| 欧美区成人在线视频| 色哟哟·www| 老熟女久久久| 大码成人一级视频| 新久久久久国产一级毛片| 777米奇影视久久| tube8黄色片| 国产真实伦视频高清在线观看| 蜜桃久久精品国产亚洲av| 成人黄色视频免费在线看| 亚洲图色成人| 色婷婷av一区二区三区视频| 男女无遮挡免费网站观看| 免费不卡的大黄色大毛片视频在线观看| 特大巨黑吊av在线直播| 亚洲第一区二区三区不卡| 91在线精品国自产拍蜜月| h视频一区二区三区| 欧美一区二区亚洲| 亚洲av中文字字幕乱码综合| 少妇的逼好多水| 91狼人影院| 精品一区在线观看国产| 蜜桃在线观看..| 黑人猛操日本美女一级片| 免费人成在线观看视频色| 欧美一区二区亚洲| 国产午夜精品一二区理论片| 国产高清有码在线观看视频| 内射极品少妇av片p| 一区二区三区免费毛片| 国产高清不卡午夜福利| 日韩 亚洲 欧美在线| 交换朋友夫妻互换小说| 亚洲va在线va天堂va国产| 国产免费福利视频在线观看| 色婷婷久久久亚洲欧美| 自拍偷自拍亚洲精品老妇| 美女高潮的动态| 一区二区av电影网| 尤物成人国产欧美一区二区三区| 街头女战士在线观看网站| 久久久久性生活片| 国产成人精品久久久久久| 日韩 亚洲 欧美在线| 男女啪啪激烈高潮av片| 国产爽快片一区二区三区| 亚洲精品日韩av片在线观看| 精品亚洲成a人片在线观看 | 久久99热这里只有精品18| av在线老鸭窝| 色婷婷av一区二区三区视频| 国产精品国产av在线观看| 国产男女内射视频| 精品亚洲成a人片在线观看 | 校园人妻丝袜中文字幕| 涩涩av久久男人的天堂| 亚洲av男天堂| 国产成人a区在线观看| 成人漫画全彩无遮挡| 老司机影院毛片| 国产成人a∨麻豆精品| 人人妻人人爽人人添夜夜欢视频 | 少妇高潮的动态图| 女性被躁到高潮视频| 亚洲电影在线观看av| 午夜激情福利司机影院| 成人毛片a级毛片在线播放| 香蕉精品网在线| 欧美亚洲 丝袜 人妻 在线| 久久精品国产亚洲网站| 美女国产视频在线观看| 亚洲av在线观看美女高潮| 麻豆国产97在线/欧美| 国产一区有黄有色的免费视频| 亚洲国产精品一区三区| 欧美高清成人免费视频www| 观看美女的网站| 欧美激情国产日韩精品一区| 噜噜噜噜噜久久久久久91| 三级经典国产精品| 一区二区三区免费毛片| 亚洲av中文av极速乱| 国产 精品1| 少妇被粗大猛烈的视频| 久久久久久久大尺度免费视频| 亚洲国产欧美在线一区| 97在线人人人人妻| 国产高清国产精品国产三级 | av视频免费观看在线观看| 久久 成人 亚洲| 国产精品一区二区性色av| 国产精品嫩草影院av在线观看| 精品人妻视频免费看| 少妇 在线观看| 久久99热这里只有精品18| 久久久欧美国产精品| 欧美成人a在线观看| 美女脱内裤让男人舔精品视频| 久久久久性生活片| 丰满人妻一区二区三区视频av| 国产一级毛片在线| 亚洲欧美精品专区久久| 国产精品嫩草影院av在线观看| 午夜福利高清视频| 中国国产av一级| 夜夜看夜夜爽夜夜摸| 国产无遮挡羞羞视频在线观看| 我的女老师完整版在线观看| 国产免费又黄又爽又色| 国产高清有码在线观看视频| 欧美精品一区二区大全| 久久久久国产精品人妻一区二区| 精品久久久久久电影网| 精品99又大又爽又粗少妇毛片| 亚洲内射少妇av| 女性生殖器流出的白浆| av国产久精品久网站免费入址| 免费高清在线观看视频在线观看| 老师上课跳d突然被开到最大视频| 国产一区亚洲一区在线观看| 一级毛片aaaaaa免费看小| 久久人人爽人人爽人人片va| 久热这里只有精品99| 久久午夜福利片| 一级毛片黄色毛片免费观看视频| 国产精品久久久久成人av| av免费观看日本| 中文字幕免费在线视频6| 男女国产视频网站| 赤兔流量卡办理| 一区二区三区乱码不卡18| 国产精品福利在线免费观看| 国产欧美日韩一区二区三区在线 | 亚洲欧美成人精品一区二区| 国产精品一区二区在线不卡| 天堂中文最新版在线下载| 97精品久久久久久久久久精品| 国产成人免费观看mmmm| av黄色大香蕉| 我要看日韩黄色一级片| 日本-黄色视频高清免费观看| 国产美女午夜福利| 一级毛片 在线播放| 国产探花极品一区二区| 午夜老司机福利剧场| 深夜a级毛片| 中文字幕制服av| 少妇的逼好多水| 日本黄色日本黄色录像| 最近中文字幕高清免费大全6| 午夜福利在线观看免费完整高清在| 日本黄色片子视频| 亚洲av欧美aⅴ国产| 欧美极品一区二区三区四区| 九色成人免费人妻av| 日韩,欧美,国产一区二区三区| 99久久综合免费| 亚洲精品自拍成人| 午夜福利影视在线免费观看| 国产有黄有色有爽视频| a级毛色黄片| 欧美精品一区二区免费开放| 中文字幕久久专区| 国内揄拍国产精品人妻在线| 男女无遮挡免费网站观看| 女人十人毛片免费观看3o分钟| 亚洲,欧美,日韩| av不卡在线播放| 久久国产乱子免费精品| 欧美高清性xxxxhd video| 免费久久久久久久精品成人欧美视频 | 高清视频免费观看一区二区| 97热精品久久久久久|