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.

Structure and dynamics of a gene network model incorporating small RNAs

As advances in molecular biology continue to reveal additional layers of complexity in gene regulation, computational models need to incorporate additional features to explore the implications of new theories and hypotheses. It has recently been suggested that eukaryotic organisms owe their phenotypic complexity and diversity to the exploitation of small RNAs as signalling molecules. Previous models of genetic systems are, for several reasons, inadequate to investigate this theory. In this study, we present an artificial genome model of genetic regulatory networks based upon previous work by Torsten Reil, and demonstrate how this model generates networks with biologically plausible structural and dynamic properties. We also extend the model to explore the implications of incorporating regulation by small RNA molecules in a gene network. We demonstrate how, using these signals, highly connected networks can display dynamics that are more stable than expected given their level of connectivity.

Dynamics of gene expression in an artificial genome

Complex systems techniques provide a powerful tool to study the emergent properties of networks of interacting genes. In this study we extract models of genetic regulatory networks from an artificial genome, represented by a sequence of nucleotides, and analyse how variations in the connectivity and degree of inhibition of the extracted networks affects the resulting classes of behaviours. For low connectivity systems were found to be very stable. Only with higher connectivity was a significant occurrence of chaos found. Most interestingly, the peak in occurrence of chaos occurs perched on the edge of a phase transition in the occurrence of attractors.

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.