引 言农田土壤水分是作物生产所必需的基本要素之一。自然降水、人工灌溉和地下水补给 ,都首先要转化成土壤水才能被作物吸收利用。在土壤水分动态模拟和预报方面国内外已做了一些工作<1～ 8> ,以前的土壤水分预报模型中多数没有考虑潜水蒸发。但对地下水埋深浅 ,雨水充沛的江淮流域 ,潜水蒸发对作物根层的土壤水分动态变化具有不可忽视的影响。确定潜水蒸发量的方法有器测法、水量平衡法和经验公式<9,1 0 > 。器测法通常用蒸渗仪<1 1 > ;水量平衡法的各项分量的确定比较复杂 ,通常采用简化的水量平衡方程式 ,即Eg=μwΔhg(式中Eg 为潜水蒸发量 ,μw 当水位下降时为给水度 ,当水位上升时为饱和差 ,Δhg 为地下水位变化量 ) ,用这个方程式可以计算过去某个时段内的潜水蒸发量 ,要预报未来某个时段内的潜水蒸发量 ,首先要建立地下水位变化量预报模式 ;常用的经验公式有怀特公式和阿维扬诺夫公式 ,也有一些通过试验建立的修正公式<1 2 > ,这些公式中的参数一般都不是常规观测资料 ,通用性不强。本文利用农业气象试验站常规观测资料 ,用统计方法和水量平衡法推导出江淮地区潜水蒸发经验计算模型 ;利用农田水分平衡原理<1 3> 、引用FAO(联合国粮农组织 )的一些结果<1 4 > ,分别在江淮地区建立了引入潜水蒸发量和没有引入潜水蒸发量的冬小麦和大豆土壤水分动态预报模型 ,并对这两种模型在地下水浅埋条件下进行预报准确度比较。本模型的特点在于 :( 1 )可与其它气候模式连接 ,利用气候模式输出的气象要素预报值 ,进行实时冬小麦和大豆土壤水分预报 ;( 2 )可逐日预报冬小麦和大豆全生育期内的土壤水分。1 资料来源及其处理研究中所用资料多数是常规观测资料 ,其中逐日降水量、蒸发量、平均气温、最高气温、最低气温、日照、水汽压和风速来自气象站地面观测资料 ;土壤湿度、作物发育期和地下水位资料来自已有的土壤湿度数据库 ;土壤水文常数来自农业气象站的农业气象报表。作物系数引用参考文献 <1 4>,可能蒸散量的计算采用FAOPenman Monteith公式<1 5 ,1 6> 。2 江淮地区潜水蒸发计算模型潜水蒸发 ,指蒸发的水源直接来自潜水。潜水蒸发量的大小 ,反映在地下水位下降的幅度上。当土壤含水量减少时 ,通过土壤的非饱和渗透作用 ,潜水补充土壤水分 ,因此 ,潜水蒸发是影响作物根层土壤水分动态变化的一个因素。潜水蒸发主要受地下水埋深、气象、植物和土壤等因素的影响。为了简化计算手续 ,本文采用在无降水补给或日降水量小于 2mm、无灌溉水下渗、无开采和无侧向水分输入情况下建立的简化的潜水蒸发量计算公式 ,即简化水量平衡方程<1 7> :Eg =μwΔhg ( 1 )式 ( 1 )中Eg 为潜水蒸发量 ,μw为给水度 (因为只考虑水位下降时段 ) ,Δhg 为地下水消退深度。地下水消退深度主要受气象条件和地下水埋藏条件的影响 ,前者可以用水面蒸发强度 (E)代表 ,后者可用地下水位 (hg)代表。因此可以建立ΔhgE 与hg 的线性相关关系<9> :ΔhgE =a+b×hg ( 2 )用最小二乘法可以求得当地的a ,b值。由公式 ( 2 )得到地下水消退深度计算式 :Δhg =E(a +b×hg) ( 3)把公式 ( 3)代入公式 ( 1 )得到下式 :Eg =μwE(a+b×hg) ( 4 )公式 ( 4 )即是潜水蒸发量计算模型。根据宿县农业气象站 1 961～ 1 992年常规逐日气象资料和逐旬地下水位资料 ,从中选出 68个无降水补给或日降水量小于 2mm、无灌溉水下渗、无开采和地下水位呈消退趋势的时段作为样本 ,用最小二乘法求得当地的a ,b值分别为 1 8.471 0和 - 0 .0 641 8,相关系数为 - 0 .638。根据相关分析法 ,一个地区的给水度( μw)可用式 ( 2 )中的a值计算得到 ,即 μw=1 /a(此式的具体推导见参考文献 <9,1 7>) ,因此得到宿县地区的给水度为 0 .0 5 4 ,并建立潜水蒸发量计算模型 :Eg =μwE( 1 8.471 0 - 0 .0 641 8hg) ( 5 )用同样的方法求得阜阳 (资料年代为 1 980～ 1 992年 )的a ,b值分别为 1 3.5 849和- 0 .0 5 95 6,样本量 5 3,相关系数 - 0 .665 9,给水度为 0 .0 74,潜水蒸发量计算模型为 :Eg =μwE( 1 3.5 849- 0 .0 5 95 6hg) ( 6)式 ( 5 )和式 ( 6)中符号的物理意义同上 ,hg 的单位为cm ,Eg 和E为mm ,μw为无量纲。本文由于资料有限 ,只建立了宿县和阜阳两个站的潜水蒸发量计算模型。其它地区可以根据本文提供的思路和方法 ,用当地常规气象和地下水位观测资料建立自己的潜水蒸发量计算模型。3 冬小麦和大豆土壤水分 动态预报模型的构建和比较影响农田土壤水分收支平衡的因素主要有降水、实际蒸散、灌溉、截留、径流、渗漏和潜水蒸发等。这些因素在自然界中或在人为干预下不断的进行着变化,因此农田土壤水分在时间和空间上始终处于动态变化之中。虽然这种变化是极其复杂的 ,但通过研究可以发现一些能为人们所掌握、并服务于农业生产的有规律性的变化机制。这也是本研究的目的。3.1 土壤水分动态预报模型为了简化土壤水分动态预报模型 ,认为作物截留量与其它土壤水分收支量比起来可以忽略不计。如果某一天的初始土壤湿度加上当天的降水量超过田间持水量 ,就把超过部分当作径流量或渗漏量去除 ,不参加下一步计算。经过上述简化后 ,根据农田水分平衡原理建立江淮地区 0～ 1 0 0cm土层土壤水分动态预报模型如下 :dWdt =-Et(u ,v ,Ea,T ,Kc,W0 ) +Eg(E ,hg) +P+I ( 7)W(i) =W(i- 1 ) -Et(i) Δt +Eg(i) Δt +P(i) Δt +I(i) Δt ( 8)式 ( 7)和式 ( 8)中 ,W为土壤含水量 (mm) ,Et 为日实际蒸散量 (mm) ,Eg 为日潜水蒸发量(mm) ,P为日降水量 (mm) ,I为日灌溉量 (mm) ,i为预报时段 (n天 )中的某一天 ,Δt为时间步长 (这里取Δt =1天 )。如果忽略不计潜水蒸发量 ,土壤水分动态预报模型为 :W(i) =W(i- 1 ) -Et(i) Δt+P(i) Δt+I(i) Δt ( 9)式 ( 9)中各符号的物理意义和量纲同式 ( 8)。式 ( 8)和式 ( 9)分别是引入潜水蒸发量和没有引入潜水蒸发量的土壤水分动态预报模型。3.2 可能蒸散子模型可能蒸散量的计算采用FAOPenman Monteith公式<1 4 ,1 5> :0 .40 8Δ(Rn -G) +γ 90 0T+ 2 73U2 (es-ea)Δ +γ( 1 + 0 .34U2 ) ( 1 0 )式中 ,E0 :参考作物蒸散量 (mm·d- 1 ) ,Rn:地表净辐射 (MJ·m- 2 ·d- 1 ) ,G :土壤热通量(MJ·m- 2 ·d- 1 ,在逐日计算公式中 ,G≈ 0 ) ,T :2m高度处平均气温 (℃ ) ,U2 :2m高度处风速 (m·s- 1 ) ,es:饱和水汽压 (kPa) ,ea:实际水汽压 (kPa) ,Δ :饱和水汽压曲线斜率 (kPa·℃ - 1 ) ,γ :干湿表常数 (kPa·℃ - 1 )。3.3 实际蒸散子模型实际蒸散量是本模型水分平衡公式中的主要支出项。计算分土壤充分供水和土壤供水不足两种情况 :Et(i) =E0 (i) Kc(i) 当W(i) ≥Wk(i) (i =1 ,… ,n) ( 1 1 )Et(i) =(W(i) -Wp(i) ) <1 -exp( - E0 (i) Kc(i)Wk(i) -Wp(i)) > 当W(i)
。3.4 潜水蒸发子模型公式 ( 4 )是潜水蒸发量计算模型 ,用公式 ( 4 )计算得到的潜水蒸发量补充给地下水水面以上的整层土壤 ,而且还有一部分通过土壤表面进入大气层<1 0 > ,0～ 1 0 0cm土层接受到的只是潜水蒸发量的一部分。为了定量表示潜水蒸发量对 0～ 1 0 0cm土层的补充量 ,本文假定 2 0 %的潜水蒸发量直接进入大气。另外 ,假定潜水蒸发量在地下水水面以上的整层土壤的分配是均匀的 ,由此建立 0～ 1 0 0cm土层潜水蒸发模型 :Eg =0 .8×1 0 0hg× μwE(a +b×hg) ( 1 3)式 ( 1 3)中a、b为待定参数 ,本文只给出了宿县和阜阳两个地区的a、b值 ,其它地区可用本文提供的方法求得 ,其它各符号的物理意义和量纲同式 ( 5 )。3.5 2种模型准确度比较用公式 ( 8)和公式 ( 9)分别建立宿县单站引入潜水蒸发量和没有引入潜水蒸发量的冬小麦和大豆土壤水分动态预报模型。为了检验 2种模型的预测准确度 ,选了安徽省宿县地区 1 980年的土壤湿度资料作为实测值。在冬小麦主要生长季选了 3个时段 ,大豆生长季选了 4个时段 ,每个时段 1 0天 ,即预报时效为 1 0天 ,地下水位呈下降趋势 ,时段内逐日实际地下水位和 0～ 1 0 0cm土壤湿度用时段初和时段末的实测值内插得到。引入潜水蒸发量预报模型输出的结果为预报值 1 ,
More abstracts about the 地下水浅埋条件下冬小麦和大豆土壤水分动态预报模型研究