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.     

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.     

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.     

Boron incorporation into mullite

Summary Boron-doped mullites were synthesized using aluminium nitrate-nonahydrate, tetraethoxysilane and boric acid in a sol–gel process with subsequent annealing at 950 and 1300 °C for five hours. Two different bulk compositions with constant Al₂O₃ contents (60 and 70 mol%, respectively) and varying SiO₂ plus B₂O₃ contents were investigated. X-ray powder diffraction analyses yielded a linear decrease of the lattice parameters with increasing bulk B₂O₃ content, which was interpreted as to be due to boron incorporation. Related to the increasing boron content, corresponding infrared spectra revealed a slight and continuous shift for most of the absorption bands. These data show that mullite is able to incorporate large amounts of boron into its structure (up to about 20 mol% B₂O₃ depending on the bulk composition of the starting materials). Infrared analyses suggest that boron is incorporated into the mullite structure in form of planar three-fold coordinated BO₃ groups. Author’s address: A. Beran, Institut für Mineralogie und Kristallographie, Universit?t Wien-Geozentrum, Althanstra?e 14, 1090 Wien, Austria     

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.     

Nonlinear Analysis of Local Site Effects on Seismic Ground Response in the Bam Earthquake

The influence of local geologic and soil conditions on the intensity of ground shaking is addressed in this study. The amplification of the ground motion due to local site effects resulted in severe damage to dwellings in the Bam area during the 2003 Bam Earthquake. A unique set of strong motion acceleration recordings was obtained at the Bam accelerograph station. Although the highest peak ground acceleration recorded was the vertical component (nearly 1 g), the longitudinal component (fault-parallel motion) clearly had the largest maximum velocity as well as maximum ground displacement. Subsurface geotechnical and geophysical (down-hole) data in two different sites have been obtained and used to estimate the local site condition on earthquake ground motion in the area. The ground response analyses have been conducted considering the nonlinear behavior of the soil deposits using both equivalent linear and nonlinear approaches. The fully nonlinear method embodied in FLAC was used to evaluate the nonlinear soil properties on earthquake wave propagation through the soil layer, and compare with the response from the equivalent linear approach. It is shown that thick alluvium deposits amplified the ground motion and resulted in significant damage in residential buildings in the earthquake stricken region. The comparison of results indicated similar response spectra of the motions for both equivalent and nonlinear analyses, showing peaks in the period range of 0.3–1.5 s. However, the amplification levels of nonlinear analysis were less than the equivalent linear method especially in long periods. The observed response spectra are shown to be above the NEHRP building code design requirements, especially at high frequencies.     

AERIAL SURVEY: COLONY OF HONG KONG

This paper describes the planning and implementation of a very extensive aerial survey contract undertaken by Hunting Surveys on behalf of the Hong Kong Government. From photography carried out during January and February 1963, photogrammetric plotting at 1 : 600 and 1 : 1200 proceeded steadily over a period of 5 years to the point where all but the highest land in the Colony had been surveyed at one or other of these two large scales. The problems arising from the rather unique terrain conditions and the urban and industrial development of Hong Kong are dealt with mainly from the photographic and photogrammetric aspect but reference is made to the contribution of the Crown Lands and Survey Office of the Public Works Department of Hong Kong, who provided all ground control and undertook field verification.