利用地下水-农艺耦合模型评价涝渍地区作物-土壤-水动态

摘要

长期以来，农田内涝一直是一个众所周知的问题，导致了不利的社会、物理、经济和环境问题。为了更好地解决这一问题，需要更好地了解复杂的植物-土壤-水动力学系统。挑战在于模拟系统中组件之间的交互。目前存在模拟植物-土壤-水系统的模型，但这些模型要么独立运行，要么使用集成模型具有很高的侵入性。本文提出了一个紧密耦合的模型，DayCent-MODFLOW，它将三维地下水流(MODFLOW)模型和一维农业生态系统模型(DayCent)联系起来。DayCent负责根区植物-土壤-水动态，而MODFLOW模拟含水层饱和区水头和地下水流动。DayCent从土壤剖面向地下水位深层渗透，在地下水位位于土壤剖面内的内涝条件下，DayCent土壤水文过程受到MODFLOW模拟地下水位存在的约束。采用并行进程间通信技术MPI (Message Passing Interface)实现了耦合。该模型应用于美国科罗拉多州北部一个涝渍农业区(22平方公里)，并对地下水水头和蒸散速率(ET)进行了测试。该模型在最大的AWS Linux服务器上与多个进程并行运行。地下水水头与实测水头匹配程度较好，蒸散发率与参考蒸散发率匹配，且与作物类型高度相关。结果表明，两个模型之间存在较强的水文相互作用。模型还估算了涝渍条件下土壤温室气体排放(N2O和CH4)。尽管该模型可用于模拟任何植物-土壤-含水层系统，但无论地下水位深度如何，本研究的结果表明，该模型可用于评估浅层地下水地区的作物生产力、补给、ET和温室气体排放。

1. 介绍

植物-土壤-水系统控制着农业景观中水分、养分和温室气体的流动。在各种水文条件下了解这一系统对粮食生产、土地管理、水管理和营养管理都很重要。植物-土壤-水系统是一个复杂的系统，包括地表水径流、入渗、土壤水分动力学、蒸散发、补给和养分淋溶、碳(C)-氮(N)循环以及随后的水化学过程，如含水层的流动和养分运输、水流排放和养分负荷以及温室气体排放。

在这样一个复杂的系统中模拟水运和养分循环的挑战主要在于“带”之间的相互作用，例如土壤剖面和含水层饱和带之间的水运动。如Alley (Alley et al .， 2002)所述，由于降水、灌溉、蒸散、土地利用、作物类型和土壤类型等因素的影响，地下水补给是最难模拟的地下水预算。植物-土壤-水相互作用的一个特殊条件是作物根区存在饱和条件(即“涝渍”)，这会降低作物产量，损害土壤健康和结构(Cannell等，1980;Cavazza和Pisa, 1988;Collaku & Harrison, n.d;Houk等人，2006,2006,2006;Kaur et al, 2017)。2010-2016年期间，美国每年因内涝造成3.6亿美元的作物生产损失，其损失甚至超过干旱(Ploschuk等，2018)。涝渍可以极大地改变土壤中碳和氮的动态。由此产生的厌氧条件降低了有机物分解的速度(Meurant和Riker, 2014)，导致土壤有机质积累，影响氮矿化和作物吸收的有效氮;也增加了反硝化作用，从而增加了甲烷(CH4)排放和氧化亚氮(N2O)排放(Bartlett和harris, 1993;Parton et al, 2001)，两者都是温室气体。

有许多基于数值物理的模型可以模拟植物-水-土壤系统中的一系列水文和化学过程。这些模型的一个子集是在土壤剖面尺度上模拟一维域水文的农艺模型。其中包括SWAP (Kroes等人，2009)、DSSAT (Jones等人，2003)和DayCent (Parton等人，1998;Zhang et al .， 2018a,b)。它们模拟灌溉、径流、渗透和土层渗透、作物蒸散发和土壤剖面底部的深层渗透。但是，由于它们没有模拟下伏含水层饱和区的地下水流动，因此没有反映地下水位的波动及其在土壤剖面和作物根区可能存在的情况。水土评估工具(SWAT) (Arnold et al .， 1998)模拟了流域尺度上的水文过程和作物产量，在整个流域景观的单个水文响应单元(hru)上发生了水平衡。该模式考虑了地下水储存和地下水向溪流的排放，但没有以物理方式模拟地下水位波动，因此不能考虑内涝对根区过程和作物产量的影响。

即使是相关的SWAT-MODFLOW模型(Bailey et al, 2016)也没有考虑到根区浅层地下水的状况——MODFLOW可以模拟HRU土壤剖面高度的地下水位，但它对SWAT的HRU土壤剖面和根区过程没有影响。链接DSSAT- modflow模型(Xiang et al .， 2020)也没有考虑浅层地下水对根区过程的影响，因为模型之间的链接是在年时间尺度上进行的，并且每天的DSSAT根区水文过程不受modflow模拟的地下水位的影响。

另一部分模型模拟水汽带水文过程和地下水位波动，但不模拟近地表水文、营养生长、根区过程和养分循环。其中包括水文和水文地质模型MODFLOW (Niswonger等人，2011)、HYDRUS (Simůnek等人，2012)、MODFLOW- surfact (Panday和Huyakorn, 2008)、STOMP (White和Oostrom, 2003)、TOUGH2 (Pruess等人，1999)、VS2DI (Healy, 2008)、VSF软件包(Thoms等人，2006)和HydroGeoSphere (Therrien等人，2010)等。目前普遍缺乏水文农艺模型，其中模拟的地下水位影响根区过程，而根区过程反过来又影响地下水位的补给。

本文的目的是提出一个耦合的农业生态系统地下水模拟系统，可以模拟水浸条件下的水运动。分别选择DayCent和MODFLOW作为农业生态系统模型和地下水模型。DayCent之所以被选中，是因为它包含了农业生态系统中的碳和氮动态，因此在应用中更加通用。之所以选择MODFLOW，是因为它是世界上使用最广泛的地下水流动模型(Langevin et al, 2017)。这些模型是紧密耦合的，使用新颖的消息传递接口(MPI) (Gropp等人，1996)，在它们之间每天的时间步上传递系统数据，避免了对DayCent和MODFLOW代码的代码修改，因此可以与DayCent和MODFLOW的更新版本一起工作。DayCent改进了MODFLOW的补给计算，MODFLOW提高了DayCent在浅水区存在下作物产量和养分循环的准确性。尽管DayCent-MODFLOW连接系统可以应用于任何植物-土壤-水系统，但在美国科罗拉多州北部具有浅层地下水的灌溉区的应用将展示其在淹水条件下的能力。在模型应用中，还将简要描述涝渍条件对温室气体的影响。