張 鸞,郭 偉,李素清,王曉軍
(1.山西大學(xué) 黃土高原研究所,太原 030006;2.山西大學(xué) 環(huán)境與資源學(xué)院,太原 030006)
晉西北丘陵風(fēng)沙區(qū)地處黃土高原東北部,降水少,氣溫低,風(fēng)沙大,原始植被殘留少,是山西省生態(tài)環(huán)境極脆弱區(qū)也是環(huán)京津塘地區(qū)防風(fēng)固沙的重要生態(tài)屏障。從20世紀(jì)70年代開始,這里陸續(xù)開展的“三北”防護(hù)林工程、退耕還林工程和京津風(fēng)沙源治理工程對生態(tài)環(huán)境的改善發(fā)揮了重要作用,也逐步形成了晉西北林業(yè)生態(tài)體系的新格局。然而隨著人工植被面積的增長,如何正確認(rèn)識和評估人工植被的生態(tài)恢復(fù)效應(yīng)受到研究人員的關(guān)注[1-4]。有研究表明:植被恢復(fù)過程中隨著演替時(shí)間的推移,土壤質(zhì)地細(xì)化,團(tuán)聚體增強(qiáng),土壤有機(jī)質(zhì)和養(yǎng)分含量增加[5-6],土壤微生物生物量提高[7],土壤質(zhì)量得到改善;植物群落演替與土壤水分之間的關(guān)系比其與土壤養(yǎng)分的關(guān)系更加密切[8],植被截留降水調(diào)節(jié)水文過程[9],但在階段性水分條件改善情況下形成的較高生物量對土壤水分消耗過大也導(dǎo)致了土壤干燥化的出現(xiàn)[10-12]。還有報(bào)道指出了植被恢復(fù)的水土保持、改善小氣候等效應(yīng)[13]。目前晉西北有關(guān)人工植被經(jīng)過多年恢復(fù)對環(huán)境響應(yīng)方面的研究相對較少[14-16]。曾經(jīng)“風(fēng)起黃沙飛,雨落洪成災(zāi)”的山西省右玉縣,經(jīng)過65年堅(jiān)持不懈造林,森林覆蓋率從0.3%上升到54%[17],是開展人工植被調(diào)查的理想?yún)^(qū)域。為此,本文選取右玉縣賈家窯流域開展人工植被分布格局和生態(tài)特征調(diào)查,利用數(shù)量生態(tài)分析方法探尋不同植被恢復(fù)模式的植被群落特征與演替動(dòng)態(tài)。
研究區(qū)位于晉西北丘陵風(fēng)沙區(qū)的右玉縣賈家窯流域,屬于典型的緩坡丘陵地貌,海拔1 386~1 489m,年均氣溫4.0℃,>10℃的年活動(dòng)積溫2 500℃,年均日照時(shí)數(shù)2 915h,無霜期平均104d,年均降水量420mm,降雨主要集中在6—9月,多年平均蒸發(fā)量1 761mm,年均風(fēng)速3m/s以上。研究區(qū)土壤類型以淡褐土和栗鈣土為主,有機(jī)質(zhì)含量不高。天然植被有針茅(Stipacapillata)、羊草(Leymuschinensis)、百里香(Thymusmongolicus)、白蓮蒿(Artenisiagmelinii)、阿爾泰狗娃花(Heteropappusaltaicus)等多年生草本,人工植被主要有油松(Pinustabulaeformis)、樟子松(Pinussylvestrismongolica)、華北落葉松(Larixprincipis-rupprechtii)、小葉楊(Populussimonii)、檸條錦雞兒(Caraganakorshinskii)和沙棘(Hippophaerhamnoides)等。
2015年7月在賈家窯流域海拔1 400~1 489m之間,采用系統(tǒng)取樣方法,根據(jù)海拔和坡向從坡底到峁頂隨機(jī)設(shè)置10m×10m喬木樣方75個(gè),5m×5m灌木樣方17個(gè),共計(jì)92個(gè)樣方。每個(gè)樣方調(diào)查記錄內(nèi)容主要包括植物種的名稱、蓋度、高度、喬木的胸徑和株數(shù),同時(shí)記錄樣方所在地的海拔高度、坡度、坡向、林齡、枯落物厚度、人類干擾程度等。剔除蓋度小于2%的植物種后92個(gè)樣方共得到32個(gè)植物種[18](表1),計(jì)算各植物種的重要值,形成32×92的重要值數(shù)據(jù)矩陣。
1) 喬木種的重要值(Pi)。計(jì)算公式[19]為:
Pi(喬木)=(相對密度+相對優(yōu)勢度+相對頻度)/300
相對密度=(某一種的個(gè)體數(shù)/全部種的個(gè)體總數(shù))×100
表1 賈家窯流域人工植被32個(gè)植物種的種類組成與生態(tài)類型Tab.1 32 species and their ecotypes of planting vegetation in Jiajiayao watershed
注:M為中生;XM為旱中生;MX為中旱生;X為旱生。
相對優(yōu)勢度=(某一種的基面積之和/全部種的基面積之和)×100
相對頻度=(某一種的頻度/全部種的頻度之和)×100
2) 灌木和草本植物種的重要值。計(jì)算公式為:
Pi(灌草)=(相對蓋度+相對高度)/200
相對蓋度=(某一種的蓋度/全部種的蓋度之和)×100
相對高度=(某一種的平均高度/全部種的平均高度之和)×100
3) 采用雙向指示種分析(TWINSPAN)和除趨勢對應(yīng)分析(DCA)[20],用PC-ORD5.0軟件完成賈家窯流域人工植被數(shù)量分類和排序分析。
4) 群落物種多樣性。計(jì)算公式[19]如下:
Patrick豐富度指數(shù),Dpa=S
Shannon-Wiener多樣性指數(shù),
式中:S為每個(gè)群落中的物種數(shù);Pi為每個(gè)物種的重要值;Pi=Ni/N,N為每個(gè)群落中全部物種的重要值之和;Ni為第i個(gè)種的重要值。
TWINSPAN等級分類法將92個(gè)樣方劃分為13個(gè)組,如圖1所示,N代表每一分組的樣方數(shù),D代表分類。依據(jù)《中國植被》的分類原則,13個(gè)組分別代表13個(gè)植物群叢,即:I小葉楊群叢(Ass.Populussimonii),Ⅱ沙棘群叢(Ass.Hippophaerhamnoides),Ⅲ檸條錦雞兒-針茅群叢(Ass.Caraganakorshinskii-Stipacapillata),Ⅳ華北落葉松-檸條錦雞兒群叢(Ass.Larixprincipis-rupprechtii-Caraganakorshinskii),V油松-檸條錦雞兒-針茅群叢(Ass.Pinustabulaeformis-Caraganakorshinskii-Stipacapillata),Ⅵ華北落葉松-針茅群叢(Ass.Larixprincipis-rupprechtii-Stipacapillata),Ⅶ華北落葉松-羊草群叢(Ass.Larixprincipis-rupprechtii-Leymuschinensis),Ⅷ樟子松-檸條錦雞兒群叢(Ass.Pinussylvestrismongolica-Caraganakorshinskii),Ⅸ檸條錦雞兒群叢(Ass.Caraganakorshinskii),X油松-檸條錦雞兒群叢(Ass.Pinustabulaeformis-Caraganakorshinskii),Ⅺ油松-羊草群叢(Ass.Pinustabulaeformis-Leymuschinensis),Ⅻ油松-白蓮蒿群叢(Ass.Pinustabulaeformis-Artenisiagmelinii),XIII油松群叢(Ass.Pinustabulaeformis)。13個(gè)群叢的主要特征如表2所示。
圖1 賈家窯流域92個(gè)樣方的TWINSPAN分類樹狀圖
表2 賈家窯流域13個(gè)群叢的主要特征Tab.2 Characteristics of 13 clusters of planting vegetation in Jiajiayao watershed
采用第一、二排序軸的特征值作樣方的二維排序圖(圖2),從圖中可以看出,92個(gè)樣方被分為13個(gè)區(qū),與TWINSPAN分類所產(chǎn)生的13個(gè)群叢相對應(yīng)。
圖2 92個(gè)樣方的DCA二維排序圖
DCA第一軸主要反映了海拔高度的梯度變化,即沿第一軸從左到右,海拔高度逐漸降低,群叢類型依次呈現(xiàn)出由群叢XIII→群叢Ⅴ→群叢Ⅲ→群叢Ⅱ過渡的趨勢,群落類型由中生的油松林向旱中生的檸條錦雞兒和沙棘灌叢過渡;DCA笫二軸主要反映了坡向的梯度變化,即沿第二軸從下到上,坡向由陰坡逐漸轉(zhuǎn)變?yōu)殛柶?群落類型依次呈現(xiàn)出由群叢Ⅶ→群叢Ⅳ、群叢Ⅴ和群叢Ⅱ→群叢Ⅸ→群叢Ⅷ過渡的趨勢。
圖3為13個(gè)群叢的物種豐富度指數(shù)和多樣性指數(shù)對比,圖中字母表示多重比較結(jié)果,同一列字母相同表示兩者差異性不顯著,字母不同,表示差異性顯著(p<0.01)。不同群叢的Patrick豐富度指數(shù)和Shannon-Wiener多樣性指數(shù)表現(xiàn)出基本一致的變化趨勢,油松-檸條錦雞兒群叢(X)、油松-檸條錦雞兒-針茅群叢(Ⅴ)、華北落葉松-檸條錦雞兒群叢(Ⅳ)、樟子松-檸條錦雞兒群叢(Ⅷ)和檸條錦雞兒群叢(Ⅸ)較大,油松-羊草群叢(Ⅺ)、華北落葉松-羊草群叢(Ⅶ)、華北落葉松-針茅群叢(Ⅵ)、小葉楊群叢(I)和油松-白蓮蒿群叢(Ⅻ)居后,沙棘群叢(Ⅱ)和檸條錦雞兒-針茅群叢(III)較小,油松群叢(XIII)最小。Simpson多樣性指數(shù)則呈現(xiàn)相反趨勢。
圖3 不同群叢物種多樣性指數(shù)對比
以Shannon-Wiener多樣性指數(shù)為例,進(jìn)行6種主要植物群落不同坡向?qū)Ρ?不難發(fā)現(xiàn)(圖3(d)),除退耕還林25a左右的華北落葉松純林在半陽坡和半陰坡Shannon-Wiener多樣性指數(shù)相差不大外,油松、樟子松、檸條錦雞兒和沙棘純林均呈現(xiàn)陰坡大于陽坡的趨勢。
1) 小葉楊、油松、樟子松、華北落葉松、檸條錦雞兒和沙棘是賈家窯流域人工植物群落的主要建群種。作者在晉北其他縣市的植被調(diào)查發(fā)現(xiàn),以上優(yōu)勢種也是晉西北丘陵風(fēng)沙區(qū)人工植物群落的主要建群種。
2) 從優(yōu)勢種的栽植時(shí)間和造林模式來看,研究區(qū)人工植被建設(shè)經(jīng)歷了從單純種植喬木、灌木到喬灌混交林的轉(zhuǎn)變。20世紀(jì)60—70年代大面積栽植小葉楊作為防風(fēng)固沙林,造林模式為純林,造林密度較大;70—80年代,油松、華北落葉松和樟子松逐漸替代小葉楊成為主要造林樹種,造林模式仍以純林為主;90年代以后造林模式有所改變,喬灌混交造林模式出現(xiàn)(群叢Ⅹ),原“小老樹”林多由油松、樟子松、檸條錦雞兒和沙棘以帶狀或塊狀方式更新改造,油松、華北落葉松純林則多由檸條錦雞兒隔行更新改造從而形成老幼異齡林結(jié)構(gòu)(群叢Ⅱ,Ⅳ,Ⅴ,Ⅷ)??傮w來看,目前研究區(qū)人工植被類型以鄉(xiāng)土種為主,群落結(jié)構(gòu)仍較為單一(13個(gè)群落中有8個(gè)為純林)。
1) DCA結(jié)果顯示,賈家窯流域人工植被群落分布主要與海拔和坡向有關(guān):檸條錦雞兒和沙棘灌叢主要分布在較低海拔特別是陡坡地段,隨著海拔升高,喬木林和喬灌混交林分布范圍擴(kuò)大。樟子松主要分布在干旱陽坡和半陽坡;華北落葉松分布在陰坡、半陰坡和半陽坡;小葉楊和沙棘主要在海拔較低的半陽坡和半陰坡生長;油松和檸條錦雞兒則在各個(gè)坡向都有種植。這表明研究區(qū)人工植被建設(shè)考慮了樹種的生物學(xué)特性和立地條件。
2) 實(shí)地調(diào)查得知,經(jīng)過多年演替,小葉楊已呈現(xiàn)出“小老樹”態(tài)勢;油松目前長勢較好,林齡在20a以上的油松林下多見幼苗更新;樟子松種植時(shí)間相對較晚,但生長速度較快,其問題在于經(jīng)過多次調(diào)查,未發(fā)現(xiàn)林下更新。華北落葉松樹高、胸徑生長都不太理想,這可能一方面因?yàn)槿A北落葉松更適應(yīng)高海拔地區(qū),另一方面與種植間距過密,后期缺乏撫育管護(hù)有關(guān)。研究區(qū)人工林林下草本植物種類單調(diào),主要優(yōu)勢種或伴生種有針茅、百里香、羊草、華北米蒿、白蓮蒿、阿爾泰狗娃花等。郁閉度較小的小葉楊林枯枝落葉層分解較快,林下草本類型較多,但多出現(xiàn)退化指征草本——瑞香狼毒;大多數(shù)油松林下無枯落物,地表裸露,土壤侵蝕較其他林型重。經(jīng)過訪問,右玉當(dāng)?shù)乩习傩沼胁墒坝退伤舍槷?dāng)燃料的習(xí)慣,交通相對較便利的油松林都不同程度地受到了干擾;樟子松林下也有人為擾動(dòng)的跡象,程度相對較輕;華北落葉松林郁閉度大,枯枝落葉層分解慢,地表形成較厚的枯枝落葉層。
1) 賈家窯流域不同植物群落的物種多樣性表現(xiàn)出較強(qiáng)的規(guī)律,群叢X,Ⅴ,Ⅳ和Ⅷ分別為油松、華北落葉松和樟子松與檸條錦雞兒混交形成的群落,林型結(jié)構(gòu)復(fù)雜、發(fā)育較好,因而物種豐富度和物種多樣性指數(shù)均較高。而群叢I,Ⅵ,Ⅶ和Ⅻ分別為小葉楊、華北落葉松和油松純林,林型結(jié)構(gòu)簡單,因而物種豐富度和物種多樣性指數(shù)較低。同種類型群落中,華北落葉松和油松純林地形條件和恢復(fù)年限一致時(shí),前者的shannon-Wiener多樣性指數(shù)較后者低,這是由于華北落葉松枯枝落葉層分解較慢,影響了林下草本發(fā)育的緣故。
2) 群叢Ⅱ和III分別是以沙棘和檸條錦雞兒占絕對優(yōu)勢的灌叢,草本層種類比較貧乏,物種豐富度和物種多樣性指數(shù)均顯著低于前述喬木純林和喬灌混交林群落。群叢Ⅸ也是檸條錦雞兒灌叢,但該群叢的物種豐富度和物種多樣性指數(shù)均較喬木純林高。比較發(fā)現(xiàn),群叢Ⅱ和III均位于較低海拔,而群叢Ⅸ的海拔高度與喬木純林的海拔高度相差不大。這說明,隨海拔高度增加,水熱條件較好,有利于物種向多樣化方向發(fā)展,因而具有較高的物種豐富度和物種多樣性。
3) 坡向使光熱和水分在微環(huán)境內(nèi)重新分配,對群落的物種豐富度和生物多樣性有著重要影響[20]。不同坡向的群落在相同恢復(fù)年限內(nèi)shannon-Wiener多樣性指數(shù)陰坡大于陽坡,這主要是由于陽坡比陰坡光照時(shí)數(shù)長,水分條件差,群落內(nèi)種間水分競爭激烈,低耐旱和喜陰物種難以存活,植被生長受到限制的緣故[21]。從林齡角度來看,對比同海拔(1 435 m)半陽坡退耕還林10a和25a左右的油松純林shannon-Wiener多樣性指數(shù)(1.41,1.33)和物種豐富度指數(shù)(7,6),前者均大于后者。這是由于早期退耕土壤水分和養(yǎng)分條件較好地為植物提供了良好的生長環(huán)境,群落內(nèi)物種種類增加明顯,隨著演替的進(jìn)行,地帶性次生草原群落形成,演替速率降低,物種多樣性指數(shù)的變化幅度較之前有所減小。
4) 群叢XIII油松純林各項(xiàng)指數(shù)均與其他群叢有顯著差異。這與右玉當(dāng)?shù)乩习傩詹墒坝退伤舍槷?dāng)燃料的習(xí)慣有關(guān),說明人類活動(dòng)過度干擾也能導(dǎo)致物種多樣性出現(xiàn)波動(dòng)。
1) 右玉縣賈家窯流域經(jīng)歷了半個(gè)多世紀(jì)的植被建設(shè),形成了以油松、樟子松、華北落葉松、小葉楊、檸條錦雞兒和沙棘為主要建群種的人工植被群落。
2) 研究區(qū)人工植被建設(shè)考慮了樹種的生物學(xué)特性和立地條件,但植被結(jié)構(gòu)存在純林比重較大、種植過密、人工植被管護(hù)不利等問題。
3) 應(yīng)加強(qiáng)現(xiàn)有植被群落結(jié)構(gòu)的優(yōu)化調(diào)控,更新改造小葉楊林;加大喬、灌、草復(fù)合配置比例,豐富群落層次結(jié)構(gòu);充分考慮土壤水分、養(yǎng)分等生態(tài)條件確定造林密度;加強(qiáng)人工植被的科學(xué)管護(hù),避免人為過度干擾影響植物群落發(fā)育;另外,對于樟子松在研究區(qū)不能實(shí)現(xiàn)自然更新的問題應(yīng)引起高度重視。
[1] Harris J.Soil microbial communities and restoration ecology:facilitators of followers?[J].Science,2009,325:573-574.
[2] Chazdon R L.Beyond deforestation:Restoring forests and ecosystem services on degraded lands[J].Science,2008,320(5882):1458-1460.
[3] 彭少麟.退化生態(tài)系統(tǒng)恢復(fù)與恢復(fù)生態(tài)學(xué)[J].中國基礎(chǔ)科學(xué)研究進(jìn)展,2001(3):18-24.
[4] 胡嬋娟,郭雷.植被恢復(fù)的生態(tài)效應(yīng)研究進(jìn)展[J].生態(tài)環(huán)境學(xué)報(bào),2012,21(9):1640-1646.
[5] Jiao F,Wen Z M,An S S.Changes in soil properties across a chronosequence of vegetation restoration on the Loess Plateau of China[J].Catena,2011,86:110-116.
[6] 王改玲,王青杵.晉北黃土丘陵區(qū)不同人工植被對土壤質(zhì)量的影響[J].生態(tài)學(xué)雜志,2014,33(6):1487-1491.
[7] 薛箑,劉國彬,戴全厚,等.侵蝕環(huán)境生態(tài)恢復(fù)過程中人工刺槐林(Robinia Pseudoacacia)土壤微生物量演變特征[J].生態(tài)學(xué)報(bào),2007,27(3):909-917.
[8] 焦峰,溫仲明,焦菊英,等.黃丘區(qū)退耕地植被與土壤水分養(yǎng)分的互動(dòng)效應(yīng)[J].草業(yè)學(xué)報(bào),2006,15(2):79-84.
[9] 馮迪.陜西省安塞縣退耕還林工程生態(tài)效益監(jiān)測與評價(jià)[D].北京:北京林業(yè)大學(xué),2010.
[10] Wang Y Q,Shao M A,Liu Z P.Large-scale spatial variability of dried soil layers and related factors across the entire Loess Plateau of China[J].Geoderma,2010,159(1-2):99-108.
[11] 劉庚,??〗?朱煒歆,等.晉西北黃土丘陵區(qū)人工林地土壤水分虧缺評價(jià)[J].水土保持通報(bào),2015,35(5):309-314.
[12] 梁海斌,史建偉,??〗?等.晉西北黃土丘陵區(qū)不同林齡檸條地土壤水分變化特征研究[J].干旱區(qū)資源與環(huán)境,2014,28(6):143-148.
[13] 李小英,段爭虎.黃土高原土壤水分與植被相互作用研究進(jìn)展[J].土壤通報(bào),2012,43(6):1508-1514.
[14] 王尚義,李素清,曹志敏,等.山西嵐縣皇姑梁小流域人工植物群落生態(tài)[J].生態(tài)學(xué)報(bào),2007,27(12):5098-5109.
[15] 趙德懷,李素清.晉西北丘陵風(fēng)沙區(qū)人工植被數(shù)量分類與排序研究[J].山西師范大學(xué)學(xué)報(bào):自然科學(xué)版,2011,25(1):103-109.
[16] 董林水.晉西基巖山地與黃土丘陵區(qū)生態(tài)過渡帶植被特征及生態(tài)恢復(fù)機(jī)理研究[D].北京:中國林業(yè)科學(xué)研究院,2005.
[17] 王濤,辛泰.綠水青山就是金山銀山——右玉縣傳承弘揚(yáng)右玉精神推進(jìn)生態(tài)文明建設(shè)紀(jì)實(shí)[N].朔州日報(bào).2017-5-9(4).
[18] 中國植被編輯委員會(huì).中國植被[M].北京:科學(xué)出版社,1983:143-156.
[19] 張金屯.數(shù)量生態(tài)學(xué)[M].北京:科學(xué)出版社,2004:20-220.
[20] 張峰,張金屯,上官鐵梁.歷山自然保護(hù)區(qū)豬尾溝森林群落植物多樣性研究[J].植物生態(tài)學(xué)報(bào),2002,26(S1):46-51.
[21] 秦偉,朱清科,張宇清,等.陜北黃土區(qū)生態(tài)修復(fù)過程中植物群落物種多樣性變化[J].應(yīng)用生態(tài)學(xué)報(bào),2009,20(2):403-409.