基于Aqua Crop 的南宁市甘蔗响应气象干旱情景模拟*

杨云川 ,程禹灏 ,梁丽青 ,廖丽萍** ,王婷艳 ,张会娅 ,杨星星,胡甲秋

(1. 广西大学土木建筑工程学院 南宁 530004;2. 广西岩溶区水安全与智慧调控工程研究中心/工程防灾与结构安全教育部重点实验室/广西防灾减灾与工程安全重点实验室 南宁 530004)

甘蔗(Saccharum officinarumL.)作为最具潜力的高产生物能源经济作物[1-2],是制糖的重要原料,其生长对光照、雨热需求较高,导致全球甘蔗种植主要分布在坡度5%～30%、纬度30°N～30°S 的地区[3]。中国是仅次于巴西和印度的世界第三大甘蔗生产国,而广西则是中国甘蔗种植和产糖量第一大省[4-5]。水分供给是高耗水甘蔗正常生长的一个主要控制因子[3-6],但广西历年降水季节分配不均,蔗区多石山丘陵,坡耕地居多,岩溶发育及土层保水能力弱,且全区甘蔗种植仍以雨养为主,致使甘蔗饱受干旱灾害侵袭而产量损失巨大。南宁市是广西甘蔗的三大种植产区之一,更是广西年降雨量最少的市区之一,区域气象干旱灾害频发,并具有多重时间尺度、旱涝并存及交替叠加的现象,对甘蔗生长及产量累积等造成显著影响[7]。为此,开展南宁市甘蔗生长过程、生物量及产量累积对各种气象干旱情景的响应机制研究,在明晰广西甘蔗旱灾响应机理,实现甘蔗旱灾风险调控及预警、田间精准智慧管理及制定有效的防灾减灾应对措施等方面具有重要科学意义。

近年来,关于甘蔗干旱的研究主要是分析甘蔗干旱时空演变特征[8-9]、基于旱灾损失资料的风险评估[10-11]、水分亏缺胁迫及抗旱的生理生化性能指标响应、干旱对其形态特征的影响等方面[12-13],而借助反映生长物理过程的作物模型开展甘蔗对各种干旱(水分亏缺胁迫)情景的响应机制研究则相对较少。作物模型基于系统科学思想,综合气候、土壤、作物生理及田间管理等因素,在揭示作物生长过程、生物量与产量累积对干旱的响应机制具有独特优势[14-15]。在众多作物模型中,联合国粮农组织(FAO)于2009年推出的Aqua Crop 作物生长模型[16-18]主要包括: 土壤水分平衡、作物生长模拟和大气组分等3 个模块,以水分为主控驱动因子,适于甘蔗这种高耗水作物的生长模拟研究。Aqua Crop 模型已被广泛应用于小麦(Triticum aestivumL.)、玉米(Zea maysL.)及水稻(Oryza sativaL.)等粮食作物的产量模拟及灌溉制度优化,证实其具有良好的模拟精度和应用潜力[19-20]。此外,针对广西气象干旱的季节性、骤发性并存及旱涝交替并发的特征[7,20-22],文中选择只需要降水资料的标准化加权平均降水指数(standard weighted average precipitation,SWAP)[23-25]进行南宁市气象干旱事件及其发生发展过程的逐日精细化识别,一方面可避免其他常用干旱指数因各气象要素资料不足或短缺而造成分析结果的不连续性;另一方面,可为结合作物模型实现甘蔗对气象干旱响应机制的业务化精确预警及旱灾风险智慧调控奠定科学与技术基础。

综上,本文基于对南宁市1978-2018 年逐日SWAP气象干旱时空特征统计分析,确定该区域甘蔗各生育期实际可能发生的气象干旱情景,在本地化Aqua Crop 模型参数基础上,实现甘蔗生长、生物量及产量累积对气象干旱的强度、历时及生育期敏感性等关键要素变化响应的精细化模拟。

1 研究区域、数据与方法

1.1 研究区域

南宁市属亚热带季风湿润气候区,有“桂中腹地”之称,地貌类型以盆地和山地丘陵为主,周边区域岩溶发育(图1a),全年气候温和,日照充足,能很好地满足甘蔗生长的热量和光照需求,是广西三大甘蔗种植市区之一。全市降雨量充沛,多年平均降水量达1304.2 mm[26-27],但其时空分布不均,加之下垫面为复杂地形地貌,农业灌溉工程不足,造成区域多年来气象干旱频发,甘蔗干旱及旱灾损失巨大[10-11]。南宁市历年甘蔗种植面积相对稳定,2018 年为141 833 hm2,占全区的16.0%。其空间分布呈总体分散、局部聚集的总体格局,密集种植区主要有武鸣区、横县、宾阳县、邕宁区以及江南区等(图1b)。多年来,南宁市甘蔗种植仍以雨养为主,气象干旱造成的土壤水分亏缺一直是影响该区域甘蔗生长及其产量的一个主要因子。

图1 南宁市数字高程模型(DEM,a)、甘蔗种植分布(2018 年)及China Meteorological Forces Dataset (CMFD) 197 个数据格点(b)Fig.1 Nanning Digital Elevation Model (DEM,a),sugarcane cultivation distribution (2018) and China Meteorological Forces Dataset (CMFD) 197 data grid points (b)

1.2 研究数据

研究中涉及气象、土壤、作物、田间试验及管理等各类数据集。其中,气象数据采用逐日尺度中国气象强迫数据集(China Meteorological Forces Dataset,CMFD),是我国学者Jie 等[28]开发的0.1°空间分辨率网格化数据,包括1979-2018 年逐日降雨量、最高气温、最低气温、平均气温、相对湿度、风速及太阳辐射量等;其通过遥感产品、再分析数据集和原位站数据融合而成,具有连续的时间覆盖和良好质量,是目前应用最广泛的气候气象数据集之一(http://data.tpdc.ac.cn/en/data/8028b944-daaa-4511-8769-965612652c49/)。本研究中,通过与地面气象站点数据对比分析,发现CFMD 日降雨、日气温与气象站数据相关性分别达0.80 和0.98 以上,总体精度良好。文中采用的土壤数据,主要是中国土壤数据库(http://vdb3.soil.csdb.cn/),并以寒区旱区科学数据中心《基于世界土壤数据库的中国土壤数据集》作为参考;具体通过Soil Water Characteristics 计算获得的站点土壤水力参数;根据Aqua Crop 模型提供的作物信息库确定甘蔗生长模拟的作物参数,并通过敏感性分析和率定参数实现本地化;验证数据采用广西大学农学院农场的甘蔗全生育期田间试验观测数据;田间管理数据按实际田间种植情况设置。

1.3 研究方法

1.3.1 标准化加权平均降雨指数(SWAP)

为描述区域旱涝的逐日变化,Lu[29]于2009 年基于降水要素,并概化考虑下垫面水分耗散的综合累积效应,提出加权平均降水指数(Weighted Average Precipitation,WAP),计算公式如下:

式中:n=0 代表当前日;N是当前日的前期影响天数;Pn为 前期第n日降水量;a为表征前一日的降水对当前日下垫面湿润程度的贡献强度经验值参数,论证时取0.9 为适宜值,此时N取44 d。此后,Lu 等[30]进一步针对WAP,提出Gamma 函数正态标准化的SWAP指数,计算如下:

SWAP 属于正态标准化变量,具有与标准化降水指数(Standardized Precipitation Index,SPI)相同的强度等级划分标准,即按正态分布曲线的数学拐点特征来划分干旱等级。但近年研究发现以出现概率(重现期)为标准来划分干旱等级更加科学合理,并且与实际旱情更吻合。因此,本文采用以出现概率为标准的干旱等级划分(表1)。

表1 标准化的加权平均降水指数(SWAP)干旱强度等级表Table 1 Drought intensity gradation according to Standard Weighted Average Precipitation (SWAP)

基于上述等级划分阈值,采用游程理论[25,31]识别并统计分析研究区的气象干旱事件特征(结果详见2.1 节)。具体设定3 个截断水平X0=-0.84、X1=-0.52和X2=0 (X表示SWAP 干旱指数值),识别步骤如下:1)当一定时间尺度下(5 d 及以上天数) SWAP 指数值小于X1(-0.52)时,则初步判定该时间段发生干旱;2)对于历时为1 个时间段(5 d)的干旱,若其干旱指数值大于X0,则认为该时间段未发生干旱;3)对于时间间隔仅为1 d 的相邻2 次干旱过程,若间隔期的干旱指数小于X2,则这2 次干旱被视为从属干旱,将其合并为1 次干旱事件,否则视为2 次独立干旱过程。

1.3.2 Aqua Crop 模型原理

根据FAO 灌溉与排水第33 号文件所述的作物产量和水分响应的转换关系如下[32]:

式中:Yx和Y0分别为作物的潜在产量和实际产量,kg·m-2;E Tx和ET0分别为作物潜在蒸散量和实际蒸散量,mm;ky为产量对水分响应的系数。

Aqua Crop 模型对上述方程进行了改进,将蒸散量进一步分为土壤蒸发量和作物蒸腾量两部分,从而避免了非生产性用水(土壤蒸发)与生产性用水(作物蒸腾)效应的混淆;最终的产量以生物量和收获指数来表示,用以突出水分胁迫对二者各自的影响。改进后的公式如下:

式中:Y为最终作物产量(t·hm-2);B为生物量(t·hm-2);HI为收获指数;WP 为生物量水分生产效率(kg·m-2·mm-1);Tr为作物蒸腾量(mm);H I0是 参考收获指数;fHI是调整系数,用来反映各种胁迫(水分胁迫、温度胁迫等)对作物产量的影响。

此外,文中涉及甘蔗各要素受旱的减小量描述,分别采用产量减少率、生物量减少率、作物总蒸腾减少率来表达,公式如下:

式中,Yw、Bw、Trw分别表示甘蔗因干旱水分胁迫的减产率、生物量减少率、总蒸腾减少率,Yi、Bi、Tri分别表示干旱情景下的产量、生物量、作物总蒸腾量,Yck、Bck、Trck表示对照组2017 年的对应要素。

1.3.3 甘蔗干旱过程模拟数据库构建

Aqua Crop 模型输入参数主要包括研究区的气象、土壤、作物生长及管理措施等4 类信息。气象数据库包括降水量(P)、参考蒸散发量(ET0)、日最高/最低气温(Tmax/Tmin)和CO2浓度4 类,ET0可通过FAO 发布的ET0 Calculator 求出。土壤数据库包括土层剖面水分和地下水两个模块,土层剖面水分参数主要有田间持水量、永久凋萎点、饱和导水率及饱和渗透系数等,地下水参数主要有地下水位和盐度等;具体采用模型中的Soil Water Characteristics 工具计算土层剖面水分参数,地下水参数影响较小,按默认值处理。作物数据库包括植物冠层、植物根系、植物蒸腾、产量、水分、盐分、肥力及温度胁迫等模块,且每个模块均包括多个参数。本研究针对甘蔗作物生长特征,综合参考Bahmani 等[33]、Zu 等[34]、阮红燕[35]学者对不同地区甘蔗作物对水分亏缺响应、潜在生产力及气候变化影响研究中所采用的作物参数值,形成本文中的甘蔗作物参数库,进而通过广西大学农学院农场开展的甘蔗田间试验观测数据(表2)进行作物参数敏感性分析及率定验证,实现Aqua Crop 模型模拟参数的本地化。管理数据库包括灌溉管理和田间管理两个模块,灌溉管理有灌溉方式、灌溉时间及灌溉量等,田间管理有施肥水平、覆盖程度及田间地表措施等。具体按表2 中的田间试验的管理措施进行设置。

表2 率定Aqua Crop 模型所用试验数据来源Table 2 Sources of experimental data for calibration of Aqua Crop model

本研究主要输出南宁市甘蔗生长过程中的冠层覆盖度(canopy cover,CC)、作物蒸腾(crop transpiration,Tr)、生物量(biomass production,B)及甘蔗产量(cane yield,Y)等4 个要素。

2 南宁气象干旱特征及甘蔗生育期干旱情景设置

2.1 南宁市气象干旱统计特征

基于南宁市甘蔗生育期时段划分(3 月10 日-12 月25 日,共291 d,详见表3),采用游程理论[25,31]识别并统计了该市1980-2018 年(197 个格点)逐日SWAP 序列表达的气象干旱;然后根据Matlab 程序计算出干旱事件的干旱强度、干旱次数和干旱频次。甘蔗全生育期时段内干旱累积历时、平均强度和平均频次空间分布如图2a-c 所示,甘蔗各生育期的历年干旱累积历时和强度变化特征如图2d-e 所示。图2a 表明,南宁市干旱历时空间分布不均匀,但总体主要在90 d·a-1以上,尤其在市区北部及横县区域达105～120 d·a-1,即南宁市平均每年约有近1/3 的时段为干旱天数。图2b 显示,南宁市所发生的气象干旱事件强度总体呈中部高南北局部低的空间分布格局,大部分区域的干旱强度等级为中旱和重旱。由图2c可知,南宁市气象干旱事件发生的频次呈中北部及横县区域高(2～3 次·a-1),其他局部区域相对较低(1～2次·a-1,如市区南部、宾阳县局部等)的分布格局。综上可知,南宁市基本每年都有气象干旱事件发生,且强度以中旱和重旱为主,年度累计干旱历时天数常在90 d 以上。由图2d、e 可知,近40 年来,南宁市气象干旱发生在甘蔗各生育期的历时和强度差异显著,总体上,成熟期＞茎伸长期＞萌芽期＞分蘖期。其中,萌芽期历年发生的气象干旱历时变化范围主要为5～40 d、强度变化范围为0～-70;茎伸长期分别为5～40 d、0～-80;成熟期分别为5～50 d、0～-70;而分蘖期时段短(40 d)且为历年雨季,基本不发生气象干旱。南宁市甘蔗全生育期历年发生的气象干旱总历时与累积强度范围分别为10～120 d 和0～-152,大部分年份干旱天数超过60 d;但目前关于气象干旱对南宁市以雨养为主的甘蔗生长造成的定量影响研究成果相对不足,因此,要制定有效的甘蔗旱灾综合防范方案,对应的甘蔗生长过程响应机制的定量解决亟待澄清。

图2 1980—2018 年南宁市气象干旱历时(a)、强度(b)和频次(c)的空间分布及历时(d)和强度(e)的年际变化Fig.2 Spatial distribution of meteorological drought duration (a),intensity (b) and frequency (c),and intrannual variations of meteorological drought duration (d) and intensity (e) in Nanning from 1980 to 2018

2.2 甘蔗生育期气象干旱情景设置

根据南宁甘蔗种植及生长成熟实际情况,将3月10 日-12 月25 日(共291 d)确定为全生育期时段,各生育期按表3 的时间进行划分。为了更好地揭示不同历时和强度气象干旱情景下,南宁市甘蔗的长势、生物量及产量累积响应机制,结合2.1 节的分析结果,将历年甘蔗气象干旱的可能情景设置如表3 所示,干旱历时变化步长设定为5 d,历年没有发生的气象干旱情景则不进行模拟分析。

表3 甘蔗各生育期历时和干旱模拟情景设置Table 3 Duration and scenario setting of different drought grades of different sugarcane growth stages

3 Aqua Crop 模型参数化及南宁市甘蔗的历史干旱响应模拟

3.1 Aqua Crop 模型参数敏感性分析

以南宁市甘蔗种植的主要品种‘新台糖16 号’进行模型主要敏感参数的分析及本地化率定。采用扩展傅里叶幅度检验法(EFAST)对Aqua Crop 模型的13 个主要参数(表4)进行全局敏感性分析[32-36],以便实现模型参数的快速本地化。分析中,各参数的变化范围设置为参考值的上下限加减30%,参数采样方法为蒙特卡罗法,采样次数3965 次。模拟分析了南宁广西大学农场试验站(表2) 2015 年甘蔗的生物量与产量响应各参数的敏感性。最终得到各参数的一阶敏感性指数(sensitivity index,Si)和全局敏感性指数(total order sensitivity index,STi) (图3),前者表示参数对模型输出的直接影响,后者则表示参数对模型输出直接和间接影响的加和。参考Dejonge 等[37]对全球农业生态系统的敏感性阈值,将Si＞0.05、STi＞0.1 的参数视为较敏感参数。上述分析过程主要通过Simlab (V2.2)和Python 工具实现。

由图3 可知,Aqua Crop 模型中的1～4、6、7、9及11 号(表4)参数敏感性均较弱,其数值变化对甘蔗产量和生物量没有显著影响。而5、10 及12 号参数的Si和STi值均显著高于敏感性指数判别阈值,其数值变化对甘蔗产量和生物量均具有较显著的影响。8 号参数Si值略高于一阶敏感性判别阈值,但低于全局敏感性判别阈值,总体敏感性有限。而13 号参数的Si和STi值显示,其变化对甘蔗产量有显著影响,而对甘蔗生物量没有影响。

图3 Aqua Crop 模型主要参数对甘蔗产量和生物量的敏感性判别Fig.3 Parameters sensitivities to sugarcane yield and biomass by Aqua Crop model

表4 Aqua Crop 模型进行敏感性分析的参数及其取值范围Table 4 Parameters and ranges of Aqua Crop model involved in sensitivity analysis

3.2 Aqua Crop 模型参数本地化

基于参数的敏感性分析结果,重点对Aqua Crop模型的敏感性参数及作物生长控制参数等进行多次调整,最终率定得到了Aqua Crop 模型模拟南宁甘蔗生长的本地化参数取值(表5)。进而对比Aqua Crop模型模拟与试验站17 个田块实测的甘蔗产量可知(图4),二者拟合R2和精确度(Pre)分别为0.92 和0.89,模拟产量的均方根误差(RMSE)及其误差百分率分别为4.05 t·hm-2和3.84%,总体模拟效果优良。由此表明,本研究中采用Aqua Crop 模型及其本地化参数开展南宁甘蔗生长过程模拟具有较高的模拟精度和应用价值。

图4 Aqua Crop 模型模拟甘蔗产量结果验证Fig.4 Verification of sugarcane yield simulation by Aqua Crop Model

表5 南宁市甘蔗生长模拟的Aqua Crop 模型本地化参数取值Table 5 Sugarcane crop localized parameters of sugarcane crop of Aqua Crop model in Nanning

3.3 南宁市历史典型年气象干旱的甘蔗响应模拟

为厘清南宁市历史SWAP 气象干旱情景下甘蔗长势、生物量及产量累积响应特征,在1978-2018年时段内依次选择典型无旱年2017 年(作为参考对照年份,CK)、典型轻旱年2016 年、典型中旱年2009 年、典型重旱年1992 年,分别对甘蔗的蒸腾量(Tr)、冠层覆盖度(CC)、生物量(B)及产量(Y)等要素进行模拟,结果如图5 所示。

从图5 可知,南宁市甘蔗Tr在分蘖期和伸长期对不同气象干旱情景具有敏感的波动响应变化特征,如该生育期时段内各典型年SWAP 与Tr序列的相关系数分别为0.38 (2017 年)、0.66 (2016 年)、0.69(2009 年)及0.70 (1992 年),均通过显著检验;Tr在伸长期中后时段受气象干旱影响呈显著下降变化,尤其在中旱(2009 年)和重旱(1992 年)年份显著小于2017 年无旱年份数值。而在萌芽期和成熟期,甘蔗Tr总体均很小(接近0 值),与气象干旱SWAP 序列无明显响应关系。

由图5 可知,甘蔗CC 在全生育期的变化主要受生理生长过程控制,即自萌芽期中后时段开始快速增大到分蘖期末的最大值90%～100%,并在伸长期能维持在最大值状态(如2017 年全时段及其他典型年的无旱时段);气象干旱对甘蔗CC 能产生一定影响,但需要水分亏缺累积到一定量值才能导致甘蔗CC的显著减小(如1992 年重旱),轻旱和中旱的影响相对较小(如2016 年和2009 年);而成熟期甘蔗CC 快速减小至0 主要是生理现象,基本不受气象干旱的影响。

图5 显示,各典型年甘蔗B均从分蘖期初时段的0 值快速增加到伸长期末时段的最大值(如2017年为66.535 t·hm-2),成熟期的增量最大约为10%;表明气象干旱能一定程度上影响甘蔗B累积的最终数量(如2016 年、2009 年、1992 年分别为58 t·hm-2、53 t·hm-2、45 t·hm-2,相 对2017 年分别减小13%、20%、33%),而不能影响其全生育期单调增加的总体规律。

图5 南宁市历时典型年甘蔗干旱响应模拟Fig.5 Simulation of sugarcane drought response in typical years in Nanning City

图5 表明,甘蔗茎产量主要在伸长期中期(约8月中旬)开始迅速累积至成熟期末达最大值,在无旱的2017 年产量达116 t·hm-2,而轻、中、重旱年分别为97 t·hm-2(2016 年)、87 t·hm-2(2009 年)及70 t·hm-2(1992 年),对应减产率分别为16%、25%及40%。由此可知,不同强度等级的气象干旱均对甘蔗产量具有显著的影响。

4 基于气象干旱情景的甘蔗响应机制模拟

3.3 节的模拟结果表明气象干旱对南宁市甘蔗的生长、生物量及产量累积等均具有显著影响,本节基于2 节中气象干旱及其在甘蔗生育期的实际可能发生情景,进一步解析各强度及历时的气象干旱发生在不同生育期对甘蔗蒸腾量(Tr)、生物量(B)及产量(Y)累积的影响机制。

4.1 不同生育期同气象干旱情景的甘蔗响应机制模拟

本节主要分析了不同生育期发生同气象干旱(强度和历时)对南宁甘蔗各生长要素的影响机制。由图6a-c 可知,甘蔗各生育期发生轻旱时,在萌芽期,甘蔗的Tr基本不受各历时轻旱的影响;甘蔗的B和Y则在轻旱历时为5～15 d 时基本不受影响,而在轻旱历时为20～40 d 情景下呈逐步减小变化,Y最终减少到110 t·hm-2(相对2017 年,减产率5%);在伸长期,甘蔗Tr、B和Y均受10 d 及以上历时的轻旱影响而显著减小,各要素最终减小到815 mm、54 t·hm-2和88 t·hm-2(减产率24%);而在成熟期,甘蔗各要素基本没有显著变化。

图6d-f 表明,甘蔗各生育期发生中旱时,在萌芽期和伸长期,甘蔗Tr、B和Y均随中旱历时的增加而呈显著减小的变化特征,其中,萌芽期三要素最小值分别为902 mm、52 t·hm-2和87 t·hm-2(减产率25%),伸长期三要素最小值分别为734 mm、47 t·hm-2和73 t·hm-2(减产率37%)。而在成熟期,各要素基本不受中旱的影响。图6g-h 显示,甘蔗各生育期发生重旱时,在伸长期,甘蔗Tr、B和Y均随重旱历时增加而显著减小至715 mm、45 t·hm-2和70 t·hm-2(减产 率40%),而在成熟期则无显著变化。

图6 不同生育期发生同一干旱强度的甘蔗产量响应情景模拟Fig.6 Scenario simulation of sugarcane yield response to the same drought intensity in different growth stages

4.2 同生育期不同气象干旱情景的甘蔗响应机制模拟

本节进一步定量化了同生育期内发生不同强度及历时气象干旱对甘蔗各生长要素的影响。由图7ab 可知,萌芽期发生轻旱时,在干旱历时为5～15 d 时,甘蔗Yw、Bw、Trw的减少量约为0 且维持不变,随干旱历时逐渐增加至15～40 d 时,三要素减少量快速减少,各自减产率最终分别达到5.9%、5.1%、6.5%;而萌芽期发生中旱时,随干旱历时的增加(15～35 d),甘蔗的Yw和Bw分别从4.0%、3.0%分别快速减少至26.8%和21.9%,Trw则总体减少量不显著。图7c-e表明,伸长期发生轻、中、重旱时,随气象干旱历时的增加,甘蔗Yw、Bw、Trw的减少量均呈显著增加的变化特征。其中,各强度气象干旱造成Yw、Bw、Trw的减少量变化范围分别为0～24%、0～18.5%及0～15.9% (轻旱历时5～35 d),25%～37%、20%～29.3%及8%～24.4% (中旱历时15～45 d),33.5%～40%、26.2%～31.7%及18.9%～25.7% (重旱历时35～50 d)。

图7 同一生育期不同强度干旱的甘蔗产量、生物量及作物总蒸腾的变化率Fig.7 Changes of sugarcane yield,biomass and crop transpiration under different drought intensities in a growth period

5 结论

本文通过分析南宁市1980-2018 年逐日SWAP气象干旱特征及在甘蔗生育期的发生情景,本地化Aqua Crop 作物模型参数,实现了该区域甘蔗生长、生物量及产量累积过程对各强度及历时气象干旱的响应机制模拟。论文研究得出以下主要结论:

1)南宁市气象干旱历时、强度及频次空间分布不均匀,年干旱历时主要在90 d 以上,强度以中旱及重旱为主,发生频次以2～3 次·a-1居多,且主要发生在甘蔗的萌芽期、伸长期和成熟期,存在长(月以上)、短(月内)历时干旱交替及并存叠加现象。

2)通过EFAST 敏感性分析及参数本地化率定,Aqua Crop 模型模拟甘蔗产量拟合精度R2达到0.92、产量均方根误差百分率为3.84%。Aqua Crop 模型应用于南宁甘蔗对气象干旱的响应机制模拟能得到良好的模拟精度。

3)历时典型干旱年的甘蔗响应模拟揭示了甘蔗蒸腾量、冠层盖度、生物量及产量等对气象干旱的响应存在滞后和累积效应,伴随干旱强度增大,南宁甘蔗最终累积的生物量及产量减少率分别为12.7%(2016 年轻旱)、20.1% (2009 年中旱)、32.7% (1992年重旱)和16.6% (轻旱)、24.8% (中旱)及40.1%(重旱)。

4)各气象干旱强度及历时遍历甘蔗各生育期的情景模拟,明晰了南宁甘蔗生长过程、生物量及产量累积对气象干旱的响应机制。萌芽期发生轻、中旱历时达到15 d 及以上时对甘蔗各生长要素开始产生显著响应,最终减产率分别达6%和27%;伸长期发生轻、中和重旱历时为5 d 及以上时甘蔗各要素即开始显著响应,最终减产率分别达24%、37%和40%,南宁甘蔗分蘖期基本不发生气象干旱,而成熟期各要素受气象干旱的影响总体微弱。