膠州灣營養鹽限制作用空間差異和長期變化的數值研究*

賈守偉， 史 潔， 高會旺

(中國海洋大學海洋環境與生態教育部重點實驗室，山東 青島 266100)

?

膠州灣營養鹽限制作用空間差異和長期變化的數值研究*

賈守偉， 史 潔**， 高會旺

(中國海洋大學海洋環境與生態教育部重點實驗室，山東 青島 266100)

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

膠州灣； 營養鹽限制； 區域差異； 長期變化

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

膠州灣是受人類活動影響顯著的半封閉海灣，外與黃海相通，水體平均存留時間為52d[8]，其營養鹽結構在北方海灣中具有較強的代表性[9]。膠州灣有十余條入海河流，其中位于西北部的大沽河占入海總流量的70%[10]，是膠州灣重要的營養鹽來源，其帶來的營養鹽直接影響河口區域。膠州灣內水產養殖業發達，養殖區多集中在東北部，是典型的對蝦和蛤蜊養殖區，餌料的過量投放和養殖生物的排泄為該區域輸送大量的氮、磷營養鹽，使該區域內這兩種營養鹽濃度遠遠高于全灣平均水平[11]。膠州灣內的營養鹽濃度常年高于外海，通過灣口處與外海的交換，大量營養鹽被輸送到黃海，因此灣口區域的營養鹽結構也有別于其他區域[12]。

現有研究表明，在膠州灣內，相較于光照和水溫，營養鹽對浮游植物生長的限制作用更為重要[9]。基于觀測資料，已有研究闡明了膠州灣內營養鹽濃度及浮游植物生物量的季節、長期變化和空間分布特征[12,14-18]。王勇等[19]在膠州灣通過營養鹽添加實驗發現，網采浮游植物對鐵或氮、磷的添加有比較明顯的響應。姚云[20]通過分析氮磷比探討了膠州灣內限制性營養鹽的種類及季節變化，認為相較于氮和磷而言，帶有硅的評價結果能夠更好地反映膠州灣的富營養化狀況。觀測和實驗能夠反映膠州灣營養鹽結構的現狀，但卻存在著時間不連續性，并難以區分各個營養鹽來源對膠州灣營養鹽結構的影響。海洋生態動力學模型是定量研究生態系統中各過程相互作用的有效方法，不僅可以反演現場觀測結果，還可以進行機制分析，模擬和預測生態系統的變化。現有的膠州灣生態系統動力學模型研究能夠模擬膠州灣營養鹽濃度和浮游植物生物量的季節變化，計算得到膠州灣內的營養鹽收支狀況，但是沒有分析不同營養鹽來源對膠州灣營養鹽限制情況的影響，缺少源的長期變化導致膠州灣生態系統演變的研究[21-24]。

本文選取膠州灣為研究海域(見圖1)，建立膠州灣箱式生態動力學模型，區分河口區、養殖區和灣口區3個典型區域，模擬并探討膠州灣內營養鹽限制作用的區域差異，分析膠州灣不同營養鹽來源對灣內營養鹽限制作用以及浮游植物生物量的影響，模擬并分析了膠州灣營養鹽結構隨營養鹽來源變化而產生的長期變化趨勢。

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

1 研究方法

1.1 模型介紹

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

模型的控制方程如下：

(1)

(2)

(3)

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

浮游植物呼吸作用分為光呼吸和基礎呼吸過程，并考慮溫度對呼吸作用的影響。

(4)

rT=er·(T-T0)

(5)

浮游植物的生長受光照條件、溫度、營養鹽的共同影響，其表達式為：

(6)

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

溫度對浮游植物生長的限制采用Q10法則，即:

(7)

光對浮游植物的生長有重要影響，其限制作用的表達函數為:

(8)

其中：Iopt為最優光強；I為水體有效太陽輻射。

本研究中采用Michaeli-Menten(M-M)關系來量化浮游植物生長受到環境中營養鹽限制作用的強弱，定義營養鹽限制因子為：

(9)

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

VN=min(VDIN，VDIP，VSIL)

(10)

模型中參數的取值(見表1)采用前人在膠州灣及中國近海的生態模型研究結果[24-26]。

表1 模式中主要參數說明及取值Table 1 Values of parameters used in the model

1.2 模型設置和數據來源

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

模型中的物理強迫包括海水溫度和海面光強，二者均用余弦函數來擬合，模型采用以下公式進行計算[23]：

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

(11)

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

(12)

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

模型采用冬季觀測的平均濃度作為初始場，無機氮、無機磷和硅酸鹽的濃度分別設為10、0.5、3.5μmol/L，葉綠素濃度為0.8mg/m3[14-15]。開始時間為1月1日，模型時間步長為1h，模型在以上強迫下連續運行3年，用第3年的結果進行分析。

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

1960年代1960s1980年代1980s21世紀初2000s單位Unit來源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]養殖排放氮⑤00.6854106mol[20]養殖排放磷⑥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 結果

2.1 敏感性分析

在一個模式中，參數的變化對模擬結果的影響程度反映了模式結果對該參數的敏感程度。由于生物過程本身的非線性特征，生態變量的計算值會對模型中一些參數的變化非常敏感。因此，參數的敏感性分析是生態模型研究中必不可少的步驟。通過敏感性實驗，找出生態系統變化敏感參數，為了比較不同參數的敏感程度，定義敏感度為：

(12)

其中：ΔF/F為生態模型中浮游植物生物量的變化率；Δα/α為模型中與生態變量F有關的某個生物參數的變化率。表1提供了10個參數的敏感性數值，浮游植物生物量對最大生長率、浮游植物基礎呼吸率、最優光強等參數較敏感，而對礦化率、氮營養鹽半飽和常數等參數相對不敏感。各參數分別變化或減小相同比例，其敏感度并不相同，表明各參數的變化是非線性變化，如浮游植物死亡率，增大50%其敏感度為0.65，減小50%其敏感度為1.91。75%的參數敏感度小于0.5，因此本文建立的模型是穩定的。

2.2 本世紀初營養鹽濃度和浮游植物生物量的季節變化

應用所建立的膠州灣箱式生態動力學模型，模擬本世紀初的膠州灣營養鹽濃度和浮游植物生物量的季節變化規律(見圖4)。

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

圖4 營養鹽濃度和浮游植物生物量的季節變化

2.3 本世紀初營養鹽限制的季節變化

3種營養鹽對浮游植物的限制作用存在強弱差異，因此探討營養鹽限制作用的季節變化對于理解營養鹽結構對浮游植物生長的影響有重要作用。通過比較模擬結果中3種營養鹽限制因子的數值大小，來確定限制浮游植物生長的營養鹽種類，即營養鹽限制因子數值最小的被認為是限制浮游植物生長的營養鹽。

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

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

圖5 營養鹽限制因子的季節變化Fig.5 Seasonal variations of nutrient limiting factors

3 討論

3.1 營養鹽限制的區域差異

膠州灣內不同區域營養鹽的來源不同，并且各個營養鹽來源的量和比例也有較大區別，因此營養鹽濃度及結構都存在著空間不均勻性，浮游植物生長受到的營養鹽限制作用也會存在較大的區域差異。

表3 不同區域的營養鹽來源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

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

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

圖6 不同區域營養鹽限制因子的季節變化

圖7 不同區域浮游植物生物量的季節變化

3.2 營養鹽限制的年際差異

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

表4 營養鹽濃度長期變化模擬結果與觀測結果對比[35]Table 4 Comparison between the simulation and observation in the 1960s and the 1980s

Note:①Simulation；②Observation

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

比較不同時代營養鹽綜合限制因子的數值可以發現，近40年來該值逐漸減小，說明膠州灣內的營養鹽限制作用有不斷增強的趨勢。圖9為模擬得到的1960年代、1980年代和本世紀初的膠州灣浮游植物生物量的季節變化。相對于1960年代，1980年代略有降低，而21世紀初浮游植物生物量大約減少23%。由于本模型只考慮一種浮游植物，只能體現浮游植物生物量的變化，無法反映浮游植物群落結構的變化。但是，自1960年代以來，隨著無機氮和無機磷輸入的增加以及硅酸鹽輸入的減少，限制性營養鹽發生轉變，可以推測膠州灣內的大型浮游植物植物不斷減少，而小型浮游植物生物量逐漸增加，并導致總體浮游植物生物量顯著減少。

圖8 營養鹽限制因子的季節變化

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

4 結語

本研究建立膠州灣箱式生態動力學模型，模擬膠州灣內不同區域及不同年代營養鹽限制浮游植物生長情況的差異。膠州灣內營養鹽限制作用有顯著的季節變化，春季最強，夏季最弱，并存在區域差異，河口區和灣口區的浮游植物生長受硅和磷的交替限制，養殖區呈現顯著的硅限制，進而導致浮游植物生物量的空間分布不同。河口區的營養鹽限制作用最弱，浮游植物生物量最高，養殖區次之，而灣口區營養鹽限制作用最強，浮游植物生物量低于全灣平均水平。

從時間演變上看，1960年代，膠州灣內營養鹽限制作用以氮限制為主；1980年代，磷和硅交替限制浮游植物的生長，夏季為磷限制，其他季節為硅限制；本世紀全年呈現顯著的硅限制。近40年來，營養鹽限制作用有逐漸增強的趨勢，受此影響，膠州灣內的浮游植物生物量也不斷下降，相對20世紀60年代，21世紀初浮游植物生物量降低約23%。營養鹽限制作用的轉變可能會導致藻類優勢種的變化，并通過食物鏈向上層營養級傳遞，最終會影響膠州灣內海洋生態結構和漁業生產。

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

[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] 張哲, 王江濤. 膠州灣營養鹽研究概述 [J]. 海洋科學, 2009, 33(11): 90-94.

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

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

[12] 任玲, 張曼平, 孫軍, 等. 膠州灣內外水體部分化學生物因子調查分析 [J]. 青島海洋大學學報： 自然科學版, 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] 李艷, 李瑞香, 王宗靈, 等. 膠州灣浮游植物群落結構及其變化的初步研究 [J]. 海洋科學進展, 2005, 23(3): 328-334.

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

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

[17] 孫曉霞， 孫松， 張永山， 等. 膠州灣葉綠素a及初級生產力的長期變化 [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] 王勇， 趙澎， 單寶田. 膠州灣營養鹽限制浮游植物生長的初步模擬現場實驗研究 [J]. 海洋科學, 2002, 26(10): 55 -58.

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

[21] 李克強， 王修林， 石曉勇， 等. 膠州灣圍隔浮游生態系統氮、磷營養鹽遷移-轉化模型研究 [J]. 海洋學報, 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] 王震勇. 膠州灣浮游生態系統四十年變化的模擬與分析 [D]. 青島： 中國海洋大學, 2007.

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

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

[28] 浦祥. 大沽河流域非點源氮輸出及其對膠州灣水質影響的數值研究 [D]. 青島: 中國海洋大學, 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] 王剛. 膠州灣入海點源、海水養殖污染物通量研究 [D]. 青島： 中國海洋大學, 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] 田恬, 魏皓, 蘇健, 等. 黃海氮磷營養鹽的循環和收支研究 [J]. 海洋科學進展, 2003, 21(1): 1-11.

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

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

[35] 沈志良. 膠州灣營養鹽結構的長期變化及其對生態環境的影響 [J]. 海洋與湖沼, 2002, 33(3): 322-331.

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

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

責任編輯 龐 旻

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

國家自然科學基金青年基金項目(41106007)；國家自然科學基金委員會—山東省人民政府聯合資助海洋科學研究中心項目(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