土壤水动力学SWD7SPAC水分传输.ppt

土壤水动力学(SoilWaterDynamics)第7章土壤植物大气连续体水分传输WaterTransportinSPAC(Soil-Plant-AtmosphereContinuum),毛晓敏中国农业大学水利与土木工程学院Telmail:maoxiaomin,第7章土壤植物大气连续体水分传输,SPAC的概念及其中的水分传输过程腾发量的估算方法SPAC水分传输模型及应用SPAC水热传输模型及应用,7.1SPAC的概念及其中的水分传输过程,水分对植物生长的作用植物生长条件下的主要水分传输过程SPAC的概念,7.1.1水分对植物生长的作用,植物细胞原生质(细胞膜、细胞质和细胞核)含水量达80%以上时，才能进行活跃的生化反映，保持代谢,光合作用光能6CO2+6H2OC6H12O6+6O2叶绿素营养物质的输送,蒸腾：植物吸收水分的绝大部分消耗于蒸腾蒸腾的作用：植物吸收、传输水分、养分的驱动力：通过蒸腾作用，维持植物体内及其与土壤的水势梯度，从而不断地吸收、传输水分、养分调节小气候，降低叶面温度，避免叶片灼伤,7.1.2植物生长条件下的主要水分传输过程,水分传输过程：土壤根表面：土壤物理学根表面根导管：植物生理学根导管叶面：植物生理学叶面大气：气象学,水流类型：液体水流：水在土壤、植物体的流动水汽扩散：干土层、大气水的相变：蒸发、蒸腾腾发,7.1.3SPAC的概念,早期对各种水分传输过程单独研究，概念不统一后来，用水势的概念综合研究水分传输，将土壤植物大气视为连续体，Philip于1966年提出SPAC的概念土壤植物大气连续体：Soil-Plant-AtmosphereContinuum,SPAC系统的三个重要环节：根系吸水过程水在植物体内的输运过程植物体内的水分向大气扩散过程,土根：q1=-1/R1根吸收：q2=-2/R2根叶：q3=-3/R3叶气：q4=-4/R4,土壤,根表面,根导管,叶面,大气,R1,R2,R3,R4,1,2,3,4,7.2腾发量的估算方法,腾发量及其影响因素腾发量的测定方法估算腾发量的理论方法估算腾发量的经验方法估算腾发量的实用方法参考作物法根据水量平衡模型估算腾发量,7.2.1腾发量及其影响因素,腾发、蒸散发：Evapotranpiration蒸发：Evaporation蒸腾(Transpiration)：主要通过叶片气孔散失水分在气孔腔和细胞间隙的叶肉表面蒸发（水面蒸发），水汽达到饱和水汽由气孔腔经气孔向叶面气层扩散水汽从叶面气层向大气扩散,叶片构造示意图,腾发的影响因素能量：辐射，潜热消耗L=2.5MJ/kg水汽输送：风速、温度梯度、湿度梯度其它影响因素：植物生长状况、土壤供水状况,7.2.2腾发量的测定方法,直接测定法：根据实测水分剖面，利用零(定位)通量法计算称重式蒸渗仪涡度相关法间接测定法：Bowen比法空气动力学法。,涡度相关法（直接测定潜热、显热通量）,涡度相关是指某种物质的垂直通量，即这种物质的浓度与其垂直速度的协方差。涡度相关法已经成为直接测定大气与植物群落气体交换通量的通用标准方法。该试验区下垫面平坦均匀，符合涡度相关法观测要求。运用涡度相关法开路系统直接测定潜热通量和显热通量，通过用实时所测的垂直风速与水汽浓度和温度的协方差所得，其湍流通量的计算式可表示为：,Bowen比法（测定两个高度处的温、湿度）,测定两个高度处的温度和湿度。由：,7.2.3估算腾发量的经验方法,主要用于土壤供水充分条件下作物最大腾发量（作物需水量）或参考作物腾发量的估算利用蒸发皿资料估算腾发量ETmi=iE0,利用单一气象因素估算腾发量辐射：LEm=a1Rs+a2积温：ETmi=HiSti需水系数积温气温：Blaney-Criddle公式：ETmi=KiPiTi/100以月为单位，耗水系数、月昼长占年比例、平均气温Thornthwaite公式：水汽压：ETmi=(es-e)=d,利用多个气象因素估算腾发量Dalton型公式：ETmi=(es-e)f(u)Christiansen-Hargreaves公式：ETm=0.324RsCTCuChCsCL辐射，修正：温度、风速、相对时度、日照百分比、高度,7.2.4估算腾发量的理论方法,紊流扩散法（空气动力学法）能量平衡法综合法,紊流扩散法：空气动力学法和质量迁移法。由于近地表的空气处于紊流状态，空气质点的流径相互不平行，且以极不规则的方式相互掺混。由于空气紊流动这种特点，能将其中的微粒及运动属性迅速弥散。其速率远大于分子扩散的速率。通过对紊流点平均速度、脉动速度、紊流输送到分析，得到：根据2个高度处的有关气象资料，确定水汽、显热的垂直输送通量风速水汽压气温,能量平衡法：达到地面的短波辐射：Rs=Ra(a+bn/N)净辐射：Rn=Rs(1-r)-RL能量平衡：Rn=C+E+GBowen比：=C/(E)潜热通量：E=(Rn-G)/(1+),综合法（1）：能量平衡与紊流扩散理论结合，只需一个高度的气象资料Penman公式：水面湿润表面其中Ea常称为干燥力，反映了空气湿度与饱和湿度差引起的蒸发量。,综合法（2）Penman-Monteith公式：引入冠层气孔阻力，可用于计算腾发其中ra空气边界层阻力，rc冠层气孔阻力应用中，气孔阻力测定问题,7.2.5估算腾发量的实用方法参考作物法,Allen,R.G.,Pereira,L.S.,Raes,D.,Smith,M.,1998.Cropevapotranspiration-Guidelinesforcomputingcropwaterrequirements.FAOIrrigationanddrainagepaper56,单作物系数法：ETm=KcET0,ET=KsETm，ET=KsKcET0ET0参考作物腾发量(ReferenceET)Kc作物系数(Cropcoefficient)Ks土壤供水系数ETm土壤供水充分时的最大腾发量（作物需水量）ET实际腾发量双作物系数法：,参考作物腾发量(ReferenceET)概念：参照作物腾发量为一种假想的参照作物冠层（作物高度为0.12m，固定的叶面阻力为70s/m，反射率为0.23，非常类似于表面开阔、高度一致、生长旺盛、完全遮盖地面而不缺水的绿色草地）的腾发速率计算公式：FAOPenman-Monteith公式,作物系数：与作物种类、品种、生长阶段、生育状况等因素有关，综合反映了实际作物与参考作物的差异FAO概化的分段作物系数曲线：越冬作物非越冬作物经验曲线：,土壤供水系数：反映土壤水分胁迫对腾发的影响，与根系层土壤含水率的大小及分布、最大腾发强度有关FAO方法：经验公式：如,7.2.6根据水量平衡模型估算腾发量,以上腾发量估算中，计算Ks需要根系层含水量资料，可根据水量平衡模型来推算：W=W2-W1=P+I-ET-Q-R,7.3SPAC水分传输模型及应用（恒温条件下土壤水分运动的模拟）,初始条件：初始含水率（基质势）分布已知，即：(z,t)=0(z)t=0边界条件:上边界为以地表土面蒸发强度控制的通量边界下边界视具体情况而定。,植被资料气象资料ET0ETmLAITmEm根系层水分状况地表水分状况s(z,t),TaEaRichards方程,7.4SPAC水热传输模型及应用,SPAC模型中冠层的处理：大叶（单源）模式：冠层与地面的温、湿状况相同，不能区分蒸发、蒸腾，可用于矮草双源模式：将冠层与地面分别考虑，可区分出蒸发、蒸腾，应用较为方便层状模型、块状模型、混合型模型多源模式：将冠层分为多层，能详细描述冠层温、湿垂向分布，但参数多、计算复杂,大叶模式,层状模式,块状模式,混合模式,7.4.1SPAC水热传输模型,能量分配Rn=Rv+RsRs=Rnexp(-kLAI)Rv=LEv+CvRs=LEs+Cs+GC=Cv+CsLE=LEv+LEs,土壤、冠层水热传输,水热传输阻力,水分胁迫条件下蒸腾量：,根区土水势修正系数：,根系吸水率S(t)：,根系吸水分布函数s(z,t)：s(z,t)=Evp(z,t)(,),土水势应力函数(,)：,根系密度分布(z,t)：(z,t)=(4m-1)/zr(t)(8m-4)z/zr(t)2(其它根系密度函数),结果：,SPAC水热传输模型：,地表以上：Rs(T1,)=LEs(T1,eb)+Cs(T1,Tb)+GRv(T1,)=LEv(Tv,eb,)+Cv(Tv,Tb)C(Tb)=Cv(Tv,Tb)+Cs(T1,Tb)LE(eb)=LEs(T1,eb)+LEv(Tv,eb,)未知量：冠层Tb、eb；叶片Tv、地表热通量G地表以下土壤层：,7