Ophiolite remnants at the eastern margin of the Bohemian Massif and their bearing on the tectonic evolution

Summary Amphibolites are widespread in the eastern part of the Bohemian Massif, east of the South Bohemian Pluton. Based on their geological situation, their metamorphic evolution and their geochemistry, they were separated into three genetically different units: the Rehberg ophiolite, the Buschandlwand amphibolite and the Raabs group. The metamorphic Rehberg ophiolite consists of a ultramafic to mafic plutonic sequence overlain by a gabbro/dike complex and a volcanic section with basaltic, andesitic and rhyolitic volcanics associated with pelitic to psammitic sediments. The entire ophiolite underwent amphibolite facies metamorphism. Rock/MORB normalisation plots and other element ratio plots, such as Ti vs. V or Ta/Yb vs. Th/Yb argue for a supra-subduction zone environment as site of the origin of the Rehberg ophiolite. The Letovice ophiolite in Moravia is structured in a similar way. Geochemically it shows a more MORB like composition but also a distinct tholeiitic island arc group. Although there is no age data for the protolith of both ophiolites, it is believed that they have formed in an oceanic basin separating the Brunovistulian-Moravian block to the east and the Moldanubian block to the west. Possible traces of the continuation of the ophiolites towards the north can be found in the Stare Mesto ophiolites and the ophiolites surrounding the Sowie Gory block in the Sudetes.

---

Ophiolitrelikte am ostrand der Böhmischen Masse und ihre bedeutung für die tektonische entwicklung

Zusammenfassung Amphibolite sind in der Böhmischen Masse, östlich des Südböhmischen Plutons weit verbreitet. Ausgehend von ihrer geologischen Position, ihrer metamorphen Entwick lung und ihrer Geochemie lassen sie sich in drei genetisch unterschiedliche Einheiten teilen: den Rehberg Ophiolith, den Buschandlwand Amphibolit und die Raabser Gruppe. Der metamorphe Rehberg Ophiolith besteht aus einer ultramatischen bis mafischen plutonischen Abfolge, die von einem Gabbro/Gangkomplex überlagert wird und einer vulkanischen Folge mit Basalten, Andesiten and Rhyoliten, die mit tonigen und sandigen Sedimenten verknüpft sind. Der gesamte Ophiolithkörper wurde unter amphibolitfaziellen Bedingungen metamorph überprägt. Rock/MORB Diagramme und andere Elementdarstellungen, wie z.B. das Ti vs. V, oder das Ta/Yb vs. Th/Yb Diagramm sprechen für eine Entstehung des Rehberg Ophioliths in einer Suprasubduction Zone. Der Letovice Ophiolith in Mähren weist eine ähnliche Struktur auf. Geochemisch zeigt er stärkere Ähnlichkeiten mit MORB, enthält aber auch Elemente von tholeiitischen Inselbogenbasalten. Auch wenn es für beide Ophiolithe noch keine Protolith-Altersdaten gibt, scheint es dennoch wahrscheinlich, daß sie in einem Ozeanbecken, das den Brunovistulisch-Moravischen Kontinentalblock im Osten vom Moldanubischen Kontinent im Westen trennt, gebildet wurden. Die mögliche Fortsetzung beider Ophiolithe könnten im Norden der Stare Mesto Ophiolith und die Ophiolithe rund um den Sowie Gory Block in den Sudeten darstellen.

---

---

With 7 Figures     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

探测冰盖和多年冻土层的融冰探头研发

The exploration of ice sheets by melting vertical holes into the ground has some tradition in terrestrial glaciology. Such probes have been used since the 1960's to investigate the vertical structure of the ice in Greenland and Antarctica and in alpine glaciers. In this paper we look into the possibility to develop similar devices for use on extraterrestrial icy bodies, like e.g. the polar areas on Mars or the icy satellites of the outer solar system. We report on some basic experiments performed in the cryo-vacuum laboratory of the Space Research Institute of the Austrian Academy of Sciences, Graz. In these experiments the penetration of a simple melting probe into compact and porous water ice (with a snow-like texture) was monitored, both under vacuum conditions and under air pressure. The observed penetration speeds for a given power supply are compared with a simple mathematical model. We conclude that a miniature melting probe with small overall dimensions and a reasonable power demand could well be part of the payload of a future planetary mission, for example to the poles of Mars. Such missions are currently under discussion in several space agencies. Moreover such probes could also e? ectively be used in terrestrial environments. A possible design is presented at the end of the paper.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

Tetrahedral Plots of the Phase Relations for Basalts

The phase relations of quaternary systems are generally represented by projections onto ternary compositional planes. Such projections often obscure relationships that would only be evident in a three-dimensional tetrahedral plot. The tetrahedral plot requires that compositions of the minerals and melts be transformed into Cartesian coordinates. It is shown here how this transformation is carried out. The application is demonstrated by tetrahedral plots of experimental melt compositions of partially molten lherzolite. Furthermore, the plot can be used to evaluate whether or not a particular basaltic composition represents a primary melt. The methods are applicable to any four-component system.     

Spatially controlled Fe and Si isotope variations: an alternative view on the formation of the Torres del Paine pluton

We present new Fe and Si isotope ratio data for the Torres del Paine igneous complex in southern Chile. The multi-composition pluton consists of an approximately 1 km vertical exposure of homogenous granite overlying a contemporaneous 250-m-thick mafic gabbro suite. This first-of-its-kind spatially dependent Fe and Si isotope investigation of a convergent margin-related pluton aims to understand the nature of granite and silicic igneous rock formation. Results collected by MC-ICP-MS show a trend of increasing δ⁵⁶Fe and δ³⁰Si with increasing silica content as well as a systematic increase in δ⁵⁶Fe away from the mafic base of the pluton. The marginal Torres del Paine granites have heavier Fe isotope signatures (δ⁵⁶Fe = +0.25 ± 0.02 2se) compared to granites found in the interior pluton (δ⁵⁶Fe = +0.17 ± 0.02 2se). Cerro Toro country rock values are isotopically light in both Fe and Si isotopic systems (δ⁵⁶Fe = +0.05 ± 0.02 ‰; δ³⁰Si = ?0.38 ± 0.07 ‰). The variations in the Fe and Si isotopic data cannot be accounted for by local assimilation of the wall rocks, in situ fractional crystallization, late-stage fluid exsolution or some combination of these processes. Instead, we conclude that thermal diffusion or source magma variation is the most likely process producing Fe isotope ratio variations in the Torres del Paine pluton.