金明清，彭月月，王 佩，馮 浪，，陳玉藍(lán)，杜宣延，李啟權(quán)，王昌全*

（1 四川農(nóng)業(yè)大學(xué)資源學(xué)院，成都 611130；2 四川中煙工業(yè)有限責(zé)任公司，成都 610000；3 四川省煙草公司涼山州公司，四川涼山 615000）

四川省鹽源縣植煙土壤氮磷鉀空間變異特征及影響因素①

金明清1，彭月月1，王 佩2，馮 浪1，3，陳玉藍(lán)3，杜宣延1，李啟權(quán)1，王昌全1*

（1 四川農(nóng)業(yè)大學(xué)資源學(xué)院，成都 611130；2 四川中煙工業(yè)有限責(zé)任公司，成都 610000；3 四川省煙草公司涼山州公司，四川涼山 615000）

利用地統(tǒng)計(jì)學(xué)和回歸分析方法分析四川省鹽源縣植煙土壤堿解氮、有效磷、速效鉀含量空間變異特征以及不同因素對(duì)其空間變異的影響，以期為鹽源烤煙種植區(qū)劃、科學(xué)施肥和產(chǎn)質(zhì)量提升提供科學(xué)依據(jù)。結(jié)果表明：①堿解氮、有效磷和速效鉀的均值分別為102.84、10.18、148.10 mg/kg，整體處于中等水平，變異系數(shù)分別為38.66%、90.08%、60.65%，達(dá)到中等變異程度，空間分布不均；②堿解氮、有效磷空間分布符合指數(shù)模型，具有中等空間自相關(guān)，速效鉀符合高斯模型，具有明顯空間自相關(guān)，空間變異受結(jié)構(gòu)性因素和隨機(jī)性因素共同影響；③鹽源植煙土壤堿解氮、有效磷和速效鉀空間變異的主控因素有所不同，堿解氮為坡度和前茬作物，有效磷為土壤類型和前茬作物，速效鉀為海拔和坡度，因此植煙區(qū)規(guī)劃和施肥需要綜合考慮不同的影響因素。

植煙土壤；堿解氮；有效磷；速效鉀；空間變異；影響因素；四川鹽源

土壤養(yǎng)分含量的高低及相互的協(xié)調(diào)性對(duì)烤煙的產(chǎn)量和質(zhì)量具有十分重要的影響［1］。氮素是影響烤煙產(chǎn)量和品質(zhì)重要元素之一［2-4］；磷素對(duì)烤煙生長發(fā)育、生理生化和產(chǎn)質(zhì)量同樣發(fā)揮著不可替代的作用［5-7］；鉀素是烤煙在生長發(fā)育過程中吸收和積累量最大的元素，供應(yīng)充足是獲得優(yōu)質(zhì)煙葉的重要條件［8-10］。土壤速效養(yǎng)分是不經(jīng)轉(zhuǎn)化能迅速被作物直接吸收利用的土壤營養(yǎng)元素［11］，對(duì)烤煙的影響更加直接和深刻。近年來，許多學(xué)者利用地統(tǒng)計(jì)學(xué)方法對(duì)不同類型植煙區(qū)的土壤養(yǎng)分空間變異做了大量的研究，如王寶順等［12］找出影響重慶彭水坡地植煙土壤速效養(yǎng)分空間變異的主要因素，并掌握坡地植煙土壤養(yǎng)分的空間分布特征及其變異規(guī)律；劉遜等［13］采用地統(tǒng)計(jì)學(xué)方法研究并掌握了湖南湘西植煙土壤有機(jī)質(zhì)的空間分布及影響因素；王育軍等［14］采用地統(tǒng)計(jì)方法建立云南羅平縣植煙土壤肥力適宜性指數(shù)的時(shí)空變異特征。

四川省鹽源縣是我省烤煙種植大縣，盛產(chǎn)以“清、甜、香”為特色的優(yōu)質(zhì)煙葉，在四川優(yōu)質(zhì)煙葉生產(chǎn)中占有重要地位。本文利用地統(tǒng)計(jì)學(xué)和回歸分析方法對(duì)鹽源植煙土壤堿解氮、有效磷、速效鉀含量進(jìn)行分析，重點(diǎn)研究其空間分布特征和影響因素，以期為鹽源烤煙種植區(qū)劃、科學(xué)施肥和產(chǎn)質(zhì)量提升提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

鹽源縣位于四川省西南部，青藏高原東南緣，雅礱江下游西岸，位于100°42′09″ ～ 102°03′44″E、27°06′31 ～ 28°16′31″N，面積為8 398.6 km2，海拔多在2 300 ～ 2 800 m之間，地形復(fù)雜，地貌類型多樣，是一個(gè)山高坡陡、河谷深切、山巒疊障、盆地居中為總特征的典型山區(qū)縣。受高空西風(fēng)氣流的南支急流和印度洋氣流的影響，全縣屬于典型大陸性氣候，四季分明，年溫差小，日溫差大，年均氣溫13.2℃，冬春干旱，夏秋雨量集中，年均降水量727.30 mm，雨熱同期，光照資源非常豐富；縣內(nèi)土壤類型主要有紅壤、水稻土、新積土、黃棕壤及紫色土等。近年來烤煙種植面積約為4 500 hm2，煙葉收購量約為 11 200 t 。

1.2 土樣采集與分析

根據(jù)鹽源縣主要土壤類型和煙草種植情況，兼顧代表性與均勻性的基礎(chǔ)上，遵循隨機(jī)、等量、同層次、多點(diǎn)混合的采樣原則，利用GPS 定位技術(shù)獲取每個(gè)采樣點(diǎn)的坐標(biāo)，采用網(wǎng)格法和多點(diǎn)混合法采集0 ～20 cm耕層的土壤樣品。樣品主要于2011年植煙前采集，采集土壤樣品317個(gè)，采集土壤樣品時(shí)詳細(xì)記錄其地形地貌、土壤類型、土地利用方式以及施肥狀況等環(huán)境信息，采樣點(diǎn)分布見圖1。土樣堿解氮、有效磷和速效鉀分別采用堿解擴(kuò)散法、鹽酸-氟化銨-鉬銻抗比色法和乙酸銨浸提-火焰光度計(jì)法測定［15］。

圖1 鹽源植煙土壤采樣點(diǎn)分布圖Fig. 1 Sampling sites of tobacco planting soil in Yanyuan

1.3 研究方法與數(shù)據(jù)處理

由于土壤本身的空間變異性，在進(jìn)行野外隨機(jī)區(qū)組設(shè)計(jì)或者隨機(jī)確定采樣點(diǎn)的過程中，可能使獲得的試驗(yàn)數(shù)據(jù)根本不具有統(tǒng)計(jì)學(xué)意義，而應(yīng)用地統(tǒng)計(jì)學(xué)方法可以有效解決這個(gè)問題［16］。本研究采用地統(tǒng)計(jì)方法對(duì)鹽源植煙土壤堿解氮、有效磷、速效鉀的空間變異特征進(jìn)行分析，擬合出堿解氮、有效磷、速效鉀各向同性時(shí)較優(yōu)的半方差函數(shù)模型。半方差函數(shù)模型主要根據(jù)擬合結(jié)果的決定系數(shù)和殘差的大小來決定?？臻g變異性強(qiáng)弱可根據(jù)塊金值與基臺(tái)值的比值，即塊金系數(shù)進(jìn)行劃分。塊金系數(shù)表示由隨機(jī)部分引起的空間變異性占總體變異的比例，若塊金系數(shù)＜25%，說明變量有很強(qiáng)的空間自相關(guān)性，25% ～ 50% 說明變量有明顯的空間自相關(guān)，50% ～ 75% 時(shí)變量有中等空間自相關(guān)，＞75% 時(shí)變量空間自相關(guān)弱，變異主要由隨機(jī)因素影響。變程是指變異函數(shù)達(dá)到基臺(tái)值所對(duì)應(yīng)的距離，它反映的是屬性因子空間自相關(guān)范圍的大?。?7-19］。

地統(tǒng)計(jì)分析首先采用GS+ 軟件進(jìn)行半方差分析，以獲取堿解氮、有效磷、速效鉀各向同性的最優(yōu)理論半方差模型和參數(shù)，再在ArcGIS軟件的地統(tǒng)計(jì)模塊中進(jìn)行Kring插值，得到研究區(qū)土壤堿解氮、有效磷、速效鉀的空間分布圖。回歸分析在SPSS軟件中完成，運(yùn)用回歸分析定量揭示各影響因素對(duì)植煙土壤堿解氮、有效磷、速效鉀空間變異的影響程度，即利用回歸方程的校正決定系數(shù)來判定各影響因素對(duì)土壤堿解氮、有效磷、速效鉀含量空間變異程度的獨(dú)立解釋能力［20］，并以此確定研究區(qū)堿解氮、有效磷、速效鉀空間分布的主控因素?？紤]到土壤類型和前茬作物為定性變量，海拔和坡度分級(jí)后也屬于定性變量，故先采用啞變量對(duì)土壤類型、海拔、坡度和前茬作物進(jìn)行賦值，再以啞變量作為自變量進(jìn)行回歸分析［21］。

1.4 植煙土壤堿解氮、有效磷、速效鉀適宜程度劃分標(biāo)準(zhǔn)

參考我國主要煙區(qū)土壤養(yǎng)分豐缺狀況評(píng)價(jià)結(jié)果［22］，綜合制定植煙土壤堿解氮、有效磷、速效鉀含量適宜程度分級(jí)標(biāo)準(zhǔn)（表1）。

表1 植煙土壤堿解氮、有效磷、速效鉀含量適宜程度分級(jí)標(biāo)準(zhǔn)（mg/kg）Table1 Suitable classification standard of available N， P and K contents of tobacco planting soil

2 結(jié)果分析與討論

2.1 植煙土壤堿解氮、有效磷、速效鉀總體特征

由表2可知鹽源植煙土壤堿解氮、有效磷和速效鉀含量范圍分別在14.3 ～ 209.4、0.4 ～ 51.0和2.0 ～409.0 mg/kg之間，含量均值分別為102.84、10.18和148.10 mg/kg，均值都處于適宜水平。從變異系數(shù)來看，堿解氮、有效磷和速效鉀的變異系數(shù)分別為38.66%、90.08% 和60.65%，均達(dá)到中等變異程度，其空間分布不均。40.38% 的樣點(diǎn)堿解氮處于適宜水平，30.60% 和16.09% 的樣點(diǎn)分別處于低和極低水平；23.03% 的樣點(diǎn)有效磷處于適宜水平，超過60%的樣點(diǎn)有效磷缺乏；僅17.03% 的樣點(diǎn)速效鉀處于適宜水平，超過60% 的樣點(diǎn)速效鉀缺乏?？傮w來看，鹽源植煙土壤堿解氮、有效磷和速效鉀含量處于中等偏下水平，部分區(qū)域需重視氮肥、磷肥和鉀肥的補(bǔ)充。

表2 鹽源植煙土壤堿解氮、有效磷、速效鉀含量描述性統(tǒng)計(jì)（mg/kg）Table2 Descriptive statistics of available N， P and K contents of tobacco planting soil in Yanyuan

2.2 植煙土壤堿解氮、有效磷、速效鉀空間結(jié)構(gòu)分析

從各向同性半方差函數(shù)擬合結(jié)果（表3、圖2）來看，鹽源植煙土壤堿解氮、有效磷含量分布符合指數(shù)模型，速效鉀符合高斯模型，擬合決定系數(shù)差異較大；速效鉀決定系數(shù)最大、堿解氮次之，有效磷最小。堿解氮、有效磷和速效鉀的塊金系數(shù)分別為56.31%、55.08% 和46.55%，說明堿解氮和有效磷具有中等空間自相關(guān)，速效鉀具有明顯空間自相關(guān)，空間變異受結(jié)構(gòu)性因素和隨機(jī)性因素共同影響。

表3 鹽源土壤堿解氮、有效磷、速效鉀含量各向同性半方差函數(shù)理論模型及相關(guān)參數(shù)Table3 Theoretical isotropic semivariogram models and corresponding parameters of soil available N， P and K contents in Yanyuan

圖2 鹽源植煙土壤堿解氮（A）、有效磷（B）、速效鉀（C）含量半方差函數(shù)圖Fig. 2 Semivariance models for available N（A）， P（B） and K（C） contents of tobacco planting soil in Yanyuan

2.3 植煙土壤堿解氮、有效磷、速效鉀空間分布特征分析

根據(jù)鹽源植煙土壤堿解氮、有效磷、速效鉀含量的空間變異特征，利用ArcGIS采用普通克里金法生成堿解氮、有效磷和速效鉀的空間分布圖（圖3）。從圖中可以看出，堿解氮、有效磷、速效鉀含量空間分布具有明顯差異性。堿解氮含量處于適宜水平 100 ～150 mg/kg之間的面積最大，主要分布在東西兩側(cè)；處于較低水平65 ～ 100 mg/kg之間的面積次之，主要分布在中部及南部；處于較高水平150 ～ 200 mg/kg之間的面積較小，主要分布西部、北部及東部的部分區(qū)域。有效磷含量處于低水平（5 ～ 10 mg/kg）范圍基本覆蓋全境，僅在西北、東北和東南有小部分區(qū)域有效磷含量達(dá)到適宜水平10 ～ 20 mg/kg。速效鉀含量處于低和適宜水平（80 ～ 150 mg/kg和150 ～ 220 mg/kg）范圍內(nèi)的面積最大，主要分布在中部、西部和南部；處于較高水平（220 ～ 350 mg/kg）范圍內(nèi)的面積較小，主要分布在北部和東北部。各養(yǎng)分含量的空間變異明顯，生產(chǎn)上可根據(jù)區(qū)域土壤養(yǎng)分實(shí)際含量進(jìn)行具有針對(duì)性的施肥，避免進(jìn)一步加劇土壤養(yǎng)分失衡。

2.4 植煙土壤堿解氮、有效磷、速效鉀影響因素分析

2.4.1 土壤類型 不同的土壤類型由不同的母質(zhì)類型、成土過程、發(fā)育程度和耕作管理措施等多種因素綜合作用而形成，不同土壤類型間養(yǎng)分元素的組成和分布具有明顯的差異性。由表4可知，鹽源植煙土壤主要以紅壤、水稻土和新積土為主。各土壤類型堿解氮含量均值為98.14 ～ 111.87 mg/kg，其中以水稻土最高、紅壤最低；各土壤類型的變異程度均為中等變異；紅壤、水稻土、新積土均有超過40% 的比例達(dá)到適宜程度，紫色土則需要進(jìn)一步補(bǔ)充氮肥。有效磷含量均值為8.43 ～ 13.32 mg/kg，以紫色土最高、紅壤最低；除黃棕壤達(dá)到強(qiáng)變異程度外，其余類型土壤達(dá)到中等變異程度；各土壤類型有效磷含量適宜度等級(jí)主要以極低和低水平為主。紫色土速效鉀含量均值以189.35 mg/kg達(dá)到最高，黃棕壤以173.79 mg/kg次之，紅壤速效鉀含量均值最低，為128.27 mg/kg；各類型土壤的變異系數(shù)為56.90% ～ 62.60%，達(dá)到中等變異程度；從適宜度等級(jí)來看，黃棕壤、紫色土達(dá)到適宜及其以上等級(jí)的比例較紅壤、水稻土和新積土高?？傮w來看，紅壤堿解氮、有效磷和速效鉀含量均較低，肥力低；水稻土堿解氮含量豐富、有效磷和速效鉀含量較高，較肥沃；紫色土有效磷、速效鉀含量豐富，土壤類型的差異可為烤煙施肥提供相應(yīng)參考。

圖3 鹽源植煙土壤堿解氮、有效磷、速效鉀含量空間分布圖Fig. 3 Spatial distribution of available N， P and K contents of tobacco planting soil in Yanyuan

表4 鹽源植煙土壤不同土壤類型堿解氮、有效磷、速效鉀含量的統(tǒng)計(jì)特征（mg/kg）Table4 Statistical characteristics of available N， P and K contents of tobacco planting soils with different types in Yanyuan

2.4.2 海拔 不同海拔的土壤環(huán)境具有不同的溫度、濕度和熱量，影響物質(zhì)的積累與分解程度，進(jìn)而影響土壤營養(yǎng)元素的含量。從表5可知，鹽源植煙區(qū)海拔主要在2 000 ～ 2 500 m范圍內(nèi)。海拔在1 500 ～2 000 m范圍內(nèi)堿解氮均值最大，為114.33 mg/kg，海拔在2 000 ～ 2 500 m范圍內(nèi)堿解氮均值最小，其值為98.38 mg/kg；各海拔變異系數(shù)在37.00% ～46.64% 之間，達(dá)到中等變異程度。有效磷含量均值以海拔2 000 m為分界線，20 00 m以下的有效磷明顯比2 000 m以上的高；從變異系數(shù)來看，除海拔＜1 500 m的變異系數(shù)達(dá)到強(qiáng)變異程度，其余的均達(dá)到中等變異強(qiáng)度。海拔＜1 500 m的速效鉀含量均值最小，為106.57 mg/kg，海拔＞2 500 m速效鉀含量均值最大，為174.73 mg/kg；各海拔的變異系數(shù)在51.63% ～82.42% 之間，均達(dá)到中等程度變異。總體來看，海拔在2 000 ～ 2 500 m范圍內(nèi)，植煙土壤堿解氮、有效磷和速效鉀含量均偏低，主要原因可能是研究區(qū)種植業(yè)主要集中在這一范圍內(nèi)，且耕作強(qiáng)度大、有機(jī)肥投入較低、無機(jī)肥施用不合理以及養(yǎng)分流失，進(jìn)而導(dǎo)致了該海拔土壤堿解氮、有效磷、速效鉀含量下降，土壤肥力降低。

表5 鹽源植煙土壤不同海拔土壤堿解氮、有效磷、速效鉀含量的統(tǒng)計(jì)特征（mg/kg）Table5 Statistical characteristics of available N， P and K contents of tobacco planting soil with different altitude in Yanyuan

表6 鹽源植煙土壤不同坡度土壤堿解氮、有效磷、速效鉀含量的統(tǒng)計(jì)特征（mg/kg）Table6 Statistical characteristics of available N， P and K contents of tobacco planting soil with different slopes in Yanyuan

2.4.3 坡度 從表6可知，鹽源植煙區(qū)坡度主要集中0 ～ 15° 范圍內(nèi)。各坡度堿解氮含量均值在90.58 ～112.62 mg/kg之間，堿解氮含量與坡度呈反比，坡度越低，堿解氮含量越大；各坡度堿解氮含量變異系數(shù)在34.43% ～ 47.18%，達(dá)到中等變異程度。有效磷含量隨坡度增高而減小，與坡度呈反比，坡度從低到高有效磷含量均值分別為12.36、10.53、9.36、8.91、 8.22 mg/kg；有效磷含量變異系數(shù)在83.87% ～115.94% 之間，除坡度 ＞25°達(dá)到強(qiáng)變異程度外，其余均為中等變異程度。速效鉀含量同樣與坡度呈反比，＜2° 均值高達(dá)170.64 mg/kg，而 ＞25° 均值僅為110.39 mg/kg；速效鉀含量變異系數(shù)在52.46% ～80.29% 之間，均達(dá)到中等變異程度。堿解氮、有效磷和速效鉀含量均與坡度呈反比，主要因?yàn)橥寥罓I養(yǎng)物質(zhì)主要以水為載體進(jìn)行遷移，坡度會(huì)促使土壤水分發(fā)生以下行為主的運(yùn)動(dòng)，進(jìn)而導(dǎo)致土壤養(yǎng)分由陡坡向緩坡富集。由于陡坡本身土壤養(yǎng)分缺乏，投入產(chǎn)出比低，耕作難度大，容易發(fā)生水土流失，所以在植煙區(qū)規(guī)劃時(shí)盡量減少陡坡植煙面積。

2.4.4 前茬作物 輪作方式不同，收獲后土壤中殘留的有效養(yǎng)分有差別，不同的前茬作物直接影響下茬作物生長的土壤養(yǎng)分［23］。由于鹽源煙田冬季以休耕為主，所以前茬作物一般指上一年度所種植的作物。從表7可知，在采樣點(diǎn)中鹽源烤煙前茬作物以玉米、烤煙、水稻、馬鈴薯和辣椒為主。堿解氮含量均值大小排序?yàn)椋核荆酒渌魑铮抉R鈴薯＞玉米＞烤煙＞辣椒，有效磷含量均值大小排序?yàn)椋嚎緹煟居衩祝抉R鈴薯＞辣椒＞其他作物＞水稻，速效鉀含量均值大小排序?yàn)椋浩渌魑铮抉R鈴薯＞玉米＞烤煙＞水稻＞辣椒，堿解氮、有效磷和速效鉀的適宜程度等級(jí)分布與各自含量的均值呈現(xiàn)出相似趨勢。烤煙、辣椒、馬鈴薯、水稻、玉米等作物在生長發(fā)育過程中對(duì)養(yǎng)分的吸收具有差異性，且農(nóng)民種植不同的作物施肥量具有明顯差異，例如這幾種作物在采樣點(diǎn)的氮肥（按N純量計(jì)算）施用量均值大小關(guān)系為：水稻（198.3 kg/hm2）＞馬鈴薯（184.5 kg/hm2）＞玉米（172.2 kg/hm2）＞辣椒（115.7 kg/hm2）＞烤煙（105.4 kg/hm2），磷肥（按P2O5純量計(jì)算）施用量均值大小關(guān)系為：烤煙（75.6 kg/hm2）＞玉米（71.0 kg/hm2）＞馬鈴薯（66.7 kg/hm2）＞辣椒（54.8 kg/hm2）＞水稻（51.1 kg/hm2），鉀肥（按K2O純量計(jì)算）施用量均值大小關(guān)系為：烤煙（397.2 kg/hm2）＞馬鈴薯（313.8 kg/hm2）＞辣椒（243.5 kg/hm2）＞玉米（157.9 kg/hm2）＞水稻（110.6 kg/hm2），因此長期單一耕作就會(huì)導(dǎo)致土壤某些營養(yǎng)元素的富集或虧缺，進(jìn)而出現(xiàn)土壤堿解氮、有效磷、速效鉀差異明顯的現(xiàn)象。所以合理的輪作制度對(duì)維持植煙土壤養(yǎng)分平衡、穩(wěn)定植煙土壤肥力、進(jìn)而保證和提高烤煙產(chǎn)質(zhì)量具有十分重要的意義。

表7 鹽源植煙土壤不同前茬作物土壤堿解氮、有效磷、速效鉀含量的統(tǒng)計(jì)特征（mg/kg）Table7 Statistical characteristics of available N， P and K contents of tobacco planting soil under different previous crops in Yanyuan

2.4.5 不同因素影響程度分析 分別對(duì)土壤類型、海拔、坡度和前茬作物進(jìn)行回歸分析，以定量揭示不同因素對(duì)鹽源植煙土壤堿解氮、有效磷、速效鉀空間變異的影響程度。由表8可知，土壤類型、海拔、坡度和前茬作物對(duì)堿解氮的獨(dú)立解釋能力分別為5.3%、4.7%、12.4%、10.7%，坡度和前茬作物的獨(dú)立解釋能力高于土壤類型和海拔的獨(dú)立解釋能力，堿解氮的空間變異主要受坡度和前茬作物的影響；各因素對(duì)有效磷的獨(dú)立解釋能力分別為13.2%、8.0%、9.1%、12.2%，土壤類型和前茬作物的獨(dú)立解釋能力高于海拔和坡度，有效磷的空間變異主要受土壤類型和前茬作物的影響；土壤類型、海拔、坡度和前茬作物對(duì)速效鉀的獨(dú)立解釋能力分別為12.7%、17.4%、16.6% 和12.0%，各因素的獨(dú)立解釋能力均較大，速效鉀的空間變異受海拔和坡度的影響更加明顯。綜上可知，結(jié)構(gòu)性因素和隨機(jī)性因素共同影響鹽源植煙土壤堿解氮、有效磷、速效鉀的空間變異，其中速效鉀的空間變異受結(jié)構(gòu)性因素影響更明顯。因此，在研究區(qū)進(jìn)行植煙區(qū)劃和科學(xué)施肥應(yīng)綜合考慮不同的影響因素。

表8 不同影響因素對(duì)鹽源植煙土壤堿解氮、有效磷、速效鉀含量的回歸分析結(jié)果Table8 Regression analysis of available N， P and K contents of tobacco planting soil with different influencing factors in Yanyuan

3 結(jié)論

1） 鹽源植煙土壤堿解氮、有效磷和速效鉀含量均值皆處于適宜水平。分別有46.69%、64.04%和60.26% 的樣點(diǎn)處于缺乏水平，且均達(dá)到中等變異程度，空間分布不均。

2） 鹽源植煙土壤堿解氮、有效磷含量分布符合指數(shù)模型，具有中等空間自相關(guān)；速效鉀符合高斯模型，具有明顯空間自相關(guān)，自相關(guān)范圍差異較大，空間變異受結(jié)構(gòu)性因素和隨機(jī)性因素共同影響。大部分煙田堿解氮含量處于適宜水平，主要分布在東西兩側(cè)；大部分煙田有效磷含量處于低和較低水平，僅在西部和東南的小部分區(qū)域達(dá)到適宜水平；大部分煙田速效鉀含量處于低和較低水平，主要分布在鹽源的西部和南部。

3） 鹽源植煙土壤堿解氮的主控因素為坡度和前茬作物、有效磷的主控因素為土壤類型和前茬作物，速效鉀的主控因素為海拔和坡度，因此在研究區(qū)進(jìn)行植煙區(qū)規(guī)劃和烤煙科學(xué)施肥需要考慮不同的影響因素。

［1］ 陳江華， 劉建利， 李志宏. 中國植煙土壤及煙草養(yǎng)分綜合管理［M］. 北京： 科學(xué)出版社， 2008

［2］ 李志， 史宏志， 劉國順， 等. 施氮量對(duì)皖南砂壤土烤煙碳氮代謝動(dòng)態(tài)變化的影響［J］. 土壤， 2010， 42（1）： 8-13

［3］ 張生杰， 黃元炯， 任慶成， 等. 氮素對(duì)不同品種烤煙葉片衰老、光合特性及產(chǎn)量和品質(zhì)的影響［J］. 應(yīng)用生態(tài)學(xué)報(bào)， 2010， 21（3）： 668-674

［4］ 劉青麗， 陳阜， 張?jiān)瀑F， 等. 我國西南煙區(qū)典型植煙土壤烤煙氮素的吸收規(guī)律［J］. 作物學(xué)報(bào)， 2013， 39（3）：486-493

［5］ 金霞， 趙正雄， 呂芬， 等. 施磷量對(duì)烤煙幾種生理生化物質(zhì)含量、赤星病發(fā)生及煙葉產(chǎn)質(zhì)量的影響［J］. 中國煙草學(xué)報(bào)， 2010， 16（3）： 53-56

［6］ 賈志紅， 易建華， 符建國， 等. 磷肥處理對(duì)烤煙生長生理及根系構(gòu)型的影響［J］. 土壤， 2011， 43（3）： 388-391

［7］ 郭燕， 許自成， 畢慶文， 等. 恩施煙區(qū)土壤磷素與烤煙磷含量的關(guān)系［J］. 土壤通報(bào)， 2010， 41（2）： 403-407

［8］ 代曉燕， 郭春燕， 王海波， 等. 鉀肥施用方式對(duì)豫西烤煙鉀含量及產(chǎn)質(zhì)量的影響［J］. 中國煙草學(xué)報(bào)， 2013， 18（3）：42-49

［9］ 鄔石根， 白厚義， 張超蘭， 等. 鉀肥和營養(yǎng)調(diào)節(jié)劑對(duì)烤煙含鉀量及重金屬含量的效應(yīng)研究［J］. 土壤通報(bào)， 2012，43（1）： 190-194

［10］ 鄭東方， 許嘉陽， 許自成， 等. 鉀肥和腐殖酸互作對(duì)烤煙有機(jī)鉀鹽指數(shù)的影響［J］. 土壤學(xué)報(bào)， 2015， 52（3）：169-177

［11］ 黃昌勇， 徐建明. 土壤學(xué)［M］. 北京： 中國農(nóng)業(yè)出版社，2010

［12］ 王寶峰， 劉國順， 葉協(xié)鋒. 坡地植煙土壤速效養(yǎng)分空間分析［J］. 西南農(nóng)業(yè)學(xué)報(bào)， 2012， 25（1）： 217-222

［13］ 劉遜， 鄧小華， 周米良， 等. 湘西植煙土壤有機(jī)質(zhì)含量分布及其影響因素［J］. 核農(nóng)學(xué)報(bào)， 2012， 26（7）： 1 037-1 042

［14］ 王育軍， 周冀衡， 孫書斌， 等. 云南省羅平縣煙區(qū)土壤肥力適宜性評(píng)價(jià)及養(yǎng)分時(shí)空變異特征［J］. 土壤， 2015，47（3）： 515-523

［15］ 張甘霖， 龔子同. 土壤調(diào)查實(shí)驗(yàn)室分析方法［M］. 北京：科學(xué)出版社， 2012

［16］ 劉愛利， 王培法， 丁園圓. 地統(tǒng)計(jì)學(xué)概論［M］. 北京： 科學(xué)出版社， 2012

［17］ 馬渝欣， 李徐生， 李德成， 等. 安徽省皖南山區(qū)農(nóng)田表層土壤有機(jī)碳空間變異及其影響因素——以宣州區(qū)為例［J］. 土壤通報(bào)， 2014， 45（6）： 1 424-1 429

［18］ 李啟權(quán)， 王昌全， 張文江， 等. 基于神經(jīng)網(wǎng)絡(luò)模型和地統(tǒng)計(jì)學(xué)方法的土壤養(yǎng)分空間分布預(yù)測［J］. 應(yīng)用生態(tài)學(xué)報(bào)，2013， 24（2）： 459-466

［19］ 臘貴曉， 顧懷勝， 劉國順， 等. 喀斯特地區(qū)煙田土壤養(yǎng)分的空間變異特征［J］. 水土保持研究， 2012， 19（3）： 48-53

［20］ 顧成軍， 史學(xué)正， 于東升， 等. 省域土壤有機(jī)碳空間分布的主控因子—土壤類型與土地利用比較［J］. 土壤學(xué)報(bào)，2013， 50（3）： 1-8

［21］ 羅由林， 李啟權(quán)， 王昌全， 等. 四川省仁壽縣土壤有機(jī)碳空間分布特征及其主控因素［J］. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)，2015， 23（1）： 34-42

［22］ 陳江華， 李志宏， 劉建利， 等. 全國主要煙區(qū)土壤養(yǎng)分豐缺狀況評(píng)價(jià)［J］. 中國煙草學(xué)報(bào)， 2004， 10（3）： 14-18

［23］ 陳小容， 雍太文， 楊文鈺， 等. 不同種植模式下茬口對(duì)小麥季土壤養(yǎng)分和酶活性的影響［J］. 中國土壤與肥料，2011（6）： 20-24

Spatial Variability Characteristic and Its Influencing Factors of Tobacco Planting Soil Available N， P and K in Yanyuan County of Sichuan， China

JIN Mingqing1， PENG Yueyue1， WANG Pei2， FENG Lang1，3， CHEN Yulan3，DU Xuanyan1， LI Qiquan1， WANG Changquan1*
（1 College of Resources， Sichuan Agricultural University， Chengdu 611130， China； 2 China Tobacco Sichuan Industrial Co. Ltd.， Chengdu 610000， China； 3. Liangshan District Tobacco Company， Liangshan， Sichuan 615000， China）

In order to provide scientific

for tobacco planting， reasonable fertilization and improvement of tobacco quality， the spatial distribution characteristics of available nutrients of tobacco planting soil in Yanyuan， Sichuan Province were investigated using geostatistics and the effects of the influencing factors on these characteristics were examined by regression analysis. The results indicated： 1） The means of soil available N （AN）， available P （AP） and available K （AK） in the study area are 102.84， 10.18 and 148.10 mg/kg， respectively， and all of them are in the middle level. Their coefficients of variation are 38.66%， 90.08% and 60.65%， respectively， and reached middle variation degree， showing uneven spatial distribution. 2） The exponential model fitted spatial distributions of soil AN and AP well with a moderate spatial dependence， while Gaussian model fitted soil AK well with a high spatial dependence. The spatial variabilities of the nutrients are affected by structural and random factors synchronously. 3） The spatial variabilities of tobacco planting soil AN， AP and AK in Yanyuan have different main controlling factors. The main controlling factors of AN are slope and previous crop. The factors for AP are soil type and previous crop and are altitude and slope for soil AK. Therefore， various influencing factors need to be considered for selection of tobacco planting soils and fertilization.

Tobacco planting soil； Available N， P， K； Spatial variability； Influencing factors； Yanyuan County of Sichuan

S158

10.13758/j.cnki.tr.2016.05.021

四川省教育廳科研項(xiàng)目（16ZB0048）、中國煙草總公司四川省公司科技項(xiàng)目（scyc201402006）、四川中煙工業(yè)有限責(zé)任公司重點(diǎn)項(xiàng)目（130235、130236）資助。

*通訊作者（wcquan@sicau.edu.cn）

金明清（1988—），男，四川郫縣人，碩士，主要從事土壤及土地資源管理方面的研究。E-mail： 448771493@qq.com