李妹洋,李龍沁,郝 廣,康瀕月,高 越,李洪遠(yuǎn),2①

(1.南開大學(xué)環(huán)境科學(xué)與工程學(xué)院,天津 300350;2.天津市跨介質(zhì)復(fù)合污染環(huán)境治理技術(shù)重點(diǎn)實(shí)驗(yàn)室,天津 300350)

生物入侵是一個(gè)全球性的環(huán)境問題,對自然生態(tài)系統(tǒng)和農(nóng)業(yè)生態(tài)系統(tǒng)有著深遠(yuǎn)的影響[1]。在國際貿(mào)易快速發(fā)展、各國交流逐漸增加的背景下,生物入侵逐漸成為物種多樣性降低的主要因素之一[2]。中國是世界上受外來入侵物種危害最嚴(yán)重的國家之一[3],作為生物多樣性大國,外來物種很容易在我國多種類型的生態(tài)系統(tǒng)中長期定殖[4],對我國生態(tài)安全、生物安全和糧食安全構(gòu)成重大威脅[5]。截至2020年8月,生態(tài)環(huán)境部發(fā)布的公報(bào)顯示,中國已出現(xiàn)660多個(gè)外來入侵物種。

土壤種子庫(soil seed bank)是植物生命周期中極為重要的環(huán)節(jié),承擔(dān)著地面植被更新的重要使命,同時(shí)也是植物入侵的重要方式和渠道[6]。多數(shù)研究表明,外來物種入侵會顯著影響土壤種子庫群落的多樣性和豐度[7]。土壤種子庫作為一種重要的種子資源,代表了一個(gè)可以發(fā)芽和建立植物群落的物種庫,它通過種子萌發(fā)和生長參與地上植被的自然更新[8],在生物多樣性保護(hù)和植被恢復(fù)方面具有重要的生態(tài)價(jià)值,所以了解土壤種子庫特征對于評估植被恢復(fù)潛力和演替軌跡具有重要意義[9]。國內(nèi)外已有部分研究涉及入侵植物對土壤種子庫的影響,如紫莖澤蘭(Ageratinaadenophora)[10]、加拿大一枝黃花(Solidagocanadensis)[11-12]和黃頂菊(Flaveriabidentis)[13]等。研究[1]表明,空心蓮子草的入侵降低了土壤種子庫中需要光照的物種和陸生物種的幼苗密度,最終影響空心蓮子草群落的自然種子恢復(fù)潛力。有關(guān)入侵植物對土壤種子庫影響的研究大多表明入侵植物對土壤種子庫物種組成、多樣性指數(shù)產(chǎn)生影響,進(jìn)而影響地上群落穩(wěn)定性及種間關(guān)系。這些研究結(jié)果為入侵植物管理和土著植被恢復(fù)提供了新的啟示。

長芒莧(Amaranthuspalmeri)作為入侵植物之一,已于2016年被納入第4批外來入侵物種名單。長芒莧原產(chǎn)于美國西部至墨西哥北部[14],是20世紀(jì)80年代中期傳入我國的一種莧科莧屬植物[15]。該植物主要生長于田園、郊野等地,入侵性和繁殖性極強(qiáng)[16],可使入侵地農(nóng)作物大量減產(chǎn),相比于同期一年生種子植物,長芒莧出苗率極高,其在適當(dāng)?shù)墓庹蘸退謼l件下就可以快速萌發(fā)。目前,已在華北、華南和華中等區(qū)域發(fā)現(xiàn)該物種的自然分布種群。國內(nèi)外針對長芒莧入侵機(jī)制的研究較少,且研究的關(guān)注點(diǎn)主要涉及其生理生態(tài)[17]及危害和防治[18]等,特別是在長芒莧除草劑的研制方面[19-20],也有研究涉及除草劑對土壤種子庫的抑制作用[21-22]。但針對長芒莧入侵機(jī)制等關(guān)鍵科學(xué)問題的研究鮮有報(bào)道,更未出現(xiàn)從土壤種子庫角度進(jìn)行的相關(guān)研究,關(guān)于入侵植物長芒莧對土壤種子庫群落特征的影響研究仍鮮見。因此,探究長芒莧對入侵地土壤種子庫組成與群落結(jié)構(gòu)的影響,有助于進(jìn)一步揭示長芒莧入侵機(jī)制,為莧屬入侵植物的科學(xué)防控提供依據(jù)。

筆者研究以入侵植物長芒莧為對象,選取天津市武清區(qū)一塊長芒莧自然入侵棄耕地為調(diào)查采樣地,通過分析入侵地土壤種子庫群落的物種組成、多樣性指數(shù)、種子密度和種間聯(lián)結(jié)性等,試圖闡明在長芒莧不同入侵壓力和不同土層深度條件下土壤種子庫物種組成及密度變化特點(diǎn),以及其對土壤種子庫物種多樣性及種間聯(lián)結(jié)性的影響。

1 材料與方法

1.1 研究區(qū)概況

調(diào)查采樣地位于天津市武清區(qū)廣賢路東側(cè),地理位置為39°20′50.97″ N,117°4′13.32″ E。該區(qū)地處華北沖積平原下段,海拔最高為13 m,最低為2.8 m。該區(qū)屬溫帶半濕潤大陸性季風(fēng)氣候區(qū),年均氣溫為11.6 ℃,年降水量為606 mm。土壤類型以潮土為主,土層深厚,土質(zhì)較黏。土壤pH偏弱堿性,土壤w(全氮)為1.7 g·kg-1,土壤w(全磷)為0.48 g·kg-1[23]。經(jīng)過前期考察,長芒莧入侵前樣地區(qū)域?yàn)樵耘嘧魑锎蠖箺壐?各樣地相鄰較近,且近6年來未開展過任何防除措施,確保無人為和地理因素干擾。代表性植物有苣荬菜(Sonchuswightianus)、狗尾草(Setariaviridis)、堿蓬(Suaedaglauca)、地膚(Kochiascoparia)和藜(Chenopodiumalbum)等常見鄉(xiāng)土植物。

1.2 樣品采集

采樣地長約400 m,寬約60 m,面積約為2.4萬m2。2021年3月初進(jìn)行土壤種子庫采樣工作,采樣時(shí)間選擇在種子萌發(fā)前,在表土采集區(qū)域依據(jù)長芒莧不同入侵壓力分別設(shè)置試驗(yàn)區(qū)組,根據(jù)長芒莧蓋度差異將試驗(yàn)區(qū)組劃分為輕度入侵壓力區(qū)組(<30%)和重度入侵壓力區(qū)組(>80%),每個(gè)區(qū)組分為2個(gè)樣地。每個(gè)樣地選取6個(gè)1 m×1 m代表性樣方,每個(gè)樣方按五點(diǎn)取樣法用取土環(huán)刀采集0～5、>5～10和>10～15 cm 3個(gè)土層深度樣品,將不同土層深度樣品分別放入塑料自封袋內(nèi)帶回實(shí)驗(yàn)室進(jìn)行萌發(fā)試驗(yàn)。

1.3 萌發(fā)試驗(yàn)

采用溫室萌發(fā)法,于土壤取樣后開始。首先,剔除土壤中雜物(砂石與植物殘?bào)w如根、莖、葉等),將過濾后的土壤樣品均勻平鋪在育苗平盤(53.5 cm × 28.5 cm × 5 cm)內(nèi),保證溫度、光照和水分等條件的一致性,定期澆水使土壤保持濕潤。溫度保持在25～30 ℃,光照為自然光,每2 d澆1次水并及時(shí)觀察種子萌發(fā)情況,待幼苗長出鑒定后,記錄并移除。定期人為翻動土壤,整個(gè)萌發(fā)試驗(yàn)開展到連續(xù)3周內(nèi)沒有新的幼苗出現(xiàn)即結(jié)束,整個(gè)萌發(fā)試驗(yàn)周期為70 d。

1.4 數(shù)據(jù)處理

采用Excel 2016對數(shù)據(jù)進(jìn)行初步處理,采用SPSS Statistics 26對數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,采用R 4.1.2的“ggplot2”軟件包制圖。采用廣義線性模型分析入侵壓力和土層深度及其交互作用對土壤種子庫的影響。利用非度量多維尺度分析(NMDS)和PERMANOVA分析,比較不同區(qū)組土壤種子庫群落結(jié)構(gòu)的空間變異性。土壤種子庫物種多樣性特征采用Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)、Simpson生態(tài)優(yōu)勢度指數(shù)和 Margalef豐富度指數(shù)表征;采用Spearman秩相關(guān)系數(shù)分析不同入侵壓力下土壤種子庫物種種間相關(guān)性。相關(guān)公式為

Shannon-Wiener多樣性指數(shù)(H):

(1)

Pielou 均勻度指數(shù)(J):

J=H/lnS。

(2)

Simpson優(yōu)勢度指數(shù)(D):

(3)

Margalef 豐富度指數(shù)(d):

d=(S-1)/lnN。

(4)

Spearman秩相關(guān)系數(shù)(r):

(5)

式(1)～(5)中,S為土壤種子庫中物種總數(shù);N為土壤種子庫中物種總個(gè)體數(shù);Pi為物種i個(gè)體數(shù)占種子庫物種總個(gè)體數(shù)的重要值;n為樣方數(shù);xij和xkj為物種i與k在樣方j(luò)中的秩。

2 結(jié)果與分析

2.1 不同入侵壓力和土層深度條件下土壤種子庫組成與密度

共調(diào)查兩個(gè)入侵壓力區(qū)組,每個(gè)區(qū)組2個(gè)樣地,共4個(gè)樣地。所有區(qū)組土壤種子庫共統(tǒng)計(jì)到30種植物,隸屬于17科30屬(表1)。

表1 所有樣地土壤種子庫密度及物種組成

表2 不同入侵壓力下土壤種子庫物種多樣性、豐富度和均勻度指數(shù)

根據(jù)廣義線性模型檢驗(yàn)結(jié)果,土壤種子庫密度在輕度入侵區(qū)組和重度入侵區(qū)組之間差異顯著。輕度入侵區(qū)組土壤種子庫有21種,隸屬于10科;禾本科植物幼苗最多,共2 285株,占幼苗總個(gè)體數(shù)的33.61%,其中,最多的為狗尾草(Setariaviridis)和馬唐(Digitariasanguinalis)。莧科植物幼苗占該區(qū)幼苗總數(shù)的15.78%,幼苗組成全部為長芒莧。重度入侵區(qū)組土壤種子庫有25種,隸屬于14科;莧科植物幼苗最多,共3 359株,占幼苗總個(gè)體數(shù)的87.27%,其中,幼苗組成全部為長芒莧。

根據(jù)植物生活型劃分,所有區(qū)組共有一、二年生草本植物21種,多年生草本植物10種。輕度入侵區(qū)組土壤種子庫種子數(shù)量占所有入侵地種子庫總量的63.85%,其中,一、二年生草本植物14種,多年生草本植物7種,兩者種子數(shù)量分別占該區(qū)總種子數(shù)量的77.15%和22.85%;重度入侵區(qū)組土壤種子庫種子數(shù)量占所有入侵地種子庫總量的36.15%,其中,一、二年生草本植物17種,多年生草本植物8種,兩者種子數(shù)量分別占該區(qū)總種子數(shù)量的99%和1%,一年生草本植物長芒莧種子數(shù)量占該區(qū)總種子數(shù)量的86.92%?？傮w來看,輕度入侵區(qū)組土壤種子庫種子數(shù)量比重度入侵區(qū)組多,兩種入侵壓力下植被類型均以一年生草本植物為主。

廣義線性模型表明,不同入侵壓力下土壤種子庫密度存在顯著差異(圖1),輕度入侵區(qū)組土壤種子庫密度顯著大于重度入侵區(qū)組(P<0.01),而長芒莧種子密度則表現(xiàn)為輕度入侵區(qū)組顯著小于重度入侵區(qū)組(P<0.01),不同土層深度間土壤種子庫中有活力的種子數(shù)量整體差異顯著(P<0.01),入侵壓力與土層深度之間無交互作用(P>0.05),所有入侵壓力下土壤種子庫中有活力的種子主要分布在上層(0～5 cm)。

就同一入侵區(qū)組而言,柱子上方英文小寫字母不同表示不同土層深度間種子密度差異顯著(P<0.05)。箱式圖上下端線表示樣本最大值和最小值,箱框線表示數(shù)據(jù)波動程度,箱內(nèi)線表示中位數(shù),箱外黑點(diǎn)表示樣本離群值。

2.2 不同入侵壓力和土層深度對土壤種子庫物種多樣性的影響

利用Shannon-Wiener多樣性指數(shù)、Simpson優(yōu)勢度指數(shù)、Pielou均勻度指數(shù)和 Margalef豐富度指數(shù)對土壤種子庫群落多樣性進(jìn)行分析。由表 2 可知,所有樣地土壤種子庫物種多樣性指數(shù)和豐富度指數(shù)均表現(xiàn)為輕度入侵區(qū)組>重度入侵區(qū)組,在入侵壓力下所有指數(shù)均表現(xiàn)為區(qū)組內(nèi)無差異(P>0.05),區(qū)組間差異顯著(P<0.01)。整體來看,該地土壤種子庫物種較為豐富。

采用廣義線性模型分析入侵壓力和土層深度對土壤種子庫物種多樣性、豐富度和均勻度指數(shù)等因變量的影響。廣義線性模型結(jié)果顯示入侵壓力與土層深度之間無顯著交互作用,表明不同土層深度土壤種子庫持久性受入侵壓力的影響不大。

圖2顯示,重度入侵壓力區(qū)組土壤種子庫Shannon-Wiener多樣性指數(shù)顯著低于輕度入侵壓力區(qū)組(P<0.05),其中,在不同土層深度間各區(qū)組土壤種子庫物種多樣性指數(shù)無顯著差異(P>0.05),且土層深度與入侵壓力之間無顯著交互作用。各區(qū)組各土層土壤種子庫Simpson優(yōu)勢度指數(shù)變化規(guī)律與Shannon-Wiener多樣性指數(shù)變化一致。

就同一幅圖中同一入侵區(qū)組而言,柱子上方英文小寫字母不同表示不同土層深度間某指標(biāo)差異顯著(P<0.05)。箱式圖上下端線表示樣本最大值和最小值,箱框線表示數(shù)據(jù)波動程度,箱內(nèi)線表示中位數(shù),箱外黑點(diǎn)表示樣本離群值。

與輕度入侵區(qū)組相比,重度入侵區(qū)組Pielou均勻度指數(shù)顯著降低(P<0.05,圖2)。不同土層深度間無顯著差異(P>0.05),土層深度與入侵壓力之間無顯著交互作用。

2.3 長芒莧入侵對土壤種子庫物種種間聯(lián)結(jié)性的影響

2.3.1不同入侵壓力下土壤種子庫群落結(jié)構(gòu)的NMDS及PERMANOVA分析

基于 Bray-Curtis 相似性系數(shù)對不同入侵壓力下土壤種子庫數(shù)據(jù)進(jìn)行非度量多維尺度分析(NMDS),結(jié)果(圖3)表明,輕度入侵區(qū)組與重度入侵區(qū)組重疊部分較少。Stress為0.114 8<0.2,說明該分析結(jié)果具有一定的解釋意義,該圖形能準(zhǔn)確地反映數(shù)據(jù)的真實(shí)分布。采用PERMANOVA分析對不同入侵壓力下土壤種子庫群落差異是否顯著進(jìn)行量化,PERMANOVA檢驗(yàn)結(jié)果為P= 0.001,表明不同入侵壓力下土壤種子庫群落結(jié)構(gòu)差異顯著,物種組成和分布相似性低,說明長芒莧入侵對土壤種子庫群落結(jié)構(gòu)影響較大。

A為輕度入侵區(qū)組,B為重度入侵區(qū)組。

2.3.2不同入侵壓力下土壤種子庫物種種間相關(guān)分析

對所有區(qū)組土壤種子庫中種子密度占優(yōu)勢(≥90%)的物種進(jìn)行Spearman秩相關(guān)分析,得到相關(guān)系數(shù)半矩陣(表3～5)。

表3 輕度入侵區(qū)組土壤種子庫物種相關(guān)系數(shù)半矩陣

表4 重度入侵區(qū)組土壤種子庫物種相關(guān)系數(shù)半矩陣

由表3～4可知,所有區(qū)組沒有與長芒莧呈顯著負(fù)相關(guān)的種對,輕度入侵區(qū)組負(fù)相關(guān)種對數(shù)占總種對數(shù)的33.0%,正相關(guān)種對數(shù)遠(yuǎn)大于負(fù)相關(guān)種對數(shù),與長芒莧相關(guān)的種對中呈不顯著負(fù)相關(guān)的植物有堿蓬(Suaedaglauca)、狗尾草(Setariaviridis)、葎草(Humulusscandens)、大豆(Glycinemax)和打碗花(Calystegiahederacea)。與輕度入侵區(qū)組相比,重度入侵區(qū)組中負(fù)相關(guān)種對數(shù)占比更高,占總種對數(shù)的39.74%,其中,大豆未在重度入侵區(qū)出現(xiàn),與長芒莧相關(guān)的種對中呈不顯著負(fù)相關(guān)的植物有堿蓬、狗尾草和早熟禾(Poaannua)。整體來看,輕度入侵區(qū)組土壤種子庫群落抵抗力和穩(wěn)定性更高。

當(dāng)將研究區(qū)整體作為分析對象時(shí),由表5可知,共有3個(gè)呈顯著和極顯著負(fù)相關(guān)的種對,這3個(gè)種對都與長芒莧有關(guān),其中,與長芒莧呈顯著負(fù)相關(guān)的植物為堿蓬,呈極顯著負(fù)相關(guān)的植物有狗尾草和大豆,進(jìn)一步說明長芒莧入侵對本地農(nóng)作物大豆、堿蓬及狗尾草有明顯的排斥作用。

表5 所有區(qū)組土壤種子庫物種相關(guān)系數(shù)半矩陣

3 討論

3.1 長芒莧入侵對土壤種子庫組成與密度的影響

研究區(qū)屬于城郊區(qū),輕度入侵區(qū)組與重度入侵區(qū)組土壤種子庫種子密度分別高達(dá)2 477和1 402?！-2,與天津城郊樣地土壤種子庫研究規(guī)模[24]類似,其中,重度入侵區(qū)組長芒莧密度是輕度入侵區(qū)組的3倍。外來物種對生態(tài)系統(tǒng)的入侵是對群落結(jié)構(gòu)的最大威脅之一[9],研究入侵種是否對群落變化產(chǎn)生影響,一般從群落組成及結(jié)構(gòu)入手進(jìn)行分析。筆者研究與張敏等[25]研究方法一致,基于Bray-Curtis距離對兩個(gè)區(qū)組土壤種子庫群落數(shù)據(jù)進(jìn)行NMDS和PERMANOVA分析,發(fā)現(xiàn)長芒莧入侵顯著干擾土壤種子庫群落結(jié)構(gòu)及物種組成,入侵壓力的不同導(dǎo)致土壤種子庫的差異性。

近地表取樣的土壤種子庫往往更加多樣化和豐富[26],筆者研究中隨著土層深度的增加,長芒莧種子密度顯著降低且主要分布在上層(0～5 cm),輕度入侵區(qū)組土壤種子庫密度為重度入侵區(qū)組的2倍,輕度入侵區(qū)組不同土層深度物種組成差異顯著且上層物種主要為長芒莧,而重度入侵區(qū)組中層(>5～10 cm)與下層(>10～15 cm)深度之間物種組成無顯著差異,認(rèn)為此與長芒莧的入侵年限有關(guān);另外,對于一些具有休眠特性的物種而言,采用誘萌法可能會帶來一些實(shí)驗(yàn)誤差。這也與長芒莧首先在上層土壤種子庫中聚集,進(jìn)而入侵整個(gè)土壤種子庫群落有關(guān)。相比于其他物種,長芒莧出苗期早、結(jié)種量大且適生性強(qiáng),種子很容易通過水進(jìn)行遠(yuǎn)距離漂浮傳播,在黏質(zhì)土中易停留在土壤表面從而更易聚集于種子庫。該植物具有較高的光合速率,可能存在較強(qiáng)的表型可塑性,因此,在無人為干擾的棄耕地土壤種子庫中更易成功定殖并快速擴(kuò)張,因此需要在土壤種子庫萌芽初期、長芒莧可能出現(xiàn)爆發(fā)性生長的時(shí)候采取一定的防控措施。

3.2 長芒莧入侵對土壤種子庫物種多樣性及種間關(guān)系的影響

筆者研究中取樣地為長芒莧入侵的大豆棄耕地,但近幾年來該棄耕地?zé)o人為干擾,是在自然演替過程中形成的土壤種子庫。這也可能是土壤種子庫中除入侵植物長芒莧外雜草占比較高而大豆占比較低的原因,DALTON等[27]研究認(rèn)為土壤種子庫的多樣性及豐富度高表明其對植物入侵具有一定的抵御能力[28]。

筆者研究中所有區(qū)組土壤種子庫中Shannon-Wiener多樣性指數(shù)、Simpson優(yōu)勢度指數(shù)和Pielou均勻度指數(shù)均表現(xiàn)為重度入侵區(qū)組小于輕度入侵區(qū)組,但Margalef豐富度卻表現(xiàn)為重度入侵區(qū)組大于輕度入侵區(qū)組?？紤]到重度入侵區(qū)組在統(tǒng)計(jì)個(gè)別物種時(shí)將其計(jì)入物種數(shù)量,但未將不同區(qū)組密度大小等本底值計(jì)入其中,個(gè)別雜草相對于整個(gè)區(qū)組而言,個(gè)體數(shù)量微乎其微,可忽略不計(jì),且周邊林下雜草有向樣地?cái)U(kuò)散的趨勢,這也可能是重度入侵區(qū)組物種豐富度大于輕度入侵區(qū)組的原因。

相較于輕度入侵區(qū)組,重度入侵區(qū)組土壤種子庫物種多樣性指數(shù)明顯降低,這與王瑞麗等[29]研究結(jié)果類似,即長芒莧與其他大部分入侵植物一樣,對土壤種子庫物種多樣性產(chǎn)生負(fù)面影響,表明長芒莧通過降低土壤種子庫物種多樣性、干擾其群落物種組成,從而快速入侵當(dāng)?shù)厣鷳B(tài)系統(tǒng)。RODRIGUEZ[30]研究認(rèn)為,將入侵物種與本土植物的種間作用納入未來的生態(tài)學(xué)研究,對于明確生態(tài)系統(tǒng)過程和生物多樣性變化至關(guān)重要。筆者研究對土壤種子庫相關(guān)物種進(jìn)行Spearman秩相關(guān)分析,結(jié)果表明大豆、堿蓬和狗尾草是受長芒莧負(fù)面影響較大的物種,考慮到長芒莧易出現(xiàn)在沿農(nóng)田邊緣的生境內(nèi)且在生長后期株高能達(dá)到2 m及以上,所以易入侵晚播作物如大豆和玉米等,BRUNDU等[14]與章茂林等[15]的研究中也提到了長芒莧的入侵易使大豆等糧食作物減產(chǎn);在花果期長芒莧植株高大、枝葉茂盛,會影響狗尾草和堿蓬等矮小植株的光合作用。

長芒莧與大豆、堿蓬和狗尾草在重度入侵區(qū)組與輕度入侵區(qū)組半矩陣系數(shù)表中呈現(xiàn)不顯著負(fù)相關(guān),當(dāng)將兩個(gè)區(qū)組作為整體分析時(shí),則呈顯著負(fù)相關(guān)。重度入侵區(qū)組沒有大豆出現(xiàn),這是因?yàn)樵摰貤壐约伴L芒莧入侵年限較長,致使大豆種子庫逐漸消亡。土壤種子庫持久性與種子大小有關(guān)[31],長芒莧種子小、易夾帶傳播,再加上大豆棄耕地?zé)o人為干擾,認(rèn)為這是其在短短幾年內(nèi)快速入侵土壤種子庫及當(dāng)?shù)厣鷳B(tài)系統(tǒng)的原因。探究長芒莧入侵農(nóng)田是否有利于當(dāng)?shù)厣锒鄻有缘姆€(wěn)定與發(fā)展以及通過何種方式作用于本土植物對于進(jìn)一步防控長芒莧十分必要[32],同時(shí)今后也要加強(qiáng)長芒莧在進(jìn)口糧谷和糧食產(chǎn)區(qū)內(nèi)的監(jiān)測預(yù)警和防范管控。