March 1986 eruptive episodes at Piton de la Fournaise volcano (Reunion Island)

From May 1985 to April 1986 five discrete eruptions have occurred at Piton de la Fournaise volcano. On March the 17th, a sixth episode began with four distinct stages. They took place along the southeast rift zone of the volcano, from the summit to the sea coast. It was the first rift zone eruption in the south since 1800 A.D. and the first ever monitored at Piton de la Fournaise volcano.Three fissural vents opened at decreasing altitude emitting about 12 to 15 × 10⁶ m³ of olivine basalts between 19th March and 1st April. Strong seismic activity was accompanied by deformation of the summit area, and large-scale variations of the magnetic field. A summital event characterized the end of the flank activity with collapse of a new pit-crater and outflow of small amounts of degassed aphyric basalt.     

Gold dispersion in a tropical rainforest weathering profile at Dondo Mobi, Gabon

At Dondo Mobi in the gold district of Eteke, south Gabon, gold mineralization was studied in a tropical forest setting. The mineralization occurs in quartz veins within amphibolites of an Archaean gneiss-amphibolite series. Gold distribution patterns were studied in the different horizons of the weathering profile and in the different grain-size fractions of the materials sampled from three pits sunk in the weathered zone: upslope (pit P1), directly overlying (pit P2) and downslope (pit P3) of the mineralization.The weathering profile consists of an upper, thick, loose sandy argillaceous horizon (H1), an irregular nodular horizon with laterite nodules or blocks (H2) and a saprolite layer (H3) up to 70 m thick. In each pit, the specific geochemical signature of the bedrock is recognized in the three horizons of the weathering profile. Some groups of elements, e.g., Cr-Ni or Ba-V-P₂O₃, characterize amphibolite or black shale, respectively. Signal attenuations in the upper horizon can be explained by a homogenization effect. That can be related to a mushroom supergene dispersion, simultaneously inducing vertical decrease and lateral enrichment in element contents. Gold also mushrooms extensively in the different horizons of the weathering profile. The dispersion pattern is somewhat anisotropic, a strong enlargement being observed in the upper (H1) and intermediate nodular (H2) horizons alike. However, the evolution of Au distribution is not the same for all grain-size fractions: (a) in the finest fraction, Au is regularly distributed in the weathering profile at the anomalous top and tends to preserve an equivalent level of concentration laterally; (b) in the coarsest fraction, the maximum Au content is found in the upper horizon just above the mineralized structure; it is rapidly decreases laterally, in the surface halo.The Au geochemical dispersion halo consists mainly of fine-grained gold developed from rather coarse-grained mineralization. Its characteristics suggest that Au dispersion was produced as a result of Au dissolution together with collapse of the weathering profile causing gold migration over short distances.The minimal volumetric reduction of saprolite in the upper horizon can be estimated to be in the range of 60%, on the basis of the content of the different grain-size fractions and the density variations in the various horizons.In exploration, geochemical techniques, even using wide sample spacing, are proving quite adequate to prospect for Au in tropical rainforest weathering environment. Indeed, significant, well-contrasted dispersion haloes facilitate detection of rather confined mineralization.     

Application of the weighted residual method to diffraction by 2-D canyons of arbitrary shape: I. Incident SH waves

A solution for the two-dimensional scattering and diffraction of plane SH waves by canyons of arbitrary shape in an elastic half space is presented. The wave field for arbitrary geometry in this paper is computed numerically by the method of weighted residues (moment method). The wave displacement field computed by the present residual method for the case of a semi-circular canyon was shown to agree analytically and numerically with that computed by the exact closed form series solution. The same observations about ground amplifications, their dependence on frequencies and orientations of the incident waves, can be stated here for canyons of arbitrary shape as previously made for circular canyons.     

Inversion for seismic anisotropy using genetic algorithms1

A general inversion scheme based on a genetic algorithm is developed to invert seismic observations for anisotropic parameters. The technique is applied to the inversion of shear-wave observations from two azimuthal VSP data sets from the Conoco test site in Oklahoma. Horizontal polarizations and time-delays are inverted for hexagonal and orthorhombic symmetries. The model solutions are consistent with previous studies using trial and error matching of full waveform synthetics. The shear-wave splitting observations suggest the presence of a shear-wave line singularity and are consistent with a dipping fracture system which is known to exist at the test site. Application of the inversion scheme prior to full waveform modelling demonstrates that a considerable saving in time is possible whilst retaining the same degree of accuracy.     

Prediction of tension crack location and riverbank erosion hazards along destabilized channels

Riverbank erosion, associated sedimentation and land loss hazards are a land management problem of global significance and many attempts to predict the onset of riverbank instability have been made. Recently, Osman and Thorne (1988) have presented a Culmann-type analysis of the stability of steep, cohesive riverbanks; this has the potential to be a considerable improvement over previous bank stability theories, which do not account for bank geometry changes due to toe scour and lateral erosion. However, in this paper it is shown that the existing Osman-Thorne model does not properly incorporate the influence of tension cracking on bank stability since the location of the tension crack on the floodplain is indirectly determined via calculation or arbitrary specification of the tension crack depth. Furthermore, accurate determination of tension crack location is essential to the calculation of the geometry of riverbank failure blocks and hence prediction of land loss and bank sediment yield associated with riverbank instability and channel widening. In this paper, a rational, physically based method to predict the location of tension cracks on the floodplain behind the eroding bank face is presented and tested. A case study is used to illustrate the computational procedure required to apply the model. Improved estimates of failure block geometry using the new method may potentially be applied to improve predictions of bank retreat and floodplain land loss along river channels destabilized as a result of environmental change.     

Bias in estimating fractal dimension with the rescaled-range (R/S) technique

Fractal geostatistics are being applied to subsurface geological data as a way of predicting the spatial distribution of hydrocarbon reservoir properties. The fractal dimension is the controlling parameter in stochastic methods to produce random fields of porosity and permeability. Rescaled range (R/S)analysis has become a popular way of estimating the fractal dimension, via determination of the Hurst exponent (H). A systematic investigation has been undertaken of the bias to be expected due to a range of factors commonly inherent in borehole data, particularly downhole wireline logs. The results are integrated with a review of previous work in this area. Small datasets. overlapping samples, drift and nonstationariry of means can produce a very large bias, and convergence of estimates of H around 0.85–0.90 regardless of original fractal dimension. Nonstationarity can also account for H>1, which has been reported in the literature but which is theoretically impossible for fractal time series. These results call into question the validity of fractal stochastic models built using fractal dimensions estimated with the R/Smethod.     

On the kinetics of textural equilibration in forsterite marbles

In the aureole of the Beinn an Dubhaich granite, Skye, the minimum observed median forsterite-calcite-calcite dihedral angle varies from 110° at the olivinein isograd to about 165° (the equilibrium value) at the granite-limestone contact. Laboratory experiments were performed to investigate the kinetics of this textural change. It was found that the rate of change of the forsterite-calcite-calcite dihedral angle followed approximately first-order kinetics with an activation energy of 48±4 kJ mol^-1 for fluid-present conditions, and 90 ±4 kJ mol^-1 for fluid-absent conditions. Scanning ion imaging demonstrated that, at least in the early stages of textural change, solution-reprecipitation of the calcite was the rate determining step in the fluid-present runs. The latter result and the value of the activation energy are both consistent with the activation energy found by previous authors for (albeit zeroth order) silicate-aqueous solution solution-reprecipitation reactions. The value of activation energy for the dry data does not correspond to those for either grain boundary or volume diffusion of oxygen or magnesium in forsterite. The mechanism for textural equilibration in the fluid-absent runs is uncertain. Application of the experimentally-derived rate equation to the Beinn an Dubhaich marbles gave activation energies higher than those obtained experimentally. It is concluded from consideration of grain growth effects that activation energies derived from the Beinn an Dubhaich marbles probably reflect textural equilibration under predominantly fluid-absent conditions.