Modeling of ore-forming and geoenvironmental systems: Roles of fluid flow and chemical reaction processes

Ore-forming and geoenviromental systems commonly involve coupled fluid flowand chemical reaction processes. The advanced numerical methods and computational modeling have become indispensable tools for simulating such processes in recent years. This enables many hitherto unsolvable geoscience problems to be addressed using numerical methods and computational modeling approaches. For example, computational modeling has been successfully used to solve ore-forming and mine site contamination/remediation problems, in which fluid flow and geochemical processes play important roles in the controlling dynamic mechanisms. The main purpose of this paper is to present a generalized overview of: (1) the various classes and models associated with fluid flow/chemically reacting systems in order to highlight possible opportunities and developments for the future; (2) some more general issues that need attention in the development of computational models and codes for simulating ore-forming and geoenviromental systems; (3) the related progresses achieved on the geochemical modeling over the past 50 years or so; (4) the general methodology for modeling of oreforming and geoenvironmental systems; and (5) the future development directions associated with modeling of ore-forming and geoenviromental systems.

Coupling of eastern and western subpolar North Atlantic: salt transport in the Irminger Current

Salt transport in the Irminger Current and thus the coupling between eastern and western subpolar North Atlantic plays an important role for climate variability across a wide range of time scales. High-resolution ocean modeling and observations indicate that 5 salinities in the eastern subpolar North Atlantic decrease with enhanced circulation of the North Atlantic subpolar gyre (SPG). This has led to the perception that a stronger SPG also transports less salt westward. In this study, we analyze a regional ocean model and a comprehensive global coupled climate model, and show that a stronger SPG transports more salt in the Irminger Current irrespective of lower salinities in its 10 source region. The additional salt converges in the Labrador Sea and the Irminger Basin by eddy transports, increases surface salinity in the western SPG, and favors more intense deep convection. This is part of a positive feedback mechanism with potentially large implications for climate variability and predictability.

New online method for water isotope analysis of speleothem fluid inclusions using laser absorption spectroscopy (WS-CRDS)

A new online method to analyse water isotopes of speleothem fluid inclusions using a wavelength scanned cavity ring down spectroscopy (WS-CRDS) instrument is presented. This novel technique allows us simultaneously to measure hydrogen and oxygen isotopes for a released aliquot of water. To do so, we designed a new simple line that allows the online water extraction and isotope analysis of speleothem samples. The specificity of the method lies in the fact that fluid inclusions release is made on a standard water background, which mainly improves the �� D robustness. To saturate the line, a peristaltic pump continuously injects standard water into the line that is permanently heated to 140 ��C and flushed with dry nitrogen gas. This permits instantaneous and complete vaporisation of the standard water, resulting in an artificial water background with well-known �� D and ��18O values. The speleothem sample is placed in a copper tube, attached to the line, and after system stabilisation it is crushed using a simple hydraulic device to liberate speleothem fluid inclusions water. The released water is carried by the nitrogen/standard water gas stream directly to a Picarro L1102-i for isotope determination. To test the accuracy and reproducibility of the line and to measure standard water during speleothem measurements, a syringe injection unit was added to the line. Peak evaluation is done similarly as in gas chromatography to obtain &delta D; and ��18O isotopic compositions of measured water aliquots. Precision is better than 1.5 ��� for �� D and 0.4 ��� for ��18O for water measurements for an extended range (���210 to 0 ��� for �� D and ���27 to 0 ��� for ��18O) primarily dependent on the amount of water released from speleothem fluid inclusions and secondarily on the isotopic composition of the sample. The results show that WS-CRDS technology is suitable for speleothem fluid inclusion measurements and gives results that are comparable to the isotope ratio mass spectrometry (IRMS) technique.

Palynological richness and evenness: insights from the taxa accumulation curve

Palynology provides the opportunity to make inferences on changes in diversity of terrestrial vegetation over long time scales. The often coarse taxonomic level achievable in pollen analysis, differences in pollen production and dispersal, and the lack of pollen source boundaries hamper the application of diversity indices to palynology. Palynological richness, the number of pollen types at a constant pollen count, is the most robust and widely used diversity indicator for pollen data. However, this index is also influenced by the abundance distribution of pollen types in sediments. In particular, where the index is calculated by rarefaction analysis, information on taxonomic richness at low abundance may be lost. Here we explore information that can be extracted from the accumulation of taxa over consecutive samples. The log-transformed taxa accumulation curve can be broken up into linear sections with different slope and intersect parameters, describing the accumulation of new taxa within the section. The breaking points may indicate changes in the species pool or in the abundance of high versus low pollen producers. Testing this concept on three pollen diagrams from different landscapes, we find that the break points in the taxa accumulation curves provide convenient zones for identifying changes in richness and evenness. The linear regressions over consecutive samples can be used to inter- and extrapolate to low or extremely high pollen counts, indicating evenness and richness in taxonomic composition within these zones. An evenness indicator, based on the rank-order-abundance is used to assist in the evaluation of the results and the interpretation of the fossil records. Two central European pollen diagrams show major changes in the taxa accumulation curves for the Lateglacial period and the time of human induced land-use changes, while they do not indicate strong changes in the species pool with the onset of the Holocene. In contrast, a central Swedish pollen diagram shows comparatively little change, but high richness during the early Holocene forest establishment. Evenness and palynological richness are related for most periods in the three diagrams, however, sections before forest establishment and after forest clearance show high evenness, which is not necessarily accompanied by high palynological richness, encouraging efforts to separate the two.

Ocean (De)oxygenation Across the Last Deglaciation: Insights for the Future

Anthropogenic warming is expected to drive oxygen out of the ocean as the water temperature rises and the rate of exchange between subsurface waters and the atmosphere slows due to enhanced upper ocean density stratification. Observations from recent decades are tantalizingly consistent with this prediction, though these changes remain subtle in the face of natural variability. Earth system model projections unanimously predict a long-term decrease in the global ocean oxygen inventory, but show regional discrepancies, particularly in the most oxygen-depleted waters, owing to the complex interplay between oxygen supply pathways and oxygen consumption. The geological record provides an orthogonal perspective, showing how the oceanic oxygen content varied in response to prior episodes of climate change. These past changes were much slower than the current, anthropogenic change, but can help to appraise sensitivities, and point toward potentially dominant mechanisms of change. Consistent with the model projections, marine sediments recorded an overall expansion of low-oxygen waters in the upper ocean as it warmed at the end of the last ice age. This expansion was not linearly related with temperature, though, but reached a deoxygenation extreme midway through the warming. Meanwhile, the deep ocean became better oxygenated, opposite the general expectation. These observations require that significant changes in apparent oxygen utilization occurred, suggesting that they will also be important in the future.

Swiss tree rings reveal warm and wet summers during medieval times

We present a 1200 year drought reconstruction for the European Alpine region based on carbon isotope variations of tree rings from living larch trees and historic timber. The carbon isotope fractionation at the study site is sensitive to summer precipitation, temperature, and irradiance, resulting in a stable and high correlation with a drought index for interannual to decadal frequencies and possibly beyond (r(2)=0.58 for 1901-2004, July/August). When combining this information with maximum latewood density-derived summer temperature, a strongly reduced occurrence of summer droughts during the warm A.D. 900-1200 period is evident, coinciding with the Medieval Climate Anomaly (MCA), with a shift to colder and drier conditions for the subsequent centuries. The warm-wet MCA contrasts strongly with the climate of the drought-prone warm phase of the recent decades, indicating different forcing mechanism for these two warm periods and pointing to beneficial conditions for agriculture and human well-being during the MCA in this region.

Iron Fertilization of the Subantarctic Ocean During the Last Ice Age

John H. Martin, who discovered widespread iron limitation of ocean productivity, proposed that dust-borne iron fertilization of Southern Ocean phytoplankton caused the ice age reduction in atmospheric carbon dioxide (CO2). In a sediment core from the Subantarctic Atlantic, we measured foraminifera-bound nitrogen isotopes to reconstruct ice age nitrate consumption, burial fluxes of iron, and proxies for productivity. Peak glacial times and millennial cold events are characterized by increases in dust flux, productivity, and the degree of nitrate consumption; this combination is uniquely consistent with Subantarctic iron fertilization. The associated strengthening of the Southern Ocean���s biological pump can explain the lowering of CO2 at the transition from mid-climate states to full ice age conditions as well as the millennial-scale CO2 oscillations.

A Lagrangian Analysis of the Northern Hemisphere Subtropical Jet

This study presents a 5-yr climatology of 7-day back trajectories started from the Northern Hemisphere subtropical jet. These trajectories provide insight into the seasonally and regionally varying angular momentum and potential vorticity characteristics of the air parcels that end up in the subtropical jet. The trajectories reveal preferred pathways of the air parcels that reach the subtropical jet from the tropics and the extratropics and allow estimation of the tropical and extratropical forcing of the subtropical jet. The back trajectories were calculated 7 days back in time and started every 6 h from December 2005 to November 2010 using the Interim European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-Interim) dataset as a basis. The trajectories were started from the 345-K isentrope in areas where the wind speed exceeded a seasonally varying threshold and where the wind shear was confined to upper levels. During winter, the South American continent, the Indian Ocean, and the Maritime Continent are preferred areas of ascent into the upper troposphere. From these areas, air parcels follow an anticyclonic pathway into the subtropical jet. During summer, the majority of air parcels ascend over the Himalayas and Southeast Asia. Angular momentum is overall well conserved for trajectories that reach the subtropical jet from the deep tropics. In winter and spring, the hemispheric-mean angular momentum loss amounts to approximately 6%; in summer, it amounts to approximately 18%; and in fall, it amounts to approximately 13%. This seasonal variability is confirmed using an independent potential vorticity���based method to estimate tropical and extratropical forcing of the subtropical jet.

Retrospective analysis of a nonforecasted rain-on-snow flood in the Alps ��� a matter of model limitations or unpredictable nature?

A rain-on-snow flood occurred in the Bernese Alps, Switzerland, on 10 October 2011, and caused significant damage. As the flood peak was unpredicted by the flood forecast system, questions were raised concerning the causes and the predictability of the event. Here, we aimed to reconstruct the anatomy of this rain-on-snow flood in the L��tschen Valley (160 km2) by analyzing meteorological data from the synoptic to the local scale and by reproducing the flood peak with the hydrological model WaSiM-ETH (Water Flow and Balance Simulation Model). This in order to gain process understanding and to evaluate the predictability. The atmospheric drivers of this rain-on-snow flood were (i) sustained snowfall followed by (ii) the passage of an atmospheric river bringing warm and moist air towards the Alps. As a result, intensive rainfall (average of 100 mm day-1) was accompanied by a temperature increase that shifted the 0�� line from 1500 to 3200 m a.s.l. (meters above sea level) in 24 h with a maximum increase of 9 K in 9 h. The south-facing slope of the valley received significantly more precipitation than the north-facing slope, leading to flooding only in tributaries along the south-facing slope. We hypothesized that the reason for this very local rainfall distribution was a cavity circulation combined with a seeder-feeder-cloud system enhancing local rainfall and snowmelt along the south-facing slope. By applying and considerably recalibrating the standard hydrological model setup, we proved that both latent and sensible heat fluxes were needed to reconstruct the snow cover dynamic, and that locally high-precipitation sums (160 mm in 12 h) were required to produce the estimated flood peak. However, to reproduce the rapid runoff responses during the event, we conceptually represent likely lateral flow dynamics within the snow cover causing the model to react "oversensitively" to meltwater. Driving the optimized model with COSMO (Consortium for Small-scale Modeling)-2 forecast data, we still failed to simulate the flood because COSMO-2 forecast data underestimated both the local precipitation peak and the temperature increase. Thus we conclude that this rain-on-snow flood was, in general, predictable, but requires a special hydrological model setup and extensive and locally precise meteorological input data. Although, this data quality may not be achieved with forecast data, an additional model with a specific rain-on-snow configuration can provide useful information when rain-on-snow events are likely to occur.

The influence of sampling design on tree-ring-based quantification of forest growth

Tree-rings offer one of the few possibilities to empirically quantify and reconstruct forest growth dynamics over years to millennia. Contemporaneously with the growing scientific community employing tree-ring parameters, recent research has suggested that commonly applied sampling designs (i.e. how and which trees are selected for dendrochronological sampling) may introduce considerable biases in quantifications of forest responses to environmental change. To date, a systematic assessment of the consequences of sampling design on dendroecological and-climatological conclusions has not yet been performed. Here, we investigate potential biases by sampling a large population of trees and replicating diverse sampling designs. This is achieved by retroactively subsetting the population and specifically testing for biases emerging for climate reconstruction, growth response to climate variability, long-term growth trends, and quantification of forest productivity. We find that commonly applied sampling designs can impart systematic biases of varying magnitude to any type of tree-ring-based investigations, independent of the total number of samples considered. Quantifications of forest growth and productivity are particularly susceptible to biases, whereas growth responses to short-term climate variability are less affected by the choice of sampling design. The world's most frequently applied sampling design, focusing on dominant trees only, can bias absolute growth rates by up to 459% and trends in excess of 200%. Our findings challenge paradigms, where a subset of samples is typically considered to be representative for the entire population. The only two sampling strategies meeting the requirements for all types of investigations are the (i) sampling of all individuals within a fixed area; and (ii) fully randomized selection of trees. This result advertises the consistent implementation of a widely applicable sampling design to simultaneously reduce uncertainties in tree-ring-based quantifications of forest growth and increase the comparability of datasets beyond individual studies, investigators, laboratories, and geographical boundaries.

The Yanque prospect (Peru): from polymetallic Zn-Pb mineralization to a nonsulfide deposit

The Yanque nonsulfide Pb-Zn deposit (inferred resources 12.5 Mt @ 3.7% Pb and @ 3.5% Zn) is located in the Andahuaylas-Yauri ore province (Cuzco, southern Peru). The deposit occurs within a base metal mineralized district, centered on the medium-sized Dolores porphyry copper. A thorough geological, mineralogical and geochemical study has carried out in order to define: the relationships between the Dolores Cu-porphyry ore and the Yanque Zn-Pb polymetallic mineralization, and the characteristics of the economic nonsulfide concentrations. Both sedimentary and igneous rocks constitute the backbone of the Yanque-Dolores area. The sedimentary lithologies belong to the Soraya, Mara and Ferrobamba Fms. (upper Jurassic-middle Cretaceous). The Yanque orebody is hosted by the Mara Fm., which prevailingly consists of a siliciclastic sedimentary breccia. The original sulfide mineralization consisted of galena, pyrite and sphalerite. The host rock has been affected by a strong hydrothermal alteration, characterized by prevailing sericite/illite, as in the typical porphyry-related phyllic-argillic alteration stage, and by minor kaolinite, dolomite and quartz. Minor element geochemistry, characterized by Sb, As, Mn, Ag and locally also by Cu, points to magmatic-hydrothermal related mineralizing fluids. The Pb isotopic compositions from Dolores and Yanque sulfides are similar, and are typical of the Tertiary magmatically-derived ores in this part of Peru. The hydrothermally altered rocks at Yanque have the same Pb isotopic compositions as the sulfides, thus confirming the hypothesis that the Yanque primary Zn-Pb mineralization may have been produced by hydrothermal circulation related to the emplacement of the Dolores Cu-porphyry, as it is the case of other porphyry Cu systems associated with polymetallic mineralization elsewhere. However, no simple genetic model for the mineralization involving just one fluid circulation episode is able to explain the data. The Yanque economic nonsulfide ore association consists of sauconite, hemimorphite, smithsonite and cerussite, which result from the weathering and alteration of the original sulfide mineralization. Zinc is allocated mainly in sauconite (Zn-smectite), rather than in carbonates: a factor strictly related to the prevailing siliciclastic character of the host rock. Distinctive features of the Yanque orebody are the comparable ore grades for both Pb and Zn (3.5% Zn and 3.7% Pb), and the inverse supergene chemical zoning. In fact, contrary to other supergene ores of this type, zinc prevails in the top zone of the Yanque deposit, whereas lead content increases with depth. Considering the different mobility of the two metals in solution, it may be assumed that most of the primary zinc that was the source for the Yanque nonsulfides was originally located far from the position occupied by the galena mineralization, whose remnants have been observed on site. Zinc sulfides may have been originally contained in the eroded rock volumes that surrounded the actual deposit: the zinc-rich solutions have possibly migrated through the siliciclastic Mara Fm. and precipitated the nonsulfide minerals by porosity filling and replacement processes. In this sense, the Yanque secondary Zn-Pb deposit could be considered as a special type of ���Exotic��� mineralization.

A re-examination of petrogenesis and 40Ar/39Ar systematics in the the Chain of Ponds K-feldspar: "diffusion domain" archetype versus polyphase hygrochronology

K-feldspar (Kfs) from the Chain of Ponds Pluton (CPP) is the archetypal reference material, on which thermochronological modeling of Ar diffusion in discrete ���domains��� was founded. We re-examine the CPP Kfs using cathodoluminescence and back-scattered electron imaging, transmission electron microscopy, and electron probe microanalysis. 40Ar/39Ar stepwise heating experiments on different sieve fractions, and on handpicked and unpicked aliquots, are compared. Our results reproduce the staircase-shaped age spectrum and the Arrhenius trajectory of the literature sample, confirming that samples collected from the same locality have an identical Ar isotope record. Even the most pristine-looking Kfs from the CPP contains successive generations of secondary, metasomatic/retrograde mineral replacements that post-date magmatic crystallization. These chemically and chronologically distinct phases are responsible for its staircase-shaped age spectra, which are modified by handpicking. While genuine within-grain diffusion gradients are not ruled out by these data, this study demonstrates that the most important control on staircase-shaped age spectra is the simultaneous presence of heterochemical, diachronous post-magmatic mineral growth. At least five distinct mineral species were identified in the Kfs separate, three of which can be traced to external fluids interacting with the CPP in a chemically open system. Sieve fractions have size-shifted Arrhenius trajectories, negating the existence of the smallest ���diffusion domains���. Heterochemical phases also play an important role in producing non-linear trajectories. In vacuo degassing rates recovered from Arrhenius plots are neither related to true Fick���s Law diffusion nor to the staircase shape of the age spectra. The CPP Kfs used to define the "diffusion domain" model demonstrates the predominance of metasomatic alteration by hydrothermal fluids and recrystallization in establishing the natural Ar distribution amongst different coexisting phases that gives rise to the staircase-shaped age spectrum. Microbeam imaging of textures is as essential for 40Ar-39Ar hygrochronology as it is for U-Pb geochronology.

Dating polygenetic metamorphic assemblages along a transect across the Western Alps

Multichronometric analyses were performed on samples from a transect in the French-Italian Western Alps crossing nappes derived from the Brian��onnais terrane and the Piemonte-Liguria Ocean, in an endeavour to constrain the high-pressure (HP) metamorphism and the retrogression history. 12 samples of white mica were analysed by 39Ar-40Ar stepwise heating, complemented by 2 samples from the Monte Rosa 100 km to the NE and also attributed to the Brian��onnais terrane. One Sm-Nd and three Lu-Hf garnet ages from eclogites were also obtained. White mica ages decrease from ca. 300 Ma in the westernmost samples (Zone Houill��re), reaching ca. 300 ��C during Alpine metamorphism, to < 48 Ma in the internal units to the East, which reached ca. 500 ��C during Alpine orogeny. The conventional ���thermochronological��� interpretation postulates Cretaceous Eo-Alpine HP metamorphism and younger ���cooling ages��� in the higher-temperature samples. However, Eocene Lu-Hf and Sm-Nd ages from the same samples cannot be interpreted as post-metamorphic cooling ages, which makes a Cretaceous eclogitization untenable. The age date from this transect require instead to replace conventional ���thermochronology��� by an approach combining age dating with detailed geochemical, petrological and microstructural investigations. Petrology reveals important mineralogical differences along the transect. Samples from the Zone Houill��re mostly contain detrital mica. White mica with Si > 6.45 atoms per formula unit becomes more abundant eastward. Across the whole traverse, HP phengitic mica forms the D1 foliation. Syn-D2 mica is Si-poorer and associated with nappe stacking, exhumation, and hydrous retrogression under greenschist facies conditions. D1 phengite is very often corroded, overgrown or intergrown by syn-D2 muscovite. Most importantly, syn-D2 recrystallization is not limited to S2 schistosity domains; microchemical fingerprinting shows that it also can form pseudomorphs after crystals that could be mistaken to have formed during D1 based on microstructural arguments alone. Thereby the Cl concentration in white mica is a useful discriminator, since D2 retrogression was associated with a less saline fluid than eclogitization. Once the petrological stage is set, geochronology is straightforward. All samples contain mixtures of detrital, syn-D1 and syn-D2 mica, and retrogression phases (D3) in greatly varying proportions according to local pressure-temperature-fluid activity-deformation conditions. The correlation of age vs. Cl/K clearly identifies 47 �� 1 Ma as the age of formation of syn-D1 mica along the entire transect, including the Monte Rosa nappe samples. The inferred age of the greenschist-facies low-Si syn-D2 mica generation ranges within 39-43 Ma, with local variations. Coexistence of D1 and D2 ages, and the constancy of non-reset D1 ages along the entire transect, are strong evidence that the D1 white mica ages are very close to formation ages. Volume diffusion of Ar in white mica (activation energy E = 250 kJ/mol; pressure-adjusted diffusion coefficient D���0 < 0.03 cm2 s-1) has a subordinate effect on mineral ages compared to both prograde and retrograde recrystallization in most samples. Eocene Lu-Hf and Sm-Nd garnet ages are prograde and predate the HP peak.

Provenance versus weathering control on the composition of tropical river mud (southern Africa)

This study presents an integrated mineralogical-geochemical data base on fine-grained sediments transported by all major rivers of southern Africa, including the Zambezi, Okavango, Limpopo, Olifants, Orange and Kunene. Clay mineralogy, bulk geochemistry, Sr and Nd isotopic signatures of river mud, considered as proxy of suspended load, are used to investigate the influence of source-rock lithology and weathering intensity on the composition of clay and silt produced in subequatorial to subtropical latitudes. Depletion in mobile alkali and alkaline-earth metals, minor in arid Namibia, is strong in the Okavango, Kwando and Upper Zambezi catchments, where recycling is also extensive. Element removal is most significant for Na, and to a lesser extent for Sr. Depletion in K, Ca and other elements, negligible in Namibia, is moderate elsewhere. The most widespread clay minerals are smectite, dominant in muds derived from Karoo or Etendeka flood basalts, or illite and chlorite, dominant in muds derived from metasedimentary rocks of the Damara Orogen or Zimbabwe Craton. Kaolinite represents 30-40% of clay minerals only in Okavango and Upper Zambezi sediments sourced in humid subequatorial Angola and Zambia. After subtracting the effects of recycling and of local accumulation of authigenic carbonates in soils, the regional distribution of clay minerals and chemical indices consistently reflect weathering intensity primarily controlled by climate. Bulk geochemistry identifies most clearly volcaniclastic sediments and mafic sources in general, but cannot discriminate the other sources of detritus in detail. Instead, Sr and Nd isotopic fingerprints are insensitive to weathering, and thus mirror faithfully the tectonic structure of the southern African continent. Isotopic tools thus represent a much firmer basis than bulk geochemistry or clay mineralogy in the provenance study of mudrocks.

Elemental and Lead Isotopic Data of Copper Finds from the Singen Cemetery, Germany ��� a Methodological Approach to Investigate Early Bronze Age Trade Networks

We report a trace element - Pb isotope analytical (LIA) database on the "Singen Copper", a peculiar type of copper found in the North Alpine realm, from its type locality, the Early Bronze Age Singen Cemetery (Germany). What distinguishes ���Singen Copper��� from other coeval copper types? (i) is it a discrete metal lot with a uniform provenance (if so, can its provenance be constrained)? (ii) was it manufactured by a special, unique metallurgical process that can be discriminated from others? Trace element concentrations can give clues on the ore types that were mined, but they can be modified (more or less intentionally) by metallurgical operations. A more robust indicator are the ratios of chemically similar elements (e.g. Co/Ni, Bi/Sb, etc.), since they should remain nearly constant during metallurgical operations, and are expected to behave homogeneously in each mineral of a given mining area, but their partition amongst the different mineral species is known to cause strong inter-element fractionations. We tested the trace element ratio pattern predicted by geochemical arguments on the Brixlegg mining area. Brixlegg itself is not compatible with the Singen Copper objects, and we only report it because it is a rare instance of a mining area for which sufficient trace element analyses are available in the literature. We observe that As/Sb in fahlerz varies by a factor 1.8 above/below median; As/Sb in enargite varies by a factor of 2.5 with a 10 times higher median. Most of the 102 analyzed metal objects from Singen are Sb-Ni-rich, corresponding to ���antimony-nickel copper��� of the literature. Other trace element concentrations vary by > 100 times, ratios by factors > 50. Pb isotopic compositions are all significantly different from each other. They do not form a single linear array and require > 3 ore batches that certainly do not derive from one single mining area. Our data suggest a heterogeneous provenance of ���Singen copper���. Archaeological information limits the scope to Central European sources. LIA requires a diverse supply network from many mining localities, including possibly Brittany. Trace element ratios show more heterogeneity than LIA; this can be explained either by deliberate selection of one particular ore mineral (from very many sources) or by processing of assorted ore minerals from a smaller number of sources, with the unintentional effect that the quality of the copper would not be constant, as the metallurgical properties of alloys would vary with trace element concentrations.