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

    中國氣候變化預(yù)估數(shù)據(jù)在河北地區(qū)的檢驗分析

    2015-03-12 02:21:04郝立生向亮張婧
    大氣科學(xué)學(xué)報 2015年3期
    關(guān)鍵詞:河北地區(qū)數(shù)值模擬

    郝立生,向亮,張婧

    (1.天津市氣象局,天津 300074;2.中國氣象局 蘭州干旱氣象研究所,甘肅 蘭州 730020;

    3.中國氣象局 干旱氣候變化與減災(zāi)重點開放實驗室,甘肅 蘭州 730020;4.河北省氣候中心,河北 石家莊 050021)

    ?

    中國氣候變化預(yù)估數(shù)據(jù)在河北地區(qū)的檢驗分析

    郝立生1,2,3,向亮4,張婧4

    (1.天津市氣象局,天津 300074;2.中國氣象局 蘭州干旱氣象研究所,甘肅 蘭州 730020;

    3.中國氣象局 干旱氣候變化與減災(zāi)重點開放實驗室,甘肅 蘭州 730020;4.河北省氣候中心,河北 石家莊 050021)

    摘要:應(yīng)用觀測資料對中國地區(qū)預(yù)估數(shù)據(jù)集進行檢驗分析。結(jié)果表明:1)7月降水量,區(qū)域模式在張家口南部桑洋河盆地和北部壩上地區(qū)的模擬值偏大,而在北京東邊興隆高山區(qū)和東南沿海地區(qū)的模擬值明顯偏小。2)7月最高氣溫,在北京東邊興隆高山區(qū)和東部平原的模擬值明顯偏高,而在張家口南部桑洋河盆地和太行山東側(cè)的模擬值偏低。3)1月最低氣溫,在北京東邊興隆高山區(qū)、太行山北段高山區(qū)和東部平原地區(qū)的模擬值偏高,而在張家口南部桑洋河盆地和太行山東側(cè)的模擬值偏低。4)無論是對降水還是氣溫,由于全球模式的空間分辨率偏低、很難描述河北地區(qū)的地形特征,所以模擬結(jié)果很差。由于區(qū)域氣候模式的分辨率提高、對河北地區(qū)地形特征的描述有了改進,所以模擬效果有明顯改善;但受分辨率所限,它對局部地形如興隆高山、桑洋河盆地、太行山北部高山等地形特征的描述不是很好,造成模擬結(jié)果在這些地區(qū)出現(xiàn)系統(tǒng)性偏差,因此應(yīng)用這些數(shù)據(jù)時需加以訂正。

    關(guān)鍵詞:氣候模式;數(shù)值模擬;預(yù)估資料;河北地區(qū);檢驗分析

    0引言

    全球氣候變化對經(jīng)濟社會發(fā)展和生態(tài)環(huán)境有深遠(yuǎn)的影響,已成為21世紀(jì)各國可持續(xù)發(fā)展中面臨的重大問題。未來氣候如何變化?對區(qū)域可持續(xù)發(fā)展會造成什么影響?迫切需要對未來氣候變化進行科學(xué)預(yù)估。為滿足氣候變化影響評估方面的科研和業(yè)務(wù)需求,中國氣象局國家氣候中心先后于2008年10月和2009年11月制作發(fā)布了“中國地區(qū)氣候變化預(yù)估數(shù)據(jù)集”(Version1.0,Version2.0)。這些數(shù)據(jù)在氣候變化影響評估中發(fā)揮了很好的作用(Gao et al.,2008;Xu et al.,2009;郝立生等,2010)。

    為更好地預(yù)估未來氣候變化,政府間氣候變化專門委員會(Intergovernmental Panel on Climate Change,IPCC)對第三次和第四次評估報告采用的排放情景進行了改進(Meinshausen et al.,2011),即“典型濃度路徑”(Representative Concentration Pathways,RCP)。這些排放情景主要包括RCP8.5、RCP6.0、RCP4.5、RCP2.6四種情景(Moss et al.,2010;王紹武等,2012)。RCP2.6情景是指輻射強迫在2100年之前達到峰值,到2100年下降至2.6 W/m2,全球平均溫度上升限制在2.0 ℃之內(nèi)。RCP4.5情景是指輻射強迫穩(wěn)定在4.5 W/m2。RCP6.0情景是指輻射強迫穩(wěn)定在6.0 W/m2。RCP8.5情景是指輻射強迫在2100年上升至8.5 W/m2。為給科研業(yè)務(wù)人員提供一套新排放情景下未來中國區(qū)域氣候變化的預(yù)估數(shù)據(jù),國家氣候中心氣候變化適應(yīng)室在2012年12月制作發(fā)布了“中國地區(qū)氣候變化預(yù)估數(shù)據(jù)集”的第三版(Version3.0),向用戶提供CMIP5(the fifth phase of the CMIP)全球氣候模式以及新版本區(qū)域氣候模式的模擬和預(yù)估數(shù)據(jù)(高學(xué)杰等,2010;石英等,2010;Gao et al.,2012a)。

    為改進氣候模擬和科學(xué)應(yīng)用模擬結(jié)果,一些學(xué)者對全球氣候模式和區(qū)域氣候模式模擬數(shù)據(jù)和預(yù)估效果作了分析評估(高學(xué)杰等,2006;許崇海等,2007;張冬峰等,2007;孫林海和劉一鳴,2008)。張宏芳和陳海山(2011)、Xu and Xu(2012)對第五次評估所用的全球模式結(jié)果作了初步評估分析,Gao et al.(2012b)評估了區(qū)域模式對季風(fēng)降水模擬存在的不確定性。無論是全球模式還是中尺度模式,下邊界地形對模擬和預(yù)測結(jié)果都有顯著影響(李書嚴(yán)和馬京津,2011;郭蕊等,2013;朱司光等,2013;敖翔宇等,2014;朱芳澤等,2014)。為更好地使用數(shù)值模式結(jié)果,對模式模擬能力和預(yù)測結(jié)果進行分析評估必不可少。如朱芳澤等(2014)對全球海氣耦合模式進行評估,發(fā)現(xiàn)許多海氣耦合模式需重點改進對南、北半球深對流海域物理過程的描述。徐敏等(2013)對IPCC AR4的25個全球氣候模式的模擬結(jié)果與觀測資料作對比分析,發(fā)現(xiàn)高分辨率模式在空間分布上模擬效果明顯要好于低分辨率模式。

    如何更好地把中國地區(qū)氣候變化預(yù)估數(shù)據(jù)集(Version3.0)在河北應(yīng)用,目前還缺乏這方面的檢驗分析。河北省地處東亞東部,北接蒙古高原,南連黃淮平原,西倚太行山,東臨渤海(王衛(wèi)等,2012)。全省地勢西北部高,東南部低,由壩上高原、太行山燕山連成的弧形山脈、東南部的平原組成(圖1)。最北部的壩上高原是內(nèi)蒙古高原東南邊緣,平均海拔1 200~1 500 m。高原南邊是東西走向的燕山山地及東北—西南走向的太行山脈。太行山脈南北長約600 km,東西寬約180 km,呈東北—西南向,平均海拔1 000~1 500 m。燕山山脈東西長約330 km,南北寬約120 km,山體孤立,排列不規(guī)則,一般海拔800~1 000 m。山區(qū)東南部為河北平原,平原海拔大多在50 m以下。這種地形對河北省氣候有非常大的影響,特別是對氣溫和降水的影響非常顯著。

    本文選擇較低排放情景RCP4.5時的1961—2050年的全球模式數(shù)據(jù)、區(qū)域模式數(shù)據(jù)與觀測數(shù)據(jù)進行對比,重點針對河北地區(qū)的氣溫、降水的模擬、預(yù)估情況進行檢驗分析。一方面為更科學(xué)地應(yīng)用這套數(shù)據(jù)作氣候變化分析提供參考意見,同時為改進氣候模式模擬技術(shù)提供參考依據(jù)。

    1資料與方法

    主要用到兩種資料:1)國家氣候中心氣候變化適應(yīng)室制作發(fā)布的第三版“中國地區(qū)氣候變化預(yù)估數(shù)據(jù)集”(Version3.0),包括全球氣候模式數(shù)據(jù)和區(qū)域氣候模式數(shù)據(jù)。全球模式數(shù)據(jù)來自于PCMDI(Program for Climate Model Diagnosis and Intercomparison)公開發(fā)布的全球23個模式組提供的全球氣候模式數(shù)值模擬結(jié)果,這些模式數(shù)據(jù)對IPCC第五次評估報告提供了重要支持。國家氣候中心將21個CMIP5全球氣候模式的模擬結(jié)果,經(jīng)過插值計算將其統(tǒng)一降尺度到1°×1°分辨率,利用簡單平均方法進行多模式集合,制作成一套包括歷史模擬數(shù)據(jù)(1901—2005年)和RCP2.6、RCP4.5、RCP8.5排放情景下預(yù)估數(shù)據(jù)(2006—2100年)的月平均資料。區(qū)域模式數(shù)據(jù)是采用國際理論物理研究中心(ICTP,the Abdus Salam International Centre for Theoretical Physics)的區(qū)域氣候模式RegCM4.0(Giorgi et al.,2012)單向嵌套BCC_CSM1.1(Beijing Climate Center_Climate System Model version 1.1)全球氣候系統(tǒng)模式(Wu et al.,2010;Xin et al.,2013),對中國地區(qū)進行模擬和預(yù)估,從而得到歷史模擬數(shù)據(jù)(1961—2005年)和RCP4.5和RCP8.5排放情景下的預(yù)估數(shù)據(jù)(2006—2050年)的月平均資料,空間分辨率0.5°×0.5°(Gao et al.,2013)。本文針對低排放情景RCP4.5下的模擬數(shù)據(jù)進行檢驗分析。2)河北省氣候中心整理的全省72站逐月氣溫、降水資料,時間1961年1月—2013年1月。

    因為社會對夏季降水量、夏季最高氣溫和冬季最低氣溫、冬季平均氣溫變化情況比較關(guān)注,這涉及到防汛、防暑、防寒、采暖工作調(diào)度和安全生產(chǎn)安排等問題。因此,本文重點選擇7月降水量、7月最高氣溫、1月最低氣溫、1月平均氣溫進行檢驗分析。

    采用的方法主要是對比檢驗分析。作區(qū)域平均時間序列分析時,首先將全球模式數(shù)據(jù)和區(qū)域模式數(shù)據(jù)插值到72個站點上(圖1),然后與河北實際觀測資料作對比分析。全省平均指72站平均,多年平均指1976—2005年平均值。序列分析時,降水是相對于1976—2005年平均值的差值百分比變化,氣溫是相對于1976—2005年平均值的距平變化??臻g分析時,偏差顯著性采用t檢驗方法。

    圖1 河北省地形和氣象站分布(彩色陰影表示地形海拔高度,單位:m)Fig.1 Topography and weather stations in Hebei province(The color shadings denote terrain height with the unit of m)

    2檢驗分析

    2.1 夏季(7月)降水量

    2.1.1 時間變化

    河北省受東亞夏季風(fēng)影響,降水高度集中在夏季。夏季(6—8月)降水量約占全年降水量的70%,其中7月降水量最多。圖2是全省平均的7月降水量距平百分比變化曲線。河北1961—2012年實際觀測到的降水距平百分比變化在-30%~40%,變化幅度較大。而全球模式模擬和預(yù)估曲線變化比較平緩,變化幅度在-10%~15%,預(yù)估的未來降水基本呈穩(wěn)定緩慢增加趨勢。模擬值和對2006—2012年預(yù)估值明顯小于實際值,不符合實際情況。因此,全球模式對河北未來降水量變化幅度的預(yù)估值偏小,可信度低。相反,區(qū)域模式模擬和對2006—2012年預(yù)估的變化幅度與觀測到的變化幅度相當(dāng),基本在-30%~40%。比較1961—2005年模擬值與實際值發(fā)現(xiàn),兩者有一致的時候,也有相反的時候,特別是20世紀(jì)90年代后期以來擬合非常接近。2006—2012年降水出現(xiàn)向增多轉(zhuǎn)折的趨勢,預(yù)估結(jié)果與實際也非常一致。因此,區(qū)域模式對未來預(yù)估結(jié)果可信度高。區(qū)域模式預(yù)估未來降水總體偏多,無明顯長期變化趨勢,但階段性明顯,2010—2020年河北降水會偏多,之后減少,2030年前后偏多,之后減少。從變化趨勢和變化幅度看,區(qū)域模式結(jié)果更接近實際,其結(jié)果明顯優(yōu)于全球模式結(jié)果。

    圖2 相對于1976—2005年平均的河北7月降水距平百分率變化(單位:%;虛線是觀測站資料,粗實線是全球模式結(jié)果,細(xì)實線是區(qū)域模式結(jié)果,均為5 a滑動平均結(jié)果)Fig.2 Variation of July precipitation anomaly percentage(%) relative to average for 1976—2005 in Hebei(Dashed line is the observed results from station data,heavy solid line is the simulated results from the global model,and thin solid line is the simulated results from the regional model.All results are 5-yr moving average)

    2.1.2空間分布

    圖3是河北省1976—2005年平均的7月降水量觀測值、模式模擬值及其與觀測值之差空間分布。從圖3a看到,東西走向的燕山地區(qū)由于與東亞夏季風(fēng)風(fēng)向垂直,使得這里成為河北夏季降水量最大值中心;太行山北段東麓由于與副高外圍東南風(fēng)垂直,使這里成為河北夏季降水次大值中心;此外,河北東南部沿海地區(qū)降水量也明顯偏多。河北平原中南部由于處于山東丘陵地形雨影區(qū),成為降水量偏少中心;北部張家口和承德壩上地區(qū),由于偏南夏季風(fēng)攜帶的水汽大量在燕山、太行山迎風(fēng)坡降落,使得其到達這里時水汽來源不足,造成降水量明顯偏少。

    圖3 1976—2005年7月河北降水的空間分布(單位:mm;陰影區(qū)表示通過0.05信度的顯著性檢驗)  a.觀測結(jié)果;b.全球模式模擬結(jié)果;c.全球模式模擬結(jié)果減去觀測結(jié)果;d.區(qū)域模式模擬結(jié)果;e.區(qū)域模式模擬結(jié)果減去觀測結(jié)果Fig.3 Spatial distribution of July rainfall(mm) in Hebei during 1976—2005(Only regions with difference exceeding 0.05 significance level are shaded)  a.observed results;b.global model simulated results;c.global model simulated results minus observed results;d.regional model simulated results;e.regional model simulated results minus observed results

    全球模式模擬的降水量值空間分布差別不大(圖3b),東北部和南部平原略偏多,而燕山南麓和張家口西北降水少。進一步分析發(fā)現(xiàn)(圖3c),對燕山南麓地區(qū)和太行山北段東麓地區(qū)模擬值明顯偏小,可能是由于全球模式空間分辨率低,不能很好地描述燕山、太行山局地地形與夏季風(fēng)的相互作用。對北部壩上地區(qū)模擬降水量值明顯偏大,通過了顯著性檢驗。可見,全球模式模擬結(jié)果的空間分布與實際偏差很大。

    區(qū)域模式模擬結(jié)果空間分布出現(xiàn)了細(xì)微結(jié)構(gòu)(圖3d),對燕山山區(qū)和太行山北段多雨區(qū)、中南部平原的少雨區(qū)基本都模擬出來了。但是,張家口南部桑洋河盆地模擬值明顯偏大,北部壩上地區(qū)模擬值也偏大(圖3e)??梢?區(qū)域模式模擬結(jié)果與實際比較接近,但空間分布形勢整體有些北移。另外,降水量值大小有些偏差,似乎過分強調(diào)地形的阻擋抬升作用而忽視平原降水機制和山脈地形的雨影作用,造成壩上地區(qū)降水模擬值明顯偏大??赡苡捎趨^(qū)域氣候模式空間分辨率還比較低,不能很好描述桑洋河盆地地形,造成這里模擬值顯著偏大,通過了顯著性檢驗。

    可見,全球模式由于空間分辨率低,不能很好描述河北地形特征,模擬效果很差。區(qū)域模式由于分辨率的提高,對河北地形特征描述有了一定改進,模擬效果明顯提高。但由于分辨率的限制,對局部地形如興隆高山、桑洋河盆地、太行山北部高山等地形特征描述還不是很好,造成區(qū)域模式在張家口南部桑洋河盆地地區(qū)、北部壩上地區(qū)存在對降水量模擬值偏大的系統(tǒng)誤差,在北京東部興隆山地、東南沿海地區(qū)存在模擬值明顯偏小的系統(tǒng)誤差,應(yīng)用這些數(shù)據(jù)時需要加以訂正。

    2.2 夏季(7月)最高氣溫

    2.2.1時間變化

    圖4是河北夏季(7月)最高氣溫變化曲線。以1976—2005年平均值為基礎(chǔ),觀測到最高氣溫距平在-1.5~2.0 ℃之間變化,總體上呈波動上升趨勢,2000年前后最高,2005年前后低,之后又升高。全球模式對1961—2005年的模擬值呈波動變化,與實際觀測值有一致的時候也有不一致的時候,變化幅度偏小,沒有模擬出2000年前后的高溫時段,但對2006—2012年的預(yù)估比較接近實際。因此,全球模式模擬值與實際比較接近,但對極端值模擬效果差。區(qū)域模式模擬和預(yù)估變化幅度明顯大于實際變化情況,尤其模擬出的極端高溫值、極端低溫值顯著偏大。因此,區(qū)域模式預(yù)估未來氣溫呈波動變化,無明顯長期升高趨勢。全球模式和區(qū)域模式都沒有模擬出2000年前后的明顯高溫,區(qū)域模式只模擬出了一個較小的高值中心??梢?無論全球模式還是區(qū)域模式,對7月最高氣溫的模擬和預(yù)估效果都不是很好。

    2.2.2空間分布

    圖5是河北省1976—2005年平均的7月最高氣溫觀測值、模式模擬值及與其觀測值之差空間分布。從圖5a看到,河北最高氣溫分布與緯度、地形分布有關(guān),高緯、高海拔地區(qū)氣溫低,低緯、低海拔氣溫高,南部平原地區(qū)最高氣溫差別不大。全省北部地區(qū)南北梯度大,中南部小,全省最高氣溫在25~32 ℃之間,南北相差僅8 ℃。

    全球模式模擬的7月最高氣溫空間分布(圖5b)與觀測(圖5a)數(shù)值大小差別不大,但空間分布形勢差別較大。對北京東邊興隆高山和太行山北段高山模擬值有些偏高(圖5c);對北部壩上地區(qū)模擬值明顯偏大,通過了顯著性檢驗;而對張家口南部桑洋河盆地和太行山東麓模擬值偏低。這些分布都明顯不符合實際,可能是由于模式空間分辨率低,不能很好地描述河北地形特征,把北京東邊興隆高山和太行山北段高山考慮過低,而對張家口南部桑洋河盆地設(shè)置過高的緣故,對太行山焚風(fēng)效應(yīng)也沒有模擬出來。因此,應(yīng)用這些模式數(shù)據(jù)時應(yīng)根據(jù)地形加以訂正。

    圖5 1976—2005年7月河北最高氣溫的空間分布(單位:℃;陰影區(qū)表示通過0.05信度的顯著性檢驗)a.觀測結(jié)果;b.全球模式模擬結(jié)果;c.全球模式模擬結(jié)果減去觀測結(jié)果;d.區(qū)域模式模擬結(jié)果;e.區(qū)域模式模擬結(jié)果減去觀測結(jié)果Fig.5 Spatial distribution of July maximum temperature(℃) in Hebei during 1976—2005(Only regions with difference exceeding 0.05 significance level are shaded)  a.observed results;b.global model simulated results;c.global model simulated results minus observed results;d.regional model simulated results;e.regional model simulated results minus observed results

    區(qū)域模式模擬的空間分布結(jié)果要好一些(圖5d),與實際氣溫空間分布特征接近。在河北北部模擬值接近實際,越往南模擬值偏高越多(圖5e)。尤其對北京東邊興隆高山和東部平原地區(qū)氣溫模擬值明顯偏高,對張家口南部桑洋河盆地和太行山東麓模擬值偏低??梢?區(qū)域氣候模式由于分辨率的提高,對河北地形特征描述有了一定改進,模擬效果顯著提高了。但是,由于模式分辨率的限制,對局地地形如興隆高山、桑洋河盆地描述還不是很好,還需要進一步改進。

    2.3 冬季(1月)最低氣溫和平均氣溫

    2.3.1時間變化

    圖6是河北冬季(1月)最低氣溫和平均氣溫變化曲線。圖6a是1月最低氣溫變化,圖6b是1月平均氣溫變化。以1976—2005年平均值為基礎(chǔ),觀測到的最低氣溫距平在-2.0~1.5 ℃之間變化,20世紀(jì)90年代、21世紀(jì)00年代偏高,2000年前后偏低,與7月最高氣溫變化趨勢相反。2010年以來顯著偏低,為近50 a來最低值。全球模式對1961—2005年的模擬值呈弱的波動變化,幅度明顯小于實際觀測值變化,對2006—2012年的預(yù)估與實際偏差很大。因此,全球模式預(yù)估至2050年逐漸升高的趨勢,可信度低。區(qū)域氣候模式對1961—2005年的模擬值呈波動變化,與實際有些偏差,但對2006—2012年最低氣溫轉(zhuǎn)折趨勢預(yù)估接近實際,區(qū)域氣候模式模擬和預(yù)估效果較好。因此,區(qū)域模式預(yù)估未來1月最低氣溫在波動中逐漸升高,可信度高。

    圖6 相對于1976—2005年平均的河北1月最低氣溫(a)和平均氣溫(b)距平的變化(單位:℃;虛線是觀測站資料,粗實線是全球模式結(jié)果,細(xì)實線是區(qū)域模式結(jié)果,均為5 a滑動平均結(jié)果)Fig.6 Variations of January (a)minimum temperature and (b)mean temperature anomalies(℃) relative to average for 1976—2005 in Hebei(Dashed line is the observed results from station data,heavy solid line is the simulated results from the global model,and thin solid line is the simulated results from the regional model.All results are 5-yr moving average)

    1月平均氣溫觀測、模擬、預(yù)估的變化趨勢同最低氣溫變化基本一致。

    2.3.2空間分布

    對比河北省1976—2005年1月最低氣溫觀測值、模式模擬值及其與觀測值之差空間分布(圖略),可以看到,河北最低氣溫分布與緯度、地形分布有關(guān),高緯、高海拔地區(qū)氣溫低,低緯、低海拔氣溫高,北部山區(qū)和壩上高原最低氣溫南北梯度很大,與7月最高氣溫分布明顯不同。桑洋河盆地在7月是暖中心,而在1月為冷中心。全省最低氣溫空間分布在-5~-24 ℃之間,南北溫差達19 ℃,比夏季南北溫差大的多。

    全球模式模擬1月最低氣溫空間分布無論是分布形勢還是數(shù)值大小都與實際偏差很大(圖略)。對燕山山區(qū)、太行山北段高山區(qū)和壩上高原模擬值明顯偏高,對桑洋河盆地、太行山南段及平原模擬值偏低(圖略)。造成這種偏差的原因可能是由于模式空間分辨率限制,不能很好描述河北地形特征,更無法模擬太行山的焚風(fēng)效應(yīng),結(jié)果在一些地區(qū)造成了模擬系統(tǒng)性偏差。

    區(qū)域模式模擬結(jié)果要好很多,空間分布形勢和數(shù)值大小比較符合實際的空間分布特征(圖略)。但對北京東邊興隆高山、太行山北段高山和張家口壩上地區(qū)模擬值顯著偏高,對張家口南部桑洋河盆地和太行山南段模擬值偏低(圖略)。這說明隨著區(qū)域氣候模式分辨率的提高,對河北地形描述有了明顯改進,模擬結(jié)果有了很大提高。但由于分辨率的限制,沒有很好地描述局部地形如北京東邊興隆高山、桑洋河盆地地形特征,也沒有模擬出太行山地形造成的焚風(fēng)效應(yīng),結(jié)果在這些局部地區(qū)出現(xiàn)了系統(tǒng)性偏差,應(yīng)用這些模式數(shù)據(jù)時需要根據(jù)地形特征加以訂正。

    1月平均氣溫的空間分布、模擬、預(yù)估空間分布形勢和存在的問題與最低氣溫分布形勢相似(圖略),存在的問題也是一樣的。

    3結(jié)論和討論

    就整體變化,區(qū)域氣候模式對河北7月降水的模擬明顯好于全球氣候模式結(jié)果,特別是對2006—2012年降水出現(xiàn)向增多轉(zhuǎn)折的趨勢預(yù)估結(jié)果非常好。區(qū)域模式和全球模式對7月最高氣溫模擬都不是很好,區(qū)域模式數(shù)據(jù)對1月最低氣溫預(yù)估效果好一些,尤其對2006—2012年最低氣溫出現(xiàn)偏低的轉(zhuǎn)折趨勢預(yù)估結(jié)果非常好。

    對于河北7月降水空間分布,區(qū)域氣候模式模擬情況明顯好于全球模式結(jié)果。由于模式空間分辨率的限制,不能很好地描述局部地形特征,造成區(qū)域模式在張家口南部桑洋河盆地地區(qū)、北部壩上地區(qū)存在模擬值偏大的系統(tǒng)偏差,在北京東邊興隆高山區(qū)、東南沿海地區(qū)存在模擬值明顯偏小的系統(tǒng)偏差。

    對于河北7月最高氣溫空間分布,全球模式由于空間分辨率偏低,很難描述河北地區(qū)地形特征,模擬結(jié)果很差。區(qū)域氣候模式由于分辨率的提高,對河北地形特征描述有了一定改進,模擬效果明顯提高。但由于模式分辨率的限制,對局地地形如興隆高山區(qū)、桑洋河盆地地形描述還不是很好,在北京東邊興隆高山和東部平原地區(qū)模擬值明顯偏高、在張家口南部桑洋河盆地和太行山東麓模擬值偏低的系統(tǒng)偏差。

    對于河北1月最低氣溫空間分布,區(qū)域氣候模式模擬情況明顯好于全球模式結(jié)果,但由于分辨率的限制,不能很好地描述局部地形特征,造成在北京東邊興隆高山區(qū)、太行山北段高山區(qū)和東部平原地區(qū)模擬值偏高、張家口南部桑洋河盆地和太行山東側(cè)模擬值偏低的系統(tǒng)偏差。

    無論是對降水還是氣溫,全球模式由于空間分辨率偏低,很難描述河北地區(qū)地形特征,模擬結(jié)果很差。區(qū)域氣候模式由于分辨率的提高,對河北地形特征描述有了一定改進,模擬效果明顯提高。但由于分辨率的限制,對局部地形如興隆高山、桑洋河盆地、太行山北部高山等地形特征描述還不是很好,造成模擬值在這些地區(qū)出現(xiàn)系統(tǒng)偏差,應(yīng)用這些數(shù)據(jù)時需要加以訂正。

    致謝:本研究所使用的全球氣候模式氣候變化預(yù)估數(shù)據(jù)是由國家氣候中心研究人員對數(shù)據(jù)進行的整理、分析和惠許使用。原始數(shù)據(jù)由各模式組提供,由WGCM(JSC/CLIVAR Working Group on Coupled Modelling)組織PCMDI(Program for Climate Model Diagnosis and Intercomparison)搜集歸類。多模式數(shù)據(jù)集的維護由美國能源部科學(xué)辦公室提供資助。感謝國家氣候中心提供了利用區(qū)域氣候模式所進行的中國區(qū)域未來氣候變化模擬結(jié)果。

    參考文獻(References):

    敖翔宇,任雪娟,湯劍平,等.2014.長江三角洲城市群對夏季降水影響機制的模擬研究[J].氣象科學(xué),34(6):684-691.Ao Xiangyu,Ren Xuejuan,Tang Jianping,et al.2014.Simulation on mechanism of urbanization effects on summer precipitation over the Yangtze River Delta[J].J Meteor Sci,34(6):684-691.(in Chinese).

    高學(xué)杰,徐影,趙宗慈,等.2006.數(shù)值模式不同分辨率和地形對東亞降水模擬影響的試驗[J].大氣科學(xué),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 J Atmos Sci,30(2):185-192.(in Chinese).

    Gao X J,Shi Y,Song R Y,et al.2008.Reduction of future monsoon precipitation over China:Comparison between a high resolution RCM simulation and the driving GCM[J].Meteor Atmos Phys,100(1):73-86.

    高學(xué)杰,石英,Giorgi F.2010.中國區(qū)域氣候變化的一個高分辨率數(shù)值模擬[J].中國科學(xué)D輯,40(7):911-922.Gao Xuejie,Shi Ying,Giorgi F.2010.A high resolution simulation of climate change over China[J].Sci China(Earth Sci),40(7):911-922.(in Chinese).

    Gao X J,Shi Y,Zhang D F,et al.2012a.Climate change in China in the 21st century as simulated by a high resolution regional climate model[J].Chin Sci Bull,57(10):1188-1195.

    Gao X J,Shi Y,Zhang D F,et al.2012b.Uncertainties in monsoon precipitation projections over China:Results from two high-resolution RCM simulations[J].Climate Res,52(1):213-226.

    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].Atmospheric and Oceanic Science Letters,6(5):381-386.

    Giorgi F,Coppola E,Solmon F,et al.2012.RegCM4:Model description and preliminary tests over multiple CORDEX domains[J].Climate Res,52(1):7-29.

    郭蕊,苗春生,張楠.2013.一次淮河流域梅雨鋒暴雨的大別山地形敏感性試驗[J].大氣科學(xué)學(xué)報,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,36(5):626-634.(in Chinese).

    郝立生,閔錦忠,劉克巖.2010.氣候變化對河北省水資源總量的影響[J].河北師范大學(xué)學(xué)報:自然科學(xué)版,34(4):491-496.Hao Lisheng,Min Jinzhong,Liu Keyan.2010.Impact of climate change on Hebei water resource[J].Journal of Hebei Normal University(Natural Science Edition),34(4):491-496.(in Chinese).

    李書嚴(yán),馬京津.2011.城市化進程對北京地區(qū)降水的影響分析[J].氣象科學(xué),31(4):414-421.Li Shuyan,Ma Jingjin.2011.Impact of urbanization on precipitation in Beijing area[J].J Meteor Sci,31(4):414-421.(in Chinese)

    Meinshausen M,Smith S J,Calvin K,et al.2011.The RCP greenhouse gas concentrations and their extensions from 1765 to 2300[J].Climatic Change,109(1):213-241.

    Moss R H,Edmonds J A,Hibbard K A,et al.2010.The next generation of scenarios for climate change research and assessment[J].Nature,463(7282):747-756.

    石英,高學(xué)杰,吳佳,等.2010.華北地區(qū)未來氣候變化的高分辨率數(shù)值模擬[J].應(yīng)用氣象學(xué)報,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].J Appl Meteor Sci,21(5):580-589.(in Chinese).

    孫林海,劉一鳴.2008.區(qū)域氣候模式對中國夏季平均氣溫和降水的評估分析[J].氣象,34(11):31-39.Sun Linhai,Liu Yiming.2008.Assessment analysis of summer temperature and rainfall over China from regional climate model[J].Meteor Mon,34(11):31-39.(in Chinese).

    王紹武,羅勇,趙宗慈,等.2012.新一代溫室氣體排放情景[J].氣候變化研究進展,8(4):305-307.Wang Shaowu,Luo Yon,Zhao Zongci,et al.2010.New generation of scenarios of greenhouse gas emission[J].Advances in Climate Change Research,8(4):305-307.(in Chinese).

    王衛(wèi),馮忠江,陳輝.2012.河北地理[M].北京:北京師范大學(xué)出版社.Wang Wei,Feng Zhongjiang,Chen Hui.2012.Hebei Geography[M].Beijing:Beijing Normal University Press.(in Chinese).

    Wu T W,Yu R C,Zhang F,et al.2010.The Beijing climate center atmospheric general circulation model:Description and its performance for the present-day climate[J].Climate Dyn,34(1):123-147.

    Xin X G,Wu T W,Li J L,et al.2013.How well does BCC_CSM1.1 reproduce the 20th century climate change over China?[J].Atmospheric and Oceanic Science Letters,6(1):21-26.

    許崇海,沈新勇,徐影.2007.IPCC AR4模式對東亞地區(qū)氣候模擬能力的分析[J].氣候變化研究進展,3(5):287-292.Xu Chonghai,Shen Xinyong,Xu Ying.2007.An analysis of climate change in East Asia by using the IPCC AR4 simulations[J].Advances in Climate Change Research,3(5):287-292.(in Chinese).

    徐敏,于庚康,徐經(jīng)緯,等.2013.不同分辨率下全球氣候模式的集合結(jié)果檢驗[J].氣象科學(xué),33(2):160-167.Xu Min,Yu Gengkang,Xu Jingwei,et al.2013.Assessment of global climate models based on different resolution sets of simulation results[J].J Meteor Sci,33(2):160-167.(in Chinese).

    Xu Y,Xu C H.2012.Preliminary assessment pf simulation of climate changes over China by CMIP5 multi-models[J].Atmospheric and Oceanic Science Letters,5(6):489-494.

    Xu Y,Gao X J,Shen Y,et al.2009.A daily temperature dataset over China and its application in validating a RCM simulation[J].Adv Atmos Sci,26(4):763-772.

    張冬峰,歐陽里程,高學(xué)杰.2007.RegCM3對東亞環(huán)流和中國氣候模擬能力的檢驗[J].熱帶氣象學(xué)報,23(5):444-452.Zhang Dongfeng,Ouyang Licheng,Gao Xuejie,et al.2007.Simulation of the atmospheric circulation over East Asia and climatr in China by RegCM [J].J Trop Meteor,23(5):444-452.(in Chinese).

    張宏芳,陳海山.2011.21個氣候模式對東亞夏季環(huán)流模擬的評估I:氣候態(tài)[J].氣象科學(xué),319(2):119-128.Zhang Hongfang,Chen Haishan.2011.Evaluation of summer circulation simulation over East Asia by 21 climate models.Part I:Climatology[J].J Meteor Sci,31(2):119-128.(in Chinese).

    朱芳澤,王召民,劉成彥,等.2014.BCC_CSM對全球海表溫度和混合層深度的模擬評估[J].大氣科學(xué)學(xué)報,37(6):684-694.Zhu Fangze,Wang Zhaomin,Liu Chengyan,et al.2014.An assessment of the global sea surface temperature and mixed layer depth simulated in BCC_CSM[J].Trans Atmos Sci,37(6):684-694.(in Chinese).

    朱司光,陳海山,周晶.2013.NCAR_CLM系列模式對全球近50 a陸面狀況的模擬及其比較分析[J].大氣科學(xué)學(xué)報,36(4):434-446.Zhu Siguang,Chen Haishan,Zhou Jing.2013.Simulations of global land surface conditions in recent 50 years with three versions of NCAR Community Land Models and their comparative analysis[J].Trans Atmos Sci,36(4):434-446.(in Chinese).

    (責(zé)任編輯:張福穎)

    A test analysis on China climate change projection data in Hebei area

    HAO Li-sheng1,2,3,XIANG Liang4,ZHANG Jing4

    (1.Tianjin Meteorological Bureau,Tianjin 300074,China;2.Institute of Arid Meteorology,CMA,Lanzhou 730020,China; 2.Key Open Laboratory of Arid Climate Change and Disaster Redution,CMA,Lanzhou 730020,China; 4.Hebei Climate Center,Shijiazhuang 050021,China)

    Abstract:Using the station data in Hebei area,a test analysis is maded on simulated results from China climate change projection data provided by National Climate Center.Results show that:1)July precipitation simulated by the regional model(RegMV) exists systematic deviation,which is larger in southern Zhangjiakou-Sangyang River basin and northern Zhangjiakou-Dam Plateau,but significantly smaller in Xinglong mountains east to Beijing and southeast coastal areas.2)July maximum temperature simulated by RegMV exists systematic deviation,which is significantly higher in Xinglong mountains east to Beijing and eastern plain region,but lower in southern Zhangjiakou-Sangyang River basin and eastern Taihang Mountains.3)January minimum temperature simulated by RegMV exists systematic deviation,which is higher in Xinglong mountains east to Beijing,North Taihang Mountains,and eastern plain region,but lower in southern Zhangjiakou-Sangyang River basin and eastern Taihang Mountains.4)The global model is difficult to describe terrain features of Hebei due to its low spatial resolution.Whether precipitation or temperature is simulated,the results are poor.The regional climate model has a certain improvement in describing terrain features of Hebei by improving the resolution,so its simulations have been improved significantly.However,due to the limit of the spatial resolution,such local terrain as Xinglong Mountains,Sangyang River basin,northern Taihang Mountains and other terrain are not well characterized,so the simulated results appear systematic deviation in these areas.It needs an improvement for the application of this set of data in Hebei area.

    Key words:climate model;numerical simulation;projection data;Hebei area;test analysis

    doi:10.13878/j.cnki.dqkxxb.20130303002

    文章編號:1674-7097(2015)03-0362-09

    中圖分類號:P467

    文獻標(biāo)志碼:A

    通信作者:郝立生,博士,正研級高工,研究方向為華北旱澇演變機理及預(yù)測技術(shù)研究,hls54515@163.com.

    基金項目:中國氣象局蘭州干旱氣象研究所干旱氣象科學(xué)研究基金項目(IAM201201);國家氣候中心短期氣候預(yù)測創(chuàng)新團隊基金

    收稿日期:2013-03-03;改回日期:2013-06-03

    郝立生,向亮,張婧.2015.中國氣候變化預(yù)估數(shù)據(jù)在河北地區(qū)的檢驗分析[J].大氣科學(xué)學(xué)報,38(3):362-370.doi:10.13878/j.cnki.dqkxxb.20130303002.

    Hao Li-sheng,Xiang Liang,Zhang Jing.2015.A test analysis on China climate change projection data in Hebei area[J].Trans Atmos Sci,38(3):362-370.(in Chinese)

    猜你喜歡
    河北地區(qū)數(shù)值模擬
    河北高校創(chuàng)業(yè)指導(dǎo)工作現(xiàn)狀及對策研究
    商情(2018年11期)2018-03-30 03:37:42
    河北地區(qū)旅游開發(fā)中的文化遺失現(xiàn)象及對策
    河北地區(qū)旅游開發(fā)中的文化遺失現(xiàn)象及對策
    河北地區(qū)玉米種植中的病蟲害防治策略探討
    張家灣煤礦巷道無支護條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    一種基于液壓緩沖的減震管卡設(shè)計與性能分析
    科技視界(2016年20期)2016-09-29 11:08:27
    超碰97精品在线观看| 一级毛片电影观看| 啦啦啦啦在线视频资源| 国产男女超爽视频在线观看| 亚洲18禁久久av| 国产日韩欧美在线精品| 在线观看一区二区三区| 最近2019中文字幕mv第一页| 美女脱内裤让男人舔精品视频| 免费大片18禁| 国产在线男女| 91久久精品国产一区二区三区| 成人午夜精彩视频在线观看| 我的女老师完整版在线观看| 免费观看在线日韩| 国产高清三级在线| 亚洲欧洲国产日韩| 人体艺术视频欧美日本| 白带黄色成豆腐渣| 欧美一级a爱片免费观看看| 别揉我奶头 嗯啊视频| 国产69精品久久久久777片| 又爽又黄无遮挡网站| 国产v大片淫在线免费观看| 免费大片18禁| 一二三四中文在线观看免费高清| 色吧在线观看| 1000部很黄的大片| 看黄色毛片网站| 偷拍熟女少妇极品色| 国产午夜精品久久久久久一区二区三区| 亚洲人成网站在线观看播放| 亚洲成人一二三区av| 欧美成人a在线观看| 观看美女的网站| 国产成人精品婷婷| 看非洲黑人一级黄片| 久久99精品国语久久久| 身体一侧抽搐| 亚洲欧美精品自产自拍| 国产白丝娇喘喷水9色精品| 成人亚洲精品一区在线观看 | 国内精品美女久久久久久| 国产免费又黄又爽又色| 精品久久久久久久人妻蜜臀av| 18+在线观看网站| 老师上课跳d突然被开到最大视频| 亚洲av中文字字幕乱码综合| 男女那种视频在线观看| 干丝袜人妻中文字幕| 麻豆av噜噜一区二区三区| 又爽又黄无遮挡网站| 欧美3d第一页| 又爽又黄a免费视频| 极品教师在线视频| 99热这里只有精品一区| 久久精品夜色国产| 美女内射精品一级片tv| 男人和女人高潮做爰伦理| 久久韩国三级中文字幕| 亚洲精品国产成人久久av| 国产黄色免费在线视频| 国产女主播在线喷水免费视频网站 | 好男人视频免费观看在线| 欧美日韩国产mv在线观看视频 | 亚洲国产日韩欧美精品在线观看| 亚洲,欧美,日韩| 伊人久久精品亚洲午夜| 黄色配什么色好看| 国产在视频线在精品| 如何舔出高潮| 国产在视频线在精品| 高清视频免费观看一区二区 | 午夜老司机福利剧场| 久久精品国产亚洲网站| 国产成人精品一,二区| 在线播放无遮挡| 国产精品不卡视频一区二区| 久久久欧美国产精品| 精品熟女少妇av免费看| 精品一区二区免费观看| 天天躁日日操中文字幕| 国产成人a∨麻豆精品| 久久久久九九精品影院| 99热这里只有是精品在线观看| 久久久久久久久久久丰满| 午夜免费男女啪啪视频观看| 日韩av不卡免费在线播放| www.色视频.com| 视频中文字幕在线观看| 美女高潮的动态| 波多野结衣巨乳人妻| 亚洲第一区二区三区不卡| 又大又黄又爽视频免费| 日日摸夜夜添夜夜爱| 久久精品久久久久久久性| 成人一区二区视频在线观看| 成人毛片a级毛片在线播放| av免费观看日本| 舔av片在线| 日韩亚洲欧美综合| 亚洲国产精品国产精品| 亚洲色图av天堂| 日韩欧美精品v在线| 久久久久久国产a免费观看| 好男人在线观看高清免费视频| videos熟女内射| 女人久久www免费人成看片| 午夜福利在线观看免费完整高清在| 3wmmmm亚洲av在线观看| 欧美性猛交╳xxx乱大交人| 久久久久网色| 美女被艹到高潮喷水动态| 人体艺术视频欧美日本| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本猛色少妇xxxxx猛交久久| 性插视频无遮挡在线免费观看| 欧美3d第一页| 特级一级黄色大片| 亚洲乱码一区二区免费版| 国产国拍精品亚洲av在线观看| 性色avwww在线观看| 亚洲精品国产av蜜桃| 中文资源天堂在线| 日韩欧美精品v在线| 精品久久久久久电影网| 能在线免费看毛片的网站| 亚洲人成网站在线播| 看黄色毛片网站| 日本黄大片高清| 亚洲欧洲国产日韩| 国精品久久久久久国模美| 免费黄色在线免费观看| 国产色爽女视频免费观看| 国产爱豆传媒在线观看| 综合色av麻豆| 欧美3d第一页| 免费大片黄手机在线观看| 亚洲精品视频女| 国产国拍精品亚洲av在线观看| 午夜精品在线福利| 日韩欧美一区视频在线观看 | 插逼视频在线观看| 国产伦理片在线播放av一区| 国产国拍精品亚洲av在线观看| 免费看美女性在线毛片视频| 婷婷色综合www| 国产精品熟女久久久久浪| 午夜激情欧美在线| 欧美丝袜亚洲另类| 亚洲伊人久久精品综合| 久久久精品94久久精品| 国产人妻一区二区三区在| 又黄又爽又刺激的免费视频.| 国产淫语在线视频| 免费无遮挡裸体视频| 2021天堂中文幕一二区在线观| 日韩不卡一区二区三区视频在线| 国产av不卡久久| 大片免费播放器 马上看| 亚洲国产最新在线播放| 成年人午夜在线观看视频 | 日本免费a在线| 午夜福利在线观看吧| 成人av在线播放网站| 久久精品国产亚洲网站| 久久久久久久国产电影| 久久久久久久亚洲中文字幕| 赤兔流量卡办理| 伊人久久精品亚洲午夜| 日本免费在线观看一区| ponron亚洲| 26uuu在线亚洲综合色| 日本wwww免费看| 欧美日韩在线观看h| 男女下面进入的视频免费午夜| 精品久久国产蜜桃| 国产精品一二三区在线看| av在线亚洲专区| freevideosex欧美| 嫩草影院精品99| 少妇熟女欧美另类| 九草在线视频观看| 99久久中文字幕三级久久日本| 日本爱情动作片www.在线观看| 晚上一个人看的免费电影| 国产午夜精品论理片| 亚洲精品日韩在线中文字幕| 偷拍熟女少妇极品色| 青春草国产在线视频| 亚洲无线观看免费| 蜜臀久久99精品久久宅男| 国产 亚洲一区二区三区 | 日本av手机在线免费观看| 非洲黑人性xxxx精品又粗又长| 国产午夜精品一二区理论片| 国产亚洲5aaaaa淫片| 国产综合精华液| 国产精品久久久久久久久免| 黄色一级大片看看| 免费av毛片视频| 日韩欧美精品免费久久| 亚洲在线观看片| 男的添女的下面高潮视频| 色综合站精品国产| 中文字幕亚洲精品专区| 欧美三级亚洲精品| 99九九线精品视频在线观看视频| 一级av片app| 免费大片黄手机在线观看| 热99在线观看视频| 人妻少妇偷人精品九色| 91狼人影院| 欧美bdsm另类| 亚洲精品影视一区二区三区av| 直男gayav资源| 大话2 男鬼变身卡| a级毛片免费高清观看在线播放| 三级毛片av免费| 一个人看视频在线观看www免费| 成人国产麻豆网| 国产亚洲av嫩草精品影院| 夜夜爽夜夜爽视频| 少妇人妻一区二区三区视频| 国模一区二区三区四区视频| 亚洲人成网站在线观看播放| 亚洲自拍偷在线| 亚洲aⅴ乱码一区二区在线播放| 亚洲美女搞黄在线观看| 插逼视频在线观看| av在线天堂中文字幕| 好男人在线观看高清免费视频| 午夜免费观看性视频| 国产av在哪里看| 日本爱情动作片www.在线观看| 久久精品国产亚洲av天美| 欧美日韩一区二区视频在线观看视频在线 | 亚洲图色成人| 18+在线观看网站| 午夜日本视频在线| 亚洲精品国产成人久久av| 嫩草影院入口| 一级毛片aaaaaa免费看小| 国产精品国产三级专区第一集| 日本爱情动作片www.在线观看| 亚洲精品一区蜜桃| 国产爱豆传媒在线观看| 一个人观看的视频www高清免费观看| 波多野结衣巨乳人妻| 一级毛片 在线播放| 超碰97精品在线观看| 日韩欧美 国产精品| 亚洲aⅴ乱码一区二区在线播放| 午夜视频国产福利| 天堂√8在线中文| 亚洲不卡免费看| 国产精品国产三级专区第一集| 日本一二三区视频观看| 一个人看的www免费观看视频| 久久久精品欧美日韩精品| 免费黄频网站在线观看国产| 亚洲国产精品专区欧美| 成人亚洲精品av一区二区| 日韩av在线免费看完整版不卡| 最近手机中文字幕大全| h日本视频在线播放| xxx大片免费视频| 亚洲四区av| 久久这里有精品视频免费| 亚洲av日韩在线播放| 超碰av人人做人人爽久久| 麻豆乱淫一区二区| 国产毛片a区久久久久| 在线免费观看的www视频| 麻豆乱淫一区二区| 成人亚洲精品一区在线观看 | 国产日韩欧美在线精品| av免费观看日本| 毛片女人毛片| 99久久九九国产精品国产免费| 国产成人免费观看mmmm| 麻豆av噜噜一区二区三区| 99视频精品全部免费 在线| 97在线视频观看| 亚洲精品国产av蜜桃| 永久网站在线| 女人十人毛片免费观看3o分钟| 最近最新中文字幕大全电影3| 亚洲av国产av综合av卡| 久久精品久久精品一区二区三区| 国产色婷婷99| 免费看a级黄色片| 久久97久久精品| 真实男女啪啪啪动态图| 国产黄频视频在线观看| 丰满少妇做爰视频| 中文字幕av在线有码专区| 黄片无遮挡物在线观看| 亚洲四区av| 亚洲在久久综合| 日本wwww免费看| 精品人妻熟女av久视频| 大香蕉久久网| 亚洲国产欧美人成| videos熟女内射| 2022亚洲国产成人精品| 毛片女人毛片| 亚洲经典国产精华液单| 嫩草影院精品99| 亚洲精品一区蜜桃| 国产亚洲5aaaaa淫片| 好男人视频免费观看在线| 波多野结衣巨乳人妻| 日韩精品有码人妻一区| 五月天丁香电影| 看黄色毛片网站| 搡女人真爽免费视频火全软件| 大片免费播放器 马上看| 少妇的逼好多水| 婷婷色av中文字幕| 日本免费在线观看一区| 人妻夜夜爽99麻豆av| 欧美日韩在线观看h| 麻豆国产97在线/欧美| 毛片一级片免费看久久久久| 国产精品蜜桃在线观看| 看黄色毛片网站| 黄色欧美视频在线观看| 身体一侧抽搐| a级毛片免费高清观看在线播放| 性插视频无遮挡在线免费观看| 99热这里只有精品一区| 国产毛片a区久久久久| 色网站视频免费| 国产精品麻豆人妻色哟哟久久 | 麻豆av噜噜一区二区三区| 97热精品久久久久久| 欧美日韩亚洲高清精品| 赤兔流量卡办理| 午夜日本视频在线| 男女国产视频网站| 国产精品久久久久久精品电影小说 | 黑人高潮一二区| 亚洲四区av| 天堂av国产一区二区熟女人妻| 日产精品乱码卡一卡2卡三| 国产精品日韩av在线免费观看| 水蜜桃什么品种好| a级毛色黄片| 真实男女啪啪啪动态图| 亚洲性久久影院| 亚洲欧美一区二区三区黑人 | 最近2019中文字幕mv第一页| 日韩在线高清观看一区二区三区| 国产 亚洲一区二区三区 | 国产精品熟女久久久久浪| 搞女人的毛片| 校园人妻丝袜中文字幕| 人妻制服诱惑在线中文字幕| 搡老妇女老女人老熟妇| 联通29元200g的流量卡| av免费观看日本| 中文乱码字字幕精品一区二区三区 | 99视频精品全部免费 在线| 精品一区在线观看国产| 久久精品国产亚洲av天美| 精品人妻偷拍中文字幕| 大片免费播放器 马上看| 丝袜喷水一区| 日韩av在线免费看完整版不卡| 国产高清有码在线观看视频| 国产高清三级在线| 亚洲欧美成人精品一区二区| 女人久久www免费人成看片| 免费高清在线观看视频在线观看| 久久精品国产自在天天线| 国产亚洲午夜精品一区二区久久 | 国产在线男女| 午夜福利在线观看免费完整高清在| 亚洲欧美精品专区久久| 精品久久久久久久久av| 成年人午夜在线观看视频 | 亚洲在线自拍视频| 亚洲国产高清在线一区二区三| 少妇的逼好多水| 非洲黑人性xxxx精品又粗又长| av天堂中文字幕网| 69av精品久久久久久| 久99久视频精品免费| 精品人妻偷拍中文字幕| 国产中年淑女户外野战色| 91午夜精品亚洲一区二区三区| 男人爽女人下面视频在线观看| 国产精品一区二区性色av| 成年人午夜在线观看视频 | www.色视频.com| 干丝袜人妻中文字幕| 欧美xxⅹ黑人| 天天躁日日操中文字幕| kizo精华| 日本猛色少妇xxxxx猛交久久| 国产淫语在线视频| 91久久精品国产一区二区三区| 免费看不卡的av| 亚洲av一区综合| 黄色配什么色好看| 一本久久精品| 国内少妇人妻偷人精品xxx网站| 亚洲精品第二区| 国产精品福利在线免费观看| 成人亚洲精品av一区二区| 九九久久精品国产亚洲av麻豆| 老师上课跳d突然被开到最大视频| 日韩大片免费观看网站| 亚洲丝袜综合中文字幕| 一级二级三级毛片免费看| 国产成人freesex在线| 少妇猛男粗大的猛烈进出视频 | 国产成人91sexporn| 中文在线观看免费www的网站| 精品午夜福利在线看| 午夜免费观看性视频| 噜噜噜噜噜久久久久久91| 日韩成人伦理影院| 久久97久久精品| 三级男女做爰猛烈吃奶摸视频| 婷婷色综合www| 纵有疾风起免费观看全集完整版 | 免费观看在线日韩| 欧美xxxx性猛交bbbb| 亚洲真实伦在线观看| 超碰97精品在线观看| 国内揄拍国产精品人妻在线| 国产综合懂色| 日韩av免费高清视频| 亚洲人成网站在线观看播放| 91在线精品国自产拍蜜月| 久久久久久久久久久免费av| 日日干狠狠操夜夜爽| 人妻一区二区av| 国产精品三级大全| 亚洲国产欧美人成| 亚洲欧美一区二区三区国产| 三级经典国产精品| 青青草视频在线视频观看| 一本久久精品| 熟妇人妻不卡中文字幕| 尤物成人国产欧美一区二区三区| 国产三级在线视频| 成人综合一区亚洲| 中文字幕制服av| 男的添女的下面高潮视频| 18禁在线无遮挡免费观看视频| 欧美变态另类bdsm刘玥| 久久久国产一区二区| 免费黄频网站在线观看国产| 热99在线观看视频| 亚洲av电影不卡..在线观看| 亚洲成人精品中文字幕电影| 97精品久久久久久久久久精品| 国产一区二区三区综合在线观看 | 精品人妻熟女av久视频| 欧美xxⅹ黑人| 亚洲精品国产av成人精品| 国产黄片美女视频| 一区二区三区四区激情视频| 久久久久免费精品人妻一区二区| 国产成人福利小说| 天堂网av新在线| 777米奇影视久久| 97超碰精品成人国产| 国产高清有码在线观看视频| 成年女人看的毛片在线观看| 99热这里只有是精品在线观看| 亚洲欧洲国产日韩| 成人一区二区视频在线观看| 一夜夜www| 欧美日韩综合久久久久久| 国产精品一区www在线观看| 只有这里有精品99| 精品人妻一区二区三区麻豆| 777米奇影视久久| 中国美白少妇内射xxxbb| 性插视频无遮挡在线免费观看| 大话2 男鬼变身卡| 婷婷色av中文字幕| 爱豆传媒免费全集在线观看| 国产伦理片在线播放av一区| 亚洲国产最新在线播放| 午夜日本视频在线| 直男gayav资源| 精品久久国产蜜桃| 99久久九九国产精品国产免费| 亚洲av免费高清在线观看| 日本与韩国留学比较| 九九爱精品视频在线观看| 国产亚洲91精品色在线| 亚洲精品一二三| 久久韩国三级中文字幕| 中文天堂在线官网| 午夜精品国产一区二区电影 | 午夜亚洲福利在线播放| 熟女电影av网| 久久人人爽人人片av| 久久这里只有精品中国| 成人午夜精彩视频在线观看| 午夜福利视频精品| 天天躁夜夜躁狠狠久久av| 中文字幕制服av| 国产成人精品婷婷| 日韩,欧美,国产一区二区三区| 亚洲va在线va天堂va国产| 亚洲精品久久久久久婷婷小说| 狠狠精品人妻久久久久久综合| 日韩三级伦理在线观看| 国产成人一区二区在线| 男插女下体视频免费在线播放| 亚洲精品乱码久久久久久按摩| 久久久久久久国产电影| 男女边吃奶边做爰视频| 女人久久www免费人成看片| 在线免费十八禁| 国产一区二区在线观看日韩| 久久精品熟女亚洲av麻豆精品 | 免费人成在线观看视频色| 国产高清有码在线观看视频| 边亲边吃奶的免费视频| 久久久久久久国产电影| 有码 亚洲区| 久久鲁丝午夜福利片| 精品一区二区三区视频在线| 看非洲黑人一级黄片| 免费观看精品视频网站| 国产av国产精品国产| 国产成人精品久久久久久| 国产69精品久久久久777片| 99久久中文字幕三级久久日本| 亚洲精华国产精华液的使用体验| 一区二区三区免费毛片| 只有这里有精品99| 亚洲人成网站高清观看| 最近中文字幕高清免费大全6| 亚洲av日韩在线播放| 晚上一个人看的免费电影| 亚洲熟妇中文字幕五十中出| 内射极品少妇av片p| 久久国内精品自在自线图片| 国产黄频视频在线观看| 国产在线一区二区三区精| 日韩强制内射视频| 国产精品久久久久久久电影| 一级av片app| 国产av不卡久久| 精品酒店卫生间| 最近2019中文字幕mv第一页| 亚洲国产精品专区欧美| 亚洲性久久影院| 久久久色成人| 3wmmmm亚洲av在线观看| 波多野结衣巨乳人妻| 日韩欧美精品免费久久| 欧美日韩在线观看h| 欧美97在线视频| 欧美zozozo另类| 欧美精品一区二区大全| av线在线观看网站| 18禁在线播放成人免费| 三级国产精品欧美在线观看| 午夜精品国产一区二区电影 | 在线观看免费高清a一片| 亚洲精品乱码久久久久久按摩| 尤物成人国产欧美一区二区三区| 免费大片18禁| 一区二区三区高清视频在线| 99久久人妻综合| 一级二级三级毛片免费看| 国产成人精品久久久久久| 亚洲av电影不卡..在线观看| a级毛片免费高清观看在线播放| 国产老妇女一区| 亚洲精品中文字幕在线视频 | 一级爰片在线观看| 一夜夜www| 精品人妻偷拍中文字幕| 中文在线观看免费www的网站| 国产黄色小视频在线观看| 国产爱豆传媒在线观看| 精品人妻一区二区三区麻豆| 观看美女的网站| 国产爱豆传媒在线观看| 国产在视频线在精品| 嫩草影院新地址| 91午夜精品亚洲一区二区三区| 国产综合懂色| 国产高清有码在线观看视频| 男人舔奶头视频| 国产成人a区在线观看| 久久韩国三级中文字幕| 国产三级在线视频| 欧美潮喷喷水| 插阴视频在线观看视频| 3wmmmm亚洲av在线观看| 免费看不卡的av| 男人和女人高潮做爰伦理| 热99在线观看视频| 欧美+日韩+精品| 国产精品久久久久久精品电影小说 | 日日摸夜夜添夜夜添av毛片| 禁无遮挡网站| 久久精品国产自在天天线| 97超碰精品成人国产| 国产片特级美女逼逼视频| 亚洲精品日韩av片在线观看| 直男gayav资源| 久久久久久国产a免费观看|