Comparative productivity of irrigated short-term ryegrass (Lolium multiflorum) pasture receiving nitrogen, grown alone or in a mixture with white (Trifolium repens) and Persian (T. resupinatum) clovers

Dairy farms in subtropical Australia use irrigated, annually sown short-term ryegrass (Lolium multiflorum) or mixtures of short-term ryegrass and white (Trifolium repens) and Persian (shaftal) (T. resupinatum) clover during the winter-spring period in all-year-round milk production systems. A series of small plot cutting experiments was conducted in 3 dairying regions (tropical upland, north Queensland, and subtropical southeast Queensland and northern New South Wales) to determine the most effective rate and frequency of application of nitrogen (N) fertiliser. The experiments were not grazed, nor was harvested material returned to the plots, after sampling. Rates up to 100 kg N/ha.month (as urea or calcium ammonium nitrate) and up to 200 kg N/ha every 2 months (as urea) were applied to pure stands of ryegrass in 1991. In 1993 and 1994, urea, at rates up to 150 kg N/ha.month and to 200 kg N/ha every 2 months, was applied to pure stands of ryegrass; urea, at rates up to 50 kg N/ha.month, was also applied to ryegrass-clover mixtures. The results indicate that applications of 50-85 kg N/ha.month can be recommended for short-term ryegrass pastures throughout the subtropics and tropical uplands of eastern Australia, irrespective of soil type. At this rate, dry matter yields will reach about 90% of their potential, forage nitrogen concentration will be increased, there is minimal risk to stock from nitrate poisoning and there will be no substantial increase in soil N. The rate of N for ryegrass-clover pastures is slightly higher than for pure ryegrass but, at these rates, the clover component will be suppressed. However, increased ryegrass yields and higher forage nitrogen concentrations will compensate for the reduced clover component. At application rates up to 100 kg N/ha.month, build-up of NO3--N and NH4+-N in soil was generally restricted to the surface layers (0-20 cm) of the soil, but there was a substantial increase throughout the soil profile at 150 kg N/ha.month. The build-up of NO3--N and NH4+-N was greater and was found at lower rates on the lighter soil compared with heavy clays. Generally, most of the soil N was in the NO3--N form and most was in the top 20 cm.

Delineating environmental envelopes to improve mapping of species distributions, via a hurdle model with CART &/or MaxEnT

Species distribution modelling (SDM) typically analyses species’ presence together with some form of absence information. Ideally absences comprise observations or are inferred from comprehensive sampling. When such information is not available, then pseudo-absences are often generated from the background locations within the study region of interest containing the presences, or else absence is implied through the comparison of presences to the whole study region, e.g. as is the case in Maximum Entropy (MaxEnt) or Poisson point process modelling. However, the choice of which absence information to include can be both challenging and highly influential on SDM predictions (e.g. Oksanen and Minchin, 2002). In practice, the use of pseudo- or implied absences often leads to an imbalance where absences far outnumber presences. This leaves analysis highly susceptible to ‘naughty-noughts’: absences that occur beyond the envelope of the species, which can exert strong influence on the model and its predictions (Austin and Meyers, 1996). Also known as ‘excess zeros’, naughty noughts can be estimated via an overall proportion in simple hurdle or mixture models (Martin et al., 2005). However, absences, especially those that occur beyond the species envelope, can often be more diverse than presences. Here we consider an extension to excess zero models. The two-staged approach first exploits the compartmentalisation provided by classification trees (CTs) (as in O’Leary, 2008) to identify multiple sources of naughty noughts and simultaneously delineate several species envelopes. Then SDMs can be fit separately within each envelope, and for this stage, we examine both CTs (as in Falk et al., 2014) and the popular MaxEnt (Elith et al., 2006). We introduce a wider range of model performance measures to improve treatment of naughty noughts in SDM. We retain an overall measure of model performance, the area under the curve (AUC) of the Receiver-Operating Curve (ROC), but focus on its constituent measures of false negative rate (FNR) and false positive rate (FPR), and how these relate to the threshold in the predicted probability of presence that delimits predicted presence from absence. We also propose error rates more relevant to users of predictions: false omission rate (FOR), the chance that a predicted absence corresponds to (and hence wastes) an observed presence, and the false discovery rate (FDR), reflecting those predicted (or potential) presences that correspond to absence. A high FDR may be desirable since it could help target future search efforts, whereas zero or low FOR is desirable since it indicates none of the (often valuable) presences have been ignored in the SDM. For illustration, we chose Bradypus variegatus, a species that has previously been published as an exemplar species for MaxEnt, proposed by Phillips et al. (2006). We used CTs to increasingly refine the species envelope, starting with the whole study region (E0), eliminating more and more potential naughty noughts (E1–E3). When combined with an SDM fit within the species envelope, the best CT SDM had similar AUC and FPR to the best MaxEnt SDM, but otherwise performed better. The FNR and FOR were greatly reduced, suggesting that CTs handle absences better. Interestingly, MaxEnt predictions showed low discriminatory performance, with the most common predicted probability of presence being in the same range (0.00-0.20) for both true absences and presences. In summary, this example shows that SDMs can be improved by introducing an initial hurdle to identify naughty noughts and partition the envelope before applying SDMs. This improvement was barely detectable via AUC and FPR yet visible in FOR, FNR, and the comparison of predicted probability of presence distribution for pres/absence.

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.

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.

Mathematical modelling of the impaction and spreading of spray droplets on leaves

This thesis concerns the development of mathematical models to describe the interactions that occur between spray droplets and leaves. Models are presented that not only provide a contribution to mathematical knowledge in the field of fluid dynamics, but are also of utility within the agrichemical industry. The thesis is presented in two parts. First, thin film models are implemented with efficient numerical schemes in order to simulate droplets on virtual leaf surfaces. Then the interception event is considered, whereby energy balance techniques are employed to instantaneously predict whether an impacting droplet will bounce, splash, or adhere to a leaf.

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.

Special Issue on Spatial Moment Techniques for Modelling Biological Processes

Over the last two decades, there has been an increasing awareness of, and interest in, the use of spatial moment techniques to provide insight into a range of biological and ecological processes. Models that incorporate spatial moments can be viewed as extensions of mean-field models. These mean-field models often consist of systems of classical ordinary differential equations and partial differential equations, whose derivation, at some point, hinges on the simplifying assumption that individuals in the underlying stochastic process encounter each other at a rate that is proportional to the average abundance of individuals. This assumption has several implications, the most striking of which is that mean-field models essentially neglect any impact of the spatial structure of individuals in the system. Moment dynamics models extend traditional mean-field descriptions by accounting for the dynamics of pairs, triples and higher n-tuples of individuals. This means that moment dynamics models can, to some extent, account for how the spatial structure affects the dynamics of the system in question.

Modelling the movement of interacting cell populations: A moment dynamics approach

Mathematical models describing the movement of multiple interacting subpopulations are relevant to many biological and ecological processes. Standard mean-field partial differential equation descriptions of these processes suffer from the limitation that they implicitly neglect to incorporate the impact of spatial correlations and clustering. To overcome this, we derive a moment dynamics description of a discrete stochastic process which describes the spreading of distinct interacting subpopulations. In particular, we motivate our model by mimicking the geometry of two typical cell biology experiments. Comparing the performance of the moment dynamics model with a traditional mean-field model confirms that the moment dynamics approach always outperforms the traditional mean-field approach. To provide more general insight we summarise the performance of the moment dynamics model and the traditional mean-field model over a wide range of parameter regimes. These results help distinguish between those situations where spatial correlation effects are sufficiently strong, such that a moment dynamics model is required, from other situations where spatial correlation effects are sufficiently weak, such that a traditional mean-field model is adequate.

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.

The performance of irrigated mixtures of tall fescue, ryegrass and white clover in subtropical Australia. 1. The effects of sowing mixture combinations, nitrogen and oversowing on establishment, productivity, botanical composition and persistence.

In the subtropics of Australia, the ryegrass component of irrigated perennial ryegrass (Lolium perenne) - white clover (Trifolium repens) pastures declines by approximately 40% in the summer following establishment, being replaced by summer-active C4 grasses. Tall fescue (Festuca arundinacea) is more persistent than perennial ryegrass and might resist this invasion, although tall fescue does not compete vigorously as a seedling. This series of experiments investigated the influence of ryegrass and tall fescue genotype, sowing time and sowing mixture as a means of improving tall fescue establishment and the productivity and persistence of tall fescue, ryegrass and white clover-based mixtures in a subtropical environment. Tall fescue frequency at the end of the establishment year decreased as the number of companion species sown in the mixture increased. Neither sowing mixture combinations nor sowing rates influenced overall pasture yield (of around 14 t/ha) in the establishment year but had a significant effect on botanical composition and component yields. Perennial ryegrass was less competitive than short-rotation ryegrass, increasing first-year yields of tall fescue by 40% in one experiment and by 10% in another but total yield was unaffected. The higher establishment-year yield (3.5 t/ha) allowed Dovey tall fescue to compete more successfully with the remaining pasture components than Vulcan (1.4 t/ha). Sowing 2 ryegrass cultivars in the mixture reduced tall fescue yields by 30% compared with a single ryegrass (1.6 t/ha), although tall fescue alone achieved higher yields (7.1 t/ha). Component sowing rate had little influence on composition or yield. Oversowing the ryegrass component into a 6-week-old sward of tall fescue and white clover improved tall fescue, white clover and overall yields in the establishment year by 83, 17 and 11%, respectively, but reduced ryegrass yields by 40%. The inclusion of red (T. pratense) and Persian (T. resupinatum) clovers and chicory (Cichorium intybus) increased first-year yields by 25% but suppressed perennial grass and clover components. Yields were generally maintained at around 12 t/ha/yr in the second and third years, with tall fescue becoming dominant in all 3 experiments. The lower tall fescue seeding rate used in the first experiment resulted in tall fescue dominance in the second year following establishment, whereas in Experiments 2 and 3 dominance occurred by the end of the first year. Invasion by the C4 grasses was relatively minor (<10%) even in the third year. As ryegrass plants died, tall fescue and, to a lesser extent, white clover increased as a proportion of the total sward. Treatment effects continued into the second, but rarely the third, year and mostly affected the yield of one of the components rather than total cumulative yield. Once tall fescue became dominant, it was difficult to re-introduce other pasture components, even following removal of foliage and moderate renovation. Severe renovation (reducing the tall fescue population by at least 30%) seems a possible option for redressing this situation.

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.