宿州市耕地土壤養(yǎng)分時(shí)空變化特征分析

丁琪洵，詹雪潔，張?zhí)於?，許諾，馬秀婷，張長(zhǎng)坤，馬友華*

（1.農(nóng)田生態(tài)保育與污染防控安徽省重點(diǎn)實(shí)驗(yàn)室，安徽農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，合肥 230036；2.宿州市農(nóng)業(yè)農(nóng)村局，安徽 宿州 234000；3.安徽華成種業(yè)股份有限公司，安徽 宿州 234000）

土壤有機(jī)質(zhì)、有效磷、速效鉀等土壤養(yǎng)分含量的豐缺程度決定著土壤肥力的高低，直接影響農(nóng)作物的生產(chǎn)和產(chǎn)量，因此深入分析土壤養(yǎng)分的時(shí)空變化規(guī)律及其產(chǎn)生的原因具有重要意義。了解掌握土壤養(yǎng)分時(shí)空變異性的驅(qū)動(dòng)因素能夠監(jiān)測(cè)土壤肥力變化、指導(dǎo)各地區(qū)改進(jìn)施肥技術(shù)、開(kāi)展科學(xué)有效養(yǎng)分管理，為構(gòu)建生態(tài)高值農(nóng)業(yè)產(chǎn)業(yè)體系奠定基礎(chǔ)。近年來(lái)，國(guó)內(nèi)外對(duì)不同尺度下耕地土壤養(yǎng)分時(shí)空變化開(kāi)展了大量研究，研究尺度主要包括省級(jí)、市級(jí)、縣域、流域、丘陵區(qū)等尺度[1-5]；種植作物類(lèi)型包括煙草、芒果、小麥、茶葉等[6-8]；影響土壤養(yǎng)分時(shí)空變化特征的因素主要包括氣候因子、坡向、土壤理化性質(zhì)、地貌等自然因素，以及土地利用、施肥、種植結(jié)構(gòu)、采樣間距等人為因素[9-13]。但目前耕地養(yǎng)分變化研究系統(tǒng)性不夠，時(shí)間尺度上長(zhǎng)期性尚欠缺。本研究基于宿州市2009 年測(cè)土配方施肥項(xiàng)目、2019 年耕地質(zhì)量等級(jí)評(píng)價(jià)項(xiàng)目中耕地質(zhì)量長(zhǎng)期監(jiān)測(cè)點(diǎn)位以及布點(diǎn)均勻合理能夠代表典型地塊的采樣點(diǎn)的養(yǎng)分?jǐn)?shù)據(jù)，建立宿州市及各縣區(qū)的耕地土壤養(yǎng)分評(píng)價(jià)數(shù)據(jù)庫(kù)，利用地統(tǒng)計(jì)學(xué)結(jié)合GIS 反距離權(quán)重空間插值的方法，探究宿州市耕地土壤有機(jī)質(zhì)、全氮、有效磷和速效鉀的時(shí)空分布特征和養(yǎng)分變化影響因素。研究結(jié)果可為宿州市土壤養(yǎng)分分區(qū)管理和農(nóng)戶(hù)精準(zhǔn)施肥提供數(shù)據(jù)支持。

1 材料與方法

1.1 研究區(qū)概況

宿州市屬于我國(guó)糧食生產(chǎn)核心區(qū)的黃淮海區(qū)，是重要的農(nóng)產(chǎn)品生產(chǎn)供應(yīng)基地，位于安徽省北部、黃淮海平原南端，東至東北與江蘇省宿遷市、徐州市接壤，南臨蚌埠市，西至西北與淮北市、河南省商丘市、山東省菏澤市毗鄰，在東經(jīng)116°09′～118°10′、北緯33°18′～34°38′之間，地勢(shì)自西北向東南傾斜，主要地形地貌為丘陵、階地、平原。宿州市年降雨量在630.3～1 040.8 mm 之間，屬暖溫帶半濕潤(rùn)季風(fēng)氣候區(qū)，轄碭山縣、蕭縣、靈璧縣、泗縣和埇橋區(qū)，主要農(nóng)作物有小麥、玉米、大豆、蔬菜等。2019 年宿州市第二次國(guó)土調(diào)查顯示總耕地面積為575 034.50 hm2，宿州市統(tǒng)計(jì)年鑒數(shù)據(jù)顯示，2020 年宿州市糧食作物播種面積約94 萬(wàn)hm2，總產(chǎn)量516.55 萬(wàn)t。宿州市主要土壤類(lèi)型包括潮土和砂姜黑土，二者占宿州市耕地總面積的97.31%。其中：潮土占耕地總面積的56.08%，主要分布在宿州市北部的碭山縣和蕭縣全縣，以及埇橋區(qū)、靈璧縣、泗縣的北部地區(qū)；砂姜黑土占耕地總面積的41.23%，主要分布在宿州市埇橋區(qū)、靈璧縣、泗縣的南部地區(qū)。此外，黃褐土、石灰（巖）土、水稻土在解集鄉(xiāng)、莊里鄉(xiāng)、官橋鎮(zhèn)等零星分布。

1.2 數(shù)據(jù)來(lái)源

數(shù)據(jù)主要包括比例尺1∶50 000的宿州市土壤圖、宿州市行政區(qū)劃圖、第二次國(guó)土調(diào)查土地利用現(xiàn)狀圖、宿州市統(tǒng)計(jì)年鑒、2005 年以來(lái)的測(cè)土配方施肥數(shù)據(jù)庫(kù)及圖件、宿州市耕地質(zhì)量長(zhǎng)期定位監(jiān)測(cè)點(diǎn)點(diǎn)位信息、2019 年宿州市耕地質(zhì)量等級(jí)評(píng)價(jià)數(shù)據(jù)庫(kù)。宿州市測(cè)土配方采樣期為2005—2009 年，于2010 年完成耕地地力評(píng)價(jià)并建立測(cè)土配方施肥數(shù)據(jù)庫(kù)，涵蓋宿州市耕地5 297 個(gè)采樣點(diǎn)的基礎(chǔ)養(yǎng)分?jǐn)?shù)據(jù)。宿州市2019 年完成耕地質(zhì)量等級(jí)評(píng)價(jià)數(shù)據(jù)庫(kù)建立，涵蓋1 577 個(gè)采樣點(diǎn)的基礎(chǔ)養(yǎng)分?jǐn)?shù)據(jù)，原始數(shù)據(jù)均由宿州市農(nóng)業(yè)農(nóng)村局提供。宿州市2010年和2019年的土壤養(yǎng)分調(diào)查點(diǎn)位分布如圖1所示。

圖1 宿州市不同時(shí)期土壤養(yǎng)分調(diào)查點(diǎn)位分布圖Figure 1 Distribution map of soil nutrient survey points in different periods of Suzhou

1.3 研究方法

將宿州市土壤圖、第二次國(guó)土調(diào)查土地利用現(xiàn)狀圖、行政區(qū)劃圖三圖疊加形成宿州市耕地土壤養(yǎng)分指標(biāo)評(píng)價(jià)單元，即相同土壤單元、土地利用現(xiàn)狀類(lèi)型和行政區(qū)劃單位的地塊組成一個(gè)評(píng)價(jià)單元。再利用ArcGIS 軟件進(jìn)行反距離權(quán)重插值，將采樣點(diǎn)的養(yǎng)分?jǐn)?shù)據(jù)通過(guò)以點(diǎn)代面的方式在耕地土壤養(yǎng)分評(píng)價(jià)單元上進(jìn)行賦值，形成不同土壤養(yǎng)分的分布圖。

1.4 耕地養(yǎng)分分級(jí)標(biāo)準(zhǔn)

參照農(nóng)業(yè)農(nóng)村部耕地質(zhì)量監(jiān)測(cè)保護(hù)中心發(fā)布的耕地監(jiān)測(cè)函〔2019〕30號(hào)文件、安徽省耕地質(zhì)量監(jiān)測(cè)指標(biāo)分級(jí)標(biāo)準(zhǔn)[14]，將耕地土壤有機(jī)質(zhì)、全氮、有效磷和速效鉀含量分為五級(jí)，具體見(jiàn)表1。

表1 安徽省耕地質(zhì)量監(jiān)測(cè)指標(biāo)養(yǎng)分分級(jí)標(biāo)準(zhǔn)Table 1 Nutrient classification standard of cultivated land quality monitoring index in Anhui Province

2 結(jié)果與分析

2.1 宿州市耕地土壤有機(jī)質(zhì)時(shí)空變化特征

有機(jī)質(zhì)是評(píng)價(jià)耕地質(zhì)量的重要養(yǎng)分指標(biāo)。從土壤有機(jī)質(zhì)含量的時(shí)空變化情況來(lái)看，2019 年宿州市耕地土壤有機(jī)質(zhì)平均值為17.95 g·kg-1，比2010 年平均值（16.91 g·kg-1）提升了1.04 g·kg-1。2010 年宿州市各縣區(qū)有機(jī)質(zhì)含量均值由高到低依次為蕭縣＞泗縣＞埇橋區(qū)＞靈璧縣＞碭山縣，2019 年由高到低依次為埇橋區(qū)＞泗縣＞蕭縣＞靈璧縣＞碭山縣。9 年來(lái)宿州市碭山縣、靈璧縣和埇橋區(qū)的耕地土壤有機(jī)質(zhì)均值分別提升0.66、1.12 g·kg-1和4.13 g·kg-1，蕭縣和泗縣有機(jī)質(zhì)均值分別降低1.60 g·kg-1和0.26 g·kg-1；不同土類(lèi)的有機(jī)質(zhì)含量均值呈砂姜黑土顯著提升、潮土和黃褐土略微提升、石灰?guī)r土和水稻土略微下降的變化規(guī)律（表2、表3）。

表2 2010年和2019年宿州市各縣區(qū)耕地養(yǎng)分含量（均值±標(biāo)準(zhǔn)差）統(tǒng)計(jì)Table 2 The nutrient content（Mean±SD）of cultivated land of different county in Suzhou in 2010 and 2019

表3 2010年和2019年宿州市各土類(lèi)耕地養(yǎng)分含量（均值±標(biāo)準(zhǔn)差）統(tǒng)計(jì)Table 3 The nutrient content（Mean±SD）of cultivated land of different soil category in Suzhou in 2010 and 2019

從空間分布特征（圖2）看，2019 年宿州市耕地土壤有機(jī)質(zhì)含量主要處于中等水平，面積為437 005.35 hm2，占宿州市總耕地面積的76.00%（表4）。碭山縣有機(jī)質(zhì)含量處于較低級(jí)別（Ⅳ級(jí)）的耕地面積占比最大，面積為24 646.98 hm2，靈璧縣、泗縣、蕭縣和埇橋區(qū)有機(jī)質(zhì)含量處于中等級(jí)別（Ⅲ級(jí)）的耕地面積占比最大。從圖2 中耕地土壤有機(jī)質(zhì)空間分布可知，有機(jī)質(zhì)含量表現(xiàn)為：碭山縣整體較缺，蕭縣酒店鄉(xiāng)較缺，埇橋區(qū)朱仙莊、灰古等中部鄉(xiāng)鎮(zhèn)較豐富，靈璧縣含量整體中等，泗縣大路口、墩集等個(gè)別鄉(xiāng)鎮(zhèn)較缺。

2.2 宿州市耕地土壤全氮時(shí)空變化特征

從土壤全氮含量的時(shí)空變化情況來(lái)看，2019 年宿州市耕地土壤全氮平均值為1.06g·kg-1，與2010 年全氮的平均值（1.07 g·kg-1）基本持平。9 年來(lái)宿州市靈璧縣、泗縣和埇橋區(qū)的耕地土壤全氮含量平均值基本穩(wěn)定，碭山縣全氮含量均值提升0.17 g·kg-1，蕭縣全氮含量均值略微降低0.10 g·kg-1，同時(shí)宿州市各土類(lèi)全氮含量平均值也基本持平（表2、表3）。2019 年宿州市全氮多為中等水平，面積為349 618.23 hm2，占宿州市耕地總面積的60.80%。2019 年宿州市耕地土壤全氮含量空間分布特征為：碭山縣全氮處于較低級(jí)別（Ⅳ級(jí)）的耕地面積占比最大，為52.74%；其余縣區(qū)全氮處于中等級(jí)別（Ⅲ級(jí)）的耕地面積占比最大（表4）。比較圖2 中耕地土壤全氮含量分布可知，全氮含量表現(xiàn)為：碭山縣整體較缺，埇橋區(qū)欄桿鎮(zhèn)等北部鄉(xiāng)鎮(zhèn)缺乏，靈璧縣澮溝、尹集鎮(zhèn)等中部鄉(xiāng)鎮(zhèn)較豐，蕭縣酒店鄉(xiāng)、趙莊鎮(zhèn)等西部鄉(xiāng)鎮(zhèn)較缺，泗縣丁湖鎮(zhèn)、大路口鄉(xiāng)等南部鄉(xiāng)鎮(zhèn)較缺。

圖2 宿州市不同時(shí)期耕地土壤養(yǎng)分含量分布圖Figure 2 Distribution map of soil nutrient contents of cultivated land in different periods in Suzhou City

續(xù)圖2 宿州市不同時(shí)期耕地土壤養(yǎng)分含量分布圖Continued figure 2 Distribution map of soil nutrient contents of cultivated land in different periods in Suzhou City

2.3 宿州市耕地土壤有效磷時(shí)空變化特征

土壤有效磷含量是衡量土壤養(yǎng)分容量和強(qiáng)度水平的重要指標(biāo)，耕地土壤中的磷一般以無(wú)機(jī)磷和有機(jī)磷兩種形態(tài)存在。從土壤有效磷含量的時(shí)空變化情況來(lái)看，2019 年耕地土壤有效磷含量平均為16.44 mg·kg-1，比2010 年（14.84 mg·kg-1）提高了1.60 mg·kg-1。2010 年宿州市各縣區(qū)有效磷含量均值由高到低依次為碭山縣＞蕭縣＞埇橋區(qū)＞靈璧縣＞泗縣，2019年由高到低依次為埇橋區(qū)＞碭山縣＞靈璧縣＞蕭縣＞泗縣。9 年來(lái)靈璧縣、泗縣和埇橋區(qū)的耕地土壤有效磷平均值分別提升0.50、1.24 mg·kg-1和5.49 mg·kg-1，碭山縣和蕭縣耕地土壤有效磷均值分別降低0.49 mg·kg-1和1.72 mg·kg-1；各土類(lèi)呈潮土和石灰?guī)r土有效磷含量基本持平，黃褐土、砂姜黑土和水稻土有效磷含量顯著提升的變化規(guī)律（表2、表3）。

從空間分布特征（圖2）看，2019 年宿州市耕地土壤有效磷含量主要處于中等水平，面積為346 931.25 hm2，占宿州市耕地總面積的60.33%。碭山縣、靈璧縣、蕭縣和埇橋區(qū)有效磷含量處于中等級(jí)別（Ⅲ級(jí)）的耕地面積占比最大，泗縣有效磷含量處于較低級(jí)別（Ⅳ級(jí)）的耕地面積占比最大，面積為55 503.99 hm2（表4）。比較圖2 中耕地土壤有效磷含量分布可知，有效磷含量表現(xiàn)為：碭山縣朱樓、程莊等南部鄉(xiāng)鎮(zhèn)較豐，埇橋區(qū)楊莊鄉(xiāng)、欄桿鎮(zhèn)等北部鄉(xiāng)鎮(zhèn)較缺，埇橋區(qū)大營(yíng)、永鎮(zhèn)等南部鄉(xiāng)鎮(zhèn)較豐，蕭縣整體含量中等，泗縣山頭、黃圩等北部鄉(xiāng)鎮(zhèn)以及草溝、丁湖等南部鄉(xiāng)鎮(zhèn)較缺。

2.4 宿州市耕地土壤速效鉀時(shí)空變化特征

土壤速效鉀指水溶性鉀和黏土礦物晶體外表面吸附的交換性鉀，可以被植物直接吸收利用，在植物生長(zhǎng)及其品質(zhì)方面具有重要作用。速效鉀含量水平反映了土壤的供鉀能力，是土壤質(zhì)量的一項(xiàng)重要指標(biāo)。從土壤速效鉀含量的時(shí)空變化情況來(lái)看，2019年宿州市耕地土壤速效鉀含量平均為149.37 mg·kg-1，比2010 年均值（135.99 mg·kg-1）顯著提高13.38 mg·kg-1。泗縣耕地土壤速效鉀含量最高，均值為178.75 mg·kg-1，碭山縣和靈璧縣耕地土壤速效鉀含量較低，均值分別為120.36 mg·kg-1和129.23 mg·kg-1。2010 年宿州市各縣區(qū)速效鉀含量均值由高到低依次為泗縣＞蕭縣＞靈璧縣＞碭山縣＞埇橋區(qū)，2019年由高到低依次為泗縣＞埇橋區(qū)＞蕭縣＞靈璧縣＞碭山縣。9 年來(lái)泗縣和埇橋區(qū)耕地土壤速效鉀均值分別顯著提升25.94 mg·kg-1和31.94 mg·kg-1，靈璧縣和蕭縣耕地土壤速效鉀均值分別降低8.29 mg·kg-1和5.81 mg·kg-1。各土類(lèi)呈潮土、黃褐土、砂姜黑土和水稻土速效鉀含量顯著提升，石灰?guī)r土速效鉀含量略微下降的變化規(guī)律（表2、表3）。

從空間分布特征（圖2）看，2019年宿州市耕地土壤速效鉀含量多處于中等水平，面積為271 995.75 hm2，占宿州市耕地總面積的47.30%。碭山縣、靈璧縣、蕭縣和埇橋區(qū)耕地土壤速效鉀含量處于中等級(jí)別（Ⅲ級(jí)）的耕地面積占比最大，泗縣耕地土壤速效鉀含量處于較高級(jí)別（Ⅱ級(jí)）的耕地面積占比最大，面積為47 884.29 hm2（表4）。比較圖2中耕地土壤速效鉀含量分布可知，速效鉀含量表現(xiàn)為：碭山縣中等，蕭縣張莊寨、大屯等西南部鄉(xiāng)鎮(zhèn)較豐，埇橋區(qū)永安、桃溝等中部鄉(xiāng)鎮(zhèn)含量豐富，楊莊、褚蘭等北部鄉(xiāng)鎮(zhèn)較缺，靈璧縣尹集、澮溝等西部鄉(xiāng)鎮(zhèn)豐富，泗縣瓦坊、劉圩等北部鄉(xiāng)鎮(zhèn)豐富。

表4 2019年宿州市耕地質(zhì)量養(yǎng)分含量分級(jí)及占比Table 4 Classification and proportion of cultivated land quality nutrient content in Suzhou City in 2019

2.5 耕地土壤養(yǎng)分指標(biāo)相關(guān)性分析

通過(guò)SPSS計(jì)算皮爾遜相關(guān)系數(shù)，進(jìn)而分析土壤養(yǎng)分指標(biāo)間的相關(guān)性，土壤有機(jī)質(zhì)、全氮、有效磷、速效鉀4個(gè)指標(biāo)之間均存在0.01水平的極顯著相關(guān)關(guān)系。其中有機(jī)質(zhì)與全氮的正相關(guān)性最強(qiáng)，2010 年相關(guān)系數(shù)為0.607，2019年相關(guān)系數(shù)為0.722；有機(jī)質(zhì)與速效鉀的正相關(guān)性次之，2010年相關(guān)系數(shù)達(dá)到0.455，2019年相關(guān)系數(shù)為0.432。2010年土壤有效磷與其他三項(xiàng)指標(biāo)呈負(fù)相關(guān)，相關(guān)系數(shù)較低，分別為-0.114、-0.114和-0.056；2019年土壤有效磷與其他三項(xiàng)指標(biāo)呈正相關(guān)，相關(guān)系數(shù)較低，分別為0.267、0.222、0.139。綜上所述，雖然養(yǎng)分指標(biāo)具體數(shù)值存在差異，但表5的結(jié)果說(shuō)明有機(jī)質(zhì)、全氮和速效鉀3個(gè)指標(biāo)間的分布趨勢(shì)可能具有一定相似性，而有效磷的分布趨勢(shì)與其他指標(biāo)可能存在一定差異。

表5 宿州市土壤養(yǎng)分指標(biāo)相關(guān)性分析Table 5 Correlation analysis of soil nutrient indexes in Suzhou City

3 討論

宿州市耕地土壤4 項(xiàng)養(yǎng)分指標(biāo)含量的整體變化情況與全國(guó)耕地主要養(yǎng)分變化情況[15]一致，即耕地土壤有機(jī)質(zhì)含量總體上升、全氮總體穩(wěn)定、有效磷顯著上升、速效鉀穩(wěn)中有升。研究表明，坡度、水文、成土母質(zhì)、土壤類(lèi)型等自然條件和耕作習(xí)慣、種植方式、機(jī)械化生產(chǎn)水平、生產(chǎn)力等人為因素的差異都會(huì)導(dǎo)致區(qū)域間土壤養(yǎng)分水平的差異。9 年來(lái)埇橋區(qū)耕地土壤有機(jī)質(zhì)、有效磷、速效鉀含量顯著提升的原因是成土母質(zhì)、地形地貌、耕層質(zhì)地等自然因素引起的空間變異逐漸減小后，精耕細(xì)作、施肥與管理、種植模式等人類(lèi)活動(dòng)對(duì)土壤養(yǎng)分的影響逐漸增強(qiáng)[16-21]。大量研究表明施用化肥和有機(jī)肥可以提高養(yǎng)分含量[22]，根據(jù)2019 年宿州市統(tǒng)計(jì)年鑒，埇橋區(qū)第一產(chǎn)業(yè)的生產(chǎn)總值是68.43億元，肥料投入80 021萬(wàn)元，平均每公頃耕地投入化肥350.8 kg。宿州市耕地土壤有機(jī)質(zhì)和全氮含量空間分布高度相似，且兩個(gè)養(yǎng)分指標(biāo)的正相關(guān)性最強(qiáng)，這與有機(jī)質(zhì)中氮素含量相對(duì)固定的研究結(jié)論較為相符。蕭縣耕地土壤養(yǎng)分水平略低的區(qū)域是酒店、黃口、馬井等鄉(xiāng)鎮(zhèn)，其主要受土壤類(lèi)型及障礙因素的影響，這些鄉(xiāng)鎮(zhèn)多分布在黃河古道和大沙河河床兩側(cè)，受大沙河影響，耕層質(zhì)地較輕，土壤類(lèi)型以潮土的飛沙土屬為主[23]。由于飛沙土顆粒較粗，保水保肥能力弱，在后續(xù)改良利用中應(yīng)做好防治鹽堿化和改土培肥工作，完善水利配套設(shè)施，并輔以增施有機(jī)肥、種植綠肥、推進(jìn)秸稈還田等措施[24-26]。泗縣耕地土壤速效鉀明顯提升的原因主要是秸稈還田和施用鉀肥。泗縣山頭鎮(zhèn)、黑塔鎮(zhèn)、屏山鎮(zhèn)土壤由黃泛黏質(zhì)沉積物發(fā)育而成，土壤類(lèi)型主要是潮土的淤土，淤土質(zhì)地黏重，對(duì)鉀的固持能力較強(qiáng)，因而速效鉀含量較高[27]。碭山縣2010 年和2019 年的耕地土壤養(yǎng)分含量基本持平，且除有效磷含量較高外，其余養(yǎng)分指標(biāo)的含量水平均排在宿州市末位，碭山縣作為中國(guó)酥梨之鄉(xiāng)，9年來(lái)更側(cè)重于果園的施肥與管理。砂姜黑土主要分布在埇橋區(qū)、靈璧縣和泗縣南部，占宿州市耕地總面積的41.23%，9 年來(lái)砂姜黑土的耕地土壤養(yǎng)分水平明顯提升，主要是由于長(zhǎng)期耕作后土壤完成旱耕熟化過(guò)程[28-29]，在一系列灌排設(shè)施、施肥耕作后土壤形態(tài)特征發(fā)生潛育程度減輕、土體通氣條件改善等變化，提高了土壤養(yǎng)分有效性。

2010—2019 年宿州市耕地土壤有機(jī)質(zhì)、有效磷、速效鉀含量明顯提升而全氮基本持平的主要原因：一是施肥水平的提高，宿州市大力發(fā)展測(cè)土配方施肥[30]，截至2020 年測(cè)土配方施肥技術(shù)應(yīng)用面積占農(nóng)作物總播種面積的91%。宿州市根據(jù)土壤類(lèi)型、土壤養(yǎng)分條件、作物需肥規(guī)律和需肥量等，在施用有機(jī)肥的基礎(chǔ)上選擇氮、磷、鉀及中微量元素等肥料，因此2019 年宿州市耕地土壤養(yǎng)分水平較2010 年略有提升；二是科學(xué)合理的種植模式，如推廣秸稈還田技術(shù)和水肥一體化措施[31]，設(shè)置小麥-玉米（大豆）輪作區(qū)，將夏收小麥秸稈粉碎后覆蓋在玉米、大豆行間，將秋收玉米秸稈粉碎后翻壓還田種植小麥，農(nóng)民勤于耕作也會(huì)使耕地速效鉀含量顯著提升；三是就地增施有機(jī)肥[32]，宿州市養(yǎng)殖業(yè)的振興為農(nóng)作物施肥提供了優(yōu)質(zhì)的有機(jī)肥源保障，為化肥減量增效、土壤培肥提供了物質(zhì)基礎(chǔ)。根據(jù)本次土壤養(yǎng)分狀況調(diào)查結(jié)果，宿州市可以進(jìn)一步優(yōu)化養(yǎng)分管理模式，在穩(wěn)定施用氮肥、磷肥、鉀肥的基礎(chǔ)上，注重中微量元素的補(bǔ)給，通過(guò)開(kāi)展土壤測(cè)試、肥效試驗(yàn)等高效利用養(yǎng)分，綜合考慮土壤養(yǎng)分的豐缺程度、合理的目標(biāo)產(chǎn)量、自然環(huán)境等因素調(diào)整養(yǎng)分投入比例，逐步建立既保證作物增產(chǎn)又保護(hù)環(huán)境的施肥指標(biāo)體系。

4 結(jié)論

（1）2010—2019 年宿州市耕地土壤養(yǎng)分呈有機(jī)質(zhì)總體上升、全氮總體穩(wěn)定、有效磷顯著上升、速效鉀穩(wěn)中有升的變化趨勢(shì)，且4 項(xiàng)養(yǎng)分指標(biāo)含量均處于中等水平，中等級(jí)別面積分別占總耕地面積的76.00%、60.80%、60.33%和47.30%。

（2）從空間分布看，土壤有機(jī)質(zhì)含量呈現(xiàn)碭山縣整體較缺、蕭縣酒店鄉(xiāng)較缺、埇橋區(qū)中部鄉(xiāng)鎮(zhèn)較豐富、靈璧縣含量整體中等、泗縣個(gè)別鄉(xiāng)鎮(zhèn)較缺的分布規(guī)律；土壤全氮含量呈現(xiàn)碭山縣整體較缺、埇橋區(qū)北部鄉(xiāng)鎮(zhèn)缺乏、靈璧縣中部鄉(xiāng)鎮(zhèn)較豐、蕭縣西部鄉(xiāng)鎮(zhèn)較缺、泗縣南部鄉(xiāng)鎮(zhèn)較缺的分布規(guī)律；有效磷含量呈現(xiàn)碭山縣南部鄉(xiāng)鎮(zhèn)較豐、埇橋區(qū)北部較缺但南部較豐、蕭縣整體含量中等、泗縣南北部鄉(xiāng)鎮(zhèn)較缺的分布規(guī)律；土壤速效鉀含量呈現(xiàn)碭山縣中等、蕭縣西南鄉(xiāng)鎮(zhèn)較豐、埇橋區(qū)中部鄉(xiāng)鎮(zhèn)豐富但北部鄉(xiāng)鎮(zhèn)較缺、靈璧縣西部鄉(xiāng)鎮(zhèn)豐富、泗縣北部鄉(xiāng)鎮(zhèn)豐富的分布規(guī)律。

（3）9 年來(lái)碭山縣耕地土壤養(yǎng)分含量基本穩(wěn)定，埇橋區(qū)耕地土壤有機(jī)質(zhì)、有效磷、速效鉀含量顯著提升，蕭縣耕地土壤有機(jī)質(zhì)、有效磷、速效鉀含量略微降低，泗縣速效鉀含量顯著提升，靈璧縣速效鉀含量略微降低。砂姜黑土有機(jī)質(zhì)、有效磷、速效鉀含量略微提升，石灰?guī)r土和水稻土有機(jī)質(zhì)含量略微下降，各土類(lèi)全氮含量基本穩(wěn)定，黃褐土和潮土4 項(xiàng)耕地養(yǎng)分含量變化差異不顯著。