陳玉婷王菁華張 慧燕怡帆朱 雅崔志鵬任學(xué)敏,2
(1. 南陽(yáng)師范學(xué)院生命科學(xué)與農(nóng)業(yè)工程學(xué)院,河南 南陽(yáng) 473061;2. 河南省南水北調(diào)中線水源區(qū)生態(tài)安全重點(diǎn)實(shí)驗(yàn)室,河南 南陽(yáng) 473061)
森林群落的物種組成、多樣性及其與環(huán)境因子的關(guān)系一直是生態(tài)學(xué)研究的熱點(diǎn)問(wèn)題,因?yàn)槠洳粌H是了解生態(tài)系統(tǒng)過(guò)程、形成及其維持機(jī)制的重要途徑,也是森林合理經(jīng)營(yíng)管理的基礎(chǔ)[1-2]。然而眾多研究表明,不同群落其主要影響因子不同[3-4]。李巧玉等[5]的研究表明,川西亞高山彩葉林群落的組成和物種多樣性受海拔和坡向調(diào)控;溫佩穎等[6]的研究顯示,坡度和坡向?qū)﹂熑~紅松(Pinus koraiensis)林物種多樣性存在顯著影響。更多研究表明,群落的物種組成和多樣性受到地形、土壤、群落等多種因素的綜合作用,如Yang等[7]發(fā)現(xiàn),海拔、坡度、坡向、土壤氮和磷含量是影響亞熱帶常綠-落葉闊葉混交林物種組成的主要環(huán)境因子;任學(xué)敏等[8]的研究顯示,巴山冷杉(Abies fargesii)-糙皮樺(Betula utilis)混交林物種分布及多樣性主要受海拔、巖石蓋度和土壤pH 值的影響。因此,在森林經(jīng)營(yíng)管理中,針對(duì)不同群落,確定影響其物種組成和多樣性的主要環(huán)境因子,能夠更有針對(duì)性地制定高效、經(jīng)濟(jì)的管理策略。
栓皮櫟(Quercus variabilis)是我國(guó)森林重要樹(shù)種之一,同時(shí)也是常用造林樹(shù)種。在伏牛山南麓較低海拔保有大面積的純林,對(duì)該區(qū)域水土保持、水源涵養(yǎng)以及經(jīng)濟(jì)發(fā)展具有重要作用,尤其是該區(qū)域作為南水北調(diào)中線工程的核心水源涵養(yǎng)區(qū),其生態(tài)功能更是不容忽視。長(zhǎng)期以來(lái),盡管學(xué)者們對(duì)該群落也進(jìn)行了大量研究,但主要涉及到群落種群空間格局[9]、種群生長(zhǎng)動(dòng)態(tài)的地理格局及驅(qū)動(dòng)因素[10]、種間關(guān)聯(lián)[11]、葉片化學(xué)計(jì)量學(xué)特征[12]、生物量和凈生產(chǎn)力[13]、群落恢復(fù)過(guò)程中植物多樣性變化[14]等,而罕見(jiàn)關(guān)于該群落物種組成和多樣性與環(huán)境因子關(guān)系的研究報(bào)道。
本研究對(duì)伏牛山南麓栓皮櫟群落物種組成、多樣性以及土壤、群落、地形等環(huán)境因子進(jìn)行了調(diào)查,探索群落物種組成和多樣性與環(huán)境因子的關(guān)系,分析影響物種組成和多樣性的主要因子,旨在為該區(qū)域栓皮櫟群落經(jīng)營(yíng)管理措施的制定提供科學(xué)依據(jù)。
研究區(qū)位于河南伏牛山國(guó)家級(jí)自然保護(hù)區(qū)(110°30′ ~ 113°05′E,32°45′ ~ 34°00′N),該區(qū)域?yàn)楸眮啛釒蚺瘻貛У倪^(guò)渡區(qū),屬典型大陸性季風(fēng)氣候。年均溫15.1 ℃,年均≥10 ℃活動(dòng)積溫4 841 ℃,無(wú)霜期224 d,年均降水量846.4 mm。該區(qū)域地形地貌復(fù)雜,包括中山、低山、丘陵和河谷平原等。土壤類型主要有山地棕壤、山地黃棕壤和山地褐土。植被類型豐富,有常綠針葉林、落葉闊葉林和半常綠闊葉林等。栓皮櫟是落葉闊葉林的主要樹(shù)種之一,在海拔600 ~ 1200 m 區(qū)域有廣泛分布,800 ~ 1 000 m 區(qū)域常形成純林,林分郁閉度81%,群落具明顯層次結(jié)構(gòu),喬木層栓皮櫟為優(yōu)勢(shì)種,平均株高10.65 m,平均胸徑10.63 cm,銳齒槲櫟(Quercus alienavar.acutiserrata)、枹櫟(Quercus serrata)等也有少量分布;灌木層主要有紅葉(Cotinus coggygriavar.cinerea)、胡枝子(Lespedeza bicolor)、綠 葉 胡 枝 子(Lespedeza buergeri)、杜鵑(Rhododendron simsii)等;草本層主要有大披針薹草(Carex lanceolata)、大油芒(Spodiopogon sibiricus)和費(fèi)菜(Phedimus aizoon)等。
野外調(diào)查于2019 年8—9 月進(jìn)行。在伏牛山南麓河南伏牛山國(guó)家級(jí)自然保護(hù)區(qū)栓皮櫟群落內(nèi)設(shè)置典型樣地(人為干擾小,生長(zhǎng)良好)30 塊(表1),每個(gè)樣地面積20 m × 20 m,用于喬木調(diào)查,在每個(gè)樣地4 角和中心各設(shè)置1 個(gè)2 m × 2 m的小樣方,用于灌木調(diào)查,在每個(gè)2 m × 2 m 的小樣方內(nèi)再設(shè)置1 個(gè)1 m × 1 m 的小樣方用于草本植物調(diào)查。記錄每塊20 m × 20 m 樣地內(nèi)株高>1.3 m的所有喬木的種名、株高和胸徑,估計(jì)林分郁閉度和各物種蓋度;分別于每個(gè)2 m × 2 m 和1 m ×1 m 小樣方內(nèi)記錄出現(xiàn)的所有灌木和草本植物的種名和數(shù)量,測(cè)量各物種高度,估計(jì)蓋度。于每個(gè)2 m × 2 m 小樣方的中心測(cè)量枯落物厚度,以5 個(gè)小樣方枯落物厚度的均值作為該樣地枯落物厚度。用手持式GPS(佳明,MAP639SC)測(cè)量每個(gè)樣地的海拔和經(jīng)緯度,用地質(zhì)羅盤測(cè)量坡度和坡向。
表 1 樣地基本概況Table 1 Basic overview of the sample plot
在每個(gè)樣地對(duì)角線上3 個(gè)1 m × 1 m 的小樣方中心,用100 cm3環(huán)刀(50.46 mm × 50 mm)分別取0 ~ 10、10 ~ 20 cm 和20 ~ 30 cm 層次的土壤原狀土,用于土壤物理性質(zhì)測(cè)定,以每樣地3 個(gè)樣點(diǎn)3 個(gè)層次的平均值作為該樣地每個(gè)土壤物理指標(biāo)的最終值。在每個(gè)樣地環(huán)刀土的取樣點(diǎn)再分別挖取30 cm 的土壤剖面,從上到下收集0 ~ 30 cm層次的土壤,將3 個(gè)樣點(diǎn)收集的土樣混合均勻,用四分法取1 kg 左右裝入密封袋,帶回實(shí)驗(yàn)室,土樣風(fēng)干后,過(guò)14 目、60 目和100 目土壤篩后進(jìn)行土壤化學(xué)性質(zhì)測(cè)定。
采用環(huán)刀法測(cè)定土壤容重、總孔隙度、毛管孔隙度、非毛管孔隙度、土壤飽和含水量和毛管最大持水量。pH 值(土∶水=1∶5)用pH 計(jì)(pHS-3CT 型,上海)測(cè)定;有機(jī)質(zhì)含量用重鉻酸鉀容量法-稀釋熱法測(cè)定,全氮含量用凱氏定氮法測(cè)定,堿解氮含量用堿解-擴(kuò)散法測(cè)定,有效磷含量用0.03 mol/L NH4F-0.025 mol/L HCl 浸提法測(cè)定,全磷含量用酸溶(H2SO4-HClO4)-鉬銻抗比色法測(cè)定,全鉀含量用NaOH 熔融-火焰光度法測(cè)定,速效鉀含量用NH4OAc 浸提-火焰光度法測(cè)定。
2.4.1 物種重要值
本研究對(duì)栓皮櫟群落的調(diào)查涉及到喬木層、灌木層和草本層,為了防止群落分層統(tǒng)計(jì)及計(jì)算公式不同造成的層間物種無(wú)可比性問(wèn)題[15],重要值(IV)計(jì)算統(tǒng)一采用以下公式[16]:
2.4.2 排序
應(yīng)用排序法探索環(huán)境因子對(duì)物種組成的影響。預(yù)先的除趨勢(shì)對(duì)應(yīng)分析(DCA)顯示,最長(zhǎng)梯度長(zhǎng)度為1.85,因此選用線性模型中的冗余分析(RDA)。排序時(shí)應(yīng)用前選法剔除P>0.05 的環(huán)境因子,最終21 個(gè)環(huán)境因子中的6 個(gè)(海拔、坡向、灌木蓋度、土壤非毛管孔隙度和堿解氮含量)進(jìn)入最后的分析。由于坡向數(shù)據(jù)為0° ~ 360°的羅盤測(cè)量值,分析時(shí)將其轉(zhuǎn)換為0 ~ 1 的數(shù)據(jù)。轉(zhuǎn)換后的數(shù)值越大說(shuō)明生境條件越熱干,越小說(shuō)明生境條件越冷濕。轉(zhuǎn)換公式[17]為:
式中:TRASP 為坡向指數(shù),A為坡向方向角度。
為了評(píng)估不同環(huán)境因子對(duì)物種組成影響的相對(duì)重要性,按照Lososová等[18]的方法,分別采用RDA 和偏RDA 排序法計(jì)算了每一環(huán)境因子對(duì)物種組成影響的總效應(yīng)和凈效應(yīng)。采用Monte Carlo 隨機(jī)置換(9 999 次)檢驗(yàn)對(duì)所有排序軸進(jìn)行顯著性檢驗(yàn)。
2.4.3 物種多樣性
選用物種豐富度指數(shù)(R)、Simpson 指數(shù)(D),Shannon-Wiener 指數(shù)(H')和Pielou 均勻度指數(shù)(J)來(lái)反映栓皮櫟群落的物種多樣性,計(jì)算公式[19]為:
式中:S為每樣地中物種總種類數(shù),Pi為樣地中物種i的相對(duì)重要值。
應(yīng)用多元線性逐步回歸分析影響物種多樣性的主要環(huán)境因子。采用Canoco for Windows5.0 軟件進(jìn)行RDA 排序和繪圖,其他統(tǒng)計(jì)分析均由IBM SPSS Statistics 24 完成。
RDA 排序結(jié)果表明,6 個(gè)環(huán)境因子共解釋了栓皮櫟群落物種組成變異的45.38%,所有典范軸顯著性檢驗(yàn)結(jié)果的極顯著(P<0.01)性說(shuō)明物種組成與環(huán)境因子之間關(guān)系密切。RDA1 軸和2 軸分別解釋了物種組成變異的34.24%和7.04%,累積解釋量占所有環(huán)境因子解釋量的90.97%,表明RDA 前2 軸在解釋物種組成變異中起主要作用(表2)。由環(huán)境因子與RDA 排序軸的相關(guān)系數(shù)(表3)可知,RDA1 軸與6 個(gè)環(huán)境因子相關(guān)性均顯著(P<0.05),表明其是一個(gè)綜合軸,但主要反映坡向和土壤堿解氮含量梯度;RDA2 軸除與堿解氮含量相關(guān)性不顯著外,與其他5 個(gè)環(huán)境因子相關(guān)均顯著,綜合反映了海拔、灌木蓋度、土壤pH 和非毛管孔隙度梯度。
表 2 栓皮櫟群落物種RDA 排序概要Table 2 Summary of RDA ordination of species inQ. variabiliscommunity
在調(diào)查的30 個(gè)樣地中,共記錄了106 種植物,隸屬于45 科85 屬,其中喬木26 種,灌木34 種,草本46 種。
物種與環(huán)境因子的RDA 排序見(jiàn)圖1。
圖 1 栓皮櫟群落植物物種與環(huán)境因子的RDA 排序圖Fig. 1 RDA plot of species and environmental factors inQ. variabiliscommunity
由圖1可知,物種組成隨環(huán)境因子表現(xiàn)出明顯的規(guī)律,栓皮櫟位于RDA1 軸左端,較干熱坡向和較低土壤堿解氮含量是其分布的適宜生境;油松(Pinus tabuliformis)、費(fèi)菜和火絨草(Leontopodium leontopodioides)等位于排序圖的左下方,其相應(yīng)的生境條件為較低海拔、土壤具有較小的非毛管孔隙度、較低酸度和堿解氮含量;位于排序圖右下方的物種,如枹櫟、紅葉、青榨槭(Acer davidii)和杜鵑等,其喜好分布的生境為較冷濕且灌木蓋度較高的坡向;大葉華北繡線菊(Spiraea fritschianavar.angulata)、 山 櫻 花(Prunus serrulata)、多腺懸鉤子(Rubus phoenicolasius)、漆(Toxi codendronvernicifluum)和絲葉薹草(Carex capilliformis)等位于RDA2 軸的上端,表明其在海拔較高、灌木蓋度較低、土壤非毛管孔隙度較大且較酸的生境有最多分布。
Monte Carlo 檢驗(yàn)顯示,6 個(gè)環(huán)境因子對(duì)栓皮櫟群落物種組成的總效應(yīng)均達(dá)顯著水平(P<0.05),其中坡向的解釋量最高,為18.41%,土壤堿解氮含量、pH、灌木蓋度和海拔的解釋量依次降低,土壤非毛管孔隙度的解釋量最?。粌粜?yīng)達(dá)到顯著水平(P<0.05)的有2 個(gè)環(huán)境因子,分別為坡向和土壤堿解氮含量,解釋量分別為16.50%和14.48%(表4)。
表 4 環(huán)境因子對(duì)栓皮櫟群落植物組成的影響Table 4 Effect of environmental factors on species composition in Q. variabiliscommunity
由表5 可知,栓皮櫟群落不同樣地間的物種豐富度存在較大差異,物種豐富度最大的樣地記錄到26 種物種,最小的記錄到8 種物種,平均15.9 種;但從Simpson 指數(shù)、Shannon-Wiener 指數(shù)和Pielou 均勻度指數(shù)來(lái)看,3 個(gè)指數(shù)的標(biāo)準(zhǔn)差均較小,分別為0.06、0.27 和0.06,說(shuō)明栓皮櫟群落物種分布較均勻,即物種豐富度較低的樣地其物種個(gè)體數(shù)較多。多元線性逐步回歸分析結(jié)果表明,物種豐富度主要受枯落物厚度的影響,Simpson 優(yōu)勢(shì)度指數(shù)和Shannon-Wiener 指數(shù)主要受坡向和堿解氮含量的影響,而影響Pielou 均勻度指數(shù)的環(huán)境因子主要是灌木蓋度、土壤pH 和堿解氮含量(表6)。說(shuō)明栓皮櫟群落物種多樣性是受枯落物厚度、坡向、土壤堿解氮含量、土壤pH 和灌木蓋度共同作用的結(jié)果。
表 5 栓皮櫟群落物種多樣性指數(shù)Table 5 Species diversity index inQ.variabiliscommunity
4.1.1 地形因子對(duì)栓皮櫟群落物種組成和多樣性的影響
以往眾多研究表明,地形因子能夠顯著地影響物種組成和多樣性[7,20-21],這種結(jié)論在本研究中再一次得到證實(shí)。本研究發(fā)現(xiàn),坡向?qū)λㄆ等郝湮锓N組成存在強(qiáng)烈影響,其總效應(yīng)和凈效應(yīng)分別為18.41%和16.50%,均達(dá)顯著水平(P<0.05);盡管海拔對(duì)物種組成的凈效應(yīng)不顯著,但其卻具有顯著的總效應(yīng)(10.38%)。多元線性逐步回歸分析表明,較干熱坡向物種多樣性較低,反之則高。一些研究對(duì)這些現(xiàn)象的解釋是,地形因子是一種間接影響因素,其主要通過(guò)對(duì)降水、光照、溫度等的再分配影響物種的組成和多樣性[6,20],本研究雖然未對(duì)地形因子與降水、光照、溫度的關(guān)系進(jìn)行直接檢驗(yàn),但對(duì)坡向與一些未進(jìn)入分析的環(huán)境因子進(jìn)行相關(guān)分析(皮爾遜,雙尾)時(shí)發(fā)現(xiàn),與土壤水分關(guān)系密切的土壤總孔隙度(r=-0.44,P<0.05)、毛管孔隙度(r= -0.23,P>0.05)、非毛管孔隙度(r= -0.32,P>0.05)、毛管最大持水 量(r= -0.34,P>0.05)和飽和含水量(r=-0.43,P<0.05)均與坡向存在負(fù)相關(guān)關(guān)系,說(shuō)明本研究坡向可能主要通過(guò)影響土壤濕度來(lái)影響物種的組成和多樣性,土壤濕度較高的坡向能為更多物種提供所需水分。
物種組成與地形因子的關(guān)聯(lián)具有重要應(yīng)用價(jià)值,因?yàn)樵谥脖换謴?fù)、營(yíng)建或改造森林群落時(shí),可根據(jù)這種關(guān)聯(lián)在不同地形條件下選擇不同的物種組合,從而更利于群落的構(gòu)建和發(fā)展[22]。本研究中,栓皮櫟總是與熱干坡向相關(guān)聯(lián),而淫羊藿(Epimedium brevicornu)、二翅六道木(Abelia macrotera)、杜鵑和青榨槭等物種總是出現(xiàn)于冷濕坡向,這些結(jié)果能為栓皮櫟群落的營(yíng)造和管理提供理論指導(dǎo)。
4.1.2 群落因子對(duì)栓皮櫟群落物種組成和多樣性的影響
群落物種組成和多樣性與植被覆蓋度通常存在密切關(guān)系。Boscutti 等[23]研究發(fā)現(xiàn),灌木蓋度顯著影響了物種組成,同時(shí)也發(fā)現(xiàn)灌木蓋度越高物種豐富度越低;萬(wàn)猛等[24]的研究表明,太行山南麓栓皮櫟群落的物種豐富度和多樣性指數(shù)主要受灌木層和草本層影響,本研究結(jié)果部分與之相同。本研究發(fā)現(xiàn)灌木蓋度對(duì)物種組成具有顯著的總效應(yīng),且物種多樣性與灌木蓋度正相關(guān),其中對(duì)物種豐富度、Simpson 優(yōu)勢(shì)度指數(shù)和Shannon-Wiener 指數(shù)的影響不顯著,對(duì)Pielou均勻度指數(shù)的影響達(dá)顯著(P<0.05)水平。出現(xiàn)這種現(xiàn)象的原因是,栓皮櫟群落生長(zhǎng)的生境相對(duì)熱干,喬木層林分郁閉度(0.81 ± 0.07)和草本層蓋度((3.15± 2.42)%)樣地間變化不大,但林下灌木蓋度具有較大的變異,最大者達(dá)45.00%,最小者僅1.20%,較高的灌木蓋度能夠提供相對(duì)蔭蔽、濕潤(rùn)的環(huán)境,為較多耐蔭物種的生長(zhǎng)提供良好條件,從而導(dǎo)致物種組成和多樣性的變異[25]。
4.1.3 土壤因子對(duì)栓皮櫟群落物種組成和多樣性的影響
土壤作為植物生長(zhǎng)的基質(zhì),其理化性質(zhì)必然會(huì)對(duì)植物組成和多樣性產(chǎn)生影響,但群落不同其主要影響因子也不同,如土壤pH 是影響巴山冷杉-糙皮樺混交林物種組成和多樣性的主要土壤因子[8],而亞熱帶常綠闊葉林物種組成則主要受土壤全磷和速效鉀含量的影響[26]。本研究發(fā)現(xiàn),影響栓皮櫟群落物種組成的主要土壤因子為堿解氮含量,其總效應(yīng)和凈效應(yīng)均達(dá)顯著(P<0.05)水平,同時(shí)也發(fā)現(xiàn)土壤堿解氮含量對(duì)Simpson 優(yōu)勢(shì)度指數(shù)、Shannon-Wiener 指數(shù)和Pielou 均勻度指數(shù)有顯著(P<0.05)的正影響,而影響物種豐富度的主要是枯落物厚度。不同群落物種組成和多樣性對(duì)土壤因子響應(yīng)的差異可能源于群落物種組成的不同,對(duì)資源需求、利用能力及方式存在的差異[26-27]。
一些研究認(rèn)為,群落內(nèi)較厚的枯落物不利于物種的生長(zhǎng)發(fā)育,但本研究發(fā)現(xiàn),栓皮櫟群落枯落物厚度對(duì)物種豐富度存在顯著的正影響,而對(duì)Simpson 優(yōu) 勢(shì) 度 指 數(shù)、Shannon-Wiener 指 數(shù)和Pielou 均勻度指數(shù)的影響雖然不顯著,但均存在正相關(guān)關(guān)系(相關(guān)系數(shù)分別為0.25、0.35 和0.12),說(shuō)明枯落物對(duì)栓皮櫟群落物種多樣性有一定的促進(jìn)作用,相似結(jié)果也出現(xiàn)于對(duì)銳齒槲櫟群落的研究報(bào)道中[28]。其原因一方面可能是栓皮櫟群落的枯落物厚度相對(duì)較?。?2.83 ± 0.77) cm),不會(huì)對(duì)物種產(chǎn)生明顯的抑制;另一方面在植物生長(zhǎng)可忍耐的枯落物厚度范圍內(nèi),其厚度的增加能夠明顯提高土壤水分含量,為更多物種提供生長(zhǎng)資源[29]。
本研究發(fā)現(xiàn),土壤堿解氮含量對(duì)栓皮櫟群落物種多樣性存在顯著的正影響,這是因?yàn)楸狙芯恐兴ㄆ等郝涞耐寥缐A解氮含量相對(duì)缺乏,僅為(66.86 ± 33.11) mg/kg,處于全國(guó)第二次土壤普查養(yǎng)分6 級(jí)分級(jí)標(biāo)準(zhǔn)中的4 級(jí)(中下,60 ~ 89 mg/kg)水平,因此,在一定范圍內(nèi),物種多樣性隨著土壤氮營(yíng)養(yǎng)含量的增加而增加。由此可見(jiàn),在管理中,適量施加氮肥可能有助于提高該群落的物種多樣性。然而,本研究也發(fā)現(xiàn),土壤堿解氮含量對(duì)作為建群種栓皮櫟的分布存在顯著(P<0.05)的負(fù)影響,因此,如果要通過(guò)施加氮肥提高群落物種多樣性,應(yīng)結(jié)合坡向這一對(duì)栓皮櫟分布有顯著影響的因子,避開(kāi)其分布較多的熱干坡向,施用于其分布較少的冷濕坡向,這樣既最大程度地保障了栓皮櫟的生長(zhǎng),也促進(jìn)了群落物種多樣性的提升。至于適宜的氮肥施用量范圍,還需要進(jìn)一步地試驗(yàn)研究。
1)環(huán)境因子能夠顯著地解釋伏牛山南麓河南伏牛山國(guó)家級(jí)自然保護(hù)區(qū)栓皮櫟群落物種組成,海拔、坡向、灌木蓋度、土壤非毛管孔隙度、pH 和堿解氮含量對(duì)物種組成的總效應(yīng)均達(dá)顯著水平,其中坡向和土壤堿解氮含量的凈效應(yīng)也達(dá)顯著水平,是最主要的影響因子。
2)栓皮櫟群落物種多樣性主要受坡向、枯落物厚度、灌木蓋度、土壤pH 和堿解氮含量的影響,較冷濕坡向、枯落物較厚、灌木蓋度較高、土壤酸度較大且堿解氮含量豐富的生境具有較高的物種多樣性。
在栓皮櫟群落經(jīng)營(yíng)管理中,應(yīng)綜合考慮環(huán)境因子對(duì)物種組成和多樣性的影響,如在較冷濕坡向適量施加氮肥,既可提高群落的物種多樣性,又可最大程度地減小對(duì)建群種栓皮櫟的影響。