大尺度生态干扰风险评估技术方法及应用研究

为实现大尺度生态干扰高风险区域的快速、有效识别,从生态易损性、干扰易达性和资源易引性三个方面选取了9个指标,构建了生态干扰风险评估指标体系和评估模型,形成了完整的大尺度生态干扰评估技术与方法.利用该技术方法,本文评估分析了全国国土空间生态干扰风险状况、空间分布格局及其成因,并结合2017~2019年全国自然保护区人类活动监测数据和《全国主体功能区规划》中生态脆弱性评价结果对生态干扰风险评估结果进行了精度分析.结果显示：超过90%的自然保护区人类活动都集中在评估得到的中高风险区域;全国生态干扰风险主要以中风险为主、低风险次之、再次是较低风险、较高和高风险区域面积最少;其中较高和高风险区域主要分布在我国中西部和东北部的秦岭、祁连山、三江源和内蒙古草原等区域,生态干扰风险空间分布格局与生态系统结构与功能状况、地形与交通条件、资源潜力等因素具有较强相关性.     

宝鸡市主要农作物冰雹灾害风险评估

冰雹灾害是宝鸡市面临的主要农业气象灾害之一,采用宝鸡市2016年的Landsat TM遥感影像,经过遥感解译得到宝鸡市主要农作物的空间分布图。根据宝鸡市1961 — 2016年55 a的降雹记录和防雹炮点分布资料,计算得到区域降雹强度、降雹频次和防雹能力值,以灾损数据确定主要农作物的承灾体敏感性指数。基于致灾因子危险性、承灾体暴露度及脆弱性三要素,采用综合评价法和自然间断分类法得出宝鸡市主要农作物的冰雹灾害风险分布图。结果表明:宝鸡市主要农作物的冰雹灾害风险等级以中级风险为主,占全部风险等级的43.5%,主要分布于千阳、凤县、岐山、扶风等地;次高风险区占25.2%,主要分布在陈仓及扶风;次低风险区占18.7%,分布于金台、眉县及太白零星散区;高风险区占9.8%,全部集中于陇县和麟游;低风险区占比最少,为2.8%,零星分布于凤县和扶风。应加强对冰雹灾害高风险区的防爆能力建设,调整防雹点分布和作物种植。     

南阳盆地东部山区土壤重金属分布特征及生态风险评价

通过分析南阳盆地东部山区表层土壤（0~20 cm）样品的As、Hg、Cd、Cr、Pb、Zn、Ni和Cu等8种重金属含量和pH值,对土壤重金属的空间分布特征、污染程度和生态风险进行研究,并对重金属的来源进行分析.结果表明,研究区内土壤Cd、Cu、Cr、Ni、Pb、Zn和As含量相比于土壤污染风险筛选值存在不同程度的超标.土壤重金属空间分布呈面状和岛状分布.含量高值区主要分布在研究区南部,且与矿区分布相吻合.除了少数样点的土壤以外,大部分样点未受到污染,Cd的污染最为严重.As、Cr、Pb、Zn、Ni和Cu在几乎全部地区的潜在生态风险系数为轻微风险等级,Cd和Hg的潜在生态风险系数为中等风险等级占比最高.研究区综合土壤重金属潜在很强和强生态风险指数占比较高,分别为58.93%和37.66%.土壤Hg、Cd和Pb主要受到矿产开采等人类活动的影响,As主要受到研究区地质背景的影响,Zn、Ni、Cr和Cu同时受到人类活动和地质背景的影响.研究区矿产资源的开采是造成土壤污染及生态风险的主要原因.     

博斯腾湖流域景观生态风险评价与时空变化

以博斯腾湖流域为研究区,利用eCognition软件对博斯腾湖流域1990年、2000年、2010年和2015年4期Landsat遥感影像分类。基于土地利用类型的生态脆弱度和景观干扰度指数构建景观生态风险评价模型,并利用空间自相关分析方法,开展博斯腾湖流域景观生态风险评价的时空变化研究。结果表明:(1)1990-2015年,博斯腾湖流域景观格局变化明显,主要特征为水域景观面积减少了64.33%,景观优势度降低;林地、耕地和城镇用地增幅显著,三者分别增加了47.96%、44.77%和85.15%;草地仍是流域主导景观类型。流域整体斑块数量减少,景观破碎化程度减弱。(2)低风险区和较低风险区主要分布在流域北部,土地利用以草地为主;高风险区主要分布在流域东南部,土地利用以未利用地为主;较高风险区主要围绕高风险区分布;中等风险区分布比较零散。(3)1990-2015年间,高风险区面积增加了386 263 hm2,低风险区面积减少了217 216 hm2,研究区总体生态风险等级增加,博斯腾湖流域生态质量存在进一步退化的趋势。     

基于乡镇尺度的地质高背景区耕地土壤重金属来源分析与风险评价

西南地区由于广泛发育碳酸盐岩导致土壤重金属具有天然的高背景属性,在重金属地质高背景区开展生态风险评价和识别重金属潜在来源对土壤污染防治、农产品安全保障具有重要意义.基于乡镇尺度,在重庆市奉节县青龙镇采集耕地表层土壤样品307件,分析了土壤As、Cd、Cr、Cu、Hg、Ni、Pb和Zn重金属含量及pH,采用地统计学,绝对因子分析-多元线性回归(APCS-MLR)受体模型探讨研究区域内8种重金属污染空间分布特征及来源,利用内梅罗指数(P_综)和潜在生态风险指数(RI)开展土壤风险评价.结果表明,耕地土壤中As、Cd、Cr、Cu、Hg、Ni和Zn的平均值和中位数超过重庆市土壤背景值,其中Cd、As累积效应较为明显.与农用地土壤筛选值相比,Cd超标率最高,达52.12%.土壤中重金属含量空间分布格局呈现出南部整体高,高值区呈条带状分布;北部整体低,高值区呈点状分布.APCS-MLR模型表明研究区耕地土壤重金属Cr、Zn、Ni和Cu以自然来源为主,主要受地质背景控制,贡献率分别为86.62%、85.46%、76.44%和64.34%;As、Pb和Hg的来源以工业源为主,受煤矿...     

新建铅蓄电池集聚区对周边土壤环境的影响:基于重金属空间特征

为探明浙北某乡镇经提升改造后一新建铅蓄电池集聚区运行7a后是否对周边土壤环境存在影响,采集该铅蓄电池集聚区周边表层土壤(0～20 cm) 76份,测定了土壤中汞(Hg)、砷(As)、铜(Cu)、锌(Zn)、铅(Pb)、镉(Cd)、镍(Ni)和铬(Cr)这8种重金属含量,并基于集聚区内、距集聚区边界50、450和850 m的空间距离,采用单因子指数法、内梅罗综合污染指数法及潜在生态风险指数法对土壤环境质量进行了评价,然后利用地统计方法分析了重金属空间分布特征,并结合相关性分析确定了对土壤环境造成影响的重金属的来源.结果表明,与当地平均背景值相比,8种重金属元素中Hg、Zn和Pb在所有空间尺度下的平均含量均高于其对应背景值,对Cd而言,除集聚区内,其余空间尺度下的Cd平均含量均大于其背景值,As只有距边界50m处的平均含量大于其背景值,而其他元素在所有空间距离下的平均含量均低于其对应背景值,其中Hg和Cd存在高度空间变异,而其他元素含量空间变异不明显,说明区域活动的影响主要集中在Hg和Cd上,且两者的含量随集聚区距离外延而增加.出现超出农用地土壤污染风险筛选值点位的元素主要为Hg和Cd,其主要分布在集聚区外450 m和850 m处,其中Hg在对应距离下超出风险筛选值的点位占33. 33%和38. 89%,Cd分别占27. 78%和55. 56%,且两者的空间分布特征与其含量一致,而其他元素中仅有Zn和Pb存在零星点位超出风险筛选值,且总体上无明显空间特征.由8种重金属元素对土壤的综合污染风险分析可知,Cd是造成土壤综合污染风险的主要来源,由于其贡献使850 m处土壤处于警戒状态(贡献率为36. 73%).土壤的生态风险主要出现在集聚区外450 m和850 m处,处于中等生态风险水平,其中生态风险主要来自Hg和Cd,Hg在对应距离下的贡献率分别为46. 30%和39. 37%,Cd分别为38. 98%和49. 30%,说明区域活动使Hg和Cd成为影响研究区土壤质量的主要元素.经地统计和多元统计分析表明,Hg和Cd含量呈现出在当地主风向(东北-西南)轴上由集聚区外围向内扩散的特征,且两者的主要来源为集聚区外围企业的燃煤活动.因此,新建的铅蓄电池集聚区运行7a后并未对集聚区及周边土壤重金属的集聚造成明显影响.     

基于PMF模型的兰州耕地土壤重金属来源解析

为了分析Pb、Cd、Cr、Hg、As、Cu、Zn、Ni等8种重金属的污染状况、空间分布及污染来源.2019年9月共采集382个兰州耕地表层土壤样品,并测定其重金属含量.结果表明：①兰州市耕地土壤重金属Pb、Cd、Cr、Hg、As、Cu、Zn、Ni的平均含量分别为20.96,0.24,0.05,6.40,8.36,48.57,101.92,34.34mg/kg,其中Hg的均值超过国家土壤污染风险筛选值,Cd、Hg、Cu、Zn、Ni的均值超过兰州市土壤背景值.内罗梅综合指数表明研究区耕地52.62%样点呈现轻度污染,37.70%样点呈现中度污染,3.66%样点呈现重度污染,无重金属污染的样点占6.02%.潜在生态风险指数表明13.09%的样点属于轻微生态风险,68.85%的样点属于中等生态风险,17.54%的样点属于强生态风险,0.52%的样点属于很强生态风险.②基于PMF模型可知8种重金属有4种主要来源,其中Pb、Cd、Zn、Ni以交通运输源为主,贡献率分别为52.4%、50.7%、56.7%和50.2%;Hg以农业活动和医疗设备源为主,贡献率为83.7%;Cr、Cu和Ni以自然源为主,贡献率分别为81.2%和61.3%和49.8%;As以工业活动源为主,贡献率为68.3%.③从空间分布来看,研究区内除Cr外,Pb、Cd、As、Cu、Zn、Ni在城关区和七里河区都有高值区出现,其中Zn和Ni的高值范围较广,可能与人为影响有关;Hg在城关区、安宁区、七里河区以外的研究区内均有高值出现.     

基于乡镇尺度的西南重金属高背景区土壤重金属生态风险评价

西南地区土壤重金属具有天然的高背景属性,在重金属高背景区开展生态风险评价和识别重金属潜在来源具有重要意义.以往工作主要以区域尺度的调查工作为主,调查精度难以满足空间规划和自然资源管理的需求,迫切需要在乡镇尺度上进行生态风险评价.云南省宣威市热水镇是典型的碳酸盐岩覆盖区,重金属来源复杂,潜在生态风险较高.本文在热水镇采集土壤表层样品(0～20 cm)1 092件,分析表层土壤中8种重金属(Cd、 Cr、 As、 Hg、 Pb、 Cu、 Zn和Ni)含量,采用统计学、地理信息系统、地累积指数、生态风险指数和正定矩阵因子分析模型(PMF)等方法开展生态风险评价和重金属源解析.研究发现,土壤中重金属Cd、 Cr、 Cu、 Hg、 Ni和Zn含量平均值和中位数均超过云南省土壤背景值,Cd、 Cr、 Cu和Ni平均含量均超过《土壤环境质量农用地土壤污染风险管控标准》(GB 15618-2018)规定的筛选值.地累积指数结果表明,研究区表层土壤污染最严重的是Cu,其次是Cr和Cd.土壤重金属Cr、 As、 Pb、 Cu、 Zn和Ni多以残渣态存在,多来自于地质背景,生物有效性较低,Hg的潜在有效组分较高,但是Hg全量较低,污染风险较小.Cd生物有效性比例较高,易于进入土壤溶液并被农作物吸收,是研究区污染风险最高的重金属元素.土壤重金属潜在生态风险指数统计结果显示,生态风险较低、中等生态危害和强生态危害比例分别为44.23%、 54.40%和1.37%,无很强生态危害样品,贡献率最高的元素为Cd和Hg,贡献率分别为34.62%和34.33%.PMF结果显示研究区重金属来源包括人类日常活动来源、自然来源、煤矿开采及交通排放引起的综合污染源和农业生产来源,各种来源的贡献率分别占9.29%、 53.67%、 11.23%和25.81%.PMF可以有效识别和量化污染物来源,为政府部门进行土地规划提供重要的参考依据.     

天山北坡经济带土壤重金属来源及污染评价

为定量识别与评价天山北坡经济带中奇台、吉木萨尔、阜康等地区土壤重金属来源与生态风险,对该区域171个表层土壤中Zn、Cu、Cr、Pb、Hg、As和Cd7种常见重金属的含量进行测定.运用统计学方法、主成分分析、正定矩阵因子分解（PMF）与潜在生态风险指数进行重金属污染程度评价以及来源分析.结果表明,研究区土壤Zn、Cu与Cr之外,Pb、Hg、As和Cd分别超过了新疆土壤背景值4.1、2.0、8.0和48.0倍;与国家土壤重金属风险筛选值相比,Zn、Cr和Hg浓度在安全范围内,Cu与Pb有少部分样点污染较严重,超出筛选值,As与Cd平均值分别超出筛选值的3.09倍与19.17倍.污染来源分析结果显示研究区土壤中的重金属元素主要来自于燃煤源、交通运输、大气降尘、农工业排放和自然因素.生态风险评价结果分析表明,Zn、Cu、Cr和Pb处于轻微风险状态;Cd处于极高风险水平;73.68%的Hg处于中等风险,18.71%处于高风险;43.86%的As生态风险处于中等风险,51.46%处于高风险.综合潜在生态指数介于472~2575.69,Cd对综合潜在生态指数贡献率达到了89.24%,其次是As与Hg,表明研究区Cd对土壤生态环境危害很大,As与Hg也需要特别引起重视.     

三门峡某铅厂遗留场地土壤重金属空间分布特征及来源解析

我国各大城市目前正在进行大规模的工业企业搬迁,并产生了众多污染搬迁遗留场地,其中重金属污染尤为严峻.为了分析三门峡某铅厂遗留场地土壤重金属的污染状况、空间分布和污染来源,采用地统计学分析场地土壤重金属的空间变异规律和分布特征,并利用PMF模型解析场地土壤重金属的主要来源.结果表明,土壤中As、Cd、Cu、Pb、Hg和Ni的平均值远超过河南省土壤环境背景值,As、Cd、Pb和Hg含量超过场地土壤污染风险筛选值,As、Pb和Hg含量超过场地土壤污染风险管控值.Cr的高值区位于废渣堆场北侧,Ni和Cd高值区位于废渣堆场北侧和场地南侧,As的高值区位于废渣堆场南侧与生活区之间,Cu和Pb高值区较为分散,主要集中于中部的原料堆场和炼炉区,Ni和Cd、Cu和Pb具有相同的空间分布特征.基于PMF模型可知,7种重金属有3种主要来源,Cd以废渣堆积源为主,贡献率为87.60%;Cu、Pb和Hg以土壤母质源为主,贡献率分别为92.50%、75.20%和95.40%;Cr、Ni和As以原料粉尘废气源为主,贡献率分别为80.80%、83.30%和62.00%.     

青藏高原生态风险及区域分异

全球变化背景下,青藏高原生态系统受到自然、人为双重影响和威胁,生态风险日益加剧。针对生态风险源、脆弱性以及风险管理能力选取30个评估指标,利用生态风险评估优化模型,综合评估了青藏高原的生态风险,并得出如下结论：青藏高原生态风险总体处于较低水平,以低和极低生态风险为主,共占研究区面积的55.84%;极高风险主要分布在北部高山和极高山地区,中等风险主要分布在高原北部以及高原的西部和西南部地区,中、高生态风险在空间分布上形成一个“C”字形结构;青藏高原生态风险整体受自然主导因子控制,人为对生态环境的影响不容忽视,协调和降低青藏高原人类活动区域人类对生态环境的影响,是今后规避生态风险的重要途径。     

中国A级旅游景点空间分布特征与可达性

以A级景点为例,基于GIS 技术,通过最近邻指数、K指数、热点聚类等方法分析中国2 424 个A级旅游景点的空间分布特征;运用栅格成本加权距离算法,计算中国A级旅游景点的空间可达性,并测算县域单元的整体可达性;利用空间关联方法分析了县域可达性的空间差异。结果表明：中国A级景点的空间分布总体上呈现出聚集分布的特点,所有景点的平均可达时间为125.88 min,60%的景点可达时间在90 min 以内;可达性在全国的分布差异十分明显,可达性空间分布具有明显的交通指向性。景点的县域可达性呈强集聚格局,可达性热点区域的分布呈现明显的热点—次热点—次冷点—冷点自东向西带状分布的格局。     

黄土丘陵沟壑区生态风险时空动态及其风险分区—以陕西省米脂县为例

目前,黄土丘陵沟壑区等生态脆弱区生态风险评估已成为地理学与生态学应对生态系统管理的研究热点之一。以黄土丘陵沟壑区米脂县为研究区,构建“风险概率—敏感性—损失度”（PSI）的三维评价框架,并以子流域为评价单元进行数据整合,分析了米脂县2009-2015年准则层与综合生态风险的时空分异及其重心转移,并基于风险主导因子给出米脂县风险防范分区及降险对策。结果如下：（1）2009-2015年风险概率分别为49.93%、52.92%,有上升趋势;生境敏感性分别为0.61、0.60,下降了1.6%,生境质量趋好;损失度分别为0.42、0.46,损失度增加。（2）生态风险呈现中间高南北低的空间分布;研究期间生态风险值分别为0.14、0.15,风险有所升高;风险重心向西南转移跃入银州川道且风险演化主体方向为西北—东南走向。（3）风险预警区、生态恢复区、预警恢复兼顾区、自然调控区面积占比分别为7.53%、6.57%、23.86%、62.04%。基于风险主导因子的风险防范分区可有效进行风险消解,促进区域生境的可持续。     

旅游景区(点)系统空间结构随机聚集分形研究——以南京市旅游景区(点)系统为例

根据旅游景区(点)系统自组织演化的特征,以南京市旅游景区(点)系统为例,采用随机聚集分形的方法对其空间结构进行了研究。认为景区(点)系统的演化类似于有限扩散集团凝聚模型,得出景区(点)系统空间结构是分形的,具有比较明显的随机聚集分形结构,在不同的地理区域上表现为不同的系统演化时段,空间结构上的向心性作用也有强弱之分,从系统的聚集中心向四周是景点密度衰减的。指出随机聚集维数的测算与分析是景区(点)系统空间结构自组织优化趋势的判据,同时提供了优化的途径。     

中国红色旅游经典景区网络关注度时空演变及影响因素

借助“百度指数”,获取2011—2019年中国31个省（市、自治区）对红色旅游经典景区的网络关注度,结合区域差异分析法、景区偏好系数分析红色旅游经典景区网络关注度的时空演变特征,利用面板数据回归模型和地理探测器揭示其影响机理。研究发现：（1）从时序上看,2011—2019年红色旅游经典景区的网络关注度持续上升,具有明显的季节性差异,呈“M”型变化形态。（2）从区域上看,东—中—西部具有明显的梯度递减特征,区域差异呈东—西—中部依次递减态势,空间结构比较均衡,东部地区多位于热点区,中部地区多处于一般区,西部地区则属于冷点区,各省（市、自治区）存在一定季节性差异。（3）景区角度方面,景区的月际变化分为单峰型、双峰型和多峰型三种,单峰型季节差异显著,双峰型和多峰型的季节差异较小,景区偏好类型呈“金字塔”结构,具有“多中心”分布特征,整体格局相对稳定。（4）影响因素方面,客源地网络发展水平的贡献>人口规模>教育水平>经济发展水平,而旅游恩格尔系数的影响不显著,目的地的星级酒店数量、旅游化水平、市外交通可达性是重要影响因素,信息化水平、空气质量、旅行社及公共管理能力是主要影响因素,其他因素影响力较小。     

基于双模型的星云湖沉积物重金属溯源及风险评价

为了解星云湖沉积物中重金属污染现状,并为湖区污染防治提供理论依据,通过分析湖区表层沉积物中Cu、Zn、Pb、Ni、Cr、Cd、As、Hg等重金属含量,研究了其风险特征及污染来源组成.结果表明:①仅Ni、Cu含量尚未超过星云湖的背景值,其余元素均出现累积,Hg含量是星云湖背景值的2.13倍,潜在生态风险最高.②星云湖沉积物中重金属的空间分布与入湖河流和湖盆地势相关,表现出由河流输送向湖心扩散的趋势,高浓度区均分布在湖心深水区和特定河口位置,Hg的分布异于其他元素,呈西北湾>南部>中部的特征.③西北湖湾潜在生态风险最高,Cd、As是主要贡献元素.④星云湖沉积物中重金属主要来自自然源（占比为43.32%）、工业源（占比为25.89%）、农业源（占比为20.52%）和交通源（占比为10.27%）.研究显示,星云湖自然源重金属污染贡献最大,人为源相对较低,但湖区重金属污染处于较高水平.星云湖的重金属污染生态风险可能会持续加剧,其中Hg是污染防控的重点.      

基于“生态水位-水质-水源地”协同作用的地下水环境风险评价方法研究

地下水生态水位、地下水水质和地下水水源地是地下水环境的关键敏感属性,而现有关于地下水环境保护的研究多侧重于污染风险评价,常忽视地下水生态水位作为环境要素综合风险的重要性.为突破以往地下水风险评价中以污染要素为主的现状,在区域尺度针对地下水污染及地下水水位变化导致地下水系统生态服务功能失衡等问题,提出以地下水生态水位、地下水水质和地下水水源地作为风险受体,综合研究地下水系统对地下水风险源的暴露途径及响应关系,采用地下水防污性能指数指征地下水环境的空间差异性,构建了基于“生态水位-水质-水源地”协同作用的GERRM模型（地下水环境相对风险模型,groundwater environmental relative risk model）,定量描述地下水污染和地下水水位突变耦合的地下水环境风险,并将建立的方法在下辽河平原进行案例研究.结果表明:①Ⅳ、Ⅴ级风险区主要位于化工企业、危险废物填埋场周围一定范围区域及沈阳市地下水水源地保护区和生态水位敏感性较高区,面积为2 107.33 km2,占总面积的8.93%.②Ⅰ、Ⅱ级风险区主要为农田种植区、林地种植区和农业城镇建设区,面积为17 704.51 km2,占总面积的75.01%.研究显示,GERRM模型适用于区域地下水环境风险评价,下辽河平原区化工企业、危废填埋场周围一定范围区域以及沈阳地下水水源地保护区相对风险最高,需采取相应的管理保护措施.      

中国红色旅游经典景点空间分布格局及其影响因素异质性

红色旅游是传承革命文化、推动爱国主义教育常态化和大众化、促进革命老区经济发展的重要方式。通过构建理论分析框架系统,综合运用多种空间分析方法与回归模型,对中国公布的300处红色旅游经典景区共763个景点的空间分布格局及影响因素异质性进行研究,并提出发展模式及对策。研究表明：（1）经典景点空间分布受资源、环境、社会和经济四大子系统综合影响;（2）核心区空间分布较其他旅游景区分散,形成7个一级核心区、7个及两大连片二级核心区、低密度与较低密度区等,空间关联性差异较大,关联效应伴随空间尺度缩小而减弱;（3）空间分布影响因素异质性明显,影响程度为红色旅游资源>海拔>政策力度>坡度>交通设施>市场需求>资源联动>消费水平;（4）提出采取网络、节点、廊道、共生—联动、跨区共促等发展模式及对策。研究可为红色旅游发展与管理、跨区联合等提供科学决策依据。     

宁夏引黄灌区农田土壤重金属生态风险评价及来源解析

农田土壤环境直接关系到农产品质量安全,为了解宁夏引黄灌区农田土壤重金属特征及主要影响因子,在2017～2021年连续5年采样分析了重金属分布特征和相关关系,评价了重金属污染状况和潜在生态风险,解析了农田重金属主要来源.结果表明,宁夏引黄灌区土壤ω(Pb)、ω(As)、ω(Zn)、ω(Ni)、ω(Cu)、ω(Hg)、ω(Cr)和ω(Cd)的平均值分别为19.74、11.67、66.88、29.09、22.55、0.03、62.27和0.19 mg·kg^-1,相对宁夏土壤环境背景值存在一定程度的富集,无中度及以上污染等级.其中,Hg和Cd有中、较高等级生态风险点位,但均未超过农用地土壤污染风险管控值,所有样点均无高风险和极高风险等级.正定矩阵因子分解模型(PMF)与相关分析结合得出的源解析结果表明,研究区农田土壤重金属来源有自然来源、工矿活动及居民生产生活的混合源、交通运输源、农业生产活动源和工业源等5个主要来源,贡献率分别为26.54%、25.59%、22.52%、15.63%和9.72%.综合来看,宁夏引黄灌区农田土壤重金属均无超标现象,且没有较高等级生态风险,农...