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

    2007年長江口鄰近海域夏季上升流演變機(jī)制研究

    2013-03-07 07:12:16曹公平宋金寶
    海洋科學(xué) 2013年1期
    關(guān)鍵詞:上升流長江口風(fēng)場

    曹公平, 宋金寶, 樊 偉

    (1. 中國科學(xué)院 海洋研究所, 山東 青島 266071; 2. 中國科學(xué)院 研究生院, 北京 100049 ; 3. 中國科學(xué)院 海洋環(huán)流與波動重點(diǎn)實(shí)驗(yàn)室, 山東 青島 266071)

    2007年長江口鄰近海域夏季上升流演變機(jī)制研究

    曹公平1,2,3, 宋金寶1,3, 樊 偉1,3

    (1. 中國科學(xué)院 海洋研究所, 山東 青島 266071; 2. 中國科學(xué)院 研究生院, 北京 100049 ; 3. 中國科學(xué)院 海洋環(huán)流與波動重點(diǎn)實(shí)驗(yàn)室, 山東 青島 266071)

    研究長江口鄰近海域夏季上升流強(qiáng)度和空間分布的變化, 對漁業(yè)生產(chǎn)和赤潮的防治具有重要的指導(dǎo)意義。采用2007年6~10月高分辨率衛(wèi)星遙感資料NGSST海表溫度和CCMP風(fēng)場, 通過經(jīng)驗(yàn)正交函數(shù)(EOF)分解和區(qū)域海洋數(shù)值模式(ROMS)研究了該海域夏季上升流的短期演變機(jī)制及其與SST異常的關(guān)系。結(jié)果表明, 夏季上升流強(qiáng)度和范圍存在明顯變化, 是引起該海域 SST異常的重要原因; 風(fēng)場對上升流短期演變起著關(guān)鍵作用, 風(fēng)應(yīng)力旋度對局地上升流變化的影響與沿岸風(fēng)應(yīng)力同等重要; 地形變化影響著上升流中心的分布, 陡而窄的海底凸起容易在順流側(cè)形成較強(qiáng)的上升流中心, 并在逆流側(cè)誘發(fā)下降流。

    長江口; 上升流; 經(jīng)驗(yàn)正交函數(shù)(EOF)分析; 區(qū)域海洋數(shù)值模式(ROMS); 演變機(jī)制

    上升流是近海最重要的海洋過程之一, 通過垂向輸運(yùn)過程對海域溫鹽垂直結(jié)構(gòu)和營養(yǎng)鹽分布產(chǎn)生重要影響, 是海洋漁場形成的重要動力環(huán)境, 也為赤潮爆發(fā)提供條件[1]。以往觀測資料和研究證實(shí), 長江口鄰近海域存在大面積的上升流現(xiàn)象, 包括浙江沿岸上升流區(qū)和長江口外海域上升流區(qū)等。因此, 研究長江口鄰近海域上升流的演變, 不僅有助于漁業(yè)工作者更好地了解漁場的分布與變遷, 而且對赤潮等生態(tài)環(huán)境問題的防治具有重要的指導(dǎo)意義。

    以往不少學(xué)者對該區(qū)域上升流現(xiàn)象進(jìn)行了研究。關(guān)于長江口外夏季上升流現(xiàn)象, 趙保仁等[2-3]認(rèn)為海底地形對臺灣暖流的抬升作用, 是長江口外海形成上升流的基本動力因素; 而朱建榮[4]則認(rèn)為鹽淡水混合產(chǎn)生的斜壓效應(yīng)和臺灣暖流入侵陸架導(dǎo)致的正壓效應(yīng)是形成長江口外水下河谷附近上升流的主要成因。對于浙江近海上升流現(xiàn)象, 胡敦欣[5]指出黑潮北上余脈沿陸架海底沿岸爬坡是其形成的主要動力, 風(fēng)不是主要因子; 而潘玉球[6]則認(rèn)為風(fēng)的作用很重要, 和臺灣暖流產(chǎn)生上升流的量級相當(dāng); 黃祖珂[7]基于數(shù)值模擬提出了潮因子誘發(fā)機(jī)制, 呂新剛[8]進(jìn)一步強(qiáng)調(diào)了潮運(yùn)動的作用。然而, 由于產(chǎn)生上升流的動力因子(臺灣暖流、風(fēng)、潮汐和長江徑流等)復(fù)雜多變, 再加上觀測資料的相對缺乏, 對于該海域上升流的形成機(jī)制至今未形成統(tǒng)一認(rèn)識。而且, 以往研究多集中在風(fēng)應(yīng)力和臺灣暖流等形成上升流的動力機(jī)制上, 對于該海域上升流的演變及其機(jī)制卻鮮有論及。

    海洋衛(wèi)星遙感能夠提供大面積長期的觀測資料,近年來在海洋研究中越來越受重視。Tang[9-10]等應(yīng)用遙感資料研究了臺灣海峽海域的上升流長周期和短周期變化, 胡明娜[11]用遙感數(shù)據(jù)研究了浙江近海夏季上升流現(xiàn)象。本文擬在前人研究的基礎(chǔ)上, 采用高分辨率的海洋遙感資料NGSST和CCMP風(fēng)場, 分析長江口鄰近海域夏季上升流的短期演變, 通過對上升流區(qū)域海表溫度異常 EOF分析, 探討上升流引起短期內(nèi)海表溫度異常的主要機(jī)制, 并借助 ROMS模式研究地形對上升流中心分布的影響。

    1 研究區(qū)域

    本文的研究區(qū)域如圖 1所示。該海域近岸海區(qū)島嶼眾多, 岸線曲折; 海底地形復(fù)雜, 在 20 m等深線與50 m等深線間有一明顯的狹長海底斜坡, 大體與等深線和岸線平行; 河口北端存在水下河谷。北上的臺灣暖流、浙江沿岸流和南下的黃海沿岸流在長江口鄰近海域交匯, 再加上從西邊界注入的長江徑流, 使該海域水系結(jié)構(gòu)異常復(fù)雜。同時, 該海域還受東亞季風(fēng)影響, 氣候復(fù)雜多變。

    圖1 研究區(qū)域海底地形(m)分布Fig. 1 Topography of the investigation area

    2 數(shù)據(jù)和方法

    海表溫度(SST)。Japan GODAE, NGSST(New Generation Sea Surface Temperature)是由日本Tohoku University開發(fā)的一種新的高分辨率海表溫度衛(wèi)星資料, 它綜合了紅外散射計AVHRR、MODIS和微波散射計 AMSR-E三種衛(wèi)星資料, 覆蓋了 116°~166°E, 13°~63°N區(qū)域。其水平分辨率為0.05°×0.05°, 時間分辨率是6 h。

    海面風(fēng)場(SSW)。多平臺交叉校正CCMP(Cross Calibrated Multi-Platform)海面風(fēng)場采用一種增強(qiáng)的變分同化分析法(VAM)融合了QuikSCAT/SeaWinds、ADEOS-II/SeaWinds、AMSR-E、TRMM/TMI和SSM/I等諸多海洋被動微波和散射計遙感平臺上采集的風(fēng)場數(shù)據(jù), 給出的是距海面10 m處風(fēng)場沿經(jīng)向和緯向的速度分量。其空間分辨率為0.25°×0.25°, 時間分辨率是6 h。

    其中U10為距海面 10 m處風(fēng)應(yīng)力矢量,τ為風(fēng)應(yīng)力矢量,x,y分別為向東和向北方向,ρa(bǔ)為空氣密度(1.29 kg/m3),Cd為拖曳系數(shù), 取×10-3[12]。

    本文采用上述遙感衛(wèi)星資料2007年6~10月期間的日平均及月平均資料, 運(yùn)用經(jīng)驗(yàn)正交函數(shù)分析SST異常的時空模態(tài), 并結(jié)合ROMS模式探討上升流的短期演變機(jī)制。

    3 結(jié)果與分析

    3.1 夏季上升流現(xiàn)象

    圖2 NGSST衛(wèi)星資料顯示的2007年月平均SST分布圖Fig. 2 Distributions of monthly average SST during 2007 from satellite data NGSST

    圖2顯示, 2007年夏季6, 7月份長江口外??诟浇S蜷_始出現(xiàn)小面積的冷水團(tuán), 7月21日到8月20日期間, 冷水團(tuán)擴(kuò)大增強(qiáng), 形成顯著的孤立的冷水中心, 8月底冷水中心開始消退。8月份SST分布顯示, 32.5°N以南海域存在兩處孤立的冷水中心, 分別位于啟東附近(32°N, 122.2°E為中心)、長江口外水下河谷和舟山群島附近(29°~31.5°N, 122°~123°E)。冷水中心平均水溫低于 27℃, 比周圍 28.5℃以上水溫低 2℃左右, 由于冷水中心被周圍的暖水包圍, 與以往觀測資料和研究證實(shí)的上升流中心[2,5]基本吻合,因此該冷水中心不應(yīng)是周圍的平流運(yùn)動引起的, 而是深層水上涌形成的, 代表著上升流中心。其中, 較大面積的水舌狀冷水中心由夏季較強(qiáng)的水下河谷附近上升流與舟山群島附近上升流連成一體形成的。水下河谷附近上升流區(qū)比趙保仁[3]等的觀測結(jié)果偏南約半個緯度, 很可能是南下的黃海沿岸流和長江沖淡水在表層驅(qū)動造成的。32.5°N以北海域存在大面積的冷水舌入侵現(xiàn)象, 可能是由于黃海沿岸流、黃海暖流與長江沖淡水匯聚產(chǎn)生氣旋型渦旋造成的[13]。從圖中可以看出, 三處上升流的中心位置基本穩(wěn)定, 上升流的強(qiáng)度和范圍發(fā)生顯著的變化。

    3.2 上升流引起海表溫度異常的EOF結(jié)果分析

    為了研究影響上升流演變的因素, 首先對由上升流引起的海表溫度異常進(jìn)行模態(tài)分解。本文選取上升流海域29°~32.5°N, 121.5°~124.2°E區(qū)域, 上升流現(xiàn)象顯著且變化明顯的2007年7月21日至8月20日期間進(jìn)行研究。首先對日平均NGSST資料進(jìn)行處理, 得到SST距平場。對SST距平場進(jìn)行EOF分解, 得到SST方差變化的空間和時間模態(tài), 其前4個模態(tài)貢獻(xiàn)率見表 1。圖 3和圖 4分別為上升流引起SST異常的EOF分解前4個模態(tài)空間分布和對應(yīng)的時間變化系數(shù)。

    表1 EOF的前4個模態(tài)的方差貢獻(xiàn)率Tab. 1 Contributions of the first four EOF modes to variance

    第一模態(tài)的貢獻(xiàn)占方差的 51.52%, 空間分布上整個區(qū)域均為負(fù)值, 這表明研究海域 SST變化趨勢空間上具有整體一致性。第一模態(tài)的時間變化系數(shù)為正值時, 表現(xiàn)為整個海域海表溫度偏低, 時間系數(shù)為負(fù)值時, 整個海域溫度偏高。上海市中心氣象臺資料顯示上海市和長江口區(qū)2007年7月23日至8月3日出現(xiàn)連續(xù)12 d的持續(xù)高溫, 29日極端高溫達(dá)39.6℃, 時間變化系數(shù)上表現(xiàn)為較大的負(fù)值, 8月 3日之后氣溫開始下降, 時間變化系數(shù)上表現(xiàn)為逐漸向正值轉(zhuǎn)變。第一模態(tài)時間變化系數(shù)很好地反映了夏季上海市和長江口區(qū)的氣溫變化特征, 表明氣溫強(qiáng)迫是SST異常的首要控制因子。SST變率較大值區(qū)域, 沿著等深線從南向北延伸到 32°N附近, 和臺灣暖流夏季入侵的路徑一致[14], 夏季穩(wěn)定的臺灣暖流沿等深線入侵可能是SST出現(xiàn)較大變率的主要原因。因此, 第一模態(tài)反映了研究海域 SST異常主要受氣溫強(qiáng)迫的控制, 具有整體一致性, 同時受臺灣暖流入侵的影響呈現(xiàn)空間變率的不均勻分布。

    第二模態(tài)的貢獻(xiàn)占方差的27.00%??臻g分布上以 31.5°N為界, 以南和以北海域分別呈現(xiàn)負(fù)異常和正異常的空間分布特征。圖5顯示, 31.5°N以北大部分海域是東南風(fēng), 以南海域則是西南風(fēng), 強(qiáng)勁的西南風(fēng)產(chǎn)生 Ekman輸運(yùn)形成風(fēng)生上升流, 使附近表層水溫度降低, 表現(xiàn)為SST負(fù)異常區(qū)域; 8月8日左右風(fēng)場由西南風(fēng)轉(zhuǎn)為東南風(fēng), 在該模態(tài)時間變化系數(shù)上表現(xiàn)為從負(fù)值轉(zhuǎn)化為正值。因此, 第二模態(tài)反映了風(fēng)場形成的上升流引起 SST異?,F(xiàn)象, 表明風(fēng)場是上升流短期演變的關(guān)鍵因子。

    第三模態(tài)的貢獻(xiàn)占方差的9.11%。SST負(fù)異常區(qū)域沿著 123°E斷面 1°左右呈帶狀分布, 恰好沿著長江沖淡水羽狀峰的外海側(cè)[15]。此分布形態(tài)不應(yīng)為長江沖淡水的溫差造成的, 而應(yīng)為長江沖淡水與海水混合產(chǎn)生的較強(qiáng)密度梯度形成的斜壓效應(yīng)在地形的作用下形成上升流, 使得該區(qū)域 SST異常變化較為顯著; 正異常區(qū)域散布在負(fù)異常的外側(cè), 為上升流附帶產(chǎn)生的下降流區(qū)域[4]。水情通報[16]顯示8月 3日和5日, 長江中下游出現(xiàn)大暴雨, 很大程度上影響了入海的實(shí)際長江徑流量, 這正是該模態(tài)時間變化系數(shù)在8月初出現(xiàn)較大正值的原因。因此, 第三模態(tài)反映了長江徑流與海水混合產(chǎn)生的斜壓效應(yīng)誘導(dǎo)的上升流引起 SST異常, 表明長江徑流量的變化是影響水下河谷附近上升流短期變化不可忽視的因子。

    第四模態(tài)的貢獻(xiàn)占方差的3.29%。SST異常區(qū)域以31°N為界呈正、負(fù)相間的分布特征, 和該海域的海底地形基本一致。SST變率最大負(fù)值區(qū)位于啟東附近, 變率最大正值區(qū)位于舟山群島附近。呂新剛等[8]研究表明長江口鄰近海域潮運(yùn)動通過斜壓和正壓機(jī)制誘導(dǎo)上升流。驗(yàn)潮站資料顯示, 長江口南、北側(cè)海域的潮汐運(yùn)動存在較大的相位差; 同時 M2分潮的同潮圖[17]顯示, 啟東和舟山附近的遲角存在60°左右的差異, 潮汐運(yùn)動的區(qū)域差異使SST異常呈現(xiàn)出如模態(tài)四所顯示的空間分布特征。同時, 對應(yīng)的時間變化系數(shù)在 1個月內(nèi)呈現(xiàn)明顯的雙峰周期性,與研究期間天文潮運(yùn)動的月周期出現(xiàn)的大、小潮形態(tài)和時間完全一致。因此, 第四模態(tài)反映了潮運(yùn)動誘導(dǎo)上升流引起的 SST異常, 表明潮運(yùn)動影響上升流短期變化, 但影響相對較小。

    圖3 EOF分解的前4個模態(tài)空間分布Fig.3 The spatial distributions of the four modes by EOF analysis

    EOF分析結(jié)果表明, 長江口鄰近海域太陽輻射和臺灣暖流入侵使 SST異常呈現(xiàn)整體一致的趨勢;風(fēng)場是上升流引起 SST短期變化的關(guān)鍵因素; 長江徑流量產(chǎn)生的斜壓效應(yīng)在水下河谷附近誘導(dǎo)的上升流是該區(qū)域的 SST異常不可忽視的因子; 潮運(yùn)動誘導(dǎo)上升流引起SST異常, 但影響有限。

    圖4 EOF分解的前4個模態(tài)時間變化系數(shù)Fig. 4 The time-varying coefficients of the four modes by EOF analysis

    3.3 風(fēng)應(yīng)力和風(fēng)應(yīng)力旋度對上升流短期演變的影響

    圖5 研究期間每5天平均的風(fēng)應(yīng)力(N/m2)分布Fig. 5 Distributions of every five days average of wind stress (N/m2)

    風(fēng)場通過風(fēng)應(yīng)力產(chǎn)生的Ekman輸運(yùn)和風(fēng)應(yīng)力旋度產(chǎn)生的 Ekman抽吸形成沿岸上升流, 是沿岸上升流形成的重要動力。Enriquez等[18]和Pickett等[19]研究認(rèn)為, 風(fēng)應(yīng)力旋度引起的 Ekman抽吸和風(fēng)應(yīng)力引起的 Ekman輸運(yùn)對于上升流的形成同等重要, 并具有相當(dāng)?shù)牧考?。研究海域夏季盛行西南季風(fēng), 受地形影響形成強(qiáng)的風(fēng)應(yīng)力旋度中心, 對上升流的形成和短期演變有重要作用。

    圖6和圖7分別為研究期間每5天平均的風(fēng)應(yīng)力旋度和SST分布。比較SST分布與同期的風(fēng)應(yīng)力和風(fēng)應(yīng)力旋度分布, 32.5°N以南海域SST低值區(qū)域很好地對應(yīng)著風(fēng)應(yīng)力和風(fēng)應(yīng)力旋度的高值區(qū)域, 同期較大的風(fēng)應(yīng)力旋度引起較冷的SST。7月23日, 風(fēng)應(yīng)力旋度高值中心旋度很小, 平均約為0.01×10-5N/m3,西南向風(fēng)應(yīng)力也較小, 平均不到 0.025 N/m2, 上升流區(qū)范圍較小, 強(qiáng)度較弱, 僅在30°N周圍有經(jīng)緯各1°左右、平均溫度約26℃的上升流中心。8月3日和8日, 風(fēng)應(yīng)力增強(qiáng), 高值中心超過 0.1N/m2, 并在32°N 和舟山群島附近存在兩處風(fēng)應(yīng)力旋度高值中心, 平均約為0.04×10-5N/m3, 分別對應(yīng)著啟東附近上升流區(qū)、舟山群島附近和水下河谷附近上升流區(qū)的SST低值區(qū), 面積較7月23日明顯擴(kuò)大。隨著 8月 13日左右隨著風(fēng)向轉(zhuǎn)向東南風(fēng), 上升流強(qiáng)度開始減弱。

    圖6 研究期間每5天平均風(fēng)應(yīng)力旋度分布Fig. 6 Distributions of the five-day average wind stress curl

    圖 8為長江口水下河谷附近和舟山附近上升流中心區(qū)域SST距平值與該區(qū)域的經(jīng)向風(fēng)應(yīng)力、緯向風(fēng)應(yīng)力和風(fēng)應(yīng)力旋度的距平值的時間序列, SST距平值與局地的經(jīng)向風(fēng)應(yīng)力、緯向風(fēng)應(yīng)力、風(fēng)應(yīng)力旋度相關(guān)性系數(shù)分別為 0.741, 0.334和 0.715, 表明經(jīng)向風(fēng)應(yīng)力和風(fēng)應(yīng)力旋度對局地上升流變化具有重要影響, 而緯向風(fēng)應(yīng)力對局地上升流變化影響不大。因此風(fēng)場對上升流變化的影響是局地的, 風(fēng)應(yīng)力和風(fēng)應(yīng)力旋度對局地上升流短周期的演變具有相同量級的貢獻(xiàn)。

    圖7 研究期間每5天平均SST分布Fig. 7 Distributions of the 5-day average SST

    圖8 上升流區(qū)域SST與局地經(jīng)向(緯向)風(fēng)應(yīng)力及風(fēng)應(yīng)力旋度距平值的時間序列圖Fig. 8 Time series of anomalies for SST, meridional wind stress, zonal wind stress and wind stress curl in upwelling areas

    3.4 潮汐作用對上升流短期演變的影響

    長江口鄰近海域潮運(yùn)動強(qiáng)烈, 以往研究表明潮運(yùn)動通過多種動力機(jī)制誘導(dǎo)上升流, 對該區(qū)域上升流的形成起著重要作用[8]。Pineda[20]研究認(rèn)為春夏季節(jié)潮汐運(yùn)動的月周期性影響著上升流強(qiáng)度的明顯變化, 較大的海表溫度異常發(fā)生在新月后的7~12 d和19~24 d。

    分別截取 ROMS模擬(氣候態(tài)模擬, 范圍27°~35°N, 119°~126°E, 水平分辨率 1/27°, 垂向分24個S層, 海面通量選用COADS05資料, 考慮M2, S2, K1和O1四個主要分潮, 開邊界由中國近海的模擬結(jié)果插值得到)的三處上升流區(qū)域中心 32°N、31°N和29°N斷面8月份月平均的溫度剖面, 和8月份10 m層垂向流速時間序列剖面圖。如圖9所示, 去除潮汐后, 32°N和31°N斷面溫度變化明顯, 29°N斷面, 溫度剖面變化很小。同時, 32°N斷面的122.5°E西側(cè)垂向流速呈現(xiàn)明顯的周期性結(jié)構(gòu), 和天文潮的大、小潮的月周期特征相符。31°N斷面垂向流速123.5°E附近也有類似的周期性, 但其他區(qū)域呈無規(guī)律變化,可能受長江徑流量、風(fēng)場等綜合作用造成的。29°N斷面流速則較弱, 只在8月初較強(qiáng), 可能是由于8月初較強(qiáng)的臺灣暖流與地形的相互作用形成的。

    數(shù)值模擬的結(jié)果和模態(tài)四反映的潮運(yùn)動引起的SST異常分布特征相吻合。因此, 潮運(yùn)動影響上升流強(qiáng)度的短周期演變, 但相對于長江徑流量和風(fēng)場的影響, 其影響較為有限。

    圖9 32°N, 31°N和29°N斷面的有潮、無潮溫度剖面和對應(yīng)10 m層垂向流速時間序列Fig. 9 Vertical distributions of temperature with tide, without tide and time series of vertical velocity in 10 meters, along section 32°N, section 31°N and section 29°N

    3.5 地形變化對上升流中心分布的影響

    地形影響著上升流形成和分布特征, 理論研究表明地形變化并不改變風(fēng)生上升流的總強(qiáng)度, 而是對其重新分配從而形成強(qiáng)的上升流中心, 影響其分布的結(jié)構(gòu)特征[21]。長江口鄰近海域地形極其復(fù)雜, 地形的變化可能是影響上升流中心分布和形態(tài)特征的重要因素。

    如圖10所示, 截取穿過上升流中心和水下河谷,且跨岸方向上地形變化較小的123°E斷面10 m層垂向流速剖面和穿過上升流中心, 且沿岸方向上地形變化較小的 30.7°N斷面 10 m層垂向流速剖面。30.7°N斷面顯示, 近岸斜坡坡度較大, 外側(cè)海底地形變化較平緩。由于底 Ekman效應(yīng), 深層水沿岸爬坡形成上升流, 上升流中心集中在近岸斜坡坡度較大處, 而較平緩的海底地形附近不易形成強(qiáng)的上升流中心。123°E斷面表明, 沿岸地形變化劇烈, 由于底層環(huán)流的爬坡使深層水上涌, 上升流中心多集中在急劇變化地形附近, 窄而陡斜坡的順流側(cè)產(chǎn)生強(qiáng)的上升流, 逆流側(cè)出現(xiàn)下降流。而 29.5°E和31.5°E附近, 底形坡度寬而緩, 垂向流速幾乎為零, 這是由于較寬地形的底摩擦作用抵消了爬坡效應(yīng)和底Ekman效應(yīng)[22]。舟山群島和長江口水下河谷附近海底地形起伏較大, 有利于形成較強(qiáng)的上升流, 是舟山漁場和長江口漁場形成的重要原因。

    圖10 30.7°N和123°E斷面垂向流速剖面圖Fig. 10 Profiles of vertical velocity distribution along section 30.7°N and section 123°E

    4 結(jié)論

    本文應(yīng)用高分辨率的海洋遙感資料NGSST海表溫度和 CCMP風(fēng)場, 研究了長江口鄰近海域上升流的短期變化。對2007年夏季上升流存在期間SST異常進(jìn)行分析, 并通過EOF分析了由上升流引起SST異常的時空分布, 最后通過 ROMS模式進(jìn)一步探討了上升流短期變化的機(jī)制。得到的主要結(jié)論如下: (1)海表溫度異常的時空特征分布較好地反映了上升流的演變。2007年夏季, 該海域上升流范圍和強(qiáng)度變化明顯, 7月下旬清晰顯現(xiàn), 8月上旬達(dá)到最強(qiáng), 8月下旬開始衰退。(2)風(fēng)場是上升流短周期變化的關(guān)鍵因子, 風(fēng)應(yīng)力旋度對局地上升流變化的影響與沿岸風(fēng)應(yīng)力同等重要; 長江徑流量變化對水下河谷附近上升流變化影響顯著; 而潮運(yùn)動對該海域上升流強(qiáng)度變化影響較為有限。(3)地形變化影響著上升流中心的分布位置和形態(tài)。較陡的跨岸斜坡有利于形成強(qiáng)的上升流中心; 起伏的沿岸地形變化, 容易在窄而陡斜坡的順流側(cè)形成強(qiáng)的上升流中心, 逆流側(cè)出現(xiàn)下降流, 緩而寬的斜坡不易形成強(qiáng)的上升流中心。

    長江口鄰近海域分布諸多島嶼和岬角, 海岸線曲折復(fù)雜, 對于風(fēng)場和環(huán)流結(jié)構(gòu)的影響較為復(fù)雜。因此, 海岸線對于上升流的影響不可忽視, 有待進(jìn)一步研究。

    致謝:衛(wèi)星遙感資料 NGSST和 CCMP風(fēng)場來自http://apdrc.soest.hawaii.edu/, 本文還引用了上海中心氣象臺氣溫資料和水位管理系統(tǒng)的長江徑流量資料, 在此一并表示感謝。

    [1] 周名江, 顏天, 鄒景忠. 長江口鄰近海域赤潮發(fā)生區(qū)基本特征初探[J]. 應(yīng)用生態(tài)學(xué)報, 2003, 14 (7): 1031-1038.

    [2] 趙保仁. 長江口外的上升流現(xiàn)象[J]. 海洋學(xué)報, 1993, 15 (2): 108-114.

    [3] 趙保仁, 李徽翡, 楊玉玲. 長江口海區(qū)上升流現(xiàn)象的數(shù)值模擬[C]//中國科學(xué)院海洋研究所. 海洋科學(xué)集刊(45).北京: 科學(xué)出版社, 2003: 64-76.

    [4] 朱建榮. 夏季長江口外水下河谷西側(cè)上升流產(chǎn)生的動力機(jī)制[J]. 科學(xué)通報, 2003, 48 (24): 2754-2758.

    [5] 胡敦欣, 呂良洪, 熊慶成. 關(guān)于浙江沿岸上升流的研究[J]. 科學(xué)通報, 1980, 25 (03): 131-133.

    [6] 潘玉球, 徐端蓉, 許建平. 浙江沿岸上升流區(qū)的鋒面結(jié)構(gòu)、變化及其原因[J]. 海洋學(xué)報, 1985, 7 (4): 401-411.

    [7] 黃祖珂, 俞光耀, 羅義勇,等.東海沿岸潮致上升流的數(shù)值模擬[J]. 青島海洋大學(xué)學(xué)報, 1996, 4: 405-412.

    [8] Lü Xingang, Qiao Fangli, Xia Changshui, et al. Upwelling off Yangtze River estuary in summer[J]. Geophys Res, 2006, 111(C11S08), doi:10.1029/2005JC003250.

    [9] Tang Danling, Kester D R, Ni I-Hsun, et al. Upwelling in the Taiwan Strait during the summer monsoon detected by satellite and shipboard measurements[J]. Remote Sensing Environ, 2002, 83:457-471.

    [10] Tang Danling, Kawamura H, Li Guan. Long-time observation of annual variation of Taiwan Strait upwelling in summer season[J]. Advances in Space Research, 2004, 33: 307-312.

    [11] 胡明娜, 趙朝方. 浙江近海夏季上升流的遙感觀測與分析[J]. 遙感學(xué)報, 2008, 12 (2): 297-304.

    [12] Wu Jin. Wind-stress coefficients over sea surface near neutral conditions a revisit[J]. Phys Oceanogr, 1980, 10: 727-740.

    [13] 毛漢禮, 胡敦欣. 東海北部的一個氣旋型渦旋[C]//中國科學(xué)院海洋研究所. 海洋科學(xué)集刊(27).北京:科學(xué)出版社, 1986: 23-99.

    [14] 翁學(xué)傳, 王從敏. 臺灣暖流深層水變化特征的分析[J]. 海洋與湖沼, 1983, 14 (4): 357-366.

    [15] 朱建榮, 丁平興, 胡敦欣. 2000年8月長江口外海區(qū)沖淡水和羽狀鋒的觀測[J].海洋與湖沼, 2003, 34(3): 249-255.

    [16] 上海市防汛信息中心. 2007年汛期(6-9月)上海地區(qū)水情總結(jié)[J]. 水情通報, 2007, 9: 5-7.

    [17] 于宜法, 劉蘭, 郭明克. 海平面上升導(dǎo)致渤、黃、東海潮波變化的數(shù)值研究Ⅱ—海平面上升后渤、黃、東海潮波的數(shù)值模擬[J]. 中國海洋大學(xué)學(xué)報, 2007, 37(1): 7-14.

    [18] Enriquez A G, Friehe C. Effects of wind stress and wind stress curl variability on coastal upwelling[J]. Phys Oceanogr, 1995, 25: 1651-1671.

    [19] Pickett M H, Paduan J D. Ekman transport and pumping in the California Current based on the U.S. Navy’s high-resolution atmospheric model (COAMPS)[J]. Geophys Res, 2003, 108(C10), 3327, doi: 10.1029/ 2003JC001902.

    [20] Pineda J. An internal tidal bore regime at nearshore stations along western U.S.A.: Predictable upwelling within the lunar cycle[J]. Continental Shelf Research, 1995, 15(8): 1023-1041.

    [21] Song Y T, Chao Y. A theoretical study of topographic effects on coastal upwelling and cross-shore exchange[J]. Ocean Modelling, 2004, 6 (2): 151-176.

    [22] Jiang L, Yan X H, Tseng Y H, et al. A numerical study on the role of wind forcing, bottom topography, and nonhydrostacy in coastal upwelling[J]. Estuarine, Coastal and Shelf Science, 2011, 95(1): 99-109.

    (本文編輯: 劉珊珊)

    Mechanism of upwelling evolvement in the Yangtze River Estuary adjacent waters in summer, 2007

    CAO Gong-ping1,2,3,SONG Jin-bao1,3,FAN Wei1,3
    (1. Institute of Oceanology, the Chinese Academy of Sciences, Qingdao 266071, China; 2. Graduate School of the Chinese Academy of Sciences, Beijing 100049, China; 3. Key Laboratory of Ocean Circulation and Waves, the Chinese Academy of Sciences, Qingdao 266071, China)

    Mar.,2,2012

    the Yangtze River Estuary; upwelling; empirical orthogonal function (EOF) analysis; regional ocean model system (ROMS); evolvement mechanism

    We utilized the high resolution satellite product, NGSST sea surface temperature (SST) and cross calibrated multi-platform (CCMP) wind data, the empirical orthogonal function (EOF) method and numerical ocean model (ROMS) to investigate the evolution of upwelling and its influence on SST anomaly in this specific region. The result indicated that the intensity and spatial distribution of upwelling in summer exhibited an evident variation, which was the key factor for the SST anomaly. The effect of the wind field was local and primarily influences its short term evolution, while the contribution of wind stress curl was on a par with the component of wind stress along the coast. The topographic variation determined the center of the upwelling, which could be induced at the strong inflow flank above the steep and narrow heave, and a corresponding downflow occurs behind the heave.

    P731.21

    A

    1000-3096(2013)01-0102-11

    2012-03-02;

    2012-10-22

    國家基金委創(chuàng)新群體項目(41121064); 國家自然科學(xué)基金(41176016)

    曹公平(1984-), 男, 河南周口人, 碩士研究生, 研究方向:物理海洋學(xué), E-mail: caoyun311@163.com; 宋金寶, 通信作者, 電話: 86-532-82898506, E-mail: songjb@qdio.ac.cn

    猜你喜歡
    上升流長江口風(fēng)場
    夏季舟山上升流長期遙感觀測與分析
    基于FLUENT的下?lián)舯┝魅S風(fēng)場建模
    遼東半島頂端海域上升流長期變化特征及影響因素*
    海洋與湖沼(2020年1期)2020-02-06 06:01:38
    “最美風(fēng)場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    側(cè)向風(fēng)場中無人機(jī)的飛行研究
    長江口鳳鱭繁殖群體的年齡結(jié)構(gòu)和生長特性
    長江口及其鄰近水域仔稚魚種類組成及分布特征
    瓊東上升流的年際變化及長期變化趨勢*
    海洋與湖沼(2016年1期)2016-01-15 03:50:49
    風(fēng)場對夏季長江口外東北海域低鹽水團(tuán)的作用
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    99国产精品99久久久久| 日韩大码丰满熟妇| 五月开心婷婷网| 波野结衣二区三区在线| 亚洲男人天堂网一区| 赤兔流量卡办理| 成年动漫av网址| 亚洲国产毛片av蜜桃av| 国产亚洲精品第一综合不卡| 亚洲精品一区蜜桃| 亚洲国产精品国产精品| 成年av动漫网址| 黑人巨大精品欧美一区二区蜜桃| 亚洲伊人色综图| 精品人妻1区二区| 青春草视频在线免费观看| 赤兔流量卡办理| 久久精品人人爽人人爽视色| 丝袜美腿诱惑在线| 一级毛片黄色毛片免费观看视频| 成人亚洲欧美一区二区av| 久久鲁丝午夜福利片| 精品视频人人做人人爽| 纯流量卡能插随身wifi吗| 欧美日韩亚洲高清精品| 中文乱码字字幕精品一区二区三区| 亚洲欧美一区二区三区久久| 首页视频小说图片口味搜索 | 可以免费在线观看a视频的电影网站| 黄色片一级片一级黄色片| 日韩大片免费观看网站| www.精华液| 18禁黄网站禁片午夜丰满| 久9热在线精品视频| 各种免费的搞黄视频| 伊人亚洲综合成人网| 性高湖久久久久久久久免费观看| 亚洲精品中文字幕在线视频| 日韩 亚洲 欧美在线| 亚洲精品av麻豆狂野| 又粗又硬又长又爽又黄的视频| 国产在线免费精品| 韩国高清视频一区二区三区| 黄色片一级片一级黄色片| 少妇裸体淫交视频免费看高清 | 丝袜美腿诱惑在线| 日韩中文字幕欧美一区二区 | 国产亚洲欧美在线一区二区| 国产欧美日韩一区二区三 | 久久久精品区二区三区| 五月天丁香电影| 中文字幕最新亚洲高清| 成人国产av品久久久| 国产亚洲一区二区精品| 丰满迷人的少妇在线观看| 国产精品av久久久久免费| 极品少妇高潮喷水抽搐| 91九色精品人成在线观看| e午夜精品久久久久久久| 国产成人精品久久二区二区91| 少妇精品久久久久久久| 国产一区二区 视频在线| 男女床上黄色一级片免费看| 国产精品一二三区在线看| 亚洲少妇的诱惑av| 99re6热这里在线精品视频| 老鸭窝网址在线观看| e午夜精品久久久久久久| 久久久久久久久免费视频了| 丝袜美足系列| 日日摸夜夜添夜夜爱| 午夜激情av网站| av线在线观看网站| 天天躁夜夜躁狠狠躁躁| 国产亚洲精品第一综合不卡| 晚上一个人看的免费电影| 女人精品久久久久毛片| 欧美黄色片欧美黄色片| av国产久精品久网站免费入址| 久久 成人 亚洲| 国产精品偷伦视频观看了| 欧美在线一区亚洲| 这个男人来自地球电影免费观看| 男人操女人黄网站| 亚洲精品久久成人aⅴ小说| 久久亚洲精品不卡| 国产深夜福利视频在线观看| 国产淫语在线视频| 久久久国产精品麻豆| 最新的欧美精品一区二区| 久久久久视频综合| 91精品国产国语对白视频| 一级黄色大片毛片| 波多野结衣av一区二区av| 国产一区二区三区综合在线观看| 亚洲国产欧美一区二区综合| av片东京热男人的天堂| 美女扒开内裤让男人捅视频| 欧美大码av| 亚洲精品av麻豆狂野| 一级a爱视频在线免费观看| 国产亚洲av片在线观看秒播厂| 亚洲国产精品999| 亚洲七黄色美女视频| 国产极品粉嫩免费观看在线| 男女国产视频网站| 天天躁日日躁夜夜躁夜夜| 欧美 日韩 精品 国产| 欧美变态另类bdsm刘玥| 国产熟女欧美一区二区| 免费观看a级毛片全部| 欧美中文综合在线视频| 午夜91福利影院| 免费看av在线观看网站| 国产精品久久久人人做人人爽| 欧美激情高清一区二区三区| 99热国产这里只有精品6| 免费久久久久久久精品成人欧美视频| 黑人欧美特级aaaaaa片| 亚洲精品久久午夜乱码| 国产极品粉嫩免费观看在线| 爱豆传媒免费全集在线观看| 精品一区在线观看国产| 在线观看免费高清a一片| 高清欧美精品videossex| 晚上一个人看的免费电影| 国产午夜精品一二区理论片| 狠狠婷婷综合久久久久久88av| 亚洲一卡2卡3卡4卡5卡精品中文| www.熟女人妻精品国产| 国产精品久久久久成人av| 亚洲午夜精品一区,二区,三区| 亚洲国产欧美日韩在线播放| 国产xxxxx性猛交| 欧美老熟妇乱子伦牲交| 久久99一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 另类精品久久| 中文字幕人妻丝袜制服| 两个人看的免费小视频| videos熟女内射| 国产精品久久久久久人妻精品电影 | 在线观看人妻少妇| 亚洲五月婷婷丁香| 中文字幕精品免费在线观看视频| 国产av一区二区精品久久| 国产精品一国产av| 超碰成人久久| 中文字幕人妻丝袜一区二区| 日韩一卡2卡3卡4卡2021年| 老司机影院成人| 亚洲国产精品999| 国产麻豆69| 久久精品亚洲熟妇少妇任你| www.999成人在线观看| 一区二区av电影网| 欧美精品亚洲一区二区| 大香蕉久久网| 黑人猛操日本美女一级片| 欧美 日韩 精品 国产| 国产熟女欧美一区二区| 9色porny在线观看| 久久国产精品男人的天堂亚洲| 又大又爽又粗| 一本综合久久免费| 亚洲成人国产一区在线观看 | 久久精品熟女亚洲av麻豆精品| 久久性视频一级片| 久久精品亚洲熟妇少妇任你| 亚洲第一青青草原| 国产一区亚洲一区在线观看| 每晚都被弄得嗷嗷叫到高潮| 交换朋友夫妻互换小说| 脱女人内裤的视频| 80岁老熟妇乱子伦牲交| 日韩伦理黄色片| 在线 av 中文字幕| 欧美日本中文国产一区发布| 免费久久久久久久精品成人欧美视频| 亚洲av成人不卡在线观看播放网 | www.精华液| 每晚都被弄得嗷嗷叫到高潮| 日韩av不卡免费在线播放| 亚洲精品av麻豆狂野| 国产xxxxx性猛交| 午夜av观看不卡| 欧美乱码精品一区二区三区| 国产不卡av网站在线观看| 女警被强在线播放| 国产精品99久久99久久久不卡| 啦啦啦 在线观看视频| 精品欧美一区二区三区在线| 嫩草影视91久久| 人成视频在线观看免费观看| 精品国产国语对白av| 国产av一区二区精品久久| 18禁黄网站禁片午夜丰满| 曰老女人黄片| 天天添夜夜摸| 搡老乐熟女国产| 老熟女久久久| 久久九九热精品免费| 咕卡用的链子| 日韩,欧美,国产一区二区三区| 亚洲精品乱久久久久久| 好男人电影高清在线观看| 性色av乱码一区二区三区2| 精品第一国产精品| 青春草视频在线免费观看| 婷婷色综合大香蕉| 亚洲精品日本国产第一区| 精品少妇黑人巨大在线播放| 精品一区在线观看国产| 妹子高潮喷水视频| 两个人免费观看高清视频| 国产精品.久久久| 最黄视频免费看| 婷婷色综合大香蕉| 免费高清在线观看视频在线观看| 久久99一区二区三区| 午夜福利在线免费观看网站| 亚洲熟女精品中文字幕| 亚洲国产日韩一区二区| 亚洲 国产 在线| a级毛片黄视频| 美女视频免费永久观看网站| 国产一卡二卡三卡精品| 精品欧美一区二区三区在线| 国产在线一区二区三区精| 亚洲精品久久久久久婷婷小说| 91九色精品人成在线观看| 国产伦理片在线播放av一区| 欧美中文综合在线视频| 波多野结衣av一区二区av| 亚洲美女黄色视频免费看| 超碰成人久久| 成在线人永久免费视频| 国产成人免费无遮挡视频| 午夜免费观看性视频| 亚洲欧洲国产日韩| 免费在线观看影片大全网站 | 一区二区三区乱码不卡18| 一本一本久久a久久精品综合妖精| 国产精品亚洲av一区麻豆| 999精品在线视频| 国产女主播在线喷水免费视频网站| 国产成人免费无遮挡视频| 久久毛片免费看一区二区三区| 精品一区二区三区四区五区乱码 | 久久久久久久久久久久大奶| av不卡在线播放| 国产三级黄色录像| www.熟女人妻精品国产| 波多野结衣av一区二区av| 美女主播在线视频| 另类精品久久| 欧美日韩成人在线一区二区| 97精品久久久久久久久久精品| 成人亚洲欧美一区二区av| 精品一区二区三区av网在线观看 | 黄色视频在线播放观看不卡| 久久精品人人爽人人爽视色| 91老司机精品| 欧美激情 高清一区二区三区| 中文精品一卡2卡3卡4更新| 91字幕亚洲| 欧美日韩成人在线一区二区| 亚洲精品第二区| 成年美女黄网站色视频大全免费| 精品久久久精品久久久| 成年人免费黄色播放视频| 国产亚洲午夜精品一区二区久久| 国产日韩欧美在线精品| 首页视频小说图片口味搜索 | 色综合欧美亚洲国产小说| 久久久久久久国产电影| 亚洲色图 男人天堂 中文字幕| 黄色片一级片一级黄色片| 侵犯人妻中文字幕一二三四区| 人体艺术视频欧美日本| 丝袜美足系列| 91麻豆精品激情在线观看国产 | 两人在一起打扑克的视频| 午夜91福利影院| 一区二区日韩欧美中文字幕| 最新的欧美精品一区二区| 日本av手机在线免费观看| 亚洲精品国产av蜜桃| 亚洲九九香蕉| 国产男人的电影天堂91| 国产片特级美女逼逼视频| 欧美精品亚洲一区二区| 搡老岳熟女国产| 一区二区三区精品91| 真人做人爱边吃奶动态| 每晚都被弄得嗷嗷叫到高潮| 久久人妻福利社区极品人妻图片 | 亚洲中文字幕日韩| 天天躁夜夜躁狠狠久久av| 精品久久久精品久久久| 99精国产麻豆久久婷婷| 欧美日韩国产mv在线观看视频| 国产在线视频一区二区| xxxhd国产人妻xxx| 欧美黑人欧美精品刺激| 久久人人爽av亚洲精品天堂| 深夜精品福利| 亚洲情色 制服丝袜| 国产av一区二区精品久久| 啦啦啦视频在线资源免费观看| 丝袜脚勾引网站| 大片免费播放器 马上看| 自线自在国产av| 精品久久久久久久毛片微露脸 | 一本久久精品| 欧美+亚洲+日韩+国产| 久久久国产精品麻豆| 狂野欧美激情性xxxx| 精品少妇一区二区三区视频日本电影| 在线天堂中文资源库| 可以免费在线观看a视频的电影网站| 一边摸一边做爽爽视频免费| 亚洲成色77777| 国产一区二区三区av在线| 精品久久久精品久久久| 赤兔流量卡办理| 一区二区av电影网| 国产精品99久久99久久久不卡| 国产欧美亚洲国产| 又粗又硬又长又爽又黄的视频| 国产无遮挡羞羞视频在线观看| 欧美日韩综合久久久久久| 国产一区二区 视频在线| 丝袜喷水一区| 免费高清在线观看视频在线观看| 亚洲国产欧美网| 国产欧美日韩一区二区三 | 亚洲精品久久久久久婷婷小说| 咕卡用的链子| 久久亚洲国产成人精品v| 十分钟在线观看高清视频www| 高清欧美精品videossex| 国产精品国产av在线观看| 国产福利在线免费观看视频| 欧美少妇被猛烈插入视频| 妹子高潮喷水视频| 肉色欧美久久久久久久蜜桃| 熟女少妇亚洲综合色aaa.| 亚洲精品在线美女| 国产欧美亚洲国产| 各种免费的搞黄视频| 国产无遮挡羞羞视频在线观看| cao死你这个sao货| 久久国产亚洲av麻豆专区| av欧美777| 国产精品秋霞免费鲁丝片| 一二三四在线观看免费中文在| 欧美日韩亚洲综合一区二区三区_| 中文字幕最新亚洲高清| 午夜影院在线不卡| 在线观看免费午夜福利视频| 狂野欧美激情性bbbbbb| 看十八女毛片水多多多| 黑人猛操日本美女一级片| 女人高潮潮喷娇喘18禁视频| 最新的欧美精品一区二区| 美国免费a级毛片| 一本一本久久a久久精品综合妖精| 18禁裸乳无遮挡动漫免费视频| 蜜桃国产av成人99| 国产成人精品久久久久久| 国产99久久九九免费精品| 久久国产精品男人的天堂亚洲| 国产精品一国产av| 亚洲精品日本国产第一区| 五月开心婷婷网| 成人18禁高潮啪啪吃奶动态图| 在线观看人妻少妇| 波多野结衣一区麻豆| 日韩 亚洲 欧美在线| 99国产综合亚洲精品| xxxhd国产人妻xxx| 赤兔流量卡办理| 日韩大片免费观看网站| 纯流量卡能插随身wifi吗| 亚洲人成电影免费在线| 久久精品国产a三级三级三级| 精品久久久久久电影网| 丝袜喷水一区| 日本猛色少妇xxxxx猛交久久| 一区二区日韩欧美中文字幕| 亚洲av电影在线观看一区二区三区| 亚洲av美国av| 真人做人爱边吃奶动态| 一级毛片黄色毛片免费观看视频| 亚洲午夜精品一区,二区,三区| 国产精品 国内视频| 亚洲欧美精品综合一区二区三区| 一区二区三区精品91| 婷婷色av中文字幕| 久久久精品国产亚洲av高清涩受| 丝袜喷水一区| 免费av中文字幕在线| 蜜桃在线观看..| 两个人看的免费小视频| 国产在线免费精品| 美国免费a级毛片| 久久精品国产亚洲av高清一级| 亚洲欧美中文字幕日韩二区| 亚洲人成电影观看| 午夜精品国产一区二区电影| 日本vs欧美在线观看视频| 亚洲欧美色中文字幕在线| 美女扒开内裤让男人捅视频| a级片在线免费高清观看视频| 黄片播放在线免费| 久久久国产精品麻豆| 一本久久精品| 电影成人av| 日韩中文字幕欧美一区二区 | 久久久国产欧美日韩av| 国产精品 欧美亚洲| 久久亚洲国产成人精品v| xxxhd国产人妻xxx| 99精国产麻豆久久婷婷| 丝袜脚勾引网站| 18在线观看网站| 美女视频免费永久观看网站| 色婷婷av一区二区三区视频| 国产精品香港三级国产av潘金莲 | 久久精品人人爽人人爽视色| 超碰成人久久| 丝袜人妻中文字幕| 国产欧美日韩一区二区三区在线| 亚洲免费av在线视频| 国产免费又黄又爽又色| 欧美精品av麻豆av| 黄色一级大片看看| 精品人妻在线不人妻| 精品一区二区三卡| 亚洲av男天堂| 久久久久久久国产电影| 超色免费av| 在现免费观看毛片| 免费高清在线观看日韩| 色视频在线一区二区三区| 欧美激情 高清一区二区三区| 日韩大片免费观看网站| 狠狠精品人妻久久久久久综合| 麻豆av在线久日| av电影中文网址| 亚洲中文av在线| 久久99热这里只频精品6学生| 亚洲人成网站在线观看播放| 久久久久视频综合| 亚洲伊人色综图| 精品国产乱码久久久久久小说| 超碰97精品在线观看| 18在线观看网站| 国产精品99久久99久久久不卡| 大话2 男鬼变身卡| 无限看片的www在线观看| 成人亚洲欧美一区二区av| 真人做人爱边吃奶动态| 国产在线视频一区二区| 国产一区二区 视频在线| 另类亚洲欧美激情| 免费av中文字幕在线| 精品人妻在线不人妻| 2018国产大陆天天弄谢| 国产伦理片在线播放av一区| 99久久99久久久精品蜜桃| 满18在线观看网站| 热re99久久精品国产66热6| 亚洲国产av新网站| 久热这里只有精品99| 一级毛片女人18水好多 | 欧美日韩av久久| 国产精品成人在线| 老熟女久久久| √禁漫天堂资源中文www| kizo精华| 午夜日韩欧美国产| 电影成人av| 国产高清视频在线播放一区 | 欧美性长视频在线观看| 最新的欧美精品一区二区| 18禁观看日本| 国产男女内射视频| 久久鲁丝午夜福利片| 免费在线观看视频国产中文字幕亚洲 | 久久精品成人免费网站| 好男人电影高清在线观看| av网站在线播放免费| 一本一本久久a久久精品综合妖精| 国产免费一区二区三区四区乱码| 老鸭窝网址在线观看| 午夜视频精品福利| 久久久久国产一级毛片高清牌| 亚洲第一青青草原| 人妻人人澡人人爽人人| 国产精品国产三级专区第一集| 国产深夜福利视频在线观看| 欧美国产精品一级二级三级| av有码第一页| 亚洲欧美中文字幕日韩二区| 美女主播在线视频| 日韩电影二区| 啦啦啦中文免费视频观看日本| 51午夜福利影视在线观看| 性色av一级| 成人国产av品久久久| 少妇粗大呻吟视频| 国产精品人妻久久久影院| 深夜精品福利| 超碰成人久久| 成人黄色视频免费在线看| 国产一区有黄有色的免费视频| 久久国产精品人妻蜜桃| 国产免费现黄频在线看| 亚洲精品日本国产第一区| 成人黄色视频免费在线看| 亚洲欧美色中文字幕在线| 天堂8中文在线网| 一级毛片黄色毛片免费观看视频| 国产精品一二三区在线看| 桃花免费在线播放| 亚洲男人天堂网一区| 国产在视频线精品| 亚洲美女黄色视频免费看| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区 视频在线| 精品国产乱码久久久久久男人| 亚洲欧美激情在线| 首页视频小说图片口味搜索 | 国语对白做爰xxxⅹ性视频网站| 视频在线观看一区二区三区| 精品亚洲乱码少妇综合久久| 制服人妻中文乱码| 精品国产乱码久久久久久男人| 国产亚洲一区二区精品| 国产精品亚洲av一区麻豆| 高清视频免费观看一区二区| 欧美 亚洲 国产 日韩一| 成人手机av| 免费少妇av软件| av在线播放精品| 亚洲av欧美aⅴ国产| 久久人人97超碰香蕉20202| 国产亚洲午夜精品一区二区久久| 人人妻人人澡人人爽人人夜夜| 97精品久久久久久久久久精品| 老司机午夜十八禁免费视频| 免费看不卡的av| 国语对白做爰xxxⅹ性视频网站| 亚洲五月婷婷丁香| 免费少妇av软件| 高清不卡的av网站| 日韩av在线免费看完整版不卡| 97在线人人人人妻| 国产成人av激情在线播放| 久久99一区二区三区| 18禁观看日本| 国产人伦9x9x在线观看| 精品国产乱码久久久久久男人| 婷婷色av中文字幕| 免费久久久久久久精品成人欧美视频| 国产免费现黄频在线看| 亚洲欧美激情在线| 爱豆传媒免费全集在线观看| 成年人免费黄色播放视频| 水蜜桃什么品种好| 操美女的视频在线观看| 91精品三级在线观看| 国产片特级美女逼逼视频| 少妇人妻 视频| a级片在线免费高清观看视频| 欧美 亚洲 国产 日韩一| 欧美日韩一级在线毛片| 国产精品国产av在线观看| 久久午夜综合久久蜜桃| 你懂的网址亚洲精品在线观看| 亚洲国产欧美日韩在线播放| 国产精品欧美亚洲77777| 国产av精品麻豆| 免费看av在线观看网站| 国产精品免费大片| 成人亚洲欧美一区二区av| 免费看av在线观看网站| 日本猛色少妇xxxxx猛交久久| 免费高清在线观看日韩| 欧美日韩黄片免| 欧美中文综合在线视频| 99国产综合亚洲精品| 免费在线观看影片大全网站 | 久久国产精品男人的天堂亚洲| 久久 成人 亚洲| 欧美亚洲 丝袜 人妻 在线| 十八禁高潮呻吟视频| 久久久久久久精品精品| 国产精品 国内视频| 日本91视频免费播放| 高清不卡的av网站| 欧美亚洲 丝袜 人妻 在线| 男人添女人高潮全过程视频| 大片免费播放器 马上看| 欧美精品一区二区免费开放| h视频一区二区三区| 国语对白做爰xxxⅹ性视频网站| 别揉我奶头~嗯~啊~动态视频 | 国产精品偷伦视频观看了| 日韩电影二区| 女人爽到高潮嗷嗷叫在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 久久毛片免费看一区二区三区|