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

    粵東及閩南近岸上升流對s局地風(fēng)場變化的響應(yīng)

    2016-09-28 01:33:44蔡尚湛靖春生許金電朱大勇
    海洋學(xué)報(bào) 2016年9期
    關(guān)鍵詞:上升流海流粵東

    蔡尚湛,靖春生,許金電,朱大勇

    (1.國家海洋局第三海洋研究所,福建 廈門 361005)

    ?

    粵東及閩南近岸上升流對s局地風(fēng)場變化的響應(yīng)

    蔡尚湛1,靖春生1,許金電1,朱大勇1

    (1.國家海洋局第三海洋研究所,福建 廈門 361005)

    本文利用2010年6-7月的實(shí)測溫鹽、水位、海流等資料,結(jié)合風(fēng)場數(shù)據(jù),討論了在臺風(fēng)影響較小的情況下,粵東及閩南近岸上升流對局地風(fēng)場變化的響應(yīng)特征,主要結(jié)論如下:(1)譜分析結(jié)果顯示,沿岸風(fēng)、水位、海流、近底層水溫均具有3.5~4.0 d、5.0~5.5 d、8.3~9.0 d的波動(dòng)周期,沿岸風(fēng)的變化引起上升流強(qiáng)度在3~9 d周期上的波動(dòng);(2)上升流對局地風(fēng)場變化的響應(yīng)過程如下:利于上升流產(chǎn)生的局地風(fēng)場發(fā)生變化時(shí),沿岸風(fēng)作用下產(chǎn)生的Ekman輸運(yùn)促使的上升流區(qū)水位的下降幅度發(fā)生改變,隨即向岸方向的壓強(qiáng)梯度力也發(fā)生變化,進(jìn)而導(dǎo)致沿岸流及近底層向岸流的增強(qiáng)或減弱,而近底層向岸流強(qiáng)度的改變又會引起近底層水溫的變化;(3)相關(guān)分析及交叉譜分析的結(jié)果表明,沿岸風(fēng)的變化將在3 d以內(nèi)影響上升流區(qū)近底層水溫。以34 m向岸流代表近底層向岸流,則“沿岸風(fēng)—水位—近底層向岸流—近底層水溫”這一過程的響應(yīng)時(shí)間依次為24 h、7 h、27 h左右。

    粵東;閩南;上升流;沿岸風(fēng);水位;向岸流

    1 引言

    沿岸上升流作為一種重要的中尺度物理現(xiàn)象,對漁業(yè)資源的分布有著顯著的影響[1]。南海北部陸架區(qū)是我國主要的季節(jié)性上升流區(qū)之一[2]。關(guān)于粵東及閩南近岸上升流,最早出現(xiàn)于管秉賢和陳上及的報(bào)道[3],之后諸多現(xiàn)場觀測資料均驗(yàn)證了該上升流的存在[4—7]。上升流區(qū)水體具有低溫、高鹽、低氧、高磷的基本特征,同時(shí)其中心位置的時(shí)空變異較大。實(shí)際上,該上升流作為南海北部海盆尺度季節(jié)性上升流的一部分,不同作者所報(bào)道的是陸架區(qū)上升流的若干中心[4]。

    隨著遙感資料應(yīng)用的日益廣泛,粵東及閩南近岸上升流的時(shí)空分布特征得以進(jìn)一步展現(xiàn)。利用SST、海表葉綠素等遙感資料以及走航CTD數(shù)據(jù),Tang等[8]、Shang等[9]、Zhang等[10]分別對1998年、2004年夏季該上升流分布范圍及強(qiáng)度的變化進(jìn)行了分析,Tang等[11]、Hong等[12]則對上升流的年際變化及對海洋生態(tài)系統(tǒng)的影響進(jìn)行了分析。

    西南季風(fēng)被普遍認(rèn)為是影響南海北部上升流強(qiáng)度的重要因素[4,13—14]。莊偉等[5]通過對衛(wèi)星遙感資料的分析指出海面風(fēng)場沿岸分量是上升流強(qiáng)度改變的重要原因。Jing等[15]利用風(fēng)場遙感數(shù)據(jù)和ECOM模型研究了西南風(fēng)對上升流的貢獻(xiàn)。Hong等[12]的研究結(jié)果也表明風(fēng)場的變異是南海北部上升流年際差異的原因。

    現(xiàn)場觀測資料和模型結(jié)果顯示,在西南風(fēng)減弱、消失甚至東北風(fēng)出現(xiàn)的情況下,只要南海北部陸架區(qū)東向的海流依然存在,上升流現(xiàn)象便可得到維持[16—17],南海北部陸架區(qū)大尺度環(huán)流以及獨(dú)特的地形對上升流的生消有著同樣重要的作用。Gan等[18]利用ROMS模型結(jié)果從動(dòng)力機(jī)制角度分析了南海北部寬陸架區(qū)地形的變化如何導(dǎo)致上升流強(qiáng)度的增大。Gan等[16]將模型結(jié)果應(yīng)用于動(dòng)量平衡方程進(jìn)行分析,指出南海北部上升流的形成和維持依賴于沿等深線方向上增強(qiáng)的壓強(qiáng)梯度力以及底摩擦作用導(dǎo)致的跨等深線方向上的海流分量。Wang等[17]利用POM模型結(jié)果分析了局地風(fēng)場、地形對南海北部上升流強(qiáng)度的相對貢獻(xiàn),指出地形的作用在近岸區(qū)域更為顯著,其在表層與風(fēng)場相當(dāng)而在底層強(qiáng)于風(fēng)場。

    地形的變化導(dǎo)致了南海北部上升流系統(tǒng)的空間分布差異,而上升流強(qiáng)度的時(shí)間變化則更多地受到局地風(fēng)場的影響。臺風(fēng)作為夏季影響南海北部的主要極端天氣事件,將對粵東上升流產(chǎn)生劇烈的影響,并引發(fā)上升流區(qū)海水的受迫振蕩和慣性振蕩[19]。本文利用2010年6-7月的實(shí)測溫鹽、水位、海流等資料,結(jié)合風(fēng)場數(shù)據(jù),從海水運(yùn)動(dòng)的角度,著重討論在臺風(fēng)影響較小的情況下,針對局地風(fēng)場的變化,粵東及閩南近岸上升流的響應(yīng)特征。

    2 數(shù)據(jù)來源與說明

    CTD溫、鹽數(shù)據(jù)取自2010年夏季在粵東及閩南近岸海域的觀測結(jié)果,站位設(shè)置見圖1,斷面J和G分別位于惠來外海和汕頭外海,站位設(shè)置從近岸至外海分別為J1-J7、G1-G8。數(shù)據(jù)觀測時(shí)間為2010年6月8-14日,觀測儀器為美國海鳥公司生產(chǎn)的SBE917 plus CTD,現(xiàn)場觀測數(shù)據(jù)經(jīng)分析、處理和質(zhì)量控制后的垂向分辨率為1 m。

    水位、海流、近底層水溫?cái)?shù)據(jù)取自2010年夏季布放在閩南近岸的座底式海床基系統(tǒng)(海床基位置見圖1),該點(diǎn)平均水深為38.7 m。海床基系統(tǒng)上所攜帶儀器包括一臺WHS-300K ADCP和一臺TGR-2050水位計(jì)。水位與近底層水溫?cái)?shù)據(jù)采樣時(shí)間間隔均為10 min;海流觀測范圍為6~34 m,共15層,垂向間隔為2 m,采樣時(shí)間間隔為1 h。本文所用海床基數(shù)據(jù)時(shí)間范圍為2010年6月10日至7月31日。

    風(fēng)場數(shù)據(jù)共有兩種。第一種數(shù)據(jù)為美國國家環(huán)境預(yù)測中心和國家大氣研究中心發(fā)布的NCEP/DOE(NCEP-2)再分析資料。數(shù)據(jù)空間分辨率約為1.875°×1.875°,時(shí)間分辨率為6 h,本文選取粵東近岸海域的P點(diǎn)(位置見圖1)風(fēng)場數(shù)據(jù)進(jìn)行分析。數(shù)據(jù)下載網(wǎng)址為http://apdrc.soest.hawaii.edu/las/v6/dataset?catitem=5870;第二種數(shù)據(jù)為ASCAT衛(wèi)星散射計(jì)(Advanced SCATterometer)所測海面上空10 m處風(fēng)場數(shù)據(jù),數(shù)據(jù)空間分辨率為0.25°×0.25°,時(shí)間分辨率為1 d。數(shù)據(jù)下載網(wǎng)址為http://apdrc.soest.hawaii.edu/las/v6/dataset?catitem=4736。

    SST(海表溫度)數(shù)據(jù)為美國遙感系統(tǒng)(Remote Sensing Systems, RSS)提供的微波和紅外SST融合數(shù)據(jù)(MISST),該數(shù)據(jù)所融合的數(shù)據(jù)源包括TMI、MODIS和AMSR-E的SST數(shù)據(jù)。數(shù)據(jù)空間分辨率為9.76 km,時(shí)間分辨率為1 d。數(shù)據(jù)下載網(wǎng)址:ftp.discover-earth.org/sst/misst/l4/mw_ir。

    臺風(fēng)路徑資料下載于中國臺風(fēng)網(wǎng)(www.typhoon.gov.cn),本文所用資料為2010年6-7月影響南海臺風(fēng)的逐6 h中心位置數(shù)據(jù)。

    本文對上升流強(qiáng)度變化的分析,主要關(guān)注波動(dòng)周期大于慣性周期的情況。以F點(diǎn)為例,該處慣性周期為Tθ=π/(Ω?sinθ)=3.141 6/(7.292×10-5×sin23.51°)=1.080×105(s)=30.00(h)。因此對時(shí)間分辨率小于1 d的數(shù)據(jù)時(shí)間序列(包括水位、海流、近底層水溫?cái)?shù)據(jù)和NCEP風(fēng)場數(shù)據(jù))均先進(jìn)行以1.5 d為截?cái)嘀芷诘牡屯V波。

    P點(diǎn)風(fēng)場數(shù)據(jù)由于時(shí)間分辨率較高,主要用于單點(diǎn)風(fēng)場時(shí)間序列分析;ASCAT風(fēng)場數(shù)據(jù),時(shí)間分辨率較低但空間分辨率較高,主要用于描述研究海域風(fēng)場特征及對P點(diǎn)風(fēng)場進(jìn)行比對驗(yàn)證。

    為了便于分析,結(jié)合岸線及等深線的走向,對P點(diǎn)風(fēng)場矢量及F點(diǎn)海流矢量均進(jìn)行沿岸方向和向岸方向的分解,沿岸方向與正北的夾角為60°,分解方向示意見圖1。

    文中應(yīng)用了相關(guān)分析和交叉譜分析方法,當(dāng)兩列數(shù)據(jù)時(shí)間分辨率不同時(shí),選擇較低時(shí)間分辨率的數(shù)據(jù)序列進(jìn)行計(jì)算。例如,水位數(shù)據(jù)和海流數(shù)據(jù)時(shí)間分辨率分別為10 min、1 h,兩者進(jìn)行相關(guān)分析計(jì)算時(shí)均取時(shí)間分辨率為1 h的數(shù)據(jù)序列。

    3 2010年6-7月南海東北部風(fēng)場特征

    由圖2可見,2010年6-7月,P點(diǎn)風(fēng)場時(shí)間序列與南海東北部(圖3中紅色虛線框)平均風(fēng)場時(shí)間序列吻合良好,兩者的沿岸風(fēng)速、向岸風(fēng)速在趨勢與量級上均較為接近,因此下文將用P點(diǎn)風(fēng)場來表征南海東北部的平均風(fēng)場特征。

    圖1 研究海域海底地形及觀測站點(diǎn)設(shè)置Fig.1 Topography of study area and observational stations

    圖2 2010年6-7月南海東北部風(fēng)場時(shí)間序列Fig.2 Time series of wind field in the northeastern South China Sea during June to July, 2010a.風(fēng)速矢量,b.沿岸風(fēng)分量,c.向岸風(fēng)分量,藍(lán)色為P點(diǎn)NCEP風(fēng)場數(shù)據(jù)時(shí)間序列,紅色為ASCAT風(fēng)場數(shù)據(jù)區(qū)域(圖3中紅色虛線框)平均時(shí)間序列a.Wind vectors, b.alongshore wind speed component, c.cross-shore wind speed component, blue line is the time series of NCEP wind data at P, red line is the time series of ASCAT wind data averaged in NSCS (indicated by a red dotted box in Fig.3)

    圖3 2010年6-7月影響南海的臺風(fēng)路徑圖Fig.3 The track of typhoons influencing the South China Sea during June to July, 2010紅色虛線框表示ASCAT風(fēng)場數(shù)據(jù)進(jìn)行平均計(jì)算的區(qū)域范圍Red dotted box is the area in which the ASCAT wind data are averaged

    由ASCAT逐日風(fēng)場圖可知(圖略),2010年5-6月初,南海東北部尚處于季風(fēng)轉(zhuǎn)換期,西南季風(fēng)僅零星出現(xiàn)于某些時(shí)段,6月9日以后,西南季風(fēng)才開始較為穩(wěn)定地出現(xiàn)。6月9日至7月11日,P點(diǎn)整體上表現(xiàn)為持續(xù)的西南風(fēng)(沿岸風(fēng)速大多為正值,正的沿岸風(fēng)速的平均值達(dá)到3.65 m/s,圖2b),僅在6月27-30日期間,出現(xiàn)微弱的東南風(fēng)(沿岸風(fēng)速的極值僅達(dá)-3.06 m/s,圖2b)。7月12-25日,受到臺風(fēng)的間接影響,P點(diǎn)風(fēng)向轉(zhuǎn)為東南向。7月26日之后,隨著臺風(fēng)影響的減弱,西南風(fēng)重新出現(xiàn)??傮w上看,除了臺風(fēng)影響期間,P點(diǎn)風(fēng)場表現(xiàn)為以西南風(fēng)為主。

    2010年6-7月,影響南海的臺風(fēng)有兩個(gè),分別為201002號臺風(fēng)“康森”和201003號臺風(fēng)“燦都”,兩者的6 h路徑圖及影響時(shí)間見圖3。這兩個(gè)臺風(fēng)都沒有直接經(jīng)過P點(diǎn)和F點(diǎn)所在海區(qū),相隔距離尚遠(yuǎn),P點(diǎn)7月中旬出現(xiàn)的東南風(fēng)只是受到臺風(fēng)外圍風(fēng)場的影響,而F點(diǎn)的水體也沒有經(jīng)歷臺風(fēng)直接過境時(shí)的那種短時(shí)強(qiáng)烈擾動(dòng)和混合,上升流過程沒有受到直接劇烈的破壞。得益于此,下文的分析較好地刻畫出了臺風(fēng)影響較小情況下的上升流短期變化規(guī)律。

    4 粵東及閩南近岸上升流的實(shí)測證據(jù)

    圖4和圖5分別為2010年6月8-14日期間研究海域溫、鹽的平面及斷面分布圖。觀測期間,南海東北部持續(xù)性的西南風(fēng)剛剛出現(xiàn)(圖2),但研究海域已可見明顯的上升流現(xiàn)象,近岸區(qū)域表現(xiàn)出低溫、高鹽特征。表層至10 m,低溫中心位于汕頭至漳浦之間近岸海域,南澳島至東山島之間表層水溫可低于24.1℃??赡苁艿浇逗恿鞯挠绊懀芯亢S虻奈鞅泵娉霈F(xiàn)低鹽水舌,高鹽區(qū)域位于南澳島與東山島之間,表層鹽度可達(dá)34.3以上。30 m層,西北面低鹽區(qū)不在,低溫、高鹽區(qū)平行岸線分布于東山島以南的近岸區(qū)域以及汕頭西南外海??傮w上看,上升流中心位于調(diào)查范圍的近岸區(qū)域,10 m以淺的鹽度分布可能受到河流淡水的影響而顯得復(fù)雜一些。

    斷面J和G均可見外海深層的低溫、高鹽水沿海底地形向近岸爬升(圖5)。斷面J的外海深層水可直接爬升至近岸區(qū)域,J2站附近10 m以淺可見水溫低于25.0℃、鹽度高于34.0的水體。由于臺灣淺灘的存在,斷面G的外海深層水在G6站附近被阻擋而未能繼續(xù)向近岸推進(jìn),但是該斷面近岸區(qū)域同樣出現(xiàn)了水溫低于24.5℃、鹽度高于34.2的水體,結(jié)合研究海域海水流動(dòng)狀況(圖6)及之前的研究成果[6-7,18]可知,汕頭以西外海深層冷水爬升至近岸,被東北向的沿岸流攜帶至汕頭及其以東的近岸區(qū)域,粵東及閩南近岸上升流的水體來源為汕頭以西的外海深層水。從研究海域30 m層溫、鹽分布圖也可以看出(圖4c、4f),低溫、高鹽水自汕頭西南外海向近岸方向推進(jìn),進(jìn)而向東北方向延伸。

    5 上升流對局地風(fēng)場變化的響應(yīng)

    5.1上升流區(qū)海流特征

    2010年6月10日至7月31日F點(diǎn)垂向海流剖面及水位、沿岸流和向岸流時(shí)間序列見圖6和圖7。由圖6可見,F(xiàn)點(diǎn)海流流速沿岸分量明顯大于向岸分量,表明東北向海流是該處海水流動(dòng)的主要特征。沿岸流速均為正值,且各層變化規(guī)律十分接近,其時(shí)間平均值介于22~32 cm/s之間,并隨水深增大而逐漸減小(圖6b)。向岸流量級雖小,卻體現(xiàn)了典型的上升流區(qū)海水運(yùn)動(dòng)特征(上層海水離岸流動(dòng),中、下層海水向岸流動(dòng)進(jìn)行補(bǔ)償)。從時(shí)間平均上看,10 m以淺為離岸流,流速值在-4 cm/s以內(nèi);10 m以深為向岸流,流速值介于0~6 cm/s之間,并隨水深增大而略微增大,22~34 m流速值相差不大(圖6d)。

    5.2局地風(fēng)場變化對上升流區(qū)海水運(yùn)動(dòng)的影響

    如圖5所示,沿岸上升流的發(fā)生是下層海水沿海底地形向近岸爬升的結(jié)果(即存在海流向岸分量),因此本文以近底層向岸流的強(qiáng)弱來表征上升流的強(qiáng)度。

    F點(diǎn)16 m以深向岸流速值始終為正(圖6c),印證了“西南風(fēng)減弱甚至消失的情況下上升流依舊可以存在”[16—17]這一現(xiàn)象。Shu等[20]也指出局地風(fēng)場的短時(shí)變化幾乎不可能改變南海北部陸架區(qū)大尺度環(huán)流的流向。將沿岸風(fēng)數(shù)據(jù)(圖2b藍(lán)線)與F點(diǎn)沿岸流、向岸流數(shù)據(jù)(圖7b、7c紅線)直接做相關(guān)分析,其結(jié)果也顯示它們之間的相關(guān)性并不顯著(圖8a)??梢姡匕讹L(fēng)變化對上升流區(qū)海流的影響可能只體現(xiàn)在某些頻帶分量上(如下文將要提到的3~9 d周期內(nèi)的變化)。

    事實(shí)上,以往研究結(jié)果表明[17,21—23],沿岸風(fēng)變化對近岸上升流區(qū)海水運(yùn)動(dòng)的影響往往遵循如下這種形式。在沿岸風(fēng)的作用下,近岸區(qū)域上Ekman層的海水離岸輸送,水位下降,導(dǎo)致垂直海岸方向出現(xiàn)壓強(qiáng)梯度力(假設(shè)其方向?yàn)槲鞅毕?,基于地轉(zhuǎn)平衡,東南向的科氏力增大,沿岸流(東北向)將得到增強(qiáng)。在底邊界層內(nèi),增強(qiáng)的沿岸流導(dǎo)致增大的摩擦力(西南向),作為Ekman平衡的結(jié)果,平衡此摩擦力的科氏力(東北向)將會增大,因而向岸流(西北向)也將得到增強(qiáng)。按照上述響應(yīng)順序,下文將對沿岸風(fēng)變化影響上升流的過程進(jìn)行分步討論與驗(yàn)證。

    圖4 2010年夏季粵東及閩南近岸海域水溫、鹽度的平面分布Fig.4 Horizontal distributions of temperature and salinity in eastern Guangdong and southern Fujian coastal seas during summer 2010a中藍(lán)框、紅框分別代表粵東及閩南近岸海域、上升流區(qū),用于SST區(qū)域平均計(jì)算The blue box and red box in a are eastern Guangdong and southern Fujian coastal seas and the upwelling zone respectively, in which the SST data are averaged

    圖5 2010年夏季粵東及閩南近岸海域J和G斷面水溫、鹽度的斷面分布Fig.5 Temperature and salinity distributions along transections J and G in eastern Guangdong and southern Fujian coastal seas during summer 2010

    圖6 2010年6月10日至7月31日F點(diǎn)流速的水深-時(shí)間剖面圖(a.沿岸流,c.向岸流)及平均流速的垂向分布(b.沿岸流,d.向岸流)Fig.6 Depth-time plots of velocities(a. alongshore component, c. cross-shore component) and vertical distributions of averaged velocities(b. alongshore component, d. cross-shore component) at F during June 10 to July 31, 2010

    圖7 F點(diǎn)水位(a)、沿岸流(b, 6~34 m平均)、向岸流(c, 26~34 m平均)隨時(shí)間變化Fig.7 Time series of sea level(a), alongshore velocity (b, 6-34 m averaged) and cross-shore velocity (c, 26-34 m averaged) at F時(shí)間序列紅線為1.5 d低通濾波結(jié)果,藍(lán)線為3~9 d帶通濾波結(jié)果Red lines are the components with periods less than 1.5 d are removed, blue lines are the 3-9 d frequency band components

    沿岸風(fēng)(圖2b藍(lán)線)與水位(圖7a紅線)的相關(guān)分析結(jié)果表明,兩者呈明顯的負(fù)相關(guān)關(guān)系,當(dāng)水位滯后沿岸風(fēng)24 h時(shí),相關(guān)系數(shù)可達(dá)-0.766 6(圖8b)。水位對沿岸風(fēng)變化的響應(yīng)在1 d左右即已完成,這與Wang等[21]的計(jì)算結(jié)果一致。

    對沿岸風(fēng)、水位、沿岸流、向岸流進(jìn)行功率譜分析(圖9),結(jié)果表明,它們在3.5~4.0 d、5.0~5.5 d以及8.3~9.0 d內(nèi)存在多個(gè)相近的譜峰,這預(yù)示著海流中3~9 d的信號與局地風(fēng)場以及水位的變化存在一定的聯(lián)系。對水位及海流數(shù)據(jù)進(jìn)行3~9 d帶通濾波,比較發(fā)現(xiàn)它們的變化趨勢較為類似(圖7中藍(lán)線)。為此,進(jìn)一步對各層沿岸流、向岸流數(shù)據(jù)以及水位數(shù)據(jù)均進(jìn)行3~9 d帶通濾波,并分別進(jìn)行相關(guān)分析,結(jié)果見圖8c、8d。

    沿岸流與水位呈明顯的負(fù)相關(guān)關(guān)系,最大相關(guān)系數(shù)均在-0.74以上,相應(yīng)的滯后時(shí)間(沿岸流滯后水位)均在5 h以內(nèi)(圖8c)。下層向岸流(26~34 m)也與水位呈較明顯的負(fù)相關(guān)關(guān)系,最大相關(guān)系數(shù)均在-0.64以上,相應(yīng)的滯后時(shí)間(向岸流滯后水位)均在7 h以內(nèi)(圖8d)。

    可見,局地風(fēng)場的變化,將首先引起水位的變化,進(jìn)而引發(fā)下層向岸流在3~9 d周期內(nèi)的波動(dòng),從而影響上升流的強(qiáng)度。這個(gè)過程可在2 d內(nèi)完成(以34 m向岸流為例,其滯后沿岸風(fēng)的時(shí)間為24 h+7 h=31 h,見圖8b、8d)。

    為了驗(yàn)證圖8b、8d的分析結(jié)果,直接對沿岸風(fēng)數(shù)據(jù)(即圖2b藍(lán)線)與26~34 m平均向岸流數(shù)據(jù)(圖7c紅線)進(jìn)行交叉譜分析,結(jié)果見圖9e、9f。兩者在3.86 d、5.40 d、9.00 d的周期上相關(guān),且向岸流的波動(dòng)晚于沿岸風(fēng)(位相差分別為157.48°、43.58°、36.83°,所對應(yīng)的滯后時(shí)間分別為40.3 h、15.7 h、22.1 h)。滯后時(shí)間從總體上看在2 d以內(nèi),這與圖8b、8d的結(jié)果類似。

    5.3上升流區(qū)水溫變化特征

    由圖10a可知,2010年5-9月,整個(gè)粵東及閩南近岸海域(圖4a中藍(lán)色虛線框),包括沿岸上升流區(qū)(圖4a中紅色虛線框),SST表現(xiàn)出極為一致的變化規(guī)律。5月1日至6月3日,研究海域尚處于季風(fēng)轉(zhuǎn)換期,SST處于相對較低的水平;6月4-23日,SST迅速增大;6月24日至8月19日,SST處于平穩(wěn)升溫期;8月20日至9月30日,SST處于平穩(wěn)降溫期。

    上升流水體未必時(shí)時(shí)在海表“露頭”,有時(shí)又與臺灣海峽西岸表層低溫水連在一起,同時(shí)SST還受到諸如太陽輻射等多種因素的影響,因此從南海東北部逐日SST分布圖(圖略)上看,并不能很好地反映出上升流強(qiáng)度的時(shí)空變化規(guī)律。

    6月10日至6月下旬,整個(gè)粵東及閩南近岸海域的海水均處于迅速升溫期,F(xiàn)點(diǎn)近底層水溫體現(xiàn)出的是和SST一樣明顯的升溫過程,上升流特征被掩蓋。此后的平穩(wěn)升溫期(6月29日至7月31日),近底層水溫的變化較好地反映了近底層向岸流的波動(dòng)特征。

    圖8 2010年6月10日至7月31日P點(diǎn)沿岸風(fēng)、F點(diǎn)水位和海流之間的相關(guān)分析結(jié)果Fig.8 Results of correlation analyses among the alongshore wind speed component at P, sea level and velocities at F during June 10 to July 31, 2010a.沿岸風(fēng)與沿岸流(6~34 m平均,紅色)、向岸流(26~34 m平均,藍(lán)色)的相關(guān)分析,b.沿岸風(fēng)與水位的相關(guān)分析,c.水位與沿岸流在3~9 d頻帶上的相關(guān)分析,d.水位與向岸流在3~9 d頻帶上的相關(guān)分析。圖c、d中藍(lán)色為最大相關(guān)系數(shù)(指絕對值最大),紅色為對應(yīng)的滯后時(shí)間a.Correlation analysis of alongshore wind speed component with alongshore velocity(6-34 m averaged, red) and cross-shore velocity(26-34 m averaged, blue), b.correlation analysis of alongshore wind speed component with sea level, c.correlation analyses of sea level with alongshore velocities in 3-9 d frequency band, d.correlation analyses of sea level with cross-shore velocities in 3-9 d frequency band. The maximum correlation coefficients (maximum in absolute value) and the corresponding time lag are shown in blue and red in c and d

    圖9 2010年6月10日至7月31日P點(diǎn)沿岸風(fēng)(a)、F點(diǎn)水位(b)、沿岸流(c,6~34 m平均)、向岸流(d,26~34 m平均)的功率譜分析以及P點(diǎn)沿岸風(fēng)與F點(diǎn)向岸流(26~34 m平均)的交叉譜分析(e.凝聚譜,f.位相譜)Fig.9 Power spectrums of the alongshore wind speed component at P(a), sea level(b), alongshore velocity(c, 6-34 m averaged) and cross-shore velocity at F(d, 26-34 m averaged), cross-spectral analysis of the alongshore wind speed component at P with cross-shore velocity(26-34 m averaged)at F(e. coherence, f. phase) during June 10 to July 31, 2010

    圖10 2010年夏季研究海域水溫時(shí)間序列(a)、6月29日至7月31日F點(diǎn)向岸流與近底層水溫在3~9 d頻帶上的相關(guān)分析(b)以及P點(diǎn)沿岸風(fēng)與F點(diǎn)近底層水溫(c)的相關(guān)分析Fig.10 Time series of water temperature in study area during summer 2010(a), correlation analysis of cross-shore velocities with near bottom temperature at F in 3-9 d frequency band (b) and correlation analysis of the alongshore wind speed component at P with near bottom temperature at F (c) during June 29 to July 31, 2010b中藍(lán)色為最大相關(guān)系數(shù)(指絕對值最大),紅色為對應(yīng)的滯后時(shí)間The maximum correlation coefficients (maximum in absolute value) and the corresponding time lag are shown in blue and red in b

    對6月29日至7月31日期間F點(diǎn)近底層水溫?cái)?shù)據(jù)進(jìn)行功率譜分析(圖略),結(jié)果表明,其具有3.62 d、5.29 d、8.59 d的譜峰,與下層向岸流具有相近的波動(dòng)周期(圖9d)。對近底層水溫、各層向岸流數(shù)據(jù)進(jìn)行3~9 d帶通濾波,并分別做相關(guān)分析,結(jié)果見圖10b。下層向岸流(26~34 m)與水溫呈明顯的負(fù)相關(guān)關(guān)系(向岸流越強(qiáng),水溫越低),最大相關(guān)系數(shù)均在-0.6以上(30~34 m更是達(dá)到-0.7以上),相應(yīng)的滯后時(shí)間(水溫滯后向岸流)介于27~32 h之間。Wang等[21]分析了粵東上升流區(qū)向岸流速與底邊界層內(nèi)溫度變化趨勢之間的關(guān)系,也得到了相似的結(jié)論。

    5.4上升流對局地風(fēng)場變化的響應(yīng)時(shí)間

    若以34 m向岸流代表近底層向岸流,那么“沿岸風(fēng)變化—導(dǎo)致水位變化—引起近底層向岸流的強(qiáng)度波動(dòng)—影響近底層水溫”整個(gè)過程,所經(jīng)歷的時(shí)間為58 h(24 h+7 h+27 h=58 h,圖8b、8d、10b)。直接對沿岸風(fēng)數(shù)據(jù)(圖2b藍(lán)線)與近底層水溫?cái)?shù)據(jù)(圖10a綠線)進(jìn)行相關(guān)分析,結(jié)果顯示兩者的最大相關(guān)系數(shù)為-0.564 2,相應(yīng)的滯后時(shí)間(水溫滯后沿岸風(fēng))為54 h(圖10c)。上述兩種計(jì)算方法得出的響應(yīng)時(shí)間較接近(58 h和54 h),說明本文所分析描述的“沿岸風(fēng)變化影響上升流強(qiáng)度”這一過程各環(huán)節(jié)的先后順序和響應(yīng)時(shí)間與實(shí)際情況是較為吻合的。從總體上看,沿岸風(fēng)的變化在3 d之內(nèi)即可引起上升流區(qū)近底層水溫的波動(dòng)。

    6 小結(jié)與討論

    本文利用2010年6-7月的實(shí)測溫鹽、水位、海流等資料,結(jié)合風(fēng)場數(shù)據(jù),討論了在臺風(fēng)影響較小的情況下,粵東及閩南近岸上升流對局地風(fēng)場變化的響應(yīng)特征,主要結(jié)論如下:

    (1) 實(shí)測溫鹽結(jié)果顯示,在持續(xù)性的西南季風(fēng)出現(xiàn)之初,上升流即已存在,上升流的水體來源為汕頭以西的外海深層水。

    (2) 譜分析結(jié)果顯示,沿岸風(fēng)、水位、海流、近底層水溫均具有3.5~4.0 d、5.0~5.5 d、8.3~9.0 d的波動(dòng)周期,沿岸風(fēng)的變化引起上升流強(qiáng)度在3~9 d周期上的波動(dòng)。

    (3) 上升流對局地風(fēng)場變化的響應(yīng)過程如下。利于上升流產(chǎn)生的局地風(fēng)場發(fā)生變化時(shí),沿岸風(fēng)作用下產(chǎn)生的Ekman輸運(yùn)促使的上升流區(qū)水位的下降幅度發(fā)生改變,隨即向岸方向的壓強(qiáng)梯度力也發(fā)生變化,進(jìn)而導(dǎo)致沿岸流及近底層向岸流的增強(qiáng)或減弱,而近底層向岸流強(qiáng)度的改變又會引起近底層水溫的變化。

    (4) 相關(guān)分析及交叉譜分析的結(jié)果表明,沿岸風(fēng)的變化將在3 d以內(nèi)影響上升流區(qū)近底層水溫。以34 m向岸流為例,“沿岸風(fēng)變化影響水位”、“水位變化影響近底層向岸流”、“近底層向岸流變化影響近底層水溫”這3個(gè)環(huán)節(jié)的響應(yīng)時(shí)間依次為24 h、7 h、27 h左右。

    本文中所涉及的水位、海流、水溫在3~9 d內(nèi)的變化,屬于亞潮頻波動(dòng),主要為局地風(fēng)場驅(qū)動(dòng)所致。這些波動(dòng)現(xiàn)象,也可能受到遠(yuǎn)地風(fēng)及大氣擾動(dòng)激發(fā)的陸架波的影響。然而,由于中國近岸海區(qū)夏季天氣系統(tǒng)的活動(dòng)不像冬季那樣頻繁而有規(guī)律[24],而作為夏季大氣系統(tǒng)“強(qiáng)擾動(dòng)源”的臺風(fēng),在2010年6-7月并未出現(xiàn)于上升流區(qū)以北的中國近岸海區(qū),因而當(dāng)年夏季(6-7月)研究海域的陸架波傳播現(xiàn)象可能并不明顯。另外,Yankovsky等[25]的研究指出新澤西近岸陸架區(qū)海流的亞慣性波動(dòng),與流經(jīng)該處的哈德遜河流量有關(guān),當(dāng)流量較大時(shí),在淡水與風(fēng)場的相互作用下,海流顯著增強(qiáng)。類似的,汕頭以東沿岸上升流強(qiáng)度的亞潮頻波動(dòng),也可能受到其“上游”珠江、韓江所帶來的淡水的影響,這方面的研究在模型試驗(yàn)中已有提及[26],尚待更多的觀測數(shù)據(jù)驗(yàn)證。

    [1]Fréon P, Barange M, Arístegui J. Eastern boundary upwelling ecosystems: Integrative and comparative approaches[J]. Progress in Oceanography, 2009, 83(1): 1-14.

    [2]洪啟明, 李立. 粵東陸架區(qū)夏季的上升流[J]. 臺灣海峽, 1991, 10(3): 271-277.

    Hong Qiming, Li Li. A study of upwelling over continental shelf off eastern Guangdong[J]. Journal of Oceanography in Taiwan Strait, 1991, 10(3): 271-277.

    [3]管秉賢, 陳上及. 中國近海的海流系統(tǒng)[R]. 青島:中國科學(xué)院海洋研究所, 1964.

    Guan Bingxian, Chen Shangji. Current System in Coastal Waters of China[R]. Qingdao: Institute of Oceanology, Chinese Academy of Sciences, 1964.

    [4]韓舞鷹, 馬克美. 粵東沿岸上升流的研究[J]. 海洋學(xué)報(bào), 1988, 10(1):52-59.

    Han Wuying, Ma Kemei. Study on coastal upwelling in eastern Guangdong[J]. Haiyang Xuebao, 1988, 10(1):52-59.

    [5]莊偉, 王東曉, 吳日升, 等. 2000年夏季福建、廣東沿海上升流的遙感與船舶觀測分析[J]. 大氣科學(xué), 2005, 29(3):438-444.

    Zhuang Wei, Wang Dongxiao, Wu Risheng, et al. Coastal upwelling off eastern Fujian-Guangdong detected by remote sensing[J]. Chinese Journal of Atmospheric Sciences, 2005, 29(3):438-444.

    [6]蔡尚湛, 吳日升, 許金電. 2006年夏季粵東至閩南近岸海域上升流的特征[J]. 臺灣海峽, 2011,30(4): 489-497.

    Cai Shangzhan, Wu Risheng, Xu Jindian. Characteristics of upwelling in eastern Guangdong and southern Fujian coastal waters during 2006 summer[J]. Journal of Oceanography in Taiwan Strait, 2011,30(4): 489-497.

    [7]許金電, 蔡尚湛, 宣莉莉, 等. 粵東至閩南沿岸海域夏季上升流的調(diào)查研究[J]. 熱帶海洋學(xué)報(bào), 2014, 33(2):1-9.

    Xu Jindian, Cai Shangzhan, Xuan Lili, et al. Observational study on summertime upwelling in coastal seas between eastern Guangdong and southern Fujian[J]. Journal of Tropical Oceanography, 2014, 33(2):1-9.

    [8]Tang D L, Kester D R, Ni I H, et al. Upwelling in the Taiwan Strait during the summer monsoon detected by satellite and shipboard measurements[J]. Remote Sensing of Environment, 2002,83(3): 457-471.

    [9]Shang S L, Zhang C Y, Hong H S, et al. Short-term variability of chlorophyll associated with upwelling events in the Taiwan Strait during the southwest monsoon of 1998[J]. Deep-Sea Research Ⅱ, 2004, 51(10): 1113-1127.

    [10]Zhang C Y, Hong H S, Hu C M, et al. Evolution of a coastal upwelling event during summer 2004 in the southern Taiwan Strait[J]. Acta Oceanologica Sinica, 2011, 30(1) :1-6.

    [11]Tang D L, Kawamura H, Guan L. Long-time observation of annual variation of Taiwan Strait upwelling in summer season[J]. Advances in Space Research, 2004, 33(3): 307-312.

    [12]Hong H S, Zhang C Y, Shang S L, et al. Interannual variability of summer coastal upwelling in the Taiwan Strait[J]. Continental Shelf Research, 2009, 29(2): 479-484.

    [13]于文泉. 南海北部上升流的初步探討[J]. 海洋科學(xué),1987(6):7-10.

    Yu Wenquan. A preliminary approach of the upwelling for the northern South China Sea[J]. Marine Sciences, 1987(6):7-10.

    [14]顏廷壯. 中國沿岸上升流成因類型的初步劃分[J]. 海洋通報(bào), 1991, 10(6):1-6.

    Yan Tingzhuang. Features and classification of coastal upwellings off China[J]. Marine Science Bulletin, 1991, 10(6):1-6.

    [15]Jing Z Y, Qi Y Q, Hua Z L, et al. Numerical study on the summer upwelling system in the northern continental shelf of the South China Sea[J]. Continental Shelf Research, 2009, 29(2):467-478.

    [16]Gan J P, Wang J J, Liang L L, et al. A modeling study of the formation, maintenance, and relaxation of upwelling circulation on the Northeastern South China Sea shelf[J]. Deep-Sea Research Ⅱ, 2015, 117:41-52.

    [17]Wang D X, Shu Y Q, Xue H J, et al. Relative contributions of local wind and topography to the coastal upwelling intensity in the northern South China Sea[J]. Journal of Geophysical Research, 2014, 119(4): 2550-2567.

    [18]Gan J P, Cheung A, Guo X G, et al. Intensified upwelling over a widened shelf in the northeastern South China Sea[J]. Journal of Geophysical Research, 2009, 114(C9): C09019.

    [19]潘愛軍, 郭小鋼, 許金電, 等. 粵東沿岸上升流對2006 年夏季臺風(fēng)響應(yīng)[J]. 中國科學(xué), 2011, 41(10):1530-1541.

    Pan Aijun, Guo Xiaogang, Xu Jindian, et al. Responses of Guangdong coastal upwelling to the summertime typhoons of 2006[J]. Science China, 2011, 41(10):1530-1541.

    [20]Shu Y Q, Wang D X, Zhu J, et al. The 4-D structure of upwelling and Pearl River plume in the northern South China Sea during summer 2008 revealed by a data assimilation model[J]. Ocean Modelling, 2011, 36(3): 228-241.

    [21]Wang D X, Zhuang W, Xie S P, et al. Coastal upwelling in summer 2000 in the northeastern South China Sea[J]. Journal of Geophysical Research, 2012, 117(C4): C04009.

    [22]Weisberg R H, Black B D, Li Z J. An upwelling case study on Florida’s west coast[J]. Journal of Geophysical Research, 2000, 105(C5):11459-11469.

    [23]Gan J P, Allen J S. A modeling study of shelf circulation off northern California in the region of the Coastal Ocean Dynamics Experiment: Response to relaxation of upwelling winds[J]. Journal of Geophysical Research, 2002, 107(C9): 3123.

    [24]Chen D K, Su J L. Continental shelf waves along the coasts of China[J]. Acta Oceanologica Sinica, 1987, 6(3): 317-334.

    [25]Yankovsky A E, Garvine R W. Subinertial dynamics on the inner New Jersey Shelf during the upwelling season[J]. Journal of Physical Oceanography, 1998, 28(12): 2444-2458.

    [26]Gan J P, Li L, Wang D X, et al. Interaction of a river plume with coastal upwelling in the northeastern South China Sea[J]. Continental Shelf Research , 2009, 29(4):728-740.

    錢洪寶,黃曉冬,田紀(jì)偉. 一個(gè)典型南海北部第二模態(tài)內(nèi)孤立波的觀測分析[J]. 海洋學(xué)報(bào), 2016, 38(9):13-20, doi:10.3969/j.issn.0253-4193.2016.09.002

    Qian Hongbao, Huang Xiaodong, Tian Jiwei. Observational study of one prototypical mode-2 internal solitary waves in the northern South China Sea[J]. Haiyang Xuebao, 2016, 38(9): 13-20,doi:10.3969/j.issn.0253-4193.2016.09.002

    Response of upwelling in eastern Guangdong and southern Fujian coastal seas to the local wind variation

    Cai Shangzhan1, Jing Chunsheng1, Xu Jindian1, Zhu Dayong1

    (1.ThirdInstituteofOceanography,StateOceanicAdministration,Xiamen361005,China)

    Based on the cruise conductivity-temperature-depth(CTD) data, the sea level, current and near bottom temperature data from seabed-based observations, and the wind field data during June to July, 2010, this paper studied the response of upwelling in eastern Guangdong and southern Fujian coastal seas to the local wind variation, while the influence of typhoon was relatively weak. The main conclusions are as follows: (1) The results of power spectrum analyses showed that, the alongshore wind speed component, sea level, current and near bottom temperature showed similar variations with the significant periods of fluctuations in 3.5-4.0 d, 5.0-5.5 d and 8.3-9.0 d. The variation of alongshore wind speed component brought about the fluctuation of the upwelling intensity in 3-9 d frequency band. (2) The response process of upwelling to the local wind variation was like this. Surface water was transported offshore within the surface Ekman layer because of the alongshore component of wind stress. Coastal sea level dropped, a cross-shore pressure gradient set up. The pressure gradient would change following the local wind variation, resulting in enhanced(or weakened) alongshore current and enhanced(or weakened) near bottom cross-shore current. Near bottom temperature changed according to the variation of near bottom cross-shore current. (3) The results of correlation analyses and cross-spectral analysis showed that, the variation of alongshore wind speed component would influence near bottom temperature in upwelling zone within 3 days. The lag time for sea level’s response to alongshore wind variation, near bottom (34 m) cross-shore current’s response to sea level variation and near bottom temperature’s response to near bottom (34 m) cross-shore current variation were 24 h, 7 h and 27 h respectively.

    eastern Guangdong; southern Fujian; upwelling; alongshore wind speed component; sea level; cross-shore velocity

    蔡尚湛,靖春生,許金電,等. 粵東及閩南近岸上升流對局地風(fēng)場變化的響應(yīng)[J]. 海洋學(xué)報(bào), 2016, 38(9): 1-12,

    10.3969/j.issn.0253-4193.2016.09.001

    Cai Shangzhan, Jing Chunsheng, Xu Jindian,et al. Response of upwelling in eastern Guangdong and southern Fujian coastal seas to the local wind variation[J]. Haiyang Xuebao, 2016, 38(9): 1-12, doi:10.3969/j.issn.0253-4193.2016.09.001

    2015-09-02;

    2016-01-23。

    國家海洋局第三海洋研究所基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(海三科2011012)。

    蔡尚湛(1983—),男,福建省南安市人,助理研究員,主要從事熱帶邊緣海海洋環(huán)境動(dòng)力學(xué)研究。E-mail:caisz@tio.org.cn

    P731.2

    A

    0253-4193(2016)09-0001-12

    猜你喜歡
    上升流海流粵東
    夏季舟山上升流長期遙感觀測與分析
    基于數(shù)據(jù)挖掘和海流要素的船舶導(dǎo)航改進(jìn)研究
    有限水深海流感應(yīng)電磁場數(shù)值模擬?
    遼東半島頂端海域上升流長期變化特征及影響因素*
    海洋與湖沼(2020年1期)2020-02-06 06:01:38
    新型海流能發(fā)電裝置控制系統(tǒng)的研究
    粵東劇界演變發(fā)展及其戲劇史意義
    戲曲研究(2018年3期)2018-03-19 08:47:24
    低速永磁同步海流發(fā)電機(jī)電磁設(shè)計(jì)與仿真分析
    瓊東上升流的年際變化及長期變化趨勢*
    海洋與湖沼(2016年1期)2016-01-15 03:50:49
    中控簽約中海油粵東LNG接收站項(xiàng)目
    用沿岸上升流指數(shù)分析中國東南沿岸風(fēng)生上升流的特征
    亚洲性夜色夜夜综合| 亚洲国产日韩欧美精品在线观看 | 亚洲熟妇中文字幕五十中出| 一区二区三区精品91| 老司机午夜福利在线观看视频| 一区二区三区国产精品乱码| 日本黄色视频三级网站网址| 亚洲激情在线av| 欧美在线一区亚洲| 校园春色视频在线观看| 国产成+人综合+亚洲专区| 好男人电影高清在线观看| 一二三四社区在线视频社区8| 久久久久国内视频| 欧美久久黑人一区二区| 欧美亚洲日本最大视频资源| 波多野结衣高清作品| 日韩高清综合在线| 久久久久精品国产欧美久久久| 国产精品国产高清国产av| 村上凉子中文字幕在线| 亚洲av电影在线进入| 日韩有码中文字幕| 免费一级毛片在线播放高清视频| 18禁观看日本| 99国产综合亚洲精品| 日本一区二区免费在线视频| 成人精品一区二区免费| 免费在线观看影片大全网站| 亚洲av日韩精品久久久久久密| 好男人在线观看高清免费视频 | 亚洲av成人av| bbb黄色大片| 脱女人内裤的视频| 国内少妇人妻偷人精品xxx网站 | 在线永久观看黄色视频| 亚洲成av人片免费观看| 丁香六月欧美| 亚洲精品一区av在线观看| 91成年电影在线观看| 精品熟女少妇八av免费久了| 免费观看精品视频网站| 伊人久久大香线蕉亚洲五| 欧美激情极品国产一区二区三区| 99久久国产精品久久久| 中文字幕精品免费在线观看视频| 一区福利在线观看| 久久久久九九精品影院| 亚洲欧洲精品一区二区精品久久久| 国产精品98久久久久久宅男小说| 一区福利在线观看| 成人手机av| 日本一本二区三区精品| 国产爱豆传媒在线观看 | 中文字幕人妻熟女乱码| 黑人操中国人逼视频| 久久久久久人人人人人| 嫩草影院精品99| 免费在线观看成人毛片| 欧美日韩瑟瑟在线播放| 精品国产美女av久久久久小说| av超薄肉色丝袜交足视频| 一卡2卡三卡四卡精品乱码亚洲| 日本熟妇午夜| 精品国产美女av久久久久小说| 免费看十八禁软件| 啦啦啦观看免费观看视频高清| 亚洲精品国产区一区二| 欧美性猛交黑人性爽| 国产成人av教育| 亚洲av中文字字幕乱码综合 | 婷婷丁香在线五月| 亚洲av五月六月丁香网| 两个人免费观看高清视频| 精品久久久久久久末码| 日日夜夜操网爽| ponron亚洲| 国产三级黄色录像| 女人被狂操c到高潮| 亚洲av片天天在线观看| 久久久久久九九精品二区国产 | 国产三级在线视频| 国产又爽黄色视频| 久久中文字幕人妻熟女| 精品久久久久久久久久免费视频| 成人av一区二区三区在线看| 天天添夜夜摸| 国产成人欧美在线观看| 国产亚洲精品久久久久久毛片| 在线免费观看的www视频| 91国产中文字幕| 怎么达到女性高潮| 日韩免费av在线播放| 免费高清在线观看日韩| 欧美丝袜亚洲另类 | 免费观看精品视频网站| 波多野结衣av一区二区av| 性色av乱码一区二区三区2| 91麻豆av在线| 亚洲熟女毛片儿| 88av欧美| 十分钟在线观看高清视频www| 国产亚洲精品一区二区www| 亚洲国产精品999在线| 69av精品久久久久久| svipshipincom国产片| 久久香蕉精品热| 亚洲精品色激情综合| 久9热在线精品视频| 中文字幕av电影在线播放| 日韩欧美国产一区二区入口| 人成视频在线观看免费观看| 观看免费一级毛片| 日韩欧美国产一区二区入口| av视频在线观看入口| 欧美国产日韩亚洲一区| 久久婷婷成人综合色麻豆| 亚洲无线在线观看| 色av中文字幕| 老汉色∧v一级毛片| 久久久国产精品麻豆| 国产亚洲欧美98| 亚洲自偷自拍图片 自拍| 国内揄拍国产精品人妻在线 | 欧美精品啪啪一区二区三区| 国产精品二区激情视频| 人人妻人人澡人人看| 国产亚洲精品av在线| 美女大奶头视频| 久久精品91无色码中文字幕| 久久久久久久久免费视频了| 日韩免费av在线播放| 亚洲九九香蕉| 最近最新中文字幕大全电影3 | 一级片免费观看大全| 日本一本二区三区精品| 97碰自拍视频| 久久亚洲真实| 一区二区三区国产精品乱码| 国产av一区二区精品久久| 嫁个100分男人电影在线观看| 怎么达到女性高潮| 中文字幕人成人乱码亚洲影| 男人操女人黄网站| 亚洲国产精品合色在线| 亚洲熟妇中文字幕五十中出| 特大巨黑吊av在线直播 | 精品一区二区三区四区五区乱码| 亚洲自偷自拍图片 自拍| 少妇被粗大的猛进出69影院| 免费在线观看亚洲国产| 不卡av一区二区三区| 日韩精品免费视频一区二区三区| 国产一区二区三区在线臀色熟女| 波多野结衣巨乳人妻| 老司机深夜福利视频在线观看| 久久国产乱子伦精品免费另类| 亚洲人成77777在线视频| xxxwww97欧美| 亚洲 欧美一区二区三区| 无限看片的www在线观看| 不卡一级毛片| 午夜福利在线在线| 亚洲精品在线观看二区| 日本 av在线| 免费女性裸体啪啪无遮挡网站| 国产亚洲精品第一综合不卡| 国产精品久久久久久亚洲av鲁大| 欧美激情极品国产一区二区三区| 精品熟女少妇八av免费久了| 夜夜躁狠狠躁天天躁| 欧美三级亚洲精品| 欧美午夜高清在线| 18美女黄网站色大片免费观看| 色综合欧美亚洲国产小说| 亚洲国产精品合色在线| 精品国产一区二区三区四区第35| 欧美成人性av电影在线观看| 久久精品国产99精品国产亚洲性色| 在线看三级毛片| 日韩大尺度精品在线看网址| 精品少妇一区二区三区视频日本电影| 国产主播在线观看一区二区| 夜夜爽天天搞| 性欧美人与动物交配| 国产亚洲精品综合一区在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 又黄又粗又硬又大视频| 日韩欧美国产一区二区入口| 黄片小视频在线播放| 黄片大片在线免费观看| 亚洲人成网站高清观看| ponron亚洲| 精品乱码久久久久久99久播| 国产真实乱freesex| 亚洲av中文字字幕乱码综合 | 精品欧美国产一区二区三| 中文字幕久久专区| 亚洲精品国产区一区二| 精品国产美女av久久久久小说| 亚洲av日韩精品久久久久久密| 日韩视频一区二区在线观看| 中亚洲国语对白在线视频| 在线观看日韩欧美| 一区二区三区国产精品乱码| 国产亚洲精品第一综合不卡| 国产亚洲精品久久久久5区| 久久午夜亚洲精品久久| 美女扒开内裤让男人捅视频| 亚洲va日本ⅴa欧美va伊人久久| 村上凉子中文字幕在线| 午夜亚洲福利在线播放| 国产视频内射| 国产高清videossex| 成人av一区二区三区在线看| 变态另类丝袜制服| 亚洲人成伊人成综合网2020| 久久 成人 亚洲| 12—13女人毛片做爰片一| 男女下面进入的视频免费午夜 | 欧美日韩亚洲国产一区二区在线观看| 国产视频一区二区在线看| 色综合亚洲欧美另类图片| 在线观看午夜福利视频| 亚洲精品av麻豆狂野| 每晚都被弄得嗷嗷叫到高潮| 亚洲久久久国产精品| 中文字幕精品亚洲无线码一区 | 中文字幕人成人乱码亚洲影| 久久久久久久精品吃奶| 久久午夜综合久久蜜桃| 国产成人欧美| 免费看美女性在线毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲熟女毛片儿| 波多野结衣高清作品| 国产区一区二久久| 一级黄色大片毛片| 国产精品九九99| 国产精品亚洲av一区麻豆| 欧美国产精品va在线观看不卡| 宅男免费午夜| bbb黄色大片| 国产熟女xx| 日韩欧美一区视频在线观看| 精品欧美国产一区二区三| 亚洲一区二区三区不卡视频| 男女做爰动态图高潮gif福利片| 一本久久中文字幕| 久久青草综合色| 久久欧美精品欧美久久欧美| 亚洲av第一区精品v没综合| 欧美久久黑人一区二区| 一进一出抽搐gif免费好疼| 久久久精品欧美日韩精品| 亚洲成国产人片在线观看| 国产精品国产高清国产av| 亚洲精品在线观看二区| 国产精品自产拍在线观看55亚洲| 99久久国产精品久久久| 91九色精品人成在线观看| 丁香欧美五月| 欧美色欧美亚洲另类二区| 久久久久免费精品人妻一区二区 | 国产精品自产拍在线观看55亚洲| 久久久国产成人精品二区| 欧美最黄视频在线播放免费| 日本黄色视频三级网站网址| 亚洲,欧美精品.| 国产熟女午夜一区二区三区| 热re99久久国产66热| 欧美av亚洲av综合av国产av| 免费在线观看成人毛片| 国产一区二区三区视频了| 亚洲av熟女| 精品久久久久久久人妻蜜臀av| 丰满的人妻完整版| √禁漫天堂资源中文www| 日本五十路高清| 国产av又大| 国产精品免费视频内射| 一夜夜www| 国产精品电影一区二区三区| 亚洲,欧美精品.| 亚洲国产中文字幕在线视频| 最好的美女福利视频网| 日本撒尿小便嘘嘘汇集6| 成人av一区二区三区在线看| 后天国语完整版免费观看| 在线观看66精品国产| 精品免费久久久久久久清纯| av电影中文网址| 日韩三级视频一区二区三区| 免费在线观看日本一区| 亚洲男人天堂网一区| 午夜久久久在线观看| 波多野结衣av一区二区av| 国产国语露脸激情在线看| 久久久久久免费高清国产稀缺| 亚洲精华国产精华精| 国产免费男女视频| 丝袜美腿诱惑在线| 搡老熟女国产l中国老女人| 日韩三级视频一区二区三区| 啦啦啦韩国在线观看视频| 黄色女人牲交| 中出人妻视频一区二区| 日韩欧美国产一区二区入口| 非洲黑人性xxxx精品又粗又长| 一区二区三区激情视频| 午夜免费激情av| 精品福利观看| 亚洲黑人精品在线| 成人国产综合亚洲| 亚洲黑人精品在线| 国语自产精品视频在线第100页| 欧美日韩乱码在线| 超碰成人久久| а√天堂www在线а√下载| 精品卡一卡二卡四卡免费| 啦啦啦 在线观看视频| 亚洲在线自拍视频| 亚洲精品色激情综合| 国产片内射在线| 巨乳人妻的诱惑在线观看| 在线看三级毛片| 好男人电影高清在线观看| 99在线视频只有这里精品首页| 一边摸一边做爽爽视频免费| 久久久久九九精品影院| 精品国产亚洲在线| 成人手机av| 亚洲三区欧美一区| 黄片大片在线免费观看| 亚洲 欧美 日韩 在线 免费| 午夜精品在线福利| 12—13女人毛片做爰片一| 香蕉国产在线看| 亚洲熟女毛片儿| 久久久久国内视频| 日本黄色视频三级网站网址| tocl精华| 少妇粗大呻吟视频| 日韩有码中文字幕| 国产成人欧美| 18禁黄网站禁片午夜丰满| 麻豆久久精品国产亚洲av| 在线观看免费日韩欧美大片| 国产一区二区三区视频了| 国产亚洲av嫩草精品影院| 欧美成人午夜精品| 制服诱惑二区| 日韩欧美 国产精品| 欧美色视频一区免费| 欧美三级亚洲精品| 麻豆国产av国片精品| 亚洲中文字幕一区二区三区有码在线看 | 性欧美人与动物交配| 午夜两性在线视频| 夜夜爽天天搞| 久久国产精品影院| 在线观看午夜福利视频| 精品一区二区三区视频在线观看免费| 99国产精品99久久久久| 亚洲国产欧洲综合997久久, | 国产精品九九99| 成人一区二区视频在线观看| 欧美性猛交黑人性爽| 成人精品一区二区免费| 国产精品精品国产色婷婷| 国产熟女xx| 丝袜美腿诱惑在线| 欧美日韩瑟瑟在线播放| 最近最新免费中文字幕在线| 国产精品1区2区在线观看.| 男女视频在线观看网站免费 | 国产高清激情床上av| 亚洲avbb在线观看| 一进一出好大好爽视频| 亚洲在线自拍视频| 波多野结衣高清作品| 十分钟在线观看高清视频www| 亚洲自拍偷在线| 久久精品国产99精品国产亚洲性色| 亚洲国产毛片av蜜桃av| 亚洲精品国产一区二区精华液| 国产不卡一卡二| 午夜福利成人在线免费观看| 色综合亚洲欧美另类图片| 国产人伦9x9x在线观看| 亚洲成av片中文字幕在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲av成人av| 精品卡一卡二卡四卡免费| 亚洲性夜色夜夜综合| 亚洲专区字幕在线| 日韩av在线大香蕉| 亚洲人成网站在线播放欧美日韩| 精品欧美国产一区二区三| ponron亚洲| 午夜激情福利司机影院| 亚洲欧美日韩无卡精品| 美女扒开内裤让男人捅视频| 他把我摸到了高潮在线观看| 午夜免费鲁丝| 一个人免费在线观看的高清视频| 超碰成人久久| 免费看十八禁软件| 悠悠久久av| 国产精品美女特级片免费视频播放器 | 老司机福利观看| 国内揄拍国产精品人妻在线 | 美女免费视频网站| 亚洲欧美激情综合另类| 久久婷婷成人综合色麻豆| 亚洲aⅴ乱码一区二区在线播放 | 亚洲人成网站高清观看| 国产男靠女视频免费网站| 男人的好看免费观看在线视频 | 亚洲国产精品合色在线| 久久精品成人免费网站| 999久久久精品免费观看国产| 99精品久久久久人妻精品| 男人舔女人的私密视频| 久久久久久亚洲精品国产蜜桃av| 精品国产一区二区三区四区第35| 欧美大码av| 哪里可以看免费的av片| 91在线观看av| 啪啪无遮挡十八禁网站| 精品国产超薄肉色丝袜足j| 久久久国产成人免费| 韩国精品一区二区三区| 国产v大片淫在线免费观看| 亚洲片人在线观看| 亚洲欧美日韩高清在线视频| videosex国产| 免费在线观看视频国产中文字幕亚洲| 精品久久久久久久末码| 12—13女人毛片做爰片一| 国内精品久久久久久久电影| 一a级毛片在线观看| 国产片内射在线| 亚洲精品美女久久久久99蜜臀| 国产1区2区3区精品| 日韩视频一区二区在线观看| 美女国产高潮福利片在线看| 十八禁人妻一区二区| 悠悠久久av| 日本撒尿小便嘘嘘汇集6| 中文字幕人成人乱码亚洲影| 男女那种视频在线观看| 这个男人来自地球电影免费观看| 999久久久精品免费观看国产| 人人妻人人澡人人看| 男女午夜视频在线观看| 精品久久蜜臀av无| 国产私拍福利视频在线观看| 国产精品电影一区二区三区| 国产伦人伦偷精品视频| 亚洲人成网站在线播放欧美日韩| 18禁黄网站禁片午夜丰满| 亚洲国产精品sss在线观看| 成人18禁在线播放| 黄色视频不卡| 熟妇人妻久久中文字幕3abv| 麻豆av在线久日| 99久久久亚洲精品蜜臀av| 黄色片一级片一级黄色片| 一级a爱片免费观看的视频| 国产精品永久免费网站| 国产av不卡久久| 精品国产国语对白av| 69av精品久久久久久| 在线观看www视频免费| 国产极品粉嫩免费观看在线| 18禁黄网站禁片免费观看直播| 亚洲天堂国产精品一区在线| videosex国产| 成年版毛片免费区| 欧美色视频一区免费| 午夜影院日韩av| 成人18禁在线播放| 女性生殖器流出的白浆| 久久久久国产精品人妻aⅴ院| 日韩高清综合在线| 神马国产精品三级电影在线观看 | 精品福利观看| 中文字幕人妻熟女乱码| www日本黄色视频网| 制服人妻中文乱码| 国产精品爽爽va在线观看网站 | 一a级毛片在线观看| 成熟少妇高潮喷水视频| 久久国产乱子伦精品免费另类| 亚洲一区高清亚洲精品| 午夜福利在线在线| 99热6这里只有精品| 午夜视频精品福利| 可以在线观看的亚洲视频| 国产午夜福利久久久久久| 后天国语完整版免费观看| 亚洲美女黄片视频| 久久国产精品人妻蜜桃| 90打野战视频偷拍视频| 欧美日韩精品网址| 国产极品粉嫩免费观看在线| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲午夜理论影院| 亚洲精华国产精华精| 国产爱豆传媒在线观看 | 美女国产高潮福利片在线看| 麻豆国产av国片精品| 国产成人系列免费观看| 男女那种视频在线观看| 日本免费a在线| 午夜两性在线视频| 国产97色在线日韩免费| 一进一出好大好爽视频| 国产亚洲精品av在线| 在线观看午夜福利视频| 女人高潮潮喷娇喘18禁视频| 中文字幕最新亚洲高清| av在线天堂中文字幕| 一级a爱视频在线免费观看| 免费在线观看黄色视频的| 美女免费视频网站| 免费观看精品视频网站| 亚洲一码二码三码区别大吗| 亚洲国产欧美日韩在线播放| 成人亚洲精品av一区二区| 国产av不卡久久| 少妇熟女aⅴ在线视频| 一本一本综合久久| 中文字幕人成人乱码亚洲影| 高潮久久久久久久久久久不卡| 日韩欧美免费精品| 亚洲精品国产区一区二| 亚洲欧美精品综合久久99| АⅤ资源中文在线天堂| 成人国产综合亚洲| 一进一出抽搐动态| 国产免费av片在线观看野外av| 亚洲中文日韩欧美视频| 1024视频免费在线观看| 国产精品影院久久| 91大片在线观看| 久久人妻福利社区极品人妻图片| 久久香蕉精品热| 国产区一区二久久| 亚洲久久久国产精品| www日本在线高清视频| 精品欧美国产一区二区三| 欧美在线黄色| 国产精品电影一区二区三区| 夜夜爽天天搞| 波多野结衣高清无吗| 中文字幕最新亚洲高清| 欧美激情高清一区二区三区| 动漫黄色视频在线观看| www国产在线视频色| 久久精品影院6| 亚洲av五月六月丁香网| 三级毛片av免费| 中文字幕人成人乱码亚洲影| 一级毛片精品| 午夜福利在线观看吧| 日本五十路高清| 国产精品一区二区免费欧美| а√天堂www在线а√下载| 日韩欧美 国产精品| 日韩欧美在线二视频| videosex国产| 中文亚洲av片在线观看爽| 少妇的丰满在线观看| 天堂√8在线中文| 国产黄片美女视频| 99久久精品国产亚洲精品| 18禁国产床啪视频网站| 日韩国内少妇激情av| 别揉我奶头~嗯~啊~动态视频| 国产精品香港三级国产av潘金莲| 久久久久久久久中文| a在线观看视频网站| 精品久久久久久,| 脱女人内裤的视频| 国产精品久久视频播放| 亚洲美女黄片视频| 脱女人内裤的视频| 久久精品91无色码中文字幕| 久久久国产欧美日韩av| 成人一区二区视频在线观看| 国产精品久久视频播放| 精品久久久久久,| 窝窝影院91人妻| 免费看十八禁软件| 美女大奶头视频| 国产精品电影一区二区三区| 麻豆av在线久日| 亚洲国产毛片av蜜桃av| 国产精品野战在线观看| 久久久久九九精品影院| 精品国产一区二区三区四区第35| 在线国产一区二区在线| 禁无遮挡网站| 一区福利在线观看| 好看av亚洲va欧美ⅴa在| 一本一本综合久久| 婷婷六月久久综合丁香| 久久精品亚洲精品国产色婷小说| 中文字幕人妻熟女乱码| 宅男免费午夜| 亚洲午夜精品一区,二区,三区| 99久久精品国产亚洲精品| АⅤ资源中文在线天堂|