Mineralogy and Magnetic Properties of Basaltic Substrate Soils: Kaho'olawe and Big Island, Hawaii

Magnetic behavior of soils can seriously hamper the performance of geophysical sensors. Currently, we have little understanding of the types of minerals responsible for the magnetic behavior, as well as their distribution in space and evolution through time. This study investigated the magnetic characteristics and mineralogy of Fe-rich soils developed on basaltic substrate in Hawaii. We measured the spatial distribution of magnetic susceptibility (χlf) and frequency dependence (χfd%) across three test areas in a well-developed eroded soil on Kaho'olawe and in two young soils on the Big Island of Hawaii. X-ray diffraction spectroscopy, x-ray fluorescence spectroscopy (XFCF), chemical dissolution, thermal analysis, and temperature-dependent magnetic studies were used to characterize soil development and mineralogy for samples from soil pits on Kaho'olawe, surface samples from all three test areas, and unweathered basalt from the Big Island of Hawaii. The measurements show a general increase in magnetic properties with increasing soil development. The XRF Fe data ranged from 13% for fresh basalt and young soils on the Big Island to 58% for material from the B horizon of Kaho'olawe soils. Dithionite-extractable and oxalate-extractable Fe percentages increase with soil development and correlate with χlf-and χfd%, respectively. Results from the temperature-dependent susceptibility measurements show that the high soil magnetic properties observed in geophysical surveys in Kaho'olawe are entirely due to neoformed minerals. The results of our studies have implications for the existing soil survey of Kaho'olawe and help identify methods to characterize magnetic minerals in tropical soils.

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.

Anthropogenic nitrogen deposition predicts local grassland primary production worldwide

Humans dominate many important Earth system processes including the nitrogen (N) cycle. Atmospheric N deposition affects fundamental processes such as carbon cycling, climate regulation, and biodiversity, and could result in changes to fundamental Earth system processes such as primary production. Both modelling and experimentation have suggested a role for anthropogenically altered N deposition in increasing productivity, nevertheless, current understanding of the relative strength of N deposition with respect to other controls on production such as edaphic conditions and climate is limited. Here we use an international multiscale data set to show that atmospheric N deposition is positively correlated to aboveground net primary production (ANPP) observed at the 1-m2 level across a wide range of herbaceous ecosystems. N deposition was a better predictor than climatic drivers and local soil conditions, explaining 16% of observed variation in ANPP globally with an increase of 1 kg N·ha-1·yr-1 increasing ANPP by 3%. Soil pH explained 8% of observed variation in ANPP while climatic drivers showed no significant relationship. Our results illustrate that the incorporation of global N deposition patterns in Earth system models are likely to substantially improve estimates of primary production in herbaceous systems. In herbaceous systems across the world, humans appear to be partially driving local ANPP through impacts on the N cycle.

Mathematical modeling for improved greenhouse gas balances, agro-ecosystems, and policy development: Lessons from the Australian experience

If the land sector is to make significant contributions to mitigating anthropogenic greenhouse gas (GHG) emissions in coming decades, it must do so while concurrently expanding production of food and fiber. In our view, mathematical modeling will be required to provide scientific guidance to meet this challenge. In order to be useful in GHG mitigation policy measures, models must simultaneously meet scientific, software engineering, and human capacity requirements. They can be used to understand GHG fluxes, to evaluate proposed GHG mitigation actions, and to predict and monitor the effects of specific actions; the latter applications require a change in mindset that has parallels with the shift from research modeling to decision support. We compare and contrast 6 agro-ecosystem models (FullCAM, DayCent, DNDC, APSIM, WNMM, and AgMod), chosen because they are used in Australian agriculture and forestry. Underlying structural similarities in the representations of carbon flows though plants and soils in these models are complemented by a diverse range of emphases and approaches to the subprocesses within the agro-ecosystem. None of these agro-ecosystem models handles all land sector GHG fluxes, and considerable model-based uncertainty exists for soil C fluxes and enteric methane emissions. The models also show diverse approaches to the initialisation of model simulations, software implementation, distribution, licensing, and software quality assurance; each of these will differentially affect their usefulness for policy-driven GHG mitigation prediction and monitoring. Specific requirements imposed on the use of models by Australian mitigation policy settings are discussed, and areas for further scientific development of agro-ecosystem models for use in GHG mitigation policy are proposed.

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.

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.

Critical Reappraisal Of Plasticity Index Of Soils

A fundamental approach, based on Gouy-Chapman theory of double layer, has been provided to micromechanistically interpret the plasticity index of soils and their relationship with liquid limit. The relationships between plasticity index and liquid limit, developed earlier, through statistical approaches and critical state concepts, have been reexamined. The statistical analysis of extensive published data has resulted in the relationship, IP = 0.74 (wL - 8). On comparison with other relationships in vogue the proposed equation has been found to give better agreement. From the reappraisal of critical state approaches consistent with the micromechanistic interpretation, the possible range of parameters have been computed and compared with those obtained by statistical means to enhance the credibility of the proposed relationship.

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.

Association of a second phase of mortality in cucumber seedlings with a rapid rate of metalaxyl biodegradation in greenhouse soils

A study was undertaken from 2004 to 2007 to investigate factors associated with decreased efficacy of metalaxyl to manage damping-off of cucumber in Oman. A survey over six growing seasons showed that growers lost up to 14.6% of seedlings following application of metalaxyl. No resistance to metalaxyl was found among Pythium isolates. Damping-off disease in the surveyed greenhouses followed two patterns. In most (69%) greenhouses, seedling mortality was found to occur shortly after transplanting and decrease thereafter (Phase-I). However, a second phase of seedling mortality (Phase-II) appeared 9-14 d after transplanting in about 31% of the surveyed greenhouses. Analysis of the rate of biodegradation of metalaxyl in six greenhouses indicated a significant increase in the rate of metalaxyl biodegradation in greenhouses, which encountered Phase-II damping-off. The half-life of metalaxyl dropped from 93 d in soil, which received no previous metalaxyl treatment to 14 d in soil, which received metalaxyl for eight consecutive seasons, indicating an enhanced rate of metalaxyl biodegradation after repeated use. Multiple applications of metalaxyl helped reduce the appearance of Phase-II damping-off. This appears to be the first report of rapid biodegradation of metalaxyl in greenhouse soils and the first report of its association with appearance of a second phase of mortality in cucumber seedlings.

Detection of Arcobacter spp. in piggery effluent and effluent-irrigated soils in southeast Queensland

Aims: To investigate the occurrence and levels of Arcobacter spp. in pig effluent ponds and effluent-treated soil. Methods and Results: A Most Probable Number (MPN) method was developed to assess the levels of Arcobacter spp. in seven pig effluent ponds and six effluent-treated soils, immediately after effluent irrigation. Arcobacter spp. levels in the effluent ponds varied from 6.5 × 105 to 1.1 × 108 MPN 100 ml-1 and in freshly irrigated soils from 9.5 × 102 to 2.8 × 104 MPN g-1 in all piggery environments tested. Eighty-three Arcobacter isolates were subjected to an abbreviated phenotypic test scheme and examined using a multiplex polymerase chain reaction (PCR). The PCR identified 35% of these isolates as Arcobacter butzleri, 49% as Arcobacter cryaerophilus while 16% gave no band. All 13 nonreactive isolates were subjected to partial 16S rDNA sequencing and showed a high similarity (>99%) to Arcobacter cibarius. Conclusions: A. butzleri, A. cryaerophilus and A. cibarius were isolated from both piggery effluent and effluent-irrigated soil, at levels suggestive of good survival in the effluent pond. Significance and Impact of the Study: This is the first study to provide quantitative information on Arcobacter spp. levels in piggery effluent and to associate A. cibarius with pigs and piggery effluent environments.

Compressibility of Partly Saturated Soils

From the considerations of the truncated diffuse double layer due to partial saturation, an equilibrium equation is written in terms of the soil state parameter viz. void ratio, e, void ratio at liquid limit, e=wG, degree of saturation, S and external applied stress, p. The type of fine grained soil is subdued by its liquid limit. The viability of the formulated approach has been experimentally verified. This approach, essentially eliminates the difficulty in the measurement of pore air pressure, u pore water pressure, u and computation of \Gv parameter.

Subsoil constraints to grain production in the cropping soils of the north-eastern region of Australia: an overview

In dryland agricultural systems of the subtropical, semi-arid region of north-eastern Australia, water is the most limiting resource. Crop productivity depends on the efficient use of rainfall and available water stored in the soil during fallow. Agronomic management practices including a period of fallow, stubble retention, and reduced tillage enhance reserves of soil water. However, access to stored water in these soils may be restricted by the presence of growth-limiting conditions in the rooting zone of the crop. These have been termed as subsoil constraints. Subsoil constraints may include compacted or gravel layers (physical), sodicity, salinity, acidity, nutrient deficiencies, presence of toxic elements (chemical) and low microbial activity (biological). Several of these constraints may occur together in some soils. Farmers have often not been able to obtain the potential yield determined by their prevailing climatic conditions in the marginal rainfall areas of the northern grains region. In the past, the adoption of soil management practices had been largely restricted to the top 100 mm soil layer. Exploitation of the subsoil as a source of water and nutrients has largely been overlooked. The key towards realising potential yields would be to gain better understanding of subsoils and their limitations, then develop options to manage them practically and economically. Due to the complex nature of the causal factors of these constraints, efforts are required for a combination of management approaches rather than individual options, with the aim to combat these constraints for sustainable crop production, managing natural resources and avoiding environmental damage.

Impact of subsoil constraints on wheat yield and gross margin on fine-textured soils of the southern Victorian Mallee

The APSIM-Wheat module was used to investigate our present capacity to simulate wheat yields in a semi-arid region of eastern Australia (the Victorian Mallee), where hostile subsoils associated with salinity, sodicity, and boron toxicity are known to limit grain yield. In this study we tested whether the effects of subsoil constraints on wheat growth and production could be modelled with APSIM-Wheat by assuming that either: (a) root exploration within a particular soil layer was reduced by the presence of toxic concentrations of salts, or (b) soil water uptake from a particular soil layer was reduced by high concentration of salts through osmotic effects. After evaluating the improved predictive capacity of the model we applied it to study the interactions between subsoil constraints and seasonal conditions, and to estimate the economic effect that subsoil constraints have on wheat farming in the Victorian Mallee under different climatic scenarios. Although the soils had high levels of salinity, sodicity, and boron, the observed variability in root abundance at different soil layers was mainly related to soil salinity. We concluded that: (i) whether the effect of subsoil limitations on growth and yield of wheat in the Victorian Mallee is driven by toxic, osmotic, or both effects acting simultaneously still requires further research, (ii) at present, the performance of APSIM-Wheat in the region can be improved either by assuming increased values of lower limit for soil water extraction, or by modifying the pattern of root exploration in the soil pro. le, both as a function of soil salinity. The effect of subsoil constraints on wheat yield and gross margin can be expected to be higher during drier than wetter seasons. In this region the interaction between climate and soil properties makes rainfall information alone, of little use for risk management and farm planning when not integrated with cropping systems models.

Models for the field-based toxicity of copper and zinc salts to wheat in 11 Australian soils and comparison to laboratory-based models

Laboratory-based relationships that model the phytotoxicity of metals using soil properties have been developed. This paper presents the first field-based phytotoxicity relationships. Wheat(Triticum aestivum L) was grown at 11 Australian field sites at which soil was spiked with copper (Cu) and zinc (Zn) salts. Toxicity was measured as inhibition of plant growth at 8 weeks and grain yield at harvest. The added Cu and Zn EC10 values for both endpoints ranged from approximately 3 to 4760 mg/kg. There were no relationships between field-based 8-week biomass and grain yield toxicity values for either metal. Cu toxicity was best modelled using pH and organic carbon content while Zn toxicity was best modelled using pH and the cation exchange capacity. The best relationships estimated toxicity within a factor of two of measured values. Laboratory-based phytotoxicity relationships could not accurately predict field-based phytotoxicity responses.

Long-term trends in fertility of soils under continuous cultivation and cereal cropping in southern Queensland. IX. Simulation of soil carbon and nitrogen pools using CENTURY model

Cultivation and cropping of soils results in a decline in soil organic carbon and soil nitrogen, and can lead to reduced crop yields. The CENTURY model was used to simulate the effects of continuous cultivation and cereal cropping on total soil organic matter (C and N), carbon pools, nitrogen mineralisation, and crop yield from 6 locations in southern Queensland. The model was calibrated for each replicate from the original datasets, allowing comparisons for each replicate rather than site averages. The CENTURY model was able to satisfactorily predict the impact of long-term cultivation and cereal cropping on total organic carbon, but was less successful in simulating the different fractions and nitrogen mineralisation. The model firstly over-predicted the initial (pre-cropping) soil carbon and nitrogen concentration of the sites. To account for the unique shrinking and swelling characteristics of the Vertosol soils, the default annual decomposition rates of the slow and passive carbon pools were doubled, and then the model accurately predicted initial conditions. The ability of the model to predict carbon pool fractions varied, demonstrating the difficulty inherent in predicting the size of these conceptual pools. The strength of the model lies in the ability to closely predict the starting soil organic matter conditions, and the ability to predict the impact of clearing, cultivation, fertiliser application, and continuous cropping on total soil carbon and nitrogen.