京津冀地区设施土壤中不同蔬菜对镉的累积特征

高 鑫，颜蒙蒙，曾希柏，白玲玉，王亚男，陈 清，赵会薇，苏世鸣*

（1.中国农业科学院农业环境与可持续发展研究所/农业部农业环境重点实验室，北京 100081；2.中国农业大学资源与环境学院，北京 100193；3.国家半干旱农业工程技术研究中心，石家庄 050051）

京津冀地区设施蔬菜产业已经发展了近30年，其在保障该地区蔬菜农产品供应、居民食品安全方面具有重要意义。该地区设施蔬菜生产水肥（尤其是粪肥和磷肥）投入大、复种指数高，特别是20世纪80年代初期推行的农田污灌等的影响，导致该地区设施土壤环境质量持续下降[1-4]，给该地区设施蔬菜安全生产带来很大影响。部分区域设施土壤出现重金属元素超标和农产品安全问题已引起较大关注[5-8]。相比其他污染元素，镉在设施土壤中的累积现象最为普遍[9]，且由于镉在土壤中有着相对较高的生物有效性，迁移性强，容易被吸收并转移到植物的不同部位，因此对蔬菜产品安全和人类健康所造成的威胁更大[10-11]。

筛选低镉吸收能力的蔬菜品种，并在镉超标设施土壤上开展替代种植是一种有效应对土壤重金属污染的修复技术。与瓜果类、根茎类蔬菜相比，叶类蔬菜对镉的累积吸收能力较强[12]，累积镉风险更高。市场上叶菜种类多，品种更替快，不同种类和品种之间对镉的累积差异亦较大，因此系统性地比较不同种类或品种叶菜对镉累积差异对于今后选择镉低吸收叶类蔬菜品种开展替代种植、阻控镉向叶类蔬菜累积具有重要意义。

本研究采用野外调查和田间小区试验相结合的方式开展，首先明确不同类型蔬菜对镉累积能力的差异，其次比较了11个不同种类蔬菜共64个蔬菜品种对镉的累积风险，并基于物种敏感度分布曲线[13]，确定了每种蔬菜对土壤中镉毒性的敏感性程度差异。相关结果对于今后通过种植替代方式减少京津冀地区设施叶菜重金属镉累积风险具有重要参考价值。

1 材料与方法

1.1 调查采样

选取河北青县、北京大兴、天津武清的设施大棚为调查区域，采用点对点采样法采集不同蔬菜可食部位样品及对应根区土壤样品。调查共采集54个大棚，其中河北青县采集大棚32个、北京大兴9个、天津武清13个。调查共采集到6个蔬菜种类，分别为黄瓜、豆角、西红柿、油麦菜、苋菜和白菜，其样本数分别为20、10、6、8、5个和5个。各大棚蔬菜和对应土壤样品均采用5点混合法采集，均匀混合后，所有样品带回实验室分析。

1.2 田间试验

试验共选择11种京津冀地区设施农业常见的叶类蔬菜，包括：菜心、油麦菜、生菜、茼蒿、小白菜、苋菜、菠菜、芥菜、苦苣、芝麻菜和咖啡菜。该11种叶类蔬菜包含64个蔬菜品种，具体编号及对应品种等信息见表1。田间小区试验于2017年在河北青县李家营村设施大棚开展。供试的两个大棚中土壤总镉含量依次为0.38 mg·kg-1和0.76 mg·kg-1，依据《温室蔬菜产地环境质量评价标准》（HJ 333—2006，当土壤pH＞7.5时土壤总镉含量限值为0.4 mg·kg-1），两个供试土壤的镉污染水平分别属于尚清洁和超标等级。大棚内各小区面积统一为15 m2（宽3 m、长5 m），各蔬菜品种处理3次重复。小区内种植密度依据各蔬菜种植特点设定。收获期采集各小区4～5株蔬菜地上部样品及对应根区土壤样品经均匀混合后带回实验室。植株样品经称取鲜质量后杀青、烘干称干质量，之后研磨，用于植株总镉含量分析；土壤样品经风干过筛后，分析土壤总镉含量。试验大棚内土壤有机质平均含量 23.7 g·kg-1，全氮 1.4 g·kg-1，有效磷11.3 mg·kg-1，速效钾103.8 mg·kg-1，土壤pH 7.6。

1.3 测定项目及方法

土壤中总镉的测定参考GB/T 17141—1997进行。称取土壤样品0.500 0 g左右，放置于聚四氟乙烯坩埚中，采用四酸消解，浓盐酸、浓硝酸、氢氟酸和高氯酸

添加比例为5∶10∶2∶1，消解液经定容后采用ICP-MS测定镉含量。植株总镉含量测定依据食品中镉的测定标准GB/T5009.15—2003进行。称取样品0.200 0 g于三角瓶中，10 mL混合酸（硝酸与高氯酸按9∶1比例混合）消解后定容，采用ICP-MS测定样品镉含量。上述样品分析过程中加入标准土壤（GBW07405、GBW07406） 和 标 准 植 物 样 品 （GBW10015、GBW10048）进行过程质量控制，标准品回收率范围为92.3%～105.6%。

表1 64种叶类蔬菜品种名称和生产商Table1 Variety names and purchasing places of 64 leaf vegetables

续表1 64种叶类蔬菜品种名称和生产商Continued Table1 Variety names and purchasing places of 64 leaf vegetables

1.4 数据统计

试验数据采用Excel 2013、SPSS Statistics进行统计作图并进行方差分析，采用最小显著性差异法（LSD）对数据进行多重比较，显著性水平为0.05。叶菜生物富集系数（Bio-concentration factors，BCFs）计算方法为叶菜可食部分镉含量与土壤总镉含量的比值；不同蔬菜品种对重金属镉的敏感性分布曲线（SSD）采用BurrliOZ 2.0软件进行，模型的具体计算方法采用Burr-Ⅲ型函数参数方程：

式中：x为叶菜1/BCF的实测值；b、c、k为函数的3个拟合参数。专业计算软件详见澳大利亚联邦科学和工业研究组织（http：//www.csiro.au/products/Burr III.ht⁃ml）。

2 结果与讨论

2.1 蔬菜类型间对镉累积能力的差异

通过调查京津冀地区主要设施蔬菜产区作物可食部位镉含量，比较了6种常见的大棚种植蔬菜镉含量差异，结果如图1所示。黄瓜、豆角、西红柿、油麦菜、苋菜和白菜可食部位镉含量分别为0.000 9、0.003 0、0.001 5、0.022 1、0.021 1、0.012 1 mg·kg-1鲜质量。总体上看，各蔬菜可食部位镉含量均明显低于国家食品卫生相关标准（0.2 mg·kg-1鲜质量）。不同蔬菜类型之间比较来看，瓜果类蔬菜（黄瓜、豆角和西红柿）可食部位镉含量显著低于叶类蔬菜（油麦菜、苋菜和白菜）。此外，6种蔬菜对土壤镉的生物富集系数具有相似规律，即瓜果类蔬菜对镉富集系数明显低于叶类蔬菜。蔬菜镉含量和生物富集系数大小顺序依次为苋菜＞油麦菜＞白菜＞豆角≈西红柿≈黄瓜。叶类蔬菜相对于瓜果类蔬菜具有更高的镉累积风险，瓜果类蔬菜应被选为京津冀地区镉污染设施土壤上推荐种植的蔬菜类型。该结果与前人研究相一致，如徐笠等[6]通过分析北京市设施菜田的大量调研数据发现，瓜果类蔬菜富集土壤镉的能力最弱，远小于叶类蔬菜和根类蔬菜。李霞等[8]的调查研究表明，不同类型蔬菜吸收镉的能力大小排序依次为叶菜类＞块茎类＞瓜果类。有研究表明，叶类蔬菜叶片表面积大，表面蒸发作用强，更容易富集吸收土壤和空气中的镉[14-15]，而瓜果类蔬菜由于土壤中重金属元素向其可食部位迁移距离较远，表现出较弱的镉富集能力[16]。此外，不同蔬菜对土壤镉的生理生化响应及其分子生物学过程差异亦是造成其富集镉能力不同的重要原因[17-19]。

图1 不同种类蔬菜镉含量和生物富集系数比较Figure1 Comparison of cadmium content and bio-concentration factor in different types of vegetables

图2 设施土壤上11种不同叶菜地上部镉含量Figure2 Cadmium content in 11 leafy vegetables grown in greenhouse soils

2.2 不同种类叶菜对镉累积能力的比较

叶类蔬菜对土壤镉的累积能力不仅存在显著的类型间差异，同时也在一定程度上表现出科属亲缘特性[20]。本研究通过田间小区试验分析了11种不同蔬菜对镉的累积能力差异。从图2中可以看到，不同叶菜种类之间以及不同镉污染土壤上叶菜地上部镉含量存在较大差异。不同种类叶菜比较来看，当土壤镉含量为0.76 mg·kg-1时，芝麻菜地上部镉含量最高，为0.13 mg·kg-1鲜质量，茼蒿最低为0.005 mg·kg-1鲜质量，两者相差约25倍。不同种类叶菜地上部镉含量大小排序依次为芝麻菜＞油麦菜＞菠菜=芥菜＞苋菜＞生菜＞菜心＞苦苣＞小白菜＞咖啡菜＞茼蒿；不同土壤镉浓度比较来看，相同种类蔬菜的镉含量与土壤中镉浓度存在一定的正相关关系，且两种镉浓度土壤上，不同叶菜地上部镉含量大小排序有很好的一致性，表明不同种类叶菜的镉累积能力差异在不同镉污染土壤上的稳定性很好。方华为[20]研究发现在盆栽土壤镉含量为2.12 mg·kg-1时，叶类蔬菜地上部植株镉含量大小排序依次为油麦菜＞生菜＞菜薹＞莴苣＞空心菜＞芹菜＞小白菜＞大白菜＞芥菜＞香菜＞苦苣＞香葱＞甘蓝＞木耳菜＞茼蒿＞苋菜＞韭菜。该结果与本论文研究存在一定差异，这可能与试验采用的供试土壤理化特征、镉浓度水平及供试蔬菜品种等方面的差异有关[21]。此外，本研究中所有蔬菜镉含量均未超过国家食品安全标准中的蔬菜镉限定值0.2 mg·kg-1鲜质量。

进一步对11种不同叶菜对土壤镉的生物富集系数进行聚类分析发现，在两种镉污染土壤上不同叶菜积累能力聚类分析具有良好的一致性。本论文选择镉超标设施土壤（0.76 mg·kg-1）上的相关结果进行分析，如图3所示。11种叶菜对镉的积累能力差异可划分为2类：第1类是菜心、苋菜、生菜、苦苣、油麦菜、小白菜、咖啡菜、茼蒿和菠菜，为镉中低累积类群；第2类是芥菜、芝麻菜，为镉高累积类群。

图3 11种叶菜类蔬菜对镉生物富集系数聚类分析Figure3 Cluster analysis of Cd bio-concentration factor of 11 leafy vegetables

图4 两种浓度土壤上64种叶菜品种地上部镉含量Figure4 Cd contents in 64 leafy vegetables grown in greenhouse soils with two levels of Cd concentration

2.3 叶类蔬菜品种对镉累积能力的差异

为明确不同蔬菜品种对镉累积能力的差异，进一步分析了上述11个蔬菜种类涉及的64个叶菜品种对镉累积能力的差异。结果如图4所示，64个叶菜品种地上部镉含量均未超过国家食品卫生标准，表明该区域设施蔬菜生产相对清洁。但不同叶菜品种间对镉的累积状况存在显著差异。如芝麻菜和紫油麦菜有着相对较高的镉累积能力，其中以芝麻菜对镉的累积量最大，其地上部镉含量为0.13 mg·kg-1；而茼蒿、汉堡速生咖啡菜、翠英256和特选板叶茼蒿等对镉的累积量较低，其地上部镉含量＜0.01 mg·kg-1，远低于国家食品安全质量标准规定的0.2 mg·kg-1。该结果表明，今后在镉污染程度严重的设施土壤上进行蔬菜栽培时，需要关注不同蔬菜品种对镉累积的差异，通过合理选择蔬菜品种可以有效减少蔬菜镉超标风险。

此外，相同种类叶菜不同品种间比较来看，仅菠菜在0.76 mg·kg-1土壤中品种间差别较小，品种间镉含量差异相对稳定，其他蔬菜品种间镉含量存在明显差异，该结果表明今后开展叶类蔬菜替代种植时应关注相同蔬菜种类不同品种间差异。另外，两种土壤镉浓度比较来看，64个叶菜品种间蔬菜镉含量大小关系排序具有一致性。通过相关性分析发现，两种镉污染设施土壤上64个叶菜品种对镉累积量具有极显著相关性（n=64，R2=0.528 1，P＜0.01），表明不同叶菜品种的镉累积能力大小在不同程度的镉污染土壤上具有较好的稳定性。本研究表明不同叶类蔬菜对镉累积能力的差异在不同镉污染设施土壤上具有很好的重现性，可作为该区域高浓度镉污染土壤开展替代种植的有效依据。

图5 64种蔬菜品种对重金属镉敏感性分布曲线（SSD）Figure5 Distribution curve of sensitivity to Cd of 64 vegetable varieties

2.4 64个叶菜品种对土壤镉毒性敏感性分布曲线

敏感性分布曲线是通过选择某个概率分布并拟合SSD曲线的方法，描述某种重金属元素对一系列物种的毒性。基于生态风险评价角度分析重金属污染对不同物种的危害程度，进而对其生态风险进行评估，评价不同物种对重金属污染的敏感程度[22]。本研究应用较新的概率分布函数模型（Burr-Ⅲ模型）[23]比较了上述京津冀地区设施生产中常见的64个叶类蔬菜品种对土壤镉的毒性敏感性分布曲线（图5），结果表明不同叶类蔬菜对土壤镉的敏感程度存在概率分布差异。具体来看，咖啡麦当菜、蔬菠2号、芝麻菜和佳园芥蓝对土壤镉的毒性最为敏感，其敏感性分布频次小于5%；其次为汉斯203紫霞、紫油麦菜、王中王尖叶油麦菜，属于较敏感品种，其敏感性分布频次大于5%，小于10%；而茼蒿属于不敏感品种。基于保护所选64种叶类蔬菜品种中95%的蔬菜产品安全考虑，叶菜最小环境有害HC5-BCFs值（HC5，Hazardous Concentration for 5%of species）为0.19。以《食品安全国家标准》（GB 2762—2017）规定的0.2 mg·kg-1鲜质量为基准，反推出土壤镉的浓度上限为1.05 mg·kg-1。应尽量避免种植咖啡麦当菜、蔬菠2号、芝麻菜和佳园芥蓝等蔬菜品种。国内外相关学者亦开展了相关研究，如李晓婧等[24]比较了10种叶菜对镉毒害的敏感程度，发现苋菜和油菜属于较为敏感品种，芹菜和茼蒿属于较不敏感品种，芹菜和茼蒿适合在镉超标土壤中开展替代种植。该结果与本论文研究具有良好的一致性。

3 结论

（1）瓜果类蔬菜对土壤镉的累积能力明显低于叶类蔬菜，11种叶类蔬菜地上部镉含量大小排序依次为芝麻菜＞油麦菜＞菠菜=芥菜＞苋菜＞生菜＞菜心＞苦苣＞小白菜＞咖啡菜＞茼蒿。

（2）不同叶类蔬菜品种对土壤镉累积能力差异较大。芝麻菜地上部镉含量最高，而茼蒿、汉堡速生咖啡菜、翠英256和特选板叶茼蒿地上部镉累积量＜0.01 mg·kg-1鲜质量，适合在京津冀镉污染设施叶菜土壤上开展替代种植。

（3）芝麻菜、蔬菠2号、咖啡麦当菜、佳园芥蓝属于土壤镉毒性敏感品种，其物种敏感性分布频次小于5%，基于保护所选64种叶类蔬菜品种中的95%安全考虑，应避免种植该类蔬菜品种。