ALBERT

Description: The final report contains a description of the results obtained within a research contract between IRE RAS and GFZ Potsdam during the period April-November 2004. The objectives of investigation included (1) the radio-holographic methods for obtaining vertical profiles of the vertical gradients of physical parameters in the atmosphere, (2) radio holographic methods for atmospheric, ionospheric and stratospheric waves, and (3) validation of the software with GPS/MET (GPS/METeorology, e.g., Rocken et al. 1997) and CHAMP (CHAllenging Minisatellite Payload, e.g., Reigber et al. 2005) data and final report with recommendations.

Description: Late Miocene to Quaternary volcanic rocks from the frontal arc to the back-arc region of the Central Volcanic Zone in the Andes show a wide range of delta 11B values (+4 to -7 ‰) and boron concentrations (6 to 60 ppm). Positive delta 11B values of samples from the volcanic front indicate involvement of a 11B-enriched slab component, most likely derived from altered oceanic crust, despite the thick Andean continental lithosphere, and rule out a pure crust-mantle origin for these lavas. The delta 11B values and B concentrations in the lavas decrease systematically with increasing depth of the Wadati-Benioff Zone. This across-arc variation in delta 11B values and decreasing B/Nb ratios from the arc to the back-arc samples are attributed to the combined effects of B-isotope fractionation during progressive dehydration in the slab and a steady decrease in slab-fluid flux towards the back arc, coupled with a relatively constant degree of crustal contamination as indicated by similar Sr, Nd and Pb isotope ratios in all samples. Modelling of fluid-mineral B-isotope fractionation as a function of temperature fits the across-arc variation in delta 11B and we conclude that the B-isotope composition of arc volcanics is dominated by changing delta 11B composition of B-rich slab-fluids during progressive dehydration. Crustal contamination becomes more important towards the back-arc due to the decrease in slab-derived fluid flux. Because of this isotope fractionation effect, high delta 11B signatures in volcanic arcs need not necessarily reflect differences in the initial composition of the subducting slab. Three-component mixing calculations for slab-derived fluid, the mantle wedge and the continental crust based on B, Sr and Nd isotope data indicate that the slab-fluid component dominates the B composition of the fertile mantle and that the primary arc magmas were contaminated by an average addition of 15 to 30 % crustal material.

Description: This manual outlines the characteristics and structure of the software and describes how to use the software. The principles and new features are outlined systematically and referred partly to existing references. Numerical examples of multi-functions and internal tests as well as external comparisons are given.

Description: We review the historical, geological, tide-gauge, GPS and gravimetric evidence advanced in favour of or against continuing land uplift around Hudson Bay, Canada. After this, we reanalyse the tide-gauge and GPS data for Churchill using longer time series than those available to previous investigators. The dependence of the mean rate of relative sea-level change obtained on the length and mid-epoch of the observation interval considered is investigated by means of the newly developed linear-trend analysis diagram. For studying the shorter-period variability of the tide-gauge record, the continuous-wavelet transform is used. The mean rate of land uplift obtained from GPS is based on a new analysis using IGS solutions of GFZ. Furthermore, sea-level indicators from the Churchill region representing the relative sea-level history during the past 8000 a are included. Finally, the four types of observable are jointly inverted in terms of mantle viscosity. The optimum values are 3×10^20 Pa s and 1.6 × 10^22 Pa s for the upper- and lower-mantle viscosities, respectively.

Description: This report is also published as Scientific Report No. 04-9 of the Danish Meteorological Institute.

Description: As observed on the Earth's surface, the magnetic field can be separated into three components: The "Main Field", the "Crustal Field" and the "External Magnetic Field". This work concentrates on the magnetospheric current systems, in order to correct the effect of large-scale magnetospheric current systems. Their exact knowledge is vital to improve the main field and secular variation models. Spherical harmonic analysis (SHA) is commonly used to describe the magnetic field. In SHA it is possible to distinguish between internal and external contributions to the Earth's magnetic field. The spherical harmonic coefficients describe dipolar, quadrupolar and higher parts in spherical coordinates. Here, CHAMP (CHAllenging Minisatellite Payload) measurements are used. Due to the inclination CHAMP is changing the orbit in local time. From time to time CHAMP is flying in resonance with the Earth's rotation, i.e. the satellite is flying over the same ground track after a couple of days. To model changes during a period of a few days these so called "repeat tracks" are very helpful for an improved main field/external field model. On the Earth's surface the magnetic effect of magnetospheric currents is often characterised by the DST-indix. In this study, it is explained how a possible constitution of the external sources of the geomagnetic field can look like and how the effect of these magnetospheric currents can be corrected in main field modelling.

Description: The Dead Sea Transform (DST) is a prominent shear zone in the Middle East. It separates the Arabian plate from the Sinai microplate and stretches from the Red Sea rift in the south via the Dead Sea to the Taurus-Zagros collision zone in the north. Formed in the Miocene »17 Ma ago and related to the breakup of the Afro-Arabian continent, the DST accommodates the left-lateral movement between the two plates. The study area is located in the Arava Valley between the Dead Sea and the Red Sea, centered across the Arava Fault (AF), which constitutes the major branch of the transform in this region. A set of seismic experiments comprising controlled sources, linear profiles across the fault, and specifically designed receiver arrays reveals the subsurface structure in the vicinity of the AF and of the fault zone itself down to about 3–4 km depth. A tomographically determined seismic P velocity model shows a pronounced velocity contrast near the fault with lower velocities on the western side than east of it. Additionally, S waves from local earthquakes provide an average P -to-S velocity ratio in the study area, and there are indications for a variations across the fault. High-resolution tomographic velocity sections and seismic reflection profiles confirm the surface trace of the AF, and observed features correlate well with fault-related geological observations. Coincident electrical resistivity sections from magnetotelluric measurements across the AF show a conductive layer west of the fault, resistive regions east of it, and a marked contrast near the trace of the AF, which seems to act as an impermeable barrier for fluid flow. The correlation of seismic velocities and electrical resistivities lead to a characterisation of subsurface lithologies from their physical properties. Whereas the western side of the fault is characterised by a layered structure, the eastern side is rather uniform. The vertical boundary between the western and the eastern units seems to be offset to the east of the AF surface trace. A modelling of fault-zone reflected waves indicates that the boundary between low and high velocities is possibly rather sharp but exhibits a rough surface on the length scale a few hundreds of metres. This gives rise to scattering of seismic waves at this boundary. The imaging (migration) method used is based on array beamforming and coherency analysis of P -to-P scattered seismic phases. Careful assessment of the resolution ensures reliable imaging results. The western low velocities correspond to the young sedimentary fill in the Arava Valley, and the high velocities in the east reflect mainly Precambrian igneous rocks. A 7 km long subvertical scattering zone (reflector) is offset about 1 km east of the AF surface trace and can be imaged from 1 km to about 4 km depth. The reflector marks the boundary between two lithological blocks juxtaposed most probably by displacement along the DST. This interpretation as a lithological boundary is supported by the combined seismic and magnetotelluric analysis. The boundary may be a strand of the AF, which is offset from the current, recently active surface trace. The total slip of the DST may be distributed spatially and in time over these two strands and possibly other faults in the area.

Description: Ein transversales Störungssystem im Nahen Osten, die Dead Sea Transform (DST), trennt die Arabische Platte von der Sinai-Mikroplatte und erstreckt sich von Süden nach Norden vom Extensionsgebiet im Roten Meer über das Tote Meer bis zur Taurus-Zagros Kollisionszone. Die sinistrale DST bildete sich im Miozän vor »17 Ma und steht mit dem Aufbrechen des Afro-Arabischen Kontinents in Verbindung. Das Untersuchungsgebiet liegt im Arava Tal zwischen Totem und Rotem Meer, mittig über der Arava Störung (Arava Fault, AF), die hier den Hauptast der DST bildet. Eine Reihe seismischer Experimente, aufgebaut aus künstlichen Quellen, linearen Profilen über die Störung und entsprechend entworfenen Empfänger-Arrays, zeigt die Untergrundstruktur in der Umgebung der AF und der Verwerfungszone selbst bis in eine Tiefe von 3–4 km. Ein tomographisch bestimmtes Modell der seismischen Geschwindigkeiten von P-Wellen zeigt einen starken Kontrast nahe der AF mit niedrigeren Geschwindigkeiten auf der westlichen Seite als im Osten. Scherwellen lokaler Erdbeben liefern ein mittleres P -zu-S Geschwindigkeitsverhältnis und es gibt Anzeichen für Änderungen über die Störung hinweg. Hoch aufgelöste tomographische Geschwindigkeitsmodelle bestätigen der Verlauf der AF und stimmen gut mit der Oberflächengeologie überein. Modelle des elektrischen Widerstands aus magnetotellurischen Messungen im selben Gebiet zeigen eine leitfähige Schicht westlich der AF, schlecht leitendes Material östlich davon und einen starken Kontrast nahe der AF, die den Fluss von Fluiden von einer Seite zur anderen zu verhindern scheint. Die Korrelation seismischer Geschwindigkeiten und elektrischer Widerstände erlaubt eine Charakterisierung verschiedener Lithologien im Untergrund aus deren physikalischen Eigenschaften. Die westliche Seite lässt sich durch eine geschichtete Struktur beschreiben, wogegen die östliche Seite eher einheitlich erscheint. Die senkrechte Grenze zwischen den westlichen Einheiten und der östlichen scheint gegenüber der Oberflächenausprägung der AF nach Osten verschoben zu sein. Eine Modellierung von seismischen Reflexionen an einer Störung deutet an, dass die Grenze zwischen niedrigen und hohen Geschwindigkeiten eher scharf ist, sich aber durch eine raue Oberfläche auf der Längenskala einiger hundert Meter auszeichnen kann, was die Streuung seismischer Wellen begünstigte. Das verwendete Abbildungsverfahren (Migrationsverfahren) für seismische Streukörper basiert auf Array Beamforming und der Kohärenzanalyse P -zu-P gestreuter seismischer Phasen. Eine sorgfältige Bestimmung der Auflösung sichert zuverlässige Abbildungsergebnisse. Die niedrigen Geschwindigkeiten im Westen entsprechen der jungen sedimentären Füllung im Arava Tal, und die hohen Geschwindigkeiten stehen mit den dortigen präkambrischen Magmatiten in Verbindung. Eine 7 km lange Zone seismischer Streuung (Reflektor) ist gegenüber der an der Oberfläche sichtbaren AF um 1 km nach Osten verschoben und lässt sich im Tiefenbereich von 1 km bis 4 km abbilden. Dieser Reflektor markiert die Grenze zwischen zwei lithologischen Blöcken, die vermutlich wegen des horizontalen Versatzes entlang der DST nebeneinander zu liegen kamen. Diese Interpretation als lithologische Grenze wird durch die gemeinsame Auswertung der seismischen und magnetotellurischen Modelle gestützt. Die Grenze ist möglicherweise ein Ast der AF, der versetzt gegenüber des heutigen, aktiven Asts verläuft. Der Gesamtversatz der DST könnte räumlich und zeitlich auf diese beiden Äste und möglicherweise auch auf andere Störungen in dem Gebiet verteilt sein.

Description: Recently released global gravity field models generated solely from CHAMP and GRACE satellite observations allow with an unprecedented accuracy and resolution the recovery of the mean sea surface topography from the difference between an altimetry-based mean sea surface height model and the gravity model's derived geoid. Here the CHAMP EIGEN-2 gravity field model, and the first GFZ GRACE gravity model, EIGEN-GRACE01S, are used. The mean sea surface height model has been compiled from four years'; worth of TOPEX altimeter data. To evaluate the accuracy and resolution limits of the CHAMP and GRACE geoids for the envisaged application, a low pass filter in the spatial domain with different cut-off wavelengths has been applied to the geoid and sea surface data before subtraction. The minimum wavelength, where noisy and erroneous features in the recovered sea surface topography are minimised, can be interpreted as an indicator for the best suited common spatial resolution. The EIGEN-2 model's geoid has been tested to have a resolution of 1800 km, which corresponds to a truncation degree of l = 22 in terms of spherical harmonics. Using the EIGEN-GRACE01S model, the resolution could be extended to 1000 km (l = 40). These boundaries can be attributed to the geoid's error, exceeding 2 cm in case of the CHAMP model, and in case of the GRACE model to spurious systematic signals increasing with increasing spherical harmonic degree. The calculated sea surface topography models have been used to derive absolute geostrophic sea surface velocities. An error propagation shows that the requirement of 1 cm/s for geoid induced velocity errors is fulfilled at the given resolutions for all latitudes excluding a narrow equatorial band. Maximum geostrophic velocities are derived in the 1000 km-resolution model for the Kuroshio region with 40 cm/s, and for the Gulf Stream east off Cape Hatteras with 25 cm/s.