基于水平衡模型的呼伦湖湖泊水量变化

针对北方寒旱区呼伦湖水位下降、水面萎缩的现象,根据气候特征,利用月水量平衡模型探究湖泊水文过程并揭示其变化规律.在此基础上,利用不同气候条件下各水平衡项对于湖泊水位的影响程度确定水位升降的直接原因.基于1963-1980年间水位的实测数据,根据水量平衡原理及其他辅助计算判断出湖泊与周边区域存在着地下水的交换,且具有一定的规律性,即历年11月至次年3月期间的累积降雪融化渗入土壤中形成浅层径流补给湖泊,而7、8月份湖泊补给周边草原.基于以上规律,根据周边坡面汇流、地下水与湖泊交换量的年内变化特征,利用水平衡方程式推算湖泊1981-2008年逐月水位变化,并与其他研究成果比较,吻合度较高.不同气候条件下,径流量对于湖泊水位的影响程度最为突出,是水位变化的主控因子.     

From rainfed agriculture to stress-avoidance irrigation: I. A generalized irrigation scheme with stochastic soil moisture

With vast regions already experiencing water shortages, it is becoming imperative to manage sustainably the available water resources. As agriculture is by far the most important user of freshwater and the role of irrigation is projected to increase in face of climate change and increased food requirements, it is particularly important to develop simple, widely applicable models of irrigation water needs for short- and long-term water resource management. Such models should synthetically provide the key irrigation quantities (volumes, frequencies, etc.) for different irrigation schemes as a function of the main soil, crop, and climatic features, including rainfall unpredictability. Here we consider often-employed irrigation methods (e.g., surface and sprinkler irrigation systems, as well as modern micro-irrigation techniques) and describe them under a unified conceptual and theoretical framework, which includes rainfed agriculture and stress-avoidance irrigation as extreme cases. We obtain mostly analytical solutions for the stochastic steady state of soil moisture probability density function with random rainfall timing and amount, and compute water requirements as a function of climate, crop, and soil parameters. These results provide the necessary starting point for a full assessment of irrigation strategies, with reference to sustainability, productivity, and profitability, developed in a companion paper [Vico G, Porporato A. From rainfed agriculture to stress-avoidance irrigation: II. Sustainability, crop yield, and net profit. Adv Water Resour 2011;34(2):272-81].     

1-g Experimental investigation of bi-layer soil response and kinematic pile bending

The effect of soil inhomogeneity and material nonlinearity on kinematic soil–pile interaction and ensuing bending under the passage of vertically propagating seismic shear waves in layered soil, is investigated by means of 1-g shaking table tests and nonlinear numerical simulations. To this end, a suite of scale model tests on a group of five piles embedded in two-layers of sand in a laminar container at the shaking table facility in BLADE Laboratory at University of Bristol, are reported. Results from white noise and sine dwell tests were obtained and interpreted by means of one-dimensional lumped parameter models, suitable for inhomogeneous soil, encompassing material nonlinearity. A frequency range from 0.1 Hz to 100 Hz and 5 Hz to 35 Hz for white noise and sine dwell tests, respectively, and an input acceleration range from 0.015 g to 0.1 g, were employed. The paper elucidates that soil nonlinearity and inhomogeneity strongly affect both site response and kinematic pile bending, so that accurate nonlinear analyses are often necessary to predict the dynamic response of pile foundations.     

Dual-porosity and kinematic wave approaches to assess the degree of preferential flow in an unsaturated soil

Abstract

The purpose of this study was to assess the degree of preferential flow in an unsaturated soil column using two different models: the dual-porosity model, MACRO, and the kinematic wave approach (KWA) based on boundary-layer flow theory. The soil column experiments consisted of six infiltrations with intensities varying from 15 to 101 mm h^?1. Bromide solution was also infiltrated at an intensity of 79 mm h^?1 and a concentration of 80 mg l^?1. Both MACRO and the KWA indicated the absence of pure preferential flow. The KWA indicated intermediate flow with dispersion of the wetting front with depth, whereas MACRO indicated flow dominated by the diffusion of capillary potential. These results shed light on the transition between flows dominated by momentum dissipation and by diffusion of capillary potential. The absence of pure macropore flow in the structured sandy soil is mainly due to efficient lateral mass exchange in this material.     

Hillslope‐scale probabilistic characterization of soil moisture dynamics and average water balance

Probabilistic water balance modelling provides a useful framework for investigating the interactions between soil, vegetation, and the atmosphere. It has been used to estimate temporal soil moisture dynamics and ecohydrological responses at a point. This study combines a nonlinear rainfall–runoff theory with probabilistic water balance model to represent varied source area runoff as a function of rainfall depth and a runoff coefficient at hillslope scale. Analytical solutions of the soil‐moisture probability density function and average water balance model are then developed. Based on a sensitivity analysis of soil moisture dynamics, we show that when varied source area runoff is incorporated, mean soil moisture is always lower and total runoff higher, compared with the original probabilistic water balance model. The increased runoff from the inclusion of varied source area runoff is mainly because of a reduction in leakage when the index of dryness is less than one and evapotranspiration when the index of dryness is greater than one. Inclusion of varied source area runoff in the model means that the actual evapotranspiration is limited by less available water (i.e. water limit), which is stricter than Budyko’s and Milly’s water limit. Application of the model to a catchment located in Western Australia showed that the method can predict annual value of actual evapotranspiration and streamflow accurately. Copyright © 2012 John Wiley & Sons, Ltd.     

Validation of a soil water balance model using soil water content and pressure head data

The validation of soil water balance models and the evaluation of the quality of the model predictions at field‐scale require time‐series of in situ measured model outputs. In our study, we have validated such a model using a 6‐year period with time‐series of automatically recorded, daily volumetric soil water contents measured with the time‐domain reflectometry with intelligent microelements (TRIME) method and daily pressure heads measured with tensiometers. The comparisons of simulated with measured soil water contents and pressure heads were analysed using the modelling efficiency index (IA) and the square root of the mean square error (RMSE) in order to evaluate the prediction quality of the model. In our study, IA and RMSE, obtained either from the comparison of simulated with measured soil water contents or the comparison of calculated with observed pressure heads, in some cases lead to different results regarding the evaluation of the simulation quality of the soil water balance model. For example, a good fit between simulated and observed soil water contents does not necessarily result in a comparably good fit between the corresponding calculated and measured pressure heads. Therefore, a combined use of both measurement techniques, which takes into account their respective advantages and disadvantages, gives a more complete overview on the simulation quality of the soil water balance model than the single use of one of those techniques. Copyright © 2004 John Wiley & Sons, Ltd.     

Coupling the Xinanjiang model to a kinematic flow model based on digital drainage networks for flood forecasting

The Xinanjiang model, which is a conceptual rainfall‐runoff model and has been successfully and widely applied in humid and semi‐humid regions in China, is coupled by the physically based kinematic wave method based on a digital drainage network. The kinematic wave Xinanjiang model (KWXAJ) uses topography and land use data to simulate runoff and overland flow routing. For the modelling, the catchment is subdivided into numerous hillslopes and consists of a raster grid of flow vectors that define the water flow directions. The Xinanjiang model simulates the runoff yield in each grid cell, and the kinematic wave approach is then applied to a ranked raster network. The grid‐based rainfall‐runoff model was applied to simulate basin‐scale water discharge from an 805‐km² catchment of the Huaihe River, China. Rainfall and discharge records were available for the years 1984, 1985, 1987, 1998 and 1999. Eight flood events were used to calibrate the model's parameters and three other flood events were used to validate the grid‐based rainfall‐runoff model. A Manning's roughness via a linear flood depth relationship was suggested in this paper for improving flood forecasting. The calibration and validation results show that this model works well. A sensitivity analysis was further performed to evaluate the variation of topography (hillslopes) and land use parameters on catchment discharge. Copyright © 2009 John Wiley & Sons, Ltd.     

The characteristics of soil water cycle and water balance on steep grassland under natural and simulated rainfall conditions in the Loess Plateau of China

Large-scale vegetation restoration has been helpful to prevent serious soil erosion, but also has aggravated water scarcity and resulted in soil desiccation below a depth of 200 cm in the Loess Plateau of China. To understand the dynamic mechanism of soil desiccation formation is very important for sustainable development of agriculture in the Loess Plateau. Based on natural and simulated rainfall, the characteristics of soil water cycle and water balance in the 0–400 cm soil layer on a steep grassland hillslope in Changwu County of Shaanxi Loess Plateau were investigated from June to November in 2002, a drought year with annual rainfall of 460 mm. It was similarly considered to represent a rainy year with annual rainfall of 850 mm under simulated rainfall conditions. The results showed that the temporal variability of water contents would decrease in the upper 0–200 cm soil layer with the increase in rainfall. The depth of soil affected by rainfall infiltration was 0–200 cm in the drought year and 0–300 cm in the rainy year. During the period of water consumption under natural conditions, the deepest layer of soil influenced by evapotranspiration (ET) rapidly reached a depth of 200 cm on July 21, 2002, and soil water storage decreased by 48 mm from the whole 0–200 cm soil layer. However, during the same investigation period under simulated rainfall conditions, soil water storage in the 0–400 cm soil layer increased by only 71 mm, although the corresponding rainfall was about 640 mm. The extra-simulated rainfall of 458 mm from May 29 to August 10 did not result in the disappearance of soil desiccation in the 200–400 cm deep soil layer. Most infiltrated rainwater retained in the top 0–200 cm soil layer, and it was subsequently depleted by ET in the rainy season. Because very little water moved below the 200 cm depth, there was desiccation in the deep soil layer in drought and normal rainfall years.     

Coupling snow accumulation and melt rate modules of monthly water balance models with the Jazim monthly water balance model

Several models for simulation of water balance processes in semi-arid mountainous basins were developed by coupling different modules of existing water balance models (WBM). Snow accumulation and snowmelt rate relationships extracted from the McCabe-Markstrom, Guo, Rao-Al Wagdany and WASMOD-M WBMs, originally developed for basins with humid climate, were coupled with the Jazim WBM, primarily developed for arid basins. Karaj Basin, central Iran, with snowy autumn–winter and dry summer periods, was selected to assess model performance. The model parameters were optimized using a genetic algorithm (GA). All coupled models performed better than the non-modified (original) WBMs in the study basin. The coupled Jazim–McCabe-Markstrom model provided the best performance in simulating low and high monthly flows. It estimated the snowmelt runoff values more accurately than other proposed coupled models because the linear relationships used in the snow module of the McCabe-Markstrom model are more compatible with snow variations in the Karaj Basin.     

A shot noise model for the generation of soil moisture data

An attempt is made to estimate the expected contribution of rainfall to soil moisture during the irrigation season. Effective rainfall and evapotranspiration are the parameters considered in the water balance carried out in the root zone. Rainfall occurrence is simulated by a Poisson process whereas evapotranspiration is described by a simple deterministic function of potential evapotranspiration and soil moisture in the root zone. Using the theory of shot noise models a closed form solution is derived from the expected soil moisture in the root zone at the end of the time interval (0,t]. For illustration purposes the proposed model is applied to a series of data from Mikra meteorological station in Greece.List of symbols x change in water storage in the root zone during the time interval t - X infiltrated rainfall of thei th storm event - ET actual evapotranspiration during thej th day - Poisson rate - number of storm events in (0,t] - t _r duration of rainfall - t _b interarrival time - h _i rainfall depth of thei th storm event - i _m mean rainfall intensity - i(t) instantaneous rainfall intensity - x(0),x(t) available soil moisture in the root zone at time 0 andt, respectively - PET potential evapotranspiration rate - x _F available soil moisture in the root zone at field capacity - soil moisture depletion rate (=PET/x _F ) - w impulse shape of filtered Poisson processes - E[·] mean value - S _i time of thei th rainfall event - N(t) time of storm events in (0,t] - estimated standard deviation The following symbols were used in this paper     

Wetlands with controlled drainage and sub‐irrigation systems—modelling of the water balance

Over the past centuries, the agricultural use of wetlands in Central Europe has required interference with the natural wetland water balance. Often this has consisted of drainage measures alone. In low‐precipitation areas, it has also involved the operation of combined drainage and sub‐irrigation systems. Model studies conducted as part of planning processes, or with a view to finding out the impact of changing climate conditions on the water balance of wetlands, must take these facts into account. For this reason, a water balance model has been devised for wetlands whose water balance is governed by water resources management systems. It is based on the WBalMo model system. Special modules were integrated into WBalMo to calculate the water balance of wetland areas (WABI module) and to regulate inflow partitioning within the wetland (REGINF module). When calculating the water balance, the WABI module takes into account precipitation and potential evapotranspiration, groundwater levels below surface, soil types, land‐use classes, inflows via the running water system, and data for target water levels. It provides actual evapotranspiration, discharge into the running water system, and groundwater levels in the area. The example of the Spreewald, a major wetland area in north‐eastern Germany, was used to design and test the WBalMo Spreewald model. The comparison of measured and calculated water balance parameters of the wetland area confirms the suitability of the model for water balance studies in wetlands with complex water resources management systems. The results reveal the strong influence of water management on the water balance of such areas. The model system has proved to be excellently suited for planning and carrying out water management measures aimed at the sustainable development of wetlands. Furthermore, scenario analyses can be used to assess the impact of global change on the water balance of wetlands. Copyright © 2006 John Wiley & Sons, Ltd.     

维波场延拓复杂表层模型校正方法及应用

如何正确地消除复杂地表对地震波场的影响，提高地下构造成像的质量一直是中国西部复杂地区地震勘探中存在的难题.本文在三维复杂表层速度模型层析反演\[1\]的基础上，研究了关于复杂地表的静校正问题，提出用三维波动方程在炮集上对地震波场进行正、反向延拓，消除复杂地表对波场的影响，实现三维复杂表层模型校正.理论和实际应用证明，该方法已超越常规静校正的含义，属时变校正方法.用本方法处理复杂地表问题，不但能消除表层对不同深度反射波产生的不同时差影响，提高叠加剖面质量，而且能使校正后的地震波场保持波动特征不发生畸变，可为建立正确的深层速度模型和波动方程叠前深度偏移奠定良好的基础.     

Development and application of a physically based landscape water balance in the SWAT model

Watershed scale hydrological and biogeochemical models rely on the correct spatial‐temporal prediction of processes governing water and contaminant movement. The Soil and Water Assessment Tool (SWAT) model, one of the most commonly used watershed scale models, uses the popular curve number (CN) method to determine the respective amounts of infiltration and surface runoff. Although appropriate for flood forecasting in temperate climates, the CN method has been shown to be less than ideal in many situations (e.g. monsoonal climates and areas dominated by variable source area hydrology). The CN model is based on the assumption that there is a unique relationship between the average moisture content and the CN for all hydrologic response units (HRUs), and that the moisture content distribution is similar for each runoff event, which is not the case in many regions. Presented here is a physically based water balance that was coded in the SWAT model to replace the CN method of runoff generation. To compare this new water balance SWAT (SWAT‐WB) to the original CN‐based SWAT (SWAT‐CN), two watersheds were initialized; one in the headwaters of the Blue Nile in Ethiopia and one in the Catskill Mountains of New York. In the Ethiopian watershed, streamflow predictions were better using SWAT‐WB than SWAT‐CN [Nash–Sutcliffe efficiencies (NSE) of 0·79 and 0·67, respectively]. In the temperate Catskills, SWAT‐WB and SWAT‐CN predictions were approximately equivalent (NSE > 0·70). The spatial distribution of runoff‐generating areas differed greatly between the two models, with SWAT‐WB reflecting the topographical controls imposed on the model. Results show that a water balance provides results equal to or better than the CN, but with a more physically based approach. Copyright © 2010 John Wiley & Sons, Ltd.     

基于水热平衡模型的青海湖水位变化趋势预测

近几十年来,随着气候干暖化,以青海湖为代表的我国内陆湖泊水位持续下降,生态环境问题日益突出,备受世人关注.运用改进的水热平衡模型预测了2050年以前青海湖逐年的湖面蒸发量,并运用多元线性回归的方法估算出流域未来径流量的变化,最终通过水量平衡的方式对2050年以前青海湖水位的变化趋势进行了定量预测.预测表明未来几十年内,青海湖水位会经历先相对稳定再继续下降的过程,2020年以前青海湖水位会相对稳定在3192.7m,之后会继续下降,到2050年约下降到3191.22m,总体上2010-2050年青海湖水位下降趋势将有所缓和.     

Comparative performance of soil water balance models in computing semi-arid aquifer recharge

Abstract

Estimating groundwater recharge is essential to ensure the sustainable use of groundwater resources, particularly in arid and semi-arid regions. Soil water balances have been frequently advocated as valuable tools to estimate groundwater recharge. This article compares the performance of three soil water balance models (Hydrobal, Visual Balan v2.0 and Thornthwaite) in the Ventós-Castellar aquifer, Spain. The models were used to simulate wet and dry years. Recharge estimates were transformed into water table fluctuations by means of a lumped groundwater model. These, in turn, were calibrated against piezometric data. Overall, the Hydrobal model shows the best fit between observed and calculated levels (r² = 0.84), highlighting the role of soil moisture and vegetation in recharge processes.

Editor D. Koutsoyiannis; Associate editor X. Chen

Citation Touhami, I., et al., 2014. Comparative performance of soil water balance models in computing semi-arid aquifer recharge. Hydrological Sciences Journal, 59 (1), 193–203.     

The impact of land use change on water balance in Zhangye city,China

Land use change has a significant effect on water balance, especially in arid region, such as Northwest China. In this paper, we analyze the effect of land use change on water balance in terms of the amount of water supply and demand from economic perspective. It's the first time to extend the basic 48 sectors input-output table to include water and land accounts that involved into multiple production processes for Zhangye city. We then perform the improved ORANI-G model, a single region Computable General Equilibrium (CGE) model, to analyze the effect of land use change on water balance under three scenarios. Subsequently, scenario-based simulation results are interpreted through selected sectors from agricultural, industrial, and service sectors respectively. Finally, the effect of land use change on water balance is analyzed through the difference between business-as-usual and land use unchanged scenarios. The results show that the extent of effect on water balance is different among sectors. Specifically, from the perspective of absolute value, service sectors are the largest, followed by industrial sectors, and the agricultural sectors are the least. Conversely, in terms of percentage change of land use, the largest extent of effect occurs in agricultural sectors. Additionally, with the rapid urbanization and the development of social economy, water balance in industrial sectors and service sectors will be stricken and reconstructed to a new high level. Simulation results also show that agricultural land shrinking will mitigate water scarcity distinctly, which indicates that balance the relationship among different stakeholders is imperative to guarantee water transformation from agricultural sectors to industrial and service sectors.     

Evaluation of the soil water balance in an alluvial flood plain with a shallow groundwater table

Abstract

Many of the hydrological and ecological functions of alluvial flood plains within watersheds depend on the water flow exchanges between the vadoze soil zone and the shallow groundwater. The water balance of the soil in the flood plain is investigated, in order to evaluate the main hydrological processes that underlie the temporal dynamics of soil moisture and groundwater levels. The soil moisture and the groundwater level in the flood plain were monitored continuously for a three-year period. These data were integrated with the results derived from applying a physically-based numerical model which simulated the variably-saturated vertical water flow in the soil. The analysis indicated that the simultaneous processes of lateral groundwater flow and the vertical recharge from the unsaturated zone caused the observed water table fluctuations. The importance of these flows in determining the rises in the water table varied, depending on soil moisture and groundwater depth before precipitation. The monitoring period included two hydrological years (September 2009–September 2011). About 13% of the precipitation vertically recharged the groundwater in the first year and about 50% in the second. The difference in the two recharge coefficients was in part due to the lower groundwater levels in the recharge season of the first hydrological year, compared to those observed in the second. In the latter year, the shallow groundwater increased the soil moisture in the unsaturated zone due to capillary rise, and so the mean hydraulic conductivity of the unsaturated soil was high. This moisture state of soil favoured a more efficient conversion of infiltrated precipitation into vertical groundwater recharge. The results show that groundwater dynamics in the flood plain are an important source of temporal variability in soil moisture and vertical recharge processes, and this variability must be properly taken into account when the water balance is investigated in shallow groundwater environments.

Citation Pirastru, M. and Niedda, M., 2013. Evaluation of the soil water balance in an alluvial flood plain with a shallow groundwater table. Hydrological Sciences Journal, 58 (4), 898–911.     

A water balance model to estimate climate change impact on groundwater recharge in Yucatan Peninsula,Mexico

ABSTRACT

The aim of this paper is to estimate the effect that climate change will have on groundwater recharge at the Yucatan Peninsula, Mexico. The groundwater recharge is calculated from a monthly water balance model considering eight methods of potential and actual evapotranspiration. Historical data from 1961–2000 and climate model outputs from five downscaled general circulation models in the near horizon (2015–2039), with representative concentration pathway (RCP) 4.5 and 8.5 are used. The results estimate a recharge of 118 ± 33 mm·year^–1 (around 10% of precipitation) in the historical period. Considering the uncertainty from GCMs under different RCP and evapotranspiration scenarios, our monthly water balance model estimates a groundwater recharge of 92 ± 40 mm·year^–1 (RCP4.5) and 94 ± 38 mm·year^–1 (RCP8.5) which represent a reduction of 23% and 20%, respectively, a result that threatens the socio-ecological balance of the region.     

The effects of initial soil moisture conditions on swale flow hydrographs

The effects of soil water content (SWC) on the formation of run‐off in grass swales draining into a storm sewer system were studied in two 30‐m test swales with trapezoidal cross sections. Swale 1 was built in a loamy fine‐sand soil, on a slope of 1.5%, and Swale 2 was built in a sandy loam soil, on a slope of 0.7%. In experimental runs, the swales were irrigated with 2 flow rates reproducing run‐off from block rainfalls with intensities approximately corresponding to 2‐month and 3‐year events. Run‐off experiments were conducted for initial SWC (SWC_ini) ranging from 0.18 to 0.43 m³/m³. For low SWC_ini, the run‐off volume was greatly reduced by up to 82%, but at high SWC_ini, the volume reduction was as low as 15%. The relative swale flow volume reductions decreased with increasing SWC_ini and, for the conditions studied, indicated a transition of the dominating swale functions from run‐off dissipation to conveyance. Run‐off flow peaks were reduced proportionally to the flow volume reductions, in the range from 4% to 55%. The swale outflow hydrograph lag times varied from 5 to 15 min, with the high values corresponding to low SWC_ini. Analysis of swale inflow/outflow hydrographs for high SWC_ini allowed estimations of the saturated hydraulic conductivities as 3.27 and 4.84 cm/hr in Swales 1 and 2, respectively. Such estimates differed from averages (N = 9) of double‐ring infiltrometer measurements (9.41 and 1.78 cm/hr). Irregularities in swale bottom slopes created bottom surface depression storage of 0.35 and 0.61 m³ for Swales 1 and 2, respectively, and functioned similarly as check berms contributing to run‐off attenuation. The experimental findings offer implications for drainage swale planning and design: (a) SWC_ini strongly affect swale functioning in run‐off dissipation and conveyance during the early phase of run‐off, which is particularly important for design storms and their antecedent moisture conditions, and (b) concerning the longevity of swale operation, Swale 1 remains fully functional even after almost 60 years of operation, as judged from its attractive appearance, good infiltration rates (3.27 cm/hr), and high flow capacity.