Structure of the upper mantle in the Northeastern Atlantic close to the azores-gibraltar ridge from surface-wave and body-wave observations

Phase velocities of teleseismic Rayleigh waves have been measured in the central North Atlantic on both sides of the Azores-Gibraltar Ridge (AGR) by means of a specially designed long-period station network. The dispersion data obtained were regionalized and then subjected to a “hedgehog” inversion, which gives a set of upper mantle models compatible with the observational data within specified error bounds.Reasonable model solutions were selected by using regional body-wave observations, such as P_n- and S_n-wave velocities determined from earthquakes along the AGR. The S_(itn) velocities measured indicate that the shear-wave velocity in the mantle part of the lithosphere is much higher on the northern side of the AGR. Strongly negative P-wave residuals in this area indicate faster seismic propagation than implied by the Jeffreys-Bullen travel-time tables, while propagation is much slower in the Gulf of Cadiz area. Furthermore the residuals show a clear difference for paths through oceanic and continental domains and suggest that the transition between these two domains extends much further into the ocean on the southern side of the AGR than on the northern side.The proposed model for the structure of the upper mantle in that region shows that there exists a pronounced velocity contrast across the AGR. Thickening of the lithospheric plate with increasing plate age is indicated to the south of the ridge. The greatest thickness is reached close to the continental margin within a zone about 500 km wide, whose velocity close to the Canary Islands and Madeira is significantly lower, probably due to the well-known volcanic activity there. These observations together with the travel time residuals reveal that this zone seems to be of a transitional nature somewhere between a continental and oceanic structure.     

Thermobarometry of the Bleikvassli Zn-Pb-(Cu) deposit, Nordland, Norway

The Bleikvassli massive sulfide ore deposit is hosted by Proterozoic pelitic, quartzofeldspathic, and amphibolitic rocks of the Uppermost Allochthon of the Scandinavian Caledonides. Staurolite-garnet-biotite and kyanite-staurolite-biotite assemblages indicate that metamorphism reached the kyanite zone of the amphibolite facies. Geothermobarometry was conducted on rocks in and around the deposit using a variety of silicate and sulfide calibrations. Temperature determinations are most reliant on the garnet-biotite exchange reaction, with analyses obtained from 259 garnet rims and adjacent biotite. Results from nine calibrations of the garnet-biotite geothermometer are considered, but compositional limitations of many calibrations involving high Ca and Mn contents in garnet and Al^VI and Ti in biotite make many of the coexisting mineral pairs unsuitable. Average temperatures calculated from the two calibrations that most closely address the garnet-biotite compositions observed at Bleikvassli are 584 °C ± 49 °C and 570 °C ± 40 °C. The application of two calibrations of the garnet-staurolite geothermometer on a limited number of samples yields 581 °C ± 2 °C and 589 °C ± 12 °C, assuming a _H2O=0.84, based upon calculations of the modal proportions of gaseous species. Pressure determinations are less constrained. Phengite and plagioclase-biotite-garnet-muscovite geobarometers give average pressures of approximately 5.0 kbar and 8.5 ± 1.2 kbar, respectively. Pressures obtained from the sphalerite-hexagonal pyrrhotite-pyrite barometer average 7.7 ± 1.0 kbar. In consideration of these results, the peak metamorphic conditions at the Bleikvassli deposit are estimated to be 580 °C and 8 kbar. Received: 18 June 1997 / Accepted: 14 May 1998     

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity I: the continuous in time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a continuous in-time setting for a strictly positive constrained specific storage coefficient. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation, and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.     

Interparticle forces and displacements in granular materials

In a micromechanics framework, the main issue is the relationship between the microscale variables and the macroscale variables. These variables are used to describe either the statics or kinematics of the system. The relationships can be classified in two ways, namely, the “averaging” relationships and the “tracking” relationships. The averaging relationships express the macroscale variable as an averaging of the microscale variables; for example, the stress as a function of contact forces. The “tracking” relationships express the microscale variable as a function of the macroscale variables; for example, the contact force at a given orientation as a function of the stress. Based on fundamental premises, a unique averaging relationship exists for either the statics or the kinematics case. However, it is generally impossible to have a unique expression of the “tracking” relationship because they are generally derived with certain assumptions. In this paper, we will present expressions of the “tracking” based on three different approaches, namely, (1) energy conservation principle, (2) representation theory, and (3) indirect scheme. The assumptions used in each approach are discussed. The results are compared among the three approaches as well as that obtained from the Discrete Element Method (DEM).     

Late thermal events based on zircon fission track ages in northeastern New South Wales and southeastern Queensland: Links to Sydney Basin seismicity?

Zircon, concentrated from basaltic terrains in northeastern New South Wales and southeastern Queensland, reveals some unexpectedly young fission track peaks. Between 2 to 13 Ma in age, these peaks are younger than known Tertiary basaltic ages from these regions which match older fission track peaks. Analysis of the fission track data suggests that the young dates are probably not reset ages due to recent heating events such as bush fires, but more likely mark thermal resetting by later volcanic eruptions.

The young ages decrease southwards from Queensland through northern New South Wales and trend toward seismic zones within the Sydney Basin in the Newcastle, Blue Mountains and Illawarra regions. A model based on northward motion of the Australian plate over a hot asthenospheric source (0.75° latitude/Ma over 12 Ma)) predicts the positions of most young zircon ages to within ± 70 km in latitude when projected from seismic sites at Newcastle and Bowral‐Robertson.

A minor hot spot source is proposed, some 250 km across, which triggered isolated basaltic and zircon‐bearing eruptions every few million years and now underlies the southern Sydney Basin. This would bring Sydney Basin seismicity into line with other seismic zones known at present hot spot positions across southeastern Australia and the Tasman Sea. It raises questions concerning activation of local seismicity, potential for volcanic risk and distribution of young uplift in the Sydney region. Similar studies are needed to test other puzzling seismic zones (e.g. the Dalton‐Gunning Zone).     

A hybrid indirect boundary element—discrete wave number method applied to simulate the seismic response of stratified alluvial valleys

A hybrid indirect boundary element – discrete wavenumber method is presented and applied to model the ground motion on stratified alluvial valleys under incident plane SH waves from an elastic half-space. The method is based on the single-layer integral representation for diffracted waves. Refracted waves in the horizontally stratified region can be expressed as a linear superposition of solutions for a set of discrete wavenumbers. These solutions are obtained in terms of the Thomson–Haskell propagators formalism. Boundary conditions of continuity of displacements and tractions along the common boundary between the half-space and the stratified region lead to a system of equations for the sources strengths and the coefficients of the plane wave expansion. Although the regions share the boundary, the discretization schemes are different for both sides: for the exterior region, it is based on the numerical and analytical integration of exact Green's functions for displacements and tractions whereas for the layered part, a collocation approach is used. In order to validate this approach results are compared for well-known cases studied in the literature. A homogeneous trapezoidal valley and a parabolic stratified valley were studied and excellent agreement with previous computations was found. An example is given for a stratified inclusion model of an alluvial deposit with an irregular interface with the half-space. Results are displayed in both frequency and time domains. These results show the significant influence of lateral heterogeneity and the emergence of locally generated surface waves in the seismic response of alluvial valleys.     

Buchbesprechungen

Ohne Zusammenfassung     

Sedimentary noble gases

This paper presents experimental data for the trapped noble gas contents in a variety of shales and related samples. These data, along with those previously available, indicate that the normal sedimentary rock pattern is one of progressive enrichment of the heavier gases in comparison with the proportions in air, from which sedimentary rock gases were presumably acquired. Within this normal pattern, however, variations considerably beyond an order of magnitude characterize both abundance and composition, so that multiple processes appear to be necessary to account for trapping of noble gases in sedimentary rocks. Relative to the normal pattern a number of cases of anomalous Ne enhancement occur. This unexplained phenomenon, previously thought to be exotic, must now be considered unexceptional.Noble gases in air occur in approximately the ‘planetary’ proportions characteristic of meteorites, except for more than an order of magnitude deficiency of Xe. It is commonly held that atmospheric noble gases are indeed planetary, the Xe deficiency in air being made up by the inventory in sedimentary rocks. We consider that while this hypothesis cannot be proved false the available data do not support it. At least equal weight should be given to the alternatives: either that (preferentially) Xe has failed ever to become part of the atmosphere or that the total terrestrial noble gas inventory is simply not planetary.