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

    全球變暖形勢(shì)下中國(guó)陸表水分的變化

    2012-01-16 01:33:08張文君譚桂容
    大氣科學(xué)學(xué)報(bào) 2012年5期
    關(guān)鍵詞:土壤濕度水汽徑流

    張文君,譚桂容

    (南京信息工程大學(xué)1.氣象災(zāi)害省部共建教育部重點(diǎn)實(shí)驗(yàn)室;2.大氣科學(xué)學(xué)院,江蘇南京210044)

    0 引言

    作為地球系統(tǒng)中緩變過(guò)程之一,海表溫度對(duì)氣候變率的重要性已經(jīng)比較清楚,并有大量的研究工作(陳海山等,2002;Wu and Zhou,2008;Zhou et al.,2008,2009;李忠賢等,2011)。有觀點(diǎn)認(rèn)為,土壤濕度在氣候中的作用僅次于海表溫度,在中高緯,其作用甚至與海溫相當(dāng)(Shukla and Mintz,1982;Etin et al.,1999;Koster et al.,2000)。短期內(nèi),土壤濕度可以通過(guò)控制地表感熱和潛熱,影響邊界層的特性和對(duì)流的觸發(fā)條件、維持時(shí)間;長(zhǎng)期內(nèi),土壤濕度可調(diào)節(jié)干旱和洪澇(馬柱國(guó)等,2001;Pan et al.,2001;左志燕和張人禾,2007;梁樂(lè)寧和陳海山,2010)。大量數(shù)值試驗(yàn)已經(jīng)證明了土壤濕度對(duì)氣候有重要影響,如Shukla and Mintz(1982)利用模式對(duì)干、濕土壤的影響做了對(duì)比試驗(yàn),結(jié)果表明干、濕土壤對(duì)后期降水和氣溫的影響差異較大;Walker and Rowntree(1977)的敏感性試驗(yàn)表明,干土壤可使未來(lái)氣溫升高,濕土壤可使后期降水持續(xù);Kim and Hong(2007)利用區(qū)域氣候模式研究了東亞地區(qū)土壤濕度對(duì)夏季降水的影響,指出東亞區(qū)域的土壤濕度對(duì)降水有正反饋?zhàn)饔?張文君等(2012)通過(guò)較好的土壤濕度強(qiáng)迫區(qū)域氣候模式發(fā)現(xiàn),合理描述土壤濕度是提高中國(guó)夏季氣候預(yù)報(bào)技巧的潛在途徑之一;類(lèi)似的工作較多(Rind,1982;Yeh et al.,1984;王萬(wàn)秋,1991;林朝暉等,2001;胡婭敏等,2009)。

    全球變暖已成為國(guó)際上大多數(shù)科學(xué)家的共識(shí)(王紹武,1994)。目前普遍認(rèn)為,過(guò)去50 a的全球增暖可能主要是人類(lèi)活動(dòng)向大氣排放溫室氣體造成的(張?zhí)m生等,2000)。全球變暖會(huì)導(dǎo)致地表蒸發(fā)的增加,從而引發(fā)全球干旱化的加劇(Mather and Feddema,1986)。隨著溫室氣體的增加,使得向下的紅外輻射增強(qiáng),到達(dá)地表的能量增加,從而地球表面蒸發(fā)增加,由于大氣的三維運(yùn)動(dòng),蒸發(fā)的加強(qiáng)不僅導(dǎo)致局地降水的增加,而且可能導(dǎo)致其他區(qū)域降水的增加;對(duì)某一區(qū)域而言,如果降水大于潛在蒸發(fā),將導(dǎo)致徑流的增加,反之,降水的增加(減少)不會(huì)導(dǎo)致徑流的增強(qiáng)(減弱),取而代之的是土壤濕度的增加(減少)(Manabe et al.,2004)。那么在全球變暖的背景下,土壤濕度究竟如何變化呢?Wetherald and Manabe(2002)、Manabe et al.(2004)運(yùn)用GFDL的耦合模式研究了在增加溫室氣體和硫酸鹽氣溶膠(相當(dāng)于IPCC(政府間氣候變化委員會(huì))中的IS92a情景;Houghton et al.,1992)以及CO2增加四倍情景下土壤濕度的變化,結(jié)果表明,半干旱區(qū)土壤變干,北半球中高緯度土壤濕度是夏季變干、冬季變濕。

    近些年,耦合模式是評(píng)估溫室氣體增加對(duì)全球和區(qū)域氣候影響的重要工具,IPCC總結(jié)報(bào)告(Lakeman,1996;Houghton,2001)指出,雖然模式模擬的氣候變化在全球和大陸尺度上大體是一致的,但在區(qū)域細(xì)節(jié)上仍存在很大的不確定,所以有必要進(jìn)行多模式分析。IPCC第四次評(píng)估報(bào)告(IPCC AR4)組織了CO2加倍的模擬試驗(yàn),其中包括了陸面模式的結(jié)果,為研究全球變暖背景下中國(guó)土壤濕度的變化提供了基礎(chǔ)。

    本文利用可獲的10個(gè)耦合模式CO2加倍試驗(yàn)和控制試驗(yàn)的模擬結(jié)果,分析在全球變暖背景下中國(guó)土壤濕度的變化。本文其他部分安排如下:第一部分為資料和方法的簡(jiǎn)單介紹,第二部分為東亞季風(fēng)的變化,第三部分為降水和降水蒸發(fā)差的變化,第四部分為土壤濕度的變化及其可能原因,第五部分為結(jié)論和討論。

    1 資料和方法

    本文所用IPCC CO2加倍試驗(yàn)的具體做法為:模式大氣的CO2是在控制試驗(yàn)的基礎(chǔ)上,以1%/a的速率遞增,大約70 a后達(dá)到加倍,隨后固定CO2濃度再積分150 a。這里通過(guò)對(duì)CO2加倍試驗(yàn)最后80 a的平均結(jié)果和控制試驗(yàn)最后80 a的平均結(jié)果進(jìn)行比較,討論全球變暖對(duì)中國(guó)土壤濕度的影響。模式介紹見(jiàn)表1,詳細(xì)信息見(jiàn)http://www-pcmdi.llnl.gov。注意NCAR CCSM3和MRI_CGCM2.3.2的CO2加倍試驗(yàn)是從當(dāng)代氣候的控制試驗(yàn)開(kāi)始積分的,而其他模式的CO2加倍試驗(yàn)是從工業(yè)革命前的控制試驗(yàn)開(kāi)始積分的。

    這里分析表層10 cm土壤濕度的變化。此外,文中的“變化”是指CO2加倍試驗(yàn)最后80 a的平均值減去控制試驗(yàn)最后80 a的平均值。

    2 東亞季風(fēng)的變化

    溫室氣體的增加使得溫度升高,CO2加倍后中國(guó)氣溫顯著上升,年均增溫約2.7℃,增溫在東南最少,向東北和西部逐漸變強(qiáng),且冬季強(qiáng)于夏季(圖1)。溫度升高將使得環(huán)流場(chǎng)發(fā)生變化,圖2為CO2加倍后夏季850 hPa風(fēng)場(chǎng)的變化??梢?jiàn),多模式集合平均在中國(guó)東部區(qū)域表現(xiàn)為西南風(fēng)異常,這意味著夏季風(fēng)增強(qiáng),不過(guò)模式間依然存在差異;模擬夏季風(fēng)增強(qiáng)的模式有CGCM3.1(T47)、GISS_ER、FGOALS_g1.0、MIROC3.2(hires)、MIROC3.2(medres)、CCSM3和UKMO_HadCM3;模擬夏季風(fēng)減弱的模式僅有MRI_CGCM2.3.2模式。

    表1 模式情況簡(jiǎn)介T(mén)able 1 The models'information

    圖3為冬季850 hPa風(fēng)場(chǎng)變化。由多模式平均結(jié)果可見(jiàn),冬季氣流異常主要是由南海、印度洋以及西太平洋向中國(guó)的氣流異常,意味著冬季風(fēng)減弱。但模式間的差異較大,模擬冬季風(fēng)減弱的模式有GISS_ER、FGOALS_g1.0、UKMO_HadCM3和UKMO_HadGEM1,而CCSM3模擬的冬季風(fēng)增強(qiáng)。

    布和朝魯(2003)指出,由于全球變暖,大陸增溫強(qiáng)于海洋,從而使得夏季海陸溫差加大,夏季風(fēng)增強(qiáng),而冬季海陸溫差減小,使得冬季風(fēng)減弱,本文結(jié)果與之一致。夏季風(fēng)增強(qiáng)使得印度西南季風(fēng)流經(jīng)孟加拉灣向中國(guó)的水汽輸送、沿105°E附近流經(jīng)我國(guó)南海的跨赤道水汽輸送、來(lái)自西太平洋菲律賓的水汽輸送均增強(qiáng)(圖略)。雖然大部分模式表現(xiàn)出冬季風(fēng)減弱,但是其水汽場(chǎng)變化與風(fēng)場(chǎng)變化并不一致。由多模式平均結(jié)果可見(jiàn),在中國(guó)30°N以南主要為西太平洋菲律賓和印度向中國(guó)的水汽輸送異常,而在中國(guó)30°N以北則為來(lái)自西北方向的弱水汽輸送異常,除了INM_CM3.0和CCSM3外,所有模式與多模式平均基本一致(圖4)。總之,不論夏季還是冬季,向中國(guó)區(qū)域輸送的水汽增加了。

    3 降水和降水蒸發(fā)差的變化

    向中國(guó)輸送的水汽增加有利于中國(guó)降水的增加,圖5為中國(guó)夏季降水的變化。由多模式平均結(jié)果看,中國(guó)夏季降水除了長(zhǎng)江流域外基本都增加了,但模式間的差別較大;CGCM3.1(T47)、GISS_ER、MIROC3.2(hires)、CCSM3和UKMO_HadCM3模式模擬的降水變化與多模式集合平均結(jié)果基本類(lèi)似,不過(guò)INM_CM3.0和MRI_CGCM2.3.2模式模擬的夏季降水在長(zhǎng)江、華北及西北基本都減小,F(xiàn)GOALS_g1.0模式模擬的降水在中部和華南減小。

    圖1 中國(guó)夏季(a)和冬季(b)溫度的變化(10個(gè)模式的平均結(jié)果;陰影區(qū)表示通過(guò)0.05信度的顯著性檢驗(yàn);單位:℃)Fig.1 Changes of the temperature over China in(a)summer(JJA)and(b)winter(DJF)(The results are the ensemble mean of 10 coupled models.Shaded regions are significant at 95%confidence level.Units:℃)

    圖6為中國(guó)冬季降水的變化。由多模式結(jié)果可見(jiàn),除了華南外,中國(guó)冬季降水都增加,且在江淮流域存在一個(gè)增加中心,除了GISS_ER和MIROC3.2(hires)模式模擬的降水在華南明顯增加外,其余模式與此結(jié)果基本一致,只不過(guò)增值中心略有不同而已。

    隨著全球增溫,到達(dá)地面的能量增加,使得蒸發(fā)增加,CO2加倍后中國(guó)蒸發(fā)不論夏季還是冬季都一致地增加,且蒸發(fā)增加量相對(duì)于降水而言空間分布較均勻(圖略),但是作為陸面水分來(lái)源之一的降水蒸發(fā)差(P-E)變化如何呢?夏季P-E變化的集合平均結(jié)果除了在東北和南方增加外,從長(zhǎng)江流域一直到西北有一帶狀減小區(qū)域,但由于模式間降水的差別較大,導(dǎo)致夏季P-E在模式間的差別也較大;MRI_CGCM2.3.2和UKMO_HadCM3模式結(jié)果與集合平均結(jié)果基本類(lèi)似,MIROC3.2(hires)和UKMO_HadGEM1模式結(jié)果通過(guò)顯著性檢驗(yàn)的區(qū)域很小,其余模式差別較大(圖7)。冬季P-E變化的集合平均結(jié)果是北方增加、南方減小,由于各模式模擬的冬季降水變化較一致,所以各模式P-E變化除了GISS_ER模式外,基本與集合平均結(jié)果類(lèi)似(圖8)。

    4 土壤濕度的變化及其可能原因

    降水的增加有利于土壤濕度的增加,而蒸發(fā)的增加則不利于土壤濕度的增加,土壤濕度的變化主要是受降水的影響還是受P-E的影響?圖9為耦合模式模擬的夏季中國(guó)土壤濕度變化。多模式集合平均結(jié)果表現(xiàn)為,除在西北和華北的東部土壤有較弱的變濕外,基本都變干;除了INM_CM3.0和UKMO_HadGEM1模式模擬的土壤基本都變干外,其余模式模擬結(jié)果與多模式結(jié)果基本類(lèi)似,此外變濕的區(qū)域基本在干旱區(qū)。圖10為耦合模式模擬的冬季中國(guó)土壤濕度的變化,除了UKMO_HadGEM1模式模擬的土壤基本都變干外,幾乎所有耦合模式都一致表現(xiàn)為干旱區(qū)土壤濕度增加,且增加的范圍和強(qiáng)度遠(yuǎn)大于夏季,其他區(qū)域土壤變干,但GISS模式模擬的南方土壤卻變濕了。

    上述分析表明,土壤濕度的變化與降水的變化并不一致。夏季,雖然多模式平均結(jié)果的土壤濕度變化與P-E的變化也不一致,但FGOALS_g1.0、INM_CM3.0、CCSM3和UKMO_HadCM3的土壤濕度變化與P-E變化基本一致,CGCM3.1(T47)和MRI_CGCM2.3.2除了西部和東北外,土壤濕度變化與P-E的變化基本一致;冬季,土壤濕度變化與P-E的變化較一致。華文劍和陳海山(2011)也分析了土壤濕度對(duì)全球變暖的響應(yīng),他們研究結(jié)果也表明亞洲中北部干旱區(qū)土壤濕度增加而南部區(qū)域土壤濕度減少。盡管本研究的結(jié)果與他們的結(jié)論定性一致,但是區(qū)域分布仍存在著較大的差異,其主要原因是所選取的模式存在差異性。這也說(shuō)明了多模式集合平均模式結(jié)果仍存在著一定的不確定性。

    圖2 各模式夏季850 hPa風(fēng)場(chǎng)變化(Ensemble為10個(gè)模式的平均結(jié)果;陰影區(qū)表示通過(guò)0.05信度的顯著性檢驗(yàn);單位:mm/s)Fig.2 The 850 hPa wind changes in summer(JJA)(“Ensemble”denotes the ensemble mean of 10 coupled models.Shaded regions are significant at 95%confidence level.Units:mm/s)

    一些學(xué)者分析了全球變暖的背景下夏季干旱的原因(Manabe et al.,1981;Manabe and Wetherald,1987;Mitchell and Warrilow,1987;Wetherald and Manabe,1999),Wetherald and Manabe(2002)在此基礎(chǔ)上分析了中高緯土壤濕度夏季變干、冬季變濕的物理機(jī)制,指出:在海洋上,溫度的季節(jié)變化小,伴隨著全球變暖,海表面飽和水汽壓增加,因此全年蒸發(fā)都增加,增加的水汽通過(guò)大氣環(huán)流使得鄰近大陸上的降水增加;與海洋上相反,在大陸上,季節(jié)循環(huán)溫度變化較大,冬季溫度很低,夏季較高,在低溫下,由于Clausius-Clapeyron方程非線性,盡管溫度增加,飽和水汽壓幾乎不增加,因此,用于到達(dá)地表紅外輻射增加的能量主要用于感熱的增加,而非潛熱的增加。另一方面,由于海洋輸送到大陸的水汽使得降水增加,所以冬季土壤濕度增加。在夏季,不管海洋還是大陸表面溫度和飽和水汽壓都很高,地表增加的能量主要用于蒸發(fā)的增加;另一方面,在中緯度,相同緯度的大陸上降水的增加小于海洋上降水的增加,這就使得土壤濕度變干。這也可能使得中國(guó)干旱區(qū)冬季變濕范圍大、強(qiáng)度強(qiáng),而夏季變濕范圍小、強(qiáng)度弱。對(duì)于干旱區(qū)土壤變濕,而其余區(qū)域變干,Manabe and Wetherald(1987)認(rèn)為由于CO2加倍導(dǎo)致紅外輻射增加使得蒸發(fā)幾乎都增加,而在半干旱區(qū)降水的增加相對(duì)較小,從而土壤變干;而在干旱區(qū)蒸發(fā)和降水的增加都很小,溫度的增加可能導(dǎo)致零星的對(duì)流運(yùn)動(dòng),從而增加降水和土壤濕度。

    圖3 各模式冬季850 hPa風(fēng)場(chǎng)變化(Ensemble為10個(gè)模式的平均結(jié)果;陰影區(qū)表示通過(guò)0.05信度的顯著性檢驗(yàn);單位:mm/s)Fig.3 The 850 hPa wind changes in winter(DJF)(“Ensemble”denotes the ensemble mean of 10 coupled models.Shaded regions are significant at 95%confidence level.Units:mm/s)

    由陸表水分平衡方程可知,土壤濕度的變化受降水、蒸發(fā)和徑流的影響,通常在此忽略積雪的影響,在緩變氣候中,此近似是成立的,但在CO2加倍的情景下,積雪的影響就不能忽略了,尤其夏季高原積雪融化給西北補(bǔ)充大量的水分,而冬季雖然溫度也升高了,但實(shí)際溫度仍較低,積雪的變化不大,所以冬季P-E的變化與土壤濕度變化較一致。

    為了研究陸表水循環(huán)的整體性,在此本文簡(jiǎn)單分析一下多模式年平均統(tǒng)計(jì)結(jié)果。表2給出了兩個(gè)典型區(qū)域的陸表水循環(huán)分量的變化。不論南方還是西北,降水、蒸發(fā)、徑流都增加了,控制試驗(yàn)中總徑流和蒸發(fā)之和與降水量幾乎是平衡的;CO2加倍試驗(yàn)中,西北地區(qū)總徑流和蒸發(fā)之和與降水量差別很小,但南方差別較大。由表3可見(jiàn),西北徑流和蒸發(fā)的增加為4.1 cm/a,而降水增加為3.9 cm/a,降水的增加不足以滿足徑流和蒸發(fā)的增加,但在西北積雪融化給該地區(qū)補(bǔ)充了大量的水分(圖略),使得土壤濕度增加;南方徑流和蒸發(fā)的增加遠(yuǎn)大于降水的增加,因?yàn)闆](méi)有別的水分補(bǔ)充,只有通過(guò)土壤水分來(lái)維持陸表水平衡,使得土壤變干。

    圖4 各模式冬季850 hPa水汽場(chǎng)變化(Ensemble為10個(gè)模式的平均結(jié)果;陰影區(qū)表示通過(guò)0.05信度的顯著性檢驗(yàn);單位:(g/kg)·(m/s))Fig.4 The water vapor transport changes at 850 hPa in winter(DJF)(“Ensemble”denotes the ensemble mean of 10 coupled models.Shaded regions are significant at 95%confidence level.Units:(g/kg)·(m/s))

    表2 多模式年平均的中國(guó)西北(75~102.5°E,35~45°N)、南方(102.5~120°E,22~30°N)降水、徑流、蒸發(fā)和土壤濕度(徑流指地表徑流,總徑流指地表徑流與地下徑流之和)Table 2 Precipitation,runoff,evaporation and soil moisture over Northwest China(35—45°N,75—102.5°E)and South China(22—30°N,102.5—120°E)(The results are the ensemble mean of 10 coupled models.Runoff denotes surface runoff and total runoff denotes the sum of surface runoff and underground runoff )

    圖5 各模式夏季降水的變化(Ensemble為10個(gè)模式的平均結(jié)果;陰影區(qū)表示通過(guò)0.05信度的顯著性檢驗(yàn);單位:mm/d)Fig.5 The precipitation changes in summer(JJA)(“Ensemble”denotes the ensemble mean of 10 coupled models.Shaded regions are significant at 95%confidence level.Units:mm/d)

    表3 表2中CO2加倍試驗(yàn)與控制試驗(yàn)的差值Table 3 Differences of the double CO2and the control experiments in Table 2cm·a-1

    8 結(jié)論和討論

    本文通過(guò)對(duì)10個(gè)耦合模式的CO2加倍試驗(yàn)與控制實(shí)驗(yàn)結(jié)果的比較,討論了全球變暖對(duì)中國(guó)土壤濕度的潛在影響,結(jié)果表明:

    1)隨著全球變暖,7個(gè)模式模擬的夏季風(fēng)增強(qiáng)(CGCM3.1(T47)、GISS_ER、FGOALS_g1.0、MIROC3.2(hires)、MIROC3.2(medres)、CCSM3和UKMO_HadCM3),而MRI_CGCM2.3.3模式模擬的夏季風(fēng)減弱;4個(gè)模式模擬的冬季風(fēng)減弱(GISS_ER、FGOALS_g1.0、UKMO_HadCM3和UKMO_HadGEM1),但CCSM3模式模擬的冬季風(fēng)增強(qiáng)。

    2)在夏季,由夏季風(fēng)帶給中國(guó)的水汽輸送增強(qiáng);在冬季,30°N以南為西太平洋菲律賓和印度向中國(guó)的水汽輸送異常,30°N以北為由西北而來(lái)的水汽輸送異常。

    圖6 各模式冬季降水的變化(Ensemble為10個(gè)模式的平均結(jié)果;陰影區(qū)表示通過(guò)0.05信度的顯著性檢驗(yàn);單位:mm/d)Fig.6 The precipitation changes in winter(DJF)(“Ensemble”denotes the ensemble mean of 10 coupled models.Shaded regions are significant at 59%confidence level.Units:mm/d)

    3)夏季降水除長(zhǎng)江流域外基本都增加,但模式間差別較大,CGCM3.1(T47)、GISS_ER、MIROC3.2(hires)、CCSM3和UKMO_HadCM3模式模擬的降水變化與多模式集合平均結(jié)果基本類(lèi)似,不過(guò)INM_CM3.0和MRI_CGCM2.3.2模式模擬的夏季降水在長(zhǎng)江、華北及西北基本都減小,而FGOALS_g1.0模式模擬的夏季降水在中國(guó)中部和華南減小;冬季,除了GISS_ER和MIROC3.2(hires)模式模擬的降水在華南為明顯增加外,其余模式模擬的降水基本一致表現(xiàn)為除華南外均增加。

    4)夏季P-E變化的集合平均結(jié)果除了在東北和南方增加外,從長(zhǎng)江流域一直到西北有一帶狀減小區(qū)域,其中FGOALS_g1.0、MRI_CGCM2.3.2和UKMO_HadCM3模式結(jié)果與集合平均結(jié)果基本一致,MIROC3.2(hires)和UKMO_HadGEM1模式結(jié)果通過(guò)顯著性檢驗(yàn)的區(qū)域很小,其余模式差別較大;冬季,幾乎所有模式的P-E均表現(xiàn)為北方增加、南方減小。

    圖7 各模式夏季P-E的變化(Ensemble為10個(gè)模式的平均結(jié)果;陰影區(qū)表示通過(guò)0.05信度的顯著性檢驗(yàn);單位:mm/d)Fig.7 The P-E changes in summer(JJA)(“Ensemble”denotes the ensemble mean of 10 coupled models.Shaded regions are significant at 95%confidence level.Units:mm/d)

    5)夏季,盡管模式間存在一些差異(如INM_CM3.0和UKMO_HadGEM1模式模擬的土壤基本都變干),但是模式基本表現(xiàn)為,除了在干旱區(qū)土壤有較弱的變濕外都變干;冬季,除了UKMO_Had-GEM1模式外,幾乎所有耦合模式都一致表現(xiàn)為干旱區(qū)土壤濕度增加,其他區(qū)域土壤變干,且增加的范圍和強(qiáng)度大于夏季。

    6)夏季,F(xiàn)GOALS_g1.0、INM_CM3.0、CCSM3和UKMO_HadCM3模式模擬的土壤濕度變化與PE變化基本一致,CGCM3.1(T47)和MRI_CGCM2.3.2模式除了西部和東北外,模擬的土壤濕度變化與P-E變化基本一致;冬季,所有模式模擬的土壤濕度變化與P-E的變化較一致。

    7)由陸表水循環(huán)量的分析可知,土壤濕度是陸表水循環(huán)量變化的綜合結(jié)果,在全球變暖的背景下,中國(guó)蒸發(fā)、降水和徑流基本都增加了,它們的綜合結(jié)果以及積雪的作用使得土壤濕度在干旱區(qū)增加。

    受復(fù)雜海陸分布和高大地形的影響,東亞氣候及其變化有其獨(dú)特性,目前的全球氣候模式對(duì)東亞氣候的模擬能力整體上還難以令人滿意(Yu et al.,2000;Zhou and Li,2002;Zhou and Yu,2006;吳蓉和張耀存,2012)。張宏芳和陳海山(2011a,2011b)對(duì)21個(gè)氣候模式模擬東亞夏季環(huán)流做了系統(tǒng)的評(píng)估,結(jié)果表明,盡管大部分模式對(duì)氣候態(tài)具有一定的模擬能力,但是對(duì)年際變率的模擬能力很弱。張文君等(2008)評(píng)估了14個(gè)全球?!憽?dú)怦詈夏J侥M的中國(guó)區(qū)域的土壤濕度,發(fā)現(xiàn)其模擬能力有限。由本文結(jié)果也可以看到,不同模式模擬的未來(lái)陸表水循環(huán)差異較大,盡管采用多模式集合平均盡量避免了模式間的誤差,但是不確定性仍較大,具有一定的模式依耐性。為了能更準(zhǔn)確地預(yù)估未來(lái)的氣候變化,提高模式模擬性能是當(dāng)務(wù)之急,一方面要不斷地提高模式分辨率,另一方面更需要關(guān)注模式相關(guān)物理過(guò)程及其參數(shù)化方案的改進(jìn)。

    圖8 各模式冬季P-E的變化(Ensemble為10個(gè)模式的平均結(jié)果;陰影區(qū)表示通過(guò)0.05信度的顯著性檢驗(yàn);單位:mm/d)Fig.8 The P-E changes in winter(DJF)(“Ensemble”denotes the ensemble mean of 10 coupled models.Shaded regions are significant at 95%confidence level.Units:mm/d)

    圖9 各模式夏季土壤濕度的變化(Ensemble為10個(gè)模式的平均結(jié)果;陰影區(qū)表示通過(guò)0.05信度的顯著性檢驗(yàn);單位:mm)Fig.9 The soil moisture changes in summer(JJA)(“Ensemble”denotes the ensemble mean of 10 coupled models.Shaded regions are significant at 95%confidence level.Units:mm)

    布和朝魯.2003.東亞季風(fēng)氣候未來(lái)變化的情景分析——基于IPCC SRES A2和B2方案的模擬結(jié)果[J].科學(xué)通報(bào),48(7):737-742.

    陳海山,孫照渤,倪東鴻,等.2002.東亞冬季風(fēng)對(duì)秋、冬季SSTA響應(yīng)的數(shù)值試驗(yàn)[J].南京氣象學(xué)院學(xué)報(bào),25(6):721-730.

    胡婭敏,丁一匯,廖菲.2009.土壤濕度資料同化對(duì)中國(guó)東部夏季區(qū)域氣候模擬的改進(jìn)[J].科學(xué)通報(bào),54(16):2388-2394.

    華文劍,陳海山.2011.陸面過(guò)程對(duì)全球變暖的響應(yīng)及可能機(jī)制——基于CMIP3的多模式集合分析[J].大氣科學(xué),35(1):121-133.

    李忠賢,陳海山,曾剛,等.2011.海溫強(qiáng)迫下的東亞夏季大氣環(huán)流潛在可預(yù)報(bào)性特征[J].大氣科學(xué)學(xué)報(bào),34(3):281-287.

    梁樂(lè)寧,陳海山.2010.春季華南土壤濕度異常與中國(guó)夏季降水的可能聯(lián)系[J].大氣科學(xué)學(xué)報(bào),33(5):536-546.

    林朝暉,楊小松,郭裕福.2001.陸面過(guò)程模式對(duì)土壤含水量初值的敏感性研究[J].氣象與環(huán)境研究,6(2):240-248.

    馬柱國(guó),符淙斌,謝力,等.2001.土壤濕度和氣候變化關(guān)系研究中的某些問(wèn)題[J].地球科學(xué)進(jìn)展,16(4):563-568.

    王紹武.1994.近百年氣候變化與變率的診斷研究[J].氣象學(xué)報(bào),52(3):261-273.

    王萬(wàn)秋.1991.土壤濕度異常對(duì)短期氣候影響的數(shù)值模擬試驗(yàn)[J].大氣科學(xué),15(5):115-123.

    吳蓉,張耀存.2012.RegCM3模式對(duì)長(zhǎng)江流域梅雨期降水和環(huán)流形勢(shì)的模擬性能檢驗(yàn)[J].氣象科學(xué),32(2):119-126.

    張宏芳,陳海山.2011a.21個(gè)氣候模式對(duì)東亞夏季環(huán)流模擬的評(píng)估I:氣候態(tài)[J].氣象科學(xué),31(2):119-128.

    張宏芳,陳海山.2011b.21個(gè)氣候模式對(duì)東亞夏季環(huán)流模擬的評(píng)估II:年際變化[J].氣象科學(xué),31(3):247-257.

    張?zhí)m生,方修琦,任國(guó)玉.2000.全球變化[M].北京:高等教育出版社.

    張文君,宇如聰,周天軍.2008.中國(guó)土壤濕度的分布與變化II.耦合模式模擬結(jié)果評(píng)估[J].大氣科學(xué),32(5):1128-1146.

    張文君,周天軍,智海.2012.土壤濕度影響中國(guó)夏季氣候的數(shù)值試驗(yàn)[J].氣象學(xué)報(bào),70(1):78-90.

    圖10 各模式冬季土壤濕度的變化(Ensemble為10個(gè)模式的平均結(jié)果;陰影區(qū)表示通過(guò)0.05信度的顯著性檢驗(yàn);單位:mm)Fig.10 The soil moisture changes in winter(DJF)(“Ensemble”denotes the ensemble mean of 10 coupled models.Shaded regions are significant at 95%confidence level.Units:mm)

    左志燕,張人禾.2007.中國(guó)東部夏季降水與春季土壤濕度的聯(lián)系[J].科學(xué)通報(bào),52(14):1722-1724.

    Cox P M,Betts R A,Bunton C B,et al.1999.The impact of new land surface physics on the GCM simulation of climate and climate sensitivity[J].Climate Dyn,15:183-203.

    Essery R,Best M,Cox P M.2001.MOSES 2.2 technical documentation[M]//Hadley Centre technical note 30.

    Etin J,Alan R,Vinnikov K Y,et al.1999.Evaluation of global soil wetness project soil moisture simulations[J].J Meteor Soc Japan,77(1B):183-198.

    Hasumi H,Emori S.2004.K-1 Coupled GCM(MIROC)description[M]//K-1 Technical Report No.1.

    Houghton J T.2001.Intergovernmental Panel on Climate Change(IPCC),climate change 2001:The science of climate change[M].New York:Cambridge University Press:572.

    Houghton J T,Callander B A,Varney S K.1992.IPCC(Intergovernmental Panel on Climate Change)climate change 1992:The supplementary report to the IPCC scientific assessment[M].New York:Cambridge University Press:865-871.

    Kim J E,Hong S Y.2007.Impact of soil moisture anomalies on summer rainfall over East Asia:A regional climate model study[J].J Climate,20:5732-5743.

    Koster R D,Suarez M J,Heiser M.2000.Variance and predictability of precipitation at seasonal-to-in terannual timescales[J].J Hydro meteor,1:26-46.

    Lakeman A.1996.Intergovernmental Panel on Climate Change(IPCC),cliamte change 1995:The science of climate change[M].New York:Cambridge University Press:572.

    Manabe S,Wetherald R T.1987.Large-scale change in soil wetness induced by an increase in atmospheric carbon dioxide[J].J Atmos Sci,44:1211-1235.

    Manabe S,Wetherald R T,Stouffer R J.1981.Summer dryness due to increase of atmospheric CO2concentration[J].Climatic Change,3:347-384.

    Manabe S,Milly P C,Wetherald R T.2004.Simulated long-term changes in river discharge and soil moisture due to global warming[J].Hydrol Sci,49(4):625-642.

    Mather J R,F(xiàn)eddema J.1986.‘Hydrologic consequence of increase in trace gases and CO2in the atmosphere’in effects of changes in stratospheric ozone and global climate,volume 3——Climate change[M].Washington D C:U.S.Environmental Protection A-gency:251-271.

    Mitchell J F B,Warrilow D A.1987.Summer dryness in northern midlatitudes due to increased CO2[J].Nature,330(19):238-240.

    Oleson K W,Dai Yongjiu,Bonan G,et al.2004.The technical description of the Community Land Model(CLM)[M]//NCAR Technical Note.

    Pan Z,Arritt R W,Gutowski W J,et al.2001.Soil moisture in a regional climate model:Simulation and projection[J].Geophys Res Lett,28(15):2947-2950.

    Radakovich J D,Wang Guiling,Chem Jiun-dar,et al.2003.Implementation of the NCAR Community Land Model(CLM)in the NASA/NCAR finite-volume climate model(FVGCM)[M]//14th symposium on global change and climate variations(compact view).

    Rind D.1982.The influence of ground moisture conditions in North A-merica on summer climate as modeled in the GISS GCM[J].Mon Wea Rev,110:1487-1494.

    Rosenzweig C,Abramopoulos F.1996.Land-surface model development for the GISS GCM[J].J Climate,10:2040-2054.

    Shukla J,Mintz Y.1982.Influence of land-surface evapotranspiration on the earth’s climate[J].Science,215:1498-1501.

    Verseghy D L.1991.CLASS:A Canadian land surface scheme for GCMs.I.Soil model[J].Int J Climatol,11:111-133.

    Verseghy D L,McFarlane N A,Lazare M.1993.CLASS:A Canadian Land Surface Scheme for GCMs.Part II:Vegetation model and coupled runs[J].Int J Climatol,13:347-370.

    Walker J,Rowntree P R.1977.The effect of soil moisture on circulation and rainfall in a tropical model[J].Quart J Roy Meteor Soc,103:29-46.

    Wetherald R T,Manabe S.1999.Detectability of summer dryness caused by greenhouse warming[J].Climatic Change,43:495-511.

    Wetherald R T,Manabe S.2002.Simulation of hydrologic changes associated with global warming[J].J Geophys Res,107:4379-4394.

    Wu B,Zhou T.2008.Oceanic origin of the interannual and interdecadal variability of the summertime western Pacific subtropical high[J].Geophys Res Lett,35,L13701.doi:10.1029/2008GL034584.

    Yeh T C,Wetherald R I,Manabe S.1984.The effect of soil moisture on the shorttern climate and hydrology change:A numerical experiment[J].Mon Wea Rev,112:474-490.

    Yu R,Li W,Zhang X,et al.2000.Climatic features related to eastern China summer rainfalls in the NCAR CCM3[J].Adv Atmos Sci,17:503-518.

    Yukimoto S,Noda A,Kitoh A,et al.2001.The new Meteorological Reasearch Institute Coupled GCM(MRI-CGCM2)——model climate and variability[J].Pap Meteor Geopyhs,51(2):27-88.

    Zhou T,Li Z.2002.Simulation of the East Asian summer monsoon by using a variable resolution atmospheric GCM[J].Climate Dyn,19:167-180.

    Zhou T,Yu R.2006.Twentieth century surface air temperature over China and the globe simulated by coupled climate models[J].J Climate,19:5843-5858.

    Zhou T,Yu R,Li H,et al.2008.Ocean forcing to changes in global monsoon precipitation over the recent half century[J].J Climate,21:3833-3852.

    Zhou T,Wu B,Wang B.2009.How well do Atmospheric General Circulation Models capture the leading modes of the in terannual variability of Asian-Australian Monsoon?[J].JClimate,22:1159-1173.

    猜你喜歡
    土壤濕度水汽徑流
    青藏高原上空平流層水汽的時(shí)空演變特征
    土壤濕度傳感器在園林綠化灌溉上的應(yīng)用初探
    基于51單片機(jī)控制花盆土壤濕度
    電子制作(2019年15期)2019-08-27 01:12:12
    1979~2011年間平流層溫度及平流層水汽的演變趨勢(shì)
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    四川盆地土壤濕度時(shí)空分布及影響因子分析
    中國(guó)不同氣候區(qū)土壤濕度特征及其氣候響應(yīng)
    深圳“5·11”特大暴雨過(guò)程的水汽輸送特征分析
    探秘“大徑流”
    攻克“大徑流”
    国产毛片a区久久久久| 欧美激情久久久久久爽电影| 国产伦精品一区二区三区四那| 联通29元200g的流量卡| 联通29元200g的流量卡| 日韩欧美在线乱码| 亚洲avbb在线观看| 亚洲avbb在线观看| 成人亚洲精品av一区二区| 亚洲一区高清亚洲精品| 亚洲精品456在线播放app | 中文字幕免费在线视频6| 成年女人毛片免费观看观看9| 热99re8久久精品国产| 麻豆精品久久久久久蜜桃| 男人狂女人下面高潮的视频| 亚洲美女搞黄在线观看 | 欧美高清性xxxxhd video| 亚洲久久久久久中文字幕| 级片在线观看| 中文字幕久久专区| 亚洲不卡免费看| 99久国产av精品| 欧美中文日本在线观看视频| 国产精品乱码一区二三区的特点| 97人妻精品一区二区三区麻豆| 亚洲av日韩精品久久久久久密| 国产精品久久视频播放| 亚洲欧美清纯卡通| 成人av在线播放网站| 久久九九热精品免费| 18禁裸乳无遮挡免费网站照片| 久久九九热精品免费| 十八禁国产超污无遮挡网站| 亚洲经典国产精华液单| 97碰自拍视频| 日本 av在线| 精品午夜福利在线看| 国模一区二区三区四区视频| 九九爱精品视频在线观看| 成人一区二区视频在线观看| 国产亚洲精品综合一区在线观看| 国产精品乱码一区二三区的特点| 久久精品久久久久久噜噜老黄 | 国产免费男女视频| 少妇丰满av| 岛国在线免费视频观看| 亚洲欧美精品综合久久99| 午夜久久久久精精品| 免费人成在线观看视频色| 国产精品亚洲美女久久久| 欧美日本亚洲视频在线播放| 日本精品一区二区三区蜜桃| 岛国在线免费视频观看| 国产69精品久久久久777片| 赤兔流量卡办理| 国产高潮美女av| 五月玫瑰六月丁香| 免费搜索国产男女视频| 女人被狂操c到高潮| 日本一本二区三区精品| 中国美白少妇内射xxxbb| 变态另类成人亚洲欧美熟女| 亚洲成人久久性| 男女啪啪激烈高潮av片| 国产 一区精品| 日本精品一区二区三区蜜桃| av在线天堂中文字幕| 亚洲av第一区精品v没综合| 真人一进一出gif抽搐免费| 18禁黄网站禁片午夜丰满| 亚洲人成网站在线播| 精品久久久久久久久久免费视频| 成人一区二区视频在线观看| 国产白丝娇喘喷水9色精品| 亚洲人成网站在线播放欧美日韩| 欧美+亚洲+日韩+国产| 久久精品国产清高在天天线| 亚洲无线在线观看| 精品午夜福利视频在线观看一区| 九九在线视频观看精品| 免费不卡的大黄色大毛片视频在线观看 | 97碰自拍视频| 91在线精品国自产拍蜜月| 欧美极品一区二区三区四区| 免费观看在线日韩| 精品乱码久久久久久99久播| 在线观看美女被高潮喷水网站| 欧美成人性av电影在线观看| 深夜a级毛片| 一级黄色大片毛片| 成人国产一区最新在线观看| 丰满人妻一区二区三区视频av| a级毛片a级免费在线| 欧美日韩乱码在线| 日本欧美国产在线视频| 人人妻人人看人人澡| 国产老妇女一区| 在线观看av片永久免费下载| 国产不卡一卡二| 国产91精品成人一区二区三区| 一本久久中文字幕| 看黄色毛片网站| 全区人妻精品视频| 国产亚洲精品久久久久久毛片| 2021天堂中文幕一二区在线观| 久久久久久久久久黄片| 国产精品日韩av在线免费观看| 嫩草影院精品99| 中国美女看黄片| 老师上课跳d突然被开到最大视频| 亚洲人成网站在线播放欧美日韩| 久久久久久大精品| 中文字幕免费在线视频6| 少妇的逼好多水| 午夜激情福利司机影院| 亚洲av五月六月丁香网| 亚洲国产欧洲综合997久久,| 干丝袜人妻中文字幕| 久久久久国产精品人妻aⅴ院| 国产高清激情床上av| 日韩欧美精品免费久久| 国产aⅴ精品一区二区三区波| 日本熟妇午夜| 日韩人妻高清精品专区| 综合色av麻豆| 亚洲国产色片| 一本精品99久久精品77| 国产伦在线观看视频一区| 91av网一区二区| 性插视频无遮挡在线免费观看| 九色国产91popny在线| 国内少妇人妻偷人精品xxx网站| 欧美性猛交黑人性爽| 久久精品国产99精品国产亚洲性色| 少妇熟女aⅴ在线视频| 国产精品av视频在线免费观看| 欧美性感艳星| 免费观看的影片在线观看| 美女 人体艺术 gogo| 日本三级黄在线观看| 国产伦在线观看视频一区| 国产一区二区三区视频了| 日本 av在线| 国内毛片毛片毛片毛片毛片| 精品一区二区三区视频在线| 91麻豆av在线| 老师上课跳d突然被开到最大视频| 国产探花极品一区二区| 国产精品一及| 非洲黑人性xxxx精品又粗又长| 特级一级黄色大片| 美女高潮喷水抽搐中文字幕| 国产精品久久久久久久久免| 亚洲中文字幕一区二区三区有码在线看| 女生性感内裤真人,穿戴方法视频| 久久国内精品自在自线图片| 久久精品国产99精品国产亚洲性色| 联通29元200g的流量卡| 在线观看av片永久免费下载| 午夜福利在线观看吧| 国内精品久久久久精免费| 亚洲精品一区av在线观看| 男女视频在线观看网站免费| 中文字幕av在线有码专区| eeuss影院久久| www.www免费av| 真人做人爱边吃奶动态| 婷婷六月久久综合丁香| 韩国av在线不卡| 黄色欧美视频在线观看| 99视频精品全部免费 在线| 九九久久精品国产亚洲av麻豆| 日本爱情动作片www.在线观看 | 午夜激情欧美在线| 亚洲三级黄色毛片| 婷婷精品国产亚洲av| 天堂动漫精品| 女的被弄到高潮叫床怎么办 | 一级黄片播放器| 婷婷亚洲欧美| 亚洲国产欧洲综合997久久,| 最近在线观看免费完整版| 在线播放国产精品三级| 不卡视频在线观看欧美| 3wmmmm亚洲av在线观看| 一夜夜www| 亚洲精品456在线播放app | 狂野欧美激情性xxxx在线观看| 亚洲最大成人手机在线| 欧洲精品卡2卡3卡4卡5卡区| 免费看光身美女| 国产一区二区三区av在线 | 成年版毛片免费区| 精品人妻偷拍中文字幕| 亚洲精品影视一区二区三区av| 男女那种视频在线观看| av在线蜜桃| 国产av不卡久久| 禁无遮挡网站| 久久国产精品人妻蜜桃| 日韩精品中文字幕看吧| 日本五十路高清| 色哟哟·www| 1024手机看黄色片| 久久久久免费精品人妻一区二区| 精品午夜福利在线看| 一级av片app| 亚洲专区国产一区二区| 99久久精品热视频| 久久午夜亚洲精品久久| 久久精品国产鲁丝片午夜精品 | 日本欧美国产在线视频| 黄色欧美视频在线观看| 久久久久国产精品人妻aⅴ院| 美女免费视频网站| 欧美一区二区亚洲| 性欧美人与动物交配| 特级一级黄色大片| 男人舔女人下体高潮全视频| 91久久精品国产一区二区成人| 亚洲avbb在线观看| 欧美日韩综合久久久久久 | 中文字幕精品亚洲无线码一区| 老师上课跳d突然被开到最大视频| 网址你懂的国产日韩在线| 精品国产三级普通话版| 欧美成人a在线观看| 直男gayav资源| 99视频精品全部免费 在线| 男女下面进入的视频免费午夜| 亚洲av熟女| 国产精品久久久久久久久免| 国产毛片a区久久久久| 男女下面进入的视频免费午夜| 亚洲精华国产精华精| 校园春色视频在线观看| 欧美不卡视频在线免费观看| 亚洲欧美日韩无卡精品| 国产精华一区二区三区| 免费av不卡在线播放| 综合色av麻豆| 中文字幕人妻熟人妻熟丝袜美| 成人综合一区亚洲| 亚洲一区二区三区色噜噜| 亚洲18禁久久av| 亚洲av免费高清在线观看| 99热精品在线国产| 身体一侧抽搐| 国产av不卡久久| 少妇人妻精品综合一区二区 | 热99在线观看视频| 一个人观看的视频www高清免费观看| 亚洲美女搞黄在线观看 | 永久网站在线| 长腿黑丝高跟| 啦啦啦啦在线视频资源| 亚洲无线观看免费| 国产精品日韩av在线免费观看| 99久久成人亚洲精品观看| av中文乱码字幕在线| 欧美成人一区二区免费高清观看| 男女那种视频在线观看| 亚洲电影在线观看av| 亚洲熟妇中文字幕五十中出| 久久国内精品自在自线图片| 大又大粗又爽又黄少妇毛片口| 少妇被粗大猛烈的视频| 啦啦啦啦在线视频资源| 日韩欧美国产一区二区入口| 最近在线观看免费完整版| 亚洲av二区三区四区| 久久久精品欧美日韩精品| 久久久精品大字幕| 五月伊人婷婷丁香| 午夜福利在线观看吧| 成人欧美大片| 久久人人精品亚洲av| 国产高清有码在线观看视频| 欧洲精品卡2卡3卡4卡5卡区| 色在线成人网| 麻豆av噜噜一区二区三区| 一进一出抽搐gif免费好疼| 精品久久久久久久末码| 国产大屁股一区二区在线视频| 亚洲第一区二区三区不卡| 亚洲内射少妇av| 国产精品永久免费网站| 欧美日韩精品成人综合77777| 少妇人妻精品综合一区二区 | 观看免费一级毛片| 91麻豆精品激情在线观看国产| 精品午夜福利视频在线观看一区| 精品人妻1区二区| 亚洲欧美日韩高清专用| 久久久久免费精品人妻一区二区| 久久久久久伊人网av| 99精品久久久久人妻精品| 搞女人的毛片| 舔av片在线| 可以在线观看毛片的网站| 免费av毛片视频| 午夜福利成人在线免费观看| av国产免费在线观看| 免费在线观看成人毛片| 国产爱豆传媒在线观看| 校园人妻丝袜中文字幕| xxxwww97欧美| 久久精品国产亚洲av涩爱 | 国产av在哪里看| 久久久久久久精品吃奶| 亚洲欧美日韩高清在线视频| 精品久久久久久久久亚洲 | 日韩人妻高清精品专区| 日本免费a在线| 日韩欧美在线二视频| 1000部很黄的大片| 听说在线观看完整版免费高清| 亚洲av不卡在线观看| 中文亚洲av片在线观看爽| 变态另类丝袜制服| 欧美丝袜亚洲另类 | 最近中文字幕高清免费大全6 | 三级男女做爰猛烈吃奶摸视频| 欧美日韩中文字幕国产精品一区二区三区| 欧美成人免费av一区二区三区| 一级黄色大片毛片| 国产一区二区亚洲精品在线观看| 精品午夜福利视频在线观看一区| 不卡一级毛片| 69av精品久久久久久| 国产色爽女视频免费观看| 999久久久精品免费观看国产| a级毛片免费高清观看在线播放| 亚洲欧美精品综合久久99| 身体一侧抽搐| 午夜亚洲福利在线播放| 欧美极品一区二区三区四区| 亚洲一区二区三区色噜噜| 免费观看人在逋| 精品久久久噜噜| 亚洲精品成人久久久久久| 一本久久中文字幕| 欧美高清成人免费视频www| 少妇人妻一区二区三区视频| 日韩 亚洲 欧美在线| 精品一区二区三区视频在线观看免费| 久久亚洲真实| 国产亚洲精品综合一区在线观看| 精品一区二区免费观看| 国产淫片久久久久久久久| 欧美人与善性xxx| 久久久久久久久大av| 国产精华一区二区三区| 小蜜桃在线观看免费完整版高清| 色噜噜av男人的天堂激情| 可以在线观看的亚洲视频| 我的老师免费观看完整版| 禁无遮挡网站| 99久久久亚洲精品蜜臀av| 国产欧美日韩一区二区精品| 国内少妇人妻偷人精品xxx网站| 麻豆成人午夜福利视频| 成人av一区二区三区在线看| 嫩草影院入口| 亚洲中文日韩欧美视频| 一a级毛片在线观看| 夜夜夜夜夜久久久久| 久久香蕉精品热| 一级毛片久久久久久久久女| 免费观看精品视频网站| 成人国产麻豆网| 国产精品美女特级片免费视频播放器| 人妻制服诱惑在线中文字幕| 国产欧美日韩精品亚洲av| 亚洲avbb在线观看| 久久久久久九九精品二区国产| 亚洲精品久久国产高清桃花| 亚洲性久久影院| 波多野结衣高清作品| 老司机福利观看| 国产精品久久久久久亚洲av鲁大| 国产成人av教育| 91久久精品电影网| 两个人视频免费观看高清| 欧美高清成人免费视频www| 免费av毛片视频| 亚洲成a人片在线一区二区| 国产白丝娇喘喷水9色精品| 久久久久久国产a免费观看| 精品99又大又爽又粗少妇毛片 | 国产欧美日韩精品亚洲av| 老司机福利观看| 男女啪啪激烈高潮av片| 亚洲国产精品成人综合色| 久久久精品大字幕| 亚洲美女视频黄频| 在线看三级毛片| 亚洲国产高清在线一区二区三| 午夜亚洲福利在线播放| 精品国内亚洲2022精品成人| 国产精品av视频在线免费观看| 桃色一区二区三区在线观看| 免费观看在线日韩| 91在线精品国自产拍蜜月| 99热精品在线国产| 久久这里只有精品中国| 麻豆av噜噜一区二区三区| 成人毛片a级毛片在线播放| 欧美成人一区二区免费高清观看| 中文资源天堂在线| 日本a在线网址| 91狼人影院| 国产精品日韩av在线免费观看| 日本在线视频免费播放| 日韩一本色道免费dvd| 性色avwww在线观看| 欧美精品啪啪一区二区三区| 亚洲av成人av| av黄色大香蕉| 国产av在哪里看| 久久久久性生活片| 国产男靠女视频免费网站| 1024手机看黄色片| 变态另类成人亚洲欧美熟女| 我要看日韩黄色一级片| av在线老鸭窝| 久久精品久久久久久噜噜老黄 | 精品久久久噜噜| www.色视频.com| 欧美另类亚洲清纯唯美| 国产一区二区三区av在线 | 我的女老师完整版在线观看| 啦啦啦观看免费观看视频高清| 国产一区二区在线av高清观看| 少妇的逼水好多| 国产精品女同一区二区软件 | 国产v大片淫在线免费观看| 午夜免费成人在线视频| 亚洲成人久久爱视频| 久久精品91蜜桃| 欧美xxxx黑人xx丫x性爽| 欧美精品啪啪一区二区三区| 亚洲va在线va天堂va国产| 日韩中文字幕欧美一区二区| 男插女下体视频免费在线播放| 小说图片视频综合网站| 久久亚洲真实| 伊人久久精品亚洲午夜| 成人午夜高清在线视频| 九九在线视频观看精品| 日韩 亚洲 欧美在线| 无遮挡黄片免费观看| 亚洲久久久久久中文字幕| 日韩欧美三级三区| 天堂网av新在线| 精品乱码久久久久久99久播| 午夜久久久久精精品| 搡女人真爽免费视频火全软件 | 午夜精品一区二区三区免费看| 丰满乱子伦码专区| 黄色配什么色好看| 日本欧美国产在线视频| av在线观看视频网站免费| 国产精品久久电影中文字幕| 国产欧美日韩精品亚洲av| 亚洲国产日韩欧美精品在线观看| 男女视频在线观看网站免费| 国产免费一级a男人的天堂| 在线免费观看不下载黄p国产 | 精品一区二区三区av网在线观看| 九九热线精品视视频播放| 少妇熟女aⅴ在线视频| 国产成人影院久久av| 中文字幕av成人在线电影| 91精品国产九色| 高清日韩中文字幕在线| 99久久精品一区二区三区| 日本爱情动作片www.在线观看 | 亚州av有码| 亚洲三级黄色毛片| 久久精品影院6| 免费在线观看日本一区| 国产男人的电影天堂91| 亚洲在线观看片| 国产精品国产三级国产av玫瑰| 国内久久婷婷六月综合欲色啪| 色噜噜av男人的天堂激情| 国产精品久久久久久久久免| 亚洲乱码一区二区免费版| 啦啦啦啦在线视频资源| 无遮挡黄片免费观看| 非洲黑人性xxxx精品又粗又长| 国产一区二区在线观看日韩| 99久久成人亚洲精品观看| 高清日韩中文字幕在线| 亚洲人成网站高清观看| 美女黄网站色视频| 中文在线观看免费www的网站| 草草在线视频免费看| 久久久久久久久久久丰满 | 给我免费播放毛片高清在线观看| 2021天堂中文幕一二区在线观| 欧美国产日韩亚洲一区| a级毛片免费高清观看在线播放| 日韩精品中文字幕看吧| 真人一进一出gif抽搐免费| 国产高清视频在线播放一区| 一进一出抽搐动态| 成人国产一区最新在线观看| 午夜激情福利司机影院| 亚洲av中文av极速乱 | 午夜福利在线在线| 亚洲欧美激情综合另类| 91av网一区二区| 亚洲图色成人| 免费av毛片视频| 99久久成人亚洲精品观看| 国产精品,欧美在线| 午夜久久久久精精品| 琪琪午夜伦伦电影理论片6080| 少妇高潮的动态图| 国产真实乱freesex| 婷婷色综合大香蕉| 99热精品在线国产| 亚洲aⅴ乱码一区二区在线播放| 悠悠久久av| 欧美成人免费av一区二区三区| 亚洲国产欧洲综合997久久,| 天天一区二区日本电影三级| 午夜老司机福利剧场| 日本-黄色视频高清免费观看| 别揉我奶头~嗯~啊~动态视频| 又粗又爽又猛毛片免费看| 亚洲人成伊人成综合网2020| 国产三级在线视频| 级片在线观看| 偷拍熟女少妇极品色| 亚洲中文字幕日韩| 精品一区二区三区视频在线| 久久久午夜欧美精品| 中出人妻视频一区二区| 热99在线观看视频| 久久国产精品人妻蜜桃| 国产av不卡久久| 免费av观看视频| 美女xxoo啪啪120秒动态图| 欧美+日韩+精品| а√天堂www在线а√下载| 亚洲av五月六月丁香网| 午夜福利在线观看吧| 日本-黄色视频高清免费观看| 欧美色视频一区免费| 嫩草影视91久久| 欧美高清性xxxxhd video| a级毛片免费高清观看在线播放| 免费av不卡在线播放| 性色avwww在线观看| 12—13女人毛片做爰片一| 热99re8久久精品国产| 久久人人精品亚洲av| 老司机福利观看| 香蕉av资源在线| 精品久久久噜噜| 色综合站精品国产| 日韩欧美在线乱码| 69人妻影院| 国产人妻一区二区三区在| 人人妻,人人澡人人爽秒播| 亚洲天堂国产精品一区在线| 偷拍熟女少妇极品色| 三级国产精品欧美在线观看| 特大巨黑吊av在线直播| 琪琪午夜伦伦电影理论片6080| av天堂中文字幕网| 国产亚洲欧美98| 十八禁网站免费在线| 白带黄色成豆腐渣| 成人性生交大片免费视频hd| 成人国产麻豆网| netflix在线观看网站| 禁无遮挡网站| 精品久久久久久,| 高清在线国产一区| 99在线人妻在线中文字幕| 日日摸夜夜添夜夜添av毛片 | 两性午夜刺激爽爽歪歪视频在线观看| 啪啪无遮挡十八禁网站| 亚洲人成网站在线播放欧美日韩| 成人av在线播放网站| 国产精品无大码| 别揉我奶头 嗯啊视频| 97超视频在线观看视频| 久久久久久久久中文| 欧美激情在线99| 国产亚洲精品久久久久久毛片| 午夜爱爱视频在线播放| 自拍偷自拍亚洲精品老妇| av在线老鸭窝| 丝袜美腿在线中文| av在线观看视频网站免费| 日韩大尺度精品在线看网址| 欧美黑人巨大hd| 自拍偷自拍亚洲精品老妇| 免费在线观看日本一区| 欧美xxxx黑人xx丫x性爽| 亚洲 国产 在线| 91麻豆av在线| 美女高潮喷水抽搐中文字幕| 国产精品人妻久久久久久| а√天堂www在线а√下载| 亚洲中文字幕日韩| 亚洲人成伊人成综合网2020| 在线观看一区二区三区| 国产综合懂色|