四川省盐源县植烟土壤氮磷钾空间变异特征及影响因素①

金明清，彭月月，王 佩，冯 浪，，陈玉蓝，杜宣延，李启权，王昌全*

（1 四川农业大学资源学院，成都 611130；2 四川中烟工业有限责任公司，成都 610000；3 四川省烟草公司凉山州公司，四川凉山 615000）

四川省盐源县植烟土壤氮磷钾空间变异特征及影响因素①

金明清1，彭月月1，王 佩2，冯 浪1，3，陈玉蓝3，杜宣延1，李启权1，王昌全1*

（1 四川农业大学资源学院，成都 611130；2 四川中烟工业有限责任公司，成都 610000；3 四川省烟草公司凉山州公司，四川凉山 615000）

利用地统计学和回归分析方法分析四川省盐源县植烟土壤碱解氮、有效磷、速效钾含量空间变异特征以及不同因素对其空间变异的影响，以期为盐源烤烟种植区划、科学施肥和产质量提升提供科学依据。结果表明：①碱解氮、有效磷和速效钾的均值分别为102.84、10.18、148.10 mg/kg，整体处于中等水平，变异系数分别为38.66%、90.08%、60.65%，达到中等变异程度，空间分布不均；②碱解氮、有效磷空间分布符合指数模型，具有中等空间自相关，速效钾符合高斯模型，具有明显空间自相关，空间变异受结构性因素和随机性因素共同影响；③盐源植烟土壤碱解氮、有效磷和速效钾空间变异的主控因素有所不同，碱解氮为坡度和前茬作物，有效磷为土壤类型和前茬作物，速效钾为海拔和坡度，因此植烟区规划和施肥需要综合考虑不同的影响因素。

植烟土壤；碱解氮；有效磷；速效钾；空间变异；影响因素；四川盐源

土壤养分含量的高低及相互的协调性对烤烟的产量和质量具有十分重要的影响［1］。氮素是影响烤烟产量和品质重要元素之一［2-4］；磷素对烤烟生长发育、生理生化和产质量同样发挥着不可替代的作用［5-7］；钾素是烤烟在生长发育过程中吸收和积累量最大的元素，供应充足是获得优质烟叶的重要条件［8-10］。土壤速效养分是不经转化能迅速被作物直接吸收利用的土壤营养元素［11］，对烤烟的影响更加直接和深刻。近年来，许多学者利用地统计学方法对不同类型植烟区的土壤养分空间变异做了大量的研究，如王宝顺等［12］找出影响重庆彭水坡地植烟土壤速效养分空间变异的主要因素，并掌握坡地植烟土壤养分的空间分布特征及其变异规律；刘逊等［13］采用地统计学方法研究并掌握了湖南湘西植烟土壤有机质的空间分布及影响因素；王育军等［14］采用地统计方法建立云南罗平县植烟土壤肥力适宜性指数的时空变异特征。

四川省盐源县是我省烤烟种植大县，盛产以“清、甜、香”为特色的优质烟叶，在四川优质烟叶生产中占有重要地位。本文利用地统计学和回归分析方法对盐源植烟土壤碱解氮、有效磷、速效钾含量进行分析，重点研究其空间分布特征和影响因素，以期为盐源烤烟种植区划、科学施肥和产质量提升提供科学依据。

1 材料与方法

1.1 研究区概况

盐源县位于四川省西南部，青藏高原东南缘，雅砻江下游西岸，位于100°42′09″ ～ 102°03′44″E、27°06′31 ～ 28°16′31″N，面积为8 398.6 km2，海拔多在2 300 ～ 2 800 m之间，地形复杂，地貌类型多样，是一个山高坡陡、河谷深切、山峦叠障、盆地居中为总特征的典型山区县。受高空西风气流的南支急流和印度洋气流的影响，全县属于典型大陆性气候，四季分明，年温差小，日温差大，年均气温13.2℃，冬春干旱，夏秋雨量集中，年均降水量727.30 mm，雨热同期，光照资源非常丰富；县内土壤类型主要有红壤、水稻土、新积土、黄棕壤及紫色土等。近年来烤烟种植面积约为4 500 hm2，烟叶收购量约为 11 200 t 。

1.2 土样采集与分析

根据盐源县主要土壤类型和烟草种植情况，兼顾代表性与均匀性的基础上，遵循随机、等量、同层次、多点混合的采样原则，利用GPS 定位技术获取每个采样点的坐标，采用网格法和多点混合法采集0 ～20 cm耕层的土壤样品。样品主要于2011年植烟前采集，采集土壤样品317个，采集土壤样品时详细记录其地形地貌、土壤类型、土地利用方式以及施肥状况等环境信息，采样点分布见图1。土样碱解氮、有效磷和速效钾分别采用碱解扩散法、盐酸-氟化铵-钼锑抗比色法和乙酸铵浸提-火焰光度计法测定［15］。

图1 盐源植烟土壤采样点分布图Fig. 1 Sampling sites of tobacco planting soil in Yanyuan

1.3 研究方法与数据处理

由于土壤本身的空间变异性，在进行野外随机区组设计或者随机确定采样点的过程中，可能使获得的试验数据根本不具有统计学意义，而应用地统计学方法可以有效解决这个问题［16］。本研究采用地统计方法对盐源植烟土壤碱解氮、有效磷、速效钾的空间变异特征进行分析，拟合出碱解氮、有效磷、速效钾各向同性时较优的半方差函数模型。半方差函数模型主要根据拟合结果的决定系数和残差的大小来决定。空间变异性强弱可根据块金值与基台值的比值，即块金系数进行划分。块金系数表示由随机部分引起的空间变异性占总体变异的比例，若块金系数＜25%，说明变量有很强的空间自相关性，25% ～ 50% 说明变量有明显的空间自相关，50% ～ 75% 时变量有中等空间自相关，＞75% 时变量空间自相关弱，变异主要由随机因素影响。变程是指变异函数达到基台值所对应的距离，它反映的是属性因子空间自相关范围的大小［17-19］。

地统计分析首先采用GS+ 软件进行半方差分析，以获取碱解氮、有效磷、速效钾各向同性的最优理论半方差模型和参数，再在ArcGIS软件的地统计模块中进行Kring插值，得到研究区土壤碱解氮、有效磷、速效钾的空间分布图。回归分析在SPSS软件中完成，运用回归分析定量揭示各影响因素对植烟土壤碱解氮、有效磷、速效钾空间变异的影响程度，即利用回归方程的校正决定系数来判定各影响因素对土壤碱解氮、有效磷、速效钾含量空间变异程度的独立解释能力［20］，并以此确定研究区碱解氮、有效磷、速效钾空间分布的主控因素。考虑到土壤类型和前茬作物为定性变量，海拔和坡度分级后也属于定性变量，故先采用哑变量对土壤类型、海拔、坡度和前茬作物进行赋值，再以哑变量作为自变量进行回归分析［21］。

1.4 植烟土壤碱解氮、有效磷、速效钾适宜程度划分标准

参考我国主要烟区土壤养分丰缺状况评价结果［22］，综合制定植烟土壤碱解氮、有效磷、速效钾含量适宜程度分级标准（表1）。

表1 植烟土壤碱解氮、有效磷、速效钾含量适宜程度分级标准（mg/kg）Table1 Suitable classification standard of available N， P and K contents of tobacco planting soil

2 结果分析与讨论

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，均值都处于适宜水平。从变异系数来看，碱解氮、有效磷和速效钾的变异系数分别为38.66%、90.08% 和60.65%，均达到中等变异程度，其空间分布不均。40.38% 的样点碱解氮处于适宜水平，30.60% 和16.09% 的样点分别处于低和极低水平；23.03% 的样点有效磷处于适宜水平，超过60%的样点有效磷缺乏；仅17.03% 的样点速效钾处于适宜水平，超过60% 的样点速效钾缺乏。总体来看，盐源植烟土壤碱解氮、有效磷和速效钾含量处于中等偏下水平，部分区域需重视氮肥、磷肥和钾肥的补充。

表2 盐源植烟土壤碱解氮、有效磷、速效钾含量描述性统计（mg/kg）Table2 Descriptive statistics of available N， P and K contents of tobacco planting soil in Yanyuan

2.2 植烟土壤碱解氮、有效磷、速效钾空间结构分析

从各向同性半方差函数拟合结果（表3、图2）来看，盐源植烟土壤碱解氮、有效磷含量分布符合指数模型，速效钾符合高斯模型，拟合决定系数差异较大；速效钾决定系数最大、碱解氮次之，有效磷最小。碱解氮、有效磷和速效钾的块金系数分别为56.31%、55.08% 和46.55%，说明碱解氮和有效磷具有中等空间自相关，速效钾具有明显空间自相关，空间变异受结构性因素和随机性因素共同影响。

表3 盐源土壤碱解氮、有效磷、速效钾含量各向同性半方差函数理论模型及相关参数Table3 Theoretical isotropic semivariogram models and corresponding parameters of soil available N， P and K contents in Yanyuan

图2 盐源植烟土壤碱解氮（A）、有效磷（B）、速效钾（C）含量半方差函数图Fig. 2 Semivariance models for available N（A）， P（B） and K（C） contents of tobacco planting soil in Yanyuan

2.3 植烟土壤碱解氮、有效磷、速效钾空间分布特征分析

根据盐源植烟土壤碱解氮、有效磷、速效钾含量的空间变异特征，利用ArcGIS采用普通克里金法生成碱解氮、有效磷和速效钾的空间分布图（图3）。从图中可以看出，碱解氮、有效磷、速效钾含量空间分布具有明显差异性。碱解氮含量处于适宜水平 100 ～150 mg/kg之间的面积最大，主要分布在东西两侧；处于较低水平65 ～ 100 mg/kg之间的面积次之，主要分布在中部及南部；处于较高水平150 ～ 200 mg/kg之间的面积较小，主要分布西部、北部及东部的部分区域。有效磷含量处于低水平（5 ～ 10 mg/kg）范围基本覆盖全境，仅在西北、东北和东南有小部分区域有效磷含量达到适宜水平10 ～ 20 mg/kg。速效钾含量处于低和适宜水平（80 ～ 150 mg/kg和150 ～ 220 mg/kg）范围内的面积最大，主要分布在中部、西部和南部；处于较高水平（220 ～ 350 mg/kg）范围内的面积较小，主要分布在北部和东北部。各养分含量的空间变异明显，生产上可根据区域土壤养分实际含量进行具有针对性的施肥，避免进一步加剧土壤养分失衡。

2.4 植烟土壤碱解氮、有效磷、速效钾影响因素分析

2.4.1 土壤类型 不同的土壤类型由不同的母质类型、成土过程、发育程度和耕作管理措施等多种因素综合作用而形成，不同土壤类型间养分元素的组成和分布具有明显的差异性。由表4可知，盐源植烟土壤主要以红壤、水稻土和新积土为主。各土壤类型碱解氮含量均值为98.14 ～ 111.87 mg/kg，其中以水稻土最高、红壤最低；各土壤类型的变异程度均为中等变异；红壤、水稻土、新积土均有超过40% 的比例达到适宜程度，紫色土则需要进一步补充氮肥。有效磷含量均值为8.43 ～ 13.32 mg/kg，以紫色土最高、红壤最低；除黄棕壤达到强变异程度外，其余类型土壤达到中等变异程度；各土壤类型有效磷含量适宜度等级主要以极低和低水平为主。紫色土速效钾含量均值以189.35 mg/kg达到最高，黄棕壤以173.79 mg/kg次之，红壤速效钾含量均值最低，为128.27 mg/kg；各类型土壤的变异系数为56.90% ～ 62.60%，达到中等变异程度；从适宜度等级来看，黄棕壤、紫色土达到适宜及其以上等级的比例较红壤、水稻土和新积土高。总体来看，红壤碱解氮、有效磷和速效钾含量均较低，肥力低；水稻土碱解氮含量丰富、有效磷和速效钾含量较高，较肥沃；紫色土有效磷、速效钾含量丰富，土壤类型的差异可为烤烟施肥提供相应参考。

图3 盐源植烟土壤碱解氮、有效磷、速效钾含量空间分布图Fig. 3 Spatial distribution of available N， P and K contents of tobacco planting soil in Yanyuan

表4 盐源植烟土壤不同土壤类型碱解氮、有效磷、速效钾含量的统计特征（mg/kg）Table4 Statistical characteristics of available N， P and K contents of tobacco planting soils with different types in Yanyuan

2.4.2 海拔 不同海拔的土壤环境具有不同的温度、湿度和热量，影响物质的积累与分解程度，进而影响土壤营养元素的含量。从表5可知，盐源植烟区海拔主要在2 000 ～ 2 500 m范围内。海拔在1 500 ～2 000 m范围内碱解氮均值最大，为114.33 mg/kg，海拔在2 000 ～ 2 500 m范围内碱解氮均值最小，其值为98.38 mg/kg；各海拔变异系数在37.00% ～46.64% 之间，达到中等变异程度。有效磷含量均值以海拔2 000 m为分界线，20 00 m以下的有效磷明显比2 000 m以上的高；从变异系数来看，除海拔＜1 500 m的变异系数达到强变异程度，其余的均达到中等变异强度。海拔＜1 500 m的速效钾含量均值最小，为106.57 mg/kg，海拔＞2 500 m速效钾含量均值最大，为174.73 mg/kg；各海拔的变异系数在51.63% ～82.42% 之间，均达到中等程度变异。总体来看，海拔在2 000 ～ 2 500 m范围内，植烟土壤碱解氮、有效磷和速效钾含量均偏低，主要原因可能是研究区种植业主要集中在这一范围内，且耕作强度大、有机肥投入较低、无机肥施用不合理以及养分流失，进而导致了该海拔土壤碱解氮、有效磷、速效钾含量下降，土壤肥力降低。

表5 盐源植烟土壤不同海拔土壤碱解氮、有效磷、速效钾含量的统计特征（mg/kg）Table5 Statistical characteristics of available N， P and K contents of tobacco planting soil with different altitude in Yanyuan

表6 盐源植烟土壤不同坡度土壤碱解氮、有效磷、速效钾含量的统计特征（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可知，盐源植烟区坡度主要集中0 ～ 15° 范围内。各坡度碱解氮含量均值在90.58 ～112.62 mg/kg之间，碱解氮含量与坡度呈反比，坡度越低，碱解氮含量越大；各坡度碱解氮含量变异系数在34.43% ～ 47.18%，达到中等变异程度。有效磷含量随坡度增高而减小，与坡度呈反比，坡度从低到高有效磷含量均值分别为12.36、10.53、9.36、8.91、 8.22 mg/kg；有效磷含量变异系数在83.87% ～115.94% 之间，除坡度 ＞25°达到强变异程度外，其余均为中等变异程度。速效钾含量同样与坡度呈反比，＜2° 均值高达170.64 mg/kg，而 ＞25° 均值仅为110.39 mg/kg；速效钾含量变异系数在52.46% ～80.29% 之间，均达到中等变异程度。碱解氮、有效磷和速效钾含量均与坡度呈反比，主要因为土壤营养物质主要以水为载体进行迁移，坡度会促使土壤水分发生以下行为主的运动，进而导致土壤养分由陡坡向缓坡富集。由于陡坡本身土壤养分缺乏，投入产出比低，耕作难度大，容易发生水土流失，所以在植烟区规划时尽量减少陡坡植烟面积。

2.4.4 前茬作物 轮作方式不同，收获后土壤中残留的有效养分有差别，不同的前茬作物直接影响下茬作物生长的土壤养分［23］。由于盐源烟田冬季以休耕为主，所以前茬作物一般指上一年度所种植的作物。从表7可知，在采样点中盐源烤烟前茬作物以玉米、烤烟、水稻、马铃薯和辣椒为主。碱解氮含量均值大小排序为：水稻＞其他作物＞马铃薯＞玉米＞烤烟＞辣椒，有效磷含量均值大小排序为：烤烟＞玉米＞马铃薯＞辣椒＞其他作物＞水稻，速效钾含量均值大小排序为：其他作物＞马铃薯＞玉米＞烤烟＞水稻＞辣椒，碱解氮、有效磷和速效钾的适宜程度等级分布与各自含量的均值呈现出相似趋势。烤烟、辣椒、马铃薯、水稻、玉米等作物在生长发育过程中对养分的吸收具有差异性，且农民种植不同的作物施肥量具有明显差异，例如这几种作物在采样点的氮肥（按N纯量计算）施用量均值大小关系为：水稻（198.3 kg/hm2）＞马铃薯（184.5 kg/hm2）＞玉米（172.2 kg/hm2）＞辣椒（115.7 kg/hm2）＞烤烟（105.4 kg/hm2），磷肥（按P2O5纯量计算）施用量均值大小关系为：烤烟（75.6 kg/hm2）＞玉米（71.0 kg/hm2）＞马铃薯（66.7 kg/hm2）＞辣椒（54.8 kg/hm2）＞水稻（51.1 kg/hm2），钾肥（按K2O纯量计算）施用量均值大小关系为：烤烟（397.2 kg/hm2）＞马铃薯（313.8 kg/hm2）＞辣椒（243.5 kg/hm2）＞玉米（157.9 kg/hm2）＞水稻（110.6 kg/hm2），因此长期单一耕作就会导致土壤某些营养元素的富集或亏缺，进而出现土壤碱解氮、有效磷、速效钾差异明显的现象。所以合理的轮作制度对维持植烟土壤养分平衡、稳定植烟土壤肥力、进而保证和提高烤烟产质量具有十分重要的意义。

表7 盐源植烟土壤不同前茬作物土壤碱解氮、有效磷、速效钾含量的统计特征（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 不同因素影响程度分析 分别对土壤类型、海拔、坡度和前茬作物进行回归分析，以定量揭示不同因素对盐源植烟土壤碱解氮、有效磷、速效钾空间变异的影响程度。由表8可知，土壤类型、海拔、坡度和前茬作物对碱解氮的独立解释能力分别为5.3%、4.7%、12.4%、10.7%，坡度和前茬作物的独立解释能力高于土壤类型和海拔的独立解释能力，碱解氮的空间变异主要受坡度和前茬作物的影响；各因素对有效磷的独立解释能力分别为13.2%、8.0%、9.1%、12.2%，土壤类型和前茬作物的独立解释能力高于海拔和坡度，有效磷的空间变异主要受土壤类型和前茬作物的影响；土壤类型、海拔、坡度和前茬作物对速效钾的独立解释能力分别为12.7%、17.4%、16.6% 和12.0%，各因素的独立解释能力均较大，速效钾的空间变异受海拔和坡度的影响更加明显。综上可知，结构性因素和随机性因素共同影响盐源植烟土壤碱解氮、有效磷、速效钾的空间变异，其中速效钾的空间变异受结构性因素影响更明显。因此，在研究区进行植烟区划和科学施肥应综合考虑不同的影响因素。

表8 不同影响因素对盐源植烟土壤碱解氮、有效磷、速效钾含量的回归分析结果Table8 Regression analysis of available N， P and K contents of tobacco planting soil with different influencing factors in Yanyuan

3 结论

1） 盐源植烟土壤碱解氮、有效磷和速效钾含量均值皆处于适宜水平。分别有46.69%、64.04%和60.26% 的样点处于缺乏水平，且均达到中等变异程度，空间分布不均。

2） 盐源植烟土壤碱解氮、有效磷含量分布符合指数模型，具有中等空间自相关；速效钾符合高斯模型，具有明显空间自相关，自相关范围差异较大，空间变异受结构性因素和随机性因素共同影响。大部分烟田碱解氮含量处于适宜水平，主要分布在东西两侧；大部分烟田有效磷含量处于低和较低水平，仅在西部和东南的小部分区域达到适宜水平；大部分烟田速效钾含量处于低和较低水平，主要分布在盐源的西部和南部。

3） 盐源植烟土壤碱解氮的主控因素为坡度和前茬作物、有效磷的主控因素为土壤类型和前茬作物，速效钾的主控因素为海拔和坡度，因此在研究区进行植烟区规划和烤烟科学施肥需要考虑不同的影响因素。

［1］ 陈江华， 刘建利， 李志宏. 中国植烟土壤及烟草养分综合管理［M］. 北京： 科学出版社， 2008

［2］ 李志， 史宏志， 刘国顺， 等. 施氮量对皖南砂壤土烤烟碳氮代谢动态变化的影响［J］. 土壤， 2010， 42（1）： 8-13

［3］ 张生杰， 黄元炯， 任庆成， 等. 氮素对不同品种烤烟叶片衰老、光合特性及产量和品质的影响［J］. 应用生态学报， 2010， 21（3）： 668-674

［4］ 刘青丽， 陈阜， 张云贵， 等. 我国西南烟区典型植烟土壤烤烟氮素的吸收规律［J］. 作物学报， 2013， 39（3）：486-493

［5］ 金霞， 赵正雄， 吕芬， 等. 施磷量对烤烟几种生理生化物质含量、赤星病发生及烟叶产质量的影响［J］. 中国烟草学报， 2010， 16（3）： 53-56

［6］ 贾志红， 易建华， 符建国， 等. 磷肥处理对烤烟生长生理及根系构型的影响［J］. 土壤， 2011， 43（3）： 388-391

［7］ 郭燕， 许自成， 毕庆文， 等. 恩施烟区土壤磷素与烤烟磷含量的关系［J］. 土壤通报， 2010， 41（2）： 403-407

［8］ 代晓燕， 郭春燕， 王海波， 等. 钾肥施用方式对豫西烤烟钾含量及产质量的影响［J］. 中国烟草学报， 2013， 18（3）：42-49

［9］ 邬石根， 白厚义， 张超兰， 等. 钾肥和营养调节剂对烤烟含钾量及重金属含量的效应研究［J］. 土壤通报， 2012，43（1）： 190-194

［10］ 郑东方， 许嘉阳， 许自成， 等. 钾肥和腐殖酸互作对烤烟有机钾盐指数的影响［J］. 土壤学报， 2015， 52（3）：169-177

［11］ 黄昌勇， 徐建明. 土壤学［M］. 北京： 中国农业出版社，2010

［12］ 王宝峰， 刘国顺， 叶协锋. 坡地植烟土壤速效养分空间分析［J］. 西南农业学报， 2012， 25（1）： 217-222

［13］ 刘逊， 邓小华， 周米良， 等. 湘西植烟土壤有机质含量分布及其影响因素［J］. 核农学报， 2012， 26（7）： 1 037-1 042

［14］ 王育军， 周冀衡， 孙书斌， 等. 云南省罗平县烟区土壤肥力适宜性评价及养分时空变异特征［J］. 土壤， 2015，47（3）： 515-523

［15］ 张甘霖， 龚子同. 土壤调查实验室分析方法［M］. 北京：科学出版社， 2012

［16］ 刘爱利， 王培法， 丁园圆. 地统计学概论［M］. 北京： 科学出版社， 2012

［17］ 马渝欣， 李徐生， 李德成， 等. 安徽省皖南山区农田表层土壤有机碳空间变异及其影响因素——以宣州区为例［J］. 土壤通报， 2014， 45（6）： 1 424-1 429

［18］ 李启权， 王昌全， 张文江， 等. 基于神经网络模型和地统计学方法的土壤养分空间分布预测［J］. 应用生态学报，2013， 24（2）： 459-466

［19］ 腊贵晓， 顾怀胜， 刘国顺， 等. 喀斯特地区烟田土壤养分的空间变异特征［J］. 水土保持研究， 2012， 19（3）： 48-53

［20］ 顾成军， 史学正， 于东升， 等. 省域土壤有机碳空间分布的主控因子—土壤类型与土地利用比较［J］. 土壤学报，2013， 50（3）： 1-8

［21］ 罗由林， 李启权， 王昌全， 等. 四川省仁寿县土壤有机碳空间分布特征及其主控因素［J］. 中国生态农业学报，2015， 23（1）： 34-42

［22］ 陈江华， 李志宏， 刘建利， 等. 全国主要烟区土壤养分丰缺状况评价［J］. 中国烟草学报， 2004， 10（3）： 14-18

［23］ 陈小容， 雍太文， 杨文钰， 等. 不同种植模式下茬口对小麦季土壤养分和酶活性的影响［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

四川省教育厅科研项目（16ZB0048）、中国烟草总公司四川省公司科技项目（scyc201402006）、四川中烟工业有限责任公司重点项目（130235、130236）资助。

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

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