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

    簡單地球模型的研究進(jìn)展

    2014-03-02 02:47:38王勇劉蘇峽邵亞平彭公炳
    關(guān)鍵詞:復(fù)雜程度概念性氣候

    王勇劉蘇峽邵亞平彭公炳,4

    (1 中國科學(xué)院地理科學(xué)與資源研究所陸地水循環(huán)及地表過程重點(diǎn)實(shí)驗(yàn)室,北京 100101;2 南京信息工程大學(xué)遙感學(xué)院,南京 210044;3 德國科隆大學(xué)地球物理與氣象研究所,科隆 50923;4 國際歐亞科學(xué)院中國科學(xué)中心,北京 100107)

    簡單地球模型的研究進(jìn)展

    王勇1,2劉蘇峽1邵亞平3彭公炳1,4

    (1 中國科學(xué)院地理科學(xué)與資源研究所陸地水循環(huán)及地表過程重點(diǎn)實(shí)驗(yàn)室,北京 100101;2 南京信息工程大學(xué)遙感學(xué)院,南京 210044;3 德國科隆大學(xué)地球物理與氣象研究所,科隆 50923;4 國際歐亞科學(xué)院中國科學(xué)中心,北京 100107)

    以簡單概念性氣候模式(SCMs)和中等復(fù)雜程度的地球系統(tǒng)模式(EMICs)為對象,介紹了簡單地球模型的研究進(jìn)展。在簡單地球模型中,簡單概念性氣候模式已經(jīng)成為了解氣候理論和氣候變化規(guī)律的有效工具;中等復(fù)雜程度的地球系統(tǒng)模式涵蓋了地球系統(tǒng)的大多數(shù)組成部分,對各部分之間的反饋和過程描述比較詳細(xì),同時(shí)計(jì)算要求比較低,使其應(yīng)用幾乎覆蓋了簡單概念性氣候模式以及耦合氣候系統(tǒng)模式的所有研究領(lǐng)域。EMICs必定會有更大的發(fā)展,從而在地球氣候模擬中起到更加重要的作用。

    簡單地球模型,簡單概念性氣候模式,地球系統(tǒng)模式(EMICs)

    1 引言

    當(dāng)今的地球,氣候變化無常,這已為世人矚目。人們一直不懈地從不同角度去探求其原因,尋找合適的應(yīng)對措施。在所有的歸因研究方法中,氣候模式和地球系統(tǒng)模式等(這里統(tǒng)稱為地球模型)是氣候變化歸因研究的重要工具。

    地球模型的共性是在地球系統(tǒng)中的動力、物理、化學(xué)和生物過程基礎(chǔ)上建立起來的數(shù)學(xué)方程組(包括動力學(xué)方程組和參數(shù)化方案)來確定其各個(gè)部分(大氣圈、水圈、冰雪圈、巖石圈、生物圈)的性狀,由此構(gòu)成地球系統(tǒng)的數(shù)學(xué)物理模型,然后用數(shù)值的方法進(jìn)行求解,編制成一種大型綜合性計(jì)算程序,并通過計(jì)算機(jī)付諸實(shí)現(xiàn)對地球系統(tǒng)復(fù)雜行為和過程的模擬與預(yù)測的科學(xué)工具[1-2]。

    20世紀(jì)60年代末,隨著計(jì)算機(jī)技術(shù)的發(fā)展,復(fù)雜耦合氣候系統(tǒng)模式出現(xiàn)了,這類模式包括了發(fā)展成熟的大氣模式、海洋模式、陸面模式、海冰模式等模塊,用以研究包括海洋狀況、冰雪過程、土壤溫濕等在內(nèi)的氣候系統(tǒng)變化規(guī)律,是目前研究大氣、海洋及陸地之間復(fù)雜相互作用的主要工具[1,3]。這類模式基本出發(fā)點(diǎn)在于完善地描述氣候系統(tǒng)的各種過程,模型極其復(fù)雜,計(jì)算時(shí)間和花費(fèi)比較大。同時(shí),根據(jù)研究的實(shí)際需要,陸續(xù)產(chǎn)生了大量不同種類的模式,主要有概念模型、能量平衡模式(EBMs)、輻射對流模式(RCMs)、統(tǒng)計(jì)動力模式(SDMs)和地球系統(tǒng)模式(EMICs)等。目前,依據(jù)復(fù)雜程度的不同,廣

    泛得到應(yīng)用的地球模式可分為3類[4-5]:(1)簡單概念性氣候模式(SCMs);(2)耦合氣候系統(tǒng)模式(CGCMs);(3)復(fù)雜程度介于前2類之間的中等復(fù)雜程度的地球系統(tǒng)模式(EMICs)。

    相對于CGCMs來說,概念模型、能量平衡模式、輻射對流模式和統(tǒng)計(jì)動力模式屬于SCMs,這些模式對氣候系統(tǒng)高度整合,強(qiáng)調(diào)氣候系統(tǒng)的主要?dú)夂蛱卣骱臀锢磉^程,空間分辨率較低,對計(jì)算機(jī)要求不高,能進(jìn)行長尺度的模擬,故在地球模擬研究中得到廣泛應(yīng)用,在研究氣候變化成因上發(fā)揮了重要的作用。在20世紀(jì)90年代末期,中等復(fù)雜程度的EMICs發(fā)展起來的,涵蓋了CGCMs中的大多數(shù)過程,不過采用相對CGCMs更加簡化的方法來表述,這里把這兩類模式統(tǒng)稱為簡單地球模型。本文以簡單概念性氣候模式、中等復(fù)雜程度的地球系統(tǒng)模式為對象,介紹它們的應(yīng)用以及存在的不足。

    2 簡單概念性氣候模式簡介

    簡單概念性氣候模式中的能量平衡模式最先由前蘇聯(lián)的Budyko提出,美國的Sellers在1969年發(fā)表了相類似的一維能量平衡模式,計(jì)算出的氣候要素隨緯度的分布與實(shí)測值吻合得很好。20世紀(jì)70年代是動力氣候研究的春天,各種類型的模式相繼出現(xiàn),如輻射對流模式、統(tǒng)計(jì)動力模式、概念模型等。下文將對部分簡單概念性氣候模式進(jìn)行詳細(xì)介紹。

    2.1 能量平衡模式

    2.1.1 零維能量平衡模式

    公式(1)未考慮大氣中的輻射過程,如果考慮大氣成分(CO2和水汽等)不僅吸收輻射能,而且發(fā)射長波輻射,使大氣氣溫增加,地球表面的氣溫Ts應(yīng)該比有效輻射溫度高。即:

    如果進(jìn)入大氣的輻射和射出大氣的輻射量不平衡,那么地球氣溫將發(fā)生變化,變化率?Ts/?t可以表示為:

    式中,c為地球的比熱容,τa為大氣紅外傳輸率,ε為比例參數(shù)。

    2.1.2 一維能量平衡模式

    一維模式不再將地球視為一點(diǎn),引入地球氣候系統(tǒng)的主要差異,即緯度變化特征,某一緯度帶(φ)能量平衡公式為:

    式中,F(Tφ)是附加項(xiàng),為射出大氣的輻射量,Sφ和Tφ都隨緯度變化。顯然,一維能量平衡模式是零維能量平衡模式的擴(kuò)展。

    將模式中附件項(xiàng)和射出大氣的輻射量經(jīng)過近似和簡化,得到如下公式:

    再利用逐次迭代法,計(jì)算出各個(gè)緯度帶的氣溫。

    除了以上兩種模式外,還可以引入二維(水平或者經(jīng)向/垂直)的能量平衡模式,另外,還有盒式能量平衡模式(Box-EBM),Box-EBM將大氣作為一層,陸面和海洋混合層為一層,海洋中間層和深海為另外兩層,和其他EBM模式相比,Box-EBM強(qiáng)調(diào)深海內(nèi)部的熱量過程以及不同深度層上的能量交換。

    2.2 輻射對流模式

    輻射對流模式[7](RCM)把大氣簡化為一個(gè)鉛直的大氣柱,詳細(xì)考慮大氣柱內(nèi)的輻射過程,通過對流調(diào)整來獲取大氣的垂直溫度分布,是一種考慮時(shí)間變化的一維模式,此模式建立在以下兩個(gè)原理上:(1)在任何高度上的太陽輻射和長波輻射通量與對流熱通量保持平衡;(2)因輻射差異引起的溫度垂直分布的不穩(wěn)定由對流調(diào)整而達(dá)到平衡。

    2.2.1 無對流調(diào)整的輻射平衡模式

    假定大氣氣溫的垂直分布是由輻射收支決定的,則溫度變化方程為:

    式中,右邊第一項(xiàng)、第二項(xiàng)分別代表長波輻射和短波輻射引起的溫度變化。

    2.2.2 有對流調(diào)整的輻射平衡模式

    基本假定:(a)在大氣層頂,凈入射短波輻射等于射出長波輻射;(b)大氣的凈輻射冷卻作用等于大氣長波輻射與短波輻射之差;(c)溫度垂直遞減率小于規(guī)定值時(shí),氣層維持局地輻射平衡;(d)當(dāng)氣溫直減率大于規(guī)定值ε時(shí),對氣溫分布進(jìn)行調(diào)整,以達(dá)到規(guī)定值。

    利用雙極晶體管在中子輻照下的硅位移損傷特性,成功研制了晶體管直流增益在線測試系統(tǒng)。該系統(tǒng)以LabVIEW為開發(fā)平臺,通過多任務(wù)并行處理的模塊化程序設(shè)計(jì),實(shí)現(xiàn)了快速響應(yīng)、高效率同步采集、實(shí)時(shí)數(shù)據(jù)處理、存儲及顯示;采用回讀技術(shù)建立了電壓補(bǔ)償方法,解決了遠(yuǎn)程監(jiān)測中因電壓損耗帶來的測量誤差,系統(tǒng)的測試精度高達(dá)0.2%。

    模式方程為:

    式中,右邊第二項(xiàng)為對流調(diào)整項(xiàng),當(dāng)不考慮對流調(diào)整時(shí)該項(xiàng)為零。在對流層質(zhì)量守恒條件下有:

    在地表面應(yīng)滿足熱量平衡方程,即:

    同時(shí)可定義各層的短波輻射加熱率和長波輻射冷卻率分別為:

    2.3 統(tǒng)計(jì)動力模式

    將大氣沿著緯度圈進(jìn)行平均,用緯度和高度組成的網(wǎng)格點(diǎn)來表示大氣,就構(gòu)成了緯向平均動力模式(ZADM),該模式包括基本的動力和物理過程,是介于一維大氣模式和三維大氣模式之間,連接一維和三維模式的橋梁,因此這類模式在氣候模擬的研究中起著重要的作用。緯向平均模式的主要困難是渦旋輸送的參數(shù)化。由于渦旋輸送的處理是建立在統(tǒng)計(jì)近似的基礎(chǔ)上,所以這類模式也常稱為統(tǒng)計(jì)動力模式[6]。

    把渦旋動量通量寫成:

    式中,KT和Km分別為熱量和動量傳輸系數(shù),可以視為常數(shù)或?yàn)闇囟忍荻鹊暮瘮?shù),緯向平均量用表示。渦旋輸送的重要性實(shí)際上在于通過渦旋熱通量和渦旋動量通量將大氣積蓄的斜壓和正壓不能穩(wěn)定能量釋放出來,這必然對大氣環(huán)流和氣候變化發(fā)生重要的作用。

    2.4 各模式之間的關(guān)系

    通過圖1來表示簡化氣候模式間及其與CGCMs之間的關(guān)系,圖中三角錐的3條邊代表模式的3個(gè)基本要素,即輻射、地面過程和動力學(xué)。垂直方向表示模式的復(fù)雜程度,越接近椎體頂端表示復(fù)雜程度越大。當(dāng)?shù)竭_(dá)頂端時(shí),模式已有了足夠的分辨率,所有因素都被正確而合理地考慮到模式里去了,達(dá)到一個(gè)逼真的氣候模式。在椎體的底部是最簡單的氣候模式,它僅僅考慮一個(gè)主要的過程。同時(shí),椎體底部被認(rèn)為是空心的,因?yàn)樵谧刁w的各邊之間基本上沒有相互作用。沿著金字塔越往上,各主要過程之間的相互作用越大。

    3 簡單概念性氣候模式的應(yīng)用現(xiàn)狀

    圖1 氣侯模式的三角錐[8](EBM:能量平衡模式;RCM:對流輻射模式;SDM:統(tǒng)計(jì)動力模式;CGCMs:耦合氣候系統(tǒng)模式)

    簡單概念性氣候模式在氣候模擬研究中應(yīng)用廣泛,在研究氣候變化成因上發(fā)揮了重要的作用。從第一、二次IPCC報(bào)告十幾個(gè)氣候模式參加評估,到第三次報(bào)告34個(gè)模式被引用,再到第四次和第五次評估報(bào)告更多模式被采用,IPCC歷次報(bào)告所引用模式的數(shù)量在不斷增加,其中,簡單氣候模式(SCM)一直未缺席(IPCC第四次和第五次報(bào)告采用EBMs簡單氣候模式的結(jié)果)。利用簡單氣候模式能夠估測由于熱膨脹導(dǎo)致的全球平均溫度的變化和海平面的上升,此類模式以一套全球或者半球框的形式代表海洋—大氣系統(tǒng),進(jìn)而利用能量平衡方程預(yù)測全球地表溫度,能夠與地球生物化學(xué)圈簡化模式耦合,而可以快速估測氣候?qū)Ω鞣N氣候情景的響應(yīng)[9-10]。Gilliland[11]用零維能量平衡模式,考慮在太陽輻射變化、溫室氣體以及火山氣溶膠強(qiáng)迫情況下,得出太陽活動變化、CO2變化可以解釋北半球氣候變化中低頻變化的93%左右。Schlesinger等[12]也用類似的模式模擬得到溫室氣體、太陽活動和平流氣溶膠三者可以解釋1886—1985年全球氣溫變化的60%。畢鳴等[13]用海氣耦合的能量平衡模式研究了130年來全球氣溫的變化,考慮溫室氣體、太陽活動、火山活動及ENSO的情況下,模擬結(jié)

    果與觀測序列的相關(guān)系數(shù)達(dá)0.884,可解釋77.8%的溫度變化。Crowley[14]利用EBM開展的研究指出,太陽輻射和火山活動可以解釋工業(yè)化以前北半球氣溫低頻變化的41%~64%,20世紀(jì)變暖主要來自溫室氣體的增加,自然變率的貢獻(xiàn)只占25%。Zhao等[15]用一維能量平衡模式研究下墊面改變對降水和溫度造成的影響,取得一定的成功。王紹武等[16]利用簡化的能量平衡模式,綜合考慮太陽、火山、溫室氣體的輻射效應(yīng),發(fā)現(xiàn)EBM能夠較好模擬近千年的全球平均氣溫序列,能夠成功再現(xiàn)中世紀(jì)暖期、小冰期和20世紀(jì)變暖等特征。Bertrand等[17]利用二維緯向平均模式所開展的研究,同樣表明太陽和火山活動可以解釋小冰期和中世紀(jì)暖期,而20世紀(jì)后期的加速變暖,則只有溫室氣體可以解釋;1850年以后氣溫的合理模擬,還需要考慮人為硫酸鹽氣溶膠強(qiáng)迫和森林采伐的影響。劉玉芝等[18]利用一維輻射對流氣候模式,詳細(xì)研究了云量、云光學(xué)厚度以及云高等要素的變化對大氣頂和地面太陽短波輻射和紅外長波輻射通量以及云的輻射強(qiáng)迫的影響。丁仲禮等[19]就晚更新世東亞古季風(fēng)變化方面建立起一個(gè)概念模型,模型考慮了太陽輻射變化、中高緯、青藏高原、邊緣海以及熱帶海洋等因素,得出東亞古季風(fēng)主要由全球冰量通過極地冰蓋擴(kuò)張、邊緣海變化以及青藏高原冰川作用3個(gè)過程來影響的,太陽輻射變化的直接控制作用相對較小,較合理地解釋了黃土所記錄的古季風(fēng)變化特征。

    從上述模擬來看,簡單概念性氣候模式涵蓋了主要?dú)夂蛱卣骱椭饕锢磉^程,如熱力動力過程、輻射對流過程或關(guān)鍵反饋過程等,同時(shí)模式對計(jì)算機(jī)要求很低,能夠進(jìn)行長時(shí)間的模擬試驗(yàn)。不過,簡單概念性氣候模式通常采用緯向平均,最多二維,只能反映大區(qū)域的一個(gè)平均過程,空間分辨率不高。隨著當(dāng)前氣候變化動力機(jī)制逐漸明確、考慮更多的物理過程的情況下,這種模式存在一定的局限性,體現(xiàn)在難以給出氣候變化的高分辨率空間特征,同時(shí)還不包含氣候系統(tǒng)內(nèi)部變率過程。如果僅考慮“簡單”背景下的研究,其研究結(jié)果具有一定的借鑒意義。

    4 中等復(fù)雜程度的地球系統(tǒng)模式簡介

    中等復(fù)雜程度的地球系統(tǒng)模式(EMICs)于20世紀(jì)90年代末期建立,它一般具有二維(兩個(gè)水平維向或垂直和一個(gè)水平維向)空間分辨率,能夠描述包括人類和自然相互作用的自然地球系統(tǒng),EMICs涵蓋了CGCMs中的大多數(shù)過程,不過采用相對CGCMs更加簡化的方法去表述,與CGCMs模擬的分辨率也有差別。EMICs的復(fù)雜程度介于簡單概念性氣候模式和CGCMs模式之間,整合化程度沒有簡單概念性氣候模式高,但比CGCMs要高,所包含的過程和描述較CGCMs更加簡化,過程比簡單概念性氣候模式要多,細(xì)節(jié)描述也較之更加豐富,兼顧了簡單概念性氣候模式和CGCMs的部分優(yōu)點(diǎn)(圖2)。不難看出,相對于CGCMs模式,EMICs也是一種簡化氣候模式。

    EMICs中的大氣模塊主要采用的是能量水汽平衡模式、準(zhǔn)地轉(zhuǎn)模式、統(tǒng)計(jì)動力學(xué)模式和大氣環(huán)流模式。大氣模塊的分辨率通常比較低,大多采用緯向平均和能量水氣平衡模式。由表1可見,EMICs在大氣模塊構(gòu)成和簡化程度方面存在明顯的差異,根據(jù)這種差異,模式可以分為兩類:一類是基于CGCMs進(jìn)行大氣或者海洋過程的簡化,使其總體的復(fù)雜程度較高,對氣候過程及反饋描述比較完善,對氣候過程描述比較細(xì)膩,如德國漢堡大學(xué)Planet Simulator模式[21-22];而另一類基于低分辨率且大量參數(shù)化的模式,這類EMICs在更長時(shí)間尺度(千年范圍乃至冰期間冰期循環(huán))氣候變化的研究中有著重要作用,如CLIMBER-2的大氣模塊是二維統(tǒng)計(jì)動力學(xué)模式。從模式的復(fù)雜程度上看,前者更像是簡化了的CGCMs,而后者則更接近于簡單概念性氣候模式。綜述所述,從整體上簡化了EMICs的復(fù)雜程度,降低了對計(jì)算能力的要求,有助于在古氣候資料的強(qiáng)迫下模擬出不同于今天的長尺度古氣候狀態(tài),充分顯示出了EMICs模式的優(yōu)點(diǎn);但EMICs的分辨率比較粗糙,對地球系統(tǒng)細(xì)節(jié)的描述也相當(dāng)簡化,因而模擬無法達(dá)到CGCMs的高度。

    5 中等復(fù)雜程度的地球系統(tǒng)模式的應(yīng)用現(xiàn)狀

    當(dāng)前EMICs在氣候研究中得到廣泛的應(yīng)用,雖然其對氣候系統(tǒng)的描繪不如CGCMs細(xì)致,但是EMICs具有與簡單概念性模式相當(dāng)?shù)挠?jì)算效率,所以EMICs主要注重長期氣候模擬,以足夠的細(xì)節(jié)來刻畫地球系統(tǒng)的長期演化。

    圖2 EMICs 的圖解定義[20]

    表1 主要EMICs及相關(guān)介紹[5]

    目前,國際古氣候模擬比較計(jì)劃(PMIP)已經(jīng)歷了PMIP1(1991—2001年)和PMIP2(2002—2007年)兩個(gè)階段,并正在開展第三階段PMIP3(2008年至今)的研究。PMIP3階段致力于利用地球系統(tǒng)模式對末次盛冰期和中全新世(距今最近的、最有代表性的寒冷和溫暖的氣候狀態(tài))的氣候進(jìn)行模擬。在古氣候模擬的研究方面,產(chǎn)生了不少研究成果。Gallcee等[23]和Berger等[24]的研究表明,氣候模式LLN-2D在變化的太陽輻射強(qiáng)迫下,當(dāng)指定的常量CO2濃度低于220×10-6后,能夠重現(xiàn)觀測到的冰期—間冰期循環(huán);而固定軌道參數(shù),讓CO2依據(jù)Vostok重建記錄變化,則LLN-2D無法再現(xiàn)冰期—間冰期循環(huán)[25]。EMICs對冰期的開始和結(jié)束也進(jìn)行了各種試驗(yàn),Yoshimori等[26]利用Uvic模式進(jìn)行了末次冰期結(jié)束對軌道參數(shù)和CO2強(qiáng)迫的敏感性試驗(yàn),Meissner等[27]使用UVic研究了植被對末次冰期開始的作用,Kageyama等[28]使用CLIMBER-2探討了末次冰期的開始特別是各種不同冰川的反饋?zhàn)饔?。Wang等[29]使用MPM-2對冰期氣候突變和D-O震蕩的機(jī)制進(jìn)行了深入探討,發(fā)現(xiàn)北大西洋經(jīng)向翻轉(zhuǎn)流(AMOC)狀態(tài)的快速轉(zhuǎn)換導(dǎo)致冰期氣候突變,這一突變被海冰面積的擴(kuò)張進(jìn)一步強(qiáng)化。

    另外,對歷史、當(dāng)今、未來氣候模擬的應(yīng)用研究,在IPCC AR4中,EMIC已較前幾次評估報(bào)告被更多人使用并得到深入的評估,在海平面上升歸因分析上,利用一個(gè)低分辨率CGCM 和數(shù)個(gè)以氣候敏感性最佳估計(jì)值3℃為基礎(chǔ)的EMICs進(jìn)行推估,得出海平面上升的促成因素是海洋熱膨脹,并對未來進(jìn)行了預(yù)估[9]。在2013年出版的IPCC AR5中,EMICs在古氣候模擬、未來氣候預(yù)估等方面發(fā)揮了重要的作用[10]。Loutre等[30]使用LLN-2D對下次冰期的開始也進(jìn)行了探討,研究了未來的氣候及其對不同CO2情景的敏感性,發(fā)現(xiàn)目前的氣候很可能經(jīng)歷較長的間冰期(大約5萬年),而在未來的6萬年和10萬年將可能出現(xiàn)較小和較大的冰期。Bertrand等[31]和Bauer等[32]分別使用MoBiDiC和CLIMBER-2對過去1000年的各種氣候強(qiáng)迫年進(jìn)行了評估,與高分辨率的重建數(shù)據(jù)進(jìn)行比較,發(fā)現(xiàn)火山和太陽輻射作為自然強(qiáng)迫,是工業(yè)革命時(shí)期以前的主要?dú)夂驈?qiáng)迫,能夠很好地解釋中世紀(jì)暖期和小冰期。但自然強(qiáng)迫不足以產(chǎn)生近代以來的全球變暖,人類活動排放的溫室氣體是造成全球變暖的主要強(qiáng)迫,而且研究結(jié)果表明所有的氣候強(qiáng)迫對模擬過去千年氣候變化都是重要的。Petoukhov等[33]對目前氣候和CO2加倍模擬的EMICs進(jìn)行比較,充分顯示EMICs對目前氣候的模擬非常接近觀測結(jié)果,而且EMICs模擬誤差均在CGCMs誤差的范圍之內(nèi)。Wang等[34]和Marsh等[35]分別使用MPM-2和GENIE對溫鹽環(huán)流進(jìn)行研究,發(fā)現(xiàn)溫鹽環(huán)流這個(gè)洋流系統(tǒng)最少有兩種穩(wěn)定的氣候狀態(tài),其中一種是目前所出現(xiàn)的,即環(huán)流非常強(qiáng)勁和向極地方向傳輸大量熱能,從而帶來溫暖氣候的狀態(tài);另一種則是環(huán)流非常微弱,甚至接近于停止的狀態(tài),從而帶來寒冷的氣候狀態(tài)。數(shù)值模擬試驗(yàn)還顯示,這兩種狀態(tài)之間是可能出現(xiàn)相互轉(zhuǎn)化的。另外,利用EMICs對陸地、海洋、全球碳循環(huán)等多個(gè)方面的問題也進(jìn)行了探討。

    結(jié)合當(dāng)前EMICs的應(yīng)用和發(fā)展,可見EMICs應(yīng)用范圍很廣,模擬能力很強(qiáng),未來如在下面幾個(gè)方面取得突破,必定會有更大的發(fā)展,從而在地球系統(tǒng)模擬中起到更加重要的作用:(1)加入人類活動、考慮固體地球過程和太陽活動的影響,使EMICs成為真正意義上的地球模式;(2)注重氣候突變機(jī)理的研究,區(qū)分自然和人為因素在突變中的相對貢獻(xiàn);(3)改進(jìn)模式算法,提高硬件設(shè)備的運(yùn)算能力,使其耗費(fèi)更少的計(jì)算時(shí)間,能夠進(jìn)行更長時(shí)間尺度上的模擬研究。

    6 結(jié)論

    本文介紹了以簡單概念性氣候模式(SCMs)和中等復(fù)雜程度的地球系統(tǒng)模式(EMICs)為代表的簡單地球模型,它們在氣候模擬方面都有不凡表現(xiàn)。簡單概念性氣候模式將對氣候變化的研究拓展到了古氣

    候領(lǐng)域,已經(jīng)成為了解氣候理論和氣候變化規(guī)律的有效工具。EMICs涵蓋了地球系統(tǒng)的大多數(shù)組成部分,對各部分之間的反饋和過程描述比較詳細(xì),同時(shí)計(jì)算要求比較低,使其應(yīng)用幾乎覆蓋了SCMs以及耦合氣候系統(tǒng)模式(CGCMs)的所有研究領(lǐng)域。EMICs必定會有更大的發(fā)展,從而在地球氣候模擬中起到更加重要的作用。

    EMICs也有其自身的缺點(diǎn),分辨率比較低,對地球系統(tǒng)細(xì)節(jié)的描述也相當(dāng)?shù)暮喕?尚不能達(dá)到CGCMs的高度。但是,EMICs能夠刻畫出氣候變化的大體特征,更能進(jìn)行CGCMs在目前計(jì)算條件下無法比擬的長期氣候變化模擬。2005年,德國漢堡大學(xué)完成了Planet Simulator模式,此模式屬于典型的EMICs模式,由全球大氣環(huán)流模式、海洋—海冰模塊和陸地土壤—生物圈模塊組成,以并行計(jì)算機(jī)為硬件支撐,支持地球和類地球行星的氣候動力學(xué)數(shù)值試驗(yàn),同時(shí)還具備結(jié)構(gòu)簡單、可視化界面和開放源代碼便于改造等特點(diǎn),依托此模式使改進(jìn)模式、模擬長尺度的古氣候和開展敏感性試驗(yàn)成為了可能。

    EMICs從出現(xiàn)到發(fā)展到今天不過短短的十幾年,尚處于發(fā)展的初始階段,隨著計(jì)算機(jī)硬件水平的提高、模式的不斷改進(jìn)和對氣候系統(tǒng)認(rèn)識的不斷加深,在未來的幾十年或者更長的時(shí)期內(nèi),對氣候時(shí)間尺度(千年甚至更長)的氣候變化的研究,EMICs很可能依然是比較理想的選擇。今后一段時(shí)間,如何利用EMICs就一些重要驅(qū)動因子(如天文和地球因子)進(jìn)行敏感性試驗(yàn),如何采用混沌理論研究其非線性問題[36],如何改進(jìn)對重點(diǎn)區(qū)域和重點(diǎn)過程的模擬,將是氣候數(shù)值模擬領(lǐng)域面臨的重要課題。

    [1]王斌, 周天軍, 俞永強(qiáng), 等. 地球系統(tǒng)模式發(fā)展展望. 氣象學(xué)報(bào), 2008, 66(6): 857-869.

    [2]《10000個(gè)科學(xué)難題》地球科學(xué)編委會. 10000個(gè)科學(xué)難題·地球科學(xué)卷. 北京: 科學(xué)出版社, 2010.

    [3]周天軍, 王在志, 宇如聰, 等. 基于LASG/IAP大氣環(huán)流譜模式的氣候系統(tǒng)模式. 氣象學(xué)報(bào), 2005, 63(5): 702-715.

    [4]Houghton J T, Ding Y, Griggs D J, et al. Climate change 2001: the scientific basis. United Kingdom and New York: Cambridge University Press, 2001.

    [5]張冉, 李力, 郭慶春, 等. 古氣候研究中氣候模式的發(fā)展與應(yīng)用. 干旱區(qū)研究, 2007, 24(5): 704-711.

    [6]李崇銀. 氣候動力學(xué)引論(第二版). 北京: 氣象出版社, 2000.

    [7]于革. 古氣候動力模擬. 北京: 高等教育出版社, 2007.

    [8]曹鴻興. 氣候模式與氣候模擬. 氣象, 1987, 13(11): 3-6.

    [9]IPCC. IPCC fourth assessment report (AR4). Cambridge: Cambridge University Press, 2007.

    [10]IPCC. IPCC fifth assessment report (AR5). Cambridge: Cambridge University Press, 2013.

    [11]Gilliland R L. Solar, volcanic and CO2forcing of recent climatic changes. Climate Change, 1982, 4: 11-13.

    [12]Schlesinger M E, Ramalakutty N. Implication for globalwarming of intercycle solar irradiance variations. Nature, 1992, 360: 330-333.

    [13]畢鳴, 王紹武, 李曉東. 130年來全球氣溫變化的模擬. 第四紀(jì)研究, 1997(1): 21-27.

    [14]Crowley T J. Causes of climate change over the past 1000 years. Science, 2000, 289: 270-277.

    [15]Zhao M, Zeng X. A theoretical analysis on the local climate change induced by the change of landuse. Adv Atm Sci, 2002, 19: 45-63.

    [16]王紹武, 謝志輝, 蔡靜寧, 等. 近千年全球平均氣溫變化的研究.自然科學(xué)進(jìn)展, 2002, 12(11): 1145-1149.

    [17]Bertrand C, Loutre M F, Crucifix M, et al. Climate of the last millennium: a sensitivity study. Tellus, 2002, 54A: 221-244.

    [18]劉玉芝, 石廣玉, 趙劍琦. 一維輻射對流模式對云—輻射強(qiáng)迫的數(shù)值模擬研究.大氣科學(xué), 2007, 21(3): 486-493.

    [19]丁仲禮, 劉東生. 晚更新世東亞古季風(fēng)變化動力機(jī)制的概念模型. 科學(xué)通報(bào), 1998, 43(2): 122-132.

    [20]Claussen M. Earth system models //Ehlers E, Krafit T, eds. Understanding the Earth system: compartments, processes and interactions. Heidelberg, Berlin, New York: Springer, 2000.

    [21]Fraedrich K, Jansen H, Kirk E, et al. The planet simulator: Green planet and desert world. Meteorologische Zeitschrift, 2005, 14(3): 305-314.

    [22]Fraedrich K, Jansen H, Kirk E, et al. The planet simulator: towards a user friendly model. Meteorologische Zeitschrift, 2005, 14(3): 299-304.

    [23]Gallcee H, Ypersele J P, Fichefet T, et al. Simulation of the last glacial cycle by a coupled, sectorially averaged climate-ice sheet model II. Response to insolation and CO2variation. Journal of Geophysical Research, 1992, 97(15):713-715.

    [24]Berger A, Loutre M F, Gallcee H. Sensitivity of the LLN climate model to the astronomical and CO2forcings over the last 200kyr. Climate Dynamics, 1998, 14: 615-629.

    [25]Loutre M F, Berger A. No glacial-interglacial cycle in the ice volume simulated under a constant astronomical forcing and avariable CO2. Geophysical Research Letters, 2000, 27(6): 783-786.

    [26]Yoshimori M, Weaver A J, MarshallS J, et al. Glacial terminations: sensitivity to orbital and CO2forcing in a coupled climate system model. Climate Dynamics, 2001, 17: 571-588.

    [27]Meissner K J, Weaver A J, Matthews H D, et al. The role of land surface dynamics in glacial inception: a study with the Uvic Earth System Model. Climate Dynamics, 2003, 21: 515-537.

    [28]Kageyama M, Charbit S, Ritz C,et al. Quantifying ice-sheet feedbacks during the last glacial inception. Geophysical Research Letters, 2004, 31, L24203, do:i 10. 1029/2004GL021339.

    [29]Wang Z, Mysak L A. Glacial abrupt climate changes and Dansgaard-Oeschger oscillations in a coupled climatemodel. Paleoceanography, 2006, 21: PA2001.

    [30]Loutre M F, Berger A. Future climatic changes: Are we entering an exceptionally long Interglacial. Climatic Change, 2000, 46: 61-90.

    [31]Bertrand C, Loutre M F, Crucifix M, et al. Climate of the last millennium: a sensitivity study. Tellus, 2002, 54A: 221-244.

    [32]Bauer E, Claussen M, Brovkin V. Assessing climate forcings of the Earth system for the past millennium. Geophysical Research Letters, 2003, 30(6): 1276-1280.

    [33]Petoukhov V, Claussen M, Berger A, et al. EMIC intercomparison project (EMIP-CO2): Comparative analysis of EMIC simulations of currentclimate and equilibrium and transient responses to atmospheric CO2doubling. Climate Dynamics, 2005, 25: 363-385.

    [34]Wang Z. Two climatic states and feedbacks on thermohaline circulation in an Earth system model of intermediate complexity. Climate Dynamics, 2005, 25: 299-314.

    [35]Marsh R, Yool A, Lenton T M, et al. Bistability of the thermohaline circulation identified through comprehensive 2-parameter sweeps of an efficient climate model. Climate Dynamics, 2004, 23: 761-777.

    [36]Shao Y, Peng G, Leslie L M. The environmental dynamic system //Peng G, Leslie LM, Shao Y. Environmental modelling and prediction. Berlin: Springer-Verlag, 2002.

    Progress in Simple Earth Models

    Wang Yong1,2, Liu Suxia1, Shao Yaping3, Peng Gongbing1,4
    (1 Key Laboratory of Water Cycle and Related Land Surface Processes, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101; 2 School of Remote Sensing, Nanjing University of Information Science and Technology, Nanijing 210044; 3 Institute for Geophysics and Meteorology, University of Cologne, Cologne 50923; 4 China Science Center, International Eurasian Academy of Sciences, Beijing 100107)

    In this paper, by taking Simple conceptual Climate Models (SCMs) and Earth system Models of Intermediate Complexity (EMICs) as research objects, simple climate models and their application in climate modeling are introduced. Simple conceptual climate models have become effective tools to understand the theory of climate and climate-change law. Earth system Models of Intermediate Complexity (EMICs) cover most of the components in the Earth system, and describe the process and feedback among the various components in more details. Because EMICs have low computational requirements, modeling applications of EMICs cover almost all fi elds of simple conceptual climate models and CGCMs simulations. EMICs will surely achieve a greater development and play a more important role in the Earth's climate simulation.

    sample earth model, simple conceptual climate models, EMICs

    10.3969/j.issn.2095-1973.2014.03.004

    2013年12月17日;

    2014年4月14日

    王勇(1979—),Email:wyong0210@igsnrr.ac.cn

    資助信息:國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目(2012CB957802)

    猜你喜歡
    復(fù)雜程度概念性氣候
    未來主義
    睿士(2020年5期)2020-05-21 09:56:35
    探究高校會計(jì)與財(cái)務(wù)的復(fù)雜性
    瞧,氣候大不同
    氣候變暖會怎樣?
    精準(zhǔn)把握錯因有效分析糾正
    初中幾何教材認(rèn)知復(fù)雜程度的比較研究
    ——以中國、新加坡教材的三角形問題為例
    夜郎古國游樂園劇院概念性方案設(shè)計(jì)
    新一代(2016年15期)2016-11-16 15:41:36
    口孜東煤礦81煤層斷裂復(fù)雜程度定量評價(jià)
    綠色科技(2015年2期)2016-01-16 01:26:27
    立冬
    概念性知識對維、漢語小數(shù)讀法的影響
    語言與翻譯(2014年2期)2014-07-12 15:49:30
    久久久久久久久久久丰满| 亚州av有码| 免费观看在线日韩| 欧美日韩国产亚洲二区| 美女被艹到高潮喷水动态| 国产日韩欧美在线精品| 国产精品久久电影中文字幕| 国产乱来视频区| 插逼视频在线观看| 国产日韩欧美在线精品| 一级毛片aaaaaa免费看小| 亚洲图色成人| 99久国产av精品| 成人av在线播放网站| 国产亚洲av嫩草精品影院| 视频中文字幕在线观看| 青春草亚洲视频在线观看| 嫩草影院精品99| 丰满乱子伦码专区| 久久久久久久午夜电影| 久久久久久久午夜电影| 最近最新中文字幕免费大全7| 最新中文字幕久久久久| 国产不卡一卡二| 亚洲真实伦在线观看| 欧美又色又爽又黄视频| 亚洲精品自拍成人| 亚洲精品aⅴ在线观看| 高清午夜精品一区二区三区| 国产高清不卡午夜福利| 国产精品一区二区在线观看99 | 中国国产av一级| 床上黄色一级片| 97超视频在线观看视频| 日日撸夜夜添| 国产精品一区二区三区四区久久| 99久久中文字幕三级久久日本| 国产一区亚洲一区在线观看| 国产男人的电影天堂91| 亚洲伊人久久精品综合 | 激情 狠狠 欧美| 亚洲精品国产成人久久av| 亚洲国产精品专区欧美| 亚洲av成人精品一二三区| 婷婷色综合大香蕉| 日本午夜av视频| 欧美精品一区二区大全| 亚洲精品久久久久久婷婷小说 | 黑人高潮一二区| 国产精品av视频在线免费观看| 亚洲欧美日韩无卡精品| 身体一侧抽搐| 久久久久久久久久久免费av| 国产毛片a区久久久久| 少妇高潮的动态图| 九九在线视频观看精品| 全区人妻精品视频| 精品一区二区三区人妻视频| 91久久精品电影网| 男的添女的下面高潮视频| 在线观看av片永久免费下载| 日日啪夜夜撸| 日日啪夜夜撸| 久久精品国产自在天天线| 久久亚洲精品不卡| 国产精品一及| www.色视频.com| 亚洲精品色激情综合| 男女边吃奶边做爰视频| 日韩强制内射视频| 免费观看性生交大片5| 亚洲五月天丁香| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品影视一区二区三区av| 亚洲五月天丁香| 一个人观看的视频www高清免费观看| 毛片一级片免费看久久久久| 国产精品人妻久久久影院| av片东京热男人的天堂| 日韩中字成人| 最近最新中文字幕大全免费视频 | 国产精品99久久99久久久不卡 | 日韩一区二区视频免费看| 天天躁夜夜躁狠狠久久av| 国产av码专区亚洲av| 午夜影院在线不卡| 少妇精品久久久久久久| 久久精品国产亚洲av天美| 狂野欧美激情性xxxx在线观看| 男女高潮啪啪啪动态图| 欧美人与性动交α欧美精品济南到 | 免费av不卡在线播放| 亚洲精品美女久久av网站| 日韩成人伦理影院| 中国美白少妇内射xxxbb| 精品人妻熟女毛片av久久网站| 精品人妻熟女毛片av久久网站| 免费看光身美女| 亚洲,欧美,日韩| 99re6热这里在线精品视频| av视频免费观看在线观看| 亚洲性久久影院| 不卡视频在线观看欧美| 久久久久精品人妻al黑| 欧美激情国产日韩精品一区| 亚洲精华国产精华液的使用体验| 国产色婷婷99| 日韩欧美一区视频在线观看| 少妇的逼水好多| 午夜视频国产福利| 亚洲,欧美,日韩| 成人亚洲精品一区在线观看| 国产有黄有色有爽视频| 最近最新中文字幕大全免费视频 | 欧美日韩一区二区视频在线观看视频在线| 香蕉丝袜av| freevideosex欧美| 日本av免费视频播放| 亚洲精品视频女| 18在线观看网站| 亚洲色图 男人天堂 中文字幕 | 国产国语露脸激情在线看| 国产国语露脸激情在线看| 亚洲伊人久久精品综合| 一区二区av电影网| 26uuu在线亚洲综合色| 日本wwww免费看| 亚洲精品乱久久久久久| 日韩av免费高清视频| 中文字幕最新亚洲高清| 美女xxoo啪啪120秒动态图| 免费观看性生交大片5| 人妻人人澡人人爽人人| 亚洲综合精品二区| 日韩免费高清中文字幕av| 一本—道久久a久久精品蜜桃钙片| 亚洲伊人色综图| 国产在视频线精品| 精品一区在线观看国产| 婷婷色综合大香蕉| 国产av精品麻豆| 少妇的逼水好多| 自线自在国产av| 中文字幕精品免费在线观看视频 | 久久精品久久久久久久性| 国产精品熟女久久久久浪| 色视频在线一区二区三区| 日本爱情动作片www.在线观看| 女性被躁到高潮视频| 少妇熟女欧美另类| 三上悠亚av全集在线观看| 日韩电影二区| 国产伦理片在线播放av一区| av.在线天堂| 久久99一区二区三区| 如日韩欧美国产精品一区二区三区| 草草在线视频免费看| 亚洲精品国产色婷婷电影| 久久精品久久精品一区二区三区| 热re99久久国产66热| 国产老妇伦熟女老妇高清| 青春草国产在线视频| 少妇 在线观看| 亚洲精品国产av蜜桃| 午夜免费观看性视频| 国产在线视频一区二区| 欧美xxⅹ黑人| 精品福利永久在线观看| 国产精品成人在线| 国产精品人妻久久久影院| 最近中文字幕高清免费大全6| 亚洲精品日本国产第一区| 另类精品久久| 国产亚洲av片在线观看秒播厂| 成年女人在线观看亚洲视频| 女性生殖器流出的白浆| 99视频精品全部免费 在线| freevideosex欧美| 最近最新中文字幕免费大全7| 亚洲精品国产色婷婷电影| 精品少妇内射三级| 51国产日韩欧美| 18禁动态无遮挡网站| av视频免费观看在线观看| 99精国产麻豆久久婷婷| 中文欧美无线码| 九色成人免费人妻av| 久久午夜福利片| 欧美成人午夜精品| 一个人免费看片子| av有码第一页| 18在线观看网站| 国产精品免费大片| 一二三四中文在线观看免费高清| 欧美日韩av久久| av线在线观看网站| 亚洲精品一二三| 免费女性裸体啪啪无遮挡网站| 日韩精品有码人妻一区| 久久影院123| 一本大道久久a久久精品| 黄色 视频免费看| av视频免费观看在线观看| 国产免费现黄频在线看| 亚洲精品久久久久久婷婷小说| 欧美少妇被猛烈插入视频| 成人手机av| 国产亚洲欧美精品永久| 有码 亚洲区| 久久午夜福利片| 久久精品国产a三级三级三级| 国产毛片在线视频| av网站免费在线观看视频| 成人国产av品久久久| 91久久精品国产一区二区三区| 在现免费观看毛片| av线在线观看网站| 捣出白浆h1v1| 天天操日日干夜夜撸| 亚洲精品乱久久久久久| 国产精品一区二区在线观看99| 久久精品熟女亚洲av麻豆精品| 国产男女超爽视频在线观看| 亚洲欧美日韩另类电影网站| 国产爽快片一区二区三区| 18禁国产床啪视频网站| 亚洲一码二码三码区别大吗| 亚洲成国产人片在线观看| 国精品久久久久久国模美| 国产乱来视频区| 久久久精品免费免费高清| 一级毛片黄色毛片免费观看视频| 女人精品久久久久毛片| 丝袜美足系列| 国产欧美日韩一区二区三区在线| 国产精品不卡视频一区二区| 国产永久视频网站| 国产亚洲欧美精品永久| 大香蕉97超碰在线| 日本91视频免费播放| 大香蕉久久网| 亚洲精品美女久久av网站| 婷婷色麻豆天堂久久| 免费久久久久久久精品成人欧美视频 | 成人国产av品久久久| 亚洲内射少妇av| 咕卡用的链子| 久热久热在线精品观看| a级片在线免费高清观看视频| 男女免费视频国产| 搡女人真爽免费视频火全软件| 日本av手机在线免费观看| 韩国高清视频一区二区三区| 午夜免费观看性视频| av网站免费在线观看视频| 一级毛片 在线播放| 丝袜美足系列| 内地一区二区视频在线| 久久亚洲国产成人精品v| 男女啪啪激烈高潮av片| 啦啦啦在线观看免费高清www| 最近手机中文字幕大全| 丝袜喷水一区| 免费观看在线日韩| 七月丁香在线播放| 国内精品宾馆在线| 亚洲第一av免费看| 成人亚洲精品一区在线观看| 搡女人真爽免费视频火全软件| 久久久久人妻精品一区果冻| 中文字幕精品免费在线观看视频 | 国产麻豆69| 大片免费播放器 马上看| 高清在线视频一区二区三区| 免费高清在线观看日韩| 亚洲国产欧美在线一区| 国产精品嫩草影院av在线观看| videos熟女内射| 肉色欧美久久久久久久蜜桃| 黑人高潮一二区| 午夜福利影视在线免费观看| 18禁观看日本| 日韩,欧美,国产一区二区三区| 看非洲黑人一级黄片| 色视频在线一区二区三区| 晚上一个人看的免费电影| 看免费成人av毛片| 26uuu在线亚洲综合色| 大片免费播放器 马上看| 日韩制服丝袜自拍偷拍| 女性生殖器流出的白浆| 日本wwww免费看| 久久精品国产综合久久久 | 久久久久精品久久久久真实原创| 看免费av毛片| 久热久热在线精品观看| 日韩,欧美,国产一区二区三区| 国产熟女午夜一区二区三区| av国产精品久久久久影院| 亚洲国产欧美日韩在线播放| 99re6热这里在线精品视频| 最黄视频免费看| 一二三四中文在线观看免费高清| 亚洲精品国产av蜜桃| 久久精品国产亚洲av涩爱| 男人操女人黄网站| 日日摸夜夜添夜夜爱| 丝袜在线中文字幕| 久久亚洲国产成人精品v| 亚洲,欧美精品.| av播播在线观看一区| 制服人妻中文乱码| 热re99久久精品国产66热6| 五月开心婷婷网| 少妇人妻 视频| 国产老妇伦熟女老妇高清| 国产伦理片在线播放av一区| 极品人妻少妇av视频| 精品一区二区三区视频在线| 亚洲在久久综合| 国产一区亚洲一区在线观看| 999精品在线视频| 久久99一区二区三区| av国产精品久久久久影院| 日韩视频在线欧美| 波多野结衣一区麻豆| 国产精品嫩草影院av在线观看| 高清视频免费观看一区二区| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品成人av观看孕妇| 在线天堂中文资源库| 又粗又硬又长又爽又黄的视频| 亚洲美女搞黄在线观看| 亚洲av国产av综合av卡| 18禁在线无遮挡免费观看视频| 久久毛片免费看一区二区三区| 美女国产视频在线观看| 激情视频va一区二区三区| 一级毛片黄色毛片免费观看视频| 制服人妻中文乱码| 日本欧美视频一区| 男人操女人黄网站| 在现免费观看毛片| 熟女av电影| 精品国产一区二区久久| 欧美日韩av久久| 久久久久久久久久久免费av| 国产国拍精品亚洲av在线观看| 欧美 亚洲 国产 日韩一| 丝袜喷水一区| 亚洲av欧美aⅴ国产| 在线观看一区二区三区激情| 亚洲av日韩在线播放| 天堂8中文在线网| 久久精品国产亚洲av天美| 男的添女的下面高潮视频| 日本爱情动作片www.在线观看| 国产亚洲av片在线观看秒播厂| freevideosex欧美| 两个人看的免费小视频| 九九爱精品视频在线观看| 男女下面插进去视频免费观看 | 午夜影院在线不卡| av在线观看视频网站免费| 黄色毛片三级朝国网站| 国产免费又黄又爽又色| 在线观看免费日韩欧美大片| 免费观看性生交大片5| 18禁裸乳无遮挡动漫免费视频| 国内精品宾馆在线| 亚洲人成网站在线观看播放| 亚洲精品视频女| 精品久久蜜臀av无| 日韩 亚洲 欧美在线| 国产精品麻豆人妻色哟哟久久| 久久久久久久久久人人人人人人| 大片电影免费在线观看免费| 亚洲五月色婷婷综合| 亚洲av电影在线进入| 九色成人免费人妻av| 国产成人欧美| 国产av码专区亚洲av| 国国产精品蜜臀av免费| 九九爱精品视频在线观看| 纵有疾风起免费观看全集完整版| 亚洲av电影在线观看一区二区三区| 成年av动漫网址| 日韩人妻精品一区2区三区| 中文字幕av电影在线播放| 夜夜爽夜夜爽视频| 亚洲美女黄色视频免费看| 大片免费播放器 马上看| 午夜视频国产福利| 日韩中字成人| 国产免费一区二区三区四区乱码| 精品国产一区二区三区四区第35| 午夜av观看不卡| 最后的刺客免费高清国语| 国产又色又爽无遮挡免| 亚洲国产av影院在线观看| 亚洲人与动物交配视频| 色网站视频免费| 男女国产视频网站| 亚洲图色成人| 丝袜人妻中文字幕| 婷婷色综合www| 少妇被粗大的猛进出69影院 | 精品一区在线观看国产| 18禁动态无遮挡网站| 欧美日韩一区二区视频在线观看视频在线| 国产高清国产精品国产三级| 九九爱精品视频在线观看| 永久免费av网站大全| 亚洲,一卡二卡三卡| 两个人免费观看高清视频| 国产一区二区在线观看日韩| 好男人视频免费观看在线| 成年女人在线观看亚洲视频| 视频区图区小说| 久久人人爽人人片av| 久久久久久久久久久免费av| 国产精品欧美亚洲77777| 久久久久久久大尺度免费视频| av网站免费在线观看视频| 毛片一级片免费看久久久久| 午夜视频国产福利| 日韩伦理黄色片| 国产老妇伦熟女老妇高清| 九草在线视频观看| 欧美另类一区| 新久久久久国产一级毛片| 最新的欧美精品一区二区| 亚洲婷婷狠狠爱综合网| 午夜视频国产福利| 精品一区二区三区四区五区乱码 | 久久97久久精品| 在线观看免费视频网站a站| 免费播放大片免费观看视频在线观看| 国产精品一区www在线观看| 色94色欧美一区二区| 亚洲情色 制服丝袜| 国产精品久久久av美女十八| 免费观看a级毛片全部| 在线 av 中文字幕| 亚洲精品美女久久av网站| √禁漫天堂资源中文www| 十八禁高潮呻吟视频| 国产精品无大码| 国产精品99久久99久久久不卡 | 日韩一区二区三区影片| www.av在线官网国产| 精品午夜福利在线看| 久久久久久人妻| 两性夫妻黄色片 | 日本91视频免费播放| 亚洲天堂av无毛| 女性被躁到高潮视频| 久久久久久久久久人人人人人人| 久久久久人妻精品一区果冻| 看非洲黑人一级黄片| 日韩一区二区视频免费看| 欧美精品一区二区大全| 2018国产大陆天天弄谢| 丝袜脚勾引网站| 少妇的逼水好多| 黄片无遮挡物在线观看| 亚洲欧美日韩卡通动漫| av线在线观看网站| av.在线天堂| 欧美精品国产亚洲| 亚洲五月色婷婷综合| 日韩成人伦理影院| 精品福利永久在线观看| 我的女老师完整版在线观看| 最近手机中文字幕大全| 女人被躁到高潮嗷嗷叫费观| 欧美激情极品国产一区二区三区 | 如何舔出高潮| 最新的欧美精品一区二区| 成年人午夜在线观看视频| 性高湖久久久久久久久免费观看| 欧美性感艳星| 在线观看免费视频网站a站| 中文欧美无线码| 亚洲欧美日韩另类电影网站| 人成视频在线观看免费观看| 欧美日韩国产mv在线观看视频| 欧美激情极品国产一区二区三区 | 尾随美女入室| 精品久久国产蜜桃| 国产色爽女视频免费观看| 久久久久久伊人网av| 男人舔女人的私密视频| 美女主播在线视频| 丰满饥渴人妻一区二区三| 日韩人妻精品一区2区三区| 大香蕉97超碰在线| 久久久久久久久久人人人人人人| 国产免费视频播放在线视频| 亚洲久久久国产精品| 黑人欧美特级aaaaaa片| 人妻 亚洲 视频| 日韩av免费高清视频| 男人添女人高潮全过程视频| 777米奇影视久久| 国产亚洲精品第一综合不卡 | 如何舔出高潮| 国产一区二区三区av在线| 欧美成人午夜精品| 久久久久国产网址| 久久人人爽av亚洲精品天堂| 亚洲精品一二三| 黑人猛操日本美女一级片| 成年动漫av网址| 自拍欧美九色日韩亚洲蝌蚪91| 久久国产精品男人的天堂亚洲 | 久久久国产欧美日韩av| 最近最新中文字幕免费大全7| 久久精品国产亚洲av天美| 精品国产乱码久久久久久小说| 久久人人爽人人片av| 亚洲精品aⅴ在线观看| 久久久久网色| 国产成人精品久久久久久| 最后的刺客免费高清国语| 美女脱内裤让男人舔精品视频| 亚洲婷婷狠狠爱综合网| 色网站视频免费| 精品卡一卡二卡四卡免费| 大陆偷拍与自拍| 卡戴珊不雅视频在线播放| 亚洲在久久综合| 日本爱情动作片www.在线观看| 99国产精品免费福利视频| 久久人人爽av亚洲精品天堂| kizo精华| 欧美性感艳星| 色5月婷婷丁香| 久久99热这里只频精品6学生| 在线精品无人区一区二区三| 激情五月婷婷亚洲| 国产伦理片在线播放av一区| 777米奇影视久久| 国产精品麻豆人妻色哟哟久久| 成人国产av品久久久| 女人久久www免费人成看片| 午夜福利视频精品| 亚洲成色77777| 久久精品国产综合久久久 | 免费大片18禁| 久久久欧美国产精品| 日韩不卡一区二区三区视频在线| 99热6这里只有精品| 久久国产精品大桥未久av| 大片电影免费在线观看免费| 黄网站色视频无遮挡免费观看| 亚洲精品美女久久久久99蜜臀 | 满18在线观看网站| 一本大道久久a久久精品| 国产免费福利视频在线观看| 高清欧美精品videossex| 寂寞人妻少妇视频99o| 国产片特级美女逼逼视频| 日韩av不卡免费在线播放| 麻豆乱淫一区二区| 国产 一区精品| 九九爱精品视频在线观看| 最黄视频免费看| 亚洲伊人色综图| 七月丁香在线播放| 黑人巨大精品欧美一区二区蜜桃 | 纯流量卡能插随身wifi吗| 亚洲熟女精品中文字幕| 九九爱精品视频在线观看| 色吧在线观看| 丁香六月天网| 五月玫瑰六月丁香| 国产精品久久久av美女十八| 有码 亚洲区| 日韩成人伦理影院| 最近中文字幕高清免费大全6| 久久久国产精品麻豆| 性色avwww在线观看| 免费女性裸体啪啪无遮挡网站| 一本久久精品| 日本黄大片高清| 97精品久久久久久久久久精品| 亚洲av欧美aⅴ国产| 人妻系列 视频| 26uuu在线亚洲综合色| 一个人免费看片子| 色哟哟·www| 免费在线观看完整版高清| 日本爱情动作片www.在线观看| 黄色一级大片看看| 午夜91福利影院| 国产免费福利视频在线观看| 亚洲精品日本国产第一区| 亚洲av国产av综合av卡| 妹子高潮喷水视频| 极品人妻少妇av视频| 在线观看一区二区三区激情| 激情五月婷婷亚洲| 久久久久精品人妻al黑| 精品久久久精品久久久| 欧美日韩综合久久久久久| 免费日韩欧美在线观看| 国产极品天堂在线| 国产精品一二三区在线看| 美女主播在线视频| 午夜福利视频在线观看免费| √禁漫天堂资源中文www| 尾随美女入室| 精品少妇内射三级| 亚洲一码二码三码区别大吗|