新疆阿克苏地区光照充足、昼夜温差大、高山雪水灌溉等独特的优势条件非常适宜苹果的生长
[1]。阿克苏苹果独特的“冰糖心”品质,以及色泽鲜红、脆爽多汁、果肉细腻等特点,深受消费者喜爱
[2]。阿克苏作为新疆苹果的核心产区之一
[3],截至2023年,阿克苏地区苹果种植面积已达3.37万hm
2,总产量为81.72 万t,苹果产业已成为当地的主要经济支柱之一。阿克苏地区苹果园灌溉常采用漫灌方式,并伴随高量肥料投入,这导致水肥利用效率较低。随着人们生活水平的不断提高,消费者对苹果品质的要求也越来越高
[4]。因此,为促进阿克苏地区苹果的优质高效生产,亟须制定科学合理的灌溉施肥制度。
研究表明,合理的水肥处理可以有效促进苹果树的生长,提升果实的产量和品质,同时提升土壤肥力
[5]。张军等
[6]的研究表明,水肥一体化技术可以提高番茄的产量,适宜的水肥配比可以提高番茄的品质。白少倩
[7]则证实,不同水肥组合处理骏枣的生长、产量和品质有显著差异,而水肥用量超过一定限值后会导致水分利用效率下降。李青青等
[8]的研究指出,水肥耦合不仅能显著提高葡萄的水分和肥料利用效率,还能够改善土壤性质,从而实现葡萄高产优质的目标。因此,合理的灌溉和施肥是提升作物产量和品质的关键。
灌溉制度优化通常分为单目标优化与多目标优化,大多以作物产量最大化或综合效益最优为目标。王凯航
[9]采用AquaCrop-FloPy和第三代非支配排序遗传算法(non-dominated sorted genetic algorithm-Ⅲ,NSGA-Ⅲ)耦合模型对玉米多目标灌溉制度进行了优化。康燕楠等
[10]通过NSGA-Ⅲ算法与FloPy模型构建了灌区水资源多目标优化模型,实现了灌区水资源多目标的优化配置。宋健等
[11]基于AquaCrop模型、NSGA-Ⅱ算法以及功效系数法,构建了适应多决策需求的灌溉制度多目标优化模型,较好地满足了不同决策情境的需求。
综上可知,水肥调控对作物的产量和品质具有积极作用,但在干旱区滴灌条件下,苹果生长的水肥耦合效应及水肥制度的优化仍需待进一步研究。本研究以阿克苏地区苹果为对象,探讨滴灌条件下水肥耦合对苹果产量、品质、水肥利用效率的影响,并通过构建包含苹果产量、灌溉水利用效率、肥料偏生产力、可溶性固形物含量和果肉硬度在内的5目标二元二次回归方程,结合NSGA-Ⅲ算法进行验证与求解,并基于熵权-Topsis法从中选择最优解,确定最优灌水量和施肥量,为干旱区滴灌苹果的水肥管理以及提质增效提供理论依据。
1 材料与方法
1.1 试验区概况
试验于2024年3月至10月在新疆农业大学国家良种基地(80°14' E,41°16' N,海拔1 133 m)进行。试验地点位于阿克苏市红旗坡十队,地处温带大陆性气候区,平均年日照时数为2 855~2 967 h,多年平均气温为11.2 ℃,极端高温高达40.9℃,极端低温达到-27.4 ℃,无霜期205~219 d,多年平均降水量为29.6~71.4 mm。
供试材料为树龄8年的苹果树,品种为‘新红一号’,种植株行距为3 m×4 m。试验地块平整,土壤质地以壤质土为主,苹果树长势一致。田间管理按照当地果园的常规操作进行,包括疏花疏果、果实套袋与摘袋、铺设反光膜以及病虫害防治等措施,确保试验条件与实际生产状况相符。
1.2 试验方案
试验设置3个灌水梯度(675,800和925 mm)及3个施肥梯度(562.5,750.0,937.5 kg/hm
2),进行完全组合设计,共9个处理,每个处理设置3次重复,每个试验小区面积60 m
2,每个小区5株苹果树,具体的灌水与施肥方案见
表1,其中3月20日进行春灌;4月22日起每10 d为1个周期进行生育期灌水,生育期灌水截至9月12日,共计灌水15次;11月20日进行冬灌。灌溉采用滴灌技术,使用内镶贴片式滴灌带,滴头流量2 L/h,滴头间距为0.3 m;每行苹果树两侧各布置一条滴灌带,距离树干0.3 m,灌水量采用水表精确计量。施肥采用水肥一体化技术,每个处理配置单独施肥罐,实现“肥随水移”,试验中所用肥料包括氮肥(尿素,N含量为46%)、磷肥(磷酸一铵,N、P
2O
5含量分别为14%,46%)、钾肥(硫酸钾,K
2O含量为52%)。
1.3 测定指标与方法
1.3.1 果实纵径、横径和果形指数
根据苹果的生长特点,将生育期划分为开花坐果期(4月15日-5月15日)、幼果期(5月16日-6月15日)、果实膨大期(6月16日-9月10日)、成熟期(9月11日-10月20日)4个阶段。开花坐果20 d后,每个处理选择3株长势相近、无病虫害的苹果树进行标记,从每株树的东、西、南、北4个方向各选取1个果实标记,并使用游标卡尺每间隔10 d测量1次果实纵径和横径,并计算果形指数(最大纵径/最大横径)。最后统计幼果期、果实膨大期、成熟期和全生育期苹果纵径、横径以及果形指数的平均值。
1.3.2 产量和品质
在每个处理中选择3株长势相近的苹果树,将所有果实采摘后称质量并计数,折算为每公顷产量。测定苹果品质时,从每个处理中选择3株树,在同一高度从东、西、南、北4个方向各选择1个果实进行分析。采用电子秤测定单果质量,用手持折射仪测定可溶性固形物含量,用GY-4型硬度计测定果肉硬度。
1.3.3 水肥利用效率
灌溉水利用效率及肥料偏生产力的计算公式为:
式中:IWUE为灌溉水利用效率(kg/(hm2·mm)),Y为总产量(kg/hm2),W为总灌水量(mm),PFP为肥料偏生产力(kg/kg),F为总施肥量(kg/hm2)。
1.4 基于NSGA-Ⅲ算法的苹果灌水量和施肥量优化
NSGA-Ⅲ算法是一种基于第二代非支配排序遗传算法(NSGA-Ⅱ)的进化多目标优化算法
[12]。NSGA-Ⅱ算法主要基于局部的种群密度进行选择,以每个个体为基准,考虑相对位置,由局部控制整体;而NSGA-Ⅲ算法则是通过分布均匀的参考点来控制种群的整体分布,以绝对均匀分布的参考点为基准,不受个体间相对位置的影响。NSGA-Ⅲ算法较NSGA-Ⅱ算法更适用于3个以上目标任务,同时计算效率更高
[13]。
为了提高数据的可对比性和精度,本研究对苹果总灌水量及总施肥量进行了标准化处理
[14],编码值见
表2。
本研究将产量、灌溉水利用效率、肥料偏生产力、可溶性固形物含量和果肉硬度设定为优化目标,将其与灌水量和施肥量的编码值进行二元二次回归方程拟合,建立回归模型,进而构建包含产量、灌溉水利用效率、肥料偏生产力、可溶性固形物含量和果肉硬度的多目标优化模型
[15],并利用NSGA-Ⅲ算法进行优化分析。
1.5 数据处理
数据处理与分析采用Excel 2021、MATLAB 2022b和Origin 2024软件进行整理、分析及绘图,显著性检验和方差分析则使用SPSS 25.0完成。
2 结果与分析
2.1 不同水肥耦合处理对苹果果实纵径和横径的影响
不同水肥耦合处理对各生育期苹果果实纵径、横径的影响见
表3和
表4。
从
表3和
表4可以看出,在W1和W3处理下,果实纵径和横径均随着施肥量的增大而增大;在W2处理下,果实纵径和横径均随着施肥量的增加先增大后减小。在相同施肥量下,果实纵径、横径均随着灌水量的增加总体呈先增大后减小的趋势,总体表现为W2>W3>W1。这表明适宜的水肥组合更有利于苹果果实的生长。9个处理中, W2F2处理幼果期、果实膨大期、成熟期的纵径和横径均最高,分别为28.57,50.29,71.15 mm和30.15,58.68,85.59 mm。与W1F1处理相比,W2F2处理苹果的纵径分别提高了11.69%,9.54%,9.87%,横径分别提高了13.69%,7.33%,8.10%。
表3和
表4显示,灌水量对各生育期及全生育期苹果果实纵径、横径有极显著影响(
P<0.01);施肥量对幼果期和果实膨大期苹果果实纵径、横径有显著影响(
P<0.05),对成熟期和全生育期苹果果实纵径、横径有极显著影响(
P<0.01);灌水量与施肥量的交互作用对各生育期及全生育期苹果果实纵径、横径的影响均不显著(
P>0.05)。
2.2 不同水肥耦合处理对苹果产量的影响
图1显示,不同水肥耦合处理苹果产量有明显差异。W1与W3处理下,苹果产量表现为F3>F2>F1;而在W2处理下,苹果产量由大到小依次为F2>F3>F1,表明高肥处理下苹果产量略有下降。在F1、F3处理下,苹果产量由大到小依次为W3>W2>W1;而在F2处理下,苹果产量由大到小顺序为W2>W3>W1。总体来看,适宜的灌水量与施肥量组合能显著提高苹果产量。9个处理中,W2F2处理苹果产量最高,达到62.24 t/hm
2,较W1F1、W3F3处理分别提高了47.87%,2.74%。
2.3 不同水肥耦合处理对苹果品质的影响
不同水肥组合对苹果品质的影响如
表5所示。
表5显示,在W1与W3处理下,单果质量由大到小依次为F3>F2>F1;在W2处理下,则表现为F2>F3>F1。在F1与F3处理下,单果质量由大到小依次为W3>W2>W1;在F2处理下,则表现为W2>W3>W1。9个处理中,W2F2处理单果质量最大,为296.68 g;W1F1处理单果质量最小,为253.60 g。相比W1F1和W3F3处理,W2F2处理的单果质量分别提高了16.99%和1.79%。可见,适宜的灌水量与施肥量组合有助于单果质量的增加。各处理果形指数差异较小,其中W2F2处理最大(0.84),W1F1处理最小(0.78)。在W1与W3处理下,可溶性固形物含量由大到小依次为F3>F2>F1;而在W2处理下,则表现为F2>F3>F1。在F1与F2处理下,可溶性固形物含量由大到小依次为W2>W3>W1;而在F3处理下,则表现为W3>W2>W1。W2F2处理的可溶性固形物含量最高,为15.30%;W1F1处理最低,为11.21%。相比W1F1、W3F3处理,W2F2处理的可溶性固形物含量分别提高了36.49%和1.59%。同一灌水量下,果肉硬度由大到小依次为F3>F2>F1;同一施肥量下,果肉硬度由大到小依次为W3>W2>W1。9个处理中,W3F3处理的果肉硬度最高,为9.15 kg/cm
2;W1F1处理的果肉硬度最低,为8.31 kg/cm
2。相比W2F2、W1F1处理,W3F3处理的果肉硬度分别提高了2.12%和10.11%。
表5表明,灌水量对单果质量、可溶性固形物含量和果肉硬度有极显著影响(
P<0.01),但对果形指数无显著影响(
P>0.05);施肥量对可溶性固形物含量和果肉硬度有极显著影响(
P<0.01),对单果质量有显著影响(
P>0.05),但对果形指数无显著影响(
P>0.05);灌水量与施肥量的交互作用对可溶性固形物有极显著影响(
P<0.01),对果肉硬度有显著影响(
P<0.05),但对单果质量和果形指数无显著影响(
P>0.05)。
2.4 不同水肥耦合处理对苹果灌溉水利用效率及肥料偏生产力的影响
不同水肥耦合处理对水肥利用效率的影响见
图2。由
图2可知,在W1与W3处理下,灌溉水利用效率由大到小依次为F3>F2>F1;在W2处理下,则表现为F2>F3>F1。在F2与F3处理下,灌溉水利用效率由大到小依次为W1>W2>W3;在F1处理下,则表现为W2>W1>W3。9个处理中,灌溉水利用效率最高的为W1F3处理,达到79.39 kg/(hm
2·mm);最低的为W3F1处理,仅为59.74 kg/(hm
2·mm)。同一灌水量下,肥料偏生产力由大到小依次为F1>F2>F3。在F1与F3处理下,肥料偏生产力由大到小依次为W3>W2>W1;在F2处理下,则表现为W2>W3>W1。肥料偏生产力最高的是W3F1处理,达到98.24 kg/kg;最低的为W1F3处理,仅为57.16 kg/kg。
2.5 基于NSGA-Ⅲ算法的苹果灌水量和施肥量优化
本研究以灌水量编码值(
w)和施肥量编码值(
f)为变量,选择苹果的产量(
Z1)、灌溉水利用效率(
Z2)、肥料偏生产力(
Z3)、可溶性固形物含量(
Z4)、果肉硬度(
Z5)5项指标作为优化目标,构建二元二次回归方程,结果(
表6)显示,回归方程的
R2>0.85,可知方程拟合度较高。
2.5.1 单因子效应分析
为进一步研究灌水量和施肥量单因子对苹果产量、水肥利用效率、品质的影响,将
表6中的二元二次回归方程进行降维处理,得到的单因子效应函数见
表7,其对应的效应函数曲线见
图3。
从
图3可以看出,单因子效应函数曲线总体呈现为开口向下的抛物线,因此在特定区间内各目标存在最大值。就产量而言,灌水量编码值(
w)的影响为:当-1.202<
w<0.368,增加灌水量有助于提高产量;当
w=0.368时,产量达到最大值;当0.368<
w<1.202时,增加灌水量会导致产量下降。施肥量编码值(
f)对产量的影响为:当-1.202<
f<0.324时,提高施肥量有利于促进产量提高;当
f=0.324时,产量达到最大;当0.324<
f<1.202时,提高施肥量则会造成产量下降。
从
图3可以看出,就灌溉水利用效率而言,灌水量编码值(
w)的影响为:当-1.202<
w<-0.429时,灌溉水利用效率与灌水量编码值呈正相关;当
w=-0.429时,灌溉水利用效率达到最大值;当
-0.429<w<1.202时,灌溉水利用效率与灌水量编码值呈负相关。施肥量编码值(f)的影响为:当
-1.202<f<0.369时,灌溉水利用效率与施肥量编码值呈正相关;当f=0.369时,灌溉水利用效率达到最大值;当0.369<f<1.202时,灌溉水利用效率与施肥量编码值呈负相关。
从
图3可以看出,就肥料偏生产力而言,灌水量编码值(
w)的影响为:当-1.202<
w<0.426时,肥料偏生产力与灌水量编码值呈正相关;当
w=0.426时,肥料偏生产力达到最大;当0.426<
w<1.202时,肥料偏生产力与灌水量编码值呈负相关。施肥量编码值(
f)的影响为:当
f=-1.973时,肥料偏生产力达到最大值;当-1.202<
f<1.202时,肥料偏生产力与施肥量编码值呈负相关。
从
图3可以看出,就苹果果实可溶性固形物含量而言,灌水量编码值(
w)的影响为:当-1.202<
w<0.451时,增加灌水量有利于提高可溶性固形物含量;当
w=0.451时,可溶性固形物含量达到最大值;当0.451<
w<1.202时,增加灌水量则导致可溶性固形物含量下降。施肥量编码值(
f)的影响为:当-1.202<
f<0.443时,增加施肥量有利于提高可溶性固形物含量;当
f=0.443时,可溶性固形物含量达到最大值;当0.443<
f<1.202时,增加施肥量抑制可溶性固形物含量的提高。
从
图3可以看出,就苹果果肉硬度而言,灌水量编码值(
w)的影响为:当-1.202<
w<0.75时,增加灌水量有利于提高果肉硬度;当
w=0.75时,果肉硬度达到最大值;当0.75<
w<1.202时,增加灌水量抑制果肉硬度提高。施肥量编码值(
f)的影响为:当-1.202<
f<0.731时,增加施肥量会促进果肉硬度的提高;当
f=0.731时,果肉硬度达到最大值;当0.731<
f<1.202时,增加施肥量抑制果肉硬度的提高。
2.5.2 基于NSGA⁃Ⅲ算法的多目标优化模型分析
为实现提质增效,通过构建Z1、Z2、Z3、Z4、Z5二元二次方程组成基于NSGA-Ⅲ算法优化模型,寻求苹果产量、水肥利用效率及品质最优时的灌水量和施肥量为:
;
;
;
;
。
通过MATLAB 2022b对上述模型进行计算,种群大小设置为50,最大遗传代数设置为50,交叉概率设置为0.5,得到了Pareto解集;通过熵权-Topsis法对Pareto解集进行分析,计算得到不同指标的权重分别为:产量0.21、灌溉水利用效率0.21、肥料偏生产力0.19、可溶性固形物含量0.13、果肉硬度0.26;最终确定的w、f最优解分别为0.203和0.213,对应的灌水量和施肥量分别为821.11 mm和783.22 kg/hm2,对应的苹果产量、灌溉水利用效率、肥料偏生产力、可溶性固形物含量及果肉硬度分别为66.68 t/hm2、79.73 kg/(hm2·mm)、84.52 kg/kg、15.47%、9.16 kg/cm2,以上指标均优于 9个处理中表现较佳的W2F2处理。
3 讨 论
灌水量和施肥量是影响作物生长发育的2个关键因素
[16-17]。合理的灌水量与施肥量有利于促进作物的健康生长
[18]。研究显示,水肥耦合效应对蓝莓首年果实纵径、横径无显著影响;但对蓝莓次年果实纵径有极显著影响,而对果实横径无显著影响
[19]。本研究中,水肥耦合效应对苹果果实纵径、横径均无显著影响,这与前人研究结果不一致,可能与研究作物种类不同有关。
适宜的灌水施肥可以协调作物养分的分配,促进光合产物的合成
[20],从而提高作物产量,改善作物品质。有研究指出,合理的灌溉量与施肥量能显著增产,但增产效果会随着灌溉量和施肥量的增加而逐渐减弱
[21]。研究发现,在适宜的水肥条件下,库尔勒梨可溶性固形物和可溶性糖含量明显增加,可滴定酸和果核细胞含量降低,从而改善果实品质
[22]。本研究发现,9个处理中,中水中肥处理(W2F2)的产量最高,显著高于低水低肥处理(W1F1)和高水高肥处理(W3F3);此外,W2F2处理的单果质量、果形指数和可溶性固形物含量均达到最高。可知过高或过低的灌水量和施肥量不仅会降低苹果产量,而且会降低果实品质,这是由于缺水缺肥抑制了苹果树的正常生长,影响了光合产物的积累,导致果实养分供应不足从而降低了苹果产量和品质;而过量的水肥则导致苹果树营养生长过剩,削弱光合产物向果实的分配能力,从而影响苹果产量和品质。
研究表明,适量施用氮肥显著提高了辣椒的产量和水肥利用效率
[23]。在同一施肥量下,马铃薯水分利用效率表现为低水>中水>高水
[24]。在同一灌水量下,施肥量对棉花肥料偏生产力的影响显著,肥料偏生产力随着施肥量的增加呈下降趋势
[25]。本研究中,当灌水量相同时,F3和F2处理的肥料偏生产力显著高于F1处理,但灌溉水利用效率显著低于F1处理,表明较高的施肥量能明显提高灌溉水利用效率;W2F1和W3F1处理的肥料偏生产力均显著高于其他处理,显示适宜的水肥组合能明显提高作物的水分吸收能力和灌溉水利用效率。
前人通过建立回归模型发现,苹果产量、灌溉水利用效率和肥料偏生产力达到最大值时,对应的最优灌水量和施肥量差异较大,难以同时兼顾最高产量与最大水肥利用效率
[20]。MWIYA等
[26]研究表明,NSGA-Ⅲ算法可优化灌溉制度,在提高作物产量的同时提高水资源利用效率。ZHAO等
[27]则通过NSGA-Ⅲ算法的优化模型,成功解决了枸杞品质与产量之间的最优灌水量与施肥量差异较大的问题,获得了自动滴灌的最优灌水量与施肥量。本研究构建了二元二次回归方程,并基于NSGA-Ⅲ算法得到了苹果产量、水肥利用效率及品质最优时的灌水量为821.11 mm,施肥量为783.22 kg/hm
2,其对应的苹果产量、水肥利用效率及品质均优于W2F2处理,不仅实现了节水节肥,而且产量及品质均有所提升,真正实现了提质增效的目标。
4 结 论
不同水肥耦合处理对苹果果实纵径、横径、产量和品质都有明显影响,其中W2F2的产量最高,同时其品质指标(单果质量、果形指数和可溶性固形物含量)也总体优于其他处理。此外,W2F2处理的水肥利用效率也表现较佳。基于NSGA-Ⅲ算法构建了苹果的最优灌水量和施肥量多目标回归模型,最终确定最优灌水量为821.11 mm,最优施肥量为783.22 kg/hm2,其对应的苹果产量、水肥利用效率及品质均优于表现较佳的W2F2处理。
新疆林果产业技术体系项目(XJLGCYJSTX04-2024)
新疆维吾尔自治区“十四五”重大专项(2023A02002-5)