高前慧 高欣梅 张静峰 姜 岩 李文忠
(1内蒙古自治区兴安盟种子管理站,乌兰浩特 137400;2内蒙古自治区兴安盟农业科学研究所,乌兰浩特 137400;3内蒙古自治区兴安盟土壤肥料站,乌兰浩特 137400)
兴农5号玉米新品种是兴安盟兴农种业有限责任公司,以改良系NS为母本、0556为父本杂交选育而成。2015年通过内蒙古自治区品种审定委员会审定(审定编号:蒙审玉2015023号),2016年在黑龙江省进行了新品种推广备案。适宜在≥10℃活动积温2100℃以上地区种植,属东北地区极早熟品种,目前在内蒙古及黑龙江推广面积超过1万hm2。该品种在实际推广中发现,有些地区的种植密度和施肥水平与审定公告推荐不符合时,仍然能获得较高产量。为了更精确地制定该品种在大田生产中的最佳农艺措施组合方案,2016年兴安盟兴农种业有限责任公司与兴安盟种子管理站、兴安盟农业科学研究所合作,进行了种植密度和施肥水平相关性试验,通过在不同密度下种肥和追肥的不同施用量对产量影响的回归方程进行研究,确定田间的最佳农艺措施组合,为该品种推广中的标准化栽培方案提供依据。
1.1 地块及材料试验在兴安盟科右前旗居力很镇万宝村的兴农种业试验园区内进行,该地块≥10℃活动积温2500℃,土质为栗钙土。试验品种为兴农5号,种肥选用复合肥(N、P2O5、K2O含量为:15%、15%、15%),追肥选用尿素(N含量46%)。
1.2 试验设计采用三因子二次正交旋转组合设计[1],各试验因子分别为密度(X1)、种肥施用量(X2)、追肥施用量(X3)。每个因子取5个水平(表1)。其中以品种审定时推荐值6.0万株/hm2作为密度0水平,步长为1.50万株/hm2;施肥量0水平是根据测土施肥配方确定。
经兴安盟土壤肥料站检测,该地块土壤中有机质含 量 3.8%,全 N 3.02g/kg、P 22.5mg/kg、K 168mg/kg,代入配方公式[2]:
N:Y=-7.0497×ln3.02+13.02=5.23kg/667m2,即78.42kg/hm2;
P2O5:Y=-3.1982×ln22.5+12.139=2.18kg/667m2,即32.72kg/hm2;
K2O:Y=-5.0566×ln168+27.955=2.05kg/667m2,即30.68kg/hm2。
为方便试验中操作,种肥的0水平设为某一养分纯量30.00kg/hm2,与测土配方中P2O5及K2O量相近,换算为复合肥量为200kg/hm2,并设步长为复合肥50kg/hm2;追肥的0水平施用量确定为纯N 46kg/hm2,与配方中纯N量减去种肥中纯N量后的数值相当,换算为尿素用量为100kg/hm2,并设步长为尿素25kg/hm2。
表1 试验因子水平表
1.3 试验方法根据试验设计中使用的正交表[1]结构,共23个处理组合,即23个小区(表2)。为了便于操作,减少相邻小区间的相互影响,小区设计为4行,行长为10.0m,行距0.6m,面积24.0m2,并设保护行。
表2 三因素回归正交设计结构及产量结果表
试验田前茬作物为玉米,有灌溉条件。2016年5月10日播种,为保证种植密度精确,播种时按设计标准等距离刨埯,抗旱坐水人工点播,每埯播种3粒。播种时施入种肥,每埯施入量精确称重。5月21日达到苗齐。6月3日定苗,定苗时严格控制每行应留株数。采用化学除草,灌溉等其他管理略高于当地生产。追肥在7月20日进行,追肥量按设计要求称重并人工施入。
在全生育期内对田间性状进行调查。收获期为10月5日,收取中间2行,面积12m2,以实际产量进行测产。
2.1 产量结果及回归方程各小区测产结果换算为含水量13%的公顷产量(表2)。通过对产量结果进行计算,最终得出回归系数(表3),并建立回归方程:
由于水平编码经过中心化处理,所以常数项要经过修正,表中X0所对应的b0为修正后值。
2.2 回归方程显著性检验回归方程是否可靠,需要对其进行F检验来评价,同时对回归系数的检验可明确各因子对产量影响的程度(表4)。通过对表4中的F值与相应的标准值比较可知:回归方程达到极显著水平,失拟检验不显著,说明回归方程可靠性强,可以依据其进行栽培方案的制作。
表3 回归系数bi计算结果表
表4 回归方程显著性检验
3.1 回归方程解析试验目的是为了指导农业生产。通过试验得到了回归方程,也就明确了试验中的各因子对产量影响的趋势和程度。于是便可通过对方程的解析[3],得到最优因子组合时的水平值,从而得到最佳农艺组合。
首先,依回归方程分别求出对应X1、X2、X3的偏导数。即:
在令其等于0后,求解可得X1=-0.2864、X2=-0.9600、X3=-1.4509把这3个值代入原回归方程,则可得到产量的一个极值Y(极值)=9486.28kg/hm2。这说明依据回归方程,当X1、X2、X3形成这一搭配时,即X1=-0.2864(相当于密度为5.57万株/hm2),X2=-0.9600(相 当 于 施 入 种 肥 151.99kg/hm2),X3=-1.4509(相当于施入追肥63.73kg/hm2),将获得产量为9486.28kg/hm2。显然,这一产量并非期望的,这是因为根据回归方程偏导求出的极值并不一定是极大值,也有可能为极小值,最简单的验证方法只需在极值周围取几个点代入方程中,当回归产量全部大于极值产量时,说明求得的很可能是最小值,本次试验便是这种情况。这时需要用另外一种方法来确定高产栽培方案。
3.2 高产栽培方案确定使用回归方程来计算极值点,会出现很多可能,除了以上出现的极小值外,还有可能出现没有极值的情况。同时,由于农业生产的特殊性,即使应用同样的栽培措施也不一定获得相同的产量。因此,在实际当中即使方程有最高值时,一般也不直接应用求极值点的对应组合作为高产栽培方案,而是使用“最佳模拟配合法”[1]来建立某一产量水平下的高产栽培组合方案。
对于本次试验,选定11000kg/hm2为目标产量,并选定在-1.5、-1、-0.5、0、0.5、1、1.5 7个量级组合,将它们分别代入回归方程中,可获得343个回归产量,其中11000kg/hm2以上回归产量139个。通过统计高于目标产量下不同因素水平出现的频数,用加权法计算平均数()、总体标准差(SX)。最后用-1.96SX~+1.96SX来确定不同农艺措施95%的置信区间,并换算为具体的农艺措施值(表5)。
表5 玉米兴农5号11000kg/hm2以上高产栽培农艺措施组合方案
从农艺措施的组合中可以看到,如果农艺措施中密度能够在6.30万~6.60万株/hm2,施用种肥在217.2~228.4kg/hm2,追肥在 110.6~116.1kg/hm2,则可获得11000kg/hm2以上的产量。因为试验中所用的种肥为复合肥,追肥为尿素,如换算成有效成分则种肥:N、P2O5、K2O 均在 32.58~34.26kg/hm2之间,追肥:N施用量在50.88~53.41kg/hm2之间。
在一般情况下,理论产量最高点在试验设计处理范围之内时,可以求出最高理论产量,但在实际中,除了无法求得极值外,还会出现理论产量最高点向外飘移的情况,这是因为试验的设计模型强烈依靠调查数值,当一个因素的观测值出现在处理水平最高或最低点时,都会出现这种情况。这也是在实际生产中,通常使用“最佳模拟配合法”的一个原因。通常情况下,制定高产栽培方案时,下一年可以在继续重复试验验证的同时进行生产验证,起到试验示范的双重目的。
[1] 陶向新,王伯伦.作物规范化栽培设计原理与应用[M].沈阳:辽宁大学出版社,1991
[2] 王学阁.兴安盟主栽作物主推技术读本[M].赤峰:内蒙古科学技术出版社,2013:229-243
[3] 张健,冯云超,余志江,等.不同密度及施肥措施对玉米新品种三峡玉9号产量的影响[J].中国种业,2016(1):51-55