Iron nodules in ferric soils of the Fraser Coast, Australia : relicts of laterisation or features of contemporary weathering and pedogenesis?

The genesis of ferruginous nodules and pisoliths in soils and weathering profiles of coastal southern and eastern Australia has long been debated. It is not clear whether iron (Fe) nodules are redox accumulations, residues of Miocene laterite duricrust, or the products of contemporary weathering of Fe-rich sedimentary rocks. This study combines a catchment-wide survey of Fe nodule distribution in Poona Creek catchment (Fraser Coast, Queensland) with detailed investigations of a representative ferric soil profile to show that Fe nodules are derived from Fe-rich sandstones. Where these crop out, they are broken down, transported downslope by colluvial processes, and redeposited. Chemical and physical weathering transforms these eroded rock fragments into non-magnetic Fe nodules. Major features of this transformation include lower hematite/goethite and kaolinite/gibbsite ratios, increased porosity, etching of quartz grains, and development of rounded morphology and a smooth outer cortex. Iron nodules are commonly concentrated in ferric horizons. We show that these horizons form as the result of differential biological mixing of the soil. Bioturbation gradually buries nodules and rock fragments deposited at the surface of the soil, resulting in a largely nodule-free 'biomantle' over a ferric 'stone line'. Maghemite-rich magnetic nodules are a prominent feature of the upper half of the profile. These are most likely formed by the thermal alteration of non-magnetic nodules located at the top of the profile during severe bushfires. They are subsequently redistributed through the soil profile by bioturbation. Iron nodules occurring in the study area are products of contemporary weathering of Fe-rich rock units. They are not laterite duricrust residues nor are they redox accumulations, although redox-controlled dissolution/re-precipitation is an important component of post-depositional modification of these Fe nodules.

Laboratory measurement of hydraulic conductivity functions of two unsaturated sandy soils during drying and wetting processes

The importance of applying unsaturated soil mechanics to geotechnical engineering design has been well understood. However, the consumption of time and the necessity for a specific laboratory testing apparatus when measuring unsaturated soil properties have limited the application of unsaturated soil mechanics theories in practice. Although methods for predicting unsaturated soil properties have been developed, the verification of these methods for a wide range of soil types is required in order to increase the confidence of practicing engineers in using these methods. In this study, a new permeameter was developed to measure the hydraulic conductivity of unsaturated soils using the steady-state method and directly measured suction (negative pore-water pressure) values. The apparatus is instrumented with two tensiometers for the direct measurement of suction during the tests. The apparatus can be used to obtain the hydraulic conductivity function of sandy soil over a low suction range (0-10 kPa). Firstly, the repeatability of the unsaturated hydraulic conductivity measurement, using the new permeameter, was verified by conducting tests on two identical sandy soil specimens and obtaining similar results. The hydraulic conductivity functions of the two sandy soils were then measured during the drying and wetting processes of the soils. A significant hysteresis was observed when the hydraulic conductivity was plotted against the suction. However, the hysteresis effects were not apparent when the conductivity was plotted against the volumetric water content. Furthermore, the measured unsaturated hydraulic conductivity functions were compared with predictions using three different predictive methods that are widely incorporated into numerical software. The results suggest that these predictive methods are capable of capturing the measured behavior with reasonable agreement.

Two-scale computational modelling of unsaturated flow in soils exhibiting small scale heterogeneities

Unsaturated water flow in soil is commonly modelled using Richards’ equation, which requires the hydraulic properties of the soil (e.g., porosity, hydraulic conductivity, etc.) to be characterised. Naturally occurring soils, however, are heterogeneous in nature, that is, they are composed of a number of interwoven homogeneous soils each with their own set of hydraulic properties. When the length scale of these soil heterogeneities is small, numerical solution of Richards’ equation is computationally impractical due to the immense effort and refinement required to mesh the actual heterogeneous geometry. A classic way forward is to use a macroscopic model, where the heterogeneous medium is replaced with a fictitious homogeneous medium, which attempts to give the average flow behaviour at the macroscopic scale (i.e., at a scale much larger than the scale of the heterogeneities). Using the homogenisation theory, a macroscopic equation can be derived that takes the form of Richards’ equation with effective parameters. A disadvantage of the macroscopic approach, however, is that it fails in cases when the assumption of local equilibrium does not hold. This limitation has seen the introduction of two-scale models that include at each point in the macroscopic domain an additional flow equation at the scale of the heterogeneities (microscopic scale). This report outlines a well-known two-scale model and contributes to the literature a number of important advances in its numerical implementation. These include the use of an unstructured control volume finite element method and image-based meshing techniques, that allow for irregular micro-scale geometries to be treated, and the use of an exponential time integration scheme that permits both scales to be resolved simultaneously in a completely coupled manner. Numerical comparisons against a classical macroscopic model confirm that only the two-scale model correctly captures the important features of the flow for a range of parameter values.

Resistance of microbial populations in DDT-contaminated and uncontaminated soils

One DDT-contaminated soil and two uncontaminated soils were used to enumerate DDT-resistant microbes (bacteria, actinomycetes and fungi) by using soil dilution agar plates in media either with 150 μg DDT ml -1 or without DDT at different temperatures (25, 37 and 55°C). Microbial populations in this study were significantly (p<0.001) affected by DDT in the growth medium. However, the numbers of microbes in long-term contaminated and uncontaminated soils were similar, presumably indicating that DDT-resistant microbes had developed over a long time exposure. The tolerance of isolated soil microbes to DDT varied in the order fungi>actinomycetes>bacteria. Bacteria from contaminated soil were more resistant to DDT than bacteria from uncontaminated soils. Microbes isolated at different temperatures also demonstrated varying degrees of DDT resistance. For example, bacteria and actinomycetes isolated at all incubation temperatures were sensitive to DDT. Conversely fungi isolated at all temperatures were unaffected by DDT.

DDT resistance and transformation by different microbial strains isolated from DDT-contaminated soils and compost materials

Bioremediation is a potential option to treat 1, 1, 1-trichloro-2, 2 bis (4-chlorophenyl) ethane (DDT) contaminated sites. In areas where suitable microbes are not present, the use of DDT resistant microbial inoculants may be necessary. It is vital that such inoculants do not produce recalcitrant breakdown products e.g. 1, 1-dichloro-2, 2-bis (4-chlorophenyl) ethylene (DDE). Therefore, this work aimed to screen DDT-contaminated soil and compost materials for the presence of DDT-resistant microbes for use as potential inoculants. Four compost amended soils, contaminated with different concentrations of DDT, were used to isolate DDT-resistant microbes in media containing 150 mg I -1 DDT at three temperatures (25, 37 and 55°C). In all soils, bacteria were more sensitive to DDT than actinomycetes and fungi. Bacteria isolated at 55°C from any source were the most DDT sensitive. However DDT-resistant bacterial strains showed more promise in degrading DDT than isolated fungal strains, as 1, 1-dichloro 2, 2-bis (4-chlorophenyl) ethane (DDD) was a major bacterial transformation product, while fungi tended to produce more DDE. Further studies on selected bacterial isolates found that the most promising bacterial strain (Bacillus sp. BHD-4) could remove 51% of DDT from liquid culture after 7 days growth. Of the amount transformed, 6% was found as DDD and 3% as DDE suggesting that further transformation of DDT and its metabolites occurred.

The relationship between N2O, NO, and N2 fluxes from fertilized and irrigated dryland soils of the Aral Sea Basin, Uzbekistan

Microbial respiratory reduction of nitrous oxide (N2O) to dinitrogen (N2) via denitrification plays a key role within the global N-cycle since it is the most important process for converting reactive nitrogen back into inert molecular N2. However, due to methodological constraints, we still lack a comprehensive, quantitative understanding of denitrification rates and controlling factors across various ecosystems. We investigated N2, N2O and NO emissions from irrigated cotton fields within the Aral Sera Basin using the He/O2 atmosphere gas flow soil core technique and an incubation assay. NH4NO3 fertilizer, equivalent to 75 kg ha−1 and irrigation water, adjusting the water holding capacity to 70, 100 and 130% were applied to the incubation vessels to assess its influence on gaseous N emissions. Under soil conditions as they are naturally found after concomitant irrigation and fertilization, denitrification was the dominant process and N2 the main end product of denitrification. The mean ratios of N2/N2O emissions increased with increasing soil moisture content. N2 emissions exceeded N2O emissions by a factor of 5 ± 2 at 70% soil water holding capacity (WHC) and a factor of 55 ± 27 at 130% WHC. The mean ratios of N2O/NO emissions varied between 1.5 ± 0.4 (70% WHC) and 644 ± 108 (130% WHC). The magnitude of N2 emissions for irrigated cotton was estimated to be in the range of 24 ± 9 to 175 ± 65 kg-N ha−1season−1, while emissions of NO were only of minor importance (between 0.1 to 0.7 kg-N ha−1 season−1). The findings demonstrate that for irrigated dryland soils in the Aral Sera Basin, denitrification is a major pathway of N-loss and that substantial amounts of N-fertilizer are lost as N2 to the atmosphere for irrigated dryland soils.

Greenhouse gas emissions from sub-tropical agricultural soils after addition of organic by-products

As the cost of mineral fertilisers increases globally, organic soil amendments (OAs) from agricultural sources are increasingly being used as substitutes for nitrogen. However, the impact of OAs on the production of greenhouse gases (CO2 and N2O) is not well understood. A 60-day laboratory incubation experiment was conducted to investigate the impacts of applying OAs (equivalent to 296 kg N ha−1 on average) on N2O and CO2 emissions and soil properties of clay and sandy loam soils from sugar cane production. The experiment included 6 treatments, one being an un-amended (UN) control with addition of five OAs being raw mill mud (MM), composted mill mud (CM), high N compost (HC), rice husk biochar (RB), and raw mill mud plus rice husk biochar (MB). These OAs were incubated at 60, 75 and 90% water-filled pore space (WFPS) at 25°C with urea (equivalent to 200 kg N ha−1) added to the soils thirty days after the incubation commenced. Results showed WFPS did not influence CO2 emissions over the 60 days but the magnitude of emissions as a proportion of C applied was RB < CM < MB < HC

Toward a model for predicting magnetic susceptibility of bedrock regolith and soils

One of the Department of Defense's most pressing environmental problems is the efficient detection and identification of unexploded ordnance (UXO). In regions of highly magnetic soils, magnetic and electromagnetic sensors often detect anomalies that are of geologic origin, adding significantly to remediation costs. In order to develop predictive models for magnetic susceptibility, it is crucial to understand modes of formation and the spatial distribution of different iron oxides. Most rock types contain iron and their magnetic susceptibility is determined by the amount and form of iron oxides present. When rocks weather, the amount and form of the oxides change, producing concomitant changes in magnetic susceptibility. The type of iron oxide found in the weathered rock or regolith is a function of the duration and intensity of weathering, as well as the original content of iron in the parent material. The rate of weathering is controlled by rainfall and temperature; thus knowing the climate zone, the amount of iron in the lithology and the age of the surface will help predict the amount and forms of iron oxide. We have compiled analyses of the types, amounts, and magnetic properties of iron oxides from soils over a wide climate range, from semi arid grasslands, to temperate regions, and tropical forests. We find there is a predictable range of iron oxide type and magnetic susceptibility according to the climate zone, the age of the soil and the amount of iron in the unweathered regolith.

Variation in parental magmas of Mt Rouse, a complex polymagmatic monogenetic volcano in the basaltic intraplate Newer Volcanics Province, southeast Australia

Monogenetic volcanoes have long been regarded as simple in nature, involving single magma batches and uncomplicated evolutions; however, recent detailed research into individual centres is challenging that assumption. Mt Rouse (Kolor) is the volumetrically largest volcano in the monogenetic Newer Volcanics Province of southeast Australia. This study presents new major, trace and Sr–Nd–Pb isotope data for samples selected on the basis of a detailed stratigraphic framework analysis of the volcanic products from Mt Rouse. The volcano is the product of three magma batches geochemically similar to Ocean–Island basalts, featuring increasing LREE enrichment with each magma batch (batches A, B and C) but no evidence of crustal contamination; the Sr–Nd–Pb isotopes define two groupings. Modelling suggests that the magmas were sourced from a zone of partial melting crossing the lithosphere–asthenosphere boundary, with batch A forming a large volume partial melt in the deep lithosphere (1.7 GPa/55.5 km); and batches B and C from similar areas within the shallow asthenosphere (1.88 GPa/61 km and 1.94 GPa/63 km, respectively). The formation and extraction of these magmas may have been due to high deformation rates in the mantle caused by edge-driven convection and asthenospheric upwelling. The lithosphere– asthenosphere boundary is important with respect to NVP volcanism. An eruption chronology involves sequential eruption of magma batches A, C and B, followed by simultaneous eruption of batches A and B. Mt Rouse is a complex polymagmatic monogenetic volcano that illustrates the complexity of monogenetic volcanism and demonstrates the importance of combining detailed stratigraphic analysis alongside systematic geochemical sampling.

Are we getting accurate measurements of Ksat for sodic clay soils?

In this paper we discuss the use of a series of column experiments to improve understanding of the effect irrigation water chemistry (saline solutions) has on measurements of saturated hydraulic conductivity (Ksat) of a sodic clay soil. We highlight in particular the use of extended leaching periods to determine whether the duration of leaching affects the results. In the experiments, mixed cation solutions of two different salinity levels, 50 meq/L and 100 meq/L, were applied under constant head to columns of a repacked sodic clay soil using three replicates for each treatment. The maximum Ksat measured during leaching with the 100 meq/L solution was approximately double the maximum Ksat measured during leaching with the 50 meq/L solution. Measured flow rates were found to increase rapidly after flow commenced then decrease gradually until flow rates became stable. The final, stable flow rate was roughly 80% less than the maximum flow rate measured. Reasons for these changes in saturated hydraulic conductivity are discussed. The key finding from these experiments is that long term leaching, involving significantly more pore volumes than is commonly reported in the literature, is required to obtain a ‘stable’ Ksat. We recommend that further studies be carried out to (1) determine whether similar behaviour in Ksat occurs in a wide range of sodic clay soils and (2) to help build a better understanding of the causes and implications of the observed behaviour in Ksat.

Two-scale computational modelling of water flow in unsaturated soils containing irregular-shaped inclusions

The focus of this paper is two-dimensional computational modelling of water flow in unsaturated soils consisting of weakly conductive disconnected inclusions embedded in a highly conductive connected matrix. When the inclusions are small, a two-scale Richards’ equation-based model has been proposed in the literature taking the form of an equation with effective parameters governing the macroscopic flow coupled with a microscopic equation, defined at each point in the macroscopic domain, governing the flow in the inclusions. This paper is devoted to a number of advances in the numerical implementation of this model. Namely, by treating the micro-scale as a two-dimensional problem, our solution approach based on a control volume finite element method can be applied to irregular inclusion geometries, and, if necessary, modified to account for additional phenomena (e.g. imposing the macroscopic gradient on the micro-scale via a linear approximation of the macroscopic variable along the microscopic boundary). This is achieved with the help of an exponential integrator for advancing the solution in time. This time integration method completely avoids generation of the Jacobian matrix of the system and hence eases the computation when solving the two-scale model in a completely coupled manner. Numerical simulations are presented for a two-dimensional infiltration problem.

Damage assessment of basaltic rock mass due to repeated blasting in a railway tunnelling project - A case study

The demand for tunnelling and underground space creation is rapidly growing due to the requirement of civil infrastructure projects and urbanisation. Blasting remains the most inexpensive method of underground excavations in hard rock. Unfortunately, there are no specific safety guidelines available for the blasted tunnels with regards to the threshold limits of vibrations caused by repeated blasting activity in the close proximity. This paper presents the results of a comprehensive study conducted to find out the effect of repeated blast loading on the damage experienced by jointed basaltic rock mass during tunnelling works. Conducting of multiple rounds of blasts for various civil excavations in a railway tunnel imparted repeated loading on rock mass of sidewall and roof of the tunnel. The blast induced damage was assessed by using vibration attenuation equations of charge weight scaling law and measured by borehole extensometers and borehole camera. Ground vibrations of each blasting round were also monitored by triaxial geophones installed near the borehole extensometers. The peak particle velocity (V-max) observations and plastic deformations from borehole extensometers were used to develop a site specific damage model. The study reveals that repeated dynamic loading imparted on the exposed tunnel from subsequent blasts, in the vicinity, resulted in rock mass damage at lesser vibration levels than the critical peak particle velocity (V-cr). It was found that, the repeated blast loading resulted in the near-field damage due to high frequency waves and far-field damage due to low frequency waves. The far field damage, after 45-50 occurrences of blast loading, was up to 55% of the near-field damage in basaltic rock mass. The findings of the study clearly indicate that the phenomena of repeated blasting with respect to number of cycles of loading should be taken into consideration for proper assessment of blast induced damage in underground excavations.

Consolidation Behavior of Soils

Based on Terzaghi's consolidation theory, percent of consolidation, U, versus the time factor, T, relationship for constant/linear excess pore water pressure distribution, it is possible to generate theoretical log10(H2/t) versus U curves where H is the length of the drainage path of a consolidating layer, and t is the time for different known values of the coefficient of consolidation, cν. A method has been developed wherein both the theoretical and experimental behavior of soils during consolidation can be simultaneously compared and studied on the same plot. The experimental log10(H2/t) versus U curves have been compared with the theoretical curves. The deviations of the experimental behavior from the theory are explained in terms of initial compression and secondary compression. Analysis of results indicates that the secondary compression essentially starts from about 60% consolidation. A simple procedure is presented for calculating the value of cv from the δ-t data using log10(H2/t) versus U plot.