Assessment of physical and hydrological properties of biological soil crusts using x-ray microtomography and modelling

Biological soil crusts (BSCs) are formed by aggregates of soil particles and communities of microbial organisms and are common in all drylands. The role of BSCs on infiltration remains uncertain due to the lack of data on their role in affecting soil physical properties such as porosity and structure. Quantitative assessment of these properties is primarily hindered by the fragile nature of the crusts. Here we show how the use of a combination of non-destructive imaging X-ray microtomography (XMT) and Lattice Boltzmann method (LBM) enables quantification of key soil physical parameters and the modeling of water flow through BSCs samples from Kalahari Sands, Botswana. We quantify porosity and flow changes as a result of mechanical disturbance of such a fragile cyanobacteria-dominated crust. Results show significant variations in porosity between different types of crusts and how they affect the flow and that disturbance of a cyanobacteria-dominated crust results in the breakdown of larger pore spaces and reduces flow rates through the surface layer. We conclude that the XMT–LBM approach is well suited for study of fragile surface crust samples where physical and hydraulic properties cannot be easily quantified using conventional methods.

Diagenesis and aquifer characteristics of the Sherwood Sandstone Group in North Nottinghamshire

The Sherwood Sandstone Group forms an important aquifer in Eastern England, which in North Nottinghamshire comprises the Nottingham Castle and Lenton Sandstone Formations. The aquifer is formed by an alluvial red-bed sequence dominated by medium-coarse grained sandstones which are texturally immature to submature and have only been subjected to shallow burial diagenesis. These sandstones reached the mature stage of the meso diagenetic regime, and four stages are recognized in their diagenetic history depending upon the physical/chemical processes prevailing and the subsequent effect on porosity and permeability. Stage "One" represents changes including dissolution of unstable silicates, clay replacement, red colouration and precipitation of authigenic minerals (quartz, feldspar, illite, l/S, kaolinite, dolomite, ferroan calcite, calcite). The net result of these changes was porosity reduction. Stage "Two" included changes due to mechanical compaction which resulted in minor porosity reduction. Stage "Three" was the main phase of secondary porosity enhancement. Stage "Four" represents changes taking place in the present groundwater where porosity and permeability may have been increased by dissolution and partly reduced by kaolinite precipitation. Porosity measured by water-resaturation and Hg-injection gave average values of 25.63% and 24.85% respectively. The results are comparable and showed marked correlation especially in highly porous/permeable rocks. Porosity measurements from photomicrographs were markedly offset from laboratory results. Horizontal Kw ranged between 1.43 x 10-5 and 1.13 x 10-1 mm/sec, with an average of  1.68 x 10-2 mm/sec. The estimated KHg ranged between 7.29 x 10-6 and 6.99 x 10-2 mm/sec with an average of 1.47 x 10-2 mm/sec. Both results are significantly correlated for highly porous/permeable rocks. The hydraulic properties are highly dependent upon the diagenetic properties (as most of the pores present are of secondary origin) as well as the pore size distribution. The chemistry of these groundwaters indicates that they are under-saturated with respect to dolomite, calcite, K-feldspar, l/S clay, and montmorillonite. The precipitation of kaolinite,and to a lesser extent illite, is favoured in the present groundwater regime.

Groundwater flow patterns within the Triassic sandstone aquifers, Birmingham, United Kingdom

This paper presents the development of a modelling study for part of the Birmingham area. Restricted access and model resolutions have limited wide applications of some of the previously developed models. The study area covers approximately 221 km2, and is underlain geologically, by a multi-layer setup with varied hydraulic properties. The basal aquifer unit is the Kidderminster sandstone Formation, overlain by the Wildmoor and Bromsgrove sandstone Formations. The presence of the Birmingham fault which acts as low permeability barrier demarcates the eastern and southern boundaries. The western boundary is defined by the presence of crystallised rocks and coal measures, while a groundwater divide defines the northern boundary. The estimated recharge flux is 112 mm/yr. The ranges of calibrated values obtained for horizontal and vertical hydraulic conductivities are 5.787x10-6 - 2.315x10-5  m/s and 5.787x10-8  - 1.157x10-7  m/s, respectively. Corresponding values obtained for the specific yield and specific storage are 0.10 - 0.12, and 1x10 -4 - 5x10 -4. The calculated numerical error is generally much less than 0.1 %. Hydraulic layering within the Permo-Triassic sandstone aquifer is thought to account for the large vertical anisotropy. Although, uncertainties are associated with the use of a simplistic delay approach to characterise the effects of the unsaturated zone, the modelled values are comparable with those obtained in the literature, and the flow pattern predictions appear to be realistic. © Research India Publications.

Mechanical properties of sheet metal under forming conditions

In the bulge test, a sheet metal specimen is clamped over a circular hole in a die and formed into a bulge by the hydraulic pressure on one side of the specirnen. As the unsupported part of the specimen is deformed in this way, its area is increased, in other words, the material is generally stretched and its thickness generally decreased. The stresses causing this stretching action are the membrane stresses in the shell generated by the hydraulic pressure, in the same way as the rubber in a toy balloon is stretched by the membrane stresses caused by the air inside it. The bulge test is a widely used sheet metal test, to determine the "formability" of sheet materials. Research on this forming process (2)-(15)* has hitherto been almost exclusively confined to predicting the behaviour of the bulged specimen through the constitutive equations (stresses and strains in relation to displacements and shapes) and empirical work hardening characteristics of the material as determined in the tension test. In the present study the approach is reversed; the stresses and strains in the specimen are measured and determined from the geometry of the deformed shell. Thus, the bulge test can be used for determining the stress-strain relationship in the material under actual conditions in sheet metal forming processes. When sheet materials are formed by fluid pressure, the work-piece assumes an approximately spherical shape, The exact nature and magnitude of the deviation from the perfect sphere can be defined and measured by an index called prolateness. The distribution of prolateness throughout the workpiece at any particular stage of the forming process is of fundamental significance, because it determines the variation of the stress ratio on which the mode of deformation depends. It is found. that, before the process becomes unstable in sheet metal, the workpiece is exactly spherical only at the pole and at an annular ring. Between the pole and this annular ring the workpiece is more pointed than a sphere, and outside this ring, it is flatter than a sphere. In the forming of sheet materials, the stresses and hence the incremental strains, are closely related to the curvatures of the workpiece. This relationship between geometry and state of stress can be formulated quantitatively through prolateness. The determination of the magnitudes of prolateness, however, requires special techniques. The success of the experimental work is due to the technique of measuring the profile inclination of the meridional section very accurately. A travelling microscope, workshop protractor and surface plate are used for measurements of circumferential and meridional tangential strains. The curvatures can be calculated from geometry. If, however, the shape of the workpiece is expressed in terms of the current radial (r) and axial ( L) coordinates, it is very difficult to calculate the curvatures within an adequate degree of accuracy, owing to the double differentiation involved. In this project, a first differentiation is, in effect, by-passed by measuring the profile inclination directly and the second differentiation is performed in a round-about way, as explained in later chapters. The variations of the stresses in the workpiece thus observed have not, to the knowledge of the author, been reported experimentally. The static strength of shells to withstand fluid pressure and their buckling strength under concentrated loads, both depend on the distribution of the thickness. Thickness distribution can be controlled to a limited extent by changing the work hardening characteristics of the work material and by imposing constraints. A technique is provided in this thesis for determining accurately the stress distribution, on which the strains associated with thinning depend. Whether a problem of controlled thickness distribution is tackled by theory, or by experiments, or by both combined, the analysis in this thesis supplies the theoretical framework and some useful experimental techniques for the research applied to particular problems. The improvement of formability by allowing draw-in can also be analysed with the same theoretical and experimental techniques. Results on stress-strain relationships are usually represented by single stress-strain curves plotted either between one stress and one strain (as in the tension or compression tests) or between the effective stress and effective strain, as in tests on tubular specimens under combined tension, torsion and internal pressure. In this study, the triaxial stresses and strains are plotted simultaneously in triangular coordinates. Thus, both stress and strain are represented by vectors and the relationship between them by the relationship between two vector functions. From the results so obtained, conclusions are drawn on both the behaviour and the properties of the material in the bulge test. The stress ratios are generally equal to the strain-rate ratios (stress vectors collinear with incremental strain vectors) and the work-hardening characteristics, which apply only to the particular strain paths are deduced. Plastic instability of the material is generally considered to have been reached when the oil pressure has attained its maximum value so that further deformation occurs under a constant or lower pressure. It is found that the instability regime of deformation has already occurred long before the maximum pressure is attained. Thus, a new concept of instability is proposed, and for this criterion, instability can occur for any type of pressure growth curves.

Development of a conceptual model of the soil-moisture-plant sub-system of the hydrological cycle

The soil-plant-moisture subsystem is an important component of the hydrological cycle. Over the last 20 or so years a number of computer models of varying complexity have represented this subsystem with differing degrees of success. The aim of this present work has been to improve and extend an existing model. The new model is less site specific thus allowing for the simulation of a wide range of soil types and profiles. Several processes, not included in the original model, are simulated by the inclusion of new algorithms, including: macropore flow; hysteresis and plant growth. Changes have also been made to the infiltration, water uptake and water flow algorithms. Using field data from various sources, regression equations have been derived which relate parameters in the suction-conductivity-moisture content relationships to easily measured soil properties such as particle-size distribution data. Independent tests have been performed on laboratory data produced by Hedges (1989). The parameters found by regression for the suction relationships were then used in equations describing the infiltration and macropore processes. An extensive literature review produced a new model for calculating plant growth from actual transpiration, which was itself partly determined by the root densities and leaf area indices derived by the plant growth model. The new infiltration model uses intensity/duration curves to disaggregate daily rainfall inputs into hourly amounts. The final model has been calibrated and tested against field data, and its performance compared to that of the original model. Simulations have also been carried out to investigate the effects of various parameters on infiltration, macropore flow, actual transpiration and plant growth. Qualitatively comparisons have been made between these results and data given in the literature.

Application of thermal neutron radiography for the mass transport of moisture through freezing soil

This thesis reports on the development of a technique to evaluate hydraulic conductivities in a soil (Snowcal) subject to freezing conditions. The technique draws on three distinctly different disciplines, Nuclear Physics, Soil Physics and Remote Sensing to provide a non-destructive and reliable evaluation of hydraulic conductivity throughout a freezing test. Thermal neutron radiography is used to provide information on local water/ice contents at anytime throughout the test. The experimental test rig is designed so that the soil matrix can be radiated by a neutron beam, from a nuclear reactor, to obtain radiographs. The radiographs can then be interpreted, following a process of remote sensing image enhancement, to yield information on relative water/ice contents. Interpretation of the radiographs is accommodated using image analysis equipment capable of distinguishing between 256 shades of grey. Remote sensing image enhancing techniques are then employed to develop false colour images which show the movement of water and development of ice lenses in the soil. Instrumentation is incorporated in the soil in the form of psychrometer/thermocouples, to record water potential, electrical resistance probes to enable ice and water to be differentiated on the radiographs and thermocouples to record the temperature gradient. Water content determinations are made from the enhanced images and plotted against potential measurements to provide the moisture characteristic for the soil. With relevant mathematical theory pore water distributions are obtained and combined with water content data to give hydraulic conductivities. The values for hydraulic conductivity in the saturated soil and at the frozen fringe are compared with established values for silts and silty-sands. The values are in general agreement and, with refinement, this non-destructive technique could afford useful information on a whole range of soils. The technique is of value over other methods because ice lenses are actually seen forming in the soil, supporting the accepted theories of frost action. There are economic and experimental restraints to the work which are associated with the use of a nuclear facility, however, the technique is versatile and has been applied to the study of moisture transfer in porous building materials and could be further developed into other research areas.

Characterization of engineered biochar for soil management

Pyrolysis is an energy conversion technology which by heating organic materials in the absence of oxygen, produces liquid, gaseous, and solid fuel products. Biochar, the solid product, can also be used as a soil amendment and, simultaneously, enables us to sequester carbon in the soil. By controlling the pyrolysis process, it is possible to engineer biochar suitable for the remediation of specific soil management problems. This research uses a characterization method more suited to producing biochar for soil amendment purposes than the existing biochar fuel characterization standards. This is the first research to use wastewater irrigated willow as a pyrolysis feedstock. The extensive characterization of biochar produced over a range of temperatures (410-810°C) yielded data on key properties relevant to soil under management: low surface area (1.4 to 5.4 m2/g), low bulk density (0.15-0.18 g/cm3), high pH values (7.8-9.4) and high water-holding capacity (1.8 to 4.3 cm3/g). Extraction experiments demonstrated low bioavailability of char nutrients (N, P, K, Ca, and Mg). This research also studied this artificial nitrogen cycle of pyrolysis: nitrogen accumulated in the wood from the wastewater and high levels of nitrogen remained in the biochar in a stable form not directly available to plants. Copyright © 2013 American Institute of Chemical Engineers Environ Prog.