Using Bayesian Mixture Models That Combine Expert Knowledge and GIS Data to Define Ecoregions

Conservation planning and management programs typically assume relatively homogeneous ecological landscapes. Such “ecoregions” serve multiple purposes: they support assessments of competing environmental values, reveal priorities for allocating scarce resources, and guide effective on-ground actions such as the acquisition of a protected area and habitat restoration. Ecoregions have evolved from a history of organism–environment interactions, and are delineated at the scale or level of detail required to support planning. Depending on the delineation method, scale, or purpose, they have been described as provinces, zones, systems, land units, classes, facets, domains, subregions, and ecological, biological, biogeographical, or environmental regions. In each case, they are essential to the development of conservation strategies and are embedded in government policies at multiple scales.

Mixture models of nucleotide sequence evolution that account for heterogeneity in the substitution process across sites and across lineages

Molecular phylogenetic studies of homologous sequences of nucleotides often assume that the underlying evolutionary process was globally stationary, reversible, and homogeneous (SRH), and that a model of evolution with one or more site-specific and time-reversible rate matrices (e.g., the GTR rate matrix) is enough to accurately model the evolution of data over the whole tree. However, an increasing body of data suggests that evolution under these conditions is an exception, rather than the norm. To address this issue, several non-SRH models of molecular evolution have been proposed, but they either ignore heterogeneity in the substitution process across sites (HAS) or assume it can be modeled accurately using the distribution. As an alternative to these models of evolution, we introduce a family of mixture models that approximate HAS without the assumption of an underlying predefined statistical distribution. This family of mixture models is combined with non-SRH models of evolution that account for heterogeneity in the substitution process across lineages (HAL). We also present two algorithms for searching model space and identifying an optimal model of evolution that is less likely to over- or underparameterize the data. The performance of the two new algorithms was evaluated using alignments of nucleotides with 10 000 sites simulated under complex non-SRH conditions on a 25-tipped tree. The algorithms were found to be very successful, identifying the correct HAL model with a 75% success rate (the average success rate for assigning rate matrices to the tree's 48 edges was 99.25%) and, for the correct HAL model, identifying the correct HAS model with a 98% success rate. Finally, parameter estimates obtained under the correct HAL-HAS model were found to be accurate and precise. The merits of our new algorithms were illustrated with an analysis of 42 337 second codon sites extracted from a concatenation of 106 alignments of orthologous genes encoded by the nuclear genomes of Saccharomyces cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, S. castellii, S. kluyveri, S. bayanus, and Candida albicans. Our results show that second codon sites in the ancestral genome of these species contained 49.1% invariable sites, 39.6% variable sites belonging to one rate category (V1), and 11.3% variable sites belonging to a second rate category (V2). The ancestral nucleotide content was found to differ markedly across these three sets of sites, and the evolutionary processes operating at the variable sites were found to be non-SRH and best modeled by a combination of eight edge-specific rate matrices (four for V1 and four for V2). The number of substitutions per site at the variable sites also differed markedly, with sites belonging to V1 evolving slower than those belonging to V2 along the lineages separating the seven species of Saccharomyces. Finally, sites belonging to V1 appeared to have ceased evolving along the lineages separating S. cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, and S. bayanus, implying that they might have become so selectively constrained that they could be considered invariable sites in these species.

Overfitting Bayesian mixture models with an unknown number of components

This paper proposes solutions to three issues pertaining to the estimation of finite mixture models with an unknown number of components: the non-identifiability induced by overfitting the number of components, the mixing limitations of standard Markov Chain Monte Carlo (MCMC) sampling techniques, and the related label switching problem. An overfitting approach is used to estimate the number of components in a finite mixture model via a Zmix algorithm. Zmix provides a bridge between multidimensional samplers and test based estimation methods, whereby priors are chosen to encourage extra groups to have weights approaching zero. MCMC sampling is made possible by the implementation of prior parallel tempering, an extension of parallel tempering. Zmix can accurately estimate the number of components, posterior parameter estimates and allocation probabilities given a sufficiently large sample size. The results will reflect uncertainty in the final model and will report the range of possible candidate models and their respective estimated probabilities from a single run. Label switching is resolved with a computationally light-weight method, Zswitch, developed for overfitted mixtures by exploiting the intuitiveness of allocation-based relabelling algorithms and the precision of label-invariant loss functions. Four simulation studies are included to illustrate Zmix and Zswitch, as well as three case studies from the literature. All methods are available as part of the R package Zmix, which can currently be applied to univariate Gaussian mixture models.

Distributions of lesser mealworm (Coleoptera: Tenebrionidae) in litter of a compacted earth floor broiler house in subtropical Queensland, Australia

Distributions of lesser mealworm, Alphitobius diaperinus (Panzer) (Coleoptera: Tenebrionidae), in litter of a compacted earth floor broiler house in southeastern Queensland, Australia, were studied over two flocks. Larvae were the predominant stage recorded. Significantly low densities occurred in open locations and under drinker cups where chickens had complete access, whereas high densities were found under feed pans and along house edges where chicken access was restricted. For each flock, lesser mealworm numbers increased at all locations over the first 14 d, especially under feed pans and along house edges, peaking at 26 d and then declining over the final 28 d. A life stage profile per flock was devised that consisted of the following: beetles emerge from the earth floor at the beginning of each flock, and females lay eggs, producing larvae that peak in numbers at 3 wk; after a further 3 to 4 wk, larvae leave litter to pupate in the earth floor, and beetles then emerge by the end of the flock time. Removing old litter from the brooder section at the end of a flock did not greatly reduce mealworm numbers over the subsequent flock, but it seemed to prevent numbers increasing, while an increase in numbers in the grow-out section was recorded after reusing litter. Areas under feed pans and along house edges accounted for 5% of the total house area, but approximately half the estimated total number of lesser mealworms in the broiler house occurred in these locations. The results of this study will be used to determine optimal deployment of site-specific treatments for lesser mealworm control.

Thermoanalysis of ammonium perchlorate and sodium nitrate mixture

A binary mixture of ammonium perchlorate-sodium nitrate in molar proportion undergoes partial fusion at 223°C and the transformation of the mixture to sodium perchlorate-ammonium nitrate occurs in the broad endothermic region. The mixture was heated and quenched at various temperatures in a differential thermal analysis assembly. Thermogravimetric analysis, X-ray diffraction, and infrared spectroscopic techniques were used to determine the composition of the quenched sample in order to explain the overall thermal phenomenon. Visual observations of the morphological changes that occur during the course of heating were made using a hot-stage microscope, 30–350°C.

Effect of feeding sorghum ergot (Claviceps africana) to sows during mid-lactation on plasma prolactin and litter performance

Diets containing 3% sorghum ergot (16 mg alkaloids/kg, including 14 mg dihydroergosine/kg) were fed to 12 sows from 14 days post-farrowing until weaning 14 days later, and their performance was compared with that of 10 control sows. Ergot-fed sows displayed a smaller weight loss during lactation of 24 kg/head vs. 29 kg/head in control sows (p > 0.05) despite feed consumption being less (61 kg/head total feed intake vs. 73 kg/head by control sows; p < 0.05). Ergot-fed sows had poorer weight gain of litters over the 14-day period (16.6 kg/litter vs. 28.3 kg/litter for controls; p < 0.05) despite an increase in consumption of creep feed by the piglets from the ergot-fed sows (1.9 kg/litter compared with 1.1 kg/litter by the control; p > 0.05). Sow plasma prolactin was reduced with ergot feeding after 7 days to 4.8 μg/l compared with 15.1 μg/l in the control sows (p < 0.01) and then at weaning was 4.9 μg/l compared with 8.0 μg/l (p < 0.01) in the control sows. Two sows fed ergot ceased lactation early, and the above sow feed intakes, body weight losses with litter weight gains and creep consumption indirectly indicate an ergot effect on milk production.

Modelling predicts positive and negative interactions between three Australian tropical tree species in monoculture and binary mixture

Computer modelling promises to be an important tool for analysing and predicting interactions between trees within mixed species forest plantations. This study explored the use of an individual-based mechanistic model as a predictive tool for designing mixed species plantations of Australian tropical trees. The `spatially explicit individually based-forest simulator' (SeXI-FS) modelling system was used to describe the spatial interaction of individual tree crowns within a binary mixed-species experiment. The three-dimensional model was developed and verified with field data from three forest tree species grown in tropical Australia. The model predicted the interactions within monocultures and binary mixtures of Flindersia brayleyana, Eucalyptus pellita and Elaeocarpus grandis, accounting for an average of 42% of the growth variation exhibited by species in different treatments. The model requires only structural dimensions and shade tolerance as species parameters. By modelling interactions in existing tree mixtures, the model predicted both increases and reductions in the growth of mixtures (up to +/-50% of stem volume at 7 years) compared to monocultures. This modelling approach may be useful for designing mixed tree plantations.

NOx Reduction and Desorption Studies Under Electric Discharge Plasma Using a Simulated Gas Mixture: A Case Study on the Effect of Corona Electrodes

In this study, reduction and desorption of oxides of nitrogen (NOx) were conducted using an electrical discharge plasma technique. The study was carried out using a simulated gas mixture to explore the possibility of re-generation of used adsorbents by a nonthermal plasma desorption technique. Three different types of corona electrodes, namely, pipe, helical wire, and straight wire, were used for analyzing their effectiveness in NOx reduction/desorption. The pipe-type corona electrode exhibited a nitric oxide (NO) conversion of 50%, which is 1.5 times that of the straight-wire-type electrode at an energy density of 175J/L. The helical-wire-type corona electrode exhibited a NOx desorption efficiency almost 4 times that of the pipe-type electrode,indicating the possibility that corona-generated species play a crucial role in desorption.

Comparative productivity of irrigated short-term ryegrass (Lolium multiflorum) pasture receiving nitrogen, grown alone or in a mixture with white (Trifolium repens) and Persian (T. resupinatum) clovers

Dairy farms in subtropical Australia use irrigated, annually sown short-term ryegrass (Lolium multiflorum) or mixtures of short-term ryegrass and white (Trifolium repens) and Persian (shaftal) (T. resupinatum) clover during the winter-spring period in all-year-round milk production systems. A series of small plot cutting experiments was conducted in 3 dairying regions (tropical upland, north Queensland, and subtropical southeast Queensland and northern New South Wales) to determine the most effective rate and frequency of application of nitrogen (N) fertiliser. The experiments were not grazed, nor was harvested material returned to the plots, after sampling. Rates up to 100 kg N/ha.month (as urea or calcium ammonium nitrate) and up to 200 kg N/ha every 2 months (as urea) were applied to pure stands of ryegrass in 1991. In 1993 and 1994, urea, at rates up to 150 kg N/ha.month and to 200 kg N/ha every 2 months, was applied to pure stands of ryegrass; urea, at rates up to 50 kg N/ha.month, was also applied to ryegrass-clover mixtures. The results indicate that applications of 50-85 kg N/ha.month can be recommended for short-term ryegrass pastures throughout the subtropics and tropical uplands of eastern Australia, irrespective of soil type. At this rate, dry matter yields will reach about 90% of their potential, forage nitrogen concentration will be increased, there is minimal risk to stock from nitrate poisoning and there will be no substantial increase in soil N. The rate of N for ryegrass-clover pastures is slightly higher than for pure ryegrass but, at these rates, the clover component will be suppressed. However, increased ryegrass yields and higher forage nitrogen concentrations will compensate for the reduced clover component. At application rates up to 100 kg N/ha.month, build-up of NO3--N and NH4+-N in soil was generally restricted to the surface layers (0-20 cm) of the soil, but there was a substantial increase throughout the soil profile at 150 kg N/ha.month. The build-up of NO3--N and NH4+-N was greater and was found at lower rates on the lighter soil compared with heavy clays. Generally, most of the soil N was in the NO3--N form and most was in the top 20 cm.

The performance of irrigated mixtures of tall fescue, ryegrass and white clover in subtropical Australia. 1. The effects of sowing mixture combinations, nitrogen and oversowing on establishment, productivity, botanical composition and persistence.

In the subtropics of Australia, the ryegrass component of irrigated perennial ryegrass (Lolium perenne) - white clover (Trifolium repens) pastures declines by approximately 40% in the summer following establishment, being replaced by summer-active C4 grasses. Tall fescue (Festuca arundinacea) is more persistent than perennial ryegrass and might resist this invasion, although tall fescue does not compete vigorously as a seedling. This series of experiments investigated the influence of ryegrass and tall fescue genotype, sowing time and sowing mixture as a means of improving tall fescue establishment and the productivity and persistence of tall fescue, ryegrass and white clover-based mixtures in a subtropical environment. Tall fescue frequency at the end of the establishment year decreased as the number of companion species sown in the mixture increased. Neither sowing mixture combinations nor sowing rates influenced overall pasture yield (of around 14 t/ha) in the establishment year but had a significant effect on botanical composition and component yields. Perennial ryegrass was less competitive than short-rotation ryegrass, increasing first-year yields of tall fescue by 40% in one experiment and by 10% in another but total yield was unaffected. The higher establishment-year yield (3.5 t/ha) allowed Dovey tall fescue to compete more successfully with the remaining pasture components than Vulcan (1.4 t/ha). Sowing 2 ryegrass cultivars in the mixture reduced tall fescue yields by 30% compared with a single ryegrass (1.6 t/ha), although tall fescue alone achieved higher yields (7.1 t/ha). Component sowing rate had little influence on composition or yield. Oversowing the ryegrass component into a 6-week-old sward of tall fescue and white clover improved tall fescue, white clover and overall yields in the establishment year by 83, 17 and 11%, respectively, but reduced ryegrass yields by 40%. The inclusion of red (T. pratense) and Persian (T. resupinatum) clovers and chicory (Cichorium intybus) increased first-year yields by 25% but suppressed perennial grass and clover components. Yields were generally maintained at around 12 t/ha/yr in the second and third years, with tall fescue becoming dominant in all 3 experiments. The lower tall fescue seeding rate used in the first experiment resulted in tall fescue dominance in the second year following establishment, whereas in Experiments 2 and 3 dominance occurred by the end of the first year. Invasion by the C4 grasses was relatively minor (<10%) even in the third year. As ryegrass plants died, tall fescue and, to a lesser extent, white clover increased as a proportion of the total sward. Treatment effects continued into the second, but rarely the third, year and mostly affected the yield of one of the components rather than total cumulative yield. Once tall fescue became dominant, it was difficult to re-introduce other pasture components, even following removal of foliage and moderate renovation. Severe renovation (reducing the tall fescue population by at least 30%) seems a possible option for redressing this situation.

Presence and incidence of food-borne pathogens in Australian chicken litter.

1. Litter samples were collected at the end of the production cycle from spread litter in a single shed from each of 28 farms distributed across the three Eastern seaboard States of Australia. 2. The geometric mean for Salmonella was 44 Most Probable Number (MPN)/g for the 20 positive samples. Five samples were between 100 and 1000 MPN/g and one at 105 MPN/g, indicating a range of factors are contributing to these varying loads of this organism in litter. 3. The geometric mean for Campylobacter was 30 MPN/g for the 10 positive samples, with 7 of these samples being 100 MPN/g. The low prevalence and incidence of Campylobacter were possibly due to the rapid die-off of this organism. 4. E. coli values were markedly higher than the two key pathogens (geometric mean 20 x 105 colony forming units (cfu)/g) with overall values being more or less within the same range across all samples in the trial, suggesting a uniform contribution pattern of these organisms in litter. 5. Listeria monocytogenes was absent in all samples and this organism appears not to be an issue in litter. 6. The dominant (70% of the isolates) Salmonella serovar was S. Sofia (a common serovar isolated from chickens in Australia) and was isolated across all regions. Other major serovars were S. Virchow and S. Chester (at 10%) and S. Bovismorbificans and S. Infantis (at 8%) with these serovars demonstrating a spatial distribution across the major regions tested. 7. There is potential to re-use litter in the environment depending on end use and the support of relevant application practices and guidelines.

Thermoanalysis of ammonium perchlorate and sodium nitrate mixture

A binary mixture of ammonium perchlorate-sodium nitrate in molar proportion undergoes partial fusion at 223°C and the transformation of the mixture to sodium perchlorate-ammonium nitrate occurs in the broad endothermic region. The mixture was heated and quenched at various temperatures in a differential thermal analysis assembly. Thermogravimetric analysis, X-ray diffraction, and infrared spectroscopic techniques were used to determine the composition of the quenched sample in order to explain the overall thermal phenomenon. Visual observations of the morphological changes that occur during the course of heating were made using a hot-stage microscope, 30–350°C.