Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

The effect of naturally acidified irrigation water on agricultural volcanic soils. The case of Asembagus, Java, Indonesia

Acid water from the Banyuputih river (pH  3.5) is used for the irrigation of agricultural land in the Asembagus coastal area (East Java, Indonesia), with harmful consequences for rice yields. The river water has an unusual composition which is caused by seepage from the acidic Kawah Ijen crater lake into the river. This unique irrigation setting allows the study of soil acidification in situ. This paper assesses the effects of volcanogenically contaminated irrigation water on the chemical properties of the agricultural soils.The changes in soil properties were evaluated by comparing samples taken from the topsoil and sub-soil (1–3 m depth) from areas irrigated with acid water and areas irrigated with neutral water. The field survey thus resulted in four soil categories. Bulk soil composition, organic matter content, moisture content and particle size distribution were determined. Reactive phases were quantified with the selective extractions 1 M KCl, 0.1 M Na-pyrophosphate and 0.2 M acid ammonium oxalate (AAO).By comparing the four soil categories it is shown that the use of the naturally polluted irrigation water has had a large influence on the chemical composition of the topsoil. The composition of the soil solution has changed over the entire investigated soil profile. Furthermore the acid irrigation water has strongly modified the composition of the reactive phases, extracted as KCl, pyrophosphate, and AAO extractable elements, and also the bulk soil composition has been significantly modified. Overall this has resulted in the net dissolution of some elements and the net precipitation of others. The changes in the reactive phases and bulk soil composition are only apparent in the topsoil (0–20 cm) but not in the deeper soil.     

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.     

Usefulness of geological,mineralogical, chemical and chemometric analytical techniques in exploitation and profitability studies of iron mines and their associated elements

Due to the high number of variables involved in mine profitability studies, it is often very difficult to establish connections among them in order to provide a blend of market saleable quality products. In this sense, analytical chemistry together with chemometry are essential and indispensable disciplines to tackle these studies. The aim of this work was to demonstrate the utility of these disciplines to carry out optimization studies of iron mines. For this purpose, one of the most important iron mines of the Iberian Peninsula was chosen, sited in the mountain range of Sierra Menera, near the location of Ojos Negros (Teruel, Spain). Geological, mineralogical and chemical composition of 148 samples was analyzed, corresponding to different depths of three drill holes (named TE1, TE2 and TE3). In particular, aspects concerning to chemical composition are very important, since the mean contents of certain elements, such as phosphorus, sodium and potassium, should be restricted to the established limits to prevent that companies can drive back the raw material if they do not fulfil the necessary requirements. On the other hand, the large number of analysed samples drove us to use a statistical processing of the data. Among other aspects, it provides a way to find possible connections among a high number of variables and classify samples into compositional groups sharing similar composition, in order to limit the mineralised area and to obtain enough information about the amount of those chemical elements associated to iron ores. Data obtained from all these analytical techniques were in good agreement and provide a methodology that can be of wide interest applied to different geological studies.     

Ground motion estimation during the Kashmir earthquake of 8th October 2005

In this article, analytical methods have been used to estimate ground motion during the 8 October 2005, Kashmir earthquake. Peak ground acceleration (PGA) values at several stations in the epicentral region have been estimated by empirical analytical source mechanism models. As an alternate analysis, PGA estimates have also been obtained using the stochastic finite fault seismological model. The estimated PGAs are compared with that obtained from damage values. A PGA contour map in the near-source region is provided. It is found that very near to the epicenter, PGA would have reached more than 1 g. It is demonstrated that empirical analytical models can be effectively used to estimate ground motion due to rupture of active faults.     

Integration and magnitude homogenization of the Egyptian earthquake catalogue

The aim of the present work is to compile and update a catalogue of the instrumentally recorded earthquakes in Egypt, with uniform and homogeneous source parameters as required for the analysis of seismicity and seismic hazard assessment. This in turn requires a detailed analysis and comparison of the properties of different available sources, including the distribution of events with time, the magnitude completeness, and the scaling relations between different kinds of magnitude reported by different agencies. The observational data cover the time interval 1900–2004 and an area between 22°–33.5° N and 25°–36° E. The linear regressions between various magnitude types have been evaluated for different magnitude ranges. Using the best linear relationship determined for each available pair of magnitudes, as well as those identified between the magnitudes and the seismic moment, we convert the different magnitude types into moment magnitudes M _W, through a multi-step conversion process. Analysis of the catalogue completeness, based on the M _W thus estimated, allows us to identify two different time intervals with homogeneous properties. The first one (1900–1984) appears to be complete for M _W ≥ 4.5, while the second one (1985–2004) can be considered complete for magnitudes M _W ≥ 3.     

Mineralogy and geochemistry of El Dorado epithermal gold deposit, El Sauce district, central-northern Chile

Summary The El Dorado Au-Cu deposit is located in an extensive intra-caldera zone of hydrothermal alteration affecting Upper Cretaceous andesites of the Los Elquinos Formation at La Serena (≈ 29°47′S Lat., 70°43′W Long., Chile). Quartz-sulfide veins of economic potential are hosted by N25W and N20E fault structures associated with quartz-illite alteration (+supergene kaolinite). The main ore minerals in the deposit are pyrite, chalcopyrite ± fahlore (As/(As + Sb): 0.06−0.98), with electrum, sphalerite, galena, bournonite-seligmanite (As/(As + Sb): 0.21−0.31), marcasite, pyrrhotite being accessory phases. Electrum, with an Ag content between 32 and 37 at.%, occurs interstitial to pyrite aggregates or along pyrite fractures. Pyrite commonly exhibits chemical zonation with some zones up to 1.96 at.% As. Electron probe microanalyses of pyrite indicate that As-rich zones do not exhibit detectable Au values. Fluid inclusion microthermometry shows homogenization temperatures between 130 and 352 °C and salinities between 1.6 and 6.9 wt.% NaCl eq. Isotope data for quartz, ankerite and phyllosilicates and estimated temperatures show that δ¹⁸O and δD for the hydrothermal fluids were between 3 and 10‰ and between −95 and −75‰, respectively. These results suggest the mineralizing fluids were a mixture of meteoric and magmatic waters. An epithermal intermediate-sulfidation model is proposed for the formation of the El Dorado deposit. Author’s present address: J. Carrillo-Rosúa, Dpto. de Didáctica de las Ciencias Experimentales, Universidad de Granada, Campus de Cartuja, 18071, Granada, Spain     

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.     

Implementation of reconstructed geomorphologic units in landslide susceptibility mapping: the Melen Gorge (NW Turkey)

In the international literature, although considerable amount of publications on the landslide susceptibility mapping exist, geomorphology as a conditioning factor is still used in limited number of studies. Considering this factor, the purpose of this article paper is to implement the geomorphologic parameters derived by reconstructed topography in landslide susceptibility mapping. According to the method employed in this study, terrain is generalized by the contours passed through the convex slopes of the valleys that were formed by fluvial erosion. Therefore, slope conditions before landsliding can be obtained. The reconstructed morphometric and geomorphologic units are taken into account as a conditioning parameter when assessing landslide susceptibility. Two different data, one of which is obtained from the reconstructed DEM, have been employed to produce two landslide susceptibility maps. The binary logistic regression is used to develop landslide susceptibility maps for the Melen Gorge in the Northwestern part of Turkey. Due to the high correct classification percentages and spatial effectiveness of the maps, the landslide susceptibility map comprised the reconstructed morphometric parameters exhibits a better performance than the other. Five different datasets are selected randomly to apply proper sampling strategy for training. As a consequence of the analyses, the most proper outcomes are obtained from the dataset of the reconstructed topographical parameters and geomorphologic units, and lithological variables that are implemented together. Correct classification percentage and root mean square error (RMSE) values of the validation dataset are calculated as 86.28% and 0.35, respectively. Prediction capacity of the different datasets reveal that the landslide susceptibility map obtained from the reconstructed parameters has a higher prediction capacity than the other. Moreover, the landslide susceptibility map obtained from the reconstructed parameters produces logical results.