首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Based on the approaches of FAO24 and FMP1, a new empirical approach simulating evapotranspiration from groundwater (ETgw) under different soil water conditions was proposed in present study. Then, through the data from volumetric lysimeter, the new approach was compared with the two approaches of FMP1, C1 and C2. The new approach gave acceptable simulation results for the treatments under different water conditions, and especially gave better results for the treatments influenced markedly by precipitation or irrigation than the other approaches. The large errors in the late stage in 2009 might be due to the effect of the early maturity of winter wheat. It also indicated that the critical soil water storages in the equation in FAO24 were not suitable to describe the relationship between ETgw and soil water storage in root zone when water table was at or above 150?cm depth, and could be replaced with the storage corresponding to the field capacity and the steady storage. The new approach was more suitable for estimating the ETgw under different soil water conditions comparing with C1 and C2.  相似文献   

2.
Richards方程在非饱和渗流模拟及其他相关领域应用广泛。在数值求解过程中,可以采用有限差分方法进行数值离散并迭代求解,为了获得较可靠的数值解,常规的均匀网格空间步长往往是较小的。在一些不利数值条件下,如入渗于干燥土壤,迭代计算费时甚至精度也不能得到很好改善。因此,文章提出Chebyshev空间网格改进方法,结合有限差分方法对Richards方程进行数值离散以获得线性方程组,并通过经典的Picard迭代方法进行迭代求解线性方程组以得到Richards方程的数值解。通过均质土和分层土2个不利情况下的非饱和渗流算例,又结合模型解析解和软件Hydrus-1D,对比研究了改进网格方法与均匀网格方法获得数值解的精度。结果表明,提出的Chebyshev网格方法相较于传统的均匀网格,可以在较少的节点数下获得较高的数值精度,又具有较小的计算开销,有较好的应用前景。  相似文献   

3.
Irrigation is critical for oasis agriculture in Alar irrigated region, Northwestern China. The reference evapotranspiration (ET0) estimation using a modified Penman-Monteith method recommended by FAO is a baseline for rational irrigation scheduling. Thus, crop coefficients (Kc) were estimated using a ratio between reference and observed evapotranspiration. Then, we analyzed the relationship between temporal crop coefficients (Kc) and synchronized normalized difference vegetation index (NDVI) derived from multi-temporal remotely sensed images. The results showed that reference evapotranspiration increased from spring to summer hereafter decreased gradually to winter with the peak of ET0 ranging from 8 to 9 mm/day. A closely linear relation between NDVI and Kc of cotton was finally acquired with a coefficient of determination of 0.94. An acceptable accuracy of this model by comparing observed and simulated evapotranspiration illustrated the credibility of this model. This paper provides a scientific approach to estimate regional crop coefficient and evapotranspiration.  相似文献   

4.
Evapotranspiration is difficult to measure and, when measured, its spatial variability is not usually taken into account. The recommended method to estimate evapotranspiration, Penman-Monteith FAO, requires variables not available in most weather stations. Simplified but less accurate methods, as Hargreaves equation, are normally used. Several approaches have been proposed to improve Hargreaves equation accuracy. In this work, 14 calibrations of the Hargreaves equation are compared. Three goodness of fit statistics were used to select the optimal, in terms of simplicity and accuracy. The best option was an annual linear regression. Its parameters were interpolated using regression-kriging combining Random Forest and Ordinary Kriging. Twelve easy to obtain ancillary variables were used as predictors. The same approach was used to interpolate Hargreaves and Penman-Monteith-FAO ET0 on a daily basis; the Hargreaves ET0 layers and the parameter layers were used to obtain calibrated ET0 estimations. To compare the spatial patterns of the three estimations the daily layers were integrated into annual layers. The results of the proposed calibration are much more similar to Penman-Monteith FAO results than those obtained with Hargreaves equation. The research was conducted in south-east Spain with 79 weather stations with data from 01/01/2003 to 31/12/2014.  相似文献   

5.
林楠  姜然哲  刘强  郭晓东  杨航  陈思 《中国地质》2021,48(5):1392-1407
分析地表蒸散发时空变化规律及驱动因素,对促进区域水资源的科学分配、做好生态系统水源保护具有重要意义。本文基于MOD16蒸散发遥感数据产品,采用趋势分析及显著性检验法,深入分析了近20年三江平原地表蒸散量的时空变化特征,根据Penman-Monteith公式选取与地表蒸散量（ET）相关的驱动因子,分析各驱动因子对地表蒸散量变化的影响,并构建岭回归统计模型,分析研究地表蒸散量变化的主要驱动因子及其相对贡献率。结果表明：近20年三江平原地表蒸散发（ET）年际起伏特征明显,整体呈上升趋势;研究区内91.53%的地区ET呈增加趋势,且ET分布的地域差异逐年缩小;年内ET呈单峰型周期性变化,季节差异性明显;研究区坡度对ET有正向影响,高程和风速对ET有负向影响;气温、日照时数、降水量及NDVI与ET均呈正相关性,其中降水量与ET相关性最为显著;构建岭回归驱动分析模型的决定系数R2为0.823,能够有效解释各驱动因素与ET的关系。模型计算结果表明：降水量和植被覆盖度对三江平原地表蒸散量影响较大,是影响地表蒸散量变化的主要驱动力。  相似文献   

6.
In semi-arid climates, phreatophytes draw on shallow aquifers, and groundwater evapotranspiration (ETG) is a principal component of groundwater budgets. Diurnal water table fluctuations, which often are a product of ETG, were monitored in the riparian zone of Red Canyon Creek, Wyoming, USA. These fluctuations were higher in a riparian wetland (2–36 mm) than a grass-covered meadow (1–6 mm). The onset and cessation of water-table fluctuations correspond to daily temperatures relative to freezing. Spatial differences were due to vegetation type and specific yield, while temporal changes were due to vegetation dormancy. Ratios of ETG to potential evapotranspiration (PET), K c,GW, were similar to ratios of actual evapotranspiration (ET) to PET, K c, in semi-arid rangelands. Before vegetation senescence, K c,GW increased between precipitation events, suggesting phreatophytes pull more water from the saturated zone as soil moisture decreases. In contrast, K c decreases with soil moisture following precipitation events as ET becomes increasingly water-limited. Error in ETG is primarily from estimates of specific yield (S y), which is difficult to quantify in heterogeneous sediments. ETG values may be more reliable because the range of acceptable S y is smaller than K c and S y does not change with vegetation type or soil moisture.  相似文献   

7.
Unsaturated flow problems in porous media often described by Richards’ equation are of great importance in many engineering applications. In this contribution, we propose a new numerical flow approach based on isogeometric analysis (IGA) for modeling the unsaturated flow problems. The non-uniform rational B-spline (NURBS) basis is utilized for spatial discretization whereas the stable implicit backward Euler method for time discretization. The nonlinear Richards’ equation is iteratively solved with the aid of the Newton–Raphson scheme. Owing to some desirable features of an efficient numerical flow approach, major advantages of the present formulation involve: (a) numerical oscillation at the wetting front can be avoided or facilitated, simply by using either an h-refinement or a lumped mass matrix technique; (b) higher-order exactness can be obtained due to the nature of the IGA features; (c) the approach is straightforward to implement and it does not need any transformation, e.g., Kirchhoff transformation or filter algorithm; and (d) in contrast to the Picard iteration scheme, which forms linear convergences, the proposed approach can however yield quadratic convergences by using the Newton–Raphson method for solving resultant nonlinear equations. Numerical model validation is analyzed by solving a three-dimensional unsaturated flow problem in soil, and its derived results are verified against analytical solutions. Numerical applications are then studied by considering three extensive examples with simple and complex configurations to further show the accuracy and applicability of the present IGA.  相似文献   

8.
Modeling of soil?Cwater, ?Cheat and ?Ccarbon (C) fluxes provides an important tool for predicting mass and energy transfers based on a hydraulic-, thermal- and C-mass balance approach. Model predictions were evaluated using measured data from two water-limited study sites, one pasture and one supporting an alfalfa crop, to indentify differences between these ecosystems. Soil water content, temperature, and evapotranspiration (ET) data were used to validate soil water dynamics components of a process-based numerical model. Soil surface CO2 efflux estimates (i.e., fluxes from soil respiration) were also made to estimate soil CO2 emissions. The results show that the Hydrus-1D numerical model can be parameterized to simulate the soil hydrodynamics and CO2 fluxes measured at both locations. Rainfall and irrigation events triggering increases in plant root and microbial respiration rates were simulated to recreate observed pulsed CO2 fluxes. There were distinct differences in ET and soil CO2 effluxes between the ecosystems and watering events significantly modified the fluxes. Differences in potential evapotranspiration and soil texture could help explain these discrepancies. The results demonstrate that numerical modeling can be a useful tool for estimating soil surface fluxes in calibrated ecosystems when micrometeorological methods may not be suitable.  相似文献   

9.
Variations in reference evapotranspiration (ET0) and drought characteristics play a key role in the effect of climate change on water cycle and associated ecohydrological patterns. The accurate estimation of ET0 is still a challenge due to the lack of meteorological data and the heterogeneity of hydrological system. Although there is an increasing trend in extreme drought events with global climate change, the relationship between ET0 and aridity index in karst areas has been poorly studied. In this study, we used the Penman–Monteith method based on a long time series of meteorological data from 1951 to 2015 to calculate ET0 in a typical karst area, Guilin, Southwest China. The temporal variations in climate variables, ET0 and aridity index (AI) were analyzed with the Mann–Kendall trend test and linear regression to determine the climatic characteristics, associated controlling factors of ET0 variations, and further to estimate the relationship between ET0 and AI. We found that the mean, maximum and minimum temperatures had increased significantly during the 65-year study period, while sunshine duration, wind speed and relative humidity exhibited significant decreasing trends. The annual ET0 showed a significant decreasing trend at the rate of ?8.02 mm/10a. However, significant increase in air temperature should have contributed to the enhancement of ET0, indicating an “evaporation paradox”. In comparison, AI showed a slightly declining trend of ?0.0005/a during 1951–2015. The change in sunshine duration was the major factor causing the decrease in ET0, followed by wind speed. AI had a higher correlation with precipitation amount, indicating that the variations of AI was more dependent on precipitation, but not substantially dependent on the ET0. Although AI was not directly related to ET0, ET0 had a major contribution to seasonal AI changes. The seasonal variations of ET0 played a critical role in dryness/wetness changes to regulate water and energy supply, which can lead to seasonal droughts or water shortages in karst areas. Overall, these findings provide an important reference for the management of agricultural production and water resources, and have an important implication for drought in karst regions of China.  相似文献   

10.
Climate change affects the environment and natural resources immensely. Rainfall, temperature and evapotranspiration are major parameters of climate affecting changes in the environment. Evapotranspiration plays a key role in crop production and water balance of a region, one of the major parameters affected by climate change. The reference evapotranspiration or ET0 is a calculated parameter used in this research. In the present study, changes in the future rainfall, minimum and maximum temperature, and ET0 have been shown by downscaling the HadCM3 (Hadley Centre Coupled Model version 3) model data. The selected study area is located in a part of the Narmada river basin area in Madhya Pradesh in central India. The downscaled outputs of projected rainfall, ET0 and temperatures have been shown for the 21st century with the HADCM3 data of A2 scenario by the Least Square Support Vector Machine (LS-SVM) model. The efficiency of the LS-SVM model was measured by different statistical methods. The selected predictors show considerable correlation with the rainfall and temperature and the application of this model has been done in a basin area which is an agriculture based region and is sensitive to the change of rainfall and temperature. Results showed an increase in the future rainfall, temperatures and ET0. The temperature increase is projected in the high rise of minimum temperature in winter time and the highest increase in maximum temperature is projected in the pre-monsoon season or from March to May. Highest increase is projected in the 2080s in 2081–2091 and 2091–2099 in maximum temperature and 2091–2099 in minimum temperature in all the stations. Winter maximum temperature has been observed to have increased in the future. High rainfall is also observed with higher ET0 in some decades. Two peaks of the increase are observed in ET0 in the April–May and in the October. Variation in these parameters due to climate change might have an impact on the future water resource of the study area, which is mainly an agricultural based region, and will help in proper planning and management.  相似文献   

11.
One of the significant problems in geo-environmental and geotechnical engineering is the unsaturated flow of soil in unsaturated soils. The model of this phenomenon in porous media is governed by the Richards equation. In this paper a new, efficient, iterative method is used to handle the Richards equation. This new technique is obtained from the variational iteration method by a simple reconstruction that is the Laplace iteration method (LIM). In order to evaluate the efficiency and accuracy of the solutions obtained by the proposed method, two representative examples were investigated. The obtained results show that the Laplace iteration method is a very effective method, simplifies the difficulty of classical techniques and is quite accurate for systems of partial differential equations.  相似文献   

12.
朱帅润  吴礼舟 《岩土力学》2022,43(3):697-707
Richards方程常用于非饱和土渗流问题,并且应用广泛。在数值求解中,对Richards方程线性化,进而采用有限差分法进行数值离散以及迭代计算。其中传统的迭代法比如Jacobi迭代、Gauss-Seidel迭代法（GS）和连续超松驰迭代法（successive over-relaxation method,简称SOR）迭代收敛率较慢,尤其在离散空间步长较小以及离散时间步长较大时。因此,采用整体校正法以及多步预处理法对传统迭代法进行改进,提出一种基于整体校正法的多步预处理Gauss-Seidel迭代法（improved Gauss-Seidel iterative method with multistep preconditioner based on the integral correction method,简称ICMP(m)-GS）求解Richards方程导出的线性方程组。通过非饱和渗流算例,并与传统迭代法和解析解对比,对改进算法的收敛率和加速效果进行了验证。结果表明,提出的ICMP(m)-GS可以很大程度地改善线性方程组的病态性,相较于常规方法GS,SOR以及单一改进方法,ICMP(m)-GS具有更快的收敛率,更高的计算效率和计算精度。该方法可以为非饱和土渗流的数值模拟提供一定参考。  相似文献   

13.
The design of a top-cover system is a very important aspect of the closure of sanitary landfills. Soil material properties, sequencing, and vegetation are critical to top-cover design. To date, very few field investigations have been undertaken to determine the effectiveness of these designs. Results from a field investigation of top-covers proposed for an active landfill site in south Florida are reported here. Measurements were made of each component of the water balance equation: S = [P +I] – [ET +R +Q] whereS, P, I, ET, R, andQ represent water storage, rainfall, irrigation, evapotranspiration, surface runoff, and deep seepage, respectively. Deep seepage, determined using soil water tensiometers, is the amount of water infiltrating the landfill top-cover and is therefore the upper limit of potential leachate generation.Results indicate that the locally available rock tailings (carbonate silt) is an effective topliner component of the top-cover system. The results of the field investigation also showed that the EPA water balance model used to estimate infiltration during landfill design was greater than the actual field measurements.  相似文献   

14.
New Approach for Estimation of Static and Seismic Active Earth Pressure   总被引：1，自引：1，他引：1  
To estimate static and seismic active earth pressure (Pad) on a rigid retaining wall, numerical analyses using different step sizes have been carried out in this paper, based on the modified Culmann line method by considering Coulomb’s planar rupture surface. Equivalent pseudo-static seismic forces are considered in the analysis. A new concept of modified unit weight by considering ground surcharge is introduced under static and seismic conditions. By numerical analysis, area of soil (A) has been estimated to obtain the ratio of A/A0 where A0 is θh2, θ is the angle between retaining structure and ground surface and h is the vertical height of the wall. This ratio remains constant for a particular type of soil and has been used to estimate the maximum active earth pressure using force diagram. Results are provided in tabular form for easy calculation of the coefficient of static and seismic active earth pressure. Present results by considering the new technique, compares well with the results obtained by earlier researchers.  相似文献   

15.
The objective of the current study is to use satellite data to assess the mutual influence between vegetation and climate. The Ismailia Governorate was selected as a case study to investigate the impact of vegetation cover expansion on both land surface and air temperature from 1983 to 2010 and vice versa. This observation site was carefully selected as a clear example of the high rate of the reclamation and vegetation expansion process in Egypt. Land surface temperature (LST) was estimated through the Advanced Very High Resolution Radiometer (a space-borne sensor embarked on the National Oceanic and Atmospheric Administration) data while air temperature (T air) was collected from ground meteorological stations in the study area. Irrigated agriculture is the largest consumer of freshwater resources. However, consistent information on irrigation water use is still lacking. Relative humidity, wind speed, solar radiation, and T air data were inserted in the Penman–Monteith equation to calculate potential evapotranspiration (ETo), while both LST and T air were used to observe the relative water status of the study area as a result of the water deficit index (WDI). Then, both WDI and ETo were used to calculate actual evepotranspiration (ETC.). The results showed that LST decreased by about 2.3 °C while T air decreased by about 1.6 °C during the study period. The results showed also that the vegetation cover expanded from 25,529.85 ha in 1985 to 63,140.49 ha in 2009 with about 147 % increase. This decrease in LST and air temperature was according to the expansion of the cultivated land that was proved through the processing of three Landsat TM and Landsat ETM+ imageries acquired in June 19, 1985, June 7, 1998, and June 29, 2009. The vegetation water consumption was affected by the decreasing surface and air temperature. The results showed that the water deficit index decreased by about 0.35, and actual evapotranspiration increased by about 2.5 mm during the study period.  相似文献   

16.
Zhu  S. R.  Wu  L. Z.  Huang  J. 《Computational Geosciences》2022,26(1):131-145

This paper studies the potential of using the successive over-relaxation iteration method with polynomial preconditioner (P(m)-SOR) to solve variably saturated flow problems described by the linearized Richards’ equation. The finite difference method is employed to numerically discretize and produce a system of linear equations. Generally, the traditional Picard method needs to re-evaluate the iterative matrix in each iteration, so it is time-consuming. And under unfavorable conditions such as infiltration into extremely dry soil, the Picard method suffers from numerical non-convergence. For linear iterative methods, the traditional Gauss-Seidel iteration method (GS) has a slow convergence rate, and it is difficult to determine the optimum value of the relaxation factor w in the successive over-relaxation iteration method (SOR). Thus, the approximate optimum value of w is obtained based on the minimum spectral radius of the iterative matrix, and the P(m)-SOR method is extended to model underground water flow in unsaturated soils. The improved method is verified using three test examples. Compared with conventional Picard iteration, GS and SOR methods, numerical results demonstrate that the P(m)-SOR has faster convergence rate, less computation cost, and good error stability. Besides, the results reveal that the convergence rate of the P(m)-SOR method is positively correlated with the parameter m. This method can serve as a reference for numerical simulation of unsaturated flow.

  相似文献   

17.
基于分层假设的Green-Ampt模型改进   总被引：2，自引：0，他引：2       下载免费PDF全文
传统Green-Ampt入渗模型由于假设过于简化而限制了计算精度,针对这一问题,在分层假设的基础上对模型进行了改进:将均质入渗土壤分为饱和层、过渡层和干土层,采用Richards模型分析了不同入渗条件下各层的厚度以及各层内含水率和导水率的变化规律,结果表明,入渗条件一定时,过渡层占湿润层(包括饱和层和过渡层)比例会随湿润层厚度的增加线性减小,积水深度和土壤初始含水率分别会影响这一线性关系的纵轴截距和斜率,而土壤饱和导水率则对其没有明显影响;过渡层内土壤含水率分布为椭圆形曲线,而导水率随深度增加线性减小。在此基础上,通过引入湿润层的等效导水率并相应给出锋面吸力的计算方法,对传统Green-Ampt模型进行了改进。结果显示,改进后模型计算精度显著提高。  相似文献   

18.
This work concerns linearization methods for efficiently solving the Richards equation, a degenerate elliptic-parabolic equation which models flow in saturated/unsaturated porous media. The discretization of Richards’ equation is based on backward Euler in time and Galerkin finite elements in space. The most valuable linearization schemes for Richards’ equation, i.e. the Newton method, the Picard method, the Picard/Newton method and the L-scheme are presented and their performance is comparatively studied. The convergence, the computational time and the condition numbers for the underlying linear systems are recorded. The convergence of the L-scheme is theoretically proved and the convergence of the other methods is discussed. A new scheme is proposed, the L-scheme/Newton method which is more robust and quadratically convergent. The linearization methods are tested on illustrative numerical examples.  相似文献   

19.
This work presents the application of a Monte Carlo simulation method to perform an statistical analysis of transient variably saturated flow in an hypothetical random porous media. For each realization of the stochastic soil parameters entering as coefficients in Richards' flow equation, the pressure head and the flow field are computed using a mixed finite element procedure for the spatial discretization combined with a backward Euler and a modified Picard iteration in time. The hybridization of the mixed method provides a novel way for evaluating hydraulic conductivity on interelement boundaries. The proposed methodology can handle both large variability and fractal structure in the hydraulic parameters. The saturated conductivity K s and the shape parameter vg in the van Genuchten model are treated as stochastic fractal functions known as fractional Brownian motion (fBm) or fractional Gaussian noise (fGn). The statistical moments of the pressure head, water content, and flow components are obtained by averaging realizations of the fractal parameters in Monte Carlo fashion. A numerical example showing the application of the proposed methodology to characterize groundwater flow in highly heterogeneous soils is presented.  相似文献   

20.
非饱和土路基毛细作用的数值与解析方法研究   总被引：2，自引：0，他引：2  
刘杰  姚海林  卢正  胡梦玲  董启朋 《岩土力学》2013,34(Z2):421-427
通过建立非饱和土毛细作用的孔隙分布的分形模型,推导获得了非饱和土毛细水的最大上升高度,同时基于Richards水分运动微分方程,引入边界及初始条件,基于Laplace变换,得到了毛细作用下非饱和土路基湿度变化的解析解;同时引入算例,将所提出的解析方法计算结果与未简化参数的数值计算结果进行了对比分析;最后考虑不同因素的影响进行了非饱和土路基毛细作用下的湿度变化分析。分析结果表明：解析求解获得的路基湿度变化趋势和未进行参数简化的数值法求解结果基本一致,证明解析解法是合理可信的;路基填筑的初始含水率越大,填土的初始吸力越小,毛细水上升的高度及湿度变化量也相应越小;透水性能较好的路基填土毛细水上升速度较快,但上升高度较小,毛细水可以在较短时间内上升到最大高度;路基的填土类型不同,路基在毛细作用下的湿度变化状态也不同,需要针对不同的填土路基进行相应的防排水措施。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号