全 文 ：　 　 倡 唵国家重点基础研究发展规划项目（２００５CB１２１１０８唱４）和黑龙江省自然科学基金重点项目（ZJN０３唱０２）资助
收稿日期 ：２００５唱０７唱１２ 　 改回日期 ：２００５唱０８唱２４
20 年来东北典型黑土地区土壤肥力质量变化 倡
汪景宽１ 　李双异１ 　张旭东２ 　魏 　丹３ 　迟风琴３
（１畅沈阳农业大学 　沈阳 　 １１０１６１ ；２畅中国科学院沈阳应用生态研究所 　沈阳 　 １１００３４ ；
３畅黑龙江农业科学院 　哈尔滨 　 １５００８６）
摘 　要 　通过大量样品分析和数据收集 ，研究了最近 ２０年来东北典型黑土地区土壤肥力质量主要指标 ——— pH 、有
机质 、速效磷 、速效钾和黏粒的变化情况 ，并在地理信息系统和地统计分析的辅助下 ，研究了该 ５项肥力指标以及
综合肥力指数的时空变异规律 。结果表明 ：２０年中该地区土壤 pH 、有机质和速效钾平均含量明显降低 ，速效磷平
均含量有较大增加 ；速效磷的变异系数变化最大 ，pH变化最小 ；２０世纪 ８０年代该地区土壤肥力综合指数以一 、二
级为主（８０％ 以上） ，但 ２１世纪初土壤肥力质量几乎被二 、三级地所占据（９８％ 以上） 。该地区土壤肥力质量明显降
低可能与长期以来重种轻养有关 。
关键词 　黑土地区 　土壤肥力 　地统计学 　空间变异
Spatial and temporal variability of soil quality in typical black soil area in Northeast China in ２０ years ．WANG Jing唱
Kuan１ ，LI Shuang唱Yi１ ，ZHANG Xu唱Dong２ ，WEI Dan３ ，CHI Feng唱Qin３ （１ ．Shenyang Agricultural University ，Shenyang
１１０１６１ ，China ；２ ．Institute of Applied Ecology ，Chinese Academy of Sciences ，Shenyang １１００３４ ，China ；３ ．Heilongjiang
Academy of Agricultural Sciences ，Harbin １５００８６ ，China） ，CJEA ，２００７ ，１５（１） ：１９ ～ ２４
Abstract 　 Geographical information system （GIS）combined with Geostatistics were applied to reveal the spatial and tem 唱
poral variability of soil fertility quality in the topsoil （０ ～ ２０cm）in the northeastern area of China during last two decades ．
Soil fertility quality indices in this study include pH ，organic matter ，available P ，available K and clay ．The results show that
soil pH ，organic matter and available K significantly decrease during the past ２０ years while the available P increases with
the biggest variation coefficient ．The areas with the first and second grade of soil fertility quality were over ８０％ in １９８０ ，
whereas nearly ９８％ belonged to the second and third grade in ２０００ ．So the soil fertility quality has declined in this area
due to the intensive cultivation in the past ２０ years ．
Key words 　 Black soil area ，Soil fertility ，Geostatistics ，Spatial variability
（Received July １２ ，２００５ ；revised Aug ．２４ ，２００５）
黑土地区土壤肥力高 ，性状好 ，但长期以来东北黑土资源未得到应有的重视 ，使水土流失加剧 、土壤退
化速度惊人 。因此黑土地区土壤质量 ，特别是黑土的肥力质量已越来越引起人们的关注 。自 ２０世纪 ８０年
代第 ２次土壤普查以后 ，科技工作者在我国东北黑土地区开展了许多研究工作［１ ～ ３］ ，已取得许多科研成果 ，
并在生产上取得明显经济效益 ，但一直没有进行较大规模黑土地区土壤肥力质量调查和评价工作 。本研究
以黑龙江省典型黑土区嫩江 、五大连池 、克山 、北安 、海伦为研究区域 ，通过野外调查 、采样分析和资料收集 ，
运用地理信息系统（GIS）和地统计分析相结合的方法 ，对 ２０年以来该地区土壤肥力质量及其各指标的变化
特征进行探讨 ，为该区土壤培肥和土壤管理提供理论依据 。
1 　研究区域概况与研究方法
研究区位于黑龙江省中部地区 ，东经 １２４°４４′１７″ ～ １２７°５４′５４″ ，北纬 ４６°５７′３７″ ～ ５０°５９′５７″ ，总面积
３９１４０km２ 。该区主要以海拔高度 ６００m以下低山 、丘陵及平原为主 。土壤类型主要包括黑土 、黑钙土 、草甸
上 、暗棕壤等 。土壤腐殖质层厚度一般在 ２０ ～ １００cm之间 ，土壤多中性 ，以黄土状母质为主 。
土壤样品的采集与分析 。根据 ２０世纪 ８０年代第 ２次土壤普查和农业肥料定点试验的样点分布情况 ，
２０００年春在研究区域对应点共采集土壤样点 ３１９个（见图 １） 。每个样品都是以 １点为中心 ，呈梅花状向四
第 １５ 儍卷第 １期 中 国 生 态 农 业 学 报 Vol 1畅１５ 　 No 厖畅１
２ ０ ０ ７ 乔年 １ 月 Chinese Journal of Eco唱Agriculture Jan ．，　 ２００７ 後
图 1 　采样点分布图 倡
Fig ．１ 　 The distribution of samples
倡 从北至南分别为嫩江 、五大连池 、克山 、北
安和海伦县 。
周辐射 ，取 ０ ～ ２０cm表层土壤 ，每 ５ ～ ６个点混合均匀 。每个中心样点均采
用 GPS定位 ，并详细记载样点周围的环境和土地利用情况 。土壤样品按
常规处理和测定方法进行基本理化性质和必要的肥力指标测定 。 １９８０年
土壤肥力情况是来自于 ２０世纪 ８０年代第 ２次土壤普查对应地点的土壤
测定结果 。根据该地区土壤资源特点 ，结合其他研究成果和专家意见 ，选
取土壤有机质 、速效磷 、速效钾 、pH值 、黏粒含量 ５个因素作为该研究区的
土壤肥力质量评价指标 。
地统计学方法是研究空间变异性的一种重要方法［４ ，５］ ，是以区域化变
量为核心和理论基础 ，以空间相关和半方差函数为基本工具的一种数学地
质方法 。地统计学在矿产地质中的应用目前已经达到了成熟阶段 ，由于自
身的特点 ，其应用范围已由当初的地质学领域逐渐推广应用到土壤学和生
态学等领域中［６ ～ ８］ 。半方差函数是地统计分析的基本工具 ，假设区域化变
量满足二阶平稳和本征假设 ，用于估计半方差函数 γ（h）的公式为 ：
γ（h） ＝ １２ n（h） 钞
n（ h）
i ＝ １
［Z（Xi ＋ h） － Z（Xi）］２ （１）
　 　式中 ，n（h）为由距离 h分隔的观测点的对数 ，h为滞后距离（Lag ） ，
Z（Xi ＋ h）和 Z（Xi）分别为区域化变量 Z（X ）在位置 Xi ＋ h和 Xi 的实测
值 。以 h为横坐标 ，γ（h）为纵坐标作图 ，即为半方差图 。
2 　结果与分析
2畅1 　土壤肥力质量单项指标的描述性统计分析
在 ２０世纪 ８０年代 ，研究区土壤有机质 、速效磷 、速效钾 、pH和黏粒这 ５种土壤肥力指标中 ，有机质 、速
效磷 、速效钾变异相对较大（表 １） ，其变异系数分别为 ４４畅８７％ 、４３畅４４％ 、４１畅５５％ ；pH和黏粒变异较小 ，其变
异系数分别为 ８畅３５％ 和 １８畅９９％ 。 ２１ 世纪初 ５ 种土壤肥力指标中 ，速效磷含量变异最大 ，其变异系数为
７３畅８３％ ；有机质和速效钾居中 ，变异系数为 ４３畅９０％ 和 ３４畅０１％ ；pH和黏粒的变异最小 ，其变异系数分别为
８畅２０％ 和 １８畅７７％ 。同 １９８０年相比 ，速效磷的变异系数有较大上升 ，说明 ２０ 年来速效磷受外界干扰较显
著 ，尤其是土地实行联产承包后 ，农户或农场的管理水平不同 ，化肥施用不均匀 ，使目前土壤速效磷含量变
化较大 。
表 1 　土壤肥力各指标的描述性统计
Tab畅１ 　 Descriptive statistics of soil fertility indices
年 　 份
Years
土壤肥力指标
Soil fertility indices
最大值
Max ．value
最小值
Min ．value
平均数 倡
Mean value
方差
Variance
标准差
Standard deviation
变异系数／％
Coefficient of variation
２０００ c
pH 　 　 　 　 　 ８ 抖畅３４ ４ 贩畅７９ ５ 後畅９８ 　 ０ 4畅２４ 　 ０ 适畅４９ 　 　 ８ n畅２０ 　 　
　 　 有机质／g·kg － １ 晻２０６ 抖畅１９ ６ 贩畅２５ ６１ 貂畅０８ ７１８ 4畅７８ 　 ２６ 适畅８１ ４３ n畅９０
　 　 速效磷／mg·kg － １ 佑１４７ 抖畅００ １ 贩畅７９ ２１ 貂畅５３ ２５２ 4畅８１ １５ 适畅９０ ７３ n畅８３
　 　 速效钾／mg·kg － １ 佑５８４ 抖畅００ ７４ 贩畅００ １９８ 貂畅４９ ４５５７ 4畅６０ ６７ 适畅５１ ３４ n畅０１
　 　 黏 　 粒／％ ４４ 抖畅９０ ６ 贩畅７４ ２１ 貂畅５９ １６ 4畅４０ ４ 适畅０５ １８ n畅７７
１９８０ c
　 　 　 　 pH ７ 抖畅３０ ４ 贩畅７０ ６ 貂畅３６ ０ 4畅２８ ０ 适畅５３ ８ n畅３５
　 　 有机质／g·kg － １ 晻１９５ 抖畅００ ７ 贩畅３１ ７８ 貂畅５７ １２４３ 4畅２７ ３５ 适畅２６ ４４ n畅８７
　 　 速效磷／mg·kg － １ 佑３０ 抖畅５７ ４ 贩畅００ １３ 貂畅０１ ３１ 4畅９２ ５ 适畅６５ ４３ n畅４４
　 　 速效钾／mg·kg － １ 佑７１２ 抖畅００ １０６ 贩畅１０ ３０２ 貂畅８８ １５８３８ 4畅２２ １２５ 适畅８５ ４１ n畅５５
　 　 黏 　 粒／％ ４６ 抖畅９０ ７ 贩畅９４ ２１ 貂畅６９ １６ 4畅９７ ４ 适畅１２ １８ n畅９９
　 　 倡 样本数为 ３１９ 。
2畅2 　土壤肥力质量单项指标的时间变异
２０００年时 ５种土壤肥力指标 ，pH 、有机质 、速效磷 、速效钾和黏粒的平均值分别为 ５畅９８ 、６１畅０８g／kg 、
２１畅５３mg／kg 、１９８畅４９mg／kg 和 ２１畅５９％ ，１９８０ 年 ５ 种肥力指标含量的平均值分别为 ６畅３６ 、７８畅５７g／kg 、
１３畅０１mg／kg 、３０２畅８８mg／kg和 ２１畅６９％ 。 ２０ 年来土壤 pH平均降低 ０畅３８ ，土壤有机质相对降低 ２２畅２６％ ，速
效钾相对降低 ３４畅４７％ ，速效磷则相对增加 ６５畅４９％ ，土壤黏粒几乎没有变化 。
２０　 中 国 生 态 农 业 学 报 第 １５ 卷
2畅3 　土壤肥力质量单项指标的空间变异特征
半方差函数曲线图表示土壤各种指标变量在距离与方向上不同的所有成对点之间的观测值的空间相
关性 。通过半方差函数及半方差函数曲线图可得 ３个极为重要的参数即变程（Range） 、基台值（Sill） 、块金值
（Nugget） 。当半方差函数 γ（h）随着间隔距离 h的增大 ，从非零值达到 １个相对稳定的常数时 ，该常数称为
基台值（C０ ＋ C） ，当间隔距离 h ＝ ０时 ，γ（０）＝ C０ ，该值称为块金常数或块金方差 。基台值是系统或系统属
性中最大的变异 ，半方差函数 γ（h）达到基台值时的间隔距离 a称为变程 。变程表示在 h ≥ a以后 ，区域化
变量 Z（ x）空间相关性消失 。块金值表示区域化变量在小于抽样尺度时非连续变异 ，由区域化变量的属性
或测量误差决定 。另外 ，块金方差／基台值［C０ ／（C０ ＋ C）］可表示空间变异程度（由随机部分引起的空间变
异占系统总变异的比例） ，若该比值较高 ，说明由随机部分引起的空间变异性程度较大 ；反之 ，则由空间自相
关部分引起的空间变异性程度较大 ；若该比值接近于 １ ，则说明该变量在整个尺度上具有恒定的变异 。从结
构性因素角度看 ，C０ ／（C０ ＋ C）可表示系统变量的空间相关性程度 ，如果比例 ＜ ２５％ ，说明变量具有强烈的
空间相关性 ；若比例在 ２５％ ～ ７５％ 之间 ，变量具有中等的空间相关性 ；比例 ＞ ７５％ 时 ，变量空间相关性很
弱［９］ 。而这些参数都可直接从半方差函数曲线图中得到 。根据各个样点各项指标的数据 ，计算出实际半方
　 　
图 2 　土壤肥力各指标的半方差函数图
Fig ．２ 　 The semivariograms of soil fertility indices
第 １ 媼期 汪景宽等 ：２０ 缮年来东北典型黑土地区土壤肥力质量变化 ２１　
差函数的值 γ（h） ，绘制实际半方差函数曲线 ，然后选择几个（各向同性）模型分别进行拟合 ，从中选出最优
的半方差函数理论模型（图 ２ 、表 ２） 。
图 ２表明 ，各指标均能用指数模型或球状模型进行拟合 ，其中 １９８０年的 pH和速效磷最为理想 ，而 １９８０
年的速效钾只是尝试性拟合 ，主要原因可能在于实际取样间距偏大 ，相关间距内的样品“配对”数目少 。由
表 ２可知 ，２０００年速效磷和速效钾的 C０ ／（C０ ＋ C）分别为 ４６畅６３％ 、４９畅９１％ ，１９８０ 年速效钾的 C０ ／（C０ ＋ C）
为 ３４畅１９％ ，说明这 ２项指标存在中等的空间相关性 ，这种分布特征可能是随机性因素和结构性因素共同作
用的结果 ；而 ２０００年的其他 ２项指标和 １９８０年的 ３项指标 C０ ／（C０ ＋ C）的范围从 １畅５６％ 到 １６畅５７％ ，说明
他们存在强的空间相关性 ，这可能是结构性因素作用的结果 。结构性因素如气候 、母质 、地形 、土壤类型等
自然因素可导致土壤养分强的空间相关性 ，而随机性因素如施肥 、耕作措施 、种植制度等各种人为活动使土
壤养分的空间相关性减弱 ，朝均一化方向发展［１０］ 。
表 2 　 2000年和 1980年土壤肥力指标半方差函数理论模型及相关参数 倡
Tab畅２ 　 Theoretical semivariogram models of soil fertility indices and corresponding parameters in ２０００ and １９８０
项目
Items
年份
Years
理论模型
Theory model
块金值
Nugget
基台值
Sill
块金值∶基台值／％
Nugget／Sill
变程／km
Range
拟合度
　 R２
pH ２０００ 指数模型 ０ M畅０２９６ ０ 览畅２７４２ １０ X畅８０ １２３ 牋畅０ ０ �畅６７３
１９８０ 球状模型 ０ a畅０５９ １ 照畅００６ ５ X畅８６ ７１１ 牋畅０ ０ �畅９５６
有机质 ２０００ 球状模型 ３１９ 苘１９２５ d１６ X畅５７ ７１１ 牋畅０ ０ �畅６８９
１９８０ 球状模型 １７０ 苘１３７９ d１２ X畅３３ ８０ 牋畅６ ０ �畅５０７
速效磷 ２０００ 球状模型 １２６ 吵畅７ ２７１ 畅７ ４６ X畅６３ １１７ 牋畅２ ０ �畅５３２
１９８０ 球状模型 １ 妸畅３ ８３畅６ １ X畅５６ ４６１ 牋畅８ ０ �畅９１６
速效钾 ２０００ 球状模型 ２８１１ 耨５６３２ d４９ X畅９１ ２４４ 牋畅２ ０ �畅７９３
１９８０ 球状模型 ５３００ 耨１５５００ y３４ X畅１９ ８３ 牋畅０ ０ �畅２１３
　 　 倡 经 K唱S 检验 ，５项指标均属正态分布 。
１９８０年各项指标的 C０ ／（C０ ＋ C）均小于 ２０００年 ，并且具有强的空间相关性 ，说明人为活动对土地的影
响程度加大 。人们加大了化肥的施用量 ，消弱了结构性因素造成的强的自相关作用 。变程表明属性因子空
间自相关范围的大小 ，它与观测尺度以及在取样尺度上影响土壤养分的各种生态过程相互作用有关 。在变
程之内 ，变量具有空间自相关性 ，反之则不存在 ，它提供了研究某种属性相似范围的一种测度 。在研究区域
内 ，pH和速效磷 ２０００年相对于 １９８０年的变程有明显减少 ，说明 pH和速效磷分布的均一性减弱 ，在小范围
变异加强 ，整体分布趋向复杂 ；有机质和速效钾则正好相反 ，其含量分布趋向于简单化 。
2畅4 　 20年来土壤肥力质量的时空变化特征
土壤肥力质量综合指数的计算 。为避免人为主观因素的干扰 ，土壤肥力质量单因子权重采用因子分析
法计算 。因子分析是寻找对观察结果起支配作用的潜在因子（潜变量）的探索性统计分析方法 。因子载荷
aij表明第 i个变量在第 j 个主因子上的负荷 ，或称为第 i个变量在第 j 个主因子上的权 ，它反映出第 i个变
量在第 j 个主因子上的相对重要性 。本研究以特征值 ＞ １为选取主因子的条件进行因子分析 ，得到各评价
因子主成分的特征值和贡献率 ，并由因子载荷矩阵计算土壤肥力质量指标的公因子方差及权重值 。经计算
得出 pH 、有机质 、速效磷 、速效钾和黏粒的 ２０００年权重值分别为 ０畅１６６ 、０畅２１７ 、０畅２０３ 、０畅１８５ 、０畅２２９ ，１９８０年
的权重值分别为 ０畅１９１ 、０畅１５４ 、０畅２１８ 、０畅１７８ 、０畅２５９ 。
采用隶属函数确定评价因子的单因子评价指数 。首先建立各评价指标的隶属函数模型 ，将肥力指标的
隶属度函数分为两种类型即 S型隶属函数和抛物线型隶属函数 ，土壤有机质 、速效磷 、速效钾可采用 S型隶
属函数 ，pH值和黏粒可采用抛物线型隶属函数［１０］ 。根据长期以来生产实践的总结以及专家评定 ，根据作物
效应曲线 ，通过确定各个函数的转折点取值而拟合出每个函数的确切函数式 ，进而得到每个评价因子的单
因子评价指数 。在获得每个样点各个指标的单因子指数后 ，根据模糊数学中的加乘法原则 ，土壤肥力质量
综合评价指标值（ IFI）采用式（２）计算 ：
IFI ＝ 钞 （qi × wi） （２）
式中 ，qi 为第 i项土壤肥力评价指标的隶属度值 ，w i 为第 i项土壤肥力评价指标的权重系数 。 IFI 取值为
０ ～ １之间 ，其值越高 ，表明土壤肥力质量越好 。
２２　 中 国 生 态 农 业 学 报 第 １５ 卷
土壤肥力质量的分级 。根据土壤取样点并结合等距划分法 ，参考各土种类型土壤的产量水平 ，提出该
地区土壤肥力质量分级参考标准 ：IFI ≥ ０畅８０ 、IFI为 ０畅８０ ～ ０畅７０ 、０畅７０ ～ ０畅６０ 、０畅６０ ～ ０畅４０和 IFI ＜ ０畅４０分
别为１ ～ ５级土壤 。
经拟合得到本研究区土壤肥力质量综合指数的半方差函数图（图 ３） 。由图 ３可知 ２０００年与 １９８０年的
土壤肥力综合指数具有相似的空间变异 ，随滞后距离 h的增加 ，它们表现出特别有规律的持续增加趋势 ，之
后趋向平稳 ，说明它们存在很好的空间结构性 。表 ３表明 ，２０世纪 ８０年代和 ２１世纪初的 C０ ／（C０ ＋ C）分别
为 １０畅８２％ 和 １２畅４５％ ，均存在强的空间相关性 ，这是由结构性因素所致 。 ２０００年的 C０ ／（C０ ＋ C）较 １９８０年
有所增加 ，人为活动可能是导致这一现象的主要因素 。由图 ４和表 ４可知 ，２１ 世纪初黑龙江 ５个县土壤肥
力质量为二 、三级的土地占 ９８％ 以上 ，而 ２０世纪 ８０年代时土壤质量为一级的土地占 １３畅８０％ 。该地区土壤
肥力质量降低的区域占据了相当大面积 ，一级质量降低 １３畅７２％ ，三级质量土壤大量增加 。说明 ２０年来该
地区土壤肥力出现明显下降 ，对目前黑土地区农业生产构成潜在威胁 。
图 3 　土壤肥力质量综合指数半方差函数图
Fig ．３ 　 The semivariograms of integrated indexes of soil fertility quality
表 3 　土壤肥力质量综合指数半方差函数的理论模型及相关参数 倡
Tab畅３ 　 Best fitted semivariogram models of integrated indexes of soil fertility quality and corresponding parameters
年份
Years
理论模型
Theory model
块金值
Nugget
基台值
Sill
块金值∶基台值／％
Nugget／sill
变程／km
Range
决定系数
R２ 趑
２０００ 貂指数模型 ０ ^畅００２３７ ０ Ζ畅０１９０４ １２ d畅４５ ５８ .畅８ ０ 亖畅６８１
１９８０ 貂球状模型 ０ ^畅０００６１ ０ Ζ畅００５６４ １０ d畅８２ ３２ .畅５ ０ 亖畅１５１
　 　 倡 经 K唱S 检验 ，１９８０ 年 、２０００ 年土壤肥力质量综合指数的分布类型均属正态分布 。
图 4 　 2000年和 1980年黑土肥力质量分级和 20年来变化图
Fig ．４ 　 The grade maps of soil fertility quality in １９８０ and ２０００ and the changing map in past ２０ years
第 １ 媼期 汪景宽等 ：２０ 缮年来东北典型黑土地区土壤肥力质量变化 ２３　
表 4 　 20年来黑土区土壤肥力质量变化情况
Tab畅４ 　 Changes of soil fertility quality in black soil area in past ２０ years
年份
Years
项 　 目
Items
等级Grade
１ 刎２ 膊３ 寣４ ５
２０００ �面积／km２ �　 　 　 ３１ 煙畅８６ 　 　 ２４１９８ ⅱ畅５５ １４２１１ j畅３３ ６９８ 蝌畅２６ ０ "
所占百分数／％ 　 　 　 ０ 煙畅０８ 　 　 ６１ ⅱ畅８３ ３６ j畅３１ １ 蝌畅７８ ０ "
１９８０ �面积／km２ �　 　 　 ５４０１ 煙畅３５ 　 　 ２５９９７ ⅱ畅４５ ７３２５ j畅２６ ４１５ 蝌畅９４ ０ "
所占百分数／％ 　 　 　 １３ 煙畅８０ 　 　 ６６ ⅱ畅４２ １８ j畅７２ １ 蝌畅０６ ０ "
变化值 变化面积／km２ 　　 　 　 － ５３６９ 煙畅４９ 　 　 － １７９８ ⅱ畅９０ ６８８６ j畅０７ ２８２ 蝌畅３２ ０ "
变化百分数／％ 　 　 　 － １３ 煙畅７２ 　 　 － ４ ⅱ畅５９ １７ j畅５９ ０ 蝌畅７２ ０ "
3 　小 　结
东北黑土区土壤肥力单项指标中 ，pH变异系数变化最小 ，由 １９８０年的 ８畅３５％ 降低到 ２０００年 ８畅２０％ ，
除低了 １畅８％ ；速效磷变异系数变化最大 ，由 １９８０年 ４３畅４４％ 增加到 ２０００年的 ７３畅８３％ ，增加了 ７０％ 。 ２０年
来研究区土壤 pH值平均降低 ０畅３８ ，土壤有机质相对降低 ２２畅２６％ ，速效钾相对降低 ３４畅４７％ ，速效磷则相对
增加 ６５畅４９％ ，土壤黏粒几乎没有变化 ，速效磷的这种明显增加与施肥密切相关 。 ２０世纪 ８０年代该地区土
壤肥力综合指数以一 、二级为主（８０％ 以上） ，但 ２１ 世纪初土壤肥力质量几乎被二 、三级地所占据（９８％ 以
上） 。目前该地区土壤肥力质量降低的区域占相当大面积 ，这可能与人们长期以来重种轻养有关 。
参 　考 　文 　献 h
１ 　 王建国 ，刘鸿翔 ，王守宇 ，等 ．黑土农田养分平衡与养分消长规律 ．土壤学报 ，２００３ ，４０（２） ：２４６ ～ ２５１
２ 　 汪景宽 ，赵永存 ，张旭东 ，等 ．海伦县土壤重金属含量的空间变异性研究 ．土壤通报 ，２００３ ，３４（５） ：３９８ ～ ４０３
３ 　 辛 　 刚 ，颜 　 丽 ，汪景宽 ，等 ．不同开垦年限黑土有机质变化的研究 ．土壤通报 ，２００２ ，３３（５） ：３３２ ～ ３３５
４ 　 郭旭东 ，傅伯杰 ，马克明 ，等 ．基于 GIS 和地统计分析的土壤养分空间变异特征研究 ．应用生态学报 ，２０００ ，１１（４） ：５５７ ～ ５６３
５ 　 李 　 艳 ，史 　 舟 ，徐建明 ，等 ．地统计学在土壤科学中的应用及展望 ．水土保持学报 ，２００３ ，１７（１） ：１７８ ～ １８２
６ 　 毛任钊 ．河海低平原盐渍区土壤养分空间变异性研究 ．中国生态农业学报 ，２００１ ，９（４） ：６３ ～ ６６
７ 　 赵玉国 ，张甘霖 ，张 　 华 ，等 ．海南岛土壤质量系统评价与区域特征探析 ．中国生态农业学报 ，２００４ ，１２（３） ：１３ ～ １５
８ 　 肖 　 玉 ，谢高地 ，安 　 凯 ．土壤速效磷含量空间插值方法比较研究 ．中国生态农业学报 ，２００３ ，１１（１） ：５６ ～ ５８
９ 　 後Cambardella C ．A ．，Moorman T ．，Novak J ．M ．，et al ．Field scale variability of soil properties in central Iowa soils ．Soil Sci ．Soc ．Am ．J ．，
１９９４ ，５８ ：１５０１ ～ １５１１
１０ 沣Chien Y ．J ．，Lee D ．Y ．，Guo H ．Y ．，et al ．Geostatistical analysis of soil properties of mid唱west Taiwan soils ．Soil Science ，１９９７ ，６２（４） ：２９１ ～
２９８
２４　 中 国 生 态 农 业 学 报 第 １５ 卷