土壤养分的空间变异规律可作为土壤资源合理利用与改良、精准施肥和环境物质修复的基本依据
[, , ]。土壤的变异来源可分为系统变异和随机变异,其中前者可认为是由母质、气候、水文、地形、生物等相对稳定的因素产生，后者可认为是由不同土地利用方式和行为，以及取样、分析等的误差产生
[, , ]。大量的研究表明地统计学方法能较好地揭示土壤异质性机理
[, , , , ]。近几年,国内外学者对不同尺度( 乡镇级、村级、田间尺度和微尺度) 的土壤性质( 土壤水分、盐分、养分、理化性质及生物学性质等)进行了广泛研究,如，连纲等通过对黄土高原县域土壤养分空间变异分析后，认为不同土地利用类型的养分含量差异显著，水田和川地的有机质和全氮含量较高，全磷含量以梯田最高。高义民等通过对黄土区村级农田土壤养分空间变异特征分析，认为施肥是导致该村土壤养分空间变异性的主要原因。 张博等通过对江苏滩涂围垦区土壤养分空间变异分析得出该区域应多施有机肥和氮肥，少施或不施钾肥和磷肥的结论。这些研究都在一定程度上揭示了土壤有机质、水解氮、速效磷、速效钾等土壤因子在各自环境区内的空间变异规律
[, , , , , , , ]。在方法论的研究上，GIS技术被引入进行景观空间异质性分析，并将景观斑块结构、形态、数量的变化特征，作为景观连通性、稳定性和破碎化的重要依据
[, , , , ]。但是，运用地统计学方法探讨作物种植前后土壤养分空间变异性的文献则鲜有发现。本文以浦江县作物土壤试验数据为基础，结合多年的农户调查数据，以地统计学方法为工具，探讨作物种植前后土壤有机质及其养分因子的空间演化规律，为不同区域作物养分合理实施与种植结构的调整进行分区指导提供理论依据。
浙江省浦江县位于浙江省中部偏西，金衢盆地的北缘，是浦阳、壶源两江之发源地。地理坐标界于东经119°42′—120° 07′、北纬29°21′—29°41′。东西宽39.25 km，南北长36.5 km。全县土地面积为907.6637km
2。境内地貌分西北山区和浦江盆地两大类型。西北部山区属浙中丘陵山地地貌，以低、高丘为主。浦江盆地属于金衢丘陵盆地地貌。以平畈、岗地、低丘为主。根据海拔高度和农业利用特点可将浦江的地形分这平原区(海拔<150 m)、半山区(150 m<海拔< 500 m)、中低山区(海拔＞500 m)，符合“七山二水一分田”中国东部地形特征。
浦江县属亚热带季风气候区，四季分明，气温适中，光照充足，雨量丰沛，自然资源丰富。由于农业用地分别邻近浦阳江、壶源江和大陈江，因而又形成3种小气候带，分别是浦江盆地温暖少雨区(Ⅰ区)、壶源江谷温和次多雨区(Ⅱ区)、中低山温凉多雨区(Ⅲ区)(见)。各气候带特征见。
浦江县区域经济发展水平略高于全国平均水平，随着产业扩张速度的不断加快，耕地保护的压力也面临前所未有的挑战如何优化区域资源是亟待解决的重要问题。本研究正是在欧盟项目IRMLA (综合农业资源管理与土地利用分析)成果的基础上，以正在执行的国家基金为支助，以连续10a的农户土地利用与管理的调查数据和土壤、作物的试验数据为基础，探讨县域农地资源的优化利用。
(1)调查方法
在浦江县选取了107个代表性农户进行调查，农户的选取采取点、线、面相结合的方式。面：尽量做到覆盖面广。将全县分3片：平原、半山区、中低山区，并兼顾到全县的每个乡镇。同时把农作物种类和种植面积、各地气候差异、典型产品和重点产业相结合综合考虑，尽量做到平衡。线：以公路为一线，池塘、沟渠为二线，山间小路为三线。点：除考虑各种不同的作物外，还兼顾好、中、差3类生产水平的农户。具体调查点分布见。在调查人员的安排上，采取乡镇从事农业技术推广专业人员和在校教师、研究生相结合的方式，根据选取的农户、统一标准、逐户调查。
(2)调查内容
在调查之前先设计详细的调查表，与本研究相关的主要内容包括：1)农户家庭土地利用方式与农业生产管理，包括养分投入，特别是N投入；2)农户耕种田块的基本信息，如：种植面积、规模化程度、地形条件、地理位置、土壤性质等，为土地利用方式的空间分析奠定基础。
(1)采样点的选取及农户投入
根据种植作物的不同从107户农户家庭的责任田中取样183个样点(即183块田的采样)，其中，水稻及其它粮油作物共69块，葡萄36块，蔬菜40块，桃形李及其它水果类31块，其它作物7块，各类作物农户平均投入调查结果见。采样地块还需考虑坡度、坡向、相邻性等特征，并做好地理特征的相关信息采集。各样点见。
(2) 土壤取样方法
1)采样时间 第1次取样为作物播种前(3月3日—12日，共10 d)，第2次取样为作物收获后、施肥前(11月25日到12月4日共10 d)，这一取样时间是综合了该区域的季节特征和农事活动的特点选择的。两次取样的田块完全相同。
2)采样部位和深度 根据耕层厚度，确定采样深度，一般取样深度0—20 cm。
3)采样方法、数量 采用多点混合土样采集方法，每个混合农化土样由20个样点组成。每个点的取土深度及重量尽量均匀一致，土样上层和下层的比例基本相同。样品处理、储存等过程没接触金属器具和橡胶制品。每个混合样品一般取1kg左右。
(3) 试验分析方法
样品主要参照鲁如坤的《土壤农业化学分析方法》中的扩散吸收法对土壤有效N进行测定
[]。本文所指的氮、磷、钾分别是指碱解氮、速效磷和速效钾。
(1)空间自相关分析可检验土壤养分空间相关性的存在与否(即是否存在区域的集聚或集中)，一般用全域Moran′s I指数、局域Moran′s I指数表达。本文利用ArcGIS空间统计模块中的全域Moran′s I，对县域内两个时点因子空间自相关性进行分析。计算模型已在先前发表的论文中作了详细阐述，不再赘述
[]。
(2)半方差函数是地统计学的基石，其块金值、基台值和变程等参数能表达区域化变量在一定尺度上的空间变异和相关程度，是Kriging精确插值的基础，也行预测和模拟的关键
[, , , , , , , , , ]。
模型是否最优，主要看拟合的理论模型与实际变异函数的残差平方和
RSS 和决定系数R
2的大小，其中，残差平方和是实际变异函数值γ(h
i)和理论模型变异函数值γ^(h
i)之和:
式中,
RSS 是回归精度的重要参数，其取值愈小，说明实际观测值与回归线靠近，拟合曲线与实际配合愈好。
但是残差平方和受因变量和自变量绝对值大小的影响，
RSS 值大的拟合模型不一定不及
RSS 值小的拟合模型，但是决定系数 R
2 就不存在这个问题。决定系数 R
2, 也称拟合度，其计算公式为:
式中,γ(h
i)为实际变异函数值，γ ^(h
i)为理论模型计算出的变异函数值，γ(h)为实际变异函数平均值。
本文利用单样本K-S检验对各样本因子分别进行检验，发现通过对数转换后各因子均较好地符合正态分布。可以进行Moran′s I指数和Kriging插值分析。
不同养分在不同时点的空间自相关性与步长的总变化趋势基本一致，在作物种植前后均呈现波动性，说明随着尺度的变化，土壤养分之间存在一定的空间关联，但又各具特点(这一点与高义民、史利江等学者得出的结果基本一致
[, , ]):(1)各因子的振幅在步长﹤1000 m时均达到最高点，并随着步长的加大而逐渐衰减，说明区域越近，作物种植结构越相似，养分投入越接近。而随着步长的加大，范围的加大，种植多样性，引起投入的多样化，但大范围的多样化投入又引致相关性程度加大。(2)作物收获后，土壤氮和磷的负振幅超过作物种植前，说明不同农户作物养分投入的差异，导致的相关性增强是由于随机变异引起的，而经过一段时间土壤表层的氮、磷含量趋于稳定的原因是由于土壤本身的吸附与消化能力作用的结果，它属系统变异。(3)作物收获后钾的自相关强度明显弱于作物种植前，而且，随着步长的加大，其相关性增强。前者与不同区域人们对钾肥投入的重视程度不一致有关，后者则说明钾的自相关性与土壤因子系统变异的关系更加密切。具体见。
为分析土壤因子间的交互相关性，论文运用SPSS软件进行Pearson相关系数分析，并在0.01 的置信水平下进行
T检验。得出：(1)作物种植前，在 置信0.01水平下，碱解氮-1与有机质-1有较强相
关性，这一点也与史利江等学者得出的结果相一致
[]。而碱解氮-2与有机质-1之间的相关性较前两者明显减弱。说明有机质的含量与碱解氮间的关系紧密，但前者转化为后者需要一定的时间。在0.05置信水平下，碱解氮-2与速效磷-2间存在一定的相关性，与速效磷-1、速效钾之间基本上不相关，说明氮、磷投入，农户有配施的意识，但对钾的配施意识不强，这与农户调查结果相符合。(2)在0.01置信水平下，速效磷-1与速效磷-2之间存在极显著的相关性，说明农户对磷的投入相对稳定且其投入磷与土壤内磷素间成较为固定的比例。(3)速效磷与速效钾之间在种植前后都具很强的相关性，这可能与调查中65.4%的农户喜欢使用复合肥有关。具体见。
作物种植前、收获后的数据表明：(1)自相关系数，种植前碱解氮、有机质、速效磷和速效钾4种因子的自相关系数均< 25%，且均在11.4% —15.8%之间，表现出强烈的空间相关。收获后，它们的自相关系数均增大，且除碱解氮外，均处于中等空间相关程度。(2)变程，种植前，碱解氮和有机质的变程最小，均为720.00 m，最大的速效磷为1830.00 m。收获后，各因子的变程普遍加大，最小的速效磷也达到了6000.00 m。(3)拟合模型，种植前各土壤因子均符合指数同性模型，拟合效果较好。收获后，拟合模型有差异，根据
R2，有机质、碱解氮、速效磷拟合效果也好于前者()。
通过对模型的精度进行分析，各土壤因子空间分布基本符合最优模型的条件，即标准平均值最接近0，均方根预测误差最接近平均标准误差，标准均方根预测误差最接近1。根据半方差分析的结果，利用ArcGIS 9. 0中的Geostatistic analysis模块采用普通Kriging法进行内插，并根据土壤养分的含量对插值图进行合理的分级，得到各养分含量等级面积的空间变异(—)。
有机质含量范围在施肥前后基本保持一致，90%以上处于20—40 g/kg范围。作物种植前，有机质含量分布较为不均，含量较高的区域主要分布在园地和林地区。作物收获后其含量分布呈现总体有差异、局部有均衡的现象，体现2点特征：(1)作物收获后，有机质含量处于30—40 g/kg范围的区域面积减少，原有林地区的有机质下降较快，主要原因：其一可能是种植区内作物在生长过程中吸收了土壤原有的碳，而新的有机质的形成需要一定时间转化，因此，总体表现为普遍性下降；其二可能是园地和林地较近，林地的有机质含量通过园地的取样点有机质含量体现，而水果在生长过程中消耗土壤中原有的有机质，而短期内又无法及时得到补充。(2)作物收获后的有机质含量以水源地周边居高，并增加较快，原因可能是有机质的形成与水有较密切的关系。一般而言同一地点土壤有机质在短时内不会有大的变化，而两次试验结果却在不同地方有不同表现，这一现象有必要针对不同地区产生的效应验证其碳收支现象()。
作物种植前后碱解氮含量的空间分布格局变化不大，作物种植前碱解氮含量75%处于50—60mg/kg，作物收获后各区域的氮含量明显增加，范围在70—90mg/kg占65%。主要特点：(1)林地和园地的土壤碱解氮的含量较高，这可能与农户长期对经济作物高施肥有较强的相关性。(2)水源地，特别是水库和水源下游地带氮的含量更高。可能过多投入的氮随着水分的流动而富集在水库、坑塘或流向地势较低的区域。(3)作物收获后，土壤碱解氮的含量普遍高于作物种植前，主要原因在于土壤中的氮降解和固化需要一定的时间。这一点从第2年作物种植前的土壤抽样结果得到证明。具体见较作物收获后速效磷的含量较种植前普遍性地有所提高，局部性增加较大，其中城郊结合部和杭坪镇的磷含量更高，超过了100 mg/kg。原因可能有3点：(1)目前农户主要使用的复合肥，导致速效磷含量普遍性有所增加，但受使用量的影响，增加的幅度不大。(2)蔬菜和水果种植区速效磷的浓度增加大。分析原因可能是蔬菜的养分投入较大，除了氮投入超过一般作物的2—3倍外，磷的投入也有明显增加，这一点可从农户养分投入调查中得到证实。从地理位置上看，速效磷含量最高的地方正是城郊蔬菜和高山蔬菜所在地()。
作物收获后速效钾的含量较种植前普遍性提高不多，但局部性增加却有超过150 mg/kg。原因可能有3点：(1)目前农户主要使用的复合肥，导致速效钾含量普遍性有所增加，但钾含量并不高，因而，增加的幅度不大。(2)从全县整体上看，钾的投入没有得到人们的重视，蔬菜区域钾浓度增大的原因与农户蔬菜养分投入量大有直接关系。此外，高山蔬菜多是直接销售到上海，其养分及其主要成份会被检测，县内有专职人员定期对养分的施用进行指导，钾的投入也正逐步引起重视。因此，对于农村大部分地区，特别是经济作物种植较为密切的地区，配备养分指导专家有针对性地进行土壤合理配方是必要的()。
(1)运用Moran′s I指数分析、Pearson相关系数分析，以及kriging插值图等方法对作物种植前和收获后的土壤因子空间变异进行分析，发现3种方法所得出的结果可以相互应证，互为补充。Moran′s I指数分析证实随着区域范围的加大，各因子间自相关程度存在差异。Pearson相关系数分析可检验出各因子间的相关性程度。kriging插值图则能直观反应各因子在不同区域的分布特征。
(2)从养分利用特征看：1)有机质和碱解氮与水的关系较为密切，水源区两因子的含量较周边高，可能是水的作用加速了有机质向氮的转化。2)养分投入后，过多的养分在土壤中降解有一段时间过程，这种过程与区域自然环境条件的关系较密切。3)由于主要是通过复合肥被施入到土壤的，测算结果反映磷、钾之间有极强的相关性，作物收获后两者在土壤中的含量均高于作物种植前，这可能与农户好用固定磷、钾比的复合肥有较密切的关系，此外，磷、钾在土壤中降解较慢也是原因之一。
