賈守偉， 史 潔， 高會(huì)旺

(中國(guó)海洋大學(xué)海洋環(huán)境與生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室，山東 青島 266100)

?

膠州灣營(yíng)養(yǎng)鹽限制作用空間差異和長(zhǎng)期變化的數(shù)值研究*

賈守偉， 史 潔**， 高會(huì)旺

(中國(guó)海洋大學(xué)海洋環(huán)境與生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室，山東 青島 266100)

以膠州灣為研究區(qū)域，建立箱式生態(tài)動(dòng)力學(xué)模型，模擬得到膠州灣浮游植物生長(zhǎng)受營(yíng)養(yǎng)鹽限制情況的季節(jié)變化；通過(guò)數(shù)值實(shí)驗(yàn)分析了膠州灣營(yíng)養(yǎng)鹽限制的區(qū)域差異，并進(jìn)一步比較了1960、1980年代及本世紀(jì)初營(yíng)養(yǎng)鹽限制情況的區(qū)別。結(jié)果表明，膠州灣內(nèi)的營(yíng)養(yǎng)鹽限制作用有顯著的季節(jié)變化。春季由于浮游植物旺發(fā)消耗大量營(yíng)養(yǎng)鹽，使得灣內(nèi)營(yíng)養(yǎng)鹽限制春末最強(qiáng)；夏季是河流豐水季，營(yíng)養(yǎng)鹽補(bǔ)充充足，灣內(nèi)營(yíng)養(yǎng)鹽限制最弱。從不同區(qū)域來(lái)看，養(yǎng)殖區(qū)全年表現(xiàn)為硅限制；河口區(qū)和灣口區(qū)夏季為磷限制，其他季節(jié)為硅限制。營(yíng)養(yǎng)鹽限制作用灣口區(qū)最強(qiáng)，養(yǎng)殖區(qū)次之，河口區(qū)最弱。受此影響，浮游植物生物量河口區(qū)最高，灣口區(qū)最低。從長(zhǎng)期變化來(lái)看，1960年代，膠州灣內(nèi)的營(yíng)養(yǎng)鹽限制作用以氮限制為主，1980年代呈硅和磷交替限制，逐漸演變到本世紀(jì)初的以硅限制為主，而且由于營(yíng)養(yǎng)鹽限制作用的逐漸增強(qiáng)，浮游植物生物量也呈現(xiàn)不斷下降的趨勢(shì)。

膠州灣； 營(yíng)養(yǎng)鹽限制； 區(qū)域差異； 長(zhǎng)期變化

海洋浮游植物是海洋生態(tài)系統(tǒng)中的初級(jí)生產(chǎn)者和食物鏈的基礎(chǔ)，在海洋生態(tài)系統(tǒng)的物質(zhì)循環(huán)和能量轉(zhuǎn)化過(guò)程中起著重要作用，其動(dòng)態(tài)變化直接影響到海洋生態(tài)系統(tǒng)的功能與結(jié)構(gòu)。浮游植物的生長(zhǎng)受到光照、溫度和營(yíng)養(yǎng)鹽等多種因素的影響。其中，營(yíng)養(yǎng)鹽水平受到人為活動(dòng)的影響，關(guān)于營(yíng)養(yǎng)鹽對(duì)浮游植物生長(zhǎng)控制作用的研究有助于污染物減排措施的制定和富營(yíng)養(yǎng)化的治理，因此得到了國(guó)內(nèi)外學(xué)者的廣泛關(guān)注[1]。若環(huán)境中某種營(yíng)養(yǎng)鹽濃度低于浮游植物生長(zhǎng)所需的最低閾值，則會(huì)產(chǎn)生該種營(yíng)養(yǎng)鹽的限制作用，營(yíng)養(yǎng)鹽限制直接影響海洋的初級(jí)生產(chǎn)力和浮游植物群落結(jié)構(gòu)[2-3]，從長(zhǎng)遠(yuǎn)來(lái)看還可能導(dǎo)致食物鏈的改變[3]。近岸海域作為從陸地到海洋的過(guò)渡區(qū)域，受到人類(lèi)活動(dòng)的影響最為顯著，而且其生產(chǎn)力水平也往往較高。隨著氣候變化和人類(lèi)活動(dòng)的影響，在近岸海域，尤其是半封閉海灣，營(yíng)養(yǎng)鹽各來(lái)源的量和比例往往具有差異，因此會(huì)導(dǎo)致?tīng)I(yíng)養(yǎng)鹽限制作用產(chǎn)生顯著的時(shí)空變化，并對(duì)浮游植物的生長(zhǎng)造成影響[4-7]。

膠州灣是受人類(lèi)活動(dòng)影響顯著的半封閉海灣，外與黃海相通，水體平均存留時(shí)間為52d[8]，其營(yíng)養(yǎng)鹽結(jié)構(gòu)在北方海灣中具有較強(qiáng)的代表性[9]。膠州灣有十余條入海河流，其中位于西北部的大沽河占入?？偭髁康?0%[10]，是膠州灣重要的營(yíng)養(yǎng)鹽來(lái)源，其帶來(lái)的營(yíng)養(yǎng)鹽直接影響河口區(qū)域。膠州灣內(nèi)水產(chǎn)養(yǎng)殖業(yè)發(fā)達(dá)，養(yǎng)殖區(qū)多集中在東北部，是典型的對(duì)蝦和蛤蜊養(yǎng)殖區(qū)，餌料的過(guò)量投放和養(yǎng)殖生物的排泄為該區(qū)域輸送大量的氮、磷營(yíng)養(yǎng)鹽，使該區(qū)域內(nèi)這兩種營(yíng)養(yǎng)鹽濃度遠(yuǎn)遠(yuǎn)高于全灣平均水平[11]。膠州灣內(nèi)的營(yíng)養(yǎng)鹽濃度常年高于外海，通過(guò)灣口處與外海的交換，大量營(yíng)養(yǎng)鹽被輸送到黃海，因此灣口區(qū)域的營(yíng)養(yǎng)鹽結(jié)構(gòu)也有別于其他區(qū)域[12]。

現(xiàn)有研究表明，在膠州灣內(nèi)，相較于光照和水溫，營(yíng)養(yǎng)鹽對(duì)浮游植物生長(zhǎng)的限制作用更為重要[9]。基于觀測(cè)資料，已有研究闡明了膠州灣內(nèi)營(yíng)養(yǎng)鹽濃度及浮游植物生物量的季節(jié)、長(zhǎng)期變化和空間分布特征[12,14-18]。王勇等[19]在膠州灣通過(guò)營(yíng)養(yǎng)鹽添加實(shí)驗(yàn)發(fā)現(xiàn)，網(wǎng)采浮游植物對(duì)鐵或氮、磷的添加有比較明顯的響應(yīng)。姚云[20]通過(guò)分析氮磷比探討了膠州灣內(nèi)限制性營(yíng)養(yǎng)鹽的種類(lèi)及季節(jié)變化，認(rèn)為相較于氮和磷而言，帶有硅的評(píng)價(jià)結(jié)果能夠更好地反映膠州灣的富營(yíng)養(yǎng)化狀況。觀測(cè)和實(shí)驗(yàn)?zāi)軌蚍从衬z州灣營(yíng)養(yǎng)鹽結(jié)構(gòu)的現(xiàn)狀，但卻存在著時(shí)間不連續(xù)性，并難以區(qū)分各個(gè)營(yíng)養(yǎng)鹽來(lái)源對(duì)膠州灣營(yíng)養(yǎng)鹽結(jié)構(gòu)的影響。海洋生態(tài)動(dòng)力學(xué)模型是定量研究生態(tài)系統(tǒng)中各過(guò)程相互作用的有效方法，不僅可以反演現(xiàn)場(chǎng)觀測(cè)結(jié)果，還可以進(jìn)行機(jī)制分析，模擬和預(yù)測(cè)生態(tài)系統(tǒng)的變化?，F(xiàn)有的膠州灣生態(tài)系統(tǒng)動(dòng)力學(xué)模型研究能夠模擬膠州灣營(yíng)養(yǎng)鹽濃度和浮游植物生物量的季節(jié)變化，計(jì)算得到膠州灣內(nèi)的營(yíng)養(yǎng)鹽收支狀況，但是沒(méi)有分析不同營(yíng)養(yǎng)鹽來(lái)源對(duì)膠州灣營(yíng)養(yǎng)鹽限制情況的影響，缺少源的長(zhǎng)期變化導(dǎo)致膠州灣生態(tài)系統(tǒng)演變的研究[21-24]。

本文選取膠州灣為研究海域(見(jiàn)圖1)，建立膠州灣箱式生態(tài)動(dòng)力學(xué)模型，區(qū)分河口區(qū)、養(yǎng)殖區(qū)和灣口區(qū)3個(gè)典型區(qū)域，模擬并探討膠州灣內(nèi)營(yíng)養(yǎng)鹽限制作用的區(qū)域差異，分析膠州灣不同營(yíng)養(yǎng)鹽來(lái)源對(duì)灣內(nèi)營(yíng)養(yǎng)鹽限制作用以及浮游植物生物量的影響，模擬并分析了膠州灣營(yíng)養(yǎng)鹽結(jié)構(gòu)隨營(yíng)養(yǎng)鹽來(lái)源變化而產(chǎn)生的長(zhǎng)期變化趨勢(shì)。

圖1 膠州灣地理位置Fig.1 Location of Jiaozhou Bay

1 研究方法

1.1 模型介紹

本文在膠州灣建立了箱式生態(tài)系統(tǒng)動(dòng)力學(xué)模型(見(jiàn)圖2)，來(lái)模擬營(yíng)養(yǎng)鹽濃度和浮游植物生物量的季節(jié)變化。模型包括無(wú)機(jī)氮(DIN： dissolved inorganic nitrogen)、無(wú)機(jī)磷(DIP：dissolved inorganic phosphorus)、硅酸鹽(SIL：silicate)3種營(yíng)養(yǎng)鹽(N)，浮游植物(P)和碎屑(D)共5個(gè)狀態(tài)變量。浮游植物利用太陽(yáng)能進(jìn)行光合作用，攝取營(yíng)養(yǎng)鹽合成有機(jī)物，同時(shí)呼吸作用釋放出營(yíng)養(yǎng)鹽，浮游植物死亡后成為碎屑的一部分，碎屑經(jīng)過(guò)礦化作用向水體補(bǔ)充營(yíng)養(yǎng)鹽。

模型的控制方程如下：

(1)

(2)

(3)

圖2 模型概念圖Fig.2 Schematic illustration of the model

浮游植物呼吸作用分為光呼吸和基礎(chǔ)呼吸過(guò)程，并考慮溫度對(duì)呼吸作用的影響。

(4)

rT=er·(T-T0)

(5)

浮游植物的生長(zhǎng)受光照條件、溫度、營(yíng)養(yǎng)鹽的共同影響，其表達(dá)式為：

(6)

其中：g(T)為溫度限制因子；VL、VN分別為光和營(yíng)養(yǎng)鹽限制因子。

溫度對(duì)浮游植物生長(zhǎng)的限制采用Q10法則，即:

(7)

光對(duì)浮游植物的生長(zhǎng)有重要影響，其限制作用的表達(dá)函數(shù)為:

(8)

其中：Iopt為最優(yōu)光強(qiáng)；I為水體有效太陽(yáng)輻射。

本研究中采用Michaeli-Menten(M-M)關(guān)系來(lái)量化浮游植物生長(zhǎng)受到環(huán)境中營(yíng)養(yǎng)鹽限制作用的強(qiáng)弱，定義營(yíng)養(yǎng)鹽限制因子為：

(9)

其中：N為營(yíng)養(yǎng)鹽的濃度；KW為浮游植物生長(zhǎng)吸收營(yíng)養(yǎng)鹽的半飽和常數(shù)。由營(yíng)養(yǎng)鹽限制因子定義可知，在KW一定的情況下，營(yíng)養(yǎng)鹽濃度越高，VN值越接近于1，對(duì)浮游植物的限制作用越小。浮游植物生長(zhǎng)受到營(yíng)養(yǎng)鹽的限制取決于VN值最小的營(yíng)養(yǎng)鹽，即：

VN=min(VDIN，VDIP，VSIL)

(10)

模型中參數(shù)的取值(見(jiàn)表1)采用前人在膠州灣及中國(guó)近海的生態(tài)模型研究結(jié)果[24-26]。

表1 模式中主要參數(shù)說(shuō)明及取值Table 1 Values of parameters used in the model

1.2 模型設(shè)置和數(shù)據(jù)來(lái)源

營(yíng)養(yǎng)鹽的源匯項(xiàng)包括河流輸入、大氣沉降、養(yǎng)殖活動(dòng)和與黃海交換。其數(shù)值大小均采用本世紀(jì)初膠州灣的模擬和觀測(cè)結(jié)果(見(jiàn)表2)。入灣河流中大沽河對(duì)膠州灣內(nèi)的營(yíng)養(yǎng)鹽貢獻(xiàn)最大，故模型中河流對(duì)膠州灣的營(yíng)養(yǎng)鹽貢獻(xiàn)僅考慮大沽河一條徑流，大沽河的營(yíng)養(yǎng)鹽輸入量由流量乘以營(yíng)養(yǎng)鹽的濃度計(jì)算得出。大沽河流量采用浦祥[28]模擬得到的2003年大沽河日均徑流量結(jié)果(見(jiàn)圖4)。大沽河營(yíng)養(yǎng)鹽濃度為2002年干、濕季的觀測(cè)結(jié)果[29]。模型計(jì)算時(shí)將養(yǎng)殖活動(dòng)和大氣沉降一年所貢獻(xiàn)的營(yíng)養(yǎng)鹽總量平均分配到各時(shí)間步長(zhǎng)。養(yǎng)殖活動(dòng)對(duì)膠州灣每年貢獻(xiàn)無(wú)機(jī)氮54×106mol、無(wú)機(jī)磷3.8×106mol[30]。大氣沉降每年為膠州灣貢獻(xiàn)無(wú)機(jī)氮56×106mol，無(wú)機(jī)磷0.42×106mol，硅酸鹽0.93×106mol[31]?？紤]膠州灣與黃海之間的營(yíng)養(yǎng)鹽交換，將黃海水中各種營(yíng)養(yǎng)鹽的濃度視為常數(shù)，其營(yíng)養(yǎng)鹽通量取決于黃海與膠州灣內(nèi)營(yíng)養(yǎng)鹽濃度的差值和海水交換率[8, 27]。黃海的營(yíng)養(yǎng)鹽濃度參考已發(fā)表的觀測(cè)資料[12,32-33]。

模型中的物理強(qiáng)迫包括海水溫度和海面光強(qiáng)，二者均用余弦函數(shù)來(lái)擬合，模型采用以下公式進(jìn)行計(jì)算[23]：

Temperature=13.5-11.0×(2π×(t-40)/365)

(11)

Rad=240-140×cos(2π×(t+10)/365)

(12)

其中：Temperature為水溫，單位為℃；Rad為海表面的光照強(qiáng)度，單位為W/m2；t的單位為d，t=1代表1月的第1天。

模型采用冬季觀測(cè)的平均濃度作為初始場(chǎng)，無(wú)機(jī)氮、無(wú)機(jī)磷和硅酸鹽的濃度分別設(shè)為10、0.5、3.5μmol/L，葉綠素濃度為0.8mg/m3[14-15]。開(kāi)始時(shí)間為1月1日，模型時(shí)間步長(zhǎng)為1h，模型在以上強(qiáng)迫下連續(xù)運(yùn)行3年，用第3年的結(jié)果進(jìn)行分析。

圖3 本世紀(jì)初大沽河徑流量季節(jié)變化Fig. 3 Seasonal variations of Dagu River discharge in the beginning of this century

1960年代1960s1980年代1980s21世紀(jì)初2000s單位Unit來(lái)源Reference河流量①7.55.03.0108m3/a[9]氮濃度②81260400μmol/L[28-29]磷濃度③1.55.37.0μmol/L[28-29]硅濃度④600420150μmol/L[9,28]養(yǎng)殖排放氮⑤00.6854106mol[20]養(yǎng)殖排放磷⑥00.053.8106mol[20]灣外氮濃度⑦2.256μmol/L[12,32-33]灣外磷濃度⑧0.170.300.30μmol/L[12,32-33]灣外硅濃度⑨65.54μmol/L[12,32-33]

Note:①Flow；②Nitrogen concentration；③Phosphorus concentration ；④Silicion concentration；⑤Nitrogen flux of aquaculture；⑥Phosphorus flux of aquaculture；⑦Nitrogen concentration of Yellow Sea；⑧Phosphorus concentration of Yellow Sea；⑨Silicion concentration of YS

2 結(jié)果

2.1 敏感性分析

在一個(gè)模式中，參數(shù)的變化對(duì)模擬結(jié)果的影響程度反映了模式結(jié)果對(duì)該參數(shù)的敏感程度。由于生物過(guò)程本身的非線性特征，生態(tài)變量的計(jì)算值會(huì)對(duì)模型中一些參數(shù)的變化非常敏感。因此，參數(shù)的敏感性分析是生態(tài)模型研究中必不可少的步驟。通過(guò)敏感性實(shí)驗(yàn)，找出生態(tài)系統(tǒng)變化敏感參數(shù)，為了比較不同參數(shù)的敏感程度，定義敏感度為：

(12)

其中：ΔF/F為生態(tài)模型中浮游植物生物量的變化率；Δα/α為模型中與生態(tài)變量F有關(guān)的某個(gè)生物參數(shù)的變化率。表1提供了10個(gè)參數(shù)的敏感性數(shù)值，浮游植物生物量對(duì)最大生長(zhǎng)率、浮游植物基礎(chǔ)呼吸率、最優(yōu)光強(qiáng)等參數(shù)較敏感，而對(duì)礦化率、氮營(yíng)養(yǎng)鹽半飽和常數(shù)等參數(shù)相對(duì)不敏感。各參數(shù)分別變化或減小相同比例，其敏感度并不相同，表明各參數(shù)的變化是非線性變化，如浮游植物死亡率，增大50%其敏感度為0.65，減小50%其敏感度為1.91。75%的參數(shù)敏感度小于0.5，因此本文建立的模型是穩(wěn)定的。

2.2 本世紀(jì)初營(yíng)養(yǎng)鹽濃度和浮游植物生物量的季節(jié)變化

應(yīng)用所建立的膠州灣箱式生態(tài)動(dòng)力學(xué)模型，模擬本世紀(jì)初的膠州灣營(yíng)養(yǎng)鹽濃度和浮游植物生物量的季節(jié)變化規(guī)律(見(jiàn)圖4)。

膠州灣內(nèi)3種營(yíng)養(yǎng)鹽的濃度均為夏季最高，春季最低。春季，由于浮游植物生長(zhǎng)攝取大量營(yíng)養(yǎng)鹽，導(dǎo)致膠州灣內(nèi)營(yíng)養(yǎng)鹽濃度一直維持在較低水平。夏季，大沽河進(jìn)入汛期，河流輸入帶來(lái)的營(yíng)養(yǎng)鹽顯著增加，使灣內(nèi)各營(yíng)養(yǎng)鹽濃度迅速升高，到8月中旬達(dá)到峰值，無(wú)機(jī)氮、無(wú)機(jī)磷和硅酸鹽的濃度峰值分別為41.6、1.0和10.3 μmol/L，出現(xiàn)時(shí)間均在第220天左右。豐水期之后，浮游植物生長(zhǎng)消耗大量的營(yíng)養(yǎng)鹽，使其濃度不斷下降。膠州灣內(nèi)浮游植物生物量的季節(jié)變化特征明顯，模擬結(jié)果顯示呈單峰狀，葉綠素年平均濃度為1.4mg/m3。冬季，太陽(yáng)輻射較弱，水溫較低，不利于浮游植物的生長(zhǎng)，葉綠素濃度比較低。夏季，營(yíng)養(yǎng)鹽充足，溫度和光照條件適宜，浮游植物大量繁殖，葉綠素濃度達(dá)到峰值4.1mg/m3，稍微滯后于營(yíng)養(yǎng)鹽濃度峰值。此后，由于營(yíng)養(yǎng)鹽大量消耗及太陽(yáng)輻射的限制作用不斷增強(qiáng)，抑制了浮游植物的生長(zhǎng)，葉綠素濃度降低，冬季達(dá)到最低。模擬得到的本世紀(jì)初膠州灣營(yíng)養(yǎng)鹽濃度和浮游植物生物量的季節(jié)變化規(guī)律與觀測(cè)結(jié)果較為吻合?，F(xiàn)場(chǎng)觀測(cè)發(fā)現(xiàn)某些年份浮游植物生物量在2月也存在峰值，有學(xué)者[14, 34]推測(cè)，這是由于某些冷水種藻類(lèi)冬季大量生長(zhǎng)造成的，但并未證實(shí)。在本模型中溫度對(duì)浮游植物生長(zhǎng)的限制采用Q10法則，因此難以體現(xiàn)冷水種浮游植物在冬季的旺發(fā)。

圖4 營(yíng)養(yǎng)鹽濃度和浮游植物生物量的季節(jié)變化

2.3 本世紀(jì)初營(yíng)養(yǎng)鹽限制的季節(jié)變化

3種營(yíng)養(yǎng)鹽對(duì)浮游植物的限制作用存在強(qiáng)弱差異，因此探討營(yíng)養(yǎng)鹽限制作用的季節(jié)變化對(duì)于理解營(yíng)養(yǎng)鹽結(jié)構(gòu)對(duì)浮游植物生長(zhǎng)的影響有重要作用。通過(guò)比較模擬結(jié)果中3種營(yíng)養(yǎng)鹽限制因子的數(shù)值大小，來(lái)確定限制浮游植物生長(zhǎng)的營(yíng)養(yǎng)鹽種類(lèi)，即營(yíng)養(yǎng)鹽限制因子數(shù)值最小的被認(rèn)為是限制浮游植物生長(zhǎng)的營(yíng)養(yǎng)鹽。

比較3種營(yíng)養(yǎng)鹽限制因子可以看出(見(jiàn)圖5)，3種營(yíng)養(yǎng)鹽限制因子的數(shù)值均表現(xiàn)為春季最小，夏季最大，這說(shuō)明膠州灣內(nèi)的營(yíng)養(yǎng)鹽限制作用夏季最弱，秋季和冬季次之，春季限制作用最強(qiáng)。夏季，大沽河徑流量達(dá)到全年最大，帶來(lái)充足的營(yíng)養(yǎng)鹽，削弱了營(yíng)養(yǎng)鹽限制作用。秋季，隨著徑流輸入的減少，加上浮游植物的生長(zhǎng)不斷消耗營(yíng)養(yǎng)鹽，營(yíng)養(yǎng)鹽限制作用逐漸增強(qiáng)。經(jīng)過(guò)冬、春季節(jié)的消耗，營(yíng)養(yǎng)鹽限制作用在春末夏初(180 d左右)達(dá)到最強(qiáng)。從數(shù)值上看，氮和磷限制因子的數(shù)值全年均在0.8以上，而硅限制因子常年小于氮、磷，僅在夏季營(yíng)養(yǎng)鹽豐富時(shí)能夠達(dá)到0.75左右，由于春季浮游植物旺發(fā)的大量吸收，5～6月份期間硅限制因子達(dá)到最小值，約為0.5，即浮游植物在硅酸鹽的限制作用下，其生長(zhǎng)速率只能達(dá)到最大生長(zhǎng)率的一半。由此可見(jiàn)，硅是限制膠州灣浮游植物生長(zhǎng)的首要營(yíng)養(yǎng)鹽，而氮和磷相對(duì)充足，不限制浮游植物的生長(zhǎng)。

膠州灣是位于黃海邊緣的一個(gè)海灣，但與黃海中部海域相比，其營(yíng)養(yǎng)鹽限制作用卻有其特性。黃海中部海域，營(yíng)養(yǎng)鹽限制作用夏季最強(qiáng)，冬季最弱，且主要受磷的潛在限制。王保棟等[36]觀測(cè)發(fā)現(xiàn)，黃海大部分海域硅氮比值大于1，可以認(rèn)為硅酸鹽不會(huì)限制浮游植物的生長(zhǎng)。Zhang[37]研究發(fā)現(xiàn)，黃河、灤河、鴨綠江、大遼河等幾條黃海主要入海河流的氮磷比值在20～5000之間，遠(yuǎn)遠(yuǎn)高于Redfield系數(shù)。因此，黃海海域水體中無(wú)機(jī)氮相對(duì)過(guò)剩，而無(wú)機(jī)磷相對(duì)不足，黃海中部海域浮游植物的生長(zhǎng)明顯受磷限制。從季節(jié)上看，黃海中部水域中營(yíng)養(yǎng)鹽限制最強(qiáng)和最弱分別出現(xiàn)在夏季和冬季，而不像膠州灣分別出現(xiàn)在春季和夏季。黃海春季出現(xiàn)較強(qiáng)的浮游植物水華現(xiàn)象，大量消耗營(yíng)養(yǎng)鹽。夏季，隨著海表面溫度升高，形成較強(qiáng)的躍層，較弱的垂直混合阻礙了底部高營(yíng)養(yǎng)鹽的水向上補(bǔ)充到真光層，進(jìn)而使得此時(shí)營(yíng)養(yǎng)鹽限制達(dá)到全年最強(qiáng)。秋季水華過(guò)后，大風(fēng)混合增強(qiáng)使得上層得到營(yíng)養(yǎng)鹽的補(bǔ)充，并且溫度降低和太陽(yáng)輻射減弱，浮游植物生長(zhǎng)速率減慢。到了冬季，低浮游植物生物量和高營(yíng)養(yǎng)鹽濃度，使得此時(shí)營(yíng)養(yǎng)鹽限制為全年最弱。黃海中部水域營(yíng)養(yǎng)鹽限制的季節(jié)變化是受到浮游植物生長(zhǎng)和躍層生消的共同作用，而膠州灣水深較淺，跨躍層補(bǔ)充的季節(jié)變化并不顯著，因此營(yíng)養(yǎng)鹽限制的季節(jié)變化主要取決于浮游植物的生長(zhǎng)和河流的豐水、枯水期變化?？梢?jiàn)，不同的營(yíng)養(yǎng)鹽補(bǔ)充機(jī)制，會(huì)導(dǎo)致不同的海域，即使是相互連通的海域，其營(yíng)養(yǎng)鹽限制情況也會(huì)呈現(xiàn)不同的季節(jié)變化規(guī)律。

圖5 營(yíng)養(yǎng)鹽限制因子的季節(jié)變化Fig.5 Seasonal variations of nutrient limiting factors

3 討論

3.1 營(yíng)養(yǎng)鹽限制的區(qū)域差異

膠州灣內(nèi)不同區(qū)域營(yíng)養(yǎng)鹽的來(lái)源不同，并且各個(gè)營(yíng)養(yǎng)鹽來(lái)源的量和比例也有較大區(qū)別，因此營(yíng)養(yǎng)鹽濃度及結(jié)構(gòu)都存在著空間不均勻性，浮游植物生長(zhǎng)受到的營(yíng)養(yǎng)鹽限制作用也會(huì)存在較大的區(qū)域差異。

表3 不同區(qū)域的營(yíng)養(yǎng)鹽來(lái)源Table 3 Nutrient sources in different areas

Note:①Aquaculture；②Atmospheric deposition；③River input；④Exchange with YS；⑤Exchange among the different areas；⑥River mouth area；⑦Aquiculture area；⑧Bay mouth area

分別將河口區(qū)、養(yǎng)殖區(qū)和灣口區(qū)作為研究區(qū)域(見(jiàn)圖1)，考慮不同區(qū)域營(yíng)養(yǎng)鹽的來(lái)源差異(見(jiàn)表3)，將生態(tài)模型應(yīng)用于3個(gè)區(qū)域，其中灣內(nèi)不同區(qū)域間營(yíng)養(yǎng)鹽的交換量取決于該區(qū)域與全灣水平的營(yíng)養(yǎng)鹽濃度差，模擬得到各區(qū)域營(yíng)養(yǎng)鹽限制因子的季節(jié)變化(見(jiàn)圖6)。氮限制因子的數(shù)值在不同區(qū)域均大于0.8，始終是3種營(yíng)養(yǎng)鹽限制因子中最大的，說(shuō)明3個(gè)區(qū)域都沒(méi)有表現(xiàn)出氮限制，膠州灣內(nèi)無(wú)機(jī)氮的現(xiàn)存量相對(duì)于浮游植物的生長(zhǎng)需要來(lái)說(shuō)過(guò)剩。河口區(qū)硅和磷限制因子的數(shù)值區(qū)間均為0.6～0.8，在3～7月份，硅限制因子小于磷限制因子，而在8～12月份，磷限制因子小于硅限制因子。在河口區(qū)由硅和磷交替限制浮游植物的生長(zhǎng)。冬、春季，大沽河處于枯水期，硅酸鹽的來(lái)源不足，浮游植物的生長(zhǎng)主要受到硅的限制；夏季，徑流急劇增加帶來(lái)豐富的營(yíng)養(yǎng)鹽，而且硅酸鹽的濃度在汛期為枯水期的12.3倍，無(wú)機(jī)磷僅為枯水期的7倍[29]，大量的硅酸鹽被輸送到膠州灣，而無(wú)機(jī)磷相對(duì)不足，使得營(yíng)養(yǎng)鹽限制由硅限制轉(zhuǎn)變?yōu)榱紫拗啤Ｔ陴B(yǎng)殖區(qū)，硅限制因子的數(shù)值全年均在0.7以下，遠(yuǎn)小于另外兩種營(yíng)養(yǎng)鹽，說(shuō)明該區(qū)域呈顯著的硅限制。養(yǎng)殖活動(dòng)為該區(qū)域貢獻(xiàn)大量氮和磷，而硅酸鹽無(wú)充足來(lái)源，此外，養(yǎng)殖區(qū)位于膠州灣的東北角，潮流作用較弱，與其他區(qū)域的水交換相對(duì)緩慢，加劇了無(wú)機(jī)氮、無(wú)機(jī)磷的富集，阻礙了其他區(qū)域向此處硅酸鹽的補(bǔ)充。從全年來(lái)看，硅對(duì)浮游植物的限制作用春末尤為突出，浮游植物的生長(zhǎng)使硅酸鹽幾乎完全被消耗。灣口區(qū)，由于與外海水交換順暢，大量營(yíng)養(yǎng)鹽被輸運(yùn)到外海，此處營(yíng)養(yǎng)鹽濃度低于全灣平均水平，無(wú)機(jī)磷和硅酸鹽交替成為限制因子，春季為硅酸鹽限制，到了秋季轉(zhuǎn)為無(wú)機(jī)磷限制。取3種營(yíng)養(yǎng)鹽限制因子中的最小值(見(jiàn)公式10)，繪制不同區(qū)域該值的全年變化曲線，可以發(fā)現(xiàn)灣口區(qū)的數(shù)值最小，全年均在0.6以下，而河口區(qū)和養(yǎng)殖區(qū)多大于0.6，其中以河口區(qū)的數(shù)值最大，因此膠州灣內(nèi)河口區(qū)的營(yíng)養(yǎng)鹽限制作用最弱，養(yǎng)殖區(qū)次之，灣口區(qū)限制作用最強(qiáng)。

通過(guò)模型結(jié)果可以看出(見(jiàn)圖7)，各區(qū)域的浮游植物生物量季節(jié)變化均呈單峰狀，峰值出現(xiàn)在夏季，但生物量存在明顯差異。營(yíng)養(yǎng)鹽限制作用較弱的河口區(qū)，浮游植物生長(zhǎng)旺盛，浮游植物生物量最大，葉綠素年平均濃度為1.7mg/m3，高出全灣平均水平20%。反之，營(yíng)養(yǎng)鹽限制作用較強(qiáng)的灣口區(qū)，其浮游植物生物量則大大低于平均水平，年平均濃度僅為1.1 mg/m3，比全灣平均值低21%。養(yǎng)殖區(qū)浮游植物生物量介于其他兩區(qū)域之間。吳玉霖等[34]觀測(cè)也發(fā)現(xiàn)，膠州灣內(nèi)葉綠素高值區(qū)多集中在西北部水域，而該水域也是大沽河影響較大的區(qū)域?？梢?jiàn)，營(yíng)養(yǎng)鹽限制作用的區(qū)域差異，最終導(dǎo)致不同區(qū)域間浮游植物生物量及初級(jí)生產(chǎn)力也有較大差別。

圖6 不同區(qū)域營(yíng)養(yǎng)鹽限制因子的季節(jié)變化

圖7 不同區(qū)域浮游植物生物量的季節(jié)變化

3.2 營(yíng)養(yǎng)鹽限制的年際差異

近40年來(lái)，膠州灣內(nèi)的浮游植物生物量不斷下降調(diào)查表明，膠州灣內(nèi)的海水溫度和鹽度在此期間沒(méi)有大的變化，而且大部分水域的水質(zhì)和沉積環(huán)境質(zhì)量較好[35]。因此，膠州灣內(nèi)浮游植物生物量的減少，可能與膠州灣內(nèi)的營(yíng)養(yǎng)鹽結(jié)構(gòu)變化密切相關(guān)。表2給出了不同時(shí)期河流輸入、養(yǎng)殖排放和外海交換對(duì)膠州灣營(yíng)養(yǎng)鹽的貢獻(xiàn)，光照和水溫仍然用余弦函數(shù)擬合，應(yīng)用所建立的箱式生態(tài)動(dòng)力學(xué)模型，分別模擬得到20世紀(jì)60年代、80年代營(yíng)養(yǎng)鹽濃度及限制因子的季節(jié)變化，結(jié)合已本世紀(jì)初的結(jié)果，比較分析40年來(lái)膠州灣營(yíng)養(yǎng)鹽限制情況的長(zhǎng)期變化特征及其與浮游植物生物量的關(guān)系。

表4 營(yíng)養(yǎng)鹽濃度長(zhǎng)期變化模擬結(jié)果與觀測(cè)結(jié)果對(duì)比[35]Table 4 Comparison between the simulation and observation in the 1960s and the 1980s

Note:①Simulation；②Observation

從長(zhǎng)期變化模擬結(jié)果來(lái)看，膠州灣營(yíng)養(yǎng)鹽濃度和限制作用變化顯著(見(jiàn)表4；圖8)。20世紀(jì)60年代，氮限制因子的數(shù)值在3種營(yíng)養(yǎng)鹽中最小，最小值約為0.7，磷限制因子略大于氮，硅的限制因子最大，其數(shù)值在夏季達(dá)到0.95，說(shuō)明該時(shí)期膠州灣內(nèi)浮游植物的生長(zhǎng)主要受到氮的限制，其次為磷，硅的限制作用最小。這是因?yàn)樵?0世紀(jì)60年代，膠州灣內(nèi)無(wú)機(jī)氮濃度相對(duì)較低，全年的氮磷比平均值僅為2.8，浮游植物的生長(zhǎng)消耗大量的營(yíng)養(yǎng)鹽，海水中的氮常常被耗盡，特別是在生長(zhǎng)高峰期，很多區(qū)域的站位檢測(cè)不到無(wú)機(jī)氮[35]。20世紀(jì)80年代，硅限制因子的數(shù)值較小(0.6～0.8)，僅在夏季磷限制因子的數(shù)值略小于硅，說(shuō)明全年內(nèi)浮游植物的生長(zhǎng)受硅限制較強(qiáng)，僅在夏季受磷限制。這是由于1980年代人類(lèi)活動(dòng)排放導(dǎo)致膠州灣內(nèi)無(wú)機(jī)氮的濃度迅速升高，氮限制作用大幅減弱，不再成為限制因子，而無(wú)機(jī)磷的來(lái)源相對(duì)較少，增速慢于無(wú)機(jī)氮。以膠州灣中東部水域?yàn)槔?960～1980年代，無(wú)機(jī)氮增加7.2倍，而無(wú)機(jī)磷僅增加2.2倍。與此同時(shí)，由于河流上游的水利工程，大沽河的徑流量和輸沙率呈下降的趨勢(shì)，大量泥沙被截留。1960年代以前，大沽河年輸沙量平均為170.03萬(wàn)t，到1970年代末僅為71.21萬(wàn)t，減少59%，導(dǎo)致硅酸鹽濃度顯著降低，成為限制浮游植物的主要因子，僅在汛期，大沽河的徑流輸入帶來(lái)大量的硅酸鹽，此時(shí)無(wú)機(jī)磷相對(duì)不足，呈現(xiàn)磷限制[33]。本世紀(jì)初，氮和磷的限制作用均大大減弱，浮游植物的生長(zhǎng)受硅限制作用明顯，近年來(lái)大沽河輸沙量減少為36.59萬(wàn)t，僅為1960年代以前的22%，導(dǎo)致膠州灣內(nèi)硅營(yíng)養(yǎng)鹽濃度的不斷降低，進(jìn)而成為浮游植物生物量的主要限制因子。

比較不同時(shí)代營(yíng)養(yǎng)鹽綜合限制因子的數(shù)值可以發(fā)現(xiàn)，近40年來(lái)該值逐漸減小，說(shuō)明膠州灣內(nèi)的營(yíng)養(yǎng)鹽限制作用有不斷增強(qiáng)的趨勢(shì)。圖9為模擬得到的1960年代、1980年代和本世紀(jì)初的膠州灣浮游植物生物量的季節(jié)變化。相對(duì)于1960年代，1980年代略有降低，而21世紀(jì)初浮游植物生物量大約減少23%。由于本模型只考慮一種浮游植物，只能體現(xiàn)浮游植物生物量的變化，無(wú)法反映浮游植物群落結(jié)構(gòu)的變化。但是，自1960年代以來(lái)，隨著無(wú)機(jī)氮和無(wú)機(jī)磷輸入的增加以及硅酸鹽輸入的減少，限制性營(yíng)養(yǎng)鹽發(fā)生轉(zhuǎn)變，可以推測(cè)膠州灣內(nèi)的大型浮游植物植物不斷減少，而小型浮游植物生物量逐漸增加，并導(dǎo)致總體浮游植物生物量顯著減少。

圖8 營(yíng)養(yǎng)鹽限制因子的季節(jié)變化

圖9 不同年代浮游植物生物量的季節(jié)變化特征

4 結(jié)語(yǔ)

本研究建立膠州灣箱式生態(tài)動(dòng)力學(xué)模型，模擬膠州灣內(nèi)不同區(qū)域及不同年代營(yíng)養(yǎng)鹽限制浮游植物生長(zhǎng)情況的差異。膠州灣內(nèi)營(yíng)養(yǎng)鹽限制作用有顯著的季節(jié)變化，春季最強(qiáng)，夏季最弱，并存在區(qū)域差異，河口區(qū)和灣口區(qū)的浮游植物生長(zhǎng)受硅和磷的交替限制，養(yǎng)殖區(qū)呈現(xiàn)顯著的硅限制，進(jìn)而導(dǎo)致浮游植物生物量的空間分布不同。河口區(qū)的營(yíng)養(yǎng)鹽限制作用最弱，浮游植物生物量最高，養(yǎng)殖區(qū)次之，而灣口區(qū)營(yíng)養(yǎng)鹽限制作用最強(qiáng)，浮游植物生物量低于全灣平均水平。

從時(shí)間演變上看，1960年代，膠州灣內(nèi)營(yíng)養(yǎng)鹽限制作用以氮限制為主；1980年代，磷和硅交替限制浮游植物的生長(zhǎng)，夏季為磷限制，其他季節(jié)為硅限制；本世紀(jì)全年呈現(xiàn)顯著的硅限制。近40年來(lái)，營(yíng)養(yǎng)鹽限制作用有逐漸增強(qiáng)的趨勢(shì)，受此影響，膠州灣內(nèi)的浮游植物生物量也不斷下降，相對(duì)20世紀(jì)60年代，21世紀(jì)初浮游植物生物量降低約23%。營(yíng)養(yǎng)鹽限制作用的轉(zhuǎn)變可能會(huì)導(dǎo)致藻類(lèi)優(yōu)勢(shì)種的變化，并通過(guò)食物鏈向上層營(yíng)養(yǎng)級(jí)傳遞，最終會(huì)影響膠州灣內(nèi)海洋生態(tài)結(jié)構(gòu)和漁業(yè)生產(chǎn)。

在近岸海域，尤其是膠州灣這樣的半封閉海灣，其營(yíng)養(yǎng)鹽限制作用會(huì)隨著營(yíng)養(yǎng)鹽來(lái)源的變化而改變，并會(huì)對(duì)生態(tài)系統(tǒng)產(chǎn)生影響。膠州灣與世界上其他類(lèi)似海域相比，有其獨(dú)特性，但也存在許多共性。Danielsson等[5]研究發(fā)現(xiàn)，在波羅的海內(nèi)，由于不同區(qū)域營(yíng)養(yǎng)鹽來(lái)源的差異，營(yíng)養(yǎng)鹽限制作用也存在顯著的區(qū)域差異，其中的芬蘭灣，其營(yíng)養(yǎng)鹽來(lái)源與膠州灣相似，內(nèi)有較大河流的輸入，營(yíng)養(yǎng)鹽限制作用在夏季、秋初表現(xiàn)為氮限制，而其余時(shí)間均處于硅限制，這也與膠州灣營(yíng)養(yǎng)鹽限制作用的季節(jié)演替規(guī)律相似。而在沒(méi)有河流來(lái)源的波的尼亞灣，全年均呈現(xiàn)顯著的磷限制，沒(méi)有季節(jié)更替。從長(zhǎng)期來(lái)看，波羅的海有逐步向硅限制轉(zhuǎn)化的潛在趨勢(shì)，而且這種轉(zhuǎn)化已經(jīng)其浮游植物群落結(jié)構(gòu)和初級(jí)生產(chǎn)力產(chǎn)生影響。同樣，在地中海的河口區(qū)域，長(zhǎng)期以來(lái)氮和磷增加顯著，而硅的變化較小，也促進(jìn)了該水域有害鞭毛藻的爆發(fā)[7]，而膠州灣的長(zhǎng)期變化趨勢(shì)也與之類(lèi)似，大量氮、磷的增加以及硅的減少，導(dǎo)致膠州灣內(nèi)硅藻等大型藻類(lèi)正在減少，取而代之的是不需要吸收硅的甲藻和鞭毛藻[18]?，F(xiàn)階段我們的模型僅能模擬出營(yíng)養(yǎng)鹽組成及浮游植物生物量的變化特征，下一步工作中，作者將在模型中增加不同浮游植物類(lèi)型，綜合討論人類(lèi)活動(dòng)和氣候變化影響下?tīng)I(yíng)養(yǎng)鹽和浮游植物群落結(jié)構(gòu)的演變。

[1] Wu J T, Chou, T L. Silicate as the limiting nutrient for phytoplankton in a subtropical eutrophic estuary of Taiwan [J]. Estuarine Coastal and Shelf Science, 2003, 58： 155-162.

[2] Redfield A, Ketchum B, Richards F. 1963. The influence of organisms on the composition of seawater [M]. //Hill M, ed. The Sea. New York： Wiley, 1963： 26-77.

[3] Tyrerell T. The relative influences of nitrogen and phosphorus on oceanic primary production [J]. Nature， 1999， 400(6744): 525-531.

[4] Wan Z, Bi H, She J, et al. Model study on horizontal variability of nutrient N/P ratio in the Baltic Sea and its impacts on primary production, nitrogen fixation and nutrient limitation [J]. Ocean Science Discussions, 2012(9): 385-419.

[5] Danielsson ?, Papush L, Rahm L. Alterations in nutrient limitations—Scenarios of a changing Baltic Sea [J]. Journal of Marine Systems, 2008, 73(3): 263-283.

[6] Fisher T, Hagy J, Boynton W, et al. Cultural eutrophication in the Choptank and Patuxent estuaries of Chesapeake Bay [J]. Limnology and Oceanography, 2006, 51(1): 435-447.

[7] Domingues R, Anselmo T, Barbosa A, et al. Nutrient limitation of phytoplankton growth in the freshwater tidal zone of a turbid, Mediterranean estuary [J]. Estuarine, Coastal and Shelf Science, 2010, 91(2): 282-297.

[8] Liu Z, Wei H, Liu G L, et al. Simulation of water exchange in Jiaozhou Bay by average residence time approach [J]. Estuarine, Coastal and shelf Science, 2004, 61(1): 5-35.

[9] 張哲, 王江濤. 膠州灣營(yíng)養(yǎng)鹽研究概述 [J]. 海洋科學(xué), 2009, 33(11): 90-94.

[10] 李莉， 梁生康， 石曉勇， 等. 2007年膠州灣入海河流污染狀況和污染物入海通量分析 [J]. 環(huán)境科學(xué)與管理, 2009, 34(6): 23-28.

[11] 宋秀賢, 俞志明. 膠州灣東北部養(yǎng)殖海域夏季營(yíng)養(yǎng)鹽分布特征及其對(duì)浮游植物生長(zhǎng)的影響 [J]. 海洋與湖沼, 2007, 38(5): 446-452.

[12] 任玲, 張曼平, 孫軍, 等. 膠州灣內(nèi)外水體部分化學(xué)生物因子調(diào)查分析 [J]. 青島海洋大學(xué)學(xué)報(bào)： 自然科學(xué)版, 2003, 33(4): 557-564.

[13] Yang D F, Gao Z H, Sun P Y, et al. Mechanism of nutrient silicon and water temperature influences on phytoplankton growth [J]. Marine Science Bulletin, 2006, 8(2): 49-59.

[14] 李艷, 李瑞香, 王宗靈, 等. 膠州灣浮游植物群落結(jié)構(gòu)及其變化的初步研究 [J]. 海洋科學(xué)進(jìn)展, 2005, 23(3): 328-334.

[15] 李超倫, 張芳, 申欣, 等. 膠州灣葉綠素的濃度、分布特征及周年變化 [J]. 海洋與湖沼, 2005, 36(6): 499-506.

[16] 孫曉霞， 孫松， 趙增霞， 等. 膠州灣營(yíng)養(yǎng)鹽濃度與結(jié)構(gòu)的長(zhǎng)期變化 [J]. 海洋與湖沼, 2011, 42(5): 662-669.

[17] 孫曉霞， 孫松， 張永山， 等. 膠州灣葉綠素a及初級(jí)生產(chǎn)力的長(zhǎng)期變化 [J]. 海洋與湖沼, 2011, 42(5): 654-661.

[18] Shen Z L. Historical changes in nutrient structure and its influences on phytoplantkon composition in Jiaozhou Bay [J]. Estuarine, Coastal and Shelf Science， 2001， 52(2): 211-224.

[19] 王勇， 趙澎， 單寶田. 膠州灣營(yíng)養(yǎng)鹽限制浮游植物生長(zhǎng)的初步模擬現(xiàn)場(chǎng)實(shí)驗(yàn)研究 [J]. 海洋科學(xué), 2002, 26(10): 55 -58.

[20] 姚云， 鄭世清， 沈志良. 膠州灣營(yíng)養(yǎng)鹽及富營(yíng)養(yǎng)化特征 [J]. 海洋通報(bào), 2007, 26(4): 91-98.

[21] 李克強(qiáng)， 王修林， 石曉勇， 等. 膠州灣圍隔浮游生態(tài)系統(tǒng)氮、磷營(yíng)養(yǎng)鹽遷移-轉(zhuǎn)化模型研究 [J]. 海洋學(xué)報(bào), 2007, 29(3): 76-83.

[22] Han H Y, Li K Q, Wang X L, et al. Environmental capacity of nitrogen and phosphorus pollutions in Jiaozhou Bay, China: Modeling and assessing [J]. Marine Pollution Bulletin, 2011, 63(5): 262-266.

[23] Ren L, Zhang M P, Brockmann, U. H, et al. Pelagic Nitrogen Cycling in Jiaozhou Bay, a model study Ⅰ: The conceptual model [J]. Chinese Journal of Oceanology and Limnology, 2003b, 21, 358-367.

[24] Liu Z, Wei H, Bai J, et al. Nutrients seasonal variation and budget in Jiaozhou Bay, China: a 3-dimensional physical-biological coupled model study [J]. Water Air & Soil Pollution: Focus, 2007, 7(6): 607-623.

[25] 王震勇. 膠州灣浮游生態(tài)系統(tǒng)四十年變化的模擬與分析 [D]. 青島： 中國(guó)海洋大學(xué), 2007.

[26] 張燕, 孫英蘭, 袁道偉, 等. 膠州灣氮、磷濃度的三維數(shù)值模擬 [J]. 中國(guó)海洋大學(xué)學(xué)報(bào)： 自然科學(xué)版, 2007, 31(1): 021-026.

[27] 趙亮, 魏皓. 膠州灣水交換的數(shù)值研究 [J]. 海洋與湖沼, 2002, 33(1): 23-29.

[28] 浦祥. 大沽河流域非點(diǎn)源氮輸出及其對(duì)膠州灣水質(zhì)影響的數(shù)值研究 [D]. 青島: 中國(guó)海洋大學(xué), 2011.

[29] Liu S M, Zhang J, Chen H T, et al. Factors influencing nutrient dynamics in the eutrophic Jiaozhou Bay, North China [J]. Progress in Oceanography, 2005, 66(1): 66-85.

[30] 王剛. 膠州灣入海點(diǎn)源、海水養(yǎng)殖污染物通量研究 [D]. 青島： 中國(guó)海洋大學(xué), 2009.

[31] Liu S M, Li X N, Zhang J, et al. Nutrient dynamics in Jiaozhou Bay [J]. Water Air Soil Pollution: Focus, 2007, 7(6): 625-643.

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

[33] 王保棟. 黃海和東海營(yíng)養(yǎng)鹽分布及其對(duì)浮游植物的限制 [J]. 青島海洋大學(xué)學(xué)報(bào)： 自然科學(xué)版, 2003, 14(7): 1122-1126.

[34] 吳玉霖， 孫松， 張永山， 等. 膠州灣浮游植物數(shù)量長(zhǎng)期動(dòng)態(tài)變化的研究 [J]. 海洋與湖沼, 2004, 35(6): 518-523.

[35] 沈志良. 膠州灣營(yíng)養(yǎng)鹽結(jié)構(gòu)的長(zhǎng)期變化及其對(duì)生態(tài)環(huán)境的影響 [J]. 海洋與湖沼, 2002, 33(3): 322-331.

[36] 王保棟. 黃海和東海營(yíng)養(yǎng)鹽分布及其對(duì)浮游植物的限制 [J]. 青島海洋大學(xué)學(xué)報(bào)： 自然科學(xué)版, 2003, 14(7): 1122-1126.

[37] Zhang J. Nutrient elements in large Chinese estuaries [J]. Continental Shelf Research, 1996, 16(8): 1023-1045.

責(zé)任編輯 龐 旻

Numerical Study on the Regional and Long-Term Variation of Nutrient Limitation in Jiaozhou Bay

JIA Shou-Wei， SHI Jie， GAO Hui-Wang

(The Key Laboratory of Marine Environment and Ecology，Ministry of Education，Ocean University of China，Qingdao 266100，China)

A box ecosystem model is developed to simulate the seasonal variations of nutrient limitation in the Jiaozhou Bay.After model validation, the regional and interannual variations of nutrient limitation are analyzed. The results show that the nutrient limitation exhibits obvious seasonal variation. In spring, the phytoplankton grows fast and exhausts the available nutrients. The nutrient limitation becomes the strongest in late spring and turns to be the weakest in summer when the river in flood season supply sufficient nutrients into the Jiaozhou Bay. Meanwhile, the regional variations of nutrient limitation are obvious. The limitation nutrient is SIL in aquaculture area, and is SIL and DIP alternatively in both the river mouth area and the mouth of the Bay. The nutrient limitation is the strongest in the mouth area of the Bay, while the weakest in the river mouth area. In 1960s, the main limiting nutrient is DIN. Under the influence of human activity and climate change, the limiting nutrient has altered from DIN to SIL until now. The DIP limitation occurs in the certain region and period. The nutrient limitation is continually enhanced as time goes on. As a result, the phytoplankton biomass at the beginning of this country decreases by 21% comparing with that in 1960s.

the Jiaozhou Bay; nutrient limitation; regional variation; interannual variation

國(guó)家自然科學(xué)基金青年基金項(xiàng)目(41106007)；國(guó)家自然科學(xué)基金委員會(huì)—山東省人民政府聯(lián)合資助海洋科學(xué)研究中心項(xiàng)目(N1406403)資助

2014-01-17；

2014-02-27

賈守偉(1987-)，男，碩士生。

** 通訊作者： E-mail：shijie@ouc.edu.cn

P735

A

1672-5174(2015)05-001-10

10.16441/j.cnki.hdxb.20140024