重庆市丰都县

重庆市丰都县紫色土养分空间变异及土壤肥力评价

陶睿;王子芳;高明;孙宇

【摘要】综合运用地统计学、GIS和模糊数学相结合的方法,分析了丰都县紫色土

养分的空间变异规律,对其土壤肥力进行综合评价,并绘制了养分空间分布图和土壤

综合肥力指数空间分布图,为该县紫色土养分分区管理及精准施肥决策等提供理论

依据.结果表明:研究区内5种肥力指标中,除有效磷属于强变异性外,pH、有机质、

碱解氮和速效钾均具有中等变异性.土壤pH、速效钾的最适模型是指数模型,块金

值与基台值之比分别为17.07%和10.65%,有较强的空间相关性;用线性模型可较好

地模拟有机质、碱解氮和有效磷,块金值与基台值之比分别为76.37%、84.21%和

87.92%,空间相关性较弱.土壤pH、有机质、碱解氮、有效磷和速效钾的变程分别

为5.43、28.02、31.24、36.48和2.28km.土壤pH在长江以北较高与较低区域

都呈带状并相间分布,在长江以南呈西低东高状态;有机质、碱解氮、有效磷及速效

钾的分布存在明显的方向不均匀性,都呈零星斑状.丰都县土壤肥力质量呈现北低南

高的趋势,5个等级地块分别占8.14%、18.86%、31.75%、28.58%、12.67%.丰

都县土壤肥力总体水平中等偏上,这主要与该地区地形地质有关,同时也受到社会经

济与人为管理措施的影响.%Thespatialvariabilityofsoilnutrientsofpurple

soilinFengduCityofChongqingwascomprehensivelyanalyzedby

geostatistics,r

diagramsofthespatialdistributionofsoilnutrientswerethen

quantitativelyconstructedonordertopotentiallyprovidestheoretical

basnotonlyfortheuofmanagementzones(MZs)withpurplesoilin

Fengdu,ultsshowedthatmoderatelevelofspatialvariabilityoccurredin4nutrients

contentsintherearchregionexceptavailableP(AP)whichbelongedto

get/stillratioswere17.07%forpHand10.65%for

availableK(AK),respectively,whichshowedstrongspatialcorrelation,and

henugget/still

ratioswere76.37%,84.21%and87.92%fororganicmatter(SOM),

availableN(AN)andavailableP(AP),respectively,whichshowedweak

spatialcorrelation,andcouldbewellmodeledwiththesphericalmodel.

TherangesofpH,SOM,AN,APandAKw学习报告怎么写 ere5.43,28.02,31.24,36.48and

2.28km,edbandingdistributionbothinthehigher

andlowerregionsinthenorthoftheYangtzeRiver,butshowedhighin

thewestregionandlowinthesouthregioninthesouthoftheYangtze

River;allthedistributionofSOM,AN,APandAKshowedapparent

lfertility

qualityofFengdudecreadfromsouthtonorth,butgenerallywasabove

themiddlelevel,theareasofthefivegradesoffertilityqualitywere8.14%,

18.86%,31.75%,28.58%and12.67%ofthetotalarea,

fertilityqualityofFengduwasaffectedbysoiltexture,socio-economicand

artificialcultivationmanagement.

【期刊名称】《土壤》

【年(卷),期】2017(049)001

【总页数】7页(P155-161)

【关键词】紫色土;养分;地统计学;空间变异;肥力

【作者】陶睿;王子芳;高明;孙宇

【作者单位】西南大学资源环境学院,重庆400715;西南大学资源环境学院,重庆

400715;西南大学资源环境学院,重庆400715;黔西南州农业委员会,贵州位于英文 兴义

562400

【正文语种】中文

【中图分类】S159

土壤资源作为一种脆弱性的非再生资源[1]，具有高度变异性的时空连续体，受气

候、母质、植被覆盖、地形等自然因素、人为作用的影响。土壤养分是土壤的基本

属性和本质特征，也是衡量土壤肥力的重要指标。土壤肥力是土地生产力的基础，

土壤肥力水平的高低直接关系到作物生长和农业生产的结构、布局及效益等方面。

对土壤养分空间变异特征的合理分析和土壤肥力的客观评价，可以有效地揭示土壤

养分的空间分布情况与各区域土壤肥力状况，为作物合理布局、土壤精准施肥、养

分管理及防治土壤退化等提供科学依据，对推动土壤科学定量化研究与精准农业的

实施也具有重要意义[2]。

地统计学是在传统统计学基础上发展起来的空间分析方法，不仅能有效地揭示属性

变量在空间上的分布、变异和相关特征，还能将空间格局与生态过程联系起来，可

有效地解释空间格局对生态过程与功能的影响[3–5]，地统计学是分析土壤特性空

间分布特征及其变异规律最为有效的方法之一[6]。GIS可以将系统变量的属性数

据同地理数据相结合，使大区域范围内进行地统计分析变得较为方便。二者的有机

结合可以取长补短，充分发挥各自的优势。近年来，在对土壤性质的空间变异特征

研究的方法上利用地统计学和GIS技术相结合已成为土壤科学的研究热点之一[7–

8]。大量学者对农田、植烟区、城市郊区等的土壤养分空间变异特征及其土壤肥力

评价作了相关研究[9–14]，但有关三峡库区紫色土的相关研究仍不全面。

紫色土主要分布于我国亚热带地区，重庆市是紫色土集中分布区之一，重庆人多地

少，土壤养分失衡、坡耕地面积大、水土流失及侵蚀等问题日益突出，特别是丰都

县地处三峡库区腹心位置。因此本研究以重庆市丰都县紫色土为对象，采用地统计

学、GIS和模糊数学相结合的方法研究该区紫色土养分的空间变异特征及肥力状况，

揭示其空间变异规律，以期为该地区农业生产合理施肥和土壤肥力的培育提供科学

依据。

1.1研究区概况

丰都县位于重庆中部，介于10728′03″~10812′37″E，2933′18″~

3016′25″N。全县地形呈南高北低、“四山”夹“三槽”的地形，县内最低海拔

175m，最高海拔2000m。丰都县地处三峡库区腹心地带，长江自西南向东北

从县境中部横穿而过，境内长度为47km，长江以南多为深丘，中低山地形地貌，

长江以北多为浅丘陵地形地貌。属亚热带润湿季风气候区，气候温和，四季分明，

年均气温18.5℃，年均降雨量823.8mm。土壤类型以紫色土为主，主要分布在

长江以北、长江沿岸及中部以东地区。丰都县地质结构较为复杂，地貌形态多样，

各种物产、资源较为丰富。

1.2样品采集与分析

土壤采样点遵循代表性、均匀性的原则，结合土地利用现状、地貌单元及植被，根

据丰都县紫色土分布情况进行布置。研究区共布置3698个取样点(图1)，其中，

山地地形设置采样点1650个，区域主要种植玉米、水稻、马铃薯、烤烟和豆类，

丘陵、平坝地形各设置采样点1872个、176个，种植玉米、水稻。因丰都县南

部的南天湖镇、三抚林场、暨龙镇、都督乡、太平坝乡及七跃山林场的土壤类型主

要是黄壤，所购买协议书范本 以采样区域及本文的相关分析均不涉及以上乡镇及林场。利用手持式

GPS定位，记录采样点的经纬度和高程，以定位点为中心，在半径10m的圆形

区域内采集4点表层(0~20cm)土样混合，四分法取大约1kg土样带回室内分

析。

样品带回实验室后经风干，研磨过筛备用。采用玻璃电极法测定土壤pH；采用重

铬酸钾氧化法测定土壤有机质含量；采用碱解扩散法测定有效氮含量；采用钼锑抗

比色法测定有效磷含量；采用火焰光度法测定速效钾含量[15]。

1.3数据处理

采用域法识别特异值，即平均值3倍标准差去除异常值，在此区间外的数据定位

特异值，然后分别用正常最大值或最小值代替特异值。利用SPSS18.0软件对样品

数据进行描述性分析，并对数据进行Kolmogorov-Smimov(K-S)检验。采用地

统计学软件GS+7.0对数据进行半方差函数的计算和理论模型拟合。在

ArcGIS9.3平台上进行Kriging插值和图形编辑，制作土壤养分空间分布图和土壤

综合肥力指数空间分布图。

1.3.1评价指标的选取及权重的确定因土壤肥力形成机制较为复杂，不同学者对土

壤肥力的理解各不相同，其评价方法与评价指标也不尽一致。作为土壤肥力评价，

应以土壤的养分含量为主[16]。本研究按照指标对土壤生产力的影响，选取pH、

有机质、碱解氮、有效磷、速效钾5个指标，建立三峡库区紫色土肥力质量评价

指标体系。各参评指标的权重(表1)采用层次分析法[17]确定，经随机一致性检验

指标CR=0.002<0.10，通过一致性检验，判断矩阵无需调整。

1.3.2评价指标的隶属函数构建根据土壤肥力指标与作物生长效应曲线将参评隶属

度函数分为两种类型，即S型和抛物线型[18]，并将曲线型函湖水简笔画 数转化为相应的折线

型函数。将每个样本的原始统计资料分别代入相应的隶属函数公式，可以把不规则

分布的、有单位的、定量或定性描述的原始数值转化为0.1~1.0分布的、无量纲

差异的隶属度值。

S型隶属度函数：属于这类函数的土壤肥力指标有土壤有机质、碱解氮、有效磷、

速效钾。相应的隶属度函数表达式为：

抛物线型：属于这类函数的土壤肥力指标有土壤pH。相应的隶属度函数表达式为：

本文根据重庆紫色土肥力特性、作物品种特点以及相关研究的结果[9,12,19]，确

定各指标的拐点值(表2)。

1.3.3综合评价指数的计算以模糊数学中权重加权求和来计算土壤综合肥力指数

(integratedfertilityindex,IFI)，其表达式如下：

式中：n为参评指标数量，fi为第i个指标的隶属度，wi为第i个指标的权重。由

于隶属度取值范围在[0.1，1.0]，因此，IFI取值也处于[0.1，1.0]，该值越接近于

1，土壤肥力越高。

2.1丰都县紫色土养分的描述性统计分析

土壤养分含量数据及K-S检验结果见表3。从表中可以看出，研究区土壤pH、有

机质、碱解氮、有效磷和速效钾的平均含量分别为6.39、13.43g/kg、75.69

mg/kg、22.06mg/kg、96.27mg/kg，其中有效磷含量较高，其他养分含量适

宜。5种土壤养分变异系数在19.71%~106.76%之间，其中土壤pH的变异系

数最小，为19.71%；有机质、碱解氮、速效钾的变异系数比较接近，在37.35%

~51.92%之间，均属于中等变异强度；有效磷的变异系数最大，达到106.76%，

属于高变异强度，其原因除了与该地区不同土地利用类型、土壤母质等有关外，更

重要与施磷肥状况及磷在土壤中的化学行为有关，当季磷肥的利用率低和土壤中的

磷不易迁移等使较多的磷残留在土壤中，最终使得土壤中磷分布不均匀[19]。各土

壤养分含量最大值和最小值差异明显，这表明各地区间土壤养分含量存在本底差异，

盲目地平均施肥，将会造成低养分地区养分继续不足和高养分地区养分过剩[20]。

经过K-S非参数检验结果显示，土壤pH近似服从正态分布，其他养分经对数转

换后均较好地服从对数正态分布，满足半方差函数分析要求。

2.2丰都县紫色土养分的空间结构分析

由土壤养分半方差函数理论模型及其拟合参数(表4)分析可以得到，除速效钾的拟

合度只有0.546外，其余养分半方差函数模型的拟合度均在0.65以上，可见半方

差函数模型的选取符合基本要求。研究区紫色土有机质、碱解氮、有效磷的最适模

型是线性模型，其中碱解氮拟合效果最好，拟合系数为0.89，而pH、速效钾用指

数模型拟合效果较好，拟合系数分别为0.757和0.546。

表4中块金值(C0)表示取样误差和小于取样尺度下的空间变异。基台值(C0+C)表

示变量在研究范围内总的空间变异强度。块金值与基台值的比值C0/(C0+C)表示

随机部分引起的空间异质性占系统灵魂出窍日遍全校 总变异的比例，表明系统变量的空间相关性的程

度[21]。有研究表明，气候、母质、地形、土壤类型等结构性因素可以导致土壤养

分强的空间相关性；而施肥、耕作方式、种植制度等人为活动使得土壤养分的空间

相关性减弱，朝均一化方向发展[22]。从表4可以看出研究区紫色土各养分的块金

值均为正值，说明存在着由采样误差、短距离的变异、随机和结构变异引起的各种

正块金效应。速效钾和pH的块金值与基台值之比均小于25%，为10.65%和

17.07%，说明要素具有强烈的空间相关性，其变异主要受结构性因素影响；有机

质、碱解氮、有效磷的块金值与基台值之比都大于75%，为76.37%、84.21%、

87.92%，要素表现出较弱的空间相关性，说明其随机性因素引起的空间异质性占

主导地位，即土壤的空间变异性主要与土地利用类型、种植制度、管理措施、污染

等人为活动有直接关系[23]。

变程反映空间最大相关距离，它表明土壤养分变量空间自相关范围的大小[24]。在

研究区域内，5种土壤养分的变程在2.28~36.48km内变化，其中有效磷、碱

解氮、有机质的变程较大，分皮开头的成语 别为36.48、31.24、28.02km，pH和速效钾的变

程仅为5.43km和2.28km。说明它们空间自相关范围具有明显的差异，除取样

间距的影响外，还可能是由于研究区生态条件、种植作物和农事操作的不同[7]。

2.3丰都县紫色土养分的空间分布特性

为了能够更深刻、全面和直观地反映土壤养分在空间上的分布特征，在GS+7.0

建立半方差函数的基础上，应用ArcGIS对土壤养分进行Kriging插值，得到各养

分含量的空间分布图(图2)。

由图2可知，在长江以北pH较高与较低区域都呈带状并相间分布；长江以南pH

西低东高，除三建军训个人总结 乡、江池镇大于7.2和栗子镇、武平镇小于5.3以外，其他区域

pH差异不明显。有机质、碱解氮、有效磷及速效钾的分布存在明显的方向不均匀

性，都呈零星斑状。有机质有多个低值中心，相邻区域含量差异不大，最高区域主

要在包鸾镇和武平镇；碱解氮较高区域分布在研究区东南角，较低区域主要在西北

角，其空间分布与有机质基本一致，说明主要受该区土壤质地影响，这与前人的研

究一致[25]。有效磷变化较强烈，斑块小而多，相邻区域含量差异明显，规律性不

强，有多个高值中心，有效磷含量低于13.4mg/kg的区域较连续，主要在北部

及长江以南中部。速效钾的含量分布在长江以北平均低于长江以南，最高区域主要

在武平镇和三元镇，较低区域主要是零星斑状散落在长江以北和研究区中部。

2.4土壤肥力综合评价

2.4.1土壤综合肥力指数(IFI)描述性统计对研究区3698个采样点的IFI值进行统

计分析表明，样点的IFI值近似服从正态分布，处于0.115~0.893之间，平均为

0.457，变异系数为32.31%，属于中等强度变异。

2.4.2土壤综合肥力指数(IFI)空间变异特征及空间分布状况利用GS+7.0对研究区

紫色土综合肥力指数进行半变异函数的计算和理论模型拟合(表5)，可知研究区综

合肥力指数的最适模型是指数模型，并且拟合效果较好，决定系数达到0.916，残

差仅为3.6910–5，块金值与基台值的比值为86.146%，表现出较弱的空间相关

性，说明其随机部分引起的空间异质性占主导地位。

对研究区紫色土综合肥力指数(IFI)进行分级，分为优(>0.8)、良好(0.7~0.8)、一

般(0.6~0.7)、中等(0.4~0.6)和差(<0.4)5个等级，采用地统计的Kriging插值

方法绘制了丰都县紫色土综合肥力指数空间分布图(图3)。从全局看研究区北部到

南部土壤肥力有逐渐升高的趋势，最高区域在西南角和东南角。利用GIS中空间

分析功能进一步对研究区各等级区域的面积进行统计，结果表明，研究区紫色土综

合肥力指数为优的区域较缺乏，仅占研究区域的8.14%，主要分布在武平镇与包

鸾镇，这两个区域土壤肥力较高主要是与该区域以浅丘低山槽坝地形为主，气候湿

润，适宜农作物生方框对勾 长，盛产水稻、玉米、茶叶、蔬菜等经济作物，还与人为管理有

关；为良的区域占18.86%，主要分布在研究区中部丰都县主城区附近，长江河谷

地带，交通便捷，土地集约化程度高，土壤养分受到社会经济因素的影响剧烈，因

此肥力整体水平较高；一般的区域占31.75%，在除崇兴乡以外的乡镇均有分布；

中等的区域占28.58%，占据了三建乡、栗子镇、龙河镇、江池镇、三合镇和保合

镇等乡镇的大部分区域，属低山丘陵区，土壤质地偏砂，容易漏肥，区域内多以林

粮牧相结合方式发展；为差的区域占12.67%，主要在研究区西北角及龙孔乡，区

域内地质结构较为复杂，地貌形态多样，养分含量较低，可因地制宜发展林木业。

1)研究区pH、有机质、碱解氮、有效磷和速效钾的平均含量分别为6.39、13.43

g/kg、75.69mg/kg、22.06mg/kg、96.27mg/kg。除有效磷属于强变异性外，

其他均属于中等变异性，其中土壤pH的变异系数最小(19.71%)，有效磷的变异

系数最大(106.76%)。

2)土壤pH、有机质、碱解氮、有效磷和速效钾的变程分别为5.43、28.02、

31.24、36.48和2.28km。土壤pH和速效钾的最适模型是指数模型，有机质、

碱解氮、有效磷的最适模型是线性模型。其中速效钾和pH具有较强的空间相关性，

有机质、碱解氮、有效磷表现出较弱的空间相关性。

3)研究区土壤pH在长江以北较高与较低区域都呈带状并相间分布，pH在长江以

南呈西低东高状态；有机质、碱解氮、有效磷及速效钾的分布存在明显的方向不均匀性，都呈零星斑状，有机质有多个低值中心，相邻区域差异不大；碱解氮含量随

着地势增高而略有增加；有效磷变化较强烈，斑块小而多，相邻区域差异明显；速

效钾含量高值中心主要分布在长江以南。

4)研究区紫色土肥力质量呈现出明显的北低南高的趋势，5个等级地块分别占

8.14%、18.86%、31.75%、28.58%、12.67%，土壤肥力总体水平中等偏上，这

主要与该地区地形地貌因素有关，同时也受到社会经济与人为管理措施的影响。

【相关文献】

[1]赵其国,孙波,张桃林.土壤质量与持续环境I.土壤质量的定义及评价方法[J].土壤,1997,29(3):

113–120

[2]张玉铭,毛任钊,胡春胜,等.华北太行山前平原农田土壤养分的空间变异性研究[J].应用生态学

报,2004,15(11):2049–2054

[3]李建辉,李晓秀,张汪寿,等.基于地统计学的北运河下游土壤养分空间分布[J].地理科学,2011,

31(8):1001–1006

[4]马渝欣,李徐生,李德成,等.江淮丘陵区农田表层土壤有机碳空间变异——以定远县为例[J].土

壤,2014,46(4):638–643

[5]张文敏,姜小三，吴明,等.杭州湾南岸土壤有机碳空间异质性研究[J].土壤学报,2014,51(5):

1087–1095

[6]江厚龙,王新中,刘国顺,等.烟田土壤质地的空间变异性研究[J].中国生态农业学报,2010,18(4):

724–729

[7]刘国顺,常摇栋,叶协锋,等.基于GIS的缓坡烟田土壤养分空间变异研究[J].生态学报,2013,

33(8):2586–2595

[8]马渝欣,李徐生,李德成,等.皖北平原蒙城县农田土壤有机碳空间变异[J].土壤学报,2014,51(5):

1153–1159[9]崔潇潇,高原,吕贻忠.北京市大兴区土壤肥力的空间变异[J].农业工程学报,2010,

26(9):327–333

[10]于洋,赵叶婷,常庆瑞.渭北台塬区耕地土壤速效养分时空变异特征[J].土壤学报,2015,52(5):

1251–1261

[11]王军,傅伯杰,邱扬,等.黄土高原小流域土壤养分的空间异质性[J].生态学报,2002,22(8):

1173–1178

[12]叶回英语对话练习 春,张世文,黄元仿,等.北京延庆盆地农田表层土壤肥力评价及其空间变异[J].中国农业科

学,2013,46(15):3151–3160

[13]于洋,赵业婷,常庆瑞.渭北台塬区耕地土壤速效养分时空变异特征[J].土壤学报,2015,52(6):

1251–1261

[14]宋效东,郭卢,赵安,等.皖北亳州烟区土壤肥力定量评价[J].土壤通报,2016,47(1):8–14

[15]鲁如坤.土壤农业化学分析方法[M].北京:中国农业科技出版社,2000

[16]周勇,张海涛,汪善勤,等.江汉平原后湖地区土壤肥力综合评价方法及其应用[J].水土保持学报,

2001,15(4):70–74

[17]曹茂林.层次分析法确定评价指标权重及Excel计算[J].江苏科技信息:信息技术,2012(2):39–40

[18]孙波,张桃林,赵其国.我国东南丘陵山区土壤肥力的综合评价[J].土壤学报,1995,32(4):362–368

[19]崔潇潇,高原,吕贻忠.北京市大兴区土壤肥力的空间变异[J].农业工程学报,2010,26(9):327–333

[20]庞夙,李廷轩,王永东,等.土壤速效氮、磷、钾含量空间变异特征及其影响因子[J].植物营养与

肥料学报,2009,15(1):114–120

[21]CambardellaCA,MoormanTB,ParkinTB,-scalevariabilityofsoil

propertiesincentral灶公 Iowasoils[J].SoilScienceSocietyofAmericaJournal,1994,58:1501–1511

[22]王子芳,高明,魏朝富,等.植烟土壤养分的空间变异特征及适宜性评价——以重庆市彭水县为

例[J].西南大学学报(自然科学版),2008,30(1):98–103

[23]胡江玲,张高.新疆精河流域土壤质量空间分异规律分析[J].干旱区资源与环境,2011,25(6):

113–121

[24]王政权.地统计学及在生态学中的应用[M].北京:科学出版社,1999:65–132

[25]王淑英,路苹,王建立,等．不同尺度下土壤有机质和全氮的空间变异特征:以北京市平谷区为

例[J]．生态学报，2008,28(10):4957–4964