Effects of straw mulching on nutrient release and nutrient balance in potato soil
-
摘要:目的
探索马铃薯种植过程中稻草覆盖还田模式下稻草N、P、K释放和土壤N、P、K养分平衡的特征。
方法采用盆栽试验的方法，对稻草覆盖还田后的腐解速率及养分释放状况进行研究。
结果稻草覆盖还田82 d后，其C、P含量变化较小，K含量降低93.96%，而N含量却提高了49.22%；覆盖稻草可为土壤提供大量养分且以K居多，每覆盖还田100 kg稻草可为马铃薯生长提供N 0.42 kg、P 0.14 kg(P2O5)和K 1.70 kg(K2O)；稻草腐解率为70.87%；各养分累积释放率表现为：K>C>P>N，C、N和P累积释放率分别为71.89%、56.85%和68.75%，K累积释放率最高，为98.22%。覆盖稻草和施用钾肥对土壤K平衡影响较大，不施钾肥和稻草土壤的K亏缺达163.80 kg/hm2，覆盖稻草或施用钾肥均可有效降低土壤钾素亏缺量，稻草与钾肥配合施用可使土壤K收支平衡，甚至盈余。
结论在稻草覆盖还田情况下，可适当减少钾肥用量，以降低生产成本，缓解我国钾肥资源匮乏的困境。
Abstract:ObjectiveTo study the characteristics of nitrogen(N), phosphorus(P) and potassium (K) nutrient release and soil N, P, K nutrient balance under straw mulching pattern during potato planting process.
MethodWe studied the decomposition rate and nutrient release status of mulching straw using pot experiment.
ResultStraw carbon (C) and P contents changed slightly after mulching for 82 days, while K content decreased by 93.96%, and N content increased by 49.22%. The mulching straw provided a large amount of nutrients (especially K) with 100 kg mulching straw contributing 0.42 kg N, 0.14 kg P (P2O5), and 1.70 kg K (K2O) to soil used for potato growth. Straw decomposition rate was 70.87%. The nutrient accumulative release rates were performed as: K (98.22%) > C (71.89%) > P (68.75%) > N (56.85%). Straw mulching and K fertilization largely affected soil K balance. The K deficient amount was 163.80 kg/hm 2 without straw mulching and K fertilizer application. K deficient amount was effectively reduced by straw mulching and K fertilization, and even K had surplus with combined application of straw and K fertilizer.
ConclusionThe straw mulching pattern can reduce K fertilizer application amount and production cost, and relieve the dilemma of K fertilizer resource deficiency in China.
-
Keywords:
- potato /
- nutrient release /
- nutrient balance /
- straw mulching /
- decomposition
-
丛枝菌根真菌(Arbuscular mycorrhizal fungi，AMF)是属于球囊霉门Glomeromycota的专性共生微生物，可以与陆地80%以上的植物根系共生。AMF自身生长所需的碳源完全依赖宿主植物，作为回馈，它们可以帮助宿主吸收矿质营养，提高宿主植物对生物及非生物胁迫的抗性[1]。研究发现AMF对植物种群、群落甚至是草地生态系统都有重要的调控作用[2-4]。
高寒草甸是青藏高原东南部主要的草地类型[5]，其主要建群植物为莎草科嵩草属的小嵩草Kobresia pygmaea、矮嵩草K. humilis、线叶嵩草K. capillifolia和藏嵩草K. tibetica等[6-7]，因此也称作高寒嵩草草甸。莎草科植物一直被认为是非丛枝菌根植物，但有研究发现在青藏高原AMF可以与嵩草属植物共生[8-9]，并对土壤团聚体的形成起重要作用[10]。基于AMF在草地生态系统中的重要性以及高寒嵩草草甸在青藏高原草地中所占的比例，研究高寒嵩草草甸生态系统中AMF的群落构建及动态对深入理解青藏高原草地土壤−植被相互作用具有重要意义。
近几十年来，受全球气候变暖和人类活动加剧的影响，青藏高原约1/3的草地发生了不同程度的退化[11]。施肥是退化草地恢复的有效方式之一，目前关于养分添加(主要是氮和磷)能否及如何影响高寒草甸AMF群落的研究相对稀少，并且仅限于海拔不超过3 500 m的区域[12-14]，不能代表海拔为3 200~5 200 m的青藏高原高寒嵩草草甸[7]。目前基于传统形态鉴别[15]和分子测序[16]的研究表明，随海拔提升，AMF侵染率、孢子数量和物种丰度下降，AMF群落组成也随之变化。不同海拔高度AMF群落对施肥的响应可能不同，为了更全面地认识高寒草甸，有必要选择海拔更高的区域进行研究。
本试验选取海拔约4 500 m的高寒嵩草草甸施肥样地，研究不同氮、磷添加对青藏高原高寒嵩草草甸AMF群落的影响及潜在驱动因子。本研究的核心科学问题为探究高寒嵩草草甸根系AMF群落对氮、磷添加的响应机制，及在添加氮、磷条件下主导AMF群落变化的因素。
1. 研究区概况和方法
1.1 研究区概况
研究样地位于中国科学院青藏高原研究所那曲生态环境观测站(E91°12ʹ~93°02ʹ，N30°31ʹ~31°55ʹ)，属于西藏那曲东南部，海拔4 480 m，年平均气温−2.1 ℃，年平均降水406.2 mm，高原山地气候，土壤类型为始成土，主要植物种类有小嵩草、矮嵩草和藏薹草Carex thibetica等。
1.2 样地设计及取样
2013年设置了48个5 m × 5 m的小区，完全随机设计，3个施氮梯度，4个施磷梯度，4次重复。施用的氮肥为尿素，各处理每年施用量分别为0、7.5、15.0 g·m−2，分别记为N0、N1和N2。施用的磷肥为过磷酸钙，有效成分为P2O5，各处理每年P2O5施用量分别为0、7.5、15.0、30.0 g·m−2，分别记为P0、P1、P2和P3。每年7月末施肥1次。
2016年8月18日，在每个小区随机挑选1个植被覆盖度相似的点作为样点。用土钻从每个样点均取深20 cm、直径8 cm的土块。将土块过2 mm筛，去除土和石块，将保留的根系装入自封袋作为1份样品，−20 ℃保存。另取部分过筛的土装入自封袋，4 ℃保存。共计取得48份根系样品和48份土壤样品，冷藏运送至北京林业大学草地资源与生态实验室。在进行后续试验分析前，将所有根系样品室温解冻，然后用清水洗净。
1.3 土壤化学成分分析
取适量土壤样品进行土壤化学成分分析。土壤有机碳用重铬酸钾滴定法[17]测定；土壤全氮用凯氏定氮法[18]测定；土壤铵态氮和硝态氮均用1 mol·L−1的氯化钾溶液浸提，然后用紫外分光光度法[19]测定；土壤全磷和有效磷分别用高氯酸−硫酸和碳酸氢钠浸提，并用钼−锑比色法[20-21]测定；pH在水土质量比为2.5︰1的条件下测量。
1.4 丛枝菌根真菌侵染率测定
根据Brundrett等[22-23]的方法，先从每份根系样品中剪取若干长1.5 cm的幼嫩根段，用台盼蓝染色。染色后在每份样品中挑出30个根段，切成1 cm根段，然后将每10个根段平行地放置在载玻片上，滴加乳酸甘油制片，每份样品制作3个玻片。在200×显微镜下根据McGonigle等[24]的方法测定AMF侵染率。
1.5 丛枝菌根真菌DNA提取、扩增及测序
在每份根系样品中，剪取40根长1 cm的幼嫩根段，用CTAB法[25]提取其中的DNA。提取出的DNA用无菌去离子水稀释20倍作为DNA原液，然后进行巢式PCR扩增。2轮PCR的反应体系均为25 μL：2 × Pfu PCR MasterMix (KP201) 12.5 μL，无菌去离子水8.5 μL，DNA模板2 μL以及2个引物各1 μL。第1轮PCR使用的引物为NS31(5ʹ-TTGGAGGGCAAGTCTGGTGCC-3ʹ)和AML2(5ʹ- GAACCCAAACACTTTGGTTTCC-3ʹ)，反应条件：94 ℃预热3 min；94 ℃变性45 s，53 ℃退火45 s，72 ℃延伸45 s，35个循环；72 ℃延伸10 min。第1轮PCR的产物稀释至1/10作为第2轮PCR的模板。第2轮PCR使用的引物为AMV4.5NF(5ʹ-AAGCTCGTAGTTGAATTTCG-3ʹ)和AMDGR(5ʹ-CCCAACTATCCCTATTAATCAT-3ʹ)[26]，其中AMDGR的5ʹ端带有由12个碱基组成的barcode。第2轮PCR反应条件与第1轮PCR相同，所得的产物用10 g·L−1 2 × TAE的琼脂糖凝胶进行电泳(电压120 V，时间30 min)检测。若样品有目标DNA条带(约280 bp)检出，取4 μL第2轮PCR模板用2倍的反应体系进行第2轮PCR。若样品无目标条带检出，用该样品的DNA原液重复巢式PCR，若2次重复均无目标条带，则去除该样品。根据此方法N1P1处理中的1个样品无条带，去除后剩余47个样品。切取含目标条带的凝胶，使用DNA胶纯化试剂盒(Axygen，美国)进行DNA纯化。纯化后的DNA使用Nanodrop 8000超微量分光光度计测定浓度，根据测定的浓度用移液枪从每个样品中吸取含100 ng DNA的溶液混合至1个离心管，送至中国科学院成都生物研究所在Illumina Miseq平台进行高通量测序。
1.6 生物信息学分析
测序得到的原始数据使用QIIME[27]进行分析。首先对数据进行质量控制，最大错误期望=0.5，最短序列长度=200，然后依次去除重复、chimeras和singletons。使用UPARSE-OTU算法和Silva数据库，依据97%的相似度阈值划分代表OTU并制作OTU表。在Excel中除去不属于球囊霉门、出现的样本数少于3以及序列数小于总序列量0.01%的OTU，并删除序列数小于总序列1%的样品以降低误差。对所有样品根据最小的样本序列量进行重新抽样，以减小样本量不同造成的差异。此外，将代表OTU序列上传至NCBI，与GenBank的已有序列进行比对，获得相近的参考序列。使用代表序列和参考序列在MEGA 7中使用p-distance模型，在bootstrap值为1 000下构建邻接树，直观查看代表OTU的系统发生情况。本研究获得的代表OTU序列均上传至欧洲核苷酸档案库，序列号为LR535994~LR536029。
1.7 统计分析
统计OTU数量作为OTU丰度，并计算AMF各科在各样品和各处理中的相对丰度。在R(3.2.2)中使用Vegan包中diversity功能计算每个样品的Shannon多样性指数。
进行方差分析前，将数据进行转换以符合数据分布的正态性和方差同质性。其中，对土壤化学成分、OTU丰度和Shannon多样性指数数据进行对数转换，侵染率和AMF科相对丰度数据作反正弦平方根转换。转换后的数据在JMP 11中进行双因素方差分析，分析氮、磷添加是否显著影响以上数据，具有显著差异(P<0.05)的结果使用Tukey's HSD检验进行多重比较。
在群落分析前，先应用R对OTU−样品矩阵进行Hellinger转换以减少稀有OTU的影响，然后用vegdist功能将矩阵转化为Bray-Curtis矩阵。在R中使用Vegan包中的adonis功能进行PerMANOVA分析，研究施氮、磷和其交互作用对AMF群落组成的影响。将土壤化学成分数据作为环境变量导入R并转化为Euclidean矩阵，与之前的群落矩阵一起进行Mantel分析，以研究土壤化学成分与AMF群落组成的关系。
最后，进行排序分析以研究不同处理AMF的群落分布情况。首先用原始的OTU−样品矩阵进行去趋势对应分析，根据排序轴长结果进一步使用典范对应分析，并用envfit功能进行Monte Carlo检验找出与AMF群落分布有显著相关的环境因子。最终使用ggplot包进行画图，直观显示各处理下AMF群落分布及其和环境因子的关系。
2. 结果与分析
2.1 土壤化学成分
方差分析得出，不同施氮处理铵态氮和硝态氮含量差异均极显著(P<0.001)；不同施磷处理全磷和有效磷含量也有极显著(P<0.001)差异；氮、磷添加的交互作用对硝态氮和有效磷含量有显著(P=0.021，P=0.040)影响(表1)。其中，N2处理铵态氮和硝态氮含量高于N0和N1处理；P3处理全磷和有效磷量高于其余施磷处理，P2处理有效磷含量高于P0处理(表2)。氮、磷添加及其交互作用未对土壤有机碳、全氮和pH产生显著影响(表1)。
表 1 不同施肥处理土壤化学成分含量差异显著性分析Table 1. Significance analyzes of soil chemical component content differences in different fertilization treatments指标 Index 氮添加
Nitrogen addition磷添加
Phosphorus addition氮、磷添加交互作用
Interaction between nitrogen and phosphorus additionsF P F P F P 有机碳 Organic carbon 0.610 0.549 0.940 0.431 1.849 0.117 全氮 Total nitrogen 0.737 0.486 0.343 0.794 0.674 0.672 铵态氮 Ammonium nitrogen 15.662 < 0.001 1.094 0.367 1.939 0.105 硝态氮 Nitrate nitrogen 21.780 < 0.001 2.711 0.061 2.964 0.021 全磷全磷 Total phosphorus 0.304 0.740 17.135 < 0.001 0.389 0.881 有效磷 Available phosphorus 0.502 0.610 33.642 < 0.001 2.509 0.040 pH 1.336 0.276 0.329 0.804 0.948 0.473 表 2 不同施肥处理试验样地的土壤化学成分含量1)Table 2. Soil chemical component contents of the experiment field in different fertilization treatments处理 Treatment w/% w/(mg·kg−1) pH N P 有机碳
Organic carbon全氮
Total nitrogen全磷
Total phosphorus铵态氮
Ammonium nitrogen硝态氮
Nitrate nitrogen有效磷
Available phosphorusN0 P0 4.79±0.45a 0.26±0.01a 0.040±0.001cd 2.53±0.77b 7.49±1.43bc 2.72±0.49d 6.74±0.08a P1 3.61±0.79a 0.21±0.06a 0.040±0.006cd 3.91±2.16b 5.23±1.25bc 6.39±2.34cd 6.77±0.11a P2 3.25±0.58a 0.23±0.04a 0.044±0.007abcd 2.06±0.27b 6.26±1.90bc 13.64±7.19bcd 6.66±0.06a P3 3.69±0.65a 0.27±0.09a 0.057±0.012abc 2.15±1.37b 4.67±1.14c 21.78±4.00abc 6.65±0.12a N1 P0 3.59±0.42a 0.22±0.01a 0.038±0.001d 1.83±0.73b 15.32±5.27bc 2.17±1.35d 6.69±0.08a P1 3.35±0.40a 0.24±0.01a 0.046±0.003abcd 2.25±1.07b 6.96±0.95bc 9.73±0.83bcd 6.47±0.15a P2 3.77±0.47a 0.22±0.04a 0.042±0.004bcd 6.07±4.38b 9.00±0.49bc 5.49±2.61d 6.73±0.19a P3 3.83±0.73a 0.24±0.04a 0.060±0.010ab 6.28±5.37b 5.79±0.91bc 26.02±13.21ab 6.67±0.23a N2 P0 2.98±0.32a 0.19±0.01a 0.038±0.003d 27.07±11.00a 35.66±15.70ab 2.28±0.12d 6.64±0.11a P1 3.57±0.94a 0.23±0.08a 0.044±0.004abcd 13.72±9.81ab 24.76±22.05bc 4.54±1.96d 6.61±0.15a P2 3.29±0.39a 0.24±0.05a 0.047±0.003abcd 7.13±0.85b 14.94±3.58bc 11.30±5.63bcd 6.61±0.05a P3 4.24±1.50a 0.23±0.03a 0.060±0.010a 15.41±11.90ab 59.55±23.70a 35.21±5.05a 6.60±0.07a 1)表中数据为平均值±标准差；同列数据后不同小写字母表示差异显著(P < 0.05，Tukey′s HSD法)；N0、N1、N2指氮施用量为0、7.5、15.0 g·m −2的处理；P0、P1、P2和P3指磷施用量为0、7.5、15.0、30.0 g·m−2的处理
1) Data in the table were mean value ± standard deviation; Different lowercase letters in the same column indicated significant differences (P<0.05, Tukey′s HSD test); N0, N1, N2 indicated nitrogen application amount 0, 7.5, 15.0 g·m−2; P0, P1, P2, P3 indicated phosphorus application amount 0, 7.5, 15.0, 30.0 g·m−22.2 丛枝菌根真菌侵染率
如图1所示，根系样品AMF侵染率平均值为73.1%。方差分析结果表明，氮、磷的施加及其交互作用对AMF侵染率无显著(P= 0.350，P = 0.119，P = 0.562)影响。
![]()
图 1 不同施肥处理下丛枝菌根真菌侵染率N0、N1、N2指氮施用量为0、7.5、15.0 g·m−2的处理；P0、P1、P2和P3指磷施用量为0、7.5、15.0、30.0 g·m−2的处理Figure 1. Arbuscular mycorrhizal fungal colonization rates in different fertilization treatmentsN0, N1, N2 indicated nitrogen application amount 0, 7.5, 15.0 g·m−2; P0, P1, P2, P3 indicated phosphorus application amount 0, 7.5, 15.0, 30.0 g·m−22.3 测序总体结果
测序数据除杂后，剩余43个样品，共获得234 041条AMF序列。每个样品序列数从297至15 397不等。将所有样品重新抽样至样本量为297后，获得36个OTU，每个样品的OTU数量从1至12个不等，平均值为5.6个。根据稀释曲线(图2)可以看出，重新抽样后N0P2、N1P0、N1P3、N2P2和N2P1对应的曲线末端趋于平滑，说明测序深度较为饱和。经比对，所有AMF的OTU分属7个科，其中球囊霉科有27个OTU，10 285条序列，占总序列的80.5%；多孢囊霉科Diversisporaceae有3个OTU，666条序列，占总序列的5.2%；碎球囊霉科Claroideoglomeraceae有2个OTU，146条序列，占总序列的1.1%；无梗囊霉科Acaulosporaceae、原囊霉科Archaeosporaceae、巨孢囊霉科Gigasporaceae和和平囊霉科Pacisporaceae各有1个OTU，序列数分别为441、532、351和350条，共占总序列的13.1%(图3)。根据OTU代表序列和参考序列构建的邻接树如图4所示。
![]()
图 2 不同施肥处理下分子测序的稀释曲线N0、N1、N2指氮施用量为0、7.5、15.0 g·m−2的处理；P0、P1、P2和P3指磷施用量为0、7.5、15.0、30.0 g·m−2的处理Figure 2. Rarefaction curve of sequencing samples in different fertilization treatmentsN0, N1, N2 indicated nitrogen application amount 0, 7.5, 15.0 g·m−2; P0, P1, P2, P3 indicated phosphorus application amount 0, 7.5, 15.0, 30.0 g·m−2![]()
图 3 不同施肥处理下丛枝菌根真菌各科的相对丰度N0、N1、N2指氮施用量为0、7.5、15.0 g·m−2的处理；P0、P1、P2和P3指磷施用量为0、7.5、15.0、30.0 g·m−2的处理Figure 3. Relative abundance of different arbuscular mycorrhizal fungal families in different fertilization treatmentsN0, N1, N2 indicated nitrogen application amount 0, 7.5, 15.0 g·m−2; P0, P1, P2, P3 indicated phosphorus application amount 0, 7.5, 15.0, 30.0 g·m−2![]()
图 4 36个丛枝菌根真菌OTU代表序列及其参考序列构建的邻接树模型：p-distance；Boostrap值：1 000；DQ846895作为outgroupFigure 4. Neighbor-joining tree constructed based on representative sequences of 36 arbuscular mycorrhizal fungal OTUs and their reference sequencesModel: p-distance;Boostrap value: 1 000; DQ846895 was used as an outgroup2.4 氮、磷添加对丛枝菌根真菌群落的影响
方差分析结果显示，氮、磷添加及其交互作用对OTU丰度和Shannon多样性指数无显著影响(图5A、5B)。但施氮对球囊霉科的相对丰度有显著影响(P<0.001)，N2处理球囊霉科的相对丰度显著低于N1(图5C)。
![]()
图 5 不同施肥处理丛枝菌根真菌OTU丰度、Shannon多样性指数及球囊霉科相对丰度不同柱子上不同大写字母表示在P<0.001水平差异显著(Tukey’ s HSD检验)Figure 5. OTU richness, Shannon diversity index of arbuscular mycorrhizal fungi and relative abundance of Glomeraceae in different fertilization treatmentsDifferent capital letters on different columns indicated significant differences at P<0.001 level (Tukey’ s HSD test)PerMANOVA结果表明氮、磷添加处理及其交互作用对AMF群落组成无显著影响(P = 0.680，P= 0.473，P= 0.589)。Mantel分析结果显示AMF群落组成与有机碳含量和硝态氮含量有显著正相关(r=0.176, P=0.110；r=0.142, P=0.041)关系。此外，Monte Carlo结果表明有机碳、硝态氮、全磷和有效磷含量均对AMF群落结构有显著影响(r=0.04, P=0.001；r=0.327, P=0.013；r=0.185, P=0.025；r=0.188, P=0.020)。这4个环境因子在典范对应分析中对排序结果的解释量为11.70%，第1轴和第2轴的解释量分别为6.13%和2.61%。(图6)。
![]()
图 6 不同施肥处理丛枝菌根真菌群落以及显著环境变量的典范对应分析图Figure 6. Canonical correspondence analysis plot of arbuscular mycorrhizal fungal community distribution and significant environmental variables among different fertilization treatments3. 讨论与结论
在农田和天然草地生态系统中，施加氮和磷会降低AMF侵染率、OTU丰度和多样性，改变根系和根际土壤中AMF的群落结构[28-29]。研究认为出现这种现象的原因是土壤养分的增加使植物更多依靠自身根系去吸收营养，降低对帮助其吸收养分的微生物的依赖[30]。本研究发现，施氮显著提高了土壤中铵态氮和硝态氮的含量，施磷显著增加了土壤中全磷和有效磷的含量，但氮、磷的施加对AMF侵染率、OTU丰度、多样性以及群落组成均无显著影响。因此，本研究的结果与此前在其他生态系统下的研究结果存在差异[28-29]，该差异可能是由以下原因造成的。
首先，差异可能是海拔不同造成的。目前相关研究[12-14]的海拔不超过3 500 m，低于本试验样地的海拔4 480 m。海拔会显著影响AMF的群落和功能，Gai等[15]和Shi等[31]通过形态学方法发现不同海拔条件下AMF侵染率、孢子密度以及菌丝密度都会有显著变化；Liu等[16]和Li等[32]通过分子测序手段发现不同海拔条件下AMF群落组成有显著差异。本研究与此前研究的海拔不同，气候、植被及土壤条件都会有不同，进而导致AMF对氮、磷添加的差异性响应。
第二，在高寒地区AMF不仅能帮助植物进行养分吸收，同时也可帮助植物应对环境胁迫。因此即使在施肥条件下植物不需要AMF来帮助其获取养分，还是要与AMF保持共生关系应对胁迫。Xiang等[13]在海拔3 220 m的青藏高原高寒草地施肥3年，发现同时添加氮和磷时AMF的OTU丰度和多样性均显著高于对照，也高于氮、磷单独添加时的水平。Chen等[33]研究认为，虽然在养分充足的条件下，植物在营养吸收上降低了对AMF的依赖，然而，植物依然需要AMF来帮助其对抗环境胁迫，尤其是球囊霉科可以提高宿主抗寒性。在本研究的高寒草甸生态系统中，一方面氮、磷的添加使得植物群落降低了对AMF吸收养分的依赖，另一方面植物群落却增加了依靠AMF来对抗胁迫的需要。因此，2个效应综合使得在氮、磷添加的情况下AMF的OTU丰度、多样性和群落组成未受显著影响。
另一个可能引起差异的原因是样品内及样品间的差异较大，影响了氮、磷添加的效应。从稀释曲线可以看出，个别样品的曲线在终点处仍在上升，说明其测序深度未达到饱和，此时增加取样量可能会有更多的OTU出现。同时，通过典范对应分析图可以直观地发现，N2P3处理3次重复的AMF群落分布距离较大，该处理3次重复的AMF群落的差异是所有处理的重复间差异最大的。基于样本内的误差，本试验结果存在一定限制性。
本研究中，PerMANOVA分析发现，施加氮、磷不影响AMF群落组成，但通过Mantel分析发现AMF群落的组成与有机碳和硝态氮含量有显著的正相关关系。典范对应分析的结果表明有机碳、硝态氮、全磷和有效磷含量与AMF群落的分布有显著的正相关关系，其中有机碳与AMF群落分布的相关性最大。有机碳主要对低磷条件(P0和P1处理)AMF群落有影响，硝态氮、全磷和有机碳主要对高磷条件(P3处理)AMF群落有影响。相关研究也表明，土壤有机质、硝态氮和有效磷均是对高寒草甸AMF群落有显著影响的土壤成分[34-35]。
Zheng等[34]发现施氮能改变AMF群落，但不是通过直接作用，而是间接地通过改变土壤其他成分(如有机碳)和植物群落来影响AMF群落。土壤成分也可以直接或间接地影响微生物群落[36-37]。因此可以解释本研究中施肥未对AMF群落产生影响，但特定的土壤因子与AMF群落分布有显著相关性。
综上所述，基于青藏高原4 500 m海拔高寒草甸的研究发现，氮、磷添加对AMF的侵染率、OTU丰度和多样性无显著影响。该发现回答了本研究要明确的第1个核心科学问题，即青藏高原高寒嵩草草甸根系中的AMF群落不受氮、磷添加的影响，与海拔较低的AMF群落对氮、磷的响应不同；对第2个核心科学问题，影响青藏高原高寒嵩草草甸根系AMF群落的因子主要是有机碳和硝态氮含量，全磷和有效磷含量是次要因子。未来研究还需要系统选取不同海拔的样地，以获取青藏高原不同海拔高寒草甸AMF群落对施肥的响应。另外，在全面剖析AMF群落变化规律及驱动机制的基础上，亟需研究AMF群落变化如何反馈影响地上植被的个体生长与群落动态，探索高寒草地退化及恢复演替中的地下生态过程及机理。
-
![]()
图 1 稻草钾肥配施对稻草C、N、P和K含量的影响
各图中，相同指标柱子上方的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
Figure 1. Influences of combined application of straw and K fertilizer on C, N, P, K contents of rice straw
In each figure, different lowercase letters on bars of the same indicator indicated significant difference(P<0.05, Duncan’ s test)
![]()
图 2 稻草钾肥配施对稻草腐解率及C、N、P和K累积释放率的影响
相同指标柱子上方的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
Figure 2. Influences of combined application of paddy-rice straw and K fertilizer on paddy-rice decomposition rates and accumulative release rates of C, N, P, K
Different lowercase letters on bars of the same indicator indicated significant difference(P<0.05, Duncan’ s test)
表 1 盆栽试验各处理的施肥和覆盖稻草量
Table 1 Amounts of fertilizer application and straw mulching in pot experiments
kg·hm−2 处理
TreatmentN P2O5 K2O 稻草
Paddy-rice strawK0S0(CK) 160.7 36.9 0 0 K0S2 160.7 36.9 0 6 000 K1S2 160.7 36.9 99.8 6 000 K2S2 160.7 36.9 199.7 6 000 K3S2 160.7 36.9 299.5 6 000 K2S0 160.7 36.9 199.7 0 K2S1 160.7 36.9 199.7 3 000 K2S3 160.7 36.9 199.7 9 000 K1S1 160.7 36.9 99.8 3 000 
下载: 导出CSV
表 2 马铃薯收获时稻草养分累积释放量1)
Table 2 Accumulative release amounts of nutrients at harvesting period of potato
kg·hm−2 处理 Treatment C N P K K0S2 1 528.69±19.86b 25.76±0.14b 3.55±0.02b 102.32±0.01b K1S2 1 508.22±140.85b 26.60±3.96b 3.65±0.31b 102.35±0.39b K2S2 1 573.74±109.86b 27.39±2.30b 3.58±0.19b 102.35±0.23b K3S2 1 456.05±75.99b 24.39±1.96b 3.15±0.28b 101.47±0.48b K2S1 764.81±19.93c 12.23±0.69c 1.70±0.07c 51.07±0.12c K2S3 2 221.43±135.76a 36.58±2.90a 4.94±0.28a 152.74±0.46a K1S1 723.99±17.39c 11.13±0.64c 1.71±0.03c 51.03±0.02c 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
表 3 稻草钾肥配施对土壤N平衡的影响1)
Table 3 Influences of combined application of paddy-rice straw and K fertilizer on soil N balance
处理
Treatment投入量/(kg·hm−2) Input amount 支出量/
(kg·hm−2)
Output amount表观盈亏量/(kg·hm−2)
Apparent surplus and deficient amount表观平衡系数
Apparent balance coefficient实际平衡率/%
Actual equilibrium rate化学N
Chemical N稻草N
Straw NK0S0 160.65 0 89.20a 71.45±1.34e 1.80 80.19 K0S2 160.65 43.35 74.43b 129.57±3.62b 2.76 175.42 K1S2 160.65 43.35 88.63a 115.37±7.56bc 2.33 133.34 K2S2 160.65 43.35 73.63ab 130.37±4.12b 2.78 178.77 K3S2 160.65 43.35 77.41ab 126.59±9.37b 2.72 171.59 K2S0 160.65 0 83.97ab 76.68±0.29e 1.91 91.33 K2S1 160.65 21.68 82.01ab 100.31±1.54d 2.23 122.47 K2S3 160.65 65.03 75.40ab 150.28±3.78a 3.01 200.91 K1S1 160.65 21.68 80.42ab 101.90±2.54cd 2.27 127.15 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
表 4 稻草钾肥配施对土壤P平衡的影响1)
Table 4 Influences of combined application of paddy-rice straw and K fertilizer on soil phosphorus balance
处理
Treatment投入量/(kg·hm−2) Input amount 支出量/
(kg·hm−2)
Output amount表观盈亏量/(kg·hm−2)
Apparent surplus and deficient amount表观平衡系数
Apparent balance coefficient实际平衡率/%
Actual equilibrium rate化学P
Chemical P稻草P
Straw PK0S0 16.11 0 13.53a 2.58±0.07d 1.19 19.09 K0S2 16.11 5.0 12.77a 8.33±1.33ab 1.69 50.60 K1S2 16.11 5.0 15.01a 6.1±1.18bc 1.42 42.32 K2S2 16.11 5.0 11.47a 9.64±1.76a 1.95 59.84 K3S2 16.11 5.0 12.77a 8.34±1.59ab 1.71 49.99 K2S0 16.11 0 12.87a 3.24±0.78cd 1.26 26.08 K2S1 16.11 2.5 13.07a 5.54±0.45bcd 1.43 42.74 K2S3 16.11 7.5 12.98a 10.63±0.48a 1.82 82.34 K1S1 16.11 2.5 12.48a 6.13±0.77bc 1.50 50.21 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
表 5 稻草钾肥配施对土壤K平衡的影响1)
Table 5 Influences of combined application of paddy-rice straw and K fertilizer on soil potassium balance
处理
Treatment投入量/(kg·hm−2) Input amount 支出量/
(kg·hm−2)
Output amount表观盈亏量/(kg·hm−2)
Apparent surplus and deficient amount表观平衡系数
Apparent balance coefficient实际平衡率/%
Actual equilibrium rate化学K
Chemical K稻草K
Straw KK0S0 0 0 163.80b −163.80±2.36f 0 −100.00 K0S2 0 86.60 177.01ab −90.40±13.68e 0.49 −50.64 K1S2 83.19 86.60 214.60a −44.81±16.21d 0.80 −20.26 K2S2 166.38 86.60 178.28ab 74.70±23.28b 1.45 45.39 K3S2 249.56 86.60 179.59ab 156.58±28.92a 1.95 94.73 K2S0 166.38 0 186.66ab −20.28±10.75cd 0.90 −10.47 K2S1 166.38 43.30 207.44ab 2.24±5.09c 1.01 1.16 K2S3 166.38 129.90 202.50ab 93.78±4.32b 1.46 46.40 K1S1 83.19 43.30 183.02ab −56.53±13.26de 0.70 −30.41 1)表中数据为3次重复的平均值±标准误，同列数据后的不同小写字母表示差异显著(P<0.05，Duncan’ s法)
1) Data were averages of triplicate samples ± standard deviations (n = 3), and different lowercase letters in the same column indicated significant difference(P<0.05, Duncan’ s test)
下载: 导出CSV
-
[1] 屈冬玉. 马铃薯产业应在我国农村区域经济发展中发挥更大作用[C]//中国马铃薯作物学会马铃薯专业委员会. 马铃薯产业与农村区域发展. 北京: 中国作物学会, 2013: 1-3. [2] 谢开云, 屈冬玉, 金黎平, 等. 中国马铃薯生产与世界先进国家的比较[J]. 世界农业, 2008, 29(5): 35-38. doi: 10.3969/j.issn.1002-4433.2008.05.012 [3] 米健, 罗其友, 高明杰. 南方冬作区马铃薯发展趋势、区域格局与增产潜力[J]. 中国农业资源与区划, 2012, 33(3): 73-79. [4] 屈冬玉, 金黎平, 谢开云. 中国马铃薯产业10年回顾[M]. 北京: 中国农业科学技术出版社, 2010: 66-69. [5] 高媛, 韦艳萍, 樊明寿. 马铃薯的养分需求[J]. 中国马铃薯, 2011, 25(3): 182-187. doi: 10.3969/j.issn.1672-3635.2011.03.018 [6] MULDER A. Potato diseases-diseases, pests and defects[M]. Holland: Aardappelwereld B V and NIVAP, 2005: 209-221.
[7] GEORGE M S, LU G, ZHOU W. Genotypic variation for potassium uptake and utilization efficiency in sweet potato (Ipomoea batatas L.)[J]. Field Crop Res, 2002, 77(1): 7-15. doi: 10.1016/S0378-4290(02)00043-6
[8] 白由路. 高价格下我国钾肥的应变策略[J]. 中国土壤与肥料, 2009, 45(3): 1-4. doi: 10.3969/j.issn.1673-6257.2009.03.001 [9] 陈杨, 樊明寿, 李斐, 等. 氮素营养诊断技术的发展及其在马铃薯生产中的应用[J]. 中国农学通报, 2009, 25(3): 66-71. [10] 张新明, 张洪秀, 李水源, 等. 恩平市典型种植户冬作马铃薯施肥状况调查分析[J]. 安徽农业科学, 2011, 39(36): 22286-22288. doi: 10.3969/j.issn.0517-6611.2011.36.034 [11] LAL R. World crop residues production and implications of its use as a biofuel[J]. Environ Int, 2005, 31(4): 575-584. doi: 10.1016/j.envint.2004.09.005
[12] 毕于运, 王红彦, 王道龙, 等. 中国稻草资源量估算及其开发利用[J]. 中国农学通报, 2011, 27(15): 137-143. [13] 李国学, 张福锁. 固体废物堆肥化和有机复混肥的生产[M]. 北京: 化学工业出版社, 2000: 1-20. [14] 曹国良, 张小曳, 王亚强, 等. 中国区域农田秸秆露天焚烧排放量的估算[J]. 科学通报, 2007, 52(15): 1826-1831. doi: 10.3321/j.issn:0023-074x.2007.15.017 [15] 谭德水, 金继运, 黄绍文, 等. 施钾和秸秆还田对栗钙土区土壤养分及小麦产量的长期效应研究[J]. 干旱地区农业研究, 2009, 27(2): 194-198. [16] 张静, 温晓霞, 廖允成, 等. 不同玉米秸秆还田量对土壤肥力及冬小麦产量的影响[J]. 植物营养与肥料学报, 2010, 16(3): 612-619. doi: 10.11674/zwyf.2010.0314 [17] CARUSO G, CARPUTO D, CONTI S, et al. Effect of mulching and plant density on out-of-season organic potato growth, yield and quality[J]. Adv Hortic Sci, 2013, 27(3): 115-121.
[18] DVORAK P, KUCHTOVA P, TOMASEK J. Response of surface mulching of potato (Solanum tuberosum) on SPAD value, colorado potato beetle and tuber yield[J]. Int J Agric Biol, 2013, 15(4): 798-800.
[19] 农业部办公厅. 农业部关于印发《测土配方施肥技术规范(2011年修订版)》的通知: 农农发[2011]3号[A].(2011-09-22)[2018-11-02]. http://jiuban.moa.gov.cn/zwllm/tzgg/tz/201109/t20110922_2293389.htm. [20] 鲍士旦. 土壤融化分析[M]. 3版. 北京: 中国农业出版社, 2000. [21] 王志勇, 白由路, 杨俐苹, 等. 低土壤肥力下施钾和秸秆还田对作物产量及土壤钾素平衡的影响[J]. 植物营养与肥料学报, 2012, 18(4): 900-906. [22] TANG Q, ZHANG C. Data processing system (DPS) software with experimental design, statistical analysis and data mining developed for use in entomological research[J]. Insect Sci, 2013, 20(2): 254-260. doi: 10.1111/j.1744-7917.2012.01519.x
[23] 闫超. 水稻秸秆还田腐解规律及对土壤养分和酶活性的影响[D]. 哈尔滨: 东北农业大学, 2012. [24] 李世清, 李生秀. 有机物料在维持土壤微生物体氮库中的作用[J]. 生态学报, 2001, 21(1): 136-142. doi: 10.3321/j.issn:1001-9332.2001.01.043 [25] ABRO S A, TIAN X H, WANG X D, et al. Decomposition characteristics of maize (Zea mays L.) straw with different carbon to nitrogen (C/N) ratios under various moisture regimes[J]. Afr J Biotechnol, 2011, 10(50): 10149. doi: 10.5897/AJB10.2261
[26] SUN B, WANG X, WANG F, et al. Assessing the relative effects of geographic location and soil type on microbial communities associated with straw decomposition[J]. Appl Environ Microbiol, 2013, 79(11): 3327-3335. doi: 10.1128/AEM.00083-13
[27] LU P, LIN Y H, YANG Z Q, et al. Effects of application of corn straw on soil microbial community structure during the maize growing season[J]. J Basic Microbiol, 2015, 55(1): 22-32. doi: 10.1002/jobm.201300744
[28] 李逢雨, 孙锡发, 冯文强, 等. 麦秆、油菜秆还田腐解速率及养分释放规律研究[J]. 植物营养与肥料学报, 2009, 15(2): 374-380. doi: 10.3321/j.issn:1008-505X.2009.02.018 [29] 戴志刚, 鲁剑巍, 李小坤, 等. 不同作物还田秸秆的养分释放特征试验[J]. 农业工程学报, 2010, 26(6): 272-276. doi: 10.3969/j.issn.1002-6819.2010.06.047 [30] 张晋京, 窦森. 玉米秸秆分解期间胡敏酸、富里酸动态变化的研究[J]. 土壤通报, 2005(1): 134-136. doi: 10.3321/j.issn:0564-3945.2005.01.037 [31] 张亮. 关中麦玉轮作区施氮对秸秆还田小麦产量和秸秆养分释放的影响[D]. 杨凌: 西北农林科技大学, 2012. [32] 鲁如坤, 时正元, 施建平. 我国南方6省农田养分平衡现状评价和动态变化研究[J]. 中国农业科学, 2000, 33(2): 63-67. doi: 10.3321/j.issn:0578-1752.2000.02.010 [33] 鲁如坤, 刘鸿翔, 闻大中, 等. 我国典型地区农业生态系统养分循环和平衡研究:Ⅳ:农田养分平衡的评价方法和原则[J]. 土壤通报, 1996, 27(5): 197-199. [34] 谭德水, 金继运, 黄绍文, 等. 不同种植制度下长期施钾与秸秆还田对作物产量和土壤钾素的影响[J]. 中国农业科学, 2007, 40(1): 133-139. doi: 10.3321/j.issn:0578-1752.2007.01.018 -
期刊类型引用(3)
1. 国显宝，韩冰，逯新民，何奕成，于露，张万通，段丽辉，耿依仪，邵新庆，俞点，韩伽怡. 氮添加对青藏高原高寒草甸丛枝菌根真菌群落的影响. 草地学报. 2024(07): 2054-2061 . 
百度学术
2. 薛壮壮，冯童禹，王超，沈仁芳. 磷肥对酸性红壤上玉米不同部位丛枝菌根真菌群落的影响. 土壤. 2023(05): 1008-1015 . 百度学术
3. 胡家欣，彭思利，张栋，葛之葳，杨楠. 氮添加对不同林龄杨树人工林丛枝菌根真菌群落的影响. 生态环境学报. 2020(09): 1768-1775 . 百度学术
其他类型引用(10)
计量
- 文章访问数: 1388
- HTML全文浏览量: 3
- PDF下载量: 1346
- 被引次数: 13
