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.     

Geological research at the Natural History Museum of Moscow University, 1759-1930 （in Russian）

In 2006, the Moscow publishing house Nauka （Science） published a 246-page book by Zoya A. Bessudnova entitled Geological Research at the Natural History. Museum of Moscow University, 1759-1930 as No. 32 in its series ＂Sketches on the History of Geological Knowledge＇. The author is a geologist and Senior Researcher in the Department for the History of Geology of the Vernadsky State Geological Museum of the Russian Academy of Sciences, and a member of the International Commission on the History of Geological Sci- ences （INHIGEO）.     

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     

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.     

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.     

Analysis of the 2003 Varunawat Landslide,Uttarkashi, India using Earth Observation data

Mass movements such as landslides in mountainous terrains are natural degradation processes and one of the most important landscape-building factors. Varunawat Parbat overlooking Uttarkashi town witnessed a series of landslides on 23 September 2003 and the debris slides and rock falls continued for 2 weeks. This landslide complex was triggered due to the incessant rainfall prior to the event, and its occurrence led to the blockage of the pilgrim route to Gangotri (source of the Ganges river) and evacuation of thousands of people to safer places. Though there was no loss of lives due to timely evacuation, heavy losses to the property were reported. High-resolution stereoscopic earth observation data were acquired after the incidence to study the landslide in detail with emphasis on the cause of the landslide and mode of failure. Areas along the road and below the Varunawat foothill region are mapped for landslide risk. It was found that the foothill region of the Varunawat Parbat was highly disturbed by man-made activities and houses are dangerously located below steep slopes. The potential zones for landslides along with the existing active and old landslides are mapped. These areas are critical and their treatment with priority is required in order to minimise further landslide occurrences.