Role of degree of saturation in denitrification of unsaturated sand specimens

Nitrate contamination of groundwater arises from anthropogenic activities, such as, fertilizer and animal manure applications and infiltration of wastewater/leachates. During migration of wastewater and leachates, the vadose zone (zone residing above the groundwater table), is considered to facilitate microbial denitrification. Particle voids in vadose zone are deficient in dissolved oxygen as the voids are partially filled by water and the remainder by air. Discontinuities in liquid phase would also restrict oxygen diffusion and therefore facilitate denitrification in the vadose/unsaturated soil zone. The degree of saturation of soil specimen (S (r)) quantifies the relative volume of voids filled with air and water. Unsaturated specimens have S (r) values ranging between 0 and 100 %. Earlier studies from naturally occurring nitrate losses in groundwater aquifers in Mulbagal town, Kolar District, Karnataka, showed that the sub-surface soils composed of residually derived sandy soil; hence, natural sand was chosen in the laboratory denitrification experiments. With a view to understand the role of vadose zone in denitrification process, experiments are performed with unsaturated sand specimens (S (r) = 73-90 %) whose pore water was spiked with nitrate and ethanol solutions. Experimental results revealed 73 % S (r) specimen facilitates nitrate reduction to 45 mg/L in relatively short durations of 5.5-7.5 h using the available natural organic matter (0.41 % on mass basis of sand); consequently, ethanol addition did not impact rate of denitrification. However, at higher S (r) values of 81 and 90 %, extraneous ethanol addition (C/N = 0.5-3) was needed to accelerate the denitrification rates.

Mechanical behaviour of unsaturated kaolin (with isotropic and anisotropic stress history). Part 2: performance under shear loading

Validation of a framework for unsaturated soil behaviour has frequently resulted in disagreement with basic propositions. A primary reason for this disparity is considered to be attributable to the anisotropic properties of the soil specimens tested as a result of preparation using one-dimensional compaction. As part of the work presented, comparison is made between tests on samples of unsaturated kaolin prepared at identical specific volumes and specific water volumes using isotropic compression and one-dimensional compression. The suctions in the samples were reduced to predefined values by wetting under low isotropic loading in a triaxial cell. The samples were then taken through various stress paths to failure, defined as the critical state strength, while the suctions were held constant. Stress path tests were also performed on samples without reducing the suction to predefined values. In the latter, constant water mass tests, the suctions were allowed to vary and were measured using a psychrometer. The results of the tests at critical state are compared with the propositions of Wheeler and Sivakumar. The shear strengths of samples with isotropic previous history are shown to be significantly greater than those of samples with one-dimensional stress history when plotted against the mean net stress. The normal compression lines, critical state lines and yield characteristics are also shown to be significantly influenced by the previous stress history and are shown to be different for isotropically and one-dimensionally prepared samples.

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.

Investigation on parameters used in warning systems for rain-induced embankment instability

A number of instrumented laboratory-scale soil embankment slopes were subjected to artificial rainfall until they failed. The factor of safety of the slope based on real-time measurements of pore-water pressure (suction) and laboratory measured soil properties were calculated as the rainfall progressed. Based on the experiment measurements and slope stability analysis, it was observed that slope displacement measurements can be used to warn the slope failure more accurately. Further, moisture content/pore-water pressure measurements near the toe of the slope and the real-time factor of safety can also be used for prediction of rainfall-induced embankment failures with adequate accuracy.

The effects of deep cracks on the rain-induced instability of slopes : A case study

Rainfall can disrupt the balance of natural soil slope. This imbalance will be accelerated by existence of cracks in soil slope, which lead to decreasing shear strength and increasing hydraulic conductivity of the soil slope. Some research works have been conducted on the effects of surface-cracks on slope stability. However, the influence of deep-cracks is yet to be investigated. Limited availability of deep crack data due to the lack of effective sub-soil investigation methods could be one of the obstacles. To emphasize the effects of deep cracks in soil slope on its rain-induced instability, a natural soil slope in Indonesia that failed in 31st October 2010 due to heavy rainfall was analyzed for stability with and without deep cracks in the slope. The slope stability analysis was conducted using SLOPE/W coupling with the results of transient seepage analysis (SEEP/W) that simulate the pore-water pressure development in the slope during the rainfall. The results of Electrical Resistivity Tomography (ERT) survey, bore-hole tests and geometrical survey conducted on the slope before its failure were used to identify the soil layers’ stratification including deep cracks, the properties of different soil layers, and geometrical parameters of the slope for the analysis. The results showed that it is vital to consider the existence of deep crack in soil slopes in analysing their instability induced by rainfalls.

Country-scale carbon accounting of the vegetation and mineral soils of Finland

The increase in global temperature has been attributed to increased atmospheric concentrations of greenhouse gases (GHG), mainly that of CO2. The threat of severe and complex socio-economic and ecological implications of climate change have initiated an international process that aims to reduce emissions, to increase C sinks, and to protect existing C reservoirs. The famous Kyoto protocol is an offspring of this process. The Kyoto protocol and its accords state that signatory countries need to monitor their forest C pools, and to follow the guidelines set by the IPCC in the preparation, reporting and quality assessment of the C pool change estimates. The aims of this thesis were i) to estimate the changes in carbon stocks vegetation and soil in the forests in Finnish forests from 1922 to 2004, ii) to evaluate the applied methodology by using empirical data, iii) to assess the reliability of the estimates by means of uncertainty analysis, iv) to assess the effect of forest C sinks on the reliability of the entire national GHG inventory, and finally, v) to present an application of model-based stratification to a large-scale sampling design of soil C stock changes. The applied methodology builds on the forest inventory measured data (or modelled stand data), and uses statistical modelling to predict biomasses and litter productions, as well as a dynamic soil C model to predict the decomposition of litter. The mean vegetation C sink of Finnish forests from 1922 to 2004 was 3.3 Tg C a-1, and in soil was 0.7 Tg C a-1. Soil is slowly accumulating C as a consequence of increased growing stock and unsaturated soil C stocks in relation to current detritus input to soil that is higher than in the beginning of the period. Annual estimates of vegetation and soil C stock changes fluctuated considerably during the period, were frequently opposite (e.g. vegetation was a sink but soil was a source). The inclusion of vegetation sinks into the national GHG inventory of 2003 increased its uncertainty from between -4% and 9% to ± 19% (95% CI), and further inclusion of upland mineral soils increased it to ± 24%. The uncertainties of annual sinks can be reduced most efficiently by concentrating on the quality of the model input data. Despite the decreased precision of the national GHG inventory, the inclusion of uncertain sinks improves its accuracy due to the larger sectoral coverage of the inventory. If the national soil sink estimates were prepared by repeated soil sampling of model-stratified sample plots, the uncertainties would be accounted for in the stratum formation and sample allocation. Otherwise, the increases of sampling efficiency by stratification remain smaller. The highly variable and frequently opposite annual changes in ecosystem C pools imply the importance of full ecosystem C accounting. If forest C sink estimates will be used in practice average sink estimates seem a more reasonable basis than the annual estimates. This is due to the fact that annual forest sinks vary considerably and annual estimates are uncertain, and they have severe consequences for the reliability of the total national GHG balance. The estimation of average sinks should still be based on annual or even more frequent data due to the non-linear decomposition process that is influenced by the annual climate. The methodology used in this study to predict forest C sinks can be transferred to other countries with some modifications. The ultimate verification of sink estimates should be based on comparison to empirical data, in which case the model-based stratification presented in this study can serve to improve the efficiency of the sampling design.

Numerical Analysis of Rainfall Infiltration in the Slope with a Fracture

With the finite volume method, a 2D numerical model for seepage in unsaturated soil has been established to study the rainfall infiltration in the fractured slope.The result shows that more rain may infiltrate into the slope due to existing fracture and then the pore pressure rises correspondingly. Very probably, it is one of the crucial factors accounting for slope failure. Furthermore a preliminary study has been conducted to investigate the influence of various fracture and rainfall factors such as the depth, width and location of a crack, surface condition, rainfall intensity and duration. Pore pressure and water volumetric content during the transient seepage are carefully examined to reveal the intrinsic mechanism.

Avaliação de proposições matemáticas para interpretação do comportamento de solos residuais não saturados

Nas últimas décadas, teorias têm sido formuladas para interpretar o comportamento de solos não saturados e estas têm se mostrado coerentes com resultados experimentais. Paralelamente, várias técnicas de campo e de laboratório têm sido desenvolvidas. No entanto, a determinação experimental dos parâmetros dos solos não saturados é cara, morosa, exige equipamentos especiais e técnicos experientes. Como resultado, essas teorias têm aplicação limitada a pesquisas acadêmicas e são pouco utilizados na prática da engenharia. Para superar este problema, vários pesquisadores propuseram equações para representar matematicamente o comportamento de solos não saturados. Estas proposições são baseadas em índices físicos, caracterização do solo, em ensaios convencionais ou simplesmente em ajustes de curvas. A relação entre a umidade e a sucção matricial, convencionalmente denominada curva característica de sucção do solo (SWCC) é também uma ferramenta útil na previsão do comportamento de engenharia de solos não saturados. Existem muitas equações para representar matematicamente a SWCC. Algumas são baseadas no pressuposto de que sua forma está diretamente relacionada com a distribuição dos poros e, portanto, com a granulometria. Nestas proposições, os parâmetros são calibrados pelo ajuste da curva de dados experimentais. Outros métodos supõem que a curva pode ser estimada diretamente a partir de propriedades físicas dos solos. Estas propostas são simples e conveniente para a utilização prática, mas são substancialmente incorretas, uma vez que ignoram a influência do teor de umidade, nível de tensões, estrutura do solo e mineralogia. Como resultado, a maioria tem sucesso limitado, dependendo do tipo de solo. Algumas tentativas têm sido feitas para prever a variação da resistência ao cisalhamento com relação a sucção matricial. Estes procedimentos usam, como uma ferramenta, direta ou indiretamente, a SWCC em conjunto com os parâmetros efetivos de resistência c e . Este trabalho discute a aplicabilidade de três equações para previsão da SWCC (Gardner, 1958; van Genuchten, 1980; Fredlund; Xing, 1994) para vinte e quatro amostras de solos residuais brasileiros. A adequação do uso da curva característica normalizada, proposta por Camapum de Carvalho e Leroueil (2004), também foi investigada. Os parâmetros dos modelos foram determinados por ajuste de curva, utilizando técnicas de problema inverso; dois métodos foram usados: algoritmo genético (AG) e Levenberq-Marquardt. Vários parâmetros que influênciam o comportamento da SWCC são discutidos. A relação entre a sucção matricial e resistência ao cisalhamento foi avaliada através de ajuste de curva utilizando as equações propostas por Öberg (1995); Sällfors (1997), Vanapalli et al., (1996), Vilar (2007); Futai (2002); oito resultados experimentais foram analisados. Os vários parâmetros que influênciam a forma da SWCC e a parcela não saturadas da resistência ao cisalhamento são discutidos.

Thermo-hydro-geomechanical simulations of methane gas production from deep sea methane hydrate formations

A new constitutive model called Methane Hydrate Critical State (MHCS) model was conducted to investigate the geomechanical response of the gas-hydrate-bearing sediments at the Nankai Trough during the wellbore construction process. The strength and dilatancy of gas-hydrate-bearing soil would gradually disappear when the bonds are destroyed because of excessively shearing, which are often observed in dense soils and also in bonded soils such as cemented soil and unsaturated soil. In this study, the MHCS model, which presents such softening features, would be incorporated into a staged-finite-element model in ABAQUS, which mainly considered the loading history of soils and the interaction between cement-casing-formation. This model shows the influence of gas-hydrate-bearing soil to the deformation and stability of a wellbore and the surrounding sediments during wellbore construction. At the same time, the conventional Mohr-Coulomb model was used in the model to show the advantages of MHCS model by comparing the results of the two models.

盐渍土的增湿剪切变形特性

盐渍土在浸水后由于盐分的溶解,对土体的强度产生较大的影响。在以往的研究中,主要以饱和与非饱和土体的强度比较来说明水对强度的影响,没有研究在部分浸水后土体的变形和强度的变化规律。针对这一点,在参考黄土增湿剪切方法的基础上研究了盐渍土的增湿剪切变形特性,探讨了垂直压力、干密度和剪应力水平对剪切变形量的影响规律。