Natural improvements of geochemical conditions of acid sulfate soils caused by free tidal inundation and its effects on the mangrove seedlings

Acid sulfate soils (ASS) are one of the stressor factors that cause many mangrove restoration projects to fail. Achieving successful rehabilitation in an ASS affected area requires an understanding of the geochemical conditions that influence the establishment and growth of mangrove seedlings. This study evaluated the effect of tidal inundation on geochemical conditions on sub layer to better understand their impacts on the density, establishment, and growth of mangrove seedlings. This study also examined the geochemical conditions under which mangrove seedlings can establish naturally, and/or be replanted in abandoned aquaculture ponds. The study area was in an area of abandoned aquaculture ponds situated in the Mare District, adjacent to Bone Bay, South Sulawesi, Indonesia.The pH, pHfox, redox potential, organic content, water soluble sulfate, SKCl, SPOS, and grain size of the soil from the sediment core at + 10 - 15 cm depth near roots were measured using. Pyrite analysis were conducted for the top and sub sediments. The density, establishment, and the relative root growth of Rhizophoraceae were also determined. Free tidal inundation at abandoned pond sites improved the sediment quality. The high density, establishment, and growth of mangrove seedlings were characterized by freely drained areas with a higher pH (field and oxidisable), lower organic content, and high proportion of silt/clay. Higher density and growth also correlated to reduced environments. Sulfur species did not influence the density, establishment, and growth of the seedlings directly. The supply of propagules from the mangrove stands, or access from good waterways were also important for seedlings to establish naturally.

Effects of rice husk biochar and sugar-mill by-products on methane consumption from two different soils

Methane (CH4) is an important greenhouse gas with a global warming potential (GWP) 25 times greater than carbon dioxide (CO2) that can be produced or consumed in soils depending on environmental conditions and other factors. Biochar application to soils has been shown to reduce CH4 emissions and to increase CH4 consumption. However, the effects of rice husk biochar (RB) have not been thoroughly investigated. Two 60-day laboratory incubation experiments were conducted to investigate the effects of amending two soil types with RB, raw mill mud (MM) and composted mill mud (CM) on soil CH4 consumption and emissions. Soil cores incubated in 1 L glass jars and gas samples were analysed for CH4 using gas chromatography. Average CH4 consumption rates varied from -0.06 to -0.68 g CH4-C( )/ha/d in sandy loam soil and -0.59 to -1.00 g CH4-C/ha/d in clay soil. Application of RB resulted in CH4 uptake of -0.52 to -0.55 g CH4-C/ha/d in sandy loam and -0.76 to -0.91 g CH4-C/ha/d in clay soil. Addition of MM showed low CH4 emissions or consumption at 60% water-filled pore space (WFPS) in both soils. However, at high water contents (>75% WFPS) the application of MM produced high rates of CH4 emissions which were significantly suppressed when RB was added. Cumulative emissions of the MM treatment produced 108.9 g CH4-C/ha at 75% WFPS and 11 459.3 g CH4-C/ha at 90% WFPS in sandy loam soil over a period of 60 days. RB can increase CH4 uptake under low soil water content (SWC) and decrease CH4 emissions under anaerobic conditions. CM expressed more potential to reduce CH4 emissions than those of MM.

Increasing NH4+ retention efficiency of sandy soils by adding a modified kaolin material: Mesolite

Water and ammonium retention by sandy soils may be low and result in leaching of applied fertiliser. To increase water and nutrient retention, zeolite is sometimes applied as a soil ameliorant for high value land uses including turf and horticulture. We have used a new modified kaolin material (MesoLite) as a soil amendment to test the efficiency of NH4+ retention and compared the results with natural zeolite. MesoLite is made by caustic reaction of kaolin at temperature between 80-95°C; although it has a moderate surface area, its cation exchange capacity is very high;(SA=13m2/g,CEC=500meq/100g). A 13cm tall sand column filled with ~450g of sandy soil homogeneously mixed with 1, 2, 4, and 8g of MesoLite or natural zeolite per 1kg of soil was prepared. After saturation with local bore water, concentrated ammonium sulfate solution was injected at the base. Then, bore water was passed from bottom to top through the column at amounts up to 6 pore volumes and at a constant flow rate of 10ml/min using a peristaltic pump. Concentrations of leached NH4+ were determined using an AutoAnalyser. The concentration of NH4+ leached from the column with 0.4% MesoLite was greatly (90%) reduced relative to unamended soil. Under these conditions NH4+ retention by the soil-MesoLite mixture was 11.5 times more efficient than the equivalent soil-natural zeolite mixture. Glasshouse experiments conducted in a separate study show that NH4+ adsorbed by MesoLite is available to plants.

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.

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 of the dry matter and C/N ratio of the OAs, WFPS, and the soil CEC. Application of RB reduced N2O emissions by as much as 42-64% depending on WFPS. The reductions in both CO2 and N2O emissions after application of RB were due to a reduced bioavailability of C and not immobilisation of N. These findings show that the effect of OAs on soil GHG emissions can vary substantially depending on their chemical properties. OAs with a high availability of labile C and N can lead to elevated emissions of CO2 and N2O, while rice husk biochar showed potential in reducing overall soil GHG emissions.

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.

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.

Direct shear testing on unsaturated silty soils to investigate the effects of drying and wetting on shear strength parameters at low suction

A modified conventional direct shear device was used to measure unsaturated shear strength of two silty soils at low suction values (0 ~ 50 kPa) that were achieved by following drying and wetting paths of soil water characteristic curves (SWCCs). The results revealed that the internal friction angle of the soils was not significantly affected by either the suction or the drying wetting SWCCs. The apparent cohesion of soil increased with a decreasing rate as suction increased. Shear stress-shear displacement curves obtained from soil specimens subjected to the same net normal stress and different suction values showed a higher initial stiffness and a greater peak stress as suction increased. A soil in wetting exhibited slightly higher peak shear stress and more contractive volume change behavior than that of soil in drying at the same net normal stress and suction.

Efficient simulation of unsaturated flow in heterogeneous soils using an exponential integrator

- Provided a practical variable-stepsize implementation of the exponential Euler method (EEM). - Introduced a new second-order variant of the scheme that enables the local error to be estimated at the cost of a single additional function evaluation. - New EEM implementation outperformed sophisticated implementations of the backward differentiation formulae (BDF) of order 2 and was competitive with BDF of order 5 for moderate to high tolerances.