论文部分内容阅读
摘要 为了确定崂山茶叶致灾的农业气象指标，把崂山茶叶多年产量资料分为“歉年”与“非歉年”2种年景，利用SPSS22.0统计分析系统，首先选取影响显著的气象因子，再进行二元logistic回归分析确定主要影响因子。结果表明，冬季日最低气温连续2 d低于-10℃是崂山茶叶遭受冻害的主要农业气象指标，低温冻害是崂山茶叶减产和品质下降的主要氣象灾害。
关键词 崂山茶;年景分析;气象因子;二元logistic回归
中图分类号：S571.1;S42 文献标识码：A 文章编号：2095-3305（2019）02-086-02
DOI： 10.19383/j.cnki.nyzhyj.2019.02.036
Abstract In order to determine the agrometeorological indicators of tea disasters in Laoshan Mountain， the multi?year yield data of Laoshan tea were divided into apologetic year and non?apologetic year. Using SPSS22.0 statistical analysis system， the meteorological factors with significant impacts were selected first， and then the main influencing factors were determined by binary logistic regression analysis. The results showed that the main agrometeorological index of Laoshan tea freezing injury was that the daily minimum temperature was lower than -10℃ for two consecutive days in winter， and low temperature freezing injury was the main meteorological disaster of Laoshan tea yield reduction and quality decline.
Key words Laoshan tea; Year’s harvest analysis; Meteorological factors; Binary logistic regression
崂山位于山东省青岛市东部，被黄海环绕，一年四季气候湿润、多云雾，土壤呈弱酸性。崂山茶自1959年开始南茶北引到20世纪80年代初引种成功，20世纪90年代初开始推广种植，随着栽培和管理技术的不断提高，崂山茶叶的氨基酸、咖啡碱等有机营养成分含量较高，崂山茶品质优良，深受人们的青睐。有学者研究了气候条件对茶叶产量与品质的影响，如李柏贞等指出茶叶的气象灾害主要是冬季的冻害、春季的霜冻和夏季的热害[1]，赖全康等调查分析认为冬季的冰冻灾害导致茶叶受损严重[2]。但目前气象灾害对崂山茶叶的影响研究较少，因此，笔者以崂山茶叶的产量资料为样本进行研究，探讨气象因素对崂山茶叶产量年景的影响。
1 资料与方法
1.1 资料来源
以崂山区统计局1994—2015年王哥庄街道平均单产为样本，气象资料为崂山区气象局国家一般气象站对应的资料。与茶叶生产相关的气象资料分别为年平均气温（TT）、年平均最低气温（TD）、年平均最高气温（TG） 、年平均气温日较差（TGD）、 年平均水汽压（EE）、年降雨量（RR）、日地面0 cm温度是否出现连续2 d低于-10℃的日数（TD0，没有出现为0，出现了为1）;日最低气温是否出现连续2 d低于-10℃的日数（TZD，没有出现为0，出现了为1）。
TD0、TZD及产量的“歉年”与“非歉年”为类别变量，没有单位;其他资料为数值型变量，单位为对应的国际单位;产量年景“歉年”与“非歉年”为因变量，气象资料为自变量。
1.2 产量资料处理
1.2.1 相对产量及产量的相对偏差 其计算公式[3]如下：
相对产量=×100%（1）
相对偏差=×100%（2）
式中，xi为某一年崂山茶叶的产量，x为22年崂山茶的平均产量。
1.2.2 茶叶产量的年景标准 相对偏差值大于-8%的年份定为“非歉年”，相对偏差值小于或等于-8%的产量年份定为“歉年”。据此计算出历史上的茶叶每年年景[4]。
2 结果与分析
2.1 气象因子的显著性检验
2.1.1 单因素检验 检验数值型变量TT、TG、TD、TGD、EE、RR在茶叶产量“歉年”与“非歉年”的显著性差异采用t检验的方法。结果表明，茶叶产量“歉年”和“非歉年”在TT、TG、TD、TGD、EE、RR各因子上不存在显著差异，即年度气象因子TT、TG、TD、TGD、EE、RR不会对崂山茶的“歉年”与“非歉年”造成影响。
2.1.2 类别变量TD0、TZD的检验 通过分析-描述性统计-交叉表分析-卡方检验-fisher精确检验，结果表明，不同的产量状况在以上的2个因子上全部存在显著差异（P<0.05），说明茶叶产量年景的“歉年”与“非歉年”2种年景最终是由TD0、TZD共同影响造成的。
2.2 回归模型建立
据计算，崂山茶叶致灾的二元logistic回归模型的函数表达式为：
P=（3）
2016年1月23、24日，青岛市崂山区出现了强寒潮天气，实测日最低气温连续2 d分别为-12.9、-15.5℃，即出现了TZD=1，因此2016年出现“歉年”的概率为96.77%。实地调查结果表明大约有30%的大田崂山茶因遭受冻害致死。
3 结论与讨论
冬季当日最低气温连续2 d出现低于-10℃以下，此时土壤处于结冰状态，土壤水分移动和上升受到阻碍，茶树体内细胞结冰，导致茶树遭受生理冻害，会使茶树整体死亡[5]。崂山区每年最冷月份为12月到翌年2月，日平均气温一般在-3.5～0℃，如果出现长时间低温，连续低于0℃的天数大于15 d，负积温总和超过-100℃，最低气温低于-10℃，崂山茶树易遭受冻害致死。崂山茶树的主要气象灾害是冻害[6]。为此，可以根据冬季寒潮强度及时发布寒潮预警，与崂山区政府各部门联合开展崂山茶树防灾气象服务，通过烟熏、灌溉、覆盖、扣棚等措施，把崂山茶树造成的损失降到最低。
参考文献
[1] 李柏贞，谢佳杏，孔萍，等.江南茶叶农业气象灾害风险区划[J].干旱气象，2015，33（6）：1017-1023.
[2] 赖全康，郭忠春，杨平仲，等.2008年三江县茶叶冰冻灾害调查报告[J].广西农学报，2008，（S1）：81-84.
[3] 欧阳海，郑步忠，王雪娥，等.农业气候学[M].北京：气象出版社，1990：37-42.
[4] 中国气象局.主要粮食作物产量年景等级.QX/T 335-2016[S]. 中国气象局，2016.
[5] 冯秀藻，陶炳炎.农业气象学原理[M].北京：气象出版社，1991：95-96，142-143.
[6] 朱良民.崂山茶[M].青岛：中国海洋大学出版社，2006：32.
责任编辑：郑丹丹
关键词 崂山茶;年景分析;气象因子;二元logistic回归
中图分类号：S571.1;S42 文献标识码：A 文章编号：2095-3305（2019）02-086-02
DOI： 10.19383/j.cnki.nyzhyj.2019.02.036
Abstract In order to determine the agrometeorological indicators of tea disasters in Laoshan Mountain， the multi?year yield data of Laoshan tea were divided into apologetic year and non?apologetic year. Using SPSS22.0 statistical analysis system， the meteorological factors with significant impacts were selected first， and then the main influencing factors were determined by binary logistic regression analysis. The results showed that the main agrometeorological index of Laoshan tea freezing injury was that the daily minimum temperature was lower than -10℃ for two consecutive days in winter， and low temperature freezing injury was the main meteorological disaster of Laoshan tea yield reduction and quality decline.
Key words Laoshan tea; Year’s harvest analysis; Meteorological factors; Binary logistic regression
崂山位于山东省青岛市东部，被黄海环绕，一年四季气候湿润、多云雾，土壤呈弱酸性。崂山茶自1959年开始南茶北引到20世纪80年代初引种成功，20世纪90年代初开始推广种植，随着栽培和管理技术的不断提高，崂山茶叶的氨基酸、咖啡碱等有机营养成分含量较高，崂山茶品质优良，深受人们的青睐。有学者研究了气候条件对茶叶产量与品质的影响，如李柏贞等指出茶叶的气象灾害主要是冬季的冻害、春季的霜冻和夏季的热害[1]，赖全康等调查分析认为冬季的冰冻灾害导致茶叶受损严重[2]。但目前气象灾害对崂山茶叶的影响研究较少，因此，笔者以崂山茶叶的产量资料为样本进行研究，探讨气象因素对崂山茶叶产量年景的影响。
1 资料与方法
1.1 资料来源
以崂山区统计局1994—2015年王哥庄街道平均单产为样本，气象资料为崂山区气象局国家一般气象站对应的资料。与茶叶生产相关的气象资料分别为年平均气温（TT）、年平均最低气温（TD）、年平均最高气温（TG） 、年平均气温日较差（TGD）、 年平均水汽压（EE）、年降雨量（RR）、日地面0 cm温度是否出现连续2 d低于-10℃的日数（TD0，没有出现为0，出现了为1）;日最低气温是否出现连续2 d低于-10℃的日数（TZD，没有出现为0，出现了为1）。
TD0、TZD及产量的“歉年”与“非歉年”为类别变量，没有单位;其他资料为数值型变量，单位为对应的国际单位;产量年景“歉年”与“非歉年”为因变量，气象资料为自变量。
1.2 产量资料处理
1.2.1 相对产量及产量的相对偏差 其计算公式[3]如下：
相对产量=×100%（1）
相对偏差=×100%（2）
式中，xi为某一年崂山茶叶的产量，x为22年崂山茶的平均产量。
1.2.2 茶叶产量的年景标准 相对偏差值大于-8%的年份定为“非歉年”，相对偏差值小于或等于-8%的产量年份定为“歉年”。据此计算出历史上的茶叶每年年景[4]。
2 结果与分析
2.1 气象因子的显著性检验
2.1.1 单因素检验 检验数值型变量TT、TG、TD、TGD、EE、RR在茶叶产量“歉年”与“非歉年”的显著性差异采用t检验的方法。结果表明，茶叶产量“歉年”和“非歉年”在TT、TG、TD、TGD、EE、RR各因子上不存在显著差异，即年度气象因子TT、TG、TD、TGD、EE、RR不会对崂山茶的“歉年”与“非歉年”造成影响。
2.1.2 类别变量TD0、TZD的检验 通过分析-描述性统计-交叉表分析-卡方检验-fisher精确检验，结果表明，不同的产量状况在以上的2个因子上全部存在显著差异（P<0.05），说明茶叶产量年景的“歉年”与“非歉年”2种年景最终是由TD0、TZD共同影响造成的。
2.2 回归模型建立
据计算，崂山茶叶致灾的二元logistic回归模型的函数表达式为：
P=（3）
2016年1月23、24日，青岛市崂山区出现了强寒潮天气，实测日最低气温连续2 d分别为-12.9、-15.5℃，即出现了TZD=1，因此2016年出现“歉年”的概率为96.77%。实地调查结果表明大约有30%的大田崂山茶因遭受冻害致死。
3 结论与讨论
冬季当日最低气温连续2 d出现低于-10℃以下，此时土壤处于结冰状态，土壤水分移动和上升受到阻碍，茶树体内细胞结冰，导致茶树遭受生理冻害，会使茶树整体死亡[5]。崂山区每年最冷月份为12月到翌年2月，日平均气温一般在-3.5～0℃，如果出现长时间低温，连续低于0℃的天数大于15 d，负积温总和超过-100℃，最低气温低于-10℃，崂山茶树易遭受冻害致死。崂山茶树的主要气象灾害是冻害[6]。为此，可以根据冬季寒潮强度及时发布寒潮预警，与崂山区政府各部门联合开展崂山茶树防灾气象服务，通过烟熏、灌溉、覆盖、扣棚等措施，把崂山茶树造成的损失降到最低。
参考文献
[1] 李柏贞，谢佳杏，孔萍，等.江南茶叶农业气象灾害风险区划[J].干旱气象，2015，33（6）：1017-1023.
[2] 赖全康，郭忠春，杨平仲，等.2008年三江县茶叶冰冻灾害调查报告[J].广西农学报，2008，（S1）：81-84.
[3] 欧阳海，郑步忠，王雪娥，等.农业气候学[M].北京：气象出版社，1990：37-42.
[4] 中国气象局.主要粮食作物产量年景等级.QX/T 335-2016[S]. 中国气象局，2016.
[5] 冯秀藻，陶炳炎.农业气象学原理[M].北京：气象出版社，1991：95-96，142-143.
[6] 朱良民.崂山茶[M].青岛：中国海洋大学出版社，2006：32.
责任编辑：郑丹丹