Modeling Transient Streaming Potentials in Falling‐Head Permeameter Tests

We present transient streaming potential data collected during falling‐head permeameter tests performed on samples of two sands with different physical and chemical properties. The objective of the work is to estimate hydraulic conductivity (K) and the electrokinetic coupling coefficient (C_l) of the sand samples. A semi‐empirical model based on the falling‐head permeameter flow model and electrokinetic coupling is used to analyze the streaming potential data and to estimate K and C_l. The values of K estimated from head data are used to validate the streaming potential method. Estimates of K from streaming potential data closely match those obtained from the associated head data, with less than 10% deviation. The electrokinetic coupling coefficient was estimated from streaming potential vs. (1) time and (2) head data for both sands. The results indicate that, within limits of experimental error, the values of C_l estimated by the two methods are essentially the same. The results of this work demonstrate that a temporal record of the streaming potential response in falling‐head permeameter tests can be used to estimate both K and C_l. They further indicate the potential for using transient streaming potential data as a proxy for hydraulic head in hydrogeology applications.     

Anisotropy and depth‐related heterogeneity of hydraulic conductivity in a bog peat. I: laboratory measurements

Anisotropy and heterogeneity of hydraulic conductivity (K) are suspected of greatly affecting rates and patterns of ground‐water seepage in peats. A new laboratory method, termed here the modified cube method, was used to measure horizontal and vertical hydraulic conductivity (K_h and K_v) of 400 samples of bog peat. The new method avoids many of the problems associated with existing field and laboratory methods, and is shown to give relatively precise measurements of K. In the majority of samples tested, K_h was much greater than K_v, indicating that the bog peat was strongly anisotropic. Log₁₀K_h, log₁₀K_v, and log₁₀ (K_h/K_v) were found to vary significantly with depth, although none of the relationships was simple. We comment on the scale dependency of our measurements. Copyright © 2002 John Wiley & Sons, Ltd.     

A Simple Field Method for Assessing Near‐Surface Saturated Hydraulic Conductivity

This technical note describes an effective and inexpensive field technique for measuring the saturated hydraulic conductivity of both undisturbed cores and repacked soil samples. The method requires no specialized equipment; everything that is required can be obtained in a hardware store. The method is a straightforward field implementation of the widely used falling‐head laboratory analysis directly derived from Darcy's law. As such, it sidesteps the need for empirical assumptions about soil texture and the relationship between saturated and unsaturated flow components which many permeameter‐based methods rely upon. The method is shown to produce results that are consistent with K values obtained elsewhere in the same homogeneous sand formation. Furthermore, the proposed method is useful for measuring hydraulic conductivity in drill cuttings obtained from direct push or auguring drill techniques, which cannot be done with any other field method. The range of hydraulic conductivity values that this test is appropriate for is on the order of 1E ? 7 m/s to 1E ? 3 m/s.     

Horizontal hydraulic conductivity of shallow streambed sediments and comparison with the grain‐size analysis results

Streambed horizontal hydraulic conductivity (K_h) has a substantial role in controlling exchanges between stream water and groundwater. We propose a new approach for determining K_h of the shallow streambed sediments. Undisturbed sediment samples were collected using tubes that were horizontally driven into streambeds. The sediment columns were analysed using a permeameter test (PT) on site. This new test approach minimizes uncertainties due to vertical flow in the vicinity of test tube and stream stage fluctuations in the computation of the K_h values. Ninety‐eight PTs using the new approach were conducted at eight sites in four tributaries of the Platte River, east‐central Nebraska, USA. The K_h values were compared with the nondirectional hydraulic conductivity values (K_g) determined from 12 empirical grain‐size analysis methods. The grain‐size analysis methods used the same sediment samples as K_h tests. Only two methods, the Terzaghi and Shepherd methods, yielded K_g values close to the K_h values. Although the Sauerbrei method produced a value relatively closer to K_h than other nine grain‐size analysis methods, the values from this method were not as reliable as the methods of Terzaghi and Shepherd due to the inconsistent fluctuation of the average estimates at each of the test sites. The Zunker, Zamarin, Hazen, Beyer, and Kozeny methods overestimated K_h, while the Slichter, US Bureau of Reclamation (USBR), Harleman, and Alyamani and Sen methods underestimated K_h. Any of these specific grain‐size methods might yield good estimates of streambed K_h at some sites, but give poor estimates at other sites, indicating that the relationship between K_g and K_h is significantly site dependent in our study. Copyright © 2011 John Wiley & Sons, Ltd.     

Depth‐dependent hydraulic conductivity distribution patterns of a streambed

Characterization of streambed hydraulic conductivity from the channel surface to a great depth below the channel surface can provide needed information for the determination of stream‐aquifer hydrologic connectedness, and it is also important to river restoration. However, knowledge on the streambed hydraulic conductivity for sediments 1 m below the channel surface is scarce. This study describes a method that was used to determine the distribution patterns of streambed hydraulic conductivity for sediments from channel surface to a depth of 15 m below. The method includes Geoprobe's direct‐push techniques and Permeameter tests. Direct‐push techniques were used to generate the electrical conductivity (EC) logs and to collect sequences of continuous sediment cores from river channels, as well as from the alluvial aquifer connected to the river. Permeameter tests on these sediment cores give the profiles of vertical hydraulic conductivity (K_v) of the channel sediments and the aquifer materials. This method was applied to produce K_v profiles for a streambed and an alluvial aquifer in the Platte River Valley of Nebraska, USA. Comparison and statistical analysis of the K_v profiles from the river channel and from the proximate alluvial aquifer indicates a special pattern of K_v in the channel sediments. This depth‐dependent pattern of K_v distribution for the channel sediments is considered to be produced by hyporheic processes. This K_v‐distribution pattern implied that the effect of hyporheic processes on streambed hydraulic conductivity can reach the sediments about 9 m below the channel surface. Copyright © 2010 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.     

A field‐scale infiltration model accounting for spatial heterogeneity of rainfall and soil saturated hydraulic conductivity

This study first explores the role of spatial heterogeneity, in both the saturated hydraulic conductivity K_s and rainfall intensity r, on the integrated hydrological response of a natural slope. On this basis, a mathematical model for estimating the expected areal‐average infiltration is then formulated. Both K_s and r are considered as random variables with assessed probability density functions. The model relies upon a semi‐analytical component, which describes the directly infiltrated rainfall, and an empirical component, which accounts further for the infiltration of surface water running downslope into pervious soils (the run‐on effect). Monte Carlo simulations over a clay loam soil and a sandy loam soil were performed for constructing the ensemble averages of field‐scale infiltration used for model validation. The model produced very accurate estimates of the expected field‐scale infiltration rate, as well as of the outflow generated by significant rainfall events. Furthermore, the two model components were found to interact appropriately for different weights of the two infiltration mechanisms involved. Copyright © 2005 John Wiley & Sons, Ltd.     

Determining hydraulic properties of a loam soil by alternative infiltrometer techniques

Testing infiltrometer techniques to determine soil hydraulic properties is necessary for specific soils. For a loam soil, the water retention and hydraulic conductivity predicted by the BEST (Beerkan Estimation of Soil Transfer parameters) procedure of soil hydraulic characterization was compared with data collected by more standard laboratory and field techniques. Six infiltrometer techniques were also compared in terms of saturated soil hydraulic conductivity, K_s. BEST yielded water retention values statistically similar to those obtained in the laboratory and K_s values practically coinciding with those determined in the field with the pressure infiltrometer (PI). The unsaturated soil hydraulic conductivity measured with the tension infiltrometer (TI) was reproduced satisfactorily by BEST only close to saturation. BEST, the PI, one‐potential experiments with both the TI and the mini disk infiltrometer (MDI), the simplified falling head (SFH) technique and the bottomless bucket (BB) method yielded statistically similar estimates of K_s, differing at the most by a factor of three. Smaller values were obtained with longer and more soil‐disturbing infiltration runs. Any of the tested infiltration techniques appears usable to obtain the order of magnitude of K_s at the field site, but the BEST, BB and PI data appear more appropriate to characterize the soil at some stage during a rainfall event. Additional investigations on both similar and different soils would allow development of more general procedures to apply infiltrometer techniques for soil hydraulic characterization. Copyright © 2015 John Wiley & Sons, Ltd.     

Transient Recharge Estimability Through Field‐Scale Groundwater Model Calibration

The estimation of recharge through groundwater model calibration is hampered by the nonuniqueness of recharge and aquifer parameter values. It has been shown recently that the estimability of spatially distributed recharge through calibration of steady‐state models for practical situations (i.e., real‐world, field‐scale aquifer settings) is limited by the need for excessive amounts of hydraulic‐parameter and groundwater‐level data. However, the extent to which temporal recharge variability can be informed through transient model calibration, which involves larger water‐level datasets, but requires the additional consideration of storage parameters, is presently unknown for practical situations. In this study, time‐varying recharge estimates, inferred through calibration of a field‐scale highly parameterized groundwater model, are systematically investigated subject to changes in (1) the degree to which hydraulic parameters including hydraulic conductivity (K) and specific yield (S_y) are constrained, (2) the number of water‐level calibration targets, and (3) the temporal resolution (up to monthly time steps) at which recharge is estimated. The analysis involves the use of a synthetic reality (a reference model) based on a groundwater model of Uley South Basin, South Australia. Identifiability statistics are used to evaluate the ability of recharge and hydraulic parameters to be estimated uniquely. Results show that reasonable estimates of monthly recharge (<30% recharge root‐mean‐squared error) require a considerable amount of transient water‐level data, and that the spatial distribution of K is known. Joint estimation of recharge, S_y and K, however, precludes reasonable inference of recharge and hydraulic parameter values. We conclude that the estimation of temporal recharge variability through calibration may be impractical for real‐world settings.     

Field Study of Subsurface Heterogeneity with Steady‐State Hydraulic Tomography

Remediation of subsurface contamination requires an understanding of the contaminant (history, source location, plume extent and concentration, etc.), and, knowledge of the spatial distribution of hydraulic conductivity (K) that governs groundwater flow and solute transport. Many methods exist for characterizing K heterogeneity, but most if not all methods require the collection of a large number of small‐scale data and its interpolation. In this study, we conduct a hydraulic tomography survey at a highly heterogeneous glaciofluvial deposit at the North Campus Research Site (NCRS) located at the University of Waterloo, Waterloo, Ontario, Canada to sequentially interpret four pumping tests using the steady‐state form of the Sequential Successive Linear Estimator (SSLE) ( Yeh and Liu 2000 ). The resulting three‐dimensional (3D) K distribution (or K‐tomogram) is compared against: ( 1 ) K distributions obtained through the inverse modeling of individual pumping tests using SSLE, and ( 2 ) effective hydraulic conductivity (K_eff) estimates obtained by automatically calibrating a groundwater flow model while treating the medium to be homogeneous. Such a K_eff is often used for designing remediation operations, and thus is used as the basis for comparison with the K‐tomogram. Our results clearly show that hydraulic tomography is superior to the inversions of single pumping tests or K_eff estimates. This is particularly significant for contaminated sites where an accurate representation of the flow field is critical for simulating contaminant transport and injection of chemical and biological agents used for active remediation of contaminant source zones and plumes.     

Synthesis of soil‐hydraulic properties and infiltration timescales in wildfire‐affected soils

We collected soil‐hydraulic property data from the literature for wildfire‐affected soils, ash, and unburned soils. These data were used to calculate metrics and timescales of hydrologic response related to infiltration and surface runoff generation. Sorptivity (S) and wetting front potential (Ψ_f) were significantly different (lower) in burned soils compared with unburned soils, whereas field‐saturated hydraulic conductivity (K_fs) was not significantly different. The magnitude and duration of the influence of capillarity during infiltration was greatly reduced in burned soils, causing faster ponding times in response to rainfall. Ash had large values of S and K_fs but moderate values of Ψ_f, compared with unburned and burned soils, indicating ash has long ponding times in response to rainfall. The ratio of S²/K_fs was nearly constant (~100 mm) for unburned soils but more variable in burned soils, suggesting that unburned soils have a balance between gravity and capillarity contributions to infiltration that may depend on soil organic matter, whereas in burned soils the gravity contribution to infiltration is greater. Changes in S and K_fs in burned soils act synergistically to reduce infiltration and accelerate and amplify surface runoff generation. Synthesis of these findings identifies three key areas for future research. First, short timescales of capillary influences on infiltration indicate the need for better measurements of infiltration at times less than 1 min to accurately characterize S in burned soils. Second, using parameter values, such as Ψ_f, from unburned areas could produce substantial errors in hydrologic modeling when used without adjustment for wildfire effects, causing parameter compensation and resulting underestimation of K_fs. Third, more thorough measurement campaigns that capture soil‐structural changes, organic matter impacts, quantitative water repellency trends, and soil‐water content along with soil‐hydraulic properties could drive the development of better techniques for numerically simulating infiltration in burned areas.     

Hydraulic properties of peat soils along a bulk density gradient—A meta study

Our understanding of hydraulic properties of peat soils is limited compared with that of mineral substrates. In this study, we aimed to deduce possible alterations of hydraulic properties of peat soils following degradation resulting from peat drainage and aeration. A data set of peat hydraulic properties (188 soil water retention curves [SWRCs], 71 unsaturated hydraulic conductivity curves [UHCs], and 256 saturated hydraulic conductivity [K_s] values) was assembled from the literature; the obtained data originated from peat samples with an organic matter (OM) content ranging from 23 to 97 wt% (weight percent; and according variation in bulk density) representing various degrees of peat degradation. The Mualem‐van Genuchten model was employed to describe the SWRCs and UHCs. The results show that the hydraulic parameters of peat soils vary over a wide range confirming the pronounced diversity of peat. Peat decomposition significantly modifies all hydraulic parameters. A bulk density of approximately 0.2 g cm^?3 was identified as a critical threshold point; above and below this value, macroporosity and hydraulic parameters follow different functions with bulk density. Pedotransfer functions based on physical peat properties (e.g., bulk density and soil depth) separately computed for bog and fen peat have significantly lower mean square errors than functions obtained from the complete data set, which indicates that not only the status of peat decomposition but also the peat‐forming plants have a large effect on hydraulic properties. The SWRCs of samples with a bulk density of less than 0.2 g cm^?3 could be grouped into two to five classes for each peat type (botanical composition). The remaining SWRCs originating from samples with a bulk density of larger than 0.2 g cm^?3 could be classified into one group. The Mualem‐van Genuchten parameter values of α can be used to estimate K_s if no K_s data are available. In conclusion, the derived pedotransfer functions provide a solid instrument to derive hydraulic parameter values from easily measurable quantities; however, additional research is required to reduce uncertainty.     

A Tube Seepage Meter for In Situ Measurement of Seepage Rate and Groundwater Sampling

We designed and evaluated a “tube seepage meter” for point measurements of vertical seepage rates (q), collecting groundwater samples, and estimating vertical hydraulic conductivity (K) in streambeds. Laboratory testing in artificial streambeds show that seepage rates from the tube seepage meter agreed well with expected values. Results of field testing of the tube seepage meter in a sandy‐bottom stream with a mean seepage rate of about 0.5 m/day agreed well with Darcian estimates (vertical hydraulic conductivity times head gradient) when averaged over multiple measurements. The uncertainties in q and K were evaluated with a Monte Carlo method and are typically 20% and 60%, respectively, for field data, and depend on the magnitude of the hydraulic gradient and the uncertainty in head measurements. The primary advantages of the tube seepage meter are its small footprint, concurrent and colocated assessments of q and K, and that it can also be configured as a self‐purging groundwater‐sampling device.     

Infiltration and runoff generation processes in fire‐affected soils

Post‐wildfire runoff was investigated by combining field measurements and modelling of infiltration into fire‐affected soils to predict time‐to‐start of runoff and peak runoff rate at the plot scale (1 m²). Time series of soil‐water content, rainfall and runoff were measured on a hillslope burned by the 2010 Fourmile Canyon Fire west of Boulder, Colorado during cyclonic and convective rainstorms in the spring and summer of 2011. Some of the field measurements and measured soil physical properties were used to calibrate a one‐dimensional post‐wildfire numerical model, which was then used as a ‘virtual instrument’ to provide estimates of the saturated hydraulic conductivity and high‐resolution (1 mm) estimates of the soil‐water profile and water fluxes within the unsaturated zone. Field and model estimates of the wetting‐front depth indicated that post‐wildfire infiltration was on average confined to shallow depths less than 30 mm. Model estimates of the effective saturated hydraulic conductivity, K_s, near the soil surface ranged from 0.1 to 5.2 mm h^?1. Because of the relatively small values of K_s, the time‐to‐start of runoff (measured from the start of rainfall), t_p, was found to depend only on the initial soil‐water saturation deficit (predicted by the model) and a measured characteristic of the rainfall profile (referred to as the average rainfall acceleration, equal to the initial rate of change in rainfall intensity). An analytical model was developed from the combined results and explained 92–97% of the variance of t_p, and the numerical infiltration model explained 74–91% of the variance of the peak runoff rates. These results are from one burned site, but they strongly suggest that t_p in fire‐affected soils (which often have low values of K_s) is probably controlled more by the storm profile and the initial soil‐water saturation deficit than by soil hydraulic properties. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Screening‐Level Sensitivity Analysis for the Design of Pump‐and‐Treat Systems

Pump‐and‐treat systems can prevent the migration of groundwater contaminants and candidate systems are typically evaluated with groundwater models. Such models should be rigorously assessed to determine predictive capabilities and numerous tools and techniques for model assessment are available. While various assessment methodologies (e.g., model calibration, uncertainty analysis, and Bayesian inference) are well‐established for groundwater modeling, this paper calls attention to an alternative assessment technique known as screening‐level sensitivity analysis (SLSA). SLSA can quickly quantify first‐order (i.e., main effects) measures of parameter influence in connection with various model outputs. Subsequent comparisons of parameter influence with respect to calibration vs. prediction outputs can suggest gaps in model structure and/or data. Thus, while SLSA has received little attention in the context of groundwater modeling and remedial system design, it can nonetheless serve as a useful and computationally efficient tool for preliminary model assessment. To illustrate the use of SLSA in the context of designing groundwater remediation systems, four SLSA techniques were applied to a hypothetical, yet realistic, pump‐and‐treat case study to determine the relative influence of six hydraulic conductivity parameters. Considered methods were: Taguchi design‐of‐experiments (TDOE); Monte Carlo statistical independence (MCSI) tests; average composite scaled sensitivities (ACSS); and elementary effects sensitivity analysis (EESA). In terms of performance, the various methods identified the same parameters as being the most influential for a given simulation output. Furthermore, results indicate that the background hydraulic conductivity is important for predicting system performance, but calibration outputs are insensitive to this parameter (K_BK). The observed insensitivity is attributed to a nonphysical specified‐head boundary condition used in the model formulation which effectively “staples” head values located within the conductivity zone. Thus, potential strategies for improving model predictive capabilities include additional data collection targeting the K_BK parameter and/or revision of model structure to reduce the influence of the specified head boundary.     

Meta‐analysis of field‐saturated hydraulic conductivity recovery following wildland fire: Applications for hydrologic model parameterization and resilience assessment

Hydrologic recovery after wildfire is critical for restoring the ecosystem services of protecting of human lives and infrastructure from hazards and delivering water supply of sufficient quality and quantity. Recovery of soil‐hydraulic properties, such as field‐saturated hydraulic conductivity (K_fs), is a key factor for assessing the duration of watershed‐scale flash flood and debris flow risks after wildfire. Despite the crucial role of K_fs in parameterizing numerical hydrologic models to predict the magnitude of postwildfire run‐off and erosion, existing quantitative relations to predict K_fs recovery with time since wildfire are lacking. Here, we conduct meta‐analyses of 5 datasets from the literature that measure or estimate K_fs with time since wildfire for longer than 3‐year duration. The meta‐analyses focus on fitting 2 quantitative relations (linear and non‐linear logistic) to explain trends in K_fs temporal recovery. The 2 relations adequately described temporal recovery except for 1 site where macropore flow dominated infiltration and K_fs recovery. This work also suggests that K_fs can have low hydrologic resistance (large postfire changes), and moderate to high hydrologic stability (recovery time relative to disturbance recurrence interval) and resilience (recovery of hydrologic function and provision of ecosystem services). Future K_fs relations could more explicitly incorporate processes such as soil‐water repellency, ground cover and soil structure regeneration, macropore recovery, and vegetation regrowth.     

Estimating field‐saturated soil hydraulic conductivity by a simplified Beerkan infiltration experiment

Field‐saturated soil hydraulic conductivity, K_fs, is highly variable. Therefore, interpreting and simulating hydrological processes, such as rainfall excess generation, need a large number of K_fs data even at the plot scale. Simple and reasonably rapid experiments should be carried out in the field. In this investigation, a simple infiltration experiment with a ring inserted shortly into the soil and the estimation of the so‐called α* parameter allowed to obtain an approximate measurement of K_fs. The theoretical approach was tested with reference to 149 sampling points established on Burundian soils. The estimated K_fs with the value of first approximation of α* for most agricultural field soils (α* = 0.012 mm^?1) differed by a practically negligible maximum factor of two from the saturated conductivity obtained by the complete Beerkan Estimation of Soil Transfer parameters (BEST) procedure for soil hydraulic characterization. The measured infiltration curve contained the necessary information to obtain a site‐specific prediction of α*. The empirically derived α* relationship gave similar results for K_fs (mean = 0.085 mm s^?1; coefficient of variation (CV) = 71%) to those obtained with BEST (mean = 0.086 mm s^?1; CV = 67%), and it was also successfully tested with reference to a few Sicilian sampling points, since it yielded a mean and a CV of K_fs (0.0094 mm s^?1 and 102%, respectively) close to the values obtained with BEST (mean = 0.0092 mm s^?1; CV = 113%). The developed method appears attractive due to the extreme simplicity of the experiment. Copyright © 2012 John Wiley & Sons, Ltd.     

Incorporating spatially heterogeneous infiltration capacity into hydrologic models with applications for simulating post‐wildfire debris flow initiation

Soils in post‐wildfire environments are often characterized by a low infiltration capacity with a high degree of spatial heterogeneity relative to unburned areas. Debris flows are frequently initiated by run‐off in recently burned steeplands, making it critical to develop and test methods for incorporating spatial variability in infiltration capacity into hydrologic models. We use Monte Carlo simulations of run‐off generation over a soil with a spatially heterogenous saturated hydraulic conductivity (K_s) to derive an expression for an aerially averaged saturated hydraulic conductivity ( ) that depends on the rainfall rate, the statistical properties of K_s, and the spatial correlation length scale associated with K_s. The proposed method for determining is tested by simulating run‐off on synthetic topography over a wide range of spatial scales. Results provide a simplified expression for an effective saturated hydraulic conductivity that can be used to relate a distribution of small‐scale K_s measurements to infiltration and run‐off generation over larger spatial scales. Finally, we use a hydrologic model based on to simulate run‐off and debris flow initiation at a recently burned catchment in the Santa Ana Mountains, CA, USA, and compare results to those obtained using an infiltration model based on the Soil Conservation Service Curve Number.     

Comparing hydraulic properties of different forest floors

The forest floor plays an important role in runoff rate, soil erosion and soil infiltration capacity by protecting mineral soils from the direct impact of falling raindrops. Forest floor consists of different kinds of litter with different hydraulic properties. In this study, the inverse method was used to estimate the hydraulic properties of three kinds of forest floor (broad‐leaved, needle‐leaved and mixed‐stand) at three replications in a completely random design. Forest floor samples were collected from the Gilan Province, Iran. The samples were piled up to make long columns 40.88 cm high with an inner diameter of 18.1 cm. Artificial rainfall experiments were conducted on top of the columns, and free drainage from the bottom of the columns was measured in the laboratory. Saturated hydraulic conductivity (K_s), saturated water content and water retention curve parameters (van Genuchten equation) were estimated by the inverse method. The results showed that the K_s of needle‐leaved samples differed significantly (p < 0.05) from those of broad‐leaved and mixed‐stand samples, whereas the latter two did not differ in this regard. No significant differences emerged in the water retention function parameters of van Genuchten (θ_r, β and α) in the three forest floor samples. The saturated water content of mixed‐stand samples was significantly different (p < 0.05) from that of broad‐leaved and needle‐leaved treatments with the latter two samples showing no significant difference. The good agreement between simulated and observed free drainage for all forest floor samples in the validation period indicates that the estimated hydraulic properties efficiently characterize the unsaturated water flow in the forest floor. Copyright © 2013 John Wiley & Sons, Ltd.     

Sorption Behaviour of Phenols on Natural Sandy Aquifer Material during Flow‐through Column Experiments: The Effect of pH

Flow‐through column experiments were carried out to investigate the influence of pH on the sorption of three phenols (2‐methyl‐4, 6‐dinitrophenol, 2, 4, 6‐trichlorophenol, pentachlorophenol) onto a natural sandy aquifer material collected from a bank filtration site of River Elbe, Germany. For the phenols investigated, an increase in sorption (retardation) with decreasing pH is observed indicating a stronger sorption of the neutral species in comparison to that of the anions formed by dissociation. The anions of 2‐methyl‐4, 6‐dinitrophenol and 2, 4, 6‐trichlorophenol do not show significant sorption. On the contrary, pentachlorophenol showed sorption not only in neutral form but also in ionic form significantly which should be taken into account while assessing the fate and transport of such compound. A linear model based on the degree of protonation (calculated from pH and pK_a) can be used to resolve the apparent (observed) sorption coefficient (K_{d, app}) into its neutral (K_{d, n}) and ionised (K_{d, i}) components. Knowing pK_a, K_{d, n}, and K_{d, i} the apparent sorption coefficient for pH values other than experimentally investigated can be predicted.