生态学报  2021, Vol. 41 Issue (24): 9815-9823

文章信息

段媛媛, 王淳忆, 张曹莉, 郭富康, 郭正刚
DUAN Yuanyuan, WANG Chunyi, ZHANG Caoli, GUO Fukang, GUO Zhenggang
植株密度和氮素互作对垂穗披碱草生物量分配的影响
Effects of interaction between plant density and nitrogen application rate on allocation of above/underground biomass of Elymus nutans
生态学报. 2021, 41(24): 9815-9823
Acta Ecologica Sinica. 2021, 41(24): 9815-9823
http://dx.doi.org/10.5846/stxb202012283289

文章历史

收稿日期: 2020-12-28
网络出版日期: 2021-08-04
植株密度和氮素互作对垂穗披碱草生物量分配的影响
段媛媛 , 王淳忆 , 张曹莉 , 郭富康 , 郭正刚     
兰州大学草地农业科技学院, 草业科学国家级实验教学示范中心, 兰州 730020
摘要: 垂穗披碱草(Elymus nutans)是高寒地区建植和改良栽培草地的首选草种。虽然合理植株密度和氮素添加量是垂穗披碱草栽培草地稳产的关键因子,但两者之间是否存在最佳互作组合仍不清楚。采用盆栽试验,通过分析不同植株密度(58、102、146株/m2)和氮素添加量(0、200、400 mg/kg)组合状态下垂穗披碱草株高、单株分蘖数、地上生物量、地下生物量、根系体积和地上地下生物量比,以确定理论上是否存在植株密度和氮素添加量的最佳组合。结果表明:随植株密度增加,垂穗披碱草株高、地上生物量和地上地下生物量比值均先增加后降低,而单株分蘖数逐渐减小,根系体积和地下生物量先增加后保持相对稳定;随氮素添加量增加,垂穗披碱草单株分蘖数、地上生物量和地上地下生物量比值均表现为先增加后降低,地下生物量逐渐降低。植株密度与氮素添加量互作虽然对垂穗披碱草的根系体积和单株分蘖数没有显著影响,但两者互作显著影响了垂穗披碱草株高、地上生物量、地下生物量、地上地下生物量比(P < 0.01),这些指标与植株密度和氮素添加量的关系均表现为一个开口向下的抛物面。当植株密度为102株/m2和氮素添加量为200 mg/kg时,垂穗披碱草栽培草地产量最大,生物量分配最优。垂穗披碱草植株密度和氮素添加互作时理论上存在最佳组合,这为垂穗披碱草栽培草地的田间管理提供了理论依据。
关键词: 植株密度    氮素添加    互作    垂穗披碱草    生物量    
Effects of interaction between plant density and nitrogen application rate on allocation of above/underground biomass of Elymus nutans
DUAN Yuanyuan , WANG Chunyi , ZHANG Caoli , GUO Fukang , GUO Zhenggang     
Lanzhou University, National Demonstration Center for Experimental Grassland Science Education, College of Pastoral Agriculture Science and Technology, Lanzhou 730020, China
Abstract: Elymus nutans is the prior crop to establish the pasture or improve the pasture in the alpine regions because it is cold-resistant tolerance with good quality. Although proper plant density and nitrogen application rate are conserved as the key factors to affect the yield stability of Elymus nutans pasture, whether there is an optimally theoretical combination of plant density and nitrogen application rate for maintaining the relative yield of Elymus nutans is not well documented. A pot experiment was conducted in a greenhouse to investigate the effect of different combination of plant densities and nitrogen application rates on the plant height, tiller per plant, aboveground biomass, underground biomass, root volume and aboveground to underground biomass ratio, which can be used to identify an optimally theoretical combination between plant density and nitrogen application rate. In this study, the plant densities consisted of 58, 102, 146 plants per square meter and nitrogen application rate was 0, 200, 400 mg/kg. This study showed that with the increase of plant density, the plant height, aboveground biomass and aboveground to underground biomass ratio of Elymus nutans increased first and then decreased, while the tiller per plant decreased gradually, the root volume and underground biomass increased first and then remained relatively stable. The tiller per plant, aboveground biomass and aboveground to underground biomass ratio of Elymus nutans increased first and then decreased with the increase of nitrogen application rates, and the underground biomass decreased gradually. Interaction between plant densities and nitrogen application rates was found to had no effect on root volume and tiller per plant, whereas it was observed to have significant effects on the plant height, aboveground biomass, underground biomass, and the aboveground to underground biomass ratio. The relationships between above-mentioned parameters (plant height, aboveground biomass, underground biomass, and the aboveground to underground biomass ratio) and nitrogen application rate and plant density showed the downward 3D response surfaces, and this demonstrated that there was an optimally theoretical combination between proper plant density and nitrogen application rate, under which plant height, aboveground biomass, underground biomass, and the aboveground to underground biomass ratio were the best. In this study, the optimal combination of plant density and nitrogen application rate were 102 plants per square meter and the nitrogen application rate was 200 mg/kg. The findings of this study suggest that there is an optimally theoretical combination of plant density and nitrogen application rate for maintaining the relative yield of Elymus nutans pasture, which can provide a scientific basis for high persistence yield and optimal allocation of above/underground biomass of Elymus nutans pasture in the field management.
Key Words: plant density    nitrogen application rate    interaction    Elymus nutans    biomass    

垂穗披碱草(Elymus nutans)是一种优质多年生禾本科牧草, 因其具有耐寒性[1]已成为我国高寒地区建植栽培草地的首选草种之一[2-3]。随着高寒地区退化草地改良和舍饲畜牧业发展, 垂穗披碱草栽培草地面积未来将会持续增加[4], 因此, 维持垂穗披碱草栽培草地稳产, 已成为目前高寒地区草牧业持续发展亟待解决的重大科学问题。

植株密度和土壤养分因影响植株间的竞争而成为栽培草地稳产的两个重要因素[1, 3, 5-6]。植株密度是植物种群内部个体特征的主要约束因素之一[7], 植株密度过低会降低栽培草地产量和品质[8], 植株密度过大则会引起植株内对光和养分资源的竞争激烈[9-10], 不仅不会增加产量[1], 反而会增加建植成本。氮是高寒地区植物生长的限制性元素[3, 11-12], 因此, 添加氮素成为维持栽培草地稳产的主要管理策略[2-3, 13]。然而氮素添加量过小并不能实现促进植物生长的功效[14], 氮素添加过量会通过影响植物营养物质代谢而降低光合速率[1, 9, 15], 而且会增加土壤非点源污染的风险[9, 13, 16-19]

垂穗披碱草产量与植株密度和氮素添加间的关系虽然已有大量研究, 但尚未形成共识。关于植株密度与垂穗披碱草生物量的关系, 有研究认为呈单峰关系, 当播种密度为237粒/m2时其产量最大[6], 然而有研究则认为垂穗披碱草产量与播种密度呈负相关关系, 表现为植株密度从55粒/m2增加到167粒/m2时, 产量逐渐降低[7], 这说明植株密度与垂穗披碱草生物量关系的分异主要来自于植株密度范围的设置。而垂穗披碱草产量与氮素添加量的关系亦出现分异, 有研究认为是正相关关系[3], 即当施氮量从0增加到50 g/m2时垂穗披碱草产量逐渐增加, 也有研究表明是单峰关系[7], 当氮素添加量为150 kg/hm2时垂穗披碱草产量最高。研究表明, 栽培草地内植株密度与土壤氮素之间存在着密切关联[7], 这种关联有时表现为协同, 有时则表现为拮抗[20]。管理栽培草地的核心目标就是将植株密度与氮素添加量调整至协同状态, 实现两者之间的耦合[9], 最终维持草地稳产和氮的最优化利用。虽然播种密度和氮素添加互作能够显著影响垂穗披碱草的生长及生殖性能[7], 但两个因素互作是否存在有助于生物量稳定的最佳理论组合, 尚需科学试验提供证据因此, 继续深入研究植株密度和氮素添加互作对垂穗披碱草产量及其地上地下生物量分配的影响, 确证植株密度和氮素添加量互作时是否存在理论的最佳组合, 是垂穗披碱草栽培草地管理时保留合理植株密度和添加适宜氮素的基础性问题, 具有重要的现实和理论意义。

1 材料与方法 1.1 试验设计

试验采用双因素设计, 两个因素为植株密度和氮素添加量。植株密度和氮素添加量参照已有研究, 植物密度设置为20、35、50株/盆[21], 即58、102、146株/m2分别用D1、D2、D3表示, 氮素添加量设置为0 mg/kg、200 mg/kg、400 mg/kg[4], 分别用N1、N2、N3表示。氮素添加量以纯氮含量为准, 氮源为含氮量为46% 的尿素。试验共9个处理(D1N1、D1N2、D1N3、D2N1、D2N2、D2N3、D3N1、D3N2、D3N3), 每个处理设置3个重复。试验以盆栽方式在兰州大学榆中校区智能温室进行, 温室温度维持在15—25 ℃, 光照强度300 lx。花盆选用聚乙烯仿瓷塑料花盆, 高28.5 cm, 口径33 cm, 底径20 cm,试验时每个花盆装土10kg。本研究采用的土壤为黑垆土, 采集于甘肃榆中县夏官营地区的耕地耕作层, 其土壤pH值为7.4, 有机质含量0.66%, 全氮含量0.11%, 速效磷57.66 mg/kg, 速效钾240 mg/kg。

2018年10月20日, 首先测定试验用土的饱和含水量[22];其次, 按照处理设计核算每个花盆需要添加的尿素量, 将尿素和10 kg土壤充分搅拌混合, 一起装入该花盆, 不再添加其他任何营养元素, 共计装盆27个, 满足每个处理3个重复的设计需求;第三, 所有花盆装土和尿素后, 缓慢浇水, 将土壤含水量调至饱和含水量的60%左右;第四, 挑选饱满均一的垂穗披碱草种子, 于10月21日播种于花盆, 每盆播种100粒;第五, 为避免土壤水分含量变化对种子萌发和幼苗生长的影响, 采用称重法每3天浇一次水, 让花盆内土壤含水量始终维持于60% 左右, 此含水量是垂穗披碱草适宜生长的水分含量[4, 21]。2018年11月20日, 垂穗披碱草完整出苗时, 采用人工间苗法将植株密度调至试验设计水平, 间苗时每个盆内尽可能的保留长势均匀的植株。此后每周随机移动花盆位置, 保证各盆植株受光均匀, 每个盆的管理措施一致。

1.2 取样与指标测定

2019年9月10日取样。首先, 每个花盆中随机选取5株, 采用直尺(mm)测定其自然株高, 在根颈处查数单株分蘖数, 5株的平均值视为该盆植株的株高和分蘖数。再将每盆中全部植株齐土壤表面刈割, 获取地上生物量鲜样。之后将花盆打碎, 用自来水洗出植物根系后用蒸馏水冲洗干净并剔除杂物, 最后用干净滤纸吸干根系表面水分后, 将根系放入盛有一定水量的量筒内, 用玻璃棒轻轻搅动排出空气, 当根系完全浸泡在水中时, 获取量筒内水增加的体积数, 即为根系体积。最后, 将植物地上部分鲜样和根系, 均放置于65 ℃的烘箱, 烘48 h至恒重, 冷却至室温后用精确度为0.001 g的电子天平称量, 获得地上生物量和地下生物量。生物量分配采用地上生物量和地下生物量比值评价。

1.3 数据处理

先采用One-Way ANOVA检验各个指标数据的方差齐性, 若不一致, 则采用Log转置数据, 再采用Two-Way ANOVA分析植株密度和氮素添加量对各个指标的影响。若植物密度和氮素添加量互作不显著, 则不做进一步分析;若植株密度和氮素添加量互作对某一个指标的影响显著, 继续采用Tukey比较法进行多重比较分析, 最后采用MATLAB拟合该指标与植株密度和氮素添加量间的关系, 建立模型。

2 结果 2.1 植株密度和氮素添加量互作对垂穗披碱草地上生物量及构成要素的影响

植株密度显著影响了垂穗披碱草的株高、单株分蘖数和地上生物量(表 1)。氮素添加量显著影响了垂穗披碱草的单株分蘖数和地上生物量(P < 0.01), 但没有显著影响株高。植株密度和氮素添加量互作显著影响了株高和地上生物量, 对单株分蘖数没有显著影响。

表 1 植株密度和氮素添加量互作对垂穗披碱草株高、单株分蘖数及地上生物量的影响(平均值±标准误差) Table 1 Effects of plant density and nitrogen application rate on plant height, tillers per plant and aboveground biomass of Elymus nutans(Mean±SE)
处理
Treatments
地上生物量
Aboveground biomass/ (g/m2)
株高/cm
Plant height
单株分蘖数
Tillers per plant/ (枝条数/株)
处理
Treatments
地上生物量
Aboveground biomass/ (g/m2)
株高/cm
Plant height
单株分蘖数
Tillers per plant/ (枝条数/株)
D1 48.02±13.92b 33.15±2.07b 7.10±0.11a D1×N1 31.20±1.46d 32.73±2.71b 6.38±0.12
D2 70.10±24.65a 35.24±2.72a 6.41±0.11b D1×N2 49.86±2.81c 34.50±1.15ab 7.61±0.44
D3 47.90±12.87b 30.27±2.76c 4.41±0.11c D1×N3 63.05±2.25b 32.21±2.05b 7.31±0.48
N1 36.35±2.72b 33.00±0.69 5.18±1.34c D2×N1 47.08±3.39c 33.63±1.73ab 5.70±0.33
N2 65.80±27.14a 33.35±0.64 6.64±1.25a D2×N2 101.80±2.25a 38.03±2.68a 7.30±0.12
N3 37.58±7.31c 32.30±0.64 6.10±1.17b D2×N3 61.44±3.48b 34.06±1.30c 6.24±0.43
D3×N1 34.45±0.79d 32.67±2.61b 3.47±0.20
D3×N2 45.68±2.69c 27.50±0.86ab 5.02±0.20
D3×N3 63.58±2.92b 30.63±1.64bc 4.75±0.26
不同小写字母代表不同处理间存在显著差异(P < 0.05); D1、D2、D3表示三种植株密度处理, 代表的植株密度分别为20、35、50株/盆(58、102、146株/m2); N1、N2、N3表示三种氮素添加量处理, 代表的氮素添加量分别为0、200、400mg/kg

随植株密度增加, 垂穗披碱草株高和地上生物量先增加后降低, 植物密度为102株/m2时最大, 其中株高较58株/m2和146株/m2时分别增加了6%和16%, 地上生物量增加了48%和46%, 而单株分蘖数随植株密度增加逐渐减小, 降幅介于61% 到45%。随氮素添加量增加, 垂穗披碱草单株分蘖数和地上生物量均表现为先增加后降低, 氮素添加量为200 mg/kg时最大, 单株分蘖数较0 mg/kg和400 mg/kg时增加28% 和8%, 地上生物量较0 mg/kg和400 mg/kg时增加81%和79%。

垂穗披碱草的株高F(X1, X2) 与植株密度(X1)、氮素添加量(X2) 之间的数学拟合方程为F(X1, X2)=19.6546 +0.9787X1+0.0042X2-0.0157X12-0.000017X22-0.00013X1X2, F 检验时, P=0.0212, R2=0.7701, 表明在一定范围内, 植株密度增加垂穗披碱草株高的效应大于氮素添加的效应(图 1), X1X2的二次项系数均为负数, 表明该面为开口向下的抛物面, 说明垂穗披碱草株高对植株密度和氮素添加互作的响应, 理论上存在植株密度和氮素添加量的最佳组合。试验设计的9个处理中, 最佳组合是植株密度为102株/m2, 氮素添加量为200 mg/kg, 此时垂穗披碱草株高较大。

图 1 植株密度和氮素添加量互作下垂穗披碱草株高曲面图 Fig. 1 3D response surface of plant height of Elymus nutans between plant density and nitrogen application rate

垂穗披碱草的地上生物量F(X1, X2) 与植株密度(X1)、氮素添加量(X2) 间的数学拟合方程为F(X1, X2)=-23.8003+2.3687 X1+0.0777 X2-0.0336 X12-0.00013 X22-0.000077 X1 X2, F 检验时, P=0.00012, R2=0.0.6788, 表明在一定范围内, 植株密度增加垂穗披碱草地上生物量的效应大于氮素添加的效应(图 2), X1X2的二次项系数均为负数, 表明该面为开口向下的抛物面, 说明理论上植株密度和氮素添加量存在最佳组合, 此时垂穗披碱草的地上生物量最大。9个处理中, 102株/m2的植株密度和200 mg/kg的氮素添加量为最佳组合, 此组合下垂穗披碱草的地上生物量较大。

图 2 植株密度和氮素添加量互作下垂穗披碱草地上生物量曲面图 Fig. 2 3D response surface of aboveground biomass of Elymus nutans between plant density and nitrogen application rate
2.2 植株密度和氮素添加量互作对垂穗披碱草根系生物量及体积的影响

植株密度极显著影响了垂穗披碱草的根系体积和地下生物量(表 2), 而氮素添加量、植株密度和氮素添加量互作仅显著影响了地下生物量(P < 0.01), 对根系体积没有显著影响。随植株密度从58株/m2增加到102株/m2时, 垂穗披碱草根系体积和地下生物量显著增加, 分别增加13% 和26%, 而当植株密度超过102株/m2后, 根系体积和地下生物量基本保持相对稳定。随氮素添加量增加, 垂穗披碱草地下生物量逐渐降低, 降幅介于6% 到25%。

表 2 植株密度和氮素添加量对垂穗披碱草根系体积和地下生物量的影响 Table 2 Effects of plant density and nitrogen application rate on root volume and underground biomass of Elymus nutans
处理
Treatments
根系体积
Root volume/ (cm3/m2)
地下生物量
Underground biomass/ (g/m2)
D1 50.80±8.01b 10.94±2.34b
D2 57.55±5.56ab 13.74±1.58a
D3 62.50±3.63a 13.16±0.73a
N1 54.51±1.87 13.80±0.56a
N2 59.95±1.87 12.98±0.82b
N3 56.35±1.87 11.08±2.34c
D1×N1 47.67±11.11 13.54±1.11bc
D1×N2 57.23±5.21 11.05±0.56e
D1×N3 47.52±3.66 8.19±0.50f
D2×N1 53.58±3.39 14.33±0.18ab
D2×N2 58.20±0.29 15.09±0.73a
D2×N3 60.86±8.19 11.79±0.70de
D3×N1 62.29±1.75 13.48±0.67bc
D3×N2 64.43±5.38 12.78±0.82cd
D3×N3 60.74±3.16 13.25±0.85bc

垂穗披碱草的地下生物量F(X1, X2) 与植株密度(X1)、氮素添加量(X2) 之间的关系为F(X1, X2)=-0.8772+0.2340X1+0.0093X2-0.0026X12-0.0000046X22-0.00014X1X2, F 检验时, P=0.00012, R2=0.8701, 表明在一定范围内, 植株密度增加垂穗披碱草地下生物量的效应要大于氮素添加的效应(图 3), 空间分布表现为一个开口向下的抛物面, 说明垂穗披碱草的地下生物量响应植株密度和氮素添加互作时, 其理论上存在最佳组合。当植株密度为102株/m2和氮素添加量为200 mg/kg时, 垂穗披碱草地下生物量较大。

图 3 植株密度和氮素添加量互作下垂穗披碱草地下生物量曲面图 Fig. 3 3D response surface of underground biomass of Elymus nutans between plant density and nitrogen application rate
2.3 植株密度和氮素添加量互作对垂穗披碱草地上地下生物量比值的影响

植株密度、氮素添加量, 以及植株密度和氮素添加互作均极显著影响了垂穗披碱草地上地下生物量比(P < 0.01)。随植株密度增加, 垂穗披碱草的地上地下生物量比值均先增加后降低, 在102株/m2时达到最大比值4.89, 较58株/m2和146株/m2的增幅为13% 和27%;随氮素添加量增加, 垂穗披碱草的地上地下生物量比值也呈现先增加后降低的趋势, 在200mg/kg时有最大值4.97, 较0 mg/kg和400 mg/kg高9% 和42%(表 3)。

表 3 植株密度和氮素添加互作对垂穗披碱草地上地下生物量比值的影响 Table 3 Effects of plant density and nitrogen application rate on aboveground biomass to underground biomass of Elymus nutans
处理
Treatment
地上地下生物量比值Aboveground biomass to underground biomass ratio
N1 N2 N3 主效应植株密度
Main effect of plant density
D1 4.65±0.13b 4.52±0.22b 3.82±0.42b 4.33±0.46b
D2 4.28±0.27b 6.39±0.54a 4.00±0.24b 4.89±1.78a
D3 4.73±0.44b 6.39±0.54b 2.65±0.22c 3.80±1.00c
主效应氮素添加量
Main effect of nitrogen application rate
4.55±0.55b 4.97±0.44a 3.49±0.22c

垂穗披碱草的地上地下生物量比F(X1, X2) 与植株密度(X1)、氮素添加量(X2) 的数学关系为F(X1, X2)=0.8745+0.2197X1+0.0085X2-0.0037X12-0.000024X22 -0.00010 X1X2, F 检验时, P=0.00022, R2=0.6586, 表明在一定范围内植株密度增加垂穗披碱草地上地下生物量比值的效应要大于氮素添加的效应(图 4), X1X2二次项系数均为负数, 则说明该面为一个开口向下的抛物面, 反映了垂穗披碱草的地上地下生物量比对植株密度和氮素添加互作的响应, 理论上存在植株密度和氮素添加量的最佳组合。本试验的9个处理中, 植株密度为102株/m2和氮素添加量为200 mg/kg处理下, 垂穗披碱草地上地下生物量比较大。

图 4 植株密度和氮素添加互作下垂穗披碱草地上地下生物量比值曲面图 Fig. 4 3D response surface of the aboveground to underground biomass ratio of Elymus nutans between plant density and nitrogen application rate
3 讨论

栽培草地产量和物质分配不仅受植株密度和氮素添加单个因素的影响, 而且受两个因素互作效应的影响, 而植物密度和氮素添加单因素的效应是分析互作效应的基础。

植株密度与垂穗披碱草地上生物量关系密切, 主要是植株密度变化时改变了植株个体获取光照资源的能力和植株个体叶面积指数[10], 这导致植株个体间的合作和竞争关系不断调整[6], 从而影响地上和地下生物量及分配。当植株密度较低时, 植株个体间间隙较大, 单位面积的叶面积指数较低, 光能利用率低[23], 地上生物量维持在较低水平;当植株密度增加时, 垂穗披碱草植株个体间的竞争关系逐渐增加[6], 植株个体间的空隙被逐渐占据, 单位面积叶面积指数逐渐增加[5, 23-25], 光合作用增强, 光合产物增多[1, 14], 此时植株间合作关系占据优势地位, 而这种合作关系增加了土壤温度、水分和氮素的有效性[6, 25], 从而促进垂穗披碱草的地上生物量逐渐增加;当植株密度过高时, 虽然植物群体的叶面积指数较大, 但植株个体间互相遮荫, 透光效果差[25-26], 下层叶片接收不到光照, 整体光合作用下降[23], 导致垂穗披碱草个体间对有限资源的高度竞争、抵消并最终压倒了促进作用[6], 从而导致垂穗披碱草的地上生物量下降, 因此垂穗披碱草只有维持在合理植株密度范围内时[5], 其地上生物量才会最优。株高和分蘖数是垂穗披碱草地上产量的主要构成要素, 两者共同决定地上生物量的高低[1]。垂穗披碱草株高随植株密度增加先增加后降低, 主要与植株个体间竞争关系的变化有关[2, 6]。当植株密度从低到高增加时, 植株个体为获取光资源维持光合作用[1], 对光资源竞争加剧[10, 27], 每个植株个体均会努力向上竞争生长[26], 以获取更多光照资源[5, 21], 客观上形成植株高度逐渐增加的局面[16, 28]。当植株密度过大后, 植株个体之间竞争极强烈, 资源的有限性导致植株个体能够获取的资源量减少[6, 24-25], 从而抑制植株生长[27], 此时则会形成植株株高反而有所降低的现象, 这与植株密度与植株高度关系的研究结果一致[10]。垂穗披碱草单株分蘖数随植物密度增加表现为逐渐降低, 这与已有垂穗披碱草密度与分蘖数的研究结果趋同[2, 7], 主要原因是随植株密度增大, 植株个体占据的资源量减少, 光截获量下降[25-27], 一方面减弱了光照对单株分蘖数的正向增加作用, 另一方面资源对营养繁殖的投入下降所致[2, 21, 27], 这与个体间种内竞争不利于单株分蘖数增加的观点一致[2]。研究结果表明, 垂穗披碱草根系体积和地下生物量随植株密度增加表现为先增加后保持相对稳定的过程, 是由于植株密度增大时, 植株获取资源时受到邻株在物理和几何上的拥挤效应[24, 26], 促使植物主动调整状态[1], 此时植物个体需要获取并利用更多地下资源, 植物根系的觅食行为会促使根系伸长生长加强, 并构建“扩散型”根系系统, 降低内部根系拥挤程度以充分利用资源, 提高自身生境适应度[1], 因此垂穗披碱草的根系体积和生物量逐渐增加。但随植株密度持续增大, 垂穗披碱草个体间竞争加剧[2, 6], 个体根系生长减弱, 但群体根系优势依然不断增大[1]。高植株密度时, 植物根系依然存在觅食效应, 但由于花盆空间限制, 植株的根系密度和根系深度达到饱和状态[19], 此时根系体积和生物量会维持在相对稳定的状态, 不再增加。研究表明, 垂穗披碱草地上地下生物量比随植株密度增加先增加后降低, 这与植物获取资源的能力和根系生物量分配的权衡有关[1, 28]。植株密度开始增加时, 由于植物的权衡特性会将有限的资源进行分配以维持各个器官的生理机能[2, 7]。垂穗披碱草为获取更多光资源生产更多有机物质, 一般优先将光合产物分配于地上生长[1, 26], 此时垂穗披碱草地上地下生物量比值增大;而当植株密度超过一定阈值后, 植物个体间地上部分相互遮荫, 整体光合产物有所减少[26], 但地下生物量却一直维持在较高水平(表 2), 此时表现为地上地下生物量比值降低。

土壤氮素变化对垂穗披碱草产量具有明显影响, 主要原因是垂穗披碱草喜嗜氮, 因此其生物量对土壤氮素变化反应敏感[1]。垂穗披碱草单株分蘖数和地上生物量随氮素添加量增加而先增加后降低, 这与玉米(Zea mays)和狼尾草(Pennisetum alopecuroides)地上生物量响应氮素变化的结果一致[18, 26]。氮素添加量较低时, 一定程度上改善了土壤整体营养状况, 植物所需养分得到及时补充[4, 13, 20, 24, 27];同时氮素刺激了植物生长和分蘖发生[1, 5, 7, 15], 有利于植物叶面积增加[1, 24, 27], 延缓叶片衰老[14, 27], 延长叶片光合作用时间增加光合产物[9-10, 24, 27], 地上生物量随之增加, 这一过程在玉米和棉花(Gossypium spp)上也到了验证[15, 18]。随氮素添加量增加, 垂穗披碱草分蘖能力继续增加, 但增加速率有所降低, 此时单株分蘖数和地上生物量依然处于缓慢增加状态。当氮素添加过量时, 垂穗披碱草生长旺盛, 植物个体捕获的光资源减少且叶片光合速率在高氮下会显著降低[1, 9, 25, 27];在拥挤胁迫下, 叶片作为吸收氮素的主要贮藏器官, 叶片氮浓度会有所下降[25], 影响植物氮素吸收利用[13-14, 16, 18, 29], 植物分蘖受到约束[5], 导致单株分蘖数降低;同时土壤高氮能够刺激植物体内硝酸盐富集, 抑制植物生长减缓叶片中光合产物的有效运输[1], 最终表现为地上生物量降低, 这一现象与狼尾草对过量氮素的响应趋同[29]。然而, 垂穗披碱地下生物量随氮素添加量增加却逐渐降低, 一方面是由于土壤氮素添加量的升高, 影响了植株光合产物的分配, 根系的分配量降低[17];一方面是高氮素使分配到植株根系的碳素主要用于转运氮素, 便于植物吸收利用, 从而导致植株根系生长受阻[30]。垂穗披碱草地上地下生物量比随氮素添加量增加表现为逐渐增加, 是由于随氮素添加量增加, 植物群体生长的限制因素从土壤矿物质资源变为光资源[30], 依据最优分配原理和植物的生活史策略[31], 植物有获得最大适合度的适应性行为[2, 7, 30], 此时植物会将更多光合产物优先分配于获取限制性资源的地上器官[1, 7, 27], 最大限度地获取光资源, 保证群体健康生长, 客观上形成了地下生物量维持在较低水平, 表现为地上生物量降低。

植株密度和氮素变化不仅明显影响垂穗披碱草的产量和物质分配, 两者的互作更会影响垂穗披碱草的产量和物质分配。研究结果表明, 植株密度和氮素添加量互作对垂穗披碱草不同指标的影响存在分异特征, 虽然植株密度与氮素添加互作对单株分蘖数和根系体积没有显著影响, 但两者互作却显著影响了垂穗披碱草株高、地上生物量、地下生物量、地上地下生物量比。不论是株高和地上生物量, 还是地下生物量和地上地下生物量比, 与植株密度和氮素添加量的关系均表现为一个开口向下的抛物面, 说明垂穗披碱草株高、地上生物量、地下生物量、地上地下生物量比响应植株密度和氮素添加量互作时, 理论上存在最佳组合, 并不是植株密度越大越好, 也不是氮素添加量越大越好。只有当植株密度和氮素添加量处于合适状态时, 两者间关系则会表现为耦合状态, 此时垂穗披碱草产量和地上地下生物量分配最佳, 既能实现高产目标, 又能节约氮素、降低建植成本, 而且适量施氮会使收获期土壤硝态氮与铵态氮含量相对较低, 减少环境污染[18], 过度施氮则会导致养分过度代谢、光合速率降低、同化物积累, 降低氮肥的生理和农艺效率[9]。本研究设计的9个组合处理中, 植株密度为102株/m2和氮素添加量为200 mg/kg时, 垂穗披碱草的地上和地下生物量最大, 且最有利于生物量向地上配置, 此组合下垂穗披碱草的株高较氮素添加量为0 mg/kg和400 mg/kg时分别增加了13% 和12%, 地上生物量增加116% 和65%, 地下生物量增加5% 和28%, 地上地下生物量比增加49% 和60%。同时最佳组合时株高(102株/m2和氮素添加量为200 mg/kg)较20株/盆和50株/盆时的株高, 分别增加了10% 和38%, 地上生物量分别增加了104% 和123%, 地下生物量分别增加37% 和18%, 地上地下生物量比分别增加41% 和0。虽然本研究证实了植株密度和氮素添加量互作理论上存在最佳组合, 此组合会维持垂穗披碱草的高产, 而且会优化生物量的配置, 然而结果仅为一年盆栽试验所得, 田间试验条件下植物密度和氮素添加量的互作可能受多种因素综合作用, 因此未来需要采用田间试验继续验证植株密度和氮素添加量互作存在最佳组合的可靠性。同时, 不同地区土壤基础氮素含量存在差异, 且适宜的植株密度可能不同, 生产实践中应该根据区域环境和气候特征, 确定适宜的植株密度和氮素添加量, 从而实现两者的互作协同, 以适宜的植株密度和氮素添加量作为当地垂穗披碱草管理的基准, 提高栽培草地的管理水平, 维持垂穗披碱草栽培草地的高产和稳产。

致谢: 感谢文雅、冯甘霖等同学帮助取样, 感谢庞晓攀、张静和王倩博士帮助数据分析。
参考文献
[1]
李雪萍, 赵成章, 任悦, 张晶, 雷蕾. 尕海湿地不同密度条件下垂穗披碱草根系分形结构. 生态学报, 2018, 38(4): 1176-1182.
[2]
朱志红, 刘建秀, 郑伟. 资源获得性和种内竞争对垂穗披碱草生长繁殖的影响. 西北植物学报, 2005(10): 2056-2061. DOI:10.3321/j.issn:1000-4025.2005.10.021
[3]
罗文蓉, 栗文瀚, 干珠扎布, 闫玉龙, 李钰, 曹旭娟, 何世丞, 旦久罗布, 高清竹, 胡国铮. 施氮对藏北垂穗披碱草人工草地叶片功能性状和种群特征的影响. 草业学报, 2018, 27(5): 51-60.
[4]
文雅, 周培, 张忠雪, 段媛媛, 冯甘霖, 邓雨森, 郭正刚. 施氮和灌溉互作对垂穗披碱草生物量及构成要素的影响. 草业科学, 2020, 37(2): 330-338.
[5]
Ghiasy-Oskoee M, AghaAlikhani M, Sefidkon F, Mokhtassi-Bidgoli A, Ayyari M. Blessed thistle growth, essential oil content, yield and composition as influenced by plant density and nitrogen fertilizer. Journal of Essential Oil Bearing Plants, 2020, 23(2): 276-291. DOI:10.1080/0972060X.2020.1747553
[6]
牛建伟, 雷占兰, 周华坤, 刘泽华, 黄瑞灵, 姚步青, 王文颖, 赵新全. 种植密度和施氮水平对垂穗披碱草生物量分配的影响. 草业科学, 2014, 31(7): 1343-1351.
[7]
雷占兰, 周华坤, 刘泽华, 黄瑞灵, 王文颖. 密度氮肥交互处理下垂穗披碱草生长与生殖特性. 中国草地学报, 2014, 36(2): 12-19. DOI:10.3969/j.issn.1673-5021.2014.02.003
[8]
Omadewu L I, Iren O B, Eneji A E. Yield of cotton cultivars as influenced by nitrogen rates and plant density in Yalingo, Nigeria. World Scientific News, 2019, 127(3): 106-122.
[9]
Zhao Y N, Huang Y F, Li S, Chu X, Ye Y L. Improving the growth, lodging and yield of different density-resistance maize by optimising planting density and nitrogen fertilisation. Plant, Soil and Environment, 2020, 66(9): 453-460. DOI:10.17221/178/2020-PSE
[10]
Labra M H, Struik P C, Calderini D F, Evers J B. Leaf nitrogen traits in response to plant density and nitrogen supply in oilseed rape. Agronomy, 2020, 10(11): 1780. DOI:10.3390/agronomy10111780
[11]
Han Y H, Dong S K, Zhao Z Z, Sha W, Li S, Shen H, Xiao J N, Zhang J, Wu X Y, Jiang X M, Zhao J B, Liu S L, Dong Q M, Zhou H K, Yeomans J C. Response of soil nutrients and stoichiometry to elevated nitrogen deposition in alpine grassland on the Qinghai-Tibetan Plateau. Geoderma, 2019, 343: 263-268. DOI:10.1016/j.geoderma.2018.12.050
[12]
Yang X X, Dong Q M, Chu H, Ding C X, Yu Y, Zhang C P, Zhang Y F, Yang Z Z. Different responses of soil element contents and their stoichiometry (C: N: P) to yak grazing and Tibetan sheep grazing in an alpine grassland on the eastern Qinghai-Tibetan Plateau'. Agriculture, Ecosystems & Environment, 2019, 285: 106628.
[13]
Maltese N E, Maddonni G A, Melchiori R J M, Ferreyra J M, Caviglia O P. Crop nitrogen status of early- and late-sown maize at different plant densities. Field Crops Research, 2020, 258: 107965. DOI:10.1016/j.fcr.2020.107965
[14]
Luo Z, Liu H, Li W P, Zhao Q, Dai J L, Tian L W, Dong H Z. Effects of reduced nitrogen rate on cotton yield and nitrogen use efficiency as mediated by application mode or plant density. Field Crops Research, 2018, 218: 150-157. DOI:10.1016/j.fcr.2018.01.003
[15]
Dai J L, Li W J, Zhang D M, Tang W, Li Z H, Lu H Q, Kong X Q, Luo Z, Xu S Z, Xin C S, Dong H Z. Competitive yield and economic benefits of cotton achieved through a combination of extensive pruning and a reduced nitrogen rate at high plant density. Field Crops Research, 2017, 209: 65-72. DOI:10.1016/j.fcr.2017.04.010
[16]
Wang X K, Wang G, Turner N C, Xing Y Y, Li M T, Guo T. Determining optimal mulching, planting density, and nitrogen application to increase maize grain yield and nitrogen translocation efficiency in Northwest China. BMC Plant Biology, 2020, 20(1): 282. DOI:10.1186/s12870-020-02477-2
[17]
Miller A J, Cramer M D. Root nitrogen acquisition and assimilation. Plant and Soil, 2005, 274(1/2): 1-36.
[18]
刘志恒, 徐开未, 王科, 解晋, 王佳锐, 赵亚妮, 陈远学. 不同施氮量对玉米产量及各器官养分积累的影响. 浙江大学学报: 农业与生命科学版, 2018, 44(5): 573-579.
[19]
Dong S X, Zhang J, Zha T, Dai X L, He M R. Increased plant density with reduced nitrogen input can improve nitrogen use efficiency in winter wheat while maintaining grain yield. Archives of Agronomy and Soil Science, 2020, 66(12): 1707-1720. DOI:10.1080/03650340.2019.1690139
[20]
Xu X L, Ouyang H, Richter A, Wanek W, Cao G M, Kuzyakov Y. Spatio-temporal variations determine plant-microbe competition for inorganic nitrogen in an alpine meadow. Journal of Ecology, 2011, 99(2): 563-571.
[21]
冯甘霖, 文雅, 段媛媛, 郭正刚. 灌溉量和密度对垂穗披碱草生长性能和物质分配的影响. 草业科学, 2019, 36(8): 2087-2095.
[22]
许桂玲, 罗强鑫, 王晓珂, 彭金凤, 冯跃华. 控制水稻盆栽土壤水分的新方法. 农业科技通讯, 2019(12): 219-221.
[23]
刘伟, 吕鹏, 苏凯, 杨今胜, 张吉旺, 董树亭, 刘鹏, 孙庆泉. 种植密度对夏玉米产量和源库特性的影响. 应用生态学报, 2010, 21(7): 1737-1743.
[24]
Kuai J, Sun Y Y, Zhou M, Zhang P P, Zuo Q S, Wu J S, Zhou G S. The effect of nitrogen application and planting density on the radiation use efficiency and the stem lignin metabolism in rapeseed (Brassica napus L.). Field Crops Research, 2016, 199: 89-98. DOI:10.1016/j.fcr.2016.09.025
[25]
Yan P, Zhang Q, Shuai X F, Pan J X, Zhang W J, Shi J F, Wang M, Chen X P, Cui Z L. Interaction between plant density and nitrogen management strategy in improving maize grain yield and nitrogen use efficiency on the North China Plain. The Journal of Agricultural Science, 2016, 154(6): 978-988. DOI:10.1017/S0021859615000854
[26]
黎磊, 周道玮, 盛连喜. 密度制约决定的植物生物量分配格局. 生态学杂志, 2011, 30(8): 1579-1589.
[27]
Ciampitti I A, Vyn T J. A comprehensive study of plant density consequences on nitrogen uptake dynamics of maize plants from vegetative to reproductive stages. Field Crops Research, 2011, 121(1): 2-18. DOI:10.1016/j.fcr.2010.10.009
[28]
Shah A N, Yang G Z, Tanveer M, Iqbal J. Leaf gas exchange, source-sink relationship, and growth response of cotton to the interactive effects of nitrogen rate and planting density. Acta Physiologiae Plantarum, 2017, 39(5): 119. DOI:10.1007/s11738-017-2402-0
[29]
李韦柳, 覃维治, 熊军, 韦民政, 唐秀桦, 闫海锋. 施氮对狼尾草在南方贫瘠旱坡地生长、能源品质及氮肥利用率的影响. 草业学报, 2018, 27(3): 144-153.
[30]
景明慧, 贾晓彤, 张运龙, 曹娟, 周伟, 王杰, 庾强. 长期氮添加对内蒙古典型草原植物地上、地下生物量及根冠比的影响. 生态学杂志, 2020, 39(10): 3185-3193.
[31]
Zheng Y, Chen L, Luo C Y, Zhang Z H, Wang S P, Guo L D. Plant identity exerts stronger effect than fertilization on soil arbuscular mycorrhizal fungi in a sown pasture. Microbial Ecology, 2016, 72(3): 647-658. DOI:10.1007/s00248-016-0817-6