• <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
    久久毛片免费看一区二区三区| 色吧在线观看| 亚洲av成人精品一二三区| 国产麻豆69| 欧美日韩一级在线毛片| 午夜日本视频在线| 国产精品偷伦视频观看了| 一区二区三区精品91| 一级,二级,三级黄色视频| 久久久久久久久久久久大奶| 伊人亚洲综合成人网| 国产精品国产三级专区第一集| 人人妻人人澡人人看| av一本久久久久| 亚洲婷婷狠狠爱综合网| 精品亚洲乱码少妇综合久久| 欧美日韩综合久久久久久| 精品国产乱码久久久久久小说| 国产男女内射视频| 一区二区三区乱码不卡18| 免费黄网站久久成人精品| 亚洲在久久综合| 老鸭窝网址在线观看| 亚洲欧洲国产日韩| 久久久国产一区二区| 国产成人精品婷婷| 久久久久精品性色| 日韩一区二区视频免费看| 久久精品国产亚洲av高清一级| 少妇被粗大猛烈的视频| 啦啦啦在线免费观看视频4| 婷婷色av中文字幕| 伦理电影大哥的女人| 69精品国产乱码久久久| 啦啦啦视频在线资源免费观看| 亚洲av综合色区一区| 国产伦理片在线播放av一区| 观看av在线不卡| 国产成人av激情在线播放| 秋霞在线观看毛片| 国产毛片在线视频| 中文字幕人妻丝袜制服| 99香蕉大伊视频| 人人澡人人妻人| 日韩中文字幕欧美一区二区 | 欧美日韩视频高清一区二区三区二| 在线观看免费高清a一片| 欧美日韩亚洲高清精品| 亚洲中文av在线| 国产极品粉嫩免费观看在线| 美女福利国产在线| 国产日韩欧美亚洲二区| 国产女主播在线喷水免费视频网站| 午夜激情久久久久久久| 韩国av在线不卡| 热99久久久久精品小说推荐| 哪个播放器可以免费观看大片| 久久久久久久精品精品| 秋霞在线观看毛片| 人妻系列 视频| 韩国精品一区二区三区| 色网站视频免费| 亚洲伊人久久精品综合| 国产精品嫩草影院av在线观看| 亚洲欧美色中文字幕在线| 你懂的网址亚洲精品在线观看| 99久久中文字幕三级久久日本| av电影中文网址| 国产日韩欧美亚洲二区| 日韩人妻精品一区2区三区| 欧美中文综合在线视频| 性色av一级| 99国产精品免费福利视频| 香蕉国产在线看| 日本av手机在线免费观看| 美女视频免费永久观看网站| 五月开心婷婷网| 少妇人妻 视频| 亚洲婷婷狠狠爱综合网| 在线观看免费视频网站a站| 日韩人妻精品一区2区三区| 欧美少妇被猛烈插入视频| 最近的中文字幕免费完整| 午夜免费男女啪啪视频观看| 婷婷色综合大香蕉| 最近手机中文字幕大全| 国产免费又黄又爽又色| 黑丝袜美女国产一区| 亚洲综合色惰| 春色校园在线视频观看| 女人被躁到高潮嗷嗷叫费观| 国产一区二区三区综合在线观看| 亚洲伊人久久精品综合| 免费在线观看视频国产中文字幕亚洲 | 久久婷婷青草| 丝袜人妻中文字幕| 久久久久久久精品精品| 国产乱人偷精品视频| 国产成人精品福利久久| 亚洲欧洲国产日韩| av女优亚洲男人天堂| 日韩一区二区视频免费看| 久久久久久久久久久久大奶| 在线观看国产h片| 久久久精品区二区三区| 夫妻性生交免费视频一级片| 日本猛色少妇xxxxx猛交久久| 亚洲欧美精品综合一区二区三区 | 久久久久视频综合| 欧美日韩av久久| av视频免费观看在线观看| 国产精品香港三级国产av潘金莲 | 国产白丝娇喘喷水9色精品| 最新中文字幕久久久久| 一本久久精品| 久久精品aⅴ一区二区三区四区 | a级毛片黄视频| 精品一区二区三卡| 国产精品一国产av| 美国免费a级毛片| 青春草视频在线免费观看| 亚洲激情五月婷婷啪啪| 在线观看美女被高潮喷水网站| 国产精品蜜桃在线观看| 在线看a的网站| 国产亚洲午夜精品一区二区久久| 激情视频va一区二区三区| 看免费成人av毛片| 亚洲av福利一区| 狂野欧美激情性bbbbbb| 久久精品国产亚洲av涩爱| 一级,二级,三级黄色视频| 免费播放大片免费观看视频在线观看| 大片电影免费在线观看免费| 啦啦啦中文免费视频观看日本| 中文乱码字字幕精品一区二区三区| 国产综合精华液| 毛片一级片免费看久久久久| 日本爱情动作片www.在线观看| 中文乱码字字幕精品一区二区三区| 三上悠亚av全集在线观看| 久久精品aⅴ一区二区三区四区 | 大码成人一级视频| 国产1区2区3区精品| 国产老妇伦熟女老妇高清| av线在线观看网站| 亚洲,欧美,日韩| av网站免费在线观看视频| 人人妻人人澡人人看| 亚洲情色 制服丝袜| 欧美精品亚洲一区二区| 日本91视频免费播放| 精品酒店卫生间| 超色免费av| 久久久国产一区二区| 国产1区2区3区精品| 在线免费观看不下载黄p国产| 91成人精品电影| 欧美激情极品国产一区二区三区| 亚洲精品视频女| 可以免费在线观看a视频的电影网站 | 国产精品成人在线| 亚洲av福利一区| 日韩欧美精品免费久久| 久久久久久人妻| 宅男免费午夜| 一级毛片黄色毛片免费观看视频| 18在线观看网站| 国产亚洲欧美精品永久| 精品国产国语对白av| 香蕉精品网在线| 熟女av电影| a级毛片在线看网站| av网站在线播放免费| 卡戴珊不雅视频在线播放| 老司机影院成人| 女人高潮潮喷娇喘18禁视频| 亚洲精品美女久久av网站| 18在线观看网站| 一级毛片电影观看| 曰老女人黄片| 日韩视频在线欧美| 精品国产乱码久久久久久男人| 视频区图区小说| 日日爽夜夜爽网站| 国产成人欧美| 欧美人与善性xxx| 黄片无遮挡物在线观看| 久久热在线av| 精品少妇一区二区三区视频日本电影 | 9色porny在线观看| 97精品久久久久久久久久精品| 日韩精品免费视频一区二区三区| 制服丝袜香蕉在线| 久久久久久伊人网av| 在线观看免费高清a一片| 欧美日韩亚洲高清精品| 99国产精品免费福利视频| 街头女战士在线观看网站| 亚洲欧洲精品一区二区精品久久久 | 日韩精品有码人妻一区| 最近中文字幕高清免费大全6| 亚洲精品在线美女| 在线精品无人区一区二区三| av在线观看视频网站免费| 搡老乐熟女国产| 一二三四在线观看免费中文在| 国产av一区二区精品久久| 亚洲精品自拍成人| 精品一区在线观看国产| 久久精品夜色国产| 1024香蕉在线观看| 亚洲精品自拍成人| 久久午夜综合久久蜜桃| 亚洲精品成人av观看孕妇| 久久毛片免费看一区二区三区| 男女啪啪激烈高潮av片| 男女高潮啪啪啪动态图| 国产精品免费视频内射| 成人午夜精彩视频在线观看| 婷婷色麻豆天堂久久| 亚洲综合色网址| 亚洲国产精品成人久久小说| 久久狼人影院| 交换朋友夫妻互换小说| av一本久久久久| 久久人人爽人人片av| 国产精品国产av在线观看| 亚洲一区二区三区欧美精品| 国产视频首页在线观看| 97精品久久久久久久久久精品| 国产成人午夜福利电影在线观看| freevideosex欧美| 国产日韩一区二区三区精品不卡| 久久久精品区二区三区| 中文字幕人妻丝袜一区二区 | 少妇被粗大猛烈的视频| av国产精品久久久久影院| 老汉色∧v一级毛片| 久久久久国产网址| 亚洲欧洲国产日韩| 9191精品国产免费久久| 伊人久久国产一区二区| 日本免费在线观看一区| 午夜激情久久久久久久| 777久久人妻少妇嫩草av网站| 欧美97在线视频| 一级毛片黄色毛片免费观看视频| 亚洲av.av天堂| 啦啦啦视频在线资源免费观看| 亚洲精品国产av蜜桃| 精品少妇黑人巨大在线播放| freevideosex欧美| 中文字幕最新亚洲高清| 97精品久久久久久久久久精品| 久久婷婷青草| 成年美女黄网站色视频大全免费| 精品国产一区二区久久| 国产免费视频播放在线视频| 精品国产超薄肉色丝袜足j| 欧美激情 高清一区二区三区| 人成视频在线观看免费观看| 国产亚洲午夜精品一区二区久久| www日本在线高清视频| 亚洲伊人色综图| 国产熟女欧美一区二区| 成人国产av品久久久| 国产精品成人在线| 欧美日韩亚洲国产一区二区在线观看 | 蜜桃国产av成人99| 99国产精品免费福利视频| 日本色播在线视频| 免费观看性生交大片5| 飞空精品影院首页| 黄色配什么色好看| 国产熟女午夜一区二区三区| 国产成人av激情在线播放| 韩国av在线不卡| 国产av一区二区精品久久| av网站免费在线观看视频| 亚洲美女视频黄频| 搡老乐熟女国产| 中文乱码字字幕精品一区二区三区| 亚洲第一av免费看| 狠狠婷婷综合久久久久久88av| 日韩伦理黄色片| 亚洲精品av麻豆狂野| 欧美亚洲日本最大视频资源| 国产伦理片在线播放av一区| 91午夜精品亚洲一区二区三区| 欧美国产精品一级二级三级| 国产精品久久久久久精品电影小说| 国产亚洲av片在线观看秒播厂| 黄色一级大片看看| 日本色播在线视频| 搡老乐熟女国产| 在线免费观看不下载黄p国产| 蜜桃在线观看..| 日韩人妻精品一区2区三区| 一级片'在线观看视频| 亚洲一区中文字幕在线| 秋霞伦理黄片| 精品福利永久在线观看| 熟女少妇亚洲综合色aaa.| 午夜激情av网站| 国产一区二区三区综合在线观看| 亚洲综合色惰| 精品少妇一区二区三区视频日本电影 | a级毛片黄视频| 国产av一区二区精品久久| 日韩av在线免费看完整版不卡| 亚洲欧洲精品一区二区精品久久久 | 极品人妻少妇av视频| 精品国产国语对白av| 巨乳人妻的诱惑在线观看| 18在线观看网站| 啦啦啦视频在线资源免费观看| 国产爽快片一区二区三区| 在线天堂最新版资源| 免费久久久久久久精品成人欧美视频| 韩国av在线不卡| 成人亚洲精品一区在线观看| 又粗又硬又长又爽又黄的视频| 国产无遮挡羞羞视频在线观看| 交换朋友夫妻互换小说| 午夜激情av网站| 精品国产乱码久久久久久男人| 美女xxoo啪啪120秒动态图| 午夜福利乱码中文字幕| 欧美国产精品一级二级三级| 亚洲人成网站在线观看播放| 国产人伦9x9x在线观看 | 亚洲,一卡二卡三卡| 伦精品一区二区三区| 午夜久久久在线观看| 成人二区视频| 国产在线免费精品| 亚洲色图 男人天堂 中文字幕| 我的亚洲天堂| 亚洲精品一区蜜桃| 亚洲一码二码三码区别大吗| 国产成人av激情在线播放| 老熟女久久久| 久久免费观看电影| 亚洲 欧美一区二区三区| 久久久国产欧美日韩av| 久久人妻熟女aⅴ| 国产激情久久老熟女| 久热久热在线精品观看| 熟女电影av网| 夫妻午夜视频| 精品亚洲乱码少妇综合久久| 亚洲av欧美aⅴ国产| 丰满少妇做爰视频| 中文字幕人妻丝袜制服| www日本在线高清视频| 久久人妻熟女aⅴ| 高清欧美精品videossex| 九草在线视频观看| 男女国产视频网站| 老熟女久久久| 麻豆av在线久日| 亚洲国产欧美网| av视频免费观看在线观看| 精品国产一区二区久久| 97人妻天天添夜夜摸| 五月开心婷婷网| 国产视频首页在线观看| 精品亚洲乱码少妇综合久久| freevideosex欧美| av电影中文网址| 亚洲成av片中文字幕在线观看 | 超碰成人久久| 精品人妻熟女毛片av久久网站| av又黄又爽大尺度在线免费看| 999精品在线视频| 国产女主播在线喷水免费视频网站| 国产成人a∨麻豆精品| 亚洲欧美色中文字幕在线| 日本av免费视频播放| 三级国产精品片| 国产精品偷伦视频观看了| 在线观看www视频免费| 街头女战士在线观看网站| 免费看av在线观看网站| 最黄视频免费看| 中文字幕制服av| 日韩精品有码人妻一区| 新久久久久国产一级毛片| 女人精品久久久久毛片| 美女午夜性视频免费| 欧美日韩亚洲国产一区二区在线观看 | 人人妻人人澡人人看| av不卡在线播放| 国产伦理片在线播放av一区| 日韩制服骚丝袜av| 欧美 亚洲 国产 日韩一| 亚洲在久久综合| 男女高潮啪啪啪动态图| 99久久人妻综合| 91国产中文字幕| 中文乱码字字幕精品一区二区三区| 国产在线免费精品| 久久ye,这里只有精品| av在线老鸭窝| 美女主播在线视频| 侵犯人妻中文字幕一二三四区| 99国产综合亚洲精品| 天美传媒精品一区二区| 成年动漫av网址| 亚洲av成人精品一二三区| 99热全是精品| 精品国产乱码久久久久久男人| 亚洲av电影在线观看一区二区三区| av在线老鸭窝| 视频区图区小说| 久久精品久久久久久久性| 交换朋友夫妻互换小说| 精品久久久精品久久久| 看免费成人av毛片| 国产精品一区二区在线不卡| 天天影视国产精品| 久久精品久久久久久噜噜老黄| 欧美成人午夜精品| 午夜福利在线观看免费完整高清在| 中文字幕av电影在线播放| 18禁国产床啪视频网站| 日韩av免费高清视频| 男女午夜视频在线观看| 国产黄频视频在线观看| 中文乱码字字幕精品一区二区三区| 久久精品人人爽人人爽视色| 国产亚洲午夜精品一区二区久久| 精品人妻一区二区三区麻豆| 日本午夜av视频| 久久人人爽人人片av| 国产福利在线免费观看视频| 亚洲国产最新在线播放| 日韩,欧美,国产一区二区三区| 国产亚洲最大av| 在线观看一区二区三区激情| 免费播放大片免费观看视频在线观看| 欧美成人午夜精品| 交换朋友夫妻互换小说| av电影中文网址| 亚洲中文av在线| 人人澡人人妻人| av视频免费观看在线观看| 亚洲国产欧美在线一区| 大片电影免费在线观看免费| 国产精品久久久久久久久免| 国产黄色免费在线视频| 久久久久网色| 亚洲成人av在线免费| 亚洲精品av麻豆狂野| 欧美黄色片欧美黄色片| 免费av中文字幕在线| 欧美97在线视频| 一本色道久久久久久精品综合| 天天躁夜夜躁狠狠躁躁| 亚洲熟女精品中文字幕| 日韩一区二区三区影片| 欧美精品人与动牲交sv欧美| 欧美日韩国产mv在线观看视频| 国产探花极品一区二区| 久久午夜综合久久蜜桃| 久久这里有精品视频免费| 大码成人一级视频| 91成人精品电影| 精品一区二区三区四区五区乱码 | 91午夜精品亚洲一区二区三区| 老司机影院成人| 欧美激情高清一区二区三区 | 欧美精品亚洲一区二区| 永久免费av网站大全| av不卡在线播放| 精品国产露脸久久av麻豆| 亚洲国产欧美网| 性高湖久久久久久久久免费观看| 久久99热这里只频精品6学生| 国产精品二区激情视频| 国产精品女同一区二区软件| 18+在线观看网站| 国产精品欧美亚洲77777| 欧美97在线视频| av又黄又爽大尺度在线免费看| 丁香六月天网| 日韩大片免费观看网站| 97在线人人人人妻| 极品人妻少妇av视频| 自线自在国产av| 国产欧美日韩综合在线一区二区| 国产一级毛片在线| 如日韩欧美国产精品一区二区三区| 三级国产精品片| 亚洲精品久久久久久婷婷小说| 一级a爱视频在线免费观看| 天天操日日干夜夜撸| 亚洲精品日本国产第一区| 天堂中文最新版在线下载| 中文字幕精品免费在线观看视频| 久久国产精品大桥未久av| 国产精品久久久久久av不卡| 女性生殖器流出的白浆| 久久国产亚洲av麻豆专区| 不卡视频在线观看欧美| 一区二区三区四区激情视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 老司机亚洲免费影院| 欧美日韩一区二区视频在线观看视频在线| 咕卡用的链子| 丝袜美腿诱惑在线| 一区二区三区乱码不卡18| 精品国产乱码久久久久久男人| 午夜福利影视在线免费观看| 一级毛片黄色毛片免费观看视频| 亚洲精品国产av成人精品| 欧美日韩一区二区视频在线观看视频在线| 人妻人人澡人人爽人人| 国产有黄有色有爽视频| 最近手机中文字幕大全| 制服丝袜香蕉在线| 中文字幕精品免费在线观看视频| 免费观看av网站的网址| 欧美 亚洲 国产 日韩一| 亚洲国产精品一区二区三区在线| 亚洲av免费高清在线观看| 国产福利在线免费观看视频| 亚洲国产色片| 一级,二级,三级黄色视频| 免费黄色在线免费观看| 最新的欧美精品一区二区| 韩国高清视频一区二区三区| 国产一区二区三区av在线| 精品99又大又爽又粗少妇毛片| 欧美日韩一级在线毛片| 女人被躁到高潮嗷嗷叫费观| 26uuu在线亚洲综合色| 丝瓜视频免费看黄片| 国精品久久久久久国模美| 五月开心婷婷网| 性色avwww在线观看| 成年美女黄网站色视频大全免费| 中文乱码字字幕精品一区二区三区| 交换朋友夫妻互换小说| 一边摸一边做爽爽视频免费| 91aial.com中文字幕在线观看| 欧美成人精品欧美一级黄| 26uuu在线亚洲综合色| 国产乱人偷精品视频| 日韩三级伦理在线观看| 精品一品国产午夜福利视频| 精品一区在线观看国产| 人人妻人人爽人人添夜夜欢视频| 亚洲 欧美一区二区三区| 国产精品一区二区在线观看99| 一区二区三区乱码不卡18| 国产福利在线免费观看视频| 一级a爱视频在线免费观看| 两性夫妻黄色片| 免费黄网站久久成人精品| 亚洲精品乱久久久久久| 人妻 亚洲 视频| 成人18禁高潮啪啪吃奶动态图| 老汉色av国产亚洲站长工具| 啦啦啦在线免费观看视频4| 久久久久久久久久人人人人人人| 超碰成人久久| 亚洲精品自拍成人| 久热久热在线精品观看| 日韩一卡2卡3卡4卡2021年| 伊人久久国产一区二区| 亚洲一级一片aⅴ在线观看| 国产精品麻豆人妻色哟哟久久| 久久久国产一区二区| 亚洲精华国产精华液的使用体验| 观看美女的网站| av有码第一页| 婷婷色综合www| 永久网站在线| 亚洲五月色婷婷综合| 国产精品一国产av| 侵犯人妻中文字幕一二三四区| 国产熟女午夜一区二区三区| 国产一级毛片在线| 久久这里只有精品19| 国产免费又黄又爽又色| 精品一区在线观看国产| 老汉色av国产亚洲站长工具| 精品少妇内射三级| 寂寞人妻少妇视频99o| 超碰成人久久| 亚洲av中文av极速乱| 飞空精品影院首页| 国产成人免费观看mmmm| 精品一区二区三卡| 9色porny在线观看| 成人亚洲欧美一区二区av| 欧美日韩精品网址| 午夜av观看不卡| 黑人巨大精品欧美一区二区蜜桃| 国产精品人妻久久久影院| 80岁老熟妇乱子伦牲交| 久久青草综合色| 精品国产露脸久久av麻豆| 午夜激情av网站| 成人黄色视频免费在线看| 国产在线一区二区三区精| 大片电影免费在线观看免费| 久久亚洲国产成人精品v| 在线观看免费高清a一片| 久久精品国产亚洲av高清一级| 秋霞在线观看毛片| 久久久久精品久久久久真实原创| 国产精品久久久av美女十八|