Scholarly article on topic 'Forensic Assessment of Metal Contaminated Rivers in the 21st Century Using Geochemical and Isotopic Tracers'

Forensic Assessment of Metal Contaminated Rivers in the 21st Century Using Geochemical and Isotopic Tracers Academic research paper on "Earth and related environmental sciences"

Share paper
Academic journal

Academic research paper on topic "Forensic Assessment of Metal Contaminated Rivers in the 21st Century Using Geochemical and Isotopic Tracers"

Minerals 2013, 3, 192-246; doi:10.3390/min3020192



ISSN 2075-163X


Forensic Assessment of Metal Contaminated Rivers in the 21st Century Using Geochemical and Isotopic Tracers

Jerry R. Miller

Department of Geosciences and Natural Resources, Western Carolina University, Cullowhee, NC 28723, USA; E-Mail:; Tel.: +1-828-227-2269; Fax:+1-828-227-7647

Received: 1 March 2013; in revised form: 25 April 2013 / Accepted: 26 April 2013 / Published: 16 May 2013

Abstract: Documentation of contaminant source and dispersal pathways in riverine environments is essential to mitigate the potentially harmful effects of contaminants on human and ecosystem health, and is required from a legal perspective (particularly where the polluter pays principle is in effect) in assessing site liability. Where multiple natural and/or anthropogenic sources exist, identification of contaminant provenance has proven problematic, and estimated contributions from a specific source are often the subject of judicial debate. The past, current, and future use of geochemical and isotopic tracers in environmental forensic investigations of contaminant provenance, transport, and fate are analyzed herein for sediment-associated trace metals in riverine environments, particularly trace metals derived from mining and refining operations. The utilized methods have evolved significantly over the past four decades. Of primary significance has been the growing integration of geomorphic and stratigraphic techniques with the use of an increasing number of geochemical tracers including stable isotopes. The isotopes of Pb have been particularly well studied, and have been applied to a wide range of environmental media. Advances in analytical chemistry since the early 1990s have allowed for the precise characterization of other non-traditional stable isotopic systems within geological materials. The potential for using these non-traditional isotopes as tracers in river systems has yet to be adequately explored, but a number of these isotopes (e.g., Cd, Cu, Cr, Hg, Sb, and Zn) show considerable promise. Moreover, some of these isotopes (e.g., those of Cu, Cr, and Hg) may provide important insights into biogeochemical cycling processes within aquatic environments. This review suggests that future environmental forensic investigations will be characterized by an interdisciplinary approach that combines the use of multiple geochemical tracers with detailed stratigraphic, geomorphic, and hydrologic data, thereby yielding results that are likely to withstand the scrutiny of judicial review.

Keywords: environmental forensics; geochemical fingerprinting; isotopic tracers; non-traditional isotopes

1. Introduction

Environmental forensics has its roots in forensic geology, a field that can be defined as the use of rocks, sediments, soils, and other geological materials in legal investigations that are likely to end up in the court of law. Forensic geology first gained notoriety in the writings of Sir Arthur Conan Doyle between 1887 and 1893 when the fictional character Sherlock Holmes effectively analyzed geological substances to help solve complex criminal cases [1]. Most of the general principles involved in the present-day application of forensic geology were laid out, however, in the early 1900s when scientists began to use geological and geochemical principles in the investigation of serious crimes. One of the early pioneers in the field was Georg Popp, a German chemist. In 1908, Popp examined three layers of soil adhering to the shoe of a suspect in the Margarethe Filbert homicide [1,2]. Popp found that the composition of the soils did not match those within the agricultural fields where the suspect said he had been. However, one of the layers on the shoe contained red sandstone fragments and other particles which were associated with soils from where the body had been found, and a second layer (containing brick, coal dust, cement and other materials) were comparable to soils from outside a castle where the suspect's gun and clothing had been discovered. Thus, Popp was able to provide evidence that supported the theorized sequence of events associated with the crime [2].

Since the early work of Popp, forensic geology has given rise to a range of forensic geosciences, including environmental forensics, in which a wide range of scientific data from such fields as mineralogy, geochemistry, sedimentology, stratigraphy, geomorphology, soil science, hydrology, and geophysics are brought together to address an environmental issue associated with civil and criminal legal investigations. In essence, multiple forms of geoscientific data must be integrated such that the resulting conclusions can withstand the scrutiny associated with the judicial process. The majority of these cases pertain to the contamination of air, water, sediments, or soils by anthropogenic activities. Perhaps one of the most frequently addressed problems in environmental forensics within riverine environments is to define the source(s) and relative contributions of contaminants to a river from human activities in general, and mining, milling, and refining operations in particular. From a legal perspective, source identification is a critical component of many, if not most, site remediation programs where the polluter pays principle is in effect and the primary responsible parties must be accurately identified. It is not uncommon for millions to tens of millions of dollars to be at stake in terms of remedial costs [3,4]. From a practical perspective, identifying contaminant source(s) is essential to site mitigation, as well as human and environmental health, because remediation strategies are likely to fail if the source(s) of contaminant influx are not identified and addressed.

Although a variety of contaminants can be released from mining and milling operations to river systems, it is fair to say that toxic trace metals and metalloids (e.g., As, Cu, Cd, Hg, Pb, Sb, and Zn) are the most common [5]. Trace metals possess a high affinity for particulate matter in aquatic environments with typically observed (unaltered) pH and Eh conditions, and numerous studies have

shown that trace metals released from a point source of contamination predominantly are dispersed as elements bound to (or sorbed to) riverine sediments [6,7]. Acknowledgement of the importance of sediments as the predominant transporting media, and as a significant trace metal sink along the drainage system, has led investigators to focus in recent years on riverine sediments for the identification of contaminant sources as well as the elucidation of the rates and processes of contaminant transport [4,8].

Within riverine environments, significant portions (10% to more than 60%) of the contaminated sediments released from a point source are deposited and stored along the river valley in channel bed, floodplain, or terrace deposits where they may pose a threat to human or ecosystem health. Environmental forensics is increasingly relying on the use of geochemical fingerprinting and tracing methods to identify the sources of contaminants to these deposits, and to assess the dispersal pathways of sediment-bound contaminants in riverine environments. Following the definitions put forth by Foster and Lees [9] and Hoefs [10], a fingerprint is defined here for sediments in river systems as a unique set of physical or chemical properties of the source material that can be used to distinguish it from other sediments within the basin and whose traits are preserved following transport such that its main characteristics are still recognizable. In Popp's case, the red sandstone fragments served as an effective fingerprint for one of the soil layers found on the suspect's shoe. More recent forensic investigations have utilized a wide range of properties for fingerprinting purposes, including sediment color [11,12], grain size [13-17], particle mineralogy [18-20], mineral magnetics [21,22], surface morphology and composition [23-26], and chemical composition as defined by bulk elemental concentrations, concentration ratios, chemical species, and isotopes [15,27-29].

The term tracer is often used synonymously with fingerprint. A tracer, however, is more strictly defined in this review as a parameter which not only exhibits the characteristics of a fingerprint in that it allows a geological material to be distinguished from other materials in the basin, but can be used to track the movement and cycling of the material through the drainage system. Tracers have been used for several different purposes in riverine environments, including: (1) the identification of contaminant provenance for sediment sampled at a given location; (2) an analysis of the contaminant transport pathways and rates; and (3) an assessment of the biogeochemical cycling of a substance within the aquatic environment. As used for the first two purposes, tracers are often most effective when they exhibit conservative behavior in which its identifiable characteristics are maintained as the sediment is cycled through the riverine environment. In the latter instance, alterations in the nature of the tracer is utilized to assess a specific biogeochemical cycling process, such as microbial or photoreduction. Thus, the nature of the alterations must be well understood.

The use of fingerprints/tracers for contaminant source and transport assessment is complicated by the fact that the origin of the geological materials that we are assessing may have been altered from that found at the source. These alterations involve two primarily processes. The first refers to a set of processes through which certain particles are separated from the total (bulk) material found at the source on the basis of particle size, shape, density, etc.; thus, the sample under investigation represents only part, or a fraction, of the original material. Fine-grained clay minerals, for instance, may be separated from the original source materials by hydraulic sorting processes during transport, and subsequently deposited in areas devoid of coarse-particles [4]. Isotopic fractionation may also

occur during which physical or chemical processes alter the relative abundance of the source material's isotopes.

The second complicating process is mixing in which sediment weathered and eroded from all of the bedrock and alluvial type areas are co-mingled with contaminated sediment from the source(s) during transport and deposition. The sediments, then, not only reflect the lithology, mineralogy, and geochemistry from the various bedrock types that underlie the basin, but contaminants from all of the sources from which they were introduced to the channel along its course. Contaminant dilution, either through the spreading of the contaminated particles over a larger area, or their amalgamation with "clean" sediments, often accompanies the mixing process. In either case, dilution can prove to be a significant problem in some river systems (see, Miller and Orbock Miller [4] for a more detailed discussion).

Tracers have been utilized to identify both point- and non-point sources of pollution in catchments. Non-point source tracers tend to encompass a wide range of properties (e.g., grain size or mineral magnetics) as the intent is often to identify the source of sediment from a particular soil or land-use type. The tracing of point-source contaminants has historically relied on geochemical parameters, including one or more of the contaminants that is of concern [15,27-35]. The use of the contaminant as a tracer allows for the direct assessment of contaminant provenance along the river, its dispersal rates, and the relative contributions from the identified sources. In many instances, however, trace metals found within the alluvial sediments are derived from multiple natural and anthropogenic sources, and the variations in metal concentrations or chemical forms between the sources may be limited. As a result, trace metal concentration data may prove to be ineffective as a tracer, particularly when the elemental concentrations are near background values. To address the problem, investigators have combined the analysis of multiple elements with multivariate statistical methods (e.g., principle component analysis, cluster analysis, or discriminate analysis) to fingerprint the source [36], an approach widely utilized to determine sediment and nutrient provenance from non-point sources of pollution [37-41]. A significant shortcoming of the approach is that it requires a large number of samples, and its effectiveness depends on the variability inherent within the source materials. Where the variability is larger within a given source than between sources, definitively fingerprinting the contaminant may be difficult even when highly sophisticated statistical analyses are utilized [36].

An alternative approach is to utilize stable isotopic data, rather than bulk elemental concentrations or concentration ratios. The potential use of stable isotopes as a tracer has been recognized since at least the 1960s and is based on the ability to characterize a material on subtle differences in the isotopic abundances of a given element [30]. The isotopes of Pb (specifically 208Pb, 207Pb, 206Pb, and 204Pb) have been particularly well studied and applied as a fingerprinting tool [29,42,43]. In fact, it has been argued that Pb isotopic fingerprinting is one of the most popular methods available for determining Pb sources throughout the world [44]. Nonetheless, in an excellent review, Bird [29] notes that Pb isotopes have only recently been applied to identify contaminant sources within river systems and "[i]n comparison to atmospheric and marine environments, the use of Pb isotopes remains relatively under-studied". In addition to Pb, significant advancements since 1990 in analytical instrumentation, particularly the multi-collector, inductively coupled plasma mass spectrometer (MC-ICP-MS) has dramatically increased the potential for using a wide range of non-traditional stable isotopes (e.g., Cd, Cu, Hg, Li, B, Mo, Sb, Zn) as environmental tracers.

The primary objectives of the remaining discussion are to: (1) provide an overview of the methods through which geochemical tracers have been used to unravel the sources, transport mechanics, and relative source contributions of the trace metal contaminants found within alluvial sediments along a river system; and (2) examine the potential use of stable isotopes, including a number of non-traditional isotopes, as an innovative method of directly documenting trace metal provenance and cycling in riverine sedimentary environments. Particular attention will be given to the use of Pb isotopes in environmental forensic investigations. We will then turn our attention to the innovative use of other, non-traditional stable isotopes that may provide a wealth of information on contaminant sources and transport pathways, but have yet to be extensively applied to riverine environments.

2. Identifying Contaminant Sources Using Spatial Patterns in Geochemical Tracers

Historically, geochemical tracing methods have: (1) utilized a contaminant of interest as a tracer; and (2) documented contaminant provenance by examining the spatial variations that occur within the contaminant/tracer along the river system. In general, the method involves documentation of an element's concentration within the source materials and its comparison to concentrations observed at multiple locations along the channel, including areas presumed to be "clean". In the case of a single point source, a distinctive downstream pattern is often observed in which low concentrations upstream of the source abruptly increase to above background values where the contaminated source materials are introduced to the channel. Concentrations then decline semi-systematically downstream from the source (Figure 1). The sampled media may include water, suspended sediment, or sediments from the channel bed or its floodplain. When the focus is on sources that have provided a long-term, semi-continuous supply of contaminants to the river, channel bed sediments are often the media of choice as they tend to show less temporal variations in concentration through time than is found within the water or suspended sediment [4].

Figure 1. Pb, Zn, and Sb concentrations measured in channel bed sediments of the Rio Chilco-Rio Tupiza drainage system of southern Bolivia. Sb mines are located within the Rio Abaro and Rio Matsu Cuya tributary basins. Polymetallic tin deposits containing limited Sb are located within Rio Tatasi. Point sources of contamination are indicated by abrupt increases in element concentration followed by downstream declines [45].

■ Pb Concentration ■ Zn Concentration ■ Sb Concentration

Rio Matsu Cuya

Rio Tatasi


Rio San Juan delOro

Successive Downstream Sites

An example is provided by the Rio Chilco-Rio Tupiza drainage system of southern Bolivia. In this case, Pb, Zn, and Sb concentrations were used as geochemical tracers of mine contamination. Figure 1 shows that sediment from the Rio Abaroa, contaminated by an Sb mine, produced a geochemical trend extending from sampling site RC-3 to RC-5 that is expected from a single point source of contamination. Concentrations abruptly increase at the tributary mouth and then decline downstream. Villarroel et al. [45] showed that the downstream declines could be attributed to a combination of hydraulic sorting of sulfide minerals from the mines, dilution, and sediment storage associated with channel bed aggradation. The downstream declines in concentration are abruptly altered at RTU-3 where Zn, and to a much lesser degree, Pb concentrations increase downstream of the Rio Tatasi. Since Sb concentrations did not increase downstream of the Rio Tatasi, the majority of the Pb and Zn at RTU-3 is unlikely to have been derived from the Sb mines. Concentrations of all three trace metals decrease downstream of the mouth of the Rio San Juan del Ore as a result of dilution by the influx of relatively "clean" sediment.

The observed downstream declines in contaminant/tracer concentrations from a point source are generally attributed to a combination of five processes, including: (1) the hydraulic sorting of sediment-borne contaminants on the basis of particle size or density which results in the deposition of larger or denser particles (e.g., sulfide minerals or mercury amalgam grains) close to the source; (2) dilution as sediment-borne contaminants are spread out over a larger area, and/or mixed with "clean" sediment from bank erosion or tributaries; (3) sediment deposition and storage within the channel or floodplain, a process that removes the contaminants from active transport and surface deposits; (4) geochemical processes involving sediment-water interactions, such a sulfide oxidation/reduction, precipitation, and or sorption on to reactive particles; and (5) contaminant uptake by biota, a process that is typically negligible along most rivers. Wolfenden and Lewin [46] illustrated that regression analyses could be used to statistically quantify the downstream changes in tracer concentrations produced by the above dispersal mechanisms. More specifically, the regression models provide a means of describing the concentrations of a contaminant/tracer from a single point source at any location along a river (i.e., between sampling sites), while departures from the trend may indicate the influx of a contaminant from an additional source. Since their seminal work other investigators have applied linear, logarithmic, polynomial, power, and other statistical models to rivers located in wide range of hydrologic and climatic environments (Figure 2) [47-49]. The application of these models has shown that interactions between the various dispersal processes can produce different trends for different trace metals along the same river, even when they originate from the same source, presumably because of differences in transport behaviors. In the case of the Rio Chilco (described above), the downstream declines in Sb differ from those of Pb indicating that the trends are metal specific. Thus, when trace metal concentrations correlate with one another along the channel, it is generally assumed that they are derived from a similar source and are transported by similar mechanisms and rates with the channel.

Figure 2. Downstream trends in Pb and Zn concentrations along the River Twymy, Rea Brook and West Allen in the UK. Solid lines represent the best-fit regression models, and show that concentrations tend to semi-systematically decrease downstream of a point source (figure modified from Lewin and Macklin [48] as shown in Miller and Orbock Miller [4]).

Not all rivers that have historically been subjected to point sources of contamination will exhibit the classic downstream declines in the concentrations of the utilized fingerprints. Perhaps the most common reason that the pattern is absent is that the river can be converted from a system dominated by point sources of elemental influx to one characterized by a line source associated with the erosion of materials from the channel banks. For example, Hg contaminated mill tailings were introduced to the Carson River of west-central Nevada primarily between 1859 and 1890 at several distinct locations. During the late 1800s, sediment enriched in Hg-Ag and Hg-Au amalgam grains were deposited along the valley floor downstream of mill processing facilities. Although there is significant variation in the data, Hg concentrations within the alluvial valley fill that were normalized according to the percent silt and clay within the sample generally declined downstream from the mills (Figure 3a). However, Hg concentrations within the modern channel bed deposits quasi-systematically increase from the original sites of mill tailings influx downstream for approximately 25-30 km before subsequently declining (Figure 3b). Miller et al. [50] argued that the observed downstream pattern in Hg concentrations resulted from the progressive erosion of dense Hg-Au and Hg-Ag amalgam grains stored within the valley fill, and their incorporation into the channel bed materials. Subsequent declines downstream were presumably related to dilution caused by the influx of clean sediment from a tributary, and the erosion of uncontaminated, pre-mining deposits from the channel banks which were more extensively exposed downstream. The point, here, is that the absence of the expected downstream trend in concentrations from a suspected point source does not necessarily rule out its potential as a contaminant source.

Figure 3. (A) Downstream variations in Hg concentration normalized for silt and clay percent measured within the alluvial valley fill of the Carson River valley, western Nevada. Increased Hg concentrations occur in response to the influx of Hg from amalgamation processing mills. Variations result from distance from mills of sampling sites, age of stratigraphic units relative to Ag/Au production and Hg influx, hydrologic conditions at time of deposition and other factors; (B) Downstream variations in Hg concentration within channel bed sediments. Concentrations within the modern channel sediments gradually increase as contaminated materials are eroded from the channel banks. The decline in concentration results from the influx of "clean" sediments from tributaries and from the erosion of uncontaminated bank materials (modified from Miller et al. [50]).

S c o O o> X "O NJ

1000 800 600 400 200

Sixmlle Canyon

Gold Canyon ■ ■

"a Lahontan Reservoir \

Carson City m _ B ■ . - s ■ \

20 40 60 80

Distance Downstream from 395 (km)

a) o c o

12 10 8 6 4 2 0

Sixmile Canyon Mineral Canyon

Gold Canyon

Y A Lahontan Reservoir

m ▼ A £ ■A • AY A \

Carson City M .1 ■ T k • " 1 k- \

20 40 60 80

Distance Downstream from 395 (km)


Point Bar T 1997 Flood-


A number of factors can complicate the use of spatial patterns in geochemical tracers to identify a contaminant source. The most significant include the following:

(1) Sedimentological Differences in the Sampled Sediments. Trace metal contaminants are predominantly associated with chemically reactive sediments that are characterized by small particle sizes with high surface areas, a high electrical charge on the particle surface, and high cation exchange capacities [7]. These sediments often consist of clay minerals, Fe and Mn hydroxides and oxides, and organic matter. Sediment composed of quartz, feldspars and carbonates, which tend to be sand sized

(>63 |im) or larger, typically act to dilute the concentrations of contaminants bound to the particles. The net result is that contaminant concentrations can vary as a function of the grain size distribution and composition of the sediments, where samples consisting of larger quantities of sand-sized materials will exhibit lower concentrations. Downstream patterns in concentration, then, may be influenced by changes in the size distribution of the sediments that occurs along the channel or within a particularly sampling location, and mask variations associated with points of contaminant influx.

To compensate for this grain-size effect, it is common to separate the coarse and fine-sediment fractions within the sample, and to analyze only the chemically reactive, fine-grained material, thereby removing the effects of dilution. Alternatively, the bulk sample may be analyzed for the contaminant/tracer and the concentrations normalized according to the amount of fine-grained sediment within the sample. While the latter method provides only an estimate of measured concentrations within the fine sediment, the method is often used because it provides the true concentration in the sediments needed for the calculation of contaminant transport rates, and to assess the potential effects of contaminants on biota. The normalization method, however, generates poor estimates of the concentrations in the fine fraction when there is limited fine-grained material present (<25%-50%) [7]. In addition, Moore et al. [51] showed that within many mine contaminated rivers, trace metals are more uniformly distributed within the sediments. Thus, analyzing the fine-material may lead to erroneous results.

Other forms of normalization, such as normalization to a "conservative" element (e.g., Al, Li, Sc, Ti, or Zr) [7,52-54], have also been used to minimize variations associated with sedimentological differences in the samples when the primary goal is to document spatial variations in the contaminant concentrations for the purpose of source identification.

(2) The Physical Sorting of Contaminants between Channel Bed Features. Closely related to #1 (above), hydraulic sorting processes, including the selective entrainment and deposition of sediment on the basis of particle size or density, often results in a non-uniform distribution of contaminants between (or even within) morphological features found on the channel bed (e.g., pools, riffles, glides, point bars, etc.) [4]. These reach-scale variations in concentration between channel units have the potential to mask downstream trends where differences in concentration along the channel are limited. To compensate for these reach-scale variations it is common to sample specific unit types at each location along the channel, to use composite sampling (where multiple samples are collected from a specific location and combined into a single sample) to reduce field variance, or both.

(3) Localized Variations in Sediment Transport and Deposition. Macklin and his colleagues, while studying the River Tyne in the U.K., found that significant variations in trace metal concentrations were present that could not solely be related to the influx of sediment at tributary confluences. Macklin and Lewin [55] found that the underlying control on the observed pattern was the organization of the valley floor into alternating zones characterized by differences in the transport and deposition of sediment and sediment-borne contaminants along the river. The transport dominated zones tended to be steeper and narrower than average, allowing sediment entering the reach to move more quickly through it with only limited storage. In contrast, reaches dominated by deposition (which they referred to as sedimentation zones) possessed wider valley floors and shallower gradients promoting deposition. The result was that most of the Pb and Zn from historic mining operations were deposited within the sedimentation zones, and metal concentrations were relatively high along these

reaches [56,57]. Since their analysis, transport and sedimentation zones have been recognized along other river systems, such as the Rio Pilcomayo [58] (Figure 4).

(4) Temporal Variations in Sediment Transport. It is not uncommon for sediment, and any contaminants/tracers that they may be carrying, to move downstream through the system as a series of sediment pulses, or slugs, complicating the pattern predicted by standard regression or mixing models [57,59]. Sediment slugs often result from the injection of sediment into the channel as a result of either natural and human activities that exceed the transport capacity of the river [59-61]. The introduced sediments result in a cycle of channel bed aggradation and degradation which propagate downstream from the point at which it was introduced. Slugs, then, are associated with periods of channel instability that range from minor, localized disturbances to basin-scale perturbations in sediment supply that may lead to a total reconfiguration of the valley floor [59,62].

Figure 4. (A) Map showing the distribution of different process zones located along the Rio Pilcomayo of Bolivia [58]. Bedrock canyons represent sedimentation zones, whereas channel and floodplain valley systems represent sedimentation zones as shown in the aerial photograph (B) [58].

Quantifying Relative Source Contributions

Geographical patterns in geochemical tracers provide only a qualitative understanding of the relative contribution of contaminants to the river from the defined sources. More quantitative estimates may be provided by statistical mixing (or unmixing) models. Although exceptions exist (e.g., [48,63,64]), these models have been most extensively used to determine the relative contribution of contaminants from non-point sources of pollution within river, wetland, and lacustrine environments [32-36]. Nonetheless, the approach may provide highly valuable information regarding the primary source(s) responsible for contaminating a specific river reach, particularly given advances in tracer technologies which are describe below.

The mixing model approach assumes that the sediment located within the channel bed or other depositional environment is a mixture of material from the contaminated sources and all other, uncontaminated sediment transported by the river. It is therefore possible to characterize the suspected contaminant sources, the background materials within the basin, and sediments obtained from along the channel for selected geochemical tracers and statistically compare their tracer characteristics to unravel the relative proportion of sediment that was derived from each source [4,40,64]. In most non-point source studies, multiple tracers are utilized to effectively fingerprint the contaminated materials; thus, the approach represents a multivariate statistical method for source identification. Constraints on the mixing model require that: (1) each source type contributes some sediment to the mixture, and thus the proportions (xj, j = 1, 2, ..., n), derived from n individual source areas must be non-negative (0 < Xj < 1); and (2) the contributions from all of the source areas must equal unity, i.e.:

T >=i (i)

In addition, some differences (error) between the values of the m measured parameters, in the source area, a^ (i = 1, ..., m; j = 1, ..., n) and the mixture, b (i = 1, ..., m) are allowed. The residual error corresponding to the ith parameter can be determined as follows:

£i = bij (2)

for i = 1, 2, .m, where a^ (i = 1, 2, 3, ., m; j = 1, 2, ., n) are a measurement on the corresponding ith parameter within the jth source area and xj is the proportion of the jth source component in the sediment mixture. When the number of measured parameters is greater than the number of sediment sources (including contaminant and background materials), the system of equations is over-determined, and a "solution" is obtained using an iterative computational method that minimizes an objective function using a gradient search, thereby obtaining a best fit solution to the entire data set [65]. There are several ways to obtain a best fit, but in previous studies, the objective function, f, has taken the form of the sum of the relative errors [65] where,

f (Xm) = T/ b'\ (3)

or [36],

f (Xl,....,x„) = £/ bi)2 (4)

A criticism of the approach is that there may be a number of solutions that are statistically equivalent, particularly where contributions from a given source approach 0 or 100% (often referred to as the equifinality problem). In other words, similar levels of model performance as measured by a likelihood function can be produced by differing sets of source contributions [66,67]. In addition, uncertainty in the modeling results may be associated with: (1) the inherent variability of the fingerprint within the source materials; (2) the source material sampling density; (3) analytical error; and (4) the changes in sediment characteristics during particle entrainment, transport and deposition which may significantly influence the chemical and physical nature of the sampled deposits, such as grain size and organic matter content [68]. As a result, recent studies have explored methods to reduce and quantify the uncertainty in mixing/unmixing models, including the use of a Bayesian and or Monte Carlo sampling framework [66-70] as well as the inclusion of grain size and organic matter correction factors (e.g., [68]). Validation of these approaches using constructed laboratory mixtures of source materials and synthetic data show the methods hold considerable promise [69].

Attempts to model complex processes generally require simplifying assumptions, and several assumptions are associated with the use of sediment mixing models to assess sediment provenance, including that: (1) the downstream trends in concentration are the product of dilution resulting from the steady-state mixing with tributary and upstream sediment; (2) all of the potential sediment sources are identified and characterized for the geochemical tracers; and (3) the trace metals behave conservatively within the channel (i.e., the geochemical tracers move with the sediment and there is no loss of elemental mass along the channel as a result of geochemical processes). In addition, the approach assumes that the sediment for which the contaminant provenance is being determined is derived directly from the ultimate sediment sources and not its proximal one. It is possible that sediment is eroded from a defined sediment source, transported downvalley and temporally deposited within the channel (or some other site) before being remobilized and transported to its current resting point where it was sampled.

Another assumption for which mixing models have been criticized is that the sediment must be assumed to leave all sources at the same time, and to be transported downstream at an equal rate, so that the sediment arrives at the sampling point simultaneously. Let's take, for example, a sample of sediment from a channel bed accumulated downstream of a contaminant source during a flood event. Assume that the contaminated sediments are derived from a site located a considerable distance from the sampling site, while background sediments are derived from the entire basin, including areas immediately adjacent to the sampling location. Also assume that equal amounts of the sediment from both the contaminant and background sources is deposited at the sampling site as the material is transported through the sampled reach. At the onset of the runoff event sediment from the background sites will reach the sampling location first; thus, the lower portions of the sampling interval will be composed of only background material. As the event continues, material from the contaminant source reaches the site, and equal proportions of the sediment from both sources are deposited. When the mixing model is applied to the sample, it will correctly indicate that a larger relative percent of sediment was derived from the background sources over the depth of material sampled. Thus, sediment

provenance with respect to the deposit has been correctly assessed within the errors inherent in the statistical analysis. However, if the intent is determine the amount of material introduced to the river from the contaminant source relative to background areas, the results will be biased such that the model will underestimate the amount of material derived from the contaminated site. A similar argument applies to differences in transport rates of different sized particles.

A variation in the use of mixing models is to determine the influx of contaminants from individual tributaries rather than distinct contaminant sources. In this case, however, the models do not apply to geochemical tracers, but rather to the changes in elemental concentrations that occur from up- to downstream of a tributary confluence. For example, Church et al. [32] applied a pair of geochemical mixing models to determine the relative contributions of Pb, Zn, and other trace metals to the Boulder River, Montana from several mine contaminated tributaries. The mixing model approach recognizes that the channel bed sediment located downstream of a tributary is a mixture of material from the tributary and the upstream reaches of the axial drainage. Thus, the relative percent of a specific trace metal (contaminant) derived from the tributary is a function of the proportion of sediment transported to the axial drainage by the tributary and the concentration of the contaminant associated with the sediment. Mathematically,

where Pct is the relative percent of the trace metal (contaminant) derived from the tributary, Pst is the relative percent (or proportion) of sediment contributed by the tributary, Cd is the concentrations of an element downstream of the tributary confluence, and Ct is the concentration of an element within the tributary sediment. In early studies, the relative percent sediment yield from the tributary (Pst) was estimated using proxy data, such as basin area, because quantitative sediment transport (yield) data are unavailable for most rivers. However, the use of proxy data is often plagued by large uncertainties. An alternative, utilized by Church et al. [32], is to estimate the percent of sediment derived from a tributary (Pst) on the basis of measured elemental concentrations in the alluvial sediments in that:

where Cu is the concentration of an element upstream from the tributary confluence, and the other parameters are as defined above. Once the relative sediment influx is determined on the basis of the changes in concentration from up- to downstream of the tributary confluence, the estimated value of Pst can be used in Equation (5) to determine the relative contribution of the trace metal by the tributary to the axial drainage.

While the approach provided insights into the relative contributions of sediment from the examined tributaries into the Boulder River, different values of Pst were calculated by Church et al. [32] for the different elements used in Equation (5). In an attempt to reduce the significance of this problem, Mackin development of a mixing model in which Pst was determined using multiple elements and minimization principles from calculus [71]. It is possible that the influx of sediment with very high metal concentrations within a tributary in comparison to the axial channel will result in an over-estimation of sediment delivered by the tributary. To avoid this problem, conservative elements derived from lithogenic sources (rather than anthropogenic sources) (e.g., rare Earth elements) may be used in the

Pct = (Pst x Ct)/ Cd

Pst = \_(Cu - Cd)/ (Cu - Ct)]

analysis. Thus, this multi-elemental method appears to hold considerable promise for reducing the uncertainties associated with the approach.

3. Use of Geochemical Tracers to Construct Pollution Histories

During the past few decades significant attention has been paid to the potential to construct a history of anthropogenic contamination by analyzing geochemical tracers within alluvial and lacustrine deposits associated with a river system [71-80]. The developed histories can form a significant component of any environmental/geological forensic analysis in that they: (1) provide an understanding of the changes in the chemical health of the river through time where water quality monitoring records are limited or absent; (2) allow short-term monitoring records to be placed into an historical context of the changes in sediment geochemistry over periods ranging from a few years to millennia, and, perhaps of most importance from a forensic perspective; (3) can identify a contaminant source and its relative importance as a contributor of contaminants to the river when a source is no longer present. Chronologies of contaminant influx may also provide highly valuable insights into whether water and sediment quality is recovering either in response to remedial actions or on its own, without any human intervention, after the contaminant source has been mitigated. These integrated stratigraphic-geochemical analyses have most frequently focused on vertically accreted overbank floodplain deposits and lacustrine deposits found within reservoirs located along the river.

3.1. Pollution Histories Based on Overbank Deposits

Catchment scale mass balance calculations have demonstrated that significant portions (ranging from 10% to more than 60%) of the sediment eroded from upland areas and transported to the axial drainage system never reaches the basin mouth [4,81]. The bulk of these sediments, and the contaminants that they carry, are deposited and stored with floodplains located along the river. Of particular interest are overbank deposits that are created when runoff events exceed the channel's bankfull capacity and inundate the floodplain surface. Sediments within the flood waters are semi-continuously deposited across its surface, and as sediments accumulate, younger materials bury older producing a trend in which the age of the sediment varies semi-systematically with depth. Assuming that the deposited sediments reflect the geochemistry of the particulates being transported along the channel, geochemical variations as a function of depth in the overbank deposits reflects changes in contaminant influx to the river. It follows, then, that: (1) the overbank deposits contain a long-term (decadal-scale) record of the variations in sediment and contaminant transport rates; and (2) depth/temporal variations in a geochemical fingerprint/tracer within the overbank deposits can be used to construct a chronology of contaminant loadings from anthropogenic sources [72-80,82,83].

The general approach for developing a chronology of pollution loadings from overbank deposits typically involves the collection of a set of sediment cores from the floodplain which are then stratigraphically described, sampled, and analyzed for selected tracers at uniform or semi-uniform intervals below the floodplain surface. In addition to the analysis of the contaminant tracers, the core materials are typically dated using one more methods, including the use of short-lived radionuclides (e.g., 137Cs, 210Pb, 10Be), radiocarbon dating (14C), or optically stimulated luminescence (see Stokes and Walling [84] for a more detailed discussion of dating methods). When the geochemical

data are combined with an understanding of deposit age, the observed variations in concentration can be linked to a history of anthropogenic activities within the basin to identify the most likely contaminant sources.

The use of overbank deposits to create a chronology of trace metal influx is well illustrated by the work of Macklin and Klimek [72] in upper Vistula basin in southwestern Poland. In this case, they analyzed overbank sediments associated with levees and inset fills confined by a dike constructed near the turn of the 19th century, and within sediments that occurred outside the dike-protected area which were deposited prior 1800. The identified deposits were dated and analyzed for Pb, Zn, and Cd, the latter of which were considered to be effective tracers of mining debris contained within the catchment. The age and geochemical data were then related to one another to create a general history of contamination that could be linked to anthropogenic activities. The results showed that sediment quality progressively declined through time within the upper Vistula basin, particularly after about 1850. Initially, elevated trace metal concentrations occurred at sites located immediately downstream of mining districts within the Przemsza River, a tributary to the Vistula, but became much more widespread, affecting the entire basin during the past 100 years.

Hudson-Edwards et al. [85] note that the investigation of pollution histories using overbank floodplain deposits fall into two general categories, those conducted at a relative coarse-scale in which general trends in pre-, peak, and post-metal contamination are determined using geochemical tracers (e.g., [32,81-88]), and analyses performed at a fine-scale during which the sampling interval is measured in centimeters, and the sediments are subdivided into discrete periods on the basis of the degree to which they are contaminated (e.g., [88-91]). In both cases, there are two important assumptions inherent in the development of pollution histories: (1) the observed changes with depth in the geochemical tracers are a reflection of the metal content within the suspended sediment at the time of deposition; and (2) once deposited, the position of the tracers within the sediment is not significantly affected by post-depositional processes. Data regarding the degree to which the overbank deposits, or those found in oxbow lakes and reservoirs for that matter, reflect the transported sediment are surprisingly limited. The available data suggest, however, that the geochemical signal may vary spatially across the floodplain as a function of distance from the channel, the microtopography of the floodplain surface, and the magnitude of flooding, among of number of other variables [4,80,92]. The geochemical signal also may be affected by changes in channel form (e.g., channel straightening), or other sedimentological process that influence the transfer of sediment from the channel to the floodplain [93,94].

Particular attention has recently focused on the role of large floods in controlling the physical and chemical composition of the overbank deposits [75,82,93-95]. Larger floods appear to deposit relatively coarse-grained sediment on the floodplain characterized by lower trace metal concentrations than finer-grained sediments deposited during more frequent, lower magnitude events [75,96,97]. To compensate for these grain size differences, it is common to analyze only the fine-grained sediment fraction, or alternatively to normalize the concentration on the basis of the percent of fine sediment in the sample or to a conservative element. These methods are not always effective. Babek et al. [75] found that short-term, flood-related variations in contaminant concentrations remained, even after normalization to a "conservative" element (Al). However, magnetic susceptibility values were found to correlate with particle grain size. More specifically, fine-grained layers exhibited higher magnetic

susceptibility values in large part because of a higher concentration of paramagnetic clay minerals and phyllosilicates. Thus, Babek et al. [75] were able to use magnetic susceptibility data as a proxy for grain size in order to minimize the effects of flood-related grain size variations on trace metal concentrations [7,63].

With regards to the second assumption, the degree of post-depositional migration of a geochemical tracer can be significant, depending on the reactivity of the tracer and the physio-chemical environment [77,85,92,94]. Hudson-Edwards et al. [85], for example, demonstrated using detailed speciation and mineralogical data from the River Tyne in the UK that trace metals often migrate over depths of a few centimeters within the floodplain. The magnitude of migration varied between elements such that the more tightly bound elements (e.g., Pb) exhibited the least mobility. It has also become clear that the degree of mobility can be enhanced by environmental change, such as channel incision which leads to a drop in water table elevations, redox conditions, and the weathering of the existing minerals (e.g., sulfides). Mobilization also tends to be more significant in areas of high water table fluctuation [76]. It follows, then, that the results from fine-resolution data are more likely to be affected by post-depositional processes than those associated with coarser-resolution data sets.

Given the difficulties of creating chronologies of contaminant influx from floodplain deposits, the results are likely to contain uncertainties which are difficult to quantify [75]. Nonetheless, when combined with other data, the results can serve as an important component of any environmental forensic analysis.

3.2. Pollution Histories Based on Lacustrine Deposits

Overbank floodplain deposits typically exhibit hiatuses and irregularities in deposition that complicate the construction of pollution histories [76,97]. In contrast, reservoirs represent depositional environments that are characterized by high rates of nearly continuous sediment accumulation [98,99]. These high rates of accumulation tend to reduce the potential for post-depositional contaminant migration by early diagenetic processes [100,101], while the continuous deposition of sediment (and the geochemical tracers associated with them) allows for the development of pollution histories of high temporal resolution. As a result, reservoir deposits are frequently targeted to document temporal variations in contaminant influx to a river upstream of the dam (e.g., [71,98,101-103]). Other "lacustrine" deposits such as oxbow lake and estuarine/coastal sediments also have been used [104,105], although the post-depositional migration of trace metals is more problematic in these environments.

The approach to construct a history of contaminant influx using lacustrine sediments is similar to that applied to overbank floodplain deposits. Audry et al. [101], for example, constructed a record of trace metal influx to the Lot River, France, a river that has been subjected to contamination by mining and smelting activities since the late 19th century. Sediment age was determined by 137Cs analysis and by extrapolating between depths of known age on the basis of average sedimentation rates. Cd and Cu were selected as geochemical tracers of the smelting wastes, and their concentrations were normalized to Sc, a conservative element, to help differentiate between lithogenic and anthropogenic sources of the two metals. The constructed record showed that prior to ~1975, Cd/Sc and Cu/Sc ratios were highly variable, reflecting high trace metal inputs from smelting activities (Figure 5). From 1975 to 1985, Cd levels increased as smelting activity grew and smelting waste piles within the catchment

expanded. Interestingly, Cu concentrations apparently declined as a result of improvements in processing technologies for metal recovery which tended to be more efficient for Cu than for Cd. An accidental seepage release of trace metal contaminants in 1986 created the peak values in concentration at that time, which were followed by a significant decrease in both Cd and Cu concentrations at the top of the core. This latter decrease was attributed by Audry et al. [101] to remediation efforts undertaken after 1987. The concentrations at the top of the core remained well above background values indicating that the remediation efforts only were partially successful.

Figure 5. Vertical (depth) variations in Cd/Sc and Cu/Sc ratios within samples from a dated sediment core extracted from the Cajarc reservoir, located along the Lot River, France. Variations reflect changes in Cu and Cd influx to the river in response to alterations in processing methods and, after 1987, remediation efforts (modified from Audry et al. [101]).

As shown for the overbank floodplain deposits, multiple contaminants/geochemical tracers can be analyzed within the lacustrine or coastal cores. It is important to note, however, that not all variations in contaminant/tracer values are necessarily related to changes in the direct influx of anthropogenic contaminants. Miller and Mackin [71], for example, used reservoir sediments from Lake Emory to document temporal changes in trace metal influx to the Little Tennessee River of western North Carolina since dam closure in 1925. They found Cu, Pb, As, and other trace metal and metalloid concentrations with the reservoir deposits frequently exceeded the consensus-based threshold effect concentrations of MacDonald et al. [106] for aquatic biota. This consensus-based threshold represents the geometric mean of threshold effect guidelines put forth in multiple studies of aquatic biota and is

likely to be applicable to a wide range of environments. The data, then, suggest that the metals have the potential to negatively impact aquatic health. In addition, Ni, Cr, and Zn concentrations were found to be relatively high near the bottom of the reservoir deposits, and decrease semi-systematically toward the surface of Lake Emory (Figure 6). The elevated concentrations were coincident with relatively high rates of sediment transport within the river, as demonstrated in Lake Emory by sedimentation rates prior to 1963 that were approximately five times higher than those measured between 1963 and 2010. Using a variety of data, they argued that the trace metals were largely derived from natural sources, specifically the erosion of subsurface soil horizons enriched in trace metals weathering from sulfide mineralized units within the underlying bedrock. The increased loading rates to the river apparently resulted from a period of extensive land-clearing during the late 1800s and early- to mid-1900s that also resulted in the deposition of post-settlement alluvial deposits along the floodplain [107]. The observed decreases in Ni, Cr, and Zn concentrations presumably resulted from improved watershed management practices which limited upland erosion in the basin during the latter part of the 20th century. In contrast to Cr, Ni, and Zn, Cu concentrations within the reservoir semi-systematically increased from the bottom to the top of the reservoir deposits (Figure 6), suggesting that Cu loadings to the channel have increased during the past several decades. The additional Cu was most likely derived from Cu-based pesticides that were extensively used on tomato fields within the watershed.

3.3. Extrapolating Site Data Using Chemostratigraphic Methods

The history of contaminant provenance and transport rates can abruptly change along a river as a function of the mechanisms and rates of sediment-born contaminant dispersal, the space and environments available for sediment deposition, the addition of new contaminant sources, and changes in the contaminant-sediment relations, among a host of other factors. As a result, a question that often arises is how to most effectively correlate the timing of specific contaminant loading events documented for one site with that of another up- or downstream site. The inability to correlate contaminant histories along a channel may introduce significant uncertainties into the extent and duration to which a specific contaminant source has impacted both sediment and water quality in the riverine environment.

Recent studies have illustrated that geochemical tracers, when combined with the concepts inherent in chemostratigraphy, can provide a powerful tool for correlating stratigraphic units, and the pollution histories based on these units, between sites. Chemostratigraphy can be defined as the detailed characterization of the spatial changes in the chemical nature of sedimentary strata, and the identification of sedimentary packages (units) with similar geochemical traits for the purpose of correlating geographically distributed units [108]. Few studies have attempted to apply these methods to recent riverine deposits; rather most chemostratigraphic studies have focused on marine, lacustrine, and lithified strata (e.g., [108-113]). However, the studies that have been conducted to date show that geochemical signatures associated within anthropogenic sources of contamination may be used to effectively correlate otherwise sedimentologically homogenous deposits over large downstream distances [56,114,115]. In addition, if the timing at which a specific geochemical fingerprint was introduced to the river is known, then the fingerprinted materials may be used to constrain the relative and absolute age of the deposits along the drainage system.

Figure 6. Vertical (depth) variations in Cu (A), Pb (B), Ni (C), and As (D) concentrations within samples collected from a sediment core extracted from Lake Emory along the Upper Little Tennessee River, western North Carolina. Pb, Ni, and As concentrations decrease from 1925 to present. High values correlate with the erosion of metal enriched subsoils during land-use change. In contrast, Cu loadings to the river have increased since dam closure in 1925 in response to the use of Cu based pesticides within the watershed (from Miller and Mackin [71]).

Ni (mgikg)

As (mgttg)

The chemostratigraphic approach is based on the identification of recognizable geochemical signatures which can be traced between sites and which allow for unit correlation. The application of chemostratigraphic methods to riverine environments is well illustrated by the work of Chillrud et al. [114] for the upper Hudson River in the northeastern U.S. The upper Hudson River was subjected to significant inputs of Cd, Sb, Pb, and Cr during the latter part of the 1900s. Chillrud et al. [114] used geochemical tracing methods on sediment cores collected upstream of three dams located a total distance of 24 miles apart to construct a history of contaminant loadings to the river. An important component of the analysis was the chemostratigraphic correlation between the sites

206 207

on the basis of Pb/ Pb ratios (Figure 7). By doing so, they were able to use the age constraints provided by 137Cs data within the upstream core to precisely date the deposits located downstream. This allowed for a detailed analysis of the influx of sediment between the sites from the early 1960s to the late 1970s on the basis of the observed differences in tracer concentrations. They found that 3 to 4 times more sediment entered the river between the upstream- and downstream-most coring sites than was transported from upstream portions of the watershed.

206 207

Figure 7. Vertical (temporal) trends in Pb/ Pb ratios within three cores obtained from reservoirs at three different locations along the upper Hudson River. Red squares represents samples from mile post (mp) 188.5 (upstream), orange circles from mp 163.6 (downstream), and open blue triangles from samples collected at mp 169 in 1995 (A).

137 206 207

Sediment ages are assigned ages on basis of Cs data. Vertical trends in


similar between the three sites, but are offset vertically. (B) Temporal trends based on both

137 206 207

Cs data and chemostratigraphic correlation on the basis of Pb/ Pb ratios between the sites so that the in 206Pb/207Pb trends match (from Chillrud el al. [114]).


o. irai O)

Q "S ê

ra "ft a) £

"SI í= ra o

1 O x '

2000 1990 1980 1970

1960 1950 1940 1930

I 1 ! 1 i 1 & □

— Va

VV"- UU .......

f I I j **

1 aA 1 1 1 1 1

CM Csl


The point to be made is that while chemostratigraphic methods are cost and effort intensive, they provide a highly underutilized approach for linking together histories of contaminant source contributions constructed for various sites along the river. In doing so, a wealth of information can be extracted from the deposits with regards to the extent and magnitude of contamination associated with each of the identified contaminant sources.

4. Elemental Isotopes as a Geochemical Tracer

As noted in the introduction, a wide range of physical and chemical parameters have been effectively utilized as contaminant tracers in river systems. The bulk concentration of trace metal contaminants, or their related ratios, has been most extensively used as a tracer in riverine environments, particularly when the contaminants are derived from mining, milling, and refining processes. In many instances, however, trace metals found within the alluvial sediments are derived from multiple natural and anthropogenic sources, and fingerprinting contaminated sediment on the basis of bulk elemental concentrations may be ineffective, particularly when the elemental concentrations are near background values. One approach to address the issue has been to use isotopes as a geochemical fingerprint/tracer. Recent studies have demonstrated that selected elements (e.g., Pb) exhibit stable isotopic abundances that fall within a narrow range of values in geological materials. Moreover, many anthropogenically derived substances exhibit isotopic values that differ from non-mineralize (lithogenic) rocks. Thus, these isotopes can serve as a highly sensitivity tracer that allow for the determination of sediment and contaminant provenance even where multiple sources exist, or where elemental concentrations are near background values [35,116]. In the following sections, the use of isotopes as a geochemical tracer of trace metal contaminants is explored.

4.1. Pb Isotopes

To date studies have focused extensively on the use of Pb isotopes as an environmental tracer in a wide range of media, including air and aerosols [117-122], snow and ice [123-127], soils [128-132], lacustrine and reservoir deposits [133-136], wetlands and peat [137-140], plants, mosses, and tree rings [141-143], human tissues and blood [121,144,145] and other biota [146-152]. Within riverine environments there is little question that the isotopes of Pb have been most extensively utilized as an isotopic tracer in forensic investigations (see Bird [29] for a detailed review of Pb as a contaminant tracer in rivers).

Pb has four stable isotopes that are used in environmental forensics, including the daughter products 206Pb, 207Pb, and 208Pb which are derived from the radioactive decay of 238U, 235U, and 232Th, respectively (Table 1). The fourth isotope, 204Pb, has no known radiogenic parent. The natural Pb isotopic composition of a geological material depends upon the relative proportion of U, Th, and Pb in the system, various mixing processes associated with metamorphism, and the age of the rocks and minerals as some U and Th will decay through time to Pb [42,145]. The isotopic abundances of Pb found within rocks, sediments, soils, etc. are typically reported as ratios (e.g., 206Pb /204Pb), but no

consistent set of ratios have been used in the literature. While one study may report the values in terms

206 207

of 206Pb/207Pb, others, particularly geologically oriented investigations, may report the same measured

207 206

abundances in terms of Pb/ Pb, making it difficult to compare the results between investigations. Moreover, many geological studies have focused on the ratio of the radiogenic isotopes to 204Pb (i.e., 206Pb /204Pb, 207Pb/204Pb, and 208Pb/204Pb) [29]. In contrast, most studies using Pb isotopes as a geochemical tracer have avoided the use of 204Pb because of its low relative abundance (Table 1), and

the difficulty of measuring it at low concentrations. In addition, Sangster et al. [152] found that 86% of

206 207 208

the discriminative (fingerprinting) power of Pb isotopes is due to

zuoPb, zu^b, and ^Pb.

Table 1. Selected properties of Pb isotopes used in environmental studies.



Half-life and Parent Isotope

1.48% 23.60% 22.60% 52.30%

204Pb (stable) 18U to 206Pb; / life of 4.5 x 109 years 15U to 207Pb; / life of 7.1 x 108 years Th to 208Pb; / life of 1.4 x 1010 years

The pronounced effectiveness of Pb isotopes as a geochemical fingerprint/tracer is related to several factors. First, the radiogenic isotopes of Pb can be measured with a high degree of precision and accuracy. Early on the analysis was primarily carried out by Thermal Ionization Mass Spectrometry (TIMS), a method that provided precise measurements (>0.005% RSD) [153], but required considerable sample preparation [29,154]. More recent investigations have relied extensively on the use of Inductively Coupled Plasma-Mass Spectrometry (ICP-MS) which requires less time and effort. The development and availability of multi-collector (MC) ICP-MS and magnetic sector (MS) ICP-MS, both of which produce highly precise (ca. > 0.002% RSD) [29,154] and accurate isotope measurements, have particularly increased the use of ICP platforms for sample analysis. The highly-precise measurements of the Pb isotopic ratios make them a particularly sensitive tracer that allows specific Pb sources to be identified with a minimal number of samples and statistical manipulation [42].

Second, effective fingerprints/tracers exhibit conservative behavior over a wide-range of environmental conditions. In the case of sediment in river systems, this means that the tracer moves with the sediment without significant loss in elemental mass, or a change in the isotopic abundances [29,155]. Pb has a high affinity for particulate matter, and once deposited, is generally immobile [132,135]. Moreover, the heavy nature of Pb, combined with the low relative atomic weight differences among its isotopes, limit the degree to which mass-dependent isotope fractionation by physical, chemical, and biological processes occur [29,42,145,153,156]. This is an extremely beneficial trait in that the Pb isotopic abundances found within the ore deposits do not change during mining, smelting, or other industrial or biological processes. Thus, the Pb released into the environment from a contaminant source retains the isotopic signature of the source material as it is dispersed through the river system. Observed differences in Pb isotopic ratios between the source materials and the alluvial sediments at any given site can therefore be attributed to its mixing with sediments from other Pb sources [42].

Third, the most effective tracers exhibit a wide range of values in geological materials. Although the isotopic ratios of Pb in ore deposits on a global scale are highly variable [29], there is a strong tendency for ore-derived Pb to be more radiogenic. Therefore, Pb ratios within anthropogenically derived materials created from those ores tend to differ from the more geogenic Pb found within non-mineralized country rocks [157]. It is often possible, then, to effectively fingerprint both natural and anthropogenic sources of Pb contained within alluvial sediments.

4.2. Applying Pb Isotopes as a Tracer in Riverine Environments

As shown earlier for other geochemical tracers, spatial variations in Pb isotopic ratios along a channel can be used to assess the sources of Pb at any given point, and the downstream changes in Pb provenances that occur. The analysis is based on the assumption that the analyzed sediment represents

a mixture of particles from all of the potential source areas in the catchment, and that the relative abundance of the Pb isotopes reflects the contributions from each source. Theoretically, however, there is some evidence to suggest that physical sorting processes may lead to differences in the observed isotopic ratios from one location to another; for example, Bird [29] notes that Pb isotopic signatures may differ between sediment sizes (e.g., fine-grained, chemically reactive silts and clays and coarser mineral fragments). This is supported by recent sequential extraction analyses. Sequential extraction is a method that treats a sample in different reagents in a stepwise manner to determine how, and with what, a trace metal is bound within the sediment. It has demonstrated that Pb isotopic ratios can vary between Pb that is associated with different physiochemical phases (e.g., the exchangeable, reducible, oxidizable, and residual phase) that are often grain size dependent [158-161]. Thus, it is possible that downstream changes in sediment size or the partitioning of sediments on the basis of size or density into specific depositional units may contribute to differences in isotopic values that occur. While this is a topic that is in need of additional research, differences related to these processes are likely to be minor along highly contaminated river systems. It follows, then, that changes in Pb isotopic signatures within alluvial sediments primarily reflect the input of Pb from different sources characterized by different Pb isotopic abundances.

An excellent example of the use of Pb isotopes to assess Pb provenance is provided by Shepherd et al. [162] for the River Wear located in northeast England. They analyzed a variety of geological materials within the catchment as well as the channel bed sediments along an 87 km reach of the river. The isotopic variations measured within the channel bed sediments are shown in Figure 8. As is commonly done, Pb concentration data are also shown on the figure to provide additional details on the levels of contamination that are associated with the various Pb sources. The data reveal that at Site 1 the 208Pb/206Pb ratio is comparable to geological materials found in the outer mineralized zone of an orefield. Thus, the high Pb concentrations found in the bed sediments are likely to reflect the influx of Pb-contaminated sediment to the river from abandoned mines, waste piles and an old processing plant associated with historic mining operations that produced about 60,000 t of lead concentrate [163]. The rapid downstream declines in concentration from Sites 1 to 5 are thought to result from dilution and the deposition of dense, finely milled Pb ore within the channel bed. Downstream of Site 5, Pb concentrations increase significantly until reaching Site 9. Interestingly, however, 208Pb/206Pb ratios continue to decline. The decline in


ratios reflects the influx of mine waste from two tributary valleys where the Pb ore was derived from the inner zone of the orefield. The rocks within the inner zone exhibit a much lower 208Pb/206Pb ratio than those found within the outer zone of the orefield. The 208Pb/206Pb ratio between Sites 10 and 11 increases abruptly, reflecting the input of Pb from Bollihope Beck, a tributary containing another cluster of mines. Pb concentrations, however, do not increase indicating that the tributary provides a relatively limited quantity of the Pb to the main stream of the River Wear (about 18% of the total). Downstream of Site 11, there are no further direct inputs of Pb from mining operations in the catchment. Without the isotopic data, decreasing Pb concentrations downstream of Site 11 could reasonably be interpreted as reflecting the dilution of sediment-born Pb by relatively "clean" sediments from the underlying bedrock. However, the increasing


ratio suggests that other source(s) of Pb characterized by a higher 208Pb/206Pb ratio are entering the channel. Shepherd et al. [162] argued that this additional Pb was derived from the influx of diffuse anthropogenic sources from tributaries draining an abandoned coalfield. The abrupt

rise and following increase in 208Pb/206Pb between sites 24 and 27 were interpreted to result from either: (1) an increased influx of diffuse Pb in the vicinity of Durham; or (2) the influx of tetra-ethyl lead from leaded gasoline. The point to be made by the analysis conducted by Shepherd et al. [162] is that a much more detailed understanding of Pb sources and their relative contributions to the contamination of the river was made possible by combining Pb concentration data with Pb isotopic fingerprinting and tracing methods.

Figure 8. Downstream changes in total Pb concentration (red circles) and 208Pb/206Pb ratios (orange circles) along the River Wear, Northern England. Sample sites 1 through 27 are shown on ordinal axial and cover a distance of 86.5 km. Error bar refers to isotopic measurements (from Shepherd et al. [162]).

In addition to examining the spatial variations in isotopic ratios along a channel, Pb sources are often identified by analyzing three-component scatter diagrams depicting differences between two isotopic ratios [64,164-170]. When the values of the analyzed samples form a strong linear trend, Pb within the samples is often interpreted to be derived from two primary sources, where the sources exhibit minimum and/or maximum values of the examined ratios at the end of the line. The trend, then, represents a mixing line between the sources. In the case of alluvial sediments, one end-member (Pb source) is typically anthropogenic, whereas the other represents Pb from the underlying bedrock or the soils developed from it [29,171]. The signature of the latter "geogenic" Pb may be determined by analyzing a number of different materials including the underlying bedrock [172], channel bed materials from uncontaminated tributaries or areas upstream of the contaminant influx [145], uncontaminated alluvial terraces [64,172,173], or old alluvial sediment found at depth within a sediment core [32].

Miller et al. [64,172], for example, used bivariate scatter diagrams for 206Pb/207Pb and 206Pb/208Pb to assess the sources and contaminant transport processes along the Rio Pilcomayo downstream of the Cerro Rico de Potosi precious metal-polymetallic tin deposits of Bolivia. Mining of the deposits has continuously occurred since 1545 and the river is severely contaminated by a wide-range of trace metals [172,174]. Miller et al. [172] found that (1) the 206Pb/207Pb and 206Pb/208Pb ratios differed

between bedrock units within the catchment, and (2) alluvial sediments collected from pre-mining terrace deposits formed a mixing line defined by the isotopic ratios of Ordovician and Mesozoic rocks. In contrast, samples collected in 2000 from the highly contaminated channel bed exhibited a significantly different trend (Figure 9) in which the isotopic end members were formed by the Ordovician rocks and mine/mill processing waste from the operations at Cerro Rico. In 2000, then, the channel bed material was dominated by natural Pb from the underlying Ordovician rocks and Pb from the ore deposits that were released anthropogenically into the river system. Pb input from the Mesozoic rocks was apparently minor and the isotopic signature was masked by the large inputs from the other two sources.

Figure 9. (A) Pb isotopic ratios measured in pre-mining (and uncontaminated) terrace deposits (purple line) vs. those measured within Ordovician, Mesozoic and Cretaceous rocks. Pre-mining terrace deposits appear to be composed of sediment from Mesozoic rocks with relatively high Pb ratios and Ordovician rocks (from Miller et al. [172]). (B) Pb Isotopic data collected in 2000 from modern channel bed (red circles and brown squares). High ratio end member is consistent with isotopic values obtained from ore deposits mined and milled from Cerro Rico, indicating Pb is derived from the mines and the underlying rocks (particularly Ordovician Rocks (From Miller et al. [64]).

(A) 0510

lXJ ,,-nß-fV illOÎ

1.13 1.15 1.17 1.19 1.21 1.23 1.25

206 Pb/ 207Pb

• Mesozoic Rocks a Ordovician Rocks + Pre-Miining Terrace Deposits

m Cretaceous Ouartzites 0 Other Rocks

^gj 0.4900

0.4850 -

n 0.4800 -

îL 0.4750

° 0.4700 -

0.4700 -


Ordovician Rocks


1.1500 1.1550 1.1600 1.1650 1.1700 1.1750 1.1800 1.1850


♦ LowChannel ■HighChannei oShaie »Terrace Deposits

One of the significant advantages of using Pb isotopic fingerprinting methods is that accurate estimates of the relative contributions of Pb from the two sources can be made where a mixing line is defined by two distinct end-members. The estimates are typically based on a simple binary model [29,168,175] where,

where %A is the contribution of source A in a sample and 20xPb/20xPb and 20xPb/20xPb are the average ratios in the sample (S), end-member for source 1 (A) and end-member for source 2 (B). In the case of the Rio Pilcomayo, Miller et al. [172] were able to use the approach to estimate that the contribution of Pb from the upstream mining operations in the samples collected in 2000 ranged from >99% to about 56% of the total Pb present (note that these values were recalculated from Miller et al. [65] on the basis of a more accurate understanding of the isotopic signatures of the ore deposits).

The data collected along the Rio Pilcomayo indicate that a significant proportion of the Pb was derived from the underlying bedrock, a conclusion that has been reached by other investigators as well [172,176]. In addition, Miller et al. [172] found that the Pb isotopic ratios measured in the channel bed sediment changed from 2000 to 2002 to 2004. The differences appear to be related, in part, to the fact the mined ores were characterized by two distinct isotopic signatures, which they referred to informally as the Porco and Cerro Rico type ores. Ore bodies characterized by more than one isotopic signature have been described for a number of locations, and are not an uncommon phenomenon for many metallic deposits. In 2000, the channel bed samples immediately downstream

206 207 206 208

of the mills exhibit 206Pb/207Pb and 206Pb/208Pb ratios similar to Porco type ores, whereas in 2002 samples collected from along this reach exhibited a signature comparable to Cerro Rico type ores (Figure 10). This led them to argue that the change in isotopic ratios from 2000 to 2002 was likely to represent the processing of differing ore types.

206 207

Another interesting observation from along the Rio Pilcomayo is that 206Pb/20'Pb values within samples collected from the lower channel tended to decrease downstream toward the values observed

206 207

within the Ordovician rocks. In contrast,


ratios measured in samples collected in 2002

206 207

showed no downstream trend; rather, the 206Pb/207Pb ratios were either indistinguishable from the ratios within the Cerro Rico type ore deposits, or fell between the two types of ore deposits (Figure 10). The change in spatial patterns of the isotopic ratios was coincident with a decrease in: (1) Pb concentrations at any given location along the channel; and (2) the downstream extent of highly elevated Pb concentrations (Figure 10a). Miller et al. [172] suggest that the observed changes in Pb concentration and Pb isotopic values reflected annual differences in system hydrology and upland

sediment inputs to the river. During wet years, such as 2000, the highly dynamic nature of the Rio Pilcomayo allowed larger quantities of Pb-contaminated sediment from the mines and mills to be transported downstream over longer distances, resulting in relatively high Pb concentrations. However, the proportion of the Pb derived from the mines decreased downstream, as indicated by the isotopic ratios, because significant quantities of uncontaminated fine-grained sediments were eroded from upland areas and delivery to the channel via tributaries as a result of significant runoff. During dry years, such as 2002, the transport of Pb-contaminated sediment from the mines and mills is reduced

(20xPb /20xPb) — (20xPb /20x % A = ±--;

decreasing the overall concentrations within the alluvial bed sediments. However, because there is much less sediment eroded from the surrounding uplands and delivered to the channel during the dry years, most of the Pb found in the channel is derived from the upstream mining operations. The study by Miller et al. [172] demonstrates that the source of Pb to the channel can vary over relatively short timeframes (i.e., a year) as a function of the hydrology and sediment transport processes that characterize dynamic catchments. To date, very few studies have examined these types of short-term variations in isotopic signatures within alluvial sediments, and the topic clearly requires additional attention.

Figure 10. (A) Changes in Pb concentrations within channel bed sediments along the Rio

206 207

Pilcomayo from 2000 to 2004. (B) Comparison of Pb/ Pb ratios measured in channel bed

206 207

sediments sampled in 2000 and 2002. In 2000, 206Pb/20'Pb ratios decreased downstream indicating that the relative contribution of Pb from mining operations at Cerro Rico declined as "clean" sediment from tributaries was introduced to the channel during relative high runoff events. In 2002, runoff events were limited during relatively low rainfall. As a result, the Pb ratios were generally indistinguishable from the Cerro Rico type ore deposits, which suggests that nearly all of the Pb was from mining operations (from Miller et al. [172]).

-E? "a>

<L> CJ

3500 — 3000 -2500 u 2000 -1500 -

1000 -

500 -0 —

S 1.18 -

1.17 -

2000 Channel 2002 Channel ~ 2004 Channel

Spatial decline in concentration from 2000 to 2004


50 100 150 200

Distance Downstream (km)

i y = -0.001Ln(x) + 1.2008 R2 = 0.16 • Range of values for Cerro Rico tailings m •

• •—

• ■ 9 • 9 •

f ■ ■ Range of values for Porco type tailings _

■ * -0.0022Ln(x) = R2 = 0.36

1 1803 ■ a ■ ■ -

Distance Downstream (km) • 2002 Channel m 2000 Channel

Within larger river systems contaminated by multiple Pb sources, or which are underlain by a wide range of bedrock types, the Pb isotopic signature of the alluvial sediments will reflect a complex mixture of material from the Pb sources [177]. As a result, well-defined mixing lines may not exist on binary scatter plots, and determining the relative contribution of Pb from the identified sources becomes much more difficult. However, where the Pb sources can be defined on the basis of isotopic ratios, statistically based mixing models may be used to determine the relative contributions of Pb to the river [168,169,172]. Bird et al. [169], for example, were able to use a two-tiered approach to assess Pb provenance within the River Maritsa catchment of Bulgaria and Turkey. Within smaller tributary basins, a binary approach was used to assess the contributions of Pb from mining operations to the alluvial sediments. However, the Pb isotopic composition of the lower catchment areas where sediment was derived from a wide range of geological materials inhibited the use of this binary method. Thus, they applied a multi-element geochemical mixing model (similar to that described earlier) to assess sediment/Pb provenance. A similar approach was used by Bird and his colleagues [168] along the lower Danube Catchment in Eastern Europe.

Miller et al. [172] used a slightly different approach by applying a mixing model directly to Pb isotopic ratios. In this case, they found that the majority of samples collected from agricultural fields

206 208 206 207

located along the Rio Pilcomayo (described above) exhibited 206Pb/208Pb and 206Pb/207Pb ratios within

the limits formed by the four identified Pb source, (including, the Porco and Cerro Rico type ores, and Mesozoic Rocks and the Ordovician Rocks) (Figure 11). They were then able to use a mixing model similar to the one described earlier to determine the relative contribution of Pb within agricultural soils located along the Rio Pilcomayo. Miller et al. [172] found that within agricultural terrace soils, both Pb concentrations and the percentage of Pb from the mining and milling operations at Cerro Rico semi-systematically decreased downstream. In upstream reaches (within 30 km of the mills) Pb from mining accounted for more than 80% of the Pb found in the examined agricultural fields. Downstream at Sotomayor, located approximately 170 km from the mills, the relative contribution of Pb from the mines was highly variable between fields, but ranged from approximately 15% to 35%.

Pb isotopic data also have been combined with bulk concentration data to differentiate changes in contaminant provenance within floodplain deposits through time [45,170]. Villarroel et al. [45], for

206 207

example, used a combination of bulk Pb concentrations and 206Pb/207Pb ratios to assess the source of Pb and other correlated trace metals contaminants to floodplains used for agriculture along the Rio Chilco-Rio Tatasi River system of southern Bolivia. Figure 12 shows, as described earlier, that Pb and Zn were delivered to the axial drainage system from Sb mines located along the Rio Chilco and Pb/Zn mines located within the Rio Tatasi basin. Since Pb was derived from both of the major mining districts in the catchment, variations in Pb concentrations alone could not be used to assess which set of mines provided the majority of the trace metal contaminants to the river, or when. However, the

206 207

incorporation of 206Pb/207Pb data into the analysis allowed the provenance of the Pb within the floodplain deposits to be determined for different depths/ages of deposition (Figure 12).

Figure 11. Soil samples from floodplain agricultural fields collected along the Rio Pilcomayo, Bolivia, exhibit Pb isotopic values that fall within the range established by four Pb sources, including Ordovician rocks, Mesozoic Rocks, and two types of ore deposits from Cerro Rico. As a result, mixing models could be used to assess the relative contribution of Pb from the underlying bedrock and mining operations (from Miller et al. [172]).

Figure 12. Pb from the local geological units, as well as from Abaroa, Chilcobija, and Rio

206 207

Tatasi mines can be fingerprinted on the basis of total Pb concentrations and Pb/ Pb ratios (A). The defined geochemical fingerprints can be used to assess changes in Pb contributions to the river system through time as recorded in floodplain deposits sampled at two different sites (B, C) (from Villarroel et al. [45]).

A different approach has occasionally been used to assess the relative contribution of Pb where three end-member sources can be identified on the basis of Pb concentration and a Pb isotopic ratio (e.g., [104,176]). In this case, sample concentrations and their associated isotopic ratios are graphically plotted. Linear mixing lines are subsequently calculated and drawn between the three sources creating a ternary diagram. The three mixing lines between the sources are calculated using two equations. The equations are:

_ (R,C, Xi ) + (R2C2 X2)

Cmix (8)

Cmix _ (ClXl) + (C2X2) (9)

where R1, R2, and Rmix are the isotopic ratios for end-member 1 and end-member 2 as well as a mixture composed exclusively of these two Pb sources, respectively, C is the concentration of Pb in the mixture and end-members, and X is the relative (fractional) contribution of Pb from the end-members. The fractional contribution from the two end-members must sum to unity (1). The relative contributions of Pb from the three sources then can be extracted from the ternary plot.

In addition to determining Pb and correlated trace metal provenance in riverine sediments, Pb isotopes have begun to be used for a wide range of other purposes, including the assessment of Pb loadings in both the dissolved and particulate form in river waters [133,166], Pb cycling, such as is involved with sediment/rock/water interactions [177,178], and the source of Pb in aquatic and terrestrial biota [147]. The use of Pb isotopes for these purposes is likely to expand significantly in the future.

4.3. Other Stable Isotopes as Geochemical Tracers

In addition to Pb, a number of other isotopes such as those of Nd, Sm, and Sr have been widely used by the geological community to determine the provenance of sediments in marine and riverine environments [109,179-182]. However, prior to the 1990s, the precision of determining the isotopic content of geological materials for many trace metal and metalloid contaminants was limited as a result of relatively small mass differences between the isotopes. This lack of precision in isotopic analysis inhibited their use as geochemical tracers as the variation in the measured isotopic composition within a given material was similar to that observed between geological materials. Recent advances in analytical capabilities have allowed for the precise determination of the relative isotopic abundances of elements with an atomic mass greater than 40 amu. These newly acquired methods have opened up the potential to utilize a wide range of non-traditional stable isotopes as geochemical fingerprints and tracers in forensic environmental studies. Given the novelty of the methods, only a handful of studies have examined the potential to effectively apply these isotopes as geochemical tracers of anthropogenic contaminants, and even fewer studies have focused on their potential use in riverine sedimentary environments. Here, we focus on six trace metal and metalloid contaminants (Cd, Cu, Cr, Hg, Sb, and Zn) that show some potential for use in river systems (Table 2).

Table 2. Summary of Selected Elemental Isotope Characteristics and Use as a Geochemical Fingerprint and Tracer.

Element Stable Isotopes * Frequently Used Isotopic Ratios 1 Primary Fractionation Processes Comments Selected References Relevant to Riverine Environments

106Cd, 108Cd, Evaporation and condensation, including Possible use to assess nutrient uptake and Gao etal. [183]; Cloquet etal. [184];

Cd 110Cd, 111Cd, 114Cd/110Cd that associated with ore processing; translocation in microorganisms and plants; Weiss et al. [185];

112Cd, 113Cd, 114Cd, 116Cd biological processes may be used to trace source, but not adequately tested Shiel et al. [186,187]

Redox reaction; biogenic accumulation; May be used to assess redox processes Weiss et al. [185];

Cu 63Cu, 65Cu 65Cu/63Cu solution speciation during biogeochemical cycling and to trace source of Cu, but not adequately tested Bird [29]; Kimball et al. [188]; Zhu et al. [189]; Petit et al. [190]

Cr 50Cr, 52Cr, 53Cr/52Cr Redox reactions; equilibrium adsorption May be used to assess redox processes Weiss et al. [185]

53Cr, 54Cr during biogeochemical cycling of elements

196Hg, 198Hg, MDF—Redox reactions, biomethylation, May be used to assess redox processes Mil-Homes et al. [105];

Hg 199Hg, 200Hg, 202Hg/198Hg evaporation, condensation; absorption; during biogeochemical cycling of elements Bergquist and Blum [191];

201Hg, 202Hg, 204Hg MIF—photoreduction, biological transformations Yin et al. [192]; Sonke et al. [193]

204Pb, 206Pb, 207Pb, 208Pb 206Pb/207Pb Negligible; Isotopic composition dependent Extensively utilized to assess Pb sources Bird [29] Villarroel et al. [45];

Pb 206Pb/208Pb 207Pb/208Pb on initial geological source Miller et al. [64,172]; Cheng [131]

Redox reactions; biological activity; Possible use to trace Sb sources; may be Rouxel et al. [194]; Wilson [195]

Sb 121Sb, 123Sb 121Sb/123Sb volatilization used to assess redox processes during biogeochemical cycling of elements

Evaporation and condensation, including Use appears promising for tracing Zn Bird [29]; Sivry et al. [196];

64Zn, 66Zn, that associated with ore processing; sources from industrial processes, Shiel et al. [187];

Zn 67Zn, 68Zn, 66Zn/64Zn biogenic uptake; diffusion, particularly from ore processing; may be Bentahila et al. [197];

70Zn adsorption/complexation on organics and Fe oxides used to assess of Zn cycling in environment Petit et al. [190]

Note: 1 Frequently used within environmental studies.

4.3.1. Zn and Cd Isotopes

Two important characteristics of Pb is that the Pb isotopic composition of an anthropogenic substance depends upon the ore deposits from which the Pb was derived, and the Pb isotopic signature of metallic ores differ from most other (geogenic) rocks and minerals. In addition, because mass-dependent fractionation of Pb by industrial or biological processes is negligible the measured Pb isotopic ratios within alluvial sediments result exclusively from the mixing of particles with different isotopic abundances. In contrast, variations in the isotopic abundances of Cd and Zn tend to be limited within natural geological materials (e.g., rocks, minerals, and ore deposits). However, measurable, but small variations in the isotopic composition of a material derived from the original rocks, minerals, ores, etc., may occur as a result of isotopic fractionation processes (Table 2) [184,186,198,199]. The created isotopic variations may then be utilized as a geochemical tracer to determine contaminant sources [199,200].

Interest in using the non-traditional isotopes of Cd and Zn as tracers is partly related to the desire to directly track the transport and fate of these potentially toxic trace metals in environmental media, and the recognition they may provide information that cannot be derived from the use of Pb isotopic fingerprinting methods alone. For example, since the Pb isotopic composition of a material depends on the initial source of the Pb, it may not be possible to separate Pb coming from different industrial processes if the utilized Pb is from the same ore deposit [131]. Geagea [201] found, for instance, that Pb from municipal waste incineration plants was indistinguishable from the Pb associated within a steel plant, inhibiting the use of Pb isotopes as a tracer of Pb provenance. In contrast, the isotopic composition of an anthropogenic substance with regards Cd and Zn depends on mass-dependent fractionation processes (e.g., evaporation and condensation). Thus, the isotope composition of industrial materials with regards to Cd and Zn are likely to differ, even if the materials are produced at the same location, as long as they are created by different processes.

The use of Cd, Zn, and other non-traditional stable isotopes as tracers, however, is not without its challenges, in part, because fractionation is produced by a number of biotic and abiotic processes that allow the isotopic composition of the source materials to change after its release into the environment [202]. Thus, the measured isotopic abundances within a given media (e.g., alluvial sediment) will not only reflect the mixing of sediments from different sources, but isotopic fractionation as it is dispersed and cycled through the environment. As a result, a sound understanding of elemental cycling within the system is required for their utilization, and their potential use on a large (watershed) scale may be limited.

Unlike the isotopic abundances of Pb which are expressed as a ratio, the isotopic abundance of the non-traditional stable isotopes including Cd and Zn are expressed per 1000 (mil, 5) or per 10,000 (s) relative to a standard, due to extremely small variations that occur in nature. Mathematically,

Ô'/ jX =

isotopic ratio of sample - isotopic ratio of standard

isotopic ratio of standard

el / jX =

isotopic ratio of sample - isotopic ratio of standard


isotopic ratio of standard

where S'ijX and £'/jX are the ó and e values of the element, respectively, for the isotopic ratio represented by i/j. In most instances, the lighter isotope is used as the denominator, such as


(often shown as 5114Cd), in which case, 5 > 0 is heavier than the standard, and 5 < 0 lighter.

To date, the use of Zn has received the most attention of the two elements, and shows the most promise as a tracer in riverine environments [29,196,203-206]. The mass-dependent fractionation of Zn is dominated by evaporation and condensation processes, although biogenic uptake, diffusion, and adsorption onto organic, Fe oxides, and other materials may also lead to fractionation (Table 2) [186,205]. Zn isotopes may therefore be used to tracer Zn derived from metal mining and refining operations, the steel industry, coal power plants, vehicle traffic, and urban waste incineration [204].

Due to the high relative abundances of 66Zn and 64Zn, and the ability to precisely measure their contents in geological materials, most studies have utilized 66Zn/64Zn ratios, reported as 566Zn. Unfortunately, studies to date have utilized a number of different materials as a standard, making the direct comparison of reported 566Zn values difficult [203].

With regards to the fingerprinting and tracing of Zn within river systems, one of the most detailed studies to date is provided by Silvry et al. [196] in an examination of Zn contamination within the Lot River basin, France by the industrial exploitation of Zn ore (the same river system described in the section above on the construction of pollution histories, Figure 5). Sivry et al. [196] found that the Zn isotopes within the ore (exhibiting 566Zn ~0.16%) was highly fractionated during metallurgical processing forming tailings materials with 566Zn values up to +1.49%. In addition, they were able to demonstrate that 566Zn values within sediments of the Riou Mort, a contaminated tributary, were significantly different than sediments found upstream of contaminant influx and the assumed regional background materials (566Zn values of +0.91+/-0.04% as compared to +0.31+/-0.06%). The changes in 566Zn values were consistent with the influx of tailings materials via the Riou Mort.

Sivry et al. [196] also found that systematic changes in 566Zn values were present within a dated sediment core extracted from a downstream hydroelectric reservoir (Figure 13). Sediments deposited between approximately 1952 and 1972 exhibited a mean 566Zn value of +0.95+/-0.08%, and are thought to be representative of the metallurgical processes used for Zn extraction at the time. Sediments deposited during the late 1970s exhibit increasing 566Zn values until reaching a maximum value in 1986, after which they remain relatively constant until the mid-1990s (Figure 13). The heavier 566Zn values were shown to correspond to a change in ore processing techniques and increases in Zn extraction efficiencies. Thus, it is likely that sediments within the reservoir contain a record of inputs to the river by different industrial processing methods.

As noted above, a limitation of using Zn isotopes to assess pollution histories is that the Zn stored within the sediments may undergo fractionation after deposition. Weiss et al. [202], for example, showed that diffusion of dissolved Zn in ombrotrophic peat bogs altered its isotopic composition, resulting in higher 566Zn values at the bottom of collected cores, exactly where diffusion related Zn fractionation was expected to be most pronounced. In spite of the potential for post-depositional fractionation, the work by Silvry et al. [196] suggests that it may not be a significant problem in all depositional environments.

Figure 13. Vertical (depth/time) variations in total Zn concentrations (red squares), Zn/SC ratios (blue circles), and S66Zn values (brown diamonds) measured in samples from a core extracted from the Cajarc reservoir along the Lot River, France (figure from Sivry et al. [196]; Zn and Zn/Sc variations from Audrey et al. [101]). The correlation between Zn concentration and isotopic trends suggest that Zn isotopes can be used to determine Zn source.

[Zn] (mg/kg) 0 2000 6000

S60Zn.iMC %o

Depth (cm)

10000 0.8

400 800 Zn/Sc


Draining water treatment Smelting waste confined

"1986 Accidental seepage release

1976-77change of waste-treatment process



Dates (year)

Although our focus in this review is on the use of geochemical tracers to identify trace metal sources in sediments, it is worth noting that Zn isotopes hold considerable promise in tracing anthropogenic contaminants in river waters [204-206]. Chen et al. [204], for example, found that spatial and temporal variations in Zn concentrations and enrichment factors varied semi-systematically with 566Zn values, showing that it was possible to use Zn isotopes as an indicator of anthropogenic contamination.

In comparison to Zn, much less work has been conducted on the use of Cd isotopes as a tracer of Cd provenance within alluvial sediments. In fact, with the exception of Gao et al. [183], the application of Cd isotopes as a tracer in riverine sedimentary systems seems to be lacking at the present time. However, recent studies have shown that Cd isotopes may prove effective at assessing Cd sources and cycling dynamics in river waters [207]. Like Zn, evaporation and condensation are the primary processes involved in the mass-dependent fractionation of Cd isotopes (Table 2). Thus, significant Cd fractionation is likely to accompany industrial evaporation and condensation processes, an assumption supported by the work of Cloquet et al. [184,208]. Other analysts suggest that fractionation also is associated with biological processes [185,209].

While generally small, the degree of fractionation, particularly by metallurgical processes (smelting, refining, etc.), appears to be significant enough to use Cd isotopes to identify specific anthropogenic Cd sources. For example, Cloquet et al. [184] argued that systematic variations in ô114/110Cd (ô114Cd) values in top soil near a smelter in northeast France could be used to trace Cd provenance. Similarly, Shiel et al. [187] suggested that differences in ô114Cd values measured in

marine bivalues was related to variations in the relative contribution of Cd from natural and anthropogenic sources.

4.3.2. Cu and Cr Isotopes

In contrast to Cd and Zn which are stable at the Earth's surface in a single valence state, Cu and Cr are stable in two oxidative states [210]. As a result, redox reactions which are not involved in the fractionation of Cd and Zn play a key role in the isotopic fractionation of both Cu and Cr. In the case of Cu, a number of studies have shown that the abiotic oxidative dissolution of Cu sulfide minerals, in which Cu(I) is the common form, will result in an aqueous Cu(II) product being heavier than in the original mineral [188,211]. In the presence of certain microbes, however, the fractionation between aqueous Cu and mineral Cu appears to decrease, apparently as a result of the preferential accumulation of the heavier 65Cu with the bacterial cell wall [210]. In the case of Cr, the reduction of highly toxic and water-soluble Cr(VI) by both abiotic and biotic processes appears to result in relatively insoluble Cr(III) that is approximately 3% to 5%o lighter than the remaining Cr(VI) [185,210,212,213].

The oxidation of Cu sulfide minerals is a common process along many mine contaminated river systems resulting in acid mine drainage. As a result, variations in the Cu isotopic values observed between sulfide minerals, river waters, and secondary mineral precipitates may be useful for assessing various biogeochemical processes that influence the dispersal and cycling of Cu within mine contaminated riverine environments [29,214,215]. The use of Cu isotopes to tracer to source of Cu in rivers is less clear. Their use in the sedimentary components of the river may be limited because: (1) Cu isotopic variations within a deposit are often greater than those observed between deposits [216]; and (2) redox reactions may alter Cu isotopic ratios both during dispersal and after deposition. Nonetheless, Bird [29] argued that while the use of Cu isotopes was not without its challenges, they have potential to effectively determine Cu provenance.

As for Cu the use of Cr isotopes to determine contaminant provenance in rivers has yet to be adequately demonstrated. However, isotopic contrasts between insoluble Cr(III) and Cr(VI) in the coexisting waters should allow for Cr isotopes to distinguish between natural and anthropogenic Cr in aquatic systems [210]. In fact, Cr isotopes have been used to identify anthropogenic sources and to assess Cr migration processes (e.g., dilution and reduction) in groundwater [217-219].

4.3.3. Mercury Isotopes

Hg is considered a priority pollutant on a global scale and once released into the environment, its biogeochemical cycling is highly complex, owing in part to the fact that it possesses a stable gaseous form (Hg0), and may undergo a wide range of redox, phase, and biologically mediated transformations. In fact, methylmercury, the form predominantly accumulated in biota, is primarily produced by sulfate reducing microbes. Until recently, understanding the primary sources and biogeochemical cycling of Hg in a particular system relied on an analysis of total Hg concentrations within the media of interest, its chemical speciation in those media, and/or the application of mass balance and geochemical models [193]. It is now clear, however, that numerous inorganic and organic reactions result in Hg isotopic fractionation (Table 2), producing large variations in the Hg isotopic composition of natural and anthropogenic materials (Figure 14) [191]. These variations in Hg isotopic values may be

used to assess Hg source(s), pollution histories, and transformations during biogeochemical cycling [105,191-193,220]. Thus, significant effort has recently been devoted to determining the mechanisms of Hg fractionation, and how Hg isotopes may be used to assess Hg transport and fate in atmospheric, aquatic, and terrestrial systems.

Figure 14. Summary of published Hg isotopic data (from Bergquist and Blum, [191]). The horizontal bars represent range in mass-dependent fractionation (A) and mass-independent fractionation (B). The vertical bars represent the estimated crust- and mantel-derived Hg. Data sources: 1, Hintelmann and Lue [221], Smith et al. [222], Smith et al. [223]; 2, Smith et al. [223]; 3, Lauretta et al. [224]; 4, Sonke et al. [225], Zambardi et al. [226];

5, Jackson et al. [227], Foucher and Hintelmann [228], Biswas et al. [229], Ghosh et al. [230], Jackson et al. [231], Fourcher et al. [232], Gehrke et al. [233];

6, Biswas et al. [229], Gosh et al. [230], Carignan et al. [234], Bergquist et al. [191].

There are seven stable isotopes of Hg. Variations in Hg isotopic values within a given media results from both mass-dependent and mass-independent processes (Table 2). Mass-independent fractionation of stable isotopes is highly unusual and, in the case of Hg has only been observed to affect the odd isotopes (199Hg, 201Hg). Thus, the ratio of 199Hg to Hg may be particularly useful for characterizing chemical pathways, such as photochemical reduction, associated with mass-independent fractionation processes [191].

Mass-dependent isotopic compositions are reported using 5 notation and the NIST Hg standard 3133 expressed mathematically as:

Ô202 X-

( 202 Hg /198 Hg)


(202 Hg / 198Hg )

-1 x(1000)

Mass-independent fractionation of the odd isotopes is reported as the difference between the measured ô199Hg and ô201Hg values to the predicted mass-dependent values as follows [235]:

A201 Hg = A199 Hg -(S202Hg x 0.2520)

A201 Hg = A201 Hg - (S202Hg x 0.2520)

As was the case for Cd and Zn, studies have suggested that metallurgical processes will result in waste products with distinct Hg isotopic signatures. As a result, the initial studies of Hg isotopes suggest that their analysis within samples from sedimentary deposits, including those within riverine environments, may be used to assess both Hg sources and loading histories [105,193,232,236-238]. Sonke et al. [193], for example, found that within the Lot River basin, France metal refining resulted in slag residues that exhibited heavier Hg isotopic values, and both mass-dependent and mass-independent Hg isotopic signatures of contaminated sediments were significantly different (heavier) than that of the local background materials. Once again a primary concern in the use of Hg isotopes within riverine environments is whether the contaminated source sediments will retain the original signature during transport and/or following deposition as a wide range of fractionation processes are available to alter its isotopic composition (Table 2) [237]. However, the work by Sonke et al. [193] suggests that Hg isotopes may behave conservatively in some sedimentary environments. Hg isotopes appear, then, to have considerable potential as a tracer of contaminant sources, transport processes, loading histories, and biogeochemical processes in river systems.

4.3.4. Sb Isotopes

Sb has two stable isotopes, 121Sb and 123Sb. Sb isotopic fractionation primarily occurs by redox changes of Sb(V) and Sb(III), the primary species found in aquatic environments, as well as by biological activity and volatilization [194]. The use of Sb isotopes as a geologic tracer has been limited to date, in part because Sb concentrations are typically low in many natural and anthropogenic materials, and the measurement of Sb isotopes is somewhat challenging. However, Rouxel et al. [194] demonstrated that precise measurements can be obtained for sea water and a variety of other geological materials using a MC-ICP-MS when Sb concentrations are at or above background values (Figure 15). Moreover, they found that the Sb isotopic composition of a variety of the studied materials varied significantly (Figure 15), leading them to suggest that Sb isotopes may serve as extremely useful tracers for the analysis of natural processes.

The only known application of Sb isotopes to riverine environments was conducted by Wilson [195] within the Rio Loa basin of northern Chile. The Rio Loa is extensively contaminated by Cu, As, Sb, and a variety of other trace metals from Cu mining and refining operations as wells by El Tatio, a natural geyser basin. Contaminated sediments from mining and refining operations and the geyser field were found to have similar Pb isotopic values, inhibiting the use of Pb as a tracer. Thus, Wilson [195] attempted to use Sb isotopes, which were hypothesized to serve as an effective tracer of sediment from geyser basin because: (1) the discharge of near boiling, silica-rich waters at El Tatio contain some of the highest reported surface water concentrations of Sb in the world; (2) the discharged Sb is incorporated at high concentrations (up to as much as 2 wt %) as Sb2O3 in silica based geyserite precipitates, and therefore should behave conservatively in sediment; and (3) concentrations

of Sb are relatively low in the local bedrock and the Cu mine related materials (ore, tailings, etc.). Wilson [195] found that the Sb isotopic composition of sediments from the geyser field varied significantly, and therefore like Pb, could not be used to distinguish between Sb from the mines and Sb from the gesyers. However, the Sb isotopic composition of the mines and geyser field differed significantly from Sb associated with other geological materials in the basin (e.g., rocks, sediments, soils), suggesting that Sb may serve as an effective tracer in mining wastes in some riverine environments.

Figure 15. Summary of s123Sb values measured in various natural materials. Variations indicate that Sb isotopes may be used as a contaminant tracer in geological media (from Rouxel et al. [194]).

5. Summary and Concluding Remarks

Environmental forensic investigations have increasingly relied on geochemical tracers to determine contaminant provenance, transport, and fate in contaminated rivers. Historically, contaminant fingerprints/tracers for point-source pollutants (e.g., trace metals from mining and refining operations) were based on bulk elemental concentrations, or related ratios, of the contaminants of interest. Provenance was determined on the basis of the geographical patterns of tracer concentration along the channel in comparison to background values. More recent investigations have combined tracer concentration data with sediment mixing models to not only determine contaminant source, but the relatively contribution of contaminants from a specific source or tributary basin. Although the modeling results are highly useful, they yield non-unique results, which when combined with the assumptions required for their use, lead to uncertainties in estimated source contributions. New

innovative statistical methods are needed to more accurately assess relative contaminant contributions when multiple natural and anthropogenic sources of trace metal contaminants exist.

Stable Pb isotopes have been shown to be particularly sensitive tracers that are highly effective at determining both the source and source contributions of Pb in rivers. Pb isotopes are most effectively used where there is a dominant source of Pb in the basin in addition to the Pb from the underlying bedrock (as is often the case in rivers contaminated by mining and/or refining operations). Where multiple Pb sources exist, fingerprinting the source materials may prove more difficult, and source contributions (when more than three sources exist) will need to rely on sediment mixing/unmixing models.

Advances in analytical chemistry, particularly the advent of the MC-ICP-MS, have created the possibility that a wide range of non-traditional stable isotopes can be used as environmental tracers, including Cd, Cu, Cr, Hg, Sb, and Zn. In contrast to Pb, differences in the isotopic composition of geological and anthropogenic materials are due to physical and chemical isotopic fractionation processes. As a result, it may not only be possible to determine a contaminant source on the basis of these isotopes, but to link trace metal concentrations to specific industrial processes used during a given time period at a single source location. In addition, some isotopic systems (e.g., Cr and Cu which are redox sensitivity) may be used to gain insights into biogeochemical cycling with river systems. The use of these non-traditional isotopes, however, is complicated by the potential for isotopic fractionation during contaminated dispersal or following deposition. Detailed information regarding the cycling of the utilized isotopes is therefore required for their effective use. Based on a relatively limited number of studies, the isotopes of Zn and Hg appear to hold the most promise in riverine environments for determining trace metal provenance in sediments or, in the case of Hg, assessing the biogeochemical cycling of trace metals in the aquatic environment. There is a dramatic need for additional studies that explore the use of these non-traditional stable isotope systems as tracers in riverine environments.

Innovative environmental forensic investigations in the 21st century will be characterized by a highly multidisciplinary approach that: (1) combines geochemical tracer studies with detailed stratigraphic, geomorphic, and hydrologic data; and (2) utilizes multiple isotopic tracers to more fully understand the temporal and spatial changes in trace metal source and source contributions along the river as well as the biogeochemical cycling of the metals within the aquatic system. In doing so, the developed conclusions will be based on multiple lines of evidence that are likely to withstand the scrutiny of judicial review.


The manuscript was significantly improved by the helpful comments of Suzanne Orbock Miller. Her efforts are greatly appreciated. The manuscript also benefited greatly from comments provided by four anonymous reviewers. All shortcomings inherent within the text are, of course, mine and do not reflect statements provide during the review process.

Conflict of Interest

The author declares no conflict of interest.


1. Murray, R.C. Evidence from the Earth; Mountain Press: Sevierville, TN, USA, 2004.

2. Murray, R.C. Forensic geology: Yesterday, today and tomorrow. In Forensic Geosciences: Principles, Techniques, and Applications; Special Publication 232; Pye, K., Croft, D.J., Eds.; The Geological Society of London: London, UK, 2004; pp. 7-9.

3. Zarull, M.A.; Hartig, J.H.; Maynard, L.; Sediment Priority Action Committee; Great Lakes Water Quality Board. Ecological Benefits of Contaminated Sediment Remediation in the Great Lakes Basin; International Joint Commission: Detroit, MI, USA, 1999. Available online: files/publications/C80.pdf (accessed on 6 May 2013).

4. Miller, J.R.; Orbock Miller, S.M. Contaminated River: A Geomorphological-Geochemical Approach to Site Assessment and Remediation; Springer: Berlin, Germany, 2007.

5. Forstner, U.; Wittmann, G.T.W. Metal Pollution in the Aquatic Environment, 2nd ed.; Springer Verlag: New York, NY, USA, 1981.

6. Gibbs, R.J. Transport phases of transition metals in the Amazon and Yukon Rivers. Geol. Soc. Am. Bull. 1977, 88, 829-843.

7. Horowitz, A.J. A Primer on Sediment-Trace Element Chemistry, 2nd ed.; Lewis Publishers: Chelsea, MI, USA, 1991.

8. Macklin, M.G.; Brewer, P.A.; Hudson-Edwards, K.A.; Bird, G.; Coulthard, T.J.; Dennis, I.A.; Lechler, P.J.; Miller, J.R.; Turner, J.N. A Geomorphological-Geochemical Approach to River Basin Management in Mining-Affected Rivers. In The Human Role in Changing Fluvial Systems, Proceedings of the 37th International Binghamton Geomorphology Symposium, Columbia, SC, USA, 20-22 October 2006; James, L.A., Marcus, W.A., Eds.; Elsevier: Amsterdam, The Netherlands, 2006; pp. 423-447.

9. Foster, I.D.L.; Lees, J.A. Tracers in geomorphology: Theory and applications in tracing fine particulate sediments. In Tracers in Geomorphology; Foster, I.D.L., Ed.; John Wiley and Sons: Hoboken, NJ, USA, 2000; pp. 3-20.

10. Hoefs, J. Geochemical fingerprints: A critical appraisal. Eur. J. Miner. 2010, 22, 3-15.

11. Bigham, J.M.; Ciolkosz, E.J. Soil Color; SSSA Special Publication No. 31; Soil Science Society of America: Madison, WI, USA, 1993.

12. Croft, D.J.; Pye, K. Colour Theory and the evaluation of an instrumental method of measurement using geological samples for forensic applications. In Forensic Geosciences: Principles, Techniques, and Applications; Special Publication 232; Pye, K., Croft, D.J., Eds.; The Geological Society of London: London, UK, 2004; pp. 49-62.

13. Dudley, R.J.; Smalldon, K.W. The objective comparison of the particle size distribution in soils with particular reference to the sand fraction. Med. Sci. Law. 1978, 18, 278-281.

14. Robertson, J.; Thomas, C.J.; Caddy, B.; Lewis, A.J.M. Particle size analysis of soils: A comparison of dry and wet sieving techniques. Forensic Sci. Int. 1984, 24, 209-217.

15. Knox, J.C. Historical valley floor sedimentation in the upper Mississippi valley. Ann. Assoc. Am. Geogr. 1987, 77, 224-244.

16. Sutherland, R.A. Selective erosion and sediment source identification, Baringo District, Kenya. Z. Geomorphol. 1991, 35, 293-304.

17. Blott, S.J.; Croft, D.J.; Pye, K.; Saye, S.E.; Wilson, H.E.. Particle size analysis by laser diffraction. In Forensic Geosciences: Principles, Techniques, and Applications; Special Publication 232; Pye, K., Croft, D.J., Eds.; The Geological Society of London: London, UK, 2004; pp. 63-73.

18. Fan, P.F. Recent silts in the Santa Clara River drainage basin, Southern California: A mineralogical investigation of their origin and evolution. J. Sediment. Petrol. 1976, 46, 802-812.

19. Graves, W.J. A mineralogical soil classification technique for the forensic scientist. J. Forensic Sci. 1979, 24, 323-337.

20. Ugolini, F.C.; Corti, G.; Agenelli, A.; Piccardi, F. Mineralogical, physical and chemical properties of rock fragments in soil. Soil Sci. 1996, 161, 521-542.

21. Oldfield, R.; Rummery, T.A.; Thomas, R.; Walling, D. Identification of suspended sediment sources by means of magnetic measurements: Some preliminary results. Water Resour. Res. 1979, 15, 211-218.

22. Walling, D.E.; Peart, M.R.; Oldfield, F.; Thompson, R. Suspended sediment sources identified by magnetic measurements. Nature 1979, 281, 110-113.

23. Fitzpatrick, F.; Thornton, J.I.M. Wide-band detector for microampere low-energy elecontron currents. J. Forensic Sci. 1974, 4, 460-475.

24. Demmelmeyer, H.; Adam, J. Forensic investigation of soil and vegetable materials. Forensic Sci. Rev. 1995, 7, 120-136.

25. Pirrie, D.; Butcher, A.R.; Power, M.R.; Gottlieb, P.; Miller, G.L. Rapid quantitative mineral and phase analysis using automated scanning electron microscopy (QuemSCAN): Potential applications in forensic geosciences. In Forensic Geosciences: Principles, Techniques, and Applications; Special Publication 232; Pye, K., Croft, D.J., Eds.; The Geological Society of London: London, UK, 2004; pp. 123-136.

26. Pye, K. Forensic Examination of rocks, sediments, soils, and dust using scanning electron microscopy and X-ray chemical microanalysis. In Forensic Geosciences: Principles, Techniques, and Applications; Special Publication 232; Pye, K., Croft, D.J., Eds.; The Geological Society of London: London, UK, 2004; pp. 103-122.

27. Lewin, J.; Wolfenden, P.J. The assessment of sediment source: A field experiment. Earth Surf. Proc. Land. 1978, 3, 171-178.

28. Knox, J.C. Rates of floodplain overbank vertical accretion. In Floodplain Evolution; Brakenridge, G.R., Hagedorn, J., Eds.; Elsevier: Amsterdam, The Netherlands, 1992.

29. Bird, G. Provenancing anthropogenic Pb within the fluvial environment: Developments and challenges in the use of Pb isotopes. Environ. Int. 2011, 37, 802-819.

30. Macklin, M.G. Floodplain sedimentation in the Upper Axe Valley, Mendip, England. Trans. Inst. Br. Geogr. 1985, 10, 235-244.

31. Passmore, D.G.; Macklin, M.G. Provenance of fine-grained alluvium and late holocene land-use change in the Tyne Basin, Northern England. Geomorphology 1994, 9, 127-142.

32. Church, S.E.; Unruh, D.M.; Fey, D.L.; Sole, T.C. Trace Element and Lead Isotopes in Streambed Sediment in Streams Affected by Historical Mining; U.S. Geological Survey Professional Paper 1652; U.S. Government Printing Office: Washington, DC, USA, 2004; Chapter D8.

33. Ulrich, A.; Moor, C.; Vonmont, H.; Jordi, H.R.; Lory, M. ICP-MS trace element analysis as a forensic tool. Anal. Bioanal. Chem. 2004, 378, 1059-1068.

34. Aggarwal, J.; Habicht-Mauche, J.; Juarez, C. Application of heavy stable isotopes in forensic geochemistry: A review. Appl. Geochem. 2008, 23, 2658-2666.

35. Zhang, H.; Yao, Q.; Zhu, Y.; Fan, S.; He, P. Review of source identification methodologies for heavy metals in solid waste. Chin. Sci. Bull. 2013, 58, 162-168.

36. Collins, A.L.; Walling, D.E.; Leeks, G.J.L. Source type ascription for fluvial suspended sediment based on a quantitative composite fingerprinting technique. Catena 1997, 29, 1-27.

37. Collins, A.L.; Walling, D.E.; Leeks, G.J.L. Use of the geochemical record preserved in floodplain deposits to reconstruct recent changes in River Basin Sediment Sources. Geomorphology 1997, 19, 151-167.

38. Collins, A.L.; Walling, D.E.; Leeks, G.J.L. Use of composite fingerprints to determine the provenance of the contemporary suspended sediment load transported by rivers. Earth Surf. Proc. Land. 1998, 23, 31-52.

39. Collins, A.L.; Walling, D.E. Selecting fingerprinting properties for discriminating potential suspended sediment sources in river basins. J. Hydrol. 2002, 261, 218-244.

40. Miller, J.; Lord, M.; Yurkovich, S.; Mackin, G.; Kolenbrander, L. Historical trends in sedimentation rates and sediment provenance, Fairfield Lake, western North Carolina. Water Resour. Bull. 2005, 41, 1053-1075.

41. Miller, J.R.; Mackin, G.; Lechler, P.; Lord, M.; Lorentz, S. Influence of basin connectivity on sediment source, transport, and storage within the Mkabela Basin, South Africa. Hydrol. Earth Syst. Sci. 2012, 17, 1-22.

42. Cheng, H.; Hu, Y. Lead (Pb) isotopic fingerprinting and its applications in lead pollution studies in China: A review. Environ Pollut. 2010, 158, 1134-1146.

43. Farmer, J.G.; Eades, L.J.; MacKenzie, A.B.; Kirika, A.; Bailey-Watts, T.E. Stable lead isotope record of lead pollution in Loch Lomond sediments since 1630 A.D. Environ. Sci. Technol. 1996,

3, 3080-3083.

44. Duzgoren-Aydin, N.S.; Weiss, A.L. Use and abuse of Pb-isotope fingerprinting technique and GIS mapping data to assess lead in environmental studies. Environ. Geochem. Health 2008, 30, 577-588.

45. Villarroel, L.F.; Miller, J.R.; Lechler, P.J.; Germanoski, D. Lead, Zinc, and Antimony contamination of the Rio Chilco-Rio Tupiza drainage system, southern Bolivia. Environ. Geol. 2006, 51, 283-299.

46. Wolfenden, P.J.; Lewin, J. Distribution of metal pollutants in floodplain sediments. Catena 1978,

4, 309-317.

47. Lewin, J.; Macklin, M.G. Metal mining and floodplain sedimentation in Britian. In International Geomorphology; Gardiner, V., Ed.; John Wiley and Sons: Hoboken, NJ, USA, 1987; pp. 1009-1027.

48. Marcus, W.A. Copper dispersion in ephemeral stream sediments. Earth Surf. Proc. Land. 1987, 12, 217-228.

49. Taylor, M.P.; Kesterton, R.G.H. Heavy metal contamination of an arid river environment: Gruben River, Namibia. Geomorphology 2002, 42, 311-327.

50. Miller, J.R.; Lechler, P.J.; Desilets, M. The role of geomorphic processing in the transport and fate of mercury in the Carson River basin, west-central Nevada. Environ. Geol. 1998, 33, 249-262.

51. Moore, J.N.; Brook, E.J.; Johns, C. Grain size partitioning of metals in contaminated, coarse-grained river floodplain sediments: Clark Fork River, Montana, USA. Environ. Geol. Water S. 1989, 12, 107-115.

52. Loring, D.H. Lithium—A new approach for the granulometric normalization of trace metal data. Mar. Chem. 1990, 29, 155-168.

53. Loring, D.H.; Rantala, R.T.T. Manual for the geochemical analysis of sediments and suspended particular matter. Earth Sci. Rev. 1992, 32, 235-283.

54. Devesa-Rey, R.; Díaz-Fierros, F.; Barral, M.T. Normalization strategies for river bed sediments: A graphical approach. Microchem. J. 2009, 91, 253-265.

55. Macklin, M.G.; Lewin, J. Sediment transfer and transformation of an alluvial valley floor: The river south Tyne, Northumbria, U.K. Earth. Surf. Proc. Land. 1989, 14, 233-246.

56. Macklin, M.G.; Dowsett, R.B. The chemical and physical speciation of trace metals in fine grained overbank flood sediments in the Tyne Basin, north-east England. Catena 1989, 16, 135-151.

57. Mackin, M.G. Fluxes and storage of sediment-associated heavy metals in floodplain systems: assessment and river basin management issues at a time of rapid environmental change. In Floodplain Processes; Anderson, M.G., Walling, D.E., Bates, P.D., Eds.; John Wiley and Sons: Hoboken, NJ, USA, 1996; pp. 441-460.

58. Miller, J.R.; Germanoski, D.; Villarroel, L.F.; Lechler, P. Spatial and temporal variations in the transport and storage of trace metal contaminants in the Upper Río Pilcomayo, Southern Bolivia. Int. J. Environ. Heal. R. 2009, 3, 334-361.

59. Nicholas, A.P.; Ashworth, P.J.; Kirkby, M.J.; Macklin, M.G.; Murray, T. Sediment slugs: Large-scale fluctuations in fluvial sediment transport rates and storage volumes. Prog. Phys. Geog. 1995, 19, 500-519.

60. Hoey, T. Temporal variation in bedload transport rates and sediment storage in gravel river beds. Prog. Phys. Geog. 1992, 16, 319-338.

61. Bartley, R.; Rutherfurd, I. Re-evaluation of the wave model as a tool for quantifying the geomorphic recovery potential of streams disturbed by sediment slugs. Geomorphology 2005, 64, 221-242.

62. Gilbert, G.K. Hydraulic-Mining Debris in the Sierra Nevada; U.S. Geological Survey Professional Paper 105; U.S. Government Printing Office: Washington, DC, USA, 1917.

63. Helgen, S.O.; Moore, J.N. Natural background determination and impact quantification in trace metal-contaminated river sediments. Environ. Sci. Technol. 1996, 30, 129-135.

64. Miller, J.R.; Lechler, P.J.; Hudson-Edwards, K.A.; Macklin, M.G. Lead isotopic fingerprinting of heavy metal contamination, Rio Pilcomayo, Bolivia. J. Geochem. Explor. Environ. Anal. 2002, 2, 225-233.

65. Yu, L.; Oldfield, F. A multivariate mixing model for identifying sediment source from magnetic measurements. Quat. Res. 1989, 32, 168-181.

66. Rowan, J.S.; Goodwill, P.; Franks, S.W. Uncertainty estimation in fingerprinting suspended sediment sources. In Tracers in Geomorphology; Foster, I.D.L., Ed.; John Wiley and Sons: Hoboken, NJ, USA, 2000; pp. 279-289.

67. Rowan, J.S.; Black, S.; Franks, S.W. Sediment fingerprinting as an environmental forensics tool explaining cyanobacteria blooms in lakes. Appl. Geogr. 2012, 32, 832-843.

68. Collins, A.L.; Walling, D.E.; Webb, L.; King, P. Apportioning catchment scale sediment sources using a modified composite fingerprinting technique incorporating property weights and prior information. Geoderma 2010, 155, 249-261.

69. Small, I.F.; Rowan, J.S.; Franks, S.W. Quantitative sediment fingerprinting using a Bayesian uncertainty estimation framework. In The Structure, Function and Management Implications of Fluvial Sedimentary Systems; Dyer, F.J., Thoms, M.C., Olley, J.M., Eds.; IAHS Publications: Wallingford, UK, 2002; pp. 443-450.

70. Small, I.F.; Rowan, J.S.; Franks, S.W.; Wyatt, A.; Duck, R.W. Bayesian sediment fingerprinting provides a robust tool for environmental forensic geosciences applications. In Forensic Geosciences: Principles, Techniques, and Applications; Special Publication 232; Pye, K., Croft, D.J., Eds.; The Geological Society of London: London, UK, 2004; pp. 207-213.

71. Miller, J.R.; Mackin, G. Concentrations, sources, and potential ecological impacts of selected trace metals on aquatic biota within the Little Tennessee River Basin, North Carolina. Water Air Soil Poll. 2013, in press.

72. Macklin, M.G.; Klimek, K. Dispersal, storage and transformation of metal-contaminated alluvium in the upper Vistula basin, southwest Poland. Appl. Geochem. 1992, 12, 7-30.

73. Hudson-Edwards, K.A.; Macklin, M.; Taylor, M. Historic metal mining input to Tees River sediment. Sci. Total Environ. 1997, 194-195, 437-445.

74. Nguygen, H.L.; Braun, M.; Szaloki, I.; Baeyens, W.; van Grieken, R.; Lerrmakers, M. Tracing the metal pollution history of the Tisza River through the analysis of a sediment depth profile. Water Air Soil Poll. 2009, 200, 119-132.

75. Bábek, O.; Faméra, M.; Hilscherová, K.; Kalvoda, J.; Dobrovolny, P.; Sedlácek, J.; Machát, J.; Holoubek, I. Geochemical traces of flood layers in the fluvial sedimentary archive: Implications for contamination history analyses. Catena 2011, 87, 281-290.

76. Matys Grygar, T.M.; Nováková, T.; Bábek, O.; Elznicová, J.; Vadinová, N. Robust assessment of moderate heavy metal contamination levels in floodplain sediments: A cast study on the Jizera River, Czech Republic. Sci. Total Environ. 2013, 452-453, 233-245.

77. Du Laing, G.; Rinklebe, J.; Vandencasteele, B.; Meers, E.; Tack, F.M.G. Trace metal behavior in estuarine and riverine floodplain sediments: A review. Sci. Total Environ. 2009, 407, 3972-3985.

78. Matys Grygar, T.; Novakoa, T.; Mihaljevic, M.; Strnad, L.; Koptikova, L. Surprisingly small increase in the sedimentation rate in the foodplain of Marova River in the Straznice area, Czech Republic, in the last 100 years. Catena 2011, 86, 461-471.

79. Matys Grygar, T.; Sedlacek, J.; Bábek, O.; Novakoa, T.; Strnad, L.; Mihaljevic, M. Regional contamination of Moravia (south-eastern Czech Republic): Temporal shift of Pb and Zn loading in fluvial sediments. Water Air Soil Pollut. 2012, 223, 739-753.

80. Ferrand, E.; Eyrolle, F.; Radakovitch, O.; Provansal, M.; Dufour, S.; Vella, C.; Raccasi, G.; Gurriaran, R. Historical levels of heavy metals and artificial radionuclides reconstructed from overbank sediment records in lower Rhône River (South-East France). Geochim. Cosmochim. Acta 2012, 82, 163-182.

81. Walling, D.E.; Owens, P.N.; Carter, J.; Leeks, G.J.L.; Lewis, S.; Meharg, A.A.; Wright, J. Storage of sediment-associated nutrients and contaminants in river channel and floodplain systems. Appl. Geochem. 2003, 18, 195-220.

82. Tobin, G.A.; Brinkmann, R.; Montz, B.E. Flooding and the distribution of selected metals in floodplain sediments in St. Maries, Idaho. Environ. Geochem. Health 2000, 22, 219-232.

83. Pease, P.; Leece, S.; Gares, P.; Rigsby, C. Heavy metal concentrations in sediment deposits on the Tary River floodplain following Hurricane Floyd. Environ. Geol. 2007, 51, 1003-1111.

84. Stokes, S.; Walling, D.E. Radiogenic and isotopic methods for the direct dating of fluvial sediments. In Tools in Fluvial Geomorphology; Kondolf, M.G., Piegay, H., Eds.; John Wiley and Sons: Chichester, UK, 2005; pp. 233-267.

85. Hudson-Edwards, K.A.; Macklin, M.G.; Curtis, C.D.; Vaughan, D.J. Chemical remobilization of contaminated metals within floodplain sediments in an incising river system: Implications for dating and chemostratigraphy. Earth Surf. Proc. Land. 1998, 23, 671-684.

86. Davies, B.E.; Lewin, J. Chronosequences in alluvial soils with special reference to historical pollution in Cardiganshire, Wales. Environ. Pollut. 1974, 6, 49-57.

87. Rang, M.C.; Schouten, C.J. Evidence for historical heavy metal pollution in floodplain soils: The Meuse. In Historical Change of Large Alluvial Rivers: Western Europe; Petts, G.E., Ed.; John Wiley and Sons: Chichester, UK, 1989; pp. 127-142.

88. Sear, D.; Carver, S. The release and dispersal of Pb and Zn contaminated sediments within an Arctic braided river system. Appl. Geochem. 1996, 11, 187-195.

89. Klimek, K.; Zawilinska, L. Trace elements in alluvia of the Upper Vistula as indicators of palaeohyology. Earth Surf. Proc. Land. 1985, 10, 273-280.

90. Taylor, M.P.; Lewin, J. River behaviour and Holocene alleviation: The river Severn at Welshpool, mid-Wales, U.K. Earth Surf. Proc. Land. 1996, 21, 77-91.

91. Harrison, J.; Heijnis, H.; Caprarelli, G. Historical pollution variability form abandoned mine sites, Greater Blue Mountains World Heritage Area, New South Wales, Australia. Environ. Geol. 2003, 43, 680-687.

92. Bradley, S.B.; Cox, J.J. Heavy metals in the Hamps and Manifold valleys, North Straffordshire, U.K.: Partitioning of metals in floodplain soils. Sci. Total Environ. 1987, 65, 135-153.

93. Knox, J.C. Responses of floods to Holocene climatic change in the Upper Mississippi Valley. Quatern. Res. 1985, 23, 287-300.

94. Berner, Z.A.; Bleeck-Schumidt, B.; Stüben, D.; Neumann, T.; Fuchs, M.; Lehmann, M. Floodplain deposits: A geochemical archive of flood history—A case study on the River Rhine, Germany. Appl. Geochem. 2012, 27, 543-561.

95. Navrâtil, T.; Rohovec, J.; Zâk, K. Floodplain sediments of the 2002 catastrophic flood at eh Vltava (Moldau) River and its tributaries: Mineralogy, chemical composition, and post-sedimentary evolution. Environ. Geol. 2008, 56, 399-412.

96. Brovo, A.G.; Loizeau, J.L.; Ancey, L.; Ungureanu, V.G.; Dominik, J. Historical record of mecury contamination in sediments from the Babeni Reservoir in the Olt River, Romania. Environ. Sci. Pollut. R. 2009, 16, 66-75.

97. Lewin, J.; Macklin, M. Perservation potential for Late Quaternary river alluvium. J. Quat. Sci. 2003, 18, 107-120.

98. Müller, J.; Ruppert, H.; Muramatsu, Y.; Schneider, J. Reservoir sediments—A witness of mining and industrial development (Malter reservoir, eastern Erzgebirge, Germany). Environ. Geol. 2000, 39, 1341-1351.

99. Arnason, J.G.; Fletcher, B.A. A 40+ year record of Cd, Hg, Pb, and U deposition in sediments of Patroon Reservoir, Albany County, NY, USA. Environ. Pollut. 2003, 123, 383-391.

100. Callender, E. Geochemical effects of rapid sedimentation in aquatic systems: Minimal diagenesis and the preservation of historical metal signatures. J. Paleolimnol. 2000, 23, 243-260.

101. Audry, S.; Schäfer, J.; Blanc, G.; Jouanneau, J.M. Fify-year sedimentary record of heavy metal pollution (CD, Zn, Cu, Pb) in the Lot River reservoirs (France). Environ .Pollut. 2004, 132, 413-426.

102. Von Gunten, H.R.; Sturm, M.; Moser, R.N. 200-year record of metals in lake sediments and natural background concentrations. Environ. Sci. Technol. 1997, 31, 2193-2197.

103. Kähkönen, M.A.; Suominen, K.P.; Manninen, P.K.; Salkinoja-Salonen, M.S. 100 years of sediment accumulation history of organic halogens and heavy metals in recipient and nonrecipient lakes of pulping industry in Finland. Environ. Sci. Technol. 1998, 32, 1741-1746.

104. Lima, A.L.; Bergquist, B.A.; Boyle, E.A.; Reuer, M.K.; Dudas, F.O.; Reddy, C.M.; Eglinton, T.I. High-resolution historical records from Pettaquamscutt River basin sediments. Geochim. Cosmochim. Acta 2005, 69, 1813-1824.

105. Mil-Homens, M.; Blum, J.; Canario, J.; Caetano, M.; Costa, A.M.; Lebreiro, S.M.; Trancoso, M.; Richter, T.; de Stigter, H.; Johnson, M.; et al. Tracing anthropogenic Hg and Pb input using stable Hg and Pb isotope ratios in sediments of the central Portuguese Margin. Chem. Geol. 2013, 336, 62-71.

106. MacDonald, D.D.; Ingersoll, C.G.; Berger, T.A. Development and evaluation of consensus-based sediment quality guidelines for freshwater ecosystems. Arch. Environ. Contaminat. Toxicol. 2000, 39, 20-31.

107. Leigh, D.S. Morphology and channel evolution of small streams in the Southern Blue Ridge Mountains of western North Carolina. Southeastern Geogr. 2010, 50, 397-421.

108. Banner, J.L. Radiogenic isotopes: Systematics and applications to earth surface processes and chemical stratigraphy. Earth Sci. Rev. 2004, 65, 141-194.

109. Calanchi, N.; Dinelli, E.; Lucchini, F.; Mordenti, A. Chemostratigraphy of late Quaternary sediments from Lake Albano and Central Adriatic Sea corres (Paliclas Project). Memoir 1st Ital. Idrobiol. 1996, 55, 247-263.

110. Boyle, J.F. Inorganic geochemical methods in palaeolimnology. In Tracking Environmental Change Using Lake Sediments: Physical and Geochemical Methods; Last, W.M., Smol, J.P., Eds., Khrwer: Dordrecht, The Netherland, 2001; pp. 85-141.

111. Lucchini, F.; Dinelli, E.; Calachi, N. Chemostratigraphy of Lago Albano sediments (central Italy): Geochemical evidence of palaeoenvironmental changes in late Quaternary. J. Paleoliminol. 2003, 29, 109-122.

112. Mabrouk, A.; Jarvis, I.; Belayouni, H.; Moody, R.T.J.; de Cabrera, S. An integrated chemostratigraphic study of the Campanian-Early Masstrichtian deposits of offshore Miskar Field in southeastern Tunisia: SIS, S13C and S18 O isotopes, and elemental geochemistry. Stratigraphy 2005, 2, 193-216.

113. Popp, C.L.; Hawley, J.W.; Love, D.W.; Dehn, M. Use of radiometric (Cs-137, Pb-210), geomorphic and stratigraphic techniques to date recent oxbow sediments in the Rio Puerco drainage, Grants Uranium region, New Mexico. Environ. Geol. Water Sci. 1988, 11, 253-269.

114. Chillrud, S.N.; Semming, S.; Shuster, E.L.; Simpson, H.J.; Bopp, R.F.; Ross, J.M.; Pederson, D.C.; Chaky, D.A.; Tolley, L.-R.; Estabrooks, F. Stable lead isotopes, contaminant metals and radionuclides in upper Hudson River sediment cores: Implications for improved time stratigraphy and transport processes. Chem. Geol. 2003, 199, 53-70.

115. Grygar, T.; Svetlik, I.; Lisa, L.; Koptikova, L.; Bajer, A.; Wray, D.; Ettler, V.; Mihaljevic, M.; Novakova, T.; Koubova, M.; et al. Geochemical tools for the stratigraphic correlation of floodplain deposits of the Morava River in Straznicke Pomoravi, Czech Republic from the last millennium. Catena 2010, 80, 106-121.

116. Miller, JR.; Hudson-Edwards, K.A.; Lechler, P.J.; Preston, D.; Macklin, M.G. Heavy metal contamination of water, soil, and produce within riverine communities of the Rio Pilcomayo basin, Bolivia. Sci. Total Environ. 2004, 320, 189-209.

117. Chow, T.J.; Johnstone, M.S. Lead isotopes in gasolina and aerosols of Los Angeles Basin, California. Science 1985, 147, 147-148.

118. Ault, W.A.; Senechal, R.G.; Erleback, W.E. Isotopic composition as a natural tracer of lead in the environment. Environ. Sci. Technol. 1970, 4, 305-313.

119. Chiaradia, M.; Cupelin, F. Behaviour of airborne lead and temporal variations of its source effects in Geneva (Switzerland): Comparison of anthropogenic versus natural processes. Atmos. Environ. 2000, 34, 959-971.

120. Gulson, B.L.; Howarth, D.; Mizon, K.J.; Waw, A.J.; Korsch, M.J.; Davis, J.J. The source of lead in humans from Broken Hill mining community. Environ. Geochem. Health 1994, 16, 19-25.

121. Gulson, B.L.; Mizon, K.J.; Korsch, M.J.; Howarth, D. Importance of monitoring family members in establishing sources and pathways of lead in blood. Sci. Total Environ. 1996, 188, 173-182.

122. Gulson, B.L.; Mizon, K.J.; Davis, J.D.; Palmer, J.M.; Vimpani, G. Identification of sources of led in children in a primary zinc-led smelter environment. Environ. Health Persp. 2004, 112, 52-60.

123. Rosman, K.J.R.; Chisolm, W.; Boutron, C.F.; Candelone, J.P.; Gorlach, U. Isotopic evidence for the source of lead in Greenland snows since the late 1960s. Nature 1993, 362, 333-334.

124. Rosman, K.J.R.; Chisolm, W.; Boutron, C.F.; Candelone, J.P.; Patterson, C.C. Anthropogenic lead isotopes in Antarctica. Geophy. Res. Lett. 1994, 21, 2669-2672.

125. Rosmann, K.J.R.; Chisolm, W.; Hong, S.; Candelone, J.P.; Boutron, C.F. Lead from Carthaginian and Roman Spanish mines isotopically identified in Greenland ice dated from 600 B.C. to 300 A.D. Environ. Sci. Technol. 1997, 31, 3413-3416.

126. Shotyk, W.; Zheng, J.C.; Krachler, M.; Zdanowicz, C.; Koerner, R.; Fisher, D. Predominance of industrial Pb in recent snow (1994-2004) and ice (1842-1996) from Devon Island, Arctic Canada. Geophys. Res. Lett. 2005, 32, L21814, doi:10.1029/2005GL023860.

127. Veysseyre, A.M.; Bollhofer, A.F.; Rosman, K.J.R.; Ferrari, C.P.; Boutron, C.F. Tracing the origin of pollution in French Alpine snow and aerosols using lead isotope ratios. Environ. Sci. Technol. 2001, 35, 4463-4469.

128. Gulson, B.L.; Tiller, K.G.; Mizon, K.J.; Merry, R.H. Use of lead isotope ratios in soils to identify the source of lead contamination near Adelaide, South Australia. Environ. Sci. Tech. 1981, 15, 691-696.

129. Steinmann, M.; Stille, P. Rare earth element behavior and Pb, Sr, Nd isotope systematics in a heavy metal contaminated soil. Appl. Geochem. 1997, 12, 607-632.

130. Hansmann, W.; Koppel, V. Lead-isotopes as tracers of pollutants in soils. Chem. Geol. 2000, 171, 123-144.

131. Mihaljevic, M.; Ettler, V.; Sebek, O.; Strnad, L.; Chrastny, V. Lead isotopic signatures of wine and vineyard soils—Tracers of lead origin. J. Geochem. Explor. 2006, 88, 130-133.

132. Reimann, C.; Flem, B.; Fabian, K.; Birke, M.; Landenberger, A.; Négrel, P.; Demetriades, A.; Hoogewerff, J.; The GEMAS Project Team. Lead and lead isotopes in agricultural soils of Europe—The continental perspective. Appl. Geochem. 2012, 27, 532-542.

133. Shirahata, H.; Elias, R.W.; Patterson, C. Chronological variations in concentrations and isotopic compositions and anthropogenic atmospheric Pb in sediments of a remote subalpine pond. Geochim. Cosmochim. Ac. 1980, 44, 149-167.

134. Petit, D.; Mennesier, J.P.; Lamberts, L. Stable Pb isotopes in pond sediments as tracer of past and present atmospheric Pb pollution in Belgium. Atmos. Environ. 1984, 18, 1189-1193.

135. Chiaradia, M.; Chenhall, B.E.; Depers, A.M.; Gulson, B.I.; Jones, B.G. Identification of historical lead sources in roof dusts and recent lake sediments from an industrialized area: Indications from lead isotopes. Sci. Total Environ. 1997, 2005, 107-128.

136. Foster, I.D.L.; Less, J.A.; Jones, A.R.; Chapman, A.S.; Turner, S.E.; Hodgkinson, R. The possible role of agricultural land drains in sediment delivery to a small reservoir, Worcestershire UK: A multiparameter tracing study. In Structure, Function and Management Implications of Fluvial Sedimentary Systems; Dyer, F.J., Thoms, M.C., Olley, J.M., Eds.; IAHS Publications: Wallingford, UK, 2002; pp. 443-442.

137. Shotyk, W.; Weiss, D.; Appleby, P.G.; Cheburkin, A.K.; Frei, R.; Gloor, M.; Kramers, J.D.; Reese, S.; Knaap, W.O.V.D. History of atmospheric lead deposition since 12,370 14C yr BP from a peat bog, Jura Mountains, Switzerland. Science 1998, 281, 1635-1640.

138. Weiss, D.; Shotyk, W.; Appleby, P.G.; Kramers, J.D.; Cheburkin, A.K. Atmospheric lead deposition since the industrial revolution recorded by five Swiss peat profiles: Enrichment factors, fluxes, isotopic composition, and sources. Environ. Sci. Technol. 1999, 33, 1340-1352.

139. Marcantonio, F.; Flowers, G.C.; Templin, N. Lead contamination in a wetland watershed: Isotopes as fingerprints of pollution. Environ. Geol. 2000, 39, 1070-1076.

140. Cloy, J.M.; Farmer, J.G.; Graham, M.C.; MacKenzie, A.B.; Cook, G.T. Historical records of atmospheric Pb deposition in four Scottish ombrotropic peat bogs: An isotopic comparison with other records from western Europe and Greenland. Glob. Biogeochem. Cycles 2008, 22, doi:10.1029/2007GB003059.

141. Bellis, D.J.; Satake, K.; McLeod, C.W. A comparison of lead isotope ratios in the bark pockets and annual rings of two beech trees collected in Derbyshire and South Yorkshire, UK. Sci. Total Environ. 2004, 32, 105-113.

142. Bindler, R.; Renberg, I.; Klaminder, J.; Emteryd, O. Tree rings as Pb pollution archives? A comparison of 206Pb/207Pb isotope ratios in pine and other environmental media. Sci. Total Environ. 2004, 319, 173-183.

143. Bi, X.; Feng, X.; Yang, Y.; Li, X.; Shin, G.P.; Li, F.; Qiu, G.; Li, G.; Liu, T.; Fu, Z. Allocation and source attribution of lead and dadmium in maize (Zea mays L.) impacted by smelting emissions. Environ. Pollut. 2009, 157, 834-939.

144. Manton, W.I. Sources of lead in blood—Identification by stable isotopes. Arch. Environ. Health 1977, 32, 149-159.

145. Keinonen, M. The isotopic composition of lead in man and the environment in Finland 1966-1987—Isotope ratios of lead as indicators of pollutant source. Sci. Total Environ. 1992, 113, 251-268.

146. Labonne, M.; Othamn, D.B.; Lee, J.M. Pb isotopes in mussels as tracers of metal sources and water movements in a lagoon (Thau Basin, S. France). Chem. Geol. 2001, 181, 181-191.

147. Miller, JR.; Anderson, J.B.; Lechler, P.J.; Kondrad, S.L.; Galbreath, P.F.; Salter, E.B. Influence of temporal variations in water chemistry on the Pb isotopic composition of rainbow trout (Oncorhynchus mykiss). Sci. Total Environ. 2005, 350, 204-224.

148. Soto-Jimenez, M.F.; Paez-Osuna, F.; Scelfo, G.; Hibdon, S.; Franks, R.; Aggarawl, J.; Flegal, A.R. Lead pollution in subtropical ecosystems on the SE Gulf of California Coast: A study of concentrations and isotopic composition. Mar. Environ. Res. 2008, 66, 451-458.

149. Soto-Jimenez, M.F.; Flegal, A.R. Origin of lead in the Gulf of California ecoregion using stable isotope analysis. J. Geochem. Explor. 2009, 101, 209-217.

150. S0ndergaard, J.; Asmund, G.; Johansen, P.; Elberling, B. Pb isotopes as tracers of mining-related Pb in lichens, seaweed and mussels near a former Pb-Zn mine in West Greenland. Environ. Pollut. 2010, 158, 1319-1326.

151. Potot, C.; Féraud, G.; Schärer, A.; Barats, A.; Durrieu, G.; Le Poupon, C.; Travi, Y.; Simler, R. Groundwater and river baseline quality using major, trace elements, organic carbon and Sr-Pb-O isotopes in a Mediterranean catchment: The case of the Lower Var Valley (south-eastern France). J. Hydrol. 2012, 472-473, 126-147.

152. Sangster, D.F.; Outridge, P.M.; Davis, W.J. Stable lead isotope characteristics of lead ore deposits of environmental significance. Environ. Rev. 2000, 8, 115-147.

153. Balcaen, L.; Moens, L.; Vanhaecke, F. Determination of isotope ratios of metals and metalloids by means of ICP mass spectrometry for provenanceing purposes. Spectrochim. Acta Part B At. Spectrosc. 2010, 65, 769-786.

154. Thirwall, M.F. Inter-laboratory and other errors in Pb isotope analyses investigated using Pb-207/Pb-204 double spike. Chem. Geol. 2000, 163, 299-322.

155. Yeager, K.M.; Santschi, P.H.; Herbert, B.E. Suspended sources and tributary effects in the lower reaches of a coastal plain stream as indicated by radionuclides, Loca Bayou, Texas. Environ. Geol. 2005, 47, 382-395.

156. Komarek, M.; Ettler, V.; Chrastny, V.; Milhaljevic, M. Lead isotopes in environmental sciences: A review. Environ. Int. 2008, 34, 562-572.

157. Hopper, J.R.; Ross, H.B.; Sturges, W.T.; Barrie, L.A. Regional source discrimination of atmsoperhic aerosols in Europe using the isotopic composition of lead. Tellus 1991, 43B, 45-60.

158. Gobiel, C.; Johnson, W.K.; MacDonald, R.W.; Wong, C.S. Sources and burden of lead in St. Lawrence Estuary sediments: Isotopic evidence. Environ. Sci. Technol. 1995, 29, 193-201.

159. Li, X.D.; Shen, Z.G.; Wai, O.W.H.; Li, Y.S. Chemical forms of Pb, Zn, and Cu in the sediment profiles of the Pearl River Estuary. Mar. Pollut. Bull. 2001, 42, 215-223.

160. Negrel, P.; Kloppman, W.; Garcin, M.; Giot, D. Lead isotope signatures of Holocene fluvial sediments from the Loire River valley. Appl. Geochem. 2004, 19, 957-972.

161. Xu, B.; Gu, Z.Y.; Zhang, Y.H.; Chen, Y.F.; Lu, Y.W. Sequence extraction and isotope analysis for discriminating the chemical forms and origins of Pb in sediment from Liaodong Bay, China. Arch. Environ. Con. Tox. 2009, 57, 230-238.

162. Shepherd, T.J.; Chenery, S.R.N.; Pashley, V.; Lord, R.A.; Ander, L.E.; Breward, N.; Hobbs, S.F.; Horstwood, M.; Klinck, B.A.; Worral, F. Regional lead isotope study of a polluted river catchment: River Wear, Northern England, UK. Sci. Total Environ. 2009, 407, 4882-4892.

163. British Geological Survey; Dunham, K.C. Geology of the Northern Pennine Orefield, Tyne to Stainmore; The Stationery Office: Norwich, UK, 1990.

164. Elbaz-Paulichet, F.; Hollinger, P.; Martin, J.M.; Petit, D. Stable lead isotope ratios in major French rivers and estuaries. Sci. Total Environ. 1986, 54, 61-76.

165. McGill, R.A.R.; Pearce, J.M.; Fortey, N.J.; Watt, J.; Ault, L.; Parrish, R.R. Contamination source apportionment by PIMMS lead isotope analysis and SEM-image analysis. Environ. Geochem. Health 2003, 25, 25-32.

166. Kurkjan, R.; Dunlap, C.; Flegal, A.R. Long-range downstream effects of urban runoff and acid mine drainage in the Debed River, Armenia: Insights from lead isotope modeling. Appl. Geochem. 2004, 19, 1567-1580.

167. Ettler, V.; Mihaljevic, M.; Sebeck, O.; Molek, M.; Grygar, T.; Zeman, J. Geochemical and Pb isotopic evidence for sources and dispersal of metal contamination in stream sediments from the mining and smelting district of Pribram, Czech Republic. Environ. Pollut. 2006, 142, 409-417.

168. Bird, G.; Brewer, P.A.; Macklin, M.G.; Nikolova, M.; Kotsev, T.; Mollov, M.; Swain, C. Pb isotope evidence from contaminant-metal dispersal in an international river system: The lower Danube catchment, Eastern Europe. Appl. Geochem. 2010, 25, 1070-1084.

169. Bird, G.; Brewer, P.A.; Macklin, M.G.; Nikolova, M.; Kotsev, T.; Mollov, M.; Swain, C. Quantifying sediment-associated dispersal using Pb isotope: Application of binary and multivariate mixing models in the Maritsa catchment, Bulgaria. Environ. Pollut. 2010, 158, 2158-2169.

170. Townsend, A.T.; Seen, A.J. Historical lead isotope record of a sediment core from the Derwent River (Tasmania, Australia): A multiple source environment. Sci. Total Environ. 2012, 424, 153-161.

171. Erel, Y.; Veron, A.; Halicz, L. Tracing the transport of anthropogenic lead in the atmosphere and in soils using isotope ratios. Geochim. Cosmochim. Acta 1997, 61, 4495-4505.

172. Miller, J.R.; Lechler, P.J.; Mackin, G.; Germanoski, D.; Villarroel, L.F. Evaluation of particle dispersal from mining and milling operations using lead isotopic fingerprinting techniques, Rio Pilcomayo Basin, Bolivia. Sci. Total Environ. 2007, 384, 355-373.

173. Church, S.E.; Fey, D.L.; Unruh, D.M. Trace elements and lead isotopes in streambed sediments in streams affected by historical mining. In Integrated Investigations of Environmental Effects of Historical Mining in the Animas River Watershed, San Juan County, Colorado; USGS Professional Paper 1651; Church, S.E., von Guerard, P., Finger, S.E., Eds.; U.S. Government Printing Office: Washington, DC, USA, 2007; pp. 283-335.

174. Hudson-Edwards, K.A.; Macklin, M.G.; Miller, J.R.; Lechler, P.J. Sources, distribution and storage of heavy metals in the Rio Pilcomayo, Bolivia. J. Geochem. Explor. 2001, 72, 229-250.

175. Monna, F.; Lancelot, J.; Croudace, I.W.; Cundy, A.B.; Lewis, J.T. Pb isotopic composition of airborne particulate materials from France and the Southern United Kingdom: Implications for Pb pollution sources in urban areas. Environ. Sci. Technol. 1997, 31, 2277-2286.

176. N'Guessan, Y.M.; Probst, J.L.; Bur, T.; Probst, A. Trace elements in stream bed sediments from agricultural catchments (Gascogne region, S-W France): Where do they come from? Sci. Total Environ. 2009, 407, 2939-2952.

177. Ip, C.C.M.; Li, X.D.; Zhang, G.; Wai, O.W.H.; Li, Y.S. Trace metal distribution in sediments of the Pearl River Estuary and the surrounding coastal area, South China. Environ. Pollut. 2007, 147, 311-323.

178. Batista, M.J.; de Oliveira, D.P.S.; Abreu, M.M.; Locutura, J.; Sheperd, T.; Matos, J.; Bel-Lan, A.; Martins, L. Sources, background and enrichments of lead and other elements: Lower Guadiana River. Geoderma 2013, 193-194, 265-274.

179. Nelson, B.K.; DePaolo, D.J. Application of Sm-Nd and Rb-Sr isotope systematics to studies of provenance and basin analysis. J. Sediment. Petrol. 1988, 58, 348-357.

180. Awwiller, D.N.; Marck, L.E. Diagenetic resetting of Sm-Nd isotope systematics in Wilcox Groups sandstones and shales, San Marcos Arch, south-central Texas. Gulf Coast Assoc. Geol. Soc. Trans. 1989, 39, 321-330.

181. McLennan, S.M. Rare earth elements in sedimentary rocks: Influence of provenance and sedimentary processes. Rev. Miner. Geochem. 1989, 21, 169-200.

182. Awwiller, D.N. Geochronology and mass transfer in Gulf Coast mudrocks (south-central Texas, U.S.A.): Rb-Sr, Sm-Nd and REE systematics. Cehm Geol. 1994, 116, 61-84.

183. Gao, B.; Liu, Y.; Sun, K.; Liang, X.R.; Peng, P.; Sheng, G.; Fu, J. Precise determination of cadmium and lead isotopic compositions in river sediments. Anal. Chim. Acta 2008, 612, 114-120.

184. Cloquet, C.; Carignan, J.; Libourel, G.; Sterckeman, T.; Perdrix, E. Tracing source pollution in soils using cadmium and lead isotopes. Environ. Sci. Technol. 2006, 40, 2525-2530.

185. Weiss, D.J.; Rehkamper, M.; Schoenberg, R.; McLaughlin, M.; Kirby, J.; Campbell, P.G.C.; Arnold, T.; Chapman, J.; Peel, K.; Gioia, S. Application of nontraditional stable-isotopes to the study of sources and fate of metals in the environment. Environ. Sci. Technol. 2008, 42, 655-664.

186. Shiel, A.E.; Weis, D.; Orians, K.J. Evaluation of zinc, cadmium, and lead isotope fractionation during smelting and refining. Sci. Total Environ. 2010, 408, 2357-2368.

187. Shiel, A.E.; Weis, D.; Orians, K.J. Tracing cadmium, zinc, and lead sources in bivalves from the coasts of western Canada and the USA using isotopes. Geochim. Cosmochim. Acta 2012, 76, 175-190.

188. Kimball, B.E.; Mathur, R.; Dohnalkova, A.C.; Wall, A.J.; Runkel, R.L.; Brantley, S.L. Copper isotope fractionation in acid mine drainage. Geochim. Cosmochim. Ac. 2009, 73, 1247-1263.

189. Zhu, X.K.; Onions, R.K.; Guo, Y.; Belshaw, N.S.; Rickard, D. Determination of natural cu-isotope variation by plasma source mass spectrometry: Implications for use as geochemical tracers. Chem. Geol. 2000, 763, 139-149.

190. Petit, J.C.J.; DeJong, J.; Chou, L.; Mattielli, N. Development of Cu and Zn isotope MC-ICP-MS measurements: Application to suspended particulate matter and sediments from the Scheldt Estuary. Geostand. Geoanal. Res. 2008, 32, 149-166.

191. Bergquist, B.; Blum, J.D. The odds and evens of mercury isotopes: Applications of mass-dependent and mass-independent isotope fractionation. Elements 2009, 5, 353-357.

192. Yi, R.; Feng, X.; Shi, W. Application of the stable-isotope system to the study of sources and fate of Hg in the environment: A review. Appl. Geochem. 2010, 25, 1467-1477.

193. Sonke, J.E.; Schäfer, J.; Chmeleff, J.; Audry, S.; Blanc, G.; Dupré, B. Sedimentary mercury stable isotope records of atmospheric and riverine pollution from two major European heavy metal refineries. Chem. Geol. 2010, 279, 90-100.

194. Rouxel, O.; Ludden, J.; Fouquet, Y. Antimony isotope variations in natural systems and implications for their use as geochemical tracers. Chem. Geol. 2003, 200, 25-40.

195. Wilson, L.R. Determination of Trace Element Provenance, Rio Loa Basin, Northern Chile. Master's Thesis, Department of Chemistry and Physics, Western Carolina University, Cullowee, NC, USA, July 2011.

196. Sivry, Y.; Riotte, J.; Sonke, J.E.; Audray, S.; Schäfer, J.; Viers, J.; Blanc, G.; Freydier, R.; Dupré, B. Zn isotopes as tracers of anthropogenic pollution from Zn-ore smelters, The Riou Mort-Lot River System. Chem. Geol. 2008, 255, 295-304.

197. Bentahila, Y.; Othman, B.; Luck, J.M. Stontium, lead, and zinc isotopes in marine cores as tracers of sedimentary provenance: A case study around Taiwan orogen. Chem. Geol. 2008, 248, 62-82.

198. Wombacher, F.; Rahkämper, M.; Mezger, K. Determination of the massdependence of cadmium isotope fractionation during evaporation. Geochim. Cosmochim. Acta 2004, 68, 2349-2357.

199. Rahkämper, M.; Wombacher, F.; Horner, T.J.; Xue, Z. Natural and anthropogenic cd isotope variations. In Handbook of Environmental Isotope Geochemistry, Advances in Isotope Geochemistry; Baskarn, M., Ed.; Springer-Verlag: Berlin, Germany, 2011; pp. 125-154.

200. Bullen, T.D. Stable isotope of transition and post-transition metals as tracers in environmental studies. In Handbook of Environmental Isotope Geochemistry; Advances in Isotope Geochemistry; Baskarn, M., Ed.; Springer-Verlag: Berlin, Germany, 2011; pp. 177-203.

201. Geagea, M.L.; Stille, P.; Gauthier-Lafaye, F.; Millet, M. Tracing of industrial aerosol sources in an urban environment using Pb, Sr, and Nd isotopes. Environ. Sci. Tech. 2008, 42, 692-698.

202. Weiss, D.J.; Rausch, N.;Mason, T.F.D.; Coles, B.J.; Wilkinson, J.J.; Ukonmaanaho, L.; Nieminen, T. Geochim. Atmospheric deposition and isotope biogeochemistry of zinc in ombrotrophic peat. Geochim. Cosmochim. Acta 2007, 71, 942-960.

203. Cloquet, C.; Carignan, J.; Lehmann, M.R.; Vanhaecke, F. Variation in the isotopic compositon of zinc in the natural environment and the use of zinc isotopes in biogeosciences: A review. Anal. Bioanal. Chem. 2008, 390, 451-463.

204. Chen, J.; Gaillardet, J.; Louvat, P. Zinc isotopes in the Seine River waters, France: A probe of anthropogenic contamination. Environ. Sci. Technol. 2008, 42, 6494-6501.

205. Aranda, S.; Borrok, D.M.; Wanty, R.B.; Balistrien, L.S. Zinc isotope investigation of surface and pore waters in a mountain watershed impacted by acid rock drainage. Sci. Total Environ. 2012, 420, 202-213.

206. Wanty, R.B.; Podda, F.; de Giudici, G.; Cidu, R.; Lattanzi, P. Zinc isotope and transition-element dynamics accompanying hydrozincite biomineralization in the Rio Naracauli, Sardinia, Italy. Chem. Geol. 2013, 337-338, 1-10.

207. Lambelet, M.; Rehkamper, M.; van de Flierdt, T.; Xue, Z.; Kreissig, K.; Coles, B.; Porcelli, D.; Andersson, P. Isotopic analysis of Cd in the mixing zone of Siberian rivers with the Arctic Ocean-New constraints on marine Cd cycling and the isotope composition of riverine Cd. Earth Plant. Sci. Lett. 2013, 361, 64-73.

208. Cloquet, C.; Rouxel, O.; Carignan, J.; Libourel, G. Natural cadmium isotopic variations in eight geological reference materials (NIST SRM 2711, BCR 176, GSS-1, GXR-1, GXR-2, GSD-12, Nod-P-1, Nod-A-1) and anthropogenic samples, measured by MC-ICP-MS. Geostand. Geoanal. Res. 2005, 29, 95-106.

209. Ripperger, S.; Rehkamper, M. Precis determination of cadmium isotope fractionation in seawater by double spike MC-ICP-MS. Geochem. Cosmochim. Acta 2007, 71, 631-642.

210. Bullen, T.D.; Walczyk, T. Environmental and Biomedical application of natural metal stable isotope variations. Elements 2009, 5, 381-385.

211. Mathur, R.; Ruiz, J.; Titley, S.; Liermann, L.; Buss, H.; Brantley, S. Cu isotopic fractionation in the supergene environment with and without bacteria. Geochim. Cosmochim. Acta 2005, 69, 5233-5246.

212. Zink, S.; Schoenberg, R.; Staubwasser, M. Isotope fractionation and reaction kinetics between Cr(III) and Cr(VI) in aqueous media. Geochim. Cosmochim. Acta 2010, 5729-5745.

213. Kitchen, J.W.; Johnson, T.M.; Bullen, T.D.; Zhu, J.; Raddatz, A. Chromium isotope fractionation factors for reduction of Cr(VI) by aqueous Fe(II) and organic molecules. Geochim. Cosmochim. Acta 2012, 89, 190-201.

214. Borrok, DM.; Wanty, R.B.; Ridley, W.I.; Lamothe, P.J.; Kimball, B.A.; Verplanck, P.L.; Runkel, R.L. Application of iron and zinc isotopes to track the sources and mechanisms of metal loading in a mountain watershed. Appl. Geochem. 2009, 24, 1270-1277.

215. Fernandez, A.; Borrok, D.M. Fractionation of Cu, Fe, and Zn isotopes during the oxidative weathering of sulphide-rich rocks. Chem. Geol. 2009, 264, 1-12.

216. Hoefs, J. Stable Isotope Geochemistry; Springer Verlag: New York, NY, USA, 2009.

217. Bullen, T.D. Chromium Stable Isotopes as a New Tool for Forensic Hydrology at Sites Contaminated with Anthropogenic Chromium. In Water-Rock Interactions, Proceedings of 12th International Symposium on Water-Rock Interactions (WRI-12), Kunming, China, 31 July-5 August 2007; Bullen, T.D., Wang, Y., Eds.; Taylor & Francis: London, UK, 2007; pp 699-702.

218. Raddatz, A.L.; Johnson, T.M.; McLing, T.L. Cr stable isotopes in Snake River plan aquifer groundwater: Evidence for natural reduction of dissolved Cr(VI). Environ. Sci. Technol. 2011, 45, 502-507.

219. Wanner, C.; Eggenberger, U.; Kurz, D.; Zink, S.; Mader, U. A chromate-contaminated site in southern Switzerland—Part I: Site characterization and the use of Cr isotopes to delineate fate and transport. App. Geochem. 2012, 27, 644-654.

220. Sherman, L.S.; Blum, J.D. Mercury stable isotopes in sediments and largemouth bass from Florida lakes, USA. Sci. Total Environ. 2013, 448, 163-175.

221. Hintelmann, H.; Lu, S.Y. High precision isotope ratio measurements of mercury isotopes in cinnabar ores using multi-collector inductively coupled plasma mass spectrometry. Analyst 2003, 128, 636-639.

222. Smith, C.N.; Kesler, S.E.; Kaue, B.; Blum, J.D. Mercury isotope fractionation in fossil hydrothermal systems. Geology 2005, 33, 825-828.

223. Smith, C.N.; Kesler, S.E.; Blum, J.D.; Rytuba, J.J. Isotope geochemistry of mercury isn source rocks, mineral depsoits and spring deposits of the California Coast Ranges, USA. Earth Planet. Sci. Lett. 2008, 269, 399-407.

224. Lauretta, D.S.; Klaue, B.; Blum, J.D.; Buseck, P.R. Mercury abundances and isotopic compositions in the Murchison (CM) and Allende (CV) carbonaceous chondrites. Geochemica Cosmochim. Acta 2001, 65, 2807-2818.

225. Sonke, J.E.; Zambardi, T.; Toutain, J.P. Indirect gold trap-MC-ICP-MS coupling for Hg stable isotope analysis using a syringe injection interface. J. Anal. Atomic Spectrom. 2008, 23, 569-573.

226. Zambardi, T.; Sonke, J.E.; Toutain, J.P.; Sortino, F.; Shinohara, H. Mercury emissions and stable isotopic compositions at Vulcano Island (Italy). Earth Planet. Sci. Lett. 2009, 227, 236-243.

227. Jackson, T.A.; Muir, D.C.G.; Vincent, W.F. Historical variations in the stable isotope composition of mercury in Arctic lake sediments. Environ. Sci. Technol. 2004, 38, 2813-2121.

228. Foucher, D.; Hintelmann, H. High-precision measurements of mercury isotope ratios in sediments using cold-vapor generation multi-collector inductively coupled plasma mass spectrometry. Anal. Bioanal. Chem. 2006, 384, 1470-1478.

229. Biswas, A.; Blum, J.D.; Bergquist, B.A.; Keeler, G.J.; Xie, Z.Q. Natural mecury isotope variation in coal deposits and organic soils. Environ. Sci. Technol. 2008, 42, 8303-8309.

230. Ghosh, S.; Xu, Y.F.; Humayun, M.; Odom, L. Mass-independent fractionation of mercury isotopes in the environment. Geochem. Geophy. Geosyst. 2008, 9, doi:10.1029/2007GC001827.

231. Jackson, T.A.; Whittle, D.M.; Evans, M.S.; Muir, D.C.G. Evidence for mass-independent and mass-dependent fractionation of the stable isotopes of mercury by natural processes in aquatic ecosystems. Appl. Geochem. 2008, 23, 547-571.

232. Foucher, D.; Ogrinc, N.; Hintelmann, H. Tracing mercury contamination from the Idraja mining region (Slovenia) to the Gulf of Trieste using Hg isotope ratio measurements. Environ. Sci. Technol. 2009, 43, 33-49.

233. Gehrke, G.E.; Blum, J.D.; Meyers, P.A. The geochemical behavior and isotopic composition of Hg in a mid-Pleistocene western Mediterranean sapropel. Geochim. Cosmochim. Acta 2009, 73, 1651-1665.

234. Carignan, J.; Estrade, N.; Sonke, J.E.; Donard, O.F.X. Odd isotope deficits in atmospheric Hg measured in lichens. Environ. Sci. Technol. 2009, 43, 5660-5664.

235. Blum, J.D.; Bergquist, B.A. Reporting of variations in natural isotopic compositions of mercury. Anal. Bioanal. Chem. 2007, 388, 353-359.

236. Feng, X.B.; Foucher, D.; Hintelmann, H.; Yan, H.Y.; He, T.R.; Qiu, G.L. Tracing mercury contamination sources in sediments using mercury isotope compositions. Environ. Sci. Technol. 2010, 44, 3363-3368.

237. Foucher, D.; Hintelmann, H.; Al, T.A.; MacQuarrie, K.T. Mercury isotope fractionation in waters and sediments of the Murray Brook mine watershed (New Brunswick, Canada): Tracing mercury contamination and transformations. Chem. Geol. 2013, 336, 87-95.

238. Yin, R.; Feng, X.; Wang, J.; Li, P.; Liu, J.; Zhang, Y.; Chen, J.; Zheng, L.; Hu, T. Mercury speciation and mercury isotope fractionation during ore roasting process and their implication to source identification of downstream sediment in the Wanshan mercury mining area, SW China. Chem. Geol. 2013, 336, 72-79.

© 2013 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article

distributed under the terms and conditions of the Creative Commons Attribution license


Copyright of Minerals (2075-163X) is the property of MDPI Publishing and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.