隆化县耕层土壤有机质的空间结构分析.docVIP

隆化县耕层土壤有机质的空间结构分析 　　【摘 要】利用ArcGIS10.0的地统计分析工具，对隆化县耕层土壤（0～20cm）有机质进行了探索性空间分析、半变异函数模型分析和地统计插值。揭示了隆化县土壤有机质的空间变异特征，和分布规律、统计了个级别分布面积。 　　【关键字】耕地；耕层土壤；有机质；空间结构分析 　　1 研究区概况 　　隆化县位于河北省北部，地处冀北山地，东经116°47′45″～118°19′ 07″，北纬41°08′48″～41°50′09″。海拔390～1800m，整个地势北高南低，属大陆性中温带季风型半湿润气候区，年均温度2～8℃；无霜期天90～140d；年降水量450～550mm。全县总面积5460km2，其中耕地面积57334hm2，土壤类型为棕壤、褐土、潮土、草甸土。土壤质地以砂壤质、轻壤质和中壤质为主。主栽农作物为玉米、水稻和杂粮，是河北省粮食生产大县之一。 　　自1981年第二次全国土壤普查以来，农业生产条件和种植结构发生了很大变化，随着农作物良种的普及和化肥施用量的增加，单位面积农作物产量大幅度提高，农田土壤养分也在发生着变化。土壤有机质是土壤肥力的重要指标，分析其含量和分布特征，可以为土壤改良和合理施肥提供科学依据。 　　2 材料与方法 　　2.1 样品采集及测试 　　根据土壤类型，采用GPS定位，于2009年采集耕层（（0～20cm））土壤样品6010个，平均每5～10hm2采取一个土壤样品。具体方法为在距每一样点10～20 m 的范围内采集耕层各5～15 份土样，混合均匀后用四分法再留取1 kg 土样风干，磨细并全部通过0.125mm 筛的土样，土壤有机质采用重铬酸钾氧化—硫酸溶液—油浴法测定。 　　2.2 数据处理 　　以ESRI 公司的ArcGIS 10.0为软件平台，利用其地统计（Geostatistics）模块，对土壤测试数据养分进行了探索性数据分析及空间插值。包括对土壤采样测试数据作直方图分析、QQ正态图分析和趋势分析，检验数据分布特征，寻找异常值；进行分析全局趋势，确定空间插值方法和参数；进行地统计插值，形成连续的预测表面。 　　3 结果与分析 　　3.1 耕层土壤有机质分布特征 　　直方图（图1）分析是对采样数据按一定的分级方案进行分级，统计采样点落入各个级别中的个数，并通过柱状图表现出来，以直观地反映采样数据的分布特征、总体规律，用来检验数据分布。 探索性数据分析的直方图窗口，还显示了耕层土壤有机质一些基本统计信息（见表1）。 　　由表1可见，全县耕层土壤有机质含量的变化范围在0.2～39.9g/kg 之间，平均含量为18.216g/kg ，标准差为6.479，变异系数为37.30% ，属中等变异强度。有机质1/4分位数为13.3； 3/4分位数为22.9，中值为17.3，略低于平均值。 　　峰度（Kurtosis）是用于描述数据分布高度的指标，正态分布的峰度等于3，隆化县耕层土壤有机质数据的峰度为 2.852，很接近正态分布。 　　偏态（Skewness）是用于描述数据分布左右对称性的指标，正态分布的偏态等于0. 隆化县耕层土壤有机质数据的偏态为0.40185，数据低值部分较集中，稍呈正偏态分布。 　　3.2 正态QQ piot图分析 　　正态QQ piot图 提供了另一种度量数据正态分布的方法，利用QQ piot图，可以将现有数据的分布与标准正态分布对比，入股数据越接近一条直线，则它越接近于服从正态分布。 　　从耕层土壤有机质的正态QQ piot图（图2）可见，隆化县耕层土壤有机质测试数据服从正态分布。 　　3.3 数据趋势面分析 　　趋势面分析是根据样点数据拟合一个数学曲面，用该数学曲面来反映空间分布的变化情况，它和可分为趋势面和偏差两部分，其中趋势面反映了土壤测试数据总体的变化趋势，受全局性、大范围的因素影响，准确识别和量化全局趋势，在ARCGIS地统计建模时可以方便的剔除全局趋势，从而更准确地模拟短程随机变异。具体做法是用XY坐标表示样点位置，Z轴表示测试值。将Z轴数值分别投影到XZ平面和YZ平面作散点图，然后用多项式拟合投影平面散点图，根据形成的曲线，识别数据是否存在全局趋势。投影面上的深色散点上的线表示南—北向的全局性趋势应变化，投影面上的浅色散点上的线表示的是东—西向全局性的趋势效应变化。 　　从图3可见，隆化县土壤有机质存在西部高于东部的趋势，旋转坐标得到拟合的趋势图（图4），显示有机质含量呈西南至东北逐渐降低的趋势的全局。 　　3.3.1 土壤养分空间变异 　　（1）土壤有机质的结构分析。地统计工具向导根据不同空间位置上有机质含量的分析数据，计算实际半方差值，并绘制变异函数曲线图，这是空间变异分析的基本步