江西是中国花生主要产区,花生常年种植面积约18.23万hm
2,位居全国第十,发展江西花生生产对保障国家油料安全至关重要
[1-2]。叶片干重和氮含量是评价花生长势质量与氮营养丰缺的重要生长参数,直接影响花生冠层光能截获利用及丰产稳产
[3-4]。因此,实时准确获取花生叶片干重和氮含量信息对花生精确诊断施肥和定向调控管理具有重要价值
[5]。
已有研究表明光谱技术可实时、准确和无损预测作物生长参数
[6-8],克服传统人工取样测定法时效性差、费时耗工及大面积快速预测难的问题
[9-10]。近年来,研究主要利用无人机搭载高光谱和多光谱相机采集作物冠层光谱影像,进而构建基于特征光谱植被指数的株高
[11]、叶面积指数
[11-12]、氮营养
[13-14]、生物量
[15-16]、叶绿素含量
[17]及产量
[18]等生长参数预测模型。这些模型具有监测范围大、分辨率高、数据采集高效、对作业环境要求低等优点。但无人机搭载高光谱和多光谱相机大多为国外进口,对用户及影像数据处理要求高、专业性强、购机成本高,在作物实际生产应用中具有一定的局限性。为此,国内已有科研机构设计开发了使用简单、精度较高及成本较低的手持式光谱仪
[19-21],建立了准确实用的作物生长参数预测诊断模型,在作物精确栽培中表现出明显的省工节肥作用
[22-24]。但是由于研究对象、作物株型、光谱仪型号及生产管理方式等不同,导致模型结构及系数存在差异、适用性较差,要提高模型预测性和适用性需开展属地化建模及检验。
江西花生主要种植在丘陵红壤旱地、土壤酸性较强,手持式光谱仪比无人机载光谱相机更适用于获取种植在小地块的花生生长参数信息。此外,江西红壤旱地花生生长参数和冠层光谱植被指数受株型影响而存在差异,花生生长季以多云寡照天气为主,被动光源的手持式和无人机载式光谱仪在南方多云寡照天气下应用具有明显的局限性,且目前基于国产手持式主动光源光谱仪的红壤旱地不同株型花生叶片干重和氮含量预测研究鲜有报道。为此,本研究拟以红壤旱地主推的松散型品种‘赣花5号’和紧凑型品种‘粤油93’为研究对象,开展不同株型花生品种与施氮量的田间试验,在花生生长过程中主要生育期采用国产手持式主动光源光谱仪测定光谱植被指数,取样测定叶片干重和氮含量,构建基于光谱植被指数的花生叶片干重和氮含量预测模型,以期为红壤旱地花生生长参数实时准确获取及精确栽培管理提供技术支撑。
1 材料与方法
1.1 试验设计
本试验于2023年和2024年4—8月在江西省南昌市进贤县开展。试验地为典型丘陵红壤旱地,耕作层含全氮1.09 g/kg、有机质15.12 g/kg、速效钾106.32 mg/kg、速效磷20.15 mg/kg。选择2个当地主推花生品种,依次为‘赣花5号’(松散型)和‘粤油93’(紧凑型),4月3日播种,8月9日收获。设置4种施氮量,依次为纯氮0、67.5、135和202.5 kg/hm2(依次用N0、N1、N2和N3表示,N2是当地高产栽培施氮量),分2次施用,70%作基肥,30%在主茎3~5片真叶时追施。另配P2O5 81 kg/hm2 和K2O 135 kg/hm2,作基肥一次施用。小区面积40 m2,条播,行距40 cm,穴距20 cm,每穴播2粒,南北行向。采用随机区组设计,3次重复。氮、磷、钾肥分别为尿素、钙镁磷肥和氯化钾。
1.2 测定项目及方法
1.2.1 光谱植被指数测定
在花生苗期、花针期、结荚期和饱果期,采用南京农业大学研制的手持式光谱仪
[20]测定花生冠层光谱植被指数。该光谱仪(CGMD402型)为主动光源,含730 和 810 nm共2个特征波段,视场角25°,开机自校,可直接测得比值植被指数、差值植被指数及归一化植被指数,计算公式详见
表1。数据采集时间为10:00—14:00,仪器传感器垂直于花生冠层,与冠层相距1 m,每个试验小区随机选取3个观测点,每个观测点测定5次,计算平均值得到该小区测定值。
1.2.2 花生叶片干重及氮含量测定
在花生苗期、花针期、结荚期和饱果期,与光谱植被指数的测定同步,在各试验小区选取代表性花生植株5株,将植株样品取回实验室按叶、茎(包含根)和果等器官分开后装入信封,在105 ℃条件下杀青30 min,在80 ℃条件下烘48 h至恒量后分别称重,根据花生种植密度计算得到单位土地面积上的叶片干重(g/m
2),采用凯氏定氮法
[25]测定花生叶片氮含量(%)。
1.3 预测模型的构建及检验
使用Excel 2021进行数据整理作图及拟合分析;使用SAS 8.0中的LSD 法进行单因素方差分析检验不同处理间差异显著性,用Shapiro-Wilk法对数据进行正态性分布检验。以试验测得的3个光谱植被指数(RVI、DVI和NDVI)作为自变量,叶片干重和氮含量作为因变量,依次进行线性方程、多项式方程、指数方程、对数方程及幂函数方程拟合比较分析(建模训练集样本量为24,测试集样本量为24),以拟合方程决定系数(Determination Coefficient, R2)值大小确定相关性高低,选择R2值最大的方程为最佳预测模型。通过对比分析均方根误差(Root Mean Square Error, RMSE)、相对均方根误差(relative Root Mean Square Error, rRMSE)和相关系数(Correlation Coefficient, r)来检验评价模型的预测效果与准确性。
2 结果与分析
2.1 红壤旱地花生叶片干重变化规律
不同生育期和施氮量下的红壤旱地花生叶片干重变化规律见
图1。可见施氮量显著影响红壤旱地花生叶片干重。其中:2个花生品种不同生育期的叶片干重随施氮量增大而增高,相同花生品种在不同施氮量间具有极显著差异(
P<0.01)(
图1(a)和(b));不施氮肥处理N0叶片干重较低,直接影响叶片光合速率及光合产物的形成;N3处理的叶片干重显著高于N0、N1和N2处理,高氮条件下易导致花生贪青晚熟(
图1(a)和(b))。在施氮量相同的条件下,花生叶片干重从苗期至饱果期呈先增后减的变化。如紧凑型品种‘粤油93’N3处理苗期、花针期、结荚期和饱果期的叶片干重依次为82.63、228.71、263.03和161.98 g/m
2(
图1(b))。松散型品种‘赣花5号’结荚期N0、N1、N2和N3的叶片干重依次为128.28、190.69、237.63和259.19 g/m
2(
图1(a));紧凑型品种‘粤油93’结荚期N0、N1、N2和N3的叶片干重依次为130.22、191.31、244.13和263.03 g/m
2(
图1(b))。株型紧凑的花生品种具有叶片厚和光能利用率高的特性,2个供试花生品种相比,‘粤油93’株型紧凑,其叶片干重高于‘赣花5号’。
2.2 红壤旱地花生叶片氮含量变化规律
不同生育期和施氮量下的红壤旱地花生叶片氮含量变化规律见
图2。可见施氮量显著影响红壤旱地花生叶片氮含量。其中:2个花生品种不同生育期的叶片氮含量随施氮量增大而增高,相同花生品种在不同施氮量间具有极显著差异(
P<0.01)(
图2(a)和(b))。不施氮肥处理N0叶片氮含量较低,直接影响叶片光合作用及物质积累;N3处理的叶片氮含量显著高于N0、N1和N2处理(
图2(a)和(b))。在施氮量相同的条件下,花生叶片氮含量从苗期至饱果期呈先增后减的变化。如松散型品种‘赣花5号’N2处理苗期、花针期、结荚期和饱果期的叶片氮含量依次为1.99%、2.83%、2.98%和2.38%(
图2(a))。松散型品种‘赣花5号’饱果期N0、N1、N2和N3的叶片氮含量依次为1.68%、2.16%、2.38%和2.49%(
图2(a));紧凑型品种‘粤油93’饱果期N0、N1、N2和N3的叶片氮含量依次为1.64%、2.12%、2.34%和2.61%(
图2(b))。总之,2个供试花生品种相比,‘赣花5号’为耐瘠薄品种,‘粤油93’为耐肥品种,‘赣花5号’N0、N1和N2处理的叶片氮含量高于‘粤油93’,其N3处理的叶片氮含量低于‘粤油93’。
2.3 红壤旱地花生叶片干重及氮含量光谱预测模型构建
利用2023年测定的试验数据构建红壤旱地花生叶片干重和氮含量光谱预测模型。基于光谱植被指数的不同生育期红壤旱地花生叶片干重和氮含量光谱预测模型见
表2。可见:不同生育期RVI与花生叶片干重相关性最高的预测模型为二次多项式模型,模型构建
R2为0.76~0.82;不同生育期DVI、NDVI与花生叶片干重相关性最高的预测模型均为幂函数模型,模型构建
R2依次为0.81~0.90和0.82~0.91;不同生育期RVI、DVI、NDVI与花生叶片氮含量相关性最高的预测模型均为二次多项式模型,模型构建
R2依次为0.70~0.87、0.74~0.86和0.84~0.90。因此,比较RVI、DVI和NDVI 3个光谱植被指数,NDVI与花生叶片干重和氮含量的建模
R2值更大,表明NDVI与花生叶片干重和氮含量的相关性更高。
2.4 红壤旱地花生叶片干重光谱预测模型检验
利用2024年试验数据对本研究构建的红壤旱地花生叶片干重光谱预测模型进行检验评价。基于NDVI的红壤旱地花生叶片干重观测值与预测值的比较见
图3。可见基于NDVI的幂函数模型预测苗期(
图3(a))、花针期(
图3(b))、结荚期(
图3(c))和饱果期(
图3(d))花生叶片干重的均方根误差(RMSE)、相对均方根误差(rRMSE)和相关系数(
r)依次为4.57~19.07 g/m
2、7.77%~11.23%和0.911 8~0.973 0。表明不同生育期的花生叶片干重观测值与预测值间具有较高的一致性,模型具有较好的预测效果与准确性。
2.5 红壤旱地花生叶片氮含量光谱预测模型检验
利用2024年试验数据对本研究构建的红壤旱地花生叶片氮含量光谱预测模型进行检验评价。基于NDVI的红壤旱地花生叶片氮含量观测值与预测值的比较见
图4。可见基于NDVI的二次多项式模型预测苗期(
图4(a))、花针期(
图4(b))、结荚期(
图4(c))和饱果期(
图4(d))花生叶片氮含量的RMSE、rRMSE和
r依次为0.09%~0.16%、3.58%~6.42%和0.917 6~0.972 7。表明不同生育期的花生叶片氮含量观测值与预测值间具有较高的一致性,模型具有较好的预测效果与准确性。
3 讨 论
丰产、优质和高效是目前江西红壤旱地花生生产的主要目标,对花生生长参数进行实时准确预测诊断是实现这个目标的基础
[5]。本研究通过实施不同株型花生品种与施氮量的小区观测试验,测得不同生育期的红壤旱地花生冠层光谱植被指数、叶片干重和氮含量数据。结果发现,红壤旱地花生叶片干重和氮含量随施氮量增大而增高,从苗期至饱果期呈先增后减的变化,这与刘佳等
[26]研究结论一致。表明花生生长前期主要进行营养生长,随着施氮量增加,养分被植株吸收转运,叶片干重和氮含量逐渐增高,结荚期后花生下部叶片出现衰亡脱落以及养分逐渐向荚果转运,叶片干重和氮含量呈递减变化。株型是表征花生冠层结构和群体光合能力的重要性状,株型紧凑的花生品种具有叶片厚及光能利用率高的特性
[27]。本研究发现,株型紧凑的花生品种叶片干重较大;耐瘠薄的花生品种在低氮条件下叶片氮含量较高,耐肥的花生品种在高氮条件下叶片氮含量较高。证明施氮量显著影响红壤旱地花生叶片干重和氮含量,进而调控群体质量与产量形成,氮肥科学定量投入是红壤旱地花生丰产稳产栽培的关键。
近年来,监测面积大、机动灵活、无损快速、高分辨率的无人机遥感技术在农作物生长参数预测反演中被广泛应用
[28-30]。利用无人机搭载高光谱和多光谱相机可实时、准确、大范围获取农作物生长参数信息,但是其使用门槛和购机成本较高,在红壤旱地花生主产区推广应用实用性不强。为此,本研究通过实施不同年份、不同花生品种及不同施氮量的小区试验,采用手持式光谱仪测定花生冠层光谱植被指数,利用2023年试验数据构建了基于RVI、DVI和NDVI的红壤旱地花生叶片干重和氮含量光谱预测模型,利用2024年试验数据对构建的叶片干重和氮含量光谱预测模型进行了检验评价。结果发现,RVI与花生叶片干重相关性最高的预测模型为二次多项式模型,DVI、NDVI与花生叶片干重相关性最高的预测模型均为幂函数模型,RVI、DVI、NDVI与花生叶片氮含量相关性最高的预测模型均为二次多项式模型,是对周晓楠等
[31]仅开展线性拟合建模研究的拓展。比较分析RVI、DVI和NDVI 3个光谱植被指数发现,NDVI与花生叶片干重和氮含量的建模
R2值更大、相关性更高。
本研究利用NDVI构建的花生叶片干重预测模型
R2为0.82~0.91,模型检验的RMSE、rRMSE和
r依次为4.57~19.07 g/m
2、7.77%~11.23%和0.911 8~0.973 0,比陈龙跃等
[32]对花生生物量的预测更准确。利用NDVI构建的花生叶片氮含量预测模型
R2为0.84~0.90,模型检验的RMSE、rRMSE和
r依次为0.09%~0.16%、3.58%~6.42%和0.917 6~0.972 7,与Zhang等
[33]对大豆冠层氮含量的预测精度相当。证明可以采用手持式光谱仪测得光谱植被指数进而实时准确反演花生叶片干重和氮含量信息,适用于江西红壤旱地花生生长参数的实时无损获取,可克服人工测定法时效性差及用工成本高的不足。此外,本研究采用国产手持式主动光源光谱仪,可全天候实时无损预测花生生长参数信息,可克服以往被动光源光谱仪
[22,31]在南方持续多云寡照天气条件下难以全天候应用的难题,减少对国外进口光谱传感器的依赖和降低成本。
当然,为提高本研究构建的花生叶片干重和氮含量光谱预测模型的预测精度及普适性,今后还需利用不同地点、年份及栽培管理措施的试验数据对模型进行验证评价,进而为花生精确栽培与智慧管理提供技术支撑。
4 结 论
为实时准确获取红壤旱地花生生长参数、实现花生精确诊断施肥和定向调控管理,以红壤旱地主推的不同株型花生品种为研究对象,构建基于光谱植被指数的花生叶片干重和氮含量预测模型。主要结论如下:红壤旱地花生叶片干重和氮含量随施氮量增大而增高,从苗期至饱果期呈先增后减的变化;比较分析手持式光谱仪测定的3个光谱植被指数发现,NDVI与叶片干重和氮含量的相关性最高;利用NDVI构建的幂函数模型和二次多项式模型可分别精确预测红壤旱地花生叶片干重和氮含量,克服人工测定法时效性差及用工成本高的不足。
江西省农机装备研发制造推广应用一体化试点项目(YCTY202405)
国家自然科学基金项目(31960364)