气象观测资料：国家站资料选取旺苍及附近4个气象站点1981-2017年逐日气温、日照时数和降水等气象资料，区域站资料选取旺苍县19个区域站点2012-2017年逐日气温、降水等气象资料，以上资料均由四川省气候中心提供(站点分布如图 1所示).

水稻数据：水稻生育期资料来自四川省水稻和玉米农业气象观测站点数据，县级水稻产量、种植面积等数据来自2017年《四川省统计年鉴》，乡镇水稻种植面积来自2017年《旺苍县统计年鉴》.

社会经济数据来自2017年《旺苍县统计年鉴》，地理信息资料来自国家基础地理信息中心提供的1∶250 000的数字高程模型(DEM)资料.

根据地理地貌类型以及种植区形成的现状，四川盆地分为盆南丘陵区、盆中浅丘区、盆西平丘区、盆周边缘山地区、盆东平行岭谷区5个种植区，其中旺苍县及本研究所选附近气象站点均属于盆周边缘山地区^[20]，片区气候具有相似性，由于乡镇水稻产量资料局限性，本研究通过提取同片区(广元、苍溪、旺苍、南江)主要气象灾害类型确定区划因子，从而对旺苍县水稻气象灾害开展评估.

本研究采用灾害风险评估原理对旺苍县水稻进行风险评估，综合考虑旺苍县水稻种植情况、地形地貌特征和经济数据等资料，基于GIS技术平台，对旺苍县水稻不同生育期气象灾害进行风险评估，采用加权综合评价法和层次分析法建立旺苍水稻综合气象灾害风险评估体系，从而对旺苍县水稻进行气象灾害风险评估，利用GIS技术中自然断点法确定旺苍县水稻气象灾害风险等级，并将其划分为低风险性、次低风险性、中风险性、次高风险性、高风险性5个等级.

灾害风险评估原理：自然灾害风险是致灾因子、孕灾环境、承灾体特征和防灾减灾能力多方面因素综合作用的结果，其中致灾因子危险性由气象灾害发生种类、频次、强度等反映. 孕灾环境是孕育灾害的外部自然环境. 承灾体特征是指暴露在自然灾害环境下的承灾体的危险程度，主要由脆弱性和暴露性两方面构成. 防灾减灾能力作为社会属性，由经济水平、文化水平、公共设施以及当地政策等因素共同构成^[21].

式中：FDRI为危险性；D为致灾因子危险性；V为承灾体脆弱性；E为承灾体暴露性；S为孕灾环境敏感性；P为防灾减灾能力.

加权综合评价法是根据各指标因子对研究对象的影响程度赋予各因子不同的数量化指标，并进行综合评价的方法^[21].

式中：f为综合评价指标；D_i为第i指标归一化值；W_i为第i指标权重；m为评价指标个数.

层次分析法是通过专家打分的方法对同一层次指标的重要性进行两两比较，并将比较结果A_ij(i，j=1，2，…，n)写成n*n阶矩阵形式，构建判断矩阵，在建立有序递阶的基础上，计算得出矩阵中各指标权重向量和特征根，计算各层指标权重系数^[22].

自然断点法是一种根据数值统计分布规律分级和分类的统计方法，利用数列中的自然转折点、特征点，将相似值进行最恰当地分组，并使各个类之间的差异最大化.

水稻属喜温、喜湿、短日照作物，研究认为水稻适宜出苗温度为12~18 ℃，日均温20~25 ℃时稻苗生长旺盛，温度过高或过低都会影响水稻的正常生长. 温度过高时，水稻无法正常开花结实，空秕粒所占比重上升，结实率下降，水稻产量降低，温度过低也会影响水稻的秧苗素质，孕穗期冷害会使得籽粒灌浆过程受阻，千粒质量下降，从而造成水稻产量大幅下降. 水稻生长过程中对水分要求也比较严格，水分不足或过多，都会影响水稻秧苗的正常生长. 综合水稻生长发育所需环境要求和旺苍县特殊地理地形，初选低温冷害(灌浆期、乳熟期)、高温(孕穗期、抽穗期)、暴雨(抽穗期、灌浆期)、连阴雨(抽穗期、灌浆期)4类8种灾害，并分析不同灾害对水稻产量的影响，具体灾害阈值如表 1所示.

采用滑动平均方法从水稻单产中分离得出各地水稻气象产量^[23]，并与各气象灾害进行相关分析(表 2). 由表 2可知，大部分站点水稻气象产量与气象灾害之间表现为负相关关系，其中乳熟期低温冷害、抽穗期连阴雨和灌浆期连阴雨对水稻产量影响较显著，综合相关分析结果和旺苍水稻实际生产，选取抽穗期高温、抽穗期暴雨、抽穗期连阴雨、灌浆期冷害、灌浆期连阴雨和乳熟期冷害6种灾害作为旺苍水稻气象灾害评估的主要影响因素.

旺苍县以山地为主，壑谷纵横，地形复杂，气象灾害与海拔、地形地貌等环境条件密切相关. 本研究提取海拔、坡度、坡向3个地形因子(图 2)，根据不同气象灾害致灾机理，选取不同地形指标来表示孕灾环境敏感性，参考专家打分法分别赋予不同地形因子对孕灾环境敏感性的影响指数(表 3).

承灾体脆弱性是指承灾体在遭受自然灾害打击时的易损程度，承灾体脆弱性越高，易损程度越严重. 本研究从作物敏感性角度分析旺苍县水稻种植环境脆弱性程度，采用历史产量的变异程度作为脆弱性评价指标^[24]. 利用自然断点法划分为5个等级，得到图 3，其中0≤V＜0.2为低脆弱性，0.2≤V＜0.39为次低脆弱性，0.39≤V＜0.57为中脆弱性，0.57≤V＜0.75为次高脆弱性，0.75≤V≤1为高脆弱性. 由图 3可知旺苍中西部地区水稻生产脆弱性较低，易损度低，檬子乡、大河乡、万山乡、五权镇、大德乡、金溪镇等东部乡镇水稻生产脆弱性较高，易损程度高.

承灾体暴露性是指暴露在自然灾害环境下的水稻种植程度，水稻种植程度越大，承灾体暴露性越高，受气象灾害威胁越严重，本研究用水稻种植比例来表征承灾体暴露性指标^[24]，绘制水稻种植暴露性分布图(图 4)，图中，0≤E＜0.21为低暴露性，0.21≤E＜0.37为次低暴露性，0.37≤E＜0.57为中暴露性，0.57≤E＜0.76为次高暴露性，0.76≤E≤1为高暴露性. 可知以白水镇至大德乡为界，南部地区水稻种植比例较大，暴露性较高，北部地区山地较多，水稻种植面积小，承灾体暴露性低.

防灾减灾能力是指风险承担者抵御自然灾害所制定的方针、政策、行动等一系列活动的总称，防灾减灾能力越强，灾害所造成的潜在损失越小. 本研究选取农村居民人均纯收入来表征当地防灾减灾能力，0≤P＜0.33为低防灾减灾能力，0.33≤P＜0.52为次低防灾减灾能力，0.52≤P＜0.65为中防灾减灾能力，0.65≤P＜0.78为次高防灾减灾能力，0.78≤P≤1为高防灾减灾能力. 可知防灾减灾能力较强地区分布在白水镇至金溪镇、燕子乡等地区，檬子乡、正源乡防灾减灾能力最弱(图 5).

采用层次分析法(AHP)，分别计算单灾种风险评估指标权重(表 4)和综合气象灾害各指标权重系数(表 5).

抽穗期是水稻高温热害的敏感期，温度过高直接影响水稻幼穗分化，导致最终产量下降. 以抽穗期高温年平均发生日数为致灾因子，综合考虑高程、坡向等地形因子和承灾体脆弱性、暴露性以及防灾减灾能力，根据自然断点法得出旺苍县水稻抽穗期高温风险分布，如图 6所示，图中5个等级分别为：低风险性(0.25~0.40)、次低风险性(0.41~0.50)、中风险性(0.51~0.60)、次高风险性(0.61~0.70)、高风险性(0.71~0.88)，由图 6可知，旺苍县水稻抽穗期高温风险分布呈现较强的层次性，西南高东北低，其中东河镇、柳溪乡、嘉川镇、黄洋镇、龙凤乡、尚武镇等地高温风险性最高.

水稻不同生育阶段耐涝性不同，据研究，水稻苗期、分蘖期耐涝性较强，抽穗开花期耐涝性最差，此时遭遇暴雨会影响水稻小穗分化和生殖细胞发育^[25-27]. 分析旺苍县抽穗期暴雨发生频率，地形因子，承灾体脆弱性、暴露性和当地防灾减灾能力，绘制旺苍县水稻抽穗期暴雨风险分布图(图 7). 分别为：低风险性(0.25~0.38)、次低风险性(0.39~0.45)、中风险性(0.46~0.52)、次高风险性(0.53~0.60)、高风险性(0.61~0.75)，由图 7可知，旺苍县水稻抽穗期暴雨风险分布自西南向东北递减，水稻抽穗期暴雨高风险区域主要分布在东河镇、黄洋镇、嘉川镇、柳溪乡、龙凤乡等地区. 中部和东北部等乡镇处于暴雨中风险性或次低风险性地区.

连阴雨是指连续多日的阴雨天气，8月下旬至9月末的秋季连阴雨正值水稻成熟收获的关键时期，连阴雨及其伴随的低温、寡照、高湿等外界环境直接影响水稻干物质积累，导致水稻减产. 旺苍县水稻抽穗期连阴雨风险分布如图 8所示，5个等级分别为：低风险性(0.13~0.29)、次低风险性(0.30~0.37)、中风险性(0.38~0.46)、次高风险性(0.47~0.56)、高风险性(0.57~0.78). 由图 8可知，旺苍县水稻抽穗期连阴雨风险性呈现东高西低的变化特征，且大部分地区处于次低、低风险性区域，仅檬子乡至三江镇以东的乡镇，风险性等级较高，其中大河乡、大德乡、金溪镇东部、五权镇东部风险等级最高.

考虑灌浆期冷害年平均日数、水稻产量变异程度以及旺苍县海拔和坡向等因素，利用自然断点法对旺苍县水稻进行灌浆期冷害风险评估，分别为：低风险性(0.27~0.40)、次低风险性(0.41~0.47)、中风险性(0.48~0.55)、次高风险性(0.56~0.61)、高风险性(0.62~0.78). 由图 9可知，水稻灌浆期冷害风险呈现自西向东递减的分布. 其中白水镇、尚武镇、枣林镇、嘉川镇、柳溪乡、东河镇、龙凤乡、九龙乡等西南乡镇水稻灌浆期冷害风险最高，次高风险性区域沿万家乡—麻英乡—龙凤乡呈倒V型分布，中风险性区域沿鼓城乡—高阳镇—三江镇分布，东部乡镇水稻灌浆期冷害风险性较低.

灌浆期是连阴雨灾害发生的主要时期，灌浆期连阴雨会影响水稻灌浆结实，增大水稻空秕率，影响水稻的成熟收获，严重时可造成田间籽粒发芽或者籽粒霉烂、病虫害滋生，水稻产量和品质严重下降. 旺苍县灌浆期连阴雨风险评估为：低风险性(0.27~0.41)、次低风险性(0.42~0.48)、中风险性(0.49~0.55)、次高风险性(0.56~0.62)、高风险性(0.63~0.75). 由图 10可知旺苍县西南乡镇大多处于高风险或次高风险区域，其中东河镇、嘉川镇、柳溪乡、尚武镇东部、枣林乡、黄洋镇、龙凤乡风险等级最高，万家乡、盐河乡、天星乡、国华镇、福庆乡等西北部乡镇水稻灌浆期连阴雨风险等级为中级，旺苍东部地区灌浆期连阴雨风险等级较低.

旺苍县水稻乳熟期一般开始于九月中旬，受大气环流控制，这一阶段旺苍县温度较低，水稻主要气象灾害为低温冷害. 综合分析水稻乳熟期低温冷害年平均发生日数、海拔、坡向、水稻产量、种植比例以及防灾减灾能力等相关因素，绘制旺苍县水稻乳熟期冷害风险分布图：低风险性(0.10~0.25)、次低风险性(0.26~0.35)、中风险性(0.36~0.43)、次高风险性(0.44~0.49)、高风险性(0.50~0.62). 由图 11可知，旺苍大部分地区水稻处于乳熟期冷害高风险或次高风险区域，主要分布在中西部乡镇，其中福庆乡—龙凤镇以西乡镇均为高风险性地区. 中风险性地区自鼓城乡至九龙乡呈南北向分布，大河乡、大德乡、金溪镇和五权镇西部水稻乳熟期冷害风险等级最低.

采用多指标综合风险评估法建立旺苍县水稻气象灾害指标体系，通过专家打分和层次分析方法设定不同气象灾害因子权重，将不同气象灾害风险均一化处理后，通过GIS平台进行栅格计算，利用自然段点分级法，将旺苍县水稻气象灾害风险划分为高风险性(1.31~1.62)、次高风险性(1.08~0.30)、中风险性(0.86~1.07)、次低风险性(0.63~0.85)、低风险性(0.34~0.62) 5个等级，得出旺苍县水稻气象灾害风险分布图(图 12). 高风险区：尚武镇南部、枣林乡东部、嘉川镇中南部、东河镇中南部、黄洋镇中南部、龙凤乡西部、柳溪乡. 次高风险区：白水镇、枣林镇西部、尚武镇北部、嘉川镇北部、东河镇北部、高阳镇西南部、黄洋镇北部、普济镇中南部、农建乡大部、化龙乡、木门镇、九龙乡. 中风险区：燕子乡、麻英乡北部、福庆乡、双汇镇、高阳镇中东部、天星乡西南部、国华镇西南部、正源乡南部、普济镇北部、三江镇中南部. 次低风险区：天星乡中北部、万家乡、盐河乡、鼓城乡西部、国华镇中北部、正源乡北部、英萃镇、大两乡大部以及三江镇北部、金溪镇. 低风险区：鼓城乡东部、檬子乡、大河乡、水磨乡、万山乡、五权镇、大德乡、大两乡东部.

本研究基于自然灾害风险评估原理，基于GIS技术综合考虑旺苍县附近国家站和区域站气象观测数据、水稻种植自然环境、旺苍地理地貌特征复杂性、水稻产量变异性、种植面积以及当地防灾减灾能力等相关信息，构建水稻气象灾害风险评估体系，对旺苍县水稻气象灾害风险开展初步评估，具体结果如下：

气象灾害风险评估是致灾因子、孕灾环境、承灾体和防灾减灾能力等因子相互作用下的综合评价的结果. 从致灾因子来看，旺苍县对水稻种植影响较大的气象灾害主要有低温冷害、连阴雨、暴雨以及高温，主要发生在水稻抽穗至乳熟阶段，其中灌浆期低温冷害和抽穗期连阴雨对水稻产量影响最大，高温对旺苍水稻生产影响最小. 从孕灾环境来看，旺苍县以山地为主，地势北高南缓，腹部低平，复杂山地下的孕灾环境对不同灾种的影响存在差异，高海拔、阴面山坡对低温冷害有促进作用，地势陡峭、海拔高地区则有利于暴雨、连阴雨等灾害的发生. 从承灾体来看，暴露性和脆弱性的大小可以用来表征孕灾环境对水稻气象灾害的促进或延缓作用，旺苍县南部乡镇地势平缓，属农业重点发展区，水稻种植面积较大，水稻暴露性较高，脆弱性较小，东北部乡镇海拔较高地势险峻，是生态重点保护区，水稻种植面积少，产量变异程度较大，脆弱性较高但暴露性较小. 防灾减灾能力是区域范围内政策、人力、物力、生产力水平等多种因素综合作用下的结果，东河、嘉川等西南乡镇是旺苍经济发展的核心地区，经济水平高，综合防灾减灾能力强，东北部乡镇经济较薄弱，综合防灾减灾能力弱.

水稻综合气象灾害风险分布受单灾种灾害风险分布影响，旺苍县水稻灾害风险呈现为西南高东北低的分布，这与抽穗期高温、抽穗期暴雨、灌浆期连阴雨、灌浆期冷害以及乳熟期冷害的风险分布类似，但综合气象灾害风险分布较单灾种风险分布层次较分明. 次高、高风险地区大多分布在西南乡镇，主要由致灾因子危险性高和承灾体暴露性高导致，应不断提高农业发展水平，选育抗性强的水稻品种，提高抵御风险能力，中风险地区大多在旺苍中部地区，该区域在生产过程中应加强水稻种植管理，改进种植方式，东北部乡镇水稻综合风险等级低，但地势较高，地形复杂，孕灾环境较为敏感，且防灾减灾能力较弱，也应加以重视.