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

    FGOALS-g3 模擬的南亞夏季風(fēng):氣候態(tài)和年際變率

    2022-10-09 08:22:54何林強(qiáng)周天軍李立娟林鵬飛陳曉龍鄒立維
    大氣科學(xué) 2022年5期
    關(guān)鍵詞:海溫南亞印度洋

    何林強(qiáng) 周天軍 李立娟 林鵬飛 陳曉龍 鄒立維

    1 中國科學(xué)院大氣物理研究所大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國家重點(diǎn)實(shí)驗(yàn)室(LASG), 北京 100029

    2 中國科學(xué)院大學(xué)地球與行星科學(xué)學(xué)院, 北京 100049

    1 引言

    南亞夏季風(fēng)主要受太陽輻射季節(jié)變化和大尺度海陸熱力差異的影響,是全球季風(fēng)系統(tǒng)的重要組成部分,每年印度半島的夏季風(fēng)降水占到全年的80%(Wang et al., 2017),對當(dāng)?shù)氐乃Y源、社會(huì)經(jīng)濟(jì)和生態(tài)環(huán)境具有顯著影響?;跉夂蛳到y(tǒng)模式,認(rèn)識(shí)和理解南亞夏季風(fēng)的歷史變化規(guī)律、預(yù)測和預(yù)估南亞夏季風(fēng)在未來全球增暖背景下的變化,是科學(xué)界高度關(guān)注的前沿問題(Huang et al., 2020a,2020b)。

    由于復(fù)雜的?!憽?dú)庀嗷プ饔?,氣候系統(tǒng)模式對于南亞夏季風(fēng)的模擬仍然是一個(gè)巨大的挑戰(zhàn)。在氣候態(tài)上,CM?P3 和CM?P5 多模式集合平均的南亞夏季降水均在赤道西印度洋存在濕偏差,季風(fēng)槽內(nèi)存在干偏差(Sperber et al., 2013)。已有的研究主要從海氣耦合、大氣物理過程以及地形和分辨率等方面討論南亞夏季風(fēng)氣候態(tài)模擬偏差的原因。在海氣耦合方面,赤道印度洋的東風(fēng)應(yīng)力偏差通過皮耶克尼斯正反饋,形成緯向海溫梯度偏差,造成西印度洋偏暖,降水增多,并進(jìn)一步減弱了南亞夏季風(fēng)環(huán)流(Annamalai et al., 2017)。耦合模式對觀測中印度洋海溫—蒸發(fā)—降水相關(guān)性模擬的偏差,也會(huì)導(dǎo)致大氣對局地海溫強(qiáng)迫的響應(yīng)偏強(qiáng),進(jìn)而影響南亞季風(fēng)區(qū)降水的模擬(Bollasina and Nigam,2009)。此外,耦合模式模擬的阿拉伯海冷偏差會(huì)減少夏季西南氣流向印度半島的水汽輸送,造成季風(fēng) 槽 降 水 偏 少(?zumo et al., 2008; Levine et al.,2013)。對比有無海洋同化的數(shù)值試驗(yàn)表明,耦合模式對印度夏季風(fēng)降水的模擬偏差與西太平洋海溫模擬偏差有關(guān)(王晨琦等, 2019)。在大氣物理過程方面,對流參數(shù)化過程導(dǎo)致的大氣瞬變偏差會(huì)直接造成赤道西印度洋降水的模擬偏差,進(jìn)而影響南亞夏季風(fēng)的模擬(Martin et al., 2010)。大氣模式中積云對流方案的夾卷系數(shù)對于云的形成和降水的分 布 有 顯 著 影 響(Neale et al., 2008; Zhao et al.,2018),通過增加NCAR-CAM4 模式中對流方案的夾卷率,敏感性試驗(yàn)結(jié)果表明赤道西印度洋的濕偏差和季風(fēng)槽的干偏差有明顯改善(Hanf and Annamalai, 2020)。在地形作用方面,觀測結(jié)果和數(shù)值試驗(yàn)都指出濕靜力能的水平梯度在印度半島北部和西側(cè)的沙漠之間達(dá)到最大,表明由于山脈隔絕了西側(cè)的干冷空氣,南亞季風(fēng)區(qū)的濕靜力能和降水才 得 以 維 持(Boos and Kuang, 2010; Boos and Hurley, 2013)。因此,當(dāng)前模式對于高原陡峭地形的平滑處理,會(huì)導(dǎo)致季風(fēng)槽降水的模擬偏少(Boos and Hurley, 2013)。此外,提高模式的水平分辨率,可以有效改善西高止山和中南半島迎風(fēng)坡地形降水的模擬偏差,但對于印度大陸的降水模擬沒有明顯提升(Johnson et al., 2016)。

    在年際尺度上,El Ni?o 與印度夏季降水(All-?ndia Rainfall,簡稱A?R)存在著顯著的負(fù)相關(guān)關(guān)系,這也是衡量氣候模式模擬南亞夏季風(fēng)性能的重要指標(biāo)(Sperber et al., 2013; Pandey et al., 2020)。在CM?P3 的眾多模式中,僅有4 個(gè)能再現(xiàn)這種熱帶太平洋海溫異常激發(fā)的大氣遙相關(guān)(Annamalai et al., 2007)。相比之下,CM?P5 模式對于El Ni?o與A?R 負(fù)相關(guān)關(guān)系的模擬有了明顯改善(Sperber et al., 2013)。基于CM?P5 多模式歷史氣候模擬的分析表明,El Ni?o—Southern Oscillation(ENSO)與季風(fēng)關(guān)系的模擬取決于模式對于氣候態(tài)季風(fēng)降水的空間分布,ENSO 相關(guān)的赤道太平洋海溫和非絕熱加熱的模擬(Sperber et al., 2013)。同時(shí),受El Ni?o 強(qiáng)迫的印度洋的海溫異常會(huì)使得熱帶太平洋海溫與印度夏季降水的遙相關(guān)關(guān)系減弱(Ashok et al., 2004),因此,ENSO 相關(guān)的印度洋和西太平洋海溫異常的模擬也會(huì)直接影響到ENSO 與季風(fēng)降水的關(guān)系(Lau and Nath, 2009; Achuthavarier et al., 2012)。

    全球季風(fēng)模式比較計(jì)劃(GMM?P)作為第六次耦合模式比較計(jì)劃(CM?P6)的子計(jì)劃之一,由中國、英國和美國學(xué)者聯(lián)合發(fā)起組織(Zhou et al.,2016)。GMM?P 通過比較不同分辨率、有無海氣耦合過程的試驗(yàn),有望明晰季風(fēng)降水模擬偏差的來源,為提升模式綜合模擬能力提供參考(周天軍等, 2019a)。中國科學(xué)院大氣物理研究所大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國家重點(diǎn)實(shí)驗(yàn)室(LASG/?AP)發(fā)展的最新版本的耦合系統(tǒng)模式FGOALSg3(Li et al., 2020b),參與了GMM?P 國際計(jì)劃。季風(fēng)的模擬,一直是FGOALS 模式研發(fā)過程中重點(diǎn)關(guān)注的問題。以FGOALS-g2 為例,該模式能夠合理地模擬出南亞夏季風(fēng)的氣候態(tài)和年際變率特征,但由于氣候態(tài)下印度洋海溫和對流層溫度等模擬偏差,導(dǎo)致印度半島的水汽輸送和動(dòng)力抬升減弱,季風(fēng)槽出現(xiàn)干偏差;El Ni?o 期間印太暖池的海溫模擬偏差造成沃克環(huán)流下沉支位置偏移,ENSO 與A?R 負(fù)相關(guān)關(guān)系偏弱(黃昕等, 2019)。從FGOALS-g2到FGOALS-g3,該模式從大氣分量到海洋分量都做了諸多改進(jìn),但是,關(guān)于該模式對南亞夏季風(fēng)的模擬能力尚未進(jìn)行系統(tǒng)評估。本文的目的是:(1)基于觀測事實(shí),評估FGOALS-g3 對于南亞夏季風(fēng)氣候態(tài)和年際變率的模擬能力;(2)通過耦合與非耦合試驗(yàn)結(jié)果的比較,理解海氣耦合過程對于南亞季風(fēng)模擬偏差的影響。

    2 資料和方法

    2.1 模式介紹

    FGOALS-g3 是由LASG/?AP 發(fā)展的耦合系統(tǒng)模式,包含了大氣、海洋、陸面和海冰四個(gè)模塊。相比于FGOALS-g2,F(xiàn)GOALS-g3 的大氣分量GAM?L3 采用了二維混合的并行剖分,在水平分辨率、水汽平流方案,物理過程和外強(qiáng)迫上進(jìn)行了更新和改進(jìn) (Li et al., 2020a)。GAM?L3 的水平分辨率由約2.8°(128×60)提高到約2°(180×80),并調(diào)整了兩步保形平流方案以提高水汽守恒(TSPAS;Yu, 1994)。在物理過程方面,GAM?L3 考慮了對流動(dòng)量傳輸過程(Wu et al., 2007)、采用了基于估計(jì)反演強(qiáng)度的積云對流參數(shù)化方案(E?S;Guo and Zhou, 2014)、人為氣溶膠效應(yīng)的簡單參數(shù)化方案(Stevens et al., 2017; Shi et al., 2019)和改進(jìn)的邊界層方案(Sun et al., 2016)。外強(qiáng)迫采用了CM?P6 推 薦 的 方 案(Nie et al., 2019),相 比FGOALS-g2 增加了火山活動(dòng)的強(qiáng)迫,同時(shí)加入了1965~2014 年人為地下水開采作為陸面模式的額外強(qiáng)迫(Li et al., 2020b)。

    FGOALS-g3 的 海 洋 分 量L?COM3 相 對 于L?COM2 的改進(jìn)主要表現(xiàn)在以下幾個(gè)方面:一是動(dòng)力框架由原來經(jīng)緯度格點(diǎn)框架改為適用于任意正交曲面坐標(biāo)并采用三極格點(diǎn),有效解決了北極點(diǎn)引起的一系列計(jì)算穩(wěn)定性問題;二是升級耦合器接口,提高了耦合效率;三是改進(jìn)了潮汐混合以及渦旋混合的物理過程(Lin et al., 2016; 俞永強(qiáng)等, 2018; 于子棚等, 2019; Lin et al., 2020)。FGOALS-g3 的陸面分量CAS-LSM(Xie et al., 2020),海冰分量Los Alamos sea ice model Version 4.0 相比FGOALSg2 均有所改進(jìn)。FGOALS-g3 采用NCAR 開發(fā)的CPL7 進(jìn)行耦合。有關(guān)FGOALS-g3 相對于FGOALSg2 模式改進(jìn)的詳細(xì)介紹和整體評估參見(Li et al.,2020b)。

    本文使用的模式數(shù)據(jù)為FGOALS-g3 的歷史試驗(yàn)和單獨(dú)大氣試驗(yàn)結(jié)果。其中,歷史試驗(yàn)使用觀測的、隨時(shí)間變化的外強(qiáng)迫驅(qū)動(dòng)耦合模式,模擬1850~2014 年的歷史氣候;單獨(dú)大氣試驗(yàn)利用1979~2014 年以來觀測的海溫和海冰作為邊界條件驅(qū)動(dòng)大氣模式GAM?L3,這是CM?P6 的標(biāo)準(zhǔn)設(shè)定(周天軍等, 2019b)。為比較FGOALS-g3 歷史試驗(yàn)與FGOALS-g2 的差異,本文還用到FGOALSg2 向CM?P5 提交的歷史試驗(yàn)數(shù)據(jù)。上述每組試驗(yàn)均選取第一個(gè)集合成員,時(shí)間為1979~2005 年。為理解內(nèi)部變率對本文結(jié)果的影響,本文在討論部分,給出了5 個(gè)樣本集合的模擬效果與單個(gè)集合成員結(jié)果的異同。為方便討論,以下分別用FGOALS-g2,F(xiàn)GOALS-g3 和GAM?L3 來表示FGOALS-g3 與g2的歷史試驗(yàn)和FGOALS-g3 的單獨(dú)大氣試驗(yàn)結(jié)果。

    2.2 觀測和再分析資料

    本文用到的觀測資料包括:(1)GPCP 2.3( Global Precipitation Climatology Project dataset version 2.3)逐月降水(Adler et al., 2003);(2)英國哈德萊中心Had?SST(Hadley Centre Global Sea ?ce and Sea Surface Temperature version 1.1)逐月海表溫度(Rayner et al., 2003);(3)日本氣象廳JRA55(Japanese 55-year reanalysis projects)逐月大氣再分析數(shù)據(jù)(Kobayashi et al., 2015)。為方便比較,所有數(shù)據(jù)均選取1979~2005 年的6~9 月并統(tǒng)一插值到2.5°×2.5°格點(diǎn)上。JRA55 再分析數(shù)據(jù)垂直方向上選取1000、925、850、700、600、500、400、300、250、200、150、100 hPa 各層。

    3 結(jié)果分析

    以下針對南亞夏季風(fēng)的氣候態(tài)及年際變率特征,首先通過比較FGOALS-g3 與觀測、FGOALS-g3與FGOALS-g2、FGOALS-g3 與GAM?L3 的試驗(yàn)結(jié)果,關(guān)注FGOALS-g3 與FGOALS-g2 和海氣耦合過程對于南亞夏季風(fēng)模擬的差異,隨后,再進(jìn)一步探討模擬偏差改進(jìn)的原因。

    3.1 FGOALS-g3 對氣候態(tài)下南亞夏季降水的模擬

    3.1.1 南亞夏季降水氣候態(tài)的模擬

    氣候態(tài)上,觀測結(jié)果(圖1a)表明南亞夏季風(fēng)期間,印度洋越赤道氣流受地轉(zhuǎn)效應(yīng)影響及東非大陸的阻擋,形成西南季風(fēng),并向印度半島和中南半島輸送水汽,降水的中心位于孟加拉灣的季風(fēng)槽內(nèi)。受地形影響,西高止山、若開山脈和喜馬拉雅山脈的迎風(fēng)坡降水較多(圖1a)。位于赤道東印度洋海上輻合帶的降水伴隨著南亞季風(fēng)的撤退而增多(圖1a、圖2a)。

    圖1 南亞夏季降水(填色,單位:mm d-1)與850 hPa 風(fēng)場(矢量,單位:m s-1)6~9 月的氣候態(tài)分布:(a)GPCP/JRA55 觀測;(b)FGOALS-g3 模 擬;(c)GAM?L3 模 擬;(d)FGOALS-g3 模 擬 結(jié) 果 減 去GPCP/JRA55 觀 測 結(jié) 果;(e)FGOALS-g2 模 擬 結(jié) 果 減 去GPCP/JRA55 觀測結(jié)果;(f)FGOALS-g3 減去GAM?L3 模擬結(jié)果。圖中白色區(qū)域表示海拔高度2000 米以上,圖(a)中三組三角形標(biāo)識(shí)自西向東分別表示西高止山、喜馬拉雅山脈和若開山脈Fig. 1 Climatology of JJAS (June–July–August–September) South Asian summer monsoon precipitation (shaded, units: mm d-1) and 850 hPa wind(vectors, units: m s-1): (a) GPCP/JRA55, (b) FGOALS-g3, (c) GAM?L3, (d) FGOALS-g3 minus GPCP/JRA55, (e) FGOALS-g2 minus GPCP/JRA55,and (f) FGOALS-g3 minus GAM?L3. White areas denote altitudes above 2000 m. Three groups of triangular marks in (a) represent the Western Ghats,the Himalayas and the Rakhine Mountains from west to east, respectively

    相比于觀測,F(xiàn)GOALS-g3 模擬的季風(fēng)環(huán)流偏弱,整個(gè)南亞季風(fēng)區(qū)存在明顯的北風(fēng)偏差,減弱了向印度半島北部的水汽輸送(圖1b 和d)。降水場上,西印度洋的濕偏差和季風(fēng)槽的干偏差是CM?P3 和CM?P5 耦合模式普遍存在的系統(tǒng)性偏差(Annamalai et al., 2017),這在FGOALS-g3 的模擬中依然存在(圖1d)。 相比于FGOALS-g2,F(xiàn)GOALS-g3 模擬的赤道印度洋東風(fēng)偏差明顯改善,赤道西印度洋和阿拉伯海的濕偏差減弱,但赤道中印度洋和孟加拉灣的濕偏差加?。▓D1d 和e)。GAM?L3 模擬的南亞夏季風(fēng)環(huán)流和降水的特征與FGOALS-g3 結(jié)果類似(圖1b 和c),表明耦合模式的模擬偏差在很大程度上源自大氣模式自身。但相比于GAM?L3,F(xiàn)GOALS-g3 模擬的北印度洋西風(fēng)偏差減弱,對應(yīng)的印度半島中部和中南半島西部降水濕偏差減弱(圖1f),這反應(yīng)了海氣耦合過程的影響。

    在年循環(huán)特征上,南亞夏季降水集中在6~9 月,雨帶中心北推至22°N 附近,赤道印度洋輻合帶的降水在夏季風(fēng)期間逐漸增強(qiáng)(圖1a、圖2a)。相比于觀測,F(xiàn)GOALS-g3 模擬的南亞夏季風(fēng)主雨帶位置偏南(8°N 附近)且強(qiáng)度偏大,同時(shí)海上輻合帶降水特征不明顯(圖2b),表現(xiàn)為由赤道印度洋延伸至印度半島的“干—濕—干”的三核型偏差(圖2d)。這與FGOALS-g2 的模擬偏差特征類似(黃昕等, 2019),但FGOALS-g3 中的三核型偏差強(qiáng)度更大(圖2d、e),表明年循環(huán)特征的改進(jìn)難度之大。此外,GAM?L3 模擬的年循環(huán)特征與FGOALS-g3 類似,但主雨帶北推位置相比于耦合模式更接近觀測(圖2c、f),表明這種偏差來自大氣模式自身,但是海氣耦合過程進(jìn)一步放大了這一偏差。

    圖2 南亞夏季降水(65°E~95°E 平均)的年循環(huán)(單位:mm d-1):(a)GPCP 觀測;(b)FGOALS-g3 模擬;(c)GAM?L3 模擬;(d)FGOALS-g3 模擬結(jié)果減GPCP 觀測結(jié)果;(e)FGOALS-g2 模擬結(jié)果減GPCP 觀測結(jié)果;(f)FGOALS-g3 減GAM?L3 模擬結(jié)果Fig. 2 Annual cycle climatology of South Asian summer monsoon precipitation averaged from 65°E to 95°E (units: mm d-1): (a) GPCP; (b) FGOALSg3; (c) GAM?L3; (d) FGOALS-g3 minus GPCP; (e) FGOALS-g2 minus GPCP; (f) FGOALS-g3 minus GAM?L3

    3.1.2 南亞夏季降水氣候態(tài)模擬偏差的原因

    南亞夏季風(fēng)的形成和變化與大尺度的經(jīng)向熱力差異有關(guān)(Webster and Yang, 1992; Jin and Wang,2017;圖3)。由于對流潛熱釋放,熱帶對流層上層溫度對海溫異常的響應(yīng)增強(qiáng),使得對流層上層經(jīng)向熱力梯度對于南亞季風(fēng)環(huán)流的驅(qū)動(dòng)作用要大于對流層低層(Dai et al., 2013)。圖3a 給出了JRA55資料計(jì)算的氣候態(tài)下300 hPa 的南亞夏季風(fēng)環(huán)流和溫度。對流層上層海洋相比于大陸更冷,其中冷中心位于西印度洋,暖中心位于印度半島北部(與整層大氣溫度水平梯度的分布類似,圖片未給出)。根據(jù)熱成風(fēng)原理,北印度洋的對流層上層盛行東風(fēng)(圖3a),低層盛行偏西風(fēng)(圖1a)。相比于觀測,F(xiàn)GOALS-g3模擬的300 hPa 溫度在南亞季風(fēng)區(qū)整體偏冷,同時(shí)冷偏差中心位于印度半島北部至中亞一帶(圖3b)。這樣的冷偏差結(jié)構(gòu)使得氣候態(tài)的經(jīng)向溫度梯度減弱,南亞季風(fēng)環(huán)流減弱(圖3b),進(jìn)一步抑制了印度半島北部及季風(fēng)槽附近的對流活動(dòng),出現(xiàn)下沉運(yùn)動(dòng)偏差(圖3f),導(dǎo)致整層水汽輻散偏差(圖4d),降水模擬偏少(圖1d)。此外,印度半島南部的抬升運(yùn)動(dòng)偏差及對應(yīng)的降水濕偏差廣泛存在于FGOALS-g2,F(xiàn)GOALSg3 以及GAM?L3 的模擬中(圖1b、e 和f),這可能與模式的分辨率難以準(zhǔn)確刻畫位于印度德干高原西部的西高止山的陡峭地形有關(guān)。

    圖3 氣候態(tài)6~9 月南亞地區(qū)300 hPa 溫度(左列,填色,單位:K)、風(fēng)場(左列,矢量,單位:m s-1)和500 hPa 垂直速度(右列,單位:10-2 Pa s-1)的空間分布:(a, e)JRA55 觀測;(b, f)FGOALS-g3 模擬結(jié)果減JRA55 觀測結(jié)果;(c, g)FGOALS-g2 模擬結(jié)果減JRA55 觀測結(jié)果;(d, h)GAM?L3 模擬結(jié)果減JRA55 觀測結(jié)果Fig. 3 Climatology of JJAS 300 hPa air temperature (left panel, shading, units: K), wind (left panel, vectors, units: m s-1), and 500 hPa vertical velocity (right panel, units: 10-2 Pa s-1): (a, e) JRA55; (b, f) FGOALS-g3 minus JRA55; (c, g) FGOALS-g2 minus JRA55; (d, h) GAM?L3 minus JRA55

    圖4 氣候態(tài)6~9 月印度洋海表溫度(左列,單位:K)和整層水汽通量積分(右列,矢量,單位:kg m-1 s-1)及其散度(右列,填色,單位:10-5 kg m-2 s-1)。(a)FGOALS-g3 模擬減Had?SST 觀測結(jié)果;(b)FGOALS-g2 模擬減Had?SST 觀測結(jié)果;(c)FGOALS-g3 減FGOALS-g2 模擬結(jié)果;(d)FGOALS-g3 減JRA55 觀測結(jié)果;(e)FGOALS-g2 模擬減JRA55 觀測結(jié)果;(f)FGOALS-g3 減GAM?L3模擬結(jié)果Fig. 4 Climatology of JJAS sea surface temperature (left panel, units: K), vertically integrated moisture flux (right panel, vectors, units: kg m-1 s-1)and its divergence (right panel, shading, units: 10-5 kg m-2 s-1): (a) FGOALS-g3 minus Had?SST; (b) FGOALS-g2 minus Had?SST; (c) GAM?L3 minus Had?SST; (d) FGOALS-g3 minus JRA55; (e) FGOALS-g2 minus JRA55; (f) GAM?L3 minus JRA55

    相比于FGOALS-g3,F(xiàn)GOALS-g2 模擬的對流層溫度同樣存在類似的冷偏差結(jié)構(gòu),并且系統(tǒng)性冷偏差的強(qiáng)度較FGOALS-g3 偏冷2 K 左右(圖3b、c),造成局地哈德萊環(huán)流減弱進(jìn)而導(dǎo)致降水的模擬偏差(圖3g;黃昕等, 2019)。GAM?L3 模擬的氣候態(tài)對流層300 hPa 溫度的水平分布上,南亞季風(fēng)區(qū)的經(jīng)向梯度同樣減弱,但相比于FGOALS-g3耦合試驗(yàn)?zāi)M的偏差強(qiáng)度較弱(圖3b、d),對應(yīng)的季風(fēng)槽的下沉運(yùn)動(dòng)偏差(圖3f、h)和降水干偏差也有所減弱(圖1a、c)。這在一定程度上說明,海氣耦合過程會(huì)放大單獨(dú)大氣模式中對流層溫度的冷偏差。

    值得注意的是,造成南亞夏季風(fēng)環(huán)流模擬偏弱的對流層溫度冷偏差不僅存在于FGOALS-g3 模式中,John and Soden(2007)指出CM?P3 耦合模式模擬的對流層溫度系統(tǒng)性地偏冷1~2 K,冷偏差隨著高度增加,并且副熱帶地區(qū)大于熱帶和高緯地區(qū)。類似的冷偏差結(jié)構(gòu)在大多數(shù)CM?P5 耦合模式中依然存在,冷偏差的大值中心位于副熱帶300 hPa附近(Tian et al., 2013)。但是,耦合模式模擬的這種系統(tǒng)性冷偏差的來源尚不清楚,而單獨(dú)大氣模式模擬的對流層溫度也存在相似的偏差結(jié)構(gòu)。因此,這種全局性的冷偏差可能與大氣模式的濕物理過程參數(shù)化方案及動(dòng)力過程的反饋有關(guān)。在南亞季風(fēng)區(qū),對流模擬偏差造成的降水干偏差導(dǎo)致對流層中上層非絕熱加熱偏弱,一方面冷偏差可以通過熱帶波動(dòng)傳播至整個(gè)熱帶副熱帶地區(qū),形成全局性的冷偏差;另一方面,冷偏差通過調(diào)整大氣環(huán)流進(jìn)而影響降水,所構(gòu)成的正反饋過程進(jìn)一步放大了降水的模擬偏差。關(guān)于對流層冷偏差影響季風(fēng)降水的機(jī)制,未來有待通過數(shù)值試驗(yàn)加以明晰。觀測證據(jù)顯示,對流層的中上層變冷能夠令南亞季風(fēng)和東亞季風(fēng)區(qū)降水減少(Yu et al., 2004; Yu and Zhou, 2007; Roxy et al.,2015),并且這種影響通過數(shù)值試驗(yàn)在東亞地區(qū)得到部分驗(yàn)證(Xin et al., 2008)。

    在對流層低層,不同于FGOALS-g2 試驗(yàn)中印度洋海溫“西暖東冷”的偶極型偏差(圖4b;Lin et al., 2013; 黃昕等, 2019),F(xiàn)GOALS-g3 模擬的海溫冷偏差整體減弱(圖4a;Li et al., 2020b),且東印度洋海溫暖偏差強(qiáng)度更大(圖4c),使得赤道印度洋海溫偏差的緯向梯度減小,緯向環(huán)流偏差減弱(圖1f),基本糾正了FGOALS-g2 中赤道印度洋的東風(fēng)偏差(圖1d)。而相比于觀測海溫,耦合模式FGOALS-g3 模擬的北印度洋海溫呈現(xiàn)“西暖東冷”的偶極型偏差,西北太平洋呈現(xiàn)冷偏差(圖4a),使得大氣模式模擬的北印度洋西風(fēng)偏差減弱(圖1f),不利于水汽向南亞季風(fēng)區(qū)輸送(圖4f),在一定程度上改善了大氣模式模擬的印度半島中部和中南半島中部的濕偏差(圖1f)?;贑M?P5 多模式集合平均的結(jié)果也表明,氣候態(tài)季風(fēng)降水模擬的改進(jìn)很大程度上來源于耦合模式模擬的海溫偏差(Yang et al., 2019)。

    3.2 FGOALS-g3 對南亞夏季降水年際變率的模擬

    3.2.1 南亞夏季降水與ENSO 關(guān)系的模擬

    在年際尺度上,ENSO 是印度半島夏季降水重要的可預(yù)報(bào)性來源(Mishra et al., 2012; Cherchi and Navarra, 2013)。觀測中,在低層環(huán)流場上(圖5a),赤道西印度洋的東北風(fēng)異常抑制了西南氣流向印度半島的水汽輸送,同時(shí)海洋性大陸的海溫冷異常通過Gill 響應(yīng)(Gill, 1980)激發(fā)低層反氣旋性環(huán)流異常,進(jìn)而驅(qū)動(dòng)邊界層的Ekman 輻散,導(dǎo)致南亞夏季風(fēng)環(huán)流減弱,印度半島降水減少。

    相比于觀測,F(xiàn)GOALS-g3 模擬的El Ni?o 期間,赤道印度洋低層輻散偏強(qiáng),對應(yīng)的赤道中東印度洋出現(xiàn)干偏差;南亞夏季風(fēng)環(huán)流偏強(qiáng),印度半島降水偏多,干異常偏弱(圖5b, d)。這與FGOALS-g2模擬的ENSO 相關(guān)的南亞夏季風(fēng)環(huán)流和降水的偏差與之類似(黃昕等, 2019),但FGOALS-g3 模擬的印度半島南端和阿拉伯海的濕異常偏差加劇(圖5d, e)。

    GAM?L3 模擬的El Ni?o 期間南亞季風(fēng)區(qū)環(huán)流和降水特征與FGOALS-g3 相似,但印度半島的干異常、印度半島南端以及中南半島的濕異常較之耦合模式偏強(qiáng)(圖5c, f),表明耦合過程減弱了上述偏差的強(qiáng)度。同時(shí),相比于GAM?L3,耦合模式模擬的赤道印度洋低層環(huán)流與觀測更加接近(圖5a–c),有效改善了GAM?L3 中西印度洋的干偏差(圖5f)。

    圖5 標(biāo)準(zhǔn)化的6~9 月Ni?o3.4 指數(shù)回歸的同期南亞降水(填色,單位:mm d-1)和850 hPa 風(fēng)場(矢量,單位:m s-1)異常的空間分布:(a)GPCP/JRA55 觀測;(b)FGOALS-g3 模擬;(c)GAM?L3 模擬;(d)FGOALS-g3 模擬減GPCP/JRA55 觀測結(jié)果;(e)g2-GPCP,(f)g3-GAM?L3。白色區(qū)域表示海拔高度2000 米以上,打點(diǎn)區(qū)域表示回歸系數(shù)通過95%的顯著性檢驗(yàn)Fig. 5 JJAS precipitation anomalies (shading, units: mm d-1) and 850 hPa wind anomalies (vectors, units: m s-1) regressed onto standardized Ni?o3.4 index: (a) GPCP/JRA55; (b) FGOALS-g3; (c) GAM?L3; (d) FGOALS-g3 minus GPCP/JRA55; (e) FGOALS-g2 minus GPCP/JRA55; (f) GAM?L3 minus GPCP/JRA55. White areas denote altitudes above 2000 m and dots areas denote the 95% confidence level

    3.2.2 南亞夏季降水與ENSO 關(guān)系模擬偏差的原因

    El Ni?o 期間,熱帶中東太平洋暖海溫異常激發(fā)的沃克環(huán)流上升支在印太暖池區(qū)下沉,進(jìn)而影響周邊地區(qū)的氣候。為探究模式對ENSO 與南亞夏季風(fēng)關(guān)系模擬偏差的原因,圖6 和圖7 給出了El Ni?o 期間海溫和環(huán)流異常的觀測和模擬結(jié)果。參照黃昕等(2019)的定義計(jì)算得出,觀測中沃克環(huán)流上升支位于155°W 附近,下沉中心位于105°E(圖6a 和圖7a)。同時(shí),El Ni?o 通過“大氣橋”影響印度洋海表的熱通量,使得印度洋海溫異常呈現(xiàn)“西暖東冷”的偶極型分布(Lau and Nath,2009),產(chǎn)生的異常環(huán)流在西印度洋上升東印度洋下沉,有利于南亞夏季風(fēng)的增強(qiáng),在一定程度上減弱了ENSO 對南亞夏季風(fēng)的遙強(qiáng)迫。

    相比于觀測,F(xiàn)GOALS-g3 模擬的赤道中東太平洋El Ni?o 海溫暖異常與觀測特征相符,變率較FGOALS-g2 有 明 顯 提 升(圖6b 和c;Li et al.,2020b; Zhang et al., 2020),但印太暖池區(qū)的海溫偏差依然較大(圖6d)。具體來說,F(xiàn)GOALS-g3試驗(yàn)?zāi)M未能模擬出熱帶西太平洋的海溫冷異常,但赤道東印度洋則出現(xiàn)冷異常,西印度洋出現(xiàn)更強(qiáng)的暖異常,呈偶極子型。這一方面造成沃克環(huán)流下沉支(圖7a 和b;92.5°E 附近)較觀測偏西,通過Gill 響應(yīng)激發(fā)的低層反氣旋異常在南亞地區(qū)也偏西(圖5b 和d),有助于水汽從阿拉伯海輸送至印度半島;同時(shí),較觀測偏強(qiáng)的印度洋偶極子進(jìn)一步抑制了ENSO 對于南亞夏季風(fēng)的遙強(qiáng)迫(Ashok et al., 2001, 2004)。這些過程都使得FGOALS-g3模擬的El Ni?o 期間南亞夏季風(fēng)環(huán)流和對流活動(dòng)較觀測偏強(qiáng),在印度半島出現(xiàn)濕偏差(圖5d)。

    FGOALS-g2 模擬的ENSO 基本特征與觀測接近(Chen et al., 2016),沃克環(huán)流下沉支位于63°E附近(圖7c; 黃昕等, 2019),相比之下,F(xiàn)GOALS-g3由于模擬的El Ni?o 期間赤道東印度洋海溫冷偏差加劇以及赤道西太平洋暖偏差減弱(圖6b, c, f),使得下沉支的模擬(92.5°E)與觀測更接近(圖7b,c),但降水模擬的提升并不明顯(圖5d, e)。這與Li et al.(2020a)和Li et al.(2020b) 結(jié)論一致,即相比于上一版本,F(xiàn)GOALS-g3 模擬的環(huán)流場、云輻射強(qiáng)迫等有所改進(jìn),但降水的改進(jìn)則不顯著。這可能是參數(shù)自動(dòng)優(yōu)化時(shí),目標(biāo)函數(shù)里包含了更多動(dòng)力、輻射場的變量,而降水場權(quán)重較小,因此總體性能雖有改進(jìn),但降水模擬的提升有限。

    圖6 標(biāo)準(zhǔn)化的6~9 月Ni?o3.4 指數(shù)回歸的同期海表溫度異常(填色,單位:K)的空間分布:(a)Had?SST 觀測;(b)FGOALS-g3 模擬;(c)FGOALS-g2 模擬;(d)FGOALS-g3 模擬減Had?SST 觀測結(jié)果;(e)FGOALS-g2 模擬減Had?SST 觀測結(jié)果(f)FGOALS-g3減FGOALS-g2 模擬結(jié)果。打點(diǎn)區(qū)域表示回歸系數(shù)通過95%的顯著性檢驗(yàn)Fig. 6 JJAS sea surface temperature anomalies (units: K) regressed onto standardized Ni?o3.4 index: (a) Had?SST; (b) FGOALS-g3; (c) FGOALSg2; (d) FGOALS-g3 minus Had?SST; (e) FGOALS-g2 minus Had?SST; (f) FGOALS-g3 minus FGOALS-g2. Dots area denote the 95% confidence level

    相比于FGOALS-g3,GAM?L3 由于受觀測海溫驅(qū)動(dòng),模擬的ENSO 相關(guān)的沃克環(huán)流下沉支中心的位置與觀測更為接近(圖7a, d),因此通過Gill 響應(yīng)激發(fā)的反氣旋位置也接近觀測(圖5a, c)。但由于缺少海洋與大氣的耦合過程,未考慮海溫—降水—云短波輻射的負(fù)反饋過程(He et al., 2020),造成模擬的赤道中東太平洋的對流活動(dòng)偏強(qiáng)(圖7d),在海洋大陸地區(qū)的下沉運(yùn)動(dòng)也偏強(qiáng),伴隨著云量的減少,使得短波輻射增加(圖7h)。在FGOALS-g3 的模擬中,短波輻射的增加使得海溫升高,進(jìn)一步促進(jìn)了對流活動(dòng),在一定程度上減弱了先前的下沉異常;而GAM?L3 由于不考慮海洋對大氣的響應(yīng),無法描述上述海氣相互作用的負(fù)反饋過程,因此模擬的動(dòng)力下沉要比觀測偏強(qiáng)(圖7d),導(dǎo)致印度半島中北部的降水響應(yīng)偏強(qiáng)(圖5c),地表接收到的短波輻射也偏強(qiáng)(圖7h)。因此,海氣耦合過程對于南亞夏季風(fēng)環(huán)流和降水年際變率的模擬有顯著影響,一方面考慮海溫—降水—云短波輻射的負(fù)反饋過程能夠減小模擬偏差的強(qiáng)度,但同時(shí)耦合模式中沃克環(huán)流下沉運(yùn)動(dòng)中心偏西也使得陸地季風(fēng)降水的負(fù)異常響應(yīng)位置偏西。因此,印度半島降水的模擬呈現(xiàn)濕偏差。

    圖7 標(biāo)準(zhǔn)化的6~9 月Ni?o3.4 指數(shù)回歸的同期850 hPa 輻散風(fēng)(矢量,單位:m s-1)與速度勢異常(填色,單位:105 m2 s-1)和地表接收到的向下的短波輻射通量異常(單位:W m-2)的空間分布:(a, e)JRA55 觀測;(b, f)FGOALS-g3 模擬;(c, g)FGOALS-g2 模擬;(d, h)GAM?L3 模擬。打點(diǎn)區(qū)域表示回歸系數(shù)通過95%的顯著性檢驗(yàn),所展示的輻散風(fēng)與速度勢均通過95%的顯著性檢驗(yàn)Fig. 7 JJAS 850 hPa divergent wind anomalies (left panel, vectors, units: m s-1), velocity potential anomalies (left panel, shading, units: 105 m2 s-1),and downward shortwave flux anomalies (right panel, units: W m-2) regressed onto standardized Ni?o3.4 index: (a, e) JRA55, (b, f) FGOALS-g3, (c, g)FGOALS-g2, and (d, h) GAM?L3. Dots areas denote the 95% confidence level

    3.3 集合平均結(jié)果的影響

    由于同一個(gè)模式不同的集合成員的差異主要體現(xiàn)在內(nèi)部變率上,因此對于氣候態(tài)南亞夏季風(fēng)的影響可以忽略。而在年際變率上, ENSO 與印度半島夏季降水在1979~2005 年間存在顯著的負(fù)相關(guān),并且這種關(guān)系受到年代際變率的影響;對于耦合模式而言,這種關(guān)系在不同集合成員之間可能存在差異,為此,我們進(jìn)一步分析了FGOALS-g2 和FGOALS-g3 歷史試驗(yàn)5 個(gè)集合成員的模擬結(jié)果。相比于單個(gè)集合成員對于南亞夏季風(fēng)的模擬(圖1和圖5),多樣本集合之后的模擬特征(圖8 和圖9)無論在氣候態(tài)還是年際變率上都與前者一致,沒有明顯差異。

    圖8 同圖1,其中FGOALS-g2 和FGOALS-g3 歷史模擬的結(jié)果為前五個(gè)集合成員的平均值Fig. 8 As in Fig 1, but the results of FGOALS-g2 and FGOALS-g3 are derived from the mean of the first five members

    圖9 同圖5,其中FGOALS-g2 和FGOALS-g3 歷史模擬的結(jié)果為前五個(gè)集合成員的平均值。(b-c)中打點(diǎn)區(qū)域表示五個(gè)集合成員同號(hào)Fig. 9 As in Fig 5, but the results of FGOALS-g2 and FGOALS-g3 are derived from the mean of the first five members. Dots in (b) and (c) denote the five members with the same sign

    4 結(jié)論

    本文通過與觀測和JRA55 再分析資料的對比,系統(tǒng)評估了FGOALS-g3 模式對于南亞夏季風(fēng)氣候態(tài)和年際變率的模擬能力,并重點(diǎn)關(guān)注FGOALSg3 與FGOALS-g2、海氣耦合模式與非耦合模式的差異。

    氣候態(tài)上,F(xiàn)GOALS-g3 模擬的南亞夏季降水由熱帶印度洋至印度半島呈現(xiàn)“干—濕—干”的三核型偏差,即海上輻合帶和陸地季風(fēng)槽降水偏少,赤道西印度洋偏多。進(jìn)一步分析表明,耦合模式模擬的對流層中上層的溫度偏低,導(dǎo)致南亞季風(fēng)區(qū)經(jīng)向溫度梯度減弱,造成季風(fēng)環(huán)流減弱,印度半島及季風(fēng)槽內(nèi)產(chǎn)生下沉運(yùn)動(dòng)偏差以及水汽輻散偏差,抑制對流活動(dòng),導(dǎo)致降水模擬偏弱。相比于FGOALS-g2,F(xiàn)GOALS-g3 對于赤道印度洋東風(fēng)偏差的模擬有明顯改善,這與海溫東西梯度的模擬偏差減小有關(guān)。相比于GAM?L3,海氣耦合過程導(dǎo)致的西北太平洋和北印度洋海溫偏差有效改善了印度半島中部和中南半島的降水濕偏差。

    在年際變率上,F(xiàn)GOALS-g3 模擬的季風(fēng)環(huán)流偏強(qiáng),ENSO—印度夏季降水負(fù)相關(guān)偏弱。分析表明,這與耦合模式模擬的El Ni?o 海溫偏差有關(guān)。一方面,模式未能模擬出熱帶西太平洋海溫冷異常造成印度半島的下沉運(yùn)動(dòng)減弱;另一方面,模式模擬的伴隨El Ni?o 發(fā)生的印度洋偶極子偏強(qiáng),有利于季風(fēng)環(huán)流的增強(qiáng),部分抵消了El Ni?o 對于南亞季風(fēng)區(qū)的遙強(qiáng)迫。相對于FGOALS-g2,F(xiàn)GOALSg3 由于模擬的El Ni?o 期間印太暖池區(qū)偏冷,沃克環(huán)流下沉支位置與觀測更加接近。相比于GAM?L3,耦合模式由于考慮了海溫—降水—云短波輻射的負(fù)反饋過程,El Ni?o 相關(guān)的環(huán)流和降水強(qiáng)度的模擬偏差顯著減小,但印太暖池區(qū)的海溫模擬偏差造成沃克環(huán)流下沉運(yùn)動(dòng)中心偏西,使得陸地季風(fēng)降水的負(fù)響應(yīng)位置偏西,印度半島降水出現(xiàn)濕異常偏差。

    綜上,相比于FGOALS-g2,由于局地海溫的模擬變化,F(xiàn)GOALS-g3 模式對于南亞夏季風(fēng)氣環(huán)流的模擬有明顯改進(jìn),但降水的模擬偏差仍然存在。海氣耦合過程對于南亞夏季風(fēng)氣候態(tài)和年際變率的模擬偏差也存在顯著影響。未來若要改進(jìn)南亞夏季風(fēng)氣候態(tài)的模擬,需要關(guān)注大氣模式的濕物理過程參數(shù)化方案及動(dòng)力過程的反饋;若要改進(jìn)年際變率的模擬,則需要重點(diǎn)關(guān)注耦合模式模擬的El Ni?o期間印太暖池區(qū)的海溫異常。此外,耦合試驗(yàn)和單獨(dú)大氣試驗(yàn)的模擬結(jié)果都顯示,在氣候態(tài)與年際變率上,印度半島南部的降水上存在著相似的濕偏差,意味著降水偏差還與模式分辨率對地形刻畫不足等因素有關(guān),這一點(diǎn)未來有待通過高分辨率模式的數(shù)值試驗(yàn)加以驗(yàn)證。

    猜你喜歡
    海溫南亞印度洋
    南亞磷肥市場需東繼續(xù)提升
    基于深度學(xué)習(xí)的海溫觀測數(shù)據(jù)質(zhì)量控制應(yīng)用研究
    基于Argo、XBT數(shù)據(jù)的蘇拉威西海溫鹽特征分析
    與南亞高壓相聯(lián)的歐亞大陸-印度洋經(jīng)向環(huán)流
    中國熱科院南亞所 研發(fā)成功菠蘿機(jī)械化種植機(jī)
    南印度洋偶極型海溫與中國西南地區(qū)初秋降水的關(guān)系
    去印度洋
    滇池(2018年1期)2018-01-17 21:34:30
    南亞煉油化工市場一覽
    南亞高壓對西藏夏季降水的影響
    西藏科技(2016年5期)2016-09-26 12:16:41
    波浪輸運(yùn)通量剩余量對上層海溫的影響分析
    天天添夜夜摸| 亚洲av美国av| 91av网站免费观看| 免费观看a级毛片全部| 香蕉国产在线看| 亚洲成国产人片在线观看| 国产99久久九九免费精品| 99国产精品一区二区三区| 动漫黄色视频在线观看| 亚洲精品一二三| 久久中文字幕一级| 18在线观看网站| 日本猛色少妇xxxxx猛交久久| 亚洲精品av麻豆狂野| 国产免费一区二区三区四区乱码| av免费在线观看网站| 一个人免费看片子| 成人亚洲精品一区在线观看| 丝袜人妻中文字幕| 欧美大码av| 日日夜夜操网爽| 国产精品久久久久久精品电影小说| 不卡一级毛片| 久热这里只有精品99| 考比视频在线观看| 欧美xxⅹ黑人| 中国美女看黄片| 欧美人与性动交α欧美精品济南到| 亚洲伊人久久精品综合| 各种免费的搞黄视频| 中文字幕最新亚洲高清| 天堂8中文在线网| 国产精品久久久久成人av| 一边摸一边做爽爽视频免费| 波多野结衣av一区二区av| 久久中文看片网| 男女之事视频高清在线观看| 一进一出抽搐动态| 91老司机精品| 又紧又爽又黄一区二区| 久久国产精品男人的天堂亚洲| 午夜福利乱码中文字幕| 90打野战视频偷拍视频| 亚洲五月婷婷丁香| 男女国产视频网站| www日本在线高清视频| 国产亚洲一区二区精品| 91麻豆精品激情在线观看国产 | 王馨瑶露胸无遮挡在线观看| 嫩草影视91久久| 国产真人三级小视频在线观看| 俄罗斯特黄特色一大片| 国产成人一区二区三区免费视频网站| 精品国内亚洲2022精品成人 | 免费高清在线观看日韩| 亚洲国产精品一区二区三区在线| 亚洲国产中文字幕在线视频| 美国免费a级毛片| 国产高清视频在线播放一区 | 国产男女内射视频| 热re99久久国产66热| 亚洲国产成人一精品久久久| 亚洲久久久国产精品| 欧美激情高清一区二区三区| 精品欧美一区二区三区在线| 咕卡用的链子| 国产精品久久久av美女十八| 日韩熟女老妇一区二区性免费视频| av网站免费在线观看视频| 777米奇影视久久| 黄色毛片三级朝国网站| 91字幕亚洲| 亚洲第一青青草原| av网站免费在线观看视频| 国产亚洲欧美在线一区二区| 国产精品成人在线| 热re99久久国产66热| 啦啦啦啦在线视频资源| 国产亚洲欧美在线一区二区| 国产精品亚洲av一区麻豆| 精品国内亚洲2022精品成人 | 香蕉国产在线看| 爱豆传媒免费全集在线观看| 精品福利观看| 男女无遮挡免费网站观看| 黄色视频不卡| 在线av久久热| 亚洲精品国产av成人精品| 国产在线一区二区三区精| 日韩三级视频一区二区三区| 最近最新中文字幕大全免费视频| 亚洲第一欧美日韩一区二区三区 | 日本猛色少妇xxxxx猛交久久| netflix在线观看网站| 建设人人有责人人尽责人人享有的| 1024香蕉在线观看| 青草久久国产| 丁香六月欧美| 国产精品久久久久久精品古装| 日本精品一区二区三区蜜桃| 亚洲人成电影观看| 91精品国产国语对白视频| av免费在线观看网站| 新久久久久国产一级毛片| 国产亚洲av片在线观看秒播厂| 精品国产一区二区三区久久久樱花| 中文欧美无线码| 十八禁人妻一区二区| 老司机福利观看| 午夜激情久久久久久久| 王馨瑶露胸无遮挡在线观看| 在线观看免费视频网站a站| 叶爱在线成人免费视频播放| 美女扒开内裤让男人捅视频| 丝袜脚勾引网站| 国产在线一区二区三区精| 五月开心婷婷网| 青草久久国产| 美国免费a级毛片| 久久精品国产a三级三级三级| 黄色视频,在线免费观看| 制服诱惑二区| av国产精品久久久久影院| 啦啦啦啦在线视频资源| 久久中文看片网| 欧美日韩视频精品一区| 黄色毛片三级朝国网站| 国产淫语在线视频| 最近最新中文字幕大全免费视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美激情久久久久久爽电影 | 亚洲全国av大片| 日韩视频一区二区在线观看| 日本wwww免费看| 久久青草综合色| 黑人巨大精品欧美一区二区蜜桃| videosex国产| 99香蕉大伊视频| 搡老熟女国产l中国老女人| 天天躁夜夜躁狠狠躁躁| 精品人妻熟女毛片av久久网站| 精品熟女少妇八av免费久了| 国产成人a∨麻豆精品| 999久久久精品免费观看国产| 欧美激情极品国产一区二区三区| 老司机午夜十八禁免费视频| 国产1区2区3区精品| 精品人妻1区二区| 岛国毛片在线播放| 午夜精品久久久久久毛片777| 日本猛色少妇xxxxx猛交久久| 日韩中文字幕欧美一区二区| 欧美激情极品国产一区二区三区| 亚洲av男天堂| 国产精品久久久久久精品电影小说| av在线老鸭窝| 黄片大片在线免费观看| 天堂8中文在线网| 啦啦啦啦在线视频资源| 亚洲中文av在线| 午夜精品久久久久久毛片777| 亚洲 欧美一区二区三区| 午夜精品国产一区二区电影| 99精国产麻豆久久婷婷| 国产成人欧美在线观看 | 在线av久久热| 色婷婷久久久亚洲欧美| 少妇被粗大的猛进出69影院| 一级,二级,三级黄色视频| 大片电影免费在线观看免费| 黄色视频在线播放观看不卡| 久久毛片免费看一区二区三区| 日韩人妻精品一区2区三区| 午夜激情久久久久久久| av免费在线观看网站| 一个人免费看片子| 久久精品国产亚洲av高清一级| 国产成人欧美在线观看 | 国产成人影院久久av| 高清在线国产一区| 狠狠精品人妻久久久久久综合| 国产精品二区激情视频| 久久久国产成人免费| 色视频在线一区二区三区| 妹子高潮喷水视频| 国产欧美亚洲国产| 精品久久蜜臀av无| 国产精品久久久久久精品电影小说| 国产真人三级小视频在线观看| 丝袜美腿诱惑在线| 高清av免费在线| 亚洲国产精品成人久久小说| 久久毛片免费看一区二区三区| 国产黄频视频在线观看| 啦啦啦中文免费视频观看日本| 久久人人97超碰香蕉20202| 久久久国产精品麻豆| 在线亚洲精品国产二区图片欧美| 叶爱在线成人免费视频播放| 高清在线国产一区| 亚洲国产欧美日韩在线播放| 日韩视频一区二区在线观看| 狂野欧美激情性xxxx| 啦啦啦中文免费视频观看日本| 久久久久久久国产电影| 少妇裸体淫交视频免费看高清 | 国产黄色免费在线视频| 一本—道久久a久久精品蜜桃钙片| 免费在线观看日本一区| 精品人妻1区二区| 狠狠狠狠99中文字幕| 国产成人啪精品午夜网站| 又黄又粗又硬又大视频| 亚洲国产看品久久| 这个男人来自地球电影免费观看| 免费av中文字幕在线| 日韩 亚洲 欧美在线| 国产一卡二卡三卡精品| 日韩一卡2卡3卡4卡2021年| 一本久久精品| 女人高潮潮喷娇喘18禁视频| 黄片大片在线免费观看| 精品高清国产在线一区| 一区二区日韩欧美中文字幕| 久久久国产成人免费| 国产老妇伦熟女老妇高清| 久久精品国产亚洲av香蕉五月 | 国产伦理片在线播放av一区| 欧美日韩福利视频一区二区| 97人妻天天添夜夜摸| 18禁国产床啪视频网站| 国产欧美日韩综合在线一区二区| 搡老熟女国产l中国老女人| 亚洲成人国产一区在线观看| 亚洲色图 男人天堂 中文字幕| 人妻人人澡人人爽人人| 欧美激情 高清一区二区三区| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美日韩另类电影网站| 国产一区二区三区在线臀色熟女 | 国产深夜福利视频在线观看| 肉色欧美久久久久久久蜜桃| 久久影院123| 新久久久久国产一级毛片| 国产精品久久久久久人妻精品电影 | 亚洲少妇的诱惑av| 1024视频免费在线观看| 啦啦啦在线免费观看视频4| av线在线观看网站| 国产熟女午夜一区二区三区| 一边摸一边抽搐一进一出视频| 亚洲精品美女久久av网站| 亚洲伊人久久精品综合| 亚洲专区国产一区二区| 国产成人精品无人区| 欧美精品av麻豆av| 巨乳人妻的诱惑在线观看| 国产精品久久久久久人妻精品电影 | 18禁裸乳无遮挡动漫免费视频| 五月天丁香电影| 午夜福利一区二区在线看| 亚洲熟女精品中文字幕| 久久久久网色| 黄频高清免费视频| 国产一区二区 视频在线| 巨乳人妻的诱惑在线观看| 桃红色精品国产亚洲av| 9热在线视频观看99| 精品人妻熟女毛片av久久网站| 欧美中文综合在线视频| 国产深夜福利视频在线观看| 69av精品久久久久久 | 爱豆传媒免费全集在线观看| 12—13女人毛片做爰片一| 亚洲av日韩在线播放| 黄片播放在线免费| 国产日韩欧美亚洲二区| 中文欧美无线码| 人妻人人澡人人爽人人| 欧美老熟妇乱子伦牲交| 欧美日韩亚洲国产一区二区在线观看 | 欧美激情 高清一区二区三区| 久久久久视频综合| 啦啦啦啦在线视频资源| 男女下面插进去视频免费观看| 色播在线永久视频| 欧美老熟妇乱子伦牲交| 在线av久久热| 狠狠精品人妻久久久久久综合| 男女下面插进去视频免费观看| 91老司机精品| 亚洲视频免费观看视频| 最新在线观看一区二区三区| 亚洲欧美激情在线| 国产福利在线免费观看视频| 日本av手机在线免费观看| 成人国产一区最新在线观看| 最黄视频免费看| 亚洲专区字幕在线| www.av在线官网国产| 国产精品亚洲av一区麻豆| 欧美成人午夜精品| 丝袜美腿诱惑在线| 91精品伊人久久大香线蕉| 国产精品久久久av美女十八| 操出白浆在线播放| 侵犯人妻中文字幕一二三四区| 亚洲中文av在线| 50天的宝宝边吃奶边哭怎么回事| 巨乳人妻的诱惑在线观看| 波多野结衣一区麻豆| 亚洲五月色婷婷综合| 五月开心婷婷网| 久久久国产一区二区| 亚洲av成人不卡在线观看播放网 | 午夜老司机福利片| 一二三四在线观看免费中文在| 亚洲精品粉嫩美女一区| 熟女少妇亚洲综合色aaa.| 国产精品影院久久| 亚洲美女黄色视频免费看| 美女国产高潮福利片在线看| videosex国产| 一级毛片女人18水好多| 91精品三级在线观看| 夜夜夜夜夜久久久久| 超碰97精品在线观看| 亚洲av电影在线进入| 国产亚洲午夜精品一区二区久久| 无遮挡黄片免费观看| 一区在线观看完整版| 又黄又粗又硬又大视频| 亚洲中文av在线| 高潮久久久久久久久久久不卡| 国产在线观看jvid| 国产精品免费大片| 女人被躁到高潮嗷嗷叫费观| 制服诱惑二区| 国产精品99久久99久久久不卡| 黄频高清免费视频| 大片免费播放器 马上看| 亚洲精品美女久久久久99蜜臀| 婷婷丁香在线五月| 一个人免费看片子| 在线看a的网站| 日韩免费高清中文字幕av| 人人妻人人爽人人添夜夜欢视频| 久久影院123| 亚洲欧美色中文字幕在线| 后天国语完整版免费观看| 国精品久久久久久国模美| 精品亚洲乱码少妇综合久久| 69精品国产乱码久久久| 亚洲中文字幕日韩| 1024视频免费在线观看| 99国产精品一区二区三区| 香蕉丝袜av| 黄片播放在线免费| 人妻 亚洲 视频| 大码成人一级视频| 狠狠精品人妻久久久久久综合| 成年美女黄网站色视频大全免费| 亚洲欧美一区二区三区久久| 亚洲成av片中文字幕在线观看| 国产一区二区 视频在线| 国产日韩欧美在线精品| 欧美日韩一级在线毛片| 久久久久久免费高清国产稀缺| 交换朋友夫妻互换小说| 欧美精品啪啪一区二区三区 | 中文字幕人妻熟女乱码| 老司机靠b影院| 欧美精品人与动牲交sv欧美| 国产成+人综合+亚洲专区| 女人久久www免费人成看片| av天堂久久9| 在线永久观看黄色视频| 大码成人一级视频| 精品人妻一区二区三区麻豆| 亚洲第一欧美日韩一区二区三区 | 欧美 日韩 精品 国产| 亚洲熟女精品中文字幕| 狂野欧美激情性bbbbbb| 国内毛片毛片毛片毛片毛片| 热re99久久国产66热| 日韩电影二区| 国产成人啪精品午夜网站| 亚洲少妇的诱惑av| 成人国语在线视频| 成人国产一区最新在线观看| 免费观看a级毛片全部| 亚洲精品第二区| 电影成人av| 欧美人与性动交α欧美软件| 日韩,欧美,国产一区二区三区| 色视频在线一区二区三区| 一区二区日韩欧美中文字幕| 国产日韩一区二区三区精品不卡| 国产精品久久久人人做人人爽| 亚洲精品乱久久久久久| 久久精品国产a三级三级三级| 制服诱惑二区| 午夜福利影视在线免费观看| 日韩视频一区二区在线观看| 久久 成人 亚洲| www.999成人在线观看| 可以免费在线观看a视频的电影网站| 久久久久久久国产电影| 在线观看免费日韩欧美大片| 久久久久国产精品人妻一区二区| 嫁个100分男人电影在线观看| 777久久人妻少妇嫩草av网站| 日韩视频一区二区在线观看| 脱女人内裤的视频| 一个人免费在线观看的高清视频 | 老熟女久久久| 国产熟女午夜一区二区三区| 蜜桃国产av成人99| 欧美日韩福利视频一区二区| 蜜桃国产av成人99| 制服诱惑二区| 欧美成人午夜精品| 精品熟女少妇八av免费久了| 午夜视频精品福利| 久久香蕉激情| 又紧又爽又黄一区二区| 嫁个100分男人电影在线观看| 九色亚洲精品在线播放| 亚洲七黄色美女视频| 亚洲av成人一区二区三| 国产一区二区三区在线臀色熟女 | 人人澡人人妻人| 一区二区日韩欧美中文字幕| 可以免费在线观看a视频的电影网站| 国产xxxxx性猛交| 在线av久久热| 无遮挡黄片免费观看| 日本wwww免费看| 久久久国产欧美日韩av| 精品国产一区二区久久| 亚洲 国产 在线| 日本av免费视频播放| 亚洲专区中文字幕在线| 在线观看免费日韩欧美大片| 欧美成人午夜精品| 男女边摸边吃奶| 动漫黄色视频在线观看| 男女边摸边吃奶| av电影中文网址| 岛国毛片在线播放| 日韩免费高清中文字幕av| 亚洲av电影在线进入| 久久久久久久精品精品| 蜜桃在线观看..| 久久国产精品人妻蜜桃| 国产伦理片在线播放av一区| 欧美精品av麻豆av| 在线永久观看黄色视频| 欧美另类一区| 午夜影院在线不卡| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久av美女十八| 少妇的丰满在线观看| 日韩精品免费视频一区二区三区| 中文精品一卡2卡3卡4更新| 99国产极品粉嫩在线观看| 国产国语露脸激情在线看| 在线观看人妻少妇| 老司机午夜十八禁免费视频| 国产深夜福利视频在线观看| 亚洲伊人久久精品综合| 91成人精品电影| 首页视频小说图片口味搜索| 国产主播在线观看一区二区| 男女国产视频网站| 19禁男女啪啪无遮挡网站| av在线播放精品| 香蕉丝袜av| 看免费av毛片| 色94色欧美一区二区| 又黄又粗又硬又大视频| 欧美精品人与动牲交sv欧美| 亚洲精品国产区一区二| 精品国产乱码久久久久久小说| 精品人妻一区二区三区麻豆| 爱豆传媒免费全集在线观看| 精品人妻一区二区三区麻豆| 成年女人毛片免费观看观看9 | 免费在线观看完整版高清| 欧美另类亚洲清纯唯美| 捣出白浆h1v1| 精品少妇久久久久久888优播| 精品国产一区二区三区久久久樱花| 99久久综合免费| 交换朋友夫妻互换小说| 欧美日韩精品网址| 成人国产一区最新在线观看| bbb黄色大片| 国产一区有黄有色的免费视频| 亚洲欧洲精品一区二区精品久久久| 大香蕉久久网| 在线观看舔阴道视频| 中文字幕另类日韩欧美亚洲嫩草| 日韩一区二区三区影片| 亚洲天堂av无毛| 少妇精品久久久久久久| 黄色片一级片一级黄色片| 精品国产乱子伦一区二区三区 | 香蕉国产在线看| 国产精品国产三级国产专区5o| 大香蕉久久成人网| 国产成人欧美| 国产成人免费观看mmmm| 男人添女人高潮全过程视频| av网站在线播放免费| 精品少妇黑人巨大在线播放| 国产av国产精品国产| 热99久久久久精品小说推荐| 日日夜夜操网爽| 日韩视频一区二区在线观看| 亚洲av片天天在线观看| 欧美激情高清一区二区三区| 丝袜美足系列| www.av在线官网国产| 中文字幕av电影在线播放| 国产在线视频一区二区| 亚洲中文日韩欧美视频| 日韩三级视频一区二区三区| 国产一区二区 视频在线| 在线天堂中文资源库| 免费一级毛片在线播放高清视频 | 欧美精品亚洲一区二区| 后天国语完整版免费观看| 亚洲国产欧美日韩在线播放| 热99re8久久精品国产| 欧美+亚洲+日韩+国产| 18禁裸乳无遮挡动漫免费视频| 国产真人三级小视频在线观看| 2018国产大陆天天弄谢| 国产欧美日韩一区二区三 | 成年av动漫网址| 国产成人一区二区三区免费视频网站| 久久人人97超碰香蕉20202| 久久久精品国产亚洲av高清涩受| 久久精品国产亚洲av香蕉五月 | 成人av一区二区三区在线看 | 高清欧美精品videossex| 肉色欧美久久久久久久蜜桃| 性色av一级| 久久性视频一级片| 少妇粗大呻吟视频| 青草久久国产| 人人妻,人人澡人人爽秒播| 狂野欧美激情性xxxx| 亚洲第一av免费看| 多毛熟女@视频| 欧美激情 高清一区二区三区| 日日夜夜操网爽| 亚洲av日韩精品久久久久久密| 精品卡一卡二卡四卡免费| 伊人亚洲综合成人网| 黄网站色视频无遮挡免费观看| 日本av免费视频播放| 免费黄频网站在线观看国产| 欧美国产精品一级二级三级| 十八禁网站免费在线| 我的亚洲天堂| 国产福利在线免费观看视频| 天堂俺去俺来也www色官网| 波多野结衣一区麻豆| 欧美精品高潮呻吟av久久| 午夜久久久在线观看| 欧美激情 高清一区二区三区| 国产亚洲精品一区二区www | 日韩制服丝袜自拍偷拍| 十分钟在线观看高清视频www| 黄片小视频在线播放| 国产成人免费无遮挡视频| 久久女婷五月综合色啪小说| 日韩免费高清中文字幕av| 欧美精品一区二区免费开放| 久久久久久久国产电影| 最新在线观看一区二区三区| 99精品欧美一区二区三区四区| 中文字幕av电影在线播放| 亚洲人成77777在线视频| 亚洲三区欧美一区| 午夜老司机福利片| 亚洲第一av免费看| 亚洲精品久久午夜乱码| 精品久久久久久久毛片微露脸 | 97在线人人人人妻| 夜夜夜夜夜久久久久| 日韩 欧美 亚洲 中文字幕| 久久人人爽av亚洲精品天堂| 91麻豆精品激情在线观看国产 | 手机成人av网站| 高清黄色对白视频在线免费看| 国产精品国产三级国产专区5o| 欧美人与性动交α欧美精品济南到| 侵犯人妻中文字幕一二三四区| 亚洲av男天堂| 久久久精品免费免费高清| 久久午夜综合久久蜜桃| 女人久久www免费人成看片| 在线 av 中文字幕| 操出白浆在线播放| 国产一区有黄有色的免费视频| 男女国产视频网站| 91av网站免费观看| 久久人人爽人人片av| 18禁黄网站禁片午夜丰满| 国产日韩欧美亚洲二区| 18禁裸乳无遮挡动漫免费视频| 国产一卡二卡三卡精品|