文章信息
- 彭继达, 马治国, 杨凯
- PENG Jida, MA Zhiguo, YANG Kai
- 基于福建省冬季低温分型的台湾青枣寒冻害特征
- Cold and freezing injury characteristics of Taiwan-based Zizyphus mauritiana based on winter low temperature classification in Fujian Province
- 亚热带农业研究, 2017, 13(4): 270-275
- Subtropical Agriculture Research, 2017, 13(4): 270-275.
- DOI: 10.13321/j.cnki.subtrop.agric.res.2017.04.010
-
文章历史
- 收稿日期: 2017-10-13
台湾青枣(Zizyphus mauritiana Lam.)别名毛叶枣,为鼠李科枣属常绿小乔木,原产于印度、北非等地,是由印度毛叶枣经台湾育种家多代改良选育而成的优良水果品种[1-2]。台湾青枣产量较高且生产周期短,当年嫁接、种植,当年挂果,效益较好且稳定[3]。此外,台湾青枣营养价值丰富,具有保健价值和生态价值,已成为我国极具发展潜力的特色热带水果。近年来,福建省将台湾青枣作为“短、丰、快”的新兴果树种类加以推广,引种的台湾青枣主要分布在漳州市,占全省总种植面积的96.2%[4]。
福建省属亚热带季风气候,气候资源优越,但气象灾害也频繁发生,其中冬季强低温过程对亚热带作物造成较大影响[5]。冬季台湾青枣处于果实成熟的关键时期,遭受寒冻害则影响其产量和品质。国内外学者对于亚热带果树的低温寒害研究,主要集中在果树受寒害影响的机理及其诊断、品种差异及抗寒性特征方面[6-8],对青枣的研究主要集中在品种选育、引种风险分析及区划、栽培方法和寒冻害指标等方面[9-13]。有关青枣寒冻害研究,主要通过寒冻害形态学差异划分低温阈值,进而确定寒冻害指标,尚未对低温的类型差异进行针对性研究[14]。在同一低温不同低温类型下,果树的寒冻害程度有所不同[15-19],因此,有必要研究不同低温类型对青枣的影响差异。本研究根据气象观测资料和台湾青枣寒冻害灾情资料,研究不同低温类型下台湾青枣的寒冻害特征差异,以期为合理评估和防御寒冻害提供参考。
1 材料与方法 1.1 台湾青枣寒冻害灾情资料及灾害等级标准根据福建省历史低温分布规律,在青枣主要种植区选取9个具有代表性天气类型的试验点,分别为长泰县、云霄县、平和县、龙海市、罗源县、永泰县、福鼎市气象局观测场以及漳州天宝热作试验站和南靖县高港村试验点。2015—2016年冬季,分别在9个试验点各放置4盆台湾青枣,共36盆,记录寒冻害情况。同时,也记录试验期间各观测点附近青枣种植园寒冻害情况。2015—2016年冬季果园小气候自动站记录数据来自于福建省气象科学研究所。将9个试验点邻近果园小气候自动站寒冻害过程气温收集整理,得到台湾青枣寒冻害过程极端最低气温值(Td)。根据台湾青枣寒冻害形态学分级标准(表 1) [20],将树体、叶片、枝条、主干受害程度划分为0~4级,分别表示无、轻度、中度、重度、严重寒冻害。2015—2016年冬季,共有23例台湾青枣寒冻害灾情样本。
寒冻害等级1) | 果实 | 树叶 | 枝梢 |
0 | 果皮、果肉均不变色 | 叶片无变化 | 枝条无变化 |
1 | 果皮变色,果肉不变色 | 嫩叶卷曲 | 枝条无变化 |
2 | 果皮、果肉均变色,果实受冻率≤20% | 外围叶片出现卷曲 | 末端嫩梢受冻干枯,次级枝梢无冻伤 |
3 | 果实受冻率>20%,外围枝条受冻干枯 | 树体中上部及外围叶片卷曲、干枯、脱落 | 主枝部分冻伤 |
4 | 果实受冻率100%,烂果 | 叶片全部干枯、脱落 | 枝条、主枝均褐变、干枯,部分主干形成层褐变、植株死亡 |
1)寒冻害等级“0~4”分别表示:无、轻度、中度、重度、严重寒冻害。 |
福建冷空气过程主要发生在12月至翌年2月。根据文献[21]~[23],将福建省寒冻害低温类型分为:辐射型、平流型和混合型。从天气过程来看:(1) 辐射型寒冻害发生于冷风过境后冷高压控制下,白天天气晴朗,夜晚辐射导致大幅度降温。辐射降温期间,静风,日均温不会很低,但气温日较差大。(2) 平流型寒冻害则是强冷平流直接导致的巨幅降温,期间天气阴冷或伴有降水,日均温低,但气温日较差较辐射降温小。(3) 混合型寒冻害则是平流降温和辐射降温接连发生导致的,先经历低温阴雨天气,然后天气急转晴[24-25],绝大多数是冷平流过后的晴夜辐射降温引起[26]。鉴于此,2015—2016年冬季台湾青枣寒冻害等级划分见表 2。
测定点 | 低温过程日期 | 过程最低气温/℃ | 过程最低气温日期 | 寒冻害等级1) | 低温类型 |
福鼎县气象局 | 12月17日—19日 | 0.2 | 12月18日 | 3 | 辐射型 |
罗源县气象局 | 1月2日 | 1.5 | 1月2日 | 2 | 辐射型 |
南靖船厂棚外 | 12月16日—19日 | 2.1 | 12月18日 | 2 | 辐射型 |
南靖高港 | 12月16日—19日 | -1.9 | 12月18日 | 3 | 辐射型 |
永泰县气象局 | 12月17日—19日 | 2.4 | 12月18日 | 2 | 辐射型 |
长泰县气象局 | 12月17日 | 4.8 | 12月17日 | 1 | 辐射型 |
长泰薛厂长青枣园 | 2月6日—7日 | 2.6 | 2月7日 | 2 | 辐射型 |
福鼎县气象局 | 1月20日—27日 | -4.2 | 1月25日 | 4 | 混合型 |
南安农科所 | 1月23日—27日 | -2.4 | 1月25日 | 3 | 混合型 |
南靖船厂棚内 | 1月22日—27日 | -2.3 | 1月25日 | 3 | 混合型 |
南靖船厂棚外 | 1月22日—27日 | -3.1 | 1月25日 | 4 | 混合型 |
南靖高港村 | 1月22日—26日 | -6.8 | 1月25日 | 4 | 混合型 |
南靖五川村 | 1月24日—25日 | -0.2 | 1月25日 | 3 | 混合型 |
南靖县气象局 | 1月22日—26日 | -0.7 | 1月25日 | 3 | 混合型 |
永泰县气象局 | 12月18日 | 1.8 | 12月18日 | 2 | 混合型 |
永泰县气象局 | 1月23日—26日 | -2.2 | 1月25日 | 3 | 混合型 |
云霄县路边 | 1月22日—26日 | 1.7 | 1月25日 | 2 | 混合型 |
长泰县气象局 | 1月23日—26日 | 0.0 | 1月25日 | 3 | 混合型 |
南安农科所 | 12月22日—23日 | 3.6 | 12月22日 | 1 | 平流型 |
南靖高港 | 12月5日 | 4.1 | 12月5日 | 1 | 平流型 |
南靖高港 | 12月22日 | -2.6 | 12月22日 | 3 | 平流型 |
龙海市气象局 | 12月22日 | 6.6 | 12月22日 | 1 | 平流型 |
长泰县气象局 | 12月30日 | 4.8 | 12月30日 | 1 | 平流型 |
1)寒冻害等级“0~4”分别表示:无、轻度、中度、重度、严重寒冻害。 |
根据寒冻害灾情样本的等级划分及低温类型分类结果(表 2),对最低气温和寒冻害等级分别用直线、对数、幂函数、指数函数方程进行拟合,以一元一次直线方程的拟合效果最佳。对辐射型、平流型、混合型和无分型拟合结果如下。
$辐射型:T_d=-2.86x+7.80(n=7, r=0.942>r_{\rm 0.01})$ | (1) |
$平流型:T_d=-3.69x+8.46(n=5, r=0.945>r_{\rm 0.01})$ | (2) |
$混合型:T_d=-3.60x+9.29(n=11, r=0.915>r_{\rm 0.01})$ | (3) |
$无分型:T_d=-3.28x+8.34(n=23, r=0.953>r_{\rm 0.01})$ | (4) |
式中,Td为最低气温,x为寒冻害等级,n为样本个数,r为相关系数。由上可知,基于低温分型的青枣寒冻害等级回归模型与无分型有一定差异。低温分类与未分类的寒冻害等级点聚图见图 1、2。由图 1、2可知,不同类型Td与寒冻害等级间均有显著的负相关关系,即最低温度越低寒冻害等级越高。辐射型、平流型、混合型和无分型的一项式相关系数依次为:0.942、0.945、0.915、0.953,相关性均较好。
2.2 3种低温类型寒冻害指标的特征分析将寒冻害等级x=0、1、2、3、4分别代入方程(1) ~(4),即可得出不同寒冻害无、轻度、中度、重度和严重等级临界低温阈值,并据此初步确定各等级的指标区间(表 3)。在无分型下,当日最低气温 < 8.34 ℃时,青枣出现寒冻害;平流型和混合型分别在 < 8.46 ℃、< 9.29 ℃呈现寒冻害,较无分型高;辐射型在 < 7.8 ℃时呈现寒冻害,临界温度较无分型低。出现寒冻害的临界温度阈值依次为:辐射型≤无分型≤平流型≤混合型。在寒冻害1级区间,同一低温强度下平流型和辐射型寒冻害程度较接近,均较混合型寒冻害轻;在寒冻害2级区间,同一低温强度下平流型寒冻害程度较混合型轻,辐射型寒冻害程度介于二者之间;在3、4级区间,同一低温强度下辐射型、混合型、平流型青枣寒冻害程度依次减轻。
类型 | 临界低温阈值/℃ | ||||
0级 | 1级 | 2级 | 3级 | 4级 | |
辐射型 | Td≥7.80 | 4.94≤Td < 7.80 | 2.08≤Td < 4.94 | -0.78≤Td < 2.08 | Td < -0.78 |
平流型 | Td≥8.46 | 4.78≤Td < 8.46 | 1.09≤Td < 4.78 | -2.60≤Td < 1.09 | Td < -2.60 |
混合型 | Td≥9.29 | 5.70≤Td < 9.29 | 2.09≤Td < 5.70 | -1.51≤Td < 2.09 | Td < -1.51 |
无分型 | Td≥8.34 | 5.06≤Td < 8.34 | 1.78≤Td < 5.06 | -1.50≤Td < 1.78 | Td < -1.50 |
2015—2016年冬季,漳州市平和县遭受了3次低温过程(12月18日、1月22日—26日和2月13日),过程降温幅度大,极端最低气温较低,导致台湾青枣遭受寒冻害。平和县山格镇和平和县气象局2个台湾青枣种植区共观测到4个寒冻害灾情样本数据,将初步得到的低温分型寒冻害等级指标与4个台湾青枣寒冻害灾情样本进行对比(表 4)。由表 4可见,用低温分型获得的初步寒冻害指标与实际寒冻害等级比较,有1个样本等级差1级外,其余3个样本等级均一致,占75%;用无低温分型获得的初步寒冻害指标与实际寒冻害等级进行比较,有2个样本等级差1级,2个样本等级一致,准确率为50%。综上所述,2种低温寒冻害指标与实际寒冻害等级比较,准确率均≥50%;低温分型寒冻害指标较无低温分型准确率高25%。
地点 | 低温过程日期 | 过程最低气温/℃ | 过程最低气温日期 | 低温类型 | 实际寒冻害等级 | 低温分型判断寒冻害等级 | 无分型判断寒冻害等级 | 低温分型验证结果 | 无分型验证结果 |
平和县气象局 | 1月22日—26日 | -1.4 | 01-25 | 混合型 | 3 | 3 | 3 | 一致 | 一致 |
平和县山格镇 | 12月18日 | 0.6 | 12-18 | 辐射型 | 3 | 3 | 3 | 一致 | 一致 |
平和县山格镇 | 1月22日—26日 | -3.0 | 01-25 | 混合型 | 3 | 4 | 4 | 差1级 | 差1级 |
平和县山格镇 | 2月13日 | 1.4 | 02-13 | 平流型 | 2 | 2 | 3 | 一致 | 差1级 |
台湾青枣寒冻害是否发生及其受害程度与气象条件及果树树龄、长势、生育期有关,但最主要取决于最低气温,即冷空气过程最低气温[27]。以往研究大多笼统利用最低气温单一指标作为阈值,本研究将福建省低温类型划分为:辐射型、平流型、混合型,分别探讨3种低温类型下台湾青枣的寒冻害特征,并建立寒冻害指标。结果表明,不同低温类型下最低气温与寒冻害等级间均有显著的负相关关系。对比分析显示,无分型下青枣的寒冻害临界温度较平流型和混合型低,较辐射型高。不同等级区间,不同低温类型下寒冻害程度有较大差异:(1) 寒冻害1级区间,在同一寒冻害等级下平流型和辐射型寒冻害致灾温度较低,混合型寒冻害致灾温度较高;(2) 寒冻害2级区间,在同一寒冻害等级下混合型寒冻害致灾温度较平流型高,辐射型介于二者之间;(3) 寒冻害3、4级区间,在同一寒冻害等级下寒冻害致灾温度依次为:辐射型>混合型>平流型。个例对比验证结果显示,低温分型寒冻害等级指标准确率≥75%,无分型寒冻害等级指标准确率≥50%,即低温分型寒冻害等级指标准确度优于无分型。
本研究对福建省冬季低温天气过程进行分类,分析台湾青枣在不同低温类型下寒冻害指标特征,通过个例对比印证有分型和无分型寒冻害指标准确率,以期为台湾青枣寒冻害评估及防御提供参考。气象灾情资料的丰富度和植物学、人为等因素亦会影响研究结果,有待今后进一步完善。
[1] | 刘代兴. 台湾青枣品种、丰产栽培及云南种植区划[J]. 热带农业科技, 2003, 26(2): 6–10, 26. |
[2] | 梁开明, 曹洪麟, 徐志防, 等. 台湾青枣及野生种的光合作用日变化及光响应特征[J]. 园艺学报, 2008, 35(6): 793–798. |
[3] | 罗华建, 罗诗, 赖永超, 等. 台湾青枣果实生长发育初探[J]. 果树学报, 2002, 19(6): 436–438. |
[4] | 陈家金, 王加义, 黄川容, 等. 福建省引种台湾青枣的寒冻害风险分析与区划[J]. 中国生态农业学报, 2013, 21(12): 1537–1544. |
[5] | 张倩, 赵艳霞, 王春乙. 我国主要农业气象灾害指标研究进展[J]. 自然灾害学报, 2010, 19(6): 40–54. |
[6] | DING C K, WANG C E, GROSS K C, et al. Jasmonate and salicylate induce the expression of pathogenesis-related-protein genes and increase resistance to chilling injury in tomato fruit[J]. Planta, 2002, 214(6): 895–901. DOI: 10.1007/s00425-001-0698-9 |
[7] | 柏秦凤. 华南寒害致灾气候因子及综合指数研究[D]. 北京: 中国气象科学研究院, 2008: 4-5. http://cdmd.cnki.com.cn/Article/CDMD-85101-2008125901.htm |
[8] | 李娜. 华南寒害气候风险区划技术研究[D]. 北京: 中国气象科学研究院, 2010: 3-5. http://www.xueweilunwen.com/doc/116512 |
[9] | 卢新坤, 林旗华, 庄文彬. 台湾青枣"蜜枣"在漳州的引种试验初报[J]. 中国南方果树, 2012, 41(2): 95–97. |
[10] | 陈惠, 王加义, 杨凯. 南亚热带主要果树寒冻害研究进展[M]. 北京: 气象出版社, 2015. |
[11] | 陈惠, 王加义, 潘卫华, 等. 南亚热带主要果树冻(寒)害低温指标的确定[J]. 中国农业气象, 2012, 33(1): 148–155. |
[12] | 杨凯, 林晶, 陈惠, 等. 福建枇杷低温害临界温度和综合气候指标[J]. 中国农业气象, 2013, 34(3): 468–473. |
[13] | 王加义, 陈家金, 李丽纯, 等. GIS在福建枇杷低温冻害分析中的应用[J]. 中国农业气象, 2011, 32(S1): 153–156. |
[14] | 杨凯, 李丽纯, 陈福梓, 等. 福建省引种台湾青枣的寒冻害等级指标初探[C]//中国气象学会. 第33届中国气象学会年会. 西安: 中国气象学会, 2016: 5. |
[15] | 吴振海, 徐宗焕, 杨晓春, 等. 番木瓜寒(冻)害评估方法探讨[J]. 福建农业科技, 2007(5): 66–68. |
[16] | 刘世业, 谭宗琨, 包辉昌. 冬季低温对广西香蕉生产及产量损失影响评估初步研究[J]. 气象研究与应用, 2012, 33(3): 38–42. |
[17] | 谭宗琨, 刘世业, 唐志鹏, 等. 香蕉寒冻害等级指标及灾损指标的初步研究[J]. 自然灾害学报, 2013, 22(4): 182–192. |
[18] | 陈国保, 夏小曼, 李永平. 两种类型的低温对冬季免耕马铃薯的影响[J]. 气象研究与应用, 2010, 31(S2): 225–237, 238. |
[19] | 霍治国, 杜尧东, 姜燕, 等. 香蕉、荔枝寒害等级: QX/T 80-2007[S]. 北京: 中国标准出版社, 2007. |
[20] | 杨凯, 陈惠, 李丽纯, 等. 引种台湾青枣的寒冻害等级指标研究[J]. 自然灾害学报, 2017, 26(4): 91–97. |
[21] | 康志明. 基于TIGGE数据的我国寒潮强冷空气自动识别预报方法研究[D]. 兰州: 兰州大学, 2008: 11-22. http://www.xueweilunwen.com/doc/230535 |
[22] | 周琳, 孙照渤. 1961-2010年我国冷空气的活动特征[J]. 大气科学学报, 2015, 38(3): 342–353. |
[23] | 鹿世瑾, 王岩. 福建气候[M]. 2版. 北京: 气象出版社, 2012. |
[24] | 谭宗琨, 黄城华, 孟翠丽, 等. 甘蔗寒冻害等级指标及灾损指标的初步研究[J]. 中国农学通报, 2014, 30(28): 169–181. |
[25] | 陈瑶, 谭志坚, 樊佳庆, 等. 橡胶树寒害气象等级研究[J]. 热带农业科技, 2013, 36(2): 7–11. |
[26] | 陈惠, 徐宗焕, 潘卫华, 等. 地形闭塞的山坡下部冬季气温特征分析[J]. 中国农业气象, 2010, 31(2): 300–304. |
[27] | 陈惠, 杨凯, 李丽纯, 等. 莲雾寒冻害低温等级指标的确定[J]. 生态学杂志, 2017, 36(4): 1047–1052. |