黄微尘，周 丰**，梁 浩，陈 磊

(1:北京大学城市与环境学院，北京 100871)(2:北京大学地表过程分析与模拟教育部重点实验室，北京 100871)(3:河海大学农业科学与工程学院，南京 210098)(4:北京师范大学环境学院，北京 100875)

稻米是全球超过一半人口的主食[1]，稻田分布在全球120多个国家，但90%的播种面积集中在东亚、南亚和东南亚[2]. 2017年，稻田消耗了全球50%的灌溉量和14%的化肥施用量[3]，该比例在中国更高，超过60%和20%[2]. 这种粗放的水肥管理，再加上气候变化，造成了过量的氮磷径流流失，严重威胁区域水安全[4]. 有研究表明，中国稻田氮磷径流流失自1990年以来分别增加了46%和30%[5]，且未来这一风险将进一步提高[6]. 然而，水稻品种类型繁多、种植制度多样、水肥管理迥异，稻田氮磷径流流失过程十分复杂，未来时空格局及其生态环境效应还存在较大的不确定性.

稻田氮磷径流流失的影响因素众多，包括降水特征(强度、雨型、持续时间等)[7]、土壤理化性质[8]、农田管理措施[9]、肥料管理(强度、类型、方法和时间)[10]、水分管理(灌溉方式、田埂排水高度)[11]等. 同时，稻田氮磷径流流失过程还与土壤养分循环的其他过程密切相关(如作物吸收、淋溶、挥发等)，完全通过观测手段来解析稻田氮磷径流流失的区域规律是不现实的. 因此，提高稻田氮磷径流流失过程模拟能力，是优化稻田水肥管理和实现气候智能型稻田的关键手段.

自1950s以来，研究者陆续开发了适用于稻田氮磷径流流失的多种模型，并利用不同区域、尺度的观测和控制试验进行了模型校准和验证. 这些模型大体上分为4类，包括统计经验模型、物理机理模型、水文过程模型和生态系统过程模型，从简单的因果分析到复杂的机理探索，从单一模块的开发到多模块的耦合，已经具备了田块尺度、单一场次降雨到区域尺度连续降雨的多情景预测功能. 目前，已有不少相关方面的综述文章，如梳理SWAT(soil and water assessment tool)、DNDC(denitrification-decomposition)等某一个具体模型的发展和应用历程[12-13]，或是聚焦面源污染、水文过程、作物生长等某一类模型的方法研发和改进[14-16]，也有着眼于降雨-径流、污染物迁移转化、氮磷循环等具体的过程机理及影响因素[17-19]. 模型开发和应用中具体的环节，如数据获取、参数率定、不确定性分析[20-22]，也是综述文章的一个重要方面. 然而，系统总结稻田氮磷径流流失模型研发历程，尤其是其中降雨引起的释放和输移等过程以及农艺管理措施等方面的综述仍较少，也尚未系统归纳模型的关键问题或技术难点.

鉴于此，本文将系统阐释稻田氮磷径流流失发生机制，梳理国内外稻田径流流失4类模型的研发历程，对比分析主流模型的性能特征，归纳模型在发生机制参数化、模型参数合理化以及农艺管理措施模块化等关键问题，在此基础上提出模型研究展望，为稻田氮磷径流流失模型的改进和应用提供科学依据.

1 稻田氮磷径流流失发生机制

与旱地不同的是，稻田在水稻生长季会经历频繁的干湿交替过程，为水稻和土壤微生物提供了独特的生境[23]，发生了作物吸收、肥料水解、有机质矿化、养分溶解态和吸附态互相转化等过程. 由于特殊的淹水环境，稻田土壤呈现还原性，有利于反硝化作用的产生. 同时，田面水溶解氧较高，水稻茎秆通气较好，使得水土界面又形成氧化层，促进硝化作用的发生[19]. 因此，水稻种植具有与稻田环境相应、区别于旱地的田间措施、水肥管理制度和种植方式[24].

稻田氮磷径流流失发生机制总结如图1. 稻田的氮磷输入主要包括干湿沉降、灌溉、施肥和生物固氮. 当发生降雨时，雨滴携带氮磷进入田面水中，改变田面水水位及氮磷浓度[5,25]. 与旱地直接产生土壤溅蚀和颗粒分离不同，稻田田面水层对于降雨的冲击力还具有缓冲作用，降低了滴溅侵蚀能力[26]. 雨水经过水稻冠层截流之后，剩余的雨滴能量在田面水表层转化成势能，产生间歇性的压力扰动，造成水土界面土壤孔隙水的释放[7]. 土壤中存在氮磷的垂向扩散，同时由于重力作用，土壤表层溶质也会向下渗透. 随着雨量的增大，土壤淋溶损失逐渐增大至饱和，当降水量大于渗透率，田面水位开始抬升并超过田埂，产生地表径流[27]. 孔隙水、田面水和降雨灌溉水混合后形成水平推流，经过田间的输移过程，以溶解态氮素和吸附态磷素为主到达排水口，排出稻田形成地表径流流失，这一过程中部分氮磷又通过吸附解析作用再次沉淀.

由于暴雨事件往往具有瞬时性和短暂性，稻田氮磷径流流失对于初始条件敏感性较高. 初始田面水位、田面水浓度、土壤浓度等会显著影响到氮磷流失量[7]. 因此，施肥量、施肥方式以及灌排制度通过改变稻田初始水量和养分含量影响流失过程. 其中，肥料在施用之后，短时间内能显著提高田面水氮磷浓度，增加流失风险[28]. 同时土壤中氮磷三态的相互转化、作物吸收、蒸散发、氨挥发、反硝化造成的气体排放、干沉降、淋溶、侧渗等过程都会通过间接途径参与整个稻田系统的径流流失. 由此可见，有效模拟稻田径流流失，需要准确量化径流流失过程、稻田初始条件和水肥管理制度.

图1 稻田氮磷径流流失发生机制与主要过程

2 模型研发历程

2.1 统计经验模型

统计经验模型是基于观测数据，利用因果分析和统计分析方法，建立的不同形态氮磷流失通量与降雨量、施肥量、灌溉量等因素的响应方程，但难以考虑稻田复杂的水文和氮磷迁移转化过程(图2). 例如，SCS模型(soil conservation service)[29]根据土壤和降雨因素来确定径流总量，包含一个代表土壤通透性的重要综合参数CN(curve number)，可适应多种土壤类型和土地利用类型下的径流量模拟. 随后，研究者通过实验手段调整不同稻区的具体数值[30-31]. Lian等[32]搜集了55个研究地点的大量降雨径流数据，建立了CN-China查找表，提供了适用于中国的参考指南. 另外，输出系数模型[33]是根据田间实验确定污染物输出系数，应用于稻田年尺度氮磷径流流失负荷估算. SPARROW模型(spatially referenced regressions on watershed attributes)[34]将统计学方法应用于污染物的质量传输过程，为大尺度应用提供了一种体现区域差异的统计方法.

由于统计经验模型只考虑输入和输出之间的关系，不涉及具体的过程和机理，对数据质量要求较低，适合于大尺度、年尺度的氮磷径流流失负荷估算. 例如，我国全国第一次污染源普查就采用输出系数法对全国尺度稻田的污染物流失系数进行核算. Hou等[25]考虑施氮强度、气候、灌溉量和土壤属性等环境因子对氮流失率的影响，建立了农田总氮径流流失量与施氮强度之间的非线性响应模式(附表Ⅰ和Ⅱ)，更新了中国第一次污染普查的结果. 同时，统计经验模型也多用来探索多输入变量之间的相互关系，来获取稻田氮磷径流流失过程的关键控制因子，识别有针对性的减排措施[5,33]. 当数据存在缺漏时，还可应用统计模型建立回归关系补全数据序列. 由于其操作简便，所需数据资料少，常纳入其他模型的模块中，或结合简单的机理模型提高其区域模拟的能力[35]，但模型的系数多为经验参数，不具备实际物理意义，所以应用于不同区域时常需要进行修正[30]，这导致模拟与试验地点自然、人为因素显著不同的区域时精度难以保障.

图2 稻田氮磷径流流失的4类模型研发历程

2.2 物理机理模型

物理机理模型比统计经验模型复杂，考虑了稻田输入、输出和中间迁移转化过程机理，属于“白箱”模型. 物理机理模型以物质守恒、能量守恒定律为基础，具有严格理论过程，模型参数有明确的物理意义，求解方法一般是基于设定的初始和边界条件进行数值求解(图2).

物理机理模型量化了水土界面氮磷迁移转化过程. 这一过程包括了降雨引起的物理侵蚀作用、化学物质的迁移转化、土壤水-径流混合推移等. 传统的模型主要包括两类，即混合层模型[36-37]和扩散模型[38-39]. 混合层模型假设土壤表层存在一个较薄的混合层，降水、径流与土壤水在该混合层完全混合，而扩散模型假定扩散过程控制土壤和径流之间的化学传递. 这两类模型自提出以来，关键公式和基本假设不断得到完善，但往往认为某一种过程起主导性的作用，且尚未考虑到这些过程的交互影响，尤其对于降雨的作用考虑不全. Gao等[40-41]结合了降雨主导的土壤水释放过程和扩散主导的深层土壤迁移到混合层(交换层)的过程，其中，降雨控制积水-径流层和交换层间的溶质迁移，扩散控制土壤层和交换层间的迁移(附表Ⅰ和Ⅱ). 理论上，这些模型可以应用于具有田面积水和田埂的稻田[42]，但田面水位对降雨溅蚀作用的缓冲作用仍然没有得到很好的模拟. Higashino等[7]从另一个角度出发，将雨滴大小、降落速度、雨滴密度(即单位面积的雨滴个数)等参数纳入模型，细化了稻田水-土界面养分释放过程. 结果表明，该释放过程对于稻田氮径流流失的贡献超过20%. 尽管物理机理模型的科学假设不断完善，机理过程逐渐完整，模拟精度和分辨率逐步提高，但输入数据和边界条件要求高，模拟过程相对复杂[43]，一般是基于少量的室内实验或田间试验进行参数率定和不确定性分析，难以在大尺度推广应用.

2.3 水文过程模型

水文过程模型在面源污染方面应用很广，涉及到降雨-径流、土壤侵蚀、污染物迁移转化等过程(图2)[44]，从而实现不同气象条件、土壤属性和管理措施的综合模拟(附表Ⅰ和Ⅱ). 水文过程模型不仅考虑不同土壤类型之间的差异，还常对土壤进行分层分区，以提高模型模拟能力. 从19世纪60年代开始，研究者陆续开发出很多适用于农田和流域尺度的水文过程模型. 农田尺度模型用于单一的土地利用和相对均一的土壤质地，是大尺度分布式参数模型建模的基础，如CREAMS(chemicals, runoff, and erosion from agricultural management systems)[45]和GLEAMS(groundwater loading effects of agricultural management systems)[46]；流域尺度模型考虑了流域内气候、土地利用、土壤类型、管理措施等的空间异质性，得到整个流域的污染物输出量，如SWAT(soil and water assessment tool)[47]和MIKE SHE(MIKE system hydrological European)[48]. 这些模型虽然为过程模型，但仍然存在很多经验的公式，如 SWAT中直接采用了SCS和USLE模型计算单个响应单元的径流和土壤侵蚀(附表Ⅰ). 同时，这些模型无法准确模拟稻田氮磷径流流失特征，如GLEAMS模型模拟渗漏时采用的Storage-Routine方法、SWAT将稻田田面水层看作蓄水池的处理等. 鉴于此，研究者不断对已有的水文模型进行改进使之具备稻田模拟能力. 例如，RICEWQ(rice water quality)[49]、PCPF-1(pesticide concentration in paddy field, v1)[50]模型考虑氮磷在田面水中的过程，CREAMS-PADDY[51]、GLEAMS-PADDY[52]模型考虑了稻田土壤中的迁移转化行为，GLEAMS-PR(GLEMAS paddy-rice)[53]和RICEWNB(RICE-water and nitrogen balance)[54]模型还探索了稻田传热和根区土壤氮循环的模拟.

计算机和3S(GIS、GPS、RS)技术的不断发展提供了大量的空间数据和统计数据，模型的数据处理能力和应用价值大大提升. 水文过程模型在此背景下从单纯的数值计算程序成为了多功能一体的专业软件系统，更多地应用于大尺度稻区的高分辨模拟，包括氮磷径流流失通量估算和减排措施优化. Yan等[55]在WALRUS模型中考虑了稻田的灌溉和排水方式，引入了新的磷模块，形成了面向圩区的水文过程模型PHPS(lowland polder hydrology and phosphorus modelling system)，对太湖流域稻区有无抽水情景进行模拟，并提出了增加抽水、频繁清淤以及提高稻区范围水体水生植物的覆盖范围等关键减排措施.

水文过程模型具备完善的水文模块和氮磷循环模块，适用于稻区水分和养分循环研究，探索流域氮磷流失规律. 这些水文过程模型全面考虑田面水的推流过程以及氮磷在土壤层中的物理化学迁移转化过程，但多用于日或月尺度模拟，且常将降雨仅仅视为水量输入而忽略了降雨对土壤的间歇性冲击扰动作用，使得稻田氮磷径流流失的模拟结果存在一定偏差. 另外，水文过程模型对农艺管理措施的参数化较为简单，对于大尺度模拟存在较大的不确定性.

2.4 生态系统模型

生态系统模型动态描述作物生长及水、氮和磷元素在土壤-作物-大气连续体中的生物地球化学循环过程(图2). 与水文过程模型不同，其更侧重于水稻的生理生态和生长发育过程，以及对氮磷输入、气象和水肥管理的响应. 农田生态系统模型开发起始于1960s-1970s[56-59]，如美国学者开发了CERES 模型并整合到DSSAT系统[60-61]. 部分农田生态系统模型针对稻田进行了改进，如CERES-Rice[62]和WHCNS-Rice[63]. 李长生等[64-65]构建了DNDC模型(denitrification-decomposition)，较好地模拟农业生态系统碳-氮-水的生物地球化学循环. 与这类模型相比，陆面过程模型引入了农作物模块[66]，更好体现作物生理生态机制[67]，例如ORCHIDEE模型(organizing carbon and hydrology in dynamic ecosystems)[68]. 王旭辉[69]在此基础上开发了ORCHIDEE-crop模型以及其参数最优化系统，实现不同稻区水稻模拟.

生态系统模型具备多时空尺度模拟能力，如田块尺度的EPIC模型[70]发展起来的APEX模型(agricultural policy/environmental extender)[71]可用于小流域尺度，该模型时空分辨率高，具有较好的水稻种植管理措施模块，已具备响应气候变化和稻田农艺管理措施的模拟能力[72]. 目前，生态系统模型已经广泛应用于稻田模拟，多用于不同情景下水稻生长、水热通量和温室气体的排放模拟，然而，氮磷径流流失模拟能力较弱. Zhao等[73]将SCS和RULSE引入DNDC模型(附表Ⅰ)[74]，实现稻田氮素径流流失的季节过程模拟，并应用于上海地区稻田最佳水肥管理措施优化. 相比物理机理模型和水文过程模型，生态系统模型更关注养分循环和作物生长发育，对稻田径流流失的模拟具有浓厚的经验性质，无法有效响应气候变化和水肥管理的影响[75]，难以实现区域普适性[76]. 同时，由于磷的迁移转化能力相对较弱，也不参与稻田温室气体的排放，生态系统模型往往忽略了磷循环过程.

稻田氮磷径流流失发生机制十分复杂，常常与许多因素(如气象水文、土地利用、土壤类型、农田管理措施等)密切相关. 以上4类模型侧重点不同，模拟氮磷径流流失过程采用的方法和公式表达存在差异，对于驱动数据的要求和适用的范围尺度等也有所区别，这导致了模拟功能、时空尺度和适应范围方面上的差别. 本文从稻田氮磷径流流失的过程环节出发，选取7个具有代表性的流失模型进行归纳和对比(附表Ⅰ和Ⅱ). 模型的开发及具体应用可参阅相关研究成果.

2.5 不确定性评估

稻田氮磷径流流失具有复杂的发生机制，受到各类随机因素的影响，而模型运用相对简单的数学公式进行表述，就会造成不确定性. 不确定性主要可以分为随机性、模糊性、灰色性和未确知性[77]. 不确定性研究旨在理清不确定性来源并试图降低不确定性，提高模型预测能力. 不确定性分析方法和理念随着稻田流失模型的改进也在不断发展. 但模型的不确定性十分复杂且相互关联，目前的研究仍偏向单一过程和案例的不确定性，缺乏推广能力[78]. 目前，稻田氮磷径流流失模型的不确定性研究包括模型和数据两个方面[79].

首先，模型不确定性包括模型结构和模型参数的不确定性. 对于稻田流失机理和过程认识的不完整，数学表达中的简化处理，模拟的尺度问题，全局最优的参数估值过程中伴随的参数敏感性、相关性、阈值等因素的影响，都会造成这种不确定性[20]. 在模型改进历程中，流失机理逐渐完善，考虑过程逐渐周全，模型也具备了模拟更大时空尺度的能力. 如Gao等[40-41]在物理模型中考虑降雨引起的释放作用后，研究者[42-43,80-81]又陆续考虑了渗透、化学转化等过程和不同土壤、水力条件下的适用性，降低了模型结构的不确定性. “试错法”是率定稻田氮磷径流流失模型参数最传统的方法，该方法费时费力，且带有很大的主观性和偶然性[47]. 随着计算机技术的发展，贝叶斯法[82]、SCE-UA[83]、SCPSO[84]等自动优选方法逐渐得到开发和应用，对于模型参数的不确定性研究具有重要意义.

其次，数据不确定性包括输入数据和验证数据(即输出数据)的不确定性. 稻田径流流失模型的输入数据常包括基础地理、降雨灌溉、土地利用和土壤属性、农艺管理措施等，而模型还需要径流通量、氮磷浓度等观测数据进行校准和验证. 受制于数据的时空尺度和分辨率，模拟结果存在不确定性. 对于基础地理、降雨等数据，目前常采用改变其输入的空间分布和尺度的方法，考察其对于最终结果的影响[20]，而氮磷浓度、径流等观测数据的测量误差受实验方法和监测手段的制约，仍没有较好的不确定性分析方法. 有研究者尝试采用随机变量[85]或加入随机干扰数据[86]来量化数据的不确定性.

3 模型关键问题

3.1 发生机制的参数化

上述4类模型的稻田氮磷径流流失发生机制还需要更好的参数化. 目前对于溶质从土壤迁移到地表径流的机理及这些过程对于自然人文相关因素的定量响应的认识仍然有限. 比如不同土壤类型、溶质、土壤渗透率、地表覆被等因素下降雨径流混合层的深度都会存在差异[87-88]，目前仍然没有有效的方法进行精确定量[89]. 土壤-水释放迁移速率、传质系数等的区域适用性仍需要研究证实. 当未来气候超过当前波动范围时，基于现在实验条件下得到的稻田氮磷径流流失响应规律是否有效也是关键问题[90].

除此之外，大量的物理机理模型往往针对坡地和旱地，导致对于具有特殊性的稻田(较高的田面水位、淹水环境下特殊的土壤分层等)的模拟效果较差，而水文过程和生态系统模型对部分复杂过程难以进行微观定量模拟，这意味着稻田氮磷流失的发生机制以及参数化的研究仍然存在诸多空白. 虽然少量研究开始关注降雨引起水土界面氮磷释放的关键作用，但释放的间歇性、孔隙水浓度、田面水中推流输移时再吸附和再沉淀等过程仍然缺乏有效参数化，无法区分常规降水和极端降水变化(频率、强度、持续时间等增加[91])的影响.

稻田氮磷径流流失是氮磷循环过程的之一，并非孤立，而与氮磷利用、氮磷损失、水-土-气界面过程等密切相关. 例如，作物生理生态过程影响蒸散发和氮磷吸收利用，从而改变田面水位和田面水氮磷浓度等径流流失的初始条件[7]. 然而，经验、物理机理和水文过程模型缺乏作物生理生态过程模拟能力，而农田生态系统模型缺乏稻田氮磷径流流失对作物生理生态的响应模拟能力. 因此，不同类型模型的耦合可能是解决这类问题的方法之一.

3.2 模型参数的合理化

模型率定和验证是保证参数合理的关键步骤[92]. 参数的选择往往被限制在预定义的范围内，再根据观测数据进行最优值的选定，以反映参数的区域适应性[93]. 传统的生态系统模型，通常在站点尺度上校准，升尺度到大尺度[94]，因此，需要高分辨率、多年多站点观测数据进行模型率定和验证，降低模型预测的不确定性. 对于水文过程模型，需要选择具有代表性的观测站点数据来反映空间分异性，但常规的观测数据难以满足要求，而人工采样难以提高时间和空间分辨率. 此外，目前率定数据多采用实际环境状态变量的观测值，即“硬约束”观测数据. 最近研究学者提出“软约束”观测数据[95-96]，即与模拟的各独立过程和平衡有关的信息[97]，这将是提高模型参数的合理化的可行途径.

稻田氮磷径流流失模型普遍缺乏针对每个发生过程(释放、径流、输移)的逐一率定，且氮磷循环过程的通量也难以同步率定，这就造成为了模拟氮磷径流流失过程时牺牲其他过程的模拟精度，从而降低模型整体的可信度. 对于高度非线性模型，很可能只在参数空间搜寻到一组最佳组合[69]. 目前，对于流失模型多采用单目标函数进行参数率定，基于有限特征指标的多目标全局优化方法尚未能考虑目标组合优化和多目标函数对模型参数集的影响，更难以解决异参同效问题[93].

目前，参数自动优选方法得到了普遍应用[82-83]，但基于启发式搜索的模型参数自动优选算法仍存在诸多问题，搜索算法的盲目随机游动性问题突出，搜索效率较低，容易陷入局部陷阱[98]. 王旭辉[69]利用粒子滤波算法构建了ORCHIDEE-crop模型的参数最优化系统，使得模型具有鲁棒性，但由于观测数据误差获取的限制，仍无法很好做到上述的多变量同步优化，且大规模计算能力也存在技术瓶颈. 目前针对稻田的模型多是采用成熟的参数率定方法，考虑到稻田的特殊性，未来仍需要采用更有效的方法来分解模型参数、模型结构、输入数据集、误差结构和概率对模型模拟结果的不确定性. 从田块尺度到流域、区域、全球尺度，主控因素也会发生变化，参数值可能不适用，需要进行建立传递函数实现参数升尺度[21].

3.3 农艺管理与措施的模块化

农艺管理措施显著改变稻田氮磷径流流失过程[99]，但相关模块开发和改进仍然存在极大的空间. 稻田的灌溉、排水方式与旱地显著不同，并直接影响着稻田氮磷径流流失. Yang等[100]发现优化灌溉相比于传统的大水漫灌能够显著减少稻田24%的氮流失，而移栽前田间被动排水也是一大贡献因素[101]. 然而，目前的模型缺乏自动灌排制度的模拟能力. 最近研制了针对稻田的灌溉排水控制方案模拟[101-102]，但过于简单，无法涵盖灌排制度的多样性和区域性. 这类模拟还受到实验观测的限制，未能进行同步率定和长期的验证，仍无法直接外推到实验区以外的稻田. 此外，目前的模型中还缺乏节水灌溉和田-沟-塘的联合调控的措施模拟. 例如，排水沟渠和塘堰系统是稻田与水体之间的一个过渡带，能够在减少稻田进入下游水体的氮磷流失方面起到重要作用[103]. 稻区沟、塘数量多且形状极其不规则，田-沟-塘的联合运行多种多样，如何识别其关键参数及其定量化是解决不同措施下稻田氮磷径流流失模拟问题的核心[102].

稻田的总氮流失量还受种植制度、施肥制度(量、类型、时间和方式)的影响[104]. 例如，深施相比表施能显著减少稻田氮径流流失的48%～72%. 不管是统计经验模型还是过程模型[5,25,69]，模型缺乏稻田氮磷径流流失对移栽和收获日期、施肥制度(量、类型、时间和方式)、秸秆还田比例等管理措施的响应能力. 另外，对于站点或小尺度而言，农艺管理措施可以通过调查直接输入，但对于大尺度而言则难以获取到调查数据. 如何建立种植制度(轮作与否、移栽和收获时间)、施肥管理(量、类型、时间、位置)，耕作(耕作、少耕、免耕)等大尺度高分辨数据集也是农艺管理与措施模块化的另一个关键.

4 展望

自1950s以来，稻田氮磷径流流失模型研究已经取得了显著的成果，在全球范围得以广泛应用，为稻田流失机理解析、通量估算及减排措施做出了重大贡献. 然而，稻田氮磷径流流失模型在发生机制、模型参数和农艺管理措施等关键问题方面尚存在不足，仍有较多的改进空间，本文总结前述，展望未来改进与发展的趋势如下：

1)提高观测数据质量和过程机理研究. 稻田氮磷径流流失过程机理复杂，模拟参数众多，不确定性大. 未来将需要改进实验设计和监测技术，获得更高时间分辨率、全通量、多站点、多种自然人文环境驱动、以及不同水稻品种和农艺管理措施下的实验观测数据，并通过机制平台共享数据，以达到模型应用的标准化和模型不确定性的降低. 另外，增强人工模拟试验与模型模拟的结合，进一步加强稻田氮磷流失形成过程各个环节的机理研究，对稻田释放、吸附解析过程和混合层深度等进行研究量化. 关注农艺管理措施的发展和改进[13,69]，实现农艺措施的模块化和参数化.

2)加强模型的开发改进与先进技术的集成. 计算机技术的发展和3S的应用给稻田氮磷流失模型的开发思路和手段带来了技术创新和革命. 近年来，人工智能和机器学习取得了飞速的发展，目前已经存在有机结合过程模型和智能模型的相关研究[105]，但尚未得到广泛的应用. 基于CPU和GPU异构并行的高性能计算技术也取得了显著的发展[106]，能够显著提高模型计算性能，增强模型可视化，推动稻田氮磷流失模型与先进技术手段的集成和融合.

3)加强氮磷流失不同过程、不同模块、不同模型类型之间的耦合. 稻田的氮磷流失是不同因素综合作用的结果，包含了复杂的物理化学生物过程的耦合和互馈，一些过程是相同的，在模型中互相借鉴，如植物生长过程；一些过程是相关的，如水文过程与生态过程，因此，不同学科的模型耦合成为必然. 目前，生态模型对于稻田氮磷迁移转化流失模拟仍较为简单，而水文模型用于稻田模拟时对于作物生长和管理措施考虑还不够全面. 最近就有研究者将生态模型APEX和水文模型SWAT模型结合起来，开发了混合模型SWAPX，来探索气候变化情景下稻区最优管理措施[107]. 因此，耦合稻田生态系统模块和流域水文模型，实现水文和稻田氮磷循环系统模拟，建立多模型集合框架[108]，有利于提高多种胁迫压力(包括极端气候事件)下模型模拟精度[109].

5 附录

附表Ⅰ和Ⅱ见电子版(DOI：10.18307/2021.0202).