广西防城港市，位于20°36′N～22°22′N、107°28′E～108°36′E之间，属亚热带季风性气候，光热条件较好。1992—2011年（20年）的年平均气温为22.8 ℃，年累积降雨量为2 427.5 mm，有气候暖干化的趋势（图1）。防城港市现有红树林面积约2 059 hm²，防城港江山半岛西岸边的石角，红树林连片面积超过1 000 hm²，是北仑河口国家级自然保护区的重要组成部分，常见树种有木榄（Bruguiera gymnorhiza）、白骨壤（Avicennia marina）、秋茄（Kandelia obovata）和红海榄（Rhizophora stylosa）等。主要灾害性天气气候事件有暴雨洪涝、台风、干旱、高温、低温冷害、大风、冰雹和雷暴等，其中，对红树林造成危害的以高温、少雨、热带气旋、大风为主^[11-12]。

图  1  防城港市1992—2011年年平均气温及年累积降雨量

1992—2011年及2014—2019年气象数据来自国家气象中心，要素包括逐日最高气温、最低气温、平均气温、降雨量、空气相对湿度等。2015—2019年广州小斑螟和柚木肖弄蝶夜蛾虫害发生面积来自广西红树林研究中心。

采用农业气候分析方法，分析红树林病虫害发生发展过程中（当年5—6月）及发展前的7个月（上年10月—当年4月）的气温、降水、相对湿度、风速等气象环境因子特征。

气温相关指标选取月（旬）平均最高温、平均最低温、平均气温和平均气温>30 ℃日数。湿度相关指标选取月（旬）平均相对湿度、相对湿度<70%日数、相对湿度>80%日数和相对湿度>90%日数。降雨相关指标选取月（旬）累积降雨及平均降雨。风速相关指标选取月（旬）平均风速。

考虑到昆虫对温度和湿度条件的综合要求，采用温湿系数表征温湿综合条件（温湿系数=相对湿度/温度）

广州小斑螟和柚木肖弄蝶夜蛾发生程度以年发生面积为依据，因子普查采用Pearson相关分析法，分析红树林病虫害发生发展过程中（当年5—6月）及发展前的7个月（上年10月—当年4月）不同月、旬尺度上的气温、相对湿度、降水、风速等气象环境因子与受灾面积的相关关系，揭示与广州小斑螟和柚木肖弄蝶夜蛾发生程度显著相关的关键气象因子及影响时段，相关系数的检验采用双尾t检验。

在旬尺度上，选取影响虫害发生的关键气象因子为自变量，以虫害面积为因变量建立多元回归模型。模型如下：

式中，S为虫害面积，$X_{1}、X_{2}、 X_{3} \cdots \cdots X_{\mathrm{k}} $为影响虫害发生的关键气象因子，b、a₁、a₂、$a_{3} \cdots \cdots a_{\mathrm{k}} $为求得的回归参数。

温度方面，2015—2019年虫害发生期间（每年5—6月）及发生的前7个月（上年 10月开始）中，6月的平均气温最高且大于30 ℃日数最多，为29.0 ℃和8 d；1月平均气温最低，为15.7 ℃。降雨方面，6月的降雨最高，为411.2 mm；当年2月的降雨最低，为32.7 mm。风速方面，上年12月的日平均风速最高，为3.1 m·s⁻¹；当年3月的日平均风速最低，为2.5 m·s⁻¹。相对湿度方面，3月的日平均相对湿度最高，为85.6%，同时相对湿度>80%和>90%的日数也相对较多，分别为24 d和11 d。上年12月相对湿度最低，同时相对湿度<70%日数也最多，为13 d（表1）。

防城港市灾情较为严重（受灾面积>66.67 hm²）的年份为2015年和2016年，较轻（受灾面积<6.67 hm²）年份为2018年，分析对比了2种受灾类型年份及历史20年（1992—2011年）降雨、风速、温度和相对湿度各项指标的差异。温度方面，受灾年份5月（虫害发生月）的日最高气温、最低气温及平均气温均高于历史40年均值。且受灾较重年份5月和6月（均为虫害发生月）日平均气温>30 ℃日数显著高于受灾较轻年份及历史水平，即受灾较重年份虫害发生月的总体温度较高（图2）。

图  2  受灾较重年份、较轻年份及历史20年虫害发生前7个月至虫害发生月每月平均最高温、最低温、平均气温及日平均气温大于30 ℃日数。相同字母表示无显著性差异，不同字母之间表示显著性差异，下同。

从图3可知，相对湿度方面，灾情重的年份4月和6月相对湿度显著高于灾情较轻年份，分别高6.5%和7.7%。从上年10月—当年6月灾情较轻年份的相对湿度>90%日数显著低于灾情较重年份，平均低98%。

图  3  受灾较重年份、较轻年份及历史20年虫害发生前7个月至虫害发生月每月平均相对湿度、相对湿度<70%日数、相对湿度>80%日数%及相对湿度>90%日数

从图4可知，降雨方面，虫害发生月（5月和6月），灾情重和轻的年份降雨量分别比历史年份低28.9%和41.6%。风速方面，受灾年份各月平均风速均低于历史年份。

图  4  受灾较重年份、较轻年份及历史20年虫害发生前7个月至虫害发生月每月累积降雨及平均风速

在温度方面，2015—2019年防城港市广州小斑螟的发生与5月最低气温和平均气温呈显著负相关且相关性最高，相关系数分别为−0.94和−0.95。与4月平均气温呈正相关，相关系数为0.63。湿度方面，虫害的发生显著受3月的湿度影响，与3月相对湿度<70%日数呈显著负相关，与相对湿度>80%日数呈显著正相关，相关系数分别为−0.98和0.98。与5月温湿系数呈正相关，但与上年10—11月的温湿系数呈负相关，说明害虫发生的不同阶段所需的环境温度和湿度配合是不同的。降雨方面，广州小斑螟虫害的发生与3月降雨相关性最高且呈显著正相关，相关系数为0.98。风速方面，显著受5月日平均风速的影响，且呈负相关关系，相关系数为−0.74（表2）。

柚木肖弄蝶夜蛾的发生与5月、6月（虫害发生月）的平均气温>30 ℃日数呈显著正相关，相关系数分别为0.84和0.85。与5月相对湿度>80%日数呈显著正相关，相关系数为0.87。降雨方面与广州小斑螟相同，同样与3月降雨相关性最高，相关系数为−0.81。风速方面显著受6月平均风速影响，相关系数为0.80（表3）。

综上，基于月尺度的分析得到，防城港市2种虫害的发生受虫害发生月（5—6月）及当年3月各气象因子影响较大。

基于月尺度的分析结果，在旬尺度上分析了虫害面积与虫害当年2、3、5和6月各旬气象要素的相关性。在温度方面，2015—2019年广州小斑螟的发生与5月下旬和6月上旬最高温呈显著负相关，相关系数分别为−0.97和−0.94。湿度方面，该虫害发生显著受2月上旬和3月下旬平均相对湿度的影响，且综合其他月份相对湿度<70%日数及相对湿度>80%日数的相关性分析，可知广州小斑螟的发生面积随相对湿度的增高而增加。其发生还显著受3月上旬和5月下旬温湿系数影响，相关系数分别为0.90和0.97。降雨方面显著受5月上旬降雨影响，且两者呈负相关，相关系数为−0.94。风速方面显著受3月下旬风速影响，相关系数为0.94（表4）。

柚木肖弄蝶夜蛾的发生在温度上显著受3月上旬和5月下旬影响，相关系数分别为−0.77和0.99。湿度方面受5月中下旬相对湿度、和2月上旬的温湿系数的影响，且随湿度的增高受虫害程度增大。该虫害发生还与2月下旬降雨呈正相关（表5）。

基于旬尺度相关性分析结果，采用多元回归方程，选取显著影响虫害发生的气象因子和影响时段，构建了广州小斑螟和柚木肖弄蝶夜蛾发生面积（S）与气象因子的关系（表6）。

本研究结果表明，广州小斑螟和柚木肖弄蝶夜蛾两种虫害均显著受温度和湿度的影响。虫害与气温、湿度变化关系密切。在自然情况下，大气的温度和湿度这两个因素是相互关联的，而且总是共同综合对虫害起作用，因而本研究引用温湿系数来描述昆虫对温湿条件的综合要求。温度可使害虫内分泌激素代谢等方面受到影响，进而影响害虫生长发育的快慢、迁入迁出的时期和繁殖的世代数^[13]。且害虫一般多为变温生物，温度对其发生及分布有着较大影响^[7]。在适宜的温度范围内，随着温度的升高，害虫的生长发育速度加快，生育期缩短，种群数量、害虫发生的世代数、危害时长和范围随之增加^[14]，如稻飞虱的发生显著受冬季最低温度的影响，冬季最低温度上升有利于提高虫卵越冬存活率^[15]。当温度超出适宜范围时，则会抑制害虫的繁殖，尤其极端低温会破坏虫体的生理结构，使虫体细胞和组织产生不可逆的变化，从而抑制虫害发生^[16]，如2001年我国北方地区5月中下旬出现了大范围晴热高温天气，抑制了麦蚜的种群增长^[17]；再如南美洲斑潜蝇种群的净增殖率遇高温会急剧下降^[18]等。本研究区域内广州小斑螟的发生与5月下旬—6月上旬温度呈负相关，可能由于当温度较高时，害虫取食较多而发育加快，寿命相对缩短。而柚木肖弄蝶夜蛾则与温度呈正相关，其原因可能是温度升高加速了害虫体内酶的活性和生长发育，延长了危害期并扩大了危害范围和程度。此外，本研究中2种害虫对温度的要求有所差异，可能是柚木肖弄蝶夜蛾生长的适宜温度高于广州小斑螟。近几十年来全球气候变化明显，主要体现为气温逐渐升高、极端天气事件增多，这样的气候变化趋势间接提高了虫害率，并增大了虫害发生频率、扩散面积和程度。统计资料显示，近年来我国林业虫害逐年递增，受中等虫害的地区高达50%^[19]，因而在全球变暖的大背景下，更加有必要做好林业虫害防治工作，减缓气候变化带来的负面影响。

湿度和降雨量可影响害虫的迁徙能力和水分平衡，引起其质量变化，通过降水强度和持续时间共同对种群数量及幼虫成活率造成影响^[7]。同时，湿度还可以引发微生物的流行，通过影响寄主植物含水量从而影响害虫的进食^[20]。广州小斑螟和柚木肖弄蝶夜蛾的发生均受当年春季降雨及相对湿度的影响，这可能由于早春湿度对越冬结束后处于复苏状态下的低龄幼虫、卵、茧蛹体内水分调节影响较大，即适宜湿度会复苏害虫体内湿度并加速害虫的繁殖^[7]。且由于雨水的适度冲刷作用，大量下沉降落的害虫也加大了初始虫源基数^[21-22]，类似的影响事实还有稻纵卷叶螟、褐飞虱、白背飞虱、黏虫等^[23]，如耒阳市2006年6月中旬的暴雨引起稻纵卷叶螟蛾数量突增^[24]，但过强和长时间的降雨则会对害虫产生较大的机械伤害，从而使虫害减少^[7]。气候变化背景下，极端天气的增多导致雨季增加，也会导致害虫发生频率和面积增大。

风会影响昆虫的传播，适宜的风速可以帮助昆虫飞翔和迁移，但风速过大则会阻碍一些昆虫的活动，如飞蝗的“小风迎风飞翔，风力稍大顺风飞翔，风力过大则停止飞翔”，因此飞蝗的分布范围常取决于风速和风向^[25]。本研究中柚木肖弄蝶夜蛾的发生受6月风速的显著影响，这可能由于6月风速与害虫取食等活动的关系十分密切，且幼虫多集中于树枝末梢，受风的干扰和吹动即吐丝下垂，飘落到其他树枝处，影响其迁飞和扩散^[26]。类型的现象还如稻水象甲^[27]及粘虫^[28]在微风条件下出现大量扩散和迁飞等。

虫害发生前的气候条件有时会影响到下一年的森林虫害的发生，即存在“滞后效应”^[29]。本研究中，广州小斑螟和柚木肖弄蝶夜蛾的发生均受当年春季（2—3月）气象因子的影响，这说明环境要素（尤其温度和湿度）的长期作用影响着红树林虫害的发生。因此，未来的研究工作可进一步扩大研究时段，逐步探究积温、干燥度等因素对虫害发生的长期效应及影响。

本研究由于数据量有限，仅考虑了一部分地面气象条件对广州小斑螟和柚木肖弄蝶夜蛾发生和发展的影响，而其他与人类活动相关的要素如林木种类、区域经济社会活动、虫害防治措施等也会对虫害的发生产生一定影响^[30]，因此，今后的研究工作可结合植被参数等遥感监测数据，综合探究各要素对虫害发展的影响，建立监测预报模型，并用于评估危害程度。同时要加强虫害防治技术，综合气候因素的影响，灵活采用生物防治、化学药剂防治、物理防治等手段，进而及时防控、高效保护有限的红树林资源。

（1）防城港市红树林受虫灾程度较重年份的当年5—6月（虫害发生月）气温和相对湿度高于虫灾程度较轻年份。

（2）月尺度上，防城港市广州小斑螟的发生受3月降雨和5月温度影响较大，柚木肖弄蝶夜蛾的发生较大程度受5月温度及相对湿度的影响。

（3）旬尺度上，在虫害发生期，防城港市广州小斑螟的发生与5月下旬的温度呈负相关，而柚木肖弄蝶夜蛾的发生与该时段的温度呈正相关。在虫害发生前期，广州小斑螟的发生受当年3月上旬和下旬的相对湿度及风速影响，柚木肖弄蝶夜蛾的发生受当年2月上旬和下旬降雨的显著影响。