我国耕地资源匮乏、整体质量不高^[1]，因长期过度施用化肥、农药造成的土壤污染等问题致使耕地不健康^[2]。耕地健康的核心在于土壤健康，土壤是农作物生产之根本，耕地土壤健康关系到粮食安全问题，对国民经济发展、社会稳定和国家自立的全局性重大战略具有重要意义^[3]。国际热点研究显示，土壤健康管理有可能是新一轮农业革命的核心^[4]。联合国粮农组织将土壤健康定义为：土壤作为一个生命系统具有的维持其功能的能力，健康的土壤能维持土壤生物群落多样性、抵抗病虫害的传播、促进基本的植物养分循环、改善土壤结构、提高作物产量，以及改善环境^[5]。曹志洪^[6]认为土壤健康质量是土壤净化容纳污染物质、维护和保障人类及动植物健康能力的量度。因此，耕地土壤健康是指在自然或管理的生态系统边界内，土壤净化容纳污染物质，维持水、气质量以及促进人类和动植物健康的能力。当前我国耕地基础地力对粮食生产的贡献仅为50%左右^[3]，在此形势背景下，开展耕地土壤健康评价，研究耕地土壤健康特征，为构建耕地健康、实现耕地数量-质量-生态三位一体保护的目标提供参考。

“健康”概念最初应用于医学领域, 随着社会与科学的发展而后逐渐被用于其他各领域。在20世纪40年代Leopold^[7]首先提出了土地健康概念，1943年新西兰在土壤与健康期刊提出了“健康土壤-健康食品-健康人群”的理念^[8]。随后国内外不同学者逐渐开展与土壤健康相关研究。20世纪80年代中期，加拿大农业部开展“土壤健康”研究项目^[9]。澳大利亚学者Pankhurst等^[10]认为土壤健康指土壤具有连续的生物生产力、促进动植物健康水平以及维持环境质量。国内单美^[11]、李强^[12]、杨希越等^[13]从耕地质量、产能与土壤环境相结合的角度对耕地健康进行评价。谭晶今^[14]选取了自然环境、社会环境以及经济环境3个方面的指标进行了研究区耕地健康评价。上述学者从不同角度构建了不同的健康评价的指标体系，评价指标中均采用了土壤重金属指标，但并未涉及有益微量元素指标。目前，基于生态环境视角的耕地评价研究较少^[11]，本研究从耕地土壤环境的角度选取评价指标体系。耕地土壤环境指标包括了土壤重金属含量、有益微量元素含量、灌溉水污染等，耕地土壤中重金属、灌溉水污染物等含量越高，耕地土壤环境受到污染越严重，耕地土壤健康状况越差，反之耕地土壤越健康；此外，在一定范围内，耕地土壤中有益微量元素含量越高，耕地土壤越健康^[11]。本文在结合已有研究的基础上，从正反两面考虑耕地土壤健康的影响因素，进行稻田土壤健康评价，探索较为合理的耕地土壤健康评价模式，为耕地土壤的保护与提升提供参考。

1 材料与方法 1.1 研究区域概况

奉新县位于江西省西北部，距省会南昌65 km，地处114°44′ ~ 115°33′E、28°34′ ~ 28°52′N之间，为东西长、南北窄的长条形地域。该区属中亚热带气候，四季分明，降雨充沛，地形以丘陵山地为主，自西北向东南逐渐倾斜，形成西高东低的地势，大小丘陵岗地分布其间。2015年全县土地总面积为160 306.86hm²，其中耕地为39 355.38 hm²，占全县总面积的24.55%。

1.2 数据采集与检测方法

本研究的稻田地块数据来源于奉新县2015年耕地质量等别更新数据，从中提取奉新县稻田图斑作为评价单元，共计16 819个，面积共计33 665.20 hm²。灌溉水环境质量数据来源于土地质量地球化学评价数据、奉新县灌溉规划，选取GB 5084数据来源于《农田灌溉水质标准》作为评价的标准值。土壤重金属数据来源于74个土壤实际采样数据样点，参照《土壤环境质量标准》作为污染评价标准值。有益微量元素基于179个土壤采样点，其中有效钼数据来自奉新县测土配方施肥数据。奉新县水稻土主要包括潴育型、潜育型、淹育型3个亚类，潴育型主要为麻沙泥田、潮沙泥田、红泥田、红砂泥田、黄泥田，潜育型、淹育型以麻沙泥田等为主。样点采集过程中充分考虑了土壤类型的分布，各土属的样点数比例与其面积比例一致，并考虑样点在空间内均匀分布的原则，各乡镇内均有代表样点。随后对样点进行普通克里格插值，土壤重金属60个样点进行插值，14个样点进行验证，检验精度R²为0.45 ~ 0.58，标准化均方根为0.72 ~ 0.86；有益微量元素150个样点进行插值，29个样点进行验证，检验精度R²为0.55 ~ 0.62，标准化均方根为0.65 ~ 0.78。克里格方法(Kriging)能够根据已知样点的属性和空间位置关系对没有采样点的局部区域进行无偏最优估值^[15]，既考虑了土壤样点间的空间结构关系，又加入了本身属性作为参考^[16]，在保证预测精度的同时也节约了采样成本。有益微量元素、有效钼与土壤重金属采样点分布见图 1。

土壤样品在室温下自然风干，磨碎，过100目尼龙筛，取适量样品进行测定。土壤重金属含量参照HJ/T 166-2004《土壤环境监测技术规范》^[17]测定，土壤As、Hg采用PERSEE原子荧光光度计(PF-7型)测定；Cd、Pb和Cr采用火焰原子吸收光谱仪(Agilent 200AA型)测定。有益微量元素测定参照《土壤农业化学分析方法》^[18]，土壤有效Zn、Cu、Mn、Fe、Mo采用DTPA浸提-原子吸收分光光度法测定，土壤有效B采用沸水浸提-姜黄素比色法测定。灌溉水环境质量测定: pH采用玻璃电极法(GB/ T6920)，六价铬采用二苯碳酰二肼分光光度法(GB/T7467)，总Hg采用冷原子吸收分光光度法(GB/T7468)，总As采用二乙基二硫代氨基甲酸银分光光度法(GB/ T7485)。

1.3 稻田土壤健康综合评价方法 1.3.1 指标体系的构建

伴随工业化的发展引起的土壤重金属污染问题日益突出^[19]，一旦其在土壤中的含量超过限定标准，并通过植物的富集作用与食物链的传播，将对人体造成严重的危害^[20]，土壤重金属指标的选取也是耕地土壤健康评价区别于传统耕地质量、耕地地力评价的一个重要方面。土壤重金属中主要对环境有危害的“五毒”元素为Hg、Cd、Pb、Cr和类金属As, 其中以Hg毒性最大, 其次为Cd、Pb、Cr、As^[21]。灌溉水环境质量对于无公害健康农产品的种植生产基地也很重要^[22]。此外，有益微量元素的滋养能使作物产品品质更加绿色、有机、无公害^[11]。

根据奉新县实际情况，本研究稻田土壤健康评价从灌溉水环境质量、土壤重金属与有益微量元素3个方面构建评价指标体系(表 1)。

1.3.2 稻田土壤健康综合评价方法

基于地力水平与健康状况关系的分析，在一定的浓度范围内，土壤健康状况的改变并不会对作物产量有影响，但是会降低或者提高作物品质。因此，稻田土壤健康状况系数的确定采用“健康状况系数的累加模式”，即以未受到污染和未含有有益微量元素的耕地为基准累加，依据灌溉水环境、土壤重金属、有益微量元素的对应分级标准进行加和。

将稻田土壤健康系数依据一定的分级标准进行健康等级划分，分级标准见表 2。

1.4 单一指标评价方法 1.4.1 内梅罗综合污染指数法

内梅罗综合污染指数法在多因子加权过程中能够兼顾极值，全面反映不同污染物的污染水平，突出主要污染物对环境质量造成的危害^[23-24]。本研究中灌溉水环境质量、土壤重金属两个指标值均采取该方法计算。计算公式如下：

$ {P_{综}} = \sqrt {\left[ {{{\left( {{{\bar P}_i}} \right)}^2} + {{\left( {{P_{i, \max }}} \right)}^2}} \right]/2} $

(1)

$ {P_i} = {C_i}/{S_i} $

(2)

式中：P_综为内梅罗综合污染指数，${\bar P_i}$为各项污染物的污染指数平均值，P_{i, max}为各污染物的污染指数最大值，P_i为污染物i的污染指数，C_i为污染物i的实测含量，S_i为污染物评价标准值。得出的综合指数依据分级标准进行分级、赋分，灌溉水环境质量、土壤重金属分级标准见表 3、表 4。

1.4.2 最小因子法

最小因子法运用木桶原理，将木桶的最大容量作为土壤有益微量元素的得分高低，最低分数的因子限制着土壤的最终得分，采取最小值定律，选取所有分值中的最小值作为评价单元的土壤有益微量元素的最终得分^[25-26]。有益微量元素分值依据峰型隶属度函数进行计算，计算公式如下：

$ {y_i} = \left\{ \begin{array}{l} 0\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{u_i} \ll {u_{t1}}或{u_i} \ge {u_{t2}}\\ 1/(1 + {a_i}{({u_i} - {c_i})^2})\;\;\;\;\;\;{u_{t1}} < {u_i} < {u_{t2}}\\ 1\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{u_i} = {C_i} \end{array} \right. $

(3)

式中：y_i为评价因子的隶属度，u_i为评价因子的值，a_i为常数，u_t1、u_t2分别为指标上、下限值。

有益微量元素隶属度值见表 5。

2 结果与分析 2.1 单一指标评价结果 2.1.1 灌溉水环境质量

依据GB5084—2005 《农田灌溉水质标准》，将水稻作物农田灌溉用水水质标准作为基础值，进行内梅罗综合污染指数计算，结合奉新县灌溉分区图，运用Arc GIS软件制作出灌溉水环境质量级别图(图 2)。

奉新县灌溉水环境质量总体良好，95.59%稻田灌溉水环境质量为1级，水质清洁；3.56%为2级，水质尚清洁，分布在澡下镇、赤岸镇、赤田镇、冯川镇等；0.85%为3级，水质受到一定污染，分布在冯川镇、赤岸镇等。

流域内土地利用情况与社会经济发展状况对该流域的水质会产生重要的影响^[27-28]。此外，汽车尾气的排放也会引起一定的污染。奉新县县城地处冯川镇，县城经济发展更快，工业化程度相对偏远乡镇较高，汽车车流量也相对较高，工业化产生的废弃物与汽车尾气通过干湿沉降，最终经河流灌溉流入农田，对县城及周边稻田耕地灌溉水质产生一定的不良影响，降低了灌溉水环境质量。

2.1.2 土壤重金属

对土壤重金属样点数据进行探索性分析，得出5种土壤重金属元素各统计值(表 6)。

根据统计结果，可以发现，Cr元素含量峰度小于3，表明其数据分布与正态分布相比较为平坦，含量分布相对均匀；Pb、Cd、As、Hg 4种元素含量峰度均大于3，表明其数据分布相对集中，较正态分布而言较为陡峭。Cr、Pb、Cd、As和Hg 5种元素的偏态均大于0，说明其数据均向右偏态分布，表明比平均值小的重金属含量数据相对较多。

对74个样点利用ArcGIS软件进行普通克里格插值计算，将属性赋到研究区所有评价单元上，随后运用内梅罗综合污染指数法进行计算，根据土壤重金属分级标准，运用ArcGIS软件制作出奉新县5种土壤重金属含量分布图与土壤重金属级别图(图 3、图 4)。

奉新县99.31%稻田土壤重金属级别为1级，只有0.69%为2级分布在西部的石溪乡等，该评价结果表明，奉新县稻田土壤重金属污染很少，大部分稻田没有受到重金属污染。

2.1.3 有益微量元素

对有益微量元素采样结果进行探索性分析，得出6种元素各统计值见表 7。

上述统计结果表明，有效Mn、有效Mo的峰度均小于3，其数据分布较为均匀，与正态分布相比较为平缓；有效Fe、有效Cu、有效B、有效Zn的峰度均大于3，其数据分布相对集中，较正态分布而言更陡峭。有效Fe偏态小于0，其数据向左偏态分布，表明比平均值大的含量值数目相对较多。有效Zn、有效B、有效Cu、有效Mn、有效Mo偏态均大于0，其数据向右偏态分布，表明比平均值小的含量值数目相对较多，其中有效Mn、有效Mo向右偏态较小，有效B、有效Cu向右偏态较大，有效Zn向右偏态最大。

利用ArcGIS软件对样点进行普通克里格插值，并将其插值结果赋到研究区所有评价单元上，得出奉新县各评价单元的土壤养分元素含量情况，根据隶属度函数求取各评价单元有益微量元素分值，之后运用最小因子法，选取隶属度最小值作为评价单元土壤有益微量元素最终得分，利用ArcGIS软件制作出奉新县6种有益微量元素含量分布图、稻田土壤健康评价有益微量元素分值图(图 5、图 6)。

由图 5可以看出，奉新县稻田有益微量元素含量中等较高，而有效Mo含量偏低，由于该单一指标评价采用最小因子法，评价结果受有效Mo含量影响较大，从而有益微量元素级别处于中等偏下。图 6表明，有益微量元素分值以0.2 ~ 0.4分居多，占73.40%，主要分布在东部山区、西部平原区的一些乡镇；其次是0.4 ~ 0.6分，为25.81%，主要分布在中部的干州镇、冯川镇、赤岸镇等；仅0.79%为0 ~ 0.2分。

2.2 稻田土壤健康综合评价结果

将单一指标评价分值结果采用耕地健康综合评价累加法进行计算，得出奉新县稻田土壤健康评价各评价单元健康系数范围为1 ~ 1.6，根据健康系数分级标准，计算得出各评价单元的稻田土壤健康级别，并统计出各级别下的稻田图斑面积、数量及比例分布状况(表 8)。

上述统计结果表明，全县非常健康的稻田共计4 475个评价单元，面积为8 527.71 hm²，占稻田总面积的25.33%；健康稻田评价单元为12 344个，面积为25 137.49 hm²，占比74.67%。

运用ArcGIS生成级别图，得出奉新县稻田土壤健康评价级别图(图 7)。

奉新县稻田土壤总体健康状况较好，健康级别在空间上呈现出中部高、东西部低的分布格局。1级主要分布在中部的干州镇、冯川镇、会埠镇、赤岸镇，少数散落在其他乡镇；2级则分布在除中部以外的东部、西部的各个乡镇。

根据单一指标评价结果，奉新县稻田灌溉水环境质量只有中部的少数部分级别较低；中部地区的土壤重金属级别与东部平原区、大多数西部山区处在同一级别；而中部的有益微量元素级别最高。综合上述评价结果，奉新县中部稻田土壤健康级别最高，耕地土壤健康状况最好。

3 讨论

土壤重金属的来源包括成土母质、工业、农业、交通、大气沉降等^[29-33]。阿吉古丽·马木提等^[34]研究表明，耕地土壤中As、Cd、Pb的来源主要受人类活动的影响。农业上施肥对土壤中重金属累积量产生直接影响^[35-36]，陈芳等^[37]通过长期定位试验表明，施用磷肥是导致研究区土壤重金属含量增加的主要原因。工业、生活污水进行农业灌溉也会导致重金属在土壤中积累，从而引发土壤重金属污染问题^[38]，影响耕地土壤健康。而根据调查统计资料，奉新县土壤自然本底较好，重金属元素含量较少；此外，县域内工业化程度不高，污染型工业较少；再者，该县注重生态农业，肥料的施用多选择不含重金属的肥料。因而奉新县稻田土壤重金属污染很少。

土壤养分的分布受气候、母质、地形、土壤类型等结构性因素与人为因素的影响^[39]。张毅等^[40]研究表明耕地土壤中有效Mo含量受到母质、海拔、地形、土壤类型的影响。江叶枫等^[41]的研究表明高程、坡度、地貌类型、土壤类型和常年耕作制度等显著影响着有效Fe、B、Cu、Zn、Mn 5种微量元素含量的空间变异。不同类型的土壤，其发育程度不同，导致微量元素的起始含量不同，因而对微量元素产生着影响^[42]。地形因素如高程主要通过影响水热条件的再分配过程影响土壤微量元素的空间分布^[43]。但在气候条件一致的区域内，生态系统经过长期自身演替，人类活动对土壤性质的影响不断增强^[44]。臧振峰等^[45]研究表明耕作、施肥以及灌溉等人为活动影响着土壤有益微量元素的空间分布。奉新县地势西高东低，由西部的山区到中部的丘陵地带再到东部的平原区，地形逐渐平缓，坡度下降。中东部低丘山地的复杂地形、气候、降雨、母质以及人类施肥行为等因素的综合作用，可能导致有益微量元素在该区的积累，使得中东部有益微量元素级别相对较高。

本研究选取了土壤重金属作为稻田土壤健康评价的指标因素之一，单美^[11]、谭晶今^[14]将土壤重金属指标作为耕地健康评价的重要因素，土壤重金属指标还被应用到与耕地土壤健康评价相关的农用地健康评价^[46-47]、农田生态系统健康评价^[48]、土地健康评价^[49]等研究中。灌溉水环境质量是本文的又一评价指标，灌溉水环境质量对耕地土壤产生的影响^[50-52]不可忽视，将灌溉水环境质量作为评价因素，丰富了评价指标体系。本文也将有益微量元素纳入到稻田土壤健康评价指标体系中，有益微量元素对作物生长发育有着重要影响^[53]，现有健康评价中鲜有涉及到有益微量元素，本文将其作为评价的指标因素之一，综合考虑了稻田土壤健康的有益与不利影响因素，从正反两面进行稻田土壤健康评价，具有一定的合理性与科学性。

本文将稻田地块作为评价单元进行小尺度领域的稻田土壤健康评价研究。评价单元小的评价精度相对越高^[11]，因而本研究相对于大尺度领域的评价而言，评价结果具有更高的精度。采用内梅罗综合污染指数法对灌溉水环境质量、土壤重金属进行评价，该方法是当前国内外计算综合污染指数最常用的方法之一，能够较全面地反映各种污染物对区域环境质量的协同作用，从而较准确地评价各类环境要素的污染程度^[54]，其计算过程简洁、易于操作^[55]，被广泛用于水质污染^[56-58]、土壤重金属污染^[59-62]等环境质量评价中。采用最小因子法对有益微量元素进行评价，常规的耕地评价多采用的是多因素综合评价法，该方法易于掩盖极值对评价结果的影响，降低了结果的精确性。与该方法相比，最小因子法易于操作，同时避免了最大值对评价结果的影响，且符合本研究的保守评价目的。最后运用累加法得出了稻田土壤健康综合评价结果，充分考虑稻田土壤健康的不利与有利因素，使得评价结果具有一定的全面性与合理性。同时，该方法简单、容易操作，便于实际工作中推广应用。

4 结论

1) 本研究单一指标评价结果表明，中部以冯川镇为主的0.85%的稻田灌溉水环境质量为最低3级，向外扩展等级逐渐升高，水质逐渐变好。冯川镇为该县县城所在地，县城较其他各乡镇而言社会工业化发展最快、人口密度最集中，人类行为导致灌溉水质受到一定污染，因而最靠近县城的灌溉水环境质量最差。全县99.31%稻田土壤重金属级别为1级，含量未超过标准值；仅0.69%以石溪乡为主的稻田受到轻微污染，级别为2级，这可能是由于自然成岩作用的影响。全县中部25.81%的稻田有益微量元素分值为最高0.4 ~ 0.6分，奉新县属亚热带湿润气候，中东部地处低山丘陵向低丘平原过渡地段，气候温暖、热量充沛、降雨量丰富，在气候、地形、母质、土壤以及人类耕作行为的综合影响下，可能导致有益微量元素级别相对最高。

2) 奉新县稻田土壤健康评价结果表明，奉新县稻田土壤健康状况良好，74.67%稻田土壤处在健康水平，25.33%稻田土壤为非常健康水平，该县自然环境条件优越，灌溉水环境质量总体中等偏下，土壤污染问题较少，综合有益微量元素的作用，该县稻田土壤总体健康状况较好。

耕地是人类之根本，为保障耕地资源的健康，在生活中，要控制农药化肥的施用，控制污水、污染物的随意排放，并逐步发展生态农业。当前，针对耕地土壤健康评价的研究尚浅，还处于探索阶段，未形成统一的评价标准，评价理论依据与方法也需要进一步探究，耕地土壤健康与其影响因子之间的关系还处于定性分析阶段，进一步的定量分析有待于在今后的研究中继续开展。