首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 703 毫秒
1.
Stone covers on loessial slopes can increase the time of infiltration by slowing the velocity of the overland flow, which reduces the transport of solutes, but few mechanistic models have been tested under water‐scouring conditions. We carried out field experiments to test a previously proposed, physically based model of water and solute transport. The area of soil infiltration was calculated from the uncovered surface area, and Richards' equation and the kinematic wave equation were used to describe water infiltration and flow along slopes with stone covers. The transport of chemicals into the run‐off from the surface soil, presumably by diffusion, and their movement in the soil profile could be described by the convection–diffusion equations of the model. The simulated and measured data correlated well. The stones on the soil surface reduced the area available for infiltration but increased the Manning coefficient, eventually leading to increased water infiltration and decreased solute loss with run‐off. Our results indicated that the traditional model of water movement and solute migration could be used to simulate water transport and solute migration for stone‐covered soil on loessial slopes.  相似文献   

2.
Jos C. van Dam 《水文研究》2000,14(6):1101-1117
Single domain models may seriously underestimate leaching of nutrients and pesticides to groundwater in clay soils with shrinkage cracks. Various two‐domain models have been developed, either empirical or physically based, which take into account the effects of cracks on water flow and solute transport. This paper presents a model concept that uses the clay shrinkage characteristics to derive crack volume and crack depth under transient field conditions. The concept has been developed to simulate field average behaviour of a field with cracks, rather than flow and transport at a small plot. Water flow and solute transport are described with basic physics, which allow process and scenario analysis. The model concept is part of the more general agrohydrological model SWAP, and is applied to a field experiment on a cracked clay soil, at which water flow and bromide transport were measured during 572 days. A single domain model was not able to mimic the field‐average water flow and solute transport. Incorporation of the crack concept considerably improved the simulation of water content and bromide leaching to the groundwater. Still deviations existed between the measured and simulated bromide concentration profiles. The model did not reproduce the observed bromide retardation in the top layer and the high bromide dispersion resulting from water infiltration at various soil depths. A sensitivity analysis showed that the amounts of bromide leached were especially sensitive to the saturated hydraulic conductivity of the top layer, the solute transfer from the soil matrix to crack water flow and the mean residence time of rapid drainage. The shrinkage characteristic and the soil hydraulic properties of the clay matrix showed a low sensitivity. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

3.
Solute plume subjected to field scale hydraulic conductivity heterogeneity shows a large dispersion/macrodispersion, which is the manifestation of existing fields scale heterogeneity on the solute plume. On the other hand, due to the scarcity of hydraulic conductivity measurements at field scale, hydraulic conductivity heterogeneity can only be defined statistically, which makes the hydraulic conductivity a random variable/function. Random hydraulic conductivity as a parameter in flow equation makes the pore flow velocity also random and the ground water solute transport equation is a stochastic differential equation now. In this study, the ensemble average of stochastic ground water solute transport equation is taken by the cumulant expansion method in order to upscale the laboratory scale transport equation to field scale by assuming pore flow velocity is a non stationary, non divergence-free and unsteady random function of space and time. Besides the stochastic explanation of macrodispersion and the velocity correction term obtained by Kavvas and Karakas (J Hydrol 179:321–351, 1996) before a new velocity correction term, which is a function of mean pore flow velocity divergence, is obtained in this study due to strict second order cumulant expansion (without omitting any term after the expansion) performed. The significance of the new velocity correction term is investigated on a one dimensional transport problem driven by a density dependent flow field.  相似文献   

4.
We studied the problem of local‐ and field‐scale infiltration over a particular class of heterogeneous soils. At the local scale, the soils are described as being vertically non‐uniform, with the saturated hydraulic conductivity continuously decreasing with depth according to a power law function. Analogous to the Green–Ampt model, analytical expressions are first developed for local‐scale infiltration using a sharp front approximation, and model results are compared with numerical solutions of the Richards equation. These results show that saturation does not occur from below in soils with such vertical non‐uniformity, thereby allowing for the use of a sharp front approximation. Because of vertical non‐uniformity, ponding conditions are achieved locally even for rainfall rates less than the surface saturated hydraulic conductivity. Furthermore, infiltration rates asymptotically approach zero at long times. To determine field‐scale infiltration properties, the spatial variability in the surface saturated hydraulic conductivity is represented by a log‐normal random field. Using cumulative infiltration as the independent variable, expressions are developed for the ensemble mean of field‐scale infiltration and the expected time for a given depth of water to infiltrate over the field. Surface horizontal heterogeneity is found to control field‐scale infiltration at small times, whereas local vertical non‐uniformity exerts a strong control at long times. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

5.
Water and solute movement in a coarse-textured water-repellent field soil   总被引：6，自引：0，他引：6  
Unstable water flow in water-repellent unsaturated soils can significantly affect the processes of infiltration and soil water redistribution. A field experiment was carried out to study the effect of water-repellency on water and bromide movement in a coarse-textured soil in the southwestern part of The Netherlands. The field data were analyzed using a relatively simple numerical model based on the standard Richards' equation for unsaturated water flow and the Fickian-based convection-dispersion equation for solute transport. Water-repellency was accounted for by multiplying the water content and the unsaturated hydraulic conductivity of the soil with F, a factor equal to the volumetric fraction of soil occupied by preferential flow paths resulting from the unstable flow process. The good comparison of simulated and measured bromide concentrations suggests that the model provides a viable method for simulating unstable water flow in water-repellent soils.  相似文献   

6.
Soil depth and soil production are highly complicated phenomena, generated from a complex interaction of physical, biological and chemical processes. It has, nevertheless, become increasingly clear that soil formation rates are closely related to chemical weathering rates. Somewhat paradoxically, it is likewise becoming apparent that such biogeochemical reactions as slowly transform rock to soil are limited by physical processes, such as flowing water and the formation of fractures. We have formulated a theoretical approach that relates soil formation rates to chemical weathering rates, and those, likewise, to solute transport rates. For such a theoretical framework to be relevant, the solute transport rates cannot equal those of the flowing water, as is the case in Gaussian solute transport. Rather, solute transport must be slowed in accordance with heavy‐tailed solute arrival time distributions. The inference is that the traditional advection–dispersion equation formulation for solute transport is inadequate in the typically heterogeneous geological media that weather to form soils. Here we examine the implications of this soil production model on the assumption of the approach to steady state. Particularly at slow erosion rates we find that many soil columns are not in equilibrium. This tendency may be accentuated in dry climates. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

7.
Macro-pores such as crab burrows are found commonly distributed in salt marsh sediments. Their disturbance on the soil structure is likely to influence both pore water flows and solute transport in salt marshes; however, the effects of crab burrows are not well understood. Here, a three-dimensional model simulated tidally driven pore water flows subject to the influence of crab burrows in a marsh system. The model, based on Richards’ equation, considered variably saturated flow in the marsh with a two-layer soil configuration, as observed at the Chongming Dongtan wetland (Shanghai, China). The simulation results showed that crab burrows distributed in the upper low-permeability soil layer, acting as preferential flow paths, affected pore water flows in the marsh particularly when the contrast of hydraulic conductivity between the lower high-permeability soil layer and the overlying low-permeability soils was high. The burrows were found to increase the volume of tidally driven water exchange between the marsh soil and the tidal creek. The simulations also showed improvement of soil aeration conditions in the presence of crab burrows. These effects may lead to increased productivity of the marsh ecosystem and enhancement of its material exchange with coastal waters.  相似文献   

8.
A one‐dimensional, two‐layer solute transport model is developed to simulate chemical transport process in an initially unsaturated soil with ponding water on the soil surface before surface runoff starts. The developed mathematical model is tested against a laboratory experiment. The infiltration and diffusion processes are mathematically lumped together and described by incomplete mixing parameters. Based on mass conservation and water balance equations, the model is developed to describe solute transport in a two‐zone layer, a ponding runoff zone and a soil mixing zone. The two‐zone layer is treated as one system to avoid describing the complicated chemical transport processes near the soil surface in the mixing zone. The proposed model was analytically solved, and the solutions agreed well with the experimental data. The developed experimental method and mathematical model were used to study the effect of the soil initial moisture saturation on chemical concentration in surface runoff. The study results indicated that, when the soil was initially saturated, chemical concentration in surface runoff was significantly (two orders of magnitude) higher than that with initially unsaturated soil, while the initial chemical concentrations at the two cases were of the same magnitude. The soil mixing depth for the initially unsaturated soil was much larger than that for the initially saturated soil, and the incomplete runoff mixing parameter was larger for the initially unsaturated soil. The higher the infiltration rate of the soil, the greater the infiltration‐related incomplete mixing parameter. According to the quantitative analysis, the soil mixing depth was found to be sensitive for both initially unsaturated and saturated soils, and the incomplete runoff mixing parameter was sensitive for initially saturated soil but not for the initially unsaturated soil; the incomplete infiltration mixing parameter behaved just the opposite. Some suggestions are made for reducing chemical loss from runoff. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

9.
Spreading‐basin methods have resulted in more than 130 million cubic metres of recharge to the unconfined Navajo Sandstone of southern Utah in the past decade, but infiltration rates have slowed in recent years because of reduced hydraulic gradients and (or) clogging. Trench infiltration is a promising alternative technique for increasing recharge and minimizing evaporation. This paper uses a variably saturated flow model to further investigate the relative importance of the following variables on rates of trench infiltration to unconfined aquifers: saturated hydraulic conductivity, trench spacing and dimensions, initial water‐table depth, alternate wet/dry periods, and number of parallel trenches. Modelling results showed (1) increased infiltration with higher hydraulic conductivity, deeper initial water tables, and larger spacing between parallel trenches, (2) deeper or wider trenches do not substantially increase infiltration, (3) alternating wet/dry periods result in less overall infiltration than keeping the trenches continuously full, and (4) larger numbers of parallel trenches within a fixed area increases infiltration but with a diminishing effect as trench spacing becomes tighter. An empirical equation for estimating expected trench infiltration rates as a function of hydraulic conductivity and initial water‐table depth was derived and can be used for evaluating feasibility of trench infiltration in other hydrogeologic settings. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

10.
Within the framework of stochastic theory and the spectral perturbation techniques, three-dimensional dispersion in partially saturated soils with fractal log hydraulic conductivity distribution is analyzed. Our analysis is focused on the impact of fractal dimension of log hydraulic conductivity distribution, local dispersivity, and unsaturated flow parameters, such as the soil poresize distribution parameter and the moisture distribution parameter, on the spreading behavior of solute plume and the concentration variance. Approximate analytical solutions to the stochastic partial differential equations are derived for the variance of asymptotic solute concentration and asymptotic macrodispersivities.  相似文献   

11.
Field determined hydraulic and chemical transport properties can be useful for the protection of groundwater resources from land-applied chemicals. Most field methods to determine flow and transport parameters are either time or energy consuming and/or they provide a single measurement for a given time period. In this study, we present a dripper-TDR field method that allows measurement of hydraulic conductivity and chemical transport parameters at multiple field locations within a short time period. Specifically, the dripper-TDR determines saturated hydraulic conductivity (Ks), macroscopic capillary length (λc), immobile water fraction (θim/θ), mass exchange coefficient (α) and dispersion coefficient (Dm). Multiple dripper lines were positioned over five crop rows in a field. Background and step solutions were applied through drippers to determine surface hydraulic conductivity parameters at 44 locations and surface transport properties at 38 locations. The hydraulic conductivity parameters (Ks, λc) were determined by application of three discharge rates from the drippers and measurements of the resultant steady-state flux densities at the soil surface beneath each dripper. Time domain reflectometry (TDR) was used to measure the bulk electrical conductivity of the soil during steady infiltration of a salt solution. Breakthrough curves (BTCs) for all sites were determined from the TDR measurements. The Ks and λc values were found to be lognormally distributed with average values of 31.4 cm h−1 and 6.0 cm, respectively. BTC analysis produced chemical properties, θim/θ, α, and Dm with average values of 0.23, 0.0036 h−1, and 1220 cm2 h−1, respectively. The estimated values of the flow and transport parameters were found to be within the ranges of values reported by previous studies conducted at nearby field locations. The dripper TDR method is a rapid and useful technique for in situ measurements of hydraulic conductivity and solute transport properties. The measurements reported in this study give clear evidence to the occurrence of non-equilibrium water and chemical movement in surface soil. The method allows for quantification of non-equilibrium model parameters and preferential flow. Quantifying the parameters is a necessary step toward determining the influences of surface properties on infiltration, runoff, and vadose zone transport.  相似文献   

12.
Understanding the dynamics and mechanisms of soil water movement and solute transport is essential for accurately estimating recharge rates and evaluating the impacts of agricultural activities on groundwater resources. In a thick vadose zone (0–15 m) under irrigated cropland in the piedmont region of the North China Plain, soil water content, matric potential, and solute concentrations were measured. Based on these data, the dynamics of soil water and solutes were analysed to investigate the mechanisms of soil water and solute transport. The study showed that the 0–15‐m vadose zone can be divided into three layers: an infiltration and evaporation layer (0–2 m), an unsteady infiltration layer (2–6 m), and a quasi‐steady infiltration layer (6–15 m). The chloride, nitrate, and sulphate concentrations all showed greater variations in the upper soil layer (0–1 m) compared to values in the deep vadose zone (below 2 m). The average concentrations of these three anions in the deep vadose zone varied insignificantly with depth and approached values of 125, 242, and 116 mg/L. The accumulated chloride, sulphate, and nitrate were 2,179 ± 113, 1,760 ± 383, and 4,074 ± 421 kg/ha, respectively. The soil water potential and solute concentrations indicated that uniform flow and preferential flow both occurred in the deep vadose zone, and uniform flow was the dominant mechanism of soil water movement in this study. The piston‐like flow velocity of solute transport was 1.14 m per year, and the average value of calculated leached nitrate nitrogen was 107 kg/ha?year below the root zone. The results can be used to better understand recharge processes and improve groundwater resources management.  相似文献   

13.
Modeling unsaturated flow in porous media requires constitutive relations that describe the soil water retention and soil hydraulic conductivity as a function of either potential or water content. Often, the hydraulic parameters that describe these relations are directly measured on small soil cores, and many cores are needed to upscale to the entire heterogeneous flow field. An alternative to the forward upscaling method using small samples are inverse upscaling methods that incorporate soft data from geophysical measurements observed directly on the larger flow field. In this paper, we demonstrate that the hydraulic parameters can be obtained from cross borehole ground penetrating radar by measuring the first arrival travel time of electromagnetic waves (represented by raypaths) from stationary antennae during a constant flux infiltration experiment. The formulation and coupling of the hydrological and geophysical models rely on a constant velocity wetting front that causes critical refraction at the edge of the front as it passes by the antennae. During this critical refraction period, the slope of the first arrival data can be used to calculate (1) the wetting velocity and (2) the hydraulic conductivity of the wet (or saturated) soil. If the soil is undersaturated during infiltration, then an estimate of the saturated water content is needed before calculating the saturated hydraulic conductivity. The hydraulic conductivity value is then used in a nonlinear global optimization scheme to estimate the remaining two parameters of a Broadbridge and White soil.  相似文献   

14.
Conservative solute injections were conducted in three first-order montane streams of different geological composition to assess the influence of parent lithology and alluvial characteristics on the hydrological retention of nutrients. Three study sites were established: (1) Aspen Creek, in a sandstone–siltstone catchment with a fine-grained alluvium of low hydraulic conductivity (1·3×10−4 cm/s), (2) Rio Calaveras, which flows through volcanic tuff with alluvium of intermediate grain size and hydraulic conductivity (1·2×10−3 cm/s), and (3) Gallina Creek, located in a granite/gneiss catchment of coarse, poorly sorted alluvium with high hydraulic conductivity (4·1×10−3 cm/s). All sites were instrumented with networks of shallow groundwater wells to monitor interstitial solute transport. The rate and extent of groundwater–surface water exchange, determined by the solute response in wells, increased with increasing hydraulic conductivity. The direction of surface water–groundwater interaction within a stream was related to local variation in vertical and horizontal hydraulic gradients. Experimental tracer responses in the surface stream were simulated with a one-dimensional solute transport model with inflow and storage components (OTIS). Model-derived measures of hydrological retention showed a corresponding increase with increasing hydraulic conductivity. To assess the temporal variability of hydrological retention, solute injection experiments were conducted in Gallina Creek under four seasonal flow regimes during which surface discharge ranged from baseflow (0·75 l/s in October) to high (75 l/s during spring snowmelt). Model-derived hydrological retention decreased with increasing discharge. The results of our intersite comparison suggest that hydrological retention is strongly influenced by the geologic setting and alluvial characteristics of the stream catchment. Temporal variation in hydrological retention at Gallina Creek is related to seasonal changes in discharge, highlighting the need for temporal resolution in studies of the dynamics of surface water–groundwater interactions in stream ecosystems. © 1997 by John Wiley & Sons, Ltd.  相似文献   

15.
The Beerkan method based on in situ single‐ring water infiltration experiments along with the relevant specific Beerkan estimation of soil transfer parameters (BEST) algorithm is attractive for simple soil hydraulic characterization. However, the BEST algorithm may lead to erroneous or null values for the saturated hydraulic conductivity and sorptivity especially when there are only few infiltration data points under the transient flow state, either for sandy soil or soils in wet conditions. This study developed an alternative algorithm for analysis of the Beerkan infiltration experiment referred to as BEST‐generalized likelihood uncertainty estimation (GLUE). The proposed method estimates the scale parameters of van Genuchten water retention and Brooks–Corey hydraulic conductivity functions through the GLUE methodology. The GLUE method is a Bayesian Monte Carlo parameter estimation technique that makes use of a likelihood function to measure the goodness‐of‐fit between modelled and observed data. The results showed that using a combination of three different likelihood measurements based on observed transient flow, steady‐state flow and experimental steady‐state infiltration rate made the BEST‐GLUE procedure capable of performing an efficient inverse analysis of Beerkan infiltration experiments. Therefore, it is more applicable for a wider range of soils with contrasting texture, structure, and initial and saturated water content. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

16.
Natural groundwater recharge is inherently difficult to quantify and predict, largely because it comprises a series of processes that are spatially distributed and temporally variable. Infiltration ponds used for managed aquifer recharge (MAR) provide an opportunity to quantify recharge processes across multiple scales under semi-controlled conditions. We instrumented a 3-ha MAR infiltration pond to measure and compare infiltration patterns determined using whole-pond and point-specific methods. Whole-pond infiltration was determined by closing a transient water budget (accounting for inputs, outputs, and changes in storage), whereas point-specific infiltration rates were determined using heat as a tracer and time series analysis at eight locations in the base of the pond. Whole-pond infiltration, normalized for wetted area, rose rapidly to more than 1.0 m/d at the start of MAR operations (increasing as pond stage rose), was sustained at high rates for the next 40 d, and then decreased to less than 0.1 m/d by the end of the recharge season. Point-specific infiltration rates indicated high spatial and temporal variability, with the mean of measured values generally being lower than rates indicated by whole-pond calculations. Colocated measurements of head gradients within saturated soils below the pond were combined with infiltration rates to calculate soil hydraulic conductivity. Observations indicate a brief period of increasing saturated hydraulic conductivity, followed by a decrease of one to two orders of magnitude during the next 50 to 75 d. Locations indicating the most rapid infiltration shifted laterally during MAR operation, and we suggest that infiltration may function as a "variable source area" processes, conceptually similar to catchment runoff.  相似文献   

17.
  相似文献   

18.
Research shows that water repellency is a key hydraulic property that results in reduced infiltration rates in burned soils. However, more work is required in order to link the hydrological behaviour of water repellent soils to observed runoff responses at the plot and hillslope scale. This study used 5 M ethanol and water in disc infiltrometers to quantify the role of macropore flow and water repellency on spatial and temporal infiltration patterns in a burned soil at plot (<10 m2) scale in a wet eucalypt forest in south‐east Australia. In the first summer and winter after wildfire, an average of 70% and 60%, respectively, of the plot area was water repellent and did not contribute to infiltration. Macropores (r > 0·5 mm), comprising just 5·5% of the soil volume, contributed to 70% and 95%, respectively, of the field‐saturated and ponded hydraulic conductivity (Kp). Because flow occurred almost entirely via macropores in non‐repellent areas, this meant that less than 2·5% of the soil surface effectively contributed to infiltration. The hydraulic conductivity increased by a factor of up to 2·5 as the hydraulic head increased from 0 to 5 mm. Due to the synergistic effect of macropore flow and water repellency, the coefficient of variation (CV) in Kp was three times higher in the water‐repellent soil (CV = 175%) than under the simulated non‐repellent conditions (CV = 66%). The high spatial variability in Kp would act to reduce the effective infiltration rate during runoff generation at plot scale. Ponding, which tend to increase with increasing scale, activates flow through macropores and would raise the effective infiltration rates at larger scales. Field experiments designed to provide representative measurements of infiltration after fire in these systems must therefore consider both the inherent variability in hydraulic conductivity and the variability in infiltration caused by interactions between surface runoff and hydraulic conductivity. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

19.
A pore‐scale model based on measured particle size distributions has been used to quantify the changes in pore space geometry of packed soil columns resulting from a dilution in electrolyte concentration from 500 to 1 mmol l?1 NaCl during leaching. This was applied to examine the effects of particle release and re‐deposition on pore structure and hydraulic properties. Two different soils, an agricultural soil and a mining residue, were investigated with respect to the change in hydraulic properties. The mining residue was much more affected by this process with the water saturated hydraulic conductivity decreasing to 0·4% of the initial value and the air‐entry value changing from 20 to 50 cm. For agricultural soil, there was little detectable shift in the water retention curve but the saturated hydraulic conductivity decreased to 8·5% of the initial value. This was attributed to localized pore clogging (similar to a surface seal) affecting hydraulic conductivity, but not the microscopically measured pore‐size distribution or water retention. We modelled the soil structure at the pore scale to explain the different responses of the two soils to the experimental conditions. The size of the pores was determined as a function of deposited clay particles. The modal pore size of the agricultural soil as indicated by the constant water retention curve was 45 µm and was not affected by the leaching process. In the case of the mining residue, the mode changed from 75 to 45 µm. This reduction of pore size corresponds to an increase of capillary forces that is related to the measured shift of the water retention curve. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

20.
Previous studies have shown that water retention curve (WRC) and the hydraulic conductivity vary because of changes of the void ratio or porosity of soil. However, limited documents pointed out the change of hydraulic properties of soil when compacted to different porosities while considering both of the drying and wetting processes of the WRC. This information is sometimes necessary for research like finger flow analysis or the occurrence of wetting and drying cycles as what would be seen in the field. Therefore, this study aims to examine the change of WRC characteristics with varied porosity considering both of the drying and wetting path in WRC by conducting a sand box experiment. Results show that the same type of sand compacted to various porosities have different hydraulic parameters. Hydraulic conductivities generally decrease with reduced porosities; shape parameter α of the van Genuchten equation (1980) linearly decreases with declining porosity and shape parameter n in a reversal manner for the sands of interest whether in the drying process or wetting process. The unsaturated properties of sand are further characterized by inspecting the variations of moisture content, matric suction and vertical displacement of soil body subject to periodic changes of the water level by another sand box experiment. The outcomes suggest that the saturated water content and residual water content are changing during the wetting–drying process, which can be an implication of the changed properties of WRC. The characteristics of volumetric deformation might be varied as well because of the observation of the dissimilar patterns of the changing vertical displacements among each wetting–drying process. Infiltration patterns of the sands also are identified through numerical modelling by introducing a constant infiltration flux from the surface followed by a no‐influx condition. Results indicate that less water accumulates in the sand near the surface for the sand compacted to higher porosity, but water can move deeper. Hydraulic conductivity is found as the prime factor dominating the evolvement of wetting fronts. However, shape parameters of water retention curves also affect the infiltration pattern to some extent. In addition, different sands with similar porosities can have quite different infiltrating characteristics. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

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

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