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.     

Channelized and hillslope sediment transport and the geomorphology of mountain belts

This paper uses the results of landscape evolution models and morphometric data from the Andes of northern Peru and the eastern Swiss Alps to illustrate how the ratio between sediment transport on hillslopes and in channels influences landscape and channel network morphologies and dynamics. The headwaters of fluvial- and debris-flow-dominated systems (channelized processes) are characterized by rough, high-relief, highly incised surfaces which contain a dense and hence a closely spaced channel network. Also, these systems tend to respond rapidly to modifications in external forcing (e.g., rock uplift and/or precipitation). This is the case because the high channel density results in a high bulk diffusivity. In contrast, headwaters where landsliding is an important sediment source are characterized by a low channel density and by rather straight and unstable channels. In addition, the topographies are generally smooth. The low channel density then results in a relatively low bulk diffusivity. As a consequence, response times are greater in headwaters of landslide-dominated systems than in highly dissected drainages. The Peruvian and Swiss case studies show how regional differences in climate and the litho-tectonic architecture potentially exert contrasting controls on the relative importance of channelized versus hillslope processes and thus on the overall geomorphometry. Specifically, the Peruvian example illustrates to what extent the storminess of climate has influenced production and transport of sediment on hillslopes and in channels, and how these differences are seen in the morphometry of the landscape. The Swiss example shows how the bedding orientation of the bedrock drives channelized and hillslope processes to contrasting extents, and how these differences are mirrored in the landscape. An erratum to this article can be found at     

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.     

Recommended Values of 239Pu, 240Pu and 239+240Pu Concentrations in Reference Material IAEA-315 (Marine Sediment) Estimated by Thermal Ionisation Mass Spectrometry, Inductively Coupled Plasma-Mass Spectrometry and Alpha Spectrometry

The recommended concentrations of ²³⁹Pu, ²⁴⁰Pu and ^239+240Pu in reference material IAEA‐315 (marine sediment) were estimated by three analytical methods: isotope dilution thermal ionisation mass spectrometry (TIMS), isotope dilution inductively coupled plasma‐mass spectrometry (ICP‐MS) and alpha spectrometry. The determination of ²³⁹Pu and ²⁴⁰Pu (^239+240Pu by alpha spectrometry) was carried out with samples from randomly selected bottles using each method. Plutonium‐238 was also measured by alpha spectrometry. A plutonium‐242 reference material was used as a spike for the quantitative analysis. The influence of ²⁴²Pu in the samples was therefore calculated; however, this contribution was less than the range of uncertainty and did not influence the final results. The obtained data were statistically analysed using variance component analysis and paired comparison. The combined standard uncertainties from “method/measurement”, “bottle” and “sub‐sample” were in the order of 3 to 6%. The main contributions to the uncertainty were from the material heterogeneity and from systematic differences between methods. Based on this study with twenty‐seven analyses using 10–14 g sample mass, concentrations of (38 ± 3) Bq kg^?1, (28 ± 3) Bq kg^?1 and (66 ± 4) Bq kg^?1 are proposed as recommended values for ²³⁹Pu, ²⁴⁰Pu and ^239+240Pu, respectively, and (9.5 ± 0.4) Bq kg^?1 for ²³⁸Pu as an information value in reference material IAEA‐315. In mass concentration units, these amount to (16.4 ± 1.2) ng kg^?1, (3.3 ± 0.4) ng kg^?1 and (0.015 ± 0.003) ng kg^?1 for ²³⁹Pu, ²⁴⁰Pu and ²³⁸Pu, respectively. The certified reference materials NIST 4350B and NIST 4354 were also analysed by TIMS for quality assurance of the method used in this study.     

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     

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.