全為民,張?jiān)茙X,齊遵利,許 敏,范瑞良,王桃妮,李楠楠,孫兆躍,周海生,李 春,張秀文

1 中國(guó)水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所,農(nóng)業(yè)農(nóng)村部東海漁業(yè)資源開(kāi)放利用重點(diǎn)實(shí)驗(yàn)室,上海 200090 2 唐山海洋牧場(chǎng)實(shí)業(yè)有限公司,唐山 063611 3 河北農(nóng)業(yè)大學(xué)海洋學(xué)院, 秦皇島 066003

牡蠣礁是由大量牡蠣固著生長(zhǎng)于硬底物表面所形成的一種生物礁系統(tǒng),被喻為溫帶地區(qū)的“珊瑚礁”[1—3]。除為人類提供大量鮮活牡蠣以供食用外,牡蠣礁具有凈化水體、維持生物多樣性、提供魚類生境、固碳和防止岸線侵蝕等生態(tài)功能[3—5]。牡蠣礁曾廣泛分布于全球亞熱帶和溫帶海區(qū)的多數(shù)海灣、河口以及部分海岸線的潮間帶和潮下帶。但過(guò)度捕撈、環(huán)境污染、生境破壞和病害侵染使全球野生牡蠣資源量急劇下降,牡蠣礁生境遭受破壞或喪失,過(guò)去100多年來(lái)全球牡蠣礁面積約下降了85%,從而改變了近岸生態(tài)系統(tǒng)的結(jié)構(gòu)與功能[4—5]。

牡蠣礁作為重要的海岸帶棲息生境之一,曾在我國(guó)溫帶和亞熱帶海區(qū)的潮間帶和淺水潮下帶有著廣泛分布[6—12],如20世紀(jì)80年代我國(guó)地質(zhì)學(xué)者曾在我國(guó)南北沿海的許多河口和海灣內(nèi)記錄到牡蠣礁分布,然而這些研究更多關(guān)注牡蠣礁的地理分布、海岸地貌發(fā)育和古環(huán)境,而很少研究牡蠣礁的生態(tài)現(xiàn)狀和生態(tài)系統(tǒng)服務(wù)價(jià)值。2000年以后,國(guó)內(nèi)才開(kāi)始關(guān)注牡蠣礁的生態(tài)功能[1],并對(duì)天津大神堂[13]和江蘇小廟洪[9—10]等少數(shù)現(xiàn)存活體牡蠣礁的生態(tài)現(xiàn)狀進(jìn)行零星的調(diào)查研究,結(jié)果顯示這些牡蠣礁退化嚴(yán)重,大部分自然牡蠣礁已經(jīng)消失[10,13—15]。截止目前,國(guó)內(nèi)有關(guān)牡蠣礁的生態(tài)學(xué)研究仍然有限,無(wú)法全面掌握目前我國(guó)自然牡蠣礁的地理分布及生態(tài)現(xiàn)狀,制約著我國(guó)海岸帶的保護(hù)與修復(fù)。

渤海灣曾是我國(guó)牡蠣礁的主要分布區(qū),但由于海岸帶歷史變遷等原因,大量牡蠣礁被埋入地下成為古牡蠣礁,而該區(qū)現(xiàn)存活體牡蠣礁的報(bào)道較少[6,13]?！稙粗葜尽肪砭拧ぁ斗庥?下)·古跡漕運(yùn)附》載:“知府孫維城復(fù)議海運(yùn)故道直捷,而所經(jīng)蠶沙綠洋各口蠣房山石臼坨……”,所描述的“蠣房山”即為牡蠣礁；據(jù)歷史資料和當(dāng)?shù)貪O民描述,該海域自然牡蠣礁分布于溯河河口以東、小青河河口以西的低潮線以下海域,后來(lái)牡蠣礁退化或消失。近年來(lái),在原有牡蠣礁分布海域又能捕撈到大量活體牡蠣,極大的野生牡蠣捕撈量引起社會(huì)的廣泛關(guān)注,拯救并保護(hù)這一珍貴的自然資源及其生境勢(shì)在必行。在實(shí)施自然牡蠣礁的保護(hù)與修復(fù)措施前,亟待解決以下3個(gè)科學(xué)問(wèn)題:(1)該自然牡蠣礁的地理分布區(qū)域及分布范圍？(2)該自然牡蠣礁中牡蠣種群及礁體動(dòng)物群落的生態(tài)現(xiàn)狀如何？(3)該自然牡蠣礁提供多大的生態(tài)系統(tǒng)服務(wù)價(jià)值？為回答以上科學(xué)問(wèn)題,本研究于2019年3月對(duì)該自然牡蠣礁的空間分布、生態(tài)環(huán)境、牡蠣生物學(xué)和礁體動(dòng)物群落開(kāi)展了調(diào)查,并評(píng)估了該自然牡蠣礁的生態(tài)系統(tǒng)服務(wù)價(jià)值,研究結(jié)果可為該自然牡蠣礁的保護(hù)與修復(fù)提供基礎(chǔ)信息。

1 材料與方法

1.1 研究地點(diǎn)

研究地點(diǎn)位于河北省唐山市曹妃甸-樂(lè)亭海域(圖1),是介于曹妃甸工業(yè)區(qū)與京唐港之間的半封閉海灣,具有離岸沙島、深槽和淺灘等多種地貌單元。海區(qū)潮汐型態(tài)系數(shù)為0.77,屬于不規(guī)則半日潮,平均高潮位為0.81 m,平均低潮位為-0.73 m,平均潮差為1.54 m[16]。根據(jù)歷史文獻(xiàn)、現(xiàn)場(chǎng)調(diào)查和牡蠣捕撈統(tǒng)計(jì)資料,初步確定溯河口(SRE)、溯河(SR)和撈魚尖(LYJ)3個(gè)自然牡蠣礁分布區(qū)。共設(shè)置12個(gè)生態(tài)調(diào)查站位,其中SR牡蠣礁3個(gè)站位、SRE牡蠣礁5個(gè)站位和LYJ牡蠣礁4個(gè)站位(圖1)。

1.2 調(diào)查與分析方法

2019年3月24—30日?qǐng)?zhí)行野外生態(tài)調(diào)查,調(diào)查內(nèi)容包括生態(tài)環(huán)境(水溫、鹽度、總懸浮顆粒物(TSS)、葉綠素a(Chla)和浮游植物)、牡蠣礁分布、牡蠣生物學(xué)(種類、密度、生物量、殼高(SH)、肥滿度(CI)和遺傳多樣性)和礁體動(dòng)物群落(種類、密度和生物量)。

圖1 研究地點(diǎn)和自然牡蠣礁分布范圍Fig.1 The study site and spatial distribution of the natural oyster reefs on the Caofeidian-Leting coast 海草床空間分布信息引自文獻(xiàn)[17]

在高潮位時(shí),采用YSI水質(zhì)分析儀現(xiàn)場(chǎng)測(cè)定表層海水水溫和鹽度；用采水器采集1000 mL表層海水,等分成兩份,現(xiàn)場(chǎng)用裝有0.45 μm微孔濾膜的便攜式真空泵抽濾,濾膜冰凍保存后帶回實(shí)驗(yàn)室分析測(cè)定,海水Chla和TSS濃度分別采用分光光度計(jì)法和重量法測(cè)定,另每站取1000 mL表層海水用魯哥式溶液現(xiàn)場(chǎng)固定后帶回室內(nèi)鏡檢計(jì)數(shù)浮游植物。

采用社會(huì)調(diào)查、潛水調(diào)查和水下攝像等方法調(diào)查自然牡蠣礁的地理分布,具體方法為:首先和當(dāng)?shù)啬迪牪稉茝臉I(yè)人員進(jìn)行訪談,初步獲取牡蠣重點(diǎn)采捕區(qū)的位置和范圍,再提取每個(gè)區(qū)域牡蠣采捕漁船的GPS導(dǎo)航航跡數(shù)據(jù),最后采用潛水調(diào)查和水下攝像進(jìn)行驗(yàn)證,綜合上述信息確定牡蠣礁地理位置和分布范圍。在SRE和SR,由潛水員下潛至牡蠣礁表面進(jìn)行采樣,每個(gè)站位采集1個(gè)0.5 m×0.5 m樣方內(nèi)(深度10cm)所有的活體牡蠣、牡蠣殼和生物樣品,裝于網(wǎng)目0.5 mm的聚乙烯網(wǎng)袋中帶回岸上分檢。在LYJ,低潮位時(shí)水深僅約0.4m,直接手工采樣。

回岸后分撿每個(gè)樣品中的活牡蠣和礁體動(dòng)物。記錄活牡蠣的個(gè)體數(shù)、測(cè)定其總重(精確至0.1 g)；從每個(gè)調(diào)查區(qū)樣品中隨機(jī)選取不少于50個(gè)牡蠣,用游標(biāo)卡尺測(cè)量殼高(SH,精確至1mm),帶部分牡蠣樣品回室內(nèi)分析種類和肥滿度。礁體動(dòng)物(不包括牡蠣)樣品用75%乙醇固定保存后,帶回室內(nèi)鑒定至最低的分類單元,并計(jì)數(shù)和稱重。

基于外殼形態(tài),選取96個(gè)牡蠣樣品進(jìn)行線粒體16S rDNA 基因序列分析,據(jù)此分析牡蠣的種類和遺傳多樣性。從每個(gè)調(diào)查區(qū)樣品中隨機(jī)選取不少于20個(gè)牡蠣,剝?nèi)∧迪牭能涹w部分,測(cè)定濕肉重(WTW,g),然后將軟體部分放入烘箱中(60℃)烘干至恒重(約24 h),并測(cè)定其干肉重(FTW,g),依據(jù)下列公式計(jì)算牡蠣的肥滿度(Condition index,CI):

網(wǎng)絡(luò)語(yǔ)言的出現(xiàn)使得人們學(xué)習(xí)漢語(yǔ)言有了束縛，主要是使人們?cè)谖淖钟洃浬洗嬖诹艘欢ǖ娜毕?。在網(wǎng)絡(luò)不斷發(fā)展的背景下，人們經(jīng)常為了便捷運(yùn)用電子設(shè)備進(jìn)行書寫，輸入法的多樣化使得人們?cè)谳斎肷辖?jīng)常出現(xiàn)錯(cuò)別字，從而在實(shí)際書寫上就會(huì)受到影響。同時(shí)，人們的交流環(huán)境也會(huì)受到網(wǎng)絡(luò)語(yǔ)言的污染，其中暗含了很多污穢的內(nèi)容，難以讓人和諧交流。由于不良信息的出現(xiàn)使人們形成了不良思想。

1.3 數(shù)據(jù)統(tǒng)計(jì)

將測(cè)得的16S rDNA序列同GenBank中登錄的牡蠣序列進(jìn)行系統(tǒng)發(fā)育分析,采用Phylip 3.5軟件的鄰接法構(gòu)建系統(tǒng)樹(shù),利用Donaxtrunculus作為外群,置信號(hào)由1000次重復(fù)檢驗(yàn)。

牡蠣或礁體動(dòng)物的密度(個(gè)/m2)統(tǒng)計(jì)為每1m2牡蠣礁中的個(gè)體數(shù),生物量(g/m2或kg/m2)為每1 m2牡蠣礁中的總鮮重。運(yùn)用一維方差分析(1-way ANOVA)檢驗(yàn)?zāi)迪犞笜?biāo)(密度、生物量和SH)或礁體動(dòng)物群落指標(biāo)(密度和生物量)在3個(gè)分布區(qū)(SR,SRE和LYJ)間的差異(Tukey后檢驗(yàn))。如數(shù)據(jù)沒(méi)有通過(guò)方差同質(zhì)性和方差齊性檢驗(yàn),則采用一維排序方差分析(1-way ANOVA on rank,Ducnn后檢驗(yàn))。另采用Pearson相關(guān)性檢驗(yàn)?zāi)迪犞笜?biāo)與礁體動(dòng)物群落指標(biāo)之間的相關(guān)關(guān)系；所有統(tǒng)計(jì)分析及制圖均在SigmaPlot 10.0軟件中完成。

1.4 生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估

根據(jù)北美自然牡蠣礁的研究結(jié)果[3],牡蠣礁的生態(tài)系統(tǒng)服務(wù)價(jià)值標(biāo)準(zhǔn)為10.6萬(wàn)美元 hm-2a-1。據(jù)此標(biāo)準(zhǔn),本研究評(píng)估了該自然牡蠣礁的生態(tài)系統(tǒng)服務(wù)價(jià)值。

2 結(jié)果與分析

2.1 環(huán)境指標(biāo)

調(diào)查期間海水溫度介于4.1—10.9 ℃,平均為8.2℃；海水鹽度介于30.8—34.2,平均為33.4；TSS濃度介于12.2—43.7 mg/L,平均為27.8 mg/L；Chla濃度介于0.84—3.45 mg/m3,平均為2.37 mg/m3。共采集到浮游植物2門21種,其中硅藻17種和甲藻4種,浮游植物豐度介于(0.1—23.94)×104個(gè)細(xì)胞/L,平均為3.29×104個(gè)細(xì)胞/L,SRE中浮游植物豐度高于SR和LYJ,主要優(yōu)勢(shì)種為旋鏈角毛藻Chaetoceroscurvisetus。

表1 3個(gè)牡蠣礁區(qū)中環(huán)境指標(biāo)比較

2.2 牡蠣礁分布

自然牡蠣礁主要分布于以下3個(gè)區(qū)域(圖1)。①SRE牡蠣礁:主要分布于老龍溝潮汐通道兩則0—5 m水深的淺水區(qū),有6個(gè)斑塊狀分布區(qū),單個(gè)牡蠣礁面積介于0.17—3.62 km2之間,總面積約10 km2,礁體發(fā)育較好,礁體高度約0.3—0.5 m,是牡蠣的主要采捕區(qū)。②SR牡蠣礁:分布于溯河河道及溯河河口段,面積約0.5 km2,尤以河床底部的牡蠣礁發(fā)育較好。③LYJ牡蠣礁:位于石臼坨和腰坨西側(cè),面積約5 km2,單個(gè)礁體斑塊面積介于0.2—1.0 m2之間,蓋度約5%,處于發(fā)育早期,牡蠣與大型海藻共生。

2.3 牡蠣指標(biāo)

經(jīng)PCR擴(kuò)增,獲得了96個(gè)牡蠣的16S rDNA基因片段擴(kuò)增產(chǎn)物(圖2),可比16S rDNA序列為450 bp左右,共檢測(cè)到20個(gè)單倍型,其中92個(gè)為長(zhǎng)牡蠣Crassostreagigas、3個(gè)為侏儒牡蠣Nanostreafluctigera、1個(gè)為巨蠣屬未知種Crassostreasp.(圖3)。長(zhǎng)牡蠣有17個(gè)單倍型、單倍型多樣性為0.641、核苷酸多樣性為0.00278；侏儒牡蠣有2個(gè)單倍型、單倍型多樣性為0.667、核苷酸多樣性為0.00147。

圖2 部分牡蠣樣本的線粒體16S rDNA擴(kuò)增電泳圖譜Fig.2 Amplified products of mitochondrial 16S rDNA gene of partial oyster specimens 1—22 為牡蠣樣本,M為DL2000分子標(biāo)記

圖3 基于96個(gè)牡蠣樣本16S rDNA數(shù)據(jù)構(gòu)建的系統(tǒng)發(fā)育樹(shù)Fig.3 Phylogenetic trees constructed from 16S rDNA sequence of 96 oyster specimens CFDLT: 曹妃甸-樂(lè)亭海域的牡蠣樣本Oyster specimens from Caofeidian-Leting coast

SRE中牡蠣密度介于2764—3912 個(gè)/m2,平均為(3251±223)個(gè)/m2；SR中牡蠣密度介于1008—1200 個(gè)/m2,平均為(1135±63)個(gè)/m2；LYJ中牡蠣密度介于104—280 個(gè)/m2,平均為(173±41)個(gè)/m2；牡蠣密度的大小排序?yàn)?SRE>SR>LYJ (圖4,Tukey后檢驗(yàn),P<0.05)；牡蠣生物量的大小排序?yàn)?SR >SRE>LYJ (圖4,Tukey后檢驗(yàn),P<0.05)。

圖5顯示了3個(gè)自然牡蠣礁分布區(qū)中牡蠣的殼高—頻率分布。SR中牡蠣SH(中值SH:58 mm)顯著大于SRE(中值SH:40 mm)和LYJ(中值SH:33 mm)(Dunn后檢驗(yàn),P<0.05),而SRE和LYJ間沒(méi)有差異(Dunn后檢驗(yàn),P>0.05)。

圖4 3個(gè)牡蠣礁區(qū)間牡蠣密度和生物量的比較Fig.4 Comparisons on oyster densities and biomass among three oyster reef zones 誤差棒為±1 標(biāo)準(zhǔn)誤,不同小寫字母表示顯著性差異(P<0.05)

圖5 3個(gè)牡蠣礁區(qū)間牡蠣殼高—頻率分布的比較Fig.5 Comparisons on oyster SH-frequency distributions among three oyster reef zones

2.4 礁體動(dòng)物群落指標(biāo)

通過(guò)定性和定量采樣,在牡蠣礁內(nèi)記錄到49種礁體動(dòng)物,其中軟體動(dòng)物16種、節(jié)肢動(dòng)物16種、環(huán)節(jié)動(dòng)物8種、棘皮動(dòng)物5種、腔腸動(dòng)物2種、星形動(dòng)物和脊索動(dòng)物各1種(表2)。

SRE中礁體動(dòng)物密度介于404—552 個(gè)/m2,平均為(493±26)個(gè)/m2；SR中礁體動(dòng)物密度介于164—320 個(gè)/m2,平均為(228±47)個(gè)/m2；LYJ中礁體動(dòng)物密度介于28—44 個(gè)/m2,平均為(35±5)個(gè)/m2。平均總密度的大小排序?yàn)?SRE>SR>LYJ (圖6,Tukey后檢驗(yàn),P<0.05)；平均總生物量的排序?yàn)?SRE顯著高于LYJ(圖6,Dunn后檢驗(yàn),P<0.05),而SR介于SRE和LYJ之間(圖6,Dunn后檢驗(yàn),P>0.05)。

2.5 牡蠣指標(biāo)與礁體動(dòng)物群落指標(biāo)之間的相關(guān)性

牡蠣密度與礁體動(dòng)物密度和生物量呈現(xiàn)顯著的正相關(guān)性(表3,圖7,P<0.05),而牡蠣生物量與礁體動(dòng)物密度和生物量之間的相關(guān)性并不顯著(表3,P>0.05)。

2.6 生態(tài)系統(tǒng)服務(wù)價(jià)值

由于SR牡蠣礁面積較小和LYJ牡蠣礁中牡蠣密度較低,僅評(píng)估SRE牡蠣礁的生態(tài)系統(tǒng)服務(wù)價(jià)值,如該自然牡蠣礁未受破壞,其生態(tài)系統(tǒng)服務(wù)價(jià)值將達(dá)到1.06×108美元/a,按當(dāng)前匯率6.51計(jì)折合成人民幣約6.90×108元/年。

3 討論

3.1 自然牡蠣礁的地理分布

渤海灣沿岸曾是我國(guó)牡蠣礁的集中分布區(qū),如在天津沿岸發(fā)現(xiàn)了許多成片的古牡蠣礁[18]；除天津大神堂有活體牡蠣礁分布外,在渤海灣其它海域未見(jiàn)活體牡蠣礁的報(bào)道。從史書記載,曹妃甸-樂(lè)亭海域歷史上分布有大量牡蠣礁,但沒(méi)有這片牡蠣礁地理分布和生態(tài)現(xiàn)狀的基礎(chǔ)信息。受人力和資金限制,本研究并未采用聲吶勘測(cè)手段精確調(diào)查自然牡蠣礁的地理分布,但通過(guò)社會(huì)統(tǒng)計(jì)和現(xiàn)場(chǎng)調(diào)查等多種方法初步描述了該海域自然牡蠣礁的地理分布,估算牡蠣礁總面積高達(dá)15 km2,表明該牡蠣礁是截至目前為止我國(guó)面積最大的自然活體牡蠣礁。這是繼曹妃甸海草場(chǎng)[17]后在該海域的又一重大發(fā)現(xiàn),進(jìn)一步補(bǔ)充了我國(guó)牡蠣礁地理分布及生態(tài)現(xiàn)狀的基礎(chǔ)信息,研究成果對(duì)于曹妃甸-樂(lè)亭海域海岸帶管理和保護(hù)具有重要意義。

本研究發(fā)現(xiàn)3個(gè)自然分布區(qū)中以溯河口牡蠣礁發(fā)育最好、牡蠣密度最高和生物多樣性最為豐富,原因可能是該牡蠣礁分布于老龍溝潮汐通道兩側(cè)的淺水區(qū),水深約1—5 m,老龍溝的漲潮流攜帶而來(lái)的豐富浮游生物為牡蠣生長(zhǎng)和牡蠣礁發(fā)育提供了充足的餌料[19]。與歷史分布區(qū)相比,該海域現(xiàn)代牡蠣礁分布區(qū)整體向南發(fā)育,原因可能是西南側(cè)實(shí)施的曹妃甸圍填海工程減弱了溯河口海域的水動(dòng)力,從而為牡蠣附著生長(zhǎng)和牡蠣礁發(fā)育提供更為穩(wěn)定的水動(dòng)力環(huán)境,減弱的水動(dòng)力環(huán)境和增加的水體滯留時(shí)間更有利于牡蠣幼體的附著與生長(zhǎng)[20—21],加快了現(xiàn)代牡蠣礁向外側(cè)發(fā)育。

3.2 牡蠣和礁體動(dòng)物多樣性

本研究采用線粒體16S rDNA 基因序列變異分析發(fā)現(xiàn),自然牡蠣礁內(nèi)分布有長(zhǎng)牡蠣和侏儒牡蠣,其中以長(zhǎng)牡蠣在數(shù)量上占絕對(duì)優(yōu)勢(shì),是造礁的牡蠣種類。長(zhǎng)牡蠣原產(chǎn)于日本、韓國(guó)和中國(guó)江蘇連云港以北沿海,為世界性養(yǎng)殖種；侏儒牡蠣分布于印度一西太平洋淺水區(qū),我國(guó)最早報(bào)道于東海北部水深47 m海區(qū)[22],有關(guān)該種的報(bào)道非常少,本研究首次在曹妃甸-樂(lè)亭海域記錄到這一物種,為闡明該種的生物地理格局提供了基礎(chǔ)信息。

牡蠣密度和大小是表征牡蠣礁生態(tài)現(xiàn)狀的通用指標(biāo)[23]。本研究發(fā)現(xiàn)該海域是以長(zhǎng)牡蠣為造礁種的活體牡蠣礁,3個(gè)自然牡蠣礁區(qū)間牡蠣密度差異很大,尤以溯河口海域中牡蠣密度最高,平均密度達(dá)到(3251 ±223)個(gè)牡蠣/m2,明顯高于國(guó)內(nèi)外自然或人工修復(fù)牡蠣礁中的牡蠣密度[24—28]；3個(gè)牡蠣礁區(qū)中牡蠣殼高均呈現(xiàn)單峰的正態(tài)分布,分布有一定數(shù)量的牡蠣成體和稚貝,表明該自然牡蠣礁擁有1個(gè)健康的、可持續(xù)的牡蠣種群；如果加以有效保護(hù),該牡蠣礁將持續(xù)發(fā)揮顯著的生態(tài)系統(tǒng)服務(wù)功能。

表2 自然牡蠣礁定居性大型底棲動(dòng)物名錄

表3 牡蠣指標(biāo)和礁體動(dòng)物群落指標(biāo)之間的相關(guān)系數(shù)(P)

圖7 3個(gè)牡蠣礁區(qū)中牡蠣密度和礁體動(dòng)物群落總密度和總生物量之間的線性回歸關(guān)系Fig.7 Linear regressions between oyster densities and benthic macroinvertebrate communities metrics (densities and biomass) across three oyster reef zones

定居性礁體動(dòng)物群落是評(píng)估牡蠣礁生態(tài)現(xiàn)狀、生物多樣性和生境價(jià)值的重要指標(biāo)之一[23]。本研究報(bào)道的礁體動(dòng)物密度明顯低于江蘇海門蠣蚜山牡蠣礁[9,28]、長(zhǎng)江口牡蠣礁[27]和美國(guó)東海岸牡蠣礁[24—26],原因可能是調(diào)查期間該海域平均水溫低于10℃,許多底棲動(dòng)物仍未經(jīng)歷春季的資源補(bǔ)充；但牡蠣礁內(nèi)動(dòng)物密度仍高于曹妃甸近岸海域底棲動(dòng)物密度[29]。另外,本研究也發(fā)現(xiàn)牡蠣密度與礁體動(dòng)物群落總密度和總生物量之間呈現(xiàn)顯著的正相關(guān)關(guān)系,原因是牡蠣密度越高,其聚集生長(zhǎng)形成的礁體結(jié)構(gòu)越復(fù)雜,而生境結(jié)構(gòu)的復(fù)雜性與其維持的礁體動(dòng)物豐度之間有顯著的正相關(guān)性[30]。

3.3 牡蠣礁保護(hù)與生態(tài)修復(fù)策略

牡蠣礁是生態(tài)系統(tǒng)服務(wù)價(jià)值最高的海洋生境之一?；诒泵雷匀荒迪牻傅脑u(píng)估結(jié)果表明[3],牡蠣礁的生態(tài)系統(tǒng)服務(wù)價(jià)值約為10.6萬(wàn)美元 hm-2a-1,高于紅樹(shù)林(8.2萬(wàn)美元 hm-2a-1)、海草場(chǎng)(3.1萬(wàn)美元 hm-2a-1)和永久濕地(2.1萬(wàn)美元 hm-2a-1)。如對(duì)該自然牡蠣礁加以保護(hù)并免遭破壞,其生態(tài)系統(tǒng)服務(wù)價(jià)值將高達(dá)6.90×108元/年；以鮮活牡蠣平均單價(jià)10元/kg計(jì)算,則該牡蠣礁產(chǎn)生的生態(tài)系統(tǒng)服務(wù)價(jià)值相當(dāng)于每年生產(chǎn)0.69×108kg牡蠣的經(jīng)濟(jì)產(chǎn)值,按牡蠣礁面積10 km2計(jì)算,則每1 m2牡蠣礁的生態(tài)系統(tǒng)服務(wù)價(jià)值相當(dāng)于6.9 kg可捕規(guī)格牡蠣的經(jīng)濟(jì)產(chǎn)值(69元)。

自2018年春季開(kāi)始,該自然牡蠣礁持續(xù)面臨較高的捕撈壓力,主要捕撈工具為牡蠣專用拖網(wǎng),該網(wǎng)具網(wǎng)口下杠有鋼制耙,拖網(wǎng)作業(yè)除捕撈活體牡蠣外,也破壞了牡蠣礁生境結(jié)構(gòu),調(diào)查期間觀察到大量牡蠣殼隨活體牡蠣一同被捕撈,生境結(jié)構(gòu)的破壞不僅降低了牡蠣幼體的附著基數(shù)量,也減少了活體牡蠣礁面積,從而對(duì)自然牡蠣礁造成長(zhǎng)期的負(fù)面影響。隨著近兩年的持續(xù)捕撈,該海域牡蠣資源量急劇下降,2018年拖網(wǎng)漁獲中活體牡蠣的比例達(dá)到70%,至2019年調(diào)查時(shí)下降至20%,表明該牡蠣礁的生態(tài)現(xiàn)狀不容樂(lè)觀,亟待加強(qiáng)保護(hù)和修復(fù)。

曹妃甸-樂(lè)亭海域水流較緩,具有廣闊的潮間帶灘涂和潮下帶淺水區(qū),生境類型異質(zhì)多樣,分布有大面積的自然牡蠣礁和海草場(chǎng),水體初級(jí)生產(chǎn)力高,生物資源十分豐富[16]。建議以溯河口海域自然牡蠣礁和海草場(chǎng)為核心區(qū)域,通過(guò)投放小型的適合牡蠣附著生長(zhǎng)的人工礁體(牡蠣殼、混凝土模塊)[31],快速修復(fù)當(dāng)?shù)氐哪迪牻干臣捌渖镔Y源,建設(shè)一個(gè)牡蠣礁國(guó)家級(jí)海洋牧場(chǎng)示范區(qū),重建魚類重要棲息生境,發(fā)揮牡蠣礁和海草場(chǎng)的魚類育幼場(chǎng)功能,養(yǎng)護(hù)海洋漁業(yè)資源,修復(fù)與優(yōu)化海洋生態(tài)系統(tǒng),探索以經(jīng)濟(jì)魚類、蟹類和螺類為產(chǎn)出的海洋牧場(chǎng)利用模式。

致謝：感謝唐山市曹妃甸蠶沙口村魚滿倉(cāng)漁業(yè)捕撈專業(yè)合作社楊連革和朱振江同志在海上調(diào)查采樣中提供的幫助,感謝上海市水產(chǎn)研究所安傳光博士鑒定礁體動(dòng)物樣品。