应用生态学报 2005年11月 第16卷 第11期 CHINESEJOURNALOFAPPLIEDECOLOGY,Nov.2005,16(11)∀2208~2212
土壤氮素循环模型及其模拟研究进展
唐国勇
3
1,2
*
黄道友 童成立
11**
张文菊
1,3
吴金水
1
(1中国科学院亚热带农业生态研究所亚热带农业生态重点实验室,长沙410125;2中国科学院研究生院,北京100039;华中农业大学资源环境学院,武汉430070)
摘要 N既是植物必需的营养元素,又是造成环境污染的重要元素.正确模拟土壤中N循环已经成为科学家共同关注的热点问题.简述了土壤N循环的基本过程,重点介绍了13种土壤N循环模型和6个土壤N循环过程的模拟,并讨论了模拟中存在的参数化问题.
关键词 土壤N N循环 模型 模拟
文章编号 1001-9332(2005)11-2208-05 中图分类号 S153.6 文献标识码 A
Researchadvancesinsoilnitrogencyclingmodelsandtheirsimulation.TANGGuoyong1,2,HUANGDaoyou1,TONGChengli1,ZHANGWenju1,3,WUJinshui1(1KeyLaboratoryofSubtropicalAgro ecology,InstituteofSubtropicalAgriculture,ChineseAcademyofSciences,Changsha410125,China;2GraduateSchoolofChineseAcademyofSciences,Beijing100039,China;3CollegeofResourcesandEnvironment,HuazhongAgriculturalUniversity,Wuhan430070,China). Chin.J.Appl.Ecol.,2005,16(11):2208~2212.
Nitrogenisoneofthenecessarynutrientsforplant,andalsoaprimaryelementleadingtoenvironmentalpollu tion.Manyresearcheshavebeenconcernedaboutthecontributionofagriculturalactivitiestoenvironmentalpollu tionbynitrogenouscompounds,andthefocusishowtosimulatesoilnitrogencyclingprocessescorrectly.Inthispaper,theprimarysoilnitrogencyclingprocesseswerereviewedinbrief,with13cyclingmodelsand6simulatedcyclingprocessesintroduced,andtheparameterizationofmodelsdiscussed.Keywords Soilnitrogen,Nitrogencycle,Model,Simulation.
1 引 言
N是植物必需的营养元素,也是评价土壤质量和土地生产力的重要指标.为了获得高产,需要施用大量的氮肥.据统计[32],仅1996年全世界氮肥(折纯N)使用总量就高达8 50!107t,但N累积利用率不高.据估计,施入土壤中的N大约有35%通过各种途径损失掉[6,32].此外,氮肥的使用还可能造成环境污染,诸如温室气体(主要是氮氧化物)和致酸雨气体(氨气)的排放、地下水硝酸盐超标、水体富营养化等[20].如2000年,比利时80%的饮用水中硝酸盐含量超标[10].目前,土壤N循环的研究已经成为土壤学家、环境学家、农学家等共同关注的热点问题之一.
土壤N循环是N生物地球化学循环中的重要环节,其模拟是作物估产、环境评价、农田管理、决策制定和长期预测的重要依据,对提高氮肥利用率、防止或减轻环境污染具有重要的理论和实践意义.20世纪60年代,就有基于单个过程的土壤N循环方面的报道[25,28].40多年来,北美和欧洲一些国家建立了大量的土壤N循环模型.我国在这方面研究还比较薄弱[3,15,24].本文拟通过简要概述土壤N循环过程,重点介绍13种土壤N循环模型和6个土壤N循环过程的模拟,并讨论模型模拟中的参数化问题,以期为深入研究土壤N循环及其模拟提供一定的参考和借鉴.
2 土壤N循环的基本过程
土壤中含N化合物种类多,理化、生物学性质各异.一般可将土壤中N划分为有机氮和无机氮,以有机氮为主.在土壤微生物等因子的作用下,N在土壤中发生一系列复杂的循环.主要循环过程有:有机氮矿化、腐殖化、硝化、反硝化、氨挥发、N沉降、硝酸盐淋失、生物固氮、铵离子晶格固定和释放、土壤粘粒吸附和解吸、植物吸收等过程.土壤N循环过程的研究是建立土壤N循环模型以及N生物地球化学循环模型的基础.
3 土壤N循环模型的研究概况
目前,农业中数学模型并无统一的分类,可从不同角度进行划分.根据建模的方法可分为经验模型和机理模型;从土壤有机氮角度可分为单组分和多组分模型;从模拟循环过程的数目方面可分为单过程和多过程模型;此外,根据模型模拟的元素也可分为独立N模型和综合模型的N子模型.
经验模型通常依据实验测定或调查的N循环分量与气
*中国科学院知识创新工程重要方向项目(KZCX3 SW 426)、国家
自然科学基金重点项目(40235057)和国家重点基础研究发展资助项目(2002CB412503).**通讯联系人.
2005-01-10收稿,2005-05-08接受.
11期 唐国勇等:土壤氮素循环模型及其模拟研究进展 候因子等已知环境变量间的数据关系,采用统计方法,通过曲线拟合建立方程,以此进行外推,如Shaffer等[29]建立的NTRM.该模型主要用于评估土壤侵蚀对土地生产力、植物竞争、作物产量和水质的影响,较详细地模拟了水分和化学物质运移、作物生长、根系发展、田间管理效应、土壤C和N动态、径流、土壤化学平衡等方面.本模型受美国农业部大平原系统研究单位支持和推广应用.在模拟N动态方面,该模型只适用于小时空尺度内N循环过程的定量研究,主要包括对矿化、腐殖化、硝化、尿素水解、植物吸收等过程的模拟.模型运行所需的参数较少,主要参数有土壤温度、土壤氧气2209
(3)
dyt/dt=-umy/(Ks+Y)
其中,k0和k1分别为零级和一级动力学常数,um为最大转化速度,Ks为半饱和常数,dyt/dt为转化速率,Y为某N成分的含量(或浓度).
20世纪60年代后,国际上建立了大量的土壤N模型,特别是随着计算机的广泛应用,土壤N循环模型有了长足的发展.其中影响较大的有:美国的NTRM、NLEAP、NC SOIL、DNDC、EPIC、CERES、CENTURY、GLEAMS、RZWQM模型;英国的SUNDIAL模型;荷兰的SOILN模型;丹麦的DAISY模型以及德国的CANDY模型.这些模型含量、土壤水分含量、土壤pH值等.在NTRM模型的基础上,经过参数调整、修正和部分理论推导,发展成为NLEAP和RZWQM等机理模型,但土壤N硝化反应由最初的回归方程演变成机理的半阶方程[7,29].一般经验模型所需的参数较少,使用比较方便,应用较为广泛,但不能准确地反映土壤N的长期动态,有较强的地域性.机理模型则从N循环机理着手,研究各种环境因子对N循环过程的影响机理,如由英国洛桑试验站开发的SUNDIAL模型
[2]
.该模型主要用于评
估农业N循环对环境的影响.模型运行以周为步长,经过洛桑实验站长达150年定位试验的检验,证明该模型能较好地预测土壤中N长期动态,不适用于模拟土壤短期(1个月内)N动态.模型主要模拟了矿化、腐殖化、硝化、反硝化、氨挥发、植物吸收、硝酸盐淋失等过程.但氨挥发过程只模拟了肥料表施和衰老植物组分中氨的挥发,硝酸盐淋失也只模拟硝酸盐淋失至深层地下水这一过程.在有机N矿化过程中,该模型考虑了粘粒含量对矿化速率的影响.模型运行所需的参数主要包括周平均气温、田间最大持水量、最大有效持水量、土壤粘粒含量、基质组分中的C/N、各过程反应速率常数、作物目标产量等.通常机理模型理论性强,延展性好,可用于模拟土壤N长期动态,但模型运行所需参数较多,要求建模者对整个N循环过程都有深入的研究.早期土壤N模型以经验模型为主,近期则多为机理模型.
单过程模型只考虑N循环中某一过程.国内N循环模型主要为单过程模型.多过程模型研究两个或多个N循环过程以及这些过程的相互关系.N循环过程多在土壤微生物驱动下进行,而微生物活动又依赖于土壤有机碳(作为微生物能源物质)的供应.因此,大多数N模型以C/N联合模型而存在,甚至是生态系统物质循环综合模型的子模型,如CERESN子模型、CENTURYN子模型等.
在N循环过程模型中,通常以零级动力学方程、一级动力学方程或米氏方程等作为模型的数学基础,根据N循环过程的影响因子,对模型进行调整和修正.零级动力学方程的微分形式可表示为:
dyt/dt=-k0
(1)
一级动力学方程的微分形式可表示为:
dyt/dt=-k1y
(2)米氏方程的微分形式可表示为:
各有其长处和不足,已有文献对此进行了介绍和比较[9,21,24],本文仅对这13种N循环模型进行简要的介绍(表1).
4 土壤N循环过程模型模拟4 1 矿化和腐殖化过程
大多数土壤N循环模型用一级动力学方程模拟土壤有机N的净矿化或净腐殖化.在NCSOIL模型[23]中,土壤基质按C/N和土壤有机质半减期(half life)分为4个组分(残体库、微生物库、中间库和稳定腐殖质库),每个库又分为易分解和难分解两部分,每组分有机N的矿化用一级动力学方程模拟,其微分形式为:
dC/dt=min(fT,fW) N[SLkLC+(1-SL)hRC](4)其中,C是土壤基质碳含量,SL和(1-SL)分别为基质中易分解和难分解部分含量,kL和hR分别为相应的分解速率常数,fT和fW分别为温度和水分对矿化的影响因子, N为基于C/N的降低因子.土壤N净矿化或净腐殖化取决于基质残体库和微生物库C/N[23].各组分的分解速率常数、土壤温度、土壤最大持水量、基质C/N、基质中易分解与难分解比例、耕作措施等参数是模型运行所必需的,如耕作能提高土壤基质中间库易分解部分的比例[23].
影响矿化和腐殖化过程的因子很多,模型大多模拟了土壤温度、土壤含水量(或土水势)、土壤通气状况、土壤有机碳含量和土壤C/N对该过程的影响[1,2,7,9,13,14,26],也有模型考虑了土壤微生物量[22,26,33,35]
、碳磷比
[11,35]
、土壤容
重[35]、土壤质地[26]、土壤离子强度[26]、木质素含量[26]等可
能影响该过程的因子.4 2 硝化过程
硝化过程可用零级、一级或米氏方程进行模拟.在NTRM模型中,最初用经验的回归方程表示硝化速率[7]:
K=a+bTCNH4
+clog10CNH3
+dlog10CNO
3
(5)
其中,K为硝化反应速率,a、b、c、d均为常数,其值分别为4 64、1 62!10-3、0 238和-2 51,T为土壤温度,CNH4为土壤铵态氮含量,CNO3
为基质硝态氮含量.经过修正,Shaffer
等
[28]
用机理性的半阶方程来表示硝化速率:
d(NH4)/dt=[k1/2(NH4)1/2O2exp(Ea/kbTs)]/H+(6)
应 用 生 态 学 报 16卷2210
表1 13个土壤N循环模型简介
Table1Introductionto13soilnitrogencyclingmodels模 型ModelsNTRM
模型简要描述Descriptionofmodels
该模型为经验模型,主要用于评估土壤侵蚀对土壤生产力、植物竞争、作物产量和水质的影响.该模型没有模拟环境因子等对N循环过程的影响,只适用于小时空尺度内N循环过程的定量研究.该模型未将土壤有机N进行分组,用回归方法建立了矿化过程、腐殖化过程、硝化过程、尿素水解、硝酸盐淋失过程的方程;植物吸收N与吸水量成线性关系[7,29]
该模型可模拟硝酸盐的淋失,用于评估水质,能与GIS结合.该模型以天为步长,适用于有机N中长期模拟预测.模型将土壤有机N分为3个组分库,其中矿化、腐殖化、反硝化、氨挥发、硝酸盐淋失过程用一级动力学方程模拟;用logistic曲线模拟作物吸N量;用零级动力学方程模拟硝化过程[30]
该模型主要用于评估土壤N循环对环境的影响.模型运行以周为步长,适用于土壤中长期N的预测.将土壤有机N分为3个组分库,用一级动力学方程模拟矿化过程、腐殖化过程、硝化过程;用土壤有机质分解产生CO2量的线性方程表示反硝化过程;氨气挥发量与铵态肥施入量呈线性关系;植物吸收过程用植物目标含N量表示[2]
该模型主要用于评估施用N对环境的影响.早期版本中将土壤有机氮分为3个组分库.用一级动力学方程模拟矿化、腐殖化、硝化过程;用零级动力学方程模拟反硝化过程;用logistic曲线模拟植物吸N量,并模拟氨挥发过程和硝酸盐淋失过程[17]
该模型主要用于评估农家肥、泥浆等含N化合物的施用对环境的影响[14,36].该模型以日为步长,适用于短期土壤N动态的模拟.早期版本中并未对土壤有机氮进行分组[14],新版本中将有机N分为6个库[36].用一级动力学方程模拟矿化和腐殖化过程;用米氏方程模拟硝化过程;植物吸N量取决于植物潜在需求和土壤供应量;也简单地模拟了硝酸盐淋失过程[14,16]
该模型主要用于6类含N有机化合物的施用对环境的影响.模型中将土壤有机氮分为3个组分库,矿化过程、腐殖化过程用一级动力学方程模拟;反硝化过程与硝酸根浓度和有机碳量成比例;也简单地模拟了氨气挥发情况;硝化过程用米氏方程模拟[8]该模型用于模拟土壤C、N动态和痕量气体的排放.该模型将土壤有机氮分为3个组分库,跟踪有机碳、氮的逐级降解;由降解作用产生的可给态C、N,被输入硝化、脱氮及发酵子模型中,进而模拟有关微生物的活动及其代谢产物,包括几种温室气体的排放[19]该模型用于模拟土壤C、N循环,与其它模型不同的是在研究N循环机制中,该模型能同时模拟标记N和非标记N的去向.模型将土壤有机氮分成4个组分库.用米氏方程模拟矿化和腐殖化过程;用零级动力学曲线模拟硝化过程、反硝化过程;用logistic曲线模拟作物生长[23]
该模型由一系列作物生长子模型构成,各子模型中的概念、模拟方法和影响因子相同.主要用来模拟作物的生长及其与之有关的过程.模型中将土壤有机氮分为腐殖质库和新鲜有机残体库.矿化和腐殖化、反硝化、氨挥发过程用一级动力学方程模拟;硝化过程用米氏方程模拟;模型详细地模拟作物对N的吸收过程[13]
该模型主要用来模拟生态系统C、N和营养物质长期动态[26,27].早期版本中以月为步长[26],新版本中某些过程以
[27]
日为步长.由于模型运行的步长较长,模型中N子模型部分较简单.矿化和腐殖化过程用一级动力学方程模拟;植物吸N量取决于植物N需求和土壤N供应;用总矿化N的5%表示氨挥发量[26];在最新版的DAYCENT模型中,硝化过程和反硝化过程用一级动力学方程模拟[27]
该模型主要特色就是能模拟不同农业管理措施及其对水质和作物产量的影响.模型将土壤有机氮分为5个组分库.用包含微生物生长的一级动力学方程模拟矿化和腐殖化过程;用一级动力学曲线模拟反硝化过程、氨挥发过程;用米氏方程模拟作物吸N量;用零级动力学过程模拟硝化过程;模型详细地模拟了硝酸盐淋失过程[1]
该模型主要用于模拟农业管理对地下水负载的影响.模型将土壤有机氮分为3个组分库.矿化和腐殖化过程、反硝化过程、氨挥发过程用一级动力学方程模拟;硝化过程用零级动力学方程模拟;简单地模拟了土壤侵蚀和地表径流中N的损失[18]
该模型主要用于评估土壤侵蚀、耕作方式对生产力的影响.土壤有机氮的分组与GLEAMS模型中的分组一样,只是其中的活性库是估算的,而非直接测定得出.矿化过程、腐殖化过程、硝化过程、反硝化过程、氨挥发过程用一级动力学方程模拟;植物吸N量取决于植物N需求和土壤N供应;也模拟了土壤侵蚀和地表径流中N的损失[35]
NLEAP
SUNDIAL
SOILNDAISY
CANDY
DNDC
NCSOIL
CERES
CENTURY
RZWQM
GLEAMS
EPIC
其中,d(NH4)/dt为硝化速率,k1/2为半阶速率常数,Ea为化学自由能,kb为波尔茨曼常数,Ts为土壤温度,NH4为土壤溶液中铵离子浓度,O2为土壤空气中氧气浓度,H+为基质氢离子浓度.在NTRM模型中,无论使用经验的回归方程还是机理的半阶方程,土壤有机氮均不分组.回归方程中所需的参数只有土壤温度、土壤硝态氮含量和铵态氮含量,而半阶方程中所需的参数包括半阶反应速率常数、土壤铵态氮含量、土壤空气中氧气含量、土壤pH值、土壤温度.
影响硝化过程的因素很多,土壤温度、土壤含水量、通气状况、土壤pH[1,2,7,9,13,14,26]是大多数模型所考虑的主要环境因子.微生物种群动态及活性[1]、硝酸根与铵根浓度比[17]等因素对硝化过程的影响也不容忽视.一些模型考虑了硝化过程中N2O的释放情况,通常用被硝化总N的一定比例来表示N2O释放量
[7,11,26]
拟
[1,2,7,9,17,26]
.在通气良好的土壤、土表薄层水体、有氧的根
际中硝化作用都比较明显[5].4 3 反硝化过程
目前反硝化过程的模拟一般采用零级或一级动力学方程,如CERES模型用一级动力学方程模拟该过程,其一级反硝化速率系数为:
kdnit=6 0!105!CW!Tf!Wfd!h!NO3
(7)
其中,kdnit为一级反硝化速率系数,CW为土壤有机质中水提取C浓度( g#g-1),Wfd和Tf分别为水分和温度影响因子,h为土层厚度.CW、Wfd和Tf又可表示为:
CW=24 5+0 0031CTf=0 1!exp(0 046Ts)
Wfd=1 0-(SAT-SW)/(SAT-DUL)
(8)(9)(10)
.也有模型用氨的相对权重而不是
[35]
其中,C为土壤有机碳含量,SAT为土壤饱和水含量,SW为土壤含水量,DUL为土壤排干时水分含量,Ts为土壤温度.该模型认为,只有当土壤含水量超过土壤排干时水分含量的情况下,反硝化反应才发生.在目前众多土壤N循环模型中,
铵离子有效浓度来表示硝化速率.此外,硝化抑制剂的使
用能显著降低硝化反应速率[18].土壤通气状况是影响硝化过程的决定性因素,常用氧分压的非线性方程来模
11期 唐国勇等:土壤氮素循环模型及其模拟研究进展 该模型是为数不多的能用于模拟稻田土壤反硝化作用的模型之一.
多数模型模拟了土壤温度、土壤含水量、厌氧状况、土壤有机质含量(尤其是土壤水溶性C含量)对反硝化过程的影响.有模型也考虑了不同反硝化产物的情况,主要是区分了N2和N2O的不同排放量以及环境因子对N2和N2O排放的影响[2,12,26,29,31].4 4 氨挥发过程
可用零级或一级动力学方程模拟氨挥发过程,如GLEAMS模型用一级动力学方程表示该过程:
2211
定量研究硝酸盐的淋失对环境的影响是许多土壤N循环模型建立的主要目的之一,如DAISY、GLEAMS、SOILN、EPIC、NLEAP、SUNDIAL、CANDY等模型.硝酸盐淋失过程属于土壤溶质迁移的范畴.在土壤N循环模型中,常用土壤溶质迁移方面的模型或方程来模拟硝酸盐淋失.土壤N循环模型主要应用对流 扩散模型[12]、两区模型[7,26,29]、两流区模型[2]来模拟该过程.硝酸盐淋失与土壤含水量和土壤水分运动关系密切.土壤水分运动包括土壤水分入渗和再分布.土壤水分入渗常用Green Ampt[1,13,26,29]来模拟,而再分布过程用Richards方程
[29]
来拟合.国内外已有大量文献对硝酸盐
r=CNH3
exp(-Kvt)
其中,r为一级氨挥发速率,CNH3为基质中氨浓度,t为时间,Kv为速率常数,该常数又可用指数方程表示:Kv=0 409!1 08ATP-20,ATP为气温.但该方程只适用于模拟表施铵态氮肥和厩肥过程中氨气挥发情况.该模型对氨挥发过程的模拟步长为1d,模拟时空尺度较小.
影响氨挥发的因素主要是土壤温度、土壤含水量(或土水势)、氨气压梯度、风速.大多数模型只模拟了这几个主要因素[1,2,13,26,29,35].也有模型考虑了土壤阳离子交换量[2]、土壤表面覆盖
[13]
、肥料(铵态氮肥和尿素)的施用方式
[13]
、土层
深度[35]等.也有模型用气温代替土温来模拟温度对氨挥发的影响[18].大多数模型认为氨挥发只发生在土表,也有模型认为整个土体中都存在氨挥发[2].
与氨挥发相反,N沉降一般不被建模者所重视.事实上,在正常大气氨浓度下,N沉降中能被植物有效利用的N占植物N需求量的10%左右[16].而在局部高浓度氨的大气中,该机制所占的份额更大,如氮肥厂、化工厂等氨高排放区域[34].
4 5 NH4+的晶格固定和释放过程
土壤NH4+的晶格固定和释放是土壤N有效性研究中一个备受关注的问题.研究者
[22,37]
认为,
NH+4
的晶格固定
和释放类似于土壤吸附和解吸过程,可用等温吸附方程来描述,如Freundlich、Langmuir方程.影响该过程的因素主要包括土壤有机质含量、土壤溶液中
NH+4
浓度、土壤溶液中钾离
子浓度、土壤含水量、土壤矿物种类、土壤粘粒含量等.但这些因素对铵晶格固定和释放过程的影响难以用单值函数描述.如钾离子在低浓度时能促进NH+4
的固定,而高浓度时又
可能会抑制
NH+4
的固定
[22,37]
.土壤中蛭石、伊利石、蒙脱石
等2∀1型层状硅酸盐矿物含量高时,能加强
NH+4
的晶格固
定和释放,也能增加土壤固定NH+4量[4,22,37]
.NH+4的释放
除了与上述因素有关外,最重要的影响因素是NH+4固定的时间.一般可将固定态NH+4分为∃新固定的%和∃旧固定的%.∃新固定的%NH+4较易重新释放,而随固定时间的推移,固定
态NH+4就难以再释放出来
[4,22,37].目前,对土壤NH4+晶格固定和释放的研究还不够深入,还很难定量模拟环境因子对该过程的影响.4 6 硝酸盐淋失过程
淋失过程进行了报道和介绍[3,6,10,17,30,31],故本文不再论述.5 问题与展望
目前只有较少的土壤N循环模型在田间条件下得到比较充分的验证.其原因之一是模型参数化过程中参数的不确定性,而且这些不确定性对模型影响程度的评估也存在一定的困难.这种不确定性包括3个方面:1)模型参数求算过程的不确定性,尤其是不可直接测量的内部参数,如分解速率常数.这些参数只能通过其他数据资料进行计算.在计算过程中大多需要设定初始值或初始状态,而初始值或初始状态的设定通常是不确定的.2)测定资料的不确定性.土壤、大气中各成分的时空变异性很大,而由于测量技术和实验仪器的限制,测定的数据并不能完全反映这种变异性.由于测定资料的不准确,通过逐级放大作用,最终影响参数的准确性.3)过程间相互作用的不确定性.影响土壤N动态的因素很多,模型并不可能考虑所有的影响因素,有时只能用∃黑箱%表示,而且N过程间也存在着相互作用.此外,土壤N循环模拟是基于单个过程的,缺乏系统性.因此,对用户来说,好的模型最重要的要素就是:模型能一步一步地指导用户如何获取模型参数.
近年来,土壤N循环模型出现了两个大的发展趋势,即理论模型趋于复杂化、细微化,应用模型趋于简单化,易于用户使用.
从理论研究方面来看,土壤N循环模型将更加突出系统性,强调各过程间的相互联系.从应用方面来看,∃3S%(地理信息系统GIS,全球定位系统GPS,遥感RS)技术与土壤N循环模型相结合将是土壤N循环模型发展的趋势.目前已有较多的模型和GIS系统相结合,如CERES GIS,NLEAP GIS,RZWQM GIS,且CERES和RZWQM已形成专家系统.参考文献
1 AhujaLR,RojasKW,HansonJD,etal.2000.RootZoneWater
QualityModel:ModellingManagementEffectsonWaterQualityandCropProduction.LLC,HighlandsRanch,CO:WaterResourcesPublication.372
2 BardburyNJ,WhitmoreAP,HartPBS,etal.1993.Modellingthefateofnitrogenincropandsoilintheyearsfollowingapplicationof15
N labelledfertilizertowinterwheat.JAgricSci,121:363~3793 CuiJ B(崔剑波),ZhuangJ P(庄季屏).1997.SoilNO3- Ntrans portunderfieldunsaturatedflowcondition&Asimulationstudy.ChinJApplEcol(应用生态学报),8(1):49~54(inChinese)
应 用 生 态 学 报 16卷2212
4 DavidsonJB,AldagRW.1982.Inorganicformsofnitrogeninsoils.In:StevensonFJ,eds.NitrogeninAgriculturalSoils.Madison,WI:AmericanSocietyofAgronomy,Inc.48
5 DeBuskWF,WhiteJR,ReddyKR.2001.Carbonandnitrogendy
namicsinwetlandsoils.In:ShafferMJ,LiwangM,HansenS,eds.ModellingCarbonandNitrogenDynamicsforSoilManagement.Bo caRaton,FL:LewisPublishers.27~54
6 DelgadoJA,FollettRJ,ShafferMJ.2000.SimulationofNO3 Ndy
namicsforcroppingsystemswithdifferentrootingdepths.SoilSciSocAmJ,64:1050~1054
7 DuttRG,ShafferMJ,MooreWJ.1972.Computersimulationmodel
ofdynamicbio physicochemicalprocessesinsoils.ArizonaAgricExpStaTechnolBull,196:128
8 FrankoU,OelschagelB,SchenkS.1995.Simulationoftempera
ture,water,andnitrogendynamicsusingthemodelCANDY.EcolMod,81:213~322
9 FrisselMJ,VanVeenJA.1978.Acritiqueofcomputersimulation
modellingfornitrogeninirrigatedcroplands.In:NielsenDR,Mac DonaldJG,eds.NitrogenintheEnvironment.NewYork:AcademicPress.145~162
10 GarcetJDP,AmauryT,JavauxM,etal.2001.Modellingnitrogen
behaviourinthesoilandvadoseenvironmentsupportingfertilizermanagementatthefarmscale.In:ShafferMJ,LiwangM,HansenS,eds.ModellingCarbonandNitrogenDynamicsforSoilManage ment.BocaRaton,FL:LewisPublishers.571~596
11 GijsmanAJ,ObersonA,TiessenH,etal.1996.Limitedapplicabili
tyoftheCENTURYmodeltohighlyweatheredtropicalsoils.A gron,88:894~903
12 GrantRF,PatteyE.1999.Mathematicalmodellingofnitrousoxide
emissionfromanagriculturalfieldduringspringthaw.GlobalBio geochemCyc,13:679~694
13 HanksJ,RitchieJT.1991.ModellingPlantandSoilSystem.Madi
son,WI:AmericanSocietyofAgronomy,Inc.
14 HansenS,JensenHE,NielsenNE,etal.1991.Simulationofnitro
gendynamicsandbiomassproductioninwinterwheatusingtheDanishsimulationmodelDAISY.FertRes,27:245~259
15 HuangY(黄 耀).2003.CarbonandNitrogenExchangeinEarth
atmosphericSystem&FromExperimentstoModels.Beijing:Meteo rologyPress.(inChinese)
16 HutchinsonGL,MillingtonRJ,PetersDB.1972.Atmosphericam
monia:Absorptionbyplantleaves.Science,175:771~772
17 JohnsonHL,BergstromPE,Jansson,etal.1987.Simulatednitro
gendynamicsandlossesinalayeredagriculturalsoil.AgricEcosystEnviron,18:333~356
18 LeonardRA,KinselWG,StillDS.1987.GLEAMS:Groundwater
loadingeffectsofagriculturalmanagementsystems.TransASAE,30:1403~141819 LiC S(李长生).2001.Biogeochemicalconceptsandmethodolo
gies:DevelopmentoftheDNDCmodel.QuatSci(第四纪研究),21(3):89~99(inChinese)
20 LiwangM,ShafferMJ.2001.Areviewofcarbonandnitrogenpro
cessesinnineU.S.soilnitrogendynamicsmodels.In:ShafferMJ,LiwangM,HansenS,eds.ModellingCarbonandNitrogenDynam icsforSoilManagement.BocaRaton,FL:LewisPublishers.55~102
21 McGechanMB,WuL.2001.Areviewofcarbonandnitrogenin
Europeansoilnitrogendynamicsmodels.In:ShafferMJ,LiwangM,HansenS,eds.ModellingCarbonandNitrogenDynamicsforSoil
Management.BocaRaton,FL:LewisPublishers.103~171
22
McGillWB,HuntHW,WoodmanseeRG,etal.1981.Phoenix,amodelofthedynamicsofcarbonandnitrogeningrasslandsoils.EcolBull,33:49~115
23
MolinaJAE,ClappCE,ShafferMJ,etal.1983.NCSOIL,amodelofnitrogenandcarbontransformationsinsoil:Description,calibra tion,andbehaviour.SoilSciSocAmJ,47:85~9124
MuX M(穆兴民),FanS L(樊小林).1999.AreviewonecologymodelsofsoilNmineralization.ChinJApplEcol(应用生态学报),10(1):114~118(inChinese)
25 OlsonJS.1963.Energystorageandthebalanceofproducersandde composersinecologicalsystem.Ecology,44:322~331
26
PartonWJ,MosierAR,OjimaDS,etal.1996.GeneralizedmodelforN2andN2Oproductionfromnitrificationanddenitrification.GlobalBiogeochemCyc,10:401~412
27
PartonWJ,OjmaDS,ColCV,etal.1994.AGeneralmodelforsoilorganicmatterdynamics:Sensitivitytolitterchemistry,texture,andmanagement.SoilSciSocAmSpecPull,39:147~167
28
SabeyBR,FrederikLR,BarthdomewWV.1969.TheformationofnitratefromNH+4 Hinsoils∋.Useofthedelayandmaximumratephaseformakingquantitativepredictions.SoilSciSocAmJ,33:276~278
29
ShafferMJ,LarsonWE.1987.NTRM,Asoil cropsimulationmod elfornitrogen,tillage,andcrop residuemanagement.USDAConsResRep,34(1):103
30
ShafferMJ,HalvorsonAD,PierceFJ.1991.Nitrateleachingande conomicanalysispackage(NLEAP):Modeldescriptionandapplica tion.In:FollettRF,ed.ManagingNitrogenforGroundwaterQuali tyandFarmProfitability.Madison,WI:SoilScienceSocietyofAmerica,Inc.285~322
31
ShafferMJ,LasnikK,OuX,etal.2001.NLEAPInternettoolsforestimatingNO3 NleachingandN2Oemission.In:ShafferMJ,Li wangM,HansenS,eds.ModellingCarbonandNitrogenDynamicsforSoilManagement.BocaRaton,FL:LewisPublishers.403~42632 SmilV.1997.Nitrogenincropproduction:Anaccountofglobalflows.GlobalBiogeochemCyc,13:647~662
33 SmithOL.1982.SoilMicrobiology:AModelofDecompositionandNutrientCycling.BocaRaton,FL:CRCPress.273
34
WhiteheadDC,LockyerDR.1987.Theinfluenceoftheconcentra tionofgaseousammoniaonitsuptakebytheleavesofItalianRye grass,withandwithoutanadequatesupplyofnitrogentotheroots.ExpBot,38:818~827
35
WilliamsJR.1995.TheEPICmodel.In:SinghVP,ed.ComputerModelsofWatershedHydrology.HighlandsRanch,CO:WaterRe sourcesPublications.909~100036
XuC,ShafferMJ,Al KaisiM.1998.Simulatingtheimpactofman agementpracticesonnitrousoxideemissions.SoilSciSocAmJ,62:736~724
37
ZhuZ L(朱兆良),WenQ X(文启孝),ChengL L(程励励),etal.1992.NitrogeninChineseSoils.Nanjing:JiangsuScienceandTechnologyPress.(inChinese)
作者简介 唐国勇,男,1980年生,硕士研究生.主要从事土壤养分循环及模拟研究,发表论文2篇.E mail:[email protected].
应用生态学报 2005年11月 第16卷 第11期 CHINESEJOURNALOFAPPLIEDECOLOGY,Nov.2005,16(11)∀2208~2212
土壤氮素循环模型及其模拟研究进展
唐国勇
3
1,2
*
黄道友 童成立
11**
张文菊
1,3
吴金水
1
(1中国科学院亚热带农业生态研究所亚热带农业生态重点实验室,长沙410125;2中国科学院研究生院,北京100039;华中农业大学资源环境学院,武汉430070)
摘要 N既是植物必需的营养元素,又是造成环境污染的重要元素.正确模拟土壤中N循环已经成为科学家共同关注的热点问题.简述了土壤N循环的基本过程,重点介绍了13种土壤N循环模型和6个土壤N循环过程的模拟,并讨论了模拟中存在的参数化问题.
关键词 土壤N N循环 模型 模拟
文章编号 1001-9332(2005)11-2208-05 中图分类号 S153.6 文献标识码 A
Researchadvancesinsoilnitrogencyclingmodelsandtheirsimulation.TANGGuoyong1,2,HUANGDaoyou1,TONGChengli1,ZHANGWenju1,3,WUJinshui1(1KeyLaboratoryofSubtropicalAgro ecology,InstituteofSubtropicalAgriculture,ChineseAcademyofSciences,Changsha410125,China;2GraduateSchoolofChineseAcademyofSciences,Beijing100039,China;3CollegeofResourcesandEnvironment,HuazhongAgriculturalUniversity,Wuhan430070,China). Chin.J.Appl.Ecol.,2005,16(11):2208~2212.
Nitrogenisoneofthenecessarynutrientsforplant,andalsoaprimaryelementleadingtoenvironmentalpollu tion.Manyresearcheshavebeenconcernedaboutthecontributionofagriculturalactivitiestoenvironmentalpollu tionbynitrogenouscompounds,andthefocusishowtosimulatesoilnitrogencyclingprocessescorrectly.Inthispaper,theprimarysoilnitrogencyclingprocesseswerereviewedinbrief,with13cyclingmodelsand6simulatedcyclingprocessesintroduced,andtheparameterizationofmodelsdiscussed.Keywords Soilnitrogen,Nitrogencycle,Model,Simulation.
1 引 言
N是植物必需的营养元素,也是评价土壤质量和土地生产力的重要指标.为了获得高产,需要施用大量的氮肥.据统计[32],仅1996年全世界氮肥(折纯N)使用总量就高达8 50!107t,但N累积利用率不高.据估计,施入土壤中的N大约有35%通过各种途径损失掉[6,32].此外,氮肥的使用还可能造成环境污染,诸如温室气体(主要是氮氧化物)和致酸雨气体(氨气)的排放、地下水硝酸盐超标、水体富营养化等[20].如2000年,比利时80%的饮用水中硝酸盐含量超标[10].目前,土壤N循环的研究已经成为土壤学家、环境学家、农学家等共同关注的热点问题之一.
土壤N循环是N生物地球化学循环中的重要环节,其模拟是作物估产、环境评价、农田管理、决策制定和长期预测的重要依据,对提高氮肥利用率、防止或减轻环境污染具有重要的理论和实践意义.20世纪60年代,就有基于单个过程的土壤N循环方面的报道[25,28].40多年来,北美和欧洲一些国家建立了大量的土壤N循环模型.我国在这方面研究还比较薄弱[3,15,24].本文拟通过简要概述土壤N循环过程,重点介绍13种土壤N循环模型和6个土壤N循环过程的模拟,并讨论模型模拟中的参数化问题,以期为深入研究土壤N循环及其模拟提供一定的参考和借鉴.
2 土壤N循环的基本过程
土壤中含N化合物种类多,理化、生物学性质各异.一般可将土壤中N划分为有机氮和无机氮,以有机氮为主.在土壤微生物等因子的作用下,N在土壤中发生一系列复杂的循环.主要循环过程有:有机氮矿化、腐殖化、硝化、反硝化、氨挥发、N沉降、硝酸盐淋失、生物固氮、铵离子晶格固定和释放、土壤粘粒吸附和解吸、植物吸收等过程.土壤N循环过程的研究是建立土壤N循环模型以及N生物地球化学循环模型的基础.
3 土壤N循环模型的研究概况
目前,农业中数学模型并无统一的分类,可从不同角度进行划分.根据建模的方法可分为经验模型和机理模型;从土壤有机氮角度可分为单组分和多组分模型;从模拟循环过程的数目方面可分为单过程和多过程模型;此外,根据模型模拟的元素也可分为独立N模型和综合模型的N子模型.
经验模型通常依据实验测定或调查的N循环分量与气
*中国科学院知识创新工程重要方向项目(KZCX3 SW 426)、国家
自然科学基金重点项目(40235057)和国家重点基础研究发展资助项目(2002CB412503).**通讯联系人.
2005-01-10收稿,2005-05-08接受.
11期 唐国勇等:土壤氮素循环模型及其模拟研究进展 候因子等已知环境变量间的数据关系,采用统计方法,通过曲线拟合建立方程,以此进行外推,如Shaffer等[29]建立的NTRM.该模型主要用于评估土壤侵蚀对土地生产力、植物竞争、作物产量和水质的影响,较详细地模拟了水分和化学物质运移、作物生长、根系发展、田间管理效应、土壤C和N动态、径流、土壤化学平衡等方面.本模型受美国农业部大平原系统研究单位支持和推广应用.在模拟N动态方面,该模型只适用于小时空尺度内N循环过程的定量研究,主要包括对矿化、腐殖化、硝化、尿素水解、植物吸收等过程的模拟.模型运行所需的参数较少,主要参数有土壤温度、土壤氧气2209
(3)
dyt/dt=-umy/(Ks+Y)
其中,k0和k1分别为零级和一级动力学常数,um为最大转化速度,Ks为半饱和常数,dyt/dt为转化速率,Y为某N成分的含量(或浓度).
20世纪60年代后,国际上建立了大量的土壤N模型,特别是随着计算机的广泛应用,土壤N循环模型有了长足的发展.其中影响较大的有:美国的NTRM、NLEAP、NC SOIL、DNDC、EPIC、CERES、CENTURY、GLEAMS、RZWQM模型;英国的SUNDIAL模型;荷兰的SOILN模型;丹麦的DAISY模型以及德国的CANDY模型.这些模型含量、土壤水分含量、土壤pH值等.在NTRM模型的基础上,经过参数调整、修正和部分理论推导,发展成为NLEAP和RZWQM等机理模型,但土壤N硝化反应由最初的回归方程演变成机理的半阶方程[7,29].一般经验模型所需的参数较少,使用比较方便,应用较为广泛,但不能准确地反映土壤N的长期动态,有较强的地域性.机理模型则从N循环机理着手,研究各种环境因子对N循环过程的影响机理,如由英国洛桑试验站开发的SUNDIAL模型
[2]
.该模型主要用于评
估农业N循环对环境的影响.模型运行以周为步长,经过洛桑实验站长达150年定位试验的检验,证明该模型能较好地预测土壤中N长期动态,不适用于模拟土壤短期(1个月内)N动态.模型主要模拟了矿化、腐殖化、硝化、反硝化、氨挥发、植物吸收、硝酸盐淋失等过程.但氨挥发过程只模拟了肥料表施和衰老植物组分中氨的挥发,硝酸盐淋失也只模拟硝酸盐淋失至深层地下水这一过程.在有机N矿化过程中,该模型考虑了粘粒含量对矿化速率的影响.模型运行所需的参数主要包括周平均气温、田间最大持水量、最大有效持水量、土壤粘粒含量、基质组分中的C/N、各过程反应速率常数、作物目标产量等.通常机理模型理论性强,延展性好,可用于模拟土壤N长期动态,但模型运行所需参数较多,要求建模者对整个N循环过程都有深入的研究.早期土壤N模型以经验模型为主,近期则多为机理模型.
单过程模型只考虑N循环中某一过程.国内N循环模型主要为单过程模型.多过程模型研究两个或多个N循环过程以及这些过程的相互关系.N循环过程多在土壤微生物驱动下进行,而微生物活动又依赖于土壤有机碳(作为微生物能源物质)的供应.因此,大多数N模型以C/N联合模型而存在,甚至是生态系统物质循环综合模型的子模型,如CERESN子模型、CENTURYN子模型等.
在N循环过程模型中,通常以零级动力学方程、一级动力学方程或米氏方程等作为模型的数学基础,根据N循环过程的影响因子,对模型进行调整和修正.零级动力学方程的微分形式可表示为:
dyt/dt=-k0
(1)
一级动力学方程的微分形式可表示为:
dyt/dt=-k1y
(2)米氏方程的微分形式可表示为:
各有其长处和不足,已有文献对此进行了介绍和比较[9,21,24],本文仅对这13种N循环模型进行简要的介绍(表1).
4 土壤N循环过程模型模拟4 1 矿化和腐殖化过程
大多数土壤N循环模型用一级动力学方程模拟土壤有机N的净矿化或净腐殖化.在NCSOIL模型[23]中,土壤基质按C/N和土壤有机质半减期(half life)分为4个组分(残体库、微生物库、中间库和稳定腐殖质库),每个库又分为易分解和难分解两部分,每组分有机N的矿化用一级动力学方程模拟,其微分形式为:
dC/dt=min(fT,fW) N[SLkLC+(1-SL)hRC](4)其中,C是土壤基质碳含量,SL和(1-SL)分别为基质中易分解和难分解部分含量,kL和hR分别为相应的分解速率常数,fT和fW分别为温度和水分对矿化的影响因子, N为基于C/N的降低因子.土壤N净矿化或净腐殖化取决于基质残体库和微生物库C/N[23].各组分的分解速率常数、土壤温度、土壤最大持水量、基质C/N、基质中易分解与难分解比例、耕作措施等参数是模型运行所必需的,如耕作能提高土壤基质中间库易分解部分的比例[23].
影响矿化和腐殖化过程的因子很多,模型大多模拟了土壤温度、土壤含水量(或土水势)、土壤通气状况、土壤有机碳含量和土壤C/N对该过程的影响[1,2,7,9,13,14,26],也有模型考虑了土壤微生物量[22,26,33,35]
、碳磷比
[11,35]
、土壤容
重[35]、土壤质地[26]、土壤离子强度[26]、木质素含量[26]等可
能影响该过程的因子.4 2 硝化过程
硝化过程可用零级、一级或米氏方程进行模拟.在NTRM模型中,最初用经验的回归方程表示硝化速率[7]:
K=a+bTCNH4
+clog10CNH3
+dlog10CNO
3
(5)
其中,K为硝化反应速率,a、b、c、d均为常数,其值分别为4 64、1 62!10-3、0 238和-2 51,T为土壤温度,CNH4为土壤铵态氮含量,CNO3
为基质硝态氮含量.经过修正,Shaffer
等
[28]
用机理性的半阶方程来表示硝化速率:
d(NH4)/dt=[k1/2(NH4)1/2O2exp(Ea/kbTs)]/H+(6)
应 用 生 态 学 报 16卷2210
表1 13个土壤N循环模型简介
Table1Introductionto13soilnitrogencyclingmodels模 型ModelsNTRM
模型简要描述Descriptionofmodels
该模型为经验模型,主要用于评估土壤侵蚀对土壤生产力、植物竞争、作物产量和水质的影响.该模型没有模拟环境因子等对N循环过程的影响,只适用于小时空尺度内N循环过程的定量研究.该模型未将土壤有机N进行分组,用回归方法建立了矿化过程、腐殖化过程、硝化过程、尿素水解、硝酸盐淋失过程的方程;植物吸收N与吸水量成线性关系[7,29]
该模型可模拟硝酸盐的淋失,用于评估水质,能与GIS结合.该模型以天为步长,适用于有机N中长期模拟预测.模型将土壤有机N分为3个组分库,其中矿化、腐殖化、反硝化、氨挥发、硝酸盐淋失过程用一级动力学方程模拟;用logistic曲线模拟作物吸N量;用零级动力学方程模拟硝化过程[30]
该模型主要用于评估土壤N循环对环境的影响.模型运行以周为步长,适用于土壤中长期N的预测.将土壤有机N分为3个组分库,用一级动力学方程模拟矿化过程、腐殖化过程、硝化过程;用土壤有机质分解产生CO2量的线性方程表示反硝化过程;氨气挥发量与铵态肥施入量呈线性关系;植物吸收过程用植物目标含N量表示[2]
该模型主要用于评估施用N对环境的影响.早期版本中将土壤有机氮分为3个组分库.用一级动力学方程模拟矿化、腐殖化、硝化过程;用零级动力学方程模拟反硝化过程;用logistic曲线模拟植物吸N量,并模拟氨挥发过程和硝酸盐淋失过程[17]
该模型主要用于评估农家肥、泥浆等含N化合物的施用对环境的影响[14,36].该模型以日为步长,适用于短期土壤N动态的模拟.早期版本中并未对土壤有机氮进行分组[14],新版本中将有机N分为6个库[36].用一级动力学方程模拟矿化和腐殖化过程;用米氏方程模拟硝化过程;植物吸N量取决于植物潜在需求和土壤供应量;也简单地模拟了硝酸盐淋失过程[14,16]
该模型主要用于6类含N有机化合物的施用对环境的影响.模型中将土壤有机氮分为3个组分库,矿化过程、腐殖化过程用一级动力学方程模拟;反硝化过程与硝酸根浓度和有机碳量成比例;也简单地模拟了氨气挥发情况;硝化过程用米氏方程模拟[8]该模型用于模拟土壤C、N动态和痕量气体的排放.该模型将土壤有机氮分为3个组分库,跟踪有机碳、氮的逐级降解;由降解作用产生的可给态C、N,被输入硝化、脱氮及发酵子模型中,进而模拟有关微生物的活动及其代谢产物,包括几种温室气体的排放[19]该模型用于模拟土壤C、N循环,与其它模型不同的是在研究N循环机制中,该模型能同时模拟标记N和非标记N的去向.模型将土壤有机氮分成4个组分库.用米氏方程模拟矿化和腐殖化过程;用零级动力学曲线模拟硝化过程、反硝化过程;用logistic曲线模拟作物生长[23]
该模型由一系列作物生长子模型构成,各子模型中的概念、模拟方法和影响因子相同.主要用来模拟作物的生长及其与之有关的过程.模型中将土壤有机氮分为腐殖质库和新鲜有机残体库.矿化和腐殖化、反硝化、氨挥发过程用一级动力学方程模拟;硝化过程用米氏方程模拟;模型详细地模拟作物对N的吸收过程[13]
该模型主要用来模拟生态系统C、N和营养物质长期动态[26,27].早期版本中以月为步长[26],新版本中某些过程以
[27]
日为步长.由于模型运行的步长较长,模型中N子模型部分较简单.矿化和腐殖化过程用一级动力学方程模拟;植物吸N量取决于植物N需求和土壤N供应;用总矿化N的5%表示氨挥发量[26];在最新版的DAYCENT模型中,硝化过程和反硝化过程用一级动力学方程模拟[27]
该模型主要特色就是能模拟不同农业管理措施及其对水质和作物产量的影响.模型将土壤有机氮分为5个组分库.用包含微生物生长的一级动力学方程模拟矿化和腐殖化过程;用一级动力学曲线模拟反硝化过程、氨挥发过程;用米氏方程模拟作物吸N量;用零级动力学过程模拟硝化过程;模型详细地模拟了硝酸盐淋失过程[1]
该模型主要用于模拟农业管理对地下水负载的影响.模型将土壤有机氮分为3个组分库.矿化和腐殖化过程、反硝化过程、氨挥发过程用一级动力学方程模拟;硝化过程用零级动力学方程模拟;简单地模拟了土壤侵蚀和地表径流中N的损失[18]
该模型主要用于评估土壤侵蚀、耕作方式对生产力的影响.土壤有机氮的分组与GLEAMS模型中的分组一样,只是其中的活性库是估算的,而非直接测定得出.矿化过程、腐殖化过程、硝化过程、反硝化过程、氨挥发过程用一级动力学方程模拟;植物吸N量取决于植物N需求和土壤N供应;也模拟了土壤侵蚀和地表径流中N的损失[35]
NLEAP
SUNDIAL
SOILNDAISY
CANDY
DNDC
NCSOIL
CERES
CENTURY
RZWQM
GLEAMS
EPIC
其中,d(NH4)/dt为硝化速率,k1/2为半阶速率常数,Ea为化学自由能,kb为波尔茨曼常数,Ts为土壤温度,NH4为土壤溶液中铵离子浓度,O2为土壤空气中氧气浓度,H+为基质氢离子浓度.在NTRM模型中,无论使用经验的回归方程还是机理的半阶方程,土壤有机氮均不分组.回归方程中所需的参数只有土壤温度、土壤硝态氮含量和铵态氮含量,而半阶方程中所需的参数包括半阶反应速率常数、土壤铵态氮含量、土壤空气中氧气含量、土壤pH值、土壤温度.
影响硝化过程的因素很多,土壤温度、土壤含水量、通气状况、土壤pH[1,2,7,9,13,14,26]是大多数模型所考虑的主要环境因子.微生物种群动态及活性[1]、硝酸根与铵根浓度比[17]等因素对硝化过程的影响也不容忽视.一些模型考虑了硝化过程中N2O的释放情况,通常用被硝化总N的一定比例来表示N2O释放量
[7,11,26]
拟
[1,2,7,9,17,26]
.在通气良好的土壤、土表薄层水体、有氧的根
际中硝化作用都比较明显[5].4 3 反硝化过程
目前反硝化过程的模拟一般采用零级或一级动力学方程,如CERES模型用一级动力学方程模拟该过程,其一级反硝化速率系数为:
kdnit=6 0!105!CW!Tf!Wfd!h!NO3
(7)
其中,kdnit为一级反硝化速率系数,CW为土壤有机质中水提取C浓度( g#g-1),Wfd和Tf分别为水分和温度影响因子,h为土层厚度.CW、Wfd和Tf又可表示为:
CW=24 5+0 0031CTf=0 1!exp(0 046Ts)
Wfd=1 0-(SAT-SW)/(SAT-DUL)
(8)(9)(10)
.也有模型用氨的相对权重而不是
[35]
其中,C为土壤有机碳含量,SAT为土壤饱和水含量,SW为土壤含水量,DUL为土壤排干时水分含量,Ts为土壤温度.该模型认为,只有当土壤含水量超过土壤排干时水分含量的情况下,反硝化反应才发生.在目前众多土壤N循环模型中,
铵离子有效浓度来表示硝化速率.此外,硝化抑制剂的使
用能显著降低硝化反应速率[18].土壤通气状况是影响硝化过程的决定性因素,常用氧分压的非线性方程来模
11期 唐国勇等:土壤氮素循环模型及其模拟研究进展 该模型是为数不多的能用于模拟稻田土壤反硝化作用的模型之一.
多数模型模拟了土壤温度、土壤含水量、厌氧状况、土壤有机质含量(尤其是土壤水溶性C含量)对反硝化过程的影响.有模型也考虑了不同反硝化产物的情况,主要是区分了N2和N2O的不同排放量以及环境因子对N2和N2O排放的影响[2,12,26,29,31].4 4 氨挥发过程
可用零级或一级动力学方程模拟氨挥发过程,如GLEAMS模型用一级动力学方程表示该过程:
2211
定量研究硝酸盐的淋失对环境的影响是许多土壤N循环模型建立的主要目的之一,如DAISY、GLEAMS、SOILN、EPIC、NLEAP、SUNDIAL、CANDY等模型.硝酸盐淋失过程属于土壤溶质迁移的范畴.在土壤N循环模型中,常用土壤溶质迁移方面的模型或方程来模拟硝酸盐淋失.土壤N循环模型主要应用对流 扩散模型[12]、两区模型[7,26,29]、两流区模型[2]来模拟该过程.硝酸盐淋失与土壤含水量和土壤水分运动关系密切.土壤水分运动包括土壤水分入渗和再分布.土壤水分入渗常用Green Ampt[1,13,26,29]来模拟,而再分布过程用Richards方程
[29]
来拟合.国内外已有大量文献对硝酸盐
r=CNH3
exp(-Kvt)
其中,r为一级氨挥发速率,CNH3为基质中氨浓度,t为时间,Kv为速率常数,该常数又可用指数方程表示:Kv=0 409!1 08ATP-20,ATP为气温.但该方程只适用于模拟表施铵态氮肥和厩肥过程中氨气挥发情况.该模型对氨挥发过程的模拟步长为1d,模拟时空尺度较小.
影响氨挥发的因素主要是土壤温度、土壤含水量(或土水势)、氨气压梯度、风速.大多数模型只模拟了这几个主要因素[1,2,13,26,29,35].也有模型考虑了土壤阳离子交换量[2]、土壤表面覆盖
[13]
、肥料(铵态氮肥和尿素)的施用方式
[13]
、土层
深度[35]等.也有模型用气温代替土温来模拟温度对氨挥发的影响[18].大多数模型认为氨挥发只发生在土表,也有模型认为整个土体中都存在氨挥发[2].
与氨挥发相反,N沉降一般不被建模者所重视.事实上,在正常大气氨浓度下,N沉降中能被植物有效利用的N占植物N需求量的10%左右[16].而在局部高浓度氨的大气中,该机制所占的份额更大,如氮肥厂、化工厂等氨高排放区域[34].
4 5 NH4+的晶格固定和释放过程
土壤NH4+的晶格固定和释放是土壤N有效性研究中一个备受关注的问题.研究者
[22,37]
认为,
NH+4
的晶格固定
和释放类似于土壤吸附和解吸过程,可用等温吸附方程来描述,如Freundlich、Langmuir方程.影响该过程的因素主要包括土壤有机质含量、土壤溶液中
NH+4
浓度、土壤溶液中钾离
子浓度、土壤含水量、土壤矿物种类、土壤粘粒含量等.但这些因素对铵晶格固定和释放过程的影响难以用单值函数描述.如钾离子在低浓度时能促进NH+4
的固定,而高浓度时又
可能会抑制
NH+4
的固定
[22,37]
.土壤中蛭石、伊利石、蒙脱石
等2∀1型层状硅酸盐矿物含量高时,能加强
NH+4
的晶格固
定和释放,也能增加土壤固定NH+4量[4,22,37]
.NH+4的释放
除了与上述因素有关外,最重要的影响因素是NH+4固定的时间.一般可将固定态NH+4分为∃新固定的%和∃旧固定的%.∃新固定的%NH+4较易重新释放,而随固定时间的推移,固定
态NH+4就难以再释放出来
[4,22,37].目前,对土壤NH4+晶格固定和释放的研究还不够深入,还很难定量模拟环境因子对该过程的影响.4 6 硝酸盐淋失过程
淋失过程进行了报道和介绍[3,6,10,17,30,31],故本文不再论述.5 问题与展望
目前只有较少的土壤N循环模型在田间条件下得到比较充分的验证.其原因之一是模型参数化过程中参数的不确定性,而且这些不确定性对模型影响程度的评估也存在一定的困难.这种不确定性包括3个方面:1)模型参数求算过程的不确定性,尤其是不可直接测量的内部参数,如分解速率常数.这些参数只能通过其他数据资料进行计算.在计算过程中大多需要设定初始值或初始状态,而初始值或初始状态的设定通常是不确定的.2)测定资料的不确定性.土壤、大气中各成分的时空变异性很大,而由于测量技术和实验仪器的限制,测定的数据并不能完全反映这种变异性.由于测定资料的不准确,通过逐级放大作用,最终影响参数的准确性.3)过程间相互作用的不确定性.影响土壤N动态的因素很多,模型并不可能考虑所有的影响因素,有时只能用∃黑箱%表示,而且N过程间也存在着相互作用.此外,土壤N循环模拟是基于单个过程的,缺乏系统性.因此,对用户来说,好的模型最重要的要素就是:模型能一步一步地指导用户如何获取模型参数.
近年来,土壤N循环模型出现了两个大的发展趋势,即理论模型趋于复杂化、细微化,应用模型趋于简单化,易于用户使用.
从理论研究方面来看,土壤N循环模型将更加突出系统性,强调各过程间的相互联系.从应用方面来看,∃3S%(地理信息系统GIS,全球定位系统GPS,遥感RS)技术与土壤N循环模型相结合将是土壤N循环模型发展的趋势.目前已有较多的模型和GIS系统相结合,如CERES GIS,NLEAP GIS,RZWQM GIS,且CERES和RZWQM已形成专家系统.参考文献
1 AhujaLR,RojasKW,HansonJD,etal.2000.RootZoneWater
QualityModel:ModellingManagementEffectsonWaterQualityandCropProduction.LLC,HighlandsRanch,CO:WaterResourcesPublication.372
2 BardburyNJ,WhitmoreAP,HartPBS,etal.1993.Modellingthefateofnitrogenincropandsoilintheyearsfollowingapplicationof15
N labelledfertilizertowinterwheat.JAgricSci,121:363~3793 CuiJ B(崔剑波),ZhuangJ P(庄季屏).1997.SoilNO3- Ntrans portunderfieldunsaturatedflowcondition&Asimulationstudy.ChinJApplEcol(应用生态学报),8(1):49~54(inChinese)
应 用 生 态 学 报 16卷2212
4 DavidsonJB,AldagRW.1982.Inorganicformsofnitrogeninsoils.In:StevensonFJ,eds.NitrogeninAgriculturalSoils.Madison,WI:AmericanSocietyofAgronomy,Inc.48
5 DeBuskWF,WhiteJR,ReddyKR.2001.Carbonandnitrogendy
namicsinwetlandsoils.In:ShafferMJ,LiwangM,HansenS,eds.ModellingCarbonandNitrogenDynamicsforSoilManagement.Bo caRaton,FL:LewisPublishers.27~54
6 DelgadoJA,FollettRJ,ShafferMJ.2000.SimulationofNO3 Ndy
namicsforcroppingsystemswithdifferentrootingdepths.SoilSciSocAmJ,64:1050~1054
7 DuttRG,ShafferMJ,MooreWJ.1972.Computersimulationmodel
ofdynamicbio physicochemicalprocessesinsoils.ArizonaAgricExpStaTechnolBull,196:128
8 FrankoU,OelschagelB,SchenkS.1995.Simulationoftempera
ture,water,andnitrogendynamicsusingthemodelCANDY.EcolMod,81:213~322
9 FrisselMJ,VanVeenJA.1978.Acritiqueofcomputersimulation
modellingfornitrogeninirrigatedcroplands.In:NielsenDR,Mac DonaldJG,eds.NitrogenintheEnvironment.NewYork:AcademicPress.145~162
10 GarcetJDP,AmauryT,JavauxM,etal.2001.Modellingnitrogen
behaviourinthesoilandvadoseenvironmentsupportingfertilizermanagementatthefarmscale.In:ShafferMJ,LiwangM,HansenS,eds.ModellingCarbonandNitrogenDynamicsforSoilManage ment.BocaRaton,FL:LewisPublishers.571~596
11 GijsmanAJ,ObersonA,TiessenH,etal.1996.Limitedapplicabili
tyoftheCENTURYmodeltohighlyweatheredtropicalsoils.A gron,88:894~903
12 GrantRF,PatteyE.1999.Mathematicalmodellingofnitrousoxide
emissionfromanagriculturalfieldduringspringthaw.GlobalBio geochemCyc,13:679~694
13 HanksJ,RitchieJT.1991.ModellingPlantandSoilSystem.Madi
son,WI:AmericanSocietyofAgronomy,Inc.
14 HansenS,JensenHE,NielsenNE,etal.1991.Simulationofnitro
gendynamicsandbiomassproductioninwinterwheatusingtheDanishsimulationmodelDAISY.FertRes,27:245~259
15 HuangY(黄 耀).2003.CarbonandNitrogenExchangeinEarth
atmosphericSystem&FromExperimentstoModels.Beijing:Meteo rologyPress.(inChinese)
16 HutchinsonGL,MillingtonRJ,PetersDB.1972.Atmosphericam
monia:Absorptionbyplantleaves.Science,175:771~772
17 JohnsonHL,BergstromPE,Jansson,etal.1987.Simulatednitro
gendynamicsandlossesinalayeredagriculturalsoil.AgricEcosystEnviron,18:333~356
18 LeonardRA,KinselWG,StillDS.1987.GLEAMS:Groundwater
loadingeffectsofagriculturalmanagementsystems.TransASAE,30:1403~141819 LiC S(李长生).2001.Biogeochemicalconceptsandmethodolo
gies:DevelopmentoftheDNDCmodel.QuatSci(第四纪研究),21(3):89~99(inChinese)
20 LiwangM,ShafferMJ.2001.Areviewofcarbonandnitrogenpro
cessesinnineU.S.soilnitrogendynamicsmodels.In:ShafferMJ,LiwangM,HansenS,eds.ModellingCarbonandNitrogenDynam icsforSoilManagement.BocaRaton,FL:LewisPublishers.55~102
21 McGechanMB,WuL.2001.Areviewofcarbonandnitrogenin
Europeansoilnitrogendynamicsmodels.In:ShafferMJ,LiwangM,HansenS,eds.ModellingCarbonandNitrogenDynamicsforSoil
Management.BocaRaton,FL:LewisPublishers.103~171
22
McGillWB,HuntHW,WoodmanseeRG,etal.1981.Phoenix,amodelofthedynamicsofcarbonandnitrogeningrasslandsoils.EcolBull,33:49~115
23
MolinaJAE,ClappCE,ShafferMJ,etal.1983.NCSOIL,amodelofnitrogenandcarbontransformationsinsoil:Description,calibra tion,andbehaviour.SoilSciSocAmJ,47:85~9124
MuX M(穆兴民),FanS L(樊小林).1999.AreviewonecologymodelsofsoilNmineralization.ChinJApplEcol(应用生态学报),10(1):114~118(inChinese)
25 OlsonJS.1963.Energystorageandthebalanceofproducersandde composersinecologicalsystem.Ecology,44:322~331
26
PartonWJ,MosierAR,OjimaDS,etal.1996.GeneralizedmodelforN2andN2Oproductionfromnitrificationanddenitrification.GlobalBiogeochemCyc,10:401~412
27
PartonWJ,OjmaDS,ColCV,etal.1994.AGeneralmodelforsoilorganicmatterdynamics:Sensitivitytolitterchemistry,texture,andmanagement.SoilSciSocAmSpecPull,39:147~167
28
SabeyBR,FrederikLR,BarthdomewWV.1969.TheformationofnitratefromNH+4 Hinsoils∋.Useofthedelayandmaximumratephaseformakingquantitativepredictions.SoilSciSocAmJ,33:276~278
29
ShafferMJ,LarsonWE.1987.NTRM,Asoil cropsimulationmod elfornitrogen,tillage,andcrop residuemanagement.USDAConsResRep,34(1):103
30
ShafferMJ,HalvorsonAD,PierceFJ.1991.Nitrateleachingande conomicanalysispackage(NLEAP):Modeldescriptionandapplica tion.In:FollettRF,ed.ManagingNitrogenforGroundwaterQuali tyandFarmProfitability.Madison,WI:SoilScienceSocietyofAmerica,Inc.285~322
31
ShafferMJ,LasnikK,OuX,etal.2001.NLEAPInternettoolsforestimatingNO3 NleachingandN2Oemission.In:ShafferMJ,Li wangM,HansenS,eds.ModellingCarbonandNitrogenDynamicsforSoilManagement.BocaRaton,FL:LewisPublishers.403~42632 SmilV.1997.Nitrogenincropproduction:Anaccountofglobalflows.GlobalBiogeochemCyc,13:647~662
33 SmithOL.1982.SoilMicrobiology:AModelofDecompositionandNutrientCycling.BocaRaton,FL:CRCPress.273
34
WhiteheadDC,LockyerDR.1987.Theinfluenceoftheconcentra tionofgaseousammoniaonitsuptakebytheleavesofItalianRye grass,withandwithoutanadequatesupplyofnitrogentotheroots.ExpBot,38:818~827
35
WilliamsJR.1995.TheEPICmodel.In:SinghVP,ed.ComputerModelsofWatershedHydrology.HighlandsRanch,CO:WaterRe sourcesPublications.909~100036
XuC,ShafferMJ,Al KaisiM.1998.Simulatingtheimpactofman agementpracticesonnitrousoxideemissions.SoilSciSocAmJ,62:736~724
37
ZhuZ L(朱兆良),WenQ X(文启孝),ChengL L(程励励),etal.1992.NitrogeninChineseSoils.Nanjing:JiangsuScienceandTechnologyPress.(inChinese)
作者简介 唐国勇,男,1980年生,硕士研究生.主要从事土壤养分循环及模拟研究,发表论文2篇.E mail:[email protected].