馬泉來，王小玉①，趙曼宇，高鳳杰，陳 浩，孫曉兵

(1.河南省資源環(huán)境調(diào)查一院，河南 鄭州 450007；2.東北農(nóng)業(yè)大學(xué)公共管理與法學(xué)院，黑龍江 哈爾濱 150030；3.中國農(nóng)業(yè)大學(xué)土地科學(xué)與技術(shù)學(xué)院，北京 100193)

土壤中的氮(N)、磷(P)元素是作物生長必需的營養(yǎng)元素，也是陸地生態(tài)系統(tǒng)植物生長過程中常見的限制性元素，在陸地生態(tài)系統(tǒng)平衡和物質(zhì)能量循環(huán)等過程中發(fā)揮著重要作用[1-2]。有關(guān)土壤中N、P元素與作物生長關(guān)系、空間分布格局的研究一直都是土壤學(xué)研究的熱點(diǎn)內(nèi)容[3-5]。有研究表明生物體自身和外界環(huán)境中的氮磷比(N/P比)影響生物的相對生長率和發(fā)育程度，與土壤N、P單一元素的變異特征相比，N/P比能夠更科學(xué)有效地闡述土壤N、P的內(nèi)部特征及平衡關(guān)系，對于認(rèn)識植物和生態(tài)系統(tǒng)養(yǎng)分供應(yīng)與需求平衡過程、養(yǎng)分耦合循環(huán)特征及其驅(qū)動機(jī)制更具重要意義[6-7]。生態(tài)化學(xué)計(jì)量學(xué)(ecological stoichiometry)作為一門研究多重生態(tài)過程化學(xué)元素及能量平衡關(guān)系的新興生態(tài)學(xué)科，能夠從元素比例的角度將土壤-作物生態(tài)過程及土壤養(yǎng)分元素循環(huán)等研究統(tǒng)一起來[8-10]。生態(tài)化學(xué)計(jì)量學(xué)的產(chǎn)生和發(fā)展為植物和生態(tài)系統(tǒng)養(yǎng)分供應(yīng)與平衡、土壤養(yǎng)分耦合循環(huán)特征及驅(qū)動機(jī)制方面的研究提供了更為有效的途徑[8,10-11]。當(dāng)前有關(guān)土壤養(yǎng)分的研究多集中于某一行政區(qū)劃范圍內(nèi)的單一或多種元素的研究，具有一定的普遍性和限制性[7-8,11-13]。因?yàn)橐孕姓^(qū)劃為單元進(jìn)行土壤養(yǎng)分元素空間分布及影響因素分析，一定程度上忽略了土壤養(yǎng)分形成及循環(huán)累積的自然過程，也不符合自然要素自然演變的機(jī)理特征，而流域作為一種集社會經(jīng)濟(jì)和自然要素于一體的復(fù)合系統(tǒng)，以流域?yàn)閱卧M(jìn)行土壤元素生態(tài)化學(xué)計(jì)量的空間變異特征分析，可以更好地揭示空間尺度上土壤養(yǎng)分的供需平衡、格局演變等復(fù)雜過程。

該研究以東北黑土區(qū)海溝河小流域表層土壤(0～20 cm)中全氮(TN)、全磷(TP)含量及N/P比為研究對象，綜合運(yùn)用地統(tǒng)計(jì)學(xué)模型和“3S”技術(shù)，對海溝河小流域內(nèi)的土壤N、P含量及其生態(tài)化學(xué)計(jì)量空間分異特征和主要影響因素進(jìn)行量化分析，以期為流域內(nèi)的耕地產(chǎn)能提升、施肥過程管理及農(nóng)業(yè)生產(chǎn)的可持續(xù)發(fā)展提供理論和技術(shù)依據(jù)，促進(jìn)黑土區(qū)糧食產(chǎn)能提升與耕地環(huán)境保護(hù)的協(xié)同發(fā)展。

1 材料與方法

1.1 研究區(qū)概況

海溝河小流域(地理坐標(biāo)為45°34′18″～45°40′50″ N，126°55′45″～127°10′05″ E)位于哈爾濱市阿城區(qū)，以料甸滿族鄉(xiāng)為主，還包括紅星鄉(xiāng)部分村屯，總面積為119.76 km2。流域內(nèi)海溝河屬松花江的二級支流，由東向西經(jīng)阿什河在哈爾濱市東南注入松花江。流域內(nèi)地形變化復(fù)雜，呈東高西低之勢，海拔在54～366 m之間，東部以低山丘陵區(qū)為主，向西逐步過渡到平原耕作區(qū)，耕地占總面積的70.97%，林地占總面積的24%。東部山地丘陵地帶土壤類型以森林暗棕壤為主；西部地勢平坦地帶主要種植水稻，土壤類型以黑土發(fā)育的水稻土為主。流域內(nèi)土地利用方式交錯復(fù)雜，自東向西分別以林地、水田、旱地為主。海溝河小流域地處城市邊緣區(qū)，社會發(fā)展和農(nóng)業(yè)生產(chǎn)活動受城市發(fā)展輻射作用強(qiáng)，人類生產(chǎn)生活導(dǎo)致的農(nóng)業(yè)生態(tài)環(huán)境問題也較為凸顯。

1.2 數(shù)據(jù)獲取及處理

1.2.1基礎(chǔ)地理數(shù)據(jù)

研究采用了Google Earth高分辨率SPOT遙感影像(拍攝日期2014年9月6日，空間分辨率2.5 m，來源于http:∥maps.google.com/)，選取地面控制點(diǎn)在ArcGIS 10.0平臺進(jìn)行幾何校正和配準(zhǔn)，并通過人工目視判讀解譯獲取海溝河小流域土地利用數(shù)據(jù)。結(jié)合野外定點(diǎn)采樣，遙感影像數(shù)據(jù)解譯準(zhǔn)確率達(dá)90%以上，可滿足研究需要。地形數(shù)據(jù)主要從中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心地理空間數(shù)據(jù)云平臺(http:∥www.gscloud.cn)中30 m精度的DEM數(shù)據(jù)提取，土壤數(shù)據(jù)由中國科學(xué)院編制的《中國東北地區(qū)土壤圖》經(jīng)矢量化處理獲得。

1.2.2樣品采集及測試分析

以海溝河小流域內(nèi)地形、土壤、土地利用現(xiàn)狀數(shù)據(jù)為輔助信息，綜合考慮地形、土地利用、居民點(diǎn)分布、交通干道、河流水系等要素影響，進(jìn)行室內(nèi)網(wǎng)格布點(diǎn)，網(wǎng)格間距為900 m。選擇2014年10月農(nóng)作物收割后進(jìn)行野外土壤樣品采集，在清除地表雜草、秸稈、落葉等干擾后，按照五點(diǎn)法和四分法收集耕作層(0～20 cm)土壤樣品120個(gè)，其中玉米地樣品58個(gè)，水稻地50個(gè)，林地10個(gè)，菜地和亞麻地各1個(gè)。同時(shí)記錄采樣點(diǎn)的土地利用類型、坡位、坡向、耕作方式、作物秸稈還田情況等信息。室內(nèi)測定土壤TN、TP含量，具體測定分析方法參照文獻(xiàn)[14]。

1.2.3異常值處理

在考慮野外采樣和實(shí)驗(yàn)室測試分析誤差的基礎(chǔ)上，采用格拉布斯(Grubbs)法進(jìn)行數(shù)據(jù)異常值的識別和處理，剔除異常數(shù)據(jù)，共得到有效數(shù)據(jù)116組，后續(xù)相關(guān)計(jì)算分析均采用處理后的數(shù)據(jù)。

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

地統(tǒng)計(jì)學(xué)中的半方差函數(shù)模型是土壤養(yǎng)分元素空間分異特征研究中的有效方法[15-16]，該研究主要通過半方差函數(shù)模型對土壤TN、TP含量和N/P比的空間變異特征進(jìn)行分析，公式為

(1)

式(1)中，γ(h)為變異函數(shù)；N(h)為分割距離為h時(shí)的樣本點(diǎn)總數(shù)；Z(xi)為Z(x)在空間位置xi處的實(shí)測值；Z(xi+h)為Z(x)在xi處距離h的實(shí)測值。

克里金插值法是自然地理要素空間插值領(lǐng)域應(yīng)用最為廣泛的方法之一。在已有研究的基礎(chǔ)上，通過SPSS 19.0軟件對樣本數(shù)據(jù)進(jìn)行描述性統(tǒng)計(jì)分析、正態(tài)分布檢驗(yàn)和方差分析。通過普通克里金插值法對研究區(qū)土壤TN、TP含量和N/P比進(jìn)行空間分布插值。半方差函數(shù)及理論模型的擬合等計(jì)算過程均采用地統(tǒng)計(jì)學(xué)軟件GS+ 9.0進(jìn)行，空間插值、距離計(jì)算及相關(guān)繪圖在ArcGIS 10.0和Origin 9.0軟件中完成。

2 結(jié)果與分析

2.1 土壤TN、TP含量和N/P比統(tǒng)計(jì)特征分析

對海溝河小流域土壤樣點(diǎn)的TN、TP含量和N/P比進(jìn)行經(jīng)典統(tǒng)計(jì)學(xué)分析，結(jié)果見表1。

表1 土壤TN、TP含量和N/P比的描述性統(tǒng)計(jì)特征

從表1可知，海溝河小流域土壤TN、TP含量和N/P比的范圍分別為0.96～5.11 g·kg-1、0.21～0.92 g·kg-1和1.23～20.19，均值分別為2.60 g·kg-1、0.55 g·kg-1和4.78，參考全國第二次土壤普查的氮素含量分級標(biāo)準(zhǔn)及黑土區(qū)相關(guān)研究[17-20]，研究區(qū)土壤TN含量和N/P比整體處于較高水平，TP含量處于較低水平，流域內(nèi)N、P含量間的平衡關(guān)系可能會成為海溝河小流域糧食產(chǎn)能提升、耕地生態(tài)環(huán)境維持、農(nóng)業(yè)面源污染控制的關(guān)鍵點(diǎn)。變異系數(shù)為<0.1、0.1～1、>1分別表示弱變異、中等變異和強(qiáng)變異[21]，研究區(qū)土壤TN、TP含量和N/P比的變異系數(shù)分別為37.69%、25.47%和56.49%，空間分布均為中等變異程度，說明研究區(qū)人類活動、耕作措施、土地利用等隨機(jī)因素對土壤TN、TP含量和N/P比的空間變異存在一定的影響。K-S非參數(shù)檢驗(yàn)(α=0.05)顯示，研究區(qū)TN含量符合正態(tài)分布，TP含量和N/P比需進(jìn)行對數(shù)轉(zhuǎn)換，然后對樣點(diǎn)數(shù)據(jù)進(jìn)行地統(tǒng)計(jì)學(xué)分析和空間分布插值。

2.2 土壤TN、TP含量和N/P比的空間變異特征分析

通過GS+ 9.0平臺對研究區(qū)土壤TN、TP含量和N/P比進(jìn)行地統(tǒng)計(jì)學(xué)分析，得到各向同性條件下的半方差函數(shù)模型(圖1)及相關(guān)擬合參數(shù)(表2)。

圖1 土壤TN、TP含量和N/P比的各向同性半方差函數(shù)Fig.1 Semi-variogram models for TN,TP contents and N/P ratio in isotropy

模型擬合結(jié)果表明研究區(qū)土壤TN含量和N/P比的半方差函數(shù)理論模型為高斯模型(Gaussian)，TP含量為指數(shù)模型(Exponential)，擬合殘差分別為4.60×10-2、3.683×10-6和2.211×10-4，擬合程度較高，能夠真實(shí)有效地反映研究區(qū)TN、TP含量和N/P比的空間變異特征。

表2 土壤TN、TP含量和N/P比的半方差函數(shù)模型及參數(shù)

塊金值(C0)能夠有效解釋區(qū)域變量在小于抽樣尺度發(fā)生變異的程度。研究區(qū)土壤TN、TP含量和N/P比的塊金值分別為0.002、0.002和0.002，處于極低水平，表明在小于當(dāng)前采樣尺度上幾乎不存在某種生態(tài)過程對其空間分布產(chǎn)生影響，采樣尺度可進(jìn)一步擴(kuò)大。基臺值(C0+C)一般用來表示系統(tǒng)內(nèi)的總變異程度，塊金效應(yīng)〔C/(C0+C)〕則用來表示系統(tǒng)內(nèi)非隨機(jī)部分引起的變異占系統(tǒng)總變異的比例。研究區(qū)土壤TN、TP含量和N/P比的塊金效應(yīng)分別為0.99、0.75和0.95，明顯高于華北等其他平原耕作區(qū)[5，22]，表明研究區(qū)人類開發(fā)利用時(shí)間短，土壤TN、TP含量和N/P比的空間變異受人類活動、農(nóng)業(yè)生產(chǎn)等隨機(jī)因素的影響弱于區(qū)域地形、成土母質(zhì)等結(jié)構(gòu)性因素，這與華北平原、西部黃土區(qū)的相關(guān)研究有較大差異[5，22]。研究區(qū)TN、TP含量和N/P比的變程分別為510、780和590 m，說明研究區(qū)TN、TP含量和N/P比在該尺度內(nèi)存在較強(qiáng)的連續(xù)性，該結(jié)果可以為該區(qū)域內(nèi)地理要素采樣尺度的選擇提供參考。

在以上研究的基礎(chǔ)上，考慮研究區(qū)地處黑土區(qū)邊緣，區(qū)域地形和土地利用情況復(fù)雜，林地、旱地、水田交錯存在，為了更好地揭示區(qū)域內(nèi)TN、TP含量和N/P比的空間分異特征，對研究區(qū)土壤TN、TP含量和N/P比做各向異性分析，繪制0°、45°、90°和 135°方向的半方差函數(shù)圖(圖2)。

圖2 土壤TN、TP含量和N/P比的各向異性半方差函數(shù)Fig.2 Semi-variogram models for TN, TP contents and N/P ratio in anisotropy

圖2表明，研究區(qū)土壤TN、TP含量和N/P比在不同方向上的變異程度存在一定的差異性。當(dāng)間隔距離小于1 km時(shí)，研究區(qū)TN含量在南北和東西方向的變異程度小，當(dāng)間隔距離大于2 km時(shí)，各方向變異程度在0.7～1.2之間，但隨著間隔距離的增加，TN含量變異程度呈上升和趨同的整體趨勢。研究區(qū)TP含量總體變異程度極小，不同方向上的變異程度在間隔距離小于2 km時(shí)存在一定差異，但隨著間隔距離增加，尤其是大于3 km時(shí)，TP含量在不同方向上的半方差函數(shù)特征表現(xiàn)出較強(qiáng)的一致性，這與研究區(qū)TP含量空間特征受地形、成土母質(zhì)等結(jié)構(gòu)性要素影響較大的結(jié)論一致。研究區(qū)N/P比的半方差函數(shù)特征隨著間隔距離的變化波動較大，當(dāng)間隔距離接近和大于6 km后，不同方向上的半方差函數(shù)值趨于一致。綜合各向異性函數(shù)分析結(jié)果，研究區(qū)土壤TN、TP含量和N/P比在不同間隔距離內(nèi)的變異程度存在一定差異性，其中TN含量和N/P比的變異程度較TP更為明顯，這對該區(qū)域后續(xù)研究中采樣點(diǎn)的布設(shè)具有重要的指導(dǎo)意義。

2.3 土壤TN、TP含量和N/P比的空間分布特征分析

在對研究區(qū)土壤TN、TP含量和N/P比進(jìn)行半方差函數(shù)分析的基礎(chǔ)上，分別選擇高斯模型和指數(shù)模型對其進(jìn)行普通克里格插值，得到研究區(qū)土壤TN、TP含量和N/P比的空間分布圖(圖3)。

圖3 土壤TN、TP含量和N/P比的空間分布Fig. 3 Spatial distribution of TN,TP contents and N/P ratio

由圖3可見，研究區(qū)土壤TN、TP含量和N/P比在空間上呈現(xiàn)自東向西階梯狀的變化趨勢，空間變化特征明顯，與研究區(qū)地形在東西方向上的變化具有較高的一致性。南北方向土壤TN、TP含量和N/P比變化較小，這與南北方向地形變化平緩密切相關(guān)。土壤TN含量的高值區(qū)分布在海拔較高的林地向平原耕作區(qū)的過渡地帶以及南部地勢較低的水田集中分布區(qū)域。土壤TP含量的高值區(qū)主要出現(xiàn)在東部山區(qū)和西部河流匯集區(qū)域，主要是土壤侵蝕和沉積作用所致，這與典型黑土區(qū)、滇池大東河流域的研究結(jié)果相近[23-24]。低值區(qū)出現(xiàn)在研究區(qū)中部耕地和林地的交錯過渡區(qū)域，說明流域內(nèi)農(nóng)林交錯地帶土壤磷含量相對不足，主要原因可能是該地帶農(nóng)作物以玉米、大豆等作物為主，地表覆蓋有限，土壤抗侵蝕能力低，易產(chǎn)生土壤侵蝕。受土壤TN和TP含量影響，土壤N/P比最大值出現(xiàn)在研究區(qū)中部的農(nóng)林交錯地帶，該區(qū)域生態(tài)環(huán)境相對脆弱，生態(tài)系統(tǒng)穩(wěn)定性差，應(yīng)作為農(nóng)業(yè)生產(chǎn)及農(nóng)業(yè)面源污染等環(huán)境問題的重點(diǎn)防控區(qū)域。N/P比與TP含量的空間分布特征整體上較為一致，可以看出N/P比的空間分布受TP含量的限制性作用較大，整體呈現(xiàn)出兩邊低、中間高的分布格局。

2.4 土壤TN、TP含量和N/P比空間分布的影響因子分析

地形因素被認(rèn)為是影響土壤理化性質(zhì)空間分布的重要因素，其作用貫穿土壤養(yǎng)分形成發(fā)育的整個(gè)過程，隨著人類對土地資源的進(jìn)一步開發(fā)與利用，在區(qū)域地形、成土母質(zhì)等結(jié)構(gòu)性要素相對一致且影響有限的情況下，人類活動及其所產(chǎn)生的環(huán)境效應(yīng)在土壤養(yǎng)分空間分布及再配置過程中的影響作用也越來越明顯，加強(qiáng)對人類活動、農(nóng)業(yè)生產(chǎn)及管理措施等非結(jié)構(gòu)性因子對小區(qū)域尺度上土壤養(yǎng)分空間分布及耕地系統(tǒng)健康等的影響研究，已經(jīng)成為土壤科學(xué)與農(nóng)業(yè)可持續(xù)發(fā)展研究的重點(diǎn)[25-26]。

筆者參考已有土壤養(yǎng)分分布及影響因子的相關(guān)分析[27-28]，結(jié)合研究區(qū)地形因素(高程、坡度)、居民點(diǎn)、河流水系及道路交通分布，對研究區(qū)土壤TN、TP含量和N/P比的空間分布與環(huán)境因子進(jìn)行相關(guān)性分析(表3)。

采樣點(diǎn)中有2處分別位于菜地及亞麻地，樣點(diǎn)過少，不參與統(tǒng)計(jì)分析。從相關(guān)性分析結(jié)果可以看出，高程、坡度、剖面曲率、平面曲率、復(fù)雜地形因子及水系距離等地形環(huán)境因子與土壤TN、TP含量和N/P比不存在顯著相關(guān)性，這可能跟流域內(nèi)人類活動較為復(fù)雜，對耕地表層土壤中N、P含量及分布產(chǎn)生了一定的干擾有關(guān)。相關(guān)性分析結(jié)果還顯示，距離村莊遠(yuǎn)近與TN含量存在一定的負(fù)相關(guān)性(r=-0.234)，距離村莊近的地方土壤TN含量較高，可能與村莊是人類活動的主要場所，日常生產(chǎn)生活(冬季取暖、廢物堆積)、畜牧生產(chǎn)等行為所產(chǎn)生的氮素在村莊周圍累積有關(guān)。流域內(nèi)土壤TP含量與道路距離呈顯著正相關(guān)(P<0.01)，隨著距離的增加，TP含量呈現(xiàn)增大趨勢。N/P比與道路距離呈顯著負(fù)相關(guān)(P<0.05)，交通干道兩側(cè)土壤N/P比較高。已有研究表明，交通干道、居民點(diǎn)分布等人類活動會直接影響農(nóng)業(yè)生產(chǎn)過程中的施肥管理、灌溉等行為，并且隨著時(shí)間的積累進(jìn)一步影響區(qū)域內(nèi)土壤養(yǎng)分的空間格局，且其影響強(qiáng)度會隨著農(nóng)業(yè)生產(chǎn)過程的延續(xù)而越來越明顯[16]。

表3 土壤TN、TP含量和N/P比與地形及環(huán)境因子的相關(guān)系數(shù)

STI為沉積物運(yùn)移指數(shù)；CTI為復(fù)合地形指數(shù)；SPI為匯流動力指數(shù)。*表示在0.05水平(雙側(cè)) 顯著相關(guān)；**表示在0.01水平(雙側(cè)) 顯著相關(guān)。

研究區(qū)土地利用復(fù)雜，坡度變化較大，土壤類型分布較為單一，以森林暗棕壤和黑土發(fā)育的水稻土為主，因此僅對不同土地利用類型和坡位的土壤 TN、TP含量和N/P比進(jìn)行均值比較分析(圖4)。

圖4 不同土地利用類型、坡位土壤TN、TP含量和 N/P比的比較分析Fig.4 Comparison analysis of mean value of soil TN, TP contents and N/P ratio in different land use and slope position

不同土地利用方式下，土壤TN含量平均值表現(xiàn)為林地>玉米地>水稻地，林地土壤TN含量明顯高于玉米地和水稻地，這可能與林地采樣點(diǎn)所處環(huán)境較為封閉，受人類活動干擾少，且林下植被落葉及微生物等長期作用有關(guān)。TP含量平均值表現(xiàn)為水稻地> 玉米地>林地，這是因?yàn)闁|北地區(qū)水稻地的磷肥施用量高于玉米地，且水稻地一般較為平坦，基本無侵蝕，TP含量通常較高，加之水稻長期處于淹水狀態(tài)，環(huán)境中pH值較高，氧化還原電位低，有利于磷的吸附和固定[16,20,23]。受土壤TN、TP含量影響，流域內(nèi)土壤N/P比平均值表現(xiàn)為林地>玉米地>水稻地。

對流域內(nèi)不同坡位土壤TN、TP含量和N/P比的均值比較分析顯示，TN含量表現(xiàn)為坡頂>坡中>坡底>平地；TP含量表現(xiàn)為平地>坡中>坡底>坡頂；N/P比表現(xiàn)為坡頂>坡中>坡底>平地，坡頂和坡中土壤TN含量整體高于平地，這可能與坡頂、坡中土壤在長期演變和形成過程中更易受到風(fēng)化作用，向坡底運(yùn)移作用有限等有關(guān)。而TP含量則表現(xiàn)為平地高于其他坡位，這是因?yàn)槠降囟嘁运?、玉米種植為主，而水稻地和玉米地受侵蝕作用小，加上農(nóng)業(yè)生產(chǎn)施肥影響，導(dǎo)致其土壤TP含量高于其他區(qū)域。平地TN含量整體上較低，而TP含量高，導(dǎo)致其N/P比達(dá)最高，因此該區(qū)域是施肥管理和農(nóng)業(yè)面源污染防控需關(guān)注的重點(diǎn)區(qū)域。

3 結(jié)論

基于經(jīng)典統(tǒng)計(jì)學(xué)和地統(tǒng)計(jì)學(xué)模型方法，綜合應(yīng)用“3S”技術(shù)，對東北黑土區(qū)海溝河小流域土壤氮磷生態(tài)化學(xué)計(jì)量空間分異特征及其影響因素進(jìn)行量化分析，得出以下結(jié)論：

(1)海溝河小流域土壤TN、TP含量和N/P比的范圍分別為0.96～5.11 g·kg-1、0.21～0.92 g·kg-1和1.23～20.19，均值分別為2.60 g·kg-1、0.55 g·kg-1和4.78，空間變異程度均為中等變異。地形、成土母質(zhì)等結(jié)構(gòu)性因素仍是影響流域內(nèi)土壤TN、TP含量和N/P比空間變異的主要因子。

(2)土壤TN、TP含量和N/P比在空間上呈現(xiàn)自東向西的階梯狀變化趨勢，南北向變化微小，這與研究區(qū)地形東西落差大、南北變化平緩具有較高的一致性。流域中部的農(nóng)林交錯地帶是N/P比高值區(qū)，N/P比整體呈現(xiàn)出兩邊低、中間高的分布格局。

(3)不同土地利用方式下土壤TN、TP含量和N/P比存在較大差異，土壤TN含量平均值為林地>玉米地>水稻地，TP含量平均值為水稻地>玉米地>林地，N/P比平均值為林地>玉米地>水稻地，耕地N/P比整體偏低。

(4)農(nóng)業(yè)管理措施、農(nóng)村居民點(diǎn)、交通道路等對土壤養(yǎng)分空間分布的影響開始凸顯，該類因素是未來土壤養(yǎng)分空間分布及耕地系統(tǒng)健康等方面的重要研究內(nèi)容。