Geochemical and Sr-Nd-Pb isotopic characteristics of Late Cenozoic leucite lamproites from the East European Alpine belt (Macedonia and Yugoslavia)

In the East European Alpine belt, leucite-sanidine-phlogopite-olivine-bearing volcanic rocks of Late Cenozoic age occur at eight localities within the Vardar suture zone and at one locality in the Southern Carpathian fold-and-thrust belt. Most of these volcanics are characterized by high Mg# (66.6–78.6), high abundances of Ni (117–373 ppm) and Cr (144–445 ppm) as well as high primary K₂O contents (5.63–7.01 %) and K₂O/Na₂O values (1.93–4.91). Rocks with more differentiated compositions are rare. A lamproite affinity of these rocks is apparent from their relatively low contents of Al₂O₃ (9.9–14.3 wt%) and CaO (6.2–8.3 wt%) in combination with high abundances of Rb (85–967 ppm), Ba (1,027–4,189 ppm), Th (18.9–76.5 ppm), Pb (19–54 ppm), Sr (774–1,712 ppm) and F (0.16–0.52 wt%), and the general lack of plagioclase. Although eruption of the magmas took place in post-collisional extensional settings, significant depletions of Nb and Ta relative to Th and La, low TiO₂ contents (0.92–2.17 %), low ratios of Rb/Cs, K/Rb and Ce/Pb as well as high ratios of Ba/La and Ba/Th suggest close genetic relationships to subduction zone processes. Whereas Sr and Nd isotope ratios show relatively large variations (⁸⁷Sr/⁸⁶Sr = 0.7078–0.7105, ¹⁴³Nd/¹⁴⁴Nd = 0.51242–0.51215), Pb isotope ratios display a very restricted range with ²⁰⁶Pb/²⁰⁴Pb = 18.68–18.88 and variable but generally high 7/4 (11–18) and 8/4 (65–95) values. The observed petrographic, geochemical and isotopic characteristics are best explained by a genetic model involving preferential melting of phlogopite-rich veins in an originally depleted lithospheric mantle source, whereby the metasomatic enrichment of the mantle source is tentatively related to the addition of components from subducted sediments during consumption of Tethyan oceanic lithosphere.Editorial responsibility: J. Hoefs     

Tropospheric corrections to SAR interferometry from GPS observations

Interferometric synthetic aperture radar (InSAR) techniques have been recognized as an ideal tool for many ground deformation monitoring applications. However, the spatially and temporally variable delay of the radar signal propagating through the atmosphere is a major limitation to accuracy. The dominant factor to be considered is the tropospheric heterogeneity, which can lead to misinterpretation of InSAR results. In this paper, a between-site (BS) and between-epoch (BE) double-differencing algorithm for the generation of tropospheric corrections to InSAR results based on GPS observations is tested. In order to correct the radar results on a pixel-by-pixel basis, the GPS-derived corrections have to be interpolated. Using experimental data it has been found that the inverse distance weighted and kriging interpolation methods are more suitable than the spline interpolation method. Differential corrections as large as several centimeters may have to be applied in order to ensure sub-centimeter accuracy for the InSAR result. The algorithm and procedures described in this paper could easily be implemented in a continuous GPS network data center. The interpolated image of BS, single-differenced tropospheric delays can be derived as a routine product to assist radar interferometry.     

A mixed-dimensional finite volume method for two-phase flow in fractured porous media

We present a vertex-centered finite volume method for the fully coupled, fully implicit discretization of two-phase flow in fractured porous media. Fractures are discretely modeled as lower dimensional elements. The method works on unstructured, locally refined grids and on parallel computers with distributed memory. An implicit time discretization is employed and the nonlinear systems of equations are solved with a parallel Newton-multigrid method. Results from two-dimensional and three-dimensional simulations are presented.     

Combined Earth orientation parameters based on homogeneous and continuous VLBI and GPS data

The CONT02 campaign is of great interest for studies combining very long baseline interferometry (VLBI) with other space-geodetic techniques, because of the continuously available VLBI observations over 2 weeks in October 2002 from a homogeneous network. Especially, the combination with the Global Positioning System (GPS) offers a broad spectrum of common parameters. We combined station coordinates, Earth orientation parameters (EOPs) and troposphere parameters consistently in one solution using technique- specific datum-free normal equation systems. In this paper, we focus on the analyses concerning the EOPs, whereas the comparison and combination of the troposphere parameters and station coordinates are covered in a companion paper in Journal of Geodesy. In order to demonstrate the potential of the VLBI and GPS space-geodetic techniques, we chose a sub-daily resolution for polar motion (PM) and universal time (UT). A consequence of this solution set-up is the presence of a one-to-one correlation between the nutation angles and a retrograde diurnal signal in PM. The Bernese GPS Software used for the combination provides a constraining approach to handle this singularity. Simulation studies involving both nutation offsets and rates helped to get a deeper understanding of this singularity. With a rigorous combination of UT1–UTC and length of day (LOD) from VLBI and GPS, we showed that such a combination works very well and does not suffer from the systematic effects present in the GPS-derived LOD values. By means of wavelet analyses and the formal errors of the estimates, we explain this important result. The same holds for the combination of nutation offsets and rates. The local geodetic ties between GPS and VLBI antennas play an essential role within the inter-technique combination. Several studies already revealed non-negligible discrepancies between the terrestrial measurements and the space-geodetic solutions. We demonstrate to what extent these discrepancies propagate into the combined EOP solution.     

Using Fuzzy Logic for the Analysis of Sea-level Indicators with Respect to Glacial-isostatic Adjustment: An Application to the Richmond-Gulf Region,Hudson Bay

An important constraint for the inference of mantle viscosity is the variation of the Holocene relative sea-level (RSL) height (with respect to today) following the last deglaciation. As a measure of this variation, sea-level indicators (SLIs) related to the RSL heights at specific past time epochs are used. For the inversion of the RSL-height change in terms of mantle viscosity, neighbouring SLIs may be grouped into an RSL diagram taken as representative for the region considered. Usually, the nominal height and age of a particular SLI are the only characteristics considered when determining the former RSL height. However, only SLIs based on isolation basins yield a narrow range for this height, whereas SLIs based on fossil samples provide a lower bound (shells), an upper bound (driftwood) or a finite interval (basal peat) for it. To also use fossil samples objectively, we develop a classification scheme of the depositional conditions based on fuzzy logic. After the definition of appropriate membership functions, this method leads to a systematic interpretation of the large number of SLIs available. We apply this method to SLIs from the Richmond-Gulf region, southeastern Hudson Bay, near the former glaciation center of Canada and derive a decay time of 5 ka for the exponential function best fitting the RSL diagram for this region.     

Genesis of itabirite-hosted Au–Pd–Pt-bearing hematite-(quartz) veins, Quadrilátero Ferrífero, Minas Gerais, Brazil: constraints from fluid inclusion infrared microthermometry, bulk crush-leach analysis and U–Pb systematics

Fluid inclusions hosted in quartz and specular hematite from auriferous (jacutinga) and barren veins in the Quadrilátero Ferrífero (QF) have been studied using conventional and near infrared microscopy, respectively. The mineralization consists of veins that cross-cut metamorphosed iron formation (itabirite) of the Paleoproterozoic Itabira Group. The sample suite comprises hematite from veins from the low-strain domain in the W and SW of the study area, as well as hematite samples from the eastern high-strain domain in the central and NE parts of the QF. Halogen ratios of fluid inclusions in quartz and hematite from all studied deposits are consistent with a fluid evolved from dissolving and reprecipitating halite that was subsequently diluted. Fluid inclusions hosted in quartz and hematite are characterized by consistent Na/K ratios and considerable SO₄ contents, and suggest similar formation conditions and, perhaps, fluid origin from a common source. Na/K and Na/Li fluid mineral geothermometers indicate water–rock interaction at approximately 340±40°C. Hematites from the high-strain domain contain fluid inclusion assemblages of high-temperature aqueous-carbonic and multiphase high-salinity, high-temperature aqueous inclusions probably due to fluid immiscibility in the system H₂O–NaCl–CO₂. Fluid inclusions hosted in hematite from barren veins in the low-strain domain, as well as in hematite from jacutinga-type mineralization from the central part of the QF, only host multiphase aqueous fluid inclusions all showing narrow ranges of salinity (7.2–11.7 wt.% NaCl equiv.) and homogenization temperatures (148 to 229°C). Lower homogenization temperatures and the absence of CO₂-rich inclusions in specular hematite from these occurrences are attributed to carbonate precipitation and/or CO₂ escape due to cooling during fluid migration from the high- to the low-strain domain. Pb–Pb and U–Pb systematics of gold, hematite and hematite-hosted fluid inclusions in combination with geochemical evidence indicate distinct sources for Pd, Au, and Pb. The formation of specular hematite veins may be related to retrograde metamorphic fluids being released during the Brazilian orogenic cycle (600–700 Ma). The Pb isotopic characteristics of all samples are readily reconciled in a simple model that involves two different Paleoproterozoic or Archean source lithologies for lead and reflects contrasting depths of fluid percolation during the Brasiliano orogeny.     

Loading fractures and Liesegang laminae: new sedimentary structures found in the north-western North Alpine Foreland Basin (Oligocene–Miocene, south-west Germany)

Abstract Spectacular sedimentary structures recently found in the Molasse Basin (Oligocene–Miocene) in southern Germany were produced by soft‐sediment deformation under highly unusual conditions. These large, apparently wedge‐like structures –‘loading fractures’– cut down into beds of marl and are filled with coarse sand and intraclasts of shale. Wrapping the sides of the structures is a thin, continuous bed of layered dark claystone – the ‘DCB’. The upper and lower layers of this bed are an organic‐rich clay; the middle layer is a laminated quartzite. The precursor of the DCB was a lacustrine gyttja rich in diatom frustules. It was supersaturated in silica as it was buried. Subsequent diffusion of oxygen into this gyttja at a burial depth of only a few metres resulted in the formation of Liesegang laminae of quartz. These laminae grew and amalgamated, forming the layer of laminated quartzite. The sediments overlying the DCB were eventually removed by erosion, probably in a high‐energy marine environment. This erosion cut down to the DCB but was unable to penetrate it. The DCB remained exposed on the sea floor until a sudden depositional event occurred – the deposition of a 2·5 metre thick bed of coarse sand with shale intraclasts. Although the DCB had been able to resist the submarine erosion, it could not support the load of this new bed. The quartzite layer in it therefore fractured, transferring that load down onto the underlying, still‐unconsolidated marl. The intraclast‐rich sands were forced down into this marl, carrying ahead of them the partly broken remains of the DCB.     

Audiomagnetotelluric investigations of shallow sedimentary basins in northern Sudan

Audiomagnetotelluric (AMT) investigations were carried out to study shallow sedimentary basins in northern Sudan and southern Egypt. Unlike recordings in settled areas, the natural variations (atmospherics) in the eastern Sahara were completely unperturbed by cultural noise, usually yielding very high data quality. Typical recordings and spectra are presented. The atmospherics show distinct daily variations, with maximum amplitudes in the local afternoon/night hours. As a sort of calibration, AMT measurements together with DC-resistivity and electromagnetic soundings were conducted at several groundwater boreholes. Under favourable conditions, a vertical anisotropy may be resolved by joint inversion of DCR and AMT data. The results of a case study in the Es Safya basin (NW Sudan) are presented. The investigations were aimed at determining basement depths and groundwater levels, as well as the direction and dip of major bounding faults. 2-D modelling and inversion revealed the overall conductivity structure of this shallow graben. The AMT method proved to be a valuable tool for a comparatively swift reconnaissance study of shallow structures.     

Local Multiscale Approximations of Geostrophic Oceanic Flow: Theoretical Background and Aspects of Scientific Computing

In modern geoscience, understanding the climate depends on the information about the oceans. Covering two-thirds of the Earth, oceans play an important role. Oceanic phenomena are, for example: oceanic circulation; water exchanges between atmosphere, land, and ocean; or temporal changes of the total water volume. All these features require new methods in constructive approximation, since they are regionally bounded and not globally observable. This article deals with a new and alternative method of handling data with locally supported basis functions, modeling them in a multiscale scheme involving a wavelet approximation and presenting the main results for the dynamic topography and the geostrophic flow, e.g., in the Northern Atlantic. Further, it is demonstrated that compressional rates of the occurring wavelet transforms can be achieved by using of locally supported wavelets and investigating the signal distribution within different frequency bands.     

Shock-induced growth and metastability of stishovite and coesite in lithic clasts from suevite of the Ries impact crater (Germany)

The microtextures of stishovite and coesite in shocked non-porous lithic clasts from suevite of the Ries impact structure were studied in transmitted light and under the scanning electron microscope. Both high-pressure silica phases were identified in situ by laser-Raman spectroscopy. They formed from silica melt as well as by solid-state transformation. In weakly shocked rocks (stage I), fine-grained stishovite (≤1.8 μm) occurs in thin pseudotachylite veins of quartz-rich rocks, where it obviously nucleated from high-pressure frictional melts. Generally no stishovite was found in planar deformation features (PDFs) within grains of rock-forming quartz. The single exception is a highly shocked quartz grain, trapped between a pseudotachylite vein and a large ilmenite grain, in which stishovite occurs within two sets of lamellae. It is assumed that in this case the small stishovite grains formed by the interplay of conductive heating and shock reverberation. In strongly shocked rocks (stages Ib–III, above ∼30 GPa), grains of former quartz typically contain abundant and variably sized stishovite (<6 μm) embedded within a dense amorphous silica phase in the interstices between PDFs. The formation of transparent diaplectic glass in adjacent domains results from the breakdown of stishovite and the transformation of the dense amorphous phase and PDFs to diaplectic glass in the solid state. Coesite formed during unloading occurs in two textural varieties. Granular micrometre-sized coesite occurs embedded in silica melt glass along former fractures and grain boundaries. These former high-pressure melt pockets are surrounded by diaplectic glass or by domains consisting of microcrystalline coesite and earlier formed stishovite. The latter is mostly replaced by amorphous silica.