谷子[
Setaria italica(L.)Beauv]是起源于中国的传统粮食作物,具有悠久的栽培历史,其耐瘠薄、抗旱性强的特点使得谷子成为旱作地区产业结构调整的主栽作物
[1-2]。宁夏位于西北内陆高原,属典型的大陆性气候,干旱少雨,谷子种植主要分布在宁夏中部干旱带和南部山区,属于西北春谷区,种植面积常年稳定在30万亩左右
[3]。但随着对谷子经济效益的最大化追求,同一地块连年种植谷子而引起的连作障碍问题已经非常突出
[4],也已成为制约当地谷子产业可持续发展的主要因素。
诸多研究表明,轮作等农业措施是缓解连作障碍的最健康和有效的途径之一,合理轮作能够促进土壤养分的均衡利用,改善土壤生态环境,并提高作物产量
[5-9],而土壤生态系统的变化也是导致谷子连作障碍形成的重要原因
[10]。牛倩云等
[11]发现大豆-马铃薯-谷子轮作能改善谷田生态环境,提高土壤有机质和速效养分含量,使谷子产量较连作显著增加37.11%,进而有效缓解连作障碍的发生。王仪明等
[12]发现,苜蓿-青贮玉米轮作中土壤有机质、速效氮和速效钾含量分别增加了60.79%、110.84%和160.02%。李洋等
[13]发现,小麦—花生和小麦—大豆轮作能够明显增加土壤中的速效氮含量,并且保持小麦产量稳定提升。张成君等
[14]研究表明,在西北半干旱区,不同轮作模式对土壤养分具有显著影响,其中玉米-豌豆和马铃薯-燕麦轮作模式提高了土壤有机质、全氮、速效氮、速效磷和速效钾含量。李丹等
[15]研究谷子不同轮作模式对土壤理化性质及细菌群落时发现,不同轮作模式对土壤质量有一定影响,其中土壤全氮、有机质含量提高程度因模式而异,玉米是避免谷子连作障碍较为适宜的作物,但在西北旱作区其它生态条件下的谷田是否适用,尚不清楚。
目前,旱作区轮作制度的研究大多数集中在单一轮作种植模式对作物生长及土壤环境影响,而西北旱作区多种主要粮食作物与谷子进行年际间轮作对土壤肥力、作物生长特性及产量等方面综合影响的相关报道仍然较少。因此,本试验以宁夏旱作区优势抗旱作物谷子为研究对象,采用谷子与不同主栽作物年际间轮作的模式,探究谷子不同轮作模式下土壤养分2年间的变化,以及轮作后对谷子生长特性、产量、经济效益的影响,旨在缓解连作障碍引起的土壤肥力失衡,并为优化谷子轮作模式提供理论依据,从而指导西北旱作区谷子产业可持续绿色发展。
1 材料与方法
1.1 试验材料
供试谷子品种为‘张杂谷13号’,由河北巡天种业公司提供;玉米品种为‘赛博159’,由河北垦丰农业科技有限公司提供;马铃薯品种为‘青薯9号’,由海原县恒兴马铃薯专业合作社提供;高粱品种为‘晋牧4号,’由山西农业大学高粱研究所提供;糜子品种为‘固糜21号’,由宁夏农林科学院固原分院提供。覆盖材料为渗水地膜,宽1650 mm、厚0.01 mm,由山西微通渗水膜生物科技有限公司提供。
1.2 试验区概况
试验在宁夏回族自治区中卫市海原县贾塘乡王塘村(36°52′ N、105°90′ E)进行,该地区属中温带半干旱大陆性季风气候,海拔1750 m,年平均降雨量为305.3 mm,平均气温15.6 ℃,≥10 ℃积温为2550~2700 ℃。试验区土质为黄绵土,前茬作物为谷子,土壤基本理化性质为:pH=8.08,有机质含量10.44 g/kg,全氮含量0.737 g/kg,速效氮含量40.89 mg/kg,全磷含量0.59 g/kg,速效磷含量13.88 mg/kg,全钾含量18.25 g/kg,速效钾含量162 mg/kg。试验期间即2021和2022年作物生育期降雨量分别为176.7、353.8 mm,作物生育期气象条件见
图1。
1.3 试验设计
以谷子连作3年为对照(M-M),设置玉米—谷子(C-M)、马铃薯—谷子(P-M)、高粱—谷子(S-M)、糜子—谷子(B-M)4种轮作模式,共计5个处理。2021年(第1年)种植作物分别为谷子、玉米、马铃薯、高粱、糜子,2022年(第2年)种植作物全部为谷子。每个处理3次重复,随机区组设计,每个小区面积6 m×5 m。
各类作物种植方式均为地膜覆盖波浪式穴播。谷子于2021年4月18日和2022年4月25日播种,播种量7.5 kg/hm2,种植密度30万株/hm2;玉米于2021年4月28日播种,播种量30 kg/hm2,种植密度8.5万株/hm2;马铃薯于2021年5月3日播种,播种量3000 kg/hm2,种植密度5.25万株/hm2;高粱于2021年4月25日播种,播种量22.5 kg/hm2,种植密度8.25万株/hm2;糜子于2021年5月15日播种,播种量22.5 kg/hm2,种植密度90万株/hm2。结合春季整地统一施完全腐熟的优质农家肥 30 000 kg/hm2、尿素180~225 kg/hm2、过磷酸钙120~150 kg/hm2、硫酸钾75~90 kg/hm2,各类作物其它田间管理采用常规管理。
1.4 测定指标与方法
1.4.1 谷子农艺性状
在种植第2年谷子成熟期,每个小区随机选取生长均匀一致的谷子20株,测定株高、茎粗(地上第2节的粗度)、穗长、倒二叶(旗叶以下数第二片叶子)叶长和叶宽。并通过计算叶面积,进一步求得叶面积指数(LAI)。
1.4.2 谷子产量及产量性状
在种植第2年谷子收获期,按照小区全部收获并脱粒、晒干、称重,每个小区随机选取生长均匀一致的谷子20株进行室内考种,按照常规方法测定千粒重、主穗粒重、主穗重。
1.4.3 土壤养分指标
在各类作物种植第1年与第2年的成熟期,每个小区采用五点采样法收集土壤样品,采样深度0~20 cm,去除大块石砾、植株残根等杂质后,将5个点的土样混合均匀,取1 kg土,风干后过1 mm筛,测定土壤有机质、全氮、全磷、全钾、速效氮、速效磷、速效钾含量。
土壤有机质含量采用重铬酸钾氧化外加热法测定;土壤全氮含量采用凯氏蒸馏法测定,速效氮含量采用碱解扩散法测定;土壤全磷含量用高氯酸消解-钼锑抗比色法测定,速效磷含量采用碳酸氢钠浸提分光光度计法测定;土壤全钾含量用氢氧化钠熔融-火焰光度计法测定,速效钾含量采用乙酸铵浸提火焰光度计法测定
[16]。
1.4.4 综合评价方法
参考王仪明等
[12]的方法将原始变量设为
p个,分别为X
1,X
2,…,X
p,将主成分看作因变量,原始变量X看作自变量,对其进行线性组合得到相应的主成分数学模型:
F1=a11X1+a12X2+…+a1p X p
F2=a21X1+a22X2+…+ a2p X p
…
F p =a p1X1+a p2X2+…+a p X p
称F1,F2,…,FP分别为第1、第2、…、第p个主成分;系数向量(a1j,a2j,…,apj ),(j=1,2,…,p)是协方差阵的特征值(λ1≥λ2≥…≥λ p )所对应的特征向量,且F1,F2,…,F p 的方差是上述的特征值,主成分的名次按特征值大小的顺序进行排列。
第i主成分的贡献率为λ i /∑λ i,表明第i个主成分综合X1,X2,…,X p 的信息的能力。
1.4.5 经济效益计算方法
经济效益分别按照作物成本、产值、总成本、总产值、纯收入、产投比计算
[17]。
作物成本(元/hm2)=种子成本(元/hm2)+肥料成本(元/hm2)+地膜成本(元/hm2)+机械成本(元/hm2)+田间管理成本(元/hm2),各类成本均为试验当年成本
作物产值(元/hm2)=作物产量(kg/hm2)×作物价格(元/kg),作物产值为试验当年产值
总成本(元/hm2)=种植第1年的作物成本(元/hm2)+种植第2年的作物成本(元/hm2)
总产值(元/hm2)=种植第1年的作物产值(元/hm2)+种植第2年的作物产值(元/hm2)
纯收入(元/hm2)=总产值(元/hm2)-总成本(元/hm2)
产投比=总产值(元/hm2)/总成本(元/hm2)
1.5 数据处理与分析
采用Excel 2007软件进行数据汇总,利用DPS 7.05软件进行方差分析和因子分析,采用新复极差法进行多重比较。利用OriginPro 2021软件进行作图。
2 结果与分析
2.1 不同轮作模式下谷子农艺性状
由
表1可见,不同处理下谷子成熟期的株高在88.03~104.50 cm,C-M和B-M处理均显著高于其它处理(
P<0.05),P-M和S-M处理与M-M处理无显著差异(
P>0.05)。轮作模式下谷子茎粗和穗长分别比M-M高6.33%~14.02%、4.69%~18.41%,但各处理间无显著差异(
P>0.05),其中B-M茎粗最高,为9.19 mm,P-M穗长最高,为24.51 cm。轮作模式下谷子叶长和叶宽均高于M-M,其中C-M叶长较M-M显著增加23.76%,C-M、S-M和B-M叶宽较M-M分别显著增加19.27%、26.45%、26.34%(
P<0.05)。叶面积指数各处理高低依次为C-M>B-M>S-M>P-M>M-M,各轮作模式较M-M显著提高24.49%~47.96%(
P<0.05)。
2.2 不同轮作模式下谷子产量性状及产量
不同轮作模式对谷子产量及产量性状的影响见
表2。各处理下谷子的产量大小依次为S-M>P-M>B-M>C-M>M-M,其中P-M和S-M处理均显著高于其它处理(
P<0.05),分别较M-M处理增产19.00%、19.27%,C-M和B-M与M-M差异不显著(
P>0.05)。从谷子产量相关性状可以看出,谷子主穗粒重和籽粒千粒重各处理间差异不显著(
P>0.05),而主穗重与产量的变化趋势表现一致,说明轮作模式下谷子主穗重可能对其产量有着显著的影响。
2.3 不同轮作模式下土壤养分含量
2.3.1 不同轮作模式下土壤有机质含量
由
图2可见,种植第1年各处理土壤有机质含量依次为S-M>C-M>P-M>B-M>M-M,其中S-M和C-M处理分别较M-M处理显著增加54.55%、28.68%(
P<0.05),P-M和B-M处理与M-M处理差异不显著(
P>0.05);第2年C-M、P-M和S-M处理土壤有机质含量均显著高于其它处理,与M-M处理相比分别显著增加34.45%、39.49%、30.25%(
P<0.05),B-M与M-M差异不显著(
P>0.05)。各轮作模式下土壤有机质含量2年平均值依次为S-M>C-M>P-M>B-M,均高于M-M处理,增加0.04%~43.02%。
2.3.2 不同轮作模式下土壤全量养分含量
由
图3可见,各处理土壤全氮含量2年平均值依次为S-M>P-M>C-M>B-M>M-M,其中C-M、P-M和S-M轮作模式较M-M分别显著提高46.89%、47.27%、59.53%(
P<0.05);种植第1年各轮作模式土壤全氮含量均显著高于M-M处理,增加42.28%~124.26%(
P<0.05),种植第2年C-M、P-M和S-M处理分别较M-M处理显著提高19.75%、30.86%、16.05%(
P<0.05)。各处理土壤全磷含量2年平均值依次为P-M>S-M>M-M>B-M>C-M,其中第1年P-M、S-M和B-M处理分别较M-M处理显著提高21.62%、10.81%、8.11%(
P<0.05),第2年P-M处理较M-M处理显著提高19.74%(
P<0.05)。各轮作模式土壤全钾含量种植第1、2年均与M-M处理差异不显著,其中S-M处理全钾含量最高,较M-M处理2年平均提高4.12%。
2.3.3 不同轮作模式下土壤速效养分含量
由
图4可见,各处理土壤速效氮含量2年平均值依次为S-M>B-M>C-M>M-M>P-M,其中种植第1年S-M和B-M处理较M-M处理显著提高56.41%、43.34%(
P<0.05),种植第2年C-M和P-M处理较M-M处理显著提高14.71%、17.74%(
P<0.05)。轮作模式土壤速效磷含量平均值较M-M处理提高4.03%~88.66%,依次为S-M>P-M>C-M>B-M>M-M,其中种植第1年P-M和S-M处理较M-M处理显著提高142.55%、58.19%(
P<0.05),第2年S-M处理较M-M处理显著提高103.66%(
P<0.05)。土壤速效钾含量2年平均值依次为S-M>P-M>B-M>C-M>M-M,轮作模式较M-M处理提高1.19%~88.07%,其中第1年S-M处理较M-M处理显著提高111.38%(
P<0.05),第2年P-M和S-M处理分别较M-M处理显著提高47.98%、54.91%(
P<0.05)。
2.4 不同轮作模式对土壤肥力影响的综合评价
采用因子分析法,将谷子不同轮作模式下土壤养分中有机质、全氮、全磷、全钾、速效氮、速效磷和速效钾含量7个指标依次作为自变量X
1~X
7,并进行因子分析。由
表3可见,7个因子中,第1和第2公因子累计贡献率达到83.87%,说明7个因子中绝大部分相关信息可由前2个公因子概括。
根据第1和第2公因子的特征值和贡献率计算综合分值并排序,由此得出谷子不同轮作模式和连作下土壤养分的综合评分。结果如
表4所示,各处理排名依次为P-M、S-M、C-M、M-M、B-M。除B-M模式外,其它轮作模式得分均高于连作,得分越高轮作模式优势更加突出,说明C-M、P-M和S-M模式改善土壤肥力的效果更佳。
2.5 不同轮作模式下经济效益分析
由
表5可以看出,种植第1年和第2年各处理投入总成本为P-M>C-M>M-M>S-M>B-M;总产值和纯收入为P-M>S-M>M-M>C-M>B-M,其中P-M和S-M处理分别较M-M处理总产值增加23.36%、10.91%,纯收入增加24.29%、15.78%;产投比为S-M>P-M>M-M>C-M>B-M,P-M和S-M处理分别较M-M处理提高1.75%、11.95%。综上,在各轮作模式中,P-M和S-M这2种轮作模式在纯收入和产投比上较连作具有明显优势。
3 讨论
土壤中有机质及氮磷钾是土壤肥力的物质基础,建立合理的轮作制度能够打破长期连作过程中土壤养分的失衡,促进土壤养分含量的积累和再循环利用
[18-19]。孙倩等
[20]研究发现,谷子与大豆轮作后,土壤中有机质、全氮、全磷和全钾含量均不同程度提高;左琴等
[21]研究水稻与不同作物轮作表明,与水稻连作相比,各轮作模式下土壤的有机质、全磷、全钾、速效氮和速效钾含量均不同程度地增加;宋佳承等
[22]在不同种植模式对土壤质量及马铃薯生长的影响研究中表明,马铃薯与藜麦、玉米轮作均增加了土壤中有机质、速效氮和有效磷含量。本试验通过研究不同轮作模式下土壤养分变化,发现轮作模式均不同程度提高了土壤的有机质含量,可能受轮作降低土壤中有机碳矿化速率的影响
[23],并且土壤中全氮、速效磷和速效钾含量也增加,这与妙佳源等
[24]研究中谷子轮作下土壤氮、磷含量较连作增加的结果基本一致。谷子连作条件下,土壤有机质、速效氮、全钾和速效钾含量均呈现下降趋势,土壤全氮和全磷含量增加,这与郝晓芬等
[25]在谷子连作2年、3年和5年的研究中,发现土壤养分呈先下降后上升趋势的结论相反,主要原因可能是试验区基本以施氮肥和磷肥为主,轮作第1年降雨量仅为第2年的二分之一,水分缺失抑制了土壤中氮磷全量养分的转换。本研究中,与玉米、马铃薯和高粱轮作后,谷田土壤的有机质和全氮含量显著增加,在高粱与谷子轮作模式下速效磷和速效钾的含量也显著增加,说明合理轮作后谷田土壤供应养分的能力明显有所提升。对土壤肥力进行综合评价,P-M、S-M、C-M轮作模式下土壤养分的综合评分均为正,而B-M轮作模式和谷子连作均为负。主要的原因可能是糜子和谷子同属一年生草本谷类作物,二者的土壤耕作层结构以及根系生长环境相似,导致糜子在改善土壤肥力环境方面有一定的局限性;而玉米和高粱具有发达的根系,马铃薯属块茎作物,这3种作物对土壤耕作层结构的改变影响着土壤养分的固定和转化,这与周旭东等
[26]关于土壤结构的改变促进了作物对养分的摄入量的结论相似。
轮作是影响土壤质量的重要措施,且直接影响作物的生长状况和产量水平
[27]。曹鑫鑫等
[28]研究发现,轮作模式和有机肥配施有效提高了谷子各生育时期内的株高、抽穗/灌浆期植株的茎粗和营养生长阶段功能叶叶面积,有利于谷子植株的形态建成。李小霞等
[29]研究发现,玉米-大豆-谷子、玉米-花生-谷子和玉米-马铃薯-谷子轮作模式下的谷子产量较连作显著增加3.07%~9.98%。本研究中,不同轮作模式不仅增加了谷子茎粗,还提升了叶面积指数,说明耕作制度的改变导致谷子农艺性状变化,由此为谷子生长后期主穗粒重与主穗重的增加提供了物质基础,进而最终实现P-M和S-M轮作模式下的谷子产量较连作显著提高。
不同轮作模式可以提升土壤质量和谷子产量,但农田整体的生产力也因为轮作作物的差异而不同。与谷子连作相比,轮作马铃薯、高粱的总体纯收入和产投比具有明显的优势,主要原因是马铃薯和高粱当年的产值高于其它作物,但生产成本差异较小,且经过轮作后,谷子生产力也明显提升。因此,在西北旱作区旱作农业的发展过程中,谷子与马铃薯和高粱短期轮作是实现谷子增产增效的重要轮作模式。
4 结论
本研究结果表明,不同轮作模式对谷子生长特性及土壤肥力有明显影响,通过谷田土壤肥力综合评价和经济效益比较可知,马铃薯—谷子和高粱—谷子2种轮作模式下,土壤肥力综合得分、经济效益和产投比均较其它轮作与连作模式明显提升,有助于谷子产量和谷田生产力的提升,可为谷子轮作制度的建立提供理论依据。由于轮作制度对作物连作障碍的缓解是多方面综合因素的影响,还需对土壤中微生物的响应等微生态环境变化进行长期的深入研究。
宁夏回族自治区优秀青年基金项目(2022AAC05048)
国家现代农业产业技术体系(CARS-06-14.5-A18)
中央引导地方科技发展专项(2024FRD05091)
宁夏农业科技自主创新专项(NNKZZCGZH-2023-02)
宁夏农业科技自主创新专项(NKYG-25-17)
小杂粮种质资源鉴评及优质抗逆新品种选育(2024FRD05091)