Comments on the Gravity and Magnetic Anomalies of Saros Bay using Wavelet Approach

In this paper, residual and regional gravity and magnetic anomalies of Saros Bay are separated using wavelet method. Wavelet is one of the modern, stochastic image approaches processing technique in electronics. The vertical, horizontal and diagonal components of wavelet output are evaluated simultaneously and an underground model is obtained by suitable cross-sections. Thus the geological and tectonic properties of Saros Bay are extracted. Our proposed model is confirmed by deep and shallow seismic researches of Turkish Petroleum Cooperation (TPAO). The South-East region of Saros Graben is formed by strike slip Ganos Fault and Anafartalar reverse Fault that lies on the east of Ganos Fault. On the Northwest region, there is strike slip component fault resulting in the Enez Graben. Here we detected two new oval type geological structure, both of them starting from the Eastern direction of Canakkale Dardanelles and but one of them continuing through Gelibolu Peninsuls, while the other tending to the Enez Graben. We think that these structures are the ruins of a very old sea in this region. This revised version was published online in November 2006 with corrections to the Cover Date.     

On the longitudinal distribution of the solar surges

After re-examination of data presented by Verma (1984), some misuse of solar parameteres published in Solar Geophysical Data is presented. Consequently some conclusions made by Verma (1984, 1986) on about the longitudinal distribution of cool solar surges based on these misuses are discussed.     

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

Mirdita Zone ophiolites and associated sediments in Albania reveal Neotethys Ocean origin

The Mirdita Ophiolite Zone in Albania is associated with widespread mélanges containing components of up to nappe-size. We dated matrix and components of the mélange by radiolarians, conodonts, and other taxa. The components consist of radiolarites, pelagic limestones and shallow-water limestones, all of Triassic age, as well as ophiolites. Triassic radiolarite as a primary cover of ophiolite material proves Middle Triassic onset of Mirdita ocean-floor formation. The mélange contains a turbiditic radiolarite-rich matrix (“radiolaritic flysch”), dated as Late Bajocian to Early Oxfordian. It formed as a synorogenic sediment during west-directed thrusting of ophiolite and sediment-cover nappes representing ocean floor and underplated fragments of the western continental margin. The tectonic structures formed during these orogenic events (“Younger Kimmeridian or Eohellenic Orogeny”) are sealed by Late Jurassic platform carbonates. The geological history conforms with that of the Inner Dinarides and adjoining areas; we therefore correlate the Mirdita-Pindos Ophiolite Zone with the Vardar Zone and explain its present position by far-distance west-directed thrusting.     

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.     

The Glarus thrust: excursion guide and report of a field trip of the Swiss Tectonic Studies Group (Swiss Geological Society, 14.–16. 09. 2006)

This excursion guide results form a field trip to the Glarus nappe complex organized by the Swiss Tectonic Studies Group in 2006. The aim of the excursion was to discuss old and recent concepts related to the evolution of the Glarus thrust. The major aspects were (i) the interplay between deformation, fluid flow and geochemical alteration, (ii) episodic versus continuous deformation and fluid flow, and (iii) the link between large-scale structures, microstructures, and geochemical aspects. Despite 150 years of research in the Glarus nappe complex and the new results discussed during the excursion, there exist controversies that still are unsolved.     

Provenance of the Bosnian Flysch

Sandwiched between the Adriatic Carbonate Platform and the Dinaride Ophiolite Zone, the Bosnian Flysch forms a c. 3000 m thick, intensely folded stack of Upper Jurassic to Cretaceous mixed carbonate and siliciclastic sediments in the Dinarides. New petrographic, heavy mineral, zircon U/Pb and fission-track data as well as biostratigraphic evidence allow us to reconstruct the palaeogeology of the source areas of the Bosnian Flysch basin in late Mesozoic times. Middle Jurassic intraoceanic subduction of the Neotethys was shortly followed by exhumation of the overriding oceanic plate. Trench sedimentation was controlled by a dual sediment supply from the sub-ophiolitic high-grade metamorphic soles and from the distal continental margin of the Adriatic plate. Following obduction onto Adria, from the Jurassic–Cretaceous transition onwards a vast clastic wedge (Vranduk Formation) was developed in front of the leading edge, fed by continental basement units of Adria that experienced Early Cretaceous synsedimentary cooling, by the overlying ophiolitic thrust sheets and by redeposited elements of coeval Urgonian facies reefs grown on the thrust wedge complex. Following mid-Cretaceous deformation and thermal overprint of the Vranduk Formation, the depozone migrated further towards SW and received increasing amounts of redeposited carbonate detritus released from the Adriatic Carbonate Platform margin (Ugar Formation). Subordinate siliciclastic source components indicate changing source rocks on the upper plate, with ophiolites becoming subordinate. The zone of the continental basement previously affected by the Late Jurassic–Early Cretaceous thermal imprint has been removed; instead, the basement mostly supplied detritus with a wide range of pre-Jurassic cooling ages. However, a c. 80 Ma, largely synsedimentary cooling event is also recorded by the Ugar Formation, that contrasts the predominantly Early Cretaceous cooling of the Adriatic basement and suggests, at least locally, a fast exhumation.     

The qualitative zoning record of minerals. A method for determining the duration of metamorphic events?

Summary The efficacy of ion exchange in petrological systems, for example the Fe-Mg exchange between garnet and biotite, is a function of grainsize, temperature, rates of temperature change and diffusion parameters. The combination of these variables determines the final zoning profile of minerals. Therefore, zoning profiles may be used to derive one of these variables if the others are known. For example, a mineral grain that experienced a short thermal event may still preserve a zoning profile characteristic of the heating path. In contrast, if that grain experienced a long thermal event it may develop a zoning profile typical of the cooling path. Conversely, for the same temperature-time cycle, large grains may not completely equilibrate at the metamorphic peak, and in smaller grains any previous record may have been erased. This is commonly observed in natural rocks where different grains within one thin section preserve often qualitatively different zoning profiles. Thus, a critical grainsize,l _crit, may exist that separates grains with qualitatively different zoning profiles so that grains of sizel > l _crit still retain zoning information about the heating path and grains of sizel <l _crit contain only information about the cooling path. If the critical grainsize can be measured and an independent estimate for the peak metamorphic temperature exists, the duration of the thermal event may, in principle, be estimated.The applicability of this method to natural garnets is hampered by a range of uncertainties. However, a parameterisation of the critical grain size may be used to illustrate the relative importance of grain size, temperature and event duration to the equilibration of minerals. In this paper, the critical grainsize is parameterized for various temperature-time cycles and its dependence on diffusion parameters is discussed. It is shown that, for Barrovian conditions,l _crit between 0.1 and 1 mm separates garnets retaining prograde information from garnets retaining cooling path information. For illustration, we compare the results with the critical grainsize of garnets from two metamorphic terranes, the Prydz Bay region (Antarctica) and the Koralm complex (Eastern Alps). Despite the large range of uncertainties attached to the method, it is shown that the critical grainsize of garnets in both terranes is consistent with a very short duration of the last thermal event that affected the two regions.

---

Die qualitative Form von Zonierungsprofilen in Mineralen. Information über die Dauer metamorpher Ereignisse?

Zusammenfassung Zonierungsprofile von Paragenesen im Ionenaustausch, zum Beispiel das Fe-Mg Austauschsgleichgewicht zwischen Granat und Biotit, sind eine Funktion von (i) Korngröße, (ii) Temperatur, (iii) Kühl- und Heizgeschwindigkeit sowie (iv) den Diffusionskonstanen. Es sollte daher möglich sein, Zonierungsprofile dazu zu verwenden, einen dieser Parameter zu bestimmen, wenn die anderen bekannt sind. Zonierungsprofile von Kristallen die ein kurzes thermisches Ereignis erfahren haben, mögen daher prograde Information beinhalten, wogegen derselbe Kristall ein retrogrades Profil aufweisen mag, wenn er einem langlebigen thermischen Ereignis unterlag (prograd und retrograd ist hier als: bei ansteigender Temperatur und bei abfallender Temperatur definiert). In Kristallen die dem gleichen Temperatur-Zeitpfad unterworfen waren, könnten kleine Körner nur die Kühlgeschichte dokumentieren, wogegen große noch Zonierungsinformation vom Heizpfad aufweisen. In natürlichen Gesteinen wird das oft dadurch beobachtet, daß verschiedene Körner auseinem Dünnschliff qualitativ verschiedene Zonierungsprofile aufweisen. Es ist daher möglich eine kritische Korngöße,l _crit, zu definieren, die Korngrößen mit qualitativ verschiedenen Zonierungsprofilen voneinander trennt. Körner mit einem Durchmesserl > l _crit haben, zumindest teilweise, noch prograde Profile, wogegen Körner mitl <l _crit nur retrograde Information dokumentieren. Wenn man diese kritische Korngröße messen kann, sollte sie dazu benutzt werden können, etwas über die Dauer des thermischen Ereignisses auszusagen.Die Anwendbarkeit dieser Methode ist durch eine Reihe von Fehlern limitiert. Nichtsdestotrotz, ist eine Parameterisierung der kritischen Korngröße nützlich um die relative Wichtigkeit von Korngröße, Temperatur und Zeit, für die Entwicklung von Zonierungsprofilen, zu illustrieren. Unser Modell kann dazu benutzt werden, um die Größenordnung von Granaten abzuschätzen, die noch prograde Information dokumentieren können. Die Abhängigkeit der kritischen Korngröße von Diffusionsparametern und der Form des Temperatur-Zeit Pfades wird ebenfalls diskutiert. Es wird gezeigt, daß für Metamorphose in der mittleren Grünschiefer und Amphibolit fazies,l _crit zwischen 0.1 and 1 mm prograd zonierte von retrograd zonierten Granaten trennt. Um das Modell zu illustrieren, werden die Ergebnisse auf die Prydz Bay Region (Antarktis) und den Koralm Komplex (Ostalpen) angewendet. Trotz der großen Fehlergrenzen der Methode glauben wir zeigen zu können, daß die kritische Korngröße beider Terrains andeutet, daß das letzte thermische Ereignis in beiden Gebieten kurzlebig war.

---

---

With 5 Figures     

Transition from eclogite to amphibolite-facies metamorphism in the Austroalpine Ulten Zone

Summary The Ulten Zone of the Austroalpine crystalline basement south-west of Meran (Italy) contains metapelitic schists and granoblastic paragneisses, leucocratic orthogneisses, migmatites (in both gneiss-lithologies), metabasites and ultramafic lenses. Metamorphic textures of the metapelitic schists and granoblastic paragneisses indicate two different metamorphic events, characterized by two mineral assemblages, which differ in mineral chemistry: (1) an eclogite facies mineral assemblage (M1) comprising Grt-Ky I-Bt. Ms-Kfs-PI-Qtz-Rt, and (2) an amphibolite facies mineral assemblage (M2) comprising Grt-KyII-Bt-Ms-PI-Qtz-Ilm±St. For the M1 event, pressures of at least 15kbar and temperatures of about 700°±50°C can be estimated. The later amphibolite facies overprint occurred at pressures of 6 to 8kbar and about 600°±50°C. The M1 and M2 assemblages belong to a continuous clockwise metamorphic evolution during the Variscan orogeny. Evidence for Alpine metamorphism can only be detected by sericite rims around kyanite and reset biotite ages. The migmatites, which contribute about 15–30vol.% of all rocks in the investigated area, were formed on the prograde path during the M1 event. Dissolution of H₂O in the melted part of the migmatites resulted in a CO₂dominated fluid, which was trapped in primary kyanite (M1) fluid inclusions. Secondary H₂O-rich fluid inclusions are found in quartz grains and may represent the fluid which enabled a pervasive equilibration during M2.

---

Übergang von eklogit-zu amphibolitfazieller Matamorphose in der austroalpinen Ultenzone

Zusammenfassung Die Ulten Zone, ein Teil des ostalpinen kristallinen basements, südwestlich von Meran, wird aus Metapeliten and granoblastischen Paragneisen, leukokraten Orthogneisen, Migmatiten (in beiden Lithologien), Metabasiten and ultramafischen Linsen aufgebaut. Metamorphe Texturen der Metapelite und granoblastischen Paragneise lassen auf zwei verschiedene metamorphe Ereignisse schließen, die durch unterschiedliche Mineral-chemismen und Paragenesen charakterisiert sind: (1) eine eklogitfazielle Paragenese (M1), bestehend aus Grt-KyI-Bt-Ms-Kfs-P1-Qtz-Rt und (2) eine amphibolitfazielle Paragenese (M2), bestehend aus Grt-KyII-Bt-Ms-P1-Qtz-Ilm±St. Für M1 konnten Minimaldrucke von 15kbar und Temperaturen von 700°±50°C abgeleitet werden. Die spätere amphibolitfazielle Überprägung fand bei 6 bis 8kbar und 600°±50°C statt. M1 und M2 gehören einer kontinuierlichen Metamorphoseentwicklung während der variszischen Orogenese an.Die Migmatite, ungefähr 15–30vol.% der Gesteine im untersuchten Gebiet, wurden am prograden Pfad während des M1 Ereignisses gebildet. Aufgrund der höheren Löslichkeit von H₂0 in der Schmelze, blieb ein CO₂, reiches Fluid zurück, das im primären Kyanit (M1) eingeschlossen wurde. Wässrige Flüssigkeitseinschlüsse können in Quarzkörnern gefunden werden. Dieses Fluid ist wahrscheinlich für die Reequilibrierung zu amphibolitfaziellen Bedingungen verantwortlich.

---

---

With 5 Figures     

Bayesian estimation of seismic hazard for two sites in Switzerland

Seismic hazard analysis is based on data and models, which both are imprecise and uncertain. Especially the interpretation of historical information into earthquake parameters, e.g. earthquake size and location, yields ambiguous and imprecise data. Models based on probability distributions have been developed in order to quantify and represent these uncertainties. Nevertheless, the majority of the procedures applied in seismic hazard assessment do not take into account these uncertainties, nor do they show the variance of the results. Therefore, a procedure based on Bayesian statistics was developed to estimate return periods for different ground motion intensities (MSK scale).Bayesian techniques provide a mathematical model to estimate the distribution of random variables in presence of uncertainties. The developed method estimates the probability distribution of the number of occurrences in a Poisson process described by the parameter . The input data are the historical occurrences of intensities for a particular site, represented by a discrete probability distribution for each earthquake. The calculation of these historical occurrences requires a careful preparation of all input parameters, i.e. a modelling of their uncertainties. The obtained results show that the variance of the recurrence rate is smaller in regions with higher seismic activity than in less active regions. It can also be demonstrated that long return periods cannot be estimated with confidence, because the time period of observation is too short. This indicates that the long return periods obtained by seismic source methods only reflects the delineated seismic sources and the chosen earthquake size distribution law.