Pesticide residues in food - acute dietary exposure

Consumer risk assessment is a crucial step in the regulatory approval of pesticide use on food crops. Recently, an additional hurdle has been added to the formal consumer risk assessment process with the introduction of short-term intake or exposure assessment and a comparable short-term toxicity reference, the acute reference dose. Exposure to residues during one meal or over one day is important for short-term or acute intake. Exposure in the short term can be substantially higher than average because the consumption of a food on a single occasion can be very large compared with typical long-term or mean consumption and the food may have a much larger residue than average. Furthermore, the residue level in a single unit of a fruit or vegetable may be higher by a factor (defined as the variability factor, which we have shown to be typically ×3 for the 97.5th percentile unit) than the average residue in the lot. Available marketplace data and supervised residue trial data are examined in an investigation of the variability of residues in units of fruit and vegetables. A method is described for estimating the 97.5th percentile value from sets of unit residue data. Variability appears to be generally independent of the pesticide, the crop, crop unit size and the residue level. The deposition of pesticide on the individual unit during application is probably the most significant factor. The diets used in the calculations ideally come from individual and household surveys with enough consumers of each specific food to determine large portion sizes. The diets should distinguish the different forms of a food consumed, eg canned, frozen or fresh, because the residue levels associated with the different forms may be quite different. Dietary intakes may be calculated by a deterministic method or a probabilistic method. In the deterministic method the intake is estimated with the assumptions of large portion consumption of a ‘high residue’ food (high residue in the sense that the pesticide was used at the highest recommended label rate, the crop was harvested at the smallest interval after treatment and the residue in the edible portion was the highest found in any of the supervised trials in line with these use conditions). The deterministic calculation also includes a variability factor for those foods consumed as units (eg apples, carrots) to allow for the elevated residue in some single units which may not be seen in composited samples. In the probabilistic method the distribution of dietary consumption and the distribution of possible residues are combined in repeated probabilistic calculations to yield a distribution of possible residue intakes. Additional information such as percentage commodity treated and combination of residues from multiple commodities may be incorporated into probabilistic calculations. The IUPAC Advisory Committee on Crop Protection Chemistry has made 11 recommendations relating to acute dietary exposure.

A new data source for fisheries resource assessment: genetic estimates of the effective number of spawners. Final Report to the Fisheries Research and Development Corporation.

The development of innovative methods of stock assessment is a priority for State and Commonwealth fisheries agencies. It is driven by the need to facilitate sustainable exploitation of naturally occurring fisheries resources for the current and future economic, social and environmental well being of Australia. This project was initiated in this context and took advantage of considerable recent achievements in genomics that are shaping our comprehension of the DNA of humans and animals. The basic idea behind this project was that genetic estimates of effective population size, which can be made from empirical measurements of genetic drift, were equivalent to estimates of the successful number of spawners that is an important parameter in process of fisheries stock assessment. The broad objectives of this study were to 1. Critically evaluate a variety of mathematical methods of calculating effective spawner numbers (Ne) by a. conducting comprehensive computer simulations, and by b. analysis of empirical data collected from the Moreton Bay population of tiger prawns (P. esculentus). 2. Lay the groundwork for the application of the technology in the northern prawn fishery (NPF). 3. Produce software for the calculation of Ne, and to make it widely available. The project pulled together a range of mathematical models for estimating current effective population size from diverse sources. Some of them had been recently implemented with the latest statistical methods (eg. Bayesian framework Berthier, Beaumont et al. 2002), while others had lower profiles (eg. Pudovkin, Zaykin et al. 1996; Rousset and Raymond 1995). Computer code and later software with a user-friendly interface (NeEstimator) was produced to implement the methods. This was used as a basis for simulation experiments to evaluate the performance of the methods with an individual-based model of a prawn population. Following the guidelines suggested by computer simulations, the tiger prawn population in Moreton Bay (south-east Queensland) was sampled for genetic analysis with eight microsatellite loci in three successive spring spawning seasons in 2001, 2002 and 2003. As predicted by the simulations, the estimates had non-infinite upper confidence limits, which is a major achievement for the application of the method to a naturally-occurring, short generation, highly fecund invertebrate species. The genetic estimate of the number of successful spawners was around 1000 individuals in two consecutive years. This contrasts with about 500,000 prawns participating in spawning. It is not possible to distinguish successful from non-successful spawners so we suggest a high level of protection for the entire spawning population. We interpret the difference in numbers between successful and non-successful spawners as a large variation in the number of offspring per family that survive – a large number of families have no surviving offspring, while a few have a large number. We explored various ways in which Ne can be useful in fisheries management. It can be a surrogate for spawning population size, assuming the ratio between Ne and spawning population size has been previously calculated for that species. Alternatively, it can be a surrogate for recruitment, again assuming that the ratio between Ne and recruitment has been previously determined. The number of species that can be analysed in this way, however, is likely to be small because of species-specific life history requirements that need to be satisfied for accuracy. The most universal approach would be to integrate Ne with spawning stock-recruitment models, so that these models are more accurate when applied to fisheries populations. A pathway to achieve this was established in this project, which we predict will significantly improve fisheries sustainability in the future. Regardless of the success of integrating Ne into spawning stock-recruitment models, Ne could be used as a fisheries monitoring tool. Declines in spawning stock size or increases in natural or harvest mortality would be reflected by a decline in Ne. This would be good for data-poor fisheries and provides fishery independent information, however, we suggest a species-by-species approach. Some species may be too numerous or experiencing too much migration for the method to work. During the project two important theoretical studies of the simultaneous estimation of effective population size and migration were published (Vitalis and Couvet 2001b; Wang and Whitlock 2003). These methods, combined with collection of preliminary genetic data from the tiger prawn population in southern Gulf of Carpentaria population and a computer simulation study that evaluated the effect of differing reproductive strategies on genetic estimates, suggest that this technology could make an important contribution to the stock assessment process in the northern prawn fishery (NPF). Advances in the genomics world are rapid and already a cheaper, more reliable substitute for microsatellite loci in this technology is available. Digital data from single nucleotide polymorphisms (SNPs) are likely to super cede ‘analogue’ microsatellite data, making it cheaper and easier to apply the method to species with large population sizes.

Odour Emissions from Tunnel Ventilated Poultry Housing

In Australia, factors such as local planning processes, urban encroachment into rural areas and intensification of the poultry industry have increased the potential for odour and dust nuisance. At present, accurate estimates of odour emissions from mechanically ventilated poultry housing systems do not exist for Australian conditions. This has made the poultry industry vulnerable to unsubstantiated criticism. Recently, the Australian poultry industry have made a significant investment in research to obtain accurate estimates of odour, dust and volatile chemical emission rates from typical poultry housing systems. This paper describes the measurement of odour emissions from tunnel ventilated poultry housing systems in different climatic zones in Queensland and Victoria, Australia (humid sub-tropical and Mediterranean respectively) during two seasons (summer and winter). Samples were collected at defined intervals over typical batch production cycles to define the odour emission profiles. These samples were analysed using dynamic olfactometry according to the Australian Standard 4323.3 to derive the odour concentration values. Ventilation rates were measured concurrently, allowing the calculation of odour emission rates. Odour concentration and emission rates were assessed in terms of ventilation rate, ambient and shed air temperature and relative humidity and litter moisture status. Odour emission rates varied with bird age. Seasonal differences in odour emission rate were also observed.

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.

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.

Fate of applied biosolids nitrogen in a cut and remove forage system on an alluvial clay loam soil

The fate of nitrogen (N) applied in biosolids was investigated in a forage production system on an alluvial clay loam soil in south-eastern Queensland, Australia. Biosolids were applied in October 2002 at rates of 6, 12, 36, and 54dryt/ha for aerobically digested biosolids (AE) and 8, 16, 48, and 72dryt/ha for anaerobically digested biosolids (AN). Rates were based on multiples of the Nitrogen Limited Biosolids Application rate (0.5, 1, 3, and 4.5NLBAR) for each type of biosolid. The experiment included an unfertilised control and a fertilised control that received multiple applications of synthetic fertiliser. Forage sorghum was planted 1 week after biosolids application and harvested 4 times between December 2002 and May 2003. Dry matter production was significantly greater from the biosolids-treated plots (21-27t/ha) than from the unfertilised (16t/ha) and fertilised (18t/ha) controls. The harvested plant material removed an extra 148-488kg N from the biosolids-treated plots. Partial N budgets were calculated for the 1NLBAR and 4.5NLBAR treatments for each biosolids type at the end of the crop season. Crop removal only accounted for 25-33% of the applied N in the 1NLBAR treatments and as low as 8-15% with 4.5NLBAR. Residual biosolids N was predominantly in the form of organic N (38-51% of applied biosolids N), although there was also a significant proportion (10-23%) as NO3-N, predominantly in the top 0.90m of the soil profile. From 12 to 29% of applied N was unaccounted for, and presumed to be lost as gaseous nitrogen and/or ammonia, as a consequence of volatilisation or denitrification, respectively. In-season mineralisation of organic N in biosolids was 43-59% of the applied organic N, which was much greater than the 15% (AN)-25% (AE) expected, based on current NLBAR calculation methods. Excessive biosolids application produced little additional biomass but led to high soil mineral N concentrations that were vulnerable to multiple loss pathways. Queensland Guidelines need to account for higher rates of mineralisation and losses via denitrification and volatilisation and should therefore encourage lower application rates to achieve optimal plant growth and minimise the potential for detrimental impacts on the environment.

A consensus genetic map of sorghum that integrates multiple component maps and high-throughput Diversity Array Technology (DArT) markers

Background: Sorghum genome mapping based on DNA markers began in the early 1990s and numerous genetic linkage maps of sorghum have been published in the last decade, based initially on RFLP markers with more recent maps including AFLPs and SSRs and very recently, Diversity Array Technology (DArT) markers. It is essential to integrate the rapidly growing body of genetic linkage data produced through DArT with the multiple genetic linkage maps for sorghum generated through other marker technologies. Here, we report on the colinearity of six independent sorghum component maps and on the integration of these component maps into a single reference resource that contains commonly utilized SSRs, AFLPs, and high-throughput DArT markers. Results: The six component maps were constructed using the MultiPoint software. The lengths of the resulting maps varied between 910 and 1528 cM. The order of the 498 markers that segregated in more than one population was highly consistent between the six individual mapping data sets. The framework consensus map was constructed using a "Neighbours" approach and contained 251 integrated bridge markers on the 10 sorghum chromosomes spanning 1355.4 cM with an average density of one marker every 5.4 cM, and were used for the projection of the remaining markers. In total, the sorghum consensus map consisted of a total of 1997 markers mapped to 2029 unique loci ( 1190 DArT loci and 839 other loci) spanning 1603.5 cM and with an average marker density of 1 marker/0.79 cM. In addition, 35 multicopy markers were identified. On average, each chromosome on the consensus map contained 203 markers of which 58.6% were DArT markers. Non-random patterns of DNA marker distribution were observed, with some clear marker-dense regions and some marker-rare regions. Conclusion: The final consensus map has allowed us to map a larger number of markers than possible in any individual map, to obtain a more complete coverage of the sorghum genome and to fill a number of gaps on individual maps. In addition to overall general consistency of marker order across individual component maps, good agreement in overall distances between common marker pairs across the component maps used in this study was determined, using a difference ratio calculation. The obtained consensus map can be used as a reference resource for genetic studies in different genetic backgrounds, in addition to providing a framework for transferring genetic information between different marker technologies and for integrating DArT markers with other genomic resources. DArT markers represent an affordable, high throughput marker system with great utility in molecular breeding programs, especially in crops such as sorghum where SNP arrays are not publicly available.

The use of heat-sums to predict harvest dates following potassium chlorate (KClO3) application in longan.

The ability to initiate and manipulate flowering with KClO3 allows flowering of longan, to be triggered outside of the normal flowering season (July-September) in Australia. Fruit maturity following normal flowering will occur approximately six-eight months (180-220 days) from flowering, depending on variety. Out of season flowering will result in differing times to maturity due to different temperature regimes during the maturity period. Knowing how long fruit will take to mature from different KClO3 application dates is potentially a valuable tool for growers to use as it would allow them to time their applications with market opportunities, e.g. Chinese New Year, periods of low volumes or periods of high prices. A simple heat-sum calculation was shown to reliably quantify fruit maturity periods, 2902 and 3432 growing degree days for Kohala and Biew Kiew respectively. Growers can use heat-sum as a predictive tool to allow for efficient planning of harvesting, packaging and freight requirements.

Lot feeding and intensive livestock manure and effluent management

Manure management emissions may present much greater opportunity for greenhouse gas mitigation in the feedlot, pig, chicken meat, egg and diary industries, than the current IPCC and DCC calculation guidelines suggest. Current literature and understanding of manure mass loss throughout the manure management system does not support these current guidelines; in which the emission rates are fixed and consequently don't allow incentives for reduced emissions.

Objective `Safe' Grazing Capacities for South-West Queensland Australia: Model Application and Evaluation

This paper describes the employment of two experienced graziers as consultants to apply and evaluate a model for calculating 'safe' long-term grazing capacities of individual properties. The model was based on ecological principles and entailed estimates of average annual forage grown (kglha) on the different land systems on each property and the calculation of the number of livestock (dry sheep equivalents, DSE) required to 'safely' utilise this forage. The grazier consultants applied and evaluated the 'safe' grazing capacity model on 20 properties of their choosing. For evaluation, model results were compared with; (a) the Department of Lands rated carrying capacities for those properties and (b) the grazing capacity assessed independently by the owners of those properties. For the 20 properties, the average 'safe' grazing capacity calculated by the model (21.0 DSE/kmZ) was 8% lighter than the average of the owner assessed capacities (22.7 DSE/kmZ), which in tum was 37% lighter than the average of the pre-1989 Department of Lands rated carrying capacity (31.0 DSE/kmZ). The grazing land management and administrative implications of these results and the role graziers played as consultants are discussed.

Linking spatial stock dynamics and economics: Evaluation of indicators and fishery management for the travelling eastern king prawn (Melicertus plebejus)

Reduced economic circumstances have moved management goals towards higher profit, rather than maximum sustainable yields in several Australian fisheries. The eastern king prawn is one such fishery, for which we have developed new methodology for stock dynamics, calculation of model-based and data-based reference points and management strategy evaluation. The fishery is notable for the northward movement of prawns in eastern Australian waters, from the State jurisdiction of New South Wales to that of Queensland, as they grow to spawning size, so that vessels fishing in the northern deeper waters harvest more large prawns. Bio-economic fishing data were standardized for calibrating a length-structured spatial operating model. Model simulations identified that reduced boat numbers and fishing effort could improve profitability while retaining viable fishing in each jurisdiction. Simulations also identified catch-rate levels that were effective for monitoring in simple within-year effort-control rules. However, favourable performance of catch-rate indicators was achieved only when a meaningful upper limit was placed on total allowed fishing effort. The methods and findings will allow improved measures for monitoring fisheries and inform decision makers on the uncertainty and assumptions affecting economic indicators.

Estimating an Aggregate Economic Genotype to Facilitate Breeding and Selection of Strawberry Genotypes in Southeast Queensland

The Queensland strawberry (Fragaria ×ananassa) breeding program in subtropical Australia aims to improve sustainable profitability for the producer. Selection must account for the relative economic importance of each trait and the genetic architecture underlying these traits in the breeding population. Our study used estimates of the influence of a trait on production costs and profitability to develop a profitability index (PI) and an economic weight (i.e., change in PI for a unit change in level of trait) for each trait. The economic weights were then combined with the breeding values for 12 plant and fruit traits on over 3000 genotypes that were represented in either the current breeding population or as progenitors in the pedigree of these individuals. The resulting linear combination (i.e., sum of economic weight × breeding value for all 12 traits) estimated the overall economic worth of each genotype as H, the aggregate economic genotype. H values were validated by comparisons among commercial cultivars and were also compared with the estimated gross margins. When the H value of ‘Festival’ was set as zero, the H values of genotypes in the pedigree ranged from –0.36 to +0.28. H was highly correlated (R2 = 0.77) with the year of selection (1945–98). The gross margins were highly linearly related (R2 > 0.98) to H values when the genotype was planted on less than 50% of available area, but the relationship was non-linear [quadratic with a maximum (R2 > 0.96)] when the planted area exceeded 50%. Additionally, with H values above zero, the variation in gross margin increased with increasing H values as the percentage of area planted to a genotype increased. High correlations among some traits allowed the omission of any one of three of the 12 traits with little or no effect on ranking (Spearman’s rank correlation 0.98 or greater). Thus, these traits may be dropped from the aggregate economic genotype, leading to either cost reductions in the breeding program or increased selection intensities for the same resources. H was efficient in identifying economically superior genotypes for breeding and deployment, but because of the non-linear relationship with gross margin, calculation of a gross margin for genotypes with high H is also necessary when cultivars are deployed across more than 50% of the available area.

Large-chamber methane and nitrous oxide measurements are comparable to the backward Lagrangian stochastic method

Measurement of individual emission sources (e.g., animals or pen manure) within intensive livestock enterprises is necessary to test emission calculation protocols and to identify targets for decreased emissions. In this study, a vented, fabric-covered large chamber (4.5 × 4.5 m, 1.5 m high; encompassing greater spatial variability than a smaller chamber) in combination with on-line analysis (nitrous oxide [N2O] and methane [CH4] via Fourier Transform Infrared Spectroscopy; 1 analysis min-1) was tested as a means to isolate and measure emissions from beef feedlot pen manure sources. An exponential model relating chamber concentrations to ambient gas concentrations, air exchange (e.g., due to poor sealing with the surface; model linear when ≈ 0 m3 s-1), and chamber dimensions allowed data to be fitted with high confidence. Alternating manure source emission measurements using the large-chamber and the backward Lagrangian stochastic (bLS) technique (5-mo period; bLS validated via tracer gas release, recovery 94-104%) produced comparable N2O and CH4 emission values (no significant difference at P < 0.05). Greater precision of individual measurements was achieved via the large chamber than for the bLS (mean ± standard error of variance components: bLS half-hour measurements, 99.5 ± 325 mg CH4 s-1 and 9.26 ± 20.6 mg N2O s-1; large-chamber measurements, 99.6 ± 64.2 mg CH4 s-1 and 8.18 ± 0.3 mg N2O s-1). The large-chamber design is suitable for measurement of emissions from manure on pen surfaces, isolating these emissions from surrounding emission sources, including enteric emissions. © © American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America.

Can soil nematode community structure be used to indicate soil carbon dynamics in horticultural systems?

There is a need to develop indicators that relate the dynamics of soil organic carbon (SOC) with changes in land management of horticultural production systems. Soil nematode communities have been shown to be sensitive to land management changes, but often do not include plant-parasites in the calculation of soil nematode community indices. The concept of nematode functional guilds was used to estimate the proportion of carbon entering the soil ecosystem through different channels, such as through decomposition of organic material, the detrital channel, through the roots of plants, the root channel or recycled through the activity of predators, a predation channel. Calculations of the indices were developed and validated using case studies in the north Queensland banana industry. Firstly, a survey of organic and conventional banana farms found a greater proportion of C entering the soil ecosystem through the detrital channel and a reduced proportion of C originating through the root channel at the organic sites relative to conventional sites. Secondly, a field experiment comparing compost amendments, found application of fresh compost significantly increased the proportion of C entering the soil ecosystem through the detrital channel and decreased proportion of C originating from the root channel. Thirdly, a field experiment comparing 'conventional' banana production to an 'alternative' system which incorporated organic matter, found the proportion of C entering the soil ecosystem through the root channel was significantly greater in the conventional systems relative to the alternative system. This research demonstrates that nematode indices can be used to assess horticultural systems, by indicating the origins of SOC.

Baseline and greenhouse-gas emissions in extensive livestock enterprises, with a case study of feeding lipid to beef cattle

For accurate calculation of reductions in greenhouse-gas (GHG) emissions, methodologies under the Australian Government's Carbon Farming Initiative (CFI) depend on a valid assessment of the baseline and project emissions. Life-cycle assessments (LCAs) clearly show that enteric methane emitted from the rumen of cattle and sheep is the major source of GHG emissions from livestock enterprises. Where a historic baseline for a CFI methodology for livestock is required, the use of simulated data for cow-calf enterprises at six sites in southern Australia demonstrated that a 5-year rolling emission average will provide an acceptable trade off in terms of accuracy and stability, but this is a much shorter time period than typically used for LCA. For many CFI livestock methodologies, comparative or pair-wise baselines are potentially more appropriate than historic baselines. A case study of lipid supplementation of beef cows over winter is presented. The case study of a control herd of 250 cows used a comparative baseline derived from simple data on livestock numbers and class of livestock to quantify the emission abatement. Compared with the control herd, lipid supplementation to cows over winter increased livestock productivity, total livestock production and enterprise GHG emissions from 990 t CO2-e to 1022 t CO2-e. Energy embodied in the supplement and extra diesel used in transporting the supplement diminished the enteric-methane abatement benefit of lipid supplementation. Reducing the cow herd to 238 cows maintained the level of livestock production of the control herd and reduced enterprise emissions to 938 t CO2-e, but was not cost effective under the assumptions of this case study.