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

    山東半島東北部海洋牧場海域冬季海水流動的時空分布特征

    2021-05-26 13:31:28胡羅煜翟方國劉子洲李培良顧艷鎮(zhèn)孫利元姜慶巖
    海洋科學(xué) 2021年4期
    關(guān)鍵詞:余流分潮經(jīng)向

    胡羅煜, 翟方國, 劉子洲, 李培良, 顧艷鎮(zhèn), 孫利元, 李 欣, 陳 棟, 姜慶巖

    山東半島東北部海洋牧場海域冬季海水流動的時空分布特征

    胡羅煜1, 翟方國1, 劉子洲1, 李培良2, 顧艷鎮(zhèn)1, 孫利元3, 李 欣1, 陳 棟2, 姜慶巖2

    (1. 中國海洋大學(xué) 海洋與大氣學(xué)院, 山東 青島 266100; 2. 浙江大學(xué) 海洋學(xué)院, 浙江 舟山 316021; 3. 山東省水生生物資源養(yǎng)護管理中心, 山東 煙臺 264003)

    為豐富山東半島近岸海洋牧場海域水動力環(huán)境研究, 本文利用2019年12月3日至2020年1月1日在山東半島東北部4個海洋牧場獲取的海流資料, 應(yīng)用功率譜分析、調(diào)和分析、余流主軸分析和相關(guān)分析, 探討冬季各海洋牧場的潮流、余流特征及其影響機制。結(jié)果表明: (1) 各海洋牧場潮流由M2分潮潮流主導(dǎo), 受地形邊界限制, 各主要分潮潮流均為往復(fù)流, 且潮流橢圓主軸平行岸線。(2) 不同海洋牧場呈現(xiàn)不同的余流特征和影響機制。煙臺安源海洋牧場余流大致垂直于岸線流向近岸, 平均流速約為0.9~1.7 cm/s; 日平均流以經(jīng)向流為主, 與經(jīng)向風(fēng)呈顯著正相關(guān), 海水受北風(fēng)強迫在近岸堆積。威海瑜泰海洋牧場余流大致垂直于岸線流向外海, 平均流速約為1.4~1.7 cm/s; 日平均流亦以經(jīng)向流為主, 與經(jīng)向風(fēng)呈顯著負相關(guān), 表層海水受北風(fēng)強迫向近岸堆積, 在近岸產(chǎn)生下降流, 海面以下存在北向的離岸流。威海西港海洋牧場余流為東南向, 平均流速約為2.5~3.0 cm/s, 日平均流具有較為顯著的正壓性。榮成楮島海洋牧場余流為東北向, 平均流流速約為5.6~9.9 cm/s, 日平均流表明海水沿桑溝灣南岸流出海灣, 推測桑溝灣海水在灣內(nèi)逆時針流動。研究成果利于進一步研究山東半島鄰近海域動力環(huán)境的多尺度時空變化特征及其影響機制。

    海洋牧場; 潮流; 余流; 山東沿岸

    隨著海產(chǎn)品市場需求的增長, 傳統(tǒng)的掠奪式的沿岸海水養(yǎng)殖業(yè)過度膨脹, 產(chǎn)生了漁業(yè)種群單一化、個體小型化以及沿海海域水質(zhì)污染等副作用, 這種水產(chǎn)模式無法保證海洋漁業(yè)的可持續(xù)發(fā)展和人們對海產(chǎn)品的需求[1]。針對這些問題, 海洋科學(xué)家們提出建設(shè)“海洋牧場”方案。自1980年來, 隨著人工魚礁、增殖放流、魚類馴化等技術(shù)的實施, 這種新型海水養(yǎng)殖產(chǎn)業(yè)模式, 已經(jīng)被證實在集魚, 提高生物多樣性, 改善水質(zhì)、營養(yǎng)鹽結(jié)構(gòu)等方面取得成效[2]。海洋牧場的建設(shè)可以提高海域的水產(chǎn)生物養(yǎng)殖容量, 這與海域營養(yǎng)物質(zhì)含量的調(diào)整、碳匯能力的提高有關(guān), 后者的改變源于海域葉綠素、溶解氧等物質(zhì)分布的變化, 而海域物質(zhì)分布與水動力環(huán)境息息相關(guān)。已有大量研究討論海流和不同海域生態(tài)系統(tǒng)的相互調(diào)控作用: 在肯尼亞的加茲灣, 潮汐不對稱性導(dǎo)致的水體交換促進河口、海灣與外海營養(yǎng)物的交換[3]; 在日本東北部海岸, 諸多海灣與北太平洋進行水交換, 此處養(yǎng)殖業(yè)、近海商業(yè)性捕撈業(yè)十分繁榮[4]; 在福建三沙灣, 箱式養(yǎng)殖設(shè)備減少養(yǎng)殖水域的流量并增強養(yǎng)殖水域之間水道的流量, 研究表明一定程度的改變網(wǎng)箱布局可以改善三沙灣水交換能力[5];在山東半島東部桑溝灣, 樊星等運用10 d海流資料, 結(jié)合雙阻力模型, 討論了夏季桑溝灣因水產(chǎn)養(yǎng)殖產(chǎn)生的海面邊界層現(xiàn)象, 以及邊界層對潮流的影響[6]。山東周邊的海洋牧場建設(shè)方興未艾, 因而研究海洋牧場海域的水動力環(huán)境能夠為了解區(qū)域物質(zhì)分布特征, 進而解決海洋牧場建設(shè)中的設(shè)施布放方案、增殖放流的選址與放流數(shù)量等難題[7]提供新思路, 可進一步提高海洋牧場之于海上糧倉建設(shè)和海洋生態(tài)可持續(xù)發(fā)展的成效。

    研究海域位于山東半島東北部, 為北黃海西部陸架淺海, 水深淺、岸線寬廣且曲折(圖1)。在研究海域, 潮波運動是海水運動的主要形式: 在山東半島東北部, 煙臺、威海附近海域潮流流速較弱, 此海域位于流場的輻合、輻散區(qū)[8]; 漲潮時, 海水自北黃海向渤海、南黃海流動, 落潮時則相反[9]。山東半島東北部不同海域潮流運動特征不同。例如, 在山東半島北部, 威海灣灣口與灣內(nèi)潮流橢圓有較大差異[10]; 在山東半島東部的成山頭近海, 潮流方向則幾乎與岸線平行[11]; 潮流與地形相互作用則可能產(chǎn)生潮汐岬角鋒等復(fù)雜的潮致流[12]。冬季山東半島東北部存在一支較為穩(wěn)定的流, 即魯北沿岸流。冬季北風(fēng)作用下, 沿岸水在渤南近岸區(qū)域堆積, 因而在山東半島北部海域, 西側(cè)水位相對較高, 驅(qū)動沿岸水向東流動, 在成山頭外海轉(zhuǎn)向南流[13](圖1a), 魯北沿岸流攜帶懸浮物質(zhì), 對沿岸地貌有一定影響[14]。山東半島東部海域動力環(huán)境受到許多學(xué)者的關(guān)注, 王志勇等運用2007年冬季大面積調(diào)查資料討論了山東半島東部余流特征, 余流均為南向, 成山頭以東外海的斷面余流流速接近5 cm/s, 桑溝灣東南處的斷面余流流速小于5 cm/s[15]。前人對山東周邊海域動力機制同樣開展了一定研究, 宋軍等對北黃海區(qū)域的潮汐過程進行動量分析和診斷, 結(jié)果表明在渤海海峽、成山頭海域, 動量方程中的壓強梯度力項與局地變化項最顯著[16]。

    圖1 研究區(qū)域與站位分布圖

    注: a: 研究區(qū)域位于山東半島東北部, 其近海沿岸流為魯北沿岸流; b: 展示4海洋牧場的站位分布, 黑色箭頭指示桑溝灣地理位置; c—f: 海洋牧場海域岸線與10 m等深線分布特征, d中黑色箭頭指示雙島港地理位置, 水深數(shù)據(jù)來自ETOPO2(https://www.ngdc.noaa. gov/mgg/global/etopo2.html), 岸線數(shù)據(jù)為GSHHG(http: //www.soest.hawaii.edu/pwessel/gshhg/)

    前人針對北黃海流場特征開展了豐富的研究, 但在近岸海域, 因船只航行、水產(chǎn)養(yǎng)殖、港工建筑等海事活動頻繁, 中長期流矢量觀測資料獲取難度大, 流場特征有待進一步研究[17]?;诖? 本文利用山東半島東北部4處海洋牧場的聲學(xué)多普勒流速剖面儀(acoustic Doppler current profiler, ADCP)的觀測數(shù)據(jù), 對近岸淺海海水流動進行分析, 綜合討論了主要分潮潮流、月平均流、日平均流的特征, 初步討論風(fēng)、地形對海水流動特征的影響。

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

    本文采用自山東省海洋牧場觀測網(wǎng)獲取的2019年12月的海流資料。海洋牧場分布如圖1所示。觀測站位均位于山東半島東北部離岸3 km以內(nèi)海域, 煙臺安源海洋牧場位于渤海海峽東部, 威海瑜泰海洋牧場位于威海市雙島港外海, 威海西港海洋牧場位于威海瑜泰海洋牧場東部、小石島西南, 榮成楮島海洋牧場位于桑溝灣南岸。海流觀測儀器為岸基有纜在線觀測系統(tǒng)搭載的坐底式ADCP, 儀器每小時觀測5個流速剖面, 水層采樣厚度為1 m。傳感器波束自海底向上發(fā)射, 底層流矢量數(shù)據(jù)距離海底1 m。本文根據(jù)觀測系統(tǒng)搭載的溫鹽深剖面儀(conductivity- temperature-depth, CTD)的同時期水深數(shù)據(jù), 選取海表面以下的流矢量數(shù)據(jù)進行功率譜分析、調(diào)和分析和余流主軸分析, 探討海洋牧場海域海水流動特征(表1)。

    表1 海洋牧場數(shù)據(jù)信息

    本文還采用了兩種輔助數(shù)據(jù): (1) 歐洲中期天氣預(yù)報中心(European centre for medium-range weather forecast, ECMWF)發(fā)布的10 m風(fēng)矢量數(shù)據(jù)(ERA5), 用于探究逐日風(fēng)矢量和日平均流的相關(guān)關(guān)系[18]。風(fēng)矢量資料時間范圍為2019年12月3日至2020年1月1日, 時間分辨率為1 h, 空間分辨率為0.25°×0.25°, 經(jīng)雙線性插值得到海洋牧場10 m風(fēng)矢量資料。(2) HYCOM (hybrid coordinate ocean model)再分析流場資料(https:// www.hycom.org/dataserver/gofs-3pt1/reanalysis), 用于討論魯北沿岸流對海洋牧場余流的影響機制, 時間范圍為2019年12月3日至2020年1月1日, 時間分辨率為3 h, 空間分辨率為1/12°×1/12°。

    2 結(jié)果與分析

    2.1 功率譜分析

    對所選水層的流矢量資料進行垂向平均, 并分別對經(jīng)向流、緯向流分量進行功率譜分析, 得到4處海洋牧場的流矢量經(jīng)緯分量的功率譜, 如圖2所示。同一海洋牧場的經(jīng)向流、緯向流具有較為一致的顯著周期, 所有海洋牧場均在半日周期和全日周期對應(yīng)頻段處譜峰最大, 在1/3、1/4、1/5、1/6 d周期對應(yīng)的高頻頻段出現(xiàn)峰值, 主要為倍潮波和復(fù)合潮波。在煙臺安源海洋牧場、威海瑜泰海洋牧場、威海西港海洋牧場, 高頻處的譜峰遠小于半日周期、全日周期譜峰, 在榮成楮島海洋牧場高頻處譜峰峰值較高。

    2.2 潮流特征

    采用MATLAB的T_TIDE程序包對實測海流資料進行潮流調(diào)和分析[19], 得到M2、S2、O1、K1主要分潮的潮流橢圓參數(shù)(表2), 圖3展示了4個海洋牧場表層原始流矢量和潮流調(diào)和分析重構(gòu)表層流矢量時間序列。

    圖2 海流資料功率譜分析

    表2 潮流調(diào)和常數(shù)表

    所有牧場表層的M2分潮潮流橢圓主軸方向基本與岸線平行(圖4)。對比相同海洋牧場、相同水層的不同分潮潮流主軸傾角, 發(fā)現(xiàn)除威海西港海洋牧場表層O1、K1分潮和榮成楮島海洋牧場底層O1分潮外, S2、O1、K1與M2潮流橢圓主軸傾角差均不超過20°, 分潮潮流最大流速方向基本與岸線平行。相同海洋牧場、相同分潮的潮流橢圓主軸垂向分布顯示, 煙臺安源海洋牧場分潮潮流最大流速方向隨深度增加順時針旋轉(zhuǎn), 威海瑜泰海洋牧場分潮潮流最大流速方向隨深度增加逆時針旋轉(zhuǎn); 威海西港海洋牧場、榮成楮島海洋牧場分潮潮流最大流速方向隨深度的變化較為復(fù)雜, 分潮角頻率不同, 最大流速方向隨深度的變化不同。

    圖3 海洋牧場表層流矢量時間序列

    注: 黑色: 原始流矢量; 紅色: 潮流調(diào)和分析重構(gòu)流矢量

    圖4 潮流橢圓分布圖

    潮流調(diào)和分析表明, 4處海洋牧場海域均屬于半日分潮潮流主導(dǎo)的海區(qū), 且M2分潮潮流最顯著, M2分潮潮流流速最大。對比4處海洋牧場海域, M2、S2分潮潮流在榮成楮島海洋牧場流速最大, 其中M2分潮流速約為7~12 cm/s; 威海瑜泰海洋牧場次之, M2分潮流速約為4~8 cm/s; 煙臺安源海洋牧場、威海西港海洋牧場M2分潮流速小于2 cm/s。O1、K1分潮潮流在威海瑜泰海洋牧場流速最大, 其中K1分潮流速約為3~6 cm/s; 榮成楮島海洋牧場次之, K1分潮潮流速約為2~4 cm/s; 煙臺安源海洋牧場、威海西港海洋牧場K1分潮流速不超過1.1 cm/s。各海洋牧場主要分潮潮流流速均隨深度增加減小。煙臺安源海洋牧場、威海西港海洋牧場M2分潮潮流表底層流速相近。威海瑜泰海洋牧場、榮成楮島海洋牧場M2分潮潮流流速中層較表層約減少10%, 底層較表層減少40%~50%, 潮流流速垂向分布表明海底與岸基的摩擦效應(yīng)使底層潮流流速顯著減小。

    潮流橢圓旋轉(zhuǎn)性如圖4所示, 煙臺安源海洋牧場潮流橢圓基本順時針旋轉(zhuǎn), 威海瑜泰海洋牧場潮流橢圓基本逆時針旋轉(zhuǎn), 威海西港海洋牧場、榮成楮島海洋牧場M2、S2半日分潮潮流基本順時針旋轉(zhuǎn), O1、K1全日分潮潮流橢圓基本逆時針旋轉(zhuǎn)。海洋牧場潮流橢圓橢率均小于0.4, 潮流以往復(fù)運動形式為主。在煙臺安源海洋牧場、威海西港海洋牧場表層, M2分潮潮流橢圓橢率分別為–0.32、–0.27, M2分潮潮流旋轉(zhuǎn)性相對較強。潮流橢圓主軸平行于岸線, 主要分潮最大流速隨深度增加減小, 往復(fù)運動的潮流特征表明4處海洋牧場海域潮流運動受地形邊界影響顯著。

    2.3 平均流特征

    原始流矢量剔除調(diào)和重構(gòu)流矢量后, 首先進行48 h低通濾波, 然后進行觀測時段平均, 得到4處海洋牧場海域觀測時段內(nèi)平均流(圖5)。煙臺安源海洋牧場平均流流速約為0.9~1.7 cm/s, 威海瑜泰海洋牧場平均流流速約為1.4~1.7 cm/s, 威海西港海洋牧場平均流流速約為2.5~3.0 cm/s, 榮成楮島海洋牧場平均流流速約為5.6~9.9 cm/s。自煙臺安源海洋牧場至榮成楮島海洋牧場平均流流速增大, 即在山東半島東北部海域, 位于東部的海洋牧場平均流流速較大。煙臺安源海洋牧場平均流流矢量表層為南向, 底層為西南向, 底層平均流與岸線幾乎垂直, 表層較底層稍向左偏, 且底層流速大于表層, 煙臺安源海洋牧場平均流流速與前人在相鄰海域觀測結(jié)果較為一致[20]。威海瑜泰海洋牧場平均流表層為東北向, 中層為北向, 底層為西北向, 平均流呈現(xiàn)海水輻散的離岸流特征, 且流速底層最大, 中層最小。在威海西港海洋牧場, 平均流流矢量表底層幾乎一致, 為東南向, 體現(xiàn)整層水體向東南岸堆積的特征。榮成楮島海洋牧場表中底層平均流幾乎同為東北向, 平均流方向隨深度增加稍向灣口中心洼地偏轉(zhuǎn), 表中層流速大小相當(dāng), 底層流速約為表層60%。

    圖5 月平均流矢量圖

    2.4 日平均流特征

    剔除實測海流中的調(diào)和重構(gòu)分潮流并進行24 h平均, 得到日平均流特征(圖6a)。煙臺安源海洋牧場日平均流東南、西南向流矢量交替出現(xiàn); 表層流矢量主要為西南向, 但東南向流速較大; 底層流矢量主要為西南向, 且底層日平均流流速強于表層。威海瑜泰海洋牧場日平均流表層主要為東北向, 中層?xùn)|北、西北流矢量交替出現(xiàn), 底層主要為西北向。威海西港海洋牧場表底層日平均流流矢量基本一致, 主要為東南向; 榮成楮島海洋牧場日平均流表中底層呈東北向流動, 隨深度增加流矢量逆時針旋轉(zhuǎn)且流速減小。

    原始流矢量剔除調(diào)和重構(gòu)流矢量后, 進行48 h低通濾波得到逐時余流[21], 對逐時余流經(jīng)向、緯向分量進行二維主成分分析得到余流主成分分量, 余流主成分方向是余流偏差量方差最大的方向, 定義為余流主軸方向[22-23], 垂直于主軸的方向定義為次軸方向][。經(jīng)統(tǒng)計余流主成分之于逐時余流的貢獻率超過75%, 且逐時余流在余流主軸方向的變化最劇烈。如圖6所示, 除煙臺安源海洋牧場底層余流主軸與岸線存在較大夾角外, 各海洋牧場余流主軸與岸線方向大致相似, 煙臺安源海洋牧場表層、威海西港海洋牧場余流主軸為東南-西北方向(圖6b, d), 威海瑜泰海洋牧場、榮成楮島海洋牧場余流主軸為東北—西南方向(圖6c, e), 說明在4處海洋牧場海域, 地形是制約余流變化的一個重要因素。

    圖6中, 橢圓長軸為主軸方向上逐時余流分量的標(biāo)準(zhǔn)差, 橢圓短軸則為逐時余流次軸分量的標(biāo)準(zhǔn)差, 余流標(biāo)準(zhǔn)差體現(xiàn)觀測時段內(nèi)余流變化的劇烈程度。煙臺安源海洋牧場、威海西港海洋牧場標(biāo)準(zhǔn)差較小, 余流變化較小; 威海瑜泰海洋牧場、榮成楮島海洋牧場標(biāo)準(zhǔn)差較大, 余流變化較劇烈。余流標(biāo)準(zhǔn)差垂向分布則說明煙臺安源海洋牧場底層余流變化較表層劇烈; 威海瑜泰海洋牧場、榮成楮島海洋牧場底層余流變化較表、中層顯著減小; 威海西港海洋牧場表底層余流變化劇烈程度基本一致。

    圖6 日平均流矢量和逐時余流標(biāo)準(zhǔn)差分布圖

    注: 子圖a: 日平均流矢量; b—e: 橢圓長軸方向: 余流主軸方向; 橢圓長軸: 逐時余流主成分分量標(biāo)準(zhǔn)差; 橢圓短軸方向: 余流次軸方向; 橢圓短軸: 逐時余流次軸分量標(biāo)準(zhǔn)差

    2.5 影響機制

    為進一步探究4處海洋牧場的余流結(jié)構(gòu), 繪制了觀測期間平均的HYCOM 0 m水深流場與ERA5 10 m高度風(fēng)場。如圖7a所示, 在西北風(fēng)背景下, 山東半島北部海域近岸存在一條自西向東流幅逐漸變窄的流動, 在成山頭轉(zhuǎn)為南向, 為魯北沿岸流[13, 17]。觀測期間, 威海西港海洋牧場平均流與此流動特征較為一致, 煙臺安源海洋牧場、威海瑜泰海洋牧場、榮成楮島海洋牧場平均流與此流動存在顯著差異。煙臺安源海洋牧場平均流指向相鄰海域風(fēng)場、流場的右側(cè)(圖7a), 推測其平均流與風(fēng)場有關(guān), 在后文討論了日平均風(fēng)與日平均流的特征。威海瑜泰海洋牧場平均流為北向、離岸方向, 相近的威海西港海洋牧場則為東向。前者與威海外海東向的HYCOM表層流場存在差異, 后者則相近, 推測威海瑜泰海洋牧場流矢量變化除受魯北沿岸流影響外, 存在其他影響機制, 在后文討論日平均風(fēng)與日平均流的特征。威海西港海洋牧場表、底層余流特征基本一致, 說明表底層余流具有較為顯著的正壓性, 平均流與臨近海域HYCOM表層流矢量方向基本一致, 為東—東南向, 推測此海域余流與魯北沿岸流密切相關(guān)。榮成楮島海洋牧場位于桑溝灣南岸灣口處(圖1f), 桑溝灣外流矢量為南向[15], 但榮成楮島海洋牧場平均流由表至底均為東北向, 與灣外流場對流, 說明海水自灣南流出, 僅考慮灣口的水交換, 則海水自灣北流入, 在灣內(nèi)逆時針旋轉(zhuǎn)。推論與以往數(shù)值模擬研究結(jié)果較為一致[24]。

    逐日風(fēng)矢量如圖7b所示。所有牧場逐日風(fēng)矢量分布基本一致, 西北風(fēng)與西南風(fēng)交替出現(xiàn)。對比圖7b與圖6, 發(fā)現(xiàn)北風(fēng)增強時煙臺安源海洋牧場日平均流流速明顯增大。為進一步討論風(fēng)速對日平均流的影響, 將同一站位風(fēng)速和不同水層余流做相關(guān)分析。結(jié)果顯示, 在煙臺安源海洋牧場, 表底層經(jīng)向余流流速與同站經(jīng)向風(fēng)速相關(guān)系數(shù)大于0.6, 表層緯向余流流速與同站經(jīng)向風(fēng)速相關(guān)系數(shù)小于–0.4, 在威海瑜泰海洋牧場, 整層水體經(jīng)向余流流速與同站經(jīng)向風(fēng)速相關(guān)系數(shù)小于–0.6。圖8展示了煙臺安源海洋牧場表層日平均經(jīng)向流異常、緯向流異常與10 m日平均經(jīng)向風(fēng)異常, 以及威海瑜泰海洋牧場表層日平均經(jīng)向流異常與10 m日平均經(jīng)向風(fēng)異常的時間序列。煙臺安源海洋牧場經(jīng)向流與經(jīng)向風(fēng)正相關(guān), 緯向流與經(jīng)向風(fēng)負相關(guān), 在12月6日至12月9日、12月15日、12月27日日平均經(jīng)向風(fēng)正異常時, 日平均經(jīng)向流正異常, 日平均緯向流負異常, 12月5日、12月11日、12月17日、12月26日、12月30日經(jīng)向風(fēng)異常為極小值時, 日平均經(jīng)向流異常為極小值, 日平均緯向流異常為極大值, 即北風(fēng)增強時, 煙臺安源海洋牧場南向流、東向流增強。威海瑜泰海洋牧場經(jīng)向余流與經(jīng)向風(fēng)速負相關(guān), 如圖8c所示, 12月5日、12月11日、12月17日、12月26日、12月30日經(jīng)向風(fēng)異常為極小值時, 日平均經(jīng)向流異常為極大值, 經(jīng)向風(fēng)異常為極大值時, 日平均經(jīng)向流異常為極小值, 即北風(fēng)增強時, 威海瑜泰海洋牧場北向流增強。

    圖7 觀測期間0 m平均流場與10 m平均風(fēng)場

    注: 子圖a: 藍色箭頭: HYCOM再分析資料0 m流場; 紅色箭頭: ERA5 10 m風(fēng)場

    圖8 煙臺安源海洋牧場、威海瑜泰海洋牧場表層日平均流異常與10 m高度經(jīng)向風(fēng)速異常時間序列圖

    煙臺安源海洋牧場平均水深5.4 m, 水深淺, 經(jīng)向風(fēng)推動表底層海水向岸流動, 因而煙臺安源海洋牧場經(jīng)向余流與經(jīng)向風(fēng)正相關(guān); 經(jīng)向風(fēng)與岸線有夾角, 北風(fēng)推動表層海水沿岸向東流動, 因而在表層, 煙臺安源海洋牧場緯向余流與經(jīng)向風(fēng)負相關(guān)(圖9a)。威海瑜泰海洋牧場平均水深10.1 m, 水深較煙臺安源海洋牧場深, 地形相對陡峭。推測北風(fēng)時, 海面有南向流, 因經(jīng)向風(fēng)幾乎與岸線垂直, 海水向岸堆積并在岸邊產(chǎn)生下降流, 在海面以下產(chǎn)生離岸流, 即北向流; 南風(fēng)時, 海面為離岸流, 海面以下則產(chǎn)生上升流, 威海瑜泰海洋牧場表面以下為南向的補償流, 因而威海瑜泰海洋牧場經(jīng)向流與經(jīng)向風(fēng)負相關(guān)(圖9b)。

    3 結(jié)論

    通過對山東半島東北部4處海洋牧場2019年12月觀測的海流數(shù)據(jù)進行分析, 得到以下結(jié)論:

    1) 威海瑜泰海洋牧場、威海西港海洋牧場海域潮流類型為不正規(guī)半日潮, 煙臺安源海洋牧場、榮成楮島海洋牧場海域潮流類型介于正規(guī)半日潮和不正規(guī)半日潮之間。M2、S2半日分潮潮流在榮成楮島海洋牧場流速最大, M2分潮流速約為7~12 cm/s; O1、K1全日分潮潮流在威海瑜泰海洋牧場流速最大, K1分潮流速約為3~6 cm/s。各海洋牧場海域的主要分潮潮流具有最大流速平行岸線, 最大流速的大小隨深度減小的往復(fù)運動特征, 表明各海洋牧場海域潮流運動受地形邊界影響顯著。

    圖9 煙臺安源海洋牧場、威海瑜泰海洋牧場流場與經(jīng)向風(fēng)場相關(guān)關(guān)系

    注: 綠色箭頭: 風(fēng)場; 藍色箭頭: 流場

    2) 煙臺安源海洋牧場平均流大致垂直于岸線流向近岸, 流速約為0.9~1.7 cm/s, 威海瑜泰海洋牧場平均流大致垂直于岸線流向外海, 流速約為1.4~ 1.7 cm/s, 威海西港海洋牧場平均流為東南方向, 流速約為2.5~3.0 cm/s, 榮成楮島海洋牧場平均流為東北方向, 流速約為5.6~9.9 cm/s。

    3) 煙臺安源海洋牧場表底層日平均流經(jīng)向流與當(dāng)?shù)亟?jīng)向風(fēng)正相關(guān), 海水受北風(fēng)強迫在近岸堆積。威海瑜泰海洋牧場日平均流以經(jīng)向分量為主, 經(jīng)向分量主要為北向, 經(jīng)向分量與當(dāng)?shù)亟?jīng)向風(fēng)負相關(guān), 北風(fēng)強迫下, 威海瑜泰海洋牧場海面為南向流, 海水向岸堆積并在近岸產(chǎn)生下降流, 在海面下為北向的離岸流。威海西港海洋牧場表底層日平均流為東南向, 余流變化較為一致, 具有較為顯著的正壓性。榮成楮島海洋牧場日平均流為東北向, 表明海水沿桑溝灣南岸流出海灣, 在桑溝灣內(nèi)逆時針流動。

    [1] 王恩辰. 海洋牧場建設(shè)及其升級問題研究[D]. 青島: 中國海洋大學(xué), 2015. WANG Enchen. Research on the construction and upgrade of marine ranching[D]. Qingdao: Ocean University of China, 2015.

    [2] 李忠義, 林群, 李嬌, 等. 中國海洋牧場研究現(xiàn)狀與發(fā)展[J]. 水產(chǎn)學(xué)報, 2019, 43(9): 1870-1880. LI Zhongyi, LIN Qun, LI Jiao, et al. Present situation and future development of marine ranching construction in China[J]. Journal of Fisheries of China, 2019, 43(9): 1870-1880.

    [3] KITHEKA J U, OHOWA B O, MWASHOTE B M, et al. Water circulation dynamics, water column nutrients and plankton productivity in a well-flushed tropical bay in Kenya[J]. Journal of Sea Research, 1996, 35(4): 257- 268.

    [4] TANAKA K, KOMATSU K, ITOH S, et al. Baroclinic circulation and its high frequency variability in Otsuchi Bay on the Sanriku ria coast, Japan[J]. Journal of Ocea-nography, 2017, 73: 25-38.

    [5] LIN H Y, CHEN Z Z, HU J Y, et al. Impact of cage aquaculture on water exchange in Sansha Bay[J]. Continental Shelf Research, 2019, 188: 1-15.

    [6] 樊星, 魏皓, 原野, 等. 近岸典型養(yǎng)殖海區(qū)的潮流垂直結(jié)構(gòu)特征[J]. 中國海洋大學(xué)學(xué)報(自然科學(xué)版), 2009, 39(2): 181-186, 192. FAN Xing, WEI Hao, YUAN Ye, et al. The features of vertical structures of tidal current in a typical coastal mariculture area of China[J]. Periodical of Ocean University of China, 2009, 39(2): 181-186, 192.

    [7] 顏慧慧, 王鳳霞. 中國海洋牧場研究文獻綜述[J]. 科技廣場, 2016(6): 162-167. YAN Huihui, WANG Fengxia. Literature review on marine ranching in China[J]. Science Mosaic, 2016(6): 162-167.

    [8] 修日晨, 李繁華, 孔祥德. 渤海及北黃海潮流場的基本特征[J]. 海洋科學(xué), 1989, 13(5): 1-7. XIU Richen, LI Fanhua, KONG Xiangde. The basic characteristic of tidal flow field in Bohai Sea and North Huanghai Sea[J]. Marine Sciences, 1989, 13(5): 1-7.

    [9] 孫長青, 王學(xué)昌, 趙可勝. 煙臺北部近岸海域潮流數(shù)值模擬[J]. 青島海洋大學(xué)學(xué)報, 1994(S1): 42-49. SUN Changqing, WANG Xuechang, ZHAO Kesheng. Numerical computation of tidal current in the offshore area of northern Yantai[J]. Journal of Ocean University of Qingdao, 1994(S1): 42-49.

    [10] 姚蘭芳. 威海灣潮流的分布特征[J]. 海岸工程, 1990(Z1): 61-65. YAO Lanfang. Distributional features of the tidal current in the Weihai Bay[J]. Coastal Engineering, 1990(Z1): 61-65.

    [11] 李慶杰. 成山頭海域波浪能和潮流能時空分布特征[D]. 青島: 中國海洋大學(xué), 2013.LI Qingjie. Space-time distributional Characteristics of wave energy and tidal current energy in Chengshantou Sea Area[D]. Qingdao: Ocean University of China, 2013.

    [12] 黃祖珂, 黃磊. 潮汐原理與計算[M]. 青島: 中國海洋大學(xué)出版社, 2005: 209-238. HUANG Zuke, HUANG Lei. Tidal Theory and Calculation[M]. Qingdao: China Ocean University Press, 2005: 209-238.

    [13] 張志欣. 中國近海沿岸流及毗鄰流系的觀測與分析研究[D]. 青島: 中國海洋大學(xué), 2014.ZHANG Zhixin. Observation and analysis of the coas-tal current and its adjacent current system in the China offshore waters[D]. Qingdao: Ocean University of China, 2014.

    [14] LU J, QIAO F L, WANG X H, et al. A numerical study of transport dynamics and seasonal variability of the Yellow River sediment in the Bohai and Yellow seas[J]. Estuarine, Coastal and Shelf Science, 2011, 95(1): 39-51.

    [15] 王勇智, 張永強, 孫惠鳳. 山東半島東部海域懸浮體分布季節(jié)變化及其冬季輸送通量研究[J]. 沉積學(xué)報, 2019, 37(3): 541-549. WANG Yongzhi, ZHANG Yongqiang, SUN Huifeng. Seasonal variation of suspended matter distribution and flux in coastal waters of eastern Shandong Peninsula[J]. Acta Sedimentologica Sinica, 2019, 37(3): 541-549.

    [16] 宋軍, 姚志剛, 郭俊如, 等. 基于潮汐逆模型技術(shù)對渤黃海正壓M2分潮開邊界條件的優(yōu)化研究: Ⅱ. 潮汐特征、潮汐動力學(xué)及潮余流[J]. 海洋通報, 2016, 35(5): 545-553.SONG Jun, YAO Zhigang, GUO Junru, et al. Refining barotropic tide simulations for the Bohai and Yellow Seas using a hybrid data assimilation approach[J]. Marine Science Bulletin, 2016, 35(5): 545-553.

    [17] 韋欽勝, 于志剛, 冉祥濱, 等. 黃海西部沿岸流系特征分析及其對物質(zhì)輸運的影響[J]. 地球科學(xué)進展, 2011, 26(2): 145-156.WEI Qinsheng, YU Zhigang, RAN Xiangbin, et al. Cha-racteristics of the Western Coastal Current of the Yellow Sea and Its Impacts on Material Transportation[J]. Advances in Earth Science, 2011, 26(2): 145-156.

    [18] HERSBACH H, BELL W, BERRISFORD P, et al. Global reanalysis: goodbye ERA-Interim, hello ERA5[J]. ECMWF Newsletter, 2019, 159: 17-24.

    [19] PAWLOWIEZ R, BEARDSLEY B, LENTZ S, et al. Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE[J]. Computers & Geosciences, 2002, 28(8): 929-937.

    [20] 皇甫雪睿, 李毅, 胡日軍, 等. 煙臺西部近岸海域冬季懸浮泥沙濃度時空變化及其輸運特征[J]. 海洋地質(zhì)前沿, 2020, 36(5): 22-33. HUANGFU Xuerui, LI Yi, HU Rijun, et al. Temporal and spatial variation of suspended sediment concentration in the nearshore waters of western Yantai and its transport characterstics in winter[J]. Marine Geology Frontiers, 2020, 36(5): 22-33.

    [21] 費岳軍, 史軍強, 堵盤軍, 等. 冬季舟山外海定點實測海流資料分析[J]. 海洋通報, 2013, 32(6): 648-656. FEI Yuejun, SHI Junqiang, DU Panjun, et al. Analysis of the mooring current data near Zhoushan Islands[J]. Marine Science Bulletin, 2013, 32(6): 648-656.

    [22] KUNDU P K, ALLEN J S. Some Three-dimensional characteristics of low-frequency current fluctuations near the Oregon coast[J]. Journal of Physical Oceanography, 1976, 6(2): 181-199.

    [23] TANIGUCHI N, HUANG C F, ARAI M, et al. Variation of residual current in the Seto Inland Sea driven by sea level difference between the Bungo and Kii Channels[J]. Journal of Geophysical Research: Oceans, 2018, 123(4): 2921-2933.

    [24] 王麗霞, 石磊, 孫長青. 桑溝灣海域的潮流數(shù)值計算[J]. 青島海洋大學(xué)學(xué)報, 1994(S1): 77-83. WANG Lixia, SHI Lei, SUN Changqing. Tidal current computation for Sanggou Bay[J]. Journal of Ocean University of Qingdao, 1994(S1): 77-83.

    Spatial and temporal characteristics of sea current within the marine ranches in the northeast of Shandong Peninsula

    HU Luo-yu1, ZHAI Fang-guo1, LIU Zi-zhou1, LI Pei-liang2, GU Yan-zhen1, SUN Li-yuan3, LI Xin1, CHEN Dong2, JIANG Qing-yan2

    (1. College of Ocean and Atmosphere Science, Ocean University of China, Qingdao 266100, China; 2. Ocean College, Zhejiang University, Zhoushan 316021, China; 3. Shandong Hydrobios Resources Conservation and Management Center, Yantai 264003, China)

    To enrich the study of the hydrodynamic environment around the near-shore area of the Shandong Peninsula in winter, the structure of currents, such as the tidal and residual components, is investigated based on current observations within four marine ranches from December 3, 2019 to January 1, 2020, employing the methods of harmonic analysis, power spectra analysis, principal component analysis, and correlation analysis. Within four marine ranches, with the M2constituent that was significant and reciprocating, the directions of the maximum velocity vector of the main tidal constituents were basically parallel to the shoreline. The characteristics and influence mechanisms of the residual current varied with the area of marine ranches. The average current within the Anyuan Marine Ranch was southeastward and shoreward, with a velocity range of 0.9–1.7 cm/s.The main component of the daily mean current was the meridional current, which had a significant positive correlation with the wind in the same direction. While the northerly wind was driving, the residual current accumulated towards the southwest bank of the Anyuan Marine Ranch synchronously. The average current within the Yutai Marine Ranch was northward and offshore, with a velocity of 1.4–1.7 cm/s. The main component of the daily mean current was the meridional current, which had a significant negative correlation with the meridional wind. The daily average current showed that the northerly wind forced the seawater to accumulate to the south bank of the Yutai Marine Ranch and produced downflow near the coast. In addition, the offshore current flowed northward below the sea surface. The average current within the Xigang Marine Ranch was southeastward, with a velocity of 2.5–3.0 cm/s. The daily average current was barotropic at both the surface and bottom layers. The average current within the Chudao Marine Ranch was northeast, with a velocity of 5.6–9.9 cm/s. The daily average current was also northeast along the south bank of Sanggou Bay, indicating that the seawater flowed counterclockwise in the bay. The research results provided a basis for the further study of multi-scale temporal and spatial characteristics and influence mechanism of the dynamic environment in the near-shore area of the Shandong Peninsula.

    marine ranch; tidal currents; residual currents; coastal area of Shandong Province

    Aug. 11, 2020

    P717

    A

    1000-3096(2021)04-0001-12

    10.11759/hykx20200811003

    2020-08-11;

    2020-09-25

    國家自然科學(xué)基金項目(41776012); 國家重點研發(fā)計劃 (2019YFD0901305); 浙江省重點研發(fā)計劃(2020C03012); 三亞崖州灣科技城管理局重大科技項目(SKJC-KJ-2019KY03)

    [National Science Foundation of China, No. 41776012; National Key Research & Development Program of China, No. 2019YFD0901305; Key Research & Development Project of Zhejiang Province, No. 2020C03012; Major Science and Technology Project of Sanya YZBSTC, No. SKJC-KJ-2019KY03]

    胡羅煜(1998—), 男, 湖北武漢人, 碩士研究生, 研究方向: 物理海洋學(xué), E-mail: huluoyu@stu.ouc.edu.cn; 劉子洲(1987—),通信作者, 實驗師, 主要從事海洋調(diào)查及近海海洋環(huán)境問題的教學(xué)與研究, 電話: 13708954171, E-mail:lzz2013@ouc.edu.cn

    (本文編輯: 叢培秀)

    猜你喜歡
    余流分潮經(jīng)向
    浸膠帆布經(jīng)向剛度影響因素分析
    橡膠科技(2022年5期)2022-07-20 02:24:04
    大亞灣雙峰水位的形成條件及準(zhǔn)調(diào)和分量應(yīng)用的分析
    與南亞高壓相聯(lián)的歐亞大陸-印度洋經(jīng)向環(huán)流
    山東鄰海長周期分潮對深度基準(zhǔn)面的影響分析
    2018年8月大氣環(huán)流中水汽經(jīng)向輸送特征
    有關(guān)副熱帶太平洋對ENSO影響研究的綜述
    基于長期觀測的遼東灣口東部海域水動力特征研究
    夏秋季泉州灣中部海域潮流和余流的變化特征
    基于走航ADCP資料的廈門內(nèi)灣東西口門海域潮流與潮致余流特征分析
    考慮內(nèi)潮耗散的南海M2分潮伴隨同化數(shù)值模擬
    少妇粗大呻吟视频| av一本久久久久| 亚洲自偷自拍图片 自拍| 久久久久久人人人人人| 精品人妻熟女毛片av久久网站| 亚洲激情五月婷婷啪啪| 欧美少妇被猛烈插入视频| 日韩免费高清中文字幕av| 免费在线观看影片大全网站 | 欧美精品人与动牲交sv欧美| 亚洲精品国产一区二区精华液| 久久精品国产综合久久久| 欧美亚洲日本最大视频资源| 岛国毛片在线播放| 人人妻人人澡人人看| 校园人妻丝袜中文字幕| 国产欧美日韩一区二区三区在线| 亚洲av日韩精品久久久久久密 | 一边摸一边抽搐一进一出视频| 精品亚洲成国产av| 99国产综合亚洲精品| 久久久精品国产亚洲av高清涩受| 好男人视频免费观看在线| 午夜视频精品福利| 日韩熟女老妇一区二区性免费视频| 麻豆国产av国片精品| 久久影院123| 一级a爱视频在线免费观看| 男人添女人高潮全过程视频| 一级,二级,三级黄色视频| 亚洲综合色网址| 亚洲国产最新在线播放| 在线观看一区二区三区激情| 一级毛片黄色毛片免费观看视频| 婷婷丁香在线五月| 侵犯人妻中文字幕一二三四区| 亚洲国产看品久久| 汤姆久久久久久久影院中文字幕| 午夜免费鲁丝| 丝瓜视频免费看黄片| 岛国毛片在线播放| 国产97色在线日韩免费| 在现免费观看毛片| 美女中出高潮动态图| 亚洲国产精品一区二区三区在线| 黑人猛操日本美女一级片| 午夜激情av网站| 亚洲精品久久成人aⅴ小说| 交换朋友夫妻互换小说| 精品一区在线观看国产| 超碰97精品在线观看| 亚洲,欧美精品.| 欧美日韩成人在线一区二区| 国精品久久久久久国模美| 国产高清视频在线播放一区 | 观看av在线不卡| 免费在线观看日本一区| 韩国精品一区二区三区| 中文字幕人妻丝袜制服| 熟女少妇亚洲综合色aaa.| 老司机影院成人| 日韩视频在线欧美| 中国美女看黄片| 激情五月婷婷亚洲| 悠悠久久av| 看免费成人av毛片| 女警被强在线播放| 国产精品久久久久成人av| 精品亚洲成国产av| 色精品久久人妻99蜜桃| 久久人人爽av亚洲精品天堂| 好男人电影高清在线观看| 久久女婷五月综合色啪小说| 国产成人欧美在线观看 | 国产精品成人在线| 高清不卡的av网站| 国产成人系列免费观看| 51午夜福利影视在线观看| 又粗又硬又长又爽又黄的视频| 国产老妇伦熟女老妇高清| 国产淫语在线视频| 在线观看免费日韩欧美大片| 少妇猛男粗大的猛烈进出视频| 成人国产一区最新在线观看 | 97精品久久久久久久久久精品| 另类亚洲欧美激情| 男人操女人黄网站| 高清不卡的av网站| 日本欧美视频一区| 国产成人影院久久av| 青春草亚洲视频在线观看| 一级黄色大片毛片| 欧美日本中文国产一区发布| 女人精品久久久久毛片| 国产精品一区二区精品视频观看| 亚洲人成网站在线观看播放| 69精品国产乱码久久久| 欧美人与性动交α欧美软件| 国产欧美日韩一区二区三 | 少妇粗大呻吟视频| 欧美日韩综合久久久久久| 亚洲av成人精品一二三区| 咕卡用的链子| 欧美变态另类bdsm刘玥| 美女主播在线视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲五月色婷婷综合| 蜜桃在线观看..| 欧美激情极品国产一区二区三区| 一区在线观看完整版| 免费一级毛片在线播放高清视频 | 悠悠久久av| 国产免费现黄频在线看| www.精华液| 亚洲熟女毛片儿| 欧美日本中文国产一区发布| 可以免费在线观看a视频的电影网站| 国产精品久久久久成人av| 日韩精品免费视频一区二区三区| 好男人在线观看高清免费视频 | cao死你这个sao货| 啦啦啦观看免费观看视频高清| 亚洲专区中文字幕在线| 99国产精品一区二区蜜桃av| 国产又黄又爽又无遮挡在线| av中文乱码字幕在线| 亚洲七黄色美女视频| 欧美国产日韩亚洲一区| 成人av一区二区三区在线看| 欧美日韩乱码在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲电影在线观看av| 哪里可以看免费的av片| 午夜福利成人在线免费观看| 亚洲一区二区三区色噜噜| videosex国产| svipshipincom国产片| 久久久久久免费高清国产稀缺| 国产亚洲欧美在线一区二区| 国产成人啪精品午夜网站| 精品不卡国产一区二区三区| 听说在线观看完整版免费高清| 老司机午夜福利在线观看视频| 黄色成人免费大全| 淫秽高清视频在线观看| 亚洲第一电影网av| 人人妻,人人澡人人爽秒播| 国产高清有码在线观看视频 | 在线国产一区二区在线| 国产免费av片在线观看野外av| 午夜精品在线福利| 国产精品国产高清国产av| 亚洲激情在线av| 麻豆成人午夜福利视频| 手机成人av网站| 欧美黑人欧美精品刺激| 国产v大片淫在线免费观看| 国产精品亚洲一级av第二区| 国内毛片毛片毛片毛片毛片| 亚洲天堂国产精品一区在线| 亚洲成人久久性| 精品国产乱码久久久久久男人| √禁漫天堂资源中文www| 国产成人欧美| 免费人成视频x8x8入口观看| 国产激情欧美一区二区| 亚洲精品久久成人aⅴ小说| 亚洲色图 男人天堂 中文字幕| 黄色毛片三级朝国网站| 久久精品夜夜夜夜夜久久蜜豆 | 十分钟在线观看高清视频www| 国产乱人伦免费视频| 长腿黑丝高跟| 欧洲精品卡2卡3卡4卡5卡区| 国产成人欧美| av超薄肉色丝袜交足视频| 国产99久久九九免费精品| 波多野结衣av一区二区av| 久久久久国内视频| 欧美绝顶高潮抽搐喷水| 国产一区二区激情短视频| 久久久久国产一级毛片高清牌| 日韩中文字幕欧美一区二区| 日日干狠狠操夜夜爽| 欧美久久黑人一区二区| 97人妻精品一区二区三区麻豆 | 国产精品99久久99久久久不卡| 国产av不卡久久| 国产单亲对白刺激| 国内少妇人妻偷人精品xxx网站 | 搡老妇女老女人老熟妇| 中文字幕av电影在线播放| 搡老岳熟女国产| 成年人黄色毛片网站| 亚洲国产日韩欧美精品在线观看 | 男男h啪啪无遮挡| 国产一区在线观看成人免费| 美女 人体艺术 gogo| 午夜激情av网站| 黄色视频不卡| 中文字幕最新亚洲高清| 2021天堂中文幕一二区在线观 | 99国产极品粉嫩在线观看| 亚洲精华国产精华精| 精品国产乱子伦一区二区三区| 国产成人系列免费观看| tocl精华| 日本成人三级电影网站| 最近最新中文字幕大全免费视频| 亚洲国产欧美网| 欧美黑人精品巨大| 国产97色在线日韩免费| 12—13女人毛片做爰片一| 亚洲aⅴ乱码一区二区在线播放 | 亚洲国产精品999在线| 18禁裸乳无遮挡免费网站照片 | 亚洲精品粉嫩美女一区| 一本精品99久久精品77| 中文字幕人成人乱码亚洲影| 国产国语露脸激情在线看| 欧美日韩精品网址| 午夜亚洲福利在线播放| 精品久久久久久久人妻蜜臀av| 日韩精品免费视频一区二区三区| 久久这里只有精品19| 亚洲av五月六月丁香网| 亚洲人成77777在线视频| 欧美中文综合在线视频| 韩国精品一区二区三区| 欧美在线一区亚洲| 在线国产一区二区在线| 精品乱码久久久久久99久播| 婷婷丁香在线五月| 99久久久亚洲精品蜜臀av| 国产黄色小视频在线观看| 国产高清videossex| ponron亚洲| 国产成人精品久久二区二区91| 欧美久久黑人一区二区| 国产不卡一卡二| 日韩有码中文字幕| 欧美成人午夜精品| а√天堂www在线а√下载| 制服人妻中文乱码| 在线视频色国产色| 久9热在线精品视频| 不卡av一区二区三区| 岛国视频午夜一区免费看| 国产日本99.免费观看| 一级黄色大片毛片| 男女之事视频高清在线观看| 99国产精品99久久久久| 欧美激情 高清一区二区三区| 国产高清有码在线观看视频 | 欧美国产日韩亚洲一区| 午夜福利18| 这个男人来自地球电影免费观看| 搞女人的毛片| 视频区欧美日本亚洲| 国产成人啪精品午夜网站| 亚洲专区中文字幕在线| 一本久久中文字幕| 国语自产精品视频在线第100页| 黄色视频不卡| 欧美中文综合在线视频| 精品久久久久久,| 国产色视频综合| 国产精品亚洲美女久久久| 亚洲男人天堂网一区| 久久精品国产99精品国产亚洲性色| 精品国产超薄肉色丝袜足j| 免费观看人在逋| 日本在线视频免费播放| 人人妻,人人澡人人爽秒播| 又大又爽又粗| 欧美不卡视频在线免费观看 | 国产精品久久久久久人妻精品电影| 丁香六月欧美| 精品乱码久久久久久99久播| 亚洲国产中文字幕在线视频| 亚洲精品一区av在线观看| 国产精华一区二区三区| 又紧又爽又黄一区二区| 69av精品久久久久久| 精品国产国语对白av| 天堂动漫精品| 黄片播放在线免费| 精品国产美女av久久久久小说| 国产成人欧美在线观看| 久久久国产欧美日韩av| 一进一出抽搐动态| 999久久久国产精品视频| 一区二区三区精品91| 激情在线观看视频在线高清| 久久精品国产亚洲av高清一级| 亚洲av成人一区二区三| 国产精品爽爽va在线观看网站 | 国产蜜桃级精品一区二区三区| www国产在线视频色| 成人国产一区最新在线观看| 黄色片一级片一级黄色片| 亚洲 欧美 日韩 在线 免费| 久久精品成人免费网站| 亚洲精品一卡2卡三卡4卡5卡| 国产黄色小视频在线观看| 色播在线永久视频| 90打野战视频偷拍视频| 亚洲国产欧美网| 亚洲精品国产精品久久久不卡| 男女那种视频在线观看| 国产精品久久久久久精品电影 | 免费高清视频大片| 一级毛片女人18水好多| 久久精品影院6| xxx96com| 亚洲av美国av| 级片在线观看| 久久中文看片网| 欧美国产日韩亚洲一区| 久久99热这里只有精品18| 特大巨黑吊av在线直播 | 在线播放国产精品三级| 色综合亚洲欧美另类图片| 成人亚洲精品av一区二区| 91成人精品电影| 老司机午夜十八禁免费视频| 久久欧美精品欧美久久欧美| 99精品久久久久人妻精品| 国产一级毛片七仙女欲春2 | 亚洲国产日韩欧美精品在线观看 | 国产免费av片在线观看野外av| 日韩高清综合在线| 成人一区二区视频在线观看| 老汉色av国产亚洲站长工具| 免费在线观看日本一区| 国产精品九九99| 我的亚洲天堂| 一二三四社区在线视频社区8| 亚洲欧洲精品一区二区精品久久久| 男女下面进入的视频免费午夜 | 久久精品人妻少妇| 午夜精品久久久久久毛片777| 精品国产亚洲在线| 国产欧美日韩一区二区精品| 久久久久久九九精品二区国产 | 99国产精品99久久久久| or卡值多少钱| 好看av亚洲va欧美ⅴa在| 黄色视频不卡| 人人妻人人看人人澡| 女性生殖器流出的白浆| 色播在线永久视频| 香蕉丝袜av| 国产成+人综合+亚洲专区| av在线天堂中文字幕| 色综合婷婷激情| 成人亚洲精品一区在线观看| 好看av亚洲va欧美ⅴa在| 黄色视频,在线免费观看| 久久精品国产综合久久久| 精品久久久久久久毛片微露脸| 麻豆一二三区av精品| 国产精品日韩av在线免费观看| 亚洲精品av麻豆狂野| 黄色视频不卡| 99久久99久久久精品蜜桃| 久久久久久久午夜电影| 中文资源天堂在线| 久久人妻福利社区极品人妻图片| 国产精品乱码一区二三区的特点| 夜夜爽天天搞| 欧美日韩精品网址| 久久中文字幕人妻熟女| 麻豆一二三区av精品| 国产亚洲欧美在线一区二区| 人人澡人人妻人| 亚洲国产精品999在线| 两人在一起打扑克的视频| 亚洲美女黄片视频| 亚洲第一欧美日韩一区二区三区| 国产不卡一卡二| 亚洲男人的天堂狠狠| 日本 av在线| 欧美国产日韩亚洲一区| 亚洲男人天堂网一区| 午夜免费鲁丝| 亚洲无线在线观看| 久久精品亚洲精品国产色婷小说| 大型av网站在线播放| 亚洲色图av天堂| 1024视频免费在线观看| 亚洲中文字幕一区二区三区有码在线看 | 丰满的人妻完整版| 国产aⅴ精品一区二区三区波| 国产精品免费一区二区三区在线| 日韩大尺度精品在线看网址| 视频区欧美日本亚洲| 色综合婷婷激情| 99国产综合亚洲精品| 国产免费男女视频| 啪啪无遮挡十八禁网站| 中文亚洲av片在线观看爽| 禁无遮挡网站| 999久久久国产精品视频| 非洲黑人性xxxx精品又粗又长| 女人高潮潮喷娇喘18禁视频| 91麻豆av在线| 国产精品久久久久久精品电影 | 国产精品精品国产色婷婷| 韩国av一区二区三区四区| 久久人妻福利社区极品人妻图片| 亚洲午夜精品一区,二区,三区| 亚洲精华国产精华精| 国产亚洲精品久久久久久毛片| 一本久久中文字幕| 日韩免费av在线播放| 黄色片一级片一级黄色片| 人人妻人人澡欧美一区二区| 成熟少妇高潮喷水视频| 国产精品,欧美在线| 50天的宝宝边吃奶边哭怎么回事| 久久久久久大精品| 欧美 亚洲 国产 日韩一| 狠狠狠狠99中文字幕| 婷婷精品国产亚洲av| 亚洲五月天丁香| 中文亚洲av片在线观看爽| 国产欧美日韩一区二区精品| 亚洲真实伦在线观看| 久久久久国产一级毛片高清牌| 免费在线观看成人毛片| 久久久久亚洲av毛片大全| 女人被狂操c到高潮| 午夜a级毛片| 美女扒开内裤让男人捅视频| 精品乱码久久久久久99久播| 精品高清国产在线一区| 免费女性裸体啪啪无遮挡网站| 老熟妇乱子伦视频在线观看| 欧美成狂野欧美在线观看| 国内少妇人妻偷人精品xxx网站 | bbb黄色大片| 最好的美女福利视频网| 法律面前人人平等表现在哪些方面| 免费av毛片视频| 九色国产91popny在线| 变态另类丝袜制服| 好看av亚洲va欧美ⅴa在| 69av精品久久久久久| 亚洲熟女毛片儿| 99精品久久久久人妻精品| 女人高潮潮喷娇喘18禁视频| 法律面前人人平等表现在哪些方面| 在线免费观看的www视频| 国产激情欧美一区二区| 国产黄色小视频在线观看| 黄色成人免费大全| 亚洲人成电影免费在线| 91成人精品电影| www.精华液| 亚洲成av片中文字幕在线观看| 成人精品一区二区免费| 国产成人精品无人区| 亚洲三区欧美一区| 国产成人av教育| 亚洲国产精品999在线| 亚洲第一电影网av| 女生性感内裤真人,穿戴方法视频| 久久久久免费精品人妻一区二区 | av电影中文网址| a级毛片在线看网站| 国产成人欧美| 久久精品国产99精品国产亚洲性色| 国产蜜桃级精品一区二区三区| 亚洲欧美日韩高清在线视频| 国内毛片毛片毛片毛片毛片| 亚洲五月色婷婷综合| av福利片在线| 精华霜和精华液先用哪个| 中文在线观看免费www的网站 | 午夜激情福利司机影院| 巨乳人妻的诱惑在线观看| 中文字幕最新亚洲高清| 一个人免费在线观看的高清视频| 成年版毛片免费区| 99精品久久久久人妻精品| 日韩成人在线观看一区二区三区| 天堂影院成人在线观看| 露出奶头的视频| 精品福利观看| 欧美日韩一级在线毛片| 国产在线观看jvid| 欧美久久黑人一区二区| 久久中文字幕一级| 午夜两性在线视频| 国产日本99.免费观看| 一区二区日韩欧美中文字幕| 国产精品电影一区二区三区| 欧美成人性av电影在线观看| 国产一级毛片七仙女欲春2 | 免费在线观看完整版高清| 亚洲av成人一区二区三| 午夜福利在线观看吧| 妹子高潮喷水视频| 国产精品亚洲av一区麻豆| 久久精品国产清高在天天线| 麻豆一二三区av精品| 大型av网站在线播放| 老汉色∧v一级毛片| 亚洲成国产人片在线观看| 国产精品,欧美在线| 国产亚洲精品综合一区在线观看 | 高清在线国产一区| a级毛片a级免费在线| 国产精品99久久99久久久不卡| 精品熟女少妇八av免费久了| 国产精品香港三级国产av潘金莲| 国产精品电影一区二区三区| cao死你这个sao货| 国内少妇人妻偷人精品xxx网站 | 狂野欧美激情性xxxx| 免费电影在线观看免费观看| 久久精品国产综合久久久| 麻豆av在线久日| 18禁观看日本| 熟女少妇亚洲综合色aaa.| 黄色视频不卡| 亚洲色图av天堂| 女警被强在线播放| 久久婷婷人人爽人人干人人爱| 国产1区2区3区精品| 满18在线观看网站| 欧美一区二区精品小视频在线| 精品日产1卡2卡| 99在线人妻在线中文字幕| 免费在线观看视频国产中文字幕亚洲| 香蕉国产在线看| 成人三级黄色视频| 中文字幕久久专区| 满18在线观看网站| 精品电影一区二区在线| 99热这里只有精品一区 | av天堂在线播放| 久久久久精品国产欧美久久久| 午夜福利免费观看在线| 女同久久另类99精品国产91| 国产视频内射| 在线免费观看的www视频| 国产成+人综合+亚洲专区| 欧美乱妇无乱码| 18禁黄网站禁片免费观看直播| 久久青草综合色| 欧美zozozo另类| 亚洲专区中文字幕在线| 国产精品永久免费网站| 亚洲狠狠婷婷综合久久图片| 夜夜看夜夜爽夜夜摸| 国产区一区二久久| 88av欧美| 女性被躁到高潮视频| 欧美日韩福利视频一区二区| 午夜福利一区二区在线看| 亚洲国产日韩欧美精品在线观看 | 国产男靠女视频免费网站| 精品第一国产精品| 亚洲欧美一区二区三区黑人| 国产爱豆传媒在线观看 | 亚洲成人免费电影在线观看| 久久精品国产99精品国产亚洲性色| 十八禁网站免费在线| 亚洲最大成人中文| 日本一区二区免费在线视频| 丰满的人妻完整版| 亚洲专区字幕在线| 麻豆国产av国片精品| 欧美黑人精品巨大| 欧美乱码精品一区二区三区| 一本一本综合久久| 巨乳人妻的诱惑在线观看| 久久国产亚洲av麻豆专区| 久久久久国内视频| 国产精品久久久av美女十八| 日日夜夜操网爽| 中国美女看黄片| 无遮挡黄片免费观看| 亚洲 欧美一区二区三区| 国产精品久久久久久人妻精品电影| 亚洲中文字幕日韩| 国产成人av教育| 国内揄拍国产精品人妻在线 | 免费av毛片视频| 国产精品电影一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 每晚都被弄得嗷嗷叫到高潮| 丰满人妻熟妇乱又伦精品不卡| 久久人妻福利社区极品人妻图片| 亚洲国产欧美日韩在线播放| 亚洲第一电影网av| 黄频高清免费视频| 日韩免费av在线播放| 欧美大码av| 日本a在线网址| www.自偷自拍.com| 黄片播放在线免费| 老鸭窝网址在线观看| 午夜免费成人在线视频| 黄色a级毛片大全视频| 国产精品 国内视频| 波多野结衣av一区二区av| 免费一级毛片在线播放高清视频| 俺也久久电影网| 久久狼人影院| 女人爽到高潮嗷嗷叫在线视频| 国产真实乱freesex| 国产91精品成人一区二区三区|