丁隆強(qiáng)，何曉輝，李新豐，方弟安，，徐東坡，

(1:上海海洋大學(xué)水產(chǎn)科學(xué)國(guó)家級(jí)實(shí)驗(yàn)教學(xué)示范中心，上海 201306)

(2:中國(guó)水產(chǎn)科學(xué)研究院淡水漁業(yè)研究中心，農(nóng)業(yè)農(nóng)村部長(zhǎng)江下游漁業(yè)資源環(huán)境科學(xué)觀測(cè)實(shí)驗(yàn)站，無(wú)錫 214081)

青魚(yú)(Mylopharyngodonpiceus)、草魚(yú)(Ctenopharynodonidellus)、鰱(Hypophthalmichthysmolitrix)和鳙(Aristichthysnobilis)統(tǒng)稱為四大家魚(yú)，是我國(guó)特有的經(jīng)濟(jì)魚(yú)類，其產(chǎn)量曾占我國(guó)淡水魚(yú)總產(chǎn)量的60%[1]. 長(zhǎng)江干流是四大家魚(yú)主要的繁殖場(chǎng)所，作為典型的產(chǎn)漂流性卵魚(yú)類，其繁殖行為及早期生長(zhǎng)與長(zhǎng)江生態(tài)水文條件相適應(yīng). 水利工程的建設(shè)改變了長(zhǎng)江原有的生態(tài)環(huán)境，使得現(xiàn)有的滿足四大家魚(yú)繁殖需求的江段減少，同時(shí)受過(guò)度捕撈等因素干擾，家魚(yú)親魚(yú)群體破壞較為嚴(yán)重，致使長(zhǎng)江流域家魚(yú)數(shù)量逐漸減少[2]. 近年來(lái)長(zhǎng)江中游監(jiān)利江段和漢江丹江口水庫(kù)上游早期資源調(diào)查結(jié)果均顯示四大家魚(yú)資源量顯著下降，2007年監(jiān)利江段四大家魚(yú)卵苗徑流量?jī)H為1981年的1.34%，2014年丹江口水庫(kù)上游家魚(yú)卵苗徑流量下降到1977年的0.25%[3-4]. 長(zhǎng)江干流四大家魚(yú)產(chǎn)卵場(chǎng)多集中在中上游江段，歷史調(diào)查資料表明下游僅在湖口和彭澤江段發(fā)現(xiàn)四大家魚(yú)產(chǎn)卵場(chǎng)[5]. 任鵬的研究表明長(zhǎng)江下游由于江湖連通效應(yīng)及支流的匯入，距離大壩越遠(yuǎn)的江段其水文形勢(shì)越趨于自然狀態(tài)，使得四大家魚(yú)仔魚(yú)資源量自中游向下游安慶江段有增多的趨勢(shì)[6].

安慶江段沿岸有眾多附屬湖泊，分布有長(zhǎng)江江豚自然生態(tài)保護(hù)區(qū)、長(zhǎng)江刀鱭國(guó)家級(jí)水產(chǎn)種質(zhì)資源保護(hù)區(qū)等重要生態(tài)功能區(qū)，是重要的漁業(yè)資源水域[7]. 位于上游的湖口江段與鄱陽(yáng)湖形成特殊的江湖復(fù)合型生態(tài)系統(tǒng)，對(duì)四大家魚(yú)的繁殖及育幼提供良好的生境，使得距離鄱陽(yáng)湖較近的安慶江段家魚(yú)仔魚(yú)豐度高于南京和靖江江段[8]. 安慶市是作為長(zhǎng)江三角洲重要的水運(yùn)樞紐和化工城市，頻繁的航運(yùn)、化工污染和過(guò)度捕撈等因素導(dǎo)致安慶江段漁業(yè)資源日益衰退. 2015年春季安慶江段成魚(yú)調(diào)查中四大家魚(yú)僅有鰱和草魚(yú)為常見(jiàn)種，且鰱個(gè)體呈現(xiàn)小型化現(xiàn)象，四大家魚(yú)資源明顯衰退[9]. 本研究以長(zhǎng)江下游安慶段的家魚(yú)仔稚魚(yú)(初孵仔魚(yú)至鱗片形成期[10])為研究對(duì)象，在繁殖高峰期對(duì)安慶江段進(jìn)行持續(xù)監(jiān)測(cè)，了解安慶江段四大家魚(yú)仔稚魚(yú)資源量現(xiàn)狀，運(yùn)用GAM模型解析該江段家魚(yú)仔稚魚(yú)豐度變化與水文環(huán)境因子之間的關(guān)系，為評(píng)估長(zhǎng)江下游四大家魚(yú)資源補(bǔ)充量、苗汛預(yù)測(cè)和保護(hù)提供基礎(chǔ)數(shù)據(jù).

1 材料與方法

1.1 采樣時(shí)間和地點(diǎn)

通過(guò)查閱長(zhǎng)江魚(yú)類早期資源的歷史文獻(xiàn)[10]及現(xiàn)場(chǎng)勘察，2016-2018年4-8月(2016年：4月17日-8月15日；2017年：4-8月逐月采樣一周；2018年：4月22日-8月10日，共計(jì)239天)于長(zhǎng)江下游安慶段開(kāi)展早期資源調(diào)查. 采樣斷面設(shè)置于長(zhǎng)江干流與皖河口交匯處上游50 m處(圖1)，斷面設(shè)置北岸(30°29′5″N， 116°59′24″E)、南岸(30°28′48″N， 116°59′59″E)和江心(30°28′56″N， 116°59′40″E)3個(gè)采樣點(diǎn)，南北兩岸采樣點(diǎn)距離岸邊10 m左右.

圖1 安慶江段魚(yú)類早期資源采樣點(diǎn)設(shè)置

1.2 調(diào)查方法和樣品的保存鑒定

四大家魚(yú)早期資源采集使用大型圓錐狀浮游生物網(wǎng)進(jìn)行主動(dòng)采集. 網(wǎng)具為直徑0.8 m、網(wǎng)深2.5 m、網(wǎng)目40目、孔徑0.25 mm的圓錐浮游生物網(wǎng)，網(wǎng)后連接圓柱形魚(yú)苗收集裝置，網(wǎng)口固定流量計(jì)用于測(cè)量網(wǎng)口過(guò)水量，采樣時(shí)使用機(jī)動(dòng)船對(duì)江水上層進(jìn)行拖拽. 調(diào)查期間，每日8：00-11：00之間對(duì)各采樣點(diǎn)進(jìn)行定量采集一次，根據(jù)采集到的魚(yú)類早期資源總捕撈數(shù)確定每次采集時(shí)間(3～10 min)，采樣同時(shí)記錄每次拖網(wǎng)持續(xù)時(shí)間. 使用Dretec電子水溫計(jì)測(cè)量水溫，采用薩氏透明度盤(pán)測(cè)量透明度，使用德國(guó)HYDRO-BIOS公司數(shù)字網(wǎng)口流量計(jì)記錄流量計(jì)轉(zhuǎn)數(shù)來(lái)?yè)Q算網(wǎng)口過(guò)水流量. 采樣斷面的水位和江水徑流量數(shù)據(jù)來(lái)分別自安徽省水文信息網(wǎng)的安慶水文站和大通水文站(由于采樣斷面和銅陵市大通鎮(zhèn)之間無(wú)大型支流匯入長(zhǎng)江干流，故可以反映該斷面的江水徑流量情況).

主動(dòng)網(wǎng)具采集的家魚(yú)仔稚魚(yú)樣品由冰水保存，樣品鑒定完成后使用75%中性酒精保存. 仔稚魚(yú)種類鑒定按照曹文宣等[10]編寫(xiě)的《長(zhǎng)江魚(yú)類早期資源》，通過(guò)形態(tài)學(xué)方法在奧林巴斯解剖鏡SZX16進(jìn)行種類鑒定，并分類計(jì)數(shù). 對(duì)于形態(tài)破損難以鑒定確認(rèn)的物種使用分子方法鑒定，具體步驟是：提取樣本的基因組DNA，基于COⅠ(Cytochorme C oxidase Subunit Ⅰ，線粒體細(xì)胞色素C氧化酶亞基Ⅰ)基因條形碼，采用PCR擴(kuò)增目的條帶并測(cè)序，將所測(cè)結(jié)果與已經(jīng)用成魚(yú)建立的長(zhǎng)江魚(yú)類條形碼數(shù)據(jù)庫(kù)進(jìn)行比對(duì)得到種類鑒定的結(jié)果[11].

1.3 數(shù)據(jù)處理與分析

1)仔稚魚(yú)密度計(jì)算方法參考易伯魯?shù)萚12]的計(jì)算方法，依據(jù)所采集四大家魚(yú)仔稚魚(yú)的數(shù)量、采集時(shí)間、流量計(jì)的始末差值和主動(dòng)網(wǎng)口面積按下列公式處理：

qi=0.3Ci·a

(1)

Di=Ni/qit

(2)

式中，qi為第i次采集網(wǎng)具網(wǎng)口的過(guò)水流量(m3/s)，Ci為第i次采集流量計(jì)的流量差，a為主動(dòng)網(wǎng)具網(wǎng)口面積(m2)，0.3為流量計(jì)的常數(shù)，t為每次采集仔稚魚(yú)時(shí)間(s)，Di為第i次采集仔稚魚(yú)的密度(ind./m3)，Ni為第i次采集仔稚魚(yú)數(shù)量(ind.).

2)采樣斷面魚(yú)苗平均密度與定點(diǎn)魚(yú)苗密度相比系數(shù)：

(3)

3) 采集期間卵苗徑流量：

Mi=Di·Qi·C

(4)

式中，Mi為第i次采集時(shí)段內(nèi)通過(guò)該江段斷面的魚(yú)苗數(shù)(ind.)，Qi為第i次采集時(shí)的斷面流量(m3/s).

非采集期間魚(yú)苗徑流量：

Mi,j+1=(Mi/ti+Mi+1/ti+1)ti,j+1/2

(5)

用相鄰兩次采集的徑流量及其間隔時(shí)間進(jìn)行插補(bǔ)計(jì)算. 式中：Mi，i+1為第i，i+1采集時(shí)間間隔內(nèi)的魚(yú)苗徑流量(ind.)，ti，i+1為第i，i+1次采集時(shí)間間隔(s).

采集江段的魚(yú)苗總徑流量為：

M=∑Mi+∑Mi,j+1

(6)

依據(jù)劉建康對(duì)長(zhǎng)江下游四大家魚(yú)魚(yú)苗垂直分布的研究結(jié)果，江水上層家魚(yú)苗平均捕獲率占總捕獲率的43.5%[13]，估算本研究中繁殖期內(nèi)通過(guò)安慶江段的四大家魚(yú)仔稚魚(yú)總徑流量.

4)廣義可加模型(generalized additive models，GAM)是一種非參數(shù)化或半?yún)?shù)化的回歸分析方法，能很好地模擬應(yīng)變量和一個(gè)或多個(gè)預(yù)測(cè)變量之間的關(guān)系[14]. 本研究運(yùn)用GAM模型分析水文因子對(duì)家魚(yú)仔稚魚(yú)豐度變化的因子進(jìn)行分析，其表達(dá)式為：

(7)

式中，g(ρ)為關(guān)聯(lián)函數(shù)；ρ為仔稚魚(yú)豐度(ind./m3)；α為適合函數(shù)中的截距；xj為解釋變量，包括水溫、透明度、水位上漲率和徑流量上漲率；fi(xj)為解釋變量關(guān)系的非參數(shù)函數(shù)；ε為誤差項(xiàng)，與解釋變量xj無(wú)關(guān)，E(ε)=0,ε=δ2模型采用樣條平滑法，對(duì)數(shù)據(jù)圖進(jìn)行平滑處理.

數(shù)據(jù)處理使用Excel 2016軟件，使用Origin 2018軟件制作直方圖等分析圖. 數(shù)據(jù)分析使用SPSS 19.0軟件，差異性分析使用單因素方差分析，相關(guān)性分析使用Pearson相關(guān)性分析. GAM模型的構(gòu)建及其檢驗(yàn)在R軟件(R 3.5.3)的mgcv包內(nèi)完成.

2 結(jié)果與分析

2.1 安慶江段四大家魚(yú)仔稚魚(yú)種類組成及徑流量

2016-2018年共調(diào)查239天，采集家魚(yú)仔稚魚(yú)3178尾，其中以鰱數(shù)量最多，共2160尾，占家魚(yú)仔稚魚(yú)總捕撈數(shù)量的67.97%；其次是草魚(yú)，共計(jì)732尾占23.03%；鳙279尾占8.78%；青魚(yú)數(shù)量最少，僅采集到7尾，所占比例為0.22%. 2016-2018年四大家魚(yú)仔稚魚(yú)變化規(guī)律呈現(xiàn)出鰱所占比例逐年上升，草魚(yú)比例不斷下降，鳙所占比例較為穩(wěn)定(表1). 根據(jù)采集到的四大家魚(yú)仔稚魚(yú)數(shù)量及斷面流量估算，2016-2018年4-8月份通過(guò)安慶江段四大家魚(yú)仔稚魚(yú)徑流量依次為21.70×108、14.62×108、12.05×108尾，占各年份仔稚魚(yú)總徑流量的2.43%、0.99%和0.65%.

表1 安慶江段四大家魚(yú)仔稚魚(yú)種類組成

2.2 四大家魚(yú)仔稚魚(yú)豐度的時(shí)空分布特征

2016年安慶江段家魚(yú)仔稚魚(yú)日平均豐度為8.57±9.03 ind./100 m3. 仔稚魚(yú)自4月中旬呈波動(dòng)性上升趨勢(shì)至7月中旬達(dá)到峰值，集中出現(xiàn)在5月下旬至7月中旬，共出現(xiàn)3次苗汛，分別為5月21日、6月15日、7月4日，對(duì)應(yīng)的仔稚魚(yú)密度為26.50、37.83、43.23 ind./100 m3. 2017年家魚(yú)日平均豐度為6.76±4.29 ind./100 m3，采集時(shí)間內(nèi)最高豐度為13.77 ind./100 m3. 對(duì)應(yīng)日期為7月8日. 2018年日平均豐度為7.49±12.52 ind./100 m3，仔稚魚(yú)豐度隨江段徑流量變化呈現(xiàn)先升后降再上升的趨勢(shì)，在5月27日和7月18日達(dá)到峰值，對(duì)應(yīng)的密度分別為13.74和55.86 ind./100 m3(圖2).

對(duì)斷面北岸、南岸和江心3個(gè)采樣點(diǎn)的家魚(yú)密度數(shù)據(jù)進(jìn)行顯著性分析，驗(yàn)證家魚(yú)密度隨斷面橫向分布的差異性. Kruskal-Wallis檢驗(yàn)表明，斷面3個(gè)采樣點(diǎn)的家魚(yú)仔稚魚(yú)漂流密度均存在顯著性差異(2016年：P=0.026<0.05，n=57；2017年：P=0.041<0.05，n=5；2018年：P=0.046<0.05，n=49). 2016-2018年斷面家魚(yú)密度橫向分布均呈現(xiàn)出北岸>南岸>江心，三年內(nèi)北岸平均漂流密度依次為14.14、17.83、16.30 ind./100 m3；南岸平均漂流密度依次為7.63、4.79、3.66 ind./100 m3；江心平均漂流密度依次為3.97、3.81、2.53 ind./100 m3.

圖2 安慶江段四大家魚(yú)仔稚魚(yú)豐度變化

2.3 四大家魚(yú)仔稚魚(yú)豐度與水文環(huán)境因子之間的關(guān)系

調(diào)查期間內(nèi)，安慶江段水溫變化范圍為18.0～29.8℃，平均值為25.30±2.44℃；透明度為16～70 cm，平均值為46.78±15.64 cm；溶解氧濃度為6.17～8.07 mg/L，平均值為7.24±2.15 mg/L；pH 7.39～8.17，平均值為7.80±0.20；江水徑流量20400～70200 m3/s，日上漲率在100～7000 m3/(s·d)；水位8.44～17.69 m，日上漲率在0.01～1.25 m/d之間. 通過(guò)對(duì)家魚(yú)仔稚魚(yú)豐度與水文環(huán)境因子進(jìn)行的相關(guān)性分析發(fā)現(xiàn)，家魚(yú)日豐度變化與水溫呈顯著正相關(guān)(P<0.05)，與透明度呈極顯著負(fù)相關(guān)(P<0.01)，與江水徑流量、水位均呈極顯著正相關(guān)(P<0.01).

GAM模型構(gòu)建過(guò)程中逐步加入各個(gè)水文因子，依據(jù)AIC準(zhǔn)則來(lái)選擇最優(yōu)模型. 模型最終結(jié)果見(jiàn)表2，其總偏差解釋率為69.8%，預(yù)測(cè)模型的校正決定系數(shù)R-sq為0.622>0.5，說(shuō)明模型較高程度上揭示了水文因子對(duì)家魚(yú)仔稚魚(yú)豐度的影響.

應(yīng)用正態(tài)Q-Q圖和Shapiro-wilk檢驗(yàn)驗(yàn)證模型殘差的正態(tài)性，診斷GAM模型. 正態(tài)Q-Q圖表明GAM模型基本符合正態(tài)分布(圖3). Shapiro-wilk檢驗(yàn)表明，模型殘差符合正態(tài)分布(w=0.941，P=0.070>0.05).

表2 安慶江段四大家魚(yú)仔稚魚(yú)GAM模型方差分析

圖3 四大家魚(yú)仔稚魚(yú)豐度的GAM模型Q-Q圖

圖4 四大家魚(yú)仔稚魚(yú)豐度與各因子之間關(guān)系的GAM分析

GAM模型結(jié)果顯示，根據(jù)偏差解釋率，各因子相對(duì)重要性依次為水位日上漲率>透明度>水溫>徑流量日上漲率. 水文因子與家魚(yú)仔稚魚(yú)豐度之間關(guān)系的分析結(jié)果圖表明(圖4)，當(dāng)水溫在19.0～27.3℃時(shí)，家魚(yú)仔稚魚(yú)大量出現(xiàn)，其豐度隨水溫升高整體呈現(xiàn)上升趨勢(shì)；在超過(guò)27.3℃后，其豐度逐漸下降. 透明度范圍主要集中在24～57 cm之間，家魚(yú)仔稚魚(yú)密度與透明度呈現(xiàn)顯著負(fù)相關(guān)性，即江水透明度下降時(shí)家魚(yú)仔稚魚(yú)密度上升. 水位和徑流量的日上漲率與家魚(yú)仔稚魚(yú)豐度均呈現(xiàn)極顯著正相關(guān)性，說(shuō)明家魚(yú)豐度隨水位和徑流量的日上漲率增大有所上升. 其中家魚(yú)高峰期對(duì)應(yīng)的水溫、水位日上漲率和徑流量日上漲率分別為27.3℃、0.23 m/d和900 m3/(s·d).

3 討論

3.1 安慶江段四大家魚(yú)仔稚魚(yú)資源現(xiàn)狀

三峽大壩運(yùn)行以來(lái)，長(zhǎng)江中上游四大家魚(yú)卵苗種類組成出現(xiàn)以濾食浮游生物的鰱魚(yú)苗數(shù)量迅速增加，草魚(yú)所占比例不斷下降，鳙和青魚(yú)處于較低水平[2，15]，本研究結(jié)果與中上游相類似. 長(zhǎng)江下游安慶江段2016-2018年四大家魚(yú)組成中鰱以67.97%的比例占絕對(duì)優(yōu)勢(shì)，且所占比例有一定上升趨勢(shì)；其次是草魚(yú)，呈現(xiàn)出逐年下降的趨勢(shì)；鳙和青魚(yú)所占比例較低. 與2013年[6]結(jié)果相比青魚(yú)資源量占比大幅度下降，可能原因是由于長(zhǎng)江下游懸移質(zhì)粒徑較小，泥沙和有機(jī)質(zhì)等營(yíng)養(yǎng)物質(zhì)在水流搬運(yùn)過(guò)程中不易沉積，不利于大型底棲動(dòng)物的生存和建群[16]. 作為底棲生物食性的青魚(yú)，餌料的匱乏直接影響其生長(zhǎng)及育幼過(guò)程中正常的攝食活動(dòng)，使得青魚(yú)早期資源量得不到有效補(bǔ)充. 長(zhǎng)江干流中青魚(yú)資源量的減少應(yīng)引起重視.

由于沿江支流的匯入及水流作用，距離大壩較遠(yuǎn)的江段其水文情勢(shì)越接近自然狀態(tài)，大壩對(duì)魚(yú)類資源的不利影響就會(huì)越小[17]. 2016年安慶江段四大家魚(yú)仔稚魚(yú)日平均密度為8.57 ind./100 m3，顯著高于同期中游黃石江段的卵苗日均密度[18]，表明四大家魚(yú)資源量自長(zhǎng)江中游部分江段向下游安慶江段有增多的趨勢(shì). 長(zhǎng)江下游距離三峽大壩較遠(yuǎn)，水文特征受影響較小，同時(shí)鄱陽(yáng)湖與長(zhǎng)江的連通關(guān)系不僅有利于家魚(yú)親魚(yú)出湖繁殖，也有利于仔魚(yú)入湖育幼，鄱陽(yáng)湖和長(zhǎng)江形成的江湖復(fù)合生態(tài)系統(tǒng)可能是安慶江段家魚(yú)資源量較大的原因.

近年來(lái)對(duì)長(zhǎng)江流域早期資源的調(diào)查結(jié)果均呈現(xiàn)出四大家魚(yú)資源量處于不斷衰退的趨勢(shì). 2014-2015年劉明典等對(duì)長(zhǎng)江中游早期資源調(diào)查結(jié)果顯示大壩下泄的不飽和水流改變了原有的水文環(huán)境，使得宜昌江段四大家魚(yú)產(chǎn)卵場(chǎng)規(guī)模呈減少趨勢(shì)[19]. 受大量涉水工程和圍湖造田的影響，安慶江段和沿江湖泊的連通關(guān)系均發(fā)生變化，導(dǎo)致江湖洄游性魚(yú)類自然繁殖受到影響. 本研究中估算2016-2018年通過(guò)安慶江段的四大家魚(yú)仔稚魚(yú)徑流量分別為21.70×108尾、14.62×108尾、12.05×108尾，呈現(xiàn)出逐年下降的趨勢(shì). 受江湖阻隔的影響，安慶江段魚(yú)類群落結(jié)構(gòu)表現(xiàn)出洄游性魚(yú)類減少的趨勢(shì)[20]，長(zhǎng)江干流四大家魚(yú)繁殖親本數(shù)量的減少，使得仔稚魚(yú)徑流量逐年下降. 此外，作為典型的產(chǎn)漂流性卵魚(yú)類，四大家魚(yú)產(chǎn)卵規(guī)模同江水漲水過(guò)程密切相關(guān)[21]. 2018年江段徑流量較往年相比上漲持續(xù)時(shí)間較短，上漲幅度較低. 漲水過(guò)程的縮減一定程度上影響四大家魚(yú)的正常繁殖活動(dòng).

3.2 影響四大家魚(yú)仔稚魚(yú)時(shí)空分布差異的因素

魚(yú)類早期資源的時(shí)間動(dòng)態(tài)能直接反映魚(yú)類繁殖時(shí)間的變化特征[22]. 安慶江段四大家魚(yú)仔稚魚(yú)集中出現(xiàn)在5-7月，高峰期多發(fā)生在6月下旬和7月，與長(zhǎng)江中游四大家魚(yú)繁殖時(shí)間大致相同. 2016-2018年安慶江段四大家魚(yú)苗高峰期豐度有所上升，但其主汛期存在推遲且持續(xù)時(shí)間有縮短的跡象. 2016和2017年主汛期均處于6月下旬至7月上旬，2018年則推遲至7月中下旬. 可能是由于2018年5月份江水漲幅較低，不滿足家魚(yú)繁殖所需的水文條件，使得四大家魚(yú)集中在7月中下旬江水上漲過(guò)程時(shí)進(jìn)行繁殖活動(dòng).

河流魚(yú)類早期階段受非生物和生物因素影響會(huì)存在一定的漂流規(guī)律. 非生物因子包括洪水的水位漲幅變化、河床形態(tài)、光照強(qiáng)度、水溫水文條件等，生物因子有魚(yú)類的發(fā)育階段、棲息地偏好及其種間關(guān)系等[23]. 徐薇等研究發(fā)現(xiàn)長(zhǎng)江魚(yú)類在早期階段存在水平分布上的顯著性差異，斷面空間卵苗密度與流速存在負(fù)相關(guān)性關(guān)系[24]. 本文在分析采樣斷面四大家魚(yú)仔稚魚(yú)豐度空間分布時(shí)，發(fā)現(xiàn)江心豐度顯著低于南北兩岸. 江心水流量大，流速較兩岸快，不利于仔稚魚(yú)的棲息和生存，使得具有一定游泳和主動(dòng)選擇能力的仔魚(yú)聚集在沿岸緩流回水區(qū). 此外其分布還受生物等其他因素的影響，采樣斷面3個(gè)點(diǎn)位的年際變化呈現(xiàn)南岸豐度均值逐年接近江心，北岸家魚(yú)仔稚魚(yú)豐度有一定的增長(zhǎng)趨勢(shì). 采樣斷面屬于順直型江段，兩岸江水流速無(wú)較大差異，其南北兩岸家魚(yú)仔稚魚(yú)豐度可能受生境差異等因素影響. 北岸擁有較高的植被覆蓋率，相比南岸航運(yùn)及農(nóng)耕等人類活動(dòng)所造成的污染可以提供更為豐富的生物餌料，有利于仔稚魚(yú)早期的攝食活動(dòng).

3.3 水文環(huán)境因子對(duì)四大家魚(yú)仔稚魚(yú)豐度的影響

四大家魚(yú)繁殖特性相似，其產(chǎn)卵所需的水文條件也相似，不同種類在相同產(chǎn)卵場(chǎng)進(jìn)行繁殖活動(dòng)，常被視作一個(gè)類群進(jìn)行研究. 作為典型的產(chǎn)漂流性卵魚(yú)類，親魚(yú)需要在特定的水溫條件下，同時(shí)有適宜水位漲幅才能促使排卵[25]. 水溫被認(rèn)為是影響魚(yú)類產(chǎn)卵的重要因素，適宜四大家魚(yú)產(chǎn)卵的水溫一般在18～24℃. 長(zhǎng)江中游監(jiān)利斷面家魚(yú)卵苗高峰期在6月下旬至7月初，水溫范圍在20～24℃[26]. 本研究中安慶江段同期水溫在22.4～25.2℃，略高于監(jiān)利江段. 推測(cè)水溫的差異是導(dǎo)致家魚(yú)資源量長(zhǎng)江中游向下游有增多趨勢(shì)的原因. GAM模型結(jié)果中安慶江段水溫范圍在19～27.3℃之間，隨著水溫的上升，家魚(yú)仔稚魚(yú)豐度整體呈上升趨勢(shì).

透明度隨江水徑流量的變化被認(rèn)為是魚(yú)類產(chǎn)卵起始的重要因素. 張曉敏等研究2004-2007年漢江中下游洪峰過(guò)程中發(fā)現(xiàn)，當(dāng)水溫和漲水條件均滿足家魚(yú)繁殖需求的情況下，僅在透明度低于14.5 cm時(shí)監(jiān)測(cè)到家魚(yú)卵[27]. 因此洪水過(guò)程中透明度的高低和家魚(yú)自然繁殖存在一定的相關(guān)性. 本文結(jié)果表明家魚(yú)仔稚魚(yú)密度與透明度呈顯著負(fù)相關(guān)關(guān)系，隨著透明度下降，家魚(yú)密度升高. 即在繁殖季節(jié)，透明度的降低提高了家魚(yú)產(chǎn)卵的可能性，同時(shí)有利于仔稚魚(yú)躲避敵害的能力.

此外，水流刺激也是四大家魚(yú)產(chǎn)卵的必備條件. 易雨君等探究長(zhǎng)江中游家魚(yú)棲息地適宜度曲線的結(jié)果表明水位的變幅是家魚(yú)繁殖活動(dòng)的必備條件，魚(yú)苗徑流量與洪峰呈正相關(guān)關(guān)系，長(zhǎng)江中游家魚(yú)產(chǎn)卵的理想水位日上漲率為0.3 m/d[28]. 本研究發(fā)現(xiàn)安慶江段四大家魚(yú)仔稚魚(yú)豐度伴隨著江水水位和徑流量上漲有顯著的上升趨勢(shì)，證明水位和徑流量的上漲是刺激四大家魚(yú)繁殖的重要因子. GAM模型結(jié)果顯示安慶江段水位日上漲率達(dá)到0.23 m/d或徑流量日上漲率達(dá)到900 m3/(s·d)時(shí)，仔稚魚(yú)豐度開(kāi)始顯著上升.

3.4 對(duì)長(zhǎng)江下游四大家魚(yú)資源保護(hù)的建議

長(zhǎng)江下游城市經(jīng)濟(jì)發(fā)展迅速，人類活動(dòng)頻繁，為保護(hù)和恢復(fù)四大家魚(yú)資源，基于以往和本文研究結(jié)果建議：(1)堅(jiān)持貫徹落實(shí)長(zhǎng)江流域禁漁制度，保證家魚(yú)親本的正常繁殖. (2)動(dòng)態(tài)監(jiān)測(cè)魚(yú)類資源，以便制定和調(diào)整增殖保護(hù)對(duì)策. 針對(duì)下游鳙和青魚(yú)數(shù)量稀少的情況，應(yīng)通過(guò)開(kāi)展科學(xué)增殖放流，增加繁殖群體數(shù)量，促進(jìn)漁業(yè)資源恢復(fù). (3)應(yīng)注意保護(hù)洲灘及岸線的自然屬性，尤其是實(shí)施航道整治和護(hù)岸等工程時(shí). (4)嚴(yán)控各類污染物的排放，保證魚(yú)類棲息地環(huán)境質(zhì)量. (5)確保通江湖泊與長(zhǎng)江干流的連通性，并逐步恢復(fù)其他沿江湖泊及支流與長(zhǎng)江干流的連通性，充分發(fā)揮江湖復(fù)合生態(tài)系統(tǒng)的作用，以此來(lái)保護(hù)和恢復(fù)長(zhǎng)江下游家魚(yú)資源.

4 小結(jié)與展望

本研究于2016-2018年對(duì)安慶江段四大家魚(yú)仔稚魚(yú)資源進(jìn)行了連續(xù)監(jiān)測(cè)，掌握了四大家魚(yú)仔稚魚(yú)資源分布情況和苗汛規(guī)律，并利用GAM模型初步揭示了水文水環(huán)境條件對(duì)家魚(yú)早期資源豐度的影響，為安慶段四大家魚(yú)資源的保護(hù)提供了數(shù)據(jù)支撐.

1)2016-2018年長(zhǎng)江安慶江段四大家魚(yú)仔稚魚(yú)種類主要是鰱和草魚(yú)；其次是鳙；青魚(yú)數(shù)量極少僅占家魚(yú)總數(shù)的0.22%. 估算3年內(nèi)通過(guò)該江段的家魚(yú)仔稚魚(yú)徑流量分別為21.70×108、14.62×108、12.05×108尾，呈逐年下降的趨勢(shì).

2)安慶江段四大家魚(yú)仔稚魚(yú)自4月中旬整體豐度呈波動(dòng)性上升趨勢(shì)，在7月中上旬達(dá)到高峰期，其主汛期有明顯推遲的現(xiàn)象. 斷面空間分布規(guī)律為北岸>南岸>江心，且北岸豐度有一定的增長(zhǎng).

3)在繁殖季節(jié)中，四大家魚(yú)仔稚魚(yú)豐度與水溫、透明度、水位日上漲率和徑流量日上漲率呈顯著相關(guān)，其中水位日上漲率影響較為顯著. 家魚(yú)仔稚魚(yú)豐度與透明度呈負(fù)相關(guān)，與水溫(24.0～27.3℃)、水位日上漲率和徑流量日上漲率均呈正相關(guān)關(guān)系.

4)鑒于仔稚魚(yú)資源的垂直空間分布對(duì)總徑流量估算的影響，后續(xù)工作將開(kāi)展仔稚魚(yú)的垂直分布情況研究，提高徑流量估算的準(zhǔn)確度.

5)利用仔稚魚(yú)耳石鑒定家魚(yú)日齡，分析其早期生長(zhǎng)特征，結(jié)合該江段水文數(shù)據(jù)推算安慶江段以上江段四大家魚(yú)產(chǎn)卵場(chǎng)位置及規(guī)模，為四大家魚(yú)保護(hù)提供更準(zhǔn)確的資料.