传递函数模型的展和农田尺度下硝态氮淋失的数值预报.pdf

摘要 氩索淋必对农业水土环境质量的影响是重要的生态农业研究课题，该研究中颇具挑战性的科 学难题之一是对农田土壤中同时受到物理、化学和生物作用影响的硝态氰淋失动态的精确模拟和 预报。为了对氮索运移和转化行为在生态环境演变过程中的作用有更明晰的认识．本文采用不同 的数学模型分别从点月度和农田尺度对此进行了研究。 本文研究的第一部分通过引入概率密度函数的修正因子，发展了一个可区分表施氮素和土壤 残留氮对硝态氮淋失贞献的传递函数模型。假设累积排水量的概率密度函数的形状不随水流状态 而改变，概率密度函数于不同水流条件下的差异体现在贮水量差的变化，将Jury等(1990)提出 的非稳定流条件下保守溶质运穆的传递函数模型推广到田间土壤硝态氨淋失的模拟中。在模型中 考虑了氮素的矿化、固持、氨挥发和作物吸氮作用，把源汇项对硝态氮淋失的影响视为对溶质运 移体中初始驻留浓度的输入，将各转化量与土壤贮水量或蒸散量建立了较简单的函数形式得到净 矿化、氨挥发和吸氮浓度。以考虑源汇项作用时的硝态氮淋失浓度的均方根差为目标函数，采用 二分法对概率密度函数的修正因子进行优化。使用1999年冬小麦一夏玉米生长条件下土壤剖面2 米埋深处硝态氮淋失动态的试验观钡0数据，检验了所发展的模型的有效性。接着．分析了模型中 各参量对模拟精度的影响，计算了模拟时段内淋失的硝态氮占肥料氰的比例。表明：作者所发展 的能够估计表施氮素和残留氮对硝态氨淋失贡献的传递函数模型的计算结果是合理的；同时考虑 流场变化、残留氮影响和氮素转化作用的传递函数模型具有更高的模拟精度，模拟的累积硝态氮 淋失量的相对误差仅为o．99％：而既不考虑流场变化又不考虑氮素转化作用影响的模拟结果的相 对误差达19．62％；仅考虑源汇项作用的硝态氮淋失量的模拟结果与发展的传递函数模型的模拟结 果相差不大，相对误差为3．03％，说明在北方旱田深排水条件下长作物生育期内硝态氮淋失的模 拟中，考虑氮素的转化作用对模拟结果的影响更大f表施氨的淋失量只占总淋失量的0．42％，占 总施氮量的O．10％，所淋失的硝态氦大部分是土壤的残留氮，占总淋失量的99．58％。表施氮素迁 移时间的概率密度函数的修正因子的取值对模拟精度的影响较大。取零时．模拟淋失量的相对误 差为9．68％，取1时，相对误差为21．61％，取值在0．28到0．32之间时，对于与本文试验条件类 似的农田土壤，应用我们所发展的数学模型可以获得满意的模拟结果。根据本文所发展的能区分 表藏氮和残留氮的传递函数模型，可为田间的1N示踪试验提供一种定量的参考工具。 本文研究的第二部分在开展大规模田间试验的基础上，获得大量具有空间变异特征的模型参 数，对农田尺度下硝态氮的淋失风险进行了数值预报。北京通州区永乐店试验站的小区面积为 27x27m2(网格间距3m)，采集30一50cm土层的样品；山东德州地区禹城试验站的小区面积为 40x40 i112(网格间距5m)，采集o-吧O和20--40cm土层的样品。根据土壤传递函数方法，生成 了各个采样点的土壤水力学参数和弥散度。假设土壤有机氮的矿化速率常数(零级动力学)与有 机质含量成正比得到矿化速率在空间上的分布。基于柱模型假设，即，土壤由一系列不发生相互 麦一夏玉米生长条件下农田尺度下土壤氮素的转化和硝态氮的淋失风险进行了数值模拟。为评价 矿化作用产生的无机氮对感兴趣深度的硝态氮淋失的影响，本文设计了两种数值模拟方案，即， 不考虑和考虑土壤有机氯矿化速率的空间变异性(分别简记为方案1和方案2)。 模拟结果表明：(1)永乐店冬小麦生育期内剖面250cm埋深处平均的土壤水渗透量、累积硝 变异系数均在1．46以上，具有很强的空间变异性；(2)禹城冬小麦生育期内剖面200cm埋深处 量相比，空间上平均的土壤水渗透量和硝态氮淋失量不高，而某些采样点处土壤水渗透量和硝态 氮淋失量太大高于均值；是农田管理和环境监测的脆弱带：(4)两种数值模拟方案对模拟的下边 氮素各转化量的影响更大，由于考虑了氮素矿化作用的空间变异性，方案2中各氮素转化量的空 间变异性大大高于方案1的模拟结果；(5)蒸发量、蒸腾量和蒸散量具有较低的空间变异性 运用地统计学软件GEOPACK对实测土壤性状、模型参数和模拟结果进行地统计学分析，表 明：(1)永乐店采样小区的大部分参数(参量)的半方差可用球状模型表示，其中，容重、饱和 分参数的半方差函数用线性无基台值模型或纯块金模型描述；(2)永乐店冬小麦和夏玉米收获时 剖面250、200、150cm埋深处土壤水遁量的半方差