Sorghum Crop Modeling and Its Utility in Agronomy and Breeding

Crop models are simplified mathematical representations of the interacting biological and environmental components of the dynamic soil–plant–environment system. Sorghum crop modeling has evolved in parallel with crop modeling capability in general, since its origins in the 1960s and 1970s. Here we briefly review the trajectory in sorghum crop modeling leading to the development of advanced models. We then (i) overview the structure and function of the sorghum model in the Agricultural Production System sIMulator (APSIM) to exemplify advanced modeling concepts that suit both agronomic and breeding applications, (ii) review an example of use of sorghum modeling in supporting agronomic management decisions, (iii) review an example of the use of sorghum modeling in plant breeding, and (iv) consider implications for future roles of sorghum crop modeling. Modeling and simulation provide an avenue to explore consequences of crop management decision options in situations confronted with risks associated with seasonal climate uncertainties. Here we consider the possibility of manipulating planting configuration and density in sorghum as a means to manipulate the productivity–risk trade-off. A simulation analysis of decision options is presented and avenues for its use with decision-makers discussed. Modeling and simulation also provide opportunities to improve breeding efficiency by either dissecting complex traits to more amenable targets for genetics and breeding, or by trait evaluation via phenotypic prediction in target production regions to help prioritize effort and assess breeding strategies. Here we consider studies on the stay-green trait in sorghum, which confers yield advantage in water-limited situations, to exemplify both aspects. The possible future roles of sorghum modeling in agronomy and breeding are discussed as are opportunities related to their synergistic interaction. The potential to add significant value to the revolution in plant breeding associated with genomic technologies is identified as the new modeling frontier.

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.

High subsoil chloride concentrations reduce soil water extraction and crop yield on Vertosols in north-eastern Australia

Salinity, sodicity, acidity, and phytotoxic levels of chloride (Cl) in subsoils are major constraints to crop production in many soils of north-eastern Australia because they reduce the ability of crop roots to extract water and nutrients from the soil. The complex interactions and correlations among soil properties result in multi-colinearity between soil properties and crop yield that makes it difficult to determine which constraint is the major limitation. We used ridge-regression analysis to overcome colinearity to evaluate the contribution of soil factors and water supply to the variation in the yields of 5 winter crops on soils with various levels and combinations of subsoil constraints in the region. Subsoil constraints measured were soil Cl, electrical conductivity of the saturation extract (ECse), and exchangeable sodium percentage (ESP). The ridge regression procedure selected several of the variables used in a descriptive model, which included in-crop rainfall, plant-available soil water at sowing in the 0.90-1.10 m soil layer, and soil Cl in the 0.90-1.10 m soil layer, and accounted for 77-85% of the variation in the grain yields of the 5 winter crops. Inclusion of ESP of the top soil (0.0-0.10 m soil layer) marginally increased the descriptive capability of the models for bread wheat, barley and durum wheat. Subsoil Cl concentration was found to be an effective substitute for subsoil water extraction. The estimates of the critical levels of subsoil Cl for a 10% reduction in the grain yield were 492 mg cl/kg for chickpea, 662 mg Cl/kg for durum wheat, 854 mg Cl/kg for bread wheat, 980 mg Cl/kg for canola, and 1012 mg Cl/kg for barley, thus suggesting that chickpea and durum wheat were more sensitive to subsoil Cl than bread wheat, barley, and canola.

Long-term response of corn to limited irrigation and crop rotations

Dwindling water supplies for irrigation are prompting alternative management choices by irrigators. Limited irrigation, where less water is applied than full crop demand, may be a viable approach. Application of limited irrigation to corn was examined in this research. Corn was grown in crop rotations with dryland, limited irrigation, or full irrigation management from 1985 to 1999. Crop rotations included corn following corn (continuous corn), corn following wheat, followed by soybean (wheat-corn-soybean), and corn following soybean (corn-soybean). Full irrigation was managed to meet crop evapotranspiration requirements (ETc). Limited irrigation was managed with a seasonal target of no more than 150 mm applied. Precipitation patterns influenced the outcomes of measured parameters. Dryland yields had the most variation, while fully irrigated yields varied the least. Limited irrigation yields were 80% to 90%> of fully irrigated yields, but the limited irrigation plots received about half the applied water. Grain yields were significantly different among irrigation treatments. Yields were not significantly different among rotation treatments for all years and water treatments. For soil water parameters, more statistical differences were detected among the water management treatments than among the crop rotation treatments. Economic projections of these management practices showed that full irrigation produced the most income if water was available. Limited irrigation increased income significantly from dryland management.

Soil properties and crop yields in a dryland Vertisol sown with cotton-based crop rotations

Information on the effects of growing cotton (Gossypium hirsutum L.)-based crop rotations on soil quality of dryland Vertisols is sparse. The objective of this study was to quantify the effects of growing cereal and leguminous crops in rotation with dryland cotton on physical and chemical properties of a grey Vertisol near Warra, SE Queensland, Australia. The experimental treatments, selected after consultations with local cotton growers, were continuous cotton (T1), cotton-sorghum (Sorghum bicolor (L.) Moench.) (T2), cotton-wheat (Triticum aestivum L.) double cropped (T3), cotton-chickpea (Cicer arietinum L.) double cropped followed by wheat (T4) and cotton-wheat (T5). From 1993 to 1996 land preparation was by chisel ploughing to about 0.2 m followed by two to four cultivations with a Gyral tyne cultivator. Thereafter all crops were sown with zero tillage except for cultivation with a chisel plough to about 0.07-0.1 m after cotton picking to control heliothis moth pupae. Soil was sampled from 1996 to 2004 and physical (air-filled porosity of oven-dried soil, an indicator of soil compaction; plastic limit; linear shrinkage; dispersion index) and chemical (pH in 0.01 M CaCl2, organic carbon, exchangeable Ca, Mg, K and Na contents) properties measured. Crop rotation affected soil properties only with respect to exchangeable Na content and air-filled porosity. In the surface 0.15 m during 2000 and 2001 lowest air-filled porosity occurred with T1 (average of 34.6 m3/100 m3) and the highest with T3 (average of 38.9 m3/100 m3). Air-filled porosity decreased in the same depth between 1997 and 1998 from 45.0 to 36.1 m3/100 m3, presumably due to smearing and compaction caused by shallow cultivation in wet soil. In the subsoil, T1 and T2 frequently had lower air-filled porosity values in comparison with T3, T4 and T5, particularly during the early stages of the experiment, although values under T1 increased subsequently. In general, compaction was less under rotations which included a wheat crop (T3, T4, T5). For example, average air-filled porosity (in m3/100 m3) in the 0.15-0.30 m depth from 1996 to 1999 was 19.8 with both T1 and T2, and 21.2 with T3, 21.1 with T4 and 21.5 with T5. From 2000 to 2004, average air-filled porosity (in m3/100 m3) in the same depth was 21.3 with T1, 19.0 with T2, 19.8 with T3, 20.0 with T4 and 20.5 with T5. The rotation which included chickpea (T4) resulted in the lowest exchangeable Na content, although differences among rotations were small. Where only a cereal crop with a fibrous root system was sown in rotation with cotton (T2, T3, T5) linear shrinkage in the 0.45-0.60 m depth was lower than in rotations, which included tap-rooted crops such as chickpea (T4) or continuous cotton (T1). Dispersion index and organic carbon decreased, and plastic limit increased with time. Soil organic carbon stocks decreased at a rate of 1.2 Mg/ha/year. Lowest average cotton lint yield occurred with T2 (0.54 Mg/ha) and highest wheat yield with T3 (2.8 Mg/ha). Rotations which include a wheat crop are more likely to result in better soil structure and cotton lint yield than cotton-sorghum or continuous cotton.

Physiological Basis for Genotypic Variation in Tolerance to and Recovery from Pre-flowering Drought in Peanut.

Drought during the pre-flowering stage can increase yield of peanut. There is limited information on genotypic variation for tolerance to and recovery from pre-flowering drought (PFD) and more importantly the physiological traits underlying genotypic variation. The objectives of this study were to determine the effects of moisture stress during the pre-flowering phase on pod yield and to understand some of the physiological responses underlying genotypic variation in response to and recovery from PFD. A glasshouse and field experiments were conducted at Khon Kaen University, Thailand. The glasshouse experiment was a randomized complete block design consisting of two watering regimes, i.e. fully-irrigated control and 1/3 available soil water from emergence to 40 days after emergence followed by adequate water supply, and 12 peanut genotypes. The field experiment was a split-plot design with two watering regimes as main-plots, and 12 peanut genotypes as sub-plots. Measurements of N-2 fixation, leaf area (LA) were made in both experiments. In addition, root growth was measured in the glasshouse experiment. Imposition of PFD followed by recovery resulted in an average increase in yield of 24 % (range from 10 % to 57 %) and 12 % (range from 2 % to 51 %) in the field and glasshouse experiments, respectively. Significant genotypic variation for N-2 fixation, LA and root growth was also observed after recovery. The study revealed that recovery growth following release of PFD had a stronger influence on final yield than tolerance to water deficits during the PFD. A combination of N-2 fixation, LA and root growth accounted for a major portion of the genotypic variation in yield (r = 0.68-0.93) suggesting that these traits could be used as selection criteria for identifying genotypes with rapid recovery from PFD. A combined analysis of glasshouse and field experiments showed that LA and N-2 fixation during the recovery had low genotype x environment interaction indicating potential for using these traits for selecting genotypes in peanut improvement programs.

Weed growth and crop yield loss in wheat as influenced by row spacing and weed emergence times

Reducing crop row spacing and delaying time of weed emergence may provide crops a competitive edge over weeds. Field experiments were conducted to evaluate the effects of crop row spacing (11, 15, and 23-cm) and weed emergence time (0, 20, 35, 45, 55, and 60 days after wheat emergence; DAWE) on Galium aparine and Lepidium sativum growth and wheat yield losses. Season-long weed-free and crop-free treatments were also established to compare wheat yield and weed growth, respectively. Row spacing and weed emergence time significantly affected the growth of both weed species and wheat grain yields. For both weed species, the maximum plant height, shoot biomass, and seed production were observed in the crop-free plots, and delayed emergence decreased these variables. In weed-crop competition plots, maximum weed growth was observed when weeds emerged simultaneously with the crop in rows spaced 23-cm apart. Less growth of both weed species was observed in narrow row spacing (11-cm) of wheat as compared with wider rows (15 and 23-cm). These weed species produced less than 5 seeds plant-1 in 11-cm wheat rows when they emerged at 60 DAWE. Presence of weeds in the crop especially at early stages was devastating for wheat yields. Therefore, maximum grain yield (4.91tha-1) was recorded in the weed-free treatment at 11-cm row spacing. Delay in time of weed emergence and narrow row spacing reduced weed growth and seed production and enhanced wheat grain yield, suggesting that these strategies could contribute to weed management in wheat.

Soybean rotation and crop residue management to reduce nitrous oxide emissions from sugarcane soils

NITROUS OXIDE (N2O) IS a potent greenhouse gas and the predominant ozone-depleting substance in the atmosphere. Agricultural nitrogenous fertiliser use is the major source of human-induced N2O emissions. A field experiment was conducted at Bundaberg from October 2012 to September 2014 to examine the impacts of legume crop (soybean) rotation as an alternative nitrogen (N) source on N2O emissions during the fallow period and to investigate low-emission soybean residue management practices. An automatic monitoring system and manual gas sampling chambers were used to measure greenhouse gas emissions from soil. Soybean cropping during the fallow period reduced N2O emissions compared to the bare fallow. Based on the N content in the soybean crop residues, the fertiliser N application rate was reduced by about 120 kg N/ha for the subsequent sugarcane crop. Consequently, emissions of N2O during the sugarcane cropping season were significantly lower from the soybean cropped soil than those from the conventionally fertilised (145 kg N/ha) soil following bare fallow. However, tillage that incorporated the soybean crop residues into soil promoted N2O emissions in the first two months. Spraying a nitrification inhibitor (DMPP) onto the soybean crop residues before tillage effectively prevented the N2O emission spikes. Compared to conventional tillage, practising no-till with or without growing a nitrogen catch crop during the time after soybean harvest and before cane planting also reduced N2O emissions substantially. These results demonstrated that soybean rotation during the fallow period followed with N conservation management practices could offer a promising N2O mitigation strategy in sugarcane farming. Further investigation is required to provide guidance on N and water management following soybean fallow to maintain sugar productivity.

Crop sequences to manage soil pathogens and increase northern grain production

This is part of a GRDC funded project led by Dr Jeremy Whish of CSIRO Ecosystem Sciences. The project aims to build a root-lesion nematode module into the crop growth simulation program APSIM (Agricultural Production Systems Simulator). This will utilise existing nematode and crop data from field, glasshouse and laboratory research led by Dr John Thompson. New data will be collected to validate and extend the model.

Subsoil constraints in Vertosols: Crop water use nutrient concentration, and grain yields of bread wheat, durum wheat, barley, chickpea and canola

Single or multiple factors implicated in subsoil constraints including salinity, sodicity, and phytotoxic concentrations of chloride (Cl) are present in many Vertosols including those occurring in Queensland, Australia. The variable distribution and the complex interactions that exist between these constraints limit the agronomic or management options available to manage the soil with these subsoil constraints. The identification of crops and cultivars adapted to these adverse subsoil conditions and/or able to exploit subsoil water may be an option to maintain productivity of these soils. We evaluated relative performance of 5 winter crop species, in terms of grain yields, nutrient concentration, and ability to extract soil water, grown on soils with various levels and combinations of subsoil constraints in 19 field experiments over 2 years. Subsoil constraints were measured by levels of soil Cl, electrical conductivity of the saturation extract (ECse), and exchangeable sodium percentage (ESP). Increasing levels of subsoil constraints significantly decreased maximum depth of water extraction, grain yield, and plant-available water capacity for all the 5 crops and more so for chickpea and durum wheat than bread wheat, barley, or canola. Increasing soil Cl levels had a greater restricting effect on water availability than did ECse and ESP. We developed empirical relationships between soil Cl, ECse, and ESP and crop lower limit (CLL) for estimating subsoil water extraction by 5 winter crops. However, the presence of gypsum influenced the ability to predict CLL based on the levels of ECse. Stronger relationships between apparent unused plant-available water (CLL - LL15; LL15 is lower limit at -1.5 MPa) and soil Cl concentrations than ESP or ECse suggested that the presence of high Cl in these soils most likely inhibited the subsoil water extraction by the crops. This was supported by increased sodium (Na) and Cl concentration with a corresponding decrease in calcium (Ca) and potassium (K) in young mature leaf of bread wheat, durum wheat, and chickpea with increasing levels of subsoil constraints. Of the 2 ions, Na and Cl, the latter appears to be more damaging than the former, resulting in plant dieback and reduced grain yields.

Assessment of functional forms of crop yield loss models of invasive plant species applied in decision support tools and bioeconomic modelling

AbstractObjectives Decision support tools (DSTs) for invasive species management have had limited success in producing convincing results and meeting users' expectations. The problems could be linked to the functional form of model which represents the dynamic relationship between the invasive species and crop yield loss in the DSTs. The objectives of this study were: a) to compile and review the models tested on field experiments and applied to DSTs; and b) to do an empirical evaluation of some popular models and alternatives. Design and methods This study surveyed the literature and documented strengths and weaknesses of the functional forms of yield loss models. Some widely used models (linear, relative yield and hyperbolic models) and two potentially useful models (the double-scaled and density-scaled models) were evaluated for a wide range of weed densities, maximum potential yield loss and maximum yield loss per weed. Results Popular functional forms include hyperbolic, sigmoid, linear, quadratic and inverse models. Many basic models were modified to account for the effect of important factors (weather, tillage and growth stage of crop at weed emergence) influencing weed–crop interaction and to improve prediction accuracy. This limited their applicability for use in DSTs as they became less generalized in nature and often were applicable to a much narrower range of conditions than would be encountered in the use of DSTs. These factors' effects could be better accounted by using other techniques. Among the model empirically assessed, the linear model is a very simple model which appears to work well at sparse weed densities, but it produces unrealistic behaviour at high densities. The relative-yield model exhibits expected behaviour at high densities and high levels of maximum yield loss per weed but probably underestimates yield loss at low to intermediate densities. The hyperbolic model demonstrated reasonable behaviour at lower weed densities, but produced biologically unreasonable behaviour at low rates of loss per weed and high yield loss at the maximum weed density. The density-scaled model is not sensitive to the yield loss at maximum weed density in terms of the number of weeds that will produce a certain proportion of that maximum yield loss. The double-scaled model appeared to produce more robust estimates of the impact of weeds under a wide range of conditions. Conclusions Previously tested functional forms exhibit problems for use in DSTs for crop yield loss modelling. Of the models evaluated, the double-scaled model exhibits desirable qualitative behaviour under most circumstances.