摘要：针对半干旱区退化紫花苜蓿人工草地饲草产量低、品质差问题，试验采用单因素随机区组设计，研究了不同深松浅旋程度：深松 30 cm+浅旋 5 cm(S1Q1 )，深松 30 cm+浅旋 10 cm(S1Q2 )，深松 40 cm+浅旋 5 cm(S2Q1 )，深松 40 cm+浅旋 10 cm(S2Q2 )，深松 50 cm+浅旋 5 cm(S3Q1 )，深松 50 cm+浅旋 10 cm(S3Q2 )以及不做处理(CK)对半干旱区退化紫花苜蓿草地土壤理化性质、生产性能及饲草品质的影响，并利用主成分分析(PCA)方法评价其改良效果。 3 年试验结果表明，深松浅旋能够不同程度降低紫花苜蓿草地土壤容重，增加土壤孔隙度，显著提高紫花苜蓿株高、分枝数和叶茎比。其中，S2Q1处理可降低 0~40 cm 土层土壤容重，S2Q2处理可提高紫花苜蓿分枝数、干草产量和粗蛋白含量，降低中性洗涤纤维、酸性洗涤纤维含量。经 PCA 综合分析，紫花苜蓿干草产量、叶茎比和粗蛋白贡献率较大，牧草相对饲喂价值和中性洗涤纤维贡献率较小;深松浅旋 S2Q2 处理能够获得较高的紫花苜蓿干草产量(6505. 44 kg·hm-2)和粗蛋白含量(20. 74%);综合性状排名由高到低依次为 S2Q2，S1Q2，S2Q1，S3Q2，S1Q1，S3Q1，CK。由此说明，深松 40 cm+浅旋 10 cm 对半干旱区退化紫花苜蓿草地改良效果最优。
关键词：深松浅旋;退化苜蓿草地;产量;品质
紫花苜蓿(Medicago sativa)作为一种优良的多年生豆科牧草，因其产草量高、营养丰富、适口性好，被誉为 “牧草之王”[1-2]。宁夏是我国重要的畜牧业基地之一，紫花苜蓿是其重要的饲草来源，目前宁夏紫花苜蓿人工草地约 40. 0 万 hm2 ，其中半干旱区达 25. 3 万 hm2 ，但受该地区自然条件、管理水平等因素的影响，6~7 年后紫花苜蓿草地土壤紧实，含水量低，苜蓿生长受阻，植株低矮，分枝数少，植被稀疏，显现退化迹象。经调查，宁夏半干旱区 7 龄以上紫花苜蓿草地超过 70%，产量不足 3000 kg·hm-2。对于退化紫花苜蓿草地，一方面可通过翻耕后重新建植，另一方面可通过耕作措施进行改良。而深松浅旋是改良退化草地的有效措施之一，可以打破犁底层，疏松土壤，改善土壤通透性，增加土壤孔隙度，提高蓄水能力，能够促进苜蓿根系对养分和水分的吸收利用，提高草地生产性能[3-6]，实现退化紫花苜蓿草地更新复壮。因此，研究深松浅旋对退化紫花苜蓿草地改良具有重要意义。
近年来，我国为使耕地可持续利用，促进农业综合生产力稳定发展，实现“藏粮于地、藏粮于技”目标，提出了要全面推进农机深松作业[7-8]。关于深松浅旋的研究大多集中在耕作方式如何改良土壤理化性状[9-11]、改善作物根系特征[12-15]以及提高作物产量[16]等方面。赵亚丽等[17]、白伟等[18]认为深松可构造出虚实相间的耕层结构，降低土壤紧实度，促进作物根系生长发育，提高水分利用效率和作物产量。张凯等[19]研究深松程度对豫北农田土壤水分与作物耗水的影响时发现，深松深度为 40 cm 时效果最好，显著增加了作物的穗数和穗粒数。魏欢欢等[20]研究表明深松能够显著提高黄土高原东南部和西北部的冬小麦(Triticum aestivum)产量和水分利用效率。孔晓民等[21]研究发现，较免耕相比，深松可降低 0~35 cm 土壤容重和紧实度，尤其以 16~25 cm 土层降幅最为明显。也有研究证明，浅旋可提高土壤含水量、氮吸收量，减少硝态氮[22]，促进作物产量提升。基于以上综述，鲜见深松浅旋在退化紫花苜蓿草地上的报道，此外，前人的改良措施多以单一的深松或浅旋技术为主，只能在一定时间或一定区域内发挥优势。鉴于此，研究深松浅旋不同耕作组合对退化紫花苜蓿草地土壤理化性质、生产性能及饲草品质的影响，探讨该区域退化紫花苜蓿草地适宜的深松浅旋程度，可为改良退化苜蓿草地提供技术支撑。
1 材料与方法
1. 1 试验地概况
试验在宁夏隆德县神林乡观音村(35°21′N，106°15′E，海拔 2100 m)进行，该区域属中温带季风气候区。试验期间 2018、2019、2020 年降水量分别为 646、576、688 mm，但多集中在夏秋 7、8 两个月。年平均气温 5. 8 ℃，年均日照时数 2303. 5 h，无霜期 125 d。试验地为“陇东”苜蓿(M. sativa cv. Longdong)草地，2011 年条播，播种量 22. 5 kg·hm-2，行距为 15 cm，播种深度 2~3 cm;2018 年紫花苜蓿的覆盖度为 60%~70%。土壤为黑垆土，pH 值 7. 7，有机质含量 5. 88 g·kg-1，碱解氮含量 29. 75 mg·kg-1，速效磷含量 18. 89 mg·kg-1，速效钾含量 118. 08 mg·kg-1。
1. 2 试验设计
试验于 2018 年 4 月至 2020 年 9 月进行，试验地无施肥灌溉措施。试验采用单因素随机区组设计，设 7 个处理，分别为：深松 30 cm+浅旋 5 cm(S1Q1 )，深松 30 cm+浅旋 10 cm(S1Q2 )，深松 40 cm+浅旋 5 cm(S2Q1 )，深松 40 cm+浅旋 10 cm(S2Q2 )，深松 50 cm+浅旋 5 cm(S3Q1 )，深松 50 cm+浅旋 10 cm(S3Q2 )，不做处理(CK);每个处理 3 个重复，总计 21 个小区，小区面积 60 m(2 10 m×6 m)。2018 年 4 月 14 日顺着苜蓿行进行深松和浅旋作业。每年刈割两茬，时间分别为：2018 年 6 月 20 日和 9 月 2 日，2019 年 6 月 3 日和 8 月 22 日，2020 年 6 月 5 日和 8 月 25 日。
试验用具：青岛鲁耕310B 深松机;1GQN-150 立轴式浅旋机(将浅旋机刀片顶端部分弧度调整为 160°，以免对紫花苜蓿根颈造成过度伤害)。
1. 3 测定指标及方法
以下各农艺性状指标均在紫花苜蓿初花期测定：
株高：每个小区随机选取 20 株紫花苜蓿测量垂直高度，取平均值。
分枝数：每个小区随机取 20 株紫花苜蓿从根颈处，数其一级分枝数，取平均值。
叶茎比：每个小区随机选取 20 株紫花苜蓿将茎和叶分离，于 105 ℃杀青 1 h，65 ℃烘干 48 h 至恒重后冷却，称量干重，叶茎比=叶片干重/茎秆干重。
鲜草产量：在每个小区垂直行向取 6 m2 样方，留茬高度 4~5 cm，刈割后称量为鲜草产量，每个小区重复 6 次，换算为每 hm2 鲜草产量。
干草产量：将所采集的紫花苜蓿鲜草风干至恒重后称量干草重量，换算为每 hm2 干草产量。
营养品质测定：在每个小区内随机取 500 g 的鲜草样品，剪为 4~5 cm 长，于烘箱中 105 ℃下杀青 0. 5 h 左右， 65 ℃烘干 48 h 至恒重。将烘干的紫花苜蓿草样粉碎，过 0. 45 mm 筛。根据《饲料分析及饲料质量检测技术》[23]测定粗蛋白(crude protein，CP)、中性洗涤纤维(neutral detergent fiber，NDF)、酸性洗涤纤维(acid detergent fiber， ADF)含量。计算饲料相对饲喂价值(relative feeding value，RFV)[24]：
RFV = ( 88.9- 0.779× ADF )×( 120/NDF ) /1.29
土壤容重：采用环刀法测定，深度分别为 0~20 cm、20~40 cm、40~60 cm，逐层取土样，每个点 6 次重复;根据土壤容重，计算土壤孔隙度(soil porosity)：总孔隙度 = ( 1- 容重 ÷ 比重 )× 100% 式中：比重取 2. 65 g·cm-3[25]。
1. 4 数据统计分析
采用 Excel 2013 软件整理数据(干草产量为每年 2 茬累加，其余指标均为 2 茬平均)，利用 SPSS Statistics 25. 0 软件进行方差分析，多重比较，用 Origin 2018 作图，并用 SAS 9. 4 软件对株高(X1 )、分枝数(X2 )、叶茎比(X3 )、干草产量(X4 )、粗蛋白(X5 )、中性洗涤纤维(X6 )、酸性洗涤纤维(X7 )、相对饲喂价值(X8 )3 年平均值进行主成分分析(principal component analysis，PCA)。
2 结果与分析
2. 1 深松浅旋对紫花苜蓿株高和一级分枝数的影响
由图 1 可知，2018-2020 年紫花苜蓿株高的变化范围分别为 57. 34~75. 23 cm、54. 17~74. 03 cm、52. 00~ 72. 71 cm;相同深松浅旋方式下，除 S2Q2处理外，2018 年株高整体高于 2019 和 2020 年;其中 2018 年株高在 S3Q2处理下最高，达到 75. 23 cm，CK 处理最低，仅为 57. 34 cm;2019 和 2020 年株高均在 S2Q2 处理下最高，分别达到 74. 03 和 72. 71 cm，与 CK 相差 19. 86 和 22. 71 cm。一级分枝数在 2018-2020 年变化范围分别是 17~26 个、19~ 22 个、18~32 个;其中，2018 年一级分枝数最多的处理是 S1Q2，达到 26 个，显著高于其他处理(P<0. 05);2019 和 2020 年一级分枝数最多的均为 S2Q2处理，分别达到 22 和 32 个，较 CK 增加了 3 和 14 个。
2. 2 深松浅旋对紫花苜蓿叶茎比和干草产量的影响
由图 2 可知，2018-2020 年叶茎比变化范围分别为 0. 65~1. 03、0. 60~1. 05、0. 62~1. 03，且在这 3 年内，深松浅旋处理的叶茎比均显著高于 CK(P<0. 05);其中 2018 年叶茎比较高的处理分别是 S3Q2、S2Q1，较 CK 增加了 0. 38 和 0. 36。2019 和 2020 年叶茎比均在 S2Q2处理下达最大值，为 1. 05 和 1. 03，较 CK 提高了 0. 35 和 0. 31。干草产量在 2018-2020 年的变化范围分别是 1921. 26~6965. 03 kg·hm-2、1889. 79~6807. 12 kg·hm-2、1863. 13~ 8064. 37 kg·hm-2，整体来看，除 S3Q2处理外，2020 年干草产量高于 2018 和 2019 年;3 年深松浅旋处理的干草产量均显著高于 CK(P<0. 05)，其中 2018 年干草产量最高的处理是 S3Q2，其次是 S1Q2，分别为 6965. 03 和 5169. 18 kg·hm-2，2019 和 2020 年干草产量最高的处理均是 S2Q2，分别达到 6807. 12 和 8064. 37 kg·hm-2，而 CK 分别仅有 1889. 79 和 1863. 13 kg·hm-2。
2. 3 深松浅旋对紫花苜蓿粗蛋白(CP)含量和相对饲喂价值(RFV)的影响
由图 3 可知，CP 含量在 2018-2020 年的变化范围分别为 17. 23%~21. 83%、16. 55%~22. 24%、17. 69%~ 19. 65%，相同深松浅旋方式下，2020 年 CP 含量整体低于 2018 和 2019 年，且 2018 和 2019 年各深松浅旋处理的 CP 含量均显著高于 CK(P<0. 05)。其中 2018 年 CP 含量在 S3Q2处理下达最高，而 2019 和 2020 年 CP 含量均在 S2Q2处理下最高，分别达到 22. 64% 和 19. 95%，较 CK 分别提高了 34. 38% 和 9. 97%。2018-2020 年 RFV 的变化范围为 150. 92~189. 74、150. 72~175. 56、122. 19~160. 53，整体来看，2018 和 2019 年 RFV 高于 2020 年，其中 2018 年 RFV 最高的处理是 S1Q2，其次是 S3Q2，分别达到 189. 74、171. 33;2019 和 2020 年的 RFV 最高的处理均是 S2Q2，分别为 175. 56、160. 53，显著高于 CK(P<0. 05)。
2. 4 深松浅旋对紫花苜蓿中性洗涤纤维(NDF)和酸性洗涤纤维(ADF)含量的影响
由 图 4 可 知 ，NDF 含 量 在 2018-2020 年 的 变 化 范 围 分 别 为 32. 81%~39. 93%、35. 45%~39. 35%、 38. 79%~46. 87%，从 3 年数据来看，2018 年 NDF 含量整体低于 2019 和 2020 年;其中 2018 年 NDF 含量最低的处理是 S1Q2，显著低于其他处理(P<0. 05)，2019 年 NDF 含量最低的是 S2Q2处理，为 35. 45%，2020 年 NDF 含量在 S1Q2 处 理 下 最 低 ，与 CK 相 差 8. 08%。 ADF 含 量 在 2018-2020 年 的 变 化 范 围 分 别 是 27. 43%~31. 64%、 28. 25%~34. 45%、27. 79%~35. 09%;其中 2018 年 ADF 含量最低的是 S3Q2处理，为 27. 43%，显著低于除 S1Q2 处理外的其他处理(P<0. 05)，2019 和 2020 年 ADF 含量最低的均是 S2Q2处理，较 CK 分别低 5. 42% 和 7. 30%。
2. 5 深松浅旋对退化紫花苜蓿草地土壤容重的影响
由图 5 可知，适宜的深松浅旋处理可以显著降低土壤容重(P<0. 05)。整体来看，S2Q2处理的 0~20 cm 土层土壤容重逐年降低，而 S1Q2、S2Q1和 S3Q2处理的 0~20 cm 土层土壤容重均逐年增加，其中 2018-2020 年的土壤容重在 S2Q1处理下均达到最低，较 CK 分别降低 0. 24、0. 24、0. 19 g·cm-3。对于 20~40 cm 土壤容重而言，2018 年土壤容重在 S3Q1处理下最低，仅有 1. 36 g·cm-3。而 2019 和 2020 年土壤容重均在 S2Q1处理下最低，与 CK 分别相差 0. 10、0. 12 g·cm-3。2018、2019 和 2020 年的 40~60 cm 土层土壤容重均在 S3Q1处理下达到最低，分别为 1. 31、 1. 42、1. 45 g·cm-3，较 CK 降低了 0. 28、0. 09、0. 07 g·cm-3。
2. 6 深松浅旋对退化紫花苜蓿草地土壤孔隙度的影响
由图 6 可知，2018 和 2019 年 0~20 cm 土层土壤孔隙度均在 S2Q1处理下最高，分别达到 45. 31%、46. 42%; 2020 年 0~20 cm 土层在 S3Q1处理下最高，达到 46. 29%，较 CK 高 3. 23%。2018 年 20~40 cm 土层土壤孔隙度在 S3Q1处理下最高，达到 46. 89%;2019 和 2020 年 20~40 cm 土层均在 S2Q1处理下最高，达到 46. 29% 和 46. 13%，较CK 分别高 2. 27% 和 4. 24%。2018 年 40~60 cm 土层土壤孔隙度在 S2Q2处理下最高，达到 50. 26%;2019 和 2020 年 40~60 cm 土层均在 S3Q1处理下最高，分别达到 46. 29% 和 45. 38%，较 CK 分 别 高 3. 27% 和 2. 86%。
2. 7 综合分析
主成分分析能够降低数据维数以消除重叠信息的不利影响，近年来，PCA 综合评价方法在作物栽培以及牧草选育方面成为热点[26]。在对不同深松浅旋处理进行评价时，不能只考虑某 1 个或几个指标的优劣，而应该科学、综合的评价其所有指标。特征值和方差贡献率如表 1，根据特征值大于 1 原则，可提取 2 个主成分，累积方差贡献率达到 97. 171%，解释了总体信息的 97. 171%。
主成分对应的特征向量和载荷矩阵如表 2 所示，第 1 主成分主要综合了干草产量(X4 )、叶茎比(X3 )、粗蛋白(X5 )、酸性洗涤纤维(X7 )，其载荷较大，权重系数分别为 0. 978、0. 954、0. 951、-0. 883，可称为复合因子。第 2 主成分主要综合了中性洗涤纤维(X6 )、相对饲喂价值(X8 )，其载荷较大，权重系数分别为 0. 740、-0. 617。
3 讨论
3. 1 深松浅旋对紫花苜蓿生产性能的影响
深松浅旋可以不同程度提高草地生产力。本研究发现，深松浅旋对紫花苜蓿株高、一级分枝数有显著影响，主要是由于深松浅旋能够构建疏松的土层结构，降低土壤紧实度和土壤容重，增加土壤孔隙度，提高土壤持水能力，为紫花苜蓿生长发育提供良好的土壤环境，有利于根系吸收水分和养分[27-28]，进而增加紫花苜蓿植株高度和促进分枝。冯倩倩等[29]在耕作方式对冬小麦-夏玉米(Zea mays)产量影响的研究中指出，深松能够有效促进耕层土壤保水性能和提高有机质积累，增加作物的产量构成要素(有效穗数、穗粒数、千粒重)。李明等[30]研究表明，浅旋能够提高土壤含水量，有利于小麦根系吸收水分和养分，提高籽粒产量。本研究表明，深松浅旋处理较对照显著提高了紫花苜蓿产量，与前人研究结果一致，主要是因为深松提高了深层土壤苜蓿根系活性，从而延缓地上部分衰老，改善了紫花苜蓿叶片光合特性，进一步提高了单位面积产量[29，31]。同时，深松浅旋还可改善土壤理化性状，活化深层土壤资源，促进紫花苜蓿根系吸收养分和水分，最终直接影响紫花苜蓿产量的形成。从本研究结果来看，随着深松浅旋后年限的延长，大部分耕作处理的紫花苜蓿产量逐年增加，仅深松 50 cm 和浅旋 10 cm 的苜蓿产量有所下降，但也显著高于对照组，由此逐步显现了耕作改良的正效应，至于改良效果持久性如何，还需进一步研究。
3. 2 深松浅旋对紫花苜蓿品质的影响
牧草中的粗蛋白、中性洗涤纤维和酸性洗涤纤维含量是评定其营养品质的重要指标。本研究发现，深松浅旋可提高紫花苜蓿的粗蛋白含量，这是因为深松浅旋对紫花苜蓿生育期有不同程度的延长，能够延缓植物叶片衰老[32-34]，增加了紫花苜蓿叶片数量;从本试验的结果来分析，深松浅旋处理显著提高了紫花苜蓿叶茎比，而紫花苜蓿中 70% 的粗蛋白来源于叶片，且叶片中粗蛋白含量远远高于茎[35]，叶茎比越大，紫花苜蓿的粗蛋白含量就越高，因此处理组紫花苜蓿的粗蛋白含量高于对照组。深松浅旋处理第 1 年(2018 年)和第 2 年(2019 年)的紫花苜蓿 粗 蛋 白 含 量 整 体 高 于 第 3 年(2020 年)，原 因 可 能 是 在 退 化 苜 蓿 草 地 有 大 量 的 裸 露 斑 块 ，给 赖 草(Leymus secalinus)等杂草提供了生长空间，在 2018 年通过深松浅旋后，大部分地表杂草减少，地下根系受到破坏，生长缓慢，在紫花苜蓿收获时，只有极少数的杂草会和苜蓿一起被收割;到第 3 年(2020 年)，杂草数量增多，收割苜蓿时会有大量的杂草掺杂，降低了收获物粗蛋白含量。
牧草中性洗涤纤维和酸性洗涤纤维含量主要影响家畜的采食量和消化率，其主要是由纤维素、半纤维素和木质素组成[36]。本研究发现，中性洗涤纤维含量表现为 2020 年高于 2018 与 2019 年，主要原因是 2020 年紫花苜蓿叶茎比低于前两年，且紫花苜蓿叶片中的中性洗涤纤维、木质素和半纤维素含量均低于茎中[37]，因此，紫花苜蓿中性洗涤纤维含量高。
3. 3 深松浅旋对土壤容重和孔隙度的影响
深松浅旋对退化紫花苜蓿草地土壤容重起到至关重要的作用。研究发现，科学合理的深松浅旋能够显著降低土壤容重，增加土壤孔隙度，促进作物生长。李小飞等[38]研究表明，深松可不同程度降低茶(Camellia sinensis)园土壤容重，尤其在深松 30 cm 时 0~30 cm 土层降幅最为明显。梁金凤等[39]研究发现，深松深度为 30、35、45 cm 时均可降低土壤容重。张凯等[19]以豫北农田为研究对象，试验发现深松深度为 35 和 40 cm 处理的土壤容重较对照显著降低。本研究表明，深松浅旋降低了 0~60 cm 土层的土壤容重，与前人研究结果一致。深松有效打破了土壤犁底层，进而降低土壤容重，改善退化苜蓿草地土壤结构。另外，在深松浅旋第 1 年(2018 年)，S1Q1、S1Q2、 S2Q1、S3Q1和 S3Q2处理下 0~20 cm 土层土壤容重整体低于 2019 和 2020 年，主要原因可能是深松处理后又对土壤进行了浅旋处理，浅层土壤受到扰动，土壤结构稳定性下降，土壤容重降低[40]。
土壤孔隙结构与土壤通气和持水性能紧密相关，土壤孔隙度越大越有利于植物根系伸展和吸收养分。 Blazewiczwozniak 等[41]研究发现，深松浅旋能降低土壤紧实度，增加土壤孔隙度。刘卫玲等[42]为探明深松方式对砂姜黑土耕层改良的影响时发现，深松可以改善土壤孔隙结构，增加孔隙体积和大孔隙数量。本研究发现，不同程度深松浅旋较对照显著提高了土壤孔隙度，与前人研究结果一致，深松浅旋打破了土壤耕作层，增加了土壤孔隙度。 ——论文作者：王斌 1 ，李满有 1 ，王欣盼 1 ，董秀 2 ，庞军宝 3 ，兰剑 1*
