Chemical contamination has been suggested as an important contributing factor to reptile population declines, but direct links are rarely reported. Population modeling provides a quantitative method to understand the long-term effects of contaminants on population persistence. We created a matrix model for Sceloporus lizards and investigated hypothetical toxic effects by reducing survival and reproductive parameters by 0 to 100% in 10% increments. We report effects on population growth rate (λ) and elasticity values for each stage due to these reductions. We then incorporated stochasticity to the model to simulate the variation seen in demographic data and quantified extinction risk. The deterministic model yielded a λ of 1.07 suggesting stability in some wild Sceloporus populations. A yearly reduction of 20 to 30% in demographic parameters was needed to push λ to decline in both our deterministic and stochastic simulations. Surprisingly, our baseline stochastic simulations had a 30% extinction probability despite a stable deterministic model. We tested three adjustments to the stochastic model, (1) increased survival/fecundity parameters, (2) higher starting densities, and (3) a density-dependent juvenile survival function. The model with density-dependent juvenile growth had the lowest extinction risk. Ultimately, 20 or 30% mortality every year is likely unrealistic, but our results provide insight in linking toxicity to population effects. Ultimately, very little reduction in demographics is needed to cause declines in these populations. Our generalized models provide important tools for screening-level risk assessment of chemical contamination, especially for taxonomic groups that tend to receive less research interest.
Keywords: Ecotoxicology; Lizard; Population modeling; Risk assessment.