环境科学  2025, Vol. 46 Issue (4): 2580-2589   PDF    
宁夏农田土壤和农产品镉含量特征、健康风险评价及其环境基准
陈林1,2,3, 杨玲1,2,3, 李虹4, 梁翔宇4, 马建军4, 王君梅4, 王金保4, 马俊花4, 贾彪5, 薛斌1,2,3, 庞丹波1,2,3, 马琨1,2,3     
1. 宁夏大学生态环境学院, 银川 750021;
2. 宁夏大学西北退化生态系统恢复与重建教育部重点实验室, 银川 750021;
3. 宁夏大学西北土地退化与生态恢复国家重点实验室培育基地, 银川 750021;
4. 宁夏农业环境保护监测站, 银川 750002;
5. 宁夏大学农学院, 银川 750021
摘要: 准确评估土壤-农产品Cd健康风险有助于提升农田污染防治决策的科学性和精确性. 通过采集宁夏农田740对土壤-农产品协同监测样品, 分析土壤-农产品Cd污染特征, 评估土壤和摄入农产品对人类健康风险等级, 推导研究区农田土壤环境基准值. 结果表明, 研究区土壤pH值> 7.3, 呈碱性. 土壤ω(Cd)在0.020~0.472 mg·kg-1之间, 低于我国农用地土壤污染风险管控值. 农产品中其他谷物和蔬菜的ω(Cd)平均值较高, 分别为0.013 mg·kg-1和0.011 mg·kg-1, 马铃薯和水果的ω(Cd)平均值最低(均为0.002 mg·kg-1). 不同农产品Cd的生物富集系数 < 1, 说明所有农产品没有积累效应. 土壤途径对人体的单项非致癌风险指数低于1, 且单项致癌风险指数也无致癌风险, 但研究区土壤Cd的非致癌和致癌风险指数均表现为儿童高于成人. 农产品摄入Cd对成人和儿童无健康风险. 由敏感性分布法推导出的玉米、小麦、水果和其他谷物地土壤安全生产临界值分别为0.503、0.336、0.234和0.342 mg·kg-1, 小于土壤环境质量标准风险筛选值, 本研究推导出的玉米、小麦、水果和其他谷物地土壤Cd阈值较标准严格;水稻和蔬菜地的土壤安全生产临界值为1.057 mg·kg-1和0.740 mg·kg-1, 高于标准中对应的风险筛选值, 推导出的水稻和蔬菜地土壤Cd阈值相较标准要宽松;马铃薯地的土壤安全生产临界值和标准基本一致. 因此, 统一用现行土壤环境质量标准的筛选值并不完全适用于研究区不同类型农产品, 需根据不同农作物中重金属的实际含量对土壤环境基准值进行调整.
关键词: 宁夏      镉(Cd)      农产品      健康风险      土壤环境基准     
Cadmium Content Characteristics, Health Risk Assessment, and Soil Environmental Benchmarks in Farmland Soil and Agricultural Products in Ningxia
CHEN Lin1,2,3 , YANG Ling1,2,3 , LI Hong4 , LIANG Xiang-yu4 , MA Jian-jun4 , WANG Jun-mei4 , WANG Jin-bao4 , MA Jun-hua4 , JIA Biao5 , XUE Bin1,2,3 , PANG Dan-bo1,2,3 , MA Kun1,2,3     
1. School of Ecology and Environment, Ningxia University, Yinchuan 750021, China;
2. Key Laboratory for Restoration and Reconstruction of Degraded Ecosystem in Northwest China of Ministry of Education, Ningxia University, Yinchuan 750021, China;
3. Breeding Base for State Key Laboratory of Land Degradation and Ecological Restoration in Northwest China, Ningxia University, Yinchuan 750021, China;
4. Agriculture Environmental Protection Monitoring Station of Ningxia, Yinchuan 750002, China;
5. School of Agriculture, Ningxia University, Yinchuan 750021, China
Abstract: A precise assessment of the health risks posed by cadmium (Cd) contamination in agricultural products from regional farmland is crucial for informed pollution prevention and control decisions. This study collected 740 pairs of soil and agricultural product samples from farmland in Ningxia, focusing on Cd content detection. The evaluation encompassed Cd pollution characteristics in soil and agricultural products and assessed the associated health risk for human ingestion. Additionally, farmland soil environmental benchmark values specific to the study area were established. The results showed that soil pH levels across the study area were consistently alkaline, with values exceeding 7.3. Cd content in the soil fell within a range of 0.020 to 0.472 mg·kg-1, lower than the pollution risk control value of agricultural soil (pH > 7.5) in China (GB 2762-2017). Average Cd content in grains and vegetables was relatively high, with values of 0.013 and 0.011 mg·kg-1, respectively. Conversely, potatoes and fruits exhibited the lowest average Cd content at 0.002 mg·kg-1 each. The biological accumulation coefficient of Cd in different agricultural products was less than 1, indicating no accumulation effect across all agricultural products. The single non-carcinogenic risk index (HQ) of the soil pathway to the human body was lower than 1, and the single carcinogenic risk index also had no carcinogenic risk (CR < 10-4), but the non-carcinogenic and carcinogenic risk indices of Cd intake via soil pathways remained below significant thresholds, yet they were higher in children compared to adults. The intake of Cd through agricultural products posed no discernible health risks to either adults or children. The critical soil Cd thresholds for corn, wheat, fruit, and other grain fields derived from the species sensitivity distribution (SSD) method were 0.503, 0.336, 0.234, and 0.342 mg·kg-1, respectively, which were lower than the risk screening values of soil environmental quality standards (GB 15618-2018) in China. Therefore, the Cd thresholds for corn, wheat, fruit, and other grain fields derived from this study were stricter than the standards. The critical values for soil Cd in rice and vegetable fields were 1.057 and 0.740 mg·kg-1, respectively, which were higher than the corresponding risk screening values in the standard. Specifically, critical values for rice and vegetable fields were found to be higher, suggesting a relaxed threshold compared to the standard. Conversely, critical values for soil Cd in potato fields aligned closely with the existing standard. This study underscores the importance of tailoring soil environmental benchmarks to specific agricultural products, advocating for adjustments based on actual content of heavy metals in different agricultural products.
Key words: Ningxia      Cd      agricultural products      health risk      soil environmental benchmark     

土壤作为农业生态系统的核心成分和物质能量循环的枢纽, 是各种重金属污染物的主要载体[1]. 据2014年发布的《全国土壤污染状况调查公报》显示, 我国土壤总点位超标率为16.1%, 主要污染物为Cd(点位超标率达7%)[2]. 土壤中Cd的活性较高、毒性较强, 且具有生物蓄积性和不可被生物降解性[3], 在低浓度条件下也可经食物链富集威胁人体健康[4], 而且长期Cd暴露将会产生毒害作用[3], 由此导致的食品和生态安全问题已受到广泛关注. 因此, 精准评估和识别农田土壤-农产品Cd污染的健康风险, 是应对粮食安全和保障人体健康的重要举措.

与土壤背景值相比, 中国农田土壤中Cd的累积量最大, 超过对应土壤背景值177.9%[5]. 在长江中游地区, Cd是潜在生态风险最大的重金属[4]. 黄河流域农田有21.69%的样点Cd含量超出《土壤环境质量农用地土壤污染风险管控标准》(GB 15618-2018)(pH > 7.5)的筛选值, 为偏中度污染[6]. 我国地域辽阔, 农田土壤环境质量差异较大, 不同作物对重金属的富集能力也有明显差异, 而现行的土壤环境质量标准风险筛选值和管控值制订时均未考虑其生态风险[2]. 用统一的标准来界定土壤中某种重金属元素的临界值缺乏准确性和科学性, 会造成实际评估中出现超标不对应的情况[7]. 加之缺乏地方性标准作为补充, 极大影响和限制了土壤环境质量定义的严肃性和科学可靠性. 因此, 开展土壤-作物系统中重金属的分布特征和富集规律研究, 针对不同作物类型制定相应的土壤安全限量标准, 对保障我国农产品安全和人体健康具有重要意义.

本研究基于宁夏全区农田土壤-农产品协同监测结果, 分析不同类型农产品和土壤的Cd含量, 开展农产品和土壤Cd污染评价, 探讨摄入农产品和土壤途径对人类健康风险, 并采用累积概率分布函数构建不同农产品的物种敏感度分布曲线, 推导不同农产品产地土壤Cd生态安全环境基准, 以期为宁夏农田土壤Cd的生态风险评估和土壤环境质量标准的制定提供数据支撑和科学参考.

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

宁夏回族自治区属中温带干旱、半干旱气候, 年平均气温为5~9℃;年平均无霜期105~163 d;年平均降水量约为200~680 mm, 主要集中在夏季和秋季;年蒸发量1 214~2 803 mm[8];年降水量在400 mm以下的干旱地区占全区总面积的80%;年日照时间约为2 194~3 082 h. 主要种植玉米、小麦、水稻、马铃薯、冷凉蔬菜和小杂粮等作物, 大部分为一年一熟.

1.2 样品采集和处理测定

考虑到宁夏土地利用特征和农作物种植情况, 同时兼顾农作物种植空间的均匀性, 基于GPS定位采样点, 于2017~2023年协同采集740对土壤(0~20 cm)和农产品(可食部分)样品, 其中, 玉米地415组、小麦地67组、水稻地50组、马铃薯地20组、水果地(包括苹果、葡萄、西瓜等)20组、蔬菜地(包括油菜、小白菜、西蓝花、芹菜、茄子、南瓜、韭菜、辣椒等)147组和其他谷物(包括荞麦、高粱、大豆等)地20组. 每个土壤样品是从田块中间位置、半径50 m范围内、由“S”形组合5个表层土样混合而成, 土壤样品重量不低于1 kg, 经自然风干, 剔除碎石、植物根系等杂物, 并用四分法取样后过10目的尼龙筛, 保存备用. 农产品采集点与土壤样品采样点一一对应(同样是740个), 采样时避开株体过大过小、遭受病虫害或机械损伤以及道路边的植株. 虽然籽粒的重金属含量与根、茎、叶中重金属含量存在一定的相关性, 但没有确定的规律[9], 因此, 本文在衡量农产品对人体的健康风险时, 选用所采集农产品的可直接食用部分进行检测. 在作物成熟期(6月或9月), 用四分法混合后获得约0.5 kg的农产品.

土壤pH值由雷磁pH计(PHS-3C, 中国上海)按土∶水(1∶2.5)进行测量. Cd的含量使用石墨炉原子吸收分光光度计(AA-7000 Shimadzu, 日本)进行测定[10,11], 样品回收率为100%. 样品制备和分析过程均使用分析试剂空白. 对10%的土壤样品同时进行了重复样品分析, 以上重复样品的相对标准偏差均小于5%.

1.3 评价方法 1.3.1 土壤重金属污染评价

采用地累积指数法(Igeo)评价宁夏农田中重金属污染程度, 该方法是由德国科学家Muller提出用于研究重金属污染程度的定量指标[12], 是目前国内外广泛应用于评价土壤重金属污染程度的一种定量评价方法, 不仅能反映自然成岩作用下重金属自然分布特征和累积程度, 而且可以判断出人为活动对重金属污染的影响, 是区分人为活动影响的重要参数[13,14], 计算公式如下:

(1)

式中, Igeo为地累积指数;CCd为实测Cd元素在土壤中的含量(mg·kg-1);Bn为Cd的地球化学背景值, 采用宁夏的土壤背景值0.112 mg·kg-1[15];1.5为考虑成岩作用可能引起背景值波动而设定的校正常数[16].

地累积指数Igeo划分标准见表 1.

表 1 地累积指数分级标准 Table 1 Classification standard of soil geo-accumulation index

1.3.2 土壤重金属生态风险评估

生态风险指数(ecological risk index, ERI)用于分析土壤重金属的污染程度和评估潜在的生态风险, Hakanson提出的潜在生态风险指数法[17]是目前应用最为广泛的重金属质量评价方法之一[18], 可从生态毒性角度进行重金属的潜在生态风险评估, 并根据结果划分单个重金属指标或综合指标的生态风险等级[19], 其计算公式为:

(2)

式中, C为重金属Cd的实测含量;B为宁夏土壤Cd的环境质量标准值(0.112 mg·kg-1[15]T为土壤重金属毒性反应因子, Cd取值为30[20]. 生态风险指数分为ERI < 40、40≤ERI < 80、80≤ERI < 160、160≤ERI < 320和ERI≥320, 分别代表轻微、中等、较高、高和极高风险等级.

1.3.3 人体健康风险评估

重金属主要通过经口摄入土壤、吸入土壤颗粒、皮肤接触和食物链摄取等暴露途径进入人体[21], 因此本研究分别从土壤途径和农产品摄入途径开展人体健康风险评估.

1.3.3.1 土壤对人类健康的风险评估

本研究采用美国环保署推荐的健康风险评价模型分析土壤Cd对不同人群的健康风险.

(3)
(4)
(5)

式中, ADD土壤_口手、ADD土壤_皮肤和ADD土壤_呼吸为不同暴露途径的长期日均暴露量, mg·kg-1·d-1;IR1为每日摄入速率, 成人取值100 mg·d-1, 儿童取值200 mg·d-1[21];IR2为土壤吸入速率, 成人取值20 mg·d-1, 儿童取值7.65 mg·d-1[22];EF为土壤暴露频率, 取值350 d·a-1[1];ED为土壤暴露持续时间, 成人取值25 a, 儿童取值6 a[1];BW为居民平均体重, 成人取值56.8 kg, 儿童取值15.9 kg[21]Cs为实测土壤Cd含量, mg·kg-1;CF为单位转化因子, 取值10-6 mg·kg-1;SA为皮肤暴露面积, 成人取值5 700 cm2, 儿童取值2 800 cm2[1,22];AF为皮肤黏附系数, 成人取值0.07 mg·cm-1·d-1, 儿童取值0.2 mg·cm-1·d-1[1,22];ABS为皮肤吸收因子, 成人取值0.003, 儿童取值0.001[22];PEF为颗粒物释放因子, 成人和儿童均取值1.36×109 mg3·kg-1[1,22];AT为平均总的暴露时间, ED×365 d(非致癌), 70×365 d(致癌)[1,22].

土壤重金属对人体产生的健康风险可分为非致癌健康风险评价和致癌风险评价, 计算公式如下:

(6)
(7)

式中, HQ为非致癌重金属Cd的健康风险指数, 当HQ < 1时, Cd对人体不存在非致癌风险. ADDi为非致癌重金属Cd的第i种暴露途径的日均暴露量, RFDi为非致癌重金属Cd的第i种暴露途径的参考剂量, 口手、皮肤和呼吸途径的参考剂量分别为1.0×10-3、2.5×10-5和1.62×10-7 mg·kg-1·d-1[23]. CR为致癌风险指数, 分级标准为[24]:CR > 10-4时, 致癌风险较高;10-6 < CR≤10-4时, 表示人体可接受水平的致癌风险;CR≤10-6时无致癌风险. SF为斜率致癌因子系数, 取值6.1 kg·d-1·mg-1[20].

1.3.3.2 农产品对人类健康的风险评估

在农田重金属健康风险评估中, 应重视食物消费暴露途径[21]. 健康风险指数(health risk index, HRI)适用于评估通过食物摄入重金属对人体健康造成的风险, 提出了目标危害方法[25]. 因此, 本文对经口摄入的农产品重金属暴露途径进行单独评估, 计算公式如下:

(8)

式中, FIR为农产品日摄入量, 为了便于统一比较不同农产品重金属的摄入风险, 将所有农产品归类为谷物、蔬菜和水果, 其中成人和儿童的谷物FIR值分别为0.400 kg·d-1和0.200 kg·d-1[26], 成人和儿童的蔬菜FIR值分别为0.345 kg·d-1和0.232 kg·d-1[26], 水果的FIR值分别为0.275 kg·d-1和0.150 kg·d-1[27]Ccrop为农产品中Cd的实测含量(mg·kg-1);EF为暴露频率, 取值为365 d·a-1ED为暴露持续时间, 成人取值25 a, 儿童取值6 a[1]WAB为平均体重, 成人取值61.75 kg, 儿童取值32.75 kg[28]TA为平均时间(d). TAEF×ED计算得出. RFD为重金属的每日允许摄入量, Cd的RFD值为0.001 mg·kg-1·d-1[25,29]. 当HRI > 1时, 可能会对人体健康造成不良影响;HRI < 1时, 表明对人体健康无不良影响.

1.3.4 土壤环境基准值推导

土壤生态环境基准是指土壤中的污染物既不影响农产品产量和生物的生存发展, 又不致使水体污染的阈值[7]. 物种敏感性分布法(species sensitivity distribution, SSD)是推导土壤环境基准的国际主流方法, 已被多个国家和机构确定为制订土壤环境质量标准的方法[2]. 该方法基于一系列物种或相同物种的不同品种间对某一种污染物敏感性差异而提出来的一种生态风险评价方法[7,22], 通过概率分布函数将不同物种间的毒理数据外推, 实现污染物在生物群落甚至生态系统水平上的风险评估[30]. 本研究采用log-logistics模型方程对数据进行拟合分析, 通过测定土壤和农产品重金属含量, 计算农产品富集系数(BCF), 用于估算作物从土壤中累积金属的能力[31], 如公式(9), 然后将其倒数作为概率分布指标, 农产品重金属的1/BCF对土壤重金属富集效应的敏感分布遵循“S”型曲线分布[22], 因此利用公式(10)拟合SSD曲线. 最后参照《食品安全国家标准食品中污染物限量》(GB 2762-2022)规定的各类农产品限量, 利用公式(11)反演出基于保护95%和5%农产品类别安全的HC5和HC95, 并将HC5和HC95定义为土壤安全生产的临界值和警戒值[32].

(9)
(10)
(11)

式中, Ccrop为农产品中的重金属含量, Cs为土壤中的重金属含量. 当BCF≤1时, 表明作物只能吸收而不能积累重金属;当BCF > 1时, 表明作物可以积累重金属[31]. x为1/BCF, y为农产品对应的累积概率, abcx0为常数. Cr为农产品Cd限量值, 分别按玉米(0.1 mg·kg-1)、小麦(0.1 mg·kg-1)、水稻(0.2 mg·kg-1)、马铃薯(0.1 mg·kg-1)、其他谷物(0.1 mg·kg-1)、其他蔬菜及豆类(0.2 mg·kg-1)和水果(0.05 mg·kg-1)进行取值.

2 结果与讨论 2.1 土壤pH统计分析

研究区农田土壤pH值和有机质含量特征统计结果见表 2. 监测点中只有兴庆区的一个西红柿和贺兰县一个西瓜样地的土壤pH值小于7.5, 分别为7.3和7.4, 呈中性偏碱性. 其他738个监测点(占所有采样点的99.73%)土壤pH值均 > 7.5, 说明研究区农田土壤呈碱性, 一定程度上会抑制重金属的活性, 有助于降低土壤重金属的环境风险[33]. 考虑到不同作物类型土壤pH的平均值> 8.3, 因此, 本研究在土壤重金属的相关分析中, 不进行土壤pH值条件下的区分, 均按照土壤pH > 7.5的相关标准值和限量值开展探讨.

表 2 研究区农田土壤pH值特征统计 Table 2 Characteristic statistics of pH value

2.2 土壤Cd污染评价 2.2.1 土壤Cd含量状况

研究区农田中不同农产品土壤(0~20 cm)样品进行描述性统计分析可知(表 3), 所有样地中土壤ω(Cd)在0.020~0.472 mg·kg-1, 平均值(0.187 mg·kg-1)高于宁夏土壤背景值(0.112 mg·kg-1)和中国表层土壤背景值(0.097 mg·kg-1), 但均低于我国农用地土壤Cd风险管控值(0.600 mg·kg-1). 其中水稻地土壤ω(Cd)平均值最高, 为(0.210±0.071)mg·kg-1, 水果地土壤ω(Cd)平均值最低[(0.144±0.051)mg·kg-1]. 从最大值和变异系数统计结果可以看出, 各农产品监测点均有不同程度的累积, 部分点位土壤Cd含量超过了我国农用地土壤污染风险筛选值(0.300 mg·kg-1), 但均低于美国(0.770 mg·kg-1)、加拿大(1.400 mg·kg-1)和芬兰(10.000 mg·kg-1)等国家的土壤Cd指导值[3].

表 3 研究区农田土壤Cd含量统计1) Table 3 Statistics of Cd content in farmland soil in study area

从变异系数来看, 除马铃薯地土壤Cd含量变异系数较低(19.88%)外, 其他农产品样地土壤Cd均具有较高程度变异(> 30%), 均属于中等程度变异[20], 说明各样地间土壤Cd含量离散程度较大, 在空间上分布不均匀, 可能受到外界人为因素不同程度的影响[34]. 但与其他地区研究结果相比[7,32,35 ~ 37], 其变异程度相对较弱. 这与采样尺度有关, 采样区域面积越小, 其空间的结构性变异越弱[35].

不同农产品土壤Cd平均值均高于宁夏土壤背景值, 说明有不同程度的累积富集效应;与全国表层土壤背景值相比, 玉米、小麦、水稻、马铃薯、水果、蔬菜和其他谷物地分别是背景值的1.95、1.96、2.16、1.75、1.48、1.90和1.65倍;而与我国农用地土壤(pH > 7.5)污染风险管控值相比, 所有农产品的土壤Cd含量均低于管控值, 没有超标管控风险的监测点.

2.2.2 土壤Cd污染评价

为了了解研究区农田土壤Cd累积可能引起的环境风险, 应用地累积指数(Igeo)开展研究区土壤Cd污染程度评价. 由图 1可知, 研究区农田土壤Cd地累积指数平均值大小顺序依次为:水稻(0.234) > 玉米(0.089) > 小麦(0.078) > 蔬菜(0.046) > 马铃薯(-0.008) > 其他谷物(-0.135) > 水果(-0.307), 平均值为0.069, 说明研究区农田土壤Cd污染程度较低, 大部分属于轻度污染及以下(无污染)等级, 仅有小部分玉米(8个)、小麦(2个)、水稻(2个)和蔬菜(4个)地的土壤Cd处于中度污染水平. 而且研究区土壤Cd的地累积指数变异系数高达735.00%, 说明各监测点空间异质性较强.

蓝色虚线为无污染和轻度污染分界线;黑色虚线为轻度污染和中度污染分界线;红色虚线为中度污染和重度污染分界线 图 1 研究区不同农产品产地土壤Cd地累积指数 Fig. 1 The Igeo of soil Cd in different agricultural product producing areas of the study area

2.2.3 土壤Cd生态风险评估(ERI)

考虑到区域土壤Cd含量的强空间异质性, 因此进一步采用能重点反映生物毒性较高的Cd对土壤环境质量影响的生态风险指数(ERI), 开展土壤Cd潜在生态风险水平的评估. 由图 2可知, 生态风险指数结果与地累积指数相似, 大部分监测点土壤Cd处于中等风险等级及以下, 玉米、小麦、水稻和蔬菜地土壤Cd处于较高风险等级的比例分别为5.78%、4.48%、16.00%和6.12%, 马铃薯、水果和其他谷物地土壤Cd没有较高风险等级, 而且研究区无高和极高风险等级的土壤点位. 综上所述, 不同农产品类型土壤Cd评价结果基本一致, 研究区农田土壤Cd生态风险普遍处于安全水平. 但仍需坚持开展相关监测工作, 避免区域出现新的污染源以加重区域土壤Cd的累积, 出现新的生态风险.

蓝色虚线为轻微和中等风险等级分界线;黑色虚线为中等和较高风险等级分界线;红色虚线为较高和高风险等级分界线 图 2 研究区不同农产品产地土壤Cd生态风险指数 Fig. 2 ERI of soil Cd in different agricultural product producing areas of the study area

2.3 农产品Cd含量和富集特征评价 2.3.1 农产品Cd含量状况

研究区不同农产品中Cd含量统计结果见表 4. 农产品ω(Cd)分布在0.000~0.150 mg·kg-1之间, 平均变异系数高达137.89%, 具有极高的离散性, 说明不同农产品Cd含量的变异程度较大. 而且同一农产品间的变异系数(CV)也较高, 马铃薯和水果Cd含量属于中等变异程度(10% < CV < 100%), 玉米、小麦、水稻、蔬菜和其他谷物属于强变异(> 100%), 说明相同农产品的Cd含量在不同监测点间也具有较大差异. 与土壤Cd含量不同, 农产品中其他谷物和蔬菜的ω(Cd)平均值较高, 分别为0.013 mg·kg-1和0.011 mg·kg-1, 其次是小麦(0.009 mg·kg-1), 水稻和玉米中的ω(Cd)平均值较低(分别为0.004 mg·kg-1和0.003 mg·kg-1), 马铃薯和水果ω(Cd)平均值最低(均为0.002 mg·kg-1), 以上均远低于国家对食品中Cd含量的限量值(GB 2762-2017).

表 4 研究区不同农产品Cd含量统计1) Table 4 Statistics of Cd content in different agricultural products in the study area

2.3.2 Cd富集特征

土壤重金属含量高, 但对应的农产品重金属含量不一定高, 这与农作物对重金属的积累能力有关[36], 而衡量农作物对重金属的积累可以通过富集因子(BCF)来判断, BCF越小, 则农产品对重金属的富集能力越弱, 抗土壤重金属污染的能力则越强. 而且BCF可在一定程度上消除区域环境因子影响差异[41]. 从研究区不同农产品Cd的生物富集系数可知(图 3), 所有农产品只是吸收Cd(BCF < 1), 远低于积累指示值(BCF > 1), 说明监测农产品对Cd并没有积累效应. 其中, 其他谷物和蔬菜的平均BCF较高, 分别为0.083和0.057, 其次是小麦(0.053), 但这3种农产品BCF的变异程度也较大, 依次为9.5%、8.5%和6.8%. BCF最小的是马铃薯(0.013), 其变异系数也最小, 仅为0.8%. 研究区水稻的BCF为0.022, 未表现出高富集性. 马铃薯和水果的BCF值较低, 玉米和蔬菜的高值离散点较多. 虽然土壤中Cd主要以可交换态或碳酸盐结合态等易于迁移的形态存在, 通过Fe2+、Ca2+和Zn2+等二价离子转运易于被作物吸收, 但这也与采集点位数和采样频次影响有关[35], 同时植物根系对Cd的吸收速率和转运方式不同、不同作物对重金属耐受性具有强弱差异以及生物量大小等, 均是导致不同作物地上部Cd积累的重要决定因素[42,43]. 此外, 和研究区土壤pH值普遍较高(> 7.5)有关, 较高的土壤pH值有利于土壤中重金属从离子交换态转变成碳酸盐沉淀, 从而形成稳定的形态, 使重金属的有效性下降[36].

红色虚线为作物不能积累和可以积累重金属分界线 图 3 研究区不同农产品Cd的生物富集系数箱型图 Fig. 3 Boxplot of biological accumulation coefficient of Cd in different agricultural products in the study area

2.4 土壤和摄入农产品对人类健康风险评价 2.4.1 土壤对人类健康风险评估

目前, 开展农田土壤重金属的人体健康风险评估时, 通常参考建设用地的评估方法, 考虑经口摄入土壤、皮肤接触和吸入土壤颗粒这3种暴露途径[21,44]. 因此, 本研究分别采用致癌和非致癌风险评价法分析土壤途径对成人和儿童的健康风险.

不同农产品产地土壤途径对人群的单项非致癌风险指数(HQ)及单项致癌风险指数(CR)分别如图 4图 5所示. 可以看出, 成人和儿童的HQ大小排序均为:水稻 > 小麦 > 玉米 > 蔬菜 > 马铃薯 > 其他谷物 > 水果, 而且HQ远低于1, 所以土壤Cd对成人和儿童均不存在非致癌风险. 而土壤Cd含量对成人和儿童的CR值也无致癌风险(CR < 10-4), 但玉米、小麦、水稻、马铃薯和蔬菜地中土壤Cd对儿童的平均CR值介于10-6和10-4之间, 即对儿童有可接受水平致癌风险等级, 而对成人则大部分为无致癌风险(CR < 10-6).

图 4 研究区不同农产品土壤途径非致癌风险指数箱型图 Fig. 4 Boxplot of non-carcinogenic risk index (HQ) for soil pathways of different agricultural products in the study area

红色虚线为人体可接受水平致癌风险和无致癌风险的分界线 图 5 研究区不同农产品产地土壤途径致癌风险指数箱型图 Fig. 5 Boxplot of carcinogenic risk index (CR) for soil pathways of different agricultural products in the study area

整体看来, 就儿童而言, 7种农产品产地土壤Cd的非致癌风险指数和致癌风险指数均高于成人, 即在土壤暴露环境中, 儿童通常比成人脆弱, 更容易受到土壤途径Cd的危害.

2.4.2 摄入农产品对人类健康风险评价

与土壤摄入和吸入等暴露途径相比, 摄入农产品途径对人类健康存在更大的直接性和不确定性. 在计算食物链摄入途径的健康风险时, 需要通过实地调查获得农田的利用类型、农产品的种类及其重金属含量等参数[21]. 从图 6可以看出, 通过农产品摄入Cd对成人和儿童的健康风险指数均 < 1, 说明研究区域农产品Cd含量对人体不存在健康风险, 即没有食源性危害. 排除极端离散值差异的影响, 由成人和儿童的健康风险平均值来看, 玉米、小麦、水稻、马铃薯、水果、蔬菜和其他谷物摄入Cd对成人和儿童的健康风险无显著差异, 而且不同农产品类别间对成人和儿童的影响同样不大.

红色虚线为对人体健康有无不良影响的分界线 图 6 研究区不同农产品Cd健康风险指数箱型图 Fig. 6 Mean index values of human health risk of Cd in different agricultural products in the study area

2.5 土壤环境基准值推导

土壤环境基准是指不考虑社会、经济和技术等人为因素, 土壤中某一物质对所在土壤生态系统的存在对象(人或其他生物等)不产生不良或有害影响的最大剂量(或称无有害效应剂量)或浓度限值[45], 是由污染物同特定对象之间的剂量-反应关系确定的. 评价基准的确定是区域重金属污染与风险评价的关键内容, 是制订土壤环境质量标准、开展土壤环境质量评价与管理的基础, 可为地方生态环境监管, 土壤环境损益评估提供借鉴, 为重金属污染环境监测预警、风险控制和生态环境修复提供参考基准和定量标尺[46].

本研究中不同农产品Cd的敏感性分布(SSD)曲线拟合结果见图 7. 若数据点集中在曲线上端, 说明该作物对重金属的吸收富集能力较弱, 敏感度较低, 反之亦然[7,47]. 玉米和小麦的1/BCF值多集中在0~400之间, 水果和蔬菜多集中在0~200, 说明这4种农产品对Cd的吸收富集能力较强, 敏感度较高. 而马铃薯和其他谷物分布相对均匀, 说明对Cd的敏感度低.

图 7 农产品中重金属敏感性分布(SSD)曲线 Fig. 7 Species sensitivity distributions (SSD) of Cd in agricultural products

由公式(11)可以推算出基于保护95%农产品安全的临界值(HC5)和保护5%农产品安全的警戒值(HC95), 也就是HC5和HC95分别为可保护监测区95%和5%农产品的Cd含量低于安全限值. 当土壤重金属含量低于标准中(GB 15618-2018)筛选值时, 农产品无安全风险;当土壤重金属含量高于管制值时, 农产品存在一定的安全风险, 需要采取管控措施. 本研究利用Origin软件拟合的SSD公式见表 5, 进而推导出的玉米、小麦、水果和其他谷物地土壤Cd安全生产临界值分别为0.503、0.336、0.234和0.342 mg·kg-1, 均小于风险筛选值, 因此本研究推导出的玉米、小麦、水果和其他谷物地土壤Cd阈值较标准严格. 水稻和蔬菜地的土壤安全生产临界值为1.057 mg·kg-1和0.740 mg·kg-1, 高于标准中对应的风险筛选值, 因此推导出的水稻和蔬菜地土壤Cd阈值相较标准要宽松. 马铃薯地的土壤安全生产临界值和标准基本一致(0.598 mg·kg-1 ≈ 0.600 mg·kg-1). 以上研究与不同农产品对土壤重金属的富集能力导致土壤健康风险值有差异的结论相一致[22].

表 5 基于log-logistic模型的拟合公式及土壤环境基准1)/mg·kg-1 Table 5 Fitting formula based on log-logistic model and soil environmental benchmark/mg·kg-1

虽然当前现行土壤环境质量标准(GB 15618-2018)作为我国农用地分类分级管理质量标准, 具有十分重要的作用[36], 但本研究结果表明, 利用该标准评价宁夏农产品土壤具有一定的不严谨性, 这与区域环境因子和不同农产品的富集能力等差异有关[48]. 因此, 统一用该标准的筛选值并不完全适用于不同土壤重金属背景值地区和不同富集特征的农产品类型, 难以准确评价土壤重金属的实际暴露水平, 需根据不同农产品中重金属的实际含量进行调整. 因此, 本研究为制订宁夏地方土壤环境质量标准提供了基础资料和参考依据.

3 结论

(1)不同农产品土壤Cd含量平均值均高于宁夏和中国表层土壤背景值, 存在累积富集效应, 但均低于管控值. 土壤Cd空间异质性较强, 污染程度较低, 大部分监测点属于轻度污染和无污染等级, 无高和极高生态风险等级的监测点位.

(2)农产品中其他谷物和蔬菜的Cd含量平均值较高, 马铃薯和水果Cd含量平均值最低, 均低于国家对食品中Cd含量的限量值, 无富集积累效应.

(3)不同农产品产地土壤Cd对儿童的非致癌和致癌风险指数均高于成人, 但对成人和儿童均不存在非致癌和致癌风险. 摄入农产品Cd含量对人体不存在健康风险, 无食源性危害.

(4)本研究推导出的玉米、小麦、水果和其他谷物地土壤Cd阈值较现行土壤环境质量标准(GB 15618-2018)严格, 而水稻和蔬菜地土壤Cd阈值相较标准要宽松.

参考文献
[1] 刘靖宇, 李若怡, 梁永春, 等. 基于特征优选和机器学习的塔里木盆地东缘绿洲土壤镉元素含量预测及健康风险评价[J]. 环境科学, 2024, 45(8): 4802-4811.
Liu J Y, Li R Y, Liang Y C, et al. Soil cadmium prediction and health risk assessment of an oasis on the eastern edge of the tarim basin based on feature optimization and machine learning[J]. Environmental Science, 2024, 45(8): 4802-4811. DOI:10.13227/j.hjkx.202308010
[2] 丁昌峰, 周志高, 王玉荣, 等. 基于生态安全的我国土壤镉环境基准研究[J]. 地学前缘, 2024, 31(2): 130-136.
Ding C F, Zhou Z G, Wang Y R, et al. Environmental criteria for cadmium in soils based on ecological safety considerations in China[J]. Earth Science Frontiers, 2024, 31(2): 130-136.
[3] 王静, 魏恒, 潘波. 中国农田土壤Cd累积分布特征及概率风险评价[J]. 环境科学, 2023, 44(7): 4006-4016.
Wang J, Wei H, Pan B. Accumulation characteristics and probabilistic risk assessment of Cd in agricultural soils across China[J]. Environmental Science, 2023, 44(7): 4006-4016. DOI:10.13227/j.hjkx.202206153
[4] 卢洪斌, 卢少勇, 李响, 等. 长江中游典型湖泊沉积物重金属分布特征、生态风险评估及溯源[J]. 环境科学, 2024, 45(3): 1402-1414.
Lu H B, Lu S Y, Li X, et al. Distribution characteristics, ecological risk assessment, and source tracing of heavy metals in the sediments of typical lakes in the middle reaches of the Yangtze river[J]. Environmental Science, 2024, 45(3): 1402-1414. DOI:10.13227/j.hjkx.202212151
[5] 杨雳, 白宗旭, 薄文浩, 等. 中国农田土壤重金属污染分析与评价[J]. 环境科学, 2024, 45(5): 2913-2925.
Yang L, Bai Z X, Bo W H, et al. Analysis and evaluation of heavy metal pollution in farmland soil in China: a meta-analysis[J]. Environmental Science, 2024, 45(5): 2913-2925. DOI:10.13227/j.hjkx.202305253
[6] 李军, 李旭, 李开明, 等. 黄河流域农田土壤重金属污染特征及其优先控制源分析[J]. 环境科学, 2024, 45(3): 1724-1738.
Li J, Li X, Li K M, et al. Characteristics and identification priority source of heavy metals pollution in farmland soils in the yellow river basin[J]. Environmental Science, 2024, 45(3): 1724-1738. DOI:10.13227/j.hjkx.202305211
[7] 刘海, 魏伟, 宋阳, 等. 皖江经济带耕地重金属健康风险评价及环境基准[J]. 环境科学, 2023, 44(6): 3531-3543.
Liu H, Wei W, Song Y, et al. Health risk assessment and environmental benchmark of heavy metals in cultivated land in wanjiang economic zone[J]. Environmental Science, 2023, 44(6): 3531-3543. DOI:10.13227/j.hjkx.202208037
[8] Wang X D, Zheng W D, Tian W, et al. Groundwater hydrogeochemical characterization and quality assessment based on integrated weight matter-element extension analysis in Ningxia, upper Yellow River, Northwest China[J]. Ecological Indicators, 2022, 135. DOI:10.1016/j.ecolind.2021.108525
[9] Ran J, Wang D J, Wang C, et al. Heavy metal contents, distribution, and prediction in a regional soil–wheat system[J]. Science of the Total Environment, 2016, 544: 422-431.
[10] Lv J S, Liu Y. An integrated approach to identify quantitative sources and hazardous areas of heavy metals in soils[J]. Science of the Total Environment, 2019, 646: 19-28. DOI:10.1016/j.scitotenv.2018.07.257
[11] 鲁如坤. 土壤农业化学分析方法[M]. 北京: 中国农业科技出版社, 2000.
[12] Zhu Y C, Wang L J, Zhao X Y, et al. Accumulation and potential sources of heavy metals in soils of the Hetao area, Inner Mongolia, China[J]. Pedosphere, 2020, 30(2): 244-252.
[13] 雷晓玲, 邱丽娜, 魏泽军, 等. 海绵城市生物滞留带重金属污染特征及风险评价[J]. 环境科学, 2022, 43(8): 4127-4135.
Lei X L, Qiu L N, Wei Z J, et al. Pollution characteristics and risk assessment of heavy metals in the bioretention systems of sponge cities[J]. Environmental Science, 2022, 43(8): 4127-4135. DOI:10.13227/j.hjkx.202110201
[14] 陈文轩, 李茜, 王珍, 等. 中国农田土壤重金属空间分布特征及污染评价[J]. 环境科学, 2020, 41(6): 2822-2833.
Chen W X, Li Q, Wang Z, et al. Spatial distribution characteristics and pollution evaluation of heavy metals in arable land soil of China[J]. Environmental Science, 2020, 41(6): 2822-2833. DOI:10.13227/j.hjkx.201910075
[15] 中国土壤元素背景值[M]. 北京: 中国环境科学出版社, 1990.
[16] 尹芳, 封凯, 尹翠景, 等. 青海典型工业区耕地土壤重金属评价及源解析[J]. 中国环境科学, 2021, 41(11): 5217-5226.
Yin F, Feng K, Yin C J, et al. Evaluation and source analysis of heavy metal in cultivated soil around typical industrial district of Qinghai Province[J]. China Environmental Science, 2021, 41(11): 5217-5226.
[17] Hakanson L. An ecological risk index for aquatic pollution control.a sedimentological approach[J]. Water Research, 1980, 14(8): 975-1001. DOI:10.1016/0043-1354(80)90143-8
[18] Liu X S, Huang D M, Zhu Y M, et al. Bioassessment of marine sediment quality using meiofaunal assemblages in a semi-enclosed bay[J]. Marine Pollution Bulletin, 2015, 100(1): 92-101. DOI:10.1016/j.marpolbul.2015.09.024
[19] 程嘉熠, 王晓萌, 杨正先, 等. 双台子河口沉积物重金属溯源及生态风险评估[J]. 中国环境科学, 2021, 41(3): 1345-1353.
Cheng J Y, Wang X M, Yang Z X, et al. Sources and assessment of heavy metals in sediments in Shuangtaizi estuary[J]. China Environmental Science, 2021, 41(3): 1345-1353.
[20] 马杰, 佘泽蕾, 王胜蓝, 等. 重庆市煤矸山周边农产品镉健康风险评价及土壤环境基准值推导[J]. 环境科学, 2023, 44(9): 5264-5274.
Ma J, She Z L, Wang S L, et al. Health risk assessment and environmental benchmark of cadmium in farmland soils around the gangue heap of coal mine, Chongqing[J]. Environmental Science, 2023, 44(9): 5264-5274. DOI:10.13227/j.hjkx.202210159
[21] 龙新宪, 刘文晶, 仇荣亮. 中国农田土壤重金属污染的人体健康风险评估: 研究进展与展望[J]. 土壤学报, 2024, 61(5).
Long X X, Liu W J, Qiu R L. Research progress and prospects of human health risk assessment of heavy metal pollution in farmland soils of China[J]. Acta Pedologica Sinica, 2024, 61(5). DOI:10.11766/trxb202306130229
[22] 刘娟, 李洋, 张敏, 等. 滇东农田土壤铅污染健康风险评价及基准研究[J]. 农业工程学报, 2021, 37(1): 241-250.
Liu J, Li Y, Zhang M, et al. Health risk assessment and benchmark of lead pollution in agricultural soils in East Yunnan, China[J]. Transactions of the Chinese Society of Agricultural Engineering, 2021, 37(1): 241-250.
[23] 李惠民, 赵刚, 文琛, 等. 赣东北典型重(类)金属污染农田土壤人体健康风险评估[J]. 生态毒理学报, 2023, 18(5): 255-265.
Li H M, Zhao G, Wen C, et al. Health risk assessment of typical heavy metal(metalloid) contaminated farmland soils in part of Northeast Jiangxi[J]. Asian Journal of Ecotoxicology, 2023, 18(5): 255-265.
[24] 陈景辉, 郭毅, 杨博, 等. 省会城市土壤重金属污染水平与健康风险评价[J]. 生态环境学报, 2022, 31(10): 2058-2069.
Chen J H, Guo Y, Yang B, et al. Pollution level of heavy metals in soil and health risk assessment in provincial capital cities of China[J]. Ecology and Environmental Sciences, 2022, 31(10): 2058-2069.
[25] Xiang M T, Li Y, Yang J Y, et al. Heavy metal contamination risk assessment and correlation analysis of heavy metal contents in soil and crops[J]. Environmental Pollution, 2021, 278. DOI:10.1016/j.envpol.2021.116911
[26] 张瑞才. 重金属在不同土壤-作物系统中的迁移与积累[D]. 昆明: 昆明理工大学, 2023.
Zhang R C. Migration and accumulation of heavy metals in different soil-crop systems[D]. Kunming: Kunming University of Science and Technology, 2023.
[27] 中国营养学会. 中国居民膳食指南-2022[M]. 北京: 人民卫生出版社, 2022.
[28] Mao C P, Song Y X, Chen L X, et al. Human health risks of heavy metals in paddy rice based on transfer characteristics of heavy metals from soil to rice[J]. CATENA, 2019, 175: 339-348.
[29] Chen Z K, Muhammad I, Zhang Y X, et al. Transfer of heavy metals in fruits and vegetables grown in greenhouse cultivation systems and their health risks in Northwest China[J]. Science of the Total Environment, 2021, 766. DOI:10.1016/j.scitotenv.2020.142663
[30] 李晓婧, 郑向群, 郑顺安. 叶菜对Cd的富集特征及敏感性分布[J]. 环境科学研究, 2017, 30(5): 720-727.
Li X J, Zheng X Q, Zheng S A. Accumulation and sensitivity distribution of cadmium in leafy vegetables[J]. Research of Environmental Sciences, 2017, 30(5): 720-727.
[31] Tong S M, Yang L S, Gong H Q, et al. Bioaccumulation characteristics, transfer model of heavy metals in soil-crop system and health assessment in plateau region, China[J]. Ecotoxicology and Environmental Safety, 2022, 241. DOI:10.1016/j.ecoenv.2022.113733
[32] 杨利玉, 吴攀, 徐梦琪, 等. 基于土壤标准适宜性评价与环境基准的黔西北高Cd农用地精准管理新思路[J]. 环境工程学报, 2022, 16(12): 4049-4057.
Yang L Y, Wu P, Xu M Q, et al. A novel approach to accurate management of agricultural land with high level Cd in Northwest Guizhou —based on current soil standards suitability evaluation and study of environmental benchmark[J]. Chinese Journal of Environmental Engineering, 2022, 16(12): 4049-4057.
[33] 宋清泉, 徐夕博, 吴泉源, 等. 基于PMF模型的土壤重金属定量源解析及环境风险评价[J]. 湖南师范大学自然科学学报, 2022, 45(1): 76-83.
Song Q Q, Xu X B, Wu Q Y, et al. Quantitative analysis of environmental risk of heavy metal sources in soil based on PMF model[J]. Journal of Natural Science of Hunan Normal University, 2022, 45(1): 76-83.
[34] 戴彬, 吕建树, 战金成, 等. 山东省典型工业城市土壤重金属来源、空间分布及潜在生态风险评价[J]. 环境科学, 2015, 36(2): 507-515.
Dai B, Lü J S, Zhan J C, et al. Assessment of sources, spatial distribution and ecological risk of heavy metals in soils in a typical industry-based city of Shandong province, eastern China[J]. Environmental Science, 2015, 36(2): 507-515. DOI:10.13227/j.hjkX.2015.02.018
[35] 李洋, 张乃明, 魏复盛. 滇东镉高背景区菜地土壤健康风险评价与基准[J]. 中国环境科学, 2020, 40(10): 4522-4530.
Li Y, Zhang N M, Wei F S. A benchmark study on soil health risks of vegetable fields in a high-cadmium background area in eastern Yunnan[J]. China Environmental Science, 2020, 40(10): 4522-4530.
[36] 徐梦琪, 杨文弢, 杨利玉, 等. 黔西北山区耕地重金属健康风险评价及环境基准[J]. 环境科学, 2022, 43(7): 3799-3810.
Xu M Q, Yang W T, Yang L Y, et al. Health risk assessment and environmental benchmark of heavy metals in cultivated land in mountainous area of northwest Guizhou province[J]. Environmental Science, 2022, 43(7): 3799-3810. DOI:10.13227/j.hjkx.202111053
[37] 冯艳红, 王国庆, 张亚, 等. 土壤-蔬菜系统中镉的生物富集效应及土壤阈值研究[J]. 地球与环境, 2019, 47(5): 653-661.
Feng Y H, Wang G Q, Zhang Y, et al. Study on bioaccumulation of cadmium in soil-vegetable system and pollution threshold in soil[J]. Earth and Environment, 2019, 47(5): 653-661.
[38] 迟清华, 鄢明才. 应用地球化学元素丰度数据手册[M]. 北京: 地质出版社, 2007.
Chi Q H, Yan M C. Handbook of elemental abundance for applied geochemistry[M]. Beijing: Geology Press, Beijing, 2007.
[39] GB 15618-2018, 土壤环境质量 农用地土壤污染风险管控标准[S].
GB 15618-2018, Soil environmental quality Risk control standard for soil contamination of agricultural land[S].
[40] GB 2762-2017, 食品安全国家标准 食品中污染物限量[S].
[41] Swartjes F A, Versluijs K W, Otte P F. A tiered approach for the human health risk assessment for consumption of vegetables from with cadmium-contaminated land in urban areas[J]. Environmental Research, 2013, 126: 223-231.
[42] 闵涛, 罗彤, 陈丽丽, 等. 溶解性有机质强化棉花修复镉污染土壤[J]. 环境科学, 2022, 43(3): 1577-1583.
Min T, Luo T, Chen L L, et al. Enhanced remediation of Cd contaminated soil by cotton with DOM[J]. Environmental Science, 2022, 43(3): 1577-1583. DOI:10.13227/j.hjkx.202107253
[43] 陈志良, 黄玲, 周存宇, 等. 广州市蔬菜中重金属污染特征研究与评价[J]. 环境科学, 2017, 38(1): 389-398.
Chen Z L, Huang L, Zhou C Y, et al. Characteristics and evaluation of heavy metal pollution in vegetables in Guangzhou[J]. Environmental Science, 2017, 38(1): 389-398. DOI:10.13227/j.hjkx.201606138
[44] Wei M C, Pan A F, Ma R Y, et al. Distribution characteristics, source analysis and health risk assessment of heavy metals in farmland soil in Shiquan County, Shaanxi Province[J]. Process Safety and Environmental Protection, 2023, 171: 225-237.
[45] 邱荟圆, 李博, 祖艳群. 土壤环境基准的研究和展望[J]. 中国农学通报, 2020, 36(18): 67-72.
Qiu H Y, Li B, Zu Y Q. Soil environmental criteria: research and prospect[J]. Chinese Agricultural Science Bulletin, 2020, 36(18): 67-72.
[46] 李婷, 吴明辉, 王越, 等. 人类扰动对重金属元素的生物地球化学过程的影响与修复研究进展[J]. 生态学报, 2020, 40(13): 4679-4688.
Li T, Wu M H, Wang Y, et al. Advances in research on the effects of human disturbance on biogeochemical processes of heavy metals and remediation[J]. Acta Ecologica Sinica, 2020, 40(13): 4679-4688.
[47] 穆德苗, 孙约兵. 西南地质高背景区蔬菜Pb的安全生产阈值与土地质量类别划分[J]. 环境科学, 2022, 43(2): 965-974.
Mu D M, Sun Y B. Safety production threshold and land quality classification of vegetable Pb in high geological background area of southwest China[J]. Environmental Science, 2022, 43(2): 965-974. DOI:10.13227/j.hjkx.202105032
[48] 冀宏伟, 徐晓航, 许志东, 等. 贵州省岩溶地质高背景地区稻田土壤-水稻系统重金属污染及健康风险评价[J/OL]. 生态学杂志, 1-15. http://kns.cnki.net/kcms/detail/21.1148.Q.20230906.1247.008.html, 2024-06-05.
Ji H W, Xu X H, Xu Z D, et al. Heavy metal pollution and health risk assessment of paddy soil-rice system in Guizhou Karst geological high background area[J/OL]. Chinese Journal of Ecology, 1-15. http://kns.cnki.net/kcms/detail/21.1148.Q.20230906.1247.008.html, 2024-06-05.