鄂西北耕層土壤全氮含量空間分異特征及其影響因素分析

木合塔爾·艾買提 周勇 吳正祥

摘要：旨在分析鄂西北耕地表層土壤全氮含量的空間變異特征及其影響因素。采用地統計學、逐步回歸方法和GIS技術。研究結果表明：（1）十堰市耕層土壤的全氮平均含量為1.12 g/kg，耕層土壤全氮含量較高，701個土壤樣點全氮含量變異系數為0.49，屬于中等變異程度。（2）土壤全氮的空間變異最佳模型為高斯模型，塊金效應值為48.09%，全氮具有較強的空間變異性特征，且受隨機性因素和結構性因素的共用影響。（3）十堰市耕層土壤全氮含量在空間分布上具有明顯的差異，全氮含量高的方向為南、西、西南和東南，全氮含量低的方向為東、東北和西北，全氮含量在南北方向的變異性最明顯，差別較大。（4）全氮與有機質、海拔、有效土層厚度有正相關關系，與土壤pH有負相關關系。在黏土、棕壤、水田、河湖沖積物、第四紀老沖積物、一年一熟、種植水稻的土壤中全氮含量較多，在砂壤、紫色土、紫色巖類風化物和種植水果、小麥、噴灌方式下的耕地的全氮含量較少。（5）10種因素對全氮的影響程度排序為土壤環境>成土母質>熟制>生物多樣性>農作物種類>耕層質地>土類>地形>耕地利用方式>灌溉方式，其中前5個因素是影響全氮空間分異的主導因素。在今后的農業生產活動中，需要在河谷和山谷帶、一年二熟區、紅砂巖類風化物和紫色巖類風化物帶、有機質含量較少區可適當地增加氮肥的施用量，需保護生物多樣性并重視人為因素對耕地質量的影響，采用農作物輪作等方式提高耕地質量和耕地生產能力。

關鍵詞：鄂西北；耕地表層；土壤全氮；半變異函數；全氮空間分異；方差分析；逐步回歸

中圖分類號：S159.2文獻標志碼：A論文編號：cjas2020-0103

Total Nitrogen in Plough Layer Soil in Northwest Hubei： Spatial Variability and Influencing Factors

Muhtar·AMAT1， ZHOU Yong2， WU Zhengxiang2，3（1College of Life and Geography Sciences， Kashi University， Kashi 844006， Xinjiang， China； 2College of Urban&Environment Sciences， Central China Normal University， Wuhan 430079， Hubei， China； 3School of Environmental Science and Tourism， Nanyang Normal University， Nanyang 473061， Henan， China）

Abstract： The objective is to understand the spatial variation of total nitrogen and its main influencing factors in plough layer soil of northwest Hubei Province. The geostatistics， stepwise regression method and GIS technology were used in the study. The results show that： （1） the average nitrogen content of Shiyan City is 1.12 g/kg， and the total nitrogen content of plough layer soil is relatively high， with a variation coefficient of 0.49 among the 701 sampling sites， the degree of variation is moderate； （2） the best model for spatial variability of soil total nitrogen is Gaussian model， and the nugget effect value is 48.09%； total nitrogen has strong spatial variability and is affected by the sharing of random factors and structural factors； （3） the total nitrogen content of the plough layer soil in Shiyan City has obvious difference in spatial distribution， areas in the south， west， southwest and southeast have high total nitrogen content， but areas in the east， northeast and northwest have low total nitrogen content， the north-south direction variation is the most obvious and the difference is large； （4） total nitrogen has a positive correlation with organic matter， elevation， and effective soil thickness， and has a negative correlation with soil pH； there is more total nitrogen in plough layer soil such as clay， brown soil， paddy fields， alluvial deposits of rivers and lakes， quaternary alluvial deposits， double cropping and rice planting soils， there is less total nitrogen in sandy soil， purple soil， purple rock weathering materials， fruit planting soil， wheat soil， and soil under sprinkler irrigation； （5） the influence degree of ten factors on total nitrogen is ranked as soil environment > soil parent material > maturity > biodiversity > crop species > soil texture > soil type > terrain > land use mode > irrigation method， and the first five factors are the dominant factors affecting the spatial differentiation of total nitrogen. In future agricultural production， it is necessary to appropriately increase the application amount of nitrogen fertilizer in the river valley and mountain valley， the double cropping area， the red sandstone weathering and the purple rock weathering zone， and the areas with less organic matter content area. Biodiversity needs to be protected and the impact of human activities on soil quality should be highlighted， and crop rotation could be adopted to improve the quality and production capacity of cultivated land.

Keywords： Northwest Hubei；Plough Layer Soil； Soil Total Nitrogen； Semi-variogram； Total Nitrogen Spatial Differentiation； Variance Analysis； Stepwise Regression

0引言

土壤的主要作用是為植物提供所需的養分，植物吸收所需養分多少或植物生長情況由土壤養分元素的含量所決定。土壤養分元素中有氮、磷、鉀為三大元素，其中土壤氮素是決定土壤質量的關鍵指標之一[1]。土壤養分含量對植物生長和糧食生產有較大作用，研究土壤養分元素的區域空間變異特征對制定科學的土壤改良和土地整治方案，為保障國家或區域糧食安全和推動測土配方施肥工作具有重要的現實意義[2]。自20世紀下葉以來，包含對全氮在內的土壤養分元素空間變異的研究取得了一定的成果，研究主要集中在土壤中養分元素含量及土壤物理化學特征空間差異和影響因素等內容[3-5]，研究涉及到的地形區不斷拓展，例如西南喀斯特、山地丘陵、黃土高原[6-10]，研究方法也在不斷地更新，如克里金、地統計學、BP神經網絡[11-15]等。耕地質量或土壤元素含量空間變異研究方法不斷更新，研究結果更符合實際，與實際問題更為貼切，在農業生產和耕地質量保護中發揮較好的作用。在土壤質量和能量的空間分布特征研究中統計學和GIS技術相結合方法已成為主流并較成熟[16-19]。

湖北省是中國重要的糧食生產基地，而鄂西北是地形較復雜的丘陵山地區，土壤所處的地形及其利用方式復雜多樣。丘陵山地區的耕地的全氮含量容易受到地形、氣候、植被、成土母質、土壤類型等自然要素和包含灌溉方式、交通、生產方式等在內的人為因素的影響。在復雜的地形條件下研究鄂西北土壤耕地全氮空間分布規律及其形成機制，在耕地施肥、農業生產規劃、生態環境平衡等方面具有較重要的實踐意義。

本研究以湖北省西北區的十堰市為研究區，以701個樣品點的原始化驗數據為基礎，運用GIS、統計學和地統計學相結合的方法，對土壤全氮含量進行統計學特征分析、空間分布變異分析和影響因素分析，并總結其空間變異規律，以期為找出影響全氮分布的因素和提高十堰市耕地質量提供決策依據。

1數據來源

結合研究區的地形地貌、土壤類型、土地利用方式等信息，遵循代表性、均衡性和可達性等原則選擇采樣點，其空間分布如圖1所示。對十堰市各縣區土壤樣品的采集在2017年農作物秋收后土壤施肥前進行。每個采樣點均在其5 m范圍內取5點制成一個混合土樣，每個土樣重量為1～1.5 kg，表層土樣的采集深度為30 cm土層。在采樣的同時記錄每個采樣點的地理坐標、土地利用類型等20余項信息，一共采樣701個。取好的土樣在登記編號后，再經過風干、研磨、過篩后取500 g樣品，用半微量凱氏定氮法來測定全氮的含量[20-21]。在土樣采集時利用手持式GPS記錄海拔數據；利用研究區的DEM數據在ArcGIS軟件平臺下計算獲得；有效土層厚度、農作物種類、灌溉方式等因子數據也是現場進行記錄；耕層質地、成土母質，土類，土地利用類型，熟制等數據在實地調研的基礎上，結合研究區實際土壤類型和氣象氣候數據、湖北省土壤分類和湖北省土種志的相關指標進行判別而獲得。

2數據處理和研究方法

利用統計學方法分析十堰市耕層土壤全氮含量的描述性特征；采用地統計方法分析其空間變異特征；利用簡單克里金法進行空間插值，并分析空間分異情況；運用方差分析法分析各分類變量的全氮含量，運用逐步回歸判別影響因子對全氮含量變異的解釋能力以及判定主導影響因子。在逐步回歸時對成土母質、生產方式等分類變量采用均值法進行賦值。

在數據處理和分析中，采用SPSS 22.0軟件對全氮含量進行描述性統計分析，用K-S檢驗法在非參數單樣本檢驗（a=0.05）下進行正態性檢驗，使用相關分析功能分析全氮與土壤有機質等定量數據關系與分類變量的關系分析時使用方差分析和逐步回歸分析法。

在描述性統計分析中得出的平均值反映全氮含量的總體情況，方差和標準差反映全氮含量離散度情況，變異系數能反映隨機變量的相對波動程度[22-23]。

利用GS+9.0軟件平臺對十堰市耕層土壤全氮含量地統計學特征進行分析。半方差函數是地統計學的一個分支，主要研究區域化變量的空間變異特征和影響因素的作用等情況。半方差函數是以區域化變量理論為基礎、以變異函數為主要工具，研究空間要素的分布情況的模型，現已經成為土壤學研究的一個重要工具。塊金系數[C0/（C0+C）]可以反映樣本間的空間相關性特征。如果塊金系數大于75%，表明空間異質性主要受隨機性因素影響；如果塊金系數小于25%，表明空間異質性主要受結構性因素影響；若塊金系數介于25%～75%之間，表明區域化變量的空間異質性是隨機因素和結構性因素共同作用的結果[24-25]。運用半方差函數時，一般考慮區域化變量的各向同性特征和各向異性特征。

利用ArcGIS軟件判別半變異理論模型和空間變異性[26-29]，進行插值分析并制作全氮空間分異圖，分析耕層土壤全氮含量空間分布情況。

在研究一個或多個分類型自變量與一個數值型應變量之間的關系時，方差分析是主要的方法之一[30-32]。他主要研究分類型自變量對數值型應變量的影響，他們之間的關系和關系強弱等情況，通過方差分析能夠說明不同類型土壤之間養分的差異情況。方差分析基本條件：各樣本是相互獨立的隨機樣本，均服從正態分布和各樣本總體方差相等，即滿足方差齊性，為此采用F檢驗法。

采用逐步回歸分析方法，分析影響因子對全氮的相對重要性，為此根據研究目的和研究區實際情況，考慮數據的獲得性和完整性等原則，選取土壤環境、耕地利用類型等十大類影響因子組成的影響因子體系。自然因素是土壤生成和發展的決定性因素[33-34]，為此在影響因素中設有土壤環境和熟制等八大類自然因子。隨著人類技術水平的提高，人對土壤環境的影響程度越來越高，除了為達到高產大量使用化肥更改土壤成分之外，不同的農作物和不同的灌溉方式也會影響土壤質量。湖北省十堰市耕層土壤全氮空間變異影響因子體系和數據類型如表1所示。

3結果與分析

3.1土壤全氮含量的統計特征

從計算結果可得知，十堰市耕層土壤全氮含量在0.126～3.740 g/kg之間，其平均值為1.12 g/kg，表明整個十堰市耕層土壤的全氮含量比較豐富，標準差和方差分別為0.55和0.31，表明每一個樣點全氮含量與平均含量之間的離散程度較小。全氮含量的變異系數為0.49，變異程度屬于中等變異。K-S正態分布檢驗結果可表明，十堰市耕地表層土壤全氮含量數據不服從正態分布，但對原數據進行對數變換后服從正態分布。對全氮地統計學分析和進行空間插值時，可對數據進行對數變換，從而保證地統計學分析的正態性。

3.2全氮半方差函數分析

3.2.1各向同性分析各向同性分析在研究區劃變量空間分異情況時不區分方向，主要考慮空間距離來判別區域化變量的空間變異情況和影響因素。利用GS+ 9.0軟件，對十堰市耕層土壤全氮進行半方差函數分析得到各向同性半方差模型（圖2a）。土壤全氮最優半變異理論模型為高斯模型，模型決定系數（R2）為0.982，擬合誤差值（RSS）為8.43′10-4，模型擬合程度較高，能夠有效地反映研究區耕層土壤全氮含量空間變異特征。隨機因素對全氮空間變異的影響程度為16.4%，全氮在研究區域內的變異程度為34.1%，隨機變量比重為48.09%，表明研究區內全氮的空間變異由隨機因素和結構性因素的共同作用下而形成，人為因素對全氮含量的空間變異性的作用小于自然因素；十堰市耕層土壤全氮的變程57 km，該值小于間隔距離（105.75 km），表明全氮含量存在較強空間自相關性。

3.2.2各向異性分析利用半方差函數判別全氮在各方向上的變異性情況。主方向0°表示正北-南（N-S）方向，45°表示東北-西南（NE-SW）方向，90°表示東-西（EW）方向，135°表示東南-西北（SE-NW）方向，15°—35°是在東北-西南（NE-WS）方向。研究區地形地貌較為復雜，地形對全氮的空間分布有一定影響，為了解不同方向上的空間變異情況，繪制在0°、45°、90°、135°方向的半方差函數圖（圖2b），土壤全氮含量在各方向上均有變異性，間隔距離越大，變異性越大，這符合區域化變量的空間相關性特性。從圖2b可知，間隔距離31.80 km之內，N-S方向的變異程度較大，變程距離在38.73～59.98 km之間NE-WS方向上的變異較大，間隔距離66.95 km后，N-S方向的變異又超過其他方向的變異。E-W和SE-WN方向的變異總低于N-S方向的變異，其中SE-WN方向的變異在間隔超過74.02 km起開始下降趨勢。各向異性分析的半方差函數和相應參數見表2。

由表2可知，在各方向上，南北方向的變異最優半方差理論模型為線性模型，東北-西南和東-西方向的最優半方差理論模型為指數模型，東南-西北方向的最優半方差理論模型為球狀模型。塊金系數在25.91%～ 28.57%之間，各方向上的變異均受隨機因素和結構性因素共同影響。影響土壤全氮空間分布情況的隨機性因素有耕作制度、施肥、種植作物、管理水平、環境污染、灌溉能力等。結構性因素對十堰市耕層土壤全氮含量空間變異影響較大，約在72%～75%，影響全氮的結構性（自然因素）因素包括母質、地形、生物多樣性、土壤類型、地質地貌、土壤質地、氣候、生物等。土壤全氮含量半方差模型是球狀模型和指數模型，半方差模型的決定系數均0.90以上，RSS值較小，表明模型可信度較高，均可以采用這些模型模擬預測全氮的空間變異情況。在表2中，正東-西方向和東南-西北方向上的主軸和亞軸比值較低，均等于1，表明該方向上全氮的變異性較小，正北-南方向上的主軸和亞軸比值為1.49，表明該方向上耕層土壤全氮的變異性較大，在不同位置上全氮含量區別較大。分析各向異性而得到的這個結果在研究區域化變量空間變異性有一定的參考價值。

3.3全氮空間分布格局

運用ArcGIS軟件的地統計功能，通過簡單克里金法對十堰市全氮進行空間插值，差值之前進行空間趨勢分析（圖3a），掌握空間變異情況。在空間插值時對全氮原始數據進行對數變換，選擇高斯函數進行插值，并按照國家土壤養分級別標準[35-36]繪制含量插值圖（圖3b）。

由圖3可知，十堰市耕層土壤全氮的分布總趨勢是南高，北低，西高，東低，西南高，東北低，東南高西北低，各方上南北方向的變異性最明顯，差別較大。

經過簡單克里金插值后，根據土壤養分級別標準，十堰市的土壤全氮可分為6個級別：全氮含量大于2 g/kg為極豐富；1.5～2 g/kg之間為豐富；1～1.5 g/kg之間為最適宜；0.75～1 g/kg之間為適宜；0.50～0.75 g/kg之間為缺乏；小于0.5 g/kg為極缺乏。由圖3b可看到，十堰市耕層土壤全氮含量不存在小于0.5 g/kg的情況，缺乏等級主要分布在西北部，全氮含量豐富和極豐富級別主要分布在南部和中部，面積大約占全區域的2/3。

3.4影響因素分析

由表1可知，影響全氮的因素有土壤環境和熟制等10大類，其中土壤環境和地形因素是定量變量，其余因素均為分類（名義）變量。采用相關分析法研究定量因素和全氮含量之間的關系，用單因子方差分析法研究名義變量與全氮之間的關系。

3.4.1土壤環境的影響環境因素和地形因素的數據均為定量數據，可直接用相關分析法得到與全氮的之間的關系，在SPSS軟件平臺下算出的全氮與環境、地形因子各指標之間的皮爾遜相關系數，如表3所示。

由表3可知，土壤全氮與有機質、有效土層厚度、土壤容重之間存在正相關性，其中與有機質和有效土層厚度相關性較顯著，但與土壤容重之間相關性不顯著。全氮與有機質的關系較密切，與有效土層厚度和土壤容重之間的關系不太密切。全氮與土壤pH之間呈弱負相關關系。

圖3和圖4進行比較可以發現，十堰市耕層土壤有機質含量空間分布與全氮空間分布基本相同，即南部高，中部北部相對較低，表明通過提高耕地土壤有機質含量可以提高全氮含量。有效土層厚度較深的南部河谷平原帶和山谷地帶土壤全氮含量較高，表明流水和風流的堆積作用對耕地土壤全氮含量有一定的積極作用。

3.4.2地形的影響利用十堰市DEM（30 m′30 m）數據，運用ArcGIS空間分析和柵格計算功能，分別計算得到701個樣點的坡度、坡向等地形特征值。由表3可知，全氮與海拔、坡度、地形起伏度、地面粗糙度呈正相關關系，與坡向之間有負相關關系，其中與海拔之間的相關性相對較強，在0.01水平下顯著，其余要素之間的相關性較弱，均未通過顯著性檢驗。

十堰市海拔較低、地形平坦的河谷帶耕層土壤全氮含量較高，其原因很可能是高山區土壤全氮受到水土流失的影響發生土壤氮素流失，堆積在地形較平淡和海拔較低的河谷帶，為此呈現出全氮和海拔呈弱正相關特征。

3.4.3生物多樣性的影響土壤環境因素中生物多樣性是定性要素（名義變量），采用單因素方差分析法計算土壤全氮與生物多樣性之間的關系。計算結果顯示，生物多樣性為豐富的107個樣點的全氮含量大約在0.150～2.946 g/kg之間，均值為1.65 g/kg左右；生物多樣性為一般的樣點個數為465個，其全氮含量在0.126～3.740 g/kg之間，均值為1.104 g/kg左右；生物多樣性不豐富的129個樣點全氮含量均值大約為0.73 g/kg左右。3種生物多樣性類型的土壤全氮含量均值排序為豐富>一般>不豐富，生物多樣性豐富時全氮含量最小值均大于其余兩類，變異程度最小。生物多樣性一般的樣點均值、最小值和最大值均大于不豐富類型的相應值，由此可得出，土壤生物多樣性與全氮含量成正比關系。

3.4.4耕層質地的影響耕層質地是土壤生成和發展的一個結果，是影響土壤全氮含量的主要結構性因素之一，對土壤全氮含量和空間分布的影響較明顯。十堰市耕地土壤樣點中有黏土和重壤等5種耕層之地類型，耕層質地與全氮之間的單因子方差分析結果如表4所示。

由表4均值可判斷，輕壤、中壤、重壤的全氮含量在1.1 g/kg左右，各類最小值和最大值之差也不大，砂壤中全氮含量最小為0.9 g/kg左右，黏土中全氮含量最多，大約在1.3 g/kg左右，最小值和最大值差別不大。通過檢驗值可判斷，土壤質地類型對全氮空間分布有明顯的影響。這些耕層質地類型中砂壤中的全氮含量變異較大為52.8%，表明砂壤對土壤全氮的作用較明顯，與其他土壤質地類型比較有顯著的差異。

3.4.5土類和成土母質的影響由表5可知，棕壤、潮土、水稻土、黃棕壤中全氮含量較高，其中黃棕壤中全氮含量大約為1.6 g/kg左右，潮土中的全氮含量比紫色土大約多1倍。棕壤的全氮含量最大，潮土、水稻土的全氮含量在1.2～1.3 g/kg左右，石灰（巖）土，黃褐土和紫色土的全氮含量較少。7種土壤類型中全氮的含量有明顯的差別，平均全氮含量高低依次排列：潮土>水稻土>棕壤>石灰巖土>黃棕壤>黃褐土>紫色土。土壤類型對全氮的影響有明顯的差異性，其中潮土與紫色土，黃棕壤與水稻土、紫色土、棕壤，水稻土與紫色土和黃褐土，石灰巖土與紫色土，紫色土與棕壤，黃褐土與棕壤之間差別較顯著。

十堰市的耕地土壤樣點所屬于的成土母質主要有第四紀老沖積物和紫色巖類風化物等8種類型。潮土的成土母質為河湖沖沉積物，主要是河流和湖泊水堆積作用而形成的，全氮含量最高1.345 g/kg左右，黃褐土的成土母質為第四紀老沖積物其全氮含量為1.003 g/kg，黃棕壤的成土母質為石英巖類風化物、泥質巖類風化物和結晶巖類風化物，其含量為1.094～1.370 g/kg之間，石灰巖土的成土母質為碳酸鹽類風化物全氮含量為1.196 g/kg，水稻土的成土母質為較復雜，主要有河湖沖（沉）積物、結晶巖類風化物、泥質巖類風化物和第四紀老沖積物，這些母質物的全氮含量均大于1 g/kg，紫色土的成土母質為紅砂巖類風化物和紫色巖類風化物，全氮含量最?。?.855 g/kg），棕壤的成土母質為泥質巖類風化物其全氮含量為1.094 g/kg。

各類成土母質對全氮含量的影響有明顯的差異性，其中河湖沖積物、結晶巖類風化物、碳酸鹽類風化物的影響較大，均含量與紫色巖類風化物比較多0.508 g/kg、0.468 g/kg、0.371 g/kg，第四紀老沖積物、紅砂巖類風化物、泥質巖類風化物與全氮的關系與前幾個類型對比較弱。各成土母質平均全氮含量高低依次排列為河湖沖沉積物>石英巖類風化物>碳酸鹽類風化物>結晶巖類風化物>泥質巖類風化物>紅砂巖類風化物>第四紀老沖積物>紫色巖類風化物。

3.4.6耕地利用方式的影響十堰市耕地利用方式有水田、水澆地、旱地3種，全氮與利用方式的單因素方差分析結果顯示，全部樣點中水田的樣點數較多為263個，全氮含量在0.147～3.740 g/kg之間，均值為1.32 g/kg左右，全氮含量最高，均大于其余2種利用方式。水澆地樣點較少為14個，全氮含量在0.690～2.50 g/kg之間，均值為1.139左右，旱地的全氮含量在0.126～ 3.310 g/kg之間，均值為1.046 g/kg。3種土壤類型全氮含量區別不大，其中水田的全氮含量較多，旱地全氮含量較少。水田和水澆地主要分布在中部和東北部地勢較平淡和山谷地帶，交通方便、生產條件便利，農民對農田的投入量較高等原因耕地土壤全氮含量較高。3種農田利用方式中水田與旱地之間的差異性顯著，水澆地是一個過渡帶，水澆地與水田和旱地之間的差異性不顯著。

3.4.7熟制的影響熟制是分類屬性要素，是名義變量，土壤全氮與熟制單因素方差分析計算結果顯示，兩種熟制對土壤全氮含量有明顯差異性（P<0.05）。研究結果表明，在一年一熟制度下的土壤全氮含量均值、最小值、最大值均大于一年二熟制度下的區耕地土壤全氮含量，表明一年一熟制度下的土壤中全氮含量較多。一年一熟制度促進土壤全氮含量的集聚，在實際農業生產過程中實施休耕生產方式有利于氮元素集累。

3.4.8生產方式的影響701個樣點包括的種植農作物種類有茶樹、果樹、水稻、玉米等，灌溉方式主要有溝灌、漫灌和噴灌類型。農作物與全氮之間的方差分析計算結果如表6所示。

由表6可知，種植水稻、茶葉、蔬菜、煙草、玉米的耕地土壤全氮含量比較高，平均值均大于1 g/kg，其中水稻地的全氮含量最高為1.382 g/kg，表明這些農作物對土壤全氮含量的影響較大。果樹、花生、小麥農作物耕地的全氮含量平均值均低于1 g/kg，這些農作物主要分布與干旱地區，受地形，交通和澆水等條件的影響，農民對農田的投入量較小，耕層土壤全氮含量較低。農作物的種類主要與地形和氣候等自然要素有關，人為因素的作用未超過自然因素的作用，不同農作物對土壤全氮含量影響有顯著的差異。

灌溉方式中溝灌和漫灌方式的土壤全氮含量平均值較均為1 g/kg左右，最大值分別為2.964 g/kg和3.740 g/kg，噴灌地的全氮含量較小，最大值1.50 g/kg，相當于漫灌地的一半。3種灌溉方式的最小和最大值之間雖然存在較大的區別，溝灌對全氮的影響高于漫灌和噴灌，但方差分析檢驗值為2.77，顯著性P>0.05，結果可信度不顯著。研究結果表明，灌溉方式對土壤全氮含量的影響不太明顯。

3.4.9影響程度差異性比較為了進一步說明各因子對全氮影響程度大小，使用逐步回歸分析方法，確定各因子對土壤全氮的影響程度[35-36]。由表7可知，各影響因子對土壤全氮含量變異的獨立解釋能力不相同，大約在2%～61.5%，其中土壤環境對全氮的空間分異的影響最大，能解釋空間變異的61.5%，灌溉方式對全氮的影響程度最小，大約為2%。這些因素對全氮變異的單獨解釋能力排序為土壤環境>成土母質>熟制>生物多樣性>農作物種類>耕層質地>土類>地形因素>耕地利用方式>灌溉方式。

由表7可知，全氮為因變量，影響因子為自變量的線性逐步回歸方程在0.0001水平顯著，表明10種因子對耕層土壤全氮的空間分異影響程度均不相同，其中土壤環境，成土母質、熟制、生物多樣性、農作物種類對全氮的影響較大。對全氮分異影響因子中土壤環境的影響最大，能夠解釋全氮變異的61.5%，對回歸方程的作用也最大，其次為成土母質21.1%，熟制和生物多樣性的解釋能力分別為10.9%和10.2%，農作物種類的解釋能力為7.5%。由此可得出，土壤環境因素對耕層土壤全氮含量的影響最大，農作物對全氮含量也有一定的顯著作用。耕層質地和灌溉方式等因素對全氮含量影響不太顯著。

設置全氮含量為因變量，10種因子為自變量，十堰市耕層土壤全氮含量與各影響因素之間的逐步回歸分析計算結果如表8所示。表12中，R2是模型決定系數，表示回歸模型的可信度，該值越高模型可信度越大，△R2是增加該變量時回歸方程判定系數的增量，該系數取值越大，說明該自變量越重要。

由表8可知，模型5的可信度最高，在0.001水平下顯著，該模型包擴土壤環境的5個自變量，表明10中影響因子中土壤環境、成土母質、熟制、生物多樣性，農作物種類為全氮含量空間變異的主導型因子，其中土壤環境因子的作用最大，空間變異的貢獻度53.3%，農作物種類的貢獻度大約為0.1%，表明各因子對全氮的影響有較大的差別。主導因素中土壤環境因子等為自然因素，農作物種類是人為因素，自然因素的貢獻度均大于人為因素，為此得出，自然因素對全氮的影響程度較高，人為因素的影響程度較小。

4討論與結論

十堰市耕層土壤全氮含量呈中等變異特征，整體上，十堰市土壤全氮含量較高，其空間分布受結構性因素和人為因素的共同影響，其中自然因素的作用大于人為因素，全氮空間分布類型為聚集型分布，在空間分異上存在較強的空間相關性，這個研究結果與于倩茹[37]和吳正祥[38]研究結果一致。在湖北西北部土壤采集點布設和野外采樣時空間距離不能超過57 km。耕層土壤全氮分布總趨勢是南高，北低，西高，東低，西南高，東北低，東南高西北低，各方向上變異性最明顯，差別較大。在北部、尤其是西北地區農田施肥時應適當增加氮肥的施用量。全氮與有機質、有效土層厚度、海拔之間均有顯著的正相關性，與土壤pH有顯著的負相關性。在黏土、棕壤、水田、河湖沖積物、第四紀老沖積物、一年一熟，生物多樣性為豐富、種植水稻的土壤中全氮含量較多，在砂壤、紫色土、紫色巖類風化物，生物多樣性為不豐富和種植水果、小麥，噴灌方式下的耕地的全氮含量較少。10種影響因子中土壤環境等自然因素的作用較大，灌溉方式等人為因素的影響程度較小。這些因子對全氮的影響程度排序為土壤環境>成土母質>熟制>生物多樣性>農作物種類>耕層質地>土類>地形因素>耕地利用方式>灌溉方式。

對湖北省西北部而言，在實際的農業生產中，應根據影響養分空間變異的主導因素的差異性，因地制宜地進行施肥。耕層土壤有機質含量較高，但區域性特征較明顯，北部，東北和西北等全氮含量較小地區在農業生產中可采取秸稈還田、農作物輪種、提高農田林網化程度等方式適當增加土壤氮素。土壤生物多樣性和有機質與全氮含量關系較密切，直接反映全氮含量變異的70%左右，根據土壤全氮含量情況進行施肥提高土地利用效率。

鄂西北地區河谷和山谷帶、一年二熟區、紅砂巖類風化物和紫色巖類風化物帶、有機質含量較少區土壤氮含量較低，在今后的農業生產活動中可適當地增加氮肥的施用量。自然因素對耕地土壤氮含量及空間分布的主要影響因素，保護生物多樣性的同時重視人為因素對耕地質量的影響，可采用農作物輪作等方式提高耕地質量和耕地生產能力。研究中對土壤氮含量影響因素中自然要素較多，人為因素較少，這對研究結果可能產生一定的影響。隨著技術水平的提高，人為因素對土壤質量的影響不斷提高，今后研究應重視交通、農民素質、政策等人為因素對耕地土壤質量及氮含量及其空間分布的影響。

參考文獻

[1]徐劍波，宋立生，彭磊，等.土壤養分空間估測方法研究綜述[J].生態環境學報，2011，20（2）：1379-1386.

[2]龐夙，李廷軒，王永東，等.土壤速效氮、磷、鉀含量空間變異特征及其影響因子[J].植物營養與肥料學報，2009，15（1）：114-120.

[3]張嬋嬋，張瑞芳，張建恒，等.高陽縣農田土壤速效養分空間變異特征研究[J].中國生態農業學報，2013，21（6）：758-764.

[4]王晶，任麗，楊聯安，等.基于云模型的西安市蔬菜區土壤肥力綜合評價[J].干旱區資源與環境，2017，31（10）：183-189.

[5]張偉，劉淑娟，葉瑩瑩，等.典型喀斯特林地土壤養分空間變異的影響因素[J].農業工程學報，2013，29（1）：93-101.

[6]劉國順，常棟，葉協鋒，等.基于GIS的緩坡煙田土壤養分空間變異研究[J].生態學報，2013，33（8）：2586-2595.

[7]吳昊.秦嶺山地松櫟混交林土壤養分空間變異及其與地形因子的關系[J].自然資源學報，2015，30（5）：858-869.

[8]董林水，宋愛云，周金星.青藏鐵路沿線土壤有機碳和速效養分空間分異特征[J].干旱區資源與環境，2016，30（11）：161-166.

[9]趙明松，李德成，張甘霖，等.江淮丘陵地區土壤養分空間變異特征——以安徽省定遠縣為例[J].土壤，2016，48（4）：762-768.

[10]陳桂香，高燈州，曾從盛，等.福州市農田土壤養分空間變異特征[J].地球信息科學學報，2017，19（2）：216-224.

[11]連綱，郭旭東，傅伯杰，等.黃土高原小流域土壤養分空間變異特征及預測[J].生態學報，2008，28（3）：946-954.

[12]萬家山，吳云志，張友華，等.基于LM-BP神經網絡的耕地土壤養分等級劃分模型─以皖南山區為例[J].中國農學通報，2015，31（26）：255-260.

[13]裴文文，黃智剛.基于BP神經網絡對蔗區土壤養分肥力等級綜合評價[J].南方農業，2013，7（10）：47-50.

[14]江葉楓，饒磊，郭熙，等.江西省耕地土壤氮素空間變異特征及其影響因素[J].長江流域資源與環境，2018，27（1）：70-79.

[15]迪里木拉提·玉蘇甫，木合塔爾·艾買提，玉蘇甫·買買提，等.溫宿縣植棉土壤表層有效性微量元素空間分布特征[J].土壤通報， 2018，49（5）：1115-1122.

[16]張澤宇，李占斌，李鵬，等.漢江小流域土壤氮素空間分布特征及影響因素[J].水土保持研究，2017，24（4）：46-52，58.

[17]張浩，李啟權，唐嘉玲，等.成都平原西部土壤氮素空間變異特征及影響因素[J].環境科學研究，2018，31（2）：295-302.

[18]叢日環，張智，鄭磊，等.基于GIS的長江中游油菜種植區土壤養分及pH狀況[J].土壤學報，2016，53（5）：1213-1224.

[19]賈振宇，張俊華，丁圣彥，等.基于GIS和地統計學的黃泛區土壤磷空間變異—以周口為例[J].應用生態學報，2016，27（4）：1211-1220.

[20]鮑士旦.土壤農化分析（第三版）[M].北京：中國農業出版社，2000： 30-34.

[21]張彬，楊聯安，馮武煥，等.基于改進TOPSIS和COK的土壤養分綜合評價[J].干旱區資源與環境，2016，30（7）：180-185.

[22]羅由林，李啟權，王昌全，等.四川省仁壽縣土壤有機碳空間分布特征及其主控因素[J].中國生態農業學報，2015，23（1）：34-42.

[23]呂真真，劉廣明，楊勁松，等.環渤海沿海區域土壤養分空間變異及分布格局[J].土壤學報，2014，51（5）：944-952.

[24]蔡明，尚濤，謝龍漢.SPSS統計分析與數據挖掘（第二版）[M].北京：電子工業出版社，2014：158-169.

[25]楊艷麗，史學正，于東升，等.區域尺度土壤養分空間變異及其影響因素研究[J].地理科學，2008，28（6）：788-792.

[26]張敏，賀鵬飛，陳偉強.基于GIS和地統計學的土壤養分空間變異分析[J].東北農業大學學報，2010，41（3）：53-58，161.

[27]張志堅，劉苑秋，吳春生，等.基于地統計學和GIS的江西省森林土壤養分空間分布特征[J].水土保持研究，2018，25（1）：38-46.

[26]楊之江，陳效民，景峰，等.基于GIS和地統計學的稻田土壤養分與重金屬空間變異[J].應用生態學報，2018，29（6）：1893-1901.

[28]賈樹海，張琦，孟維軍，等.基于GIS與地統計學土壤養分空間變異特征研究——以遼寧省凌源市6鄉鎮為例[J].水土保持通報， 2009，29（3）：197-201.

[29]徐建華.地理建模方法[M].北京：科學出版社，2010.

[30]王金哲，張光輝，嚴明疆，等.環渤海平原區土壤鹽分分布特征及影響因素分析[J].干旱區資源與環境，2012，26（11）：104-109.

[31]王永壯，陳欣，史奕.農田土壤中磷素有效性及影響因素[J].應用生態學報，2013，24（1）：260-268.

[32]杜佩穎，張海濤，郭龍，等.平原丘陵過渡區土壤有機質空間變異及其影響因素[J].土壤學報，2018，55（5）：1286-1295.

[33]李友宏，董莉麗，王芳，等.寧夏銀北灌區灌淤土營養元素空間變異性研究[J].干旱地區農業研究，2006（6）：68-72.

[34]趙莉敏，史學正，黃耀，等.太湖地區表層土壤養分空間變異的影響因素研究[J].土壤，2008，40（6）：1008-1012.

[35]康文慧，徐進才，張武文，等.內蒙古達拉特旗土地利用方式對土壤養分的影響[J].干旱區資源與環境，2014，28（8）：161-166.

[36]趙明松，張甘霖，李德成，等.江蘇省土壤有機質變異及其主要影響因素[J].生態學報，2013，33（16）：5058-5066.

[37]于倩茹，徐煖銀，孫思琦，等.鄂西北地區土地利用格局及其生態系統服務價值的時空變化[J].生態科學，2019，38（4）：119-128.

[38]吳正祥，周勇木，木合塔爾·艾買提，等.鄂西北山區耕層土壤pH值空間變異特征及其影響因素研究[J].長江流域資源與環境，2020， 29（2）：488-498.