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

    夏季煙臺(tái)—威海北部近海溶解氧濃度垂向分布的最小值

    2021-12-09 08:23:46孫利元劉子洲丁金強(qiáng)翟方國(guó)顧艷鎮(zhèn)
    海洋科學(xué) 2021年11期
    關(guān)鍵詞:耗氧量溶解氧表觀

    孫利元, 劉子洲, 丁金強(qiáng), 翟方國(guó), 顧艷鎮(zhèn)

    夏季煙臺(tái)—威海北部近海溶解氧濃度垂向分布的最小值

    孫利元1, 劉子洲2, 丁金強(qiáng)1, 翟方國(guó)2, 顧艷鎮(zhèn)2

    (1. 山東省漁業(yè)發(fā)展和資源養(yǎng)護(hù)總站, 山東 煙臺(tái) 264003; 2. 中國(guó)海洋大學(xué) 海洋與大氣學(xué)院, 山東 青島 266100)

    基于2020年夏季的大面航次觀測(cè)數(shù)據(jù), 分析了煙臺(tái)—威海北部海洋牧場(chǎng)及鄰近海域海水溶解氧濃度垂向分布最小值(氧最小值層)的空間分布特征, 并探討了影響因素。從6月至8月, 海水溶解氧濃度不斷減小, 垂向結(jié)構(gòu)亦存在顯著變化。海水溶解氧濃度垂向分布的最小值主要集中于7月的近岸海域, 最小值大致從外海向近岸方向減小, 其距離海底高度及與底層溶解氧濃度之差的絕對(duì)值均于雙島灣鄰近海域?yàn)樽畲?。海水溶解氧濃度垂向分布的最小值位于最?qiáng)密度層結(jié)以下。但是海水溶解氧濃度垂向分布最小值的強(qiáng)度向北減小, 而密度層結(jié)向北增大, 兩者的空間分布基本相反, 說明密度層結(jié)抑制垂向湍流擴(kuò)散可極大減少深層海水溶解氧的來源, 是海水溶解氧濃度垂向分布最小值形成的必要條件, 但不是主導(dǎo)因素。在海水溶解氧濃度垂向分布的最小值層, 表觀耗氧量存在垂向分布的最大值, 大部分站點(diǎn)的pH存在垂向分布的最小值, 說明局地增強(qiáng)、持續(xù)的生物地球化學(xué)耗氧是控制海水溶解氧濃度垂向分布最小值形成和空間分布的一個(gè)重要過程。研究結(jié)果表明氧最小值層是夏季煙臺(tái)—威海北部近岸海水溶解氧垂向結(jié)構(gòu)的典型特征之一。

    煙臺(tái)—威海北部近海; 海水溶解氧濃度; 垂向分布的最小值; 表觀耗氧量; pH

    隨著世界人口對(duì)水產(chǎn)品需求的不斷增長(zhǎng)和海洋主要天然漁場(chǎng)的衰退甚至消失[1], 越來越多的沿海國(guó)家開始將海洋牧場(chǎng)作為漁業(yè)資源增殖和海洋生態(tài)環(huán)境修復(fù)與保護(hù)的重要手段[2-4]。近年來, 我國(guó)海洋牧場(chǎng)發(fā)展迅速[5-6], 尤其是在山東省近海海域[7]。山東近海海洋牧場(chǎng)無論是發(fā)展速度還是總體數(shù)量均位于全國(guó)前列[7]。隨著近海海洋牧場(chǎng)數(shù)量的不斷增多、面積的不斷擴(kuò)大, 認(rèn)識(shí)海洋牧場(chǎng)及周邊海域生態(tài)環(huán)境的時(shí)空變化規(guī)律及其影響因素已成為海洋牧場(chǎng)高質(zhì)量發(fā)展的重要保障之一。

    煙臺(tái)—威海北部(圖1)是山東省海洋牧場(chǎng)分布最為密集的區(qū)域之一, 該海域生態(tài)環(huán)境的一個(gè)顯著特點(diǎn)是在夏季底層海水溶解氧濃度偏低, 個(gè)別年份甚至出現(xiàn)低氧[7-8]。低氧是指海水溶解氧濃度低于2.0 mg/L的狀態(tài)[9]。因?yàn)楹K芙庋鯘舛容^低, 低氧會(huì)導(dǎo)致生活在底層且活動(dòng)能力較弱的海洋生物大量死亡, 嚴(yán)重破壞海洋生態(tài)系統(tǒng), 從而引起海洋災(zāi)害[10]。煙臺(tái)—威海北部是山東省重要的貝、藻生態(tài)園區(qū), 其中主要經(jīng)濟(jì)品種行動(dòng)能力弱, 對(duì)海水溶解氧濃度的變化較為敏感。歷史上, 前人專門針對(duì)該海域海水溶解氧濃度時(shí)空變化特征及影響機(jī)制的研究較少。Zhang等[8]和Yang等[11]分別展示了2016年和2017年夏季海水溶解氧濃度的空間分布, 但并未分析影響機(jī)制。

    為了更好地認(rèn)識(shí)煙臺(tái)—威海北部近海水文-生態(tài)環(huán)境要素的變化動(dòng)態(tài), 理解夏季底層海水溶解氧虧損的時(shí)空分布及影響機(jī)理, 我們于2020年6月、7月和8月對(duì)該海域的水文、生態(tài)環(huán)境要素進(jìn)行了大面航次觀測(cè)。觀測(cè)發(fā)現(xiàn), 該海域海水溶解氧濃度垂向剖面中的最小值發(fā)生的深度存在顯著空間差異。在外海海域, 海水溶解氧濃度的最小值均發(fā)生在海洋表層或者底層。但是在近岸海域, 海水溶解氧濃度的垂向分布存在氧最小值層, 即溶解氧濃度最小值并非出現(xiàn)在表層或者底層, 而是出現(xiàn)在表層和底層之間的水體層, 與以往觀測(cè)不同[8]。本研究將基于2020年夏季的大面觀測(cè)數(shù)據(jù), 分析該海域海水溶解氧濃度垂向分布最小值(氧最小值層)的變化特征, 并探討可能影響因素。

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

    1.1 觀測(cè)數(shù)據(jù)

    在2020年夏季, 我們搭乘威海市瑜泰海洋牧場(chǎng)的工作船只分別于6月13號(hào)、7月15號(hào)和16號(hào)、8月17號(hào)和18號(hào)對(duì)煙臺(tái)—威海北部近海的水文、生態(tài)環(huán)境要素進(jìn)行了大面觀測(cè), 觀測(cè)站點(diǎn)分布如圖1所示。其中, 6月航次僅包含雙島灣外海的經(jīng)向斷面, 而7月和8月航次則包含了雙島灣外海的經(jīng)向斷面和雙島灣西側(cè)的沿岸斷面和金山港外海的經(jīng)向斷面。本文中, 雙島灣外海的經(jīng)向斷面稱為東斷面, 雙島灣西側(cè)的沿岸斷面與金山港外海的經(jīng)向斷面相連并統(tǒng)稱為西斷面。所有航次均使用日本亞力克公司生產(chǎn)的多參數(shù)水質(zhì)儀(型號(hào): AAQ171)對(duì)海水溫度、電導(dǎo)率、壓強(qiáng)、溶解氧濃度、pH和葉綠素(以下簡(jiǎn)稱葉綠素)濃度等的垂向剖面進(jìn)行觀測(cè), 儀器采樣頻率為4 Hz。本文中, 葉綠素和海水溶解氧濃度為質(zhì)量濃度。觀測(cè)深度和鹽度由儀器自帶軟件根據(jù)觀測(cè)的壓強(qiáng)、電導(dǎo)率、溫度等計(jì)算得到。

    圖1 2020年夏季煙臺(tái)—威海北部航次觀測(cè)站位分布圖

    1.2 方法

    本文根據(jù)Garcia等[12]得到的公式、利用觀測(cè)資料計(jì)算了海水的飽和溶解氧濃度和表觀耗氧量(apparent oxygen utilization, AOU)。表觀耗氧量等于飽和溶解氧濃度與實(shí)測(cè)海水溶解氧濃度之差, 量化了生物地球化學(xué)過程耗氧導(dǎo)致的海水溶解氧濃度的變化[13]。

    2 觀測(cè)結(jié)果

    2.1 海水溶解氧濃度的典型垂向剖面

    觀測(cè)結(jié)果顯示, 從6月至8月海水溶解氧濃度不斷減小, 其垂向結(jié)構(gòu)亦存在顯著變化。圖2展示了威海北部近海(37°34′N, 121°57′E)海水溶解氧濃度垂向剖面的時(shí)間變化。6月, 所有站點(diǎn)處海水溶解氧濃度的垂向結(jié)構(gòu)基本一致, 如圖2中紅色線所示, 在表層和底層之間存在濃度最大值, 而在表層或底層出現(xiàn)最小值。海水溶解氧濃度在最大值之上存在較為顯著的躍層。關(guān)于溶解氧濃度垂向分布中的最大值, 前人在黃海海域已有較多觀測(cè), 并認(rèn)為是夏季海水溶解氧濃度垂向結(jié)構(gòu)中一個(gè)較為常見的現(xiàn)象[14-17]。7月, 部分站點(diǎn)海水溶解氧濃度存在氧最小值層, 即溶解氧濃度最小值并非出現(xiàn)在表層或者底層, 而是出現(xiàn)在表層和底層之間的水體層, 如圖2中綠色線所示。在該站點(diǎn), 海水溶解氧濃度最小值約為6.4 mg/L, 其深度約為11.7 m, 距離海底高度約為10.0 m。在8月份,海水溶解氧濃度的分層現(xiàn)象更為明顯, 如圖2中藍(lán)色線所示。海水溶解氧濃度在深層隨深度增大而迅速減小, 并于底層達(dá)到最小值。根據(jù)上述分析可以看出, 觀測(cè)期間海水溶解氧濃度僅在7月份存在氧最小值層。

    2.2 海水溶解氧濃度垂向分布最小值的空間分布

    圖3展示了2020年7月海水溶解氧濃度沿東斷面和西斷面的空間分布。沿東斷面, 海水溶解氧濃度約為5.29~8.54 mg/L; 而沿西斷面, 海水溶解氧濃度約為4.81~9.04 mg/L。在兩個(gè)斷面的表層, 海水溶解氧濃度大致從外海向近岸區(qū)域減小, 與以往觀測(cè)一致[8]。而在兩個(gè)斷面的底層, 海水溶解氧濃度的空間分布較為復(fù)雜, 于外海和37°35′N處存在極大值。底層海水溶解氧濃度的空間分布與以往觀測(cè)稍有不同, 可能存在年際等不同時(shí)間尺度的變化[18]。海水溶解氧濃度垂向結(jié)構(gòu)的最顯著特征是在外海區(qū)域存在氧最大值層, 而在近岸區(qū)域存在氧最小值層。觀測(cè)期間海水溶解氧濃度的最大值和最小值亦均出現(xiàn)在表層和底層之間的深度層, 其中最大值出現(xiàn)在西斷面外海站位約6 m深度處, 最小值出現(xiàn)在離岸最近站位約12 m深度處。

    圖2 2020年夏季威海北部近海觀測(cè)站點(diǎn)(37°34′N, 121°57′E)海水溶解氧濃度垂向剖面的時(shí)間變化

    沿東、西兩斷面, 海水溶解氧濃度垂向分布最小值的空間變化稍有不同。在東斷面, 海水溶解氧濃度垂向分布的最小值在空間上較為集中, 存在于37°37′N以南的近岸站位。而在西斷面, 海水溶解氧濃度垂向分布的最小值在空間上較為分散, 存在于近岸的4個(gè)站位、37°35′N附近的3個(gè)站位和外海的1個(gè)站位。圖3(c)和圖3(d)分別比較了海水溶解氧濃度垂向分布最小值及其與底層海水溶解氧濃度之差沿東斷面和西斷面的空間變化。在夏季, 由于密度層結(jié)和生物地球化學(xué)耗氧等過程的影響, 深層海水溶解氧濃度一般隨深度增大而減小, 如圖2所示。因此, 可以使用底層海水溶解氧濃度減去海水溶解氧濃度最小值的差表征氧最小值層的強(qiáng)度。在東斷面, 海水溶解氧濃度垂向分布的最小值大致向近岸方向單調(diào)減小, 而其強(qiáng)度大致向近岸增大, 強(qiáng)度最大值(0.91 mg/L)出現(xiàn)在從岸邊至外海第2個(gè)觀測(cè)站點(diǎn)。在西斷面, 海水溶解氧濃度垂向分布最小值在37°37′N以南亦向近岸方向減小, 其強(qiáng)度的空間變化較為復(fù)雜, 強(qiáng)度最大值(0.47 mg/L)出現(xiàn)在37°34′N附近。

    圖3 2020年7月海水溶解氧濃度觀測(cè)

    注: a和b中, 黑色三角形表示觀測(cè)站點(diǎn), 紅色星號(hào)標(biāo)注海水溶解氧濃度垂向分布的最小值

    總體而言, 在2020年7月觀測(cè)期間, 海水溶解氧濃度垂向分布的最小值主要集中于近岸海域, 且最小值從外海向近岸方向減小。底層海水溶解氧濃度與溶解氧濃度最小值的差于雙島灣鄰近海域最大。

    3 影響因素分析與討論

    海水溶解氧濃度的空間分布和時(shí)間變化受各種海洋動(dòng)力過程和生物地球化學(xué)過程的共同影響[13, 19]。動(dòng)力因素和過程主要包括海水層結(jié)、垂向湍流擴(kuò)散、潮流或者海流引起的平流輸運(yùn)等, 生物地球化學(xué)過程則主要包括浮游植物等光合作用生成氧氣、魚類等呼吸作用消耗氧氣、有機(jī)物分解消耗氧氣等[19]。對(duì)于海水溶解氧濃度在不同區(qū)域、不同時(shí)間尺度上的變化, 主要影響因素和過程亦不相同。下面將基于觀測(cè)資料分析討論海水溶解氧濃度垂向分布最小值形成和空間變化的影響因素。

    3.1 各觀測(cè)要素的垂向結(jié)構(gòu)

    首先以位于威海北部近海(37°34′N, 121°57′E)的觀測(cè)站點(diǎn)為例, 分析2020年夏季各觀測(cè)要素垂向結(jié)構(gòu)的時(shí)間變化。如圖4(a)和(b)所示, 海水溫度和鹽度呈現(xiàn)相反的時(shí)間變化和垂向結(jié)構(gòu)。從6月至8月, 溫度升高、鹽度降低; 從表層至底層, 溫度減小、鹽度增大。溫度和鹽度的垂向結(jié)構(gòu)存在較為顯著的躍層, 均于7月最強(qiáng), 最大值和最小值均出現(xiàn)在表層或者底層。受溫度和鹽度的影響, 海水密度從6月至8月隨時(shí)間減小、從表層至底層隨深度增大, 垂向結(jié)構(gòu)存在較為顯著的密度躍層, 并于7月最強(qiáng)[圖4(c)]。溫度和鹽度可以影響海水的飽和溶解氧濃度。根據(jù)觀測(cè)資料計(jì)算表明, 海水密度和飽和溶解氧濃度的時(shí)空變化主要由溫度控制[18]。受溫度的影響, 海水飽和溶解氧濃度從6月至8月隨時(shí)間減小、從表層至底層隨深度增大, 最小值出現(xiàn)在表層, 與7月溶解氧濃度的垂向分布不一致。

    相比于溫度和鹽度, 葉綠素濃度的垂向結(jié)構(gòu)存在顯著的時(shí)間變化, 在6月隨深度增大, 而在7月和8月則在水體內(nèi)部存在最大值層[圖4(d)]。在7月的觀測(cè)中, 葉綠素濃度垂向分布最大值的深度和溶解氧濃度垂向分布最小值的深度較為接近。

    時(shí)間變化方面, 表觀耗氧量與觀測(cè)溶解氧濃度和飽和溶解氧濃度呈完全相反的趨勢(shì)。從6月至8月, 表觀耗氧量隨時(shí)間增大, 且在底層增大速度最快。在6月, 表觀耗氧量約為0, 隨深度增大略有增加。而在8月, 整層水體的表觀耗氧量均大于0, 且隨深度增大而顯著增加, 最大可達(dá)3.20 mg/L。該結(jié)果說明從6月開始, 生物地球化學(xué)耗氧過程對(duì)觀測(cè)海域深層海水溶解氧濃度變化的影響顯著增強(qiáng)。在7月, 表觀耗氧量的垂向結(jié)構(gòu)與觀測(cè)溶解氧濃度的垂向結(jié)構(gòu)相反, 在深度約11.7 m處存在垂向分布的最大值(1.37 mg/L)。

    海水pH呈現(xiàn)與海水溶解氧濃度相似的時(shí)間變化和垂向結(jié)構(gòu), 從6月至8月隨時(shí)間減小、且減小速度在底層最大。在6月, pH約為8.47~8.48, 垂向變化較小。在8月, 表層pH減小至8.41, 而底層pH減小至8.18。深層海水pH的減小間接說明生物地球化學(xué)耗氧過程對(duì)深層海水溶解氧濃度的變化具有重要影響。在7月, pH的垂向結(jié)構(gòu)與溶解氧濃度的垂向結(jié)構(gòu)類似, 在深度約12.3 m處存在最小值(8.34)。

    上述分析說明, 觀測(cè)站點(diǎn)處海水密度層結(jié)和生物地球化學(xué)耗氧過程對(duì)海水溶解氧濃度垂向分布最小值的形成有重要影響。下面將繼續(xù)結(jié)合7月各要素的斷面分布探討海水溶解氧濃度垂向分布最小值形成和空間變化的影響因素。

    3.2 影響要素的斷面分布

    圖5(a)和(b)展示了2020年7月浮力頻率平方沿東、西斷面的分布。沿東斷面從近岸向外海方向, 密度躍層加深、強(qiáng)度增大。而沿西斷面從近岸向外海方向, 密度躍層先加深后變淺、強(qiáng)度先增大后減小。海水溶解氧濃度垂向分布的最小值均出現(xiàn)在浮力頻率平方最大值、即最強(qiáng)密度層結(jié)以下。這說明夏季密度層結(jié)抑制溶解氧的垂向湍流擴(kuò)散是影響海水溶解氧濃度垂向分布最小值形成的重要因素。但是在近岸海域, 密度層結(jié)從外海向近岸方向減弱, 與海水溶解氧濃度垂向分布最小值的空間變化基本相反, 與后者強(qiáng)度的空間變化亦無顯著相關(guān)關(guān)系(圖3)。因此, 密度層結(jié)抑制垂向湍流擴(kuò)散可極大減少深層海水溶解氧的來源, 是海水溶解氧濃度垂向分布最小值發(fā)生的必要條件, 但不是主導(dǎo)因素。

    圖4 2020年夏季煙臺(tái)—威海北部近海觀測(cè)站點(diǎn)(37°34′N, 121°57′E)各觀測(cè)要素垂向剖面的時(shí)間變化

    圖5(c)和(d)展示了2020年7月葉綠素濃度沿東、西斷面的空間分布。水平上, 海水葉綠素濃度大致從近岸向外海減小, 與衛(wèi)星觀測(cè)結(jié)果一致[20]。垂向上, 各站點(diǎn)處葉綠素濃度的垂向結(jié)構(gòu)均存在最大值層, 與夏季黃海其他海域的觀測(cè)結(jié)果一致[17, 21-22]。海水溶解氧濃度垂向分布最小值的深度基本位于葉綠素濃度垂向分布最大值的深度之下, 且兩者較為接近。海水葉綠素濃度越高, 浮游植物光合作用生成的氧氣就越多, 可導(dǎo)致海水溶解氧濃度越高。因此, 觀測(cè)結(jié)果表明海水溶解氧濃度垂向分布最小值的形成和空間變化并不是由浮游植物等光合作用生成氧氣過程的減弱引起的, 與葉綠素濃度的時(shí)空變化無關(guān)。

    為了了解生物地球化學(xué)耗氧過程的影響, 本文首先分析了2020年7月觀測(cè)期間海水表觀耗氧量的空間分布。如圖5(e)和(f)所示, 在外海海域上層, 海水表觀耗氧量為負(fù)值, 說明海水溶解氧處于過飽和狀態(tài); 而在外海海域深層和近岸海域, 海水表觀耗氧量為正值, 說明生物地球化學(xué)耗氧過程強(qiáng)烈、海水溶解氧存在顯著虧損??傮w而言, 海水表觀耗氧量的空間分布和溶解氧濃度的空間分布呈顯著負(fù)相關(guān), 兩者沿東、西斷面的空間相關(guān)系數(shù)分別為–0.92、–0.97。對(duì)應(yīng)于海水溶解氧濃度垂向分布的最小值, 表觀耗氧量垂向分布存在最大值。此外, 表觀耗氧量的最小值(–1.58 mg/L)和最大值(2.82 mg/L)分別出現(xiàn)在西斷面的外海觀測(cè)站位和離岸最近的觀測(cè)站位。根據(jù)表觀耗氧量的空間分布, 初步推測(cè)夏季近海深層生物地球化學(xué)耗氧過程的增強(qiáng)和持續(xù)是該海域深層海水溶解氧顯著虧損和溶解氧濃度垂向分布最小值形成的重要原因。

    圖5 2020年7月海水浮力頻率平方、葉綠素濃度、表觀耗氧量和pH的空間分布

    注: 黑色三角形為觀測(cè)站點(diǎn), 星號(hào)為海水溶解氧濃度垂向分布最小值的深度; e和f中, 黑色實(shí)線為表觀耗氧量為0的等值線

    為了進(jìn)一步證實(shí)上述觀點(diǎn), 圖5(g)和(h)展示了海水pH沿東、西斷面的空間分布。深層海水的pH明顯小于上層海水的pH, 且在近岸海域向岸邊減小。海水表觀耗氧量的空間分布與pH的空間分布呈顯著負(fù)相關(guān), 兩者沿東、西斷面的空間相關(guān)系數(shù)分別為–0.96、–0.97。表觀耗氧量與pH空間分布的負(fù)相關(guān)關(guān)系在深層尤為顯著, 如圖6(a)和(b)所示。基于最小二乘法可得兩者之間的線性關(guān)系式: 東斷面為AOU= –18.5·pH+155.8, 西斷面為AOU= –17.4·pH+ 146.2, 其中AOU和pH分別為海水的表觀耗氧量和pH。

    圖6 2020年7月密度躍層以下海水表觀耗氧量與pH的比較[(a)和(b)]; 溶解氧濃度垂向分布最小值對(duì)應(yīng)的pH及其與底層pH之差的空間變化[(c)和(d)]

    圖6(c)和(d)比較了海水溶解氧濃度垂向分布最小值對(duì)應(yīng)的pH及其與底層海水pH之差沿東、西斷面的空間分布。利用底層海水pH減去海水溶解氧濃度最小值深度處的pH表征氧最小值層pH的強(qiáng)度??梢钥闯? 在氧最小值層中, pH及其強(qiáng)度與溶解氧質(zhì)量濃度及其強(qiáng)度的空間分布基本一致。在海洋中, pH的變化主要受CO2含量的影響。海洋-大氣界面的氣體交換和海洋生物呼吸、有機(jī)物分解等過程是海洋CO2的主要來源[23]。在觀測(cè)期間, 研究海域存在顯著的密度層結(jié), 可同步抑制溶解氧和CO2的垂向湍流擴(kuò)散, 減弱海洋-大氣界面氣體交換對(duì)深層溶解氧和CO2含量的影響。因此對(duì)于深層水體而言, 表觀耗氧量與pH之間顯著的負(fù)相關(guān)關(guān)系進(jìn)一步證明局地生物呼吸和有機(jī)物分解等生物地球化學(xué)過程是影響海水溶解氧和CO2含量同步變化的重要過程。

    根據(jù)上述結(jié)果, 溶解氧濃度垂向分布最小值的形成機(jī)制可能為: 一方面, 觀測(cè)海域存在密度層結(jié), 可極大抑制溶解氧的垂向湍流擴(kuò)散, 從而減少密度躍層以下海水溶解氧的來源; 另一方面, 密度躍層以下水體中生物地球化學(xué)過程持續(xù)耗氧, 且耗氧過程在海底之上的水體層最強(qiáng), 從而形成溶解氧濃度垂向分布的最小值。

    此外, 在部分海水溶解氧濃度存在最小值層的站位, pH垂向分布中的最小值出現(xiàn)在底層。這說明海水溶解氧濃度垂向分布最小值的形成還受其他過程或者因素的影響。在觀測(cè)海域, 底層海水溶解氧濃度和pH從外海向近岸方向減小。Zhai等[18]根據(jù)定點(diǎn)連續(xù)觀測(cè)數(shù)據(jù)指出, 該海域流場(chǎng)主要受海面風(fēng)的控制[24]。在夏季, 離岸南風(fēng)可引起底層向岸流。因此, 底層流的向岸平流輸運(yùn)也可能導(dǎo)致近岸海域海水溶解氧濃度的最小值出現(xiàn)在海底之上的水體層而非底層[18]。由于研究區(qū)域岸線復(fù)雜、缺少高分辨率的流場(chǎng)觀測(cè)數(shù)據(jù)或數(shù)值模擬結(jié)果, 本文無法量化底層流對(duì)海水溶解氧濃度最小值層的形成和時(shí)空變化的貢獻(xiàn), 未來還需要進(jìn)一步研究。

    4 結(jié)論與展望

    作為海洋生態(tài)環(huán)境的重要參數(shù), 海水溶解氧濃度的時(shí)空變化對(duì)海洋牧場(chǎng)的可持續(xù)健康發(fā)展極為重要?;?020年夏季的大面航次觀測(cè)數(shù)據(jù), 本文首次分析了煙臺(tái)—威海北部海洋牧場(chǎng)及鄰近海域海水溶解氧濃度垂向分布最小值(即氧最小值層)的空間變化特征, 并探討了影響因素。主要結(jié)論如下:

    1) 從6月至8月, 觀測(cè)海域海水溶解氧濃度不斷減小, 垂向結(jié)構(gòu)亦存在顯著變化。海水溶解氧濃度垂向分布的最小值僅出現(xiàn)在7月。沿觀測(cè)斷面, 海水溶解氧濃度垂向分布的最小值主要出現(xiàn)在近岸海域, 最小值大致從外海向近岸方向減小, 其距離海底高度和與底層溶解氧濃度之差的絕對(duì)值均于雙島灣鄰近海域?yàn)樽畲蟆?/p>

    2) 海水溶解氧濃度垂向分布的最小值位于最強(qiáng)密度層結(jié)以下, 但兩者的空間分布基本相反, 說明密度層結(jié)抑制垂向湍流擴(kuò)散可極大減少深層海水溶解氧的來源, 是海水溶解氧濃度垂向分布最小值形成的必要條件, 但不是主導(dǎo)因素。

    3) 葉綠素濃度的垂向結(jié)構(gòu)在表層和底層之間存在最大值層, 其深度與海水溶解氧濃度最小值的深度相近, 表明葉綠素濃度的時(shí)空變化并不是海水溶解氧濃度垂向分布最小值的形成原因。

    4) 從6月至8月, 密度層結(jié)以下海水表觀耗氧量為正值且顯著增大, pH顯著減小。在大部分站點(diǎn)處, 對(duì)應(yīng)于溶解氧濃度垂向分布的最小值, 表觀耗氧量存在垂向分布的最大值, pH存在垂向分布的最小值, 說明局地增強(qiáng)、持續(xù)的生物地球化學(xué)耗氧過程是控制海水溶解氧濃度垂向分布最小值形成和空間變化的一個(gè)重要過程。

    煙臺(tái)—威海北部近海海水溶解氧濃度的最小值層與大洋中常見的次表層氧最小值層(oxygen minimum zone)[25-26]在時(shí)空變化方面不同。大洋中的次表層氧最小值層存在于太平洋、大西洋和印度洋, 且常年存在。但是兩者在形成機(jī)制方面是類似的, 其中一個(gè)重要機(jī)制是在密度層結(jié)抑制溶解氧垂向擴(kuò)散的同時(shí), 生物地球化學(xué)過程持續(xù)耗氧導(dǎo)致溶解氧濃度不斷降低。未來工作將繼續(xù)結(jié)合高分辨率觀測(cè)或者水動(dòng)力-生物地球化學(xué)耦合數(shù)值模擬結(jié)果研究該海域生物地球化學(xué)耗氧過程的時(shí)空變化和影響機(jī)制, 同時(shí)利用溶解氧的收支方程探討水動(dòng)力導(dǎo)致的平流輸運(yùn)等不同過程對(duì)海水溶解氧濃度最小值層形成和時(shí)空變化的貢獻(xiàn)。

    [1] 劉子飛. 中國(guó)近海漁場(chǎng)荒漠化: 評(píng)價(jià)、原因與治理[J]. 農(nóng)業(yè)經(jīng)濟(jì)問題, 2019(6): 105-116.

    LIU Zifei. Research on the status, causes and governance of marine fishing ground desertification[J]. Issue in Agricultural Economy, 2019(6): 105-116.

    [2] MUSTAFA S. Stock enhancement and sea ranching: objectives and potential[J]. Reviews in Fish Biology and Fisheries, 2003, 13: 141-149.

    [3] BELL J. A new era for restocking, stock enhancement and sea ranching of coastal fisheries resources[J]. Reviews in Fisheries Science, 2008, 16: 1-9.

    [4] 陳丕茂, 舒黎明, 袁華榮, 等. 國(guó)內(nèi)外海洋牧場(chǎng)發(fā)展歷程與定義分類概述[J]. 水產(chǎn)學(xué)報(bào), 2019, 43(9): 1851- 1869.

    CHEN Pimao, SHU Liming, YUAN Huarong, et al. Review on development, definition and classification of marine ranching in domestic and overseas[J]. Journal of Fisheries of China, 2019, 43(9): 1851-1869.

    [5] 楊紅生. 我國(guó)海洋牧場(chǎng)建設(shè)回顧與展望[J]. 水產(chǎn)學(xué)報(bào), 2016, 40(7): 1133-1140.

    YANG Hongsheng. Construction of marine ranching in China: reviews and prospects[J]. Journal of Fisheries of China, 2016, 40(7): 1133-1140.

    [6] YU Jinkai, ZHANG Linlin. Evolution of marine ranching policies in China: Review, performance and prospects[J]. Science of the Total Environment, 2020, 737: 139782.

    [7] 翟方國(guó), 顧艷鎮(zhèn), 李培良, 等. 山東省海洋牧場(chǎng)觀測(cè)網(wǎng)的建設(shè)與發(fā)展[J]. 海洋科學(xué), 2020, 44(12): 93-106.

    ZHAI Fangguo, GU Yanzhen, LI Peiliang, et al. Construction and development of marine ranch observation network in Shandong Province[J]. Marine Sciences, 2020, 44(12): 93-106.

    [8] ZHANG Yong, GAO Xuelu, GUO Weidong, et al. Origin and dynamics of dissolved organic matter in mari-cultured area suffering from summertime hypoxia and acidification[J]. Frontiers in Marine Science, 2018, 5(325): 1-14.

    [9] DIAZ R. Overview of hypoxia around the world[J]. Journal of Environmental Quality, 30(2): 275-281.

    [10] BREITBURG D, HONDORP D, DAVIAS L et al. Hypoxia, nitrogen, and fisheries: Integrating effects across local and global landscapes[J]. Annual Review of Marine Science, 2009, 1: 329-349.

    [11] YANG Bo, GAO Xuelu. Chromophoric dissolved organic matter in summer in a coastal mariculture region of northern Shandong Peninsula, North Yellow Sea[J]. Continental Shelf Research, 2019, 176: 19-35.

    [12] GARCIA H, Gordon I. Oxygen solubility in seawater: Better fitting equations[J]. Limnology and Oceanology, 1992, 37(6): 1307-1312.

    [13] FENNEL K, TESTA J. Biogeochemical controls on coastal hypoxia[J]. Annual Review of Marine Science, 2019, 11: 105-130.

    [14] 顧宏堪. 海水溶解氧夏季垂直分布中的最大值[J]. 海洋與湖沼, 1966, 8(2): 85-91.

    GU Hongkan. On the maximum value of dissolved oxygen in its vertical distribution in the sea[J]. Oceano-logia et Limnologia Sinica, 1966, 8(2): 85-91.

    [15] 藍(lán)淑芳. 南黃海斜斷面溶解氧分布特征及其微細(xì)結(jié)構(gòu)[J]. 海洋科學(xué), 1988, 12(3): 1-5.

    LAN Shufang. Distribution characteristics and microstructure of dissolved oxygen along an oblique section in the south Huanghai Sea[J]. Marine Sciences, 1988, 12(3): 1-5.

    [16] 王保棟, 王桂云, 鄭昌洙, 等. 南黃海溶解氧的垂直分布特性[J]. 海洋學(xué)報(bào), 1999, 21(5): 72-77.

    WANG Baodong, WANG Guiyun, ZHENG Changzhu, et al. Features in vertical profiles of dissolved oxygen in the southern Huanghai Sea[J]. Acta Oceanologica Sinica, 1999, 21(5): 72-77.

    [17] 韋欽勝, 傅明珠, 李艷, 等. 南黃海冷水團(tuán)海域溶解氧和葉綠素最大值現(xiàn)象及營(yíng)養(yǎng)鹽累積的季節(jié)演變[J]. 海洋學(xué)報(bào), 2013, 35(4): 142-154.

    WEI Qinsheng, FU Mingzhu, LI Yan, et al. Observation of the seasonal evolution of DO, Chlorophyll a maximum phenomena and nutrient accumulating in the southern Huanghai (Yellow) Sea Cold Water Mass area[J]. Acta Oceanologica Sinica, 2013, 35(4): 142-154.

    [18] ZHAI Fangguo, LIU Zizhou, LI Peiliang, et al. Physical controls of summer variations in bottom layer oxygen concentrations in the coastal hypoxic region off the northeastern Shandong Peninsula in the Yellow Sea[J]. Journal of Geophysical Research Oceans, 2021, 126(5): e2021JC017299.

    [19] TESTA J, KEMP W. Oxygen – Dynamics and biogeochemical consequences[M]//WOLANSKI E, MCLUSKY D. Treatise on estuarine and coastal science. Waltham: Academic Press, 2011: 163-199.

    [20] SONG Xiukai, GU Yanzhen, ZHAI Fangguo, et al. Climatology and seasonal variability of satellite-derived chlorophyllaround the Shandong Peninsula[J]. Journal of Oceanology and Limnology, 2021, 39(4): 1222-1244.

    [21] 傅明珠, 王宗靈, 孫萍, 等. 2006年夏季南黃海浮游植物葉綠素a分布特征及其環(huán)境調(diào)控機(jī)制[J]. 生態(tài)學(xué)報(bào), 2009, 29(10): 5366-5375.

    FU Mingzhu, WANG Zongling, SUN Ping, et al. Spatial distribution characteristics and the environmental regulation mechanisms of phytoplankton chlorophyll a in the Southern Yellow Sea during summer 2006[J]. Acta Ecologica Sinica, 2009, 29(10): 5366-5375.

    [22] FU Mingzhu, SUN Ping, WANG Zongling, et al. Structure, characteristics and possible formation me-chanisms of the subsurface chlorophyll maximum in the Yellow Sea Cold Water Mass[J]. Continental Shelf Research, 2018, 165: 93-105.

    [23] CAI Weijun, FEELY R, TESTA J, et al. Natural and anthropogenic drivers of acidification in large estuaries[J]. Annual Review of Marine Science, 2021, 13: 23-55.

    [24] 胡羅煜, 翟方國(guó), 劉子洲, 等. 山東半島東北部海洋牧場(chǎng)海域冬季海水流動(dòng)的時(shí)空分布特征[J]. 海洋科學(xué), 2021, 45(4): 1-12.

    HU Luoyu, ZHAI Fangguo, LIU Zizhou, et al. Spatial and temporal characteristics of sea current within the marine ranches in the northeast of Shandong Peninsula[J]. Marine Sciences, 2021, 45(4): 1-12.

    [25] KARSTENSEN J, STRAMMA L, VISBECK M. Oxygen minimum zones in the eastern tropical Atlantic and Pacific Oceans[J]. Progress in Oceanography, 2008, 77: 331-350.

    [26] PAULMIER A, RUIZ-PINO D. Oxygen minimum zones (OMZs) in the modern ocean[J]. Progress in Oceanography, 2009, 80: 113-128.

    Minima in vertical distributions of dissolved oxygen concentration in the northern coastal oceans of Yantai-Weihai

    SUN Li-yuan1, LIU Zi-zhou2, DING Jin-qiang1, ZHAI Fang-guo2, GU Yan-zhen2

    (1. Shandong Fisheries Development and Resources Conservation Center, Yantai 264003, China; 2. College of Oceanic and Atmospheric Sciences, Ocean University of China, Qingdao 266100, China)

    The present study investigates the spatial distribution of minima of dissolved oxygen (DO) concentration in vertical distributions (oxygen minimum layer) and the influencing factors in the northern coastal oceans of Yantai-Weihai based on cruise observations recorded in 2020’s summer. From June to August, the DO concentration decreased, and its vertical distribution changed significantly. In July, DO concentration minima in the vertical distributions were observed primarily in nearshore regions and decreased shoreward; its height above the sea bottom and difference from the bottom layer DO concentration were the largest in Shuangdao Bay. Besides, DO concentration minima in the vertical distributions were located below the strongest density stratification. However, the strength of DO concentration minima decreased northward, and this was opposite to density stratification, which increased northward. This finding indicated that density stratification could significantly restrict oxygen from entering the bottom layer of water by inhibiting turbulence-induced vertical diffusion. Thus, it was necessary but not a dominant factor for the formation of the oxygen minimum layer. Apparent oxygen utilization had maxima and pH had minima in the DO minimum layers at most observation sites. This high degree of correlation suggested that locally enhanced and continuous biogeochemical oxygen consumption were important in controlling the formation and spatial distribution of the oxygen minimum layer. The present study indicated that the oxygen minimum layer is one of the prominent features in vertical distributions of DO in nearshore waters of northern Yantai-Weihai.

    northern coastal oceans of Yantai-Weihai; dissolved oxygen concentration; minimum in vertical distribution; apparent oxygen utilization; pH

    Apr. 12, 2021

    P715

    A

    1000-3096(2021)11-0020-10

    10.11759/hykx20210412002

    2021-04-12;

    2021-05-25

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

    [Shandong Provincial Natural Science Foundation, No. ZR2020MD059; National Science Foundation of China, No. 42176016; Key Research and Development Project of Zhejiang Province, No. 2020C03012; Key Research and Development Project of Guangdong Province, No. 2020B1111030002; Major Science and Technology Project of Sanya YZBSTC, No. SKJC-KJ-2019KY03; National Key Research and Development Program of China, No. 2019YFD0901305]

    孫利元(1980—), 男, 山東招遠(yuǎn)人, 高級(jí)工程師, 主要研究方向?yàn)楹Q竽翀?chǎng)、漁業(yè)資源增殖, E-mail: heroland80@ 163.com; 翟方國(guó)(1983—),通信作者, 男, 山東濰坊人, 副教授, 主要從事海洋動(dòng)力、生態(tài)與氣候變化的研究, E-mail: gfzhai@ ouc.edu.cn

    (本文編輯: 叢培秀)

    猜你喜歡
    耗氧量溶解氧表觀
    綠盲蝽為害與赤霞珠葡萄防御互作中的表觀響應(yīng)
    河北果樹(2021年4期)2021-12-02 01:14:50
    淺析水中溶解氧的測(cè)定
    鋼結(jié)構(gòu)表觀裂紋監(jiān)測(cè)技術(shù)對(duì)比與展望
    上海公路(2019年3期)2019-11-25 07:39:28
    例析對(duì)高中表觀遺傳學(xué)的認(rèn)識(shí)
    影響飲用水水質(zhì)耗氧量檢測(cè)的因素探討
    尤溪河某河段表觀耗氧量與營(yíng)養(yǎng)鹽相關(guān)性研究*
    污水活性污泥處理過程的溶解氧增益調(diào)度控制
    城市河道洲灘對(duì)水流溶解氧分布的影響
    水中耗氧量的準(zhǔn)確測(cè)定影響因素分析
    表觀遺傳修飾在糖脂代謝中的作用
    遺傳(2014年3期)2014-02-28 20:58:52
    9热在线视频观看99| 国产亚洲一区二区精品| 人人妻人人澡人人看| 母亲3免费完整高清在线观看| 国产精品二区激情视频| 国产高清激情床上av| 99国产极品粉嫩在线观看| 欧美日韩av久久| 夜夜骑夜夜射夜夜干| 怎么达到女性高潮| 啦啦啦中文免费视频观看日本| av天堂久久9| 一级片'在线观看视频| 欧美黄色淫秽网站| 最黄视频免费看| 久久久久国内视频| 啪啪无遮挡十八禁网站| 日本a在线网址| 亚洲性夜色夜夜综合| 人妻一区二区av| 欧美激情 高清一区二区三区| 日韩视频在线欧美| 男女无遮挡免费网站观看| 国产高清激情床上av| 日韩精品免费视频一区二区三区| 色在线成人网| 男女午夜视频在线观看| 99久久人妻综合| 人人妻,人人澡人人爽秒播| 嫩草影视91久久| 亚洲少妇的诱惑av| 丁香六月天网| 大型黄色视频在线免费观看| 亚洲国产欧美日韩在线播放| 丰满饥渴人妻一区二区三| 天堂俺去俺来也www色官网| 手机成人av网站| 99久久99久久久精品蜜桃| 99国产精品99久久久久| 一区二区三区乱码不卡18| 午夜日韩欧美国产| 国产一区二区 视频在线| 免费观看a级毛片全部| 丝袜人妻中文字幕| 国产在视频线精品| 超色免费av| 十八禁人妻一区二区| 欧美日韩一级在线毛片| 无人区码免费观看不卡 | 十八禁人妻一区二区| 精品国产乱码久久久久久小说| 久久久国产成人免费| 免费人妻精品一区二区三区视频| 一区在线观看完整版| 女警被强在线播放| 亚洲第一欧美日韩一区二区三区 | 亚洲国产欧美网| 人人妻人人添人人爽欧美一区卜| 黄色怎么调成土黄色| 99国产精品99久久久久| 久久精品亚洲精品国产色婷小说| 美国免费a级毛片| 在线观看舔阴道视频| 757午夜福利合集在线观看| 老熟妇乱子伦视频在线观看| 777久久人妻少妇嫩草av网站| 下体分泌物呈黄色| 国产精品亚洲一级av第二区| 精品一区二区三区视频在线观看免费 | bbb黄色大片| 久久久国产成人免费| 国产精品秋霞免费鲁丝片| 精品久久久久久电影网| 人妻久久中文字幕网| 91麻豆av在线| 淫妇啪啪啪对白视频| 国产成人免费无遮挡视频| 热re99久久国产66热| netflix在线观看网站| 日本av手机在线免费观看| 中文字幕av电影在线播放| 亚洲精品av麻豆狂野| 一区二区三区国产精品乱码| 91精品国产国语对白视频| 亚洲中文日韩欧美视频| 免费日韩欧美在线观看| 亚洲精品一卡2卡三卡4卡5卡| 麻豆成人av在线观看| 国产无遮挡羞羞视频在线观看| 久久婷婷成人综合色麻豆| 日韩欧美国产一区二区入口| 精品亚洲成a人片在线观看| 国产又爽黄色视频| 男女无遮挡免费网站观看| 欧美激情 高清一区二区三区| 欧美日本中文国产一区发布| 黄片小视频在线播放| 亚洲男人天堂网一区| 美女高潮喷水抽搐中文字幕| 亚洲精品国产精品久久久不卡| 99re6热这里在线精品视频| 捣出白浆h1v1| 黑丝袜美女国产一区| 国产av精品麻豆| 亚洲黑人精品在线| 一区二区三区精品91| 国产欧美日韩综合在线一区二区| 可以免费在线观看a视频的电影网站| 婷婷丁香在线五月| 亚洲第一欧美日韩一区二区三区 | 老司机亚洲免费影院| 亚洲成人免费电影在线观看| 黄色 视频免费看| 精品亚洲成国产av| 午夜福利一区二区在线看| 99国产精品免费福利视频| 亚洲av日韩精品久久久久久密| 性色av乱码一区二区三区2| 麻豆av在线久日| 久久这里只有精品19| 成年女人毛片免费观看观看9 | 高清毛片免费观看视频网站 | 热99国产精品久久久久久7| 老熟妇乱子伦视频在线观看| 成人特级黄色片久久久久久久 | 午夜日韩欧美国产| 少妇的丰满在线观看| 国产精品 欧美亚洲| 制服人妻中文乱码| 黑人欧美特级aaaaaa片| 日日摸夜夜添夜夜添小说| 搡老熟女国产l中国老女人| 下体分泌物呈黄色| 亚洲专区中文字幕在线| 久久毛片免费看一区二区三区| 久热这里只有精品99| a级片在线免费高清观看视频| 欧美日韩中文字幕国产精品一区二区三区 | 免费日韩欧美在线观看| 亚洲成人免费电影在线观看| 97在线人人人人妻| 精品一品国产午夜福利视频| 亚洲成av片中文字幕在线观看| 国产伦人伦偷精品视频| av网站免费在线观看视频| 国产欧美日韩一区二区精品| 国产97色在线日韩免费| 夫妻午夜视频| 亚洲情色 制服丝袜| 美国免费a级毛片| 另类精品久久| 自线自在国产av| 国产精品香港三级国产av潘金莲| 高清视频免费观看一区二区| 国产成人精品久久二区二区91| 欧美日韩一级在线毛片| 大片电影免费在线观看免费| 午夜免费鲁丝| 美女国产高潮福利片在线看| 不卡一级毛片| 亚洲国产毛片av蜜桃av| 久久天堂一区二区三区四区| 乱人伦中国视频| 国产xxxxx性猛交| 国产91精品成人一区二区三区 | 久久精品91无色码中文字幕| 99riav亚洲国产免费| 久久久久久久久久久久大奶| 欧美日韩成人在线一区二区| e午夜精品久久久久久久| 国产亚洲午夜精品一区二区久久| 国产精品国产高清国产av | 一边摸一边做爽爽视频免费| tube8黄色片| 久久久国产精品麻豆| 午夜福利在线观看吧| 亚洲综合色网址| 丰满迷人的少妇在线观看| 午夜91福利影院| 另类精品久久| 精品一品国产午夜福利视频| 中文字幕人妻丝袜制服| 欧美人与性动交α欧美精品济南到| 丰满人妻熟妇乱又伦精品不卡| 亚洲自偷自拍图片 自拍| 18禁观看日本| 午夜福利,免费看| 久久亚洲真实| 最新在线观看一区二区三区| 免费一级毛片在线播放高清视频 | 欧美精品啪啪一区二区三区| 999精品在线视频| 日本欧美视频一区| 一个人免费在线观看的高清视频| a在线观看视频网站| 亚洲第一欧美日韩一区二区三区 | 99国产极品粉嫩在线观看| 免费在线观看完整版高清| 亚洲黑人精品在线| 亚洲国产欧美在线一区| 精品国内亚洲2022精品成人 | 正在播放国产对白刺激| 在线十欧美十亚洲十日本专区| 欧美亚洲日本最大视频资源| 成年动漫av网址| 俄罗斯特黄特色一大片| av福利片在线| 国产片内射在线| 岛国在线观看网站| 国产亚洲精品第一综合不卡| 国产欧美日韩一区二区三| 国产xxxxx性猛交| 国产精品影院久久| 99久久人妻综合| 亚洲av成人不卡在线观看播放网| 悠悠久久av| 自线自在国产av| av线在线观看网站| 日韩制服丝袜自拍偷拍| 另类精品久久| 日本黄色视频三级网站网址 | www.熟女人妻精品国产| 亚洲avbb在线观看| 婷婷丁香在线五月| 午夜久久久在线观看| 桃花免费在线播放| bbb黄色大片| 午夜老司机福利片| 精品国产乱码久久久久久小说| 后天国语完整版免费观看| 久久久久国内视频| 中国美女看黄片| 电影成人av| 夜夜爽天天搞| 亚洲五月婷婷丁香| 夜夜夜夜夜久久久久| 菩萨蛮人人尽说江南好唐韦庄| 50天的宝宝边吃奶边哭怎么回事| 精品少妇一区二区三区视频日本电影| 国产xxxxx性猛交| 日韩欧美免费精品| bbb黄色大片| 看免费av毛片| 久久久精品区二区三区| h视频一区二区三区| 精品免费久久久久久久清纯 | 视频在线观看一区二区三区| 成人国语在线视频| 美女福利国产在线| 精品一区二区三区视频在线观看免费 | 成年动漫av网址| 欧美黑人欧美精品刺激| 人成视频在线观看免费观看| 精品久久久久久电影网| 中文字幕av电影在线播放| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利视频在线观看免费| 国产亚洲欧美在线一区二区| 欧美精品一区二区免费开放| 色尼玛亚洲综合影院| 我要看黄色一级片免费的| 精品亚洲乱码少妇综合久久| 亚洲精品美女久久久久99蜜臀| 久久亚洲真实| 精品欧美一区二区三区在线| 女人爽到高潮嗷嗷叫在线视频| 日韩欧美国产一区二区入口| 亚洲成av片中文字幕在线观看| 一边摸一边抽搐一进一小说 | 一区二区三区精品91| 老鸭窝网址在线观看| 欧美黑人欧美精品刺激| 国产视频一区二区在线看| 中文字幕人妻丝袜制服| 精品午夜福利视频在线观看一区 | 熟女少妇亚洲综合色aaa.| av欧美777| 中文字幕制服av| 日韩免费av在线播放| 日韩欧美国产一区二区入口| 乱人伦中国视频| 国产又爽黄色视频| 久久天堂一区二区三区四区| 99国产精品免费福利视频| 欧美日本中文国产一区发布| 久久精品亚洲av国产电影网| 两个人看的免费小视频| 动漫黄色视频在线观看| 在线天堂中文资源库| 日韩有码中文字幕| 国产1区2区3区精品| 午夜视频精品福利| 91精品三级在线观看| 亚洲色图 男人天堂 中文字幕| 久久99热这里只频精品6学生| 精品一区二区三区视频在线观看免费 | 亚洲男人天堂网一区| 制服诱惑二区| 国产精品电影一区二区三区 | 成人av一区二区三区在线看| 啦啦啦 在线观看视频| 99香蕉大伊视频| 真人做人爱边吃奶动态| 韩国精品一区二区三区| 国产精品.久久久| 国产成人一区二区三区免费视频网站| 午夜福利视频精品| 亚洲少妇的诱惑av| 夜夜爽天天搞| 国产成人av激情在线播放| 热re99久久精品国产66热6| 成人亚洲精品一区在线观看| 久久久精品国产亚洲av高清涩受| 久久国产亚洲av麻豆专区| 精品亚洲成a人片在线观看| 亚洲国产精品一区二区三区在线| 法律面前人人平等表现在哪些方面| 国产日韩欧美视频二区| 亚洲久久久国产精品| 视频区图区小说| 在线 av 中文字幕| 人妻 亚洲 视频| 汤姆久久久久久久影院中文字幕| 啦啦啦视频在线资源免费观看| 国产精品一区二区精品视频观看| 99精品欧美一区二区三区四区| 欧美日韩一级在线毛片| 久久精品亚洲精品国产色婷小说| 国产主播在线观看一区二区| 欧美成狂野欧美在线观看| 亚洲成人手机| 色在线成人网| 日韩人妻精品一区2区三区| 91大片在线观看| 蜜桃国产av成人99| 天天添夜夜摸| 日本vs欧美在线观看视频| 欧美 亚洲 国产 日韩一| 人妻 亚洲 视频| 国产成人影院久久av| 日日夜夜操网爽| 亚洲av成人不卡在线观看播放网| 亚洲国产精品一区二区三区在线| 91大片在线观看| 99久久人妻综合| 国产成人精品无人区| 久久影院123| 国产欧美日韩一区二区三区在线| 精品久久久久久久毛片微露脸| 少妇的丰满在线观看| 国产日韩欧美视频二区| 啦啦啦 在线观看视频| av网站在线播放免费| 纯流量卡能插随身wifi吗| 男女之事视频高清在线观看| 成人18禁高潮啪啪吃奶动态图| 丝袜在线中文字幕| 国产一区二区三区在线臀色熟女 | 国产在线观看jvid| 久久亚洲真实| svipshipincom国产片| 91字幕亚洲| 国产伦人伦偷精品视频| 制服人妻中文乱码| 18禁观看日本| 国产伦人伦偷精品视频| 国产高清国产精品国产三级| 狠狠精品人妻久久久久久综合| 中文字幕另类日韩欧美亚洲嫩草| 黄色丝袜av网址大全| 动漫黄色视频在线观看| 国产精品一区二区在线观看99| av又黄又爽大尺度在线免费看| 亚洲精品成人av观看孕妇| 日本黄色日本黄色录像| 精品福利永久在线观看| 久久 成人 亚洲| 一二三四在线观看免费中文在| 99在线人妻在线中文字幕 | 欧美日韩亚洲国产一区二区在线观看 | 视频在线观看一区二区三区| 看免费av毛片| 免费少妇av软件| 欧美日韩视频精品一区| 成人永久免费在线观看视频 | 国产精品亚洲av一区麻豆| 侵犯人妻中文字幕一二三四区| 成人免费观看视频高清| 美女扒开内裤让男人捅视频| 不卡一级毛片| 久久中文字幕人妻熟女| 午夜激情久久久久久久| 中文亚洲av片在线观看爽 | 成人免费观看视频高清| e午夜精品久久久久久久| 亚洲欧美日韩另类电影网站| 国产欧美日韩精品亚洲av| 欧美一级毛片孕妇| 激情视频va一区二区三区| 亚洲av电影在线进入| 法律面前人人平等表现在哪些方面| 丝袜在线中文字幕| 亚洲精品美女久久久久99蜜臀| 亚洲天堂av无毛| 不卡av一区二区三区| 亚洲成av片中文字幕在线观看| 日韩大片免费观看网站| 动漫黄色视频在线观看| 欧美激情高清一区二区三区| 最近最新中文字幕大全电影3 | 日韩有码中文字幕| 亚洲午夜精品一区,二区,三区| 99在线人妻在线中文字幕 | 一级片'在线观看视频| 欧美亚洲日本最大视频资源| 日韩欧美免费精品| 国产精品一区二区免费欧美| 久久精品熟女亚洲av麻豆精品| 一区在线观看完整版| 变态另类成人亚洲欧美熟女 | 99精国产麻豆久久婷婷| 久久影院123| 97在线人人人人妻| 大香蕉久久成人网| 亚洲va日本ⅴa欧美va伊人久久| 99在线人妻在线中文字幕 | 法律面前人人平等表现在哪些方面| 久久久久久久国产电影| 久久久久久久久久久久大奶| 午夜福利视频精品| 两个人免费观看高清视频| 性少妇av在线| 老汉色∧v一级毛片| 国产福利在线免费观看视频| 少妇猛男粗大的猛烈进出视频| 久久人妻熟女aⅴ| 捣出白浆h1v1| 国产激情久久老熟女| 国产在线一区二区三区精| 丝袜人妻中文字幕| 亚洲专区字幕在线| 亚洲久久久国产精品| 2018国产大陆天天弄谢| 欧美+亚洲+日韩+国产| 久久性视频一级片| 久久久久网色| 成在线人永久免费视频| 50天的宝宝边吃奶边哭怎么回事| 欧美在线黄色| 99精品久久久久人妻精品| 麻豆乱淫一区二区| 丝袜美腿诱惑在线| 日日夜夜操网爽| 久久久欧美国产精品| 两人在一起打扑克的视频| 成人特级黄色片久久久久久久 | 日本一区二区免费在线视频| 激情在线观看视频在线高清 | 亚洲熟妇熟女久久| www.熟女人妻精品国产| 男女午夜视频在线观看| 久久香蕉激情| 麻豆乱淫一区二区| 母亲3免费完整高清在线观看| 热re99久久国产66热| 国产精品免费一区二区三区在线 | 亚洲熟女精品中文字幕| 黄色成人免费大全| 狠狠狠狠99中文字幕| 欧美日韩国产mv在线观看视频| 国产在线免费精品| 免费在线观看日本一区| 国产精品久久久人人做人人爽| 国产成人欧美在线观看 | 国产野战对白在线观看| 欧美日韩视频精品一区| 变态另类成人亚洲欧美熟女 | 精品少妇黑人巨大在线播放| 午夜激情久久久久久久| 午夜福利视频在线观看免费| 极品人妻少妇av视频| 天天躁夜夜躁狠狠躁躁| 成人手机av| 午夜激情av网站| 国产伦人伦偷精品视频| 王馨瑶露胸无遮挡在线观看| 国产亚洲欧美在线一区二区| 免费在线观看日本一区| 亚洲第一青青草原| 色尼玛亚洲综合影院| av欧美777| 狂野欧美激情性xxxx| 亚洲美女黄片视频| 日韩一卡2卡3卡4卡2021年| 狠狠精品人妻久久久久久综合| 成人黄色视频免费在线看| 亚洲三区欧美一区| 男女边摸边吃奶| 老司机在亚洲福利影院| 国产有黄有色有爽视频| 亚洲欧美一区二区三区久久| 大香蕉久久网| 99精品在免费线老司机午夜| 亚洲久久久国产精品| av一本久久久久| 国产av一区二区精品久久| 免费久久久久久久精品成人欧美视频| 国产野战对白在线观看| 欧美日韩一级在线毛片| 美女视频免费永久观看网站| 久久久久国内视频| 九色亚洲精品在线播放| 老鸭窝网址在线观看| 欧美精品亚洲一区二区| 亚洲黑人精品在线| 国产在视频线精品| 久久久国产欧美日韩av| 这个男人来自地球电影免费观看| 老司机在亚洲福利影院| 夜夜爽天天搞| 曰老女人黄片| 亚洲av日韩在线播放| 美女国产高潮福利片在线看| 中文字幕另类日韩欧美亚洲嫩草| 王馨瑶露胸无遮挡在线观看| 国产av精品麻豆| 99九九在线精品视频| 亚洲性夜色夜夜综合| 美女福利国产在线| 久久久国产精品麻豆| 99九九在线精品视频| 欧美激情高清一区二区三区| 欧美人与性动交α欧美软件| 少妇裸体淫交视频免费看高清 | 亚洲欧美一区二区三区久久| 人人妻人人澡人人爽人人夜夜| 黄色毛片三级朝国网站| 五月开心婷婷网| 啦啦啦视频在线资源免费观看| 国产一区二区激情短视频| 黄色毛片三级朝国网站| 亚洲av日韩精品久久久久久密| 丁香六月欧美| 国产成人精品在线电影| 王馨瑶露胸无遮挡在线观看| 老司机影院毛片| 亚洲专区字幕在线| 18禁美女被吸乳视频| 下体分泌物呈黄色| 久久天躁狠狠躁夜夜2o2o| 久久国产精品影院| 久久精品成人免费网站| 亚洲熟女毛片儿| 成年人黄色毛片网站| 亚洲精品成人av观看孕妇| 黄片大片在线免费观看| 黄色 视频免费看| 一区在线观看完整版| 狠狠精品人妻久久久久久综合| 精品国产乱码久久久久久小说| 日韩制服丝袜自拍偷拍| 日本撒尿小便嘘嘘汇集6| 国产一区有黄有色的免费视频| 国产区一区二久久| 国产一区有黄有色的免费视频| 熟女少妇亚洲综合色aaa.| 国产高清激情床上av| 亚洲精品国产色婷婷电影| 久久av网站| 怎么达到女性高潮| 日韩三级视频一区二区三区| 免费一级毛片在线播放高清视频 | 交换朋友夫妻互换小说| 欧美激情久久久久久爽电影 | 久久久久视频综合| 午夜福利乱码中文字幕| 黄色丝袜av网址大全| 欧美激情 高清一区二区三区| 国产av国产精品国产| 老司机影院毛片| 老司机亚洲免费影院| 亚洲精品自拍成人| 少妇 在线观看| 1024视频免费在线观看| 国产日韩欧美视频二区| 成人黄色视频免费在线看| 国产欧美日韩一区二区三区在线| 国产精品秋霞免费鲁丝片| 免费少妇av软件| 亚洲精品国产色婷婷电影| 99精品欧美一区二区三区四区| 窝窝影院91人妻| 大陆偷拍与自拍| 又黄又粗又硬又大视频| 国产精品av久久久久免费| 在线十欧美十亚洲十日本专区| 老鸭窝网址在线观看| 丰满饥渴人妻一区二区三| 久久久国产成人免费| videos熟女内射| 大片电影免费在线观看免费| 欧美午夜高清在线| 国精品久久久久久国模美| 亚洲九九香蕉| 午夜激情av网站| 两性夫妻黄色片| 丝瓜视频免费看黄片| 日韩大片免费观看网站| 久久精品aⅴ一区二区三区四区| 9色porny在线观看| 视频区欧美日本亚洲| 一本大道久久a久久精品| cao死你这个sao货| 建设人人有责人人尽责人人享有的| 69精品国产乱码久久久|