06精準养分管理-基于GIS的土壤养分分区管理模型研究.docVIP

精准养分管理 基于GIS的土壤养分分区管理模型研究 白由路 金继运 杨俐苹 梁鸣早 中国农科院土壤肥料研究所 　　施用化学肥料是农业生产中投入资金最多的项目之一，但是，化肥的当季利用率还不足50%，氮素化肥的利用率更低[1]。其原因除施用方法不当外，也与土壤养分的空间变异密切相关 [2~4]。国内外对于测土推荐施肥都做了大量的研究[5~6]， 但是，大部分是基于大田平均值，或基于大田混合样品进行的[5]。近年来，随着GPS、GIS和地统计学等方法应用于土壤领域[7]，土壤特性的空间变异越来越受到人们的重视，特别在土壤养分管理上，以精准农业（precision agriculture）为中心的养分管理正成为土壤养分管理的热点[8~9]。但是，我国目前农业生产管理主要以农户为单位，地块面积较小。再由于农业生产特别是粮食生产的比较效益低下，以农户地块为单位的测土推荐施肥难以进行。以大区域或大面积混合样品的测定结果进行施肥推荐时，推荐量不能适合所有地点，造成推荐偏差[6]。　　本研究结合我国农村的实际，应用地理信息系统，研究区域土壤养分分区管理的模式，旨在为充分发挥土壤的生产潜力、提高肥料利用率、保护生态环境提供有关的科学理论依据和发展模式。　　1 材料与方法　　1.1 地点及取样　　试验区在河北省辛集市马兰村进行，面积为173.4ha，地势平坦，土壤为潮土，质地为中壤，主要种植制度为小麦-玉米轮作。2000年6月初在小麦刚收获后，用差分式GPS在试验区用100m的网格采取土壤耕层样品，风干备用[10]。　　1.2 分析项目与方法　　用土壤养分状况系统研究法（ASI）[11]分析土壤养分状况。土壤速效磷(P)、钾(K)、铜(Cu)、铁(Fe)、锰(Mn)和锌(Zn)采用ASI联合浸提剂（0.25mol/L NaHCO3-0.01 mol/L EDTA-0.01mol/L NH4F）浸提；速效硼(B)和硫(S)用0.08mol/L的过磷酸钙溶液浸提；速效钙 (Ca)、镁(Mg)和铵态氮(N)用1mol/L的KCl浸提；有机质用0.2 mol/L NaOH - 0.01 mol/L EDTA - 2% 甲醇溶液提取。有机质和非金属元素用比色法测定，金属元素用原子吸收分光光度计测定。　　1.3 数据处理平台　　GIS平台为ESRI公司的Arc/Info7.0和Arcview3.0[12]。　　2 结果与分析　　2.1 土壤养分分区管理的概念模型　　在我国目前以农户为单位的农业生产体制下，地块既小又不连片，每个地块面积仅0.1~0.2ha，若以该面积为取样单位，会大幅度提高单位面积测土推荐施肥的费用，难以有效推广。当以较大面积为测土施肥推荐单元时，由于空间变异，测土施肥的推荐量不具有较强的代表性，所以，在建立土壤养分分区管理模型时，必须克服上述问题。本研究以土壤特性相对均匀的村为单位，建立以地块或农户为单位的土壤养分分区管理模型(图1)。该模型以土壤网格取样和ASI的土壤养分推荐模型为基础，经地理信息系统的一系列处理，形成土壤养分管理图，由于地理信息系统中的数据库管理功能，也可将土壤养分管理方法以表的形式输出。在该管理模型中，管理单元是人为划分的，在较小的规模下，管理单元可以是农户的地块，也可以是生产小组的地块。在比较宏观的情况下，也可以以村或县为单位。 (图：图1 土壤养分管理的概念模型) 　　2.2 模型实现的基本步骤　　为了实现土壤养分的分区管理，以Arc/Info 7.1为基本平台[12]，对土壤养分管理进行了探索。在该模型中，所需的基本数据与图层为取样点图、区域边界图、分区边界图和基于ASI养分管理的作物类型、产量目标及土壤养分分级。在Arc/Info平台上需建立11个文件，进行9个步骤的处理。同时，还需要在数据库平台或EXCEL等软件上进行必要的数据处理过程。如果没有GPS定位工具，可采用相对坐标。在较大比例尺的地图上进行数字化也可以得到相同的效果。有高精度GPS定位工具时，田间图包括地界图都可用GPS定位，然后导入GIS中，不必进行数字化。Kriging插值是该过程的关键步骤之一。可根据养分管理的精度进行插值[13]。养分分类是施肥推荐准确程度的重要依据，养分测定方法要与施肥推荐方法一致。为了便于管理和实施，在一个管理单元内，只能有一种施肥量。但是，由于土壤养分的空间变异性，在大多数情况下，一个管理单元内会有多种养分等级和施肥量，需要进行必要的归类，可供选择的归类方法有：　　①面积优先法：即在一个管理单元内，用面积最大的某一类养分的值代表该单元的养分值。　　②加权平均法：即在一个管理单元内，用面积作权重，对养分进行加权平均，用平均值代表该单元的养分值进行施肥