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

    2013年春季W(wǎng)yrtki急流變異及成因分析

    2017-06-22 14:36:14韓國慶段永亮王關(guān)鎖王輝武劉延亮王海員于衛(wèi)東
    海洋科學(xué)進展 2017年2期
    關(guān)鍵詞:海表赤道西風(fēng)

    韓國慶,劉 琳,2*,段永亮,2,王關(guān)鎖,王輝武,2,劉延亮,2,王海員,馮 琳,2,于衛(wèi)東,2

    (1.國家海洋局第一海洋研究所海洋與氣候研究中心,山東青島266061; 2.青島海洋科學(xué)與技術(shù)國家實驗室區(qū)域海洋動力學(xué)與數(shù)值模擬功能實驗室,山東青島266200; 3.國家海洋局第一海洋研究所海洋環(huán)境數(shù)值模擬實驗室,山東青島266061)

    2013年春季W(wǎng)yrtki急流變異及成因分析

    韓國慶1,劉 琳1,2*,段永亮1,2,王關(guān)鎖2,3,王輝武1,2,劉延亮1,2,王海員1,馮 琳1,2,于衛(wèi)東1,2

    (1.國家海洋局第一海洋研究所海洋與氣候研究中心,山東青島266061; 2.青島海洋科學(xué)與技術(shù)國家實驗室區(qū)域海洋動力學(xué)與數(shù)值模擬功能實驗室,山東青島266200; 3.國家海洋局第一海洋研究所海洋環(huán)境數(shù)值模擬實驗室,山東青島266061)

    通過分析深海潛標(biāo)系統(tǒng)的實際觀測海流數(shù)據(jù)對2013年赤道印度洋上層Wyrtki急流(WJ)的變異進行了研究,發(fā)現(xiàn)WJ在2013年呈現(xiàn)明顯異常變化,主要表現(xiàn)為春季W(wǎng)J強度略強于同年秋季W(wǎng)J。再分析資料的分析結(jié)果也同樣驗證此結(jié)論,而這與氣候態(tài)WJ的演化特征相反。在氣候平均態(tài)下,印度洋春季W(wǎng)J強度遠弱于秋季W(wǎng)J強度。進一步數(shù)據(jù)分析與數(shù)值模式結(jié)果表明,2013年春季W(wǎng)J異常加強與赤道印度洋海表風(fēng)場變化密切相關(guān),而春季海表風(fēng)場的變化主要歸因于2013年春季異常增強的季節(jié)內(nèi)振蕩(ISO)事件,發(fā)生在2013-05的ISO對流位相是導(dǎo)致春季印度洋海表西風(fēng)異常的直接原因。本研究揭示了中小時間尺度海氣相互作用可以影響更長時間尺度海洋環(huán)流系統(tǒng)的年際變化。

    Wyrtki急流;季節(jié)內(nèi)振蕩;印度洋;年際變化

    印度洋是世界第三大洋,其獨特的地理特征決定了它別具特色的季風(fēng)性氣候。印度洋北側(cè)被亞洲大陸阻擋,南部為廣袤開闊的南大洋,獨特的地形分布以及太陽輻射的年變化造就了印度洋處于強大的亞洲季風(fēng)系統(tǒng)影響下。在亞洲季風(fēng)系統(tǒng)控制下,印度洋特別是熱帶印度洋地區(qū)存在形式多樣的局地氣候及海洋環(huán)流系統(tǒng)[1-10]。印度洋冬季被東北季風(fēng)控制,夏季被強盛的西南季風(fēng)影響,春季和秋季分是冬季風(fēng)和夏季風(fēng)的轉(zhuǎn)換期。赤道印度洋每年春季(4—5月)和秋季(10—翌年1月)出現(xiàn)兩次東向Wyrtki急流(WJ)[11],它發(fā)生在表層且流幅窄、流速強,根據(jù)漂流浮標(biāo)觀測,2°S~2°N平均的氣候態(tài)流速在春、秋季分別為50和70 cm/s[12-13]。該海流系統(tǒng)主要發(fā)生在(50°~80°E,3°S~3°N)的赤道印度洋區(qū)域內(nèi),其深度主要為海洋上層,從100 m深度至海面[14]。WJ影響熱帶東西印度洋之間上層海洋的水體、鹽度和熱量的東西向輸送,并對印度洋海盆尺度海氣相互作用產(chǎn)生影響[15-16]。WJ存在顯著年際變化[17-18],WJ年際變化受印度洋偶極子事件(IOD)和ENSO影響[19-20]并對IOD產(chǎn)生作用[15]。

    相對于太平洋與大西洋,印度洋海洋環(huán)境監(jiān)測及觀測均處于起步階段。前人只能通過有限的衛(wèi)星遙感海面資料以及稀少的現(xiàn)場觀測資料對WJ演化特征及影響機制進行研究,更多則是利用數(shù)值模式對WJ的演化特征、發(fā)生和發(fā)展規(guī)律及控制因素進行分析。2000年后印度洋海洋觀測系統(tǒng)的發(fā)展,特別是它的重要組成部分“非洲-亞洲-澳大利亞季風(fēng)分析和預(yù)測研究錨系浮標(biāo)陣列(RAMA)”的建立為研究WJ的三維流速、溫-鹽結(jié)構(gòu)和多時間尺度變化等諸多細結(jié)構(gòu)提供了第一手資料?;赗AMA浮標(biāo)數(shù)據(jù),目前對WJ流量[21]、季節(jié)變化[22]和年際變化動力學(xué)[23]有了進一步認識。本文中,我們利用多普勒聲學(xué)剖面儀(ADCP)觀測的印度洋上層海洋海流剖面資料,結(jié)合多個海洋再分析數(shù)據(jù)集,對2013年春季W(wǎng)J進行了分析,發(fā)現(xiàn)2013年春季W(wǎng)J強度與秋季W(wǎng)J強度相比略強,這與氣候態(tài)WJ強度的特征不符。我們進一步使用了POM數(shù)值海洋模式對2013年春季W(wǎng)J異常產(chǎn)生原因進行了分析和探討。

    (李 燕 編輯)

    1 數(shù)據(jù)及模式介紹

    1.1 數(shù)據(jù)簡介

    本研究使用了國家海洋局“全球變化與海氣相互作用專項”2013年印度洋南部水體環(huán)境綜合調(diào)查春季航次所布放的深水潛標(biāo)觀測數(shù)據(jù)。該航次于2013-03—05執(zhí)行,并于2013-04在赤道中印度洋(85°E,0°)處布放了一套深水潛標(biāo)系統(tǒng)(圖1),現(xiàn)場水深4 200 m,潛標(biāo)主浮體設(shè)計深度350 m。主浮體上安裝的一套ADCP(頻率為150 K)可以實現(xiàn)對海洋上層海流剖面進行連續(xù)觀測。ADCP觀測設(shè)置為垂直間隔16 m,時間間隔1 h。我們對ADCP原始資料進行了插值處理,海流數(shù)據(jù)垂向間隔插值為10 m,時間分辨率為日平均。由于ADCP在海面附近受海表反射聲信號影響較大,本文中我們忽略了上層40 m內(nèi)數(shù)據(jù),只采用40~150 m深度觀測結(jié)果。ADCP從2013-04-05開始記錄數(shù)據(jù),2014-04-18終止記錄,時間范圍完整覆蓋了2013年春季、秋季W(wǎng)J過程。

    圖1 深水潛標(biāo)位置示意圖Fig.1 Location of mooring system

    本文使用了Ocean Surface Current Analyses Real Time(OSCAR,http:∥www.oscar.noaa.gov/index. html)、Global Ocean Data Assimilation System(GODAS,http:∥www.cpc.ncep.noaa.gov/products/GODAS/)和Geophysical Fluid Dynamics Laboratory(GFDL,http:∥www.gfdl.noaa.gov/)三套不同海洋再分析資料以及RAMA中位于(80°30'E,0°)點處的浮標(biāo)觀測數(shù)據(jù)(http:∥www.pmel.noaa.gov/tao/rama/)。OSCAR,GFDL和GODAS時間分辨率為月平均,RAMA數(shù)據(jù)時間分辨率為日平均。OSCAR資料選取海洋上層15 m處水平流速,空間水平分辨率為1°×1°,使用資料時間跨度為2013-01—2013-12,覆蓋區(qū)域為30°30'~119°30'E,29°30'S~29°30'N。GODAS資料選取海洋上層15 m處水平流速,空間水平分辨率為1°×(1/3)°,使用資料時間跨度為2013-01—12,覆蓋區(qū)域為1°~360°E,74°S~65°N。GFDL資料選取海洋上層15 m處水平流速,空間水平分辨率為1°×1°,使用資料時間跨度為2013-01—12,覆蓋區(qū)域為279°~80° E,81°S~90°N。本研究使用的風(fēng)場數(shù)據(jù)為NCEP/NCAR再分析資料6 h平均三維風(fēng)場數(shù)據(jù)(http:∥www. esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.html)。NCEP/NCAR再分析資料是由美國國家環(huán)境預(yù)報中心和國家大氣研究中心聯(lián)合開發(fā)的數(shù)據(jù)產(chǎn)品,水平空間分辨率為2°30'×2°30'。使用資料時間跨度為1958-01—2015-12,覆蓋區(qū)域為0°~357°30'E,90°S~90°N。向外長波輻射(OLR)數(shù)據(jù)來自于美國國家海洋和大氣管理局(NOAA)(http:∥www.esrl.noaa.gov/psd/data/gridded/data.interp_OLR.html),空間水平分辨率為2°30'×2°30',使用資料時間跨度為1974-01—2013-12,覆蓋區(qū)域為0°00'~357°30'E,90°S~90°N。

    1.2 POM模式簡介

    Princeton Ocean Model(POM)是由美國普林斯頓大學(xué)Blumberg和Mellor在1977年共同建立起來的一個三維斜壓原始方程數(shù)值海洋模式,被當(dāng)今國內(nèi)外應(yīng)用較為廣泛的河口、近岸海洋模式,模式模擬內(nèi)容主要有:潮流、風(fēng)生流、混合層和躍層、熱鹽環(huán)流、海洋環(huán)流和輸運及與AGCM(大氣環(huán)流模式)和冰模式的耦合。垂向混合系數(shù)由二階湍流閉合模型確定,這在一定程度上擺脫了人為因素的困擾。垂直方向采用σ坐標(biāo),水平網(wǎng)格采用的是曲線正交坐標(biāo)系統(tǒng),變量空間配置使用“Arakawa C”網(wǎng)格,可以較好的匹配岸界。水平時間差分是顯式的,而垂直時間差分是隱式的,這可以保證垂向不受CFL條件限制,從而保證模式在海洋表層和底層可以有很高的垂向分辨率。模式區(qū)域設(shè)定在0°~360°E,89°45'S~90°15'N范圍內(nèi)。采用直角網(wǎng)格,東西劃分為721個網(wǎng)格,南北劃分為361個網(wǎng)格,空間水平分辨率為30'×30',垂向分為21層。模式采用實際水深。

    利用POM模式,我們設(shè)計了兩組數(shù)值實驗,控制實驗和敏感性實驗??刂茖嶒灷脷夂驊B(tài)大氣強迫場數(shù)據(jù)對POM模式進行強迫,敏感性實驗利用2013年大氣強迫場對模式進行強迫。兩組數(shù)值實驗前均采用NCEP風(fēng)場對模式進行強迫,從2008-01-01—2012-12-31,利用當(dāng)日輸出值作為兩組數(shù)值實驗的初始條件。

    2 觀測資料分析

    圖2a給出了RAMA浮標(biāo)(80°30'E,0°)處緯向海流的氣候態(tài)時間演化圖。氣候態(tài)上來看,每年4月和11月赤道印度洋均出現(xiàn)WJ,且春季W(wǎng)J強度低于秋季W(wǎng)J。圖2b給出2013年ADCP觀測的日平均緯向流時間演化圖。如圖所示,ADCP完整觀測到2013年發(fā)生的2次WJ事件:春季W(wǎng)J盛期出現(xiàn)在2013-05,最大值達1.8 m/s,流速最大區(qū)覆蓋范圍最深處達110 m,強流區(qū)一直延伸到海洋內(nèi)部140 m處;春季W(wǎng)J自2013-05-25起開始減弱,強度由強盛期的1.8 m/s逐漸減弱為0.6 m/s;2013-06-04春季W(wǎng)J消亡。秋季W(wǎng)J于2013-11-26出現(xiàn),其東向流速大于0.5 m/s。2013-12月初,WJ迅速發(fā)展,強度達1 m/s;2013-12-20,秋季W(wǎng)J開始減弱,強度減弱為0.8 m/s;2013-12-25,秋季W(wǎng)J消亡;2014-01-04,赤道印度洋重新被西向流控制。ADCP觀測結(jié)果顯示,2013年WJ發(fā)展與WJ在氣候平均態(tài)情況下的春季弱、秋季強顯著不同,并且這一變化也與最近研究發(fā)現(xiàn)的春季W(wǎng)J減弱趨勢[24]不同,因此值得進一步關(guān)注其發(fā)生和發(fā)展原因。

    為了進一步了解ADCP實測資料揭示的2013年異常WJ現(xiàn)象,接下來我們分析了不同來源的3套海洋再分析資料產(chǎn)品,從而對2013年WJ變化的空間分布特征進行全面把握。

    圖2 緯向海流時間-深度剖面圖Fig.2 Temporal evolution of zonal velocity

    圖3 經(jīng)向平均的赤道印度洋緯向流時間-經(jīng)度分布圖Fig.3 Temporal evolution of zonal velocity along the equator

    圖3給出了2013年OSCAR,GFDL,GODAS資料刻畫的赤道印度洋緯向流的逐月變化。通過與ADCP實測資料比對(圖3d),我們發(fā)現(xiàn)3套再分析資料具有一致性的季節(jié)變化,即可以揭示春季W(wǎng)J與秋季W(wǎng)J以及春季W(wǎng)J強度略強于秋季W(wǎng)J,因此可以認為3套再分析資料可以較好再現(xiàn)2013年春季和秋季W(wǎng)J事件。雖然3套再分析資料能夠再現(xiàn)2013年春季和秋季W(wǎng)J現(xiàn)象,但是三者之間也存在顯著差異。與ADCP實測結(jié)果相比較,GODAS資料對春季W(wǎng)J刻畫強于觀測約0.3 m/s,GFDL與OSCAR對春季W(wǎng)J現(xiàn)象的描述比觀測弱。對秋季W(wǎng)J刻畫上,3套再分析資料產(chǎn)品均表現(xiàn)出弱于觀測結(jié)果,幅度大約為0.5 m/s。雖然3套再分析資料與觀測結(jié)果存在一定偏差,但是3套再分析資料產(chǎn)品均可以較好把握2013年春季W(wǎng)J事件。因此我們可以從3套不同再分析資料產(chǎn)品出發(fā),對2013年春季W(wǎng)J現(xiàn)象的空間分布及時間演化特征進行分析。圖3a、圖3b、圖3c分別是OSCAR,GFDL,GODAS三套再分析資料產(chǎn)品中緯向流在2°30'S~2°30'N范圍內(nèi)經(jīng)向平均的結(jié)果。結(jié)果表明,2013年春季W(wǎng)J在5月出現(xiàn),且強度強于秋季W(wǎng)J。其中GODAS資料中WJ現(xiàn)象最為明顯,強度最強,最大值達1.5 m/s。GFDL資料中WJ強度最弱,海流最大值只有0.8 m/s,并且在GFDL資料中,2013-10—11沒有出現(xiàn)WJ,僅在12月出現(xiàn),WJ主要覆蓋60°~85°E之間區(qū)域。從3套不同的再分析資料產(chǎn)品可以看出,2013年春季W(wǎng)J強度均明顯強于秋季W(wǎng)J,這與平均意義下的結(jié)果有所差別(圖2a)。

    Wyrtki最早提出WJ受季風(fēng)轉(zhuǎn)換期的赤道緯向風(fēng)影響,是西風(fēng)強迫下海洋的響應(yīng)[11]。Han等通過不同復(fù)雜性的模式進一步確認海表風(fēng)場強迫是WJ的最主要驅(qū)動機制[25]。圖4給出ADCP位置處2013年春夏季緯向風(fēng)場與海表緯向流的時間演化結(jié)果。結(jié)果顯示在該站位附近,春季海表緯向風(fēng)場與海洋上層流場存在較強相關(guān),兩者相關(guān)系數(shù)為0.39,通過99%信度檢驗標(biāo)準(zhǔn)。伴隨海表風(fēng)場在4月底增強,WJ開始出現(xiàn),并逐漸發(fā)展增強。鑒于海表風(fēng)場是WJ產(chǎn)生的重要因素[26-27],那么我們提出一個假設(shè):2013年春季強WJ的主要原因是2013年春季海表風(fēng)場異常增強所導(dǎo)致。

    圖4 2013年觀測點處NCEP緯向風(fēng)場和ADCP觀測的緯向流時間序列圖Fig.4 Time series of zonal current and zonal surface wind in the year of 2013 at observation station

    3 數(shù)值模式模擬結(jié)果

    為理解2013年春季W(wǎng)J異常增強原因,我們通過海洋環(huán)流模式POM開展不同組別的數(shù)值實驗進行進一步分析。首先,我們利用氣候態(tài)風(fēng)場對POM進行強迫,開展控制實驗?zāi)M,檢驗POM模式對WJ事件的模擬能力。圖5為控制實驗結(jié)果。控制實驗中,在東北季風(fēng)影響下,2013-01—03,赤道印度洋被西向流控制。2013-04,西南季風(fēng)開始在西南印度洋建立,赤道印度洋上空逐漸被緯向西風(fēng)控制,春季W(wǎng)J開始在赤道印度洋出現(xiàn)。2013-05,赤道印度洋海表西風(fēng)達到最強,伴隨風(fēng)場的改變,春季W(wǎng)J達到盛期,最大值達到0.8 m/s。2013-06,隨著西南季風(fēng)北推,赤道印度洋上空西風(fēng)開始減弱,春季W(wǎng)J亦逐漸衰退以至消失。2013-07—09,赤道印度洋仍舊被弱西向流控制。2013-10,亞洲夏季風(fēng)開始衰退,東北季風(fēng)逐漸建立,此時在赤道印度洋上空重新出現(xiàn)了強盛的西風(fēng)環(huán)流,秋季W(wǎng)J開始出現(xiàn)。2013-11,東北季風(fēng)開始建立,秋季W(wǎng)J亦達到一年中的第二個盛期。2013-12,隨著赤道印度洋上空的西風(fēng)開始衰退,秋季W(wǎng)J開始消失。

    控制實驗表明,POM模式可以較好把握氣候態(tài)WJ事件,特別是春季W(wǎng)J的發(fā)生和發(fā)展變化規(guī)律。秋季W(wǎng)J雖然發(fā)生在2013-11,但春季W(wǎng)J比秋季W(wǎng)J強度略偏強。接下來,我們利用2013年風(fēng)場資料對POM模式進行強迫,來檢查POM模式對2013年WJ事件模擬情況。

    圖5 控制實驗?zāi)M的氣候態(tài)月平均印度洋緯向流(顏色填充)和NCEP風(fēng)場(矢量箭頭)Fig.5 Climatological monthly zonal surface current(shaded)simulated by control run and NCEP surface wind fields(arrows)

    圖6 是敏感性實驗得到的2013年印度洋海表月平均緯向流及海表風(fēng)場分布。在敏感性實驗中,2013-04,印度洋海表風(fēng)場由東北季風(fēng)開始轉(zhuǎn)變?yōu)槲髂霞撅L(fēng),赤道印度洋逐漸被西風(fēng)控制,春季W(wǎng)J開始出現(xiàn)。2013-05,赤道印度洋上空的西風(fēng)達到最強,同時春季W(wǎng)J亦達到最強,最大值達1.2 m/s,顯著強于控制實驗中的春季W(wǎng)J強度。2013-06,春季W(wǎng)J開始消失。2013-11,赤道印度洋上空重新被西風(fēng)所控制,秋季W(wǎng)J開始出現(xiàn)。2013-12,秋季W(wǎng)J達到最強,最強處流速達到0.8 m/s。從垂向結(jié)構(gòu)來看,敏感性實驗對WJ垂向結(jié)構(gòu)模擬也與觀測符合較好(圖2c):春季W(wǎng)J峰值出現(xiàn)在5月并且最大深度可以達到水面以下150 m。模式模擬得到的秋季W(wǎng)J主流區(qū)范圍也小于觀測結(jié)果,在100 m水深處流速已不足0.1 m/s。雖然POM模式對2013年秋季W(wǎng)J模擬存在一定偏差,但是對春季W(wǎng)J的模擬與觀測較符合。

    圖7給出敏感性實驗與控制實驗的差別。圖7表明,2013年春季赤道印度洋出現(xiàn)異常增強的緯向西風(fēng),2013-05的海表風(fēng)場強度比氣候態(tài)強約5 m/s,而敏感性實驗得到的赤道緯向流場比控制實驗強約0.45 m/s。以上結(jié)果表明2013年春季W(wǎng)J確實受到海表風(fēng)場影響而產(chǎn)生及發(fā)展。圖8為區(qū)域平均(60°~85°E, 2°S~2°N)的海表風(fēng)場及模擬得到的海表流場在敏感性實驗與控制實驗中的差別。圖8表明,模式模擬得到的海表流場隨著強迫風(fēng)場的改變而進行變化,二者具有較好相關(guān),風(fēng)場提前1周時兩者相關(guān)系數(shù)0.46,通過了99%顯著性檢驗。數(shù)值實驗結(jié)果證實WJ強度及變化與赤道印度洋風(fēng)場密切相關(guān)。

    圖6 敏感實驗?zāi)M的2013年印度洋緯向流(陰影填充)和NCEP風(fēng)場(矢量箭頭)Fig.6 Simulated zonal surface current(shaded)by sensitive run and NCEP surface wind fields(arrows)

    圖7 控制實驗和敏感性實驗之差(敏感性實驗減控制實驗)Fig.7 Differences between sensitive run and control run

    圖8 控制實驗和敏感性實驗之差在(60°~85°E,2°S~2°N)7 d滑動平均區(qū)域平均時間序列圖Fig.8 Variation of wind and current differences between sensitive run and control run at(60°~85°E,2°S~2°N) in a 7-day moving window

    4 討論和分析

    前人研究表明WJ強度受印度洋海表風(fēng)場控制。上一節(jié)的數(shù)值實驗結(jié)果也進一步證實2013年春季W(wǎng)J異常增強受赤道上空海表風(fēng)場影響。那么是什么原因?qū)е铝?013年赤道印度洋風(fēng)場異常增強?經(jīng)過分析, 2013年既非ENSO事件發(fā)生年份,也非IOD顯著發(fā)生年份,因此,可以排除這兩類太平洋、印度洋顯著年際時間尺度海氣相互作用對WJ的影響。

    熱帶印度洋是大氣季節(jié)內(nèi)震蕩(intraseasonal oscillations,ISO)發(fā)生和發(fā)展最完善的區(qū)域,ISO與印度洋地區(qū)的夏季風(fēng)暴發(fā)有著密切聯(lián)系[28-29]。而大氣低空西風(fēng)異常,是ISO對流位相的一個顯著特征,因此存在著2013年春季ISO影響赤道印度洋低空風(fēng)場的可能性。接下來,我們進一步檢查在2013-05,是否存在著ISO的對流位相。圖9給出了30~90 d帶通濾波后OLR與海表風(fēng)場的侯平均結(jié)果。圖9表明,2013-04中旬,ISO對流位相控制著熱帶中東印度洋,赤道印度洋海表西風(fēng)為正距平。同時,ISO抑制對流位相在西印度洋開始出現(xiàn),并逐漸東傳,于4月底傳到中東印度洋,赤道印度洋海表風(fēng)場亦呈現(xiàn)西風(fēng)負異常,該抑制對流位相產(chǎn)生的低空東風(fēng)距平不利于春季W(wǎng)J的產(chǎn)生與發(fā)展。與此同時,ISO對流位相在西南印度洋生成并于5月上旬控制熱帶中東印度洋,使得赤道印度洋低空被強大的西風(fēng)控制。因此,2013-05月初的ISO對流位相是5月赤道印度洋低空西風(fēng)異常的直接原因,并進一步使得2013年春季W(wǎng)J異常增強。

    以上分析可以發(fā)現(xiàn),2013-05在赤道印度洋出現(xiàn)了異常增強ISO對流位相,進而進一步激發(fā)了赤道印度洋上空的強烈西風(fēng)距平,從而誘使了春季W(wǎng)J的異常增強。ISO對WJ的影響不僅限于2013年,2011年開展的國際大型ISO聯(lián)合研究計劃SINDY/DYNAMO執(zhí)行期間,聯(lián)合觀測團隊同樣發(fā)現(xiàn)了類似現(xiàn)象,2011年秋季的一次ISO對流位相事件,同樣引發(fā)了赤道印度洋上層出現(xiàn)了強烈的東向流[30]。

    圖9 30~90 d帶通濾波后OLR距平(陰影填充)和緯向風(fēng)距平時空分布圖(2013-04-01—05-15,侯平均)Fig.9 Pented spatial distribution of 30~90 days band-pass filtered outgoing longwave radiation(shading) and zonal surface wind(contour)from April 1 to May 15,2013

    5 結(jié) 論

    利用“全球變化與海氣相互作用”專項布放的深海潛標(biāo)觀測資料及3套不同來源海洋再分析資料,本文對2013年赤道東印度洋WJ現(xiàn)象進行分析,發(fā)現(xiàn)2013年春季W(wǎng)J呈現(xiàn)異常變化特征,主要表現(xiàn)為春季W(wǎng)J強度略強于秋季W(wǎng)J,這與氣候態(tài)平均意義下的WJ年變化不同。進一步分析表明,2013年赤道印度洋海表風(fēng)場異常變化是導(dǎo)致春季W(wǎng)J異常的主要原因。而2013-05出現(xiàn)在赤道印度洋區(qū)域的ISO對流位相是赤道印度洋低空西風(fēng)距平產(chǎn)生的直接影響因子。最后,我們利用POM數(shù)值模式對2013年WJ現(xiàn)象進行了模擬,模擬結(jié)果可以反映上述特征,特別是2013年春季W(wǎng)J強于秋季W(wǎng)J,并進一步證實2013-05出現(xiàn)在赤道印度洋區(qū)域的ISO對流位相所誘導(dǎo)的西風(fēng)距平是春季W(wǎng)J變化的主要原因。本研究表明,作為季節(jié)內(nèi)時間尺度變化的ISO事件可以對WJ事件在年際時間尺度上進行調(diào)制,從一個方面反映了小尺度海氣相互作用過程對大尺度海洋環(huán)流事件的影響。

    [1] CHAO J P,YUAN S Y.Concerted development of atmosphere-ocean interaction events in the Tropical Indian Ocean and Pacific Ocean [J].Advances in Marine Science,2004,22(3):247-252.巢紀(jì)平,袁紹宇.熱帶印度洋和太平洋海氣相互作用事件的協(xié)調(diào)發(fā)展[J].海洋科學(xué)進展,2004,22(3):247-252.

    [2] LIU Y L,LüL G,SHU Q,et al.Seasonal variation of heat content and its influence factors in the upper northeastern Indian Ocean[J]. Advances in Marine Science,2009,27(4):434-443.劉延亮,呂連港,舒啟,等.東北印度洋上層海水熱含量及影響要素的季節(jié)變化[J].海洋科學(xué)進展,2009,27(4):434-443.

    [3] QIAO B,LIU Z Z,ZHANG S Y,et al.Equatorial current system structure and hydrologic characteristics in monsoonal wind transition period[J].Advances in Marine Science,2014,32(3):301-305.喬彬,劉子洲,張書穎,等.季風(fēng)轉(zhuǎn)換期東印度洋的赤道流系結(jié)構(gòu)和水文特征[J].海洋科學(xué)進展,2014,32(3):301-305.

    [4] QIU Y,LI L.Review of study on upper ocean circulation in the Bay of Bengal[J].Advances in Marine Science,2006,24(4):593-603.邱云,李立.孟加拉灣上層環(huán)流研究綜述[J].海洋科學(xué)進展,2006,24(4):593-603.

    [5] LIU L,YU W D,DIAO X Y,et al.Variability of Tropical Indian Ocean atmosphere during Tropical Indian Ocean Dipole mode event-The analysis of the stream function and the potential velocity of the wind field[J].Chinese Journal of Atmospheric Sciences,2008,32(5): 1803-1093.劉琳,于衛(wèi)東,刁新源,等.熱帶印度洋偶極子事件期間印度洋大氣的變化——大氣流函數(shù)場和勢函數(shù)場的分析[J].大氣科學(xué),2008,32(5):1803-1093.

    [6] LIU L,YU W D,LI T.Dynamic and thermodynamic air-sea coupling associated with the Indian Ocean Dipole diagnosed from 23 WCRP CMIP3 models[J].Journal of Climate,2011,24(18):4941-4958.

    [7] SU B,LIU L,LI K P,et al.The impacts of Equatorial Indian Ocean kelvin wave on the depth of thermocline in eastern coastal Andaman Sea[J].Advances in Marine Science,2013,31(3):343-350.蘇博,劉琳,李奎平,等.印度洋赤道Kelvin波對安達曼海東部近岸溫躍層深度的影響[J].海洋科學(xué)進展,2013,31(3):343-350.

    [8] YANG Y,YU W D,SUN J L,et al.Impact of Madden-Julian Oscillation on diurnal variation of sea surface temperature in the Tropical Indian Ocean[J].Advances in Marine Science,2013,31(3):351-359.楊洋,于衛(wèi)東,孫即霖,等.熱帶印度洋季節(jié)內(nèi)振蕩對海表溫度日變化的影響[J].海洋科學(xué)進展,2013,31(3):351-359.

    [9] LIU L,FENG L,YU W D,et al.The distribution and variability of simulated chlorophyll concentration over the Tropical Indian Ocean from five CMIP5 models[J].Journal of Ocean University of China,2013,12(2):253-259.

    [10] LIU H F,ZHANG X M,TANG Y M,et al.A Review of the studies on Indian Ocean Dipole and its predictability[J].Advances in Marine Science,2014,32(3):405-414.劉華鋒,章向明,唐佑民,等.印度洋偶極子及其可預(yù)報性研究進展[J].海洋科學(xué)進展,2014,32 (3):405-414.

    [11] WYRTKI K.An equatorial jet in the Indion Ocean[J].Science,1973,181(4096):262-264.

    [12] QIU Y,LI L,YU W.Behavior of the Wyrtki jet observed with surface drifting buoys and satellite altimeter[J].Geopgysical Research Letters,2009,36(18):120-131.

    [13] XUAN L L,QIU Y,XU J D,et al.Seasonal variation of surface-layer circulation in the eastern Tropical Indian Ocean[J].Journal of Tropical Oceanography,2014,33(1):26-35.宣莉莉,邱云,徐金電,等.熱帶東印度洋表層環(huán)流季節(jié)變化特征研究[J].熱帶海洋學(xué)報, 2014,33(1):26-35.

    [14] WANG Y,CUI F J.The structure and seasonal variation of upper-layer currents at central Equatorial Indian Ocean[J].Oceanologia et Limnologia Sinica,2015,46(2):241-247.王毅,崔鳳娟.赤道中印度洋上層環(huán)流結(jié)構(gòu)與季節(jié)變化特征分析[J].海洋與湖沼,2015,46 (2):241-247.

    [15] ZHANG Y,DU Y,ZHANG Y H,et al.Asymmetric influences of positive and negative IOD events on salinity transport by the fall Wyrtki jet along the Equatorial Indian Ocean[J].Journal of Tropical Oceanography,2015,34(5):1-10.張瑩,杜巖,張玉紅,等.赤道印度洋秋季W(wǎng)yrtki急流鹽度輸運對正、負IOD事件的不對稱響應(yīng)[J].熱帶海洋學(xué)報,2015,34(5):1-10.

    [16] ZHANG Y H,DU Y,XU H M.Seasonal variability of the west-east water mass exchange on the section of central Equatorial Indian Ocean and its regional difference[J].Haiyang Xuebao,2012,34(2):30-38.張玉紅,杜巖,徐海明.赤道印度洋中部斷面東西水交換的季節(jié)變化及其區(qū)域差異[J].海洋學(xué)報,2012,34(2):30-38.

    [17] SREENIVAS P,CHOWDARY J S,GNANASEELAN C.Impact of tropical cyclones on the intensity and phase propagation of fall Wyrtki jets[J].Geopgysical Research Letters,2012,39(22):L22603.

    [18] NYADJRO E S,MCPHADEN M J.Variability of zonal currents in the eastern Equatorial Indian Ocean on seasonal to interannual time scales[J].Journal of Geophysical Research Oceans,2014,119(11):7969-7986.

    [19] CHU P C.Observational studies on association between eastward equatorial jet and Indian Ocean Dipole[J].Journal of Oceanography, 2010,66(3):429-434.

    [20] GNANASEELAN C,ATIDI D,MCPHADEN M J.Impact of Indian Ocean Dipole and El Ni?o/Southern Oscillation wind-forcing on the Wyrtki jets[J].Journal of Geophysical Research Oceans,2012,117(C8):C08005.

    [21] MCPHADEN,M J,WANG Y,RAVICHANDRAN M.Volume transports of the Wyrtki jets and their relationship to the Indian Ocean Dipole[J].Journal of Geopgysical Research,2015,120(8):5302-5317.

    [22] NAGURA M,MCPHADEN M J.Wyrtki jet dynamics:seasonal variability[J].Journal of Geophysical Research Oceans,2010a,117 (C7):319-336.

    [23] NAGURA M,MCPHADEN M J.Dynamics of zonal current variations associated with the Indian Ocean Dipole[J].Journal of Geophysical Research Oceans,2010b,115(C11):2027-2042.

    [24] JOSEPH S,WALLCRAFT A J,JENSEN T G,et al.Weakening of spring Wyrtki jets in the Indian Ocean during 2006-2011[J].Journal of Geopgysical Research,2012,117(C4):76-85.

    [25] HAN W,MCCREARY J P,ANDERSON D L T,et al.Dynamics of the eastern surface jets in the Equatorial Indian Ocean[J].Journal of Physical Oceanography,1999,29(9):2191-2209.

    [26] NAGURA M,MCPHADEN M J.The dynamics of zonal current variations in the central Equatorial Indian Ocean[J].Geophysical Research Letters,2008,35(35):186-203.

    [27] GRODSKY S A,CARTON J A,MURTUGUDDE R.Anomalous surface currents in the Tropical Indian Ocean[J].Geophysical Research Letters,2001,28(22):4207-4210.

    [28] LI K P,YU W D,LI L L.Physical structure of Madden-Julian oscillation in the Tropical Indian Ocean during winter[J].Advances in Marine Science,2009,27(2):133-139.李奎平,于衛(wèi)東,李俐俐.冬季熱帶印度洋季節(jié)內(nèi)振蕩的物理結(jié)構(gòu)[J].海洋科學(xué)進展,2009,27 (2):133-139.

    [29] LI C Y,PAN J,SONG J.Progress on the MJO research in recent years[J].Chinese Journal of Atmospheric Sciences,2013,37(2):229-252.李崇銀,潘靜,宋潔.MJO研究新進展[J].大氣科學(xué),2013,37(2):229-252.

    [30] SHINODA T,JENSEN T G,FLATAU M,et al.Large-scale oceanic variability associated with the Madden-Julian oscillation during the CINDY/DYNAMO field campaign from satellite observations[J].Romote Sensing,2013,5(5):2072-2092.

    Anomalous Behavior of Spring Wyrtki Jet in Equatorial Indian Ocean During 2013

    HAN Guo-qing1,LIU Lin1,2,DUAN Yong-liang1,2,WANG Guan-suo2,3,WANG Hui-wu1,2,
    LIU Yan-liang1,2,WANG Hai-yuan1,FENG Lin1,2,YU Wei-dong1,2(1.Center for Ocean and Climate Research,The First Institute of Oceanography,SOA,Qingdao 266061,China; 2.Laboratory for Regional Oceanography and Numerical Modeling,Qingdao National Laboratory for Marine Science and Technology,Qingdao 266200,China; 3.Key Laboratory of Marine Science and Numerical Modeling,The First Institute of Oceanography, SOA,Qingdao 266061,China)

    In-situ measurements of the upper ocean currents reveal remarkable abnormal behavior of Wyrtki Jet in boreal spring in tropical Indian Ocean in 2013.The Wyrtki Jet in boreal spring was unusually stronger than its counterpart in fall,clearly against previous understanding,and this phenomenon is also demonstrated by reanalysis data.Further analysis and numerical experiments show that the anomalously enhanced Wyrtki Jet is related to the surface wind anomaly in equatorial Indian Ocean,which is resulted from the strong intra-seasonal oscillation(ISO)event in that season.This study show that mesoscale airsea interaction can influence relatively large scale ocean current on interannual timescale.

    Wyrtki Jet;intra-seasonal oscillation;Indian Ocean;interannual variability

    March 29,2016

    P732

    A

    1671-6647(2017)02-0189-11

    10.3969/j.issn.1671-6647.2017.02.004

    2016-03-29

    全球變化研究國家重大科學(xué)研究計劃項目——太平洋印度洋對全球變暖的響應(yīng)及其對氣候變化的調(diào)控作用(2012CB955601);國家自然科學(xué)基金項目——全球變暖背景下印度洋年際時間尺度海氣相互作用對季節(jié)內(nèi)振蕩的影響(41376037);國家自然科學(xué)基金委員會項目-山東省人民政府聯(lián)合資助海洋科學(xué)研究中心項目——海洋環(huán)境動力學(xué)和數(shù)值模擬(U1606405);全球變化與海氣相互作用專項項目——季風(fēng)變化對印度洋東部和南海南部上升流生態(tài)系統(tǒng)的影響(GASI-03-01-03-03),熱帶印度洋海洋過程與海氣相互作用(GASIIPOVAI-02),太平洋-印度洋與亞洲季風(fēng)的相互影響(GASI-IPOVAI-03)和東印度洋南部水體綜合調(diào)查春、秋航次(GASI-02-IND-STSspr,GASI-02-IND-STSaut)

    韓國慶(1990-),男,山東濰坊人,碩士研究生,主要從事海氣相互作用方面研究.E-mail:gqhan@fio.org.cn

    *通訊作者:劉 琳(1978-),男,山西陽泉人,副研究員,博士,主要從事海氣相互作用方面研究.E-mail:liul@fio.org.cn

    猜你喜歡
    海表赤道西風(fēng)
    “雪龍2”號過赤道,救下一船人!
    軍事文摘(2024年2期)2024-01-10 01:58:44
    古道西風(fēng)(水彩)
    基于無人機的海表環(huán)境智能監(jiān)測系統(tǒng)設(shè)計與應(yīng)用
    2016與1998年春季北大西洋海表溫度異常的差異及成因
    距赤道越遠越喜愛黃色
    融合海表溫度產(chǎn)品在渤黃東海的對比分析及初步驗證
    海洋通報(2020年6期)2020-03-19 02:10:08
    秋漸濃,分付西風(fēng)此夜涼
    文苑(2019年20期)2019-11-16 08:52:46
    浣溪沙 誰念西風(fēng)獨自涼
    飛魔幻A(2019年3期)2019-06-12 09:22:35
    太陽總輻照度對熱帶中太平洋海表溫度年代際變化的可能影響
    最熱的地方不在赤道
    亚洲在线自拍视频| 亚洲精品色激情综合| 看十八女毛片水多多多| 尤物成人国产欧美一区二区三区| 国产精品嫩草影院av在线观看| 亚洲精品色激情综合| 国产精品久久久久久亚洲av鲁大| 黄色欧美视频在线观看| 男插女下体视频免费在线播放| 免费高清视频大片| 在线观看免费视频日本深夜| 中文在线观看免费www的网站| 99久久精品国产国产毛片| 丰满乱子伦码专区| 深夜a级毛片| av天堂中文字幕网| 欧美一级a爱片免费观看看| 欧美色视频一区免费| 亚洲自拍偷在线| 精品久久久噜噜| 国产v大片淫在线免费观看| 亚洲高清免费不卡视频| 黑人高潮一二区| 亚洲精品亚洲一区二区| 亚洲精华国产精华液的使用体验 | 国产精品一区二区三区四区久久| 亚洲精品在线观看二区| 国产伦在线观看视频一区| 直男gayav资源| 一级黄片播放器| 亚洲五月天丁香| 亚洲av美国av| 最好的美女福利视频网| 成年版毛片免费区| 99热精品在线国产| 性欧美人与动物交配| 久久久色成人| 色播亚洲综合网| 久久精品综合一区二区三区| 色吧在线观看| 国内揄拍国产精品人妻在线| 国产高清视频在线播放一区| 中文字幕av在线有码专区| 人妻久久中文字幕网| 国内精品宾馆在线| 亚洲七黄色美女视频| 国产一区亚洲一区在线观看| 久久久午夜欧美精品| 久久久久国产网址| 99精品在免费线老司机午夜| 亚洲一区高清亚洲精品| 99riav亚洲国产免费| 国产精品女同一区二区软件| 国产精品久久久久久精品电影| 日本免费一区二区三区高清不卡| 欧美不卡视频在线免费观看| 亚洲精华国产精华液的使用体验 | 免费大片18禁| 精品99又大又爽又粗少妇毛片| 成人漫画全彩无遮挡| 天堂√8在线中文| 丰满的人妻完整版| 亚洲成人中文字幕在线播放| 精品人妻偷拍中文字幕| 99久久九九国产精品国产免费| 成年版毛片免费区| 亚洲在线观看片| 久久草成人影院| 99热这里只有是精品在线观看| 日本欧美国产在线视频| 欧美zozozo另类| 听说在线观看完整版免费高清| 在线免费观看的www视频| 一区福利在线观看| 免费看a级黄色片| 一区二区三区高清视频在线| .国产精品久久| 中国美白少妇内射xxxbb| 又爽又黄无遮挡网站| 久久久久九九精品影院| 男女视频在线观看网站免费| 国产一区二区三区在线臀色熟女| 国产精品一区二区三区四区久久| 国产蜜桃级精品一区二区三区| 日本色播在线视频| 欧美一区二区精品小视频在线| av天堂中文字幕网| 国产伦在线观看视频一区| 日韩大尺度精品在线看网址| 日韩欧美国产在线观看| 欧美高清性xxxxhd video| 干丝袜人妻中文字幕| 精品无人区乱码1区二区| 精品欧美国产一区二区三| 一边摸一边抽搐一进一小说| 男女之事视频高清在线观看| 欧美激情在线99| 免费观看人在逋| 少妇裸体淫交视频免费看高清| 黄色视频,在线免费观看| 99在线人妻在线中文字幕| 美女黄网站色视频| 日本 av在线| 一区福利在线观看| 国产精品av视频在线免费观看| 少妇人妻一区二区三区视频| 亚洲美女黄片视频| 一区二区三区高清视频在线| 亚洲中文字幕一区二区三区有码在线看| 欧美区成人在线视频| 免费av不卡在线播放| 少妇熟女aⅴ在线视频| 国产成人一区二区在线| 日韩在线高清观看一区二区三区| 少妇的逼水好多| 成人一区二区视频在线观看| 青春草视频在线免费观看| 色尼玛亚洲综合影院| 国产av一区在线观看免费| 麻豆av噜噜一区二区三区| 亚洲国产色片| 1024手机看黄色片| 久久午夜亚洲精品久久| 最近在线观看免费完整版| 久久午夜亚洲精品久久| 亚洲欧美日韩高清专用| 美女高潮的动态| 欧美色欧美亚洲另类二区| 午夜爱爱视频在线播放| 午夜精品在线福利| 亚洲精品国产av成人精品 | 成人美女网站在线观看视频| 极品教师在线视频| 国产中年淑女户外野战色| 精品无人区乱码1区二区| 国产欧美日韩一区二区精品| 亚洲精品久久国产高清桃花| 亚洲欧美日韩卡通动漫| 嫩草影院新地址| av在线天堂中文字幕| 中国国产av一级| 国产精品女同一区二区软件| 一个人免费在线观看电影| 欧美性感艳星| 亚洲国产精品成人综合色| 91久久精品国产一区二区成人| 亚洲国产精品成人综合色| 嫩草影院入口| 性色avwww在线观看| 给我免费播放毛片高清在线观看| 俺也久久电影网| 国产成人影院久久av| 久久午夜亚洲精品久久| 亚洲va在线va天堂va国产| 久久九九热精品免费| 啦啦啦观看免费观看视频高清| av黄色大香蕉| 亚洲在线观看片| 中文字幕av在线有码专区| 久久久久久久久久成人| 天天一区二区日本电影三级| 三级经典国产精品| 又粗又爽又猛毛片免费看| av免费在线看不卡| 麻豆国产97在线/欧美| 舔av片在线| 狠狠狠狠99中文字幕| 男女视频在线观看网站免费| 日本一二三区视频观看| 国产精品美女特级片免费视频播放器| 99热全是精品| 老司机福利观看| 日本一本二区三区精品| 赤兔流量卡办理| 日韩,欧美,国产一区二区三区 | 久久精品综合一区二区三区| 偷拍熟女少妇极品色| 搡老妇女老女人老熟妇| 国产成人a∨麻豆精品| 久久午夜亚洲精品久久| 最近的中文字幕免费完整| 在线看三级毛片| 午夜精品国产一区二区电影 | 日韩精品中文字幕看吧| av在线老鸭窝| 成人漫画全彩无遮挡| 日本爱情动作片www.在线观看 | 亚洲中文字幕日韩| 国语自产精品视频在线第100页| 三级经典国产精品| 亚洲一区二区三区色噜噜| 又粗又爽又猛毛片免费看| 97人妻精品一区二区三区麻豆| 久久亚洲精品不卡| 深夜精品福利| 国产一区二区激情短视频| 最近视频中文字幕2019在线8| 乱人视频在线观看| 免费大片18禁| 亚洲国产精品成人综合色| 亚洲精品国产av成人精品 | 又黄又爽又免费观看的视频| 伦理电影大哥的女人| 欧美三级亚洲精品| 国产av在哪里看| 午夜a级毛片| АⅤ资源中文在线天堂| 久久久久久国产a免费观看| 国产av不卡久久| 最近视频中文字幕2019在线8| 男女之事视频高清在线观看| 女人十人毛片免费观看3o分钟| 亚洲av熟女| 亚洲精品久久国产高清桃花| 欧美xxxx黑人xx丫x性爽| 日韩亚洲欧美综合| 成年免费大片在线观看| 成人亚洲精品av一区二区| 特大巨黑吊av在线直播| 久久草成人影院| 特大巨黑吊av在线直播| 色av中文字幕| 亚洲天堂国产精品一区在线| 精品久久久久久久久久久久久| 久久久久精品国产欧美久久久| 级片在线观看| 亚洲在线观看片| 精品久久久噜噜| 成人精品一区二区免费| 欧美日韩国产亚洲二区| 看免费成人av毛片| 亚洲成人av在线免费| 精品熟女少妇av免费看| 两个人的视频大全免费| 亚洲精品日韩在线中文字幕 | 天堂av国产一区二区熟女人妻| 黄色欧美视频在线观看| 真人做人爱边吃奶动态| 午夜精品一区二区三区免费看| 国产精华一区二区三区| av在线播放精品| 久久人人爽人人爽人人片va| 久久久精品欧美日韩精品| 成人午夜高清在线视频| 一区二区三区免费毛片| 亚洲人与动物交配视频| 91精品国产九色| 国产aⅴ精品一区二区三区波| 麻豆一二三区av精品| 欧美激情久久久久久爽电影| av在线蜜桃| 波多野结衣高清无吗| 国产av在哪里看| 九九久久精品国产亚洲av麻豆| 免费高清视频大片| 18禁在线播放成人免费| 夜夜爽天天搞| 日韩亚洲欧美综合| 最近最新中文字幕大全电影3| 日韩成人伦理影院| 插逼视频在线观看| 两个人的视频大全免费| 亚洲欧美日韩无卡精品| 久久久久久伊人网av| 国产高清视频在线观看网站| 最近中文字幕高清免费大全6| 久久午夜福利片| 国产麻豆成人av免费视频| 中文字幕免费在线视频6| av在线老鸭窝| 国产伦一二天堂av在线观看| 国产精品久久久久久久电影| 国产女主播在线喷水免费视频网站 | 麻豆久久精品国产亚洲av| 国产三级中文精品| 婷婷六月久久综合丁香| 亚洲最大成人av| 免费观看的影片在线观看| 国产69精品久久久久777片| 亚洲第一电影网av| 成人永久免费在线观看视频| 一区二区三区高清视频在线| avwww免费| 听说在线观看完整版免费高清| 最近手机中文字幕大全| 男人狂女人下面高潮的视频| 国内精品久久久久精免费| 亚洲国产精品合色在线| 人妻夜夜爽99麻豆av| 九九久久精品国产亚洲av麻豆| 我要看日韩黄色一级片| 男人的好看免费观看在线视频| 啦啦啦啦在线视频资源| 亚洲三级黄色毛片| 欧美最新免费一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 老熟妇仑乱视频hdxx| 亚洲综合色惰| 我要搜黄色片| 干丝袜人妻中文字幕| 免费观看精品视频网站| 大香蕉久久网| 自拍偷自拍亚洲精品老妇| 十八禁网站免费在线| 久久韩国三级中文字幕| 久久久国产成人精品二区| 99精品在免费线老司机午夜| 插逼视频在线观看| 少妇被粗大猛烈的视频| 噜噜噜噜噜久久久久久91| 春色校园在线视频观看| 亚洲av不卡在线观看| 久久久久久久久中文| 麻豆久久精品国产亚洲av| 日韩欧美一区二区三区在线观看| 综合色丁香网| 欧美xxxx性猛交bbbb| 精品久久久噜噜| 国产精品久久久久久久电影| 男女下面进入的视频免费午夜| a级毛片a级免费在线| 最近手机中文字幕大全| 亚洲av不卡在线观看| a级一级毛片免费在线观看| 永久网站在线| 国产精品国产三级国产av玫瑰| av视频在线观看入口| 久久久久性生活片| 亚洲av成人精品一区久久| 国产白丝娇喘喷水9色精品| 深夜a级毛片| 免费看美女性在线毛片视频| 简卡轻食公司| 精品一区二区三区av网在线观看| 国产熟女欧美一区二区| 男人舔女人下体高潮全视频| 亚洲性夜色夜夜综合| 少妇人妻精品综合一区二区 | 男女边吃奶边做爰视频| 嫩草影院新地址| 97超级碰碰碰精品色视频在线观看| 国产在线男女| 一区二区三区高清视频在线| 亚洲一级一片aⅴ在线观看| 久久韩国三级中文字幕| 欧美3d第一页| 国产亚洲精品av在线| 中文字幕av成人在线电影| 久久中文看片网| 性插视频无遮挡在线免费观看| 成人鲁丝片一二三区免费| 联通29元200g的流量卡| av天堂中文字幕网| 少妇丰满av| 欧美一区二区国产精品久久精品| 免费观看的影片在线观看| 色尼玛亚洲综合影院| 久久久久国产精品人妻aⅴ院| 一进一出抽搐gif免费好疼| 搡女人真爽免费视频火全软件 | 99久国产av精品| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久| 日韩欧美三级三区| 麻豆国产av国片精品| 天堂网av新在线| 我的老师免费观看完整版| 日韩欧美三级三区| videossex国产| 亚洲国产精品久久男人天堂| 国产老妇女一区| 丰满乱子伦码专区| 99久国产av精品国产电影| 亚洲久久久久久中文字幕| 色5月婷婷丁香| 婷婷精品国产亚洲av在线| 乱人视频在线观看| 久久精品国产亚洲网站| 日韩一区二区视频免费看| 两个人的视频大全免费| www日本黄色视频网| 嫩草影院精品99| 少妇熟女欧美另类| 国产精品亚洲美女久久久| 国产一区二区三区av在线 | 日韩精品有码人妻一区| 舔av片在线| 美女高潮的动态| 国产成人freesex在线 | 中文字幕精品亚洲无线码一区| 亚洲丝袜综合中文字幕| 你懂的网址亚洲精品在线观看 | 51国产日韩欧美| 一级毛片aaaaaa免费看小| 欧美一区二区精品小视频在线| 麻豆国产97在线/欧美| 国产精品伦人一区二区| 日韩,欧美,国产一区二区三区 | 狂野欧美白嫩少妇大欣赏| 天堂av国产一区二区熟女人妻| 波多野结衣高清无吗| 九九热线精品视视频播放| 少妇人妻精品综合一区二区 | 国产片特级美女逼逼视频| 国产av麻豆久久久久久久| 日韩精品青青久久久久久| 嫩草影院精品99| 99久国产av精品| 女同久久另类99精品国产91| 亚洲自偷自拍三级| 亚洲最大成人手机在线| 精品一区二区三区av网在线观看| 不卡一级毛片| 变态另类丝袜制服| 别揉我奶头~嗯~啊~动态视频| 人妻丰满熟妇av一区二区三区| 人妻夜夜爽99麻豆av| 色哟哟·www| 国产精品一二三区在线看| 国产午夜精品论理片| 国产一区二区三区av在线 | 国产精品一区二区三区四区久久| 男女视频在线观看网站免费| 亚洲图色成人| 99视频精品全部免费 在线| 99久久中文字幕三级久久日本| 男人和女人高潮做爰伦理| 99热精品在线国产| 国内少妇人妻偷人精品xxx网站| 黑人高潮一二区| 最近中文字幕高清免费大全6| 男女之事视频高清在线观看| 亚洲aⅴ乱码一区二区在线播放| 综合色av麻豆| 欧美一区二区国产精品久久精品| 最近视频中文字幕2019在线8| 一级毛片我不卡| 国产色爽女视频免费观看| 我的老师免费观看完整版| a级毛片免费高清观看在线播放| 男人狂女人下面高潮的视频| 色综合色国产| 亚洲国产欧洲综合997久久,| 桃色一区二区三区在线观看| 欧美不卡视频在线免费观看| 亚洲av免费在线观看| 国产中年淑女户外野战色| 日本熟妇午夜| 亚州av有码| 亚洲成a人片在线一区二区| 国产精品国产高清国产av| 精品乱码久久久久久99久播| 可以在线观看的亚洲视频| 免费看光身美女| 久久久成人免费电影| 别揉我奶头~嗯~啊~动态视频| 熟妇人妻久久中文字幕3abv| АⅤ资源中文在线天堂| 亚洲欧美清纯卡通| 中国美白少妇内射xxxbb| 热99在线观看视频| 精品人妻熟女av久视频| 少妇熟女aⅴ在线视频| 国产高清视频在线观看网站| 欧美丝袜亚洲另类| 日韩精品中文字幕看吧| 国产成年人精品一区二区| 国产精品一二三区在线看| 在线观看免费视频日本深夜| 一级毛片电影观看 | 黄色配什么色好看| 一进一出抽搐gif免费好疼| 国产精品一区二区三区四区久久| 中文字幕精品亚洲无线码一区| 人人妻人人澡人人爽人人夜夜 | 在线国产一区二区在线| 久久这里只有精品中国| 韩国av在线不卡| 亚洲中文字幕一区二区三区有码在线看| 国产精品爽爽va在线观看网站| 蜜桃久久精品国产亚洲av| 日本黄大片高清| 俄罗斯特黄特色一大片| 蜜桃久久精品国产亚洲av| 两个人的视频大全免费| 99国产精品一区二区蜜桃av| 成人精品一区二区免费| 91麻豆精品激情在线观看国产| 国产亚洲欧美98| 色哟哟哟哟哟哟| 亚洲婷婷狠狠爱综合网| 麻豆一二三区av精品| 亚洲一级一片aⅴ在线观看| 久久草成人影院| 一区二区三区高清视频在线| 一个人观看的视频www高清免费观看| 美女被艹到高潮喷水动态| 菩萨蛮人人尽说江南好唐韦庄 | 久久久欧美国产精品| 亚洲欧美日韩无卡精品| 在线观看午夜福利视频| 最新中文字幕久久久久| 亚洲国产欧洲综合997久久,| 六月丁香七月| 中文在线观看免费www的网站| 免费观看在线日韩| 精品久久国产蜜桃| 在线播放国产精品三级| 精品一区二区三区视频在线观看免费| 久久人人爽人人爽人人片va| 成人鲁丝片一二三区免费| 黄色一级大片看看| 99热网站在线观看| 在线观看午夜福利视频| 女人被狂操c到高潮| 国产伦精品一区二区三区四那| 内射极品少妇av片p| 亚洲久久久久久中文字幕| 国产成人freesex在线 | 免费看a级黄色片| 中文字幕熟女人妻在线| 久久久国产成人精品二区| 天天躁夜夜躁狠狠久久av| 日韩av在线大香蕉| av在线老鸭窝| 精华霜和精华液先用哪个| 一边摸一边抽搐一进一小说| 国产精品久久电影中文字幕| 久久久久国产网址| 亚洲熟妇中文字幕五十中出| 欧美xxxx性猛交bbbb| 国产精品亚洲美女久久久| 免费av毛片视频| 免费在线观看成人毛片| 亚洲国产高清在线一区二区三| 亚洲最大成人av| 色尼玛亚洲综合影院| 一进一出抽搐gif免费好疼| 一夜夜www| 免费观看精品视频网站| 日本黄大片高清| 亚洲无线在线观看| 九九爱精品视频在线观看| 亚洲在线观看片| 日本黄色视频三级网站网址| av在线播放精品| 精品久久久噜噜| 精品人妻视频免费看| 美女内射精品一级片tv| 日本熟妇午夜| 丰满人妻一区二区三区视频av| 欧美色视频一区免费| 伊人久久精品亚洲午夜| 久久久精品94久久精品| 干丝袜人妻中文字幕| .国产精品久久| 国产一区二区在线av高清观看| 成人亚洲精品av一区二区| 最好的美女福利视频网| 99riav亚洲国产免费| 国产精品不卡视频一区二区| 国产精品久久久久久精品电影| 中出人妻视频一区二区| 婷婷精品国产亚洲av| 色哟哟·www| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美 国产精品| 亚洲欧美清纯卡通| 亚洲人成网站在线观看播放| 久久久久久伊人网av| 99久久无色码亚洲精品果冻| 国内精品一区二区在线观看| 国产精品一区二区三区四区久久| 亚洲成a人片在线一区二区| 国产熟女欧美一区二区| 精品久久久久久久末码| 国产av麻豆久久久久久久| 国产 一区精品| 日韩高清综合在线| 亚洲欧美日韩高清在线视频| 97人妻精品一区二区三区麻豆| 午夜免费激情av| 尾随美女入室| 日韩欧美一区二区三区在线观看| 成人亚洲精品av一区二区| 国产高清不卡午夜福利| 国产黄a三级三级三级人| 免费观看精品视频网站| 日本与韩国留学比较| 精品久久久久久久末码| 日韩一区二区视频免费看| 男女边吃奶边做爰视频| 在线播放国产精品三级| 婷婷六月久久综合丁香| 国产精品一区www在线观看| 全区人妻精品视频| 日本一二三区视频观看| 看免费成人av毛片| а√天堂www在线а√下载| 人人妻人人看人人澡| 中国国产av一级| 啦啦啦观看免费观看视频高清| 少妇被粗大猛烈的视频| 欧美另类亚洲清纯唯美| 欧美又色又爽又黄视频| 欧美xxxx性猛交bbbb| 免费观看在线日韩| 亚洲欧美日韩高清专用| 免费av不卡在线播放| 精品久久国产蜜桃| 午夜免费男女啪啪视频观看 | aaaaa片日本免费| 婷婷精品国产亚洲av| 深爱激情五月婷婷| 亚洲人成网站高清观看|