基于Meta分析中國(guó)水稻產(chǎn)量對(duì)施肥的響應(yīng)特征

韓天富，馬常寶，黃晶,3，柳開樓,4，薛彥東，李冬初,3，劉立生,3，張璐,3，劉淑軍,3，張會(huì)民,3

?

基于Meta分析中國(guó)水稻產(chǎn)量對(duì)施肥的響應(yīng)特征

韓天富1，馬常寶2，黃晶1,3，柳開樓1,4，薛彥東2，李冬初1,3，劉立生1,3，張璐1,3，劉淑軍1,3，張會(huì)民1,3

（1中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所/耕地培育技術(shù)國(guó)家工程實(shí)驗(yàn)室，北京 100081；2農(nóng)業(yè)農(nóng)村部耕地質(zhì)量監(jiān)測(cè)保護(hù)中心，北京 100125；3中國(guó)農(nóng)業(yè)科學(xué)院衡陽紅壤實(shí)驗(yàn)站/祁陽農(nóng)田生態(tài)系統(tǒng)國(guó)家野外試驗(yàn)站，湖南祁陽 426182；4江西省紅壤研究所/國(guó)家紅壤改良工程技術(shù)研究中心，南昌 330046）

【】定量分析近30年施肥對(duì)中國(guó)水稻產(chǎn)量的綜合效應(yīng)和影響機(jī)制，為水稻種植區(qū)域肥料的科學(xué)施用提供依據(jù)。以全國(guó)水稻土長(zhǎng)期監(jiān)測(cè)點(diǎn)為平臺(tái)，將相應(yīng)的監(jiān)測(cè)數(shù)據(jù)按照種植區(qū)域、試驗(yàn)時(shí)間、種植制度、作物類型、施肥類型、土壤質(zhì)地、土壤pH、土壤有機(jī)質(zhì)含量、土壤全氮含量、土壤有效磷含量、土壤速效鉀含量、土壤緩效鉀含量進(jìn)行分組，以不施肥處理作為對(duì)照，利用Meta-analysis方法探究施肥對(duì)水稻產(chǎn)量的綜合效應(yīng)及其影響因素。近10年（2008—2017）以來，無論施肥與否，水稻產(chǎn)量均顯著高于1988—1997和1998—2007年對(duì)應(yīng)的水稻產(chǎn)量。與不施肥相比，施肥顯著提高水稻產(chǎn)量，其提高幅度平均為80.8%。在西南地區(qū)施肥對(duì)水稻產(chǎn)量的提高幅度最高（98.5%），顯著高于華北地區(qū)（70.3%）。不同試驗(yàn)時(shí)間下，施肥比不施肥處理在1988—1997年對(duì)水稻產(chǎn)量提高的幅度（99.1%）高于1998—2007年（84.2%）和 2008-2017年（78.1%）。不同種植制度下，施肥較不施肥處理能顯著提高一年三熟水稻產(chǎn)量（92.0%），且提高幅度均高于一年一熟（76.2%）和一年兩熟（81.9%）。與不施肥相比，雙季稻施肥對(duì)水稻產(chǎn)量的提高幅度（85.9%）高于單季稻區(qū)（75.9%）和水稻-其他作物（79.5%）。與不施肥相比，有機(jī)肥與無機(jī)肥配合施用對(duì)水稻產(chǎn)量提高幅度（88.3%）高于化肥單施處理（76.6%）。施肥較不施肥處理能顯著提高黏質(zhì)土壤水稻產(chǎn)量（92.0%），提高幅度顯著高于砂質(zhì)土壤（58.0%）和壤質(zhì)土壤（77.5%）。隨著土壤有機(jī)質(zhì)和有效磷含量的增加，施肥較不施肥處理水稻產(chǎn)量提高的幅度呈降低趨勢(shì)。在較高的土壤pH（＞7.5）、較低土壤全氮（＜1.5 g·kg-1）和緩效鉀（＜150 mg·kg-1）情況下，施肥較不施肥處理水稻產(chǎn)量提高的幅度較高。隨機(jī)森林分析結(jié)果表明：施肥對(duì)水稻產(chǎn)量提高幅度主要受水稻種植區(qū)域、土壤全氮和種植制度的影響。此外，肥料的農(nóng)學(xué)效率與施肥對(duì)水稻產(chǎn)量增產(chǎn)幅度呈極顯著正相關(guān)。雖然當(dāng)前施肥對(duì)水稻產(chǎn)量增加的趨勢(shì)在降低，但是適量的肥料投入（尤其是西南地區(qū)）是提高和維持水稻高產(chǎn)的重要措施，尤其是有機(jī)肥與無機(jī)肥配合施用增產(chǎn)效果更加顯著。同時(shí)，在種植制度的基礎(chǔ)上，各水稻種植區(qū)域應(yīng)結(jié)合土壤質(zhì)地、土壤氮素和鉀素等方面作為肥料投入的主要依據(jù)。

施肥；水稻產(chǎn)量；Meta-analysis；響應(yīng)比；農(nóng)學(xué)效率

0 引言

【研究意義】全球超過60%的人口以水稻 (L.) 為主食[1]，預(yù)計(jì)到2035年大米需求量將從2010年的6.76億噸增加到8.52億噸[2]。中國(guó)的水稻種植面積和單位產(chǎn)量均位居世界各國(guó)之首，對(duì)世界的糧食安全做出了重要貢獻(xiàn)[3]。肥料的投入是保證作物高產(chǎn)穩(wěn)產(chǎn)的重要手段之一[4]，不合理的肥料投入不僅抑制作物正常的生長(zhǎng)發(fā)育[5]，同時(shí)還可能造成土壤板結(jié)、土壤酸化、肥料利用率下降、土壤酶活性降低等一系列負(fù)面效應(yīng)，最終影響作物產(chǎn)量的提高[6]。因此，研究施肥對(duì)水稻產(chǎn)量的影響及其關(guān)鍵作用因子對(duì)于水稻高產(chǎn)穩(wěn)產(chǎn)、土壤培肥等具有重要意義。【前人研究進(jìn)展】盡管國(guó)內(nèi)外在稻田培肥方面進(jìn)行了大量的研究，但大部分主要集中在施肥對(duì)水稻產(chǎn)量或肥料利用率等方面的影響[7-9]。比如，苑俊麗等[10]運(yùn)用整合分析方法研究了高效氮肥較常規(guī)化肥施用對(duì)中國(guó)水稻產(chǎn)量和氮素吸收量的影響，很少從不同種植區(qū)域、管理措施、土壤理化性質(zhì)等方面探討水稻產(chǎn)量對(duì)施肥響應(yīng)的差異特征及影響因素。首先，土壤基礎(chǔ)肥力水平是決定作物能否高效利用肥料的關(guān)鍵因子[11]。其次，在品種和其他管理措施相對(duì)穩(wěn)定的情況下，施肥處理的產(chǎn)量高低主要決定于施肥處理本身與環(huán)境互作效應(yīng)[12]。然而，由于我國(guó)水稻種植區(qū)域遼闊，土壤肥力差異較大，導(dǎo)致水稻產(chǎn)量對(duì)施肥的響應(yīng)特征各不相同[13-14]。方暢宇等[15]研究表明，基礎(chǔ)地力較低的土壤上優(yōu)先施用化肥，輔助施用有機(jī)肥；肥力較高的土壤上輕施化肥，多施有機(jī)肥以達(dá)到水稻高產(chǎn)穩(wěn)產(chǎn)的目的。因此，探明不同稻田土壤肥力水平下長(zhǎng)期施肥對(duì)水稻產(chǎn)量的影響，進(jìn)而為各水稻種植區(qū)域不同水稻種植制度下合理施肥提供依據(jù)顯得尤為重要。【本研究切入點(diǎn)】首先，獨(dú)立田間試驗(yàn)研究難以準(zhǔn)確回答較大區(qū)域上施肥對(duì)產(chǎn)量的影響，且在全國(guó)尺度上探究施肥對(duì)水稻產(chǎn)量的影響及其關(guān)鍵因素的研究到目前為止還鮮見報(bào)道。農(nóng)業(yè)農(nóng)村部在1988—2017年期間先后設(shè)置了一批田間施肥試驗(yàn)（常規(guī)施肥和不施肥處理），涵蓋了中國(guó)水稻主要的種植區(qū)域，主要記錄水稻產(chǎn)量、土壤肥力、管理措施等數(shù)據(jù)，結(jié)合Meta分析方法——對(duì)同一主題下多個(gè)研究結(jié)果進(jìn)行綜合定量分析的方法[16]，可從全國(guó)尺度回答近30年不同水稻種植區(qū)域、管理措施、土壤理化性質(zhì)等條件下施肥對(duì)水稻產(chǎn)量的影響。其次，基于隨機(jī)森林算法——一種高效的組合分類方法，根據(jù)變量重要性程度進(jìn)而提取特征變量[17]，進(jìn)而探究施肥對(duì)水稻產(chǎn)量影響的關(guān)鍵因素。【擬解決的關(guān)鍵問題】本研究運(yùn)用Meta分析方法，以農(nóng)業(yè)農(nóng)村部設(shè)置的不施肥處理為對(duì)照，從全國(guó)尺度上分析近30年水稻主產(chǎn)區(qū)的水稻產(chǎn)量對(duì)施肥的響應(yīng)特征及其關(guān)鍵的影響因素，旨在為肥料的合理施用并實(shí)現(xiàn)水稻的高產(chǎn)穩(wěn)產(chǎn)提供理論依據(jù)。

1 材料與方法

1.1 數(shù)據(jù)來源

本研究所用的數(shù)據(jù)均來源于1988—2017年“農(nóng)業(yè)農(nóng)村部耕地土壤質(zhì)量監(jiān)測(cè)”工作開展期間所收集的水稻田間試驗(yàn)數(shù)據(jù)。水稻土監(jiān)測(cè)點(diǎn)的布局主要分為以下幾個(gè)部分：東北地區(qū)（黑龍江、吉林和遼寧）、華北地區(qū)（河南、山東和陜西）、西南地區(qū)（云南、四川、重慶和貴州）、長(zhǎng)江中游（湖北、湖南和江西）、長(zhǎng)江下游（安徽、江蘇、上海和浙江）、華南地區(qū)（福建、廣東、廣西和海南）。監(jiān)測(cè)地塊的地理位置、耕作制度、土壤類型、作物類型、分布面積、管理水平等在各區(qū)域均有較好的代表性。各監(jiān)測(cè)點(diǎn)試驗(yàn)設(shè)置不施肥區(qū)（空白區(qū)）和常規(guī)施肥區(qū)（農(nóng)民習(xí)慣施肥管理）兩個(gè)處理，進(jìn)而探究與不施肥相比，常規(guī)施肥處理對(duì)水稻產(chǎn)量差異的影響。不施肥處理中，小區(qū)面積為32—67 m2，用水泥板或者其他材料做擋板，防止水肥橫向的轉(zhuǎn)移；施肥處理按照當(dāng)?shù)剞r(nóng)田常年肥料用量均勻撒施，小區(qū)面積不低于334 m2，鑒于小區(qū)面積較大能夠基本反映施肥和不施肥的差異，再加上經(jīng)費(fèi)等方面的限制，因此，所有處理均為1次重復(fù)。不施肥處理和施肥處理除了施肥用量不一致以外，其他措施均相同。監(jiān)測(cè)內(nèi)容主要包括：作物產(chǎn)量、施肥用量和土壤理化性狀等。作物產(chǎn)量包括每一季作物的實(shí)際產(chǎn)量，采用去邊行后實(shí)打、實(shí)收的方法測(cè)定；施肥情況主要包括每一季作物有機(jī)肥和化肥的施用日期、肥料品種、施肥次數(shù)和施肥用量等，秸稈還田按照實(shí)際用量以有機(jī)肥形式記錄；土壤理化性狀主要包括土壤有機(jī)質(zhì)、全氮、有效磷、速效鉀和緩效鉀、pH、質(zhì)地等，在每年最后一季作物收獲以后，按照“隨機(jī)”“等量”“多點(diǎn)混合”的原則，立即采集耕層土壤樣品進(jìn)行上述指標(biāo)的測(cè)定，按照《土壤分析技術(shù)規(guī)范》[18]進(jìn)行樣品的測(cè)定。但是由于部分點(diǎn)位管理不是特別完善，導(dǎo)致某些年份水稻產(chǎn)量數(shù)據(jù)出現(xiàn)缺失，為了彌補(bǔ)某些年份產(chǎn)量數(shù)據(jù)的缺失和各點(diǎn)位試驗(yàn)設(shè)置缺少重復(fù)，且綜合考慮監(jiān)測(cè)時(shí)間的跨度后將所有數(shù)據(jù)分為：1988—1997、1998—2007和2008—2017年3個(gè)時(shí)間段，將各時(shí)間段內(nèi)的產(chǎn)量數(shù)據(jù)作為重復(fù)進(jìn)行下一步的數(shù)據(jù)分析。通過對(duì)所有數(shù)據(jù)的篩選計(jì)算，最終得到462組數(shù)據(jù)。

1.2 研究方法

本研究中的數(shù)據(jù)均來自農(nóng)業(yè)農(nóng)村部設(shè)置的定位監(jiān)測(cè)數(shù)據(jù)。統(tǒng)計(jì)學(xué)指標(biāo)采用權(quán)重響應(yīng)比（response ratios，）表示，并計(jì)算其95%的置信區(qū)間（95%）。其計(jì)算公式為：

1.3 數(shù)據(jù)分析

首先，通過卡方檢驗(yàn)（Chi-square test）明確試驗(yàn)處理之間及各試驗(yàn)結(jié)果是否存在異質(zhì)性（處理間或不同研究結(jié)果間的變異是否由隨機(jī)誤差引起）。若納入的各研究結(jié)果無異質(zhì)性（＞0.05），采用固定效應(yīng)模型進(jìn)行分析（fixed effect model，F(xiàn)EM），相反，則采用隨機(jī)效應(yīng)模型（random effect model，REM）[23]。表1為本研究數(shù)據(jù)的異質(zhì)性檢驗(yàn)結(jié)果，故采用REM進(jìn)行數(shù)據(jù)分析。其次，采用Meta Win 2.1軟件進(jìn)行分析合并計(jì)數(shù)資料的響應(yīng)比得出加權(quán)平均響應(yīng)。

表1 樣本量描述性統(tǒng)計(jì)

Q 為異質(zhì)性檢驗(yàn)的統(tǒng)計(jì)量；Q 為Q統(tǒng)計(jì)量顯著性

Q is the statistic of heterogeneity;Q is the significant value of Q

本研究綜合考慮了水稻種植區(qū)域（東北地區(qū)、華北地區(qū)、西南地區(qū)、長(zhǎng)江中游、長(zhǎng)江下游、華南地區(qū)）、試驗(yàn)時(shí)間（1988—1997、1998—2007、2008—2017年）、種植制度（一年一熟、一年兩熟、一年三熟）、作物類型（單季稻、雙季稻、水稻-其他作物）、施肥類型（化肥單施、化肥+有機(jī)肥）、土壤質(zhì)地（黏土、壤土、砂土）、土壤pH（＜6.5、6.5—7.5、＞7.5）、土壤有機(jī)質(zhì)含量（＜20、20—30、＞30 g·kg-1）、土壤全氮含量（＜1.5、1.5—2、＞2 g·kg-1）、土壤有效磷含量（＜10、10—20、＞20 mg·kg-1）、土壤速效鉀含量（＜50、50—100、＞100 mg·kg-1）、土壤緩效鉀含量（＜150、150—300、＞300 mg·kg-1）對(duì)水稻產(chǎn)量響應(yīng)施肥的影響，其中，種植區(qū)域、種植制度、作物類型、施肥類型沒有時(shí)間上的差異，而土壤理化性質(zhì)是指3個(gè)時(shí)間階段相對(duì)應(yīng)的初始年份土壤理化性質(zhì)，根據(jù)已發(fā)表的文獻(xiàn)[24]和全國(guó)第二次土壤普查得到的土壤養(yǎng)分含量分級(jí)標(biāo)準(zhǔn)，并結(jié)合本研究數(shù)據(jù)整體分布范圍進(jìn)行相應(yīng)的等級(jí)劃分。采用單因素方差分析和最小顯著性差（LSD）評(píng)價(jià)不同區(qū)域、土地利用類型和土壤理化性質(zhì)下水稻產(chǎn)量數(shù)據(jù)異質(zhì)性和差異性[25]。

肥料農(nóng)學(xué)效率[26]的計(jì)算公式如下：

肥料農(nóng)學(xué)效率=（施肥處理水稻產(chǎn)量-不施肥處理水稻產(chǎn)量）/肥料總施用量 （7）

數(shù)據(jù)處理和部分圖形制作采用Excel 2003；利用R語言（3.4.4）進(jìn)行數(shù)據(jù)的隨機(jī)森林重要性分析；采用SPSS軟件11.0（SPSS Inc.，Chicago，IL，USA）進(jìn)行數(shù)據(jù)統(tǒng)計(jì)分析和顯著性檢驗(yàn)；利用Sigmaplot 10.0軟件進(jìn)行圖形的制作。

2 結(jié)果

2.1 不同施肥處理下水稻產(chǎn)量的差異及其響應(yīng)比的分布

近30年以來，施肥處理較不施肥處理能顯著增加水稻產(chǎn)量（圖1-A），兩者的平均產(chǎn)量分別為9.7和5.2 t·hm-2。近10年（2008—2017）以來，無論施肥與否，水稻產(chǎn)量（9.9、5.4 t·hm-2）均顯著高于1988—1997和1998—2007年對(duì)應(yīng)的水稻產(chǎn)量（8.5、4.5 t·hm-2和9.5、5.1 t·hm-2）。利用Meta軟件分析了462組水稻產(chǎn)量對(duì)施肥的響應(yīng)比（本研究水稻的響應(yīng)比是指施肥對(duì)產(chǎn)量增加的幅度），且分布檢驗(yàn)表明全部響應(yīng)比符合正態(tài)分布（＜0.001），平均值為0.60±0.31（圖1-B），因此不用進(jìn)行數(shù)據(jù)轉(zhuǎn)化而直接進(jìn)行因素分析。

2.2 水稻產(chǎn)量對(duì)施肥的響應(yīng)及因素分析

Meta分析結(jié)果表明，與不施肥相比，施肥能顯著提高水稻產(chǎn)量，其提高幅度為80.8%（置信區(qū)間為76.4%—83.4%）（圖2）。不同種植區(qū)域下，施肥較不施肥處理對(duì)水稻產(chǎn)量的提高幅度各不相同，各區(qū)域增加幅度分別為：西南區(qū)為98.5%（置信區(qū)間為85.7%—112.1%）、華南區(qū)為80.0%（置信區(qū)間為69.6%—90.9%）、長(zhǎng)江下游為78.1%（置信區(qū)間為69.4%—87.3%）、東北區(qū)為78.0%（置信區(qū)間為59.1%—99.1%）、長(zhǎng)江中游為77.8%（置信區(qū)間為70.6%—85.3%）、華北區(qū)為70.3%（置信區(qū)間為45.0%—100.1%）。不同試驗(yàn)時(shí)間下，施肥較不施肥處理對(duì)水稻產(chǎn)量的提高幅度具體表現(xiàn)為：1988—1997年為99.1%（置信區(qū)間為79.4%—121.0%）、1998—2007年為84.2%（置信區(qū)間為75.4%—93.6%）和 2008—2017年為78.1%（置信區(qū)間為72.9%—83.5%）。不同種植制度下，施肥較不施肥能顯著提高一年三熟水稻產(chǎn)量（92.0%，置信區(qū)間為74.2%—111.7%），且提高幅度高于一年一熟（76.2%，置信區(qū)間為68.5%—84.2%）和一年兩熟（81.9%，置信區(qū)間為76.5%—87.4%）。

CK：無肥處理，F(xiàn)：常規(guī)施肥處理。不同的小寫字母表示差異顯著（P＜0.05）；箱形框中間的實(shí)線代表中位數(shù)，虛線代表平均值。箱形框上下邊緣分別代表全部數(shù)據(jù)的75%和25%。箱形框上下兩條線的邊緣分別代表全部數(shù)據(jù)的95%和5%，上下兩個(gè)實(shí)心點(diǎn)為異常值。括號(hào)里的數(shù)字表示該處理有水稻產(chǎn)量數(shù)據(jù)的個(gè)數(shù)

點(diǎn)和誤差線分別代表響應(yīng)比及其95%的置信區(qū)間，如果誤差線沒有跨越零線表示處理與對(duì)照存在顯著差異；括號(hào)內(nèi)的數(shù)值代表樣本數(shù)。下同

由圖3可知，施肥對(duì)水稻產(chǎn)量的提高與作物類型、施肥措施和土壤質(zhì)地密切相關(guān)。與不施肥相比，施肥在雙季稻對(duì)水稻產(chǎn)量的提高幅度（85.9%，置信區(qū)間為78.5%—93.5%）高于單季稻區(qū)（75.9%，置信區(qū)間為67.9%—84.2%）和水稻—其他作物（79.5%，置信區(qū)間為71.9%—87.5%）與不施肥相比，化肥與有機(jī)肥配施對(duì)水稻產(chǎn)量提高的幅度（88.3%，置信區(qū)間為81.1%—95.8%）略高于化肥單施（76.6%，置信區(qū)間為71.6%—81.7%）。不同土壤質(zhì)地條件下，施肥較不施肥處理在黏土上對(duì)水稻產(chǎn)量的提高幅度（96.6%，置信區(qū)間為86.1%—107.7%）顯著高于砂土（57.9%，置信區(qū)間為25.8%—98.5%），但是與在壤土上的增幅（77.5%，置信區(qū)間為72.7%—82.4%）相比無顯著差異。

施肥對(duì)產(chǎn)量的提高程度與土壤的理化性質(zhì)也密切相關(guān)（圖4和5）。就土壤pH而言，在偏中性（pH為6.5—7.5）的土壤上，施肥較不施肥處理對(duì)水稻產(chǎn)量的提高幅度為77.2%（置信區(qū)間為67.7%—88.4%），在pH＞7.5和pH＜6.5的土壤上提高幅度較高，分別為87.1%（置信區(qū)間為72.3%—103.0%）和80.7%（置信區(qū)間為75.8%—85.8%）。隨著土壤有機(jī)質(zhì)（SOM）含量的增加，施肥較不施肥處理對(duì)水稻產(chǎn)量提高的幅度呈降低趨勢(shì)，具體為：SOM＜20 g·kg-1時(shí)提高的幅度為83.1%（置信區(qū)間為70.8%—96.3%）、20 g·kg-1＜SOM＜30 g·kg-1時(shí)提高的幅度為81.6%（置信區(qū)間為74.4%—89.1%）、SOM＞30 g·kg-1時(shí)提高的幅度為79.7%（置信區(qū)間為73.7%%—85.8%）。就土壤有效磷（AP）含量而言，在AP較低（＜10 mg·kg-1）情況下，施肥較不施肥處理對(duì)水稻產(chǎn)量提高幅度（87.4%，置信區(qū)間為79.7%—95.4%）顯著高于AP較高（＞20 mg·kg-1）情況下的提高幅度（74.1%，置信區(qū)間為67.1%—81.5%）。隨著土壤全氮含量（STN）的增加，施肥較不施肥處理對(duì)水稻產(chǎn)量提高的幅度呈先降低后增加趨勢(shì)，具體表現(xiàn)為：當(dāng)STN＜1.5 g·kg-1時(shí)提高的幅度為86.8%（置信區(qū)間為78.4%95.5%）、1.5 g·kg-1＜STN＜2 g·kg-1時(shí)提高的幅度為76.2%（置信區(qū)間為69.4%—83.3%）、STN＞2 g·kg-1時(shí)提高的幅度為81.2%（置信區(qū)間為73.4%—89.4%）。施肥對(duì)水稻產(chǎn)量提高的幅度隨土壤速效鉀含量的增加呈先升高后降低趨勢(shì)。就土壤緩效鉀（SAK）含量而言，在SAK ＜150 mg·kg-1情況下，施肥較不施肥處理對(duì)水稻產(chǎn)量提高幅度（87.5%，置信區(qū)間為77.5%—98.1%）高于SAK≥150 mg·kg-1（79.0%，置信區(qū)間為71.1%—86.1%）水平下的提高幅度。

圖3 不同作物、肥料類型和土壤質(zhì)地下水稻產(chǎn)量對(duì)施肥的權(quán)重響應(yīng)比

SOM：土壤有機(jī)質(zhì)；STN：土壤全氮。下同

AP：土壤有效磷；AK：土壤速效鉀，SAK：土壤緩效鉀。下同

2.3 不同因素對(duì)水稻產(chǎn)量響應(yīng)比的重要性

利用隨機(jī)森林對(duì)水稻產(chǎn)量的影響因素進(jìn)行重要性分析，結(jié)果如圖6所示。Mean decrease accuracy是指預(yù)測(cè)誤差準(zhǔn)確性降低的程度，該值越大表示該變量的重要性越大[27]。各指標(biāo)對(duì)水稻產(chǎn)量響應(yīng)比均有一定的影響，通過比較各變量因素的重要性可知，其中種植區(qū)域（Region）、土壤全氮（STN）、種植制度（Cropping system）、土壤質(zhì)地（Soil texture）和緩效鉀（SAK）5個(gè)因素的重要程度較大，速效鉀（AK）、土壤有機(jī)質(zhì)（SOM）、pH、施肥類型（Fertilization type）、時(shí)間（Time）和土壤有效磷（AP）影響較小。

圖6 變量的重要性

2.4 水稻產(chǎn)量響應(yīng)比與肥料農(nóng)學(xué)效率的關(guān)系

進(jìn)一步分析了水稻產(chǎn)量響應(yīng)比與肥料農(nóng)學(xué)效率之間的關(guān)系（圖7），發(fā)現(xiàn)兩者呈極顯著正相關(guān)關(guān)系（＜0.001），通過線性擬合可得，每增加1個(gè)單位的肥料農(nóng)學(xué)效率，水稻產(chǎn)量的響應(yīng)比相應(yīng)地提高0.05個(gè)單位。就目前較為普遍的兩種施肥類型而言（化肥單施、化肥與有機(jī)肥配施），在增加相同單位的肥料農(nóng)學(xué)效率情況下，化肥與有機(jī)肥配施處理對(duì)水稻產(chǎn)量響應(yīng)比的增加速率高于化肥單施處理。

CF：化肥單施；CMF：化肥+有機(jī)肥

3 討論

通過施肥為作物提供足夠的養(yǎng)分是維持作物高產(chǎn)穩(wěn)產(chǎn)、保障世界糧食安全的重要途徑之一[28-29]。本研究通過整合農(nóng)業(yè)農(nóng)村部耕地質(zhì)量監(jiān)測(cè)的水稻產(chǎn)量數(shù)據(jù)發(fā)現(xiàn)，與不施肥相比，近30年以來常規(guī)施肥條件下水稻產(chǎn)量平均提高80.8%（置信區(qū)間為76.4%—83.4%）。這與CHIANU等[30]研究得到施肥能夠增加作物產(chǎn)量50%—100%的結(jié)果相符。作物生長(zhǎng)過程中需要大量的營(yíng)養(yǎng)物質(zhì)，而化肥中含有高濃度的養(yǎng)分物質(zhì)，施用以后不僅能夠供給植物吸收，還能補(bǔ)充因作物吸收而帶走或隨水分流失的養(yǎng)分[31]。不施肥處理由于土壤養(yǎng)分常年被消耗且供應(yīng)不足導(dǎo)致作物產(chǎn)量顯著低于常規(guī)施肥處理（圖1），也間接反映施肥的重要性。研究表明，通過合理的水肥養(yǎng)分管理，不僅能消除由于營(yíng)養(yǎng)過剩造成的負(fù)面環(huán)境效應(yīng)，同時(shí)可實(shí)現(xiàn)近30%的增產(chǎn)潛力[32]。近10年（2008—2017年）水稻產(chǎn)量顯著高于其他年間施肥處理，與張福鎖等[11]研究結(jié)果類似，這主要與近年來高產(chǎn)品種的大力推廣[33]和土壤肥力總體普遍提高[34]密切相關(guān)。據(jù)聯(lián)合國(guó)糧農(nóng)組織（FAO）的數(shù)據(jù)統(tǒng)計(jì)，化肥的增產(chǎn)作用占到農(nóng)作物產(chǎn)量的50%[35]，這主要是因?yàn)楸狙芯渴窃诔耸┓剩渌蛩鼐３忠恢碌那闆r下開展的田間試驗(yàn)，相對(duì)提高了施肥的作用，弱化了品種、管理措施等因素對(duì)產(chǎn)量的影響。

不同區(qū)域因氣候特征、自然環(huán)境和耕作習(xí)慣等因素不同導(dǎo)致施肥對(duì)產(chǎn)量的影響差異較大。本研究發(fā)現(xiàn)相比其他水稻種植區(qū)域，施肥在西南地區(qū)增產(chǎn)效果較好。研究表明，在四川省施用氮肥最高可增產(chǎn)90%，且相對(duì)貧瘠的山區(qū)施肥增產(chǎn)效果更加顯著[36]。本研究的西南區(qū)域包括云南、四川、重慶和貴州，屬于云貴高原區(qū)域，土壤肥力相對(duì)較低[37]，因此施肥增產(chǎn)效果高于其他地區(qū)。近30年間，施肥對(duì)產(chǎn)量的增加趨勢(shì)在下降，這主要與土壤質(zhì)量的改善和土壤肥力普遍提高有關(guān)[11]；其次是長(zhǎng)期連續(xù)種植同一作物容易產(chǎn)生“連作障礙”，如土壤養(yǎng)分異常累積，微生物種群結(jié)構(gòu)失衡等，最終影響作物產(chǎn)量的提高[38]。一年一熟（單季稻）水稻種植區(qū)域主要分布在東北地區(qū)，受氣候條件的影響，導(dǎo)致該地區(qū)水稻產(chǎn)量均較低[39]，另外就是該地區(qū)土壤肥力較高，相對(duì)弱化了施肥增產(chǎn)的效應(yīng)，施肥較不施肥處理對(duì)單季稻產(chǎn)量提高的幅度較低也印證這一結(jié)果。一年三熟地區(qū)主要分布在南方高溫多雨的紅壤地區(qū)，土壤養(yǎng)分貧瘠，酸化嚴(yán)重[40]，施肥能夠快速供應(yīng)大量的氮磷鉀養(yǎng)分，滿足作物的生長(zhǎng)需求，相比不施肥能顯著增加作物產(chǎn)量。本研究結(jié)果表明：有機(jī)無機(jī)肥配施在水稻上的增產(chǎn)效果高于單施化肥，這與國(guó)內(nèi)外大部分長(zhǎng)期定位試驗(yàn)結(jié)果相似[41-42]。有機(jī)無機(jī)肥配施有利于水稻中后期干物質(zhì)累積和養(yǎng)分吸收，提高了單位面積總穗數(shù)和穗粒數(shù)[43]。隨機(jī)森林結(jié)果表明：水稻種植區(qū)域和種植制度對(duì)水稻產(chǎn)量響應(yīng)比的影響較大，因此，針對(duì)不同的種植區(qū)域和種植制度，需要指定相應(yīng)的施肥方案，與前人提出的“大配方、小調(diào)整”區(qū)域施肥方案相符[14]。

土壤是作物生長(zhǎng)的主要載體, 也是我們當(dāng)前培肥管理的主要對(duì)象。土壤理化性質(zhì)的優(yōu)劣不僅決定作物根系生長(zhǎng)是否良好，同時(shí)與施入土壤中的養(yǎng)分能否高效利用密切相關(guān)。本研究結(jié)果表明，在砂質(zhì)土壤上，施肥對(duì)水稻增產(chǎn)的效果較差，主要是因?yàn)樯巴令w粒較大，養(yǎng)分的固定位點(diǎn)較少，導(dǎo)致施入的養(yǎng)分易隨水分淋失[44]，其次是砂質(zhì)土壤中養(yǎng)分含量少，且保水保肥能力較差[45]。土壤有機(jī)質(zhì)含量是衡量土壤綜合肥力的一個(gè)重要指標(biāo)[46]，在低有機(jī)質(zhì)含量下土壤肥力較低，施肥增產(chǎn)的效應(yīng)較高（圖4），這與曾明祥等[3]研究結(jié)果一致。本研究結(jié)果表明：隨著土壤全氮含量的增加，施肥較不施肥處理對(duì)水稻產(chǎn)量提高的幅度呈先降低后增加趨勢(shì)。原因主要為：（1）氮素對(duì)水稻生產(chǎn)的影響僅次于水分，增施氮肥能夠顯著的提高水稻產(chǎn)量[47]，且在土壤低氮水平下，增施氮肥的增產(chǎn)效果較明顯[48]；（2）隨著土壤全氮含量的增加，土壤能夠供給水稻較多的氮素，進(jìn)而弱化施肥增產(chǎn)的效應(yīng)，為氮肥減施提供依據(jù)；（3）本研究在高氮情況下，施肥處理水稻產(chǎn)量均較高（平均為10.6 t·hm-2，未在結(jié)果部分展示），此時(shí)土壤中的氮素并不能充分滿足水稻對(duì)氮素的吸收，進(jìn)而能夠凸顯肥料的增產(chǎn)效應(yīng)，也說明高產(chǎn)稻田同樣需要培肥。其他各化學(xué)指標(biāo)不同水平下施肥對(duì)產(chǎn)量的影響結(jié)果各不相同，通過隨機(jī)森林分析結(jié)果可知：土壤質(zhì)地、全氮、緩效鉀和速效鉀對(duì)水稻產(chǎn)量響應(yīng)比影響相對(duì)較大。這也表明土壤一致的情況下，要注重氮肥的合理施用；其次是鉀素，因?yàn)殁浰匾詿o機(jī)形態(tài)存在土壤中，易隨水遷移，因此也要注重鉀肥的投入，尤其是南方缺鉀地區(qū)，結(jié)果顯示緩效鉀含量較低情況下的增產(chǎn)效果較高也印證這一結(jié)果。

水稻產(chǎn)量響應(yīng)比與肥料農(nóng)學(xué)效率之間呈極顯著正相關(guān)關(guān)系，這說明施肥在增產(chǎn)的同時(shí)也能增加肥效。在增加相同單位的肥料農(nóng)學(xué)效率情況下，化肥與有機(jī)肥配施處理對(duì)水稻產(chǎn)量的響應(yīng)比提高幅度高于單施化肥處理。這說明在等養(yǎng)分投入條件下，合理配施養(yǎng)分對(duì)產(chǎn)量的提高至關(guān)重要，該結(jié)果不僅與前人的研究結(jié)果一致[49]，同時(shí)也為我們當(dāng)前提倡的化肥減施的情況下保持增效提供理論依據(jù)[50]。另外，還應(yīng)當(dāng)結(jié)合施肥對(duì)產(chǎn)量響應(yīng)比的影響因素，綜合考慮水稻種植區(qū)域、種植制度和土壤理化性質(zhì)后進(jìn)行合理施肥，最終提高肥料的農(nóng)學(xué)效率，達(dá)到增產(chǎn)增效的目的。

4 結(jié)論

4.1 1988—2017年間，施肥能夠顯著提高水稻的產(chǎn)量，相比不施肥能增產(chǎn)80.8%，但是增產(chǎn)的效應(yīng)呈逐漸減弱趨勢(shì)。與化肥單施相比，化肥與有機(jī)肥配施是提高和維持水稻高產(chǎn)的重要措施。

4.2 施肥對(duì)水稻產(chǎn)量的提高效應(yīng)與水稻種植區(qū)域、土壤全氮、種植制度和土壤鉀素含量等密切相關(guān)。建議在指導(dǎo)農(nóng)民施肥時(shí)應(yīng)結(jié)合上述指標(biāo)進(jìn)行合理推薦。

4.3 水稻產(chǎn)量響應(yīng)比與肥料農(nóng)學(xué)效率之間呈極顯著正相關(guān)關(guān)系，增加相同單位的肥料農(nóng)學(xué)效率情況下，化肥與有機(jī)肥配施處理對(duì)水稻產(chǎn)量的響應(yīng)比提高量高于單施化肥處理，為我國(guó)當(dāng)前提倡化肥減施情況下如何實(shí)現(xiàn)水稻增產(chǎn)增效提供理論依據(jù)。

致謝：Muhammad Qaswar對(duì)本文英文摘要的修改給予了指導(dǎo)，在此致以誠(chéng)摯的感謝！

[1] PATEL D P, DAS A, MUNDA G C, GHOSH P K, BORDOLOI J S, KUMAR M. Evaluation of yield and physiological attributes of high-yielding rice varieties under aerobic and flood-irrigated management practices in mid-hills ecosystem., 2010, 97: 1269-1276.

[2] ZHANG H M, XU M G, SHI X J, LI Z Z, HUANG Q H, WANG X J. Rice yield, potassium uptake and apparent balance under long-term fertilization in rice-based cropping systems in southern China., 2010, 88: 341-349.

[3] 曾祥明, 韓寶吉, 徐芳森, 黃見良, 蔡紅梅, 石磊. 不同基礎(chǔ)地力土壤優(yōu)化施肥對(duì)水稻產(chǎn)量和氮肥利用率的影響. 中國(guó)農(nóng)業(yè)科學(xué), 2012, 45(14): 2886-2894.

ZENG X M, HAN B J, XU F S, HUANG J L, CAI H M, SHI L. Effect of optimized fertilization on grain yield of rice and nitrogen use efficiency in paddy fields with different basic soil fertilities., 2012, 45(14): 2886-2894. (in Chinese)

[4] HAIDER I K. Appraisal of biofertilizers in rice: To supplement inorganic chemical fertilizer., 2018, 25(6): 357-362.

[5] 孫浩燕, 李小坤, 任濤, 叢日環(huán), 魯劍巍. 淺層施肥對(duì)水稻苗期根系生長(zhǎng)及分布的影響. 中國(guó)農(nóng)業(yè)科學(xué), 2014, 47(12), 2476-2484.

SUN H Y, LI X K, REN T, CONG R H, LU J W. Effects of fertilizer in shallow soils on growth and distribution of rice roots at seedling stage., 2014, 47(12): 2476-2484. (in Chinese)

[6] 張國(guó)榮, 李菊梅, 徐明崗, 高菊生, 谷思玉. 長(zhǎng)期不同施肥對(duì)水稻產(chǎn)量及土壤肥力的影響. 中國(guó)農(nóng)業(yè)科學(xué), 2009, 42(2): 543-551.

ZHANG G R, LI J M, XU M G, GAO J S, GU Si Y. Effects of chemical fertilizer and organic manure on rice yield and soil fertility., 2009, 42(2): 543-551. (in Chinese)

[7] YADAV R L. Assessing on-farm efficiency and economics of fertilizer N, P and K in rice wheat systems of India., 2003, 81(1): 39-51.

[8] HAEFELE S M, WOPEREIS M C S. Spatial variability of indigenous supplies for N, P and K and its impact on fertilizer strategies for irrigated rice in West Africa., 2005, 270(1): 57-72.

[9] 吳萍萍, 劉金劍, 周毅, 謝小立, 沈其榮, 郭世偉. 長(zhǎng)期不同施肥制度對(duì)紅壤稻田肥料利用率的影響. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2008, 14(2): 277-283.

WU P P, LIU J J, ZHOU Y, XIE X L, SHEN Q R, GUO S W. Effects of different long term fertilizing systems on fertilizer use efficiency in red paddy soil., 2008, 14(2): 277-283. (in Chinese)

[10] 苑俊麗梁新強(qiáng)李亮葉玉適傅朝棟宋清川.中國(guó)水稻產(chǎn)量和氮素吸收量對(duì)高效氮肥響應(yīng)的整合分析. 中國(guó)農(nóng)業(yè)科學(xué), 2014, 47(17): 3414-3423.

YUAN J L, LIANG X Q, LI L, YE Y S, FU C D, SONG Q C. Response of rice yield and nitrogen uptake to enhanced efficiency nitrogen fertilizer in China: A meta-analysis., 2014, 47(17): 3414-3423. (in Chinese)

[11] 張福鎖, 王激清, 張衛(wèi)峰, 崔振嶺, 馬文奇, 陳新平, 江榮風(fēng). 中國(guó)主要糧食作物肥料利用率現(xiàn)狀與提高途徑. 土壤學(xué)報(bào), 2008, 45(5): 915-924.

ZHANG F S, WANG J Q, ZHANG W F, CUI Z L, MA W Q, CHEN X P, JIANG R F. Nutrient use efficiencies of major cereal crops in China and measures for improvement., 2008, 45(5): 915-924. (in Chinese)

[12] 胡建利, 王德建, 王燦, 孫瑞娟. 不同施肥方式對(duì)水稻產(chǎn)量構(gòu)成及其穩(wěn)定性的影響. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2009, 17(1): 48-53.

HU J L, WANG D J, WANG C, SUN R J. Effect of different fertilization systems on rice yield components and their stability., 2009, 17(1): 48-53. (in Chinese)

[13] 包耀賢, 徐明崗, 呂粉桃, 黃慶海, 聶軍, 張會(huì)民, 于寒青. 長(zhǎng)期施肥下土壤肥力變化的評(píng)價(jià)方法. 中國(guó)農(nóng)業(yè)科學(xué), 2012, 45(20): 4197-4204.

BAO Y X, XU M G, Lü F T, HUANG Q H, NIE J, ZHANG H M, YU H Q. Evaluation method on soil fertility under long-term fertilization., 2012, 45(20): 4197-4204. (in Chinese)

[14] 吳良泉, 武良, 崔振嶺, 陳新平, 張福鎖. 中國(guó)水稻區(qū)域氮磷鉀肥推薦用量及肥料配方研究. 中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào), 2016, 21(9): 1-13.

WU L Q, Wu L, CUI Z L, CHEN X P, ZHANG F S. Studies on recommended nitrogen, phosphorus and potassium application rates and special fertilizer formulae for different rice production regions in china., 2016, 21(9): 1-13. (in Chinese)

[15] 方暢宇, 屠乃美, 張清壯, 易鎮(zhèn)邪. 不同施肥模式對(duì)稻田土壤速效養(yǎng)分含量及水稻產(chǎn)量的影響. 土壤, 2018, 50(3): 462-468.

FANG C Y, TU N M, ZHANG Q Z, YI Z X. Effects of fertilization modes on available nutrient contents of reddish paddy soils and rice yields., 2018, 50(3): 462-468. (in Chinese)

[16] 銀敏華, 李援農(nóng), 陳朋朋, 徐路全, 申勝龍, 王星壵. 基于Meta-analysis的中國(guó)北方地區(qū)免耕玉米產(chǎn)量效應(yīng)研究. 中國(guó)農(nóng)業(yè)科學(xué), 2018, 51(5): 843-854.

YIN M H, LI Y N, CHEN P P, XU L Q, SHEN S L, WANG X Y. Effect of no-tillage on maize yield in northern region of China: A meta-analysis., 2018, 51(5): 843-854. (in Chinese)

[17] MELLOR A, HAYWOOD A, STONE C, JONES S. The performance of random forests in an operational setting for large area sclerophyll forest classification., 2013, 5(6): 2838-2856.

[18] 全國(guó)農(nóng)業(yè)技術(shù)推廣服務(wù)中心. 土壤分析技術(shù)規(guī)范. 2版. 北京: 中國(guó)農(nóng)業(yè)出版社, 2006.

National Agricultural Technology Extension and Service Center.. Beijing: China Agriculture Press, 2006. (in Chinese)

[19] CURTIS P S, WANG X. A meta-analysis of elevated CO2effects on woody plant mass, form, and physiology., 1998, 113(3): 299-313.

[20] 任鳳玲, 張旭博, 孫楠, 徐明崗, 柳開樓. 施用有機(jī)肥對(duì)中國(guó)農(nóng)田土壤微生物量影響的整合分析. 中國(guó)農(nóng)業(yè)科學(xué), 2018, 51(1): 119-128.

REN F L, ZHANG X B, SUN N, XU M G, LIU K L. A meta-analysis of manure application impact on soil microbial biomass across China's croplands., 2018, 51(1): 119-128. (in Chinese)

[21] LUO Y Q, HUI D F, ZHANG D Q. Elevated CO2stimulates net accumulations of carbon and nitrogen in land ecosystems: A meta-analysis., 2006, 87(1): 53-63.

[22] HEDGES L V, GUREVITCH J, CURTIS P S. The meta-analysis of response ratios in experimental ecology., 1999, 80(4): 1150-1156.

[23] 蔡岸冬, 張文菊, 楊品品, 韓天富, 徐明崗. 基于meta-analysis研究施肥對(duì)中國(guó)農(nóng)田土壤有機(jī)碳及其組分的影響. 中國(guó)農(nóng)業(yè)科學(xué), 2015, 48(15): 2995-3004.

CAI A D, ZHANG W J, YANG P P, HAN T F, XU M G. Effect degree of fertilization practices on soil organic carbon and fraction of croplands in China-based on meta-analysis., 2015, 48(15): 2995-3004. (in Chinese)

[24] CHEN Y, CAMPS-ARBESTAIN M, SHEN Q H, SINGH B, CAYUELA M L. The long-term role of organic amendments in building soil nutrient fertility: A meta-analysis and review., 2018, 111(2/3): 103-125.

[25] WANG J Z, WANG X J, XU M G, FENG G, ZHANG W J, LU C A. Crop yield and soil organic matter after long-term straw return to soil in China., 2015, 102(3): 371-381.

[26] 徐霞, 趙亞南, 黃玉芳, 閆軍營(yíng), 葉優(yōu)良. 不同地力水平下的小麥?zhǔn)┓市?yīng). 中國(guó)農(nóng)業(yè)科學(xué), 2018, 51(21): 4076-4086.

XU X, ZHAO Y L, HUANG Y F, YAN J Y, YE Y L. Fertilization effect of wheat under different soil fertilities., 2018, 51(21): 4076-4086. (in Chinese)

[27] RODRIGUEZ-GALIANO V F, GHIMIRE B, ROGAN J, CHICA- OLMO M, RIGLI-SANCHEZ J P. An assessment of the effectiveness of a random forest classifier for land-cover classification., 2012, 67(1): 93-104.

[28] LARSON B A, FRISVOLD G B. Fertilizers to support agricultural development in sub-Saharan Africa: What is needed and why., 1996, 21(6): 509-525.

[29] MUELLER N D, GERBER J S, JOHNSTON M, Ray D K, RAMANKUTTY N, FOLEY J A. Closing yield gaps through nutrient and water management., 2012, 490(7419): 254-257.

[30] CHIANU J N, CHIANU J N, MAIRURA F. Mineral fertilizers in the farming systems of sub-Saharan Africa. A review., 2012, 32(2): 545-566.

[31] 石元亮, 王玲莉, 劉世彬, 聶鴻光. 中國(guó)化學(xué)肥料發(fā)展及其對(duì)農(nóng)業(yè)的作用. 土壤學(xué)報(bào), 2008, 45(5): 852-864.

SHI Y L, WANG L L, LIU S B, NIE H G. Development of chemical fertilizer industry and its effect on agriculture of China., 2008, 45(5): 852-864. (in Chinese)

[32] ALLEY M M, VANLAUWE B. The role of fertilizers in integraded plant nutrient management., 2009, 15: 1-59.

[33] ZHANG W F, CAO G X, LI X L, ZHANG H Y, WANG C, LIU Q Q, CHEN X P, CUI Z L, SHEN J B, JIANG R F, MI G H, MIAO Y X, ZHANG F S, DOU Z X. Closing yield gaps in China by empowering smallholder farmers., 2016, 537(7622): 671-674.

[34] CHEUNG F. Yield: The search for the rice of the future., 2014, 514(7524): S60-S61.

[35] 黃耀, 孫文娟. 近20年來中國(guó)大陸農(nóng)田表土有機(jī)碳含量的變化趨勢(shì). 科學(xué)通報(bào), 2007, 51: 750-763.

HUANG Y, SUN W J. Change trend of organic carbon in farmland surface soil of main land China during the past 20 years., 2007, 51: 750-763. (in Chinese)

[36] 張智, 王偉妮, 李昆, 馬紅菊, 茍曦, 魯劍巍, 李小坤, 叢日環(huán). 四川省不同區(qū)域水稻氮肥施用效果研究. 土壤學(xué)報(bào), 2015, 52(1): 234-241.

ZHANG Z, WANG W N, LI K, MA H J, GOU X, LU J W, LI X K, CONG R H. Effects of nitrogen fertilization on rice in different regions of Sichuan Province., 2015, 52(1): 234-241. (in Chinese)

[37] 任意, 張淑香, 穆蘭, 田有國(guó), 盧昌艾. 我國(guó)不同地區(qū)土壤養(yǎng)分的差異及變化趨勢(shì). 中國(guó)土壤與肥料, 2009(6): 13-17.

REN Y, ZHANG S X, MU L, TIAN Y G, LU C A. Change and difference of soil nutrients for various regions in China., 2009(6): 13-17 . (in Chinese)

[38] YANG X L, GAO W S, ZHANG M, CHEN Y Q, SUI P. Reducing agricultural carbon footprint through diversified crop rotation systems in the north China plain., 2014, 76: 131-139.

[39] 湯勇華, 黃耀. 中國(guó)大陸主要糧食作物地力貢獻(xiàn)率和基礎(chǔ)產(chǎn)量的空間分布特征. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2009, 28(5): 1070-1078.

TANG Y H, HUANG Y. Spatial distribution characteristics of the percentage of soil fertility contribution and its associated basic crop yield in mainland China., 2009, 28(5): 1070-1078. (in Chinese)

[40] GUO J H, LIU X J, ZHANG Y, SHEN J L, HAN W X, ZHANG W F, CHEISTIE P, GOULDING K W T, VITOUSEK P M, ZHANG F S. Significant acidification in major Chinese croplands., 2010, 327 (5968): 1008-1010.

[41] YOUSAF M, LI J F, LU J W, REN T, CONG R H, FAHAD S, LI X K. Effects of fertilization on crop production and nutrient-supplying capacity under rice-oilseed rape rotation system., 2017, 7(1): 1-9.

[42] 黃欠如, 胡鋒, 李輝信, 賴濤, 袁穎紅. 紅壤性水稻土施肥的產(chǎn)量效應(yīng)及與氣候、地力的關(guān)系. 土壤學(xué)報(bào), 2006, 43(6): 926-929.

HUANG Q R, HU F, LI H X, LAI T, YUAN Y H. Crop yield response to fertilization and its relations with climate and soil fertility in red paddy soil., 2006, 43(6): 926-929. (in Chinese)

[43] 徐明崗, 李冬初, 李菊梅, 秦道珠, 八木一行, 寶川靖和. 化肥有機(jī)肥配施對(duì)水稻養(yǎng)分吸收和產(chǎn)量的影響. 中國(guó)農(nóng)業(yè)科學(xué), 2008, 41(10): 3133-3139.

XU M G, LI D C, LI J M, QIN D Z, KAZUYUKI Y, YASUKAZU H. Effects of organic manure application with chemical fertilizers on nutrient absorption and yield of rice in Hunan of Southern China., 2008, 41(10): 3133-3139. (in Chinese)

[44] LIN Y R, WATTS D B, SANTEN E V, CAO G Q. Influence of poultry litter on crop productivity under different field conditions: A meta-analysis., 2018, 110(3): 807-818.

[45] 呂貽忠. 土壤學(xué). 北京: 中國(guó)農(nóng)業(yè)出版社, 2006.

Lü Y Z.Beijing: China Agricultural Press, 2006. (in Chinese)

[46] TIROLPADRE A, LADHA J K, REGMI A P, BHANDARI A L, INUBUSHI K. Organic amendments affect soil parameters in two long-term rice-wheat experiments., 2007, 71(2): 442-452.

[47] 彭少兵, 黃見良, 鐘旭華, 楊建昌, 王光火, 鄒應(yīng)斌, 張福鎖, 朱慶森, ROLAND B, CHRISTIAN W. 提高中國(guó)稻田氮肥利用率的研究策略. 中國(guó)農(nóng)業(yè)科學(xué), 2002, 35(9): 1095-1103.

PENG S B, HUANG J L, ZHONG X H, YANG J C, WANG G H, ZHOU Y B, ZHANG F S, ZHU Q S, ROLAND B, CHRISTIAN W. Research strategy in improving fertilizer-nitrogen use efficiency of irrigated rice in China., 2002, 35(9): 1095-1103. (in Chinese)

[48] 馮洋, 陳海飛, 胡孝明, 蔡紅梅, 徐芳森. 高、中、低產(chǎn)田水稻適宜施氮量和氮肥利用率的研究. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2014, 20(1): 7-16.

FENG Y, CHEN H F, HU X M, CAI H M, XU F S. Optimal nitrogen application rates on rice grain yield and nitrogen use efficiency in high, middle and low-yield paddy fields., 2014, 20(1): 7-16. (in Chinese)

[49] 冀建華, 侯紅乾, 劉益仁, 劉秀梅, 馮兆濱, 劉光榮, 楊濤, 李文娟. 長(zhǎng)期施肥對(duì)雙季稻產(chǎn)量變化趨勢(shì)、穩(wěn)定性和可持續(xù)性的影響. 土壤學(xué)報(bào), 2015, 52(3): 607-619.

JI J H, HOU H Q, LIU Y R, LIU X M, FENG Z B, LIU G R, YANG T, LI W J. Effects of long-term fertilization on yield variation trend, yield stability and sustainability in the double cropping rice system., 2015, 52(3): 607-619. (in Chinese)

[50] 白由路. 高效施肥技術(shù)研究的現(xiàn)狀與展望. 中國(guó)農(nóng)業(yè)科學(xué), 2018, 51(11): 2116-2125.

BAI Y L. The situation and prospect of research on efficient fertilization., 2018, 51(11): 2116-2125. (in Chinese)

Variation in Rice Yield Response to Fertilization in China: Meta-analysis

HAN TianFu1, MA ChangBao2, HUANG Jing1, 3, LIU KaiLou1, 4, XUE YanDong2, LI DongChu1, 3, LIU LiSheng1, 3, ZHANG Lu1, 3, LIU ShuJun1, 3, ZHANG HuiMin1, 3

(1Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences/National Engineering Laboratory for Improving Quality of Arable Land, Beijing 100081;2Center of Arable Land Quality Monitoring and Protection, Ministry of Agriculture and Rural Affairs, Beijing 100125;3Red Soil Experimental Station of CAAS in Hengyang/National Observation and Research Station of Farmland Ecosystem in Qiyang, Qiyang 426182, Hunan;4Jiangxi Institute of Red Soil/ National Engineering and Technology Research Center for Red Soil Improvement, Nanchang 330046)

【】A meta study was conducted to investigate the comprehensive effect of fertilization on rice yield in Chinese paddy soils during the past 30 years, and to provide a theoretical basis for the scientific correct application of fertilizers in rice cultivation areas.【】Based on the long-term paddy soil monitoring sites from Ministry of Agriculture and Rural Affairs, we conducted meta-analysis to investigate the rice yield response to no fertilization versus fertilization in different agro-climatic regions.】Rice yield in the past 10 years (2008-2017) was significantly higher than the corresponding rice yield in 1988-1997 and 1998-2007, regardless of fertilization. The increase of rice yield with fertilization in southwest of China was by 98.5%, which was significantly higher than that of in north of China (70.3%). Fertilization increased rice yield by 99.1%, 84.2% and 78.1% during 1988-1997, 1998-2007 and 2008-2017, respectively. For the cropping system, the increase of rice yield under triple cropping system (92.0%) was significantly higher than that under single cropping system (76.2%) and double cropping system (81.9%). Fertilization increased rice yield by 85.9% under double rice cropping system, by 75.9% under single cropping system, and by 79.5% under other cropping system. Compared with no fertilizer, chemical plus organic fertilizer application increased rice yield by 88.3%, which was higher than that of single chemical fertilizer application (76.6%). Fertilization significantly increased rice yield in clay soil by 92.0%, compared with no fertilization, which significantly higher than that in sandy soil (58.0%) and loam soil (77.5%). With the increase of soil organic matter and available phosphorus, the increasing trend of fertilization on rice yield was decreased compared with no fertilization. Under higher soil pH (＞7.5) and lower soil total nitrogen (＜1.5 g·kg-1) and slow available potassium (＜150 mg·kg-1), the rice yield increasing was more than that of corresponding the rest of level. Random forest analysis showed that the region, soil total nitrogen and cropping system had greater impact on the response ratio () of rice yield. In addition, the agronomic efficiency of fertilizer was positively correlated with rice yield. 【】Although the trend of increasing rice yield by fertilization was decreasing at present, but combined appropriate chemical plus organic fertilizer, especially in southwest of China, were important measures to improve and maintain high rice yield. Base on the cropping system, combining soil texture, soil nitrogen and potassium content should be the main basis for fertilizer input in different rice cultivation areas.

fertilization; rice yield; meta-analysis; response ratio; agronomic efficiency

10.3864/j.issn.0578-1752.2019.11.007

2018-12-04；

2019-01-18

國(guó)家重點(diǎn)研發(fā)計(jì)劃“糧食豐產(chǎn)增效科技創(chuàng)新”重點(diǎn)專項(xiàng)（2016YFD0300901）、國(guó)家自然科學(xué)基金（41671301，41371293）

韓天富，E-mail：hantianfu123@126.com。通信作者張會(huì)民，E-mail：zhanghuimin@caas.cn

（責(zé)任編輯 李云霞）