• <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)鹽分布及結構特征
    国产麻豆成人av免费视频| 一级黄色大片毛片| 婷婷色综合大香蕉| 超碰av人人做人人爽久久| 给我免费播放毛片高清在线观看| 久久久久久久亚洲中文字幕 | 国产三级黄色录像| 午夜精品久久久久久毛片777| 老女人水多毛片| 白带黄色成豆腐渣| 久久精品久久久久久噜噜老黄 | 成年人黄色毛片网站| 99久国产av精品| 免费人成在线观看视频色| 精品国产亚洲在线| 五月玫瑰六月丁香| 午夜老司机福利剧场| 欧美日韩亚洲国产一区二区在线观看| 老司机深夜福利视频在线观看| 日韩欧美国产一区二区入口| 国产精品爽爽va在线观看网站| 最近在线观看免费完整版| 精品久久久久久成人av| 国内精品一区二区在线观看| 久久精品久久久久久噜噜老黄 | 亚洲第一欧美日韩一区二区三区| 人妻制服诱惑在线中文字幕| 色在线成人网| 免费观看人在逋| 国产中年淑女户外野战色| 国产探花在线观看一区二区| 国产精品一区二区性色av| 91麻豆精品激情在线观看国产| 免费人成视频x8x8入口观看| 久久精品国产亚洲av天美| 日本熟妇午夜| 色噜噜av男人的天堂激情| 免费看美女性在线毛片视频| 欧美3d第一页| 欧美+日韩+精品| 午夜免费成人在线视频| 听说在线观看完整版免费高清| 国产美女午夜福利| 亚洲成av人片在线播放无| ponron亚洲| 国产欧美日韩精品亚洲av| 午夜免费激情av| 高清在线国产一区| 欧美三级亚洲精品| 国产一区二区在线av高清观看| 亚洲最大成人手机在线| 国产伦精品一区二区三区四那| 欧美日韩黄片免| 夜夜躁狠狠躁天天躁| 亚洲专区国产一区二区| 丰满人妻一区二区三区视频av| 久久午夜福利片| 国产v大片淫在线免费观看| 成人欧美大片| 婷婷精品国产亚洲av在线| 久久精品综合一区二区三区| 99热这里只有是精品在线观看 | 日日干狠狠操夜夜爽| 欧美黄色淫秽网站| 日韩国内少妇激情av| 一级a爱片免费观看的视频| 美女xxoo啪啪120秒动态图 | 宅男免费午夜| 亚洲五月婷婷丁香| a级一级毛片免费在线观看| 少妇高潮的动态图| 欧美日韩国产亚洲二区| 国产成人a区在线观看| 怎么达到女性高潮| 啦啦啦观看免费观看视频高清| 尤物成人国产欧美一区二区三区| 99久久精品国产亚洲精品| 精品午夜福利视频在线观看一区| 国产综合懂色| 91九色精品人成在线观看| 国产精品久久久久久久电影| 日韩欧美国产一区二区入口| 国产蜜桃级精品一区二区三区| 色综合站精品国产| 毛片一级片免费看久久久久 | 久久6这里有精品| 精品熟女少妇八av免费久了| 日本在线视频免费播放| 在线观看免费视频日本深夜| 亚洲激情在线av| 国产成人影院久久av| 精品久久久久久久久久免费视频| 亚洲 国产 在线| 久久久久久大精品| 亚洲久久久久久中文字幕| 少妇高潮的动态图| 免费高清视频大片| 亚洲人与动物交配视频| 免费看光身美女| 精品国产三级普通话版| 午夜精品在线福利| 91在线观看av| 亚洲片人在线观看| 日韩亚洲欧美综合| 国产熟女xx| 亚洲国产精品sss在线观看| 嫩草影院新地址| 精品日产1卡2卡| 亚洲色图av天堂| 久久久久久久久大av| 亚洲最大成人手机在线| 国产午夜精品论理片| 亚洲国产精品sss在线观看| 亚洲aⅴ乱码一区二区在线播放| 精品久久久久久成人av| 国产一区二区三区视频了| 天堂动漫精品| 美女高潮的动态| 伊人久久精品亚洲午夜| 久久中文看片网| 在现免费观看毛片| 成人永久免费在线观看视频| 免费av毛片视频| 国产亚洲av嫩草精品影院| 啦啦啦观看免费观看视频高清| 此物有八面人人有两片| 免费人成视频x8x8入口观看| 日本一本二区三区精品| 亚洲精品一区av在线观看| 精品免费久久久久久久清纯| 色尼玛亚洲综合影院| 熟女电影av网| 一个人免费在线观看的高清视频| 欧美另类亚洲清纯唯美| 亚洲一区二区三区色噜噜| 精品久久久久久久久av| 国产精品久久久久久久久免 | 极品教师在线免费播放| 日本在线视频免费播放| 亚洲 国产 在线| 男女之事视频高清在线观看| 观看免费一级毛片| 国内少妇人妻偷人精品xxx网站| 黄色视频,在线免费观看| 欧美日韩亚洲国产一区二区在线观看| av在线蜜桃| 亚洲av日韩精品久久久久久密| 国产视频内射| 婷婷精品国产亚洲av| 美女高潮喷水抽搐中文字幕| 国产麻豆成人av免费视频| 熟妇人妻久久中文字幕3abv| 亚洲18禁久久av| 国产成人av教育| 国产免费男女视频| 亚洲av免费在线观看| 成人性生交大片免费视频hd| 人妻制服诱惑在线中文字幕| 床上黄色一级片| 久久久国产成人精品二区| .国产精品久久| 变态另类丝袜制服| 欧美3d第一页| 成人特级黄色片久久久久久久| 亚洲内射少妇av| 97超级碰碰碰精品色视频在线观看| 少妇熟女aⅴ在线视频| 丰满的人妻完整版| 久9热在线精品视频| 99久久99久久久精品蜜桃| 午夜a级毛片| 久久性视频一级片| 熟妇人妻久久中文字幕3abv| 波多野结衣高清作品| 中文字幕精品亚洲无线码一区| 欧美绝顶高潮抽搐喷水| 深夜精品福利| 男人狂女人下面高潮的视频| 亚洲精品在线美女| 18禁在线播放成人免费| 亚洲人成网站高清观看| 无遮挡黄片免费观看| 天天躁日日操中文字幕| 国产伦精品一区二区三区四那| 精品久久久久久久久亚洲 | 乱码一卡2卡4卡精品| 国产高清视频在线播放一区| 亚洲七黄色美女视频| 亚洲人成网站高清观看| 亚洲av电影不卡..在线观看| 久久精品91蜜桃| 成人特级av手机在线观看| 色综合欧美亚洲国产小说| 亚洲国产高清在线一区二区三| 成人特级av手机在线观看| 精品久久久久久成人av| 国内少妇人妻偷人精品xxx网站| 国产高清有码在线观看视频| 天堂网av新在线| 亚洲精品一区av在线观看| 在线观看午夜福利视频| 国产精品亚洲美女久久久| 国产精品,欧美在线| 少妇人妻精品综合一区二区 | 十八禁网站免费在线| 欧美色欧美亚洲另类二区| 一个人看的www免费观看视频| 亚洲精品粉嫩美女一区| 亚洲第一电影网av| 午夜福利18| 久久性视频一级片| 国语自产精品视频在线第100页| 国产精品一区二区三区四区久久| 亚洲人成网站在线播| 人妻夜夜爽99麻豆av| 免费观看人在逋| 亚洲在线自拍视频| 日本 av在线| 别揉我奶头~嗯~啊~动态视频| 久久精品国产亚洲av涩爱 | 人人妻人人看人人澡| 最新中文字幕久久久久| 亚洲人成伊人成综合网2020| av在线蜜桃| 亚洲欧美日韩卡通动漫| 亚洲中文日韩欧美视频| 欧美精品国产亚洲| 亚洲人成网站在线播| 日韩人妻高清精品专区| 欧美乱妇无乱码| 一个人看视频在线观看www免费| 亚洲国产欧洲综合997久久,| av黄色大香蕉| 婷婷色综合大香蕉| 亚洲人与动物交配视频| 身体一侧抽搐| 高清毛片免费观看视频网站| 757午夜福利合集在线观看| 国产精品三级大全| 国产成人啪精品午夜网站| 波多野结衣高清无吗| 欧美xxxx黑人xx丫x性爽| 中文字幕高清在线视频| 亚洲18禁久久av| 校园春色视频在线观看| 自拍偷自拍亚洲精品老妇| 高清毛片免费观看视频网站| 国产探花极品一区二区| 久久久久久大精品| 亚洲三级黄色毛片| 久久精品国产清高在天天线| 色综合欧美亚洲国产小说| 日本免费一区二区三区高清不卡| 亚洲成av人片免费观看| 亚洲av二区三区四区| 久久久久久国产a免费观看| 欧洲精品卡2卡3卡4卡5卡区| 色在线成人网| 在线看三级毛片| 亚洲中文字幕一区二区三区有码在线看| 观看美女的网站| 免费观看精品视频网站| 高清日韩中文字幕在线| 极品教师在线视频| 免费看a级黄色片| 亚洲最大成人手机在线| 国产亚洲精品久久久久久毛片| 美女xxoo啪啪120秒动态图 | av视频在线观看入口| 午夜视频国产福利| www.999成人在线观看| 欧洲精品卡2卡3卡4卡5卡区| 3wmmmm亚洲av在线观看| 国产单亲对白刺激| 欧美日本亚洲视频在线播放| 国产人妻一区二区三区在| 亚洲精品一卡2卡三卡4卡5卡| 国内毛片毛片毛片毛片毛片| 在线播放无遮挡| 亚洲av免费在线观看| 免费黄网站久久成人精品 | 日本在线视频免费播放| 在线观看一区二区三区| 亚洲av不卡在线观看| 91午夜精品亚洲一区二区三区 | 高清在线国产一区| 久久亚洲精品不卡| 婷婷丁香在线五月| 亚洲精品日韩av片在线观看| 97超级碰碰碰精品色视频在线观看| 两个人的视频大全免费| 日本一二三区视频观看| 一级黄色大片毛片| 中国美女看黄片| 成人高潮视频无遮挡免费网站| 国产精品日韩av在线免费观看| 国产在视频线在精品| 久久国产乱子免费精品| 男人狂女人下面高潮的视频| 亚洲美女黄片视频| 亚洲欧美激情综合另类| 好男人电影高清在线观看| 桃红色精品国产亚洲av| 在线a可以看的网站| 制服丝袜大香蕉在线| 日日摸夜夜添夜夜添小说| 免费av毛片视频| 欧美日韩亚洲国产一区二区在线观看| 3wmmmm亚洲av在线观看| 精品日产1卡2卡| 2021天堂中文幕一二区在线观| 麻豆国产av国片精品| av天堂中文字幕网| 男女视频在线观看网站免费| 国产伦一二天堂av在线观看| 亚洲性夜色夜夜综合| 亚洲av中文字字幕乱码综合| 久久久久久久午夜电影| 中文字幕久久专区| 亚洲熟妇熟女久久| 国产色婷婷99| 国产成年人精品一区二区| 亚洲国产精品久久男人天堂| 在线观看一区二区三区| 老司机福利观看| 午夜福利在线在线| x7x7x7水蜜桃| 久久久久久九九精品二区国产| 欧美xxxx黑人xx丫x性爽| 黄色一级大片看看| 深爱激情五月婷婷| 亚洲无线观看免费| 免费看美女性在线毛片视频| 熟妇人妻久久中文字幕3abv| 日本a在线网址| 欧美日本视频| a级毛片免费高清观看在线播放| 精品一区二区三区av网在线观看| 国产精品一区二区三区四区久久| 亚洲三级黄色毛片| 精品久久国产蜜桃| 亚洲中文日韩欧美视频| 舔av片在线| 国产欧美日韩精品一区二区| 99在线人妻在线中文字幕| 淫妇啪啪啪对白视频| 97超视频在线观看视频| 又黄又爽又刺激的免费视频.| 欧美午夜高清在线| 国内久久婷婷六月综合欲色啪| 亚洲精品成人久久久久久| 99热这里只有是精品在线观看 | 别揉我奶头 嗯啊视频| 美女免费视频网站| 精品日产1卡2卡| 精华霜和精华液先用哪个| 两人在一起打扑克的视频| 午夜免费成人在线视频| 国产精品亚洲av一区麻豆| 精品熟女少妇八av免费久了| 每晚都被弄得嗷嗷叫到高潮| 国产真实乱freesex| 两个人的视频大全免费| 狂野欧美白嫩少妇大欣赏| 一区二区三区激情视频| 99久久久亚洲精品蜜臀av| 亚洲狠狠婷婷综合久久图片| 欧美成人免费av一区二区三区| 亚洲 国产 在线| 欧美一区二区亚洲| 中文字幕av在线有码专区| 有码 亚洲区| 国产欧美日韩精品一区二区| 国产精品爽爽va在线观看网站| 成人精品一区二区免费| 国产黄a三级三级三级人| 天堂网av新在线| 国产精品一区二区三区四区久久| 美女被艹到高潮喷水动态| 麻豆国产97在线/欧美| 日日摸夜夜添夜夜添小说| 女生性感内裤真人,穿戴方法视频| av黄色大香蕉| 五月伊人婷婷丁香| 老司机午夜十八禁免费视频| 夜夜看夜夜爽夜夜摸| 国产精品人妻久久久久久| 亚洲 欧美 日韩 在线 免费| 露出奶头的视频| 1024手机看黄色片| 男女那种视频在线观看| 日韩中字成人| 日本一二三区视频观看| 少妇熟女aⅴ在线视频| 少妇的逼好多水| 麻豆国产av国片精品| 伊人久久精品亚洲午夜| 亚州av有码| 久久精品国产99精品国产亚洲性色| 亚洲第一区二区三区不卡| 亚洲av电影在线进入| www.熟女人妻精品国产| 欧美+亚洲+日韩+国产| 久久天躁狠狠躁夜夜2o2o| 赤兔流量卡办理| 国产免费av片在线观看野外av| 亚洲人成电影免费在线| 亚洲色图av天堂| 国产一区二区三区在线臀色熟女| 乱码一卡2卡4卡精品| 亚洲欧美日韩无卡精品| 99精品久久久久人妻精品| 国语自产精品视频在线第100页| 我要看日韩黄色一级片| 国产成人a区在线观看| 国产高清激情床上av| 精品无人区乱码1区二区| 全区人妻精品视频| 性色av乱码一区二区三区2| 久久久久久久久久成人| 性插视频无遮挡在线免费观看| 国产亚洲精品久久久久久毛片| 男女视频在线观看网站免费| 一个人看视频在线观看www免费| 亚洲狠狠婷婷综合久久图片| 色在线成人网| 亚洲精品久久国产高清桃花| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品99久久久久久久久| 国产精品免费一区二区三区在线| 真人一进一出gif抽搐免费| 欧美日韩瑟瑟在线播放| 欧美黑人欧美精品刺激| 99热这里只有精品一区| 久久久久亚洲av毛片大全| 亚洲欧美清纯卡通| 熟妇人妻久久中文字幕3abv| 97超视频在线观看视频| 成人亚洲精品av一区二区| av天堂在线播放| 久久性视频一级片| 成年版毛片免费区| 九九热线精品视视频播放| 成人欧美大片| 我的女老师完整版在线观看| 久久久精品大字幕| 一个人免费在线观看电影| 伦理电影大哥的女人| 一区二区三区免费毛片| 亚洲一区二区三区色噜噜| 少妇人妻精品综合一区二区 | 成人av在线播放网站| 日韩人妻高清精品专区| 久久草成人影院| 久久99热这里只有精品18| 国产 一区 欧美 日韩| 有码 亚洲区| 国产精品亚洲美女久久久| 一个人观看的视频www高清免费观看| 亚洲中文字幕日韩| 中文字幕精品亚洲无线码一区| 97人妻精品一区二区三区麻豆| 丰满人妻熟妇乱又伦精品不卡| 亚洲熟妇中文字幕五十中出| 中文资源天堂在线| 男女视频在线观看网站免费| 三级国产精品欧美在线观看| 亚洲人成伊人成综合网2020| 国产不卡一卡二| 观看免费一级毛片| a级毛片免费高清观看在线播放| 久久久久久久久久成人| 国模一区二区三区四区视频| 身体一侧抽搐| 国产69精品久久久久777片| 五月伊人婷婷丁香| 国产高清三级在线| 亚洲狠狠婷婷综合久久图片| 成年版毛片免费区| 亚洲成av人片免费观看| 国产av在哪里看| 国产精品久久电影中文字幕| 欧美精品国产亚洲| 成人精品一区二区免费| 亚洲国产欧美人成| 亚洲av成人精品一区久久| 黄片小视频在线播放| 亚洲天堂国产精品一区在线| 国产一区二区在线观看日韩| 99久久成人亚洲精品观看| 亚洲av成人av| 日日摸夜夜添夜夜添av毛片 | 国产v大片淫在线免费观看| 高清日韩中文字幕在线| 国产精品人妻久久久久久| 久久亚洲真实| 在线观看av片永久免费下载| 成人亚洲精品av一区二区| 日本精品一区二区三区蜜桃| 国产亚洲精品av在线| 日日干狠狠操夜夜爽| 精品一区二区三区视频在线观看免费| x7x7x7水蜜桃| 日韩欧美一区二区三区在线观看| 中文字幕高清在线视频| 嫩草影院入口| 成人av一区二区三区在线看| 一级av片app| 一区福利在线观看| 欧美国产日韩亚洲一区| av欧美777| 日本三级黄在线观看| 我的女老师完整版在线观看| 久久久久久久亚洲中文字幕 | 日本成人三级电影网站| 中文字幕精品亚洲无线码一区| 亚洲片人在线观看| 亚洲av中文字字幕乱码综合| 毛片一级片免费看久久久久 | 3wmmmm亚洲av在线观看| 免费在线观看日本一区| 性色av乱码一区二区三区2| 精品国产亚洲在线| 一级av片app| 国产成年人精品一区二区| 国产伦一二天堂av在线观看| 禁无遮挡网站| 亚洲av第一区精品v没综合| av福利片在线观看| 精品无人区乱码1区二区| 九九久久精品国产亚洲av麻豆| 欧美高清性xxxxhd video| 男女做爰动态图高潮gif福利片| 别揉我奶头 嗯啊视频| 久久精品国产清高在天天线| 女人被狂操c到高潮| 亚洲第一区二区三区不卡| 狠狠狠狠99中文字幕| 国产午夜精品论理片| 久久久久久国产a免费观看| 免费观看精品视频网站| 床上黄色一级片| 日韩大尺度精品在线看网址| 日本精品一区二区三区蜜桃| 99精品在免费线老司机午夜| 九九热线精品视视频播放| 好看av亚洲va欧美ⅴa在| 日韩欧美免费精品| 免费看光身美女| 免费人成视频x8x8入口观看| 三级男女做爰猛烈吃奶摸视频| 久久久久九九精品影院| 熟女电影av网| 欧美色视频一区免费| 免费在线观看影片大全网站| 麻豆国产av国片精品| 91午夜精品亚洲一区二区三区 | 精品99又大又爽又粗少妇毛片 | 天堂√8在线中文| 中文字幕高清在线视频| 超碰av人人做人人爽久久| 国产精品国产高清国产av| 日韩高清综合在线| 精品一区二区三区视频在线观看免费| 精品人妻熟女av久视频| 99久久久亚洲精品蜜臀av| 桃色一区二区三区在线观看| 男女下面进入的视频免费午夜| 亚洲人成网站高清观看| 欧美一区二区国产精品久久精品| 在线十欧美十亚洲十日本专区| 99精品久久久久人妻精品| 久久这里只有精品中国| 久久国产精品人妻蜜桃| а√天堂www在线а√下载| av在线蜜桃| 精品久久久久久成人av| or卡值多少钱| 全区人妻精品视频| 精品午夜福利视频在线观看一区| 久久久久久国产a免费观看| 又紧又爽又黄一区二区| xxxwww97欧美| 国产高清视频在线观看网站| 白带黄色成豆腐渣| 欧美色欧美亚洲另类二区| 五月玫瑰六月丁香| 99久久九九国产精品国产免费| 欧美日韩福利视频一区二区| 3wmmmm亚洲av在线观看| 中文字幕av在线有码专区| 嫁个100分男人电影在线观看| 18禁黄网站禁片免费观看直播| 免费一级毛片在线播放高清视频| 久久久久国内视频| 一级a爱片免费观看的视频| 免费一级毛片在线播放高清视频| 性色avwww在线观看| 一级a爱片免费观看的视频| 亚洲av中文字字幕乱码综合| 十八禁人妻一区二区| 精华霜和精华液先用哪个| 麻豆一二三区av精品| 亚洲精品久久国产高清桃花| 久久国产乱子免费精品| 真人做人爱边吃奶动态| 成人无遮挡网站| 无遮挡黄片免费观看| 国产精品野战在线观看| 久久久久九九精品影院| 欧美bdsm另类| 人人妻人人看人人澡| 在现免费观看毛片| 九九久久精品国产亚洲av麻豆| 久久久久久久久久成人|