• <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é)的初步比較
    成年动漫av网址| www日本在线高清视频| 人妻人人澡人人爽人人| 欧美日韩亚洲高清精品| 人人妻人人爽人人添夜夜欢视频| 免费人妻精品一区二区三区视频| av天堂久久9| 老司机深夜福利视频在线观看 | 国产成人一区二区三区免费视频网站| 在线观看免费高清a一片| 午夜两性在线视频| 欧美黄色片欧美黄色片| 亚洲天堂av无毛| 高清在线国产一区| 国产精品一二三区在线看| 精品一区二区三卡| 成年美女黄网站色视频大全免费| 亚洲精品乱久久久久久| 色老头精品视频在线观看| 在线观看人妻少妇| 成年美女黄网站色视频大全免费| 免费久久久久久久精品成人欧美视频| 97人妻天天添夜夜摸| 男女无遮挡免费网站观看| 欧美xxⅹ黑人| 高清av免费在线| av不卡在线播放| 精品亚洲乱码少妇综合久久| 国产有黄有色有爽视频| 精品久久蜜臀av无| 久久热在线av| 国产激情久久老熟女| 建设人人有责人人尽责人人享有的| 大片免费播放器 马上看| 国产精品久久久久成人av| 麻豆乱淫一区二区| 亚洲欧美精品综合一区二区三区| 飞空精品影院首页| 国产高清视频在线播放一区 | 91麻豆精品激情在线观看国产 | 久久精品国产亚洲av香蕉五月 | 9191精品国产免费久久| 国产精品影院久久| 亚洲九九香蕉| 十分钟在线观看高清视频www| 国产亚洲欧美精品永久| 日日爽夜夜爽网站| 国产高清视频在线播放一区 | 美女扒开内裤让男人捅视频| 一二三四在线观看免费中文在| 国产淫语在线视频| 十八禁人妻一区二区| 中文字幕最新亚洲高清| 亚洲精品第二区| 免费黄频网站在线观看国产| 国产亚洲av片在线观看秒播厂| 视频在线观看一区二区三区| av福利片在线| 色婷婷av一区二区三区视频| 正在播放国产对白刺激| 淫妇啪啪啪对白视频 | 久久久久久久国产电影| 少妇被粗大的猛进出69影院| 日韩熟女老妇一区二区性免费视频| 国产av精品麻豆| 亚洲一码二码三码区别大吗| 18在线观看网站| 精品国产乱码久久久久久小说| 正在播放国产对白刺激| 精品国产国语对白av| 亚洲熟女毛片儿| tube8黄色片| 国产欧美日韩精品亚洲av| 国产有黄有色有爽视频| 最近最新免费中文字幕在线| 男人操女人黄网站| 日本vs欧美在线观看视频| 多毛熟女@视频| 99久久精品国产亚洲精品| 在线观看免费高清a一片| 国产精品99久久99久久久不卡| 午夜精品国产一区二区电影| h视频一区二区三区| 国产三级黄色录像| 最近最新中文字幕大全免费视频| 麻豆国产av国片精品| 亚洲精品国产一区二区精华液| 久久午夜综合久久蜜桃| 美女福利国产在线| 久久精品亚洲熟妇少妇任你| 一级a爱视频在线免费观看| 国产野战对白在线观看| 久久av网站| 91成年电影在线观看| 男女免费视频国产| 午夜精品久久久久久毛片777| 日韩中文字幕视频在线看片| 欧美日韩亚洲国产一区二区在线观看 | 叶爱在线成人免费视频播放| 亚洲av欧美aⅴ国产| 69av精品久久久久久 | 日韩有码中文字幕| 最新在线观看一区二区三区| 欧美一级毛片孕妇| 免费在线观看黄色视频的| 欧美少妇被猛烈插入视频| 国产在视频线精品| 日韩 亚洲 欧美在线| 麻豆乱淫一区二区| avwww免费| 女人久久www免费人成看片| 国产伦人伦偷精品视频| 一级,二级,三级黄色视频| 十分钟在线观看高清视频www| 国产高清国产精品国产三级| 中文字幕最新亚洲高清| av在线播放精品| 国产精品99久久99久久久不卡| 亚洲伊人色综图| 欧美黄色片欧美黄色片| 人成视频在线观看免费观看| 久久这里只有精品19| 丰满少妇做爰视频| 制服人妻中文乱码| 99久久99久久久精品蜜桃| 国产精品国产三级国产专区5o| 天天躁夜夜躁狠狠躁躁| 日韩精品免费视频一区二区三区| 免费高清在线观看日韩| 交换朋友夫妻互换小说| 最黄视频免费看| 爱豆传媒免费全集在线观看| 啦啦啦在线免费观看视频4| 国产精品1区2区在线观看. | 久久精品国产亚洲av高清一级| 纵有疾风起免费观看全集完整版| 丰满饥渴人妻一区二区三| 亚洲成人免费av在线播放| 少妇被粗大的猛进出69影院| 一个人免费在线观看的高清视频 | 中文字幕av电影在线播放| 亚洲 国产 在线| 亚洲精品国产精品久久久不卡| 一进一出抽搐动态| 在线观看舔阴道视频| 欧美精品高潮呻吟av久久| 交换朋友夫妻互换小说| 久久ye,这里只有精品| 妹子高潮喷水视频| 久久久久精品人妻al黑| 女警被强在线播放| 美国免费a级毛片| 午夜影院在线不卡| 国产精品一区二区在线观看99| 纵有疾风起免费观看全集完整版| 久久性视频一级片| 一级片'在线观看视频| 午夜91福利影院| 999久久久精品免费观看国产| videos熟女内射| 美女脱内裤让男人舔精品视频| www.精华液| 999久久久精品免费观看国产| 精品久久久久久久毛片微露脸 | 亚洲少妇的诱惑av| 大香蕉久久成人网| 黑人操中国人逼视频| 日韩欧美一区视频在线观看| 欧美人与性动交α欧美精品济南到| 王馨瑶露胸无遮挡在线观看| 国产一卡二卡三卡精品| 久久久久视频综合| 欧美日韩亚洲高清精品| 国产福利在线免费观看视频| 窝窝影院91人妻| 精品亚洲乱码少妇综合久久| 两性夫妻黄色片| 91九色精品人成在线观看| 一区福利在线观看| 男人操女人黄网站| 一本综合久久免费| 亚洲久久久国产精品| av片东京热男人的天堂| 性色av一级| 最近中文字幕2019免费版| 黑人猛操日本美女一级片| 性少妇av在线| 黄色毛片三级朝国网站| 一区二区日韩欧美中文字幕| 桃花免费在线播放| 妹子高潮喷水视频| 国产精品 国内视频| 久久久久久久精品精品| 考比视频在线观看| 日本撒尿小便嘘嘘汇集6| 久热爱精品视频在线9| 女性生殖器流出的白浆| 欧美大码av| 十八禁高潮呻吟视频| 免费看十八禁软件| 欧美国产精品va在线观看不卡| a级片在线免费高清观看视频| 人人妻人人澡人人看| 老司机影院成人| 一二三四社区在线视频社区8| 热99re8久久精品国产| 久久久久网色| 国产精品.久久久| 欧美另类亚洲清纯唯美| 成年av动漫网址| 69av精品久久久久久 | 男人舔女人的私密视频| 免费高清在线观看视频在线观看| 成年人午夜在线观看视频| 99国产精品免费福利视频| 各种免费的搞黄视频| 激情视频va一区二区三区| 黑丝袜美女国产一区| 国产片内射在线| 99国产极品粉嫩在线观看| 午夜福利在线观看吧| 亚洲 欧美一区二区三区| 国产真人三级小视频在线观看| 国产一区二区三区av在线| 不卡av一区二区三区| 精品少妇黑人巨大在线播放| 多毛熟女@视频| 乱人伦中国视频| 午夜福利在线观看吧| 99九九在线精品视频| 国产精品国产av在线观看| 三级毛片av免费| 亚洲专区国产一区二区| 国产免费福利视频在线观看| 国产免费一区二区三区四区乱码| 欧美日本中文国产一区发布| 中国国产av一级| 国产精品麻豆人妻色哟哟久久| 两性夫妻黄色片| 五月天丁香电影| 在线天堂中文资源库| 夫妻午夜视频| 国产精品久久久av美女十八| 一区二区av电影网| 蜜桃国产av成人99| 嫁个100分男人电影在线观看| 一区二区三区乱码不卡18| 在线精品无人区一区二区三| 国产高清国产精品国产三级| 在线观看免费视频网站a站| 18禁黄网站禁片午夜丰满| 亚洲精品成人av观看孕妇| 亚洲九九香蕉| 久久 成人 亚洲| 国产一卡二卡三卡精品| 法律面前人人平等表现在哪些方面 | 成年美女黄网站色视频大全免费| 国产成人免费无遮挡视频| 午夜福利一区二区在线看| 少妇人妻久久综合中文| 久久久久久人人人人人| 国产精品 国内视频| 一区二区三区乱码不卡18| 人妻一区二区av| 大片电影免费在线观看免费| 欧美精品高潮呻吟av久久| 99热国产这里只有精品6| 久久国产精品影院| 欧美xxⅹ黑人| 欧美另类一区| 精品久久久精品久久久| 午夜激情久久久久久久| 久久 成人 亚洲| 美女视频免费永久观看网站| 91精品三级在线观看| 别揉我奶头~嗯~啊~动态视频 | 超碰97精品在线观看| 岛国毛片在线播放| 久久久久久亚洲精品国产蜜桃av| 国产精品久久久久成人av| 久久精品aⅴ一区二区三区四区| 一级片免费观看大全| 人人澡人人妻人| 青春草视频在线免费观看| 精品国产一区二区久久| 美女视频免费永久观看网站| 国产色视频综合| 久久精品亚洲熟妇少妇任你| 黄色视频,在线免费观看| 精品少妇久久久久久888优播| 日韩欧美一区二区三区在线观看 | 黄频高清免费视频| 丰满少妇做爰视频| 女人久久www免费人成看片| 成人影院久久| 一区福利在线观看| 熟女少妇亚洲综合色aaa.| 国产精品影院久久| 日韩制服骚丝袜av| 国产精品影院久久| 亚洲精品中文字幕一二三四区 | 欧美一级毛片孕妇| 色播在线永久视频| 亚洲专区国产一区二区| 国产不卡av网站在线观看| 69av精品久久久久久 | 亚洲欧美清纯卡通| 他把我摸到了高潮在线观看 | 女警被强在线播放| 日日摸夜夜添夜夜添小说| 精品人妻一区二区三区麻豆| 亚洲人成电影免费在线| 成人av一区二区三区在线看 | 黑人巨大精品欧美一区二区蜜桃| 亚洲人成电影免费在线| 婷婷色av中文字幕| 欧美日韩成人在线一区二区| 黄网站色视频无遮挡免费观看| 大型av网站在线播放| 黑人操中国人逼视频| 女人久久www免费人成看片| 黄色怎么调成土黄色| 在线永久观看黄色视频| 欧美另类一区| 久久天堂一区二区三区四区| 精品一区在线观看国产| 国产一卡二卡三卡精品| 欧美大码av| 少妇被粗大的猛进出69影院| 成人影院久久| 亚洲国产精品999| 丰满人妻熟妇乱又伦精品不卡| 欧美在线黄色| 欧美日韩亚洲国产一区二区在线观看 | 亚洲综合色网址| tube8黄色片| www.999成人在线观看| 夫妻午夜视频| 蜜桃在线观看..| 免费久久久久久久精品成人欧美视频| 国产又爽黄色视频| 久久久久久亚洲精品国产蜜桃av| 免费黄频网站在线观看国产| 亚洲第一欧美日韩一区二区三区 | 国产欧美日韩一区二区三 | 丁香六月天网| 中文字幕色久视频| 亚洲国产看品久久| 欧美中文综合在线视频| 亚洲国产看品久久| 黑人猛操日本美女一级片| 日韩精品免费视频一区二区三区| 多毛熟女@视频| 一区福利在线观看| 伦理电影免费视频| 免费在线观看完整版高清| 成人影院久久| 中文字幕人妻丝袜一区二区| 青青草视频在线视频观看| 久久精品aⅴ一区二区三区四区| 高清av免费在线| 国产淫语在线视频| 法律面前人人平等表现在哪些方面 | 国产欧美日韩综合在线一区二区| 男人舔女人的私密视频| 日韩一卡2卡3卡4卡2021年| 91字幕亚洲| 18禁黄网站禁片午夜丰满| 久久久久久久国产电影| 国产精品秋霞免费鲁丝片| 久久精品国产a三级三级三级| 午夜老司机福利片| 亚洲成人免费av在线播放| 成人国语在线视频| 丁香六月天网| 各种免费的搞黄视频| av在线老鸭窝| 国产男女内射视频| 国产免费av片在线观看野外av| 黑人猛操日本美女一级片| 另类亚洲欧美激情| 精品第一国产精品| 亚洲精品粉嫩美女一区| 亚洲中文日韩欧美视频| 免费观看人在逋| 亚洲va日本ⅴa欧美va伊人久久 | 一本综合久久免费| 亚洲一码二码三码区别大吗| 一级黄色大片毛片| tocl精华| 亚洲全国av大片| 视频区图区小说| av天堂在线播放| 亚洲av男天堂| 免费少妇av软件| 老司机福利观看| 美女中出高潮动态图| 国产激情久久老熟女| 成年动漫av网址| 欧美另类一区| 日本精品一区二区三区蜜桃| 啦啦啦视频在线资源免费观看| 亚洲av日韩精品久久久久久密| 国产免费一区二区三区四区乱码| 欧美激情高清一区二区三区| 女性被躁到高潮视频| 国产成人精品久久二区二区91| 精品福利观看| 久久九九热精品免费| 亚洲精品中文字幕在线视频| 午夜老司机福利片| 国产一区有黄有色的免费视频| 大片电影免费在线观看免费| a级片在线免费高清观看视频| 黄色a级毛片大全视频| 久久久久精品人妻al黑| 精品国产国语对白av| 免费观看av网站的网址| 欧美日韩精品网址| 亚洲美女黄色视频免费看| 老熟妇仑乱视频hdxx| 国产成人av教育| 日本91视频免费播放| 国产亚洲欧美精品永久| 涩涩av久久男人的天堂| 人人妻人人澡人人看| 制服人妻中文乱码| xxxhd国产人妻xxx| 国产真人三级小视频在线观看| tocl精华| 亚洲欧美日韩另类电影网站| av超薄肉色丝袜交足视频| 色老头精品视频在线观看| 永久免费av网站大全| 色老头精品视频在线观看| 免费观看人在逋| 高清欧美精品videossex| 少妇裸体淫交视频免费看高清 | 韩国高清视频一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 一级毛片女人18水好多| 日韩欧美一区视频在线观看| e午夜精品久久久久久久| 老司机靠b影院| a在线观看视频网站| 九色亚洲精品在线播放| 免费一级毛片在线播放高清视频 | 蜜桃在线观看..| 欧美老熟妇乱子伦牲交| 2018国产大陆天天弄谢| 久久狼人影院| 精品久久蜜臀av无| 极品少妇高潮喷水抽搐| 视频区图区小说| 亚洲av成人不卡在线观看播放网 | av免费在线观看网站| 日韩熟女老妇一区二区性免费视频| 久久久久精品国产欧美久久久 | 女人久久www免费人成看片| 黑人欧美特级aaaaaa片| 久久精品熟女亚洲av麻豆精品| e午夜精品久久久久久久| 在线亚洲精品国产二区图片欧美| 成人手机av| 亚洲专区国产一区二区| 国产福利在线免费观看视频| 一区二区三区乱码不卡18| 五月开心婷婷网| 50天的宝宝边吃奶边哭怎么回事| av视频免费观看在线观看| 色综合欧美亚洲国产小说| 在线观看人妻少妇| 夜夜夜夜夜久久久久| 国产精品影院久久| 国产国语露脸激情在线看| 飞空精品影院首页| 亚洲色图 男人天堂 中文字幕| 又紧又爽又黄一区二区| 老鸭窝网址在线观看| 欧美日韩av久久| 亚洲国产欧美一区二区综合| www日本在线高清视频| 免费观看人在逋| 少妇 在线观看| 亚洲人成电影免费在线| 国产真人三级小视频在线观看| 亚洲国产精品一区三区| 亚洲性夜色夜夜综合| 九色亚洲精品在线播放| 亚洲人成电影免费在线| 亚洲国产看品久久| 久久99热这里只频精品6学生| av有码第一页| 中文字幕精品免费在线观看视频| 亚洲人成电影观看| 亚洲av成人不卡在线观看播放网 | 国产主播在线观看一区二区| 欧美+亚洲+日韩+国产| 菩萨蛮人人尽说江南好唐韦庄| 新久久久久国产一级毛片| 夜夜夜夜夜久久久久| 亚洲成人免费av在线播放| 久久久久久亚洲精品国产蜜桃av| 涩涩av久久男人的天堂| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久久国产电影| 国产极品粉嫩免费观看在线| 12—13女人毛片做爰片一| 久久久国产成人免费| 悠悠久久av| 久久香蕉激情| 免费在线观看影片大全网站| 别揉我奶头~嗯~啊~动态视频 | 午夜精品久久久久久毛片777| 成人av一区二区三区在线看 | 三上悠亚av全集在线观看| 女警被强在线播放| 乱人伦中国视频| e午夜精品久久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 成人国产av品久久久| 另类精品久久| 后天国语完整版免费观看| 国产高清视频在线播放一区 | 嫩草影视91久久| 亚洲综合色网址| 捣出白浆h1v1| 亚洲午夜精品一区,二区,三区| 欧美人与性动交α欧美精品济南到| 久久久国产欧美日韩av| 久久国产精品人妻蜜桃| 久久中文字幕一级| 亚洲欧美激情在线| 成人国语在线视频| 中文字幕人妻丝袜一区二区| 欧美国产精品va在线观看不卡| 男人操女人黄网站| 青春草视频在线免费观看| av国产精品久久久久影院| 国产成人精品无人区| 成人av一区二区三区在线看 | 精品人妻熟女毛片av久久网站| 不卡av一区二区三区| 成人黄色视频免费在线看| 国产野战对白在线观看| 人人澡人人妻人| 欧美xxⅹ黑人| 久久久久久人人人人人| √禁漫天堂资源中文www| 午夜福利免费观看在线| 侵犯人妻中文字幕一二三四区| 日韩一卡2卡3卡4卡2021年| 亚洲精品一二三| 在线观看舔阴道视频| 国产黄频视频在线观看| 国产日韩欧美在线精品| 窝窝影院91人妻| av网站免费在线观看视频| 91字幕亚洲| 亚洲欧洲精品一区二区精品久久久| 亚洲性夜色夜夜综合| 动漫黄色视频在线观看| 亚洲综合色网址| 国产三级黄色录像| 又紧又爽又黄一区二区| 欧美黄色淫秽网站| 男女床上黄色一级片免费看| 女性被躁到高潮视频| 日本一区二区免费在线视频| 久久久欧美国产精品| 在线观看舔阴道视频| 欧美精品一区二区免费开放| 国产精品久久久av美女十八| 日日爽夜夜爽网站| 亚洲精品自拍成人| 亚洲少妇的诱惑av| 精品国产乱子伦一区二区三区 | 亚洲国产欧美日韩在线播放| 欧美另类亚洲清纯唯美| 人人妻人人澡人人爽人人夜夜| 黑人巨大精品欧美一区二区mp4| 日韩大片免费观看网站| 成年人免费黄色播放视频| 下体分泌物呈黄色| 99久久精品国产亚洲精品| 久久久欧美国产精品| 亚洲自偷自拍图片 自拍| av福利片在线| 久久天堂一区二区三区四区| 欧美 亚洲 国产 日韩一| 69精品国产乱码久久久| 一区福利在线观看| 黑人操中国人逼视频| 欧美成狂野欧美在线观看| 成年动漫av网址| 一区二区三区四区激情视频| 国产有黄有色有爽视频| 亚洲欧美精品综合一区二区三区| 亚洲av欧美aⅴ国产| 丰满饥渴人妻一区二区三| 99精品久久久久人妻精品| 免费在线观看影片大全网站| 精品国产乱码久久久久久小说| www.自偷自拍.com| 十分钟在线观看高清视频www| 午夜福利乱码中文字幕| 久久热在线av| 日韩熟女老妇一区二区性免费视频| 我要看黄色一级片免费的| 色精品久久人妻99蜜桃| 欧美日韩黄片免| 国产一区二区在线观看av| 日韩电影二区| 国产精品欧美亚洲77777| 在线精品无人区一区二区三|