Modelling predicts positive and negative interactions between three Australian tropical tree species in monoculture and binary mixture

Computer modelling promises to be an important tool for analysing and predicting interactions between trees within mixed species forest plantations. This study explored the use of an individual-based mechanistic model as a predictive tool for designing mixed species plantations of Australian tropical trees. The `spatially explicit individually based-forest simulator' (SeXI-FS) modelling system was used to describe the spatial interaction of individual tree crowns within a binary mixed-species experiment. The three-dimensional model was developed and verified with field data from three forest tree species grown in tropical Australia. The model predicted the interactions within monocultures and binary mixtures of Flindersia brayleyana, Eucalyptus pellita and Elaeocarpus grandis, accounting for an average of 42% of the growth variation exhibited by species in different treatments. The model requires only structural dimensions and shade tolerance as species parameters. By modelling interactions in existing tree mixtures, the model predicted both increases and reductions in the growth of mixtures (up to +/-50% of stem volume at 7 years) compared to monocultures. This modelling approach may be useful for designing mixed tree plantations.

Simulating the effects of saline and sodic subsoils on wheat crops growing on Vertosols

Soils with high levels of chloride and/or sodium in their subsurface layers are often referred to as having subsoil constraints (SSCs). There is growing evidence that SSCs affect wheat yields by increasing the lower limit of a crop's available soil water (CLL) and thus reducing the soil's plant-available water capacity (PAWC). This proposal was tested by simulation of 33 farmers' paddocks in south-western Queensland and north-western New South Wales. The simulated results accounted for 79% of observed variation in grain yield, with a root mean squared deviation (RMSD) of 0.50 t/ha. This result was as close as any achieved from sites without SSCs, thus providing strong support for the proposed mechanism that SSCs affect wheat yields by increasing the CLL and thus reducing the soil's PAWC. In order to reduce the need to measure CLL of every paddock or management zone, two additional approaches to simulating the effects of SSCs were tested. In the first approach the CLL of soils was predicted from the 0.3-0.5 m soil layer, which was taken as the reference CLL of a soil regardless of its level of SSCs, while the CLL values of soil layers below 0.5 m depth were calculated as a function of these soils' 0.3-0.5 m CLL values as well as of soil depth plus one of the SSC indices EC, Cl, ESP, or Na. The best estimates of subsoil CLL values were obtained when the effects of SSCs were described by an ESP-dependent function. In the second approach, depth-dependent CLL values were also derived from the CLL values of the 0.3-0.5 m soil layer. However, instead of using SSC indices to further modify CLL, the default values of the water-extraction coefficient (kl) of each depth layer were modified as a function of the SSC indices. The strength of this approach was evaluated on the basis of correlation of observed and simulated grain yields. In this approach the best estimates were obtained when the default kl values were multiplied by a Cl-determined function. The kl approach was also evaluated with respect to simulated soil moisture at anthesis and at grain maturity. Results using this approach were highly correlated with soil moisture results obtained from simulations based on the measured CLL values. This research provides strong evidence that the effects of SSCs on wheat yields are accounted for by the effects of these constraints on wheat CLL values. The study also produced two satisfactory methods for simulating the effects of SSCs on CLL and on grain yield. While Cl and ESP proved to be effective indices of SSCs, EC was not effective due to the confounding effect of the presence of gypsum in some of these soils. This study provides the tools necessary for investigating the effects of SSCs on wheat crop yields and natural resource management (NRM) issues such as runoff, recharge, and nutrient loss through simulation studies. It also facilitates investigation of suggested agronomic adaptations to SSCs.

Non-specific conducting polymer-based array capable of monitoring odour emissions from a biofiltration system in a piggery building

A commercial non-specific gas sensor array system was evaluated in terms of its capability to monitor the odour abatement performance of a biofiltration system developed for treating emissions from a commercial piggery building. The biofiltration system was a modular system comprising an inlet ducting system, humidifier and closed-bed biofilter. It also included a gravimetric moisture monitoring and water application system for precise control of moisture content of an organic woodchip medium. Principal component analysis (PCA) of the sensor array measurements indicated that the biofilter outlet air was significantly different to both inlet air of the system and post-humidifier air. Data pre-processing techniques including normalising and outlier handling were applied to improve the odour discrimination performance of the non-specific gas sensor array. To develop an odour quantification model using the sensor array responses of the non-specific sensor array, PCA regression, artificial neural network (ANN) and partial least squares (PLS) modelling techniques were applied. The correlation coefficient (r(2)) values of the PCA, ANN, and PLS models were 0.44, 0.62 and 0.79, respectively.

Ocean dispersal modelling for propagules of pond apple (Annona glabra L.)

Pond apple invades riparian and coastal environments with water acting as the main vector for dispersal. As seeds float and can reach the ocean, a seed tracking model driven by near surface ocean currents was used to develop maps of potential seed dispersal. Seeds were ‘released’ in the model from sites near the mouths of major North Queensland rivers. Most seeds reach land within three months of release, settling predominately on windward-facing locations. During calm and monsoonal conditions, seeds were generally swept in a southerly direction, however movement turns northward during south easterly trade winds. Seeds released in February from the Johnstone River were capable of being moved anywhere from 100 km north to 150 km south depending on prevailing conditions. Although wind driven currents are the primary mechanism influencing seed dispersal, tidal currents, the East Australian Current, and other factors such as coastline orientation, release location and time also play an important role in determining dispersal patterns. In extreme events such as tropical cyclone Justin in 1997, north east coast rivers could potentially transport seed over 1300 km to the Torres Strait, Papua New Guinea and beyond.

Applications of CLIMEX modelling leading to improved biological control

Wilmot Senaratne, Bill Palmer and Bob Sutherst recently published their paper 'Applications of CLIMEX modelling leading to improved biological control' in Proceedings of the 16th Australian Weeds Conference. They looked at three examples where modern climate matching techniques using computer software produces decisions and results than might happen using previous techniques such as climadiagrams. Assessment of climatic suitability is important at various stages of a biological control project; from initial foreign exploration, to risk assessment in preparation for the release of a particular agent, through to selection of release sites that maximise the agent´s chances of initial establishment. It is now also necessary to predict potential future distributions of both target weeds and agents under climate change.

Estimation in a multiplicative mixed model involving a genetic relationship matrix.

Genetic models partitioning additive and non-additive genetic effects for populations tested in replicated multi-environment trials (METs) in a plant breeding program have recently been presented in the literature. For these data, the variance model involves the direct product of a large numerator relationship matrix A, and a complex structure for the genotype by environment interaction effects, generally of a factor analytic (FA) form. With MET data, we expect a high correlation in genotype rankings between environments, leading to non-positive definite covariance matrices. Estimation methods for reduced rank models have been derived for the FA formulation with independent genotypes, and we employ these estimation methods for the more complex case involving the numerator relationship matrix. We examine the performance of differing genetic models for MET data with an embedded pedigree structure, and consider the magnitude of the non-additive variance. The capacity of existing software packages to fit these complex models is largely due to the use of the sparse matrix methodology and the average information algorithm. Here, we present an extension to the standard formulation necessary for estimation with a factor analytic structure across multiple environments.

Effect of timing of a deficit-irrigation allocation on corn evapotranspiration, yield, water use efficiency and dry mass

Water regulations have decreased irrigation water supplies in Nebraska and some other areas of the USA Great Plains. When available water is not enough to meet crop water requirements during the entire growing cycle, it becomes critical to know the proper irrigation timing that would maximize yields and profits. This study evaluated the effect of timing of a deficit-irrigation allocation (150 mm) on crop evapotranspiration (ETc), yield, water use efficiency (WUE = yield/ETc), irrigation water use efficiency (IWUE = yield/irrigation), and dry mass (DM) of corn (Zea mays L.) irrigated with subsurface drip irrigation in the semiarid climate of North Platte, NE. During 2005 and 2006, a total of sixteen irrigation treatments (eight each year) were evaluated, which received different percentages of the water allocation during July, August, and September. During both years, all treatments resulted in no crop stress during the vegetative period and stress during the reproductive stages, which affected ETc, DM, yield, WUE and IWUE. Among treatments, ETc varied by 7.2 and 18.8%; yield by 17 and 33%; WUE by 12 and 22%, and IWUE by 18 and 33% in 2005 and 2006, respectively. Yield and WUE both increased linearly with ETc and with ETc/ETp (ETp = seasonal ETc with no water stress), and WUE increased linearly with yield. The yield response factor (ky) averaged 1.50 over the two seasons. Irrigation timing affected the DM of the plant, grain, and cob, but not that of the stover. It also affected the percent of DM partitioned to the grain (harvest index), which increased linearly with ETc and averaged 56.2% over the two seasons, but did not affect the percent allocated to the cob or stover. Irrigation applied in July had the highest positive coefficient of determination (R2) with yield. This high positive correlation decreased considerably for irrigation applied in August, and became negative for irrigation applied in September. The best positive correlation between the soil water deficit factor (Ks) and yield occurred during weeks 12-14 from crop emergence, during the "milk" and "dough" growth stages. Yield was poorly correlated to stress during weeks 15 and 16, and the correlation became negative after week 17. Dividing the 150 mm allocation about evenly among July, August and September was a good strategy resulting in the highest yields in 2005, but not in 2006. Applying a larger proportion of the allocation in July was a good strategy during both years, and the opposite resulted when applying a large proportion of the allocation in September. The different results obtained between years indicate that flexible irrigation scheduling techniques should be adopted, rather than relying on fixed timing strategies.

Surveillance protocols for management of invasive plants: modelling Chilean needle grass (Nassella neesiana) in Australia

Aim: To develop a surveillance support model that enables prediction of areas susceptible to invasion, comparative analysis of surveillance methods and intensity and assessment of eradication feasibility. To apply the model to identify surveillance protocols for generalized invasion scenarios and for evaluating surveillance and control for a context-specific plant invasion. Location: Australia. Methods: We integrate a spatially explicit simulation model, including plant demography and dispersal vectors, within a Geographical Information System. We use the model to identify effective surveillance protocols using simulations of generalized plant life-forms spreading via different dispersal mechanisms in real landscapes. We then parameterize the surveillance support model for Chilean needle grass [CNG; Nassella neesiana (Trin. & Rupr.) Barkworth], a highly invasive tussock grass, which is an eradication target in south-eastern Queensland, Australia. Results: General surveillance protocols that can guide rapid response surveillance were identified; suitable habitat that is susceptible to invasion through particular dispersal syndromes should be targeted for surveillance using an adaptive seek-and-destroy method. The search radius of the adaptive method should be based on maximum expected dispersal distances. Protocols were used to define a surveillance strategy for CNG, but simulations indicated that despite effective and targeted surveillance, eradication is implausible at current intensities. Main conclusions: Several important surveillance protocols emerged and simulations indicated that effectiveness can be increased if they are followed in rapid response surveillance. If sufficient data are available, the surveillance support model should be parameterized to target areas susceptible to invasion and determine whether surveillance is effective and eradication is feasible. We discovered that for CNG, regardless of a carefully designed surveillance strategy, eradication is implausible at current intensities of surveillance and control and these efforts should be doubled if they are to be successful. This is crucial information in the face of environmentally and economically damaging invasive species and large, expensive and potentially ineffective control programmes.

Improving the agent selection, release and evaluation process: the role of bioclimatic modelling

The value of CLIMEX models to inform biocontrol programs was assessed, including predicting the potential distribution of biocontrol agents and their subsequent population dynamics, using bioclimatic models for the weed Parkinsonia aculeata, two Lantana camara biocontrol agents, and five Mimosa pigra biocontrol agents. The results showed the contribution of data types to CLIMEX models and the capacity of these models to inform and improve the selection, release and post release evaluation of biocontrol agents. Foremost among these was the quality of spatial and temporal information as well as the extent to which overseas range data samples the species’ climatic envelope. Post hoc evaluation and refinement of these models requires improved long-term monitoring of introduced agents and their dynamics at well selected study sites. The authors described the findings of these case studies, highlighted their implications, and considered how to incorporate models effectively into biocontrol programs.

Modelling management options for management of feral camels in central Australia

Camels (Camelus dromedarius) were introduced into Australia from the 1840s to the early 1900s for transport and hauling cargo in arid regions. Feral populations remained small until the 1930s when many were released after they were superseded for transport by trucks and rail. Although camels have a relatively slow population growth (<10% per annum), the population has not reached carrying capacity and therefore, requires control to reduce the increasing impacts on central Australia. The model developed for the Northern Territory suggested that currently there are insufficient numbers being removed. The model also investigated which control options would have greatest impacts and found harvesting to be most important. The extent to which commercial harvesting can feasibly reduce camel populations requires further analysis. Due to the wide dispersal of camels in Australia, fertility control, even if technically feasible, will not target adults, the most important age class of the population. Habitat preferences were also investigated in the model but more validation is required as the population is still under range expansion. Immediate action is suggested to alleviate future costs as camel populations and their impacts rise.

Predicting risk and benefit a priori in biological control of invasive plant species: a systems modelling approach

While the method using specialist herbivores in managing invasive plants (classical biological control) is regarded as relatively safe and cost-effective in comparison to other methods of management, the rarity of strict monophagy among insect herbivores illustrates that, like any management option, biological control is not risk-free. The challenge for classical biological control is therefore to predict risks and benefits a priori. In this study we develop a simulation model that may aid in this process. We use this model to predict the risks and benefits of introducing the chrysomelid beetle Charidotis auroguttata to manage the invasive liana Macfadyena unguis-cati in Australia. Preliminary host-specificity testing of this herbivore indicated that there was limited feeding on a non-target plant, although the non-target was only able to sustain some transitions of the life cycle of the herbivore. The model includes herbivore, target and non-target life history and incorporates spillover dynamics of populations of this herbivore from the target to the non-target under a variety of scenarios. Data from studies of this herbivore in the native range and under quarantine were used to parameterize the model and predict the relative risks and benefits of this herbivore when the target and non-target plants co-occur. Key model outputs include population dynamics on target (apparent benefit) and non-target (apparent risk) and fitness consequences to the target (actual benefit) and non-target plant (actual risk) of herbivore damage. The model predicted that risk to the non-target became unacceptable (i.e. significant negative effects on fitness) when the ratio of target to non-target in a given patch ranged from 1:1 to 3:2. By comparing the current known distribution of the non-target and the predicted distribution of the target we were able to identify regions in Australia where the agent may be pose an unacceptable risk. By considering risk and benefit simultaneously, we highlight how such a simulation modelling approach can assist scientists and regulators in making more objective decisions a priori, on the value of releasing specialist herbivores as biological control agents.

Use of CLIMEX modelling to identify prospective areas for exploration to find new biological control agents for prickly acacia

Selection of biocontrol agents that are adapted to the climates in areas of intended release demands a thorough analysis of the climates of the source and release sites. We present a case study that demonstrates how use of the CLIMEX software can improve decision making in relation to the identification of prospective areas for exploration for agents to control the woody weed, prickly acacia Acacia nilotica ssp. indica in the arid areas of north Queensland.

Analysis and Modelling of the Effects of Water Stress on Maize Growth and Yield in Dryland Conditions.

It is essential to provide experimental evidence and reliable predictions of the effects of water stress on crop production in the drier, less predictable environments. A field experiment undertaken in southeast Queensland, Australia with three water regimes (fully irrigated, rainfed and irrigated until late canopy expansion followed by rainfed) was used to compare effects of water stress on crop production in two maize (Zea mays L.) cultivars (Pioneer 34N43 and Pioneer 31H50). Water stress affected growth and yield more in Pioneer 34N43 than in Pioneer 31H50. A crop model APSIM-Maize, after having been calibrated for the two cultivars, was used to simulate maize growth and development under water stress. The predictions on leaf area index (LAI) dynamics, biomass growth and grain yield under rain fed and irrigated followed by rain fed treatments was reasonable, indicating that stress indices used by APSIM-Maize produced appropriate adjustments to crop growth and development in response to water stress. This study shows that Pioneer 31H50 is less sensitive to water stress and thus a preferred cultivar in dryland conditions, and that it is feasible to provide sound predictions and risk assessment for crop production in drier, more variable conditions using the APSIM-Maize model.

Modelling the distribution and relative abundance of feral camels in the Northern Territory using count data

The objectives of this study were to predict the potential distribution, relative abundance and probability of habitat use by feral camels in southern Northern Territory. Aerial survey data were used to model habitat association. The characteristics of ‘used’ (where camels were observed) v. ‘unused’ (pseudo-absence) sites were compared. Habitat association and abundance were modelled using generalised additive model (GAM) methods. The models predicted habitat suitability and the relative abundance of camels in southern Northern Territory. The habitat suitability maps derived in the present study indicate that camels have suitable habitat in most areas of southern Northern Territory. The index of abundance model identified areas of relatively high camel abundance. Identifying preferred habitats and areas of high abundance can help focus control efforts.

Best use modelling for sustainable Australia sports field surfaces

The Best Use Modelling for Sustainable Australian Sports Field Surfaces project has achieved significant success. The project has attracted participation from councils throughout Australia, with in excess of 300 sports fields evaluated from 18 councils to date. An important project component is the derivation of a recommended standard procedure for specifying the performance of playing surfaces. An associated step has been to establish recommended playing surface performance standards for community level sports fields. The derived modelling also provides information on the expected usage and associated costs of different sports surface development options. This is expected to assist the Australian turf production industry through demonstrating to councils that cost effective natural turf options exist that can meet higher usage expectation (as a viable alternative to synthetic turf). A web-accessed data base system will be made available to councils from January 2010 on (reference to www.passturf.com). This system will enable participating councils to record and analyse field performance over time. The system is considered world-leading, and will help keep the Australian parks industry to the international forefront. Tools developed as part of the project offer councils the opportunity to internally assess the performance of their current sports field provision, to identify any deficiencies and to determine the best corrective measure if any deficiency is identified. This is expected to offer community benefits to both sports facility providers and facility user groups. In turn this will aid the provision of affordable community access to safe and good quality playing surfaces. Tools and associated information material will be made available to councils throughout Australia by the end of this year, via the Parks and Leisure Aust. web site. The Best Use Modelling Project is work in progress. On-going input will be needed to ensure the web-accessed database software is as user friendly as possible, new performance testing data will need to be inputted, and tools provided to participating councils updated. Through the support of HAL there is now a well-structured, nationally-supported system in place for benchmarking playing surfaces and for assisting councils to optimise their resource allocation to sports field upgrade or maintenance work.