武强， 王旭， 方丽， 江姣，孙恩虹， 韩旭， 陈思英

1.重庆市气象科学研究所，重庆 401147；2.重庆市巴南区气象局，重庆 401320；3.重庆市长寿区气象局，重庆 长寿 401220；4.重庆市江津现代农业气象试验站，重庆 江津 402260；5.安徽省科技研究开发中心，合肥 230088

茶树是典型的亚热带常绿植物，具有喜温、喜湿，且耐阴、喜散射光的生长特性.气象因子是决定茶树的物候期、鲜叶适采期的重要因素[1-2]，同时影响着茶叶品质形成的生理生化过程.在茶叶生产中，对茶叶物候期尤其是开采期的准确预判具有重要的意义，新梢嫩芽若不及时采摘，会叶片老化，品质下降，适制性范围变小，严重影响成茶经济效益.

作物物候期的气象预测方法中，基于阶段积温学说的积温预报法应用十分广泛，即作物发育主要是受温度的影响，完成某一发育阶段所需的积温基本一定[3].基于不同下限温度的活动积温和有效积温在阶段积温与作物物候期的相关关系研究中[4]，以及作物模型如林果类花期、成熟期的物候期预报[5-6]中有广泛的应用.逐步回归统计法通过分析天气气候条件与作物生育期相关关系，筛选关键气象因子，建立预测模型，实现作物物候期预测[7-8]，在观赏类植物的花期气象预报中得到应用，进而为旅游景区的打造和规划提供参考依据[9-10].集成预报法对多种预报方法预测结论进行综合集成，构建新的预报模型，可以优化单一预报方法预报结果不稳定的问题，达到提高预测准确率的目的.该建模思想在天气模式、环境气象以及农业气象研究中均有应用，并取得较好的预报效果[11-13].

对于茶叶开采期的预测研究，孙秀邦等[14]基于积温法对安徽泾县春茶开采期进行预测，朱兰娟等[15]研究了西湖龙井开采期的影响因子与预报方法，较为准确地对西湖龙井开采期进行了预报.特早种茶树春茶开采期较常规品系明显偏早，因而开采期对气象条件的需求与常规品系也有一定的差异.本研究以特早种茶树品种“巴渝特早”为研究对象，以开采期关键气象影响因子以及预测方法研究为切入点，应用阶段积温预报法、逐步回归预报法、集成预报法分别建立“巴渝特早”开采期气象预报模型，比较不同预报方法对“巴渝特早”开采期的预报效果，旨在为特早种茶树的开采期预报和关于春茶开采的农业管理决策提供依据.

1 材料与方法

1.1 研究区域及品种

本研究区域为“巴渝特早”品系比较试验所在地重庆市巴南区二圣镇白象山(29°30′N，106°42′E，海拔400～700 m)，气候条件优越，年≥10 ℃活动积温为5 600 ℃·d左右，茶树生长期内常年无冻害发生，且春季气温回暖早，土壤水分适宜，空气湿度大，气候条件适宜茶树生长.

研究品种为“巴渝特早”，鉴定编号为国品鉴茶2014001[16-17]，具有开采期特早的优点.巴南白象山种植区的多年观测资料显示，开采期在每年3月上旬至中旬.

1.2 资料来源及处理

茶叶开采期观测资料(2009-2020年)来自于重庆茶叶集团有限公司.将春茶开采日期转换为年日序值，即1月1日记为日序1，1月2日记为日序2，以此类推.

气象资料来源于巴南区国家基本气象站、区域自动气象站的逐日气象数据，包括：平均气温、极端最高气温、极端最低气温、降雨量、日照时数、空气相对湿度、最小空气相对湿度.资料年限结合开采期观测资料，选用时段为2009年1月-2020年4月.将茶叶开采期观测资料与气象资料按年限分为两类数据样本集，2009-2018年为校正样本，用于建立茶叶开采期预报模型；2019-2020年为验证样本，用于验证所建立茶叶开采期预报模型的预报效果.

1.3 预报方法

阶段积温预报法：计算历年春茶开采期≥0 ℃积温以及不同下限温度的活动积温、有效积温，比较得出变异系数最小的积温类型，计算多年平均值作为预报标准.待预测年份从当年1月1日，即日序1开始统计，结合前期气象条件以及后期温度预报，以达到积温预报标准的日序作为茶叶预测开采期.

逐步回归预报法：通过相关分析法分别计算开采期前各类气象因子与开采期年日序相关系数，筛选出相关系数在0.05水平有统计学意义的气象因子作为预报初选气象因子.将初选气象因子与开采期年日序通过逐步回归构建开采期气象预报模型，在有统计学意义的所有预报方案中，选择预报效果最优的作为茶叶开采期预报模型.利用自变量气象因子计算得到茶叶预测开采期.

集成预报法：对阶段积温预报法与逐步回归预报法得到的日序值通过多元线性回归得到集成预报模型.通过利用阶段积温预报法和逐步回归预报法得到待预测年份日序值，代入集成预报模型获得最终的茶叶预测开采期.

2 结果与分析

2.1 “巴渝特早”开采期气象条件分析

2.1.1 热量条件

“巴渝特早”属于特早种茶树品种，在茶树遗传特性基础上，热量、水分、光照等气象条件共同决定了茶芽春季萌动的日期，其中，热量条件是茶芽萌动的主要因素.一般认为茶树春季生长温度为日平均气温稳定通过10 ℃，新梢开始缓慢生长，随着活动积温的增加，春季第一片新叶展开，进入春茶开采期[18].

下限温度指标是积温计算过程中的必要要素，且会随作物品种和生育期不同而不同，相同下限温度以活动积温和有效积温进行作物生育期的积温统计，或者使用单一积温统计方法而选取不同下限温度得到积温统计值[19-20].为明确“巴渝特早”春梢生长的热量需求，结合2009-2018年“巴渝特早”开采期，除统计当年≥0 ℃积温以外，分别设定5，6，7，8，9，10 ℃的不同下限温度，进行阶段活动积温和有效积温的统计，筛选确定“巴渝特早”芽叶生长的下限温度[18].此外，将采摘期前一段时间的平均气温作为热量条件研究的指标，设置“巴渝特早”开采前15，10，5 d的平均气温.统计并计算各热量条件指标的多年平均值以及变异系数.经计算，2009-2018年，“巴渝特早”开采期前平均气温作为热量条件指标，以开采期前15 d平均气温变异系数较小，为0.138；有效积温与活动积温中，以5 ℃为下限温度的活动积温变异系数最小，为0.125；≥0 ℃积温平均值为564.4 ℃·d，变异系数为0.098，是所有热量条件指标中最小变异系数.随着不同年份间的差异，开采期前15 d平均气温年际差异在5 ℃左右，≥0 ℃积温、≥5 ℃活动积温总体变化较为平稳，个别年份(2011，2016，2018年)较多年平均值偏低50～100 ℃·d(图1).2015-2017年，开采期前15 d的平均气温变化与积温变化呈现相反的趋势，2015年和2017年积温值较高，但开采期前15 d气温明显偏低，主要是因为当年气温条件整体较好，前期积温值偏高，茶芽达到萌发的热量要求，但是在开采期前出现一段降温过程，导致开采期前15 d气温平均值偏低，但积温总量已经满足茶树芽叶生长需求，达到开采标准.2016年的情况与之相反，当年气温总体偏低，积温条件偏差，但是在开采期前遇升温过程，导致开采期前15 d平均气温偏高.

图1 “巴渝特早”开采期热量条件

2.1.2 影响气象因子

因“巴渝特早”开采期一般在3月上旬至中旬，故选取2009-2018年1-3月各月、旬的平均气温、平均日最高气温、平均日最低气温、平均空气相对湿度、平均日最低空气相对湿度、总降水量、总日照时数，分别与开采期日序进行相关分析，筛选出在0.01水平有统计学意义的气象因子作为初选因子.影响“巴渝特早”开采期的气象因子主要集中于2月份，温度与开采期日序的相关性高于空气相对湿度和降水，日照时数与开采期日序无统计学意义(表1).白象山茶叶产区常年2-5月气温资料显示，进入2月中旬以后，日平均气温在8 ℃左右，与多年观测到的独芽开采期日平均气温一致，但是早春气温回升以及年际差异波动较大.进入3月上旬，日平均气温升高到10 ℃以上，此时气温一般较为稳定，新梢开始缓慢生长，随着活动积温的增加，春季第一片新叶展开，进入春茶采摘期.统计历年2-5月0～50 cm土层平均土壤相对含水量以及空气相对湿度，发现白象山茶叶产区水分条件较为稳定，且满足茶叶生长的水分需求，不构成春茶开采限制因素.而茶树原生环境为亚热带丛林，为高大乔木遮盖下的灌木群落，具有耐阴的生理特性，故光照条件对开采期影响并不显著.

表1 与“巴渝特早”开采期日序显著相关的气象因子筛选结果

2.2 “巴渝特早”开采期预报模型的建立与验证

2.2.1 阶段积温预报法

统计2009-2018年春茶开采前各热量指标，变异系数最小的因子为≥0 ℃积温，变异系数为0.098.多年平均值为564.4 ℃·d.因此，采用≥0 ℃积温值作为“巴渝特早”开采期预报标准.从1月1日开始将≥0 ℃的日平均气温进行累加，累加值达到564.4 ℃·d的日期作为“巴渝特早”开采期.

以2009-2018年开采期日序作为校正数据集，以2019，2020年开采期日序作为验证数据集，利用积温预报开采期，验证样本预测值与调查值偏差均为1 d，预测效果较好，但在校正样本中，预测值与调查值普遍偏差为5 d以内，个别年份偏差较大，如2011年的预测结果比实际情况偏晚14 d(表2).经调查，2011年1月气温持续偏低，月平均气温2.3 ℃，较常年同期偏低4.1 ℃，2011年1月积温71.2 ℃·d，较其余年份1月积温平均值偏低了127.5 ℃·d，此时茶芽尚在休眠期，对芽体生长影响不大.进入2月，气温迅速升高，茶芽生长较快，进入开采期，此时积温为453.0 ℃·d，较历年开采期积温偏低了111.4 ℃·d，可见2011年开采期积温偏低主要是1月的持续低温导致.说明利用≥0 ℃积温预报茶叶开采期，可能受到冬季茶树休眠期的异常气温影响，增加预测结果的不确定性.因此，虽然在验证样本中，对“巴渝特早”开采期的预测偏差较小，但是在未来茶树休眠期气温异常的年份，预测结果仍可能出现较大偏差.

表2 阶段积温预报结果与误差

2.2.2 逐步回归预报模型

以表1初选因子作为自变量，开采期年日序作为因变量进行逐步回归，建立预报模型，预报模型在0.05水平有统计学意义.分别采用2月、2月中旬、2月下旬的平均气温、平均日最高气温、平均日最低气温、平均空气相对湿度、平均日最低空气相对湿度、总降水量、总日照时数作为预报因子，以2009-2018年开采期日序作为校正数据集，以2019，2020年开采期日序作为验证数据集，建立逐步回归预报模型，经比较分析预测效果，筛选得到“巴渝特早”逐步回归开采期预报模型：

Y=124.883-7.212X1+11.223X2-5.053X3-0.544X4

式中：X1为2月中旬平均日最低气温；X2为2月中旬平均气温；X3为2月中旬平均日最高气温；X4为2月平均空气相对湿度.预报模型在0.05水平有统计学意义.

逐步回归预报模型在验证样本预测值与调查值时偏差分别为1.9，0.5 d，预测效果较好，在校正样本的预报偏差为3.5 d以内，但个别年份偏差较大，如2016年预报偏差偏晚7 d(表3).调查显示，2016年1月中上旬气温较高且稳定，平均值为7.3 ℃，利于特早种茶树萌发，在茶芽萌发以后，2月上旬升温较快，到2月上旬末平均气温大于10 ℃，气温升高促进新芽快速伸展，故使得开采期整体提前，以2月中旬气温为主要预报因子的预测结果出现一定偏差.

表3 逐步回归模型预报结果与误差

2.2.3 开采期集成预报模型

以2009-2018年积温预报和逐步回归预报模型预测结果的日序作为自变量，以开采期实际调查值作为因变量，采用多元回归方法确定积温预报法与逐步回归预报法的权重系数，构建集成预报模型，得到：

Y=1.63+0.127Y1+0.848Y2

式中：Y1为积温预报日序值；Y2为逐步回归预报模型预测日序值.

以2009-2018年开采期日序作为校正数据集，以2019，2020年开采期日序作为验证数据集，集成预报模型在验证样本中，预测值与调查值偏差分别为1.6，0.5 d，预测效果较好.校正样本的预报偏差为3.4 d以内，个别年份，如2016年预报偏差偏晚6.4 d，预报偏差也较逐步回归预报模型变小，预报精度有所提高(表4).

表4 “巴渝特早”开采期集成预报模型预报结果与误差

2.3 不同预报方法性能比较

比较3种“巴渝特早”开采期预报方法的预测结果与实际值的差异(图2)，阶段积温预报值与实际值拟合优度R2为0.539，与调查值一致性较好，说明积温预报方法能够较好地对开采期进行预测，结合实际预报结果与实际值的比较，发现预测结果可能存在一定的偏差；逐步回归模型预测值与实际值拟合优度R2为0.749，拟合效果较积温预报方法有明显的提升，能够对“巴渝特早”开采期进行更加准确地预报；集成预报模型预测值与实际值拟合优度R2为0.765，较逐步回归预测模型的准确性有所提升.

图2 预报模型预测值与实际值拟合比较

比较3种开采期预报模型性能指标(表5)，阶段积温预报的预报性能表现相对较差，集成预报模型预报性能表现最好.阶段积温预报结果最大偏差可达14 d，集成预报模型最大偏差为6.4 d，偏差明显缩小，排除掉个别异常偏差值以后，积温预报结果偏差在5 d以内，集成预报模型偏差在3.4 d以内，预测偏差的缩小在实际生产中具有重要意义，当春季新芽叶达到采摘标准以后，及时采摘才能保证茶叶的品质以及经济效益.利用阶段积温预报的结果残差平方和为273.0，表明阶段积温预报法对茶叶开采期预测的准确率的不确定性较大，而集成预报模型预测结果残差平方和仅为74.5，集成预报法对茶叶开采期预测的准确率的不确定性明显减小.阶段积温预报的结果剩余均方差为5.225，集成预报模型预测结果剩余均方差为2.729，能够更加准确地对“巴渝特早”开采期进行预测.

表5 预报模型性能指标

3 结论与讨论

比较3种“巴渝特早”开采期预报方法的预测结果与实际值的差异，阶段积温法、逐步回归法、集成预报法的预报值与实际值拟合优度R2分别为0.539，0.749，0.765，阶段积温预报结果最大偏差可达14 d，集成预报模型最大偏差为6.4 d，排除掉个别异常偏差值，积温预报结果偏差在5 d以内，集成预报模型偏差在3.4 d以内.阶段积温预报的预报性能表现相对较差，集成预报模型预报性能表现最好.

阶段积温预报要结合前期气温以及后期气温预报来做统计，对特早芽种的茶树来说，当年冬季异常偏暖可能导致茶芽提前萌发，1月1日前的偏高气温在积温统计中就有所遗漏，导致结果预测不准确；此外，积温预报仅考虑了热量条件对茶芽发育与叶片展开的影响，忽略了降水、光照等气象要素对生育期的综合影响，也会为预测结果增加一定的不确定性.

利用逐步回归方法建立“巴渝特早”开采期预报模型，相较于阶段积温预报法预测结果的准确性明显提升.选择多个关键气象因子进行逐步回归建模，对开采期真值的拟合效果更好，但是对于关键气象因子以外的气象因子的开采期影响考虑不足，会导致逐步回归模型预报的局限性.

集成预报模型融合阶段积温预报与逐步回归预报结果优化了阶段积温预报法与逐步回归预报法预测结果的不稳定性，一定程度规避了积温预报法个别年份较大的不确定性，同时较逐步回归模型预测结果的基础上有所提升，预测精度较高，适用于特早种茶叶的开采期气象预报.

本研究的开采期观测资料依然较少，随着观测资料的积累，可对各预测模型进行进一步优化.