深松耕对石灰岩红壤物理性状和玉米生长发育的影响
Effect of Sub-soiling on the Physical Properties of Limestone Red Soil and Maize Growth
-
Keywords:
- red soil /
- sloping land /
- sub-soiling /
- soil physical properties /
- maize
-
细胞分类学(Cytotaxonomy)是20世纪初在细胞形态学的基础上发展起来的一门新兴学科,其实质是将染色体的数目、形态结构、行为、DNA含量以及核苷酸重复序列等细胞学资料运用在分类学上,进而探讨和研究类群之间的亲缘关系、分类地位、物种形成机制和进化关系等问题[1-5]。在细胞分类学的基础上,通过与其他学科的交叉融合,还发展出了物种生物学(Biosystematics)、细胞地理学(Cytogeography)等衍生学科[3, 6]。中国西南地区植物类群丰富,且成分复杂,仅云南就具有几乎所有的分布型和亚型,特有现象明显[7];细胞分类学可以解决常规形态分类难以解决的问题,是现代综合分类的重要组成部分[8],对于该地区的植物资源利用及系统分类学具有重要意义。
蝇子草属(Silene L.)为石竹科(Caryophyllaceae)一年或多年生草本植物,全球约400种。中国有112种2亚种17变种,广布长江流域和北部各省区,以西北和西南地区较多[9]。横断山地区共有49种18个种下单位[10]。在国外,蝇子草属植物的染色体数目有较多研究,且大部分种已开展细胞学研究[11-12]。但是在西南横断山地区,截至目前有染色体报道的蝇子草属植物只有7种,分别是细蝇子草(S. gracilicaulis,2n=2x=24=22m+2sm)、垫状蝇子草(S. davidii,2n=4x=48=48m)、尼泊尔蝇子草(S. nepalensis,2n=2x=24=16m+8sm)、变黑蝇子草(S. nigrescens,2n=4x=48=40m+8sm)、内蒙古女娄菜(S. orientalimongolica,2n=2x=24=24m)、糙叶蝇子草(S. trachyphylla,2n=2x=24=8m+6sm+10st)和腺毛蝇子草(S. yetii,2n=6x=72=64m+8sm)[13]。对于该地区染色体的报道十分有限,需要更多的数据补充。
本研究对3种蝇子草属植物进行核型的首次报道,分别为(1)喜马拉雅蝇子草(S. himalayensis),分布于中国河北、湖北、陕西、四川、云南和西藏等省区,以及印度锡金地区,生于海拔2 000~5 000 m的灌丛间或高山草甸;(2)喇嘛蝇子草(S. lamarum)为多年生草本,分布于四川(盐源、西昌、雷波、木里)和云南(德钦、中甸),生于海拔3 400~4 000 m的高山草地或灌丛中;(3)墨脱蝇子草(S. namlaensis),为多年生草本,高30~40 cm,叶片狭椭圆形或披针形,花单生或2~3朵,微俯垂,后期直立蒴果卵形;花期8—9月。墨脱蝇子草分布于中国西藏(林芝、波密、墨脱、加查、措那、聂木拉),生长在海拔3 600~4 500 m的高山草地[9]。
1. 材料与方法
1.1 材料
本试验供试材料的采集具体地点见表1。凭证标本存放在中国科学院昆明植物研究所标本馆(KUN)及西南林业大学林学院标本室,永久装片保存于云南师范大学生命科学学院植物细胞学实验室。
表 1 蝇子草属3种植物材料的来源信息Table 1. The original information of three species of Silene种名species 采集地location 海拔/m altitude 经纬度position 标本号sample No. 喜马拉雅蝇子草 S. himalayensis 四川稻城Daocheng, Sichuan 4 420 E100°16′42.25′′, N28°52′54.96′′ MS17-689 喇嘛蝇子草 S. lamarum 云南昆明Kunming, Yunnan 4 101 E102°55′40.09′′, N26°09′21.48′′ MSC-058 墨脱蝇子草 S. namlaensis 西藏林芝Linzhi, Tibet 4 578 E94°38′50.47′′, N29°36′59.08′′ Tsui-645 1.2 方法
把野外收集的种子放置于4 ℃冰箱,然后将种子置于固体培养基,24 ℃恒温培养箱萌发,待其根尖长到1~2 cm时取下,用0.003 mol/L的8-羟基喹啉常温下避光处理3 h,之后用蒸馏水清洗2~3次,用卡诺氏固定液(无水乙醇∶冰乙酸= 3∶1)在4 ℃冰箱固定24 h,用蒸馏水清洗之后,加入1 mol/L的HCl,在60 ℃ 水浴锅中解离10~15 min,随后清洗,用改良苯酚品红染液染色24 h,使用常规压片法镜检观察,挑选处于分裂中期的分散良好的细胞拍照并制作永久装片。
染色体核型分析标准按照李懋学等[14]的标准,按照ARANO[15]的方法计算核型不对称系数:As.K =(染色体组长臂总和/染色体组总长度和)× 100%;参考PASZKO[16]、PERUZZI等[17]和ASTUTI等[18]对核型参数使用的建议,得出相关核型参数及核型公式;最后取5~6个染色体分散良好的细胞照片进行数据分析。
2. 结果与分析
2.1 喜马拉雅蝇子草[S. himalayensis (Rohrb.) Majumdar]
喜马拉雅蝇子草染色体数目为2n=48,核型公式为2n=4x=48m,为四倍体,核型类型为1A型;染色体长度范围为3.78~2.24 μm,臂比范围为1.51~1.26 μm,核型不对称指数(As.K)为56.85%,臂比>2∶1的为0,染色体内不对称程度(MCA)为13.68。染色体参数如表2、3所示,染色体核型模式图如图1所示。
表 2 蝇子草属3个种的染色体参数Table 2. Chromosome parameters of three species of Silene种名 species 染色体序号 chromosome No. 长臂 L 短臂 S 总长 L+S 臂比 L/S 染色体类型 chromosome type 喜马拉雅蝇子草
S. himalayensis1 2.15 1.41 3.55 1.54 m 2 1.96 1.59 3.55 1.25 m 3 1.75 1.41 3.16 1.24 m 4 1.86 1.25 3.11 1.50 m 5 1.66 1.29 2.94 1.29 m 6 1.67 1.26 2.93 1.33 m 7 1.48 1.27 2.75 1.17 m 8 1.56 1.19 2.74 1.30 m 9 1.44 1.19 2.62 1.21 m 10 1.46 1.03 2.49 1.41 m 11 1.40 1.02 2.41 1.38 m 12 1.30 1.03 2.33 1.27 m 喇嘛蝇子草
S. lamarum1 3.40 2.91 6.31 1.17 m 2 3.41 2.45 5.86 1.39 m 3 3.03 2.32 5.35 1.31 m 4 2.95 2.23 5.18 1.32 m 5 2.71 2.20 4.91 1.23 m 6 3.01 1.84 4.84 1.64 m 7 2.45 2.15 4.61 1.14 m 8 2.82 1.84 4.66 1.53 m 9 2.97 1.40 4.37 2.12 m 10 2.53 1.80 4.33 1.41 m 11 2.96 1.11 4.07 2.67 m 12 2.33 1.75 4.08 1.33 sm 墨脱蝇子草
S. namlaensis1 1.67 1.24 2.90 1.34 m 2 1.63 1.36 3.00 1.20 m 3 1.56 1.26 2.81 1.24 m 4 1.62 1.03 2.66 1.58 m 5 1.43 1.17 2.59 1.22 m 6 1.50 1.07 2.57 1.40 m 7 1.43 1.00 2.43 1.42 m 8 1.29 1.11 2.40 1.16 m 9 1.24 1.02 2.27 1.22 m 10 1.30 0.9 2.19 1.44 m 11 1.14 0.97 1.12 1.17 m 12 1.1 0.81 1.91 1.36 m 注:L. 长臂;S. 短臂;L+S. 染色体长度;L/S. 臂比;m. 中部着丝粒区;sm.近中部着丝粒区。
Note: L. long arm; S. short arm; L+S. total length; L/S. arms ratio; m. median centeomere chromosome; sm. submedian centromere chromosome.表 3 蝇子草属3个种的染色体核型参数Table 3. Karyotype and chromosomal data of three species of Silene种名species 染色体数目/基数/倍性
chromosome No./basic No./polyploidy核型karyotype CVCL CVCI AI THL MCA As.K/% 核型公式
karyotype formula喜马拉雅蝇子草
S. himalayensis48/12/4x 1A 14.16 5.06 0.72 34.58 13.68 56.94 2n=4x=48m 喇嘛蝇子草
S. lamarum24/12/2x 2A 14.21 14.18 2.01 58.57 18.85 59.02 2n=2x=19m+5sm 墨脱蝇子草
S. namlaensis72/12/6x 1A 13.39 5.79 0.76 29.85 13.23 56.65 2n=6x=72m 注:CVCL. 染色体长度变异系数; CVCI. 着丝粒指数变异系数; AI. 染色体不对称指数; THL. 核型单倍体染色体总长; MCA. 染色体内不对称程度; As.K. 核型不对称指数。
Note:CVCL. coefficient of variation of chromosome length; CVCI. coefficient of variation of centromere index; AI. chromosome asymmetry index; THL. total chromosome length of karyotype haploid; MCA. degree of intrachromosomal asymmetry; As.K. karyotype asymmetry index.图 1 3种蝇子草的核型细胞及模式图注:A1~A5. 喇嘛蝇子草;B1~B5. 喜马拉雅蝇子草;C1~C5. 墨脱蝇子草;1~5分别代表间期核、前期核、中期核、核型模式图和核型图。Figure 1. The metaphase chromosomes and idiograms of SileneNote: A1-A5. S. lamarum; B1-B5. S. himalayensis; C1-C5. S. namlaensis; 1-5 interphase nucleus, prophase nucleus, metaphase nucleus, idiogram and karyogram.2.2 喇嘛蝇子草(S. lamarum C. Y. Wu)
喇嘛蝇子草染色体数目为2n=24,为二倍体,核型公式为2n=2x=19 m+5 sm,核型类型为2A型;染色体长度范围为3.4~2.33 μm,臂比范围为1.14~2.67 μm,核型不对称指数(As.K)为59.02%,臂比>2∶1的占16.7%,染色体内不对称程度(MCA)为18.85。染色体参数如表2、3所示,染色体核型模式图如图1所示。
2.3 墨脱蝇子草[S. namlaensis (Marquand) Bocquet]
墨脱蝇子草染色体数目为2n=72,核型公式为2n=6x=72,核型类型为1A型,染色体长度范围为2.9~1.91 μm,臂比范围为1.58~1.16 μm, 核型不对称指数(As.K)为56.65%,染色体内不对称程度(MCA)为13.23。染色体参数如表2、3所示,染色体核型模式图如图1所示。
3. 讨论
本研究首次报道了蝇子草属3个种的核型及染色体数目,3个种的间期核形态都相似,染色质比较均匀地分布在细胞核内,周围有染色较浅的染色质丝和纤维状结构,根据TANAKA[19]的间期核划分标准为扩散型(diffuse type);有丝分裂前期核染色体上着色较深且比较均匀,隶属于中间型前期核(interstitial type)。有丝分裂中期的染色体数目为:喇嘛蝇子草(S. lamarum)为24条,是二倍体;喜马拉雅蝇子草(S. himalayensis)为48条,是四倍体;墨脱蝇子草(S. namlaensis)为72条,是六倍体。3个种的染色体基数都是12,但是在倍性上存在不同,分别是二倍、四倍和六倍,这可能是由于它们生长的地理位置导致的倍性变化,物种为了适应恶劣的环境,改变自身的倍性,产生多倍化现象,这也是物种分化变异的主要来源之一。在本研究中,墨脱蝇子草生长的海拔总体较高,其次是喜马拉雅蝇子草、喇嘛蝇子草,从这个角度来看,随着海拔的升高,倍性也在增加,但是由于试验数据有限,更普遍的结论需要更进一步的试验验证。
从核型类型来看,喜马拉雅蝇子草和墨脱蝇子草核型都属于1A型,喇嘛蝇子草则为2A型,根据STEBBINS[20]的理论,在生物的进化过程中,染色体核型是由对称性向非对称性演化的。理论认为:核型对称程度越高的生物,其染色体变异越小,进化程度也越低;而非对称程度越高的生物,其染色体变异越大,进化程度越高。按照这一理论及本试验结果可推测:在系统演化上,喜马拉雅蝇子草和墨脱蝇子草属于进化程度相对较低的种类,而喇嘛蝇子草的染色体则属于进化程度比较高的类型。本研究结果可为蝇子草属染色体提供更多的数据补充,并为蝇子草属的系统分类学及探讨多倍化现象、多倍化物种形成等研究提供有价值的参考。
致谢:感谢西南林业大学林学院徐波、中国科学院昆明植物研究所马祥光和孙文光提供试验种子以及云南师范大学生命科学学院苏富明、饶培瑜和胡晶晶在试验上的帮助。
-
表 1 深松耕对玉米产量性状的影响
Table 1 Effect of sub-soiling on the yield traits of maize
年份 year 处理 treatment 穗长/cm ear length 穗粗/cm ear diameter 秃尖长/cm length of bald 穗粒数 grain number 百粒质量/g 100-kernel weight 2016 旋耕 rotary tillage 18.63±0.71 a 16.12±0.19 a 3.42±0.62 a 502.24±41.74 a 28.77±2.17 a 深松耕 sub-soiling 19.90±0.32 b 16.51±0.22 b 2.44±1.37 a 552.54±34.20 a 31.27±1.61 a 2017 旋耕 rotary tillage 22.08±2.33 a 16.79±0.81 a 3.53±0.53 a 538.80±39.76 a 40.98±3.94 a 深松耕 sub-soiling 22.38±0.87 a 16.72±0.49 a 2.67±1.37 a 562.54±34.20 a 39.65±1.78 a 2018 旋耕 rotary tillage 16.92±1.74 a 15.27±0.61 a 0.61±0.29 a 529.90±89.48 b 31.17±2.82 a 深松耕 sub-soiling 16.20±2.03 a 15.28±0.65 a 1.12±0.80 b 474.15±92.61 a 33.33±3.49 b 注:不同字母代表差异显著(P<0.05);下同。
Note: Different letters represent significant diffences (P<0.05); the same below.表 2 深松耕对大喇叭口时期玉米生长的影响
Table 2 Effect of sub-soiling on the maize growth at the big flare stage
年份 year 处理 treatment 株高/cm plant height 茎粗/cm stem diameter 叶面积指数 leaf area index 2016 旋耕 rotary tillage 187.30±13.26 a 8.22±0.84 a 3.75±0.24 a 深松耕 sub-soiling 193.55±16.89 a 8.38±0.15 a 3.80±0.14 a 2017 旋耕 rotary tillage 210.06±5.75 a 7.84±0.39 a 3.95±0.13 a 深松耕 sub-soiling 212.85±11.22 a 7.95±0.27 a 4.18±0.13 b 2018 旋耕 rotary tillage 190.33±22.35 a 8.39±0.91 a 3.30±0.34 b 深松耕 sub-soiling 175.80±21.77 a 8.06±0.72 a 3.08±0.45 a 表 3 深松耕对收获期玉米生长的影响
Table 3 Effect of sub-soiling on the maize growth at the harvest time
年份 year 处理 treatment 株高/cm plant height 茎粗/cm stem diameter 穗位高/cm ear height 2016 旋耕 rotary tillage 228.10±9.30 a 6.44±0.39 a 67.85±5.17 a 深松耕 sub-soiling 236.58±9.74 a 6.78±0.45 a 68.15±3.45 a 2017 旋耕 rotary tillage 250.85±3.01 a 7.84±0.15 a 93.20±1.60 a 深松耕 sub-soiling 255.30±10.07 a 7.95±0.27 a 91.10±1.80 a 2018 旋耕 rotary tillage 217.90±15.57 a 6.73±0.93 a 73.15±12.21 a 深松耕 sub-soiling 210.38±18.51 a 6.75±0.79 a 70.10±10.30 a 表 4 深松耕对土壤容重的影响
Table 4 Effect of sub-soiling on the soil bulk density
年份 year 处理 treatment 土壤容重/(g·cm−3) soil bulk density 0~10 cm >10~20 cm >20~30 cm 2016 旋耕 rotary tillage 1.28±0.07 b 1.26±0.05 a 1.34±0.11 a 深松耕 sub-soiling 1.15±0.06 a 1.22±0.17 a 1.31±0.10 a 2017 旋耕 rotary tillage 1.21±0.01 b 1.33±0.05 b 1.49±0.11 b 深松耕 sub-soiling 1.12±0.02 a 1.19±0.03 a 1.31±0.02 a 2018 旋耕 rotary tillage 1.31±0.09 a 1.43±0.07 a 1.56±0.06 a 深松耕 sub-soiling 1.26±0.07 a 1.38±0.11 a 1.45±0.15 a 表 5 深松耕对土壤含水量的影响
Table 5 Effect of sub-soiling on the soil moisture content
年份 year 处理 treatment 土壤含水量/% soil moisture content 0~10 cm >10~20 cm >20~30 cm 2016 旋耕 rotary tillage 22.18±1.34 a 21.68±2.10 a 21.76±1.14 a 深松耕 sub-soiling 24.57±2.84 a 24.44±1.19 a 26.58±2.38 a 2017 旋耕 rotary tillage 24.71±0.83 a 23.94±0.94 a 21.37±0.72 a 深松耕 sub-soiling 27.26±1.97 a 28.42±1.51 b 27.84±1.96 b 2018 旋耕 rotary tillage 20.62±2.05 a 20.52±0.99 a 21.00±0.88 a 深松耕 sub-soiling 23.13±2.84 a 24.00±0.65 b 24.50±0.62 b 表 6 深松耕对土壤总孔隙度的影响
Table 6 Effect of sub-soiling on the soil total porosity
年份 year 处理 treatment 土壤总孔隙度/% soil total porosity 0~10 cm >10~20 cm >20~30 cm 2016 旋耕 rotary tillage 51.55±2.61 a 52.48±1.71 a 49.28±4.14 a 深松耕 sub-soiling 56.47±2.27 b 53.78±6.32 a 50.75±3.70 a 2017 旋耕 rotary tillage 54.30±0.43 a 49.79±1.72 a 43.78±4.20 a 深松耕 sub-soiling 57.88±0.61 b 55.19±1.12 b 50.39±0.90 b 2018 旋耕 rotary tillage 50.54±3.34 a 45.98±2.75 a 41.07±2.24 a 深松耕 sub-soiling 52.38±2.78 a 47.85±4.06 a 45.12±5.60 a 表 7 深松耕对土壤紧实度的影响
Table 7 Effect of sub-soiling on the soil compactness
年份 year 处理 treatment 土壤紧实度/kPa soil compactness 0~10 cm >10~20 cm 2016 旋耕 rotary tillage 561.54±54.09 a 795.80±48.23 b 深松耕 sub-soiling 496.08±47.40 a 694.17±22.74 a 2017 旋耕 rotary tillage 544.31±28.66 a 666.61±34.38 a 深松耕 sub-soiling 500.70±47.33 a 626.99±23.22 a 2018 旋耕 rotary tillage 881.92±149.03 a 1422.79±237.57 a 深松耕 sub-soiling 961.50±264.37 a 1528.34±335.20 a -
[1] 胡钧铭, 陈胜男, 韦翔华, 等. 耕作对健康耕层结构的影响及发展趋势[J]. 农业资源与环境学报, 2018, 35(2): 95. DOI: 10.13254/j.jare.2017.0242. [2] YANG X, PANG H C, LI Y B, et al. Effects of deep rotary sub-soiling tillage on the physical properties and crop growth of the sticky loamy soil in North China[J]. Scientia Agricultura Sinica, 2013, 46(16): 3401. DOI: 10.3864/j.issn.0578-1752.2013.16.011.
[3] SENE M, VEPRASKAS M J, NADERMAN G C, et al. Relationships of soil texture and structure to corn yield response to subsoiling[J]. Soil Science Society of America Journal, 1984, 49(2): 422. DOI: 10.2136/sssaj1985.03615995004900020030x.
[4] 王万宁, 强小嫚, 刘浩, 等. 麦前深松对夏玉米土壤物理性状和生长特性的影响[J]. 水土保持学报, 2017, 31(6): 229. DOI: 10.13870/j.cnki.stbcxb.2017.06.036. [5] 薛永杰, 姚影, 赵文生, 等. 通辽地区玉米机械化深松改土效果分析[J]. 北方农业学报, 2017, 45(4): 6. DOI: 10.3969/j.issn.2096-1197.2017.04.02. [6] 郭家萌, 刘振朝, 高强, 等. 深松对玉米产量和养分吸收的影响[J]. 水土保持学报, 2016, 30(2): 249. DOI: 10.13870/j.cnki.stbcxb.2016.02.043. [7] 张丽娜, ASENSO E, 张陆勇, 等. 耕作方式对旱地红壤物理特性的影响[J]. 水土保持研究, 2018, 25(3): 46. [8] 温美娟, 王成宝, 霍琳, 等. 深松和秸秆还田对甘肃引黄灌区土壤物理性状和玉米生产的影响[J]. 应用生态学报, 2019, 30(1): 224. DOI: 10.13287/j.1001-9332.201901.030. [9] 吕巨智, 闫飞燕, 程伟东, 等. 不同耕作方式对土壤理化性状及玉米产量的影响[J]. 江苏农业科学, 2015, 43(11): 118. DOI: 10.15889/j.issn.1002-1302.2015.11.035. [10] 王岩, 刘玉华, 张立峰, 等. 耕作方式对冀西北栗钙土土壤物理性状及莜麦生长的影响[J]. 农业工程学报, 2014, 30(4): 109. DOI: 10.3969/j.issn.1002-6819.2014.04.014. [11] 曾维军, 张建生, 刘淑霞, 等. 云南省耕地自然质量等流域空间分布特征研究[J]. 水土保持研究, 2016, 23(5): 267. DOI: 10.13869/j.cnki.rswc.2016.05.034. [12] 郎凤莲, 张晓云, 李永贤, 等. 耕作方式对红壤坡地土壤物理性质和玉米产量的影响[J]. 云南农业大学学报(自然科学), 2019, 34(3): 377. DOI: 10.12101/j.issn.1004-390X(n).201804026. [13] HE J, MCHUGH A D, LI H W, et al. Permanent raised beds improved soil structure and yield of spring wheat in arid North-western China[J]. Soil Use and Management, 2012, 28(4): 536. DOI: 10.1111/j.1475-2743.2012.00445.x.
[14] 汤文光, 肖小平, 唐海明, 等. 长期不同耕作与秸秆还田对土壤养分库容及重金属Cd的影响[J]. 应用生态学报, 2015, 26(1): 168. DOI: 10.13287/j.1001-9332.20141029.011. [15] 张合云, 汤丽芬. 土地深耕深松机械化技术推广[J]. 云南农业, 2016(8): 87. DOI: 10.3969/j.issn.1005-1627.2016.08.046. [16] 张晓云, 郎凤莲, 李永贤, 等. 坡地玉米马铃薯间作群体地上部垂直层化对水土流失的影响[J]. 云南农业大学学报(自然科学), 2017, 32(5): 904. DOI: 10.16211/j.issn.1004-390X(n).2017.05.023. [17] 鲍士旦. 土壤农化分析[M]. 3版. 北京: 中国农业出版社, 2000. [18] LI X, TANG M J, ZHANG D X, et al. Effects of sub-soiling on soil physical quality and corn yield[J]. Transactions of the Chinese Society of Agricultural Engineering, 2014, 30(23): 65. DOI: 10.3969/j.issn.1002-6819.2014.23.009.
[19] 周华, 张居敏, 祝英豪, 等. 秸秆还田深松旋埋联合耕整机设计与试验[J]. 农业工程学报, 2017, 33(22): 17. DOI: 10.11975/j.issn.1002-6819.2017.22.003. [20] 王宪良, 王庆杰, 李洪文, 等. 免耕条件下轮胎压实对土壤物理特性和作物根系的影响[J]. 农业机械学报, 2017, 48(6): 168. DOI: 10.6041/j.issn.1000-1298.2017.06.022. [21] HE J, LI H W, WANG X Y, et al. The adoption of annual subsoiling as conservation tillage in dryland maize and wheat cultivation in Northern China[J]. Soil and Tillage Research, 2006, 94(2): 493. DOI: 10.1016/j.still.2006.10.005.
[22] TIAN S Z, WANG Y, NING T Y, et al. Continued no-till and subsoiling improved soil organic carbon and soil aggregation levels[J]. Agronomy Journal, 2014, 106(1): 212. DOI: 10.2134/agronj2013.0288.
[23] 杨志忠. 机械深松深耕增产技术及机理探讨[J]. 当代农机, 2014(6): 65. DOI: 10.3969/J.ISSN.1673-632X.2014.06.08. [24] 杨艳. 不同耕作措施对农田土壤理化性质和作物产量的影响[D]. 杨凌: 西北农林科技大学, 2017. [25] 郑侃, 何进, 李洪文, 等. 中国北方地区深松对小麦玉米产量影响的Meta分析[J]. 农业工程学报, 2015, 31(22): 7. DOI: 10.11975/j.issn.1002-6819.2015.22.002. [26] 何进, 李洪文, 高焕文. 中国北方保护性耕作条件下深松效应与经济效益研究[J]. 农业工程学报, 2006, 22(10): 62. DOI: 10.3321/j.issn:1002-6819.2006.10.013. [27] LOPEZ F C, DORADO J, PARDO M T. Effects of zone-tillage in rotation with no-tillage on soil properties and crop yields in a semi-arid soil from central Spain[J]. Soil and Tillage Research, 2007, 95(1): 266. DOI: 10.1016/j.still.2007.01.005.
[28] SONG S Y, JIN Q M, ZHANG W. Effects of occur tendency of maize diseases by pathogenetic fungi distribution in subsoiling fields[J]. Journal of Maize Sciences, 2007, 15(3): 130. DOI: 10.3969/j.issn.1005-0906.2007.03.032.