• <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天堂| 午夜av观看不卡| 国产色婷婷99| 午夜激情福利司机影院| 国产视频首页在线观看| 日本色播在线视频| 欧美日韩综合久久久久久| 国产一区有黄有色的免费视频| 国产亚洲精品久久久com| 中文字幕最新亚洲高清| 哪个播放器可以免费观看大片| 女性生殖器流出的白浆| 99热这里只有精品一区| 免费观看的影片在线观看| 国产男人的电影天堂91| 国产一区二区在线观看日韩| 精品熟女少妇av免费看| 国产av一区二区精品久久| 狂野欧美白嫩少妇大欣赏| 天堂8中文在线网| 国产综合精华液| 女人精品久久久久毛片| 免费高清在线观看视频在线观看| 国产精品久久久久久久电影| 十八禁网站网址无遮挡| 乱码一卡2卡4卡精品| 性色av一级| 高清欧美精品videossex| av天堂久久9| 日本wwww免费看| 天天躁夜夜躁狠狠久久av| 久久久国产精品麻豆| 免费少妇av软件| 嘟嘟电影网在线观看| 精品久久久精品久久久| 久久ye,这里只有精品| 国产欧美日韩综合在线一区二区| 中文字幕人妻丝袜制服| 欧美激情国产日韩精品一区| 国产日韩欧美在线精品| 国产视频内射| 中文字幕精品免费在线观看视频 | 午夜福利在线观看免费完整高清在| 日韩欧美一区视频在线观看| 高清在线视频一区二区三区| 如何舔出高潮| 一区二区三区乱码不卡18| 国产精品女同一区二区软件| 成人国产麻豆网| 日韩av在线免费看完整版不卡| 免费av不卡在线播放| 国产伦理片在线播放av一区| 亚洲精品国产av蜜桃| 亚洲国产最新在线播放| 国产精品欧美亚洲77777| 精品国产乱码久久久久久小说| 亚洲av日韩在线播放| a级片在线免费高清观看视频| av黄色大香蕉| 国产一区有黄有色的免费视频| 老司机亚洲免费影院| 免费大片18禁| 欧美另类一区| 少妇被粗大猛烈的视频| 亚洲av国产av综合av卡| 最近中文字幕2019免费版| 国产探花极品一区二区| 人人妻人人爽人人添夜夜欢视频| 少妇被粗大的猛进出69影院 | 亚洲欧洲国产日韩| 免费观看在线日韩| 一个人看视频在线观看www免费| 一级爰片在线观看| 久久久久久久国产电影| 少妇 在线观看| 美女cb高潮喷水在线观看| 久久久久久久精品精品| 少妇的逼好多水| 国产精品秋霞免费鲁丝片| 日韩av不卡免费在线播放| 国产精品一区二区在线观看99| 午夜免费男女啪啪视频观看| 91成人精品电影| 国产日韩欧美亚洲二区| 久久人人爽人人片av| 亚洲四区av| 性色avwww在线观看| 国产成人a∨麻豆精品| 观看美女的网站| 国产亚洲最大av| 少妇的逼好多水| 亚洲一级一片aⅴ在线观看| 视频区图区小说| 麻豆成人av视频| 熟妇人妻不卡中文字幕| 91精品国产九色| 一本一本综合久久| 91精品国产国语对白视频| 人人澡人人妻人| 有码 亚洲区| 精品国产一区二区久久| 国产黄频视频在线观看| 日产精品乱码卡一卡2卡三| 天堂中文最新版在线下载| 日本黄色日本黄色录像| 成人亚洲精品一区在线观看| 人妻制服诱惑在线中文字幕| 亚洲欧美成人精品一区二区| 久久婷婷青草| 十八禁高潮呻吟视频| 香蕉精品网在线| 在线观看www视频免费| 简卡轻食公司| 美女cb高潮喷水在线观看| av卡一久久| 国产高清有码在线观看视频| 制服人妻中文乱码| 免费大片18禁| 国国产精品蜜臀av免费| 精品久久国产蜜桃| 少妇人妻久久综合中文| 国产精品无大码| 蜜臀久久99精品久久宅男| 成人毛片60女人毛片免费| 欧美xxⅹ黑人| 久久99一区二区三区| 亚洲欧美精品自产自拍| 22中文网久久字幕| 视频在线观看一区二区三区| 少妇被粗大猛烈的视频| 免费不卡的大黄色大毛片视频在线观看| 男女高潮啪啪啪动态图| 我要看黄色一级片免费的| 99热全是精品| 久久精品久久久久久久性| 国产精品一二三区在线看| 久久精品夜色国产| 欧美xxxx性猛交bbbb| 亚洲国产精品成人久久小说| 丝袜美足系列| 国产成人a∨麻豆精品| 啦啦啦中文免费视频观看日本| 黄色怎么调成土黄色| 亚洲国产精品一区三区| 国产亚洲最大av| 久久久久人妻精品一区果冻| √禁漫天堂资源中文www| 黑人高潮一二区| 中国三级夫妇交换| 亚洲国产av新网站| 久久久久精品久久久久真实原创| 国产精品无大码| 人成视频在线观看免费观看| 中文字幕人妻丝袜制服| 成人免费观看视频高清| 久久久久久久久久久免费av| 麻豆成人av视频| 色5月婷婷丁香| 少妇的逼水好多| 精品人妻熟女毛片av久久网站| 国产av码专区亚洲av| 人妻夜夜爽99麻豆av| 成人国产麻豆网| 777米奇影视久久| 丝袜在线中文字幕| 大香蕉97超碰在线| 国产精品久久久久久精品电影小说| 国产老妇伦熟女老妇高清| 中文乱码字字幕精品一区二区三区| 在现免费观看毛片| 精品亚洲成国产av| 精品酒店卫生间| 国产精品成人在线| 免费久久久久久久精品成人欧美视频 | 高清在线视频一区二区三区| 精品酒店卫生间| 成人黄色视频免费在线看| 在线观看国产h片| 欧美xxxx性猛交bbbb| 久久精品熟女亚洲av麻豆精品| 天堂中文最新版在线下载| 一级毛片 在线播放| 少妇猛男粗大的猛烈进出视频| 嘟嘟电影网在线观看| 高清视频免费观看一区二区| 国产精品无大码| 成人漫画全彩无遮挡| 在线亚洲精品国产二区图片欧美 | 亚洲精品成人av观看孕妇| 国产精品一区二区三区四区免费观看| 久久久久视频综合| 成年av动漫网址| 男女免费视频国产| 丝袜美足系列| 国产精品国产av在线观看| 最黄视频免费看| 亚洲精品日本国产第一区| 久久精品国产亚洲av涩爱| 爱豆传媒免费全集在线观看| 亚洲精品中文字幕在线视频| 天天躁夜夜躁狠狠久久av| 卡戴珊不雅视频在线播放| 国产免费又黄又爽又色| 国产一区二区在线观看av| 午夜激情av网站| 女人久久www免费人成看片| 午夜久久久在线观看| 特大巨黑吊av在线直播| 亚洲av国产av综合av卡| 亚洲精品国产av成人精品| 日韩成人伦理影院| 99热网站在线观看| 嘟嘟电影网在线观看| 三级国产精品欧美在线观看| 毛片一级片免费看久久久久| 蜜桃在线观看..| 七月丁香在线播放| 国产精品蜜桃在线观看| 蜜臀久久99精品久久宅男| 午夜福利影视在线免费观看| 女人久久www免费人成看片| 有码 亚洲区| 高清毛片免费看| 天天躁夜夜躁狠狠久久av| 18禁在线无遮挡免费观看视频| 男女无遮挡免费网站观看| 看非洲黑人一级黄片| 亚洲国产精品成人久久小说| 熟女电影av网| 国产精品久久久久久av不卡| 伦精品一区二区三区| 五月玫瑰六月丁香| 18禁观看日本| 91国产中文字幕| 亚洲av免费高清在线观看| 日本黄色日本黄色录像| 日韩强制内射视频| 精品人妻一区二区三区麻豆| 亚洲熟女精品中文字幕| 麻豆乱淫一区二区| av网站免费在线观看视频| 大片免费播放器 马上看| 考比视频在线观看| 啦啦啦在线观看免费高清www| 婷婷成人精品国产| 久久99精品国语久久久| 在线 av 中文字幕| 欧美日韩综合久久久久久| 9色porny在线观看| 国产成人午夜福利电影在线观看| 高清av免费在线| 久久精品国产a三级三级三级| av在线观看视频网站免费| 2021少妇久久久久久久久久久| 久久鲁丝午夜福利片| 夫妻午夜视频| 青春草视频在线免费观看| 日韩精品有码人妻一区| 美女视频免费永久观看网站| freevideosex欧美| 免费日韩欧美在线观看| 高清午夜精品一区二区三区| 久久 成人 亚洲| 考比视频在线观看| 国产男女内射视频| 国产免费又黄又爽又色| 一本色道久久久久久精品综合| 99热6这里只有精品| 欧美三级亚洲精品| 一区二区三区乱码不卡18| 日本黄色日本黄色录像| 免费大片18禁| 热99久久久久精品小说推荐| 国产精品欧美亚洲77777| 日本91视频免费播放| 日韩伦理黄色片| 亚洲五月色婷婷综合| 一级毛片黄色毛片免费观看视频| 91精品国产国语对白视频| 尾随美女入室| 亚洲欧洲精品一区二区精品久久久 | 波野结衣二区三区在线| 亚洲精品久久午夜乱码| 人妻一区二区av| 夫妻午夜视频| 久久这里有精品视频免费| 少妇人妻久久综合中文| 久久精品国产亚洲av天美| 午夜av观看不卡| 免费黄色在线免费观看| 欧美精品一区二区大全| 亚洲av免费高清在线观看| 夫妻性生交免费视频一级片| 卡戴珊不雅视频在线播放| 国产毛片在线视频| 18禁在线播放成人免费| 尾随美女入室| 成人18禁高潮啪啪吃奶动态图 | 99热网站在线观看| 精品亚洲乱码少妇综合久久| 搡女人真爽免费视频火全软件| 国产精品久久久久久精品古装| a级毛片在线看网站| 日本欧美国产在线视频| 欧美xxⅹ黑人| 日本欧美国产在线视频| 制服丝袜香蕉在线| 99久国产av精品国产电影| 免费高清在线观看日韩| videossex国产| www.av在线官网国产| 18禁在线播放成人免费| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 色网站视频免费| 搡老乐熟女国产| 亚洲丝袜综合中文字幕| 精品久久久久久久久av| 秋霞伦理黄片| 黄色配什么色好看| 黄色怎么调成土黄色| 97超碰精品成人国产| 男女边吃奶边做爰视频| 夜夜爽夜夜爽视频| 亚洲综合色网址| 国产乱人偷精品视频| 又大又黄又爽视频免费| 欧美成人午夜免费资源| 草草在线视频免费看| 亚洲精品456在线播放app| 久久99精品国语久久久| 欧美丝袜亚洲另类| 只有这里有精品99| 青春草国产在线视频| 久久久亚洲精品成人影院| av不卡在线播放| 你懂的网址亚洲精品在线观看| av国产久精品久网站免费入址| 日韩熟女老妇一区二区性免费视频| 久久人人爽av亚洲精品天堂| 尾随美女入室| 久久久久久久精品精品| 久久久久国产网址| 美女福利国产在线| .国产精品久久| 在线观看免费高清a一片| 亚洲国产日韩一区二区| 少妇猛男粗大的猛烈进出视频| 久久鲁丝午夜福利片| 日本av手机在线免费观看| a级片在线免费高清观看视频| av一本久久久久| 欧美变态另类bdsm刘玥| 久久午夜综合久久蜜桃| 欧美人与善性xxx| 黑人欧美特级aaaaaa片| 伊人久久精品亚洲午夜| 高清不卡的av网站| 国产精品一区二区在线观看99| 美女福利国产在线| 99精国产麻豆久久婷婷| 99久久精品国产国产毛片| 国产熟女欧美一区二区| 免费观看的影片在线观看| 你懂的网址亚洲精品在线观看| 人妻人人澡人人爽人人| 亚洲av日韩在线播放| 国产精品嫩草影院av在线观看| 制服丝袜香蕉在线| 国产 精品1| 国产永久视频网站| 国产乱来视频区| av在线观看视频网站免费| 国产精品国产三级国产专区5o| 18在线观看网站| 一级黄片播放器| 汤姆久久久久久久影院中文字幕| 另类精品久久| 日韩伦理黄色片| 国产爽快片一区二区三区| 视频区图区小说| av线在线观看网站| 亚洲一级一片aⅴ在线观看| 国产色爽女视频免费观看| 黑人高潮一二区| 黄色欧美视频在线观看| 久久久精品94久久精品| 午夜福利,免费看| 色视频在线一区二区三区| 人妻一区二区av| 视频中文字幕在线观看| 精品酒店卫生间| 亚洲成人手机| 天堂中文最新版在线下载| 老司机影院成人| 亚洲av欧美aⅴ国产| 狂野欧美激情性bbbbbb| 国产午夜精品久久久久久一区二区三区| 妹子高潮喷水视频| 久久久久精品久久久久真实原创| videos熟女内射| 日韩三级伦理在线观看| 80岁老熟妇乱子伦牲交| 国产成人av激情在线播放 | 黑人猛操日本美女一级片| 成人午夜精彩视频在线观看| 黄片播放在线免费| 久久久国产精品麻豆| 欧美日韩国产mv在线观看视频| 肉色欧美久久久久久久蜜桃| 国产 精品1| 这个男人来自地球电影免费观看 | 久久久久久久久大av| 亚洲少妇的诱惑av| 国产亚洲欧美精品永久| 亚洲精品一二三| 色视频在线一区二区三区| 黄色一级大片看看| 亚洲精品国产av成人精品| 亚洲国产精品一区二区三区在线| 亚洲精品第二区| 黄片无遮挡物在线观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲四区av| av网站免费在线观看视频| 少妇熟女欧美另类| 欧美人与性动交α欧美精品济南到 | 一本—道久久a久久精品蜜桃钙片| 免费少妇av软件| 欧美另类一区| 国产成人freesex在线| 美女福利国产在线| 国产成人精品久久久久久| 简卡轻食公司| 国产精品无大码| 蜜桃久久精品国产亚洲av| 99热网站在线观看| 午夜免费观看性视频| 国产精品免费大片| 欧美人与善性xxx| 少妇猛男粗大的猛烈进出视频| 99国产综合亚洲精品| 国产精品嫩草影院av在线观看| 亚洲激情五月婷婷啪啪| 亚洲精品456在线播放app| 亚洲怡红院男人天堂| 亚洲av成人精品一区久久| 多毛熟女@视频| 欧美+日韩+精品| 男男h啪啪无遮挡| 欧美三级亚洲精品| 韩国av在线不卡| 观看av在线不卡| 亚洲图色成人| 精品久久久久久久久av| 精品少妇黑人巨大在线播放| 一区二区三区乱码不卡18| 91久久精品电影网| 亚洲精品中文字幕在线视频| 久久久久久伊人网av| 日韩电影二区| 久久久久人妻精品一区果冻| 成人免费观看视频高清| 在线 av 中文字幕| 亚洲美女黄色视频免费看| 蜜桃久久精品国产亚洲av| av又黄又爽大尺度在线免费看| 成人漫画全彩无遮挡| 亚洲在久久综合| videossex国产| 永久免费av网站大全| 久久精品国产亚洲网站| 国产精品秋霞免费鲁丝片| 夜夜爽夜夜爽视频| 国产成人av激情在线播放 | 最近中文字幕高清免费大全6| 夫妻性生交免费视频一级片| 午夜激情av网站| 成人黄色视频免费在线看| 极品少妇高潮喷水抽搐| 精品一区在线观看国产| 性高湖久久久久久久久免费观看| 日韩亚洲欧美综合| 91国产中文字幕| 国产在线一区二区三区精| 免费久久久久久久精品成人欧美视频 | 欧美一级a爱片免费观看看| 中文字幕最新亚洲高清| 色94色欧美一区二区| 中文乱码字字幕精品一区二区三区| 亚洲精品,欧美精品| 亚洲图色成人| 欧美国产精品一级二级三级| 精品酒店卫生间| 国产精品熟女久久久久浪| 欧美亚洲 丝袜 人妻 在线| 自线自在国产av| 一级片'在线观看视频| 国产精品99久久久久久久久| 欧美日韩国产mv在线观看视频| 天天躁夜夜躁狠狠久久av| 特大巨黑吊av在线直播| 亚洲综合色网址| 久久久久久久久久成人| 春色校园在线视频观看| 亚洲精品,欧美精品| 十八禁高潮呻吟视频| 欧美另类一区| 色94色欧美一区二区| 人妻少妇偷人精品九色| 3wmmmm亚洲av在线观看| 亚洲三级黄色毛片| 国产av国产精品国产| 男女免费视频国产| 久热久热在线精品观看| 国产男人的电影天堂91| 日本黄色日本黄色录像| 女性被躁到高潮视频| 亚洲精品日韩av片在线观看| 午夜影院在线不卡| 日本av手机在线免费观看| 国产午夜精品久久久久久一区二区三区| 精品人妻一区二区三区麻豆| 久久午夜综合久久蜜桃| 精品国产露脸久久av麻豆| 国产乱来视频区| 99热网站在线观看| 欧美激情 高清一区二区三区| 王馨瑶露胸无遮挡在线观看| 这个男人来自地球电影免费观看 | 夫妻午夜视频| 成人免费观看视频高清| 亚洲人与动物交配视频| 菩萨蛮人人尽说江南好唐韦庄| 少妇被粗大的猛进出69影院 | 欧美日韩av久久| 26uuu在线亚洲综合色| 一级毛片aaaaaa免费看小| 日日摸夜夜添夜夜爱| 成人毛片a级毛片在线播放| 亚洲国产成人一精品久久久| 午夜久久久在线观看| 女人久久www免费人成看片| 精品人妻熟女av久视频| 亚洲av成人精品一二三区| 欧美变态另类bdsm刘玥| 一区二区三区免费毛片| 欧美日韩成人在线一区二区| 少妇精品久久久久久久| 最近中文字幕高清免费大全6| 人妻夜夜爽99麻豆av| 亚洲av不卡在线观看| 午夜久久久在线观看| 久久久国产精品麻豆| 午夜激情福利司机影院| 亚洲av.av天堂| 麻豆成人av视频| 一本一本综合久久| 精品国产乱码久久久久久小说| 日韩一区二区视频免费看| 亚洲成人手机| 欧美日韩成人在线一区二区| 国产高清三级在线| xxxhd国产人妻xxx| 伊人亚洲综合成人网| 亚洲图色成人| 成人漫画全彩无遮挡| 天堂中文最新版在线下载| 飞空精品影院首页| 国产亚洲午夜精品一区二区久久| 亚洲人成77777在线视频| 建设人人有责人人尽责人人享有的| 我的女老师完整版在线观看| 久久精品国产亚洲av涩爱| av女优亚洲男人天堂| 少妇人妻 视频| 国产精品国产av在线观看| 日韩免费高清中文字幕av| 日本色播在线视频| 一区二区日韩欧美中文字幕 | 久久久久精品性色| 国产伦精品一区二区三区视频9| 日韩在线高清观看一区二区三区| 免费观看的影片在线观看| 3wmmmm亚洲av在线观看| 少妇被粗大的猛进出69影院 | 国产日韩欧美视频二区| 色吧在线观看| 成年人免费黄色播放视频| 国产精品一区www在线观看| 少妇人妻久久综合中文| 91精品三级在线观看| a 毛片基地| 久久久久久久久大av| 久久精品国产亚洲av涩爱| 欧美日韩av久久| 亚洲综合色惰| 国产成人freesex在线| 狂野欧美激情性bbbbbb| 亚洲一级一片aⅴ在线观看| 黄色配什么色好看| 少妇熟女欧美另类| 在线观看三级黄色| 国产成人freesex在线| 我要看黄色一级片免费的| 亚洲精品乱码久久久久久按摩| 性色av一级| 少妇人妻久久综合中文| 亚洲精品久久成人aⅴ小说 | 国产高清三级在线| 国精品久久久久久国模美| 最新的欧美精品一区二区| 国产日韩欧美在线精品| 嫩草影院入口| 亚洲色图 男人天堂 中文字幕 | 国产精品一区二区三区四区免费观看|