应用反求方法模拟冬小麦土柱实验中硝态氮运移.pdfVIP

http :// 应用反求方法模拟冬小麦土柱实验中的硝态氮运移1 石建初，左强* 中国农业大学资源与环境学院土壤与水科学系，北京 (100094) E-mail ：shijianchu@ 摘 要：硝态氮在土壤中的运移转化受到了众多因素的综合影响，其动态过程的模拟一直是 相关研究领域专家学者关注的焦点。本文应用反求方法迭代求解硝态氮运移方程（CDE 方 程）估算其源汇项的平均分布，并进一步优化获得了源汇项中难以直接测定的根系吸氮因子， 从而建立了室内砂培土柱实验中冬小麦的根系吸氮模型。应用该模型模拟冬小麦砂培土柱实 验中各种水分与氮素供应条件下土壤硝态氮动态变化过程的结果表明：冬小麦各生长阶段土 壤硝态氮浓度剖面的模拟值与实测值吻合较好，冬小麦根系吸氮总量模拟值与实测值之间的 相对误差均在 10%以内。 关键词：反求方法 冬小麦 硝态氮 运移 模拟 根系吸氮 中图分类号：2-5-2 1. 引言 为满足日益增长的人口对粮食的需求，全球范围内每年都有大量化肥被投放到土壤-作 物系统中。化肥施用量的不断增长在提高农作物产量的同时，也增加了对土壤、地下水、大 气和地表水环境的污染，其中以硝态氮对地下水的污染最为严重[1]-[3] 。一般植物每年的吸氮 -2 [3] 量介于 50-400Kg hm 之间，而全球范围内氮肥的施用量却远远高于这个水平 。过量使用 氮肥加之不合理的灌溉制度所引起的面源污染是地下水硝酸盐污染的首要原因[4, 5] 。因此如 何在不严重影响农作物产量的情形下减少化肥（尤其是氮肥）的施用量已成为土壤、农业、 环境等领域专家学者们致力研究的焦点问题，准确地模拟预测作物生长条件下土壤硝态氮的 动态变化过程在其中具有举足轻重的作用。 对流-弥散方程（Convection-Dispersion Equation，CDE ）可以有效地描述土壤中硝态氮 的运移转化过程，已被国内外专家学者广泛采用[6]-[9] 。在运用CDE方程模拟土壤硝态氮运移 时，氮素在土壤中的诸多转化过程（包括作物根系吸收、硝化、固持、反硝化等）通常在方 程的源汇项中予以考虑。由于受土壤微环境条件的综合影响，各种土壤氮素转化速率都很难 实际测定，对其准确地定量描述则更为困难[10] 。迄今为止，尽管已有众多的模拟模型用于 描述土壤中氮素的运移转化过程，但建模目标和方式却千差万别，即使对于同一氮素转化过 程，不同学者所建立的模型也大相径庭（如零级反应、一级反应等），由于受到土壤含水量、 氮素浓度、有机质含量、根系生长分布、土壤温度等众多因素的影响，各模型中的参数（如 零级或一级动力学常数、根系吸氮因子等）往往难以准确获得[10]-[12]，目前，多采用试差法 [12] 来优化获取模型中的这些参数，所得参数组合常常不唯一或缺乏严格的数学、物理意义 。 近年来，反求方法越来越受到国内外专家学者的青睐，在估算一些难以准确测定且动态变化 较弱的土壤水力参数（如土壤饱和导水率、水分特征曲线等）方面已得到广泛应用[2, 13, 14] 。 对于动态变化较强的参数（如根系吸水速率、氮素转化速率等），适时适地的瞬时值几乎无 1本课题得到国家自然科学基金项目（项目编号）与高等学校博士学科点专项科研 基金项目（项目编号：20040019001 ）的资助。 -1- http :// 法实测，但与瞬时值相比，其在某一时间段内的平均值可能更有实际意义。基于以上思路， Zuo和Zhang[15]提出了通过两个连续测量的土壤含水量剖面估算平均根系吸水速