13精准养分管理-土壤速效养分空间变异特征及评价.docVIP

精准养分管理 土壤速效养分空间变异特征及评价 姜城 杨俐苹 金继运 中国农业科学院土壤肥料研究所 　　空间变异是土壤本身存在的一种自然特性，认识土壤空间变异对于评价和有效地利用土壤都是十分重要的[1]。从理论上说，如果一个观测点能为其邻近的点位提供一定的信息，我们就可以说在它们之间存在着一定的空间依赖性（Spatially Dependent）[2]。处理具有空间依赖性的变量的方法主要是基于60年代初法国学者G.Matheron所创立区域性变量理论[3,4,5]。70年代后，国内外学者开始将该方法应用于土壤性质的研究，但对于农田土壤速效养分空间变异的研究却鲜有报道[6~9]。另外，计算技术的落后也极大地限制了这一研究的进一步深入[10]。　　本文拟在前人研究的基础上，结合地理信息系统和土壤养分系统研究法，探讨在一定农业生产条件下的农田土壤养分空间变异特征，以揭示土壤速效养分的空间变异规律，为实现土壤养分的精确管理提供理论基础。　　1 材料与方法　　1.1 样品的采集和分析 　　试验在邯郸陈刘营村54.8ha棉花田上进行，土壤为褐土，质地砂壤，前作为玉米，取土时间1999年2月10日，取土前对当地的基本农业生产条件（包括前作作物品种、产量水平、施肥习惯、水利条件等）进行详细调查。　　利用网格取样法采集土壤样本98个，取样间隔80m。方法为在以网格点为圆心10m为半径的范围内采集10钻0~20cm的耕层土壤组成代表该点的混合样本（Composite Sample）。取样点分布如图1所示。　　土壤样品的分析采用国际农化服务中心的ASI方法测定[11]。土壤样品用容量法取样。　　 (图：图1　陈刘营棉田网格取样示意图) 　　1.2 数据处理　　本文的数据分析采用Fisher统计学和地统计学结合的方法，涉及地统计学的主要有半方差函数、半方差模型和半方差图。　　1.2.1 半方差函数　　是描述土壤性质空间变异的一个函数，即：　　 　　式中：区域化变量Z(x)是一个随机变量，它是某区域内x位置上记录的不同的Z值。h为样本的位置间距矢量，或称为位差(lag)。r(h)即为只依赖于点与点之间距离h的半方差。　　1.2.2 半方差模型　　半方差的理论模型有多种，常用的有球形模型、线性模型、指数模型及Gauss模型等[12]，本文主要采用球形模型，其模型如下：　　 　　 　　式中：C0为块金系数（间距为0时的半方差），C0+C是基台值（半方差函数随间距递增到一定程度出现的平稳值），C为结构方差（基台值与块金系数的差值），a为变程（半方差达到基台值的样本间距）。基台值近似于采样方差S2，所以间距小于变程的样点是空间相关的，而间隔大于变程的样点则无空间相关性。因此，变程就理所当然地又称为空间相关域。　　区域化变量和土壤特性几乎总是呈现出空间上的依赖性，这种空间依赖性可以用半方差图定量化。半方差图是地统计学解释土壤空间变异性结构的基础，它的精确估计是成功的空间内插和栅格地图制作的关键[7]。　　2 结果与分析　　2.1 土壤速效养分的一般分布特征　　图2是棉田土壤六种元素速效养分的频数分布图。从图中可以看出不同土壤养分含量的变异规律。虽然在各种元素之间存在着一定的差异，但几乎所有元素的频数分布都呈对称或不完全对称的单峰分布，它们的土壤测定结果在平均值附近分布较多，远离平均值出现的频率较低。另外，Fe、Zn和N的分布范围相对较宽，从频数分布状况来说，这些元素变异较大，不容易统一管理。K、P和S的分布则相对集中，K大部分集中在58.648mg/L~93.834mg/L之间，S集中在41.389mg/L~57.805mg/L之间，而P则明显地集中在15.800mg/L~47.300mg/L之间，在这个范围之外分布的样点很少。　　 (图：图2　土壤速效养分频数分布图) 　　对所获得的养分数据赋予相对的地理坐标，使之具有地理属性，然后利用地理信息系统MapGIS稀疏网格加密和Kriging空间内插的功能作出各种养分的分布图，再与地块图进行图层叠加后就可以明确各种养分在地块上的分布状况。图3为地块土壤速效磷和速效钾的分布图。通过地理信息系统的统计模块进行面积分析，结果表明，土壤速效磷含量〈25mg/L的占总面积的26.97%，25 mg/L ~40 mg/L的占62.44%，40 mg/L ~60 mg/L的占5.70%，＞60 mg/L的占4.89%，说明绝大多数土壤速效磷处于中等偏上水平，极高或极低区域所占的比例相当小。土壤速效钾含量〈78.2 mg/L的占总面积37.59%，78.2 mg/L ~117.3 mg/L的占55.27%，＞117.3 m