全 文 :中国生态农业学报 2013年 12月 第 21卷 第 12期
Chinese Journal of Eco-Agriculture, Dec. 2013, 21(12): 1537−1544
* 中国气象局气象关键技术集成与应用项目(CAMGJ2012M27)、福建省自然科学基金项目(2012J01160)、中国气象局小型业务项目
(2012209)和福建省气象局基金项目(2011K03)资助
陈家金(1964—), 男, 本科, 正研级高级工程师, 主要从事农业气象研究。E-mail: cjj8284@163.com
收稿日期: 2013−04−24 接受日期: 2013−09−05
DOI: 10.3724/SP.J.1011.2013.30403
福建省引种台湾青枣的寒冻害风险分析与区划*
陈家金 1 王加义 1 黄川容 1 郑东旗 2 林 晶 1 李丽纯 1 马治国 1
(1. 福建省气象科学研究所 福州 350001; 2. 福州市气象局 福州 350014)
摘 要 为了摸清台湾青枣在福建复杂地形下引种的寒冻害风险, 规避在高风险区引种, 本文依据灾害风险
系统理论、气象和农业经济资料, 在分析寒冻害危险性、青枣脆弱性和引种区防寒防冻能力的基础上, 应用层
次分析法、熵权系数法、加权综合评价法等构建了台湾青枣寒冻害风险评价指标和模型, 对福建省引种台湾
青枣的寒冻害风险程度进行等级评价和区域划分。结果表明: 在风险因素构成上, 寒冻害危险性在综合风险中
起决定性作用, 引种区防寒防冻能力只能略微减轻寒冻害带来的风险。惠安县、厦门市以南沿海县市的青枣
引种区为轻度寒害风险区域, 该区域寒冻害致灾危险性指数低, 中度以上的寒冻害发生频次少或没有出现过;
漳州市内陆的较低海拔地区, 福州市沿海低海拔地区, 同安区、惠安县及永春县的部分区域为中度寒冻害风险
区域, 该区域寒冻害危险性指数相对较低, 中度以上寒冻害发生频次相对较少; 其余青枣引种区的冻害风险
在重度以上, 其中龙岩市、漳州市西部山区、福州市西部和北部山区、永春县、宁德市辖区及霞浦县有严重
冻害风险。结合近年来福建出现的青枣寒冻害个例分析, 青枣寒冻害形态表征程度与该地域风险评估结果较
为吻合, 风险区划结果可为福建引种台湾青枣提供布局决策参考。
关键词 青枣 引种 福建省 台湾省 寒冻害 风险分析 风险区划
中图分类号: S166 文献标识码: A 文章编号: 1671-3990(2013)12-1537-08
Risk analysis and regionalization of cold and freezing damage to
Taiwan-based Zizyphus mauritiana in Fujian Province
CHEN Jia-Jin1, WANG Jia-Yi1, HUANG Chuan-Rong1, ZHENG Dong-Qi2,
LIN Jing1, LI Li-Chun1, MA Zhi-Guo1
(1. Fujian Institute of Meteorological Sciences, Fuzhou 350001, China; 2. Fuzhou Meteorology Bureau, Fuzhou 350014, China)
Abstract This study determined the risk of cold and freezing damage to Zizyphus mauritiana introduced to the complex terrains of
Fujian Province from Taiwan Province in order to avoid planting the crop in high risk areas. Based on disaster risk system theory, the
study used meteorological, agricultural and economic data to analyze potential hazards of cold and freezing damage and vulnerability
of Z. mauritiana. The study also investigated the prevention measures of cold and freezing damage in the cultivated regions and
developed a model-driven index system of evaluation based on the AHP, entropy-weight coefficient and comprehensive weighting
methods. The study further determined the degree and region of risk of cold and freezing damage to Z. mauritiana in Fujian. In terms
of risk factor composition, the results showed that potential hazard was the determinant factor of cold and freezing damage and
prevention measures only mitigated this potential. For risk region, the areas of lighter risk of cold and freezing damage to Z.
mauritiana were Hui’an County of Xiamen Region and the coastal counties to the south of Xiamen Region. The potential hazard of
cold and freezing damage was low in these regions and frequency occurrence of moderate or severe damages were less or even
non-existent. The regions with moderate risk were in the low elevation areas of Zhangzhou Inland Region, low elevation areas of the
coastal counties of Fuzhou Region, and parts of Yongchun County, Tong’an District and Hui’an County. The occurrence frequency of
moderate or severe damages was relatively low in the areas. Other cultivated regions of Z. mauritiana had severe risk of cold and
freezing damage. Among the regions, cold and freezing damage risk was severe in Longyan City, the western mountain areas of
Zhangzhou Region, the western and northern mountain areas of Fuzhou Region, Yongchun County, Ningde City and Xiapu County.
1538 中国生态农业学报 2013 第 21卷
Based on the risk evaluation results and Z. mauritiana cold and freezing damage, the morphological characterization of cold and
freezing damage to Z. mauritiana was consistent with the risk evaluation results. The results of the study offered the basis for
decisions on Z. mauritiana cultivation and arrangement in Fujian Province.
Key words Zizyphus mauritiana, Introduce, Fujian Province, Taiwan Province, Cold and freezing damage, Risk analysis,
Risk regionalization
(Received Apr. 24, 2013; accepted Sep. 5, 2013)
台湾青枣 , 学名为毛叶枣(Zizyphsu mauritiana
Lam.), 属热带、亚热带常绿小乔木果树。台湾青枣
适应性强, 当年种植当年挂果, 且果实成熟期正值
冬春季节, 是水果淡季, 种植效益高, 因此, 近年来
福建将其作为短、丰、快的新兴果树种类加以引种
推广。但也存在不合理布局导致的寒冻害问题, 导
致植株受冻、果实受害并影响产量[1–3]。因此, 开展
台湾青枣在福建引种的灾害风险区划, 对实现该区
域台湾青枣的安全引种和可持续发展至关重要。
果树气象灾害风险区划的目的正是将种植灾害
风险控制在引种前, 找到无果树气象灾害风险或风
险低的区域。目前国内外研究者针对热带果树的寒
冻害研究, 更多地停留在寒冻害指标以及影响分析
上, 而风险区划方面的研究则较少。霍治国等[4−8]针
对华南地区香蕉 (Musa nana Lour.)、荔枝 (Litchi
chinensis Sonn.)开展过寒冻害风险分析与区划研究。
针对热带果树台湾青枣的研究大多侧重于栽培技术
的介绍, 还未见有关青枣气象灾害风险区划方面的
研究。因此, 亟需针对台湾青枣引种过程中最主要
的农业气象灾害开展风险区划, 在福建省复杂地形
下甄别出青枣气象灾害低风险的精细分布区域, 把
灾害风险规避在种植前。本文针对福建省果业结构调
整的需要, 以福建省引种台湾青枣的区域为背景, 针
对台湾青枣引种过程中存在的最主要农业气象灾害
即寒冻害风险, 采用“多指标综合风险评估法”, 构建
风险区划指标体系和评估模型, 对台湾青枣在福建
省引种的寒冻害风险进行分析和区划, 以期为台湾
青枣在福建省的合理布局、安全生产提供决策依据。
1 材料与方法
1.1 研究区概况
福建省地处我国东南沿海 , 位居北纬 23°33′~
28°19′, 东经 115°50′~120°43′之间, 跨中、南亚热带,
境内地形复杂, 多山、丘陵起伏, 地形气候多样。福
建省各县市的年平均气温为 15.0~21.7 ℃; 平均稳
定通过 0 ℃的积温为 5 473.6~7 948.8 ℃, 除鹫峰山
区的屏南县、周宁县、寿宁县以及长汀县以北至浦
城县的西北部外, 其他县市的积温在 6 500 ℃以上;
年极端最低气温为−12.8~3.8 ℃ , 并呈现由东南向
西北递减的趋势 , 极端最低气温≤ 0 ℃县市占
92%。近年来, 福建不断推进农业结构调整, 出台果
业发展规划, 以促进果业调整, 尤其是南亚热带区域
从台湾引种了大量的热带果树, 一定程度优化了果
业结构并取得了成效。2011 年福建省引种的台湾青
枣面积已达 2 109 hm2, 主要分布在漳州市, 占全省
总种植面积的 96.2%。但是在亚热带区域内引种热带
果树, 热量条件成为主要限制因子, 研究区域内的寒
冻害成为台湾青枣引种和扩种最主要的制约要素。
1.2 数据来源
青枣引种区各县 1971—2011年气温、降水、风
速的资料来源于国家气象信息中心和福建省气象局
的气象观测数据; 农业统计资料来源于历年福建省
农村统计年鉴; 地理信息资料来源于“数字福建”提
供的 1︰250 000基础地理背景资料, 包括 DEM、行
政边界等数据。
1.3 区划方法
采用“多指标综合风险评估法”对台湾青枣在福
建省引种的寒冻害综合风险进行区划和评估。其主
要步骤[9]是: (1)构建指标体系; (2)对诸指标进行归一
化处理, 消除量纲差异; (3)对诸指标用主、客观赋权
等手段综合确定权重; (4)计算评估单元的风险指数;
(5)根据权重构建评估模型; (6)对计算出的区域风险
指数值进行指标分级; (7)应用GIS技术制作风险区划
图; (8)评估不同区域的风险大小。限于篇幅, 本文重
点对区划方法中指标体系构建、指标权重确定、风险
指数计算和风险区划模型构建的关键技术进行阐述。
1.3.1 指标体系构建方法
基于“自然灾害风险指数法”的思路构建风险区
划指标体系。灾害风险受致灾因子危险性、承险体
脆弱性和防灾减灾能力 3 个因子共同影响, 是这 3
个因子整体作用的结果[9–10]。通常采用自然灾害风
险指数表征风险程度, 可表示为: 自然灾害风险指
数=f(危险性、脆弱性、防灾减灾能力)。
通过分析台湾青枣的生物学特性以及在福建省
引种的寒冻害情况, 确定台湾青枣寒冻害致灾危险
性区划指标; 分析台湾青枣自身的脆弱性, 即种植
的暴露性和对寒冻害的敏感性, 确定青枣脆弱性指
标; 分析台湾青枣引种区防灾减灾能力情况, 确定
引种区防灾减灾能力大小的指标; 构建出由寒冻害
致灾因子危险性、青枣脆弱性以及引种区防冻防寒
能力这 3方面组成的风险区划指标体系(图 1)。
第 12期 陈家金等: 福建省引种台湾青枣的寒冻害风险分析与区划 1539
图 1 台湾青枣引种寒冻害风险区划指标体系
Fig. 1 Index system of the risk regionalization of cold and freezing damage for Zizyphus mauritiana
图 1中 C1~C4为极端最低气温, C5为种植面积占全省种植面积比, C6为歉年平均减产率, C7为歉年减产率变异系数, C8为减产率发生概
率, C9为地膜覆盖面积占农作物播种面积比, C10为有效灌溉面积占耕地面积比, C11为农民人均纯收入。C1−C4: extreme minimum temperature;
C5: percentage of planting area to total planting area; C6: average yield-decreasing rate during disaster years; C7: coefficient of variance of
yield-decreasing rate during disaster years; C8: occurrence probability of yield-decreasing rate; C9: percentage of films mulching area to planting area;
C10: percentage of effective irrigation area to cultivated area; C11: net income per farmer.
1.3.2 指标权重确定方法
由于台湾青枣引种寒冻害风险指标体系中每个
指标对最终风险构成的贡献率不同, 因此, 必须采
取合适的指标赋权方法来反映每个指标的风险重要
性。本文综合应用层次分析法(AHP)、熵权系数法确
定风险区划指标权重。
层次分析法是一种对指标进行定性定量分析的
方法[11–13]。主要步骤: (1)构建风险区划指标体系框
架; (2)邀请不同领域的专家, 按照1~9比例标度法进
行打分, 对所选指标进行重要性比较; (3)构建判断
矩阵; (4)层次单排序及一致性检验; (5)层次总排序
及组合一致性检验; (6)确定各指标的权重(ωi)。
熵权系数法是一种基于分析各评价指标特征值
之间变异程度的客观赋权方法[14–15]。主要步骤: (1)
构建区划指标特征值矩阵; (2)对特征值进行标准化处
理, 使其值处于0和1之间, 以消除量纲的影响; (3)计
算区划指标特征值比重; (4)计算区划指标的熵权(ai)。
将上述2种方法确定的主观权重和客观权重相
融合得出主客观综合权重值, 作为风险评估指标权
重值(σi)。
1
/
m
i i i i i
i
a aσ ω ω
=
= ∑ (1)
式中, i为评价指标, m为评价指标个数, ai为第 i个指
标熵权, ωi为第 i个指标层次分析权重。
1.3.3 风险指数计算方法
采用“加权综合评价法(WCA)”计算风险指数[9–10]。
该方法是依据评价指标对评价总目标影响的重要程
度, 预先分配一个响应的权重系数, 然后与相应的
被评价对象各指标的特征量化值相乘后再相加。计
算公式为:
1
m
j i ij
i
pλ σ
=
= ∑ (2)
式中, λj为某评估对象(单元)计算所得的风险指数值,
σi为第 i项指标的综合权重系数(σi>0,
1
1
m
i
i
σ
=
=∑ ), pij
为每个指标 i特征值在评估单元 j所占比重(0≤pij≤1),
i为评价指标, j为评估单元, m为评价指标个数。
1.3.4 寒冻害风险模型构建方法
地理模型构建。由于青枣寒冻害危险性指标评
估指数是基于评估单元(县)地面气象观测站单点观
测值计算所得, 属点状数据, 未能反映评估单元(县)
内其他点的风险情况。因此, 通过建立寒冻害危险
性指数与测站经度、纬度和海拔的多元回归模型 ,
利用GIS技术对其进行地理插值推算 , 将点状数据
转化为面状格网数据, 解决无气象观测资料地区的
风险指数计算问题, 实现寒冻害危险性评估指数的
空间连续分布; 对与地理因子无关的青枣脆弱性、
引种区防寒防冻能力的评估指标采用基本单元值计
算, 以面上数据代表评估单元内的任一点数据。
区划模型构建。按照“加权综合评价法”, 对风险
区划指标体系中各指标构成因子权重进行加权, 得
出青枣寒冻害危险性、青枣脆弱性、种植区防寒防
冻能力和综合风险指数的区划模型。
2 结果与分析
2.1 台湾青枣的寒冻害风险分析
2.1.1 寒冻害危险性分析
台湾青枣属于热带水果, 引种至福建亚热带区
域, 其主要限制因子是温度; 冬季极端低温决定引
种的台湾青枣露地栽培能否越冬成活, 同时冬季恰
是青枣果实发育至成熟期, 寒冻害还会造成青枣果
实落果或烂果。据研究, 台湾青枣在极端最低气温
达到 8 ℃以下超过 5 d, 幼果会受寒害; 气温下降至
3 ℃时, 叶片枯死; −1~−2 ℃时枝条、花果冻死; −4 ℃
1540 中国生态农业学报 2013 第 21卷
时成年枣树地上部分被冻死[16–21]。虽然台湾青枣树
体本身较为耐寒, 但由于以果越冬且果实耐寒性较
树体差 , 因此 , 在分析台湾青枣寒冻害危险性时 ,
着重考虑青枣以果越冬带来的果实寒冻害危险, 同
时兼顾树体本身的寒冻害危险, 确定以极端最低气
温(Td)作为寒冻害危险性表征指标, 并对其进行分
级, 设轻度、中度、重度、严重 4 级寒冻害危险性
程度。表 1 列出了青枣不同程度寒冻害危险性表征
指标及权重, 可见极端最低气温越低, 青枣的致灾
危险性越大。
表 1 福建省引种台湾青枣的寒冻害危险性表征指标及权重
Table 1 Index and weight of potential hazards of cold and freezing damage of Zizyphus mauritiana introduced from Taiwan
Province to Fujian Province
危险性等级 Potential hazards grade 指标 Index 指标权重 Weight of the index
轻度 Light 极端最低气温(1.0≤Td<5.0) Extreme minimum temperature ( )℃ 0.010 2
中度 Moderate 极端最低气温(−1.0≤Td<1.0) Extreme minimum temperature (℃) 0.054 7
重度 Severe 极端最低气温(−3.0≤Td<−1.0) Extreme minimum temperature (℃) 0.163 9
严重 Very severe 极端最低气温(Td<−3.0) Extreme minimum temperature (℃) 0.771 2
图 2a给出了青枣寒冻害危险性区划结果。结果
表明泉州市以南沿海县市台湾青枣引种区为轻度寒
害危险性区域; 漳州市内陆(长泰县、南靖县、平和
县、云霄县)的低海拔地区、福州市沿海低海拔地区、
同安区及惠安县的部分区域为中度寒冻害危险性区
域; 其余引种区(龙岩市、漳州市西部高海拔区域、
福州市高海拔区域、宁德市)的冻害危险性在重度以
上, 其中龙岩市、漳州市西部山区、福州市西部和
北部(闽侯县及连江县的山区)、永春县、宁德市辖区
及霞浦县有严重冻害危险性。
2.1.2 青枣脆弱性分析
对于台湾青枣这个承险体而言, 采用最小评估
单元(县)的青枣种植面积占全省总种植面积来反映
青枣暴露性, 种植面积越大, 其所承受的灾害风险
越大; 此外, 青枣果实寒冻害是主要农业气象灾害,
常导致产量损失, 因此, 歉年减产率相关指标可一
定程度反映台湾青枣寒冻害灾损的敏感性。表 3 列
出了青枣脆弱性表征指标及权重, 其重要性从大到
小的顺序是: 种植面积比例>减产发生概率>歉年平
均减产率>歉年减产率变异系数。
图 2b给出了青枣脆弱性区划结果。结果表明长
泰县和云霄县的脆弱性低, 主要源于这两县台湾青
枣歉年平均减产率和歉年发生概率均比较低, 此外
云霄县青枣种植面积小; 高脆弱性区域分布在漳浦
县、南靖县、平和县, 主要源于漳浦县、南靖县和
平和县的台湾青枣种植面积大, 歉年平均减产率、
歉年变异系数和发生概率均比较高的缘故; 其余引
种区属中度脆弱性, 该区域种植面积也不大, 歉年
平均减产率、歉年变异系数和歉年发生概率相对属
中等。
2.1.3 引种区防寒防冻能力分析
对于台湾青枣引种区而言, 为防止或减轻引种
区寒冻害对台湾青枣生长发育的影响, 在基础防灾
减灾能力方面, 主要考虑对防御寒冻害的资金投入
能力; 而在专项应灾能力方面, 主要考虑预防寒冻
害而配备的灌溉工程及地膜覆盖情况。表 3 列出了
青枣引种区防寒防冻能力的表征指标及权重, 可见
地膜覆盖对防寒防冻起主要作用。
图 2c给出了引种区防寒防冻能力区划结果。结
果表明厦门市的防寒防冻能力高, 主要源于该市地
膜覆盖比例高, 同时农民人均纯收入也高; 霞浦县、
福清市、惠安县、漳浦县、云霄县和东山县的防寒
防冻能力低, 主要是源于地膜覆盖比例不高, 部分
县市农民人均纯收入较低, 其余引种区的防寒防冻
能力属中等。
2.2 寒冻害综合风险区划
2.2.1 区划指标和等级划分
台湾青枣寒冻害综合风险区划指标包括寒冻害
危险性指标、青枣脆弱性指标和引种区防寒防冻能
力指标。表 2 列出了各指标的权重计算结果, 可以
看出寒冻害危险性在综合风险中起决定性作用。
在确定各风险指标的权重后, 根据式(2)计算各
个评估单元的风险评估指数(λj)。表 3列出了各评估
单元台湾青枣引种寒冻害风险指标归一化评估指数
的计算结果。
从表 3 看出, 台湾青枣引种区评估单元气象观
测点寒冻害致灾因子危险性归一化指数为 0~1.000 0,
东山县的寒冻害危险指数为 0; 福清市以南沿海县市
的寒冻害危险指数在 0.040 0 以下; 指数为 0.040 0~
0.1000 0 的县市有闽侯县、宁德市辖区、南靖县、
平和县; 指数为 0.1~0.5的县市有霞浦县、连江县、
上杭县、龙岩市辖区、永春县; 大于 0.500 0的县市
有永定县和武平县, 以武平县为最大。
第 12期 陈家金等: 福建省引种台湾青枣的寒冻害风险分析与区划 1541
图 2 福建省引种台湾青枣的寒冻害危险性(a)、脆弱性(b)、防寒防冻能力(c)及综合风险(d)区划图
Fig. 2 Regionalization maps of potential hazards (a), vulnerability (b), prevention capacity (c), integrated risk (d) of cold and
freezing damage of Zizyphus mauritiana introduced from Taiwan Province to Fujian Province
青枣脆弱性归一化指数为 0.296 4~1.000 0。脆
弱性风险指数低于 0.400 0 的县市有长泰县和云霄
县, 以云霄县(0.296 4)最低; 指数大于 0.800 0的县
市有南靖县、漳浦县和平和县, 以漳浦县为最大; 其
余种植区的脆弱性风险指数为 0.400 0~0.800 0。
青枣引种区防寒防冻能力归一化指数为 0.180 0~
1.000 0, 其中指数小于 0.300 0的县市有霞浦县、福
清市、惠安县、漳浦县、东山县、云霄县; 指数大
于 0.700 0的县市有厦门市; 其余引种区的防寒防冻
能力归一化指数为 0.300 0~0.700 0。
青枣引种区评估单元气象观测点各综合风险归
一化指数为 0.020 0~1.000 0, 其中指数小于 0.100 0
的县市有长泰县、龙海市、漳浦县、厦门市、惠安
县、诏安县、东山县; 指数为 0.100 0~0.200 0的县
市有闽侯县、宁德市辖区、连江县、永春县、南靖
县和平和县; 指数为 0.200 0~0.500 0 的县市有霞浦
县、连江县、上杭县、龙岩市辖区; 指数大于 0.500 0
的县市有永定县和武平县, 以武平县为最大。
采用自然断点法(属不等值分级方法)对风险指数
进行分级, 即通过评估区域内所有数据值与平均值之
差的原则来寻找特征点, 结合实地风险大小的调查,
得出台湾青枣寒冻害综合风险等级划分标准(表4)。
1542 中国生态农业学报 2013 第 21卷
表 2 福建省引种台湾青枣的脆弱性、防寒防冻能力及综合风险表征指标及权重
Table 2 Index and weight of vulnerability, prevention capacity and integrated risk of cold and freezing damage of Zizyphus mau-
ritiana introduced from Taiwan Province to Fujian Province
项目 Item 指标 Index 指标权重 Index weight
种植面积占全省种植面积比 Percentage of planting area to total planting area 0.386 5
歉年平均减产率 Average yield-decreasing rate during disaster years 0.153 3
歉年减产率变异系数
Coefficient of variance of yield-decreasing rate during disaster years
0.097 1
青枣脆弱性
Vulnerability of Z. mauritiana
减产率发生概率 Occurrence probability of yield-decreasing rate 0.363 1
地膜覆盖面积占农作物播种面积比
Percentage of films mulching area to planting area
0.809 6
有效灌溉面积占耕地面积比
Percentage of effective irrigation area to cultivated area
0.137 4
防寒防冻能力
Prevention capacity of cold and
freezing damage
农民人均纯收入 Net income per farmer 0.053 0
寒冻害危险性 Potential hazards of cold and freezing damage 0.926 2
青枣脆弱性 Vulnerability of Z. mauritiana 0.061 1
综合风险
Integrated risk
引种区防寒防冻能力
Prevention capacity of cold and freezing damage of introduction areas
0.012 7
表 3 台湾青枣引种寒冻害归一化风险评估指数
Table 3 Normalized indexes of comprehensive risk assessment on cold and freezing damage of Zizyphus mauritiana introduction in
Fujian Province
评估单元 Evaluation unit 危险性
Hazard potential
脆弱性
Vulnerability
防寒防冻能力
Prevention capacity
综合风险
Comprehensive risk
霞浦县 Xiapu County 0.194 8 0.776 1 0.205 8 0.241 2
闽侯县 Minhou County 0.076 9 0.679 3 0.353 2 0.117 8
宁德市 Ningde City 0.081 2 0.765 9 0.313 5 0.130 3
连江县 Lianjiang County 0.189 6 0.411 0 0.374 3 0.200 9
武平县 Wuping County 1.000 0 0.722 2 0.403 9 1.000 0
上杭县 Shanghang County 0.364 0 0.643 5 0.484 1 0.386 0
龙岩市 Longyan City 0.256 8 0.614 0 0.351 4 0.283 4
永春县 Yongchun County 0.113 2 0.446 2 0.306 9 0.132 4
福清市 Fuqing City 0.020 6 0.691 1 0.298 0 0.066 2
永定县 Yongding County 0.672 2 0.713 2 0.455 2 0.686 2
长泰县 Changtai County 0.030 9 0.384 8 0.459 0 0.046 0
南靖县 Nanjing County 0.092 8 0.868 7 0.547 3 0.146 8
平和县 Pinghe County 0.081 7 0.989 7 0.591 0 0.146 2
龙海市 Longhai County 0.015 3 0.626 7 0.353 5 0.054 5
漳浦县 Zhangpu County 0.010 8 1.000 0 0.291 5 0.084 5
同安区 Tong’an District 0.020 1 0.678 3 1.000 0 0.053 2
惠安县 Hui’an County 0.006 1 0.433 3 0.213 2 0.030 7
厦门市 Xiamen City 0.005 0 0.678 3 1.000 0 0.038 9
诏安县 Zhao’an County 0.021 0 0.423 7 0.305 0 0.042 6
东山县 Dongshan County 0.000 0 0.611 6 0.183 3 0.041 3
云霄县 Yunxiao County 0.009 4 0.296 4 0.203 4 0.021 7
表 4 台湾青枣寒冻害综合风险等级划分标准
Table 4 Grading standard for comprehensive risk assessment of cold and freezing damage of Zizyphus mauritiana
风险等级 Risk grade 轻度 Mild 中度 Moderate 重度 Severe 严重 Severity
综合风险归一化指数
Normalized index of comprehensive risk assessment
0≤I<0.1 0.1≤I<0.3 0.3≤I<0.5 0.5≤I≤1
根据式(3)台湾青枣引种寒冻害综合风险指数区
划模型, 将寒冻害危险性、青枣脆弱性、引种区防
寒防冻能力 3 个图层叠加, 结合综合风险等级划分
标准, 得到青枣寒冻害综合风险分级区划图(图 2d)。
I =0.926 2Ih +0.061 1Iv−0.012 7Ic (3)
式(3)中, Ih、Iv、Ic、I分别代表寒冻害危险性、青枣
第 12期 陈家金等: 福建省引种台湾青枣的寒冻害风险分析与区划 1543
脆弱性、引种区防寒防冻能力和综合风险的归一化
评估指数。
2.2.2 综合风险区划结果
从台湾青枣引种区寒冻害综合风险区划图(图
2d)可以看出, 惠安县、厦门市以南沿海县市台湾青
枣引种区为轻度寒害风险区域, 该区域台湾青枣寒
害危险性轻, 区域内中度以上的寒冻害危险发生频
次少或没有出现过(厦门市、东山县), 虽然区域脆弱
性大部属中等, 防寒防冻能力各异, 但脆弱性和防
寒防冻能力的影响权重较小, 寒冻害危险性对综合
风险起主导作用, 因此, 寒冻害对台湾青枣影响的
综合风险属轻度, 是福建省引种台湾青枣灾害风险
最低的区域; 漳州市内陆(长泰县、南靖县、平和县、
漳浦县、云霄县)的较低海拔地区, 福州市沿海低海
拔地区, 同安区、惠安县及永春县的部分区域为中
度寒冻害风险区域, 该区域中度寒冻害发生的频次
也相对较少, 寒冻害危险性属中等, 引种区台湾青
枣的脆弱性和防寒防冻能力虽有不同, 但危险性占
决定性作用, 因此, 寒冻害对台湾青枣影响的综合
风险属中度, 是福建省引种台湾青枣灾害风险较低
的区域; 其余青枣引种区(龙岩市、漳州市西部高海
拔区域、永春县大部、福州市高海拔区域、宁德市)
的冻害危险性在重度以上, 其中龙岩市、漳州市西
部山区、福州市西部和北部(闽侯县及连江县的山
区)、永春县、宁德市辖区及霞浦县有严重冻害风险,
引种风险高。
2.2.3 风险区划结果的验证
以近年来福建省引种的台湾青枣的寒冻害个例
验证风险区划结果的可靠性。
1999年12月21日 , 永定县极端最低气温达到
−1.7 , 23℃ 日极端低温达−5.1 , ℃ 露天栽培的青枣
大部分被冻死, 但树体较大的植株于次年重新萌发
结果 [1]; 而风险区划结果显示 , 永定县青枣寒冻害
的综合风险属重度。
2005年1月1日, 霞浦县极端最低气温达−0.8 , ℃
台湾青枣遭受冻害。作者曾对该过程造成的霞浦县
青枣冻害进行调查, 调查显示该引种地的最低气温
低于霞浦县地面观测台站的气温, 地处该县的水潮
村引种的台湾青枣枝叶和果实均遭受重度以上冻害,
出现落果和烂果。而风险区划结果显示, 霞浦县青
枣寒冻害的综合风险属重度以上。
2010年12月17日, 平和县极端最低气温达−0.4 , ℃
青枣出现冻害。作者曾对该过程造成的平和县青枣
寒冻害进行调查, 调查显示平和县郊区露天栽培的
台湾青枣出现冻果, 部分脱落, 留在树上的果实长
势较差, 影响产量, 但不冻树, 枝叶无冻, 属于中度
冻害形态特征。而风险区划结果显示, 平和县青枣
寒冻害的综合风险属中度至重度, 其中调查区域(平
和县郊区)属于中度风险区域。
上述事例表明, 台湾青枣寒冻害的形态表征程
度与区划结果较为吻合, 某种程度验证了风险区划
结果的可靠性。
3 结论与讨论
(1)惠安县、厦门市以南沿海县市的青枣引种区
为轻度寒害风险区域; 漳州市内陆的较低海拔地区,
福州市沿海低海拔地区, 同安区、惠安县及永春县
的部分区域为中度寒冻害风险区域; 其余引种区的
青枣冻害风险在重度以上。近年来福建省区域出现
的青枣寒冻害实况个例与区划结果较为吻合, 一定
程度验证了风险区划结果的可靠性, 可为福建省引
种台湾青枣提供布局决策参考。
(2)在农业气象灾害风险区划方法当中, 指标体
系构建、指标权重的确定以及风险区划指标的分级
直接影响到评估结果的可靠性。本研究在考虑数据
来源的基础上, 从寒冻害危险性、青枣脆弱性和引
种区防寒防冻能力 3 个方面确定的风险表征指标能
较充分反映寒冻害风险构成, 当然还可尝试其他表
征指标, 尤其是青枣脆弱性表征指标; 而在指标权
重确定方面, 采用的融合主客观赋权方法克服了以
往一些研究单纯依赖专家主观打分带来的局限性 ,
消除了单纯依赖专家经验或基于统计指标发生概率
计算指标权重的片面性, 当然还可以尝试不同赋权
方法, 比较择优确定权重; 在区划图的不同风险区
域分级上, 采用指标自然断点法并结合对各地历年
灾害程度的判断, 划分出风险等级区划指标, 基本
可反映各地的风险大小情况, 但应随着数据时间序
列的推移以及灾害极端个例的出现, 不断进行动态
分析和评估, 以修订和完善风险区划结果。
参考文献
[1] 卢扬荣 . 台湾青枣栽培技术总结[J]. 中国南方果树 , 2001,
30(4): 62–63
Lu Y R. The summary on cultivation techniques of Zizyphus
mauritiana[J]. South China Fruits, 2001, 30(4): 62–63
[2] 郑小琴 . 闽南地区毛叶枣低温伤害分析 [J]. 福建果树 ,
2006(1): 33–34
Zheng X Q. Analysis of low temperature and freezing injury
of Zizyphus mauritiana Lam. in Southern Fujian[J]. Fujian
Fruits, 2006(1): 33–34
[3] 范光南 . 浅析寒害对龙岩市引种南亚热带果树的影响[J].
福建热作科技, 2004, 29(1): 45–46
Fan G N. Analysis on influence of cold damage to subtropical
fruit introduction in Longyan[J]. Fujian Science & Technol-
1544 中国生态农业学报 2013 第 21卷
ogy of Tropical Crops, 2004, 29(1): 45–46
[4] 霍治国, 李世奎, 王素艳, 等. 主要农业气象灾害风险评估
技术及其应用研究[J]. 自然资源学报, 2003, 18(6): 692–702
Huo Z G, Li S K, Wang S Y, et al. Study on the risk evalua-
tion technologies of main agrometeorological disasters and
their application[J]. Journal of Natural Resources, 2003, 18(6):
692–702
[5] 植石群 , 刘锦銮 , 杜尧东 , 等 . 广东省香蕉寒害风险分析
[J]. 自然灾害学报, 2003, 12(2): 113–116
Zhi S Q, Liu J L, Du Y D, et al. Risk analysis of cold damage
to banana in Guangdong Province[J]. Journal of Natural Dis-
asters, 2003, 12(2): 113–116
[6] 杜尧东 , 李春梅 , 毛慧琴 . 广东省香蕉与荔枝寒害致灾因
子和综合气候指标研究 [J]. 生态学杂志 , 2006, 25(2):
225–230
Du Y D, Li C M, Mao H Q. Disaster inducing factors and in-
tegrated climatic index for banana and litchi chilling injuries
in Guangdong Province[J]. Chinese Journal of Ecology, 2006,
25(2): 225–230
[7] 刘锦銮 , 杜尧东 , 毛慧琴 . 华南地区荔枝寒害风险分析与
区划[J]. 自然灾害学报, 2003, 12(3): 126–130
Liu J L, Du Y D, Mao H Q. Risk analysis and zonation of cold
damage to litchi in South China[J]. Journal of Natural Disas-
ters, 2003, 12(3): 126–130
[8] 李娜, 霍治国, 贺楠, 等. 华南地区香蕉、荔枝寒害的气候
风险区划[J]. 应用生态学报, 2010, 21(5): 1244–1251
Li N, Huo Z G, He N, et al. Climatic risk zoning for banana
and litchi’s chilling in jury in South China[J]. Chinese Journal
of Applied Ecology, 2010, 21(5): 1244–1251
[9] 陈家金, 王加义, 李丽纯, 等. 影响福建省龙眼产量的多灾
种综合风险评估[J]. 应用生态学报, 2012, 23(3): 819–826
Chen J J, Wang J Y, Li L C, et al. Integrated risk evaluation of
multiple disasters affecting Longyan yield in Fujian Province,
East China[J]. Chinese Journal of Applied Ecology, 2012,
23(3): 819–826
[10] 张会, 张继权, 韩俊山. 基于 GIS 技术的洪涝灾害风险评
估与区划研究——以辽河中下游地区为例[J]. 自然灾害学
报, 2005, 14(6): 141–146
Zhang H, Zhang J Q, Han J S. GIS-based assessment and
zoning of flood /waterlogging disaster risk: a case study on
middle and lower reaches of Liaohe River[J]. Journal of
Natural Disasters, 2005, 14(6): 141–146
[11] 徐建华. 现代地理学中的数学方法[M]. 北京: 高等教育出
版社, 2004: 224–250
Xu J H. Mathematical methods in contemporary geography[M].
Beijing: Higher Education Press, 2004: 224–250
[12] 文世勇, 赵冬至, 陈艳拢, 等. 基于 AHP 法的赤潮灾害风
险评估指标权重研究[J]. 灾害学, 2007, 22(2): 9–14
Wen S Y, Zhao D Z, Chen Y L, et al. Research on index
weighting of risk assessment of harmful algal bloom hazard
based on AHP[J]. Journal of Catastrophology, 2007, 22(2):
9–14
[13] 王以彭 , 李结松 , 刘立元 . 层次分析法在确定评价指标权
重系数中的应用 [J]. 第一军医大学学报 , 1999, 19(4):
377–379
Wang Y P, Li J S, Liu L Y. Application of AHP in the
determination of the weights coefficient of evaluation
index[J]. Journal of First Military Medical University, 1999,
19(4): 377–379
[14] 汤瑞凉, 王 . 农作物品种综合评判的熵权系数法研究[J].
资源开发与市场, 2002, 18(5): 3–4
Tang R L, Wang Y. Method of entropic coefficients for com-
prehensive evaluation of crop varieties[J]. Resource Devel-
opment and Market, 2002, 18(5): 3–4
[15] 付强. 数据处理方法及其农业应用[M]. 北京: 科学出版社,
2006: 466–469
Fu Q. Data processing method and its application in agricul-
ture[M]. Beijing: Science Press, 2006: 466–469
[16] 孙浩元, 续九如, 王玉柱, 等. 应用灰色聚类和灰色关联分
析对台湾青枣引种适宜区的选择[J]. 中国农学通报, 2006,
22(4): 143–146
Sun H Y, Xu J R, Wang Y Z, et al. Selection of suitable cli-
matic regions for introduction of Ziziphus mauritiana Lam.
using the methods of gray clustering and gray connecting
analysis[J]. Chinese Agricultural Science Bulletin, 2006,
22(4): 143–146
[17] 姚志文 , 黄巧敏 . 新丰地区种植台湾青枣的农业气候条件
分析[J]. 广东气象, 2004(4): 26–27, 29
Yao Z W, Huang Q M. Analysis of agricultural climatic con-
dition for Ziziphus mauritiana planting in Xinfeng region[J].
Guangdong Meteorology, 2004(4): 26–27, 29
[18] 卢新坤 , 林旗华 , 庄文彬 . 台湾青枣“蜜枣”在漳州的引种
试验初报[J]. 中国南方果树, 2012, 41(2): 95–97
Lu X K, Lin Q H, Zhuang W B. Introduction experiment re-
port of Ziziphus mauritiana “candied jujube” in Zhangzhou[J].
South China Fruits, 2012, 41(2): 95–97
[19] 徐宗焕, 陈家金, 林俩法, 等. GIS 在漳州青枣低温害分析
中的应用[J]. 广西气象, 2005, 26(增刊): 12–13
Xu Z H, Chen J J, Lin L F, et al. GIS application to the anal-
ysis of low-temperature damage of Ziziphus mauritiana in
Zhangzhou[J]. Journal of Guangxi Meteorology, 2005, 26
(Supplement): 12–13
[20] 陈清西. 台湾青枣无公害栽培[M]. 福州: 福建科学技术出
版社, 2006: 18–19
Chen Q X. Pollutant-free cultivation for Ziziphus mauritiana
[M]. Fuzhou: Fujian Science & Technology Press, 2006:
18–19
[21] 黄光民, 刘尉, 王广伦. 低温对广东果树、花卉形态及花期
的影响[J]. 广东气象, 2010, 32(5): 39–41
Huang G M, Liu W, Wang G L. Influences of low-temperature
on fruit tree, flower morphology and blooming period[J].
Guangdong Meteorology, 2010, 32(5): 39–41