广西中东部9县区农田土壤Se输入通量研究

2020-07-26陈雪杨忠芳陈岳龙杨琼王磊韦雪姬

物探与化探 2020年4期

陈雪，杨忠芳，陈岳龙，杨琼，王磊，韦雪姬

(1. 中国地质大学(北京) 地球科学与资源学院，北京 100083； 2.中国人民武装警察部队黄金地质研究所，河北廊坊 065000； 3.广西壮族自治区第四地质队，广西南宁 530031； 4.广西土地质量地球化学评价办公室，广西南宁 530023； 5.广西地球物理勘察院，广西柳州 545005)

0 引言

硒(Se)是与人和动物健康密切相关的微量生命元素[1-2]。人体获取Se的主要途径为食物摄入，土壤是基础媒介。土壤中Se的含量和赋存状态对植物中Se含量水平具有决定性作用，因此土壤中Se的分布特征、赋存状态及来源成为当前研究的热点[3-7]。研究证明，土壤和植物中过量的Se会导致人和动物患“碱性病”和“盲珊症”等慢性中毒症[8]；而土壤和植物中Se含量过低会引起人和动物发生白肌病、克山病、大骨节病等多种缺Se病[9]。

全世界有40多个国家和地区缺Se，中国是世界上缺Se最严重的国家之一，低Se或缺Se地区约占全国总面积的72%，其中30%为严重缺Se地区[10-11]。然而，在这样的大背景下， 1∶250 000 多目标区域地球化学调查成果表明，广西富Se土壤面积多达21 200 km2，为截止到2013年我国圈定出的最大面积连片富Se土壤，较为罕见。自治区地矿局对南宁、北海、贵港3个地区4 350万亩土地的调查发现，富Se土壤达3 181万亩，其中富Se水田592.65万亩，旱地563.1万亩，林地1 553.85万亩，草地269.1万亩[12]。根据黄碧燕等对广西农产品产地土壤Se资源普查及分级评价研究的结果[13]，广西农产品产地土壤28 568个普查点位Se含量范围为0.02×10-6～14.4×10-6，中位值为0.45×10-6，平均含量达0.49×10-6，是全国土壤Se平均含量(0.29×10-6[14])的1.7 倍。广西14 个地级市农产品产地土壤Se含量范围为0.34×10-6～0.57×10-6，均高于全国的平均水平，是全国土壤Se平均含量的1.2～1.9倍。

针对富Se土壤成因，已经有许多专家学者进行了探讨，总的认为，土壤Se有各种来源，如成土母质、大气沉降、灌溉水、污泥、农用石灰、化学肥料、人工加Se等[15]。由于地表土壤中Se与作物中Se含量具有直接的联系，因此研究农田生态系统中的Se输入通量，不但对认识Se的地球化学循环规律具有重要意义，而且对科学利用富Se土地资源具有重要的现实意义。通量通常定义为单位时间内，通过单位面积的物质的量[16]。在本次研究中，输入通量主要考虑面源影响，涉及到大气干湿沉降输入通量(D)、灌溉水输入通量(I)、施肥输入通量(F)。前人针对不同性质农田土壤Se输入通量做了相关研究[17-19]，但没有进行Se与农田土壤的重金属元素相关性方面的研究。土壤Se资源的开发利用必须注意安全问题，生物补Se伴随有重金属潜在危害性， Se在土壤—生物系统迁移的同时，进入土壤的重金属元素也会沿着食物链传递，最终进入人体，危害健康。因此，要明确富Se地区是否存在重金属的污染与危害，并进行相关性分析与评价，这对开发安全高效、生态高值的富Se农产品极其重要[20]。笔者选取了广西中东部9县区，依托《广西中东部地区土壤硒元素和重金属元素地球化学研究》项目收集的广西 “农业办”1∶5 万土地质量地球化学调查资料，开展大气沉降、灌溉和施肥外源输入对农田区土壤Se含量的影响研究，以及Se与As、Cd、Cr6+、Hg、Pb 共5种重金属元素的相关性研究，以期对Se的不同来源能有更清晰的认识，对合理利用研究区农田富Se土壤及控制重金属污染防治具有一定意义。

1 研究区概况

研究区位于广西中东部，区域面积23 322 km2。区内气候为亚热带季风气候和南亚热带季风气候，年平均气温21.5 ℃，年平均降水量1 599.4 mm。地形以孤峰平原地形和残丘平原地形为主。土壤类型以赤红壤为主，也有少量砖红壤和红壤。研究区9个县区包括临桂县、鹿寨县、象州县、兴宾区、平南县、港南区、兴业县、北流市、博白县，其中临桂县属于桂北经济区，鹿寨县、象州县、兴宾区属于桂中经济区，平南县、港南区、兴业县、北流市、博白县属于桂东经济区。本研究区9县区所属5个市的土壤Se平均含量由高到低依次为来宾市、桂林市、柳州市、玉林市、贵港市，其中来宾市、桂林市、柳州市高于或等于广西农产品产地土壤Se平均含量水平，玉林市、贵港市低于广西农产品产地土壤平均Se含量水平[13]。前人对广西部分地区土壤Se地球化学特征及其富集来源进行了探讨。北流地区、兴业县和贵港地区的Se富集物质来源都以内源性为主[21-23]。不在本次研究区范围内的北部湾经济区[24]、全州县[25]、武鸣县[26]、柳城县[27]、横县[28]这几个地区，土壤Se也均主要来源于富 Se岩层。由此可见，成土母质是土壤形成的物质基础，是制约土壤中元素含量的最重要因素[29]。前人研究发现，不同时代成土母质，其Se含量存在一定的差异性，如变质岩>岩浆岩>沉积岩和酸性岩>基性岩[30]。同时，多数研究表明土壤Se含量不仅受成土母岩的影响，还受铁锰氧化物、硫化物成矿元素、有机质及pH值的影响[28]。另外，在不同土地利用方式下，土壤Se经过水热交替、物理化学及生物过程等作用下发生迁移贫化或富集，是影响土壤Se含量的重要因素之一[27]。在不同土地利用方式中，林地和草地中的Se由于树叶、草叶等在土壤表层腐烂堆积而富集于表层；在耕地中，长期耕作过程中由于Se大量支出而导致土壤Se含量减少，形成低Se土壤[31-32]。与此同时，工农业生产及环境污染状况也会对土壤Se的含量产生影响。桂北和桂中经济区(临桂县、鹿寨县、象州县、兴宾区)工业发展处于整个广西领先地位，重点发展以汽车为主的机械工业、电子工业、冶金工业、化学医药工业等。而桂东经济区(平南县、港南区、兴业县、北流市、博白县)主要开发本地有资源优势的轻化工、林化工业，发展技术密集型工业群[33]。

2 样品采集与分析

2.1 布点与采样

按照DZ/T 0289-2015 《区域生态地球化学评价规范》[34]和DZ/T 0295-2016 《土地质量地球化学评价规范》[35]的要求，采集了大气干湿沉降、灌溉水、肥料样品。所有样品都分析Se元素以及As、Cd、Cr6+、Hg、Pb共5种重金属元素含量。图1为大气干湿沉降、肥料、灌溉水采样点位分布图。

图1 广西中东部地区3种外源的采样点位分布Fig.1 Map of the sampling locations of 3 outside-sourced in the middle east of Guangxi

大气干湿沉降样品：从2014～2015年，共布设190个采样点。由于各县区项目承担单位不同，采样接尘缸口径各有差异，主要为 24.5、25和 30 cm的塑料桶，使用前用盐酸溶液(φHCl=10%)浸泡24 h后，再用纯水洗净。每个点放置3个接尘缸，位置为距地面 5～10 m的屋顶开阔平台上，采样口距平台1.0～1.5 m，以避免平台扬尘的影响，最后固定好接尘缸。大气干湿沉降样品每半年接收一次。

灌溉水样品：按照不同灌溉水系布设样点，于灌溉季节在灌溉口处采样，共采集样品500件。采集前用采样点处的水洗涤采样瓶和瓶塞2～3次。根据测试指标不同，每个样点需制备3瓶水样：聚乙烯瓶装原水1瓶，1 500 mL；玻璃瓶装酸化水1瓶，500 mL；聚乙烯瓶装加重铬酸钾水1瓶，500 mL。

肥料样品：共采集复合肥125件，氮肥35件，磷肥10件，钾肥26件，有机肥12件，共计208 件样品，同时记录不同种类肥料施用量，调查肥料的生产地等信息。

2.2 样品分析和质量监控

大气干湿沉降样品、灌溉水样品以及北流、博白、兴业地区的肥料样品由广西壮族自治区地质矿产测试研究中心分析测试，港南、平南、临桂地区的肥料样品由自然资源部南京矿产资源监督检测中心分析测试，鹿寨、象州、兴宾地区的肥料样品分别由自然资源部武汉矿产资源监督检测中心、自然资源部福州矿产资源监督检测中心以及自然资源部合肥矿产资源监督检测中心分析测试。

在分析过程,通过标准参考样、加标回收和室内外重复样、密码样的检验控制分析精密度和准确度。所有测试数据的精密度和准确度合格率都达到100%，符合DD2005-03 《生态地球化学评价样品分析技术要求(试行) 》[36](表1)。

表1 样品各元素分析方法及检出限

3 结果与讨论

大气干湿沉降、施肥、灌溉水等外源因素，不仅为土壤带来有益元素Se，还会带来As、Cd、Cr6+、Hg、Pb 等有害重金属元素，对土壤造成污染。土壤中的重金属，进入植物体内，能诱导其体内产生某些对酶和代谢具有毒害作用和不利影响的物质，如H2O2、C2H2等，会对植物产生一定的毒害作用，引起株高、主根长度、叶面积等一系列生理特征的改变[37]。土壤中的重金属通过食物链进入人体，直接对人体健康造成威胁[38-39]。因此笔者将Se及As、Cd、Cr6+、Hg、Pb共5种重金属元素输入通量进行对比分析，以期对有益元素及有害元素的输入途径、输入量大小等获得更加清晰地认识，有助于富Se土壤的有效利用及预防和控制重金属污染。

3.1 大气干湿沉降年输入通量

大气干湿沉降是Se输入土壤的重要途径。Pa-cyna[40]等对全球人为源每年向大气的平均释Se量定，肥料中Cr6+在合肥采用ICP-AES测定，肥料和干沉降中Cr6+在南宁、福建、武汉、南京采用ICP-MS测定。

进行过估算；Wang等[41]根据降雨Se浓度估算中国每年Se沉降为600～1 200 t。

在本研究中，大气沉降样品分为溶液和沉淀两部分，分别测得Se及 As、Cd、Cr6+、Hg、Pb浓度，另根据实际样品溶液的总体积和沉淀物的总质量分别求得各元素的沉降量，再由沉降收集桶的接收面积计算得到每公顷土壤上各元素的年沉降通量(D)：

D=(CsV+CiM)/(S×100)，

式中，D为年沉降通量(g/(hm2·a))，S为接尘面积，即采样桶桶口面积(采样桶口径不同，其接尘面积也不同，分别为0.047、0.049、0.07 m2)，Cs为溶液样品中某元素的浓度(μg/L)，V为年溶液总体积(L/a)，Ci为沉淀样品中某元素的质量分数(μg/g)，M为年沉淀总质量(g/a) 。

根据190个大气干湿沉降接收点全年实测数据，计算得出广西中东部9县区Se及 As、Cd、Cr6+、Hg、Pb的大气干湿沉降输入通量(表2)。所有样点中， Se年沉降通量最大值为66.10 g/(hm2·a)，最小值为0.93 g/(hm2·a)，平均值为6.36 g/(hm2·a)。Se在各地区的沉降通量由大到小依次为临桂>鹿寨>兴宾>平南>北流>博白>港南>兴业>象州(图2)。Se年沉降通量表现为明显的地区差异，桂北和桂中经济区(临桂县、鹿寨县、兴宾区)高于桂东经济区(港南区、平南县、兴业县、北流市、博白县)，而桂中经济区的象州县最低。

表2 研究区Se及重金属元素3外源输入通量 g/(hm2·a)

Se与5 种重金属元素的沉降通量最大值基本都出现在鹿寨县和临桂区，其中Se与As的空间分布趋势最为相似，说明两种元素的来源有一定共性，可能与煤的燃烧、颜料的生产、冶金以及化工生产有关。桂林和柳州是广西传统的工业区，柳州是我国工业重镇，制造中心；桂林的电子工业占广西的70%以上，食品工业、制药工业以及客车制造业位居广西前列。因此这两个地区的大气污染较其他区域严重，大气中的Se以及重金属元素沉降通量较大。而Cd的沉降最大通量出现在象州县，这可能是因为象州县拥有广西面积最大的工业园，园区内产生Cd污染的企业类型较多，因此造成大气污染物中Cd含量高于其他地区。Se与As、Cr6+的沉降量低值出现在象州县，这可能是因为象州县的化工类产业最少。而Cd、Hg、Pb的沉降量低值出现在桂东经济区，这可能是因为该区的支柱产业为轻化工和林化工，这些企业产生的Cd、Hg、Pb污染相对较少。玉林市(北流、博白、兴业)沉降通量整体最低，北流、博白、兴业三地的差别很小，这与当地支柱产业类型以及污染治理程度紧密相关。总体来看，桂北和桂中经济区重金属沉降通量大于桂东经济区，5种重金属元素在各地区沉降输入通量总和由大到小依次为临桂>兴宾>鹿寨>象州>北流>平南>兴业>博白>港南(图2)，这与干湿降尘中Se的沉降通量分布规律基本一致。由此可见，工业较为发达的桂北和桂中经济区污染排放更高，因此重金属元素和Se沉降通量均较高；而发展技术密集型工业群的桂东经济区污染较小，因此大气质量相对较好，重金属元素和Se排放量相对也较少。

图2 广西中东部9县区农田土壤Se及重金属大气干湿沉降平均输入通量Fig.2 The average deposition flux of Se and heavy metals in 9 counties of the middle east of Guangxi

3.2 灌溉水年输入通量

农田灌溉对于提高农作物的产量具有重要的意义，同时也是Se及重金属元素进入土壤的途径之一。根据《广西壮族自治区2014年水资源公报》和《广西壮族自治区2015年水资源公报》公布的数据，2014年全省农田实灌亩均用水量916 m3，有效利用系数为0.446，折合6 128 m3/(hm2·a)，2015年全省农田实灌亩均用水量873 m3，有效利用系数为0.465，折合6 089 m3/(hm2·a)，计算每公顷每年Se及重金属元素的输入通量(I)，即：

I=CwV/1000，

式中，I为年输入通量(g/(hm2·a))，Cw为样品中某元素的浓度(μg/L)，V为年灌溉水量(m3/(hm2·a))。

研究区灌溉水Se输入通量平均值为2.45 g/(hm2·a)，9个县区输入通量呈现明显的地区差异，位于河流下游的北流、博白、兴业的输入通量要高于位于河流上游的其他6个县区(图3)。5种重金属元素的输入通量在县区间的差异和Se的特征基本一致。

图3 广西中东部9县区农田土壤Se及重金属灌溉水平均输入通量Fig.3 The average irrigation flux of Se and heavy metals in 9 counties of the middle east of Guangxi

3.3 施肥年输入通量

肥料是植物的粮食，农作物的高产离不开肥料的施用，不论发达国家还是发展中国家，施肥都是粮食增产最有效的措施[42-43]，但其在补充作物生长必需营养元素的同时，不可避免地会将一些有毒有害物质(如重金属)带入到土壤中[44-45]，一旦这些有害物质积累到一定程度，势必会影响植物的营养生长和生殖生长，从而影响产品品质[46-47]。土壤中Se及重金属元素的累积都是漫长的过程，了解长期施肥对土壤和作物产品Se及重金属含量的影响，对肥料产品的安全生产和科学施用，合理利用富Se土壤，防止农产品污染，保证农产品安全和农业可持续发展都具有十分重要的意义。

研究区复合肥、氮肥、磷肥、钾肥、有机肥中Se及 As、Cd、Cr6+、Hg、Pb重金属元素含量统计结果见表3。肥料在带入Se的同时，也造成了重金属的积累。5种主要肥料中磷肥的Se及 As、Cd、Cr6+、Hg、Pb含量明显最高，其次为有机肥。磷肥的主要原料是磷矿石，磷矿石加工过程中会有部分重金属随原料进入磷肥[48-49]。欧洲12个国家196 种磷肥重金属平均质量分数为：w(Ni)=14.8×10-6、w(Cd)=7.4×10-6、w(Zn)=166×10-6、w(Pb)=2.9×10-6、w(As) =7.6×10-6和w(Cr)=89.5 ×10-6，进入土壤的量和含磷水平紧密相关[50]。大量研究表明，施用有机肥影响作物体内重金属含量。有机肥对作物可食部位重金属含量的影响，与有机肥种类、用量、土壤类型和pH以及作物种类等有很大关系[51-53]。氮肥和钾肥对土壤重金属累积的影响不大，主要是因为氮、钾肥含有的重金属较少[54]。

表3 广西中东部9县区肥料中Se及重金属元素含量平均值 10-6

因为各种肥料的元素含量差异较大，且年施用量各不相同，所以不能把不同种类肥料的元素含量简单相加。本文采用的方法是，根据不同种类肥料年使用量，分别计算各种肥料的输入量，然后相加。肥料输入通量(F) 公式计算为：

式中，F为某点位的年输入通量(g/hm2·a)，i={复合肥，氮肥，磷肥，钾肥，有机肥}，Ci为某元素在肥料i中的含量(10-6),Mi为某点位肥料i的年施用量(g/(hm2·a) )。

根据不同种类肥料的施肥量计算出研究区各地肥料中Se及 As、Cd、Cr6+、Hg、Pb的年输入通量(图4) 。Se的施肥平均输入通量为0.15 g/(hm2·a)，输入通量呈现较为明显的地区差异，兴业、港南、临桂地区肥料中Se的年输入量普遍较高(平均通量为0.22～0.50 g/(hm2·a))。As、Cd、Cr6+、Hg、Pb的年输入通量也表现出明显的地区差异，兴业、港南、鹿寨、临桂地区肥料中重金属总量输入通量较高(平均通量为20.47～75.87 g/(hm2·a))，这取决于肥料中元素的含量以及施肥量。Se的输入通量地区差异与5种重金属元素总和的变化趋势基本一致。

图4 广西中东部9县区农田土壤Se及重金属施肥平均输入通量Fig.4 The average fertilization flux of Se and heavy metals in 9 counties of the middle east of Guangxi

3.4 年输入总量

研究区土壤中 Se 的外源输入总量(IP)由 3 种途径的通量相加而得，即：

IP=D+I+F。

研究区不同地区Se输入通量存在明显差异(图5)，9个县区按通量大小可分为 3 组，临桂区、鹿寨县的输入通量最高(>16.70 g/(hm2·a))；博白县、北流市、兴宾区输入通量大小分布在平均值(9.11 g/(hm2·a))附近；兴业县、港南区、平南县、象州县输入通量低于平均值。综合分析本文的3种输入途径，研究区大气干湿沉降、灌溉水和施肥对农田土壤Se的贡献率分别为70.93%、27.35%和1.72%，Se的主要来源是大气干湿沉降，大部分地区的大气干湿沉降输入比例在50%以上，灌溉水的输入量次之，而施肥所占输入比例最小(图6)。因此Se输入总量的地区差异也主要是由大气干湿沉降的差异所致，经济和工业更发达的临桂区、鹿寨县Se输入总量高于轻工业为主的桂东经济区，这与重金属元素输入通量特征一致(大气干湿沉降、灌溉水和施肥3种输入途径对5种重金属输入通量总和的贡献率分别为74.21%、15.67%、10.12%)，反映了研究区冶金、机械、化工等各项活动与大气中Se和As、Cd、Cr6+、Hg、Pb重金属元素增加有关。5种重金属元素，除As是灌溉水输入通量占主导地位，其他元素都是大气干湿沉降的输入通量占主要比例，这也和Se一致(图7)。前人对江西鄱阳湖流域、黑龙江松嫩平原以及四川绵阳的农田土壤Se的地球化学循环进行了详细研究[17-19]，本研究区Se输入通量与这些地区进行比较，大气干湿沉降是所有地区的主要输入途径。