Simulating crop-parasitic weed interactions using APSIM: Model evaluation and application

The parasitic weed Orobanche crenata inflicts major damage on faba bean, lentil, pea and other crops in Mediterranean environments. The development of methods to control O. crenata is to a large extent hampered by the complexity of host-parasite systems. Using a model of host-parasite interactions can help to explain and understand this intricacy. This paper reports on the evaluation and application of a model simulating host-parasite competition as affected by environment and management that was implemented in the framework of the Agricultural Production Systems Simulator (APSIM). Model-predicted faba bean and O. crenata growth and development were evaluated against independent data. The APSIM-Fababean and -Parasite modules displayed a good capability to reproduce effects of pedoclimatic conditions, faba bean sowing date and O. crenata infestation on host-parasite competition. The r(2) values throughout exceeded 0.84 (RMSD: 5.36 days) for phenological, 0.85 (RMSD: 223.00 g m(-2)) for host growth and 0.78 (RMSD: 99.82 g m(-2)) for parasite growth parameters. Inaccuracies of simulated faba bean root growth that caused some bias of predicted parasite number and host yield loss may be dealt with by more flexibly simulating vertical root distribution. The model was applied in simulation experiments to determine optimum sowing windows for infected and non-infected faba bean in Mediterranean environments. Simulation results proved realistic and testified to the capability of APSIM to contribute to the development of tactical approaches in parasitic weed control.

The introduction and release of Chiasmia inconspicua and C. assimilis (Lepidoptera: Geometridae) for the biological control of Acacia nilotica in Australia

Two geometrid moths Chiasmia inconspicua and Chiasmia assimilis, identified as potential biological control agents for prickly acacia Acacia nilotica subsp. indica, were collected in Kenya and imported into quarantine facilities in Australia where laboratory cultures were established. Aspects of the biologies of both insects were studied and CLIMEX® models indicating the climatically favourable areas of Australia were developed. Host range tests were conducted using an approved test list of 74 plant species and no-choice tests of neonate larvae placed on both cut foliage and potted plants. C. inconspicua developed through to adult on prickly acacia and, in small numbers, Acacia pulchella. C. assimilis developed through to adult on prickly acacia and also in very small numbers on A. pulchella, A. deanei, A. decurrens, and A. mearnsii. In all experiments, the response on prickly acacia could be clearly differentiated from the responses on the non-target species. Both insects were approved for release in Australia. Over a three-year period releases were made at multiple sites in north Queensland, almost all in inland areas. There was no evidence of either insect's establishment and both colonies were terminated. A new colony of C. assimilis was subsequently established from insects collected in South Africa and releases of C. assimilis from this new colony were made into coastal and inland infestations of prickly acacia. Establishment was rapid at one coastal site and the insect quickly spread to other infestations. Establishment at one inland area was also confirmed in early 2006. The establishment in coastal areas supported a CLIMEX model that indicated that the climate of coastal areas was more suitable than inland areas.

Agonosoma trilineatum (Heteroptera: Scutelleridae) a biological control agent of the weed bellyache bush, Jatropha gossypiifolia (Euphorbiaceae)

Bellyache bush, Jatropha gossypiifolia L., is a serious weed of northern Australia. Agonosoma trilineatum (F.) is an insect from tropical America released in Australia in 2003 as a biological control agent against bellyache bush. It feeds on seeds and has the potential to reduce seed production, thereby potentially reducing the rate of spread and recruitment. To test the host specificity of A. trilineatum, four biological responses to host plant species were determined: development of nymphs, oviposition preferences, adult feeding and frequency of mating. Development of nymphs to adults and adult feeding only occurred on three Jatropha spp. These species also supported mating and oogenesis but only J. gossypiifolia was accepted for oviposition. Mating did not occur in the presence of other plant species. The evidence indicates that there is little risk associated with the release of this insect species in Australia and probably other countries where this weed is a problem. The probability of this insect expanding its host range is low because multiple aspects of the biology would need to change simultaneously. A. trilineatum was released in Australia between 2003 and 2007. A Climex model indicated that coastal areas of Queensland and the Northern Territory would be climatically most suitable for this insect.

New efforts at biological control of Mikania micrantha H.B.K. (Asteraceae) in Papua New Guinea and Fiji.

New efforts at biological control of Miconia calvescens (Melastomataceae) is a serious invader in the tropical Pacific, including the Hawaiian and Tahitian Islands, and currently poses a major threat to native biodiversity in the Wet Tropics of Australia. The species is fleshy-fruited, small-seeded and shade tolerant, and thus has the potential to be dispersed widely and recruit in relatively intact rainforest habitats, displacing native species. Understanding and predicting the rate of spread is critical for the design and implementation of effective management actions. We used an individual-based model incorporating a dispersal function derived from dispersal curves for similar berry-fruited native species, and life-history parameters of fecundity and mortality to predict the spatial structure of a Miconia population after a 30 year time period. We compared the modelled population spatial structure to that of an actual infestation in the rainforests of north Queensland. Our goal was to assess how well the model predicts actual dispersion and to identify potential barriers and conduits to seed movement and seedling establishment. The model overpredicts overall population size and the spatial extent of the actual infestation, predicting individuals to occur at a maximum 1,750 m from the source compared with the maximum distance of any detected individual in the actual infestation of 1,191 m. We identify several characteristic features of managed invasive populations that make comparisons between modelled outcomes and actual infestations difficult. Our results suggest that the model’s ability to predict both spatial structure and spread of the population will be improved by incorporating a spatially explicit element, with dispersal and recruitment probabilities that reflect the relative suitability of different parts of the landscape for these processes. Mikania micrantha H.B.K. (Asteraceae) in Papua New Guinea and Fiji.

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.

Population dynamics of the European rabbit (Oryctolagus cuniculus) in north eastern Australia: simulated responses to control

Wild European rabbits are a serious problem to agriculture in Australia, with an estimated annual cost of A$ 113 million. Biological control agents (myxomatosis and rabbit haemorrhagic disease virus) have caused large and sustained declines in rabbit populations throughout Australia. A simulation model incorporates these diseases as well as warren destruction as methods of controlling rabbit populations in Queensland, north eastern Australia. These diseases reduced populations by 90-99% and the combination of these and warren destruction led to 100% control in simulations at six sites across southern Queensland. Increasing monthly pasture growth by 15% had little effect on simulated populations whereas a 15% decrease reduced populations by 0-50%. An increase in temperature of 2.5 °C would lead to a 15-60% decrease in populations. These effects suggest that climate change will lead to a decrease in the population of rabbits in Queensland and a retraction in the northern limit of their distribution in Australia.

How widespread is Parthenium hysterophorus and its biological control agent Zygogramma bicolorata in South Asia?

Parthenium hysterophorus is a weed of global significance causing severe economic, environmental, human and animal health problems in Asia, Africa, Australia and the Pacific. In South Asia, P. hysterophorus occurs in India, Pakistan, Sri Lanka, Bangladesh and Nepal. A host-specific leaf-feeding beetle Zygogramma bicolorata from Mexico was introduced into India in 1984, as a biological control agent for P. hysterophorus. In this study, a GIS-based distribution map of P. hysterophorus and its biological control agent Z. bicolorata in South Asia based on meta-analysis is presented. The map highlights the limited published information on P. hysterophorus incidence in many of the states and territories in India, as well as in neighbouring Bangladesh, Bhutan, Nepal and Pakistan. Incidence of Z. bicolorata was recorded as three distinct clusters, covering many states in India. In Pakistan, Z. bicolorata was recorded in the Punjab region bordering India. A CLIMEX model based on the current distribution of Z. bicolorata in India suggests that the geographic range of this agent in India and Pakistan can extend to other P. hysterophorus-infested areas in the region. The CLIMEX model also suggests that all of Bangladesh and Sri Lanka, and parts of Nepal are climatically suitable for Z. bicolorata.

Applying a simulation model to the management of an infestation of Miconia calvescens in the wet tropics of Australia

A simulation model that combines biological, search and economic components is applied to the eradication of a Miconia calvescens infestation at El Arish in tropical Queensland, Australia. Information on the year M. calvescens was introduced to the site, the number of plants controlled and the timing of control, is used to show that currently there could be M. calvescens plants remaining undetected at the site, including some mature plants. Modelling results indicate that the eradication programme has had a significant impact on the population of M. calvescens, as shown by simulated results for uncontrolled and controlled populations. The model was also used to investigate the effect of changing search effort on the cost of and time to eradication. Control costs were found to be negligible over all levels of search effort tested. Importantly, results suggest eradication may be achieved within several decades, if resources are increased slightly from their current levels and if there is a long-term commitment to funding the eradication programme.

Demography of feral camels in central Australia and its relevance to population control

Since their release over 100 years ago, camels have spread across central Australia and increased in number. Increasingly, they are being seen as a pest, with observed impacts from overgrazing and damage to infrastructure such as fences. Irregular aerial surveys since 1983 and an interview-based survey in 1966 suggest that camels have been increasing at close to their maximum rate. A comparison of three models of population growth fitted to these, albeit limited, data suggests that the Northern Territory population has indeed been growing at an annual exponential rate of r = 0.074, or 8% per year, with little evidence of a density-dependent brake. A stage-structured model using life history data from a central Australian camel population suggests that this rate approximates the theoretical maximum. Elasticity analysis indicates that adult survival is by far the biggest influence on rate of increase and that a 9% reduction in survival from 96% is needed to stop the population growing. In contrast, at least 70% of mature females need to be sterilised to have a similar effect. In a benign environment, a population of large mammals such as camels is expected to grow exponentially until close to carrying capacity. This will frustrate control programs, because an ever-increasing number of animals will need to be removed for zero growth the longer that culling or harvesting effort is delayed. A population projection for 2008 suggests ~10 500 animals need to be harvested across the Northern Territory. Current harvests are well short of this. The ability of commercial harvesting to control camel populations in central Australia will depend on the value of animals, access to animals and the presence of alternative species to harvest when camels are at low density.

Simple software processes and tests improve the reliability and usefulness of a model.

Models are abstractions of reality that have predetermined limits (often not consciously thought through) on what problem domains the models can be used to explore. These limits are determined by the range of observed data used to construct and validate the model. However, it is important to remember that operating the model beyond these limits, one of the reasons for building the model in the first place, potentially brings unwanted behaviour and thus reduces the usefulness of the model. Our experience with the Agricultural Production Systems Simulator (APSIM), a farming systems model, has led us to adapt techniques from the disciplines of modelling and software development to create a model development process. This process is simple, easy to follow, and brings a much higher level of stability to the development effort, which then delivers a much more useful model. A major part of the process relies on having a range of detailed model tests (unit, simulation, sensibility, validation) that exercise a model at various levels (sub-model, model and simulation). To underline the usefulness of testing, we examine several case studies where simulated output can be compared with simple relationships. For example, output is compared with crop water use efficiency relationships gleaned from the literature to check that the model reproduces the expected function. Similarly, another case study attempts to reproduce generalised hydrological relationships found in the literature. This paper then describes a simple model development process (using version control, automated testing and differencing tools), that will enhance the reliability and usefulness of a model.

Validation and development of PigBal model

This project involves validating and upgrading the PigBal model to improve the accuracy of manure production predictions from intensive piggery operations.

Validation and development of PigBal model - Stage 2

The proposed project involves validate and upgrading the PigBal model to improve the accuracy of manure and GHG production predictions.

Introduction, establishment, and potential geographic range of Carmenta sp. nr ithacae (Lepidoptera: Sesiidae), a biological control agent for Parthenium hysterophorus (Asteraceae) in Australia

Parthenium (Parthenium hysterophorus L.), a major weed causing economic, environmental, and human and animal health problems in Australia and several countries in Asia, Africa, and the Pacific, has been a target for biological control in Australia since the mid-1970s. Nine species of insects and two rust fungi have been introduced as biological control agents into Australia. These include Carmenta sp. nr ithacae, a root feeding agent from Mexico. The larvae of C. sp. nr ithacae bore through the stem-base into the root where they feed on the cortical tissue of the taproot. During 1998-2002, 2,816 larval-infested plants and 387 adults were released at 31 sites across Queensland, Australia. Evidence of field establishment was first observed in two of the release sites in central Queensland in 2004. Annual surveys at these sites and nonrelease sites during 2006-2011 showed wide variations in the incidence and abundance of C. sp. nr ithacae between years and sites. Surveys at three of the nine release sites in northern Queensland and 16 of the 22 release sites in central Queensland confirmed the field establishment of C. sp. nr ithacae in four release sites and four nonrelease sites, all in central Queensland. No field establishment was evident in the inland region or in northern Queensland. A CLIMEX model based on the native range distribution of C. sp. nr ithacae predicts that areas east of the dividing range along the coast are more suitable for field establishment than inland areas. Future efforts to redistribute this agent should be restricted to areas identified as climatically favorable by the CLIMEX model.

A review of the biology, ecology, distribution and control of Mozambique tilapia, Oreochromis mossambicus (Peters 1852) (Pisces: Cichlidae) with particular emphasis on invasive Australian populations

Oreochromis mossambicus (Peters 1852) are native to the eastward flowing rivers of central and southern Africa but from the early 1930s they have been widely distributed around the world for aquaculture and for biological control of weeds and insects. While O. mossambicus are now not commonly used as an aquaculture species, the biological traits that made them a popular culture species including tolerance to wide ranging ecological conditions, generalist dietary requirements and rapid reproduction with maternal care have also made them a 'model' invader. Self-sustaining populations now exist in almost every region to which they have been imported. In Australia, since their introduction in the 1970s, O. mossambicus have become established in catchments along the east and west coasts and have the potential to colonise other adjacent drainages. It is thought that intentional translocations are likely to be the most significant factor in their spread in Australia. The ecological and physical tolerances and preferences, reproductive behaviour, hybridization and the high degree of plasticity in the life history traits of O. mossambicus are reviewed. Impacts of O. mossambicus on natural ecosystems including competitive displacement of native species, habitat alteration, predation and as a vector in the spread of diseases are discussed. Potential methods for eradicating or controlling invasive populations of O. mossambicus including physical removal, piscicides, screens, environmental management and genetic technologies are outlined.

Individual-based modelling of the efficacy of fumigation tactics to control lesser grain borer (Rhyzopertha dominica) in stored grain

Increasing resistance to phosphine (PH 3) in insect pests, including lesser grain borer (Rhyzopertha dominica) has become a critical issue, and development of effective and sustainable strategies to manage resistance is crucial. In practice, the same grain store may be fumigated multiple times, but usually for the same exposure period and concentration. Simulating a single fumigation allows us to look more closely at the effects of this standard treatment.We used an individual-based, two-locus model to investigate three key questions about the use of phosphine fumigant in relation to the development of PH 3 resistance. First, which is more effective for insect control; long exposure time with a low concentration or short exposure period with a high concentration? Our results showed that extending exposure duration is a much more efficient control tactic than increasing the phosphine concentration. Second, how long should the fumigation period be extended to deal with higher frequencies of resistant insects in the grain? Our results indicated that if the original frequency of resistant insects is increased n times, then the fumigation needs to be extended, at most, n days to achieve the same level of insect control. The third question is how does the presence of varying numbers of insects inside grain storages impact the effectiveness of phosphine fumigation? We found that, for a given fumigation, as the initial population number was increased, the final survival of resistant insects increased proportionally. To control initial populations of insects that were n times larger, it was necessary to increase the fumigation time by about n days. Our results indicate that, in a 2-gene mediated resistance where dilution of resistance gene frequencies through immigration of susceptibles has greater effect, extending fumigation times to reduce survival of homozygous resistant insects will have a significant impact on delaying the development of resistance. © 2012 Elsevier Ltd.