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

传递函数模型的发展和农田尺度下硝态氮淋失的数值预报氮素淋失对农业水土环境质量的影响是重要的生态农业研究课题，该研究中颇具挑战性的科学难题之一是对农田土壤中同时受到物理、化学和生物作用影响的硝态氮淋失动态的精确模拟和预报。为了对氮素运移和转化行为在生态环境演变过程中的作用有更明晰的认识，本文采用不同的数学模型分别从点尺度和农田尺度对此进行了研究。本文研究的第一部分通过引入概率密度函数的修正因子，发展了一个可区分表施氮素和土壤残留氮对硝态氮淋失贡献的传递函数模型。假设累积排水量的概率密度函数的形状不随水流状态而改变，概率密度函数于不同水流条件下的差异体现在贮水量差的变化，将Jury等（1990）提出的非稳定流条件下保守溶质运移的传递函数模型推广到田间土壤硝态氮淋失的模拟中。在模型中考虑了氮素的矿化、固持、氨挥发和作物吸氮作用，把源汇项对硝态氮淋失的影响视为对溶质运移体中初始驻留浓度的输入，将各转化量与土壤贮水量或蒸散量建立了较简单的函数形式得到净矿化、氨挥发和吸氮浓度。以考虑源汇项作用时的硝态氮淋失浓度的均方根差为目标函数，采用二分法对概率密度函数的修正因子进行优化。使用1999年冬小麦—夏玉米生长条件下土壤剖面2米埋深处硝态氮淋失动态的试验观测数据，检验了所发展的模型的有效性。接着，分析了模型中各参量对模拟精度的影响，计算了模拟时段内淋失的硝态氮占肥料氮的比例。表明：作者所发展的能够估计表施氮素和残留氮对硝态氮淋失贡献的传递函数模型的计算结果是合理的；同时考虑流场变化、残留氮影响和氮素转化作用的传递函数模型具有更高的模拟精度，模拟的累积硝态氮淋失量的相对误差仅为0.99％；而既不考虑流场变化又不考虑氮素转化作用影响的模拟结果的相对误差达19.62％；仅考虑源汇项作用的硝态氮淋失量的模拟结果与发展的传递函数模型的模拟结果相差不大，相对误差为3.03％，说明在北方旱田深排水条件下长作物生育期内硝态氮淋失的模拟中，考虑氮素的转化作用对模拟结果的影响更大；表施氮的淋失量只占总淋失量的0.42％，占总施氮量的0.10％，所淋失的硝态氮大部分是土壤的残留氮，占总淋失量的99.58％。表施氮素迁移时间的概率密度函数的修正因子的取值对模拟精度的影响较大，取零时，模拟淋失量的相对误差为9.68％，取1时，相对误差为21.61％，取值在0.28到0.32之间时，对于与本文试验条件类似的农田土壤，应用我们所发展的数学模型可以获得满意的模拟结果。根据本文所发展的能区分表施氮和残留氮的传递函数模型，可为田间的¹⁵N示踪试验提供一种定量的参考工具。本文研究的第二部分在开展大规模田间试验的基础上，获得大量具有空间变异特征的模型参数，对农田尺度下硝态氮的淋失风险进行了数值预报。北京通州区永乐店试验站的小区面积为27×27m²（网格间距3m），采集30—50cm土层的样品；山东德州地区禹城试验站的小区面积为40×40m²（网格间距5m），采集0—20和20—40cm土层的样品。根据土壤传递函数方法，生成了各个采样点的土壤水力学参数和弥散度。假设土壤有机氮的矿化速率常数（零级动力学）与有机质含量成正比得到矿化速率在空间上的分布。基于柱模型假设，即，土壤由一系列不发生相互作用的一维土柱组成，运用HYDRUS-1D软件，对永乐店2000—2001年和禹城2001—2003冬小麦—夏玉米生长条件下农田尺度下土壤氮素的转化和硝态氮的淋失风险进行了数值模拟。为评价矿化作用产生的无机氮对感兴趣深度的硝态氮淋失的影响，本文设计了两种数值模拟方案，即，不考虑和考虑土壤有机氮矿化速率的空间变异性（分别简记为方案1和方案2）。模拟结果表明：（1）永乐店冬小麦生育期内剖面250cm埋深处平均的土壤水渗透量、累积硝态氮淋失量分别为2.25mm、0.0098m沙m，，夏玉米生育期内的值分别为l.35mm、o.oo65m沙m，，变异系数均在1.46以上，具有很强的空间变异性;（2）禹城冬小麦生育期内剖面ZO0cm埋深处平均的土壤水渗透量、累积硝态氮淋失量分别为1930mm、0.0790m沙mZ，夏玉米生育期内的值分别为24.48mm、0.0948m沙m，，变异系数均在0.50以上，表现为中等变异强度;（3）与总输入量相比，空间上平均的土壤水渗透量和硝态氮淋失量不高，而某些采样点处土壤水渗透量和硝态氮淋失量大大高于均值，是农田管理和环境监测的脆弱带;（4）两种数值模拟方案对模拟的下边界处（永乐店:250cm;禹城:ZO0cm）硝态氮淋失的影响很小，其差异主要在于方案2对土壤氮素各转化量的影响更大，由于考虑了氮素矿化作用的空间变异性，方案2中各氮素转化量的空间变异性大大高于方案1的模拟结果;（5）蒸发量、蒸腾量和蒸散量具有较低的空间变异性（CV<0.06），而剖面不同埋深处的水通量（250、200、150、IO0cm）呈现较强的空间变异性。运用地统计学软件GEOPACK对实测土壤性状、模型参数和模拟结果进行地统计学分析，表明:（l）永乐店采