索炎炎 張 翔 司賢宗 余 瓊 李 亮 余 輝

(1河南省農(nóng)業(yè)科學(xué)院植物營養(yǎng)與資源環(huán)境研究所，河南 鄭州 450000；2正陽縣花生研究所，河南 正陽 463600)

土壤是一個不均一變化的連續(xù)體，受氣候、生物、母質(zhì)、地形等自然因素和人為活動的影響，具有較高的時空變異性[1-2]。同一時期，即使土壤類型和質(zhì)地相同的土壤也存在明顯的空間變異[3-4]。土壤養(yǎng)分是植物生長發(fā)育所必需的營養(yǎng)元素，是土壤肥力的物質(zhì)基礎(chǔ)，了解土壤肥力尤其是土壤養(yǎng)分的空間變異特征，可有效揭示土壤養(yǎng)分空間分布狀況，為土壤養(yǎng)分分區(qū)管理和精準(zhǔn)施肥提供基本依據(jù)。

地理信息系統(tǒng)(geographic information system，GIS)和地統(tǒng)計(jì)學(xué)技術(shù)相結(jié)合是研究土壤空間異質(zhì)性最有效的方法之一。近幾年，國內(nèi)外學(xué)者采用該方法對土壤的空間異質(zhì)性特征進(jìn)行了大量研究，如徐新朋等[5]應(yīng)用GIS 和地統(tǒng)計(jì)學(xué)技術(shù)研究了吉林省玉米種植區(qū)土壤養(yǎng)分空間變異特征；楊之江等[3]基于該方法研究了湖南省長沙縣北山鎮(zhèn)稻田土壤養(yǎng)分與重金屬的空間分布特征；Zhu 等[6]運(yùn)用該方法對水稻-小麥輪作農(nóng)田土壤氮素的時空變異進(jìn)行了研究。前人開展的多為單項(xiàng)土壤因子的空間分布特征研究，但土壤肥力是多種物理、化學(xué)和生物因子的綜合反映[7-8]，單個肥力指標(biāo)不能代表土壤肥力水平，更不能以此作為制定調(diào)節(jié)土壤肥力綜合措施的依據(jù)。因此，需要對土壤肥力進(jìn)行綜合評價。

模糊隸屬度函數(shù)是重要的統(tǒng)計(jì)和評價方法[9-11]，國內(nèi)外學(xué)者利用該方法開展了大量土壤肥力的量化綜合評價研究，如Bijanzadeh 等[12]利用隸屬度函數(shù)法研究了伊朗設(shè)拉子?xùn)|部麥田肥力水平及其與土壤鹽分的關(guān)系；李強(qiáng)等[13]和呂小娜等[14]采用GIS 和模糊數(shù)學(xué)相結(jié)合的方法，對植煙土壤的肥力空間變異進(jìn)行了分析與綜合評價；徐志超等[15]基于GIS 和土壤采樣數(shù)據(jù)，利用隸屬度函數(shù)法評價并分析了長三角典型區(qū)占補(bǔ)耕地的土壤肥力狀況及其空間特征。但是，土壤肥力綜合指標(biāo)的空間分布特征在作物種植前后的變化研究鮮有報(bào)道。

花生(Arachis hypogaeaL.)是我國重要的油料作物，種植面積達(dá)461.97 萬hm2，年產(chǎn)超過1 733 萬t，總產(chǎn)居全國油料作物之首，2018年河南省花生種植面積達(dá)120.32 萬hm2，已經(jīng)成為河南省第三大作物，總產(chǎn)和單產(chǎn)均居全國首位[16]。砂姜黑土是黃淮海平原的主要耕作土壤類型之一，也是河南省花生種植區(qū)的主要土壤類型。本研究以河南省正陽縣為研究區(qū)域，利用GIS 和地統(tǒng)計(jì)學(xué)相結(jié)合的方法，對該區(qū)花生種植前和成熟期土壤養(yǎng)分狀況及其空間變異特征進(jìn)行分析，并采用模糊數(shù)學(xué)的隸屬度函數(shù)法對該區(qū)土壤肥力水平進(jìn)行量化綜合評價，以明確花生生產(chǎn)中的施肥狀況，旨在為該地區(qū)花生優(yōu)質(zhì)高產(chǎn)的合理施肥和培肥地力提供科學(xué)指導(dǎo)。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)設(shè)置在河南省駐馬店正陽縣河南省南部，位于(32°16′～32°47′N，114°12′～114°53′E)，南鄰淮河，屬大陸性濕潤季風(fēng)氣候，年均降雨量938 mm，年平均氣溫15℃，四季分明，氣候溫和，適合多種農(nóng)作物生長。該研究區(qū)花生種植面積超過11.34萬hm2，占全縣總耕地面積的74%，每年6—10月為花生生育期，花生生育期內(nèi)平均氣溫和總積溫分別為25.36℃和3 220.63℃；花生生育期內(nèi)降水量為538.44 mm。該研究區(qū)土壤類型以砂姜黑土為主，花生種植方式為起壟種植，種植花生品種主要為遠(yuǎn)雜9102。研究區(qū)種植方式主要為小麥-花生兩熟制，根據(jù)調(diào)查發(fā)現(xiàn)當(dāng)?shù)剞r(nóng)民的耕作管理方式、周年施肥習(xí)慣等基本一致。

1.2 土壤樣品采集

于2017年5月22 —26日(花生種植前)和9月12—15日(花生成熟期)，在前茬為小麥后季種植花生的同一田塊采集土壤。田間取樣時，進(jìn)行GPS 定位，記錄采樣點(diǎn)的經(jīng)緯度，為保證樣品的代表性，在定位的田塊地里采用“S”法采集0～20 cm 的耕層土壤樣品10～15 個，制成1 個混合土樣，四分法保留1 kg，風(fēng)干后，過20 目和100 目篩，備用。花生種植前和收獲后均定位采集48 個土壤樣品，最終獲得96 個土壤樣品。根據(jù)花生種植面積，每個種植花生的鄉(xiāng)鎮(zhèn)均設(shè)有采樣點(diǎn)，土壤采樣點(diǎn)分布見圖1。

1.3 測定項(xiàng)目與方法

參照文獻(xiàn)[17]測定土壤pH 值以及有機(jī)質(zhì)、全氮、有效磷、速效鉀、有效硫、有效鐵、有效鋅含量。其中pH 值采用FE20-K pH 計(jì)(梅特勒-托利多儀器，上海)測定(水∶土=2.5∶1)；有機(jī)質(zhì)含量采用重鉻酸鉀容量法測定；全氮含量利用2300 全自動凱氏定氮儀(Foss公司，瑞典)測定；有效磷含量采用鉬銻抗顯色法測定；速效鉀含量利用FP6410 火焰光度計(jì)(上海欣益儀器儀表有限公司)測定；有效硫含量采用硫酸鋇比濁法測定；有效鐵和有效鋅含量采用二乙烯三胺五乙酸(diethylene triamine pentaacetic acid，DTPA)溶液浸提，利用HitachiZ-2000 原子吸收分光光度計(jì)(北京天美科學(xué)儀器有限公司)測定。

1.4 土壤肥力綜合評價方法

本研究采用模糊數(shù)學(xué)隸屬度函數(shù)模型進(jìn)行土壤肥力評價。由于評價因素指標(biāo)之間缺乏可比性，因此在評價土壤肥力之前先利用隸屬度函數(shù)將各評價指標(biāo)進(jìn)行歸一化處理。常用的隸屬函數(shù)主要有S 型隸屬函數(shù)和拋物線型隸屬函數(shù)，分別見公式(1)和公式(2)[18]:

根據(jù)前人研究成果[15，18-19]，結(jié)合研究區(qū)土壤中養(yǎng)分的實(shí)際情況，土壤養(yǎng)分指標(biāo)中有機(jī)質(zhì)含量、全氮含量、有效磷含量、速效鉀含量、有效硫含量、有效鐵含量、有效鋅含量隸屬度函數(shù)采用S 型隸屬函數(shù)，pH 值采用拋物線型隸屬度函數(shù)。土壤肥力各評價指標(biāo)的隸屬度函數(shù)及閾值見表1。

采用相關(guān)系數(shù)法確定評價因子權(quán)重系數(shù)[20]。首先計(jì)算單項(xiàng)評價指標(biāo)之間的相關(guān)系數(shù)，然后計(jì)算某一評價指標(biāo)與其他指標(biāo)之間相關(guān)系數(shù)的平均值，用該平均值與所有評價指標(biāo)相關(guān)系數(shù)平均值總和的比值，作為該評價指標(biāo)的權(quán)重。

表1 土壤肥力評價指標(biāo)隸屬度函數(shù)類型及閾值Table 1 Membership function types and related parameters of soil fertility index

根據(jù)模糊數(shù)學(xué)中的加成法原則，用公式(3)計(jì)算土壤肥力綜合評價指數(shù)(integrated fertility index，IFI):

式中，Wi和Ni分別表示第i個指標(biāo)的權(quán)重值和隸屬度值。

1.5 數(shù)據(jù)統(tǒng)計(jì)與分析

采用Microsoft Excel 2013 進(jìn)行數(shù)據(jù)處理。異常值識別采用樣本平均值±3 倍標(biāo)準(zhǔn)差法，大于平均值＋3倍標(biāo)準(zhǔn)差和小于平均值-3 倍標(biāo)準(zhǔn)差的異常值分別用數(shù)據(jù)正常值的最大值、最小值進(jìn)行替換處理。采用軟件SPSS 16.0 對數(shù)據(jù)進(jìn)行描述性分析，采用軟件GS＋7.0 進(jìn)行半方差函數(shù)計(jì)算和理論模型擬合，并利用軟件ArcGIS 10.2.2 地統(tǒng)計(jì)模塊的普通克里格(Kriging)插值方法，繪制土壤綜合肥力指標(biāo)的空間分布圖。

2 結(jié)果與分析

2.1 土壤養(yǎng)分含量描述性統(tǒng)計(jì)分析

由表2可知，花生種植前土壤pH 值范圍為4.73～6.22，平均值5.34，呈酸性；有機(jī)質(zhì)、全氮、有效磷、速效鉀、有效硫、有效鐵和有效鋅含量范圍分別為12.80 mg·kg-1～24.16 g·kg-1、0.83～1.71 g·kg-1、11.52～127.81 mg·kg-1、73.01～236.00 mg·kg-1、19.60～161.00 mg·kg-1、73.10～344.00 mg·kg-1和0.51～1.33 mg·kg-1，平均值分別為17.15 g·kg-1、1.18 g·kg-1、53.89 mg·kg-1、125.13 mg·kg-1、74.33 mg·kg-1、178.00 mg·kg-1和0.78 mg·kg-1?；ㄉ墒炱谕寥纏H值范圍為4.84～6.12，平均值5.31，也呈酸性；有機(jī)質(zhì)、全氮、有效磷、速效鉀、有效硫、有效鐵和有效鋅含量范圍分別為12.30～24.16 g·kg-1和0.89～1.66 g·kg-1、15.60～137.40 mg·kg-1、62.04～218.00 mg·kg-1、35.30～121.10 mg·kg-1、105.00～329.00 mg·kg-1和0.51～1.47 mg·kg-1，平均值分別為17.23 g·kg-1、1.30 g·kg-1、62.71 mg·kg-1、128.23 mg·kg-1、65.89 mg·kg-1、190.22 mg·kg-1和0.85 mg·kg-1。與種植前相比，花生成熟期土壤全氮、有效磷、有效鐵和有效鋅含量明顯增加，增幅分別為10.17%、16.37%、6.87%和8.97%，有機(jī)質(zhì)和速效鉀含量也有所升高，但增幅較小，而土壤有效硫含量則出現(xiàn)明顯下降，土壤pH 值也有所降低。

根據(jù)全國第二次土壤普查養(yǎng)分分級標(biāo)準(zhǔn)[21](表3)可知，試驗(yàn)區(qū)土壤pH 值呈酸性，有效磷、有效硫和有效鐵含量屬于一級水平，全氮和速效鉀含量屬于三級水平，有機(jī)質(zhì)和有效鋅含量屬于四級水平。從變異系數(shù)來看，花生種植前和成熟期土壤各養(yǎng)分各指標(biāo)變異系數(shù)變化不大，花生種植前和成熟期土壤pH 值變異系數(shù)分別為6.18%和6.36%，屬弱變異；其余養(yǎng)分指標(biāo)均呈中等強(qiáng)度變異，其中花生種植前和成熟期有效磷的變異系數(shù)最大，分別為40.55%和44.94%，而有機(jī)質(zhì)的變異系數(shù)最小，分別為12.71%和12.58%。

表2 花生種植前和成熟期土壤pH 值和養(yǎng)分含量描述性統(tǒng)計(jì)Table 2 Descriptive statistics of soil pH value and nutrient contents at pre-planting and maturity of peanut

表3 土壤養(yǎng)分及pH 值分級標(biāo)準(zhǔn)Table 3 Classification criterions for soil nutrients and pH value

2.2 花生種植前和成熟期土壤養(yǎng)分狀況綜合評價

2.2.1 指標(biāo)權(quán)重值的確定 土壤肥力評價指標(biāo)的相關(guān)系數(shù)平均值及權(quán)重見表4。不同土壤肥力評價指標(biāo)對土壤肥力的貢獻(xiàn)存在差異，其中有機(jī)質(zhì)含量的權(quán)重最大，為0.216，說明有機(jī)質(zhì)含量在決定土壤肥力質(zhì)量變化中起最主要作用；速效鉀含量、全氮含量、有效磷含量的權(quán)重居中，分別為0.181、0.163 和0.164，pH值、有效鋅含量、有效硫含量和有效鐵含量的權(quán)重較小，分別為0.104、0.075、0.065 和0.033，說明這幾個指標(biāo)對土壤肥力的貢獻(xiàn)較小。

2.2.2 花生田土壤肥力綜合評價指標(biāo)的描述性分析 土壤肥力綜合評價指數(shù)的取值范圍為0～1，值越接近1，土壤肥力水平越高[20]。由圖2-A 可知，花生種植前和成熟期土壤肥力綜合指數(shù)(IFI)多集中于0.58～0.72 之間，土壤肥力屬于中等偏上水平?；ㄉN植前土壤IFI 介于0.44～0.83 之間，平均值為0.64；花生成熟期土壤IFI 介于0.42～0.83 之間，平均值為0.66。可見，花生成熟期土壤肥力綜合指標(biāo)的平均值略高于花生種植前。

表4 土壤肥力評價指標(biāo)相關(guān)系數(shù)平均值及權(quán)重Table 4 Mean of correlation coefficients and the weight of soil fertility indices

2.2.3 花生田土壤肥力綜合評價指標(biāo)的計(jì)算及等級劃分 采用累積曲線法[15，22]對正陽縣花生田土壤肥力進(jìn)行等級劃分。為使研究結(jié)果具有可比性，本研究對花生種植前和成熟期的IFI 采用統(tǒng)一的劃分標(biāo)準(zhǔn)，具體是將花生種植前和成熟期的值合并在一起，繪制累積分布圖(圖2-B)，然后根據(jù)累積曲線的突變點(diǎn)劃分IFI 的分界點(diǎn)，最終將研究區(qū)花生田土壤肥力劃分為4 個等級(表5)。

表5 研究區(qū)土壤肥力等級劃分Table 5 Grades of soil integrated fertility index of study area

2.2.4 花生田土壤肥力綜合評價指標(biāo)的半方差分析 半方差分析結(jié)果如表6所示，花生種植前IFI 的最佳擬合模型為球狀模型，擬合的決定系數(shù)為0.674；花生成熟期IFI 的最佳擬合模型為高斯模型，擬合的決定系數(shù)為0.718，擬合效果優(yōu)于花生種植前IFI。塊金效應(yīng)為塊金值/基臺值，即C0/(C＋C0)表示隨機(jī)部分引起的空間變異占總變異的比例?；ㄉN植前和成熟期C0/(C＋C0)分別為17.09%和55.92%，說明花生成熟期土壤肥力空間分布比種植前受隨機(jī)因素(如施肥、耕作管理等人為因素)的影響大。

2.3 花生田土壤肥力綜合評價指標(biāo)的空間分布特征

花生種植前和成熟期土壤肥力綜合評價指標(biāo)(IFI)的空間分布特征如圖3所示，根據(jù)累積曲線獲取的分級標(biāo)準(zhǔn)，利用ArcGIS 10.2.2 軟件中空間分析(Spatial analyst)模塊進(jìn)行不同土壤肥力等級的面積分析，發(fā)現(xiàn)在花生種植前，土壤肥力為Ⅰ級(＜0.52)的土壤稀少，僅出現(xiàn)在汝南埠鎮(zhèn)；正陽縣大部分區(qū)域土壤肥力為Ⅱ級水平(0.52＜IFI＜0.67)，占正陽縣面積的84.22%；皮店鄉(xiāng)、陡溝鎮(zhèn)和彭橋鄉(xiāng)及蘭青鄉(xiāng)東南部花生田的土壤肥力為Ⅲ級水平(0.67＜IFI＜0.78)，占正陽縣面積的15.76%；未發(fā)現(xiàn)土壤肥力為Ⅳ級(＞0.78)的土壤。與花生種植前相比，花生成熟期正陽縣土壤肥力為Ⅱ級的土壤面積縮小，土壤肥力為Ⅲ級的土壤面積增大，分別占正陽縣面積的43.73%和56.10%；土壤肥力為Ⅰ級的土壤零星分布于王勿橋鄉(xiāng)和袁寨鄉(xiāng)；土壤肥力為Ⅳ級的土壤零星分布于陡溝鎮(zhèn)和慎水鄉(xiāng)。總之，成熟期花生田土壤肥力水平總體有所提高，且呈現(xiàn)出由南向北推移的現(xiàn)象。

表6 花生種植前和成熟期土壤肥力綜合指標(biāo)的半方差函數(shù)理論模型及其相關(guān)參數(shù)Table 6 Semi-variogram models of soil integrated fertility index and their parameters

3 討論

根據(jù)全國第二次土壤普查養(yǎng)分分級標(biāo)準(zhǔn)[21]，正陽縣花生田土壤有機(jī)質(zhì)含量屬于中等偏下水平，有效磷、有效硫和有效鐵含量屬于豐富水平，全氮和速效鉀含量屬于中等水平，有效鋅含量屬于缺鋅水平。與2005—2014年《測土配方施肥土壤基礎(chǔ)養(yǎng)分?jǐn)?shù)據(jù)集》[23]中正陽縣土壤pH 值(6.1)以及有機(jī)質(zhì)(15.02 g·kg-1)、全氮(0.838 g·kg-1)、有效磷(16.6 mg·kg-1)和速效鉀(88.4 mg·kg-1)含量平均值相比，本研究pH 值降低約0.8 個單位，其余指標(biāo)的測定結(jié)果均明顯增加。與種植前相比，花生成熟期土壤有效硫含量出現(xiàn)明顯下降，其余土壤養(yǎng)分指標(biāo)均有不同程度的提高，尤其是土壤全氮和有效磷含量增加明顯，增幅分別為10.17%和16.37%，但土壤有效鋅和有機(jī)質(zhì)含量仍處于較低水平。綜合分析，正陽縣花生種植戶在施肥上應(yīng)適當(dāng)控制氮磷用量，增加有機(jī)肥和碎肥的施用。

變異系數(shù)反映數(shù)據(jù)間的離散程度，可直觀反映出樣本的空間變異性大小。一般認(rèn)為，CV＜10%為弱變異，10%≤CV≤100%為中等強(qiáng)度變異，CV≥100%為強(qiáng)變異[24]?；ㄉN植前和成熟期土壤養(yǎng)分各指標(biāo)變異系數(shù)相差不大，pH 值變異系數(shù)分別為6.18%和6.36%，屬弱變異，其余養(yǎng)分指標(biāo)變異系數(shù)在12.58%～44.97%之間，屬中等強(qiáng)度變異，與鄭昊楠等[25]和舒田等[26]研究結(jié)果一致。pH 值變異強(qiáng)度較弱，與廖超林等[27]研究結(jié)果一致，這主要受土壤母質(zhì)影響，其余指標(biāo)中等強(qiáng)度變異可能是由農(nóng)戶的施肥水平、耕作栽培及管理措施不同所致，塊金效應(yīng)(55.92%)也進(jìn)一步證實(shí)了本結(jié)論。

土壤肥力是土壤為植物生長供應(yīng)和協(xié)調(diào)水、肥、氣、熱的能力。土壤肥力高低不僅受土壤養(yǎng)分和植物吸收能力的獨(dú)立作用，更取決于各因子間的協(xié)調(diào)程度[14]。目前鮮有選擇土壤肥力評價指標(biāo)劃分肥力等級以及確定權(quán)重統(tǒng)一的方法。以往研究采用各肥力因素得分之和表示土壤肥力高低，主觀人為因素對評價結(jié)果影響較大[28-29]。本研究采用模糊數(shù)學(xué)的隸屬度函數(shù)對土壤肥力進(jìn)行評價，利用相關(guān)系數(shù)法確定各評價指標(biāo)的權(quán)重，有效地減少人為因素對土壤肥力評價的影響，使評價結(jié)果更加客觀。本研究所用8 個土壤因子的權(quán)重系數(shù)在0.033～0.216 之間，其中有機(jī)質(zhì)含量所占權(quán)重較大，表明有機(jī)質(zhì)是影響土壤肥力質(zhì)量變異的主要因子，這與以往研究結(jié)果一致[29-30]；土壤肥力是土壤多項(xiàng)肥力因素綜合的結(jié)果，本研究花生種植前和成熟期土壤肥力綜合指標(biāo)(IFI)多集中于0.58～0.72 之間，與楊麗萍等[31]得出的結(jié)論相吻合，但不同于以往其他研究結(jié)果[14-15，18]，這可能是由研究區(qū)域和模型所選擇的指標(biāo)不同造成的。因此，為提高土壤肥力評價精度，未來需將更多的土壤物理、化學(xué)和生物因子加入到模糊數(shù)學(xué)隸屬度模型中；同時需要用作物產(chǎn)量對土壤肥力水平進(jìn)行檢驗(yàn)，以驗(yàn)證評價結(jié)果的準(zhǔn)確性。

塊金值與基臺值比值C0/(C＋C0)，表示隨機(jī)部分引起的空間異質(zhì)性占總變異的比例，一般認(rèn)為，若比值＜25%，表明系統(tǒng)空間相關(guān)性強(qiáng)烈，空間變異主要由結(jié)構(gòu)性因素起作用；若比值在25%～75%之間，表明系統(tǒng)的空間相關(guān)性中等，空間變異由結(jié)構(gòu)性因素和隨機(jī)因素共同起作用；若比值＞75%說明系統(tǒng)空間相關(guān)性較弱，空間變異主要由隨機(jī)因素起作用[32]。本研究中花生種植前IFI 的C0/(C＋C0)＜25%，具有較強(qiáng)的空間相關(guān)性，表明正陽縣花生田種植前土壤肥力的空間分布主要取決于結(jié)構(gòu)性因素(如土壤母質(zhì)、氣候、地形等自然因素)；花生成熟期IFI 的C0/(C＋C0)在25%～75%之間，空間相關(guān)性中等，說明正陽縣花生田成熟期土壤肥力空間分布由結(jié)構(gòu)性因素和隨機(jī)性因素(施肥、耕作管理等人為因素)共同作用的結(jié)果。另外，花生成熟期的平均IFI 高于花生種植前，其原因可能有以下幾點(diǎn):1)花生生產(chǎn)上農(nóng)戶的施肥量高于花生植株需求量；2)花生為豆科作物，可進(jìn)行根瘤固氮，增加了土壤氮含量；3)土壤溫度和濕度是影響土壤微生物活動和活性季節(jié)性變化的重要因素[33-34]，花生生長季為每年的6—9月份，該季節(jié)氣溫較高，雨水較充沛，土壤溫度和水分均適宜微生物的生長和繁殖，利于土壤中有機(jī)養(yǎng)分礦化。比較花生種植前和成熟期土壤肥力空間分布圖發(fā)現(xiàn)，正陽縣的土壤肥力南部高于北部，可能是由于正陽縣南鄰淮河，水熱條件利于土壤養(yǎng)分的積累。綜上，在未來花生生產(chǎn)上應(yīng)著重開展氮磷減施、生物固氮方面、以及水熱條件下土壤微生物的肥力效應(yīng)等方面的研究。

4 結(jié)論

本研究結(jié)果表明，與種植前相比，花生成熟期土壤pH 值和土壤有效硫含量出現(xiàn)下降，其余土壤養(yǎng)分指標(biāo)均有不同程度的提高，尤其是土壤全氮和有效磷含量增加明顯，但有機(jī)質(zhì)和有效鋅含量仍處于較低水平。從變異系數(shù)來看，8 個土壤肥力指標(biāo)中，pH 值屬弱變異，其余指標(biāo)均屬中等變異。從各肥力指標(biāo)權(quán)重來看，不同的土壤肥力指標(biāo)對IFI 的貢獻(xiàn)率存在差異，有機(jī)質(zhì)含量的權(quán)重最大，說明有機(jī)質(zhì)含量對土壤肥力的貢獻(xiàn)最大。半方差分析結(jié)果表明，花生種植前土壤肥力空間相關(guān)性較強(qiáng)，變異函數(shù)以球狀模型擬合效果最佳，空間分布特征主要受自然因素的影響；花生成熟期土壤肥力空間相關(guān)性中等，空間分布特征受自然因素和人為活動共同影響。此外，花生種植前IFI 低于花生成熟期，土壤肥力處于中等偏上(Ⅱ級和Ⅲ級)水平，IFI 空間分布特征呈南高北低。綜合分析，在正陽縣花生生產(chǎn)上，建議控制氮磷用量，增加有機(jī)肥和鋅肥的施用，尤其是在其南部地區(qū)。