桑溝灣養(yǎng)殖生態(tài)系統(tǒng)健康綜合評價

傅明珠，蒲新明，* ，王宗靈，劉新杰

(1.國家海洋局第一海洋研究所，青島 266061;2.海洋生態(tài)環(huán)境科學與工程國家海洋局重點實驗室，青島 266061;3.山東省榮成市海洋與漁業(yè)局，榮成 264300)

桑溝灣是以筏式貝藻大規(guī)模養(yǎng)殖為主要特征的典型海灣，養(yǎng)殖種類、方式與布局具有自身的特征，年產(chǎn)貝類(扇貝、牡蠣)、大型藻數(shù)萬噸，是我國北方最具有代表性的重點養(yǎng)殖海區(qū)和海產(chǎn)品基地。在桑溝灣養(yǎng)殖業(yè)的發(fā)展過程中，養(yǎng)殖企業(yè)為追求高產(chǎn)量、高產(chǎn)出，在提高養(yǎng)殖密度的同時，亦大規(guī)模擴展養(yǎng)殖面積，海帶養(yǎng)殖已擴展到灣外20 m水深海域。隨著養(yǎng)殖規(guī)模的不斷擴大，該海灣近來出現(xiàn)了海帶在收獲季節(jié)開始腐爛，櫛孔扇貝養(yǎng)殖個體小型化、養(yǎng)殖周期延長、死亡率升高、產(chǎn)品質(zhì)量下降等問題，嚴重影響了該灣海水養(yǎng)殖業(yè)的可持續(xù)發(fā)展［1］。

關于桑溝灣的生態(tài)環(huán)境問題已有較多研究，例如水體的營養(yǎng)狀況［2-4］、浮游植物多樣性的長期變化［5］、貝類養(yǎng)殖對水動力條件［6-7］和浮游與底棲環(huán)境的影響［8-9］等，然而將生態(tài)系統(tǒng)作為一個整體進行綜合健康評價的研究還未見報道。為了量化桑溝灣養(yǎng)殖生態(tài)系統(tǒng)受損和退化的程度，本文應用《海水養(yǎng)殖生態(tài)系統(tǒng)健康綜合評價:方法與模式》一文所建立的方法，對桑溝灣的環(huán)境壓力、生物群落結(jié)構、生態(tài)系統(tǒng)功能等方面進行系統(tǒng)地定性與定量評價，以期真實、客觀的反映該海灣養(yǎng)殖生態(tài)系統(tǒng)的現(xiàn)狀及可能的變化趨勢，并提出可持續(xù)發(fā)展管理的建議，為桑溝灣養(yǎng)殖生態(tài)系統(tǒng)退化機制與修復技術的研究提供科學依據(jù)和技術支持。

1 材料和方法

1.1 桑溝灣海水養(yǎng)殖概況

桑溝灣位于山東半島東端，榮成市境內(nèi)(122°24'—122°35'E，37°01'—37°09'N)，北、西、南三面為陸地環(huán)抱，灣口朝東面向黃海(圖1)，為半封閉海灣。該灣南北口長約11.5 km，東西寬約7.5 km，灣內(nèi)面積約133.33 km2，灣內(nèi)平均水深7—8 m，最大水深15—17 m，呈自西向東遞增的趨勢。

桑溝灣養(yǎng)殖業(yè)十分發(fā)達，養(yǎng)殖面積幾乎覆蓋全灣，成為我國黃海沿岸的重要水產(chǎn)養(yǎng)殖水域。該灣的主要養(yǎng)殖類型由灣口向灣內(nèi)依次為海帶養(yǎng)殖區(qū)、貝藻混養(yǎng)區(qū)和貝類養(yǎng)殖區(qū)(圖1)，傳統(tǒng)的魚類網(wǎng)箱養(yǎng)殖在灣西南角也有分布。桑溝灣養(yǎng)殖業(yè)始于20世紀60年代末，當時海帶是唯一的養(yǎng)殖品種，70年代曾一度養(yǎng)殖貽貝，80年代隨著扇貝人工育苗的成功，櫛孔扇貝養(yǎng)殖得到迅速發(fā)展，90年代末期由于櫛孔扇貝夏季大規(guī)模死亡，養(yǎng)殖者轉(zhuǎn)養(yǎng)牡蠣，目前桑溝灣以牡蠣和海帶養(yǎng)殖為主［1，10］(圖2)，魚類的網(wǎng)箱養(yǎng)殖規(guī)模近年也有所擴大。

1.2 數(shù)據(jù)來源

對桑溝灣健康現(xiàn)狀的評價主要采用了2006—2007年度908專項的調(diào)查數(shù)據(jù)，桑溝灣養(yǎng)殖數(shù)據(jù)來自榮成市海洋與漁業(yè)局提供的內(nèi)部資料。歷史變化趨勢的比較使用了國家海洋局第一海洋所的 1983—1984、1988—1989、1999—2000年的歷史調(diào)查資料，其它過程參數(shù)以及908生態(tài)調(diào)查沒有涉及的指標引用了文獻中發(fā)表的數(shù)據(jù)，詳見下文計算與評價過程。

1.3 評價方法與數(shù)據(jù)處理

對桑溝灣生態(tài)系統(tǒng)健康狀況的評價，采用《海水養(yǎng)殖生態(tài)系統(tǒng)健康綜合評價Ⅰ:方法與模式》一文中建立的指標框架和評分標準。對于評價過程中遇到的幾種特殊情況按以下方式處理:

(1)若某項指標沒有同期調(diào)查數(shù)據(jù)，則采用歷史上該參數(shù)的多年平均值進行評價得分;(2)若某項指標沒有歷史調(diào)查數(shù)據(jù)，則將該指標的權重按比例分配至與其同層次的其他指標中;(3)對每個站的指標分別進行評價得分，然后將各站點的得分進行平均得到整個研究區(qū)域的評價得分;(4)若同一指標有多個參數(shù)，則取得分最低的參數(shù)作為該項的得分。

圖1 桑溝灣2006—2007年生態(tài)調(diào)查站位及養(yǎng)殖布局示意圖Fig.1 Diagram of aquaculture distribution and investigation stations in Sanggou Bay during 2006—2007

2 結(jié)果與討論

2.1 壓力分析

2.1.1 外源污染壓力

桑溝灣入灣河流有桑干河、崖頭河、沽河、小落河等，陸源物質(zhì)來源較少。陸源物質(zhì)主要隨工業(yè)廢水和生活污水排入海中。位于灣西北部的沽河口為城市污水和工業(yè)廢水進入桑溝灣的入口(圖1)，約占該海域等標污染負荷比的99.03%［11］。因此本文以沽河對桑溝灣的污染物輸入為基礎，計算該灣的外源污染壓力。

根據(jù)海洋養(yǎng)殖生態(tài)系統(tǒng)評價方法，營養(yǎng)鹽壓力采用以下污染指數(shù)進行計算:

各參數(shù)的選擇和計算方法如下:沽河平均徑流量(6.47×107m3/a)引自張朝暉等［12］的實測數(shù)據(jù)，沽河DIN平均濃度(0.36 mol/m3)與PO4-P平均濃度(0.015 mol/m3)采用榮成市環(huán)保局提供的數(shù)據(jù);沽河營養(yǎng)鹽入海通量Fi用徑流平均濃度與沽河平均徑流量乘積得到;桑溝灣水體中營養(yǎng)鹽濃度Ci采用2006—2007年908專項調(diào)查數(shù)據(jù)的年度平均值;桑溝灣水體體積Vi用海域面積(133.33 km2)與平均水深(8 m)的乘積得到，計算結(jié)果見表1。

表1 桑溝灣營養(yǎng)鹽污染壓力評價參數(shù)及結(jié)果Table 1 Parameters in the assessment of nutrient pressure

根據(jù)海水養(yǎng)殖生態(tài)系統(tǒng)外源污染壓力評價得分標準，DIN的污染指數(shù)Pi在0.5—2之間，得分為4，PO4-P的污染指數(shù)Pi在2—5之間，得分為3，取最低評價得分3作為營養(yǎng)鹽污染壓力指數(shù)。評價結(jié)果表明桑溝灣營養(yǎng)鹽污染壓力較大。

由于沒有查到桑溝灣重金屬和石油烴陸源輸入通量數(shù)據(jù)，本次評價將二者的權重分配至營養(yǎng)鹽污染壓力指標中，即營養(yǎng)鹽污染壓力指標權重變?yōu)?4.5，因此本模塊綜合評價得分即為營養(yǎng)鹽壓力得分=3分。

2.1.2 養(yǎng)殖壓力

本文以桑溝灣的主要養(yǎng)殖品種海帶、扇貝、鮑魚、牡蠣進行養(yǎng)殖壓力評價。

(1)養(yǎng)殖密度/規(guī)模

根據(jù)榮成市海洋與漁業(yè)局提供的資料，桑溝灣1999—2009年主要養(yǎng)殖品種的養(yǎng)殖面積以及不同養(yǎng)殖品種的比例見圖2。近十年來，海帶始終是桑溝灣的主要養(yǎng)殖品種，養(yǎng)殖面積變化不大，扇貝的養(yǎng)殖面積出現(xiàn)下降趨勢，而牡蠣養(yǎng)殖面積呈增加趨勢。為了與環(huán)境質(zhì)量模塊各指標時間一致(908調(diào)查為2006—2007年)，本文采用2006年和2007年桑溝灣養(yǎng)殖面積的平均值進行養(yǎng)殖密度的評價。

對養(yǎng)殖密度/規(guī)模的評價采用養(yǎng)殖面積占海灣面積的比例進行計算，2006年和2007年桑溝灣平均養(yǎng)殖面積為155km2。需要注意的是，對于筏式養(yǎng)殖，1養(yǎng)殖畝實際占用水面面積并不等于1實測畝。在榮成，每666.7m25 道梗，筏間距為4—5 m，每養(yǎng)殖 666.7m2占水面實測2.5×666.7m2［13］。本文一律將養(yǎng)殖畝換算成實測面積61.78×2.5=155 km2，超過桑溝灣水域面積133.33 km2。這可能是由于桑溝灣的養(yǎng)殖海域已經(jīng)擴展到灣外的原因，這從桑溝灣的衛(wèi)星圖片上也可以清晰的看出。

張繼紅［14］論文中指出，桑溝灣現(xiàn)在的養(yǎng)殖面積約占海域的70%，另外根據(jù)2008年和2009年山東省海洋環(huán)境質(zhì)量公報［15-16］，桑溝灣養(yǎng)殖面積為1萬hm2(約占海域的75%)，可以認為桑溝灣養(yǎng)殖面積占水域面積在70%—80%之間，根據(jù)評價標準，養(yǎng)殖密度/規(guī)模評價得分為2分。

圖2 桑溝灣養(yǎng)殖面積(1999—2009)與產(chǎn)量(1980—2009)變化趨勢Fig，2 Variations trend of aquaculture area(1999—2009)and yield(1980—2009)in Sanggou Bay

(2)養(yǎng)殖自身污染

桑溝灣的海帶養(yǎng)殖和貝類養(yǎng)殖均屬于自養(yǎng)生態(tài)系統(tǒng)，不需要投放餌料。貝類的生物沉積作用和養(yǎng)殖筏架對海流的阻礙作用，是形成養(yǎng)殖自身污染的兩個重要因素。本文主要考慮養(yǎng)殖生物對海區(qū)所產(chǎn)生的營養(yǎng)鹽負荷的污染(主要是影響較大的無機氮)，針對桑溝灣的養(yǎng)殖狀況，主要包括以下3個過程:海帶對氮的吸收，貝類對氮的排泄以及沉積物中氮的溶出量。

養(yǎng)殖自身污染壓力也采用污染壓力指數(shù)法進行評價，污染壓力指數(shù)等于該污染物通過養(yǎng)殖生物的入海通量除以評價海區(qū)中該污染物總量。

式中，Qj為評價海區(qū)中j養(yǎng)殖生物每年的養(yǎng)殖產(chǎn)量，單位為kg;Eij為在一個養(yǎng)殖周期中j養(yǎng)殖生物向水體中輸出的i污染物的量，單位為kg/kg。具體計算過程如下:

①海帶吸收氮的量 已有研究資料表明，海帶的無機氮含量為1.2%［17］，桑溝灣2006年和2007年海帶平均產(chǎn)量干重為94454 t［18］，則桑溝灣每年海帶吸收N的量為1133 t;

②貝類對氮的排泄量 周毅等［19］對幾種雙殼類對氮、磷的排泄速率進行了研究，結(jié)果表明櫛孔扇貝和牡蠣對 NH4-N 的排泄速率差異不大，分別為 3.69 μmol h－1個體－1和 3.46 μmol h－1個體－1，本文取二者的平均值(3.58 μmol h－1個體－1)估算養(yǎng)殖貝類對NH4-N 的釋放量。另外，根據(jù)Yuan等［20］提供的扇貝產(chǎn)量與扇貝個體數(shù)的換算(93 800 t=1.59×109個體)，桑溝灣2006年和2007年貝類的個體數(shù)約為9.0×108個體。

貝類釋放 NH4-N 的量=9.0×108個體×3.58μmol h－1個體－1×24 h×365 d×14×10－12t=395t。

劉學海［7］在對桑溝灣養(yǎng)殖海域生態(tài)系統(tǒng)模型中，參考“海岸帶復合系統(tǒng)中的生態(tài)海水養(yǎng)殖研究”項目的調(diào)查結(jié)果及有關文獻［21-22］，對桑溝灣櫛孔扇貝和牡蠣的養(yǎng)殖密度均采用59個體/m2。桑溝灣2006和2007年貝類平均養(yǎng)殖面積為15.55 km2，按照59個體/m2這樣的養(yǎng)殖密度計算的話，桑溝灣貝類的個體數(shù)=15.55×106m2×59 個體/m2=9.2×108個體，與用 Yuan 等［20］換算公式結(jié)果一致。

③沉積物氮的溶出量:根據(jù)孫珊等［23］的估算研究，桑溝灣2006—2007年沉積物-水體界面NH4-N的交換通量為 0.78—4.13 mg m－2d－1。利用其平均值(2.78 mg m－2d－1)和桑溝灣的總面積計算，估算得到桑溝灣全年(365 d)由沉積物擴散進入上覆水的總無機氮的量約為135 t。

近20年來對桑溝灣沉積物營養(yǎng)鹽溶出動力學的研究結(jié)果列于表2。盡管桑溝灣的貝類養(yǎng)殖已有30a的歷史，但與其他淺海環(huán)境相比，自1999年至今沉積物-水界面的營養(yǎng)鹽通量一直處于較低或中等水平，甚至表現(xiàn)出下降的趨勢(表2)，當然這可能與測定方法和取樣站點的不同有一定的關系。

綜合上述①—③的計算結(jié)果，桑溝灣每年氮凈釋放量=395+135－1133＜0。因此，自身污染壓力指數(shù)Pi＜0，桑溝灣自身污染壓力指標得分為5分。

從養(yǎng)殖規(guī)模與品種比例長期變化趨勢看(圖2)，由于海帶一直是桑溝灣的主要養(yǎng)殖品種，對營養(yǎng)鹽吸收量較大，因此盡管貝類對氮的排泄量和沉積物中氮的溶出量是外源污染壓力的兩倍，桑溝灣養(yǎng)殖自身污染壓力仍然較小。

表2 桑溝灣沉積物-海水界面NH4-N擴散通量時間變化Table 2 Diffusion fluxes of NH4-N between sediment and seawater in Sanggou Bay

(3)養(yǎng)殖品種多樣性/互補性

采用專家評判法，得分為4.25分。

綜上所述，桑溝灣養(yǎng)殖壓力模塊的得分為3.0(具體見表3)，總的壓力相對較高。其中養(yǎng)殖自身污染壓力得分最高(5分)，養(yǎng)殖品種多樣性得分也較高(4.25分)，養(yǎng)殖品種較多，不同養(yǎng)殖品種之間在空間和季節(jié)上的搭配比較合理;養(yǎng)殖密度/規(guī)模的得分為2分，說明桑溝灣養(yǎng)殖密度/規(guī)模壓力較大，擴展空間有限。總的來說，桑溝灣養(yǎng)殖業(yè)自身污染程度不大，養(yǎng)殖品種搭配也比較合理，但較大的養(yǎng)殖規(guī)模還是給生態(tài)系統(tǒng)造成了較高壓力。

表3 桑溝灣生態(tài)系統(tǒng)壓力指標評價結(jié)果Table 3 Assessment results of ecosystem pressure in Sanggou Bay

2.2 生態(tài)環(huán)境現(xiàn)狀與變化趨勢

本文從水文環(huán)境、水體環(huán)境、底質(zhì)環(huán)境以及生物群落質(zhì)量4個方面，對桑溝灣的生態(tài)環(huán)境現(xiàn)狀及可能的變化趨勢進行評價，評價結(jié)果總結(jié)于表4。

2.2.1 水文環(huán)境現(xiàn)狀

劉學海［7］通過海上觀測和數(shù)值模擬研究了桑溝灣水動力學特征。觀測表明養(yǎng)殖活動明顯降低了所在水層的流速。模型結(jié)果顯示，不考慮養(yǎng)殖影響該灣的半交換周期為7 d，僅考慮貝類養(yǎng)殖為12 d，同時考慮貝類和海帶養(yǎng)殖半交換周期達16 d。灣內(nèi)西北、西南部水交換能力最差。

綜合考慮桑溝灣水交換狀況的得分為4分。

2.2.2 水體環(huán)境現(xiàn)狀與變化趨勢

桑溝灣2006—2007年水體環(huán)境指標的評價結(jié)果顯示(表4)，海灣水體處于程度較輕的富營養(yǎng)化狀態(tài)，這主要是由較高的DIN濃度所致，底層水體溶解氧含量較高，重金屬和石油烴污染不嚴重，水體懸浮物評價得分最低，濃度處于中等水平。

將各評價指標與自1983年以來可得的歷史數(shù)據(jù)進行了比較:桑溝灣底層水體溶解氧含量一直較高(＞5 mg/L)，在2000年之后表現(xiàn)出一定程度的升高;水體營養(yǎng)鹽含量和結(jié)構變化較明顯(圖3)，DIN濃度在過去的20a表現(xiàn)出明顯的上升趨勢，PO4-P濃度一直維持在較低水平，N/P比值顯著上升，自2003年以來N/P比值已遠高于Redfield比值，富營養(yǎng)化指數(shù)也表現(xiàn)出逐漸上升的趨勢，雖然目前仍處于較輕富營養(yǎng)化水平，但隨著DIN濃度的進一步升高，水體可能會出現(xiàn)富營養(yǎng)化程度加重的風險。

圖3 桑溝灣水體營養(yǎng)鹽濃度與結(jié)構及富營養(yǎng)化指數(shù)的長期變化趨勢Fig.3 Long-term variations of nutrient concentrations and ratios and EI index in seawater of Sanggou Bay

3.2.3 底質(zhì)環(huán)境現(xiàn)狀

由于桑溝灣沉積物歷史調(diào)查數(shù)據(jù)較少，本文僅對底質(zhì)環(huán)境的現(xiàn)狀進行了評價。桑溝灣沉積物質(zhì)量總體狀況良好(表4):重金屬生態(tài)風險污染指數(shù)很低;氧化還原電位值較高，與底層水體的溶解氧含量的情況是一致的;有機碳平均含量中等為1.10%［26］，處于養(yǎng)殖海域的正常水平，表明桑溝灣底質(zhì)沒有受到貝類養(yǎng)殖生物沉積作用的明顯影響。

3.2.4 生物群落質(zhì)量狀況

多樣性是生物群落的重要特征，是評價生態(tài)系統(tǒng)功能與穩(wěn)定性的重要指標。本文采用Shannon多樣性指數(shù)對桑溝灣浮游植物、浮游動物和底棲生物群落進行了評價。結(jié)果表明，2006—2007年桑溝灣浮游植物(0.36—3.58)與浮游動物(0.82—2.26)多樣性指數(shù)變化范圍較大，平均得分分別為3.1和2.7分，大型底棲生物多樣性水平相對較高(1.5—3.73)，平均得分為4.2分。另外，根據(jù)專家對桑溝灣重要生境的評判結(jié)果，由于桑溝灣的養(yǎng)殖占有了原海灣部分重要生境，大葉藻的分布區(qū)顯著縮小，對原有生境造成中等程度干擾。

關于海水養(yǎng)殖對生態(tài)環(huán)境的影響已有較多研究報道，但對于養(yǎng)殖對生物群落多樣性的研究則較少。少量的研究表明海水養(yǎng)殖會對生物多樣性產(chǎn)生一定的負面影響，例如貝類養(yǎng)殖會通過對浮游植物的選擇性攝食而影響浮游植物的數(shù)量豐度和群落結(jié)構［8，14，27-28］。從桑溝灣網(wǎng)采浮游植物多樣性指數(shù)的長期變化看(圖4)，2000年之后出現(xiàn)下降趨勢。

圖4 桑溝灣浮游植物物種多樣性長期變化Fig.4 Long-term variation of phytoplankton biodiversity(H')in Sanggou Bay

2.3 生態(tài)系統(tǒng)響應

生態(tài)系統(tǒng)對外源壓力和環(huán)境變化的響應是目前生態(tài)系統(tǒng)健康評價中較少涉及的問題，本文采用專家評判法，對海水養(yǎng)殖生態(tài)系統(tǒng)中常見的生態(tài)異常問題進行了半定量評價(表5)。桑溝灣的有害赤潮和外來生物入侵問題不嚴重，但是由于隨著養(yǎng)殖規(guī)模的不斷擴大，該灣近來出現(xiàn)了海帶在收獲季節(jié)開始腐爛，櫛孔扇貝養(yǎng)殖個體小型化、養(yǎng)殖周期延長、死亡率升高、產(chǎn)品質(zhì)量下降等問題，養(yǎng)殖病害和產(chǎn)品質(zhì)量下降是桑溝灣海水養(yǎng)殖生態(tài)系統(tǒng)的主要威脅之一。

表5 桑溝灣生態(tài)系統(tǒng)生態(tài)響應評價結(jié)果Table 5 Assessment result of the ecosystem response in Sanggou Bay

2.4 綜合評價結(jié)果

綜合3.1—3.3中對桑溝灣海水養(yǎng)殖生態(tài)系統(tǒng)的壓力、狀態(tài)和響應的評價，得到桑溝灣海水養(yǎng)殖生態(tài)系統(tǒng)健康綜合評價結(jié)果(表6)。桑溝灣生態(tài)系統(tǒng)壓力等級為中等，系統(tǒng)狀態(tài)等級為較好，系統(tǒng)響應等級為中等;總的生態(tài)系統(tǒng)健康指數(shù)為3.4，為較好健康水平。同時給出了有關管理建議，見表6。為了方便識別，采用不同的顏色對得分等級進行了標識。

3 結(jié)論

本文利用《海水養(yǎng)殖系統(tǒng)健康綜合評價:方法與模式》一文中建立的方法，對我國北方重要養(yǎng)殖型海灣桑溝灣的生態(tài)系統(tǒng)健康狀況進行了綜合評價，并結(jié)合可得的歷史數(shù)據(jù)對其時間變化趨勢進行了探討，主要結(jié)論如下:

(1)沽河輸入的營養(yǎng)鹽負荷和過高的養(yǎng)殖規(guī)模/密度是影響桑溝灣生態(tài)系統(tǒng)健康的主要壓力因素，盡管貝類排泄的無機氮與河流輸入量相當，但由于大規(guī)模海帶養(yǎng)殖對營養(yǎng)鹽的吸收，桑溝灣的養(yǎng)殖自身污染程度并不高;

(2)由于桑溝灣較好的水交換能力，以及貝藻混養(yǎng)的多元養(yǎng)殖模式，水體富營養(yǎng)化程度較輕，底層水體溶氧含量較高，大規(guī)模養(yǎng)殖活動對沉積物和底棲生物未造成明顯影響，但由于桑溝灣的養(yǎng)殖占有了原海灣部分重要生境，大葉藻的分布區(qū)顯著縮小，對原有生境造成中等程度干擾;

(3)養(yǎng)殖病害問題和養(yǎng)殖產(chǎn)品質(zhì)量下降，是桑溝灣生態(tài)系統(tǒng)對外源壓力和環(huán)境變化的主要響應方式，對生態(tài)系統(tǒng)健康具有中等程度威脅;

表6 桑溝灣海水養(yǎng)殖生態(tài)系統(tǒng)健康綜合評價表Table 6 Integrated assessment of Sanggou Bay marine aquaculture ecosystem health

(4)從評價指標的歷史變化趨勢看，桑溝灣的養(yǎng)殖規(guī)模在1990年后有大幅增加，近10a來養(yǎng)殖品種和規(guī)模變化不大，海帶養(yǎng)殖始終占有很大比例，能夠有效減輕貝類養(yǎng)殖的自身污染，水體無機氮濃度和氮磷比值顯著增加，具有潛在富營養(yǎng)化風險。

(5)桑溝灣生態(tài)系統(tǒng)壓力等級為中等，系統(tǒng)狀態(tài)等級為較好，系統(tǒng)響應等級為中等;總的生態(tài)系統(tǒng)健康指數(shù)為3.4，為較好健康水平。

［1］ Fang J G，Sun H L，Kang S H，Sun Y，Zhou S L，Song Y L，Cui Y，Zhao J，Yang Q F，Li F，Zhang A J，Wang X Z，Tang T Y.Mariculture status and optimising measurements for the culture of scallopColalys farreriand kelpLaminaria japonicain Sanggou Bay.Marine Fisheries Research，1996，17(2):95-102.

［2］ Zhang J H，Jiang Z J，Wang W，Zou J，Xue S Y，F(xiàn)ang J G，Lian Y，Zhang X J，Liu X J，Zhou Y M.Seasonal distribution and variation of nutrients and nutrients limitation in Sanggou Bay.Progress in Fishery Sciences，2010，31(4):16-25.

［3］ Sun P X，Zhang Z H，Hao L H，Wang B，Wang Z L，Liu P，Lian Y，Chang Z Y，Xie L P.Analysis of nutrient distributions and potential eutrophication in seawater of the Sanggou Bay.Advances in Marine Science，2007，25(4):436-445.

［4］ Song Y L，Cui Y，Sun Y，F(xiàn)ang J G，Sun H L，Yang Q F，Kuang S H.Study on nutrient state and influencing factors in Sanggou Bay.Marine Fisheries Research，1996，17(2):41-52.

［5］ Song H J，Li R X，Wang Z L，Zhang X L，Liu P.Interannual variations in phytoplankton diversity in the Sanggou Bay.Advances in Marine Science，2007，25(3):332-339.

［6］ Shi J.Numerical Study on the Influences of Physical Processes on the Aquaculture Carrying Capacity in a Semi-enclosed Bay［D］.Qingdao:Ocean University of China，2009.

［7］ Liu X H.Ecosystem Dynamic Models for South Yellow Sea and An Aquaculture Bay［D］.Qingdao:Ocean University of China，2008.

［8］ Ji R B，Mao X H，Zhu M Y.Impacts of coastal shellfish aquaculture on bay ecosystem.Journal of Oceanography of Huanghai and Bohai Seas，1998，16(1):21-27.

［9］ Zhang L H，Zhang X L，Li R X，Wang Z L，Li Y，Wang L C，Lian Y，Liu Y.Impact of scallop culture on dinoflagellate abundance in the Sanggou Bay.Advances in Marine Science，2005，23(3):342-346.

［10］ Zhang W X.Marine Aquaculture in Sanggou Bay.Beijing:Ocean Press，1992.

［11］ Wu J X.Development and Applications of Dynamic Model for Dissolved Nitrogen and Phosphate Budgets in Sungo Bay［D］.Qingdao:Ocean University of China，2005.

［12］ Zhang Z H，Lü J B，Ye S F，Zhu M Y.Values of marine ecosystem services in Sanggou Bay.Chinese Journal of Applied Ecology，2007，18(11):2540-2547.

［13］ Liu H，F(xiàn)ang J G，Dong S L，Wang L C，Lian Y.Annual variation of major nutrients and limiting factors in Laizhou Bay and Sanggou Bay.Journal of Fishery Sciences of China，2003，10(3):227-234.

［14］ Zhang J H.Effect of Filter Feeding Shellfish Mariculture on the Ecosystem and the Evaluation of Ecology Carrying Capacity［D］.Qingdao:Institute of Oceanology Chinese Academy of Sciences，2008.

［15］ Marine Environment Quality Bulletin of Shandong Province，2008.Ji'nan:Ocean and Fisheries Bureau of Shandong Province，2009.

［16］ Marine Environment Quality Bulletin of Shandong Province，2009.Ji'nan:Ocean and Fisheries Bureau of Shandong Province，2010.

［17］ Fang J G，Sun H L，Kang S H，Sun Y，Zhou S H，Song Y L，Cui Y，Zhao J，Yang Q F，Li F，Wang X Z，Tang T Y.Assessing the carrying capacity of Sanggou Bay for culture of kelp Laminaria japonica.Marine Fisheries Research，1996，17(2):7-17.

［18］ Zhang J H，Hansen P K，F(xiàn)ang J G，Wang W，Jiang Z J.Assessment of the local environmental impact of intensive marine shellfish and seaweed farming—Application of the MOM system in the Sungo Bay，China.Aquaculture，2009，287(3/4):304-310.

［19］ Zhou Y，Yang H S，He Y C，Zhang F S.Nitrogen and phosphorus excretion and its ecological effect by several bivalves and fouling animals.Oceanologia Et Limnologia Sinica，2002，33(4):424-431.

［20］ Yuan X T，Zhang M J，Liang Y B，Liu D，Guan D M.Self-pollutant loading from a suspension aquaculture system of Japanese scallop(Patinopecten yessoensis)in the Changhai sea area，Northern Yellow Sea of China.Aquaculture，2010，304(1/4):79-87.

［21］ Nunes J P，F(xiàn)erreira J G，Gazeau F，Lencart-Silva J，Zhang X L，Zhu M Y，F(xiàn)ang J G.A model for sustainable management of shellfish polyculture in coastal bays.Aquaculture，2003，219(1/4):257-277.

［22］ Zhu M Y，Zhang X L，Tang T Y，F(xiàn)erreira J G，F(xiàn)ang J G，Wang X Z.Application of ecological model in studing the sustainable development of coastal shellfish culture.Advances in Marine Science，2002，20(4):34-42.

［23］ Sun S，Liu S M，Ren J L，Zhang J H，Jiang Z J.Distribution features of nutrients and flux across the sediment-water interface in the Sanggou Bay.Acta Oceanologica Sinica，2010，32(6):108-117.

［24］ Zhang X L，Zhu M Y，Tang T Y，Wang Y，Liu G，Martin J L M.Fluxes of nutrients at sediment-water in Sanggou Bayand Jiaozhou Bay in summer.Marine Environmental Science，2004，23(1):1-4.

［25］ Cai L S，F(xiàn)ang J G，Dong S L.Preliminary studies on nitrogen and phosphorus fluxes between seawater and sediment in Sungo Bay.Marine Fisheries Research，2004，25(4):57-64.

［26］ Wu J X，Sun Y，Zhang Q Q，Wang X L.Research on the exchange rates of TOC，TN，TP at the sediment-water interface in aquaculture water areas of Sungo Bay.Marine Fisheries Research，2005，26(2):62-67.

［27］ Wu L J，Yang L.Advances on impact of mariculture on biodiversity.Marine Environmental Science，2008，27(4):397-400.

［28］ Song H J，Li R X，Wang Z L，Zhang X L，Liu P.Interannual variations in phytoplankton diversity in the Sanggou Bay.Advances in Marine Science，2007，25(3):332-339.

參考文獻:

［1］ 方建光，孫慧玲，匡世煥，孫耀，周詩賚，宋云利，崔毅，趙俊，楊琴芳，李鋒，張愛君，王興章，湯庭耀.桑溝灣海水養(yǎng)殖現(xiàn)狀評估及優(yōu)化措施.海洋水產(chǎn)研究，1996，17(2):95-102.

［2］ 張繼紅，蔣增杰，王巍，鄒建，薛素艷，方建光，連巖，張新軍，劉新杰，周穎敏.桑溝灣營養(yǎng)鹽時空分布及營養(yǎng)鹽限制分析.漁業(yè)科學進展，2010，31(4):16-25.

［3］ 孫丕喜，張朝暉，郝林華，王波，王宗靈，劉萍，連巖，常忠岳，謝琳萍.桑溝灣海水中營養(yǎng)鹽分布及潛在性富營養(yǎng)化分析.海洋科學進展，2007，25(4):436-445.

［4］ 宋云利，崔毅，孫耀，方建光，孫慧玲，楊琴芳，匡世煥.桑溝灣養(yǎng)殖海域營養(yǎng)狀況及其影響因素分析.海洋水產(chǎn)研究，1996，17(2):41-52.

［5］ 宋洪軍，李瑞香，王宗靈，張學雷，劉萍.桑溝灣浮游植物多樣性年際變化.海洋科學進展，2007，25(3):332-339.

［6］ 史潔.物理過程對半封閉海灣養(yǎng)殖容量影響的數(shù)值研究［D］.青島:中國海洋大學，2009.

［7］ 劉學海.南黃海及養(yǎng)殖功能海域生態(tài)動力學模型研究［D］.青島:中國海洋大學，2009.

［8］ 季如寶，毛興華，朱明遠.貝類養(yǎng)殖對海灣生態(tài)系統(tǒng)的影響.黃渤海海洋，1998，16(1):21-27.

［9］ 張莉紅，張學雷，李瑞香，王宗靈，李艷，王立超，連巖，劉瑤.桑溝灣扇貝養(yǎng)殖對甲藻數(shù)量的影響.海洋科學進展，2005，23(3):342-346.

［10］ 張為先.桑溝灣增養(yǎng)殖.北京:海洋出版社，1992.

［11］ 武晉宣.桑溝灣養(yǎng)殖海域氮、磷收支及環(huán)境容量模型［D］.青島:中國海洋大學，2005.

［12］ 張朝暉，呂吉斌，葉屬峰，朱明遠.桑溝灣海洋生態(tài)系統(tǒng)的服務價值.應用生態(tài)學報，2007，18(11):2540-2547.

［13］ 劉慧，方建光，董雙林，王立超，連巖.萊州灣和桑溝灣養(yǎng)殖海區(qū)主要營養(yǎng)鹽的周年變動及限制因子.中國水產(chǎn)科學，2003，10(3):227-234.

［14］ 張繼紅.濾食性貝類養(yǎng)殖活動對海域生態(tài)系統(tǒng)的影響及生態(tài)容量評估［D］.青島:中科院海洋研究所，2008.

［15］ 2008年山東省海洋環(huán)境質(zhì)量公報.濟南:山東省海洋與漁業(yè)廳，2009.

［16］ 2009年山東省海洋環(huán)境質(zhì)量公報.濟南:山東省海洋與漁業(yè)廳，2010.

［17］ 方建光，孫慧玲，匡世煥，孫耀，周詩賚，宋云利，崔毅，趙俊，楊琴芳，李鋒，王興章，湯庭耀.桑溝灣海帶養(yǎng)殖容量的研究.海洋水產(chǎn)研究，1996，17(2):7-17.

［19］ 周毅，楊紅生，何義朝，張福綏.四十里灣幾種雙殼貝類及污損動物的氮、磷排泄及其生態(tài)效應.海洋與湖沼，2002，33(4):424-431.

［22］ 朱明遠，張學雷，湯庭耀，F(xiàn)erreira J G，方建光，王興章.應用生態(tài)模型研究近海貝類養(yǎng)殖的可持續(xù)發(fā)展.海洋科學進展，2002，20(4):34-42.

［23］ 孫珊，劉素美，任景玲，張繼紅，蔣增杰.桑溝灣養(yǎng)殖海域營養(yǎng)鹽和沉積物-水界面擴散通量研究.海洋學報，2010，32(6):108-117.

［24］ 張學雷，朱明遠，湯庭耀，王英，劉崗，Martin J L M.桑溝灣和膠州灣夏季的沉積物-水界面營養(yǎng)鹽通量研究.海洋環(huán)境科學，2004，23(1):1-4.

［25］ 蔡立勝，方建光，董雙林.桑溝灣養(yǎng)殖海區(qū)沉積物-海水界面氮、磷營養(yǎng)鹽的通量.海洋水產(chǎn)研究，2004，25(4):57-64.

［26］ 武晉宣，孫耀，張前前，王曉麗.桑溝灣養(yǎng)殖水域沉積物中營養(yǎng)要素(TOC、TN和TP)溶出動力學特性.海洋水產(chǎn)研究，2005，26(2):62-67.

［27］ 吳隆杰，楊林.海水養(yǎng)殖對生物多樣性的影響研究進展.海洋環(huán)境科學，2008，27(4):397-400.

［28］ 宋洪軍，李瑞香，王宗靈，張學雷，劉萍.桑溝灣浮游植物多樣性年際變化.海洋科學進展，2007，25(3):332-339.