张 喜 霍 达 张佐玉 王莉莉 崔迎春
(1 贵州省林业科学研究院 贵阳 550011;2 贵州财经大学 贵阳 550005)
赤水河下游紫色土不同形态钾含量与地形要素及毛竹林植物多样性指数的关系
张 喜1霍 达1张佐玉1王莉莉2崔迎春1
(1 贵州省林业科学研究院 贵阳 550011;2 贵州财经大学 贵阳 550005)
研究了赤水河下游紫色土不同形态钾含量与地形要素及毛竹林植物多样性指数的相关规律。结果表明:1)紫色土不同形态钾含量在毛竹林不同类型及土层间差异不显著,随土层厚度增加,全钾含量增加、速效钾含量降低。全钾含量与土壤主要物理指标相关性不显著,速效钾含量与A(B)层土壤密度及孔隙度相关性显著,不同形态钾含量与不同土层土壤阳离子交换量、以及表层土壤有机质含量、全氮含量及有效磷含量相关性显著。2)毛竹林乔木层植物Margalef指数、Simpson指数、Shannon-Wiener指数和Pielou指数值随毛竹重要值降低而增加,灌木及草本层植物多样性指数值变化的趋势性各异。不同层内植物多样性指数相关性显著,不同层间植物Margalef指数、乔木和灌木层植物Simpson指数相关性显著。3)不同土层土壤全钾含量与乔木及草本植物多样性因子相关性显著,速效钾含量与地形及草本植物多样性因子相关性显著。土壤全钾含量与毛竹林乔木、灌木及草本层植物多样性指数间相关显著的直线和二次多项式指标对数量分别为50.00%、50.00%,土壤速效钾含量的相应值为13.89%、44.44%,表明土壤不同形态钾含量与植物多样性指数关系存在无关型、直线型和曲线型3种类型,印证了毛竹林植物多样性的土壤生态功能多样性。
毛竹林;紫色土;土壤全钾含量;土壤速效钾含量;植物多样性指数
毛竹(Phyllostachyspubescens)广泛分布于秦岭、汉水流域至长江以南各省,是我国栽培历史悠久、面积最大、经济价值最高的竹种之一[1]。位于赤水河下游的贵州省赤水市及相邻的习水县属毛竹中心分布区[2],境内特殊的地质地貌及紫色土壤和丰富的水热资源有利于毛竹生长[3],分布面积达2.68 万hm2(2006年),占贵州省毛竹林面积的58.51%,其中赤水市被誉为“中国十大竹子之乡”[4]。毛竹林科学经营对该区域的社会经济发展和长江上游生态环境保护有重要意义。
对已有文献[5-6]分析发现,该区域内野竹坪林区紫色土毛竹林经营30年后,平均胸径由9.9~13.6 cm降至7.0~10.8 cm,这一现象也出现在山地黄壤毛竹林[7]中,表明毛竹纯林长期经营可致林分质量下降的现象具有普遍性。另一方面,石灰土壤养分指标与淡竹林[8]以及紫色土水源涵养指标和毛竹林[9]的植物多样性指数间有一定关联性,同其他植被类型[10-11]的相似研究结论相近。采用乔木层物种调控与林下套种措施可提高红壤毛竹林林地自肥能力[12],表明通过植被层植物多样性能调控土壤肥力。紫色砂页岩易于崩解、矿质养分丰富,发育的紫色土较一般土壤富含钙、磷和钾元素,肥力较高。分布在赤水河下游的酸性紫色土[3]除具有紫色土的一般特性外,其有机质和全氮含量相对较高,磷和钾含量相对较低。土壤不同形态钾含量与地形要素及毛竹林植物多样性指数的关系如何[13-14]?调控毛竹林植物多样性指数能否成为紫色土不同形态钾含量的改善途径?这些问题均需要科学的试验研究。
国内学者对紫色土毛竹林不同径级类型的土壤理化指标[15]和林分结构指标[16]变化规律,以及水源涵养指标和毛竹林植物多样性指数的关系[9]都有过研究,但未见土壤不同形态钾含量与地形要素及毛竹林植物多样性指数相关性的报道。本文主要研究赤水河下游紫色土不同形态钾含量与地形要素及毛竹林植物多样性指数的关系,揭示其相互作用机理,旨在为毛竹林的可持续经营及通过植物多样性途径调控林地土壤不同形态钾含量提供理论参考。
1.1 研究区概况
研究区位于赤水河下游的赤水市和习水县交界区域。两县(市)位于贵州省西北部,地处东经105°36′25″—106°44′30″、北纬28°06′35″—28°50′15″,海拔变幅221~1 871.9 m。区域地质属杨子区四川盆地分区,主要地层有侏罗系、白垩系和第四系,成土母岩有石英砂岩、紫色砂(页)岩和泥岩等,土壤主要有紫色土、黄壤和黄棕壤等。该区域处于娄山山脉由东南至西北向的延伸带,属四川盆地边缘河谷和中山峡谷地貌。多年平均温度(海拔293 m)18.1 ℃、变幅17.5~19.1 ℃;1月平均温度7.9 ℃、变幅5.8~10.8 ℃;7月平均温度28 ℃、变幅23.9~33.0 ℃;≥0 ℃年积温计365 d和6 621.7 ℃、≥10 ℃年积温计279 d和5 888.3 ℃。年降水量1 268.8 mm,年蒸发量1 307.1 mm。年平均日照时数1 297.7 h、变幅999.4~1 529.2 h。立体气候明显,区域差异显著。主要森林类型包括中亚热带常绿阔叶林、针阔混交林、针叶林和竹林等。
1.2 样地设置与植被调查
在研究区内,利用贵州省森林资源二类调查小班数据库(2006年),筛选紫色土的毛竹林小班建立新的数据库。依据毛竹林小班分布图的集散状态和交通便利程度,实地考查、建立毛竹林调查样地43个,其中赤水市33个、习水县10个。调查样地间海拔高度变幅500~1 183 m,坡度变幅0~45°,土壤A层厚度变幅3~40 cm、A(+B)层厚度变幅17~85 cm;乔木层毛竹平均胸径变幅6.21~12.02 cm、平均高度变幅9.9~17.6 m、平均密度变幅850~4 650株/hm2;杂木的相应值为0~38.42 cm、0~16.3 m、0~1 850株/hm2;灌木层植物平均地径变幅0.15~1.38 cm、平均高度变幅0.13~1.43 m、平均密度变幅0.20万~8.51万株/hm2,草本层植物平均高度变幅0.16~1.29 m、平均密度变幅2.7万~98.0万株/hm2。
调查样地面积20 m × 20 m,分别乔木、灌木及草本3层调查。乔木层样方面积10 m×10 m,计4个;灌木层与草本层样方面积分别为5 m×5 m、1 m×1 m,选择各层片植物发育较好、一般和较差的3种类型各设置1个。记录样地空间位置,记载乔木层物种、胸径、株数、高度,灌木层物种、地径、株数、平均高度,草本层物种、株数、平均高度。在灌木层与草本层植物调查样方内,选择有代表性地段挖掘土壤剖面,按土壤发生层记录剖面特征,分层提取环刀,密封后带回室内测定土壤主要物理指标;分层取土,去除石砾和根系后,等量混合,保留质量约1 kg,带回室内风干,分析土壤主要化学指标。植被调查与土壤取样同时进行,在2008年8月一次完成。
1.3 土壤理化指标测定
土壤理化指标分析参见《森林土壤分析方法》[17]。其中土壤密度、孔隙度和持水量指标测定采用环刀法,pH值采用电位法,有机质含量采用重铬酸钾氧化-外加热法,全氮含量采用半微量凯氏扩散法,水解氮含量采用碱解-扩散法,全磷含量采用碱溶-钼锑抗比色法,有效磷含量采用盐酸-硫酸浸提法,全钾含量采用碱溶-火焰光度法,速效钾含量采用乙酸铵浸提-火焰光度法,阳离子交换量采用乙酸铵交换法,交换性盐基量采用乙酸铵交换-中和滴定法。
1.4 数据分析
1.4.1 植物多样性指数[18]
其中:RD为相对优势度,RF为相对频度,RT为相对密度。
其中:S为物种数,Ni和N分别为某物种和总物种株数,Pi=Ni/N。
1.4.2 数据统计
数据分析运用SPSS17.0和EXCEL 2003软件。
2.1 毛竹林类型及植物多样性指数变化
将毛竹林调查样地乔木层物种重要值(IV)采用离差平方和法及欧氏距离平方聚类划分成5个类型。Ⅰ为毛竹(IV=0.921 3,下同。)-丝栎栲(Castanopsifargesii,0.029 0)-杉木(Cunninghamialanceolata,0.012 9)林,乔木层组成物种20种,包括20个样地;Ⅱ为毛竹(0.389 6)-杉木(0.096 1)-白栎(Quercusfabri,0.091 8)林,乔木层组成物种20种,包括2个样地;Ⅲ为毛竹(0.545 4)-丝栎栲(0.286 0)-西南木荷(Schimawallichii,0.078 0)林,乔木层组成物种9种,包括4个样地;Ⅳ为毛竹(0.700 5)-杉木(0.275 3)-棕榈(Trachycarpusfortunei,0.015 8)林,乔木层组成物种4种,包括5个样地;Ⅴ为毛竹(0.710 7)-丝栎栲(0.079 8)-木荷(Schimasuperba,0.024 3)林,乔木层组成物种33种,包括12个样地。乔木层植物中毛竹重要值大小呈Ⅰ类>Ⅴ类>Ⅳ类>Ⅲ类>Ⅱ类的变化趋势。
毛竹林不同类型间乔木层植物Margalef指数、Shannon-Wiener指数和Pielou指数值差异显著(P<0.05),Simpson指数值在部分类型间差异显著(P<0.05);在灌木层中,类型Ⅰ和Ⅴ的植物Pielou指数值差异显著(P<0.05),其他毛竹林类型的灌木层植物多样性指数值差异不显著(P>0.05);草本层不同植物多样性指数值差异不显著(P>0.05)(表1)。毛竹林不同类型乔木层植物多样性指数值较低,灌木层植物Margalef指数和Simpson指数值较高,草本层植物Shannon-Wiener指数和Pielou指数值较高。
表1 不同类型毛竹林植被层植物多样性指标
注:Mar.,Sim.,Sha.和Pie.分别代表Margalef指数、Simpson指数、Shannon-Wiener指数和Pielou指数,下同。表内数值后的字母代表相应层植物多样性指数在毛竹林不同类型间S检验差异显著(P<0.05)的符号,空格者表示差异不显著。Ⅰ、Ⅱ、Ⅲ、Ⅳ和Ⅴ类型的S检验差异显著符号为a、b、c、d和e。样本量=43。
毛竹林不同层内部植物Margalef指数、Simpson指数、Shannon-Wiener指数和Pielou指数值相关显著(P<0.05)、趋势性各异,其中乔木层植物各指数间为正相关,灌木及草本层植物除Margalef指数和Simpson指数、Shannon-Wiener指数和Pielou指数值为正相关外,其他指数值间均为负相关(表2)。乔木层和灌木层及草本层植物Margalef指数、以及草本层和乔木层及灌木层植物Simpson指数值正相关显著(P<0.05),其他植物多样性指数值间相关不显著(P>0.05)、趋势性各异。
表2 毛竹林植被层植物多样性指数的相关性
注:表内数值右上角*P<0.05、**P<0.01,右上角空格者表示相关性不显著。下同。样本量=43。
2.2 毛竹林土壤不同形态钾含量变化以及与土壤理化指标的相关性
S检验表明,毛竹林不同类型间相同土层的土壤全钾含量和速效钾含量差异不显著(P>0.05),相同类型不同土层的土壤全钾含量和速效钾含量差异也不显著(P>0.05)(表3)。土壤全钾含量呈A层C层、Ⅱ、Ⅲ及Ⅳ类型呈B层
表3 不同类型毛竹林土壤钾含量
注:样本量在A、B和C层的相应值分别为43、41和39。
分析发现,在土壤A层(P<0.01)、B层(P<0.01)及C层(P<0.05)全钾含量和速效钾含量的相关性显著,二者相互影响。毛竹林不同土层土壤密度与全钾含量及速效钾含量呈负相关;土壤最大持水量、毛管持水量、田间持水量、非毛管孔隙度、毛管孔隙度及总孔隙度与全钾含量及速效钾含量呈正相关,其中A层及B层土壤速效钾含量与土壤密度(P<0.05)、最大持水量(P<0.05)、毛管持水量(P<0.01)、田间持水量(P<0.01)及总孔隙度(P<0.05)的相关性显著,B层土壤速效钾含量受土壤物理指标的影响大于A层(表4)。不同土层土壤有机质含量、全氮含量、水解氮含量及阳离子交换量与全钾含量及速效钾含量呈正相关;全磷含量、有效磷含量及pH值与全钾含量及速效钾含量呈负相关;交换性盐基量与全钾含量呈负相关、与速效钾含量呈正相关。其中不同层次土壤阳离子交换量与全钾含量(P<0.01)及速效钾含量(P<0.01)相关性显著,全钾含量与全磷含量(P<0.01)负相关性显著;A层土壤全钾含量与有机质含量(P<0.05)、全氮含量(P<0.05)及有效磷含量(P<0.01)相关性显著,速效钾含量与有机质含量(P<0.01)、全氮含量(P<0.01)、水解氮含量(P<0.05)、有效磷含量(P<0.01)及pH值(P<0.05)相关性显著,速效钾含量受表层土壤化学指标的影响大于全钾含量。
2.3 土壤不同形态钾含量与地形要素及毛竹林植物多样性指数的因子分析
选择地形要素中的海拔高度、坡向、坡度和土壤(A+B)层厚度,乔木层、灌木层及草本层植物Margalef指数、Simpson指数、Shannon-Wiener指数及Pielou指数进行因子分析(BTS=684.209、P<0.000,KMO=0.716)。第1主分量由乔木层植物多样性指数构成,方差贡献率29.92%;第2主分量由灌木层植物多样性指数构成,方差贡献率18.42%;第3主分量由草本层植物多样性指数构成,方差贡献率15.22%;第4主分量由地形要素构成,方差贡献率9.06%;第5主分量由地形要素和植物多样性指数构成,方差贡献率8.17%,5个主分量总贡献率达80.80%。相关分析表明(表5),不同土层土壤全钾含量与乔木植物多样性因子(P<0.05)及草本植物多样性因子(P<0.01)相关性显著,速效钾含量与地形因子(P<0.05)、以及A层土壤速效钾含量和草本植物多样性因子(P<0.05)相关性显著。
表4 毛竹林土壤不同形态钾含量与土壤主要理化指标的相关性
注:物理指标样本量在A、B和C层的相应值为32、30和29,化学指标样本量同表3。
表5 毛竹林土壤不同形态钾含量与地形要素及植物多样性指数主分量的相关性
注:不同土层样本量同表3。
2.4 土壤不同形态钾含量与毛竹林植物多样性指数值的相关性
毛竹林不同层次土壤全钾含量与乔木层植物Margalef指数、Simpson指数及Shannon-Wiener指数、灌木层植物Margalef指数、草本层植物Margalef指数及Simpson指数的直线相关性显著(P<0.05)(表6),C层土壤速效钾含量与乔木层植物Simpson指数、A层及B层土壤速效钾含量与草本层植物Margalef指数及Simpson指数的直线相关性显著(P<0.05)。不同层次土壤全钾含量及速效钾含量与毛竹林植物多样性指数的二次多项式相关系数≥直线相关系数;其中直线相关性不显著、二次多项式相关显著的有:A层土壤速效钾含量与乔木层植物Margalef指数及Simpson指数、草本层植物Shannon-Wiener指数及Pielou指数(P<0.05),B层土壤速效钾含量与乔木层植物Simpson指数、灌木层植物Shannon-Wiener指数及Pielou指数(P<0.05),C层土壤速效钾含量与乔木层植物Margalef指数、灌木层植物Simpson指数、Shannon-Wiener指数及Pielou指数(P<0.05)。表明二次多项式更能反映土壤不同形态钾含量与毛竹林植物多样性的关系,植物多样性指数对土壤全钾含量的影响高于速效钾含量,呈乔木层>草本层>灌木层的变化趋势。
表6 毛竹林土壤不同形态钾含量与植被层植物多样性指数的相关性
注:表中数据分子和分母表示直线和二次多项式的相关系数。
3.1 土壤钾的赋存形式与影响因素
土壤钾主要以离子、化合物和矿物态存在,速效钾含量包括溶于土壤溶液的液态钾和部分土壤胶体呈物理化学吸附的离子态钾,植物残体、根系和微生物钾有效性较高[14]。山地黄壤灌草坡—马尾松林—常绿阔叶林[13]和毛竹林[7]的土壤全钾及速效钾含量均低于本研究紫色土的相应值,尤以速效钾含量差异较大。紫色土不同土层内全钾含量与速效钾含量相关性显著,二者相互影响;毛竹林不同类型及相同类型不同土层间的相应值差异不显著,随土层厚度增加,全钾含量呈增加趋势、速效钾含量呈降低趋势[15]。相关性显著的指标对数量在土壤A层、B层及C层的分布中,全钾含量为33.33%、20.00%、13.33%,速效钾含量为80.00%、46.67%、13.33%,表明土壤全钾含量及速效钾含量受土壤其他理化指标的影响程度随土层厚度增加而降低。
3.2 植物多样性指数值变化
毛竹林乔木层植物多样性指数随毛竹重要值降低而增加,同已有研究结论相似[9,16];灌木层及草本层植物多样性指数随毛竹重要值降低而趋势性各异。同一层内植物不同多样性指数间相关性显著,乔木层的趋势性相异于灌木与草本层。乔木、灌木与草本层植物间Margalef指数、灌木与草本层植物间Margalef指数、乔木和草本层植物间Simpson指数正相关性显著。
3.3 土壤不同形态钾含量与毛竹林植物多样性指数的关系
因子分析表明,紫色土全钾含量主要受毛竹林乔木层及草本层植物多样性因子影响,速效钾含量主要受地形和草本层植物多样性因子的影响。土壤全钾含量随毛竹重要值降低而减少,呈低重要值类型<中重要值类型<高重要值类型的变化趋势;速效钾含量在毛竹低重要值类型中较高。土壤全钾含量与毛竹林乔木、灌木与草本层植物多样性指数间相关显著的直线和二次多项式指标对数量分别为50.00%和50.00%,土壤速效钾含量的相应值为13.89%和44.44%,表明土壤不同形态钾含量与植物多样性指数关系存在无关型、直线型和曲线型3种类型,印证了植物多样性的土壤生态功能多样性[9-11]。紫色土速效钾含量与乔木层及草本层植物多样性指数的相关性趋势相异于山地黄壤不同植被类型[13]、也不同于石灰岩山地淡竹林[8],反映了紫色土速效钾含量与植物多样性指数关系的特异性。大径级毛竹林土壤的全钾含量和速效钾含量较高[15]、植物多样性指数也较高[16]。
3.4 毛竹林是一个干扰性生态系统
年度间采笋或周期性伐竹均能引起毛竹林植物多样性指数变化,导致紫色土不同形态钾含量与毛竹林植物多样性指数关系的不确定性,周期内持续跟踪和周期间对比分析将有助于提高此研究结论的科学性。另一方面,紫色土钾[3,14]源及赋存形式有别于土壤其他养分元素[7-8,12-13,15],其流动性强,提高土壤孔隙度和林地光热水气条件的营林措施有利于增加不同形态钾含量,这也是土壤不同形态钾含量和毛竹林植物多样性指数值的直线型负相关显著的原因之一。由于二次多项式较直线型更能反映土壤不同形态钾含量与毛竹林植物多样性指数值的相关性,部分相关显著模式的拐点值可成为紫色土毛竹林经营的植物多样性指数管理依据[9-12]。
[1] 中国植物志编辑委员会. 中国植物志: 第9卷第1分册[M]. 北京: 科学出版社, 1996.
[2] 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会. 毛竹林丰产技术: GB/T20391-2006[S]. 北京: 中国标准出版社, 2006.
[3] 《赤水县综合农业区划》编写组. 赤水县综合农业区划[M]. 贵阳: 贵州人民出版社, 1989.
[4] 国家林业局. 国家林业局关于命名中国竹子之乡的决定(林造发[2006]190号) [EB/OL]. [2017-03-09]. http://www.forestry.gov.cn/portal/main/s/495/content-24205.html.
[5] 巫启新. 贵州毛竹林类型与林分结构的研究[J]. 竹子研究汇刊, 1983, 2(1): 112-124.
[6] 陈家洪, 张喜, 尹洁, 等. 抚育和采伐组合措施对毛竹林生产力的影响[J]. 安徽农业科学, 2010, 38(4): 2110-2111, 2152.
[7] 楼一平, 吴良如, 邵大方, 等. 毛竹纯林长期经营对林地土壤肥力的影响[J]. 林业科学研究, 1997, 10(2): 125-129.
[8] 施建敏, 王龙凤, 郭春兰, 等. 石灰岩山地淡竹林演替序列土壤养分及其与群落特征的关系[J]. 水土保持学报, 2016, 30(1): 153-158.
[9] 张佐玉, 张喜. 毛竹林紫色土壤主要水源涵养指标和地形要素及植物多样性指数的关系[J]. 竹子研究汇刊, 2015, 34(2): 1-7.
[10] Hector A, Bagchi R. Biodiversity and ecosystem multifunctionality[J]. Nature, 2007(448): 188-191.[11] Donald R Z, William E H, David C W, et al. Plant diversity, soil microbial communities, and ecosystem function: are there any links? [J]. Ecology, 2003, 84(8): 2042-2050.
[12] 李正才, 傅懋毅, 谢锦忠, 等. 毛竹竹阔混交林群落地力保持研究[J]. 竹子研究汇刊, 2003, 22(1): 32-37.
[13] 杨万勤, 钟章成, 陶建平, 等. 缙云山森林土壤速效K的分布特征及其与物种多样性的关系[J]. 生态学杂志, 2001, 20(6): 1-3.
[14] 占丽平, 李小坤, 鲁剑巍, 等. 土壤钾素运移的影响因素研究进展[J]. 土壤, 2012, 44(4): 548-553.
[15] 张佐玉, 张喜, 崔迎春, 等. 贵州毛竹林不同径级类型的土壤理化特性变化[J]. 竹子研究汇刊, 2012, 31(2): 6-10.
[16] 张喜, 张佐玉, 崔迎春, 等. 贵州赤水河流域毛竹林不同径级类型林分结构指标变化[J]. 世界竹藤通讯, 2014, 12(4): 1-5.
[17] 国家林业局. 森林土壤分析方法[M]. 北京: 中国标准出版社, 2000.
[18] 张金屯. 数量生态学[M]. 北京: 科学出版社, 2004.
Correlation of Purple Soil Potassium Contents with Terrain Factors and Plant Diversity Indices of Moso Bamboo (Phyllostachyspubescens)Forest at Downstream of Chishui River, Guizhou Province, China
Zhang Xi1Huo Da1Zhang Zuoyu1Wang Lili2Cui Yingchun1
(1 Guizhou Provincial Academy of Forestry, Guiyang 550011, China;2 Guizhou University of Finance and Economics, Guiyang 550005, China)
Correlation of the contents of potassium at different forms in purple soil with terrain factors and plant diversity indices of moso bamboo (Phyllostachyspubescens) forest were studied in this paper at the downstream of Chishui river, Guizhou Province, China. The results showed that: (1) No significant variation has been ever found between the purple soil’s contents of different forms of potassium in different types of moso bamboo forest and soil stratum, and as the soil thickness increases from surface to mother rock, the total potassium content is increasing and the available potassium content is reducing. The total potassium content is not significantly correlated with main physical indicators of soil while the correlation coefficient of the available potassium content and soil density and porosity of A (or B) stratum was significantly. The contents of different forms of potassium is significantly correlated with the cation exchange capacity of different soil stratum, and the content of organic matter, the total nitrogen content and the effective phosphorus content in surface soil; (2) Values of Margalef index, Simpson index, Shannon-Wiener index and Pielou index of tree layer in the moso bamboo forest increase with the reducing bamboo important values, and the plant diversity index values of shrub and herb layers are different. The plant diversity indices values of different layers are significantly related, and Margalef index values of different layers is correlated significantly with Simpson index of the tree and shrub layers; (3) The total potassium contents in different soil layers are significantly correlated with the plant diversity factors of tree and herbaceous layers, while the available potassium contents are significantly correlated with the factors of terrain and herb plant diversity. The linear and quadratic polynomial numbers of significant correlation between the soil total potassium content and plant diversity indices of tree, shrub and herb layers are 50.00% and 50.00%, respectively, and the corresponding values of the soil available potassium content are 13.89% and 44.44% in the moso bamboo forest. It showed that the relationship between different form of soil potassium contents and plant diversity indices includes 3 types, i.e., no correlation, straight correlation and curve correlation, which confirms the diversity of soil ecological function, which is embodied in moso bamboo forest plant diversity.
moso bamboo forest, purple soil, soil total potassium content, soil available potassium content, plant diversity index
贵州省科技厅项目(编号:黔科合R字[2014]2034);贵州省科技厅重大专项(编号:[2007]6002030201);黔科合院士站项目(编号:[2014]4006)。
张喜,博士,研究员,主要从事森林生态学和竹林经营研究。E-mail: 1003434130@qq.com。
10.13640/j.cnki.wbr.2017.03.001