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

    全球變暖不同階段熱帶輻合帶的移動(dòng)及其與大氣能量輸送的關(guān)系?

    2020-02-19 03:04:58鄭小童
    關(guān)鍵詞:北半球赤道短波

    楊 靜, 鄭小童

    (1.中國海洋大學(xué)物理海洋教育部重點(diǎn)實(shí)驗(yàn)室,海洋與大氣學(xué)院 山東 青島 266100;2. 青島海洋科學(xué)與技術(shù)試點(diǎn)國家實(shí)驗(yàn)室, 山東 青島 266237)

    作為熱帶海區(qū)經(jīng)向上集中的強(qiáng)降雨帶,同時(shí)也是哈德萊(Hadley)環(huán)流圈的上升中心,熱帶輻合帶(Intertropical Convergence Zone, ITCZ)是熱帶氣候系統(tǒng)的重要組成部分[1]。ITCZ位置和強(qiáng)度的微小變化會(huì)引起熱帶地區(qū)顯著的降水和環(huán)流改變,進(jìn)而造成整個(gè)熱帶甚至全球的氣候異常,理解并預(yù)測ITCZ的變化是氣候研究中的重要問題。

    在年平均氣候態(tài)上,ITCZ并非與接受太陽短波輻射通量最多的赤道重合,而是常年處于赤道以北8°N附近,并隨著季節(jié)南北移動(dòng)[2]。傳統(tǒng)觀點(diǎn)認(rèn)為是由于南北半球海陸分布和大陸岸界走向的影響,南北半球海陸分布差異使得赤道以北海表面溫度(Sea Surface Temperature, SST)略高于赤道以南,在風(fēng)-蒸發(fā)-SST反饋[3]、上升流-風(fēng)-SST反饋[4]和層云-SST反饋[5]等熱帶局地海氣相互作用下維持ITCZ一直位于赤道以北。

    考慮到ITCZ對應(yīng)Hadley環(huán)流的上升支,由于濕靜力能隨著高度增加而增大,因此高空支主要與能量輸送有關(guān),相應(yīng)的低空支主要與水汽輸送有關(guān),如此將大氣能量與ITCZ聯(lián)系起來[6]。近年來,一些學(xué)者從大氣能量的角度出發(fā),指出ITCZ的位置不僅受熱帶地區(qū)局地作用的影響,熱帶外地區(qū)對大氣的能量強(qiáng)迫也會(huì)對其產(chǎn)生影響[7-8]。Kang等從能量的角度提出了一個(gè)新的理論框架來解釋ITCZ對于熱帶外地區(qū)熱力強(qiáng)迫的響應(yīng),通過冷卻(加熱)40°N(S)以北(南)地區(qū)來改變南北半球大氣能量差異,發(fā)現(xiàn)北半球副熱帶地區(qū)向極地的渦旋通量增加來補(bǔ)償北溫帶地區(qū)的失熱,使得北半球熱帶地區(qū)熱量減少,熱帶地區(qū)內(nèi)部通過Hadley環(huán)流重新調(diào)整,大氣高層產(chǎn)生向北的跨赤道能量輸送,低層產(chǎn)生向南的水汽輸送,引起ITCZ向南移動(dòng)[9]。

    基于以上能量框架理論,人們對ITCZ的位置提出新的解釋。例如Frierson等指出了海洋經(jīng)向翻轉(zhuǎn)流對于ITCZ位于赤道以北的貢獻(xiàn):由于海洋向北輸送熱量,北半球熱帶外地區(qū)大氣得到更多的熱量,因而北半球副熱帶地區(qū)向高緯度輸送的能量較南半球副熱帶地區(qū)少,熱帶高層大氣產(chǎn)生向南的跨赤道能量輸送,導(dǎo)致ITCZ位于赤道以北[10]。同時(shí),氣候模式中雙ITCZ的系統(tǒng)性模擬偏差也可以由南半球中高緯度海區(qū)低云模擬偏差導(dǎo)致大氣得到過多的能量來解釋[11]。

    工業(yè)革命以來溫室氣體排放導(dǎo)致的全球變暖對全球氣候有重要影響。在全球變暖下,前人的研究也發(fā)現(xiàn)了跨赤道的能量輸送,進(jìn)而影響ITCZ的位置。例如全球變暖下大西洋經(jīng)向翻轉(zhuǎn)流(Atlantic Meridional Overturning Circulation, AMOC)將會(huì)減弱,引起異常的大氣向北跨赤道能量輸送,造成ITCZ南移[12]。Frierson和Hwang利用CMIP3中9個(gè)混合層海洋-大氣耦合模式探究了全球變暖下ITCZ的響應(yīng),雖然各個(gè)模式中ITCZ的位移與大氣跨赤道能量輸送差異較大,但兩者呈顯著的負(fù)相關(guān),通過歸因分析,溫帶地區(qū)的云和冰的響應(yīng)可以解釋大部分的ITCZ位移,而熱帶地區(qū)由水汽、云引起的正反饋會(huì)放大溫帶地區(qū)的熱力作用[13]。McFarlane和Frierson利用CMIP5模式在RCP8.5排放情景下研究了全球變暖下ITCZ的變化,發(fā)現(xiàn)大氣跨赤道能量輸送變化與ITCZ的移動(dòng)之間呈很好的負(fù)相關(guān),不同模式中,各個(gè)能量項(xiàng)引起的跨赤道能量輸送的方向基本一致,但數(shù)值上有較大的差異,多模式平均下跨赤道能量輸送幾乎為零,氣溶膠濃度、海洋熱含量、大氣層頂凈的短波輻射的變化使得ITCZ北移,而大氣層頂凈的長波輻射、海洋熱輸送的變化會(huì)造成ITCZ南移,ITCZ整體移動(dòng)不大[14]。

    由于海洋熱容量遠(yuǎn)遠(yuǎn)大于地球其他圈層,海洋在氣候系統(tǒng)中扮演著“熱容器”的角色,延長了氣候系統(tǒng)(特別是地球表面溫度)對外強(qiáng)迫響應(yīng)的時(shí)間尺度[15]。因此,全球變暖特別是全球平均溫度(Global-Mean Surface Temperature, GMST)對外輻射強(qiáng)迫變化的響應(yīng)存在快(幾年)、慢(幾百到上千年)兩種尺度,分別對應(yīng)于混合層的快速響應(yīng)、深層海洋的緩慢響應(yīng)[16]。在溫室氣體快速增長的過程中,海洋混合層因吸收大量熱量而快速增暖,并且向深層海洋傳輸,將熱量儲(chǔ)存在深海中,因此延緩了GMST的增加。而在溫室氣體濃度達(dá)到穩(wěn)定后時(shí),儲(chǔ)存在深海里的熱量會(huì)對混合層海洋產(chǎn)生反饋,使GMST在輻射強(qiáng)迫保持不變時(shí)仍呈一個(gè)緩慢增長的趨勢。相對應(yīng)的SST增暖形態(tài)對溫室氣體先增長后保持不變的響應(yīng)同樣也分為快響應(yīng)和慢響應(yīng)兩個(gè)階段[17]。在快響應(yīng)下,CO2快速增加,深層海洋的增暖延緩了SST的增暖,由于南半球海洋面積遠(yuǎn)大于北半球,因而北半球增暖大于南半球;而在慢響應(yīng)下,CO2趨于穩(wěn)定不變,深層海洋反過來加熱上層海洋,南半球增暖大于北半球。由于SST和能量有著密切的關(guān)系,大氣能量傳輸方向是從增暖的半球傳到冷卻的半球[18-19],快慢響應(yīng)中不同的SST增暖型必然會(huì)導(dǎo)致ITCZ不同的變化。前人的研究都是針對溫室氣體增加造成的跨赤道能量輸送變化和ITCZ的移動(dòng)。對于溫室氣體達(dá)到穩(wěn)定后,雖然有工作研究了局地海溫和降水變化的關(guān)系[17]以及大氣環(huán)流的響應(yīng)過程[20],但這一階段ITCZ的位置如何變化,以及與南北半球能量傳輸?shù)年P(guān)系,還未有人給出明確的答案。

    綜上所述,本文在前人研究的基礎(chǔ)上,對CMIP5多模式數(shù)據(jù)集進(jìn)行分析,給出了全球變暖不同階段ITCZ的位置變化及其與能量的關(guān)系,并對快慢響應(yīng)階段進(jìn)行比較,指出了海洋對于氣候變化的重要調(diào)控作用,這對理解全球變暖背景下全球的大氣環(huán)流和水循環(huán)變化有重要意義,為通過大氣能量收支預(yù)估未來熱帶降水變化提供了理論基礎(chǔ)。

    1 資料和方法

    本文使用了第五次國際間耦合模式比較計(jì)劃(CMIP5)中12個(gè)模式的歷史模擬(Historical)和典型濃度路徑2100年達(dá)到4.5 W/m2(Representative Concentration Pathways, RCP4.5)模擬試驗(yàn)的月平均資料結(jié)果(所用模式見表格1)。在RCP4.5試驗(yàn)中,輻射強(qiáng)迫穩(wěn)定增加,于2070年左右基本達(dá)到4.5 W/m2,然后保持穩(wěn)定不變,并且包含氣溶膠濃度的減少。如圖1所示,大多數(shù)模式中GMST在2070年之前快速增加(多模式平均變化速率為0.023 ℃/a),2070年之后增長速率顯著變慢(0.003 ℃/a)。雖然GMST增長速率的拐點(diǎn)為2070年,本文仍沿用Long等文章中的定義:快響應(yīng)為RCP4.5試驗(yàn)中21世紀(jì)后50年(2051—2100年)平均與歷史試驗(yàn)后50年(1956—2005年)平均之差,慢響應(yīng)為RCP4.5在23世紀(jì)后50年(2251—2300年)平均與22世紀(jì)前50年(2101—2150年)平均之差[17]。

    表1 選取的CMIP5模式

    (采用11年滑動(dòng)平均。11 years moving average.)

    大氣向北跨赤道能量輸送使用McFarlane和Frierson文章中的定義:

    AHTcross-eq=

    即為半球間(南半球減北半球)能量不對稱性。其中:QA為大氣凈的能量收支,包括地球表面的感熱、潛熱通量、以及地球表面和大氣層頂?shù)拈L短波輻射通量;a為地球半徑;λ為經(jīng)度;φ為緯度[14]。為了方便計(jì)算,本文將能量的方向統(tǒng)一為進(jìn)入大氣為正,即地球表面通量(Surface flux)向上為正,大氣層頂通量(Top of the Atmosphere flux, TOA flux)向下為正。

    此外,ITCZ位移用熱帶南北降水不對稱性[21]、熱帶降水中心、熱帶降水重心[22]的變化來表征:

    熱帶降水重心φcent=

    為探究全球變暖不同階段下ITCZ位移與能量的關(guān)系,本文還使用近似部分輻射擾動(dòng)法(Approximate Partial Radiative Perturbation, APRP)對大氣層頂短波輻射進(jìn)行分解,得到云輻射、非云輻射、表面反照率輻射響應(yīng)對短波輻射的貢獻(xiàn),其中云輻射響應(yīng)可分為云反射、散射、吸收輻射響應(yīng);非云輻射響應(yīng)反映大氣的作用,可分為大氣散射(主要與硫酸鹽氣溶膠的散射作用有關(guān))、大氣吸收輻射響應(yīng)(主要與黑炭、水汽、臭氧等物質(zhì)的吸收作用有關(guān))[23]。此外還運(yùn)用了多模式集成分析、相關(guān)分析等統(tǒng)計(jì)方法。由于各個(gè)模式分辨率不一致,本文在進(jìn)行所有計(jì)算之前將數(shù)據(jù)插值到1°×1°網(wǎng)格點(diǎn)上。

    2 全球變暖快響應(yīng)階段ITCZ的移動(dòng)及其與能量的關(guān)系

    2.1 全球變暖快響應(yīng)階段跨赤道AHT的變化

    在RCP4.5試驗(yàn)中全球變暖快響應(yīng)下,大多數(shù)模式顯示大氣存在向南的跨赤道能量輸送(見圖2(a)),多模式平均值為-0.03 PW(以南半球多為正)。本文將大氣跨赤道AHT進(jìn)行分解診斷,考察surface flux和TOA flux對跨赤道AHT的貢獻(xiàn),發(fā)現(xiàn)在大多數(shù)模式中,南半球surface flux大于北半球,多模式平均值為0.07 PW,而TOA flux則相反,顯示為北半球大于南半球,多模式平均值為-0.10 PW,說明快響應(yīng)下大氣跨赤道AHT由TOA flux決定。本文再將TOA flux分為TOA長波與TOA短波,發(fā)現(xiàn)所有模式都模擬北半球得到較多的TOA 短波輻射,多模式平均值為-0.23 PW,是跨赤道向北AHT的主要貢獻(xiàn)因素;相反多模式平均TOA長波輻射的南北差異為0.13 PW,這與快響應(yīng)下北半球增暖更多,大氣層頂向外釋放更多長波輻射有關(guān)。為進(jìn)一步探究TOA短波輻射對跨赤道AHT的調(diào)控機(jī)制,本文利用APRP方法將TOA短波輻射分解為表面反照率、云和非云輻射響應(yīng)(見圖2(a)右側(cè)部分)。通過分解發(fā)現(xiàn)大多數(shù)模式都模擬北半球云、表面反照率、非云輻射響應(yīng)大于南半球,對多模式平均跨赤道AHT的貢獻(xiàn)值分別為-0.08、-0.08和-0.05 PW,可見多模式平均TOA短波輻射通量的南北差異主要與云、表面反照率輻射響應(yīng)有關(guān)。表面反照率輻射響應(yīng)與全球變暖下北極放大效應(yīng)有關(guān),其造成北半球冰雪融化大于南半球,北半球反照率減少使得從TOA進(jìn)來的短波輻射增多,而云輻射響應(yīng)可進(jìn)一步歸因至云反射、散射、吸收輻射響應(yīng)的變化,三者多模式平均值分別為-0.03、-0.05和-0.007 PW,因此云輻射響應(yīng)的半球間差異主要與云短波反射、散射輻射響應(yīng)有關(guān)。以上結(jié)果與前人的研究結(jié)果相似[14]。

    圖2 RCP4.5試驗(yàn)中由輻射引起的大氣向北跨赤道能量輸送在(a) 快響應(yīng)階段和(b)慢響應(yīng)階段的變化(單位:PW)

    2.2 全球變暖快響應(yīng)階段ITCZ的移動(dòng)及其與大氣能量收支的模式間關(guān)系

    對于溫室氣體增加的快響應(yīng)階段ITCZ的移動(dòng),前人研究主要聚焦于模式間平均結(jié)果,對于模式間差異的原因以及對應(yīng)的大氣能量空間分布研究較少,本文針對這一問題進(jìn)行了分析。在RCP4.5試驗(yàn)中全球變暖快響應(yīng)下,12個(gè)模式中有8個(gè)模式顯示ITCZ北移。ITCZ的移動(dòng)與跨赤道AHT變化之間呈顯著的負(fù)相關(guān)(r=-0.79,見圖3(a)),即模式中大氣存在向北(南)的異常能量輸送時(shí),ITCZ南(北)移,能量變化與ITCZ移動(dòng)的對應(yīng)關(guān)系與前人的研究一致[14]。本文對跨赤道AHT進(jìn)行分解,發(fā)現(xiàn)ITCZ的移動(dòng)與TOA flux變化之間呈顯著的負(fù)相關(guān)(r=-0.81,見圖3(c)),與surface flux變化之間呈顯著的正相關(guān)(r=0.70,見圖3(b)),因此快響應(yīng)下ITCZ的移動(dòng)主要由TOA flux驅(qū)動(dòng)。同時(shí)我們看到所有模式都模擬南半球表面溫度(Surface Temperature, TS)增暖小于北半球(見圖3(d)),并且模式間南北半球TS差異的變化與ITCZ位置移動(dòng)和跨赤道AHT的變化顯著相關(guān),相關(guān)系數(shù)分別達(dá)到-0.76和0.85。本文進(jìn)一步將TOA flux分解為TOA 短波輻射通量和長波輻射通量兩項(xiàng),發(fā)現(xiàn)短波輻射通量的南北差異變化與ITCZ的移動(dòng)呈顯著的負(fù)相關(guān)(r=-0.78,見圖4(a),而長波相反,其南北半球差異變化與ITCZ移動(dòng)呈正相關(guān)(r=0.34)。為了進(jìn)一步考察短波輻射通量變化影響ITCZ移動(dòng)的具體物理原因,本文使用APRP方法對TOA短波輻射通量變化進(jìn)行診斷分析,發(fā)現(xiàn)ITCZ移動(dòng)的模式間差異主要由云、表面反照率輻射響應(yīng)貢獻(xiàn),相關(guān)系數(shù)分別為-0.55和-0.68(見圖4(b),(c)),而云輻射響應(yīng)的貢獻(xiàn)主要與云短波反射、散射輻射響應(yīng)有關(guān)(見圖4(d),(e))。因此快響應(yīng)下云反射、云散射、表面反照率輻射響應(yīng)是引起ITCZ模式間差異的主要因子。盡管云、非云短波吸收輻射響應(yīng)也與ITCZ移動(dòng)存在模式間相關(guān),但其半球間差異變化過小,對ITCZ移動(dòng)的貢獻(xiàn)可以忽略不計(jì)(見圖4(f)和(h))。

    ((a) 總能量;(b)表面通量;(c)大氣層頂能量;(d)表面溫度。實(shí)線為線性擬合結(jié)果。 (a) Total energy; (b) Surface flux; (c) TOA flux; (d) Surface temperature. The solid line denotes the linear regression.)

    圖3 快響應(yīng)下ITCZ移動(dòng)與半球間能量以及溫度不對稱性(南半球減北半球)變化散點(diǎn)圖

    Fig.3 Scatter plots of ITCZ shift with interhemispheric energy and temperature asymmetry (southern minus northern hemisphere) change in the fast response

    影響ITCZ移動(dòng)的大氣能量變化具有明顯的空間分布特征。在快響應(yīng)下,調(diào)控ITCZ南北移動(dòng)的TOA flux存在顯著的南北半球差異(見圖5(a)),北半球大多數(shù)地區(qū)為正相關(guān)(陸地上最顯著),南半球大多數(shù)地區(qū)為負(fù)相關(guān),即北半球從TOA進(jìn)入大氣的能量較南半球多時(shí),ITCZ北移。ITCZ與水平平均TOA flux的相關(guān)分別在南北半球低緯度和中高緯度較高(見圖5(b)),其中低緯度的相關(guān)反映了ITCZ位置的南北移動(dòng):當(dāng)ITCZ北移時(shí),對流也隨之北移,南(北)半球熱帶地區(qū)向外長波輻射(Outgoing longwave radiation, OLR)增強(qiáng)(減弱),TOA flux減小(增大)。而中高緯度TOA flux的變化是ITCZ移動(dòng)的驅(qū)動(dòng)因素。本文對TOA flux進(jìn)行分解,主要考察與ITCZ南北移動(dòng)模式間差異相關(guān)的云反射、云散射、表面反照率輻射響應(yīng)的空間分布特征。發(fā)現(xiàn)南大洋地區(qū)顯著的負(fù)值區(qū)與云反射輻射響應(yīng)有關(guān)(見圖6(a)),此地區(qū)海溫較低,上空以低云為主,對短波輻射有較強(qiáng)的阻擋作用,其云反射輻射響應(yīng)的變化對ITCZ的移動(dòng)有重要貢獻(xiàn)。而北半球中高緯度顯著的正值區(qū)與表面反照率和云散射輻射響應(yīng)都有關(guān)(見圖6(b),(c)),其中與云散射輻射響應(yīng)有關(guān)的信號主要出現(xiàn)在大陸上,作者猜測這與RCP4.5試驗(yàn)中在21世紀(jì)南北半球不對稱的氣溶膠排放減少有關(guān)[24],而表面反照率輻射響應(yīng)與北半球高緯度地區(qū)的冰蓋變化有一定聯(lián)系[12]。

    ((a) TOA短波;(b)云輻射響應(yīng);(c)表面反照率輻射響應(yīng);(d)云反射輻射響應(yīng);(e)云散射輻射響應(yīng);(f)云吸收輻射響應(yīng);(g)非云散射輻射響應(yīng);(h)非云吸收輻射響應(yīng)。實(shí)線為線性擬合結(jié)果。 (a) TOA shortwave radiation; (b) Cloud radiative response; (c) Surface albedo radiative response; (d) Cloud reflection radiative response; (e) Cloud scattering radiative response; (f) Cloud absorption radiative response; (g) Noncloud scattering radiative response; (h) Noncloud absorption radiative response. The solid line denotes the linear regression.)

    圖4 快響應(yīng)下ITCZ移動(dòng)與半球間TOA短波輻射不對稱性變化以及APRP分解的散點(diǎn)圖

    Fig.4 Scatter plots of ITCZ shift with interhemispheric TOA shortwave asymmetry change and APRP decomposition in the fast response

    ( (a) 空間分布;(b)水平平均。黑點(diǎn)表示過95%置信度檢驗(yàn)的區(qū)域,虛線為95%置信度的臨界值。 (a) Spatial pattern; (b) Zonal mean. Black dots in (a) denote areas exceeding the 95% confidence level; Dashed lines in (b) mean the 95% confidence level.)

    圖5 快響應(yīng)下ITCZ移動(dòng)與TOA flux變化相關(guān)圖

    Fig.5 Correlation of ITCZ shift with TOA flux change in the fast response

    ( (a) 云反射輻射響應(yīng);(b)云散射輻射響應(yīng);(c)表面反照率輻射響應(yīng)。黑點(diǎn)表示過95%置信度檢驗(yàn)的區(qū)域。(a) Cloud reflection radiative response; (b) Cloud scattering radiative response; (c) Surface albedo radiative response. Black dots denote areas exceeding the 95% confidence level.)

    圖6 快響應(yīng)下ITCZ移動(dòng)與能量變化相關(guān)圖

    Fig.6 Correlation of ITCZ shift with energy change in the fast response

    3 全球變暖慢響應(yīng)階段ITCZ的移動(dòng)及其與能量的關(guān)系

    3.1 全球變暖慢響應(yīng)階段跨赤道AHT的變化

    針對慢響應(yīng)階段的氣候特征,前人就局地海溫和降水變化的關(guān)系[17]以及大氣環(huán)流的響應(yīng)過程[20]開展了研究工作,但這一階段ITCZ的移動(dòng)及其與大氣能量傳輸?shù)年P(guān)系仍不明確,這也是本文要解決的另一個(gè)問題。通過分析,作者發(fā)現(xiàn),在RCP4.5試驗(yàn)中全球變暖慢響應(yīng)下(見圖2(b)),由于GMST增暖幅度較小,其造成的半球間能量差異也小于快響應(yīng)。與快響應(yīng)不同的是,大多數(shù)模式都模擬大氣向北輸送能量,多模式平均值為0.02 PW。將慢響應(yīng)下跨赤道AHT進(jìn)行分解,發(fā)現(xiàn)在大多數(shù)模式中南半球surface flux大于北半球,多模式平均值為0.03 PW,而TOA flux則相反,多模式平均值為-0.01 PW,說明慢響應(yīng)階段南北半球能量差異主要是surface flux變化的貢獻(xiàn)。進(jìn)一步將surface flux分解,發(fā)現(xiàn)多模式平均表面通量的南北不對稱性主要由表面短波輻射通量、潛熱通量引起(多模式平均值均為0.01 PW),其中潛熱通量存在很大的模式間不確定性。

    3.2 全球變暖慢響應(yīng)階段ITCZ的移動(dòng)及其與大氣能量收支的模式間關(guān)系

    在RCP4.5試驗(yàn)中全球變暖慢響應(yīng)下,大多數(shù)模式模擬ITCZ南移,由于溫室氣體濃度已經(jīng)趨于穩(wěn)定,ITCZ移動(dòng)的幅度較快響應(yīng)顯著減小,但依然受跨赤道AHT的調(diào)控(r=-0.74,見圖7(a))。與快響應(yīng)階段不同的是,這種調(diào)控關(guān)系主要由surface flux貢獻(xiàn),ITCZ與surface flux的變化呈顯著的模式間負(fù)相關(guān)(r=-0.56,見圖7(b)),而與TOA flux的變化無模式間相關(guān)性(見圖7(c))。大多數(shù)模式都模擬南半球增暖大于北半球增暖(見圖7(d)),跨赤道AHT的分布仍然與TS的分布有關(guān)(r=0.75)。本文將surface flux分解為潛熱通量、感熱通量、短波輻射和長波輻射通量,分別考察與ITCZ變化的模式間關(guān)系(見圖8),發(fā)現(xiàn)ITCZ與surface flux的模式間相關(guān)性主要由潛熱通量貢獻(xiàn)(見圖8(a)),因而潛熱通量是慢響應(yīng)下影響ITCZ移動(dòng)的主要因子。

    慢響應(yīng)下調(diào)控ITCZ南北移動(dòng)的surface flux在水平平均上存在一定的南北半球差異(見圖9(a)),南半球大多數(shù)地區(qū)為負(fù)相關(guān),北半球大多數(shù)地區(qū)為正相關(guān),這樣的特征在表面潛熱通量上也有所體現(xiàn),并且中緯度45°S~55°S地區(qū)貢獻(xiàn)最大(見圖9(b)),作者猜測這與南大洋西風(fēng)急流的強(qiáng)度變化有關(guān)。由于風(fēng)速變化帶來的潛熱通量變化及其對應(yīng)的海洋吸熱改變,是慢響應(yīng)下決定不同模式ITCZ南北移動(dòng)的主要因素。后續(xù)作者計(jì)劃對這一部分工作進(jìn)行專門研究。

    4 結(jié)論與討論

    4.1 結(jié)論

    本文在前人研究的基礎(chǔ)上利用CMIP5多模式數(shù)據(jù)分析了全球變暖不同階段跨赤道AHT、ITCZ位置的變化及其物理機(jī)制,主要得到以下結(jié)論:

    (1)快響應(yīng)下,大多數(shù)模式模擬大氣向南輸送能量,ITCZ向北移動(dòng),與跨赤道AHT的變化呈顯著的模式間負(fù)相關(guān)。這種負(fù)相關(guān)關(guān)系主要由TOA flux南北差別的模式間差異貢獻(xiàn),而surface flux與ITCZ的移動(dòng)在模式間呈正相關(guān)。進(jìn)一步對TOA flux各項(xiàng)進(jìn)行診斷分析,發(fā)現(xiàn)這一階段ITCZ移動(dòng)的模式間差異主要受南大洋地區(qū)云短波反射輻射響應(yīng)以及北半球中高緯度地區(qū)云短波散射、表面短波反照率輻射響應(yīng)的影響。

    (2)慢響應(yīng)下,大多數(shù)模式模擬大氣向北輸送能量,對應(yīng)ITCZ向南移動(dòng),兩者在模式間仍呈顯著的負(fù)相關(guān),但此時(shí)的負(fù)相關(guān)主要由surface flux南北差別的模式間差異貢獻(xiàn)。進(jìn)一步對surface flux各項(xiàng)進(jìn)行診斷分析,發(fā)現(xiàn)這一階段ITCZ模式間差異主要與潛熱通量的南北半球差異有關(guān),其中45°S~55°S的潛熱通量變化貢獻(xiàn)最大。

    ( (a) 總能量;(b)表面通量;(c)大氣層頂能量;(d)表面溫度。實(shí)線為線性擬合結(jié)果。(a) Total energy; (b) Surface flux; (c) TOA flux; (d) Surface temperature. The solid line denotes the linear regression.)

    圖7 慢響應(yīng)下ITCZ移動(dòng)與半球間能量以及溫度不對稱性變化散點(diǎn)圖

    Fig.7 Scatter plots of ITCZ shift with interhemispheric energy and temperature asymmetry change in the slow response

    ((a) 表面潛熱通量;(b)表面感熱通量;(c)表面短波輻射通量;(d)表面長波輻射通量。實(shí)線為線性擬合結(jié)果。 (a) Surface latent heat; (b) Surface sensible heat; (c) Surface shortwave radiation flux; (d) Surface longwave radiation flux. The solid line denotes the linear regression.)

    圖8 慢響應(yīng)下ITCZ移動(dòng)與半球間能量不對稱性變化散點(diǎn)圖

    Fig.8 Scatter plots of ITCZ shift with interhemispheric energy asymmetry change in the slow response

    ((a)表面通量;(b)表面潛熱通量。虛線為95%置信度的臨界值。(a) Surface flux; (b) Surface latent heat flux. The dashed lines mean the 95% confidence level.)

    圖9 慢響應(yīng)下ITCZ移動(dòng)與水平平均能量變化相關(guān)圖

    Fig.9 Correlation of ITCZ shift with zonal mean energy change in the slow response

    4.2 討論

    根據(jù)對南北半球大氣能量收支差異和ITCZ南北移動(dòng)關(guān)系的分析,本文發(fā)現(xiàn)在全球變暖的快慢響應(yīng)兩個(gè)階段,ITCZ移動(dòng)都和跨赤道AHT的變化顯著相關(guān)(見圖3(a)和圖7(a)),但兩者變化的原因在兩個(gè)階段中是不同的。在快響應(yīng)階段,氣候系統(tǒng)主要受外強(qiáng)迫和氣候反饋的影響發(fā)生變化,模式間的ITCZ位置和跨赤道AHT受TOA flux調(diào)控(見圖3(c)),同時(shí)surface flux變化的南北半球差異與ITCZ的移動(dòng)呈正相關(guān)(見圖3(b))。當(dāng)TOA 向下的能量在北(南)半球較多(少)時(shí),會(huì)導(dǎo)致向南的跨赤道AHT以及ITCZ的向北移動(dòng),同時(shí)北(南)半球海洋(即surface flux變化)也得到更多(少)的能量,造成海洋的跨赤道能量向南的輸運(yùn)。這體現(xiàn)了海洋和大氣對氣候強(qiáng)迫和反饋南北半球不對稱信號同向的跨赤道能量輸送響應(yīng)。同時(shí),在快響應(yīng)階段ITCZ的南北移動(dòng)與TS增暖的南北不對稱顯著相關(guān)(見圖3(d)),也體現(xiàn)了surface flux對表面加熱的不對稱特征。

    相較而言,在慢響應(yīng)階段,ITCZ的位置以及跨赤道AHT主要受surface flux的南北半球不對稱影響(見圖7(b))。在這個(gè)時(shí)期,氣候強(qiáng)迫信號(溫室氣體和氣溶膠)達(dá)到穩(wěn)定,因此海表增暖在模式中的差異會(huì)通過影響surface flux來影響ITCZ的移動(dòng)。在這一階段,大氣通過surface flux在南(北)半球得到更多(少)的能量,造成ITCZ的南移以及向北的跨赤道AHT,同時(shí)相應(yīng)的海洋在北(南)半球得熱更多(少)。在慢響應(yīng)階段,大氣能量在南半球表面得到更多能量的原因主要是在快響應(yīng)階段,南大洋會(huì)吸收大部分進(jìn)入海洋的熱量[25],在溫室氣體穩(wěn)定后會(huì)反饋到海表,造成了大氣能量收支的南北不對稱。因此在慢響應(yīng)階段,南北半球的表面溫度變化差異也與跨赤道AHT呈顯著的正相關(guān)。

    為了進(jìn)一步理解快慢響應(yīng)過程下ITCZ移動(dòng)的原因,本文對不同階段TOA flux和surface flux的能量收支進(jìn)行了診斷,并研究了與ITCZ移動(dòng)相關(guān)的空間分布特征。在快響應(yīng)階段,與ITCZ移動(dòng)相關(guān)的TOA flux主要體現(xiàn)在北半球中高緯度區(qū)域以及南半球45°S位置(見圖5(a),(b))。其中北半球中高緯度地區(qū)與表面短波反照率和云短波散射輻射響應(yīng)有關(guān),而南半球45°S地區(qū)與云短波反射輻射響應(yīng)有關(guān)(見圖6)。并且北半球的輻射響應(yīng)受氣溶膠的影響,由于歷史試驗(yàn)和RCP4.5情景試驗(yàn)中存在人類活動(dòng)排放氣溶膠的作用,且氣溶膠的減少主要發(fā)生在北半球,引起北半球增暖,北極冰雪融化反照率減小,同時(shí)氣溶膠影響云的生成,使得北半球云對短波的散射作用減弱,進(jìn)入大氣的TOA短波輻射較南半球多,引起ITCZ北移。在慢響應(yīng)階段,對surface flux進(jìn)一步診斷發(fā)現(xiàn),ITCZ的南北移動(dòng)主要與潛熱通量的南北半球差異有關(guān),其中45°S~55°S的潛熱通量變化貢獻(xiàn)最大(見圖9(b))。

    通過本文的分析,初步明確了快慢響應(yīng)下跨赤道AHT、ITCZ位置的變化,并且發(fā)現(xiàn)快響應(yīng)階段主要由TOA flux代表的氣候強(qiáng)迫和反饋信號決定,而慢反應(yīng)下主要由表征海洋吸熱變化的surface flux驅(qū)動(dòng)。兩個(gè)階段ITCZ移動(dòng)驅(qū)動(dòng)原因的差異體現(xiàn)了不同階段氣候變化的物理本質(zhì),即快響應(yīng)階段的氣候變化信號是對輻射強(qiáng)迫的直接響應(yīng),而慢響應(yīng)階段的氣候變化信號是對海洋吸熱調(diào)整的響應(yīng)。以上結(jié)果為預(yù)估未來熱帶降水和大氣能量平衡的關(guān)系提供了理論依據(jù)。

    猜你喜歡
    北半球赤道短波
    “雪龍2”號過赤道,救下一船人!
    軍事文摘(2024年2期)2024-01-10 01:58:44
    北半球最強(qiáng)“星空攝影師”開工啦
    軍事文摘(2023年24期)2023-12-19 06:50:06
    清涼一夏
    距赤道越遠(yuǎn)越喜愛黃色
    南北半球天象
    軍事文摘(2019年18期)2019-09-25 08:09:22
    最熱的地方不在赤道
    樂海短波
    人民音樂(2016年1期)2016-11-07 10:02:42
    工運(yùn)短波
    工運(yùn)短波
    綠野短波
    国产精品一区二区三区四区久久| 亚洲成人中文字幕在线播放| 亚洲,欧美,日韩| 亚洲人成网站在线播| 一边摸一边抽搐一进一小说| 十八禁网站免费在线| 男女啪啪激烈高潮av片| 亚洲av五月六月丁香网| 午夜精品一区二区三区免费看| 一进一出抽搐动态| 免费观看人在逋| 亚洲av.av天堂| 久久热精品热| 国产精品乱码一区二三区的特点| 亚洲国产精品久久男人天堂| 最新中文字幕久久久久| 日韩大尺度精品在线看网址| 美女高潮的动态| 欧美色欧美亚洲另类二区| 在线观看午夜福利视频| 亚洲av美国av| 亚洲精品日韩av片在线观看| 久久久久久九九精品二区国产| 99视频精品全部免费 在线| 国产精品不卡视频一区二区| 99久久成人亚洲精品观看| 精品无人区乱码1区二区| 欧美成人a在线观看| 精品久久国产蜜桃| 婷婷色综合大香蕉| 免费在线观看影片大全网站| 欧美精品啪啪一区二区三区| 成人美女网站在线观看视频| 国产精品伦人一区二区| 久久久国产成人免费| 亚洲欧美日韩卡通动漫| 天堂√8在线中文| 国产高清激情床上av| av天堂在线播放| 99视频精品全部免费 在线| 欧美日本视频| 中文在线观看免费www的网站| aaaaa片日本免费| av视频在线观看入口| 国产欧美日韩精品亚洲av| 3wmmmm亚洲av在线观看| 国产三级中文精品| 久久久久国产精品人妻aⅴ院| 国产在线精品亚洲第一网站| 亚洲成a人片在线一区二区| 天天躁日日操中文字幕| or卡值多少钱| 久久这里只有精品中国| 此物有八面人人有两片| 亚洲成人免费电影在线观看| 国产成人一区二区在线| 日日夜夜操网爽| 一本精品99久久精品77| 精品午夜福利视频在线观看一区| 亚洲五月天丁香| 中文字幕熟女人妻在线| 国产91精品成人一区二区三区| 亚洲精品在线观看二区| 国产乱人伦免费视频| 在线观看舔阴道视频| 日韩大尺度精品在线看网址| 亚洲美女搞黄在线观看 | eeuss影院久久| 国产精品久久久久久亚洲av鲁大| 一区二区三区激情视频| 一个人看的www免费观看视频| 久久久久免费精品人妻一区二区| 免费观看精品视频网站| 久久精品国产鲁丝片午夜精品 | 午夜精品一区二区三区免费看| 国产精华一区二区三区| 精品一区二区免费观看| 国产三级在线视频| 老熟妇仑乱视频hdxx| 亚洲专区国产一区二区| netflix在线观看网站| 欧美一级a爱片免费观看看| 亚洲美女搞黄在线观看 | 男插女下体视频免费在线播放| 91在线精品国自产拍蜜月| 亚洲最大成人av| 97超视频在线观看视频| eeuss影院久久| 日本与韩国留学比较| 日本免费一区二区三区高清不卡| 看十八女毛片水多多多| 91狼人影院| 看片在线看免费视频| 免费观看精品视频网站| 国产在线男女| 精品一区二区三区av网在线观看| 欧美日韩亚洲国产一区二区在线观看| 精品无人区乱码1区二区| 国产精品一区www在线观看 | 成人高潮视频无遮挡免费网站| 国产精品野战在线观看| 亚洲电影在线观看av| 国内精品一区二区在线观看| 欧美性猛交╳xxx乱大交人| 美女大奶头视频| 国产主播在线观看一区二区| 国产精品免费一区二区三区在线| 欧美又色又爽又黄视频| 亚洲欧美日韩东京热| 国产精品久久电影中文字幕| 97热精品久久久久久| 亚洲av日韩精品久久久久久密| 国产精品av视频在线免费观看| 波多野结衣高清无吗| 国产精品亚洲一级av第二区| 99久久精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 99国产精品一区二区蜜桃av| 精品久久久久久久久av| 少妇被粗大猛烈的视频| 人人妻,人人澡人人爽秒播| 欧美精品啪啪一区二区三区| 国产 一区 欧美 日韩| 久久精品国产亚洲网站| 免费搜索国产男女视频| 欧美性猛交╳xxx乱大交人| 国产高清视频在线观看网站| 日韩精品有码人妻一区| 又粗又爽又猛毛片免费看| 伊人久久精品亚洲午夜| 免费观看精品视频网站| 亚洲男人的天堂狠狠| 久久久久性生活片| 校园人妻丝袜中文字幕| 九九热线精品视视频播放| 一级a爱片免费观看的视频| 永久网站在线| 欧美激情久久久久久爽电影| 欧美最黄视频在线播放免费| 中文字幕免费在线视频6| 少妇人妻精品综合一区二区 | 嫁个100分男人电影在线观看| 亚洲性夜色夜夜综合| 十八禁网站免费在线| 99热网站在线观看| avwww免费| 国产精品无大码| 一进一出抽搐gif免费好疼| 麻豆成人av在线观看| 色综合色国产| 亚洲欧美日韩东京热| 亚洲国产日韩欧美精品在线观看| 日日夜夜操网爽| 久久人妻av系列| 午夜福利成人在线免费观看| 国产精品一区www在线观看 | 亚洲成av人片在线播放无| 亚洲在线自拍视频| 久久精品国产鲁丝片午夜精品 | 国产精品亚洲美女久久久| 精品一区二区三区av网在线观看| 国产人妻一区二区三区在| 很黄的视频免费| 国产精品久久久久久久久免| 蜜桃亚洲精品一区二区三区| 国产av麻豆久久久久久久| 天天一区二区日本电影三级| 两个人的视频大全免费| 久久精品国产清高在天天线| 国产亚洲91精品色在线| 久久久精品欧美日韩精品| 亚洲aⅴ乱码一区二区在线播放| 一本久久中文字幕| 熟女人妻精品中文字幕| 91精品国产九色| 国产精品永久免费网站| 联通29元200g的流量卡| 亚洲午夜理论影院| 日韩精品中文字幕看吧| 九九在线视频观看精品| 麻豆久久精品国产亚洲av| 男人舔奶头视频| 亚洲乱码一区二区免费版| 国产一区二区在线观看日韩| 88av欧美| 观看美女的网站| 丝袜美腿在线中文| 久久久色成人| 女人被狂操c到高潮| 又爽又黄a免费视频| 午夜免费男女啪啪视频观看 | 校园春色视频在线观看| 精品人妻一区二区三区麻豆 | 最近最新免费中文字幕在线| 日韩亚洲欧美综合| 能在线免费观看的黄片| 亚洲人成网站高清观看| 亚洲性久久影院| 熟妇人妻久久中文字幕3abv| 99久久九九国产精品国产免费| 亚洲中文字幕日韩| 18禁黄网站禁片免费观看直播| 色吧在线观看| 少妇裸体淫交视频免费看高清| 午夜激情福利司机影院| 久久久精品大字幕| 婷婷精品国产亚洲av| 91av网一区二区| 成人特级av手机在线观看| 成人av在线播放网站| 亚洲人成网站高清观看| 日本成人三级电影网站| 97人妻精品一区二区三区麻豆| 88av欧美| 嫩草影院入口| 特大巨黑吊av在线直播| 日韩在线高清观看一区二区三区 | 国产精品久久久久久久电影| 国产一区二区三区av在线 | 韩国av在线不卡| 尤物成人国产欧美一区二区三区| 少妇高潮的动态图| 亚洲精品一区av在线观看| 琪琪午夜伦伦电影理论片6080| 成人av在线播放网站| 中文字幕人妻熟人妻熟丝袜美| 日日摸夜夜添夜夜添小说| 女同久久另类99精品国产91| 一a级毛片在线观看| 免费看美女性在线毛片视频| 日韩亚洲欧美综合| 欧美极品一区二区三区四区| 伦精品一区二区三区| 午夜福利在线观看免费完整高清在 | 在线播放无遮挡| 色综合色国产| 亚洲av二区三区四区| ponron亚洲| 婷婷六月久久综合丁香| 亚洲狠狠婷婷综合久久图片| 嫁个100分男人电影在线观看| 日本撒尿小便嘘嘘汇集6| 日本三级黄在线观看| 欧美高清成人免费视频www| 亚洲专区国产一区二区| 婷婷精品国产亚洲av在线| 日本黄大片高清| 国产真实伦视频高清在线观看 | 99久久精品国产国产毛片| 国产精品一区www在线观看 | 两个人的视频大全免费| 亚洲精品日韩av片在线观看| 亚洲18禁久久av| 一夜夜www| 别揉我奶头~嗯~啊~动态视频| 一卡2卡三卡四卡精品乱码亚洲| 午夜精品一区二区三区免费看| 亚洲欧美日韩东京热| 三级国产精品欧美在线观看| 一级a爱片免费观看的视频| 亚洲av美国av| 精品免费久久久久久久清纯| 美女cb高潮喷水在线观看| 日本黄色片子视频| 欧美一级a爱片免费观看看| 欧美+日韩+精品| 欧美日韩综合久久久久久 | 国产一区二区在线观看日韩| 成年女人永久免费观看视频| www.www免费av| 两人在一起打扑克的视频| 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲精品综合一区在线观看| 国产熟女欧美一区二区| 一级a爱片免费观看的视频| 女生性感内裤真人,穿戴方法视频| 色精品久久人妻99蜜桃| 欧美3d第一页| 日日摸夜夜添夜夜添小说| 免费av不卡在线播放| 精品人妻1区二区| 欧美黑人欧美精品刺激| 69av精品久久久久久| 欧美色欧美亚洲另类二区| 免费在线观看日本一区| 在线观看66精品国产| 亚洲美女黄片视频| 天堂√8在线中文| 欧美xxxx黑人xx丫x性爽| 一级av片app| 赤兔流量卡办理| 国产一区二区三区视频了| 久久精品夜夜夜夜夜久久蜜豆| 乱人视频在线观看| 国产真实乱freesex| 中文字幕人妻熟人妻熟丝袜美| 欧美激情久久久久久爽电影| 人人妻,人人澡人人爽秒播| 国产精品久久久久久精品电影| 极品教师在线免费播放| 欧美色视频一区免费| 中文字幕人妻熟人妻熟丝袜美| 国产毛片a区久久久久| 一区二区三区四区激情视频 | 成年免费大片在线观看| 国产真实乱freesex| 黄色丝袜av网址大全| 老司机深夜福利视频在线观看| 91午夜精品亚洲一区二区三区 | 免费黄网站久久成人精品| 91在线观看av| 成人毛片a级毛片在线播放| 我的女老师完整版在线观看| 免费看美女性在线毛片视频| 午夜激情福利司机影院| 亚洲三级黄色毛片| 国产精品亚洲一级av第二区| 精品国产三级普通话版| 观看美女的网站| 桃色一区二区三区在线观看| av专区在线播放| 亚洲国产高清在线一区二区三| 午夜a级毛片| 制服丝袜大香蕉在线| 国产真实乱freesex| 欧美一级a爱片免费观看看| 亚洲人成网站在线播放欧美日韩| 婷婷色综合大香蕉| 亚洲av一区综合| bbb黄色大片| 成人欧美大片| 久久久久久九九精品二区国产| 亚洲18禁久久av| 给我免费播放毛片高清在线观看| 久久久久久久久久久丰满 | 一区二区三区高清视频在线| 深夜精品福利| 一本一本综合久久| 欧美一级a爱片免费观看看| 久久香蕉精品热| 少妇丰满av| 一夜夜www| 在线免费观看不下载黄p国产 | 久久久久久伊人网av| 国产精品久久视频播放| 欧美高清性xxxxhd video| 白带黄色成豆腐渣| 一a级毛片在线观看| 欧美zozozo另类| 亚洲人成网站在线播| 看十八女毛片水多多多| 成人国产麻豆网| 中国美女看黄片| 直男gayav资源| 日本成人三级电影网站| 日韩欧美一区二区三区在线观看| 久久午夜福利片| av黄色大香蕉| 日本成人三级电影网站| 亚洲综合色惰| 春色校园在线视频观看| 久久久精品大字幕| 国产高潮美女av| 国产一区二区激情短视频| 人妻少妇偷人精品九色| 国产激情偷乱视频一区二区| 此物有八面人人有两片| 亚洲狠狠婷婷综合久久图片| xxxwww97欧美| av在线亚洲专区| 亚洲乱码一区二区免费版| 制服丝袜大香蕉在线| 亚洲av不卡在线观看| 成年版毛片免费区| 国产精品免费一区二区三区在线| 国产亚洲精品久久久com| 99热精品在线国产| 日韩欧美 国产精品| 琪琪午夜伦伦电影理论片6080| 午夜激情欧美在线| 色5月婷婷丁香| 在线观看一区二区三区| 久久午夜福利片| 国产高清三级在线| 欧美一区二区精品小视频在线| 亚洲精品色激情综合| 在线免费十八禁| 免费av观看视频| 91av网一区二区| 一边摸一边抽搐一进一小说| 九色国产91popny在线| 久久精品91蜜桃| 亚洲国产欧美人成| 桃色一区二区三区在线观看| 免费看光身美女| 国产免费av片在线观看野外av| 国产亚洲精品久久久com| 午夜视频国产福利| 国产真实伦视频高清在线观看 | 午夜免费男女啪啪视频观看 | 亚洲专区中文字幕在线| 男女之事视频高清在线观看| 老司机福利观看| 国产精品亚洲一级av第二区| 天堂√8在线中文| 97超级碰碰碰精品色视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 久久精品夜夜夜夜夜久久蜜豆| 桃红色精品国产亚洲av| 成人欧美大片| 无遮挡黄片免费观看| 最近最新免费中文字幕在线| 中文字幕久久专区| 精品人妻1区二区| 亚洲aⅴ乱码一区二区在线播放| 国产综合懂色| 午夜福利高清视频| 日韩欧美免费精品| 色尼玛亚洲综合影院| 欧美日韩黄片免| or卡值多少钱| 国产精品久久电影中文字幕| 午夜福利在线观看免费完整高清在 | 亚洲av第一区精品v没综合| 亚洲精品乱码久久久v下载方式| 亚洲第一区二区三区不卡| 性色avwww在线观看| 国内精品久久久久精免费| 国产av在哪里看| 国产精品久久久久久av不卡| 免费人成在线观看视频色| 日本 av在线| 久久九九热精品免费| 大又大粗又爽又黄少妇毛片口| 日本五十路高清| 嫩草影院入口| 窝窝影院91人妻| 18禁黄网站禁片免费观看直播| 男人的好看免费观看在线视频| 午夜视频国产福利| 国产不卡一卡二| 久久午夜亚洲精品久久| 国产亚洲精品综合一区在线观看| 精品人妻偷拍中文字幕| 最近视频中文字幕2019在线8| 免费黄网站久久成人精品| 国内精品美女久久久久久| 狠狠狠狠99中文字幕| 99热这里只有是精品50| 国产精品福利在线免费观看| 中亚洲国语对白在线视频| 黄色欧美视频在线观看| 免费看a级黄色片| 最后的刺客免费高清国语| 日韩一区二区视频免费看| 一a级毛片在线观看| 99在线视频只有这里精品首页| 欧美激情国产日韩精品一区| 免费电影在线观看免费观看| 精品无人区乱码1区二区| av天堂中文字幕网| 久久久精品欧美日韩精品| 国产精品永久免费网站| 免费看a级黄色片| 春色校园在线视频观看| 91狼人影院| 美女高潮喷水抽搐中文字幕| 久久久国产成人精品二区| 夜夜看夜夜爽夜夜摸| 美女xxoo啪啪120秒动态图| 中文在线观看免费www的网站| 搡老岳熟女国产| 亚洲人与动物交配视频| 免费av不卡在线播放| 精品久久国产蜜桃| 国产精品永久免费网站| 国产精品久久久久久亚洲av鲁大| 最近最新免费中文字幕在线| 波多野结衣巨乳人妻| 久久久成人免费电影| 国产视频内射| 动漫黄色视频在线观看| 日韩欧美免费精品| 亚洲av成人精品一区久久| 精华霜和精华液先用哪个| 国产淫片久久久久久久久| 三级男女做爰猛烈吃奶摸视频| 一级毛片久久久久久久久女| 淫妇啪啪啪对白视频| 国产 一区精品| 99精品久久久久人妻精品| 一个人看的www免费观看视频| 露出奶头的视频| 久久人妻av系列| 国产不卡一卡二| 午夜激情福利司机影院| 国产视频一区二区在线看| 99久久无色码亚洲精品果冻| 看黄色毛片网站| 日韩欧美一区二区三区在线观看| 又爽又黄无遮挡网站| 久久热精品热| 高清毛片免费观看视频网站| 夜夜看夜夜爽夜夜摸| 午夜福利视频1000在线观看| 欧美绝顶高潮抽搐喷水| 日韩欧美免费精品| 日韩大尺度精品在线看网址| 日韩中字成人| 日韩欧美精品免费久久| 亚洲国产日韩欧美精品在线观看| 深夜精品福利| 91麻豆av在线| 一区二区三区高清视频在线| 久久久久久久精品吃奶| 联通29元200g的流量卡| 黄片wwwwww| 能在线免费观看的黄片| av福利片在线观看| av在线天堂中文字幕| 人妻丰满熟妇av一区二区三区| 精品久久久久久久久亚洲 | 国产成人aa在线观看| 人妻丰满熟妇av一区二区三区| 成年女人看的毛片在线观看| 日本a在线网址| 免费一级毛片在线播放高清视频| 亚洲美女视频黄频| 中文字幕精品亚洲无线码一区| 九九在线视频观看精品| 欧美日韩乱码在线| 99精品在免费线老司机午夜| 男女做爰动态图高潮gif福利片| xxxwww97欧美| 麻豆久久精品国产亚洲av| 国产精品久久久久久av不卡| 国产亚洲精品av在线| 国产精品人妻久久久影院| 欧美另类亚洲清纯唯美| 国产女主播在线喷水免费视频网站 | 日本免费a在线| 99久久精品一区二区三区| 91麻豆精品激情在线观看国产| 国产日本99.免费观看| 精品一区二区三区人妻视频| 一个人观看的视频www高清免费观看| 国产激情偷乱视频一区二区| 久久99热这里只有精品18| 久久久国产成人精品二区| 国产精品综合久久久久久久免费| 欧美成人a在线观看| 亚洲,欧美,日韩| 少妇的逼好多水| 老女人水多毛片| 日韩精品有码人妻一区| 99热6这里只有精品| 三级国产精品欧美在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲成av人片在线播放无| 国产国拍精品亚洲av在线观看| 欧美成人免费av一区二区三区| 人人妻,人人澡人人爽秒播| 最近在线观看免费完整版| av天堂中文字幕网| 亚洲七黄色美女视频| 久久久久免费精品人妻一区二区| 亚洲一级一片aⅴ在线观看| 成人亚洲精品av一区二区| 男人舔奶头视频| 国产亚洲av嫩草精品影院| 亚洲国产精品成人综合色| 看片在线看免费视频| 美女黄网站色视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品一区av在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲真实伦在线观看| 日韩 亚洲 欧美在线| 日本黄大片高清| 日本一二三区视频观看| 在线a可以看的网站| 国产精品久久久久久av不卡| 少妇熟女aⅴ在线视频| 舔av片在线| 窝窝影院91人妻| 国产欧美日韩精品亚洲av| av在线老鸭窝| 国产精品一及| 尤物成人国产欧美一区二区三区| avwww免费| 九九爱精品视频在线观看| 一a级毛片在线观看| 91麻豆av在线| av国产免费在线观看| 国产单亲对白刺激| 午夜福利18| 女人被狂操c到高潮| 国产精品福利在线免费观看| 久久久久性生活片| 人人妻,人人澡人人爽秒播| 51国产日韩欧美| 亚洲精品亚洲一区二区| 美女 人体艺术 gogo| 亚洲 国产 在线| 欧美xxxx性猛交bbbb| 97碰自拍视频| 国产麻豆成人av免费视频| 在线观看av片永久免费下载| 亚洲成人久久爱视频| 久久精品人妻少妇| 国产麻豆成人av免费视频| 亚洲 国产 在线| 在线免费观看不下载黄p国产 | 免费人成视频x8x8入口观看| 色综合色国产|