Effects of intercropping soybeans with different root architecture and sweet maize on crop growth and phosphorus uptake
-
摘要:目的
研究不同磷水平下不同根构型大豆与甜玉米间作对作物磷吸收与生长的影响，探究间作体系内根系形态、构型及根际土壤磷有效性等的关系。
方法田间试验于2022年8月在广东省广州市增城区华南农业大学教学试验基地进行。以2个不同根构型的大豆品种‘本地2号’(深根型)和‘粤春03-3’(浅根型)为材料，采用裂区设计，主区为施磷(+P：大豆40 kg·hm−2、甜玉米120 kg· hm−2)和不施磷(−P：大豆0 kg·hm−2、甜玉米0 kg·hm−2) 2种施磷水平，副区为甜玉米单作、‘本地2号’‖甜玉米、‘粤春03-3’‖甜玉米、‘本地2号’单作和‘粤春03-3’单作5种种植模式。测定甜玉米和大豆的产量、生物量、磷吸收量以及根系形态和构型的相关指标，计算间作系统的土地当量比和种间竞争力。
结果无论施磷还是不施磷条件下，间作大豆显著增加甜玉米产量，平均增产33.4%；不施磷条件下，间作大豆显著增加甜玉米生物量，平均增加62.7%。甜玉米和大豆间作具有间作优势，且受施磷水平影响；施磷条件下甜玉米与大豆间作的土地当量比为1.08，而不施磷条件下为1.21。种间竞争力分析表明，间作体系中，甜玉米竞争力显著强于大豆(种间竞争力>0)，这种竞争优势在不施磷条件下更明显。此外，不施磷条件下间作深根型大豆‘本地2号’显著促进甜玉米磷吸收，平均增加40.6%。进一步分析发现，不施磷间作改变深根型大豆细根(直径≤0.5 mm)占比，同时诱导其根系拓宽，并显著提升甜玉米根际土壤有效磷含量。对间作体系内大豆和甜玉米性状的相关性分析发现，甜玉米的生物量与其总根长、大豆根宽和大豆根际土壤有效磷含量呈显著正相关；大豆的生物量与其磷吸收量、总根长、根宽以及根际土壤有效磷含量呈显著正相关，与其细根占比(直径≤0.5 mm)呈显著负相关。
结论不施磷条件下，间作深根型大豆可有效促进甜玉米生长与磷吸收，提高甜玉米产量。研究结果为充分挖掘甜玉米‖大豆间作系统磷素利用潜力、筛选适合间作的大豆品种、实现磷肥减施增效提供了重要理论依据。
Abstract:ObjectiveTo study the effects of intercropping soybeans with different root architecture and sweet maize on crop phosphorus (P) uptake and growth under different P levels, and explore the relationship among root morphology, root architecture, P availability in rhizospheric soil and so on within the intercropping system.
MethodField experiments were conducted in August 2022 at the Teaching and Experimental Base of South China Agricultural University in Zengcheng District, Guangzhou City, Guangdong Province. Two different soybean cultivars, namely ‘Bendi 2’ (deep-rooted genotype) and ‘Yuechun 03-3’ (shallow-rooted genotype) were selected as experimental materials. The field experiment employed a split-plot design, with the main plots comprising two P application levels: With P (+P: 40 kg·hm−2 for soybean, 120 kg·hm−2 for sweet maize) and without P (−P: 0 kg·hm−2 for both soybean and sweet maize). The subplots consisted of five different planting patterns as follows: Sweet maize monocropping, soybean ‘Bendi 2’ ‖ sweet maize, soybean ‘Yuechun 03-3’ ‖ sweet maize, soybean ‘Bendi 2’ monocropping and soybean ‘Yuechun 03-3’ monocropping. The yield, biomass, P uptake, root morphology and architecture of sweet maize and soybean were measured, and the land equivalent ratio (LER) and interspecific competitiveness of the intercropping system were calculated.
ResultIntercropping soybean significantly increased the yield of sweet maize by average of 33.4% under both P levels. Under −P condition, intercropping soybean significantly increased the biomass of sweet maize by average of 62.7%. Intercropping sweet maize and soybean exhibited advantages, influencing by P application levels. Under +P condition, the average LER for intercropped sweet maize and soybean was 1.08, while under −P condition, the average LER increased to 1.21. Interspecific competition analysis indicated that in the intercropping system, sweet maize showed significantly stronger competitiveness than soybean (interspecitic competitiveness > 0). This competitive advantage was more pronounced under −P condition. In addition, intercropping with deep-rooted soybean genotype (‘Bendi 2’) significantly promoted P uptake in sweet maize under −P condition by average of 40.6%. Further analysis revealed that for deep-rooted soybean ‘Bendi 2’, intercropping with sweet maize under −P condition widened its root width, altered the proportion of fine roots (diameter ≤ 0.5 mm), and significantly increased P availability in the rhizosphere of intercropping maize. The correlation analysis of traits of soybean and sweet maize in the intercropping system revealed a significant positive correlation between the biomass of sweet maize and its total root length, root width of soybean, and available P content in the soybean rhizosphere soil. The biomass of soybean showed a significant positive correlation with its P uptake, total root length, root width, and available P content in the rhizosphere soil. However, it exhibited a significant negative correlation with its proportion of fine roots (diameter ≤ 0.5 mm).
ConclusionIntercropping with deep rooting soybean under −P condition can effectively promote the growth and P absorption of sweet maize, and increase its yield. These findings provide important theoretical basis for fully exploring the P utilization potential in the sweet maize‖soybean intercropping system, screening suitable soybean varieties for crop intercropping, and achieving reduced P fertilization rate with increased P use efficiency.
-
Keywords:
- Phosphorus level /
- Intercropping /
- Soybean /
- Sweet maize /
- Root configuration /
- Phosphorus uptake
-
研究发现，马铃薯70%的块茎损伤出现在收获环节[1]，机械损伤严重影响了马铃薯的品质、产量和经济效益[2-5]。随着马铃薯收获机械化、加工自动化的快速发展，马铃薯由碰撞导致的损伤问题亟待解决。有研究根据马铃薯果肉的变色来度量其损伤程度，提出了损伤的评价标准[6-8]。吴永根等[9]将经冲击的果实在自然环境中放置24 h，损伤的果肉会呈现褐色，从而获得损伤体积、损伤深度等试验指标。桑永英等[10]根据淀粉遇碘变色原理，对碰撞后的果实涂抹碘，根据颜色深浅划分损伤程度。康璟等[11]提出通过多因素试验研究马铃薯机械损伤的建议，如从挖掘深度、下落高度等角度分析马铃薯的机械损伤，从应力应变及其临界极限值方面研究马铃薯的机械损伤等。王咏梅等[12]研究了马铃薯收获中的机械损伤，分析了机械损伤产生的过程和机理，建议对马铃薯动态力学特性进行深入研究。李辉等[13]通过对比国内外的马铃薯机械损伤研究，鼓励多从实际作业中发现问题、解决问题。魏忠彩等[14]系统研究了收获、清选和分级对马铃薯机械损伤的影响。冯斌等[15]通过试验及参数的数学关系获得马铃薯块茎碰撞恢复系数，并分析跌落高度、碰撞材料等跌落过程的参数与块茎碰撞恢复系数的关系。吕金庆等[16]研究了跌落高度、升运链倾角和线速度对马铃薯损伤综合指数和伤薯率的影响。Ito等[17]研究了跌落高度和碰撞材料对马铃薯表面损伤的影响。Dwelle等[18]研究了温度、放置时间、气体压力对马铃薯内部黑斑瘀伤形成的影响。Thomson等[19]研究探讨了马铃薯块茎的特性及影响马铃薯表面损伤的处理技术。Nikara等[20]采用扫描电镜(Scanning electron microscope, SEM)分析了马铃薯碰撞后组织的微观力学变化。卢立新等[21]通过冲击力传感器研究苹果跌落的冲击力学特性，获得加速度与时间、冲击力与变形量、应力与应变以及跌落冲击弹性恢复系数等动态关系。李晓娟等[22]通过悬摆式碰撞试验，采用加速度传感器，得到苹果碰撞过程的加速度−时间曲线及数学表达式。Brusewitz等[23]使用力学传感器研究苹果的跌落，同时测得最大撞击力、撞击持续时间、到达最大撞击力所需要的时间、冲量等冲击特性的参数。刘治震[24]通过冲击力传感器对马铃薯进行跌落碰撞试验，得到马铃薯的临界损伤高度。卢琦[25]通过平膜盒测力传感器获得马铃薯的瞬时冲击力和冲击应力；通过与马铃薯的坚实度比较，得到马铃薯的损伤情况。胡奔[26]通过力学传感器获得马铃薯跌落碰撞的冲击力，经过转换研究得到马铃薯损伤的极限冲击应力。冯斌等[27]通过压力传感器测量马铃薯的冲击特性参数与损伤综合指数的关系。但是传统的测量方法不精确，普通的力学传感器只能测定马铃薯与其他材料发生碰撞时冲击力的大小，无法测定碰撞时的应力分布情况。对马铃薯与各作业部件的接触力学分布特性进行研究，可为马铃薯的机械化收割和自动化减损提供参考，具有重要的意义。
有学者用超声波技术[28]或TekScan5051柔性薄膜网格压力传感器[29]测量了苹果准静态下的压缩接触应力分布。美国和日本通过微胶囊颗粒缓释控制技术研制出了压力感应胶片，实现了动载下接触应力的精确测量，并成功应用于医学[30-31]。采用感压胶片对苹果挤压和碰撞的接触应力进行分析，并获得损伤面积和体积[32-33]。采用感压胶片对香梨跌落碰撞的接触应力分布进行研究，确定接触应力分布和香梨果实损伤间的关系[34-35]。
本文采用Prescale感压胶片及配套的压力图像数字解析系统对马铃薯跌落碰撞的接触应力分布进行研究，精确测量马铃薯在动载下的接触应力分布特性，研究马铃薯在动载下的力学特性和损伤规律，以实现对其机械损伤的可靠性评估和预测。
1. 材料与方法
1.1 材料
马铃薯试样为华中农业大学农业农村部马铃薯生物学与生物技术国家重点试验室试验田培育的‘中薯5号’，早熟品种，从出苗到收获约60~80 d，株高约50 cm，休眠期短。‘中薯5号’具有较好的增产性，薯块大而整齐、表面光滑、芽眼极浅、综合性状好，结薯集中，商品薯率高[36-37]。
收获期为2019年6月，在确保马铃薯无虫害、无损伤后迅速冷藏。马铃薯的含水率(w)为77.36%，贮藏温度为−5~5 ℃，相对湿度为80%~90%。
试验材料为土块、橡胶、塑料和钢板。土块采集自马铃薯收获时的试验田，含水率(w)为22.58%；橡胶为工业橡胶中的丁晴橡胶；塑料的材料是ABS。钢板采用65Mn钢，为我国薯类收获机上升运链建议采用的材料[38]。
高速摄像机（Pco.dimaxHD）为德国PCO公司生产。高速摄像机放置在距离试验台架水平2 m处，记录马铃薯跌落碰撞时的运动情况，选定的拍摄帧率为2 000帧/s，曝光时间为942 ns，像素分辨率为1 920 ×1 080 。
Fuji公司研发的感压胶片(Prescale®双层特超低压LLLW型)测量范围为0.2~0.6 MPa，测试精度≤±10%，配套的压力图像数字解析系统V2.0(FPD-8010E型)。
EPSON平板扫描仪(V370)为爱普生公司生产。
马铃薯跌落碰撞测试系统示意图见图1。
![]()
图 1 马铃薯跌落碰撞测试系统1：计算机；2：高速摄影机；3：跌落孔；4：调节阀；5：直尺；6：高度调节手柄；7：支架；8：碰撞材料；9：感压胶片；10：摄影范围Figure 1. Schematic diagram of potato drop crash test system and physical picture1：Computer；2：High-speed photographic instrument；3：Down hole；4：Regulating valve；5：Vertical ruler；6：Height adjusting handle；7：Stent；8：Collision material；9：Pressure sensitive film；10：Scope of photography1.2 方法
1.2.1 马铃薯跌落试验
马铃薯试样通过调节阀固定在跌落孔中，松开调节阀，马铃薯试样自由下落，保证其长轴扁平部位与不同的材料碰撞。反弹后在再次落地前接住马铃薯试样，保证马铃薯试样经过反弹达最大高度，也避免二次碰撞。高速摄影机记录全过程。
1.2.2 接触应力分布的测定
将感压胶片放置在碰撞接触材料表面，使感压胶片能够获取完整的碰撞接触部位。碰撞后的感压胶片，通过扫描仪进行平滑处理，获得接触应力分布的特征图，然后经相配套的压力图像数字解析系统V2.0读取特征图并进行数值解析，最后获得接触面积、碰撞冲击力、接触应力及其分布等参数。
1.2.3 马铃薯冲击压缩变形量的测定
通过马铃薯的冲击压缩变形量来衡量马铃薯的损伤。马铃薯的冲击压缩变形量越大，损伤越大。
冲击压缩变形量(D)表征马铃薯跌落碰撞时弹性形变的最大值，定义公式如下：
$$ D=\frac{\Delta s}{d} ,$$ (1) 式中，
$\Delta {{s}}$ 是指形变量，马铃薯从碰撞开始到压缩到最大值（下一秒会发生反弹）(图2)；d是与接触材料发生碰撞时，马铃薯竖直方向的直径。形变量
$\Delta s={\int }_{0}^{t}{v}_{t}{\rm d}t$ ，${v}_{t}={v}_{0}+{\int }_{0}^{T}a{\rm d}T$ ，T为碰撞弹性变形时间，通过高速摄影技术获得，t=2T。1.3 马铃薯试验方案设计
通过Design-expert设计了3因素3水平的组合正交试验方案，研究因素及组合正交试验设计方案的编码值和实际值的对应关系如表1所示。
表 1 因素水平编码表Table 1. Coding table of factor and level水平
Level碰撞材料(A)
Collision
material跌落高度(B)/mm
Dropping
height马铃薯质量(C)/g
Potato
mass−1 65Mn钢 65Mn steel 200 0~50 0 马铃薯 Potato 500 100~150 1 土块 Clods 800 150~200 跌落高度是根据马铃薯收获机振动筛与地面的高度进行筛选的。碰撞材料是根据马铃薯收获时可能发生碰撞的材料选取的。马铃薯质量是将收获后的马铃薯进行分级所得。最后通过分析较显著的因素，补充详细完整的因素水平，进行单因素试验。
2. 结果与分析
2.1 方差分析及响应曲面分析
组合正交试验的方差分析见表2。跌落高度和碰撞材料的P<0.000 1，对冲击压缩变形量的影响极显著；马铃薯质量的P= 0.026 5<0.05，对冲击压缩变形量的影响是显著的。而两两因素间的交互作用影响不显著，因素的二次方影响也不显著。通过方差分析获得冲击压缩变形量(D)与3个因素的二次回归模型。
表 2 冲击压缩变形量方差分析Table 2. Variance analysis of impact compression deformation amount方差来源1)
Source of variance离差平方和
Sum of deviation square自由度
Degree of freedom平均离差平方和
Mean sum of deviation squareF P 模型 Model 160.25 9 17.81 12.69 0.001 5 A 54.03 1 54.03 38.52 0.000 4 B 87.62 1 87.62 62.46 < 0.000 1 C 11.01 1 11.01 7.85 0.026 5 AB 0.83 1 0.83 0.59 0.467 4 AC 0.24 1 0.24 0.17 0.691 5 BC 0.19 1 0.19 0.13 0.725 8 A2 0.11 1 0.11 0.08 0.789 2 B2 6.30 1 6.30 4.49 0.071 8 C2 0.03 1 0.03 0.02 0.891 5 残差 Residual 9.82 7 1.40 失拟 Lack of fit 2.36 3 0.79 0.42 0.747 5 误差 Error 7.46 4 1.86 总计 Total 170.07 16 1)A为碰撞材料，B为跌落高度，C为马铃薯质量
1)A: Collision material; B: Dropping height; C: Potato mass$$ \begin{split} D= & 11.27-2.60A+3.31B-1.17C+0.45AB-0.25AC-\\ & 0.22BC+0.16{A}^{2}-1.22{B}^{2}+0.082{C}^{2}{\text{。}} \end{split} $$ (2) 失拟P=0.747 5>0.05，影响不显著，说明方差分析模型是可信的。
对跌落高度、碰撞材料和马铃薯质量的响应曲面分析如图3所示，二维的投影曲线反映因素间的关系，曲线在不同坐标轴具有不同疏密程度，疏密程度反映因素对冲击压缩变形量的不同贡献程度，曲线越密集说明对冲击压缩变形量的贡献越大。
![]()
图 3 冲击压缩变形量的响应曲面图Figure 3. Response surface diagram of impact compression deformation amount由图3a、3c可看出，随着跌落高度增加，曲线斜率均较陡峭，说明跌落高度对冲击压缩变形量的影响较显著，冲击压缩变形量随着跌落高度的增加而增大。由二维的投影可知，跌落高度的曲线密集程度比马铃薯质量的高，比碰撞材料的低。由图3a、3b可看出，不同的碰撞材料中马铃薯的冲击压缩变形量差异较大，碰撞材料的曲线密集程度比跌落高度、马铃薯质量的高。由图3b、3c可看出，当马铃薯质量较小时，曲线较陡峭，说明此时马铃薯质量对冲击压缩变形量影响较显著。马铃薯质量的曲线密集程度比跌落高度、碰撞材料都低。综上可得各因素对冲击压缩变形量贡献的程度为碰撞材料>跌落高度>马铃薯质量。因此后续着重分析碰撞材料和跌落高度对马铃薯跌落碰撞损伤的影响。
2.2 马铃薯跌落碰撞的损伤规律
为减少马铃薯外形差异对试验结果的影响，保证试验因素以外的环境和条件一致，选取的马铃薯为100~150 g，以减小其质量对试验的误差影响。不同跌落高度下马铃薯与不同材料碰撞的冲击压缩变形量见图4。
![]()
图 4 各接触材料的跌落高度与冲击压缩变形量的关系Figure 4. Relation between dropping height and impact compression deformation amount for each collision material由图4可知，随着跌落高度的增加，不同碰撞材料的冲击压缩变形量均呈上升趋势：土块和马铃薯的曲线交织上升，65Mn钢、塑料ABS和丁晴橡胶曲线上升趋势比较清晰。65Mn钢的冲击压缩变形量随着跌落高度增加而增加，在跌落高度>300 mm后，曲线呈近似直线上升趋势，冲击压缩变形量比其他材料的大。由于硬度比较大，65Mn钢碰撞时的弹性变形较小，使马铃薯的形变较大，因而65Mn钢的冲击压缩变形量比其他材料大。
塑料ABS的冲击压缩变形量曲线，首先是逐渐上升，在跌落高度>400 mm后趋于平缓。塑料ABS具有一定的弹性，对马铃薯的碰撞变形起到缓冲作用，因而塑料ABS的冲击压缩变形量处于中间水平，与土块和马铃薯的冲击压缩变形量相差不大。
马铃薯和土块的冲击压缩变形量曲线在跌落高度>400 mm后上升，在跌落高度>600 mm之后趋于平缓。马铃薯具有弹性缓冲性能，因此两两碰撞时，能够减少变形，而使冲击压缩变形量不大。土块本身具有一定的硬度，但由于收获时的土块含水量较高，能够对马铃薯的碰撞有一定的缓冲作用，加之土块具有较多的孔隙，跌落碰撞时的接触面积较大，从而减少了形变，因而土块的冲击压缩变形量处于中间水平。
丁晴橡胶的冲击压缩变形量在跌落高度为300 mm时最小，在跌落高度>400 mm后，曲线呈近似直线上升，在跌落高度为600~800 mm，冲击压缩变形量的差异很小。丁晴橡胶的冲击压缩变形量整体上比其他材料低。丁晴橡胶是缓冲材料，工业设计材料内部具有较多孔隙，能够对马铃薯起到缓冲减震的作用，使马铃薯的冲击变形较小，因而丁晴橡胶的冲击压缩变形量最小。
综上，不同碰撞材料造成冲击压缩变形量的差异与其接触应力、冲击力、接触面积等有关。
2.3 马铃薯碰撞接触应力的分布特性
2.3.1 不同碰撞材料的接触应力分布特性
马铃薯与5种接触材料进行碰撞试验，通过感压胶片及其压力数字解析系统获得马铃薯碰撞的接触应力及其分布。马铃薯与各种材料在不同高度碰撞时的接触应力与接触面积的关系见图5，不同接触材料在不同高度碰撞后的应力分布如图6所示。
![]()
图 5 不同跌落高度下马铃薯与不同材料碰撞的接触应力分布Figure 5. Contact stress distribution of potato in collision with different materials at different dropping height![]()
图 6 马铃薯与不同材料的碰撞接触应力典型分布图Figure 6. Typical distribution of contact stress on potato in collision with different materials马铃薯与各种材料在不同高度碰撞时，接触应力≤0.50 MPa占主要的接触面积。接触应力为0.50~0.60 MPa应是应力的峰值，因为接触应力>0.60 MPa的接触面积非常小，约为0.01~1.00 mm2，几乎可以忽略不计。接触应力在≤0.50 MPa的接触面积最大，对马铃薯的损伤起主要贡献作用；而接触应力为0.50~0.60 MPa的接触面积很小，不是造成马铃薯损伤的主要应力。
当接触应力≤0.20 MPa时，不同跌落高度间的接触面积相差不大。当接触应力>0.20 MPa后，跌落高度为200和300 mm的接触面积相差很小，而跌落高度≥400 mm的接触面积比200和300 mm的接触面积大很多。说明随着跌落高度增加，接触应力增大，接触面积也增大，即高应力占主要的接触面积，是导致马铃薯碰撞损伤的主要原因。
由图6可以看出，马铃薯与65Mn钢碰撞，当跌落高度为200 mm时，低应力区域(绿色区域)占主要面积；当跌落高度≥300 mm后，低应力区域分布在边缘，高应力区域占主要面积。马铃薯与塑料ABS、马铃薯、土块和丁晴橡胶碰撞，当跌落高度≤300 mm时，低应力区域占主要面积；当跌落高度≥400 mm后，低应力区域分布在边缘，而高应力区域占主要面积。从图6可看出，马铃薯与土块碰撞时的接触应力分布的轮廓非常分散，存在较多孔隙，同时接触面积也比其他材料大；马铃薯与65Mn钢、塑料ABS、马铃薯、丁晴橡胶碰撞的轮廓均呈较规则的椭圆形。
以上分析说明，与65Mn钢碰撞，跌落高度≥300 mm后，马铃薯开始出现损伤；与塑料ABS、土块、马铃薯、丁晴橡胶碰撞，跌落高度≥400 mm后，马铃薯开始产生损伤。
与65Mn钢碰撞，跌落高度<300 mm时，≤0.20 MPa的区域(低应力区域)占据主要的接触面积 (74.57%)；在跌落高度≥300 mm后，>0.20 MPa的区域(高应力区域)为主要接触面积 (74.37%)。与塑料ABS、马铃薯、土块、丁晴橡胶碰撞时，马铃薯在跌落高度<400 mm时，≤0.20 MPa的区域(低应力区域)占主要的接触面积，分别为200 mm的74.37%、60.38%、61.73%和65.68%，300 mm的58.23%、76.98%、71.42%和55.41%；在跌落高度≥400 mm时，>0.20 MPa的区域(高应力区域)占主要接触面积，分别为83.52%，73.31%，66.22%和83.87%。因此可以将马铃薯受到0.20 MPa的应力作为其损伤的临界应力。
2.3.2 马铃薯接触面积和接触应力
从图7可以看出，随着跌落高度的增加，不同材料的接触面积随之上升，且接触面积(Ac)随着跌落高度(H)上升呈高度线性正相关，决定系数(R2)均大于0.95。
![]()
图 7 马铃薯与不同材料碰撞时跌落高度与接触面积的关系Figure 7. Relation between dropping height and contact area of potato in collision with different materials不同材料的接触面积由高到低依次为：土块、65Mn钢、丁晴橡胶、马铃薯、塑料ABS。从图7可以看出，土块的接触面积十分分散，是因为土块本身具有较多孔隙，碰撞时通过增大接触面积，降低了马铃薯的变形和损伤。马铃薯与65Mn钢碰撞时的接触面积较大，是由于碰撞时65Mn钢塑性变形很小，形变集中于马铃薯果肉，使马铃薯碰撞变形和接触面积较大。马铃薯与马铃薯、丁晴橡胶碰撞时的接触面积较大，且马铃薯和丁晴橡胶均具有缓冲减震的作用，因此接触面积较大能够减少对马铃薯的损伤。马铃薯与塑料ABS碰撞时接触面积较小，应力集中，是马铃薯碰撞损伤严重的主要原因。
由此可见，马铃薯接触面积在不同碰撞材料上变化的原理不同，导致马铃薯损伤程度的不同。
如图8所示，马铃薯与各材料碰撞的接触应力随着高度增加呈正态分布。65Mn钢板和塑料ABS的接触应力较大，范围为0.20~0.31 MPa，随着跌落高度增大略微起伏；而丁晴橡胶、马铃薯土块的接触应力明显小于钢板和塑料ABS，为0.18~0.26 MPa。
![]()
图 8 马铃薯与不同材料碰撞的接触应力与跌落高度的关系Figure 8. Relation between contact stress and dropping height of potato in collision with different materials不同材料的接触应力具有显著差异，说明接触应力与马铃薯跌落的接触面积存在关联。
2.3.3 马铃薯碰撞接触应力分布与损伤的关系
马铃薯的接触应力（P）与接触面积（Ac）存在关联，而且接触应力与接触面积的乘积，即为碰撞冲击力（F），即：
$$ {{F}} = {{P}} \times {{{A}}_{\rm{c}}}{\text{。}} $$ (3) 马铃薯冲击压缩变形量与碰撞冲击力的线性拟合见图9，两者高度线性相关。冲击压缩变形量的回归方程见表3，除了土块(R2=0.963 4)，其他材料的R2均大于0.99。说明通过感压胶片测量的接触应力和接触面积能够准确地预测和评估马铃薯的损伤。
表 3 马铃薯与不同材料碰撞的冲击压缩变形量的回归方程Table 3. Regression equation of impact compression deformation amount of potato in collision with different materials碰撞材料
Collision material回归方程1)
Regression equationR2
Determination coefficient65Mn钢 65Mn Steel D=0.088 21AcP+5.370 0.993 1 塑料ABS Plastic ABS D=0.054 63AcP+6.932 0.992 9 马铃薯 Potato D=0.083 92AcP+5.308 0.992 1 土块 Clods D=0.034 42AcP+9.026 0.963 4 丁晴橡胶 Nitrile rubber D=0.027 18AcP+3.296 0.991 1 1) D表示冲击压缩变形量，Ac表示接触面积，P表示接触应力
1) D represents impact compression deformation amount, Ac represents contact area, P represents contact stress![]()
图 9 马铃薯与不同材料碰撞时冲击压缩变形量与冲击力的关系Figure 9. Relation between impact compression deformation amount and impact force of potato in collision with different materials3. 结论
本研究采用感压胶片及配套的压力图像数字解析系统分析马铃薯跌落碰撞的接触应力及其分布特性，获得马铃薯跌落碰撞的损伤预测模型。研究了碰撞材料和跌落高度对马铃薯碰撞损伤影响的数学关系，建立了相关模型。在方差分析中，马铃薯质量对其碰撞损伤同样具有显著影响，马铃薯质量对马铃薯碰撞损伤机理有待进一步的研究。
马铃薯组合正交试验的响应曲面分析表明，跌落高度、碰撞材料、马铃薯质量对马铃薯碰撞的冲击压缩变形量影响显著，且影响程度由大到小为碰撞材料>跌落高度>马铃薯质量。
跌落高度为200~800 mm时，接触应力≤0.50 MPa占主要的接触面积，对马铃薯的损伤起主要贡献作用。马铃薯接触应力0.50~0.60 MPa是接触应力峰值，因为马铃薯接触应力>0.60 MPa的接触面积约为0.01~1.00 mm2，可忽略不计。不同碰撞材料马铃薯损伤的跌落高度不一致，碰撞材料为65Mn钢时，在跌落高度≥300 mm后，马铃薯开始出现损伤，而塑料ABS、土块、马铃薯、丁晴橡胶在跌落高度≥400 mm后，马铃薯开始发生损伤。当碰撞材料的跌落高度<400 mm(65Mn钢<300 mm)时，≤0.20 MPa的区域占据主要的接触面积，当跌落高度≥400 mm(65Mn钢≥300 mm)时，(0.20，0.60] MPa的区域为主要接触面积。因此可以将0.20 MPa的接触应力作为马铃薯损伤的临界应力。跌落高度与接触面积呈高度线性正相关，R2均大于0.95。随着跌落高度增大，65Mn钢和塑料ABS的接触应力有略微起伏，为0.20~0.31 MPa；马铃薯、土块和丁晴橡胶的接触应力明显低于65Mn钢和塑料ABS，为0.18~0.26 MPa。冲击压缩变形量与接触应力和接触面积的乘积即碰撞冲击力呈高度线性相关，R2大于0.96，说明冲击压缩变形量与碰撞冲击力的回归方程模型能够准确地预测和评估马铃薯的损伤。
致谢：感谢宋波涛教授和陈惠兰教授对试验给予的帮助和鼓励！
-
![]()
图 1 大豆、甜玉米间作和单作种植模式
Figure 1. Layouts of soybean, sweet maize intercropping system and monocultures
![]()
图 3 不同根构型大豆与甜玉米间作对作物产量、生物量及磷吸收量的影响
+P：施磷处理，−P：不施磷处理；BD：‘本地2号’，YC：‘粤春03-3’；箱体上不同大写和小写字母分别表示施磷和不施磷条件下不同种植模式之间差异显著(P<0.05，Duncan’s法)
Figure 3. Effects of intercropping of soybean with different root architecture and sweet maize on crop yield, biomass and phosphorus uptake
+P: Treatment with phosphate fertilizer, −P: Treatment without phosphate fertilizer; BD: ‘Bendi 2’, YC: ‘Yuechun 03-3’; Different capital and lowercase letters on the boxes indicate significant differences among different planting modes under +P and −P conditions respectively (P<0.05, Duncan’s method)
![]()
图 4 不同根构型大豆与甜玉米间作的间作优势(A)和种间竞争力(B)
+P：施磷处理，−P：不施磷处理；BD：‘本地2号’，YC：‘粤春03-3’
Figure 4. Intercropping advantages (A) and interspecific competitiveness (B) of soybean with different root architecture intercropping with sweet maize
+P: Treatment with phosphate fertilizer, −P: Treatment without phosphate fertilizer; BD: ‘Bendi 2’, YC: ‘Yuechun 03-3 ’
![]()
图 5 不同根构型大豆与甜玉米间作对大豆根宽的影响
+P：施磷处理，−P：不施磷处理；BD：‘本地2号’，YC：‘粤春03-3’；箱体上不同大写和小写字母分别表示施磷和不施磷条件下不同种植模式之间差异显著(P<0.05，Duncan’s法)
Figure 5. Effects of intercropping soybean with different root architecture and sweet maize on root width of soybeans
+P: Treatment with phosphate fertilizer, −P: Treatment without phosphate fertilizer; BD: ‘Bendi 2’, YC: ‘Yuechun 03-3’; Different capital and lowercase letters on the boxes indicate significant differences among different planting modes under +P and −P conditions respectively (P<0.05, Duncan’s method)
![]()
图 6 不同根构型大豆与甜玉米间作对植株根系形态的影响
+P：施磷处理，−P：不施磷处理；BD：‘本地2号’，YC：‘粤春03-3’；箱体上不同大写和小写字母分别表示施磷和不施磷条件下不同种植模式之间差异显著(P<0.05，Duncan’s法)
Figure 6. Effects of intercropping soybean with different root architecture and sweet maize on root morphology of plants
+P: Treatment with phosphate fertilizer, −P: Treatment without phosphate fertilizer; BD: ‘Bendi 2’, YC: ‘Yuechun 03-3’; Different capital and lowercase letters on the boxes indicate significant differences among different planting modes under +P and −P conditions respectively (P<0.05, Duncan’s method)
![]()
图 7 不同根构型大豆与甜玉米间作对植株根际土壤有效磷含量的影响
+P：施磷处理，−P：不施磷处理；BD：‘本地2号’，YC：‘粤春03-3’；箱体上不同大写和小写字母分别表示施磷和不施磷条件下不同种植模式之间差异显著(P<0.05，Duncan’s法)
Figure 7. Effects of intercropping soybean with different root architecture and sweet maize on available phosphorus content in rhizospheric soil of plants
+P: Treatment with phosphate fertilizer, −P: Treatment without phosphate fertilizer; BD: ‘Bendi 2’, YC: ‘Yuechun 03-3’; Different capital and lowercase letters on the boxes indicate significant differences among different planting modes under +P and −P conditions respectively (P<0.05, Duncan’s method)
![]()
图 8 间作大豆和甜玉米表型参数间的相关性分析
1：大豆总生物量，2：大豆总磷吸收量，3：大豆总根长，4：大豆直径≤0.5 mm细根占比，5：大豆根际土壤有效磷含量，6：大豆根宽，7：甜玉米总生物量，8：甜玉米总磷吸收量，9：甜玉米总根长，10：甜玉米直径≤0.5 mm细根占比，11：甜玉米根际土壤有效磷含量；相关系数(R)在图中以不同颜色展示，右侧图例是不同R值的颜色区间，红色表示正相关，蓝色表示负相关；“*”表示在P<0.05水平显著相关(Spearman法)
Figure 8. Correlation analysis of phenotypic parameters in intercropped soybean and sweet maize
1: Total biomass of soybean, 2: Total phosphorus uptake of soybean, 3: Total root length of soybean, 4: Proportion of diameter ≤0.05 mm fine roots to total roots in soybean, 5: Available phosphorus content in soybean rhizospheric soil, 6: Root width of soybean, 7: Total biomass of sweet maize, 8: Total phosphorus uptake of sweet maize, 9: Total root length of sweet maize, 10: Proportion of diameter ≤0.05 mm fine roots to total roots in sweet maize, 11: Available phosphorus content in sweet maize rhizospheric soil; Colors in the heatmap indicate different correlation coefficients (R) and the color interval is shown in the right legend, red indicates positive correlation, blue indicates negative correlation; “*” indicates significant correlation at P<0.05 level (Spearman method)
-
[1] SHEN J, YUAN L, ZHANG J, et al. Phosphorus dynamics: From soil to plant[J]. Plant Physiology, 2011, 156(3): 997-1005. doi: 10.1104/pp.111.175232
[2] 刘瑾, 杨建军. 近三十年农田土壤磷分子形态的研究进展[J]. 土壤学报, 2021, 58(3): 558-567. [3] 赵学强, 潘贤章, 马海艺, 等. 中国酸性土壤利用的科学问题与策略[J]. 土壤学报, 2023, 60(5): 1248-1263. [4] REID J B, TAN Y, GRIFFIN W B, et al. Root system plasticity as a factor in crop breeding for improved nutrient uptake[C]//Agronomy New Zealand. 40th Agronomy Society Conference. Lincoln, 2010, 40: 51-64.
[5] KARLOVA R, BOER D, HAYES S, et al. Root plasticity under abiotic stress[J]. Plant Physiology, 2021, 187(3): 1057-1070. doi: 10.1093/plphys/kiab392
[6] HALING R E, BROWN L K, BENGOUGH A G, et al. Root hairs improve root penetration, root-soil contact, and phosphorus acquisition in soils of different strength[J]. Journal of Experimental Botany, 2013, 64(12): 3711-3721. doi: 10.1093/jxb/ert200
[7] SOUMYA P R, SHARMA S, MEENA M K, et al. Response of diverse bread wheat genotypes in terms of root architectural traits at seedling stage in response to low phosphorus stress[J]. Plant Physiology Reports, 2021, 26(1): 152-161. doi: 10.1007/s40502-020-00540-6
[8] HALING R E, YANG Z, SHADWELL N, et al. Root morphological traits that determine phosphorus-acquisition efficiency and critical external phosphorus requirement in pasture species[J]. Functional Plant Biology, 2016, 43(9): 815-826. doi: 10.1071/FP16037
[9] 刘鹏, 区伟贞, 王金祥, 等. 磷有效性与植物侧根的发生发育[J]. 植物生理学通讯, 2006(3): 395-400. [10] MIGUEL M A, POSTMA J A, LYNCH J P. Phene synergism between root hair length and basal root growth angle for phosphorus acquisition[J]. Plant Physiology, 2015, 167(4): 1430-1439. doi: 10.1104/pp.15.00145
[11] YI K, LI X, CHEN D, et al. Shallower root spatial distribution induced by phosphorus deficiency contributes to topsoil foraging and low phosphorus adaption in sugarcane (Saccharum officinarum L. )[J]. Frontiers in Plant Science, 2022, 12: 797635. doi: 10.3389/fpls.2021.797635
[12] XIA H, ZHAO J, SUN J, et al. Dynamics of root length and distribution and shoot biomass of maize as affected by intercropping with different companion crops and phosphorus application rates[J]. Field Crops Research, 2013, 150: 52-62. doi: 10.1016/j.fcr.2013.05.027
[13] LI L, LIU Y X, LI X F. Intercropping to maximize root-root interactions in agricultural plants: Agronomic aspects[M]//The root systems in sustainable agricultural intensification. 2021: 309-328. doi: 10.1002/9781119525417.ch12.
[14] 覃潇敏, 潘浩男, 肖靖秀, 等. 不同磷水平下玉米−大豆间作系统根系形态变化[J]. 应用生态学报, 2021, 32(9): 3223-3230. [15] AN R, YU R, XING Y, et al. Enhanced phosphorus-fertilizer-use efficiency and sustainable phosphorus management with intercropping[J]. Agronomy for Sustainable Development, 2023, 43(5): 57. doi: 10.1007/s13593-023-00916-6
[16] 雍太文, 刘小明, 宋春, 等. 种植方式对玉米−大豆套作体系中作物产量、养分吸收和种间竞争的影响[J]. 中国生态农业学报, 2015, 23(6): 659-667. [17] 唐劲驰, ISMAEL A M, 佘丽娜, 等. 大豆根构型在玉米/大豆间作系统中的营养作用[J]. 中国农业科学, 2005(6): 1196-1203. [18] 赵静, 付家兵, 廖红, 等. 大豆磷效率应用核心种质的根构型性状评价[J]. 科学通报, 2004, 49(13): 1249-1257. [19] 朱同林, 方素琴, 李志垣, 等. 基于图像重建的根系三维构型定量分析及其在大豆磷吸收研究中的应用[J]. 科学通报, 2006, 51(16): 1885-1893. [20] 郑林生, 刘志强, 陈康, 等. 基因型影响磷镁互作下大豆生长及根瘤和菌根性状[J]. 植物营养与肥料学报, 2022, 28(9): 1685-1698. [21] BURRIDGE J, JOCHUA C N, BUCKSCH A, et al. Legume shovelomics: High-throughput phenotyping of common bean (Phaseolus vulgaris L.) and cowpea (Vigna unguiculata subsp, unguiculata) root architecture in the field[J]. Field Crops Research, 2016, 192: 21-32.
[22] WALK T C, JARAMILLO R, LYNCH J P. Architectural tradeoffs between adventitious and basal roots for phosphorus acquisition[J]. Plant and Soil, 2006, 279(1/2): 347-366.
[23] AO J, FU J, TIAN J, et al. Genetic variability for root morph-architecture traits and root growth dynamics as related to phosphorus efficiency in soybean[J]. Functional Plant Biology, 2010, 37(4): 304-312. doi: 10.1071/FP09215
[24] 阮文亮, 彭松, 祝晓慧, 等. 减量施氮与间作大豆对甜玉米土壤团聚体及有机碳含量的影响[J]. 四川农业大学学报, 2023, 41(5): 811-819. [25] 白伟, 孙占祥, 郑家明, 等. 东北风沙半干旱区仁用杏作物间作对作物产量和土地生产力的影响[J]. 生态学杂志, 2017, 36(9): 2521-2528. [26] 董宛麟, 张立祯, 于洋, 等. 向日葵和马铃薯间作模式的生产力及水分利用[J]. 农业工程学报, 2012, 28(18): 127-133. [27] YANG W, LI Z, WANG J, et al. Crop yield, nitrogen acquisition and sugarcane quality as affected by interspecific competition and nitrogen application[J]. Field Crops Research, 2013, 146: 44-50. doi: 10.1016/j.fcr.2013.03.008
[28] BARRY K E, MOMMER L, VAN RUIJVEN J, et al. The future of complementarity: Disentangling causes from consequences[J]. Trends in Ecology & Evolution, 2019, 34(2): 167-180.
[29] LI C, HOFFLAND E, KUYPER T W, et al. Syndromes of production in intercropping impact yield gains[J]. Nature Plants, 2020, 6(6): 653-660 doi: 10.1038/s41477-020-0680-9
[30] YU X, XIAO S, YAN T, et al. Interspecific competition as affected by nitrogen application in sweet corn-soybean intercropping system[J]. Agronomy, 2023, 13(9): 2268. doi: 10.3390/agronomy13092268
[31] 刘天学, 张绍芬, 赵霞, 等. 我国玉米主要间作技术研究进展[J]. 河南农业科学, 2008(5): 14-17. [32] BETENCOURT E, DUPUTEL M, COLOMB B, et al. Intercropping promotes the ability of durum wheat and chickpea to increase rhizosphere phosphorus availability in a low P soil[J]. Soil Biology and Biochemistry, 2012, 46: 181-190. doi: 10.1016/j.soilbio.2011.11.015
[33] TANG X, PLACELLA S A, DAYDÉ F, et al. Phosphorus availability and microbial community in the rhizosphere of intercropped cereal and legume along a P-fertilizer gradient[J]. Plant and Soil, 2016, 407(1/2): 119-134.
[34] BEDOUSSAC L, JUSTES E. The efficiency of a durum wheat-winter pea intercrop to improve yield and wheat grain protein concentration depends on N availability during early growth[J]. Plant and Soil, 2010, 330(1/2): 19-35.
[35] ADAMS A E, BESOZZI E M, SHAHROKHI G, et al. A case for associational resistance: Apparent support for the stress gradient hypothesis varies with study system[J]. Ecology Letters, 2022, 25(1): 202-217. doi: 10.1111/ele.13917
[36] LI L, LI S, SUN J, et al. Diversity enhances agricultural productivity via rhizosphere phosphorus facilitation on phosphorus-deficient soils[J]. Proceedings of the National Academy of Sciences of the United States of America, 2007, 104(27): 11192-11196. doi: 10.1073/pnas.0704591104
[37] 安迪, 杨令, 王冠达, 等. 磷在土壤中的固定机制和磷肥的高效利用[J]. 化工进展, 2013, 32(8): 1967-1973. [38] HODGE A, BERTA G, DOUSSAN C, et al. Plant root growth, architecture and function[J]. Plant and Soil, 2009, 321(1/2): 153-187.
[39] 任媛媛, 李朝怡, 闫敏飞, 等. 黄土塬区玉米||大豆间作对根系形态的影响[J]. 草业科学, 2021, 38(12): 2449-2459. [40] 李隆. 间套作强化农田生态系统服务功能的研究进展与应用展望[J]. 中国生态农业学报, 2016, 24(4): 403-415. [41] LUO W, NI M, WANG Y, et al. Limited evidence of vertical fine-root segregation in a subtropical forest[J]. New Phytologist, 2021, 231(6): 2308-2318. doi: 10.1111/nph.17546
[42] ZHANG D, LI H, SHEN J. Strategies for root’s foraging and acquiring soil nutrient in high efficiency under intensive cropping systems[J]. Journal of Plant Nutrition and Fertilizers, 2017, 23(6): 1547-1555.
-
期刊类型引用(4)
1. 黎家蓉，劳秀娟，戴平凤，黄汝红，尤明双. 气象因子对于灵山县晚稻产量的影响分析. 农业灾害研究. 2025(01): 193-195 . 
百度学术
2. 黎家蓉，黄汝红，尤明双，戴平凤，劳秀娟. 气象因子对灵山县早稻产量的影响研究. 新农民. 2024(22): 81-83 . 百度学术
3. 赵平，柴铭，郝鑫炎，张亚青，耿纪允，孙晓杰. 衡水市典型作物生育期降水特征分析. 湖北农业科学. 2024(S1): 84-88 . 百度学术
4. 徐相明，顾品强，姚寅秋，王正大，汤晨阳，尹荔阳. 上海奉贤1960—2019年冷暖冬气候特征分析. 陕西气象. 2020(04): 30-33 . 百度学术
其他类型引用(2)
下载:
计量
- 文章访问数: 536
- HTML全文浏览量: 8
- PDF下载量: 28
- 被引次数: 6
