摘要：目前在樂安河流域內(nèi)展開的河流重金屬的污染研究亟待發(fā)掘。為充分闡明樂安河河源區(qū)水庫環(huán)境特征及其環(huán)境變化驅(qū)動過程，以婺源縣境內(nèi)不同生態(tài)環(huán)境與開發(fā)背景的清華和石井水庫為研究對象。通過實驗記錄沉積柱巖芯中硅藻種類及理化指標，并結合現(xiàn)代監(jiān)測資料及多元統(tǒng)計方法，分析影響水庫中沉積硅藻組合群落變化及沉積環(huán)境變化的環(huán)境特征驅(qū)動因子。結果表明，水庫營養(yǎng)水平隨深度至表層逐漸升高，沉積物污染評價也揭示了營養(yǎng)物質(zhì)的匯集過程。沉積硅藻物種耐污種和富營養(yǎng)化指示種取代由普生種和清潔指示種成為優(yōu)勢，代表性物種為Achnanthidium minutissima、Navicula minima。多元統(tǒng)計分析結果表明，TP、TN及TP、TOC分別為清華和石井水庫硅藻組合群落結構分布的主要驅(qū)動因子（p＜0.05）。18項硅藻指數(shù)的箱線圖分析結果顯示，CEE、LOBO指數(shù)被篩選并表現(xiàn)出良好的評價適宜性及相關性。不同人類活動方式影響庫區(qū)內(nèi)輸入的外源物質(zhì)，從而驅(qū)動硅藻群落的演替，2座水庫皆有向富營養(yǎng)化邁進的趨勢。

關鍵詞：硅藻；沉積物；富營養(yǎng)化；指數(shù)；樂安河

中圖分類號：TV21 文獻標識碼：A 文章編號：1001-9235（2024）04-0077-10

Distribution Characteristics and Environmental Response of Diatoms Deposited in Reservoirs in the Source Region of Le’an River

YIN Hongwei，ZHANG Jie^*，YAN Ruyu，JI Yong，XIAO Weiqi，CAI Gaotang

（School of Hydraulic and Ecological Engineering，Nanchang Institute of Technology，Nanchang 330099，China）

Abstract： At present，the study of heavy metal pollution in the Le’an River Basin should be explored.To fully elucidate the environmental characteristics and change-driven processes of reservoirs in the source region of the Le’an River，this paper employs Qinghua and Shijing reservoirs with different ecological environments and development backgrounds as research objects.Experimental records of diatom species and physicochemical indicators in sedimentary column cores are conducted，and meanwhile combined with modern monitoring data and multivariate statistical methods，environmental characteristics driving factors affecting changes in sedimentary diatom assemblages and sedimentary environment in reservoirs are analyzed.The results show that nutrient levels in the reservoirs gradually increase with the depth to the surface，and the sediment pollution evaluation also reveals the collection process of nutrients.Universal and cleaning indicator species have been replaced by fouling-tolerant and eutrophication indicator species of sedimentary diatoms，with representative species being Achnanthidium minutissima and Navicula minima.The multivariate statistical analysis results show that TP，TN，TP，and TOC are the main driving factors of diatom community structure distribution in Qinghua and Shijing reservoirs respectively （plt;0.05）.Box plot analysis results of 18 diatom indices reveal that CEE and LOBO indices are screened and show sound evaluation suitability and correlation.The input of exogenous materials into the reservoirs is affected by different anthropogenic practices to drive the succession of the diatom community，with a trend towards eutrophication in both reservoirs.

Keywords：diatom;sediment;eutrophication;index;Le’an River

在經(jīng)濟快速發(fā)展條件下，水庫正面臨環(huán)境質(zhì)量和生態(tài)系統(tǒng)健康帶來挑戰(zhàn)^［1］。其中，人類活動造成的營養(yǎng)物質(zhì)輸入，會使得水域環(huán)境在短時間內(nèi)急劇惡化，加速水體富營養(yǎng)化過程。硅藻作為指示物種主要作用為評估水域生態(tài)環(huán)境的一個有力手段^［2-3］。硅藻分布廣泛、分類系統(tǒng)成熟、數(shù)量龐大，同時對環(huán)境微小變化反應明顯，通常用來解釋水域營養(yǎng)狀態(tài)的變化^［4］。硅藻死亡后，其遺骸大部分能夠長期完整保留下來。硅藻作為指示生物，國際上目前研究主要在以下幾個領域：①構建適合于特定水域的監(jiān)測與評價指數(shù)體系^［5］；②建立近代流域環(huán)境變化過程^［6］；③反演流域古環(huán)境變化過程^［7］；④評估與預測氣候變化過程^［8］。因此，在確保采樣點沉積記錄的代表性、連續(xù)性和完整性的前提下，應用硅藻群落分析水庫生態(tài)系統(tǒng)短期的發(fā)展和演變歷程，是一種新型和極具潛力的研究方法。

樂安河為中國第一淡水湖鄱陽湖支流饒河的上游主干流，受上游分支德興銅礦、銀山鉛鋅礦和冶煉廠的影響，樂安河重金屬污染問題已受到廣泛關注并在區(qū)域中開展了諸多研究^［9-12］。前期眾多學者主要聚焦于樂安河流域重金屬分布及其影響因素和風險評價，未能充分了解調(diào)查正源區(qū)所處環(huán)境狀況，對流經(jīng)婺源的正源區(qū)域進行硅藻群落分析及環(huán)境調(diào)查較少^［13-15］。樂安河正源所處婺源縣具有良好的生態(tài)環(huán)境，區(qū)域水資源豐富，水庫眾多，是地區(qū)主要水資源也是流域內(nèi)重要的航道，其環(huán)境質(zhì)量也將影響下游輸入鄱陽湖的水環(huán)境質(zhì)量^［16］。為了解短時間尺度上樂安河源區(qū)水庫環(huán)境狀況及沉積物污染情況，本文通過研究河源區(qū)分布于城鎮(zhèn)周邊的石井水庫和河道型清華水庫沉積硅藻群落結構特征結合沉積物理化指標，探求河源區(qū)水庫沉積環(huán)境特征及關鍵驅(qū)動因素，以期為保護樂安河流域水質(zhì)安全及鄱陽湖生態(tài)系統(tǒng)健康提供前期數(shù)據(jù)。

樂安河婺源以上流域面積2 404.4 km²，生態(tài)系統(tǒng)完整，地貌復雜多樣，屬饒河一級支流，具備從源頭至河口不同區(qū)段的河流特征，同時各種人類活動（水利工程、重金屬污染、湖庫富營養(yǎng)化及各類土地開發(fā)利用等）干擾顯著，目前在樂安河流域內(nèi)展開的河流主要針對流經(jīng)德興的上游分支的重金屬污染研究，針對較少污染的樂安河上游婺源段所開展的生態(tài)研究亟待發(fā)掘。根據(jù)不同區(qū)段自然環(huán)境特征及人類干擾行為，本研究開展區(qū)域自然生態(tài)環(huán)境調(diào)查、硅藻種群分布特征調(diào)查、指示硅藻對不同脅迫因子的指示效應。通過柱狀沉積物樣本分析、底棲硅藻監(jiān)測、水環(huán)境與水生態(tài)分析，獲得研究區(qū)域完整的硅藻、水文與水生態(tài)歷史過程數(shù)據(jù)。利用對沉積物物理、化學特征的分析，可以進一步識別不同的環(huán)境因子的生態(tài)影響程度，分析沉積物理化指標對環(huán)境因素的響應模式差異，建立有效的硅藻與不同環(huán)境因子之間的統(tǒng)計學關系，闡明不同指標的響應模式差異和可能的有效指標組合，研究結果具有較好的理論價值與區(qū)域應用前景，為探明樂安河上游區(qū)生態(tài)環(huán)境具有一定的意義。

1 研究區(qū)域概況

樂安河正源段莘水發(fā)源于婺源縣北部大庾山、五龍山南麓。流域?qū)僦衼啛釒駶櫦撅L區(qū)，氣候溫和，雨量充沛。流域內(nèi)多年平均氣溫 16.7℃，多年平均蒸發(fā)量和降水量為1 342.5、1 845.19 mm。該區(qū)域?qū)儆诮系湫偷那鹆晟絽^(qū)，地勢東南高，西北低，流域內(nèi)林木蔥郁、峰巒疊嶂。本研究區(qū)域位于江西省上饒市婺源縣境內(nèi)清華水上清華水庫以及石井水庫。據(jù)前期調(diào)研以及當?shù)貙嶋H情況，研究聚焦于樂安河上游支流集發(fā)電、灌溉、防洪和養(yǎng)殖等綜合利用的水庫，分別選取樂安河支流高砂水上中云鎮(zhèn)曉林村的石井水庫，支流清華水上清華鎮(zhèn)的清華水庫，主支段莘水上的段莘水庫以及支流賦春水上的鴛鴦湖景區(qū)的大塘塢水庫，這些水庫代表著樂安河上游區(qū)域。在調(diào)研過程中發(fā)現(xiàn)其中石井水庫、段莘水庫存在網(wǎng)箱養(yǎng)殖的現(xiàn)象，清華水庫主要流經(jīng)城鎮(zhèn)區(qū)域，大塘塢水庫為鴛鴦湖景區(qū)主要發(fā)展旅游業(yè)?？紤]現(xiàn)場實際采樣情況，故安排石井水庫和大塘塢水庫在庫尾區(qū)設置采樣，清華水庫和段莘水庫在庫中進行采樣。2座水庫均具備不同環(huán)境特征，清華水庫為流經(jīng)清華鎮(zhèn)的河道型水庫，受城鎮(zhèn)居民生活干擾；石井水庫為婺源縣城周邊封閉型水庫，存在網(wǎng)箱養(yǎng)殖現(xiàn)象（圖1）。

2 材料與方法

2.1 樣品采集

清華水庫及石井水庫分別位于婺源縣城上下游區(qū)域，其受到不同生態(tài)特征及人類活動類型影響。2020年12月利用柱狀采樣器采集了25 cm（清華水庫）、43 cm（石井水庫）的沉積柱狀樣品，利用1 cm分樣器經(jīng)現(xiàn)場分層后立即裝入塑封袋標號記錄，共獲取有效樣品Q1—Q25和S1—S37，柱狀樣采集點詳見圖1。并使用GPS（佳明eTrex-30X）記錄采樣點區(qū)域的海拔與經(jīng)緯度信息。

2.2 化學分析

沉積物樣品自然風干后去除石子樹葉等雜質(zhì)后經(jīng)100 mm篩分后以備后續(xù)測定TN、TP、有機碳（TOC）等指標。其中，沉積物中TOC采用重鉻酸鉀氧化-分光光度法^［17］，TN采用凱氏法^［18］利用石墨消解儀（海能SH220F）消解樣品后利用凱氏定氮儀（海能K9840）進行蒸餾滴定測得，依據(jù)鉬銻抗比色法^［19］使用紫外分光光度計在700 nm波長環(huán)境下測定TP濃度。

2.3 硅藻鑒定

利用標準化硅藻處理方法^［20-21］處理樣品中沉積硅藻，稱取0.5 g過篩沉積樣，經(jīng)80℃水浴去除碳酸鹽及有機質(zhì)等雜質(zhì)后用水洗至中性，使用中性樹膠制成硅藻永久片。硅藻鏡檢使用奧林巴斯顯微鏡（型號CX31）在1 000倍油鏡下進行鑒定和計數(shù)，視野內(nèi)超過1/2完好程度均進入計數(shù)，每個樣點計數(shù)不小于300個，硅藻種屬的鑒定依據(jù)《中國淡水藻志》及《歐洲硅藻鑒定系統(tǒng)》^［22-23］。

2.4 分析評價

2.4.1 沉積物污染情況評價

利用有機指數(shù)和有機氮評價法進行對沉積物營養(yǎng)鹽污染情況的評判^［24］。有機指數(shù)一般用作反映水域沉積物環(huán)境狀況，有機氮為表征在氮水平上沉積物污染程度的指標。計算公式如下：

有機指數(shù)=有機碳（%）×有機氮（%）（1）

有機碳=有機質(zhì)÷1.724（2）

有機氮（%）=總氮（%）×95%（3）

2.4.2 硅藻相對豐度及優(yōu)勢度

計算硅藻豐度（P）并采用Mcnaughton優(yōu)勢度指數(shù)（Y）^［25］將優(yōu)勢度大于等于 0.02 且至少出現(xiàn)于 2 個樣品中的藻類列為樂安河源區(qū)水庫的優(yōu)勢種。

式中 N_sp——某個樣本中某種硅藻的個數(shù)；N——某個樣本中的總硅藻數(shù)量；f_i——某樣點中某種硅藻在所有樣點中出現(xiàn)的頻率。

2.4.3 硅藻指數(shù)

使用OMNIDIA 5.0（硅藻指數(shù)軟件）計算各樣點在18個硅藻指數(shù)的得分，分別為阿圖瓦-皮卡第硅藻指數(shù)（IDAP）、希佛爾和施賴納指數(shù)（SHE）、特定污染敏感指數(shù)（IPS）、勒克萊爾指數(shù)（IDSE）、硅藻富營養(yǎng)化指數(shù)（TDI）、南美大草原硅藻指數(shù)（IDP）、富營養(yǎng)污染硅藻指數(shù)（EPI-D）、瑞士硅藻指數(shù)（DI-CH）、斯雷德切克指數(shù)（SLA）、硅藻屬指數(shù)（IDG）、污染耐受百分比（PT）、羅特營養(yǎng)指數(shù)（ROTTtroph）、硅藻生物指數(shù)（IBD）、渡邊指數(shù)（WAT）、戴斯指數(shù)（DES）、歐盟硅藻指數(shù)（CEE）、洛博指數(shù)（LOBO）、羅特指數(shù)（ROTTsap）。除LOBO指數(shù)外，其余指數(shù)得分越高，水質(zhì)越好^［26-27］。

2.5 數(shù)據(jù)分析

對硅藻豐度及環(huán)境理化數(shù)據(jù)進行l(wèi)g（x+1）標準化處理以消除極值影響。首先對硅藻優(yōu)勢種數(shù)據(jù)進行主成分分析（principal components analysis，PCA）以判定硅藻群落變化的主要方向梯度，之后利用冗余分析（Redundancy Analysis，RDA）以此識別影響硅藻群落分布的驅(qū)動因子。本研究利用Excel 2010進行數(shù)據(jù)整理，利用Spass 22.0和 Canoco 5.0軟件對整理數(shù)據(jù)進行統(tǒng)計分析，最終結果通過Origin 2020、CorelDRAW及ArcGIS10.2軟件進行圖像繪制。

3 結果與討論

3.1 沉積環(huán)境分析

3.1.1 沉積柱沉積特征

由圖2可知，清華和石井水庫沉積物各理化指標在垂向深度上總體呈現(xiàn)逐步下降的趨勢。清華水庫沉積物TOC含量處于10～20.85 g/kg，隨深度含量逐步上升至18～11 cm處快速累積后并于10～0 cm區(qū)間逐步增加，有機質(zhì)含量高指示著清華水庫在這一時期的生產(chǎn)力也相應較高^［28］。沉積物TN的含量在0.97～1.65 g/kg，與TOC含量的變化趨勢相似。石井水庫中沉積物的TOC和TN含量分別為10.02～20.71、0.78～1.51 g/kg，于6～0 cm區(qū)間迅速增加，體現(xiàn)了短時間尺度上營養(yǎng)鹽的大量輸入。清華和石井水庫中TP變化范圍分別為0.22～0.73、0.25～0.71 g/kg。在靠近表層處，石井水庫各理化指標都顯現(xiàn)出較高于清華水庫的含量，可能意味著石井水庫受到較多方面營養(yǎng)物質(zhì)的累積及沉降。通過沉積物C/N值可以較好判定水域有機質(zhì)和營養(yǎng)鹽的來源，2座水庫的C/N值在9.7～14.5、11.4～14.0，表明水庫沉積物內(nèi)有機質(zhì)及營養(yǎng)鹽來源的復雜性，主要來源于水庫中藻類和水生植物，既存在水庫種內(nèi)源因素的影響也包括局部外源因素的作用。

3.1.2 沉積物污染情況評價

有機指數(shù)（指示沉積物環(huán)境狀況）和有機氮評價法（氮水平上沉積物污染程度）結果見圖3。有機指數(shù)結果顯示，清華水庫在垂直深度上呈現(xiàn)較清潔—尚清潔的變化趨勢，其中表層與中部（12～19 cm）顯示尚清潔程度，與理化指標反映趨勢一致。石井水庫在垂直深度上與清華水庫變化趨勢一致為較清潔-尚清潔，與之不同的是，石井水庫只在表層部分呈現(xiàn)了尚清潔的狀態(tài)。而有機氮指數(shù)結果顯示清華水庫經(jīng)歷了尚清潔-有機污染-尚清潔的梯度變化，石井水庫也經(jīng)歷了尚清潔-有機污染的變化趨勢。將2種評價指標結果進行對比后發(fā)現(xiàn)有機氮污染評價等級普遍高于有機污染指數(shù)，由此結果推論兩座水庫受含氮元素污染程度普遍較高，其中清華和石井水庫的部分點位中甚至指示了有機污染，表明其水域中營養(yǎng)狀態(tài)逐漸升高。

樂安河上游水庫遭到強烈的流域開發(fā)，主要以農(nóng)業(yè)生產(chǎn)和旅游業(yè)為主，同時也承擔魚類的網(wǎng)箱養(yǎng)殖的壓力影響。其中清華水庫毗鄰清華鎮(zhèn)，屬于河道型水庫，而石井水庫靠近婺源縣中心，四周農(nóng)田覆蓋。但近年來，隨著水庫周圍禁止養(yǎng)殖以及治理政策的展開，兩水庫中有機物污染評價指數(shù)較之前有所下降。從清華和石井水庫沉積物狀況及硅藻組合分布可知，兩座水庫皆存在先升后降的富營養(yǎng)化趨勢。

3.2 沉積硅藻特征

3.2.1 硅藻分布

清華水庫25個沉積樣品共鑒定出硅藻2綱6目10科23屬162種，石井水庫37個沉積樣品共鑒定出19屬123種。將相對豐度大于等于2%的屬種定義為優(yōu)勢種，其中清華水庫沉積柱中以附生和底棲類硅藻占主要優(yōu)勢，石井水庫由附生和底棲類硅藻逐漸過渡至以浮游類硅藻為主。根據(jù)沉積硅藻豐度在垂向深度上的變化，利用聚類分析將清華和石井水庫整個沉積柱巖芯劃分為2個組合帶（圖4）。

a）組合帶Ⅰ（Q25—Q15；S37—S9）。清華水庫中Achnanthidium minutissima的含量在這一階段急劇上升后下降，含量最高達到20.6%。Navicula cryptocephala和Navicula sp.2含量逐漸下降，降幅變化相似最終豐度都降至3.7%。Meridion circulare常見于水位較低的淡水環(huán)境，其含量逐步上升從2%到8.3%。耐污染種Synedra ulna以低含量偶有出現(xiàn)，這一階段表明清華水庫受污染較少，水體較為清潔。石井水庫中Aulacoseira granulata豐度迅速增加以平均含量12.2%占據(jù)絕對優(yōu)勢，這一階段主要以附生和底棲種為主。富營養(yǎng)指示種Cyclotella meneghiniana豐度逐漸增長，高營養(yǎng)鹽指示種Aulacoseira amigua豐度逐步減少。Navicula minima豐度在這一階段在優(yōu)勢種中占比也較高，平均達到6.1%。Achnanthidium minutissima、Navicula sp.2、Navicula palea、Cymbella affinis、Eunotia minor廣泛分布。

b）組合帶Ⅱ（Q14—Q1；S8—S1）。清華水庫以附生營養(yǎng)種Cocconeis placentula呈現(xiàn)波動性起伏迅速上升占據(jù)絕對優(yōu)勢，平均豐度12.3%。耐污染種Gomphonema parvulum和Navicula cryptocephala含量逐步上升，表明清華水庫正受到多環(huán)境因素的污染輸入。營養(yǎng)物質(zhì)隨地表徑流匯入水庫中，攜帶著流域生活污水及農(nóng)業(yè)肥料使用殘渣，導致了水域的營養(yǎng)富集、水質(zhì)惡化及藻類的爆發(fā)。石井水庫在組合帶Ⅱ中以Aulacoseira granulata迅速上升成為絕對優(yōu)勢，喜富營養(yǎng)環(huán)境的Gomphonema parvulum、Synedra ulna僅在表層處于優(yōu)勢。在這一組合帶石井水庫沉積硅藻群落主要以富營養(yǎng)化及耐污染指示種為主，也側面反映了石井水庫的富營養(yǎng)化趨勢。其中石井水庫受網(wǎng)箱養(yǎng)殖影響較為嚴重，餌料投放過多以及魚類成長過程中的糞便增加，都會直接導致大量的營養(yǎng)物質(zhì)進入水域內(nèi)，使得營養(yǎng)化現(xiàn)象加劇。推斷每年網(wǎng)箱養(yǎng)魚投放大量飼料等現(xiàn)象很可能是石井水庫營養(yǎng)物質(zhì)積聚的原因，同時也為藻類生長提供了先決條件。

3.2.2 硅藻指數(shù)分析

聚類分析結果顯示，清華和石井水庫皆被分為2個類群，故所有樣點共被分為4個類群，記為QⅠ（Q15—Q25）、QⅡ（Q1—Q14）、SⅠ（S9—S37）、SⅡ（S1—S8）。由沉積物污染情況評價結果及硅藻組合優(yōu)勢種變化情況可知，各個組所指示的水生態(tài)系統(tǒng)狀況為組合帶Ⅰ優(yōu)于組合帶Ⅱ，即QⅠgt;QⅡ，SⅠgt;SⅡ。除LOBO指數(shù)外，其余指數(shù)得分隨水質(zhì)變差得分降低并劃分為5個等級：健康（大于17）、好（大于 15 且小于 17）、中等（大于 12 且小于 15）、較差（大于 9 且小于12）、極差（小于 9）。計算 4 個類群在 18 種硅藻指數(shù)下的得分，并利用箱型圖剔除與研究

區(qū)水環(huán)境狀況變化不相符的指標，篩選出適宜評價樂安河源區(qū)水庫水環(huán)境狀況的硅藻指數(shù)。

硅藻指數(shù)作為有效的水環(huán)境狀況評價方法之一，被廣泛應用于水域生態(tài)環(huán)境的水質(zhì)評價^［29-30］。然而，硅藻指數(shù)也存在著環(huán)境異質(zhì)性，不同地理環(huán)境下所適宜的硅藻指數(shù)也不盡相同^［31］。為了篩選出適宜樂安河源區(qū)進一步研究和評價的硅藻指數(shù)，選取了18個歐洲硅藻指數(shù)對清華和石井水庫的樣點進行計算以獲取最適水環(huán)境狀況的硅藻指數(shù)。圖5所示，IDAP、SHE、CEE、LOBO指數(shù)與環(huán)境狀況匹配度較高，呈現(xiàn)出了相似的水環(huán)境趨勢。IDAP與大多數(shù)指數(shù)得分、SHE與環(huán)境因子沒有顯著相關性，CEE、LOBO指數(shù)與其他硅藻指數(shù)和環(huán)境變量在plt;0.05上顯著相關。因此，CEE、LOBO指數(shù)可作為樂安河源區(qū)水環(huán)境狀況評價的有效硅藻指標。

3.2.3 硅藻驅(qū)動因子分析

利用主成分分析探討清華及石井水庫硅藻組合群落變化的主要方向與潛在環(huán)境變化梯度，清華水庫由廣布種Achnanthidium minutissima、Gomphonema patricki、Meridion circulare等迅速轉向以混合型富營養(yǎng)化指示種Gomphonema parvulum、Cocconeis placentula等為主的硅藻群落結構。石井水庫從混合型以附生底棲種為主的Navicula minima、Achnanthidium minutissima等轉變到以耐污種和富營養(yǎng)化種為主的Cyclotella meneghiniana、Cymbella affinis等群落結構。耐污種及富營養(yǎng)化指示種的數(shù)量激增指示著水庫營養(yǎng)狀況的發(fā)生突變，表示在硅藻組合Ⅱ帶的富營養(yǎng)化趨勢，同時也側面反映水庫向富營養(yǎng)化水體轉變的趨勢。組合帶Ⅰ和Ⅱ分別與環(huán)境因子具有一定的相關性，隨著組合帶的變化與環(huán)境因子的相關性愈發(fā)強烈，表明水庫受到的污染輸入逐漸增加，也側面反映出與硅藻組合帶指向相似的富營養(yǎng)化趨勢。

為確定影響水庫沉積硅藻分布的關鍵驅(qū)動因素，通過DCA和RDA分析手段探求硅藻組合群落分布的趨勢，提取影響硅藻分布的環(huán)境因子^［32］。圖6結果表明，營養(yǎng)指標TP、TN是清華水庫硅藻組合群落結構分布的主要驅(qū)動因子（p＜0.05），清華水庫主要接收了來自沿途城鎮(zhèn)居民生活的各項營養(yǎng)物質(zhì)的輸入，人口和經(jīng)濟的快速增長導致農(nóng)業(yè)產(chǎn)量及生活污水直接流入水體，化肥的施用和人畜糞便攜帶大量營養(yǎng)鹽，持續(xù)性的外源營養(yǎng)物質(zhì)輸入造成水庫沉積物中 TP 含量不斷升高，加快水體富營養(yǎng)化進程。TOC、TP是石井水庫巖芯中硅藻群落分布的主要驅(qū)動因子（p＜0.05），可知因過剩的營養(yǎng)物質(zhì)，藻類得以大量繁衍生長促進水域的初級生產(chǎn)力^［33］。石井水庫相對封閉的水庫條件圍網(wǎng)養(yǎng)殖、開墾種植等人類活動致使富營養(yǎng)化情況較為明顯，也可解釋水域內(nèi)TP所反映的生態(tài)響應。大量的研究也揭示了大量的人類活動對硅藻的物種組成豐度及硅藻群落變化具有重大影響^［34-35］。

在人為活動影響下，伴隨水庫中營養(yǎng)鹽的持續(xù)富集，水庫經(jīng)歷2個階段：①以附生和底棲類為主要的優(yōu)勢時期，Achnanthidium minutissima、Navicula cryptocephala、Fragilaria capucina等相對清潔種豐度較高指示水域內(nèi)以較清潔為主的天然水體狀況；②以快速增長的硅藻富營養(yǎng)指示種時期，水域內(nèi)呈現(xiàn)營養(yǎng)物質(zhì)急劇輸入的富營養(yǎng)趨勢。在短期尺度上，清華和石井水庫中Cyclotella meneghiniana、Gomphonema parvulum、Cocconeis placentula等富營養(yǎng)種的頻繁出現(xiàn)表明庫區(qū)內(nèi)生態(tài)系統(tǒng)受到內(nèi)外雙重影響的過程作用，結果顯示在近代沉積物含量較高和日益增長的農(nóng)作需求及網(wǎng)箱養(yǎng)殖等人類活動下加重了水質(zhì)向富營養(yǎng)化狀態(tài)邁進的趨勢。在所計算的硅藻指數(shù)中，數(shù)據(jù)表明CEE、LOBO指數(shù)是最適宜樂安河源區(qū)水庫的有效硅藻評價指標。以靠近婺源縣城市的樂安河源區(qū)域2座不同生態(tài)類型的水庫為例，樂安河源頭區(qū)應注重水庫中營養(yǎng)物質(zhì)輸入和人類活動影響問題。

3.3 討論

本研究利用生物環(huán)境代用指標硅藻，探求樂安河上游區(qū)水庫的環(huán)境特征及硅藻與之的相關性和驅(qū)動因素。利用底棲硅藻物種多樣性及水環(huán)境理化性質(zhì)，橫向?qū)钒埠由嫌螀^(qū)4座水庫短時間尺度的水環(huán)境質(zhì)量做出評價，得出底棲硅藻群落結構變化的關鍵因素；再利用沉積硅藻優(yōu)勢種分布特征及沉積環(huán)境理化指標，縱向分析各水庫長時間尺度下環(huán)境質(zhì)量，得出影響沉積硅藻群落結構的驅(qū)動因素。并結合背景、氣候數(shù)據(jù)及經(jīng)濟指標數(shù)據(jù)等對樂安河上游區(qū)水庫的變化及現(xiàn)狀探討成因，為了解樂安河上游區(qū)婺源縣生態(tài)環(huán)境狀況及后續(xù)研究提供科學依據(jù)，也為下游鄱陽湖流域的污染控制與治理研究打好上游基礎。但是，在研究過程中本論文受客觀條件制約，仍有不足之處。

a）在實際采樣過程中，由于水庫水位較深、采樣地形復雜及采樣條件工具的限制，在枯水期采樣時只能采集靠近庫尾處的樣品，未能在水庫中心進行樣品采集，故部分樣品只能表征近現(xiàn)代的環(huán)境特征，從而無法獲得更長尺度下更為精確的沉積及硅藻數(shù)據(jù)。

b）研究通過垂直深度下以1 cm為單位進行樣品分割對沉積硅藻及環(huán)境指標進行分析。但所得的數(shù)據(jù)只能定性反映其變化趨勢和關鍵影響因素，無法結合實際年代從時間角度上定量分析環(huán)境演變過程，無法更為精確且明晰樂安河上游區(qū)水庫的演變過程?？梢栽诤罄m(xù)的研究中結合測年手段以完善定量分析的缺憾。

c）本文主要選取水質(zhì)理化指標（TP、TN、NH₃-N、Chl-a、DO、COD_Mn、SD）7項水質(zhì)數(shù)據(jù)及沉積理化指標（TN、TP、TOC、C/N）4項沉積數(shù)據(jù)與硅藻進行綜合分析，探討影響其群落變化的驅(qū)動因素。但同時沒有考慮到其他因素（水溫、水深、總鉀、總硫等）對硅藻群落結構變化的影響。因此，在分析探討中可能存在一定程度的局限性而導致結果過于片面?；诖瞬蛔?，可以在今后的研究中考慮多環(huán)境因子綜合對硅藻群落結構變化的驅(qū)動機制。

參考文獻：

［1］MORELLóN M，VALERO-GARCéS B，GONZáLEZ-SAMPéRIZ P，et al.Climate changes and human activities recorded in the sediments of Lake Estanya （NE Spain） during the Medieval Warm Period and Little Ice Age［J］.Journal of Paleolimnology，2011，46（3）：423-452.

［2］LAN B，ZHANG D L，YANG Y P，et al.Diatom-based reconstructions of hydrological variations and the underlying mechanisms during the past 520 years in the central Tianshan Mountains［J］.Journal of Hydrology，2019，575：945-954.

［3］YANG J F，JI Y，YAN R Y，et al.Applicability of Benthic Diatom Indices Combined with Water Quality Valuation for Dish Lake from Nanjishan Nature Reserve，Lake Poyang［J］.Water，2020，12（10）：2732-2748.

［4］康文剛，陳光杰，王教元，等.大理西湖流域開發(fā)歷史與硅藻群落變化的模式識別［J］.應用生態(tài)學報，2017，28（3）：1001-1012.

［5］SOLAK C N，PESZEK U，YILMAZ E，et al.Use of Diatoms in Monitoring the Sakarya River Basin，Turkey［J］.Water，2020，12（3）.DOI：10.3390/w12030703.

［6］SOKAL M A，HALL R I，WOLFE B B.Relationships between hydrological and limnological conditions in lakes of the Slave River Delta （NWT，Canada） and quantification of their roles on sedimentary diatom assemblages［J］.Journal of Paleolimnology，2008，39（4）：533-550.

［7］HADLEY K R，PATERSON A M，HALL R I，et al.Effects of multiple stressors on lakes in south-central Ontario：15years of change in lakewater chemistry and sedimentary diatom assemblages［J］.Aquatic Sciences，2013，75（3）：349-360.

［8］CHIBA T，ENDO K，SUGAI T，et al.Reconstruction of Lake Balkhash levels and precipitation / evaporation changes during the last 2000 years from fossil diatom assemblages［J］.Quaternary International，2016，397：330-341.

［9］JI Y，WU P J，ZHANG J，et al.Heavy metal accumulation，risk assessment and integrated biomarker responses of local vegetables：a case study along the Le’an river［J］.Chemosphere，2018 ，199：361-371.

［10］ZHANG J，HAN L，JI Y，et al.Heavy metal investigation and risk assessment along the Le’An River from non-ferrous metal mining and smelting activities in Poyang，China［J］.Journal of Environmental Biology，2018，39（4）：536-545.

［11］JI Y，ZHANG J，LI X L，et al.Biomarker responses of rice plants growing in a potentially toxic element polluted region：A case study in the Le’An Region［J］.Chemosphere，2017，187：97-105.

［12］夏楨妍.樂安河水域重金屬污染下的硅藻標志物響應分析［D］.南昌：南昌工程學院，2019.

［13］王越.樂安河流域土壤重金屬污染評價及影響因素分析［D］.濟南：山東科技大學，2020.

［14］何柱錕，陳海洋，陳瑞暉，等.樂安河沉積物重金屬污染評價與來源解析［J］.北京師范大學學報（自然科學版），2020，56（1）：78-85.

［15］余楊，呂雅寧，王偉杰，等.樂安河中下游重金屬時空分布特征及風險評價［J］.環(huán)境科學，2020，41（2）：691-701.

［16］舒旺，王鵬，肖漢玉，等.鄱陽湖流域樂安河水化學特征及影響因素［J］.長江流域資源與環(huán)境，2019，28（3）：681-690.

［17］田云飛，安俊芳，陶士敏.重鉻酸鉀氧化-分光光度法測定土壤有機碳含量的研究［J］.現(xiàn)代化工，2020，40（4）：231-235.

［18］袁靜，季平揚，袁藝.修正的凱氏法測定土壤中的總氮［J］.環(huán)境科學與管理，2009，34（2）：143-146.

［19］陳元松，侯明韜，馬丹，等.堿熔-鉬銻抗分光光度法測定土壤中總磷［J］.中國標準化，2018（S1）：193-196.

［20］劉娟，段昌兵，唐紅渠，等.沉積物硅藻提純方法的優(yōu)化［J］.生態(tài)科學，2013，32（6）：763-768.

［21］周娟.近60年來仙女湖硅藻群落變化特征及富營養(yǎng)化驅(qū)動因素研究［D］.南昌：南昌工程學院，2019.

［22］克拉默，蘭格-貝爾塔洛.歐洲硅藻鑒定系統(tǒng)［M］.廣州：中山大學出版社，2012.

［23］魏印心.中國淡水藻類：系統(tǒng)，分類及生態(tài)［M］.北京：科學出版社，2006.

［24］向愛農(nóng)，蘇甜，劉方平，等.撫河故支河道沉積物污染物分布特征及污染評價［J］.中國農(nóng)村水利水電，2020（10）：117-120.

［25］薛浩，鄭丙輝，孟凡生，等.甘河著生藻類群落結構及其與環(huán)境因子的關系［J］.生態(tài)環(huán)境學報，2020，29（2）：328-336.

［26］KALYONCU H，ERBETCI B.Applicability of Diatom-Based Water Quality Assessment Indices in Dari Stream，Isparta Turkey［J］.World Academy of Science，Engineering and Technology，International Journal of Environmental，Chemical，Ecological，Geological and Geophysical Engineering，2013，7（6）：386-394.

［27］劉麟菲，徐宗學，殷旭旺，等.應用硅藻指數(shù)評價渭河流域水生態(tài)健康狀況［J］.北京師范大學學報（自然科學版），2016，52（3）：317-321.

［28］丁宗艷，張家武，楊盼盼，等.柴達木盆地托素湖不同位置巖芯沉積物指標對比及其環(huán)境意義［J］.湖泊科學，2020，32（1）：259-270.

［29］趙楊，陸欣鑫，馬煜，等.阿什河流域底棲硅藻群落特征及水環(huán)境健康評價［J］.海洋與湖沼，2020，51（2）：307-317.

［30］張昆，蔡德所，林金城，等.硅藻指數(shù)在龍江與刁江流域生態(tài)健康評價中的應用［J］.中國農(nóng)村水利水電，2021（4）：98-106.

［31］TAYLOR J C，VUUREN M ，PIETERSE A.The application and testing of diatom-based indices in the Vaal and Wilge Rivers，South Africa［J］.Water Sa，2007，33（1）：51-59.

［32］馬露莎.三江平原沼澤濕地沉積硅藻分布規(guī)律及其環(huán)境相關性研究［D］.長春：中國科學院大學（中國科學院東北地理與農(nóng)業(yè)生態(tài)研究所），2018.

［33］NOGUEIRA M G.Phytoplankton composition，dominance and abundance as indicators of environmental compartmentalization in Jurumirim Reservoir （Paranapanema River），So Paulo，Brazil［J］.Hydrobiologia，2000，431（2/3）：115-128.

［34］The diatoms：applications for the environmental and earth sciences［M］.London：Cambridge University Press，2010.

［35］RüHLAND K M，PATERSON A M ，SMOL J P .Lake diatom responses to warming：reviewing the evidence［J］.Journal of Paleolimnology，2015，54（1）：1-35.

（責任編輯：高天揚）