賈朋遠(yuǎn) 林治家 陳珍寶 宋江濤 賀晨騁
(湖南省地質(zhì)調(diào)查所,湖南 長(zhǎng)沙 410014)
土壤作為歷史自然體,受氣候、生物、母質(zhì)、地形、成土?xí)r間等成土因素的影響,具有復(fù)雜性和時(shí)空變異性,其隨空間位置發(fā)生的變化,被稱為土壤的空間變異性。傳統(tǒng)土壤空間變異的研究方法多是定性的描述,易受人為因素干擾,隨著地統(tǒng)計(jì)學(xué)和GIS技術(shù)的出現(xiàn),定量性的研究被廣泛應(yīng)用到土壤特性時(shí)空變化的研究中。國內(nèi)外許多學(xué)者研究表明,地統(tǒng)計(jì)學(xué)法和多元統(tǒng)計(jì)分析法是目前研究土壤特性空間變異最為有效的方法之一,其在土壤物理性質(zhì)和化學(xué)性質(zhì)空間變異、土壤采樣設(shè)計(jì)以及土壤質(zhì)量管理等方面得到了廣泛應(yīng)用。在土壤物理性質(zhì)方面,顧鑫研究了松嫩平原鹽堿化土壤的電導(dǎo)率、容重、孔隙度、含水量等的空間變異特征;楊士凱以杉木林地為研究對(duì)象,分析了不同土壤深度土壤容重的空間變異特征。在土壤化學(xué)性質(zhì)方面,賈魯凈分析了寶雞市農(nóng)耕區(qū)土壤有機(jī)質(zhì)全氮、速效磷、速效鉀的空間分布規(guī)律,通過冗余分析法,從土壤類型、地形、人為因素等3方面探討了自然和人類活動(dòng)對(duì)土壤養(yǎng)分空間變異的影響程度;趙明松運(yùn)用地統(tǒng)計(jì)學(xué)法和GIS技術(shù)研究了徐淮黃泛平原區(qū)表層土壤有機(jī)質(zhì)含量的空間變異特征,利用方差分析和回歸分析定量研究了區(qū)域內(nèi)土壤有機(jī)質(zhì)空間變異的影響因素;肖凱琦以汨羅市1∶5萬土地質(zhì)量地球化學(xué)調(diào)查為基礎(chǔ),分析了土壤的pH值、有機(jī)質(zhì)、全氮、全磷、全鉀等養(yǎng)分要素,采用地統(tǒng)計(jì)學(xué)與GIS相結(jié)合的方法對(duì)土壤養(yǎng)分空間分布及變異規(guī)律進(jìn)行了研究。在土壤采樣設(shè)計(jì)方面,李俊穎研究了宜興市土壤重金屬的空間變異以及縣域尺度的合理采樣數(shù)。在土壤質(zhì)量管理方面,劉文全研究了山東省萊州灣廢棄鹽田復(fù)墾區(qū)空間變異特征,探討了田間管理分區(qū)的精準(zhǔn)劃分方法。
土壤養(yǎng)分是反映耕地生產(chǎn)力和土壤環(huán)境質(zhì)量的重要指標(biāo),其豐缺程度直接決定農(nóng)作物長(zhǎng)勢(shì),進(jìn)而影響糧食產(chǎn)量。新化縣是湖南省重要的糧食產(chǎn)出地,在1∶5萬土地質(zhì)量調(diào)查之前,對(duì)于區(qū)內(nèi)土地地球化學(xué)質(zhì)量現(xiàn)狀缺乏系統(tǒng)的了解。本次研究基于該項(xiàng)目成果,采用傳統(tǒng)統(tǒng)計(jì)學(xué)、地統(tǒng)計(jì)學(xué)和GIS分析方法對(duì)土壤養(yǎng)分元素的空間分布規(guī)律及影響因素進(jìn)行研究,以期為土壤養(yǎng)分分區(qū)管理及合理施肥提供科學(xué)依據(jù)。
新化縣位于湖南省中部偏西,婁底市西部,在雪峰山南麓,跨資水中游兩岸。地理坐標(biāo)為N27°31′35″~28°14′04″,E110°45′58″~111°41′59″。新化縣四周為低山、低中山,中部為資水及其支流河谷,靠近資水河谷海拔較低,地形相對(duì)平緩,遠(yuǎn)離資江河谷,地勢(shì)逐漸抬升,地形高陡。新化縣屬于亞熱帶季風(fēng)氣候區(qū),四季分明,春暖夏熱,年平均氣溫16.8℃,年平均降水量1432.2mm·a-1。
根據(jù)新化縣第二次國土調(diào)查成果,見圖1,新化縣現(xiàn)有水田590.50km2、旱地169.79km2、園地255.25km2、林地2036.63km2,其他地類567.84km2。其中,水田、旱地按照不少于9點(diǎn)·km-2進(jìn)行網(wǎng)格化布樣,由于地塊較為分散,實(shí)際采樣4776個(gè),采樣密度為6.28個(gè)·km-2。園地按照不少于4點(diǎn)·km-2進(jìn)行布設(shè),實(shí)際采樣595個(gè),采樣密度2.33個(gè)·km-2。林地按照1點(diǎn)·4km-2的進(jìn)行網(wǎng)格化布設(shè),實(shí)際采樣393個(gè),采樣密度0.19個(gè)·km-2。
圖1 新化縣土地利用現(xiàn)狀及采樣點(diǎn)位圖
樣品采集按照《土地質(zhì)量地球化學(xué)評(píng)價(jià)規(guī)范(DZ/T 0295-2016)》執(zhí)行,在布設(shè)采樣點(diǎn)周圍20m范圍內(nèi)采用梅花形布點(diǎn)法(平緩區(qū))或蛇形布點(diǎn)法(低山丘陵區(qū))等量采集5個(gè)子樣組成1個(gè)樣品,水田、旱地、林地采集地表至20cm深處的土柱,園地采集深度0~60cm,采集時(shí)去除雜草、草根、礫石、磚塊、肥料團(tuán)塊等雜物。
“湖南省新化縣土地質(zhì)量調(diào)查”項(xiàng)目共采集測(cè)試了5764個(gè)面積性土壤樣品。樣品加工按照《土地質(zhì)量地球化學(xué)評(píng)價(jià)規(guī)范(DZ/T 0295-2016)》執(zhí)行,樣品檢測(cè)由湖南省地質(zhì)調(diào)查所測(cè)試中心完成,檢測(cè)過程依據(jù)《區(qū)域地球化學(xué)樣品分析方法(DZ/T 0279-2016)》執(zhí)行。該項(xiàng)目測(cè)定了pH、有機(jī)質(zhì)、全氮、全磷、全鉀、堿解氮、有效磷、速效鉀、B、Mn、Zn、Cu、Se、Mo、V、Ni、Cd、Hg、As、Pb、Cr等指標(biāo)。
土壤中的有機(jī)質(zhì)和有效態(tài)測(cè)定方法如下。有機(jī)質(zhì):稱取0.20g樣品,加重鉻酸鉀-硫酸溶液、油浴加熱消煮,采用氧化還原容量法(VOL)測(cè)定有機(jī)質(zhì)。堿解氮:稱取1.00g樣均勻地平鋪于擴(kuò)散外室,在土壤外室內(nèi)加1g鋅-硫酸亞鐵還原劑平鋪土樣上,加3mL 20g·L-1硼酸指示劑溶液于擴(kuò)散皿內(nèi)室;在擴(kuò)散皿外室邊緣上方涂堿性膠液,蓋好毛玻璃,慢慢轉(zhuǎn)開毛玻璃一邊,使擴(kuò)散皿的一邊露出一條狹縫,在此缺口加入10.0mL 1.8mol·L-1氫氧化鈉溶液于皿的外室,立即把毛玻璃蓋嚴(yán);將擴(kuò)散皿放在恒溫箱中,于40℃保溫24h;用0.01mol·L-1鹽酸標(biāo)準(zhǔn)溶液滴定內(nèi)室中吸收的氨量,顏色由藍(lán)變紫紅即達(dá)終點(diǎn),計(jì)算堿解氮結(jié)果。速效鉀:稱取5.00g樣品于浸提瓶中,加50mL 1mol·L-1乙酸銨溶液,加塞震蕩30min,用干濾紙過濾,電感耦合等離子體發(fā)射光譜法(ICP-OES)測(cè)定濾液中鉀。有效磷:酸性、中性土壤中,稱取5.00g樣品于浸提瓶中,加50mL氟化銨-鹽酸浸提劑,震蕩5min后過濾,電感耦合等離子體發(fā)射光譜法(ICP-OES)測(cè)定濾液中磷;堿性土壤中,稱取5.00g樣品于浸提瓶中,加100mL 0.5mol·L-1碳酸氫鈉溶液,在25℃溫度中震蕩30min后過濾,電感耦合等離子體發(fā)射光譜法(ICP-OES)測(cè)定濾液中磷。
由于特異值的存在會(huì)對(duì)變異函數(shù)具有顯著的影響,因此計(jì)算變異函數(shù)前需要剔除特異值。本研究采用拉依達(dá)準(zhǔn)則識(shí)別特異值,即樣本平均值±3倍標(biāo)準(zhǔn)差,在此區(qū)間以外的數(shù)據(jù)均定為特異值,分別用正常的最大和最小值代替特異值。
本次研究對(duì)處理后的數(shù)據(jù)進(jìn)行描述性統(tǒng)計(jì)分析,以地統(tǒng)計(jì)學(xué)基本原理和方法為基礎(chǔ),運(yùn)用GS+9.0軟件進(jìn)行半方差函數(shù)分析,以半方差分析結(jié)果中擬合決定系數(shù)R2最大、殘差平方和RSS最小為原則,選取最優(yōu)插值理論模型和參數(shù),運(yùn)用ArcGIS 10.3軟件進(jìn)行克里金插值,獲得研究區(qū)耕作層土壤養(yǎng)分含量空間分布圖。采用SPSS 28.0軟件對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析和相關(guān)性分析,探討研究區(qū)耕作層土壤養(yǎng)分含量的影響因素。
新化縣4種土壤養(yǎng)分含量及評(píng)價(jià)結(jié)果見表1。有機(jī)質(zhì)含量平均值為34.06g·kg-1,最大值72.85g·kg-1,最小值1.95g·kg-1,變異系數(shù)34.60%,有機(jī)質(zhì)一、二級(jí)地塊面積占了全區(qū)評(píng)價(jià)面積的72.76%。堿解氮含量平均值為154.96mg·kg-1,最大值296.24mg·kg-1,最小值15.26mg·kg-1,變異系數(shù)29.62%,堿解氮一、二級(jí)地塊面積占了全區(qū)評(píng)價(jià)面積的92.32%。有效磷含量平均值21.02mg·kg-1,最大值73.77mg·kg-1,最小值0.13mg·kg-1,變異系數(shù)77.74%,有效磷含量中等及以上的地區(qū)占全區(qū)評(píng)價(jià)面積的92.88%。速效鉀含量平均值為81.84mg·kg-1,最大值205.79mg·kg-1,最小值8.68mg·kg-1,變異系數(shù)48.83%,速效鉀含量中等及以下的地區(qū)占全區(qū)評(píng)價(jià)面積的95.68%,新化縣速效鉀大部分地區(qū)都是中等或較缺乏,僅在南部的石沖口鎮(zhèn)局地較為富集。新化縣土壤鉀素含量整體偏低,應(yīng)通過平衡施肥、輪作和休耕等方式改良和培肥土壤。
表1 土壤養(yǎng)分統(tǒng)計(jì)結(jié)果
變異系數(shù)可以反映土壤元素的離散程度,進(jìn)而表示土壤元素分布的均勻程度。變異系數(shù)越大,元素的分布越不均勻,離散程度越大。本次研究將變異程度分為4種類型:CV<30%為均勻分布,30%≤CV<60%為弱分異(中等起伏),60%≤CV<100%為較強(qiáng)分異(較大起伏),CV≥100%為強(qiáng)分異(很大起伏)。參照這一分類標(biāo)準(zhǔn)顯示堿解氮為均勻分布,有機(jī)質(zhì)、速效鉀為弱分異,有效磷為較強(qiáng)分異。
本研究采用偏度—峰度法檢驗(yàn)數(shù)據(jù)是否服從正態(tài)分布。一般而言,偏度值越接近0,數(shù)據(jù)越服從正態(tài)分布;峰度值越接近于3,數(shù)據(jù)越服從正態(tài)分布。從表1可知,研究區(qū)耕作層土壤養(yǎng)分含量均服從正態(tài)分布(或近似服從正態(tài)分布)。
對(duì)土壤養(yǎng)分狀況的常規(guī)統(tǒng)計(jì)分析能夠概括土壤養(yǎng)分的全貌和整體特征,但不能反映局部的變化特征,即只在一定程度上反映樣本全體,而不能定量刻畫土壤養(yǎng)分的結(jié)構(gòu)性和隨機(jī)性。半方差分析能較好地判定土壤元素的空間自相關(guān)性,并通過空間自相關(guān)性的強(qiáng)弱區(qū)分土壤元素結(jié)構(gòu)性來源(如地質(zhì)背景、地形地勢(shì)等)和隨機(jī)來源(各種人為活動(dòng))。
本文運(yùn)用GS+9.0軟件對(duì)4種養(yǎng)分元素進(jìn)行半變異函數(shù)分析。要合理設(shè)置有效滯后距離(Active Lag Distance)和步長(zhǎng)(Lag Class Distance Interval)。有效滯后距離是用來指定計(jì)算自相關(guān)范圍的,一般設(shè)置為最大采樣點(diǎn)間距的1/2~1/3。由于新化縣地形地貌較為復(fù)雜,土地利用現(xiàn)狀底圖的圖斑比較分散,有效滯后距離設(shè)置太小,會(huì)導(dǎo)致一小部分離得較遠(yuǎn)的圖斑搜不到點(diǎn),降低插值的準(zhǔn)確度和可靠度。本次工作有效滯后距離設(shè)置為最大采樣點(diǎn)間距的1/2,通過Arcgis點(diǎn)距離分析功能得知新化縣最大采樣點(diǎn)間距是86373.51m,故有效滯后距離設(shè)置為43186.75m。步長(zhǎng)是用來定義將多少成對(duì)采樣點(diǎn)分組的,本次分析擬設(shè)置15對(duì)采樣點(diǎn),故步長(zhǎng)設(shè)置為有效滯后距離的1/15,即2879.12m。在設(shè)置好有效滯后距離和步長(zhǎng)后,根據(jù)決定系數(shù)R2最大、殘差平方和RSS最小的原則,選出最優(yōu)模型,見表2、圖2。分析結(jié)果表明,土壤有機(jī)質(zhì)、堿解氮、有效磷、速效鉀的最佳擬合模型均為指數(shù)模型,有機(jī)質(zhì)對(duì)應(yīng)理論模型擬合度最高。
表2 土壤養(yǎng)分半方差函數(shù)模型
圖2 新化縣表層土壤養(yǎng)分半方差圖
塊金值用C0表示,反映的是最小抽樣尺度以下變量的變異性及測(cè)量誤差,表示隨機(jī)部分的空間變異性。新化縣土壤養(yǎng)分塊金值均大于0,表明存在最小抽樣尺度下土壤性質(zhì)的變異或儀器測(cè)量誤差引起的變異。塊金系數(shù)是塊金值與基臺(tái)值的比值,表示由隨機(jī)因素引起的空間變異占系統(tǒng)內(nèi)總變異的比例,反映了空間自相關(guān)性的強(qiáng)弱。塊基比越小,說明受結(jié)構(gòu)性因子影響越顯著,也就是空間自相關(guān)性越強(qiáng)。塊基比小于0.25,說明系統(tǒng)具有強(qiáng)烈的空間自相關(guān)性;塊基比大于0.75,說明系統(tǒng)空間自相關(guān)性很弱;塊基比在0.25~0.75,說明系統(tǒng)具有中等的空間自相關(guān)性。土壤有機(jī)質(zhì)、堿解氮、有效磷的塊金系數(shù)分別是0.13、0.14、0.11,呈現(xiàn)極強(qiáng)的空間相關(guān)性,說明這3種養(yǎng)分元素主要受結(jié)構(gòu)因素(如地形地貌、成土母質(zhì)、土壤類型、氣候等)影響。土壤速效鉀的塊金系數(shù)為0.50,呈現(xiàn)中等空間相關(guān)性,說明速效鉀的空間變異受結(jié)構(gòu)性和隨機(jī)性因素共同作用。
變程是由塊金值到達(dá)基臺(tái)值時(shí)的距離,表示在某種觀測(cè)尺度下空間相關(guān)性的作用范圍,在變程范圍內(nèi),采樣點(diǎn)的距離越近,其相似性越大,即空間相關(guān)性越大。當(dāng)已知點(diǎn)與未知點(diǎn)間距大于最大變程時(shí),該點(diǎn)數(shù)據(jù)不能用于內(nèi)插或外推。土壤有機(jī)質(zhì)、堿解氮、有效磷的變程分別為2850m、3540m、3450m,說明控制其空間變異的影響因素作用范圍相近。速效鉀的變程171810m,說明控制其空間變異的影響因素作用范圍較大,可能是由于速效養(yǎng)分受施肥影響大、遷移性強(qiáng)的原因。
根據(jù)半方差函數(shù)模型及參數(shù),通過克里金插值對(duì)4種土壤養(yǎng)分進(jìn)行擬合,獲取研究區(qū)土壤養(yǎng)分空間分布圖,見圖3。結(jié)果顯示,有機(jī)質(zhì)變異性主要體現(xiàn)在南北方向上,呈南高北低的分布格局;堿解氮變異性主要體現(xiàn)在北東—南西方向上,呈現(xiàn)由南西向北東逐漸遞減的趨勢(shì);有效磷變異性主要體現(xiàn)在東西方向上,呈西高東低的分布格局;速效鉀呈中部低、東部高、西部稍高的分布格局。具體分布特征如下。
圖3 新化縣表層土壤養(yǎng)分空間分布圖
有機(jī)質(zhì)分布的總體特征見圖3a。有機(jī)質(zhì)含量的高值區(qū)主要分布于桑梓鎮(zhèn)、曹家鎮(zhèn),主要位于車前江向斜核部,與二疊系龍?zhí)督M、大隆組等煤系地層關(guān)系密切,可以看出明顯伴隨著黑色巖系地層分布。有機(jī)質(zhì)含量次高值區(qū)主要分布在海拔200~300m的低山、丘陵及山前階地等部位,包括石沖口鎮(zhèn)、維山鄉(xiāng)、西河鎮(zhèn)西部等地,與震旦系金家洞組、留茶坡組和寒武系牛蹄塘組等黑色巖系分布密切相關(guān)。另外,有機(jī)質(zhì)含量次高值區(qū)在爐觀鎮(zhèn)、上梅鎮(zhèn)、奉家鎮(zhèn)(紫鵲界梯田)等地有分布,該區(qū)域多為水稻種植區(qū)域,有機(jī)肥投入較多,且當(dāng)?shù)囟嘤薪斩捇靥锏牧?xí)慣,有助于土壤肥力的積累。有機(jī)質(zhì)含量的低值區(qū)主要分布于水車、西河、吉慶、田坪鎮(zhèn)一帶,這些地區(qū)以山地為主,農(nóng)田面積較少,出露巖性為花崗巖、碳酸鹽巖,表層土壤養(yǎng)分水平較低。
堿解氮分布的總體特征見圖3b。堿解氮含量高值區(qū)主要分布在奉家鎮(zhèn)、西河鎮(zhèn)、瑯塘鎮(zhèn)、金鳳鄉(xiāng)。堿解氮含量次高值區(qū)主要分布在石沖口鎮(zhèn)、維山鄉(xiāng)、科頭鄉(xiāng)等地。高值區(qū)、次高值區(qū)地勢(shì)平坦,是主要的農(nóng)作物產(chǎn)區(qū)。氮含量的低值區(qū)主要分布于水車、文田、吉慶、坐石一帶,這些地區(qū)以山地為主,農(nóng)田面積較少,水土流失嚴(yán)重,出露巖性為花崗巖、碳酸鹽巖,土壤養(yǎng)分水平較低。值得一提的是,資水兩岸普遍為低值區(qū),可能跟河流沖刷侵蝕、養(yǎng)分流失有關(guān)。
有效磷分布的總體特征見圖3c。有效磷含量高值區(qū)分布在白馬山巖體周邊以及大熊山林場(chǎng)北部、榮華鄉(xiāng)西北部,可以看出明顯伴隨著震旦系地層分布。這是因?yàn)楹鲜≈饕哿讜r(shí)期為震旦紀(jì)、寒武紀(jì),以震旦紀(jì)陡山沱組、寒武紀(jì)牛蹄塘組為主要含礦層。在東部的泥盆紀(jì)—石炭紀(jì)碳酸鹽巖分布區(qū),有效磷含量明顯低于其他地層出露區(qū)。
速效鉀分布的總體特征見圖3d。速效鉀含量高值區(qū)主要分布在東部吉慶鎮(zhèn)、坐石鄉(xiāng)、溫塘鎮(zhèn)等,主要為泥盆紀(jì)—石炭紀(jì)碳酸鹽巖區(qū)。速效鉀含量次高值區(qū)主要分布在西部奉家鎮(zhèn)、文田鎮(zhèn)、水車鎮(zhèn)等,為花崗巖出露區(qū),主要巖體有水車、望云山巖體等。巖性以酸性二長(zhǎng)花崗巖及中酸性閃長(zhǎng)巖為主,成土母質(zhì)鉀元素含量較高。
通過以上分析,對(duì)研究區(qū)表層土壤4種主要養(yǎng)分元素含量的空間變異特征和分布規(guī)律有了一定了解,但是對(duì)其來源影響因素尚缺乏具體認(rèn)識(shí)。本文參考以往研究結(jié)果,結(jié)合新化縣實(shí)際情況,從成土母質(zhì)、土壤類型、海拔、坡度、pH值等方面分析其對(duì)土壤養(yǎng)分空間變異的影響。研究思路:根據(jù)插值結(jié)果,提取每一個(gè)農(nóng)用地圖斑(共30106個(gè))的養(yǎng)分元素含量,與圖4a(成土母質(zhì)圖)、圖4b(土壤類型圖)、圖4c(海拔圖)、圖4d(坡度圖)、圖4e(pH值圖)疊加,提取相應(yīng)的屬性值,在SPSS 28.0軟件中進(jìn)行統(tǒng)計(jì)分析。
圖4 新化縣土壤養(yǎng)分相關(guān)因子系列圖
成土母質(zhì)是巖石經(jīng)風(fēng)化作用后物理性質(zhì)改變,在地殼表面形成的疏松的風(fēng)化物,是形成土壤的物質(zhì)基礎(chǔ)。因此,成土母質(zhì)對(duì)土壤養(yǎng)分的含量有著或多或少的影響。研究區(qū)土地資源肥沃,受地層巖性、地形地貌以及沉積環(huán)境等因素影響,可以把研究區(qū)表層土壤成土母質(zhì)劃分為石灰?guī)r風(fēng)化物、砂巖風(fēng)化物、板頁巖風(fēng)化物、花崗巖風(fēng)化物、白堊紀(jì)紅土、河流沖積物等6種類型,見圖4a。
研究區(qū)不同成土母質(zhì)的養(yǎng)分元素統(tǒng)計(jì)表明,板頁巖風(fēng)化物有機(jī)質(zhì)、堿解氮、有效磷含量均高于其他成土母質(zhì)類型;河流沖積物有機(jī)質(zhì)、堿解氮含量均低于其他成土母質(zhì)類型,河流沖積物中有效磷、速效鉀含量也相對(duì)較低,見表3。
表3 不同成土母質(zhì)養(yǎng)分元素含量平均值
對(duì)比圖3和圖4a也能得出相似規(guī)律,有機(jī)質(zhì)在石灰?guī)r類母質(zhì)、花崗巖類母質(zhì)上的含量較低,而在板頁巖類母質(zhì)(主體是早寒武世的黑色巖系、晚二疊世的煤系地層)上的含量較高。堿解氮分布規(guī)律與有機(jī)質(zhì)相似。有效磷在碎屑巖(板頁巖、花崗巖)母質(zhì)上的含量較高,在石灰?guī)r類母質(zhì)上的含量較低。速效鉀在河流沖積母質(zhì)上的含量較低,在石灰?guī)r類母質(zhì)上的含量較高,在碎屑巖母質(zhì)上的含量中等。
第二次全國土壤普查顯示,新化縣共有潮土、黑色石灰土、紅壤、紅色石灰土、黃壤、水稻土、紫色土等7大類,細(xì)分為潮土、黑色石灰土、紅壤、紅壤性土、黃紅壤、紅色石灰土、黃壤、漂白型水稻土、潛育型水稻土、淹育型水稻土、潴育型水稻土、石灰性紫色土、酸性紫色土等13個(gè)亞類。
研究區(qū)不同土壤類型的養(yǎng)分元素統(tǒng)計(jì)表明,黃壤有機(jī)質(zhì)、有效磷、速效鉀含量均高于其他土壤類型;黑色石灰土有機(jī)質(zhì)、堿解氮、速效鉀含量均低于其他土壤類型,有效磷含量也很低,見表4。
表4 不同土壤類型養(yǎng)分元素含量平均值
對(duì)比各元素的空間分布圖和圖4b,發(fā)現(xiàn)不同土壤類型對(duì)各養(yǎng)分元素的空間分布影響不同。在水稻土上,有機(jī)質(zhì)、有效磷含量比較低,堿解氮、速效鉀含量中等。在黃壤上,有機(jī)質(zhì)、堿解氮、有效磷含量比較高。紅壤分布在黃壤周圍,對(duì)養(yǎng)分元素的影響規(guī)律不明顯。其他土壤類型分布面積相對(duì)較小,對(duì)養(yǎng)分元素的空間分布影響較小。
本文使用的是ALOS DEM數(shù)據(jù),是日本宇宙航空研究所使用ALOS(Advanced Land Observing Satellite)衛(wèi)星相控陣型L波段合成孔徑雷達(dá)生產(chǎn)的高程數(shù)據(jù),水平及垂直精度可達(dá)12.5m。利用Arcgis 10.3軟件的空間分析工具提取研究區(qū)的海拔分布圖和坡度分布圖,見圖4c、圖4d。
對(duì)海拔與養(yǎng)分含量做Spearman相關(guān)性分析可以發(fā)現(xiàn),見表5,有機(jī)質(zhì)、堿解氮與海拔高度呈負(fù)相關(guān)(R海拔-有機(jī)質(zhì)=-0.194,R海拔-堿解氮=-0.158),有效磷、速效鉀與海拔高度呈正相關(guān)(R海拔-有效磷=0.168,R海拔-速效鉀=0.271)。將圖4c與圖3做對(duì)比分析也可以看出,有機(jī)質(zhì)、堿解氮的低值區(qū)位于海拔較高的水車、吉慶一帶,而這些地區(qū)的有效磷、速效鉀明顯含量較高。4種養(yǎng)分中海拔對(duì)速效鉀的影響最大。
表5 土壤養(yǎng)分含量與海拔、坡度、pH值的相關(guān)性
對(duì)坡度與養(yǎng)分含量做Spearman相關(guān)性分析可以發(fā)現(xiàn),見表5,有機(jī)質(zhì)、堿解氮與坡度呈負(fù)相關(guān)(R坡度-有機(jī)質(zhì)=-0.289,R坡度-堿解氮=-0.240),有效磷、速效鉀與坡度呈正相關(guān)(R坡度-有效磷=0.112,R坡度-速效鉀=0.098)。將圖4d與圖3做對(duì)比分析可以看出,有機(jī)質(zhì)、堿解氮的低值區(qū)位于坡度較高的天門、大熊山林場(chǎng)一帶,而這些地區(qū)的有效磷、速效鉀含量較高。4種養(yǎng)分中坡度對(duì)有機(jī)質(zhì)的影響最大。
新化縣土壤pH值分布的總體特征見圖4e。以圳上鎮(zhèn)—白溪鎮(zhèn)—孟公鎮(zhèn)—爐觀鎮(zhèn)—石沖口鎮(zhèn)形成的弧形為界,弧西為強(qiáng)酸性環(huán)境,弧東主要為酸性—微酸性環(huán)境,局部有弱堿性環(huán)境?;∥饕栽绻派樾紟r、花崗巖為主,弧東以晚古生代碳酸鹽巖夾碎屑巖為主,說明pH值分布與巖性組合密切相關(guān)。這是因?yàn)樵绻派貙影殡S酸性巖漿活動(dòng),土壤多呈強(qiáng)酸性。而晚古生代大面積出露的是灰?guī)r、白云巖等碳酸鹽巖,表層土壤常伴有石灰反應(yīng),多呈中性甚至弱堿性反應(yīng)。
對(duì)pH值與養(yǎng)分含量做Spearman相關(guān)性分析可以發(fā)現(xiàn),見表5,有機(jī)質(zhì)、堿解氮、速效鉀與pH值呈正相關(guān)(RpH-有機(jī)質(zhì)=0.223,RpH-堿解氮=0.017,RpH-速效鉀=0.204),有效磷與pH值呈負(fù)相關(guān)(RpH-有效磷=-0.379)。將pH值分布圖與4種養(yǎng)分元素含量的空間分布圖做對(duì)比分析可以看出,孟公、游家、田坪等pH值較高的地方,有機(jī)質(zhì)、堿解氮、速效鉀含量也較高。而有效磷在強(qiáng)酸性環(huán)境區(qū)值普遍較高,在酸性—微酸性環(huán)境區(qū)值普遍較低,有效磷與pH值呈明顯負(fù)相關(guān)性。
綜合采用傳統(tǒng)統(tǒng)計(jì)學(xué)、地統(tǒng)計(jì)學(xué)和GIS分析方法,研究了新化縣表層土壤中有機(jī)質(zhì)、堿解氮、有效磷、速效鉀等養(yǎng)分指標(biāo)含量空間變異特征及影響因素。結(jié)果如下。
新化縣土壤有機(jī)質(zhì)、堿解氮、有效磷一級(jí)、二級(jí)地塊占比較高,速效鉀中等及以下的地塊占全區(qū)評(píng)價(jià)面積的95.68%,速效鉀大部分地區(qū)是中等或較缺乏。堿解氮變異系數(shù)29.62%,為均勻分布;有機(jī)質(zhì)變異系數(shù)34.60%、速效鉀變異系數(shù)48.83%,為弱分異;有效磷變異系數(shù)77.74%,為較強(qiáng)分異。
4種養(yǎng)分的最佳擬合模型均為指數(shù)模型,速效鉀對(duì)應(yīng)理論模型擬合度最高。有機(jī)質(zhì)、堿解氮、有效磷塊金系數(shù)分別為0.13、0.14、0.11,呈極強(qiáng)的空間相關(guān)性,速效鉀塊金系數(shù)0.5,呈中等空間相關(guān)性。有機(jī)質(zhì)、堿解氮、有效磷的變程分別為2850m、3540m、3450m,控制其空間變異的影響因素作用范圍較小,速效鉀的變程171810m,控制其空間變異的影響因素范圍較大。
有機(jī)質(zhì)變異性主要體現(xiàn)在南北方向上,呈南高北低的分布格局;堿解氮變異性主要體現(xiàn)在北東—南西方向上,呈現(xiàn)由南西向北東逐漸遞減的趨勢(shì);有效磷變異性主要體現(xiàn)在東西方向上,呈西高東低的分布格局;速效鉀呈中部低、東部高、西部稍高的分布格局。
板頁巖風(fēng)化物有機(jī)質(zhì)、堿解氮、有效磷含量均高于其他成土母質(zhì)類型;河流沖積物有機(jī)質(zhì)、堿解氮含量均低于其他成土母質(zhì)類型,有效磷、速效鉀含量也相對(duì)較低。黃壤的有機(jī)質(zhì)、有效磷、速效鉀含量均高于其他土壤類型;黑色石灰土有機(jī)質(zhì)、堿解氮、速效鉀含量均低于其他土壤類型,有效磷含量也很低。有機(jī)質(zhì)、堿解氮與海拔高度、坡度呈負(fù)相關(guān),有效磷、速效鉀與海拔高度、坡度呈正相關(guān)。有機(jī)質(zhì)、堿解氮、速效鉀與pH值呈正相關(guān),有效磷與pH值呈明顯負(fù)相關(guān)性。