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.

Domestication to crop improvement: Genetic resources for Sorghum and Saccharum (Andropogoneae).

Background: Both sorghum (Sorghum bicolor) and sugarcane (Saccharum officinarum) are members of the Andropogoneae tribe in the Poaceae and are each other's closest relatives amongst cultivated plants. Both are relatively recent domesticates and comparatively little of the genetic potential of these taxa and their wild relatives has been captured by breeding programmes to date. This review assesses the genetic gains made by plant breeders since domestication and the progress in the characterization of genetic resources and their utilization in crop improvement for these two related species. Genetic Resources: The genome of sorghum has recently been sequenced providing a great boost to our knowledge of the evolution of grass genomes and the wealth of diversity within S. bicolor taxa. Molecular analysis of the Sorghum genus has identified close relatives of S. bicolor with novel traits, endosperm structure and composition that may be used to expand the cultivated gene pool. Mutant populations (including TILLING populations) provide a useful addition to genetic resources for this species. Sugarcane is a complex polyploid with a large and variable number of copies of each gene. The wild relatives of sugarcane represent a reservoir of genetic diversity for use in sugarcane improvement. Techniques for quantitative molecular analysis of gene or allele copy number in this genetically complex crop have been developed. SNP discovery and mapping in sugarcane has been advanced by the development of high-throughput techniques for ecoTILLING in sugarcane. Genetic linkage maps of the sugarcane genome are being improved for use in breeding selection. The improvement of both sorghum and sugarcane will be accelerated by the incorporation of more diverse germplasm into the domesticated gene pools using molecular tools and the improved knowledge of these genomes.

Amelioration of soil compaction can take 5 years on a Vertisol under no till in the semi-arid subtropics

Heavy wheel traffic causes soil compaction, which adversely affects crop production and may persist for several years. We applied known compaction forces to entire plots annually for 5 years, and then determined the duration of the adverse effects on the properties of a Vertisol and the performance of crops under no-till dryland cropping with residue retention. For up to 5 years after a final treatment with a 10 Mg axle load on wet soil, soil shear strength at 70-100 mm and cone index at 180-360 mm were significantly (P < 0.05) higher than in a control treatment, and soil water storage and grain yield were lower. We conclude that compaction effects persisted because (1) there were insufficient wet-dry cycles to swell and shrink the entire compacted layer, (2) soil loosening by tillage was absent and (3) there were fewer earthworms in the compacted soil. Compaction of dry soil with 6 Mg had little effect at any time, indicating that by using wheel traffic only when the soil is dry, problems can be avoided. Unfortunately such a restriction is not always possible because sowing, tillage and harvest operations often need to be done when the soil is wet. A more generally applicable solution, which also ensures timely operations, is the permanent separation of wheel zones and crop zones in the field--the practice known as controlled traffic farming. Where a compacted layer already exists, even on a clay soil, management options to hasten repair should be considered, e.g. tillage, deep ripping, sowing a ley pasture or sowing crop species more effective at repairing compacted soil.

Long-term impacts of harvest residue management on nutrition, growth and productivity of an exotic pine plantation of sub-tropical Australia

Residue retention is an important issue in evaluating the sustainability of production forestry. However, its long-term impacts have not been studied extensively, especially in sub-tropical environments. This study investigated the long-term impact of harvest residue retention on tree nutrition, growth and productivity of a F1 hybrid (Pinus elliottii var. elliottii × Pinus caribaea var. hondurensis) exotic pine plantation in sub-tropical Australia, under three harvest residue management regimes: (1) residue removal, RR0; (2) single residue retention, RR1; and (3) double residue retention, RR2. The experiment, established in 1996, is a randomised complete block design with 4 replicates. Tree growth measurements in this study were carried out at ages 2, 4, 6, 8 and 10 years, while foliar nutrient analyses were carried out at ages 2, 4, 6 and 10 years. Litter production and litter nitrogen (N) and phosphorus (P) measurements were carried out quarterly over a 15-month period between ages 9 and 10 years. Results showed that total tree growth was still greater in residue-retained treatments compared to the RR0 treatment. However, mean annual increments of diameter at breast height (MAID) and basal area (MAIB) declined significantly after age 4 years to about 68-78% at age 10 years. Declining foliar N and P concentrations accounted for 62% (p < 0.05) of the variation of growth rates after age 4 years, and foliar N and P concentrations were either marginal or below critical concentrations. In addition, litter production, and litter N and P contents were not significantly different among the treatments. This study suggests that the impact of residue retention on tree nutrition and growth rates might be limited over a longer period, and that the integration of alternative forest management practices is necessary to sustain the benefits of harvest residues until the end of the rotation.

Crop and environmental attributes underpinning genotype by environment interaction in synthetic-derived bread wheat evaluated in Mexico and Australia

Synthetic backcrossed-derived bread wheats (SBWs) from CIMMYT were grown in the north-west of Mexico (CIANO) and sites across Australia during 3 seasons. A different set of lines was evaluated each season, as new materials became available from the CIMMYT crop enhancement program. Previously, we have evaluated both the performance of genotypes across environments and the genotype x environment interaction (G x E). The objective of this study was to interpret the G x E for yield in terms of crop attributes measured at individual sites and to identify the potential environmental drivers of this interaction. Groups of SBWs with consistent yield performance were identified, often comprising closely related lines. However, contrasting performance was also relatively common among sister lines or between a recurrent parent and its SBWs. Early flowering was a common feature among lines with broad adaptation and/or high yield in the northern Australian wheatbelt, while yields in the southern region did not show any association with the maturity type. Lines with high yields in the southern and northern regions had cooler canopies during flowering and early grain filling. Among the SBWs with Australian genetic backgrounds, lines best adapted to CIANO were tall (>100 cm), with a slightly higher ground cover. These lines also displayed a higher concentration of water-soluble carbohydrates in the stem at flowering, which was negatively correlated with stem number per unit area when evaluated in southern Australia (Horsham). Possible reasons for these patterns are discussed. Selection for yield at CIANO did not specifically identify the lines best adapted to northern Australia, although they were not the most poorly adapted either. In addition, groups of lines with specific adaptation to the south would not have been selected by choosing the highest yielding lines at CIANO. These findings suggest that selection at CIMMYT for Australian environments may be improved by either trait based selection or yield data combined with trait information. Flowering date, canopy temperature around flowering, tiller density, and water-soluble carbohydrate concentration in the stem at flowering seem likely candidates.

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.

Analysis of high yielding maize production - A study based on a commercial crop

This paper reports on the use of APSIM - Maize for retrospective analysis of performance of a high input, high yielding maize crop and analysis of predicted performance of maize grown with high inputs over the long-term (>100 years) for specified scenarios of environmental conditions (temperature and radiation) and agronomic inputs (sowing date, plant population, nitrogen fertiliser and irrigation) at Boort, Victoria, Australia. It uses a high yielding (17 400 kg/ha dry grain, 20 500 kg/ha at 15% water) commercial crop grown in 2004-05 as the basis of the study. Yield for the agronomic and environmental conditions of 2004-05 was predicted accurately, giving confidence that the model could be used for the detailed analyses undertaken. The analysis showed that the yield achieved was close to that possible with the conditions and agronomic inputs of 2004-05. Sowing dates during 21 September to 26 October had little effect on predicted yield, except when combined with reduced temperature. Single year and long-term analyses concluded that a higher plant population (11 plants/m2) is needed to optimise yield, but that slightly lower N and irrigation inputs are appropriate for the plant population used commercially (8.4 plants/m2). Also, compared with changes in agronomic inputs increases in temperature and/or radiation had relatively minor effects, except that reduced temperature reduces predicted yield substantially. This study provides an approach for the use of models for both retrospective analysis of crop performance and assessment of long-term variability of crop yield under a wide range of agronomic and environmental conditions.

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.

Detection of nitrogen deficiency in wheat from spectral reflectance indices and basic crop eco-physiological concepts

We tested the capacity of several published multispectral indices to estimate the nitrogen nutrition of wheat canopies grown under different levels of water supply and plant density and derived a simple canopy reflectance index that is greatly independent of those factors. Planar domain geometry was used to account for mixed signals from the canopy and soil when the ground cover was low. A nitrogen stress index was developed, which adjusts shoot %N for plant biomass and area, thereby accounting for environmental conditions that affect growth, such as crop water status. The canopy chlorophyll content index (CCCi) and the modified spectral ratio planar index (mSRPi) could explain 68 and 69% of the observed variability in the nitrogen nutrition of the crop as early as Zadoks 33, irrespective of water status or ground cover. The CCCi was derived from the combination of 3 wavebands 670, 720 and 790 nm, and the mSRPi from 445, 705 and 750 nm, together with broader bands in the NIR and RED. The potential for their spatial application over large fields/paddocks is discussed.

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.

Detecting the attributes of a wheat crop using digital imagery acquired from a low-altitude platform

A low-altitude platform utilising a 1.8-m diameter tethered helium balloon was used to position a multispectral sensor, consisting of two digital cameras, above a fertiliser trial plot where wheat (Triticum spp.) was being grown. Located in Cecil Plains, Queensland, Australia, the plot was a long-term fertiliser trial being conducted by a fertiliser company to monitor the response of crops to various levels of nutrition. The different levels of nutrition were achieved by varying nitrogen application rates between 0 and 120 units of N at 40 unit increments. Each plot had received the same application rate for 10 years. Colour and near-infrared images were acquired that captured the whole 2 ha plot. These images were examined and relationships sought between the captured digital information and the crop parameters imaged at anthesis and the at-harvest quality and quantity parameters. The statistical analysis techniques used were correlation analysis, discriminant analysis and partial least squares regression. A high correlation was found between the image and yield (R2 = 0.91) and a moderate correlation between the image and grain protein content (R2 = 0.66). The utility of the system could be extended by choosing a more mobile platform. This would increase the potential for the system to be used to diagnose the causes of the variability and allow remediation, and/or to segregate the crop at harvest to meet certain quality parameters.

Annual forecasting of the Australian macadamia crop - integrating tree census data with statistical climate-adjustment models

To facilitate marketing and export, the Australian macadamia industry requires accurate crop forecasts. Each year, two levels of crop predictions are produced for this industry. The first is an overall longer-term forecast based on tree census data of growers in the Australian Macadamia Society (AMS). This data set currently accounts for around 70% of total production, and is supplemented by our best estimates of non-AMS orchards. Given these total tree numbers, average yields per tree are needed to complete the long-term forecasts. Yields from regional variety trials were initially used, but were found to be consistently higher than the average yields that growers were obtaining. Hence, a statistical model was developed using growers' historical yields, also taken from the AMS database. This model accounted for the effects of tree age, variety, year, region and tree spacing, and explained 65% of the total variation in the yield per tree data. The second level of crop prediction is an annual climate adjustment of these overall long-term estimates, taking into account the expected effects on production of the previous year's climate. This adjustment is based on relative historical yields, measured as the percentage deviance between expected and actual production. The dominant climatic variables are observed temperature, evaporation, solar radiation and modelled water stress. Initially, a number of alternate statistical models showed good agreement within the historical data, with jack-knife cross-validation R2 values of 96% or better. However, forecasts varied quite widely between these alternate models. Exploratory multivariate analyses and nearest-neighbour methods were used to investigate these differences. For 2001-2003, the overall forecasts were in the right direction (when compared with the long-term expected values), but were over-estimates. In 2004 the forecast was well under the observed production, and in 2005 the revised models produced a forecast within 5.1% of the actual production. Over the first five years of forecasting, the absolute deviance for the climate-adjustment models averaged 10.1%, just outside the targeted objective of 10%.

Laboratory information management software for genotyping workflows: Applications in high throughput crop genotyping

Background: With the advances in DNA sequencer-based technologies, it has become possible to automate several steps of the genotyping process leading to increased throughput. To efficiently handle the large amounts of genotypic data generated and help with quality control, there is a strong need for a software system that can help with the tracking of samples and capture and management of data at different steps of the process. Such systems, while serving to manage the workflow precisely, also encourage good laboratory practice by standardizing protocols, recording and annotating data from every step of the workflow Results: A laboratory information management system (LIMS) has been designed and implemented at the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT) that meets the requirements of a moderately high throughput molecular genotyping facility. The application is designed as modules and is simple to learn and use. The application leads the user through each step of the process from starting an experiment to the storing of output data from the genotype detection step with auto-binning of alleles; thus ensuring that every DNA sample is handled in an identical manner and all the necessary data are captured. The application keeps track of DNA samples and generated data. Data entry into the system is through the use of forms for file uploads. The LIMS provides functions to trace back to the electrophoresis gel files or sample source for any genotypic data and for repeating experiments. The LIMS is being presently used for the capture of high throughput SSR (simple-sequence repeat) genotyping data from the legume (chickpea, groundnut and pigeonpea) and cereal (sorghum and millets) crops of importance in the semi-arid tropics. Conclusions: A laboratory information management system is available that has been found useful in the management of microsatellite genotype data in a moderately high throughput genotyping laboratory. The application with source code is freely available for academic users and can be downloaded from http://www.icrisat.org/bt-software-d-lims.htm

Estimating crop area using seasonal time series of enhanced vegetation index from MODIS satellite imagery

Cereal grain is one of the main export commodities of Australian agriculture. Over the past decade, crop yield forecasts for wheat and sorghum have shown appreciable utility for industry planning at shire, state, and national scales. There is now an increasing drive from industry for more accurate and cost-effective crop production forecasts. In order to generate production estimates, accurate crop area estimates are needed by the end of the cropping season. Multivariate methods for analysing remotely sensed Enhanced Vegetation Index (EVI) from 16-day Moderate Resolution Imaging Spectroradiometer (MODIS) satellite imagery within the cropping period (i.e. April-November) were investigated to estimate crop area for wheat, barley, chickpea, and total winter cropped area for a case study region in NE Australia. Each pixel classification method was trained on ground truth data collected from the study region. Three approaches to pixel classification were examined: (i) cluster analysis of trajectories of EVI values from consecutive multi-date imagery during the crop growth period; (ii) harmonic analysis of the time series (HANTS) of the EVI values; and (iii) principal component analysis (PCA) of the time series of EVI values. Images classified using these three approaches were compared with each other, and with a classification based on the single MODIS image taken at peak EVI. Imagery for the 2003 and 2004 seasons was used to assess the ability of the methods to determine wheat, barley, chickpea, and total cropped area estimates. The accuracy at pixel scale was determined by the percent correct classification metric by contrasting all pixel scale samples with independent pixel observations. At a shire level, aggregated total crop area estimates were compared with surveyed estimates. All multi-temporal methods showed significant overall capability to estimate total winter crop area. There was high accuracy at pixel scale (>98% correct classification) for identifying overall winter cropping. However, discrimination among crops was less accurate. Although the use of single-date EVI data produced high accuracy for estimates of wheat area at shire scale, the result contradicted the poor pixel-scale accuracy associated with this approach, due to fortuitous compensating errors. Further studies are needed to extrapolate the multi-temporal approaches to other geographical areas and to improve the lead time for deriving cropped-area estimates before harvest.