Geostatistical estimation of chemical contamination in stream sediments: The case study of Vale das Gatas mine (northern Portugal)

Based on an environmental geochemistry case study carried out in the neighbourhood of a W–Sn abandoned mine, the pollution in stream sediments was modelled through a Global Contamination Index. Such an index permits one to summarize the combination of deleterious elements in a single variable, obtained by the projection of samples onto the first axis of a PCASD (Principal Components Analysis of Standardized Data) applied to the entire n × p matrix containing the available concentrations of p = 16 elements in the set of n = 220 collected samples.In order to provide a sound basis for a coherent planning of the remediation process which will be put in operation in the affected area, it is necessary to balance the costs of reclaiming with the probabilities of exceeding the upper limits accepted for concentrations of environmentally harmful elements in sediments. Given these limits, they are back-transformed in the index values, providing a practical threshold between ‘clean’ and ‘contaminated’ samples. On the other hand, the minimum dimension of the cell to be reclaimed is restrained by the selected remediation process to be applied in the affected area. Hence, to meet the constraints of such a remediation process, it is required to estimate the probabilities of exceeding the index threshold in technologically meaningful sub-areas. For this end, the Indicator Block Kriging technique was applied, producing a series of maps where sub-areas to be reclaimed can be spotted for different probability levels. These maps, on which the decision making remediation agency can rely for its cost-benefit analysis, take into account both the spatial structure of ‘clean’ vs. ‘contaminated’ samples and the constraints of the reclaiming process.     

Methodology of hydrogeological characterization of deep carbonate aquifers as potential reservoirs of groundwater. Case of study: the Jurassic aquifer of <Emphasis Type="Italic">El Maestrazgo</Emphasis> (Castellón,Spain)

A methodology for the characterization of deep carbonate aquifers has been developed and applied to El Maestrazgo Jurassic aquifer in Castellón, Spain. Characterization of these aquifer formations, located at more than 300 m deep, consisted of a previous phase of compilation, analysis and synthesis of the existing information about the area, followed by a coordinated combination of different speciality studies: geology, stratigraphy, structural analysis, hydrogeology, hydrochemistry, geophysics and remote sensing. Geological studies included geological mapping, definition of stratigraphical units and facies and structural analysis. The aim of the hydrogeology study was to define aquifer formations, recharge area, aquifer points inventory and groundwater flow directions for the establishment of piezometric and water quality observation nets. Special techniques were applied, like thermal infrared aerial images and the evaluation of submarine groundwater discharge by means of natural radium isotopes. Hydrochemical techniques, including majority elements characterization and stable isotopes (¹⁸O, ²H and ³H) determination, allowed classifying hydrochemical facies and establishing a renewal pattern for water within the system. Geophysics was useful in determining the aquifer geometry, the features of the basement and the petrophysical characteristics of the geological formations. Preliminary results show an important tectonic complexity and the possibilities for groundwater uses in the area of study.     

Spatial Variability of Rock Depth in Bangalore Using Geostatistical, Neural Network and Support Vector Machine Models

Geospatial technology is increasing in demand for many applications in geosciences. Spatial variability of the bed/hard rock is vital for many applications in geotechnical and earthquake engineering problems such as design of deep foundations, site amplification, ground response studies, liquefaction, microzonation etc. In this paper, reduced level of rock at Bangalore, India is arrived from the 652 boreholes data in the area covering 220 km². In the context of prediction of reduced level of rock in the subsurface of Bangalore and to study the spatial variability of the rock depth, Geostatistical model based on Ordinary Kriging technique, Artificial Neural Network (ANN) and Support Vector Machine (SVM) models have been developed. In Ordinary Kriging, the knowledge of the semi-variogram of the reduced level of rock from 652 points in Bangalore is used to predict the reduced level of rock at any point in the subsurface of the Bangalore, where field measurements are not available. A new type of cross-validation analysis developed proves the robustness of the Ordinary Kriging model. ANN model based on multi layer perceptrons (MLPs) that are trained with Levenberg–Marquardt backpropagation algorithm has been adopted to train the model with 90% of the data available. The SVM is a novel type of learning machine based on statistical learning theory, uses regression technique by introducing loss function has been used to predict the reduced level of rock from a large set of data. In this study, a comparative study of three numerical models to predict reduced level of rock has been presented and discussed.     

SOME EXPERIMENTS CONCERNING THE PRIMARY SEISMIC PULSE*

Recordings were made with three types of detector of the primary compressional (P) and shear (S) wave pulses generated by explosions in boreholes. Charge weights varied from 0.08 kg to 9.5 kg and detector distances varied from about 3 m to about 80 m. Scaling by the simple factor W^1/3 where W is the charge weight, enabled observations from different sized charges to be fitted to a single expression. Experiments were carried out in the Bunter sandstone and the London clay and both fluid and solid tamping were used. This variation in tamping had no significant effect on the P-waves but it may have affected the generation of SV-waves. In both media the P-wave energy carried at 30 m from the shot by frequencies less than 100 Hz decreased rapidly with depth and was usually 1–2 % of the available chemical energy for a shot depth of 15 m. The S-wave energy was much less than this, but was highly directional. The P-wave pulse had the appearance of a damped sinusoid in very good agreement with the predictions of the ‘equivalent radiator’ hypothesis. However, the surface of this radiator should be identified not with the blown cavity but with the surface at which the tensile stresses associated with the stress wave become less than the tensile strength of the rock. The predominant frequency for a 1 kg charge at a depth of 15 m was 24 Hz in the clay and 52 Hz in the sandstone. In these and similar media, therefore, an effort should be made to keep individual charges less than 1 kg in reflection shooting and less than 10 kg in refraction shooting. The value of Q was about 50 in clay and about 25 in the sandstone. These estimates are rather uncertain because of the small distances over which the pulses were observed. The Z-transforms of the sampled pulses indicated that they were all of minimum phase, or very near to it.     

THE ELECTRICAL CURRENT PATTERN INDUCED BY AN OSCILLATING MAGNETIC DIPOLE IN A SEMI INFINITE CONDUCTIVE THIN PLATE*

In a paper by Koefoed and Kegge (1968), which was based on previous work of Wesley (1958), the electrical current pattern has been derived that is induced by an oscillating magnetic dipole in a semi-infinite thin plate of infinitesimal resistivity. In the present paper, the range of validity of the assumptions, on which the work of Wesley is based, is subjected to a theoretical analysis. It is found that the decisive factor for the validity of Wesley's derivation is the quotient of the square of the penetration depth of the electrical current over the product of the thickness of the plate and a distance that is indicative of the size of the current loops in the plate. Wesley's derivation is shown to be valid only when this factor is negligible. It is also shown that in this condition the imaginary component of the anomaly must be negligible. Model experiments are described in which the electrical current pattern is studied also in the range in which the derivation of Wesley is not valid. The procedure used in these model experiments was to measure the tangential component of the magnetic field strength very close to the metal plate that simulated the conductive dyke. In order to express the results of the measurements in terms of the imaginary to real ratio, these results are compared with an interpretation graph for field measurements that was published by Hedström and Parasnis (1958). It is found that the current pattern in the plate is essentially the same as that which follows from Wesley's derivation, provided that the imaginary to real ratio is less than one third. The measurements do not permit to draw conclusions regarding the current pattern in the plate in conditions that correspond to larger values of the imaginary to real ratio.     

PHOTOGRAMMETRIC MEASUREMENTS OF GLACIAL LANDFORMS

Photogrammetric methods using sequential aerial photography can provide the geomorphologist with comparative measurements of rapidly changing glacial landforms. Examples are given of eskers, kame and kettle areas, an ice-dammed lake and coastal features in south-east Iceland.     

ANALYTICAL AERIAL TRIANGULATION BLOCK ADJUSTMENT; THE DIRECT HEIGHT SOLUTION INCORPORATING TIE-STRIPS

The historical development of analytical aerial triangulation block adjustment is mentioned by way of introduction. The advantages of separating the plan and height adjustments are given and the observation equations of a typical adjustment section are analysed to show that this convenient practice is justified in normal circumstances. The principles of a direct solution of the normal equations suitable for a moderate-sized computer, are explained and developed in detail for the height adjustment. The advantages of tie-strips in the height adjustment are demonstrated and the difficulties which then arise are discussed. An approximate method which overcomes these difficulties is developed. Details of the computer programmes used are given in Appendix B.     

On the basic principles of geophysical radioactive measurements

Summary The process of radon diffusion in space and time through rocks is treated in detail from the point of view of radiometric prospecting for oil. Taking into account the depth and diffusivity conditions encountered in nature it is stated that vertical radon diffusion reaching from the oil deposits up to the surface cannot be expected. Time intervals as required for setting in of a stationary concentration state and decay conditions of the radon lead to the conclusion that the penetration of emanometric measurements seldom exceeds the depth of 5–10 m.