Near‐saturated surface soil hydraulic properties under different land uses in the St Denis National Wildlife Area,Saskatchewan, Canada

Surface soil hydraulic properties are key factors controlling the partition of rainfall and snowmelt into runoff and soil water storage, and their knowledge is needed for sound land management. The objective of this study was to evaluate the effects of three land uses (native grass, brome grass and cultivated) on surface soil hydraulic properties under near‐saturated conditions at the St Denis National Wildlife Area, Saskatchewan, Canada. For each land use, water infiltration rates were measured using double‐ring and tension infiltrometers at ?0·3, ?0·7, ?1·5 and ?2·2 kPa pressure heads. Macroporosity and unsaturated hydraulic properties of the surface soil were estimated. Mean field‐saturated hydraulic conductivity (K_fs), unsaturated hydraulic conductivity at ?0·3 kPa pressure head, inverse capillary length scale (α) and water‐conducting macroporosity were compared for different land uses. These parameters of the native grass and brome grass sites were significantly (p < 0·1) higher than that of the cultivated sites. At the ?0·3 kPa pressure head, hydraulic conductivity of grasslands was two to three times greater than that of cultivated lands. Values of α were about two times and values of K_fs about four times greater in grasslands than in cultivated fields. Water‐conducting macroporosity of grasslands and cultivated fields were 0·04% and 0·01% of the total soil volume, respectively. Over 90% of the total water flux at ?0·06 kPa pressure head was transmitted through pores > 1·36 × 10^?4 m in diameter in the three land uses. Land use modified near‐saturated hydraulic properties of surface soil and consequently may alter the water balance of the area by changing the amount of surface runoff and soil water storage. Copyright © 2004 John Wiley & Sons, Ltd.     

Influence of Degree of Saturation in the Electric Resistivity–Hydraulic Conductivity Relationship

The relationship between aquifer hydraulic conductivity and aquifer resistivity, either measured on the ground surface by vertical electrical sounding (VES) or from resistivity logs, or measured in core samples have been published for different types of aquifers in different locations. Generally, these relationships are empirical and semi-empirical, and confined in few locations. This relation has a positive correlation in some studies and negative in others. So far, there is no potentially physical law controlling this relation, which is not completely understood. Electric current follows the path of least resistance, as does water. Within and around pores, the model of conduction of electricity is ionic and thus the resistivity of the medium is controlled more by porosity and water conductivity than by the resistivity of the rock matrix. Thus, at the pore level, the electrical path is similar to the hydraulic path and the resistivity should reflect hydraulic conductivity. We tried in this paper to study the effect of degree of groundwater saturation in the relation between hydraulic conductivity and bulk resistivity via a simple numerical analysis of Archie’s second law and a simplified Kozeny-Carmen equation. The study reached three characteristic non-linear relations between hydraulic conductivity and resistivity depending on the degree of saturation. These relations are: (1) An inverse power relation in fully saturated aquifers and when porosity equals water saturation, (2) An inverse polynomial relation in unsaturated aquifers, when water saturation is higher than 50%, higher than porosity, and (3) A direct polynomial relation in poorly saturated aquifers, when water saturation is lower than 50%, lower than porosity. These results are supported by some field scale relationships.     

Solute movement through intact columns of cryoturbated Upper Chalk

Cryoturbated Upper Chalk is a dichotomous porous medium wherein the intra‐fragment porosity provides water storage and the inter‐fragment porosity provides potential pathways for relatively rapid flow near saturation. Chloride tracer movement through 43 cm long and 45 cm diameter undisturbed chalk columns was studied at water application rates of 0·3, 1·0, and 1·5 cm h^?1. Microscale heterogeneity in effluent was recorded using a grid collection system consisting of 98 funnel‐shaped cells each 3·5 cm in diameter. The total porosity of the columns was 0·47 ± 0·02 m³ m^?3, approximately 13% of pores were ≥ 15 µm diameter, and the saturated hydraulic conductivity was 12·66 ± 1·31 m day^?1. Although the column remained unsaturated during the leaching even at all application rates, proportionate flow through macropores increased as the application rate decreased. The number of dry cells (with 0 ml of effluent) increased as application rate decreased. Half of the leachate was collected from 15, 19 and 22 cells at 0·3, 1·0, 1·5 cm h^?1 application rates respectively. Similar breakthrough curves (BTCs) were obtained at all three application rates when plotted as a function of cumulative drainage, but they were distinctly different when plotted as a function of time. The BTCs indicate that the columns have similar drainage requirement irrespective of application rates, as the rise to the maxima (C/C_o) is almost similar. However, the time required to achieve that leaching requirement varies with application rates, and residence time was less in the case of a higher application rate. A two‐region convection–dispersion model was used to describe the BTCs and fitted well (r² = 0·97–0·99). There was a linear relationship between dispersion coefficient and pore water velocity (correlation coefficient r = 0·95). The results demonstrate the microscale heterogeneity of hydrodynamic properties in the Upper Chalk. Copyright © 2007 John Wiley & Sons, Ltd.     

Investigating peat hydrological properties using field and laboratory methods: application to the Lanoraie peatland complex (southern Quebec,Canada)

Flow dynamics within a peatland are governed by hydraulic parameters such as hydraulic conductivity, dispersivity and specific yield, as well as by anisotropy and heterogeneity. The aim of this study is to investigate hydraulic parameters variability in peat through the use of different field and laboratory methods. An experimental site located in the Lanoraie peatland complex (southern Quebec, Canada) was used to test the different approaches. Slug and bail tests were performed in piezometer standpipes to investigate catotelm hydraulic conductivity. Combined Darcy tests and tracer experiments were conducted on cubic samples using the modified cube method (MCM) to assess catotelm hydraulic conductivity, anisotropy and dispersivity. A new laboratory method is proposed for assessing acrotelm hydraulic conductivity and gravity drainage using a laboratory experimental tank. Most of slug tests' recovery curves were characteristic of compressible media, and important variability was observed depending on the initial head difference. The Darcy experiments on cubic samples provided reproducible results, and anisotropy (K_h > K_v) was observed for most of samples. All tracer experiments displayed asymmetrical breakthrough curves, suggesting the presence of retardation and/or dual porosity. Hydraulic conductivity estimates performed using the experimental tank showed K variations over a factor of 44 within the upper 40 cm of the acrotelm. The results demonstrate that the intrinsic variability associated with the different field and laboratory methods is small compared with the spatial variability of hydraulic parameters. It is suggested that a comprehensive assessment of peat hydrological properties can be obtained through the combined use of complementary field and laboratory investigations. Copyright © 2007 John Wiley & Sons, Ltd.     

TDR pressure cell for monitoring water content retention and bulk electrical conductivity curves in undisturbed soil samples

The water retention curve (θ(ψ)), which defines the relationship between soil volumetric water content (θ) and matric potential (ψ), is of paramount importance in characterizing the hydraulic behaviour of soils. However, few methods are so far available for estimating θ(ψ) in undisturbed soil samples. We present a new design of TDR‐pressure cell (TDR‐Cell) for estimating θ(ψ) in undisturbed soil samples. The TDR‐Cell consists of a 50‐mm‐long and 50‐mm internal diameter stainless steel cylinder (which constitutes the outer frame of a coaxial line) attached to a porous ceramic disc and closed at the ends with two aluminium lids. A 49‐mm‐long and 3‐mm‐diameter stainless steel rod, which runs longitudinally through the centre of the cylinder, constitutes the inner rod of a coaxial TDR probe. The TDR‐Cell was used to determine the θ(ψ) curves of a packed sand and seven undisturbed soil samples from three profiles of agricultural soils. These θ(ψ) curves were subsequently compared to those obtained from the corresponding 2‐mm sieved soils using the pressure plate method. Measurements of bulk electrical conductivity, σ_a, as a function of the water content, σ_a(θ), of the undisturbed soil samples were also performed. An excellent correlation (R² = 0·988) was found between the θ values measured by TDR on the different undisturbed soils and the corresponding θ obtained from the soil gravimetric water content. A typical bimodal θ(ψ) function was found for most of the undisturbed soil samples. Comparison between the θ(ψ) curves measured with the TDR‐Cell and those obtained from the 2‐mm sieved soils showed that the pressure plate method overestimates θ at low ψ values. The σ_a(θ) relationship was well described by a simple power expression (R² > 0·95), in which the power factor, defined as tortuosity, ranged between 1·18 and 3·75. Copyright © 2011 John Wiley & Sons, Ltd.     

Uncertainties in DRAINMOD predictions of subsurface drain flow for an Indiana silt loam using the GLUE methodology

Good modelling practice requires the incorporation of uncertainty analysis into hydrologic/water quality models. The generalized likelihood uncertainty estimation procedure was used to evaluate the uncertainty in DRAINMOD predictions of daily, monthly, and yearly subsurface drain flow. A variance‐based sensitivity analysis technique, the extended Fourier amplitude sensitivity test, was used to identify the main sources of prediction uncertainty. The analysis was conducted for the experimental drainage field at the Southeast Purdue Agricultural Center in Indiana. Six years of data were used and the uncertainties in eight model parameters were considered to analyse how uncertainties in input parameters propagate to model outputs. The width of 90% confidence interval bands of drain flow ranged from 0 to 0·6 cm day^?1 for daily predictions, from 0 to 3·1 cm month^?1 for the monthly predictions, and from 7·6 to 12·4 cm year^?1 for yearly predictions. Annual drain flow predicted by DRAINMOD fell well within the 90% confidence bounds. Model results were most sensitive to the vertical saturated hydraulic conductivity of the restrictive layer and the lateral hydraulic conductivity of the deepest soil layer, followed by the lateral hydraulic conductivity of the top soil layer and surface micro‐storage. Parameter interactions also contributed to the prediction uncertainty. Copyright © 2006 John Wiley & Sons, Ltd.     

Heterogeneous characteristics of streambed saturated hydraulic conductivity of the Touchet River,south eastern Washington,USA

Traditionally a streambed is treated as a layer of uniform thickness and low saturated hydraulic conductivity (K) in surface‐ and ground‐water studies. Recent findings have shown a high level of spatial heterogeneity within a streambed and such heterogeneity directly affects surface‐ and ground‐water exchange and can have ecological implications for biogeochemical transformations, nutrient cycling, organic matter decomposition, and reproduction of gravel spawning fish. In this study a detailed field investigation of K was conducted in two selected sites in Touchet River, a typical salmon spawning stream in arid south eastern Washington, USA. In‐stream slug tests were conducted to determine K following the Bouwer and Rice method. For the upper and lower sites, each 50 m long and 9 m wide and roughly 20 m apart, a sampling grid of 5 m longitudinally and 3 m transversely was used. The slug tests were performed for each horizontal coordinate at 0·3–0·45, 0·6–0·75, 0·9–1·05 and 1·2–1·35 m depth intervals unless a shallower impenetrable obstruction was encountered. Additionally, water levels were measured to obtain vertical hydraulic gradient (VHG) between each two adjacent depth intervals. Results indicated that K ranged over three orders of magnitude at both the upper and lower sites and differed between the two sites. At the upper site, K did not differ significantly among different depth intervals based on nonparametric statistical tests for mean, median, and empirical cumulative distribution, but the spatial pattern of K varied among different depth intervals. At the lower site, K for the 0·3–0·45 m depth interval differed statistically from those at other depth intervals, and no similar spatial pattern was found among different depth intervals. Zones of upward and downward water flow based on VHG also varied among different depth intervals, reflecting the complexities of the water flow regime. Detailed characterization of the streambed as attempted in this study should be helpful in providing information on spatial variations of streambed hydraulic properties as well as surface‐ and ground‐water interaction. Copyright © 2009 John Wiley & Sons, Ltd.     

Effects of meteorological and hydrogeological factors on gross recharge percentage at unconfined sandy aquifers with an equatorial climate

A thorough understanding of rainfall recharge processes and their controlling factors is essential for management of groundwater systems. This study investigates the effects of various meteorological and hydrogeological factors on the gross recharge percentages, the rainfall–recharge relationships and the recharge threshold values for unconfined sandy aquifers under an equatorial climate. Among the meteorological factors investigated, rainfall intensity was found to have the most significant impact on the gross recharge rate. The effects of potential evaporation rate, relative humidity and air temperature on the gross recharge percentage were significant when the vadose zone thickness is larger than 2·5 m. The recharge threshold values were found to depend strongly on the vadose zone thickness. The rainfall–recharge relationships could generally be well defined by a normal–log relationship. The rainfall–recharge relationships derived here are applicable to yield estimates of gross recharge percentages for unconfined sandy aquifers under an equatorial climate, using rainfall intensity and vadose zone thickness as input variables. In this study, a theory was developed and validated to provide physical explanations for the observations, based on the residence time of the percolated rainwater within the vadose zone. Among the soil hydraulic parameters tested, porosity and saturated hydraulic conductivity were found to have the most pronounced effects on the gross recharge percentage. Utilizing the sensitivity results and the theory derived, an approach was developed for extending the application of the derived rainfall–recharge relationships to other sand textures. The approach was found to be capable of producing rough and fast estimations of gross recharge percentage for other sand textures. Copyright © 2007 John Wiley & Sons, Ltd.     

Correlating specific conductivity with total hardness in gypsum karst waters

Modern digital conductivity meters are readily portable, robust, cheap, and give precisely reproducible values of specific electrical conductivity (SpC, in µS cm^?1). Here we investigate the accuracy of their estimates of the amounts of gypsum dissolved in waters collected in gypsum karst terrains, expressed as total hardness (TH) in mg L^?1 of CaSO₄·2H₂O (GYP). Total dissolved solid concentrations (TDS) are also considered. Curves obtained with the program PHREEQC, for the dissolution of pure gypsum in water at 25 C, are compared with 574 comprehensive water chemical analyses selected from gypsum karst studies in Europe and the Americas. Principal common and foreign ions encountered are the BNC group (bicarbonates, nitrates, chlorides). It is found that GYP = 1·12^·SpC + 62 where BNC < 33% (Cl^? < 5%), with one standard error <5% for waters with SpC > 2400 µS cm^?1; GYP = 0·74^·SpC + 777 where BNC < 33% (5% ≤ Cl^? < 15%), with one standard error <10% for waters with SpC > 3100 µS cm^?1; GYP = 0·97^·SpC ? 209 where BNC < 33% and Cl^? ≥ 15%, with one standard error <10% for samples with SpC > 4300 µS cm^?1. There are similar results for the more complex waters found in gypsum karsts where much carbonate rock or salt is also present, to the limit of BNC < 50% for what may reasonably be defined as ‘gypsum waters’. Values of R² for linear correlations of different subsets of the water samples range from 0·69 to 0·96, the majority being >0·8. Copyright © 2006 John Wiley & Sons, Ltd.     

Hydraulic conductivity in frozen soils

An apparatus has been developed for investigation of hydraulic conductivity of frozen soils. The test procedure is isothermal and involves the passage of water from one reservoir into the frozen sample and out of the frozen sample into a second reservoir. The water in the reservoirs remains unfrozen because it contains dissolved lactose. The concentration of lactose is such that, initially, the water in the reservoirs is in thermodynamic equilibrium with the water in the soil. On application of pressure to one reservoir a known hydraulic gradient is established and flow takes place. Flow is shown to vary linearly with hydraulic gradient. The hydraulic conductivity coefficient depends on soil type and temperature and is related to the unfrozen water content. At temperatures within a few tenths of 0°C the coefficient apparently ranges from 10^?5 to 10^?9 cm sec^?1, and decreases only slowly below about ?0·5°C. Soils known to be susceptible to frost heave are shown to have significant hydraulic conductivities well below 0°C.     

Synergistic effects of water repellency and macropore flow on the hydraulic conductivity of a burned forest soil,south‐east Australia

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 m²) 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 (K_p). 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 K_p 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 K_p 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.     

Determination of Hydraulic Conductivity from Grain‐Size Distribution for Different Depositional Environments

Over 400 unlithified sediment samples were collected from four different depositional environments in global locations and the grain‐size distribution, porosity, and hydraulic conductivity were measured using standard methods. The measured hydraulic conductivity values were then compared to values calculated using 20 different empirical equations (e.g., Hazen, Carman‐Kozeny) commonly used to estimate hydraulic conductivity from grain‐size distribution. It was found that most of the hydraulic conductivity values estimated from the empirical equations correlated very poorly to the measured hydraulic conductivity values with errors ranging to over 500%. To improve the empirical estimation methodology, the samples were grouped by depositional environment and subdivided into subgroups based on lithology and mud percentage. The empirical methods were then analyzed to assess which methods best estimated the measured values. Modifications of the empirical equations, including changes to special coefficients and addition of offsets, were made to produce modified equations that considerably improve the hydraulic conductivity estimates from grain size data for beach, dune, offshore marine, and river sediments. Estimated hydraulic conductivity errors were reduced to 6 to 7.1 m/day for the beach subgroups, 3.4 to 7.1 m/day for dune subgroups, and 2.2 to 11 m/day for offshore sediments subgroups. Improvements were made for river environments, but still produced high errors between 13 and 23 m/day.     

ALLUVIAL CHARACTERISTICS,GROUNDWATER–SURFACE WATER EXCHANGE AND HYDROLOGICAL RETENTION IN HEADWATER STREAMS

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⁻⁴ cm/s), (2) Rio Calaveras, which flows through volcanic tuff with alluvium of intermediate grain size and hydraulic conductivity (1·2×10⁻³ cm/s), and (3) Gallina Creek, located in a granite/gneiss catchment of coarse, poorly sorted alluvium with high hydraulic conductivity (4·1×10⁻³ 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.     

Using a pore‐scale model to quantify the effect of particle re‐arrangement on pore structure and hydraulic properties

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.     

Dynamics of deformation and water flow in heterogeneous porous media and its impact on soil properties

The interaction of geomechanics and flow within a soil body induces deformation and pore pressure change. Deformation may change hydrogeological and elastic properties, which alters the mechanical behaviour and results in non‐linearity. To investigate this interaction effect in a heterogeneous porous medium, a stochastic poroelastic model is proposed. Monte Carlo simulations are performed to determine the mean and uncertainty of the parameter changes, displacement, and change in pore water pressure. Hydraulic conductivity is treated as the only random variable in the coupled geomechanics‐flow system due to its large variation compared to other mechanical and hydrogeological properties in natural environments. The three considered non‐linear models for the interaction between parameters and deformation are those that consider (1) porosity and hydraulic conductivity; (2) porosity and Young's modulus; and (3) a combined effect that includes porosity, hydraulic conductivity, and Young's modulus. Boundary effects on the coupled system are also explored. The relationships between changes of porosity, hydraulic conductivity, and Young's modulus are analytically shown to be non‐linear. Among the considered parameters, the deformation effect induces the largest reduction in hydraulic conductivity. The deformation‐induced change in hydraulic conductivity shows the most significant effect on the mean and variance of the change in pore water pressure and displacement, while changes in Young's modulus have the least effect. When the deformation effect is considered, the superposition relationship does not exist in the mean displacement and mean change in pore water pressure for the three scenarios considered; it exists for the case without deformation effects. Deformation also causes a reduction in the effective hydraulic conductivity for the whole domain. The scenario that considers both loading and discharge boundaries has larger changes in hydrogeological and geo‐mechanical parameters than those in scenarios that consider loading and discharge boundaries separately. The results indicate that the interaction between deformation and changes in parameters has a profound effect on the poroelastic system. The effect of deformation should thus be considered in modelling and practice. Copyright © 2009 John Wiley & Sons, Ltd.     

Estimation of anisotropic hydraulic conductivity using geophysical data in a coastal aquifer of Karnataka,India

Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.     

Characteristics of preferential flow and groundwater discharge to Shingobee Lake,Minnesota, USA

Small‐scale heterogeneities and large changes in hydraulic gradient over short distances can create preferential groundwater flow paths that discharge to lakes. A 170 m² grid within an area of springs and seeps along the shore of Shingobee Lake, Minnesota, was intensively instrumented to characterize groundwater‐lake interaction within underlying organic‐rich soil and sandy glacial sediments. Seepage meters in the lake and piezometer nests, installed at depths of 0·5 and 1·0 m below the ground surface and lakebed, were used to estimate groundwater flow. Statistical analysis of hydraulic conductivity estimated from slug tests indicated a range from 21 to 4·8 × 10^?3 m day^?1 and small spatial correlation. Although hydraulic gradients are overall upward and toward the lake, surface water that flows onto an area about 2 m onshore results in downward flow and localized recharge. Most flow occurred within 3 m of the shore through more permeable pathways. Seepage meter and Darcy law estimates of groundwater discharge agreed well within error limits. In the small area examined, discharge decreases irregularly with distance into the lake, indicating that sediment heterogeneity plays an important role in the distribution of groundwater discharge. Temperature gradients showed some relationship to discharge, but neither temperature profiles nor specific electrical conductance could provide a more convenient method to map groundwater–lake interaction. These results suggest that site‐specific data may be needed to evaluate local water budget and to protect the water quality and quantity of discharge‐dominated lakes. Copyright © 2002 John Wiley & Sons, Ltd.     

Assessment of runoff and sediment yield in the Miyun Reservoir catchment by using SWAT model

The Soil and Water Assessment Tool (SWAT) was used to simulate the transport of runoff and sediment into the Miyun Reservoir, Beijing in this study. The main objective was to validate the performance of SWAT and the feasibility of using this model as a simulator of runoff and sediment transport processes at a catchment scale in arid and semi‐arid area in North China, and related processes affecting water quantity and soil erosion in the catchment were simulated. The investigation was conducted using a 6‐year historical streamflow and sediment record from 1986 to 1991; the data from 1986 to 1988 was used for calibration and that from 1989 to 1991 for validation. The SWAT generally performs well and could accurately simulate both daily and monthly runoff and sediment yield. The simulated daily and monthly runoff matched the observed values satisfactorily, with a Nash‐Sutcliffe coefficient of greater than 0·6, 0·9 and a coefficient of determination 0·75, 0·9 at two outlet stations (Xiahui and Zhangjiafen stations) during calibration. These values were 0·6, 0·85 and 0·6, 0·9 during validation. For sediment simulation, the efficiency is lower than that for runoff. Even so, the Nash‐Sutcliffe coefficient and coefficient of determination were greater than 0·48 and 0·6 for monthly sediment yield during calibration, and these values were greater than 0·84 and 0·95 during validation. Sensitivity analysis shows that sensitive parameters for the simulation of discharge and sediment yield include curve number, base flow alpha factor, soil evaporation compensation factor, soil available water capacity, soil profile depth, surface flow lag time and channel re‐entrained linear parameter, etc. Copyright © 2009 John Wiley & Sons, Ltd.     

Influence of pore size and geometry on peat unsaturated hydraulic conductivity computed from 3D computed tomography image analysis

In organic soils, hydraulic conductivity is related to the degree of decomposition and soil compression, which reduce the effective pore diameter and consequently restrict water flow. This study investigates how the size distribution and geometry of air‐filled pores control the unsaturated hydraulic conductivity of peat soils using high‐resolution (45 µm) three‐dimensional (3D) X‐ray computed tomography (CT) and digital image processing of four peat sub‐samples from varying depths under a constant soil water pressure head. Pore structure and configuration in peat were found to be irregular, with volume and cross‐sectional area showing fractal behaviour that suggests pores having smaller values of the fractal dimension in deeper, more decomposed peat, have higher tortuosity and lower connectivity, which influences hydraulic conductivity. The image analysis showed that the large reduction of unsaturated hydraulic conductivity with depth is essentially controlled by air‐filled pore hydraulic radius, tortuosity, air‐filled pore density and the fractal dimension due to degree of decomposition and compression of the organic matter. The comparisons between unsaturated hydraulic conductivity computed from the air‐filled pore size and geometric distribution showed satisfactory agreement with direct measurements using the permeameter method. This understanding is important in characterizing peat properties and its heterogeneity for monitoring the progress of complex flow processes at the field scale in peatlands. Copyright © 2010 John Wiley & Sons, Ltd.     

Water flow path interactions with soil hydraulic properties in till soil at Gårdsjön, Sweden

In 1989, in a hydrological research programme within a deacidification project in the Gårdsjön area in southwest Sweden, flow paths and residence times of soil water and groundwater in microcatchments were examined to support the interpretation of the hydrochemical changes. Saturated hydraulic conductivity and soil water retention were analysed on more than 100 cylinder samples. The catchments have shallow sandy-silty till soil with a mean depth in the main catchment of 43 cm. Porosity of the mineral soil in the main catchment was high and ranged from 38 to 85%. The samples from the B-horizon had generally higher porosity. Porosity and the content of organic matter were correlated. The soil water retention was relatively high at all tensions, likely owing to the high content of organic matter. Dissolved organic substances were most probably transported from the shallow soil on the steep sides of the catchment down to the valley where it precipitated. The high porosities could be a consequence of long-term weathering, provided that the organic substances present have increased the leaching of the weathering products. Measured values of saturated hydraulic conductivity were close to log-normally distributed with a mean for all samples of 3 × 10⁻⁵ m s⁻¹. There was a significant increase in conductivity toward the ground surface with the mean conductivity of the samples in the uppermost 10 cm of the mineral soil of 4 × 10⁻⁵ m s⁻¹, which was about 13 times higher than the conductivity of 3 × 10⁻⁶ m s⁻¹ at 1 m depth. From the relationship between runoff at the catchment outlet and groundwater levels, the conductivity was estimated to be 15–200 times higher in the upper soil layer than in the deeper ones. In one profile, 44–64% of the yearly lateral flow was estimated to occur above 30 cm depth. The conductivity was correlated with the content of drainable water, which indicated the importance of the largest pores for the saturated hydraulic conductivity.