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

    南海海面風(fēng)場、高度場和溫度場近十年的變化趨勢及年際變化特征

    2017-06-22 14:36:14王佳瑩方國洪王永剛
    海洋科學(xué)進(jìn)展 2017年2期
    關(guān)鍵詞:海表年際風(fēng)場

    王佳瑩,方國洪,2,3*,王永剛,2,3

    (1.國家海洋局第一海洋研究所,山東青島266061; 2.海洋環(huán)境科學(xué)和數(shù)值模擬國家海洋局重點(diǎn)實(shí)驗(yàn)室,山東青島266061; 3.青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室區(qū)域海洋動(dòng)力學(xué)和數(shù)值模擬功能實(shí)驗(yàn)室,山東青島266237)

    研究論文

    南海海面風(fēng)場、高度場和溫度場近十年的變化趨勢及年際變化特征

    王佳瑩1,方國洪1,2,3*,王永剛1,2,3

    (1.國家海洋局第一海洋研究所,山東青島266061; 2.海洋環(huán)境科學(xué)和數(shù)值模擬國家海洋局重點(diǎn)實(shí)驗(yàn)室,山東青島266061; 3.青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室區(qū)域海洋動(dòng)力學(xué)和數(shù)值模擬功能實(shí)驗(yàn)室,山東青島266237)

    利用2002—2011年南海海面風(fēng)場、高度場和溫度場的月平均資料分析了這些要素的變化趨勢及年際變化特征。計(jì)算得到區(qū)域平均海面風(fēng)東分量的線性趨勢為(-0.012±0.014)(m·s-1)·a-1,北分量的線性趨勢為(-0.014±0.019)(m·s-1)·a-1,其年際變化與NIN~O3指數(shù)相關(guān)系數(shù)可達(dá)0.63,滯后1個(gè)月;區(qū)域平均海面風(fēng)應(yīng)力旋度的線性趨勢為(0.099±0.330)×10-9(N·m-3)·a-1,其年際變化與中太平洋厄爾尼諾(EMI)指數(shù)呈高度相關(guān),達(dá)到-0.85,滯后3個(gè)月;區(qū)域平均海面高度的線性趨勢為(0.665±0.200)cm·a-1,其年際變化與EMI指數(shù)相關(guān)系數(shù)可達(dá)-0.80,滯后2個(gè)月;區(qū)域平均海表溫度的線性趨勢為(-0.016±0.017)℃·a-1,其年際變化與NIN~O3指數(shù)相關(guān)系數(shù)可達(dá)0.68,滯后4個(gè)月。其中海表溫度的負(fù)趨勢與近期全球變暖的停滯相一致。采用經(jīng)驗(yàn)正交函數(shù)分解方法分析了各要素的年際變化特征,結(jié)果顯示海面風(fēng)場的第一模態(tài)呈現(xiàn)海盆尺度的反氣旋格局,其對(duì)應(yīng)的時(shí)間序列與NIN~O3.4指數(shù)相關(guān)系數(shù)最大可達(dá)0.68,滯后5個(gè)月。海面高度場的第一模態(tài)沿南海東邊界和越南東部出現(xiàn)較高值,其對(duì)應(yīng)的時(shí)間序列的變化滯后于EMI指數(shù)3個(gè)月,呈負(fù)相關(guān),但系數(shù)最大為-0.32。海面溫度場的第一模態(tài)呈西高東低的分布狀態(tài),其對(duì)應(yīng)的時(shí)間序列與NIN~O3指數(shù)的相關(guān)系數(shù)最大可達(dá)0.62,滯后4個(gè)月。

    南海;風(fēng)場;高度場;溫度場;變化趨勢;年際變化

    南海是太平洋最大的一個(gè)邊緣海,水域面積約為365萬km2,平均水深約為1 200 m,其中最深處位于南海中部偏北,約為5 000 m。南海北靠中國大陸和臺(tái)灣島,其間有臺(tái)灣海峽與東海相連;南臨加里曼丹和蘇門答臘島,經(jīng)卡里馬塔海峽與爪哇海相連;東臨菲律賓群島,經(jīng)呂宋海峽和明都洛海峽與太平洋相連;西臨中南半島和馬來半島,其南端有馬六甲海峽與安達(dá)曼海相通(圖1)。我們?nèi)?8°~124°E,2°S~24°N的區(qū)域作為研究范圍。

    南海的地理位置決定了其環(huán)境受東亞季風(fēng)的影響而表現(xiàn)出強(qiáng)烈的季節(jié)性變化。同時(shí),作為熱帶太平洋的邊緣海,南海環(huán)境也展示了十分顯著的年際變化特征。

    對(duì)于南海各要素的研究一直是物理海洋學(xué)家關(guān)注的熱點(diǎn)。Chao等[1]利用美國氣候中心COADS(Comprehensive Ocean-Atmosphere Data Set)風(fēng)應(yīng)力資料分析了1982—1983年ENSO事件與南海海面風(fēng)場異常的相關(guān)性。Wu等[2]利用美國國家環(huán)境預(yù)報(bào)中心(NCEP)的再分析風(fēng)應(yīng)力資料分析了1992—1995年間南海的季節(jié)和年際變化特征。Liang等[3]利用歐洲中期天氣預(yù)報(bào)的再分析結(jié)果研究了1997—1998年厄爾尼諾事件期間南海風(fēng)場的異?,F(xiàn)象。Hwang和Chen[4]利用ERS-1/2衛(wèi)星觀測5 a海面風(fēng)應(yīng)力資料研究了南海從季節(jié)到年際的變化特征。

    (李 燕 編輯)

    在很多前人利用衛(wèi)星高度計(jì)資料的研究中都表明,南海海面高度存在顯著的季節(jié)、年際和長期變化。例如Ho等[5]的研究中提到南海海面高度場對(duì)厄爾尼諾現(xiàn)象的響應(yīng);Li等[6]發(fā)現(xiàn)在1993—1999年間南海的海面高度以平均每年1 cm的速度升高。

    Chu等[7]利用1982—1994年NCEP的海表溫度月平均數(shù)據(jù)研究南海海表溫度的時(shí)空變化特征,給出了海表溫度的季節(jié)分布特征并定義了南海北部強(qiáng)烈的冷暖異常。Klein等[8]發(fā)現(xiàn)海表溫度變化滯后ENSO指數(shù)5個(gè)月發(fā)生。Wang等[9]揭示了1997—1998年南海暖事件及其與厄爾尼諾事件的密切聯(lián)系。Liu等[10]發(fā)現(xiàn)南海冬季冷舌的變化特征與NIN~O3指數(shù)有很好的相關(guān)性。

    以上的研究都涉及南海環(huán)境要素與ENSO的關(guān)系,但是大部分都只涉及一種要素。此外,除了海面高度場,其他要素的整體趨勢都沒被提及。Fang等[11]對(duì)南海海面風(fēng)場、高度場和溫度場的變化趨勢和年際變化同時(shí)進(jìn)行了分析,揭示了這些環(huán)境參量重要的變化特征。由于當(dāng)時(shí)可利用的資料只到2003年,而現(xiàn)在可獲取的海面風(fēng)場觀測資料序列最長可到2011-12,海面高度場觀測資料可到2014-04,海表溫度場資料可到2015年。為統(tǒng)一,我們將對(duì)2002—2011年的這些要素進(jìn)行分析。

    圖1 南海地圖Fig.1 Map of the South China Sea

    1 資料來源與數(shù)據(jù)處理

    我們采用的海面高度資料是法國衛(wèi)星中心(Aviso)的Ssalto/Duacs多任務(wù)高度計(jì)(包括Saral,Cryosat-2,Jason-1/2,T/P,Envisat,GFO,ERS1/2和Geosat)產(chǎn)品,是由周平均延遲海平面距平資料進(jìn)行每月平均計(jì)算得到的月平均資料,空間分辨率為(1/4)°×(1/4)°,當(dāng)前可獲取的時(shí)間序列為1993-01—2014-04;海表溫度資料是美國國家環(huán)境預(yù)報(bào)中心(NCEP)的Reynolds最優(yōu)插值海面溫度產(chǎn)品,是由美國國家海洋和大氣管理局(National Oceanic and Atmospheric Administration,NOAA)的高級(jí)甚高分辨率輻射計(jì)(AVHRR)得到的衛(wèi)星SST數(shù)據(jù)和由現(xiàn)場SST觀測值通過最優(yōu)插值得到的月平均資料,空間分辨率為1°×1°,可獲取的時(shí)間序列為1981-12至今;在Fang等[11]的文章中采用的風(fēng)場資料來自法國ERS-1數(shù)據(jù)處理與歸檔中心(Centre ERS d'Archivage et de Traitement)的CERSAT資料,由于CERSAT的時(shí)間序列較短(1992-02—2007-12),因而在本研究中采用的海面風(fēng)場資料改用CCMP(cross-calibrated multi-platform,多平臺(tái)交互校正)月平均風(fēng)場資料,該產(chǎn)品是通過變分分析方法對(duì)由遙感系統(tǒng)得到的交互校準(zhǔn)衛(wèi)星風(fēng)場資料進(jìn)行處理,得到高分辨率(1/4)°×(1/4)°的分析結(jié)果,其中交互校準(zhǔn)衛(wèi)星風(fēng)場資料來自SSM/I,SSMIS,AMSR-E,TRMM TMI,QuikSCAT,SeaWinds,WindSat和一些其他的衛(wèi)星儀器,可獲取的時(shí)間序列為1988-01—2011-12。此外,我們還利用CCMP風(fēng)場資料計(jì)算得出海面風(fēng)應(yīng)力場。為統(tǒng)一資料時(shí)間長度,本研究選取2002-01—2011-12作為4組數(shù)據(jù)的統(tǒng)一時(shí)間長度。

    我們著重研究這十年變化趨勢和年際變化特征,因此在進(jìn)行資料處理時(shí)對(duì)所有格點(diǎn)數(shù)據(jù)進(jìn)行了12個(gè)月的滑動(dòng)平均,濾掉頻率高于每年一周的信號(hào),特別也去掉了季節(jié)信號(hào)。在對(duì)資料進(jìn)行經(jīng)驗(yàn)正交函數(shù)(Empirical Orthogonal Function,EOF)分析時(shí),我們先對(duì)資料進(jìn)行去線性趨勢的處理,即用原始資料減去平均值和線性趨勢值,以濾掉背景場的影響。為了便于比較,我們還計(jì)算了氣候態(tài)月平均場和年平均場。

    2 背景場分析

    2.1 風(fēng) 場

    南海地處東亞季風(fēng)區(qū),受季風(fēng)影響較大。圖2a為1月氣候態(tài)平均風(fēng)場,由圖可以看出,冬季整個(gè)南海海域都被較強(qiáng)的東北季風(fēng)控制著。圖2b為1月氣候態(tài)平均風(fēng)應(yīng)力旋度場,由圖可知呂宋島西側(cè)和南海南部出現(xiàn)了氣旋性風(fēng)應(yīng)力旋度。圖2c為7月氣候態(tài)平均風(fēng)場,由圖可知夏季南海主要被西南風(fēng)控制,且南海南部風(fēng)速大于北部風(fēng)速,最大風(fēng)速出現(xiàn)在中南半島的東南外海。圖2d為7月氣候態(tài)平均風(fēng)應(yīng)力旋度場,對(duì)應(yīng)圖2b中風(fēng)速最大的位置北部有一個(gè)氣旋性風(fēng)應(yīng)力旋度結(jié)構(gòu),其南部有一個(gè)反氣旋性風(fēng)應(yīng)力旋度結(jié)構(gòu)。這種形式的風(fēng)應(yīng)力場是形成東南越南離岸流的原因。圖2e為氣候態(tài)年平均風(fēng)場,圖中南海主要被偏北風(fēng)控制,北部風(fēng)速大于南部風(fēng)速,說明冬季季風(fēng)要強(qiáng)于夏季季風(fēng),且兩者之差在南海北部更大。圖2f為氣候態(tài)年平均風(fēng)應(yīng)力旋度場,圖中在呂宋島西側(cè)和中南半島東側(cè)各有一個(gè)氣旋性風(fēng)應(yīng)力旋度結(jié)構(gòu),前者大于后者。本文得到的氣候態(tài)平均風(fēng)場圖和風(fēng)應(yīng)力旋度場圖與Fang等[11]利用1993-2003年CERSAT風(fēng)場資料計(jì)算得到的結(jié)果大體一致。Lian等[21]給出了由9種風(fēng)應(yīng)力資料計(jì)算得到的氣候態(tài)平均風(fēng)應(yīng)力旋度場圖,表明高分辨率資料得到的結(jié)果更為合理。我們得到的分布特征與Lian等[21]所得結(jié)果基本一致。

    2.2 高度場

    圖3a為南海1月氣候態(tài)平均海面高度場(SSH),總體上呈西北高、東南低,顯然與冬季東北季風(fēng)有關(guān)。此外,在呂宋島西側(cè)和中南半島東南部各有一個(gè)低SSH區(qū)域,與圖2d中的氣旋性風(fēng)應(yīng)力旋度對(duì)應(yīng),表明該海面高度異常是由風(fēng)應(yīng)力旋度引起的。圖3b為南海7月氣候態(tài)平均海面高度場,總體上呈西北低、東南高,顯然與夏季西南季風(fēng)有關(guān)。此外,在(110°E,10°N)附近還有一個(gè)較弱的高SSH區(qū)域,這個(gè)高SSH區(qū)會(huì)在8,9月進(jìn)一步發(fā)展變強(qiáng),形成一個(gè)反氣旋式環(huán)流。圖3c為海面高度的氣候態(tài)年平均,總體上為正值,數(shù)值在3~4 cm,系南海海面高度呈升高趨勢的結(jié)果。Fang等[11]利用1993—2003年AVISO海面高度資料計(jì)算得到的氣候態(tài)平均海面高度場圖與本文中分布相似。Ho等[5]利用1992-11—1997-10的T/P海面高度資料計(jì)算各月氣候態(tài)平均海面高度場分布,其中1月和7月分布結(jié)構(gòu)與本文的結(jié)果大體一致。與前人結(jié)果的主要差別是,我們得到的平均海面高度要比他們的略高,這是由于海平面在不斷上升而我們采用了更新近的資料。

    圖2 2002-2011年期間南海海面風(fēng)場及南海海面風(fēng)應(yīng)力旋度場Fig.2 Climatological mean fields of the surface winds and the wind stress curl during 2002-2011

    圖3 2002-2011年期間南海海面高度場(cm)Fig.3 Climatological mean fields of the sea surface height anomalies(cm)during 2002-2011

    2.3 溫度場

    圖4a為1月氣候態(tài)平均海表溫度場,可以明顯看出海面溫度呈北低南高分布,且南海北部等溫線更為密集,變化梯度更大。這是由于冬季東北季風(fēng)導(dǎo)致的,東北季風(fēng)對(duì)南海北部的影響有2種作用方式,一是季風(fēng)直接帶來的干冷空氣使海面水溫降低,二是季風(fēng)在北部形成西南向沿岸流,將較冷的沿岸水由東海經(jīng)臺(tái)灣海峽帶入南海。圖4b為海面溫度場7月氣候態(tài)平均,等溫線分布仍是北高南低,只是較冬季分布來得均勻,南北梯度比冬季要小。同時(shí),在越南東南沿岸出現(xiàn)了一個(gè)冷水舌,這個(gè)冷水舌是由上升流引起的。圖4c為海表溫度場氣候態(tài)年平均,等溫線分布與圖4a類似,只是南北梯度相對(duì)變小。Fang等[11]利用1993—2003年NCEP海表溫度資料分析所得的氣候態(tài)平均海表溫度場圖分布與本文所得結(jié)果相近。Xie等[22]利用1985—1999年AVHRR海表溫度資料計(jì)算得到的7月氣候態(tài)平均分布圖與圖4b中分布基本一致。Chu等[7]利用1982—1994年NCEP的海表溫度資料計(jì)算得到的氣候態(tài)年平均分布圖與圖4c也基本一致。與前人結(jié)果的主要差別是,本文得到的平均海表溫度要比他們的略高,這是由于海表溫度在不斷上升而我們采用了更新近的資料。例如與Xie等[22]比較,我們得到的7月氣候態(tài)平均SST(圖4b)在越南東南外海低于29℃的范圍要明顯小于他們的結(jié)果[22](圖3b)。

    圖4 2002-2011年期間南海海表溫度場(℃)Fig.4 Climatological mean fields of the sea surface temperature during 2002-2011(℃)

    3 線性趨勢

    本文利用最小二乘法求取各參量的線性趨勢,最小二乘法是通過最小化偏差的平方和來尋找數(shù)據(jù)的回歸系數(shù),即趨勢值。圖5中曲線為南海海面風(fēng)速東分量和北分量、海面高度和海表溫度區(qū)域平均值的時(shí)間序列,直線代表利用線性回歸分析方法得到的線性趨勢。圖6為南海海面風(fēng)場、高度場和溫度場線性趨勢的空間分布。本文對(duì)所有資料都事先進(jìn)行了12個(gè)月滑動(dòng)平均,因此結(jié)果中不存在季節(jié)信號(hào)。

    表1 2002—2011年南海海面風(fēng)場、海面風(fēng)應(yīng)力旋度場、海面高度場和海表溫度場的線性趨勢值和標(biāo)準(zhǔn)差Table 1 Linear trends and standard deviations of the Surface Wind,Sea Surface Wind Stress Curl,Sea Surface Height, and Sea Surface Temperature of the South China Sea during 2002-2011

    表1給出了各要素線性趨勢及其標(biāo)準(zhǔn)差,其中標(biāo)準(zhǔn)差根據(jù)Fang等[11](Appendix A)提供的方法計(jì)算。圖5a是南海區(qū)域平均的海面風(fēng)場東分量和北分量的時(shí)間序列和線性趨勢圖(2002—2011年),由圖可知,風(fēng)場東分量的區(qū)域平均趨勢為(-0.012±0.014)(m·s-1)·a-1(表1),風(fēng)場北分量的區(qū)域平均趨勢為(-0. 014±0.019)(m·s-1)·a-1。圖6a是海面風(fēng)場線性趨勢的空間分布圖,南海區(qū)域大體為東風(fēng),在南海海盆西部和南海南部分別出現(xiàn)了輻聚和輻散,同時(shí)南海北部風(fēng)速增長率要小于南海南部風(fēng)速增長率,最大風(fēng)速增長率出現(xiàn)在巴拉望島西側(cè)。

    圖5b是南海海域平均海面風(fēng)應(yīng)力旋度的時(shí)間序列和線性趨勢圖(2002—2011年),海面風(fēng)應(yīng)力旋度場的區(qū)域平均趨勢為(0.099±0.33)×10-9(N·m-3)·a-1,呈上升的趨勢。圖6b是海面風(fēng)應(yīng)力旋度場線性趨勢的空間分布圖,由圖可知,南海海域風(fēng)應(yīng)力旋度整體呈上升趨勢分布,在南海海盆中存在一條東北-西南方向的條帶,這里風(fēng)應(yīng)力旋度呈下降趨勢。

    圖5c是南海區(qū)域平均的海面高度場的時(shí)間序列和線性趨勢圖(2002—2011年),海面高度區(qū)域平均的線性趨勢為(0.665±0.200)cm·a-1,呈上升趨勢,與Fang等[11]得到的1993-2003年上升速率基本相同。Church等[12]給出的全球海面高度上升率(1950—2000年)為(0.180±0.030)cm·a-1,馮偉等[13]估計(jì)所得的全球海面高度上升率(2005—2013年)為(0.272±0.070)cm·a-1,這表明近20 a南海海平面上升速率顯著高于全球速率。通過對(duì)比圖5b和5c可知,海面高度場的時(shí)間序列變化與海面風(fēng)應(yīng)力旋度場的時(shí)間序列變化具有一定相關(guān)性,說明海面高度的年際變化受到風(fēng)應(yīng)力旋度變化的影響。圖6c是海面高度線性趨勢的空間分布圖,整體趨勢值在0.1~1.3 cm·a-1范圍內(nèi),其中深海盆區(qū)有相對(duì)較高的上升速率,最大的上升速率位于呂宋島的西側(cè)。

    圖5d是南海區(qū)域平均的海表溫度的時(shí)間序列和線性趨勢圖。2002—2011年期間線性趨勢為(-0.016± 0.017)℃·a-1,呈下降趨勢。鑒于Reynolds最優(yōu)插值海表溫度產(chǎn)品的時(shí)間范圍為1982—2015年,故作為比較,我們還計(jì)算了1982—2015年南海區(qū)域平均海表溫度的線性趨勢,得到的數(shù)值為0.006℃·a-1,比Fang等[11]計(jì)算得到1982—2004年線性趨勢0.026℃·a-1要小。這個(gè)結(jié)果與前面得到的近十年負(fù)增溫率一致。從20世紀(jì)末以來,全球變暖出現(xiàn)了停滯現(xiàn)象[14],南海海表溫度的負(fù)增溫應(yīng)當(dāng)是全球變暖停滯的一個(gè)反映。圖6d是海面溫度場線性趨勢的空間分布圖,整體趨勢值的范圍為-0.05~0.03℃·a-1,南海區(qū)域內(nèi)大部分增溫率為負(fù)值,在南中國陸架和巽他陸架附近負(fù)增溫率最大,而在呂宋海峽附近卻出現(xiàn)正增溫。

    圖5 2002—2011年期間南海各參量區(qū)域平均值的時(shí)間序列和線性趨勢Fig.5 Time series and linear trends of South China Sea regional means during 2002-2011

    圖6 2002—2011年期間南海各參量線性趨勢的空間分布Fig.6 Spatial patterns of the linear trends during 2002-2011

    4 年際變化

    4.1 區(qū)域平均值的年際變化

    為了研究各參量的年際變化,我們將對(duì)參量的時(shí)間序列與代表氣候年際變化的主要指數(shù)進(jìn)行超前/滯后相關(guān)性分析。這些指數(shù)包括NIN~O3,NIN~O3.4,IOD及EMI。其中EMI是Ashok等[15]提出的一個(gè)新的指數(shù),稱作類厄爾尼諾指數(shù)(El Ni?o Modoki Index),是一種不同但類似于傳統(tǒng)的厄爾尼諾現(xiàn)象,這種現(xiàn)象在赤道太平洋中部發(fā)生,也稱中太平洋El Ni?o,發(fā)生時(shí)熱帶太平洋海表溫度呈現(xiàn)一種馬蹄形分布,即沿赤道兩側(cè)出現(xiàn)相對(duì)較低的海表溫度異常(SSTA)夾著赤道太平洋中部較高的SSTA。EMI指數(shù)與NIN~O3指數(shù)基本上不相關(guān),但與NIN~O3.4有一定的相關(guān)性[15]。

    在2002—2011年期間南海各參量區(qū)域平均值時(shí)間序列圖(圖5)中黑色虛線給出了對(duì)應(yīng)期間相關(guān)性最高的氣候年際變化指數(shù)的標(biāo)準(zhǔn)化后時(shí)間序列。圖5a為風(fēng)分量與NIN~O3指數(shù)的時(shí)間序列對(duì)比圖。相關(guān)分析表明,風(fēng)東分量與NIN~O3指數(shù)相關(guān)系數(shù)最大為0.56,超前6個(gè)月;風(fēng)北分量與NIN~O3指數(shù)相關(guān)系數(shù)最大為0.63,滯后1個(gè)月。圖5b為海面風(fēng)應(yīng)力旋度與負(fù)EMI指數(shù)的相近序列對(duì)比圖,海面風(fēng)應(yīng)力旋度與EMI指數(shù)相關(guān)系數(shù)的負(fù)最大值為-0.85,滯后3個(gè)月。圖5c為海面高度與負(fù)EMI指數(shù)的時(shí)間序列對(duì)比圖,海面高度場與EMI指數(shù)相關(guān)系數(shù)的負(fù)最大值為-0.80,滯后2個(gè)月。圖5d為海表溫度與NIN~O3指數(shù)的時(shí)間序列對(duì)比圖,海面溫度場與NIN~O3指數(shù)相關(guān)系數(shù)最大為0.68,滯后4個(gè)月。我們根據(jù)Fang等[11](Appendix A)提供的方法對(duì)以上相關(guān)關(guān)系進(jìn)行顯著性檢驗(yàn),結(jié)果顯示其顯著性水平都低于1%。

    4.2 主成分的年際變化

    4.2.1 主成分分解

    主成分分析(principal component analysis,PCA)也稱特征向量分析或經(jīng)驗(yàn)正交函數(shù)(empirical orthogonal function,EOF)分析,是一種分析矩陣數(shù)據(jù)中的結(jié)構(gòu)特征,提取數(shù)據(jù)主要特征量的方法,最早由Lorenz[16]引入氣象和氣候研究,并在地球科學(xué)及其他學(xué)科中得到了廣泛應(yīng)用。本研究中我們利用EOF方法對(duì)各參量的觀測值進(jìn)行成份分解,將觀測值分解成空間形態(tài)(spatial pattern)與時(shí)間函數(shù)(time coefficient function,TCF)乘積的線性組合,以便進(jìn)行進(jìn)一步的分析。

    為研究各參量的年際變化情況,我們對(duì)原始資料進(jìn)行了12個(gè)月滑動(dòng)平均處理,去掉平均值和線性趨勢值。假定位置為m(m=1,2,…,M),時(shí)間為t(t=1,2,…,N)的參量為x(m,t),則經(jīng)過EOF處理后就可將它分解為一系列的時(shí)間函數(shù)αi(t)和空間函數(shù)Fi(m)的乘積之和,即

    式中,K取(M,N)中較小值;i代表第i個(gè)模態(tài);αi即代表第i個(gè)模態(tài)的時(shí)間函數(shù);Fi代表第i個(gè)模態(tài)的空間形態(tài)。模態(tài)順序是按照x(m,t)的協(xié)方差矩陣特征值λi大小排列的,λi的大小也代表了各模態(tài)的方差。為了方便對(duì)比,本文中將TCF進(jìn)行標(biāo)準(zhǔn)化處理,即滿足

    則相應(yīng)的方差

    在對(duì)具有東分量和北分量的海面風(fēng)場資料進(jìn)行EOF分析時(shí),我們選取了實(shí)向量方法[17](real-vector method),它是多變量EOF分析[18](mutlivariate EOF analysis)的一種特殊情況。

    表2 南海風(fēng)場、風(fēng)應(yīng)力旋度場、高度場及溫度場EOF前3個(gè)模態(tài)的方差、方差貢獻(xiàn)率和標(biāo)準(zhǔn)差Table 2 Variances,contribution percentage and standard deviations of the first three EOFs of the surface wind, wind stress curl,sea surface height,and sea surface temperature of the South China Sea

    4.2.2 風(fēng) 場

    圖7為海面風(fēng)場和海面風(fēng)應(yīng)力旋度場的第一模態(tài)的空間分布和時(shí)間序列,由表2可知,風(fēng)場的第一模態(tài)方差貢獻(xiàn)率為44%,風(fēng)應(yīng)力旋度的第一模態(tài)方差貢獻(xiàn)率為33%。由圖7a可以看出,在南海海盆北部有一個(gè)明顯的反氣旋結(jié)構(gòu),中心位于呂宋島西側(cè)。由圖7b可以看出,對(duì)應(yīng)反氣旋結(jié)構(gòu)的位置存在一個(gè)反氣旋中心,同時(shí)在南海北部大陸沿岸和呂宋海峽以西存在氣旋性中心。圖7c是南海海面風(fēng)場EOF第一模態(tài)的時(shí)間序列與NIN~O3.4指數(shù)的對(duì)比圖。通過對(duì)該時(shí)間序列與NIN~O3.4指數(shù)的超前/滯后相關(guān)分析得出,當(dāng)風(fēng)場第一模態(tài)的時(shí)間序列滯后于NIN~O3.4指數(shù)5個(gè)月時(shí),出現(xiàn)最大的相關(guān)系數(shù)0.68,置信水平超過99%。風(fēng)場第一模態(tài)時(shí)間序列中較大的4個(gè)峰值基本上分別與2002—2003年、2004—2005年、2006—2007年和2009—2010年四次厄爾尼諾事件相對(duì)應(yīng)。Wang等[19]認(rèn)為,在厄爾尼諾事件發(fā)生期間,西北太平洋會(huì)出現(xiàn)一個(gè)反氣旋異常,同時(shí)東亞冬季季風(fēng)會(huì)相對(duì)減弱。南海這一模態(tài)中出現(xiàn)的反氣旋異常實(shí)際上是西北太平洋反氣旋異常的一部分,印證了Wang等[19]的觀點(diǎn)。當(dāng)出現(xiàn)厄爾尼諾暖位相時(shí),南海北部的東北季風(fēng)將會(huì)減弱,南海南部的東北季風(fēng)則所受影響較小。我們還對(duì)南海海面風(fēng)場第一模態(tài)的時(shí)間序列與EMI指數(shù)進(jìn)行超前/滯后相關(guān)計(jì)算,得出當(dāng)時(shí)間序列滯后4個(gè)月時(shí)出現(xiàn)相關(guān)系數(shù)最大值0.58,低于與NIN~O3.4的相關(guān)系數(shù)。南海海面風(fēng)應(yīng)力第一模態(tài)的分布情況與風(fēng)場第一模態(tài)相近。圖7d為相應(yīng)的南海海面風(fēng)應(yīng)力旋度第一模態(tài)的時(shí)間序列與NIN~O3指數(shù)的對(duì)比圖,當(dāng)NIN~O3指數(shù)滯后風(fēng)應(yīng)力旋度第一模態(tài)時(shí)間序列3個(gè)月時(shí)出現(xiàn)相關(guān)系數(shù)最大值0.74,置信水平也超過99%。此外,南海海面風(fēng)應(yīng)力旋度場第一模態(tài)的時(shí)間序列與NIN~O3.4滯后2個(gè)月的相關(guān)系數(shù)達(dá)到最大值0.73,低于與NIN~O3的相關(guān)系數(shù)。對(duì)以上的相關(guān)關(guān)系進(jìn)行顯著性檢驗(yàn),由結(jié)果可知其顯著性水平都低于1%。

    圖7 南海海面風(fēng)場和風(fēng)應(yīng)力第一EOF模態(tài)的空間形態(tài)和時(shí)間變化Fig.7 The first EOF modes of the surface wind and wind stress variability

    圖8 為南海海面風(fēng)場和海面風(fēng)應(yīng)力旋度場的第二模態(tài)的空間形態(tài)和時(shí)間序列,由表2可知,風(fēng)場的第二模態(tài)方差貢獻(xiàn)率為31%,風(fēng)應(yīng)力的第二模態(tài)方差貢獻(xiàn)率為14%。圖8a為南海海面風(fēng)場第二模態(tài)空間分布情況,不同于海面風(fēng)場的第一模態(tài),由圖可以看出整個(gè)南海海域并無明顯渦旋,整體呈現(xiàn)西南風(fēng)趨勢,且中部風(fēng)速要大于南部和北部。這表明在厄爾尼諾正位相期間東北季風(fēng)會(huì)減弱、西南季風(fēng)會(huì)增強(qiáng)。圖8b為南海海面風(fēng)應(yīng)力旋度場第二模態(tài)空間分布情況,由圖可知在呂宋海峽以西和越南沿海存在一定的氣旋性結(jié)構(gòu),在西沙群島一帶出現(xiàn)一個(gè)反氣旋中心,可能有利于西沙暖渦[20]的形成。圖8c為南海海面風(fēng)場第二模態(tài)時(shí)間序列與NIN~O3的對(duì)比圖,由圖可見,當(dāng)風(fēng)場超前時(shí)間為1個(gè)月時(shí)具有最大的相關(guān)系數(shù)0.52,置信水平超過99%。圖8d為相應(yīng)的南海海面風(fēng)應(yīng)力旋度第二模態(tài)的時(shí)間序列與負(fù)IOD指數(shù)的對(duì)比圖,當(dāng)滯后5個(gè)月時(shí)出現(xiàn)相關(guān)系數(shù)負(fù)最大值-0.63,置信水平也超過99%。對(duì)上述兩個(gè)相關(guān)關(guān)系進(jìn)行顯著性檢驗(yàn),由結(jié)果可知二者的顯著性水平都低于1%。

    圖8 南海海面風(fēng)場和風(fēng)應(yīng)力第二EOF模態(tài)的空間形態(tài)和時(shí)間變化Fig.8 The second EOF modes of the surface wind and wind stress variability

    4.2.3 高度場

    南海海面高度場第一模態(tài)的方差貢獻(xiàn)率為42%,第二、第三模態(tài)的方差貢獻(xiàn)率分別為17%和14%(表2)。圖9a和圖9c分別為南海海面高度場第一模態(tài)的空間分布和時(shí)間序列。由圖9a可以看出,第一模態(tài)的空間形態(tài)基本上全海域?yàn)檎?極大值區(qū)域位于呂宋海峽以西和越南南部海岸,沿南中國大陸架海面高度異常值較低,分別與圖7b中反氣旋性、氣旋性結(jié)構(gòu)對(duì)應(yīng)。由圖9c可知,當(dāng)海面高度場第一模態(tài)的時(shí)間序列滯后EMI指數(shù)3個(gè)月時(shí),相關(guān)系數(shù)達(dá)負(fù)最大值-0.32,相關(guān)程度較低,其顯著性水平超過10%。圖9b和圖9d分別為南海海面高度場第二模特的空間形態(tài)和時(shí)間序列。由圖9b可知,海面高度異常極大值區(qū)域與圖7b中的反氣旋性結(jié)構(gòu)對(duì)應(yīng),出現(xiàn)在西沙群島附近和呂宋島西側(cè),此外,以極大值區(qū)域?yàn)榻?海面高度異常向南北兩側(cè)逐漸減小。由圖9d可得,當(dāng)海面高度場第二模態(tài)的時(shí)間序列超前NIN~O3指數(shù)5個(gè)月時(shí),相關(guān)系數(shù)達(dá)負(fù)最大值,為-0.37,相關(guān)程度也較低,顯著性水平在1%~5%。由此可見,南海海面高度變化受海面風(fēng)應(yīng)力旋度影響,反氣旋性中心對(duì)應(yīng)海面高度較低。海面高度與厄爾尼諾事件的關(guān)聯(lián)主要表現(xiàn)在區(qū)域性整體的變化,并且與中太平洋厄爾尼諾相關(guān)性最高,且為負(fù)相關(guān);亦即,當(dāng)中太平洋厄爾尼諾處于正位相時(shí),南海海平面會(huì)降低。

    圖9 南海海面高度第一和第二EOF模態(tài)的空間形態(tài)和時(shí)間變化Fig.9 The first and second EOF modes of the sea surface height variability

    4.2.4 溫度場

    南海海表溫度場第一、二模態(tài)的方差貢獻(xiàn)率分別為68%和14%。圖10a和圖10c為南海海表溫度場EOF第一模態(tài)的空間形態(tài)圖和時(shí)間序列。由圖10a可以看出,第一模態(tài)的空間分布整體上為西高東低,最大值出現(xiàn)在中南半島東南外海,反映出厄爾尼諾年第二年夏季這里的上升流會(huì)減弱。南海海表溫度場第一模態(tài)的時(shí)間序列與NIN~O3指數(shù)相比(圖10c)可知,當(dāng)時(shí)間序列滯后NIN~O3指數(shù)4個(gè)月時(shí)相關(guān)系數(shù)達(dá)到最大值0.62,置信水平超過99%。圖10b為南海海表溫度第二模態(tài)的空間分布圖,與第一模態(tài)不同,第二模態(tài)的分布為北高南低,量值較小。由圖10d可知,第二模態(tài)的時(shí)間序列滯后IOD指數(shù)1個(gè)月時(shí)相關(guān)系數(shù)達(dá)到最大值0.69,置信水平也超過99%。這兩個(gè)模態(tài)的顯著性水平都低于1%。由此可見,南海海表溫度年際變化主要受控于太平洋的變異,印度洋的變異對(duì)其也有一定的影響。

    圖10 南海海表溫度第一和第二EOF模態(tài)的空間形態(tài)和時(shí)間變化Fig.10 The first and second EOF modes of the sea surface temperature variability

    5 結(jié) 語

    本文基于衛(wèi)星觀測資料分析了2002-2011年期間南海海面風(fēng)場、高度場和溫度場的變化趨勢和年際變化特征。其中,南海區(qū)域平均海面風(fēng)速東分量的線性趨勢值為(-0.012±0.014)(m·s-1)·a-1,北分量為(-0.014±0.019)(m·s-1)·a-1。南海區(qū)域平均海面風(fēng)應(yīng)力旋度的線性趨勢值為(0.099±0.33)×10-9(N·m-3)·a-1。南海區(qū)域平均海面高度的線性趨勢值為(0.665±0.200)cm·a-1,大約是2005-2013年間全球海面高度線性趨勢的2.4倍。南海北部海盆為南海海面高度上升速率較大的區(qū)域。南海區(qū)域平均海表溫度線性趨勢值為(-0.016±0.017)℃·a-1,表明近十年南海海表溫度呈下降趨勢,與20世紀(jì)末以來全球變暖的停滯相一致。南海海面溫度趨勢分布為東高西低,特別呂宋海峽附近出現(xiàn)正增溫。

    從本文的分析結(jié)果來看,南海海面風(fēng)場、高度場和溫度場的年際變化與ENSO變化具有一定相關(guān)性。南海海面風(fēng)場EOF第一模態(tài)時(shí)間序列的變化比NIN~O3.4指數(shù)滯后5個(gè)月,二者相關(guān)系數(shù)為0.68;比EMI滯后4個(gè)月,相關(guān)系數(shù)為0.58。海面風(fēng)應(yīng)力旋度場第一模態(tài)時(shí)間序列的變化比NIN~O3指數(shù)滯后3個(gè)月,相關(guān)系數(shù)為0.74。海面風(fēng)場第二模態(tài)時(shí)間序列的變化比NIN~O3指數(shù)超前1個(gè)月,相關(guān)系數(shù)為0.52。海面風(fēng)應(yīng)力旋度場第二模態(tài)時(shí)間序列的變化比IOD指數(shù)滯后5個(gè)月,相關(guān)系數(shù)為-0.63。南海海面高度場EOF第一模態(tài)時(shí)間序列的變化比EMI滯后3個(gè)月,相關(guān)系數(shù)為-0.32。高度場第一模態(tài)空間分布為東高西低,在越南東南外海出現(xiàn)一個(gè)海面高度較高的區(qū)域。高度場第二模態(tài)時(shí)間序列的變化比NIN~O3指數(shù)超前5個(gè)月,相關(guān)系數(shù)為-0.37。高度場第二模態(tài)空間分布為中部高,南北低。南海海表溫度場第一模態(tài)時(shí)間序列的變化比NIN~O3指數(shù)滯后4個(gè)月,相關(guān)系數(shù)為0.62。溫度場第一模態(tài)空間分布為東低西高,最大值出現(xiàn)在中南半島東南部。溫度場第二模態(tài)的空間分布為北高南低,時(shí)間序列的變化比IOD指數(shù)滯后1個(gè)月,相關(guān)系數(shù)為0.69。

    南海海面風(fēng)場、高度場和溫度場在2002—2011年期間的變化趨勢與1993—2003年期間[11]相比,一個(gè)最重大的差別是海表溫度由上升轉(zhuǎn)為下降,這與20世紀(jì)末以來全球變暖的停滯[14]相一致。然而在這兩個(gè)時(shí)期南海海面高度卻在持續(xù)上升。這說明在2002-2011年期間雖然海表溫度在下降,但海洋中的熱含量仍然在增加。這與新近的一些意見,認(rèn)為熱量進(jìn)入了海洋深層[14]相一致。從海面風(fēng)場來看,在1993—2003年期間有西向風(fēng)增強(qiáng)的趨勢,而2002—2011年期間有東向風(fēng)增強(qiáng)的趨勢,這是否與海表溫度的變化趨勢改變有關(guān),尚待進(jìn)一步研究。

    在年際變化方面,Fang等[11]和我們的研究都顯示了南海的變異與熱帶太平洋的變異即使在不同時(shí)期都密切相關(guān)。Fang等[15]的研究僅采用NIN~O3.4指數(shù)代表熱帶太平洋的變異,本文則更進(jìn)一步研究了南海的變異與NIN~O3.4、NIN~O3及EMI指數(shù)的關(guān)系,其中EMI指數(shù)是新近提出的。結(jié)果發(fā)現(xiàn),對(duì)于南海海表溫度的年際變化,無論是區(qū)域平均還是EOF第一主成份,均與以NIN~O3指數(shù)為代表的東太平洋厄爾尼諾相關(guān)程度最高;而對(duì)于南海海面高度的年際變化,無論是區(qū)域平均還是EOF第一主成份,均與以EMI指數(shù)為代表的中太平洋厄爾尼諾相關(guān)程度最高,而且呈負(fù)相關(guān)關(guān)系。這說明控制南海海表溫度和海面高度的主要熱力和動(dòng)力機(jī)制應(yīng)當(dāng)有所不同,值得進(jìn)一步研究。

    [1] CHAO S Y,SHAW P T,WU S Y.El Ni?o modulation of the South China Sea circulation[J].Progress in Oceanography,1996,38(1): 51-93.

    [2] WU C R,SHAW P T,CHAO S Y.Seasonal and interannual variations in the velocity field of the South China Sea[J].Journal of Oceanography,1998,54(4):361-372.

    [3] LIANG W D,JAN J C,TANG T Y.Climatological wind and upper ocean heat content in the South China Sea[J].Acta Oceanography, 2000,38:91-114.

    [4] HWANG C,CHEN S A.Circulation and eddies over the South China Sea derived from TOPEX/Poseidon altimetry[J].Journal of Geophysical Research:Oceans,2000,105(C10):23943-23965.

    [5] HO C R,KUO N J,ZHENG Q,et al.Dynamically active areas in the South China Sea detected from TOPEX/POSEIDON satellite altimeter data[J].Remote Sensing of Environment,2000,71(3):320-328.

    [6] LI L,XU J D,CAI R S.Trends of sea level rise in the South China Sea during the 1990s:An altimetry result[J].Chinese Science Bulletin,2002,47(7):582-585.

    [7] CHU P C,LU S H,CHEN Y C.Temporal and spatial variabilities of the South China Sea surface temperature anomaly[J].Journal of Geophysical Research:Oceans,1997,102(C9):20937-20955.

    [8] KLEIN S A,SODEN B J,LAO N C.Remote sea surface temperature variations during ENSO:Evidence for a tropical atmospheric bridge [J].Journal of Climate,1999,12(4):917-932.

    [9] WANG D X,XIE Q,DU Y,et al.The 1997-1998 warm event in the South China Sea[J].Chinese Science Bulletin,2002,47(14): 1221-1227.

    [10] LIU Q Y,JIANG X,XIE S P,et al.A gap in the Indo-Pacific warm pool over the South China Sea in boreal winter:Seasonal development and interannual variability[J].Journal of Geophysical Research:Oceans,2004,109(C7):101-111.

    [11] FANG G H,CHEN H Y,WEI Z X,et al.Trends and interannual variability of the South China Sea surface winds,surface height,andsurface temperature in the recent decade[J].Journal of Geophysical Research Part 2006,111(C11):1-16.

    [12] CHURCH J A,WHITE N J,COLEMAN R,et al.Estimates of the regional distribution of sea level rise over the 1950-2000 period [J].Journal of Climate,2004,17(13):2609-2625.

    [13] FENG W,ZHONG M,XU H Z.Global sea level changes estimated from satellite altimetry,satellite gravimetry and Argo data during 2005-2013[J].Progress in Geophysics,2014,29(2):471-477.馮偉,鐘敏,許厚澤.聯(lián)合衛(wèi)星測高、衛(wèi)星重力和海洋浮標(biāo)資料研究2005-2013年的全球海平面變化[J].地球物理學(xué)進(jìn)展,2014,29(2):471-477.

    [14] MEEHL G A,ARBLASTER J M,FASULLO J T,et al.Model-based evidence of deep-ocean heat uptake during surface-temperature hiatus periods[J].Nature Climate Change,2011,1(7):360-364.

    [15] ASHOK K,BEHERA S K,RAO S A,et al.El Ni?o Modoki and its possible teleconnection[J].Journal of Geophysical Research:O-ceans,2007,112(C11):C11007.

    [16] LORENZ E N.Empirical orthogonal functions and statistical weather prediction[J].Science Report,1956,409(2):997-999.

    [17] KAIHATU J M,HANDLER R A,MARMORINO G O,et al.Empirical orthogonal function analysis of ocean surface currents using complex and real-vector methods[J].Journal of Atmospheric and Oceanic Technology,1998,15(4):927-941.

    [18] PICAUT J,HACHKERT E,BUSALACCHI A J,et al.Mechanisms of the 1997-1998 El Ni?o-La Ni?a,as inferred from space-based observations[J].Journal of Geophysical Research:Oceans,2002,107(C5):3037.

    [19] WANG B,WU R,FU X.Pacific-East Asian teleconnection:How does ENSO affect East Asian climate?[J].Journal of Climate,2000, 13(9):1517-1536.

    [20] CHU X,XUE H,QI Y,et al.An exceptional anticyclonic eddy in the South China Sea in 2010[J].Journal of Geophysical Research:O-ceans,2014,119(2):881-897.

    [21] LIAN Z,FANG G H,WEI Z X,et al.A comparison of wind stress datasets for the South China Sea[J].Ocean Dynamics,2015,65(5): 721-734.

    [22] XIE S P,XIE Q,WANG D,et al.Summer upwelling in the South China Sea and its role in regional climate variations[J].Journal of Geophysical Research:Oceans,2003,108(C8):343-367.

    T
    rends and Interannual Variability of the South China Sea Surface Winds,Surface Height and Surface Temperature in the Recent Decade

    WANG Jia-ying1,FANG Guo-hong1,2,3,WANG Yong-gang1,2,3
    (1.The First Institute of Oceanology,SOA,Qingdao 266061,China; 2.Key Lab of Marine Science and Numerical Modeling,SOA,Qingdao 266061,China; 3.Laboratory for Regional Oceanography and Numerical Modeling,Qingdao National Laboratory for Marine Science and Technology,Qingdao 266237,China)

    The monthly products of the surface winds(SW),sea surface height(SSH)and sea surface temperature(SST)of the South China Sea(SCS)from satellite observations during 2002-2011 are used to analyze trends and interannual variability.Time series are smoothed with a 12-month running mean filter.The linear trends of the east and north components of the regional mean SW,wind stress curl,SSH and SST are(-0.012±0.014)(m·s-1)·a-1,(-0.014±0.019)(m·s-1)·a-1,(0.099±0.33)×10-9(N·m-3)·a-1,(0.665±0.200)cm·a-1,(-0.016±0.017)℃·a-1,respectively.The interannual variability of regional mean SW,wind stress curl,SSH and SST correlates with NIN~O3,EMI(El Ni?o Modoki Index),EMI and NIN~O3 with coefficient of about 0.6,-0.85,-0.80 and 0.68,and with lag of 1,3, 2 and 4 months,respectively.The negative warming rate of the SCSSST is consistent with the global warming hiatus since the end of the last century.Empirical Orthogonal Function(EOF)analysis shows that the first EOF of the SW is characterized by a basin-wide anticyclonic pattern.The corresponding time coefficient function(TCF)correlates with the NIN~O3.4 index at the 99%confidence level(correlation coefficient=0.68),with a lag of 5 months.The first EOF of the SSH is characterized by high sea level along the east boundary and east of Vietnam.The corresponding TCF correlates with the EMI with a coefficient of -0.32 and a lag of 3 month.The first EOF of the SST is characterized by high values in the west and low values in the east.The corresponding TCF correlates with NIN~O3 index at 99%confidence level(correlation coefficient=0.63),with a lag of 4 months.

    South China Sea;trends;interannual variability

    February 29,2016

    P731.11

    A

    1671-6647(2017)02-0159-17

    10.3969/j.issn.1671-6647.2017.02.002

    2016-02-29

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目——南海環(huán)流和海峽水交換對(duì)海氣相互作用的影響(2011CB403502);國家自然科學(xué)基金委員會(huì)-山東省人民政府聯(lián)合資助海洋科學(xué)研究中心項(xiàng)目——海洋環(huán)境動(dòng)力學(xué)和數(shù)值模擬(U1406404)

    王佳瑩(1991-),女,遼寧鞍山人,碩士研究生,主要從事物理海洋學(xué)方面研究.E-mail:wangjy@fio.org.cn

    *通訊作者:方國洪(1939-),男,浙江瑞安人,研究員,博士,主要從事海洋潮汐、海洋環(huán)流和海洋數(shù)值模擬方面研究. E-mail:fanggh@fio.org.cn

    猜你喜歡
    海表年際風(fēng)場
    基于FLUENT的下?lián)舯┝魅S風(fēng)場建模
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    基于無人機(jī)的海表環(huán)境智能監(jiān)測系統(tǒng)設(shè)計(jì)與應(yīng)用
    2016與1998年春季北大西洋海表溫度異常的差異及成因
    融合海表溫度產(chǎn)品在渤黃東海的對(duì)比分析及初步驗(yàn)證
    太陽總輻照度對(duì)熱帶中太平洋海表溫度年代際變化的可能影響
    “最美風(fēng)場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    側(cè)向風(fēng)場中無人機(jī)的飛行研究
    亞洲夏季風(fēng)的年際和年代際變化及其未來預(yù)測
    與北大西洋接壤的北極海冰和年際氣候變化
    avwww免费| 高清在线国产一区| 国产精品国产高清国产av| 好男人电影高清在线观看| 男女床上黄色一级片免费看| 一区二区三区精品91| 亚洲avbb在线观看| 美女国产高潮福利片在线看| 精品国产一区二区三区四区第35| 大陆偷拍与自拍| 在线观看免费日韩欧美大片| 亚洲国产精品久久男人天堂| 欧美日韩精品网址| 婷婷六月久久综合丁香| www.999成人在线观看| 亚洲av第一区精品v没综合| 日韩av在线大香蕉| 最新在线观看一区二区三区| 国产亚洲精品一区二区www| 亚洲专区中文字幕在线| 亚洲精品粉嫩美女一区| bbb黄色大片| 人妻丰满熟妇av一区二区三区| 国产精品亚洲一级av第二区| 亚洲片人在线观看| 男女下面插进去视频免费观看| 亚洲免费av在线视频| 色播在线永久视频| 午夜福利18| videosex国产| 女人精品久久久久毛片| 在线十欧美十亚洲十日本专区| 涩涩av久久男人的天堂| 大型av网站在线播放| 亚洲成国产人片在线观看| 9191精品国产免费久久| 国产在线观看jvid| 亚洲av熟女| www国产在线视频色| 咕卡用的链子| 国产男靠女视频免费网站| 色综合站精品国产| 精品国产一区二区久久| 韩国av一区二区三区四区| 亚洲 国产 在线| 午夜精品国产一区二区电影| 中文字幕人妻丝袜一区二区| 亚洲人成77777在线视频| 老鸭窝网址在线观看| 成熟少妇高潮喷水视频| 亚洲五月色婷婷综合| 国产主播在线观看一区二区| 十八禁人妻一区二区| 成人18禁在线播放| 亚洲成人久久性| 国产精品一区二区三区四区久久 | 欧美激情高清一区二区三区| aaaaa片日本免费| 中文字幕高清在线视频| 亚洲自拍偷在线| 99国产综合亚洲精品| 欧美色欧美亚洲另类二区 | 国产av精品麻豆| 极品人妻少妇av视频| 国产99久久九九免费精品| 在线观看免费视频日本深夜| 国产欧美日韩精品亚洲av| 国产欧美日韩一区二区精品| 欧美最黄视频在线播放免费| 人人妻,人人澡人人爽秒播| 国产精品爽爽va在线观看网站 | 午夜成年电影在线免费观看| 激情视频va一区二区三区| 精品乱码久久久久久99久播| 在线观看免费午夜福利视频| 午夜福利影视在线免费观看| 在线观看免费视频网站a站| 亚洲国产精品成人综合色| 一卡2卡三卡四卡精品乱码亚洲| 日韩精品中文字幕看吧| 免费在线观看完整版高清| 18禁裸乳无遮挡免费网站照片 | 精品卡一卡二卡四卡免费| 亚洲少妇的诱惑av| 免费搜索国产男女视频| 成人18禁在线播放| 亚洲 欧美 日韩 在线 免费| 午夜福利免费观看在线| 久久久精品国产亚洲av高清涩受| 一本综合久久免费| 国产av一区二区精品久久| 1024视频免费在线观看| 窝窝影院91人妻| 91字幕亚洲| 免费高清在线观看日韩| 熟女少妇亚洲综合色aaa.| 精品福利观看| 88av欧美| 亚洲免费av在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 看免费av毛片| 久久国产乱子伦精品免费另类| 51午夜福利影视在线观看| 精品日产1卡2卡| 国产精品国产高清国产av| 久久精品国产亚洲av香蕉五月| 视频在线观看一区二区三区| 黄色毛片三级朝国网站| 日韩成人在线观看一区二区三区| 色播亚洲综合网| 高清在线国产一区| 国产高清视频在线播放一区| or卡值多少钱| 校园春色视频在线观看| 亚洲三区欧美一区| 久久精品人人爽人人爽视色| 法律面前人人平等表现在哪些方面| 久久国产精品影院| 国产黄a三级三级三级人| 黄色片一级片一级黄色片| 香蕉久久夜色| 男人舔女人下体高潮全视频| 这个男人来自地球电影免费观看| 淫妇啪啪啪对白视频| 纯流量卡能插随身wifi吗| 国产精品自产拍在线观看55亚洲| 午夜成年电影在线免费观看| 少妇粗大呻吟视频| 高清黄色对白视频在线免费看| 丰满人妻熟妇乱又伦精品不卡| 黄网站色视频无遮挡免费观看| 19禁男女啪啪无遮挡网站| 国产精品秋霞免费鲁丝片| 国产伦一二天堂av在线观看| 国产精品久久久久久精品电影 | 一级毛片女人18水好多| 美女免费视频网站| 成年人黄色毛片网站| 麻豆一二三区av精品| 亚洲美女黄片视频| 十八禁人妻一区二区| 成人欧美大片| 精品国产一区二区三区四区第35| 久久久国产成人免费| 国产亚洲精品av在线| 99re在线观看精品视频| 两性夫妻黄色片| 免费一级毛片在线播放高清视频 | 亚洲精品国产色婷婷电影| 亚洲中文字幕一区二区三区有码在线看 | 国产精品永久免费网站| 亚洲 国产 在线| 日本免费a在线| 日本在线视频免费播放| 天天躁夜夜躁狠狠躁躁| 久久亚洲真实| 亚洲色图综合在线观看| 搡老岳熟女国产| 这个男人来自地球电影免费观看| 香蕉久久夜色| 国产欧美日韩一区二区三| 国产av又大| 亚洲国产日韩欧美精品在线观看 | 最近最新免费中文字幕在线| 国产亚洲精品综合一区在线观看 | 国产野战对白在线观看| 亚洲九九香蕉| 操美女的视频在线观看| 啦啦啦观看免费观看视频高清 | 极品教师在线免费播放| 国产极品粉嫩免费观看在线| 最近最新中文字幕大全电影3 | 国产精品久久久久久亚洲av鲁大| 日韩欧美免费精品| 欧美人与性动交α欧美精品济南到| 久久久久精品国产欧美久久久| 国产伦一二天堂av在线观看| 欧美+亚洲+日韩+国产| 国产成年人精品一区二区| 久久久久久久久久久久大奶| 最新在线观看一区二区三区| 精品久久久久久久人妻蜜臀av | 99国产精品免费福利视频| 人成视频在线观看免费观看| bbb黄色大片| 啦啦啦观看免费观看视频高清 | 真人一进一出gif抽搐免费| 99精品在免费线老司机午夜| 久久人人97超碰香蕉20202| 国产精品日韩av在线免费观看 | 久久中文看片网| 亚洲av日韩精品久久久久久密| 真人一进一出gif抽搐免费| 日本在线视频免费播放| 久久久久国产精品人妻aⅴ院| 久久中文看片网| 操美女的视频在线观看| 免费在线观看影片大全网站| 久久精品91无色码中文字幕| 美国免费a级毛片| 精品国产美女av久久久久小说| av福利片在线| 变态另类成人亚洲欧美熟女 | 两性午夜刺激爽爽歪歪视频在线观看 | 在线观看免费午夜福利视频| 国产一区二区三区综合在线观看| 这个男人来自地球电影免费观看| 成人国语在线视频| 黄色a级毛片大全视频| 国产精品98久久久久久宅男小说| 亚洲国产精品sss在线观看| 成人亚洲精品一区在线观看| 久久午夜亚洲精品久久| 成人18禁高潮啪啪吃奶动态图| 美女国产高潮福利片在线看| 久久国产精品男人的天堂亚洲| 欧美日韩亚洲综合一区二区三区_| 午夜老司机福利片| 国产成人av激情在线播放| 不卡一级毛片| 精品乱码久久久久久99久播| 一个人观看的视频www高清免费观看 | 又黄又粗又硬又大视频| 国产精品日韩av在线免费观看 | 亚洲av成人av| 日本三级黄在线观看| 91字幕亚洲| 午夜福利在线观看吧| 久久精品国产综合久久久| 午夜免费观看网址| 大码成人一级视频| 此物有八面人人有两片| 一区在线观看完整版| 欧美激情久久久久久爽电影 | 女性被躁到高潮视频| 国产一区二区在线av高清观看| 免费久久久久久久精品成人欧美视频| 色在线成人网| 精品一区二区三区四区五区乱码| 亚洲国产欧美日韩在线播放| 午夜免费激情av| 午夜福利影视在线免费观看| 久久久久久久午夜电影| 亚洲在线自拍视频| 丁香欧美五月| 国产97色在线日韩免费| 好男人电影高清在线观看| 又大又爽又粗| 一边摸一边抽搐一进一出视频| 久久久久久国产a免费观看| 国产精品永久免费网站| 精品国产一区二区三区四区第35| 国产区一区二久久| 自线自在国产av| a在线观看视频网站| 欧美 亚洲 国产 日韩一| 免费高清视频大片| 亚洲情色 制服丝袜| 一级a爱片免费观看的视频| 免费在线观看黄色视频的| 97人妻天天添夜夜摸| 亚洲人成77777在线视频| 久久人人爽av亚洲精品天堂| 午夜免费观看网址| 久久伊人香网站| 99国产极品粉嫩在线观看| 久久九九热精品免费| 精品福利观看| 亚洲aⅴ乱码一区二区在线播放 | 99精品在免费线老司机午夜| 男女床上黄色一级片免费看| 国产日韩一区二区三区精品不卡| 国产精品久久久久久精品电影 | 曰老女人黄片| 国产精品自产拍在线观看55亚洲| 99国产精品99久久久久| 亚洲熟妇熟女久久| 国产精品一区二区免费欧美| 精品人妻在线不人妻| aaaaa片日本免费| 身体一侧抽搐| 午夜福利视频1000在线观看 | 制服诱惑二区| 妹子高潮喷水视频| 一级,二级,三级黄色视频| 亚洲自偷自拍图片 自拍| 日韩三级视频一区二区三区| 看免费av毛片| АⅤ资源中文在线天堂| 淫妇啪啪啪对白视频| 久热这里只有精品99| 可以免费在线观看a视频的电影网站| 99久久久亚洲精品蜜臀av| 一级,二级,三级黄色视频| 不卡av一区二区三区| 90打野战视频偷拍视频| 精品高清国产在线一区| 亚洲精品国产精品久久久不卡| 制服人妻中文乱码| 在线永久观看黄色视频| 亚洲国产精品久久男人天堂| avwww免费| 久久精品国产99精品国产亚洲性色 | 国产欧美日韩精品亚洲av| 亚洲欧洲精品一区二区精品久久久| 巨乳人妻的诱惑在线观看| 无人区码免费观看不卡| 国产欧美日韩一区二区三区在线| 男人操女人黄网站| 激情视频va一区二区三区| 色精品久久人妻99蜜桃| 亚洲欧美精品综合一区二区三区| 国产精品免费一区二区三区在线| 成人亚洲精品av一区二区| 久久人人精品亚洲av| 日韩有码中文字幕| 久久香蕉精品热| 精品免费久久久久久久清纯| 色播在线永久视频| 嫩草影院精品99| 国产1区2区3区精品| 中文字幕人妻丝袜一区二区| 久久人妻av系列| 久久婷婷人人爽人人干人人爱 | www.999成人在线观看| av天堂久久9| 国产成人系列免费观看| 亚洲成国产人片在线观看| av中文乱码字幕在线| 国产国语露脸激情在线看| 亚洲一区二区三区色噜噜| 国产一区二区三区视频了| 两个人视频免费观看高清| 搡老熟女国产l中国老女人| 长腿黑丝高跟| 日日夜夜操网爽| 欧美黑人欧美精品刺激| 黄片播放在线免费| 深夜精品福利| 国产成人av激情在线播放| 久久中文看片网| 国产精品 欧美亚洲| 99re在线观看精品视频| 亚洲男人的天堂狠狠| 国产在线观看jvid| 欧美日韩亚洲国产一区二区在线观看| 精品国产一区二区三区四区第35| 亚洲国产中文字幕在线视频| 在线永久观看黄色视频| 久久影院123| av视频免费观看在线观看| 免费久久久久久久精品成人欧美视频| 婷婷六月久久综合丁香| 亚洲欧洲精品一区二区精品久久久| 久久精品aⅴ一区二区三区四区| 精品久久久久久成人av| 亚洲国产看品久久| 色尼玛亚洲综合影院| 在线观看一区二区三区| 法律面前人人平等表现在哪些方面| 久久狼人影院| 黄色丝袜av网址大全| 国产精品久久久人人做人人爽| 国产成人啪精品午夜网站| 一a级毛片在线观看| 精品欧美一区二区三区在线| 精品一品国产午夜福利视频| 正在播放国产对白刺激| 黄片小视频在线播放| 天堂影院成人在线观看| 黄片小视频在线播放| 久久国产精品人妻蜜桃| 久久久水蜜桃国产精品网| 欧美最黄视频在线播放免费| 久久久国产精品麻豆| 真人做人爱边吃奶动态| 国产精品av久久久久免费| 亚洲午夜理论影院| 18禁美女被吸乳视频| 在线十欧美十亚洲十日本专区| 日日干狠狠操夜夜爽| 欧美性长视频在线观看| 啪啪无遮挡十八禁网站| 精品少妇一区二区三区视频日本电影| 给我免费播放毛片高清在线观看| 国产高清视频在线播放一区| 亚洲欧美激情综合另类| 搡老岳熟女国产| 色av中文字幕| 国产免费av片在线观看野外av| 色尼玛亚洲综合影院| 国产欧美日韩一区二区三| 日本三级黄在线观看| 免费av毛片视频| 精品第一国产精品| 啦啦啦 在线观看视频| 91精品国产国语对白视频| or卡值多少钱| 在线十欧美十亚洲十日本专区| 久久久久久久精品吃奶| 国产亚洲精品久久久久久毛片| 青草久久国产| 少妇粗大呻吟视频| 免费不卡黄色视频| 美女免费视频网站| 久久久久精品国产欧美久久久| 欧美老熟妇乱子伦牲交| 免费在线观看日本一区| 日韩视频一区二区在线观看| 熟妇人妻久久中文字幕3abv| 69av精品久久久久久| 亚洲色图av天堂| 日韩高清综合在线| 在线观看免费视频日本深夜| 久9热在线精品视频| 久久人人97超碰香蕉20202| 日本vs欧美在线观看视频| 午夜影院日韩av| 国产成人欧美在线观看| 亚洲精品久久成人aⅴ小说| 国产精品 欧美亚洲| 亚洲第一青青草原| 国产成人免费无遮挡视频| 天天添夜夜摸| 在线十欧美十亚洲十日本专区| 日本 欧美在线| 久久精品91蜜桃| 别揉我奶头~嗯~啊~动态视频| 成人国产一区最新在线观看| 亚洲第一欧美日韩一区二区三区| 一区福利在线观看| 黑人欧美特级aaaaaa片| 国产精品日韩av在线免费观看 | 高清黄色对白视频在线免费看| 亚洲精品一卡2卡三卡4卡5卡| 两个人看的免费小视频| 国产私拍福利视频在线观看| 国产精品野战在线观看| 久久人妻熟女aⅴ| 成年版毛片免费区| 如日韩欧美国产精品一区二区三区| 一区二区三区高清视频在线| 国产精品亚洲av一区麻豆| 国产真人三级小视频在线观看| 多毛熟女@视频| 岛国视频午夜一区免费看| 日日夜夜操网爽| 国产av精品麻豆| 亚洲午夜精品一区,二区,三区| 久久久久国内视频| 91字幕亚洲| 成年女人毛片免费观看观看9| 午夜a级毛片| 精品久久久久久久久久免费视频| 亚洲狠狠婷婷综合久久图片| 9热在线视频观看99| 青草久久国产| 国产av在哪里看| 自拍欧美九色日韩亚洲蝌蚪91| 国产又色又爽无遮挡免费看| 人人澡人人妻人| 久久久久久久午夜电影| 欧美日韩瑟瑟在线播放| 青草久久国产| 最好的美女福利视频网| www.精华液| 日本精品一区二区三区蜜桃| 午夜日韩欧美国产| 电影成人av| 麻豆久久精品国产亚洲av| 欧美色欧美亚洲另类二区 | 精品第一国产精品| 国产成人精品在线电影| 国产精品国产高清国产av| 国产精华一区二区三区| 精品人妻1区二区| 欧美成人免费av一区二区三区| 欧美老熟妇乱子伦牲交| 99久久99久久久精品蜜桃| 日韩av在线大香蕉| 99久久久亚洲精品蜜臀av| 成人国语在线视频| 久久中文看片网| 丝袜美足系列| 一级,二级,三级黄色视频| 午夜视频精品福利| 婷婷丁香在线五月| 午夜老司机福利片| 久久久久亚洲av毛片大全| 自拍欧美九色日韩亚洲蝌蚪91| 在线免费观看的www视频| 亚洲av片天天在线观看| 色播亚洲综合网| 黄色 视频免费看| 夜夜夜夜夜久久久久| 婷婷六月久久综合丁香| 一级毛片女人18水好多| 国产一区二区在线av高清观看| 亚洲 欧美一区二区三区| www.999成人在线观看| 欧美一级a爱片免费观看看 | 亚洲av片天天在线观看| 免费高清在线观看日韩| 欧美日本视频| 俄罗斯特黄特色一大片| 一级毛片高清免费大全| 女性生殖器流出的白浆| 麻豆一二三区av精品| 久久久国产精品麻豆| 中文字幕高清在线视频| 久久久久久大精品| 一区二区三区国产精品乱码| 久久香蕉国产精品| 黄色毛片三级朝国网站| 亚洲中文av在线| 午夜日韩欧美国产| 嫩草影院精品99| 欧美日本视频| 精品国产乱码久久久久久男人| 一本综合久久免费| 少妇裸体淫交视频免费看高清 | 国产精品爽爽va在线观看网站 | 级片在线观看| 免费av毛片视频| 免费在线观看黄色视频的| ponron亚洲| 亚洲男人的天堂狠狠| 欧美激情极品国产一区二区三区| 亚洲午夜精品一区,二区,三区| 别揉我奶头~嗯~啊~动态视频| 啦啦啦免费观看视频1| 国内精品久久久久精免费| 亚洲av第一区精品v没综合| 亚洲一区二区三区色噜噜| 成年女人毛片免费观看观看9| 人人妻人人爽人人添夜夜欢视频| 欧美绝顶高潮抽搐喷水| 波多野结衣一区麻豆| 国产精品久久久av美女十八| 啦啦啦观看免费观看视频高清 | 久久久久久免费高清国产稀缺| 国产精品亚洲av一区麻豆| 深夜精品福利| 欧美精品啪啪一区二区三区| 别揉我奶头~嗯~啊~动态视频| 国产又色又爽无遮挡免费看| 黄色毛片三级朝国网站| 国产精品久久视频播放| 国产av在哪里看| 精品第一国产精品| 18禁美女被吸乳视频| 国产成人欧美| 免费久久久久久久精品成人欧美视频| 琪琪午夜伦伦电影理论片6080| 亚洲av五月六月丁香网| 美女 人体艺术 gogo| 亚洲激情在线av| 美国免费a级毛片| 欧美乱妇无乱码| 欧美在线黄色| 一区二区三区国产精品乱码| 这个男人来自地球电影免费观看| 午夜影院日韩av| 夜夜躁狠狠躁天天躁| 国产高清有码在线观看视频 | 亚洲av成人一区二区三| 又大又爽又粗| 亚洲国产精品sss在线观看| 国产99久久九九免费精品| 久久亚洲真实| av中文乱码字幕在线| 精品久久久久久久久久免费视频| 国产成人精品久久二区二区91| 999久久久国产精品视频| 黄色丝袜av网址大全| 欧美日韩亚洲综合一区二区三区_| 国产一级毛片七仙女欲春2 | 在线观看66精品国产| 很黄的视频免费| 人人妻人人爽人人添夜夜欢视频| 欧美激情久久久久久爽电影 | avwww免费| 亚洲国产欧美网| 青草久久国产| av欧美777| 色综合婷婷激情| 久久性视频一级片| 国产不卡一卡二| 亚洲av成人一区二区三| 青草久久国产| 精品一品国产午夜福利视频| 国产亚洲精品av在线| 制服诱惑二区| 两个人免费观看高清视频| 亚洲一区二区三区不卡视频| 如日韩欧美国产精品一区二区三区| 国产精品香港三级国产av潘金莲| 波多野结衣巨乳人妻| 大香蕉久久成人网| 午夜福利18| 日日爽夜夜爽网站| 99riav亚洲国产免费| 国产xxxxx性猛交| 热re99久久国产66热| 给我免费播放毛片高清在线观看| 国产精品香港三级国产av潘金莲| 亚洲国产日韩欧美精品在线观看 | 日日摸夜夜添夜夜添小说| 久久天躁狠狠躁夜夜2o2o| 在线视频色国产色| 成人欧美大片| 久久久久九九精品影院| av在线播放免费不卡| 国产熟女xx| 操出白浆在线播放|