汪宗蘭，王春霞，汪思佳，楊躍發(fā)，張景瑞

（1. 石河子大學(xué)水利建筑工程學(xué)院，新疆 石河子 832000； 2. 現(xiàn)代節(jié)水灌溉兵團(tuán)重點(diǎn)實(shí)驗(yàn)室，新疆 石河子 832000）

0 引 言

土壤肥力是土地生產(chǎn)力的基礎(chǔ)，是農(nóng)業(yè)可持續(xù)發(fā)展的基本要素［1，2］。隨著可持續(xù)農(nóng)業(yè)的提出與發(fā)展，現(xiàn)代農(nóng)業(yè)生產(chǎn)向著優(yōu)質(zhì)高產(chǎn)的方向轉(zhuǎn)變。長期大量元素化肥的施用，導(dǎo)致土壤養(yǎng)分失衡，因此科學(xué)、合理、實(shí)用地評(píng)價(jià)土壤肥力水平對(duì)指導(dǎo)農(nóng)民合理施肥和促進(jìn)農(nóng)業(yè)可持續(xù)發(fā)展有重要作用［3-5］。

目前，地統(tǒng)計(jì)學(xué)和地理空間信息技術(shù)的發(fā)展提供了更為精準(zhǔn)、便捷的土壤質(zhì)量評(píng)價(jià)方法。近些年來國內(nèi)外學(xué)者逐步將土壤質(zhì)量評(píng)價(jià)與空間信息技術(shù)相結(jié)合研究［6，7］，尤其對(duì)大中尺度土壤養(yǎng)分的空間變異特征及評(píng)價(jià)等研究較多［8-10］，一般而言大中尺度的土壤養(yǎng)分空間變異性分析對(duì)提高地力、建立土壤數(shù)據(jù)庫有重要意義［11］。而評(píng)價(jià)單元的大小直接影響到土壤質(zhì)量評(píng)價(jià)的精度，評(píng)價(jià)單元越小評(píng)價(jià)結(jié)果精度越高［12］。因此，國內(nèi)外學(xué)者對(duì)小尺度土壤肥力空間分布特征進(jìn)行了相關(guān)研究，李亞莉等［13］對(duì)新疆奇臺(tái)縣小麥土壤進(jìn)行分析，發(fā)現(xiàn)隨種植年限增加，土壤速效養(yǎng)分減少，空間分布受人為因素和隨機(jī)因素共同影響。展秀麗等［14］發(fā)現(xiàn)土壤有機(jī)碳和全氮在小尺度上表現(xiàn)出明顯的空間變異性，土壤機(jī)械組成對(duì)有機(jī)碳和全氮有明顯影響。連玉珍等［15］分析色季拉山西坡表層土壤的有機(jī)碳，結(jié)果表明小尺度下苔草高寒草甸和林芝杜鵑灌叢的土壤有機(jī)碳及其密度主要受自然因素的影響，具有較強(qiáng)的空間自相關(guān)性。鄒心雨等［16］對(duì)河北壩上地區(qū)土壤機(jī)械組成進(jìn)行研究，發(fā)現(xiàn)坡面小尺度土壤各粒級(jí)含量具有較強(qiáng)的空間自相關(guān)性，沙粒的空間變異性最強(qiáng)。喬學(xué)瑾等［17］對(duì)新疆安集海膜下滴灌棉田土壤進(jìn)行研究，發(fā)現(xiàn)土壤含鹽量受結(jié)構(gòu)性因素與人為活動(dòng)等隨機(jī)性因素共同影響呈強(qiáng)空間自相關(guān)性，采樣尺度越小，空間分布更加清晰明確。由此可見小尺度土壤養(yǎng)分具有明顯空間變異性，小尺度的研究更有利于對(duì)土壤實(shí)施科學(xué)合理的管理措施。

我國膜下滴灌技術(shù)的發(fā)展開始于1996 年在新疆生產(chǎn)建設(shè)兵團(tuán)第八師石河子121 團(tuán)試驗(yàn)的成功，并在此后推廣應(yīng)用于棉田生產(chǎn)。至2020 年新疆膜下滴灌面積突破128 萬hm2，成為我國大田應(yīng)用膜下滴灌技術(shù)面積最大的地區(qū)［18，19］。但由于滴灌“少灌勤灌”的特性，加上該區(qū)域蒸發(fā)強(qiáng)烈，土壤鹽分容易在土壤表層積聚，潛在的土壤次生鹽漬化風(fēng)險(xiǎn)增大，導(dǎo)致耕地土壤理化性質(zhì)退化、水肥保持能力降低。因此準(zhǔn)確、客觀的評(píng)價(jià)長期膜下滴灌農(nóng)田的肥力情況是防治膜下滴灌耕地土壤退化和提高生產(chǎn)力的關(guān)鍵措施之一。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于新疆石河子市石河子大學(xué)現(xiàn)代節(jié)水灌溉兵團(tuán)重點(diǎn)實(shí)驗(yàn)室基地（85°57′49＂E，44°19'28＂N）。該地區(qū)是典型的溫帶大陸性氣候。夏季相對(duì)濕度在30%～50%之間，年蒸發(fā)量為1 000～1 500 mm。年日照時(shí)長2 300～2 700 h。無霜期182 d。地下水埋深9 m 以下，灌溉水含鹽量約為1.35 g/L。田塊平整，種植作物主要為棉花和番茄，種植方式為膜下滴灌，滴灌年限為13 a。

1.2 樣品采集與分析

在2021 年10 月下旬整地前選擇膜下滴灌棉花和番茄地進(jìn)行土壤樣品采集。0～100 cm 土層平均干容重為1.62 g/cm3，土壤質(zhì)地為壤土，土壤砂粒、粉粒、黏粒含量分別為46.32%、38.53%、15.15%，土壤飽和含水量（質(zhì)量含水量）為26.11%，田間持水量為20.73%。該田塊種植作物、土壤質(zhì)地及種植模式在研究區(qū)具有代表性，采樣區(qū)域種植作物情況見圖1。按照20 m×20 m網(wǎng)格采樣方式采集表層（0～20 cm）土壤樣品88個(gè)。每個(gè)采樣點(diǎn)以定位點(diǎn)為中心在1 m 半徑內(nèi)采集土壤樣品4 個(gè)，混合均勻后采用四分法留存約1 kg土樣。樣品室內(nèi)自然風(fēng)干研磨過2 mm篩。采樣田塊取樣點(diǎn)位置分布詳見圖1。

圖1 研究區(qū)種植作物類型及采樣點(diǎn)分布圖Fig.1 Distribution of crop types and sampling points location

測(cè)定的指標(biāo)有土壤pH、電導(dǎo)率、有機(jī)質(zhì)、堿解氮、速效磷和速效鉀。指標(biāo)測(cè)定方法:pH 用電位法（土水比為5:1）所用pH計(jì)型號(hào)為哈希LA-pH10 型；電導(dǎo)率用電導(dǎo)率儀（土水比為5:1）所用電導(dǎo)率儀型號(hào)為梅特勒-托利多FE38 型；有機(jī)質(zhì)用重鉻酸鉀外加熱法；堿解氮用堿解擴(kuò)散法；速效磷用0.5 mol/L NaHCO3浸提鉬藍(lán)比色法；速效鉀用火焰光度法。

1.3 數(shù)據(jù)處理

1.3.1 權(quán)重計(jì)算

本研究采用相關(guān)系數(shù)法確定土壤各肥力指標(biāo)的權(quán)重。首先用統(tǒng)計(jì)分析軟件計(jì)算出各單項(xiàng)指標(biāo)間相關(guān)系數(shù)（表1），相關(guān)系數(shù)為負(fù)的取絕對(duì)值，進(jìn)而求出某單項(xiàng)指標(biāo)與其他指標(biāo)間相關(guān)系數(shù)的平均值，并將該平均值除以所有肥力指標(biāo)相關(guān)系數(shù)平均值總和的結(jié)果，作為該單項(xiàng)肥力指標(biāo)在土壤綜合肥力中的權(quán)重［20］（表2）。

表1 指標(biāo)間相關(guān)系數(shù)Tab.1 Correlation coefficients between indicators

表2 土壤肥力評(píng)價(jià)指標(biāo)權(quán)重Tab.2 Average correlation coefficients and index weights

1.3.2 隸屬度函數(shù)模型確定

確定隸屬度函數(shù)計(jì)算各個(gè)肥力指標(biāo)的隸屬度值。土壤電導(dǎo)率、有機(jī)質(zhì)、速效磷、速效鉀和堿解氮采用S型隸屬函數(shù)［公式（1）］，pH 采用拋物線型隸屬函數(shù)［公式（2）］［21］。根據(jù)研究區(qū)膜下滴灌土壤特征、種植作物及相關(guān)研究結(jié)果［22-24］，確定隸屬函數(shù)曲線轉(zhuǎn)折點(diǎn)的取值［24］，見表3。

表3 隸屬度函數(shù)轉(zhuǎn)折點(diǎn)取值Tab.3 Turning point values of evaluation indicators in the affiliation function curve

1.3.3 肥力評(píng)價(jià)指數(shù)計(jì)算

計(jì)算土壤綜合肥力指數(shù)（IFI），具體計(jì)算公式如下:

式中:n為參評(píng)因子數(shù)；Wi為第i個(gè)因子的權(quán)重；Fi為第i個(gè)因子的隸屬度值。IFI的結(jié)果為0-1，其值越接近1 代表土壤質(zhì)量越好。

1.3.4 土壤養(yǎng)分綜合肥力分級(jí)

將膜下滴灌土壤IFI值根據(jù)等距法分為5 個(gè)等級(jí)［25］，具體見表4。

表4 土壤綜合肥力指數(shù)分級(jí)Tab.4 Soil integrated fertility index classification

1.3.5 空間變異性分析

采用地統(tǒng)計(jì)學(xué)分析各項(xiàng)指標(biāo)的空間變異性其半變異函數(shù)表達(dá)式（4）［20］:

式中:r(h)為半變異函數(shù)；h為樣本點(diǎn)間距；N(h)以h為間距時(shí)點(diǎn)的成對(duì)數(shù)；Z(xi)和Z(xi+h)為區(qū)域化變量在xi與xi+h處的實(shí)測(cè)值。在GS+7.0 中完成各土壤指標(biāo)半方差函數(shù)理論模型的構(gòu)建與擬合。

2 結(jié)果與分析

2.1 土壤各項(xiàng)肥力指標(biāo)特征分析

2.1.1 各項(xiàng)肥力指標(biāo)描述性統(tǒng)計(jì)分析

在SPSS 25.0中對(duì)各項(xiàng)肥力指標(biāo)進(jìn)行描述性統(tǒng)計(jì)分析，從表5 可以看出，pH 的均值為8.43，土壤呈堿性。有機(jī)質(zhì)、堿解氮和速效磷的均值分別為22.32、108.46 和19.61 g/kg 含量適中。速效鉀均值為423.25 g/kg，含量處于較高水平。土壤各肥力指標(biāo)經(jīng)過單一樣本Kolmogorov-Smirnov（K-S）正態(tài)分布檢驗(yàn)，電導(dǎo)率的p值小于0.05為非正態(tài)分布，對(duì)其進(jìn)行對(duì)數(shù)轉(zhuǎn)化后，結(jié)果符合正態(tài)分布特征。pH、有機(jī)質(zhì)、堿解氮、速效磷和速效鉀的p值分別為0.20、0.20、0.07、0.06和0.08，均大于0.05，為正態(tài)分布。pH的變異系數(shù)最小，僅為0.02，表現(xiàn)為弱變異。其他肥力指標(biāo)的變異系數(shù)在0.13～0.40之間，屬于中等變異強(qiáng)度。各肥力指標(biāo)極差有大（比如電導(dǎo)率）有?。ū热鏿H），可以看出田間尺度上各點(diǎn)間土壤養(yǎng)分含量差異明顯。有機(jī)質(zhì)和pH 的偏度系數(shù)小于零，其分布峰為左偏。電導(dǎo)率、堿解氮、速效磷和速效鉀的偏度值大于零，分布峰為右偏。有機(jī)質(zhì)、堿解氮、速效磷和速效鉀的峰度系數(shù)為正，說明這些指標(biāo)的分布曲線相較于正態(tài)分布更加陡峭，反之pH 峰度系數(shù)為負(fù)，說明該組數(shù)據(jù)的分布曲線相較于正態(tài)分布更加平緩。

表5 各項(xiàng)肥力指標(biāo)含量描述性統(tǒng)計(jì)Tab.5 Descriptive statistics of the content of each index

2.1.2 各項(xiàng)肥力指標(biāo)空間變異特征分析

在GS+7.0中對(duì)各項(xiàng)指標(biāo)進(jìn)行地統(tǒng)計(jì)學(xué)分析，確定最適的半方差函數(shù)模型，確定原則是殘差平方和最小、決定系數(shù)最大。pH 和電導(dǎo)率符合高斯模型，有機(jī)質(zhì)符合線性模型，速效磷和速效鉀符合球形模型，堿解氮符合指數(shù)模型。除堿解氮的擬合系數(shù)為0.13 外，有機(jī)質(zhì)和速效磷擬合效果最好，決定系數(shù)分別為0.88 和0.86，pH、電導(dǎo)率和速效鉀擬合效果較好，擬合系數(shù)分別為0.75、0.62 和0.65。膜下滴灌土壤各項(xiàng)肥力指標(biāo)半方差函數(shù)模型見表6。

表6 土壤各項(xiàng)肥力指標(biāo)半方差函數(shù)模型及參數(shù)Tab.6 Semi-covariance function model and parameters for each soil fertility index

田塊尺度上土壤pH、電導(dǎo)率、堿解氮和速效磷的塊金系數(shù)大于75%分別為97.6%、98.5%、97.2%和99.8%，表現(xiàn)出較弱的空間自相關(guān)性，表明這些肥力指標(biāo)主要受人為因素的影響，即施肥、灌溉和管理措施等人為因素對(duì)上述指標(biāo)的影響較大。有機(jī)質(zhì)和速效鉀指標(biāo)的塊金系數(shù)在25%～75%之間，分別為51.2%和62.5%，表明二者具有中等強(qiáng)度的空間自相關(guān)性，空間變異受隨機(jī)因素和結(jié)構(gòu)性因素共同影響，并且兩者的塊基比都大于0.5，說明更多的受隨機(jī)因素（人為因素）影響。各項(xiàng)肥力指標(biāo)的變程在23.89～120.20 m 之間，說明各項(xiàng)肥力指標(biāo)的空間自相關(guān)距離存在較大差異。原因在于“隨水滴肥”的灌溉方式及種植作物不同。長期隨水滴施肥料，肥料大部分僅在土壤濕潤體內(nèi)分布，所以田塊內(nèi)肥料空間分布不均勻，使土壤肥力指標(biāo)空間變異性增強(qiáng)。各項(xiàng)土壤肥力指標(biāo)的變程均大于取樣間距，所以在以后分析膜下滴灌土壤肥力狀況時(shí)可以適當(dāng)?shù)脑黾硬蓸娱g距。

在Surfer 15.0 中繪制各項(xiàng)肥力指標(biāo)空間分布圖（圖2），圖2中可以看出:①pH 范圍在8.0～9.1 之間，全區(qū)土壤偏堿性，在空間上呈帶狀分布特征。整體東側(cè)pH 值要高于西側(cè)，其中北部區(qū)域pH 值達(dá)到8.5 以上，堿性程度強(qiáng)烈，東北部的防護(hù)林土壤堿性程度較高。②電導(dǎo)率值成斑塊狀分布且含量差異巨大，北部電導(dǎo)率值低于南部，尤其是東北側(cè)防護(hù)林所在區(qū)域明顯偏低，其在研究區(qū)南部有島狀高含量點(diǎn)存在，原因在于此區(qū)域有過微咸水灌溉試驗(yàn)處理。③有機(jī)質(zhì)田塊尺度上分布較為均勻，含量在20～30 g/kg 之間，在肥力等級(jí)劃分［26］（表7）中處于中等水平，在研究區(qū)東北部有島狀高點(diǎn)存在，原因可能是胡楊樹落葉長期積累導(dǎo)致該區(qū)域有機(jī)質(zhì)含量偏高。④堿解氮含量斑塊狀分布明顯，呈現(xiàn)四周高中間低，北部高南部低的分布格局，同時(shí)在研究區(qū)東北部有島狀高點(diǎn)存在，原因在于枯枝落葉等凋落物的分解顯著提高表層土壤堿解氮的含量。⑤速效鉀含量普遍較高，根據(jù)表7整個(gè)區(qū)域最低值遠(yuǎn)大于豐富水平，且其含量在研究區(qū)分布不均勻，東部番茄地速效鉀含量普遍高于西部棉花地，東北部胡楊林種植區(qū)有島狀高點(diǎn)存在。⑥速效磷含量分布呈現(xiàn)東高西低南高北低的特點(diǎn)，在研究區(qū)北部和南部有兩個(gè)集中高點(diǎn)，其余區(qū)域分布較為均勻?？傮w來說，研究區(qū)土壤各指標(biāo)在田塊尺度上分布復(fù)雜、差異性顯著，說明種植制度、灌溉和作物等因素的影響，導(dǎo)致膜下滴灌土壤養(yǎng)分在田間尺度上存在明顯的差異性。

表7 土壤養(yǎng)分分級(jí)標(biāo)準(zhǔn)Tab.7 Soil nutrient grading criteria

圖2 土壤養(yǎng)分空間分布圖Fig.2 Soil nutrient spatial distribution map

2.2 土壤綜合肥力指數(shù)（IFI）特征分析

2.2.1IFI描述性統(tǒng)計(jì)分析

對(duì)88個(gè)樣點(diǎn)的土壤綜合肥力指數(shù)進(jìn)行描述性統(tǒng)計(jì)分析（表8），研究區(qū)IFI值在0.51～0.74 之間，平均值為0.62，變異系數(shù)為0.08，屬于弱變異。IFI值的偏度系數(shù)大于零，表明分布峰為左偏。IFI值采用單一樣本K-S 檢驗(yàn)法（表8），其p值為0.17（大于0.05）為正態(tài)分布。將各項(xiàng)土壤各項(xiàng)肥力指標(biāo)的平均隸屬度值繪制雷達(dá)圖（圖3）。根據(jù)雷達(dá)圖的幾何含義，坐標(biāo)軸上的點(diǎn)越靠近原點(diǎn)，代表它所反映的指標(biāo)的肥力水平越低，反之越高。可以看出pH、速效磷和電導(dǎo)率平均隸屬度值最小，約為0.05；有機(jī)質(zhì)次之，約為0.15；速效鉀和堿解氮約為0.2。

表8 IFI值描述性統(tǒng)計(jì)Tab.8 Descriptive statistics of IFI values

圖3 各項(xiàng)肥力指標(biāo)平均隸屬度雷達(dá)圖Fig.3 Radar map of the average affiliation of each fertility index

2.2.2IFI空間分布特征分析

由圖4可知，該研究區(qū)膜下滴灌土壤全被Ⅱ類（良好肥力水平）和Ⅲ類（中等肥力水平）肥力（見表4）占據(jù)，其中Ⅱ類（良好肥力水平）地分布最廣，主要分布在研究區(qū)中部，占整個(gè)區(qū)域面積71.11%。田塊南部和北部以Ⅲ類（中等肥力水平）地為主，占整個(gè)區(qū)域面積28.89%。因此，按照模糊綜合評(píng)價(jià)法，研究區(qū)土壤的肥力水平整體處于中等偏上水平。研究區(qū)北部因接近地塊邊緣，肥力水平低于其他區(qū)域；中部區(qū)域以種植棉花和番茄為主，農(nóng)民注重投入和管理，故總體肥力水平偏高。南部地區(qū)由于局部地塊做過微咸水灌溉處理因此土壤肥力偏低。

圖4 土壤IFI分布圖Fig.4 Soil IFI distribution

3 討 論

土壤是農(nóng)業(yè)生產(chǎn)的最基本資料，土壤肥力水平的高低直接影響到農(nóng)作物的生產(chǎn)效益。我國精準(zhǔn)農(nóng)業(yè)快速發(fā)展，作為其中的一部分，小尺度土壤養(yǎng)分空間變異的研究十分重要。進(jìn)行土壤養(yǎng)分綜合評(píng)價(jià)有助于了解研究區(qū)土地利用現(xiàn)狀，明確土壤生產(chǎn)力限制因子，為土壤質(zhì)量改善以及土地可持續(xù)利用提供基礎(chǔ)數(shù)據(jù)支持［27，28］。

描述性統(tǒng)計(jì)分析表明，土壤pH 和堿解氮、速效鉀、速效磷含量呈明顯的負(fù)相關(guān)，說明適當(dāng)降低土壤pH 值，可以提高化肥的利用率；電導(dǎo)率的極差很大，pH的極差很小，這表明土壤屬性在田間尺度存在明顯差異。pH的變異系數(shù)最小為0.02，可能原因在于研究區(qū)范圍小，pH 值改變不大，也可以看出土壤本身的屬性具有更好的穩(wěn)定特征這與劉月嬌等［29］的研究結(jié)果一致。有機(jī)質(zhì)、堿解氮、速效磷、速效鉀均屬于中等變異強(qiáng)度，主要因?yàn)樗傩юB(yǎng)分受施肥和耕作方式、作物的種類的影響較大陶睿等［27］在對(duì)重慶豐都紫色土壤的評(píng)價(jià)也得出了相似結(jié)論。

地統(tǒng)計(jì)學(xué)分析表明，在田塊尺度上各項(xiàng)肥力指標(biāo)的空間自相關(guān)距離有大（如有機(jī)質(zhì)）有?。ㄈ鐗A解氮），說明有機(jī)質(zhì)受隨機(jī)因素影響較小，空間自相關(guān)的范圍較大，堿解氮受施肥等人為因素影響較大，空間自相關(guān)距離較小。土壤各肥力指標(biāo)在田塊尺度上斑塊狀分布特征明顯，是因?yàn)椤半S水滴肥”的灌溉方式，肥料大部分僅在土壤濕潤體內(nèi)分布，所以田塊內(nèi)肥力指標(biāo)含量空間分布不均勻，受隨機(jī)因素的影響，使土壤肥力指標(biāo)空間分布的差異性增大，周英霞［30］對(duì)棉花生長以及農(nóng)田尺度上土壤水肥鹽時(shí)空分布特征研究中也得出了相似結(jié)論。

根據(jù)劃分的土壤肥力綜合指標(biāo)值（IFI）等級(jí)，土壤肥力特征在空間上呈現(xiàn)出斑塊狀。主要原因在于滴灌使局部土壤濕潤，養(yǎng)分則在作物耕作層滴灌帶下分布較高，周圍裸地上分布較少，由此使肥料在土壤中缺乏連續(xù)性這與李亞莉等［13］的研究結(jié)果類似。同時(shí)田塊尺度上膜下滴灌土壤肥力整體較高，也說明膜下滴灌技術(shù)的水肥一體化實(shí)施可以提高作物的肥料利用效率［31］。通過平均隸屬度雷達(dá)圖可知研究區(qū)土壤速效磷含量低，是土壤肥力的主要制約因素。

本研究選取了滴灌13 年田塊土壤耕層的6 種化學(xué)指標(biāo)分析土壤肥力狀況，但缺少生物等指標(biāo)，需要在后續(xù)研究中完善。此外本研究所選取得膜下滴灌田塊種植作物、滴灌的應(yīng)用年限等都具有代表性，因此所得出的土壤各肥力指標(biāo)的評(píng)價(jià)結(jié)果對(duì)于類似膜下滴灌土壤具有一定的參考作用。

4 結(jié) 論

利用模糊評(píng)判法和地統(tǒng)計(jì)學(xué)方法對(duì)長期膜下滴灌（13 年）農(nóng)田土壤養(yǎng)分空間分布特征及肥力進(jìn)行綜合評(píng)價(jià)。主要結(jié)論如下:

（1）6 項(xiàng)肥力指標(biāo)的經(jīng)典統(tǒng)計(jì)表明，研究區(qū)土壤pH、電導(dǎo)率、有機(jī)質(zhì)、堿解氮、速效磷、速效鉀的平均含量為8.43、5 726.57 μS/cm、22.32 g/kg、108.46 mg/kg、19.61 mg/kg、423.25 mg/kg。pH 的變異系數(shù)為0.02，屬于弱變異，均值為8.43，土壤呈堿性。其他肥力指標(biāo)變異系數(shù)在0.13～0.40 之間，屬于中等變異。

（2）土壤pH、電導(dǎo)率、堿解氮和速效磷指標(biāo)的變程分別為32.9、31.0、23.9、36.7m且它們的塊金系數(shù)大于75%，表明在空間上呈現(xiàn)弱空間自相關(guān)性。有機(jī)質(zhì)和速效鉀指標(biāo)的變程分別為120.2和45.1 m，它的塊金系數(shù)位于25%～75%之間，呈現(xiàn)出中等強(qiáng)度的空間自相關(guān)性。

（3）各土壤屬性在田塊尺度上空間分布復(fù)雜土壤pH 呈帶狀分布特征；電導(dǎo)率斑塊狀分布，且含量差異很大；有機(jī)質(zhì)分布較為均勻；堿解氮呈現(xiàn)北高南低、四周高中間低的趨勢(shì)；速效鉀分布不均勻；速效磷東高西低、南高北低，南北部有兩個(gè)集中高點(diǎn)。

（4）研究區(qū)土壤肥力綜合指標(biāo)值（IFI）值在0.51～0.74 之間，平均值為0.62。研究區(qū)土壤全被Ⅱ、Ⅲ類土壤占據(jù)，占比分別為71.11%和28.89%，土壤肥力質(zhì)量處于中上水平，但有效磷含量偏低應(yīng)該通過增施磷肥方式改良土壤。