The Othris Ophiolite, Greece: A snapshot of subduction initiation at a mid-ocean ridge

The mantle section of the Tethyan-type Othris Ophiolite, Greece, records tectono-magmatic processes characteristic of both mid-ocean ridges and supra-subduction zones. The Othris Ophiolite is a remnant of the Jurassic Neotethys Ocean, which existed between Eurasia and Gondwanaland. Othris peridotites range from fertile plagioclase lherzolites to depleted harzburgites. Abundances of Al₂O₃ and CaO show well-defined inverse linear correlations with MgO, suggesting that the Othris peridotites formed as residua from variable degrees of partial melting.

Peridotites from the Fournos Kaïtsa and western Katáchloron sub-massifs are similar to abyssal peridotites and can be explained by a multistage model with some melting in the garnet stability field followed by moderate degrees of anhydrous near-fractional melting in the spinel stability field. In contrast, the peridotites from the Metalleio, Eretria, and eastern Katáchloron sub-massifs, and the Vourinos ophiolite are highly depleted and have extremely low concentrations of Al₂O₃ and heavy rare earth elements. These peridotites have enriched light REE contents compared to the middle REE. These residua are best modelled by hydrous melting due to a flux of slab-derived fluid to the mantle wedge during melting.

The occurrence of both styles of melting regimes within close spatial and temporal association in the same ophiolite is explained by intra-oceanic thrusting and forced subduction initiation at (or near) a mid-ocean ridge. Thus, the Othris Ophiolite, and probably Tethyan-type ophiolites in general, represent a transient phase of plate tectonic reorganisation rather than quasi-steady state plate tectonics.     

Garnet growth at high- and ultra-high pressure conditions and the effect of element fractionation on mineral modes and composition

In this paper we show that thermodynamic forward modelling, using Gibbs energy minimisation with consideration of element fractionation into refractory phases and/or liberated fluids, is able to extract information about the complex physical and chemical evolution of a deeply subducted rock volume. By comparing complex compositional growth zonations in garnets from high-and ultra-high pressure samples with those derived from thermodynamic forward modelling, we yield an insight into the effects of element fractionation on composition and modes of the co-genetic metamorphic phase assemblage. Our results demonstrate that fractionation effects cause discontinuous growth and re-crystallisation of metamorphic minerals in high pressure rocks. Reduced or hindered mineral growth at UHP conditions can control the inclusion and preservation of minerals indicative for UHP metamorphism, such as coesite, thus masking peak pressure conditions reached in subducted rocks.Further, our results demonstrate that fractional garnet crystallisation leads to strong compositional gradients and step-like zonation patterns in garnet, a feature often observed in high-and ultra-high pressure rocks. Thermodynamic forward modelling allows the interpretation of commonly observed garnet growth zonation patterns in terms of garnet forming reactions and the relative timing of garnet growth with respect to the rock's pressure–temperature path. Such a correlation is essential for the determination of tectonic and metamorphic rates in subduction zones as well as for the understanding of trace element signatures in subduction related rocks. It therefore should be commonplace in the investigation of metamorphic processes in subduction zones.     

Evaluation of Probe-Induced Flow Distortion of Campbell CSAT3 Sonic Anemometers by Numerical Simulation

The Campbell CSAT3 sonic anemometer is one of the most popular instruments for turbulence measurements in basic micrometeorological research and ecological applications. While measurement uncertainty has been characterized by field experiments and wind-tunnel studies in the past, there are conflicting estimates, which motivated us to conduct a numerical experiment using large-eddy simulation to evaluate the probe-induced flow distortion of the CSAT3 anemometer under controlled conditions, and with exact knowledge of the undisturbed flow. As opposed to wind-tunnel studies, we imposed oscillations in both the vertical and horizontal velocity components at the distinct frequencies and amplitudes found in typical turbulence spectra in the surface layer. The resulting flow-distortion errors for the standard deviations of the vertical velocity component range from 3 to 7%, and from 1 to 3% for the horizontal velocity component, depending on the azimuth angle. The magnitude of these errors is almost independent of the frequency of wind speed fluctuations, provided the amplitude is typical for surface-layer turbulence. A comparison of the corrections for transducer shadowing proposed by both Kaimal et al. (Proc Dyn Flow Conf, 551–565, 1978) and Horst et al. (Boundary-Layer Meteorol 155:371–395, 2015) show that both methods compensate for a larger part of the observed error, but do not sufficiently account for the azimuth dependency. Further numerical simulations could be conducted in the future to characterize the flow distortion induced by other existing types of sonic anemometers for the purposes of optimizing their geometry.     

Force-directed layout of origin-destination flow maps

This paper introduces a force-directed layout method for creating origin-destination flow maps. Design principles derived from manual cartography and automated graph drawing to increase readability of flow maps and graph layouts are taken into account. The origin-destination flow maps produced with our algorithm show flows with quadratic Bézier curves that reduce flow-on-flow and flow-on-node overlaps, and avoid sharp or irregular bends in flow lines. A survey of expert cartographers found that flow maps created with our automated method are similar in quality to manually produced flow maps.     

Phenomenon-based Specification of the Digital Representation of Terrain Surfaces

The phenomenon known as 'terrain' is a continuous surface. However, when a digital terrain representation is based on a regular raster (i.e. a DEM) the digital surface is commonly not continuous. This is the case for the derivation of variables such as slope, aspect, and curvature values as performed in today's Geographic Information Systems (GIS). Often, there is no surface specified at all, as, for instance, when flow lines or watersheds are constructed. The discrepancy between the phenomenon to be modelled and its digital representation causes the terrain analysis results to be less accurate than they could be. Furthermore, if more than one type of terrain information is derived the results are likely to be based on different specifications of the seemingly same terrain surface. The combined application of the derivation results will likely introduce inconsistencies. This paper suggests founding the specification of digital terrain representations on a careful analysis of the properties of the phenomenon. The paper details the reasons for, and advantages of, continuous surface representations and emphasises the importance of a comprehensive documentation of the conceptual models underlying digital terrain representations. A review of suitable interpolation approaches for the specification of terrain surfaces is given. The paper discusses how the resulting digital surfaces are analysed and how measurement uncertainty may be accounted for.     

The maximum Ice Age (LGM) glaciation of the Central- and South Karakorum: an investigation of the heights of its glacier levels and ice thicknesses as well as lowest prehistoric ice margin positions in the Hindukush, Himalaya and in East-Tibet on the Minya Konka-massif

A continuing prehistoric ice stream network between the Karakorum main crest and the Nanga Parbat massive has been evidenced, which, flowing down from the current Baltoro- and Chogolungma glaciers and filling the Shigar valley as well as the Skardu Basin, has flowed together with the Gilgit valley glacier to a joint Indus parent glacier through the Indus gorge. The ice stream network received an influx by a plateau glacier covering the Deosai plateau, which was connected through outlet glaciers to the ice filling of the Skardu Basin and the Astor glacier at the Nanga Parbat, as well as to the lower Indus glacier. The field observations introduced here in part confirm the results as to the Ice Age glacier surface area of Lydekker, Oestreich and Dainelli, but go beyond it. In additon, a reconstruction of the surface level of this ice stream network and its glacier thicknesses up to the highest regions of the present-day Karakorum valley glaciers has been carried out for the first time. In the area under investigation the Karakorum ice stream network showed three ice cupolas, culminating at an altitude of 6200–6400 m. Between the mountain groups towering 1000–2000 m higher up, they communicated with each other over the transfluence passes in a continuous glacier surface without breaks in slope. In the Braldu- and Basna valley ice thicknesses of 2400–2900 m have been reached. In the Skardu Basin, where the glacier thickness had decreased to c. 1500–1000 m, the ELA at an ice level of 3500-3200 m asl had fallen short to the extent that from here on down the Indus glacier a surface moraine cover has to be suggested. However, 80% of the surface of the ice stream network was devoid of debris and had an albedo of 75-90%. The lowest joint glacier terminus of the ice stream network was situated - as has already been published in 1988 – in the lower Indus valley at 850–800 m asl. The reconstructed maximum extension of the ice stream network has been classified as belonging to the LGM in the wider sense (60–18 Ka BP). Four Late Glacial glacier positions (I–IV), with a decreasing ice filling of the valleys, have been differentiated, which can be locally recognized through polish lines and lateral moraine ledges. The valley (trough-) flanks with their ground moraine covers, oversteepened by glacier abrasion, have been gravitationally destroyed by crumblings, slides and rock avalanches since the deglaciation, so that an interglacial fluvial-, i.e. V-shaped valley relief has been developed from the in part preserved glacial relief. The contrast of the current morphodynamics with regard to the preserved forms is seen as an indication of the prehistorically completely different - namely glacigenic – valley development and the obvious rapidity of this reshaping at still clearly preserved glacial forms provides evidence of their LGM-age. In an additional chapter the lowest ice margin positions, so far unpublished, are introduced, which have been reconstructed for the Hindukush, Central Himalaya and on the eastern margin of Tibet.     

Age,origin and geodynamic significance of a polymetamorphic felsic intrusion in the Ötztal Crystalline Basement,Tirol, Austria

Summary A multi-method approach was applied to derive the age and origin of an orthogneiss body located in the central Kaunertal, western Ötztal Crystalline Basement (ÖCB). The Tieftal orthogneiss body is an internally differentiated, polymetamorphosed epizonal intrusion, embedded in amphibolites. It comprises leucocratic hedenbergite-hornblende-, hornblende- and biotite-hornblende-gneisses, but also some melanocratic rock types. The leucocratic Tieftal gneisses are granitic, have a near eutectic melt composition and share some features of A-type granites, such as high Na₂O+K₂O(8.07 to 8.58wt%), Zr (379 to 554ppm) and Y (58 to 79ppm) contents. The REE-patterns are rather flat ((La/Yb)_N=2.4 to 3.7), with distinct negative Eu anomalies. Single zircon evaporation dating of two samples and Sm-Nd dating of relict magmatic titanite resulted in ages of 487±7, 484±3 and 487±5Ma, respectively. The weighted mean of 485±3Ma is interpreted as the primary crystallization age of the Tieftal orthogneiss body. Rb-Sr whole rock dating results in a well defined regression line, corresponding to an age of 411±9Ma. This age clearly documents at least a partial resetting of the whole rock Rb-Sr system, which is most probably due to subsequent metamorphic overprint. The leucocratic Tieftal gneisses are isotopically rather primitive with an Nd _CHUR ^{485 Ma} value of +1.7 and a calculated magmatic initial⁸⁷Sr/⁸⁶Sr ratio of 0.7047. These data suggest a major mantle contribution. Most probably, they originated through fractionation of the magmatic precursors of the accompanying tholeiitic metabasites. The more primitive isotopic composition of ÖCB metabasites and some late Archean/early Proterozoic and Cambrian inheritance in Tieftal gneiss zircons suggest some involvement of old crustal rocks, too. The amount of crustal contamination can be calculated to be in the range of 10 to 40%. The Tieftal gneisses and the accompanying metabasites are interpreted as remnants of igneous rocks related to an early Ordovician rifting and incipient formation of new oceanic crust, an event which can be traced throughout the central and western European Variscan and Alpine terranes.

---

Alter, Genese und geologische Bedeutung einer polymetamorphen felsischen Intrusion im Ötztalkristallin, Tirol, Österreich

Zusammenfassung Ein Vielzahl von Methoden wurde angewandt, um das Alter und die Genese eines Orthogneiskörpers im mittleren Kaunertal, westliches Ötztalkristallin, abzuleiten. Der Tieftal-Orthogneiskörper ist eine in Amphiboliten eingeschaltete, intern differenzierte, polymetamorph überprägte, epizonale Intrusion. Er umfaßt sowohl leukokrate Hedenbergit-Hornblende-, Hornblende- und Biotit-Hornblende-Gneise als auch untergeordnet melanokrate Gesteine. Die leukokraten Tieftal-Gneise besitzen einen granitischen, beinahe einer eutektischen Schmelze entsprechenden Chemismus; einige Parameter wie hohe Na₂O+K₂O(8.07 bis 8.58Gew%), Zr(379 bis 554ppm) und Y(58 bis 79ppm) Gehalte weisen auf eine A-Typ Affinität hin. Die SEE-Spektren sind nur gering fraktioniert ((La/Yb)_N=2.3 bis 3.7) und weisen eine markante negative Eu-Anomalie auf. Einzelzirkon-Evaporationsdatierungen an 2 Proben und eine Sm-Nd Datierung von reliktischem magmatischem Titanit ergeben Alter von 487±7, 484±3 und 487±5Ma. Der gewichtete Mittelwert von 485±3Ma wird als das primäre magmatische Kristallisationsalter des Tieftal-Orthogneiskörpers interpretiert. Eine Rb-Sr Gesamtgesteinsdatierung ergibt eine gut definierte Regressionsgerade mit einem Alter von 411±9Ma. Dieses Alter beweist eine postmagmatische Störung des Rb-Sr Gesamtgesteinssystems, die durch die metamorphen überprägungen verursacht wurde. Die leukokraten Tieftal-Gneise besitzen eine relative primitive isotopische Zusammensetzung mit einem Nd _CHUR ^{485 Ma} a Wert von +1.7 und einem zurückgerechneten magmatischen⁸⁷Sr/⁸⁶Sr Initialverhältnis von 0.7047. Diese Daten machen eine große Beteiligung von Mantelmaterial wahrscheinlich. Am ehesten entstanden die leukokraten Tieftal-Gneise durch magmatische Fraktionierungsprozesse aus den Ausgangsgesteinen der begleitenden tholeiitischen Metabasite. Die noch primitivere isotopische Zusammensetzung der Metabasite im Ötztalkristallin und spätarchaische/frühproterozoische sowie kambrische Komponenten in den Zirkonen der leukokraten Tieftal-Gneise weisen aber auch auf die Beteiligung alten krustalen Materials hin. Der Anteil der krustalen Komponente liegt im Bereich von 10 bis 40%. Der Tieftal-Orthogneiskörper und die begleitenden Metabasite werden als Relikte magmatischer Gesteine, die während eines frühordovizischen Riftings und der beginnenden Bildung neuer ozeanischer Kruste entstanden sind, gedeutet. Zeugen dieses Vorganges sind in allen variszisch und alpidisch geprägten Gebieten Westund Mitteleuropas zu finden.

---

---

With 13 Figures     

Modular development of an inline monitoring system for waterborne pathogens in raw and drinking water

The state-of-the-art monitoring of drinking water hygiene is based on the cultivation and enumeration of indicator bacteria. Despite its proven reliability, this approach has the disadvantages of being (a) relatively slow and (b) limited to a small number of indicator bacteria. Ideally, alternative methods would be less time-consuming while providing information about a larger set of hygienically relevant microorganisms including viruses. In this paper, we present insights into the design of a modular concentration and detection system for bacteria, bacteriophages and viruses. Following further validation, this or similar techniques have the potential to extend and speed up the monitoring of raw and drinking water hygiene in the future. The system consists of different modules for the concentration of microorganisms, an amplification and detection unit that includes a module for the differentiation between live and dead microorganisms, and an automated system for decision support and self-diagnosis. The ongoing testing under controlled laboratory conditions and real-life conditions in the water supply industry yields further system improvements. Moreover, the increased sensitivity and broader range of microbiological parameters emphasize the need for a reconsideration of the currently used criteria for the assessment of (drinking) water hygiene.     

Real-time raster projection for web maps

The virtual globe is the default visualization for Digital Earth applications, but it can only show one half of the Earth. This article introduces user-adjustable, on-the-fly projection of georeferenced raster images for web mapping and web GIS applications. This technique allows users to center a map on arbitrary locations, while still seeing the entire Earth surface. Modern web mapping libraries can apply map projection transformations to vector data in the browser, but they do not currently support the projection of raster images at interactive frame rates. We use the cross-platform WebGL (Web Graphics Library) graphics pipeline for hardware-accelerated projection of raster data in web browsers. Two algorithmic techniques – inverse per-fragment projection, and combined forward per-triangle and inverse per-fragment projection – for georeferenced raster imagery are documented. The resulting raster maps combine the ease of use and flexibility of interactive virtual globes with the ability to show the entire Earth. The adjustable projection of raster imagery allows users to observe global phenomena that are difficult to visualize on virtual globes or traditional maps with static projections.