Meta Analysis on Effects of N-fertilization on Yield and N-uptake of Cotton Plants
-
摘要:目的 定量评价棉花产量和氮素吸收量对施氮的响应及其影响因素，为科学施用氮肥提供参考。方法 采用Meta分析的方法，对2002–2019年符合研究要求的43篇文献数据按照种植时段、种植区域、土壤类型、种植密度、施氮量、土壤有机质含量、土壤全氮含量、土壤碱解氮含量、土壤速效磷含量、土壤速效钾含量等进行分组，以不施氮肥为对照，探讨施氮对棉花产量和氮素吸收量影响的综合效应。结果 在2008–2013年棉花产量增量幅度（36%）与氮素吸收量提升幅度（51%）均显著高于2002–2007年（20%、34%）和2014–2019年（28%、33%）。施氮量为300–450 kg·hm−2时棉花产量和氮素吸收量（38%、58%）显著高于施氮量>450 kg·hm−2（30%、54%）和 施氮量<300 kg·hm−2（27%、31%） 。西北区施氮肥对棉花产量提高幅度最高（45%），但氮素吸收量提高幅度最低（35%）；华中区对棉花产量提高幅度最低（25%），但氮素吸收量提高幅度最高（65%）；华东区居中（33%、38%）。棉花产量和氮素吸收量在种植密度＜5 万株·hm−2（40%、62%）时显著高于5万～10万株·hm−2（33%、35%）和＞10万株·hm−2（18%、37%）。随着种植密度的增加，棉花产量和氮素吸收量有下降趋势。黏土施用氮肥后棉花产量和氮素吸收量提高幅度（40%、56%）显著高于砂壤土（25%、28%）和壤土（36%、39%）。在土壤有机质＞15 g·kg−1、土壤碱解氮含量50～100 g·kg−1、土壤速效磷含量15～30 mg·kg−1、土壤速效钾含量＞300 mg·kg−1情况下，施氮较不施氮棉花产量提高幅度较高。在土壤有机质10～15 g·kg−1、土壤全氮含量＞1.6 g·kg−1、土壤碱解氮含量＞100 g·kg−1、土壤速效磷含量15～30 mg·kg−1、土壤速效钾含量＜200 mg·kg−1情况下，施氮较不施氮棉花氮素吸收量幅度较高。结论 施入300～450 kg·hm−2氮肥是提高和维持棉花高产的重要措施。Abstract:Objective Yield and N-uptake of cotton plants in response to N applications and factors affecting the reaction were quantitatively analyzed with documented database.Method Meta analysis was conducted on 43 relevant articles published from 2002 to 2019 to compare the effects of the planting time, area, and N fertilization as well as the soil type, density, organic matters, total N, alkali-hydrolysable N, available phosphorus, and available potassium on the yield and N-uptake of cotton plants with those without any N application.Result The increase of 36% on yield and 51% on N-uptake of the cotton plants in 2008–2013 were significantly higher than 20% and 34%, respectively, in 2002–2007 or 28% and 33%, respectively, in 2014–2019. Increased N application promoted both yield and N-uptake of cotton plants more from the first data collection period to the second period than from the second to the last 5 years. Cotton yield and N-uptake were significantly higher with N application 300–450 kg·hm−2(38%, 58%) than N application more than 450 kg·hm−2(30%, 54%) or less than 300 kg·hm−2 (27%, 31%). There were regional variations on the effects as well. For instance, in northwest China, the fertilization resulted in the greatest increase on yield at 45% but the least on N-uptake at 35%, whereas the least on yield at 25% and the greatest on N-uptake at 65% in central China and the medium on yield at 33% and N-uptake at 38% in eastern China were observed. Increases on planting density tended to reduce the cotton yield and N-uptake. Soil type also affected the yield and N uptake as they increased significantly more in clay soil (40% and 56%, respectively) than in sandy loam (25% and 28%, respectively) and loam (36% and 39%, respectively). Grown on soil containing organic matters more than 15 g·kg−1, alkali-hydrolyzable N in the range of 50–100 g·kg−1, available phosphorus between 15–30 mg·kg−1, and available potassium greater than 300 mg·kg−1, the cotton yield was higher than without N application. And the soil containing organic matters in the range of 10–15 g·kg−1, total N more than 1.6 g·kg−1, alkali-hydrolyzable N higher than 100 g·kg−1, available phosphorus between 15–30 mg·kg−1, and available potassium less than 200 mg·kg−1 produced higher cotton N-uptake than without N fertilization.Conclusion Application of 300-450 kg·hm−2 of N fertilizer on the filed appeared to be critical for high cotton yield.
-
Keywords:
- N application /
- cotton /
- yield /
- N-uptake /
- meta-analysis
-
0. 引言
【研究意义】棉花是我国的重要经济作物，2020年棉花产量位列全球第二，为591.0万t，种植面积约占全国耕地总面积的2.64%。中国是棉花消费大国也是进口大国，每年对棉花的需求量很大。氮肥对提高农作物产量起到重要作用，其产量贡献率为30%～50%[1]。大田生产投入过多的氮肥，部分氮肥通过NH3挥发损失和反硝化成为主要的氮素损失途径，造成资源浪费、农田土壤酸化[2]和地表水体富营养化[3]等负面影响。若施入较少的氮肥则影响作物后期的养分吸收与干物质积累，较低的生物量积累使作物产量也下降。因此，明确施氮对棉花产量的影响，分析不同农艺措施与不同土壤条件下施氮的产量效应和氮素吸收量效应的变化趋势，对于提高棉花产量具有重要意义。【前人研究进展】目前，我国已进行了大量有关施氮影响棉花产量的大田试验。李鹏程等[4]的研究表明，增施氮肥360 kg·hm−2时，产量提升698.9 kg·hm−2。同时，施氮的增产效益受多方面因素的影响。李景慧等[5]在新疆库尔勒市（试验土壤有机质含量为11.89 g·kg−1）试验比较了不同施氮量对棉花的产量影响，发现施氮量在384.9～430.4 kg·hm−2时，棉花产量最高。在河南安阳（试验土壤有机质含量为11.84 g·kg−1）种植密度与施氮量互作研究表明，种植密度为5.25 万株·hm−2、施氮量为225.0 kg·hm−2时棉花产量最高[6]。薛晓萍等[7]在两种不同土壤类型下种植棉花，在黏土下施氮量为240 kg·hm−2时，氮素吸收量最高；在砂壤土下施氮量为480 kg·hm−2时，氮素吸收量最高。这说明氮素合理施用可显著提高棉花产量[8-10]，但是不同的土壤类型和不同土壤肥力下，施氮对棉花增产和氮素吸收量影响并不一致。Meta分析是整合一个主题对其独立研究进行定量结合和总结的统计方法[11]，能够有效分析不同因素对试验目标的影响。【本研究切入点】独立田间试验研究难以解读较大区域施肥对产量和氮素吸收量的影响，且在全国尺度上研究施肥对棉花产量和氮素吸收量的影响及其关键因素的研究鲜有报道。【拟解决的关键问题】本研究以Meta分析为手段，从全国尺度上解读我国棉花种植时段、种植区域、种植密度、土壤类型及其土壤理化性质等条件下施氮对棉花产量和氮素吸收量的影响，以期为棉花氮肥高效利用与管理提供借鉴。
1. 材料与方法
1.1 数据来源
本研究所采用的文献数据来源于Web of Science、中国知网（CNKI）、万方和维普数据库，文献检索所用的关键词主要包括氮（Nitrogen）、施氮量（N application）、棉花（Cotton）、产量（Yield）、氮素吸收量（N uptake）等。为了提高整合分析结果的准确性和代表性，选用的文献必须满足以下标准：（1）试验地点必须在中国大陆区域范围内；（2）试验必须为田间试验，试验材料为棉花，田间试验必须有重复；（3）试验必须包含不施氮和施氮处理；（4）文献报道中至少包含产量和氮素吸收量中的一项；（5）文献报道中应包括试验组和对照组的平均值。文献中的数据有3种呈现方式：表格、折线图和柱状图；表格中的数据直接提取，折线图和柱状图用软件Web Plot Digitizer 4.4[12]将图形数值化后再提取。按照以上标准，筛选出184组棉花产量和103组氮素吸收量相关数据（产量均值、标准差和样本量）。
1.2 Meta分析
1.2.1 响应比及其方差计算
本研究选用响应比（Response ration，RR）为统计学指标，用响应比自然对数（lnRR）来评价施氮对棉花产量和氮素吸收量的影响程度[13]，即施氮产量效应值和氮素吸收量效应。如果lnRR大于0，说明施氮对棉花产量和氮素吸收量有正效应，且lnRR值越大说明作用越明显；反之，则表示产生负效应。每项独立研究中响应比及其对数化的计算公式为：
RR=Xt/Xc (1) lnRR=ln(XtXc) (2) 式中，Xt和Xc分别表示处理组（施氮）和对照组（不施氮）的棉花产量和氮素吸收量平均值。lnRR的方差（V）计算公式为：
V=SD2tntX2t+SD2cncX2c (3) 式中，
SD2t 和SD2c 分别为处理组和对照组棉花产量和氮素吸收量的标准差，文献中SD缺失值，根据Andreas Gattingera的方法，采取相对应的试验组和对照组的1/10作为替代[14]。nt和nc分别为处理组和对照组的样本数。权重（Wij）的计算公式为：
Wij=1V (4) 加权响应比（RR++）的计算公式为：
RR++=m∑i=1ki∑j=1WijRRijm∑i=1ki∑j=1Wij (5) 式中， m是分组数（例如，不同的种植区域或种植密度等），ki是第i分组的总比较对数。
总体标准误计算公式为：
S(RR++)=√1m∑i=1ki∑j=1Wij (6) 式中，S（RR++）为总体标准误。
加权响应比（RR++）的95%置信区间计算公式为：
95%CI=RR++±1.96S(RR++) (7) 1.2.2 异质性与发表偏倚检验
本研究通过卡方检验（Chi-square test）来确定不同试验结果间是否存在异质性，并根据异质性情况来选择分析模型。若显著性检验结果为P＞0.05，说明不同试验结果间无明显差异，可选择固定效应模型（Fixed effect model，FEM）；若P＜0.05，则说明不同试验结果存在明显异质性，应选择随机效应模型（Random effect model，REM）。
发表偏倚检验采用失安全系数法（Fail-safe number，Nfs），若Nfs＞5n+10（n为数据量），则说明不存在发表偏倚，结论可信[15]。
1.2.3 加权综合响应比及其置信区间
因为纳入本研究的试验数据间存在明显的异质性，因此本研究选择随机效应模型，运用R软件中metafor包[16]来计算各研究的权重、加权综合响应比（RR++）及其95%置信区间（95%CI）。95%CI可用于RR++的假设检验：若95%CI下限＞0，说明施氮对棉花的增产效应和氮素吸收效应显著；若95%CI上限＜0，说明施氮对棉花产量和氮素吸收量存在显著的负效应；如果95%CI与0重叠，则说明施氮对棉花产量和氮素吸收无明显影响[13]。为更好地描述施氮对产量和氮素吸收量变化的影响，通过以下公式将RR++转换为变化百分率(E)：
E=(eRR++−1)×100% (8) 1.2.4 亚组分析与Meta回归分析
为明确异质性来源与大小，分析不同条件下的施氮产量效应和氮素吸收量及其影响因素，本研究根据基础产量（不施氮处理的产量）、氮素吸收量（不施氮处理的氮素吸收量）、种植区域、种植年限、种植密度、土壤条件等特征因素进行亚组分析和Meta回归分析。根据我国种植区域（华中区、西北区、华东区）、种植时段（2002–2007、2008–2013、2014–2019）、土壤类型（砂壤土、壤土、黏土）、种植密度（＜5万、5万～10万和＞10万株·hm−2）、施氮量（＜300、300～450、＞450 kg·hm−2）、土壤有机质含量（＜10、10～15、＞15 g·kg−1）、土壤全氮含量（＜1、1～1.6、＞1.6 g·kg−1）、土壤碱解氮（＜50、50～100、＞100 g·kg−1）、土壤速效磷含量（＜15、15～30、＞30 mg·kg−1）、土壤速效钾含量（＜200、200～300、＞300 mg·kg−1）来探究其对棉花产量和氮素吸收量响应氮肥的影响。
1.3 统计分析
使用Excel 2019 进行数据合并、分类，使用R语言（4.0）计算相关数据，采用GraphPad Prism 8.0 软件进行图形绘制。
2. 结果与分析
2.1 施氮对棉花产量和氮素吸收量影响的Meta分析
由表1可知，棉花产量和氮素吸收量的失安全系数为196793、121726，5n+10为930、525，失安全系数远大于5n+10，说明本研究结果受发表偏倚的影响不明显，结论具有较高的可信度。由图1可知，施氮使棉花产量平均增加了32.22%（29.30%～38.80%），氮素吸收量平均增加了43.29%（38.80%～47.78%），说明氮肥的施用在不同程度上显著提高了棉花产量和氮素吸收量。
表 1 施氮影响棉花产量和氮素吸收量的Meta分析Table 1. Meta analysis on effects of N-fertilization on yield and N-uptake of cotton plants项目
Item样本量
Sample size显著性检验
Significance test95%置信区间下限
Bottom of 95% confidence interval95%置信区间上限
Top of 95% confidence interval失安全系数
Fail-safe number5n+10 产量 Yield 184 ＜0.01 0.2930 0.3514 196793 930 氮素吸收量 N uptake 103 ＜0.01 0.3880 0.4778 121726 525 ![]()
图 1 施用氮肥对棉花产量和氮素吸收量的响应比（RR++）注：红色实线为整体效应值；红色虚线为0。Figure 1. Response ratio (RR++) of N-fertilization to yield and N-uptake of cotton plantsNote：The solid red line shows the overall response ratio; the red dotted line is 0.2.2 棉花产量和氮素吸收量的影响因素
图2所示为不同种植时段、种植地区施氮对棉花产量和氮素吸收量的效应值。由图2可以看出，种植时段对施肥增加棉花产量的提高幅度具体体现为：2002–2007年提高了20%，2008—2013年提高了36%，2014–2019年提高了28%；施肥较不施肥处理对氮素吸收量的提高幅度具体体现为：2002—2007年提高了34%，2008–2013年提高了51%，2014–2019年提高了33%。
![]()
图 2 不同种植时段、种植地区施氮对棉花产量和氮素吸收量的效应值Figure 2. Effects of N-fertilization on yield and N-uptake of cotton in different periods and regions of data collection施氮显著提高了各个区域棉花产量（图2），且在不同区域间的增产效应和氮素吸收效应存在明显差异。其中，西北区的施氮产量增幅最高为45%，但西北区的氮素吸收量增幅最低为35%；华东区增幅分别为33%和38%，华中区产量增幅较低为25%，但华中区的氮素吸收量增幅最高为67%。
图3所示为不同土壤类型、种植密度对施氮棉花产量和氮素吸收量效应值影响。由图3可知，黏土施氮增产效果最高（40%），其次为壤土（36%）和砂壤土（25%）。就氮素吸收量而言，黏土氮素吸收效应最高（56%），壤土次之（39%），砂壤土最低（28%）。
![]()
图 3 不同土壤类型、种植密度施氮对棉花产量和氮素吸收量的效应值Figure 3. Effects of N-fertilization on yield and N-uptake of cotton planted on different soil types and varied planting densities种植密度显著影响棉花产量和氮素吸收量（图3）。当种植密度小于5 万株·hm−2时，施氮的产量增幅为40%，种植密度越大，棉花产量越低。对氮素吸收量而言，种植密度小于5 万株·hm−2时提高了62%，明显高于5万～10万株·hm−2（35%）和＞10万株·hm−2（37%）的密度条件处理。
图4所示为不同施氮量、土壤有机质含量对棉花产量和氮素吸收量的效应值。由图4可知，施氮量在300～450 kg·hm−2时，施氮效果最显著，增产效应为38%，氮素吸收量提升了58%；其次为施氮量＞450 kg·hm−2时，棉花产量涨幅为30%，氮素吸收量提升了54%；最后为施氮量＜300 kg·hm−2时，棉花产量增长了27%，氮素吸收量增长了31%。
![]()
图 4 不同施氮量、土壤有机质含量施氮对棉花产量和氮素吸收量的效应值Figure 4. Effects of N-fertilization rate and soil organic matters on cotton yield and N-uptake随着土壤有机质含量的增加，施肥较不施肥处理对棉花产量提高的幅度呈增高趋势，但氮素吸收量呈先升高再降低的趋势（图4），在有机质含量为10～15 g·kg−1时，氮素吸收量增幅为48%，明显高于有机质含量为＞15 g·kg−1（40%）和＜10 g·kg−1（21%）。
图5所示为不同土壤碱解氮含量、土壤全氮含量下施氮对棉花产量和氮素吸收量的效应值。由图5可知，在土壤碱解氮含量50～100 g·kg−1时施氮增产效应最高（37%），土壤碱解氮含量＞100 g·kg−1时施氮的增产效应次之（34%），土壤碱解氮含量＜50 g·kg−1时施氮的增产效应最低（32%）。不同于产量的增加趋势，在土壤碱解氮含量＞100 g·kg−1（46%）和＜50 g·kg−1（44%）时氮素吸收效应明显高于土壤速效氮含量50～100 g·kg−1（19%）时氮素吸收效应。
![]()
图 5 不同土壤碱解氮含量、土壤全氮含量施氮对棉花产量和氮素吸收量的效应值Figure 5. Effects of N-fertilization on yield and N-uptake of cotton planted on soil with different contents of alkali hydrolyzable and total N随着土壤全氮含量的增加，施肥较不施肥处理对棉花产量提高幅度差异不大，但氮素吸收量提高幅度差异较大（图5），土壤全氮含量为＞1.6 g·kg−1、1～1.6 g·kg−1和＜1 g·kg−1时增幅为51%、28%和44%。
图6所示为不同土壤速效磷、速效钾含量对施氮棉花产量和氮素吸收量的效应值。土壤速效磷含量15～30 mg·kg−1时增产效应最高（37%），＞30 mg·kg−1的增产效应次之（36%），＜15 mg·kg−1时增产效益最低（27%）；与产量的增加趋势相同，土壤速效磷含量15～30 mg·kg−1时氮素吸收量增效最高（50%），土壤速效磷含量＞30 mg·kg−1次之（44%），＜15 mg·kg−1仅为21%。
![]()
图 6 不同土壤速效磷含量、土壤速效钾含量施氮对棉花产量和氮素吸收量的效应值Figure 6. Effects of N-fertilization on yield and N uptake of cotton planted on soil with different contents of available phosphorus and potassium土壤速效钾含量＞300 mg·kg−1时增产效应为45%，明显高于速效钾含量200～300 mg·kg−1（22%）及＜200 mg·kg−1（25%）的增产效应（图6）。随着土壤速效钾含量的升高，氮素吸收量呈先降低后增加的趋势。具体表现为：土壤速效钾含量为＜200 mg·kg−1时提高的幅度为49%，200～300 mg·kg−1提高的幅度为26%，＞300 mg·kg−1提高的幅度为39%。
2.3 棉花产量和氮素吸收量影响因素的异质性分析
表2所示为产量主要影响因素的Meta回归分析。从表2可看出，种植时间、种植区域、土壤类型、种植密度、施氮量、速效钾含量等因素引起的异质性较大，说明它们对施氮的增产效应存在明显差异且有显著影响；而土壤有机质含量、全氮含量、碱解氮含量、速效磷含量引起的异质性比较小且未达到显著水平，说明这4个因素对施氮的增产效果影响较小。
表 2 产量主要影响因素的Meta回归分析Table 2. Meta regression analysis on main factors affecting cotton yield产量影响因素
Factors affecting yield自由度
df异质性统计量
Qm显著性检验
PQM种植时段 Planting time 2 13.7413 ＜0.01 种植区域 Planting region 2 26.1813 ＜0.01 土壤类型 Soil type 2 21.3611 ＜0.01 种植密度 Planting density 2 35.0722 ＜0.01 施氮量 N application 2 11.8923 ＜0.01 有机质含量 SOM content 2 2.8682 0.2383 全氮含量 Total N content 2 0.0427 0.9789 碱解氮含量 Available N content 2 0.7656 0.6819 速效磷含量 Available P content 2 8.8255 0.0121 速效钾含量 Available K content 2 60.2263 ＜0.01 表3所示为氮素吸收量主要影响因素的Meta回归分析。从表3可看出，只有土壤全氮含量引起的异质性比较小且未达到显著水平，说明土壤全氮量对施氮效果影响较小。其余因素引起的异质性较大，说明它们对施氮的氮素吸收量效应密切相关且有显著影响。
表 3 氮素吸收量主要影响因素的Meta回归分析Table 3. Meta regression analysis on main factors affecting N-uptake of cotton plantse氮素吸收量影响因素
Factors affecting N uptake自由度
df异质性统计量
Qm显著性检验
PQM种植时段 Planting time 2 17.5224 ＜0.01 种植区域 Planting region 2 48.5255 ＜0.01 土壤类型 Soil type 2 13.0509 ＜0.01 种植密度 Planting density 2 36.2180 ＜0.01 施氮量 N application 2 33.7026 ＜0.01 有机质含量 SOM content 2 11.3464 ＜0.01 全氮含量Total N content 2 7.9870 0.0184 碱解氮含量 Available N content 2 10.1454 ＜0.01 速效磷含量 Available P content 2 9.9292 ＜0.01 速效钾含量 Available K content 2 10.4530 ＜0.01 3. 讨 论
3.1 施氮对棉花产量和氮素吸收量的整体影响分析
Meta分析结果表明，氮肥的施用可以显著提高我国棉花产量和氮素吸收量，增幅分别为32.22%和43.29%。董合林[17]研究认为棉田施用氮肥产量增长幅度为6.0%～114.7%。施肥增加作物产量的原因主要是施肥补充土壤养分的不足。棉花种植因生态区域、生态气候、栽培习惯不一致，产量差异较大[9]，施氮显著提高了各种植区域棉花产量和氮素吸收量，棉花产量和氮素吸收量表现出不同的显著性。西北地区土壤较为贫瘠，土壤氮素供应能力较弱，施氮增加氮素养分供给，增产幅度较高。任意等[18]研究也表明西北地区因水土流失、降水稀少、气候干旱、土壤水分含量低等环境因素限定了养分的可利用性，所以施用的氮肥较难吸收并转化为产量，进而产量提升幅度较大，这与本文研究结果一致。
3.2 各类影响因素分析
就种植时段而言，不同试验时间的增产效应和氮素吸收量存在显著差异，近十年由于棉花品种改良、生产技术更新，较2008年以前棉花产量和氮素吸收量增幅较大。
就种植密度而言，娄善伟等[10]研究发现种植密度影响了棉株养分累积进而影响产量的形成。氮肥的施用与密度存在密切的交互作用[6]。本研究表明，种植密度对棉花施氮产量效应和氮素吸收量有显著影响，随着密度的增加，增产效应和氮素吸收量存在下降趋势。这与广大研究者研究结果一致[19]。
就施氮量而言，赵新华等[20]提出施氮量过低以至于棉铃干物质和氮的累积也较低。薛晓萍等[21]研究表明，棉花取得了较高的生物量和氮素累积量，但因为营养体内氮素积累量较多，营养生长较旺盛，干物质与养分向生殖器官转移量少，生殖器官比例下降，由此导致经济产量降低。这与本文研究结果一致，棉花产量随着施氮量的增加先升高后降低，氮素吸收与产量趋势相同。本文研究结果为当施氮量大于450 kg·hm−2时，施氮量对棉花产量增长幅度较小。养分含量较高易使棉花营养生长和生殖生长失调，同化物更多的去向茎叶分配，导致茎叶徒长，进而降低产量[22]。
施入土壤的养分能否被植物高效吸收利用与土壤类型和土壤理化性质关系密切[23]。本研究结果表明，壤土、砂壤土和黏土施氮均显著增加了棉花产量，其黏土施肥的增产效果高于砂壤土和壤土，这与徐国春等[24]的研究结果一致。在黏土条件下黏土的氮素吸收量的响应最高，这可能是由于黏土具有较好的保肥能力[25]，且颗粒较小、养分吸附位点多[26]等原因。本研究结果表明在砂壤土条件下施氮对棉花产量和氮素吸收量的效应都较小，主要是因为砂壤土养分含量少，养分固定位点较少，施入的养分易随着水分淋失[27]。
土壤有机质作为衡量土壤肥力的重要指标，有机质含有土壤中80%以上的氮素含量，其影响氮素的供求关系[28]。本研究结果表明，有机质含量越高，施氮对棉花产量的效应越大。所以在有机质含量较低的土壤条件下，尽可能补充更多的外源氮素才能满足棉花高产的需求。刘德平等[29]研究指出氮磷有明显的互作关系，适宜的氮磷含量更好地提升作物产量。本研究结果表明土壤速效磷含量为15～30 mg·kg−1时，棉花产量最高。
3.3 Meta不确定性分析
因为Meta分析对数据要求较高，不足以纳入所有的相关文献，使本研究数据库的规模受到限制[25]。本研究没有考虑水肥一体化、测土配方等特殊施肥技术，也没有考虑到土壤质地、施肥方式、年降雨量和积温等对棉花产量和氮素吸收量的影响，没有涉及水分利用率，也未涉及水分管理模式对棉花产量和氮素吸收量的影响，同时由于不同文献中存在不同差异，这些都增加分析的不确定性。尽管如此，本研究结果依然在一定程度上阐明了施氮对我国棉花产量和氮素吸收量的影响，可以为棉花施氮效果提高和氮素营养综合管理提供借鉴。
4. 结论
利用Meta分析定量地研究了施氮肥对中国棉花产量和氮素吸收量的影响，结果显示施入氮肥对提高我国棉花产量和氮素吸收量有重要作用，施氮的平均增产效应为32.22%，氮素吸收量效应值平均增幅为43.29%。施氮肥对棉花产量的增产效应与种植时段、种植区域、土壤类型、种植密度、施氮量和土壤速效钾含量密切相关；施氮肥对棉花氮素吸收量的增长效应与种植时段、种植区域、土壤类型、种植密度、施氮量、土壤有机质含量、土壤碱解氮含量、土壤速效磷含量和土壤速效钾含量密切相关。
-
![]()
图 1 施用氮肥对棉花产量和氮素吸收量的响应比（RR++）
注：红色实线为整体效应值；红色虚线为0。
Figure 1. Response ratio (RR++) of N-fertilization to yield and N-uptake of cotton plants
Note：The solid red line shows the overall response ratio; the red dotted line is 0.
![]()
图 2 不同种植时段、种植地区施氮对棉花产量和氮素吸收量的效应值
Figure 2. Effects of N-fertilization on yield and N-uptake of cotton in different periods and regions of data collection
![]()
图 3 不同土壤类型、种植密度施氮对棉花产量和氮素吸收量的效应值
Figure 3. Effects of N-fertilization on yield and N-uptake of cotton planted on different soil types and varied planting densities
![]()
图 4 不同施氮量、土壤有机质含量施氮对棉花产量和氮素吸收量的效应值
Figure 4. Effects of N-fertilization rate and soil organic matters on cotton yield and N-uptake
![]()
图 5 不同土壤碱解氮含量、土壤全氮含量施氮对棉花产量和氮素吸收量的效应值
Figure 5. Effects of N-fertilization on yield and N-uptake of cotton planted on soil with different contents of alkali hydrolyzable and total N
![]()
图 6 不同土壤速效磷含量、土壤速效钾含量施氮对棉花产量和氮素吸收量的效应值
Figure 6. Effects of N-fertilization on yield and N uptake of cotton planted on soil with different contents of available phosphorus and potassium
表 1 施氮影响棉花产量和氮素吸收量的Meta分析
Table 1 Meta analysis on effects of N-fertilization on yield and N-uptake of cotton plants
项目
Item样本量
Sample size显著性检验
Significance test95%置信区间下限
Bottom of 95% confidence interval95%置信区间上限
Top of 95% confidence interval失安全系数
Fail-safe number5n+10 产量 Yield 184 ＜0.01 0.2930 0.3514 196793 930 氮素吸收量 N uptake 103 ＜0.01 0.3880 0.4778 121726 525 
下载: 导出CSV
表 2 产量主要影响因素的Meta回归分析
Table 2 Meta regression analysis on main factors affecting cotton yield
产量影响因素
Factors affecting yield自由度
df异质性统计量
Qm显著性检验
PQM种植时段 Planting time 2 13.7413 ＜0.01 种植区域 Planting region 2 26.1813 ＜0.01 土壤类型 Soil type 2 21.3611 ＜0.01 种植密度 Planting density 2 35.0722 ＜0.01 施氮量 N application 2 11.8923 ＜0.01 有机质含量 SOM content 2 2.8682 0.2383 全氮含量 Total N content 2 0.0427 0.9789 碱解氮含量 Available N content 2 0.7656 0.6819 速效磷含量 Available P content 2 8.8255 0.0121 速效钾含量 Available K content 2 60.2263 ＜0.01
下载: 导出CSV
表 3 氮素吸收量主要影响因素的Meta回归分析
Table 3 Meta regression analysis on main factors affecting N-uptake of cotton plantse
氮素吸收量影响因素
Factors affecting N uptake自由度
df异质性统计量
Qm显著性检验
PQM种植时段 Planting time 2 17.5224 ＜0.01 种植区域 Planting region 2 48.5255 ＜0.01 土壤类型 Soil type 2 13.0509 ＜0.01 种植密度 Planting density 2 36.2180 ＜0.01 施氮量 N application 2 33.7026 ＜0.01 有机质含量 SOM content 2 11.3464 ＜0.01 全氮含量Total N content 2 7.9870 0.0184 碱解氮含量 Available N content 2 10.1454 ＜0.01 速效磷含量 Available P content 2 9.9292 ＜0.01 速效钾含量 Available K content 2 10.4530 ＜0.01
下载: 导出CSV
-
[1] ERISMAN J W, SUTTON M A, GALLOWAY J, et al. How a century of ammonia synthesis changed the world [J]. Nature Geoscience, 2008, 1(10): 636−639. DOI: 10.1038/ngeo325
[2] GUO J H, LIU X J, ZHANG Y, et al. Significant acidification in major Chinese croplands [J]. Science, 2010, 327(5968): 1008−1010. DOI: 10.1126/science.1182570
[3] CONLEY D J, PAERL H W, HOWARTH R W, et al. Controlling eutrophication: Nitrogen and phosphorus [J]. Science, 2009, 323(5917): 1014−1015. DOI: 10.1126/science.1167755
[4] 李鹏程, 董合林, 刘爱忠, 等. 施氮量对棉花功能叶片生理特性、氮素利用效率及产量的影响 [J]. 植物营养与肥料学报, 2015, 21(1):81−91. DOI: 10.11674/zwyf.2015.0109 LI P C, DONG H L, LIU A Z, et al. Effects of nitrogen application rates on physiological characteristics of functional leaves, nitrogen use efficiency and yield of cotton [J]. Journal of Plant Nutrition and Fertilizer, 2015, 21(1): 81−91.(in Chinese) DOI: 10.11674/zwyf.2015.0109
[5] 李景慧, 韩焕勇, 王友华, 等. 氮素水平对高产杂交棉功能叶生理特性及产量的影响 [J]. 江苏农业科学, 2012, 40(4):80−83. DOI: 10.3969/j.issn.1002-1302.2012.04.025 LI J H, HAN H Y, WANG Y H, et al. Effects of nitrogen level on physiological characteristics of functional leaves and yield of high-yield hybrid cotton [J]. Jiangsu Agricultural Sciences, 2012, 40(4): 80−83.(in Chinese) DOI: 10.3969/j.issn.1002-1302.2012.04.025
[6] 李鹏程, 董合林, 刘爱忠, 等. 种植密度氮肥互作对棉花产量及氮素利用效率的影响 [J]. 农业工程学报, 2015, 31(23):122−130. DOI: 10.11975/j.issn.1002-6819.2015.23.016 LI P C, DONG H L, LIU A Z, et al. Effects of planting density and nitrogen fertilizer interaction on yield and nitrogen use efficiency of cotton [J]. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(23): 122−130.(in Chinese) DOI: 10.11975/j.issn.1002-6819.2015.23.016
[7] 薛晓萍, 沙奕卓, 郭文琦, 等. 棉花蕾花铃生物量、氮累积特征及临界氮浓度稀释模型 [J]. 生态学报, 2008, 28(12):6204−6211. DOI: 10.3321/j.issn:1000-0933.2008.12.051 XUE X P, SHA Y Z, GUO W Q, et al. Accumulation characteristics of biomass and nitrogen and critical nitrogen concentration dilution model of cotton reproductive organ [J]. Acta Ecologica Sinica, 2008, 28(12): 6204−6211.(in Chinese) DOI: 10.3321/j.issn:1000-0933.2008.12.051
[8] 凌启鸿. 作物群体质量[M]. 上海: 上海科学技术出版社, 2000. [9] WATT M S, CLINTON P W, WHITEHEAD D, et al. Above-ground biomass accumulation and nitrogen fixation of broom (Cytisus scoparius L. ) growing with juvenile Pinus radiata on a dryland site [J]. Forest Ecology and Management, 2003, 184(1/2/3): 93−104.
[10] 娄善伟, 高云光, 郭仁松, 等. 不同栽培密度对棉花植株养分特征及产量的影响 [J]. 植物营养与肥料学报, 2010, 16(4):953−958. DOI: 10.11674/zwyf.2010.0426 LOU S W, GAO Y G, GUO R S, et al. Effects of planting density on nutrition characteristics and yield of cotton [J]. Plant Nutrition and Fertilizer Science, 2010, 16(4): 953−958.(in Chinese) DOI: 10.11674/zwyf.2010.0426
[11] 李长卓. 基于Meta分析研究不同耕作方式对麦田N2O排放的影响[D]. 泰安: 山东农业大学, 2020. LI C Z. Effect of different tillage methods on N2O emissions from wheat fields based on meta-analysis[D]. Taian: Shandong Agricultural University, 2020. (in Chinese)
[12] BURDA B U, O'CONNOR E A, WEBBER E M, et al. Estimating data from figures with a Web-based program: Considerations for a systematic review [J]. Research Synthesis Methods, 2017, 8(3): 258−262. DOI: 10.1002/jrsm.1232
[13] HEDGES L V, GUREVITCH J, CURTIS P S. The meta-analysis of response ratios in experimental ecology [J]. Ecology, 1999, 80(4): 1150−1156. DOI: 10.1890/0012-9658(1999)080[1150:TMAORR]2.0.CO;2
[14] GATTINGER A, MULLER A, HAENI M, et al. Enhanced top soil carbon stocks under organic farming [J]. PNAS, 2012, 109(44): 18226−18231. DOI: 10.1073/pnas.1209429109
[15] ROSENTHAL R. The file drawer problem and tolerance for null results [J]. Psychological Bulletin, 1979, 86(3): 638−641. DOI: 10.1037/0033-2909.86.3.638
[16] VIECHTBAUER W. Conducting meta-analyses in R with the metaphor Package [J]. Journal of Statistical Software, 2010, 36(3): 1−48.
[17] 董合林. 我国棉花施肥研究进展 [J]. 棉花学报, 2007, 19(5):378−384. DOI: 10.3969/j.issn.1002-7807.2007.05.008 DONG H L. Research progress on fertilization technology of cotton [J]. Cotton Science, 2007, 19(5): 378−384.(in Chinese) DOI: 10.3969/j.issn.1002-7807.2007.05.008
[18] 任意, 张淑香, 穆兰, 等. 我国不同地区土壤养分的差异及变化趋势 [J]. 中国土壤与肥料, 2009(6):13−17. DOI: 10.3969/j.issn.1673-6257.2009.06.003 REN Y, ZHANG S X, MU L, et al. Change and difference of soil nutrients for various regions in China [J]. Soil and Fertilizer Sciences in China, 2009(6): 13−17.(in Chinese) DOI: 10.3969/j.issn.1673-6257.2009.06.003
[19] 赖奕英, 郭承君, 占东霞, 等. 不同种植密度对新疆棉花产量及纤维品质的影响 [J]. 中国棉花, 2019, 46(9):16−18. LAI Y Y, GUO C J, ZHAN D X, et al. Effect of different planting density on yield and fiber quality of cotton in Xinjiang [J]. China Cotton, 2019, 46(9): 16−18.(in Chinese)
[20] 赵新华, 束红梅, 王友华, 等. 施氮量对棉铃干物质和氮累积及分配的影响 [J]. 植物营养与肥料学报, 2011, 17(4):888−897. DOI: 10.11674/zwyf.2011.0422 ZHAO X H, SHU H M, WANG Y H, et al. Effects of nitrogen fertilization on accumulation and distribution of dry weight and nitrogen of cotton bolls [J]. Plant Nutrition and Fertilizer Science, 2011, 17(4): 888−897.(in Chinese) DOI: 10.11674/zwyf.2011.0422
[21] 薛晓萍, 王建国, 郭文琦, 等. 氮素水平对初花后棉株生物量、氮素累积特征及氮素利用率动态变化的影响 [J]. 生态学报, 2006, 26(11):3631−3640. DOI: 10.3321/j.issn:1000-0933.2006.11.015 XUE X P, WANG J G, GUO W Q, et al. Effect of nitrogen applied levels on the dynamics of biomass, nitrogen accumulation and nitrogen fertilization recovery rate of cotton after initial flowering [J]. Acta Ecologica Sinica, 2006, 26(11): 3631−3640.(in Chinese) DOI: 10.3321/j.issn:1000-0933.2006.11.015
[22] 郭金强, 危常州, 侯振安, 等. 施氮量对膜下滴灌棉花氮素吸收、积累及其产量的影响 [J]. 干旱区资源与环境, 2008, 22(9):139−142. DOI: 10.3969/j.issn.1003-7578.2008.09.027 GUO J Q, WEI C Z, HOU Z N, et al. Effect of N rates on N uptake, accumulation and yield of cotton under drip irrigation and mulch [J]. Journal of Arid Land Resources and Environment, 2008, 22(9): 139−142.(in Chinese) DOI: 10.3969/j.issn.1003-7578.2008.09.027
[23] 冯媛媛, 申艳, 徐明岗, 等. 施磷量与小麦产量的关系及其对土壤、气候因素的响应 [J]. 植物营养与肥料学报, 2019, 25(4):683−691. DOI: 10.11674/zwyf.18171 FENG Y Y, SHEN Y, XU M G, et al. Relationship between phosphorus application amount and grain yield of wheat and its response to soil and climate factors [J]. Journal of Plant Nutrition and Fertilizers, 2019, 25(4): 683−691.(in Chinese) DOI: 10.11674/zwyf.18171
[24] 许国春, 纪荣昌, 邱永祥, 等. 我国马铃薯产量对施氮的响应及其影响因素分析 [J]. 植物营养与肥料学报, 2020, 26(4):727−737. DOI: 10.11674/zwyf.19278 XU G C, JI R C, QIU Y X, et al. Responses of potato yields to nitrogen application and associated driving factors in China [J]. Journal of Plant Nutrition and Fertilizers, 2020, 26(4): 727−737.(in Chinese) DOI: 10.11674/zwyf.19278
[25] 关连珠. 普通土壤学[M]. 2版. 北京: 中国农业大学出版社, 2016. [26] 葛楠楠, 石芸, 杨宪龙, 等. 黄土高原不同土壤质地农田土壤碳、氮、磷及团聚体分布特征 [J]. 应用生态学报, 2017, 28(5):1626−1632. GE N N, SHI Y, YANG X L, et al. Distribution of soil organic carbon, total nitrogen, total phosphorus and water stable aggregates of cropland with different soil textures on the Loess Plateau, Northwest China [J]. Chinese Journal of Applied Ecology, 2017, 28(5): 1626−1632.(in Chinese)
[27] LIN Y R, WATTS D B, VAN SANTEN E, et al. Influence of poultry litter on crop productivity under different field conditions: A meta-analysis [J]. Agronomy Journal, 2018, 110(3): 807−818. DOI: 10.2134/agronj2017.09.0513
[28] 朱洪芬, 南锋, 徐占军, 等. 黄土高原盆地土壤有机质与影响因子的空间多尺度关系 [J]. 生态学报, 2017, 37(24):8348−8360. ZHU H F, NAN F, XU Z J, et al. Multi-scale spatial relationships between soil organic matter and influencing factors in basins of the Chinese Loess Plateau [J]. Acta Ecologica Sinica, 2017, 37(24): 8348−8360.(in Chinese)
[29] 刘德平, 杨树青, 史海滨, 等. 氮磷配施条件下作物产量及水肥利用效率 [J]. 生态学杂志, 2014, 33(4):902−909. LIU D P, YANG S Q, SHI H B, et al. Crop yield and water-fertilizer utilization efficiency under combined application of nitrogen and phosphorous [J]. Chinese Journal of Ecology, 2014, 33(4): 902−909.(in Chinese)
-
期刊类型引用(3)
1. 王嘉琳，张世柏，王涛，赵竹青，石磊. 我国油菜、大豆和烤烟硼肥增产提质的整合分析. 中国土壤与肥料. 2024(12): 176-184 . 
百度学术
2. 朱盼盼，石磊. 我国主要粮食作物施用锌肥增产效果的整合分析. 中国土壤与肥料. 2023(06): 213-218 . 百度学术
3. 程艳莉，张芬，刘发波，方林发，孙铭，梁涛，刘敦一，陈新平，王孝忠. 滴灌施肥对不同类型蔬菜和果树产量、水氮利用效率和品质的影响. 植物营养与肥料学报. 2023(09): 1677-1688 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 531
- HTML全文浏览量: 405
- PDF下载量: 28
- 被引次数: 6
