Organic Matter Dynamics in the Wet Tropics Under Different Sugarcane Residue Management Regimes

Retention of sugarcane leaves and tops on the soil surface after harvesting has almost completely replaced pre- and post-harvest burning of crop residues in the Australian sugar industry. Since its introduction around 25 years ago, residue retention has increased soil organic matter to improve soil fertility as well as improve harvest flexibility and reduce erosion. However, in the wet tropics residue retention also poses potential problems of prolonged waterlogging, and late-season release of nitrogen which can reduce sugar content of the crop. The objective of this project is to examine the management of sugarcane residues in the wet tropics using a systems approach. Subsidiary objectives are (a) to improve understanding of nitrogen cycling in Australian sugarcane soils in the wet tropics, and (b) to identify ways to manage crop residues to retain their advantages and limit their disadvantages. Project objectives will be addressed using several approaches. Historic farm production data recorded by sugar mills in the wet tropics will be analysed to determine the effect of residue burning or retention on crop yield and sugar content. The impact of climate on soil processes will be highlighed by development of an index of nitrogen mineralisation using the Agricultural Production Systems Simulator (APSIM) model. Increased understanding of nitrogen cycling in Australian sugarcane soils and management of crop residues will be gained through a field experiment recently established in the Australian wet tropics. From this experiment the decomposition and nitrogen dynamics of residues placed on the soil surface and incorporated will be compared. The effect of differences in temperature, soil water content and pH will be further examined on these soils under glasshouse conditions. Preliminary results show a high ammonium to nitrate ratio in tropics soils, which may be due to low rates of nitrification that increase the retention of nitrogen in a form (ammonium) that is less subject to leaching. Further results will be presented at Congress.

An anisotropic viscous representation of Mohr-Coulomb failure for use is modelling coupled mantle-continent dynamics

In mantle convection models it has become common to make use of a modified (pressure sensitive, Boussinesq) von Mises yield criterion to limit the maximum stress the lithosphere can support. This approach allows the viscous, cool thermal boundary layer to deform in a relatively plate-like mode even in a fully Eulerian representation. In large-scale models with embedded continental crust where the mobile boundary layer represents the oceanic lithosphere, the von Mises yield criterion for the oceans ensures that the continents experience a realistic broad-scale stress regime. In detailed models of crustal deformation it is, however, more appropriate to choose a Mohr-Coulomb yield criterion based upon the idea that frictional slip occurs on whichever one of many randomly oriented planes happens to be favorably oriented with respect to the stress field. As coupled crust/mantle models become more sophisticated it is important to be able to use whichever failure model is appropriate to a given part of the system. We have therefore developed a way to represent Mohr-Coulomb failure within a code which is suited to mantle convection problems coupled to large-scale crustal deformation. Our approach uses an orthotropic viscous rheology (a different viscosity for pure shear to that for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation - neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane which is formally equivalent in any symmetric formulation. It is also important to allow strain-weakening of the material. The material should remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-Integration-point finite element code and will show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion, and comparisons between mantle convection models using the von Mises versus the Mohr-Coulomb yield criteria. The formulation itself is general and applies to 2D and 3D problems, although it is somewhat more complicated to identify the slip plane in 3D.