Estimation of light interception in research environments: A joint approach using directional light sensors and 3D virtual plants applied to sunflower (Helianthus annuus) and Arabidopsis thaliana in natural and artificial conditions

Light interception is a major factor influencing plant development and biomass production. Several methods have been proposed to determine this variable, but its calculation remains difficult in artificial environments with heterogeneous light. We propose a method that uses 3D virtual plant modelling and directional light characterisation to estimate light interception in highly heterogeneous light environments such as growth chambers and glasshouses. Intercepted light was estimated by coupling an architectural model and a light model for different genotypes of the rosette species Arabidopsis thaliana (L.) Heynh and a sunflower crop. The model was applied to plants of contrasting architectures, cultivated in isolation or in canopy, in natural or artificial environments, and under contrasting light conditions. The model gave satisfactory results when compared with observed data and enabled calculation of light interception in situations where direct measurements or classical methods were inefficient, such as young crops, isolated plants or artificial conditions. Furthermore, the model revealed that A. thaliana increased its light interception efficiency when shaded. To conclude, the method can be used to calculate intercepted light at organ, plant and plot levels, in natural and artificial environments, and should be useful in the investigation of genotype-environment interactions for plant architecture and light interception efficiency. This paper originates from a presentation at the 5th International Workshop on Functional–Structural Plant Models, Napier, New Zealand, November 2007.

The role of root architectural traits in adaptation of wheat to water-limited environments

Better understanding of root system structure and function is critical to crop improvement in water-limited environments. The aims of this study were to examine root system characteristics of two wheat genotypes contrasting in tolerance to water limitation and to assess the functional implications on adaptation to water-limited environments of any differences found. The drought tolerant barley variety, Mackay, was also included to allow inter-species comparison. Single plants were grown in large, soil-filled root-observation chambers. Root growth was monitored by digital imaging and water extraction was measured. Root architecture differed markedly among the genotypes. The drought-tolerant wheat (cv. SeriM82) had a compact root system, while roots of barley cv. Mackay occupied the largest soil volume. Relative to the standard wheat variety (Hartog), SeriM82 had a more uniform rooting pattern and greater root length at depth. Despite the more compact root architecture of SeriM82, total water extracted did not differ between wheat genotypes. To quantify the value of these adaptive traits, a simulation analysis was conducted with the cropping system model APSIM, for a wide range of environments in southern Queensland, Australia. The analysis indicated a mean relative yield benefit of 14.5% in water-deficit seasons. Each additional millimetre of water extracted during grain filling generated an extra 55 kg ha-1 of grain yield. The functional implications of root traits on temporal patterns and total amount of water capture, and their importance in crop adaptation to specific water-limited environments, are discussed.

Detection and use of QTL for complex traits in multiple environments.

QTL mapping methods for complex traits are challenged by new developments in marker technology, phenotyping platforms, and breeding methods. In meeting these challenges, QTL mapping approaches will need to also acknowledge the central roles of QTL by environment interactions (QEI) and QTL by trait interactions in the expression of complex traits like yield. This paper presents an overview of mixed model QTL methodology that is suitable for many types of populations and that allows predictive modeling of QEI, both for environmental and developmental gradients. Attention is also given to multi-trait QTL models which are essential to interpret the genetic basis of trait correlations. Biophysical (crop growth) model simulations are proposed as a complement to statistical QTL mapping for the interpretation of the nature of QEI and to investigate better methods for the dissection of complex traits into component traits and their genetic controls.

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.

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.

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.

Using simulation techniques to investigate methods to determine resistance of helminths to anthelmintic treatment

The widespread and increasing resistance of internal parasites to anthelmintic control is a serious problem for the Australian sheep and wool industry. As part of control programmes, laboratories use the Faecal Egg Count Reduction Test (FECRT) to determine resistance to anthelmintics. It is important to have confidence in the measure of resistance, not only for the producer planning a drenching programme but also for companies investigating the efficacy of their products. The determination of resistance and corresponding confidence limits as given in anthelmintic efficacy guidelines of the Standing Committee on Agriculture (SCA) is based on a number of assumptions. This study evaluated the appropriateness of these assumptions for typical data and compared the effectiveness of the standard FECRT procedure with the effectiveness of alternative procedures. Several sets of historical experimental data from sheep and goats were analysed to determine that a negative binomial distribution was a more appropriate distribution to describe pre-treatment helminth egg counts in faeces than a normal distribution. Simulated egg counts for control animals were generated stochastically from negative binomial distributions and those for treated animals from negative binomial and binomial distributions. Three methods for determining resistance when percent reduction is based on arithmetic means were applied. The first was that advocated in the SCA guidelines, the second similar to the first but basing the variance estimates on negative binomial distributions, and the third using Wadley’s method with the distribution of the response variate assumed negative binomial and a logit link transformation. These were also compared with a fourth method recommended by the International Co-operation on Harmonisation of Technical Requirements for Registration of Veterinary Medicinal Products (VICH) programme, in which percent reduction is based on the geometric means. A wide selection of parameters was investigated and for each set 1000 simulations run. Percent reduction and confidence limits were then calculated for the methods, together with the number of times in each set of 1000 simulations the theoretical percent reduction fell within the estimated confidence limits and the number of times resistance would have been said to occur. These simulations provide the basis for setting conditions under which the methods could be recommended. The authors show that given the distribution of helminth egg counts found in Queensland flocks, the method based on arithmetic not geometric means should be used and suggest that resistance be redefined as occurring when the upper level of percent reduction is less than 95%. At least ten animals per group are required in most circumstances, though even 20 may be insufficient where effectiveness of the product is close to the cut off point for defining resistance.

Spectral and thermal sensing for nitrogen and water status in rainfed and irrigated wheat environments

Variable-rate technologies and site-specific crop nutrient management require real-time spatial information about the potential for response to in-season crop management interventions. Thermal and spectral properties of canopies can provide relevant information for non-destructive measurement of crop water and nitrogen stresses. In previous studies, foliage temperature was successfully estimated from canopy-scale (mixed foliage and soil) temperatures and the multispectral Canopy Chlorophyll Content Index (CCCI) was effective in measuring canopy-scale N status in rainfed wheat (Triticum aestivum L.) systems in Horsham, Victoria, Australia. In the present study, results showed that under irrigated wheat systems in Maricopa, Arizona, USA, the theoretical derivation of foliage temperature unmixing produced relationships similar to those in Horsham. Derivation of the CCCI led to an r2 relationship with chlorophyll a of 0.53 after Zadoks stage 43. This was later than the relationship (r2 = 0.68) developed for Horsham after Zadoks stage 33 but early enough to be used for potential mid-season N fertilizer recommendations. Additionally, ground-based hyperspectral data estimated plant N (g kg)1) in Horsham with an r2 = 0.86 but was confounded by water supply and N interactions. By combining canopy thermal and spectral properties, varying water and N status can potentially be identified eventually permitting targeted N applications to those parts of a field where N can be used most efficiently by the crop.

Yield of synthetic backcross-derived lines in rainfed environments of Australia

Wheat is one of the major food crops in the world. It is Australia's largest crop and most important agricultural commodity. In Australia the crop is grown under rainfed conditions with inherently important regional environmental differences; wheat growing areas are characterized by winter dominant rainfall in southern and western Australia and summer rainfall in northern Australia. Maximizing yield potential across these diverse regions is dependent upon managing, either genetically or agronomically, those factors in the environment that limit yield. The potential of synthetic backcross lines (SBLs) to increase yield in the diverse agroecological zones of Australia was investigated. Significant yield advantages were found for many of the SBLs across diverse environments. Depending on the environment, the yield of the SBLs ranged from 8% to 30% higher than the best local check in Australia. Apart from adaptation to semiarid water stressed conditions, some SBLs were also found to be significantly higher yielding under more optimal (irrigated) conditions. The four testing environments were classified into two groups, with the northern and southern environments being in separate groups. An elite group of SBLs was identified that exhibited broad adaptation across all diverse Australian environments included in this study. Other SBLs showed specific adaptation to either northern or southern Australia. This study showed that SBLs are likely to provide breeders with the opportunity to significantly improve wheat yield beyond what was previously possible in a number of diverse production environments.

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.

CIMMYT-selected derived synthetic bread wheats for rainfed environments: Yield evaluation in Mexico and Australia.

Synthetic backcrossed-derived bread wheats (SBWs) from CIMMYT were grown in the Northwest of Mexico at Centro de Investigaciones Agrícolas del Noroeste (CIANO) and sites across Australia during three seasons. During three consecutive years Australia received “shipments” of different SBWs from CIMMYT for evaluation. A different set of lines was evaluated each season, as new materials became available from the CIMMYT crop enhancement program. These consisted of approximately 100 advanced lines (F7) per year. SBWs had been top and backcrossed to CIMMYT cultivars in the first two shipments and to Australian wheat cultivars in the third one. At CIANO, the SBWs were trialled under receding soil moisture conditions. We evaluated both the performance of each line across all environments and the genotype-by-environment interaction using an analysis that fits a multiplicative mixed model, adjusted for spatial field trends. Data were organised in three groups of multienvironment trials (MET) containing germplasm from shipment 1 (METShip1), 2 (METShip2), and 3 (METShip3), respectively. Large components of variance for the genotype × environment interaction were found for each MET analysis, due to the diversity of environments included and the limited replication over years (only in METShip2, lines were tested over 2 years). The average percentage of genetic variance explained by the factor analytic models with two factors was 50.3% for METShip1, 46.7% for METShip2, and 48.7% for METShip3. Yield comparison focused only on lines that were present in all locations within a METShip, or “core” SBWs. A number of core SBWs, crossed to both Australian and CIMMYT backgrounds, outperformed the local benchmark checks at sites from the northern end of the Australian wheat belt, with reduced success at more southern locations. In general, lines that succeeded in the north were different from those in the south. The moderate positive genetic correlation between CIANO and locations in the northern wheat growing region likely reflects similarities in average temperature during flowering, high evaporative demand, and a short flowering interval. We are currently studying attributes of this germplasm that may contribute to adaptation, with the aim of improving the selection process in both Mexico and Australia.

A comparison of plotless density estimators using Monte Carlo simulation on totally enumerated field data sets

Background: Plotless density estimators are those that are based on distance measures rather than counts per unit area (quadrats or plots) to estimate the density of some usually stationary event, e.g. burrow openings, damage to plant stems, etc. These estimators typically use distance measures between events and from random points to events to derive an estimate of density. The error and bias of these estimators for the various spatial patterns found in nature have been examined using simulated populations only. In this study we investigated eight plotless density estimators to determine which were robust across a wide range of data sets from fully mapped field sites. They covered a wide range of situations including animal damage to rice and corn, nest locations, active rodent burrows and distribution of plants. Monte Carlo simulations were applied to sample the data sets, and in all cases the error of the estimate (measured as relative root mean square error) was reduced with increasing sample size. The method of calculation and ease of use in the field were also used to judge the usefulness of the estimator. Estimators were evaluated in their original published forms, although the variable area transect (VAT) and ordered distance methods have been the subjects of optimization studies. Results: An estimator that was a compound of three basic distance estimators was found to be robust across all spatial patterns for sample sizes of 25 or greater. The same field methodology can be used either with the basic distance formula or the formula used with the Kendall-Moran estimator in which case a reduction in error may be gained for sample sizes less than 25, however, there is no improvement for larger sample sizes. The variable area transect (VAT) method performed moderately well, is easy to use in the field, and its calculations easy to undertake. Conclusion: Plotless density estimators can provide an estimate of density in situations where it would not be practical to layout a plot or quadrat and can in many cases reduce the workload in the field.

Simulating pasture growth rates in Australian and New Zealand grazing systems

DairyMod, EcoMod, and the SGS Pasture Model are mechanistic biophysical models developed to explore scenarios in grazing systems. The aim of this manuscript was to test the ability of the models to simulate net herbage accumulation rates of ryegrass-based pastures across a range of environments and pasture management systems in Australia and New Zealand. Measured monthly net herbage accumulation rate and accumulated yield data were collated from ten grazing system experiments at eight sites ranging from cool temperate to subtropical environments. The local climate, soil, pasture species, and management (N fertiliser, irrigation, and grazing or cutting pattern) were described in the model for each site, and net herbage accumulation rates modelled. The model adequately simulated the monthly net herbage accumulation rates across the range of environments, based on the summary statistics and observed patterns of seasonal growth, particularly when the variability in measured herbage accumulation rates was taken into account. Agreement between modelled and observed growth rates was more accurate and precise in temperate than in subtropical environments, and in winter and summer than in autumn and spring. Similarly, agreement between predicted and observed accumulated yields was more accurate than monthly net herbage accumulation. Different temperature parameters were used to describe the growth of perennial ryegrass cultivars and annual ryegrass; these differences were in line with observed growth patterns and breeding objectives. Results are discussed in the context of the difficulties in measuring pasture growth rates and model limitations.

Relative contributions of light interception and radiation use efficiency to the reduction of maize productivity under cold temperatures

Maize (Zea mays L.) is a chill-susceptible crop cultivated in northern latitude environments. The detrimental effects of cold on growth and photosynthetic activity have long been established. However, a general overview of how important these processes are with respect to the reduction of productivity reported in the field is still lacking. In this study, a model-assisted approach was used to dissect variations in productivity under suboptimal temperatures and quantify the relative contributions of light interception (PARc) and radiation use efficiency (RUE) from emergence to flowering. A combination of architectural and light transfer models was used to calculate light interception in three field experiments with two cold-tolerant lines and at two sowing dates. Model assessment confirmed that the approach was suitable to infer light interception. Biomass production was strongly affected by early sowings. RUE was identified as the main cause of biomass reduction during cold events. Furthermore, PARc explained most of the variability observed at flowering, its relative contributions being more or less important according to the climate experienced. Cold temperatures resulted in lower PARc, mainly because final leaf length and width were significantly reduced for all leaves emerging after the first cold occurrence. These results confirm that virtual plants can be useful as fine phenotyping tools. A scheme of action of cold on leaf expansion, light interception and radiation use efficiency is discussed with a view towards helping breeders define relevant selection criteria. This paper originates from a presentation at the 5th International Workshop on Functional–Structural Plant Models, Napier, New Zealand, November 2007.

The accuracy of varietal selection using factor analytic models for multi-environment plant breeding trials

Modeling of cultivar x trial effects for multienvironment trials (METs) within a mixed model framework is now common practice in many plant breeding programs. The factor analytic (FA) model is a parsimonious form used to approximate the fully unstructured form of the genetic variance-covariance matrix in the model for MET data. In this study, we demonstrate that the FA model is generally the model of best fit across a range of data sets taken from early generation trials in a breeding program. In addition, we demonstrate the superiority of the FA model in achieving the most common aim of METs, namely the selection of superior genotypes. Selection is achieved using best linear unbiased predictions (BLUPs) of cultivar effects at each environment, considered either individually or as a weighted average across environments. In practice, empirical BLUPs (E-BLUPs) of cultivar effects must be used instead of BLUPs since variance parameters in the model must be estimated rather than assumed known. While the optimal properties of minimum mean squared error of prediction (MSEP) and maximum correlation between true and predicted effects possessed by BLUPs do not hold for E-BLUPs, a simulation study shows that E-BLUPs perform well in terms of MSEP.