[发明专利]一种基于区域尺度旱地作物氮投入环境阈值的确定方法

摘要

本发明公开一种基于区域尺度旱地作物氮投入环境阈值的确定方法。该方法包括区域内最高相对产量下施氮量的确定、初步确定区域尺度旱地作物氮投入的环境阈值、初步确定的氮投入环境阈值下的区域产量风险和环境风险评估，最终确定区域尺度旱地作物氮投入环境阈值。本发明确定区域尺度氮投入环境阈值为98％最高相对产量对应的施氮量，其应用十分简单方便。本发明方法确定的氮投入环境阈值是基于区域尺度的历史数据和区域尺度土壤类型的当前试验数据获得，具有区域广泛适应性和可行性，采用试验点进行产量和环境两个风险评估才最终确定，使得实施效果安全可靠。不仅可以维持作物高产，而且可以有效地、大幅度地降低氮肥过量施用导致的环境污染风险。

主权项

1.一种基于黄淮海平原区区域尺度小麦氮投入环境阈值的确定方法，其特征在于包括以下步骤：(1)区域内最高相对产量下施氮量的确定①查阅区域内某一旱地作物的施氮量及其产量的历史文献数据，根据查阅的历史文献施氮量及其作物产量数据，计算每一施氮量对应的作物相对产量，作物相对产量计算为公式1：Y相对产量＝Yi/Ymax×100％   (公式1)，公式1中，Yi为某一施氮梯度下某一施氮量的作物产量，Ymax为该施氮梯度下最大作物产量，Y相对产量为该施氮量下的作物相对产量；②根据区域内的土壤类型设置施氮梯度试验，每个土壤类型至少设置一个施氮梯度试验，每个施氮梯度试验至少设置一个年度的试验，实测各施氮梯度试验中各施氮量对应的作物产量，并计算各施氮梯度试验中每一施氮量对应的作物相对产量，该作物相对产量的计算采用公式2计算；Y试验相对产量＝Y试验i/Y试验max×100％   (公式2)，公式2中，Y试验i为某一年内某一施氮梯度试验某一施氮量下的作物产量，Y试验max为该年内该施氮梯度试验最大的作物产量，Y试验相对产量为该年内该施氮梯度试验该施氮量下的作物相对产量；③以步骤(1)①和步骤(1)②所述的各施氮量为横坐标，以步骤(1)①和步骤(1)②所述的各施氮量对应的作物相对产量为纵坐标，用数据统计软件，建立区域尺度上施氮量与作物相对产量的产量效应方程，所建立的区域尺度上施氮量与作物相对产量的产量效应方程为公式3：Y相对产量＝a3x2+b3x+c3   (公式3)，公式3中，Y相对产量为作物相对产量，x为施氮量，a3和b3为系数，c3为常数；④根据公式3的顶点坐标公式：〔(‑b3/2a3),(4a3c3‑b32)/4a3)〕，求得区域内最高相对产量下施氮量X最高相对产量下施氮量和最高相对产量Y最高相对产量，X最高相对产量下施氮量＝(‑b3/2a3)，Y最高相对产量＝(4a3c3‑b32)/4a3；(2)初步确定区域尺度旱地作物氮投入的环境阈值①以步骤(1)④计算的最高相对产量Y最高相对产量为基准，通过公式4，计算从99％至94％每降低一个百分点的相对产量，将以最高相对产量Y最高相对产量为基准，从99％至94％每降低一个百分点的相对产量称为不同产量保证率，所述不同产量保证率分别为99％最高相对产量、98％最高相对产量、97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量，W＝Y最高相对产量×E   (公式4)，公式4中E分别为0.99、0.98、0.97、0.96、0.95、0.94，W为不同产量保证率，E分别为0.99、0.98、0.97、0.96、0.95、0.94时，W分别对应地为99％最高相对产量、98％最高相对产量、97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量；将计算得到的99％最高相对产量、98％最高相对产量、97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量分别代入公式3中计算得到不同产量保证率下对应的施氮量；②绘制不同产量保证率下边际土壤无机氮储量图A、建立施氮量与土壤无机氮储量的关系方程实测步骤(1)②各施氮梯度试验中的每一个施氮量下在0‑100cm土壤中的土壤无机氮储量，以各施氮梯度试验的每一个施氮量为横坐标，每一个施氮量对应的土壤无机氮储量为纵坐标，用数据统计软件，建立施氮量与土壤无机氮储量的关系方程，如果施氮梯度试验为两年以上的试验，则纵坐标为每一个施氮量对应的两年以上的平均土壤无机氮储量，施氮量与土壤无机氮储量的关系方程为公式5：Y土壤无机氮储量＝a5eb′x   (公式5)，公式5中，Y土壤无机氮储量为0‑100cm土壤中的土壤无机氮储量，单位为kg·hm‑2；x为施氮量、单位为kg，a5为常数、b′为系数；分别将步骤(2)①计算得到的99％最高相对产量、98％最高相对产量、97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量对应的施氮量代入公式5，计算得到对应的土壤无机氮储量Y土壤无机氮储量；B、绘制不同产量保证率下边际土壤无机氮储量图边际土壤无机氮储量采用公式6计算：边际土壤无机氮储量y＝土壤无机氮储量(y+1)‑土壤无机氮储量y   (公式6)，公式6中，土壤无机氮储量y分别为99％最高相对产量、98％最高相对产量、97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量对应的土壤无机氮储量，土壤无机氮储量(y+1)为比土壤无机氮储量y对应的相对产量高一个百分点的相对产量对应的土壤无机氮储量；在横坐标上从左至右依次取如下不同产量保证率：99％最高相对产量，98％最高相对产量，97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量，纵坐标为不同产量保证率对应的边际土壤无机氮储量，绘制不同产量保证率下边际土壤无机氮储量图；③绘制不同产量保证率下边际氨挥发量图A、建立施氮量与氨挥发量的关系方程实测步骤(1)②各施氮梯度试验中的每一个施氮量下的氨挥发量，以各施氮梯度试验的每一个施氮量为横坐标，每一个施氮量对应的氨挥发量为纵坐标，用数据统计软件，建立施氮量与氨挥发量的关系方程，如果施氮梯度试验为两年以上的试验，则纵坐标为每一个施氮量对应的两年以上的平均氨挥发量，施氮量与氨挥发量的关系方程为公式7：Y氨挥发量＝b7+k7x   (公式7)，公式7中，Y氨挥发量为氨挥发量、单位为kg N·hm‑2；x为施氮量、单位为kg，k7为系数、b7为参数；分别将步骤(2)①计算得到的99％最高相对产量、98％最高相对产量、97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量对应的施氮量代入公式7，计算得到对应的氨挥发量Y氨挥发量；B、绘制不同产量保证率下边际氨挥发量图边际氨挥发量采用公式8计算：边际氨挥发量y＝氨挥发量(y+1)‑氨挥发量y   (公式8)，公式8中，氨挥发量y分别为99％最高相对产量、98％最高相对产量、97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量对应的氨挥发量，氨挥发量(y+1)为比氨挥发量y对应的相对产量高一个百分点的相对产量对应的氨挥发量；在横坐标上从左至右依次取如下不同产量保证率：99％最高相对产量，98％最高相对产量，97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量，纵坐标为不同产量保证率对应的边际氨挥发量，绘制不同产量保证率下边际氨挥发量图；④绘制不同产量保证率下边际氮淋溶流失量图A、建立施氮量与氮淋溶流失量的关系方程实测步骤(1)②各施氮梯度试验中的每一个施氮量下氮淋溶流失量，以各施氮梯度试验的每一个施氮量为横坐标，每一个施氮量对应的氮淋溶流失量为纵坐标，用数据统计软件，建立施氮量与氮淋溶流失量的关系方程，如果施氮梯度试验为两年以上的试验，则纵坐标为每一个施氮量对应的两年以上的平均氮淋溶流失量，施氮量与氮淋溶流失量的关系方程为公式9：Y氮淋溶流失量＝a9x2+b9x+c9(公式9)，公式9中，Y氮淋溶流失量为氮淋溶流失量，单位为kg N·hm‑2；x为施氮量，单位为kg，a9和b9为系数，c9为常数；分别将步骤(2)①计算得到的99％最高相对产量、98％最高相对产量、97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量对应的施氮量代入公式9，计算得到对应的氮淋溶流失量Y氮淋溶流失量；B、绘制不同产量保证率下边际氮淋溶流失量图边际氮淋溶流失量采用公式10计算：边际氮淋溶流失量y＝氮淋溶流失量(y+1)‑氮淋溶流失量y(公式10)，公式10中，氮淋溶流失量y分别为99％最高相对产量、98％最高相对产量、97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量对应的氮淋溶流失量，氮淋溶流失量(y+1)为比氮淋溶流失量y对应的相对产量高一个百分点的相对产量对应的氮淋溶流失量；在横坐标上从左至右依次取如下不同产量保证率：99％最高相对产量，98％最高相对产量，97％最高相对产量、96％最高相对产量、95％最高相对产量、94％最高相对产量，纵坐标为不同产量保证率对应的边际氮淋溶流失量，绘制不同产量保证率下边际氮淋溶流失量图；⑤比较分析步骤(2)②绘制的不同产量保证率下边际土壤无机氮储量图、步骤(2)③绘制的不同产量保证率下边际氨挥发量图、步骤(2)④绘制不同产量保证率下边际氮淋溶流失量图，该三幅图中99％最高相对产量至98％最高相对产量这一段各环境指标下降幅度均大于其它任意相邻的两个不同产量保证率下各环境指标的下降幅度，初步确定98％最高相对产量对应的施氮量为初步确定的区域尺度旱地作物氮投入环境阈值；所述初步确定的区域尺度旱地作物氮投入环境阈值为将98％最高相对产量代入步骤(1)③所述的区域尺度上施氮量与作物相对产量的产量效应方程公式3：Y相对产量＝a3x2+b3x+c3中计算得到的x值；所述环境指标分别为边际土壤无机氮储量、边际氨挥发量、边际氮淋溶流失量；(3)初步确定的氮投入环境阈值下的区域产量风险和环境风险评估，最终确定区域尺度旱地作物氮投入环境阈值①产量风险评估A、分别对步骤(1)②每个施氮梯度试验建立每一年的施氮梯度试验施氮量与作物实际产量的产量效应方程，每个施氮梯度试验每一年的施氮梯度试验的施氮量与作物实际产量的产量效应方程均按以下方法建立：以步骤(1)②所述某一施氮梯度试验某一年的各施氮量为横坐标，以实测的该各施氮量对应的作物实际产量为纵坐标，用数据统计软件，建立某一施氮梯度试验某一年的施氮梯度试验的施氮量与作物实际产量的产量效应方程为公式11：Y梯度试验实际产量＝a11x2+b11x+c11   (公式11)，公式11中，Y梯度试验实际产量为某一施氮梯度试验某一年的作物实际产量，x为该施氮梯度试验该年的施氮量，a11和b11为系数，c11为常数；B、将步骤(2)⑤初步确定的区域尺度旱地作物氮投入环境阈值代入公式11，计算出各施氮梯度试验各年度在初步确定的区域尺度旱地作物氮投入环境阈值下的作物实际产量Y梯度试验实际产量，根据公式12计算各施氮梯度试验各年度作物产量波动率：公式12中，y_i为任一施氮梯度试验某一年在初步确定的区域尺度旱地作物氮投入阈值下的作物实际产量，为该年各施氮梯度试验在初步确定的区域尺度旱地作物氮投入环境阈值下的平均作物实际产量；C、绘制作物产量波动率柱状图在横坐标上设置各梯度试验，纵坐标为作物产量波动率，绘制各施氮梯度试验各年度作物产量波动率柱状图，比较各施氮梯度试验的作物产量波动率和减产2％的作物产量波动率的大小，当绘制的各施氮梯度试验的产量波动率大于减产2％时的作物产量波动时，确定步骤(2)⑤以98％最高相对产量对应的施氮量为初步确定的区域尺度旱地作物氮投入环境阈值为无产量风险；②环境风险评估A、将步骤(2)⑤初步确定的区域尺度旱地作物氮投入环境阈值分别代入施氮量与土壤无机氮储量的关系方程公式5：Y土壤无机氮储量＝a5eb′x、施氮量与氨挥发量的关系方程公式7：Y氨挥发量＝b7+k7x、施氮量与氮淋溶流失量的关系方程公式9：Y氮淋溶流失量＝a9x2+b9x+c9，分别计算出初步确定的区域尺度旱地作物氮投入环境阈值下的土壤无机氮储量、氨挥发量、氮淋溶流失量；B、将步骤(1)④通过公式3的顶点坐标公式求得的区域内最高相对产量下施氮量X最高相对产量下施氮量分别代入公式5：Y土壤无机氮储量＝a5eb′x、公式7：Y氨挥发量＝b7+k7x、公式9：Y氮淋溶流失量＝a9x2+b9x+c9，分别计算出区域内最高相对产量下施氮量对应的土壤无机氮储量、氨挥发量、氮淋溶流失量；C、比较氮的潜在流失风险，当初步确定的区域尺度旱地作物氮投入环境阈值下的土壤无机氮储量与区域内最高相对产量下施氮量对应的土壤无机氮储量差异显著，同时初步确定的区域尺度旱地作物氮投入环境阈值下的氮流失总量与区域内最高相对产量下施氮量对应的氮流失总量差异显著时，确定步骤(2)⑤以98％最高相对产量对应的施氮量为初步确定的区域尺度旱地作物氮投入环境阈值为无环境风险；所述氮流失总量＝氨挥发量+氮淋溶流失量；③当产量风险评估确定初步确定的区域尺度旱地作物氮投入环境阈值为无产量风险和环境风险评估确定初步确定的区域尺度旱地作物氮投入环境阈值为无环境风险时，则步骤(2)⑤以98％最高相对产量对应的施氮量即为基于区域尺度旱地作物氮投入环境阈值；所述旱地作物为小麦。