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.     

Biaxial test simulations using a packing of polygonal particles

The mechanical response of cohesionless granular materials under monotonic loading is studied by performing molecular dynamic simulations. The diversity of shapes of soil grains is modelled by using randomly generated convex polygons as granular particles. Results of the biaxial test obtained for dense and loose media show that samples achieve the same void ratio at large strains independent of their initial density state. This limit state resembles the so‐called critical state of soil mechanics, except for some stress fluctuations, which remain for large deformations. These fluctuations are studied at the micro‐mechanical level, by following the evolution of the co‐ordination number, force chains and the fraction of the sliding contacts of the sample. Copyright © 2007 John Wiley & Sons, Ltd.     

Methodology for preliminary assessment of Natech risk in urban areas

Concern for natural hazard-triggered technological disasters (Natech disasters) in densely populated and industrialized areas is growing. Residents living in urban areas subject to high natural hazard risk are often unaware of the potential for secondary disasters such as hazardous materials releases from neighboring industrial facilities, chemical storage warehouses or other establishments housing hazardous materials. Lessons from previous disasters, such as the Natech disaster during the Kocaeli earthquake in Turkey in 1999 call for the need to manage low frequency/high consequence events, particularly in today’s densely populated areas. However, there is little guidance available on how local governments and communities can assess Natech risk. To add to the problem, local governments often do not have the human or economic resources or expertise to carry out detailed risk assessments. In this article, we propose a methodology for preliminary assessment of Natech risk in urban areas. The proposed methodology is intended for use by local government officials in consultation with the public. The methodology considers possible interactions between the various systems in the urban environment: the physical infrastructure (e.g., industrial plants, lifeline systems, critical facilities), the community (e.g., population exposed), the natural environment (e.g., delicate ecosystems, river basins), and the risk and emergency management systems (e.g., structural and nonstructural measures). Factors related to vulnerability and hazard are analyzed and qualitative measures are recommended. Data from hazardous materials releases during the Kocaeli, Turkey earthquake of August 17, 1999 are used as a case study to demonstrate the applicability of the methodology. Limitations of the proposed methodology are discussed as well as future research needs.

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     

Mineralogy and geochemistry of platinum-group elements in the Aguablanca Ni-Cu deposit (SW Spain)

Summary The Aguablanca Ni-Cu-(PGE) magmatic sulphide deposit is associated with a magmatic breccia located in the northern part of the Aguablanca gabbro (SW, Iberia). Three types of ores are present: semi-massive, disseminated, and chalcopyrite-rich veined ore. The principal ore minerals are pyrrhotite, pentlandite and chalcopyrite. A relatively abundant platinum-group mineral (PGM) assemblage is present and includes merenskyite, melonite, michenerite, moncheite and sperrylite. Moreover, concentrations of base and precious metals and micro-PIXE analyses were obtained for the three ore-types. The mineralogy and the mantle-normalised chalcophile element profiles strongly suggest that semi-massive ore represents mss crystallisation, whereas the disseminated ore represents an unfractionated sulphide liquid and the chalcopyrite-rich veined ore a Cu-rich sulphide liquid. Palladium-bearing minerals occur commonly enclosed within sulphides, indicating a magmatic origin rather than hydrothermal. The occurrences and the composition of these minerals suggest that Pd was initially dissolved in the sulphides and subsequently exsolved at low temperatures to form bismutotellurides. Negative Pt and Au anomalies in the mantle-normalised chalcophile element profiles, a lack of Cu-S correlation and textural observations (such as sperrylite losing its euhedral shape when in contact with altered minerals) suggest partial remobilisation of Pt, Au and Cu by postmagmatic hydrothermal fluids after the sulphide crystallisation. Authors’ addresses: R. Pi?a, L. Ortega, R. Lunar, Departamento de Cristalografía y Mineralogía, Facultad de Geología, Universidad Complutense de Madrid, ES-28040 Madrid, Spain; F. Gervilla, Facultad de Ciencias, Instituto Andaluz de Ciencias de la Tierra, Universidad de Granada-CSIC, Avda. Fuentenueva, s/n, ES-18002 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.     

Flexural isostatic response of the Alps to increased Quaternary erosion recorded by foreland basin remnants,SE France

We test the hypothesis that flexural isostatic compensation of the mass removed by enhanced Quaternary erosion is responsible for uplift of the Western European Alps and their forelands. We use two well‐preserved and well‐dated (1.8 Ma) abandonment surfaces of foreland basin remnants in SE France (the Chambaran and Valensole plateaux) as passive benchmarks for tilting of the foreland. Estimating their initial slope from morphometric scaling relationships, we determine bulk post‐depositional tilting of 0.5–0.8% for these surfaces. The calculated isostatic response of the Alpine lithosphere to erosional unloading, using the method recently proposed by Champagnac et al. [Geology 35 (2007) 195–198] , yields a predicted tilting of 0.3–0.4% in the considered areas, explaining approximately half of the determined post‐depositional tilting. Such long‐term deformation being insensitive to cyclic loading/unloading because of glaciations, we suspect the other half to be related to as yet undetermined long‐wavelength and long‐lived tectonic process(es).     

Ermittlung der QNH-werte aus dem luftdruckfeld der synoptischen bodenkarte

Резюме Дабление воздуха, переснитанное иа уровень моря по стандартной атмосфере в Q-коде обозчачается через QNH. Давление воздуха пересчитанное на уровень моря по высотной барометрической формуле обознаеается через QFF. Для целей авиационной службы погоды должны быть известны значения QNH, однако Зе барическое поле на синоптических картах выражается через QFF. С помощью рис. 1 для соответствующей температуры воздуха на станцин и ее высоты н. у. м. можно определить разность значений QFF—QNH при давленин QFF=1000мб. Далее по табл. 2 можно определить поправку для каждого значения QFF отличного от значения QFF при ином давлении, чем 1000мб путем умножения табулированного значения ва разность QFF—1000мб и его алгебраического сложения со значением, полученным по рис. 1.