Refine
Year of publication
Document Type
- Doctoral Thesis (87) (remove)
Language
- English (87) (remove)
Has Fulltext
- yes (87)
Is part of the Bibliography
- no (87)
Keywords
- Boden (2)
- Caribbean (2)
- China (2)
- Roman (2)
- bank filtration (2)
- degradation (2)
- groundwater (2)
- 1,4-dioxane (1)
- 3D-Modellierung (1)
- Abbau (1)
Institute
- Geowissenschaften (87) (remove)
A graph theoretical approach to the analysis, comparison, and enumeration of crystal structures
(2008)
As an alternative approach to lattices and space groups, this work explores graph theory as a means to model crystal structures. The approach uses quotient graphs and nets - the graph theoretical equivalent of cells and lattices - to represent crystal structures. After a short review of related work, new classes of cycles in nets are introduced and their ability to distinguish between non-isomorphic nets and their computational complexity are evaluated. Then, two methods to estimate a structure’s density from the corresponding net are proposed. The first uses coordination sequences to estimate the number of nodes in a sphere, whereas the second method determines the maximal volume of a unit cell. Based on the quotient graph only, methods are proposed to determine whether nets consist of islands, chains, planes, or penetrating, disconnected sub-nets. An algorithm for the enumeration of crystal structures is revised and extended to a search for structures possessing certain properties. Particular attention is given to the exclusion of redundant nets and those, which, by the nature of their connectivity, cannot correspond to a crystal structure. Nets with four four-coordinated nodes, corresponding to sp3 hybridised carbon polymorphs with four atoms per unit cell, are completely enumerated in order to demonstrate the approach. In order to render quotient graphs and nets independent from crystal structures, they are reintroduced in a purely graph-theoretical way. Based on this, the issue of iso- and automorphism of nets is reexamined. It is shown that the topology of a net (that is the bonds in a crystal) constrains severely the symmetry of the embedding (that is the crystal), and in the case of connected nets the space group except for the setting. Several examples are studied and conclusions on phases are drawn (pseudo-cubic FeS2 versus pyrite; α- versus β- quartz; marcasite- versus rutile-like phases). As the automorphisms of certain quotient graphs stipulate a translational symmetry higher than an arbitrary embedding of the corresponding net would show, they are examined in more detail and a method to reduce the size of such quotient graphs is proposed. Besides two instructional examples with 2-dimensional graphs, the halite, calcite, magnesite, barytocalcite, and a strontium feldspar structures are discussed. For some of the structures it is shown that the quotient graph which is equivalent to a centred cell is reduced to a quotient graph equivalent to the primitive cell. For the partially disordered strontium feldspar, it is shown that even if it could be annealed to an ordered structure, the unit cell would likely remain unchanged. For the calcite and barytocalcite structures it is shown that the equivalent nets are not isomorphic.
In dieser Arbeit wurde der chemische Ozonverlust in der arktischen Stratosphäre über elf Jahre hinweg, zwischen 1991 und 2002, mit Hilfe der so genannten "Ozon-Tracer Korrelationstechnik" (TRAC), untersucht. Bei dieser Methode werden Korrelationen zwischen Ozon und langlebigen Spurenstoffen im Verlauf des Winters im Polarwirbels beobachtet und so der jährliche akkumulierte Ozonverlust berechnet. Die Ergebnisse dieser Arbeit basieren im wesentlichen auf Messdaten der Satelliteninstrumente: HALOE (Halogen Occultation Experiment) auf UARS (Upper Atmosphere Research Satellite) und ILAS (Improved Limb Atmospheric Spectrometer) Instrument auf ADEOS (Advanced Earth Observing Satellite). Das HALOE Instrument misst seit Oktober 1991 kontinuierlich alle zwei bis drei Monate für einige Tage in höheren nördlichen Breiten. ILAS lieferte ausschließlich für den Winter 1996-97 Messungen, die über sieben Monate hinweg in hohen Breiten aufgenommen wurden. Aufgrund der eingeführten Erweiterungen und Verbesserungen der Methode in dieser Arbeit, konnte die Methode anhand einer detaillierten Studie für den Winter 1996-97 validiert werden. Die ILAS Messreihe wurde dazu verwendet, erstmals die Untersuchung der zeitlichen Entwicklung von Ozon-Tracer Korrelationen kontinuierlich für die gesamte Lebensdauer des Polarwirbels durchzuführen. Dabei wurden auch Korrelationen während der Bildung des Wirbels untersucht und im Besonderen mögliche Mischungsvorgänge zwischen Wirbelluft und Luftmassen außerhalb des Wirbels. Ausserdem wurde ein Vergleich der Ergebnisse von ILAS und HALOE Messdaten durchgeführt und Unterschiede in den Ergebnissen tiefgreifend analysiert. Basierend auf HALOE Messungen konnte die erweiterte TRAC Methode über elf Jahren hinweg angewendet werden. Damit war erstmals eine konsistente Analyse von Ozonverlust und Chloraktivierung über diesen Zeitraum möglich. Die Erweiterungen führten zu einer Verringerung und genauen Quantifizierung von Unsicherheiten der Ergebnisse. Ein deutlicher Zusammenhang zwischen meteorologischen Bedingungen, Chloraktivierung und dem chemischen Ozonverlust wurde deutlich. Weiterhin zeigte sich eine Abhängigkeit zwischen den meteorologischen Bedingungen und der Homogenität des Ozonverlustes innerhalb eines Winters, sowie der mögliche Einfluss von horizontaler Mischung auf Luftmassen in einem schwach ausgeprägten Polarwirbel. In dieser Arbeit wurde eine positive Korrelation zwischen den über die gesamte Lebensdauer des Wirbels auftretenden möglichen PSC-Flächen und den akkumulierten Ozonverlusten für die elf untersuchten Jahre deutlich. Es konnte darüber hinaus gezeigt werden, dass der Ozonverlust von deutlich mehr Einflüssen als nur von der Fläche möglichen PSC Auftretens bestimmt wird, sondern zum Beispiel von der Stärke der Sonneneinstrahlung abhängt. Außerdem lassen sich Auswirkungen von Vulkanausbrüchen, wie zum Beispiel im Jahr 1991 der des Mount Pinatubo, identifizieren.
Derivation and characterization of a new filter for nonlinear high-dimensional data assimilation
(2015)
Data assimilation (DA) combines model forecasts with real-world observations to achieve an optimal estimate of the state of a dynamical system. The quality of predictions in nonlinear and chaotic systems such as atmospheric or oceanic circulation is strongly sensitive to the initial conditions. Therefore, beyond the consistent reconstruction of past states, a primary relevance of advanced DA methods concerns the proper model initialization. The ensemble Kalman filter (EnKF) and its deterministic variants, mostly square root filters such as the ensemble transform Kalman filter (ETKF), represent a popular alternative to variational DA schemes. They are applied in a wide range of research and operations. Their forecast step employs an ensemble integration that fully respects the nonlinear nature of the analyzed system. In the analysis step, they implicitly assume the prior state and observation errors to be Gaussian. Consequently, in nonlinear systems, the mean and covariance of the analysis ensemble are biased and these filters remain suboptimal. In contrast, the fully nonlinear, non-Gaussian particle filter (PF) relies on Bayes' theorem without further assumptions, which guarantees an exact asymptotic behavior. However, it is exposed to weight collapse, particularly in higher-dimensional settings, known as the curse of dimensionality.
This work presents a new method to obtain an analysis ensemble with mean and covariance that exactly match the corresponding Bayesian estimates. This is achieved by a deterministic matrix square root transformation of the forecast ensemble, and subsequently a suitable random rotation that significantly contributes to filter stability while preserving the required second-order statistics. The forecast step remains as in the ETKF. The algorithm, which is fairly easy to implement and computationally efficient, is referred to as the nonlinear ensemble transform filter (NETF). The limitation with respect to fully-nonlinear filtering is that the NETF only considers the mean and covariance of the Bayesian analysis density, neglecting higher-order moments.
The properties and performance of the proposed algorithm are investigated via a set of experiments. The results indicate that such a filter formulation can increase the analysis quality, even for relatively small ensemble sizes, compared to other ensemble filters in nonlinear, non-Gaussian scenarios. They also confirm that localization enhances the applicability of this PF-inspired scheme in larger-dimensional systems. Finally, the novel filter is coupled to a large-scale ocean general circulation model with a realistic observation scenario. The NETF remains stable with a small ensemble size and shows a consistent behavior. Additionally, its analyses exhibit low estimation errors, as revealed by a comparison with a free ensemble integration and the ETKF. The results confirm that, in principle, the filter can be applied successfully and as simple as the ETKF in high-dimensional problems. No further modifications are needed, even though the algorithm is only based on the particle weights. Thus, it is able to overcome the curse of dimensionality, even in deterministic systems. This proves that the NETF constitutes a promising and user-friendly method for nonlinear high-dimensional DA.
Over the last several decades, spinel-structured minerals with the chemical formula AB2O4 (where A and B stand for divalent and trivalent cations, respectively) have attracted more and more attention, particularly with regards to their breakdown at high pressures and temperatures and the nature of the so-called "post-spinel" phases. Spinel-structured phases with different endmember compositions, like magnetite (Fe3O4), hercynite (FeAl2O4) or spinel (MgAl2O4), are known to breakdown differently at high pressure-temperature conditions (e.g., Akaogi et al. 1999; Schollenbruch et al. 2010; Woodland et al. 2012). Such phases are of particular interest when they incorporate ferric (Fe3+) and ferrous (Fe2+) cations as this makes their stability sensitive to redox conditions. Since magnetite and magnesioferrite (MgFe3+ 2O4) have been found as inclusions in diamond (e.g., Stachel et al. 1998; Harte et al. 1999; Wirth et al. 2014; Palot et al. 2016; Jacob et al. 2016), understanding their phase relations is important for setting constraints on the conditions of their formation.
This study aimed to experimentally investigate the phase relations of Fe-Mg spinel-structured phases at conditions of the deep upper mantle and transition zone. Exploring the stability of new post-spinel phases and their characterization were also major goals of this study. Approaching a pyrolitic mantle composition by adding amounts of SiO2 in the system allowed constraints on the relevance of Fe-Mg post-spinel phases coexisting with mantle silicates to be made. ...
This thesis deals with the analysis of “presolar” silicates and oxides by high resolution mass spectrometry and electron microscopy techniques. This “stardust” was identified by its extreme oxygen isotopic anomalies, which point to nucleosynthetic reactions in stellar interiors, in the carbonaceous chondrite Acfer 094. Isotopic, chemical and mineralogical studies on these stardust grains therefore allow the testing of astrophysical questions on Earth, which are otherwise only accessible by spectroscopy and theoretical models. The class of presolar silicates has been identified only six years ago in 2002, although it was known already from spectroscopic observations that silicates represent the most abundant type of dust in the galaxy. The development of the “NanoSIMS” was a crucial step in this respect, because this ion probe with its superior spatial resolution of only 50 nm allowed the detection of the typically 300 nm sized presolar silicates. A total of 142 presolar silicates and 20 presolar oxides were identified within Acfer 094, whose matrix therefore contains 163 ± 14 ppm presolar silicates and 26 ± 6 ppm presolar oxides. This is among the highest amounts reported so far for any primitive solar system material. The majority of detected stardust grains derive from asymptotic giant branch stars of 1 – 2.5 Msun and close-to-solar or slightly lower-than-solar metallicity. However, by measuring the Si isotopic compositions of some enigmatic grains, it could be shown that there is a sub-class of presolar silicates characterized by an extreme enrichment of 17O and a moderate enhancement of 30Si relative to solar, whose origins might be explained by formation in binary stellar systems. About 10% of all grains exhibit an enrichment in 18O and some of them also of 28Si relative to solar, which most likely point to an origin in type II supernova explosions. The Si isotopic measurements also allowed to quantify the effect of the s-process on the Si isotopes in low-mass asymptotic giant branch stars. The results agree well with theoretical predictions. The grains were furthermore characterized by SEM and the chemistries of about half of the grains were determined by Auger electron spectroscopy. The majority of grain morphologies are consistent with what is expected from condensation experiments. However, a lot of grains are altered by Fe-rich minerals, which are either of primary condensation or of secondary ISM or solar nebula origin. Furthermore, complex presolar grains consisting of refractory Al-rich grains attached to silicate material could be identified, which have been predicted by condensation theory and observational evidence. Nine presolar silicates were analyzed by combined NanoSIMS/TEM studies. The majority of grains are Mg-rich and amorphous, which is in contrast to astrophysical evidence, which mainly postulate crystalline Mg-rich and amorphous Fe-rich circumstellar condensates. However, the grains might have been rendered amorphous by secondary processes in the ISM or could have condensed under non-equilibrium, low-temperature conditions in the circumstellar outflow. The grains are more likely characterized by a variable, pyroxene-like chemistry, which could be a result of sputtering in the ISM, which preferentially removes Mg. The detected crystalline presolar silicates in this study and in other work are all olivines, whereas grains with a pyroxene stoichiometry are all amorphous except one. This supports astrophysical models which point to different formation pathways for these two types of grains and therefore different crystallinity. However, the relatively high Fe content of three detected presolar olivines in this study and in other work is in contrast to astrophysical evidence and theoretical considerations, which predict essentially Fe-free crystalline grains. It is therefore possible that the infrared spectra might also be compatible with less Mg-rich olivines. The only crystalline presolar silicate with a pyroxene-like stoichiometry is the unusual grain 1_07: although it is chemically enstatite, the electron diffraction pattern could only be indexed to silicate perovskite, which is stable above ~23 GPa. The discovery of a high-pressure phase of presolar origin shows that dust grains encountering interstellar shocks might not necessarily be completely destroyed. In astrophysical models it is in principle also possible that a fraction of larger grains might survive such a shock wave encounter as a high-pressure modification, which is supported by this discovery.
During this study clumped isotope analysis of carbonates was established at the Goethe University of Frankfurt, Germany. Therefore, preparation protocols and analytical parameters were elaborated to obtain precise and accurate Δ47 data. Briefly, analyte CO2 was cleaned cryogenically using glass extraction lines to remove traces of water that enable re-equilibration of C–O bonds in the gases. Furthermore, analyte CO2 was passed through a gas chromatograph (GC) to clean it from contaminants that produce isobaric interferences with m/z 47. Initially, phosphoric acid digestions of carbonates was conducted at 25 °C in McCrea-type reaction vessels. Afterwards samples were reacted at 90 °C using a common acid bath. Mass spectrometric analyses were performed using a MAT 253 equipped with a dual inlet system. Δ47 values were directly projected to the absolute scale using CO2 gases equilibrated at distinct temperatures.
In cooperation with Stefano Bernasconi and his research group at ETH Zurich we studied the non-linearity that occurs for the measurement of m/z 47. This effect results from secondary electrons created by the m/z 44 beam. These electrons cause a negative background on the m/z 47 collector. A correction procedure was proposed that relies on the determination of the negative background on the m/z 47 Faraday cup. This approach might reduce time-consuming analyses of heated gases which were used so far to account for the observed non-linearity. However, the suggested correction of the negative background on the m/z 47 cup is only applicable if the slit width of the m/z 44 beam is significantly wider than that of the m/z 47 beam.
This thesis, furthermore, presents a comparison of the different phosphoric acid digestion techniques which are commonly used for carbonate clumped isotope analysis. For calcitic and aragonitic material digested at 25 °C in McCrea-type vessels we observed that the sample size has an effect on Δ47 data: higher mean Δ47 values and a larger scatter of data were received for samples <7 mg than for larger aliquots. For carbonate samples digested at 90 °C in a common acid bath no sample size effect was determined. We assume that secondary re-equilibration of CO2 with water preferentially occurs at 25 °C producing the observed differences. However, a sample size effect can be avoided if reaction temperature is increased to 90 °C.
In order to make carbonate Δ47 data obtained from acid digestions at 90 °C comparable to Δ47 data received from reactions at 25 °C the difference of the acid fractionation factores (Δ47*25-90) between both temperatures has to be known. For the determination of the Δ47*25-90 value we have considered Δ47 data made at 25 °C from samples >7 mg only. For calicte and aragonite we obtained differences in fractionation factores of 0.075‰ and 0.066‰, respectively. These Δ47*25-90 values are coincident with the theoretical prediction of 0.069‰ proposed for calcite (Guo et al., 2009).
Moreover, this dissertation comprises a calibration study of the clumped isotope thermometer based on various natural calcites that grew between 9 and 38 °C. The samples include a brachiopod shell, a bivalve shell, an eggshell of an ostrich and foraminifera tests which formed from distinct biomineralizing processes. Furthermore we included an authigenic carbonate crystallized from biological-induced precipitation. The following linear relationship between 1/T2 and Δ47 was determined (with Δ47 in ‰ and T in K):
Δ47 = 0.0327 (± 0.0026) x 106 / T2 + 0.3030 (± 0.0308) (R2 = 0.9915)
This equation differs from the pioneering Ghosh et al. (2006a) calibration. However, our regression line is statistically indistinguishable from that of Henkes et al. (2013) which is based on aragonitic mollusks and calcitic brachiopod shells. Both studies have in common that calibration data were, at first, directly referenced to the absolute scale. In addition, both datasets rely on similar digestion techniques. Furthermore, the two calibrations are conform with the theoretical prediction of Guo et al. (2009).
The calcite calibration of the clumped isotope paleothermometer received in this study was applied to Δ47 data measured for Silurian brachiopods shells from Gotland/Sweden. Prior to isotopic analysis the fossils were intensively investigated for their preservation state (CL, SEM, trace elements). The lowest T(Δ47) values of ca. 28 to 33 °C were estimated from ultrastructurally well-preserved regions of some shells. For these samples also the lowest δ18Ow values of Silurian seawater were determined. These estimates of ca. −1‰ confirm the assumption that the δ18O value of the Silurian ocean was buffered to (0 ± 1)‰.
Nevertheless, most studied shells were characterized by a patchwork of pristine and altered shell portions resulting in elevated T(Δ47) values which plot mostly between 40 and 60 °C. Our results indicate that the clumped isotopic composition of the shells were altered at low water-rock ratios, not affecting the δ18O values. Δ47 and δ18O data of associated diagenetic phases (sparitic and micritic phases of the inner fillings of the fossils) provide evidence that the sparitic cements grew during several diagenetic events which occurred at different temperatures in fluid-buffered systems. We, furthermore, conclude that the micritic phases lithified at a very early diagenetic stage with the δ18O values being most probably close to a Silurian seawater composition
Atmospheric nanoaerosols have extensive effects on the Earth’s climate and human health. This cumulative work focuses on the development and characterization of instrumentation for measuring various parameters of atmospheric nanoaerosols, and its use to understand new particle formation from organic precursors. The principal research question is, how the chemical composition of nanoaerosol particles can be measured and how atmospheric chemistry influences aerosol processes, especially new particle formation and growth. Therefore, nanoaerosols are investigated under various aspects. More specifically, an instrument is developed to analyze nanoparticles, and field as well as chamber studies are conducted.
The main project is the instrument development of the Thermal Desorption Differential Mobility Analyzer (TD-DMA, project 1, Wagner et al. (2018)). This instrument analyzes the chemical composition of small aerosol particles. By characterization and testing in chamber experiments, it is proven to be suitable for the analysis of freshly nucleated particles.
The second project (Wagner et al. (2017)) applies a broad spectrum of aerosol measurement instruments for the characterization of aerosol particles produced by a skyscraper blasting. A comprehensive picture of the particle population emitted by the demolition is obtained.
Project 3 (K¨urten et al. (2016)) is also an ambient aerosol measurement, focusing of new particle formation in a rural area in central Germany, and the ability of a negative nitrate CI-APi-TOF to detect various substances in atmosphere. Project 4 (Heinritzi et al. (2016)) is a characterization of the negative nitrate CI-APi-TOF used in projects 1, 3, 5, 6, 7 and 8. The following projects focus on understanding new particle formation from atmospherically abundant organic precursors. Key instruments comprise the negative nitrate CI-APiTOF for gas-phase measurements of the nucleating species, and various sizing and counting instruments for quantifying the particle formation and growth. Project 5 (Kirkby et al. (2016)) shows that biogenic organic compounds formed from alpha-pinene can nucleate on their own without the influence of e.g. sulfuric acid. Project 6 (Tr¨ostl et al. (2016)) describe the subsequent growth of these particles. Project 7 (Stolzenburg et al. (2018)) covers the temperature dependence of this growth and in project 8 (Heinritzi et al. (2018)), the suppressing influence of isoprene on the new particle formation is assessed.
Since the discovery of the ozone hole [Farman et al., 1985], the dynamics of the stratosphere and the transport of anthropogenic trace gases from the surface to the higher atmosphere has come into the focus of interest. In the tropics, air rises high into the stratosphere and is transported poleward. Trace gases like the CFCs1, photochemically stable in the troposphere, are thus transported into regions where they are photolyzed. The products of the photolysis reactions (eg. Cl, Br) interact in the catalytic ozone cycles and lead to enhanced ozone depletion. Regarding the transport of trace gases, the so-called lowermost stratosphere (LMS) is a very interesting region, where the troposphere and the stratosphere directly interact and air masses out of both regions are mixed. It is the lowest part of the stratosphere between the tropopause and the 380 K isentrope. Tropospheric air can enter this region directly via isentropic transport across the extra-tropical tropopause whereas stratospheric air descends across the 380 K isentrope via the mean meridional circulation. Stratosphere-troposphere exchange (STE) controls the chemical composition of the LMS as well as of the tropopause region and thus has an important effect on the radiative and chemical balance of these regions and the climate system. STE exhibits a strong seasonality [Holton et al., 1995]. While downwelling of stratospheric air across the 380 K isentrope is the dominant process in winter, troposphere-to-stratosphere transport (TST) gains importance in summer, when the downwelling from the stratosphere is only weak. Isentropic transport across the extra-tropical tropopause occurs in regions where the tropopause is strongly disturbed and is connected to tropopause folds, streamer events, frontal zones, polar and subtropical jets, warm conveyor belts and cut-off low systems. A short introduction into STE, the LMS region, and methods to study atmospheric transport is given in Chapter 1. One useful tool to analyse the motions of air and transport processes are longlived trace gases. Since the lifetimes of these tracers are longer than the time scale of the transport processes they are involved in, the distribution of tracers in the atmosphere is mostly determined by dynamics. In the context of this thesis, measurements of such long-lived tracers were performed and used to study transport into the LMS region in the northern hemisphere. During the Vintersol/EuPLEx and ENVISAT validation campaigns in winter 2003, long-lived tracers such as N2O, CH4, CFC-12, CFC-11, H-1211, H2, SF6 and CO2 were measured with the High Altitude Gas Analyser (HAGAR), a two channel in-situ gas chromatograph combined with a CO2 instrument, based on nondispersive infrared absorption. Combined with measurements taken during campaigns in Forli/Italy (ENVISAT validation) in July and October 2002, tracer data were gathered from the tropopause up to altitudes around 20 km during 25 flights on board the Russian high-altitude aircraft M55 Geophysica. Thus, a substantial set of high quality tracer data has been obtained covering the polar vortex region as well as the mid latitudes of the northern hemisphere. Chapter 2 gives an overview of the HAGAR instrument and necessary improvements of the instrumental set up (implementing a CH4 channel) that were performed in the context of this thesis, and review data processing, the measurement campaigns. In order to study transport into the LMS it is assumed that air basically enters the LMS via three different pathways: a) quasi-isentropic transport from the troposphere, b) downward advection from the middle stratosphere through the 380 K surface and c) in the polar vortex region subsidence of air from of the polar vortex. Fractions of air originating in each of these source regions are determined with a simple mass balance calculation by using observations of a subset of the above species with distinct lifetimes (N2O, CH4, CFC-11, H-1211, H2 and O3) yielding complementary constraints on transport from each region. Details of the mass balance calculation and the results are presented in Chapter 3. During the mid-latitude measurement campaigns in Forlí the passing of a cut-off low system associated with an elongated streamer over Europe was observed. The impacts of this event on the trace gas mixing ratios in the LMS are examined in Chapter 4. Finally, a summary is given in Chapter 5.
This PhD thesis has been carried out within an interdisciplinary cooperational project between the Deutsches Bergbau-Museum Bochum and the Goethe-Universität Frankfurt, which is dedicated to ancient Pb-Ag mining and metal production in the hinterland of the municipium Ulpiana in central Kosovo. Geochemical analysis (OM, XRD, EMP, MC-ICP-MS) of ores, metallurgical (by-) products and metal artefacts allowed to reconstruct the local chaîne opératoire and to decipher significant chronological differences between presumably Roman/late antique and medieval/early modern metallurgical processing. Pb isotope provenance studies documented the relevance of local metal production within the Roman Empire and confirmed the actual existence of a Metalla Dardanica district, which until now solely has been suspected on basis of epigraphy.
The predominant abundance of the by-products matte (Cu, Pb, Fe and Zn sulphides) and speiss (ferrous speiss: Fe-As compounds; base metal speiss: ~(Cu,Ni,Fe,Ag )x(Sb,Sn,As )y ) at smelting sites with a preliminary Roman/late antique dating points to treatment of complex polymetallic ore. Pb isotope analysis demonstrated that the mining district of Shashkoc-Janjevo (partially) supplied six of the ten investigated metallurgical sites. In this mineralisation, parageneses with elevated Cu, As and Sb abundances comprise significant proportions of particularly tennantite-tetrahedrite minerals, chalcopyrite, arsenopyrite and were generated during the early and main stages of ore formation. Later precipitated ore in contrast is marked by a significantly less versatile mineralogy and consists almost exclusively of galena, sphalerite and pyrite/marcasite. Besides increased Cu, As and Sb contents, ore from the main formation stage also exhibits generally higher Ag abundances, which are mainly hosted by fahlore and locally abundant secondary Cu sulphides (chalcocite, digenite and covellite) and oxidised phases (e.g. malachite, azurite). The higher precious metal grades of this ore type, whose geochemical signature (i.e. higher proportions of Cu, As and Sb) is mirrored by the abundance of the metallurgical by-products matte and speiss (almost exclusively found at potentially Roman/late antique smelting sites; see above), presumably were a pivotal factor leading to its preferential exploitation in earlier times. Matte and base metal-rich speiss contain notable amounts of Ag, which are mainly present in Cu-(Fe) sulphides and particularly antimonides ((Cu,Ni)2Sb, Ag3Sb), respectively. While the speiss compounds due to their close association with Pb bullion presumably were cupelled automatically, the metallurgical treatment of matte could not have been proven unambiguously, but overall certainly is highly likely.
The beneficiated ore (i.e. crushed and sorted, potentially also treated by more lavish techniques such as grinding, sieving or wet-mechanical methods) possibly was partially roasted and subsequently together with fluxes and charcoal submitted to the furnaces. The working temperatures approximately ranged between 1100 and 1400 °C. Slags from all presumably Roman/late antique dated and few of their potentially medieval/early modern analogues were produced from smelting of (partially roasted) ore with charcoal and added siliceous material, thus resulting in fayalite-dominant phase assemblages or rarely observed glassy parageneses. Even though several subtypes of fayalite slags have been established on basis of the abundance of Fe-rich oxide phases (i.e. spinel ss and wüstite), late clinopyroxene and the general solidification sequence of the slags, the process conditions (i.e. temperature, fO2, added fluxing agents) must have been widely similar; chemical variations could be explained by varying degrees of interaction of the slag melt with charcoal ash and furnace material. The other investigated metallurgical remains indicate employment of a calcareous flux, which led to formation of Ca-rich olivine-, olivine+clinopyroxene-, clinopyroxene- or melilite-type slags. These types as well as glassy slags were generated at more oxidising conditions outside the fayalite stability field (FMQ buffer equilibrium, cf. Lindsley, 1976) than their olivine-dominant analogues. Conclusions on the furnace construction could be drawn on basis of the typology of the slags, which mostly were tapped into a basin located outside the furnace, but partially (at two presumably medieval/early modern sites) also accumulated in a reservoir within the smelter.
Lead artefacts excavated in Ulpiana could be isotopically related to ores from mineralisations in its vicinity and demonstrate that the resources were at least utilised for local metal production. However, also ship wreck cargo from Israel - including several lead ingots with the inscription 'MET DARD' (Raban, 1999) - and late antique lead-glazed pottery from Serbia and Romania (Walton & Tite, 2010) could be related to a possible Kosovarian/Serbian provenance of the raw material and thus indicate flourishing trade of metal from the Metalla Dardanica district within the Roman Empire.
References:
Lindsley, D. H. (1976). Experimental studies of oxide minerals. In D. Rumble, III (Hrsg.), Oxide minerals (61-88). Reviews in Mineralogy, Volume 3. Washington, DC: Mineralogical Society of America.
Raban, A. (1999). The lead ingots from the wreck site (area K8). Journal of Roman Archaeology, Supplementary Series, 35, 179-188.
Walton, M. S., & Tite, M. S. (2010). Production technology of Roman lead-glazed pottery and its continuance into late antiquity. Archaeometry, 52(5), 733-759.
Forty two samples of the Late Eocene Kiliran oil shale, Central Sumatra Basin, Indonesia were collected from a 102 m long drill core. The oil shale core represents the deposition time of about 240.000 years. Palynofacies and geochemical analyses have been carried out to reconstruct the paleoenvironmental conditions and paleoecology during deposition of the oil shale. Amorphous organic matter (AOM) is very abundant (>76%). B. braunii palynomorphs are present (3-16%) as the only autochtonous structured organic matter and generally more abundant in the middle part of the profile. The stable carbon isotopic composition of bulk organic matter (13C) varies from -27.0 to -30.5‰ and is generally more depleted in the middle part of the profile. The ratio of total organic carbon to sulfur (TOC/S), used as salinity indicator, ranges from 2.5 to 15.8 and shows variations along the profile. Slightly less saline environments are observed in the middle part of the profile. Fungal remains are generally present only in this part with a distinct peak of abundance. The presence of fungal remains is regarded as an indication for a relatively warmer climate during deposition of the middle part of the profile. The warmer climate is thought to influence the establishment of a thermocline, limiting the supply of recycled nutrients to epilimnion. Consequently, the primary productivity in the Kiliran lake decreased during deposition of the middle part of the profile as indicated by the relatively depleted 13C values and the blooming of B. braunii. The chemocline was also shoaling during the deposition according to the higher abundance of total isorenieratane and its derivatives originated from green sulfur bacteria dwelling in the photic zone euxinia. The warmer climate is also thought to influence the slightly decrease of water salinity during deposition of the middle part of the profile. The occurrence of B. braunii in Kiliran lake is also recognized from organic geochemical data. The distribution of n-alkanes is characterized by the unusual high amount of C27 n-alkane relative to the other long-chain n-alkanes. The concentrations of C27 n-alkane vary from 30.1 to 393.7 μg/g TOC and are generally in parallel with the abundances of B. braunii palynomorphs along the profile. The 13C values of this compound are about -31‰ and up to 2‰ enriched relative to those of the adjacent long-chain n-alkanes. B. braunii race A can thus be regarded as the significant biological source of the C27 n-alkane. Lower amounts of lycopane are observed in many oil shale samples (0 to 54.7 μg/g TOC). The 13C value of this compound is 17.2‰. This strong enrichment of 13C suggests that the lycopane was derived from B. braunii race L. The concentrations of lycopane develop generally in opposite with those of C27 μalkane. It is likely that both B. braunii races bloomed in alternation in the lake, probably due to changes on specific water chemistry. Norneohop-13(18)-ene and neohop-13(18)-ene derived from methanotrophic bacteria are the dominant hopanoid hydrocarbons. The sum of their concentrations varies from 40.6 to 360.0 μg/g TOC. The 13C of these compounds are extremely depleted (-45.2 to -50.2‰). The occurrence of abundant bacteria including methanotrophic bacteria was responsible for the recycling of carbon below the chemocline of the lake. The effect of the recycling of carbon is observed by the presence of a concomitant depletion (about 7-9‰) in 13C of some specific biomarkers derived from organisms dwelling in the whole phototrophic zone. 4-Methylsterane and 4-methyldiasterene homologues occur in the oil shale as the predominant biomarkers. The sum of the concentrations of all homologues are about 40.3-1,009.2 μg/g TOC with generally higher values in the uppermost and lower parts of the profile. Calcium (Ca) accounts as the predominant element in the oil shale, ranging from 5.0 to 16.7%. This element shows generally parallel variation with the 4-methylsterane and 4-methyldiasterene homologues along the profile. This suggests that these compounds were derived from biological sources favoring more alkaline and more trophic environments. On the other hand, these compounds were less abundant in the middle part of the profile which is consistent with less alkaline and less trophic environments promoting B. braunii to bloom. Alternation between Dinoflagellates and B. braunii in ancient lacustrine environments due to water chemistry changes have been known from previous studies. In the present case, distinct alternation between B. braunii abundances and concentrations of 4-methylsterane and 4-methyldiasterene homologues along the studied oil shale profile suggest a hypothesis that these compounds were derived from freshwater Dinoflagellates although dinosterane is not present in the sediment extracts. Water alkalinity and trophic level changes were most likely responsible for the alternation of Dinoflagellates and B. braunii blooming.