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.

Unmodified silver nanoparticles for rapid analysis of the organophosphorus pesticide, dipterex, often found in different waters

A novel, uncomplicated and rapid method of analysis for organophosphorus (OP) pesticides was researched and developed using the important, common OP, dipterex, as a typical example. The basis of the method involved the citrate-capped silver nanoparticles (citrate-capped AgNPs) and Acetylthiocholine (ATCh). The latter compound can be catalyzed by Acetylcholinesterase (AChE) to form thiocholine (TCh), which induces the aggregation of AgNPs. Correspondingly, the color of AgNPs in solution changes from bright yellow to pink, and the UV–vis characteristic absorption peak of AgNPs at about 400 nm decreases; simultaneously, a new absorption band appears at about 520 nm. Irreversible inhibition of AChE activity caused by dipterex, prevents the aggregation of AgNPs. Thus, a UV–vis spectrophotometric method was developed for the analysis of dipterex. The absorbance ratio A396 nm/A520 nm was found to be linearly related to the concentration of dipterex in the range of 0.25–37.5 ng mL−1 with a detection limit of 0.18 ng mL−1. This method was used successfully to analyse dipterex in spiked, different water samples.

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.

Temporal trend of pesticide concentrations in the Chikugo River (Japan) with changes in environmental regulation and field infrastructure

Contamination of pesticides, which are applied to rice paddy fields, in river water has been a major problem in Japan for decades. A prolonged water holding period after pesticide application in paddy fields is expected to reduce the concentration of rice pesticides in river water. Therefore, a long monitoring campaign was conducted from 2004 to 2010 to measure the concentrations of pesticides in water samples collected from several points along the Chikugo River (Japan) including tributaries and the main stream to see if there was any reduction in the level of pesticide contamination after the extension of the water holding period (from 3–4 days to 7 days) was introduced in 2007 by the new water management regulation. No significant difference (p > 0.05) was found in pesticide concentrations between the periods before and after 2007 in all monitoring points, except in one tributary where the pesticide concentrations after 2007 were even higher than that of the previous period. A detailed study in one of the tributaries also revealed that the renovated infrastructure did not reduce the pesticide concentrations in the drainage canals. Neither the introduction of the new regulation nor the improved infrastructure had any significant effect on reducing the contamination of pesticides in water of the Chikugo River. It is probably because most farmers did not properly implement the new requirement of holding paddy water within the field for 7 days after the application of pesticides. Only tightening the regulation would not be sufficient and more actions should be taken to enforce/provide extension support for the new water management regulation in order to reduce the level of residual pesticides in river water in Japan.

Effect of water management practice on pesticide behavior in paddy water

The fate and transport of three herbicides commonly used in rice production in Japan were compared using two water management practices. The herbicides were simetryn, thiobencarb and mefenacet. The first management practice was an intermittent irrigation scheme using an automatic irrigation system (AI) with a high drainage gate and the second one was a continuous irrigation and overflow drainage scheme (CI) in experimental paddy fields. Dissipation of the herbicides appeared to follow first order kinetics with the half-lives (DT₅₀) of 1.6-3.4 days and the DT₉₀ (90% dissipation) of 7.4-9.8 days. The AI scheme had little drainage even during large rainfall events thus resulting in losses of less than 4% of each applied herbicide through runoff. Meanwhile the CI scheme resulted in losses of about 37%, 12% and 35% of the applied masses of simetryn, thiobencarb and mefenacet, respectively. The intermittent irrigation scheme using an automatic irrigation system with a high drainage gate saved irrigation water and prevented herbicide runoff whereas the continuous irrigation and overflow scheme resulted in significant losses of water as well as the herbicides. Maintaining the excess water storage is important for preventing paddy water runoff during significant rainfall events. The organic carbon partition coefficient K_oc seems to be a strong indicator of the aquatic fate of the herbicide as compared to the water solubility (S_W). However, further investigations are required to understand the relation between K_oc and the agricultural practices upon the pesticide fate and transport. An extension of the water holding period up to 10 days after herbicide application based on the DT₉₀ from the currently specified period of 3-4 days in Japan is recommended to be a good agricultural practice for controlling the herbicide runoff from paddy fields. Also, the best water management practice, which can be recommended for use during the water holding period, is the intermittent irrigation scheme using an automatic irrigation system with a high drainage gate. © 2006 Elsevier B.V. All rights reserved.

Applicability of ELISA in pesticide monitoring to control runoff of bensulfuron-methyl and simetryn from paddy fields

The applicability of ELISA kits was evaluated as an alternative to monitor bensulfuron-methyl and simetryn behavior in paddy water under intermittent (Plot 1) and continuous (Plot 2) irrigation schemes. Simetryn concentrations in both plots decreased exponentially from the peak of the first day. However, the simetryn kit systematically underestimated by a factor of 0.79 as compared to the GC method. Bensulfuron-methyl concentrations exhibited similar dissipation kinetics in paddy water and the drainage water. The bensulfuron-methyl kit was capable of distinguishing spatial variations of concentrations in the paddy field. The ELISA kits clearly indicated differences in the loss of both herbicides between the two plots and therefore may be useful for evaluating the water management practice of pesticide runoff control in paddy fields.

Simulation of pesticide behavior in a paddy block by a pesticide fate and transport model

A simulation model (PCPF-B) was developed based on the PCPF-1 model to predict the runoff of pesticides from paddy plots to a drainage canal in a paddy block. The block-scale model now comprises three modules: (1) a module for pesticide application, (2) a module for pesticide behavior in paddy fields, and (3) a module for pesticide concentration in the drainage canal. The PCPF-B model was first evaluated by published data in a single plot and then was applied to predict the concentration of bensulfuron-methyl in one paddy block in the Sakura river basin, Ibaraki, Japan, where a detailed field survey was conducted. The PCPF-B model simulated well the behavior of bensulfuron-methyl in individual paddy plots. It also reflected the runoff pattern of bensulfuron-methyl at the block outlet, although overestimation of bensulfuronmethyl concentrations occurred due to uncertainty in water balance estimation. Application of water management practice such as water-holding period and seepage control also affected the performance of the model. A probabilistic approach may be necessary for a comprehensive risk assessment in large-scale paddy areas.

Pesticide discharge and water management in a paddy catchment in Japan

Concentrations of several pesticides were monitored in a paddy block and in the Kose river, which drains a paddy catchment in Fukuoka prefecture, Japan. Detailed water management in the block was also monitored to evaluate its effect on the pesticide contamination. The concentrations of applied pesticides in both block irrigation channel and drainage canal increased to tens of μg/L shortly after their applications. The increase in pesticide concentrations was well correlated with the open of irrigation and drainage gates in the pesticide-applied paddy plots only 1–3 days after pesticide application. High concentration of other pesticides, mainly herbicides, was also observed in the inflow irrigation and drainage waters, confirming the popularity of early irrigation and drainage after pesticide application in the area. The requirement of holding water after pesticide application (as a best management practice) issued by the authority was thus not properly followed. In a larger scale of the paddy catchment, the concentration of pesticides also increased significantly to several μg/L in the water of the Kose river shortly after the start of the pesticide application period either in downstream or mid–upstream areas, confirming the effect of current water management to the water quality. More extension and enforcement on water management should be done in order to control pesticide pollution from rice cultivation in Japan.

Exposure risk assessment and evaluation of the best management practice for controlling pesticide runoff from paddy fields. Part 2: Model simulation for the herbicide pretilachlor

BACKGROUND: Monitoring studies revealed high concentrations of pesticides in the drainage canal of paddy fields. It is important to have a way to predict these concentrations in different management scenarios as an assessment tool. A simulation model for predicting the pesticide concentration in a paddy block (PCPF-B) was evaluated and then used to assess the effect of water management practices for controlling pesticide runoff from paddy fields. RESULTS: The PCPF-B model achieved an acceptable performance. The model was applied to a constrained probabilistic approach using the Monte Carlo technique to evaluate the best management practices for reducing runoff of pretilachlor into the canal. The probabilistic model predictions using actual data of pesticide use and hydrological data in the canal showed that the water holding period (WHP) and the excess water storage depth (EWSD) effectively reduced the loss and concentration of pretilachlor from paddy fields to the drainage canal. The WHP also reduced the timespan of pesticide exposure in the drainage canal. CONCLUSIONS: It is recommended that: (1) the WHP be applied for as long as possible, but for at least 7 days, depending on the pesticide and field conditions; (2) an EWSD greater than 2 cm be maintained to store substantial rainfall in order to prevent paddy runoff, especially during the WHP.

Pesticide exposure assessment in rice paddy areas: A Japanese perspective

This chapter provides an overview of the Japanese regulatory issues regarding pesticide use in rice paddies and an introduction of the new pesticide registration program. In addition, the experience of the environmental monitoring of pesticides and the modeling approaches used for the calculation of predicted environmental concentrations (PECs) in surface water and ground water systems adjacent to rice paddies in Japan are also discussed. Japan has been one of the major pesticide users in the world. Although having a long history in rice cultivation, the pesticide exposure assessment for paddy rice production received less attention compared with EU and US. Applications of up-to-date techniques and the development of realistic assessment procedures under specific climatic conditions as well as mitigation management practices for controlling pesticide contamination are important for an environmental-friendly rice production. Through the international cooperation and research exchanges, advances in pesticide risk assessment for rice paddies in Asian region and other rice-growing areas in the world would contribute to sustainable rice production. Transplanting of rice seedlings grows almost all rice paddies in Japan. The land preparation starts around April and June, and the harvest season lasts from August to October depending on the region and the rice varieties. Most of the rice paddies are treated with herbicides and other crop protection products, such as fungicides and insecticides that are applied during the crop season accordingly. Newly developed insecticides and fungicides are also applied during seedbed preparation.

Modeling approaches for pesticide exposure assessment in rice paddies

Pesticide use in paddy rice production may contribute to adverse ecological effects in surface waters. Risk assessments conducted for regulatory purposes depend on the use of simulation models to determine predicted environment concentrations (PEC) of pesticides. Often tiered approaches are used, in which assessments at lower tiers are based on relatively simple models with conservative scenarios, while those at higher tiers have more realistic representations of physical and biochemical processes. This chapter reviews models commonly used for predicting the environmental fate of pesticides in rice paddies. Theoretical considerations, unique features, and applications are discussed. This review is expected to provide information to guide model selection for pesticide registration, regulation, and mitigation in rice production areas.

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.

Bayesian mixture model estimation of aerosol particle size distributions

In this paper, we examine approaches to estimate a Bayesian mixture model at both single and multiple time points for a sample of actual and simulated aerosol particle size distribution (PSD) data. For estimation of a mixture model at a single time point, we use Reversible Jump Markov Chain Monte Carlo (RJMCMC) to estimate mixture model parameters including the number of components which is assumed to be unknown. We compare the results of this approach to a commonly used estimation method in the aerosol physics literature. As PSD data is often measured over time, often at small time intervals, we also examine the use of an informative prior for estimation of the mixture parameters which takes into account the correlated nature of the parameters. The Bayesian mixture model offers a promising approach, providing advantages both in estimation and inference.

Modeling transport through an environment crowded by a mixture of obstacles of different shapes and sizes

Many biological environments are crowded by macromolecules, organelles and cells which can impede the transport of other cells and molecules. Previous studies have sought to describe these effects using either random walk models or fractional order diffusion equations. Here we examine the transport of both a single agent and a population of agents through an environment containing obstacles of varying size and shape, whose relative densities are drawn from a specified distribution. Our simulation results for a single agent indicate that smaller obstacles are more effective at retarding transport than larger obstacles; these findings are consistent with our simulations of the collective motion of populations of agents. In an attempt to explore whether these kinds of stochastic random walk simulations can be described using a fractional order diffusion equation framework, we calibrate the solution of such a differential equation to our averaged agent density information. Our approach suggests that these kinds of commonly used differential equation models ought to be used with care since we are unable to match the solution of a fractional order diffusion equation to our data in a consistent fashion over a finite time period.