高 琳，陳曉遠(yuǎn)，林昌華，王 衛(wèi)，張宇鵬(韶關(guān)學(xué)院英東農(nóng)業(yè)科學(xué)與工程學(xué)院，廣東 韶關(guān) 512005)

【研究意義】土壤肥力是影響農(nóng)作物生長(zhǎng)發(fā)育和產(chǎn)量最為重要的因素之一，而土壤養(yǎng)分的含量分布及形態(tài)特征是土壤肥力的重要標(biāo)志[1]。土壤養(yǎng)分本身存在著復(fù)雜的空間變異特征，在同一時(shí)間不同點(diǎn)存在著明顯的差異性，因此了解和掌握農(nóng)田土壤養(yǎng)分的空間變異規(guī)律是提高農(nóng)田土壤肥力的重要手段，也是科學(xué)施肥工作的重要依據(jù)[1]。隨著地理信息系統(tǒng)的廣泛應(yīng)用和土壤科學(xué)的發(fā)展，運(yùn)用地統(tǒng)計(jì)學(xué)對(duì)土壤養(yǎng)分的空間變異的影響層出不窮，而土壤理化性質(zhì)的時(shí)空變異特征也成為土壤學(xué)的研究熱點(diǎn)[2]。因此在保證穩(wěn)產(chǎn)、高產(chǎn)、優(yōu)質(zhì)的前提下，將GIS與地統(tǒng)計(jì)學(xué)相結(jié)合起來(lái)研究農(nóng)田土壤養(yǎng)分空間變異特性具有重要的意義。【前人研究進(jìn)展】國(guó)內(nèi)學(xué)者對(duì)作物的土壤養(yǎng)分的空間特征進(jìn)行了諸多的研究與分析，張彬等研究得出陜西省禮泉縣蘋果主產(chǎn)區(qū)的土壤有機(jī)質(zhì)、堿解氮、有效磷和速效鉀的養(yǎng)分空間分布特征[3]，叢日環(huán)等研究了長(zhǎng)江中游油菜種植區(qū)的土壤養(yǎng)分的等級(jí)水平和空間分布圖[4]，王永豪等研究了川北煙區(qū)的土壤有效磷的空間變異特征及其主控因素，為煙區(qū)的平衡施肥和種植區(qū)劃提供依據(jù)[2]，金明清等對(duì)四川省的植煙土壤養(yǎng)分指標(biāo)的空間變異規(guī)律和地質(zhì)影響因素進(jìn)行了研究[5]。但香芋作為我國(guó)南方重要的經(jīng)濟(jì)作物之一，目前主要是在香芋的栽培技術(shù)、加工工藝和品種等方面進(jìn)行研究，僅有楊飛等對(duì)韶關(guān)香芋產(chǎn)區(qū)的土壤條件進(jìn)行了研究[6]，而對(duì)香芋作物產(chǎn)區(qū)的土壤空間變異性特征鮮有報(bào)道。韶關(guān)市位于廣東省的粵北地區(qū)，而粵北地區(qū)具有優(yōu)越的生態(tài)環(huán)境條件，這為韶關(guān)的農(nóng)業(yè)發(fā)展提供了得天獨(dú)厚的條件。近年來(lái)，韶關(guān)在穩(wěn)定糧食生產(chǎn)的同時(shí)，大力發(fā)展特色經(jīng)濟(jì)作物。其中，香芋是韶關(guān)重要的經(jīng)濟(jì)作物之一，更成為出口創(chuàng)匯的支柱產(chǎn)業(yè)，在全國(guó)也有一定的知名度。隨著近年來(lái)化肥生產(chǎn)和消費(fèi)量持續(xù)增長(zhǎng)，香芋的產(chǎn)量有所提升，但是根據(jù)田間試驗(yàn)的統(tǒng)計(jì)結(jié)果，目前化肥利用效率低，這既浪費(fèi)資源，又增加生產(chǎn)成本，而且不可避免地導(dǎo)致土壤中某些區(qū)域肥料用量過(guò)高或過(guò)低，造成土壤養(yǎng)分供應(yīng)與作物養(yǎng)分需求在空間上不同步，同時(shí)大部分未被利用的養(yǎng)分會(huì)通過(guò)滲漏和徑流直接進(jìn)入環(huán)境造成污染[8]?！颈狙芯壳腥朦c(diǎn)】因此本文通過(guò)室外取樣和室內(nèi)試驗(yàn)，運(yùn)用傳統(tǒng)統(tǒng)計(jì)學(xué)SPSS軟件、地統(tǒng)計(jì)學(xué)GS+軟件的半方差函數(shù)和GIS中克里金插值的研究韶關(guān)香芋主產(chǎn)區(qū)土壤養(yǎng)分含量、空間變異特征和空間分布?！緮M解決的關(guān)鍵問題】以期為韶關(guān)市香芋種植區(qū)劃、科學(xué)施肥和產(chǎn)量品質(zhì)的提升提供依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)域位于韶關(guān)市乳源縣的桂頭鄉(xiāng)、游溪鎮(zhèn)、一六鄉(xiāng)、乳城鄉(xiāng)、大橋鎮(zhèn)、大布鄉(xiāng)，和樂昌市的樂城鎮(zhèn)、北鄉(xiāng)鎮(zhèn)、長(zhǎng)來(lái)鎮(zhèn)、廊田鎮(zhèn)，地處東經(jīng)113°03′18″～113°04′27″ ，北緯24°56′24″～24°57′29″，面積約為 2600 km2(圖1左)。研究區(qū)地勢(shì)由西北向東南傾斜，西北部、西部峰巒環(huán)峙，屬高山地帶，溶蝕高原地貌顯著，是韶關(guān)市主要石灰?guī)r地區(qū)之一；東北部屬丘陵地帶，河流兩岸地勢(shì)平緩，海拔區(qū)間54～1833 m(圖1右)。氣候?qū)僦衼啛釒Ъ撅L(fēng)氣候，區(qū)間氣候懸殊，全年平均氣溫16～20 ℃。降雨量東南部及西部山區(qū)偏多，全年平均有2000 mm以上；南部和北部降雨量偏少，全年降雨量平均1400～1500 mm。研究區(qū)母質(zhì)主要為含碳酸鹽河流沖積物，土地利用方式以旱田和水田為主，主要是香芋-水稻-蔥等輪作。

1.2 土樣采集與測(cè)定

2016年5月-6月在韶關(guān)市乳源縣的桂頭鄉(xiāng)、游溪鎮(zhèn)、一六鄉(xiāng)、乳城鄉(xiāng)、大橋鎮(zhèn)、大布鄉(xiāng)，和樂昌市的樂城鎮(zhèn)、北鄉(xiāng)鎮(zhèn)、長(zhǎng)來(lái)鎮(zhèn)、廊田鎮(zhèn)選取具有代表性的采樣點(diǎn)進(jìn)行調(diào)查取樣。結(jié)合實(shí)際情況，用GPS定位采集204份香芋土壤樣品，保證每個(gè)土壤樣品具體到村、組、戶、地塊。采樣間距為1～2 km，深度為耕作(0～30 cm)，采用“S”型布點(diǎn)法，用四分法留取500 g左右裝入貼好標(biāo)簽的塑料袋中用于室內(nèi)指標(biāo)的測(cè)定。在取樣的同時(shí)采用實(shí)地測(cè)量和農(nóng)戶走訪等多種形式，詳細(xì)調(diào)查每個(gè)樣點(diǎn)的坡度、坡向、土層厚度、土地利用方式、耕作制度、施肥等信息，填寫土壤樣本調(diào)查表[9]。

圖1 研究區(qū)樣點(diǎn)分布圖(左)及地形圖(右)Fig.1 Sample points distribution (left) and topographic map (right) of studied areas

樣品在實(shí)驗(yàn)室自然風(fēng)干后，研磨過(guò)1和0.25 mm的網(wǎng)篩，備用。根據(jù)《土壤農(nóng)業(yè)化學(xué)分析方法》對(duì)樣品進(jìn)行常規(guī)測(cè)定，在測(cè)試方法上，土壤有機(jī)質(zhì)采用重鉻酸鉀氧化-滴定法測(cè)定，土壤堿解氮采用堿解擴(kuò)散法測(cè)定，有效磷采用碳酸氫鈉浸提-鉬銻抗比色法測(cè)定，速效鉀采用乙酸銨-火焰光度計(jì)發(fā)進(jìn)行測(cè)定。

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

利用SPSS19.0對(duì)土壤有機(jī)質(zhì)、堿解氮、有效磷和速效鉀的養(yǎng)分含量進(jìn)行經(jīng)典統(tǒng)計(jì)學(xué)分析；地統(tǒng)計(jì)學(xué)分析采用GS+9.0軟件對(duì)有機(jī)質(zhì)、堿解氮、有效磷和速效鉀進(jìn)行半方差函數(shù)分析，獲取其最優(yōu)理論半方差模型和參數(shù)；利用ArcGIS10.2地統(tǒng)計(jì)學(xué)模塊的正態(tài)檢驗(yàn)、趨勢(shì)分析和克里金插值繪制土壤有機(jī)質(zhì)、堿解氮、有效磷和速效鉀的養(yǎng)分空間分布圖。

2 結(jié)果與分析

2.1 土壤養(yǎng)分含量的總體特征

從表1可以看出，韶關(guān)香芋土壤有機(jī)質(zhì)、堿解氮、有效磷和速效鉀含量范圍分別為0.89 %～6.65 %，27.79～257.25 mg/kg，5.79～108.44 mg/kg，15.03～270.18 mg/kg，變異系數(shù)為35.2 %，46 %，64.7 %，62.4 %，其中有效磷和速效鉀的含量變化范圍較大。有機(jī)質(zhì)和堿解氮含量均值為2.98 %，109.79 mg/kg，均處于中等水平，有效磷均值為40.30 mg/kg，處于豐富水平，速效鉀含量均值為81.60 mg/kg，處于缺乏水平。其中86.3 %的樣點(diǎn)有機(jī)質(zhì)處于中等以上等級(jí)，57.3 %的樣點(diǎn)堿解氮處于中等以上，92.1 %的樣點(diǎn)有效磷處于中等以上，但是76.7 %的樣點(diǎn)速效鉀處于中等以下。這與之前研究的樂昌張溪村土壤肥力情況相接近，張溪村土壤有機(jī)質(zhì)處于中等水平，速效鉀處于缺乏水平，堿解氮和有效磷處于豐富水平。

按不同鄉(xiāng)鎮(zhèn)區(qū)域(表2)來(lái)看，有機(jī)質(zhì)含量依次乳城鄉(xiāng) (3.5 %)> 桂頭鄉(xiāng)(3.12 %) > 大橋鎮(zhèn)(2.79 %) > 大布鄉(xiāng)(2.69 %) > 一六鄉(xiāng)(2.55 %) > 長(zhǎng)來(lái)鎮(zhèn)(2.15 %) >游溪鎮(zhèn)(2.07 %) > 廊田鎮(zhèn)(1.89 %) > 樂城鎮(zhèn) (1.79 %)> 北鄉(xiāng)鎮(zhèn)(1.59 %)，堿解氮含量依次游溪鎮(zhèn)(134.75 mg/kg) > 一六鄉(xiāng) (129.2 mg/kg)> 大橋鎮(zhèn)(121.89 mg/kg) > 長(zhǎng)來(lái)鎮(zhèn)(120.75 mg/kg) > 廊田鎮(zhèn)(113.75 mg/kg) > 桂頭鄉(xiāng)(111.41 mg/kg) > 樂城鎮(zhèn)(105 mg/kg) > 大布鄉(xiāng)(88.97 mg/kg) > 北鄉(xiāng)鎮(zhèn)(85.75 mg/kg) > 乳城鄉(xiāng)(79.19 mg/kg)，有效磷含量依次樂城鎮(zhèn)(102 mg/kg) >北鄉(xiāng)鎮(zhèn)(82.3 mg/kg) > 游溪鎮(zhèn)(81.92 mg/kg) > 大布鄉(xiāng)(80.72 mg/kg) > 長(zhǎng)來(lái)鎮(zhèn)(62.09 mg/kg) > 一六鄉(xiāng)(43.56 mg/kg) > 大橋鎮(zhèn)(31.73 mg/kg) > 桂頭鄉(xiāng)(29.39 mg/kg) > 乳城鄉(xiāng)(27.32 mg/kg) > 廊田鎮(zhèn)(5.79 mg/kg)，速效鉀的含量依次為游溪鎮(zhèn)(270.18 mg/kg) > 北鄉(xiāng)鎮(zhèn)(250.48 mg/kg) > 樂城鎮(zhèn)(148.52 mg/kg) > 大橋鎮(zhèn)(112.14 mg/kg) > 桂頭鄉(xiāng)(74.09 mg/kg) > 一六鄉(xiāng)(68.54 mg/kg) > 乳城鄉(xiāng)(67.4 mg/kg) > 長(zhǎng)來(lái)鎮(zhèn)(65.63 mg/kg) > 廊田鎮(zhèn)(63.81 mg/kg) > 大布鄉(xiāng)(45.1 mg/kg)。

表2 不同區(qū)域香芋土壤養(yǎng)分含量及級(jí)別Table 2 Soil nutrient content and level of fragrant taro in different areas

2.2 土壤養(yǎng)分含量的空間變異性分析

半方差函數(shù)分析是地統(tǒng)計(jì)學(xué)分析的基礎(chǔ)，運(yùn)用地統(tǒng)計(jì)學(xué)探索性分析對(duì)香芋土壤養(yǎng)分指標(biāo)進(jìn)行預(yù)處理，有機(jī)質(zhì)、堿解氮、有效磷和速效鉀數(shù)據(jù)經(jīng)過(guò)對(duì)數(shù)轉(zhuǎn)換后均符合正態(tài)分布，均存在全局趨勢(shì)采用一次剔除。半方差函數(shù)通過(guò)“變程”反映變量的影響范圍，變程A的大小反映了區(qū)域化變量影響范圍的大小，即區(qū)域化變量空間變異尺度或空間自相關(guān)尺度。從結(jié)構(gòu)性因素角度看，C0/(C0+C)的比例可表示系統(tǒng)變量空間相關(guān)程度，如果比值<25 %，則說(shuō)明系統(tǒng)具有強(qiáng)烈的空間相關(guān)性，如果比值在25 %～75 %之間表明系統(tǒng)具有中等的空間相關(guān)性；如果比值大于75 %說(shuō)明系統(tǒng)在空間相關(guān)中占據(jù)的比值很小，而隨機(jī)性在空間相關(guān)中比值較大。

從半方差函數(shù)擬合結(jié)果(表3、圖2)來(lái)看，香芋主產(chǎn)區(qū)的土壤養(yǎng)分的空間相關(guān)性依次為 速效鉀 >有機(jī)質(zhì) > 堿解氮 > 有效磷。其中土壤有機(jī)質(zhì)符合球狀模型，變程28.3 km，說(shuō)明有機(jī)質(zhì)的空間變異尺度為28.3 km，塊金系數(shù)47.8 %，比值在25 %～50 %，屬于中等空間相關(guān)性；堿解氮符合指數(shù)模型，變程22.47 km，說(shuō)明堿解氮的空間變異尺度為22.47 km，塊金系數(shù)是 37.6 %，比值在25 %～50 %，屬于中等空間相關(guān)性；有效磷符合高斯模型，變程是61.38 km，塊金系數(shù)是32.5 %，屬于中等空間相關(guān)性，空間變異尺度為61.38 km；速效鉀符合高斯模型，變程是56.29 km，塊金系數(shù)是49.98 %，屬于中等空間相關(guān)性，空間變異尺度為56.29 km。

表3 土壤養(yǎng)分指標(biāo)的變異函數(shù)理論模型及相關(guān)參數(shù)Table 3 Oretical isotropic semivariogram model and corresponding parameters of soil nutrient

圖2 香芋土壤有機(jī)質(zhì)、堿解氮、有效磷和速效鉀的半方差函數(shù)Fig.2 Isotropic semivariogram of organic matter, available N, P and K contents of fragrant taro

表4 土壤養(yǎng)分分級(jí)標(biāo)準(zhǔn)Table 4 Grading standard of soil nutrients

2.3 土壤養(yǎng)分含量的空間分布特征

根據(jù)香芋土壤有機(jī)質(zhì)、堿解氮、有效磷和速效鉀含量的空間變異性，利用ArcGIS采用普通克里金法插值形成土壤有機(jī)質(zhì)、堿解氮、有效磷和速效鉀的空間分布圖(圖3)。根據(jù)土壤養(yǎng)分分級(jí)標(biāo)準(zhǔn)可以看出(表4)，香芋主產(chǎn)區(qū)的有機(jī)質(zhì)含量主要在2 %～4 %的范圍內(nèi)，處于養(yǎng)分級(jí)別的中等和豐富水平，其中乳源縣的桂頭鄉(xiāng)、一六鄉(xiāng)、大布鄉(xiāng)和游溪鎮(zhèn)，樂昌的北鄉(xiāng)鎮(zhèn)、樂城鎮(zhèn)、廊田鎮(zhèn)和長(zhǎng)來(lái)鎮(zhèn)的有機(jī)質(zhì)含量處于中等水平，而乳源縣的大橋鎮(zhèn)、乳城鎮(zhèn)有機(jī)質(zhì)含量處于豐富水平。堿解氮的含量主要在60～150 mg/kg的范圍內(nèi)，其中樂昌的北鄉(xiāng)鎮(zhèn)、樂城鎮(zhèn)、廊田鎮(zhèn)和長(zhǎng)來(lái)鎮(zhèn)和乳源縣的桂頭鄉(xiāng)、一六鄉(xiāng)堿解氮含量處于中等水平，而乳源縣的大橋鎮(zhèn)、乳城鎮(zhèn)堿解氮含量處于豐富水平，而乳源縣的桂頭鄉(xiāng)、一六鄉(xiāng)和大布鄉(xiāng)的東南部堿解氮處于缺乏水平。有效磷的含量在20 mg/kg以上，其中乳源縣的桂頭鄉(xiāng)、一六鄉(xiāng)、大橋鎮(zhèn)、乳城鎮(zhèn)、游溪鎮(zhèn)，和樂昌的北鄉(xiāng)鎮(zhèn)、樂城鎮(zhèn)、廊田鎮(zhèn)和長(zhǎng)來(lái)鎮(zhèn)的有效磷含量處于豐富水平，而乳源縣的大布鄉(xiāng)的有效磷的含量處于很豐富水平。香芋主產(chǎn)區(qū)的速效鉀主要在50～100 mg/kg的范圍內(nèi)，主要處于缺乏的水平。

圖3 香芋土壤有機(jī)質(zhì)、堿解氮、有效磷和速效鉀含量空間分布Fig.3 Spatial distributions of organic matter, available N, P and K contents of fragrant taro

3 討 論

近些年隨著地理信息系統(tǒng)的快速發(fā)展，許多學(xué)者運(yùn)用GIS和地統(tǒng)計(jì)學(xué)對(duì)不同作物產(chǎn)區(qū)的土壤養(yǎng)分的空間變異特征進(jìn)行研究，揭示了陜西省禮泉縣蘋果主產(chǎn)區(qū)[3]、長(zhǎng)江中游油菜主產(chǎn)區(qū)[4]、川北和四川省植煙區(qū)[2,5]等不同背景下的土壤養(yǎng)分的空間變異特征。相關(guān)研究表明，不同作物產(chǎn)區(qū)土壤養(yǎng)分的空間變異性差異明顯，分布特征即受到土壤類型、地形地貌、海拔等結(jié)構(gòu)因素的影響，又會(huì)受到土地利用方式、耕作制度、施肥等隨機(jī)因素的影響。這與王巖、張超和羅新寧等人的研究一致，土壤養(yǎng)分元素的含量和形態(tài)分布不僅與氣候、母質(zhì)、成土過(guò)程、地形有關(guān)，還受到耕作制度、施用肥料和作物類型等人為因素的影響[10-12]。但是目前對(duì)香芋產(chǎn)區(qū)的土壤養(yǎng)分空間變異特征方面研究甚少，因此，立足于廣東韶關(guān)市典型香芋產(chǎn)區(qū)，在保證穩(wěn)產(chǎn)、高產(chǎn)、優(yōu)質(zhì)的前提下，將GIS與地統(tǒng)計(jì)學(xué)相結(jié)合起來(lái)研究農(nóng)田土壤養(yǎng)分空間變異特性具有重要的意義，也為進(jìn)一步探討韶關(guān)香芋產(chǎn)區(qū)水肥資源的利用效率、肥料的精準(zhǔn)施用和種植規(guī)劃打下基礎(chǔ)。

目前較多的研究證明，農(nóng)田土壤耕層有效磷含量水平大幅提高，這主要與磷肥的廣泛施用有關(guān)[2]，本研究根據(jù)農(nóng)戶走訪調(diào)查，農(nóng)民在施肥方面，主要無(wú)機(jī)肥和有機(jī)肥聯(lián)合施用，尤其重視磷肥施用，這可能是種植區(qū)土壤有效磷有機(jī)質(zhì)含量豐富的原因之一，在對(duì)廣東樂昌香芋產(chǎn)區(qū)的土壤養(yǎng)分特征與香芋品質(zhì)的相關(guān)性研究中，也得到土壤有效磷含量與香芋總糖含量正相關(guān)，有機(jī)質(zhì)和香芋Vc含量正相關(guān)，而速效鉀和香芋香芋Vc含量負(fù)相關(guān)[13]，因此農(nóng)戶在保證產(chǎn)量和品質(zhì)的前提下，不重視鉀肥的施入，這可能也是種植區(qū)速效鉀含量偏低的原因之一。 研究區(qū)地勢(shì)起伏度變化大，地勢(shì)西高東低，種植區(qū)主要分布在東部河谷地帶，母質(zhì)以沖積物為主，土壤質(zhì)地多為粘壤和粉壤為主，具有一定的保水保肥能力，相比其他作物，產(chǎn)區(qū)種植制度比較單一，主要是香芋-蔥-水稻輪作，總體來(lái)看肥料補(bǔ)給有限而且人為影響較小。土壤有機(jī)質(zhì)、堿解氮、有效磷和速效鉀的塊金系數(shù)主要在30 %～50 %之間，更多是受到地形和海拔等結(jié)構(gòu)性因素的影響，因此，在后續(xù)研究中，要進(jìn)一步研究香芋產(chǎn)區(qū)土壤養(yǎng)分空間變異特征與地形因素之前的關(guān)系。除此之外，通過(guò)對(duì)韶關(guān)香芋香芋主產(chǎn)區(qū)的土壤養(yǎng)分級(jí)別和空間變異特征研究，在此基礎(chǔ)上可進(jìn)一步研究空間和時(shí)間上的差異性信息，分析影響產(chǎn)量差異的原因，采取技術(shù)上可行、經(jīng)濟(jì)上有效的施肥調(diào)控措施[14]。

4 結(jié) 論

(1) 韶關(guān)香芋土壤有機(jī)質(zhì)和堿解氮含量均值為2.98 %，109.79 mg/kg，均處于中等水平，有效磷均值為40.30 mg/kg，處于豐富水平，速效鉀含量均值為81.60 mg/kg，處于缺乏水平。其中不同鄉(xiāng)鎮(zhèn)區(qū)域有機(jī)質(zhì)含量依次乳城鄉(xiāng) > 桂頭鄉(xiāng) > 大橋鎮(zhèn) > 大布鄉(xiāng) > 一六鄉(xiāng) > 長(zhǎng)來(lái)鎮(zhèn) >游溪鎮(zhèn) > 廊田鎮(zhèn) > 樂城鎮(zhèn) > 北鄉(xiāng)鎮(zhèn)，堿解氮含量依次游溪鎮(zhèn)> 一六鄉(xiāng) > 大橋鎮(zhèn) > 長(zhǎng)來(lái)鎮(zhèn)> 廊田鎮(zhèn) > 桂頭鄉(xiāng) > 樂城鎮(zhèn) > 大布鄉(xiāng) > 北鄉(xiāng)鎮(zhèn)> 乳城鄉(xiāng)，有效磷含量依次樂城鎮(zhèn) >北鄉(xiāng)鎮(zhèn) > 游溪鎮(zhèn) > 大布鄉(xiāng) > 長(zhǎng)來(lái)鎮(zhèn) > 一六鄉(xiāng) > 大橋鎮(zhèn)> 桂頭鄉(xiāng) > 乳城鄉(xiāng) > 廊田鎮(zhèn)，速效鉀的含量依次為游溪鎮(zhèn) > 北鄉(xiāng)鎮(zhèn)> 樂城鎮(zhèn) > 大橋鎮(zhèn) > 桂頭鄉(xiāng)> 一六鄉(xiāng) > 乳城鄉(xiāng) > 長(zhǎng)來(lái)鎮(zhèn)> 廊田鎮(zhèn) > 大布鄉(xiāng)。

(2) 香芋主產(chǎn)區(qū)的土壤養(yǎng)分的空間相關(guān)性程度依次為速效鉀> 有機(jī)質(zhì) > 堿解氮 > 有效磷 ，均屬于中等空間相關(guān)性，結(jié)構(gòu)性變異為主，隨機(jī)性變異次之。其中土壤有機(jī)質(zhì)符合球狀模型，在28.3 km內(nèi)存在空間相關(guān)性；堿解氮符合指數(shù)模型，在22.47 km內(nèi)存在空間相關(guān)性；有效磷符合高斯模型，在61.38 m內(nèi)存在空間相關(guān)性；速效鉀符合高斯模型，變56.29 km內(nèi)存在空間相關(guān)性。

(3) 香芋主產(chǎn)區(qū)的土壤養(yǎng)分總體上呈現(xiàn)北低南高的趨勢(shì)，其中有機(jī)質(zhì)含量主要在2 %～4 %的范圍內(nèi)，處于養(yǎng)分級(jí)別的中等和豐富水平，堿解氮的含量主要在60～150 mg/kg的范圍內(nèi)，其中樂昌的北鄉(xiāng)鎮(zhèn)、樂城鎮(zhèn)、廊田鎮(zhèn)和長(zhǎng)來(lái)鎮(zhèn)和乳源縣的桂頭鄉(xiāng)、一六鄉(xiāng)堿解氮含量處于中等水平，而乳源縣的大橋鎮(zhèn)、乳城鎮(zhèn)堿解氮含量處于豐富水平，而乳源縣的桂頭鄉(xiāng)、一六鄉(xiāng)和大布鄉(xiāng)的東南部堿解氮處于缺乏水平。有效磷的含量在20 mg/kg以上，大部分地區(qū)處于豐富水平，而乳源縣的大布鄉(xiāng)的有效磷的含量處于很豐富水平。香芋主產(chǎn)區(qū)的速效鉀主要在50～100 mg/kg的范圍內(nèi)，大部分地區(qū)處于缺乏的水平。

參考文獻(xiàn)：

[1]王雪梅，柴仲平，武紅旗.典型干旱荒漠綠洲區(qū)耕層土壤養(yǎng)分空間變異[J].水土保持通報(bào)，2016，36(1)： 51-56.

[2]王永豪，王昌全，李啟全.川北煙區(qū)土壤有效磷空間變異特征及主控因素[J].核農(nóng)學(xué)報(bào)，2016，30(12): 2425-2433.

[3]張 彬，楊聯(lián)安，粉 莉，等.蘋果主產(chǎn)區(qū)土壤養(yǎng)分空間分布特征及其影響因素—以陜西省禮泉縣為例[J].土壤，2016，48(4)：777-784.

[4]叢日環(huán)，張 智，鄭 磊，等. 基于GIS的長(zhǎng)江中游油菜種植區(qū)土壤養(yǎng)分及pH狀況[J].土壤學(xué)報(bào)，2016，53(5)：1213-1224.

[5]金明清，彭月月，王 佩，等. 四川省鹽源縣植煙土壤氮磷鉀空間變異特征及影響因素[J]. 土壤，2016，48(5)：984-991.

[6]楊 飛. 土壤條件對(duì)香芋品質(zhì)影響的研究[D]. 內(nèi)蒙古：內(nèi)蒙古農(nóng)業(yè)大學(xué)，2011：5-7.

[7]黃立梅，黃紹文，劉雙全，等. 規(guī)模經(jīng)營(yíng)稻田土壤養(yǎng)分空間變異特征與高效平衡施肥效應(yīng)[J]. 華北農(nóng)學(xué)報(bào)，2010(2)：205-211.

[8]農(nóng)業(yè)部.中華人民共和國(guó)農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn):測(cè)土配方施肥技術(shù)規(guī)范[S]. 北京:中國(guó)農(nóng)業(yè)出版社，2006：27-59.

[9]王 巖，陳永金，劉加珍. 黃河三角洲濕地土壤養(yǎng)分空間分布特征[J]. 人民黃河，2013，35(2)：72-74.

[10]張 超，劉國(guó)彬. 黃土丘陵區(qū)撂荒農(nóng)耕地土壤有效態(tài)微量元素演變特征[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2013，46(18)：3809-3817.

[11]羅新寧，朱友娟，張宏勇，等. 塔里木綠洲種植制度對(duì)棉田土壤養(yǎng)分性狀的影響[J]. 干旱地區(qū)農(nóng)業(yè)研究，2012，30(3)：114-118.

[12]高 琳. 廣東樂昌市香芋產(chǎn)區(qū)土壤養(yǎng)分特征與香芋品質(zhì)的相關(guān)性[J]. 貴州農(nóng)業(yè)科學(xué)，2014，42(8)：179-181.

[13]高祥照，胡克林，郭 焱，等. 土壤養(yǎng)分與作物產(chǎn)量的空間變異特征與精確施肥[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2002，35(6)：660-666.