Examining the information content of time-lapse crosshole GPR data collected under different infiltration conditions to estimate unsaturated soil hydraulic properties

Time-lapse geophysical data acquired during transient hydrological experiments are being increasingly employed to estimate subsurface hydraulic properties at the field scale. In particular, crosshole ground-penetrating radar (GPR) data, collected while water infiltrates into the subsurface either by natural or artificial means, have been demonstrated in a number of studies to contain valuable information concerning the hydraulic properties of the unsaturated zone. Previous work in this domain has considered a variety of infiltration conditions and different amounts of time-lapse GPR data in the estimation procedure. However, the particular benefits and drawbacks of these different strategies as well as the impact of a variety of key and common assumptions remain unclear. Using a Bayesian Markov-chain-Monte-Carlo stochastic inversion methodology, we examine in this paper the information content of time-lapse zero-offset-profile (ZOP) GPR traveltime data, collected under three different infiltration conditions, for the estimation of van Genuchten-Mualem (VGM) parameters in a layered subsurface medium. Specifically, we systematically analyze synthetic and field GPR data acquired under natural loading and two rates of forced infiltration, and we consider the value of incorporating different amounts of time-lapse measurements into the estimation procedure. Our results confirm that, for all infiltration scenarios considered, the ZOP GPR traveltime data contain important information about subsurface hydraulic properties as a function of depth, with forced infiltration offering the greatest potential for VGM parameter refinement because of the higher stressing of the hydrological system. Considering greater amounts of time-lapse data in the inversion procedure is also found to help refine VGM parameter estimates. Quite importantly, however, inconsistencies observed in the field results point to the strong possibility that posterior uncertainties are being influenced by model structural errors, which in turn underlines the fundamental importance of a systematic analysis of such errors in future related studies.

Estimate of Groundwater Recharge Based on Water Balance in The Unsaturated Soil Zone

ABSTRACT Groundwater management depends on the knowledge on recharge rates and water fluxes within aquifers. The recharge is one of the water cycle components most difficult to estimate. As a result, despite the chosen method, the estimates are subject to uncertainties that can be identified by means of comparison with other approaches. In this study, groundwater recharge estimates based on the water balance in the unsaturated zone is assessed. Firstly, the approach is evaluated by comparing the results with those of another method. Then, the estimates are used as inputs in a transient groundwater flow model in order to assess how the water table would respond to the obtained recharges rates compared to measured levels. The results suggest a good performance of the adopted approach and, despite some inherent limitations, it has advantages over other methods since the data required are easier to obtain.

Aspects of Mechanical Behavior and Modeling of a Tropical Unsaturated Soil

The research studies the applicability of two elastoplastic models for the collapse prediction of the lateritic soil profile from Southeastern Brazil. These tropical soils have peculiar geotechnical behavior, due to their mineralogical composition and porous structure coming from intense process of formation. Two elastoplastic models were analyzed: the Barcelona Basic Model (BBM) and another one based on BBM, however developed for tropical soils. Oedometric tests with suction control were performed at three distinct depths of the soil profile. The BBM was not suitable for the upper layer of the soil profile, because BBM considers the compressible behavior of the soil in function of the reduction of the elastoplastic compressibility index with the increase of the matric suction. The model developed for tropical soils showed better suited to the compressible behavior of the soil profile, resulting in good prediction of the collapse potential, mainly by accepting increasing values of the elastoplastic compressibility index of the soil profile with the matric suction rise. © 2013 Springer Science+Business Media Dordrecht.

Modeling field load tests in lateritic unsaturated soil

Many models for unsaturated soil have been developed in the last years, accompanying the development of experimental techniques to deal with such soils. The benchmark of the models for unsaturated soil can be assigned to the Barcelona Basic Model (BBM) now incorporated in some codes such as the CODE_BRIGHT. Most of those models were validated considering limited laboratory test results and not much validation is available considering real field problems. This paper presents modeling results of field plate load tests performed under known suction on a lateritic unsaturated soil. The required input data were taken from laboratory tests performed under suction control. The modeling nicely reproduces field tests allowing appreciating the influence of soil suction on the stress-settlement curve. In addition, wetting induced or collapse settlements were calculated from field tests and were nicely duplicated by the numerical analysis performed.

Small strain elastic behaviour of unsaturated soil investigated by bender/extender element testing

The aim of this project was to investigate very small strain elastic behaviour of soils under unsaturated conditions, using bender/extender element (BEE) testing. The behaviour of soils at very small strains has been widely studied under saturated conditions, whereas much less work has been performed on very small strain behaviour under unsaturated conditions. A suction-controlled double wall triaxial apparatus for unsaturated soil testing was modified to incorporate three pairs of BEEs transmitting both shear and compression waves with vertical and horizontal directions of wave transmission and wave polarisation. Various different techniques for measuring wave travel time were investigated in both the time domain and the frequency domain and it was concluded that, at least for the current experimental testing programme, peak-to-first-peak in the time domain was the most reliable technique for determining wave travel time. An experimental test programme was performed on samples of compacted speswhite kaolin clay. Two different forms of compaction were employed (i.e. isotropic and anisotropic). Compacted kaolin soil samples were subjected to constant suction loading and unloading stages at three different values of suction, covering both unsaturated conditions (s= 50kPa and s= 300kPa) and saturated conditions (s=0). Loading and unloading stages were performed at three different values of stress ratio (η=0, η=1 and η=-1 ). In some tests a wetting-drying cycle was performed before or within the loading stage, with the wetting-drying cycles including both wetting-induced swelling and wetting-induced collapse compression. BEE tests were performed at regular intervals throughout all test stages, to measure shear wave velocity Vs and compression wave velocity Vp and hence to determine values of shear modulus G and constrained modulus M. The experimental test programme was designed to investigate how very small strain shear modulus G and constrained modulus M varied with unsaturated state variables, including how anisotropy of these parameters developed either with stress state (stress-induced anisotropy) or with previous straining (strain-induced anisotropy). A new expression has been proposed for the very small strain shear modulus G of an isotropic soil under saturated and unsaturated conditions. This expression relates the variation of G to only mean Bishop’s stress p* and specific volume v, and it converges to a well-established expression for saturated soils as degree of saturation approaches 1. The proposed expression for G is able to predict the variation of G under saturated and unsaturated conditions at least as well as existing expressions from the literature and it is considerably simpler (employing fewer state variables and fewer soil constants). In addition, unlike existing expressions from the literature, the values of soil constants in the proposed new expression can be determined from a saturated test. It appeared that, in the current project at least, any strain-induced anisotropy of very small strain elastic behaviour was relatively modest, with the possible exception of loading in triaxial extension. It was therefore difficult to draw any firm conclusion about evolution of strain-induced anisotropy and whether it depended upon the same aspects of soil fabric as evolution of anisotropy of large strain plastic behaviour. Stress-induced anisotropy of very small strain elastic behaviour was apparent in the experimental test programme. An attempt was made to extend the proposed expression for G to include the effect of stress-induced anisotropy. Interpretation of the experimental results indicated that the value of shear modulus was affected by the values of all three principal Bishop’s stresses (in the direction of wave transmission, the direction of wave polarisation and the third mutually perpendicular direction). However, prediction of stress-induced anisotropy was only partially successful, and it was concluded that the effect of Lode angle was also significant.

A software to calculate soil hydraulic conductivity in internal drainage experiments (SHC, Version 2.00)

A software for the calculation of unsaturated soil hydraulic conductivity K(theta) is presented for commonly used methods found in the literature, based on field experiments in which a soil profile is submitted to water infiltration followed by internal drainage. The software is available at: dourado@esalq.usp.br.

Fractal concepts in relation to soil water diffusivity

Fractal geometry would appear to offer promise for new insight on water transport in unsaturated soils, This study was conducted to evaluate possible fractal influence on soil water diffusivity, and/or the relationships from which it arises, for several different soils, Fractal manifestations, consisting of a time-dependent diffusion coefficient and anomalous diffusion arising out of fractional Brownian motion, along with the notion of space-filling curves were gleaned from the literature, It was found necessary to replace the classical Boltzmann variable and its time t(1/2) factor with the basic fractal power function and its t(n) factor, For distinctly unsaturated soil water content theta, exponent n was found to be less than 1/2, but it approached 1/2 as theta approached its sated value, This function n = n(theta), in giving rise to a time-dependent, anomalous soil water diffusivity D, was identified with the Hurst exponent H of fractal geometry, Also, n approaching 1/2 at high water content is a behavior that makes it possible to associate factal space filling with soil that approaches water saturation, Finally, based on the fractally interpreted n = n(theta), the coalescence of both D and 8 data is greatly improved when compared with the coalescence provided by the classical Boltzmann variable.

Evaluation of Electrical Resistivity in a Tropical Sandy Soil Compacted

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The Effect of Contact on the Filter Paper Method for Measuring Soil Suction

The filter paper method is one of the most commonly used and critiqued techniques for measuring soil suction. However, many aspects related to its use still require some clarification. The results of a comprehensive study on the effect of the contact between the soil grains and soil water and the filter paper are presented herein. We investigated the influence of the equilibration time, the texture of the porous material and the degree of contact, or lack thereof, between the soil grains and the filter paper using Miamian #42 and three different types of porous material. To enhance the difference between the total suction and the matrix suction, osmotic suction was induced by saturating the specimens with a sodium chloride solution.

Unconfined shear strength of compacted unsaturated plastic soils

A central goal in unsaturated soil mechanics research is to create a smooth transition between traditional soil mechanics approaches and an approach that is applicable to unsaturated soils. Undrained shear strength and the liquidity index of reconstituted or remoulded saturated soils are consistently correlated, which has been demonstrated by many studies. In the liquidity index range from 1 (at w(l)) to 0 (at w(p)), the shear strength ranges from approximately 2 kPa to 200 kPa. Similarly, for compacted soil, the shear strength at the plastic limit ranges from 150 kPa to 250 kPa. When compacted at their optimum water content, most soils have a suction that ranges from 20 kPa to 500 kPa; however, in the field, compacted materials are subjected to drying and wetting, which affect their initial suction and as a consequence their shear strength. Unconfined shear tests were performed on five compacted tropical soils and kaolin. Specimens were tested in the as-compacted condition, and also after undergoing drying or wetting. The test results and data from prior literature were examined, taking into account the roles of void ratio, suction, and relative water content. An interpretation of the phenomena that are involved in the development of the undrained shear strength of unsaturated soils in the contexts of soil water retention and Atterberg limits is presented, providing a practical view of the behaviour of compacted soil based on the concept of unsaturated soil. Finally, an empirical correlation is presented that relates the unsaturated state of compacted soils to the unconfined shear strength.

Transport of iodide in structured clay–loam soil under maize during irrigation experiments analyzed using HYDRUS model

Transport of radioactive iodide 131I− in a structured clay loam soil under maize in a final growing phase was monitored during five consecutive irrigation experiments under ponding. Each time, 27 mm of water were applied. The water of the second experiment was spiked with 200 MBq of 131I− tracer. Its activity was monitored as functions of depth and time with Geiger-Müller (G-M) detectors in 11 vertically installed access tubes. The aim of the study was to widen our current knowledge of water and solute transport in unsaturated soil under different agriculturally cultivated settings. It was supposed that the change in 131I− activity (or counting rate) is proportional to the change in soil water content. Rapid increase followed by a gradual decrease in 131I− activity occurred at all depths and was attributed to preferential flow. The iodide transport through structured soil profile was simulated by the HYDRUS 1D model. The model predicted relatively deep percolation of iodide within a short time, in a good agreement with the observed vertical iodide distribution in soil. We found that the top 30 cm of the soil profile is the most vulnerable layer in terms of water and solute movement, which is the same depth where the root structure of maize can extend.

Long-term flow rates and biomat zone hydrology in soil columns receiving septic tank effluent

Soil absorption systems (SAS) are used commonly to treat and disperse septic tank effluent (STE). SAS can hydraulically fail as a result of the low permeable biomat zone that develops on the infiltrative surface. The objectives of this experiment were to compare the hydraulic properties of biomats grown in soils of different textures, to investigate the long-term acceptance rates (LTAR) from prolonged application of STE, and to assess if soils were of major importance in determining LTAR. The STE was applied to repacked sand, Oxisol and Vertisol soil columns over a period of 16 months, at equivalent hydraulic loading rates of 50, 35 and 8 L/m(2)/d, respectively Infiltration rates, soil matric potentials, and biomat hydraulic properties were measured either directly from the soil columns or calculated using established soil physics theory. Biomats 1 to 2 cm thick developed in all soils columns with hydraulic resistances of 27 to 39 d. These biomats reduced a 4 order of magnitude variation in saturated hydraulic conductivity (K.) between the soils to a one order of magnitude variation in LTAR. A relationship between biomat resistance and organic loading rate was observed in all soils. Saturated hydraulic conductivity influenced the rate and extent of biomat development. However, once the biomat was established, the LTAR was governed by the resistance of the biomat and the sub-biomat soil unsaturated flow regime induced by the biomat. Results show that whilst initial soil K. is likely to be important in the establishment of the biomat zone in a trench, LTAR is determined by the biomat resistance and the unsaturated soil hydraulic conductivity, not the K, of a soil. The results call into question the commonly used approach of basing the LTAR, and ultimately trench length in SAS, on the initial K, of soils. (c) 2006 Elsevier Ltd. All rights reserved.

Alternative Analytical Expressions for the General van Genuchten-Mualem and van Genuchten-Burdine Hydraulic Conductivity Models

The van Genuchten expressions for the unsaturated soil hydraulic properties, first published in 1980, are used frequently in various vadose zone flow and transport applications assuming a specific relationship between the m and n soil hydraulic parameters. By comparison, probably because of the complexity of the hydraulic conductivity equations, the more general solutions with independent m and n values are rarely used. We expressed the general van Genuchten-Mualem and van Genuchten-Burdine hydraulic conductivity equations in terms of hypergeometric functions, which can be approximated by infinite series that converge rapidly for relatively large values of the van Genuchten-Mualem parameter n but only very slowly when n is close to one. Alternative equations were derived that provide very close approximations of the analytical results. The newly proposed equations allow the use of independent values of the parameters m and n in the soil water retention model of van Genuchten for subsequent prediction of the van Genuchten-Mualem and van Genuchten-Burdine hydraulic conductivity models, thus providing more flexibility in fitting experimental pressure-head-dependent water content, theta(h), and hydraulic conductivity, K(h), or K(theta) data.

Measuring Hydraulic Conductivity to Wilting Point Using Polymer Tensiometers in an Evaporation Experiment

The polymer tensiometer is a novel instrument to measure soil water pressure heads from saturation to permanent wilting conditions. We used tensiometers of this type in an experiment to determine the hydraulic properties of evaporating soil samples in the laboratory. Relative errors in the hydraulic conductivity function in the wet part were high due to the relatively low accuracy of the pressure transducers, resulting in a large uncertainty in the hydraulic gradient and therefore in the calculated hydraulic conductivity. In the dry part, the error related to this accuracy was on the same order of magnitude as the error related to balance accuracy. Therefore, the method can be assumed adequate for measuring soil hydraulic properties except under very wet conditions. In our experiments, relative error and bias increased significantly at pressure heads less negative than -1 m.