-
由于邻近植物个体对有限资源的共同需求,群落中植物种间存在相互竞争,从而使竞争双方或某一方植物的生长和发育受到阻碍[1],而植物间的竞争作用会对森林群落物种的天然更新、分布格局、多样性、结构等方面产生影响[2-3],进而影响群落的生产力和稳定性。生态位能有效反映物种在其所处的生态系统中的地位关系,可以将种间关系定量化描述和研究,其理论成为解释群落中物种共存与竞争机制的基本理论[4-6],而物种多样性和群落结构可以由物种生态位宽度和种间生态位重叠决定,能定量反映物种对资源的利用能力、在群落或生态系统中的功能位置及其所在群落的稳定性[7-9]。所以,基于生态位理论研究群落内物种间的相互关系,对深入揭示群落物种分布、共存及群落稳定性机制至关重要。
油松(Pinus tabuliformis Carr.)为我国特有树种,是华北植物区系的代表种之一,其根系发达、枝叶繁茂、具有较强的抗旱性和天然更新能力[10-11],对区域水土保持起到非常重要的作用[12]。本研究通过对晋南油松人工林林下植被的调查发现,油松天然更新存在障碍。牛健[13]通过综合前人研究发现,球果和种子质量、枯落物、坡向、坡位、郁闭度、干扰等是导致油松天然更新障碍的主要原因;但近年来,中条山林业主管部门对其辖区内不同类型油松人工林采取了人工促进天然更新(图1),如清除枯枝落叶、间伐等措施,发现当年春季油松萌发数量较多,但第2年存活较少。所以,种实、枯落物、非生物环境因素并不是影响该地区油松更新障碍的主要原因,基于此,本研究提出以下假设:物种之间的竞争可能是影响该地区油松天然更新的主要原因。为了证明该假设,本研究基于中条山油松人工林群落样地调查数据,通过分析油松人工林下主要植物的种间关联性、生态位重叠及群落稳定性,探讨种间关联和生态位理论在油松人工林天然更新障碍及群落稳定性形成机制中的作用,为油松人工林经营管理提供理论依据。
-
本研究位于山西省中条山十河林场(35°31′~35°51′ N,111°50′~112°00′ E),属太行山向南延伸的余脉,处于黄土高原向华北平原过渡的东南边缘区,土壤类型随海拔由低到高分别为山地褐土、淋溶褐土、棕壤、亚高山草甸土。该区属于暖温带半湿润大陆季风气候,年均气温10.3~14.3℃,≥10℃的年积温为3 400~4 700℃,年均降水量为524.3~654.3 mm,无霜期199~229 d[14-16]。植被属暖温带落叶阔叶林带,主要树种有油松、华山松(P. armandii Franch.)、辽东栎(Quercus wutaishanica Mary)、橿子栎(Q. baronii Skan)、侧柏(Platycladus orientalis Linn.)、鹅耳枥(Carpinus turczaninowii Hance)、栓皮栎(Q. variabilis Bl.)和白桦(Betula platyphylla Suk.)等,林下植被主要由辽东栎、华山松、油松、连翘(Forsythia suspensa Thunb.)、土庄绣线菊(Spiraea pubescens Turcz.)、披针薹草(Carex lancifolia C. B. Clarke)、早开堇菜(Viola prionantha Bunge)、茜草(Rubia cordifolia L.)和金线草(R. membranace Thunb.)等物种组成。
-
在对中条山油松人工林全面踏查的基础上,在油松人工林分布的核心地区进行样地调查,样地概况见表1。样地面积为675 m2或1 200 m2,利用网格法将每个样地分割成3个20 m×20 m或15 m×15 m的样方,对样方内胸径(DBH)≥5 cm的木本植物进行每木检尺,记录每个个体的物种名、胸径、树高等,根据二类调查数据查询林分年龄。在每个20 m×20 m(15 m×15 m)样方内设置16(9)个5 m×5 m和1 m×1 m的灌木和草本样方,记录物种名、多度、株高、盖度等。共计录207种植物,其中,乔木29种、灌木47种、草本112种、藤本19种;同时记录每个样地的经纬度、海拔、郁闭度、坡度、坡向等环境因子。
序号
No.经度(E)
Longitude纬度(N)
Latitude树种组成
Species composition林龄
Age/a海拔
Elevation/
m样地面积
Area/
m2坡向
Slope aspect坡度
Slope/
(°)郁闭度
Coverage/
%平均胸径
Dq/
cm平均树高
Height/
m密度
Density/
(株∙ hm−2)1 111.92° 35.50° 6 油+4 华 40 1 916 1 200 西南 SW 22 80 24.08 14.89 808 2 111.92° 35.50° 6 油+4 辽 46 1 924 675 东南 SE 34 75 17.56 9.31 1 511 3 111.88° 35.50° 10 油 28 1 880 1 200 西 W 19 60 14.77 12.43 1 600 4 111.94° 35.51° 8 油+2 辽 38 1 750 1 200 西南 SW 36 55 18.66 11.20 925 5 111.93° 35.51° 8 油+2 华 48 1 722 1 200 东 E 21 85 25.64 15.46 508 6 111.90° 35.49° 9 油+1 华 48 1 933 1 200 东北 NE 13 80 27.84 16.31 583 7 111.90° 35.49° 10 油 41 1 917 1 200 西南 SW 20 60 27.79 14.74 500 8 111.94° 35.48° 8 油+2 辽 36 1 509 1 200 南 S 26 80 17.39 11.60 1 383 9 111.92° 35.51° 8 油+1 辽+1 华 39 1 997 1 200 西南 SW 25 75 14.79 9.07 2 475 10 111.89° 35.49° 9 油+1 落 32 1 872 1 200 东南 SE 11 70 17.39 12.95 1 250 11 111.90° 35.49° 10 油 48 1 928 1 200 西 W 8 60 29.32 15.94 467 12 111.90° 35.50° 7 油+2 华+1 辽 46 1 960 1 200 西南 SW 23 75 21.76 10.90 558 Table 1. The information of community plots of P. tabuliformis Carr. plantation
-
基于样地调查数据,选取油松林下植物物种为研究对象,建立“样方—物种多度”矩阵,共计36个样方数据用于种间关联分析。采用方差比率(VR)检验群落物种总体联结性,若VR>1则表示群落物种总体正关联,VR<1则表示负关联。采用W统计量来检验方差比率的显著程度,若W落入95%的置信区间内,则群落物种总体关联不显著,反之则表示显著关联[7]。采用经过Yates[17]矫正的χ2检验物种之间是否存在显著关联;基于联结系数(AC)得到物种之间联结的紧密程度,若AC值越趋向于1.00,表示种间正联结越紧密,AC值越趋向于−1.00,表示种间负联结越紧密,AC=0则表示种间相对独立;采用Pearson相关系数判断种对之间的相关程度[8];采用Levins指数度量物种的生态位宽度(Bi),同时计算种对之间的Levins生态位重叠指数(Oik)[9],并绘制物种生态位重叠指数半矩阵图。由于物种数较多,群落中常见种种间关系更容易监测到,且常见种决定群落物种分布[18-19],所以,选择群落中频度最高的30个常见物种进行分析。
-
采用Godron稳定性测定方法评价群落的稳定性。以36个样方为单位,将群落中不同物种的频度由大到小排序,将物种的频度换算成相对频度,按相对频度由大到小逐步累积相加,同时,一一对应逐步累积总物种数倒数,采用散点平滑曲线模拟两次方程和两个坐标轴100处连接的对角直线来判断累积总物种数倒数与累积相对频度比值的交点坐标,20/80是群落的稳定点,该交点坐标越接近20/80,群落就越稳定[20]。
-
数据整理采用Excel 2016(Microsoft,USA)完成,分析及作图采用R 4.0.3软件中“vegan”、“spaa”程序包完成。
-
中条山油松人工林下最常见的30个物种的总体联结性分析结果表明,30个物种的方差比率为1.98,大于1,表明种间存在正联结,其W统计量为71.22,高于χ2分布临界值(17.71),表明物种间总体上显著正联结。
表2表明:30个物种共形成406个种对,其中,辽东栎因频度太大在每个样方中都有分布而无法分析其与共处群落中其他种之间的关联性,即偏离种。χ2检验结果显示:无联结种对共计372对,占全部种对的91.63%;有联结种对共计34对,占全部种对的8.37%,其中,极显著和显著联结的种对分别为15对和19对,分别占有联结种对的44.12%和55.88%,说明常见种间联结较松散。
种间关联 Interspecific association 类型 Types 数值范围 Value range 种对数 Species pair number 百分比/% χ2 检验
Chi−square test无联结 No association p≥0.05 372 91.63 有联结 Association p<0.01 15 3.69 0.01≤p<0.05 19 4.68 总数 Total 406 100.00 联结系数
Association coefficient正联结 Positive association AC=1 26 6.40 0.60≤AC<1.00 23 5.67 0.20≤AC<0.60 96 23.65 0<AC<0.20 66 16.26 无联结 No association AC=0 1 0.25 负联结 Negative association −0.20<AC<0 26 6.40 −0.60<AC≤−0.20 73 17.98 −1.00<AC≤−0.60 23 5.67 AC=−1 72 17.73 总数 Total 406 100.00 Pearson相关
Pearson correlation正相关 Positive correlation p<0.01 50 11.49 0.01≤p<0.05 25 5.75 p≥0.05 139 31.95 负相关 Negative correlation p<0.01 1 0.23 0.01≤p<0.05 10 2.30 p≥0.05 210 48.28 总数 Total 435 100.00 Table 2. Results of interspecific association analysis of common species
联结系数分析结果(表2)显示:406个种对中(存在1个偏离种),正联结211个种对,负联结194个种对,无联结1个种对,分别占全部种对的51.97%、47.78%、0.25%。正联结显著的种对(AC≥0.20,145对)小于负联结显著的种对(AC≤−0.20,168对),其中,极显著正联结(AC≥0.60)种对共计49对,显著正联结(0.20≤AC<0.60)种对共计96对,极显著负联结(AC≤−0.60)种对共计95对,显著负联结(−0.60<AC≤−0.20)种对73对。此外,AC=1.00和AC=−1.00的种对分别为26对和72对,其中,油松天然更新竞争最激烈的物种主要为土庄绣线菊、五角枫、粗齿铁线莲、青榨槭、金线草、披碱草,均为AC=−1.00。
Pearson相关分析结果(表2)显示:435个种对中显著相关(p<0.05)的种对为86对,无显著相关(p>0.05)的种对为349对,分别占全部种对的19.77%和80.23%,其中,极显著(p<0.01)和显著(0.01≤p<0.05)正相关的种对分别为50对和25对,极显著和显著负相关的种对分别为1对和10对。常见种整体正相关与负相关比值为0.97,其中,显著(p<0.05)正相关与显著负相关比值为6.82。
-
基于Levins指数计算常见种生态位宽度,结果(表3)表明:中条山油松人工林下常见种披针薹草的生态位最宽(23.47),辽东栎次之(20.60),而红椋子的生态位最窄(3.87),三裂绣线菊次之(4.86),且披针薹草的生态位宽度是红椋子的6.06倍,灌木、草本、藤本、乔木更新的平均生态位宽度分别为10.13、10.25、11.97、12.79,油松天然更新的生态位宽度为14.69,相对较高。
物种
Species生态位宽度
Niche width物种
Species生态位宽度
Niche width物种
Species生态位宽度
Niche width披针薹草 Carex lancifolia 23.47 青榨槭 Acer davidii 12.89 高山露珠草 Circaea alpina 8.09 辽东栎 Quercus wutaishansea 20.60 茜草 Rubia cordifolia 11.84 费菜 Sedum aizoon 7.99 华山松 Pinus armandii 18.79 五角枫 Acer pictum subsp. mono 11.69 南蛇藤 Celastrus orbiculatus 7.98 金线草 Rubia membranacea 17.15 土庄绣线菊 Spiraea pubescens 11.55 赤瓟 Thladiantha dubia 7.85 粗齿铁线莲 Clematis argentilucida 15.03 香薷 Elsholtzia ciliata 10.83 青麸杨 Rhus potaninii 7.07 油松 Pinus tabuliformis 14.69 早开堇菜 Viola prionantha 10.72 筋骨草 Ajuga ciliata 5.52 木梨 Pyrus xerophila 14.04 接骨木 Sambucus williamsii 10.41 鹅耳枥 Carpinus turczaninowii 5.31 太平花 Philadelphus pekinensis 14.02 毛丁香 Syringa tomentella 10.06 披碱草 Elymus dahuricus 5.18 榆叶梅 Amygdalus triloba 13.87 灰栒子 Cotoneaster acutifolius 9.82 三裂绣线菊 Spiraea trilobata 4.86 鞘柄菝葜 Smilax stans 13.08 唐古特忍冬 Lonicera tangutica 9.73 红椋子 Swida hemsleyi 3.87 注:表中乔木均为油松人工林下天然更新。
Notes: The trees in the table are natural regeneration under the P. tabuliformis Carr. plantation.Table 3. Calculating the niche width of common species based on Levins index
-
基于Levins生态位重叠指数分析了30个常见种种对间的生态位重叠情况,并绘制半矩阵图(图2)。由图2和表4可知:生态位重叠指数Oik<0.59的种对有338对(77.70%),而Oik≥0.59的种对有97对(22.30%),30个常见种总体生态位重叠较低,但有11个种对生态位重叠指数大于0.98,如筋骨草与辽东栎、披针薹草、青榨槭、木梨、五角枫等物种的生态位重叠均较大,说明部分种之间的生态相似性较大。油松天然更新与费菜、披针薹草、华山松、赤瓟等物种的生态位重叠较大,说明油松更新与以上几种植物资源利用相似或占据了相近的生态空间。
项目 Item Oik≥0.98 0.78≤Oik<0.98 0.59≤Oik<0.78 0.40≤Oik<0.59 0.21≤Oik<0.40 Oik≤0.21 种对数 Species pair number 11 36 50 87 115 136 百分比/% 2.53 8.28 11.49 20.00 26.44 31.26 Table 4. Statistics of Levins niche overlap
-
将油松人工林内所有物种(207)的累积相对频度与所有物种数的累积倒数一一对应,绘制散点图(图3)。根据种总数倒数累计-相对频度累计的散点图,完成平滑曲线的模拟,其Godron稳定性曲线为y=0.009 6x2−0.038 3x+5.006 9(R2=0.993 7),求得交点坐标(61.28,38.72),距离群落稳定点(20,80)相对较远。
Effects of Interspecific Relationship on Natural Regeneration and Community Stability of Pinus tabuliformis Plantation in Zhongtiao Mountain
- Received Date: 2021-03-30
- Accepted Date: 2021-10-14
- Available Online: 2022-02-20
Abstract: