Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

A GIS Based Assessment of Land Suitable for Growing Hoop Pine in the Atherton, Eacham and Herberton Shires of North Queensland

The area of private land suitable and available for growing hoop pine (Araucaria cunninghamii) on the Atherton Tablelands in North Queensland was modelled using a geographic information system (GIS). In Atherton, Eacham and Herberton shires, approximately 64,700 ha of privately owned land were identified as having a mean annual rainfall and soil type similar to Forestry Plantations Queensland (FPQ) hoop pine growth plots with an approximate growth rate of 20 m3 per annum. Land with slope of over 25° and land covered with native vegetation were excluded in the estimation. If land which is currently used for high-value agriculture is also excluded, the net area of land potentially suitable and available for expansion of hoop pine plantations is approximately 22,900 ha. Expert silvicultural advice emphasized the role of site preparation and weed control in affecting the long-term growth rate of hoop pine. Hence, sites with less than optimal fertility and rainfall may be considered as being potentially suitable for growing hoop pine at a lower growth rate. The datasets had been prepared at various scales and differing precision for their description of land attributes. Therefore, the results of this investigation have limited applicability for planning at the individual farm level but are useful at the regional level to target areas for plantation expansion.

Effect of ionic strength and clay mineralogy on Na-Ca exchange and the SAR-ESP relationship

The relationship between sodium adsorption ratio (SAR) and exchangeable sodium percentage (ESP) for all soils has traditionally been assumed to be similar to that developed by the United States Salinity Laboratory (USSL) in 1954. However, under certain conditions, this relationship has been shown not to be constant, but to vary with both ionic strength and clay mineralogy. We conducted a detailed experiment to determine the effect of ionic strength on the Na+-Ca2+ exchange of four clay minerals (kaolinite, illite, pyrophyllite, and montmorillonite), with results related to the diffuse double-layer (DDL) model. Clays in which external exchange sites dominated (kaolinite and pyrophyllite) tended to show an overall preference for Na+, with the magnitude of this preference increasing with decreasing ESP. For these external surfaces, increases in ionic strength were found to increase preference for Na+. Although illite (2:1 non-expanding mineral) was expected to be dominated by external surfaces, this clay displayed an overall preference for Ca2+, possibly indicating the opening of quasicrystals and the formation of internal exchange surfaces. For the expanding 2:1 clay, montmorillonite, Na+-Ca2+ exchange varied due to the formation of quasicrystals (and internal exchange surfaces) from individual clay platelets. At small ionic strength and large ESP, the clay platelets dispersed and were dominated by external exchange surfaces (displaying preference for Na+). However, as ionic strength increased and ESP decreased, quasicrystals (and internal exchange surfaces) formed, and preference for Ca2+ increased. Therefore, the relationship between SAR and ESP is not constant and should be determined directly for the soil of interest.

A review of the use of the basic cation saturation ratio and the "ideal" soil

The use of 'balanced' Ca, Mg, and K ratios, as prescribed by the basic cation saturation ratio (BCSR) concept, is still used by some private soil-testing laboratories for the interpretation of soil analytical data. This review aims to examine the suitability of the BCSR concept as a method for the interpretation of soil analytical data. According to the BCSR concept, maximum plant growth will be achieved only when the soil’s exchangeable Ca, Mg, and K concentrations are approximately 65 % Ca, 10 % Mg, and 5 % K (termed the ‘ideal soil’). This ‘ideal soil’ was originally proposed by Firman Bear and co-workers in New Jersey (USA) during the 1940s as a method of reducing luxury K uptake by alfalfa (Medicago sativa L.). At about the same time, William Albrecht, working in Missouri (USA), concluded through his own investigations that plants require a soil with a high Ca saturation for optimal growth. Whilst it now appears that several of Albrecht’s experiments were fundamentally flawed, the BCSR (‘balanced soil’) concept has been widely promoted, suggesting that the prescribed cationic ratios provide optimum chemical, physical, and biological soil properties. Our examination of data from numerous studies (particularly those of Albrecht and Bear, themselves) would suggest that, within the ranges commonly found in soils, the chemical, physical, and biological fertility of a soil is generally not influenced by the ratios of Ca, Mg, and K. The data do not support the claims of the BCSR, and continued promotion of the BCSR will result in the inefficient use of resources in agriculture and horticulture.

Parliament House Competition Model, Canberra

View of model for competition entry.

How to unravel and solve soil fertility problems

This book provides a way for farmers in developing countries to benefit from scientific knowledge on plant nutrition and soil fertility. Specifically, it will help farmers recognise and deal with shortages or excesses of chemical elements.

Structural analysis and molecular model of a self-incompatibility RNase from wild tomato

Self-incompatibility RNases (S-RNases) are an allelic series of style glycoproteins associated with rejection of self-pollen in solanaceous plants. The nucleotide sequences of S-RNase alleles from several genera have been determined, but the structure of the gene products has only been described for those from Nicotiana alata. We report on the N-glycan structures and the disulfide bonding of the S-3-RNase from wild tomato (Lycopersicon peruvianum) and use this and other information to construct a model of this molecule. The S-3-RNase has a single N-glycosylation site (Asn-28) to which one of three N-glycans is attached. S-3-RNase has seven Cys residues; six are involved in disulfide linkages (Cys-16-Cys-21, Cys-46-Cys-91, and Cys-166-Cys-177), and one has a free thiol group (Cys-150). The disulfide-bonding pattern is consistent with that observed in RNase Rh, a related RNase for which radiographic-crystallographic information is available. A molecular model of the S-3-RNase shows that four of the most variable regions of the S-RNases are clustered on one surface of the molecule. This is discussed in the context of recent experiments that set out to determine the regions of the S-RNase important for recognition during the self-incompatibility response.

Improving wheat simulation capabilities in Australia from a cropping systems perspective - II. Testing simulation capabilities of wheat growth

To simulate cropping systems, crop models must not only give reliable predictions of yield across a wide range of environmental conditions, they must also quantify water and nutrient use well, so that the status of the soil at maturity is a good representation of the starting conditions for the next cropping sequence. To assess the suitability for this task a range of crop models, currently used in Australia, were tested. The models differed in their design objectives, complexity and structure and were (i) tested on diverse, independent data sets from a wide range of environments and (ii) model components were further evaluated with one detailed data set from a semi-arid environment. All models were coded into the cropping systems shell APSIM, which provides a common soil water and nitrogen balance. Crop development was input, thus differences between simulations were caused entirely by difference in simulating crop growth. Under nitrogen non-limiting conditions between 73 and 85% of the observed kernel yield variation across environments was explained by the models. This ranged from 51 to 77% under varying nitrogen supply. Water and nitrogen effects on leaf area index were predicted poorly by all models resulting in erroneous predictions of dry matter accumulation and water use. When measured light interception was used as input, most models improved in their prediction of dry matter and yield. This test highlighted a range of compensating errors in all modelling approaches. Time course and final amount of water extraction was simulated well by two models, while others left up to 25% of potentially available soil water in the profile. Kernel nitrogen percentage was predicted poorly by all models due to its sensitivity to small dry matter changes. Yield and dry matter could be estimated adequately for a range of environmental conditions using the general concepts of radiation use efficiency and transpiration efficiency. However, leaf area and kernel nitrogen dynamics need to be improved to achieve better estimates of water and nitrogen use if such models are to be use to evaluate cropping systems. (C) 1998 Elsevier Science B.V.