A Bayesian approach to parameter estimation for a crayfish (Procambarus spp) bioaccumulation model

Bioaccumulation models are used to describe chemical uptake and clearances by organisms. Averaged input parameter values are traditionally used and yield point estimates of model outputs. Hence, the uncertainty and variability of model predictions are ignored. Probabilistic modeling approaches, such as Monte Carlo simulation and the Bayesian method, have been recommended by the U.S. Environmental Protection Agency to provide a quantitative description of the degree of uncertainty and/or variability in risk estimates in ecological hazards and human health effects. In this study, a Bayesian analysis was conducted to account for the combined uncertainty and variability of model parameters in a crayfish bioaccumulation model. After a 5-d exposure in the LaBranche Wetlands (LA, USA), crayfish were analyzed for polycyclic aromatic hydrocarbon concentrations and lipid fractions. The posterior distribution of model parameters were derived from the joint posterior parameter distributions using a Markov chain Monte Carlo approach and the experimental data. The results were then used to predict the distribution of chrysene concentration versus time in the crayfish to compare the predicted ranges at the different study sites.