Abstract
This study presents a new 3600-year record of past metal contamination from a bog located close to the Leadhills and Wanlockhead orefield of southwest Scotland. A peat core, collected from Toddle Moss, was radiocarbon (14C) dated and analysed for trace metal concentrations (by EMMA) and lead isotopes (by ICP-MS) to reconstruct the atmospheric deposition history of trace metal contamination, in particular, lead. The results show good agreement with documented historical and archaeological records of mining and metallurgy in the region: the peak in metal mining during the 18th century, the decline of lead mining during the Anglo-Scottish war and lead smelting during the early medieval period. There may also have been earlier workings during the Late Bronze and Iron Ages indicated by slight increases in lead concentrations, the Pb/Ti ratio and a shift in 206Pb/207Pb ratios, which compare favourably to the signatures of a galena ore from Leadhills and Wanlockhead. In contrast to other records across Europe, no sizeable lead enrichment was recorded during the Roman Iron Age, suggesting that the orefield was not a significant part of the Roman lead extraction industry in Britain. These findings add to the various strands of archaeological evidence that hint at an early lead extraction and metallurgical industry based in southern Scotland. The results also provide further evidence for specific regional variations in the evolution of mining and metallurgy and an associated contamination signal during prehistoric and Roman times across Europe.
Introduction
Since its first use approximately seven millennia ago, lead has played an important role in human history, including aspects of art, medicine and technology (Hong et al., 1994; Nriagu, 1998). Lead became particularly important in the 5th millennium
Lead in Bronze Age artefacts confirms that it was being used in Britain prior to 1500
Archaeological evidence for lead extraction, objects and its wider use in later prehistory is also scant. Lead was a rare addition to northern British alloys during the Iron Age (Dungworth, 1996), and only a few examples of lead use in ornaments, or as solder to secure iron objects, exist in Scotland (Toolis, 2007). Roman lead extraction and smelting is more abundant, with known mines and/or smelting sites occurring at Charterhouse in the Mendips, SW Britain, Alderley Edge in Cheshire and sites in Wales (Timberlake et al., in press), but so far, there are no known sites in Scotland. One region that has been not fully investigated is the Leadhills and Wanlockhead orefield in SW Scotland. A study is therefore warranted, especially given the recent discovery of a stone hammer in Wanlockhead, which is indicative of early mining (Pickin, 2008).
The aim of this study was to reconstruct the history of exploitation of insular ore sources in the Leadhills/Wanlockhead orefield from prehistory to the present. To do so, we present an atmospheric metal contamination history from an ombrotrophic mire, Toddle Moss, for the last 3600 years using total concentrations, Pb/Ti ratios and lead isotope ratios (206Pb/207Pb) contained in the peat.
Materials and methods
Location, sampling and sub-sampling strategy
The orefield of Leadhills and Wanlockhead, on the border of Dumfriesshire and Lanarkshire, SW Scotland (Figure 1), is rich in metalliferous deposits as a result of two phases of mineralization in Ordovician sediments: a quartz vein mineral phase of Carboniferous age and a possibly metallic one of lead and zinc, both of the Carboniferous era. This was followed by a later phase of secondary enrichment (Pattrick and Polkya, 1993).

(a) Location of study area in Britain, (b) location of Leadhills, Wanlockhead and Toddle Moss and (c) evidence of early lead working after Pickin (2010). The place name ‘bail’ is considered to indicate a lead smelting site.
Toddle Moss is located approximately 4 km northwest of the village of Leadhills and 0.5 km east of Elvanfoot in the Elvan Water river valley (Figure 1). This area has a rich history of mining. Alluvial sediments have been worked in this valley for gold, and lodes rich in lead and copper have also been exploited, particularly in the
A 7.5-m-deep Sphagnum–Eriophorum-rich peat core was taken from Toddle Moss using a Russian corer with 30 cm × 10 cm chamber in 2004. The samples were placed in plastic guttering, wrapped in polythene and placed in cold storage. The top 3.5 m was analysed in this study.
Chronology
Radiocarbon dates were determined at Beta Analytic Ltd (Miami) and the Poznań Radiocarbon Laboratory using conventional and accelerator mass spectrometry (AMS) methods, respectively. The total carbon of one sample of fresh peat (Beta-15142) was dated using conventional radiometric dating. For AMS, a 1-cm-thick slice of peat (Poz-19215) and Sphagnum macrofossils (Poz-56748) were selected for dating (cf. Nilsson et al., 2001; Piotrowska et al., 2011).
To provide a highly resolved chronology for the last 100–150 years, the unsupported 210Pbun activity within samples towards the peat surface was ascertained by subtraction of the supported component (measured as 214Pb at 295.22 and 351.93 keV) from the total 210Pb activity measured at 46.54 keV (Wallbrink et al., 2002). 210Pb and 214Pb activities were measured using EG&G ORTEC hyper-pure Germanium detectors in a well configuration (11 mm diameter, 40 mm depth) housed at Coventry University. The method for calculating the depth–age relationship follows procedures described by Appleby and Oldfield (1978), Appleby (2001) and Walling et al. (2002). The CRS dating model was used to calculate ages as accumulation rates varied down core (Appleby, 2001; Appleby et al., 1988). The CLAM software package (Blaauw, 2010) was used to create an age–depth model, combining the 14C and 210Pb ages, to infer approximate ages for all levels.
Geochemical analyses
The core was cut into contiguous 1-cm-thick slices, oven dried at 40°C and homogenized to improve the efficiency of the chemical digest and to provide better representation of the total metal concentration within the samples. Calcium and magnesium were determined by inductively coupled plasma–optical emission spectrometry (ICP-OES). An estimation of the efficiency of the digestion method and of the accuracy of the analytical measurements was obtained through the use of replicate sub-samples, spiked blanks and certified reference materials (Ebdon et al., 1998; Fifield and Kealey, 2000). Spiked samples of known concentration (10 mg/L) were used to test the efficiency of the acid microwave digestion. Two certified reference materials were also used: Sphagnum energy peat (NJV 94-2) and Carex energy peat (NJV 94-1). The reference materials were imported from the Swedish University of Agricultural Sciences, Department of Agricultural Research for Northern Sweden, Laboratory for Chemistry and Biomass. Standards of known metal concentrations were used to calibrate the instrument and to ensure it was performing at its optimum efficiency (Holler et al., 1996). Total metal concentrations in each sample are expressed in units of microgram per gram.
Recovery of calcium and magnesium from the CRMs for Sphagnum and Carex was between 101% and 130%. Three spiked samples yielded recovery of between 101% and 114%. These results suggest that metal recovery using the microwave digestion was very efficient and that the ICP-OES provided reliable data.
The elemental composition of 89 dried, milled and homogenized samples between 0 and 300 cm depth were obtained by EMMA-XRF analyses (Cheburkin and Shotyk, 1996; Weiss and Shotyk, 1998) including concentrations of major and trace lithogenic elements (silicon, aluminium, titanium, gallium, yttrium and zirconium) and trace metals and metalloids (lead, chromium and arsenic). The instruments are hosted at the RIAIDT (Infrastructure Network for the Support of Research and Technological Development) facility of the University of Santiago de Compostela, Spain. Standard reference materials were used for the calibration of the instruments. Quantification limits were 0.001% for Ti, 0.01% for Al, 0.05% for Si, 0.5 µg/g for Pb and 1 µg/g for other trace elements. Replicate measurements were taken for one of every five samples in order to account for reproducibility; all replicates were within 5% agreement.
A total of 28 sub-samples of peat from the same core were selected for lead isotope analysis at the School of Geosciences, University of Edinburgh. Sub-samples (~0.25 g) were air-dried, then washed at 450°C for 4 h and finally digested using a modified US EPA Method 3052 Protocol microwave-assisted HF-HNO3 digestion method (Yafa and Farmer, 2006; Yafa et al., 2004). Digests were evaporated to 1 mL on a hotplate and then made up to 25 mL with 2% (v/v) HNO3. All reagents used in sample preparation were of the highest analytical quality available, that is, Aristar nitric acid (69%) and hydrofluoric acid (48%) and high purity water (18.2 MΩ cm) from a Milli-Q water system (Millipore, Watford, UK). Lead isotopic ratios were determined in the prepared 2% v/v HNO3 solutions using a PlasmaQuad (PQ) 3 ICP-MS instrument (Thermo Electron, Winsford, UK), equipped with a Meinhard nebulizer, nickel sampler and skimmer cones, Gilson autosampler and a Gilson Minipuls 3 peristaltic pump (Anachem, Luton, UK). A solution of the National Institute of Standards and Technology (NIST) common lead isotopic reference standard SRM 981 (206Pb/207Pb = 1.093, 208Pb/206Pb = 2.168, 208Pb/207Pb = 2.370) was used for calibration and mass bias correction (Farmer et al., 2000). Analytical precision on these ratios was typically <±0.2%.
To ensure the quality of analytical procedures and data, an ombrotrophic peat reference material (NIMT/UOE/FM001) (Yafa et al., 2004) was analysed along with the samples. The mean values (n = 5) of 1.177 ± 0.001, 2.093 ± 0.002 and 2.463 ± 0.004 determined for the isotope ratios 206Pb/207Pb, 208Pb/206Pb and 208Pb/207Pb, respectively, in the reference material were in good agreement with corresponding ‘information only’ values of 1.176 ± 0.001, 2.092 ± 0.002 and 2.461 ± 0.003 reported in Yafa et al. (2004).
Statistics
Following the procedure described by Martínez Cortizas et al. (2013) and Hermanns and Biester (2013), we used factor analysis by principal component analysis (PCA) to identify sources and processes related to the distribution of the measured elements using the SPSS 20 software package. PCA of compositional data is usually undertaken on transformed variables (Baxter, 1995), particularly when the values cover several orders of magnitude and there are outliers (Baxter, 1999). Transformation also avoids any scaling effects (Eriksson et al., 1999). Thus, the PCA was done on log-transformed and standardized (z-scores) data, using varimax rotation to maximize the variance of the elements in the principal components (Eriksson et al., 1999). Similarly distributed elements will load on to the same principal component and are most likely to be controlled by the same environmental factor(s). Hence, interpretation of the signals with regard to the underlying cause or causes of variation of a group of elements should be more evident.
Results
Age–depth modelling
The results are shown in Table 1 with 2σ calibrated age ranges (in calibrated years
Radiocarbon dates from Toddle Moss.

An age–depth model for Toddle Moss using Clam (after Blaauw, 2010).
Geochemistry
The calcium/magnesium ratio (Figure 3) is consistently below 1 and much lower than the measured value for rainwater (approximately 1.9) at Raeburn Flow, which is located approximately 74 km to the southeast of Toddle Moss (Küttner et al., 2014). These low values are sufficient to infer ombrotrophic conditions for the bog (cf. Shotyk, 1996). Notwithstanding a series of short-lived peaks, titanium concentrations remain relatively constant from the base of the profile to 101 cm. They then rise gradually with a sustained increase from 43 cm to the surface of the bog. Lead concentrations are generally low from the base of the core to 35 cm depth: thereafter, they increase dramatically to peak at 11 cm before decreasing to much lower concentrations at the bog surface. The Pb/Ti ratio follows a similar pattern. The 206Pb/207Pb profile shows more radiogenic ratios from 342 to 308 cm (SI Table 1, available online). A shift to less radiogenic values occurs between 308 and 222 cm (Figure 3). Thereafter, the ratios gradually rise from 222 to 24 cm before they become less radiogenic towards the bog surface.

Calcium/magnesium ratio (by ICP-OES), titanium concentrations, lead concentrations (by EMMA; dashed line on a log scale), Pb/Ti ratio and 206Pb/207Pb ratios (by ICP-MS) from Toddle Moss.
Concentrations of arsenic, chromium, gallium, yttrium, zinc and zirconium, determined using EMMA, are shown in Figure 4. Arsenic, gallium and yttrium have similar trends to lead: low concentrations from the base of the core until 35 cm, a sharp rise to peak in the top 10 cm, followed by a decline to much lower concentrations. Zinc and chromium are characterized by low but highly fluctuating concentrations from the base of the core until approximately 80 cm. They also peak in the uppermost 10 cm with other elements. Chromium concentrations also decline close to the peat surface, but zinc remains relatively high. Zirconium concentrations are more erratic with a series of peaks throughout the profile (e.g. 265, 237, 179, 135, 113, 63, 33 and 17 cm) but show a gradual increase in the upper metre of the profile.

Concentrations of lead, arsenic, gallium, yttrium, zinc and chromium as determined by EMMA.
PCA
Two components explain 78% of the total variance (Table 2). The first component (Cp1, 44% of the total variance) is characterized by large-to-moderate positive loadings of metals typically associated with mining/metallurgy (namely, arsenic, zinc, lead and chromium; Table 2). Gallium and yttrium, which are usually considered to be lithogenic elements, show large loadings in Cp1 (Table 2) and thus can also be associated with atmospheric metal pollution, probably derived from dust emissions during mining. This is also supported by the extremely high metal concentrations in the upper section of the peat, which are only comparable with those found close to pollution sources (e.g. in the Harz mountains, Germany, where Pb concentrations exceed 1000 µg/g during medieval times (Kempter and Frenzel, 2000) and equivalent to Pb concentrations determined within several kilometres of a lead smelter (Mihaljevič et al., 2006)). Cp1 scores show a typical record of atmospheric metal pollution, with a large peak in the upper 30 cm of the core and a sharp decrease in the upper 8 cm (Figure 5). It also shows two minor increases in scores during the early medieval period: between 120 and 140 cm (5th–7th centuries
Loadings of the variables in the components extracted by PCA on the chemical composition of the peat.
Cp1 and Cp2: components extracted using the whole data set; Cp1-PI to Cp3-PI: components extracted using data for pre-industrial peat sections; Eigv: eigenvalues; Var: percentage of total variance; Com: communality (proportion of variance of each element explained by the two principal components); PCA: principal component analysis.

Records of scores of the extracted principal components.
The second component (Cp2, 33.7% of the total variance) is characterized by large positive loadings of the lithogenic elements (Ti, Si, Al and Zr; Table 2). This chemical association reflects the mineral content of the peat because of deposition of dust, probably derived from soil erosion. The record of scores shows a ‘see-saw’ pattern, with eight peaks in dust deposition (Figure 5): at 265 cm (c. 870 cal.
Due to the expected large effect on the PCA of the metal concentrations of the peat sections with ages younger than
Moreover, in the second analysis, the metal signal is divided into two components: Cp2-PI with zinc, arsenic and chromium which has a record of scores similar to that of Cp2 (r = 0.68) and Cp3-PI with gallium and lead (Table 2). Thus, the metal signature of the peat for pre-industrial times seems to indicate that there were differences in the accumulation of the metals in Toddle Moss. The record of Cp2-PI scores suggests that the history of zinc, arsenic and chromium enrichment is quite similar during late prehistory and early Middle Ages. However, the elevated scores of Cp2-PI between 250 and 236 cm, corresponding to the period c. 700–550 cal.
Interpretation and discussion
Lead is essentially immobile once it becomes incorporated into ombrotrophic peat (MacKenzie et al., 1997; Shotyk et al., 1997), and there is a plethora of studies that have demonstrated that the pattern of lead is faithfully preserved in peat bogs which can be reliably matched with other archaeological and historical documentary records (e.g. Mighall et al., 2002b; Shotyk et al., 1997). Lead concentrations, lead–titanium ratios and isotopic ratios are now regularly used to identify evidence of anthropogenic forcing on the lead biogeochemical cycle. Lead–titanium ratios are used to identify non-silicate sources of lead (Görres and Frenzel, 1997; Shotyk, 1996), whereas isotopic ratios are also considered to reflect accurately anthropogenic lead emissions especially when the isotopic signature of potential sources is well known (Martínez Cortizas et al., 2002). The record of lead derived from Toddle Moss presented here should therefore provide a reliable chronological record for past lead deposition onto the bog surface. Notwithstanding the numerous complicating variables that can influence the dispersion of gaseous and particulate pollution from source, bogs located close to industrial sites should provide robust records of emissions from these sites as pollutants are deposited onto the mire surface (Mighall et al., 2002a, 2002b).
Slightly elevated lead concentrations (Figure 3) might well be attributed to Middle–Late Bronze Age metallurgical activities: centring on 272 and 247 cm (c. 940 to 670 cal.
For purposes of comparison, the lead isotope ratios of the Toddle Moss peat samples are plotted in Figure 6 besides selected lead ores from other locations in the British Isles and the major Spanish mines of Rio Tinto and Murcia (Rohl, 1996; Stos-Gale et al., 1995). Lead isotopic signatures from Flanders Moss and Lindow Moss (cf. Cloy et al., 2005; Le Roux et al., 2004) fall within the cluster of British ores from the Mendips, Alderley Edge, NE Wales and the mines at Leadhills and Wanlockhead. Because of the overlap of the isotopic ratios, it is not possible to attribute the origin of this lead to a particular British ore source (cf. Cloy et al., 2005; Le Roux et al., 2004). Nevertheless, the results clearly show that the lead is likely to be of British origin, as the isotopic values for the British ores are clearly separable from those of the heavily exploited Spanish sources (Shotyk et al., 1998).

Plot of 208Pb/206Pb versus 206Pb/207Pb ratios from samples from the Toddle Moss core and galena from various British and Spanish ores.
The peat samples from Toddle Moss of the section between 296 and 224 cm clearly fall outside all of the clusters shown in Figure 6. An analysis of galena samples from Wanlockhead has established the isotopic signature for lead ore at this location. Cloy et al. (2005) reported a 206Pb/207Pb ratio of 1.172 ± 0.003, which is in close agreement with a value of 1.170 ± 0.003 for the Leadhills and Wanlockhead lead ore reported by Sugden et al. (1993). Rohl (1996) calculated the mean 206Pb/207Pb ratio from five ore samples from Wanlockhead plus six from Leadhills as 1.171 ± 0.001. All these values plot within the ‘British ore’ cluster (Figure 6). This suggests that the source of lead determined within the Toddle Moss peat samples between 296 and 224 cm does not originate from the main galena bearing lodes from the Leadhills or Wanlockhead orefield. However, one sample of galena from Wanlockhead has a 206Pb/207Pb ratio of 1.142, and it is plotted in Figure 6 as the ‘Wanlockhead outlier’. This ratio is in much better agreement with the Toddle Moss samples between 296 and 224 cm. If early miners and metallurgists did exploit lead at Wanlockhead, then they appear to have targeted lodes bearing a similar lead isotopic signature.
A more extensive analysis of the isotopic signatures of the mineralized zones in the orefield could resolve the apparent signature discrepancies between the peat samples, the Wanlockhead ‘outlier’ and the other Leadhills ore samples. Whether the two distinct phases of lead formation at Leadhills resulted in the different isotopic signatures is unknown. The natural lead isotope composition of rocks depends upon the age of the lithogenic system, the U/Pb and Th/Pb ratios of the system and mixing during remobilization and metamorphism (Keinonen, 1992). During the formation of lead ore deposits, lead is separated from the parent uranium and thorium isotopes, with the lead isotopic composition of hydrothermal fluids being ‘frozen’ into lead-bearing minerals (Church et al., 1993). Thus, the isotopic composition of a given ore deposit is a function of four parameters: (1) the decay rate of parent isotopes, (2) the initial ratio of the abundance of the parent to the abundance of Pb (238Pb/204Pb, 232Th/204Pb) in the source reservoir (e.g. mantle or continental crust), (3) the initial isotopic composition of reservoir Pb and (4) the duration of reservoir evolution prior to separation of Pb by geological processes (Sangster et al., 2000).
A second phase of increasing lead enrichment and relatively higher Pb/Ti ratios at Toddle Moss occurs during the late Iron Age c. cal. 365

Comparison of lead records from selected peat bogs in Scotland. (a) Leadhills lead concentrations and Pb/Ti ratio; dashed line represents exaggerated values down the profile (scale at base of the graph). (b) Raeburn Flow lead concentrations and lead enrichment factor. Enrichment factor is calculated using the equation proposed by Shotyk (1996). Titanium was used as the reference lithogenic element (from, and full details in, Küttner et al., 2014). Dashed lines represent exaggerated values down the profiles (scale at base of the graph). (c) Flanders Moss lead concentrations (from, and full details in, Cloy et al., 2005).
Shifts in the 206Pb/207Pb ratios and/or increased lead have been regularly recorded in peat profiles dated to the Late Iron and Roman times (De Vleeschouwer et al., 2010; Martínez Cortizas et al., 1997; Renberg et al., 2001). This includes sites in the British Isles: northwest and southwest England (Le Roux et al., 2004; Meharg et al., 2012), central Wales (Mighall et al., 2002b, 2009), at Flanders Moss in central Scotland (Cloy et al., 2005, 2008) and Raeburn Flow in southern Scotland (Küttner et al., 2014). The results suggest that British ores were exploited at least two centuries before the Roman occupation (Cloy et al., 2005) and that Roman exploitation always followed an earlier indigenous (British) lead extraction industry. The equivalent time frame at Toddle Moss is contained between approximately 182 and 148 cm (
The Roman occupation of Scotland was short and intermittent, spanning approximately 150 years between the late 1st century and early 3rd century
Martínez Cortizas et al. (2013) note that records from mires, lakes and lagoons do not always show metal enrichment during the Late Iron Age and Roman period. For example, some studies, such as those from southern France (Labonne et al., 1998), the Eifel area and Ireland (Schettler and Romer, 2006) do not show any evidence of contamination during the Iron Age while mires from Bavaria (Küster and Rehfuess, 1997), central and southeastern France (Baron et al., 2005; Monna et al., 2004) also show no metal enrichment during the Roman period. The Toddle Moss record provides additional evidence for specific regional variations in the evolution of mining and metallurgy during prehistoric and Roman times across Europe (cf. Martínez Cortizas et al., 2013).
The results suggest that phases of enhanced chromium accumulation occurred between 191 and 161 cm (c. 60 cal.
Despite the low concentrations, the lead record between
After
Lead, arsenic, chromium and zinc concentrations all increase above 38 cm (c.
Increased concentrations of zinc, arsenic and, to a lesser extent, gallium are most likely to be associated with lead mining, coal combustion (Oremland and Stolz, 2003; Rothwell et al., 2009; Shotyk et al., 1996) and/or plant uptake (Zaccone et al., 2008). Coal was used as fuel for lead smelting from
The marked reduction of lead, arsenic, chromium and gallium concentrations in the top 10 cm reflects the demise of the Leadhills and Wanlockhead mines in the 1930s and the subsequent phasing out of leaded gasoline. Less radiogenic 206Pb/207Pb ratios once again fall close to the isotopic ratio of the Wanlockhead outlier (Figure 6).
These values reflect a change in the source of the lead deposited onto the bog, with the final closure of the mines and the loss of the main Leadhills isotopic signature, the increasing influence of imported Australian lead (206Pb/207Pb ratio = 1.04) and other alkyl lead additives in petrol which are also phased out in the recent past. These sources would dilute any remnant lead deposition from the Leadhills galena ores. These trends are commonly recorded in bogs across the British Isles (e.g. Cloy et al., 2008; Farmer et al., 1997; Le Roux et al., 2004; Mighall et al., 2002b, 2004, 2009; West et al., 1997). The acrotelm is also affected by ongoing peat forming processes such as decomposition, plant uptake/recycling, element mobility and fluctuations at the acrotelm/catotelm transition (e.g. Martínez Cortizas et al., 2007). These processes may also have played a role in influencing the distribution of elements through the peat such as zinc (Biester et al., 2012; Espi et al., 1997; Jones, 1987; Shotyk, 1988).
Conclusion
The Toddle Moss record provides evidence of three major phases of atmospheric metal contamination which accord well with historical and archaeological records of mining and metallurgy from medieval times to present: 5th–7th centuries
Patterns of lead might reflect earlier activity during the Bronze and Iron Ages, but further analysis will be required to confirm whether the Leadhills/Wanlockhead was an early source of metals as part of an insular metal mining industry.
Footnotes
Acknowledgements
Thanks to Nathan Pittam, Claire Devonport and Alice Humphries for help with fieldwork. We thank Joanna Cloy and John Farmer for provision of lead isotope data and for helpful comments. Jason Jordan prepared some of the material for analysis. Thanks to David Cranstone for information about lead sources. We thank Alison Sandison and Jenny Johnson at the University of Aberdeen for cartographical support. Peter Gayford (Hopetoun Estates) and James C. E. Hodge kindly gave permission to take a core from Toddle Moss. We would like to thank the two anonymous reviewers for their considered and useful comments.
Funding
Coventry University and Society of Antiquaries of Scotland kindly provided funding for this research.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
