张 蕾 郭安红 何 亮 侯英雨 赵秀兰 钱永兰 蔡 哲
1 国家气象中心,北京 100081
2 江西省农业气象中心,南昌 330046
提 要: 基于高分辨率的土地利用/土地覆被遥感监测数据,融合下垫面信息、格点化寒露风等级和不同尺度寒露风指数,分析了2020年我国南方双季晚稻区的寒露风特点。结果表明:2020年晚稻寒露风发生时间早,持续时间长,普遍出现轻度至重度寒露风。寒露风发生面积占稻田总面积的24.7%,为2000年以来第三高,发生轻度、中度、重度寒露风的面积比例分别为9.7%、12.8%和4.9%;湖南、江西、浙江发生面积比例超过90.0%,其余省(自治区)不足30.0%;2020年湖南、湖北、安徽、江苏发生面积比例为2000年以来最高,江西发生重度寒露风面积比例为2000年以来最高,浙江、福建、广西发生寒露风以轻度为主,广东寒露风发生面积较小;晚稻区寒露风指数达4.21,为2000年以来第二高;2020年湖南、江西、湖北寒露风指数较高,分别达7.16、7.16、12.59,相应分别为2000年以来第二高、最高、最高。
农业生产与气象、气候条件息息相关,气候变化,尤其是气候变暖背景下极端气候事件的变化会对农业生产产生较大影响。随着全球气候不断变暖,极端温度和极端温度事件的发生特点有所改变(李金洁等,2019);总体而言近年来极端高温发生频率明显增加,持续时间延长(聂羽等,2018),而极端低温的发生总体呈降低趋势,随之造成的农业气象灾害发生频率有减有增。作为我国主要粮食作物之一,水稻生产有其适宜的气候条件范围(易灵伟等,2016;徐敏等,2018),而伴随气候变化所引起的极端温度条件改变会引起灾害胁迫变化(徐敏等,2015;许莹等,2020;杨建莹等,2020;Zhang et al,2018)。双季晚稻是我国南方稻区主栽作物,生产中遭受的主要温度灾害之一是寒露风(吕晓敏和周广胜,2018)。寒露风是指晚稻抽穗杨花期因低温造成扬花受阻、空壳率增加的灾害性天气。受寒露风影响,晚稻生理特征(如叶生理活性下降)、品质构成(如出糙率、精米率、整精米率、垩白度、直链淀粉含量下降,垩白粒率、胶稠度上升)和产量构成(如结实率降低,千粒重减轻,产量下降)等方面会受到不同程度的影响(林洪鑫等,2016;余焰文等,2014)。
对南方双季晚稻寒露风的发生的气候特点,目前已有不少学者从寒露风气象指标、气候特征、风险评估和影响评估等方面开展研究,包括:基于大多数农业气象业务和研究工作者认可的寒露风指标,分析了南方双季晚稻区和省(自治区)等不同尺度下晚稻不同等级寒露风指标发生频率、危害日数、发生等级等时空分布特征(戴剑波等,2012;黄晚华等,2011;刘文英等,2009;刘丹等,2019;罗伯良和李易芝,2015;彭莉莉等,2014;苏荣瑞等,2012);基于寒露风发生特点,构建寒露风风险评价指标,开展了不同区域寒露风风险和影响评估(黄晚华等,2011;田俊和崔海建,2016;王华等,2018;吴立等,2014)。上述针对寒露风的研究较好地总结了寒露风的发生特点和风险水平,对开展寒露风防御工作具有重要指导作用。以往研究多用气象站点的表征指数来代表整个县或者区域的寒露风情况,较少从下垫面角度分析,也很少能反映出寒露风发生面积。近年来,随着多源数据的应用,从面尺度开展寒露风的评估工作逐步开展,Chen et al(2019)利用中国气象局陆面数据同化系统(CLDAS)的格点数据开展了寒露风的监测,统计分析了不同等级寒露风的受影响面积;这项工作拓展了寒露风影响范围的评估,但没有考虑下垫面信息从而难以获取精确的影响范围。2020年9月中下旬,南方双季晚稻区出现了明显的寒露风天气,针对此次典型的寒露风天气过程结合历史寒露风发生特点,本文融合高分辨率土地利用/土地覆被遥感监测数据和气象数据,从寒露风的发生时间、面积、等级、强度等角度精细化分析晚稻寒露风发生特点。
研究所用土地利用数据为中国多时期土地利用/土地覆被遥感监测数据集(CNLUCC),来源于中国科学院资源环境科学数据中心,空间分辨率为30 m×30 m(徐新良等,2018)。该数据集采用中国科学院土地资源分类系统:1级分类划分为6类(耕地、林地、草地、水域、建设用地和未利用土地);2级分类主要根据土地资源的自然属性分为25类(如耕地划分为水田和旱地)。该数据集经过遥感解译,具备较好的应用精度(刘纪远等,2018)。在处理过程中,考虑到目前国家级农业气象灾害监测评估业务服务中产品的空间分辨率为5 km×5 km,因此以5 km×5 km格点为对象,判断该格点内划分的30 m×30 m格点对应CNLUCC数据集中的土地利用类型,如果是水田则计为1,最后水田比例即为格点内水田数与格点包含的30 m×30 m总数的比值(图1)。
图1 双季晚稻区水田分布Fig.1 Location of double-season late rice
研究所用的气象资料包含南方双季晚稻区(湖北南部、江苏南部、安徽南部、上海、浙江、江西、湖南、福建、广东、广西)的789个气象站点2000年以来逐日平均温度、最低温度等。晚稻发育期数据包括晚稻区82个农业气象观测站2000年以来逐年观测的发育期(播种、出苗、三叶、移栽、返青、分蘖、拔节、孕穗、抽穗、乳熟、成熟期)日期。
寒露风等级计算:首先,将计算时段内晚稻区气象站点的逐日观测数据插值到5 km×5 km格点,然后判断格点所处的发育期(发育期用其最邻近的农业气象观测站观测的晚稻发育期代替)是否进入抽穗扬花阶段;若进入,则参考《寒露风等级》(QX/T 94—2008,中国气象局,2008)和《大宗作物气象服务手册》(毛留喜和魏丽,2015)确定的业务服务中的晚稻寒露风等级划分指标(表1),即可判定该格点的寒露风等级。
表1 寒露风等级划分Table 1 Classification of degree of chilling dew wind
寒露风面积提取:在格点化寒露风等级的基础上,与水田面积比例图层相叠加,即格点内寒露风等级发生面积=水田面积比例×寒露风等级×格点面积;省级/区域寒露风面积为省/区域内格点判定的寒露风面积的和。
不同空间尺度寒露风指数构建:在分析时段内,通过寒露风等级划分计算轻度、中度、重度寒露风的天数,考虑不同等级寒露风的影响程度,进一步构建站点寒露风指数(H):
Hi=0.2LlDl+0.3LmDm+0.5LhDh
(1)
式中:Hi为i站的寒露风指数;Ll、Lm、Lh分别为轻度、中度、重度寒露风等级,即对应1、2、3;Dl、Dm、Dh分别为轻度、中度、重度寒露风等级出现的天数。
通过对不同站点寒露风指数进行集成,构建省级寒露风指数(HP):
(2)
式中:Hi为i站的寒露风指数;wi为i站的权重,即为该i站点所在县的水田面积占全省水田总面积的比例,m为省内站点数。
对不同省寒露风指数进行集成,构建晚稻区寒露风指数(HS):
(3)
式中:权重wj为该j省水田面积占晚稻区水田总面积的比例,n为省份个数。
2020年9月13日至月底,南方晚稻区出现寒露风天气,与2000—2019年寒露风的平均发生日期对比来看(图2),湖南北部、湖北南部普遍偏早5 d以上,湖南中部—东南部偏早3~5 d,湖南南部、江西中北部、安徽南部、江苏南部、浙江、广西西北部等地偏早1~2 d。从南方双季晚稻区寒露风发生天数来看,湖北南部、安徽南部、江苏南部、浙江、江西、湖南、福建西北部和广西北部等地出现3 d以上寒露风天气,湖北南部、安徽南部、江西北部、湖南中北部和浙江中部等地达9 d以上。依据寒露风等级判断标准,上述地区出现轻度及以上寒露风,其中湖南中部、江西北部、湖北东南部、安徽南部等地出现重度寒露风。
图2 2020年双季晚稻区寒露风(a)发生日期与2000—2019年平均发生日期对比,(b)持续天数,(c)等级Fig.2 (a) Comparison of occurrence date between chilling dew wind in 2020 and the multi-year average during 2000-2019, (b) lasting days and (c) degree of chilling dew wind in 2020 across double-season late rice region
从不同等级双季晚稻区寒露风发生面积来看(图3),2020年双季晚稻区寒露风发生面积占水田总面积的24.7%,低于2011年的28.3%和2002年的27.5%;其中,轻度、中度、重度寒露风面积比例分别为9.7%、12.8% 和4.9%,分别为2000年以来第七、第三、第三高。
图3 2000—2020年双季晚稻区寒露风发生面积比例Fig.3 The percent of area exposed to chilling dew wind across double-season rice region during 2000-2020
从分省寒露风发生面积来看(图4),2020年江苏、安徽、湖北、湖南、江西、浙江、福建、广东、广西该面积占比分别为11.1%、18.9%、17.7%、99.9%、99.4%、91.8%、23.1%、2.8%、22.2%。2020年湖南寒露风发生面积比例为2000年以来最高,轻度、中度、重度寒露风面积比例分别为16.9%、63.6%和19.4%;江西寒露风发生面积比例为2000年以来第二高,轻度、中度、重度寒露风面积比例分别为36.3%、37.2%和25.9%,重度寒露风面积比例为2000年以来最高;湖北寒露风发生面积比例为2000年以来最高,轻度、中度、重度寒露风面积比例分别为0.0%、14.6%和3.0%;安徽寒露风发生面积比例为2000年以来最高,轻度、中度、重度寒露风面积比例分别为2.4%、12.0%和4.4%;江苏寒露风发生面积比例为2000年以来最高,轻度、中度、重度寒露风面积比例分别为8.8%、2.3%和0.0%;浙江寒露风发生面积比例为2000年以来第四高,轻度、中度、重度寒露风面积比例分别为53.9%、25.9%和12.0%;福建寒露风发生面积比例为2000年以来第五高,轻度、中度、重度寒露风面积比例分别为20.8%、0.4%和1.9%;广东寒露风发生面积较小,仅为轻度寒露风;广西寒露风发生面积比例为2000年以来第三高,轻度、中度、重度寒露风面积比例分别为20.6%、1.6%和0.0%。
图4 2000—2020年(a)江苏,(b)安徽,(c)湖北,(d)湖南,(e)江西,(f)浙江,(g)福建,(h)广东,(i)广西不同等级寒露风发生面积比例Fig.4 The percent of area exposed to chilling dew wind in Jiangsu (a), Anhui (b), Hubei (c), Hunan (d), Jiangxi (e), Zhejiang (f), Fujian (g), Guangdong (h) and Guangxi (i) during 2000-2020
从寒露风指数来看(图5),湖北南部、安徽西南部、江西北部、湖南中西部、浙江中部寒露风指数在5以上,江西中部、湖南东南部、福建北部寒露风指数在1~5,其余大部晚稻区寒露风指数不足1。
图5 2020年双季晚稻区寒露风指数Fig.5 The index of chilling dew wind across double-season late rice region in 2020
双季晚稻区寒露风指数2020年为4.21,为2000年以来第二高,仅次于2011年的5.16(图6)。从省(自治区)尺度寒露风指数来看(图6),2020年湖南、江西、湖北、安徽、广西、广东、福建、浙江寒露风指数分别为7.16、7.16、12.59、3.83、0.16、0.00、1.06、1.72。湖南寒露风指数为2000年以来第二高,低于2011年的10.33;江西、湖北寒露风指数为2000年以来最高;安徽、广西、广东、福建、浙江寒露风指数分别为2000年以来第六、第五、第七、第五、第五高。
图6 2000—2020年各省(自治区)和双季晚稻区寒露风指数Fig.6 The integrated indexes of chilling dew wind across double-season late rice region and provinces during 2000-2020
利用土地利用/土地覆被遥感监测数据提取下垫面信息,在此基础上构建点-省-晚稻区的寒露风指数,分析了2020年南方双季晚稻区寒露风特点,主要结论如下:
(1)发生时间早、持续时间长:湖北南部、安徽南部、江苏南部、江西、湖南、浙江和广西北部等地寒露风天气出现日期普遍偏早,寒露风持续天数在3 d以上,大部出现轻度至重度寒露风。
(2)影响范围大:双季晚稻区寒露风发生面积比例为24.7%,是2000年以来第三高,轻度、中度、重度寒露风面积比例分别为9.7%、12.8%和4.9%;在各省(自治区)中,湖南、江西、浙江发生面积比例较大,分别为99.9%、99.4%、91.8%,其余省(自治区)不足30.0%。湖南、湖北、安徽、江苏寒露风发生面积比例为2000年以来最高;江西寒露风发生面积比例为2000年以来第二高,重度寒露风面积比例为2000年以来最高;浙江、福建、广西寒露风发生面积比例分别为2000年以来的第四、第五、第三高,以轻度为主;广东寒露风发生面积较小。
(3)强度高:湖北南部、安徽西南部、江西北部、湖南中西部、浙江中部寒露风指数较高,普遍在5以上。双季晚稻区寒露风指数为4.21,为2000年以来第二高;湖南、江西、湖北寒露风指数分别为 7.16、7.16、12.59,其余省(自治区)不足4.00;湖南寒露风指数为2000年以来第二高,江西、湖北为2000年以来最高。
本文通过构建寒露风指数,分析寒露风时空多尺度发生特点。对比历史寒露风发生特点,1981年以来寒露风指数整体呈降低趋势,与1981—2010年相比(平均指数为1.53),近十年寒露风强度减弱(平均指数为1.39);但寒露风指数高值出现频率并没有降低,近十年中寒露风指数较高的年份出现了两年(2011年和2020年为历史前五位)。可见,气候变暖背景下强寒露风的发生并没有减弱,这与已有研究结果较为一致(刘丹等,2019)。在这种背景下,寒露风指数较高的年份(2011年和2020年),晚稻一旦遭受寒露风天气会受到明显影响,正如田间调查结果显示:2011年湖南南部、江西南部、浙江南部、福建西部和北部以及两广北部局部等地晚稻遭受轻至中度寒露风危害,晚稻无法正常抽穗开花和授粉、结实率降低,部分处于抽穗扬花始期的晚稻空壳率高达60%~70%,处于抽穗扬花普遍期的空壳率为40%~50%,处于抽穗末期的空壳率为20%~30%;2020年湖南中北部、江西北部等地区晚稻遭受寒露风天气影响较重,出现抽穗开花缓慢、花粉败育、授粉不良、结实率下降、空壳率上升、灌浆不充分,后期千粒重偏低,影响晚稻产量形成。此外,胡磊等(2020)研究表明2006年属于中度寒露风,这也与计算的寒露风指数属于中等水平相一致。当然,本研究建立的寒露风指数虽然能反映出寒露风年际间的差异,但是要量化寒露风的具体影响程度(如对产量的影响),还需要通过进一步获取有效的试验数据样本从而构建综合性的寒露风影响指数(如低温积寒等)。此外,本研究在数据方法上也存在不确定性。首先,受土地利用数据的限制,计算过程中采用了统一的分类数据,即粗略地认为近年来水田面积一致,忽略了年际间的差异,这与实际水田面积存在一定的差异。其次,在土地利用数据分类解译时把水田粗略地认为是晚稻田,这里也必然存在一定的误差,因为水田并不等于晚稻田,下一步将尝试利用遥感反演和地面实测结合的方式确定晚稻田的具体分布情况,从而从水田中进一步分离出晚稻田,提高精细化服务水平。
根据寒露风天气特点,尤其是区域发生特点,可以为采取区别化、针对性的生产措施调整和防御提供参考。在产前,根据区域寒露风发生和风险特点,尤其是在江西、湖南等寒露风风险较高的区域,需要选择合适的能在寒露风来临前大概率完成齐穗的品种,在充分考虑品种生育特性的基础上尽可能早播早插(王华等,2011)。在产中,根据实时天气和晚稻生长发育进程,对寒露风进行提前预测(罗伯良和李易芝,2015;吴立等,2016a;2016b),从而采取相应的管理技术(李超等,2018;苏荣瑞等,2012),如在寒露风到来时,对于尚未抽穗晚稻采取灌深水增温,对已进入始穗的晚稻采取结合叶面适量喷施磷酸二氢钾加调节剂等促进晚稻早抽穗;对正处于抽穗扬花的晚稻应灌深水保温御寒、喷施增温剂以减轻低温危害确保安全齐穗。寒露风过后,受影响水稻对病虫害的抵抗力下降,需预防白叶枯病和细菌性条斑病以及稻飞虱等危害,避免或减轻病虫害对于水稻产量的影响。