高志英,杨永辉,张兆沛
(1.山西运城农业职业技术学院 农林与工程系,山西运城 044000; 2.河南省农业科学院 植物营养与资源环境研究所,郑州 450002;3.农业部作物高效用水原阳科学观测实验站,河南原阳 453514;4.河南科技学院,河南新乡 453000)
磷作为主要营养元素之一,以磷化合物和磷酸离子参与许多重要物质组成和生命代谢活动,并对作物的生长发育和产量形成有重要影响[1-2]。大豆是主要的油料作物,缺磷会显著影响植株对其他养分的吸收,抑制大豆的生长,从而降低籽粒产量和品质[3-5],而过量施磷,易造成植株呼吸作用加强,养分大量消耗,不利于大豆干物质积累和产量的形成[2,6-7]。适量磷的供应可促进大豆对养分的吸收,碳水化合物的运输和产量的提高[8-9]。玉米-大豆间作是一种高效种植模式,在与玉米共同生长期内,低位作物大豆受高位玉米遮阴,减弱了大豆可接受的光照辐射[10-11],从而对大豆的形态及生理特征产生显著影响,如,茎秆变细,叶片变小变薄 ,光合速率降低[12-13],造成物质生产显著降低[14-16],玉米收获后,恢复自然光照,大豆进入单独生长期,在正常光照条件下,光合生产能力逐渐恢复[17]。关于遮阴[17]、磷肥单因素对大豆生长发育[18]、养分吸收[19]、物质积累分配[20-21]等方面已有大量研究,但在这两个因素共同作用下,干物质及磷素的积累动态变化尤其是积累过程各阶段特征如何,则鲜见报道。据此,本研究通过模拟玉米-大豆间作前期遮阴后期恢复自然光的特点,研究在不同磷用量下,苗期遮光复光后大豆干物质及磷素积累动态的变化特征,以期进一步阐明大豆干物质及磷素的积累规律,从而为生产实践中磷肥的合理施用提供科学依据。
试验在河南省新乡市河南科技学院园艺实训基地(35°18′ N,113°54′ E)进行,供试土壤为黄潮土,含有机质17.2 g·kg-1、碱解氮41.8 g· kg-1、速效磷 12. 3 mg· kg-1、速效钾110.6 mg ·kg-1。试验采用塑料盆栽,盆高30 cm,内径38 cm, 每盆装混匀干土10 kg,供试大豆品种为‘中黄61’。试验于2015-06-07播种,9月26日收获。播种前浇透水(以底部排水孔排出水为准),播种时,每盆4穴,每穴2粒大豆种子,出苗1周后间苗,每盆留长势均匀的4株苗株。试验采用两因素随机区组试验设计,二因素分别为磷肥用量和光照条件。磷肥用量(P2O5)设置0 mg·kg-1(P0)、75 mg·kg-1(P1)、125 mg·kg-1(P2)和175mg·kg-1(P3)4个水平;光照条件设自然光( L1)和遮光( L2)2个水平,自然光处理即全生育期自然光照,遮光处理即在一层遮阳网罩住的遮阳棚中进行,从播种开始进行遮阴处理,出苗后30 d结束遮阴,2个光照条件下4个施磷处理随机排列,每个处理10盆。
玉米各器官干物质量与磷素质量分数测定 分别于复光当天、复光后15 d(开花期)、30 d(结荚期)、45 d(鼓粒期)、75 d(成熟期)进行取样,各处理分别于不同盆随机取地上植株4株,样品于105 ℃杀青30 min,75 ℃烘至恒量,称量获得单株地上部干物质量(g),并在成熟期每处理取8株, 风干后获得单株产量(g)。地上部干物质称量后,干样品粉碎过筛,参照鲍士旦法[22],测定植株的磷素质量分数(mg/g),并依据公式:单株吸磷总量=植株干质量×植株磷素质量分数,算出植株磷积累量。
干物质及磷素积累特征各参数计算采用吴雨珊等[23]方法,用Logisitc方程对不同磷处理下大豆干物质与磷素积累的进程进行拟合,并计算相应的模型参数进行积累动态分析。其方程为:A=K/[1+e(a-bt)],其中,A为干物质(磷素)积累量;t为大豆复光后时间;K为大豆干物质(磷素)理论最大积累量;a、b为待定系数。对方程进行一阶和二阶推导,获得相关积累特征数,具体公式如表1。
表1 干物质与磷素积累过程特征数计算公式Table 1 Calculation formula of characteristic parameters on dry matter and phosphorus accumulation process
采用SPSS 10.0进行Logisitc方程拟合,利用Excel 2016和origin 8.0进行数据整理和作图。
由图1可知,随着生育期的推进,大豆干物质和植株磷素的积累量逐渐增加,两者均在复光取样前期(0~15 d)呈缓慢增加,随后快速增加(15~30 d),到生育后期(30~75 d)又呈缓慢增加。在各取样期,相同磷用量下,植株干物质及磷素积累量在2种光处理下,均呈现L1>L2,而在相同光照条件下,总干物质表现为P2最大,P0最小,增施磷肥促进了干物质及磷素的积累。在收获期,干物质及磷素积累量表现为P2显著高于P1、P0。
图中不同小写字母表示处理间差异显著(P<0.05),下同
Different lowercase letters mean significant difference among different treatment(P<0.05),The same below
图1大豆干物质和磷素积累动态变化
Fig.1Dynamicchangesofdrymatterandphosphorusaccumulationinsoybean
2.2.1 模型特征参数 大豆干物质积累量和磷素积累均可用Logistic 方程进行描述(表2),且各积累方程的决定系数R2都在0.97以上,F检验均达极显著水平,可见大豆干物质积累和磷素积累符合S型增长。不同光照条件下,随着施磷量的增加,干物质Rmax先增后降,当施磷量在P2处理,干物质Rmax最大,L1和L2分别为1.00 g·d-1和0.97 g·d-1。大豆的磷素Rmax在L1和L2光照环境下,分别在P3和P2处理下最大。干物质积累和磷素积累的Rmean在两光照条件下均与两者的Rmax变化趋势相同。此外,由表2可知,L1光照下,P0、P1、P2处理下干物质Tmax均早于磷素的Tmax,而L2光照下,干物质Tmax在各磷用量处理下也均早于磷素Tmax,这与干物质积累速率峰值比磷素峰值出现的早的结论相吻合(图2)。干物质与磷素积累速率最大时生长量Amax在L1、L2下均随施肥量增加先增后降,在P2处理达到最大。
2.2.2 积累速率 对Logistic积累方程进行求导,可得到各处理条件下干物质积累速率和磷素积累速率分别与时间关系的方程,以复光取样时间为横坐标,可获得不同处理条件下干物质及磷素积累速率分别与复光后时间的关系曲线图。由图2可知,大豆干物质积累速率与磷素积累速率在整个生育期均呈单峰曲线变化,但两者变化又有区别,表现为磷素的积累速率在L1、L2 光条件下对磷用量水平的响应较干物质敏感。另外可由图2可知,L1、L2下施磷可提高干物质积累和磷素积累的速率。
2.3.1 干物质积累阶段 由Logistic 曲线进一步将大豆干物质与磷素积累过程分为渐增期、快增期和缓增期(表3)。2光照条件下,干物质积累量在积累的3个阶段均随磷用量先升后降,并均在P2处理,即施磷量为125 mg·kg-1时达到最大,在L1、L2条件下,干物质积累量最大值在渐增期、快增期和缓增期分别为8.75 g、23.91 g、8.34 g和8.56 g、23.40 g、8.16 g。积累持续期在L1、L2条件下3个阶段分别为8.13~10.97 d、24.39~27.73 d、30.35~33.88 d和7.79~9.41 d、24.85~27.44 d、28.71~32.14 d,即缓增期>快增期>渐增期。而积累平均速率在L1、L2条件下分别为0.17~0.22 g·d-1、0.72~0.88 g·d-1、0.20~0.25 g·d-1和0.11~0.22 g·d-1、0.44~0.85 g·d-1、0.12~0.28 g·d-1,即快增期>缓增期>渐增期;2光环境下的3个积累阶段,积累平均速率均随施肥量先升后降,在施磷量为125 mg·kg-1时达到最大,说明适度的增加施肥量可以促进整个生育期积累速率增加,从而促进干物质的积累。
表2 苗期遮复光后不同磷肥处理大豆干物质与磷素积累特征参数Table 2 Characteristic parameters of dry matter accumulation and phosphorus accumulation in soybean underdifferent phosphate application rates after light recovery of seedling-shading
2.3.2 磷素积累阶段 由表3可知,L1光条件下,各阶段磷素积累持续期为缓增期>快增期>渐增期,分别为24.83~30.92 d、19.95~24.85 d、18.13~23.26 d;L2条件下,各阶段磷素积累持续期为缓增期>渐增期>快增期,分别为23.41~29.56 d、20.02~25.27 d、18.81~23.73 d。磷素积累量在L1、L2条件下,渐增期、快增期和缓增期分别为47.91~70.25 mg、130.91~191.94 mg、45.65~66.93 mg和40.80~70.83 mg、111.47~193.53 mg、38.87~67.48 mg,均表现为快增期>渐增期>缓增期,且在各阶段,磷素积累量随施磷量的增加呈先升后降趋势,均在P2处理条件,即施磷量为125 mg·kg-1下积累量最大。L1条件下,磷素积累平均速率随磷用量增加而增加,在P3处理,即磷用量为175 mg·kg-1下最大,而在L2条件下,磷素积累平均速率在P2处理达到最大,在整个生育期,平均速率表现为快增期>缓增期>渐增期。
2.4.1 大豆产量 由图3可看出,2种光照环境下的大豆籽粒产量均呈现P2>P3>P1>P0。L1环境下, P1、P2、P3产量分别较对照P0增加7.73%、12.16%、10.60%,而L2环境下,P1、P2、P3产量分别较对照P0增加7.80%、12.85%、11.44%。相同磷处理下,L2较L1降低幅度表现为P2 2.4.2 相关分析 由表4可知,大豆籽粒产量与干物质及磷素积累的最大速率(Rmax)、平均速率(Rmean)、最大速率时的物质量(Amax)均呈显著正相关,与干物质最大速率出现时间(Tmax)相关不显著,而与磷素的Tmax呈显著负相关,相关系数为-0.90。籽粒产量与渐增阶段的干物质及磷素积累量(A渐)、平均积累速率(R渐)、积累持续期(T渐)呈显著相关,并与A渐、R渐达极显著水平;与快增期、缓增期的干物质量(A快、A缓)、平均积累速率(R快、R缓)呈极显著相关,而与两阶段的持续期相关不显著;对与两阶段的磷素积累因子关系,与磷素A渐、A缓呈极显著相关,而与其他磷素积累因子相关不显著。 表3 苗期遮复光后不同磷肥处理大豆干物质与磷素积累 3 个阶段的特征参数Table 3 Characteristics parameters of dry matter and phosphorus accumulation under different phosphorusapplication rates after light recovery of seedling-shading in three stages of soybean 注:A、R和T分别代表3个积累阶段的干物质与磷素积累量、平均积累速率、持续时间。 Note:A、RandTrepresent the increased accumulation amount of dry matter and phosphorus, mean accumulation rate and accumulation duration at 3 stages. 图3 不同处理单株籽粒产量Fig.3 Seed yield per plant of soybeanunder different treatments 在玉米一大豆间作模式下,对大豆营养生长期进行遮阴处理会影响大豆形态特征及生物量的积累[12-14]。有研究对大豆不同生育时期遮荫发现,在出苗-盛花期间遮荫,主要影响主茎形态,而对产量影响不大,可能是在花期前复光,有利于叶片光合能力的恢复[24-25]。本研究发现,在大豆苗期进行遮光,同一施磷量的2种光照环境下,产量差异虽不显著,但增施磷肥缓解了因遮光可能引起对大豆的不利生长,提高最终的大豆产量,说明磷素对于大豆苗期遮光后进行复光具有一定的补偿效应。 表4 大豆产量与干物质、磷素积累因子的相关性Table 4 Correlation between yield and dry matter, phosphorus accumulation factors of soybeans 注Note:**P<0.01;*P<0.05. 大豆是需磷量较大的作物,磷肥的合理施用是调控大豆生长发育与产量形成的最重要措施之一[26]。乔振江等[27]研究发现,无论地上干物质量还是产量,磷高效基因型和磷低效基因型两个大豆品种均表现出高磷、中磷>低磷的趋势。本研究在4个不同施磷水平的试验表明,在各生育期,两种光照环境下的大豆干物质和产量表现为P2>P3>P1>P0,且作为干物质和产量形成基础的磷素积累量也表现相同趋势。因此在生产上,合理施用磷肥才能取得较好的经济效益。 作为描述作物生长特征的Logistic方程,在小麦[28]、大豆[29]、玉米[30]等作物上均有相关报道。此前, 有关大豆干物质生产特征研究主要侧重单因素,如品种[23],耕作方式[31]等,而对于遮复光、施磷量双因素对干物质及磷素积累特征的研究鲜见报道。相关分析表明,大豆籽粒产量与干物质及磷素积累的Rmax、Rmean、Amax、各阶段的积累量和T渐均呈显著相关。本研究发现,2种光照条件下,干物质及磷素积累的Rmax、Rmean及Amax均在P2处理下达最大。干物质及磷素的S曲线积累过程分为渐增期、快增期和缓增期[23],在磷肥用量0~125 mg·kg-1范围内,2光环境下各阶段的干物质及磷素积累量均随供磷量增施而增加。在渐增阶段,干物质积累的持续时间T渐在L1、L2光条件下,均在P2处理下最长。L1环境下, P1、P2、P3产量分别较对照P0增加7.73%、12.16%、10.60%,而L2环境下,P1、P2、P3产量分别较对照P0增加7.80%、12.85%、11.44%,说明增施磷对大豆产量的促进作用在L2环境下更大,且随着施磷量增加,促进作用增强。可见适当增施磷肥可补偿苗期因弱光胁迫对磷素与干物质积累造成的不利影响,从而促进了大豆产量的增加。 因此,本研究可以阐释在大田条件下,大豆与玉米间作时,增施磷肥有利于缓解因玉米对大豆的前期遮光而导致的减产,实现大豆稳产,甚至增产[32],从而为阐明间作条件下增施磷肥促进大豆增产的作用机制提供了有利依据。而增施磷肥对大豆不同生育期遮复光的补偿作用,及全生育期遮光的影响特征如何,尚需进一步研究。 Reference: [1] 孙 羲.植物营养与肥料[M].北京:中国农业出版社,1988:74-77. SUN X.Plant Nutrition and Fertilizer[M].Beijing:China Agriculture Press,1988:74-77. [2] 吴明才,肖昌珍,郑普英.大豆磷素营养研究[J].中国农业科学,1999,32(3):59-65. WU M C,XIAO CH ZH,ZHENG P Y. Study on the physiological function of phosphorus to soybean[J].ScientiaAgriculturalSinica,1999,32(3):59-65. [3] 吴冬婷,张晓雪,龚振平,等.磷素营养对大豆磷素吸收及产量的影响[J].植物营养与肥料学报,2012,18(3):670-677. WU D T,ZHANG X X,LONG ZH P,etal. Effects of phosphorus nutrition on P absorption and yields of soybean[J].PlantNutritionandFertilizerScience,2012,18(3):670-677. [4] ZARRIN F,HUHAMMAD Z,MFAYYAZ C. Effect of rhizobium strains and phosphorus on growth of soybean(Glycinemax) and survival of rhizobium and P solubilizing bacteria[J].PakistanJournalofBotany,2006,38(2):459 - 464. [5] 王建国,李兆林,李文斌,等.磷肥与大豆产量及品质的关系[J].农业系统科学与综合研究,2006,22(1):55-57. WANG J G,LI ZH L,LI W B,etal. Application of phosphorus in relation to soybean yield and quality [J].SystemSciencesandComprehensiveStudiesinAgriculture,2006,22(1):55-57. [6] 张小明,曾宪楠,孙 羽.磷素对大豆生长发育影响的研究进展[J].大豆科学,2016,35(1):176-180. ZHANG X M,ZENG X N,SUN Y. The progress of phosphorus fertilizer effect on soybean[J].SoybeanScience,2016,35(1):176-180. [7] 向达兵,郭 凯,杨文钰.不同磷钾处理下套作大豆干物质积累及钾肥利用率的动态变化[J].中国油料作物学报,2012,34(2):163-167. XIANG D B,GUO K,YANG W Y. Dynamics of dry matter accumulation and potassium utilization in relay strip intercropped soybean under different phosphorus and potassium levels[J].ChineseJournalofOilCropSciences,2012,34(2):163-167. [8] 刘克礼,高聚林,王立刚.大豆对氮、磷、钾的平衡吸收动态的研究[J].中国油料作物,2004,16(1):51-54. LIU K L,GAO J L,WANG L G.Study on dynamic balance assimilation of N,P and K in soybean[J].ChineseJournalofOilCropSciences,2004,16(1):51-54. [9] 曹景勤,陈碧云.黄淮海地区施肥技术对夏大豆增产效果的研究[J].土壤通报,1997,28(4):158-159. CAO J Q,CHEN B Y.Study on fertilization technology study of summer soybean in Huang-huai-hai region[J].ChineseJournalofSoilScience,1997,28(4):158-159. [10] 高 阳,段爱旺,刘祖贵,等.玉米和大豆条带间作模式下的光环境特性[J].应用生态学报,2008,19(6):1248-1254. GAO Y,DUAN A W,LIU Z G,etal. Light environment characteristics in maize soybeans tripinter cropping system[J].ChineseJournalofAppliedEcology,2008,19(6):1248-1254. [11] 范元芳,杨 峰,何知舟,等. 套作大豆形态、光合特征对玉米荫蔽及光照恢复的响应[J].中国生态农业学报,2016,24(5):608-617. FAN Y F,YANG F,HE ZH ZH,etal. Effects of shading and light recovery on soybean morphology and photosynthetic characteristics in soybean-maize intercropping system[J].ChineseJournalofEco-Agriculture,2016,24(5):608-617. [12] 任梦露,刘卫国,刘小明,等. 荫蔽信号对大豆幼苗生长和光合特性的影响[J].中国生态农业学报,2016,24(4):499-505. REN M L,LIU W G,LIU X M,etal.Effect of shading signal on growth and photosynthetic characteristics of soybean seedlings[J].ChineseJournalofEco-Agriculture,2016,24(4):499-505. [13] 李 瑞,文 涛,唐艳萍,等.遮阴对大豆幼苗光合和荧光特性的影响[J].草业学报,2014,23(6):198-206. LI R,WEN T,TANG Y P,etal. Effect of shading on photosynthetic and chlorophⅡ fluorescence characteristics of soybean[J].PrataculturaeSinica,2014,23(6):198-206. [14] GONG W Z,JING C D,WU Y S,etal. Tolerance vs. avoidance:two strategies of soybean(Glycinemax) seedlings in response to shade in intercropping[J].Photosynthetica,2015,53(2):259-268 . [15] SU B Y,SONG Y X,SONG C,etal. Growth and photosynthetic responses of soybean seedlings to maize shading in relay intercropping system in Southwest China[J].Photosynthetica,2014,52(3):332-340. [16] YANG F,HUANG S,GAO R C,etal.Growth of soybean seedlings in relay strip intercropping systems in relation to light quantity and red:far-red ratio[J].FieldCropsResearch,2014,155(155):245-253. [17] 吴雨珊,龚万灼,刘 鑫,等.苗期解除弱光胁迫后大豆恢复能力的生理评价[J].中国油料作物学报,2016,38 (6):804 -810. WU Y SH,GONG W ZH,LIU X,etal.Physiological evaluation index of recovery ability of soybean seedling after removing shade stress[J].ChineseJournalofOilCropSciences,2016,38(6):804 -810. [18] 丁玉川,陈明昌,程 滨,等.磷营养对不同大豆品种生长和磷吸收利用效率的影响[J].华北农学报,2006,21(1):121-124. DING Y CH,CHEN M CH,CHENG B,etal. Effect of phosphorus on plant growth and phosphorus uptake and use efficiency in different soybean cultivars[J].ActaAgricultureBoreali-Sinica,2006,21(1):121-124. [19] 乔振江,蔡昆争,骆世明.低磷和干旱胁迫对大豆植株干物质积累及磷效率的影响[J].生态学报,2011,31(19):5578-5587. QIAO ZH J,CAI K Z,LUO SH M.Interactive effects of low phosphorus and drought stress on dry matter accumulation and phosphorus efficiency of soybean plants[J].ActaEcologicaSinica,2011,31(19):5578-5587. [20] 唐 梅,李伏生,张富仓,等.不同磷钾条件下苗期适度水分亏缺对大豆生长及干物质积累的影响[J].干旱地区农业研究,2006,24(5):109-114. TANG M,LI F SH,ZHANG F C,etal. Effect of moderate water deficit during seedling stage on growth and dry matter accumulation of soybean under different P and K levels[J].AgriculturalResearchintheAridAreas,2006,24(5):109-114. [21] 刘 君,蔡昆争,骆世明,等.磷处理对不同磷效率基因型大豆根系性状和干物质积累及产量的影响[J].青岛农业大学学报(自然科学版),2008,25(3):211-214. LIU J,CAI K Z,LUO SH M,etal. Effects of phosphorous treatment on root traits and dry matter accumulation and yield of soybean genotypes with different P efficiency[J].JournalofQingdaoAgriculturalUniversity(NaturalScience),2008,25(3):211-214. [22] 鲍士旦.土壤农化分析[M].北京:中国农业出版社,2000:268-270. BAO SH D. Soil Agro-chemistrical Analysis[M].Beijing:China Agriculture Press,2000:268-270. [23] 吴雨珊,龚万灼,杨文钰,等.带状套作复光后不同大豆品种干物质积累模型与特征分析[J].中国生态农业学报,2017,25(4):572-580. WU Y SH,GONG W ZH,YANG W Y,etal. Dynamic model and characteristics analysis of dry matter production after light recovery of different soybean varieties in relay strip intercropping systems[J].ChineseJournalofEco-Agriculture,2017,25(4):572-580. [24] JIANG H F,EGLI D B.Soybean seed number and crop growth rate during flowering[J].AgronomyJournal,1995,87(2):264-267. [25] 王 一,杨文钰,张 霞,等. 不同生育时期遮阴对大豆形态性状和产量的影响[J].作物学报,2013,39(10):1871-1879. WANG Y,YANG W Y,ZHANG X,etal.Effects of shading at different growth stages on different traits and yield of soybean[J].ActaAgronomicaSinica,2013,39(10):1871-1879. [26] 王维军.大豆的磷素营养与施肥[J].中国农业科学,1963(11):41- 44. WANG W J. Phosphorus nutrition and fertilization of soybean[J].ScientiaAgriculturalSinica,1963(11):41-44. [27] 乔振江,蔡昆争,骆世明.低磷和干旱胁迫对大豆植株干物质积累及磷效率的影响[J].生态学报,2011,31(19):5578-5587. QIAO ZH J,CAI K ZH,LUO SH M. Interactive effects of low phosphorus and drought stress on dry matter accumulation and phosphorus efficiency of soybean plants[J].ActaEcologicaSinica,2011,31(19):5578-5587. [28] 李 娜,张保军,张正茂,等.不同施氮量和播量对‘普冰151’干物质积累特征及籽粒灌浆特性的影响[J].西北农业学报,2017,26(5):693-701. LI N,ZHANG B J,ZHANG ZH M,etal. Effects of different N application rate and seeding rate on dry matter accumulation and grain filling characteristics of ‘Pubing 51’ [J].ActaAgriculturaeBoreali-occidentalisSinica,2017,26(5):693-701. [29] 姜丽霞,那济海,朱海霞,等.基于 Logistic 方程的大豆产量预报方法[J].大豆科学,2008,27(3):414-419. JIANG L X,NA J H,ZHU H X,etal.Method for forecasting soybean yield based on logistic equation[J].SoybeanScience,2008,27(3):414-419. [30] 李文娟,何 萍,金继运.钾素营养对玉米生育后期干物质和养分积累与转运的影响[J].植物营养与肥料学报,2009,15(4):799-807. LI W J,HE P,JIN J Y. Potassium nutrition on dry matter and nutrients accumulation and translocation at reproductive stage of maize[J].PlantNutritionandFertilizerScience,2009,15(4):799-807. [31] 唐江华,苏丽丽,罗家祥,等. 不同耕作方式对夏大豆干物质积累及转运特性的影响 [J].核农学报,2015,29(10):2026-2032. TANG J H,SU L L,LUO J X,etal. Effects of different tillage methods on dry matter accumulation and transfer characteristics of summer soybean[J].JournalofNuclearAgrculturalSciences,2015,29(10):2026-2032. [32] 陈远学,周 涛,黄 蔚,等.麦 /玉米 /大豆间套作体系中小麦施磷后效对大豆产量、营养状况的影响[J].植物营养与肥料学报 2013,19(2):331-339. CHEN Y X,ZHOU T,HUANG W,etal.Phosphrous after effects on soybean yield and nutrition status in wheat / maize / soybean intercropping system[J].PlantNutritionandFertilizerScience,2013,19(2):331 -339.3 讨论与结论