土壤结构和形态相当复杂，人类活动进一步增强了土壤养分的空间变异性和不确定性，农田土壤养分的空间变异性具有普遍性，它主要受土壤类型、水文地貌、气候条件、耕作方式、作物类型、施肥方式等因素的影响。加强土壤养分空间变异及其影响因素的研究，对环境预测、精准农业和资源管理起着重要的作用，将会直接影响土壤养分的循环和作物的产量。国内外学者对农田土壤养分空间变异性及保护性耕作对土壤养分含量的影响进行了大量研究，但保护性耕作对土壤养分空间变异性的影响研究相对较少。以吉林省梨树县的金山传统耕作和高家免耕土壤养分的空间变异性为例，利用地统计学模型，结合 GIS 技术，研究免耕土壤碱解氮、速效磷、速效钾和有机质的空间变异性。本研究有利于科学合理地制定农田施肥方案，提高养分资源的利用率，促进变量施肥技术的发展，为实现精准施肥奠定基础。
1、 研究区域与研究方法
1.1 研究区域
梨树县位于吉林省西南部，地处东经 123°45′～124°53′，北纬 42°49′～43°46′之间，总面积达 4 209km2，属中型县。梨树县居于东北大陆腹地，属于温带半湿润大陆性季风气候，冬季漫长寒冷，夏季炎热多雨，春季风沙较大，秋季气候凉爽。年平均气温为5.8℃，年湿润度为 0.7。0℃以上积温为 3 477.9℃·d，平均全年日照时数为 2 698.5 h，无霜期平均为152 d，年降水量平均为 577.2 mm，主要集中在 6、7、8 三个月。热量比较充足，雨量比较适中，且热量与雨量同期，可满足一年一熟农作物的生长需要。梨树县地处东辽河下游，东、北、西北为东辽河所环绕，地势东南高西北低。东南部为低山丘陵，系长白山支脉大黑山余脉。中部为波状平原，北部系低洼平原，地势平坦，土质肥沃。全县土地面积42.1 万 hm2， 其中耕地面积23.4万hm2，占55.6%，属于以农业生产为主的平原县。全县 27个乡镇，337 个行政村，总人口 86 万人，农业人口64 万人。
本文在梨树县境内选择高家村和金山村两块试验田作为研究对象，其中高家试验田是玉米覆盖保护性耕作方式(2007 年开始)，金山村试验田是传统耕作方式(图 1)。两块地的自然属性相似，都是多年种植玉米的农田，因此可用来对比保护性耕作与传统耕作对土壤养分空间变异性的影响。
1.2 数据采集与测试
2009 年 10 月(玉米收获之后)分别在高家村、金山村用手动脚踏式直压取土器进行土壤样品的采集。运用“S”型 5 点采样法采样，网格采样间距是40 m×40 m，采样深度为 0～20 cm。高家和金山均采集 60 个土壤样品，将从野外取回的土样进行风干、磨细、过筛、混匀、装袋处理，并对土壤养分中的碱解氮、速效磷、速效钾、有机质进行测试。采用重铬酸钾容量法测定土壤中有机质的含量,碱解氮的测定采用碱解扩散法，速效磷的测定采用 0.5mol/L NaHCO3浸提 - 钼锑抗比色法，速效钾的测定采用 1 mol/L NH4OAc 浸提 - 火焰光度法。
1.3 研究方法
运用 SPSS15.0 对每个地块的土壤养分数据进行常规统计分析，描述地统计学分析的变量包括标准差、变异系数、偏度系数和峰度系数。变异系数(C.V)表示土壤特性空间变异性的大小，变异系数 C.V≤0.1 为弱变异性，0.1<C.V<1.0 为中等变异性，C.V≥1.0 为强变异性。
利用 GS+ 软件计算出高家和金山土壤养分的半方差模型，分析其决定系数，块金值(Nugget)与基台值(Sill)的比值和变程，进而分析其空间变异性。块金值与基台值的比值可以表明土壤养分的空间相关程度。该比值 <25%，说明系统具有较强的空间相关性，由自然因素引起的空间变异性较高，具有稳定的结构性；比值在 25%～75%之间，表明系统具有中等相关性，说明由随机因素引起的空间变异性较高；比值 >75%，说明系统相关性很弱；如果该比值接近 1，则说明该变量在整个尺度上具有随机性。最大相关距离即变程(Range)是指变异函数达到基台值时，采样点的间隔距离，它表明了土壤要素的空间自相关范围，表示在某种观测尺度下，空间相关性的作用范围。在变程范围内，样点间的距离越小，其相似性越强，即空间相关性越大。
常用的空间插值方法有反距离加权插值、最近邻点插值法、克里格插值方法、多元回归法、径向基函数法、线性插值三角网法、局部多项式法等，其中克里格插值方法是地学中常用的一种插值方法，是根据变异函数模型发展起来的一系列地统计的空间插值方法。本文用 ArcGIS9.3 的普通克里格插值对土壤养分进行插值，寻找土壤养分的空间变异规律。
2、 结果与分析
2.1 描述性统计特征值
土壤养分的描述性统计分析是建立空间变异性模型的前提，运用 SPSS15.0 对高家和金山土壤养分数据进行描述性统计分析(表 1)。
从变异系数对比可以看出，除金山传统性耕作模式下土壤速效钾的变异属于弱变异强度外，其他农田土壤养分的变异都属于中等变异强度。
不管是保护性耕作还是传统性耕作方式，速效磷的空间变异最大，其次是有机质，再其次是碱解氮，速效钾的空间变异最小。整体上看，保护性耕作土壤养分空间变异性较大，而传统性耕作土壤养分空间变异性较小，主要是因为传统性耕作方式受人类干扰较强，人类通过对农田大量施用化肥、进行秸秆回收等措施控制土壤养分，经过长期的干扰与控制，土壤养分的空间变异性逐渐减小，而保护性耕作主要利用秸秆还田、复耕等措施使土壤接近于自然状态，土壤养分的空间变异较大。
2.2 空间结构性和相关性分析
利用 GS+ 软件计算得到高家保护性耕作和金山传统性耕作土壤养分的半方差模型及其拟合参数(图 2，表 2)。
高家和金山土壤养分的块金 / 基台值均小于25%，说明土壤养分的空间变异都具有较强的空间相关性。保护性耕作方式下，速效钾和碱解氮的空间结构性最强，有机质居中，速效磷的空间结构性最差。传统耕作方式下，碱解氮的空间结构性最强，有机质和速效磷居中，速效钾的空间结构性最差。
从保护性耕作和传统性耕作土壤养分半方差模型的对比可以看出(图 2，表 2)，保护性耕作模式下土壤养分的块金 / 基台值均小于传统性耕作模式下土壤养分的块金 / 基台值，表明保护性耕作模式下土壤养分比传统性耕作模式下土壤养分的结构性强，具有较强的空间自相关性，说明保护性耕作模式下土壤养分受成土母质、地形地貌、气候等自然因素影响较大。传统性耕作模式下，受施肥、作物品种、管理方式等人为随机因素引起的空间变异性程度较大，结构性相对较差，空间自相关较弱，这也可以从最大变程看出，除速效磷外，传统性耕作模式下的土壤养分的最大变程均小于保护性耕作模式下土壤养分的最大变程。
2.3 土壤养分的空间变异性
精准农业土壤养分空间变异图的作用是根据土壤养分空间差异，划分管理小区，指导精准施肥。
本文使用 ArcGIS 9.3 软件，对高家保护性耕作模式和金山传统性耕作模式的土壤养分进行克里格插值，得到土壤养分的空间变异图(图 3，图 4)。本文将土壤养分含量分为 10 个等级，分析对比不同耕作方式下农田土壤养分空间变异性的差异性。
对比分析图 3、图 4 可以看出，高家保护性耕作模式下土壤养分的空间变异结构性较强，空间变异程度大，呈团块状分布。而金山传统性耕作模式下土壤养分的空间变异结构性相对较差，空间变异程度小，呈垂直于垄向的条带状分布，这是由于这块地的垄向是西北 - 东南走向，不同农户间施用化肥量、作物品种、管理方式等不同造成的。
总之，保护性耕作受自然因素影响较大，而传统性耕作受人为随机因素影响较大。
3、 结论与讨论
3.1 结论
(1)保护性耕作农田土壤养分的空间变异都属于中等变异强度，保护性耕作土壤养分空间变异性明显大于传统性耕作土壤养分空间变异性。
(2)保护性耕作农田土壤养分具有较强的空间相关性，其中速效钾和碱解氮的空间结构性最强，有机质居中，速效磷的空间结构性较差。保护性耕作模式下土壤养分比传统性耕作模式下土壤养分的结构性强。除速效磷外，保护性耕作模式下土壤养分的最大变程大于传统性耕作模式下土壤养分的最大变程。
(3)保护性耕作模式下土壤养分的空间变异程度大，呈团块状分布。而传统性耕作模式下土壤养分的空间变异程度小，呈垂直于垄向的条带状分布。保护性耕作受自然因素影响较大，而传统性耕作受人为随机因素影响较大。
3.2 讨论
本文采用描述性统计特征值的变异系数和地统计学的半方差模型对比分析了不同耕作方式下土壤养分的空间变异性。保护性耕作土壤养分的变异系数都大于传统性耕作的变异系数，而半方差模型中的块金 / 基台值正相反，说明变异系数和块金 / 基台值有一定的相关性，变异系数越大，半方差模型的块金 / 基台值越小，即空间变异性越大，空间结构性越强。但这种相关性具有一定的适用范围，如保护性耕作变异系数表明，速效磷的空间变异最大，其次是有机质，再次是碱解氮，速效钾的空间变异最小；而半方差模型表明，速效钾的空间结构性最强，其次是碱解氮，再次是有机质，速效磷的空间结构性最差，得出的结论应是空间变异性越大，空间结构性越差，得出完全相反的两种结论。
古典的 Fisher 统计法和地统计学方法是研究土壤养分空间变异的两种主要方法。Fisher 统计法的统计原理是假设研究的变量为纯随机变量，样本之间是完全独立。通过计算土壤养分的均值、标准差、变异系数以及进行显著性检验来描述土壤特性的空间变异。该方法在土壤学研究工作中取得了一定的成功，但由于其基本上是定性的描述，只能概括土壤特性变化的全貌，不能把土壤局部的空间变异特性反映出来。国外的研究表明土壤特性也不像 Fisher 原理假设的那样属于纯随机变量，而是在一定的范围内存在空间上的相关性。地统计学是既考虑样本间的结构性，又考虑样本间的随机性的一种空间统计方法，具有较强的优越性，弥补了经典统计方法上的缺陷。因此地统计学的半方差模型比变异系数更能说明土壤养分的空间变异性，是研究土壤养分空间变异性的一种主要方法。
保护性耕作模式下土壤养分的空间变异大于传统性耕作模式，而且结构性和相关性较强，说明精准施肥更适用于保护性耕作模式。保护性耕作和精准农业都是依靠先进的现代科学技术，通过农业机械化、规模化生产，最大限度地优化使用各项农业投入(如化肥、农药、水、种子和其他方面的投入量)以获取最高产量和最大经济效益，同时减少化学物质使用，保护农业生态环境，保护土地等自然资源，实现农业的可持续发展，因此二者有集成的基础。保护性耕作和精准农业的有效集成是未来农业的主要发展趋势之一。
