昌黎海域文昌魚時空分布及其與棲息地底質關系的探究＊

周衛彪，吳 凱，高偉明，韓曉慶

（1.河北師范大學 資源與環境科學學院，河北 石家莊，050024；2.河北省國土資源利用規劃院，河北 石家莊，050051）

文昌魚（Branchiostomalanceolatum）為國家二級野生保護動物，是一種珍貴的海洋生物，青島文昌魚是廈門文昌魚的一個變種，是黃渤海唯一一種文昌魚，分布在秦皇島近岸、煙臺、青島膠州灣及大公島附近海域。灤河口外昌黎海域為國內青島文昌魚棲息密度最高的海域之一。文昌魚在砂質中營半穴居生活，其分布與沉積物粒度有密切的關系，研究清楚文昌魚的分布與棲息地底質結構的關系，對加強文昌魚的保護有積極意義［1－7］。曹玉萍等曾對該海域的文昌魚種群結構和生物學特征進行研究［3］，馬明輝曾于2000－08對該海域進行調查，并對文昌魚與底質的關系做出分析［4］。但是，關于該海域多年間的文昌魚分布及棲息地底質的時間、空間變化及文昌魚與棲息地底質之間的關系少有研究。運用該海域2005年、2008年及2011年3期文昌魚分布及底質調查數據，對該海域的文昌魚分布及底質結構時空變化進行了研究，分析了該研究海區的文昌魚分布和棲息地底質特征存在的關系，并對該海域威脅文昌魚生存的主要環境壓力進行了分析。

1 研究區概況

昌黎縣位于河北省東北部，隸屬秦皇島市。昌黎縣海域面積757.88km2，其中淺海面積715.86km2，占海域總面積的94.46%。灤河口外昌黎海域北起大蒲河河口南岸，南至灤河河口北岸（圖1），屬暖溫帶濕潤大陸性季風氣候，具有光照充足、降水豐沛、雨熱同期等特點。研究區內主要河流為灤河，其余河流均為季節性河流。

圖1 研究區位置圖Fig.1 Location of the study area

2 研究數據及方法

2.1 調查范圍及采樣站位

調查范圍為灤河口以北至大蒲河口以南，高潮線至20m等深線區域內近350km2的海域。在上述區域內共設站位19個（圖2）。采用2005－08，2008－08及2011－08三期文昌魚和底質監測數據。

2.2 分析測試方法

分別 于 2005－08－15，2008－08－25 和 2011－08－04 利 用0.05m2抓斗式采泥器，進行文昌魚及沉積物樣品采集。文昌魚樣本的分析測試方法采用計數法，樣本處理按照相關規范［8］進行。沉積物處理采用篩分法結合沉析法即綜合法，取2份分析樣本，其中一份供作篩分法，分析粒徑大于0.063mm的粒級配比，另一份供作沉析法，分析粒徑小于0.063mm的粒級配比。最后，2種方法所求的體積分數，統一平差，求出校正值。粒度參數統計采用Folk等公式計算，沉積類型劃分采用Shepard三角圖式［9－10］。

圖2 監測站位分布圖Fig.2 Locations of the monitoring stations

3 結果與分析

3.1 沉積物類型

監測海區的底質沉積物粒度范圍在0.001～2.000mm，以砂為主。其中，2005年監測海區底質沉積物類型：砂－粉砂－粘土（STY）、極細砂（VFS）、粘土質砂（YS）、粉砂質砂（TS）、粘土質粉砂（YT）、細砂（FS）、中細砂（MFS）、細中砂（FMS）、中粗砂（MCS）9種。2008年監測海區底質沉積物類型：粘土質粉砂（YT）、砂質粉砂（ST）、粉砂質砂（TS）、細砂（FS）、中細砂（MFS）、細中砂（FMS）、中砂（MS）7種。2011年監測海區底質沉積物類型：砂（S）、粉砂質砂（TS）、粘土質粉砂（YT）三種。3期監測期間，各粒組平均比例變化（表1）。

表1 監測期各粒組體積分數變化Table 1 Changes in the proportions of particle factions in the monitoring period

3.2 粒度特征

監測海區沉積物粒度特征變化見表2。其中，2005年沉積物粒度分選性中等，2008年分選性差，2011年分選性好。

表2 監測海區沉積物粒度特征Table 2 Grain－size characteristics of the sediments in the monitoring sea area

3.3 文昌魚的分布變化

2005年文昌魚的分布主要分布在以大蒲河口外18號站與新開口外7號站連線為軸心分布，其中高密度區為11和10號站，11號站位的文昌魚分布達到475尾/m2；2008年文昌魚的分布分為3個部分，灤河口近岸區域、大蒲河口至新開口的近岸區域及遠岸區域，其中高密度區為11和14號站；2011年文昌魚的分布范圍及密度均下降明顯，只在灤河口近岸區域及大蒲河口至新開口的遠岸區域少量分布。3期監測中各站位文昌魚密度及生物量見表3。

表3 各站位文昌魚密度及生物量Table 3 The density and biomass of Amphioxus at each stations

3.4 文昌魚的分布與底質的關系

文昌魚對底質的要求相對嚴格，本調查區域的溫度、鹽度、深度等均在文昌魚的適宜范圍內［11］。

3.4.1 底質類型對文昌魚的影響

監測數據顯示，文昌魚集中分布的底質類型為細砂、中細砂、細中砂及中粗砂。利用各站點的底質粒組體積分數及文昌魚密度數據，運用SPSS 19.0軟件將各粒組體積分數與文昌魚密度和生物量做聚類分析，發現0.125～1.000mm粒組比例與文昌魚的密度和生物量關系密切（圖3）。研究區底質的粒徑主要為0.125～0.500mm，2005年該粒組的體積分數為72.95%～92.95%，2008年為65.66%～95.49%。2011年為83.60%～97.08%，而0.500～1.000mm的粒組值較小。做相關性分析發現0.500～1.000 mm粒組與文昌魚分布呈顯著相關（r＝0.360，p＜0.01，n＝57）。數據表明，0.500～1.000mm粒組值均高于平均值，以文昌魚分布最多的站位（2005年11號站）為例，該站位0.500～1.000mm粒組的體積分數為5.52%，遠高于2005年該粒組在研究區分布的平均值1.62%。因此，初步推斷0.500～1.000mm粒組也是影響該海區文昌魚分布的一個重要粒組組成。文昌魚選擇這種粗細混摻的底質，是因為此類底質較疏松，通氣性好，便于其潛居生活。

圖3 文昌魚密度、生物量與各粒組比例聚類分析圖Fig.3 Hierarchical clustering analysis of the density and the biomass of Amphioxus and the individual particle fractions

文昌魚主要分布在0.125～1.000mm粒組大于90%的區域，2005年該組分體積分數大于90%的區域主要分布在大蒲河口外；2008年該組分體積分數高于90%的區域主要分布在大蒲河口外，呈條帶狀分布；2011年該組分體積分數高于90%的區域主要分布在大蒲河口至新開口的遠岸區域及灤河口近岸區域（圖4）。圖4中曲線表示0.125～1.000mm粒組的體積分數（%），該粒組體積分數＞90%的分布變化與文昌魚的分布變化具有一致性。

圖4 研究區0.125～1mm粒組分布與文昌魚分布關系示意圖Fig.4 Relationship between the distribution of 0.125～1mm fraction and that of Amphioxus in the study area

在2000年馬明輝探討該海域文昌魚棲息地底質特征時曾指出，影響本地區的文昌魚分布的主要粒組為0.125～0.500mm［4］。本研究發現，灤河口外昌黎海域文昌魚的適宜棲息粒徑范圍有所擴大，這顯示出文昌魚對環境變化的適應性。同時與黃厝海域文昌魚適宜粒徑［5］、閩江口至南澳海岸帶的文昌魚適宜粒徑［6］、大放雞海域文昌魚適宜粒徑［7］相比較（表4），灤河口外昌黎海域的文昌魚適宜粒徑相對較小，這很可能是由昌黎海域文昌魚群體與黃厝海域、閩江口至南澳海域、大放雞海域文昌魚群體差異造成的，也顯示了不同海域文昌魚對不同底質粒度的適應。

表4 各地文昌魚適宜粒組比較Table 4 Comparison of the particle fractions that favor the inhabitance of Amphioxus at different locations

3.4.2 底質粒度特征對文昌魚分布的影響

沉積物的粒度特征是沉積物的基本性質之一，是沉積物物源、沉積動力、搬運距離等綜合作用的結果，在沉積物粒度中包含了水動力條件和物源的信息［12］。中值粒徑指示了沉積物粒度分布的集中趨勢，反映了沉積介質的平均動能。研究區的文昌魚集中分布區中值粒徑變化見表5。文昌魚集中分布區沉積顆粒粗細狀況出現先細化后粗化的波動，即在2008年偏細，其中文昌魚的分布密度大于50尾/m2的站位平均粒徑1.38～2.98φ（圖5），在2011年偏粗，而這種波動主要受研究海域水動力的增減和沉積物物質來源的變化的影響。圖5中曲線表示中值粒徑，單位為φ。

表5 文昌魚集中分布海區底質粒度特征Table 5 The grain－size characteristics of the bottom sediments in the areas where Amphioxus is concentrated

圖5 研究區中值粒徑與文昌魚分布關系示意圖Fig.5 Relationship between the medium diameters of sediments and the distribution of Amphioxus in the study area

分選系數代表了沉積物顆粒大小的均勻程度，若粒級少，主要粒級突出，則分選就好，分選系數的數值小；反之，粒級分布范圍很廣，主要粒級不突出，則分選就差，分選系數的數值大。研究區的文昌魚集中分布區分選系數變化（表5）：2005年沉積物粒度分選性好；2008年分選性中等；2011年分選性很好。除2008年1、8、10號站位及2005年15號站位分選系數大于2，分選性差，其余站位分選系數均小于2，其中文昌魚密度大于50尾/m2的站位，分選系數均小于1，分選性為好或很好（圖6）。圖6中曲線表示分選系數。

圖6 研究區分選系數與文昌魚分布關系示意圖Fig.6 Relationship between the sorting coefficients of sediments and the distribution of Amphioxus in the study area

2008年研究海域分選系數為2.17，分選性差，說明研究海域水動力過程波動幅度較大，沉積環境不穩定，受此影響，在2008年文昌魚的集中分布區，分選系數的數值也較2005年、2011年大，分選性中等。偏態度用來量度沉積物顆粒頻率分布的對稱程度，反應沉積環境能量的變異。研究海域中的文昌魚集中分布區偏態度變化見表5。四分位偏態均近似對稱分布，呈正偏。

3.5 文昌魚種群數量下降原因分析

首先，文昌魚生境退化的主要原因是灤河入海懸沙量明顯減少，沙源補充嚴重不足所致。文昌魚對其生活沙質環境要求很高。灤河徑流量和輸沙量的季節變化極為明顯，隨之作用營力也發生顯著的季節變化。在洪水季節，以河流作用為主，表現為三角洲的堆積；在枯水季節，以波浪作用為主，表現為三角洲被破壞和改造。在上述條件的作用下，灤河口以北海域，沉積特征表現為近岸為細砂、中細砂，向外出現大片中砂，其范圍可至水深10～15m的平緩地帶，形成了文昌魚適宜的棲息環境。

灤河水利工程（包括干流的潘家口水庫、大黑汀水庫和支流上的桃林口水庫及引灤工程）的實施控制了灤河流域89.9%的匯水面積，調蓄了灤河下游63.6%的徑流量，也使其泥沙量減少了99.3%［13］，從而引起灤河下游直至灤河三角洲的快速響應，由此導致灤河口的衰亡和灤河三角洲由淤漲轉變為侵蝕后退。

近幾年來灤河三角洲的形態發生了嚴重改變，已近衰亡［14］。三角洲在衰亡過程中前緣在侵蝕狀態下逐漸變得平直，濱外壩向陸遷移，外廓線曲率愈來愈小；三角洲兩側外廓線曲率較頂部相對較大。因文昌魚棲息地漲潮流速大于落潮流速，中砂及中細砂被不斷搬運到潮間帶，再在由風動力搬運補充給沿岸沙丘，導致文昌魚棲息地中砂和中細砂含量下降，棲息地發生嚴重改變，文昌魚種群數量下降。

同時，昌黎海域近年來海洋養殖規模面積猛增，規模不斷擴大（表6），對文昌魚群的生存和繁衍有一定的威脅。養殖密度過大，養殖總量超過水域容納量，養殖病害在逐漸增多，養殖貝類糞大量沉降到海底也加劇了文昌魚棲息環境的改變。另外由于水產養殖大量的圍墾，造成水動力條件下降，設置的大量養殖筏及吊籠會改變水動力條件［15］，也對文昌魚的棲息地造成一定的影響。育苗室、儲水池等建設對海灘造成嚴重的破壞，向海內排放養殖廢水，水體中N、P等含量升高，臨近水域營養鹽含量不斷增加，加速了藻類在春夏季的大量繁殖和水體的富營養化［16］，導致文昌魚的密度和生物量呈現下滑趨勢。

表6 昌黎縣歷年海水養殖面積統計Table 6 The statistics of the ocean aquaculture areas over the years in the Changli County

此外，灤河口昌黎海域南為灤河口，北為大蒲河口，每年都有大量污染物由河口排入。灤河與大蒲河排入海區的主要污染物為有機污染物（以COD計）和石油類。大蒲河口入海有機污染物占總污染物的97%，石油烴占2.8%，灤河入海有機污染物占總污染物的97.5%，石油烴為2.2%。灤河水系，收納的廢污水占河北省濱海三市排污口進入入海河流的污廢水量的4.14%，污染物占8.50%。由于文昌魚的濾食活動有助于凈化水體，其數量多的水域通常是清潔的，但文昌魚對環境污染極其敏感，當水環境污染的程度大于文昌魚凈化水體的能力時，不利于文昌魚的生長。

4 結論

1）灤河口外昌黎海域文昌魚主要分布在細砂、中細砂、細中砂及中粗砂類型的底質中，中值粒徑范圍為1.38～2.98φ，分選性好，四分位偏態近似對稱，呈正偏；

2）影響文昌魚分布的主要粒組為0.125～1.000mm，文昌魚主要分布在該粒組體積分數＞90%的海域，該粒組體積分數＞90%的底質分布變化與文昌魚的分布變化具有一致性，同時發現該海域文昌魚的適宜棲息粒組范圍擴大；

3）與其他海域的適宜文昌魚棲息底質對比，發現昌黎海域的適宜文昌魚棲息底質的顆粒較細；

4）文昌魚的空間分布由片狀、條帶狀分布向斑點狀分布演變，文昌魚密度及生物量也呈下降趨勢，特別是在2008～2011年監測期，下降趨勢明顯；

5）該海域文昌魚種群退化的主要原因是是灤河入海懸沙量明顯減少，沙源補充嚴重不足，同時該海域海洋養殖產業的迅猛發展，改變了文昌魚適宜生存和繁衍的底質和水動力條件，成為文昌魚種群退化的一個重要原因，此外，該海域的環境污染也對文昌魚的生存造成一定威脅。

灤河口外昌黎海域文昌魚的棲息環境已遭受破壞，文昌魚種群在該海域出現嚴重的退化，我們研究了該海域底質與文昌魚的時空變化及相關關系，但是，對文昌魚的棲息環境僅僅分析底質這單一因素，顯然不足以了解文昌魚的生命生態全貌，仍需要在今后的研究對其它因素諸如水質環境、餌料等加以綜合考慮。

致謝：河北省昌黎黃金海岸國家級自然保護區管理處為本研究提供相關數據。

（References）：

［1］YUAN W，WANG J，LIN Q，et al.The distribution and habitat use of Branchiostoma belcheri at the mouth of Jiaozhou Bay，Qingdao in autumn［J］.Biodiversity Science，2011，19（4）：470－475.袁偉，王俊，林群，等.青島膠州灣口海域秋季文昌魚的分布及棲息環境的特征［J］.生物多樣性，2011，19（4）：470－475.

［2］YANG J W，REN Y P，XU B D.Study on living resources and amphioxus resources in Qingdao natural reserve area abstract［J］.Journal of Hydroecology，2008，28（3）：112－114.楊建威，任一平，徐賓鐸.青島自然保護區文昌魚的資源調查與生長研究［J］.水利漁業，2008，28（3）：112－114.

［3］CAO Y P，YAN L N，XIE S，et al.Preliminary investigation of Amphioxus in Changli［J］.Chinese Journal of Zoology，2001，36（3）：10－13.曹玉萍，閆路娜，謝松，等.昌黎海區文昌魚初步調查［J］.動物學雜志，2001，36（3）：10－13.

［4］MA M H，ZHANG Z N，FENG Z Q，et al.The characteristics of sediments and distribution of Branchiostoma belcheritsingtaneuse in Luanhe Estuary［J］.Marine Environmental Science，2005，24（2）：39－42.馬明輝，張志南，馮志權，等.灤河口青島文昌魚分布與棲息地底質特征［J］.海洋環境科學，2005，24（2）：39－42.

［5］ZHOU Q L，HE M H，SHAO H D，et al.Distribution of lancelet and its relationship to bottom soil in Huangcuo waters of Xiamen［J］.Journal of Oceanography in Taiwan Strait，1990，9（2）：185－189.周秋麟，何明海，邵合道，等.廈門黃厝海區文昌魚的分布及其與底質的關系［J］.臺灣海峽，1990，9（2）：185－189.

［6］WANG W Y，CHEN B Z，YAO L T，et al.Distribution investigation on amphioxus living in Minjiang Estuary to Nan′ao Island Coastal Zone［J］.Journal of Fujian Fisheries，1989，（1）：14－16.汪偉洋，陳必哲，姚聯騰，等.閩江口至南澳海岸帶文昌魚分布調查［J］.福建水產，1989，（1）：14－16.

［7］FENG B，LU H S，CHEN W H.Distribution of Branchiostoma belcheni and its relation with bottom sediments in waters around Dafangji Island of Maoming，Guangdong［J］.Journal of Tropical Oceanography，2006，25（2）：18－23.馮波，盧伙勝，陳文河.廣東茂名大放雞海區文昌魚分布與底質關系［J］.熱帶海洋學報，2006，25（2）：18－23.

［8］GB17378.7－2007The Specification for Marine Monitoring［S］.Beijing：China Ocean Press，2008.GB 17378.7－2007海洋監測規范［S］.北京：海洋出版社，2008.

［9］SHEPARD F P.Nomenclature based on sand－silt－clay rations［J］.J.Sed.Petrol.，1954，24：151－158.

［10］FOLK R L，WARD W C.Brazos river Bar：A study in the significance of grain size parameters［J］.J.Sed.Petrol.，1957，27（1）：3－26.

［11］FANG Y Q.The ecology habits of Amphioxus and the protection about its resource［J］.Chinese Journal of Zoology，1987，22（2）：41－45.方永強.文昌魚生態習性及其資源的保護［J］.動物學雜志，1987，22（2）：41－45.

［12］YANG X H，FENG X L，CHU Z X，et al.Surface sediment distribution and sedimentary environment on the East China continental shelf［J］.Periodical of Ocean University of China：Natural Science edition，2012，42（7－8）：126－134.楊旭輝，馮秀麗，褚忠信，等.中國東部陸架表層沉積物粒度特征及其沉積環境淺析［J］.中國海洋大學學報：自然科學版，2012，42（7－8）：126－134.

［13］ZHANG L Y，CAI Z，LI Q C，et al.The Runoff－silt Relationship in Luanhe River Basin and Effects of the Water Conservancy Projects on the Lower Reaches［J］.Hubei Agricultural Sciences，2012，51（15）：3222－3225.張麗云，蔡湛，李慶辰，等.灤河的水沙關系及水利工程對下游河道的影響［J］.湖北農業科學，2012，51（15）：3222－3225.

［14］FENG J L，ZHANG W.Geometric Features of Modern Luanhe River Delta Evolutton［J］.Adv.Mar.Sci.，1997，15（3）：22－25.馮金良，張穩.灤河現代三角洲演變的幾何學特征［J］.黃渤海海洋，1997，15（3）：22－25.

［15］SU Y，LIU J，HAN X Q，et al.Effects of Mariculture on Coastal ecological Environment－a case study of Changli county in Hebei province［J］.Jiangsu Agricultural Sciences，2012，40（3）：306－309.蘇藝，劉佳，韓曉慶，等.海水養殖對海洋生態環境的影響——以河北省昌黎縣為例［J］.江蘇農業科學，2012，40（3）：306－309.

［16］DONG S L，PAN K H，UWE B.Review on Effects of Mariculture on Coastal Environment［J］.Journal of Ocean University of Qingdao：Natural Science edition，2000，30（4）：575－582.董雙林，潘克厚，Uwe Brockmann.海水養殖對沿岸生態環境影響的研究進展［J］.青島海洋大學學報：自然科學版，2000，30（4）：575－582.