• ISSN 1008-505X
  • CN 11-3996/S

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于GIS和丰缺指标法的区域施肥管理体系的构建

高宏艳 索全义 郑海春 郜翻身

引用本文:
Citation:

基于GIS和丰缺指标法的区域施肥管理体系的构建

    作者简介: 高宏艳 E-mail:905868022@qq.com;
    通讯作者: 索全义, E-mail:paul98@sina.com
  • 基金项目: 农业农村部测土配方施肥项目

Construction of regional fertilization system based on GIS and nutrient abundance index

    Corresponding author: SUO Quan-yi, E-mail:paul98@sina.com ;
  • 摘要:   【目的】  采用GIS定位和养分丰缺指标相结合的方法,研究西辽河平原耕地土壤碱解氮、有效磷、速效钾的空间分布特点,建立直观、精准的推荐施肥系统。  【方法】  西辽河平原包括科尔沁区、开鲁县、奈曼旗、科尔沁左翼后旗、科尔沁左翼中旗。以该区域的15421个土壤测试数据和143个玉米“3414试验”数据为基础,运用GIS方法计算耕地土壤碱解氮、有效磷、速效钾的空间插值结果,结合土壤丰缺指标,将西辽河平原耕地土壤碱解氮、有效磷、速效钾养分分区,确定不同土壤氮磷钾养分组合下的施肥量,并制作玉米氮磷钾分区施肥图。  【结果】  西辽河平原耕地土壤养分空间分布不均,大部分土壤碱解氮含量较低,有效磷、速效钾含量中等。耕地土壤有效氮、磷、钾养分含量呈低氮–中磷–中钾 (碱解氮 ≤ 69.5 mg/kg、有效磷4.18~10.40 mg/kg、速效钾90.3~152.5 mg/kg)、中氮–中磷–中钾 (碱解氮69.5~107.9 mg/kg、有效磷4.18~10.40 mg/kg、速效钾90.3~152.5 mg/kg)、中氮–中磷–高钾 (碱解氮69.5~107.9 mg/kg、有效磷4.18~10.40 mg/kg、速效钾 > 152.5 mg/kg) 的面积占比分别为42.4%、14.1%、13.5%。土壤碱解氮、有效磷、速效钾含量 (x) 与经济最佳施肥量 (y) 之间存在较好的对数关系,氮施肥模型y = –102.5ln(x) + 617.22,磷施肥模型y = –36.11ln(x) + 174.1,钾施肥模型为y = –25.89ln(x) + 180.96。不同土壤养分组合种植玉米所需的经济最佳N–P2O5–K2O用量 (kg/hm2) 组合分别为202–110–58、168–107–57、164–102–47。  【结论】  土壤碱解氮、有效磷、速效钾含量与相对产量存在较好的对数函数关系。西辽河灌区土壤有效氮磷钾丰缺指标的划分以75%和85%作为相对产量的分级区间较为合适。在土壤有效氮磷钾养分大数据的支撑下,采用GIS和丰缺指标法相结合的方法,可以准确快捷地构建区域施肥管理体系。
  • 图 1  土壤采样点分布图

    Figure 1.  Distribution map of soil sampling points

    图 2  试验点分布图

    Figure 2.  Distribution map of test points

    图 3  玉米最佳氮磷钾肥经济施用量与土壤养分含量的函数关系图

    Figure 3.  The function relationship between optimal nutrient rate and soil available nutrient content

    图 4  西辽河平原土壤养分分区及氮磷钾肥料最佳用量

    Figure 4.  Soil nutrient zoning and the optimal fertilizer rate in Xiliao River Plain

    表 1  不同产量目标施肥量等级的试验数量占比

    Table 1.  Proportion of test number at different levels of fertilizer

    产量目标
    Yield target
    养分
    Nutrient
    总试验数量
    Total test
    number
    平均施肥量
    Mean rate
    (kg/hm2)
    各级施肥量试验数量占比 (%)
    The proportion of test number under each nutrient rate
    ≤ 5050~100100~150150~200 ≥ 200
    最佳经济产量
    Optimal economical yield
    N137178.002.232.831.433.6
    P2O13098.98.543.038.57.72.3
    K2O12464.022.670.26.40.80
    最高产量
    Maximun yield
    N137209.30019.730.749.6
    P2O130127.70.830.049.211.58.5
    K2O12474.48.979.09.72.40
    下载: 导出CSV

    表 2  不同相对产量划分的土壤养分丰缺指标

    Table 2.  Soil nutrient abundance index of different relative yields

    土壤养分
    Soil nutrient
    相对产量类型
    Relative
    yield type
    丰缺指标 Abundance index (mg/kg)相对产量与土壤
    养分关系式
    Relative yield and soil
    nutrient function
    R2FP
    低 Low
    ( ≤ 75%)
    中 Medium
    (75%~85%)
    高 High
    (>85%)
    碱解氮
    Alkeline-N
    RY1 ≤ 68.268.2~107.2>107.2y = 22.123ln (x) – 18.4180.49095.057<0.001
    RY2 ≤ 69.269.2~109.3>109.3y = 21.902ln (x) – 17.8010.49295.631<0.001
    RY3 ≤ 69.569.5~107.9>107.9y = 22.716ln (x) – 21.3320.553122.237<0.001
    有效磷
    Available P
    RY1 ≤ 4.094.09~10.14>10.14y = 11.022ln (x) + 59.4670.610140.775<0.001
    RY2 ≤ 4.184.18~10.40>10.40y = 10.985ln (x) + 59.2770.630151.759<0.001
    RY3 ≤ 3.803.80~10.32>10.32y = 10.019ln (x) + 61.6140.611141.495<0.001
    速效钾
    Available K
    RY1 ≤ 90.390.3~152.5>152.5y = 19.062ln (x)–10.8290.606143.137<0.001
    RY2 ≤ 90.990.9~153.9>153.9y = 18.964ln (x)–10.5130.604141.805<0.001
    RY3 ≤ 78.378.3~141.7>141.7y = 16.858ln (x) + 1.48570.542111.111<0.001
    注(Note):RY1 = 无肥区产量/最经济产量; RY2 = 无肥区产量/最高产量; RY3 = 无肥区产量/N2P2K2 处理产量。 函数关系方程中,y 代表相对产量,x 代表土壤养分含量。RY1 = no-fertilizer yield/optimum yield; RY2 = no-fertilizer yield/maximum yield; RY3 = no-fertilizer yield/N2P2K2 yield. In the function, y and x represent relative yield and nutrient content, respectively.
    下载: 导出CSV

    表 3  不同丰缺指标下玉米最佳N、P2O5、K2O肥料用量

    Table 3.  Maize optimal fertilizer use under different ranges of soil nutrient content

    土壤养分
    Soil nutrient
    丰缺程度
    Grade
    相对产量
    Relative yield
    养分含量 (mg/kg)
    Nutrient content
    经济施肥范围 (kg/hm2)
    Optimal fertilizerr ate
    碱解氮 Alkeline-N低 Low ≤ 75% ≤ 69.5182.5~235.6
    中 Medium75%~85%69.5~107.9137.4~182.5
    高 High > 85%>107.944.4~137.4
    有效磷 Available P低 Low ≤ 75% ≤ 4.18122.4~151.1
    中 Medium75%~85%4.18~10.4089.5~122.4
    高 High > 85%>10.4067.6~89.5
    速效钾 Available K低 Low ≤ 75% ≤ 90.364.4~77.2
    中 Medium75%~85%90.3~152.550.8~64.4
    高 High > 85%>152.539.6~50.8
    下载: 导出CSV
  • [1] 唐春闺, 李强, 李帆, 等. 基于GIS和FCM的宁乡植烟土壤养分管理分区[J]. 湖北农业科学, 2017, 56(6): 1029–1032, 1039. Tang C G, Li Q, Li F, et al. Definition of management zones of tobacco soil nutrients based on GIS and FCM in Ningxiang[J]. Hubei Agricultural Sciences, 2017, 56(6): 1029–1032, 1039.
    [2] 王宜伦, 苏瑞光, 刘举, 等. 养分专家系统推荐施肥对潮土夏玉米产量及肥料效率的影响[J]. 作物学报, 2014, 40(3): 563–569. Wang Y L, Su R G, Liu J, et al. Effects of nutrient expert recommend fertilization on yield and fertilizer efficiency of summer maize in fluvo-aquic soil[J]. Acta Agronomica Sinica, 2014, 40(3): 563–569. doi:  10.3724/SP.J.1006.2014.00563
    [3] Bongiovanni R G, Robledo C W, Lambent D M. Economics of site-specific nitrogen management for protein content in wheat[J]. Computers and Electronics in Agriculture, 2007, 58: 13–24. doi:  10.1016/j.compag.2007.01.018
    [4] Fridgen J J, Kitchen N R, Sudduth K A, et al. Management zone analyst (MZA): Software for subfield management zone delineation[J]. Agronomy Journal, 2004, 96(1): 100–108. doi:  10.2134/agronj2004.0100
    [5] 武德传, 罗红香, 宋泽民, 等. 黔南山地植烟土壤主要养分空间变异和管理分区[J]. 应用生态学报, 2014, 25(6): 1701–1707. Wu D C, Luo H X, Song Z M, et al. Spatial variability and management zone of soil major nutrients in tobacco fields in Qiannan mountainous region[J]. Chinese Journal of Applied Ecology, 2014, 25(6): 1701–1707.
    [6] 郭澎涛, 李茂芬, 林钊沐, 等. 基于多源环境变量的橡胶园土壤管理分区[J]. 农业工程学报, 2014, 30(12): 96–104. Guo P T, Li M F, Lin Z M, et al. Delineating soil management zones in rubber plantation using multisource data of environmental variables[J]. Transactions of the Chinese Society of Agricultural Engineering, 2014, 30(12): 96–104. doi:  10.3969/j.issn.1002-6819.2014.12.012
    [7] 车升国. 区域作物专用复合(混)肥料配方制定方法与应用[D]. 北京: 中国农业大学博士学位论文, 2015.

    Che S G. Design method and application of formula of regional crop-based compound fertilize[D]. Beijing: PhD Dissertaion of China Agricultural University, 2015.
    [8] 张树清. 甘肃农田氮磷钾养分资源特征研究[D]. 兰州: 甘肃农业大学博士学位论文, 2001.

    Zhang S Q. Characteristics of nitrogen, phosphorus and potassium nutrient resources in farmlands of Gansu[D]. Lanzhou: PhD Dissertaion of Gansu Agricultural University, 2001.
    [9] Mzuku M, Khosla R, Reich R, et al. Spatial variability of measured soil properties across site-specific management zones[J]. Soil Science Society of America Journal, 2005, 69(5): 1572–1579. doi:  10.2136/sssaj2005.0062
    [10] 王子龙, 付强, 姜秋香, 等. 基于粒子群优化算法的土壤养分管理分区[J]. 农业工程学报, 2008, 24(10): 80–84. Wang Z L, Fu Q, Jiang Q X, et al. Soil nutrient management zones based on particle swarm optimization algorithm[J]. Transactions of the Chinese Society of Agricultural Engineering, 2008, 24(10): 80–84. doi:  10.3321/j.issn:1002-6819.2008.10.016
    [11] 串丽敏, 何萍, 赵同科. 作物推荐施肥方法研究进展[J]. 中国农业科技导报, 2016, 18(1): 95–102. Chuan L M, He P, Zhao T K. Research advance on recommendation for crop fertilization methodology[J]. Journal of Agricultural Science and Technology, 2016, 18(1): 95–102.
    [12] 杨俐苹, 白由路, 王贺, 等. 测土配方施肥指标体系建立中“3414”试验方案应用探讨——以内蒙古海拉尔地区油菜“3414”试验为例[J]. 植物营养与肥料学报, 2011, 17(4): 1018–1023. Yang L P, Bai Y L, Wang H, et al. Application of “3414”field trial design for establishing soil testing and fertilizer recommendation index[J]. Journal of Plant Nutrition and Fertilizers, 2011, 17(4): 1018–1023.
    [13] 章明清, 李娟, 孔庆波, 等. 作物肥料效应函数模型研究进展与展望[J]. 土壤学报, 2016, 53(6): 1343–1356. Zhang M Q, Li J, Kong Q B, et al. Progress and prospect of the study on crop-response-to-fertilization function model[J]. Acta Pedologica Sinica, 2016, 53(6): 1343–1356.
    [14] 周伟, 周亚星, 马日亮, 等. 西辽河平原灌区玉米种植密度调查分析[J]. 天津农林科技, 2015, (S1): 29–32. Zhou W, Zhou Y X, Ma R L, et al. Investigation and analysis of maize planting density in Xiliaohe plain irrigation area[J]. Science and Technology of Tianjin Agriculture and Forest, 2015, (S1): 29–32.
    [15] 童倩倩, 李莉婕, 赵泽英, 等. 基于GIS的贵州稻田土壤养分管理分区[J]. 西南农业学报, 2017, 30(12): 2727–2731. Tong Q Q, Li L J, Zhao Z Y, et al. Management subarea of paddy soil nutrients based on GIS in Guizhou[J]. Southwest China Journal of Agricultural Sciences, 2017, 30(12): 2727–2731.
    [16] 张南冰, 潘艳花. 河西内陆灌区盐渍土甜菜氮磷钾肥最大利润施用量研究[J]. 中国糖料, 2017, 39(4): 14–17. Zhang N B, Pan Y H. NPK fertilizer application amount for maximum profit rate of sugarbeet in irrigation area of Hexi inland saline soil[J]. Sugar Crops of China, 2017, 39(4): 14–17.
    [17] Yao X, Yu K Y, Deng Y B, et al. Spatial variability of soil organic carbon and total nitrogen in the hilly red soil region of Southern China[J]. Journal of Forestry Research, 2020, 31(6): 2385–2394. doi:  10.1007/s11676-019-01014-8
    [18] 潘峰, 陶杰, 冯建军, 等. 内蒙古扎鲁特旗春玉米推荐施肥指标体系研究[J]. 内蒙古民族大学学报(自然科学版), 2011, 26(4): 424–428. Pan F, Tao J, Feng J J, et al. Study on index system of optimal fertilizer treatment on spring corn in Zhalute Banner[J]. Journal of Inner Mongolia University for Nationalities, 2011, 26(4): 424–428.
    [19] 向德明, 李强, 黎娟, 等. 湘西植烟土壤养分管理分区研究[J]. 安徽农学通报, 2019, 25(13): 119–122. Xiang D M, Li Q, Li J, et al. Study on soil nutrient management zoning of tobacco planting in Xiangxi[J]. Anhui Agricultural Science Bulletin, 2019, 25(13): 119–122. doi:  10.3969/j.issn.1007-7731.2019.13.046
    [20] 张云贵, 刘青丽, 李志宏, 等. 分散经营条件下烤烟精准养分管理研究[J]. 植物营养与肥料学报, 2015, 21(4): 1066–1074. Zhang Y G, Liu Q L, Li Z H, et al. Research on the field-specific nutrient management of flue-cured tobacco in condition of decentralized operation[J]. Journal of Plant Nutrition and Fertilizers, 2015, 21(4): 1066–1074. doi:  10.11674/zwyf.2015.0427
    [21] 李文彪, 刘荣乐, 郑海春, 等. 内蒙古河套灌区春玉米推荐施肥指标体系研究[J]. 中国农业科学, 2012, 45(1): 93–101. Li W B, Liu R L, Zheng H C, et al. Study on index system of optimal fertilizer recommendation for spring corn in Hetao irrigation area of Inner Mongolia[J]. Scientia Agricultura Sinica, 2012, 45(1): 93–101. doi:  10.3864/j.issn.0578-1752.2012.01.011
    [22] 李文彪, 郑海春, 郜翻身, 等. 内蒙古河套灌区春小麦推荐施肥指标体系研究[J]. 植物营养与肥料学报, 2011, 17(6): 1327–1334. Li W B, Zheng H C, Gao F S, et al. Study on index of fertilizer recommendation for spring wheat in Hetao irrigated area[J]. Journal of Plant Nutrition and Fertilizers, 2011, 17(6): 1327–1334. doi:  10.11674/zwyf.2011.1089
    [23] 刘疆, 翟丙年, 殷振江. 咸阳灌区冬小麦推荐施肥指标体系研究[J]. 西北农林科技大学学报(自然科学版), 2014, 42(6): 55–62. Liu J, Zhai B N, Yin Z J. Index system of fertilizer recommendation for winter wheat in Xianyang irrigation area[J]. Journal of Northwest A& F University (Natural Science Edition), 2014, 42(6): 55–62.
    [24] 付莹莹, 同延安, 赵佐平, 等. 陕西关中灌区夏玉米土壤养分丰缺及推荐施肥指标体系的建立[J]. 干旱地区农业研究, 2010, 28(1): 88–93. Fu Y Y, Tong Y A, Zhao Z P, et al. The establishment of soil nutrient abundance and fertilizer recommendation index system for summer maize in Guanzhong irrigation areas, Shaanxi Province[J]. Agricultural Research in the Arid Areas, 2010, 28(1): 88–93.
    [25] 刘芬, 同延安, 王小英, 等. 陕西关中灌区冬小麦施肥指标研究[J]. 土壤学报, 2013, 50(3): 556–563. Liu F, Tong Y A, Wang X Y, et al. Indices for fertilization of winter wheat in Guanzhong irrigation area of Shaanxi[J]. Acta Pedologica Sinica, 2013, 50(3): 556–563.
    [26] Fageria N K, Baligar V C, Jones C A. Growth and mineral nutrition of field]crops[M]. Boca Raton, USA: CRC Press/Taylor & Francis Group, 2011.
    [27] Troeh F R, Thompson L M. Soil and soil fertility[M]. Lowa, America: Blackwell Publishing, 2005. 251–268.
    [28] 张福锁. 测土配方施肥技术要览[M]. 北京: 中国农业大学出版社, 2006. 93–110.

    Zhang F S. Technical overview of soil test and formula fertilization[M]. Beijing: China Agricultural University Press, 2006. 93–110.
    [29] 陈新平, 张福锁. 通过“3414”试验建立测土推荐施肥技术指标体系[J]. 中国农技推广, 2006, 22(4): 36–39. Chen X P, Zhang F S. The technical index system of recommended fertilization for soil measurement was established through "3414" test[J]. China Agricultural Technology Extension, 2006, 22(4): 36–39. doi:  10.3969/j.issn.1002-381X.2006.04.026
    [30] 黄亿, 李廷轩, 张锡洲, 等. 基于“3414”试验的川中丘陵区油菜施肥指标体系构建[J]. 中国农业科学, 2013, 46(10): 2058–2066. Huang Y, Li T X, Zhang X Z, et al. Establishment of fertilization recommendation indexes of rapeseed soil based on the “3414” field experiments in the middle of Sichuan hilly regions[J]. Scientia Agricultura Sinica, 2013, 46(10): 2058–2066. doi:  10.3864/j.issn.0578-1752.2013.10.011
    [31] 戢林, 张锡洲, 李廷轩. 基于“3414”试验的川中丘陵区水稻测土配方施肥指标体系构建[J]. 中国农业科学, 2011, 44(1): 84–92. Ji L, Zhang X Z, Li T X. Establishing fertilization recommendation index of paddy soil based on the “3414” field experiments in the middle of Sichuan hilly regions[J]. Scientia Agricultura Sinica, 2011, 44(1): 84–92. doi:  10.3864/j.issn.0578-1752.2011.01.010
    [32] 宋朝玉, 高峻岭, 张清霞, 等. “3414”肥料试验结果统计分析方法的探讨[J]. 山东农业科学, 2009, (9): 93–96. Song C Y, Gao J L, Zhang Q X, et al. Discussion on statistical analysis method of "3414" fertilizer test results[J]. Shandong Agricultural Sciences, 2009, (9): 93–96. doi:  10.3969/j.issn.1001-4942.2009.09.031
    [33] 孙义祥, 郭跃升, 于舜章, 等. 应用“3414”试验建立冬小麦测土配方施肥指标体系[J]. 植物营养与肥料学报, 2009, 15(1): 197–203. Sun Y X, Guo Y S, Yu S Z, et al. Establish phosphorus and potassium fertilization recommendation index based on the “3414” field experiments[J]. Journal of Plant Nutrition and Fertilizers, 2009, 15(1): 197–203. doi:  10.3321/j.issn:1008-505X.2009.01.029
  • [1] 陈日升康文星吕中诚黄志宏赵仲辉邓湘雯 . 杉木林物质生产中养分利用特征分析. 植物营养与肥料学报, 2019, 25(9): 1588-1599. doi: 10.11674/zwyf.18361
    [2] 徐倩宋佳田汇张洋胥婷婷张荣左亚杰马爱生姚致远高亚军 . 青海省春油菜区土壤养分状况及施肥策略. 植物营养与肥料学报, 2019, 25(1): 157-166. doi: 10.11674/zwyf.18040
    [3] 王鹏祝丽香陈香香冯惠孙文帅秦宁 . 桔梗与大葱间作对土壤养分、微生物区系和酶活性的影响. 植物营养与肥料学报, 2018, 24(3): 668-675. doi: 10.11674/zwyf.17325
    [4] 王乐张淑香马常宝李春花 . 潮土区29年来土壤肥力和作物产量演变特征. 植物营养与肥料学报, 2018, 24(6): 1435-1444. doi: 10.11674/zwyf.18187
    [5] 李大明柳开楼叶会财胡志华余喜初徐小林杨旭初周利军胡秋萍胡惠文黄庆海 . 长期不同施肥处理红壤旱地剖面养分分布差异. 植物营养与肥料学报, 2018, 24(3): 633-640. doi: 10.11674/zwyf.17265
    [6] 魏猛张爱君诸葛玉平李洪民唐忠厚陈晓光 . 长期不同施肥对黄潮土区冬小麦产量及土壤养分的影响. 植物营养与肥料学报, 2017, 23(2): 304-312. doi: 10.11674/zwyf.16275
    [7] 郭腾飞梁国庆周卫刘东海王秀斌孙静文李双来胡诚 . 施肥对稻田温室气体排放及土壤养分的影响. 植物营养与肥料学报, 2016, 22(2): 337-345. doi: 10.11674/zwyf.14557
    [8] 潘俊峰万开元李祖章陶勇刘益仁刘毅李志国张过师陈防 . 施肥模式对晚稻田杂草群落的影响. 植物营养与肥料学报, 2015, 21(1): 200-210. doi: 10.11674/zwyf.2015.0122
    [9] 赵亮张贺翠廉小平陆广涛朱利泉 . 喀斯特地形区水稻测土配方施肥指标体系研究. 植物营养与肥料学报, 2015, 21(4): 1056-1065. doi: 10.11674/zwyf.2015.0426
    [10] 徐新朋赵士诚张云贵何萍高强 . 吉林省玉米种植区土壤养分空间变异特征研究. 植物营养与肥料学报, 2011, 17(6): 1342-1350. doi: 10.11674/zwyf.2011.1044
    [11] 邓阳春黄建国 . 长期连作对烤烟产量和土壤养分的影响. 植物营养与肥料学报, 2010, 16(4): 840-845. doi: 10.11674/zwyf.2010.0409
    [12] . 黄土塬面果园土壤养分特征及演变. 植物营养与肥料学报, 2010, 16(5): 1170-1175. doi: 10.11674/zwyf.2010.0518
    [13] 庞夙李廷轩王永东余海英吴德勇 . 土壤速效氮、磷、钾含量空间变异特征及其影响因子 . 植物营养与肥料学报, 2009, 15(1): 114-120. doi: 10.11674/zwyf.2009.0116
    [14] 陈磊郝明德戚龙海 . 长期施肥对黄土旱塬区土壤—植物系统中氮、磷养分的影响. 植物营养与肥料学报, 2007, 13(6): 1006-1012. doi: 10.11674/zwyf.2007.0604
    [15] 吕卫光沈其荣余廷园诸海涛 . 酚酸化合物对土壤酶活性和土壤养分的影响. 植物营养与肥料学报, 2006, 12(6): 845-849. doi: 10.11674/zwyf.2006.0615
    [16] 郭晓敏牛德奎郭熙陈防张过师张斌胡冬南 . 奉新毛竹林土壤养分空间变异性研究. 植物营养与肥料学报, 2006, 12(3): 420-425. doi: 10.11674/zwyf.2006.0322
    [17] 许红卫高克异王珂周斌 . 稻田土壤养分空间变异与合理取样数研究. 植物营养与肥料学报, 2006, 12(1): 37-40. doi: 10.11674/zwyf.2006.0107
    [18] 李新举胡振琪刘宁赵庚星田素峰 . 黄河三角洲土壤肥力质量的时空演变——以垦利县为例. 植物营养与肥料学报, 2006, 12(6): 778-783. doi: 10.11674/zwyf.2006.0604
    [19] 路鹏黄道友宋变兰黄敏刘守龙苏以荣肖和艾吴金水 . 亚热带红壤丘陵区典型景观单元土壤养分的空间变异. 植物营养与肥料学报, 2005, 11(6): 717-723. doi: 10.11674/zwyf.2005.0602
    [20] 金继运 . “精准农业”及其在我国的应用前景. 植物营养与肥料学报, 1998, 4(1): 1-7. doi: 10.11674/zwyf.1998.0101
  • 加载中
图(4)表(3)
计量
  • 文章访问数:  417
  • HTML全文浏览量:  317
  • PDF下载量:  20
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-04-21
  • 网络出版日期:  2021-10-09
  • 刊出日期:  2021-09-25

基于GIS和丰缺指标法的区域施肥管理体系的构建

    作者简介:高宏艳 E-mail:905868022@qq.com
    通讯作者: 索全义, paul98@sina.com
  • 1. 内蒙古农业大学草原与资源环境学院,内蒙古呼和浩特 010018
  • 2. 内蒙古自治区土壤肥料和节水农业工作站,内蒙古呼和浩特 010010
  • 3. 内蒙古自治区农牧业科学院,内蒙古呼和浩特 010037
  • 基金项目: 农业农村部测土配方施肥项目
  • 摘要:   【目的】  采用GIS定位和养分丰缺指标相结合的方法,研究西辽河平原耕地土壤碱解氮、有效磷、速效钾的空间分布特点,建立直观、精准的推荐施肥系统。  【方法】  西辽河平原包括科尔沁区、开鲁县、奈曼旗、科尔沁左翼后旗、科尔沁左翼中旗。以该区域的15421个土壤测试数据和143个玉米“3414试验”数据为基础,运用GIS方法计算耕地土壤碱解氮、有效磷、速效钾的空间插值结果,结合土壤丰缺指标,将西辽河平原耕地土壤碱解氮、有效磷、速效钾养分分区,确定不同土壤氮磷钾养分组合下的施肥量,并制作玉米氮磷钾分区施肥图。  【结果】  西辽河平原耕地土壤养分空间分布不均,大部分土壤碱解氮含量较低,有效磷、速效钾含量中等。耕地土壤有效氮、磷、钾养分含量呈低氮–中磷–中钾 (碱解氮 ≤ 69.5 mg/kg、有效磷4.18~10.40 mg/kg、速效钾90.3~152.5 mg/kg)、中氮–中磷–中钾 (碱解氮69.5~107.9 mg/kg、有效磷4.18~10.40 mg/kg、速效钾90.3~152.5 mg/kg)、中氮–中磷–高钾 (碱解氮69.5~107.9 mg/kg、有效磷4.18~10.40 mg/kg、速效钾 > 152.5 mg/kg) 的面积占比分别为42.4%、14.1%、13.5%。土壤碱解氮、有效磷、速效钾含量 (x) 与经济最佳施肥量 (y) 之间存在较好的对数关系,氮施肥模型y = –102.5ln(x) + 617.22,磷施肥模型y = –36.11ln(x) + 174.1,钾施肥模型为y = –25.89ln(x) + 180.96。不同土壤养分组合种植玉米所需的经济最佳N–P2O5–K2O用量 (kg/hm2) 组合分别为202–110–58、168–107–57、164–102–47。  【结论】  土壤碱解氮、有效磷、速效钾含量与相对产量存在较好的对数函数关系。西辽河灌区土壤有效氮磷钾丰缺指标的划分以75%和85%作为相对产量的分级区间较为合适。在土壤有效氮磷钾养分大数据的支撑下,采用GIS和丰缺指标法相结合的方法,可以准确快捷地构建区域施肥管理体系。

    English Abstract

    • 构建区域施肥管理体系是根据不同地域土壤养分状况进行精确施肥管理的重要途径。根据土壤养分管理分区的情况,科学合理地确定肥料的投入量,不仅可以有效地提高肥料的利用率,促进土壤生产潜力的发挥,而且还能够提高作物产量和改善作物品质[1-4]。近年来,国内外学者针对不同研究区域进行了大量的分区施肥管理研究,在土壤养分管理分区时,采用的方法和分区标准各有不同,未形成统一的方法。以中小尺度为研究对象时,通常采用模糊c-均值聚类法、K均值聚类法等定量分析法[5-6],而针对大尺度的施肥区划多数是基于定性分析[7-8],综合考虑土壤、气候等环境因素,并参考生态分区或种植业分区等区划类型,而缺乏定量分析的理论支撑。研究表明,土壤养分管理分区可以有效区分土壤养分的空间变异性,可以将不同管理分区土壤养分的均值作为变量施肥的依据[9-10]。土壤养分丰缺指标法是区域养分分级和分区的经典方法,适合应用于有大量田间试验数据和土壤测试数据支撑的区域。根据不同作物、不同土壤养分条件确定最佳肥料用量的方法也多种多样,大致可以分为基于土壤养分的推荐施肥方法和基于作物的推荐施肥方法两大类[11]。最经典的肥料推荐方法主要是肥料效应函数法和目标产量法,而肥料效应函数法又是平衡施肥中广泛应用的一种方法,也是养分管理中确定推荐施肥量的常用方法[12-13]

      内蒙古西辽河平原土壤肥沃,耕地面积广,种植的主要粮食作物是玉米,素有“北方粮仓”之称[14],对该区域土壤养分进行分区管理是玉米科学施肥的基础,也是摸清区域土壤养分特征的重要手段,对于农业的可持续发展有重要的意义[15]。但是,由于西辽河平原耕地面积大且地块分散,不能实现每个地块土壤养分的测试。土壤养分在空间分布上表现出一定的空间异质性,而传统的区域推荐施肥方法不能满足空间数据分析的需要,存在一定的局限性,且缺少空间上的直观展示效果。基于这个原因,将GIS方法与土壤养分丰缺指标法、肥料效应函数法等经典方法结合起来,划分土壤养分的分级和制定准确的施肥方案。本研究借助地统计学和GIS手段,将分散的土壤点状养分测试结果,赋值到面域,依托土壤养分的大数据基础,对土壤养分进行分区管理,并以大量的玉米“3414”肥料试验为基础,采用土壤养分丰缺指标法结合GIS方法,构建了不同土壤养分分区的玉米施肥体系,旨在探索科学的区域施肥推荐方法和先进的养分管理技术,为农业生产和科学施肥提供理论依据。

      • 研究中数据为内蒙古西辽河平原2009年15421个土壤样品的碱解氮、有效磷、速效钾测试数据和143个玉米“3414”试验数据,分布于开鲁县 (Kailu Country,KL)、科尔沁区 (Horqin District,HQ)、科尔沁左翼后旗 (简称科左后旗,Horqin Left Back Banner,HLB)、科尔沁左翼中旗 (简称科左中旗,Horqin Left Middle Banner,HLM)、奈曼旗 (Naiman Banner,NM)。以上数据均来自农业农村部测土配方施肥项目 (土壤采样点分布见图1),“3414”试验点分布见图2

        图  1  土壤采样点分布图

        Figure 1.  Distribution map of soil sampling points

        图  2  试验点分布图

        Figure 2.  Distribution map of test points

      • 用Excel软件作回归分析,分别模拟143个玉米“3414”试验的氮、磷、钾肥料用量 (纯量) 和相应处理玉米产量的一元二次效应函数。将符合肥料报酬递减律和一元二次函数曲线特点的方程,进行边际效应分析,采用微分求导法[16]计算经济最佳施肥量和最高产量施肥量。

      • 相对产量是指在不施某种肥料时,作物单位面积的籽粒产量与施用所有肥料时单位面积产量的比值,常用百分数来表示。研究中采用3种不同的方法计算相对产量,通过对比得出最优相对产量。相对产量1 (RY1) 是肥料效应函数中无肥区产量与最经济产量的比值;相对产量2 (RY2) 是肥料效应函数中无肥区产量与最高产量的比值;相对产量3 (RY3) 是无肥区产量与N2P2K2处理区产量的比值。

      • 采用普通克里格插值法 (Ordinary Kriging)[17]对土壤碱解氮、有效磷、速效钾含量进行空间插值,该方法是以变异函数理论和结构分析为基础,通过对已知样本点赋权重来求得未知样本点 (区域) 的值,其计算公式为:$Z\left ( {{\chi _0}} \right) = $$ \displaystyle\sum\limits_{i = {\rm{1}}}^n {{\lambda _i}Z\left ( {{\chi _i}} \right)} $,式中,Z(χ0) 是未知样点χ0处的估计值,Z(χi) 是已知采样点χi处的实测值,λi为第i个已知采样点对未知样点的影响权重系数,n为搜索半径内已知采样点的个数。

      • 应用SPSS软件和Excel软件进行数据的统计分析和散点图的制作。应用ArcMap 10.5软件进行GIS图件制作。

      • 通过玉米“3414”试验点数据,计算经济最佳施肥量和最高产量施肥量,用于对比不同方法计算相对产量的差异以及为推荐区域施肥量提供点状数据。按照等距离法将经济最佳施肥量和最高产量施肥量分别划分为5个等级区间,统计不同等级区间分布的试验数量占比 (表1)。

        表 1  不同产量目标施肥量等级的试验数量占比

        Table 1.  Proportion of test number at different levels of fertilizer

        产量目标
        Yield target
        养分
        Nutrient
        总试验数量
        Total test
        number
        平均施肥量
        Mean rate
        (kg/hm2)
        各级施肥量试验数量占比 (%)
        The proportion of test number under each nutrient rate
        ≤ 5050~100100~150150~200 ≥ 200
        最佳经济产量
        Optimal economical yield
        N137178.002.232.831.433.6
        P2O13098.98.543.038.57.72.3
        K2O12464.022.670.26.40.80
        最高产量
        Maximun yield
        N137209.30019.730.749.6
        P2O130127.70.830.049.211.58.5
        K2O12474.48.979.09.72.40

        表1可知,西辽河平原玉米的经济最佳施N、P2O5、K2O量和最高产量施N、P2O5、K2O量均存在较大变幅,在不同的施肥量等级区间均有分布且分布不均,因此,需要对西辽河平原地区进行分区施肥管理。

      • 以往内蒙古划分土壤养分丰缺指标,将相对产量取值65%定为低级,相对产量取值90%定为高级[18]。经过计算统计,发现相对产量小于70%和大于90%时,所对应的土壤养分含量在实际中分布很少,没有实际划分意义。因此,我们将西辽河平原土壤有效养分的丰缺指标调整为相对产量85%和75%。

      • 分别模拟3种不同方法确定的相对产量与各试验点的土壤碱解氮、有效磷、速效钾含量的对数函数关系。

        选择F值和R2较大的函数关系式用于计算各养分的丰缺指标。将75%、85%分别代入相对产量和土壤养分的函数关系式中,即得到了土壤碱解氮、有效磷、速效钾的丰缺指标 (表2)。由表2可知,碱解氮 ≤ 69.5 mg/kg为低级,69.5~107.9 mg/kg为中级,>107.9 mg/kg为高级。有效磷 ≤ 4.18 mg/kg为低级,4.18~10.40 mg/kg为中级,>10.40 mg/kg为高级。速效钾 ≤ 90.3 mg/kg为低级,90.3~152.5 mg/kg为中级,>152.5 mg/kg为高级。

        表 2  不同相对产量划分的土壤养分丰缺指标

        Table 2.  Soil nutrient abundance index of different relative yields

        土壤养分
        Soil nutrient
        相对产量类型
        Relative
        yield type
        丰缺指标 Abundance index (mg/kg)相对产量与土壤
        养分关系式
        Relative yield and soil
        nutrient function
        R2FP
        低 Low
        ( ≤ 75%)
        中 Medium
        (75%~85%)
        高 High
        (>85%)
        碱解氮
        Alkeline-N
        RY1 ≤ 68.268.2~107.2>107.2y = 22.123ln (x) – 18.4180.49095.057<0.001
        RY2 ≤ 69.269.2~109.3>109.3y = 21.902ln (x) – 17.8010.49295.631<0.001
        RY3 ≤ 69.569.5~107.9>107.9y = 22.716ln (x) – 21.3320.553122.237<0.001
        有效磷
        Available P
        RY1 ≤ 4.094.09~10.14>10.14y = 11.022ln (x) + 59.4670.610140.775<0.001
        RY2 ≤ 4.184.18~10.40>10.40y = 10.985ln (x) + 59.2770.630151.759<0.001
        RY3 ≤ 3.803.80~10.32>10.32y = 10.019ln (x) + 61.6140.611141.495<0.001
        速效钾
        Available K
        RY1 ≤ 90.390.3~152.5>152.5y = 19.062ln (x)–10.8290.606143.137<0.001
        RY2 ≤ 90.990.9~153.9>153.9y = 18.964ln (x)–10.5130.604141.805<0.001
        RY3 ≤ 78.378.3~141.7>141.7y = 16.858ln (x) + 1.48570.542111.111<0.001
        注(Note):RY1 = 无肥区产量/最经济产量; RY2 = 无肥区产量/最高产量; RY3 = 无肥区产量/N2P2K2 处理产量。 函数关系方程中,y 代表相对产量,x 代表土壤养分含量。RY1 = no-fertilizer yield/optimum yield; RY2 = no-fertilizer yield/maximum yield; RY3 = no-fertilizer yield/N2P2K2 yield. In the function, y and x represent relative yield and nutrient content, respectively.
      • 各试验点经济最佳肥料用量与土壤碱解氮、有效磷、速效钾含量存在对数函数关系 (图3)。

        图  3  玉米最佳氮磷钾肥经济施用量与土壤养分含量的函数关系图

        Figure 3.  The function relationship between optimal nutrient rate and soil available nutrient content

        将土壤碱解氮、有效磷、速效钾的丰缺指标所对应的土测值分别代入到图3 N、P2O5、K2O最佳施肥量的施肥模型中,可以得到相应的N、P2O5、K2O的经济最佳肥料用量,根据内蒙古西辽河平原土壤养分的极大值和极小值确定经济合理施肥的区间取值范围 (表3)。

        表 3  不同丰缺指标下玉米最佳N、P2O5、K2O肥料用量

        Table 3.  Maize optimal fertilizer use under different ranges of soil nutrient content

        土壤养分
        Soil nutrient
        丰缺程度
        Grade
        相对产量
        Relative yield
        养分含量 (mg/kg)
        Nutrient content
        经济施肥范围 (kg/hm2)
        Optimal fertilizerr ate
        碱解氮 Alkeline-N低 Low ≤ 75% ≤ 69.5182.5~235.6
        中 Medium75%~85%69.5~107.9137.4~182.5
        高 High > 85%>107.944.4~137.4
        有效磷 Available P低 Low ≤ 75% ≤ 4.18122.4~151.1
        中 Medium75%~85%4.18~10.4089.5~122.4
        高 High > 85%>10.4067.6~89.5
        速效钾 Available K低 Low ≤ 75% ≤ 90.364.4~77.2
        中 Medium75%~85%90.3~152.550.8~64.4
        高 High > 85%>152.539.6~50.8
      • 将西辽河平原15421个土壤样品的碱解氮、有效磷、速效钾测试值进行空间插值。按照土壤碱解氮、有效磷、速效钾的丰缺指标,对西辽河平原耕地土壤的碱解氮、有效磷、速效钾含量水平进行高、中、低等级的赋值。将耕地图与行政区划图叠加,形成耕地土壤碱解氮、有效磷、速效钾含量水平分布图,并在Arcgis10.5软件中进行叠加处理,形成碱解氮、有效磷、速效钾的丰缺水平组合图。根据西辽河平原土壤氮磷钾养分丰缺水平组合分布图,将每个管理单元的土壤碱解氮、有效磷、速效钾含量代入土壤测试值与最佳施肥量的函数关系式,得到每个管理单元玉米的最佳N、P2O5、K2O肥料用量。统计分析内蒙古西辽河平原每种土壤氮磷钾养分丰缺水平组合的平均N、P2O5、K2O肥料用量,将平均值作为该组合的N、P2O5、K2O肥料用量,并用不同的颜色显示,生成了西辽河平原不同土壤养分分区玉米的N、P2O5、K2O最佳肥料用量图 (图4)。

        图  4  西辽河平原土壤养分分区及氮磷钾肥料最佳用量

        Figure 4.  Soil nutrient zoning and the optimal fertilizer rate in Xiliao River Plain

        图4可以看出,西辽河平原耕地土壤养分状况以低氮–中磷–中钾的面积占比最高,为42.4%,即土壤碱解氮含量 ≤ 69.5 mg/kg、有效磷4.18~10.4 mg/kg、速效钾90.3~152.5 mg/kg;其次是中氮–中磷–中钾组合、中氮–中磷–高钾组合、低氮–低磷–低钾组合,耕地面积占比分别为14.1%、13.5%、10.2%,以上4种组合的面积占比之和为80.2%,其他14种组合的面积占比只有19.8%,说明内蒙古西辽河平原耕地的整体养分状况以这4种组合方式为主,且耕地土壤碱解氮、有效磷、速效钾养分含量空间分布不均匀,需要根据耕地土壤的不同养分状况进行合理的推荐施肥。由图4可知,该4种土壤养分组合模式下,种植玉米所需的最佳N、P2O5、K2O用量 (kg/hm2) 组合分别为202-110-58、168-107-57、164-102-47、218-130-68。

      • 采用科学的方法对土壤养分进行分区是养分管理最重要的一个环节。向德明等[19]采用GPS和模糊C-均值聚类相结合的方法,将湘西植烟区划分为3个土壤养分管理区。张云贵等[20]将土壤无机氮、有效磷、速效钾作为数据源,运用聚类分析的方法对研究地块进行了养分管理分区。近年来,我国学者也有针对灌区农作物使用丰缺指标法进行推荐施肥的研究[21-23]。付莹莹等[24]、刘芬等[25]采用丰缺指标法和肥料效应函数法建立了关中灌区夏玉米、冬小麦的丰缺指标体系,并确定了不同肥力水平下的平均推荐施肥量。土壤养分丰缺指标法是国际上最通用的测土配方施肥方法[26-27],也是中国进行测土配方施肥工作所形成的较为成熟的方法之一。本研究中将传统方法与GIS方法相结合,从空间上全方位的将耕地土壤碱解氮、有效磷、速效钾的养分进行了区域划分,并根据养分的空间分布状况确定了氮、磷、钾肥料的用量,具有较好的空间展示效果,实现了土壤养分和推荐施肥的区域化管理。

        在土壤养分丰缺指标的相关研究中,学者们确定丰缺指标分级时采用的分级标准并不统一。在测土配方施肥技术要览中,相对产量小于50%时所对应的养分含量为极低,50%~75%为低,75%~90%为中等,90%~95%为高,大于95%为极高,推荐使用该分级标准进行分级[28]。由于研究区域土壤养分含量水平不同,在划分丰缺指标时会出现缺少极低或极高等级的情况,学者对分级标准 (临界值) 做了不同程度的调整。陈新平等[29]以相对产量50%、75%和95%作为划分土壤养分的丰缺标准,将我国农田的大量元素养分划分为极低、低、中、高4个等级。黄亿等[30]在构建川中丘陵区油菜施肥指标体系时,以相对产量75%、80%、85%、90%为标准,将土壤有效磷的丰缺程度划分为低、较低、中、较高、高5个等级。戢林等[31]以水稻相对产量75%、80%、85%、90%和95%为标准,将川中丘陵区土壤有效养分划分为低、较低、中、较高和高5个等级。本研究依据第二次土壤普查的标准初步判断西辽河平原灌区耕地土壤养分的水平属于中氮中磷中钾,缺少较低和较高含量水平的类型。经统计,在西辽河平原的143个试验点中,与相对产量小于70%和大于90%相对应的土壤养分在实际中的分布很少,因此,把相对产量的分级标准定为75%、85%,将土壤养分丰缺指标划分为低、中、高三级。由于不能充分覆盖土壤养分肥力的分布区间范围,这可能是因为西辽河平原常年大量施用化学肥料,致使耕地土壤养分肥力水平有所提高,玉米的相对产量也有所提高。但是化学肥料的长期施用破坏了土壤结构,使得土壤养分没有出现极高的含量,因此,大于90%的作物相对产量也极少。这个分级标准与其他学者的标准不同,经过验证表明该分级标准适合西辽河平原的土壤养分分级。

        学者们通过相对产量分级建立了不同作物的土壤养分丰缺指标,但是在计算相对产量时采用的方法也不同。大部分学者用无肥区产量与全肥区产量的比值来计算相对产量[30, 32]。也有学者用“3414”试验中无肥区产量与4个处理中最高产量的比值来计算相对产量[33]。本研究中在划分土壤养分丰缺指标时,采用3种不同的方法计算了相对产量,不同方法得到的相对产量与土壤养分的函数模型的模拟效果不同,氮肥、磷肥、钾肥的相对产量分别以RY3、RY2、RY1的拟合效果最优。由于研究对象的不同,试验设计的不同,没有形成统一的相对产量的计算方法。因此,需要根据不同研究对象的特点来确定适宜的计算方法。丰缺指标的高中低分级标准也不统一,应该根据研究区域土壤养分的高低和具体的分布特征,来确定合理的分级标准。

      • 在西辽河平原,以75%和85%的相对产量作为划分土壤有效氮磷钾高中低丰缺指标的临界值更为合理。按照此标准,土壤碱解氮、有效磷、速效钾最主要的组合方式为低氮–中磷–中钾、中氮–中磷–中钾、中氮–中磷–高钾组合,面积占比分别为42.4%、14.1%、13.5%,对应的玉米所需的最佳N、P2O5、K2O用量 (kg/hm2) 组合分别为202–110–58、168–107–57、164–102–47。

    参考文献 (33)

    目录

      /

      返回文章
      返回