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

    21世紀(jì)海上絲綢之路海表溫度異常與氣候變率的相關(guān)性初探

    2017-06-06 14:09:21齊慶華蔡榕碩
    海洋開發(fā)與管理 2017年4期
    關(guān)鍵詞:變率海表海區(qū)

    齊慶華,蔡榕碩

    (1.國(guó)家海洋局第三海洋研究所 廈門 361005;2.國(guó)家海洋局海洋-大氣化學(xué)與全球變化重點(diǎn)實(shí)驗(yàn)室 廈門 361005)

    ?

    21世紀(jì)海上絲綢之路海表溫度異常與氣候變率的相關(guān)性初探

    齊慶華1,2,蔡榕碩1,2

    (1.國(guó)家海洋局第三海洋研究所 廈門 361005;2.國(guó)家海洋局海洋-大氣化學(xué)與全球變化重點(diǎn)實(shí)驗(yàn)室 廈門 361005)

    區(qū)域性海洋系統(tǒng)及其環(huán)境要素(尤其是海表溫度等)受到局地和其他不同時(shí)空尺度氣候變率的交互影響,加之全球變暖背景下,區(qū)域氣候和海洋變化、變率及災(zāi)害風(fēng)險(xiǎn)的不確定性會(huì)進(jìn)一步加強(qiáng),因而分析其間的凈相關(guān)性是研究海洋環(huán)境要素變化機(jī)制機(jī)理和預(yù)測(cè)預(yù)報(bào)的重要前提。文章選取主要的海洋、大氣變率和遙相關(guān)型,在一定程度濾除干擾信號(hào)后,分析探討各氣候變率對(duì)21世紀(jì)海上絲綢之路沿途海區(qū)(簡(jiǎn)稱絲路海區(qū))海表溫度異常變化的相關(guān)性特征和差異。通過(guò)分析可為我國(guó)有關(guān)氣候變化下區(qū)域海洋環(huán)境響應(yīng)、預(yù)報(bào)減災(zāi)研究和決策以及海洋環(huán)境安全保障提供新的依據(jù)和科技支撐。

    海表溫度;氣候變率;溫度梯度;海洋環(huán)境安全;21世紀(jì)海上絲綢之路

    圖1 1950—2015年間絲路海區(qū)海表溫度線性改變量 數(shù)據(jù)來(lái)源:HadISST.

    海洋作為氣候系統(tǒng)的主要組成部分,是氣候變化的重要調(diào)節(jié)器,也是人類生存和可持續(xù)發(fā)展的戰(zhàn)略資源基地。我國(guó)地處東亞大陸東南隅,氣候受來(lái)自印度洋,我國(guó)渤海、黃海、東海和南海以及西太平洋等周邊海域熱力狀況的巨大影響。該海域是我國(guó)21世紀(jì)海上絲綢之路戰(zhàn)略實(shí)施的主戰(zhàn)場(chǎng),即“絲路海區(qū)”(圖1)。絲路海區(qū)不僅是影響我國(guó)氣候異常的關(guān)鍵區(qū)[1],也是深入理解和認(rèn)識(shí)海洋環(huán)境預(yù)報(bào)減災(zāi)以及氣候變化區(qū)域響應(yīng)和應(yīng)對(duì)研究的重點(diǎn)海域,主要任務(wù)涉及海洋與氣候變化研究與預(yù)測(cè)、海洋災(zāi)害風(fēng)險(xiǎn)分析評(píng)估以及減災(zāi)防災(zāi)決策服務(wù)和海洋環(huán)境安全保障體系建立等[2]。上述科學(xué)與實(shí)際問題的解決是保障21世紀(jì)海上絲綢之路建設(shè)實(shí)施和發(fā)展的必然要求,迫切需要開展相關(guān)的基礎(chǔ)應(yīng)用研究。絲路海區(qū)的海洋-大氣-陸地相互作用強(qiáng)烈而復(fù)雜,這里有世界上范圍最廣、海表溫度最高的西太平洋暖池區(qū)[3-5],太平洋強(qiáng)西邊界流系將低緯度能量輸送到中高緯度[6-7],亞澳3大季風(fēng)系統(tǒng)也在此交匯[8-9],這些海洋大氣系統(tǒng)不僅交互影響,而且還受其他氣候變率遙相關(guān)過(guò)程的調(diào)配,進(jìn)而對(duì)周邊地區(qū)氣候及異常變化起著關(guān)鍵性作用。其中,印度洋偶極子(Indian Ocean Dipole,IOD)、ENSO(El Nio-South Oscillation)和太平洋10年期振蕩(PDO)是印度洋和太平洋(印-太)海域主要的海洋變率[10-12],更是我國(guó)海洋預(yù)報(bào)減災(zāi)研究的重點(diǎn)。此外,Nitta[13]發(fā)現(xiàn)存在于菲律賓附近的對(duì)流活動(dòng)激發(fā)的東亞-太平洋相關(guān)型(EAP型),其對(duì)太平洋和我國(guó)近海海洋環(huán)境均有明顯影響[14-15]。而如北大西洋濤動(dòng)(NAO)和北極濤動(dòng)(AO)等大氣遙相關(guān)型與印-太海洋變率間可能存在相互作用關(guān)系[16],進(jìn)而對(duì)絲路海區(qū)海洋環(huán)境的變化產(chǎn)生聯(lián)合作用。由于氣候變化背景下局地氣候變率和遙相關(guān)型的交互影響,絲路海區(qū)海洋環(huán)境變化具有一定的復(fù)雜性,尤其是影響海表溫度變化的熱動(dòng)力過(guò)程更為復(fù)雜,而認(rèn)清各種氣候變率與絲路海區(qū)海洋環(huán)境,尤其是熱力環(huán)境變化的凈相關(guān)性,是研究海洋環(huán)境變化及其成因機(jī)理(包括強(qiáng)迫因子、影響途徑和作用貢獻(xiàn)等)乃至最終達(dá)到預(yù)報(bào)減災(zāi)目的的根本前提。目前關(guān)于絲路海區(qū)海洋環(huán)境要素與氣候變率的相關(guān)性和作用關(guān)系的研究已較多,本研究不再贅述。但有關(guān)研究大多缺乏對(duì)各時(shí)間尺度氣候變率之間相關(guān)關(guān)系和貢獻(xiàn)的辨識(shí),由于各變率間可能存在一定的關(guān)聯(lián)性,如不從一定程度剔除其他變率信號(hào)的干擾,所得的關(guān)聯(lián)性可能更多是一種綜合疊加效果,而非純粹相關(guān)關(guān)系(如偏相關(guān)分析所指的凈相關(guān)關(guān)系),這給歸因分析帶來(lái)一定影響,即無(wú)法判斷關(guān)聯(lián)性來(lái)自于哪種變率或各自的作用貢獻(xiàn)如何。從研究對(duì)象和尺度配置而言,作為背景場(chǎng)的氣候變率可以不考慮其與其他變率信號(hào)的關(guān)系,而以某個(gè)變率為主體的診斷研究則需剔除其他變率信號(hào)的影響,這對(duì)于厘清海洋系統(tǒng)內(nèi)相互作用以及與外界強(qiáng)迫的關(guān)系至關(guān)重要,也是當(dāng)前提升氣候模式預(yù)測(cè)預(yù)報(bào)水平的難點(diǎn)之一[17]。眾所周知,海表溫度是氣候和海洋環(huán)境變化預(yù)測(cè)預(yù)報(bào)和減災(zāi)防災(zāi)工作中的關(guān)鍵指標(biāo),作為初步嘗試,本研究擬側(cè)重相關(guān)性分析的技術(shù)方法及應(yīng)用,從季風(fēng)、海洋變率和大氣遙相關(guān)型入手,基于氣候變化背景下氣候系統(tǒng)內(nèi)部變率的交互影響,在一定程度(有限但必要)濾除其他干擾信號(hào)后,重新審視和探討其與絲路海區(qū)海表溫度年際和年代際異常變化的相關(guān)性特征等,以期為我國(guó)有關(guān)海洋環(huán)境變化的診斷預(yù)報(bào)和減災(zāi)防災(zāi)提供必要的科學(xué)依據(jù)。

    1 資料和方法

    1950年1月至2015年12月的月平均地表溫度主要取自英國(guó)氣象局Hadley中心的海表溫度再分析數(shù)據(jù)和英國(guó)環(huán)境數(shù)據(jù)分析中心的陸表氣溫再分析數(shù)據(jù)?;谟绊懡z路海區(qū)氣候和海洋環(huán)境預(yù)測(cè)預(yù)報(bào)的主要?dú)夂蛞蜃樱x取各變率指標(biāo)和指數(shù)序列(表1),包括主要的年際、年代際海洋變率,大氣遙相關(guān)型和季風(fēng)變化等,其中ENSO指標(biāo)以及AO、NAO和PNA指數(shù)均取自NOAA氣候預(yù)測(cè)中心網(wǎng)站。需要指出的是,目前針對(duì)各氣候變率的指標(biāo)指數(shù)較多,文中氣候變率的選擇并非全面,其指標(biāo)指數(shù)的選擇也沒有嚴(yán)格標(biāo)準(zhǔn),主要從數(shù)據(jù)源(包括時(shí)間序列的連續(xù)性等)、權(quán)威性和普適程度等方面考慮;同時(shí),由于主要關(guān)注氣候變率之間的交互作用,包括聯(lián)合協(xié)同作用以及抵消作用對(duì)環(huán)境要素變化相關(guān)性的影響,從技術(shù)方法角度考慮,指標(biāo)指數(shù)的選擇范圍和條件并不影響本研究分析結(jié)論。本研究主要利用統(tǒng)計(jì)相關(guān)分析方法探討主要海洋、大氣變率與絲路海區(qū)海表溫度的年際和年代際尺度變化上的線性相關(guān)性,同時(shí)采用回歸分析方法濾除其他變率信號(hào)。

    表1 1950—2015年各類氣候變率指標(biāo)、指數(shù)相關(guān)系數(shù)

    注:**表示各相關(guān)系數(shù)通過(guò)99%的顯著性檢驗(yàn);*表示各相關(guān)系數(shù)通過(guò)95%的顯著性檢驗(yàn)。

    相關(guān)分析表明(表1),本研究選取的所有變率指標(biāo)均與其他某個(gè)或某些變率之間存在顯著相關(guān),其中ENSO指標(biāo)作為我國(guó)預(yù)報(bào)減災(zāi)的重點(diǎn),其相關(guān)范疇最廣,有4個(gè)相關(guān)指標(biāo);而印度洋偶極子和北大西洋濤動(dòng)的相關(guān)范圍則較小,主要分別與ENSO和北極濤動(dòng)相關(guān),因此在探討其與絲路海區(qū)海表溫度異常變化的相關(guān)性時(shí)對(duì)其他變率信號(hào)的干擾應(yīng)予以濾除。下文均以上述得到的各變率的相關(guān)性為基礎(chǔ),先對(duì)各指標(biāo)序列進(jìn)行標(biāo)準(zhǔn)化處理,然后就某一變率的原始指標(biāo)減去通過(guò)回歸分析得到的其他變率在該變率中的信號(hào),以達(dá)到對(duì)其他干擾信號(hào)的濾除;同時(shí)為消除全球變化的影響(圖1),首先對(duì)各物理量和指標(biāo)序列作去傾處理,然后進(jìn)行相關(guān)分析。此外,南海夏季風(fēng)與東亞夏季風(fēng)具有高度相關(guān)(相關(guān)系數(shù)為0.898),而南海夏季風(fēng)作為東亞夏季風(fēng)的肇始,與我國(guó)雨型及旱澇等災(zāi)害分布密切相關(guān)[9,24]。限于篇幅,在季風(fēng)指數(shù)、局地海洋變率和遙相關(guān)型等氣候變率中,本研究主要圍繞南海夏季風(fēng)、ENSO、PDO和AO與絲路海區(qū)海表溫度的相關(guān)關(guān)系特征進(jìn)行探討。

    需要指出的是,相關(guān)關(guān)系是變量(要素)間不確定的依存關(guān)系,必須在一定的物理等意義下統(tǒng)籌考慮,才能認(rèn)識(shí)和理解相關(guān)關(guān)系背后的機(jī)制機(jī)理;相關(guān)關(guān)系和機(jī)制機(jī)理本身具有復(fù)雜性,而認(rèn)識(shí)相關(guān)性是研究事物變化機(jī)制機(jī)理的根本前提,尤其本研究以相關(guān)性分析的技術(shù)方法及應(yīng)用為側(cè)重點(diǎn),因此關(guān)于相關(guān)性首先提出以下概念。在涉及第三方C(可能為要素或變率集合)的影響下,分析要素A與變率B的相關(guān)性主要考慮3種類型(圖2):①直接相關(guān);②間接(遞接)相關(guān);③協(xié)應(yīng)(對(duì)接)相關(guān),即A和B沒有直接和間接關(guān)系,而是共同對(duì)C作出響應(yīng)。需要注意的是,實(shí)際的相關(guān)要復(fù)雜得多,一般A和B二者的相關(guān)主要看作3種相關(guān)的疊加,而上述對(duì)各要素的去傾處理則是消除第三種相關(guān)的影響。通常,研究中采用的偏相關(guān)分析是針對(duì)如圖2(a)所示的情況,即如果A和C也存在相關(guān)(②存在的情況下),對(duì)A和B的相關(guān)性則可用偏相關(guān)分析。本研究主要關(guān)注第一種和第二種相關(guān),目的在于找到要素和變率間的直接關(guān)系或作用通道,從而為分析探討要素變化背后的機(jī)理以及對(duì)其進(jìn)行預(yù)測(cè)預(yù)報(bào)提供必要的科學(xué)依據(jù)。

    圖2 相關(guān)類型示意(箭頭和連線表示作用關(guān)系和途徑)

    2 南海夏季風(fēng)與絲路海區(qū)海表溫度異常變化的相關(guān)性

    首先對(duì)南海夏季風(fēng)標(biāo)準(zhǔn)化指數(shù)與絲路海區(qū)海表溫度異常作相關(guān)分析(簡(jiǎn)稱原始相關(guān)),然后將濾除南亞夏季風(fēng)信號(hào)的南海夏季風(fēng)指數(shù)與絲路海區(qū)海表溫度異常作相關(guān)分析(簡(jiǎn)稱濾波相關(guān)),最后將濾除南亞夏季風(fēng)信號(hào)的南海夏季風(fēng)指數(shù)與絲路海區(qū)海表溫度異常作時(shí)滯相關(guān)分析(簡(jiǎn)稱時(shí)滯相關(guān))(圖3)。原始相關(guān)分析表明,南海夏季風(fēng)主要與熱帶的西南印度洋、中太平洋以及包括我國(guó)近海、日本海以及日本以東的北太平洋3個(gè)海域的海表溫度呈正相關(guān)關(guān)系,而與南海中部、珊瑚海等海域存在負(fù)相關(guān)關(guān)系(圖3(a))。比較而言,濾除信號(hào)后,正相關(guān)區(qū)域的范圍均有明顯縮小,其中我國(guó)近海海域的正相關(guān)關(guān)系已不明顯,而日本以東的北太平洋海域正相關(guān)范圍變化較?。欢叩呢?fù)相關(guān)區(qū)域范圍擴(kuò)大,其中珊瑚海的負(fù)相關(guān)區(qū)域擴(kuò)大至澳洲北部東西兩側(cè)的大范圍海域??梢钥闯?,由于南亞夏季風(fēng)參與影響,南海夏季風(fēng)與絲路海區(qū)海表溫度正相關(guān)區(qū)域范圍擴(kuò)大,而負(fù)相關(guān)區(qū)域則相對(duì)減少(圖3(b))。時(shí)滯相關(guān)分析表明,前期主要是20°N以北的北太平洋海域海表溫度對(duì)南海夏季風(fēng)可能存在正相關(guān)影響(圖3(c)),而后期南海夏季風(fēng)則可能對(duì)我國(guó)近海及澳洲北部等大范圍海域存在負(fù)相關(guān)影響(圖3(d)),這可能與季風(fēng)的增強(qiáng)會(huì)增加海氣熱通量的釋放有關(guān)。以上差異應(yīng)在南海夏季風(fēng)變化及影響的預(yù)測(cè)預(yù)報(bào)中予以關(guān)注。

    圖3 南海夏季風(fēng)與絲路海區(qū)海表溫度相關(guān)場(chǎng)

    3 海洋年際、年代際變率與絲路海區(qū)海表溫度異常變化的相關(guān)性

    ENSO(主要指暖事件)與絲路海區(qū)海表溫度異常的原始相關(guān)性分析表明,其主要與我國(guó)南海和印度洋海表溫度呈正相關(guān)關(guān)系,而與熱帶西太平洋海域海表溫度呈負(fù)相關(guān)關(guān)系(圖4(a))。濾除信號(hào)后,除印度洋外,二者的正相關(guān)區(qū)域擴(kuò)大到包括我國(guó)近海及以北太范圍的西北太平洋海域;同時(shí)赤道以北的熱帶西太平洋負(fù)相關(guān)區(qū)域相關(guān)性更加明顯,赤道以南有所減弱(圖4(b))。時(shí)滯相關(guān)分析表明,前期印度洋和熱帶(20°S~20°N)西太平洋分別與ENSO呈反位相相關(guān)關(guān)系(圖4(c)),大尺度的環(huán)流系統(tǒng)和局地反饋可能是ENSO爆發(fā)的前提,這方面的研究仍有待深入,如中東太平洋的西風(fēng)異常(爆發(fā))與ENSO的關(guān)聯(lián)性等;ENSO對(duì)后期絲路海區(qū)海表溫度的影響主要體現(xiàn)在ENSO發(fā)生后可能對(duì)包括我國(guó)南海以及東印度洋海域海表溫度產(chǎn)生負(fù)相關(guān)影響,而對(duì)西印度洋,我國(guó)渤海、黃海、東海以及日本以東的北太平洋產(chǎn)生正相關(guān)影響(圖4(d))。

    圖4 ENSO與絲路海區(qū)海表溫度異常相關(guān)場(chǎng)

    從PDO與絲路海區(qū)海表溫度異常的原始相關(guān)分析結(jié)果可以看出(圖5(a)),在熱帶海域存在一定程度的ENSO信號(hào),在濾除包括ENSO等變率的影響后,PDO對(duì)絲路海區(qū)海表溫度的相關(guān)關(guān)系呈現(xiàn)約以20°N和40°N為界的南北三極反位相相關(guān)關(guān)系(圖5(b))。時(shí)滯相關(guān)分析表明,前期西太平洋海域海表溫度異常與PDO呈正負(fù)相間的帶狀相關(guān)性分布(圖5(c)),而PDO對(duì)后期包括我國(guó)南海、印度洋和赤道中太平洋海表溫度呈現(xiàn)正相關(guān)影響,對(duì)20°N~40°N北太平洋海域海表溫度則呈負(fù)相關(guān)影響;我國(guó)渤海、黃海、東海對(duì)PDO的響應(yīng)不明顯,也可能與信號(hào)傳播和影響的時(shí)間尺度有關(guān)(圖5(d))。

    圖5 PDO與絲路海區(qū)海表溫度異常相關(guān)場(chǎng)

    4 北極濤動(dòng)與絲路海區(qū)海表溫度異常變化的相關(guān)性

    北極濤動(dòng)是中緯度地區(qū)與極地間的一種遙相關(guān)振蕩,除日本海外,其與絲路海區(qū)海表溫度異常變化的相關(guān)性并不明顯(圖6(b));而由于北大西洋濤動(dòng)和PDO等其他變率的影響,AO對(duì)包括我國(guó)渤海、黃海、東海和30°N以北的北太平洋海域海表溫度異常呈現(xiàn)顯著的正相關(guān)關(guān)系(圖6(a))。時(shí)滯相關(guān)分析表明,前期40°N以北的北太平洋海域?qū)O可能存在主要的正相關(guān)影響,同時(shí)赤道以南的印度洋及鄰近海域海表溫度異常對(duì)AO也可能存在較弱的正相關(guān)影響(圖6(c));AO的存在主要對(duì)我國(guó)南海、印度洋全海域海表溫度異常存在正相關(guān)影響(圖6(d)),這可能通過(guò)副熱帶高壓和西風(fēng)急流等大氣波列產(chǎn)生作用,具體影響機(jī)理需進(jìn)一步分析探討。

    圖6 AO與絲路海區(qū)海表溫度異常相關(guān)場(chǎng)

    5 討論

    近幾十年來(lái),全球變化以增溫為顯著特征[2],而不同區(qū)域?qū)θ蜃兣捻憫?yīng)并非一致(圖7(a))。一般來(lái)說(shuō),南極、北極、青藏高原(由于其隆起成為地球中低緯度地區(qū)唯一大面積進(jìn)入冰凍圈的地區(qū))以及赤道附近暖池區(qū)域被認(rèn)為是地球的4個(gè)極,主要對(duì)應(yīng)冷極和暖極[5];在氣候變暖背景下,由于冰凍圈及其快速變化對(duì)全球和區(qū)域的氣候變化有至關(guān)重要的影響,作為冷極的極地地區(qū)和青藏高原被學(xué)者稱為全球和區(qū)域氣候變化的驅(qū)動(dòng)器和放大器,這里有關(guān)極地地區(qū)氣候變化的討論從略。另一方面,就全球海洋熱力狀況的平均態(tài)及其要素本身變化而言(圖7(a)、圖7(b)),熱力要素的梯度則體現(xiàn)為氣候變化響應(yīng)的放大器,梯度的大值區(qū)(圖7(c)中方格填充區(qū))成為對(duì)氣候變暖響應(yīng)的最顯著區(qū),在海洋中主要位于我國(guó)東海和日本海、墨西哥灣流、南美洲及澳洲的東南部沿海等中高緯度區(qū)域;同時(shí)由于溫度梯度的存在,在同等響應(yīng)下(變化幅度相同),中高緯度地區(qū)的氣候變化速度較小,而低緯度地區(qū)的變化速度較大(圖7(c)中斜線填充區(qū)),進(jìn)而表明溫度梯度在全球氣候變化一致性和氣候變化區(qū)域響應(yīng)的非均勻性上起著關(guān)鍵性作用,尤其是在氣候變化區(qū)域響應(yīng)上充分體現(xiàn)放大器的效用。更重要的是,區(qū)域性和局地氣候因子及其變率不僅交互影響,而且還制約著氣候平均態(tài)的構(gòu)成,從而進(jìn)一步調(diào)控全球變化和氣候區(qū)域響應(yīng)特征。

    圖7 1950—2015年間全球(不包括極地及 鄰近區(qū)域)地表溫度氣候特征

    鑒于全球氣候變化的一致性和氣候變化區(qū)域響應(yīng)的非均勻性,區(qū)域性氣候和海洋環(huán)境變化研究及其預(yù)測(cè)預(yù)報(bào)更具現(xiàn)實(shí)意義。加之全球變暖背景下,區(qū)域氣候和海洋變化變率(尤其是急轉(zhuǎn)和突變等)及災(zāi)害風(fēng)險(xiǎn)的不確定性(主要指災(zāi)害本身發(fā)生頻次的不確定性,尤其重大災(zāi)害的突發(fā)性和時(shí)空密集性)會(huì)進(jìn)一步加強(qiáng),從而使氣候和環(huán)境要素的變化、轉(zhuǎn)變更為復(fù)雜,極端天氣氣候事件突發(fā),這會(huì)給預(yù)測(cè)預(yù)報(bào)和減災(zāi)防災(zāi)工作帶來(lái)很大挑戰(zhàn),而厘清相關(guān)要素與不同時(shí)空尺度氣候變率的凈相關(guān)性是科學(xué)問題解決及應(yīng)用的重要前提。因此,本研究選取一些典型的海洋和大氣中的氣候變率,從線性系統(tǒng)的角度(不考慮非線性相關(guān)),基于區(qū)域氣候因子(變率)間一定物理意義上的相互影響和作用關(guān)系的認(rèn)識(shí),在暫不考慮時(shí)滯效應(yīng)的前提下,分析其與絲路海區(qū)海表溫度異常變化的凈相關(guān)性特征,經(jīng)濾除干擾信號(hào)后的主要變率與海表溫度異常變化之間存在顯著相關(guān)。尤其值得提出的是,分析表明各氣候變率對(duì)海表溫度變化可能存在交互或聯(lián)合影響,是否濾除其他變率的干擾對(duì)相關(guān)性有時(shí)會(huì)存在較大程度的影響,而其間的差異很可能是氣候變率與環(huán)境要素之間相關(guān)性時(shí)空和時(shí)頻改變(如二者相關(guān)性空間分布變化尤其是顯著性區(qū)域的轉(zhuǎn)移、強(qiáng)度的改變乃至位相轉(zhuǎn)換和時(shí)間突變等)的重要原因之一,更是氣候因子及其變率指標(biāo)體系構(gòu)建的根本依據(jù)。這一點(diǎn)在研究氣候變率影響機(jī)理及其氣候因子和海洋環(huán)境變化的預(yù)測(cè)預(yù)報(bào)時(shí)應(yīng)予以關(guān)注,需作必要的干擾信號(hào)濾除處理,從而以更清晰的物理圖景認(rèn)識(shí)海表溫度變化和未來(lái)趨勢(shì),更好地為相關(guān)氣候變化分析評(píng)估和未來(lái)情景預(yù)估提供科學(xué)依據(jù),并逐步提升科學(xué)應(yīng)對(duì)氣候變化和海洋環(huán)境安全保障水平,最終為我國(guó)海洋強(qiáng)國(guó)戰(zhàn)略布局和實(shí)施服務(wù)。

    同時(shí),本研究對(duì)海表溫度與氣候變率之間相互作用關(guān)系的認(rèn)識(shí)還不是很透徹,尚待進(jìn)一步圍繞某一變率對(duì)海洋環(huán)境要素的影響進(jìn)行細(xì)致分析,進(jìn)而明確其間的關(guān)聯(lián)性以及相互作用關(guān)系和影響過(guò)程機(jī)理等。此外,本研究沒有采用偏相關(guān)分析,主要考慮要素空間分布及與變率相關(guān)的區(qū)域特征和復(fù)雜度,其可操作性較差;與控制變量的偏相關(guān)分析相比較,本研究首先提出分類型的相關(guān)分析和信號(hào)濾除更能了解變量間的相關(guān)性,尤其是相互作用關(guān)系、影響途徑和通道及其貢獻(xiàn)等,同時(shí)結(jié)合一定的物理意義可更清楚地認(rèn)識(shí)要素變化的機(jī)制機(jī)理;即使使用偏相關(guān)分析,仍需對(duì)偏相關(guān)中的相關(guān)性進(jìn)行信號(hào)濾除處理,而一般的偏相關(guān)缺乏這一步驟。

    6 結(jié)語(yǔ)

    本研究側(cè)重于相關(guān)性分析技術(shù)方法與應(yīng)用,旨在辨識(shí)氣候變化背景下各變率對(duì)包括我國(guó)及周邊等絲路海區(qū)海洋環(huán)境變化的交互影響,經(jīng)一定程度地(有限但必要)濾除干擾信號(hào),探討絲路海區(qū)海表溫度異常變化與各氣候變率之間的凈相關(guān)特征和顯著差異,以期為提高我國(guó)海洋環(huán)境變化的診斷和預(yù)報(bào)水平提供必要的參考依據(jù),進(jìn)而提升我國(guó)海洋領(lǐng)域科學(xué)應(yīng)對(duì)氣候變化以及海洋環(huán)境安全保障水平。主要結(jié)論如下:

    (1)南海夏季風(fēng)主要與中國(guó)近海(尤其是南海)和澳洲北部海域呈負(fù)相關(guān)關(guān)系,而與日本及以東的北太平洋海域呈正相關(guān)關(guān)系。由于南亞季風(fēng)的影響,可能會(huì)增強(qiáng)南海夏季風(fēng)與我國(guó)近海海表溫度的正相關(guān)關(guān)系(主要包括相關(guān)關(guān)系的范圍、強(qiáng)度和位相等),并可能使我國(guó)海域與南海夏季風(fēng)呈現(xiàn)反位相的相關(guān)關(guān)系。前期大致以20°N為界,絲路海區(qū)海表溫度對(duì)南海夏季風(fēng)呈南北反位相的相關(guān)影響,尤其是20°N以北海域存在明顯的正相關(guān)關(guān)系;而南海夏季風(fēng)對(duì)后期我國(guó)近海和澳洲北部海域海表溫度的降低可能存在明顯影響,這可能與季風(fēng)引起的海氣相互作用及熱通量交換增強(qiáng)有關(guān),尚待進(jìn)一步深入分析。

    (2)ENSO(暖事件)主要與印度洋和包括我國(guó)近海的西北太平洋呈正相關(guān)關(guān)系,而與熱帶西太平洋呈現(xiàn)負(fù)相關(guān)關(guān)系;由于其他變率的影響,可能會(huì)明顯減弱ENSO對(duì)包括我國(guó)近海的西北太平洋海域海表溫度的正相關(guān)關(guān)系。前期印度洋和熱帶西太平洋對(duì)ENSO可能主要呈反位相的相關(guān)影響,而ENSO主要對(duì)后期印度洋偶極子式的海表溫度異常變化以及包括我國(guó)近海的西太平洋海域海表溫度異常反位相變化產(chǎn)生影響。PDO與絲路海區(qū)海表溫度主要呈現(xiàn)南北正負(fù)相間的帶狀相關(guān)性分布,ENSO信號(hào)在PDO與絲路海區(qū)海表溫度相關(guān)性中的影響可能較大。前期西太平洋海表溫度南北正負(fù)相間的帶狀分布對(duì)PDO可能存在主要影響,PDO對(duì)后期印度洋海盆尺度的海表溫度異常變化可能主要存在一致的正相關(guān)影響,同樣對(duì)我國(guó)南海也存在正相關(guān)影響,而西太平洋仍以正負(fù)相間的海表溫度變化作為對(duì)PDO的主要響應(yīng)特征。此外,我國(guó)渤海、黃海和東海對(duì)PDO的后期響應(yīng)不明顯,也可能與信號(hào)未被捕捉有關(guān)。

    (3)北極濤動(dòng)與絲路海區(qū)海表溫度的關(guān)系較不明顯,由于其他變率的聯(lián)合影響,包括我國(guó)渤海、黃海和東海以及30°N以北的北太平洋海域海表溫度異常與AO呈現(xiàn)顯著的正相關(guān)關(guān)系。前期40°N以北的北太平洋海域?qū)O可能存在主要的正相關(guān)影響,而AO對(duì)我國(guó)南海、印度洋全海域海表溫度后期的異常變化可能存在正相關(guān)影響。

    (4)除極地區(qū)域外,就全球熱力狀況的平均態(tài)及其基礎(chǔ)要素本身變化而言,熱力要素的空間梯度則體現(xiàn)為其對(duì)氣候變化區(qū)域響應(yīng)的放大器,其中溫度梯度在全球氣候變化一致性和氣候變化區(qū)域響應(yīng)的非均勻性上很可能起著關(guān)鍵性作用。尤其是梯度的大值區(qū)成為對(duì)氣候變暖響應(yīng)的最顯著區(qū),在海洋中主要位于我國(guó)東海和日本海、墨西哥灣流、南美洲和非洲以至澳洲的東南部沿海等中高緯度區(qū)域。而溫度梯度的分布和變化與區(qū)域性和局地氣候因子及變率具有直接關(guān)系(如海溫梯度與季風(fēng)等),具體機(jī)制機(jī)理將進(jìn)一步分析研究。

    鑒于海-陸-氣系統(tǒng)內(nèi)各要素和變率相關(guān)性異常復(fù)雜,本研究并非驗(yàn)證氣候變率相關(guān)指標(biāo)指數(shù)的適用性和正確性以及其被用于有關(guān)機(jī)制機(jī)理研究的合理性,今后將致力于從聚類分析的角度,重點(diǎn)辨識(shí)全球和區(qū)域性海洋、陸地和大氣中的系統(tǒng)性結(jié)構(gòu),著眼于純粹相關(guān)性構(gòu)建氣候變率指標(biāo)體系,綜合分析研究各體系間的(包括非線性)相互作用關(guān)系和機(jī)制,以豐富海-陸-氣相互作用與氣候環(huán)境變化和異常演變的機(jī)理研究,尤其為氣候和海洋變化的預(yù)測(cè)預(yù)報(bào)和災(zāi)害風(fēng)險(xiǎn)評(píng)估及適應(yīng)和對(duì)策研究提供必要的科學(xué)依據(jù)和技術(shù)支撐,也為我國(guó)21世紀(jì)海上絲綢之路建設(shè)和海洋強(qiáng)國(guó)戰(zhàn)略實(shí)施提供環(huán)境安全保障服務(wù)。

    [1] 吳國(guó)雄,李建平,周天軍,等.影響我國(guó)短期氣候異常的關(guān)鍵區(qū):亞印太交匯區(qū)[J].地球科學(xué)進(jìn)展,2006,21(11):1109-1117.

    [2] 蔡榕碩,齊慶華.氣候變化與全球海洋:影響、適應(yīng)和脆弱性評(píng)估之解讀[J].氣候變化研究進(jìn)展,2014,10(3):185-190.

    [3] 翁學(xué)傳,張啟龍,顏廷壯.熱帶西太平洋暖池域次表層水熱含量變化及其與我國(guó)東部汛期降水和副高的相關(guān)關(guān)系[J].海洋科學(xué)集刊,1996,37:1-9.

    [4] 張啟龍,翁學(xué)傳.熱帶西太平洋暖池的某些海洋學(xué)特征分析[J].海洋科學(xué)集刊,1997,38:31-38.

    [5] 齊慶華,蔡榕碩.西太平洋暖池水體振蕩與赤道中東太平洋海表溫度異常的關(guān)聯(lián)性初探[J].海洋科學(xué),2015,39(9):77-85.

    [6] 翁學(xué)傳,張啟龍,楊玉玲,等.東海黑潮熱輸送及其與黃淮平原區(qū)汛期降水的關(guān)系[J].海洋與湖沼,1996,27(3):237-245.

    [7] 齊慶華,蔡榕碩,張啟龍.臺(tái)灣以東黑潮熱輸送變異及可能的氣候效應(yīng)[J].海洋學(xué)報(bào),2012,34(5):31-38.

    [8] 黃榮輝,顧雷,徐予紅,等.東亞夏季風(fēng)爆發(fā)和北進(jìn)的年際變化特征及其與熱帶西太平洋熱狀況的關(guān)系[J].大氣科學(xué),2005,29(1):20-36.

    [9] 陳際龍,黃榮輝.亞澳季風(fēng)各子系統(tǒng)氣候?qū)W特征的異同研究I.夏季風(fēng)流場(chǎng)結(jié)構(gòu)[J].大氣科學(xué),2006,30(6):1 091-1 102.

    [10] 李崇印,穆明權(quán).赤道印度洋海溫偶極子型振蕩及其氣候影響[J].大氣科學(xué),2001,25(4):433-443.

    [11] 李若鈍,武建平,泮惠周,等.ENSO與中國(guó)氣候[J].海洋通報(bào),1988,7(3):7-12.

    [12] FRANCIS R C,HARE S R.Decadal-scale regime shifts in the large marine ecosystems of the Northeast Pacific:a case for historical science[J].Fish Oceanogr,1994,3:279-291.

    [13] NITTA T.Convective activities in the tropical western Pacific and their impact on the Northern Hemisphere summer circulation[J].J.Meteor.Soc.Japan,1987,65:373-390.

    [14] CAI R,TAN H,QI Q.Impacts of and adaptation to inter-decadal marine climate change in coastal China seas[J].International Journal of Climatology,2015,DOI:10.1002/joc.4591.

    [15] 宗海鋒,張慶云,陳烈庭.東亞-太平洋遙相關(guān)型形成過(guò)程與ENSO盛期海溫關(guān)系的研究[J].大氣科學(xué),2008,32(2):220-230.

    [16] 孫建奇,王會(huì)軍.北極濤動(dòng)與太平洋年代際振蕩的關(guān)系[J].科學(xué)通報(bào),2005,50(15):1648-1653.

    [17] BARROS V R,F(xiàn)IELD C B,DOKKEN D J,et al.Climate Change 2014:Impacts,Adaptation,and Vulnerability.Part B:Regional Aspects[M].Cambridge,United Kingdom and New York,USA:Cambridge University Press,2014:1133-1197.

    [18] LI J P,ZENG Q C.A new monsoon index,its interannual variability and relation with monsoon precipitation[J].Climatic and Environmental Research,2005,10(3):351-365.

    [19] BARNSTON A G,ROPELEWSKI C F.Prediction of ENSO episodes using canonical correlation analysis[J].J.Climate,1992,5:1316-1345.

    [20] MANTUA N J,HARE S R,ZHANG Y,et al.A Pacific interdecadal climate oscillation with impacts on salmon production[J].Bull.Amer.Meteor.Soc.,1997,78:1069-1079.

    [21] SAJI N H,GOSWAMI B N,VINAYACHANDRAN P N,et al.A dipole mode in the tropical Indian Ocean[J].Nature,1999,401:360-363.

    [22] ZHOU S A,MILLER J,WANG J,et al.Trends of NAO and AO and their associations with stratospheric processes[J].Geophys.Res.Lett.,2001,28:4107-4110.

    [23] BARNSTON A G,LIVEZEY R E.Classification,seasonality and persistence of low-frequency atmospheric circulation patterns[J].Mon.Wea.Rev.,1987,115:1083-1126.

    [24] 齊慶華,蔡榕碩.南海海表溫時(shí)空演變與南海夏季風(fēng)爆發(fā)早晚相關(guān)性初探[J].海洋學(xué)報(bào),2014,36(3):94-103.

    Spatio-Temporal Change of Sea Surface Temperature Anomalies in Seas of 21st Century Maritime Silk Road and Its Net Correlation to Climate Variability

    QI Qinghua1,2,CAI Rongshuo1,2

    (1.Third Institute of Oceanography,SOA,Xiamen 361005,China;2.Key Laboratory of Global Change and Marine-Atmospheric Chemistry,SOA,Xiamen 361005,China)

    Regional ocean system and its thermal environmental elements are affected reciprocally by local and other teleconnected climate variability on different spatio-temporal scales.Analysis of the net relationship between them is one of the important aspects for the research of mechanisms and predictability on changes of marine environmental elements.This paper chose the main oceanic and atmospheric variability and teleconnection patterns with filtering interference signals to a certain degree,and analyzed the main correlation characteristics and difference between the climate variability and sea surface temperature (SST) anomalies in seas of 21st Century Maritime Silk Road (MSRS).The analysis provides a new basis and technical support services for the research on the ocean and climate changes and the marine environment security of China.

    Sea surface temperature,Climate variability,Temperature gradient,Marine environment security,21st Century Maritime Silk Road

    2016-08-31;

    2017-03-11

    中國(guó)清潔生產(chǎn)發(fā)展機(jī)制基金項(xiàng)目(2014112);福建省自然科學(xué)基金面上項(xiàng)目(2017J01076);國(guó)家海洋局第三海洋研究所基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(海三科20015030).

    齊慶華,助理研究員,博士,研究方向?yàn)楹j憵庀嗷プ饔门c氣候、環(huán)境變化及災(zāi)害風(fēng)險(xiǎn)評(píng)估與適應(yīng)對(duì)策,電子信箱:qqh_2002@163.com

    P7

    A

    1005-9857(2017)04-0041-09

    猜你喜歡
    變率海表海區(qū)
    內(nèi)部變率和全球變暖對(duì)春季北太平洋維多利亞模態(tài)增強(qiáng)的相對(duì)貢獻(xiàn)
    研究顯示降水變率將隨氣候增暖而增強(qiáng)
    不正規(guī)半日潮海區(qū)高(低)潮選取的探討
    基于無(wú)人機(jī)的海表環(huán)境智能監(jiān)測(cè)系統(tǒng)設(shè)計(jì)與應(yīng)用
    2016與1998年春季北大西洋海表溫度異常的差異及成因
    融合海表溫度產(chǎn)品在渤黃東海的對(duì)比分析及初步驗(yàn)證
    太陽(yáng)總輻照度對(duì)熱帶中太平洋海表溫度年代際變化的可能影響
    地鐵長(zhǎng)大過(guò)海區(qū)間隧道人員疏散模擬分析
    Does a monsoon circulation exist in the upper troposphere over the central and eastern tropical Pacifc?
    3個(gè)不同海區(qū)鳶烏賊漁業(yè)生物學(xué)的初步比較
    女人爽到高潮嗷嗷叫在线视频| 啦啦啦观看免费观看视频高清| 国产欧美日韩一区二区精品| 亚洲成人久久爱视频| e午夜精品久久久久久久| 日韩欧美在线二视频| 黑人欧美特级aaaaaa片| 天天一区二区日本电影三级| 妹子高潮喷水视频| 99re在线观看精品视频| 国产精品自产拍在线观看55亚洲| 在线观看午夜福利视频| 久久草成人影院| 国产高清videossex| 精品久久久久久久人妻蜜臀av| 国产黄片美女视频| 亚洲乱码一区二区免费版| 1024香蕉在线观看| 极品教师在线免费播放| 嫁个100分男人电影在线观看| 99久久国产精品久久久| 狂野欧美激情性xxxx| 午夜精品一区二区三区免费看| 日本熟妇午夜| 人人妻人人澡欧美一区二区| 久久这里只有精品中国| 老汉色av国产亚洲站长工具| 亚洲av熟女| 两个人看的免费小视频| 久久 成人 亚洲| 成人手机av| 欧美色欧美亚洲另类二区| 精品久久久久久久久久免费视频| 婷婷亚洲欧美| 五月伊人婷婷丁香| 久久人妻福利社区极品人妻图片| 很黄的视频免费| 正在播放国产对白刺激| 岛国在线观看网站| 午夜福利18| 十八禁网站免费在线| 99国产精品一区二区蜜桃av| 国产69精品久久久久777片 | 成熟少妇高潮喷水视频| 久久精品国产综合久久久| 老司机在亚洲福利影院| av国产免费在线观看| 国产一区二区激情短视频| 日日干狠狠操夜夜爽| 欧美日韩福利视频一区二区| 两个人看的免费小视频| 一区二区三区激情视频| 午夜福利欧美成人| 国产高清有码在线观看视频 | 精品不卡国产一区二区三区| 日韩成人在线观看一区二区三区| 久久中文字幕一级| 欧美乱妇无乱码| 老司机福利观看| 日本在线视频免费播放| 看黄色毛片网站| av福利片在线观看| 国产蜜桃级精品一区二区三区| 一个人免费在线观看电影 | 亚洲av成人一区二区三| 成人av一区二区三区在线看| 亚洲精品美女久久久久99蜜臀| 亚洲国产精品成人综合色| 亚洲成a人片在线一区二区| 精品国内亚洲2022精品成人| 三级国产精品欧美在线观看 | 1024视频免费在线观看| svipshipincom国产片| 日本在线视频免费播放| 一二三四在线观看免费中文在| 美女大奶头视频| 亚洲狠狠婷婷综合久久图片| aaaaa片日本免费| 一级毛片精品| 欧美日本亚洲视频在线播放| 制服丝袜大香蕉在线| 久久这里只有精品19| 国产精品电影一区二区三区| 人人妻,人人澡人人爽秒播| 午夜福利免费观看在线| 最新在线观看一区二区三区| 亚洲最大成人中文| 少妇人妻一区二区三区视频| 天天添夜夜摸| 亚洲欧美精品综合久久99| 久久精品亚洲精品国产色婷小说| 色噜噜av男人的天堂激情| 最近最新中文字幕大全免费视频| 国产片内射在线| 在线观看免费视频日本深夜| 亚洲熟妇熟女久久| 国产成人av教育| 久久久精品欧美日韩精品| 一a级毛片在线观看| 国产黄色小视频在线观看| 久久精品国产99精品国产亚洲性色| 韩国av一区二区三区四区| 亚洲av成人不卡在线观看播放网| а√天堂www在线а√下载| 成熟少妇高潮喷水视频| 亚洲人成网站在线播放欧美日韩| 19禁男女啪啪无遮挡网站| 欧美日韩一级在线毛片| 国产精品亚洲美女久久久| 无遮挡黄片免费观看| 免费观看人在逋| 狂野欧美激情性xxxx| 亚洲va日本ⅴa欧美va伊人久久| 无遮挡黄片免费观看| 亚洲av熟女| 日本在线视频免费播放| 99久久久亚洲精品蜜臀av| 欧美日韩中文字幕国产精品一区二区三区| 国产单亲对白刺激| 免费无遮挡裸体视频| 日本 av在线| 又大又爽又粗| 波多野结衣巨乳人妻| 人人妻,人人澡人人爽秒播| 丝袜美腿诱惑在线| 日韩欧美在线二视频| 亚洲中文字幕日韩| 嫩草影院精品99| 91麻豆av在线| 国产亚洲精品久久久久久毛片| 男人舔女人下体高潮全视频| 亚洲九九香蕉| 18禁黄网站禁片午夜丰满| 人妻久久中文字幕网| 人妻丰满熟妇av一区二区三区| 天天一区二区日本电影三级| 亚洲精品av麻豆狂野| 可以在线观看毛片的网站| 国产精品1区2区在线观看.| 国产精品久久久久久亚洲av鲁大| 欧美黑人欧美精品刺激| 国产成人av教育| 国产成+人综合+亚洲专区| 国产在线精品亚洲第一网站| 波多野结衣高清无吗| 亚洲真实伦在线观看| 午夜亚洲福利在线播放| 99久久精品国产亚洲精品| 国产熟女xx| 亚洲精品在线美女| 黑人欧美特级aaaaaa片| 国内揄拍国产精品人妻在线| 国产亚洲av高清不卡| 精品一区二区三区四区五区乱码| 99国产综合亚洲精品| 91国产中文字幕| 午夜福利视频1000在线观看| 亚洲国产精品成人综合色| 亚洲av日韩精品久久久久久密| 国产伦人伦偷精品视频| 村上凉子中文字幕在线| 动漫黄色视频在线观看| 国产免费av片在线观看野外av| 亚洲美女视频黄频| 国产野战对白在线观看| 波多野结衣高清作品| 亚洲国产欧洲综合997久久,| 一进一出好大好爽视频| 他把我摸到了高潮在线观看| 亚洲专区中文字幕在线| 在线十欧美十亚洲十日本专区| 精品一区二区三区视频在线观看免费| 日韩 欧美 亚洲 中文字幕| 亚洲性夜色夜夜综合| 国产成人精品久久二区二区91| x7x7x7水蜜桃| 脱女人内裤的视频| 亚洲国产精品sss在线观看| 中文字幕熟女人妻在线| 99久久无色码亚洲精品果冻| 可以在线观看毛片的网站| 久久国产乱子伦精品免费另类| 中文字幕熟女人妻在线| 中文字幕精品亚洲无线码一区| 中文资源天堂在线| 亚洲乱码一区二区免费版| 一个人免费在线观看的高清视频| 国产亚洲欧美98| 人人妻人人看人人澡| 日韩欧美 国产精品| 久久精品夜夜夜夜夜久久蜜豆 | 一级毛片女人18水好多| 俺也久久电影网| 变态另类丝袜制服| av欧美777| 久久久久久久精品吃奶| a级毛片在线看网站| 美女 人体艺术 gogo| 午夜精品一区二区三区免费看| 精品久久蜜臀av无| 天堂影院成人在线观看| 看片在线看免费视频| 脱女人内裤的视频| 欧美人与性动交α欧美精品济南到| 一边摸一边抽搐一进一小说| cao死你这个sao货| 日本一区二区免费在线视频| 午夜精品一区二区三区免费看| 一本综合久久免费| 国产精品1区2区在线观看.| 亚洲电影在线观看av| 国产精品久久久人人做人人爽| 香蕉丝袜av| 丰满人妻熟妇乱又伦精品不卡| 欧洲精品卡2卡3卡4卡5卡区| 国产99白浆流出| 亚洲,欧美精品.| 久久久久久亚洲精品国产蜜桃av| 亚洲av五月六月丁香网| 欧美在线黄色| 午夜福利在线观看吧| 精品久久久久久成人av| 老司机在亚洲福利影院| 我要搜黄色片| 精品第一国产精品| 欧美日本亚洲视频在线播放| 麻豆国产av国片精品| 夜夜躁狠狠躁天天躁| 一级黄色大片毛片| 精品国产乱码久久久久久男人| 久久九九热精品免费| 一个人观看的视频www高清免费观看 | 日本在线视频免费播放| 人人妻人人澡欧美一区二区| cao死你这个sao货| 久久午夜综合久久蜜桃| 国产野战对白在线观看| 人成视频在线观看免费观看| 欧美大码av| 亚洲av成人不卡在线观看播放网| 日韩有码中文字幕| 国产一区二区三区在线臀色熟女| 精品福利观看| 国产精品影院久久| 美女免费视频网站| 亚洲中文字幕一区二区三区有码在线看 | 国产成人精品久久二区二区免费| 观看免费一级毛片| 精品国内亚洲2022精品成人| 亚洲国产日韩欧美精品在线观看 | 亚洲av中文字字幕乱码综合| 国内精品一区二区在线观看| 国产精品免费视频内射| 又紧又爽又黄一区二区| 九九热线精品视视频播放| 国产不卡一卡二| 亚洲18禁久久av| 久久久久久国产a免费观看| 听说在线观看完整版免费高清| 亚洲国产欧美一区二区综合| 国产成人欧美在线观看| 欧美一级a爱片免费观看看 | 亚洲电影在线观看av| 亚洲中文字幕一区二区三区有码在线看 | 国产伦人伦偷精品视频| 亚洲精品国产精品久久久不卡| 草草在线视频免费看| 精品乱码久久久久久99久播| 丰满人妻一区二区三区视频av | 人妻丰满熟妇av一区二区三区| 一级毛片精品| 国产亚洲精品一区二区www| 亚洲 国产 在线| 我的老师免费观看完整版| www日本黄色视频网| 国产精品爽爽va在线观看网站| 国产视频一区二区在线看| 欧美日本视频| 日本五十路高清| 国产精品久久久人人做人人爽| 淫秽高清视频在线观看| 欧美一区二区国产精品久久精品 | 三级男女做爰猛烈吃奶摸视频| 久久国产精品影院| 啦啦啦韩国在线观看视频| 久久人人精品亚洲av| 亚洲成人中文字幕在线播放| 无遮挡黄片免费观看| 欧美黄色片欧美黄色片| 好男人在线观看高清免费视频| 国产69精品久久久久777片 | 我的老师免费观看完整版| 国产不卡一卡二| 叶爱在线成人免费视频播放| 国产真实乱freesex| 日本撒尿小便嘘嘘汇集6| 国产av麻豆久久久久久久| 亚洲一区中文字幕在线| 在线观看66精品国产| 亚洲国产高清在线一区二区三| 亚洲欧美日韩无卡精品| 19禁男女啪啪无遮挡网站| 一本久久中文字幕| 国产v大片淫在线免费观看| 小说图片视频综合网站| 亚洲精品久久成人aⅴ小说| 日韩欧美一区二区三区在线观看| 一个人免费在线观看的高清视频| 啪啪无遮挡十八禁网站| 制服诱惑二区| 九色国产91popny在线| 老司机深夜福利视频在线观看| 一级黄色大片毛片| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲精品美女久久av网站| 香蕉丝袜av| 一级作爱视频免费观看| 无人区码免费观看不卡| 国产v大片淫在线免费观看| 窝窝影院91人妻| 精品少妇一区二区三区视频日本电影| 久久欧美精品欧美久久欧美| or卡值多少钱| 麻豆成人午夜福利视频| 欧美另类亚洲清纯唯美| 男男h啪啪无遮挡| 国产一级毛片七仙女欲春2| 久久久国产精品麻豆| 在线观看免费日韩欧美大片| 亚洲精品一卡2卡三卡4卡5卡| 国产av不卡久久| 亚洲精品久久国产高清桃花| 国产成+人综合+亚洲专区| 国产av一区二区精品久久| 日本免费a在线| 9191精品国产免费久久| 91麻豆精品激情在线观看国产| 婷婷丁香在线五月| 欧美日韩精品网址| 国产精品 国内视频| 久久久久性生活片| 久久欧美精品欧美久久欧美| 黄频高清免费视频| 国产激情欧美一区二区| 亚洲国产欧洲综合997久久,| 亚洲第一电影网av| 99久久精品热视频| 女警被强在线播放| 麻豆国产av国片精品| 全区人妻精品视频| а√天堂www在线а√下载| 久久久久久久久免费视频了| 天天一区二区日本电影三级| av中文乱码字幕在线| 1024手机看黄色片| 成人特级黄色片久久久久久久| 国产精品精品国产色婷婷| 熟女少妇亚洲综合色aaa.| 国产精品电影一区二区三区| av欧美777| 99久久无色码亚洲精品果冻| 国产精品电影一区二区三区| 老司机在亚洲福利影院| 亚洲国产欧洲综合997久久,| www国产在线视频色| 国产精品久久久久久亚洲av鲁大| 俄罗斯特黄特色一大片| 亚洲午夜理论影院| 欧美黑人欧美精品刺激| 久久午夜综合久久蜜桃| 伊人久久大香线蕉亚洲五| 日本一区二区免费在线视频| 国产精品免费视频内射| 久久久久九九精品影院| 国产精品国产高清国产av| 搡老妇女老女人老熟妇| 一级毛片高清免费大全| 国产精品爽爽va在线观看网站| 日本三级黄在线观看| 国产野战对白在线观看| 18美女黄网站色大片免费观看| 精品第一国产精品| 国产欧美日韩精品亚洲av| 给我免费播放毛片高清在线观看| 韩国av一区二区三区四区| 久久精品夜夜夜夜夜久久蜜豆 | 韩国av一区二区三区四区| 深夜精品福利| 日本一二三区视频观看| 制服人妻中文乱码| 大型黄色视频在线免费观看| 日韩欧美 国产精品| АⅤ资源中文在线天堂| 真人一进一出gif抽搐免费| 精品欧美一区二区三区在线| 亚洲中文日韩欧美视频| 老司机在亚洲福利影院| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品自产拍在线观看55亚洲| 在线播放国产精品三级| 性欧美人与动物交配| 老司机靠b影院| 国产高清视频在线观看网站| 在线观看午夜福利视频| 成人三级做爰电影| 国产视频一区二区在线看| 免费电影在线观看免费观看| 日本一二三区视频观看| 亚洲第一电影网av| 日韩免费av在线播放| 欧美日韩乱码在线| 国产探花在线观看一区二区| 日韩 欧美 亚洲 中文字幕| 国内精品久久久久久久电影| 不卡一级毛片| 亚洲中文av在线| 亚洲av中文字字幕乱码综合| 国产欧美日韩精品亚洲av| 国产v大片淫在线免费观看| 男女下面进入的视频免费午夜| 三级毛片av免费| 变态另类成人亚洲欧美熟女| 色综合亚洲欧美另类图片| 婷婷丁香在线五月| 国语自产精品视频在线第100页| АⅤ资源中文在线天堂| 日韩欧美三级三区| 午夜免费观看网址| 亚洲狠狠婷婷综合久久图片| 久久热在线av| 国产亚洲精品av在线| 亚洲男人天堂网一区| 日本一二三区视频观看| 成年版毛片免费区| 亚洲全国av大片| cao死你这个sao货| 国产精品乱码一区二三区的特点| 久久久久亚洲av毛片大全| 床上黄色一级片| 久久香蕉激情| 国产激情偷乱视频一区二区| 久久99热这里只有精品18| 色精品久久人妻99蜜桃| 久久草成人影院| 特级一级黄色大片| 国产一区二区激情短视频| 1024手机看黄色片| 国产亚洲精品av在线| 男女下面进入的视频免费午夜| a在线观看视频网站| 亚洲av中文字字幕乱码综合| 成人三级黄色视频| 国产精华一区二区三区| 国产亚洲精品一区二区www| 国产亚洲欧美在线一区二区| 亚洲在线自拍视频| 美女大奶头视频| 成年女人毛片免费观看观看9| 欧美又色又爽又黄视频| a级毛片在线看网站| 欧美精品啪啪一区二区三区| 精品久久久久久久毛片微露脸| 精品久久久久久久久久免费视频| 亚洲精品美女久久久久99蜜臀| 亚洲国产中文字幕在线视频| 啦啦啦免费观看视频1| 一进一出抽搐动态| ponron亚洲| 午夜精品在线福利| 特级一级黄色大片| 在线观看午夜福利视频| 婷婷精品国产亚洲av| 亚洲七黄色美女视频| 亚洲av成人精品一区久久| 精品国产美女av久久久久小说| 老司机福利观看| 久久久久久免费高清国产稀缺| 女警被强在线播放| 欧美一区二区精品小视频在线| 俺也久久电影网| 嫩草影视91久久| 亚洲欧美日韩高清在线视频| www.自偷自拍.com| 黄片大片在线免费观看| 在线观看www视频免费| 变态另类丝袜制服| 欧美大码av| 国产av一区二区精品久久| aaaaa片日本免费| 久久精品国产综合久久久| 亚洲成av人片在线播放无| 久久久水蜜桃国产精品网| 又黄又粗又硬又大视频| 国产精品 国内视频| 亚洲av五月六月丁香网| 亚洲欧美日韩高清专用| 国产精品国产高清国产av| a在线观看视频网站| 在线观看免费视频日本深夜| a级毛片a级免费在线| 欧美绝顶高潮抽搐喷水| 亚洲国产欧美人成| 亚洲av电影不卡..在线观看| 波多野结衣高清作品| 熟女少妇亚洲综合色aaa.| 999久久久国产精品视频| 国产亚洲欧美98| 日本免费一区二区三区高清不卡| 可以在线观看的亚洲视频| 国内少妇人妻偷人精品xxx网站 | 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久亚洲av毛片大全| 亚洲成人久久爱视频| 最近最新中文字幕大全电影3| 久久久久久国产a免费观看| 午夜精品一区二区三区免费看| 男女做爰动态图高潮gif福利片| 黄色丝袜av网址大全| 久久久久性生活片| 天天躁夜夜躁狠狠躁躁| 色播亚洲综合网| 亚洲欧美精品综合一区二区三区| 一级a爱片免费观看的视频| 成人国语在线视频| 亚洲全国av大片| 18禁裸乳无遮挡免费网站照片| 99国产极品粉嫩在线观看| 亚洲黑人精品在线| 亚洲中文字幕日韩| 两性夫妻黄色片| 天堂av国产一区二区熟女人妻 | 国产亚洲精品第一综合不卡| 真人一进一出gif抽搐免费| 婷婷亚洲欧美| 五月玫瑰六月丁香| 色精品久久人妻99蜜桃| 国产精品亚洲一级av第二区| 久久精品亚洲精品国产色婷小说| 岛国视频午夜一区免费看| 午夜激情福利司机影院| 一区二区三区激情视频| 一级作爱视频免费观看| 国产高清videossex| 免费电影在线观看免费观看| xxx96com| 色噜噜av男人的天堂激情| 在线观看免费午夜福利视频| 99热这里只有精品一区 | 国产欧美日韩一区二区三| 日本在线视频免费播放| 午夜免费成人在线视频| 日韩av在线大香蕉| 国产亚洲精品综合一区在线观看 | 岛国视频午夜一区免费看| 淫妇啪啪啪对白视频| 国产av不卡久久| 蜜桃久久精品国产亚洲av| 巨乳人妻的诱惑在线观看| 窝窝影院91人妻| 啪啪无遮挡十八禁网站| 成人一区二区视频在线观看| 一级片免费观看大全| 又黄又粗又硬又大视频| 久久久国产成人免费| 天天躁夜夜躁狠狠躁躁| 国产精品亚洲av一区麻豆| 国产一区在线观看成人免费| 精品久久蜜臀av无| 日本免费一区二区三区高清不卡| 国产日本99.免费观看| 少妇粗大呻吟视频| 亚洲人成伊人成综合网2020| 1024香蕉在线观看| 亚洲人成77777在线视频| 国产视频一区二区在线看| 俺也久久电影网| 成人高潮视频无遮挡免费网站| 国产伦人伦偷精品视频| 老熟妇仑乱视频hdxx| 性欧美人与动物交配| 一a级毛片在线观看| 黄色毛片三级朝国网站| 亚洲18禁久久av| 一级毛片高清免费大全| 91av网站免费观看| 精品日产1卡2卡| 国产在线观看jvid| 天堂av国产一区二区熟女人妻 | 中文字幕人成人乱码亚洲影| 制服人妻中文乱码| 好男人电影高清在线观看| 亚洲欧美日韩高清专用| 国产三级在线视频| 国产片内射在线| 久久精品影院6| 国产私拍福利视频在线观看| 怎么达到女性高潮| 久久久国产欧美日韩av| 国产片内射在线| 草草在线视频免费看| 窝窝影院91人妻| 一区二区三区激情视频| 一边摸一边做爽爽视频免费| 久久亚洲真实| 国产精品九九99| 丁香六月欧美| 亚洲片人在线观看| 制服诱惑二区| 国产av麻豆久久久久久久| 99久久99久久久精品蜜桃| 久久精品国产亚洲av香蕉五月| 午夜久久久久精精品| 精品电影一区二区在线| 精品一区二区三区四区五区乱码| 日韩欧美 国产精品|