• <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| 免费女性裸体啪啪无遮挡网站| 麻豆成人av在线观看| 亚洲成人免费电影在线观看| 男女午夜视频在线观看| av欧美777| 97在线人人人人妻| 精品一区二区三卡| 两个人免费观看高清视频| 91精品三级在线观看| 国产精品国产av在线观看| 好男人电影高清在线观看| 一区二区三区激情视频| av免费在线观看网站| 女人精品久久久久毛片| 丰满饥渴人妻一区二区三| 黄色视频在线播放观看不卡| 亚洲欧美一区二区三区久久| 国产精品1区2区在线观看. | 久久av网站| 91成年电影在线观看| 女性生殖器流出的白浆| 成人免费观看视频高清| 亚洲精品一二三| 精品乱码久久久久久99久播| 久久精品国产a三级三级三级| 一区二区三区乱码不卡18| 久久av网站| 国产精品久久久av美女十八| 中文字幕制服av| 婷婷成人精品国产| 婷婷成人精品国产| 亚洲av成人不卡在线观看播放网| 在线观看免费视频日本深夜| 亚洲色图av天堂| 欧美国产精品va在线观看不卡| 精品一区二区三区av网在线观看 | 大码成人一级视频| 制服人妻中文乱码| 高潮久久久久久久久久久不卡| 香蕉久久夜色| 一边摸一边做爽爽视频免费| 女人高潮潮喷娇喘18禁视频| 欧美午夜高清在线| 人人妻人人澡人人爽人人夜夜| 久久中文看片网| 12—13女人毛片做爰片一| 亚洲专区国产一区二区| 日本黄色视频三级网站网址 | 一区二区三区乱码不卡18| 日韩欧美国产一区二区入口| 久久国产精品影院| 中亚洲国语对白在线视频| av不卡在线播放| 成人永久免费在线观看视频 | 超碰97精品在线观看| 国产精品99久久99久久久不卡| 一级黄色大片毛片| 久久久精品区二区三区| 国产精品国产高清国产av | 一级a爱视频在线免费观看| 一边摸一边抽搐一进一出视频| 亚洲国产成人一精品久久久| 国产日韩欧美在线精品| 久久久国产精品麻豆| 香蕉丝袜av| 欧美精品一区二区免费开放| 色综合欧美亚洲国产小说| 日本av手机在线免费观看| 中文字幕av电影在线播放| 亚洲 欧美一区二区三区| 欧美国产精品va在线观看不卡| 香蕉久久夜色| 超碰97精品在线观看| 久久国产精品影院| 丝瓜视频免费看黄片| 日韩精品免费视频一区二区三区| 久久久久精品人妻al黑| 99精国产麻豆久久婷婷| 成人手机av| 日韩一卡2卡3卡4卡2021年| 制服人妻中文乱码| 午夜福利影视在线免费观看| 热re99久久精品国产66热6| 极品少妇高潮喷水抽搐| 久久中文看片网| 最近最新中文字幕大全免费视频| 麻豆国产av国片精品| 欧美另类亚洲清纯唯美| 91精品国产国语对白视频| 亚洲欧美日韩另类电影网站| 国产在线一区二区三区精| 飞空精品影院首页| 亚洲视频免费观看视频| 精品午夜福利视频在线观看一区 | 国产精品 欧美亚洲| 人妻久久中文字幕网| 国产男女内射视频| 国产精品 欧美亚洲| 久久久久国产一级毛片高清牌| 久久久久国产一级毛片高清牌| 三级毛片av免费| 日本wwww免费看| 亚洲av国产av综合av卡| 亚洲欧美一区二区三区黑人| 日本a在线网址| 女人久久www免费人成看片| 日韩成人在线观看一区二区三区| 下体分泌物呈黄色| 一本综合久久免费| 啦啦啦在线免费观看视频4| 激情在线观看视频在线高清 | 在线观看免费视频日本深夜| 欧美激情久久久久久爽电影 | 在线观看免费视频网站a站| 91麻豆精品激情在线观看国产 | 国产成人免费无遮挡视频| 黑人操中国人逼视频| 97人妻天天添夜夜摸| 女同久久另类99精品国产91| 王馨瑶露胸无遮挡在线观看| 2018国产大陆天天弄谢| 国产成人啪精品午夜网站| 亚洲精品国产一区二区精华液| 国产高清videossex| 国产高清videossex| 亚洲国产欧美在线一区| 正在播放国产对白刺激| 亚洲成人手机| 一本综合久久免费| 大片免费播放器 马上看| 在线观看舔阴道视频| 热99国产精品久久久久久7| 水蜜桃什么品种好| 婷婷成人精品国产| 日韩成人在线观看一区二区三区| 丁香欧美五月| 两个人免费观看高清视频| 中文字幕人妻丝袜一区二区| 19禁男女啪啪无遮挡网站| 亚洲七黄色美女视频| 久久精品国产亚洲av高清一级| 99国产精品一区二区蜜桃av | 十八禁高潮呻吟视频| 国产精品免费视频内射| 两个人看的免费小视频| 精品国产国语对白av| 色播在线永久视频| 亚洲国产中文字幕在线视频| 国产成人免费无遮挡视频| 窝窝影院91人妻| 国产精品98久久久久久宅男小说| 黄色视频在线播放观看不卡| 日韩免费高清中文字幕av| 欧美国产精品va在线观看不卡| 亚洲五月色婷婷综合| videos熟女内射| 宅男免费午夜| 国产三级黄色录像| av超薄肉色丝袜交足视频| av片东京热男人的天堂| 国产精品一区二区在线不卡| 久久国产亚洲av麻豆专区| 丝瓜视频免费看黄片| 欧美av亚洲av综合av国产av| 高清视频免费观看一区二区| 久久亚洲真实| 一级毛片电影观看| 91九色精品人成在线观看| 国产在视频线精品| 国产男靠女视频免费网站| 99在线人妻在线中文字幕 | 最近最新中文字幕大全免费视频| 色婷婷av一区二区三区视频| 亚洲色图综合在线观看| 亚洲精品中文字幕在线视频| 国产主播在线观看一区二区| 人成视频在线观看免费观看| 高清av免费在线| av免费在线观看网站| 考比视频在线观看| cao死你这个sao货| 高清av免费在线| 亚洲色图av天堂| 十八禁高潮呻吟视频| 国产av精品麻豆| 波多野结衣一区麻豆| 人人澡人人妻人| 国产一区二区三区视频了| 老汉色av国产亚洲站长工具| 国产黄频视频在线观看| 久久国产精品大桥未久av| 日韩人妻精品一区2区三区| 中文字幕色久视频| www日本在线高清视频| 欧美在线黄色| 天堂中文最新版在线下载| 国产成人欧美在线观看 | av免费在线观看网站| 国产免费av片在线观看野外av| 国产成人精品无人区| 91麻豆av在线| 黄色视频在线播放观看不卡| 另类亚洲欧美激情| 麻豆av在线久日| 丰满人妻熟妇乱又伦精品不卡| 99国产精品一区二区三区| h视频一区二区三区| 中国美女看黄片| 淫妇啪啪啪对白视频| 亚洲国产欧美日韩在线播放| 女性生殖器流出的白浆| 丝袜人妻中文字幕| 日本av手机在线免费观看| 纯流量卡能插随身wifi吗| 国产亚洲欧美在线一区二区| 热99re8久久精品国产| 国产在视频线精品| 久久人人97超碰香蕉20202| 蜜桃在线观看..| 91国产中文字幕| 日日摸夜夜添夜夜添小说| 欧美精品高潮呻吟av久久| 男女下面插进去视频免费观看| 法律面前人人平等表现在哪些方面| 国产亚洲av高清不卡| 少妇被粗大的猛进出69影院| 亚洲少妇的诱惑av| 国产有黄有色有爽视频| 777米奇影视久久| 日韩免费高清中文字幕av| 少妇精品久久久久久久| 99精国产麻豆久久婷婷| 亚洲国产av影院在线观看| 国产亚洲欧美精品永久| 色综合欧美亚洲国产小说| 亚洲全国av大片| 久久精品人人爽人人爽视色| 日韩有码中文字幕| 婷婷丁香在线五月| 狂野欧美激情性xxxx| 精品高清国产在线一区| 国产在线观看jvid| 久久性视频一级片| 一区二区三区国产精品乱码| av线在线观看网站| 满18在线观看网站| svipshipincom国产片| 美女国产高潮福利片在线看| 在线观看一区二区三区激情| 亚洲精品美女久久av网站| 国产麻豆69| 精品国产国语对白av| 国产日韩欧美在线精品| 中亚洲国语对白在线视频| 国产成人av激情在线播放| 悠悠久久av| 国产成人精品在线电影| 777米奇影视久久| 日韩一区二区三区影片| 天堂动漫精品| 女警被强在线播放| 亚洲自偷自拍图片 自拍| 另类亚洲欧美激情| 久久九九热精品免费| 久久人人爽av亚洲精品天堂| 黑人巨大精品欧美一区二区蜜桃| 国产福利在线免费观看视频| 久久久久精品国产欧美久久久| 12—13女人毛片做爰片一| 黑丝袜美女国产一区| 日韩中文字幕视频在线看片| 欧美成人午夜精品| 日韩三级视频一区二区三区| 可以免费在线观看a视频的电影网站| 午夜两性在线视频| 久久精品人人爽人人爽视色| 亚洲第一欧美日韩一区二区三区 | 国产一区二区 视频在线| 超碰成人久久| av在线播放免费不卡| 露出奶头的视频| av一本久久久久| 999久久久国产精品视频| 天堂中文最新版在线下载| 久久国产精品人妻蜜桃| 伊人久久大香线蕉亚洲五| 男女免费视频国产| 高潮久久久久久久久久久不卡| 男人舔女人的私密视频| 久久国产精品大桥未久av| 色婷婷av一区二区三区视频| 午夜激情av网站| 女人久久www免费人成看片| 性少妇av在线| 搡老熟女国产l中国老女人| 国产不卡一卡二| 久久性视频一级片| 色综合婷婷激情| 精品视频人人做人人爽| 国产精品偷伦视频观看了| xxxhd国产人妻xxx| 国产在线观看jvid| 一边摸一边抽搐一进一小说 | 欧美日韩av久久| 国产成人精品在线电影| 亚洲成人手机| 亚洲av电影在线进入| 亚洲情色 制服丝袜| 国产免费av片在线观看野外av| 老司机福利观看| 97在线人人人人妻| e午夜精品久久久久久久| 亚洲全国av大片| 久9热在线精品视频| 日韩中文字幕欧美一区二区| 亚洲国产精品一区二区三区在线| 欧美精品啪啪一区二区三区| 最近最新中文字幕大全免费视频| 757午夜福利合集在线观看| 亚洲九九香蕉| 一区二区三区乱码不卡18| 大香蕉久久网| 99国产精品99久久久久| 90打野战视频偷拍视频| 亚洲午夜理论影院| 国产av精品麻豆| 午夜日韩欧美国产| 成人亚洲精品一区在线观看| 国产1区2区3区精品| 天天影视国产精品| 高清av免费在线| 亚洲专区中文字幕在线| 午夜老司机福利片| 国产欧美日韩一区二区三区在线| 国产精品美女特级片免费视频播放器 | 老司机亚洲免费影院| 亚洲性夜色夜夜综合| 精品高清国产在线一区| 99精品欧美一区二区三区四区| 一级,二级,三级黄色视频| 十分钟在线观看高清视频www| 99国产精品99久久久久| 老熟女久久久| 色精品久久人妻99蜜桃| 精品少妇久久久久久888优播| 国产亚洲精品第一综合不卡| 亚洲一区二区三区欧美精品| 国产免费视频播放在线视频| 久久狼人影院| 69精品国产乱码久久久| 欧美性长视频在线观看| 热re99久久精品国产66热6| 高清av免费在线| 日本五十路高清| 国产欧美亚洲国产| 黑人猛操日本美女一级片| 国产亚洲av高清不卡| 久久中文字幕一级| 国产精品免费大片| 午夜精品国产一区二区电影| 少妇裸体淫交视频免费看高清 | 成人国产av品久久久| 99国产精品免费福利视频| 菩萨蛮人人尽说江南好唐韦庄| 水蜜桃什么品种好| 99riav亚洲国产免费| 黑人巨大精品欧美一区二区mp4| 母亲3免费完整高清在线观看| 9热在线视频观看99| 在线看a的网站| 在线 av 中文字幕| 麻豆av在线久日| 久久国产精品影院| 欧美精品一区二区免费开放| 国产精品久久久久久精品古装| 成年人黄色毛片网站| 久久婷婷成人综合色麻豆| 国产欧美日韩一区二区精品| 成人国产一区最新在线观看| 两人在一起打扑克的视频| 在线观看人妻少妇| 91麻豆av在线| 国产精品国产av在线观看| 搡老乐熟女国产| 久久中文看片网| 国产aⅴ精品一区二区三区波| kizo精华| 黄色怎么调成土黄色| 99在线人妻在线中文字幕 | 欧美+亚洲+日韩+国产| 亚洲国产av影院在线观看| 免费看a级黄色片| 日韩大片免费观看网站| av欧美777| 国产片内射在线| 1024香蕉在线观看| 亚洲 欧美一区二区三区| 国产精品.久久久| 精品国产国语对白av| 叶爱在线成人免费视频播放| 免费日韩欧美在线观看| 一级片免费观看大全| 亚洲精品成人av观看孕妇| 高清毛片免费观看视频网站 | 亚洲 欧美一区二区三区| 女人久久www免费人成看片| 午夜日韩欧美国产| 满18在线观看网站| 一边摸一边做爽爽视频免费| 中文字幕人妻熟女乱码| 男女午夜视频在线观看| 精品久久蜜臀av无| 一进一出好大好爽视频| 亚洲 欧美一区二区三区| 90打野战视频偷拍视频| 国产在视频线精品| av超薄肉色丝袜交足视频| 欧美成人午夜精品| 免费在线观看黄色视频的| 极品少妇高潮喷水抽搐| 国产精品久久电影中文字幕 | 国产午夜精品久久久久久| 亚洲成人手机| 满18在线观看网站| 日韩中文字幕欧美一区二区| 视频区图区小说| 午夜福利视频在线观看免费| 国产精品久久电影中文字幕 | 午夜福利在线免费观看网站| 亚洲va日本ⅴa欧美va伊人久久| 99国产精品一区二区蜜桃av | 欧美黄色片欧美黄色片| 少妇猛男粗大的猛烈进出视频| 免费高清在线观看日韩| a级毛片黄视频| 十八禁人妻一区二区| 每晚都被弄得嗷嗷叫到高潮| 夜夜夜夜夜久久久久| 国产不卡一卡二| 国产aⅴ精品一区二区三区波| 成年人免费黄色播放视频| 亚洲中文日韩欧美视频| 久久人妻福利社区极品人妻图片| 九色亚洲精品在线播放| 后天国语完整版免费观看| 老熟女久久久| 2018国产大陆天天弄谢| 亚洲第一青青草原| 久久ye,这里只有精品| 一二三四社区在线视频社区8| 色播在线永久视频| 亚洲av欧美aⅴ国产| 动漫黄色视频在线观看| 极品人妻少妇av视频| 日韩 欧美 亚洲 中文字幕| 男女高潮啪啪啪动态图| 高潮久久久久久久久久久不卡| 成年人免费黄色播放视频| 美国免费a级毛片| 精品福利永久在线观看| 美女视频免费永久观看网站| 最新在线观看一区二区三区| 国产精品久久久人人做人人爽| 青青草视频在线视频观看| 久久精品人人爽人人爽视色| 国产精品久久久久久精品电影小说| 国产黄色免费在线视频| 亚洲 国产 在线| 欧美一级毛片孕妇| 深夜精品福利| 91九色精品人成在线观看| av网站免费在线观看视频| 日本黄色视频三级网站网址 | 91九色精品人成在线观看| 久久亚洲真实| 纯流量卡能插随身wifi吗| 国产成人免费观看mmmm| 久久久欧美国产精品| 脱女人内裤的视频| 日本a在线网址| 日韩中文字幕欧美一区二区| aaaaa片日本免费| 久久亚洲精品不卡| 中文字幕av电影在线播放| 性色av乱码一区二区三区2| 一二三四在线观看免费中文在| 亚洲一卡2卡3卡4卡5卡精品中文| 精品亚洲成国产av| av一本久久久久| 在线观看免费视频网站a站| 高清av免费在线| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲色图综合在线观看| 天堂中文最新版在线下载| 午夜两性在线视频| 日本vs欧美在线观看视频| 精品人妻在线不人妻| 国产av国产精品国产| 18禁美女被吸乳视频| 精品国产亚洲在线| 操美女的视频在线观看| 80岁老熟妇乱子伦牲交| 久久国产亚洲av麻豆专区| 高清黄色对白视频在线免费看| 精品一区二区三区视频在线观看免费 | 天天躁狠狠躁夜夜躁狠狠躁| av一本久久久久| 日韩熟女老妇一区二区性免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日韩成人在线一区二区| 色婷婷久久久亚洲欧美| 肉色欧美久久久久久久蜜桃| 国产亚洲精品一区二区www | 亚洲av欧美aⅴ国产| 久久 成人 亚洲| 水蜜桃什么品种好| 日韩视频一区二区在线观看| 在线观看免费高清a一片| 91精品三级在线观看| 黑人欧美特级aaaaaa片| 欧美变态另类bdsm刘玥| 午夜久久久在线观看| 人妻 亚洲 视频| 国产精品 欧美亚洲| 欧美成狂野欧美在线观看| 精品国产超薄肉色丝袜足j| 男人操女人黄网站| 午夜视频精品福利| 日韩人妻精品一区2区三区| 精品久久蜜臀av无| 国产又色又爽无遮挡免费看| 国产精品欧美亚洲77777| 少妇被粗大的猛进出69影院| 最新美女视频免费是黄的| 99re6热这里在线精品视频| 人人妻人人添人人爽欧美一区卜| www.999成人在线观看| 久久国产亚洲av麻豆专区| 90打野战视频偷拍视频| 香蕉久久夜色| 人成视频在线观看免费观看| 日韩一区二区三区影片| 欧美人与性动交α欧美精品济南到| 久久久久久人人人人人| 视频区图区小说| 天堂俺去俺来也www色官网| 国产成人影院久久av| 色老头精品视频在线观看| 亚洲国产欧美一区二区综合| 亚洲成a人片在线一区二区| 18禁观看日本| 99热网站在线观看| 涩涩av久久男人的天堂| 午夜激情久久久久久久| h视频一区二区三区| 中文字幕最新亚洲高清| 亚洲精品美女久久av网站| 9191精品国产免费久久| 亚洲av第一区精品v没综合| 一二三四在线观看免费中文在| 久久精品亚洲av国产电影网| 91字幕亚洲| 亚洲免费av在线视频| tube8黄色片| 国产精品 欧美亚洲| www.999成人在线观看| 精品久久久精品久久久| 国产精品av久久久久免费| 我要看黄色一级片免费的| 最近最新免费中文字幕在线| 每晚都被弄得嗷嗷叫到高潮| 亚洲天堂av无毛| 丝瓜视频免费看黄片| 国产精品亚洲av一区麻豆| 考比视频在线观看| 国产aⅴ精品一区二区三区波| 中文字幕人妻丝袜一区二区| 国产激情久久老熟女| 在线观看舔阴道视频| 精品国产一区二区三区久久久樱花| 日韩免费高清中文字幕av| av网站免费在线观看视频| 亚洲中文字幕日韩| 一区在线观看完整版| 亚洲成人免费av在线播放| 多毛熟女@视频| 亚洲精品粉嫩美女一区| 99国产精品一区二区三区| 在线观看免费日韩欧美大片| 人人妻,人人澡人人爽秒播| 夜夜夜夜夜久久久久| 少妇被粗大的猛进出69影院| 久久国产亚洲av麻豆专区| 久久久久久久精品吃奶| 老司机靠b影院| 老司机午夜福利在线观看视频 | 亚洲精品中文字幕一二三四区 | 成年人黄色毛片网站| 国产男靠女视频免费网站| 国产精品 欧美亚洲| 久久久水蜜桃国产精品网| 欧美性长视频在线观看| 成年人黄色毛片网站| 日韩成人在线观看一区二区三区| 亚洲欧洲精品一区二区精品久久久| 欧美日韩中文字幕国产精品一区二区三区 | 久久热在线av| 国产男女超爽视频在线观看| 亚洲欧美精品综合一区二区三区|