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.

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.

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.

Effects of fishing on tropical reef associated shark populations on the Great Barrier Reef

Three data sets were examined to define the level of interaction of reef associated sharks with the commercial Coral Reef Fin Fish Fishery within the Great Barrier Reef (GBR). Data were examined from fishery logbooks, an observer program within the fishery and a fishery-independent survey conducted as part of the Effects of Line Fishing (ELF) Experiment. The majority of the identified catch was comprised of grey reef (62-72%), whitetip reef (16-29%) and blacktip reef (6-13%) sharks. Logbook data revealed spatially and temporally variable landings of shark from the GBR. Catch per unit effort (CPUE) through time was stable for the period from 1989 to 2006 with no evidence of increase or decline. Data from observer and ELF data sets indicated no differences in CPUE among regions. The ELF data set demonstrated that CPUE was higher in Marine National Park zones (no fishing) when compared to General Use zones (open to fishing). The ongoing and consistent catches of reef sharks in the fishery and effectiveness of no-fishing zones suggest that management zones within the GBR Marine Park are effective at protecting a portion of the reef shark population from exploitation.

Refined Global Analysis of Bemisia tabaci (Hemiptera: Sternorrhyncha: Aleyrodoidea: Aleyrodidae) Mitochondrial Cytochrome Oxidase 1 to Identify Species Level Genetic Boundaries.

Identifying species boundaries within morphologically indistinguishable cryptic species complexes is often contentious. For the whitefly Bemisia tabaci (Gennadius) (Hemiptera: Sternorrhyncha: Aleyrodoidea: Aleyrodidae), the lack of a clear understanding about the genetic limits of the numerous genetic groups and biotypes so far identified has resulted in a lack of consistency in the application of the terms, the approaches use to apply them and in our understanding of what genetic structure within B. tabaci means. Our response has been to use mitochondrial gene cytochrome oxidase one to consider how to clearly and consistently define genetic separation. Using Bayesian phylogenetic analysis and analysis of sequence pairwise divergence we found a considerably higher to number of genetic groups than had been previously determined with two breaks in the distribution, one at 11% and another at 3.5%. At >11% divergence, 11 distinct groups were resolved, whereas at >3.5% divergence 24 groups were identified. Consensus sequences for each of these groups were determined and were shown to be useful in the correct assignment of sequences of unknown origin. The 3.5% divergence bound is consistent with species level separations in other insect taxa and Suggests that B. tabaci is it cryptic species composed of at least 24 distinct species. We further show that the placement of Bemesia atriplex (Froggatt) within the B. tabaci in, group adds further weight to the argument for species level separation within B. tabaci. This new analysis, which constructs consensus sequences and uses these its a standard against which unknown sequences call be compared, provides for the first time it consistent means of identifying the genetic hounds of each species with it high degree of certainty.

MOE and MOR assessment technologies for improving graded recovery of exotic pines in Australia.

This project was designed to provide the structural softwood processing industry with the basis for improved green and dry grading to allow maximise MGP grade yields, consistent product performance and reduced processing costs. To achieve this, advanced statistical techniques were used in conjunction with state-of-the-art property measurement systems. Specifically, the project aimed to make two significant steps forward for the Australian structural softwood industry: • assessment of technologies, both existing and novel, that may lead to selection of a consistent, reliable and accurate device for the log yard and green mill. The purpose is to more accurately identify and reject material that will not make a minimum grade of MGP10 downstream; • improved correlation of grading MOE and MOR parameters in the dry mill using new analytical methods and a combination of devices. The three populations tested were stiffness-limited radiata pine, strength-limited radiata pine and Caribbean pine. Resonance tests were conducted on logs prior to sawmilling, and on boards. Raw data from existing in-line systems were captured for the green and dry boards. The dataset was analysed using classical and advanced statistical tools to provide correlations between data sets and to develop efficient strength and stiffness prediction equations. Stiffness and strength prediction algorithms were developed from raw and combined parameters. Parameters were analysed for comparison of prediction capabilities using in-line parameters, off-line parameters and a combination of in-line and off-line parameters. The results show that acoustic resonance techniques have potential for log assessment, to sort for low stiffness and/or low strength, depending on the resource. From the log measurements, a strong correlation was found between the average static MOE of the dried boards within a log and the predicted value. These results have application in segregating logs into structural and non-structural uses. Some commercial technologies are already available for this application such as Hitman LG640. For green boards it was found that in-line and laboratory acoustic devices can provide a good prediction of dry static MOE and moderate prediction for MOR.There is high potential for segregating boards at this stage of processing. Grading after the log breakdown can improve significantly the effectiveness of the mill. Subsequently, reductions in non-structural volumes can be achieved. Depending on the resource it can be expected that a 5 to 8 % reduction in non structural boards won’t be dried with an associated saving of $70 to 85/m3. For dry boards, vibration and a standard Metriguard CLT/HCLT provided a similar level of prediction on stiffness limited resource. However, Metriguard provides a better strength prediction in strength limited resources (due to this equipment’s ability to measure local characteristics). The combination of grading equipment specifically for stiffness related predictors (Metriguard or vibration) with defect detection systems (optical or X-ray scanner) provides a higher level of prediction, especially for MOR. Several commercial technologies are already available for acoustic grading on board such those from Microtec, Luxscan, Falcon engineering or Dynalyse AB for example. Differing combinations of equipment, and their strategic location within the processing chain, can dramatically improve the efficiency of the mill, the level of which will vary depending of the resource. For example, an initial acoustic sorting on green boards combined with an optical scanner associated with an acoustic system for grading dry board can result in a large reduction of the proportion of low value low non-structural produced. The application of classical MLR on several predictors proved to be effective, in particular for MOR predictions. However, the usage of a modern statistics approach(chemometrics tools) such as PLS proved to be more efficient for improving the level of prediction. Compared to existing technologies, the results of the project indicate a good improvement potential for grading in the green mill, ahead of kiln drying and subsequent cost-adding processes. The next stage is the development and refinement of systems for this purpose.

AG15 (Sustainable Agricultural State-level Investment Program) natural resource management extension in irrigated cotton and grain

Natural Resource Management project developing reources and supporting best practice management for irrigated cotton and grain growers in Queensland.

High-level coordination and strategy in the 2007 equine influenza outbreak response.

The equine influenza (EI) outbreak presented many challenges that required high-level coordination and decision making, as well as the development of new approaches for satisfactory and consistent resolution. This paper outlines the elements of the national coordination arrangements, preparatory arrangements in place prior to the outbreak that facilitated national coordination, and some of the issues faced and resolved in the response.

Animal-Level Factors Affecting Ovarian Function in Bos indicus Heifers Treated to Synchronize Ovulation with Intravaginal Progesterone-Releasing Devices and Oestradiol Benzoate

The primary objective of this study was to investigate the impact of animal-level factors including energy balance and environmental/management stress, on the ovarian function of Bos indicus heifers treated to synchronize ovulation. Two-year-old Brahman (BN) (n = 30) and BN-cross (n = 34) heifers were randomly allocated to three intravaginal progesterone-releasing device (IPRD) treatment groups: (i) standard-dose IPRD [Cue-Mate (R) (CM) 1.56 g; n = 17]; (ii) half-dose IPRD [0.78 g progesterone (P4); CM 0.78 g; n = 15]; (iii) half-dose IPRD + 300 IU equine chorionic gonadotrophin at IPRD removal (CM 0.78 g + G; n = 14); (iv) and a control group, 2x PGF2a [500 mu g prostaglandin F2a (PGF2a)] on Day -16 and -2 (n = 18). Intravaginal progesterone-releasing device-treated heifers received 250 mu g PGF2a at IPRD insertion (Day -10) and IPRD removal (Day -2) and 1 mg oestradiol benzoate on Day -10 and -1. Heifers were managed in a small feedlot and fed a defined ration. Ovarian function was evaluated by ultrasonography and plasma P4 throughout the synchronized and return cycles. Energy balance was evaluated using plasma insulin-like growth factor 1 (IGF-I) and glucose concentrations. The impact of environmental stressors was evaluated using plasma cortisol concentration. Heifers that had normal ovarian function had significantly higher IGF-I concentrations at commencement of the experiment (p = 0.008) and significantly higher plasma glucose concentrations at Day -2 (p = 0.040) and Day 4 (p = 0.043), than heifers with abnormal ovarian function. There was no difference between the mean pre-ovulatory cortisol concentrations of heifers that ovulated or did not ovulate. However, heifers that ovulated had higher cortisol concentrations at Day 4 (p = 0.056) and 6 (p = 0.026) after ovulation than heifers that did not ovulate.

The rph2 Gene Is Responsible for High Level Resistance to Phosphine in Independent Field Strains of Rhyzopertha dominica

The lesser grain borer Rhyzopertha dominica (F.) is one of the most destructive insect pests of stored grain. This pest has been controlled successfully by fumigation with phosphine for the last several decades, though strong resistance to phosphine in many countries has raised concern about the long term usefulness of this control method. Previous genetic analysis of strongly resistant (SR) R. dominica from three widely geographically dispersed regions of Australia, Queensland (SRQLD), New South Wales (SRNSW) and South Australia (SRSA), revealed a resistance allele in the rph1 gene in all three strains. The present study confirms that the rph1 gene contributes to resistance in a fourth strongly resistant strain, SR2(QLD), also from Queensland. The previously described rph2 gene, which interacts synergistically with rph1 gene, confers strong resistance on SRQLD and SRNSW. We now provide strong circumstantial evidence that weak alleles of rph2, together with rph1, contribute to the strong resistance phenotypes of SRSA and SR2(QLD). To test the notion that rph1 and rph2 are solely responsible for the strong resistance phenotype of all resistant R. dominica, we created a strain derived by hybridising the four strongly resistant lines. Following repeated selection for survival at extreme rates of phosphine exposure, we found only slightly enhanced resistance. This suggests that a single sequence of genetic changes was responsible for the development of resistance in these insects.

Modelling interactions between farm-level structural adjustment and a regional economy: A case of the Australian rice industry

Climate change and on-going water policy reforms will likely contribute to on-farm and regional structural adjustment in Australia. This paper gathers empirical evidence of farm-level structural adjustments and integrates these with a regional equilibrium model to investigate sectoral and regional impacts of climate change and recent water use policy on rice industry. We find strong evidence of adjustments to the farming system, enabled by existing diversity in on-farm production. A further loss of water with additional pressures to adopt less intensive and larger-scale farming, will however reduce the net number of farm businesses, which may affect regional rice production. The results from a regional CGE model show impacts on the regional economy over and above the direct cost of the environmental water, although a net reduction in real economic output and real income is partially offset by gains in rest of the Australia through the reallocation or resources. There is some interest within the industry and from potential new corporate entrants in the relocation of some rice production to the north. However, strong government support would be crucial to implement such relocation.

Impacts of level of utilisation by grazing on an Astrebla (Mitchell grass) grassland in north-western Queensland between 1984 and 2010. 1. Herbage mass and population dynamics of Astrebla spp

Managing large variations in herbage production, resulting from highly variable seasonal rainfall, provides a major challenge for the sustainable management of Astrebla (Mitchell grass) grasslands in Australia. A grazing study with sheep was conducted between 1984 and 2010 on an Astrebla grassland in northern Queensland to describe the effects of a range of levels of utilisation of the herbage at the end of the summer growing season (April–May in northern Australia) on the sustainability of these grasslands. In unreplicated paddocks, sheep numbers were adjusted annually to achieve 0, 10, 20, 30, 50 and 80% utilisation of the herbage mass at the end of the summer over the ensuing 12 months. Higher levels of utilisation reduced both total and Astrebla spp. herbage mass because of the effects of higher utilisation on Astrebla spp. and this effect was accentuated by drought. The tussock density of Astrebla spp. varied widely among years but with few treatment differences until 2005 when density was reduced at the 50% level of utilisation. A major change in density resulted from a large recruitment of Astrebla spp. in 1989 that influenced its density for the remainder of the study. Basal area of the tussocks fluctuated among years, with increases due to rainfall and decreases during droughts. Seasonal rainfall was more influential than level of utilisation in changes to the basal area of perennial grasses. Drought resulted in the death of Astrebla spp. tussocks and this effect was accentuated at higher levels of utilisation. A series of three grazing exclosures were used to examine the recovery of the density and basal area of Astrebla spp. after it had been reduced by 80% utilisation over the preceding 9 years. This recovery study indicated that, although grazing exclusion was useful in the recovery of Astrebla spp., above-average rainfall was the major factor driving increases in the basal area of perennial grasses. Spring values of the Southern Oscillation Index and associated rainfall probabilities were considered to have potential for understanding the dynamics of Astrebla spp. It was concluded that Astrebla grassland remained sustainable after 26 years when grazed at up to 30% utilisation, while, at 50% utilisation, they became unsustainable after 20 years. Results from this study emphasised the need to maintain the population of Astrebla spp. tussocks.

Impacts of level of utilisation by grazing on an Astrebla (Mitchell grass) grassland in north-western Queensland between 1984 and 2010. 2. Plant species richness and abundance

The occurrence of interstitial species in Astrebla grasslands in Australia are influenced by grazing and seasonal rainfall but the interactions of these two influences are complex. This paper describes three studies aimed at determining and explaining the changes in plant species richness and abundance of the interstitial species in a long-term sheep utilisation experiment in an Astrebla grassland in northern Queensland. In the first study, increasing utilisation increased the frequency of Dactyloctenium radulans (Button grass) and Brachyachne convergens (Downs couch) and reduced that of Streptoglossa adscendens (Mint bush). In the second study, seasonal rainfall variation between 1984 and 2009 resulted in large annual differences in the size of the seed banks of many species, but increasing utilisation consistently reduced the seed bank of species such as Astrebla spp. and S. adscendens and increased that of species such as B. convergens, D. radulans, Amaranthus mitchellii (Boggabri) and Boerhavia sp. (Tar vine). In the third study, the highest species richness occurred at the lightest utilisation because of the presence of a range of palatable forbs, especially legumes. Species richness was reduced as utilisation increased. Species richness in the grazing exclosure was low and similar to that at the heaviest utilisation where there was a reduction in the presence of palatable forb species. The pattern of highest species richness at the lightest grazing treatment was maintained across three sampling times, even with different amounts of seasonal rainfall, but there was a large yearly variation in both the density and frequency of many species. It was concluded that the maintenance of highest species richness at the lightest utilisation was not aligned with other data from this grazing experiment which indicated that the maximum sustainable wool production occurred at moderate utilisation.

Complementary forages – integration at a whole-farm level

A high proportion of the Australian and New Zealand dairy industry is based on a relatively simple, low input and low cost pasture feedbase. These factors enable this type of production system to remain internationally competitive. However, a key limitation of pasture-based dairy systems is periodic imbalances between herd intake requirements and pasture DM production, caused by strong seasonality and high inter-annual variation in feed supply. This disparity can be moderated to a certain degree through the strategic management of the herd through altering calving dates and stocking rates, and the feedbase by conserving excess forage and irrigating to flatten seasonal forage availability. Australasian dairy systems are experiencing emerging market and environmental challenges, which includes increased competition for land and water resources, decreasing terms of trade, a changing and variable climate, an increasing environmental focus that requires improved nutrient and water-use efficiency and lower greenhouse gas emissions. The integration of complementary forages has long been viewed as a means to manipulate the home-grown feed supply, to improve the nutritive value and DM intake of the diet, and to increase the efficiency of inputs utilised. Only recently has integrating complementary forages at the whole-farm system level received the significant attention and investment required to examine their potential benefit. Recent whole-of-farm research undertaken in both Australia and New Zealand has highlighted the importance of understanding the challenges of the current feedbase and the level of complementarity between forage types required to improve profit, manage risk and/or alleviate/mitigate against adverse outcomes. This paper reviews the most recent systems-level research into complementary forages, discusses approaches to modelling their integration at the whole-farm level and highlights the potential of complementary forages to address the major challenges currently facing pasture-based dairy systems.

A participatory systems approach to understanding climate adaptation needs

Emerging literature on climate adaptation suggests the need for effective ways of engaging or activating communities and supporting community roles, coupled with whole-of-system approaches to understanding climate change and adaptation needs. We have developed and evaluated a participatory approach to elicit community and stakeholder understanding of climate change adaptation needs, and connect diverse community members and local office bearers towards potential action. The approach was trialed in a series of connected social-ecological systems along a transect from a rural area to the coast and islands of ecologically sensitive Moreton Bay in Queensland, Australia. We conducted ‘climate roundtables’ in each of three areas along the transect, then a fourth roundtable reviewed and extended the results to the region as a whole. Influence diagrams produced through the process show how each climate variable forecast to affect this region (heat, storm, flood, sea-level rise, fire, drought) affects the natural environment, infrastructure, economic and social behaviour patterns, and psychosocial responses, and how sets of people, species and ecosystems are affected, and act, differentially. The participatory process proved effective as a way of building local empathy, a local knowledge base and empowering participants to join towards future climate adaptation action. Key principles are highlighted to assist in adapting the process for use elsewhere.