Applicability of the relationship between stand DBH and height to evaluate the site quality of major coniferous stands in China
-
摘要:目的
建立科学的立地质量评价体系，指导林业生产实践。
方法根据全国主要针叶林分样地数量情况，划分16个针叶树种组，采用Richards、Logistic和Korf模型拟合导线曲线，构建立地形指数模型，并进行落点检验，运用1999—2018年连续4期清查数据进行立地形等级动态变化分析。
结果Richards、Logistic和Korf模型拟合导向曲线决定系数均值均大于0.95，建立的立地形指数模型落点检验值均大于90.00%，落点检验均值达96.59%，适用于实际生产。20年间，针叶林分Ⅰ级和Ⅱ级均值合计增长了7.60个百分点，Ⅲ级均值合计减少了3.50个百分点，Ⅳ级和Ⅴ级均值合计减少了4.10个百分点，立地质量表现为较好的改善趋势。
结论基于胸径−树高关系建立全国统一的立地质量评价模型具有可行性和合理性，通过减少气候差异导致基于树龄的生长速率对立地质量评价的影响偏差，使不同地区间相同林分的评价结果具有可比性，在大尺度水平具有较好的适用性，但仍然需要警惕经营措施和小样本数据导致的评价结果不确定。
Abstract:ObjectiveTo establish a scientific site quality evaluation system, guide forestry production practice.
MethodBased on the number of major coniferous forest stand plots across the country, 16 coniferous tree groups were classified. The Richards, Logistic, and Korf models were applied to fit the guiding curves and establish a site form index model, which was then validated through a scatter plot test. Data from four consecutive forest inventory periods from 1999 to 2018 were used to analyze dynamic changes in site form levels.
ResultThe mean coefficients of determination for guiding curves fitted with the Richards, Logistic, and Korf models were all above 0.95. The scatter plot validation of the established site form index model showed test values exceeding 90.00%, with an average value of 96.59%, indicating its feasibility for practical use. Over 20 years, site quality for coniferous stands improved, with a combined proportion increase of 7.60 percent in grades I and II, a 3.50 percent decrease in grade III, and a combined reduction of 4.10 percent in grades IV and V.
ConclusionDeveloping a national unified site quality assessment model based on DBH-height relationships is feasible and reasonable. This approach reduces the impact deviations of climate-induced and age-based growth rate on site quality assessment, allowing for comparable site quality evaluations across regions. The model demonstrates good applicability on a large scale, though it is essential to consider potential uncertainties due to management practices and limited sample data.
-
Keywords:
- Site quality /
- DBH-height /
- Dynamic change /
- Coniferous forest stand
-
猪格拉瑟病是由副猪格拉瑟菌Glaesserella parasuis引起的以猪的多发性浆膜炎、脑膜炎、关节炎、心包炎为特征的细菌性传染病[1]。G. parasuis可感染断奶至屠宰任何阶段的猪群，并且与多种细菌、病毒、支原体等病原的混合感染广泛发生，给养猪业带来严重损失[2-3]。G. parasuis具有很强的基因变异性[4]，目前已知有15种血清型，不同血清型在基因型构成、毒力方面差异较大，这极大地制约了人们对G. parasuis感染的控制[5-7]，使得该病成为世界范围内引起仔猪死亡的重要因素之一。在临床治疗中，抗生素的广泛使用加速了G. parasuis耐药性的产生，为该病的防控带来巨大挑战。疫苗免疫是G. parasuis防控的主要手段，当前中国市场上的G. parasuis可用商品化疫苗主要是本地分离株灭活苗类。但灭活苗的保护效果有限，不同血清型之间缺少交叉保护作用[8]，而不同猪群流行的G. parasuis血清型可能不同，甚至同一头病猪的不同部位也可能存在不止一种血清型[9]。因此，研究具有交叉保护力的新型疫苗成为G. parasuis防控研究的重点[10]。
磷酸甘油醛脱氢酶(Glyceraldehyde-3-phosphate dehydrogenase，GAPDH)是广泛存在于真核生物以及原核生物中的一种酶，可与免疫细胞表面的受体结合激活相关信号通路。研究发现，GAPDH是G. parasuis的新型免疫原[11]，Fu等[12]通过基因测序与生物信息学分析发现，GAPDH广泛存在于各种血清型的G. parasuis中并且高度保守。但GAPDH作为G. parasuis普遍存在的免疫原性蛋白，却较少有针对该蛋白重组疫苗株的研究。外膜蛋白(Outer membrane proteins，OMPs)是细菌表面重要的抗原成分，能够参与决定宿主的特异性免疫应答。有学者发现在针对G. parasuis的免疫应答中存在OMPs的抗体，但没有脂多糖等成分的抗体，这表明OMPs比该生物体其他成分或许更具有免疫原性[13]。多个团队的研究者[14-16]从G. parasuis基因组中筛选出不同的OMPs作为潜在毒力相关因素，并证实其在豚鼠、小鼠模型中对致死量G. parasuis的攻击具有显著的保护作用。已有研究[17]显示，编码G. parasuis GAPDH和OmP26蛋白的基因在G. parasuis不同血清型菌株中普遍存在，且能够介导免疫反应，这意味着它们有可能成为广谱的免疫原。
理想的疫苗载体应具备多种优良特性，包括安全性、免疫原性等。赵战勤等[18]通过基因工程方法构建猪霍乱沙门氏菌Salmonella choleraesuis asd−缺失株C500Δasd，并证明该缺失株相较亲本菌株C500毒力降低；同时具有良好的生长特性、稳定的遗传背景及高效的表达能力。随后，徐引弟等[19]和赵战勤等[20]对C500Δasd缺失株展开更进一步研究，证实在C500Δasd缺失株中插入不同外源抗原如支气管败血波氏杆菌免疫原性基因、巴氏杆菌毒素基因后重组菌株仍具有良好的遗传稳定性，并高效表达外源抗原，将重组菌株通过口服免疫小鼠和仔猪，外源抗原有良好的表达活性，且接种动物无任何不良反应。上述研究成果表明，C500Δasd缺失株是构建重组菌株的有力候选载体。
因此，本研究选择GAPDH基因和OmP26基因作为免疫原性基因，将上述基因同时导入S. choleraesuis C500Δasd缺失株，以求构建能够高效表达G. parasuis免疫基因，同时能有效刺激机体产生针对S. choleraesuis特异性免疫保护作用的重组菌株，并对其基本生物学特性及免疫效果进行评估；旨在探索一种更为有效的猪格拉瑟病防控策略，同时为G. parasuis和S. choleraesuis混合感染的防控提供新的思路和方法，以应对这2种疾病给养猪业带来的挑战。
1. 材料与方法
1.1 主要试剂和试验动物
S. choleraesuis C500弱毒疫苗株的asd基因缺失株、大肠埃希菌Escherichia coli x6097及无抗性原核表达质粒pYA3493由华南农业大学兽医学院动物传染病教研室保存；限制性内切酶EcoR I + Sal I、Sal I + Hind III等各种内切酶购自New England Biolabs公司；T4 DNA连接酶购自TaKaRa公司；His标签蛋白纯化试剂盒(P2226)购自碧云天公司；4周龄雌性BALB/c小鼠购自广东省实验动物中心。
1.2 引物设计和基因克隆
参考Genebank中G. parasuis血清5型SH02165菌株的膜蛋白基因GAPDH(GenBank ID：
72785966 )和外膜蛋白基因OmP26(GenBank ID：219691637 )的序列，利用Primer 5.1设计引物(GAPDH-F：CGGAATTCATGGCAATTAAAATTG，GAPDH-R：GCGTCGACTTAGCCTTTGTAGTTG；OmP26-F：GCGTCGACATGAAAAATTTATTTAAACTTGC，OmP26-R：CCAAGCTTTTATTTTTTCACTTCTTCTGG)，并以G. parasuis SH02165菌株的基因组为模板扩增引物。1.3 重组质粒的构建
分别利用限制性内切酶EcoR I + Sal I、Sal I + Hind III对“1.2”回收获得的基因片段GAPDH、OmP26进行双酶切，原核表达质粒pYA3493使用对应限制性内切酶进行双酶切，然后经琼脂糖凝胶电泳回收GAPDH、OmP26基因片段与相应的线性质粒pYA3493，以T4 DNA连接酶进行连接。将质粒转入x6097，涂布在无抗性培养基上，挑取单菌落，在LB液体培养基中于37 ℃培养12 h后进行质粒小提，以酶切鉴定方法获取阳性克隆pYA-GAPDH和pYA-OmP26。用Sal I + Hind III双酶切基因片段OmP26与pYA-GAPDH载体，回收、纯化、连接后转入x6097，重复上述步骤，获得pYA-GAPDH-OmP26质粒，并将其送至上海生工生物工程有限公司进行测序鉴定。
1.4 重组菌株的构建
将重组质粒pYA-GAPDH、pYA-OmP26和pYA-GAPDH-OmP26电转入C501，完成电击后迅速加入37 ℃预热培养基，并将菌液转移到无菌管中；于37 ℃、120 r/min震荡1 h培养后离心收集菌体，重悬于单纯培养基并涂布于LB平皿表面。使用目的基因引物对平皿表面生长的菌落进行扩增鉴定。
1.5 重组菌株的遗传与生长稳定性
重组菌株C501(pYA-GAPDH)、C501(pYA-OmP26)及C501(pYA-GAPDH-OmP26)在无抗性LB液体培养基中连续传代培养，每25代进行PCR鉴定，测定其遗传特性。于37 ℃条件下培养重组菌株，每隔1 h测定D600 nm，以菌液培养时间为横坐标，菌液D600 nm为纵坐标绘制重组菌株的生长曲线。
1.6 重组菌株的生化特性与表达特性
将重组菌株C501(pYA-GAPDH)、C501(pYA-OmP26)及C501(pYA-GAPDH-OmP26)培养于无抗性LB液体培养基中，转接至以阿拉伯糖、乳糖、棉子糖、山梨醇、淀粉、半乳糖、葡萄糖、甘露糖、果糖、鼠李糖和麦芽糖等为唯一碳源的细菌微量生化鉴定管进行生化试验。对重组菌株全菌裂解物进行Western-Blot鉴定。
1.7 重组菌株C501(pYA-GAPDH-OmP26)对小鼠的免疫试验
选取40只4周龄雌性BALB/c小鼠，分3组：C501(pYA-GAPDH-OmP26)组(20只小鼠)、C501(pYA3493)组(10只小鼠)、PBS阴性对照组(10只小鼠)，按每只5×108 CFU/mL剂量皮下多点位注射。首次免疫后第14天进行加强免疫，操作程序同首次免疫。加强免疫14 d后各组分别按照每只2.5×109 CFU/mL剂量腹腔攻毒G. parasuis 5型强毒株SH0165和每只5×107 CFU/mL剂量灌胃攻毒S. choleraesuis强毒株C78-1，各小组中2种菌株攻毒相同数量小鼠，连续观察14 d，记录发病及死亡情况。
2. 结果与分析
2.1 重组菌株的构建
目的基因GAPDH和OmP26扩增结果分别为1 020、798 bp，与预期片段长度一致，测序结果正确(图1)。重组质粒pYA-GAPDH、pYA-OmP26酶切鉴定均能成功扩增出约3 000 bp的外源基因片段(图2)。pYA-GAPDH、pYA-OmP26和pYA-GAPDH-OmP26测序结果均符合预期。重组菌株C501(pYA-GAPDH)、C501(pYA-OmP26)及C501(pYA-GAPDH-OmP26)用GAPDH、OmP26引物扩增，测序结果与目的外源基因一致。
![]()
图 1 目的基因PCR扩增片段Figure 1. PCR amplified fragments of target genesM: DNA Marker, 1: GAPDH, 2: OmP26.![]()
图 2 重组质粒pYA-GAPDH和pYA-OmP26的酶切鉴定Figure 2. Enzymatic identification of recombinant plasmids pYA-GAPDH and pYA-OmP26M: DNA Marker, 1: pYA-GAPDH, 2: pYA-OmP26.2.2 重组菌株遗传与生长稳定性
重组菌株C501(pYA-GAPDH)、C501(pYA-OmP26)及C501(pYA-GAPDH-OmP26)连续传代100代，每25代进行1次PCR鉴定，结果显示，不同代次均能扩增出外源基因条带(图3)，表明重组质粒pYA-GAPDH、pYA-OmP26和pYA-GAPDH-OmP26均能在宿主菌菌株C501中稳定遗传。重组菌株的生长曲线与对照菌株C501(pYA3493)非常接近(图4)，表明外源基因基本不影响宿主菌株生长趋势。
![]()
图 3 重组菌株遗传稳定性M：DNA Marker；1~4：C501(pYA-GAPDH) 25、50、75、100代；6~9：C501(pYA-OmP26) 25、50、75、100代；11~14：C501(pYA-GAPDH-OmP26) 25、50、75、100代；5、10：pYA3493。Figure 3. Genetic stability of recombinant strainsM: DNA Marker; 1−4: C501(pYA-GAPDH) 25, 50, 75, 100 generations; 6−9: C501(pYA-OmP26) 25, 50, 75, 100 generations; 11−14: C501(pYA-GAPDH-OmP26) 25, 50, 75, 100 generations; 5, 10: pYA3493.![]()
图 4 重组菌株C501(pYA-GAPDH-OmP26)与对照菌株C501(pYA3493)的生长曲线Figure 4. Growth curves of recombinant strain C501 (pYA-GAPDH-OmP26) and control strain C501 (pYA3493)2.3 重组菌株生化特性
重组菌株C501(pYA-GAPDH)、C501(pYA-OmP26)及C501(pYA-GAPDH-OmP26)与对照菌株C501(pYA3493)均不能利用阿拉伯糖、乳糖、棉子糖、山梨醇、淀粉和半乳糖，但可以利用葡萄糖、甘露醇、果糖、麦芽糖和鼠李糖(表1)，表明外源片段引入不影响Salmonella自身的生化特性。
表 1 重组沙门氏菌对不同碳源的利用情况1)Table 1. Utilization of different carbon sources by recombinant Salmonella spp.菌株
Strain阿拉伯糖Arabinose 乳糖Lactose 棉子糖Raffinose 山梨醇Sorbitol 淀粉Starch 半乳糖Galactose 葡萄糖Glucose 甘露醇Mannitol 果糖Fructose 鼠李糖Rhamnose 麦芽糖Maltose C501(pYA3493) − − − − − − + + + + + C501(pYA-GAPDH) − − − − − − + + + + + C501(pYA-OmP26) − − − − − − + + + + + C501(pYA-GAPDH-OmP26) − − − − − − + + + + + 1) −：阴性；+：阳性。
1) −: Negative; +: Positive.2.4 重组菌株表达特性
重组菌株C501(pYA-GAPDH)、C501(pYA-OmP26)全菌裂解物Western-Blot结果显示，与空质粒对照相比，重组疫苗菌株能够正确表达外源蛋白，蛋白相对分子质量分别为37 400、37 000(图5)。
![]()
图 5 重组菌株Western-Blot分析Mr：相对分子质量，M：Marker，1：C501(pYA3493)，2：C501(pYA-GAPDH)，3：C501(pYA-OmP26)。Figure 5. Western-Blot analysis of recombinant strainsMr: Relative molecular mass, M: Marker, 1: C501(pYA3493), 2: C501(pYA-GAPDH), 3: C501(pYA-OmP26).2.5 免疫小鼠攻毒保护试验
C501(pYA-GAPDH-OmP26)组免疫小鼠对S. choleraesuis强毒株C78-1攻毒的保护率为62.5%，对G. parasuis 5型强毒株SH0165攻毒的保护率为50.0%；C501(pYA3493)组、PBS对照组对G. parasuis 5型强毒株SH0165攻毒的保护率均为0(图6)。
3. 讨论与结论
3.1 讨论
接种疫苗是预防猪格拉瑟病的有效措施。但现有的商业疫苗对异源菌株的交叉保护作用有限，不能提供广泛的异源保护[21]。如基于血清型2和5组合、血清型1和6组合的市售疫苗，对其他血清型甚至同一血清型不同菌株的保护力常常无法达到预期[22]。因此，研发对不同血清型具有交叉保护作用的新型疫苗对该病的防治具有重要意义。诸多研究表明，S. choleraesuis C500Δasd缺失株作为活疫苗载体能够同时表达多种不同病原的抗原，且具有极高的安全性和可行性，作为胞内侵袭细菌，Salmonella能够有效递呈抗原，激发机体产生抗Salmonella免疫反应，同时诱导产生针对外源基因的特异性细胞免疫和体液免疫[23-24]。S. choleraesuisC500Δasd缺失株导入含asd质粒后才能在无外源二氨基庚二酸条件下存活，这使得质粒携带的外源基因不易丢失，同时避免Salmonella过度繁殖[25-26]。本研究选取的外源抗原OmP26是G. parasuis的外膜蛋白，本身具有良好的免疫活性；GAPDH作为G. parasuis多种血清型普遍存在的蛋白，具有免疫调节功能。将二者结合同时克隆，可增加机体可识别的抗原位点，从而提高机体免疫应答的效率。同时，GAPDH和OmP26二者同时表达可能存在协同刺激免疫系统的效果，产生更强的免疫反应。此外，同时包含GAPDH和OmP26基因的表达载体对于血清型众多的G. parasuis而言，能够提供更广泛的保护。
基于此，本研究从G. parasuis中筛选出具有良好免疫效应的蛋白GAPDH、OmP26，将其同时插入S. choleraesuis C500Δasd缺失株，成功构建二联基因工程弱毒株C501(pYA-GAPDH-OmP26)。试验结果显示，C501(pYA-GAPDH-OmP26)能稳定携带、表达外源基因，并与对照菌株C501(pYA3493)具有相近的生物学特性，表明外源片段的插入基本不影响宿主菌株的生长与代谢。通过皮下注射C501(pYA-GAPDH-OmP26)免疫的小鼠对致死剂量的G. parasuis强毒株表现出50.0%的保护率，与对照组C501(pYA3493)组、PBS组(保护率为0)均存在明显差异。同时，C501(pYA-GAPDH-OmP26)、C501(pYA3493)对S. choleraesuis C78-1的保护率分别为62.5%、50.0%，说明外源GAPDH和OmP26抗原的表达没有降低C500Δasd缺失株针对Salmonella感染的免疫保护力。但该重组菌株在小鼠与在猪体内引发的免疫反应可能存在差异，目前针对G. parasuis重组疫苗菌株的研究也大都在小鼠试验模型上进行，少有在猪模型上检验[27]。在后续研究中，有必要在本体动物仔猪动物模型中进一步评估该重组菌株的免疫效力。
3.2 结论
本研究选取G. parasuis的GAPDH和OmP26为外源性抗原并成功研制G. parasuis-S. choleraesuis二联基因工程弱毒株C501(pYA-GAPDH-OmP26)，与空质粒组C501(pYA3493)相比，该重组菌株对G. parasuis 5型强毒株SH0165的保护效果更好且存在明显差异，同时保留了对S. choleraesuis C78-1的高效保护，为进一步研发G. parasuis新型疫苗提供了参考。
-
![]()
图 1 主要针叶林分胸径和树高拟合曲线图
Figure 1. Fitted curve plots of DBH and height for major coniferous stands
![]()
图 2 主要针叶林分立地形指数模型曲线簇
Figure 2. Curve families of site form index model for major coniferous stands
![]()
图 3 2003、2008、2013和2018年主要针叶林分立地形等级占比动态直方图
Figure 3. Proportion dynamic histogram of site form class for major coniferous stands in 2003, 2008, 2013 and 2018
表 1 立地形指数模型研建数据描述性统计
Table 1 Descriptive statistics of model establishment data for site form index
林分类型 Stand types 样地数量
Number of plots林分胸径/cm Stand DBH 林分树高/m Stand height 区间 Range 均值 Mean 标准差 SD 区间 Range 均值 Mean 标准差 SD 冷杉 Abies fabri 4363 6.5~77.1 32.8 11.05 2.8~43.7 19.8 5.95 云杉 Picea asperata 11869 5.0~60.0 27.2 10.61 1.7~41.5 17.0 6.11 铁杉 Tsuga chinensis 182 6.9~58.0 29.9 9.61 4.0~29.0 16.6 4.42 油杉 Keteleeria fortunei 217 6.9~27.1 13.0 4.18 2.8~18.3 7.2 2.97 落叶松 Larix gmelinii 9105 5.0~69.7 16.4 9.19 2.0~36.2 13.5 5.00 红松 Pinus koraiensis 290 5.0~58.9 17.5 10.52 1.5~29.2 12.1 5.83 樟子松 Pinus sylvestris 563 5.4~43.2 16.9 7.13 2.8~26.7 11.6 4.78 赤松 Pinus densiflora 296 5.3~23.1 10.6 3.55 1.8~16.7 6.1 3.06 黑松 Pinus thunbergii 346 5.3~20.0 10.8 3.13 2.2~16.0 6.3 2.46 油松 Pinus tabuliformis 4525 5.0~35.0 13.2 5.16 1.5~23.0 7.7 3.27 华山松 Pinus armandii 1029 5.0~34.3 13.9 5.60 1.5~25.0 9.3 3.83 马尾松 Pinus massoniana 15430 5.0~39.9 12.9 4.90 1.5~28.5 9.2 3.58 云南松 Pinus yunnanensis 3510 5.0~43.0 14.6 6.58 2.2~30.0 9.6 4.48 思茅松 Pinus kesiya var. langbianensis 478 5.7~33.4 16.5 5.14 2.9~27.2 12.4 4.13 高山松 Pinus densata 5125 5.3~40.0 27.9 6.53 2.0~28.0 17.8 4.41 其他松类1) Other pines 1350 5~27.4 12.8 4.44 2.5~18.7 8.1 3.06 1) 其他松类指样地数量较少的针叶林分类型。
1) Other pines indicate stand types with less sample plots.
下载: 导出CSV
表 2 主要针叶林分胸径和树高模型参数
Table 2 Model parameters of DBH and height for major coniferous stands
林分
Stand typesRichards Logistic Korf a b c R2 SEE a b c R2 SEE a b c R2 SEE 冷杉 34.9304 0.0278 1.2032 0.94 2.08 29.7918 6.2168 − 0.0711 0.92 2.25 71.0181 − 8.2024 0.5228 0.94 2.07 云杉 35.0515 0.0292 1.2862 0.99 0.83 27.6166 7.9591 − 0.0872 0.98 1.12 104.1793 − 7.7122 0.4304 0.99 0.82 铁杉 42.7731 0.0100 0.7756 0.84 2.36 25.0287 4.3132 − 0.0612 0.81 2.52 228.4099 − 6.1365 0.2395 0.84 2.35 油杉 49.4001 0.0217 1.5328 0.96 0.94 17.6832 15.4880 − 0.1514 0.95 0.96 356.5069 − 9.6653 0.3327 0.96 0.94 落叶松 21.7763 0.0825 1.6565 0.92 1.54 21.4231 5.7592 − 0.1286 0.91 1.63 24.4672 − 16.0920 1.1943 0.91 1.63 红松 23.7699 0.0706 1.9944 0.95 1.54 22.4409 9.9858 − 0.1324 0.95 1.64 31.3234 − 15.3821 0.9800 0.95 1.60 樟子松 23.1622 0.0573 1.6375 0.97 1.07 19.8216 9.2197 − 0.1350 0.96 1.27 69.4325 − 7.1841 0.4718 0.97 1.06 赤松 11.2997 0.1945 5.7997 0.89 1.32 10.6562 48.3629 − 0.3444 0.92 1.14 15.2157 − 32.2581 1.4462 0.88 1.43 黑松 32.0888 0.0315 1.4875 0.97 0.60 13.1164 14.9074 − 0.2020 0.98 0.49 646.3492 − 8.9301 0.2569 0.97 0.62 油松 34.1642 0.0217 1.2051 0.99 0.28 17.6943 9.4751 − 0.1219 0.99 0.49 474.2441 − 8.0463 0.2431 0.99 0.26 华山松 23.2458 0.0527 1.5916 0.99 0.59 17.9922 10.4218 − 0.1496 0.98 0.68 69.1080 − 7.7007 0.4897 0.99 0.59 马尾松 20.7730 0.0585 1.4939 0.99 0.28 17.5935 8.1581 − 0.1427 0.99 0.55 41.5400 − 7.4772 0.5971 0.99 0.24 云南松 32.6336 0.0345 1.4617 0.99 0.38 22.9034 10.6239 − 0.1183 0.99 0.67 182.4129 − 7.9799 0.3567 0.99 0.38 思茅松 23.3830 0.0688 1.8721 0.97 0.87 19.5171 10.3107 − 0.1589 0.96 1.01 42.3634 − 10.2905 0.7343 0.98 0.84 高山松 47.8737 0.0196 1.2483 0.99 0.71 25.0354 10.2222 − 0.1076 0.98 0.89 378.0811 − 8.2583 0.2902 0.99 0.71 其他松类1) 15.6588 0.0804 1.7907 0.99 0.35 13.0084 10.8083 − 0.1960 0.99 0.28 36.7931 − 7.3685 0.5885 0.99 0.40 1) 其他松类指样地数量较少的针叶林分类型。
1) Other pines indicate stand types with less sample plots.
下载: 导出CSV
-
[1] WEISKITTEL A R, HANN D W, KERSHAW J A, et al. Forest growth and yield modeling[M]. New York: John Wiley & Sons, 2011.
[2] AGUIRRE A, MORENO-FERN´ANDEZ D, ALBERDI I, et al. Mapping forest site quality at national level[J]. Forest Ecology and Management, 2022, 508: 120043. doi: 10.1016/jforece.2022.120043.
[3] 夏洪涛, 郭晓斌, 张珍, 等. 基于不同立地质量评价指标的杉木大径材林分树高−胸径模型[J]. 中南林业科技大学学报, 2023, 43(10): 80-88. [4] 李贤伟, 李守剑, 张健, 等. 四川盆周西缘水杉人工林林地立地质量评价研究[J]. 四川农业大学学报, 2002, 20(2): 106-109. doi: 10.3969/j.issn.1000-2650.2002.02.010 [5] 刘智军, 朱丽艳, 吴恒, 等. 基于清查数据的立地质量评价方法准确性分析[J]. 西南林业大学学报(自然科学), 2019, 39(5): 127-135. [6] VANCLAY J K, HENRY N B. Assessing site productivity of indigenous cypress pine forest in southern Queensland[J]. The Commonwealth Forestry Review, 1988, 67(1): 53-64.
[7] 吴恒, 党坤良, 田相林, 等. 秦岭林区天然次生林与人工林立地质量评价[J]. 林业科学, 2015, 51(4): 78-88. [8] 沈剑波, 雷相东, 雷渊才, 等. 长白落叶松人工林地位指数及立地形的比较研究[J]. 北京林业大学学报, 2018, 40(6): 1-8. [9] PADILLA-MARTÍNEZ J R, PAUL C, CORRAL-RIVAS J J, et al. Evaluation of the site form as a site productive indicator in temperate uneven-aged multispecies forests in Durango, Mexico[J]. Plants, 2022, 11(20): 2764. doi: 10.3990/plants11202764.
[10] 张会儒, 雷相东, 张春雨, 等. 森林质量评价及精准提升理论与技术研究[J]. 北京林业大学学报, 2019, 41(5): 1-18. [11] 胡庭兴, 李贤伟, 杨祯禄. 立地质量综合评价及其应用的研究[J]. 四川农业大学学报, 1993,11(3): 397-403. [12] PERIN J, HÉBERT J, BROSTAUX Y, et al. Modeling the top-height growth and site index of Norway spruce in Southern Belgium[J]. Forest Ecology and Management, 2013, 298: 62-70.
[13] 付晓, 曹霖, 王雪军, 等. 吉林省林地立地质量评价及生产潜力研究[J]. 中南林业科技大学学报, 2019, 39(5): 1-9. [14] 黄国胜, 马炜, 王雪军, 等. 基于一类清查数据的福建省立地质量评价技术[J]. 北京林业大学学报, 2014, 36(3): 1-8. [15] MACFARLANE D W, GREEN E J, BURKHART H E. Population density influences assessment and application of site index[J]. Canadian Journal of Forest Research, 2000, 30(9): 1472-1475.
[16] 雷相东, 符利勇, 李海奎, 等. 基于林分潜在生长量的立地质量评价方法与应用[J]. 林业科学, 2018, 54(12): 116-126. doi: 10.11707/j.1001-7488.20181213 [17] SHARMA M, PARTON J. Height–diameter equations for boreal tree species in Ontario using a mixed-effects modeling approach[J]. Forest Ecology and Management, 2007, 249(3): 187-198.
[18] KROON J, ANDERSSON B, MULLIN T J. Genetic variation in the diameter-height relationship in Scots pine (Pinus sylvestris)[J]. Canadian Journal of Forest Research, 2008, 38(6): 1493-1503.
[19] 王彬. 祁连山青海云杉树木生长及其对气候变化的响应[D]. 北京: 中国林业科学研究院, 2020. [20] MOTALLEBI A, KANGUR A. Are allometric relationships between tree height and diameter dependent on environmental conditions and management[J]. Trees, 2016, 30(4): 1429-1443.
[21] 郭小阳, 吴恒, 田相林, 等. 基于优势高模型分析多源数据对立地质量评价的影响[J]. 西北林学院学报, 2017, 32(6): 184-189. doi: 10.3969/j.issn.1001-7461.2017.06.28 [22] MOLINA-VALERO J A, DIÉGUEZ-ARANDA1 U, ÁLVAREZ-GONZÁLEZ1 J G, et al. Assessing site form as an indicator of site quality in even-aged Pinus radiata D. Don stands in north-western Spain[J]. Annals of Forest Science, 2019, 76(4): 113. doi: 10.1007/s13595-019-0904-1
[23] 刘俊, 孟雪, 高媛赟, 等. 基于不同立地质量的松树林分生物量遥感估测[J]. 中南林业科技大学学报, 2016, 36(5): 41-46. [24] GARCIA O. Comparing and combining stem analysis and permanent sample plot data in site index models[J]. Forest Science, 2005, 51(4): 277-283. doi: 10.1093/forestscience/51.4.277
计量
- 文章访问数: 25
- HTML全文浏览量: 7
- PDF下载量: 7
