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

    大沽河口水位變化特征分析

    2021-02-03 13:22:16賴秋宇劉子洲翟方國(guó)顧艷鎮(zhèn)姜慶巖
    海洋科學(xué) 2021年1期
    關(guān)鍵詞:大沽河口水深

    賴秋宇, 劉子洲, 翟方國(guó), 顧艷鎮(zhèn), 陳 棟, 姜慶巖

    大沽河口水位變化特征分析

    賴秋宇1, 劉子洲1, 翟方國(guó)1, 顧艷鎮(zhèn)1, 陳 棟2, 姜慶巖2

    (1. 中國(guó)海洋大學(xué) 海洋與大氣學(xué)院, 山東 青島 266100; 2. 浙江大學(xué) 海洋學(xué)院, 浙江 舟山 316021)

    本文基于2016年11月8日—2019年6月29日大沽河入??谒钣^測(cè)數(shù)據(jù)分析了大沽河河口水位變化特征, 并結(jié)合風(fēng)場(chǎng)、降水量、衛(wèi)星高度計(jì)融合產(chǎn)品資料對(duì)其影響因素展開了討論。結(jié)果表明: 1) 大沽河口水位變化由潮汐過(guò)程主導(dǎo), 潮汐類型為正規(guī)半日潮, M2分潮占主導(dǎo); 2) 余水位在2017年7月—2019年1月存在周期約為110—150天的顯著季節(jié)內(nèi)變化, 主要受到緯向風(fēng)的影響, 監(jiān)測(cè)系統(tǒng)處在大沽河入??谖靼? 東向(西向)風(fēng)將驅(qū)動(dòng)水體向東(西)輸運(yùn), 導(dǎo)致西岸監(jiān)測(cè)系統(tǒng)處水量減少(增加), 從而觀測(cè)到余水位下降(上升); 3) 觀測(cè)期間, 余水位存在顯著下降趨勢(shì), 約為–0.53×10–2m/月, 主要受到大沽河流域降水量減少的影響。

    大沽河口; 水位變化; 緯向風(fēng); 降水

    在全球氣候變暖這一氣候變化背景下衍生出的全球海平面上升問(wèn)題[1], 其將增加感潮河段(特別是河口三角洲和地勢(shì)低洼區(qū)域)潮位, 導(dǎo)致在極端天氣條件(臺(tái)風(fēng)、強(qiáng)溫帶氣旋、寒潮等)下沿海城市、港口與航道對(duì)風(fēng)暴潮災(zāi)害抵御能力下降, 直接影響沿海地區(qū)人民生活水平和制約社會(huì)經(jīng)濟(jì)發(fā)展[2-5]。青島市海岸帶區(qū)域人口密集導(dǎo)致其遭受風(fēng)暴潮災(zāi)害的脆弱性較高[6-7], 故有必要建設(shè)完善的海洋環(huán)境監(jiān)測(cè)系統(tǒng), 加強(qiáng)對(duì)海洋災(zāi)害性事件的預(yù)警預(yù)報(bào), 利用觀測(cè)水位時(shí)間序列是以往學(xué)者研究中國(guó)近海海平面變化的主要研究手段之一[8-10]。

    大沽河, 位于36°10′—37°12′N, 120°03′—120°25′E, 發(fā)源于山東省招遠(yuǎn)市境內(nèi), 于膠州灣馬頭村南注入膠州灣, 干流全長(zhǎng)179.9 km, 是膠州灣入海徑流最大的河流[11], 被稱作青島市的“母親河”。大沽河流域處在華北暖溫帶季風(fēng)氣候區(qū), 氣候具有夏季炎熱多雨, 冬季寒冷干燥的特點(diǎn), 流域降水集中于每年6—9月的汛期, 7、8月的降水量在全年總降水中占主導(dǎo), 且呈現(xiàn)下降趨勢(shì)[12], 以往學(xué)者對(duì)大沽河口的研究多集中在生態(tài)環(huán)境[13-16]、水質(zhì)[17-18]、水沙輸運(yùn)[19-20]、地下水交換[21-22]等方面。關(guān)于大沽河口水動(dòng)力過(guò)程的研究相對(duì)較少。尹則高等[23]通過(guò)建立潮流和徑流作用下大沽河口水動(dòng)力數(shù)學(xué)模型, 研究了潮位、潮流變化在不同徑流條件下的響應(yīng)。李兆欽等[24]在討論潮汐過(guò)程對(duì)大沽河口底層溶解氧濃度變化時(shí), 基于2016年11月6日—2017年6月20日水深觀測(cè)資料指出了大沽河口潮汐類型為正規(guī)半日潮。受限于缺乏觀測(cè)資料和觀測(cè)資料的時(shí)間長(zhǎng)度, 以往學(xué)者未能揭示大沽河口水位較低頻時(shí)間尺度變化, 對(duì)其變化趨勢(shì)和影響因素討論也較少, 故有必要利用更長(zhǎng)觀測(cè)資料做進(jìn)一步研究討論。

    本文基于放置于大沽河口的監(jiān)測(cè)系統(tǒng)在2016年11月8日—2019年6月29日期間的水深觀測(cè)數(shù)據(jù), 分析了大沽河口水位變化特征及其變化趨勢(shì), 并結(jié)合風(fēng)場(chǎng)、降雨量、外海水位等數(shù)據(jù)對(duì)影響水位變化的原因做出了討論。

    1 監(jiān)測(cè)系統(tǒng)與數(shù)據(jù)集

    1.1 監(jiān)測(cè)系統(tǒng)概況

    2016年11月在青島市海洋環(huán)境監(jiān)測(cè)系統(tǒng)項(xiàng)目下建設(shè)的一套海洋環(huán)境有纜在線監(jiān)測(cè)系統(tǒng), 位于大沽河入??跂|北—西南向河道西岸, 坐標(biāo)為36°11′44.15″N,120°7′0.57″E, 具體位置如圖1中紅色五角星所示。溫度、鹽度、pH、溶解氧、濁度、葉綠素等生態(tài)要素?cái)?shù)據(jù)由系統(tǒng)搭載的加拿大AML公司Metric X 50160水質(zhì)儀觀測(cè)得到, 觀測(cè)間隔為1 min。壓力、流速、流向、波浪數(shù)據(jù)由聲學(xué)波浪流速剖面儀(Acoustic Wave and Current Profiler, AWAC)觀測(cè)得到, 觀測(cè)間隔為1 min, 壓力的精度為0.01%, 分辨率為10–8。

    1.2 數(shù)據(jù)

    用于分析水位變化的水深數(shù)據(jù)由壓力觀測(cè)數(shù)據(jù)根據(jù)Fofonoff和Millard公式[25]換算得到, 時(shí)間間隔為1 min, 本文選取2016年11月8日—2019年6月29日期間數(shù)據(jù), 剔除缺測(cè)天數(shù), 有效數(shù)據(jù)長(zhǎng)度為931 d。

    為了討論影響大沽河口水位變化的影響因素, 本文還采用以下數(shù)據(jù): 由歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium-Range Weather Forecasts, ECMWF)提供的天平均ERA-Interim再分析資料[26], 主要使用經(jīng)向風(fēng)應(yīng)力(Meridional Wind Stress, MWS)和緯向風(fēng)應(yīng)力(Zonal Wind Stress, ZWS)數(shù)據(jù), 空間分辨率為0.125°×0.125°, 時(shí)間范圍為2016年1月1日—2019年7月1日, 利用雙線性插值法, 可以獲得觀測(cè)系統(tǒng)處的風(fēng)應(yīng)力時(shí)間序列; 歐盟哥白尼海洋環(huán)境監(jiān)測(cè)服務(wù)中心 (E.U. Copernicus Marine Environment Monitoring Se-rvice, CMEMS) Unidata Dataset Discovery v1.0 (http://marine.copernicus.eu/web/27-service-commitments- and-licence.php)多衛(wèi)星(Altika Drifting Phase, Cryosat-2, Haiyang-2A, OSTM/Jason-2)融合格點(diǎn)數(shù)據(jù)集。主要結(jié)果采用海平面異常(Sea Level Anomaly, SLA)數(shù)據(jù), 時(shí)間分辨率為1天, 空間分辨率0.125°×0.125°, 時(shí)間范圍為2016年1月1日—2019年5月13日; 降水?dāng)?shù)據(jù)來(lái)自亞洲—太平洋數(shù)據(jù)研究中心(Asia Pacific Data Resea-rch Center, APDRC)(http://apdrc.soest.hawaii. edu/data/data.php)的Climate Prediction Center (CPC)全球月氣候格點(diǎn)數(shù)據(jù)集[27-28], 該數(shù)據(jù)集由美國(guó)國(guó)家海洋大氣管理局(National Oceanic and Atmos-phere Administration station, NOAA)的氣候預(yù)報(bào)中心提供。數(shù)據(jù)通過(guò)插值第2版全球歷史氣候數(shù)據(jù)網(wǎng)和異常氣候監(jiān)測(cè)系統(tǒng)站點(diǎn)觀測(cè)數(shù)據(jù), 并考慮時(shí)空變化和地形影響, 空間分辨率為0.5°×0.5°, 能夠反映我國(guó)降水的時(shí)空特征[29], 時(shí)間跨度為2016年1月—2019年12月。本文選取包括山東半島在內(nèi)大范圍降水?dāng)?shù)據(jù), 以獲取降水分布特征及長(zhǎng)期降水趨勢(shì)。

    2 研究結(jié)果

    觀測(cè)期間, 大沽河口觀測(cè)系統(tǒng)處水深變化范圍為2.50~7.17 m, 平均水深為4.62 m。如圖2(a)所示, 水深呈現(xiàn)出多種時(shí)間尺度的變化, 小時(shí)平均水深變化顯示出明顯的高頻潮汐信號(hào), 大小潮過(guò)程顯著; 天平均水深變化也存在多時(shí)間尺度變化; 在季節(jié)尺度上, 水深自每年2月份開始增大, 至同年8月份上升至峰值, 于次年1月下降至谷值, 這表明水深的季節(jié)變化與流域豐水期和枯水期的交替變化相一致, 豐水期水深較大, 枯水期水深較小; 通過(guò)比較各月平均水深可以發(fā)現(xiàn)水深還呈現(xiàn)一定下降趨勢(shì)。綜上所述, 河口處水深變化在一定程度上受到入海徑流量和膠州灣潮汐的影響。水深時(shí)間序列減去觀測(cè)期間內(nèi)平均水深可得到水深異常, 近似等于水位異常。為獲取水位異常的顯著變化信號(hào), 本文選取整個(gè)觀測(cè)期間小時(shí)平均水位異常時(shí)間序列進(jìn)行功率譜分析[30], 結(jié)果如圖2(b)所示, 水位異常功率譜存在多個(gè)顯著峰值, 其中最大譜值對(duì)應(yīng)周期為12.42小時(shí), 對(duì)應(yīng)M2分潮的周期; 在4、6、12、24小時(shí)附近譜值較顯著, 分別與淺水分潮族、半日分潮族、全日分潮族的周期相近。上述結(jié)果表明, 大沽河口水位變化以潮信號(hào)為主。

    圖1 大沽河口外地形

    注: 圖中紅色五角星為監(jiān)測(cè)系統(tǒng)的位置; 藍(lán)色實(shí)線為大沽河河道; 水深圖像由海圖數(shù)據(jù)繪制; P點(diǎn)為CPC降水?dāng)?shù)據(jù)其中一個(gè)站點(diǎn)

    圖2 2016年11月8日至2019年6月29日小時(shí)平均、天平均、月平均水深時(shí)間序列(a), 小時(shí)平均水位異常功率譜密度(b)

    2.1 潮汐

    2.2 季節(jié)內(nèi)變化

    觀測(cè)水位異常與回報(bào)水位異常的差異可以表征余水位(Residual Water Level, RWL)異常變化[31], 圖4(a)展示了觀測(cè)期間RWL異常的天平均時(shí)間序列, 其變化范圍為–0.62~0.40 m, 標(biāo)準(zhǔn)差為0.07 m。本文通過(guò)Morlet小波變換方法[32], 分析RWL異常時(shí)間序列的顯著變化周期, 由圖4(b)展示結(jié)果可知, RWL異常在2017年7月—2019年1月存在110~150天的顯著季節(jié)內(nèi)變化。前人研究[31, 33]指出: 河口水位除受到外海潮汐影響外, 還受到風(fēng)速、風(fēng)向的影響。

    圖3 觀測(cè)水位異常與回報(bào)水位異常時(shí)間序列

    表1 大沽河口主要分潮調(diào)和常數(shù)

    注: 表中所有結(jié)果均由T_TIDE程序運(yùn)行得到

    通過(guò)相關(guān)性分析ZWS異常與RWL異常存在顯著的負(fù)相關(guān), 相關(guān)系數(shù)為–0.54,<0.05; RWL與MWS和風(fēng)應(yīng)力旋度間相關(guān)性較差。通過(guò)對(duì)ZWS的Morlet小波變換分析結(jié)果[圖4(c)]可知其存在約5個(gè)月左右的顯著變化周期, 與RWL的顯著周期對(duì)應(yīng)。以上分析說(shuō)明RWL對(duì)緯向風(fēng)的影響較為敏感, 本文通過(guò)帶通濾波分別獲取RWL和ZWS異常的110~150天季節(jié)內(nèi)變化信號(hào), 結(jié)果如圖4(d)所示, 二者相關(guān)系數(shù)為–0.6,<0.05。綜上所述, RWL季節(jié)內(nèi)變化主要受緯向風(fēng)的影響, 監(jiān)測(cè)系統(tǒng)處在大沽河入??谖靼? 西(東)風(fēng)將驅(qū)動(dòng)水體向河口東(西)岸輸運(yùn), 西岸出水量減少(增加), 致使監(jiān)測(cè)系統(tǒng)觀測(cè)到RWL下降(上升)。

    圖4 天平均RWL時(shí)間序列(a), 天平均RWL小波系數(shù)實(shí)部(b), ZWS小波系數(shù)實(shí)部(c)和110~150天RWL和ZWS的帶通濾波信號(hào)(d)

    本節(jié)中還計(jì)算了RWL的季節(jié)內(nèi)變化同流域內(nèi)P點(diǎn)(36.75°N, 121.25°E, 圖1)降水量異常的相關(guān)性, 二者間相關(guān)系數(shù)為–0.13, 表明降水量的變化并非影響RWL季節(jié)內(nèi)變化的主要原因。

    2.3 線性趨勢(shì)

    圖5(a)展示了將觀測(cè)期間內(nèi)天平均RWL進(jìn)行逐月平均處理后得到的月平均RWL時(shí)間序列, 其呈現(xiàn)出顯著的下降趨勢(shì), 約為–0.53×10–2m/月, 通過(guò)顯著性水平為0.05的Mann-Kendall趨勢(shì)檢驗(yàn)[34]。前人研究指出風(fēng)場(chǎng)、外海相對(duì)海平面、流域徑流等因素影響河口水位的變化[31, 35-38], 高宗軍等[11]曾利用歷史資料分析指出大沽河降水與徑流間存在較好的正相關(guān)關(guān)系, 故在缺少?gòu)搅鲾?shù)據(jù)的情況下, 本節(jié)將研究風(fēng)場(chǎng)、外海水位、大范圍降水等因素對(duì)RWL趨勢(shì)的影響。

    圖5 觀測(cè)期間, 月平均RWL及其趨勢(shì)(a), 月平均ZWS及其趨勢(shì)(b), P站降水量異常時(shí)間序列(c)

    圖5(b)展示了ZWS異常的時(shí)間序列及其趨勢(shì), 其下降趨勢(shì)約為–11.29 Ns/(m2月), 對(duì)應(yīng)東風(fēng)在觀測(cè)期間有增強(qiáng)的趨勢(shì), 在3.2節(jié)中分析可知東風(fēng)對(duì)RWL的貢獻(xiàn)主要是使得RWL升高, 與觀測(cè)RWL變化趨勢(shì)不符, 所以在觀測(cè)期間內(nèi), 緯向風(fēng)的變化趨勢(shì)(東風(fēng)增強(qiáng))并非是RWL呈下降趨勢(shì)的主要原因。圖6(a)展示了觀測(cè)期間外海SLA的趨勢(shì)分布, 結(jié)果顯示膠州灣外海水位在觀測(cè)期間呈現(xiàn)上升趨勢(shì), 外海水位上升將使得更多外海海水進(jìn)入膠州灣, 致使灣內(nèi)包括大沽河口區(qū)水位整體上升, 這與觀測(cè)結(jié)果不符。通過(guò)討論局地風(fēng)場(chǎng)(緯向風(fēng))和外海SLA的變化對(duì)RWL的影響, 二者并非觀測(cè)期間大沽河口RWL呈下降趨勢(shì)的主要影響因素, 故作者推測(cè)致使RWL呈下降趨勢(shì)的影響因素更多來(lái)源于大沽河口上游。于是作者計(jì)算了觀測(cè)期間中國(guó)東部大范圍降水趨勢(shì)分布, 如圖6(b)所示, 觀測(cè)期間山東半島降水量普遍呈現(xiàn)下降趨勢(shì), 圖5(c)以大沽河周邊的P測(cè)站降水量異常時(shí)間序列為例展示了觀測(cè)期間的降水變化及其趨勢(shì), 同前人[40]研究結(jié)果一致, 降水極值集中于汛期(7~8月), 其趨勢(shì)約為–0.037 mm/月, 相似的變化的趨勢(shì)在大沽河流域周邊站點(diǎn)同樣存在。流域降水減少使得進(jìn)入大沽河補(bǔ)充徑流的水體減少, 另外降水減少還將使得流域內(nèi)人類生產(chǎn)生活對(duì)大沽河河水汲取量增加, 以上兩個(gè)過(guò)程導(dǎo)致大沽河流域徑流的減少, 致使大沽河口監(jiān)測(cè)系統(tǒng)觀測(cè)到RWL的下降。

    3 結(jié)論與展望

    本文通過(guò)分析大沽河口監(jiān)測(cè)系統(tǒng)2016年11月8日—2019年6月29日的水深觀測(cè)數(shù)據(jù), 對(duì)水位變化的潮汐過(guò)程、季節(jié)內(nèi)變化、觀測(cè)期間變化趨勢(shì)等特征, 以及影響因素展開了分析討論, 得到了以下結(jié)論:

    1) 大沽河口水位變化由潮汐過(guò)程主導(dǎo), 每月都存在兩個(gè)大潮過(guò)程和兩個(gè)小潮過(guò)程; 潮汐類型為正規(guī)半日潮, 其中M2分潮的振幅最大, 達(dá)1.23 m, 淺水分潮(K1、O1)和全日分潮(M4、MS4)的振幅相當(dāng), 二者振幅在0.15~0.23 m范圍內(nèi)。

    2) 緯向風(fēng)變化與日平均RWL呈顯著相關(guān), 具體表現(xiàn)為ZWS異常與RWL異常間存在顯著負(fù)相關(guān), 相關(guān)系數(shù)為–0.54,<0.05; 通過(guò)Morlet小波變化分析方法得知RWL異常在2017年7月—2019年1月存在110~150天左右的顯著季節(jié)內(nèi)變化, 用同樣的方法可以發(fā)現(xiàn)ZWS在2017年10月—2018年10月期間存在約5個(gè)月左右的顯著變化周期, 與RWL的變化對(duì)應(yīng), 進(jìn)一步通過(guò)帶通濾波分別獲取RWL和ZWS異常110~150 d信號(hào), 二者相關(guān)系數(shù)為–0.60,<0.05。RWL的季節(jié)內(nèi)變化主要受緯向風(fēng)的影響, 西(東)風(fēng)將驅(qū)動(dòng)水體向河口東(西)岸輸運(yùn), 西岸出水量減少(增加), 致使監(jiān)測(cè)系統(tǒng)觀測(cè)到RWL下降(上升)。

    圖 6 觀測(cè)期間, 海平面異常趨勢(shì)分布(a), 地面降水量趨勢(shì)分布(b)

    注: 子圖a、b中紅色五角星為監(jiān)測(cè)系統(tǒng)位置, 子圖b中品紅色虛線框內(nèi)區(qū)域代表大沽河流域

    3) 在觀測(cè)期間, 月平均RWL還呈現(xiàn)顯著下降趨勢(shì), 約為–0.53×10–2m/月, 主要受到大沽河流域內(nèi)降水減少的影響, 本文推測(cè)降水量的減少對(duì)RWL下降的影響可以分為兩方面: 其一是降水減少直接導(dǎo)致流域內(nèi)徑流減少, 另外降水減少還將使得流域內(nèi)人類生產(chǎn)生活對(duì)大沽河河水汲取量增加。

    本文通過(guò)對(duì)大沽河口長(zhǎng)期、連續(xù)水深觀測(cè)的研究, 揭示了河口的潮汐過(guò)程, 發(fā)現(xiàn)了觀測(cè)期間水位呈現(xiàn)的季節(jié)內(nèi)尺度變化以及水位的變化趨勢(shì), 針對(duì)影響水位變化的動(dòng)力因素得到了一些初步結(jié)論, 受到觀測(cè)資料時(shí)間長(zhǎng)度和分析數(shù)據(jù)種類的限制, 需進(jìn)一步結(jié)合觀測(cè)資料分析探究。

    [1] Alley R B, Clark P U, Huybrechts P, et al. Ice-sheet and sea-level changes[J]. Science, 2005, 310(5747): 456-460.

    [2] 王芳. 海平面上升的影響及損失預(yù)測(cè)[J]. 上海環(huán)境科學(xué), 1998, 17(10): 9-11. Wang Fang. Impact of sea level rise and loss forecast[J]. Shanghai Environmental Sciences, 1998, 17(10): 9-11.

    [3] 韓小燕, 潘曉東, 馬林芳, 等. 溫州沿岸平均海平面變化特征分析[J]. 海洋預(yù)報(bào), 2011, 28(1): 66-71.Han Xiaoyan, Pan Xiaodong, Ma Linfang, et al. Analysis on the characteristics of the average sea level change along the coast of Wenzhou[J]. Marine Forecasts, 2011, 28(1): 66-71.

    [4] 于華. 氣候變暖背景下風(fēng)場(chǎng)對(duì)江蘇近海海平面的影響研究[D]. 江蘇南京: 南京信息工程大學(xué), 2014. Yu Hua. Research on the influence of wind field on the sea level of Jiangsu coast under the background of climate warming[D]. Nanjing, Jiangsu- Province: Nanjing University of Information Science & Technology, 2014.

    [5] 宗虎城, 章衛(wèi)勝, 張金善. 海平面上升對(duì)黃浦江風(fēng)暴潮水位影響研究[J]. 人民長(zhǎng)江, 2014, 45(9): 1-3.Zong Hucheng, Zhang Weisheng, Zhang Jinshan. Resea-rch on the impact of sea level rising on the storm surge level of the Huangpu River[J]. Yangtze River, 2014, 45(9): 1-3.

    [6] 徐曉曉. 海岸帶地區(qū)風(fēng)暴潮災(zāi)害風(fēng)險(xiǎn)研究[D]. 山東青島: 中國(guó)海洋大學(xué), 2015. Xu Xiaoxiao. Research on the disaster risk of storm surge in coastal zone[D]. Qingdao, Shandong Province: Ocean University of China, 2015.

    [7] 王舒鴻, 郭越, 李逸超, 等. 城市海岸帶環(huán)境管理模式優(yōu)化研究——以青島市膠州灣為例[J]. 中國(guó)環(huán)境管理, 2019, 11(2): 69-75, 79.Wang Shuhong, Guo Yue, Li Yichao, et al. Research on optimization of urban coastal zone environmental management mode—Taking Jiaozhou Bay in Qingdao as an example[J]. Chinese Journal of Environmental Management, 2019, 11(2): 69-75, 79.

    [8] 周天華, 陳宗鏞, 田暉, 等. 近幾十年來(lái)中國(guó)沿岸海面變化趨勢(shì)的研究[J]. 海洋學(xué)報(bào)(中文版), 1992, 14(2): 1-8. Zhou Tianhua, Chen Zongyong, Tian Hui, et al. Research on the changing trend of China’s coastal sea surface in recent decades[J]. Haiyang Xuebao, 1992, 14(2): 1-8.

    [9] 周中剛, 蔣國(guó)榮, 沈春, 等. 中國(guó)近海平均海平面計(jì)算方法的分析與比較[J]. 海洋預(yù)報(bào), 2013, 30(5): 37-43. Zhou Zhonggang, Jiang Guorong, Shen Chun, et al. Analysis and comparison of the calculation methods of the Coast Areas in China mean sea level[J]. Marine Forecasts, 2013, 30(5): 37-43.

    [10] 劉聚, 暴景陽(yáng), 許軍, 等. 中國(guó)香港驗(yàn)潮站1962-2017年水位相對(duì)變化分析[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2020, 45(7): 1065-1072.Liu Ju, Bao Jingyang, Xu Jun, et al. Analysis of the relative change of water level at tide gauge station in Hong Kong, China from 1962 to 2017[J]. Geomatics and Information Science of Wuhan University, 2020, 45(7): 1065-1072.

    [11] 高宗軍, 宋翠玉, 蔡玉林, 等. 大沽河流域水文要素監(jiān)測(cè)體系建設(shè)與實(shí)踐[M]. 北京: 水利水電出版社. 2017: 1-133.Gao Zongjun, Song Cuiyu, Cai Yulin, et al. Construction and practice of hydrological elements monitoring system in Dagu River basin[M]. Beijing: Water Resources and Hydropower Press. 2017: 1-133.

    [12] 盛茂剛, 黃修東, 左林遠(yuǎn), 等. 大沽河流域近60年降水量時(shí)空變化特征分析[J]. 水資源與水工程學(xué)報(bào), 2016, 27(2): 65-68. Sheng Maogang, Huang Xiudong, Zuo Linyuan, et al. Analysis of temporal and spatial variation characteristics of precipitation in Dagu River basin for recent 60 years[J]. Journal of Water Resources & Water Engineering, 2016, 27(2): 65-68.

    [13] 王艷玲, 崔文連, 劉峰, 等. 青島市大沽河河口區(qū)生態(tài)環(huán)境現(xiàn)狀研究[J]. 中國(guó)環(huán)境監(jiān)測(cè), 2007, 23(3): 77-81. Wang Yanling, Cui Wenlian, Liu Feng, et al. Study on the ecological environment status of Dagu River Estuary in Qingdao[J]. Environmental Monitoring in China, 2007, 23(3): 77-81.

    [14] 柳新偉, 秦志華. 大沽河河口灘涂植被群落動(dòng)態(tài)特征研究[J]. 青島農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010, 27(3): 241-243. Liu Xinwei, Qin Zhihua. Study on dynamic characteris-tics of tidal flat vegetation communities in Dagu River Estuary[J]. Journal of Qingdao Agricultural Universi-ty(Natural Science), 2010, 27(3): 241-243.

    [15] 張曉紅, 汪進(jìn)生, 劉旭東. 2014年秋季大沽河河口浮游植物多樣性研究[J]. 環(huán)境科學(xué)與管理, 2015, 40(7): 129-131. Zhang Xiaohong, Wang Jinsheng, Liu Xudong. Study on phytoplankton diversity of Dagu River estuary in autumn 2014[J]. Environmental Science and Management, 2015, 40(7): 129-131.

    [16] 劉娜, 印萍, 朱志剛, 等. 膠州灣大沽河河口表層沉積物中多環(huán)芳烴分布特征、來(lái)源及生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J]. 2016, 35(6): 831-837. Liu Na, Yin Ping, Zhu Zhigang, et al. Distribution characteristics, sources and ecological risk assessment of polycyclic aromatic hydrocarbons in surface sediments of Dagu River estuary, Jiaozhou Bay[J]. 2016, 35(6): 831-837.

    [17] 婁安剛, 王學(xué)昌, 吳德星, 等. 膠州灣大沽河口鄰近海域海水水質(zhì)預(yù)測(cè)[J]. 海洋環(huán)境科學(xué), 2002, 21(1): 54-56. Lou Angang, Wang Xuechang, Wu Dexing, et al. Predi-ction of sea water quality in the adjacent sea area of Dagu River estuary in Jiaozhou Bay[J]. Marine Environmental Science, 2002, 21(1): 54-56.

    [18] 韓彬, 曹磊, 李培昌, 等. 膠州灣大沽河河口及鄰近海域海水水質(zhì)狀況與評(píng)價(jià)[J]. 海洋科學(xué), 2010, 34(8): 46-49.Han Bin, Cao Lei, Li Peichang, et al. Status and evalua-tion of sea water quality at Dagu River estuary in Jiao-zhou Bay and adjacent seas[J]. Marine Sciences, 2010, 34(8): 46-49.

    [19] 韓樹宗, 趙瑾, 魏福寶, 等. 膠州灣大沽河口洪水期三維水沙數(shù)值模擬研究[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2007, 37(5): 689-694.Han Shuzong, Zhao Jin, Wei Fubao, et al. Three-dimen-sional numerical simulation of water and sediment in Dagu estuary of Jiaozhou Bay in flood period[J]. Periodi-cal of Ocean University of China, 2007, 37(5): 689-694.

    [20] 徐統(tǒng). 徑流、潮流耦合作用下大沽河河口懸沙輸移擴(kuò)散數(shù)值研究[D]. 山東青島: 中國(guó)海洋大學(xué), 2014.Xu Tong. Numerical study on suspended sediment transport and diffusion in the Dagu River estuary under the coupling action of runoff and tidal current[D]. Qingdao, Shandong Province: Ocean University of China, 2014.

    [21] 郭迢. 海灣河口泥質(zhì)潮灘海水—地下水相互交換定量化研究—以膠州灣大沽河河口某剖面為例[D]. 北京: 中國(guó)地質(zhì)大學(xué), 2017.Guo Tiao. Quantitative study on the exchange of seawater and groundwater in the muddy tidal flat of the gulf estuary-taking a section at the mouth of Dagu River in Jiaozhou Bay as an example[D]. Beijing: Geosicence Uni-versity of China, 2017.

    [22] Xiong G, An Q, Fu T, et al. Evolution analysis and environmental management of intruded aquifers of the Dagu River Basin of China[J]. Science of The Total Environment, 2020, 719: 137260.

    [23] 尹則高, 徐統(tǒng), 王振魯, 等. 潮流和徑流作用下大沽河河口區(qū)水動(dòng)力特性研究[J]. 中國(guó)海洋大學(xué)學(xué)報(bào), 2015, 45(7): 119-124. Yin Zegao, Xu Tong, Wang Zhenlu, et al. Dagu estuary hydrodynamic characteristic research under the coupled action of runoff and tide[J]. Periodical of Ocean University of China, 2015, 45(7): 119-124.

    [24] 李兆欽, 李欣, 孫利元, 等. 大沽河口底層海水溶解氧濃度分析[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2019, 49(9): 23-33. Li Zhaoqin, Li Xin, Sun Liyuan, et al. Analysis of the dissolved oxygen concentration in bottom of seawater in the estuary of Dagu River[J]. Periodical of Ocean University of China, 2019, 49(9): 23-33.

    [25] Fofonoff N P, Millard R C. Algorithms for computation of fundamental properties of seawater[J]. Unesco Technical Papers in Marine Science, 1983, 44: 1-53.

    [26] Dee D P, Uppala S M, Simmons A J, et al. The ERA- Interim reanalysis: configuration and performance of the data assimilation system[J]. Quarterly Journal of the Royal Meteorological Society, 2011, 137(656): 553-597.

    [27] Chen M Y, Xie P P, Janowiak J E, et al. Global land precipitation: a 50-yr monthly analysis based on gauge observations[J]. Journal of Hydrometeorology, 2002, 3(3): 249-266.

    [28] Fan Y, Dool H. A global monthly land surface air temperature analysis for 1948—present[J]. Journal of Geophysical Research, 2008, 113(D1): JD008470.1?JD008470.7.

    [29] 吳嫻, 黃偉, 陳發(fā)虎. 1951—2012年中國(guó)大陸0.025°× 0.025°高分辨率月氣溫和降水量格點(diǎn)數(shù)據(jù)集的建立及其初步應(yīng)用[J]. 蘭州大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014, 50(2): 213-220.Wu Xian, Huang Wei, Chen Fahu. Construction and application of monthly air temperature and precipitation gridded datasets with high resolution (0.025°×0.025°) over China during 1951—2012[J]. Journal of Lanzhou University (Natural Sciences), 2014, 50(2): 213-220.

    [30] 姚武川, 姚天任. 經(jīng)典譜估計(jì)方法的MATLAB分析[J]. 華中理工大學(xué)學(xué)報(bào), 2000, 28(4): 45-47. Yao Wuchuan, Yao Tianren. Analyzing classical spectral estimation by MATLAB[J]. Journal of Huazhong University of Science and Technology(Natural Science Edition), 2000, 28(4): 45-47.

    [31] 裴文斌, 牛桂芝, 董海軍. 余水位及潮汐差分方法[J]. 水道港口, 2007, 28(6): 439-443. Pei Wenbin, Niu Guizhi, Dong Haijun. Residual water level and tide difference method[J]. Journal of Waterway and Harbor, 2007, 28(6): 439-443.

    [32] Torrence C, Compo G P. A practical guide to wavelet analysis[J]. Bulletin of the American Meteorological Society, 1998, 79(1): 61-78.

    [33] 鮑道陽(yáng), 朱建榮. 近60年來(lái)長(zhǎng)江河口河勢(shì)變化及其對(duì)水動(dòng)力和鹽水入侵的影響Ⅱ. 水動(dòng)力[J]. 海洋學(xué)報(bào), 2017, 39(2): 1-15. Bao Daoyang, Zhu Jianrong. The effects of river regime changes in the Changjiang Estuary on hydrodynamics and salinity intrusion in the past 60 years Ⅱ. Hydrodynamics[J]. Acta Oceanologica Sinica, 2017, 39(2): 1-15.

    [34] Sneyers R. On the statistical analysis of series of observations, WMO Technical Note No. 143[R]. Geneva, Switzerland: Secretariat of the World Meteorological Organization, 1990.

    [35] Piecuch C G, Bittermann K, Kemp A C, et al. River- discharge effects on United States Atlantic and Gulf coast sea-level changes[J]. Proceedings of the National Academy of Sciences, 2018, 115(30): 7729-7734.

    [36] Jay D A, Leffler K, Degens S. Long-term evolution of Columbia River tides[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 2011, 137(4): 182-191.

    [37] Ralston D K, Talke S, Geyer W R, et al. Bigger tides, less flooding: Effects of dredging on water level in the Hudson River estuary[J]. Journal of Geophysical Research: Oceans, 2019, 124(1): 196-211.

    [38] 楊昊, 歐素英, 姚鵬, 等. 河口區(qū)斜壓梯度對(duì)余水位的累積影響及其機(jī)制探討[J]. 海洋學(xué)報(bào), 2019, 41(1): 21-31. Yang Hao, Ou Suying, Yao Peng, et al. The cumulative effect of baroclinic pressure gradient on residual water level in estuary area and its mechanism[J]. Haiyang Xuebao, 2019, 41(1): 21-31.

    [39] 姜德娟, 王曉利. 膠東半島大沽河流域徑流變化特征[J]. 干旱區(qū)研究, 2013, 30(6): 965-972. Jiang Dejuan, Wang Xiaoli. Variation of runoff volume in the Dagu River Basin in the Jiaodong peninsula[J]. Arid Zone Research, 2013, 30(6): 965-972.

    Analysis of the characteristics of water level variation in Dagu River estuary

    LAI Qiu-yu1, LIU Zi-zhou1, ZHAI Fang-guo1, GU Yan-zhen1, CHEN Dong2, JIANG Qing-yan2

    (1. College of Ocean and Atmospheric Sciences, Ocean University of China, Qingdao 266100, China; 2. Ocean College, Zhejiang University, Zhoushan 316021, China)

    This study analyzed the characteristics of water level variation in the Dagu River estuary and its influencing factors combined with the wind field, precipitation, and satellite altimetry data based on the depth observation data of Dagu River estuary from November 8, 2016 to June 29, 2019. Results showed that (1) the water level variation in Dagu River estuary was dominated by a normal semidiurnal tidal process and the M2 subtide was dominant. (2) The residual water level (RWL) had a significant intraseasonal period of approximately 110–150 days, which was mainly affected by zonal wind, i.e., the eastward (westward) wind would drive water transport to the east (west), decreasing (increasing) the water volume at the west bank at which the monitoring system was located, leading to the decrease (increase) in RWL. (3) During the observation period, the RWL exhibited a significant downward trend, approximately ?0.53×10?2m/month, which was mainly affected by the decrease in precipitation in the Dagu River basin.

    Dagu River estuary; water level variation; zonal wind; precipitation

    Aug. 21, 2020

    P731.2

    A

    1000-3096(2021)01-0022-10

    10.11759/hykx20200821002

    2020-08-21;

    2020-09-26

    國(guó)家自然科學(xué)基金項(xiàng)目(41776012); 國(guó)家重點(diǎn)研發(fā)計(jì)劃 (2019YFD0901305); 浙江省重點(diǎn)研發(fā)計(jì)劃(2020C03012); 三亞崖州灣科技城管理局重大科技項(xiàng)目(YZ2019ZD0X)

    [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. YZ2019ZD0X]

    賴秋宇(1995-), 男, 廣西桂林市人, 碩士研究生, 主要從事海洋水動(dòng)力研究, E-mail: QYordinary@163.com; 劉子洲,通信作者, 男, 山東省淄博人, 實(shí)驗(yàn)師, 主要從事海洋調(diào)查及近海海洋環(huán)境問(wèn)題的教學(xué)和研究, E-mail: lzz2013@ouc.edu.cn

    (本文編輯: 趙衛(wèi)紅)

    猜你喜歡
    大沽河口水深
    書法靜水深流
    河北水利(2022年10期)2022-12-29 11:48:12
    基于水深分段選擇因子的多光譜影像反演水深
    大沽河
    大沽河
    北極光(2018年5期)2018-11-19 02:23:56
    他們?yōu)槭裁催x擇河口
    河口,我們的家
    特殊的河口水
    河口
    GPS RTK技術(shù)在水深測(cè)量中的應(yīng)用
    浸入式水深監(jiān)測(cè)儀器的設(shè)計(jì)
    亚洲av欧美aⅴ国产| 大陆偷拍与自拍| 欧美激情国产日韩精品一区| 国产乱人偷精品视频| 国产淫片久久久久久久久| 一区二区三区精品91| 午夜免费观看性视频| 国产 一区 欧美 日韩| 在线亚洲精品国产二区图片欧美 | 欧美日韩视频精品一区| 日韩中文字幕视频在线看片 | 丝瓜视频免费看黄片| 高清不卡的av网站| 最近最新中文字幕免费大全7| 亚洲av国产av综合av卡| 日韩人妻高清精品专区| 日本黄大片高清| av国产免费在线观看| 中国三级夫妇交换| 亚洲国产欧美人成| 国产精品国产三级专区第一集| 色综合色国产| 99久久精品一区二区三区| 中文字幕制服av| 日日啪夜夜爽| 久久久久国产网址| 免费观看在线日韩| 夫妻性生交免费视频一级片| 亚洲av在线观看美女高潮| 简卡轻食公司| 久久精品久久精品一区二区三区| 久久综合国产亚洲精品| 黄色怎么调成土黄色| 一二三四中文在线观看免费高清| 国产男女内射视频| 国产精品国产三级专区第一集| 久久av网站| 91精品一卡2卡3卡4卡| 精品人妻一区二区三区麻豆| 久久精品国产亚洲av涩爱| 国内少妇人妻偷人精品xxx网站| 在线观看免费日韩欧美大片 | 一区二区三区免费毛片| 国产免费福利视频在线观看| 国产日韩欧美亚洲二区| 中文字幕人妻熟人妻熟丝袜美| 国产成人aa在线观看| 免费看av在线观看网站| 看免费成人av毛片| 国产精品一及| 日本爱情动作片www.在线观看| 久久人人爽av亚洲精品天堂 | 亚洲婷婷狠狠爱综合网| 日韩中文字幕视频在线看片 | 一本一本综合久久| 亚洲成人手机| 黄色视频在线播放观看不卡| 国产在线一区二区三区精| 亚洲av中文字字幕乱码综合| 欧美高清性xxxxhd video| 国产在视频线精品| 校园人妻丝袜中文字幕| 男男h啪啪无遮挡| 亚洲欧美日韩卡通动漫| 中文字幕精品免费在线观看视频 | 夜夜看夜夜爽夜夜摸| 亚洲电影在线观看av| 在线观看免费视频网站a站| 国产深夜福利视频在线观看| 国产精品久久久久久精品古装| 狂野欧美激情性bbbbbb| 秋霞在线观看毛片| 国产午夜精品一二区理论片| 国产在线免费精品| 大话2 男鬼变身卡| 天堂中文最新版在线下载| 日韩,欧美,国产一区二区三区| 日韩一区二区视频免费看| 内射极品少妇av片p| 人妻制服诱惑在线中文字幕| 黄色日韩在线| 亚洲av欧美aⅴ国产| 一区在线观看完整版| 国产精品国产三级专区第一集| 国产高清不卡午夜福利| 日日摸夜夜添夜夜添av毛片| 看免费成人av毛片| 夜夜看夜夜爽夜夜摸| 深夜a级毛片| 国产欧美另类精品又又久久亚洲欧美| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品亚洲一区二区| 在线看a的网站| 一个人看的www免费观看视频| 亚洲经典国产精华液单| 国产精品一区二区在线观看99| 国产成人精品福利久久| 午夜福利在线观看免费完整高清在| 成人国产麻豆网| 日韩欧美精品免费久久| 老熟女久久久| 熟女av电影| 国产白丝娇喘喷水9色精品| 国产在线免费精品| 精品少妇黑人巨大在线播放| 欧美高清性xxxxhd video| 亚洲精品中文字幕在线视频 | 多毛熟女@视频| 美女cb高潮喷水在线观看| 如何舔出高潮| 中文字幕免费在线视频6| 国产无遮挡羞羞视频在线观看| 国产69精品久久久久777片| 在线免费观看不下载黄p国产| 建设人人有责人人尽责人人享有的 | 人妻少妇偷人精品九色| 国产av一区二区精品久久 | 国产极品天堂在线| 成人亚洲欧美一区二区av| 国产伦理片在线播放av一区| 人人妻人人看人人澡| 欧美激情极品国产一区二区三区 | 我的老师免费观看完整版| 人体艺术视频欧美日本| 欧美人与善性xxx| 亚洲精品国产av蜜桃| 久久久久久人妻| 成人一区二区视频在线观看| 亚洲色图av天堂| 国产黄频视频在线观看| 欧美少妇被猛烈插入视频| 一级毛片久久久久久久久女| 国产无遮挡羞羞视频在线观看| 国产精品三级大全| 亚洲熟女精品中文字幕| 精品人妻视频免费看| 九九久久精品国产亚洲av麻豆| 国产视频首页在线观看| 下体分泌物呈黄色| 老熟女久久久| 久久鲁丝午夜福利片| 国产精品无大码| 国产精品欧美亚洲77777| 激情五月婷婷亚洲| 亚洲国产欧美人成| 亚洲精品国产色婷婷电影| 王馨瑶露胸无遮挡在线观看| 久久99热6这里只有精品| 免费少妇av软件| 亚洲欧美日韩东京热| 天天躁夜夜躁狠狠久久av| 18+在线观看网站| 国产精品嫩草影院av在线观看| 久久久久性生活片| videos熟女内射| 日本黄色日本黄色录像| 亚洲一级一片aⅴ在线观看| 国产人妻一区二区三区在| 尾随美女入室| av又黄又爽大尺度在线免费看| 亚洲欧美一区二区三区国产| 99精国产麻豆久久婷婷| 亚洲欧洲国产日韩| 欧美区成人在线视频| 国产伦精品一区二区三区四那| 麻豆成人av视频| 人体艺术视频欧美日本| 亚洲电影在线观看av| 国产午夜精品久久久久久一区二区三区| 午夜日本视频在线| 日韩大片免费观看网站| 美女xxoo啪啪120秒动态图| av不卡在线播放| 91精品国产九色| 赤兔流量卡办理| 欧美日韩精品成人综合77777| 日韩免费高清中文字幕av| 久久这里有精品视频免费| 观看av在线不卡| 一区二区三区乱码不卡18| 小蜜桃在线观看免费完整版高清| 国产成人91sexporn| 欧美xxxx性猛交bbbb| h视频一区二区三区| 国产白丝娇喘喷水9色精品| 一级毛片黄色毛片免费观看视频| 99国产精品免费福利视频| 久久久久久久久久人人人人人人| 国产免费福利视频在线观看| 蜜臀久久99精品久久宅男| 18禁裸乳无遮挡动漫免费视频| 久久国产亚洲av麻豆专区| 国产亚洲一区二区精品| 久久婷婷青草| 看免费成人av毛片| 一个人看的www免费观看视频| 欧美最新免费一区二区三区| 久久精品人妻少妇| 97超视频在线观看视频| 嫩草影院新地址| 嘟嘟电影网在线观看| 久久婷婷青草| av天堂中文字幕网| 久久女婷五月综合色啪小说| 亚洲av男天堂| 少妇裸体淫交视频免费看高清| 成人漫画全彩无遮挡| 国产精品久久久久久av不卡| 最近手机中文字幕大全| 水蜜桃什么品种好| 老师上课跳d突然被开到最大视频| 十八禁网站网址无遮挡 | 男女下面进入的视频免费午夜| 国产老妇伦熟女老妇高清| 在线精品无人区一区二区三 | 日韩成人伦理影院| 九色成人免费人妻av| 狠狠精品人妻久久久久久综合| 大话2 男鬼变身卡| 成人无遮挡网站| 小蜜桃在线观看免费完整版高清| xxx大片免费视频| 超碰97精品在线观看| 亚洲第一区二区三区不卡| 日韩成人av中文字幕在线观看| 亚洲第一av免费看| 99热国产这里只有精品6| 99久久综合免费| 99热全是精品| 成人一区二区视频在线观看| 国产伦在线观看视频一区| 国产黄色免费在线视频| 狂野欧美激情性xxxx在线观看| 免费黄色在线免费观看| 丝袜喷水一区| 国产黄片美女视频| 大话2 男鬼变身卡| 十分钟在线观看高清视频www | 国产精品国产三级专区第一集| 国产精品免费大片| 乱系列少妇在线播放| 男人添女人高潮全过程视频| 久久鲁丝午夜福利片| 国产一区二区三区综合在线观看 | 十八禁网站网址无遮挡 | 最黄视频免费看| 国产精品99久久99久久久不卡 | 久久ye,这里只有精品| 青春草亚洲视频在线观看| 亚洲精品久久午夜乱码| 激情五月婷婷亚洲| av女优亚洲男人天堂| 99热这里只有是精品50| 国产精品一区二区在线观看99| 国产亚洲欧美精品永久| 尤物成人国产欧美一区二区三区| 最新中文字幕久久久久| 香蕉精品网在线| 国产成人精品久久久久久| 欧美性感艳星| 亚洲精品中文字幕在线视频 | 丰满迷人的少妇在线观看| 国语对白做爰xxxⅹ性视频网站| 日韩 亚洲 欧美在线| 免费少妇av软件| 久久国产精品大桥未久av | 一本—道久久a久久精品蜜桃钙片| 午夜老司机福利剧场| 亚洲欧美中文字幕日韩二区| 日本黄色片子视频| 国产免费又黄又爽又色| 日日摸夜夜添夜夜添av毛片| 少妇被粗大猛烈的视频| 国产在视频线精品| 老女人水多毛片| 午夜福利网站1000一区二区三区| 亚洲精品色激情综合| 亚洲va在线va天堂va国产| 亚洲,欧美,日韩| 日本猛色少妇xxxxx猛交久久| 18禁动态无遮挡网站| 黄色配什么色好看| 亚洲欧美精品专区久久| 少妇被粗大猛烈的视频| 亚洲自偷自拍三级| 成人18禁高潮啪啪吃奶动态图 | 婷婷色av中文字幕| 久久久久久久久久成人| 欧美另类一区| a级一级毛片免费在线观看| 最近最新中文字幕免费大全7| 亚洲国产成人一精品久久久| 观看av在线不卡| 麻豆精品久久久久久蜜桃| 精品酒店卫生间| 精品久久久久久电影网| av在线播放精品| 最近手机中文字幕大全| 国产69精品久久久久777片| 在线观看人妻少妇| 网址你懂的国产日韩在线| 又大又黄又爽视频免费| 久久韩国三级中文字幕| 丝袜喷水一区| 中国三级夫妇交换| 最近最新中文字幕大全电影3| 欧美日韩精品成人综合77777| 嫩草影院入口| 成人一区二区视频在线观看| 制服丝袜香蕉在线| 久久久久久久久大av| 伦理电影免费视频| 欧美xxxx性猛交bbbb| 午夜免费男女啪啪视频观看| 91aial.com中文字幕在线观看| 卡戴珊不雅视频在线播放| 日韩视频在线欧美| 国产真实伦视频高清在线观看| 亚洲精品久久午夜乱码| 国产精品蜜桃在线观看| 欧美老熟妇乱子伦牲交| 男女免费视频国产| 国产一区二区三区综合在线观看 | 国产熟女欧美一区二区| 亚洲av福利一区| 免费黄色在线免费观看| 你懂的网址亚洲精品在线观看| 欧美日韩精品成人综合77777| 在线天堂最新版资源| 少妇人妻一区二区三区视频| 99九九线精品视频在线观看视频| 水蜜桃什么品种好| av专区在线播放| 午夜日本视频在线| 国产亚洲最大av| 美女cb高潮喷水在线观看| 夜夜爽夜夜爽视频| 中文天堂在线官网| 男女下面进入的视频免费午夜| 91精品国产国语对白视频| 国产高清三级在线| 国产精品国产三级专区第一集| 久久99热这里只有精品18| 国产欧美另类精品又又久久亚洲欧美| 亚洲怡红院男人天堂| 成人18禁高潮啪啪吃奶动态图 | 亚洲欧美日韩另类电影网站 | 久久久久人妻精品一区果冻| 亚洲欧美成人精品一区二区| 久久久久人妻精品一区果冻| 国产伦在线观看视频一区| 亚洲av福利一区| 免费看日本二区| 日本免费在线观看一区| av一本久久久久| 国产精品99久久99久久久不卡 | 亚洲三级黄色毛片| 精品熟女少妇av免费看| 啦啦啦在线观看免费高清www| 高清毛片免费看| 国内揄拍国产精品人妻在线| 蜜桃久久精品国产亚洲av| 欧美bdsm另类| 欧美高清性xxxxhd video| 国产精品无大码| 亚洲成人一二三区av| 黄片无遮挡物在线观看| 18禁动态无遮挡网站| 最黄视频免费看| 中文字幕制服av| 欧美bdsm另类| 精品亚洲乱码少妇综合久久| 99久久中文字幕三级久久日本| 亚洲人成网站在线播| 久久精品人妻少妇| 91久久精品国产一区二区成人| 少妇精品久久久久久久| 欧美xxⅹ黑人| 日韩强制内射视频| 中文乱码字字幕精品一区二区三区| 极品少妇高潮喷水抽搐| 亚洲精品,欧美精品| 亚洲国产精品一区三区| 蜜桃久久精品国产亚洲av| 欧美精品一区二区大全| 伦理电影免费视频| 女的被弄到高潮叫床怎么办| 少妇的逼水好多| 亚洲国产精品999| 美女主播在线视频| 精品一区二区三区视频在线| 成人亚洲精品一区在线观看 | a 毛片基地| 夫妻性生交免费视频一级片| 在线观看国产h片| 成人综合一区亚洲| 边亲边吃奶的免费视频| 久久久久久久精品精品| 少妇的逼水好多| 自拍欧美九色日韩亚洲蝌蚪91 | 国产午夜精品一二区理论片| 亚洲伊人久久精品综合| 国产v大片淫在线免费观看| 我的老师免费观看完整版| 天天躁日日操中文字幕| 久久毛片免费看一区二区三区| 99久久精品热视频| 超碰97精品在线观看| 久久女婷五月综合色啪小说| 国产av国产精品国产| 国产精品国产av在线观看| 亚洲精品乱久久久久久| 日韩欧美精品免费久久| 欧美变态另类bdsm刘玥| 国产极品天堂在线| 五月伊人婷婷丁香| 国产免费一区二区三区四区乱码| a级一级毛片免费在线观看| 精品久久久久久久久亚洲| 国产在视频线精品| 尾随美女入室| 精品久久久久久电影网| 秋霞在线观看毛片| 国产av精品麻豆| 一级片'在线观看视频| 免费黄网站久久成人精品| 亚洲av日韩在线播放| 国产综合精华液| 国产在线一区二区三区精| 精品人妻视频免费看| 水蜜桃什么品种好| 街头女战士在线观看网站| 国产一级毛片在线| 亚洲怡红院男人天堂| 亚洲人成网站高清观看| xxx大片免费视频| 精品亚洲乱码少妇综合久久| 汤姆久久久久久久影院中文字幕| 美女xxoo啪啪120秒动态图| 欧美日韩一区二区视频在线观看视频在线| 亚洲美女搞黄在线观看| 22中文网久久字幕| 亚洲综合精品二区| 91狼人影院| av又黄又爽大尺度在线免费看| 在线亚洲精品国产二区图片欧美 | 啦啦啦视频在线资源免费观看| 亚洲av在线观看美女高潮| 午夜日本视频在线| 国产淫片久久久久久久久| 欧美成人精品欧美一级黄| 免费大片黄手机在线观看| 多毛熟女@视频| 国产深夜福利视频在线观看| 免费观看a级毛片全部| 内地一区二区视频在线| av又黄又爽大尺度在线免费看| 又黄又爽又刺激的免费视频.| 日韩伦理黄色片| 日韩精品有码人妻一区| 亚洲国产色片| 99精国产麻豆久久婷婷| 天美传媒精品一区二区| 久久这里有精品视频免费| 亚洲欧洲国产日韩| 亚洲欧美一区二区三区国产| 九九在线视频观看精品| 欧美一区二区亚洲| 亚洲不卡免费看| 大陆偷拍与自拍| 我要看日韩黄色一级片| 欧美97在线视频| 少妇丰满av| 香蕉精品网在线| 你懂的网址亚洲精品在线观看| 国产精品三级大全| 久久久久久久亚洲中文字幕| 免费人妻精品一区二区三区视频| 99九九线精品视频在线观看视频| 亚洲精品乱码久久久v下载方式| 这个男人来自地球电影免费观看 | 日韩免费高清中文字幕av| 黄色日韩在线| 久久韩国三级中文字幕| 亚洲最大成人中文| 青春草亚洲视频在线观看| 尤物成人国产欧美一区二区三区| 国产黄色免费在线视频| 婷婷色综合www| 亚洲国产日韩一区二区| 欧美另类一区| 亚洲第一av免费看| 久久久久精品久久久久真实原创| 婷婷色av中文字幕| 亚洲精品日韩在线中文字幕| 久久精品夜色国产| 日产精品乱码卡一卡2卡三| 国产爽快片一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 国产有黄有色有爽视频| 99久国产av精品国产电影| 亚洲电影在线观看av| 国产美女午夜福利| 夜夜看夜夜爽夜夜摸| 丰满少妇做爰视频| 波野结衣二区三区在线| 久久 成人 亚洲| 日韩,欧美,国产一区二区三区| 国产亚洲最大av| 亚洲av在线观看美女高潮| 亚洲av成人精品一二三区| 成人二区视频| 午夜日本视频在线| 我要看黄色一级片免费的| 女人久久www免费人成看片| 青青草视频在线视频观看| 国产精品人妻久久久久久| 激情五月婷婷亚洲| 精品视频人人做人人爽| 十分钟在线观看高清视频www | 少妇 在线观看| 在线观看一区二区三区| 亚洲丝袜综合中文字幕| 久久精品久久精品一区二区三区| 午夜福利在线观看免费完整高清在| 美女高潮的动态| 中国三级夫妇交换| 国产在线男女| 伦理电影免费视频| 国产精品不卡视频一区二区| 久久久精品94久久精品| 成人综合一区亚洲| 我的女老师完整版在线观看| 中文字幕免费在线视频6| av在线观看视频网站免费| 国内精品宾馆在线| 国产成人精品福利久久| 日本av手机在线免费观看| 人妻制服诱惑在线中文字幕| 黄色怎么调成土黄色| 亚洲无线观看免费| 亚洲国产精品一区三区| 亚洲无线观看免费| av国产久精品久网站免费入址| 成人毛片a级毛片在线播放| 国产精品99久久99久久久不卡 | 欧美高清性xxxxhd video| 日韩一本色道免费dvd| 国产视频内射| 五月伊人婷婷丁香| 97在线人人人人妻| 2018国产大陆天天弄谢| 日日摸夜夜添夜夜添av毛片| 一级二级三级毛片免费看| 国产成人午夜福利电影在线观看| 午夜免费鲁丝| 亚洲av二区三区四区| 国模一区二区三区四区视频| 国产免费一级a男人的天堂| 亚洲欧美中文字幕日韩二区| 免费人妻精品一区二区三区视频| 午夜精品国产一区二区电影| av免费在线看不卡| 伊人久久国产一区二区| 九色成人免费人妻av| 国产成人一区二区在线| 欧美成人a在线观看| 1000部很黄的大片| 少妇熟女欧美另类| 日韩精品有码人妻一区| 成人美女网站在线观看视频| 亚洲精品国产成人久久av| 免费大片黄手机在线观看| 日本一二三区视频观看| av网站免费在线观看视频| 久热久热在线精品观看| 成人国产av品久久久| 午夜免费观看性视频| 搡老乐熟女国产| 亚洲性久久影院| 一区二区三区免费毛片| 精品熟女少妇av免费看| 亚洲国产精品成人久久小说| freevideosex欧美| 免费看日本二区| 精品久久久久久电影网| 大码成人一级视频| 国产成人a区在线观看| 一边亲一边摸免费视频| 97超碰精品成人国产| 国产黄片视频在线免费观看| 三级经典国产精品| 国产在视频线精品| 国产高潮美女av| 国产精品99久久99久久久不卡 | 寂寞人妻少妇视频99o| av女优亚洲男人天堂| 亚洲一区二区三区欧美精品| 婷婷色av中文字幕| 国产伦精品一区二区三区视频9| 在线观看人妻少妇| 国产精品免费大片| 尾随美女入室| 一区二区三区四区激情视频| 一个人看的www免费观看视频| 最近手机中文字幕大全| 国产精品一区二区性色av| 国产人妻一区二区三区在| 国产黄色免费在线视频| 国产无遮挡羞羞视频在线观看| 亚洲欧美精品专区久久| 一级av片app| av免费观看日本| 97在线人人人人妻| 亚洲自偷自拍三级|