GEOLOGIJA 46/1, 129–134, Ljubljana 2003
Organic geochemical records of hydrothermal alteration at Idrija
mercury deposit, Slovenia
Organski geokemi~ni zapisi hidrotermalnih sprememb v rudi{~u ‘ivega srebra
Idrija, Slovenija
Jo{t V. LAVRI^ 1, Jorge E. SPANGENBERG1 & Bojan RE@UN2
1Institut de Minéralogie et Géochimie, Université de Lausanne, BFSH2, CH-1015 Lausanne,
Switzerland 2Rudnik `ivega srebra Idrija v zapiranju d.o.o., Arkova 43, SI-5280 Idrija, Slovenia
Key words: stable isotopes, PAH, organic matter, hydrothermal alteration, Idrija, Slovenia
Klju~ne besede: stabilni izotopi, policikli~ni aromatski ogljikovodiki, organska snov, hidrotermalne spremembe, Idrija, Slovenija
Abstract
A combined molecular and stable isotopes (C, N) investigation of the organic matter associated to the Idrija mercury deposit is being used to track the pathway of the mineralizing hydrothermal fluids and their interaction with the organic matter. The studied samples include regional barren rocks, and host rocks and ore from the mine. The Rock-Eval parameters indicate that the organic matter disseminated in regional and mine rocks is mature to post mature. The bitumens from mineralized samples are depleted in aliphatic hydrocarbons and enriched in polycyclic aromatic hydrocarbons (PAH), aromatic sulfur compounds (S-PAH) and hydrogenated PAH. The isotopic compositions of kerogens, bitumens, and individual hydrocarbons point to thermal and oxidative degradation of indigenous and migrated bitumens during mineralization and later evolution of the Idrija deposit.
Kratka vsebina
Kombinirana molekularna in stabilno izotopska (C, N) raziskava organske snovi Idrijskega ‘ivosrebrovega rudi{~a je bila uporabljena za sledenje poti hidrotermalnih raztopin in njihovega vpliva na organsko snov. Preiskani vzorci zajemajo regionalne kamnine ter prikamnino in rudo iz Idrijskega rudnika. Parametri Rock-Eval ka‘ejo na zrelost do prezrelost organske snovi, ki je razpr{ena v idrijskih in okolnih kamninah. Bitumni iz mineraliziranih vzorcev nakazujejo zmanj{ane vsebnosti nasi~enih ogljikovodikov ter po-ve~ane koncentracije policikli~nih aromatskih ogljikovodikov (PAH), aromatskih ‘veplo-vih spojin (S-PAH) in hidrogeniranih PAH. Izotopske sestave kerogenov, bitumnov in posameznih ogljikovodikov pri~ajo o termi~ni in oksidativni degradaciji prvotnih in mi-griranih bitumnov med mineralizacijo in poznej{im razvojem idrijskega rudi{~a.
Introduction
The Idrija world-class mercury deposit is located in western Slovenia, about 50 km west of Ljubljana (Fig. 1). As the world second largest mercury mine after Almadén (Spain), Idrija has produced more than 12.7 Mt ore with 145,000 t Hg since 1490 (Mlak a r , 1974). In 1988 the Idrija mine stopped production and initiated a shutdown program, which shall finish in 2006.
The geology and the genesis of the Idrija deposit were described in numerous studies (e.g. Mlakar, 1967; Mlakar & Drove-n i k , 1971; Placer , 1982; ^ ar , 1975, 1990). The geochemical investigations include trace element concentrations in host rocks (^ a -d e ` et al., 1981) and cinnabar (Berce, 1958;
Drovenik et al., 1980), mercury contents at deposit scale (Berce, 1965), stable isoto-pic studies of ore and host rocks (Ozerova et al., 1973; Drovenik et al., 1976, 1991), and mineralogical, molecular and isotopic characterization of the polycyclic aromatic hydrocarbon (PAH) mineral idrialite associated to the Idrija ore (S t r u n z & C o n t a g , 1965; Blumer, 1975; Wise et al., 1986; Span genberg et al., 1999). However, in order to better understand the mineralization process and in its relationship with the organic matter, further geochemical investigations are necessary. The ongoing organic and inorganic geochemical study includes trace elements of the host rocks, the isotopic composition of carbonate and sulfur minerals (?13C, ?18O, ?34S), and molecular and iso-
130
Jo{t V. Lavri~, Jorge E. Spangenberg & Bojan Re`un
Figure 1. Gas chromatograms of (a) aliphatic and (b) aromatic hydrocarbons extracted from the Upper Scythian dolostone. The n-alkanes are indicated by their carbon
numbers. Pr = pristane; Ph =
phytane; UCM = unresolved complex
mixture. Compounds corresponding
to numbered peaks are listed in
Table 1
Slika 1. Plinski kromatogrami (a)
nasi~enih in (b) aromatskih ogljikovodikov iz zgornjeskitskega
dolomita. Normalni alkani so
ozna~eni s {tevilkami, ki pomenijo
{tevilo ogljikovih atomov Pr =
pristan; Ph = fitan; UCM =
nerazlo~ena kompleksna me{anica.
Spojine, ki ustrezajo o{tevil~enim
kromatogramskim vrhovom, so
navedene v tabeli 1
topic (?D, ?13C, ?15N) composition of the associated organic matter (e.g. Lavri~ & Spangenberg, 2001, 2002). Here we present Rock-Eval pyrolysis parameters, ?13C and ?15N of the kerogens (insoluble organic matter), ?13C and the molecular composition of the bitumens (extractable organic matter), ?13C of the individual n-alkanes, and ?13C of pyrobitumen. These data provide insights into the nature of the hydrothermal alteration and remobilization of the Idrija organic matter, and help to constrain the ore fluids pathway.
Geological setting and mineralization
The stratigraphic succession at the Idrija deposit comprises about 5,500 m of sedimentary rocks of Permocarboniferous to Eocene age, from which the lowermost 800 m host the mineralization (Mlakar & Dro-venik, 1971). The Idrija structure developed as a part of an E-W extending failed rift (Idrija graben) during the Middle Triassic intra-continental rifting (Placer & ^ar, 1977). The host rocks were affected by two main post-ore deformation episodes: (1) an Early Tertiary folding and thrusting of about 30 km toward SSW, and (2) a Late Tertiary dextral NW-SE strike-slip faulting, with a horizontal displacement of up to 2.5 km (Mlakar, 1969; Placer, 1982).
The mineralization took place in two phases during the Scythian-Ladinian rifting and bimodal volcanism from near neutral mer-curiferous hydrothermal fluids channeled by a system of deep subvertical faults within the Idrija graben (Mlakar & Drovenik , 1971). The ore occurs as syngenetic stratiform bodies in the Upper Ladinian Skonca beds and the overlying tuffs, and as epigene-tic open space filling and replacement in fault zones in Permocarboniferous to Upper Ladinian beds. The ore consists of cinnabar and native mercury, with minor pyrite, mar-casite and metacinnabar. Up to 1.5 m thick lenses of evaporites (gypsum and anhydrite) occur in the Upper Permian and Lower Scythian dolostones (^ade`, 1977). The main gangue minerals are quartz, calcite and dolomite, with rare barite and fluorite (M l a-kar & Drovenik , 1971). Three distinct types of organic matter occur at Idrija: (1) kerogen and bitumen in the host rocks and ore, (2) open-space filling black solid pyro-bitumen, and (3) idrialite intergrown with the mercury ore.
Sampling
Permocarboniferous to Upper Ladinian barren and mineralized lithologies were sampled at outcrops up to 6 km from the deposit (n = 56), and from the mine walls at
Organic geochemical records of hydrothermal alteration at Idrija mercury deposit, Slovenia            131
different levels of the Idrija mine (n = 125). The hydrothermal alteration and flow direction of the mineralizing fluids on mine scale were investigated in profiles across mineralized and barren zones. The largest part of the host rock and ore samples (n = 111) was collected in seven profiles at level IV.
Analyses
Rock powders of selected regional (n = 18) and mine (n = 45) samples were submitted to the determination of total organic carbon (TOC) and Rock-Eval pyrolysis. The organic matter was studied following the procedures described in Spangenberg & Macko (1998). The C and N stable isotopic composition of the kerogens, and the ?13C of bulk bitumens and pyrobitumens were determined by using a Carlo Erba 1108 elemental analyzer connected to a Thermo Finnigan Delta S isotope ratio mass spectrometer (EA/IRMS). The isotopic data are reported in the delta (d) notation as the per mil (‰) deviations relative to the VPDB and AIR for carbon and nitrogen, respectively. The reproducibility of the EA/ IRMS measurements for carbon and nitrogen is better than ±0.1 and ±0.3‰ (1s), respectively. The aliphatic and aromatic fractions of the bitumens were chemically characterized by an Agilent Technologies 6890 gas chromatograph coupled to a 5973N mass selective detector (GC/MSD). The C stable isotopic composition of the individual n-alkanes was determined in triplicate by using the GC coupled to the Delta S IRMS by a combustion (C) interface III (GC/C/IRMS). The reproducibility of the GC/C/IRMS analyses of the Idrija n-alkanes ranges between 0.1 and 1.0‰.
Results
Bulk organic geochemical and isotopic data The TOC content for mine samples ranges from 0.01 to 1.47 wt%, with median values of 0.45 wt% for the Permocarboniferous shales, 0.25 wt% for the Upper Permian dolostone, 0.13 wt% for the Lower Scythian to Anisian dolostones, and 0.85 wt% for the Upper Ladinian Skonca beds. For regional samples the TOC ranges from 0 to 1.75 wt%, with median values of 0.61 wt% for the Permocarboniferous shales, 0.31 for the Up-
per Permian dolostone, 0.03 wt% for the Lower Scythian to Anisian dolostones, and 0.98 wt% for the Upper Ladinian Skonca beds. The Rock-Eval S1 and S2 peaks are generally small (< 0.2 mg HC/g) or absent, which indicates a low petroleum-generative potential, and makes the S2-derived temperatures (Tmax) unreliable (Peters, 1986). Exceptions are one regional and one mine sample of the Upper Ladinian Skonca shale, having S2 peaks of 0.27 and 0.41 mg HC/g, and Tmax values of 524 and 504°C, respectively. The samples have low hydrogen and high oxygen indices (HI < 29 mg HC/g TOC, OI < 267 mg CO2/g TOC), which is typical of highly recycled and oxidized kerogens. An additional decrease of the HI, and increase of the OI due to the mineral matrix effect have to be taken into account for samples with TOC < 1.5 wt% (e.g. Hunt, 1996). The C and N isotopic compositions of the kerogens (n = 12) range from –29.0 to –23.1‰ and from –1.2 to +4.6‰, respectively. The kerogens of mineralized samples are enriched in 15N (up to ~1‰) and in 13C (up to 3.6‰) compared to barren or slightly mineralized samples. In both mineralized and barren mine samples the kerogens are isoto-pically lighter (–29.0 to –23.1‰, median = – 26.3‰) compared to the associated bitumens (–27.7 to –22.8‰, median = –25.2‰). The pyro-bitumens are isotopically heavier (?13C ? – 23.0‰), reflecting their origin from thermally altered migrated hydrocarbons.
Molecular organic geochemistry
The main resolvable compounds in the aliphatic hydrocarbons (HC) fraction of the Idrija bitumens are unimodally-distributed n-alkanes in the C12 to C30 range and the acyclic isoprenoids pristane (Pr) and phyta-ne (Ph). Trace amounts of the biomarker hydrocarbons hopanes and steranes are present. Regional samples and barren mine samples have higher concentrations of aliphatic HC as the mineralized samples. The mineralized samples show a larger hump of unresolved complex mixture (UCM) of light hydrocarbons moieties, and higher concentrations of lower molecular weight n-alka-nes compared to barren samples (Fig. 1). The molecular parameters Pr/Ph (0.36 to 3.94), Pr/n-C17 (0.10 to 1.83), and Ph/n-C18 (0.06 to 3.31) are highly variable, and independent of the degree of mineralization. Two distinct trends of the d13C values of the individual n-
132
Jo{t V. Lavri~, Jorge E. Spangenberg & Bojan Re`un
Peak   Compound
Peak     Compound
1           Methylnaphtalene
2           Biphenyl
3           Methylbiphenyl
4           Bibenzyl
5           Fluorene
6           Dimethylbiphenyl
7           Dimethylmethylphenylbenzene
8           Dihydroanthracene
9           Methylfluorene
10         Dihydrotrimethylphenylindene
11         Dibenzothiophene
12         Phenanthrene/Anthracene
13         Methyldibenzothiophene
14         Methylphenanthrene
15         Tetrahydrofluoranthene
16         Fluoranthene
17           Pyrene
18           Tetrahydrobenzofluorene
19           Benzofluorene
20           Benzonaphtothiophene
21           Tetrahydrobenzanthracene
22           Tetrahydrochrysene
23           Chrysene/Triphenylene
24           Methylbenzonaphtothiophene
25           Benzofluoranthene
26           Perylene
27           Cyclohexylbenzene
28           Ethenylnaphtalene
29           Cyclohexylmethylbenzene
30           Dimethyldibenzothiophene
31           Dimethylphenanthrene
Table 1. Aromatic hydrocarbons identified in Figures 1 and 2 Tabela 1. Aromatski ogljikovodiki, identificirani v slikah 1 in 2
alkanes were observed. In mineralized samples the C=17 n-alkanes are up to 3‰ lighter compared to the longer C-chain homologu-es, whereas in barren and slightly mineralized samples short and long-chain n-alkanes are enriched in 13C (1 –27‰), with the lowest values at n-C18 and n-C19 (~ –33‰).
The aromatic fraction of bitumens from mineralized samples is characterized by up to 5-rings polycyclic aromatic hydrocarbons (PAH), aromatic sulfur compounds (S-PAH; e.g. di-benzothiophene, benzonaphtothiophene) and their alkylated homologues (Table 1). Furthermore, the mineralized samples contain substantial concentrations of hydrogenated PAH (e.g., dihydroanthracene, tetrahydrochrysene), which occur in smaller concentrations or are not detected in barren mine samples and regional samples (Fig. 1). The PAH content increases with the degree of mineralization. Regional and barren mine samples have lower content in aromatic hydrocarbons which are mainly represented by lower molecular weight compounds (e.g. benzenes, naphtalenes). Exceptions are the barren samples from the Permocarboniferous shale and the Upper Permian dolostone which show substantial PAH concentrations (Fig. 2).
Discussion
The TOC and Rock-Eval data indicate that the Idrija host rocks, which were sub-
jected to prolonged hydrothermal activity, burial and extensive tectonic deformation, are organic-lean (TOC < 1 wt.%). The Tmax of the Upper Ladinian Skonca shale samples (504 and 557°C) suggest that they are post mature and may have been heated at temperatures of up to ~ 190°C (e.g. Hunt , 1996). This is in line with the homogenization temperatures of fluid inclusions in single quartz crystals from the Gruebler orebody of Idrija deposit (160 to 218°C; Palinka{ et al., 2001).
The relatively wide scatter of the kerogen C and N isotopic compositions can be attributed to (1) distinct primary organic sources in the different lithostratigraphic units, and (2) different degrees of alteration. The overall 13C enrichment of bitumens, associated to isotopically lighter kerogens, can be explained by a combination of thermal maturation and oxidative degradation of the indigenous HC on deposit scale during and after mineralization. A contribution of migrated solid bitumen (pyrobitumen) explains the additional 13C enrichment in the mineralized samples of the different litho-stratigraphic units. The variations of the hydrocarbons distribution (Pr/Ph, Pr/n-C17, Ph/n-C18) and d13C values of the individual n-alkanes indicate different maturation degrees and enhanced thermal alteration (e.g. cracking of longer chain n-alkanes, water washing) along the hydrothermal fluids pathways.
Organic geochemical records of hydrothermal alteration at Idrija mercury deposit, Slovenia            133
Figure 2. Gas chromatograms of the aromatic
hydrocarbons extracted from Upper Permian
dolostone. Compounds corresponding to
numbered peaks are listed in Table 1
Slika 2. Plinski kromatogrami aromatskih
ogljikovodikov iz zgornjepermskega dolomita.
Spojine, ki ustrezajo o{tevil~enim
kromatogramskim vrhovom, so navedene v
tabeli 1
PAH at Idrija hydrothermal system may have been formed through: (1) pyrolytic fragmentation of organic compounds followed by reformation (cyclization, aromatization, annelation), and (2) transformation of biologic precursors (e.g. saturated cyclic compounds) through dehydrogenation and dealk-ylation (e.g. M cCollom et al., 1999). The efficiency of both processes increases with temperature. Further aromatization of Idrija hydrocarbons may be related to S-cataly-zed reactions (e.g. Hunt , 1996) during ther-mochemical sulfate reduction which was one of the processes involved in the precipitation of cinnabar (Lavri~ & Spangenberg, 2002). The high content of PAH in the barren Permocarboniferous shales and Upper Permian dolostones suggest that these rocks most likely were the source of the aromatic petroleum staining Idrija ore and host rocks. Hydrous pyrolysis experiments showed that PAH could be hydrogenated at elevated temperatures (330°C) in an aqueous environment (M c C ollom et al., 1999). Thus, the occurrence of hydrogenated PAH in the Idrija ore samples indicates the degree of alteration and pathway of the hydrothermal fluids. The S-PAH are incorporated into the sedimentary organic matter during early diagenesis (e.g. O rr & Sinninghe Damsté, 1990). In hydrothermal deposits a part of the sulfur can be incorporated in the PAH structure during mineralization (e.g. Landais & Gize , 1997). Thus, the high concentrations of S-PAH in Idrija mineralized samples were formed during mineralization, probably related to thermally mediated reduction of sulfate by organic matter (Fig. 2)
Conclusions
The molecular and isotopic characteristics of the Idrija organic matter reflect migration, thermal maturation, and oxidation, which were enhanced by the mineralizing hydrothermal and post-ore fluids in the mineralized zones. The abundance and distribution patterns of hydrocarbons indicate that the organic matter in the Permocarbo-niferous shales and Upper Permian dolosto-nes was the main source of the aromatic hydrothermal petroleum associated to the Idrija ore. The organic geochemistry confirms the fracture-controlled nature of the
134
Jo{t V. Lavri~, Jorge E. Spangenberg & Bojan Re`un
Idrija hydrothermal system, previously established by the inorganic geochemical and isotopic patterns.
Acknowledgements
This research was supported by the Swiss National Science Foundation (grant 2100-059198.99/1) and the University of Lausanne. We thank J. ^ ar (University of Ljubljana) for field assistance and discussions on the geology of the Idrija deposit.
References
Berce , B. 1958: Geology of the Idrija mercury deposit.– Geologija, 4, 5–62, Ljubljana.
Berce , B. 1965: The use of mercury in geoche-mical prospecting for mercury .– Econ. Geol., 60, 1516–1528, Laucaster.
Blumer , M. 1975: Curtisite, idrialite and pen-dletonite, polycyclic aromatic hydrocarbon minerals: their composition and origin.– Chem. Geol., 16, 245–256, Amsterdam
^ade`, F. 1977: Gypsum and anhydrite occurrences in Idria region. – Geologija, 20, 289– 301, Ljubljana.
^ a d e ` , F., H u d n i k , V. & G o g a l a , A. 1981: Trace elements in the Idrija ore deposit country rocks.– Mining and metallurgy quarterly, 28, 33– 43, Ljubljana.
^ a r , J. 1975: Olisthostromes in the Idrija Middle Triassic trough–fault – Geologija, 18, 157–183, Ljubljana.
^ a r , J. 1990: Angular tectonic-erosional unconformity in the deposit part of the Idrija Middle Triassic tectonic structure. – Geologija, 31–32, 267–284, Ljubljana.
D r o v e n i k , M., D u h o v n i k , J. & P e z d i ~ , J. 1976: The sulfur isotope composition of sulfides from ore deposits in Slovenia. – Mining and metallurgy quarterly, 24, 193–46, Ljubljana.
D r o v e n i k , M., P l e n i ~ a r , M. & D r o v e n i k , F. 1980: The origin of Slovenian ore deposits.– Geologija, 23, 1–157, Ljubljana.
D r o v e n i k , M., D o l e n e c , T., Re`un, B. & P e z d i ~ , J. 1991: On the mercury ore from the Grübler orebody, Idrija – Geologija, 33, 397–446, Ljubljana.
H u n t , J. M. 1996: Petroleum geochemistry and geology. – 2nd edn. Freeman, 743 pp., New York.
L a n d a i s , P. & G i z e , A. P. 1997: Organic matter in hydrothermal ore deposits. In: H. L. Barnes (ed.), Geochemistry of hydrothermal ore deposits. – 3rd edn. Wiley, 613–655,
New York.
L a v r i ~ , J. V. & S p a n g e n b e r g , J. E. 2001: Geochemistry of the Idrija mercury deposit, Slovenia: insights from d13C, d18O, d34S and organic geochemical data. In: A. Piestrzynski et al. (eds.), Mineral deposits at the beginning of the 21st century. – Swets & Zeitlinger, 55–58, Lisse.
L a v r i ~ , J. V. & S p a n g e n b e r g , J. E. 2002: Stable isotope (C, O, S) systematics of the mer-
cury mineralization at Idrija, Slovenia: constraints on ore fluid source and alteration processes.– Miner. Deposita, submitted.
McCollom, T. M., Simoneit, B. R. T. & Shock, E. L. 1999: Hydrous pyrolysis of polyc-yclic aromatic hydrocarbons and implications for the origin of PAH in hydrothermal petroleum.– Energy & Fuels 13, 401–410, Washington.
M l a k a r , I. 1967: Relations between the lower and the upper structure of the Idrija ore deposit.– Geologija, 10, 87–126, Ljubljana.
M l a k a r , I. 1969: Nappe structure of the Idri-ja-@iri region.– Geologija, 12, 5–72, Ljubljana.
M l a k a r , I. & D r o v e n i k , M. 1971: Structural and genetic particularities of the Idrija mercury ore deposit. – Geologija, 14, 67–126, Ljubljana.
M l a k a r , I. 1974: An outline of production of the Idrija mercury mine through centuries.– Idrijski razgledi, 3–4, XIX, 1–40, Idrija.
Mlakar, I. 1996: On the Marija Reka mercury deposit and its comparison with the Litija and Idrija deposits from the aspect of plate tectonics. – Geologija, 37–38, 321–376, Ljubljana.
Orr, W. L. & Sinninghe Damsté, J. S. (1990) Geochemistry of sulfur in petroleum systems. In: W. L. Orr & C. M. White (eds.), Geochemistry of sulfur in fossil fuels. – ACS Symp. Ser. 429, 2–29, Washington.
O z e r o v a , N. A., V i n o g r a d o v , V. I., M l a -k a r , I., F e d o r c h u k , V. P. & T i t o v , I. N. 1973: Isotopic composition of the sulfur of some ores from ore deposits in the western part of the Mediterranean mercury belt (in Russian). In: F. K. Shipulin & K. M. Feodotiev (eds.) Abstract of the Geochemistry of Selected Elements, Izdat Nauka, 275–310, Moscow.
P a l i n k a { , L. A., Strmi }, S. & Herlec, U. 2001: The ore-forming fluids in the Idria mercury mine, Slovenia. In: A. Piestrzynski et al. (eds.), Mineral deposits at the beginning of the 21st cen-tury.– Swets & Zeitlinger, 321–324, Lisse.
P e t e r s , K. E. 1986: Guidelines for evaluating petroleum source rock using programmed pyrolysis. Am. Assoc. Petrol. Geol. Bull, 70, 318–329, Tulsa.
P l a c e r , L. 1982: Structural history of the Idrija mercury deposit. Geologija, 25, 7–94, Ljubljana.
P l a c e r , L. & ^ a r , J. 1977: The middle trias-sic structure of the Idrija region. – Geologija, 20, 141–166, Ljubljana.
S p a n g e n b e r g , J. E. & M a c k o , S. A. 1998: Organic geochemistry of the San Vicente zinc-lead district, eastern Pucara Basin, Peru.– Chem. Geol., 146, 1–23, Amsterdam.
S p a n g e n b e r g , J. E., H u n z i k e r , J., M e i s -s e r , N. & H e r l e c , U. 1999: Compound specific isotope analysis of the organic minerals hatchet-tite and idrialite in geodes, coal, and mineral deposits. In: C. J. Stanley et al. (eds.), Mineral Deposits: Processes to Processing.– Balkema, 275– 278, Rotterdam.
S t r u n z , H. & C o n t a g , B. 1965: Evenkit, Flagstaffit, Idrialin und Reficit.– Neues Jahrb. Mineral. Monatsh. 1, 19–25, Stuttgart.
W i s e , J. A., C a m p b e l l , R. M., W e s t , W. R., L e e , M. L. & B a r t l e , K. D. 1986: Characterization of polycyclic aromatic hydrocarbon minerals curtisite, idrialite, and pendeltonite using highperformance liquid chromatography, gas chroma-tography, mass spectrometry and nuclear resonance spectroscopy. – Chem. Geol., 54, 339–357, Amsterdam.