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

    黃渤海氮磷營養(yǎng)鹽的分布、收支與生態(tài)環(huán)境效應

    2016-10-13 21:37:25趙晨英臧家業(yè)冉祥濱
    中國環(huán)境科學 2016年7期
    關鍵詞:黃渤海營養(yǎng)鹽黃海

    趙晨英,臧家業(yè),劉 軍,2,孫 濤,冉祥濱*

    ?

    趙晨英1,臧家業(yè)1,劉 軍1,2,孫 濤1,冉祥濱1*

    (1.國家海洋局第一海洋研究所生態(tài)研究中心,山東 青島 266061;2.中國海洋大學,海洋化學理論與工程技術教育部重點實驗室,山東 青島 266100)

    基于在黃渤海的綜合調查結果,分析了水體和沉積物間隙水中溶解無機氮(DIN)和溶解無機磷(DIP)的分布;結合歷史數(shù)據(jù)構建了黃渤海DIN和DIP的收支模型,并分析了陸源輸入變化對研究區(qū)域生態(tài)環(huán)境的影響.結果表明,黃渤海DIN和DIP的含量受季節(jié)、河流輸入和沉積物界面擴散作用的影響,具有秋季高于春季和近岸高于離岸的時空分布特征.收支模型計算結果表明,底界面擴散是黃渤海水體DIN的主要來源,其次是大氣、周邊河流、地下水和東海的輸入;黃渤海水體DIN的支出主要是通過沉積埋藏和反硝化.黃渤海水體DIP的來源主要是磷酸鹽吸附解吸,占91%,底界面擴散和大氣輸入為其次,河流和地下水的輸入貢獻較小.DIP的支出主要是通過沉積埋藏和向東海的輸出.黃渤海每年有11Gmol的氮在水體積累,并導致其濃度提高約0.6μmol/(L·a).近些年來陸地向黃渤海輸入氮的持續(xù)增加,加劇了氮營養(yǎng)鹽的積累,導致非硅藻類浮游植物比例以及赤潮發(fā)生頻率和面積顯著增加,同時還提高了水體初級生產(chǎn)力和海洋磷的埋藏量以及加劇了磷限制的趨勢,并可能威脅生態(tài)系統(tǒng)的穩(wěn)定.

    黃渤海;溶解無機氮;溶解無機磷;收支

    營養(yǎng)鹽的濃度和結構深刻影響著海洋浮游植物初級生產(chǎn)力和種群結構[1-2].近年來,人類活動改變了陸架邊緣海營養(yǎng)鹽的水平和循環(huán)模導致近海富營養(yǎng)化加重,有害赤潮頻發(fā),以及浮游生物群落改變等[3-7],這嚴重威脅甚至破壞了海洋生態(tài)系統(tǒng)的穩(wěn)定.因此,陸架邊緣海營養(yǎng)鹽的生物地球化學循環(huán)和收支過程的研究是海洋環(huán)境變化研究和管理以及海洋環(huán)境治理的基礎,也是碳循環(huán)研究中非常重要的一個方面[1,8],為當前全球環(huán)境變化領域研究的熱點.

    盡管陸架邊緣海區(qū)僅占全球海洋面積的8%左右[9],但該區(qū)域初級生產(chǎn)水平顯著高于大洋[3],在物質循環(huán)和全球氣候變化方面發(fā)揮著重要的作用[10-12].然而,近年來陸地地表過程變化和人類活動增強顯著改變了邊緣海獲取營養(yǎng)鹽的總量和生態(tài)功能[13-15],使之成為全球環(huán)境變化研究熱點區(qū)域.渤海和黃海作為典型的半封閉型陸架邊緣海[16],是人類活動和地表變化影響最為突出的海區(qū)之一.近年來日益增強的人類活動顯著改變了黃渤海入海河流營養(yǎng)鹽的含量和結構[13,17-21],這必然改變該海域營養(yǎng)鹽的水平和結構[2,22-24],從而打破水體原有的營養(yǎng)鹽平衡,造成黃渤海的浮游植物群落組成和初級生產(chǎn)的變化[2,22,25],產(chǎn)生深遠的生態(tài)環(huán)境效應[2,18,22,25].更為關鍵的是,由于黃渤海水體停留時間較長[22,26],其內部營養(yǎng)鹽水平和結構的變化對生態(tài)系統(tǒng)的影響可能比其它開放水域更深遠.因此,我們有必要對黃渤海營養(yǎng)鹽的來源、循環(huán)過程和生態(tài)效應等進行深入的研究.但目前已有的收支研究多局限于有限的時間或空間范圍內,缺少營養(yǎng)鹽積累和反硝化作用等關鍵過程以及陸源輸入在近海生態(tài)環(huán)境中所發(fā)揮的作用等方面的量化分析.在陸地營養(yǎng)鹽輸入變化和海洋營養(yǎng)鹽限制日益突出的今天,有關黃渤海氮(N)與磷(P)營養(yǎng)鹽的收支及控制機制方面的研究則顯得尤為重要.

    本文根據(jù)黃渤海、黃河和長江的綜合調查以及相關的歷史資料,分析了黃渤海氮和磷的組成與分布,建立了黃渤海氮-磷的收支模型,討論了影響該區(qū)域氮和磷分布、界面通量和埋藏的主要生物地球化學循環(huán)過程,探討了河流物質輸入變化對黃渤海氮-磷循環(huán)和生態(tài)環(huán)境造成的影響,以期為全面理解黃渤海乃至中國東部陸架海域環(huán)境演化預測提供依據(jù).

    1 材料與方法

    1.1 樣品采集與分析

    于2012年春季(5月3~24日)和秋季(11月2~20日)搭載“東方紅2”號調查船在黃渤海進行了海洋科學綜合調查,采樣站位如圖1所示.使用溫鹽深儀(Seabird 911CTD Plus)采集水樣,同時獲取水體溫度與鹽度等水文參數(shù);部分站位采集了表層沉積物和柱狀沉積物樣品.

    取500mL的水樣用孔徑0.45μm聚醚砜(預先用體積比1:1000HCl溶液浸泡24h,并以Milli-Q水洗至中性,45℃烘干稱重)過濾,濾液分裝于聚乙烯樣品瓶中(樣品瓶預先在1:5HCl中浸泡48h以上,用Milli-Q水清洗數(shù)遍)中,20℃冷凍保存,用于硝酸鹽(NO3-)、亞硝酸鹽(NO2-)、銨鹽(NH4+)與無機磷酸鹽(DIP)營養(yǎng)鹽的測定.現(xiàn)場用箱式采泥器采集表層沉積物,上覆水采集后,取表層沉積物(0~1cm)于密封袋中,冷凍保存.部分站位在采集表層沉積物的同時采集了柱狀沉積物樣品(采樣管內徑=9cm),現(xiàn)場分割取樣,取樣間隔為1cm,并取50g新鮮沉積物離心過濾制得間隙水.上覆水和間隙水樣品處理和保存方法同水樣,用于分析溶解無機氮(DIN)和溶解無機磷(DIP),沉積物烘干后用于有機氮(ON)的分析.

    DIN和DIP利用QUTRRAO(德國SEAL公司)自動分析儀進行測定,相對標準偏差<5%. DIN=NO3-+NO2-+NH4+.沉積物有機氮分析采用元素分析儀(利曼公司,Euro Vector EA3000)進行測定;分析前沉積物先用4mol/L的HCl清洗后,洗至中性并低溫烘干,上機測試;有機氮標準為C8H9NO(E1101,Euro Vector S.p.A,Italy),相對偏差<10%.

    1.2 收支計算

    黃渤海的氮和磷收支評估為穩(wěn)態(tài)零維箱式模型,即以水體為收支計算的核心,考慮主要的水文、化學與生物過程,將和DIP的外源輸入(河流輸入、大氣干濕沉降、東海的輸入、底界面的擴散和地下水輸入等)以及支出(黃海向東海的輸出、反硝化作用和沉積埋藏)以及內部過程(初級生產(chǎn)消耗和降解作用)等納入到模型中,氮和磷在初級生產(chǎn)、沉積和內部循環(huán)過程中遵循Redfield比值等基本原則[28].水的收支采用劉軍等[27]給出的黃渤海水量收支計算.主要過程的計算方法如下:

    1.2.1 初級生產(chǎn)(P) 由于缺少近些年系統(tǒng)的初級生產(chǎn)的實測數(shù)據(jù),黃渤海全年的初級生產(chǎn)力(以碳計)根據(jù)1984~1985年黃渤海的季節(jié)性(2月、5月、8月和11月)調查資料[29],以及利用分級初級生產(chǎn)力模式反演2003~2005年黃渤海年平均初級生產(chǎn)力[30],綜合確定并結合Redfield比值(C:N:P= 106:16:1)估算出浮游植物從海水中吸收DIN和DIP的速率,再根據(jù)黃渤海的水域面積(黃海: 38.0×1010m2,渤海:7.7×1010m2[16])計算因初級生產(chǎn)力所消耗的DIN和DIP的總量.

    1.2.2 河流輸入(R) 黃渤海沿岸陸源輸入主要考慮影響較大的12條河流,根據(jù)多年的水文數(shù)據(jù)和營養(yǎng)鹽監(jiān)測資料,通過累加得到黃海與渤海主要入海河流的DIN和DIP的入海通量.其中,黃河和長江的營養(yǎng)鹽來源于本課題組2013~ 2014年在長江(觀測點:長江江陰大橋)和黃河(觀測點:黃河勝利浮橋)現(xiàn)場調查,其它河流數(shù)據(jù)均引自歷史文獻資料(表1).考慮到長江并非主要輸入到黃海,按其入海徑流總量的14.1%[31]計算.

    1.2.3 大氣輸入(A) 黃渤海DIN和DIP的大氣干濕沉降依據(jù)黃?;厩Ю飵r島[32-33]降水和氣溶膠中和DIP的濃度,以及黃渤海降水量[34],干沉降速率[33]和黃渤海水域面積計算得到.

    1.2.4 黃海(Y)與東海(E)水交換 黃海與東海之間的水交換主要考慮東側北上的黃海暖流、黑潮支流和西側南下的蘇北沿岸流.根據(jù)本研究和東海DIN和DIP的相關資料(表1)以及水交換通量為基礎估算黃渤海和東海DIN和DIP的交換通量.

    1.2.5 沉積物—水界面釋放(BD) DIN和DIP在沉積物—水界面的擴散通量根據(jù)Fick第一擴散定律[35]計算:

    式中:F為擴散速率mol/(m2d),負號表示擴散是由沉積物向水體方向進行的;為沉積物孔隙率,引自文獻[36]相同或相近站位數(shù)據(jù),0.55~0.83;s為沉積物中DIN和DIP的擴散系數(shù);?/?為沉積物水界面處DIN和DIP的濃度梯度;0為無限稀釋溶液中溶質的分子擴散系數(shù)[37];為經(jīng)驗系數(shù),£0.7,=2;>0.7,=2.5~3.0[38].為了減小誤差,除了本研究計算得到的結果,DIN和DIP的擴散速率還引用了渤海[23]和黃海[26]的相關結果,以求覆蓋整個研究海域.

    1.2.6 氮和磷的凈埋藏(NB)與沉積量(S) 黃渤海有機氮(ON)和有機磷(OP)的沉積通量由其平均沉積速率和黃渤海面積所得,由于氮磷的快速再生[1,14],將計算結果視為凈沉積通量,即沉積量的一部分;其為沉積物中的過程,不參與水體N和P平衡過程計算.凈沉積速率計算公式如下:

    式中:N和P為沉積速率,mol/(m2·a);N和P為沉積物表層有機氮和有機磷含量,%;R為沉積速率,g/(m2·a);w為含水率,%,相關變量見表1.

    鑒于氮的快速再生,初級生產(chǎn)轉化的有機氮由水體向沉積物-水界面沉積量應為凈埋藏量與沉積物-水界面釋放量之和.磷的沉積量由Redfield比值根據(jù)氮的沉積量計算得出;該計算基于顆粒物在水深較淺的近海沉積所需時間也相對較短,水柱內物理化學過程對氮磷循環(huán)的速率差異化改造有限的假設[1].

    1.2.7 地下水輸入(GW) 黃渤海地下水DIN和DIP的輸入主要是通過226Ra、228Ra示蹤和端元模型估算地下水交換量[39-40],并根據(jù)地下水中DIN和DIP的含量[23,41]求出其通量.

    1.2.8 反硝化作用(Den) 反硝化作用作為水體N去除的重要途徑.采用如下經(jīng)驗公式[42]計算:

    Den=16.1×0.30(5)

    Den=Den×N(6)

    其中:Den為反硝化效率,%;為水體停留時間, month,以月計;Den為反硝化量,Gmol/a;N為外部輸入氮的總量,Gmol/a.該公式基于多個河口和近海水體的觀測結果回歸得到,可表征近岸水體反硝化的特征.另外,中國邊緣海的固氮作用并不強烈,多小于河流輸入和大氣沉降輸入[43],且應被海洋內部的反硝化作用(非外源部分)抵消,因此在本研究中可予以忽略.

    1.2.9 磷的吸附解吸與含磷礦物的風化作用(O) 與氮不同的是,磷與顆粒物之間存在強烈的界面作用[44-45],并對水體中磷的含量起到緩沖作用,如顆粒物中生物可利用磷等,其可在磷限制情況下提供生物生長所需要的磷;另外,輸入到海洋的含磷顆粒物以及海底磷酸鹽礦物也會對磷起到補充作用,其負荷采用收支計算的余值求出,正值表示移出,負值則表示添加.在水體氮相比于磷過剩的情況下,這一計算過程具有合理性.計算公式如下:

    O=R+GW+A+E+BD-S-Y(7)

    1.2.10 內部循環(huán)(Rc) 初級生產(chǎn)合成的有機顆粒在向沉積物界面輸送的過程中會有部分氮磷釋放出來,其應為初級生產(chǎn)合成的氮磷總量與沉積通量的差值;與沉積過程一樣,水柱內氮磷的再循環(huán)也遵循Redfield比值.

    1.2.11 水體氮的累積(Ac) 近年來黃渤海水體氮的濃度逐漸升高,這表明氮存在過剩和累積的現(xiàn)象,其負荷也用收支計算的余值給出,正值表示積累,負值則表示消耗.計算公式如下:

    Ac=R+GW+A+E+BD-S-Y-Den(8)

    2 結果與討論

    2.1 水體中營養(yǎng)鹽分布

    2.1.1 水體中DIN的分布特征 研究海域水體DIN的含量變化范圍為0.09~43.6μmol/L,平均值為8.86μmol/L.春季水體DIN的含量介于0.09~ 38.5μmol/L之間,平均值為5.62μmol/L;秋季DIN的含量介于0.21~43.6μmol/L之間,平均值為12.8μmol/L.硝酸鹽貢獻了77%的DIN,是水體DIN的主要形態(tài),亞硝酸鹽和銨鹽的比例分別為17%和6%.

    圖2為黃渤海DIN的分布.春季表層水體中DIN呈近岸高、離岸低的變化趨勢,高值區(qū)位于黃河口和長江口附近,渤海中部和黃海北部相對較低,表明陸源輸入對于DIN含量和分布的影響較大,尤其是在河口區(qū)域,DIN和DIP分布均與鹽度呈相反的梯度變化;底層水體中的DIN在河口區(qū)含量也較高,深水區(qū)DIN含量較沿岸區(qū)域高,這與表層的DIN分布略有不同,可能是受到底界面擴散和有機物降解釋放的影響.秋季水體中DIN的分布與春季相似,但在黃海中部底層出現(xiàn)DIN高值區(qū)域,這與黃海冷水團積累效應有關[46].

    2.1.2 水體中DIP的分布特征 黃渤海水體DIP的含量變化范圍為0.01~1.71μmol/L,平均值為0.24μmol/L.春季研究區(qū)域水體DIP的含量介于0.01~0.71μmol/L之間,平均值為0.13μmol/L;秋季DIP的含量介于0.01~1.71μmol/L之間,平均值為0.37μmol/L.

    由圖3可見,春季表底層的DIP分布同DIN的分布相似;秋季底層水體中的DIP在黃海中部海域較高,這是底界面擴散、生物擾動以及冷水團累積效應共同作用的結果.秋季表層DIP的含量在渤海、南黃海和長江口區(qū)域含量較高,長江沖淡水作用是長江口DIP出現(xiàn)高值的主要原因;渤海底層水體DIP含量呈沿岸高,離岸低的變化趨勢,尤其是萊州灣海域DIP高出其他海域數(shù)倍,這可能是因為底層物質再懸浮和礦化作用和海水養(yǎng)殖共同作用的結果.在真光層,浮游植物光合作用對真光層營養(yǎng)鹽的消耗和底層的營養(yǎng)鹽再生應是表層DIP含量低于底層的主要原因.

    2.2 沉積物間隙水中氮磷營養(yǎng)鹽的分布特征

    間隙水中氮營養(yǎng)鹽的垂直分布如圖4所示.間隙水中氨氮的含量遠遠高于硝酸鹽和亞硝酸鹽的含量,是間隙水中氮營養(yǎng)鹽的主要存在形態(tài),其濃度變化范圍在17~1080μmol/L,硝酸鹽和亞硝酸鹽的濃度變化分別為1~20μmol/L和0.1~15μmol/L,這主要是受微生物降解作用以及沉積物中氧化還原條件改變的影響.上覆水-沉積物界面DIN的濃度梯度表明調查站位的氨氮和亞硝酸鹽存在從沉積物向上覆水擴散的現(xiàn)象,其擴散速率分別為3.18mmol/(m2·d)和0.013mmol/(m2·d);DIN擴散速率為0.012mmol/ (m2·d).

    間隙水中DIP的垂直分布如圖4所示.各站位間隙水中DIP的空間分布差異顯著,其與沉積環(huán)境、沉積物類型和早期成巖礦化速率等密切相關.上覆水-沉積物界面DIP的濃度梯度顯示調查站位均存在DIP從沉積物向上覆海水擴散的現(xiàn)象,表明黃渤海沉積物是水體DIP的源,其平均擴散速率0.006mmol/(m2d),與Liu等[23,26]的結果相近.然而,間隙水DIP與DIN的比值(<0.01)遠小于Redfield比值[28],這表明磷酸鹽在沉積物-水界面存在吸附現(xiàn)象,并對磷由沉積物向水體釋放起到緩沖作用.

    2.3 黃渤海氮磷營養(yǎng)鹽的收支

    2.3.1 初級生產(chǎn)力 海洋中的DIN和DIP主要是通過浮游植物初級生產(chǎn)消耗.經(jīng)統(tǒng)計黃海和渤海真光層平均初級生產(chǎn)力分別為(112±47)g/ (m2·a)和(155±16)g/(m2·a)(以碳記),換算成黃海和渤海初級生產(chǎn)力所消耗的DIN分別為(1.42±0.59)mol/(m2·a)和(1.95±0.20)mol/(m2·a), DIP分別為(0.09±0.04)mol/(m2·a)和(0.12±0.01) mol/(m·a).由此估算得到黃渤海消耗的DIN和DIP的通量分別為(850±120)mol/(m2·a)和(53± 7.6)mol/(m2·a).

    3.3.2 河流輸入 黃渤海主要入海河流徑流量、DIN和DIP的含量及通量如表1所示.沿岸入海河流DIN和DIP的通量為(30±4.1)Gmol/a和(0.44±0.07)Gmol/a.對于DIN而言,入海的12條主要河流中黃河和長江的輸入量占總的河流輸入的15%和49%,鴨綠江、淮河和韓國的漢河(Han River)的輸入量則分別為10%、9%和9%,其余的小河流貢獻了剩余8%;長江(14.1%的徑流量進入黃海)、淮河和漢河是提供DIP的主要入海河流,分別為總的河流輸入量的51%、26%和15%.沿岸河流輸入到黃渤海的DIN和DIP的通量分別占中國入海河流總輸入量[47]的25%和15%.

    2.3.3 大氣輸入 根據(jù)表1,估算得到黃渤海的DIN的干濕沉降通量為62Gmol/a,黃渤海的DIP的干濕沉降通量為0.49Gmol/a.黃渤海DIN的大氣輸入的量約為河流的2倍,DIP的輸入的量與河流相當,這與黃海的結果相似[26].可見,大氣輸入是黃渤海氮磷營養(yǎng)鹽的重要來源.

    表1 黃渤海氮和磷的主要過程 Table 1 Main processes of N and P budget in the Bohai Sea and Yellow Sea

    續(xù)表1

    沉積埋藏海區(qū)沉積速率含水率ON含量OP含量ON沉積速率OP沉積速率ON沉積通量OP沉積通量 [g/(m2·a)](%)(%)(%)[mol/(m2·a)][mol/(m2·a)](Gmol/a)(Gmol/a) 渤海0.2~2.8 [58]41[59]0.02~0.080.01 [60]0.100.0037.40.23 黃海0.11~0.66 [61]33 [60]0.02~0.100.006~0.014 [62]0.150.004452.68 反硝化海區(qū)反硝化速率反硝化效率反硝化量 [kmol/(km2·a)](%)(Gmol/a) 渤海2436319 黃海1334350

    注:1墾利黃河浮橋2013~2015年數(shù)據(jù)(未發(fā)表數(shù)據(jù));2長江江陰大橋2013~2014年數(shù)據(jù)(未發(fā)表數(shù)據(jù));3韓國的河流;其他未標明的數(shù)據(jù)均來源于本研究或公式計算.考慮到系統(tǒng)誤差,各過程計算的最終結果保留1~2位有效數(shù)字.

    2.3.4 黃海與東海水交換 黃海西南部水體中DIN和DIP的濃度分別為3.5μmol/L和0.08μmol/L,東海東北部水體中DIN和DIP的濃度分別為5.65,0.05μmol/L.根據(jù)表1統(tǒng)計,得到黃海向東海輸出的DIN通量為25Gmol/a,DIP的通量為0.56Gmol/a;東海向黃海輸入的DIN通量為38Gmol/a,DIP的通量為0.33Gmol/a.由此可知,東海向黃海凈輸入的DIN的通量為13Gmol/a,黃海向東海凈輸出的DIP的通量為0.23Gmol/a.

    2.3.5 沉積物—水界面釋放 渤海和黃海沉積物—水界面的氨氮、硝酸鹽、亞硝酸鹽和磷酸鹽的擴散速率分別為2.81,-0.005,0.013, 0.0062mmol/(m2·d),再根據(jù)黃渤海的面積計算得到底界面DIN和DIP的擴散通量分別為350, 0.96Gmol/a.黃渤海相對較高的底界面DIN和DIP擴散通量對于維持水體較高的初級生產(chǎn)力具有重要的意義,同時也是底層水體DIN和DIP濃度高于表層水體的主要原因.

    2.3.6 地下水輸入 地下水的輸入對黃渤海氮磷的收支貢獻見表1,其DIN和DIP輸入通量分別為24,0.25Gmol/a.地下水氮磷營養(yǎng)鹽輸入通量與河流的貢獻在同一個數(shù)量級上,盡管黃渤海地下水的研究有限,但地下水輸入顯然是陸架邊緣海外部輸入氮磷的重要途徑[48],其對區(qū)域海洋物質循環(huán)的影響不可忽視.

    2.3.7 氮和磷的埋藏通量 由式3和4計算得到黃渤海氮和磷的埋藏通量分別為53, 2.9Gmol/a (表1).結合氮和磷的溶解擴散通量可知,黃渤海氮和磷的沉積通量約為400,25Gmol/a.有機氮和磷的凈埋藏效率分別為沉積通量的13%和12%.上述結果表明黃渤海是潛在的氮和磷的匯,對區(qū)域物質循環(huán)和碳的保存意義重大.

    2.3.8 反硝化作用 渤海和黃海水體停留時間分別為8年[22]和2.2年[26],根據(jù)式5計算可知,輸入到渤海和黃海的氮分別有63%和43%由反硝化作用被去除,這一結果與全球輸入到河口海灣的無機氮約50%通過反硝化作用從水體清除[49]的結論相符.根據(jù)上述結果計算得到黃渤海反硝化量為69Gmol/a,黃海和渤海的反硝化速率分別為240×103,130×103mol/(km2·a),這一計算結果與世界近岸陸架海域反硝化速率的平均值[610× 103mol/(km2·a)[42]]處于同一量級.

    2.3.9 磷的吸附解吸與含磷礦物的風化 根據(jù)黃渤海磷收支的分析,得出黃渤海磷的吸附解吸與含磷礦物的風化作用是水體初級生產(chǎn)所需磷的主要供應者,可向水體貢獻約23Gmol/a的活性磷酸鹽.黃渤海表層沉積物中生物可利用磷總量[60,62]與磷的吸附解吸的釋放量相當,表明解吸等過程具備向近海水體貢獻磷的基本的條件以及收支計算的合理性.同時,與Redfield比值[28]相比,沉積物間隙水中DIP含量(圖4)相對于DIN而言普遍較低,也表明沉積物-水界面磷的吸附對磷向水體釋放的影響.

    2.3.10 氮的累積 根據(jù)氮收支平衡(式8)大約為11Gmol/a的氮在水體累積;若以黃渤海水體體積(1810×1010m3,其中渤海平均水深18m,黃海平均水深44m),黃渤海水體11Gmol/a的氮盈余會提高水體氮的濃度,提高速率約為0.6μmol/(L·a),這與近10年來渤海氮濃度升高超5μmol/L的實測數(shù)據(jù)相符[22-23].

    2.3.11 收支計算 根據(jù)黃渤海氮收支的分析,得出黃渤海溶解無機氮收支各過程通量(圖5).對黃渤海而言,初級生產(chǎn)者消耗的DIN的總量為850Gmol/a,這與王保棟等[63]和田恬等[64]的計算相近;估算黃渤海內部循環(huán)的DIN的通量為450Gmol/a,內部循環(huán)的DIN貢獻了初級生產(chǎn)的53%.收支計算表明河流和大氣輸入,黃、東海水交換,地下水輸入以及底界面擴散釋放等過程提供的DIN的量不足以支撐黃渤海初級生產(chǎn)的量,大部分維持黃渤海初級生產(chǎn)的DIN和DIP更多的來源有顆粒態(tài)營養(yǎng)鹽的轉化和有機態(tài)營養(yǎng)鹽的降解.從氮的收支過程來看,底界面擴散是黃渤海水體DIN的主要來源,占總輸入的73%;其次是大氣輸入,占13%;河流、地下水和東海輸入的貢獻分別為6%、5%和3%;黃渤海水體DIN的支出主要是通過沉積埋藏和反硝化,其比例分別為85%和15%.黃渤海水體每年氮的積累占河流輸入的37%.

    由磷的收支過程得出黃渤海溶解無機磷的收支各過程通量見圖5.年初級生產(chǎn)是磷主要的匯,消耗的DIP的總量為53Gmol/a,該值與田恬等[64]的估算相近;黃渤海水體磷通過生物的吸收與隨后的沉積埋藏,其通量為25Gmol/a(磷依據(jù)Redfield比值[28]和氮的沉積通量求出);在黃渤海內部溶解無機磷循環(huán)平衡的前提下,我們估算黃渤海內部循環(huán)的DIP的通量為28Gmol/a,內部循環(huán)的DIP的量貢獻了初級生產(chǎn)的53%,這說明河流和大氣輸入,黃、東海水交換,地下水輸入以及底界面擴散釋放等過程提供的DIP的量不足以支撐黃渤海初級生產(chǎn)的量,其必然存在內部補充作用.在磷收支平衡的前提下,黃渤海水體DIP外源貢獻主要是磷的吸附解吸與含磷礦物的風化作用,其負荷為23Gmol/a,為水體獲取外部磷總量的91%,底界面擴散和大氣輸入,分別占總輸入的量分別為3.8%和1.9%,其次是河流輸入,占1.8%,地下水輸入貢獻為1.0%.DIP的支出主要是通過沉積埋藏和向東海的輸出,其比例分別為97%和3%.磷的埋藏效率((S-BD)/S)約為97%,這與江輝煌和劉素美[65]在渤海的研究結果相符.收支計算表明,黃渤海內部緩沖作用是磷獲得補充的主要方式,超過了河流等外部負荷,這在水體磷限制情況下具有合理性.

    2.4 N與P的歷史變化與長期生態(tài)效應

    由于人類活動影響,黃渤海的營養(yǎng)鹽結構自1980年代至今發(fā)生了很大的變化,黃渤海區(qū)域水體中DIN含量逐漸增加,DIP卻逐年降低;渤海的DIN由1985年的0.53μmol/L[22]上升到2000年5.12μmol/L再到2012年的10.6μmol/L;DIP則由1960年的0.8μmol/L[2]下降到2012年的0.3μmol/ L.同樣的黃海的DIN由1984年的0.79μmol/L[66]上升到2006~2007年的5.8~9.3μmol/L[24];DIP則由1960年的0.39μmol/L[66]下降到2012年的0.17μmol/L.水體氮磷濃度的相對變化表明磷的埋藏量在增加,而氮與磷非平衡的變化,必將對近海生態(tài)系統(tǒng)產(chǎn)生深遠的影響.

    更為重要的一點,由于周邊河流河水的氮磷比[13,31,47]明顯高于渤海和黃海的比值,而且DIN的通量也在逐漸增加[13,47],這對黃渤海營養(yǎng)鹽結構(特別是N/P嚴重比偏離Redfield比值)和區(qū)域海洋生態(tài)系統(tǒng)造成了很大的影響;如圖6所示,黃渤海近幾十年來的赤潮爆發(fā)的頻率逐年增加,年赤潮累積發(fā)生面積也在逐漸擴大(圖6a和6b),同時渤海初級生產(chǎn)相對于1990年代也有較大幅度的提高,這些很大程度上是由于入海氮營養(yǎng)鹽通量的增加(圖6c,與河流磷負荷無顯著相關關系)提高了近海浮游植物生產(chǎn)力.在收支計算中,沉積物埋藏的氮和磷表明近海沉積物匯的作用將影響包括氮磷在內的物質循環(huán)的總量.沉積物中凈埋藏的N/P比(18)遠小于其它外源負荷輸入的N/P比(44~130),這意味著沉積物埋藏磷的相對量比氮多,即表現(xiàn)為水體DIP濃度的逐漸降低,并會進一步加劇磷限制的趨勢,同時也會對近海碳的埋藏產(chǎn)生重要的影響.

    陸地氮負荷的增加提高了海洋初級生產(chǎn)力水平(圖6c).若目前的初級生產(chǎn)力水平得以維持一段時間,磷由于初級生產(chǎn)力水平的提高而產(chǎn)生的高的埋藏量會導致其水體濃度進一步的降低,并使得黃渤海水體由潛在的比值磷限制[28]向絕對濃度上的磷限制轉變(即<0.1μmol/L)[74],這對于生態(tài)系統(tǒng)而言將是災難性的,尤其在氮累積量逐年增加的當前.因此,我們應該關注黃渤海氮磷營養(yǎng)鹽的結構變化及生態(tài)環(huán)境效應,并進行陸海統(tǒng)籌,逐步減少陸地氮的排放,以維持海洋生態(tài)系統(tǒng)的健康發(fā)展.

    3 結論

    3.1 春季黃渤海水體DIN和DIP的分布有明顯的差異性,高值基本都出現(xiàn)在近岸和河口區(qū)域.河流輸入對于其含量與分布影響較為明顯.沉積物間隙水中DIN和DIP濃度顯著高于上覆水體,上覆水-沉積物界面存在DIN和DIP從沉積物向上覆水擴散的現(xiàn)象.

    3.2 黃渤海氮的收支評估表明,底界面擴散是黃渤海水體DIN的主要來源,占總輸入的73%;其次是大氣輸入,占13%;河流、地下水和東海輸入分別貢獻了6%、5%和3%;黃渤海水體DIN的支出主要是通過沉積埋藏和反硝化,其比例分別為85%和15%.

    3.3 黃渤海水體DIP的來源主要是磷的吸附解吸與含磷礦物的風化作用,其貢獻為91%,底界面擴散和大氣輸入,分別占總輸入的量分別為3.8%和1.9%;其次是河流輸入,占1.8%;地下水輸入貢獻為1.0%.DIP的支出主要是通過沉積埋藏和向東海的輸出,其比例分別為97%和3%.

    3.4 黃渤海每年有11Gmol/a的氮在水體積累,并導致氮濃度提高約0.6μmol/(L·a).近幾十年來陸源氮負荷持續(xù)增加,提高了磷的埋藏率,并致使黃渤海的營養(yǎng)鹽結構發(fā)生了很大的變化,使得水體中氮含量逐漸增加,磷含量逐漸減少,N/P遠大于Redfield值,赤潮發(fā)生頻次和累計面積也有較大幅度的上升.我們應該進行陸海統(tǒng)籌,逐步減少陸地氮的排放,以維持海洋生態(tài)系統(tǒng)的結構和功能.

    Tyrrell, T. The relative influences of nitrogen and phosphorus on oceanic primary production [J]. Nature, 1999,400:525–531.

    Tang Q, Jin X, Wang J, et al. Decadal-scale variations of ecosystem productivity and control mechanisms in the Bohai Sea [J]. Fisheries Oceanography, 2003,12(4/5):223-233.

    高 嵩,范士亮,韓秀榮,等.滸苔綠潮與南黃海近岸海域水質的關系 [J]. 中國環(huán)境科學, 2014,34(1):213-218.

    Diaz R J, Rosenberg R. Spreading dead zones and consequences for marine ecosystems [J]. Science, 2008,321(5891):926-929.

    Wei G F, Tang D L, Wang S. Distribution of chlorophyll and harmful algal blooms (HABs): A review on space based studies in the coastal environments of Chinese marginal seas [J]. Advances in Space Research, 2008,41(1):12-19.

    Wu Z, Yu Z, Song X, et al. The spatial and temporal characteristics of harmful algal blooms in the southwest Bohai Sea [J]. Continental Shelf Research, 2013,59:10-17.

    Rousseaux C S, Gregg W W. Recent decadal trends in global phytoplankton composition [J]. Global Biogeochemical Cycles, 2015,29(10):1674-1688.

    Bauer J E, Cai W J, Raymond P A, et al. The changing carbon cycle of the coastal ocean [J]. Nature, 2013,504:61-70.

    Berner R A. Burial of organic carbon and pyrite sulfur in the modern ocean: its geochemical and environmental significance [J]. American Journal of Science, 1982,282(4):451-473.

    Hedges J I, Keil R G. Sedimentary organic matter preservation: an assessment and speculative synthesis [J]. Marine chemistry, 1995,49(2):81-115.

    Mantoura R F C, Martin J M, Wollast R. Ocean margin processes in global change [M]. John Wiley and Sons Ltd, 1991.

    Muller-Karger F E, Varela R, Thunell R, et al. The importance of continental margins in the global carbon cycle [J]. Geophysical Research Letters, 2005,32(1), L01602.

    Dai Z, Du J, Zhang X, et al. Variation of Riverine Material Loads and Environmental Consequences on the Changjiang (Yangtze) Estuary in Recent Decades (1955?2008) [J]. Environmental Science & Technology, 2010,45(1):223-227.

    Andrews J E, Burgess D, Cave R R, et al. Biogeochemical value of managed realignment, Humber estuary, UK [J]. Science of the Total Environment, 2006,371:19–30.

    Smith S V, Swaney D P, Buddemeier R W, et al. River nutrient loads and catchment size [J]. Biogeochemistry, 2005,75(1):83- 107.

    曾呈奎,鴻 儒,春 林.中國海洋志 [M]. 鄭州:大象出版社, 2003.

    鞏 瑤.黃河下游利津站營養(yǎng)鹽輸送規(guī)律及影響因素研究 [D]. 青島:中國海洋大學, 2012,75-80.

    Liu S M. Response of nutrient transports to water–sediment regulation events in the Huanghe basin and its impact on the biogeochemistry of the Bohai [J]. Journal of Marine Systems, 2015,141:59-70.

    Ran X B, Che H, Zang J Y, et al. Variability in the composition and export of silica in the Huanghe River Basin [J]. Science China: Earth Sciences, 2015,58(11):2078-2089.

    Wang H, Yang Z, Saito Y, et al. Stepwise decreases of the Huanghe (Yellow River) sediment load (1950–2005): Impacts of climate change and human activities [J]. Global and Planetary Change, 2007,57(3):331-354.

    王 婷.2002~2004年及調水調沙期間黃河下游營養(yǎng)鹽的變化特征 [D]. 青島:中國海洋大學, 2007:14-59.

    Ning X, Lin C, Su J, et al. Long-term environmental changes and the responses of the ecosystems in the Bohai Sea during 1960~ 1996 [J]. Deep Sea Research Part II: Topical Studies in Oceanography, 2010,57(11):1079-1091.

    Liu S M, Li L W, Zhang Z. Inventory of nutrients in the Bohai [J]. Continental Shelf Research, 2011,31(16):1790-1797.

    Wei Q, Yao Q, Wang B, et al. Long-term variation of nutrients in the southern Yellow Sea [J]. Continental Shelf Research, 2015, 111:184-196.

    孫 萍,李瑞香,李 艷,等.2005年夏末渤海網(wǎng)采浮游植物群落結構 [J]. 海洋科學進展, 2008,26(3):354-363.

    Liu S M, Zhang J, Chen S Z, et al. Inventory of nutrient compounds in the Yellow Sea [J]. Continental Shelf Research, 2003,23(11):1161-1174.

    劉 軍,臧家業(yè),張麗君,等.黃海硅的分布與收支研究 [J]. 中國環(huán)境科學, 2016,36(1):157-166.

    Redfield A C. The biological control of chemical factors in the environment [J]. American scientist, 1958:205-221.

    朱明遠,毛興華,呂瑞華,等.黃海海區(qū)的葉綠素a和初級生產(chǎn)力 [J]. 黃渤海海洋, 1993,11(3):38-51.

    檀賽春,石廣玉.中國近海初級生產(chǎn)力的遙感研究及其時空演化 [J]. 地理學報, 2007,61(11):1189-1199.

    Liu S M, Zhang J, Chen H T, et al. Nutrients in the Changjiang and its tributaries [J]. Biogeochemistry, 2003,62(1):1-18.

    Zhang J, Yu Z G, Raabe T, et al. Dynamics of inorganic nutrient species in the Bohai seawaters [J]. Journal of Marine Systems, 2004,44(3):189-212.

    Zhang G, Zhang J, Liu S. Characterization of nutrients in the atmospheric wet and dry deposition observed at the two monitoring sites over Yellow Sea and East China Sea [J]. Journal of Atmospheric Chemistry, 2007,57(1):41-57.

    Martin J M, Zhang J, Shi M C, et al. Actual flux of the Huanghe (Yellow River) sediment to the western Pacific Ocean [J]. Netherlands Journal of Sea Research, 1993,31(3):243-254.

    Berner R A. Early diagenesis: A theoretical approach [M]. Princeton University Press, Princeton, 1980:241-242.

    李玲偉.沉積物—水界面交換和黃河輸入對渤海營養(yǎng)鹽的影響 [D]. 青島:中國海洋大學, 2010:47-55.

    Li Y H, Gregory S. Diffusion of ions in sea water and in deep-sea sediments [J]. Geochimica et Cosmochimica Acta, 1974,38(5): 703-714.

    Ullman W J, Aller R C. Diffusion coefficients in nearshore marine sediments [J]. Limnology and Oceanography, 1982,27(3):552- 556.

    Kim G, Ryu J W, Yang H S, et al. Submarine groundwater discharge (SGD) into the Yellow Sea revealed by228Ra and226Ra isotopes: Implications for global silicate fluxes [J]. Earth and Planetary Science Letters, 2005,237(1):156-166.

    Peterson R N, Burnett W C, Taniguchi M, et al. Radon and radium isotope assessment of submarine groundwater discharge in the Yellow River delta, China [J]. Journal of Geophysical Research, 2008,113,C09021.

    Waska H, Kim G. Submarine groundwater discharge (SGD) as a main nutrient source for benthic and water-column primary production in a large intertidal environment of the Yellow Sea [J]. Journal of Sea Research, 2011,65(1):103-113.

    Seitzinger S, Harrison J A, B?hlke J K, et al. Denitrification across landscapes and waterscapes: a synthesis [J]. Ecological Applications, 2006,16(6):2064-2090.

    張 潤.中國邊緣海生物固氮作用的研究 [D]. 廈門:廈門大學, 2011.

    Slomp C P, Malschaert J F P, Van Raaphorst W. The role of adsorption in sediment‐water exchange of phosphate in North Sea continental margin sediments [J]. Limnology and Oceanography, 1998,43(5):832-846.

    Slomp C P, Van Cappellen P. The global marine phosphorus cycle: sensitivity to oceanic circulation [J]. Biogeosciences, 2007,4: 155-171.

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

    Liu S M, Hong G H, Zhang J, et al. Nutrient budgets for large Chinese estuaries [J]. Biogeosciences, 2009,6(10):2245-2263.

    Moore W S. The Effect of Submarine Groundwater Discharge on the Ocean [J]. Annual Review of Marine Science, 2010,2:59-88.

    Seitzinger S P, Kroeze C. Global distribution of nitrous oxide production and N inputs in freshwater and coastal marine ecosystems [J]. Global biogeochemical cycles, 1998,12(1):93- 113.

    中華人民共和國水利部.中國河流泥沙公報2007-2013 [DB/OL]. http://www.mwr.gov.cn/.

    水利部海河水利委員會.海河流域水資源公報2005年[EB/OL]. [2014-12-01] http://www.hwcc.gov.cn/pub/hwcc/static/szygb/ gongbao2005/index.htm.

    卞少偉,孫 韌,梅鵬蔚,等.2013年春夏季天津近岸海域浮游植物的群落結構 [J]. 水生態(tài)學雜志, 2015,36(1):47-52.

    斌.2005年夏季環(huán)渤海16條主要河流的污染狀況及入海通量 [D]. 青島:中國海洋大學, 2007:13-16.

    劉敬偉,李富祥,劉 月,等.鴨綠江下游徑流時序變化特征研究 [J]. 人民黃河, 2011,33(10):34-36.

    Kim J K, Jung S, Eom J, et al. Dissolved and particulate organic carbon concentrations in stream water and relationships with land use in multiple-use watersheds of the Han River (Korea) [J]. Water International, 2013,38(3):326-339.

    楊守業(yè),Jung Hoi-soo,李從先,等.黃河,長江與韓國Keum, Yeongsan江沉積物常量元素地球化學特征 [J]. 地球化學, 2004,33(1):99-105.

    米鐵柱,姚慶禎,孟 佳,等.年春、夏季黃海、東海營養(yǎng)鹽分布特征的研究 [J]. 海洋與湖沼, 2012,43(3):678-688.

    董太祿.渤?,F(xiàn)代沉積作用與模式的研究 [J]. 海洋地質與第四紀地質, 1996,16(4):43-53.

    盧 博,李趕先,黃韶健,等.中國黃海,東海和南海北部海底淺層沉積物聲學物理性質之比較 [J]. 海洋技術, 2005,24(2):28-33.

    趙志梅.渤海灣沉積物磷形態(tài)及營養(yǎng)鹽在沉積物-水界面交換的研究 [D]. 楊凌:西北農(nóng)林科技大學, 2005,27-28.

    李鳳業(yè),高 抒,賈建軍,等,黃、渤海泥質沉積區(qū)現(xiàn)代沉積速率 [J]. 海洋與湖沼, 2002,33(4):364-369.

    宋國棟.黃東海沉積物中磷的形態(tài)及其影響因素研究 [D]. 青島:中國海洋大學, 2014,22-25.

    王保棟,單寶田,戰(zhàn) 閆.黃、渤海無機氮的收支模式初探 [J]. 海洋科學, 2002,(2):33-36.

    田 恬,魏 皓,蘇 健,等.黃海氮磷營養(yǎng)鹽的循環(huán)和收支研究 [J]. 海洋科學進展, 2003,21(1):1-11.

    江輝煌,劉素美.渤海沉積物中磷的分布與埋藏通量 [J]. 環(huán)境科學學報, 2013,33(1):125-132.

    Li H M, Zhang C S, Han X R, et al. Changes in concentrations of oxygen, dissolved nitrogen, phosphate, and silicate in the southern Yellow Sea, 1980–2012: Sources and seaward gradients [J]. Estuarine, Coastal and Shelf Science, 2014,163(20):44-55.

    國家海洋局.中國海洋災害公報1995~2011 [EB/OL]. http://www.soa.gov.cn/zwgk/hygb/zghyzhgb/.

    Li Z Y, Bai J, Shi J H, et al. Distributions of inorganic nutrients in the Bohai Sea of China [J]. Journal of Ocean University of Qingdao, 2003,2(1):112-116.

    孫 軍,劉東艷,徐 俊,等.1999年春季渤海中部及其鄰近海域的網(wǎng)采浮游植物群落 [J]. 生態(tài)學報, 2004(9):179-192.

    苑明莉,孫 軍,翟惟東.2012年秋季渤海和北黃海浮游植物群落 [J]. 天津科技大學學報, 2014,29(6):56-64.

    張志鋒,賀 欣,張 哲,等.渤海富營養(yǎng)化現(xiàn)狀,機制及其與赤潮的時空耦合性 [J]. 海洋環(huán)境科學, 2012,31(4):465-468.

    Wang X, Cui Z, Guo Q, et al. Distribution of nutrients and eutrophication assessment in the Bohai Sea of China [J]. Chinese Journal of Oceanology and Limnology, 2009,27:177-183.

    孫 軍,劉東艷,柴心玉,等.1998~1999年春秋季渤海中部及其鄰近海域葉綠素a濃度及初級生產(chǎn)力估算 [J]. 生態(tài)學報, 2003,23(3):517-526.

    Nelson D M, Brzezinski M A. Kinetics of silicate acid uptake by natural diatom assemblages in two Gulf and Stream warm-core rings [J]. Marine Ecology Progress Series, 1990,62:283-292.

    致謝:本研究的現(xiàn)場采樣工作由國家海洋局第一海洋研究所車宏和鄭莉莉等研究生協(xié)助完成,在此表示感謝.

    * 責任作者, 副研究員, rxb@fio.org.cn

    Distribution and budget of nitrogen and phosphorus and their influence on the ecosystem in the Bohai Sea and Yellow Sea

    ZHAO Chen-ying1, ZANG Jia-ye1, LIU Jun1,2, SUN Tao1, RAN Xiang-bin1*

    (1.Research Center for Marine Ecology, First Institute of Oceanography, State Oceanic Administration, Qingdao 266061, China;2.Key Laboratory of Marine Chemistry Theory and Technology, Ministry of Education, Ocean University of China, Qingdao 266100, China)., 2016,36(7):2115~2127

    Based on results obtained from the comprehensive investigation in the Bohai Sea and the Yellow Sea (BYS), distributions and sources of dissolved inorganic nitrogen (DIN) and dissolved inorganic phosphorus (DIP) in the sea water and sediment pore water in the study area were discussed; and a budget model of nitrogen (N) and phosphorus (P) in the BYS was established. Concentrations of DIN and DIP in the study area are affected by riverine input and benthic flux. The major sources of the DIN in the water column of the BYS are benthic diffusion, followed by the atmospheric deposition and inputs from adjacent rivers the groundwater and the East China Sea; while the dominant removals of DIN from water column in the BYS are phytoplankton utilization and sedimentation. The key sources of the DIP are phosphorus minerals desorption and weathering, accounting for 91%, followed by benthic flux and atmospheric deposition, then riverine flux and groundwater input; while the dominant outputs of DIP are sedimentation and the outflow to the East China Sea. There was a 11Gmol/a of nitrogen accumulated in the BYS, which would result in an increasing of DIN concentration by 0.6 μmol/(L·a). DIN concentrations and N/P ratio in the Bohai Sea and the Yellow Sea BYS kept increasing in recent decades due to the increasing of nitrogen N accumulation and phosphorus P sedimentation. Terrestrial nitrogen loadings might enhance the primary production and alter the structure and function of ecosystem, which had would have significant effects on the ecosystem of Bohai Sea and Yellow Sea the BYS in the long-time scale. We, therefore, should reduce the riverine nitrogen N loads to maintain the ecosystem stability in the BYS Bohai Sea and Yellow Sea.

    Bohai Sea and Yellow Sea;dissolved inorganic nitrogen;dissolved inorganic phosphorus phosphate;flux and budget

    X55

    A

    1000-6923(2016)07-2115-13

    趙晨英(1991- ),女,陜西寶雞人,國家海洋局第一海洋研究所碩士研究生,主要從事海洋生物地球化學研究.發(fā)表論文1篇.

    2015-12-07

    國家自然科學基金(41376093,41106072)

    猜你喜歡
    黃渤海營養(yǎng)鹽黃海
    黃渤海不同水層中浮游植物對灰霾添加的響應
    東方濕地 黃海明珠
    華人時刊(2020年21期)2021-01-14 01:32:28
    黃海簡介
    出海
    涼水河子河營養(yǎng)鹽滯留能力評估
    休漁補網(wǎng)忙
    三角恒等變換去哪兒了
    黃渤海區(qū)域衛(wèi)星云圖反演海面溫度設計初探夏光濱
    瓊東海域冬季、夏季營養(yǎng)鹽結構特征及其對浮游植物生長的影響
    2012年冬季南海西北部營養(yǎng)鹽分布及結構特征
    亚洲avbb在线观看| 中文字幕人妻丝袜一区二区| 一个人看视频在线观看www免费 | 午夜精品一区二区三区免费看| 真人一进一出gif抽搐免费| 亚洲乱码一区二区免费版| 国产免费男女视频| 久久久久性生活片| 两个人视频免费观看高清| 国产精品一及| 精品一区二区三区视频在线观看免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 9191精品国产免费久久| 色综合欧美亚洲国产小说| 亚洲国产精品999在线| 搞女人的毛片| 色av中文字幕| 18美女黄网站色大片免费观看| 99热只有精品国产| 色在线成人网| 国产精品乱码一区二三区的特点| 老司机福利观看| 亚洲欧美日韩高清专用| 欧美日韩中文字幕国产精品一区二区三区| 18禁黄网站禁片免费观看直播| 午夜两性在线视频| 国产精品亚洲美女久久久| 美女 人体艺术 gogo| 两个人视频免费观看高清| 国产私拍福利视频在线观看| 午夜福利在线观看免费完整高清在 | 9191精品国产免费久久| 亚洲av电影不卡..在线观看| a在线观看视频网站| 草草在线视频免费看| 国产精品国产高清国产av| 日本在线视频免费播放| 我要搜黄色片| 成人国产一区最新在线观看| 亚洲成av人片在线播放无| 国产亚洲精品综合一区在线观看| 一级作爱视频免费观看| 偷拍熟女少妇极品色| 色在线成人网| 成年女人看的毛片在线观看| 亚洲人成网站在线播| 精品欧美国产一区二区三| 免费人成视频x8x8入口观看| 观看美女的网站| 亚洲精品一区av在线观看| 日韩大尺度精品在线看网址| 免费搜索国产男女视频| 久久久久久人人人人人| 日本熟妇午夜| 亚洲欧美日韩高清在线视频| 女警被强在线播放| 性色av乱码一区二区三区2| 国产欧美日韩一区二区三| 午夜福利在线在线| 成人国产综合亚洲| 亚洲人成电影免费在线| 88av欧美| 成人国产综合亚洲| 中文字幕人妻熟人妻熟丝袜美 | 久久久久精品国产欧美久久久| 午夜福利在线观看吧| 天堂动漫精品| 国产爱豆传媒在线观看| 日韩欧美在线二视频| 男女做爰动态图高潮gif福利片| 波多野结衣巨乳人妻| 精品熟女少妇八av免费久了| 悠悠久久av| 在线国产一区二区在线| 99热6这里只有精品| 国产精品精品国产色婷婷| www日本黄色视频网| 欧美不卡视频在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 久久精品综合一区二区三区| 成人三级黄色视频| 最好的美女福利视频网| 久久久久久九九精品二区国产| 一区二区三区激情视频| 又紧又爽又黄一区二区| 国内精品一区二区在线观看| 韩国av一区二区三区四区| 窝窝影院91人妻| 国产在线精品亚洲第一网站| 亚洲天堂国产精品一区在线| 欧洲精品卡2卡3卡4卡5卡区| 欧美成人免费av一区二区三区| 亚洲成a人片在线一区二区| 欧美日韩一级在线毛片| 国产精品嫩草影院av在线观看 | 内地一区二区视频在线| 又爽又黄无遮挡网站| 国产日本99.免费观看| 很黄的视频免费| 成人性生交大片免费视频hd| 禁无遮挡网站| 男人的好看免费观看在线视频| 一边摸一边抽搐一进一小说| 亚洲五月天丁香| 久久亚洲真实| 天天躁日日操中文字幕| 69av精品久久久久久| 亚洲七黄色美女视频| 亚洲av一区综合| 俄罗斯特黄特色一大片| 偷拍熟女少妇极品色| 国产精品自产拍在线观看55亚洲| 3wmmmm亚洲av在线观看| 亚洲一区二区三区色噜噜| 成人永久免费在线观看视频| 国产又黄又爽又无遮挡在线| 97超级碰碰碰精品色视频在线观看| 国产精品 国内视频| 国产精品精品国产色婷婷| 草草在线视频免费看| 久久久久亚洲av毛片大全| www日本在线高清视频| 成人鲁丝片一二三区免费| 免费观看的影片在线观看| 久久久精品欧美日韩精品| 国产主播在线观看一区二区| 精品人妻偷拍中文字幕| 伊人久久精品亚洲午夜| 久久久久久久亚洲中文字幕 | 免费av毛片视频| 亚洲真实伦在线观看| 亚洲专区国产一区二区| 看片在线看免费视频| 成人18禁在线播放| 非洲黑人性xxxx精品又粗又长| 18美女黄网站色大片免费观看| 丰满乱子伦码专区| 国产三级中文精品| 亚洲人成网站在线播放欧美日韩| 亚洲国产欧美网| 久久久久精品国产欧美久久久| 无人区码免费观看不卡| 日本 欧美在线| 国产精品自产拍在线观看55亚洲| 国产私拍福利视频在线观看| av中文乱码字幕在线| 国产精品一区二区三区四区免费观看 | 国产探花在线观看一区二区| 两人在一起打扑克的视频| 男人的好看免费观看在线视频| 偷拍熟女少妇极品色| 变态另类丝袜制服| 高清日韩中文字幕在线| 亚洲国产精品999在线| 麻豆久久精品国产亚洲av| 波野结衣二区三区在线 | 丰满人妻熟妇乱又伦精品不卡| 日本撒尿小便嘘嘘汇集6| 亚洲欧美日韩高清专用| 在线观看午夜福利视频| 给我免费播放毛片高清在线观看| 亚洲欧美日韩卡通动漫| 久久中文看片网| www日本黄色视频网| 亚洲 国产 在线| 国产精品1区2区在线观看.| 亚洲中文日韩欧美视频| 久久久久久久亚洲中文字幕 | 欧美三级亚洲精品| 又黄又粗又硬又大视频| 男女那种视频在线观看| 99久国产av精品| 成年免费大片在线观看| 国产亚洲av嫩草精品影院| 狠狠狠狠99中文字幕| 精品久久久久久久毛片微露脸| 深夜精品福利| 日韩 欧美 亚洲 中文字幕| 女人高潮潮喷娇喘18禁视频| 12—13女人毛片做爰片一| 免费看光身美女| 脱女人内裤的视频| 一本精品99久久精品77| 国产亚洲av嫩草精品影院| 69av精品久久久久久| 亚洲av电影在线进入| 他把我摸到了高潮在线观看| av国产免费在线观看| 国产精品综合久久久久久久免费| av福利片在线观看| 757午夜福利合集在线观看| 91九色精品人成在线观看| 免费在线观看日本一区| 99热这里只有是精品50| 欧美成狂野欧美在线观看| 亚洲av成人av| 男女做爰动态图高潮gif福利片| 午夜影院日韩av| 不卡一级毛片| 成人国产综合亚洲| 国产激情偷乱视频一区二区| 国产精品久久电影中文字幕| 日本撒尿小便嘘嘘汇集6| tocl精华| 俄罗斯特黄特色一大片| 亚洲熟妇熟女久久| 高清日韩中文字幕在线| 精品国产三级普通话版| 欧美av亚洲av综合av国产av| 日韩精品青青久久久久久| ponron亚洲| 国内精品美女久久久久久| 欧美日韩综合久久久久久 | 少妇裸体淫交视频免费看高清| 少妇的逼好多水| 国产精品,欧美在线| 最新中文字幕久久久久| 午夜免费成人在线视频| 国产精品综合久久久久久久免费| 丁香欧美五月| 久久久久久九九精品二区国产| 亚洲精品一卡2卡三卡4卡5卡| 一个人看视频在线观看www免费 | 怎么达到女性高潮| 中文资源天堂在线| 搡老熟女国产l中国老女人| 成人一区二区视频在线观看| 一本精品99久久精品77| 国产亚洲精品久久久com| 亚洲国产欧美人成| 少妇的逼好多水| 日韩 欧美 亚洲 中文字幕| 国内少妇人妻偷人精品xxx网站| 成人高潮视频无遮挡免费网站| 国产在线精品亚洲第一网站| 国产成人啪精品午夜网站| 97超级碰碰碰精品色视频在线观看| 欧美黑人巨大hd| 99久久无色码亚洲精品果冻| 日本a在线网址| 日本一二三区视频观看| aaaaa片日本免费| 国产日本99.免费观看| 长腿黑丝高跟| 我要搜黄色片| 国产日本99.免费观看| 精品久久久久久久人妻蜜臀av| 国产老妇女一区| 欧美+亚洲+日韩+国产| 人人妻人人澡欧美一区二区| 精品人妻偷拍中文字幕| 老司机在亚洲福利影院| 一个人看视频在线观看www免费 | 国产伦在线观看视频一区| 国产老妇女一区| 久久久久亚洲av毛片大全| 精品国产美女av久久久久小说| 天天添夜夜摸| 69av精品久久久久久| 免费在线观看亚洲国产| 亚洲不卡免费看| 免费一级毛片在线播放高清视频| 天堂av国产一区二区熟女人妻| 欧美日韩国产亚洲二区| 成人永久免费在线观看视频| 日韩国内少妇激情av| 757午夜福利合集在线观看| 国产精品98久久久久久宅男小说| 亚洲成人免费电影在线观看| 亚洲 国产 在线| 热99在线观看视频| 国产一区在线观看成人免费| 网址你懂的国产日韩在线| 亚洲av五月六月丁香网| 国语自产精品视频在线第100页| 亚洲av二区三区四区| 欧美日韩瑟瑟在线播放| 欧美最新免费一区二区三区| 熟女电影av网| 国产午夜福利久久久久久| 中文字幕av成人在线电影| 中文欧美无线码| 亚洲av成人精品一二三区| 免费看光身美女| 亚洲国产精品国产精品| 啦啦啦啦在线视频资源| 成人亚洲精品av一区二区| 蜜桃亚洲精品一区二区三区| 一级av片app| 午夜爱爱视频在线播放| 国产亚洲一区二区精品| 女的被弄到高潮叫床怎么办| 欧美成人午夜免费资源| 亚洲在线自拍视频| 成人综合一区亚洲| 久久这里有精品视频免费| 亚洲欧美成人精品一区二区| 亚洲在久久综合| 人妻制服诱惑在线中文字幕| 大陆偷拍与自拍| 身体一侧抽搐| 午夜福利高清视频| 亚洲精品日韩在线中文字幕| 人人妻人人澡欧美一区二区| 免费播放大片免费观看视频在线观看| 国产成人91sexporn| 亚洲成人中文字幕在线播放| 99热网站在线观看| 2018国产大陆天天弄谢| 国产精品精品国产色婷婷| 欧美性感艳星| 午夜福利在线观看吧| 日韩亚洲欧美综合| a级毛色黄片| 97精品久久久久久久久久精品| 成人鲁丝片一二三区免费| 男人爽女人下面视频在线观看| 国产69精品久久久久777片| 又爽又黄无遮挡网站| 久久精品国产亚洲av天美| 日本一本二区三区精品| 日韩在线高清观看一区二区三区| 国产精品一区二区三区四区免费观看| 国产熟女欧美一区二区| 伦理电影大哥的女人| 午夜福利在线观看免费完整高清在| 成人亚洲精品一区在线观看 | 中文欧美无线码| 久久这里只有精品中国| 99热网站在线观看| 国产一区二区亚洲精品在线观看| 观看美女的网站| 精品久久久精品久久久| 国产黄片视频在线免费观看| av女优亚洲男人天堂| 欧美潮喷喷水| 中文字幕av成人在线电影| 一级毛片电影观看| 国产成人aa在线观看| 赤兔流量卡办理| 2022亚洲国产成人精品| 国产精品精品国产色婷婷| 久久精品久久久久久久性| 亚洲国产日韩欧美精品在线观看| 久久久午夜欧美精品| 嫩草影院精品99| 日韩一区二区三区影片| 麻豆成人av视频| 日本欧美国产在线视频| 联通29元200g的流量卡| 成人毛片a级毛片在线播放| 国产一区二区亚洲精品在线观看| 最近的中文字幕免费完整| 在线播放无遮挡| av在线观看视频网站免费| 午夜亚洲福利在线播放| 欧美成人精品欧美一级黄| 免费无遮挡裸体视频| 国产成人91sexporn| 日本一二三区视频观看| 夜夜爽夜夜爽视频| 男女啪啪激烈高潮av片| 国产亚洲最大av| 成人二区视频| 亚洲色图av天堂| 一级毛片 在线播放| 国产老妇伦熟女老妇高清| 成人一区二区视频在线观看| 日韩精品青青久久久久久| 青青草视频在线视频观看| 久久久久久久国产电影| 大陆偷拍与自拍| 亚洲熟妇中文字幕五十中出| 日韩 亚洲 欧美在线| 亚洲性久久影院| 欧美日韩精品成人综合77777| 大片免费播放器 马上看| 国产精品麻豆人妻色哟哟久久 | 观看美女的网站| 欧美成人午夜免费资源| 亚洲电影在线观看av| 国产精品日韩av在线免费观看| 麻豆av噜噜一区二区三区| 少妇高潮的动态图| 高清日韩中文字幕在线| 国产淫片久久久久久久久| 国产亚洲精品av在线| 欧美bdsm另类| 禁无遮挡网站| 中文字幕av在线有码专区| 亚洲精品日韩av片在线观看| 天天躁日日操中文字幕| 成人美女网站在线观看视频| 日韩精品青青久久久久久| 亚洲精品久久午夜乱码| 亚洲精品乱码久久久v下载方式| 天堂中文最新版在线下载 | 亚洲欧美精品专区久久| 毛片一级片免费看久久久久| 老师上课跳d突然被开到最大视频| 边亲边吃奶的免费视频| 国产综合精华液| 国产乱人偷精品视频| 边亲边吃奶的免费视频| 国产黄a三级三级三级人| 免费看av在线观看网站| 老司机影院毛片| 亚洲欧洲日产国产| 精品人妻熟女av久视频| 日韩一区二区视频免费看| 日韩欧美精品v在线| 性色avwww在线观看| 亚洲人成网站在线观看播放| 国产 一区精品| 婷婷色综合www| 亚洲国产日韩欧美精品在线观看| 日韩精品有码人妻一区| 国产老妇伦熟女老妇高清| 欧美人与善性xxx| 色吧在线观看| 麻豆乱淫一区二区| 黄色日韩在线| 色综合站精品国产| 熟妇人妻不卡中文字幕| 成人毛片60女人毛片免费| 偷拍熟女少妇极品色| 天堂影院成人在线观看| av网站免费在线观看视频 | 午夜福利成人在线免费观看| 国产免费福利视频在线观看| 国产精品女同一区二区软件| 国产在线一区二区三区精| 91久久精品电影网| 成人毛片a级毛片在线播放| 亚洲最大成人av| 欧美极品一区二区三区四区| 午夜福利网站1000一区二区三区| 校园人妻丝袜中文字幕| 精品久久久噜噜| 免费观看a级毛片全部| 久久久色成人| 天天一区二区日本电影三级| 国产男人的电影天堂91| 干丝袜人妻中文字幕| 成年女人看的毛片在线观看| 超碰av人人做人人爽久久| 在线观看人妻少妇| 自拍偷自拍亚洲精品老妇| 高清欧美精品videossex| 欧美人与善性xxx| 国产精品一区二区三区四区久久| 日本黄色片子视频| 国产午夜福利久久久久久| 国产黄片视频在线免费观看| 成人综合一区亚洲| 久久精品国产自在天天线| 国产片特级美女逼逼视频| 欧美激情久久久久久爽电影| 午夜激情久久久久久久| 成年av动漫网址| 日日啪夜夜撸| 一区二区三区乱码不卡18| 欧美日本视频| 五月玫瑰六月丁香| 欧美3d第一页| 日韩欧美 国产精品| 亚洲婷婷狠狠爱综合网| 午夜免费男女啪啪视频观看| 日韩制服骚丝袜av| 久久精品久久久久久久性| 你懂的网址亚洲精品在线观看| 一级毛片久久久久久久久女| 国产综合懂色| 久久热精品热| 天天躁日日操中文字幕| 国产精品不卡视频一区二区| 国产精品人妻久久久影院| 精品久久久久久成人av| freevideosex欧美| 亚洲av免费高清在线观看| 精品酒店卫生间| 亚洲av电影在线观看一区二区三区 | 高清视频免费观看一区二区 | 亚洲天堂国产精品一区在线| 亚洲乱码一区二区免费版| 欧美xxxx黑人xx丫x性爽| 久久久久久久久中文| 国产 亚洲一区二区三区 | 在线免费观看不下载黄p国产| 色综合色国产| 成人特级av手机在线观看| 高清在线视频一区二区三区| 日产精品乱码卡一卡2卡三| 最后的刺客免费高清国语| 久久久久免费精品人妻一区二区| 美女国产视频在线观看| 18禁动态无遮挡网站| 国产男女超爽视频在线观看| 久久久久久久午夜电影| 国产在线一区二区三区精| 国产黄色免费在线视频| 久久久久网色| 一级毛片黄色毛片免费观看视频| 丝袜美腿在线中文| 久久久久免费精品人妻一区二区| 高清av免费在线| 十八禁国产超污无遮挡网站| 日韩伦理黄色片| 国产69精品久久久久777片| 国产乱人视频| 好男人在线观看高清免费视频| 国产白丝娇喘喷水9色精品| 一级毛片电影观看| 中文字幕人妻熟人妻熟丝袜美| 男人舔奶头视频| 性色avwww在线观看| 搞女人的毛片| 乱人视频在线观看| 午夜免费观看性视频| 亚州av有码| 日日啪夜夜爽| 美女主播在线视频| 亚洲一区高清亚洲精品| 免费观看av网站的网址| 午夜亚洲福利在线播放| 久久久精品免费免费高清| 一边亲一边摸免费视频| 亚洲av成人精品一区久久| 丰满人妻一区二区三区视频av| 午夜久久久久精精品| 免费av毛片视频| 99久久人妻综合| 免费播放大片免费观看视频在线观看| 欧美成人精品欧美一级黄| 哪个播放器可以免费观看大片| 婷婷色av中文字幕| 亚洲欧美精品专区久久| 午夜爱爱视频在线播放| 免费无遮挡裸体视频| 欧美另类一区| 男女那种视频在线观看| 97人妻精品一区二区三区麻豆| 乱系列少妇在线播放| 亚洲国产高清在线一区二区三| 丝袜美腿在线中文| 少妇高潮的动态图| 久久精品久久精品一区二区三区| 男女那种视频在线观看| 国产 一区精品| 一级av片app| 日韩欧美 国产精品| 日日摸夜夜添夜夜爱| 国产成人一区二区在线| av在线天堂中文字幕| 亚洲国产日韩欧美精品在线观看| or卡值多少钱| 日日摸夜夜添夜夜添av毛片| 青春草亚洲视频在线观看| 国产又色又爽无遮挡免| 一边亲一边摸免费视频| 久久久a久久爽久久v久久| 少妇丰满av| 国产亚洲一区二区精品| 男女边摸边吃奶| 国产麻豆成人av免费视频| 成人毛片60女人毛片免费| 在现免费观看毛片| 午夜福利在线观看免费完整高清在| av女优亚洲男人天堂| 国产爱豆传媒在线观看| 亚洲精品乱码久久久久久按摩| 五月天丁香电影| 成人亚洲欧美一区二区av| 欧美成人一区二区免费高清观看| 69人妻影院| 女人久久www免费人成看片| 久久午夜福利片| 熟女电影av网| 国产亚洲av嫩草精品影院| 免费高清在线观看视频在线观看| 成人鲁丝片一二三区免费| av免费在线看不卡| 欧美变态另类bdsm刘玥| 久久精品综合一区二区三区| 18+在线观看网站| 一区二区三区四区激情视频| 欧美激情国产日韩精品一区| 啦啦啦中文免费视频观看日本| 夜夜看夜夜爽夜夜摸| 国产视频内射| 国产成人免费观看mmmm| 国产欧美日韩精品一区二区| 18禁在线播放成人免费| 欧美不卡视频在线免费观看| 亚洲高清免费不卡视频| 国产精品一区二区性色av| 国产视频内射| 99热网站在线观看| 中文字幕av在线有码专区| 美女xxoo啪啪120秒动态图| 九草在线视频观看| 午夜福利在线观看吧| 亚洲欧美一区二区三区黑人 | 日韩成人av中文字幕在线观看| 大话2 男鬼变身卡| 久久精品国产亚洲av天美| 国产综合精华液| 午夜福利高清视频| 最近最新中文字幕大全电影3| 久久97久久精品| 久久久午夜欧美精品| 成人综合一区亚洲| 国产精品久久视频播放| 91av网一区二区| 成人午夜精彩视频在线观看|