王 祺,周荣飞,李宝鑫,张俊佩,张 强,裴 东,白永超*
(1. 林木遗传育种全国重点实验室,国家林业和草原局林木培育重点实验室,中国林业科学研究院林业研究所,北京 100091;2. 喀什地区林果产业工作站,新疆 喀什 844000;3. 新疆林业科学院,新疆 乌鲁木齐 830000)
核桃(Juglans regiaL. ) 为胡桃科(Juglandaceae)核桃属(JuglansL.)多年生落叶经济林木,是我国重要木本油料树种之一,被誉为“四大坚果”之首[1]。当前,我国大部分核桃产区规模扩张结束、面积趋于稳定,产业发展势头良好。然而,随着人们对核桃高营养价值的青睐以及我国粮油结构性短缺问题的严峻,在我国主栽区种植扩张面积趋于稳定的形势下[2],化学肥料因具有提高单位面积产量的功能而成为了农林业生产中最基础且最重要的物质投入[3]。基于此,种植户常以“大水大肥”的管理模式来提高产量,对土地、水资源进行长期的掠夺性经营[4],这样不仅忽视了土壤养分状况和树体自身的营养规律,亦没有考虑树体与土壤养分之间的平衡关系,致使施用的化肥仅有10%~40%可被植物吸收利用,残留的大部分化肥加剧了土壤盐渍化水平,植物所需营养元素严重失衡,对土壤环境和果实品质造成严重的负面影响[5-8],极大地限制了我国农林业的可持续健康发展。然而,有机肥作为补充或替代化肥的一种高有机质含量肥料,对土壤结构改善、土壤团粒形成、土壤肥力保持和有益土壤微生物丰度具有积极的促进作用。但是,由于有机肥的速效养分含量、养分释放率普遍低于化肥,短期内见效慢,加之有机肥施用时费时耗力,使得有机肥施用量大幅减少[9]。因此,本研究拟开展有机-无机肥配施,并探究其对核桃园土壤养分资源高效管理和树体健康的影响。
近年来,有机-无机肥配施被广泛用于探究其对改善土壤酶活性、根系分泌物、根际微生物群落组成结构及其功能特性的影响[10-11]。张倩[12]等研究表明,‘香玲’核桃园中化肥、有机肥和生物有机肥混合配施可显著改善土壤肥力,提高核桃产量及改善坚果品质;杜天宇[13]等研究结果表明,核桃园中化肥、有机肥和生物有机肥配施能增加土壤有机质、提高土壤肥力和核桃产量,改善坚果品质;Stamatiadis 等[14]研究认为,施用有机肥可有效调节土壤酸碱环境,改善植物根系微环境的水、肥、气、热及土壤结构等条件。矿源黄腐酸钾作为一种新兴的腐殖酸有机肥,富含大量腐殖酸、黄腐酸钾、有机质等有机组分,既能调控土壤酸碱度、提高土壤肥力[15],又能作为调节剂增强植物抗性、促进植物生长[16]。同时,矿源黄腐酸钾与大中微量元素混合施用,既可沟施也可冲施,可降低人工成本。张亚飞等[17]研究表明,黄腐酸钾与化肥控释可提高氮肥利用率,提高根系活力;庄振东[18]研究结果表明,腐植酸氮肥可减少土壤氮素淋溶损失,提高作物氮素利用率;此外,黄腐酸有机肥还可缓解盐碱胁迫下刺槐生理干旱[19],提高平邑甜茶和八棱海棠的耐盐性[20]。上述研究为核桃园合理施肥提供了参考。然而,新疆核桃主要分布在干旱和半干旱地区,土壤类型以砂土为主,常年气候干旱、降水稀少、蒸发量大,是典型的灌溉农业区与盐渍化严重发生区[21]。当前,新疆核桃主栽区以单施化肥、少施或不施有机肥的土肥管理措施为主,致使核桃园盐渍化不断加重,盐离子不断增加,土壤有机质严重匮乏,植物所需的营养元素严重失衡[22-23],植物逆境胁迫程度加剧,生长发育受阻,核桃焦叶症现象呈现逐年扩张且加重的趋势[24-25]。土壤理化指标与叶片矿质元素含量丰缺状况作为核桃园科学施肥的重要依据,针对南疆核桃主栽区有机-无机肥配施对核桃园的土壤化学特性和叶片矿质元素的影响研究鲜见报道。因此,明确有机-无机肥配施对土壤化学特性与叶片矿质元素的影响及其相关性对指导新疆核桃合理施肥和树体健康生长具有重要的理论意义。
本研究以新疆主栽品种10 年生‘温185’核桃为试材,在喀什地区叶城县核桃主栽区选取树势基本一致的核桃园进行有机-无机肥配施,分析有机-无机肥配施对土壤化学特性和叶片矿质元素的变化,筛选影响叶片养分的关键土壤化学因子,将为新疆核桃园土壤改良、合理施肥、叶片养分平衡等高效养分资源和树体健康管理提供理论基础和指导方案。
研究区位于新疆喀什地区叶城县洛克乡(37°53′39″ N,77°37′57″ E),海拔1 360 m,紧连塔克拉玛干大沙漠,属于暖温带大陆性干旱气候,年均气温在12.0 ℃以下,年均降雨量约50.0 mm,年均无霜期约235 d,土壤为砂质土壤,地势平坦,无防护林带干扰(图1)。
图1 研究区分布位置Fig. 1 Study area distribution location
供试植物材料为10 年生‘温185’核桃树,长势均一,树势中庸,树高(6.0 ± 0.5)m、地径(10.4 ± 0.7) cm、东西冠幅(3.9 ± 0.6)m、南北冠幅(4.4 ± 0.5)m,树形为小冠疏散分层形,东西向栽植,株行距5 × 7 m,灌溉方式以大水漫灌为主。在研究区以放射状施肥方式进行2 个施肥处理YZ 和YD,其中,YZ 是有机-无机肥处理,YD 措施是当地主要的化肥施用类型。YZ 为矿源黄腐酸钾(黄腐酸钾≥33%、腐殖酸≥35%、有机质≥60%、N≥10%、P2O5+ K2O≥8%)和中量元素水溶肥( Ca + Mg≥10%,EDTA-Zn≥0.1%,EDTA-Fe≥0.02%,B≥0.05%),YD 为磷酸氢二铵(N + P2O5≥64%)和硫酸钾(K2O≥52%,Cl-≤1.5%,S≥17%),每种施肥处理间设置防水畦,防止两个处理间发生灌溉水交换。每种处理设置5 亩试验地,3 次施肥位置示意图如图2 所示。施肥时通过4 条放射状沟施(长、宽、深分别为100 cm、 50 cm、30 cm),施肥时间见表1。
表1 3 次施肥时间及每棵树的施肥量Table 1 Time and amount of three rounds of fertilization
图2 3 次放射状施肥区域示意图Fig. 2 three diagrams of radial fertilization areas
分别在第1 次、第2 次、第3 次施肥前及核桃成熟采收前采集土壤和叶片样品。按照“S”型取样方式在试验区选取5 个点(每个点随机选取3 棵树,每棵树采集4 个方位混合为一个样本)进行根系集中分布层土壤样品采集(N=5),在采集土壤对应的树体上随机采集4 个方位的当年生营养枝上的成熟复叶(每个点随机选取3 棵树,每棵树采集4 个方位混合为一个样本,N=5),采集土壤样品时避开施肥区域。将采集的叶片一部分经干冰运输至实验室储存于-80 ℃冰箱用于叶片Cl--L 的测定,其余叶片经70 ℃ 恒温烘干、粉碎、过筛后,用于叶片矿质元素含量分析[26];将采集的土壤带回实验室后除去根系、砾石等杂质后,一部分保存在-20 ℃冰箱中用于分析铵态氮和硝态氮,另一部分经自然晾干、分级过筛后用于土壤化学特性分析。
土壤铵态氮和硝态氮采用紫外分光光度法测定;土壤有效磷、钾采用 Mehlich(M3)浸提剂浸提后,有效磷采用紫外分光光度法测定,有效钾采用原子吸收分光光度计测定,钠离子(Na+)采用NH4OAc-NH4OH 火焰光度法[27]。土壤有机质(OM)、pH、电导率(EC)、土壤氯离子(Cl--S)和叶片N、P、K、Ca、Mg、Fe、Mn、Zn、B、Na 元素的测定参照《土壤农化分析》并使用电感耦合等离子光谱仪测定[27-28],叶片氯离子Cl--L 含量测定采用离子色谱法。
采用 Excel 2019、 SPSS 19.0 和Origin 9.1 软件进行数据处理与统计学分析,采用单因素方差分析(p<0.05),Duncan 法进行多重比较。
应用典型相关分析方法,探究叶片矿质元素(氮(N)、磷(P)、钾(K)、钙(Ca)、镁(Mg)、钠(Na)、铁(Fe)、锰(Mn)、锌(Zn)、硼(B)、叶片氯离子(Cl--L))和土壤化学因子(pH 、 电导率(EC )、 有机质(OM)、土壤氯离子(Cl--S)、硝态氮(NO3--N)、铵态氮(NH4+-N)、有效磷(AP)、速效钾(AK)、钠离子(Na+))中各指标之间及叶片(氮(N)、磷(P)、钾(K)、钙(Ca)、镁(Mg )、 钠(Na )、 氯离子(Cl--L )、铁(Fe)、锰(Mn)、锌(Zn)、硼(B))与土壤之间的相关性;应用多元线性回归分析方法筛选影响叶片矿质元素的主要土壤化学因子,以土壤pH(x1)、EC(x2)、OM(x3)、Cl--L(x4)、NO3--N(x5)、NH4+-N (x6)、AP (x7)、AK(x8)、Na+(x9)为一个总体,核桃叶片矿质元素N(y1)、P(y2)、K(y3)、Ca(y4)、Mg(y5)、Na(y6)、Cl--L(y7)、Fe(y8)、Mn(y9)、Zn(y10)、B(y11)为另一个总体,建立多元线性回归方程,通过F检验,土壤化学因子中呈显著水平的指标被确定为影响核桃叶片矿质元素含量的重要土壤因子。同时对所建立的线性回归方程进行显著性检验,当F值呈显著(p<0.05)或极显著(p<0.01)水平说明建立的回归方程可靠。
2.1.1 有机-无机肥配施对核桃园土壤化学特性的影响 有机-无机肥配施对核桃园土壤化学特性影响的分析结果显示(图3),在5 月,YZ 处理中土壤NO3--N 和AK 含量显著高于YD 处理(分别高出460.22%、46.46%,p<0.05,下同),相反,土壤EC、Cl--S、NH4+-N、AP、Na+含量显著低于YD 处理(27.02%、22.86%、64.19%、53.82%、11.12%,p<0.05);在7 月,YZ 处理中土壤AK 含量显著高于YD 处理(42.35%,p<0.05),而土壤pH、EC、Cl--S、NH4+-N、AP、Na+含量则显著低于YD 处理(2.21%、27.25%、23.37%、61.09%、30.03%、13.24%,p<0.05);在9 月,YZ 处理中土壤OM、NO3--N、AK 含量显著高于YD 处理(22.17%、59.09%、56.47%,p<0.05),而土壤EC、NH4+-N、AP、Na+含量显著低于YD 处理(20.45%、37.31%、21.74%、11.36%,p<0.05)。综上,有机-无机肥配施(YZ)可显著降低核桃园土壤pH、Cl--S、Na+等因子含量,提高土壤NO3--N、OM、AK 等有效养分含量,由此说明YZ 处理可有效降低核桃园次生盐渍化水平。
图3 不同施肥处理下土壤化学特性的年变化Fig. 3 Annual changes of soil chemical factors under different fertilization treatments
2.1.2 土壤化学因子间的相关性分析 土壤化学特性是决定土壤肥力的重要组成部分,研究土壤化学因子间的相互关系,可为核桃园配方或精准施肥技术提供依据。核桃园土壤化学因子间的相关性分析结果如图4 所示,土壤pH 与OM,土壤EC 与AK,土壤Cl--S 与NO3--N、AK,土壤NO3--N 与NH4+-N,土壤NH4+-N 与AK 均呈显著或极显著负相关;土壤Cl--S 与NH4+-N、AP,土壤NO3--N与AK,土壤NH4+-N 与AP,土壤AP 与Na+均呈显著或极显著正相关,其他因子间的相关性较弱。由此可知,土壤化学因子的各指标之间可能存在协同或拮抗,在核桃园中开展配方或精准施肥时应考虑土壤化学因子间的相互作用。
图4 核桃园土壤化学因子间的相关性Fig. 4 Correlation coefficients between soil chemical factors in walnut orchards
叶片矿质元素含量丰缺状况可间接反映植株生长发育和果实品质的优良程度,有机-无机肥配施对核桃叶片矿质元素含量影响的分析结果如图5所示。从5 月—9 月,在YZ 处理下,叶片中N(2.90% ~3.32% )、K (1.76% ~2.29% )、Na(0.01% ~0.06% )、Fe (345.33 ~622.67 mg·kg-1)、Zn(7.72~12.70 mg·kg-1)元素含量均呈下降趋势,Cl--L(0.46%~0.88%)、Ca(2.45% ~3.28% )、Mn (152.00 ~280.67 mg·kg-1)元素含量均呈上升趋势,P(0.21%~0.28%)、B(98.60~187.67 mg·kg-1)元素含量均呈先上升后下降的趋势,7 月时达最高值;在YD处理下,叶片中N(2.75%~2.96%)、P(0.20%~0.25%)、K(1.24%~1.95%)、Na(0.01%~0.04%)、Fe(337.67 ~548.10 mg·kg-1)、Mn(97.43 ~128.33 mg·kg-1)、Zn(7.01 ~11.23 mg·kg-1)元素含量均呈下降趋势,Ca(2.12%~3.07%)元素含量呈上升趋势,Cl--L(0.55%~0.86%)、Mg(0.55%~0.75%)、B(112.33~142.33 mg·kg-1)元素含量呈波动趋势,其中,Cl--L(0.55%)、Mg(0.55%)元素含量在7 月时最低,而B(142.33 mg·kg-1)元素含量在7 月达到最高值。YZ 与YD 处理相比,在5 月,YZ 施肥处理下的叶片N、K、Zn 元素含量均显著高于YD(分别高出12.16%、17.46%、21.88%,p<0.05,下同),相反,Cl--L 含量则显著低于YD(28.20%,p<0.05);在7 月,YZ 施肥处理下的叶片N、P、Ca、Mg、Na、Fe、Mn、B 元素含量均显著高于YD(4.16%、13.25%、20.73%、44.04%、66.67%、24.62%、52.47%、31.85%,p<0.05),相反,Cl--L 含量则显著低于YD(15.94%,p<0.05);在9 月,YZ 施肥处理下的叶片K、Na、Mn 元素含量均显著高于YD(42.28%、75.00%、188.06%,p<0.05)。整体而言,YZ比YD 处理具有更丰富的叶片矿质元素含量(Cl--L 除外),由此推测YZ 处理下的叶片矿质元素含量水平更均衡。
图5 不同施肥处理下叶片矿质元素含量的年变化Fig. 5 Annual changes of leaf mineral element content under different fertilization treatments
2.3.1 核桃叶片矿质元素与土壤化学因子间的相关性分析 土壤养分是核桃树生长发育中必需矿质元素的主要来源,研究土壤化学因子与叶片矿质元素之间的相关性,可以更好地了解施肥对核桃叶片矿质元素的影响。土壤化学因子与叶片矿质元素之间的相关性分析结果如图6 所示,叶片N、K、Na、Fe、Zn 与土壤pH,叶片B 与土壤EC,叶片Cl--L 与土壤OM,叶片Ca、Mg、Mn 与土壤Cl--S, 叶片Ca 、Mn 与NH4+-N , 叶片Ca、Mg 与土壤AP,叶片Ca、Mg、Cl--L 与土壤Na+均呈显著或极显著负相关;反之,叶片Cl--L 与土壤pH,叶片K 与土壤EC,叶片N、K、Na、Fe、Zn 与土壤OM,叶片Mn 与土壤NO3--N,叶片B 与土壤AP,叶片Ca、Mn、B 与土壤AK,叶片P、Fe、Zn、B 与土壤Na+均呈显著或极显著正相关。然而,核桃叶片中的其他元素与土壤化学因子间的关系较为复杂,为进一步明确土壤化学因子与核桃叶片矿质元素间的相互关系,需借助多元统计分析方法研究其相关性。
图6 核桃果园土壤理化因子与叶片矿质元素相关性分析Fig. 6 Correlation coefficients between soil physicochemical factors and leaf mineral elements in walnut orchards
2.3.2 核桃叶片矿质元素与土壤化学因子间的多元线性回归分析 通过建立回归方程并对其进行显著性检验,方程均达到显著差异水平,说明建立的回归方程可靠。由表2 可知,土壤pH、Na+和Cl--S 是影响核桃叶片矿质元素含量的重要土壤化学因子。其中,土壤pH 显著影响叶片中N、P、K、Na、Fe、Zn 含量(p<0.05),降低土壤pH 有利于叶片N、P、K、Na、Fe、Zn 元素的积累;叶片Fe、Zn、B 含量受到土壤Na+的显著影响(p<0.05),并且随土壤Na+含量的增加而升高;降低土壤Cl--S 含量可能促进叶片Ca 元素的积累。
表2 影响核桃叶片矿质元素含量的土壤化学因子筛选及回归方程建立Table 2 Screening of soil chemical factors affecting mineral element content in walnut leaves and establishment of regression equation
与第二次全国土壤养分分级标准相比[29],试验区的土壤在施肥前pH、EC、Cl--S 和Na+含量普遍偏高(图3, 3 月),属于中度盐碱土,而土壤OM 和碱解N(NO3--N 和NH4+-N)含量则偏低,分别处于5 级和6 级,属于缺乏状态,由此可知,本研究区域的核桃园土壤肥力低,有效养分不足,盐碱化较严重,可能引起植物出现盐碱危害现象[30],在该区域进行合理的土壤养分资源管理对核桃健康生长发育非常必要。本研究结果表明,有机-无机肥配施可显著改善核桃园的土壤化学特性(图3),提高叶片矿质元素含量水平(图5)。根据前人的研究结果,施用有机肥可以显著提高土壤有机质含量[30-32],在本试验中,春季施用腐熟油渣后,YZ 和YD 处理(5 月)的有机质含量相比施肥前(3 月)分别提高了64.39% 和52.92%,说明腐熟油渣可作为有机肥替代农家肥施用提高土壤有机质含量,降低施肥成本,具有较好的应用潜力。然而,在植物生长季,YZ 和YD 处理对土壤化学特性的影响各异,YZ 处理可显著提高核桃园土壤N、P、K 等有效养分,而YD 处理提高了EC、pH、Na+、Cl--S、NH4+-N 等可能加剧土壤次生盐渍化致使核桃树发生盐害现象的盐碱因子的含量,两种施肥处理产生的差异可能是由于YZ 和YD 中的主成分和元素间的相互作用所致(图4)。对于YZ 处理显著提高土壤肥力的现象,推测可能是因为YZ 中的矿源黄腐酸钾含有丰富的腐殖酸、黄腐酸钾、有机质等有机组分,可改善土壤微生物群落结构及其功能特性和土壤酶活性[15,33],富集根际硝化细菌群落[34]或是降低土壤盐碱化程度而促进了硝化作用[35],从而提高了土壤中的NO3--N 含量,也可能是由于黄腐酸钾具有较强的螯合能力和腐熟油渣缓慢的氮素释放速率,改善了土壤团聚体和阳离子交换能力,增强了对N 素的固持作用,减少了NO3--N 的淋溶。此外,YZ 处理中的矿源黄腐酸钾中含有丰富的K2O,可能对盐碱土壤中的Na+毒害有一定的缓解作用[36]。更重要的是YZ 处理显著降低了Cl--S、Na+等土壤盐离子和土壤pH(5 月),这与前人研究含黄腐酸钾的肥料可减轻盐碱对植物根系胁迫的研究结果一致[37];在7 月,YZ 处理中土壤pH、EC、Cl--S、NH4+-N、AP、Na+因子含量分别比YD 处理低2.21%、27.25% 、23.37% 、61.09% 、30.03%、13.24%(图3,p<0.05),并且这些盐离子含量处于核桃树可耐受的最大范围内[38],说明有机-无机肥配施可显著降低土壤盐性离子含量,推测YZ 措施的应用将对新疆盐碱地核桃园土壤改善具有较好的促进作用,但在不同季节间表现出不同的效果,建议在7 月前进行2 次有机-无机肥配施。
土壤养分亏缺或过量对果树生长发育或产量形成均有负面影响,合适的土壤肥力水平是植物健康生长的基础,是经济林木稳产、高产与产品优质的重要前提。对于多年生经济林木而言,童期至初果期阶段,土壤肥力是树体长势、冠幅形成及枝干增粗拉长等营养生长期重要的果园影响因子,进入盛果期后,矿质元素的丰缺与平衡是影响产量形成和品质提高的重要因子。本研究结果表明,有机-无机肥配施显著改善了核桃叶片矿质元素含量水平,尤其YZ 处理,显著提高了5 月至7 月期间的叶片N、K、Zn 等元素含量,而叶片中的Cl--L 含量显著低于YD 处理,这可能是因为YZ 处理中土壤NO3--N、AK 含量显著高于YD 处理,其土壤养分易被核桃根系吸收利用。同时,YZ 处理中土壤EC、Cl--S 等盐离子显著低于YD 处理,减轻了盐性离子对其他矿质离子吸收的影响(图3),因此,YZ 处理中叶片矿质元素含量相对YD 处理较高。氯作为植物必需生长元素之一,约0.1%即可满足植物正常生长发育需要[39],前人[38]研究表明,在美国加州地区核桃叶片中,Cl--L 含量超过0.3%时,核桃树体会出现氯毒害现象,然而,本研究中YZ 和YD 处理下的叶片Cl--L 含量均高于0.3%,但均未出现氯毒害现象,这可能是因为对Cl--L 耐受强度的差异可能与品种、树龄、环境等因素有关[40]。但是,在新疆半干旱沙漠地区,盐碱化较严重,合理的水肥管理措施降低叶片Cl--L 含量对树体的正常生长发育极其重要[41],YZ 措施显著提高了土壤NO3--N 含量,可能抑制根系对Cl--S 的吸收[42]。进一步研究发现(表2),根系分布层土壤养分显著影响核桃叶片N、P、K、Fe、Zn 元素和Cl--L 含量,但叶片中多数元素与对应的土壤化学因子间无显著相关性,与前人对笃斯越橘(Vaccinium uliginosumL.)[43]、甜橙(Citrus sinensis(L.) Osbeck)[44]等研究结果较一致。多元线性回归分析结果表明(表2),叶片矿质元素含量主要受根系层土壤pH、Na+和Cl--S 的影响,尤其土壤pH 显著影响叶片中N、P、K、Na、Fe、Zn 含量(p<0.05),降低土壤pH 有利于叶片N、P、K、Na、Fe、Zn 元素的积累,这可能是因为在新疆半干旱荒漠化地区的耕地土壤以高盐碱为主[21],而盐离子(Na+和Cl--S)和高pH 是制约植物健康生长的主要限制因子[45]。本研究中有机-无机肥配施后,黄腐酸钾中的腐殖酸、有机质等有机组分对土壤pH 具有较强的缓冲能力和阳离子交换能力[46],除了可以显著改善土壤肥力外,还能减少氮肥施用和促进植物对其矿质养分的吸收利用,平衡树体营养水平[47],但需进一步研究明确增施矿源黄腐酸钾等高碳含量的有机肥对土壤微生物过程和养分循环的影响,为可持续生产系统下提高林果产量提供施肥参考。
有机-无机肥配施显著影响新疆核桃园土壤和叶片养分,YZ 处理显著提高了土壤养分含量,降低了盐性离子含量,同时平衡了叶片矿质元素含量,相反,YD 处理显著提高了土壤盐性离子含量,表明YZ 处理对新疆核桃园土壤改良和叶片矿质元素平衡具有积极的促进作用。此外,提高土壤NO3--N、OM、AK 含量,降低土壤pH、Na+、Cl--S 含量是核桃生长发育期的关键土肥管理技术,在今后的核桃园养分资源管理中,尤其在新疆核桃主栽区,应避免长期单一施用化肥,有机-无机肥配施将对解决南疆核桃主栽区次生盐渍化问题和缓解盐害发生提供有效保障。