海南两系杂交稻制种关键发育期气候风险区划

吕 润，车秀芬，吴 慧，陈小敏，张亚杰，邹海平，白 蕤

（1.海南省气候中心，海口 570203;2.海南省南海气象防灾减灾重点实验室，海口 570203）

相比于三系法水稻，两系杂交稻具有简化繁殖、杂种优势更明显、降低种子生产成本等优势[1-2]，但由于两系杂交水稻制种的育性转换受光温条件的影响较大，光温条件的变化容易导致杂交稻不育系自交结实，从而对制种产量造成严重的影响[3]，海南两系杂交稻制种对气候条件有着极为严格的要求[4-6]，气候条件对制种产量和杂交种子的纯度均有一定的制约[7]。开展两系杂交稻制种气候风险区划将为海南岛南繁水稻制种布局的优化提供重要依据[8]。国内对两系杂交水稻育种区划研究颇多，主要集中在气候条件、气象安全和资源利用等方面。殷剑敏等[9]、汪扩军等[10]、帅细强等[11]均通过对气候资料的系统分析，分别制作了江西省、湖南省杂交稻制种气候风险区划图。

海南作为我国重要的南繁制种基地，在杂交水稻种子生产中起着重要的作用。但是影响海南水稻南繁制种的气象灾害复杂[12-13]，制约着南繁制种面积的扩大。目前大多风险区划研究所使用的气象站数据主要为国家气象观测站数据，站点较少、气象资料稀疏，随着地面区域自动气象观测站在全国的建设和发展，海南岛全部乡镇均已覆盖区域自动气象观测站，相比只利用国家气象站进行插值，区域气象站插值具有高精准、高空间分辨率等优势[14]。因此，本研究使用海南岛南部8个市县的144个区域自动气象观测站点历史资料，在数据质量控制的基础上，利用农业气候区划技术与方法[15-17]，对海南两系杂交稻制种关键发育期进行精细化气候风险区划技术研究，旨在为海南两系杂交稻制种产业区域布局提供参考依据。

1 研究区域、数据与方法

1.1 研究区域海南岛南部地区，包括昌江县、保亭县、乐东县、东方市、五指山市、三亚市、陵水县和万宁市。

1.2 气象资料气象资料为海南岛南部地区144个区域自动气象观测站2006—2020年逐日降水量资料、最低气温、最高气温、平均气温；地理信息资料为海拔高度、经度和纬度等数据。

1.3 地理信息资料1∶25万海南岛地理数据由国家基础地理信息中心提供；数字高程模型数据（DEM）、海南岛行政边界等资料，空间分辨率30 m。

1.4 区域自动气象观测站气象数据质量控制方法为筛选检测区域自动气象观测数据中隐藏的错误数据，从而保证气象资料的可靠性，采用数据质量控制方法[18]。（1）内部一致性检查：检查不同要素之间的物理联系是否符合，如是否符合日最大风速>日平均风速；（2）时间一致性检查：检查气象要素随时间的变化是否符合一定的规律；（3）空间一致性检查：检查气象参数在空间上的分布特点是否满足一定的相关性，本研究采用Madsen-Allerupt方法、空间极值比较法和空间平均值比较法；（4）界限值检查：包括区域界限值、气候学界限值及台站界限值检查。

1.5 精细化气候风险区划指标与方法海南岛南部地区光照充足，两系杂交稻制种的气象限制性因素主要为温度和降水，参考杂交稻制种国家标准[19]，如雷东阳等[20]、汪扩军等[21]、左钦月等[22]提出的两系杂交水稻制种气候适宜性判断标准，并通过咨询海南当地两系杂交稻制种专家[23]，结合海南岛的气象条件和生产实际，将两系杂交稻南繁分为育性敏感期和扬花授粉期2个关键时期，育性敏感期大部分在4月上旬至4月中旬，扬花授粉期大部分在4月下旬至5月上旬，制定安全性与适宜性判断标准（表1）。

表1 两系法杂交稻南繁育性敏感期和扬花授粉期气候风险等级

育性转换气候风险由不育系育性敏感期间出现超过临界温度指标气象事件的概率确定。

式中，P为育性敏感期期间出现超过不育临界温度气象事件的概率；N为育性敏感期期间出现超过不育临界温度气象事件的天数；L为气候资料总天数。

扬花授粉期气候风险由气候风险指数确定。

式中，H为扬花授粉期气候风险指数；L为气象资料总天数；i为第i个年份（i=1,2,…L）；M为气象灾害种类的数量；j为第j种气象灾害（j=1,2,…M）；Hij为第i年第j种气象灾害的强度。在扬花授粉期内，第i年第j种灾害的持续时间不足3 d时，认为灾害没有发生，记Hij=0；当持续时间等于3 d时，记Hij=1；当持续时间大于3 d时，每超过1 d则Hij加上1；如果多个灾害同时发生，则先分别用上述方法计算各自的强度值，再累加求得Hij值。本研究为扬花授粉期，主要气象灾害指标为连续3 d日平均气温≤24℃、3 d日最高气温≥35℃和连续3 d阴雨天。

两系杂交稻制种综合气候风险等级由育性敏感期气候风险等级和扬花授粉期气候风险等级共同决定，当等级不一致时，取最低气候风险等级。等级指标如表2所示。

表2 两系法杂交稻南繁制种气候风险等级

利用ArcGIS以反距离权重法（Inverse Distance Weighting，IDW）将气候风险指数进行空间插值计算，从实际情况考虑，由于海拔较高的山区不适宜大面积制种，故将海拔在200～400 m之间地区的气候风险值增加一倍，海拔大于400 m的地区则划为高风险区。最终得到两系杂交稻制种关键发育期气候风险区划。

2 结果与分析

2.1 育性转换安全期气候风险区划当不育临界温度22℃时，除部分高海拔地区，绝大部分地区均为较低风险区或极低风险区。其中极低风险区主要在东方市西南部地区，乐东县西部、南部、东北部，三亚市西部、南部、东部地区，保亭县东部、中部、南部地区，陵水县中部、南部地区。较低风险区主要在昌江县西部、中部、北部地区，东方市北部、东部地区，乐东县东部地区，三亚市北部地区，保亭县北部地区，万宁市绝大部分地区。中高风险区主要分布在海拔较高的山区，如五指山市、昌江县东南部、乐东县北部、三亚市北部、保亭县西部、陵水县北部和万宁市西部部分地区。不育临界温度23℃时，极低风险区主要在乐东县西南部、三亚市南部和东部、陵水县南部，较低风险区主要在东方市西部和南部、乐东县西北部和东部、三亚市西部、保亭县中部和东部、陵水县中部、万宁市东南部。中高风险区主要分布在昌江县、东方市东北部、万宁市中北部和海拔较高的山区。当不育临界温度23.5℃时，仅有乐东县和三亚市小部分地区为极低风险区，较低风险区主要在东方市西南部、乐东县西部和南部、三亚市南部和东部、陵水县中部和南部地区。昌江县、五指山市、保亭县、万宁市以及东方市、乐东县、三亚市、陵水县的高海拔山区为中高风险区。当不育临界温度24℃时，仅有乐东县西部、三亚南部和东部、陵水南部为较低风险区，其余地区均为中高风险地区，因此对临界温度为24℃的不育系制种，须慎重考虑，尽量选择在乐东、三亚、陵水海拔较低的地区。

2.2 扬花授粉期气候风险区划连续3 d日平均气温≤24℃结果表明，绝大部分地区为极低风险区，中高风险区最要集中在五指山市、昌江县东南部、东方市南部、乐东县西北部和三亚市北部的高海拔山区。连续 3 d日最高气温≥35℃结果表明，极低风险区主要分布在乐东县南部、保亭县中部、三亚市大部、陵水县大部和万宁市大部，较低风险区主要分布在昌江县西部、东方市西部和东部、乐东县东北部，中高风险区主要分布在昌江县大部、东方市大部、乐东县北部、五指山市、保亭县西部和海拔较高的山区。连续3 d阴雨天结果表明，极低风险区主要分布在昌江县西部、东方市西部、乐东县西部和南部、三亚市大部、陵水县南部，较低风险区主要分布在昌江县中部、东方市东部、乐东县东北部、保亭县东部、陵水县北部、万宁市大部，中高风险区主要分布在五指山市、昌江县东部、东方市东部和南部、乐东县北部和东部、三亚市北部、保亭县西部和北部、陵水县北部和海拔较高的山区。综合所得区划的结果表明，极低风险区主要分布在乐东县南部、三亚市南部和东部、陵水县南部，较低风险区主要分布在昌江县西部、东方市西北部、乐东县西北部和东北部、保亭县东部、陵水县北部、万宁市部分地区，中高风险区主要分布在五指山市、昌江县中部和东部、东方市大部、乐东县北部和东部、三亚市北部、保亭县西部、陵水县北部、万宁市北部和中部。

2.3 两系杂交稻制种综合气候风险区划当不育临界温度22℃时，极低风险区主要在乐东县大部、三亚市大部、保亭县大部和陵水县中南部地区，面积约为区划区域的26%。较低风险区主要在昌江县大部、东方市大部、乐东县东部、三亚市北部、保亭县北部和万宁市绝大部分地区，面积约为区划区域的39%。中高风险区主要分布在海拔较高的山区，如五指山市、昌江县东南部、乐东县北部、三亚市北部、保亭县西部、陵水县北部和万宁市西部部分地区，面积约为区划区域的35%。不育临界温度23℃时，极低风险区主要在乐东县西南部、三亚市东南部部、陵水县南部，面积约为区划区域的9%，较低风险区主要在东方市西南部、乐东县西北部和东部、三亚市西部、保亭县中东部、陵水县中部、万宁市东南部，面积约为区划区域的33%。中高风险区主要分布在昌江县、东方市东北部、万宁市中北部和海拔较高的山区，面积约为区划区域的58%。当不育临界温度23.5℃时，仅有乐东县、三亚市小部分地区为极低风险区，面积约为区划区域的2%，较低风险区主要在东方市西南部、乐东县西南部、三亚市东南部、陵水县中南部地区，面积约为区划区域的24%。昌江县、五指山市、保亭县、万宁市及东方市、乐东县、三亚市、陵水县的高海拔山区为中高风险区，面积约为区划区域的74%。当不育临界温度24℃时，仅有乐东县西部、三亚东南部、陵水南部为较低风险区，面积约为区划区域的14%，其余地区均为中高风险地区。

3 讨 论

根据两系杂交稻制种关键发育期对气象条件的要求，采用海南岛南部地区区域自动气象观测站的观测资料，开展两系杂交稻制种关键发育期气候风险区划，结果表明，两系杂交稻制种综合气候风险随不育临界温度的升高而增加。不育临界温度22℃时，两系杂交稻制种综合气候低风险区分布于昌江、东方、乐东、三亚、保亭、陵水和万宁；不育临界温度23℃时，两系杂交稻制种综合气候低风险区分布于东方、乐东、三亚、保亭和陵水大部分地区；不育临界温度23.5℃时，两系杂交稻制种综合气候低风险区分布于乐东、三亚和陵水大部分地区；不育临界温度24℃时，两系杂交稻制种综合气候低风险区分布于乐东、三亚和陵水沿海地区。目前，在杂交稻制种上广泛使用的两系不育系中起点温度大多数在23.5℃左右，不育起点温度低于23℃的较少。从上述分析结果可知，海南两系杂交水稻制种育性转换安全期的最适宜地区主要在东方市西部、乐东县西部南部东北部、三亚市西南部东部、陵水县南部中部地区。据资料[24-25]显示，目前海南南繁制种水稻主要分布于三亚、陵水、乐东和东方地区，李东秀等[26]通过在三亚、乐东、昌江、临高进行试制试验，结果表明三亚、乐东的产量较高，由此可见，本区划结果与实际水稻制种情况较吻合。

本研究采用订正后的区域自动气象观测站观测资料开展区划工作，克服了因观测站点稀疏，采用插值方法开展精细化区划工作存在准确性差、精度低等问题。海南岛三亚市、陵水县和乐东县是全国最大的南繁育种研究试验基地和制种生产基地，从区划结果来看，与目前海南岛杂交稻制种的情况基本吻合，且南繁为海南自由贸易港建设的重点工作之一，本区划将为海南两系杂交稻制种种植规划提供及时科学的技术支撑。本研究杂交稻制种的关键期时段和气候风险指标的选择主要参考前人研究及专家经验，未经过大田对比试验及田间观测的验证，且未考虑种植技术、土地利用率、土壤等因素。特别是海南水稻制种多在冬季，而海南冬季易发生干旱[27]，因此在实际水稻制种地选择过程中需要关注灌溉用水保障。今后的研究应通过田间试验更新完善制种区划指标，充分考虑上述因子对海南两系杂交稻制种的综合影响开展气候风险区划。