CN112595834B - 一种土壤重金属溯源及污染途径确定方法 - Google Patents
一种土壤重金属溯源及污染途径确定方法 Download PDFInfo
- Publication number
- CN112595834B CN112595834B CN202011406026.8A CN202011406026A CN112595834B CN 112595834 B CN112595834 B CN 112595834B CN 202011406026 A CN202011406026 A CN 202011406026A CN 112595834 B CN112595834 B CN 112595834B
- Authority
- CN
- China
- Prior art keywords
- pollution
- area
- farmland
- determining
- source
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000002689 soil Substances 0.000 title claims abstract description 33
- 229910001385 heavy metal Inorganic materials 0.000 title claims abstract description 24
- 238000011160 research Methods 0.000 claims abstract description 31
- 230000002285 radioactive effect Effects 0.000 claims abstract description 12
- 230000035772 mutation Effects 0.000 claims abstract description 10
- 238000012544 monitoring process Methods 0.000 claims description 44
- 239000011159 matrix material Substances 0.000 claims description 25
- 238000011109 contamination Methods 0.000 claims description 16
- 238000001514 detection method Methods 0.000 claims description 14
- 230000037361 pathway Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 3
- 238000000513 principal component analysis Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000002159 abnormal effect Effects 0.000 claims description 2
- 229910052793 cadmium Inorganic materials 0.000 claims description 2
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 238000013459 approach Methods 0.000 description 5
- 230000002265 prevention Effects 0.000 description 2
- 238000003900 soil pollution Methods 0.000 description 2
- 238000009360 aquaculture Methods 0.000 description 1
- 244000144974 aquaculture Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000491 multivariate analysis Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/245—Earth materials for agricultural purposes
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Food Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Processing Of Solid Wastes (AREA)
Abstract
本申请涉及一种农田土壤重金属溯源及污染途径确定方法,包括研究区域确定及数据获取、污染差异度计算、污染途径确定和污染源贡献率,根据污染差异度确定污染差异度突变地块位置,并在所述突变地块进行辐射性地块布设,进一步确定潜在污染源位置及污染途径,大大降低了污染源溯源的工作量,并提高了农田土壤重金属溯源效率。
Description
技术领域
本发明属于农业环境技术领域,具体涉及一种农田土壤重金属溯源方法及污染途径确定方法。
背景技术
农田土壤重金属污染严重且来源复杂,直接影响农产品质量安全与人体健康。目前对于农田污染土壤的研究工作多局限于重金属分布、赋存形态危害等方面,污染源溯源研究较少且多采用同位素示踪、空间分析及多元统计结合分析。此外,关于污染途径确认及潜在污染源确定的相关研究缺乏。
上述技术方法存在的主要问题有:(1)现有污染源溯源方法复杂,实施过程耗时长成本高;(2)现有方法难以根据已有土壤检测数据直接确定潜在污染源;(3)农田土壤污染防治措施缺乏针对性,人力资金等资源过度投入。因此,需要一种能够对土壤重金属污染溯源并能够确定污染途径的方法。
发明内容
本发明利用研究区域农田及已知污染源区域农田土壤数据,明确已知污染源是否为实际污染源,若为实际污染源则明确其污染途径;此外,通过污染途径,明确潜在污染源位置。同时通过计算所有污染源的贡献率,明确主要污染源,为区域性土壤污染防治工作提供决策支撑,制定针对性方案措施。
本发明的目的是提供一种农田土壤重金属溯源及污染途径确定方法。
为了解决上述技术问题,本发明公开了一种农田土壤重金属溯源及污染途径确定方法,所述方法包括:
(1)研究区域确定及数据获取;
(1.1)选取研究区域并明确背景信息,
(1.2)污染区域面积确定,
具体公式如下:
S=max(S研究农田,S污染区域)
其中,S为研究区域内农田面积,S研究农田为研究农田面积,S污染区域是以污染源为中心,向四周辐射1-5km的区域内农田面积,面积单位为亩;
(1.3)点位布设,
(1.4)获取点位数据
获取所述点位的土壤,并对土壤进行检测,获取监测指标对应检测结果;
(1.5)对研究农田的土壤检测结果进行主成分分析,获取载荷矩阵X,所述载荷矩阵X的行为指标,列为主成分,获取得分系数矩阵Y,所述得分系数矩阵Y的行为监测指标,列为主成分,根据所述载荷矩阵X可得到任意单一监测指标Ii:
Ii=a1F1+a2F2+…aiFi
其中,Ii为单一监测指标,仅示意无取值;Fi为单一监测指标对应主成分,仅示意无取值;ai指成分矩阵中单一监测指标Ii对应的数值;
同时根据Y得到任意Fi:
Fi=m1I1+m2I2+…+miIi
其中:mi指得分系数矩阵中单一监测指标Ii对应的数值;
(2)污染差异度计算
(2.1)对于污染源区域内农田,所述区域内农田点位的监测指标均存在以下关系式,并可通过最小二乘法求解wi:
Ei(Fi)=∑Cimi(Fi)
Ci=w1(Fi)E1+w2(Fi)E2+…+wi(Fi)Ei
其中,Ei(Fi)为污染源地块单一监测指标Ii(如Cd)对应Fi下的成分因子,Ci为监测指标对应检测结果,mi(Fi)为得分系数矩阵中单一监测指标Ii对应Fi下的数值,wi(Fi)为单一监测指标Ii对应Fi下的值;
(2.2)依据单一监测指标Ii的ai(Fi)和wi(Fi)分别构建矩阵,求单一监测指标污染差异度,确定已知污染源是否为实际污染源,具体公式如下:
其中cosθi为单一监测指标污染差异度,a为ai(Fi)组成的矩阵,且ai(Fi)为成分矩阵中单一监测指标Ii对应Fi下的数值,w为wi(Fi)组成的矩阵;所述差异度≤0.85认为已知污染源为研究区域实际污染源;>0.85则认为已知污染源不是研究区域实际污染源;
(3)污染途径确定;
污染途径确定,其特征在于:
(3.1)选取污染差异度≤0.85的污染源,依据区域背景信息预设污染途径,并布设污染途径地块,所述污染途径地块内点位数量、监测指标信息与污染区域保持一致,
(3.2)潜在污染源确定
依据步骤(1)-(2)的步骤计算步骤(3.1)确定的所述污染途径地块的污染差异度,确定差异度突变地块,针对所述突变地块进行辐射性地块布设,获取辐射性地块的污染差异度,循环此操作,直至差异度稳定变化,
所述辐射性地块内点位数量、监测指标等信息与污染区域保持一致,每次循环布设一个或多个辐射性地块;
(4)污染源贡献率
依据上述已知污染源及潜在污染源,获取污染地块所有监测指标的w值,所有污染地块与研究地块存在关系如下式;联立所有监测指标关系式,通过最小二乘法求解ki,ki即为各污染源地块的贡献率:
ai(Fi)=k1wi (1) (Fi)+k2wi (2) (Fi)+…+kiwi (n) (Fi)
其中,wi (1) (Fi)、wi (2) (Fi)、wi (n) (Fi)与wi求法一致,仅用于区分不同污染源地块;
所述ki和为1。
进一步地,所述辐射性地块差异度稳定出现变化——越接近突变地块差异度越小,则最终布设的辐射性地块为潜在污染源所在位置;
进一步地,所述背景信息包括农田位置、类型、面积,已知污染源数量、类型、位置、面积、时期,研究区域水系流向、风向、道路,研究农田到各污染源间的距离;
进一步地,采取均匀布点法计算各区域内的点位数量,计算公式为:
其中,N为点位数量,计算结果为小数时向上取整;n为监测指标数量,且n≥5;
所述研究农田及所述污染区域内农田中点位布设数量一致;
进一步地,所述监测的检测指标包括但不限于Cd、As、Pb、Cr、Cu、Zn、Ni、Mg、Fe、Mn、Ca、Si、Se;
进一步地,所述预设污染途径依据区域背景信息,取研究农田与各污染源连线,并在连线上设置污染途径地块,所述污染途径地块中心点设置在所述连线上;
进一步地,每条所述污染源连线上设置的所述污染途径地块布设数量≥1;
进一步地,所述确定差异度突变地块为所述研究农田、污染途径地块及污染源的污染差异度呈单调递增或递减变化过程中的异常污染差异度对应的地块;
本发明的农田土壤重金属溯源及污染途径确定方法,具有以下优点:
1.点位数据获取时间快,土壤检测结果实时获取;点位数据获取成本低:传统判别是否为污染源,需进行同位素分析,成本较高,本发明单个点位数据获取成本仅为常规成本的10-20%;
2.依据已有土壤数据,可推断潜在污染源位置:现有的方法只能分析已知污染源,已经关停的企业或停用的灌渠等造成的重金属等污染,现有方法无法对其分析;
3.已知污染源,可同时推测污染源贡献率及污染途径,更利于制定针对性决策措施建议,减少现有大区域措施的人力物资等资源浪费;
4.根据污染差异度确定污染差异度突变地块位置,并在所述突变地块进行辐射性地块布设进一步确定潜在污染源,大大降低了工作量,为土壤重金属污染的工作提供了技术支持。
附图说明
图1为农田土壤重金属溯源及污染途径确定方法的技术流程图;
图2为某研究区域示意图;
图3为研究区域点位分布图;
图4为各地块污染差异度示意图;
图5为污染途径示意图;
图6为污染区域W值示意图;
具体实施方式
下面通过实施例对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
应当理解,本文所使用的诸如“具有”,“包含”以及“包括”术语并不排除一个或多个其它元件或其组合的存在或添加。
1.研究区域确定及数据获取
①选取某地内农田M为研究对象,确定该乡镇存在两个污染源A(工矿业)和B(养殖业);M到A距离约为8.5km,M到B距离约为3km(图2);
②确定农田M的面积为780亩,以污染源为中心划分相同面积农田区域;
③点位布设:依据点位布设公式,各区域地块均匀布设41个点位,获取点位地理信息(图3);
④点位数据获取:A、B、M三个区域所有点位通过X荧光快速检测仪器检测Cd、As、Pb、Cr、Cu、Zn、Ni、Mg、Fe、Mn、Ca、Si、Se等指标,现场获取检测结果。
2.污染差异度计算
①以Cd为研究指标,对农田M点位信息进行主成分分析,得到载荷矩阵和得分系数矩阵;
②依据差异度计算公式,计算污染区域污染差异度,其中A、B区域差异度均小于0.85,见图4。
③依据区域水系流向、风向等信息,在农田M到污染源A、B间均预设污染途径,并根据污染途径距离布设污染途径地块,其中M到A布设3个地块(C、D、E),M到B布设1个地块(F);
3.污染途径确定
①依据污染差异度,可知区域D污染差异度突增(由0.014突增至0.029),说明其与农田M差异度增强,可能受其他污染源影响,在其两边预设相反污染途径并布设地块(G、H),布设点位并获取点位地理信息;
②通过X荧光快速检测仪器获取G、H点位检测数据,计算区域污染差异度,确定三条污染途径(MA、MB、DH),见图5。
4.污染源贡献率确定
①依据最小二乘法,计算出污染区域A、B所有指标wi值,见图6;
②通过农田的成分矩阵中的ai值及污染区域wi值,计算得到A的贡献率为0.865,B的贡献率为0.135,确定A与B均为农田M的污染源,且A为农田M的主要污染源,判定结果与实际认知相符。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域,对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的实施例。
Claims (10)
1.一种土壤重金属溯源及污染途径确定方法,其特征在于,所述方法包括:
(1)研究区域确定及数据获取;
(1.1)选取研究区域并明确背景信息,
(1.2)污染区域面积确定,
(1.3)点位布设,
(1.4)获取点位监测数据
获取所述点位的土壤,并对土壤进行检测,获取监测指标对应检测结果;
(1.5)对研究农田的土壤检测结果进行主成分分析,获取载荷矩阵X,所述载荷矩阵X的行为监测指标,列为主成分,获取得分系数矩阵Y,所述得分系数矩阵Y的行为监测指标,列为主成分,根据所述载荷矩阵X可得到单一监测指标Ii:
Ii=a1F1+a2F2+…aiFi
其中,Ii为单一监测指标,仅示意无取值;Fi为单一监测指标对应主成分,仅示意无取值;ai指成分矩阵中单一监测指标Ii对应的数值;
同时根据Y得到任意Fi:
Fi=m1I1+m2I2+…+miIi
其中:mi指得分系数矩阵中单一监测指标Ii对应的数值;
(2)污染差异度计算
(2.1)对于污染源区域内农田,所述农田内点位的监测指标均
存在以下关系式,并可通过最小二乘法求解wi:
Ei(Fi)=∑Cimi(Fi)
Ci=w1(Fi)E1+w2(Fi)E2+…+wi(Fi)Ei
其中,Ei(Fi)为污染源地块单一监测指标Ii对应Fi下的成分因子,Ci为监测指标Ii对应检测结果,mi(Fi)为得分系数矩阵中单一监测指标Ii对应Fi下的数值,wi(Fi)为单一监测指标Ii对应Fi下的值;
(2.2)依据单一监测指标Ii的ai(Fi)和wi(Fi)分别构建矩阵,求单一监测指标污染差异度,确定已知污染源是否为实际污染源,具体公式如下:
其中cosθi为单一监测指标污染差异度,a为ai(Fi)组成的矩阵,且ai(Fi)为成分矩阵中单一监测指标Ii对应Fi下的数值,w为wi(Fi)组成的矩阵;所述差异度≤0.85认为已知污染源为研究区域实际污染源;>0.85则认为已知污染源不是研究区域实际污染源;
(3)污染途径确定;
污染途径确定,其特征在于:
(3.1)选取污染差异度≤0.85的污染源,依据区域背景信息预设污染途径,并布设污染途径地块,所述污染途径地块内点位数量、监测指标信息与污染区域保持一致,
(3.2)潜在污染源确定
依据步骤(1)-(2)的步骤计算步骤(3.1)确定的所述污染途径地块的污染差异度,确定差异度突变地块,针对所述突变地块进行辐射性地块布设,获取辐射性地块的污染差异度,循环此操作,直至差异度稳定变化,
所述辐射性地块点位数量、监测指标信息与污染区域保持一致,每次循环布设一个或多个辐射性地块;
(4)污染源贡献率
依据上述已知污染源及潜在污染源,获取污染地块所有监测指标的w值,所有污染地块与研究地块存在关系如下式;联立所有监测指标关系式,通过最小二乘法求解ki,ki即为各污染源地块的贡献率:
ai(Fi)=k1wi (1) (Fi)+k2wi (2) (Fi)+…+kiwi (n) (Fi)
其中,wi (1) (Fi)、wi (2) (Fi)、wi (n) (Fi)与wi求法一致,仅用于区分不同污染源地块;
所述ki和为1。
2.如权利要求1所述的一种土壤重金属溯源及污染途径确定方法,其特征在于,所述背景信息包括农田位置、类型、面积,已知污染源数量、类型、位置、面积、时期,研究区域水系流向、风向、道路,研究农田到各污染源间的距离。
4.如权利要求1所述的一种土壤重金属溯源及污染途径确定方法,其特征在于,所述监测的检测指标包括但不限于Cd、As、Pb、Cr、Cu、Zn、Ni、Mg、Fe、Mn、Ca、Si、Se。
5.如权利要求1所述的一种土壤重金属溯源及污染途径确定方法,其特征在于,所述预设污染途径依据区域背景信息,取研究农田与各污染源连线,并在连线上设置污染途径地块,所述污染途径地块中心点设置在所述连线上。
6.如权利要求5所述的一种土壤重金属溯源及污染途径确定方法,其特征在于,每条所述污染源连线上设置的所述污染途径地块布设数量≥1。
7.如权利要求6所述的一种土壤重金属溯源及污染途径确定方法,其特征在于,所述污染源连线上设置的所述污染途径地块布设数量依据所述连线的距离确定。
8.如权利要求6所述的一种土壤重金属溯源及污染途径确定方法,其特征在于,所述污染途径地块在所述污染源连线上均匀设置。
9.如权利要求1所述的一种土壤重金属溯源及污染途径确定方法,其特征在于,所述确定差异度突变地块为所述研究农田、污染途径地块及污染源地块的污染差异度呈单调递增或递减变化过程中的异常污染差异度对应的地块。
10.如权利要求1所述的一种土壤重金属溯源及污染途径确定方法,其特征在于,所述污染区域面积的确定公式如下:
S=max(S研究农田,S污染区域)
其中,S为研究区域内农田面积,S研究农田为研究农田面积,S污染区域是以污染源为中心,向四周辐射1-5km的区域内农田面积,面积单位为亩。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011406026.8A CN112595834B (zh) | 2020-12-03 | 2020-12-03 | 一种土壤重金属溯源及污染途径确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011406026.8A CN112595834B (zh) | 2020-12-03 | 2020-12-03 | 一种土壤重金属溯源及污染途径确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112595834A CN112595834A (zh) | 2021-04-02 |
CN112595834B true CN112595834B (zh) | 2022-11-04 |
Family
ID=75188285
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011406026.8A Active CN112595834B (zh) | 2020-12-03 | 2020-12-03 | 一种土壤重金属溯源及污染途径确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112595834B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114210718B (zh) * | 2021-12-09 | 2023-03-24 | 农业农村部环境保护科研监测所 | 一种农产品污染治理修复方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107767032A (zh) * | 2017-09-27 | 2018-03-06 | 北京农业信息技术研究中心 | 一种农田土壤重金属污染决策***及方法 |
CN111474307A (zh) * | 2020-03-26 | 2020-07-31 | 南方科技大学 | 污染物溯源方法、装置、计算机设备及存储介质 |
AU2020101615A4 (en) * | 2020-08-02 | 2020-09-10 | Beijing Normal University | A Method for Source Apportionment of PAHs in Roadway Sediments Coupled with Transport and Transformation Process |
CN111678969A (zh) * | 2020-06-05 | 2020-09-18 | 农业农村部环境保护科研监测所 | 利用土壤剖面表层重金属累积比例解析重金属污染来源的方法 |
CN111751513A (zh) * | 2020-06-04 | 2020-10-09 | 农业农村部环境保护科研监测所 | 一种稻田剖面土壤重金属累积及形态特征解析重金属主要污染源的分析方法 |
CN111766368A (zh) * | 2020-06-30 | 2020-10-13 | 中国环境科学研究院 | 一种铅同位素的重金属源解析方法 |
CN111797918A (zh) * | 2020-06-29 | 2020-10-20 | 平安国际智慧城市科技股份有限公司 | 大气污染源识别方法、装置、计算机设备及存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101349306B1 (ko) * | 2013-09-13 | 2014-01-09 | 한국지질자원연구원 | 연속 추출 방법과 동위 원소 분석에 의한 중금속 오염원 분석 방법 |
KR101511292B1 (ko) * | 2014-12-30 | 2015-04-13 | 한국지질자원연구원 | 연속 추출 방법과 안정 동위 원소 분석을 통한 다수 원인자에 의한 토양 오염 기여도 배분 방법 |
-
2020
- 2020-12-03 CN CN202011406026.8A patent/CN112595834B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107767032A (zh) * | 2017-09-27 | 2018-03-06 | 北京农业信息技术研究中心 | 一种农田土壤重金属污染决策***及方法 |
CN111474307A (zh) * | 2020-03-26 | 2020-07-31 | 南方科技大学 | 污染物溯源方法、装置、计算机设备及存储介质 |
CN111751513A (zh) * | 2020-06-04 | 2020-10-09 | 农业农村部环境保护科研监测所 | 一种稻田剖面土壤重金属累积及形态特征解析重金属主要污染源的分析方法 |
CN111678969A (zh) * | 2020-06-05 | 2020-09-18 | 农业农村部环境保护科研监测所 | 利用土壤剖面表层重金属累积比例解析重金属污染来源的方法 |
CN111797918A (zh) * | 2020-06-29 | 2020-10-20 | 平安国际智慧城市科技股份有限公司 | 大气污染源识别方法、装置、计算机设备及存储介质 |
CN111766368A (zh) * | 2020-06-30 | 2020-10-13 | 中国环境科学研究院 | 一种铅同位素的重金属源解析方法 |
AU2020101615A4 (en) * | 2020-08-02 | 2020-09-10 | Beijing Normal University | A Method for Source Apportionment of PAHs in Roadway Sediments Coupled with Transport and Transformation Process |
Non-Patent Citations (4)
Title |
---|
Heavy metals and organic compounds contamination in Heavy metals and organic compounds contamination in soil from an e-waste region in South China;Ming Liu等;《Environmental Science Processes & Impacts》;20131231;第15卷;第919–929页 * |
Tracing source and distribution of heavy metals in road dust, soil and soakaway sediment through speciation and isotopic fingerprinting;Manish Kumar等;《Geoderma》;20130730;第211-212卷;第8–17页 * |
丹江口水库新增淹没区农田土壤重金属源解析;韩培培等;《中国环境科学》;20161231;第36卷(第8期);第2437-2443页 * |
湘潭县农田土壤重金属污染及生态风险评价;刘瑞雪等;《农业环境科学学报》;20190731;第38卷(第7期);第1523-1530页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112595834A (zh) | 2021-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wen et al. | Heavy metals in agricultural soil in China: A systematic review and meta-analysis | |
Lu et al. | Multivariate and geostatistical analyses of the spatial distribution and origin of heavy metals in the agricultural soils in Shunyi, Beijing, China | |
Zhang et al. | Multivariate and geostatistical analyses of the sources and spatial distribution of heavy metals in agricultural soil in Gongzhuling, Northeast China | |
Li et al. | Regional simulation of nitrate leaching potential from winter wheat-summer maize rotation croplands on the North China Plain using the NLEAP-GIS model | |
CN103106347A (zh) | 一种基于土壤属性空间分布的农业面源磷污染估算方法 | |
CN112595834B (zh) | 一种土壤重金属溯源及污染途径确定方法 | |
CN113020232B (zh) | 一种受污染耕地综合治理及动态调控方法 | |
CN114210718B (zh) | 一种农产品污染治理修复方法 | |
Qiao et al. | Identifying the relationships and drivers of agro-ecosystem services using a constraint line approach in the agro-pastoral transitional zone of China | |
Ghiglieri et al. | Potential risks of nitrate pollution in aquifers from agricultural practices in the Nurra region, northwestern Sardinia, Italy | |
Djodjic et al. | Optimizing placement of constructed wetlands at landscape scale in order to reduce phosphorus losses | |
CN114219262A (zh) | 一种区域农产品超标概率预测及误差分析方法 | |
Xiao et al. | Water quality assessment and pollution source apportionment using multivariate statistical techniques: a case study of the Laixi River Basin, China | |
Bi et al. | The impact of land use and socio-economic factors on ammonia nitrogen pollution in Weihe River watershed, China | |
Song et al. | High-throughput sequencing clarifies the spatial structures of microbial communities in cadmium-polluted rice soils | |
Qu et al. | Spatial-temporal characteristics and driving factors of the chemical fertilizer supply/demand correlation network in China | |
Agyeman et al. | Combination of enrichment factor and positive matrix factorization in the estimation of potentially toxic element source distribution in agricultural soil | |
Anthony et al. | Soil nitrogen and phosphorus behavior in a long‐term fertilization experiment | |
CN117648793A (zh) | 基于rusle模型的流域土壤侵蚀分析*** | |
CN105975647A (zh) | 一种高耐磨铸造磨球成分设计与热处理制度的可视化方法 | |
Tian et al. | The effect of the urban wastewater treatment ratio on agricultural water productivity: based on provincial data of China in 2004–2010 | |
Feng et al. | The impacts of spatio-temporal variation of natural and agricultural influences on the environmental water quality in a fluvial-lacustrine watershed in China | |
Zou et al. | Analysis of soil fertility and optimal nitrogen application of brown earth (luvisols) in China | |
Yan et al. | A knowledge model system for wheat production management | |
Li et al. | Suitability evaluation of remediation technology for polluted farmland |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |