漳江口紅樹植物群落周轉(zhuǎn)對大型底棲動物群落結(jié)構(gòu)的影響

陳國貴,王文卿,谷 宣, 洪 偉, 林宇楓,王 瑁 ,*

1 廈門大學(xué)環(huán)境與生態(tài)學(xué)院 廈門大學(xué)濱海濕地生態(tài)系統(tǒng)教育部重點實驗室, 廈門 361102 2 臺灣海峽海洋生態(tài)系統(tǒng)國家野外觀測研究站, 廈門 361102

目前,地上與地下的關(guān)系已經(jīng)成為生態(tài)學(xué)領(lǐng)域的研究熱點[1]。底棲動物是濕地生態(tài)系統(tǒng)中重要的消費者和分解者,對生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量流動的維持及濕地生態(tài)系統(tǒng)結(jié)構(gòu)和功能的穩(wěn)定都具有重要作用[2- 4]。底棲動物群落受到植物群落組成和沉積物理化性質(zhì)的影響[3,5-6],同時底棲動物群落也通過反饋作用影響植物群落[7-8]。因此,底棲動物是濕地生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量流動的中心環(huán)節(jié),也是濕地生態(tài)系統(tǒng)演化的重要驅(qū)動因子。

紅樹林濕地素有“地球之腎”的美稱[9],并以高生產(chǎn)力、高歸還率和高生物多樣性著稱[9],是存儲“藍(lán)碳”的寶庫[10]。紅樹林植被通常會促進(jìn)潮間帶生境的異質(zhì)性進(jìn)而影響林下底棲動物的分布和生存[5-6]。不同紅樹植被類型和林型對大型底棲動物群落密度、生物量以及功能多樣性具有顯著影響[5-6]。另外,由于紅樹林植被變化引起的生境變化也將導(dǎo)致大型底棲群落的變化[3,6]。受城市化和人類活動的影響,紅樹林濕地嚴(yán)重退化,表現(xiàn)為濕地面積嚴(yán)重萎縮、植物群落退化、功能明顯變化[9,11-12],但是目前對紅樹林濕地生態(tài)系統(tǒng)的地下結(jié)構(gòu)和功能變化趨勢并不清楚。因此,對紅樹林濕地生態(tài)系統(tǒng)變化過程中不同大型底棲動物群落的研究具有重要的生態(tài)學(xué)意義。

本研究以福建漳江口國家級紅樹林自然保護(hù)區(qū)中主要的紅樹林大型底棲動物群落為研究對象,旨在探究紅樹林濕地生態(tài)系統(tǒng)中植物群落周轉(zhuǎn)(即,物種替換和豐度變化)對大型底棲動物群落周轉(zhuǎn)的影響。研究結(jié)果不僅能夠提供紅樹林濕地底棲動物群落特征的基礎(chǔ)資料,還能揭示城市化背景下紅樹林濕地生態(tài)系統(tǒng)地下過程的變化趨勢,為紅樹林生態(tài)系統(tǒng)健康評價和退化紅樹林濕地的生態(tài)修復(fù)提供科學(xué)依據(jù)。

1 研究區(qū)域概況

福建漳江口國家級紅樹林自然保護(hù)區(qū)位于福建省南部的云霄縣境內(nèi),地處東經(jīng)117°24′—117°30′,北緯23°53′—23°56′,總面積2360 hm2,是國家級紅樹林濕地保護(hù)區(qū)的自然地理分布北界。保護(hù)區(qū)屬亞熱帶海洋季風(fēng)氣候,年均氣溫21.3 ℃,年均降雨量1714.5 mm,年均蒸散發(fā)量1718.4 mm。紅樹植物的優(yōu)勢種有秋茄(Kandeliaobovata)、桐花樹(Aegicerascorniculatum)、白骨壤(Avicenniamarina)和老鼠簕(Acanthusilicifolius)等。近幾十年來,在城市化和人類活動等因素的影響下,漳江口紅樹林濕地不斷退化,如水質(zhì)惡化、生境破碎化、植物群落結(jié)構(gòu)變化等問題層出不窮。

2 研究方法

2.1 樣地設(shè)置

2017年12月—2018年1月及2018年6月—10月,本研究在福建漳江口國家級紅樹林自然保護(hù)區(qū)設(shè)置了6條斷面(圖1)。每條斷面設(shè)置高潮帶、低潮帶和中潮帶3個采樣點,每個采樣點設(shè)置3個10 m×10 m的植物樣方進(jìn)行植物群落調(diào)查,記錄樹種和植株數(shù)等指標(biāo)。每個植物樣方設(shè)置三個重復(fù),并對沉積物水質(zhì)的理化參數(shù)進(jìn)行原位測定,包括鹽度、pH、溶解氧、氧化還原電位、電導(dǎo)率和溫度。每個植物樣方采集3—4個底泥樣品進(jìn)行大型底棲動物調(diào)查,樣方規(guī)格為25 cm × 25 cm,深度為30 cm。泥樣用0.5 mm篩網(wǎng)淘洗,挑選出軟體動物和蟹類樣本固定在5%福爾馬林溶液中,帶回實驗室分類鑒定。所收集軟體動物樣品的處理按照《海洋調(diào)查規(guī)范第6部分: 海洋生物調(diào)查》(GB/T 12763.6—2007)進(jìn)行。

圖1 福建漳江口紅樹林采樣點分布圖Fig.1 Sampling sites in Zhangjiang estuary, Fujian province

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

用基于蒙特卡洛的零模型模擬分別分析軟體動物和蟹類的共存格局(非隨機(jī)/隨機(jī))[13-14]。采用R軟件v4.0.2版本“EcoSimR”包隨機(jī)模擬運行50000次,計算群落的C-score指標(biāo)[14]。采用置換檢驗判斷觀測值和模擬值之間的顯著差異(P<0.05)。如果C-score的觀測值比模擬值大,說明群落為非隨機(jī)分散結(jié)構(gòu)；如果C-score的觀測值比模擬值小,說明群落為非隨機(jī)聚集結(jié)構(gòu)；如果C-score的觀測值等于模擬值,說明群落為隨機(jī)結(jié)構(gòu)[13-14]。另外,基于Spearman相關(guān)性,使用“igraph”包對軟體動物共存網(wǎng)絡(luò)結(jié)構(gòu)進(jìn)行模塊化分析[15]。模塊化程度越高,表示群落結(jié)構(gòu)越復(fù)雜,抵抗干擾的能力越強。

采用R軟件v4.0.2版本“Biodiversity”和“geosphere”包分別計算軟體動物、蟹類和紅樹植物群落的Bray-Curtis不相似性,環(huán)境的歐式距離(Euclidean distance)以及采樣點的地理距離[16],以探測群落周轉(zhuǎn)、環(huán)境變化程度和地理尺度。群落的不相似性越高表示群落周轉(zhuǎn)程度越高。歐式距離越大表示環(huán)境異質(zhì)性和變化程度越大。采用廣義加性模型(Generalized additive models,GAM)分析底棲動物群落與植物群落、環(huán)境距離和地理距離的關(guān)系,以探究植物群落周轉(zhuǎn)、環(huán)境過程和空間過程對底棲動物群落周轉(zhuǎn)的影響。用Spearman相關(guān)系數(shù)進(jìn)行Mantel檢驗進(jìn)一步量化各個因子(環(huán)境、植被和空間)對底棲動物群落周轉(zhuǎn)的影響。

3 結(jié)果與分析

3.1 蟹類、軟體動物、植物群落的周轉(zhuǎn)和環(huán)境變化程度

本研究采樣點間的最小地理距離為8.5 m,最大距離為1628 m,平均地理距離為669.3 m(圖2)。漳江口紅樹林濕地的環(huán)境距離最小值為0.24,最大值為7.19,平均值為3.28(圖2)。表明,在1.6 km的地理尺度上,漳江口紅樹林濕地的環(huán)境變化較大,出現(xiàn)7.19個環(huán)境梯度。本研究共采集20種蟹類、12種軟體動物并記錄了5種紅樹植物。其中,蟹類群落的平均不相似性為0.89,軟體動物群落的不相似性為0.83,植物群落的不相似性為0.76(圖2)。說明漳江口蟹類、軟體動物和紅樹植物群落的周轉(zhuǎn)率較高。漳江口紅樹林濕地是研究群落隨環(huán)境梯度變化的理想生態(tài)系統(tǒng)。

圖2 漳江口采樣地理距離、環(huán)境距離和生物群落不相似性Fig.2 Geographical distance of sampling, environmental distance and dissimilarity of biological community in Zhangjiang Estuary

3.2 蟹類、軟體動物的共存格局

漳江口軟體動物群落的實測C-score值小于零模型的模擬值(置換檢驗P=0.043),表現(xiàn)為非隨機(jī)聚集群落格局。蟹類群落的實測C-score值大于零模型的模擬值(置換檢驗P=0.038),表現(xiàn)為非隨機(jī)分散群落格局。軟體動物和蟹類的共存網(wǎng)絡(luò)都很復(fù)雜,但蟹類的網(wǎng)絡(luò)模塊化程度比軟體群落高(5:3；圖3)。說明軟體動物群落和蟹類群落可能受到不同生態(tài)過程地驅(qū)動。

圖3 蟹類和軟體動物群落的物種共存格局Fig.3 Co-occurrence pattern of crabs and mollusks

3.3 蟹類和軟體動物群落與植被、環(huán)境和空間因子的關(guān)系

廣義加性模型結(jié)果顯示,環(huán)境距離、地理距離和植物群落不相似性及其兩兩交互作用均與蟹類和軟體動物群落不相似性存在顯著的正相關(guān)關(guān)系(P<0.05；圖4)。單因素效應(yīng)分析顯示,環(huán)境距離對軟體動物群落的影響效應(yīng)為線性效應(yīng)(自由度等于1；圖4)。地理距離和植物群落不相似性對軟體動物群落的影響則呈非線

圖4 基于廣義加性模型,環(huán)境距離、地理距離和植物群落不相似性以及其兩兩交互作用與底棲動物群落不相似性的相關(guān)性Fig.4 Based on Generalized additive models, the correlations between the dissimilarity of macrobenthic communities and environmental distance/geographical distance/dissimilarity of mangrove community/their interactions

性效應(yīng)(自由度為6.88和5.77,大于1；圖4)。環(huán)境距離、地理距離和植物群落不相似性對蟹類群落的影響均為線性效應(yīng)(自由度為1；圖4)。雙因素效應(yīng)分析顯示,環(huán)境距離、地理距離和植物群落不相似性的兩兩交互作用均對蟹類和軟體動物群落產(chǎn)生非線性效應(yīng)(自由度大于1；圖4)。說明蟹類和軟體動物群落均受到環(huán)境過程、空間過程和植物群落周轉(zhuǎn)過程的影響,并以非線性效應(yīng)為主。

Mantel檢驗的分析結(jié)果顯示,溫度、pH、氧化還原電位、秋茄、白骨壤和空間因子顯著影響蟹類群落(P<0.05；圖5)。溫度、鹽度、氧化還原電位、秋茄、老鼠簕和空間因子則顯著影響軟體動物群落(P<0.05；圖5)。溫度、氧化還原電位、秋茄和空間距離是蟹類群落和軟體動物群落的共同影響因子,pH和白骨壤顯著影響蟹類群落,鹽度和老鼠簕則顯著影響軟體動物群落。

圖5 基于Mantel檢驗,分析蟹類和軟體動物群落與各個因子的相關(guān)性Fig.5 Based on the Mantel test, the correlations between crab and mollusk community and each factor

4 討論和結(jié)論

基于距離方法(distance-based method)的beta多樣性表示沿著某一環(huán)境梯度物種替代的程度或速率、物種周轉(zhuǎn)率、生物變化速度,可以靈敏地探測群落某個物種的得失,是研究群落隨環(huán)境變化的有效方法[17-18]。本研究結(jié)果表明,在漳江口紅樹林濕地1.6 km地理尺度的7.19個環(huán)境梯度上(圖2),植物群落平均物種周轉(zhuǎn)率為76%(圖2),蟹類和軟體動物群落平均物種周轉(zhuǎn)率分別為89%和83% (圖2)。高群落周轉(zhuǎn)率跟紅樹林濕地生態(tài)系統(tǒng)的特點有關(guān)。首先,紅樹林濕地是潮汐周期性淹沒的開放生態(tài)系統(tǒng),這不僅為紅樹林濕地提供物種補充還可以為在此定殖的海洋物種提供源源不斷的食物[8]。其次,紅樹林濕地生態(tài)系統(tǒng)生境異質(zhì)性很高,導(dǎo)致物種的空間分布差異巨大[2]。此外,有研究表明,不同生物類群存在功能差異,因此這些生物類群的群落結(jié)構(gòu)以及周轉(zhuǎn)率差異明顯[19-20]。例如：本研究發(fā)現(xiàn)漳江口紅樹林濕地蟹類群落結(jié)構(gòu)為非隨機(jī)分散格局,而軟體動物群落結(jié)構(gòu)為非隨機(jī)聚集格局(圖3)；蟹類群落的網(wǎng)絡(luò)模塊化程度比軟體動物群落的網(wǎng)絡(luò)模塊化程度高(圖3)；活動性較高的蟹類群落周轉(zhuǎn)率比活動性較低的軟體動物群落周轉(zhuǎn)率高6%(圖2)。因此,將蟹類和軟體動物整合在大型底棲動物的框架下研究必將導(dǎo)致一些重要生態(tài)過程被忽略。

在本研究中,隨著植物群落不相似性的增加,軟體動物和蟹類群落組成也分別發(fā)生明顯的線性和非線性變化(P<0.05；圖4)。其中,秋茄、白骨壤和老鼠簕種群周轉(zhuǎn)是紅樹林大型底棲動物群落的驅(qū)動因子(P<0.05；圖5)。這些結(jié)果表明紅樹植物群落周轉(zhuǎn)促進(jìn)底棲軟體動物和蟹類群落的周轉(zhuǎn)。有研究表明,底棲動物種類分布與紅樹植物種類關(guān)系密切[3],例如在泰國紅樹林中,擬沼螺種群密度與紅樹植物群落密度呈顯著正相關(guān),相手蟹的分布與紅樹植物群落類型關(guān)系明顯[21]；在廣東湛江紅樹林不同的植物群落(白骨壤+桐花樹群落、桐花樹群落和木欖(Bruguieragymnorhiza)+桐花樹群落)中,大型底棲動物群落的物種數(shù)、種群密度、豐富度指數(shù)和多樣性指數(shù)均存在顯著差異[7]；在福建九龍江口,秋茄和桐花樹恢復(fù)推動了大型底棲動物群落發(fā)育和演替[22-23]；紅樹林植被恢復(fù)初期改變底棲濱螺科動物群落結(jié)構(gòu)[24]。因此,植物群落是影響底棲動物群落結(jié)構(gòu)的重要因素[6-7,25]。植物群落為底棲動物提供充足的食物來源,滿足其生長發(fā)育需求,還能為底棲動物提供安全、適宜的棲息生境,由此形成對底棲動物自下而上的調(diào)控機(jī)制[6,26]。有研究發(fā)現(xiàn),高植被密度導(dǎo)致低地下動物群落多樣性,而中等和低植被密度促進(jìn)林下動物多樣性[27]。在漳江口紅樹林中,桐花樹、秋茄、白骨壤和老鼠簕為優(yōu)勢樹種。桐花樹的種群密度非常高(618株/100 m2),不利于大型底棲動物活動；而種群密度相對較低的秋茄(36株/100 m2)、白骨壤(54株/100 m2)和老鼠簕(74株/100 m2)生境則吸引更多大型底棲動物定殖。因此,超高密度的紅樹林并不利于維持底棲動物群落多樣性。種群密度較低的秋茄、白骨壤和老鼠簕不僅可以為底棲動物提供充足的生存空間,還可以為底棲動物提供多種微生境以維持底棲動物多樣性。例如,秋茄的板狀根是蟹類棲息的場所,白骨壤的密集呼吸根可以營造豐富的微生境以吸引底棲動物棲息和覓食。本研究強調(diào)紅樹植物群落變化顯著影響底棲軟體動物和蟹類群落結(jié)構(gòu)。

環(huán)境過濾是影響紅樹林底棲生物群落的重要生態(tài)過程[13,20]。底棲蟹類和軟體動物群落的不相似性均與環(huán)境距離呈顯著正相關(guān)非線性關(guān)系(圖4),其中,溫度、鹽度、pH、氧化還原電位等影響顯著(P<0.05；圖5),表明漳江口紅樹林大型底棲群落受環(huán)境過濾的影響顯著[13,20]。在全球變化的背景下,環(huán)境壓力是紅樹林濕地退化的主要問題之一[11]。例如,養(yǎng)殖污水和生活污水任意排放等。有研究表明,紅樹林大型底棲動物群落受環(huán)境的影響較大,主要包括溫度、海水鹽度、pH、潮汐的影響[3]。在東寨港和浙江西門島紅樹林中,大型底棲動物分布與溫度明顯相關(guān)[28]。在河口尺度上,溫度對底棲群落的影響主要表現(xiàn)在季節(jié)差異上。溫度對漳江口紅樹林大型底棲動物的影響跟本研究進(jìn)行夏季(濕季)和冬季(旱季)的采樣調(diào)查有關(guān)[29]。在漳江口紅樹林濕地,pH和鹽度與底棲動物群落呈高度相關(guān)(圖5)。但是,植物群落變化也會通過改變環(huán)境因素而間接影響底棲動物的生態(tài)過程[6,20,30]。紅樹植物群落結(jié)構(gòu)退化已然成為一個全球性的問題[11-12]。紅樹植物群落在變化過程中,會產(chǎn)生諸如水文條件變化等環(huán)境因素的改變,影響著與之伴生的海洋生物[5]。一方面,紅樹林通常導(dǎo)致生境沉積物酸化和鹽漬化,植物群落變化往往改變沉積物酸化和鹽漬化[5]。另一方面,漳江口紅樹林濕地的上游有大量的養(yǎng)殖塘和云霄縣居住區(qū),養(yǎng)殖廢水和生活污水的排放會降低水體鹽度和酸度。另外,最近的研究表明紅樹林大型底棲動物群落結(jié)構(gòu)受擴(kuò)散限制等隨機(jī)過程的影響[30]。在漳江口紅樹林濕地,底棲蟹類和軟體動物均與地理距離呈顯著正相關(guān)關(guān)系(圖4),且與空間因子關(guān)系顯著(圖5),表明漳江口紅樹林大型底棲動物群落也顯著受到擴(kuò)散限制的影響。一方面,紅樹植物群落結(jié)構(gòu)變化導(dǎo)致水文過程改變,進(jìn)而影響大型底棲動物的擴(kuò)散[30]。另一方面,水文條件的變化(如,酸化、鹽漬化或富營養(yǎng)化)不僅影響底棲動物的擴(kuò)散過程而且還影響底棲動物擴(kuò)散后的定殖[30]。因此,紅樹植物群落變化同樣影響底棲群落的擴(kuò)散限制等隨機(jī)生態(tài)過程。本研究強調(diào)紅樹植物群落變化對底棲群落生態(tài)過程的影響。

因此,紅樹植物群落變化(即物種更替和植株數(shù)量增減),導(dǎo)致理化環(huán)境因子和大型底棲動物的內(nèi)稟特征變化,進(jìn)而影響環(huán)境過濾和擴(kuò)散限制等生態(tài)過程,使大型底棲動物群落結(jié)構(gòu)發(fā)生變化。紅樹植物群落結(jié)構(gòu)變化會驅(qū)動地下底棲生物群落的結(jié)構(gòu),進(jìn)而改變底棲群落功能,最終影響整個紅樹林生態(tài)系統(tǒng)功能。