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

    RegCM4對中國東部區(qū)域氣候模擬的輻射收支分析

    2016-11-01 08:31:54韓振宇王宇星聶羽
    大氣科學(xué)學(xué)報 2016年5期
    關(guān)鍵詞:云量長波短波

    韓振宇,王宇星,聶羽

    ?

    RegCM4對中國東部區(qū)域氣候模擬的輻射收支分析

    韓振宇①*,王宇星②,聶羽①

    ① 中國氣象局 國家氣候中心,北京 100081;

    ② 國家海洋局 海洋減災(zāi)中心,北京 100194

    2015-11-04收稿,2016-05-05接受

    中國氣象局氣候變化專項項目(CCSF201626;CCSF201509);國家自然科學(xué)基金資助項目(41405101);公益性行業(yè)(氣象)科研專項(GYHY201306019)

    利用衛(wèi)星和再分析數(shù)據(jù),評估了區(qū)域氣候模式RegCM4對中國東部地區(qū)輻射收支的基本模擬能力,重點關(guān)注地表凈短波(SNS)、地表凈長波(SNL)、大氣頂凈短波(TNS)、大氣頂凈長波(TNL)4個輻射分量。結(jié)果表明:1)短波輻射的誤差值在夏季較大,而長波輻射的誤差值在冬季較大。但各輻射分量模擬誤差的空間分布在冬、夏季都有較好的一致性。2)對于地表輻射通量,SNS表現(xiàn)為正偏差(向下凈短波偏多),在各分量中誤差最大,區(qū)域平均誤差值近50 W/m2;SNL表現(xiàn)為負(fù)偏差(向上凈長波偏多);對于大氣頂輻射通量,TNS和TNL分別表現(xiàn)為“北負(fù)南正”的誤差分布和整體正偏差。3)利用空間相關(guān)和散點線性回歸方法對4個輻射分量的模擬誤差進(jìn)行歸因分析,發(fā)現(xiàn)在云量、地表反照率、地表溫度三個直接影響因子中,云量模擬誤差的貢獻(xiàn)最大,中國東部地區(qū)云量模擬顯著偏少。

    輻射收支

    模擬評估

    區(qū)域氣候模式

    云量

    地表反照率

    大氣頂?shù)膬糨椛涫钦麄€氣候系統(tǒng)的主要能量來源,這種輻射能量能轉(zhuǎn)化為其他形式的能量重新分布,如動能、潛熱能和感熱能等,最終起驅(qū)動氣候系統(tǒng)各圈層內(nèi)的運動、維持地球生命活動等作用。太陽輻射的季節(jié)變化及其導(dǎo)致的海陸溫差變化是東亞季風(fēng)的重要驅(qū)動因子之一。19世紀(jì)70年代以來,利用衛(wèi)星數(shù)據(jù)來精確測定全球輻射收支的工作陸續(xù)開展,提供了大量有關(guān)大氣頂輻射收支的信息,這些數(shù)據(jù)成為了評估氣候模式的重要依據(jù)。但僅根據(jù)大氣頂?shù)妮椛涫罩Р⒉荒苋娣从炒髿獾妮椛涮卣?輻射在大氣和地表之間的分配也是非常重要的(Gleckler,2005;Stephens,2005)。因此大氣頂(Top Of Atmosphere,TOA)和地表(SurFaCe,SFC)的輻射收支模擬都是區(qū)域氣候模式(Regional Climate Model,RCM)性能評估的重要內(nèi)容(Jaeger et al.,2008;Kothe et al.,2011;Güttler et al.,2014;Pessacg et al.,2014),可為RCM在氣候預(yù)測、氣候變化等方面的研究提供基礎(chǔ),但東亞地區(qū)的相關(guān)研究還較少。

    臺站觀測能夠提供精確的SFC輻射收支數(shù)據(jù),常被用來進(jìn)行RCM等模式的評估(閻宇平等,2001;劉樹華等,2008;于濤,2010;宋耀明等,2014)。但有限的臺站觀測無法提供時空連續(xù)的長期觀測資料,同時也缺少TOA的輻射資料;因此再分析資料和衛(wèi)星資料也常用來進(jìn)行模擬評估(Li et al.,2009;Kothe et al.,2011;Güttler et al.,2014;Wang et al.,2014)。

    國際聯(lián)合區(qū)域氣候降尺度試驗CORDEX(COordinated Regional climate Downscaling EXperiment),旨在基于全球模式的多模式、多情景預(yù)估結(jié)果,利用RCM等降尺度方法獲得高分辨率區(qū)域氣候變化信息(Giorgi et al.,2009)。目前,CORDEX-EastAsia第二階段已在開展,RCM分辨率由50 km提高到25 km,能更好地適應(yīng)東亞復(fù)雜的下墊面特征。本文的研究目的是,以再分析和衛(wèi)星資料的輻射數(shù)據(jù)為參考,對區(qū)域氣候模式RegCM4的CORDEX試驗在中國東部的模擬結(jié)果進(jìn)行初步評估。

    1 模式、試驗和資料

    1.1模式簡介及試驗設(shè)計

    國際理論物理中心研制的RegCM系列區(qū)域氣候模式有著廣泛的應(yīng)用,該系列模式對東亞和中國地區(qū)氣候有較好的模擬能力(Gao et al.,2003,2008;李巧萍和丁一匯,2004;鄒靖和謝正輝,2012;Zou and Zhou,2013)。本研究使用RegCM4.4版本,該模式有多種物理參數(shù)化方案選擇(Giorgi et al.,2012)。經(jīng)過在再分析資料驅(qū)動下的大量試驗的對比分析(王美麗,2015;Gao et al.,2016),本研究選擇了對中國氣候有較好模擬效果的參數(shù)化組合:輻射采用CCM3方案,行星邊界層使用Holtslag方案,大尺度降水采用SUBEX方案,積云對流選擇Emanuel方案,陸面使用CLM3.5方案。試驗使用的土地覆蓋資料在中國區(qū)域內(nèi)基于中國1:100萬植被圖得到(韓振宇等,2015)。該版本已經(jīng)被應(yīng)用在短期氣候預(yù)測(張冬峰等,2015)和CORDEX-East Asia第二階段的系列試驗。

    模擬區(qū)域采用CORDEX-EastAsia第二階段的推薦區(qū)域,覆蓋整個中國大陸及周邊地區(qū)(http://www.cordex.org/index.php?option=com_content&view=article&id=88&Itemid=625)。模式垂直方向分18層,水平分辨率取25 km。模式的初始和側(cè)邊界場數(shù)據(jù)、海溫數(shù)據(jù)都采用ERA-Interim再分析資料。試驗積分從1990年1月1日到2010年12月31日(世界時,下同),其中第1年作為模式初始化時段。

    1.2資料

    衛(wèi)星資料和再分析資料中的輻射數(shù)據(jù)都存在一定的不確定性,因此本文選用兩套資料,以考察觀測資料不確定性對評估結(jié)果的影響。

    1)CERES(Clouds and the Earth’s Radiant Energy System)EBAF逐月衛(wèi)星反演輻射資料(簡稱CERES),水平分辨率為1°×1°,時間跨度為2000年3月—2014年3月(Loeb et al.,2009;Kato et al.,2013)。CERES資料因其較高的精度已被廣泛用于氣候模式評估(吳春強(qiáng)和周天軍,2011)。

    2)ERA-Interim再分析資料(簡稱ERAI)中的輻射通量、地表溫度、地表反照率、云量,水平分辨率為1.5°×1.5°,時間跨度為1979年1月—2014年5月(Dee et al.,2011)。受觀測資料時段的限制,本文選取2001—2010年的平均值作為氣候態(tài)進(jìn)行評估分析。

    文中將不同分辨率的觀測和模式結(jié)果都水平插值到1°×1°分辨率,選取中國東部地區(qū)(108~123°E,20~40°N)為東亞季風(fēng)區(qū)的主要陸地區(qū)域進(jìn)行分析,文中的區(qū)域平均等統(tǒng)計分析全部基于該區(qū)域。需要注意,文中所有能量通量都是以向下為正。

    2 輻射收支的模擬評估

    2.1地表凈短波輻射(SNS)

    相比CERES和ERAI,模式模擬的區(qū)域平均SNS都偏高,其中以ERAI為參考得到的誤差相對較低(圖1)。在觀測中,冬季云量較少,地表短波輻射以緯向分布為主,由南部沿海向北逐漸減小,數(shù)值分布在50~150 W/m2;夏季受夏季風(fēng)相關(guān)云系的影響,江淮地區(qū)有明顯的低值中心,數(shù)值在150 W/m2附近,而華北地區(qū)的高值能超過200 W/m2(圖略)。

    圖1 相對于觀測基準(zhǔn)(ERAI+CERES)/2,CERES、ERAI和RegCM4中輻射分量SNS、TNS、SNL、TNL的區(qū)域平均偏差(單位:W/m2)Fig.1 Biases of CERES,ERA-Interim and RegCM4 relative to (ERA-Interim+CERES)/2 for SNS,TNS,SNL and TNL(for the whole investigation area and time period;units:W·m-2)

    圖2 相對于CERES(a,c,e,g)和ERAI(b,d,f,h),RegCM4模擬SNS和SNL的誤差(單位:W/m2)a,b.冬季SNS;c,d.夏季SNS;e,f.冬季SNL;g,h.夏季SNLFig.2 Biases of RegCM4 relative to (a,c,e,g)CERES and (b,d,f,h)ERA-Interim for SNS and SNL(units:W·m-2):(a,b)SNS in winter;(c,d)SNS in summer;(e,f)SNL in winter;(g,h)SNL in summer

    圖2a、c分別給出了模式多年平均SNS相對于CERES的差值,從空間分布圖中可以看出,冬季該區(qū)域以模擬偏多為主,尤其在30°N以南的地區(qū),SNS模擬得明顯偏多,中心偏高可超過80 W/m2;夏季模擬誤差明顯大于冬季,最高超過100 W/m2,但空間分布類似。如果以ERAI為參考,得到的誤差普遍較低,但空間分布一致。

    2.2地表凈長波輻射(SNL)

    相比CERES和ERAI,模式模擬的區(qū)域平均SNL都偏低,即模式模擬的地表長波熱量釋放偏大,其中以ERAI為參考得到的誤差相對較低(圖1)。在觀測中,冬季SNL的分布與云量有關(guān),-50 W/m2的等值線大致位于32°N,其以南區(qū)域的云量大都超過五成。湖南及以西的云量超過六成,對應(yīng)SNL的高值中心,可達(dá)-20 W/m2;夏季SNL的分布也與云量有很好的對應(yīng)關(guān)系,除華北個別區(qū)域外,多數(shù)區(qū)域的云量都超過五成,SNL大于-50 W/m2(圖略)。

    圖2e、g分別給出了模式多年平均SNL相對于CERES的差值,從空間分布圖中可以看出,冬季該區(qū)域以模擬偏少為主,尤其在南部沿海地區(qū),SNL模擬得明顯偏少,中心偏低可超過-80 W/m2;夏季,模擬誤差明顯低于冬季,但空間分布類似,誤差值大都分布在-5~-15 W/m2。如果以ERAI為參考,得到的誤差普遍較低,但空間分布一致。

    2.3大氣頂凈短波輻射(TNS)

    相比CERES和ERAI,模式模擬的區(qū)域平均TNS都偏低,其中以ERAI為參考得到的誤差相對較高(圖1)。在觀測中,冬、夏季,TNS的空間分布與SNS的相似,數(shù)值分別分布在120~220 W/m2和270~350 W/m2(圖略)。圖3a、c分別給出了模式多年平均TNS相對于CERES的差值,從空間分布圖中可以看出,冬季,以30°N為界,模擬誤差空間分布為“北負(fù)南正”,誤差值遠(yuǎn)低于SNS,分布在-30~40 W/m2。夏季的空間分布與冬季相似,北部負(fù)偏差增強(qiáng),而南部正偏差減弱。如果以ERAI為參考,得到的誤差以負(fù)偏差為主,但空間和季節(jié)分布變化不大。

    圖3 相對于CERES(a,c,e,g)和ERAI(b,d,f,h),RegCM4模擬TNS和TNL的誤差(單位:W/m2)a,b.冬季TNS;c,d.夏季TNS;e,f.冬季TNL;g,h.夏季TNLFig.3 Biases of RegCM4 relative to (a,c,e,g)CERES and (b,d,f,h)ERA-Interim for TNS and TNL(units:W·m-2):(a,b)TNS in winter;(c,d)TNS in summer;(e,f)TNL in winter;(g,h)TNL in summer

    2.4大氣頂凈長波輻射(TNL)

    相比CERES和ERAI,模式模擬的區(qū)域平均TNL都偏高,即模式模擬的大氣頂長波熱量釋放偏小,其中以ERAI為參考得到的誤差相對較高(圖1)。在觀測中,冬季云量較少,TNL以緯向分布為主,由南向北逐漸增大,數(shù)值分布在-270~-210 W/m2;夏季,受夏季風(fēng)相關(guān)云系的影響,江淮和華南地區(qū)有明顯的高值中心,數(shù)值大于-230 W/m2,而華北地區(qū)的低值超過-260 W/m2(圖略)。

    圖3e、g分別給出了模式多年平均TNL相對于CERES的差值,從空間分布圖中可以看出,冬季該區(qū)域以正偏差為主,尤其在30°N以北的地區(qū),TNL模擬誤差較大,多數(shù)區(qū)域超過20 W/m2;夏季的模擬誤差總體小于冬季,且在南部沿海有較小的負(fù)偏差,但是北部部分區(qū)域的誤差較大,最高超過40 W/m2。如果以ERAI為參考,得到的誤差普遍較高,但空間分布一致。

    總體來看,分別以CERES和ERAI為參考,得到的模式誤差的季節(jié)和空間分布都較為一致,所以能基本排除觀測資料的不確定性對模式評估結(jié)果的影響。模式模擬的SFC輻射通量更接近ERAI,而TOA輻射通量更接近CERES,其中SNS的模擬誤差最大。

    3 影響因子分析

    模式顯著高估了SFC的凈入射短波輻射和凈出射長波輻射,而在TOA是相反的誤差。對于總的輻射收支,這些分量的模擬誤差會部分抵消。模式模擬SNS偏高通常有兩種可能的原因,一是云量或者其他大氣短波吸收量的模擬偏低,導(dǎo)致模擬的向下短波輻射偏大;二是由地表反照率的模擬偏低所引起。SNL的模擬誤差可能來自云量或者地表溫度模擬誤差引起的向外長波高估。TNS、TNL模擬誤差的影響因子分別與SNS、SNL類似,但具體的物理過程有所區(qū)別。將在本節(jié)詳細(xì)分析誤差的來源(如無特別說明,以下的模擬誤差都是基于CERES)。

    圖4 相對于CERES,RegCM4模擬總云量(a—c)和地表反照率(d—f)的誤差  a,d.年平均;b,e.冬季;c,f.夏季Fig.4 Biases of RegCM4 relative to CERES for (a—c)total cloud fraction and (d—f)surface albedo:(a,d)annual mean;(b,e)winter mean;(c,f)summer mean

    推測云量(CFR)、地表反照率(ALB)、地表溫度(TS)是直接影響輻射收支各分量模擬的三個主要因素。如果CFR和ALB模擬偏低,可能會起到增加SNS和TNS的作用;TS模擬偏低,可能會起到增加SNL和TNL的作用。對比這些要素模擬誤差的空間分布,可以發(fā)現(xiàn)一些定性的對應(yīng)關(guān)系:CFR的年平均模擬誤差以偏少為主,在30°N以南的部分地區(qū)偏少超過二成,而在北部部分地區(qū)存在正偏差(圖4a),正偏差在冬季更為顯著(圖4b)。這樣的空間分布與SNS和TNS模擬誤差有很好的空間對應(yīng)關(guān)系,且在冬季相關(guān)更強(qiáng),其相關(guān)系數(shù)達(dá)-0.87和-0.86(表1);而ALB模擬以偏低為主,誤差值的空間分布較為均勻,沒有明顯的高低值中心(圖4d),且季節(jié)差異較小(圖4e—f)。ALB與SNS或者TNS模擬誤差之間的空間相關(guān)較弱,且不是推測的負(fù)相關(guān)(表1),因此判斷ALB模擬誤差對短波輻射模擬的影響較小;出乎預(yù)料的是,TS和SNL模擬誤差的相關(guān)在冬、夏季符號相反,且TS和TNL模擬誤差之間并沒有顯著的相關(guān)(表1),因此推測TS模擬誤差對長波輻射模擬的影響較小。

    表1輻射分量模擬誤差與CFR、ALB、TS模擬誤差的空間相關(guān)系數(shù)

    Table 1The spatial correlation coefficients between biases of radiation components and that of CFR,ALB,and TS

    CFRALB年冬季夏季年冬季夏季SNS-080-087-049023010008TNS-081-086-056040019026CFRTS年冬季夏季年冬季夏季SNL080077062005044-049TNL072065064-018001-029

    下面進(jìn)一步給出CFR、ALB、TS這三個要素綜合影響各輻射分量模擬的定量分析。圖5a給出的是SNS模擬誤差對CFR和ALB誤差(dCFR和dALB)的依賴度。當(dāng)ALB的誤差為零時,隨著CFR誤差值由正到負(fù),SNS的模擬誤差會從-14 W/m2增加到78 W/m2。從圖5c可以看出,多數(shù)樣本點(約84%)都位于坐標(biāo)系的第三象限,dCFR和dALB都是負(fù)值,即模擬結(jié)果中云量偏少、反照率偏弱,對應(yīng)圖5a中SNS模擬偏高。并且SNS誤差僅在縱軸方向(dCFR變化的方向)上的變化較為顯著,再次說明ALB與SNS模擬誤差之間的相關(guān)較弱。對TNS的分析也是類似的(圖5b)。

    圖5 模式中SNS(a)和TNS(b)模擬誤差(單位:W/m2)與ALB和CFR模擬誤差的對應(yīng)關(guān)系以及ALB和CFR模擬誤差的聯(lián)合分布(c,用樣本點數(shù)表示)Fig.5 Monthly means of model biases in (a)SNS and (b)TNS (units:W·m-2) against biases in ALB and CFR and (c)the number of occurrence in each ALB-CFR class

    圖6a給出的是SNL模擬誤差對CFR和TS誤差(dCFR和dTS)的依賴度。當(dāng)TS的誤差為零時,隨著CFR誤差值由正到負(fù),SNL的模擬誤差會從15 W/m2降低到-42 W/m2。從圖6c可以看出,近半數(shù)樣本點(約46%)都位于坐標(biāo)系第二、三象限的零值附近。dCFR是負(fù)值,說明模擬結(jié)果中云量偏少,對應(yīng)圖6a中SNL模擬偏低,地表長波熱量釋放偏大。并且SNL誤差在橫軸方向(dTS變化的方向)上幾乎沒有變化,說明TS與SNL模擬誤差之間的相關(guān)沒有CFR強(qiáng)。對于TNL(圖6b),其模擬誤差都是正值,但隨著CFR誤差值由正到負(fù)變化時,正偏差也在逐漸變小。這是因為,CFR偏小時大氣頂向外長波熱量釋放較大,TNL數(shù)值偏小。同樣,TS與TNL模擬誤差之間的相關(guān)較弱。

    假設(shè)不同誤差之間的關(guān)系是線性的,線性回歸的決定系數(shù)R2可以用來衡量誤差貢獻(xiàn)率。基于CERES和ERAI的計算都顯示,dCFR對各個輻射分量模擬誤差的貢獻(xiàn)率都遠(yuǎn)大于dALB或者dTS,且對各個輻射分量的影響也是相當(dāng)?shù)?圖略)。

    4 結(jié)論和討論

    利用衛(wèi)星和再分析數(shù)據(jù),初步評估了中國東部區(qū)域模式RegCM4模擬的輻射收支過程,主要結(jié)論有:

    1)在中國東部區(qū)域,RegCM4顯著高估了區(qū)域平均的SFC凈入射短波輻射和凈出射長波輻射,而在TOA是相反的誤差。這一結(jié)論不依賴于觀測資料的選擇,且誤差的空間分布隨季節(jié)變化不大。短波輻射的誤差值在夏季較大,季節(jié)差異超過20 W/m2;而長波輻射的誤差值在冬季較大,季節(jié)差異在10 W/m2附近。輻射分量SNS的模擬誤差最大,區(qū)域年平均誤差值將近50 W/m2。

    2)選擇CFR、ALB、TS三個直接影響因子,對四個輻射分量的模擬誤差進(jìn)行線性歸因分析。中國東部地區(qū)CFR模擬顯著偏少,年平均模擬值在30°N以南的部分地區(qū)偏少超過二成。其空間和散點分布都與四個輻射分量的誤差有很好的相關(guān)關(guān)系??臻g相關(guān)、散點線性回歸的分析結(jié)果都顯示,三個因子中CFR模擬誤差的貢獻(xiàn)最大,貢獻(xiàn)率大都超過25%,空間相關(guān)系數(shù)也超過0.5,高于ALB和TS低于10%的貢獻(xiàn)率和較低的空間相關(guān)系數(shù)。

    圖6 模式中SNL(a)和TNL(b)模擬誤差(單位:W/m2)與TS和CFR模擬誤差的對應(yīng)關(guān)系及TS和CFR模擬誤差的聯(lián)合分布(c,用樣本點數(shù)表示)Fig.6 Monthly means of model biases in (a)SNL and (b)TNL(units:W·m-2) against biases in TS and CFR and (c)the number of occurrence in each TS-CFR class

    RegCM4中SNS的模擬誤差最為突出,以CERES為參考,RegCM4模擬SNS偏高,與眾多CMIP3/5模式的模擬結(jié)果相似(Li et al.,2013)。ERAI作為再分析系統(tǒng)的直接模擬輸出也表現(xiàn)出類似的差異(圖2a—d),且RegCM4模擬的云量也與ERAI更加接近(圖略),可能是RegCM4與這些全球氣候模式及再分析系統(tǒng)存在相似的缺陷。例如,Li et al.(2013)指出全球模式在計算輻射傳輸時普遍未考慮云中雨(雪)水或?qū)α髟圃扑妮椛鋸?qiáng)迫作用,RegCM4是否有類似的影響還需進(jìn)一步分析。目前,RegCM4已引入新的大尺度降水方案,是包括五種水物質(zhì)的混合相云方案(Nogherotto et al.,2012;Nogherotto,2015),其對云輻射效應(yīng)的描述更加準(zhǔn)確,未來需要進(jìn)一步分析新的方案對RegCM4輻射收支模擬的影響。此外,本文僅對比分析了對輻射通量模擬誤差產(chǎn)生直接影響的三個物理因子,未來對RegCM4模式的進(jìn)一步改進(jìn)時,還需要深入分析環(huán)流等偏差或者物理過程參數(shù)化方案對云量、地表反照率等模擬的影響。

    References)

    Dee D,Uppala S,Simmons A,et al.,2011.The ERA-Interim reanalysis:Configuration and performance of the data assimilation system[J].Quart J Roy Meteor Soc,137(656):553-597.

    Gao X,Luo Y,Lin W,et al.,2003.Simulation of effects of land use change on climate in China by a regional climate model[J].Adv Atmos Sci,20(4):583-592.

    Gao X,Shi Y,Song R,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.

    Gao X J,Shi Y,Giorgi F,2016.Comparison of convective parameterizations in RegCM4 experiments over China with CLM as the land surface model[J].Atmospheric & Oceanic Science Letters,doi:10.1080/16742834.2016.1172938.

    Giorgi F,Jones C,Asrar G R,2009.Addressing climate information needs at the regional level:The CORDEX framework[J].WMO Bulletin,58(3):175-183.

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

    Gleckler P J,2005.Surface energy balance errors in AGCMs:Implications for ocean-atmosphere model coupling[J].Geophys Res Lett,32(15):L15708.

    韓振宇,高學(xué)杰,石英,等,2015.中國高精度土地覆蓋數(shù)據(jù)在RegCM4/CLM模式中的引入及其對區(qū)域氣候模擬影響的分析[J].冰川凍土,37(4):857-866.Han Z Y,Gao X J,Shi Y,et al.,2015.Development of Chinese high resolution land cover for the RegCM4/CLM and its impact on regional climate simulation[J].Journal of Glaciology and Geocryology,37(4):857-866.(in Chinese).

    Jaeger E B,Anders I,Lüthi D,et al.,2008.Analysis of ERA40-driven CLM simulations for Europe[J].Meteor Z,17(4):349-367.

    Kato S,Loeb N G,Rose F G,et al.,2013.Surface irradiances consistent with CERES-derived top-of-atmosphere shortwave and longwave irradiances[J].J Climate,26(9):2719-2740.

    Kothe S,Dobler A,Beck A,et al.,2011.The radiation budget in a regional climate model[J].Climate Dyn,36(5/6):1023-1036.

    Li J D,Liu Y M,Wu G X,2009.Cloud radiative forcing in Asian monsoon region simulated by IPCC AR4 AMIP models[J].Adv Atmos Sci,26(5):923-939.

    Li J L,Waliser D E,Stephens G,et al.,2013.Characterizing and understanding radiation budget biases in CMIP3/CMIP5 GCMs,contemporary GCM,and reanalysis[J].J Geophys Res:Atmospheres,118(15):8166-8184.

    李巧萍,丁一匯,2004.區(qū)域氣候模式對東亞季風(fēng)和中國降水的多年模擬與性能檢驗[J].氣象學(xué)報,62(2):140-153.Li Q P,Ding Y H,2004.Multi-year simulation of the east Asian monsoon and precipitation in China using a regional climate model and evaluation[J].Acta Meteorologica Sinica,62(2):140-153.(in Chinese).

    劉樹華,蔣浩宇,胡非,等,2008.利用區(qū)域尺度氣象模式模擬黑河地區(qū)地表能量通量的研究[J].大氣科學(xué),32(6):1392-1400.Liu S H,Jiang H Y,Hu F,et al.,2008.A study of surface energy fluxes in Heihe region simulated with a mesoscale atmospheric model[J].Chin J Atmos Sci,32(6):1392-1400.(in Chinese).

    Loeb N G,Wielicki B A,Doelling D R,et al.,2009.Toward optimal closure of the Earth's top-of-atmosphere radiation budget[J].J Climate,22(3):748-766.

    Nogherotto R,2015.A numerical framework for multiple phase cloud microphysics in regional and global atmospheric models[D].Trieste:University of Trieste.

    Nogherotto R,Tompkins A M,Giorgi F,et al.,2012.The new microphysic cloud scheme implemented in RegCM4[C]//EGU General Assembly Conference.Abstracts,8017.

    Pessacg N L,Solman S A,Samuelsson P,et al.,2014.The surface radiation budget over South America in a set of regional climate models from the CLARIS-LPB project[J].Climate Dyn,43(5-6):1221-1239.

    宋耀明,范軼,馬天嬌,2014.陸面過程模式CLM4.5在半干旱區(qū)退化草原站的模擬性能評估[J].大氣科學(xué)學(xué)報,37(6):794-803.Song Y M,Fan Y,Ma T J,2014.Evaluation of simulation performance of land surface model NCAR_CLM4.5 at a degrated glassland station in semi-arid area[J].Trans Atmos Sci,37(6):794-803.(in Chinese).

    Stephens G L,2005.Cloud feedbacks in the climate system:A critical review[J].J Climate,18(2):237-273.

    Wang F,Yang S,Wu T,2014.Radiation budget biases in AMIP5 models over the East Asian monsoon region[J].J Geophys Res:Atmospheres,119(23):13,400-413,426.

    王美麗,2015.RegCM4.4區(qū)域氣候模式對中國當(dāng)代氣候及極端事件的高分辨率數(shù)值模擬[D].北京:中國氣象科學(xué)研究院.Wang M L,2015.A high resolution simulation of present climate and extreme events over East Asia by RegCM4.4[D].Beijing:Chinese Academy of Meteorological Sciences.(in Chinese).

    吳春強(qiáng),周天軍,2011.CFMIP大氣環(huán)流模式模擬的東亞云輻射強(qiáng)迫特征[J].氣象學(xué)報,69(3):381-399.Wu C Q,Zhou T J,2011.Characteristics of cloud radiative forcings over East Asia as simulated by the AGCMs in the CFMIP[J].Acta Meteorologica Sinica,69(3):381-399.(in Chinese).

    閻宇平,王介民,Menenti M,et al.,2001.黑河實驗區(qū)非均勻地表能量通量的數(shù)值模擬[J].高原氣象,20(2):132-139.Yan Y P,Wang J M,Menenti M,et al.,2001.Numerical simulation of surface energy fluxes over heterogeneous land surfaces in HEIFE area[J].Plateau Meteor,20(2):132-139.(in Chinese).

    于濤,2010.綠洲—沙漠系統(tǒng)地表輻射收支的模擬研究[J].中國沙漠,30(3):686-690.Yu T,2010.Simulation of surface radiation balance over the oasis-desert system[J].Journal of Desert Research,30(3):686-690.(in Chinese).

    張冬峰,高學(xué)杰,馬潔華,2015.CCSM4.0模式及其驅(qū)動下RegCM4.4模式對中國夏季氣候的回報分析[J].氣候與環(huán)境研究,20(3):307-318.Zhang D F,Gao X J,Ma J H,2015.Analysis of summer climate over China from hindcast experiments by CCSM4.0 and RegCM4.4 models[J].Climatic Environ Res,20(3):307-318.(in Chinese).

    鄒靖,謝正輝,2012.RegCM4 中陸面過程參數(shù)化方案對東亞區(qū)域氣候模擬的影響[J].氣象學(xué)報,70(6):1312-1326.Zou J,Xie Z H,2012.The effects of the land-surface process parameterization of the RegCM4 on climate simulation in East Asia[J].Acta Meteorologica Sinica,70(6):1312-1326.(in Chinese).

    Zou L,Zhou T,2013.Near future(2016—40) summer precipitation changes over China as projected by a regional climate model(RCM) under the RCP8.5 emissions scenario:Comparison between RCM downscaling and the driving GCM[J].Adv Atmos Sci,30(3):806-818.

    The net radiation at the top of the atmosphere(TOA) is one of the most important driving forces of the Earth’s climate system,while the radiation budget of the climate system directly influences atmospheric and oceanic motion,as well as the atmospheric water cycle.In this paper,the performance of the regional climate model RegCM4 in simulating the radiation budget over eastern China is evaluated.A long-term simulation at a grid spacing of 25 km driven by ERA-Interim data for the period 1990—2010 is conducted.The East Asia domain of phase II of the international COordinated Regional climate Downscaling EXperiment(CORDEX) is used as the model domain.The simulation is compared with a satellite dataset and reanalysis over eastern China during 2001—10 in terms of surface net shortwave(SNS),surface net longwave(SNL),TOA net shortwave(TNS),and TOA net longwave(TNL) radiation.Results show that biases in shortwave fluxes are larger in summer,and biases in longwave fluxes are larger in winter.However,spatial distributions of the biases in summer and winter are similar.For surface fluxes,biases of SNS are positive(overestimation of net downward shortwave flux),with a regional average of nearly 50 W/m2,which is the maximum among these radiation flux biases.Biases of SNL are negative(overestimation of net upward longwave flux).For TOA fluxes,biases of TNS are negative in the north but positive in the south,and biases of TNL are overall positive.Attribution analysis on the biases in radiation fluxes is conducted by two methods—spatial correlation and linear regression.The biases in the cloud fraction,surface albedo,and surface temperature are considered.We find that biases in the radiation fluxes can be explained,to a large degree,by significant negative biases in cloud fraction.

    radiation budget;model evaluation;regional climate model;cloud fraction;surface albedo

    (責(zé)任編輯:孫寧)

    The radiation budget in a regional climate simulation by RegCM4 for eastern China

    HAN Zhenyu1,WANG Yuxing2,NIE Yu1

    1NationalClimateCenter,ChinaMeteorologicalAdministration,Beijing100081,China;2NationalMarineHazardMitigationService,Beijing100194,China

    10.13878/j.cnki.dqkxxb.20151104001

    引用格式:韓振宇,王宇星,聶羽,2016.RegCM4對中國東部區(qū)域氣候模擬的輻射收支分析[J].大氣科學(xué)學(xué)報,39(5):683-691.

    Han Z Y,Wang Y X,Nie Y,2016.The radiation budget in a regional climate simulation by RegCM4 for eastern China[J].Trans Atmos Sci,39(5):683-691.doi:10.13878/j.cnki.dqkxxb.20151104001.(in Chinese).

    *聯(lián)系人,E-mail:hanzy@cma.gov.cn

    猜你喜歡
    云量長波短波
    贛州地區(qū)云量變化特征及其與降水的關(guān)系
    ASC200型地基雙波段全天空云量自動觀測儀云量比對分析
    氣象科技(2022年2期)2022-04-28 09:35:46
    廣義對稱正則長波方程的孤波解和周期波解及它們與Hamilton能量的關(guān)系
    1971—2010年虎林市云量與氣溫、降水的年際變化特征分析
    樂海短波
    人民音樂(2016年1期)2016-11-07 10:02:42
    工運短波
    工運短波
    綠野短波
    基于構(gòu)架點頭角速度的軌道垂向長波不平順在線檢測
    掃描型長波紅外連續(xù)變焦光學(xué)系統(tǒng)
    一级毛片久久久久久久久女| 国产在线视频一区二区| 欧美3d第一页| 国产黄片视频在线免费观看| 国产精品免费大片| 日韩伦理黄色片| 成人毛片a级毛片在线播放| 久久久精品免费免费高清| 建设人人有责人人尽责人人享有的| 国产伦精品一区二区三区视频9| 欧美日韩视频精品一区| 两个人的视频大全免费| 亚洲国产日韩一区二区| 精品酒店卫生间| 少妇被粗大的猛进出69影院 | 欧美精品高潮呻吟av久久| 亚洲精品第二区| videossex国产| 精品久久久精品久久久| 国产片特级美女逼逼视频| 简卡轻食公司| 国产精品久久久久久精品古装| 国产精品蜜桃在线观看| videos熟女内射| av在线app专区| 啦啦啦啦在线视频资源| 欧美日韩av久久| av国产久精品久网站免费入址| 精品国产国语对白av| 精品少妇内射三级| 国产免费一区二区三区四区乱码| 性色avwww在线观看| 卡戴珊不雅视频在线播放| 久久99热这里只频精品6学生| 我要看黄色一级片免费的| av国产久精品久网站免费入址| 一本久久精品| 久久韩国三级中文字幕| 亚洲av成人精品一二三区| 国产欧美日韩综合在线一区二区 | 久久人人爽av亚洲精品天堂| 美女xxoo啪啪120秒动态图| 欧美日韩国产mv在线观看视频| 亚洲成人手机| 一本久久精品| 美女cb高潮喷水在线观看| 91在线精品国自产拍蜜月| 青春草亚洲视频在线观看| 亚洲国产av新网站| 大香蕉久久网| 三级国产精品欧美在线观看| 高清视频免费观看一区二区| 少妇精品久久久久久久| 免费人成在线观看视频色| 精品久久久噜噜| 黄色一级大片看看| 一个人看视频在线观看www免费| 午夜免费男女啪啪视频观看| 国产欧美日韩精品一区二区| 亚洲av在线观看美女高潮| 啦啦啦啦在线视频资源| 美女cb高潮喷水在线观看| 国产精品一区www在线观看| 2021少妇久久久久久久久久久| 午夜老司机福利剧场| 久久97久久精品| av在线播放精品| 99九九在线精品视频 | 欧美少妇被猛烈插入视频| 啦啦啦在线观看免费高清www| 丰满乱子伦码专区| 亚洲欧洲日产国产| 国产国拍精品亚洲av在线观看| 美女国产视频在线观看| 久久精品久久久久久噜噜老黄| 美女国产视频在线观看| 日韩精品有码人妻一区| 亚洲av成人精品一二三区| 国产伦精品一区二区三区视频9| 中文精品一卡2卡3卡4更新| 一本—道久久a久久精品蜜桃钙片| 国产精品三级大全| 日韩中文字幕视频在线看片| 自拍欧美九色日韩亚洲蝌蚪91 | videos熟女内射| 国产高清三级在线| 精品久久久噜噜| 丝袜喷水一区| 国产亚洲av片在线观看秒播厂| 成人综合一区亚洲| 成人综合一区亚洲| 十八禁网站网址无遮挡 | 少妇裸体淫交视频免费看高清| 日本av手机在线免费观看| 久久97久久精品| 欧美日韩视频高清一区二区三区二| 多毛熟女@视频| 在线观看免费高清a一片| 校园人妻丝袜中文字幕| 色视频在线一区二区三区| 久久韩国三级中文字幕| 欧美日韩国产mv在线观看视频| 性色avwww在线观看| 黄色视频在线播放观看不卡| 精品一区二区三区视频在线| 亚洲国产色片| 最近2019中文字幕mv第一页| 性色av一级| 妹子高潮喷水视频| 免费看日本二区| 精品亚洲成a人片在线观看| 国产爽快片一区二区三区| 97超碰精品成人国产| 大又大粗又爽又黄少妇毛片口| 亚洲人成网站在线观看播放| 人妻一区二区av| 亚洲电影在线观看av| 精品人妻偷拍中文字幕| 亚洲av二区三区四区| 久久久国产精品麻豆| 你懂的网址亚洲精品在线观看| 久久热精品热| 免费高清在线观看视频在线观看| 欧美亚洲 丝袜 人妻 在线| 欧美日韩在线观看h| 久久国产乱子免费精品| 国产一区二区三区av在线| 日日撸夜夜添| 自拍欧美九色日韩亚洲蝌蚪91 | videos熟女内射| 一级爰片在线观看| 精品视频人人做人人爽| 亚洲中文av在线| 久久女婷五月综合色啪小说| av.在线天堂| a级毛片在线看网站| 只有这里有精品99| 国产午夜精品久久久久久一区二区三区| 免费观看a级毛片全部| 亚洲怡红院男人天堂| 黑人高潮一二区| 亚洲欧洲日产国产| 国产亚洲最大av| 青青草视频在线视频观看| 国语对白做爰xxxⅹ性视频网站| 午夜91福利影院| av又黄又爽大尺度在线免费看| 亚洲一级一片aⅴ在线观看| 2018国产大陆天天弄谢| 日本猛色少妇xxxxx猛交久久| 久久综合国产亚洲精品| 美女内射精品一级片tv| 中文字幕精品免费在线观看视频 | av网站免费在线观看视频| 如日韩欧美国产精品一区二区三区 | 亚洲欧美清纯卡通| 亚洲四区av| 亚洲av不卡在线观看| 欧美精品高潮呻吟av久久| 午夜视频国产福利| 性高湖久久久久久久久免费观看| 国产精品久久久久久久电影| 永久网站在线| 欧美精品一区二区免费开放| 日日撸夜夜添| 国产亚洲午夜精品一区二区久久| 能在线免费看毛片的网站| 亚洲精品日本国产第一区| 国产精品一二三区在线看| 欧美成人午夜免费资源| 成人美女网站在线观看视频| a级片在线免费高清观看视频| 国产精品一区二区在线观看99| 国产精品久久久久久久电影| 国产精品国产三级专区第一集| 欧美日韩精品成人综合77777| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲最大av| 毛片一级片免费看久久久久| 久久国内精品自在自线图片| 亚洲精品中文字幕在线视频 | 成年人免费黄色播放视频 | 中国三级夫妇交换| 午夜久久久在线观看| 中文字幕av电影在线播放| 亚洲av免费高清在线观看| 日韩一区二区三区影片| 国产淫语在线视频| 不卡视频在线观看欧美| 国产黄色免费在线视频| 观看av在线不卡| 亚洲自偷自拍三级| av卡一久久| 热99国产精品久久久久久7| 韩国高清视频一区二区三区| 精品一区二区三卡| 国产在线视频一区二区| 91aial.com中文字幕在线观看| 成人特级av手机在线观看| 亚洲欧美清纯卡通| 亚洲精品第二区| 九九爱精品视频在线观看| 国产精品人妻久久久久久| 欧美最新免费一区二区三区| 这个男人来自地球电影免费观看 | 国产黄片美女视频| 亚洲精品视频女| 亚洲伊人久久精品综合| a级毛片免费高清观看在线播放| 美女福利国产在线| 欧美日韩视频精品一区| 这个男人来自地球电影免费观看 | 亚洲国产精品999| 久久99一区二区三区| 日本av免费视频播放| 国产男人的电影天堂91| 嫩草影院新地址| 欧美日韩av久久| 黑人巨大精品欧美一区二区蜜桃 | 观看美女的网站| 女性被躁到高潮视频| 九九爱精品视频在线观看| 欧美变态另类bdsm刘玥| 精品亚洲成国产av| 久久亚洲国产成人精品v| 久久久国产精品麻豆| 日日啪夜夜撸| 国产精品人妻久久久影院| 黄色毛片三级朝国网站 | 我的老师免费观看完整版| 亚洲美女视频黄频| 激情五月婷婷亚洲| 简卡轻食公司| 97在线人人人人妻| 欧美日韩一区二区视频在线观看视频在线| 精品久久久久久久久亚洲| 一级毛片 在线播放| 国产淫语在线视频| 青春草视频在线免费观看| 国产色婷婷99| 国产午夜精品一二区理论片| 日韩一区二区视频免费看| 亚洲精品一二三| 大话2 男鬼变身卡| 日韩不卡一区二区三区视频在线| 18禁在线播放成人免费| 亚洲精品国产av成人精品| 国产深夜福利视频在线观看| 国产精品国产av在线观看| 日韩欧美 国产精品| 婷婷色综合大香蕉| 天美传媒精品一区二区| 麻豆成人午夜福利视频| 国产色爽女视频免费观看| 久久精品国产亚洲av涩爱| 成人黄色视频免费在线看| 在线观看三级黄色| av不卡在线播放| 久久久久久久大尺度免费视频| 亚洲国产日韩一区二区| 久久热精品热| 亚洲国产欧美在线一区| 亚洲国产最新在线播放| 久久久久久久国产电影| 久久精品国产亚洲av涩爱| 国产黄频视频在线观看| 狂野欧美激情性xxxx在线观看| 一本一本综合久久| a级一级毛片免费在线观看| 高清午夜精品一区二区三区| 国产精品国产三级国产av玫瑰| 啦啦啦视频在线资源免费观看| 中文字幕制服av| 内射极品少妇av片p| 男人舔奶头视频| 九色成人免费人妻av| 国产精品久久久久成人av| 亚洲精品一区蜜桃| 成人免费观看视频高清| av免费观看日本| av又黄又爽大尺度在线免费看| 久久6这里有精品| 在线观看人妻少妇| 成人亚洲精品一区在线观看| 亚洲av成人精品一二三区| 久久99热这里只频精品6学生| 久久99热6这里只有精品| 91精品一卡2卡3卡4卡| 欧美精品人与动牲交sv欧美| 国产亚洲最大av| av女优亚洲男人天堂| 日韩一区二区视频免费看| 久久久亚洲精品成人影院| 黑人猛操日本美女一级片| 国产精品一区www在线观看| 日韩电影二区| 寂寞人妻少妇视频99o| 精品久久国产蜜桃| 日韩亚洲欧美综合| 国产淫语在线视频| 三级经典国产精品| 草草在线视频免费看| 午夜老司机福利剧场| 精品亚洲乱码少妇综合久久| 久久女婷五月综合色啪小说| 精品午夜福利在线看| 亚洲四区av| a级一级毛片免费在线观看| 亚洲激情五月婷婷啪啪| 国产亚洲91精品色在线| 老司机亚洲免费影院| 国产日韩欧美在线精品| 久久99蜜桃精品久久| 天美传媒精品一区二区| 男人狂女人下面高潮的视频| 国产精品99久久久久久久久| 街头女战士在线观看网站| 亚洲婷婷狠狠爱综合网| 99热这里只有精品一区| 高清午夜精品一区二区三区| 国产成人91sexporn| 亚洲av免费高清在线观看| 99热这里只有是精品50| 亚洲精品aⅴ在线观看| 啦啦啦视频在线资源免费观看| 91成人精品电影| 精品国产一区二区三区久久久樱花| 丰满人妻一区二区三区视频av| 我要看黄色一级片免费的| 秋霞伦理黄片| 久久久午夜欧美精品| av国产精品久久久久影院| 欧美 日韩 精品 国产| 黄色一级大片看看| 大码成人一级视频| 男人和女人高潮做爰伦理| 乱码一卡2卡4卡精品| 夜夜看夜夜爽夜夜摸| 午夜激情久久久久久久| 精品一区二区三区视频在线| 国产淫语在线视频| freevideosex欧美| 一级二级三级毛片免费看| 亚洲成人一二三区av| 最近中文字幕2019免费版| av免费在线看不卡| 成年人午夜在线观看视频| 午夜免费观看性视频| 3wmmmm亚洲av在线观看| 亚洲成人手机| tube8黄色片| 国产免费一级a男人的天堂| 高清av免费在线| 老司机亚洲免费影院| 日本欧美视频一区| 天堂中文最新版在线下载| 色哟哟·www| 国产亚洲一区二区精品| 女性生殖器流出的白浆| 97超视频在线观看视频| 老司机影院成人| 久久婷婷青草| 日日啪夜夜撸| 国精品久久久久久国模美| 涩涩av久久男人的天堂| 日韩av在线免费看完整版不卡| 成人漫画全彩无遮挡| 特大巨黑吊av在线直播| 日本免费在线观看一区| 多毛熟女@视频| 国产69精品久久久久777片| 十八禁网站网址无遮挡 | 在线免费观看不下载黄p国产| 国产在线视频一区二区| 精品久久久久久电影网| 国产片特级美女逼逼视频| 18+在线观看网站| 日本av免费视频播放| 热re99久久国产66热| 成人美女网站在线观看视频| 男人舔奶头视频| av国产精品久久久久影院| 精品99又大又爽又粗少妇毛片| 老司机影院毛片| 久久久国产一区二区| 午夜激情福利司机影院| 久久国产精品大桥未久av | 久久精品国产亚洲av天美| 成人免费观看视频高清| 我的女老师完整版在线观看| 高清欧美精品videossex| 国产精品麻豆人妻色哟哟久久| 国产精品一区二区性色av| 久久久久人妻精品一区果冻| 91精品国产九色| 成人亚洲欧美一区二区av| 国产国拍精品亚洲av在线观看| 亚洲国产欧美日韩在线播放 | 久久6这里有精品| 一级爰片在线观看| 免费高清在线观看视频在线观看| 好男人视频免费观看在线| 久久国产精品男人的天堂亚洲 | 久久人人爽人人片av| 欧美人与善性xxx| 最新中文字幕久久久久| 国产在视频线精品| 久久这里有精品视频免费| 国产午夜精品久久久久久一区二区三区| 毛片一级片免费看久久久久| 亚洲av在线观看美女高潮| 国产黄色免费在线视频| 黄色毛片三级朝国网站 | 亚洲综合色惰| 国产精品久久久久成人av| 99国产精品免费福利视频| 菩萨蛮人人尽说江南好唐韦庄| 国产成人aa在线观看| 成人亚洲精品一区在线观看| 伦理电影免费视频| 亚洲人成网站在线播| 成人毛片60女人毛片免费| av天堂久久9| 乱码一卡2卡4卡精品| 91aial.com中文字幕在线观看| 精品一区二区三区视频在线| 日韩成人伦理影院| 一级片'在线观看视频| 成人二区视频| 午夜福利网站1000一区二区三区| 亚洲色图综合在线观看| 欧美一级a爱片免费观看看| 国产熟女午夜一区二区三区 | 春色校园在线视频观看| 国产成人午夜福利电影在线观看| 精品久久久久久久久亚洲| 最近的中文字幕免费完整| 久久久久久久久久久免费av| 色网站视频免费| 男人舔奶头视频| 黄色毛片三级朝国网站 | 日本猛色少妇xxxxx猛交久久| 久久久a久久爽久久v久久| 自拍偷自拍亚洲精品老妇| videossex国产| 国产成人一区二区在线| 国产精品人妻久久久影院| 中文字幕亚洲精品专区| 日韩人妻高清精品专区| 毛片一级片免费看久久久久| 国产又色又爽无遮挡免| 日本猛色少妇xxxxx猛交久久| 国产成人精品一,二区| 国产成人免费观看mmmm| 丝瓜视频免费看黄片| 两个人的视频大全免费| 免费观看性生交大片5| 欧美高清成人免费视频www| 亚洲怡红院男人天堂| 超碰97精品在线观看| 日本黄大片高清| 啦啦啦视频在线资源免费观看| 久热这里只有精品99| 国产成人91sexporn| 涩涩av久久男人的天堂| 日韩中文字幕视频在线看片| 女的被弄到高潮叫床怎么办| 五月伊人婷婷丁香| 亚洲av不卡在线观看| 搡女人真爽免费视频火全软件| 欧美日韩亚洲高清精品| a级毛片在线看网站| 国产乱来视频区| 国产精品福利在线免费观看| 晚上一个人看的免费电影| 国产精品无大码| 国产精品久久久久久精品古装| 如日韩欧美国产精品一区二区三区 | 久久国产乱子免费精品| 欧美成人精品欧美一级黄| 如日韩欧美国产精品一区二区三区 | 日本wwww免费看| 亚洲精品456在线播放app| 亚洲真实伦在线观看| 国产精品嫩草影院av在线观看| .国产精品久久| 亚洲精品久久午夜乱码| 人人妻人人看人人澡| 有码 亚洲区| 男女无遮挡免费网站观看| 99热网站在线观看| 亚洲不卡免费看| 人人妻人人澡人人看| 在线观看www视频免费| freevideosex欧美| 少妇人妻精品综合一区二区| 美女大奶头黄色视频| 99九九在线精品视频 | 亚洲成人手机| 看十八女毛片水多多多| 久久久欧美国产精品| 亚洲国产av新网站| 亚洲人成网站在线播| 日韩伦理黄色片| 欧美成人午夜免费资源| 国产伦精品一区二区三区视频9| 久久精品国产a三级三级三级| 国产伦精品一区二区三区视频9| 极品少妇高潮喷水抽搐| 看非洲黑人一级黄片| 久热久热在线精品观看| 国产欧美日韩综合在线一区二区 | 日本黄大片高清| 亚洲精品日本国产第一区| 国产免费一级a男人的天堂| 国产精品熟女久久久久浪| 午夜福利在线观看免费完整高清在| 一本一本综合久久| 午夜福利在线观看免费完整高清在| 一本久久精品| 国产黄片美女视频| 大片免费播放器 马上看| 黄色日韩在线| 好男人视频免费观看在线| 久久久国产精品麻豆| 久久韩国三级中文字幕| 欧美日韩综合久久久久久| 在线观看免费高清a一片| 大香蕉久久网| 99九九线精品视频在线观看视频| 欧美高清成人免费视频www| 亚洲情色 制服丝袜| 欧美激情国产日韩精品一区| 国产av国产精品国产| 黄色日韩在线| av国产久精品久网站免费入址| 久久女婷五月综合色啪小说| 亚洲成人手机| 色网站视频免费| 男男h啪啪无遮挡| 精品酒店卫生间| 亚洲高清免费不卡视频| 成人亚洲精品一区在线观看| 亚洲欧美日韩卡通动漫| 一个人看视频在线观看www免费| 一级毛片久久久久久久久女| 国产在线视频一区二区| 99久久综合免费| 色94色欧美一区二区| 热re99久久国产66热| 色网站视频免费| 午夜91福利影院| 成人亚洲欧美一区二区av| 精品国产国语对白av| 波野结衣二区三区在线| 毛片一级片免费看久久久久| 人人妻人人澡人人爽人人夜夜| 国产精品久久久久久av不卡| 日韩,欧美,国产一区二区三区| 天堂8中文在线网| 欧美高清成人免费视频www| 噜噜噜噜噜久久久久久91| 久久久欧美国产精品| 我要看日韩黄色一级片| 国产精品不卡视频一区二区| 色哟哟·www| 久热这里只有精品99| 欧美另类一区| 午夜免费男女啪啪视频观看| 激情五月婷婷亚洲| 亚洲欧美成人精品一区二区| 视频中文字幕在线观看| 人妻一区二区av| 免费观看性生交大片5| 欧美激情极品国产一区二区三区 | 久久人妻熟女aⅴ| 女性被躁到高潮视频| 亚洲激情五月婷婷啪啪| 在线观看www视频免费| 中国三级夫妇交换| 97超碰精品成人国产| 插阴视频在线观看视频| 91aial.com中文字幕在线观看| 大又大粗又爽又黄少妇毛片口| 精品少妇久久久久久888优播| 久久国产精品大桥未久av | 人体艺术视频欧美日本| 精品人妻熟女av久视频| 观看av在线不卡| 久久精品国产鲁丝片午夜精品| 啦啦啦在线观看免费高清www| 久久久欧美国产精品| 国产av国产精品国产| 少妇的逼水好多| 天美传媒精品一区二区| 国产高清不卡午夜福利| 日韩人妻高清精品专区| 曰老女人黄片| 成人黄色视频免费在线看| 国产亚洲av片在线观看秒播厂| 国产亚洲欧美精品永久| 午夜视频国产福利| 免费大片18禁| 国产av精品麻豆| 亚洲欧美一区二区三区国产| 国产精品成人在线| 国产乱人偷精品视频| 97超视频在线观看视频| 在线观看免费视频网站a站| 久久人人爽av亚洲精品天堂| 人人妻人人澡人人爽人人夜夜| 精品少妇内射三级| 久久久亚洲精品成人影院| 国产无遮挡羞羞视频在线观看| 亚洲精品,欧美精品|