Processing math: 6%
  • Scopus收录期刊
  • CSCD(核心库)来源期刊
  • 中文核心期刊
  • 中国科技核心期刊
  • F5000顶尖学术来源期刊
  • RCCSE中国核心学术期刊
高级检索

84份核桃种质资源萌芽期和雌雄花期调查

刘警, 宋雅萍, 于秋香

刘警, 宋雅萍, 于秋香. 84份核桃种质资源萌芽期和雌雄花期调查[J]. 北京林业大学学报, 2025, 47(2): 49-57. DOI: 10.12171/j.1000-1522.20240271
引用本文: 刘警, 宋雅萍, 于秋香. 84份核桃种质资源萌芽期和雌雄花期调查[J]. 北京林业大学学报, 2025, 47(2): 49-57. DOI: 10.12171/j.1000-1522.20240271
Liu Jing, Song Yaping, Yu Qiuxiang. Investigation on germination period and male and female flowering period of 84 walnut germplasm resources[J]. Journal of Beijing Forestry University, 2025, 47(2): 49-57. DOI: 10.12171/j.1000-1522.20240271
Citation: Liu Jing, Song Yaping, Yu Qiuxiang. Investigation on germination period and male and female flowering period of 84 walnut germplasm resources[J]. Journal of Beijing Forestry University, 2025, 47(2): 49-57. DOI: 10.12171/j.1000-1522.20240271

84份核桃种质资源萌芽期和雌雄花期调查

基金项目: 河北省农林科学院创新工程课题(2022KJCXZX-CGS-9)。
详细信息
    作者简介:

    刘警,助理研究员。主要研究方向:核桃种质资源搜集及分子辅助育种。Email:1550957472@qq.com 地址:066600河北省秦皇岛市昌黎县碣阳大街东段39号昌黎果树研究所

    责任作者:

    于秋香,副研究员。主要研究方向:核桃种质资源创新及丰产栽培技术研究。Email:75708282@qq.com 地址:同上。

  • 中图分类号: S664.1

Investigation on germination period and male and female flowering period of 84 walnut germplasm resources

  • 摘要:
    目的 

    针对核桃生产中春季频发晚霜和雌雄花期不相遇等问题,探究河北省秦皇岛地区84份核桃种质资源萌芽期、雌雄花盛期起始时间,旨在为选育避晚霜品种、合理配置授粉树、提高核桃产量提供参考。

    方法 

    2024年对河北省秦皇岛市昌黎县84份核桃种质资源进行萌芽期、雌雄花盛期起始时间调查,并记录此时期温度,计算萌芽和开花期起始所需≥10 ℃有效积温。

    结果 

    (1)84份核桃种质资源萌芽期起始时间为3月24日—3月31日,所需平均≥10 ℃有效积温为6.12 ℃·d,雄花盛期起始时间为4月17日—5月10日,所需平均≥10 ℃有效积温为103.69 ℃·d;雌花盛期起始时间为4月18日—5月15日,所需平均≥10 ℃有效积温为109.25 ℃·d。(2)84份核桃种质资源雌雄花期可划分为3大类:早花期(4月17日—4月26日),中花期(4月27日—5月6日),晚花期(5月7日及之后)。(3)在未来核桃育种工作中,可选择萌芽较晚的品种(‘中核3号’‘西洛7号’‘金薄香2号’‘金薄香4号’‘西林1号’‘绿香’‘晋香’‘花生1号’‘中核2号’等),以避开春季晚霜危害(≤0 ℃);可合理搭配早中晚花资源,‘礼品2号’(雌花早花)搭配‘西岭’(雄花早花),‘绿岭’(雌花中花)搭配‘西林1号’(雄花中花),‘小果黑核桃’(雌花晚花)搭配‘黑核桃’(雄花晚花),提高授粉率和坐果率。

    结论 

    本研究明确了84份核桃种质资源在河北省秦皇岛地区萌芽期和雌雄花盛期起始时间,研究结果为84份核桃种质资源的有效开发及利用,合理有效配置授粉品种,筛选和选育避晚霜品种,预测萌芽期和花期以提前采取防寒措施,产业增产增效提供了重要依据。

    Abstract:
    Objective 

    To solve frequent late frosts in spring and the mismatch of male and female flowering periods in walnut production, this study investigated the starting time of germination and flower period of 84 walnut germplasm resources in Qinhuangdao, Hebei Province of northern China, aiming to provide a reference for breeding late frost-resistant varieties, optimizing the arrangement of pollinator trees, and enhancing walnut yield.

    Method 

    In 2024, the starting time of germination and flowering period of 84 walnut germplasm resources in Changli County of Qinhuangdao City were investigated, and the temperature during this period were recorded. The effective accumulated temperature of ≥10 ℃ required for germination and flowering was calculated.

    Result 

    (1) The germination period of 84 walnut germplasm resources started from March 24 to March 31, requiring an average ≥10 ℃ effective accumulated temperature of 6.12 ℃·d; the peak period for male flowers started from April 17 to May 10, requiring an average ≥10 ℃ effective accumulated temperature of 103.69 ℃·d; the peak period for female flowers started from April 18 to May 15, requiring an average ≥10 ℃ effective accumulated temperature of 109.25 ℃·d. (2) The 84 germplasm resources were categorized into three groups: early-flowering group (April 17 − April 26), mid-flowering group (April 27 − May 6), and late-flowering group (May 7 and later). (3) In future walnut breeding work, we can select cultivars with late germination period such as ‘Zhonghe3hao’ ‘Xiluo7hao’ ‘Jinboxiang2hao’ ‘Jinboxiang4hao’ ‘Xilin1hao’ ‘Lvxiang’ ‘Jinxiang’ ‘Huasheng1hao’ ‘Zhonghe2hao’ to avoid the late frost in spring (≤0 ℃); and combine reasonably of cultivars with early, mid, and late flowers, such as ‘Lipin2hao’ (female early-flowering group) with ‘Xiling’ (male early-flowering group), ‘Lvling’ (female mid-flowering group) with ‘Xilin1hao’ (male mid-flowering group), ‘Small fruit black walnut’ (female late-flowering group) with ‘Black walnut’ (male late-flowering group) to improve the pollination rates and fruit setting rates.

    Conclusion 

    The findings of this research provide important basis for the effective development and utilization of 84 walnut germplasm resources, optimize the arrangement of pollinator varieties, screen and breed late frost-resistant varieties, predict germination and flower period to implement preemptive frost prevention measures, and promote increased production and efficiency in the walnut industry.

  • 树高和胸径是森林资源调查和经营效果评价的主要林分调查因子,也是林木生长发育和林分材积数量研究的重要变量。树高与胸径关系一直是林学研究的热点。通过量化特定树种的树高与胸径关系拟合树高−胸径预测模型,在计算树干材积、出材率表编制、生长和收获模型、森林资源调查、立地质量评价以及森林碳储量估算方面具有广泛应用。

    分位数回归方法具有较长的统计学历史和理论基础,最早于1978年由Koenker等[1]根据中位数回归推导提出。由于对调查数据模型方差没有严格要求,利用自变量的多个分位数求得因变量条件分布的相应分位数方程,能够描述任意发散点处变量的变化趋势。对具有尖峰、厚尾、异方差显著的数据拟合效果更加稳健[2],相对最小二乘回归具有更大优势。已广泛应用于林分自疏边界线、直径分布规律、林分密度指数、林木干形、树冠轮廓以及森林病虫害等领域。国外学者利用分位数回归建立了树高−胸径关系模型[3]。国内学者利用分位数回归拟合了树冠轮廓[4]、树冠表面积及树冠体积[5],结合哑变量方法构建了大兴安岭落叶松(Larix gmelinii)的树高与胸径关系最优模型[6]。Cao等[7]提出利用组合分位数法预测火炬松 (Pinus taeda)干形曲线,并通过插值系数修正削度方程以提高预测精度。

    树木生长受气候、立地条件和经营措施等多种因素的影响和制约,不同样地间的树高−胸径关系存在较大差异,仅含胸径一个变量的树高模型不能描述不同林况下树高−胸径关系的差异,因此很多引入林分变量或者单木变量构建的广义模型也被用于描述不同林分的模型差异。除此之外,近年来混合效应模型在树高模型中也得到了一定的应用,可以反映总体平均变化趋势及个体之间的变异,不需要假设观测值相互独立[8],显著提高了模型的预测精度。Zhang等[9]在以样地为随机效应的混合模型中引入林分密度和立地指数,结果表明该模型可以显著提高长白落叶松 (Larix olgensis)树高−胸径模型的预测精度。燕云飞等[10]基于广义混合效应构建了拟合效果最优的红松(Pinus koraiensis)人工林枝下高模型。

    以往树高−胸径关系研究均选定某种方程进行非线性回归,且对多层次、多树种的混交林树高−胸径关系的描述相比纯林要困难和复杂很多[11]。同时随着方程通用程度的提高,用于局部地区的误差也会增大[12]。因此,本研究以吉林蛟河针阔混交林典型树种为研究对象,构建含有林分变量的广义树高−胸径模型,并基于基础模型和广义模型分布拟合不同树种的分位数回归和以样方效应作为分类变量的混合模型,对比分析各树种不同树高−胸径模型的预测效果和拟合差异,为温带针阔混交林资源调查、收获、经营与管理提供依据。

    吉林省林业实验区国有林保护中心位于吉林省东部、长白山西麓,隶属于长白山系张广才岭山脉。行政区域处于吉林省蛟河市前进乡境内,属温带大陆性季风气候,年平均气温3.4 ℃,降水量708.8 mm,无霜期为120 ~ 130 d。地貌类型以低山丘陵区为主,海拔459 ~ 517 m。森林资源和植物种类丰富,植被类型为次生针阔混交林,主要乔木树种包括红松、水曲柳(Fraxinus mandshurica)、蒙古栎(Quercus mongolica)、紫椴(Tilia amurensis)、色木槭(Acer mono)、白牛槭(Acer mandshuricum)和胡桃楸(Juglans mandshurica)等。

    在研究区设置总面积大小为21.12 hm2监测样地,经纬度为43°57′35″ ~ 43°57′54″N,127°43′37″ ~ 127°44′09″E。2009年对植被开展初次调查,2014年、2019年对样地进行复测。调查时将样地划分为528个20 m × 20 m的连续调查样方,记录所有样方中胸径大于1 cm的木本植物种类,测量胸径、树高、冠幅、坐标、林分郁闭程度等。基于样地内树种组成(表1)和各树种树高−胸径模型的初步拟合结果,选择重要值最高的乔木树种作为研究样本,其中筛选出1 531株红松、3 066株色木槭、1 822株紫椴和1 856株水曲柳。建模时将各组数据随机分成两部分,其中80%作为训练数据,20%作为测试数据(表2)。根据调查数据,计算各样方的林分密度、胸高断面积和优势木高(表3)。在每个20 m × 20 m样方内,分别选择1 ~ 3株最高的林木计算平均值,作为该树种在该样方的优势木高。利用R软件的quantreg包和nlme包进行非线性分位数回归和混合效应模型分析。

    表  1  优势树种主要指标统计
    Table  1.  Main indicator statistics of dominant tree species
    树种
    Tree species
    重要值
    Important value
    多度
    Abundance
    频度
    Frequency
    胸高断面积/(m2·hm−2
    Basal area/(m2·ha−1
    色木槭 Acer mono 11.42 3 071 505 87.58
    水曲柳 Fraxinus mandshurica 11.13 1 861 420 119.81
    紫椴 Tilia amurensis 9.46 1 825 472 84.91
    红松 Pinus koraiensis 8.21 1 535 411 74.72
    暴马丁香 Syringa reticulata 7.28 3 483 466 4.66
    春榆 Ulmus davidiana 6.79 2 279 446 26.94
    拧筋槭 Acer triflorum 6.51 1 474 444 41.76
    蒙古栎 Quercus mongolica 5.62 646 337 55.80
    白牛槭 Acer mandshuricum 4.67 1 492 358 15.56
    白桦 Betula platyphylla 4.44 876 249 36.81
    下载: 导出CSV 
    | 显示表格
    表  2  单木树高−胸径模型拟合及检验数据
    Table  2.  Fitting and testing data of individual tree height-DBH model
    数据分类
    Dataset classification
    树种
    Tree species
    样本数
    Number of
    sample
    胸径 DBH/cm 树高 Tree height/m
    平均值
    Mean
    最小值
    Min. value
    最大值
    Max. value
    标准差
    SD
    平均值
    Mean
    最小值
    Min. value
    最大值
    Max. value
    标准差
    SD
    拟合数据
    Fitting data
    红松 Pinus koraiensis 1 226 19.14 1.60 75.80 15.90 12.25 1.00 33.40 6.75
    色木槭 Acer mono 2 452 14.13 1.20 88.80 12.62 11.22 2.00 37.20 6.10
    紫椴 Tilia amurensis 1 458 19.86 1.20 70.40 13.76 14.42 1.60 32.70 6.36
    水曲柳 Fraxinus mandshurica 1 484 26.17 2.00 75.00 12.32 17.82 2.70 27.30 4.15
    检验数据
    Testing data
    红松 Pinus koraiensis 305 19.28 2.00 65.30 16.12 12.11 1.90 34.40 6.98
    色木槭 Acer mono 614 14.32 1.40 62.00 12.55 11.49 1.20 29.40 6.42
    紫椴 Tilia amurensis 364 20.12 1.70 65.80 14.60 14.45 1.80 34.00 6.68
    水曲柳 Fraxinus mandshurica 372 24.84 1.80 62.40 11.76 17.75 4.50 27.00 4.21
    下载: 导出CSV 
    | 显示表格
    表  3  林分调查因子统计表
    Table  3.  Statistical table of stand description factors
    数据分类
    Dataset classification
    树种
    Tree species
    胸高断面积/(m2·hm−2
    Basal area/(m2·ha−1
    林分密度/(株·hm−2
    Stand density/(tree·ha−1
    优势木高
    Dominant tree height/m
    平均值
    Mean
    最小值
    Min.
    value
    最大值
    Max.
    value
    标准差
    SD
    平均值
    Mean
    最小值
    Min.
    value
    最大值
    Max.
    value
    标准差
    SD
    平均值
    Mean
    最小值
    Min.
    value
    最大值
    Max.
    value
    标准差
    SD
    拟合数据
    Fitting data
    红松 Pinus koraiensis 32.91 15.10 257.89 13.62 1 223.63 350 2 900 384.76 16.30 1.80 33.40 6.97
    色木槭 Acer mono 32.66 15.02 257.89 17.34 1 244.02 350 3 100 400.18 16.98 4.00 37.20 5.56
    紫椴 Tilia amurensis 34.00 15.02 257.89 14.38 1 284.93 350 2 900 387.05 19.62 3.40 34.00 4.79
    水曲柳 Fraxinus mandshurica 31.44 15.02 257.89 15.41 1 420.70 350 3 100 447.84 24.21 2.80 35.30 4.49
    检验数据
    Testing data
    红松 Pinus koraiensis 33.49 19.21 58.18 7.35 1 215.16 450 2 575 346.34 16.88 2.00 33.40 7.14
    色木槭 Acer mono 33.46 16.19 257.89 16.51 1 258.15 450 2 900 374.19 17.23 3.20 37.20 5.31
    紫椴 Tilia amurensis 33.82 15.02 257.89 16.92 1 310.16 600 2 900 412.81 19.67 6.20 34.00 5.15
    水曲柳 Fraxinus mandshurica 31.91 15.10 257.89 21.16 1 488.58 600 3 100 443.16 24.95 15.20 34.80 3.97
    下载: 导出CSV 
    | 显示表格

    以往大多选定某种非线性方程描述树种的树高−胸径关系,而混交林相对于纯林或年龄结构和林层单一的林分来说,其树高与胸径的关系更为复杂[13]。根据以往的研究,基于树高−胸径散点图和树种的分布特征,选择以下11个常用的树高模型作为候选基础模型(表4),对树高−胸径分布数据进行非线性拟合,比较各模型的拟合精度,筛选预测精度最高的作为研究区树高−胸径的最优基础模型。

    表  4  候选树高−胸径关系模型
    Table  4.  Candidate models of tree height-DBH relation
    模型编号 Model No. 表达式 Equation 参考文献 Reference
    M1 H=1.3+β1Dβ2 [14]
    M2 H=1.3+β1[D1+D]β2 [15]
    M3 H=1.3+exp(β1+β2D+1) [16]
    M4 H=1.3+β1[1exp(β2D)] [17]
    M5 H=1.3+β1exp(β2D) [18]
    M6 H=1.3+β1[1exp(β2D)]β3 [19]
    M7 H=1.3+β1D[β2D(β3)] [20]
    M8 H=1.3+exp(β1+β2D+β3) [21]
    M9 H=1.3+β1exp[β2exp(β3D)] [22]
    M10 H=1.3+D2β1+β2D+β3D2 [23]
    M11 H=1.3+β11+β2exp(β3D) [24]
    注:H为树高预测值;D为胸径;β1β2β3为模型参数。Notes: H is the predicted tree height; D is the DBH; β1, β2 and β3 are the model parameters.
    下载: 导出CSV 
    | 显示表格

    但在结构更加复杂的针阔混交林中,仅仅将胸径作为预测变量,无法反映其他林分变量对树高的影响。本研究采用再参数化的方法建立树高−胸径模型,依据林分变量与树高之间的关系,适当融合林分密度、胸高断面积和优势木高度等变量。具体做法如下:基于所选择的最优基础模型,对林分变量与树高进行相关性分析;比较不同变量引入对模型的影响,筛选模型中的其他变量,确定与树高相关性较大且能显著提高模型预测能力的变量作为协变量;基于最优基础模型构建各树种预测精度最高的广义模型。

    分位数回归是利用解释变量估计响应变量的条件分位数分布,全面刻画整个条件分布的综合回归方法[2526]。不仅可以度量响应变量对分布中心的影响,而且可以度量对分布上尾和下尾的影响。最小二乘法假定误差服从正态分布,而分位数回归的响应变量的平均预测值对残差分布没有特殊要求[27],从而确保参数估计的稳健性和模型拟合的灵活性。最小二乘回归的目标是最小化误差平方和,分位数回归的核心思想是从均值推广到分位数,也是通过最小化残差绝对值的非对称损失函数来估计分为点处参数值。

    min (1)

    式中:τ表示要估计的分位数, {y}_{\tau } {\hat{y}}_{\tau } 分别为响应变量在分位数τ处的观测值和估计值。

    分位数组合法利用抽样数据拟合不同组合的分位数曲线,基于相邻两条分位数回归曲线插值得到修正的分位数曲线。假设第i个样方第j个树木的胸径值为 {x}_{ij} ,树高观测值为 {y}_{ij} 分布在第m和(m + 1)分位数估计值之间,即 {\hat{y}}_{m}\left({x}_{ij}\right) {y}_{ij} {\hat{y}}_{m + 1}\left({x}_{ij}\right) 。使用插值法计算通过该点的修正曲线。

    {\hat{y}}_{ij}=\alpha {\hat{y}}_{m}\left({x}_{ij}\right) + (1-\alpha ){\hat{y}}_{m + 1}\left({x}_{ij}\right) (2)

    式中: \alpha = \dfrac{{\hat{y}}_{m + 1}\left({x}_{ij}\right)-{y}_{ij}}{{\hat{y}}_{m + 1}\left({x}_{ij}\right)-{\hat{y}}_{m}\left({x}_{ij}\right)} 表示插值系数, {\hat{y}}_{ij} 为树高估计值, {\hat{y}}_{m}\left({x}_{ij}\right)\mathrm{和}{\hat{y}}_{m + 1}\left({x}_{ij}\right) 分别为胸径为 {x}_{ij} 的树木在第m和(m + 1)分位数处的树高估计值。

    当抽样数量大于1时,在各分位数组合中选择使平均误差值符号发生改变的相邻曲线。当胸径值大于最大分位数回归的预测值时,将 \hat{y} m \hat{y} m + 1分别定义为 \hat{y} m − 1 \hat{y} m,之后利用公式(2)得到树高曲线。

    由于基础模型只有一个预测变量缺乏其他变量信息,很难完整描述响应变量在不同分位点的分布。本研究选择9个分位点(τ = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9),基于所选的最优基础模型拟合树高−胸径曲线,同时使用三分位数、五分位数和九分位数组合法得到树高修正曲线,通过模型评价和比较得到基于最佳分位数的基础模型,进而在该模型中融入林分变量构建广义分位数模型。

    在筛选出的最优基础模型和广义模型加入随机效应构建非线性混合模型,本研究中选择样方为随机效应,混合效应模型一般形式可表示为

    \hat{\boldsymbol{H}}_{ij}=f\left(\hat{\boldsymbol{\phi}}_i,\hat{\boldsymbol{D}}_{ij}\right)+\boldsymbol{\varepsilon}_{ij},\ {\boldsymbol{\varepsilon}}_{ij}\sim{{{{\rm{N}}}}}(0,\ \sigma^2) (3)

    式中: \hat{\boldsymbol{H}}_{ij} \hat{\boldsymbol{D}}_{ij} 分别为第i个样方第j株树的树高和胸径预测向量,f为包含树高和胸径变量的函数表达式, \hat{\boldsymbol{\phi}}_i 为混合效应参数向量,εij为误差项。

    参数向量 \hat{\boldsymbol{\phi}}_i 的固定效应参数和随机效应参数可表示为

    {\boldsymbol{\phi}}_i=\hat{\boldsymbol{X}}_i \hat{\boldsymbol{\beta}}+ {\hat{\boldsymbol{Z}}_i} {\hat{\boldsymbol{u}}_i},\; {\hat{\boldsymbol{u}}_i} \sim {\text{N}}(0,\; {\text{φ}}) (4)

    式中: \hat{\boldsymbol{\beta }} 为固定效应参数向量, {\hat{\boldsymbol{u}}}_{i} 为第i个样方的随机效应参数向量,遵循均值为0、协方差矩阵为{\text{φ}} 的正态分布, {\hat{\boldsymbol{X}}}_{i} {\hat{\boldsymbol{{Z}}}}_{i} 分别是 \hat{\boldsymbol{\beta }} {\hat{\boldsymbol{u}}}_{i} 的设计矩阵。

    非线性混合效应模型中,方差−协方差结构由包括相关因子和加权因子的矩阵Ri来表示。

    {\boldsymbol{R}}_{i}={\sigma }^{2}{\boldsymbol{G}}_{i}^{0.5}{\boldsymbol{\varGamma }}_{i}{\boldsymbol{G}}_{i}^{0.5} (5)

    式中:σ2指模型的残差方差值,Γi为组内误差相关性结构,Gi为描述样方内方差异质性的对角矩阵,其对角元素由方差函数提供[28]

    随机效应的方差−协方差结构(D)反映了样方间的差异,矩阵因随机参数个数的不同也会有不同的形式,本研究中选用广义正定矩阵结构,当混合模型中具有两个随机参数时表达式如下

    \boldsymbol{D}=\left[\begin{array}{cc}{\sigma }_{{u}_{1}}^{2}& {\sigma }_{{u}_{1}{u}_{2}}\\ {\sigma }_{{u}_{1}{u}_{2}}& {\sigma }_{{u}_{2}}^{2}\end{array}\right] (6)

    式中: {u}_{1} {u}_{2} 为随机参数, {\sigma }_{{u}_{1}}^{2} {\sigma }_{{u}_{2}}^{2} 分别为 {u}_{1} {u}_{2} 的方差, {\sigma }_{{u}_{1}{u}_{2}} {u}_{1} {u}_{2} 的协方差。

    混合效应模型中固定效应部分的检验方法与传统模型相同,对于随机效应的检验则需要计算随机参数值,本研究需计算在拟合数据中不涉及的样方的随机参数,利用式7[29]进行计算。

    \hat{u}_i\approx\hat{D}\hat{Z}_i^{\mathrm{T}}\left(\hat{Z}_i\hat{D}\hat{Z}_i^{\mathrm{T}}+\hat{R}_i\right)^{-1}\hat{e}_i (7)

    式中: \hat{\boldsymbol{D}} 为随机效应参数的方差协方差矩阵, {\hat{\boldsymbol{R}}}_{i} 为样方内方差协方差矩阵, {\hat{\boldsymbol{Z}}}_{i} 为设计矩阵, {\hat{\boldsymbol{Z}}}_{i}^{\mathrm{T}} 为转置矩阵, {\hat{e}}_{i} 为实际值减去用固定效应参数计算的预测值。

    研究中通过决定系数(R2)、平均绝对误差(MAE)和均方根误差(RMSE)评价和比较模型拟合结果,并利用未参与建模的20%独立样本检验11个候选非线性模型,上述评价指标的计算公式。

    {R}^{2}=1-{\sum} _{i,j=1}^{n}{\left({H}_{ij}-{\hat{H}}_{ij}\right)}^{2}/{\sum} _{i,j=1}^{n}{\left({H}_{ij}-{\overline{H}}_{ij}\right)}^{2} (8)
    \mathrm{M}\mathrm{A}\mathrm{E}={\sum} _{i,j=1}^{n}\left|{H}_{ij}-{\hat{H}}_{ij}\right|/n (9)
    \mathrm{R}\mathrm{M}\mathrm{S}\mathrm{E}=\sqrt{{\sum} _{i,j=1}^{n}{\left({H}_{ij}-{\hat{H}}_{ij}\right)}^{2}/n} (10)

    式中: {H}_{ij} {\hat{H}}_{ij} 为第i个样方第j个样本的树高实际值和预测值, {\overline{H}}_{ij} 为预测值的均值;n 为样本数。

    除考虑以上指标外,还选择赤池信息准则(AIC)、贝叶斯信息准则(BIC)和似然比检验(LRT)对混合效应模型进行评价,AIC和BIC数值越小则表明模型拟合精度越高,LRT得到的P值若 < 0.000 1,则认为两个模型差异显著。

    \mathrm{AIC}=-2\mathrm{\mathit{L}}+2\lambda (11)
    \mathrm{BIC}=-2\mathrm{\mathit{L}}+\lambda\mathrm{ln}\ n (12)

    式中:L为模型的对数似然值,λ 为固定效应参数的个数。

    基于红松、色木槭、紫椴和水曲柳的建模数据,对表4中的11个基础模型进行拟合,比较其MAE、RMSE和R2值(表5)。结果表明,对于4个树种的树高−胸径关系拟合,各模型均具有良好表现。MAE范围为1.887 ~ 2.818,RMSE范围为2.557 ~ 3.582,R2则为0.552 ~ 0.811,其中红松、紫椴水曲柳以M6即Richards模型拟合结果最优,色木槭的Richards模型拟合结果与其最优模型相差极小。使用检验数据对拟合模型进行评价比较,得到相同结果。且Richards方程生物学意义明确、灵活性较高易收敛,被认为很适合描述树高与胸径之间的关系[3031],因此选择Richards作为4个树种的基础模型。

    表  5  候选模型拟合与评价
    Table  5.  Fitting and evaluation of candidate models
    模型
    Model
    红松 Pinus koraiensis 色木槭 Acer mono 紫椴 Tilia amurensis 水曲柳 Fraxinus mandshurica
    MAE RMSE R2 MAE RMSE R2 MAE RMSE R2 MAE RMSE R2
    M1 2.381 3.153 0.780 2.118 2.846 0.781 2.818 3.582 0.684 2.248 2.774 0.552
    M2 2.361 3.128 0.784 2.080 2.806 0.787 2.521 3.357 0.722 2.079 2.571 0.615
    M3 2.211 2.982 0.803 1.926 2.697 0.803 2.503 3.344 0.724 2.083 2.575 0.614
    M4 2.179 2.945 0.808 1.895 2.654 0.809 2.530 3.343 0.724 2.088 2.575 0.614
    M5 2.279 3.032 0.797 2.030 2.772 0.793 2.547 3.374 0.719 2.076 2.568 0.616
    M6 2.140 2.925 0.811 1.890 2.654 0.809 2.491 3.332 0.726 2.063 2.557 0.619
    M7 2.149 2.936 0.809 1.887 2.653 0.809 2.495 3.337 0.725 2.069 2.562 0.617
    M8 2.144 2.932 0.810 1.888 2.652 0.810 2.496 3.335 0.725 2.074 2.565 0.617
    M9 2.174 2.945 0.808 1.931 2.691 0.804 2.513 3.348 0.723 2.077 2.570 0.615
    M10 2.146 2.935 0.809 1.887 2.652 0.810 2.493 3.336 0.725 2.068 2.560 0.618
    M11 2.231 2.991 0.802 1.986 2.742 0.796 2.553 3.381 0.718 2.091 2.587 0.610
    注:MAE.平均绝对误差;RMSE.均方根误差;R2.决定系数。下同。Notes: MAE, mean absolute error; RMSE, root mean squared error; R2, coefficient of determination. The same below.
    下载: 导出CSV 
    | 显示表格

    表3的林分特征因子进行相关性分析,结果如表6所示。4个树种的树高和胸径与优势木高呈显著正相关关系(P < 0.01),色木槭和紫椴树高、胸径与林分密度呈显著负相关关系(P < 0.01)。因此将优势木高和林分密度两个变量分别引入基础模型构建广义模型,但包含林分密度因子的色木槭和紫椴树高方程不收敛,而引入胸高断面积的水曲柳树高方程也不收敛,这表明添加优势木高度可显著提高模型精度。因此本研究最优广义模型的表达式为

    表  6  树高与林分变量相关性分析
    Table  6.  Correlation analysis between tree height and stand variables
    树种 Tree species因子 Factor树高
    Tree height
    胸径
    DBH
    胸高断面积
    Basal area
    林分密度
    Stand density
    优势木高
    Dominate tree height
    红松 Pinus koraiensis树高 Tree height1.000
    胸径 DBH0.855**1.000
    胸高断面积 Basal area0.0200.0401.000
    林分密度 Stand density0.0100.0060.092**1.000
    优势木高 Dominate tree height0.587**0.440**0.069*0.0181.000
    色木槭 Acer mono树高Tree height1.000
    胸径 DBH0.845**1.000
    胸高断面积 Basal area0.0310.0301.000
    林分密度 Stand density−0.284**−0.270**0.089**1.000
    优势木高 Dominate tree height0.484**0.353**0.139**−0.261**1.000
    紫椴 Tilia amurensis树高 Tree height1.000
    胸径 DBH0.776**1.000
    胸高断面积 Basal area0.068**0.064**1.000
    林分密度 Stand density−0.207**−0.200**0.105**1.000
    优势木高 Dominate tree height0.537**0.294**0.178**−0.175**1.000
    水曲柳 Fraxinus mandshurica树高Tree height1.000
    胸径 DBH0.684**1.000
    胸高断面积 Basal area0.067**0.133**1.000
    林分密度 Stand density−0.037−0.081**0.084**1.000
    优势木高 Dominate tree height0.256**0.264**−0.021**0.318**1.000
    注:** 表示两个变量在 P < 0.01 水平上显著相关,* 表示两个变量在 P < 0.05 水平上显著相关。Notes: ** means significant correlation at P < 0.01 level between two variables, * means significant correlation at P < 0.05 level between two variables.
    下载: 导出CSV 
    | 显示表格
    H=1.3+(\beta_1+a_1H_{\mathrm{D},i})\left[1-\exp\left(-\beta_2D\right)\right]^{\beta_3} (13)

    式中:HD为某个树种的优势木高,HD,i为第i个样方的优势木高, {\beta }_{1} {\beta }_{2} {\beta }_{3} {a}_{1} 为模型参数。

    在不同分位点处(τ = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9),基于Richards方程利用分位数回归分别拟合树高−胸径关系模型(表5)。利用回归评价指标检验不同分位数模型对观测值的拟合程度(表7),基础模型及分位数回归模型的决定系数R2值均在0.6以上,表明各分位点的树高方程均具有较好的拟合效果。其中τ = 0.5的中位数回归与非线性回归的拟合效果基本一致,反映了分位数回归模型的灵活性。

    表  7  分位数回归模型拟合精度
    Table  7.  Goodness-of-fit tests for quantile regression models
    分位数 Quantile红松 Pinus koraiensis色木槭 Acer mono紫椴 Tilia amurensis水曲柳 Fraxinus mandshurica
    MAERMSER2MAERMSER2MAERMSER2MAERMSER2
    Richards2.1402.9250.8111.8902.6540.8092.4913.3320.7262.0632.5570.619
    0.13.4154.4970.7803.1034.1050.7704.0295.1660.6963.4264.1560.614
    0.22.7793.8260.7932.4603.4060.7933.1404.2380.7152.6803.3870.615
    0.32.3683.3220.8052.0812.9630.8042.7423.7910.7202.3232.9760.615
    0.42.1793.0630.8091.9242.7640.8072.5203.4930.7232.1162.6930.615
    0.52.1222.9420.8111.8732.6660.8092.4613.3740.7242.0562.5830.617
    0.62.1912.9530.8101.9182.6750.8082.5423.3630.7242.1192.5960.616
    0.72.4403.1810.8072.1612.8970.8032.8303.6300.7192.3722.8460.614
    0.82.9123.6900.7972.5773.3540.7913.3884.2680.6972.8873.4220.615
    0.94.0815.0270.7533.4514.2930.7624.6295.5190.6793.9654.5820.602
    下载: 导出CSV 
    | 显示表格

    利用未参与建模的20%测试数据检验模型的预测能力(表8)。分位数组合法分别采用三分位数组合(τ = 0.1,0.5,0.9)、五分位数组合(τ = 0.1,0.3,0.5,0.7,0.9)和九分位数组合(τ = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9)。随机抽样或抽取平均木时幼树会带来较大预测误差,基于王君杰等[32]的研究结果,选取10株最大树木计算插值系数进行预测。结果表明:对于分位数组合法,三分位数和五分位数组合的预测精度高于九分位数组合;其中红松树高的拟合效果略优于中位数回归,但其他树种模型预测精度均低于中位数回归。因此,本研究选择中位数为4个树种的最优分位数树高曲线式(13)构建广义分位数模型。

    表  8  分位数模型验证评价结果
    Table  8.  Quantile model validation evaluation results
    方法 Method红松 Pinus koraiensis色木槭 Acer mono紫椴 Tilia amurensis水曲柳 Fraxinus mandshurica
    RMSEMAERMSEMAERMSEMAERMSEMAE
    中位数回归 Median regression2.9632.0942.8982.0203.5402.5492.5602.050
    九分位数组合 Nine quantile combination2.9562.1182.9972.1853.5102.7612.6012.150
    五分位数组合 Five quantile combination2.9272.0952.9712.1553.5072.7582.6022.152
    三分位数组合 Triquartile combination2.9262.0922.9832.1623.5962.8512.5922.139
    下载: 导出CSV 
    | 显示表格
    表  11  树高−胸径模型拟合评价统计量
    Table  11.  Fitting and evaluation statistics of tree height-DBH models
    树种
    Tree species
    指标
    Index
    基础模型
    Base model
    广义模型
    Generalized
    model
    中位数回归 Median quantile regression 非线性混合效应 Nonlinear mixed-effects
    简单模型
    Simple model
    广义模型
    Generalized model
    简单模型
    Simple model
    广义模型
    Generalized model
    红松 Pinus koraiensis MAE 2.140 1.669 2.122 1.666 1.259 1.358
    RMSE 2.925 2.225 2.942 2.228 1.687 1.826
    R2 0.811 0.891 0.811 0.891 0.937 0.927
    色木槭 Acer mono MAE 1.890 1.575 1.873 1.568 1.258 1.359
    RMSE 2.654 2.159 2.666 2.166 1.740 1.875
    R2 0.809 0.875 0.809 0.874 0.919 0.906
    紫椴 Tilia amurensis MAE 2.491 1.722 2.461 1.706 1.399 1.611
    RMSE 3.332 3.379 3.374 2.401 1.955 2.235
    R2 0.726 0.860 0.724 0.858 0.906 0.877
    水曲柳 Fraxinus mandshurica MAE 2.063 2.041 2.056 2.031 2.004 2.527
    RMSE 2.557 2.524 2.583 2.536 2.478 2.038
    R2 0.619 0.629 0.617 0.627 0.643 0.629
    下载: 导出CSV 
    | 显示表格

    在基础模型和广义模型基础上考虑样方效应对树高的影响,在4个树种的基础模型上增加样方为随机效应,共8种情况收敛(表9),红松、紫椴和水曲柳增加1个随机效应参数时拟合效果最好,而色木槭引入2个随机参数的混合模型预测精度更高,模型表达式分别为

    表  9  混合效应模型验证评价结果
    Table  9.  Verification and evaluation results of mixed effect model
    树种
    Tree species
    基础混合效应模型
    Base mixed-effect model
    广义混合效应模型
    Generalized mixed-effect model
    随机参数
    Random
    parameter
    AIC BIC L LRT P 随机参数
    Random
    parameter
    AIC BIC L LRT P
    红松
    Pinus koraiensis
    {u}_{1} 5 619.3 5 644.8 −2 804.6 {u}_{1} 5 346.9 5 377.5 −2 667.4
    {u}_{1} , {u}_{2} 5 643.2 5 678.9 −2 814.6 19.878 2 < 0.001 {u}_{3} 5 354.2 5 384.8 −2 671.1
    {u}_{1} , {u}_{3} 5 350.9 5 391.8 −2 667.5 7.238 0 0.026 8
    色木槭
    Acer mono
    {u}_{2} 10 966.2 10 995.2 −5 478.1 {u}_{1} 10 597.5 10 632.3 −5 292.7
    {u}_{1} , {u}_{2} 10 962.2 11 002.9 −5 474.1 7.964 9 0.018 6 {u}_{3} 10 602.8 10 637.6 −5 295.4
    {u}_{1} , {u}_{3} 10 594.2 10 640.6 −5 289.1 12.571 9 0.001 9
    {u}_{1} , {u}_{2} 10 601.4 10 647.8 −5 292.8
    紫椴
    Tilia amurensis
    {u}_{1} 7 109.6 7 136.0 −3 549.8 {u}_{1} 6 662.5 6 694.2 −3 325.2
    {u}_{3} 6 661.8 6 693.5 −3 324.9
    {u}_{1},{u}_{3} 6 665.8 6 708.0 −3 324.9 4.415 0 0.998 7
    水曲柳
    Fraxinus
    mandshurica
    {u}_{1} 7 007.8 7 034.4 −3 498.9 {u}_{1} 6 966.0 6 997.8 −3 477.0
    {u}_{2} 7 009.0 7 035.5 −3 499.5 {u}_{3} 6 964.3 6 996.1 −3 476.1
    {u}_{1} , {u}_{2} 7 011.8 7 049.0 −3 498.9 1.126 9 0.569 2 {u}_{1},{u}_{3} 6 964.9 7 007.4 −3 474.5 3.322 9 0.189 9
    注:AIC.赤池信息准则;BIC.贝叶斯信息准则;L.对数似然值;LRT.似然比检验。Notes: AIC, Akaike information criterion; BIC, Bayesian information criterion; L, log-likelihood value; LRT, likelihood ratio test.
    下载: 导出CSV 
    | 显示表格
    H=1.3 + {({\beta }_{1} + {u}_{1i})[1-\exp(-{\beta }_{2}D\left)\right]}^{{\beta }_{3}} (14)
    H=1.3 + {({\beta }_{1} + {u}_{1i})[1-\exp(-{\beta }_{2}D\left)\right]}^{{(\beta }_{3 + }{u}_{2i})} (15)

    式中: {u}_{1i} {u}_{2i} 为样方水平的随机参数。

    在广义模型基础上,引入不同组合形式的随机效应作为自变量构建广义混合效应模型,4个树种共有13种情况收敛(表9),当引入3个以上随机参数时,模型不收敛;仅引入参数 {u}_{1} 的红松树高模型和仅引入参数 {u}_{3} 的紫椴及水曲柳树高模型AIC、BIC值最小,且L值达到最大,而色木槭以随机效应作用于 {u}_{1} {u}_{3} 两个参数时模型预测精度最高,一般表达式分别为

    H=1.3+(\beta_1+u_{1i}+a_1H_{\mathrm{D},\mathit{\mathrm{\mathit{i}}}})\left[1-\exp\left(-\beta_2D\right)\right]^{\beta_3} (16)
    H=1.3+\left[\beta_1+u_{1i}+(a_1+u_{3i})H_{\mathrm{D},\mathrm{\mathit{i}}}\right]\left[1-\exp\left(-\beta_2D\right)\right]^{\beta_3} (17)
    H=1.3+\left[\beta_1+(a_1+u_{3i})H_{\mathrm{D},\rm{\mathit{i}}}\right]\left[1-\exp\left(-\beta_2D\right)\right]^{\beta_3} (18)

    式中: {u}_{1i} {u}_{3i} 为样方水平的随机参数。

    选择各树种拟合效果最好的混合效应模型,参数估计值见表10

    表  10  混合效应模型参数估计
    Table  10.  Parameter estimates of mixed-effect models
    参数
    Parameter
    基础混合效应模型
    Base mixed-effect model
    广义混合效应模型
    Generalized mixed-effect model
    红松
    Pinus
    koraiensis
    色木槭
    Acer mono
    紫椴
    Tilia amurensis
    水曲柳
    Fraxinus
    mandshurica
    红松
    Pinus koraiensis
    色木槭
    Acer mono
    紫椴
    Tilia amurensis
    水曲柳
    Fraxinus
    mandshurica
    公式编号
    Equation No.
    (14) (15) (14) (14) (16) (17) (18) (18)
    a 21.192(0.387*** 18.626(0.251*** 19.366(0.263*** 20.021(0.190*** 9.186(0.535*** 7.776(0.467*** 4.522(0.468*** 16.513(0.523***
    b 0.065(0.004*** 0.085(0.004*** 0.100(0.005*** 0.101(0.007*** 0.071(0.004*** 0.092(0.004* 0.104(0.005*** 0.110(0.007***
    c 1.366(0.056*** 1.142(0.033*** 1.367(0.066*** 1.445(0.126*** 1.305(0.057*** 1.142(0.034*** 1.322(0.066*** 1.567(0.142***
    a1 0.577(0.024*** 0.554(0.023*** 0.715(0.022*** 0.132(0.019***
    {\mathrm{\sigma }}_{{u}_{1}}^{2} 14.927 2 10.920 5 11.410 7 0.201 9 2.877 0 14.850 5
    {\mathrm{\sigma }}_{{u}_{2}}^{2} 0.011 9
    {{\mathrm{\sigma }}^{2}}_{{u}_{3}} 0.038 0 0.001 7 0.000 4
    {\mathrm{\sigma }}_{{u}_{1}{u}_{2}} 0.162 1
    {\mathrm{\sigma }}_{{u}_{1}{u}_{3}} −0.718 2
    {\mathrm{\sigma }}^{2} 3.637 0 3.638 8 4.908 0 6.422 3 3.811 3 3.858 9 5.275 9 6.201 5
    注:***表示在 P < 0.001 水平上显著。Note: *** means significant at P < 0.001 level.
    下载: 导出CSV 
    | 显示表格

    4个树种的各模型拟合统计量见表11。结果表明:基于混合效应的模型精度明显高于不含随机参数的模型;在不含随机效应的模型中,中位数回归形式和一般形式的广义模型精度差异较小;广义中位数模型在预测精度上高于基础模型的中位数回归,红松、色木槭、紫椴和水曲柳的模型R2分别提高了9.9%、8.0%、18.5%和1.6%,MAE值分别降低了27.4%、19.5%、44.3%和1.2%,RMSE值分别下降了32.0%、23.1%、40.5%和1.9%。在增加随机效应的模型中,基础混合模型的拟合结果略优于广义混合模型,4个树种R2仅分别提高了1.1%、1.4%、3.3%和2.2%,但两者并未有显著差异。红松、色木槭、紫椴和水曲柳的基础混合效应和广义混合效应模型的残差分布图如图1所示,残差分布较为均匀,没有出现明显的不规则形状,这说明各树种的混合模型均没有明显异方差现象。

    图  1  基础混合效应模型(A1、B1、C1、D1)和广义混合效应模型(A2、B2、C2、D2)的残差分布图
    Figure  1.  Residual plots of base mixed-effect models (A1, B1, C1, D1) and generalized mixed-effect models (A2, B2, C2, D2)

    混合效应模型中随机效应部分的检验需要估算随机参数值,本研究以20%的样木作为检验数据,利用公式(7)分别计算未包含在拟合数据的样方的随机参数;同时与基础模型、广义模型和分位数模型的预测精度进行比较,检验结果如表12。可以看出,混合效应模型的MAE和RMSE值均小于其他模型,拟合效果相较分位数回归得到了显著提高,其中基础混合模型相比于广义混合模型预测精度更高。

    表  12  树高−胸径模型统计检验
    Table  12.  Goodness-of-fit statistics of tree height-DBH models
    树种
    Tree species
    指标
    Index
    基础模型
    Base model
    广义模型
    Generalized
    model
    中位数回归 Median quantile regression非线性混合效应 Nonlinear mixed-effects
    简单模型
    Simple model
    广义模型
    Generalized model
    简单模型
    Simple model
    广义模型
    Generalized model
    红松 Pinus koraiensisMAE2.1282.1192.0942.1631.6001.856
    RMSE2.9633.0242.9633.0632.2162.635
    色木槭 Acer monoMAE2.0181.9062.0202.0521.6201.712
    RMSE2.8672.7822.8982.9442.3412.529
    紫椴 Tilia amurensisMAE2.5711.6852.5492.6971.6091.635
    RMSE3.4712.3493.5403.7542.0582.206
    水曲柳
    Fraxinus mandshurica
    MAE2.0412.0372.0502.0512.0392.049
    RMSE2.5262.5172.5602.5612.5282.534
    下载: 导出CSV 
    | 显示表格

    分位数回归既能描述均值条件下树高的分布趋势,也能描述特定分位点树高−胸径曲线的形状变化[33]。传统的最小二乘法易受残差分布限制,与之相比分位数回归具有更加稳健的参数估计值,模型稳定性和灵活性也更强,能够避免选取数据的人为影响。在不同分位点拟合树高−胸径模型,同时基于参数估计和建模样本计算三分位数、五分位数和九分位数组合的插值系数,王君杰等[34]最近发现中位数回归优于非线性回归;本文与之不同发现中位数回归与非线性回归的拟合效果基本一致。经过对11个非线性模型的评价比较,选择Richards方程作为基础模型进行分位数回归,其中中位数回归拟合效果最好,4个树种的R2值均大于0.6,对树高的预测能力不逊于非线性回归,这表明Richards方程能广泛适用树高−胸径模型的拟合。

    三分位数组合的预测效果均优于五分位数和九分位数组合,尤其在红松树高的预测中,两类分位数组合的预测精度相较中位数和非线性回归均有提高;以往研究显示五分位数组合对树高曲线和削度方程的拟合要优于三分位数组合,但都没有考虑抽取方案对分位数组合的影响[35]。因此抽取方案的不同会导致分位数组合之间的差异不同,总体上三分位数组合(τ = 0.1,0.5,0.9)的预测能力较好且效果稳定,本研究中选取多个分位数组合也并未使预测结果有明显提高[32]

    在林分生长过程中,有学者利用含有优势木胸径、树高和林分密度等因子的模型来描述树高−胸径关系,模型拟合结果得到提升[36]。为了比较不同样方间树高−胸径的变化规律,本研究对胸高断面积、林分密度和优势木高3个林分变量进行Pearson相关性分析,发现引入变量的广义模型的拟合效果均优于基础模型,当引入1个变量时,红松、色木槭、紫椴和水曲柳均仅增加优势木高变量的模型可以收敛,部分其他林分变量与树高有显著关系,但含优势木高的模型相较含林分密度或胸高断面积的广义模型对树高−胸径的关系影响更大;4个树种中水曲柳的广义模型相较基础模型的拟合效果提升幅度最小,R2仅提高0.01,这可能与水曲柳的直径分布有关,幼树数量较少,成熟林木较多,在林分中多位于林冠层,内部资源和空间竞争相对稳定,处于优势地位,这类大个体的树高不易受林分变量较大影响。这说明即使是在具有相似地形和气候条件的林分中,不同树种的树高−胸径关系可能不同。森林群落中特定的环境条件和树木之间的复杂联系可能会导致树高−胸径生长关系不同[37]。优势木高反映了立地质量以及林分生长和收获,优势木高的参数估计值为正值,这表明由于立地质量的不同,优势木高与树高之间呈正相关关系,引入优势木高的广义模型更适用于树高−胸径模拟的拟合。

    基于基础模型和各树种的最优广义模型,对应构建以样方水平为随机效应的混合模型,比较不同随机参数下各树种树高模型的拟合效果。研究发现,红松、紫椴和水曲柳基础混合效应模型均以增加随机参数 {u}_{1} ,广义混合效应模型中红松增加参数 {u}_{1} ,紫椴和水曲柳将 {u}_{3} 同时作为随机参数,拟合效果最优;色木槭基础混合效应模型增加随机参数 {u}_{1} {u}_{2} ,广义混合效应模型则是将 {u}_{1} {u}_{3} 同时作为随机参数,预测精度最高。通过对各模型进行检验评价,结果表明混合效应模型对树高的预测效果优于固定效应模型,这与王冬至等[38]的研究结果一致,表明考虑样方水平的随机影响能更好地描述树高生长规律。此外,本研究的模型评价和检验结果显示基础混合效应模型的预测能力略优于广义混合效应模型,但未有显著差异,同样的情况也出现在陈浩等[39]对马尾松和周晏平等[40]对樟子松的研究结果中,说明在增加样方效应后,林分变量对树高−胸径关系的影响被样方水平所涵盖,对混合效应模型的影响较弱。在Zang等[41]对中国北部和东北落叶松人工林树高−胸径关系的研究中,基础混合效应模型比其他模型提供了更好的拟合和更稳定的估计,表明不需要额外的树木和林分变量。

    本研究基于分位数回归和混合效应模型对阔叶红松林中优势树种红松、色木槭、紫椴和水曲柳的树高−胸径关系进行模型拟合,在11个非线性候选模型中确定Richards方程为基础模型,并引入优势木高构建广义模型,模型拟合结果得到较大提高;4个树种基于中位数回归的基础模型和广义模型拟合精度与基础模型的非线性回归相近。

    加入样方水平为随机效应后,混合效应模型相较固定效应模型更具优越性,模型预测能力得到显著提高,而基础混合模型和广义混合模型的拟合优度差异较小。这表明林分变量和随机效应对树高模型的影响是显著的,在不考虑样方效应的情况下,可选择广义模型作为东北地区针阔混交林树高模型;若考虑样方效应,基础混合效应模型能提供更稳定的预测趋势和拟合结果。

  • 图  1   昌黎县2024年3—5月份日平均温度变化及≥10 ℃有效积温

    Figure  1.   Daily average temperature changes and effective accumulated temperature (≥ 10 ℃) in Changli County from March to May in 2024

    图  2   核桃萌芽及雌雄花开放进程

    A. 休眠期;B. 萌动期;C. 膨大期;D. 展叶期;E. 叶片生长期;F. 雄花膨大期;G. 雄花伸长期;H. 雄花散粉初期;I. 雄花散粉盛期;J. 雄花散粉末期;K. 雌花花蕾期;L、M. 雌花初期;N. 雌花盛期;O. 雌花末期

    Figure  2.   Walnut germination and male and female flower blooming process

    图  3   84份核桃种质资源2024年萌芽期统计

    Figure  3.   Statistics of germination period of 84 walnut germplasm resources in 2024

    图  4   84份核桃种质资源2024年雌雄花盛期统计

    Figure  4.   Statistics of blooming period of 84 walnut germplasm resources in 2024

    表  2   雌雄花期起始所需平均≥10 ℃有效积温及种质资源数量统计

    Table  2   Statistics of average effective accumulated temperature of ≥ 10 °C required for the starting of flowering period and the number of germplasm resources

    日期 ≥10 ℃有效
    积温/(℃·d)
    当日雌花盛期起始
    种质资源数量
    当日雄花盛期起始
    种质资源数量
    划分时期
    04−17 67.0 0 2 早花期
    04−18 75.5 13 18
    04−19 80.0 4 8
    04−20 83.5 1 4
    04−21 88.0 4 3
    04−22 92.5 9 10
    04−23 92.5 1 2
    04−24 99.0 3 0
    04−25 107.5 10 5
    04−26 113.0 10 2
    04−27 119.0 8 6 中花期
    04−28 122.5 3 1
    04−29 123.0 3 4
    04−30 130.0 8 11
    05−01 139.5 1 1
    05−02 149.0 2 1
    05−03 158.5 2 3
    05−04 164.0 0 1
    05−09 204.5 0 1 晚花期
    05−10 214.5 0 1
    05−15 265.0 2 0
    注:经计算,雌花盛期起始(04−18—05−15)所需平均 ≥ 10 ℃有效积温为109.25 ℃·d,雄花盛期起始(04−17—05−10)所需平均 ≥ 10 ℃有效积温为103.69 ℃·d。
    下载: 导出CSV

    表  1   萌芽期起始所需平均≥10 ℃有效积温和萌芽资源统计

    Table  1   Average effective accumulated temperature of ≥ 10 °C required for the starting of germination period and statistics of germination resources

    日期 ≥10 ℃有效积温/(℃·d) 当日萌芽期起始种质资源
    03−24 4.5 ‘绿岭’‘中核1号’‘鲁果5号’‘实生4’‘辽宁1号’‘核桃楸’‘辽瑞丰’
    03−25 4.5 ‘绿波’‘金薄丰’‘礼品1号’‘早硕’‘实生6’‘实生10’
    03−26 4.5 ‘实生1’‘纸皮3号’‘西岭’‘丰香’‘礼品2号’‘辽宁4号’‘魁香’‘扎木台5号’‘安边2号’‘新新2’‘鲁果6号’‘晋丰’‘辽宁6号’‘岱香’‘实生12’‘实生15’‘清香’‘硕香’‘麻核桃’
    03−27 7.0 ‘新光’‘元宝’‘西洛9号’‘龙珠’‘硕宝’‘鲁果7号’‘A53’‘西林5号’‘Pedro’‘美香’‘寒丰’‘陇南11号’‘极早丰’‘实生7’‘庆阳11号’‘B76’‘辽宁5号’‘扎木台3号’‘磨盘’‘实生9’‘实生14’‘2-34’‘实生16’
    03−28 7.0 ‘扎343’‘实生2’‘辽宁7号’‘新早丰’‘红光’‘华山5号’‘青林’‘实生5’‘黑核桃’‘A19’‘10901’‘辽宁10号’‘实生8’‘纸皮扎01’‘旮旯铺’‘实生11’‘实生13’
    03−29 7.5 ‘中核3号’‘西洛7号’‘金薄香2号’‘金薄香4号’‘西林1号’‘小果黑核桃’‘绿香’‘晋香’‘花生1号’‘4-22’
    03−30 7.5 ‘中核2号’
    03−31 7.5 ‘实生3’
    注:经计算,萌芽期起始(03−24—03−31)所需平均≥10 ℃有效积温为6.12 ℃·d。
    下载: 导出CSV

    表  3   84份核桃种质资源盛花期类群划分

    Table  3   Classification of blooming period groups of 84 walnut germplasm resources

    花期类群 种质资源 数量
    雄花早花
    (4月17日—4月26日)
    ‘A19’‘10901’‘扎343’‘绿岭’‘鲁果5号’‘硕宝’‘新早丰’‘丰香’‘鲁果7号’‘辽宁4号’‘魁香’‘辽宁1号’‘礼品1号’‘寒丰’‘新新2’‘辽宁6号’‘实生10’‘实生14’‘清香’‘硕香’‘新光’‘实生2’‘金薄丰’‘中核2号’‘西岭’‘A53’‘旮旯铺’‘实生13’‘花生1号’‘岱香’‘纸皮扎01’‘麻核桃’‘红光’‘安边2号’‘2-34’‘实生1’‘辽宁7号’‘龙珠’‘实生4’‘Pedro’‘B76’‘实生8’‘辽瑞丰’‘实生12’‘实生15’‘西洛7号’‘磨盘’‘中核1号’‘实生6’‘辽宁5号’‘扎木台3号’‘辽宁10号’‘纸皮3号’‘鲁果6号’ 54
    雄花中花
    (4月27日—5月6日)
    ‘西洛9号’‘礼品2号’‘美香’‘早硕’‘陇南11号’‘极早丰’‘西林1号’‘中核3号’‘绿波’‘元宝’‘晋香’‘金薄香2号’‘金薄香4号’‘西林5号’‘华山5号’‘青林’‘实生5’‘庆阳11号’‘晋丰’‘实生9’‘4-22’‘实生16’‘实生11’‘实生3’‘绿香’‘扎木台5号’‘实生7’‘核桃楸’ 28
    雄花晚花
    (5月7日及之后)
    ‘小果黑核桃’‘黑核桃’ 2
    雌花早花
    (4月17日—4月26日)
    ‘中核3号’‘西林1号’‘中核1号’‘晋香’‘美香’‘早硕’‘实生6’‘陇南11号’‘极早丰’‘鲁果6号’‘辽宁5号’‘辽宁10号’‘辽瑞丰’‘绿波’‘西洛9号’‘金薄香2号’‘礼品2号’‘纸皮3号’‘青林’‘实生7’‘庆阳11号’‘扎木台3号’‘元宝’‘金薄香4号’‘绿香’‘扎木台5号’‘西林5号’‘实生5’‘晋丰’‘实生9’‘实生16’‘实生11’‘辽宁4号’‘10901’‘实生15’‘中核2号’‘西岭’‘鲁果7号’‘辽宁1号’‘新新2’‘A19’‘辽宁6号’‘纸皮扎01’‘实生12’‘硕香’‘实生3’‘鲁果5号’‘新早丰’‘丰香’‘红光’‘礼品1号’‘B76’‘花生1号’‘核桃楸’‘岱香’ 55
    雌花中花
    (4月27日—5月6日)
    ‘新光’‘金薄丰’‘辽宁7号’‘A53’‘华山5号’‘寒丰’‘实生10’‘实生14’‘绿岭’‘4-22’‘实生13’‘实生1’‘魁香’‘Pedro’‘扎343’‘实生2’‘西洛7号’‘龙珠’‘硕宝’‘安边2号’‘磨盘’‘旮旯铺’‘实生8’‘实生4’‘清香’‘2-34’‘麻核桃’ 27
    雌花晚花
    (5月7日及之后)
    ‘小果黑核桃’‘黑核桃’ 2
    下载: 导出CSV
  • [1] 黄梅, 任华东, 姚小华, 等. 薄壳山核桃主要生物活性成分及其作用研究进展[J]. 中国油脂, 2023, 48(6): 99−104.

    Huang M, Ren H D, Yao X H, et al. Advances in studies on main bioactive components and their effects of pecans[J]. China Oils and Fats, 2023, 48(6): 99−104.

    [2] 傅本重, 邹路路, 朱洁倩, 等. 中国核桃生产现状与发展思路[J]. 江苏农业科学, 2018, 46(18): 5−8.

    Fu B Z, Zou L L, Zhu J Q, et al. Production status and development of Chinese walnut[J]. Jiangsu Agricultural Sciences, 2018, 46(18): 5−8.

    [3] 怀婷婷, 卫伟, 刘春晓, 等. 核桃产业和贸易现状及发展建议[J]. 安徽农业科学, 2023, 51(18): 203−208. doi: 10.3969/j.issn.0517-6611.2023.18.047

    Huai T T, Wei W, Liu C X, et al. Walnut industry and trade status and development suggestions[J]. Journal of Anhui Agricultural Sciences, 2023, 51(18): 203−208. doi: 10.3969/j.issn.0517-6611.2023.18.047

    [4] 刘朝斌, 李荣, 陈诗婷, 等. 核桃开花相关基因JrSOC1的筛选及其在开花中的作用分析[J]. 西北林学院学报, 2023, 38(4): 97−103. doi: 10.3969/j.issn.1001-7461.2023.04.12

    Liu C B, Li R, Chen S T, et al. Identification of flowering related JrSOC1 gene in walnut and its role in flowering[J]. Journal of Northwest Forestry University, 2023, 38(4): 97−103. doi: 10.3969/j.issn.1001-7461.2023.04.12

    [5] 崔元博, 储国林, 赵娟, 等. 薄壳山核桃主栽品种花期观察及品种配置[J]. 植物资源与环境学报, 2024, 33(3): 113−116. doi: 10.3969/j.issn.1674-7895.2024.03.12

    Cui Y B, Chu G L, Zhao J, et al. Flowering period observation and cultivar configuration of main pecan (Carya illinoinensis) cultivars[J]. Journal of Plant Resources and Environment, 2024, 33(3): 113−116. doi: 10.3969/j.issn.1674-7895.2024.03.12

    [6] 王彦平, 赵子威, 李洪峰. 大兴安岭地区不同海拔高度木本植物物候期及其对气候变化的响应[J]. 江西农业学报, 2022, 34(4): 172−179.

    Wang Y P, Zhao Z W, Li H F. Phenological characteristics of representative woody plants at different altitudes in Daxinganling region and their response to climate change[J]. Acta Agriculturae Jiangxi, 2022, 34(4): 172−179.

    [7]

    Bujdoso G, Soelimani A, Illes B, et al. Base temperature comparisons for leafing date, pistillate flower receptivity, and pollen shedding in Persian walnut[J]. Erwerbs-Obstbau, 2022, 64(4): 549−558. doi: 10.1007/s10341-022-00724-4

    [8] 郑景云, 葛全胜, 赵会霞. 近40年中国植物物候对气候变化的响应研究[J]. 中国农业气象, 2003(1): 29−33. doi: 10.3969/j.issn.1000-6362.2003.01.009

    Zheng J Y, Ge Q S, Zhao H X. Changes of plant phenological period and its response to climate change for the last 40 years in China[J]. Chinese Journal of Agrometeorology, 2003(1): 29−33. doi: 10.3969/j.issn.1000-6362.2003.01.009

    [9] 李恒, 张春花, 顾国栋, 等. 大泡核桃与圆菠萝核桃在攀枝花的引种迁移适应性[J]. 中国南方果树, 2023, 52(6): 196−199.

    Li H, Zhang C H, Gu G D, et al. Adaptability of introduction and migration of big bubble walnut and Yuanboluo walnut in Panzhihua[J]. South China Fruits, 2023, 52(6): 196−199.

    [10]

    Zaitchik B F, Evans J P, Smith R B. Regional impact of an elevated heat source: the Zagros Plateau of Iran[J]. Journal of Climate, 2007, 20(16): 4133−4146. doi: 10.1175/JCLI4248.1

    [11]

    Cosmulescu S, Baciu A, Botu M, et al. Environmental factors’ influence on walnut flowering[J]. Acta Horticulturae, 2010(861): 83−88.

    [12]

    Hájková L, Možný M, Oušková V, et al. Evaluation of walnut tree flowering and frost occurrence probability during 1961−2012[J]. Acta Universitatis Agriculturae Et Silviculturae Mendelianae Brunensis, 2022, 70(3): 235−248. doi: 10.11118/actaun.2022.018

    [13]

    Menzel A, Sparks T H, Estrella N, et al. European phenological response to climate change matches the warming pattern[J]. Global Change Biology, 2006, 12(10): 1969−1976. doi: 10.1111/j.1365-2486.2006.01193.x

    [14] 陈金海, 罗仁仙. 晚霜与冰雪灾害对核桃产业的影响及应对措施[J]. 经济林研究, 2011, 29(2): 119−123. doi: 10.3969/j.issn.1003-8981.2011.02.022

    Chen J H, Luo R X. Influence of late frost and freezing damage on walnut industry and some countermeasures[J]. Non-wood Forest Research, 2011, 29(2): 119−123. doi: 10.3969/j.issn.1003-8981.2011.02.022

    [15] 王凤妞, 姜涛, 李明媛, 等. 早实核桃‘西岭’成花调控相关基因和miRNA的发掘[J/OL]. 分子植物育种, 2024. [2024−04−16]. https://link.cnki.net/urlid/46.1068.S.20240415.0839.002.

    Wang F N, Jiang T, Li M Y, et al. Isolation of flowering-related genes and miRNAs in early bearing walnut cultivar ‘Xiling’[J/OL]. Molecular Plant Breeding, 2024. [2024−04−16]. https://link.cnki.net/urlid/46.1068.S.20240415.0839.002.

    [16]

    Bernard A, Marrano A, Donkpegan A, et al. Association and linkage mapping to unravel genetic architecture of phenological traits and lateral bearing in Persian walnut (Juglans regia L.)[J]. BMC Genomics, 2020(21): 203−227.

    [17] 王磊, 王红梅, 贾君, 等. 铁线莲属植物品种及野生种的物候期观察分析[J]. 浙江农业科学, 2016, 57(2): 198−201.

    Wang L, Wang H M, Jia J, et al. Phenophase observation and analysis of Clematis varieties and wild species[J]. Journal of Zhejiang Agricultural Sciences, 2016, 57(2): 198−201.

    [18] 李阳, 史高川, 王晋峰, 等. 山西晋中地区避晚霜核桃品种的筛选[J]. 果树资源学报, 2023, 4(6): 1−5.

    Li Y, Shi G C, Wang J F, et al. Selection of walnut cultivars avoiding late frost in Jinzhong Area of Shanxi Province[J]. Journal of Fruit Resources, 2023, 4(6): 1−5.

    [19] 张绥林, 李洋, 李琰, 等. 核桃晚霜危害特性及影响机制研究进展[J]. 中国农业科技导报, 2024, 26(4): 18−26.

    Zhang S L, Li Y, Li Y, et al. Research progress on harmful characteristics and mechanism of walnut late frost[J]. Journal of Agricultural Science and Technology, 2024, 26(4): 18−26.

    [20] 王企, 曹尚银, 陈利娜, 等. 核桃雌雄异熟性研究进展[J]. 中国果树, 2018(3): 65−71.

    Wang Q, Cao S Y, Chen L N, et al. Research progress on dichogamy of walnut[J]. China Fruits, 2018(3): 65−71.

    [21] 陈虹, 朱小虎, 黄学芹, 等. 不同早实核桃品种物候期观察与低温抗性评价[J]. 新疆农业大学学报, 2010, 33(6): 479−483. doi: 10.3969/j.issn.1007-8614.2010.06.003

    Chen H, Zhu X H, Huang X Q, et al. Observation on phenological period of different precocious walnut cultivars and evaluation of its cold tolerance in Aksu area of Xinjiang[J]. Journal of Xinjiang Agricultural University, 2010, 33(6): 479−483. doi: 10.3969/j.issn.1007-8614.2010.06.003

    [22] 高露璐, 李林芳, 马育珠, 等. 铁线莲品种群的花期观赏性状分析[J]. 园艺学报, 2017, 44(5): 921−932.

    Gao L L, Li L F, Ma Y Z, et al. Analysis of florescence characteristics on Clematis cultivars group[J]. Acta Horticulturae Sinica, 2017, 44(5): 921−932.

    [23] 何的明, 傅松玲, 毕慧慧, 等. 江淮地区薄壳山核桃主栽品种开花物候观测[J]. 经济林研究, 2022, 40(4): 42−51.

    He D M, Fu S L, Bi H H, et al. Observation on flowering phenology of main varieties of pecan in the area between Yangtze River and Huaihe River[J]. Non-wood Forest Research, 2022, 40(4): 42−51.

    [24] 缑锋利. 3种核桃花粉生活力的比较研究[J]. 山西农业科学, 2014, 42(2): 140−142. doi: 10.3969/j.issn.1002-2481.2014.02.10

    Gou F L. Pollen viability comparison between three walnuts[J]. Journal of Shanxi Agricultural Sciences, 2014, 42(2): 140−142. doi: 10.3969/j.issn.1002-2481.2014.02.10

    [25] 魏海林, 欧阳群文, 蒋瑶. 薄壳山核桃成花坐果的主要机理、影响因素及其促进措施[J]. 湖南农业科学, 2021(8): 106−111.

    Wei H L, Ouyang Q W, Jiang Y. Review on the main mechanism and factors affecting the flowering and fruit-setting of Carya illinoensis and its technical measures of high yield cultivation[J]. Hunan Agricultural Sciences, 2021(8): 106−111.

    [26] 刘晓芬, 凌晓祺, 向理理, 等. 温度和赤霉素对春兰开花的调控[J]. 浙江农业学报, 2023, 35(2): 9. doi: 10.3969/j.issn.1004-1524.2023.02.13

    Liu X F, Ling X Q, Xiang L L, et al. Effect of temperature and gibberellin on flowering regulation of Cymbidium goeringii[J]. Acta Agriculturae Zhejiangensis, 2023, 35(2): 9. doi: 10.3969/j.issn.1004-1524.2023.02.13

  • 期刊类型引用(13)

    1. 陈仕媛,马岚,陈佩岩. 北方土石山区不同粒径团聚体特征及其对坡面侵蚀过程的影响. 水土保持研究. 2025(02): 102-110+139 . 百度学术
    2. 张紫薇,赵文俊,李奇,马一淳,田亮,杨广勇,李忠武,刘窑军. 表层土壤结构对红壤坡耕地产流及产沙过程的影响. 土壤学报. 2024(02): 434-444 . 百度学术
    3. 王琳涛 ,王云琦 ,王玉杰 ,张建聪 ,尹乐 ,张洪珊 . 火干扰对缙云山土壤团聚体稳定性的影响. 中国水土保持科学(中英文). 2023(03): 25-33 . 百度学术
    4. 孙旭,马岚,张金阁,王福星,薛梦华,宋天翼,张旖璇,李浩玉. 模拟降雨条件下枯落物覆盖对坡面地表和壤中流产流过程的影响. 陕西师范大学学报(自然科学版). 2022(04): 136-144 . 百度学术
    5. 刘文景,许文年,夏露,陈剑刚,程虎,余飞,夏栋,刘大翔,余绍鹏. 鄂西长江小流域土地利用方式对土壤团聚体稳定性的影响. 长江流域资源与环境. 2021(04): 925-935 . 百度学术
    6. 郭岩松,毕曼. 土壤侵蚀影响因素研究进展. 河南科技. 2021(20): 136-140 . 百度学术
    7. 黎俊佑,马岚,刘京晶,沈雨珺,张金阁,孙占薇,燕琳,张栋. 雨型对华北土石山区坡面土壤侵蚀的影响. 水土保持学报. 2020(04): 1-6 . 百度学术
    8. 胡昭,张懿,郭建. 黄土丘陵区典型林地土壤团聚体稳定性评价. 人民黄河. 2020(08): 104-108 . 百度学术
    9. 郭威震,胡斐南,谭滔滔,马任甜,刘婧芳,李喆,赵世伟. 土壤表面电场对黑土团聚体破碎和侵蚀的影响. 应用生态学报. 2020(08): 2644-2652 . 百度学术
    10. 张金阁,马岚,张栋,刘京晶,黎俊佑,孙占薇,燕琳,孙旭,王福星. 华北土石山区草地过滤带布设方式对坡面产流产沙的影响. 水土保持学报. 2020(05): 76-80 . 百度学术
    11. 毛旭芮,王明力,杨建军,丁肇龙,曹月娥. 采煤对露天煤矿土壤理化性质及可蚀性影响. 西南农业学报. 2020(11): 2537-2544 . 百度学术
    12. 李肖,陈晨,林杰,朱茜,董波,丁鸣鸣. 侵蚀强度对淮北土石山区土壤团聚体组成及稳定性的影响. 水土保持研究. 2019(04): 56-61+67 . 百度学术
    13. 王婷,李建平,张翼,井乐,张茹. 不同降水下天然草地土壤水稳定性团聚体分布特征. 草业科学. 2019(08): 1935-1943 . 百度学术

    其他类型引用(13)

图(4)  /  表(3)
计量
  • 文章访问数:  113
  • HTML全文浏览量:  17
  • PDF下载量:  29
  • 被引次数: 26
出版历程
  • 收稿日期:  2024-08-20
  • 修回日期:  2024-11-13
  • 录用日期:  2025-01-06
  • 网络出版日期:  2025-01-08
  • 刊出日期:  2025-02-24

目录

/

返回文章
返回