Development of a compatible biomass model system for individual trees of major tree species in Hunan Province of central China
-
摘要:目的
建立杉木、马尾松和栎类3个树种的单木各组分相容性生物量模型系统，为实现湖南主要树种模型的科学化和标准化以及精准预估森林生物量提供技术支撑。
方法基于湖南省3个主要树种468株伐倒木数据（杉木160株、马尾松153株和栎类155株），分别比较单木6个组分生物量（地上、地下、树干、树皮、树枝和树叶生物量）基础模型并寻求拟合指标最优模型，进一步比较不同起源（天然林和人工林）的生物量模型，并选择引起参数差异显著的起源作为哑变量，构建比例控制的非线性联立方程组预估单木各组分生物量。
结果3个树种单木各组分生物量基础模型的拟合结果表明，地上B1、地下B2、树干B3和树皮B4这4个组分生物量模型的自变量选择胸径和树高，树枝B5和树叶B6生物量模型选择胸径、树高和冠幅；以起源为哑变量的模型显示，杉木的地上生物量B1（平均预估误差MPE较基础模型降低1.16%）、马尾松的地下生物量B2（MPE降低4.65%）、以及栎类的地上生物量B1（MPE降低3.24%）和地下生物量B2（MPE降低29.13%）均表现出天然林的生物量要显著高于人工林的生物量；以起源为哑变量的单木生物量模型系统的拟合效果较佳，3个树种地上生物量模型的决定系数R²均超过0.90，平均预估误差MPE低于4%；而地下生物量模型的R²同样均在0.90以上，其MPE则低于14%。
结论以起源为哑变量的单木各组分生物量模型拟合效果最好，其次是非线性联立方程组相容性模型系统，拟合效果最差的是基础模型。虽然模型系统的拟合效果略劣于哑变量模型，但模型系统能解决单木各组分生物量的相容性和可加性，且它们参数之间的差异不大，因此优先推荐采用相容性模型系统来预测湖南省杉木、马尾松和栎类3个树种的单木各组分生物量。
Abstract:ObjectiveThis paper aims to establish a compatible biomass model system for individual tree components of Cunninghamia lanceolata, Pinus massoniana, and Quercus species, so as to provide technical support for the scientific and standardized modeling of main tree species in Hunan Province of central China and the accurate prediction of forest biomass.
MethodBased on the felling data of 468 individual trees of three main tree species in Hunan Province (160 of Cunninghamia lanceolata, 153 of Pinus massoniana, and 155 of Quercus species), the basic models of biomass for six individual tree components (aboveground, underground, trunk, bark, branch, and leaf biomass) were compared respectively, and the optimal model with the best fitting indicators was sought. Further, the biomass models of different origins (natural forests and plantations) were compared. Dummy variable models were established by selecting parameters adding origin as variable with significant differences, and the biomass of individual tree components was predicted by constructing nonlinear simultaneous equations with proportion control.
ResultThe fitting results of basic biomass models for individual tree components of three tree species showed that for four components of aboveground biomass B1, underground biomass B2, trunk biomass B3, and bark biomass B4, the independent variables selected were DBH and tree height, while for the biomass models of branches B5 and leaves B6, DBH, tree height (H), and crown width (CW) were selected. The models with origin as a dummy variable indicated that for aboveground biomass B1 (the mean prediction error MPE was reduced by 1.16% compared with the base model) of Cunninghamia lanceolata, underground biomass B2 (MPE was reduced by 4.65%) of Pinus massoniana, and aboveground biomass B1 (MPE was reduced by 3.24%) and underground biomass B2 (MPE was reduced by 29.13%) of Quercus species, the biomass in natural forests was significantly higher than that in plantations. The fitting effect of individual tree biomass model system with origin as a dummy variable was better, with the determination coefficients (R2) of aboveground biomass models of three tree species all above 0.90, and the average predicted error MPE was below 4%, while the underground biomass models also showed R² above 0.90 and MPE below 14%.
ConclusionThe individual tree component biomass model with origin as a dummy variable has the best fitting effect, followed by the nonlinear simultaneous equation compatibility model system, and the basic model has the worst fitting effect. Although the fitting effects of model systems are slightly worse than the dummy variable model, these model systems can solve the compatibility and additivity of individual tree component biomass, and the differences among their parameters are not significant. Therefore, it is recommended to use compatibility model systems to predict the biomass of individual tree components of Cunninghamia lanceolata, Pinus massoniana, and Quercus species in Hunan Province of central China.
-
Keywords:
- biomass model /
- origin /
- dummy variable /
- nonlinear simultaneous system of equations /
- compatibility
-
化肥过量施用、土壤侵蚀、河岸植被破坏等易使面源污染物通过地表径流和地下渗流进入河湖水体中[1−2]。全球范围内近50%的地表水资源已经受到面源污染的影响。然而，面源污染具有很大的随机性，分布面积广，时间和空间变化幅度较大，所以在监测、控制和治理上存在诸多困难和复杂性，这使得减轻农业面源污染成为一项艰巨的任务，也使面源污染治理逐渐成为水环境保护与水资源利用中的重点和难点[3]。
河岸植被缓冲带是一个由土壤、水、植被等构成的生态系统，通过植物吸收、土壤吸附等一系列物理、化学和生物过程，截留、吸收和转化污染物，阻止径流水中的污染物进入水体，被认为是面源污染防治的有效管理措施[4]。黄斌斌等[5]研究发现，森林的水质净化率最高，达91%，其次是灌丛（87.80%）和草地（82.28%）。Haukos等[6]研究发现，40 ~ 60 m宽的植被缓冲带可以显著减少地表径流中的养分、重金属和其他污染物的负荷。农业区河流失去河岸植被，可能会对河流生态系统的结构和功能产生深远的影响[3，6]。因此，恢复和重建河岸缓冲带以减少农业活动对河流的污染负荷，应受到重视并成为流域治理的重要途径。
生物炭是由生物质在低氧或无氧、高温条件下裂解产生的一种富碳物质[7−8]。由于生物炭化学性质稳定，具有碱性、多孔的结构、良好的阳离子交换能力以及大的比表面积等特点，在土壤污染修复、提高土壤肥力、增加土壤碳封存等方面都有重要意义[9]。特别是生物炭通过改善土壤物理性质，改变土壤微生物群落结构和功能，影响土壤养分循环和有机质分解，并可能减少氮的淋溶，提高氮的滞留[10]。因此，生物炭作为一种环境功能性材料在改善土壤、净化水质、物质循环、土壤固碳等方面得到了广泛的应用。但生物炭的生态效益因植被类型、生物炭性质、土壤类型和土地利用类型而异[10]。以往的研究表明[11]，河岸缓冲带的植被类型对污染物截留去除效果差异显著，但也有研究[6]表明缓冲带的植被类型对污染物去除差异不大。生物炭通过影响土壤微生物群落和土壤酶活性对土壤微生物产生影响，但生物炭对微生物群落的影响仍存在争议。一些研究认为生物炭对微生物群落有促进作用或抑制作用，而另一些研究认为生物炭对微生物群落没有影响。因此，生物炭对污染物截留去除的效果，对土壤微生物群落变化的影响尚不完全清楚，有待进一步的研究。同时，现有的研究更多关注生物炭对农田的影响，且多室内盆栽试验，野外研究相对匮乏，生物炭添加对河岸林地影响的研究相对较少。
白洋淀的生态环境是雄安新区健康发展的重要保障，而入淀河流的水质影响着白洋淀的水环境质量。近年来，白洋淀流域通过减少外源污染排入、建立人工湿地、清淤和生态补水等方式进行了综合治理，但入淀河流仍然存在热点污染区域[12]。研究发现，白洋淀流域水体和沉积物的氮污染一直处于较高水平，是影响白洋淀流域水质的主要污染因子[13]。由于白洋淀流域河岸缓冲带存在生态功能脆弱、破碎化程度高的特点，且流域农业种植面积广、化肥施用量高，农业面源污染物很容易通过地表径流进入河流生态系统，影响流域污染治理效果[14]。因此，本研究选择白洋淀流域毛白杨林地为研究对象，探究不同生物炭添加量对河岸林土壤化学因子、微生物群落组成的影响，以及生物炭添加对河岸林地氮污染的阻控机制。研究旨在通过添加生物炭，强化河岸植被缓冲带的截留去污能力，减轻农业面源污染对流域水环境的威胁，为流域面源污染阻控及生态修复技术提供依据。
1. 研究区概况与研究方法
1.1 研究区域
白洋淀流域属海河流域大清河水系（113°39′ ~ 116°11′ E，38 °12′ ~ 39 °54′ N），跨山西和河北两省以及北京市，流域面积31 200 km2。属典型的温带大陆性季风气候，夏季高温多雨，冬季寒冷干燥，雨热同期，年平均气温7 ℃，多年平均降水量564 mm，年内降水量分配不均，70% ~ 80%集中在6—8月份，且多以暴雨形式出现[15]。流域地形地貌复杂，地势西高东低，其中山区占 43%，平原占57%，山区的土地利用类型主要是林地和草地等，土壤为褐土，平原地区主要是农田、城镇用地，土壤为潮土[16]。该流域农业种植面积广、化肥施用量高，氮肥和磷肥的平均施用量分别为100 ~ 300 kg/hm2和50 ~ 150 kg/hm2[17]。2017年雄安新区成立，围绕白洋淀而建，建设定位是中国未来生态与经济协调发展的先进新区，这对白洋淀流域生态环境提出了更高要求。由于其独特的自然条件和社会因素，白洋淀流域对维护华北地区生态环境具有不可替代的作用[18]。
1.2 样地概况和样地布设
白洋淀流域平原地区的河岸林基本都为毛白杨（Populus tomentosa），因此选择毛白杨河岸林为研究对象。于2023年5月18日，在白洋淀流域府河沿岸选择毛白杨人工林布设试验（图1a）。毛白杨行距为4 m，株距为3 m，密度为840株/hm2，林地郁闭度为0.75，林下植被稀少，偶见平车前（Plantago depressa）、艾（Artemisia argyi）、藜（Chenopodium album）等草本植物。林下稀疏生长毛白杨更新幼苗。
1.3 试验设计
所用的生物炭原材料为板栗壳。管式炉（型号KJ-T1200-S80-800S，规格为80 mm × 800 mm）通入高纯氮气30 min，以排净管中的空气。当管式炉温度升到指定温度（500 ℃）时，将装有板栗壳的石英管放入管式炉中进行炭化，炭化停留时间为2 h，自然冷却获得本试验所用生物炭。其表面积15.6 m2/g，容重0.23 g/cm3，pH8.56，全碳485.30 g/kg，全氮4.56 g/kg，全磷0.82 g/kg。
设置3种处理：未添加生物炭的对照（BC0）、添加生物炭3 t/hm2（BC3）、添加生物炭6 t/hm2（BC6）。每种处理3次重复，共设置9个样地，各处理组样地条件基本一致。样地设置在河流和农田之间的林地内（图1b），起始段和末端分别距农田、河流3 ~ 5 m，呈长方形布局，面积为2 m × 10 m。样地土壤为潮土，坡度为25°，坡长28 m。9块样地平行排列，间隔5 m。
在样地内开横沟施加生物炭，横沟长2 m，宽20 cm，深20 cm。每块样地开10条沟，沟与沟之间间距1 m。生物炭沟渠与地表径流方向垂直。在每块样地起始端（靠近农田处）添加氮肥（尿素，含N量为46%），添加带长3 m，宽0.5 m。模拟农业面源氮污染（图2），氮肥按照1 t/hm2，30 m2面积添加，每块样地的总添加量为3 kg。每月添加一次，连续添加4个月（5—8月），每次添加量为0.75 kg。
1.4 样品采集及指标测定
2023年9月中旬，采用五点取样法，清除地表枯枝落叶，采集样地中部0 ~ 20 cm表层土壤，装入塑封袋后放入保温箱中带回实验室。样品分为两部分，一份置于4 ℃冰箱用于土壤化学性质和土壤酶指标的测定；一份置于−80 ℃冰箱用于土壤微生物多样性的测定。
参照鲍士旦[19]提到的方法，进行土壤化学性质和酶的测定。土壤养分指标全氮的测定采用凯氏定氮法（Kjeltec 8400）；用重铬酸钾氧化加热法测定土壤有机质。土壤pH值使用玻璃电极pH计进行测定（水土比5∶1）；铵态氮和硝态氮采用流动分析仪测定（Smartchem 200）。土壤脲酶活性采用苯酚−次氯酸钠法测定，土壤硝酸还原酶活性采用酚二磺酸比色法测定，土壤亚硝酸还原酶活性用格里试剂显色与酶促反应前后NO−2-N变化表征。
1.5 细菌群落结构分析
采用高通量测序分析土壤细菌群落结构，在北京百迈克生物科技有限公司进行测序分析。采用E.Z.N.A. ® Soil DNA Kit提取土壤总DNA，以土壤总DNA为模板，采用341F（5′-ACTCCTACGGGAGGCAGCAG-3′）和806R（5′-GGACTACHVGGGTWTCTAAT-3′）对细菌16S rRNA基因的V3-V4高变区片段进行PCR扩增，扩增后的PCR产物用2%琼脂糖凝胶回收PCR产物，使用 AxyPrep DNA Gel Extraction kit进行纯化，Tris-HCl洗脱，2%琼脂糖电泳检测。利用QuantiFluorTM-ST检测定量，然后根据测序量需求将各样品按相应比例混合，构建Miseq文库，使用Illumina MiSeq平台测序。
1.6 数据处理和分析
采用SPSS26.0进行数据处理和分析，数据由平均值 ± 标准误表示，对土壤化学性质、土壤酶活性进行单因素方差分析和相关性分析。优势菌门与环境因子之间进行相关性热图分析，确定影响土壤细菌群落的关键环境因子。使用R、Origin 2021 软件进行绘图。
2. 结果与分析
2.1 不同处理组毛白杨林地土壤化学性质和酶活性
不同生物炭添加量的毛白杨林地土壤化学性质具有显著性差异（P < 0.05）（表1），其中BC6的土壤养分含量多高于BC3。与对照相比，BC3生物炭添加组全氮（TN）、铵态氮（NH+4-N）、硝态氮（NO−3-N）、有机质（SOM）、有机碳（SOC）分别提高10.39%、46.15%、8.57%、38.03%和 51.23%。BC6生物炭添加组TN、NH+4-N、NO−3-N、SOM、SOC分别提高 26.62%、92.31%、65.00%、78.05%和135.47%。
表 1 生物炭添加量对土壤养分和酶活性的影响Table 1. Effects of biochar addition on soil nutrients and enzyme activities指标 处理 BC0 BC3 BC6 全氮/(g·kg−1) 1.54 ± 0.05c 1.70 ± 0.1b 1.95 ± 0.02a 铵态氮/(mg·kg−1) 0.13 ± 0.04b 0.19 ± 0.03a 0.25 ± 0.07a 硝态氮/(mg·kg−1) 1.40 ± 0.21b 1.52 ± 0.23b 2.31 ± 0.60a 有机质/(g·kg−1) 8.02 ± 0.66c 11.07 ± 1.23b 14.28 ± 1.11a 有机碳/(g·kg−1) 4.06 ± 0.30c 6.14 ± 0.46b 9.56 ± 0.43a pH值 6.63 ± 0.21a 6.61 ± 0.29a 6.72 ± 0.19a 脲酶/(mg·g−1) 0.60 ± 0.03b 0.81 ± 0.06a 0.98 ± 0.07a 硝酸还原酶/(μg·g−1) 2.55 ± 0.64c 4.32 ± 1.01b 5.21 ± 1.35a 亚硝酸还原酶/(mg·g−1) 0.77 ± 0.15b 1.43 ± 0.19a 1.37 ± 0.25a 注：不同处理BC0、BC3、BC6的生物炭添加量分别为0、3、6 t/hm2。不同小写字母表示同一指标不同处理间差异显著（P<0.05）。下同。 不同生物炭添加量的毛白杨林地土壤酶活性具有显著性差异（P < 0.05），BC6生物炭添加组土壤酶活性高于BC3。同对照相比，BC3生物炭添加组脲酶（S-UE）、硝酸还原酶（S-NR）、亚硝酸还原酶（S-NIR）活性分别提高35.00%、69.41%和85.71%。BC6生物炭添加组S-UE、S-NR、S-NIR活性分别提高63.33%、104.31%和77.92%（表1）。
2.2 土壤酶活性与理化因子的相关性分析
表2显示：NO−3-N与NH+4-N、SOC、SOM之间呈现出极其显著正相关（P < 0.001）。TN与NO−3-N、SOC、SOM之间呈现极显著或显著正相关（P < 0.01或P < 0.05）。S-NR与S-NIR呈极其显著正相关（P < 0.001）。S-UE与pH呈极其显著正相关（P < 0.001）。S-NR、S-NIR分别与NH+4-N、NO−3-N、SOC、SOM呈极显著正相关（P < 0.01）。
表 2 土壤酶活性与养分因子相关分析Table 2. Correlation analysis of soil enzyme activity and nutrient factors指标 TN NO−3-N NH+4-N SOC SOM S-NR S-NIR S-UE pH TN 1 0.65** 0.21 0.55* 0.55* 0.41 0.41 0.42 0.42 NO−3-N 1 0.78*** 0.72*** 0.72*** 0.69** 0.69** 0.18 0.18 NH+4-N 1 0.37 0.37 0.60** 0.60** −0.10 −0.10 SOC 1 1.00*** 0.64** 0.64** −0.06 −0.06 SOM 1 0.64** 0.64** −0.06 −0.06 S-NR 1 1.00*** −0.21 −0.21 S-NIR 1 −0.21 −0.21 S-UE 1 1.00*** pH 1 注：*表示相关性在0.05水平显著，**表示相关性在0.01水平极显著，***表示相关性在0.001水平极其显著。TN.全氮；NO−3-N.硝态氮；NH+4-N.氨态氮；SOC.有机碳；SOM.有机质；S-NR.硝酸还原酶；S-NIR.亚硝酸还原酶；S-UE.脲酶。下同。 2.3 生物炭对毛白杨林地土壤细菌群落结构的影响
2.3.1 土壤细菌OTU数量分析
对不同处理组细菌操作性分类单元（operational taxonomic unit，OTU）进行韦恩分析。3个处理9个样品得到812 706个细菌有效序列，对土壤样品序列通过拆分、去冗余后，在97%相似度下进行OTU聚类共得到36 720个细菌OTU，不同处理土壤样品中共有的细菌OTU总数有1 023个。其中BC0处理14 923个OTU，BC3处理12 469个 OTU，BC6处理13 344个 OTU。BC0、BC3、BC6所特有的 OTU 数目分别为12 616、10 130、10 981个，占比分别为 34.36%、27.59%、29.91%。与BC0相比较，2种生物炭添加量都降低了土壤细菌中的OTU数量。
2.3.2 微生物群落α多样性分析
表3显示：土壤样本的覆盖率均在99%以上，且无显著差异（P > 0.05），这表明各组测得的数据都具有代表性，测序结果能够准确反应供试土样真实情况。各处理间Ace、Chao指数、Shannon多样性指数和Simpson多样性指数均无显著性差异（P > 0.05），生物炭添加对土壤细菌的α多样性影响较小。
表 3 不同处理生物覆盖度和多样性指数Table 3. Biological coverage and diversity indices under different treatments处理 Ace Chao指数 Shannon多样性指数 Simpson多样性指数 覆盖度/% BC0 2 802.2 ± 45.62 2 796.4 ± 45.58 10.12 ± 0.04 0.99 ± 0.00 99.96 ± 0.00 BC3 2 093.6 ± 125.69 2 089.1 ± 123.69 9.76 ± 0.39 0.99 ± 0.00 99.97 ± 0.00 BC6 2 587.6 ± 228.70 2 581.2 ± 226.47 10.14 ± 0.02 0.99 ± 0.00 99.96 ± 0.00 注：数值为均值 ± 标准误差。 2.3.3 土壤微生物门水平群落结构组成
在相似性水平97%的条件下，对OTU的代表序列做分类分析，共得到42个门，105纲，314个目，624个科，
1105 个属，1343 个种的土壤细菌（图3）。不同处理下土壤中主要的细菌群落组成相似，但相对丰度却有差异。在门分类水平上，相对丰度前10的细菌门包括变形杆菌门（Proteobacteria）、酸杆菌门（Acidobacteriota）、放线菌门（Actinobacteriota）、未分类菌门、绿弯菌门（Chloroflexi）、拟杆菌门（Bacteroidota）、芽单胞菌门（Gemmatimonadota）、黏球菌门（Myxococcota）、甲基微菌门（Methylomirabilota）、疣微菌门（Verrucomicrobiota）。![]()
图 3 土壤细菌门水平群落结构组成Figure 3. Composition of community structure of soil bacterial phylum level前六大优势菌门的相对丰度占总群落的比重接近80%，各处理变形杆菌门的相对丰度分别为BC6（27.86%） > BC3（27.52%） > BC0（26.76%）；酸杆菌门：BC3（23.84%） > BC0（22.28%） > BC6（21.37%）；放线菌门为BC0（13.09%） > BC3（9.61%） > BC6（8.12%）；未分类菌门为BC6（9.15%） > BC0（6.64%） > BC3（4.86%）；绿弯菌门为BC0（7.06%） > BC3（6.29%） > BC6（5.22%）；拟杆菌门为BC6（7.19%） > BC3（4.84%） > BC0（4.36%）。
总体来看，生物炭添加提高了变形杆菌门、未分类菌门和拟杆菌门的相对丰度；降低了酸杆菌门的相对丰度。变形杆菌门在生物炭添加组的相对丰度占比较高，其中含有大量的反硝化细菌，是影响反硝化效果的重要因素。反硝化作用是土壤中氮素转化与去除的主要方式之一。因此，生物炭添加通过影响土壤主要细菌门的相对丰度（如变形杆菌门），进而影响土壤中氮素的转化与去除。
2.4 环境因子与微生物群落结构相关性
土壤化学性质、酶活性与相对丰度前10的细菌门的相关性分析如图4所示。NO−3-N、NH+4-N、SOC、SOM、S-NR、S-NIR分别与变形杆菌门、拟杆菌门、疣微菌门正相关，其中SOC、SOM、S-NR、S-NIR分别与变形杆菌门、拟杆菌门、疣微菌门显著正相关（P < 0.05），NO−3-N、NH+4-N与变形杆菌门显著正相关（P < 0.05）。
![]()
图 4 土壤细菌群落结构和土壤化学性质相关性分析Figure 4. Correlation analysis of soil bacterial community structure and soil chemical propertiesNO−3-N、NH+4-N、SOC、SOM、S-NR、S-NIR分别与芽单胞菌门、酸杆菌门、绿弯菌门负相关，其中SOC、SOM、S-NR、S-NIR分别与芽单胞菌门、酸杆菌门、绿弯菌门显著负相关（P < 0.05），NO−3-N、NH+4-N与芽单胞菌门显著负相关（P < 0.05），SOC、SOM还与放线菌门显著负相关（P < 0.05）。
可见，SOC、SOM、S-NR、S-NIR是影响土壤微生物群落门水平的主要因子。
3. 讨 论
3.1 生物炭添加对毛白杨林地土壤养分和酶活性的影响
河流和河岸林地联系密切，通过物质循环、能量流动、动植物栖息、洪水和光照等相互作用高度交织在一起[3]。根据流域治理中的“源头减排−过程阻断−末端治理”模式，通过对污染物迁移、转化、扩散过程的拦截和阻断，能够有效减少污染物进入河湖水体[20]。研究[21−22]表明，生物炭对土壤水分和营养元素有较强吸持能力，能提高土壤保水能力，改良土壤质量，减少土壤养分淋失。本研究发现，6 t/hm2（BC6）处理组对氮素的截持能力更强，这可能是因为生物炭能够通过改善土壤结构来增强土壤的持水能力和吸附点位[23]。生物炭富含碳元素，能够提升土壤中的无机碳含量，从而增强土壤的养分保持能力和水分保持能力。此外，它还可以促进林地表面及土壤内部木质素的分解[24]。张晗芝[25]研究发现，生物炭施于土壤中能显著提高土壤碳含量。与对照组相比，6 t/hm2（BC6）处理组林地SOC、SOM分别提高了78.05%、135.47%。生物炭不仅增加了林地土壤有机质含量，还有助于降低土壤压实和侵蚀、增强土壤持水能力，并维持立地生产力[26]。同时，乔木叶片木质素含量高，分解有限[27]，导致微生物可利用碳源不足，生物炭添加使林地土壤碳含量提高，为微生物提供充足的碳源。
微生物群落及其产生的胞外酶在有机质分解和养分循环等过程中扮演着重要的调节角色，能够反映生物炭等土壤添加剂对土壤养分循环的潜在长期影响[28]。与赵蕊蕊等[24]研究结果一致，本研究中生物炭处理组林地土壤酶活性显著高于对照（P < 0.05）。土壤硝酸还原酶能够将硝酸还原成亚硝酸，而亚硝酸还原酶能将亚硝酸还原成氮气，进而影响土壤中氮素的转化与去除。因此，随生物炭添加量增加，毛白杨林地土壤氮循环相关酶活性也随之提高，有助于土壤中氮素被还原成氮气并从土壤中去除。
3.2 生物炭添加对毛白杨林地土壤细菌群落结构的影响
生物炭添加在改变土壤化学性质的同时，也可影响土壤微生物群落结构和多样性[29]。有研究[30]表明，随着稻壳生物炭添加浓度的增加，樟子松（Pinus sylvestris var. mongolica）的细菌群落总数和丰度出现了下降。本研究发现，生物炭添加不同处理的Ace、Chao指数、Shannon多样性指数和Simpson多样性指数均无显著性差异（P > 0.05），生物炭添加对毛白杨林地土壤细菌群落α多样性的影响小。这可能与生物炭添加时间较短有关，有待进一步长期观测研究。本研究中随着生物炭添加浓度的增高，毛白杨林地土壤细菌OTU总数有所下降，且微生物群落结构发生了明显变化。优势门水平上，生物炭添加提高了变形菌门、拟杆菌门的相对丰度；降低了放线菌门的相对丰度。
河岸植被缓冲带的土壤干湿交替频繁，有利于土壤中微生物的生长，也有利于好氧硝化作用和厌氧反硝化作用的进行。在微生物作用下，铵态氮被氧化成亚硝态氮和硝态氮，然后再被反硝化细菌还原成N2或N2O等气体排入空气中，反硝化作用被认为是清除土壤里氮素的最佳途径[31]。变形菌门在反硝化过程中起着主导作用，其相对丰度的多少是影响反硝化效果的重要因素。本研究中，变形菌门在生物炭添加组的相对丰度占比较高。Zhang等[8]研究也发现生物炭的的添加能显著提高系统中的变形杆菌门的比例，强化系统脱氮性能，并通过优化生物炭的添加量，使系统中全氮的最终去除率达67%。与此同时，BC6和BC3的拟杆菌门的相对丰度较高，这是由于该菌门是一种典型的兼性厌氧菌，能适应交替的好氧缺氧环境，将亚硝酸盐还原为NO，有利于高效脱氮[32]。
生物炭比表面积大于土壤，为微生物提供附着表面和栖息环境，从而影响土壤微生物活性和群落结构，进而影响胞外酶的活性[33−34]。本研究中生物炭添加对毛白杨林地土壤养分提高和土壤酶活性增加也具有显著的促进作用，这可能与相关土壤微生物的活性增大和生物量增加有关。相关性分析表明，SOC、SOM、S-NR、S-NIR是影响土壤微生物群落门水平的主要因子，SOC、SOM与S-NR、S-NIR极显著正相关（P < 0.01），这与土壤微生物多为异养有关，生物炭添加会为微生物提供碳源和能量，从而提高了土壤氮循环相关酶活性[34]。
4. 结 论
（1）生物炭添加能够显著提升毛白杨林地土壤对氮素的截留能力以及氮循环相关酶活性。生物炭处理组土壤养分均高于对照，其中6 t/hm2处理效果更好。
（2）高通量测序结果表明，随着生物炭添加量的增大，毛白杨林地土壤细菌OTU总数有所下降，但对土壤细菌多样性影响较小。生物炭添加影响了毛白杨林地土壤优势细菌门的相对丰度，促进了河岸林地氮素的转化与去除。
（3）相关性分析表明，土壤化学性质（SOC、SOM）以及酶活性（S-NR、S-NIR）是影响土壤微生物群落门水平结构的主要因素。
-
![]()
图 1 湖南省3个树种样木伐倒木分布
Figure 1. Distribution of felled trees for three tree species in Hunan Province
![]()
图 2 单木生物量基础模型变量间的相关性
B1 ~ B6为地上、地下、树干、树皮、树枝和树叶的生物量。D. 胸径；H. 树高；T. 年龄；CW. 冠幅；CL. 冠长；CS. 冠形率；CR. 树冠率。下同。
Figure 2. Correlation between variables of individual tree biomass base model
表 1 湖南省3个树种生物量建模数据统计
Table 1 Statistics of biomass modeling variables for three tree species in Hunan Province
树种 变量 样本数 平均值 标准差 最小值 最大值 杉木
Cunninghamia
lanceolata胸径/cm 160 19.88 7.77 2.20 35.30 树高/m 160 14.27 5.37 2.50 25.60 冠幅/m 160 4.65 1.47 0.80 9.45 地上生物量/kg 160 122.71 103.81 0.60 456.00 地下生物量/kg 74 20.55 16.65 0.20 77.14 树干/kg 160 85.95 74.43 0.27 320.04 树皮/kg 160 14.03 11.41 0.03 55.51 树枝/kg 160 14.63 13.49 0.10 59.19 树叶/kg 160 8.09 8.18 0.10 39.25 马尾松
Pinus
massoniana胸径/cm 153 17.81 8.72 3.50 38.20 树高/m 153 14.19 5.29 3.90 27.60 冠幅/m 153 4.07 1.98 0.90 11.70 地上生物量/kg 153 148.95 145.15 1.45 693.14 地下生物量/kg 54 32.16 40.69 0.43 201.43 树干/kg 153 110.04 113.89 0.76 489.65 树皮/kg 153 10.42 9.37 0.23 48.22 树枝/kg 153 23.81 24.09 0.16 143.22 树叶/kg 153 4.69 5.37 0.06 36.24 栎类
Quercus spp.胸径/cm 155 16.72 8.89 3.50 38.90 树高/m 155 12.97 4.69 3.70 24.70 冠幅/m 155 4.79 2.06 1.15 10.50 地上生物量/kg 155 189.62 201.10 2.22 940.20 地下生物量/kg 49 48.15 49.42 0.84 189.40 树干/kg 155 121.30 133.81 1.25 667.73 树皮/kg 155 21.08 25.44 0.25 121.89 树枝/kg 155 39.38 41.35 0.38 184.03 树叶/kg 155 7.85 7.74 0.06 41.01 
下载: 导出CSV
表 2 单木生物量基础模型的参数估计与评价指标
Table 2 Parameter estimation and evaluation indicators for the basic model of individual tree biomass
树种 组分 变量 参数 评价指标 β1 β2 β3 β4 RMSE R2 TRE MPE/% AIC 杉木
Cunninghamia
lanceolataB1 D 0.049 6 2.531 3 27.618 0 0.928 8 −0.605 7 3.536 6 1 521.970 5 D + H 0.037 2 2.169 8 0.507 9 25.437 6 0.939 6 −0.381 9 3.267 9 1 497.653 5 D + T 0.051 1 2.494 8 0.025 9 27.600 2 0.928 9 −0.658 2 3.545 7 1 523.763 4 D + H + CW 0.041 0 2.186 7 0.498 3 −0.079 0 25.230 9 0.940 6 −0.429 2 3.251 9 1 497.042 3 B2 D 0.013 7 2.422 0 4.816 9 0.899 8 0.458 3 7.963 4 963.143 3 D + H 0.013 6 2.627 4 −0.229 7 4.787 7 0.901 0 0.406 1 7.940 7 963.198 8 D + T 0.014 2 1.811 0 0.591 7 3.975 5 0.931 7 −0.370 7 6.593 6 903.711 8 D + H + CW 0.013 9 2.410 7 −0.089 5 0.183 9 4.705 4 0.904 4 0.280 1 7.829 6 959.650 8 B3 D 0.029 4 2.583 8 20.818 5 0.921 3 −0.732 6 3.806 0 1 431.529 9 D + H 0.019 4 2.051 4 0.743 5 17.708 3 0.943 0 −0.502 4 3.247 8 1 381.750 5 D + T 0.033 7 2.413 8 0.119 9 20.567 7 0.923 2 −0.987 6 3.772 3 1 429.651 3 D + H + CW 0.019 7 2.053 6 0.742 4 −0.010 9 17.705 3 0.943 1 −0.508 2 3.257 8 1 383.697 5 B4 D 0.010 1 2.349 6 3.995 9 0.876 5 −0.200 6 4.474 0 903.346 3 D + H 0.008 5 2.112 1 0.330 8 3.914 5 0.881 5 −0.028 7 4.397 0 898.762 8 D + T 0.011 9 2.138 8 0.152 1 3.943 0 0.879 7 −0.536 8 4.429 0 901.081 2 D + H + CW 0.009 1 2.126 7 0.322 9 −0.064 8 3.902 8 0.882 2 −0.065 3 4.398 1 899.803 6 B5 D 0.010 9 2.341 4 7.457 2 0.692 7 −0.622 7 8.009 7 1 103.000 0 D + H 0.012 1 2.489 4 −0.204 1 7.437 7 0.694 3 −0.743 6 8.014 5 1 104.162 1 D + T 0.007 0 2.947 5 −0.445 2 7.161 6 0.716 5 0.038 1 7.716 9 1 092.054 0 D + H + CW 0.019 8 2.567 7 −0.249 7 −0.388 1 7.243 5 0.710 0 −1.052 0 7.830 5 1 097.692 2 B6 D 0.002 4 2.631 0 4.655 7 0.673 8 −0.058 5 9.038 7 952.248 4 D + H 0.002 9 2.810 9 −0.272 9 4.639 4 0.676 0 −0.395 5 9.036 2 953.127 6 D + T 0.001 6 3.097 1 −0.329 3 4.574 1 0.685 1 0.448 5 8.909 0 948.593 2 D + H + CW 0.004 8 2.993 0 −0.394 0 −0.472 2 4.466 1 0.699 8 −0.390 7 8.727 0 942.948 8 马尾松
Pinus
massonianaB1 D 0.135 1 2.325 4 42.301 2 0.914 5 −0.209 1 4.566 2 1 586.108 4 D + H 0.051 3 2.001 6 0.699 8 31.849 7 0.951 5 0.079 6 3.449 6 1 501.269 5 D + T 0.125 3 2.440 0 −0.084 6 42.115 4 0.915 3 −0.026 2 4.561 5 1 586.761 4 D + H + CW 0.055 8 1.950 4 0.700 7 0.046 9 31.783 2 0.951 7 −0.001 7 3.454 1 1 502.630 3 B2 D 0.001 4 3.225 4 10.248 5 0.870 3 8.799 7 14.516 3 1 152.297 7 D + H 0.000 8 2.714 2 0.794 6 9.774 9 0.882 0 8.821 6 13.892 2 1 139.818 6 D + T 0.001 4 3.306 3 −0.063 6 10.241 9 0.870 5 8.806 1 14.556 0 1 154.100 4 D + H + CW 0.002 3 2.197 4 0.820 9 0.295 8 9.576 4 0.886 8 7.047 1 13.656 4 1 135.540 7 B3 D 0.084 2 2.379 1 44.683 2 0.845 1 −0.396 0 6.529 0 1 602.871 9 D + H 0.018 1 1.834 6 1.142 4 29.328 5 0.933 2 −0.224 9 4.299 9 1 476.034 9 D + T 0.080 6 2.446 7 −0.049 9 44.648 7 0.845 3 −0.295 5 6.546 0 1 604.635 5 D + H + CW 0.014 2 1.976 2 1.138 3 −0.122 8 29.051 0 0.934 5 0.025 2 4.273 7 1 475.125 1 B4 D 0.024 3 2.029 0 3.281 1 0.876 6 −0.131 4 5.061 0 803.778 2 D + H 0.021 1 1.976 1 0.107 6 3.268 4 0.877 5 −0.107 3 5.058 5 804.591 1 D + T 0.021 1 2.326 5 −0.235 1 3.201 4 0.882 5 0.254 5 4.954 7 798.249 0 D + H + CW 0.018 5 2.069 4 0.100 3 −0.087 3 3.257 0 0.878 4 0.029 1 5.058 1 805.526 2 B5 D 0.027 5 2.249 3 10.666 4 0.802 7 −0.073 4 7.204 4 1 164.526 0 D + H 0.079 7 2.554 3 −0.717 7 9.471 8 0.844 4 −0.939 9 6.419 2 1 130.181 2 D + T 0.024 2 2.436 9 −0.137 8 10.618 0 0.804 4 0.263 7 7.196 0 1 165.136 4 D + H + CW 0.168 5 1.890 3 −0.672 1 0.739 1 7.950 5 0.890 4 −0.655 7 5.406 5 1 078.606 2 B6 D 0.011 2 2.015 9 3.814 6 0.492 1 1.515 1 13.091 4 849.878 8 D + H 0.053 4 2.247 5 −0.806 8 3.666 5 0.530 8 −1.405 0 12.625 7 839.761 5 D + T 0.006 0 2.621 1 −0.387 9 3.785 4 0.499 9 3.898 5 13.035 3 849.531 4 D + H + CW 0.127 9 1.054 9 −0.698 6 1.526 8 3.159 6 0.651 6 2.564 0 10.917 3 796.234 6 栎类
Quercus spp.B1 D 0.285 0 2.200 1 63.745 6 0.898 3 −1.034 2 5.364 7 1 745.037 7 D + H 0.030 5 1.978 4 1.046 6 34.912 3 0.969 5 1.520 7 2.947 9 1 559.194 4 D + T 0.139 4 2.906 3 −0.412 4 60.649 8 0.907 9 0.616 0 5.121 1 1 731.505 2 D + H + CW 0.032 7 1.925 2 1.033 2 0.075 8 34.683 0 0.969 9 1.436 1 2.938 3 1 559.138 5 B2 D 0.017 4 2.569 3 11.623 0 0.891 0 1.516 9 12.579 5 1 214.041 6 D + H 0.015 7 2.557 1 0.051 3 11.618 1 0.891 1 1.558 4 12.615 8 1 215.910 1 D + T 0.025 7 1.703 8 0.624 9 11.178 3 0.899 2 1.142 0 12.138 3 1 203.869 8 D + H + CW 0.006 6 3.124 2 0.066 8 −0.555 1 11.340 5 0.896 2 2.921 9 12.355 5 1 210.364 3 B3 D 0.170 4 2.222 4 50.787 1 0.854 2 −1.386 8 6.682 4 1 674.133 2 D + H 0.008 6 1.955 4 1.362 4 23.420 2 0.969 0 2.090 3 3.091 8 1 434.631 8 D + T 0.057 5 3.245 6 −0.585 4 47.522 2 0.872 4 1.103 6 6.273 6 1 655.402 5 D + H + CW 0.008 7 1.946 3 1.360 2 0.012 7 23.416 1 0.969 0 2.076 5 3.101 6 1 436.577 3 B4 D 0.023 5 2.298 9 11.701 8 0.786 3 −2.784 8 8.879 5 1 216.148 4 D + H 0.000 8 1.881 2 1.690 6 6.652 7 0.930 9 0.342 0 5.064 9 1 041.954 5 D + T 0.007 0 3.446 1 −0.656 6 11.136 5 0.806 4 −0.282 6 8.478 5 1 202.698 7 D + H + CW 0.000 8 1.832 1 1.680 3 0.067 5 6.638 5 0.931 2 0.283 8 5.071 0 1 043.291 5 B5 D 0.063 7 2.172 9 15.842 1 0.851 3 0.393 5 6.411 6 1 310.662 4 D + H 0.058 8 2.164 3 0.038 3 15.837 1 0.851 4 0.485 8 6.430 8 1 312.564 2 D + T 0.069 9 2.070 0 0.062 8 15.830 3 0.851 6 0.179 9 6.428 0 1 312.429 9 D + H + CW 0.073 0 1.970 5 −0.012 4 0.298 1 15.537 4 0.857 0 0.232 2 6.330 2 1 308.603 4 B6 D 0.072 5 1.618 7 4.604 2 0.642 0 −0.344 6 9.327 5 925.124 4 D + H 0.244 1 1.886 5 −0.748 8 4.366 5 0.678 0 −1.467 4 8.875 1 910.581 5 D + T 0.112 2 0.785 7 0.587 7 4.474 1 0.661 9 −1.495 7 9.094 0 918.183 5 D + H + CW 0.270 9 1.723 5 −0.780 9 0.275 2 4.336 9 0.682 3 −1.435 1 8.844 4 910.460 1
下载: 导出CSV
表 3 以起源为哑变量的单木生物量模型的参数估计与评价指标
Table 3 Parameter estimation and evaluation indicators for individual tree biomass model with origin as dummy variable
树种 组分 参数 评价指标 β1 β2 α1 β3 β4 RMSE R2 TRE MPE/% F P 杉木
Cunninghamia
lanceolataB1 0.037 9 2.200 9 −0.052 0 0.523 4 25.061 2 0.941 4 −0.380 2 3.230 0 4.721 3 0.031 3 B2 0.013 6 2.458 8 0.168 3 −0.230 0 4.787 6 0.901 0 0.434 5 7.966 3 0.006 7 0.935 1 B3 0.019 7 2.076 2 −0.042 3 0.757 6 17.539 3 0.944 1 −0.497 3 3.227 3 3.019 5 0.084 2 B4 0.008 6 2.143 7 −0.051 0 0.343 8 3.886 0 0.883 2 −0.037 1 4.379 2 2.299 8 0.131 4 B5 0.022 5 2.688 7 −0.165 6 −0.223 2 −0.433 2 7.005 7 0.728 7 −1.046 0 7.598 3 10.700 2 0.001 3 B6 0.004 7 2.986 4 0.008 0 −0.394 6 −0.470 7 4.466 0 0.699 8 −0.391 7 8.755 2 0.012 0 0.913 1 马尾松
Pinus
massonianaB1 0.052 9 1.991 6 −0.003 8 0.703 7 31.836 1 0.951 6 0.016 3 3.459 9 0.126 7 0.722 4 B2 0.003 2 2.123 5 −0.084 2 1.050 9 9.289 0 0.895 5 4.425 6 13.246 6 15.995 0 0.000 1 B3 0.016 3 1.870 9 0.013 1 1.125 7 29.233 8 0.933 7 −0.014 5 4.300 6 0.967 0 0.327 0 B4 0.020 2 1.987 3 0.006 2 0.105 0 3.266 8 0.877 6 −0.010 1 5.073 3 0.143 6 0.705 3 B5 0.205 4 1.869 7 −0.024 1 −0.677 8 0.711 2 7.884 3 0.892 2 −1.072 1 5.379 9 2.497 0 0.116 2 B6 0.296 1 1.038 1 −0.122 5 −0.733 9 1.324 2 2.965 9 0.693 0 0.079 0 10.283 2 19.963 3 0.000 0 栎类
Quercus spp.B1 0.034 0 1.987 2 −0.018 1 1.003 8 33.637 9 0.971 8 1.338 6 2.852 3 6.400 2 0.012 4 B2 0.170 2 2.373 8 −0.132 0 −0.516 3 8.463 9 0.943 0 0.382 9 8.941 0 100.782 2 0.000 0 B3 0.008 4 1.960 0 −0.001 5 1.362 2 23.048 7 0.970 1 1.898 9 3.055 2 0.039 7 0.842 3 B4 0.001 1 1.919 5 −0.058 8 1.524 3 6.207 9 0.940 1 0.436 6 4.734 3 23.367 6 0.000 0 B5 0.094 2 1.988 8 −0.052 6 −0.092 2 0.278 8 14.942 8 0.868 6 0.258 4 6.121 4 12.170 0 0.000 6 B6 0.211 4 1.735 8 0.030 6 −0.705 7 0.254 8 4.323 5 0.686 3 −1.333 8 8.882 3 1.769 9 0.185 4
下载: 导出CSV
表 4 以起源为哑变量的单木生物量模型系统的参数估计与评价指标
Table 4 Parameter estimation and evaluation indicators for individual tree biomass aggregation model systems with origin as dummy variable
树种 组分 参数 评价指标 β1 β2 α1 β3 β4 RMSE R2 TRE MPE/% 杉木
Cunninghamia
lanceolataB1 0.038 5 2.203 8 −0.052 8 0.516 1 25.061 9 0.941 4 −0.414 8 3.230 1 B2 0.013 6 2.627 4 / −0.229 7 4.787 7 0.901 0 0.406 3 7.940 7 B3 / / / / 17.522 1 0.944 2 −0.360 2 3.358 0 B4 0.437 2 0.058 1 / −0.411 1 3.892 4 0.882 8 0.076 0 4.568 4 B5 0.972 8 0.627 2 −0.104 8 −1.095 6 −0.188 0 7.110 9 0.720 5 −1.078 3 8.006 4 B6 0.224 3 0.958 3 / −1.303 7 −0.144 1 4.595 2 0.682 2 −0.633 5 9.351 9 马尾松
Pinus
massonianaB1 0.051 4 2.000 5 / 0.700 2 31.849 7 0.951 5 0.072 1 3.449 6 B2 0.003 2 2.123 5 −0.084 2 1.050 9 9.289 0 0.895 5 4.425 5 13.246 6 B3 / / / / 28.797 9 0.935 6 0.405 0 4.405 0 B4 1.836 9 0.104 1 / −1.146 6 3.261 6 0.878 0 −0.553 6 5.266 6 B5 17.109 7 −0.005 9 / −2.006 5 0.827 5 7.857 6 0.892 9 −0.586 1 5.555 9 B6 25.778 4 −0.590 7 −0.157 0 −1.932 1 0.850 0 3.101 0 0.664 4 −2.864 2 11.140 9 栎类
Quercus spp.B1 0.034 0 1.987 2 −0.018 1 1.003 8 33.637 9 0.971 8 1.338 5 2.852 3 B2 0.170 2 2.373 8 −0.132 0 −0.516 3 8.463 9 0.943 0 0.382 9 8.941 0 B3 / / / / 23.375 9 0.969 3 1.653 7 3.243 7 B4 0.227 8 −0.104 8 −0.054 5 0.054 4 6.360 5 0.937 1 −1.562 0 5.077 9 B5 8.657 5 0.100 0 −0.052 1 −1.394 4 0.199 3 14.997 4 0.867 6 0.882 5 6.409 9 B6 3.188 2 0.454 6 / −1.628 9 −0.474 7 4.723 6 0.625 6 7.111 5 10.124 6
下载: 导出CSV
-
[1] 董利虎, 李凤日. 大兴安岭东部主要林分类型乔木层生物量估算模型[J]. 应用生态学报, 2018, 29(9): 2825−2834. Dong L H, Li F R. Stand-level biomass estimation models for the tree layer of main forest types in East Daxing'an Mountains, China[J]. Chinese Journal of Applied Ecology, 2018, 29(9): 2825−2834.
[2] Asrat Z, Eid T, Gobakken T, et al. Aboveground tree biomass prediction options for the dry Afromontane forests in south-central Ethiopia[J]. Forest Ecology and Management, 2020, 473: 118335. doi: 10.1016/j.foreco.2020.118335
[3] Fehrmann L, Kleinn C. General consideration about the use of allometric equations for biomass estimation on the example of Norway spruce in central Europe[J]. Forest Ecology and Management, 2006, 236(2): 412−421.
[4] 罗云建, 张小全, 王效科, 等. 森林生物量的估算方法及其研究进展[J]. 林业科学, 2009, 45(8): 129−134. doi: 10.3321/j.issn:1001-7488.2009.08.023 Luo Y J, Zhang X Q, Wang X K, et al. Forest biomass estimation methods and their prospects[J]. Scientia Silvae Sinicae, 2009, 45(8): 129−134. doi: 10.3321/j.issn:1001-7488.2009.08.023
[5] 赵菡, 雷渊才, 符利勇. 江西省不同立地等级的马尾松林生物量估计和不确定性度量[J]. 林业科学, 2017, 53(8): 81−93. doi: 10.11707/j.1001-7488.20170810 Zhao H, Lei Y C, Fu L Y. Biomass and uncertainty estimates of Pinus massoniana forest for different site classes in Jiangxi Province[J]. Scientia Silvae Sinicae, 2017, 53(8): 81−93. doi: 10.11707/j.1001-7488.20170810
[6] 曾伟生, 唐守正. 东北落叶松和南方马尾松地下生物量模型研建[J]. 北京林业大学学报, 2011, 33(2): 1−6. Zeng W S, Tang S Z. Establishment of below-ground biomass equations for larch in northeastern and Masson pine in southern China[J]. Journal of Beijing Forestry University, 2011, 33(2): 1−6.
[7] 曹磊, 刘晓彤, 李海奎, 等. 广东省常绿阔叶林生物量生长模型[J]. 林业科学研究, 2020, 33(5): 61−67. Cao L, Liu X T, Li H K, et al. Biomass growth models for evergreen broad-leaved forests in Guangdong[J]. Forest Research, 2020, 33(5): 61−67.
[8] 王佳慧, 李凤日, 董利虎. 基于不同预测变量的天然椴树可加性地上生物量模型构建[J]. 应用生态学报, 2018, 29(11): 3685−3695. Wang J H, Li F R, Dong L H. Additive aboveground biomass equations based on different predictors for natural Tilia Linn[J]. Chinese Journal of Applied Ecology, 2018, 29(11): 3685−3695.
[9] 王金池, 邓华锋, 黄国胜, 等. 天然云杉相容性生物量估算模型[J]. 应用生态学报, 2017, 28(10): 3189−3196. Wang J C, Deng H F, Huang G S, et al. Compatible biomass models of natural spruce (Picea asperata)[J]. Chinese Journal of Applied Ecology, 2017, 28(10): 3189−3196.
[10] Fu L Y, Lei Y C, Zeng W S, et al. Uncertainty assessment in aboveground biomass estimation at the regional scale using a new method considering both sampling error and model error[J]. Canadian Journal of Forest Research, 2017, 47(8): 1095−1103. doi: 10.1139/cjfr-2016-0436
[11] 李海奎, 宁金魁. 基于树木起源、立地分级和龄组的单木生物量模型[J]. 生态学报, 2012, 2(3): 740−757. Li H K, Ning J K. Individual tree biomass model by tree origin, site classes, and age groups[J]. Acta Ecologica Sinica, 2012, 2(3): 740−757.
[12] 曾伟生, 王雪军, 陈振雄, 等. 林分起源对立木生物量模型的影响分析[J]. 林业资源管理, 2014(2): 40−45. Zeng W S, Wang X J, Chen Z X, et al. Analysis of impacts of forest origin on single tree biomass models[J]. Forest Resources Management, 2014(2): 40−45.
[13] 薛春泉, 徐期瑚, 林丽平, 等. 广东主要乡土阔叶树种含年龄和胸径的单木生物量模型[J]. 林业科学, 2019, 55(2): 97−108. doi: 10.11707/j.1001-7488.20190210 Xue C Q, Xu Q H, Lin L P, et al. Biomass models with breast height diameter and age for main native tree Species in Guangdong Province[J]. Scientia Silvae Sinicae, 2019, 55(2): 97−108. doi: 10.11707/j.1001-7488.20190210
[14] 欧建德, 欧家琳, 康永武. 基于树冠形态特征因子的云南红豆杉单木生物量模型拟合[J]. 西南林业大学学报(自然科学), 2022, 42(4): 117−124. Ou J D, Ou J L, Kang Y W. Single tree biomass simulation of Taxus yunnanensis plantation based on crown morphological index[J]. Journal of Southwest Forestry University (Natural Science), 2022, 42(4): 117−124.
[15] 符利勇, 雷渊才, 曾伟生. 几种相容性生物量模型及估计方法的比较[J]. 林业科学, 2014, 50(6): 42−54. Fu L Y, Lei Y C, Zeng W S. Comparison of several compatible biomass models and estimation approaches[J]. Scientia Silvae Sinicae, 2014, 50(6): 42−54.
[16] Fu L Y, Lei Y C, Wang G X, et al. Comparison of seemingly unrelated regressions with error in-variable models for developing a system of nonlinear additive biomass equations[J]. Trees, 2015, 30: 839−857.
[17] Zeng W S, Zhang L J, Chen X Y, et al. Construction of compatible and additive individual-tree biomass models for Pinus tabulaeformis in China[J]. Canadian Journal of Forest Research, 2017, 47: 467−475. doi: 10.1139/cjfr-2016-0342
[18] Dong L, Zhang Y, Zhang Z, et al. Comparison of tree biomass modeling approaches for larch (Larix olgensis Henry) trees in Northeast China[J]. Forests, 2020, 11: 202. doi: 10.3390/f11020202
[19] 孟延山, 孟俐君, 王静洁, 等. 青海省2种主要树种的生物量分配格局和单木生物量模型[J]. 西部林业科学, 2019, 48(6): 21−28. Meng Y S, Meng L J, Wang J J, et al. Models for estimating biomass and its allocation patterns in organs of two major tree species in Qinghai Province[J]. Journal of West China Forestry Science, 2019, 48(6): 21−28.
[20] 蔡会德, 卢峰, 徐占勇, 等. 桉树相容性可加性立木生物量模型系统研建[J]. 林业资源管理, 2023(1): 87−93. Cai H D, Lu F, Xu Z Y, et al. Research and development of compatible and additive individual tree biomass model systems for eucalyptus[J]. Forest Resources Management, 2023(1): 87−93.
[21] 梁瑞婷, 王轶夫, 邱思玉, 等. 人工神经网络与相容性生物量模型预测单木地上生物量的比较[J]. 应用生态学报, 2022, 33(1): 9−16. Liang R T, Wang Y F, Qiu S Y, et al. Comparison of artificial neural network with compatible biomass model for predicting aboveground biomass of individual tree[J]. Chinese Journal of Applied Ecology, 2022, 33(1): 9−16.
[22] Zhu G , Hu S , Chhin S , et al. Modelling site index of Chinese fir plantations using a random effects model across regional site types in Hunan Province, China[J]. Forest Ecology and Management, 2019, 446: 143−150.
[23] Parresol B R. Assessing tree and stand biomass: a review with examples and, critical comparisons[J]. Forest Science, 1999, 45(4): 573−593. doi: 10.1093/forestscience/45.4.573
[24] 曾伟生. 全国立木生物量方程建模方法研究[D]. 北京:中国林业科学研究院, 2011. Zeng W S. Methodology on modeling of single-tree biomass equations for national biomass estimation in China[D]. Beijing: Chinese Academy of Forestry, 2011.
[25] R Core Team. R: A language and environment for statistical computing[M].Vienna: R Foundation for Statistical Computing, 2024.
[26] 国家林业局. 立木生物量建模方法技术规程: LY/T 2258-2014[S]. 北京: 中国标准出版社, 2014. State Forestry Administration. Technical regulations for tree biomass modeling: LY/T 2258-2014[S]. Beijing: Standards Press of China, 2014.
[27] 欧光龙, 胥辉. 森林生物量模型研究综述[J]. 西南林业大学学报(自然科学), 2020, 40(1): 1−10. Ou G L, Xu H. A review on forest biomass models[J]. Journal of Southwest Forestry University, 2020, 40(1): 1−10.
[28] 曾伟生, 唐守正. 利用度量误差模型方法建立相容性立木生物量方程系统[J]. 林业科学研究, 2010, 23(6): 797−803. Zeng W S, Tang S Z. Using measurement error modeling method to establish compatible single-tree biomass equations system[J]. Forest Research, 2010, 23(6): 797−803.
[29] 刘秀红, 姜春前, 徐睿, 等. 相容性单木生物量模型估计方法的比较-以青冈栎为例[J]. 林业科学, 2020, 56(9): 164−173. Liu X H, Jiang C Q, Xu R, et al. Comparison of methods to construct compatible individual tree biomass models: a case study of Cyclobalanopsis glauca[J]. Scientia Silvae Sinicae, 2020, 56(9): 164−173.
[30] Wang X P, Fang J Y, Zhu B. Forest biomass and root shoot allocation in northeast China[J]. Forest Ecology and Management, 2008, 255: 4007−4020. doi: 10.1016/j.foreco.2008.03.055
[31] Luo Y J, Wang X K, Zhang X Q. Root: shoot ratios across China’s forests: Forest type and climatic effects[J]. Forest Ecology and Management, 2012, 269: 19−25. doi: 10.1016/j.foreco.2012.01.005
[32] Widagdo F, Li F, Zhang L, et al. Aggregated biomass model systems and carbon concentration variations for tree carbon quantification of natural Mongolian oak in Northeast China[J]. Forests, 2020, 11: 397. doi: 10.3390/f11040397
计量
- 文章访问数: 33
- HTML全文浏览量: 3
- PDF下载量: 7
