Applying HYDRUS to flow in a sodic clay soil with solution composition–dependent hydraulic conductivity

Management of sodic soils under irrigation often requires application of chemical ameliorants to improve permeability combined with leaching of excess salts. Modeling irrigation, soil treatments, and leaching in these sodic soils requires a model that can adequately represent the physical and chemical changes in the soil associated with the amelioration process. While there are a number of models that simulate reactive solute transport, UNSATCHEM and HYDRUS-1D are currently the only models that also include an ability to simulate the impacts of soil chemistry on hydraulic conductivity. Previous researchers have successfully applied these models to simulate amelioration experiments on a sodic loam soil. To further gauge their applicability, we extended the previous work by comparing HYDRUS simulations of sodic soil amelioration with the results from recently published laboratory experiments on a more reactive, repacked sodic clay soil. The general trends observed in the laboratory experiments were able to be simulated using HYDRUS. Differences between measured and simulated results were attributed to the limited flexibility of the function that represents chemistry-dependent hydraulic conductivity in HYDRUS. While improvements in the function could be made, the present work indicates that HYDRUS-UNSATCHEM captures the key changes in soil hydraulic properties that occur during sodic clay soil amelioration and thus extends the findings of previous researchers studying sodic loams.

Induced swelling of kaolinitic soil in alkali solution

Unexpected swelling induced in foundation soils can cause distress to structures founded on them. In this paper, the swelling of kaolinitic soils due to interaction with alkali solution has been reported. The induced swelling is attributed to the formation of new minerals, which has been confirmed by X-ray diffraction patters and SEM studies. To understand the effect of alkali concentration and duration of interaction, two series of consolidation experiments have been carried out. In series 1, the specimen were remoulded with water and inundated with alkali solutions and in series 2, the specimen were remoulded and inundated with same alkali solutions. A steep compression during loading cycle and no abnormal swelling during unloading cycle has been noticed for the specimen remoulded with water and inundated with 1 N NaOH solutions. The steep compression is due to the segregation or break down of clay minerals due to alkali interactions. In case of specimen inundated with 4 N NaOH solutions, abnormal swelling has been observed during unloading cycle of the consolidation test. New minerals are formed on interaction of soil with 4 N solution as confirmed by X-ray diffraction patterns. These minerals are known to have very fine pores and possess high water holding capacity. The differences in the amount of swelling of samples remoulded with water and remoulded with alkali solution are due to variations in the concentration of alkali and duration of interaction.

Mineralogical changes and geotechnical properties of an expansive soil interacted with caustic solution

The type and amount of clay mineral plays an important role in the behaviour of fine-grained soils. Clay minerals are the primary source and moisture is often the external agent of swelling in soils. Also soils may exhibit increased/reduced swelling due to interaction with chemicals. Alkalis used in industrial operations are one such example. Concentrations of alkali and mineral type are the key factors in such interactions. The present paper reports the changes in the properties of an expansive Black Cotton soil containing a mixed layer mineral, rectorite upon interaction with high concentration caustic solutions. X-ray diffraction studies have shown that the rectorite present in the soil undergoes changes with increase in the concentration of alkali. Saponite gets transformed to nantronite. Small amount of kaolinitic mineral present in the soil also reacts with alkali producing some changes in its mineralogy. Many hydroxides are produced. Differential thermal analysis studies have been supportive of these changes. Consequent of these changes, the soil-specific surface increases, changes its Atterberg limits and free swell volume increases. The results have been supported by the characteristics and behaviour of samples contaminated in the field with alkali from an alumina extraction plant.

An Analytical Solution for Three-Dimensional Diffraction of Plane P-Waves by a Hemispherical Alluvial Valley with Saturated Soil Deposits

An analytical solution for the three-dimensional scattering and diffraction of plane P-waves by a hemispherical alluvial valley with saturated soil deposits is developed by employing Fourier-Bessel series expansion technique. Unlike previous studies, in which the saturated soil deposits were simulated with the single-phase elastic theory, in this paper, they are simulated with Biot's dynamic theory for saturated porous media, and the half space is assumed as a single-phase elastic medium. The effects of the dimensionless frequency, the incidence angle of P-wave and the porosity of soil deposits on the surface displacement magnifications of the hemispherical alluvial valley are investigated. Numerical results show that the existence of a saturated hemispherical alluvial valley has much influence on the surface displacement magnifications. It is more reasonable to simulate soil deposits with Biot's dynamic theory when evaluating the displacement responses of a hemispherical alluvial valley with an incidence of P-waves.

A finite element solution scheme for an elastic–viscoplastic soil model

A Newton–Raphson solution scheme with a stress point algorithm is presented for the implementation of an elastic–viscoplastic soilmodel in a finite element program. Viscoplastic strain rates are calculated using the stress and volumetric states of the soil. Sub-incrementsof time are defined for each iterative calculation of elastic–viscoplastic stress changes so that their sum adds up to the time incrementfor the load step. This carefully defined ‘iterative time’ ensures that the correct amount of viscoplastic straining is accumulated overthe applied load step. The algorithms and assumptions required to implement the solution scheme are provided. Verification of the solutionscheme is achieved by using it to analyze typical boundary value problems.

Sequential exhaustive extraction of a Mollisol soil, and characterizations of humic components, including humin, by solid and solution state NMR.

A comprehensive sequential extraction procedure was applied to isolate soil organic components using aqueous solvents at different pH values, base plus urea (base-urea), and finally dimethylsulfoxide (DMSO) plus concentrated H2SO4 (DMSO-acid) for the humin-enriched clay separates. The extracts from base-urea and DMSO-acid would be regarded as 'humin' in the classical definitions. The fractions isolated from aqueous base, base-urea and DMSO-acid were characterized by solid and solution state NMR spectroscopy. The base-urea solvent system isolated ca. 10% (by mass) additional humic substances. The combined base-urea and DMSO-acid solvents isolated ca. 93% of total organic carbon from the humin-enriched fine clay fraction (<2 ?m). Characterization of the humic fractions by solid-state NMR spectroscopy showed that oxidized char materials were concentrated in humic acids isolated at pH 7, and in the base-urea extract. Lignin-derived materials were in considerable abundance in the humic acids isolated at pH 12.6. Only very small amounts of char-derived structures were contained in the fulvic acids and fulvic acids-like material isolated from the base-urea solvent. After extraction with base-urea, the 0.5 m NaOH extract from the humin-enriched clay was predominantly composed of aliphatic hydrocarbon groups, and with lesser amounts of aromatic carbon (probably including some char material), and carbohydrates and peptides. From the combination of solid and solution-state NMR spectroscopy, it is clear that the major components of humin materials, from the DMSO-acid solvent, after the exhaustive extraction sequence, were composed of microbial and plant derived components, mainly long-chain aliphatic species (including fatty acids/ester, waxes, lipids and cuticular material), carbohydrate, peptides/proteins, lignin derivatives, lipoprotein and peptidoglycan (major structural components in bacteria cell walls). Black carbon or char materials were enriched in humic acids isolated at pH 7 and humic acids-like material isolated in the base-urea medium, indicating that urea can liberate char-derived material hydrogen bonded or trapped within the humin matrix.

Oxygen transport in soil and the vertical distribution of roots

An analytical solution for steady-state oxygen transport in soils including 2 sink terms, viz roots and microbes with the corresponding vertical distribution scaling lengths forming a ratio p, showed p governed the critical air-filled porosity, θ_c, needed by most plants. For low temperature and p, θ_c was <0.1 but at higher temperatures and p = 1, θ_c was >0.15 m³/m³. When root length density at the surface was 10⁴ m/m³ and p > 3, θ_c was 0.25 m³/m³, more than half the pore space. Few combinations of soil and climate regularly meet this condition. However, for sandy soils and seasonally warm, arid regions, the theory is consistent with observation, in that plants may have some deep roots. Critical θ_c values are used to formulate theoretical solutions in a forward mode, so different levels of oxygen uptake by roots may be compared to microbial activity. The proportion of respiration by plant roots increases rapidly with p up to p ≈2.

Controls on iron in soils and soil waters of a forested, coastal catchment in subtropical Australia

Soluble organic matter derived from exotic Pinus vegetation forms stronger complexes with iron (Fe) than the soluble organic matter derived from most native Australian species. This has lead to concern about the environmental impacts related to the establishment of extensive exotic Pinus plantations in coastal southeast Queensland, Australia. It has been suggested that the Pinus plantations may enhance the solubility of Fe in soils by increasing the amount of organically complexed Fe. While this remains inconclusive, the environmental impacts of an increased flux of dissolved, organically complexed Fe from soils to the fluvial system and then to sensitive coastal ecosystems are potentially damaging. Previous work investigated a small number of samples, was largely laboratory based and had limited application to field conditions. These assessments lacked field-based studies, including the comparison of the soil water chemistry of sites associated with Pinus vegetation and undisturbed native vegetation. In addition, the main controls on the distribution and mobilisation of Fe in soils of this subtropical coastal region have not been determined. This information is required in order to better understand the relative significance of any Pinus enhanced solubility of Fe. The main aim of this thesis is to determine the controls on Fe distribution and mobilisation in soils and soil waters of a representative coastal catchment in southeast Queensland (Poona Creek catchment, Fraser Coast) and to test the effect of Pinus vegetation on the solubility and speciation of Fe. The thesis is structured around three individual papers. The first paper identifies the main processes responsible for the distribution and mobilisation of labile Fe in the study area and takes a catchment scale approach. Physicochemical attributes of 120 soil samples distributed throughout the catchment are analysed, and a new multivariate data analysis approach (Kohonen’s self organising maps) is used to identify the conditions associated with high labile Fe. The second paper establishes whether Fe nodules play a major role as an iron source in the catchment, by determining the genetic mechanism responsible for their formation. The nodules are a major pool of Fe in much of the region and previous studies have implied that they may be involved in redox-controlled mobilisation and redistribution of Fe. This is achieved by combining a detailed study of a ferric soil profile (morphology, mineralogy and micromorphology) with the distribution of Fe nodules on a catchment scale. The third component of the thesis tests whether the concentration and speciation of Fe in soil solutions from Pinus plantations differs significantly from native vegetation soil solutions. Microlysimeters are employed to collect unaltered, in situ soil water samples. The redox speciation of Fe is determined spectrophotometrically and the interaction between Fe and dissolved organic matter (DOM) is modelled with the Stockholm Humic Model. The thesis provides a better understanding of the controls on the distribution, concentration and speciation of Fe in the soils and soil waters of southeast Queensland. Reductive dissolution is the main mechanism by which mobilisation of Fe occurs in the study area. Labile Fe concentrations are low overall, particularly in the sandy soils of the coastal plain. However, high labile Fe is common in seasonally waterlogged and clay-rich soils which are exposed to fluctuating redox conditions and in organic-rich soils adjacent to streams. Clay-rich soils are most common in the upper parts of the catchment. Fe nodules were shown to have a negligible role in the redistribution of dissolved iron in the catchment. They are formed by the erosion, colluvial transport and chemical weathering of iron-rich sandstones. The ferric horizons, in which nodules are commonly concentrated, subsequently form through differential biological mixing of the soil. Whereas dissolution/ reprecipitation of the Fe cements is an important component of nodule formation, mobilised Fe reprecipitates locally. Dissolved Fe in the soil waters is almost entirely in the ferrous form. Vegetation type does not affect the concentration and speciation of Fe in soil waters, although Pinus DOM has greater acidic functional group site densities than DOM from native vegetation. Iron concentrations are highest in the high DOM soil waters collected from sandy podosols, where they are controlled by redox potential. Iron concentrations are low in soil solutions from clay and iron oxide rich soils, in spite of similar redox potentials. This is related to stronger sorption to the reactive clay and iron oxide mineral surfaces in these soils, which reduces the amount of DOM available for microbial metabolisation and reductive dissolution of Fe. Modelling suggests that Pinus DOM can significantly increase the amount of truly dissolved ferric iron remaining in solution in oxidising conditions. Thus, inputs of ferrous iron together with Pinus DOM to surface waters may reduce precipitation of hydrous ferric oxides and increase the flux of dissolved iron out of the catchment. Such inputs are most likely from the lower catchment, where podosols planted with Pinus are most widely distributed. Significant outcomes other than the main aims were also achieved. It is shown that mobilisation of Fe in podosols can occur as dissolved Fe(II) rather than as Fe(III)-organic complexes. This has implications for the large body of work which assumes that Fe(II) plays a minor role. Also, the first paper demonstrates that a data analysis approach based on Kohonen’s self organising maps can facilitate the interpretation of complex datasets and can help identify geochemical processes operating on a catchment scale.

The influence of vegetation and soil type on the speciation of iron in soil water

Soluble organic matter derived from exotic Pinus species has been shown to form stronger complexes with iron (Fe) than that derived from most native Australian species. It has also been proposed that the establishment of exotic Pinus plantations in coastal southeast Queensland may have enhanced the solubility of Fe in soils by increasing the amount of organically complexed Fe, but this remains inconclusive. In this study we test whether the concentration and speciation of Fe in soil water from Pinus plantations differs significantly from soil water from native vegetation areas. Both Fe redox speciation and the interaction between Fe and dissolved organic matter (DOM) were considered; Fe - DOM interaction was assessed using the Stockholm Humic Model. Iron concentrations (mainly Fe 2+) were greatest in the soil waters with the greatest DOM content collected from sandy podosols (Podzols), where they are largely controlled by redox potential. Iron concentrations were small in soil waters from clay and iron oxide-rich soils, in spite of similar redox potentials. This condition is related to stronger sorption on to the reactive clay and iron oxide mineral surfaces in these soils, which reduces the amount of DOM available for electron shuttling and microbial metabolism, restricting reductive dissolution of Fe. Vegetation type had no significant influence on the concentration and speciation of iron in soil waters, although DOM from Pinus sites had greater acidic functional group site densities than DOM from native vegetation sites. This is because Fe is mainly in the ferrous form, even in samples from the relatively well-drained podosols. However, modelling suggests that Pinus DOM can significantly increase the amount of truly dissolved ferric iron remaining in solution in oxic conditions. Therefore, the input of ferrous iron together with Pinus DOM to surface waters may reduce precipitation of hydrous ferric oxides (ferrihydrite) and increase the flux of dissolved Fe out of the catchment. Such inputs of iron are most probably derived from podosols planted with Pinus.

Soil moisture monitoring at the field scale using neutron probe

Measurement of the moisture variation in soils is required for geotechnical design and research because soil properties and behavior can vary as moisture content changes. The neutron probe, which was developed more than 40 years ago, is commonly used to monitor soil moisture variation in the field. This study reports a full-scale field monitoring of soil moisture using a neutron moisture probe for a period of more than 2 years in the Melbourne (Australia) region. On the basis of soil types available in the Melbourne region, 23 sites were chosen for moisture monitoring down to a depth of 1500 mm. The field calibration method was used to develop correlations relating the volumetric moisture content and neutron counts. Observed results showed that the deepest “wetting front” during the wet season was limited to the top 800 to 1000 mm of soil whilst the top soil layer down to about 550mmresponded almost immediately to the rainfall events. At greater depths (550 to 800mmand below 800 mm), the moisture variations were relatively low and displayed predominantly periodic fluctuations. This periodic nature was captured with Fourier analysis to develop a cyclic moisture model on the basis of an analytical solution of a one-dimensional moisture flow equation for homogeneous soils. It is argued that the model developed can be used to predict the soil moisture variations as applicable to buried structures such as pipes.