金洪宇，李 雷，張澤鵬，馬 波

(1.國(guó)家農(nóng)業(yè)科學(xué)漁業(yè)資源環(huán)境撫遠(yuǎn)觀測(cè)試驗(yàn)站，黑龍江哈爾濱 150070；2.農(nóng)業(yè)農(nóng)村部黑龍江流域漁業(yè)資源環(huán)境科學(xué)觀測(cè)實(shí)驗(yàn)站，黑龍江哈爾濱 150070；3.中國(guó)水產(chǎn)科學(xué)研究院黑龍江水產(chǎn)研究所，黑龍江哈爾濱 150070；4.中國(guó)水產(chǎn)科學(xué)研究院營(yíng)口增殖實(shí)驗(yàn)站，遼寧營(yíng)口 115004)

生物完整性是生物集合群落所需的一種穩(wěn)定能力，其具體是指與自身生存環(huán)境區(qū)域相適應(yīng)的生物群落具有維持自身平衡、保持結(jié)構(gòu)完整和適應(yīng)環(huán)境變化的能力[1-2]。魚類作為水生生態(tài)系統(tǒng)的重要組成部分，具有生命周期長(zhǎng)和活動(dòng)范圍廣等特點(diǎn)，能夠準(zhǔn)確地反映長(zhǎng)時(shí)間序列下大尺度范圍內(nèi)水生生態(tài)系統(tǒng)的環(huán)境變化特征[3-5]。1981 年，有學(xué)者提出了生物完整性指數(shù)(index of biotic integrity，IBI)概念[6]，并以魚類作為指示生物，評(píng)價(jià)了美國(guó)中西部河流的健康狀況。目前，生物完整性指數(shù)在底棲生物[7]、藻類[8]和水生植物[9]研究中得到了應(yīng)用。利用生物完整性指數(shù)，以魚類或者大型底棲動(dòng)物作為指示生物，評(píng)價(jià)了長(zhǎng)江中、上游[10-11]、黃河[12]和大清河[13]的健康狀況。由于自然狀況和周邊人類活動(dòng)干擾程度的差異性，生物完整性指數(shù)評(píng)價(jià)方法并不適用于所有河流的健康狀況評(píng)估。因此，需要根據(jù)所研究河流中的生物群落組成和自然環(huán)境特點(diǎn)，構(gòu)建各條河流獨(dú)有的生物完整性指數(shù)評(píng)價(jià)體系。

牡丹江是松花江第二大支流，發(fā)源于長(zhǎng)白山牡丹嶺，全長(zhǎng)726 km，流經(jīng)敦化市、寧安市、海林市、牡丹江市、林口縣和依蘭縣，最終匯入松花江。牡丹江的主要支流有海浪河、珠爾多河、五虎林河和北安河等[14]。牡丹江流經(jīng)地區(qū)的氣候都屬于溫帶大陸性濕潤(rùn)季風(fēng)氣候。在牡丹江中，生活著哲羅鮭(Hucho taimen)和細(xì)鱗鮭(Brachymystaxlenok)等國(guó)家二級(jí)重點(diǎn)保護(hù)野生魚類。牡丹江是北方冷水性瀕危魚類的重要棲息地。

于2021 年5 月和10 月，本研究對(duì)牡丹江及其支流中的魚類資源開展了調(diào)查，構(gòu)建了魚類完整性指數(shù)評(píng)價(jià)體系，對(duì)牡丹江及其支流的健康狀況進(jìn)行了評(píng)價(jià)，旨在豐富北方寒冷地區(qū)河流健康狀況的生物完整性指數(shù)評(píng)價(jià)體系研究，為制定牡丹江的魚類資源保護(hù)和水環(huán)境修復(fù)的相關(guān)政策提供參考依據(jù)。

1 材料與方法

1.1 研究對(duì)象

本研究以牡丹江及其支流(43°43′5′′N至46°17′24′′N，128°22′26′′E至129°57′14′′E)為研究對(duì)象。

1.2 采樣點(diǎn)設(shè)置、魚類樣本采集和種類鑒定

參照《內(nèi)陸水域漁業(yè)自然資源調(diào)查手冊(cè)》[15]，依據(jù)牡丹江流經(jīng)地區(qū)的地形、地貌和水文狀況，共設(shè)置了16 個(gè)采樣點(diǎn)(采樣點(diǎn)S1 至采樣點(diǎn)S16)，其分別位于依蘭鎮(zhèn)、蘭崗鎮(zhèn)、三道通鎮(zhèn)、三道河子鎮(zhèn)、二道河子鎮(zhèn)、鐵嶺鎮(zhèn)、土城子鄉(xiāng)、大孤山、老孤砬子、小石咀子和林江村江段以及鐵嶺河、蛤蟆河、海浪河、雙河和北安河(圖1)。本研究將受人類活動(dòng)干擾很小的采樣點(diǎn)作為參照采樣點(diǎn)，將受人類活動(dòng)干擾很大的采樣點(diǎn)作為受損采樣點(diǎn)。在16個(gè)采樣點(diǎn)中，依蘭鎮(zhèn)江段(采樣點(diǎn)S1)、蘭崗鎮(zhèn)江段(采樣點(diǎn)S2)、老孤砬子江段(采樣點(diǎn)S9)和林江村江段(采樣點(diǎn)S11)的采樣點(diǎn)為參照采樣點(diǎn)，其他采樣點(diǎn)都為受損采樣點(diǎn)。

圖1 牡丹江及其支流中的采樣點(diǎn)分布示意圖Fig.1 The schematic map of sampling points in Mudanjiang River and its tributaries

根據(jù)牡丹江中魚類的種類、生活習(xí)性和季節(jié)分布規(guī)律，采用三層流刺網(wǎng)(用于采集中、上層魚類樣本)和地籠(用于采集底層魚類樣本)相結(jié)合的采捕方法，調(diào)查魚類物種和種群數(shù)量，采集魚類樣本。

在每個(gè)采樣點(diǎn)，布設(shè)4 組三層刺網(wǎng)和2 個(gè)地籠。刺網(wǎng)的規(guī)格為50 m 長(zhǎng)、1 m 高、網(wǎng)眼孔徑為1～4 cm。地籠的規(guī)格為20 m 長(zhǎng)、網(wǎng)眼直徑為5 mm、面積為5 m2。

采樣時(shí)間為2021年5月1日至10日和10月3日至12日的白天，時(shí)長(zhǎng)為12 h。

依據(jù)《東北地區(qū)淡水魚類》[16]、《黑龍江省魚類》[17]和《黑龍江流域(中國(guó))魚類識(shí)別手冊(cè)》[18]，鑒定捕獲到的魚類標(biāo)本的種類，統(tǒng)計(jì)各種魚類物種的數(shù)量。利用量魚板，測(cè)量魚類樣本的體長(zhǎng)(精確到1 mm)。利用電子天平，稱量魚類樣本的體質(zhì)量(精確到0.1 g)。

對(duì)于在采用現(xiàn)場(chǎng)未能鑒別出種類的魚類樣本，使用于澤標(biāo)本保存液進(jìn)行保存，運(yùn)回實(shí)驗(yàn)室后，再進(jìn)行其種類鑒定。

1.3 魚類群落的優(yōu)勢(shì)物種確定和多樣性評(píng)價(jià)方法

利用相對(duì)重要指數(shù)，確定魚類的優(yōu)勢(shì)物種。將相對(duì)重要指數(shù)值大于100 的魚類物種作為優(yōu)勢(shì)物種。

采用Margalef豐富度指數(shù)、Shannon-Wiener多樣性指數(shù)和Pielou 均勻度指數(shù)，研究魚類群落的物種多樣性。

1.4 魚類完整性指數(shù)評(píng)價(jià)方法

1.4.1 候選指標(biāo)及其篩選原則

依據(jù)牡丹江及其支流中特有的魚類組成和生態(tài)系統(tǒng)特征，參照文獻(xiàn)[19]和文獻(xiàn)[20]，將魚類物種組成和豐度、營(yíng)養(yǎng)結(jié)構(gòu)、耐受性、繁殖共位群、魚類數(shù)量和健康狀況共5類30種指標(biāo)作為候選指標(biāo)(表1)，用于構(gòu)建牡丹江的魚類完整性指數(shù)評(píng)價(jià)體系。

表1 魚類完整性指數(shù)評(píng)價(jià)體系構(gòu)建的候選指標(biāo)Table 1 Candidate indicators for the construction of fish integrity index evaluation system

按照以下3方面要求，篩選出牡丹江及其支流的魚類完整性指數(shù)評(píng)價(jià)體系構(gòu)建指標(biāo)。

第一，統(tǒng)計(jì)各采樣點(diǎn)的魚類物種數(shù)量，若都少于3種，則刪除該指標(biāo)；統(tǒng)計(jì)各采樣點(diǎn)魚類物種數(shù)量占物種總數(shù)量的比例，若都小于10%，則刪除該指標(biāo)；對(duì)于各采樣點(diǎn)的任意1 個(gè)指標(biāo)值，若90%以上都為0，則刪除該指標(biāo)。

第二，利用箱體圖，對(duì)參照采樣點(diǎn)與受損采樣點(diǎn)不同指標(biāo)的25%～75%分位數(shù)范圍進(jìn)行判讀，將不存在重疊顯現(xiàn)或者各自中位數(shù)未在對(duì)方箱體之內(nèi)出現(xiàn)部分重疊現(xiàn)象的指標(biāo)，予以保留。

第三，采用判別分析方法，篩選保留后的指標(biāo)，進(jìn)行指標(biāo)值間的Pearson 相關(guān)分析，剔除顯著相關(guān)(相關(guān)系數(shù)的絕對(duì)值大于0.9)的指標(biāo)，保留最能反映結(jié)果的指標(biāo)。

利用Excel 2016軟件和R語(yǔ)言軟件，對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。

1.4.2 指標(biāo)賦值和評(píng)價(jià)標(biāo)準(zhǔn)

采用比值法[21]，統(tǒng)一各指標(biāo)的量綱。在外界干擾影響下數(shù)值減小的指標(biāo)，其最佳期望值為95%分位數(shù)，指標(biāo)的分值等于實(shí)際值/最佳期望值；反之，其最佳期望值為5%分位數(shù)，指標(biāo)的分值等于最大值與實(shí)際值之差除以最大值與最佳期望值之差。若計(jì)算出的指標(biāo)分值大于1，則都計(jì)為1。

當(dāng)魚類完整性指數(shù)值大于或者等于參照采樣點(diǎn)25%分位數(shù)值時(shí)，其所在河流的健康狀況等級(jí)為“健康”；采用四分法，將河流的健康狀況向下依次劃分為“亞健康”、“一般”、“差”和“極差”。

2 結(jié)果與分析

2.1 魚類群落的物種組成

在調(diào)查中，共采集了2 559 尾魚類標(biāo)本，其總質(zhì)量為78.55 kg；記錄了7 目13 科37 屬49 種魚類(表2)。其中，鯉形目、鱸形目、鮭形目、鲇形目、鲉形目、鱈形目和七鰓鰻目的科、屬、種數(shù)量分別為2科25屬36種、3科4屬4種、2科2屬3種、3科3屬3種、1科1屬1種、1科1屬1種和1科1屬1種。

表2 調(diào)查期間牡丹江及其支流中的魚類物種名錄Table 2 List of fish species in the Mudanjiang River and its tributaries during the investigation

在49 種魚類中，鯽的分布范圍最廣，在16 個(gè)采樣點(diǎn)分布(表3)；其次為麥穗魚，在12 個(gè)采樣點(diǎn)分布；再其次為犬首鮈，在11 個(gè)采樣點(diǎn)分布；?和大鰭鱊的分布也比較廣泛；而花羔紅點(diǎn)鮭、細(xì)鱗鲴、條紋似白鮈、平口鮈、東北鳈、高體鮈、棒花魚和花斑副沙鰍都僅在1個(gè)采樣點(diǎn)分布，屬于牡丹江及其支流中的少見物種。

表3 調(diào)查期間出現(xiàn)魚類物種的采樣點(diǎn)Table 3 The sampling points of fish species occurred during the investigation

在牡丹江及其支流中，有15 種魚類為優(yōu)勢(shì)物種，其分別為鯽、麥穗魚、鲇、光澤黃顙魚、大鰭鱊、?、犬首鮈、黑龍江泥鰍、黑龍江鳑鲏、黃顙魚、高體鳑鲏、細(xì)體鮈、蒙古鲌、北方須鰍和大銀魚(表4)。其中，鯽的相對(duì)重要指數(shù)值、種群數(shù)量所占比例和質(zhì)量所占比例都最大。

表4 牡丹江及其支流中魚類優(yōu)勢(shì)物種的相對(duì)重要指數(shù)值及其種群數(shù)量和質(zhì)量所占比例Table 4 The values of relative importance index of dominant fish species in the Mudanjiang River and its tributaries and the proportions of their population quantities and weights

2.2 魚類群落的物種多樣性

在16個(gè)采樣點(diǎn)中，參照采樣點(diǎn)S1魚類群落的Shannon-Wiener多樣性指數(shù)值(4.06)和Margalef豐富度指數(shù)值(3.02)都最大，其Pielou 均勻度指數(shù)值(0.84)也較大，說(shuō)明參照采樣點(diǎn)S1魚類群落的多樣性水平較高；受損采樣點(diǎn)S6 魚類群落的Pielou 均勻度指數(shù)值(0.89)最大，其次是受損采樣點(diǎn)S7魚類群落的Pielou均勻度指數(shù)值(0.88)(表5)。

表5 各采樣點(diǎn)魚類群落的多樣性指數(shù)值Table 5 The values of diversity indexes of the fish communities in the sampling points

2.3 牡丹江及其支流的健康狀況

2.3.1 篩選出的指標(biāo)

根據(jù)候選指標(biāo)篩選標(biāo)準(zhǔn)，將指標(biāo)M24、指標(biāo)M25和指標(biāo)M30從候選指標(biāo)中剔除，還剩下27個(gè)指標(biāo)。對(duì)參照采樣點(diǎn)的27 個(gè)指標(biāo)(表6)的逐步判別分析結(jié)果顯示，指標(biāo)M1、指標(biāo)M2、指標(biāo)M3、指標(biāo)M10、指標(biāo)M13、指標(biāo)M14、指標(biāo)M15、指標(biāo)M17、指標(biāo)M18、指標(biāo)M20、指標(biāo)M23、指標(biāo)M27和指標(biāo)M28 的分值都大于2，這13 個(gè)候選指標(biāo)為篩選出的進(jìn)行相關(guān)分析的指標(biāo)。

表6 參照采樣點(diǎn)候選指標(biāo)的度量值Table 6 The metric values of candidate indicators for the reference sampling points

相關(guān)分析結(jié)果顯示，指標(biāo)M1、指標(biāo)M2、指標(biāo)M3、指標(biāo)M18、指標(biāo)M20 和指標(biāo)M28 兩兩之間顯著相關(guān)(表7)，指標(biāo)M20為敏感性魚類物種數(shù)量所占比例，其所包含的信息量相對(duì)最多，因此只保留指標(biāo)M20，而將其他5 個(gè)指標(biāo)剔除。候選指標(biāo)M13、指標(biāo)M14、指標(biāo)M15、指標(biāo)M17、指標(biāo)M23和指標(biāo)M27 兩兩之間顯著相關(guān)，指標(biāo)M13 為植食性魚類種群數(shù)量所占比例，其所包含的信息量相對(duì)最多，因此只保留指標(biāo)M13，而將其他5 個(gè)指標(biāo)剔除。指標(biāo)M10(中下層魚類物種數(shù)量所占比例)與任何候選指標(biāo)都不相關(guān)。最終，選取指標(biāo)M10、指標(biāo)M13 和指標(biāo)M20 作為魚類完整性指標(biāo)，構(gòu)建魚類完整性指數(shù)，對(duì)牡丹江及其支流的健康狀況進(jìn)行評(píng)價(jià)。

表7 13個(gè)指標(biāo)間的相關(guān)系數(shù)Table 7 The correlation coefficients between the 13 indicators

2.3.2 魚類完整性指標(biāo)的賦值

指標(biāo)M10的分值等于中下層魚類物種數(shù)量所占比例除以51.39。指標(biāo)M13的分值等于66.41減去植食性魚類種群數(shù)量所占比例再除以66.41。指標(biāo)M20的分值等于敏感性魚類物種數(shù)量所占比例除以28.19。各采樣點(diǎn)所在江段或支流健康狀況的魚類完整性指數(shù)評(píng)價(jià)分值等于上述3 個(gè)魚類完整性指標(biāo)分值之和。

2.3.3 各采樣點(diǎn)所在江段和支流的健康狀況

將參照采樣點(diǎn)魚類完整性指數(shù)值的25%分位數(shù)值(1.97)作為劃分河流健康狀況等級(jí)的閾值。當(dāng)河流的魚類完整性指數(shù)值≥1.97時(shí)，表示河流處于健康狀態(tài)；將小于25%分位數(shù)的魚類完整性指數(shù)值進(jìn)行四分法分級(jí)，當(dāng)1.48≤魚類完整性指數(shù)值＜1.97 時(shí)，表示河流處于亞健康狀態(tài)；當(dāng)0.99≤魚類完整性指數(shù)值＜1.48 時(shí)，表示河流的健康狀況一般；當(dāng)0.49≤魚類完整性指數(shù)值＜0.99 時(shí)，表示河流的健康狀況差；當(dāng)魚類完整性指數(shù)值小于0.49時(shí)，表示河流的健康狀況極差。

由表8可知，蘭崗鎮(zhèn)江段、土城子鄉(xiāng)江段、大孤山江段、老孤砬子江段和林江村江段處于健康狀態(tài)；依蘭鎮(zhèn)江段、三道通鎮(zhèn)江段和蛤蟆河處于亞健康狀態(tài)；三道河子鎮(zhèn)江段、二道河子鎮(zhèn)江段、鐵嶺鎮(zhèn)江段、小石咀子江段和海浪河的健康狀態(tài)一般；鐵嶺河、雙河和北安河的健康狀況差。

表8 牡丹江及其支流的健康狀況Table 8 The health status of the Mudanjiang River and its tributaries

3 討 論

以往的研究表明，牡丹江中魚類的物種數(shù)量為37 種[22]。在本研究中，在牡丹江及其支流中發(fā)現(xiàn)了49種魚類，但是，其中有11種為外來(lái)物種，未發(fā)現(xiàn)哲羅鮭、細(xì)鱗鮭、黑龍江茴魚(Thymallus grubii)和烏蘇里白鮭(Coregonus ussuriensis)等珍稀土著物種，說(shuō)明牡丹江及其支流中的魚類完整性水平有所下降。近年來(lái)，隨著人們對(duì)水產(chǎn)品需求量的增大，大量外來(lái)魚類物種被引入當(dāng)?shù)仞B(yǎng)殖，在魚類養(yǎng)殖過(guò)程中，因偶然逃逸和直接放生等原因，使大銀魚、花羔紅點(diǎn)鮭等外來(lái)物種進(jìn)入了牡丹江及其支流水域，使魚類物種數(shù)量增多。另外，通過(guò)對(duì)比發(fā)現(xiàn)，以往相關(guān)研究所采用的調(diào)查方法與本研究的調(diào)查方法存在差異，這也許是導(dǎo)致本研究記錄的牡丹江及其支流中魚類物種數(shù)量增多的原因之一。牡丹江及其支流周邊縣市的產(chǎn)業(yè)主要是工業(yè)和畜牧業(yè)，有大量未經(jīng)有效處理的工業(yè)和畜牧業(yè)廢水被直接排入牡丹江及其支流中，導(dǎo)致水體中化學(xué)需氧量和氨氮含量等超標(biāo)，污染了哲羅鮭、細(xì)鱗鮭、黑龍江茴魚和烏蘇里白鮭等冷水性魚類的棲息地，使這些物種逐漸絕跡[23-24]。大銀魚、花羔紅點(diǎn)鮭等外來(lái)物種與同屬于鮭科的哲羅鮭、細(xì)鱗鮭、黑龍江茴魚和烏蘇里白鮭之間存在生態(tài)位競(jìng)爭(zhēng)關(guān)系，也導(dǎo)致土著物種資源量減少甚至消失。

魚類多樣性能直接體現(xiàn)河流的健康狀況，其與可以被魚類物種利用的生態(tài)位、生存環(huán)境密切相關(guān)[25]。本研究的牡丹江中魚類群落的Shannon-Wiener多樣性指數(shù)值(1.69～4.06)小于以往研究的相關(guān)數(shù)據(jù)(2.16～4.29[22])，說(shuō)明牡丹江中魚類群落的穩(wěn)定性有變差的趨勢(shì)。在本研究的16個(gè)采樣點(diǎn)中，牡丹江依蘭鎮(zhèn)江段魚類群落的Shannon-Wiener多樣性指數(shù)值最大，而鐵嶺河、雙河和北安河中魚類群落的Shannon-Wiener多樣性指數(shù)值都較小。總體上，牡丹江干流中魚類群落的Margalef豐富度指數(shù)、Shannon-Wiener多樣性指數(shù)和Pielou均勻度指數(shù)值較大。Margalef 豐富度指數(shù)和Pielou均勻度指數(shù)能夠反映群落結(jié)構(gòu)的穩(wěn)定性[26]，指數(shù)值越小，表明群落結(jié)構(gòu)的復(fù)雜度和穩(wěn)定性越差，說(shuō)明牡丹江干流中的魚類群落結(jié)構(gòu)相對(duì)更合理。鐵嶺河、雙河和北安河等牡丹江支流水域狹窄，容易出現(xiàn)斷流現(xiàn)象，導(dǎo)致魚類物種可以利用的生態(tài)位和生存環(huán)境萎縮，使魚類群落的多樣性水平降低。魚類多樣性與餌料生物豐度密切相關(guān)[27-28]，而牡丹江支流中的餌料生物豐度遠(yuǎn)低于牡丹江干流，導(dǎo)致牡丹江支流中魚類的攝食、生長(zhǎng)和繁衍受到限制。

利用各采樣點(diǎn)的魚類完整性指數(shù)值，評(píng)價(jià)了牡丹江及其支流的健康狀況，結(jié)果顯示，鐵嶺河、雙河和北安河的健康狀況差，這些河流中的魚類群落承受著很大的環(huán)境壓力。北安河附近分布著大量的工廠，部分未達(dá)標(biāo)的工廠廢水被排入河流中，導(dǎo)致河流水體被污染。蛤蟆河的水體主要被隨雨水流入河中的農(nóng)業(yè)和畜牧業(yè)污染物和廢水所污染。近年來(lái)，牡丹江上游的三道河子鎮(zhèn)江段的旅游業(yè)發(fā)展迅速，已建造了蓮花湖等旅游景點(diǎn)，該江段的河道被水渠化，破壞了產(chǎn)黏性卵的魚類物種的生存環(huán)境。牡丹江上游的二道河子鎮(zhèn)江段和鐵嶺鎮(zhèn)江段的水環(huán)境也存在污染加劇的現(xiàn)象。調(diào)查發(fā)現(xiàn)，這兩處江段的沿岸共有27 家重點(diǎn)工業(yè)企業(yè)和37 家畜禽養(yǎng)殖場(chǎng)，每日還有大量生活污水產(chǎn)生，這些工業(yè)企業(yè)、畜禽養(yǎng)殖場(chǎng)和居民的生活污水被排入了這兩處江段，導(dǎo)致江水的水質(zhì)下降，魚類的生存環(huán)境遭到破壞。綜上所述，污水排放、開發(fā)旅游業(yè)等人類活動(dòng)是造成牡丹江部分江段和支流的健康狀況較差的主要原因，本研究對(duì)牡丹江及其支流健康狀況的評(píng)價(jià)結(jié)果與實(shí)際情況相符合。

4 結(jié) 論

2021年，牡丹江干流上游的少數(shù)江段(蘭崗鎮(zhèn)江段)、寧江市以南的絕大多數(shù)干流江段(土城子鄉(xiāng)江段、大孤山江段、老孤砬子江段和林江村江段)處于健康狀態(tài)；上游的多數(shù)江段(依蘭鎮(zhèn)江段、三道通鎮(zhèn)江段)和牡丹江支流蛤蟆河處于亞健康狀態(tài)；牡丹江干流中游的江段(三道河子鎮(zhèn)江段、二道河子鎮(zhèn)江段、鐵嶺鎮(zhèn)江段)和牡丹江支流海浪河的健康狀況一般；牡丹江支流鐵嶺河、雙河和北安河的健康狀況差。

建議嚴(yán)格控制污染牡丹江及其支流水體的污染源，加大對(duì)生活污水和企業(yè)廢水排放的監(jiān)管力度，使牡丹江及其支流的水環(huán)境恢復(fù)健康；禁止開發(fā)和利用魚類產(chǎn)卵場(chǎng)所在江段的水資源；采取人與自然和諧相處的旅游開發(fā)模式，避免旅游業(yè)污染和破壞牡丹江及其支流的水環(huán)境。