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

    未來30 年亞洲降水情景預(yù)估及偏差訂正

    2022-01-26 06:21:54楊陽戴新剛汪萍
    大氣科學(xué) 2022年1期
    關(guān)鍵詞:距平對數(shù)增量

    楊陽 戴新剛 汪萍

    1 中國科學(xué)院東亞區(qū)域氣候與環(huán)境重點實驗室,大氣物理研究所,北京 100029

    2 中山大學(xué)大氣科學(xué)學(xué)院,珠海 519082

    3 中國氣象科學(xué)研究院大氣成分與環(huán)境氣象研究所,北京 100081

    4 蘭州大學(xué)大氣科學(xué)學(xué)院,蘭州 730000

    1 引言

    季風(fēng)氣候區(qū)是地球上植被最繁茂、人口最密集、耕地最多、農(nóng)業(yè)生產(chǎn)最活躍的區(qū)域(Webster,1981; Dai and Wang, 2017; Chang et al., 2018; 戴新剛和汪萍, 2020; Wang et al., 2021)。亞洲南部和東部是地球上最著名的季風(fēng)區(qū)之一,其中南亞屬于熱帶季風(fēng)區(qū),東亞季風(fēng)屬于副熱帶季風(fēng)(Dai et al.,2004b; Ding and Chan, 2005; Wang, 2006; 朱志偉和何金海, 2013)。南亞地區(qū)的旱澇主要與季風(fēng)的濤動有關(guān),在強(qiáng)季風(fēng)年,南亞降水偏多,洪澇災(zāi)害頻發(fā);弱季風(fēng)年南亞降水偏少,部分地區(qū)干旱,影響作物收成(Revadekar and Preethi, 2012; Lu et al.,2020)。對于東亞而言,在強(qiáng)季風(fēng)年西太平洋副熱帶高壓(簡稱西太副高)往往偏北偏東,汛期中國東南沿海及黃淮流域降水偏多,長江中下游地區(qū)降水偏少;弱季風(fēng)年西太副高偏南偏西,長江及江南等地多雨(丁一匯, 1993; 陶詩言, 1980; 陶詩言等,2001),而北方干旱(Dai et al., 2003, 2005),造成糧食產(chǎn)量下降(Terjung et al., 1989)。影響亞洲季風(fēng)降水的因子眾多,其中ENSO(El Ni?o-Southern Oscillation)是亞洲季風(fēng)年際變化的主要驅(qū)動因子( Krishnamurthy and Goswami, 2000;Fasullo and Webster, 2002)。這種影響既可以是間接的,也可以是直接的,受海溫年代尺度/多年代尺度遙相關(guān)型調(diào)制(Dai et al., 2004a; Wang et al.,2013; Zhu et al., 2014; Dong et al., 2018),如太平洋年代尺度濤動(PDO/IPO)(Feng et al., 2014;Dai et al., 2015a; Dong and Dai, 2015),北大西洋多年代尺度濤動(AMO; Naidu et al., 2020)等。此外,全球季風(fēng)系統(tǒng)也會受到全球氣候暖化趨勢的明顯影響(Wang et al., 2015)。

    通常,El Ni?o 年亞洲季風(fēng)偏弱、西太平洋副熱帶高壓偏南偏西、東亞季風(fēng)雨帶隨之偏南,汛期中國大陸易出現(xiàn)“南澇北旱”格局,冬季氣溫偏高,且易出現(xiàn)暖冬;反之,在La Ni?a 年亞洲季風(fēng)偏強(qiáng),西太平洋副熱帶高壓偏北偏東,汛期易出現(xiàn)所謂“Ⅰ類雨型”,即中國大陸南北降水多,中間少的格局(符淙斌, 1978; 黃榮輝等, 2003; 吳國雄等,2003; Wu et al., 2003; 王會軍等, 2020)。在PDO的正位相,北半球中緯度太平洋東部海溫偏高,西北太平洋北部海溫偏低,亞洲季風(fēng)偏弱,中國大陸易出現(xiàn)年代尺度“南澇北旱”格局(Yang et al.,2017a, 2017b);與此相反,PDO 的負(fù)位相使得亞洲季風(fēng)偏強(qiáng),冬季氣溫偏低,易出現(xiàn)冷冬(Wang et al., 2008; Chen et al., 2013; Yun et al., 2015; Huang et al., 2020a)。已有數(shù)值試驗證實,PDO 的位相變化還能明顯影響印度洋異常海溫分布,從而對南亞及東亞季風(fēng)均產(chǎn)生影響(Krishnamurthy and Krishnamurthy, 2014; Wang et al., 2014)。在AMO正位相,北大西洋海溫偏暖,負(fù)位相海溫偏冷。在多年代尺度上,它不僅能影響中緯度西風(fēng)帶的濕度和降水,也與亞洲季風(fēng)存在關(guān)聯(lián)(Goswami et al.,2006; Krishnamurthy and Krishnamurthy, 2016; Li et al., 2017)。

    此外,過去百年的全球氣候暖化已經(jīng)逐漸影響到地球氣候系統(tǒng)的多個方面(IPCC, 2013; Dai and Wang, 2018),例如,北極海冰減少(Kumar et al., 2020)、山岳冰川融化(張宏杰等, 2018)、中太平洋暖事件增多(Ashok et al., 2007)、副熱帶高壓面積擴(kuò)展(Seidel et al., 2008)等。氣候系統(tǒng)熱機(jī)效率的提高導(dǎo)致全球能量與水循環(huán)速率加快(Laliberté et al., 2015),致使極端天氣氣候事件的出現(xiàn)頻率呈現(xiàn)上升趨勢(Li et al., 2020; Zhong et al., 2020)。氣候系統(tǒng)的這一系列變化已經(jīng)影響到區(qū)域尺度氣候及亞洲季風(fēng)氣候(Sooraj et al., 2015;Duan et al., 2018; Cai et al., 2021)。這增加了氣候變化機(jī)制的復(fù)雜性,并且造成預(yù)測難度上升。因此,研究或預(yù)估未來氣候暖化背景下幾十年或百年以上區(qū)域尺度氣候的可能變化趨勢,必須得考慮未來社會經(jīng)濟(jì)發(fā)展情景(IPCC, 2013)。

    氣候變化趨勢是氣候年際變化及年代際變化的背景。前者往往屬于外部強(qiáng)迫的結(jié)果,如全球氣候暖化等,后者屬于氣候系統(tǒng)的內(nèi)部變率,如?!獨庀嗷プ饔玫龋℉uang et al., 2020b)。未來全球氣候暖化趨勢在很大程度上取決于溫室氣體的排放情景(Dai et al., 2015b; Li et al., 2020)。目前,科學(xué)家對于未來幾十年或百年氣候趨勢預(yù)估都基于不同溫室氣體排放情景(RCPs/SSPs)下的國際耦合模式比較計劃第五/六階段(CMIP5/CMIP6)的多模式輸出數(shù)據(jù)(Taylor et al., 2012; Eyring et al., 2016)。例如,Lee and Wang (2014)基于CMIP5 模式資料對未來全球季風(fēng)變化進(jìn)行預(yù)估; Cook et al.(2020)分析了CMIP6 模式輸出的預(yù)估資料后指出,21 世紀(jì)全球多個地區(qū)可能存在干旱趨勢,如美國西部、中美洲、歐洲和地中海、亞馬遜、南部非洲、中國部分地區(qū)、東南亞及澳大利亞等,其中與高排放情景相比,低排放情景下某些地區(qū)的嚴(yán)重干旱會減弱甚至消失。這說明氣候暖化會在一定程度上改變?nèi)蚝禎撤植几窬?,盡管區(qū)域性降水還會受到年際或年代際氣候變化的影響,且與多個海區(qū)的海溫遙相關(guān)型有關(guān)。

    Xu et al.(2018)利用CMIP5 數(shù)據(jù)對未來50年中國極端氣溫和降水進(jìn)行了預(yù)估,指出區(qū)域性極端氣候事件出現(xiàn)的潛在風(fēng)險。Jiang et al.(2020)利用CMIP6 多模式集合數(shù)據(jù)預(yù)估了未來百年中亞干旱區(qū)降水的變化,發(fā)現(xiàn)在不同排放情景下中亞降水幾乎都增加,且在高排放情景(SSP5-8.5)下,降水增加更加明顯。例如,在中亞北部表現(xiàn)為夏季降水減少,春季降水明顯增加。張蓓等(2019)、楊陽等(2019)利用CMIP5 多模式集合數(shù)據(jù)對中國降水的預(yù)估研究表明,2016~2035/2045 年間,從中國西部的三江源區(qū)到淮河流域一帶的降水呈增加趨勢,而西南地區(qū)可能會持續(xù)干旱。

    然而,一些預(yù)估檢驗也指出,CMIP5 模式輸出數(shù)據(jù)存在系統(tǒng)性偏差,尤其是降水,因此氣候預(yù)估時必須進(jìn)行偏差訂正(張蓓和戴新剛, 2016)。模式的氣候漂移是氣候預(yù)估偏差的主要部分,大多數(shù)氣候預(yù)估研究都選擇通過兩個時段平均之差或只計算趨勢的方式消除或回避模式集合的氣候漂移(Chen et al., 2020)。此外,模式降水偏差還存在非定常振蕩,會對某些地區(qū)的幾十年降水預(yù)估產(chǎn)生明顯影響,對此也需要訂正。張蓓等(2019)利用線性回歸方程法訂正了2016~2035 年中國大陸情景預(yù)估降水的定常及非定常偏差。楊陽等(2019)研究發(fā)現(xiàn),回歸方程偏差訂正效果還與緯度有關(guān),提出所謂“區(qū)域組合回歸”訂正法,即在30°N 以北地區(qū)使用一種回歸方案,以南地區(qū)使用另一種,如此可以得到更高的預(yù)估降水距平符號一致率。

    雖然已有一系列的研究用CMIP5/CMIP6 對未來百年或21 世紀(jì)后三十年的降水進(jìn)行預(yù)估,但幾乎均未考慮其中的非定常偏差訂正問題(Li et al.,2019; Cook et al., 2020; Jiang et al., 2020; Wang et al.,2020),增加了預(yù)估結(jié)果的不確定性。本文基于CMIP5 多模式集合輸出資料,在中等排放情景(RCP4.5)下,對未來30 年(2021~2050)亞洲降水進(jìn)行預(yù)估并對其做定常和非定常偏差訂正,給出未來30 年亞洲可能出現(xiàn)的旱澇分布格局或變化趨勢。下文的第二節(jié)介紹了所用數(shù)據(jù)和方法;第三節(jié)是CMIP5 歷史氣候模擬的偏差分析;第四節(jié)對降水預(yù)估偏差進(jìn)行訂正檢驗,并建立統(tǒng)計回歸訂正模型;第五節(jié)給出未來30 年亞洲降水預(yù)估結(jié)果;最后一節(jié)是總結(jié)和討論。

    2 數(shù)據(jù)和方法

    本文使用耦合模式比較計劃第五階段(CMIP5)中多模式輸出的1960~2005 年歷史氣候模擬及未來RCP4.5 排放情景下(2006~2015 及2021~2050 年)的逐月降水?dāng)?shù)據(jù)(Taylor et al., 2012),以及英國氣候研究所全球陸地格點分析數(shù)據(jù)集(CRU TS v4.0)中的格點月降水量(Harris et al.,2020),其分辨率為0.5°×0.5°,視為觀測降水??紤]到過多模式的集合平均可能會使結(jié)果更接近于氣候的強(qiáng)迫場結(jié)果(Dai et al., 2015a),而模式過少又會造成降水集合平均的統(tǒng)計穩(wěn)定性變差。另一方面,在降水演變中的趨勢項一般都很弱,其年代際變化就顯得更為重要,這就要求在選取模式數(shù)量上既不能太多,也不能太少。本文選用了20 個常用的CMIP5 模式,所選模式參數(shù)見表1。

    由于CMIP5 各模式分辨率不一致,本文使用雙線性插值法將每個模式的降水和觀測格點分析降水都統(tǒng)一插值到1.0°×1.0°經(jīng)緯網(wǎng)格點上再進(jìn)行計算。本文使用三種模式偏差訂正法,第一種是去除模式氣候漂移,即通過去除模式和觀測降水的平均偏差來實現(xiàn);后兩種分別為一元對數(shù)回歸和一元年際增量回歸(或一元差分回歸)(Fan et al., 2009;張凱靜等, 2017)。根據(jù)這三種偏差訂正法在地理區(qū)域上的優(yōu)劣又進(jìn)一步可以設(shè)計區(qū)域組合回歸偏差訂正法,即在30°N 以南/北使用一元年際增量或一元對數(shù)回歸法對偏差進(jìn)行訂正(楊陽等, 2019),并對檢驗期中的盲區(qū)(距平符號相反區(qū))再用去除模式氣候漂移的結(jié)果進(jìn)行補(bǔ)充訂正填補(bǔ)。兩種回歸模型參數(shù)均以1960~2005 年為訓(xùn)練期,通過20個CMIP5 模式集合平均降水和觀測降水(CRU TS v4.0)確定。

    3 歷史氣候降水模擬

    CMIP5 歷史氣候模擬中的多模式集合平均對亞洲降水氣候分布型具有一定的模擬能力(楊陽等, 2019)。圖1a 是20 個CMIP5 模式集合平均模擬的1960~2005 年亞洲區(qū)域的年降水(簡稱為模式降水,下同)。與觀測降水(圖1b)比較可以看到,模式集合能基本再現(xiàn)1960~2005 年亞洲氣候平均年降水分布的主要特征。例如,在亞洲季風(fēng)區(qū),包括南亞、東亞及東南亞等地降水較多,中國西北、中亞和阿拉伯半島等干旱氣候區(qū)降水稀少等。然而,若仔細(xì)比較二者可以發(fā)現(xiàn)模式降水依舊存在偏差。例如,模式降水在北方偏多,致使干旱或半干旱區(qū)的范圍縮小,這在中國西北和中亞最為明顯,而模式降水在青藏高原東部明顯偏多,但在印度次大陸偏少。

    圖1 1960~2005 年亞洲氣候平均年降水量(單位:mm):(a) 20 個CMIP5 模式集合模擬;(b)觀測(CRU TS v4.0)Fig. 1 Climate mean annual precipitation (units: mm) from 1960 to 2005: (a) Climate mean simulated using 20 CMIP5 models; (b) observations(CRU TS v4.0)

    暖季(5~10 月)和冷季(11 月至次年4 月)的模式降水偏差分布不同。用模式降水減去觀測降水得到模式降水的偏差,對其再除以觀測降水便可以得到偏差百分率。圖2a 表示暖季即亞洲南部季風(fēng)季節(jié)的模式降水偏差百分率,可見20 個CMIP5模式集合對亞洲季風(fēng)區(qū)大部分地區(qū)降水模擬偏弱,如中國東南沿海和華南沿海、南亞、朝鮮半島、日本關(guān)西地區(qū)、越南大部等地,其中模式降水在印度半島中北部偏少超過40%,在越南中北部、緬甸西南部、朝鮮半島和日本關(guān)西地區(qū)偏少約30%,在中國華南沿海偏少10%~20%,在我國東南沿海偏少5%~10%。相反,在中國西部、北部,中亞、遠(yuǎn)東和阿拉伯半島等地模式降水明顯偏多,其中青藏高原、華北、中亞南部、阿拉伯半島大部模式降水偏多40%以上,在中國東北、中亞北部、印度南部、中南半島西南部模式降水偏少20%~30%。

    在冷季,亞洲大部分地區(qū)模式降水偏多且超過40%,尤其在高緯度地區(qū)降水偏多更為明顯(圖2b)。與模式集合模擬的夏季風(fēng)降水偏少相反,冷季印度半島大部分地區(qū)的模式降水偏多,北部偏多30%~40%,僅在中國華南沿海、中南半島、中亞南部、阿拉伯半島等地模式降水偏少,其中在華南沿海和中南半島偏少10%~20%,在中亞南部及伊朗高原偏少20%~30%,在阿拉伯半島東南部偏少超過40%。可見CMIP5 歷史氣候模擬降水中存在明顯偏差,在使用前必須做不確定性分析并進(jìn)行偏差訂正。

    4 模式降水預(yù)估偏差

    20 個CMIP5 模式集合平均預(yù)估的2006~2015 年亞洲暖季降水偏差分布型(圖3a)與1960~2005年歷史氣候模擬降水偏差分布相似(圖2a)。例如,預(yù)估降水在東亞北部偏多,其中在中國西部、蒙古國及中亞北部等地偏多達(dá)40%及以上;在亞洲低緯度及阿拉伯半島等地降水偏少,其中印度中北部和孟加拉國大部偏少40%及以上,巴基斯坦、阿拉伯半島部分地區(qū)、越南中北部、中國的華南沿海和朝鮮半島等地偏少20%~30%。同樣,模式集合預(yù)估的2006~2015 年亞洲冷季降水偏差分布型(圖3b)也與1960~2005 年平均歷史氣候模擬亞洲冷季降水偏差(圖2b)相似。這種相似性說明二者偏差分布型幾乎不隨時間變化,說明CMIP5氣候模式集合存在系統(tǒng)性偏差即模式氣候漂移。張蓓和戴新剛(2016)對CMIP5 模式集合偏差的分析表明,模式降水的氣候漂移約占總偏差的80%或以上??梢?006~2015 年平均降水預(yù)估偏差的絕大部分應(yīng)該屬于模式集合氣候漂移,具有準(zhǔn)定常性。因此,在CMIP5 模式集合的氣候預(yù)估中直接扣除模式氣候漂移可以減少降水預(yù)估的不確定性。通常的做法是用模式集合預(yù)估降水減去過去一段歷史氣候模擬降水來消除模式氣候漂移。這里假設(shè)模式氣候漂移是定常的。然而,比較圖3 和圖2 中偏差分布的細(xì)節(jié)就會發(fā)現(xiàn),二者之間還存在一些小的差異。例如,在中國東部沿海地區(qū)預(yù)估的降水偏少率幅度要大于歷史氣候模擬;在我國的華北和東北、蒙古國、中亞等地預(yù)估降水偏多率要高于歷史氣候模擬(圖3a,圖2a)。同樣冷季的預(yù)估降水與歷史氣候模擬降水在我國華南沿海和南亞等地也存在明顯差異(圖3b,圖2b)。這些差異說明,模式預(yù)估降水偏差還受到非定常偏差的影響,同樣需要加以訂正。

    圖2 1960~2005 年20 個CMIP5 模式集合平均亞洲降水的偏差百分率:(a)暖季(5~10 月)降水;(b)冷季(11 月至次年4 月)降水Fig. 2 Bias percentage of Asian precipitation between the ensemble mean of 20 CMIP5 models during 1960-2005: (a) Warm season (May-October);(b) cold season (November throughout next April)

    圖3 RCP4.5 情景下2006~2015 年20 個CMIP5 模式集合平均與觀測(CRU TS v4.0)亞洲降水的偏差百分率:(a)暖季;(b)冷季Fig. 3 Bias percentage of Asian precipitation between the ensemble mean of 20 CMIP5 models and the observation (CRU TS v 4.0) during 2006-2015 under the RCP4.5 scenario: (a) Warm season; (b) cold season

    5 模式降水偏差訂正

    5.1 模式氣候漂移

    在CMIP5 模式集合降水預(yù)估中扣除模式氣候漂移有多種方式,例如,可以在模式降水預(yù)估中減去歷史氣候模擬中一個時期平均的降水偏差,或常用的做法是只計算預(yù)估降水的相對變化,即計算一段時間平均的預(yù)估降水與另一時間段平均的歷史氣候模擬降水之差,以此扣除所謂模式氣候漂移,其中假定模式氣候漂移不隨時間段變化。這在很大程度上可以削弱模式降水的氣候漂移偏差,減小降水預(yù)估的不確定性。圖4 給出了2006~2015 年CMIP5多模式集合平均降水預(yù)估相對于1976~2005 年CMIP5 歷史氣候模擬平均值(參考態(tài))的距平百分率。圖4a 顯示,在中等排放情景RCP4.5 下,2006~2015 年暖季預(yù)估降水在亞洲大陸大部分地區(qū)都有所增多,其中,我國西北中東部和南亞西北部降水增加率為5%~10%,阿拉伯半島降水增加10%~20%,在我國江南、四川盆地、華北西部、中南半島中部、印度的部分地區(qū)、青藏高原西部、蒙古國東部、里海周邊及東非東北海岸等地降水減少。在冷季,2006~2015 年模式預(yù)估降水在北方增多,南方減少,其中降水增多比較明顯的區(qū)域位于我國西北東部、印度西北部、巴基斯坦、貝加爾湖附近及遠(yuǎn)東等地,增加幅度在5%~10%之間;而降水減少的區(qū)域主要位于中國中東部、南部、西南部、青藏高原東南部、印度東北部、伊朗、及阿拉伯半島等地(圖4b)。

    圖4 RCP4.5 情景下2006~2015 年20 個CMIP5 模式集合預(yù)估降水相對于1976~2005 年歷史氣候模擬降水的距平百分率(上)及其與觀測降水(CRU TS v4.0)的距平符號一致格點(“+”)分布(下):(a、c)暖季;(b、d)冷季Fig. 4 Percentage of the precipitation anomaly projected using 20 CMIP5 models ensemble mean vs the counterpart of its historical simulations for 1976-2005 (upper panel) and distributions of grid points with the same sign (marked as “+”) in the precipitation anomaly with the observation (CRU TS v4.0,lower panel) in Asia during 2006-2015 under the RCP4.5 scenario: (a, c) Warm season; (b, d) cold season

    然而,進(jìn)一步分析扣除氣候漂移后預(yù)估降水與觀測距平符號的一致性,發(fā)現(xiàn)在許多地區(qū)都存在不一致的情形,即存在預(yù)估盲區(qū)。圖4c 顯示幾乎有一半格點的預(yù)估降水距平與觀測距平相反,即沒有預(yù)估能力,距平符號一致率僅為49%。這說明圖4a的預(yù)估降水距平在這些“盲區(qū)”(距平反號區(qū))無效,其中包括中國的東北和華北北部、中東部大部分地區(qū)、西北大部及青藏高于南部等地。但是,冷季模式預(yù)估降水的距平符號一致率較高,達(dá)到68%,距平符號不一致的格點主要集中在40°N 附近緯度帶及中亞等地。這說明,僅扣除模式氣候漂移存在一定局限性,模式降水預(yù)估的非定常偏差亦很重要,需要對其做進(jìn)一步訂正,以降低模式降水預(yù)估的不確定性。

    除距平符號一致率外,模式降水預(yù)估偏差還存在距平幅度的偏差,即盡管模式降水距平符號與觀測一致,若距平幅度與觀測相比太小也沒有實際意義。對比圖5a、b 與圖4a、b 發(fā)現(xiàn),大部分地區(qū)模式預(yù)估降水距平都遠(yuǎn)小于觀測。例如,在模式預(yù)估的暖季區(qū)域性降水增多幅度明顯較小的區(qū)域有我國西北東部或三江源區(qū)、印度次大陸、巴基斯坦和阿富汗、中南半島西部,以及高緯度地帶;模式預(yù)估降水負(fù)距平幅度明顯偏小區(qū)域為我國西南及東北部分地區(qū)、伊朗南部、里海附近及印度東北部等(圖5a,圖4a)。比較圖4b 與圖5b 可知,冷季模式預(yù)估降水距平幅度明顯低于觀測,盡管模式預(yù)估的冷季降水距平與觀測之間的符號一致率遠(yuǎn)高于暖季。例如,模式冷季降水增多幅度在我國華北和東北、印度次大陸中部及西北部、中南半島大部、菲律賓群島、俄羅斯遠(yuǎn)東及西西伯利亞等地明顯偏?。辉谀J筋A(yù)估冷季降水偏少的區(qū)域距平幅度也明顯偏弱,如我國西南部、青藏高原東部、印度西北部、伊朗高原及阿拉伯半島東部等地。因此,模式降水偏差訂正不僅要提高其距平符號一致率,也要考慮解決模式預(yù)估降水距平幅度過低的問題。此外,為進(jìn)一步減少模式偏差,除了模式氣候漂移還應(yīng)考慮對模式降水非定常偏差的訂正。

    圖5 2006~2015 年亞洲觀測(CRU TS v4.0)降水距平百分率:(a)暖季;(b)冷季。參考態(tài):1976~2005 年Fig. 5 Percentage of precipitation anomaly in Asia from CRU TS v4.0 for 2006-2015: (a) Warm season; (b) cold season. Reference years: 1976-2005

    5.2 回歸偏差訂正

    模式降水的定常和非定常偏差可以通過回歸模型等方式統(tǒng)一訂正,但統(tǒng)計方法可以訂正的僅為非定常偏差中的統(tǒng)計平穩(wěn)部分。根據(jù)楊陽等(2019)對中國大陸降水預(yù)估偏差訂正的檢驗表明,一元對數(shù)回歸和一元年際增量回歸等模型對模式降水偏差訂正的效果較好。因此,本文選擇一元對數(shù)回歸和一元年際增量回歸模型,檢驗其對亞洲模式降水預(yù)估偏差的訂正效果。這兩種統(tǒng)計模型的建立都基于1961~2005 年CMIP5 歷史氣候模擬降水?dāng)?shù)據(jù)和降水觀測數(shù)據(jù)(CRU TS v4.0)。

    圖6a 表示,使用一元對數(shù)回歸訂正偏差后的2006~2015 年暖季模式預(yù)估降水相對于1976~2005 年觀測氣候平均值的距平及其與觀測降水距平同號格點分布(“+”),距平一致格點率為56.6%,明顯高于圖4c 的49%,而且降水距平幅度與圖4a相當(dāng)。但圖6c 顯示,一元對數(shù)回歸偏差訂正在南亞次大陸幾乎無效,其在30°N 以南距平同號率僅為46.7%,以北達(dá)到60%。這說明一元對數(shù)回歸訂正法在低緯度存在較多盲區(qū),比較適合訂正中高緯度暖季模式降水偏差。

    圖6 RCP4.5 情景下2006~2015 年線性回歸訂正后亞洲暖季模式預(yù)估降水與觀測降水距平百分率(上)及其距平同號格點分布(下):(a、c)一元對數(shù)回歸;(b、d)一元年際增量回歸。參考態(tài):1976~2005 年Fig. 6 Anomaly percentage of the warm season precipitation projected using 20 CMIP5 model ensemble mean under RCP4.5 for 2006-2015 (upper)and the corresponding grid distribution with the correct sign (“+”) with the observation precipitation anomaly (lower): (a, c) Logarithmic regression; (b,d) year-to-year increment regression. Reference years: 1976-2005

    比較圖6b 與圖5b 可知,一元年際增量回歸訂正后預(yù)估降水距平幅度高于一元對數(shù)回歸,與觀測比較接近,其距平符號一致率為62%,高于對數(shù)回歸,其30°N 以南距平同號率甚至達(dá)到約67%,也明顯高于對數(shù)回歸,但在30°N 以北與對數(shù)回歸相當(dāng)(圖6d)。因此,在暖季可以用所謂“組合訂正法”,即在30°N 以南選用一元年際增量回歸進(jìn)行偏差訂正,而30°N 以北用一元對數(shù)回歸訂正法。鑒于一元對數(shù)回歸訂正后降水距平在北方太小,在暖季本文還是使用一元年際增量回歸訂正,盡管其降水距平在北方幅度超過觀測降水距平(見圖5a)。

    在冷季,圖7a,b 顯示,一元對數(shù)回歸訂正后降水距平幅度依舊小于一元年際增量回歸,后者的幅度接近于觀測值(圖5b)。而圖7c,d 中一元對數(shù)回歸的距平符號一致率僅為48%,一元年際增量回歸的略高,為51%。根據(jù)圖7c 中距平符號不一致的區(qū)域即盲區(qū)(空白處),可以推斷使用一元對數(shù)回歸偏差訂正后冷季在華北、南亞西北部及中南半島大部等地的降水偏少無效(圖7a)。同樣,使用一元年際增量回歸偏差訂正后冷季降水在我國華北、中南半島等地偏少,在我國三江源區(qū)、中亞南部等地偏多的估計也均無效。比較圖7c,d不難看出,二者在南方和北方的距平符號一致率不同,距平同號格點的分布亦不一樣。統(tǒng)計結(jié)果表明,一元對數(shù)回歸在30°N 以南為56%,遠(yuǎn)高于一元年際增量回歸的48%;相反,前者的距平同號率在北方只有45%,低于后者的53%。因此,對于2006~2015 年冷季降水預(yù)估偏差訂正,可以用所謂“組合訂正法”,即在30°N 以南選用一元對數(shù)回歸訂正,以北用一元年際增量回歸訂正。

    圖7 同圖6,但為冷季降水Fig. 7 Same as Fig. 6, but for cold season precipitation

    顯然,單純地使用一元對數(shù)回歸或一元年際增量回歸偏差訂正法都不能得到冷季預(yù)估降水的最優(yōu)偏差訂正效果,若將二者組合起來可以得到更優(yōu)結(jié)果。圖8 即是用區(qū)域組合回歸法進(jìn)行偏差訂正后得到的2006~2015 年亞洲暖季和冷季降水距平百分率分布,其中暖季在30°N 以南用一元年際增量回歸訂正,以北用一元對數(shù)回歸訂正(圖8a)。相反,冷季在在30°N 以南用一元對數(shù)訂正,以北用一元年際增量訂正(圖8b)。顯然,組合訂正后圖8 與圖5 的距平幅度大小更為接近,距平符號一致率也高于單純的一元年際增量回歸訂正或一元對數(shù)回歸訂正的結(jié)果。不過,如果參考圖6c,d 可知,圖8a 中的暖季降水預(yù)估結(jié)果在中南半島中南部、南亞西北部、西北東部等地的預(yù)估均無效,因為是盲區(qū)(圖6c,d)。但些盲區(qū)在圖4c 中卻是有效區(qū),可以用其中的預(yù)估降水距平填補(bǔ)這些區(qū)域,以降低預(yù)估的不確定性。同樣,參考圖7c,d 可知,在圖8b 中的印度中南部、南亞西北部、中南半島中北部、中國西南部、華北大部、西北東部及新疆南部等地的區(qū)域組合降水預(yù)估亦均無效。也可以參考圖4d,用其降水距平填補(bǔ)這些盲區(qū)。因此,區(qū)域組合訂正法也需要結(jié)合去模式氣候漂移的結(jié)果進(jìn)行補(bǔ)充訂正,以進(jìn)一步減少亞洲降水預(yù)估的不確定性。

    圖8 在RCP4.5 情景下20 個CMIP5 模式集合預(yù)估并進(jìn)行組合回歸訂正后2006~2015 年亞洲降水距平百分率: (a)在暖季30°N 以北用一元對數(shù)回歸訂正,以南用一元年際增量回歸訂正;(b)冷季與暖季相反。參考態(tài):1976~2005 年Fig. 8 Asian precipitation anomaly percentage projected using CMIP5 20 model ensemble mean under the RCP4.5 scenario with bias correction for 2006-2015 versus that of the 1976-2005 observation: (a) Logarithm regression in the north of 30°N with year-to-year increment regression in the south for the warm seasons; (b) opposite combination in the regressions for the cold seasons

    6 未來30 年亞洲降水預(yù)估

    在RCP4.5 情景下用區(qū)域組合回歸訂正法對20 個CMIP5 模式集合平均預(yù)估的未來30 年(2021~2050)亞洲降水進(jìn)行偏差訂正,其中暖季用一元年際增量回歸訂正,冷季在30°N 以南用一元對數(shù)回歸訂正,30°N 以北用一元年際增量回歸訂正。對于一元對數(shù)或年際增量回歸訂正法可能的盲區(qū)用單純?nèi)ツJ狡平Y(jié)果填補(bǔ),如果后者在檢驗期屬于有效區(qū),否則不予填補(bǔ)訂正。圖9a 顯示,相對于1976~2005 年氣候平均,2021~2050 年暖季降水距平分布與2006~2015 年類似(見圖6b),即中國西南部和江南中部、南亞東北部、環(huán)里海南部、阿拉伯半島東北部等地降水偏少10%~20%;從中國的三江源區(qū)到淮河流域降水增加20%左右,東北南部降水增加約10%;新疆北部降水增加約10%,南部增加約20%,華北和東北大部降水減少10%~20%;中南半島北部降水增加約10%,南部變化不大(圖9c);亞洲高緯度地帶降水略有增加。

    圖9b,d 顯示,在2021~2050 年亞洲冷季預(yù)估降水呈現(xiàn)北方增加、南方減少的格局。其中南亞降水減少最明顯,達(dá)-10%左右,中國西南部約為~5%;西北和青藏高原降水增加幅度約為20%~40%,華北和東北增加約5%,西伯利亞及遠(yuǎn)東降水增加約為10%~40%。這里考慮到一元對數(shù)回歸與年際增量回歸在冷季降水預(yù)估的盲區(qū)較多,例如在南亞、中南半島、中國北方大部等(圖7c,d),因此在降水距平符號選擇上主要參考圖9d 及圖7d;在降水距平幅度估計上,如果在圖9b 與圖9d 中同號,則選擇前者,否則選擇后者,若遇到二者共同的盲區(qū)則不再做預(yù)估推薦。

    圖9 RCP4.5 情景下經(jīng)偏差訂正后20 個CMIP5 模式集合預(yù)估的2021~2050 年亞洲年暖季(左)和冷季(右)降水距平百分率:(a)一元年際增量回歸訂正;(b)組合回歸訂正即一元年際增量回歸(30°N 以北)和一元對數(shù)回歸(30°N 以南);(c,d)去除模式氣候漂移后暖季和冷季降水距平百分率。圖a 和b 的參考態(tài):1976~2005 年觀測降水;圖c 和d 的參考態(tài):1976~2005 年 CMIP5 歷史模擬降水Fig. 9 Bias-corrected precipitation anomaly percentage projected using CMIP5 20 model ensemble mean under RCP4.5 scenario for 2021-2050 in the warm season (left) and cold season (right): (a) Year-to-year increment regression; (b) year-to-year increment regression in the north of 30°N and logarithm in the south; (c, d) model drift removed for warm and cold seasons, respectively. Figures a and b are with reference to the 1976-2005 observation; figures c and d are with reference to the precipitation of the CMIP5 historical climate simulation for 1976-2005

    7 討論與結(jié)論

    本文分析了CMIP5 多模式集合歷史模擬及預(yù)估降水偏差在亞洲的分布特征,檢驗了三種偏差訂正法,即去模式氣候漂移、一元對數(shù)回歸及年際增量回歸在2006~2015 年亞洲降水預(yù)估中的偏差訂正效果,構(gòu)造了區(qū)域組合回歸偏差訂正法,預(yù)估了2021~2050 年暖季和冷季亞洲降水的可能變化。主要結(jié)論如下:

    (1)在CMIP5 歷史氣候模擬中,多模式集合平均降水在亞洲存在明顯偏差,在中國華南沿海、南亞、越南、朝鮮半島、伊朗南部及阿拉伯半島等地均偏少,其中在南亞明顯偏少,達(dá)-20%~-40%,而在青藏高原、內(nèi)蒙古、蒙古國等地偏多達(dá)30%~40%。

    (2)2006~2015 年CMIP5 多模式集合預(yù)估降水偏差分布與歷史氣候模擬類似,說明模式的氣候漂移偏差具有準(zhǔn)定常性,可以通過歷史氣候模擬確定并在預(yù)估中扣除。

    (3)用CMIP5 預(yù)估時段的平均降水模式減去歷史氣候模擬時段平均降水可以在很大程度上消去模式氣候漂移,但預(yù)估的降水距平幅度明顯偏弱,盡管距平符號一致率較高。

    (4)對于2006~2015 年降水預(yù)估檢驗了兩種偏差訂正方法,即一元對數(shù)和年際增量回歸訂正法,發(fā)現(xiàn)其各有優(yōu)勢區(qū)域及盲區(qū),偏差訂正后前者在北方的距平幅度較小,后者偏大。以此構(gòu)造了南北區(qū)域組合回歸訂正法,即暖季在30°N 以南使用一元年際增量回歸,以北用一元對數(shù)或年際增量回歸訂正;相反,冷季在30°N 以南用一元對數(shù)回歸訂正,以北用一元年際增量回歸訂正,可以獲得更優(yōu)的偏差訂正效果。

    (5)對于組合回歸訂正法的盲區(qū),用去模式氣候漂移后的結(jié)果填補(bǔ),若也遇到盲區(qū)則放棄預(yù)估推薦。

    (6)經(jīng)區(qū)域組合回歸偏差訂正及去模式氣候漂移的附加訂正后,在RCP4.5 情景下對2021~2050 年的亞洲降水進(jìn)行了預(yù)估。結(jié)果表明,相對于1976~2005 年氣候平均,在暖季,中國南方、南亞東北部、中亞南部、阿拉伯半島東北部等地降水減少10%~20%,從中國的三江源區(qū)到淮河流域帶降水增加約20%,東北南部降水增加約10%,新疆北部降水增加約10%,南部約20%;華北和東北大部降水減少約10%~20%;中南半島北部降水增加約10%,南部變化不大;亞洲高緯度地帶降水略有增加。在冷季,亞洲降水呈現(xiàn)北方增加,南方減少的格局,其中南亞降水減少最明顯,達(dá)-10%左右,中國西南部約為-5%,中國西部降水增加幅度約為20%~40%,華北和東北增加約5%,亞洲高緯度降水增加約為10%~40%。

    降水預(yù)估結(jié)果表明,隨著全球氣候暖化未來30 年的亞洲的降水格局會發(fā)生或多或少的變化,中國西部洪澇可能會增多,暖季西南和華北降水可能會減少,冷季南亞北部和中亞南部等地可能出現(xiàn)干旱災(zāi)害。不同偏差訂正結(jié)果的檢驗還表明,可以設(shè)計更加復(fù)雜的偏差訂正方法,更細(xì)分辨率的偏差訂正組合方案,以期進(jìn)一步減小模式降水預(yù)估的不確定性。

    猜你喜歡
    距平對數(shù)增量
    颶風(fēng)Edouard(2014)暖心結(jié)構(gòu)的多資料對比分析
    含有對數(shù)非線性項Kirchhoff方程多解的存在性
    提質(zhì)和增量之間的“辯證”
    指數(shù)與對數(shù)
    指數(shù)與對數(shù)
    “價增量減”型應(yīng)用題點撥
    對數(shù)簡史
    近40年阿里地區(qū)云量和氣溫的年際變化
    西藏科技(2018年9期)2018-10-17 05:51:30
    甘肅省降水和冰雹天氣氣候分析
    基于均衡增量近鄰查詢的位置隱私保護(hù)方法
    最近最新中文字幕大全电影3 | 一进一出抽搐动态| 亚洲国产看品久久| a级毛片黄视频| 欧美成人免费av一区二区三区 | av网站免费在线观看视频| 亚洲欧美色中文字幕在线| 亚洲精品一卡2卡三卡4卡5卡| 很黄的视频免费| 51午夜福利影视在线观看| 久久久水蜜桃国产精品网| 丁香欧美五月| 久久久精品区二区三区| 热re99久久国产66热| 99热网站在线观看| av网站免费在线观看视频| 国产免费现黄频在线看| 国产一卡二卡三卡精品| 黄色 视频免费看| 精品久久久精品久久久| 欧美另类亚洲清纯唯美| 一边摸一边抽搐一进一小说 | 91国产中文字幕| 99精品欧美一区二区三区四区| 久久国产精品男人的天堂亚洲| 少妇 在线观看| 久久精品亚洲精品国产色婷小说| 9色porny在线观看| 欧美日韩成人在线一区二区| 中文字幕人妻丝袜制服| 亚洲av日韩精品久久久久久密| 亚洲全国av大片| 一边摸一边抽搐一进一出视频| 黄色怎么调成土黄色| 久久午夜综合久久蜜桃| 操美女的视频在线观看| 如日韩欧美国产精品一区二区三区| 欧美日韩乱码在线| 国产精品免费一区二区三区在线 | 亚洲国产毛片av蜜桃av| 欧美日韩精品网址| 久久久久久久久免费视频了| 免费久久久久久久精品成人欧美视频| 极品人妻少妇av视频| 18在线观看网站| 69av精品久久久久久| 视频区欧美日本亚洲| 日韩欧美一区视频在线观看| 搡老岳熟女国产| 侵犯人妻中文字幕一二三四区| 午夜激情av网站| 欧美激情久久久久久爽电影 | 免费少妇av软件| 涩涩av久久男人的天堂| 久久国产精品影院| 免费黄频网站在线观看国产| 精品一区二区三区av网在线观看| 国产亚洲欧美98| 性少妇av在线| 精品熟女少妇八av免费久了| 国产成人影院久久av| 51午夜福利影视在线观看| 亚洲七黄色美女视频| 国产一区二区三区在线臀色熟女 | 午夜影院日韩av| 天天躁夜夜躁狠狠躁躁| 少妇粗大呻吟视频| 日本黄色视频三级网站网址 | 久久久久久久久久久久大奶| 51午夜福利影视在线观看| 精品国产一区二区三区久久久樱花| 久久香蕉国产精品| 乱人伦中国视频| 亚洲人成伊人成综合网2020| 一级a爱视频在线免费观看| 久久午夜综合久久蜜桃| 国产成人精品无人区| 丝袜美腿诱惑在线| 国产av精品麻豆| 叶爱在线成人免费视频播放| 深夜精品福利| 久久天堂一区二区三区四区| 99re在线观看精品视频| 色在线成人网| 嫩草影视91久久| 日韩制服丝袜自拍偷拍| 悠悠久久av| 手机成人av网站| 欧美乱色亚洲激情| 国产精品电影一区二区三区 | 亚洲精品国产区一区二| 老司机在亚洲福利影院| 老司机深夜福利视频在线观看| 中文字幕高清在线视频| 亚洲在线自拍视频| 国产精品免费大片| 国产亚洲欧美精品永久| 久久草成人影院| 亚洲avbb在线观看| 亚洲成国产人片在线观看| 亚洲成国产人片在线观看| bbb黄色大片| 男人操女人黄网站| 久久午夜亚洲精品久久| 久久午夜亚洲精品久久| 大香蕉久久成人网| 欧美日韩视频精品一区| 久久精品亚洲av国产电影网| 中国美女看黄片| 中亚洲国语对白在线视频| 51午夜福利影视在线观看| 国产伦人伦偷精品视频| 99国产综合亚洲精品| 国产精品av久久久久免费| 日本五十路高清| 久久性视频一级片| 国产熟女午夜一区二区三区| 免费不卡黄色视频| 婷婷丁香在线五月| 国产aⅴ精品一区二区三区波| 国内久久婷婷六月综合欲色啪| 国产精品久久电影中文字幕 | 久久 成人 亚洲| 在线观看免费高清a一片| 久久精品国产综合久久久| 午夜91福利影院| 国产精品免费视频内射| 黄色成人免费大全| 亚洲av片天天在线观看| 亚洲av第一区精品v没综合| 亚洲精品av麻豆狂野| 人人澡人人妻人| 脱女人内裤的视频| 国产成人精品久久二区二区免费| 丰满迷人的少妇在线观看| 中文字幕色久视频| 女人精品久久久久毛片| 亚洲精品乱久久久久久| 丰满人妻熟妇乱又伦精品不卡| 国内久久婷婷六月综合欲色啪| xxxhd国产人妻xxx| 91精品国产国语对白视频| 免费在线观看视频国产中文字幕亚洲| 国产精品综合久久久久久久免费 | ponron亚洲| 少妇粗大呻吟视频| 色在线成人网| 国产高清激情床上av| 在线永久观看黄色视频| 亚洲av成人av| 免费高清在线观看日韩| 亚洲九九香蕉| 麻豆av在线久日| 成人黄色视频免费在线看| 美女 人体艺术 gogo| 高清av免费在线| 黑人欧美特级aaaaaa片| 看片在线看免费视频| 91九色精品人成在线观看| 手机成人av网站| 麻豆乱淫一区二区| 91字幕亚洲| 精品一区二区三卡| 欧美黄色片欧美黄色片| 国产三级黄色录像| 免费人成视频x8x8入口观看| 久久久国产成人免费| 精品久久久久久电影网| 欧美中文综合在线视频| 黄色怎么调成土黄色| 欧美乱码精品一区二区三区| 中文字幕精品免费在线观看视频| 欧美不卡视频在线免费观看 | 亚洲av欧美aⅴ国产| 日韩大码丰满熟妇| 精品一区二区三区四区五区乱码| 久久久久视频综合| 欧洲精品卡2卡3卡4卡5卡区| 91精品国产国语对白视频| 激情在线观看视频在线高清 | 中亚洲国语对白在线视频| 国产精品美女特级片免费视频播放器 | 51午夜福利影视在线观看| 亚洲 国产 在线| 日韩成人在线观看一区二区三区| 免费少妇av软件| a级片在线免费高清观看视频| 亚洲熟妇中文字幕五十中出 | 久久精品国产清高在天天线| 欧美激情极品国产一区二区三区| 亚洲黑人精品在线| 美女 人体艺术 gogo| 久久热在线av| 黑人操中国人逼视频| 亚洲中文字幕日韩| 巨乳人妻的诱惑在线观看| 亚洲精品国产精品久久久不卡| 夜夜爽天天搞| 欧美黄色淫秽网站| 最新美女视频免费是黄的| 国产亚洲精品一区二区www | 亚洲黑人精品在线| 久久香蕉激情| 狂野欧美激情性xxxx| 亚洲色图综合在线观看| 老熟妇乱子伦视频在线观看| 精品国产超薄肉色丝袜足j| 一区二区日韩欧美中文字幕| 亚洲性夜色夜夜综合| 99精品久久久久人妻精品| 午夜成年电影在线免费观看| 在线观看免费高清a一片| 久热爱精品视频在线9| 村上凉子中文字幕在线| 黄色视频不卡| 亚洲精品自拍成人| 美女 人体艺术 gogo| 欧美精品av麻豆av| 欧美乱色亚洲激情| 国产高清视频在线播放一区| 亚洲精品美女久久久久99蜜臀| 欧美乱妇无乱码| 淫妇啪啪啪对白视频| 日韩中文字幕欧美一区二区| av中文乱码字幕在线| 日韩欧美一区二区三区在线观看 | 欧美日韩福利视频一区二区| 大陆偷拍与自拍| 99国产精品一区二区蜜桃av | 欧美日韩黄片免| 在线永久观看黄色视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲片人在线观看| 国产高清视频在线播放一区| 国产精品久久久久成人av| 国产精品1区2区在线观看. | 老司机深夜福利视频在线观看| 校园春色视频在线观看| 两个人看的免费小视频| 国产成人免费观看mmmm| 亚洲综合色网址| 亚洲一区高清亚洲精品| 无人区码免费观看不卡| 免费观看a级毛片全部| 少妇的丰满在线观看| 久久久久久亚洲精品国产蜜桃av| 伦理电影免费视频| 亚洲自偷自拍图片 自拍| 成人特级黄色片久久久久久久| 涩涩av久久男人的天堂| 他把我摸到了高潮在线观看| 视频区欧美日本亚洲| 视频区欧美日本亚洲| 国产av精品麻豆| 欧美黄色片欧美黄色片| 成人免费观看视频高清| 国产精品久久视频播放| 免费看十八禁软件| 国产精品久久视频播放| 香蕉国产在线看| 18在线观看网站| 亚洲精品自拍成人| 夜夜爽天天搞| 热99re8久久精品国产| 国产欧美日韩精品亚洲av| 国产精品二区激情视频| 免费少妇av软件| 黄色女人牲交| 岛国毛片在线播放| 电影成人av| 国产av又大| 久久久久视频综合| 99精品在免费线老司机午夜| 又黄又爽又免费观看的视频| 亚洲全国av大片| 国产主播在线观看一区二区| 亚洲人成77777在线视频| 中亚洲国语对白在线视频| 黄色片一级片一级黄色片| 最近最新中文字幕大全电影3 | a级毛片在线看网站| 操美女的视频在线观看| 中文字幕最新亚洲高清| 欧美日韩瑟瑟在线播放| 女人被躁到高潮嗷嗷叫费观| 国产高清激情床上av| 国产伦人伦偷精品视频| 国产精品偷伦视频观看了| 熟女少妇亚洲综合色aaa.| 精品卡一卡二卡四卡免费| 久久精品亚洲精品国产色婷小说| 亚洲 国产 在线| 在线免费观看的www视频| 亚洲五月天丁香| 亚洲第一欧美日韩一区二区三区| 精品国产国语对白av| 中文字幕制服av| 国产亚洲欧美在线一区二区| 99精品欧美一区二区三区四区| 日韩大码丰满熟妇| 精品国产乱码久久久久久男人| 精品一区二区三区四区五区乱码| 日韩欧美在线二视频 | 国产精品香港三级国产av潘金莲| 欧美激情极品国产一区二区三区| 亚洲一区二区三区欧美精品| 国产精品免费大片| 精品国内亚洲2022精品成人 | 婷婷精品国产亚洲av在线 | 久久热在线av| 久久国产精品影院| 午夜福利一区二区在线看| 久久99一区二区三区| 亚洲成a人片在线一区二区| 怎么达到女性高潮| 亚洲精品中文字幕在线视频| 超碰97精品在线观看| 韩国精品一区二区三区| 国产主播在线观看一区二区| 一二三四社区在线视频社区8| 国内毛片毛片毛片毛片毛片| 久久性视频一级片| 欧美黄色淫秽网站| 脱女人内裤的视频| 久久久久久免费高清国产稀缺| 午夜福利一区二区在线看| av网站免费在线观看视频| 国产成人欧美在线观看 | 欧美黄色片欧美黄色片| 大香蕉久久成人网| 婷婷丁香在线五月| 久久精品亚洲熟妇少妇任你| 天天影视国产精品| 亚洲国产欧美一区二区综合| 国产成人啪精品午夜网站| 国产精品欧美亚洲77777| 老司机深夜福利视频在线观看| 伊人久久大香线蕉亚洲五| 他把我摸到了高潮在线观看| 亚洲av日韩在线播放| 欧美黄色淫秽网站| 午夜福利影视在线免费观看| 女人被躁到高潮嗷嗷叫费观| 日本黄色日本黄色录像| 男人的好看免费观看在线视频 | 一二三四社区在线视频社区8| 在线观看免费日韩欧美大片| 精品第一国产精品| 中文字幕制服av| 五月开心婷婷网| 99热国产这里只有精品6| 人人妻人人添人人爽欧美一区卜| 热re99久久精品国产66热6| 国产精品二区激情视频| 欧美成狂野欧美在线观看| 亚洲情色 制服丝袜| 久久久久久免费高清国产稀缺| 18禁国产床啪视频网站| 交换朋友夫妻互换小说| 可以免费在线观看a视频的电影网站| 色播在线永久视频| 老熟妇仑乱视频hdxx| 欧美黑人欧美精品刺激| 一级片免费观看大全| 亚洲avbb在线观看| 成年版毛片免费区| 精品视频人人做人人爽| 国产伦人伦偷精品视频| 美女扒开内裤让男人捅视频| 久久久久国产精品人妻aⅴ院 | 国产成人影院久久av| 最新美女视频免费是黄的| 一夜夜www| 久久久国产精品麻豆| 亚洲五月婷婷丁香| 国产日韩一区二区三区精品不卡| 久久婷婷成人综合色麻豆| 一个人免费在线观看的高清视频| av片东京热男人的天堂| 两性午夜刺激爽爽歪歪视频在线观看 | 岛国在线观看网站| 欧美黑人精品巨大| 亚洲国产精品一区二区三区在线| 成熟少妇高潮喷水视频| 水蜜桃什么品种好| a级片在线免费高清观看视频| 丝袜人妻中文字幕| 两个人看的免费小视频| 亚洲全国av大片| 日韩欧美免费精品| 法律面前人人平等表现在哪些方面| 亚洲欧美日韩另类电影网站| 亚洲精品美女久久av网站| 精品国产一区二区久久| 国产免费男女视频| 日韩精品免费视频一区二区三区| 视频在线观看一区二区三区| 久久人人爽av亚洲精品天堂| 日本黄色视频三级网站网址 | 成年女人毛片免费观看观看9 | 国产欧美日韩一区二区三区在线| 欧美午夜高清在线| 后天国语完整版免费观看| 色尼玛亚洲综合影院| 国产精品99久久99久久久不卡| 国产精品二区激情视频| 无人区码免费观看不卡| 国产免费男女视频| 91麻豆精品激情在线观看国产 | a在线观看视频网站| 亚洲色图综合在线观看| 在线视频色国产色| 一级,二级,三级黄色视频| 我的亚洲天堂| 国产xxxxx性猛交| 午夜精品在线福利| 亚洲中文日韩欧美视频| 久久精品国产a三级三级三级| 91av网站免费观看| 亚洲成国产人片在线观看| 香蕉久久夜色| 亚洲精品一二三| 一级作爱视频免费观看| 女人爽到高潮嗷嗷叫在线视频| 精品人妻在线不人妻| 国产黄色免费在线视频| 国产日韩一区二区三区精品不卡| 日韩免费av在线播放| 精品免费久久久久久久清纯 | 日韩欧美国产一区二区入口| 国产成人精品久久二区二区91| 9色porny在线观看| 精品人妻1区二区| 夜夜爽天天搞| 欧美最黄视频在线播放免费 | 91老司机精品| 欧美av亚洲av综合av国产av| 亚洲五月婷婷丁香| 色播在线永久视频| 亚洲国产看品久久| 91国产中文字幕| 美女福利国产在线| 国产又爽黄色视频| 青草久久国产| 咕卡用的链子| 少妇被粗大的猛进出69影院| 国产成人一区二区三区免费视频网站| 亚洲专区字幕在线| 宅男免费午夜| 亚洲成a人片在线一区二区| 咕卡用的链子| 午夜福利一区二区在线看| 99热国产这里只有精品6| 不卡av一区二区三区| 国产精品九九99| 两人在一起打扑克的视频| 亚洲欧美日韩另类电影网站| 亚洲精品美女久久av网站| 国产亚洲欧美精品永久| 午夜激情av网站| 欧美日韩中文字幕国产精品一区二区三区 | 99riav亚洲国产免费| 咕卡用的链子| 国产一区二区三区综合在线观看| 精品久久蜜臀av无| 国产精品98久久久久久宅男小说| 国产高清激情床上av| 操出白浆在线播放| av不卡在线播放| 天堂动漫精品| 欧美成人免费av一区二区三区 | 国产xxxxx性猛交| 国产精华一区二区三区| 国产精品久久久久成人av| 日本精品一区二区三区蜜桃| 国产成人欧美在线观看 | 母亲3免费完整高清在线观看| 亚洲在线自拍视频| 亚洲欧美色中文字幕在线| 久久久国产成人免费| 欧美日韩精品网址| 90打野战视频偷拍视频| svipshipincom国产片| 日韩欧美国产一区二区入口| 俄罗斯特黄特色一大片| 国产激情久久老熟女| 亚洲av成人不卡在线观看播放网| 成在线人永久免费视频| 亚洲综合色网址| 超碰97精品在线观看| 成在线人永久免费视频| 久久久国产欧美日韩av| 欧美日韩亚洲综合一区二区三区_| 国产男女内射视频| 男女午夜视频在线观看| 老鸭窝网址在线观看| av电影中文网址| 少妇被粗大的猛进出69影院| 欧美黄色淫秽网站| 好男人电影高清在线观看| 麻豆av在线久日| 在线国产一区二区在线| 久久久久精品人妻al黑| 法律面前人人平等表现在哪些方面| 日本五十路高清| 日韩有码中文字幕| 久久国产精品影院| 国产精品一区二区在线观看99| 在线av久久热| 91字幕亚洲| 久久午夜亚洲精品久久| 在线观看www视频免费| 欧美精品一区二区免费开放| 水蜜桃什么品种好| 天堂中文最新版在线下载| 亚洲 欧美一区二区三区| 热99久久久久精品小说推荐| 中文字幕色久视频| 每晚都被弄得嗷嗷叫到高潮| av超薄肉色丝袜交足视频| 国产成+人综合+亚洲专区| 国产一区二区三区综合在线观看| av一本久久久久| 午夜福利欧美成人| 国产区一区二久久| 三上悠亚av全集在线观看| 91国产中文字幕| 亚洲专区国产一区二区| 久久人妻福利社区极品人妻图片| 50天的宝宝边吃奶边哭怎么回事| av在线播放免费不卡| 操美女的视频在线观看| 欧美成人午夜精品| 亚洲国产精品一区二区三区在线| 国产在视频线精品| videosex国产| 国精品久久久久久国模美| 深夜精品福利| 国产99久久九九免费精品| 高清视频免费观看一区二区| 99热网站在线观看| 欧美日韩亚洲高清精品| 成年人黄色毛片网站| 国产精品99久久99久久久不卡| 欧美激情高清一区二区三区| 国产1区2区3区精品| 一二三四社区在线视频社区8| 国产高清国产精品国产三级| 国产av一区二区精品久久| 久久久久久亚洲精品国产蜜桃av| 成人av一区二区三区在线看| 又大又爽又粗| 欧美日韩av久久| 老鸭窝网址在线观看| 99国产精品一区二区三区| 久久人妻福利社区极品人妻图片| 亚洲欧美一区二区三区久久| 日韩视频一区二区在线观看| 亚洲精品久久午夜乱码| 99热只有精品国产| 搡老乐熟女国产| 久久精品国产清高在天天线| 美女扒开内裤让男人捅视频| 波多野结衣av一区二区av| 两人在一起打扑克的视频| 老熟女久久久| 三级毛片av免费| 纯流量卡能插随身wifi吗| 一区在线观看完整版| 国产无遮挡羞羞视频在线观看| 这个男人来自地球电影免费观看| 女同久久另类99精品国产91| 国产成人精品无人区| 两性午夜刺激爽爽歪歪视频在线观看 | 啦啦啦 在线观看视频| 50天的宝宝边吃奶边哭怎么回事| 国产一卡二卡三卡精品| 欧美日韩亚洲国产一区二区在线观看 | 九色亚洲精品在线播放| 大型av网站在线播放| 国产免费男女视频| 欧美日韩乱码在线| av不卡在线播放| 18禁美女被吸乳视频| 欧美性长视频在线观看| 日韩欧美免费精品| 亚洲一区二区三区不卡视频| 99国产综合亚洲精品| 丝袜美足系列| 亚洲精品粉嫩美女一区| 亚洲在线自拍视频| 亚洲精品一二三| www.熟女人妻精品国产| 另类亚洲欧美激情| 午夜成年电影在线免费观看| 欧美黑人精品巨大| 成人18禁在线播放| 一进一出好大好爽视频| 9色porny在线观看| 99re6热这里在线精品视频| 欧美久久黑人一区二区| 亚洲成人国产一区在线观看| 欧美久久黑人一区二区| 91精品三级在线观看| 日日摸夜夜添夜夜添小说| 国产免费现黄频在线看| 好看av亚洲va欧美ⅴa在| 亚洲av第一区精品v没综合| 色老头精品视频在线观看| 色精品久久人妻99蜜桃| 一边摸一边抽搐一进一出视频| 人成视频在线观看免费观看| 久热爱精品视频在线9| 中文欧美无线码| 嫩草影视91久久|