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

    應(yīng)用渦動相關(guān)技術(shù)的底棲溶解氧通量原位觀測
    --以青島匯泉灣礫石海灘短周期變化為例

    2017-12-27 02:01:23崔尚公于新生趙廣濤
    關(guān)鍵詞:渦動漲潮溶解氧

    崔尚公,于新生,趙廣濤

    (1. 國家海洋信息中心 天津市 300171;2. 中國海洋大學(xué)海洋地球科學(xué)學(xué)院 海底科學(xué)與探測技術(shù)重點實驗室 山東青島 266100)

    應(yīng)用渦動相關(guān)技術(shù)的底棲溶解氧通量原位觀測
    --以青島匯泉灣礫石海灘短周期變化為例

    崔尚公1,2,于新生2,趙廣濤2

    (1. 國家海洋信息中心 天津市 300171;2. 中國海洋大學(xué)海洋地球科學(xué)學(xué)院 海底科學(xué)與探測技術(shù)重點實驗室 山東青島 266100)

    海岸帶區(qū)域內(nèi),穿過沉積物-海水界面的溶解氧通量是研究水體環(huán)境內(nèi)生物地球化學(xué)循環(huán)周轉(zhuǎn)的重要指標。由于傳統(tǒng)采樣方法存在局限性,原位觀測儀器布放難度較大,因此礫石海灘環(huán)境下的溶解氧通量研究還較少。利用非侵入性的渦動相關(guān)技術(shù)對青島匯泉灣潮間帶礫石基底的底棲溶解氧通量做出評估。渦動相關(guān)原位觀測結(jié)果表明,在一個漲潮期內(nèi)溶解氧通量變換范圍為-16.7888 ±5.0~+49.3344 ±3.7mmol O2/m2/d。不穩(wěn)定的通量結(jié)果表明近海潮間帶底棲溶解氧通量在漲潮時段主要受波浪作用影響呈現(xiàn)出極為復(fù)雜的變化形式。頻譜分析表明匯泉灣潮間帶底棲溶解氧主要受東南向波浪作用控制,0.093~0.279 Hz(對應(yīng)時間間隔為3.58~10.75s)頻帶為其主要貢獻頻段。通過對不同時刻累加共譜的比較分析,可以推斷波浪作用和破碎波作用分別在漲潮期前后對底棲溶解氧通量變化起著主導(dǎo)推進作用。此項研究也表明渦動相關(guān)技術(shù)可應(yīng)用于近海潮間帶底棲溶解氧觀測,為從宏觀到微觀研究海岸帶生物地球化學(xué)進程研究提供了新的觀測手段。

    沉積物-海水界面;氧通量;渦動相關(guān);礫石海灘;原位觀測

    海床是水體與沉積物間地球化學(xué)和生物進程緊密耦合的區(qū)域[1]。沿海近岸的沉積物-海水界面既是有機質(zhì)成巖作用的重要場所[2],也是溶質(zhì)和顆粒物交換循環(huán)的活躍地帶。因此,界面間的溶解氧通量作為地球化學(xué)反應(yīng)中的重要參數(shù),是評估海底邊界層碳礦化和初級生產(chǎn)的關(guān)鍵指標之一[3,4]。準確測定溶解氧通量對研究不同時空尺度上生物地球化學(xué)循環(huán)與生態(tài)系統(tǒng)活性至關(guān)重要[5]。

    傳統(tǒng)上,一般使用底棲通量室或微電極剖面法原位測量底邊界層溶解氧通量[6-9]。底棲培養(yǎng)室技術(shù)是將開口箱體嵌入到表層沉積物以下,使得已知面積區(qū)域的沉積物及其上覆水體在海底呈封閉狀態(tài),同時通過嵌入傳感器對封閉上覆水中的如O2、H2S等化學(xué)元素濃度進行實時跟蹤采樣或測量。該方法會封閉區(qū)域沉積物與上覆水,影響水體自然流動和水結(jié)組分的自然交換,使外界動力環(huán)境變化無法在封閉的箱室系統(tǒng)中再現(xiàn)[5]。微電極剖面技術(shù)是采用步進馬達,將溶解氧微電極以微米為單位從上覆水開始緩慢穿過沉積物-水界面逐步刺入沉積物中,檢測不同深度位置的氧濃度,獲得垂直方向氧濃度剖面分布信息,實現(xiàn)沉積物-水界面擴散邊界層中的變化過程測量。這一方法雖然垂直分辨率很高,但僅能獲取一維梯度溶解氧通量信息,無法反映遠離氧微電極穿透區(qū)域的生物擾動和灌洗作用的影響;并且由于微電極尖端非常細,在實際觀測應(yīng)用中極易導(dǎo)致微電極的損壞。近年來,日本及歐美許多發(fā)達國家相繼將大氣科學(xué)中的渦動相關(guān)概念應(yīng)用于水底邊界層界面,利用溶解氧微電極測量距海底15~30cm固定一點的氧濃度,同時使用多普勒矢量流速儀 (Acoustic Doppler Velocimeter, ADV) 對溶解氧濃度測點位置的水流垂向流速進行測量,之后將二者以極高的頻率進行擬合計算氧通量[10]。渦動相關(guān)技術(shù)既沒有封閉沉積物水體,也沒有將微電極插入沉積物,因此可以獲得無人為干擾的原位數(shù)據(jù)。

    美國弗吉尼亞大學(xué)利用渦動相關(guān)系統(tǒng)分別在海草區(qū)、海冰區(qū)和珊瑚礁區(qū)實施了測量,對特殊底質(zhì)區(qū)域的溶解氧通量釋放過程進行了研究[5,11,12]; 瑞 士 聯(lián) 邦水體科學(xué) 與 技 術(shù) 研 究 所 科 學(xué)家使用渦動相關(guān)方法探討了阿爾卑納赫湖受湖震影響的水體氧輸運機理[13]。目前,國內(nèi)外底棲氧通量研究領(lǐng)域的重點主要集中在河口或湖泊有機物沉積區(qū),對堅硬底質(zhì)區(qū)氧通量的研究還未見報道,這主要是由于在堅硬底質(zhì)區(qū)難以進行采樣和儀器布放操作。青島市沿岸具有典型的礫石海灘,遍布藻類及底棲生物,氧交換速率是評估該區(qū)域生物代謝活動的重要參數(shù)。利用渦動相關(guān)技術(shù)在青島匯泉灣礫石海灘對底棲溶解氧通量進行了測量和評估。通過對不同時刻累加共譜的比較分析,推斷出波浪作用和破碎波作用分別在漲潮期前后對底棲溶解氧通量變化起著主導(dǎo)推進作用。

    1 方 法

    1.1 研究區(qū)概況

    研究區(qū)位于膠州灣以東的匯泉灣 (見圖1)。匯泉灣西側(cè)、北側(cè)、東側(cè)三面環(huán)山,以喇叭狀向西南方向開口與黃海相連。海灣高潮時海岸線寬約1.5km,岸灘走向西北-南東,最大面積為1.24km2,平均水深7m。由于山體阻擋了北向來風(fēng),因此灣內(nèi)主要受SSE向風(fēng)浪影響[14]。匯泉灣為典型的砂質(zhì)海岸,自西北向東南沉積物類型由細礫、粗砂、中砂、細砂依次分布變化[15]。海灣潮汐類型屬正規(guī)半日潮,潮流半日周期往復(fù)運動,漲潮時流向為順時針方向,退潮時流向為逆時針方向[16,17]。匯泉灣內(nèi)0~2cm表層沉積物中葉綠素a含量在1.73~5.71μg/g,2~5cm沉積物中葉綠素a含量在0.51~4.32μg/g范圍內(nèi)[18]。

    圖1 研究區(qū)域地形圖(黑色星號為測點區(qū)域)

    圖2 (A)渦動相關(guān)系統(tǒng)現(xiàn)場布設(shè)情況;(B)測點位置情況;(C)長有貽貝類生物的礫石;(D)長有藻類生物的礫石

    測點位置位于匯泉灣東側(cè)中潮帶 (36°03′ 06.38" N, 120°20′23.36"E), 滿 潮 水 深 將 近3.5m。測點附近沉積物類型主要為細砂,砂體之上布有長滿貽貝類和藻類生物的礫石 (見圖2)。

    1.2 渦動相關(guān)系統(tǒng)集成與布設(shè)

    從使用便捷、拆裝簡易、運行安全的角度出發(fā),我們將渦動相關(guān)系統(tǒng)設(shè)計為小型可拆卸式坐底觀測平臺,平臺主要包含溶解氧微電極、聲學(xué)多普勒流速儀、溫鹽深儀三個部分??蚣懿捎媒Y(jié)構(gòu)穩(wěn)定的三角框架,集成了不同的傳感器進行組合檢測。選用挪威Nortek公司生產(chǎn)的聲學(xué)多普勒流速儀 (Acoustic Doppler Velocimeter, ADV)測量固定點三維流速;丹麥Unisense公司生產(chǎn)的溶解氧微電極作為溶解氧傳感器 (型號 OX-100 fast),尖端直徑 90~110μm,90%響應(yīng)時間 <2s,測量精度可達0.3μmol/L;并裝配了一臺意大利IDRONAUT公司生產(chǎn)的CTD多參數(shù)水質(zhì)儀 (型號:OCEAN SEVEN 304Plus)獲取海水深度、溫度和電導(dǎo)率信息,為后期分析氧通量數(shù)據(jù)提供依據(jù)和參考。

    2016年6月4日上午11:00,在匯泉灣東側(cè)潮間帶進行了渦動相關(guān)系統(tǒng)觀測實驗。海流數(shù)據(jù)模塊、溶解氧數(shù)據(jù)模塊和溫鹽深數(shù)據(jù)模塊都被固定在三角支架上 (見圖2),微電極尖端指向ADV測量區(qū)域,測點距沉積物表面24cm,CTD安裝位置稍高,距沉積物表面55cm。觀測過程中,ADV以16Hz的預(yù)設(shè)頻率和20min的預(yù)設(shè)周期進行間歇采樣,循環(huán)周期包括15min持續(xù)采樣時間和5min間歇時間;微電極以16Hz頻率進行連續(xù)采樣;CTD以8Hz的預(yù)設(shè)頻率進行連續(xù)采樣。

    布設(shè)按照原定計劃于低潮時 (上午10:20)開始進行,安放三角支架時特別地將其開口方向正對海流主流速方向以避免支架腳對海流的干擾。數(shù)據(jù)采集初始階段,由于潮水并未完全淹沒測點及傳感器,因此所測數(shù)據(jù)無效,將這部分數(shù)據(jù)歸為 “空采數(shù)據(jù)”,在后期處理中予以剔除。實驗原計劃水下工作10~12h,以獲取完整潮汐漲落周期內(nèi)溶解氧通量變化數(shù)據(jù),但由于中途數(shù)據(jù)記錄發(fā)生故障,系統(tǒng)采集被迫中止。綜合以上兩點,本次觀測采集到了從上午11:20~下午15:00約4h的有效數(shù)據(jù)。據(jù)國家海洋局北海預(yù)報中心當日發(fā)布的潮汐預(yù)報,當日漲潮期為上午10:18~下午15:36,由此可見本次采集數(shù)據(jù)基本覆蓋了漲潮期。

    1.3 氧通量的計算方法

    1.3.1 基本原理

    根據(jù)Reynolds分解理論,在流體力學(xué)研究中常常對溶質(zhì)濃度和流速做出如下分解[19,20]:

    經(jīng)過幾十年來科學(xué)家對水底邊界層中懸浮顆粒和相關(guān)物質(zhì)遷移的研究,發(fā)現(xiàn)可以通過三維溶質(zhì)質(zhì)量守恒方程來確定其分布狀態(tài)[21],即溶質(zhì) (例如溶解氧)在均勻水底邊界層中滿足:

    公式 (4)表示在高度h處測得的溶解氧垂向通量等于沉積物-海水界面垂向通量與厚度為h的水體內(nèi)溶解氧源匯效應(yīng)之和 (w為垂向流速)。當測量點位距沉積物-海水界面很近 (即h很小)時,源匯效應(yīng)項也可以忽略不計。那么就有:

    公式 (5)表示可以用沉積物-海水界面以上可忽略高度的一點的溶解氧垂向通量近似代表界面通量。

    1.3.2 數(shù)據(jù)處理分析

    為簡化數(shù)據(jù)處理,首先將采集到的16Hz頻率數(shù)據(jù)轉(zhuǎn)換為8Hz數(shù)據(jù),并對溶解氧和海流流速數(shù)據(jù)進行了去尖峰濾波、去噪濾波處理,對流速數(shù)據(jù)進行了傾斜校正。之后,使用滑動平均方法進行湍流波動計算。為選擇出合適的濾波窗口長度Nr,采用累加計算方法,將初始窗口長度設(shè)為1,并逐漸增大窗口長度,反復(fù)計算滑動平均值、波動值和通量值,當窗口長度增大到一定程度后,通量值達到近乎穩(wěn)定 (見圖3),以上午12:00~12:15這一數(shù)據(jù)段為例,將3000個數(shù)據(jù)點 (對應(yīng)時間為6.25min)作為最終的濾波窗口長度,以此計算湍流波動值C′和u′。隨后對每個數(shù)據(jù)段分別計算兩個波動值的互相關(guān)函數(shù)(見圖4),以下午13:00~13:15這一數(shù)據(jù)段為例,湍流波動值C′與u′的最大互相關(guān)出現(xiàn)在第7175點位置,距7200點的中心位置偏差了25個數(shù)據(jù)點 (對應(yīng)時間為3.125s),因此將溶解氧濃度和垂向波動速度時間序列相對移動25個數(shù)據(jù)點,使最大互相關(guān)值位于中心位置。隨后根據(jù)式(6)計算功率譜密度并進行頻率響應(yīng)校正。

    最后計算出每個15min時間段內(nèi)的溶解氧通量。

    2 結(jié) 果

    圖3 12:00~12:15數(shù)據(jù)段滑動平均濾波窗口選取示意圖

    圖4 13:00~13:15數(shù)據(jù)段互相關(guān)校正圖

    經(jīng)過處理和計算得到的海流垂向流速、溶解氧濃度、水深、溫度及溶解氧通量數(shù)據(jù)見圖5和表1。由于CTD安裝位置稍高,導(dǎo)致被海水淹沒前的"空采期"較長,因此舍去了水深和水溫的前20min數(shù)據(jù);另外CTD供電系統(tǒng)在13:40左右發(fā)生故障,因此水溫的后1h數(shù)據(jù)未在圖2-8中顯示,水深數(shù)據(jù)由ADV測量值代替補充。測量期間隨著潮水的上漲,水深變化范圍為0.3~2.8m,平均水深1.68m;水溫在近1.5h內(nèi)從17.58℃持續(xù)下降到16.82℃,據(jù)ADV溫度傳感器記錄數(shù)據(jù)顯示,當日下午14:00~17:00期間水溫依然在下降,但變化幅度極小,17:00時刻水溫為16.48℃。海流水平流速變化范圍為10.380~25.031cm/s,平均水平流速 19.43cm/s;垂直流速變化范圍為0.012~0.686cm/s,平均水平流速為0.36cm/s。

    圖5 渦動相關(guān)數(shù)據(jù)時序變化特征:(A)海流垂向速度;(B)溶解氧濃度;(C)水深;(D)水溫;(E)計算得到的溶解氧通量(正通量表示通量方向為從沉積物指向上覆水)(誤差條代表標準差)

    測量過程中,溶解氧濃度變化范圍為292.80 umol/L到248.13 ?mol/L。在測量初始階段,由于潮水上漲,測點逐漸被海水淹沒,溶解氧濃度快速下降。在12:10左右,溶解氧濃度出現(xiàn)小幅度回升,這可能是由于潮間帶沉積物在被海水淹沒之前暴露于空氣中,沉積物中的間隙水較少,因此溶解態(tài)的氧也相對較少,而漲潮后沉積物被海水覆蓋,沉積物間隙中的氧氣快速溶于海水,使溶解氧濃度增加。在此之后,溶解氧濃度隨著水深的增長持續(xù)下降,這主要是由于測點位置固定不變,測點與海水-大氣界面距離逐漸擴大,因此溶解氧濃度呈緩慢下降趨勢。在下午14:30后,水深達到2.5m,深度對溶解氧濃度的影響逐漸減弱,溶解氧濃度基本保持穩(wěn)定,穩(wěn)定值約為252.89 ?mol/L。

    表1給出了實測流速、水深、水溫、溶解氧濃度及溶解氧通量計算結(jié)果的具體數(shù)據(jù)。從中可以看出溶解氧通量存在明顯的規(guī)律性。在上午11:40到11:55的觀測初始階段,溶解氧通量為負值,這是由于在漲潮初期水深較淺,測點位置主要受到海水-大氣界面作用影響,近海面氧濃度顯著高于沉積物表層,造成了溶解氧自上而下向沉積物運動。在12:00到12:35這一時段,溶解氧通量轉(zhuǎn)為正值。如上文所述,漲潮后潮間帶沉積物被海水覆蓋,沉積物間隙中的氧氣快速溶于海水,使沉積物間隙溶解氧濃度增加,高于海水中氧濃度,形成了下高上低的濃度梯度分布。此后,氧通量緩慢下降,逐漸由正值轉(zhuǎn)為負值,并在基本達到滿潮時開始回升,重新變?yōu)檎?,最終達到+49.3 ±3.7 mmol/m2/d。這可能由幾種機制造成:首先,隨著潮水上漲,水深逐漸增大,測點與海水-大氣界面的距離擴大,形成了從海水表面向下指向沉積物的垂向通量;其次,海流的作用造成了上高下低的氧濃度梯度。據(jù)Holtappels M等人研究表明,近岸海域中海流流速隨著距海底距離的增加而增長,這種流速分布會破環(huán)近岸水體中原有的氧濃度梯度,產(chǎn)生或增強自上而下方向的溶解氧通量 (見圖6)[23]。第三,當水深達到2.5m左右時,海水-大氣界面作用對溶解氧濃度的影響逐漸減弱消失,測點位置水體環(huán)境基本保持穩(wěn)定。本次研究測區(qū)海底礫石和沉積物上長有大量貽貝類和底棲藻類生物,底部藻類光合作用會生成大量的氧,形成自下而上的溶解氧通量。另外,滲透性砂質(zhì)沉積物在海流的反復(fù)沖刷過程中,不斷地被重置與翻新。氧和固態(tài)、溶解態(tài)有機質(zhì)被向下帶入滲透性沉積物,增強了微生物的礦化作用,同時缺氧水被沖出沉積物混入上覆水[24,25]。此時的砂質(zhì)沉積物便成為了生物催化作用的過濾器和催化劑,刺激了碳循環(huán)和氧循環(huán)的周轉(zhuǎn)[26],造成了較高的底

    棲氧通量值。

    表1 原位渦動相關(guān)實驗詳細數(shù)據(jù)(括號內(nèi)為平均值)(正通量表示通量方向為從沉積物指向上覆水)

    圖6 近岸海流流速分布引起的垂向濃度梯度示意

    為進一步評估渦動相關(guān)數(shù)據(jù)和所得通量的質(zhì)量,分別對每個15min數(shù)據(jù)段 (burst)進行累加通量計算,計算結(jié)果見圖7。累加通量的斜率代表了所測時間段內(nèi)的平均氧通量值,可以根據(jù)其升降趨勢評估計算的氧通量的細微變化。從圖中可以看出,在漲潮期的大部分時間內(nèi),累加通量都沒有顯示出完美的線性趨勢,這表明研究站位的水動力條件極度不穩(wěn)定[5]。隨著海流流速的增長,在不平坦的砂質(zhì)沉積物表面反復(fù)沖刷的海流引發(fā)的壓力差會在上層砂質(zhì)層中產(chǎn)生強勁的孔隙水流[28]。下午14:30后接近滿潮時,累加通量開始呈現(xiàn)較為線性的走勢,表明此時水動力條件趨于穩(wěn)定,所測氧通量值對該時段真實氧通量具有很好的代表性。

    圖7 觀測過程中波動的累加通量

    3 分析與討論

    3.1 氧通量質(zhì)量評估

    根據(jù)匯泉灣測點附近底質(zhì)條件、測點固定高度估算了本次實驗指示區(qū)域范圍,區(qū)域長約為55.275m;寬1.567m;最大貢獻點位于海流上游距測點2.512m位置。本次實驗由于觀測時間處于漲潮期階段,因此測得的溶解氧通量極度不穩(wěn)定,變化范圍為 -16.7888±5.0~+49.3344±3.7mmol O2/m2/d, 平均值為+11.3876 mmol O2/m2/d。這與美國科學(xué)家AC Hume利用渦動相關(guān)方法在海草生長區(qū)測得的+121 mmol O2/m2/d和美國科學(xué)家Berg在西法爾茅斯河口區(qū)域測得的+77.5 mmol O2/m2/d的溶解氧通量相比數(shù)值較低。這可能有以下幾個原因。首先,研究區(qū)域存在較大差異。本次實驗區(qū)域位于匯泉灣潮間帶,匯泉灣內(nèi)無河流入海,砂源貧乏,加之測點區(qū)域的砂質(zhì)沉積物上滿布礫石,雖然存在貽貝類和藻類生物,但沉積物內(nèi)有機質(zhì)含量和指示區(qū)域內(nèi)海草密度都無法與海草和河口區(qū)域相比。其次,2016年6月4日當天陰有小雨,較弱的太陽輻照度減弱了底棲植物的光合作用。第三,波浪和海流引發(fā)的沉積物再懸浮使海水濁度增加,也會對底棲植物光合作用產(chǎn)生負面影響。

    另外,通過對所得漲潮期時段內(nèi)溶解氧通量和水平、垂向流速序列的相關(guān)系數(shù)分析,發(fā)現(xiàn)氧通量和水平流速序列相關(guān)系數(shù)只有0.084,與垂向流速序列相關(guān)系數(shù)為0.398,表明本次實驗測得的溶解氧通量主要受垂向流速影響,這與渦動相關(guān)理論依據(jù)是相符的。

    3.2 頻譜分析

    為進一步分析海流和波浪作用對溶解氧通量的影響,分別計算了溶解氧濃度、垂向速度的功率譜密度及二者的共功率譜密度和累加共譜。功率譜密度函數(shù)指示了數(shù)據(jù)中不同頻率成分所占百分比。典型的功率譜結(jié)果見圖8。對于觀測起始的11:40~11:55數(shù)據(jù)段 (burst 1),可以看出0.093~0.279 Hz頻段 (對應(yīng)時間間隔為3.58~10.75s)是總功率譜的主要貢獻頻段,這一時間間隔與匯泉灣歷史記錄的5.4~11s的平均波浪周期極為吻合[14]。據(jù)Lorrai C等人研究發(fā)現(xiàn),當沉積物之上1m水平流速達到20cm/s,該區(qū)域可被視為高度湍流區(qū)域,其渦動時間尺度約為0.2~10 s[29]。這也與此次實驗結(jié)果極為一致。不過也可以發(fā)現(xiàn)這一時段內(nèi)頻率大于1 Hz的高頻成分占比很低,僅占15.59%。對于觀測過程中間階段的13:20~13:35數(shù)據(jù)段 (burst 6),雖然0.1~0.3Hz頻段貢獻依然突出,但可以看出其大于1 Hz的高頻成分占比顯著增長,達到63.57%。另外,0.4~0.62 Hz頻段 (對應(yīng)時間間隔為1.61~2.50s)成為一新的主要貢獻頻段。累加共譜表示了大于其頻率的頻段貢獻所占百分比,因此隨著頻率逐漸減小,其百分比值由零逐漸趨近于一。從圖8(D)可以看出,漲潮期間隨著海流和波浪作用的加劇,高頻成分貢獻逐漸擴大。根據(jù)匯泉灣波浪周期的歷史統(tǒng)計信息,可以推斷夏季東南向波浪作用是該測區(qū)溶解氧動態(tài)變化的主導(dǎo)因素。

    圖8 burst1和burst6功率譜密度計算結(jié)果(A)溶解氧濃度歸一化功率譜;(B)垂向流速歸一化功率譜;(C)共功率譜密度;(D) 累加共譜

    通過比較每個時間段的累加共譜曲線,可以看出隨著時間推移,數(shù)據(jù)主要貢獻頻帶逐漸向高頻移動 (見圖9)。對于burst1時間段,共譜區(qū)間為一條凹型曲線,主要貢獻頻段為0.093~0.279 Hz(對應(yīng)時間間隔為3.58~10.75s)。對于burst3、7、10時間段來說,曲線形狀逐漸變?yōu)橥剐颓€,這表明在漲潮期開始的2h內(nèi)高頻成分的影響快速擴大并在接近滿潮時逐漸趨于穩(wěn)定。大于0.8Hz的頻帶占比從20.17% 增長到71.40%,湍流引發(fā)的極度不穩(wěn)定的水動力條件導(dǎo)致了溶解氧通量的大幅度波動。

    圖9 burst1,burst3,burst7和burst10四個時段累加共譜圖

    5 結(jié) 論

    基于Reynolds分解理論和Taylor凍結(jié)流假說,利用聲學(xué)多普勒流速儀 (ADV)、溶解氧微電極和溫鹽深儀 (CTD)構(gòu)建了渦動相關(guān)原位觀測系統(tǒng),并將其應(yīng)用在青島近海潮間帶礫石海灘,對溶解氧通量進行了原位觀測。觀測結(jié)果表明,在一個漲潮期內(nèi)溶解氧通量變換范圍為 -16.7888±5.0~+49.3344±3.7mmol O2/m2/d。不穩(wěn)定的通量結(jié)果表明近海潮間帶底棲溶解氧通量在漲潮時段會受到海水上漲、海流流速加快等物理因素和間隙海水溶解、植物光合作用等生物化學(xué)因素的共同作用,呈現(xiàn)出極為復(fù)雜的變化形式。頻譜分析表明匯泉灣潮間帶底棲溶解氧主要受東南向波浪作用控制,0.093~0.279 Hz(對應(yīng)時間間隔為3.58~10.75s)頻帶為其主要貢獻頻段。通過對不同時刻累加共譜的比較分析,可以推斷波浪作用和破碎波作用分別在漲潮期前后對底棲溶解氧通量變化起著主導(dǎo)推進作用。本次實驗也驗證了渦動相關(guān)系統(tǒng)應(yīng)用于近海潮間帶底棲溶解氧觀測的可行性,為進一步了解海底邊界層有機質(zhì)的生產(chǎn)分解,開展早期成巖規(guī)律和生物地球化學(xué)進程研究提供了準確的科學(xué)依據(jù)和可靠的技術(shù)支持。

    [1] HOPKINSON CS, WETZEL RL. In situ Measurements of Nutrient and Oxygen Fluxes in a Coastal Marine Benthic Community[J]. Marine Ecology Progress, 1982, 10,29-35.

    [2] MIDDELBURG JJ, CARLOS CM, GATTUSO JJ. Respiration in coastal benthic communities[M]//Williams Peter, JleB, del G, Paul A. Respiration in aquatic ecosystems. Oxford University Press, New York, 2008: 206-224.

    [3] CANFIELD DE, THAMDRUP B, HANSEN JW. The anaerobic degradation of organic matter in Danish coastal sediments: iron reduction, manganese reduction, and sulfate reduction[J]. Geochimica Et Cosmochimica Acta,1993, 57, 3867-83.

    [4] GLUD RN. Oxygen dynamics of marine sediments[J].Marine Biology Research, 2008, 4, 243-289.

    [5] BERG P, GLUD RN, HUME A, et al. Eddy correlation measurements of oxygen uptake in deep ocean sediments[J]. Limnology and Oceanography Methods, 2009, 7,576-584.

    [6] JORGENSEN BB, BOUDREAU BP. Diagenesis and sediment-water exchange[M]//Boudreau B, J rgensen BB.The benthic boundary layer. Oxford University Press,New York, 2001: 221-244.

    [7] REIMERS CE, JAHNKE RA, THOMSEN L. In situ sampling in the benthic boundary layer[M]//Boudreau B,J0rgensen BB. The benthic boundary layer. Oxford University Press, New York, 2001: 245-268.

    [8] WENZHOFER F, GLUD RN. Small-scale spatial and temporal variability in coastal benthic O2dynamics: effects of fauna activity[J]. Limnology and Oceanography, 2004,49, 1471-1481.

    [9] BEER DD, WENZHOFER F, FERDELMAN TG et al.Transport and mineralization rates in North Sea sandy intertidal sediments, Sylt-R?m? Basin, Wadden Sea[J].Limnology and Oceanography, 2005, 50, 113-127.

    [10] BERG P, ROY H, JANSSEN F, et al. Oxygen uptake by aquatic sediments measured with a novel non-invasive eddy-correlation technique[J]. Marine Ecology Progress,2003, 261, 75-83.

    [11] HUME AC, BERG P, MCGLATHERY KJ. Dissolved oxygen fluxes and ecosystem metabolism in an eelgrass(Zostera marina) meadow measured with the eddy correlation technique[J]. Limnology and Oceanography, 2011,56, 86-96.

    [12]LONG MH, KOOPMANS D, BERG P, et al. Oxygen exchange and ice melt measured at the ice-water interface by eddy correlation[J]. Biogeosciences Discussions, 2012,8, 11255-11284.

    [13] BRAND A, MCGINNIS DF, WEHRLI B,et al. Intermittent oxygen flux from the interior into the bottom boundary of lakes as observed by eddy correlation[J]. Limnology and Oceanography, 2008, 53, 1997-2006.

    [14] 常瑞芳,劉鐳,范元炳. 波浪對青島匯泉灣潮間沙壩的塑造作用[J]. 中國海洋大學(xué)學(xué)報自然科學(xué)版, 1992,(4):61-70.

    [15] 王偉偉, 莊麗華, 閻軍, 等. 青島市匯泉灣海水 浴場 表層沉積物粒度特征及輸運趨勢[J]. 中國石油大學(xué)學(xué)報自然科學(xué)版, 2007, 31(3): 13-17.

    [16] 鄭建瑜, 且鐘禹, 李學(xué)倫. 青島南海岸海水浴場的旅游環(huán)境質(zhì)量評價[J]. 海洋環(huán)境科學(xué), 1998, (1): 66-72.

    [17] 李春柱, 徐宏明, 侯國本. 青島匯泉灣浴場部分海灘潰蝕及治理的初步探討[J]. 海岸工程, 1992, (1): 26-31.

    [18] 吳以平, 劉曉收. 青島灣潮間帶沉積物中葉綠素的分析[J]. 海洋科學(xué), 2005, 29 (11): 8-12.

    [19] REYNOLDS O. On the Dynamical Theory of Incompressible Viscous Fluids and the Determination of the Criterion[J]. Philosophical Transactions of the Royal Society of London A, 1995, 186 (1941): 123-164.

    [20] STANISIC MM. The mathematical theory of turbulence[J]. Universitext, 1988, 41 (1): 68.

    [21] HILL PS,NOWELL ARM. Comparison of two models of aggregation in continental-shelf bottom boundary layers[J]. Journal of Geophysical Research Oceans, 1995,100 (C11): 22749-22763.

    [22] TAYLOR GI. The spectrum of turbulence[J]. Proceedings of the Royal Society of London, 1938, 164 (919):476-409.

    [23] HOLTAPPELS M, GLUD RN, DONIS D, et al. Effects of transient bottom water currents and oxygen concentrations on benthic exchange rates as assessed by eddy correlation measurements[J]. Journal of Geophysical Research Oceans, 2013, 118 (3): 1157-1169.

    [24] PRECHT E., HUETTEL M. Advective pore-water exchange driven by surface gravity waves and its ecological implications[J]. Limnology & Oceanography, 2003, 48:1674-1684.

    [25] PRECHT E, FRANKE U, POLERECKY L, et al. Oxygen dynamics in permeable sediments with wave-driven pore water exchange[J]. Limnology & Oceanography,2004, 49 (3): 693-705.

    [26] HUETTEL M, WEBSTER IT. Porewater flow in permeablesediments[M]//Boudreau B, J?rgensen BB. The benthic boundary layer. Oxford University Press, New York,2001: 144-179.

    [27] HOLTAPPELS M, GLUD RN, DONIS D, et al. Effects of transient bottom water currents and oxygen concentrations on benthic exchange rates as assessed by eddy correlation measurements[J]. Journal of Geophysical Research Oceans, 2013, 118 (3): 1157-1169.

    [28] BERG P, LONG MH, HUETTEL M, et al. Eddy correlation measurements of oxygen fluxes in permeable sediments exposed to varying current flow and light[J]. Limnology & Oceanography, 2013, 58 (4): 1329-1343.

    [29] LORRAI C, MCGINNIS DF, BERG P, et al. Application of Oxygen Eddy Correlation in Aquatic Systems[J]. Journal of Atmospheric & Oceanic Technology, 2010, 27(9):1533-1546.

    2017-11-21

    國家自然科學(xué)基金資助項目(41276089;41176078)。

    猜你喜歡
    渦動漲潮溶解氧
    漲潮
    散文詩世界(2023年3期)2023-04-15 17:12:36
    淺析水中溶解氧的測定
    “水城”被淹
    污水活性污泥處理過程的溶解氧增益調(diào)度控制
    城市河道洲灘對水流溶解氧分布的影響
    BTA鉆桿渦動數(shù)學(xué)建模及實驗研究
    神奇的潮汐
    理想條件下BTA鉆鉆桿的渦動分析
    漲潮啦
    基于LabVIEW的油膜渦動及油膜振蕩故障特征分析研究
    亚洲最大成人中文| 成人影院久久| 毛片女人毛片| 成年av动漫网址| 国产伦理片在线播放av一区| freevideosex欧美| 在线观看三级黄色| 一区二区三区免费毛片| 伊人久久国产一区二区| 看十八女毛片水多多多| 日韩不卡一区二区三区视频在线| 99久久综合免费| 国产v大片淫在线免费观看| 亚洲精品一二三| 两个人的视频大全免费| 成人特级av手机在线观看| 在现免费观看毛片| 精品一区二区三卡| 九九在线视频观看精品| 精品久久久久久电影网| 久久 成人 亚洲| 国产v大片淫在线免费观看| 麻豆成人av视频| 婷婷色麻豆天堂久久| 久久久国产一区二区| 日韩免费高清中文字幕av| 亚洲图色成人| 国产熟女欧美一区二区| 我的女老师完整版在线观看| 高清欧美精品videossex| 欧美极品一区二区三区四区| 伦理电影大哥的女人| 欧美日韩精品成人综合77777| 深爱激情五月婷婷| 国产视频内射| 伦精品一区二区三区| 蜜臀久久99精品久久宅男| 日本猛色少妇xxxxx猛交久久| 亚洲四区av| 免费黄频网站在线观看国产| 制服丝袜香蕉在线| 久久久久久久久大av| 一区二区三区四区激情视频| 久久久久久久久久久丰满| 国产日韩欧美在线精品| 午夜福利网站1000一区二区三区| 少妇的逼好多水| 一本一本综合久久| 精品亚洲乱码少妇综合久久| 亚洲欧美精品自产自拍| 久久久精品94久久精品| 婷婷色av中文字幕| 一个人看的www免费观看视频| 久久精品久久久久久噜噜老黄| 国产高清有码在线观看视频| 亚洲精品国产av蜜桃| 国产成人freesex在线| 日本av免费视频播放| 一级毛片久久久久久久久女| 午夜精品国产一区二区电影| 亚洲成人av在线免费| 亚洲欧美成人综合另类久久久| 女人十人毛片免费观看3o分钟| 久热久热在线精品观看| 国产成人freesex在线| 婷婷色综合大香蕉| 免费观看性生交大片5| 男人狂女人下面高潮的视频| 成人免费观看视频高清| 亚洲精品乱久久久久久| 最近的中文字幕免费完整| av视频免费观看在线观看| 日韩av在线免费看完整版不卡| 丰满乱子伦码专区| 嘟嘟电影网在线观看| 不卡视频在线观看欧美| 国内揄拍国产精品人妻在线| 我要看黄色一级片免费的| 亚洲人成网站高清观看| 亚洲精品自拍成人| 女人十人毛片免费观看3o分钟| 熟妇人妻不卡中文字幕| 欧美xxⅹ黑人| 亚洲内射少妇av| 欧美精品亚洲一区二区| 尾随美女入室| 亚洲伊人久久精品综合| 99久久精品国产国产毛片| 亚洲av成人精品一区久久| 亚洲第一区二区三区不卡| 国产成人精品婷婷| 国产69精品久久久久777片| 亚洲精华国产精华液的使用体验| 最黄视频免费看| 国产免费视频播放在线视频| 观看免费一级毛片| 国产精品国产三级专区第一集| 成人漫画全彩无遮挡| 一级av片app| 大香蕉久久网| 中文精品一卡2卡3卡4更新| 久久久久网色| 精品少妇黑人巨大在线播放| 日日撸夜夜添| 久久久久性生活片| 网址你懂的国产日韩在线| 国产成人91sexporn| 一本久久精品| 日韩在线高清观看一区二区三区| 日韩制服骚丝袜av| 熟女人妻精品中文字幕| 色综合色国产| 少妇的逼好多水| 少妇熟女欧美另类| 免费av不卡在线播放| 色吧在线观看| 精品一区二区三区视频在线| 少妇的逼水好多| 国产精品无大码| 一二三四中文在线观看免费高清| 久久国产精品男人的天堂亚洲 | videossex国产| 国产一区有黄有色的免费视频| av在线蜜桃| 国产精品欧美亚洲77777| 日韩一区二区三区影片| 人人妻人人澡人人爽人人夜夜| 成人国产av品久久久| 国产精品精品国产色婷婷| 极品教师在线视频| 国产成人精品久久久久久| 国国产精品蜜臀av免费| 内地一区二区视频在线| 天天躁夜夜躁狠狠久久av| 中国国产av一级| 91精品国产九色| 亚洲色图av天堂| 成人国产av品久久久| av卡一久久| 又黄又爽又刺激的免费视频.| 成人毛片a级毛片在线播放| 日韩强制内射视频| 97在线人人人人妻| 男的添女的下面高潮视频| 丰满迷人的少妇在线观看| 狠狠精品人妻久久久久久综合| 黑丝袜美女国产一区| 亚洲精品色激情综合| 精品人妻熟女av久视频| 成人亚洲精品一区在线观看 | 国产深夜福利视频在线观看| 国产精品99久久久久久久久| 国产永久视频网站| 日本色播在线视频| 久久久久久久久久人人人人人人| 国产毛片在线视频| 我的老师免费观看完整版| 国内揄拍国产精品人妻在线| 校园人妻丝袜中文字幕| 少妇熟女欧美另类| 亚洲国产毛片av蜜桃av| 亚洲在久久综合| 天天躁夜夜躁狠狠久久av| 777米奇影视久久| 国模一区二区三区四区视频| 18禁动态无遮挡网站| 97超碰精品成人国产| av网站免费在线观看视频| 日韩大片免费观看网站| 欧美成人a在线观看| 久久久久久久大尺度免费视频| 男人狂女人下面高潮的视频| 免费看不卡的av| 亚洲激情五月婷婷啪啪| 九色成人免费人妻av| 一区二区三区乱码不卡18| 成年av动漫网址| 国产成人freesex在线| 午夜福利在线在线| 国产亚洲午夜精品一区二区久久| 欧美成人午夜免费资源| 91久久精品电影网| 91精品伊人久久大香线蕉| 免费人成在线观看视频色| 少妇人妻一区二区三区视频| 国产精品一区www在线观看| 免费看光身美女| 国产精品国产三级国产av玫瑰| 人妻系列 视频| 亚洲av二区三区四区| 一级片'在线观看视频| 韩国高清视频一区二区三区| 制服丝袜香蕉在线| 国产精品久久久久久av不卡| 亚洲精品一二三| 国产成人91sexporn| 午夜激情福利司机影院| 国产av精品麻豆| 国产亚洲最大av| 内地一区二区视频在线| 黄片无遮挡物在线观看| 久久久久网色| 色网站视频免费| 97在线视频观看| 九九爱精品视频在线观看| 精品人妻偷拍中文字幕| 春色校园在线视频观看| 看免费成人av毛片| 精品一区在线观看国产| 亚洲av成人精品一区久久| 日本av免费视频播放| 日本一二三区视频观看| 国产白丝娇喘喷水9色精品| 国产精品欧美亚洲77777| 免费黄网站久久成人精品| 边亲边吃奶的免费视频| 亚洲va在线va天堂va国产| 亚洲精华国产精华液的使用体验| 亚洲精品色激情综合| 美女视频免费永久观看网站| 国产乱人视频| 能在线免费看毛片的网站| 尾随美女入室| 国产精品欧美亚洲77777| 亚洲精品乱码久久久v下载方式| 久久久久久九九精品二区国产| 久久韩国三级中文字幕| 国产av一区二区精品久久 | 99久久综合免费| 精品一品国产午夜福利视频| 亚洲综合色惰| 青春草国产在线视频| 欧美区成人在线视频| 久久综合国产亚洲精品| 欧美97在线视频| 婷婷色综合大香蕉| 美女高潮的动态| 亚洲第一区二区三区不卡| 日本爱情动作片www.在线观看| 97在线人人人人妻| 在线观看一区二区三区| 久久精品久久久久久久性| 精品熟女少妇av免费看| 亚洲色图综合在线观看| 亚洲av成人精品一区久久| 日韩中字成人| 麻豆乱淫一区二区| 日韩欧美 国产精品| 日日摸夜夜添夜夜爱| 日韩欧美一区视频在线观看 | 国模一区二区三区四区视频| 国产在线免费精品| 青青草视频在线视频观看| 一级爰片在线观看| 亚洲伊人久久精品综合| 成年av动漫网址| av在线观看视频网站免费| 亚洲欧美清纯卡通| 一个人看的www免费观看视频| 我要看日韩黄色一级片| 精品一品国产午夜福利视频| 美女国产视频在线观看| 国产精品一二三区在线看| 亚洲av二区三区四区| 久久ye,这里只有精品| 精品亚洲乱码少妇综合久久| 又大又黄又爽视频免费| 极品教师在线视频| 亚洲欧美日韩无卡精品| 亚洲欧美日韩另类电影网站 | av卡一久久| 国产成人一区二区在线| 成人18禁高潮啪啪吃奶动态图 | 狂野欧美激情性bbbbbb| 在线看a的网站| 美女内射精品一级片tv| 女性生殖器流出的白浆| 视频中文字幕在线观看| 欧美日本视频| 亚洲精品国产色婷婷电影| 亚洲av成人精品一区久久| 免费播放大片免费观看视频在线观看| 男人狂女人下面高潮的视频| 欧美精品一区二区免费开放| 成年免费大片在线观看| 国产极品天堂在线| 97精品久久久久久久久久精品| 国产有黄有色有爽视频| 卡戴珊不雅视频在线播放| 亚洲精品一二三| 久久人人爽人人爽人人片va| 久久国产精品男人的天堂亚洲 | 免费看不卡的av| 又粗又硬又长又爽又黄的视频| 久久精品国产a三级三级三级| 亚洲精品,欧美精品| 国产人妻一区二区三区在| 99精国产麻豆久久婷婷| 成人无遮挡网站| 国产精品一区二区三区四区免费观看| 啦啦啦中文免费视频观看日本| 午夜免费男女啪啪视频观看| 丰满乱子伦码专区| 久久毛片免费看一区二区三区| 欧美精品国产亚洲| 亚洲人成网站高清观看| 国产一区亚洲一区在线观看| 少妇精品久久久久久久| 久久久久视频综合| av在线app专区| 久久精品国产自在天天线| 在线 av 中文字幕| 永久免费av网站大全| 中文字幕精品免费在线观看视频 | 国产 精品1| 日本-黄色视频高清免费观看| 六月丁香七月| 一个人看的www免费观看视频| 最近中文字幕高清免费大全6| 亚洲欧美成人精品一区二区| 亚洲av电影在线观看一区二区三区| av专区在线播放| 国产成人精品婷婷| 黄色欧美视频在线观看| 欧美精品国产亚洲| 国产淫语在线视频| 日韩av免费高清视频| 老师上课跳d突然被开到最大视频| 免费播放大片免费观看视频在线观看| 免费黄网站久久成人精品| 免费观看a级毛片全部| 欧美精品亚洲一区二区| 精品久久久久久电影网| 免费黄网站久久成人精品| 日韩中字成人| 九九爱精品视频在线观看| 3wmmmm亚洲av在线观看| 国产欧美日韩一区二区三区在线 | 91午夜精品亚洲一区二区三区| 老司机影院成人| 欧美xxⅹ黑人| 国产精品久久久久久精品电影小说 | 一级毛片电影观看| 久久久欧美国产精品| 久久久精品94久久精品| 亚洲欧美成人精品一区二区| 国产精品女同一区二区软件| 国产成人一区二区在线| av在线播放精品| 99国产精品免费福利视频| 偷拍熟女少妇极品色| 91精品一卡2卡3卡4卡| 国产精品爽爽va在线观看网站| 秋霞伦理黄片| 狂野欧美激情性xxxx在线观看| 午夜视频国产福利| 激情 狠狠 欧美| 午夜免费男女啪啪视频观看| 丰满迷人的少妇在线观看| 亚洲精品一二三| 日韩强制内射视频| 在线精品无人区一区二区三 | 美女脱内裤让男人舔精品视频| 嫩草影院入口| 一级毛片黄色毛片免费观看视频| 日本欧美国产在线视频| 久久热精品热| 建设人人有责人人尽责人人享有的 | 男女边吃奶边做爰视频| 亚洲人成网站在线观看播放| 一级二级三级毛片免费看| 日韩成人伦理影院| 在线观看美女被高潮喷水网站| 尾随美女入室| 一区二区三区四区激情视频| 波野结衣二区三区在线| 免费观看的影片在线观看| 亚洲电影在线观看av| 我要看黄色一级片免费的| 天堂中文最新版在线下载| 毛片女人毛片| 午夜福利高清视频| 99热6这里只有精品| 91午夜精品亚洲一区二区三区| 少妇 在线观看| 婷婷色麻豆天堂久久| 人妻系列 视频| 黄片无遮挡物在线观看| 国产爽快片一区二区三区| 草草在线视频免费看| 亚洲精品乱码久久久v下载方式| 国产精品秋霞免费鲁丝片| 熟妇人妻不卡中文字幕| av免费在线看不卡| 日韩电影二区| 舔av片在线| 久久久久久久久久成人| 欧美zozozo另类| .国产精品久久| 日韩免费高清中文字幕av| 永久网站在线| 久久久久国产精品人妻一区二区| 国产黄频视频在线观看| 国产91av在线免费观看| 日韩大片免费观看网站| 日韩一区二区视频免费看| 成人特级av手机在线观看| 国产精品不卡视频一区二区| 久久6这里有精品| 欧美人与善性xxx| 校园人妻丝袜中文字幕| 夫妻性生交免费视频一级片| 亚洲四区av| 午夜福利高清视频| 免费看av在线观看网站| 亚洲国产av新网站| 一级爰片在线观看| 亚洲精品国产色婷婷电影| 人人妻人人爽人人添夜夜欢视频 | 国产亚洲精品久久久com| 制服丝袜香蕉在线| 秋霞在线观看毛片| av国产免费在线观看| 免费黄网站久久成人精品| 日日撸夜夜添| av专区在线播放| 成人18禁高潮啪啪吃奶动态图 | 亚洲图色成人| 自拍欧美九色日韩亚洲蝌蚪91 | 汤姆久久久久久久影院中文字幕| 久热久热在线精品观看| 99热网站在线观看| 大又大粗又爽又黄少妇毛片口| h日本视频在线播放| 国产精品免费大片| 久久久久久久久久成人| 王馨瑶露胸无遮挡在线观看| 在线看a的网站| 午夜激情福利司机影院| 欧美日韩精品成人综合77777| av线在线观看网站| 久久久国产一区二区| 久久久欧美国产精品| 最近中文字幕高清免费大全6| 亚洲av二区三区四区| 毛片女人毛片| 汤姆久久久久久久影院中文字幕| 亚洲精品国产av蜜桃| 成人国产av品久久久| 99热这里只有是精品在线观看| 五月伊人婷婷丁香| 国产色爽女视频免费观看| 色哟哟·www| 亚洲,欧美,日韩| 啦啦啦视频在线资源免费观看| 99精国产麻豆久久婷婷| 国产成人aa在线观看| 精品人妻视频免费看| 蜜桃久久精品国产亚洲av| 我的女老师完整版在线观看| av视频免费观看在线观看| 色视频在线一区二区三区| 亚洲成色77777| 大香蕉97超碰在线| 午夜日本视频在线| 内地一区二区视频在线| 2018国产大陆天天弄谢| 日本与韩国留学比较| 国产乱来视频区| 精华霜和精华液先用哪个| 尤物成人国产欧美一区二区三区| 国产男女内射视频| 日韩,欧美,国产一区二区三区| 精品久久久精品久久久| 日韩一本色道免费dvd| 久热久热在线精品观看| 亚洲精品乱久久久久久| 久久亚洲国产成人精品v| 91精品伊人久久大香线蕉| 国产免费福利视频在线观看| 天堂俺去俺来也www色官网| 久久久久精品久久久久真实原创| 精品一区二区三卡| 校园人妻丝袜中文字幕| 欧美日韩一区二区视频在线观看视频在线| 久久久久国产精品人妻一区二区| 久久久久精品久久久久真实原创| 男女无遮挡免费网站观看| 偷拍熟女少妇极品色| 免费av不卡在线播放| 亚洲欧美中文字幕日韩二区| 国产精品无大码| 精品亚洲成a人片在线观看 | 欧美最新免费一区二区三区| 精品一区二区三区视频在线| 亚洲精品日韩av片在线观看| 99热这里只有是精品50| 天堂8中文在线网| 波野结衣二区三区在线| 欧美xxxx性猛交bbbb| 久久精品国产自在天天线| 成人美女网站在线观看视频| 美女xxoo啪啪120秒动态图| 成人亚洲精品一区在线观看 | 一级毛片黄色毛片免费观看视频| 99热全是精品| 美女国产视频在线观看| 欧美高清成人免费视频www| 毛片一级片免费看久久久久| 人人妻人人爽人人添夜夜欢视频 | 亚洲精华国产精华液的使用体验| 99久久中文字幕三级久久日本| 极品少妇高潮喷水抽搐| 日日撸夜夜添| av在线播放精品| 成人美女网站在线观看视频| 亚洲av日韩在线播放| 少妇裸体淫交视频免费看高清| 五月开心婷婷网| av国产久精品久网站免费入址| 九九在线视频观看精品| 久久久久网色| 亚洲av国产av综合av卡| 晚上一个人看的免费电影| 草草在线视频免费看| 亚洲av综合色区一区| 亚洲精品亚洲一区二区| 日韩一区二区三区影片| 成年美女黄网站色视频大全免费 | 成人毛片a级毛片在线播放| 最近2019中文字幕mv第一页| 久久久午夜欧美精品| 亚洲av综合色区一区| 99热国产这里只有精品6| 老熟女久久久| 久久久久视频综合| 国产欧美日韩精品一区二区| 少妇丰满av| 美女视频免费永久观看网站| 最黄视频免费看| 老熟女久久久| 各种免费的搞黄视频| 精品一区二区三卡| 欧美另类一区| 欧美高清成人免费视频www| 久久人人爽av亚洲精品天堂 | 久久国产亚洲av麻豆专区| 网址你懂的国产日韩在线| 18禁裸乳无遮挡免费网站照片| 国产免费福利视频在线观看| 熟女人妻精品中文字幕| 一区二区三区精品91| 亚洲精品456在线播放app| 日韩不卡一区二区三区视频在线| av黄色大香蕉| 国产国拍精品亚洲av在线观看| 久久婷婷青草| 国产综合精华液| 麻豆成人av视频| 国产成人免费无遮挡视频| av天堂中文字幕网| 日韩欧美精品免费久久| 亚洲国产高清在线一区二区三| 中文乱码字字幕精品一区二区三区| 久久久久久九九精品二区国产| 国产乱来视频区| 精品久久久久久久久亚洲| 蜜桃久久精品国产亚洲av| 97超视频在线观看视频| 在线观看美女被高潮喷水网站| 青春草国产在线视频| 亚洲美女视频黄频| 日日啪夜夜爽| xxx大片免费视频| 老熟女久久久| 搡女人真爽免费视频火全软件| 亚洲精品一区蜜桃| 少妇的逼好多水| 欧美人与善性xxx| 成人无遮挡网站| 亚洲精品亚洲一区二区| 狠狠精品人妻久久久久久综合| 大码成人一级视频| 少妇的逼水好多| 亚洲av电影在线观看一区二区三区| 亚洲国产日韩一区二区| 大片免费播放器 马上看| 我的老师免费观看完整版| 婷婷色av中文字幕| 18禁裸乳无遮挡动漫免费视频| 免费看不卡的av| 欧美成人a在线观看| 国产女主播在线喷水免费视频网站| 亚洲欧美日韩另类电影网站 | 中国三级夫妇交换| 精品99又大又爽又粗少妇毛片| 日韩不卡一区二区三区视频在线| 高清在线视频一区二区三区| 国产男人的电影天堂91| 观看免费一级毛片| 精品久久国产蜜桃| 80岁老熟妇乱子伦牲交| 激情 狠狠 欧美| 黄色日韩在线| av女优亚洲男人天堂| 亚洲欧美精品专区久久| 免费黄网站久久成人精品| 久久久久视频综合| 亚洲精华国产精华液的使用体验| 亚洲,一卡二卡三卡| 波野结衣二区三区在线| 卡戴珊不雅视频在线播放| 色网站视频免费| 人妻少妇偷人精品九色| 国产女主播在线喷水免费视频网站|