High-resolution saturated hydraulic conductivity logging of borehole cores using air permeability measurements

Saturated hydraulic conductivity (K _s) is one of the most important parameters determining groundwater flow and contaminant transport in both unsaturated and saturated porous media. The hand-held air permeameter technique was investigated for high-resolution hydraulic conductivity determination on borehole cores using a spatial resolution of ～0.05 m. The suitability of such air permeameter measurements on friable to poorly indurated sediments was tested to improve the spatial prediction of classical laboratory-based K _s measurements obtained at a much lower spatial resolution (～2 m). In total, 368 K _s measurements were made on ～350 m of borehole cores originating from the Campine basin, northern Belgium, while ～5,230 air permeability measurements were performed on the same cores, resulting in a K _s range of seven orders of magnitude. Cross-validation demonstrated that, using air permeameter data as the secondary variable for laboratory based K _s measurements, the performance increased from R ²?=?0.35 for ordinary kriging (laboratory K _s only) to R ²?=?0.61 for co-kriging. The separate treatment of horizontal and vertical hydraulic conductivity revealed considerable anisotropy in certain lithostratigraphical units, while others were clearly isotropic at the sample scale. Air permeameter measurements on borehole cores provide a cost-effective way to improve spatial predictions of traditional laboratory based K _s.     

Measurement of streambed hydraulic conductivity and its anisotropy

 A method is described for the measurement of streambed hydraulic conductivity. Unlike permeameter methods, this method applies straight and l-shaped standpipes directly to streambeds for measurements of in-situ hydraulic conductivity in the vertical (K _v ) and horizontal (K _h ) directions, as well as in other oblique directions (K _s ). This method has advantages in determination of K _v values over grain-size analysis, permeameter tests, or slug test methods. Also unique to this method is that it provides K _s values of a streambed. The measured results can be used to construct a hydraulic conductivity ellipse and to evaluate the anisotropy of streambed sediments. Field examples from the Republican River, Nebraska, demonstrated the usefulness of this method in the determination of streambed hydraulic conductivity and anisotropy along or across a river channel. Results indicate that the K _h is about three to four times larger than K _v , whereas K _s values are larger than K _v but smaller than K _h . Received: 6 March 2000 · Accepted: 18 April 2000     

Uncertainty of rainfall-induced landslides considering spatial variability of parameters

A cross-correlation analysis is conducted to determine the impacts of the heterogeneity of hydraulic conductivity K_s, soil cohesion c′ and soil friction angle (tan φ′) on the uncertainty of slope stability in time and space during rainfall. We find the relative importance of tan φ′ and c′ depends on the effective stress. While the sensitivity of the stability to the variability of K_s is small, the large coefficient of variation of K_s may exacerbate the variability of pore-water pressure. Therefore, characterizing the heterogeneity of hydraulic properties and pore-water distribution in the field is critical to the stability analysis.     

Estimating regional-scale fractured bedrock hydraulic conductivity using discrete fracture network (DFN) modeling

Estimating bedrock hydraulic conductivity of regional fractured aquifers is challenging due to a lack of aquifer testing data and the presence of small and large-scale heterogeneity. This study provides a novel approach for estimating the bedrock hydraulic conductivity of a regional-scale fractured bedrock aquifer using discrete fracture network (DFN) modeling. The methodology is tested in the mountainous Okanagan Basin, British Columbia, Canada. Discrete fractures were mapped in outcrops, and larger-scale fracture zones (corresponding to lineaments) were mapped from orthophotos and LANDSAT imagery. Outcrop fracture data were used to generate DFN models for estimating hydraulic conductivity for the fractured matrix (K _m). The mountain block hydraulic conductivity (K _mb) was estimated using larger-scale DFN models. Lineament properties were estimated by best fit parameters for a simulated pumping test influenced by a fracture zone. Unknown dip angles and directions for lineaments were estimated from the small-scale fracture sets. Simulated K _m and K _mb values range from 10^–8 to 10^–7?m/s and are greatest in a N–S direction, coinciding with the main strike direction of Okanagan Valley Fault Zone. K _mb values also decrease away from the fault, consistent with the decrease in lineament density. Simulated hydraulic conductivity values compare well with those estimated from pumping tests.     

Impact of pore water conductivity and porosity on the electrical conductivity of kaolinite

The purpose of this study is to quantify the magnitudes of surface conduction and pore water conduction from the measured electrical conductivity of kaolinite, with the ultimate goal of estimating the electrical conductivity of kaolinite with a wide range of pore water conductivities (σ _w = 0.013–3.356 S/m) and porosities (n = 0.368–1.0). Therefore, the theoretical background of the electrical conductivity in soils was reviewed, and electrical conductivity measurements on kaolinite were performed using both slurry and consolidation tests in this study. The results of this study demonstrate that the variations of measured electrical conductivity (σ _mix) with n are debatable according to the values of σ _w, because a decrease in n results in both an increase in surface conduction (K _s) and a decrease in pore water conduction (K _w); this causes the relative magnitude of K _s compared to that of K _w to vary with σ _w and n. Consequently, this study develops the relation between the porosity-normalized K _s/K _w and 1/σ _w. Additionally, the surface conductivity of the tested kaolinite is back-calculated and compared with the previous relationship between K _s and zeta potential of kaolinite. The measured and estimated σ _mix values are compared with the varying pore water conductivity and porosity values.     

Heterogeneity of hydraulic conductivity and Darcian flux in the submerged streambed and adjacent exposed stream bank of the Beiluo River,northwest China

Recognizing the heterogeneity of hydraulic conductivity and hyporheic flow is critical for understanding contaminant transfer and biogeochemical and hydrological processes involving streams and aquifers. In this study, the heterogeneity of hydraulic conductivity and Darcian flux in a submerged streambed and its adjacent exposed stream banks were investigated in the Beiluo River, northwest China. In the submerged streambed, Darcian flux was estimated by measurement of vertical hydraulic conductivity (K _v) and vertical head gradient (VHG) using in-situ permeameter tests. On exposed stream banks, both horizontal hydraulic conductivity (K _h) and K _v were measured by on-site permeameter tests. In the submerged streambed, K _v values gradually decreased with depth and the higher values were concentrated in the center and close to the erosional bank. Compared to the exposed stream banks, the K _v values were higher in the streambed. From stream stage to the topmost layer of tested sediment, through increasing elevation, the K _h values increased on the erosional bank, while they decreased on the depositional bank. The values of VHG along the thalweg illustrate that downwelling flux occurred in the deepest area while upwelling flux appeared in the other areas, which might result from the change of streambed elevation. The higher value of the Darcian flux in the submerged streambed existed near the erosional bank.     

Estimation of Hydraulic Conductivity and Its Uncertainty from Grain-Size Data Using GLUE and Artificial Neural Networks

Various approaches exist to relate saturated hydraulic conductivity (K _s) to grain-size data. Most methods use a single grain-size parameter and hence omit the information encompassed by the entire grain-size distribution. This study compares two data-driven modelling methods??multiple linear regression and artificial neural networks??that use the entire grain-size distribution data as input for K _s prediction. Besides the predictive capacity of the methods, the uncertainty associated with the model predictions is also evaluated, since such information is important for stochastic groundwater flow and contaminant transport modelling. Artificial neural networks (ANNs) are combined with a generalised likelihood uncertainty estimation (GLUE) approach to predict K _s from grain-size data. The resulting GLUE-ANN hydraulic conductivity predictions and associated uncertainty estimates are compared with those obtained from the multiple linear regression models by a leave-one-out cross-validation. The GLUE-ANN ensemble prediction proved to be slightly better than multiple linear regression. The prediction uncertainty, however, was reduced by half an order of magnitude on average, and decreased at most by an order of magnitude. This demonstrates that the proposed method outperforms classical data-driven modelling techniques. Moreover, a comparison with methods from the literature demonstrates the importance of site-specific calibration. The data set used for this purpose originates mainly from unconsolidated sandy sediments of the Neogene aquifer, northern Belgium. The proposed predictive models are developed for 173 grain-size K _s-pairs. Finally, an application with the optimised models is presented for a borehole lacking K _s data.     

Mesures in situ par infiltrométrie des propriétés hydrodynamiques des sols de Mnasra (Maroc)

The aim of this work is the field hydraulic characterisation of Mnasra soils in northern Morocco, which represents an essential step to study the hydraulic and chemical transports through the vadose zone. We have used a tension infiltrometer associated with a transient axisymmetric infiltration method to determine the hydraulic conductivity, which reduces the duration of measurements. This allows us to characterise a large area with many measurements. Parameters of the characteristic functions K(h) and θ(h) are estimated for six different soils belonging to two geomorphologically different domains: a sandy zone and an alluvial plain. To cite this article: K. Tamoh, A. Maslouhi, C. R. Geoscience 336 (2004).     

Predicting saturated hydraulic conductivity using particle swarm optimization and genetic algorithm

Soil saturated hydraulic conductivity (K_s) is considered as soil basic hydraulic property, and its precision estimation is a key element in modeling water flow and solute transport processes both in the saturated and vadose zones. Although some predictive methods (e.g., pedotransfer functions, PTFs) have been proposed to indirectly predict K_s, the accuracy of these methods still needs to be improved. In this study, some easily available soil properties (e.g., particle size distribution, organic carbon, calcium carbonate content, electrical conductivity, and soil bulk density) are employed as input variables to predict K_s using a fuzzy inference system (FIS) trained by two different optimization techniques: particle swarm optimization (PSO) and genetic algorithm (GA). To verify the derived FIS, 113 soil samples were taken, and their required physical properties were measured (113 sample points?×?7 factors?=?791 input data). The initial FIS is compared with two methods: FIS trained by PSO (PSO-FIS) and FIS trained by GA (GA-FIS). Based on experimental results, all three methods are compared according to some evaluation criteria including correlation coefficient (r), modeling efficiency (EF), coefficient of determination (CD), root mean square error (RMSE), and maximum error (ME) statistics. The results showed that the PSO-FIS model achieved a higher level of modeling efficiency and coefficient of determination (R²) in comparison with the initial FIS and the GA-FIS model. EF and R² values obtained by the developed PSO-FIS model were 0.69 and 0.72, whereas they were 0.63 and 0.54 for the GA-FIS model. Moreover, the results of ME and RMSE indices showed that the PSO-FIS model can estimate soil saturated hydraulic conductivity more accurate than the GA-FIS model with ME?=?10.4 versus 11.5 and RMSE?=?5.2 versus 5.5 for PSO-FIS and GA-FIS, respectively.     

Determination of the anisotropy of an upper streambed layer in east-central Nebraska, USA

Information on the anisotropy of streambed hydraulic conductivity (K) is a necessity for analyses of water exchange and solute transport in the hyporheic zone. An approach is proposed for the determination of K, developed from existing in-situ permeameter test methods. The approach is based on determination of vertical and horizontal hydraulic conductivity of streambed sediments on-site and eliminates the effects of vertical flow in the hyporheic zone and stream-stage fluctuation, which normally influence in situ permeameter tests. The approach was applied to seven study sites on four tributaries of the Platte River in east-central Nebraska, USA. On-site permeameter tests conducted on about 172 streambed cores for the determination of vertical hydraulic conductivity (K _v) and horizontal hydraulic conductivity (K _h) at the study sites indicate that the study sites have relatively small anisotropic ratios, ranging from 0.74 to 2.40. The ratios of K _h to K _v from individual locations within a study site show greater variation than the anisotropic ratios from the mean or median K at each of the study sites.     

Time-Variant Reliability Analysis of Unsaturated Soil Slopes Under Rainfall

Time-variant reliability analysis for a typical unsaturated soil slope is performed. Eight rainfall conditions are considered, and three slope models are set up for studying the influence of shear strength parameters, hydraulic conductivity parameters, rainfall intensity and duration on the reliability of the soil slope. Sensitivity analysis shows that when the saturated hydraulic conductivity (k _s) is very small, the variation of hydraulic conductivity has little effect on the reliability index (β). For saving the computation effort, only the shear strength parameters are needed in performing the reliability analysis in this condition. With the increase of k _s, the importance of hydraulic conductivity becomes large. The reliability index of the soil slope is changing with time (t), and the shape of β–t curves for different slope model is quite different for they depend on the value of k _s. When k _s is very small, β keeps decreasing for all the 18 simulation days. With the increase of k _s, β decreases to its minimum value at about the cessation day of rainfall events, and it then increases gradually due to the redistribution of suction in the soil slope.     

Estimating the hydraulic properties of an aquitard from in situ pore pressure measurements

A workflow is described to estimate specific storage (S _s) and hydraulic conductivity (K) from a profile of vibrating wire piezometers embedded into a regional aquitard in Australia. The loading efficiency, compressibility and S _s were estimated from pore pressure response to atmospheric pressure changes, and K was estimated from the earliest part of the measurement record following grouting. Results indicate that S _s and K were, respectively, 8.8?×?10^?6 to 1.2?×?10^?5 m^?1 and 2?×?10^?12 m s^?1 for a claystone/siltstone, and 4.3?×?10^?6 to 9.6?×?10^?6 m^?1 and 1?×?10^?12 to 5?×?10^?12 m s^?1 for a thick mudstone. K estimates from the pore pressure response are within one order of magnitude when compared to direct measurement in a laboratory and inverse modelled flux rates determined from natural tracer profiles. Further analysis of the evolution and longevity of the properties of borehole grout (e.g. thermal and chemical effects) may help refine the estimation of formation hydraulic properties using this workflow. However, the convergence of K values illustrates the benefit of multiple lines of evidence to support aquitard characterization. An additional benefit of in situ pore pressure measurement is the generation of long-term data to constrain groundwater flow models, which provides a link between laboratory scale data and the formation scale.     

Statistical and geostatistical features of streambed hydraulic conductivities in the Platte River, Nebraska

This paper presents streambed hydraulic conductivities of the Platte River from south-central to eastern Nebraska. The hydraulic conductivities were determined from river channels using permeameter tests. The vertical hydraulic conductivities (K _v ) from seven test sites along this river in south-central Nebraska belong to one statistical population. Its mean value is 40.2 m/d. However, the vertical hydraulic conductivities along four transects of the Ashland test site in eastern Nebraska have lower mean values, are statistically different from the K _v values in south-central Nebraska, and belong to two different populations with mean values of 20.7 and 9.1 m/d, respectively. Finer sediments carried from the Loup River and Elkhorn River watersheds to the eastern reach of the Platte River lowers the vertical hydraulic conductivity of the streambed. Correlation coefficients between water depth and K _v values along a test transect indicates a positive correlation – a larger K _v usually occurs in the part of channel with deeper water. Experimental variograms derived from the vertical hydraulic conductivities for several transects across the channels of the Platte River show periodicity of spatial correlation, which likely result from periodic variation of water depth across the channels. The sandy to gravelly streambed contains very local silt and clay layers; spatially continuous low-permeability streambed was not observed in the river channels. The horizontal hydraulic conductivities were larger than the vertical hydraulic conductivities for the same test locations.     

Comparison of in situ uranium KD values with a laboratory determined surface complexation model

Reactive solute transport simulations in groundwater require a large number of parameters to describe hydrologic and chemical reaction processes. Appropriate methods for determining chemical reaction parameters required for reactive solute transport simulations are still under investigation. This work compares U(VI) distribution coefficients (i.e. K_D values) measured under field conditions with K_D values calculated from a surface complexation model developed in the laboratory. Field studies were conducted in an alluvial aquifer at a former U mill tailings site near the town of Naturita, CO, USA, by suspending approximately 10 g samples of Naturita aquifer background sediments (NABS) in 17-5.1-cm diameter wells for periods of 3 to 15 months. Adsorbed U(VI) on these samples was determined by extraction with a pH 9.45 NaHCO₃/Na₂CO₃ solution. In wells where the chemical conditions in groundwater were nearly constant, adsorbed U concentrations for samples taken after 3 months of exposure to groundwater were indistinguishable from samples taken after 15 months. Measured in situ K_D values calculated from the measurements of adsorbed and dissolved U(VI) ranged from 0.50 to 10.6 mL/g and the K_D values decreased with increasing groundwater alkalinity, consistent with increased formation of soluble U(VI)-carbonate complexes at higher alkalinities. The in situ K_D values were compared with K_D values predicted from a surface complexation model (SCM) developed under laboratory conditions in a separate study. A good agreement between the predicted and measured in situ K_D values was observed. The demonstration that the laboratory derived SCM can predict U(VI) adsorption in the field provides a critical independent test of a submodel used in a reactive transport model.     

The spatial variability of streambed vertical hydraulic conductivity in an intermittent river, northwestern China

Streambed vertical hydraulic conductivity (K) plays an important role in river water and groundwater interaction. The K at the ten transects (Ts1–Ts10) at the Donghe River (an intermittent river) in the Ejina Basin, northwestern China, was measured to investigate its spatial variation. Based on the sediment characteristics and vertical hydraulic conductivity of the riverbed, the entire riverbed of the Donghe River could be divided arbitrarily into two parts: an upper part (starting at Ts1 and ending at Ts9, without an obvious and continuous clogging layer) and a lower part (the remaining riverbed, with an obvious and continuous clogging layer). In the upper part, although the K varied with depth within the 0–30 cm layer, the variability with depth could be ignored in practice. The arithmetic mean K of the upper part ranged from 12 to 27.6 m/day, three orders of magnitude larger than that of the lower part (0.06 m/day). The change of K along the river cross section was significant, and larger values of K often occurred in the parts of the channels with greater water depth. However, there were no consistent patterns of the variability of K at transects across the river, which was influenced by the variation in streambed characteristics. The results could be useful for the estimation of groundwater recharge from river and groundwater resources evaluation in the Ejina Basin.     

Soil-hydrological properties response to grazing exclusion in a steppe grassland of the Loess Plateau

Soil–water characteristics are necessary for water quality monitoring, solute migration and plant growth. Soil–water characteristic curve (SWCC) is a relationship between suction and water content or degree of saturation. However, little information is available concerning the impacts of grazing exclusion management on soil–water characteristics. Here, the soil–water characteristics of grasslands, which were excluded grazing for 5 (GE5) and 15 years (GE15), were studied. The saturated hydraulic conductivity (K _s), SWCC, particle composition, field capacity and some other indexes were determined. Results showed that the clay content and K _s of grassland soil were higher for GE15 than GE5. For both treatments, in low suction condition (≤100 kPa), the water holding capacity of 0–10 cm soil was the best. Water holding capacity of topsoil decreased gradually with the increasing of suction, and it reached the strongest when the suction reached 600 kPa. In all soil water suction, the water holding capacity of subsoil was the weakest. The van Genuchten expression was applicable for most of the samples, except 20–30 cm of GE5 and 10–20 cm of GE15. Dual porosity equation was applicable for all the samples. The soil–water capability and soil structure of which was fenced for 15 years is superior to that of 5 years. This study suggests that the enclosure management improved the soil structure and soil–water capability.     

Hydraulic,textural and geochemical characteristics of the Ajali Formation,Anambra Basin,Nigeria: implication for groundwater quality

This study highlights the distribution of hydraulic conductivity (K) in the regional aquiferous Ajali Formation of SE-Nigeria on one hand and assesses the possible influences of textural and geochemical characteristics on the hydraulic conductivity on the other hand. The investigation approach involved field sampling and collection of 12 sandstone samples from different outcrop locations, followed by laboratory studies such as grain-size analysis (GSA), constant head permeameter test and geochemical analysis of major and trace elements using X-ray fluorescence method. GSA and textural studies show that the sandstones range from fine to medium sands, constituting about <75–99% sand fraction, with graphic mean grain size of 0.23–0.53 mm. Other parameters such as coefficient of uniformity (C_u) range from 1.58 to 5.25 (av. 2.75), while standard deviation (sorting) values of 0.56Ø–1.24Ø imply moderately well sorted materials. In addition, the order of the estimated K values is K_permeameter>K_Beyer>K_Hazen>K_{Kozeny-Carmen}>K_Fair-Hatch with average values of 1.4×10^?3, 4.4×10^?4, 3.8×10^?4, 2.2×10^?4 and 8.1×10^?5m/s, respectively. These values fall within the range of 10^?5 and 10^?3m/s for fine to medium sands. However, multivariate factor analysis of the data revealed significant positive dependence of the empirically determined K values on graphic mean grain size and percentage sand content and much less dependence on sorting and total porosity. Geochemical profiles of the fresh samples are dominated by quartz with corresponding SiO₂ content of 76.1–98.2% (av. 89.7%) while other major oxides are generally below 1.0wt.% in the fresh samples. However, the ferruginized samples exhibited elevated concentrations of Al₂O₃ (3.50–11.60wt.%) and Fe₂O₃ (1.80–3.60wt.%), which are clear indications of weathering/ferruginization processes with attendant trace metal release/enrichment (2.5mg/l Cu, 7.5mg/l Pb, 6.5mg/l Zn, 3.9mg/l Ni and 19.6mg/l Cr) call for concern in respect of the chemical quality of the groundwater system. The associated groundwater is generally soft, slightly acidic, and with low dissolved solids (EC=14–134μs/cm) dominated by silica; implying water from clean sandy aquifer devoid of labile and weatherable minerals. Nonetheless, most of the metals (with exception of Si, Fe and Mn) exhibited higher degree of mobility (2–12 folds), which can be attributed to reduction of Fe-/Mn-oxyhydroxides as sinks/hosts for trace metals. Consequently, infiltration-induced geochemical reactions (redox, ferruginization and leaching processes) signify potential environmental impact in terms of groundwater quality as well as borehole/aquifer management, especially under humid tropical environment of the study area.     

Hydraulic conductivity prediction based on grain-size distribution using M5 model tree

The hydraulic conductivity, K_s, is one of the most important hydraulic properties which controls the water and solute movement into the soil. It is measured on soil specimens in the laboratory. On the other hand, sometimes it is obtained by tests carried out in the field by a number of researchers. Therefore, several experimental formulas have developed to predict it. Recently, soft computing tools have been used to evaluate the hydraulic conductivity. However, these tools are not as transparent as empirical formulas. In this study, another soft computing approach, i.e. model trees, have been used for predicting the hydraulic conductivity. The main advantage of model trees is that, unlike the other data learning tools, they are easier to use and represent understandable mathematical rules more clearly. In this paper, a new formula that includes some parameters is derived to estimate the hydraulic conductivity. To develop the new formulas, experimental data sets of hydraulic conductivity were used. A comparison is made between the estimated hydraulic conductivity by this new formula and formulas given by other’s researches.     

Soil hydraulic properties on the steep karst hillslopes in northwest Guangxi, China

Soil hydraulic properties such as soil infiltration rate and hydraulic conductivity are closely linked to runoff generation and infiltration processes but little is known about them on karst hillslopes. The objectives of this paper were to investigate the change in soil stable infiltration rate (q _s) and near-saturated hydraulic conductivity (K _ns) in different slope positions and to understand their relationship with rock fragment content and soil texture within the topsoil in subtropical karst regions of southwest China. Tension infiltrometers (20 cm in diameter) were used to measure q _s and K _ns at pressure head of −20 mm on hillslopes 1 (a disintegrated landslide failure) and 2 (an avalanche slope). The change of q _s and K _ns was great and they mostly had a moderate variability with coefficient of variations (CV) between 0.1 and 1.0 in the different slope positions. On average, q _s ranged from 0.43 to 4.25 mm/min and K _ns varied from 0.75 to 11.00 mm/min. These rates exceed those of most natural rainfall events, confirming that overland flow is rare on karst hillslopes. From bottom to top, q _s and K _ns had a decrease–increase–decrease trend due to the presence of large rock outcrops (>2 m in height) on hillslope 1 but had an increasing trend on hillslope 2 with less complex landform. They tended to increase with increase in total rock fragment content (5–250 mm) within the topsoil as well as slope gradient on both hillslopes. Pearson correlation analysis suggested that higher coarse pebble (20–75 mm), cobble (75–250 mm), and sand (2–0.05 mm) contents as well as total rock fragment content could significantly facilitate water infiltration into soils, but higher clay (<0.002 mm) content could restrict water movement. This result indicated that rock fragment, sand, and clay contents may remarkably affect water flow in the topsoil layers, and should be considered in hydrological modeling on karst hillslopes in subtropical regions.