Calculating actual crop evapotranspiration under soil water stress conditions with appropriate numerical methods and time step     

Songhao Shang 《水文研究》2012,26(22):3338-3343

Calculation of actual crop evapotranspiration under soil water stress conditions is crucial for hydrological modeling and irrigation water management. Results of actual evapotranspiration depend on the estimation of water stress coefficient from soil water storage in the root zone, which varies with numerical methods and time step used. During soil water depletion periods without irrigation or precipitation, the actual crop evapotranspiration can be calculated by an analytical method and various numerical methods. We compared the results from several commonly used numerical methods, including the explicit, implicit and modified Euler methods, the midpoint method, and the Heun's third‐order method, with results of the analytical method as the bench mark. Results indicate that relative errors of actual crop evapotranspiration calculated with numerical methods in one time step are independent of the initial soil water storage in the range of soil water stress. Absolute values of relative error decrease with the order of numerical methods. They also decrease with the number of time step, which can ensure the numerical stability of successive simulation of soil water balance. Considering the calculation complexity and calculation errors caused by numerical approximation for different time step and maximum crop evapotranspiration, the explicit Euler method is recommended for the time step of 1 day (d) or 2 d for maximum crop evapotranspiration less than 5 mm/d, the midpoint method or the modified Euler method for the time step of up to one week or 10 d for maximum crop evapotranspiration less than 5 mm/d, and the Heun's third‐order method for the time step of up to 15 d. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

内蒙正镶白旗岩体中长石的初步研究     

姚培琳  饶明辉 《华东地质学院学报》1989,12(1):7-17

  相似文献   

利用位场及其高阶导数的功率谱计算岩层的深度与厚度   总被引：3，自引：2，他引：3       下载免费PDF全文

侯重初  李保国 《地球物理学报》1988,31(1):90-98

本文导出了利用重力与磁异常及其高阶导数的功率谱计算密度体和磁性体上顶与下底埋深的近似公式,还提供了利用垂向或水平导数谱计算深度的方法。理论模型的试验结果表明,由近似公式计算的精度可以满足实际需要。本文最后根据一条重力剖面资料,计算了密度层上、下界面。  相似文献   

我国月降水和气温网格点资料的处理和分析   总被引：14，自引：2，他引：14  

魏凤英  曹鸿兴 《气象》1994,20(10):26-30

用逐步搜索逼近法的客观分析技术计算出１９５１－１９９２年中国月降水量和月平均气温的格点资料，并利用图象识别原理对格点值和站点观测值进行了比较。  相似文献   

Oscillatory Rossby Solitary Waves in the Atmosphere     

Chen Zhongming  Liu Fuming  Li Xiaoping  Tao Jie 《大气科学进展》1994,11(1):65-73

The linear Rossby wave frequency expression is expanded at higher accuracy based on the scale difference char?acteristics of atmospheric long waves in the and directions. That the nature of the waves represented by the expan?sion is identical to that of the original ones is demonstrated both in phase velocity and wave energy dispersion speed , followed by the derivation of the nonlinear expression describing atmospheric long wave behaviors with the associated approximate analytic solution obtained. Then, for the first time atmospheric’ oscillatory Rossby solitary wave’ with its dispersion relation is obtained by numerical calculation with the aid of physical parameters of the real atmosphere. The solitary wave is found to be very close to such longwave systems as blocking highs and cut-off de?pressions in the actual atmosphere.  相似文献   

纬向切变基流中的非线性Rossby波     

何建中  郭品文 《高原气象》1993,12(4):361-366

本文在半地转近似下讨论了纬向切变基流中的非线性正压Ｒｏｓｓｂｙ波，给出了存在波解的切变基流条件。  相似文献   

Waveforms of Long-Period Body Waves In A Slightly Aspherical Earth Model     

Xiang-Dong Li  Toshiro Tanimoto 《Geophysical Journal International》1993,112(1):92-102

  相似文献   

Research　of　pattern　dynamics　parameters　of　crustal　deformation　field　in　seismogenic　process     

周硕愚  吴云  王若柏  杨国华 《地震学报(英文版)》1994,7(3):427-432

ＲｅｓｅａｒｃｈｏｆｐａｔｔｅｒｎｄｙｎａｍｉｃｓｐａｒａｍｅｔｅｒｓｏｆｃｒｕｓｔａｌｄｅｆｏｒｍａｔｉｏｎｆｉｅｌｄｉｎｓｅｉｓｍｏｇｅｎｉｃｐｒｏｃｅｓｓＳｈｕｏ－ＹｕＺＨＯＵ（周硕愚），ＹｕｎＷＵ（吴云），Ｒｕｏ－ＢｏＷＡＮＧ（王若柏）ａｎｄ...  相似文献   

一种新的计算冰面饱和水汽压方法   总被引：1，自引：0，他引：1  

刘公波  胡志晋 《气象》1994,20(1):24-26

综述了国内外近似计算平面饱和水汽压公式。给出了较为精确、简单的平冰面饱和水汽压计算式，并与Ｗｅｘｌｅｒ等的计算结果进行了比较。  相似文献   

计算最小走时和射线路径的界面网全局方法   总被引：27，自引：16，他引：27       下载免费PDF全文

刘洪  孟凡林 《地球物理学报》1995,38(6):823-832

用慢度分块均匀正方形模型将介质参数化，仅在正方形单元的边界上设置计算结点，这些结点构成界面网．根据Ｈｕｖｓｅｎｓ和Ｆｅｒｍａｔ原理，由不断扩张、收缩的波前点扫描代替波前面搜索，在波前点附近点的局部最小走时计算中对波前点之间的走时使用双曲线近似，通过比较确定最小走时和相应的次级源位置，记录在以界面网点位置为指针的３个一维数组中．借助这些数组通过向源搜索可计算任意点（包括界面网以外的点）上的全局最小走时和射线路径．这一方法不受介质慢度差异大小限制，占内存少，计算速度较快，适于走时反演和以Ｍａｓｌｏｖ射线理论为基础的波场计算．  相似文献