劉海平,葉少文,楊雪峰,張良松,鐘國輝,何勇平,李鐘杰**
(1:西藏大學(xué)農(nóng)牧學(xué)院,林芝860000)
(2:中國科學(xué)院水生生物研究所,淡水生態(tài)與生物技術(shù)國家重點實驗室,武漢430072)
(3:中國科學(xué)院大學(xué),北京100049)
(4:錫林浩特市第六中學(xué),錫林浩特026000)
(5:福建省水產(chǎn)技術(shù)推廣總站,福州350003)
(6:西藏自治區(qū)林芝地區(qū)農(nóng)牧局,林芝860000)
著生藻類(periphytic algae)又稱為周叢藻類,是一個微型級、復(fù)合的群落,在物質(zhì)循環(huán)過程中發(fā)揮著重要的作用[1].著生藻類種類和數(shù)量與其粘附的介質(zhì)以及水體深度有很大的關(guān)系[2],一般來說,附植藻類和附石藻類的種類較附泥藻類豐富[1],硅藻門著生藻類與大型水生植物存在著密切的聯(lián)系,非硅藻門著生藻類與大型水生植物聯(lián)系較為微弱[3].著生藻類群落會隨著光線以及水質(zhì)理化(pH值、水溫、離子組成、氯離子、總氮、總磷等)的變動表現(xiàn)出較為明顯的時空異質(zhì)性[4],因此可以利用著生藻類豐度來判別河流的水質(zhì)狀況[5].馬沛明[6]探索性證明了以絲狀藻類為主的著生藻類對氮磷具有較高的吸收率和良好的去除效果;孫巍[7]指出藻類的初始投加量與N、P營養(yǎng)物的去除呈正相關(guān),磷元素可以成為著生藻類生長的限制性因子,總磷值過低時藻類生長速率降低;況琪軍等[8]指出剛毛藻(Cladophora oligoclona)在凈化污染水體、修復(fù)受損湖泊及防治水體富營養(yǎng)化等方面具有潛在的應(yīng)用前景.鑒于著生藻類對水體生態(tài)環(huán)境的響應(yīng)的差異性和敏感性,不同水域須篩選出不同的生態(tài)敏感指示著生藻類[9],從而有效推動水域生態(tài)系統(tǒng)健康發(fā)展.
關(guān)于研究區(qū)域概況、采樣點設(shè)置、采樣時間、環(huán)境因子的采集和分析和數(shù)據(jù)統(tǒng)計和分析等內(nèi)容在文獻[10]中進行了介紹,這里就不再贅述.
著生藻類的采集和分析:預(yù)先放置20 cm×20 cm的玻璃板于采樣點的河道中,14 d后,用毛刷清洗玻璃板,收集附著在其上的著生藻類,用魯哥試劑固定,帶回實驗室進行定量分析,每個樣品觀察2片,每片30個視野,鑒定到種或?qū)伲?1-12].按照如下公式對著生藻類進行定量分析:
使用Shannon-Wiener多樣性指數(shù)(H')、物種豐富度、總豐度、出現(xiàn)頻率和相對豐度來判別尼洋河著生藻類的多樣性.
Shannon-Wiener多樣性指數(shù) H'[13]計算公式為:
式中,S為著生藻類種類數(shù),pi為著生藻類i占所有著生藻類的比例(%).
物種豐富度指某一采樣點或者季節(jié)出現(xiàn)的物種個數(shù).
總豐度指某一采樣點或者季節(jié)單位面積出現(xiàn)的著生藻類個體數(shù)量(cells/cm2).
出現(xiàn)頻率指某一著生藻類出現(xiàn)的次數(shù)占所有調(diào)查樣點數(shù)的百分比.
相對豐度指某一著生藻類個體數(shù)占所有物種個體數(shù)的百分比.
尼洋河共有著生藻類6門30科49屬,包括硅藻門、藍藻門、綠藻門、隱藻門、黃藻門和裸藻門,其中硅藻門9科21屬,綠藻門12科17屬,藍藻門5科7屬,裸藻門2科2屬,黃藻門和隱藻門各出現(xiàn)1科1屬(表1).
表1 尼洋河著生藻類出現(xiàn)頻率和相對豐度Tab.1 Occurrence frequency and relative abundance of periphytic algae species in Niyang River
將著生藻類出現(xiàn)頻率在90%以上的定義為優(yōu)勢種,以此可知,尼洋河著生藻類優(yōu)勢類群以硅藻為主(表1),脆桿藻科有4屬,包括脆桿藻屬、等片藻屬、針桿藻屬、蛾眉藻屬,舟形藻科的舟形藻屬,橋彎藻科的橋彎藻屬,異極藻科的異極藻屬,曲殼藻科的曲殼藻屬,菱板藻科的菱形藻屬,如上幾種硅藻相對豐度之和達97.35%,其中以曲殼藻屬相對豐度最大,達49.60%,其實是橋彎藻屬和菱形藻屬,分別占18.38%和6.65%.
用PCA方法分析尼洋河著生藻類時空特征(圖1),采樣點Ⅳ囊括的著生藻類最多,其次依次為采樣點Ⅲ、Ⅱ和Ⅰ.可以初步判斷,隨著海拔的升高,著生藻類的物種豐富度和總豐度呈遞減趨勢.季節(jié)方面,以冬季囊括的著生藻類最多,接下來是春季、夏季和秋季.在底質(zhì)方面,砂石底質(zhì)較黏土底質(zhì)囊括的著生藻類較多.
圖1 基于PCA分析尼洋河著生藻類時空特征(圖中數(shù)字為樣點標(biāo)記,其中1~4分別代表采樣點Ⅰ的春、夏、秋、冬;5~8分別代表采樣點Ⅱ的春、夏、秋、冬;9~12分別代表采樣點Ⅲ的春、夏、秋、冬;13~16分別代表采樣點Ⅳ的春、夏、秋、冬.S1為采樣點Ⅰ,S2為采樣點Ⅱ,S3為采樣點Ⅲ,S4為采樣點Ⅳ.第一主成分解釋率為34.6%(圖中左上圖第一條黑色柱所示),第二主成分解釋率為11.5%(圖中左上圖第二條黑色柱所示))Fig.1 Spatio-temporal characteristics for periphytic algae in Niyang River based on PCA
為了更好地區(qū)分尼洋河各個采樣點和季節(jié)著生藻類群落的時空特征以及差異性,采用Duncan檢驗法對各采樣點以及季節(jié)之間著生藻類物種豐富度、總豐度和Shannon-Wiener多樣性指數(shù)進行差異性分析.結(jié)果顯示,尼洋河著生藻類物種豐富度、總豐度和Shannon-Wiener多樣性指數(shù)在各個采樣點之間均無顯著差異(P>0.05),尼洋河中下游較中上游著生藻類物種豐富度大,交匯處(采樣點Ⅱ為巴河和尼洋河交匯處,采樣點Ⅳ為尼洋河和雅魯藏布江交匯處)的著生藻類總豐度較其他河段(采樣點Ⅰ和Ⅲ)大,以尼洋河中游為分界,中上游河段和中下游河段著生藻類Shannon-Wiener多樣性指數(shù)均呈下降趨勢(圖2).
圖2 尼洋河著生藻類時空變化特征(用Duncan法檢驗各采樣點以及季節(jié)之間著生藻類相關(guān)參數(shù)差異性,包括總豐度、物種豐富度和Shannon-Wiener多樣性指數(shù),不同字母表示差異達顯著水平(P<0.05))Fig.2 Spatio-temporal dynamic characteristics for periphytic algae in Niyang River
秋、冬季著生藻類物種豐富度存在顯著差異(P<0.05),并且以秋、春、夏和冬季順序物種豐富度逐漸遞增,著生藻類總豐度在各個季節(jié)之間不存在顯著差異(P>0.05),以冬、春季的總豐度最高,夏、秋季與冬、春季著生藻類Shannon-Wiener多樣性指數(shù)存在顯著差異(P<0.05)(圖2).夏、秋季的著生藻類Shannon-Wiener多樣性指數(shù)較高,在1.8以上,冬、春季的著生藻類Shannon-Wiener多樣性指數(shù)較低,在1.8以下.
圖3 基于CCA方法分析尼洋河著生硅藻(a)、綠藻(b)和其它藻類(c)的豐度、環(huán)境因子、樣點之間的關(guān)系Fig.3 Relationships among abundance of periphytic algae(a:Bacillariophyta;b:Chlorophyta;C:others),environment factors and sampling sites in Niyang River based on CCA
為了更加清晰地闡釋尼洋河著生藻類與環(huán)境因子之間的關(guān)系,我們將著生藻類分為3個類群與9項環(huán)境指標(biāo)進行CCA分析,這3個類群是硅藻、綠藻以及以藍藻、隱藻、裸藻、黃藻為一類的藻類,9項環(huán)境指標(biāo)包括表層水溫、表層pH、硬度、礦化度、表層溶解氧、總氮、銨態(tài)氮、總磷、總堿度(圖3),對于硅藻門,尼洋河部分硅藻與理化因子相關(guān)聯(lián),如:雙壁藻屬(Diploneis)的豐度與總磷相關(guān)聯(lián),窗紋藻屬(Epithemia)的豐度與銨態(tài)氮相關(guān)聯(lián),雙菱藻屬(Surirella)的豐度與pH值相關(guān)聯(lián);綠藻門部分藻類豐度與理化因子相關(guān)聯(lián),如:小球藻屬(Chlorella)、柵藻屬(Scenedesmus)、溪菜屬(Prasiola)的豐度與水溫相關(guān)聯(lián),小椿藻屬(Characium)的豐度與硬度相關(guān)聯(lián),轉(zhuǎn)板藻屬(Mougeotia)的豐度與溶解氧、礦化度相關(guān)聯(lián),新月藻屬(Closterium)的豐度與堿度相關(guān)聯(lián).藍藻門、隱藻門、裸藻門、黃藻門的藻類分布與尼洋河各河段以及理化因子3者之間關(guān)聯(lián)性均較弱.
用CART預(yù)測了尼洋河著生藻類總豐度與環(huán)境因子之間的關(guān)系(圖4a).當(dāng)水溫≥12.98℃時,尼洋河著生藻類總豐度為266.66 cells/cm2;當(dāng)水溫介于12.4~12.98℃之間時,尼洋河著生藻類總豐度為967.28 cells/cm2;當(dāng)水溫 <12.4℃且為冬季時,尼洋河著生藻類總豐度為4442.41 cells/cm2;當(dāng)水溫<12.4℃且非冬季,同時海拔<3068 m時,尼洋河著生藻類總豐度為29127.25 cells/cm2;當(dāng)水溫 <12.4℃且非冬季,同時海拔≥3068 m、礦化度≥92.06 mg/L時,尼洋河著生藻類總豐度為363.63 cells/cm2;當(dāng)水溫<12.4℃且非冬季,同時海拔≥3068 m、礦化度 <92.06 mg/L時,尼洋河著生藻類總豐度為203.02 cells/cm2.綜上所述,在一定的水溫范圍之內(nèi)(12.98℃≥水溫≥12.4℃),較低水溫時的著生藻類總豐度比較高水溫的高;低海拔水域的著生藻類總豐度較高海拔高,高礦化度水域的著生藻類總豐度較低礦化度的高.
用CART預(yù)測了尼洋河著生藻類Shannon-Wiener多樣性指數(shù)與環(huán)境因子之間的關(guān)系(圖4b).當(dāng)水溫<8.12℃且非冬季時,尼洋河著生藻類Shannon-Wiener多樣性指數(shù)為1.529;當(dāng)水溫 <8.12℃且非冬季時,尼洋河著生藻類Shannon-Wiener多樣性指數(shù)為1.471;當(dāng)水溫≥8.12℃且礦化度 <83.35 mg/L時,尼洋河著生藻類Shannon-Wiener多樣性指數(shù)為1.878;當(dāng)水溫≥8.12℃且礦化度≥83.35 mg/L,同時底質(zhì)為黏土?xí)r,尼洋河著生藻類Shannon-Wiener多樣性指數(shù)為1.918;當(dāng)水溫≥8.12℃且礦化度≥83.35 mg/L,同時底質(zhì)為砂石、硬度≥2.57°DH 時,尼洋河著生藻類Shannon-Wiener多樣性指數(shù)為1.780;當(dāng)水溫≥8.12℃且礦化度≥83.35 mg/L,同時底質(zhì)為砂石、硬度<2.57°DH時,尼洋河著生藻類Shannon-Wiener多樣性指數(shù)為1.391.綜上所述,在水溫低于8.12℃時,尼洋河著生藻類冬季的Shannon-Wiener多樣性指數(shù)較春季的高;另外,尼洋河底質(zhì)為黏土的著生藻類Shannon-Wiener多樣性指數(shù)較砂石的高;尼洋河高硬度河段的著生藻類Shannon-Wiener多樣性指數(shù)較低硬度河段高.
圖4 基于分類回歸樹分析尼洋河著生藻類總豐度(a)、Shannon-Wiener多樣性指數(shù)(b)與環(huán)境因子之間的關(guān)系(分支左側(cè)代表條件成立,右側(cè)代表相反條件成立)Fig.4 Relationships between total abundance(a)and Shannon-Wiener diversity index(b)of periphytic algae and environment factors in Niyang River based on CART
海拔決定了一個地區(qū)的溫度和光照等環(huán)境因素的變化,屬于宏觀尺度的環(huán)境因子,通常情況下,由于海拔的升高、水域溫度降低、冰凍期延長[14-15],物種的豐富度也隨之降低[16],優(yōu)勢種群隨之發(fā)生改變,并且隨著海拔的上升,水溫對著生藻類群落結(jié)構(gòu)的影響愈發(fā)重要[17].通常情況下,著生藻類物種豐富度與水溫呈正相關(guān)[18],不同著生藻類群落對溫度的響應(yīng)機制有所區(qū)別,如運動型(motile ecological guild)的硅藻對溫度較為敏感,而喜靜水型(low profile guild)以及喜流水型(high profile guild)的硅藻則受總氮、總磷等理化因子驅(qū)動作用更多一些[4].本研究中,尼洋河中下游較中上游(低海拔)著生藻類物種豐富度大,尼洋河交匯處(低海拔)的著生藻類總豐度較其他河段大.
隨著水溫的降低,尤其到了冬季,較雨季而言,著生藻類生物量呈現(xiàn)增加趨勢[19],加之在這段時間,降雨量或者融雪量極大地減少,河流處于處于枯水期,流速降低,流速會直接影響到著生藻類的群落結(jié)構(gòu).由于水體的沖刷作用,流速快的地方,著生藻類生物量較低[1],某些硅藻門著生藻類,如Didymosphenia geminata更偏好于生存在具有穩(wěn)定水流的環(huán)境之中[20].另外水流速度與河床的底質(zhì)存在著相關(guān)性.就尼洋河而言,但凡流速快的地方,底質(zhì)為砂石,流速慢的地方,底質(zhì)則為黏土或者細砂,往往粗糙的基質(zhì)上著生的著生藻類物種豐富度較光滑基質(zhì)上的大,而總豐度則沒有明顯的區(qū)別[21].因此,本研究中尼洋河底質(zhì)為黏土河段的著生藻類物種豐富度和總豐度較底質(zhì)為砂石的河段大.
值得注意的是,魚類的習(xí)性以及攝食行為也會對著生藻類群落結(jié)構(gòu)產(chǎn)生影響.研究資料表明,西藏不乏以著生藻類為食的本土魚類,如異齒裂腹魚[22]、拉薩裸裂尻[23]等,以及以著生藻類為食的外來魚,如鯽魚[24],從一定程度上講,著生藻類的豐度是影響以其為食的西藏魚類種群數(shù)量較為關(guān)鍵的生物因子[25].不容忽視的是,隨著水溫的降低,魚類的攝食行為也在減弱,在以著生藻類為食的魚類所生存的水域,著生藻類群落結(jié)構(gòu)在冬季也會因此產(chǎn)生響應(yīng).
本實驗中,礦化度是根據(jù)水體的殘渣量得來的,總氮和總磷是殘渣量重要的組成成分,也是影響著生藻類群落結(jié)構(gòu)的重要水環(huán)境因子[26],總磷濃度迅速上升會導(dǎo)致葉綠素a含量和著生藻類數(shù)量急劇增加[27],在較高的照度,較高的氮、磷濃度的水體中,著生藻類總豐度也比較高[28],部分非硅藻門著生藻類與總氮或者總磷存在負相關(guān)[18],當(dāng)總氮和總磷濃度均較高時,無法明確判斷其與著生藻類生物量的關(guān)系[19],鑒于一些著生藻類自身具有固氮作用,磷元素可以成為著生藻類生長的限制性因子,總磷值過低時藻類生長速率會降低[7].本研究中,尼洋河高礦化度河段的著生藻類總豐度較低礦化度的高.
殘渣量還與水體的渾濁度有著密切的聯(lián)系,4-7月為尼洋河漲水期,水體較其他時間渾濁,光線變?nèi)?,有資料表明,光線對著生藻類的生物量有著重要的影響,隨著光線的減弱,著生藻類的生物量也在減弱[29].本文中,尼洋河冬季以著生藻類為食的西藏本土魚類以及外來魚類攝食行為減弱,同時,此階段尼洋河水體能夠沐浴充足的光線,著生藻類總豐度和物種豐富度有所增加.
影響著生藻類多樣性的因子包括水流、水深、濁度、溶解氧和水溫等[30],值得強調(diào)的是,由于空間的異質(zhì)性,導(dǎo)致了環(huán)境條件的多樣化,生態(tài)系統(tǒng)中著生藻類多樣性主要依存于空間的異質(zhì)性[31],流域上游的著生藻類Shannon-Wiener多樣性指數(shù)較其他河段低[32].中間高度膨脹(mid-altitude bulge)理論指出,往往在中間海拔區(qū)域內(nèi),物種的Shannon-Wiener多樣性指數(shù)較高,在低海拔或者高海拔區(qū)域內(nèi),物種的Shannon-Wiener多樣性指數(shù)較低[33].本研究中,尼洋河著生藻類Shannon-Wiener多樣性指數(shù)符合中間高度膨脹理論,即:尼洋河中游著生藻類Shannon-Wiener多樣性指數(shù)最大,中上游河段和中下游河段呈下降趨勢.
圖5 尼洋河著生藻類Shannon-Wiener多樣性指數(shù)與總豐度關(guān)系擬合Fig.5 Relationship fitting between Shannon-Wiener diversity index and total abundance of periphytic algae in Niyang River
另外,著生藻類的Shannon-Wiener多樣性指數(shù)也可以折射出流域的污染情況以及季節(jié)的差異性等信息,如夏季著生藻類的Shannon-Wiener多樣性指數(shù)較低[34],往往隨著生藻類生物量增加,著生藻類的Shannon-Wiener多樣性指數(shù)有降低的趨勢[31].鑒于此,為了更清楚地闡釋Shannon-Wiener多樣性指數(shù)和總豐度之間的關(guān)系,擬合了對數(shù)函數(shù),即:
式中,a、b為常量,Y和X是變量.結(jié)果如圖5所示,擬合效果較好.同時,因擬合的函數(shù)為對數(shù)函數(shù),尼洋河著生藻類總豐度與Shannon-Wiener多樣性指數(shù)之間不存在線性相關(guān),而是與Shannon-Wiener多樣性指數(shù)的自然對數(shù)呈線性負相關(guān),說明不能用影響著生藻類種類和數(shù)量的因子來解釋著生藻類的多樣性.
[1]劉建康等.高級水生生物學(xué).北京:科學(xué)出版社,1999:260-277.
[2]Yang H,F(xiàn)lower RJ.Effects of light and substrate on the benthic diatoms in an oligotrophic lake:a comparison between natural and artificial substrates.J Phycol,2012,48(5):1166-1177.
[3]Schneider SC,Lawniczak AE,Picinska-Faltynowicz J et al.Diatoms and non-diatom benthic algae give redundant information?Results from a case study in Poland.Limnologica,2012,42(3):204-211.
[4]Stenger-Kovacs C,Lengyel E,Crossetti LO et al.Diatom ecological guilds as indicators of temporally changing stressors and disturbances in the small Torna-stream,Hungary.Ecol Indic,2013,24:138-147.
[5]凌旌瑾,顧詠潔,許春梅等.黃浦江和蘇州河的著生藻類與水質(zhì)因子關(guān)系的多元分析.環(huán)境科學(xué)研究,2008,21(5):184-189.
[6]馬沛明.利用著生藻類去除N、P營養(yǎng)物質(zhì)的研究[學(xué)位論文].武漢:中國科學(xué)院水生生物研究所,2005.
[7]孫 巍.著生藻類去除微污染水體中氮磷的實驗研究[學(xué)位論文].成都:西南交通大學(xué),2008.
[8]況琪軍,凌曉歡,馬沛明等.著生剛毛藻處理富營養(yǎng)化湖泊水.武漢大學(xué)學(xué)報:理學(xué)版,2007,53(2):213-218.
[9]袁信芳,施華宏,王曉蓉.太湖著生藻類的時空分布特征.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2006,25(4):1035-1040.
[10]劉海平,葉少文,楊雪峰等.西藏尼洋河水生生物群落時空動態(tài)及與環(huán)境因子關(guān)系:1.浮游植物.湖泊科學(xué),2013,25(5):695-706.
[11]朱蕙忠,陳嘉佑.中國西藏硅藻.北京:科學(xué)出版社,2000.
[12]胡鴻鈞,李堯英,魏心印等.中國淡水藻類.上海:上??茖W(xué)技術(shù)出版社,1980.
[13]Shannon EE,Weaver W.The mathematical theory of communication.London:University Illinois Press,1949:125.
[14]Rundle SD,Jenkins A,Ormerod SJ.Macroinvertebrate communities in streams inHimalaya,Nepal.Freshwater Biol,1993,30(1):169-180.
[15]Ao M,Alfred JRB,Gupta A.Studies on some lotic systems in the north-eastern hill regions ofIndia.Limnologica,1984,15(1):135-141.
[16]Suren AM.Macroinvertebrate communities of streams in western Nepal:Effects of altitude and land use.Freshwater Biol,1994,32(2):323-336.
[17]Eli-Anne L,Stein W,Johansen et al.Periphyton in running waters-long-term studies of natural variation.Hydrobiologia,2004,521(1/2/3):63-86.
[18]Stancheva R,F(xiàn)etscher AE,Sheath RG.A novel quantification method for stream-inhabiting,non-diatom benthic algae,and its application in bioassessment.Hydrobiologia,2012,684(1):225-239.
[19]Kohler TJ,Heatherly TN,El-Sabaawi RW et al.Flow,nutrients,and light availability influence Neotropical epilithon biomass and stoichiometry.Freshwater Science,2012,31(4):1019-1034.
[20]Kirkwood AE,Shea T,Jackson LJ et al.Didymosphenia geminata in two Alberta headwater rivers:an emerging invasive species that challenges conventional views on algal bloom development.Can J Fish Aquat Sci,2007,64(12):1703-1709.
[21]Schneck F,Schwarzbold A,Melo A.Substrate roughness affects stream benthic algal diversity,assemblage composition,and nestedness.J N Am Benthol Soc,2011,30(4):1049-1056.
[22]季 強.異齒裂腹魚食性的初步研究.水利漁業(yè),2008,28(3):51-53,82.
[23]楊學(xué)峰.拉薩裸裂尻魚的食性及食物選擇的研究[學(xué)位論文].武漢:華中農(nóng)業(yè)大學(xué),2011.
[24]陳 鋒,陳毅峰.拉薩河魚類調(diào)查及保護.水生生物學(xué)報,2010,34(2):278-285.
[25]季 強.六種裂腹魚類攝食消化器官形態(tài)學(xué)與食性的研究[學(xué)位論文].武漢:華中農(nóng)業(yè)大學(xué),2008.
[26]文 航,蔡佳亮,蘇 玉等.滇池流域入湖河流豐水期著生藻類群落特征及其與水環(huán)境因子的關(guān)系.湖泊科學(xué),2011,23(1):40-48.
[27]王朝暉,胡 韌,谷陽光等.珠江廣州河段著生藻類的群落結(jié)構(gòu)及其與水質(zhì)的關(guān)系.環(huán)境科學(xué)學(xué)報,2009,29(7):1510-1516.
[28]Sanches LF,Guariento RD,Caliman A et al.Effects of nutrients and light on periphytic biomass and nutrient stoichiometry in a tropical black-water aquatic ecosystem.Hydrobiologia,2011,669(1):35-44.
[29]Bowes MJ,Ings NL,McCall SJ et al.Nutrient and light limitation of periphyton in the River Thames:Implications for catchment management.Sci of the Total Environ,2012,434(SI):201-212.
[30]Sharma RC,Bahuguna M,Chauhan P.Periphytonic diversity in Bhagirathi:Preimpoundment study of Tehri Dam Reservoir.J Environ Sci Eng,2008,50(4):255-262.
[31]Saravia LA,Giorgi A,Momo F.Multifractal spatial patterns and diversity in an ecological succession.PloS One,2012,7(3):1-8.
[32]Mihalic KC,Vilicic D,Ahel M et al.Periphytic algae development in the upper reach of the Zrmanja Estuary(Eastern A-driatic Coast).Vie Milieu,2008,58(3/4):203-213.
[33]賀金生,陳偉烈.陸地植物群落物種多樣性的梯度變化特征.生態(tài)學(xué)報,1996,7(l):91-99.
[34]Morin S,Duong TT,Dabrin A et al.Long-term survey of heavy-metal pollution,biofilm contamination and diatom community structure in the Riou Mort watershed,South-West France.Environ Pollut,2008,151(3):532-542.