申桐，王恒飛，杜文波，周懷平，王瑞，張建杰，靳東升，徐明崗

褐土區(qū)典型縣域耕地土壤肥力時空演變特征及主控因素

申桐1,2，王恒飛1,2，杜文波3，周懷平1，王瑞3，張建杰1，靳東升1,2，徐明崗1,2

1山西農(nóng)業(yè)大學資源環(huán)境學院，山西太谷 030800；2山西農(nóng)業(yè)大學生態(tài)環(huán)境產(chǎn)業(yè)技術(shù)研究院/土壤環(huán)境與養(yǎng)分資源山西省重點實驗室，太原 030031；3山西省耕地質(zhì)量監(jiān)測保護中心，太原 030001

【目的】探明褐土區(qū)典型縣域耕地土壤肥力時空演變特征及主控因素，為合理施肥及土壤培肥提供科學依據(jù)?！痉椒ā炕谏轿魇坳柨h1983、2007和2017年的土壤肥力（包括土壤有機質(zhì)、全氮、有效磷、速效鉀含量和pH值）數(shù)據(jù)，采用內(nèi)梅羅指數(shù)法計算土壤肥力綜合指數(shù)；通過隨機森林（RF）方法定量土壤肥力演變的主控因素；利用GIS和地統(tǒng)計學相結(jié)合的方法，確定土壤肥力綜合指數(shù)及其主控因素的時空演變特征和分布格局?！窘Y(jié)果】（1）1983—2017年間壽陽縣耕地土壤肥力水平總體呈上升趨勢，肥力綜合指數(shù)由1.16增至1.42，增加了0.26。土壤肥力指標變化具有明顯的階段性，1983—2007年，有機質(zhì)、全氮和速效鉀增幅明顯，年均增加量分別為0.09 g·kg-1、0.0021 g·kg-1和1.61 mg·kg-1，而有效磷增幅不顯著；2007—2017年，有機質(zhì)、全氮和有效磷增幅明顯，年均增加量分別為0.20 g·kg-1、0.01 g·kg-1、0.24mg·kg-1，而速效鉀增幅不顯著。（2）RF分析結(jié)果顯示，1983年壽陽縣耕地土壤肥力的主控因素為全氮和有機質(zhì)，貢獻率分別為75.3%和17.8%；2007年土壤肥力的主控因素則為速效鉀、有效磷和全氮，貢獻率分別為31.8%、27.1%和26.8%；2017年變?yōu)槿?、速效鉀和有效磷，貢獻率分別為31.8%、27.1%和26.8%。（3）壽陽縣耕地土壤肥力主控因素的變化存在顯著的空間差異特征。1983—2007年，全域土壤有機質(zhì)整體呈現(xiàn)增加趨勢；全氮除在西北部和中部降低外，在其他區(qū)域均呈增加趨勢；有效磷則呈現(xiàn)北部降低、南部增加的變化特征；速效鉀在全域呈增加趨勢。2007—2017年，有機質(zhì)在東南部增加較快，其他區(qū)域增加較慢；全氮則呈現(xiàn)在縣域中部減少、其他區(qū)域增加的變化趨勢；有效磷在縣域東部下降，其他區(qū)域上升；速效鉀則呈現(xiàn)在東部和西部下降、中部上升的空間變化特征?！窘Y(jié)論】34年以來，壽陽縣耕地土壤肥力總體呈上升趨勢，有機質(zhì)、全氮、有效磷和速效鉀是影響土壤肥力演變的主控因素，有機質(zhì)、全氮和有效磷在南部增長較快，速效鉀則在東部增長較快。未來建議壽陽縣采取全區(qū)尤其是中部增氮，全區(qū)穩(wěn)磷控鉀的養(yǎng)分調(diào)控策略。

土壤肥力；時空演變；褐土區(qū)；主控因素；壽陽縣

0 引言

【研究意義】褐土是我國主要土壤類型之一，廣泛分布于我國北方地區(qū)，包括內(nèi)蒙古、北京、山西、河南、陜西和甘肅等地，其中，山西省是我國褐土分布最集中和典型的區(qū)域[1-2]。褐土區(qū)作為我國玉米和小麥等糧食作物的主產(chǎn)區(qū)之一，揭示其肥力時空演變特征及主控因素對于今后提升耕地質(zhì)量至關(guān)重要。壽陽縣地處褐土典型區(qū)域，農(nóng)業(yè)條件相對較好，發(fā)展?jié)摿Υ?，是山西省重要的農(nóng)業(yè)主產(chǎn)區(qū)。科學評價壽陽縣耕地土壤肥力及其主控因素的時空演變特征，可為褐土農(nóng)田肥力提升、實現(xiàn)農(nóng)業(yè)可持續(xù)發(fā)展提供理論依據(jù)?！厩叭搜芯窟M展】自然因素與人為因素共同影響土壤肥力的演變，已有學者針對山西省褐土區(qū)土壤肥力指標的時空變化特征做了大量的研究工作。但大多都局限于分析土壤肥力指標的靜態(tài)變化或單一指標的動態(tài)變化。例如，張建杰等[3-4]分別研究了太原市和臨汾盆地全氮和有機質(zhì)的空間變化特征，并闡釋了其影響因素；申若禹等[5]對山西省不同類型土壤有機質(zhì)含量研究發(fā)現(xiàn)，褐土中有機質(zhì)含量最高，栗鈣土有機質(zhì)含量最低；解文艷等[6]研究發(fā)現(xiàn)，瀟河流域有機質(zhì)含量1982—2017年呈增加趨勢，施肥是影響其時空分異的主要因素。上述研究對某一肥力指標的時空變化特征進行了定量表征，并且分析了主要影響因素，然而未能準確刻畫區(qū)域尺度整體的肥力時空演變規(guī)律。目前，針對褐土多個肥力指標多時間段的時空演變特征研究鮮見報道，已有研究也僅僅局限于服務(wù)耕地質(zhì)量的分等定級[7]，不能夠從科學層面反映耕地質(zhì)量變化的原因，更不能夠為耕地土壤肥力的提高提供依據(jù)。此外，對于土壤肥力演變主控因素的探究，國內(nèi)學者在其他省市已經(jīng)開展了相應(yīng)工作。例如，陳延華等[8]基于29年間全國28個褐土區(qū)長期試驗點數(shù)據(jù)，發(fā)現(xiàn)褐土區(qū)土壤生產(chǎn)力受有機肥用量和土壤全氮、有機質(zhì)含量的影響最大；李建軍等[9]對長江中下游糧食主產(chǎn)區(qū)研究發(fā)現(xiàn)，初期土壤肥力主控因素為有機質(zhì)、全氮和堿解氮，后期主控因素變?yōu)槿?、堿解氮和速效鉀；周波等[10]研究指出，深圳市土壤肥力的限制因子由初期的速效養(yǎng)分變?yōu)楹笃诘牡兀N植模式和施肥的改變是導(dǎo)致土壤肥力限制因子發(fā)生轉(zhuǎn)變的主要原因。綜上發(fā)現(xiàn)，耕地土壤肥力的主控因素為全氮、有機質(zhì)、速效養(yǎng)分等土壤肥力指標。目前在山西省褐土區(qū)耕地土壤肥力研究中，針對肥力演變主控因素的探究還鮮見報道?！颈狙芯壳腥朦c】近年來，隨著土地利用方式、農(nóng)田管理措施等的改變，褐土區(qū)肥力發(fā)生了明顯變化，然而其土壤肥力指數(shù)及其主控因素的時空演變特征尚不明確。山西省壽陽縣是典型的褐土區(qū)，明確該縣耕地土壤肥力時空演變特征及其主控因素，對于分析全國褐土肥力演變也具有參考價值。【擬解決的關(guān)鍵問題】基于壽陽縣耕地土壤肥力長期監(jiān)測和調(diào)查數(shù)據(jù)，分析該縣耕地土壤肥力綜合指數(shù)及其主控因素的時空變異特征，為褐土區(qū)耕地土壤培肥和養(yǎng)分管理提供支撐。

1 材料與方法

1.1 區(qū)域概況

壽陽縣位于山西省東部，居黃土高原東部瀟河流域中上游地區(qū)（37°34′—38°5′N， 112°6′—113°28′E），氣候類型為溫帶大陸性氣候，年平均氣溫7.4 ℃，年降水量518.3 mm，無霜期140 d左右。土壤類型主要為褐土，占全縣土壤總面積的96.8%，成土母質(zhì)主要有黃土及黃土狀母質(zhì)、殘積、坡積母質(zhì)和沖積母質(zhì)等。農(nóng)業(yè)種植方式為一年一季，主要農(nóng)作物以玉米、馬鈴薯、小雜糧為主。

1.2 樣品采集與分析

1983年耕地土壤屬性數(shù)據(jù)收集于第二次土壤普查資料，樣本數(shù)量為726個，空間分布均勻，有較強代表性。2007年前后和2017年前后的土壤肥力數(shù)據(jù)均為采樣實測數(shù)據(jù)，數(shù)量分別為5 667和1 401個。采樣時，以第二次土壤普查形成的土壤類型分布圖為參考，依據(jù)相對均勻的空間分布原則，在播種前按照五點法采集耕層（0—20 cm）土壤樣品，并利用GPS記錄每個采樣點的經(jīng)緯度，采樣點分布如圖1所示。將采集的土樣帶回實驗室風干，剔除動植物殘體及石塊等雜質(zhì)并磨細，分別過2和0.149 mm篩，用于土壤理化性質(zhì)的測定。其中，土壤有機質(zhì)含量采用重鉻酸鉀氧化外加熱法測定；全氮含量采用凱氏定氮法測定；有效磷含量采用碳酸氫鈉浸提-鉬銻抗比色法測定；速效鉀含量采用NH4Ac浸提，火焰光度法測定；pH采用水土比為2.5﹕1的玻璃電極法進行測定。

圖1 1983、2007和2017年壽陽縣耕地土壤采樣點分布圖

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

利用Excel 2016進行數(shù)據(jù)整理，采用閾值法[11]對異常值進行剔除。采用SPSS 23.0完成描述性統(tǒng)計分析、顯著性分析（LSD差異分析）及正態(tài)分布檢驗。利用R語言中的“Random-Forest”軟件包進行隨機森林的重要性分析。在GS+9.0軟件中完成最優(yōu)半方差參數(shù)及模型擬合。在ArcGIS10.8的地統(tǒng)計分析模塊（Geostatistical Analyst）中，運用克里金插值繪制土壤肥力主要指標及綜合指數(shù)的空間分布圖。

1.4 評價方法

采用修正后的內(nèi)梅羅指數(shù)法進行土壤肥力評價[12]。首先，按照標準化處理方法對采樣點的有機質(zhì)、全氮、有效磷、速效鉀和pH進行標準化，以消除各評價因子之間的量綱差別[13-14]，參評因子分級標準見表1。之后，采取較為客觀的修正后的內(nèi)梅羅（Nemerow N. C.）公式評價土壤肥力[14]。

最后，根據(jù)土壤肥力指數(shù)，將其等級分為4級：一級（≥2）、二級（：＜2.0—1.5）、三級（：＜1.5—1.0）、四級（＜1.0）[14]。土壤肥力主要指標（有機質(zhì)、全氮、有效磷、速效鉀）等級劃分參照第二次土壤普查標準[15]。

表1 土壤質(zhì)量評價因子分級標準

2 結(jié)果

2.1 壽陽縣耕地土壤肥力指標時間演變特征

1983—2017年間，壽陽縣耕地有機質(zhì)、全氮、有效磷、速效鉀與pH整體均呈現(xiàn)不同程度的上升趨勢（表2）。2007和2017年，壽陽縣耕地土壤有機質(zhì)均值為12.8和14.8 g·kg-1，分別較1983年（10.7 g·kg-1）增加了19.6%和38.3%。1983—2007年有機質(zhì)年均增加量為0.1 g·kg-1；2007—2017年均增加量為0.2 g·kg-1，增速是前者的2倍。

2007和2017年壽陽縣耕地全氮均值為0.71和0.83 g·kg-1，分別較1983年（0.66 g·kg-1）增加了7.6%和25.8%。1983—2007年全氮年均增加量為0.0021 g·kg-1；2007—2017年均增加量為0.01 g·kg-1，其增速是前者的4.8倍。

壽陽縣耕地土壤有效磷含量，1983年（10.4 mg·kg-1）與2007年（10.9 mg·kg-1）無顯著差異；而2017年（13.3 mg·kg-1）則較2007和1983年分別顯著提高了22.0%和27.9%。土壤有效磷年均變化量也呈現(xiàn)為階段性變化特征，1983—2007年無顯著變化，2007—2017年均增加量為0.2 mg·kg-1，增幅明顯。

表2 1983、2007和2017年壽陽縣耕地土壤肥力指標描述性統(tǒng)計結(jié)果

同一欄內(nèi)均值后不同字母表示不同年份同一土壤肥力指標存在顯著差異（＜0.05）

Different letters after the mean values in the same column indicate significant differences in the same soil fertility index in different years (＜0.05)

2007和2017年土壤速效鉀均值為135和136 mg·kg-1，分別較1983年（96 mg·kg-1）增加了40.6%和41.7%；1983—2007年速效鉀年均增加量為1.6 mg·kg-1，增幅顯著，而2007—2017年速效鉀無顯著變化。

2007和2017年土壤pH均為8.28，較1983年（8.14）略有增加，但均在作物適宜的pH范圍內(nèi)。

2.2 不同時期壽陽縣耕地土壤肥力演變主控因素

34年間壽陽縣耕地土壤肥力的主控因素發(fā)生了較大變化（圖2）。1983年全氮和有機質(zhì)是土壤肥力的主要貢獻指標，貢獻率分別為75.3%和17.8%，其他指標的重要性相對較弱；2007年，速效鉀、有效磷和全氮是土壤肥力的主要貢獻指標，貢獻率分別為34.8%、27.9%和27.8%；2017年，對土壤肥力的主要貢獻轉(zhuǎn)變?yōu)槿⑺傩р浐陀行Я?，貢獻率分別為31.8%、27.1%和26.8%。速效養(yǎng)分在后期發(fā)揮的作用越來越重要。

圖2 1983、2007和2017年壽陽縣各指標對土壤肥力的貢獻率

2.3 壽陽縣耕地土壤肥力主控因素及綜合指數(shù)的時空演變特征

2.3.1 耕地土壤肥力主控因素的時空分布 1983、2007和2017年壽陽縣耕地土壤肥力主控因素（有機質(zhì)、全氮、有效磷和速效鉀）存在明顯的時空差異性（圖3、圖4、表3）。

1983年壽陽縣耕地土壤有機質(zhì)含量呈現(xiàn)北部高，南部低的分布格局，有機質(zhì)含量整體偏低，以四、五等級為主。2007和2017年土壤有機質(zhì)空間分布相對均勻，且與1983年相比，四等級的比例顯著增加，分別增至接近94%（2007年）和98%（2017年）（表3）。34年來該縣南部土壤有機質(zhì)含量逐漸上升，而北部區(qū)域土壤有機質(zhì)的變化不顯著。

1983年壽陽縣耕地土壤全氮含量的分布相對均勻，全氮含量整體偏低，且以五等級為主。2007年該縣土壤全氮含量以四、五等級為主，其中，在平頭鎮(zhèn)和馬首鄉(xiāng)土壤全氮降低，等級由五級降為六級；北部和南部土壤全氮增加，等級由五級增長為四級。至2017年，該縣土壤全氮含量整體呈增加趨勢，其中在南部區(qū)域增加最為顯著。34年來，該縣西北部和中部的土壤全氮含量先減少后增加，而在其他區(qū)域土壤全氮均有不同程度的增加，其中在北部和南部的邊緣地區(qū)，土壤全氮的增加最為明顯。整體來看，該縣土壤全氮含量處于四等級的中等偏低水平。

1983年壽陽縣土壤有效磷含量從北到南呈現(xiàn)遞減的趨勢，且以三和四等級為主。2007年，該縣土壤有效磷含量呈東南部上升趨勢，依舊以三和四等級為主。至2017年，該縣土壤有效磷含量整體增加，其中在西南部土壤有效磷由原來的四級上升至三級。34年來，該縣南部區(qū)域土壤有效磷含量增加明顯，而北部區(qū)域變化不顯著。整體來看，該縣土壤有效磷含量處于三、四等級的中等偏上水平。

1983年壽陽縣土壤速效鉀含量呈現(xiàn)西南部和東部較高，北部和中部較低的分布格局，且以四等級為主。至2007年，該縣土壤速效鉀含量整體上升，以三等級為主，占比達90%（表3），西南部部分區(qū)域達二級水平。至2017年，該縣耕地土壤速效鉀含量變化不顯著，依舊以三等級為主。整體看來，34年來該縣耕地土壤速效鉀含量顯著增加，且在北部和中部地區(qū)增加最為顯著。

圖3 1983、2007和2017年壽陽縣土壤有機質(zhì)和全氮的空間分布

圖4 1983、2007和2017年壽陽縣土壤有效磷和速效鉀的空間分布

2.3.2 耕地土壤肥力主控因素的年均變化率 土壤肥力主控因素的年均變化率（圖5和圖6）呈現(xiàn)階段性特征。1983—2007年土壤有機質(zhì)年均增加率在0—5%的占比最高，為95.9%；2007—2017年該縣東南部景尚鄉(xiāng)、羊頭崖鄉(xiāng)和松塔鎮(zhèn)年均增加速率最大，達5%—10%。整體看來，2007—2017年有機質(zhì)含量年均增加率較1983—2007年更大。

1983—2007年，壽陽縣西北部和中部全氮含量呈現(xiàn)下降趨勢，占比達33.8%，其他區(qū)域全氮含量增加0—5%，占比達52.0%。2007—2017年，該縣土壤全氮年均變化率整體呈現(xiàn)中北部下降，其他區(qū)域增加的趨勢，尤其在平頭鎮(zhèn)和馬首鄉(xiāng)增速最大，年均增加率達5%—10%。總體看來，2007—2017年全氮含量年均增加率較1983—2007年更大。

1983—2007年，壽陽縣土壤有效磷增加0—5%的占比最高，達70.0%。2007—2017年間增加0—5%占比為62.5%，其中，在中部區(qū)域土壤有效磷增加顯著，上升幅度達10%—20%。整體看來，2007—2017年有效磷含量年均增加率較1983—2007年更大。

1983—2007年，壽陽縣土壤速效鉀增加在0—5%占比最高，達94.4%。2007—2017年，速效鉀含量整體呈現(xiàn)東西部下降，中部上升的趨勢，其中速效鉀年均減少占比為51.2%，年均增加占比為48.8%。整體來看，1983—2007年速效鉀含量年均增加率較2007— 2017年更大。

2.3.3 耕地土壤肥力綜合指數(shù)的時空分布 圖7顯示，1983年壽陽縣土壤肥力綜合指數(shù)平均值為1.16，且以三級為主，占比為64%，其次為四級，占比為36%。2007年和2017年土壤肥力綜合指數(shù)平均值分別為1.31和1.42，2007年三級占比達92%，2017年三級占比為82%，且在中部和北部區(qū)域出現(xiàn)二級水平，占比為13%。整體來看，1983—2017年，該縣土壤肥力綜合指數(shù)增加0.26，其中第三等級占比不斷增大，而第四等級占比逐漸減少，且部分區(qū)域土壤肥力指數(shù)出現(xiàn)二級水平。

圖5 1983—2007和2007—2017年壽陽縣土壤有機質(zhì)和全氮的年均變化率

表3 1983、2007和2017年壽陽縣土壤肥力等級面積比例統(tǒng)計結(jié)果

圖6 1983—2007和2007—2017年壽陽縣土壤有效磷和速效鉀的年均變化率

圖7 1983、2007和2017年壽陽縣土壤肥力綜合指數(shù)的空間分布

3 討論

3.1 不同時期土壤肥力主控因素的變化

土壤肥力主控因素的改變表示著各指標對綜合肥力指數(shù)的貢獻率發(fā)生變化。本研究表明，在1983—2017年間，壽陽縣耕地土壤肥力演變的主控因素由最初的全氮和有機質(zhì)（1983年），變?yōu)槿?、有效磷和速效鉀?007和2017年）。34年間，全氮始終是土壤肥力的主要貢獻因子，但其相對重要性逐漸降低，而有效磷、速效鉀養(yǎng)分對壽陽縣土壤肥力的貢獻逐漸增加。這可能是因為壽陽縣土壤全氮初始含量較低，僅為0.66 g·kg-1，且1983年還沒有大規(guī)模施用化肥，壽陽縣氮肥年平均施用量僅為47 kg·hm-2，之后，2007和2017年，肥料施用大幅度增加，氮肥施用量分別達133和246 kg·hm-2（表4）。隨著氮肥施用量的增加，全氮含量提升顯著，全氮一直為土壤肥力主控因素，而土壤有機質(zhì)在2007和2017年分布相對均勻，空間差異小，不再是土壤肥力的主控因素。此外，隨著化肥的投入，磷肥和鉀肥施用量也不斷增加，2007和2017年磷肥施用量分別為71和91 kg·hm-2，鉀肥施用量分別為32和50 kg·hm-2（表4），導(dǎo)致有效磷、速效鉀養(yǎng)分也成為土壤肥力的主控因素。與我們的結(jié)果相類似，李建軍等[9]指出，主要肥力貢獻因子的變化與化肥投入量的持續(xù)增加有關(guān)，在監(jiān)測初期長江中下游地區(qū)土壤肥力的主控因素為有機質(zhì)、全氮和堿解氮，在監(jiān)測后期主控因素變?yōu)槿?、堿解氮和速效鉀。此外，王遠鵬等[16]指出，進賢縣土壤肥力主控因素由1982年的堿解氮和有效磷變?yōu)?017年的速效鉀和有效磷。整體來看，土壤全氮和速效養(yǎng)分仍然是壽陽縣土壤肥力的主控因素。因此，未來壽陽縣農(nóng)田土壤管理中應(yīng)注意氮、磷、鉀肥料的合理施用。

表4 1983、2007和2017年壽陽縣農(nóng)田施肥情況和農(nóng)作物秸稈總量

3.2 不同時期土壤肥力主控因素及綜合指數(shù)時空變化

3.2.1 土壤有機質(zhì)和全氮時空變化 土壤有機質(zhì)影響?zhàn)B分儲存、持水能力以及團聚體穩(wěn)定性，是土壤質(zhì)量和肥力的核心[17]。34年來，壽陽縣土壤有機質(zhì)整體呈現(xiàn)上升趨勢，這與解文艷等[6]的研究結(jié)果相似。通常，土壤有機質(zhì)主要來源于作物根茬、秸稈還田、有機肥施用以及綠肥的翻壓等，其受自然因素和人為因素的共同影響，且人為因素（農(nóng)田管理措施如施肥等）對土壤有機質(zhì)含量的影響更大[17]。前人研究表明，有機肥含有大量氮磷鉀等營養(yǎng)物質(zhì)和有益菌種，與化肥配施后，可以促進土壤微生物代謝，加快養(yǎng)分循環(huán)，改善土壤結(jié)構(gòu)，提高土壤肥力[18-20]。此外，與秸稈不還田相比，秸稈還田可以明顯提高土壤有機質(zhì)含量[21]。本研究表明，與1983—2007年相比，2007—2017年間壽陽縣土壤有機質(zhì)含量增加更加明顯。這可能與2007年以來，壽陽縣政府大力推廣秸稈還田和有機肥施用有關(guān)。自2005年以來，壽陽縣政府大力實施沃土計劃推行玉米秸稈還田，使得大量秸稈進入土壤，其中2008年壽陽縣實現(xiàn)80%以上秸稈還田，且2016年有機肥增施量達1 000 kg。表4表明，1983—2017年壽陽縣秸稈總量逐年增長，這也意味著自2005年以來，壽陽縣秸稈還田量在增加。因此，壽陽縣耕地土壤有機質(zhì)顯著增加顯然與政府大力推廣秸稈還田和有機肥增施策略有關(guān)。此外，從空間變化來看，壽陽縣南部區(qū)域有機質(zhì)增長更快，這是因為南部地區(qū)是農(nóng)業(yè)重點發(fā)展區(qū)域，耕地多，主要種植作物為玉米，秸稈量大，還田較多；而中部呈現(xiàn)降低趨勢，這是因為中部為鄉(xiāng)鎮(zhèn)工業(yè)發(fā)展區(qū)，其工業(yè)產(chǎn)值占總產(chǎn)值的43%—56%，耕地面積較少，秸稈還田量少。綜上，壽陽縣增施有機肥和秸稈還田的措施是導(dǎo)致該縣34年來土壤有機質(zhì)含量增加的重要措施。

土壤全氮作為土壤肥力的重要內(nèi)容，是農(nóng)田土壤氮肥施用的重要依據(jù)。通常，土壤全氮和土壤有機質(zhì)密不可分。34年來，壽陽縣土壤全氮含量呈現(xiàn)整體上升趨勢。代子俊等[22]指出，自然因素和人為因素共同導(dǎo)致了氮的時空分異，人類活動如施肥對其干預(yù)作用更強。前人研究表明，土壤氮素主要以有機態(tài)氮和硝態(tài)氮等的形式存在于土壤有機質(zhì)中[17]，增施有機肥能明顯增加土壤有機氮庫的儲量，提高土壤供氮能力[23]。馬力等[24]研究也表明，長期施用有機肥以及秸稈還田可以明顯提高耕層土壤氮含量。1983—2007年間，壽陽縣土壤全氮含量變化速度較慢，這可能是因為這一時期氮肥施用量相對較少，僅由1983年的47 kg·hm-2增加至2007年的133 kg·hm-2。在2007—2017年，土壤全氮含量顯著增加，這與氮肥和有機肥施用量增加以及秸稈還田有關(guān)。表4結(jié)果顯示，從2007到2017年，壽陽縣氮肥施用量由133 kg·hm-2增加至246 kg·hm-2，且農(nóng)作物秸稈量增加，也意味著秸稈還田量的增加。從空間變化看，全氮變化與有機質(zhì)相似，均呈現(xiàn)中部降低的趨勢，這也與中部重點發(fā)展工業(yè)有關(guān)。因此，壽陽縣增施氮肥和有機肥以及秸稈還田等措施是導(dǎo)致該縣34年來土壤全氮含量增加的主要原因。

3.2.2 土壤速效磷和速效鉀時空變化 土壤有效磷是反映土壤供磷性狀的重要指標，其含量高低直接影響到作物對土壤磷素的吸收利用[25]。土壤速效鉀含量對于判斷土壤中鉀素供應(yīng)狀況具有重要意義[26]。34年來，壽陽縣耕地土壤有效磷和速效鉀含量整體呈現(xiàn)上升趨勢，表現(xiàn)為四級面積占比減少，三級面積占比增加。王齊齊等[27]的研究表明，土壤速效養(yǎng)分含量的增加與化肥施用量增加有關(guān)。前人研究表明，有機無機配施有利于土壤中活性磷的轉(zhuǎn)化[28]，進而提高土壤有效磷含量[23,29]，而秸稈還田后分解產(chǎn)生的大量有機酸也可以為土壤提供陰離子，減少無機磷的固定，提高磷的有效性[30-31]。1983—2007年間，壽陽縣耕地土壤有效磷沒有顯著變化而速效鉀含量顯著增加，這與施肥量密切相關(guān)。2007年施磷量（71 kg·hm-2）較1983年（60 kg·hm-2）沒有顯著變化，而施鉀量由1983的4 kg·hm-2提高至2007年32 kg·hm-2。因此，這一時期鉀肥施用量的增大是導(dǎo)致速效鉀含量顯著增加的重要因素。2007—2017年間，土壤有效磷顯著增加而速效鉀沒有顯著變化。這是因為2007—2017年間，壽陽縣開始進行測土配方施肥，推廣有機無機配施和秸稈還田，該時期施磷量增加，2017年是2007年的1.9倍，而施鉀量變化較小。因此，磷肥和鉀肥施用量的增加、推廣測土配方施肥以及秸稈還田等措施可能是導(dǎo)致34年來壽陽縣耕地土壤有效磷鉀含量增加的重要原因。

3.2.3 土壤肥力綜合指數(shù)時空變化 土壤肥力綜合指數(shù)可以較為全面地反映土壤的肥力特征。34年間壽陽縣褐土肥力綜合指數(shù)逐漸增加，這與趙秀娟[31]的研究結(jié)果一致。此外，1983年壽陽縣邊緣地區(qū)土壤肥力較差，而到2017年，壽陽縣土壤肥力水平整體到達三級，且空間分布均勻，尤其邊緣地區(qū)土壤肥力提升明顯。這可能與壽陽縣自2007年以來實行的測土配方施肥、增施有機肥以及推行秸稈還田等政策有關(guān)。

綜上分析，34年來壽陽縣褐土區(qū)耕地土壤肥力的提高與化肥和有機肥施用以及農(nóng)作物秸稈還田等人為措施密切相關(guān)。由于大面積采樣工作量大，且時間跨度較長，本研究只有3次采樣數(shù)據(jù)結(jié)果，因此可能并未能十分精準反映出34年以來壽陽縣耕地土壤肥力時空演變的情況，但其大致演變趨勢應(yīng)對未來壽陽縣褐土耕地質(zhì)量提升和作物施肥管理具有一定的科學指導(dǎo)意義。此外，受制于數(shù)據(jù)的局限性，本文只關(guān)注了34年間壽陽縣施肥與秸稈還田等人為措施對土壤肥力綜合指數(shù)和基礎(chǔ)5項肥力因子（有機質(zhì)、全氮以及速效養(yǎng)分等）時空演變的影響，而其他因素，比如土壤物理性質(zhì)、生物性質(zhì)、氣候、地形地貌等影響土壤肥力的可能因素均未能包括。因此，今后研究中仍需加強長時間尺度下，自然和人為管理等因素對褐土肥力以及質(zhì)量的綜合影響研究。

4 結(jié)論

1983—2017年間壽陽縣褐土肥力水平整體提高，綜合指數(shù)增加了0.26，當前耕地土壤肥力整體處于三級中等水平。土壤肥力的主控因素前期以全氮和有機質(zhì)為主導(dǎo)，而后期有效磷和速效鉀養(yǎng)分的作用愈加重要。

全縣有機質(zhì)、全氮、有效磷和速效鉀均呈現(xiàn)不同程度的上升趨勢，其中有機質(zhì)、全氮和有效磷在南部增長較快，速效鉀在中部增長較快，2017年有機質(zhì)、全氮、有效磷和速效鉀含量均值分別為14.8 g·kg-1、0.83 g·kg-1、13.3 mg·kg-1和136 mg·kg-1。

全縣大部分區(qū)域全氮處于四等級及以下的中等偏低水平，尤其中部地區(qū)全氮含量偏低。有效磷和速效鉀含量較高，大部分處于二和三等級，中等偏上水平。因此，在壽陽縣全域尤其是中部區(qū)域需注意增施氮肥，而全區(qū)需穩(wěn)磷控鉀。

[1] 陳延華, 王樂, 張淑香, 郭寧, 馬常寶, 李春花, 徐明崗, 鄒國元. 我國褐土耕地質(zhì)量的演變及對生產(chǎn)力的影響. 中國農(nóng)業(yè)科學, 2019, 52(24): 4540-4554. doi: 10.3864/j.issn.0578-1752.2019.24.009.

CHEN Y H, WANG L, ZHANG S X, GUO N, MA C B, LI C H, XU M G, ZOU G Y. Quality change of cinnamon soil cultivated land and its effect on soil productivity. Scientia Agricultura Sinica, 2019, 52(24): 4540-4554. doi: 10.3864/j.issn.0578-1752.2019.24.009. (in Chinese)

[2] 于婧文. 褐土區(qū)鄉(xiāng)鎮(zhèn)級土壤養(yǎng)分的時空變異研究[D]. 太原: 山西大學, 2010.

YU J W. Spatiotemporal variation research of soil nutrient on cinnamon soil in town[D]. Taiyuan: Shanxi University, 2010. (in Chinese)

[3] 張建杰, 李富忠, 胡克林, 張強, 郭彩霞, 黃元仿. 太原市農(nóng)業(yè)土壤全氮和有機質(zhì)的空間分布特征及其影響因素. 生態(tài)學報, 2009, 29(6): 3163-3172.

ZHANG J J, LI F Z, HU K L, ZHANG Q, GUO C X, HUANG Y F. Spatial characteristics and impact factors of soil total nitrogen and soil organic matter in Taiyuan. Acta Ecologica Sinica, 2009, 29(6): 3163-3172. (in Chinese)

[4] 張建杰, 張強, 楊治平, 陳明昌, 李勇, 陳德利. 山西臨汾盆地土壤有機質(zhì)和全氮的空間變異特征及其影響因素. 土壤通報, 2010, 41(4): 839-844.

ZHANG J J, ZHANG Q, YANG Z P, CHEN M C, LI Y, CHEN D L. Spatial variability of topsoil organic matter and total nitrogen in Linfen Basin, Shanxi and its influencing factors. Chinese Journal of Soil Science, 2010, 41(4): 839-844. (in Chinese)

[5] 申若禹, 張吳平, 王國芳, 張茜, 喬磊. 基于地統(tǒng)計學的山西省不同類型土壤有機質(zhì)空間變異分析研究. 山東農(nóng)業(yè)科學, 2019, 51(4): 110-116.

SHEN R Y, ZHANG W P, WANG G F, ZHANG Q, QIAO L. Spatial variability analysis of organic matter in different types of soil based on geostatistics in Shanxi Province. Shandong Agricultural Sciences, 2019, 51(4): 110-116. (in Chinese)

[6] 解文艷, 周懷平, 楊振興, 馮悅晨, 白雪, 杜艷玲. 黃土高原東部瀟河流域農(nóng)田土壤有機質(zhì)時空變異及影響因素. 農(nóng)業(yè)資源與環(huán)境學報, 2019, 36(1): 96-104.

XIE W Y, ZHOU H P, YANG Z X, FENG Y C, BAI X, DU Y L. The spatial-temporal variation of soil organic matter and its influencing factors in Xiaohe River Basin in eastern Loess Plateau, China. Journal of Agricultural Resources and Environment, 2019, 36(1): 96-104. (in Chinese)

[7] 趙亮. 山西省2010—2013年度耕地質(zhì)量等別分布情況分析. 華北國土資源, 2016(1): 120-123.

ZHAO L. Analysis on the graded distribution of cultivated land quality in Shanxi Province from 2010 to 2013. Huabei Land and Resources, 2016(1): 120-123. (in Chinese)

[8] 陳延華, 王樂, 張淑香, 任意, 李春花, 徐明崗, 趙同科. 長期施肥下褐土生產(chǎn)力的演變及其影響因素. 植物營養(yǎng)與肥料學報, 2018, 24(6): 1445-1455.

CHEN Y H, WANG L, ZHANG S X, REN Y, LI C H, XU M G, ZHAO T K. Evolution of cinnamon soil productivity and the main influencing factors under long-term fertilization. Journal of Plant Nutrition and Fertilizers, 2018, 24(6): 1445-1455. (in Chinese)

[9] 李建軍, 辛景樹, 張會民, 段建軍, 任意, 孫楠, 徐明崗. 長江中下游糧食主產(chǎn)區(qū)25年來稻田土壤養(yǎng)分演變特征. 植物營養(yǎng)與肥料學報, 2015, 21(1): 92-103.

LI J J, XIN J S, ZHANG H M, DUAN J J, REN Y, SUN N, XU M G. Evolution characteristics of soil nutrients in the main rice production regions, the middle-lower reach of Yangtze River of China. Journal of Plant Nutrition and Fertilizer, 2015, 21(1): 92-103. (in Chinese)

[10] 周波, 陳旭飛, 張聰俐, 李靜娟, 張靜, 任宗玲, 張池, 戴軍. 深圳市光明新區(qū)土壤肥力時空演變的主成分分析. 華南農(nóng)業(yè)大學學報, 2012, 33(4): 448-452.

ZHOU B, CHEN X F, ZHANG C L, LI J J, ZHANG J, REN Z L, ZHANG C, DAI J. Principal component analysis of spatial and temporal evolution of soil fertility in Guangming new district, Shenzhen city. Journal of South China Agricultural University, 2012, 33(4): 448-452. (in Chinese)

[11] 康日峰, 任意, 吳會軍, 張淑香. 26年來東北黑土區(qū)土壤養(yǎng)分演變特征. 中國農(nóng)業(yè)科學, 2016, 49(11): 2113-2125. doi: 10.3864/j.issn. 0578-1752.2016.11.008.

KANG R F, REN Y, WU H J, ZHANG S X. Changes in the nutrients and fertility of black soil over 26 years in northeast China. Scientia Agricultura Sinica, 2016, 49(11): 2113-2125. doi: 10.3864/j.issn.0578- 1752.2016.11.008. ( in Chinese)

[12] 毛偉, 李文西, 陳明, 龔鑫鑫, 高暉, 張月平, 褚宏遠, 趙海濤, 柏彥超, 錢曉晴. 近30年揚州市耕地土壤肥力變異特征及其驅(qū)動因素分析. 植物營養(yǎng)與肥料學報, 2020, 26(11): 1998-2009.

MAO W, LI W X, CHEN M, GONG X X, GAO H, ZHANG Y P, CHU H Y, ZHAO H T, BAI Y C, QIAN X Q. Evolution of soil fertility and the main driving factors in Yangzhou city during the past 30 years. Journal of Plant Nutrition and Fertilizers, 2020, 26(11): 1998-2009. (in Chinese)

[13] 王緒奎, 徐茂, 汪吉東, 王建明, 張永春. 太湖地區(qū)典型水稻土大時間尺度下的肥力質(zhì)量演變. 中國生態(tài)農(nóng)業(yè)學報, 2009, 17(2): 220-224.

WANG X K, XU M, WANG J D, WANG J M, ZHANG Y C. Response of typical paddy soil fertility to long-term cultivation in Taihu Lake Region. Chinese Journal of Eco-Agriculture, 2009, 17(2): 220-224. (in Chinese)

[14] 秦明周, 趙杰. 城鄉(xiāng)結(jié)合部土壤質(zhì)量變化特點與可持續(xù)性利用對策: 以開封市為例. 地理學報, 2000, 55(5): 545-554.

QIN M Z, ZHAO J. Strategies for sustainable use and characteristics of soil quality changes in urban-rural marginal area: A case study of Kaifeng. Acta Geographica Sinica, 2000, 55(5): 545-554. (in Chinese)

[15] 全國土壤普查辦公室. 中國土壤普查技術(shù). 北京: 農(nóng)業(yè)出版社, 1992: 111-112.

National Soil Census Office. Soil General Survey Technology in China. Beijing: Agricultural Press, 1992: 111-112. (in Chinese)

[16] 王遠鵬, 黃晶, 孫鈺翔, 柳開樓, 周虎, 韓天富, 都江雪, 蔣先軍, 陳金, 張會民. 近35年紅壤稻區(qū)土壤肥力時空演變特征: 以進賢縣為例. 中國農(nóng)業(yè)科學, 2020, 53(16): 3294-3306. doi: 10.3864/j.issn. 0578-1752.2020.16.008.

WANG Y P, HUANG J, SUN Y X, LIU K L, ZHOU H, HAN T F, DU J X, JIANG X J, CHEN J, ZHANG H M. Spatiotemporal variability characteristics of soil fertility in red soil paddy region in the past 35 years—A case study of Jinxian County. Scientia Agricultura Sinica, 2020, 53(16): 3294-3306. doi: 10.3864/j.issn.0578-1752.2020.16.008. (in Chinese)

[17] CARDOSO E J B N, VASCONCELLOS R L F, BINI D, MIYAUCHI M Y H, DOS SANTOS C A, ALVES P R L, DE PAULA A M, NAKATANI A S, DE MORAES PEREIRA J, NOGUEIRA M A. Soil health: Looking for suitable indicators. What should be considered to assess the effects of use and management on soil health? Scientia Agricola, 2013, 70(4): 274-289.

[18] DARILEK J L, HUANG B, WANG Z G, QI Y B, ZHAO Y C, SUN W X, GU Z Q, SHI X Z. Changes in soil fertility parameters and the environmental effects in a rapidly developing region of China. Agriculture, Ecosystems & Environment, 2009, 129(1/2/3): 286-292.

[19] 張祺, 孫萬春, 俞巧鋼, 周家昊, 林輝, 周建利, 馬軍偉. 長期施用有機肥對稻田土壤理化性狀的影響. 浙江農(nóng)業(yè)科學, 2022, 63(7): 1441-1444, 1448.

ZHANG Q, SUN W C, YU Q G, ZHOU J H, LIN H, ZHOU J L, MA J W. Effect of long-term application of organic fertilizer on physical and chemical properties of paddy soil. Journal of Zhejiang Agricultural Sciences, 2022, 63(7): 1441-1444, 1448. (in Chinese)

[20] LUAN H A, GAO W, HUANG S W, TANG J W, LI M Y, ZHANG H Z, CHEN X P. Partial substitution of chemical fertilizer with organic amendments affects soil organic carbon composition and stability in a greenhouse vegetable production system. Soil and Tillage Research, 2019, 191: 185-196.

[21] 成臣, 汪建軍, 程慧煌, 羅亢, 曾勇軍, 石慶華, 商慶銀. 秸稈還田與耕作方式對雙季稻產(chǎn)量及土壤肥力質(zhì)量的影響. 土壤學報, 2018, 55(1): 247-257.

CHENG C, WANG J J, CHENG H H, LUO K, ZENG Y J, SHI Q H, SHANG Q Y. Effects of straw returning and tillage system on crop yield and soil fertility quality in paddy field under double- cropping-rice system. Acta Pedologica Sinica, 2018, 55(1): 247-257. (in Chinese)

[22] 代子俊, 趙霞, 李德成, 劉峰, 石平超, 龐龍輝. 近30年湟水流域土壤全氮時空變異及影響因素. 土壤學報, 2018, 55(2): 338-350.

DAI Z J, ZHAO X, LI D C, LIU F, SHI P C, PANG L H. Spatio-temporal variation of soil total nitrogen in Huangshui River Basin and its affecting factors in the past 30 years. Acta Pedologica Sinica, 2018, 55(2): 338-350. (in Chinese)

[23] 高菊生, 黃晶, 董春華, 徐明崗, 曾希柏, 文石林. 長期有機無機肥配施對水稻產(chǎn)量及土壤有效養(yǎng)分影響. 土壤學報, 2014, 51(2): 314-324.

GAO J S, HUANG J, DONG C H, XU M G, ZENG X B, WEN S L. Effects of long-term combined application of organic and chemical fertilizers on rice yield and soil available nutrients. Acta Pedologica Sinica, 2014, 51(2): 314-324. (in Chinese)

[24] 馬力, 楊林章, 顏廷梅, 王建國, 李斗爭, 殷士學. 長期施肥水稻土氮素剖面分布及溫度對土壤氮素礦化特性的影響. 土壤學報, 2010, 47(2): 286-294.

MA L, YANG L Z, YAN T M, WANG J G, LI D Z, YIN S X. Profile distribution and mineralization characteristics of nitrogen in relation to temperature in paddy soil under long-term fertilization. Acta Pedologica Sinica, 2010, 47(2): 286-294. (in Chinese)

[25] 黃晶, 張楊珠, 徐明崗, 高菊生. 長期施肥下紅壤性水稻土有效磷的演變特征及對磷平衡的響應(yīng). 中國農(nóng)業(yè)科學, 2016, 49(6): 1132-1141. doi: 10.3864/j.issn.0578-1752.2016.06.009.

HUANG J, ZHANG Y Z, XU M G, GAO J S. Evolution characteristics of soil available phosphorus and its response to soil phosphorus balance in paddy soil derived from red earth under long-term fertilization. Scientia Agricultura Sinica, 2016, 49(6): 1132-1141. doi: 10.3864/j.issn.0578-1752.2016.06.009. (in Chinese)

[26] 焉莉, 王寅, 馮國忠, 高強. 吉林省農(nóng)田土壤肥力現(xiàn)狀及變化特征. 中國農(nóng)業(yè)科學, 2015, 48(23): 4800-4810. doi: 10.3864/j.issn.0578- 1752.2015.23.021.

YAN L, WANG Y, FENG G Z, GAO Q. Status and change characteristics of farmland soil fertility in Jilin Province. Scientia Agricultura Sinica, 2015, 48(23): 4800-4810. doi: 10.3864/j.issn. 0578-1752.2015.23.021. (in Chinese)

[27] 王齊齊, 徐虎, 馬常寶, 薛彥東, 王傳杰, 徐明崗, 張文菊. 西部地區(qū)紫色土近30年來土壤肥力與生產(chǎn)力演變趨勢分析. 植物營養(yǎng)與肥料學報, 2018, 24(6):1492-1499.

WANG Q Q, XU H, MA C B, XUE Y D, WANG C J, XU M G, ZHANG W J. Change of soil fertility and productivity of purple soil in Western China in recent 30 years. Journal of Plant Nutrition and Fertilizers, 2018, 24(6): 1492-1499. (in Chinese)

[28] HU W G, JIAO Z F, WU F S, LIU Y J, DONG M X, MA X J, FAN T L, AN L Z, FENG H Y. Long-term effects of fertilizer on soil enzymatic activity of wheat field soil in Loess Plateau, China. Ecotoxicology, 2014, 23(10): 2069-2080.

[29] DUAN Y, CHEN L, LI Y M, WANG Q Y, ZHANG C Z, MA D H, LI J Y, ZHANG J B. N, P and straw return influence the accrual of organic carbon fractions and microbial traits in a Mollisol. Geoderma, 2021, 403: 115373.

[30] 解文艷, 周懷平, 楊振興, 路慧英, 關(guān)春林, 武文麗. 秸稈還田方式對褐土鉀素平衡與鉀庫容量的影響. 植物營養(yǎng)與肥料學報, 2015, 21(4): 936-942.

XIE W Y, ZHOU H P, YANG Z X, LU H Y, GUAN C L, WU W L. Effect of different straw return modes on potassium balance and potassium pool in cinnamon soil. Journal of Plant Nutrition and Fertilizer, 2015, 21(4): 936-942. (in Chinese)

[31] 趙秀娟, 任意, 張淑香. 25年來褐土區(qū)土壤養(yǎng)分演變特征. 核農(nóng)學報, 2017, 31(8): 1647-1655.

ZHAO X J, REN Y, ZHANG S X. Evolution characteristics of cinnamon soil nutrients in 25 years. Journal of Nuclear Agricultural Sciences, 2017, 31(8): 1647-1655. (in Chinese)

Spatial-Temporal Variability Characteristics and Main Controlling Factors of Soil Fertility in Shouyang County of Cinnamon Soil Area

SHEN Tong1, 2, WANG HengFei1, 2, DU WenBo3, ZHOU HuaiPing1, WANG Rui3, ZHANG JianJie1, JIN DongSheng1,2, XU Minggang1, 2

1College of Resources and Environment, Shanxi Agricultural University, Taigu 030800, Shanxi;2Academy of Eco-Environment and Industrial Technology, Shanxi Agricultural University/Shanxi Province Key Laboratory of Soil Environment and Nutrient Resources, Taiyuan 030031;3Shanxi Farmland Quality Monitoring and Protection Center, Taiyuan 030001

【Objective】The main controlling factors of soil fertility and temporal and spatial evolution characteristics in typical county of cinnamon soil region were investigated, so as to provide the basis for scientific management of nutrients and soil fertilization. 【Method】Based on the data of soil fertility, including soil organic matter (SOM), total nitrogen (TN), available phosphorus (AP), available potassium (AK) and pH, at Shouyang County in 1983, 2007 and 2017, the Nemerow index method was used to calculate the comprehensive index of soil fertility, the random forest method was used to explore the main controlling factors of soil fertility, and the combined method of GIS and geostatistics was used to determine the temporal and spatial variability characteristics and the distribution patterns of the comprehensive index of soil fertility and its main factors. 【Result】(1) The overall level of cultivated land soil fertility showed an upward trend in Shouyang County during the past 34 years, and the comprehensive index of soil fertility increased by 0.26, which changed from 1.16 to 1.42. The changes of soil fertility index were characterized by stages. From 1983 to 2007, the average annual increases of SOM, TN and AK were 0.09 g·kg-1, 0.0021 g·kg-1and 1.61 mg·kg-1, respectively, while the changes of AP were not significant. From 2007 to 2017, SOM, TN and AP increased significantly, with the average annual increases of 0.25 g·kg-1, 0.01 g·kg-1and 0.31 mg·kg-1, respectively, while the changes of AK were not significant. (2) The analysis results of the random forest model showed that, the main controlling factors of soil fertility were TN and SOM in 1983, with the importance of 75.3% and 17.8%, respectively. In 2007, the main controlling factors of soil fertility became AK, AP and TN, with the importance of 31.8%, 27.1% and 26.8%, respectively. In 2017, the main controlling factors were TN, AK and AP, and the importance for the three factors were 31.8%, 27.1% and 26.8%, respectively. (3) There were certain spatial differences in the main controlling factors of soil fertility. From 1983 to 2007, SOM increased in the whole county; TN decreased in the northwest and central regions, but increased in other regions; AP decreased in the north regions, but increased in the south regions; AK increased in the whole county. From 2007 to 2017, SOM increased fast in the southeast regions but slow in other regions; TN decreased in the middle regions but increased in other regions; AP decreased in the east regions, but increased in other regions; AK decreased in the east and west regions, but increased in the central region. 【Conclusion】After 34 years, the cultivated land soil fertility of Shouyang County has been improved, and SOM, TN, AP and AK were the main controlling factors affecting the soil fertility variability. SOM, TN and AP increased fast in the south regions, while AK increased fast in the east regions. It was suggested that the whole Shouyang County still needed to increase the application of nitrogen fertilizer moderately, especially in the central region, and stabilize the application of phosphate fertilizer, while control the application of potassium fertilizer in the future.

soil fertility; spatial-temporal variability; cinnamon soil area; main controlling factors; Shouyang County

2022-12-18；

2023-03-28

科技基礎(chǔ)資源調(diào)查專項（2021FY100501）、山西省科技合作交流項目（202104041101002）、山西省博士畢業(yè)生來晉獎勵科研項目（SXBYKY2022083）

申桐，E-mail：3211419826@qq.com。通信作者徐明崗，E-mail：xuminggang@caas.cn

（責任編輯 李云霞）