APSIM – Evolution towards a new generation of agricultural systems simulation

Agricultural systems models worldwide are increasingly being used to explore options and solutions for the food security, climate change adaptation and mitigation and carbon trading problem domains. APSIM (Agricultural Production Systems sIMulator) is one such model that continues to be applied and adapted to this challenging research agenda. From its inception twenty years ago, APSIM has evolved into a framework containing many of the key models required to explore changes in agricultural landscapes with capability ranging from simulation of gene expression through to multi-field farms and beyond. Keating et al. (2003) described many of the fundamental attributes of APSIM in detail. Much has changed in the last decade, and the APSIM community has been exploring novel scientific domains and utilising software developments in social media, web and mobile applications to provide simulation tools adapted to new demands. This paper updates the earlier work by Keating et al. (2003) and chronicles the changing external challenges and opportunities being placed on APSIM during the last decade. It also explores and discusses how APSIM has been evolving to a “next generation” framework with improved features and capabilities that allow its use in many diverse topics.

Quantifying the response of cotton production in eastern Australia to climate change

Abstract The paper evaluates the effect of future climate change (as per the CSIRO Mk3.5 A1FI future climate projection) on cotton yield in Southern Queensland and Northern NSW, eastern Australia by using of the biophysical simulation model APSIM (Agricultural Production Systems sIMulator). The simulations of cotton production show that changes in the influential meteorological parameters caused by climate change would lead to decreased future cotton yields without the effect of CO2 fertilisation. By 2050 the yields would decrease by 17 %. Including the effects of CO2 fertilisation ameliorates the effect of decreased water availability and yields increase by 5.9 % by 2030, but then decrease by 3.6 % in 2050. Importantly, it was necessary to increase irrigation amounts by almost 50 % to maintain adequate soil moisture levels. The effect of CO2 was found to have an important positive impact of the yield in spite of deleterious climate change. This implies that the physiological response of plants to climate change needs to be thoroughly understood to avoid making erroneous projections of yield and potentially stifling investment or increasing risk.

Integrating pest population models with biophysical crop models to better represent the farming system

Farming systems frameworks such as the Agricultural Production Systems simulator (APSIM) represent fluxes through the soil, plant and atmosphere of the system well, but do not generally consider the biotic constraints that function within the system. We designed a method that allowed population models built in DYMEX to interact with APSIM. The simulator engine component of the DYMEX population-modelling platform was wrapped within an APSIM module allowing it to get and set variable values in other APSIM models running in the simulation. A rust model developed in DYMEX is used to demonstrate how the developing rust population reduces the crop's green leaf area. The success of the linking process is seen in the interaction of the two models and how changes in rust population on the crop's leaves feedback to the APSIM crop modifying the growth and development of the crop's leaf area. This linking of population models to simulate pest populations and biophysical models to simulate crop growth and development increases the complexity of the simulation, but provides a tool to investigate biotic constraints within farming systems and further moves APSIM towards being an agro-ecological framework.

The shifting influence of drought and heat stress for crops in Northeast Australia

Characterization of drought environment types (ETs) has proven useful for breeding crops for drought-prone regions. Here we consider how changes in climate and atmospheric carbon dioxide (CO2) concentrations will affect drought ET frequencies in sorghum and wheat systems of Northeast Australia. We also modify APSIM (the Agricultural Production Systems Simulator) to incorporate extreme heat effects on grain number and weight, and then evaluate changes in the occurrence of heat-induced yield losses of more than 10, as well as the co-occurrence of drought and heat. More than six million simulations spanning representative locations, soil types, management systems, and 33 climate projections led to three key findings. First, the projected frequency of drought decreased slightly for most climate projections for both sorghum and wheat, but for different reasons. In sorghum, warming exacerbated drought stresses by raising the atmospheric vapor pressure deficit and reducing transpiration efficiency (TE), but an increase in TE due to elevated CO2 more than offset these effects. In wheat, warming reduced drought stress during spring by hastening development through winter and reducing exposure to terminal drought. Elevated CO2 increased TE but also raised radiation use efficiency and overall growth rates and water use, thereby offsetting much of the drought reduction from warming. Second, adding explicit effects of heat on grain number and grain size often switched projected yield impacts from positive to negative. Finally, although average yield losses associated with drought will remain generally higher than for heat stress for the next half century, the relative importance of heat is steadily growing. This trend, as well as the likely high degree of genetic variability in heat tolerance, suggests that more emphasis on heat tolerance is warranted in breeding programs. At the same time, work on drought tolerance should continue with an emphasis on drought that co-occurs with extreme heat. This article is protected by copyright. All rights reserved.

Physiological basis of yield variation in response to row spacing and plant density of mungbean grown in subtropical environments

In this study, we investigated the extent and physiological bases of yield variation due to row spacing and plant density configuration in the mungbean Vigna radiata (L.) Wilczek variety “Crystal” grown in different subtropical environments. Field trials were conducted in six production environments; one rain-fed and one irrigated trial each at Biloela and Emerald, and one rain-fed trial each at Hermitage and Kingaroy sites in Queensland, Australia. In each trial, six combinations of spatial arrangement of plants, achieved through two inter-row spacings of 1 m or 0.9 m (wide row), 0.5 m or 0.3 m (narrow row), with three plant densities, 20, 30 and 40 plants/m2, were compared. The narrow row spacing resulted in 22% higher shoot dry matter and 14% more yield compared to the wide rows. The yield advantage of narrow rows ranged from 10% to 36% in the two irrigated and three rain-fed trials. However, yield loss of up to 10% was also recorded from narrow rows at Emerald where the crop suffered severe drought. Neither the effects of plant density, nor the interaction between plant density and row spacing, however, were significant in any trial. The yield advantage of narrow rows was related to 22% more intercepted radiation. In addition, simulations by the Agricultural Production Systems Simulator model, using site-specific agronomy, soil and weather information, suggested that narrow rows had proportionately greater use of soil water through transpiration, compared to evaporation resulting in higher yield per mm of soil water. The long-term simulation of yield probabilities over 123 years for the two row configurations showed that the mungbean crop planted in narrow rows could produce up to 30% higher grain yield compared to wide rows in 95% of the seasons.

The shifting influence of drought and heat stress for crops in Northeast Australia

Characterization of drought environment types (ETs) has proven useful for breeding crops for drought-prone regions. Here we consider how changes in climate and atmospheric carbon dioxide (CO2) concentrations will affect drought ET frequencies in sorghum and wheat systems of Northeast Australia. We also modify APSIM (the Agricultural Production Systems Simulator) to incorporate extreme heat effects on grain number and weight, and then evaluate changes in the occurrence of heat-induced yield losses of more than 10%, as well as the co-occurrence of drought and heat. More than six million simulations spanning representative locations, soil types, management systems, and 33 climate projections led to three key findings. First, the projected frequency of drought decreased slightly for most climate projections for both sorghum and wheat, but for different reasons. In sorghum, warming exacerbated drought stresses by raising the atmospheric vapor pressure deficit and reducing transpiration efficiency (TE), but an increase in TE due to elevated CO2 more than offset these effects. In wheat, warming reduced drought stress during spring by hastening development through winter and reducing exposure to terminal drought. Elevated CO2 increased TE but also raised radiation use efficiency and overall growth rates and water use, thereby offsetting much of the drought reduction from warming. Second, adding explicit effects of heat on grain number and grain size often switched projected yield impacts from positive to negative. Finally, although average yield losses associated with drought will remain generally higher than for heat stress for the next half century, the relative importance of heat is steadily growing. This trend, as well as the likely high degree of genetic variability in heat tolerance, suggests that more emphasis on heat tolerance is warranted in breeding programs. At the same time, work on drought tolerance should continue with an emphasis on drought that co-occurs with extreme heat. This article is protected by copyright. All rights reserved.

An improved simulation model to predict pre-harvest aflatoxin risk in maize

Aflatoxin is a potent carcinogen produced by Aspergillus flavus, which frequently contaminates maize (Zea mays L.) in the field between 40° north and 40° south latitudes. A mechanistic model to predict risk of pre-harvest contamination could assist in management of this very harmful mycotoxin. In this study we describe an aflatoxin risk prediction model which is integrated with the Agricultural Production Systems Simulator (APSIM) modelling framework. The model computes a temperature function for A. flavus growth and aflatoxin production using a set of three cardinal temperatures determined in the laboratory using culture medium and intact grains. These cardinal temperatures were 11.5 °C as base, 32.5 °C as optimum and 42.5 °C as maximum. The model used a low (≤0.2) crop water supply to demand ratio—an index of drought during the grain filling stage to simulate maize crop's susceptibility to A. flavus growth and aflatoxin production. When this low threshold of the index was reached the model converted the temperature function into an aflatoxin risk index (ARI) to represent the risk of aflatoxin contamination. The model was applied to simulate ARI for two commercial maize hybrids, H513 and H614D, grown in five multi-location field trials in Kenya using site specific agronomy, weather and soil parameters. The observed mean aflatoxin contamination in these trials varied from <1 to 7143 ppb. ARI simulated by the model explained 99% of the variation (p ≤ 0.001) in a linear relationship with the mean observed aflatoxin contamination. The strong relationship between ARI and aflatoxin contamination suggests that the model could be applied to map risk prone areas and to monitor in-season risk for genotypes and soils parameterized for APSIM.

An improved simulation model to predict pre-harvest aflatoxin risk in maize

Aflatoxin is a potent carcinogen produced by Aspergillus flavus, which frequently contaminates maize (Zea mays L.) in the field between 40° north and 40° south latitudes. A mechanistic model to predict risk of pre-harvest contamination could assist in management of this very harmful mycotoxin. In this study we describe an aflatoxin risk prediction model which is integrated with the Agricultural Production Systems Simulator (APSIM) modelling framework. The model computes a temperature function for A. flavus growth and aflatoxin production using a set of three cardinal temperatures determined in the laboratory using culture medium and intact grains. These cardinal temperatures were 11.5 °C as base, 32.5 °C as optimum and 42.5 °C as maximum. The model used a low (≤0.2) crop water supply to demand ratio—an index of drought during the grain filling stage to simulate maize crop's susceptibility to A. flavus growth and aflatoxin production. When this low threshold of the index was reached the model converted the temperature function into an aflatoxin risk index (ARI) to represent the risk of aflatoxin contamination. The model was applied to simulate ARI for two commercial maize hybrids, H513 and H614D, grown in five multi-location field trials in Kenya using site specific agronomy, weather and soil parameters. The observed mean aflatoxin contamination in these trials varied from <1 to 7143 ppb. ARI simulated by the model explained 99% of the variation (p ≤ 0.001) in a linear relationship with the mean observed aflatoxin contamination. The strong relationship between ARI and aflatoxin contamination suggests that the model could be applied to map risk prone areas and to monitor in-season risk for genotypes and soils parameterized for APSIM.