• <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)鹽分布及結構特征
    真人做人爱边吃奶动态| 黄色视频不卡| 成人午夜精彩视频在线观看| 乱人伦中国视频| 美女视频免费永久观看网站| a级毛片黄视频| 又黄又粗又硬又大视频| 国产老妇伦熟女老妇高清| 久久综合国产亚洲精品| 在线观看免费高清a一片| 国产成人av激情在线播放| 黑人欧美特级aaaaaa片| 9色porny在线观看| 久久久久久久精品精品| 各种免费的搞黄视频| 国产黄色免费在线视频| 婷婷色综合大香蕉| cao死你这个sao货| 国产亚洲精品久久久久5区| 99re6热这里在线精品视频| 国产不卡av网站在线观看| 啦啦啦啦在线视频资源| 国产精品一区二区在线不卡| 少妇人妻久久综合中文| 又大又爽又粗| 99香蕉大伊视频| 国产黄色免费在线视频| 天天操日日干夜夜撸| 久久久国产精品麻豆| 啦啦啦 在线观看视频| 午夜福利,免费看| 国产淫语在线视频| 美女中出高潮动态图| 亚洲精品一卡2卡三卡4卡5卡 | 成人午夜精彩视频在线观看| 九草在线视频观看| 成年美女黄网站色视频大全免费| 亚洲av欧美aⅴ国产| av福利片在线| av网站在线播放免费| 欧美少妇被猛烈插入视频| 91精品国产国语对白视频| 午夜福利影视在线免费观看| 97精品久久久久久久久久精品| 韩国高清视频一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 每晚都被弄得嗷嗷叫到高潮| 黑人巨大精品欧美一区二区蜜桃| 男人添女人高潮全过程视频| 欧美变态另类bdsm刘玥| 欧美老熟妇乱子伦牲交| 精品一品国产午夜福利视频| 亚洲精品美女久久久久99蜜臀 | 日日摸夜夜添夜夜爱| 久久 成人 亚洲| 国产精品 欧美亚洲| 99久久精品国产亚洲精品| 啦啦啦 在线观看视频| 国产日韩欧美视频二区| 久久人人爽人人片av| 精品少妇内射三级| 国产免费视频播放在线视频| 亚洲国产最新在线播放| 亚洲av电影在线观看一区二区三区| √禁漫天堂资源中文www| 天天躁夜夜躁狠狠久久av| 国产精品一二三区在线看| 一二三四在线观看免费中文在| 免费看av在线观看网站| 亚洲国产看品久久| 国产深夜福利视频在线观看| 精品人妻1区二区| 婷婷色av中文字幕| 在线观看www视频免费| 一级黄色大片毛片| 狠狠精品人妻久久久久久综合| 丝袜美足系列| 日韩中文字幕视频在线看片| 国产精品 欧美亚洲| 精品少妇久久久久久888优播| 在线观看免费日韩欧美大片| 丝袜人妻中文字幕| 国产女主播在线喷水免费视频网站| 一级片免费观看大全| 亚洲av欧美aⅴ国产| 蜜桃在线观看..| 亚洲第一青青草原| 亚洲视频免费观看视频| 91成人精品电影| 日韩电影二区| 免费看不卡的av| 叶爱在线成人免费视频播放| 国产精品免费大片| 亚洲成人免费电影在线观看 | 女警被强在线播放| 丰满少妇做爰视频| avwww免费| 国产精品九九99| 熟女av电影| 亚洲国产精品一区二区三区在线| 大码成人一级视频| 亚洲精品一卡2卡三卡4卡5卡 | 欧美黄色片欧美黄色片| 日本午夜av视频| 在线观看国产h片| 一级黄色大片毛片| 九草在线视频观看| 精品国产一区二区久久| 中文字幕精品免费在线观看视频| 男人添女人高潮全过程视频| 1024视频免费在线观看| 国产精品三级大全| 亚洲情色 制服丝袜| 欧美+亚洲+日韩+国产| tube8黄色片| 亚洲国产毛片av蜜桃av| 欧美成人午夜精品| 91字幕亚洲| 久久精品国产a三级三级三级| 国产一区亚洲一区在线观看| 亚洲国产av新网站| 黄色片一级片一级黄色片| 操美女的视频在线观看| 婷婷色综合大香蕉| 亚洲,欧美精品.| 亚洲精品一区蜜桃| 国产精品久久久久久精品古装| 老司机影院成人| 秋霞在线观看毛片| 亚洲精品在线美女| 欧美日韩精品网址| 久久精品久久久久久久性| 激情视频va一区二区三区| 亚洲熟女毛片儿| 亚洲少妇的诱惑av| 亚洲久久久国产精品| 满18在线观看网站| 欧美人与性动交α欧美精品济南到| 国产亚洲欧美精品永久| 亚洲精品成人av观看孕妇| 久久亚洲国产成人精品v| 一区在线观看完整版| 国产免费视频播放在线视频| 精品亚洲成国产av| 中国美女看黄片| 黄色怎么调成土黄色| 亚洲 国产 在线| 丝袜美腿诱惑在线| 午夜久久久在线观看| 亚洲三区欧美一区| 国产av精品麻豆| 精品久久蜜臀av无| 一区二区三区激情视频| av有码第一页| 在线av久久热| 狂野欧美激情性bbbbbb| 九草在线视频观看| 一边亲一边摸免费视频| 久久国产精品影院| av在线app专区| 色婷婷久久久亚洲欧美| 最新在线观看一区二区三区 | 美女大奶头黄色视频| 人成视频在线观看免费观看| 欧美激情极品国产一区二区三区| 国产深夜福利视频在线观看| 汤姆久久久久久久影院中文字幕| 国产精品麻豆人妻色哟哟久久| av福利片在线| 男人操女人黄网站| 国产精品秋霞免费鲁丝片| 婷婷色综合大香蕉| 精品国产乱码久久久久久小说| 电影成人av| 久久免费观看电影| 欧美亚洲 丝袜 人妻 在线| 久久久精品国产亚洲av高清涩受| 中国国产av一级| 一本久久精品| 免费在线观看视频国产中文字幕亚洲 | 性色av乱码一区二区三区2| 中文欧美无线码| 日本色播在线视频| 国产视频首页在线观看| 91字幕亚洲| 国产精品二区激情视频| 高清视频免费观看一区二区| 18禁国产床啪视频网站| 亚洲七黄色美女视频| 欧美日韩亚洲国产一区二区在线观看 | 男女午夜视频在线观看| 国产精品久久久久久精品古装| 熟女av电影| 亚洲五月婷婷丁香| 777米奇影视久久| 精品国产国语对白av| 晚上一个人看的免费电影| 亚洲欧美色中文字幕在线| 十八禁网站网址无遮挡| 精品免费久久久久久久清纯 | 久久精品成人免费网站| 十八禁人妻一区二区| 成人国产av品久久久| 人妻 亚洲 视频| 黄片小视频在线播放| 大话2 男鬼变身卡| 精品熟女少妇八av免费久了| 欧美精品av麻豆av| 韩国精品一区二区三区| 亚洲激情五月婷婷啪啪| 一级毛片我不卡| 亚洲第一av免费看| 操出白浆在线播放| 一边摸一边抽搐一进一出视频| 亚洲,欧美精品.| 极品人妻少妇av视频| 手机成人av网站| 99re6热这里在线精品视频| 欧美黄色淫秽网站| 亚洲欧美激情在线| 男女免费视频国产| 99久久综合免费| 欧美黑人精品巨大| 十分钟在线观看高清视频www| 秋霞在线观看毛片| 精品国产乱码久久久久久男人| 亚洲欧洲日产国产| 新久久久久国产一级毛片| 老鸭窝网址在线观看| 色视频在线一区二区三区| 新久久久久国产一级毛片| 中国美女看黄片| 校园人妻丝袜中文字幕| 国产淫语在线视频| 国产成人av教育| 国产精品免费视频内射| 90打野战视频偷拍视频| 久久综合国产亚洲精品| 精品一品国产午夜福利视频| 亚洲国产av影院在线观看| 美女国产高潮福利片在线看| 国产欧美亚洲国产| 欧美少妇被猛烈插入视频| 18禁观看日本| 亚洲久久久国产精品| 国产成人免费观看mmmm| 超碰97精品在线观看| 亚洲欧美精品自产自拍| 日韩一卡2卡3卡4卡2021年| 在现免费观看毛片| 看免费成人av毛片| av天堂久久9| 日本91视频免费播放| 国产精品偷伦视频观看了| 中国国产av一级| 国产亚洲av片在线观看秒播厂| 在线观看免费日韩欧美大片| 我要看黄色一级片免费的| 99re6热这里在线精品视频| 极品人妻少妇av视频| 一级黄片播放器| 久久精品人人爽人人爽视色| 久久久久久久久免费视频了| 一区福利在线观看| 成人国语在线视频| 乱人伦中国视频| 只有这里有精品99| 97人妻天天添夜夜摸| 国产成人精品在线电影| 天堂8中文在线网| 精品人妻一区二区三区麻豆| 久9热在线精品视频| 亚洲午夜精品一区,二区,三区| 黄色视频不卡| 亚洲精品国产色婷婷电影| 国产日韩欧美亚洲二区| 69精品国产乱码久久久| 99热全是精品| 人成视频在线观看免费观看| 中文字幕精品免费在线观看视频| 午夜福利在线免费观看网站| av在线播放精品| 80岁老熟妇乱子伦牲交| 成人国语在线视频| 老司机在亚洲福利影院| 免费观看av网站的网址| 中文字幕最新亚洲高清| 老汉色∧v一级毛片| 成人亚洲精品一区在线观看| 看免费成人av毛片| 韩国高清视频一区二区三区| 男女午夜视频在线观看| 久久久久久久精品精品| 在线亚洲精品国产二区图片欧美| 亚洲人成网站在线观看播放| 亚洲av电影在线进入| 国产成人一区二区三区免费视频网站 | 自拍欧美九色日韩亚洲蝌蚪91| 欧美人与性动交α欧美软件| 欧美日韩亚洲高清精品| 亚洲国产欧美网| 涩涩av久久男人的天堂| 欧美日韩综合久久久久久| 免费av中文字幕在线| 又紧又爽又黄一区二区| 国产日韩欧美在线精品| a级毛片在线看网站| 啦啦啦在线观看免费高清www| 热re99久久国产66热| 国产日韩欧美亚洲二区| 天天添夜夜摸| 精品久久久久久电影网| 亚洲自偷自拍图片 自拍| 丝袜喷水一区| 少妇裸体淫交视频免费看高清 | 亚洲欧美一区二区三区久久| 一个人免费看片子| 一区二区av电影网| 婷婷成人精品国产| 国产免费现黄频在线看| 99国产精品一区二区蜜桃av | 国产主播在线观看一区二区 | 中文字幕高清在线视频| tube8黄色片| 99精国产麻豆久久婷婷| √禁漫天堂资源中文www| 国产精品国产三级国产专区5o| 在线观看免费视频网站a站| 中文欧美无线码| 欧美日韩一级在线毛片| 国产精品一区二区精品视频观看| 午夜免费男女啪啪视频观看| 99精品久久久久人妻精品| 丝袜美足系列| av视频免费观看在线观看| 免费观看a级毛片全部| 不卡av一区二区三区| 国产又爽黄色视频| 亚洲五月婷婷丁香| 99国产精品99久久久久| 丁香六月天网| 青草久久国产| 欧美 亚洲 国产 日韩一| 五月开心婷婷网| 两个人看的免费小视频| avwww免费| 如日韩欧美国产精品一区二区三区| 老司机影院成人| av有码第一页| 伦理电影免费视频| 久久狼人影院| 久久这里只有精品19| 侵犯人妻中文字幕一二三四区| 国产有黄有色有爽视频| 黑人欧美特级aaaaaa片| 美女午夜性视频免费| 日韩大片免费观看网站| 天天躁夜夜躁狠狠躁躁| 1024视频免费在线观看| 女人爽到高潮嗷嗷叫在线视频| a级片在线免费高清观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久人人97超碰香蕉20202| 丝袜人妻中文字幕| 欧美日韩av久久| 99国产综合亚洲精品| 一级毛片女人18水好多 | 久9热在线精品视频| 久久国产精品影院| www.av在线官网国产| 国产男女内射视频| 亚洲精品国产av蜜桃| 纵有疾风起免费观看全集完整版| 男女国产视频网站| 国产欧美日韩一区二区三区在线| 免费女性裸体啪啪无遮挡网站| 欧美国产精品一级二级三级| 国产亚洲一区二区精品| av有码第一页| 我的亚洲天堂| 无遮挡黄片免费观看| 欧美人与性动交α欧美精品济南到| 少妇 在线观看| 丝袜美腿诱惑在线| 高清不卡的av网站| 一区二区日韩欧美中文字幕| 各种免费的搞黄视频| 久久人人97超碰香蕉20202| 国产一区二区三区av在线| 99国产精品一区二区三区| 国产精品一区二区在线观看99| 久久综合国产亚洲精品| 在线观看免费视频网站a站| 菩萨蛮人人尽说江南好唐韦庄| 在线观看www视频免费| 免费女性裸体啪啪无遮挡网站| 免费在线观看日本一区| 欧美97在线视频| 少妇人妻 视频| 亚洲av综合色区一区| 一级,二级,三级黄色视频| 一本综合久久免费| 美女高潮到喷水免费观看| 亚洲av电影在线观看一区二区三区| 国产日韩欧美亚洲二区| 别揉我奶头~嗯~啊~动态视频 | 国产精品久久久久久精品古装| 亚洲国产av影院在线观看| 青草久久国产| 国产国语露脸激情在线看| 久久ye,这里只有精品| 久久精品久久久久久久性| 国产成人精品在线电影| 亚洲国产中文字幕在线视频| 亚洲精品国产av蜜桃| 久久性视频一级片| 97精品久久久久久久久久精品| 亚洲一码二码三码区别大吗| 人妻人人澡人人爽人人| www.av在线官网国产| 91老司机精品| 午夜福利乱码中文字幕| 一本色道久久久久久精品综合| 久久久久久人人人人人| 天天躁夜夜躁狠狠久久av| 免费日韩欧美在线观看| av有码第一页| 欧美日韩视频高清一区二区三区二| 国产淫语在线视频| 午夜激情久久久久久久| 97在线人人人人妻| 视频区图区小说| 建设人人有责人人尽责人人享有的| 国产精品二区激情视频| 少妇的丰满在线观看| 亚洲国产成人一精品久久久| 免费看十八禁软件| 久久综合国产亚洲精品| 国产精品久久久av美女十八| 汤姆久久久久久久影院中文字幕| 亚洲欧美清纯卡通| 久久精品人人爽人人爽视色| 我的亚洲天堂| 免费在线观看完整版高清| 精品欧美一区二区三区在线| 国产爽快片一区二区三区| 久久免费观看电影| 亚洲欧美一区二区三区久久| 久久亚洲国产成人精品v| 可以免费在线观看a视频的电影网站| 日韩制服丝袜自拍偷拍| 91国产中文字幕| 丝袜脚勾引网站| 男女之事视频高清在线观看 | 美国免费a级毛片| 好男人电影高清在线观看| 成年人黄色毛片网站| 一级黄片播放器| 亚洲av欧美aⅴ国产| 日日夜夜操网爽| 国产欧美日韩精品亚洲av| 搡老乐熟女国产| 麻豆乱淫一区二区| 最近最新中文字幕大全免费视频 | 久久国产精品大桥未久av| 午夜福利在线免费观看网站| 亚洲成国产人片在线观看| 飞空精品影院首页| videosex国产| 高潮久久久久久久久久久不卡| 赤兔流量卡办理| 美女扒开内裤让男人捅视频| 国产亚洲av高清不卡| 亚洲激情五月婷婷啪啪| 色婷婷久久久亚洲欧美| 黄片小视频在线播放| 国产激情久久老熟女| 免费在线观看影片大全网站 | 国产黄频视频在线观看| 狂野欧美激情性bbbbbb| 日韩熟女老妇一区二区性免费视频| 欧美日韩国产mv在线观看视频| 亚洲国产日韩一区二区| 不卡av一区二区三区| 国产成人影院久久av| 亚洲精品第二区| 99精品久久久久人妻精品| 国产成人精品久久二区二区免费| 国产一区二区三区综合在线观看| 亚洲自偷自拍图片 自拍| 日韩 欧美 亚洲 中文字幕| 亚洲情色 制服丝袜| 老司机影院毛片| 高清不卡的av网站| 99国产精品一区二区蜜桃av | 国产成人精品久久久久久| 两人在一起打扑克的视频| 久久久久网色| 亚洲中文av在线| 成年女人毛片免费观看观看9 | 亚洲,欧美,日韩| 亚洲欧美成人综合另类久久久| 少妇粗大呻吟视频| 天堂中文最新版在线下载| 亚洲午夜精品一区,二区,三区| 国产免费福利视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 欧美变态另类bdsm刘玥| bbb黄色大片| 中文字幕制服av| 天天操日日干夜夜撸| 美女国产高潮福利片在线看| 亚洲人成网站在线观看播放| www.精华液| 丰满饥渴人妻一区二区三| 一级片免费观看大全| 少妇猛男粗大的猛烈进出视频| 久久久国产精品麻豆| 在线观看免费高清a一片| 亚洲第一av免费看| 婷婷丁香在线五月| 国产亚洲精品第一综合不卡| 一本色道久久久久久精品综合| 一边摸一边抽搐一进一出视频| 肉色欧美久久久久久久蜜桃| 伊人亚洲综合成人网| 日韩中文字幕视频在线看片| 国产黄色视频一区二区在线观看| 精品一区二区三卡| 汤姆久久久久久久影院中文字幕| 亚洲欧美一区二区三区国产| 看免费成人av毛片| 久久 成人 亚洲| 久久久精品免费免费高清| 夫妻午夜视频| 两个人看的免费小视频| 精品一品国产午夜福利视频| 国产亚洲欧美精品永久| 国产精品一二三区在线看| 国产精品偷伦视频观看了| 午夜福利视频精品| 成年动漫av网址| 青草久久国产| 女人被躁到高潮嗷嗷叫费观| 如日韩欧美国产精品一区二区三区| 亚洲情色 制服丝袜| 女警被强在线播放| 97人妻天天添夜夜摸| 国产亚洲欧美在线一区二区| 一级毛片黄色毛片免费观看视频| 波野结衣二区三区在线| 成人国产av品久久久| 久久精品成人免费网站| 国产成人精品久久二区二区免费| 国产免费现黄频在线看| 国精品久久久久久国模美| 母亲3免费完整高清在线观看| 美女视频免费永久观看网站| 精品国产一区二区久久| av视频免费观看在线观看| 王馨瑶露胸无遮挡在线观看| 女人久久www免费人成看片| 男人爽女人下面视频在线观看| 99热国产这里只有精品6| 日本五十路高清| 日本色播在线视频| 天天影视国产精品| 国产高清videossex| 国产有黄有色有爽视频| 亚洲av综合色区一区| 婷婷色综合大香蕉| 两个人免费观看高清视频| 丝袜美足系列| 九草在线视频观看| 欧美 亚洲 国产 日韩一| 午夜av观看不卡| 国产精品香港三级国产av潘金莲 | 免费日韩欧美在线观看| 操出白浆在线播放| 国产欧美日韩一区二区三 | 18在线观看网站| 国产精品国产三级国产专区5o| 狂野欧美激情性xxxx| 亚洲精品美女久久久久99蜜臀 | 亚洲男人天堂网一区| 老司机亚洲免费影院| 18禁黄网站禁片午夜丰满| 国产av一区二区精品久久| 色婷婷av一区二区三区视频| 男人添女人高潮全过程视频| 国产成人精品久久二区二区免费| 一区二区三区精品91| 亚洲国产毛片av蜜桃av| 中国美女看黄片| 日本五十路高清| 日本av手机在线免费观看| 纵有疾风起免费观看全集完整版| 午夜精品国产一区二区电影| 精品欧美一区二区三区在线| 午夜福利,免费看| av国产久精品久网站免费入址| 国产精品麻豆人妻色哟哟久久| 亚洲欧洲精品一区二区精品久久久| 免费女性裸体啪啪无遮挡网站| 狠狠精品人妻久久久久久综合| av天堂久久9| 久久人妻熟女aⅴ| 亚洲精品第二区| 久久久精品国产亚洲av高清涩受| 国产主播在线观看一区二区 | 99久久人妻综合| 国产一区二区 视频在线| 交换朋友夫妻互换小说| 日韩中文字幕欧美一区二区 | 亚洲av在线观看美女高潮|