Assessing the sustainability of wheat-based cropping systems using simulation modelling: sustainability

Concepts of agricultural sustainability and possible roles of simulation modelling for characterising sustainability were explored by conducting, and reflecting on, a sustainability assessment of rain-fed wheat-based systems in the Middle East and North Africa region. We designed a goal-oriented, model-based framework using the cropping systems model Agricultural Production Systems sIMulator (APSIM). For the assessment, valid (rather than true or false) sustainability goals and indicators were identified for the target system. System-specific vagueness was depicted in sustainability polygons-a system property derived from highly quantitative data-and denoted using descriptive quantifiers. Diagnostic evaluations of alternative tillage practices demonstrated the utility of the framework to quantify key bio-physical and chemical constraints to sustainability. Here, we argue that sustainability is a vague, emergent system property of often wicked complexity that arises out of more fundamental elements and processes. A 'wicked concept of sustainability' acknowledges the breadth of the human experience of sustainability, which cannot be internalised in a model. To achieve socially desirable sustainability goals, our model-based approach can inform reflective evaluation processes that connect with the needs and values of agricultural decision-makers. Hence, it can help to frame meaningful discussions, from which actions might emerge.

Resting pastures to improve land condition in northern Australia: guidelines based on the literature and simulation modelling

Pasture rest is a possible strategy for improving land condition in the extensive grazing lands of northern Australia. If pastures currently in poor condition could be improved, then overall animal productivity and the sustainability of grazing could be increased. The scientific literature is examined to assess the strength of the experimental information to support and guide the use of pasture rest, and simulation modelling is undertaken to extend this information to a broader range of resting practices, growing conditions and initial pasture condition. From this, guidelines are developed that can be applied in the management of northern Australia’s grazing lands and also serve as hypotheses for further field experiments. The literature on pasture rest is diverse but there is a paucity of data from much of northern Australia as most experiments have been conducted in southern and central parts of Queensland. Despite this, the limited experimental information and the results from modelling were used to formulate the following guidelines. Rest during the growing season gives the most rapid improvement in the proportion of perennial grasses in pastures; rest during the dormant winter period is ineffective in increasing perennial grasses in a pasture but may have other benefits. Appropriate stocking rates are essential to gain the greatest benefit from rest: if stocking rates are too high, then pasture rest will not lead to improvement; if stocking rates are low, pastures will tend to improve without rest. The lower the initial percentage of perennial grasses, the more frequent the rests should be to give a major improvement within a reasonable management timeframe. Conditions during the growing season also have an impact on responses with the greatest improvement likely to be in years of good growing conditions. The duration and frequency of rest periods can be combined into a single value expressed as the proportion of time during which resting occurs; when this is done the modelling suggests the greater the proportion of time that a pasture is rested, the greater is the improvement but this needs to be tested experimentally. These guidelines should assist land managers to use pasture resting but the challenge remains to integrate pasture rest with other pasture and animal management practices at the whole-property scale.

Boosting the productivity and profitability of northern Australian beef enterprises: Exploring innovation options using simulation modelling and systems analysis

The financial health of beef cattle enterprises in northern Australia has declined markedly over the last decade due to an escalation in production and marketing costs and a real decline in beef prices. Historically, gains in animal productivity have offset the effect of declining terms of trade on farm incomes. This raises the question of whether future productivity improvements can remain a key path for lifting enterprise profitability sufficient to ensure that the industry remains economically viable over the longer term. The key objective of this study was to assess the production and financial implications for north Australian beef enterprises of a range of technology interventions (development scenarios), including genetic gain in cattle, nutrient supplementation, and alteration of the feed base through introduced pastures and forage crops, across a variety of natural environments. To achieve this objective a beef systems model was developed that is capable of simulating livestock production at the enterprise level, including reproduction, growth and mortality, based on energy and protein supply from natural C4 pastures that are subject to high inter-annual climate variability. Comparisons between simulation outputs and enterprise performance data in three case study regions suggested that the simulation model (the Northern Australia Beef Systems Analyser) can adequately represent the performance beef cattle enterprises in northern Australia. Testing of a range of development scenarios suggested that the application of individual technologies can substantially lift productivity and profitability, especially where the entire feedbase was altered through legume augmentation. The simultaneous implementation of multiple technologies that provide benefits to different aspects of animal productivity resulted in the greatest increases in cattle productivity and enterprise profitability, with projected weaning rates increasing by 25%, liveweight gain by 40% and net profit by 150% above current baseline levels, although gains of this magnitude might not necessarily be realised in practice. While there were slight increases in total methane output from these development scenarios, the methane emissions per kg of beef produced were reduced by 20% in scenarios with higher productivity gain. Combinations of technologies or innovative practices applied in a systematic and integrated fashion thus offer scope for providing the productivity and profitability gains necessary to maintain viable beef enterprises in northern Australia into the future.

Modelling the manager: Representing rule-based management in farming systems simulation models

We trace the evolution of the representation of management in cropping and grazing systems models, from fixed annual schedules of identical actions in single paddocks toward flexible scripts of rules. Attempts to define higher-level organizing concepts in management policies, and to analyse them to identify optimal plans, have focussed on questions relating to grazing management owing to its inherent complexity. “Rule templates” assist the re-use of complex management scripts by bundling commonly-used collections of rules with an interface through which key parameters can be input by a simulation builder. Standard issues relating to parameter estimation and uncertainty apply to management sub-models and need to be addressed. Techniques for embodying farmers' expectations and plans for the future within modelling analyses need to be further developed, especially better linking planning- and rule-based approaches to farm management and analysing the ways that managers can learn.

The effect of row configuration on yield reliability in grain sorghum: II. Modelling the effects of row configuration

In recent years many sorghum producers in the more marginal (<600 mm annual rainfall) cropping areas of Qld and northern NSW have utilised skip row configurations in an attempt to improve yield reliability and reduce sorghum production risk. But will this work in the long run? What are the trade-offs between productivity and risk of crop failure? This paper describes a modelling and simulation approach to study the long-term effects of skip row configurations. Detailed measurements of light interception and water extraction from sorghum crops grown in solid, single and double skip row configurations were collected from three on-farm participatory research trials established in southern Qld and northern NSW. These measurements resulted in changes to the model that accounted for the elliptical water uptake pattern below the crop row and reduced total light interception associated with the leaf area reduction of the skip configuration. Following validation of the model, long-term simulation runs using historical weather data were used to determine the value of skip row sorghum production as a means of maintaining yield reliability in the dryland cropping regions of southern Qld and northern NSW.

Simulation of the supply capability of sheep on the Mitchell grass downs of north west Queensland

Sheep in western Queensland have been predominantly reared for wool. When wool prices became depressed interest in the sheep meat industry, increased. For north west Queensland producers, opportunities may exist to participate in live sheep and meat export to Asia. A simulation model was developed to determine whether this sheep producing area has the capability to provide sufficient numbers of sheep under variable climatic conditions while sustaining the land resources. Maximum capacity for sustainability of resources (as described by stock numbers) was derived from an in-depth study of the agricultural and pastoral potential of Queensland. Decades of sheep production and climatic data spanning differing seasonal conditions were collated for analysis. A ruminant biology model adapted from Grazplan was used to simulate pregnancy rate. Empirical equations predict mortalities, marking rates, and weight characteristics of sheep of various ages from simple climatic measures, stocking rate and reproductive status. The initial age structure of flocks was determined by running the model for several years with historical climatic conditions. Drought management strategies such as selling a proportion of wethers progressively down to two-tooth and oldest ewes were incorporated. Management decisions such as time of joining, age at which ewes were cast-for-age, wether turn-off age and turning-off rate of lambs vary with geographical area and can be specified at run time. The model is run for sequences of climatic conditions generated stochastically from distributions based on historical climatic data correlated in some instances. The model highlights the difficulties of sustaining a consistent supply of sheep under variable climatic conditions.

Modelling Mortality and Reproduction Rates for Management of Sheep Flocks in Northern Australia

By quantifying the effects of climatic variability in the sheep grazing lands of north western and western Queensland, the key biological rates of mortality and reproduction can be predicted for sheep. These rates are essential components of a decision support package which can prove a useful management tool for producers, especially if they can easily obtain the necessary predictors. When the sub-models of the GRAZPLAN ruminant biology process model were re-parameterised from Queensland data along with an empirical equation predicting the probability of ewes mating added, the process model predicted the probability of pregnancy well (86% variation explained). Predicting mortality from GRAZPLAN was less successful but an empirical equation based on relative condition of the animal (a measure based on liveweight), pregnancy status and age explained 78% of the variation in mortalities. A crucial predictor in these models was liveweight which is not often recorded on producer properties. Empirical models based on climatic and pasture conditions estimated from the pasture production model GRASP, predicted marking and mortality rates for Mitchell grass (Astrebla sp.) pastures (81% and 63% of the variation explained). These prediction equations were tested against independent data from producer properties and the model successfully validated for Mitchell grass communities.

Modelling Leaf Production and Crop Development in Maize (Zea mays L.) After Tassel Initiation Under Diverse Conditions of Temperature and Photoperiod

Prediction of the initiation, appearance and emergence of leaves is critically important to the success of simulation models of crop canopy development and some aspects of crop ontogeny. Data on leaf number and crop ontogeny were collected on five cultivars of maize differing widely in maturity and genetic background grown under natural and extended photoperiods, and planted on seven sowing dates from October 1993 to March 1994 at Gatton, South-east Queensland. The same temperature coefficients were established for crop ontogeny before silking, and the rates of leaf initiation, leaf tip appearance and full leaf expansion, the base, optimum and maximum temperatures for each being 8, 34 and 40 degrees C. After silking, the base temperature for ontogeny was 0 degrees C, but the optimum and maximum temperatures remained unchanged. The rates of leaf initiation, appearance of leaf tips and full leaf expansion varied in a relatively narrow range across sowing times and photoperiod treatments, with average values of 0.040 leaves (degrees Cd)-1, 0.021 leaves (degrees Cd)-1, and 0.019 leaves (degrees Cd)-1, respectively. The relationships developed in this study provided satisfactory predictions of leaf number and crop ontogeny (tassel initiation to silking, emergence to silking and silking to physiological maturity) when assessed using independent data from Gatton (South eastern Queensland), Katherine and Douglas Daly (Northern Territory), Walkamin (North Queensland) and Kununurra (Western Australia).

Short-term responses of leaf growth rate to water deficit scale up to whole-plant and crop levels: An integrated modelling approach in maize.

Physiological and genetic studies of leaf growth often focus on short-term responses, leaving a gap to whole-plant models that predict biomass accumulation, transpiration and yield at crop scale. To bridge this gap, we developed a model that combines an existing model of leaf 6 expansion in response to short-term environmental variations with a model coordinating the development of all leaves of a plant. The latter was based on: (1) rates of leaf initiation, appearance and end of elongation measured in field experiments; and (2) the hypothesis of an independence of the growth between leaves. The resulting whole-plant leaf model was integrated into the generic crop model APSIM which provided dynamic feedback of environmental conditions to the leaf model and allowed simulation of crop growth at canopy level. The model was tested in 12 field situations with contrasting temperature, evaporative demand and soil water status. In observed and simulated data, high evaporative demand reduced leaf area at the whole-plant level, and short water deficits affected only leaves developing during the stress, either visible or still hidden in the whorl. The model adequately simulated whole-plant profiles of leaf area with a single set of parameters that applied to the same hybrid in all experiments. It was also suitable to predict biomass accumulation and yield of a similar hybrid grown in different conditions. This model extends to field conditions existing knowledge of the environmental controls of leaf elongation, and can be used to simulate how their genetic controls flow through to yield.

Modelling the potential spread of Solenopsis invicta Buren (Hymenoptera: Formicidae) (red imported fire ant) in Australia

Quantifying the potential spread and density of an invading organism enables decision-makers to determine the most appropriate response to incursions. We present two linked models that estimate the spread of Solenopsis invicta Buren (red imported fire ant) in Australia based on limited data gathered after its discovery in Brisbane in 2001. A stochastic cellular automaton determines spread within a location (100 km by 100 km) and this is coupled with a model that simulates human-mediated movement of S. invicta to new locations. In the absence of any control measures, the models predict that S. invicta could cover 763 000–4 066 000 km2 by the year 2035 and be found at 200 separate locations around Australia by 2017–2027, depending on the rate of spread. These estimated rates of expansion (assuming no control efforts were in place) are higher than those experienced in the USA in the 1940s during the early invasion phases in that country. Active control efforts and quarantine controls in the USA (including a concerted eradication attempt in the 1960s) may have slowed spread. Further, milder winters, the presence of the polygynous social form, increased trade and human mobility in Australia in 2000s compared with the USA in 1940s could contribute to faster range expansion.

Modelling climatic risks of aflatoxin contamination in maize

Aflatoxins are highly carcinogenic mycotoxins produced by two fungi, Aspergillus flavus and A. parasiticus, under specific moisture and temperature conditions before harvest and/or during storage of a wide range of crops including maize. Modelling of interactions between host plant and environment during the season can enable quantification of preharvest aflatoxin risk and its potential management. A model was developed to quantify climatic risks of aflatoxin contamination in maize using principles previously used for peanuts. The model outputs an aflatoxin risk index in response to seasonal temperature and soil moisture during the maize grain filling period using the APSIM's maize module. The model performed well in simulating climatic risk of aflatoxin contamination in maize as indicated by a significant R2 (P ≤ 0.01) between aflatoxin risk index and the measured aflatoxin B1 in crop samples, which was 0.69 for a range of rainfed Australian locations and 0.62 when irrigated locations were also included in the analysis. The model was further applied to determine probabilities of exceeding a given aflatoxin risk in four non-irrigated maize growing locations of Queensland using 106 years of historical climatic data. Locations with both dry and hot climates had a much higher probability of higher aflatoxin risk compared with locations having either dry or hot conditions alone. Scenario analysis suggested that under non-irrigated conditions the risk of aflatoxin contamination could be minimised by adjusting sowing time or selecting an appropriate hybrid to better match the grain filling period to coincide with lower temperature and water stress conditions.

Architectural modelling of maize under water stress

Two field experiments using maize (Pioneer 31H50) and three watering regimes [(i) irrigated for the whole crop cycle, until anthesis, (ii) not at all (experiment 1) and (iii) fully irrigated and rain grown for the whole crop cycle (experiment 2)] were conducted at Gatton, Australia, during the 2003-04 season. Data on crop ontogeny, leaf, sheath and internode lengths and leaf width, and senescence were collected at 1- to 3-day intervals. A glasshouse experiment during 2003 quantified the responses of leaf shape and leaf presentation to various levels of water stress. Data from experiment 1 were used to modify and parameterise an architectural model of maize (ADEL-Maize) to incorporate the impact of water stress on maize canopy characteristics. The modified model produced accurate fitted values for experiment 1 for final leaf area and plant height, but values during development for leaf area were lower than observed data. Crop duration was reasonably well fitted and differences between the fully irrigated and rain-grown crops were accurately predicted. Final representations of maize crop canopies were realistic. Possible explanations for low values of leaf area are provided. The model requires further development using data from the glasshouse study and before being validated using data from experiment 2 and other independent data. It will then be used to extend functionality in architectural models of maize. With further research and development, the model should be particularly useful in examining the response of maize production to water stress including improved prediction of total biomass and grain yield. This will facilitate improved simulation of plant growth and development processes allowing investigation of genotype by environment interactions under conditions of suboptimal water supply.

Graze to grain - measuring and modelling the effects of grazed pasture leys on soil nitrogen and sorghum yield on a Vertosol soil in the Australian subtropics

Highly productive sown pasture systems can result in high growth rates of beef cattle and lead to increases in soil nitrogen and the production of subsequent crops. The nitrogen dynamics and growth of grain sorghum following grazed annual legume leys or a grass pasture were investigated in a no-till system in the South Burnett district of Queensland. Two years of the tropical legumes Macrotyloma daltonii and Vigna trilobata (both self regenerating annual legumes) and Lablab purpureus (a resown annual legume) resulted in soil nitrate N (0-0.9 m depth), at sorghum sowing, ranging from 35 to 86 kg/ha compared with 4 kg/ha after pure grass pastures. Average grain sorghum production in the 4 cropping seasons following the grazed legume leys ranged from 2651 to 4012 kg/ha. Following the grass pasture, grain sorghum production in the first and second year was < 1900 kg/ha and by the third year grain yield was comparable to the legume systems. Simulation studies utilising the farming systems model APSIM indicated that the soil N and water dynamics following 2-year ley phases could be closely represented over 4 years and the prediction of sorghum growth during this time was reasonable. In simulated unfertilised sorghum crops grown from 1954 to 2004, grain yield did not exceed 1500 kg/ha in 50% of seasons following a grass pasture, while following 2-year legume leys, grain exceeded 3000 kg/ha in 80% of seasons. It was concluded that mixed farming systems that utilise short term legume-based pastures for beef production in rotation with crop production enterprises can be highly productive.

Modelling environmental effects on phenology and canopy development of diverse sorghum genotypes

The ability to predict phenology and canopy development is critical in crop models used for simulating likely consequences of alternative crop management and cultivar choice strategies. Here we quantify and contrast the temperature and photoperiod responses for phenology and canopy development of a diverse range of elite Indian and Australian sorghum genotypes (hybrid and landrace). Detailed field experiments were undertaken in Australia and India using a range of genotypes, sowing dates, and photoperiod extension treatments. Measurements of timing of developmental stages and leaf appearance were taken. The generality of photo-thermal approaches to modelling phenological and canopy development was tested. Environmental and genotypic effects on rate of progression from emergence to floral initiation (E-FI) were explained well using a multiplicative model, which combined the intrinsic development rate (Ropt), with responses to temperature and photoperiod. Differences in Ropt and extent of the photoperiod response explained most genotypic effects. Average leaf initiation rate (LIR), leaf appearance rate and duration of the phase from anthesis to physiological maturity differed among genotypes. The association of total leaf number (TLN) with photoperiod found for all genotypes could not be fully explained by effects on development and LIRs. While a putative effect of photoperiod on LIR would explain the observations, other possible confounding factors, such as air-soil temperature differential and the nature of model structure were considered and discussed. This study found a generally robust predictive capacity of photo-thermal development models across diverse ranges of both genotypes and environments. Hence, they remain the most appropriate models for simulation analysis of genotype-by-management scenarios in environments varying broadly in temperature and photoperiod.

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.