Download A framework for assessing the biological risks of increasing salinity

Document related concepts

Statistics wikipedia , lookup

History of statistics wikipedia , lookup

Transcript
A Framework for Assessing the
Biological Risks of Increasing
Salinity in Victoria
W.J. Dixon
June 2007
Arthur Rylah Institute for Environmental Research
Technical Report Series No. 162
Arthur Rylah Institute for Environmental Research
Technical Report Series No. 162
A framework for assessing the biological risks
of increasing salinity in Victoria
William J. Dixon
June 2007
Published by: Arthur Rylah Institute for Environmental Research
Department of Sustainability and Environment
PO Box 137
Heidelberg, Victoria 3094
Australia
Telephone: (03) 9450 8600
www.dse.vic.gov.au/ari/reports
This publication may be cited as:
Dixon, W.J. (2007) A framework for assessing the biological risks of increasing
salinity in Victoria. Arthur Rylah Institute for Environmental Research, Technical
Report Series No. 162, Department of Sustainability and Environment, Heidelberg,
Victoria.
© The State of Victoria
Department of Sustainability and Environment 2007
This publication is copyright. Apart from any fair dealing for private study, research,
criticism or review allowed under the Copyright Act 1968, no part of this publication
may be reproduced, stored in a retrieval system or transmitted in any forms or by
any means, electronic, photocopying or other, without the prior permission of the
copyright holder.
ISBN 978-1-74208-052-9 (Print)
ISBN 978-1-74208-053-6 (Online)
ISSN 1835-3827 (Print)
ISSN 1835-3835 (Online)
Disclaimer
This publication may be of assistance to you but the State of Victoria and its
employees do not guarantee that the publication is without flaw of any kind or is
wholly appropriate for your particular purposes and therefore disclaims all liability
for any error, loss or other consequence, which may arise from you relying on any
information in this publication.
Authorised by the Victorian Government, Melbourne.
Abbreviations
ANZECC – Australian and New Zealand Environment and Conservation Council
ARI - Arthur Rylah Institute for Environmental Research
ARMCANZ - Agriculture and Resource Management Council of Australia and New
Zealand
CDF – Cumulative Density Function
DPI – Department of Primary Industries (Victoria)
DSE - Department of Sustainability and Environment (Victoria)
EC05 - Effective Concentration 5%: The concentration of a toxicant that results in a
5% effect
ECD - Environmental Concentration Distribution
FFG – The Flora and Fauna Guarantee Act 1988 (Victoria)
GLM - Generalised Linear Model
GLMM - Generalised Linear Mixed Model
LC50 - Lethal Concentration 50%: The concentration of a toxicant that results in a
50% lethal effect
LD50 - Lethal Dose 50%
LOEC - Lowest Observed Effect Concentration
LOEL - Lowest Observed Effect Level
LR - Likelihood Ratio
MCMC – Markov Chain Monte Carlo
ML – Maximum Likelihood
NOEC - No Observed Effect Concentration
PBA – Probability Bounds Analysis
PDF - Probability Density Function
SSD - Species Sensitivity Distribution
I
II
Abstract
Increasing salinity presents a significant threat to a wide range of native plants and
animals in Victoria. In order to ensure adequate conservation of our biodiversity
assets, management action is often required to mitigate or ameliorate the effects of
salinity. In many situations, however, little is known about which species, areas or
ecosystems might be most affected and how. Furthermore, management efforts must
often be prioritised in order to make the best use of limited resources for the
protection, enhancement and restoration of our biodiversity assets. Such decision
making processes can often be better informed through quantitative risk and
scenario assessments. In particular, explicit recognition of the uncertainties involved
in risk assessments can help to avoid erroneous decisions, increase transparency
and highlight critical information gaps. The integration of more objective and
defensible scientific approaches into the decision making process has many other
benefits, including satisfying the aims of various legislative policies and guidelines,
ensuring efficiency in government investment, facilitating a broader landscape
approach to land and water management, and maximising the potential for
biodiversity protection.
This report examines and develops risk assessment methods for predicting the
effects of increasing salinity on biodiversity assets. A detailed investigation of
various approaches to modelling the risks of salinity is presented and the relative
merits and limitations of each are discussed. In summary, techniques that have
previously been applied in this area of risk modelling have focused on the use of
scalar summaries of effect, ‘area under the curve’ formulations of risk and Species
Sensitivity Distributions (SSDs). However, these techniques involve assumptions that
are potentially incorrect and may limit their applications in risk assessment. This
report describes the development of a number of alternative methods for risk
modelling that require fewer statistical assumptions, correctly characterise and
propagate uncertainty and can be used to develop a greater ecological understanding
of
potential
effects.
The
methods
presented
range
from
basic
graphical
representations of tolerance data to complex probabilistic models. In particular,
Probability Bounds Analysis and Bayesian Markov Chain Monte Carlo (MCMC)
techniques are recommended as the basis for improved risk modelling approaches.
The quality and availability of data on which to base risk assessments is also
identified as an important factor contributing to uncertainty in risk management
decisions. It is recommended that these issues be given careful consideration before
deciding on a course of action. Details and examples of suggested risk assessment
techniques are given in the Appendix and also in a series of associated case studies
(Dixon 2007a, b; Dixon & Scroggie 2007).
3
4
Contents
ABBREVIATIONS............................................................................................................ I
ABSTRACT ............................................................................................................... III
ACKNOWLEDGEMENTS ............................................................................................. 1
1 INTRODUCTION AND BACKGROUND ................................................................. 2
AIM OF SALT IMPACT MODELLING .......................................................................................... 3
RECOMMENDATIONS OF STAGE ONE REVIEW ............................................................................. 4
2 REVIEW OF METHODS FOR SALT IMPACT MODELLING. ..................................... 6
OVERVIEW OF CURRENT METHODS FOR MODELLING SPECIES SENSITIVITY DISTRIBUTIONS .................... 6
Parametric Modelling – Probability Distributions ............................................................. 6
Bootstrap Techniques.................................................................................................. 7
Non-parametric Methods ............................................................................................. 8
SPECIES SENSITIVITY DISTRIBUTIONS IN PROBABILISTIC RISKS ASSESSMENTS ................................... 9
3
SOURCES OF UNCERTAINTY IN SPECIES SENSITIVITY DISTRIBUTIONS .......... 11
DATA ASSOCIATED UNCERTAINTY ........................................................................................ 11
PREVIOUS APPROACHES TO THE TREATMENT OF UNCERTAINTY IN TOXICITY BASED RISK ASSESSMENTS. 13
LIMITATIONS OF PREVIOUS APPROACHES TO UNCERTAINTY ANALYSIS AND SSD MODELLING IN
ECOTOXICOLOGY ............................................................................................................. 14
4
REFINED APPROACHES FOR SALT IMPACT MODELLING ....................................18
METHODS FOR UNCERTAINTY PROPAGATION IN RISK ASSESSMENTS .............................................. 18
Interval Arithmetic and Probability Bounds Analysis ..................................................... 18
Bayesian Markov Chain Monte Carlo Estimation Techniques .......................................... 19
UNCERTAINTY CHARACTERISATION AND PROPAGATION IN SSD MODELLING ...................................... 20
5 A GENERALISED LINEAR MODEL APPROACH TO PREDICTIVE EXPOSUREEFFECTS MODELLING ............................................................................................. 23
Probability Bounds Analysis........................................................................................ 24
Bayesian Markov Chain Monte Carlo Methods .............................................................. 26
6 SUMMARY AND CONCLUSIONS. ........................................................................... 29
APPENDIX 1: BASIC METHODS FOR SALT IMPACT MODELLING ............................. 32
A.1 - LOOKUP TABLES AND PLOTS ..................................................................................... 32
A.2 – BURRLIOZ SOFTWARE EXAMPLE ................................................................................ 35
A.3 - PARAMETRIC SSD MODELLING .................................................................................. 37
REFERENCES ........................................................................................................... 39
5
Acknowledgements
Particular thanks are extended to Dr Robin Hale for his editorial and formatting
assistance. A number of other people contributed to the production of this report
and their efforts are greatly appreciated: Prof. Mark Burgman (University of
Melbourne), Dr Scott Ferson (Applied Biomathematics Corp.) Phil Papas (ARI), Tom
Ryan (Environous), Dr Micheal Scroggie (ARI), Dr Charles Todd (ARI), and Dr Jenny
Wilson (DSE). This research was funded by the National Action Plan for Salinity and
Water Quality (NAP) and the National Heritage Trust (NHT).
1
1
Introduction and Background
Salinity presents a significant threat to much of Victoria’s biodiversity, yet little is
known about how effects might occur and what can be done to alleviate them. In
order to ensure suitable protection, enhancement and restoration of our natural
resources, land and water managers are often tasked with deciding how best to
allocate limited resources. It is increasingly recognised that such decisions may be
better informed through quantitative risk and scenario assessments (Thompson &
Graham 1996; Thompson 2002; Burgman 2005). The integration of these more
objective and defensible scientific approaches into the decision making process has
many benefits, including satisfying aims within relevant state strategies, ensuring
efficiency in government investment, using a broader landscape approach to land
and water management, and maximising the potential for biodiversity protection.
The creation of the Salt Sensitivity Database, available through Land and Water
Australia (LWA)1, has provided an important source of collated information on the
salinity tolerance of a range of Australian species (Bailey et al. 2002). However, few
tools are available that can be used to combine this sensitivity data with projections
of increases in salinity in order to determine the risk to species, communities or
ecosystems. In particular, approaches for evaluating uncertainties in these types of
risk assessments have, to date, not been adequately described. Consequently,
decisions must currently be based on risk estimates involving unknown degrees of
confidence.
This document describes the refinement process and outcomes for the Salt Impact
Modelling approach developed as part of Stage 2 of the National Action Plan for
Salinity and Water Quality and the National Heritage Trust funded project: Strategic
Frameworks to Indicate Biodiversity thresholds to Salinity (hereafter ‘The Project’).
Stage 1 of The Project consisted of an initial exploration of different methods for
predicting the effects of salinity on biodiversity, including basic modelling
techniques and exploration of possible vital attributes or aggregation models that
could be used to extrapolate between or provide analogues of different organisms.
The approach sought to utilise data from the LWA database to develop models that
could guide managers in making planning decisions to protect biodiversity assets
from salinity (Cant et al. 2003). The development and use of Salt Impact Models
aimed to “identify which biodiversity assets are most at threat from salinity in order
to target research and management resources to protect those assets” (Cant et al.
2003 p7). An initial trial of the approach, and subsequent review (Webb & Hart 2004),
1
www.rivers.gov.au/Tools_and_Techniques/Salt_Sensitivity_Database/index.aspx
2
revealed that such tools are likely to be useful to managers but that the modelling
techniques
required
further
development.
A
series
of
regional
stakeholder
workshops, consultation with a number of relevant experts and case study work to
test the utility of the refined SIM
approach also contributed to the process
of
refining the modelling approach.
This report describes the outcomes of Stage 2 of The Project on developing ‘refined’
Salt Impact Modelling approaches. A review of methods that are currently applied to
modelling the environmental risks of hazardous agents is provided first. Second, the
sources of uncertainty in ecotoxicological risk assessments are outlined and the
problems and limitations of previous approaches to analysing this uncertainty
discussed. Third, a number of new, alternative techniques for risk modelling are
described and presented as recommended methods for salt impact modelling. These
new approaches focus on the use of Probability Bounds Analysis (PBA) and Bayesian
Markov Chain Monte Carlo (MCMC) methods of uncertainty propagation and are
intended to provide a valuable addition to the suite of assessment and management
tools available for mitigating and ameliorating the risks of salinity to key
biodiversity assets. To aid in understanding the approaches developed, a series of
associated case studies have been performed that describe the methods in further
detail and give examples of their application (Dixon 2007a, b; Dixon & Scroggie
2007)2. Further detail on some basic approaches to modelling salinity effects are also
given in Appendix 1.
Aim of Salt Impact Modelling
The aim of Salt Impact Modelling is to predict the likely effects of changes in levels
of salinity at a particular site for a certain species or group of species. Regional scale
spatial predictions are also possible using the techniques presented here, however
these will invariably involve greater uncertainties due to the lack of suitable
information on the distribution of organisms and the effects of salt at these scales.
The predictions made through Salt Impact Modelling are intended to provide useful
information to decision makers in managing the effects of salinity and can also form
a basis for quantitative risk comparisons between places or situations in order to
assist in the prioritisation and allocation of resources and management works. Due
to the uncertain nature of much of the available data, estimates of confidence and/or
uncertainty in risk estimates will be critical in describing and communicating the
reliability that can be ascribed to predictions, and therefore subsequent management
2
Associated material can be found at www.dse.vic.gov.au/ari/reports
3
decisions. Further refinement of predictions can occur through monitoring, the
integration of more complex ecological models, and additional field research.
The revised Salt Impact Models (SIMs) will also aid in a more strategic approach to
conservation management by providing quantitative information about which
biodiversity assets are most at risk from increasing salinity and the relative
uncertainty surrounding the data on which management decisions can be made. The
SIM approach will enable more informed and objective decisions to be made based
on the likely outcomes of a range of different management scenarios and the likely
impact of any management changes on key biodiversity assets. The SIMs help to
achieve State and Regional conservation goals by providing a rigorous method of
assessing risks, in addition to helping manage threatened species and communities
because specific species can be targeted for assessment.
Recommendations of Stage One Review
A review of Stage 1 of the Project was conducted by Dr Angus Webb & Professor
Barry Hart of the Water Studies Centre, Monash University. This review made a
number of suggestions about refining the modelling process. These included:

The general Probabilistic Risk Assessment (PRA) approach used was highly
supported but it was recommended that some changes be made to the methods
to ensure credibility and suitability of the Project;

Further clarification of the advantages and disadvantages of the approach used
was required, in particular, discussion of relative and absolute risks, as well as
the limitations and uncertainties of the methods used;

Consideration and discussion of the relative merits and inferential consequences
of using different data sources. Specifically, issues associated with using field
and laboratory tolerance information and the ecological relevance of exceeding
particular thresholds or likelihoods of effect;

Validation and statistical testing of assumptions associated with parametric or
other modelling approaches;

Consideration and development of approaches to link together the results of the
SIMs with salinity models such as the Catchment Assessment Tool (CAT); and

Use of a case study approach to test models and validate predictions to develop
credibility and acceptability of the approach.
4
These points formed the starting point for this review and the inclusion of the
alternative modelling approaches described in this report.
5
2
Review of Methods for Salt Impact Modelling
Overview of Current Methods for Modelling Species Sensitivity
Distributions
The designation of ‘safe’ levels of exposure to hazardous agents has become a
fundamental component of the regulation, control and management of both human
and environmental health around the world (Posthuma et al. 2002b). In human and
occupational health situations these (safe) levels are commonly described as
exposure standards, such as maximum allowable dietary or inhalation doses. In
environmental situations such levels are routinely expressed as ‘environmental
quality criteria’, which specify the amount of a hazardous agent below which adverse
effects are considered unlikely to occur (or considered acceptable) (Suter 1993; Rand
et al. 1995; Cairns Jr 1999). Although commonly used to describe synthetic
chemicals, the term ‘hazardous agents’ is used here to refer to a wide range of
potentially toxic natural and anthropogenic phenomena, including salinity.
The most common approach currently used to evaluate the environmental
‘acceptability’ of levels of hazardous agents for environmental applications involves
the construction of a Species Sensitivity Distribution (SSD). These are statistical or
mathematical models that relate an amount of hazardous agent with an effect,
expressed as a proportion of a ‘group of species’ responding (Posthuma et al. 2002a).
This group of species may be defined a number of ways, such as those present at a
specific site or all species of a particular taxonomic type. While empirical models are
probably the most commonly used, a limited number of parametric and nonparametric approaches have been explored, with fewer still of these incorporating
uncertainty analysis (Verdonck et al. 2001b; Posthuma et al. 2002a). In the SSD
approach the risk to the defined group of species is considered to be the overlap of
the SSD (the likelihood of effect from exposure) and Environmental Exposure
Distribution (ECD, the likelihood of exposure to the hazardous agent).
Parametric Modelling – Probability Distributions
Species Sensitivity Distributions, as popularised by Aldenberg and Slob (1993), are
predominately modelled by fitting parametric relationships, usually as a Cumulative
Distribution Function (CDF), to sets of tolerance data. Techniques for performing
parametric distribution fitting are included as standard in many statistical packages
and this method has been widely applied in regulatory situations in Australia and
around the world (Hart et al. 1999; ANZECC & ARMCANZ 2000; van Straalen & van
Leeuwen 2002; Van Sprang et al. 2004).
6
Parametric SSD methods have a number of advantages, including: (1) relatively little
data is required to develop a model, (2) the fitting techniques are included as
standard in many statistical packages and (3) the methods of testing and checking
model adequacy and fit are relatively well understood. Unfortunately, these methods
often require a number of assumptions to be made that may limit their application.
These issues are discussed further in Section 3.
In Australia, parametric methods for SSD modelling were applied by Fox (1999) and
Shao (2000) in developing the modelling approach used in the Australian and New
Zealand Guidelines for Fresh and Marine Water Quality (ANZECC & ARMCANZ 2000).
This method involved fitting a Burr Type III distribution, a three parameter
exponential distribution of which the log-log distribution is a special case (Fox 1999;
ANZECC & ARMCANZ 2000; Shao 2000). This distribution was considered to have
the advantage of being able to encompass a wide range of shapes (although there
may be some problems of colinearity with three parameter distributions). As part of
the ANZECC guidelines, Shao (2000) developed a software
program,
“BurrliOZ”,
which fits the Burr Type III distribution to scalar values and calculates confidence
intervals. Unfortunately, the Burr Type III distribution is not included in many
statistical packages and the BurrliOZ program has limited goodness of fit
applications, making determination of model adequacy difficult. An example using
the BurrliOZ software is given in Appendix 1 (Section A2).
In other Australian examples, Webb and Hart (2004) used log-normal models to
describe the effects of salinity at a number of sites in the Goulburn-Broken
catchment and applied graphical methods for propagating uncertainty, while Kefford
et al. (2005) and Kefford et al. (2006) investigated the use of empirical SSD models
for macro-invertebrates.
Bootstrap Techniques
Bootstrap and other re-sampling techniques are increasingly being applied in risk
modelling because of their ability to cope with multi-modal data sets, lack of
assumptions about distribution type and shape, and the relative ease of computation
of associated statistics (Manly 1997; Frey & Burmaster 1999; Bixio et al. 2002). Jagoe
and Newman (1997) and Newman et al. (2000) have suggested applying bootstrap
techniques to modelling SSDs. These authors propose using bootstrap estimation
methods to overcome the problems and assumptions associated with many
parametric distributional analyses. In an example of the application of boot strap
techniques, 30 toxicity data sets were used to compare between the utility of
lognormal, log-logistic and Gompertz distributions for modelling SSDs (Jagoe &
7
Newman 1997; Newman et al. 2000). The number of values in each data set ranged
from 20 to 91, with an average of 32 and 80% of data sets having less than 40 values.
The authors found that half of the 30 data sets failed the test of lognormality
(alpha=0.05). Goodness of fit tests using Pearson’s X2 statistic indicated that a loglogistic model was a better fit to the data twice as often as the lognormal
distribution, but that the Gompertz distribution was the best fit overall. Most data
sets were found to be multimodal, which was considered to be the result of
similarities between taxonomic groups within the data sets. A strong correlation was
found between the bootstrap and log-normal percentiles (r2 = 0.954). The authors
concluded that assumptions of log-normality were not supported, that the SSD
approach was useful given the limited data available, and that it provided a direct
means of comparing cumulative exposure concentration distribution to cumulative
SSD.
In a response to Newman et al. (2000), Verdonck et al. (2001a) highlight a number of
methodological flaws in the bootstrap approach that was applied. In particular,
because of the small number of observations in each dataset, the sample size for
each bootstrap sample was necessarily larger than those in the original data set. This
meant that all data values were used in the construction of the boot-strap
distribution, not a random sub-sample, and that some values were sampled multiple
times. Although Newman (2000) refers to the use of this approach in the first (1991)
edition of Manly’s work on bootstrap and randomisation methods, this section has
been removed from the second edition as it is not a particularly rigorous practice
(Manly 1997). For many SSD modelling situations the number of data points may be
very few, creating problems with both bootstrap and randomisation approaches.
Newman et al. (2000) further applied this re-sampling approach to calculating
confidence limits for percentiles of the SSD, although the reasoning involved in this
method is again circular, in that the more re-samples taken, the narrower the
confidence intervals were, as also noted by Verdonck et al. (2001a). It is thus evident
that these particular applications of bootstrap and re-sampling methods to SSD
modelling involve inferential and reasoning problems that should be avoided.
Non-parametric Methods
A number of other approaches to modelling SSDs have been proposed. These include
such methods as the ‘Hazen Plotting System’ (Verdonck et al. 2001b), fitting
empirical CDFs, and other non-parametric methods (van der Hoeven 2001; Verdonck
et al. 2001b). While some of
these methods do overcome particular limitations
associated with parametric techniques, they unfortunately tend to require nonstandard statistical packages and have not been widely applied. Further, many of the
8
applications of these methods to date have focused on calculating particular
percentiles of the distribution of effect and not on characterising the whole
distribution or its associated uncertainties. The use of empirical distributions, where
possible, does have the advantage that assumptions of underlying shape are not
required.
Species Sensitivity Distributions in Probabilistic Risks Assessments
The predominant approach to deriving probabilistic measures of the environmental
risk of hazardous agents involves specifying both the Species Sensitivity Distribution
(SSD)
and
Environmental
Concentration
Distribution
(ECD)
as
probability
distributions (Aldenberg et al. 2002; Forbes & Calow 2002; Posthuma et al. 2002a; de
Zwart 2003). In this formulation, the risk is defined as the overlap between these two
distributions, which is also the integral or Area Under the Curve (AUC) of the joint
probability distribution (van Straalen & Denneman 1989; Aldenberg et al. 2002;
Verdonck et al. 2002). Graphical interpretations of the AUC have been described
(Soloman et al. 2000; Soloman & Sibley 2002; Webb & Hart 2004), however these
approaches have been criticised as mathematically incorrect methods for uncertainty
propagation
(Verdonck
et
al.
2003).
A
number
of
different
mathematical
formulations of the AUC have also been described (van Straalen 1990; Cardwell et al.
1993; Aldenberg & Jaworska 2000; Aldenberg et al. 2002; van Straalen 2002a;
Warren-Hicks et al. 2002). Aldenberg et al. (2002) reviewed these different
approaches and found them to be mathematically equivalent expressions of the
integral of the joint probability distribution, which they also describe as the chance
of one distribution exceeding the other, expressed as:

p(ECD SSD)   PDFECD (x).CDFSSD (x)dx
(1)
or

p(ECD SSD)   (1CDFECD (x)).PDFSSD (x)dx
9
(2)
Where:
-
ECD is the Environmental Concentration Distribution,
-
SSD is the Species Sensitivity Distribution,
-
PDF is a Probability Density Function,
-
CDF is a Cumulative Distribution Function, and
-
and x is the level of a hazardous agent.
Aldenberg et al (2002) further describe a simplification of this risk formulation,
termed the ‘Probabilistic Risk Quotient’:
Risk = p(ECD > SSD)
= p(ECD / SSD > 1)
= p(log(ECD) – log(SSD) > 0)
If both the ECD and SSD are assumed to be random (log)normally distributed
variables, then the result is also a normal distribution with mean, µ = µ
standard deviation, σ=√(σ
σ
SDD
) + (σ
2
ECD
) , where µ
2
SDD
ECD
and µ
SDD
ECD
-µ
are the means, and σ
SDD
ECD
and
and
the standard deviations, of the log-transformed ECD and SSD respectively
(Aldenberg et al. 2002; Verdonck et al. 2003). Although this solution provides a
probabilistic measure of risk, it does not allow for propagation of measures of
uncertainty and the assumption of (log)normality is generally untested.
10
3
Sources of Uncertainty in Species Sensitivity
Distributions
The evaluation of uncertainty in risk estimates is of central importance to ensuring
rigor, transparency and success in the risk management process (Morgan & Henrion
1990; Winkler 1996; Suter et al. 2000; Andrews et al. 2004; Burgman 2005). Despite
widespread
acknowledgment
of
the
need
for
uncertainty
analysis
in
risk
assessments (NRC 1993, 1994; Pastorok 2002; Burgman 2005), to date there have
been very few applications of such techniques in toxicological risk assessments
(Aldenberg & Jaworska 2000; Posthuma et al. 2002a; Verdonck et al. 2002; Van
Sprang et al. 2004). Furthermore, it has been recognised that current probabilistic
formulations involve a number of untested assumptions and may not be amenable to
uncertainty analysis (Suter 1998; Forbes & Calow 2002; Newman et al. 2002; Suter et
al. 2002; van Straalen 2002a; van der Hoeven 2004).
Uncertainty in risk assessments arises from a multitude of sources including natural
variability in the sensitivity of organisms to toxicants, confounding and modifying
factors of toxicity, systematic measurement error, lack of knowledge about
mechanisms determining toxicity, and extrapolation across spatial, temporal and
organisational scales (Chapman et al. 1996; Chapman et al. 1998; Cairns Jr 1999). A
distinction is often made between aleatory uncertainty, or the non-reducible inherent
variability in a parameter, and epistemic uncertainty, or incertitude about the
process or situation that should, at least in part, be reducible (Ferson & Ginzburg
1996; Regan et al. 2002a; Ferson et al. 2003).
Data Associated Uncertainty
Perhaps the most significant source of uncertainty in risk assessments involves the
data (or lack there of) on which they are based. Issues surrounding the comparability
of differing spatial and temporal scales of observation and exposure, test endpoints,
taxonomic resolution and a general lack of information on the tolerance of many
species to hazards contribute to uncertainty in the predictions made from risk
models.
Information about the effects of salinity generally comes from either toxicity studies
or field observations. Toxicity studies performed in laboratories, glass houses or
field trials are designed to determine a particular magnitude of effect at different
levels of exposure to a hazardous agent. These studies can also provide
more
detailed information of use in risk assessments, such as the toxic mode of action or
11
the presence of acute responses thresholds. However, because many of these studies
are necessarily performed over short time periods and away from the complicating
factors found in an organism’s natural environment, their applicability to field
situations is often unknown. A multitude of factors may affect an organisms’
sensitivity to a toxicant, and hence the relevance of laboratory/toxicity data to real
world situations. These include:

The length of the exposure period;

The rate at which exposure is increased (i.e. direct transfer vs. gradual change);

The frequency of exposure;

The season or life history stage at which exposure occurs;

The endpoint (i.e. acute vs chronic);

The presence of other stressors, such as predators, chemical pollutants, drought,
food availability; and

Water quality associated modifying factors such as carbonate concentration, pH
or temperature.
Despite these limitations, toxicity experiments do provide direct evidence of the
effect of a toxicant and a means by which to test hypotheses about its effects.
Observations of field tolerance, on the other hand, can potentially provide a more
realistic and encompassing indication of likely effects, although the relevance of
such measures of ‘toxic effect’ may also be affected by the factors listed above.
Unfortunately, few field studies are designed to determine the maximum tolerance
of an organism or community to salinity. Instead, much of the data on the field
tolerance of organisms comes from associated research, such as distributional
surveys, where salinity levels have been recorded. Many of these surveys have only
been conducted once over limited areas and thus only provide a snap-shot of the
range of variable environmental conditions that may be experienced over a wider
spatial and temporal extent. Consequently, this information is often uninformative
about where the observed salinities lie within the range of a species’ occurrence and
how long that level can be tolerated.
Issues of which value to use for a species when there are multiple values available
are also currently unresolved. Arguments might be made for using values from the
closest location or the lowest values to ensure that risk assessments are
conservative. These and other approaches to modelling and propagating uncertainty
are discussed in more detail in the next section.
A final and particularly important source of uncertainty in predicting the effects of
salinity relates to the consideration of ‘salinity’ as a toxicant. This is particularly
12
problematic as ‘salinity’ comprises a complex mixture of many different chemical
ions whose relative composition may vary considerably from place to place. Very
little is known about the toxicity of these different ions on native species, let alone
their combined long-term effects, therefore careful consideration should be given to
the relevance and comparability of different mixtures of salts.
Previous Approaches to the Treatment of Uncertainty in Toxicity
Based Risk Assessments
In the risk formulations described above (Section 2.2), only variability in the mean
response is considered and not other sources of knowable uncertainty associated
with the variables,
parameterized
by
such as
probability
inherent
variability that
distributions
or
is greater than that
uncertainty
associated
with
assumptions of model structure. Although the potential effects of other sources of
uncertainty in SSDs have been widely recognised (Suter 1998; Newman et al. 2000;
Forbes & Calow 2002; Newman et al. 2002; Suter et al. 2002), relatively few
approaches have been suggested that include these in risk estimates (Aldenberg &
Jaworska 2000; Aldenberg et al. 2002; Verdonck et al. 2002; Verdonck et al. 2003;
Van Sprang et al. 2004).
Verdonck et al. (2002; 2003) and Van Sprang et al. (2004) describe a Monte Carlo
approach to estimating uncertainty in probabilistic risk quotients, which attempts to
account for uncertainty in both the SSD and ECD. In this approach, successive
samples are drawn from both distributions and used to create a density of values for
a distribution of risk quotients. This density can then be used to obtain the
likelihood that the risk quotient is equal to one. Similar Monte Carlo approaches
have been widely applied in human health risk assessments (Finkel 1995; Brattin et
al. 1996; Thompson 1999; Andrews et al. 2004).
Aldenberg et al. (2000) and Aldenberg et al (2002a) describe a sampling theory
extrapolation approach to
estimating
confidence intervals for percentiles of the
normal distribution. This uses extrapolation factors derived from the non-central tdistribution that relate the number of observations in the sample with estimates of
confidence in the ‘true’ population percentile. Lookup tables of the extrapolation
factors for 90%, 2-sided confidence intervals are provided for various percentiles of
interest (Aldenberg & Jaworska 2000). Although calculations were made for both
deriving the concentration hazardous to a percentile of a population (the ‘forward’
situation) and the fraction affected at a particular concentration (the ‘reverse’
situations), only uncertainty in the SSD was considered.
13
There are also an increasing number of studies that have investigated uncertainties
in larger and more complex models of exposure and ecological risk (see: Ferson et al.
1996; Fogarty et al. 1996; Burmaster & Wilson 1998; Moore et al. 1999; Jager et al.
2001; Linkov et al. 2001; Regan et al. 2002b; Regan et al. 2002c; Bartell et al. 2003;
Bates et al. 2003; Breton et al. 2003; Crane et al. 2003; Landis 2003; Pastorok et al.
2003; Merrick et al. 2005). Many of these studies begin to address the ecological
aspects of the effects of chemicals, however to date there has been little
consideration of uncertainty in concentration-response relationships within these
models or its impact in overall measures of risk. Instead, toxicity information has
generally been represented as a scalar toxicity reference value or dose.
Limitations of Previous Approaches to Uncertainty Analysis and SSD
Modelling in Ecotoxicology
While the approaches outlined above have provided a useful grounding for SSD
based risk assessments (Forbes & Calow 2002; Posthuma et al. 2002a), they involve a
number of assumptions that may be incorrect and limit their application in risk
management (Suter 1998; Suter et al. 2002; van der Hoeven 2004):
First, the proposition that there is a ‘safe level’ of exposure has been the subject of
some criticism. It has been argued that for some hazardous agents there is no safe
level of exposure (Newman et al. 2000). While this criticism is probably valid for
many anthropogenic toxicants, in cases involving salinity it would be reasonable to
assume that there is some safe level of exposure, given that the compounds involved
occur naturally and many are essential to survival.
Second, the aim of SSD modelling is to estimate a certain percentile of species to be
protected, such as the lower 5th percentile (i.e. assumes that 95% of species will be
protected). While this may appear to be an honourable goal, the basic SSD modelling
process is not informative about which species are lost or the ecological
consequences of the loss of those particular species. This would
require
identification of the expected species to be lost and consideration of the ecological
effects of that loss within the community. However, very little is known about the
behaviour of many of the systems of interest, further adding to the uncertainty in
decisions made.
A third limitation of the SSD approach is that models might only be informative for
species for which there is data. If a community contains 100 species but there is only
information for about 50 of these, predictions that describe the proportion
of
species affected will only be based on the 50 species for which there is data (so there
14
is potentially another 50 species that could be affected that were not included in the
assessment). While the process of fitting probability distributions to the data is used
to extrapolate from a sample of species to a larger number, the test species used are
generally not an unbiased random sample from a larger population of potentially
affected species. Additionally, for many situations there is a high degree of
taxonomic uncertainty, both in the description of different species and their
distribution.
Fourth, a number of authors have noted that the treatment of an exposure-response
relationship as a random variable (i.e. as represented as a probability distribution) is
mathematically incorrect (Suter 1998; Forbes & Calow 2002; Suter et al. 2002; van
Straalen 2002b; van der Hoeven 2004). Particular problems with this probabilistic
description of the SSD include:
-
test species are not randomly selected;
-
interpretation of the underlying variability generating mechanisms is
ambiguous, as is the interpretation of the differentiable relationship between
the PDF and CDF of the SSD;
-
different species are treated as equally important and non-dependant on each
other in a community; and
-
interactions between species are not considered.
A fifth problem is the use of scalar summaries of exposure-response distributions.
The most common of these is referred to as the NOEC or No Observed Effect
Concentration, and is often defined as the highest exposure concentration in a
toxicity test that did not produce a significantly different response from the control
(Rand 1995). A number of other similar ‘no-effect’ summary measures are commonly
used, however these values have been found by a number of authors to be statistical
and experimental artefacts, and to lead to misleading interpretations (Laskowski
1995; Chapman et al. 1996; van der Hoeven 1997; Yanagawa & Kikuchi 2001; Suter et
al. 2002; van der Hoeven 2004). Although in practice it has often been necessary to
base risk assessments on scalar summaries due to a lack of previously published
whole datasets or other measures of effect, the use of scalar summaries is
uninformative with regards to levels of uncertainty and variability in the estimates
produced (Hoekstra & van Ewijk 1993; Chapman et al. 1996; van der Hoeven 1997;
Newman et al. 2000). This has the flow-on effect that uncertainty and variability in
the response to a toxicant is ignored in the calculation of risk estimates, even though
it is obviously a fundamentally important consideration if the intention is to predict
the true range of values over which a species or community could be susceptible to a
toxicant. A number of authors have suggested the use of Effect Concentration (EC)
15
values and whole exposure-response distributions as more informative alternatives
(Chapman et al. 1996; de Bruijn & Hof 1997; van der Hoeven 1997; Chapman et al.
1998), although there have been only a few applications of these in risk assessments
to date (Nayak & Kundu 2001; Englehardt 2004).
A sixth limitation of previously applied SSD approaches involves statistical and
mathematical assumptions required in the formulations of risk described above
(Section 2.2). For example, it has been shown that assumptions of one underlying
distributional type, such as the normal or log-normal, rarely hold when applied to
multiple data sets (Newman et al. 2000; van der Hoeven 2001, 2004). Importantly,
discrepancies caused by these assumptions often have the greatest implications for
values in the tails of a distribution, which is generally the focus of risk assessments
(Hattis 1990; Newman et al. 2000). The assumption of ‘independence’ between
parameters in probabilistic calculations is also often untested or invalid and can lead
to an underestimation of risks (Ferson 2002; Ferson & Hajagos 2004). This can cause
significant problems with many applications of Monte Carlo techniques, where
parameters in a model are assumed to be independent of each other and the
inclusion of correlations between them can be difficult (Ferson 1996; Nayak & Kundu
2001; Ades & Lu 2003; Ferson et al. 2003; Berleant & Zhang 2004; Englehardt 2004;
Ferson & Hajagos 2004). Further limitations arise from the formulation of the risk
model as the closed form joint-integral of two probability distributions, as it is
difficult to expand the model to incorporate other factors, such as modifiers of
toxicity or exposure, random effects or population growth rates. Also, exact
solutions to the integral are only available for certain combinations of distributions.
Finally, methods for propagating uncertainty, such as those discussed below, are in
practice mostly limited to arithmetic operations and may not be appropriate for
many graphical and integration based model formulations (Ferson & Ginzburg 1996;
Ferson et al. 2003; Ferson & Hajagos 2004).
Further to these points, it is important to recognise that:
-
Sensitivity Distributions are really a function that describes the expected
‘map’ or effect for increasing levels of exposure. Standard portrayals of an
SSD show its values over the full range of those possible, whereas the risk is
the result of the SSD function for particular levels of exposure;
-
CDFs must be increasing and thus cannot incorporate alternative exposureresponse functions (such as hormesis or avoidance);
16
-
the use of percentile standards inherent in these approaches is based on the
assumption that the loss of some proportion of species is acceptable, and is
non-informative about which species may be lost;
-
the focus on changes in the number of species within a community assumes
that effects on populations are of no consequence.
While it has been argued that the probabilistic interpretation of species sensitivity
has some utility in representing uncertainty in responses across species and that it is
in the interpretation of the results of such processes that problems arise (Forbes &
Calow 2002; van Straalen 2002a), changing the interpretation of the SSD does not
alter the mathematical validity.
17
4
Refined Approaches for Salt Impact Modelling
This section describes three refined approaches for Salt Impact Modelling that
provide probabilistic measures of risk and incorporate uncertainty propagation.
Following on from methods developed in Dixon (2005) the first of these concerns the
application of Probability Bounds Analysis (PBA) to uncertainty characterisation in
SSD models. The other two methods involve the application of PBA and Bayesian
MCMC techniques for propagating uncertainty in generalised linear risk models of
population effects (Dixon 2005). While the two uncertainty propagation approaches
described here are not the only ones that may be appropriate for modelling salinity
risks, they do offer a number of advantages over the standard Monte Carlo based
approaches previously applied to these problems. A general background and
overview of the techniques is first given, before the particular applications of these
methods are outlined in the context of risk modelling. More detail on the application
of these methods is provided in a series of associated case-studies (Dixon 2007a, b;
Dixon & Scroggie 2007).
Methods for Uncertainty Propagation in Risk Assessments
Interval Arithmetic and Probability Bounds Analysis
Probability Bounds Analysis (PBA) is a combination of interval arithmetic and
probability theory that provides a rigorous method for propagating uncertainty
through calculations (Ferson 2002). In particular, PBA provides solutions to problems
involving unknown dependencies between variables and uncertainties in the exact
nature of distributions (Ferson 2002; Ferson & Hajagos 2004). The solving of some
key algorithms used in the computation of probability bounds (Yager 1986;
Williamson & Downs 1990; Ferson et al. 2003) and the recent development of
software that makes use of these techniques (RAMAS RiskCalc) 3, has enabled interval
and PBA approaches to become viable techniques for use in risk assessments (Ferson
2002; Regan et al. 2002b; Regan et al. 2002c; Ferson & Hajagos 2004; Regan et al.
2004).
Conceptually, PBA could be considered an extension of interval arithmetic. If two
uncertain numbers are considered, A and B, defined by the respective intervals [a ,1
a ] and [b , b ], then their sum must lie within the bounds [a +b , a +b ], and their
2
1
2
1
1
2
2
product within the bounds [min(a b , a b , a b , a b ), max(a b , a b , a b , a b )] (Ferson
1
1
1
2
2
1
2
2
1
1
1
2
2
1
2
2
2002). Similar rules exist for other arithmetic operations, and analogous solutions
3
RiskCalc is available from www.ramas.com
18
have been found for convolutions of bounds on whole distributions, referred to as a
probability box (p-box) (Yager 1986; Williamson & Downs 1990; Ferson et al. 2003).
PBA has the particular advantage over other uncertain number techniques that it is
not necessary to assume independence, or to know the dependency, between the
parameters in a model because the bounded solution is guaranteed to encompass
the true value (Ferson & Hajagos 2004). Because p-boxes describe a non-parametric
region of probability space and knowledge of model structure within the box is not
required, various aspects of model and parameter uncertainty can be investigated
simultaneously. The benefit of this is that different types or sources of uncertainty
may be propagated in the same calculation and that the resultant bounds on a
variable are the best possible, given the information available (Ferson et al. 2003;
Ferson & Hajagos 2004). PBA does have the limitation, however, that the relative
likelihood of values occurring inside or outside the bounds specified is not resolved.
To date there have been very few applications of intervals and PBA techniques in risk
assessments of hazardous agents and even fewer in ecotoxicology. Regan et al.
(2002c) applied PBA with an exposure model used to determine soil screening levels
for terrestrial wildlife. Risk was expressed as hazard quotients and the sensitivity
distribution was represented as a single, minimum toxicity reference value. This
study compared PBA with Monte Carlo approaches and showed that significant
discrepancies could result from unjustified assumptions of independence between
parameters (Regan et al. 2002c). In another application, Regan et al. (2002b)
compared two-dimensional Monte Carlo and PBA of a food web exposure model,
again using a minimum value toxicity reference value, and found similar
underestimation of uncertainties even with two-dimensional Monte Carlo techniques
(Regan et al. 2002b).
Bayesian Markov Chain Monte Carlo Estimation Techniques
With the assistance of increased computing power and the refinement of algorithms
used in Markov Chain Monte Carlo (MCMC) estimation procedures, Bayesian
approaches to data analysis and predictive modelling are becoming a powerful
alternative to frequentist or Maximum Likelihood Estimation (MLE) approaches
(Godsill 2001; Gelman et al. 2004; Lee 2004). The particular utility of MCMC methods
in risk assessment applications lies in the relative ease with which correlation
structures between parameters can be captured and propagated through models
(Ades & Lu 2003).
Bayesian approaches to data analysis have been well documented elsewhere
(Congdon 2001, 2003; Gelman et al. 2004; Lee 2004) and so shall only be covered
19
briefly here. Bayes’ theorem states that the posterior probability density p(θ|y) of a
parameter θ, given data y, is related to the prior probability of the parameter
independent of the data p(θ), and the likelihood of the data given the parameter
p(y|θ):
p(θ|y)  p(θ) x p(y|θ)
This general form can be extended to apply to each of the parameters in a model and
to make predictions for y(ỹ):
p(ỹ| y) = ∫p(ỹ|θ)p(θ|y)dθ
Advocates of Bayesian approaches argue that this presents a sensible framework
within which to be explicit about the prior assumptions of the distribution of
parameters within a model and also the structure of the model itself in multiparameter situations (Lee 2004).
Much of the controversy surrounding Bayesian
techniques has stemmed from the use of subjective priors or distributions of p(θ).
However, in many situations, the use of ‘flat’ or non-informative prior distributions
produces results identical to classical methods (Gelman et al. 2004; Lee 2004).
In toxicity risk assessments, a basic Bayesian method has been applied to calculating
confidence intervals for the fraction of species affected at specific exposures
(Aldenberg & Jaworska 2000). A similar approach has also been applied to the
determination of compliance with regulatory percentile standards (McBride & Ellis
2001). Nayak and Kundu (2001) discuss the use of Bayesian methodology to
propagate uncertainty in risk assessments and apply these methods to calculating
hazard quotients for ingestion exposure in humans. In another example Englehardt
(2004) described a Bayesian dose-response assessment for human health risk
assessment of viruses.
Uncertainty characterisation and propagation in SSD modelling
The first of the refined methods presented here involves the application of
Probability Bounds Analysis (PBA) to SSD modelling. In this method, PBA is used to
construct bounds on a SSD from commonly available Effect Concentration (EC)
estimates and their confidence intervals. In a similar manner to fitting probability
distributions to groups of scalar tolerance values, in order to construct a standard
parametric SSD, PBA can be used to combine the upper and lower confidence
intervals into a probability box that describes the uncertainty bounds on the SSD.
Specifically, the ‘histogram’ and ‘mixture’ functions in RAMAS RiskCalc enable p-
20
boxes to be constructed from groups of intervals (min-max) data. This operation is
performed empirically, without assumptions of underlying distributional shapes
(although parametric distributions may also be used), and can be used in
calculations of the (bounded) likelihood that various values are exceeded for a given
exposure distribution.
Using this PBA method, bounds on the SDD can also be constructed using other data
sources, such as the highest and lowest observed field tolerance of a species. In this
way, taxonomic uncertainty can also be incorporated into the model by including the
highest and lowest observed tolerance for the lowest level of taxonomic resolution
available. For example, if an SSD is to be constructed for the salinity tolerance of a
group of organisms that occur in a particular area of interest, except that many of
them have only been identified to Family level, bounds on the SSD could be defined
using the minimum and maximum reported tolerance values for that family. This
carries with it the consequence that the bounds would generally be much wider than
they would for a particular species within that family, encapsulating the increased
uncertainty caused by low taxonomic resolution data. As further research is
conducted and taxonomic resolution increases, the bounds on the SSD model will
generally narrow, reducing uncertainty in the SSD. Where a number of data points
are available for a family or species, the SSD model can be constructed using the
median of all the reported tolerance values for that species (i.e. the medial response)
in order to provide a less biased estimator of effects.
An example of a bounded SSD is given in Figure 1, which compares a normal
distribution fitted to log-transformed tolerance data (blue line), a non-parametric
distribution (ECDF) of the median tolerance values (black line) and a p-box of a
mixture distribution of the corresponding 95% confidence intervals (red lines) (Dixon
2005). The application of this approach is described in more detail in the associated
case study, Dixon (2007a).
21
Proportion of Species Affected
Lower
Bound
1
Upper Bound
0.5
Log-normal
Mean
Tolerance
0
0
1
2
3
4
5
6
Concentration (mg/L)
Figure 1 Example of a bounded Species Sensitivity Distribution model. The blue
line is a normal distribution fitted to log-transformed tolerance data, the black line is
an empirical CDF of the median tolerance values (medial response), and the red
lines are a p-box (mixture distribution) of the corresponding 95% confidence
intervals.
22
5
A Generalised Linear Model Approach to Predictive
Exposure-Effects Modelling
While the uncertain SSD approach described above improves on previous methods, it
still relies on a probabilistic interpretation of species’ response. In order to avoid
this assumption the second and third refined methods presented for Salt Impact
Modelling involve the construction of parametric functions to model the response of
each species’ population. Generalised Linear Models (GLM) and Generalised Linear
Mixed Models (GLMM) form the basis for these parametric risk models and are a
versatile, informative and natural approach to describing exposure effects data.
These approaches are also amenable to uncertainty propagation
using
the
techniques outlined below.
A binomially distributed logit model is the basis for the population level risk model
described here (although other formulations may be appropriate, depending on the
situation):
logit() = B + B *X
0
(3)
1
where: logit = log(/1-),  is the proportion of organisms responding and X is the
concentration of a hazardous agent. This model describes a wide range of sigmoid
shaped curves commonly used to describe toxicological responses where a binomial
outcome (i.e. survival) is concerned. In this case, a logit model is fit to tolerance data
for each species. This is then used to either consider the impact of changes on each
species separately or the individual logit model might be combined into a larger
meta-model of community response (discussed later).
In both the PBA and Bayesian techniques, uncertainty in the model is estimated
directly from the experimental data and described as distributions on the
parameters B and B . However, the methods used to estimate and propagate this
0
1
uncertainty differ. Also, dependencies between the parameters are accounted for in
the PBA approach by effectively including all possible dependencies between the
parameters, while in the Bayesian method correlation structures are estimated
directly from the data and preserved through the calculation of predictions.
Unlike previous ‘joint probability’ approaches to calculating ecotoxicological risk, the
method described here involves arithmetic operations of probabilistic terms within a
parametric model. This approach has four main advantages: First, techniques for
propagating uncertainty are more amenable to arithmetic formulation and therefore
are more easily applied within a parametric function. Second, and perhaps more
importantly, this approach models whole exposure-response distributions, rather
23
than scalar summary values, and is easily expanded to include more terms such as
other sources of uncertainty or modifying factors. Third, it provides an alternative to
the SSD approach and does not require the associated assumptions. Fourth, this
general formulation of risk model can handle a wide range of parameterisations and
underlying error-distributions, including higher-order, non-linear effects and models
with mixtures of error distributions.
Probability Bounds Analysis
In the PBA approach to uncertainty characterisation in GLMs, a logit or other model
is fit to experimental data using maximum likelihood techniques found as standard
in many statistical packages (such as SAS, S-Plus or R). Uncertainty in the fit of the
model to the data is determined using the standard error (SE) values for the
parameter estimates (i.e. 95% intervals on the mean parameter estimates are ±
1.96*SE) and these are used to describe lower and upper bounds on the parameters,
B and B . Model predictions are then computed with RiskCalc, Version 4.0 (Ferson
0
1
2002), using the logistic equation:
~
y 
1
1 exp(B0 B1X )
(4)
Or more specifically with bounds of the parameters:
~
y 
1
1 exp
([B0(MIN) , B0(MAX) ] [B1(MIN) , B1(MAX) ] *(dist(X) ~ (mean[MIN,MAX] , SD[MIN,MAX] )))
(5)
where dist(X) (the Environmental Concentration Distribution) can be any probability
distribution, scalar, interval or p-box. In this example the distribution is shown as a
p-box with bounds on its parameters, shown as minimum and maximum values for
the mean and standard deviation. A mean predicted response representing
the
results of a Monte Carlo simulation can also be calculated to enable comparisons to
be made with the results of PBA. An example of the output from this type of
PBA/GLM risk assessment is given in Figure 2.
24
1
Lab
Pr
0.5
Field
Both
Mean
0
0
1
Proportion Affected
Figure 2 Example of risk distribution produced from Probability Bounds Analysis
uncertainty propagation in a GLM risk model showing the relative contribution of
various sources of uncertainty. The y-axis gives the likelihood that an exposure level
has been exceeded, while the x-axis gives the corresponding proportion of the
population expected to be affected. The mean response is given in black the
uncertainty in the exposure distribution in blue, the uncertainty in the tolerance
relationship in red and the combined uncertainties in green.
25
Bayesian Markov Chain Monte Carlo Methods
In
the
Bayesian
approach,
model
fitting
and
generating
predictions
occur
simultaneously. Although computationally intensive, this has the advantage that
correlation structures between parameters are preserved and incorporated into the
predictions (Ades & Lu 2003). The Bayesian techniques also produce a continuous
density of likelihoods, so that the probability of observing any effect size can be
determined, not just the min-max bounds. However, this comes at the cost
of
limiting the range of dependencies that could occur. Binomial logit models can be fit
using MCMC techniques in programs such as WinBugs or OpenBugs (Spiegelhalter et
al. 2003). The algorithms applied in MCMC estimation, such as the Gipps sampling
method,
are
iterative
and
involve
sequential
sampling
from
approximate
distributions for the parameters in a model, θ , θ ,…θ , in order to find the best fitting
1
2
n
distribution for the parameters, conditional on the data. At each successive iteration,
the ‘random walk’ sampling chain is corrected such that the next sample draw
depends only on the previous value, with the aim that the chain of samples
eventually converges on the target distribution, p(θ|y). Recall that this is the posterior
probability density of a parameter given the data, conditional on the other
parameters in the model.
The use of Bayesian techniques to estimate logistic models for standard bioassay
data is discussed in Gelman et al. (2004) and Goodman (2004).
The general risk model is of the form of Equation 3:
logit (θ ) = B + B *X
i
0
1
i
where θ is the probability of an effect (y |n ), given concentration X , and logit(θ ) =
i
i
i
i
i
log(θ |1- θ ). For each replicate, the number of organisms responding y out of the total
i
i
i
n , is binomially distributed (Bin), such that:
i
y |θ  Bin(n , θ )
i
i
i
i
which gives the posterior distribution:
p(B , B |y, n, X )  p(B , B |n, X ) x p(y |B , B , n, X )
0
1
0
1
0
1
And the likelihood function:
l(B0 , B1 | y)  i
yi
(1 i )
ni yi
i.e. the likelihood l, of the parameters B and B , conditional on the data y, is given by
0
1
the product () of the binomial probabilities (Gelman et al. 2004; Lee 2004).
26
In the Bayesian MCMC fitting procedure, the parameters to be estimated are defined
as nodes within an overall model that specifies the full joint probability distribution
of all quantities. Estimation then seeks to find the joint posterior density of the
parameters, conditional on the observed quantities (y, the data), the ‘prior’
assumptions and the constraints of the model. As with the PBA approach,
predictions can be made to find the likely effects for different exposure distributions
or to determine which exposure regimes may have minimal effect on the population.
An example of the output from a Bayesian risk assessment is given in Figure 3 which
compares the risks of exposure to two different populations of a species. Further
details and examples of the application of these techniques to population level risk
assessment are given in an associated case study, Dixon (2007b).
27
1
Cumulative Probability (Exccedence)
0.9
0.8
W o o rine n
0.7
R o und
0.6
0.5
0.4
0.3
0.2
0.1
0
0
0 .1
0.2
0.3
0 .4
0 .5
0.6
0 .7
0 .8
0 .9
P ropo rtion A ffe c te d
Figure 3 Example of uncertainty risk distributions produced using Bayesian MCMC techniques. Compares the effects of a hazardous agent
(salinity) on two different populations of a species. The x-axis gives the proportion of the population affected, while the y-axis is the
probability of exceeding these values, given the exposure distributions. The spread of points describes the uncertainty in each distribution.
28
1
6
Summary and Conclusions
While risk assessment tools, such as probabilistic modelling, can provide the basis
for a more objective and transparent approach to decision making, it should be
noted that predicting the ecological changes from increasing salinity is a complex
and difficult task. The tolerance of a species to a hazardous agent may not be
constant and can vary between locations and times. The effects of a hazardous agent
on other species and ecological processes can lead to intricate indirect effects that
may not be predictable or well understood. Additionally, the limited understanding
of the behaviour of most ecosystems restricts the ability to predict their response to
disturbance.
For these reasons, risk assessment approaches are focused on using the information
available to inform decision making in the best way possible. Many models used in
risk assessment may not be correct or accurate, yet still provide adequate
comparisons of the likely outcomes of different management actions. In risk
assessments of the effects of salinity, no particular set of models or data will be
suitable for all situations and thus risk assessments must have some element of
localisation or specificity.
In refining the approach for Salt Impact Modelling, an investigation into methods for
uncertainty analysis in risk models was conducted. This found that one of the major
limitations of previous approaches has been describing uncertainty as a lack of fit of
scalar values to an assumed model, rather than the propagation of uncertainty
estimated from the original data or additional information.
Unknowable sources of incertitude, such as extrapolation from laboratory to field
scales of exposure and effect can introduce considerable uncertainty into risk
assessments. While accurate prediction of in situ effects may not be the goal of
many risk assessments, quantification of uncertainty is just as important when
comparing risks under management scenarios that involve tradeoffs and relative
risks.
Species Sensitivity Distribution based approaches to risk assessment have many
attractive properties (e.g. being able to extrapolate to community level effects) and in
the past have provided a useful framework within which to consider the
environmental risk of hazardous agents (Forbes & Calow 2002). However, these and
other current approaches used in ecotoxicological risk assessments are limited by
their underlying assumptions and involve unknown uncertainties. In particular, the
notion that the loss of a proportion of species is acceptable, without specific regard
to either the community structure, which species are expected to be lost or the effect
29
on species’ populations, obviously has limited ecological relevance (Suter et al. 2002;
van der Hoeven 2004).
An approach was described for characterising uncertainty in SSDs using Probability
Bounds Analysis to construct bounds on the model, based on Effect Concentration
(EC%) data and associated confidence intervals. This uncertain SSD method can be
easily implemented in the calculation of both the forward and inverse situations
(predicting effects and setting exposure standards) using currently available
software, although it should be noted that some of the assumptions involved in this
approach may be incorrect (as discussed above).
In comparison to the SSD approach, a more detailed method for risk modelling was
described that focused on each species’ response. These species level responses can
then be built into a larger model of community effects that incorporates the whole
response and uncertainties of each population and may also include interactions and
complexities between species. Single species risk assessments are also important in
the management of threatened or target species and are focused on predicting
effects on populations, rather than communities (although the methods used may be
expanded to apply to communities). The application of Bayesian MCMC and PBA
techniques to these models enables quantification of distributions of risk and
associated uncertainty, as well as direct quantitative comparisons between species,
situations or scenarios.
In extending these methods to multi-species or community level risk assessments, a
number of approaches are possible depending on the information available, the
outcomes required and the assumptions that can be justified. If the aim of a risk
assessment is to minimise effects on a particular species and/or those it depends on
then separate modelling of the risks to the individual species would provide
information relevant to the management of the situation. If a larger model of the
expected community sensitivities to exposure was required, a hierarchical approach
that formed a (community) distribution from underlying (species) population
tolerance models could be undertaken (e.g. Dixon & Scroggie 2007). Bayesian MCMC
techniques are highly amenable to hierarchical modelling and it would be possible to
construct a model of species within a community that might be affected, where each
underlying population has a cut off for inclusion (e.g. 95% chance of 50% of the
population affected). Mixture models of different underlying distributions might also
be suitable for community analysis and treatment of overdispersion in larger
datasets. Finally, perhaps the preferred approach would be to include the ‘uncertain’
risk distributions as a parameter in a larger model of exposure and ecological effect,
where the interdependencies between species and their populations is considered.
30
The methods described in this report enable quantification of risk distributions and
associated uncertainties from the original data, propagation of this uncertainty
through the risk models and quantitative comparisons between risk distributions. In
developing these alternative risk models few assumptions were required about the
underlying distribution of parameters and those assumptions that were made are
easily
tested.
The
approaches
developed
are
also
amenable
to
alternative
distributional structures, re-parameterisation, inclusion of other estimates of
uncertainty and do not assume independence between the parameters. Although the
examples provided are focused on single species populations, the approaches used
can be extended to apply to communities or to feed into larger models of ecosystem
behaviour and risk.
The integration of more objective and defensible scientific approaches into the
decision making process helps achieve regional and State conservation goals by
providing a rigorous and transparent method of assessing risks and prioritising
resources. This ensures efficiency in government investment, facilitates a landscape
level approach to land and water management, and helps to maximise the potential
for biodiversity protection.
31
Appendix 1: Basic Methods for Salt Impact Modelling
This section describes a number of basic methods for predicting the effects of
salinity on biodiversity. The first approach described is to use tolerance or look-up
tables to compare between measures of sensitivity for different species of interest.
Although simplistic, this approach has the advantage that it does not make the
statistical assumptions of SSD based techniques.
The other two methods described here are for SSD approaches to effects modelling.
The first uses a program called BurrliOZ, which was developed for the ANZECC
guidelines, while the second applies more standard approaches using commonly
available software packages.
A.1 - Lookup Tables and Plots
This basic approach to describing tolerance data provides a way to visualise the
information and indicate which species are likely to be affected at various levels of
salinity. The method involves constructing a lookup table or ranked-plot of the
tolerances of the different species of concern. The limitations of this approach are
that it does not enable extrapolations or estimates of confidence/uncertainty to be
made, and thus the results should be considered cautiously.
As with the other methods described here, the initial step is to define the subject
area, the species of interest contained within it, and to collate the relevant data on
sensitivity. This data should then be entered into a spreadsheet, such as in Microsoft
Excel, with species name in one column and the corresponding tolerance value/s in
the next. The next stage is to ‘sort’ the data in ascending order by tolerance. This is
done by highlighting both columns, clicking ‘tools’ and then ‘sort’ tabs, selecting the
relevant column for the tolerance values and checking the ‘ascending’ box. A plot
can then be then created (i.e. in Excel select the graph wizard and choose column
graph) with species on the x-axis and tolerance on the y-axis. The salinity level of
concern is then read across the graph from the y-axis. Any bars that fall below this
line represent species that are likely to be affected at this level and subsequent
attention can then be given to their conservation status or other specialist concerns.
Visual comparisons may also be made between groups within the plot, such as
between native and introduced species.
This method can only provide a basic indication of likely effects of salinity on
species of concern, and does not provide estimates of confidence or uncertainty.
32
Careful consideration should be given to the relevance and inferential consequences
of using different data sources. An example of the tolerance plot approach is given
in Figure A1.
33
70000
Salinity (us/cm)
60000
50000
Native
Introduced
40000
30000
20000
10000
0
Figure A1 Compares the salt tolerance of native and introduced fishes. It can be seen here that at a level of 20,000 EC (red line) at least 5 native
species would be expected to be affected, a number of which are listed and threatened or endangered. It can also been seen that this level is close to
the tolerance of three other native species. The effect on introduced species can also be compared.
34
A.2 – BurrliOZ Software Example
BurrliOZ is a free program that performs parametric modelling of SSDs. Developed
by Shao (2000), BurrliOZ fits a Burr Type III distribution to sets of tolerance data and
computes confidence intervals around a specified percentile. It has the advantage of
being easy to use and has been used in the ANZECC water quality guidelines (Fox
1999; ANZECC & ARMCANZ 2000; Shao 2000). However, the limitations of this
method are that goodness of fit and model checking statistics are limited and the
distribution itself is not characterised in a way that can be used in risk assessments.
The program is available from CSIRO4 and includes instructional documentation.
Once the species list and tolerance data have been obtained (from sources such as
the Land and Water database), the individual tolerance points are entered into the
program,
along
with
the
effect-percentile
of
interest.
Resultant
hazardous
concentrations and confidence intervals are then computed.
This is an easy way to perform parametric SSD modelling and compute confidence
intervals. However, care should be given to the influence of the data used and the
lack of statistical checking diagnostics. An example of this approach is given in
Figure A2.
4
www.cmis.csiro.au/Envir/burrlioz/
35
Figure A2 Display of typical output from the BurrliOZ program. The solid blue line represents
the Burr Type III distribution that is used to calculate percentiles. A log-log (orange line) and
log-normal (green line) fit are shown for comparison. Red dots represent data points and the
purple line indicates the percentile of choice.
36
A.3 - Parametric SSD Modelling
While a number of statistical packages can fit parametric probability distributions to
data, only a few estimate confidence intervals for percentiles. In particular, recent
versions of Minitab (i.e. 13 and 14) have included this function as standard. The
more specific and powerful packages such as SAS and S-plus are also able to produce
confidence intervals, although they are more expensive and difficult to use.
Many of the published studies involving SSDs assume that data is log-normally
distributed and relatively few provide tests of normality or comparisons with other
parametric models. A few studies have compared the fit of different parametric
models, generally finding that the assumption of (log)normality was not supported
in over half of cases and that log-logistic was a better fit (Newman et al. 2000). There
seems to be little justification for the a priori assumption that species sensitivities
are normally or log-normally distributed (Newman et al. 2000). However, statistical
tests of normality (such as Shapiro-Wilks) should provide some basis for the choice
of distribution. Comparison between the fit of different distributions should be
made using a goodness of fit statistic such as Anderson-Darling, KolmogorovSmirnov or the various likelihood information criterion (such as
Akaike’s
Information Criterion) (Buckland et al. 1997; Burnham & Anderson 2001), and may
be aided by examining plots of the fitted distribution. The choice of distribution is a
much debated issue and one particular model will not be advocated over another
here, except to note that in published studies the log-logistic has been found to be a
better fit more often than other common models (Newman et al. 2000; Posthuma et
al. 2002a).
To fit the model, tolerance or sensitivity data sets are entered into the packages
spread-sheet and the distributions fit (Minitab has a menu command for this
whereas SAS and other programs require command line instructions). Minitab
automatically produces a plot and table of confidence intervals and an Anderson
Darling statistic. These percentiles, confidence intervals and plots can then be used
to assess the fit of the model, choose better fitting models and make inferences
about the likely effects of salinity on species or groups of species. A Minitab example
is given in Figure A3.
37
Figure A3 Compares the fit of a log-logistic (parametric) distribution to data for two different
lakes (shown in red and black). The straight, middle, lines depict the mean response, with the
outer lines representing the upper and lower 95% confidence intervals (a table of these is
also produced by the software). The solid dots are the data points. The parameters for each
distribution are given as ‘loc’ (location) and ‘scale’ in the key, as are the number of data points
(N), Anderson-Darling statistic (AD) and p-value.
38
References
Ades A. E. & Lu G. (2003) Correlations between parameters in risk models: estimation
and propagation of uncertainty by Markov Chain Monte Carlo. Risk Analysis 23:
1165-1172.
Aldenberg T. & Jaworska J. S. (2000) Uncertainty of the hazardous concentration and
fraction affected for normal species sensitivity distributions. Ecotoxicology and
Environmental Safety 46: 1-18.
Aldenberg T., Jaworska J. S. & Traas T. P. (2002) Normal species sensitivity
distributions and probabilistic ecological risk assessment. In: Species Sensitivity
Distributions in Ecotoxicology (eds. L. Posthuma, G. W. Suter II & T. P. Traas) pp. 49102. Lewis Publishers, Boca Raton.
Aldenberg T. & Slob W. (1993) Confidence limits for hazardous concentrations based
on logistically distributed NOEC toxicity data. Ecotoxicology and Environmental
Safety 25: 48-63.
Andrews C. J., Hassenzahl D. M. & Johnson B. B. (2004) Accommodating uncertainty
in comparative risk. Risk Analysis 24: 1323-1335.
ANZECC & ARMCANZ (2000) National Water Quality Management Strategy Australian and New Zealand Guidelines for Fresh and Marine Water Quality.
Australian and New Zealand Environment and Conservation Council (ANZECC) &
Agriculture and Resource Management Council of Australia and New Zealand
(ARMCANZ), Canberra, Australia.
Bailey
P.,
Boon
P.
&
Morris
K.
(2002)
Salt
Sensitivity
Database.
www.rivers.gov.au/Tools_and_Techniques/Salt_Sensitivity_Database/index.aspx.
Land and Water Austalia.
Bartell S. M., Pastorok R. A., Akçakaya H. R., Regan H., Ferson S. & Mackay C. (2003)
Realism and relevance of ecological models used in chemical risk assessment.
Human and Ecological Risk Assessment 9: 907-938.
Bates S. C., Cullen A. & Raftery A. E. (2003) Bayesian uncertainty assessment in
multicompartment deterministic simulation models for environmental risk
assessment. Environmetrics 14: 35-371.
Berleant D. & Zhang J. (2004) Using pearson correlation to improve envelopes around
the distributions of functions. Reliable Computing 10: 139-161.
Bixio D., Parmentier G., Rousseau D., Verdonck F. A. M., Meirlaen J., Vanrolleghem P.
A. & Thoeye C. (2002) A quantitative risk analysis tool for design/simulation of
wastewater treatment plants. Water Science and Technology 46: 301-307.
Brattin W. J., Barry T. M. & Chiu N. (1996) Monte Carlo modeling with uncertain
probability density functions. Human and Ecological Risk Assessment 2: 820-840.
Breton R. L., Teed R. S. & Mooreb D. R. J. (2003) An ecological risk assessment of
phenol in the aquatic environment. Human and Ecological Risk Assessment 9: 549568.
Buckland S. T., Burnham K. P. & Augustin N. H. (1997) Model selection: An
intergrated part of inference. Biometrics 53: 603-618.
Burgman M. A. (2005) Risks and Decisions for Conservation and Environmental
Management. Cambridge University Press, London, UK.
Burmaster D. E. & Wilson J. C. (1998) Risk assessment for chemicals in the
environment. In: The Encyclopedia of Biostatistics, First Edition. John Wiley & Sons,
Inc.
Burnham K. P. & Anderson D. R. (2001) Kullback–Leibler information as a basis for
strong inference in ecological studies. Wildlife Research 28: 111-119.
39
Cairns Jr J. (1999) Absence of certainty is not synonymous with absence of risk.
Environmental Health Perspectives 107: A56-A66.
Cant B., James K. R. & Ryan T. (2003) Salt impact model for strategic framework and
information to indicate biodiversity thresholds to salinity pp. 82. Arthur Rylah
Institute, Department of Sustainability & Environment, Heidelberg, Victoria.
Cardwell R. D., Parkhurst B. R., Warren-Hicks W. & Volosin J. S. (1993) Aquatic
ecological risk. Water Environment and Technology 4: 47–51.
Chapman P. M., Caldwell R. S. & Chapman P. F. (1996) A warning: NOECs are
inappropriate for regulatory use. Environmental Toxicology and Chemistry 15: 77-79.
Chapman P. M., Fairbrother A. & Brown D. (1998) A critical evaluation of safety
(uncertainty) factors for ecological risk assessment. Environmental Toxicology and
Chemistry 17: 99-108.
Congdon P. (2001) Bayesian Statistical Modelling. John Wiley & Sons., Chichester, UK.
Congdon P. (2003) Applied Bayesian Modelling. John Wiley & Sons, Ltd., Chichester,
UK.
Crane M., Whitehouse P., Comber S., Watts C., Giddings J. M., Moore D. R. J. & Grist E.
(2003) Evaluation of probabilistic risk assessment of pesticides in the UK:
chlorpyrifos use on top fruit. Pest Management Science 59: 512-526.
de Bruijn J. H. M. & Hof M. (1997) How to measure no effect. Part IV: How acceptable
is the ECx from an environmental policy point of view? Environmetrics 8: 263-267.
de Zwart D. (2003) Ecological effects of pesticide use in the Netherlands. RIVM
Report 500002003/2003, Bilthoven.
Dixon W. J. (2005) Uncertainty in aquatic toxicological exposure-effect models: the
toxicity of 2,4-Dichlorophenoxyacetic acid and 4-Chlorophenol to Daphnia carinata.
PhD Thesis. RMIT University, Melbourne, Australia.
Dixon W. J. (2007a) The use of probability bounds analysis for characterising and
propagating uncertainty in species sensitivity distributions. Arthur Rylah Institute
for Environmental Research, Technical Report Series No. 163. Department of
Sustainability and Environment, Heidelberg, Victoria.
Dixon W. J. (2007b) Uncertainty propagation in population level salinity risk models.
Arthur Rylah Institute for Environmental Research, Technical Report Series No. 164.
Department of Sustainability and Environment, Heidelberg, Victoria.
Dixon W. J. & Scroggie M. P. (2007) A Bayesian approach to multi-species risk
assessments using whole exposure-response distributions. Arthur Rylah Institute for
Environmental Research, Technical Report Series No. 165. Department of
Sustainability and Environment, Heidelberg, Victoria.
Englehardt J. D. (2004) Predictive Bayesian dose-response assessment for appraising
absolute health risk from available information. Human and Ecological Risk
Assessment 10: 69 - 78.
Ferson S. (1996) What Monte Carlo methods cannot do. Human and Ecological Risk
Assessment 2: 990-1007.
Ferson S. (2002) RAMAS Risk Calc 4.0 Software: Risk Assessment with Uncertain
Numbers. Lewis Publishers,, Boca Raton, Florida.
Ferson S. & Ginzburg L. R. (1996) Different methods are needed to propagate
ignorance and variability. Reliability Engineering & System Safety 54: 133-144.
Ferson S., Ginzburg L. R. & Goldstein R. A. (1996) Inferring ecological risk from
toxicity bioassays. Water, Air and Soil Pollution 90: 71-82.
Ferson S. & Hajagos J. G. (2004) Arithmetic with uncertain numbers: rigorous and
(often) best possible answers. Reliability Engineering & System Safety 85: 135-152.
40
Ferson S., Kreinovich V., Ginzburg L. R., Myers D. S. & Sentz K. (2003) Constructing
Probability Boxes and Dempster-Shafer Structures. Report Number: SAND2002-4015.
Sandia National Laboratories, Albuquerque, New Mexico.
Finkel A. M. (1995) Toward less misleading comparisons of uncertain risks: the
example of Aflatoxin and Alar. Environmental Health Perspectives 103: 376-385.
Fogarty M. J., Mayo R. K., O'Brien L., Serchuk F. M. & Rosenberg A. A. (1996) Assessing
uncertainty and risk in exploited marine populations. Reliability Engineering &
System Safety 54: 183-195.
Forbes V. E. & Calow P. (2002) Species sensitivity distributions revisited: a critical
appraisal. Human and Ecological Risk Assessment 8: 473-492.
Fox D. R. (1999) Setting water quality guidelines - A statistician's perspective. SETAC
News May: 17-18.
Frey H. C. & Burmaster D. E. (1999) Methods for characterizing variability and
uncertainty: comparison of bootstrap simulation and likelihood-based approaches.
Risk Analysis 19: 109-130.
Gelman A., Carlin J. B., Stern H. S. & Rubin D. B. (2004) Bayesian Data Analysis.
Second Edition. Chapman & Hall/CRC, Boca Ranton, Florida, USA.
Godsill S. J. (2001) On the relationship between markov chain monte carlo methods
for model uncertainty. Journal of Computational and Graphical Statistics 10: 1-19.
Hart B. T., Maher B. & Lawrence I. (1999) New generation water quality guidelines for
ecosystem protection. Freshwater Biology 41: 347-359.
Hattis D. (1990) Three candidate "laws" of uncertainty analysis. Risk Analysis 10: 11.
Hoekstra J. A. & van Ewijk P. H. (1993) Alternatives for the no-observed-effect level.
Environmental Toxicology and Chemistry 12: 187-194.
Jager T., den Hollander H. A., van der Poel P., Rikken M. G. J. & Vermeire T. (2001)
Probabilistic environmental risk assessment for Dibutylphthalate (DBP). Human and
Ecological Risk Assessment 7: 1681-1697.
Jagoe R. H. & Newman M. C. (1997) Bootstrap estimation of community NOEC values.
Ecotoxicology 6: 293-306.
Kefford B. J., Nugegoda D., Metzeling L. & Fields E. J. (2006) Validating species
sensitivity distributions using salinity tolerance of riverine macroinvertebrates in the
southern Murray-Darling Basin (Victoria, Australia). Canadian Journal of Fisheries
and Aquatic Sciences 63: 1865-1877.
Kefford B. J., Palmer C. G., Warne M. S. & Nugegoda D. T. (2005) What is meant by
"95% of species"? An argument for the inclusion of rapid tolerance testing. Human
and Ecological Risk Assessment 11: 1025-1046.
Landis W. G. (2003) The frontiers in ecological risk assessment at expanding spatial
and temporal scales. Human and Ecological Risk Assessment 9: 1415 - 1424.
Laskowski R. (1995) Some good reasons to ban the use of NOEC, LOEC and related
concepts in ecology. Oikos 73: 140-144.
Lee P. M. (2004) Bayesian Statistics. Third Edition. Hodder Arnold, London, UK.
Linkov I., von Stackelberg K. E., Burmistrov D. & Bridges T. S. (2001) Uncertainty and
variability in risk from trophic transfer of contaminants in dredged sediments. The
Science of the total environment 274: 255-269.
Manly B. F. J. (1997) Randomization, Bootstrap and Monte Carlo Methods in Biology,
2nd Edition. Chapman & Hall, London, UK.
McBride G. B. & Ellis J. C. (2001) Confidence of compliance: a bayesian approach for
percentiles standards. Water Research 35: 1117-1124.
41
Merrick J. R. W., van Dorp J. R. & Dinesh V. (2005) Assessing uncertainty in
simulation-based maritime risk assessment. Risk Analysis 25: 731-743.
Moore D. R. J., Sample B. E., Suter G. W., Parkhurst B. R. & Teed R. S. (1999) A
probabilistic risk assessment of the effects of methylmercury and PCBs on mink and
kingfishers along east fork Poplar Creek, Oak Ridge, Tennessee, USA. Environmental
Toxicology and Chemistry 18: 2941-2953.
Morgan M. G. & Henrion M. eds. (1990) Uncertainty: A Guide to Dealing with
Uncertainty in Quantitative Risk and Policy Analysis. Cambridge University Press,
New York, NY USA.
Nayak T. K. & Kundu S. (2001) Calculating and describing uncertainty in risk
assessment: the Bayesian approach. Human and Ecological Risk Assessment 7: 307328.
Newman M. C., Ownby D. R., Mezin L. C. A., Powell D. C., Chiristensen T. R. L., Lerberg
S. B. & Anderson B. A. (2000) Applying species-sensitivity distributions in ecological
risk assessment: assumptions of distribution type and sufficient numbers of species.
Environ. Toxicol. Chem. 19: 508-515.
Newman M. C., Ownby D. R., Mézin L. C. A., Powell D. C., Christensen T. R. L., Lerberg
S. B., Anderson B.-A. & Padma T. V. (2002) Species sensitivity distributions in
ecological risk assessment: distributional assumptions, alternate bootstrap
techniques, and estimation of adequate number of species. In: Species Sensitivity
Distributions in Ecotoxicology (eds. L. Posthuma, G. W. Suter II & T. P. Traas) pp. 109132. Lewis Publishers, Boca Raton.
NRC (1993) Issues in Risk Assessment: a Paradigm for Ecological Risk Assessment.
NRC (National Research Council), Committee on Risk Assessment Methodology,
Board on Environmental Studies and Toxicology, Commission on Life Sciences, and
National Research Council, National Academy Press, Washington D.C. USA.
NRC (1994) Science and Judgment in Risk Assessment. NRC (National Research
Council), Committee on Risk Assessment of Hazardous Air Pollutants, Board on
Environmental Studies and Toxicology, Commission on Life Sciences, and National
Research Council, National Academy Press, Washington D.C. USA.
Pastorok R. A., Akçakaya H. R., Regan H., Ferson S. & Bartell S. M. (2003) Role of
ecological modeling in risk assessment. Human and Ecological Risk Assessment 9:
939-972.
Pastorok R. A., Bartell, S.M., Ferson, S., Ginzburg, L.R. ed. (2002) Ecological Modeling
in Risk Assessment: Chemical Effects on Populations, Ecosystems and Landscapes.
Lewis Publishers, Boca Raton, FL.
Posthuma L., Suter II G. W. & Traas T. P. eds. (2002a) Species Sensitivity Distributions
in Ecotoxicology. Lewis Publishers, Boca Raton.
Posthuma L., Traas T. P., de Zwart D. & Suter II G. W. (2002b) Conceptual and
technical outlook on species sensitivity distributions. In: Species Sensitivity
Distributions in Ecotoxicology (eds. L. Posthuma, G. W. Suter II & T. P. Traas) pp. 475508. Lewis Publishers, Boca Raton.
Rand G. M. ed. (1995) Fundamentals of Aquatic Toxicology, Second Edition. Taylor and
Francis, Washington, D.C. USA.
Regan H. M., Colyvan M. & Burgman M. A. (2002a) A taxonomy and treatment of
uncertainty for ecology and conservation biology. Ecological Applications 12: 618628.
Regan H. M., Ferson S. & Berleant D. (2004) Equivalence of methods for uncertainty
propagation of real-valued random variables. International Journal of Approximate
Reasoning 36: 1-30.
42
Regan H. M., Hope B. K. & Ferson S. (2002b) Analysis and portrayal of uncertainty in a
food-web exposure model. Human and Ecological Risk Assessment 8: 1757 - 1777.
Regan H. M., Sample B. E. & Ferson S. (2002c) Comparison of deterministic and
probabilistic calculation of ecological soil screening levels. Environmental Toxicology
and Chemistry 21: 882-890.
Shao Q. (2000) Estimation for hazardous concentrations based on NOEC toxicity
data: an alternative approach. Environmetrics 11: 583-595.
Soloman K. R., Giesy J. & Jones P. (2000) Probabilistic risk assessment of
agrochemicals in the environment. Crop Protection 19: 649-655.
Soloman K. R. & Sibley P. (2002) New concepts in ecological risk assessment: where
do we go from here? Marine Pollution Bulletin 44: 279-285.
Spiegelhalter D. J., Thomas A., Best N. G. & Gilks W. R. (2003) BUGS: Bayesian
inference Using Gibbs Sampling, Version 1.4. MRC Biostatistics Unit, Cambridge.
Suter G. W. (1998) Comments on the interpretation of distributions in “Overview of
recent developments in ecological risk assessment.” Risk Analysis 18: 3-4.
Suter G. W., II, Efroymson R. A., Sample B. E. & Jones D. S. (2000) Ecological Risk
Assessment for Contaminated Sites. Lewis Pub., Boca Raton, USA.
Suter G. W., II, Traas T. P. & Posthuma L. (2002) Issues and practices in the derivation
and use of species sensitivity distributions. In: Species Sensitivity Distributions in
Ecotoxicology (eds. L. Posthuma, G. W. Suter II & T. P. Traas) pp. 437-474. Lewis
Publishers, Boca Raton.
Thompson K. M. (1999) Developing univariate distributions from data for risk
analysis. Human and Ecological Risk Assessment 5: 755-783.
Thompson K. M. (2002) Variability and uncertainty meet risk management and risk
communication. Risk Analysis 22: 647-654.
Thompson K. M. & Graham J. D. (1996) Going beyond the single number: using
probabilistic risk assessment to improve risk management. Human and Ecological
Risk Assessment 2: 1008-1034.
van der Hoeven N. (1997) How to measure no effect. Part III: Statistical aspects of
NOEC, ECx and NEC estimates. Environmetrics 8: 255-261.
van der Hoeven N. (2001) Estimating the 5-percentile of the species sensitivity
distribution without any assumptions about the distribution. Ecotoxicology 10: 25-34.
van der Hoeven N. (2004) Current Issues in Statistics and Models for Ecotoxicological
Risk Assessment. Acta Biotheoretica 52: 201–217.
Van Sprang P. A., Verdonck F. A. M., Vanrolleghem P. A., Vangheluwe M. L. & Janssen
C. R. (2004) Probabilistic environmental risk assessment of Zinc in Dutch surface
water. Environmental Toxicology and Chemistry 23: 2993-3002.
van Straalen N. M. (1990) New methodologies for estimating the ecological risk of
chemicals in the environment. In: Proceedings of the 6th Congress of the International
Association of Engineering Geology (ed. D. G. Price) pp. pp. 165–173. A.A. Balkema,
Rotterdam, the Netherlands.
van Straalen N. M. (2002a) Theory of ecological risk assessment based on species
sensitivity distributions. In: Species Sensitivity Distributions in Ecotoxicology (eds. L.
Posthuma, G. W. Suter II & T. P. Traas) pp. 37-48. Lewis Publishers, Boca Raton.
van Straalen N. M. (2002b) Threshold models for species sensitivity distributions
applied to aquatic risk assessment for zinc. Environmental Toxicology and
Pharmacology 11: 167-172.
van Straalen N. M. & Denneman C. A. J. (1989) Ecotoxicological evaluation of soil
quality criteria. Ecotoxicology and Environmental Safety 18: 241-251.
43
van Straalen N. M. & van Leeuwen C. J. (2002) European history of species sensitivity
distributions. In: Species Sensitivity Distributions in Ecotoxicology (eds. L. Posthuma,
G. W. Suter II & T. P. Traas) pp. 19-34. Lewis Publishers, Boca Raton.
Verdonck F., Jaworska J. & Vanrolleghem P. A. (2001a) Comments on species
sensitivity distributions sample size determination based on Newman et al. (2001).
SETAC Globe, Learned Discourse 2: 22-24.
Verdonck F. A. M., Aldenberg T., Jaworska J. S. & Vanrolleghem P. A. (2003)
Limitations of current risk characterization methods in probablistic environmental
risk assessment. Environmental Toxicology and Chemistry 22: 2209-2213.
Verdonck F. A. M., Jaworska J., Janssen C. R. & Vanrolleghem P. A. (2002)
Probabilistic ecological risk assessment framework for chemical substances. In:
Integrated Assessment and Decision Support, Proceedings of the First Biennial Meeting
of the International Environmental Modelling and Software Society (IEMSS), Volume 1
(eds. A. E. Rizzoli & A. J. Jakeman) pp. 144- 149. IEMSS.
Verdonck F. A. M., Jaworska J., Thas O. & Vanrolleghem P. A. (2001b) Determining
environmental standards using bootstrapping, Bayesian and maximum likelihood
techniques: a comparative study. Analytica Chimica Acta 446: 429-438.
Warren-Hicks W. J., Parkhurst B. R. & Butcher J. B. (2002) Methodology for Aquatic
Ecological Risk Assessment. In: Species Sensitivity Distributions in Ecotoxicology (eds.
L. Posthuma, G. W. Suter II & T. P. Traas) pp. 345-382. Lewis Publishers, Boca Raton.
Webb J. A. & Hart B. T. (2004) Environmental Risks from Salinity Increases in the
Goulburn-Broken Catchment - Draft. Water Studies Centre, Monash University,
Melbourne, Australia.
Williamson R. C. & Downs T. (1990) Probabilistic arithmetic I: numerical methods for
calculating convolutions and dependency bounds. International Journal of
Approximate Reasoning 4: 89–158.
Winkler R. L. (1996) Uncertainty in probabilistic risk assessment. Reliability
Engineering & System Safety 54: 127-132.
Yager R. R. (1986) Arithmetic and other operations on Dempster-Shafer structures.
International Journal of Man-Machine Studies 25: 357-366.
Yanagawa T. & Kikuchi Y. (2001) Statistical issues on the determination of the noobserved-adverse-effect levels in toxicology. Environmetrics 12: 319-325.
44