Experimental estimation of velocities and anisotropy of a series of Swedish crystalline rocks and ores

To provide a guide for future deep (<1.5 km) seismic mineral exploration and to better understand the nature of reflections imaged by surface reflection seismic data in two mining camps and a carbonatite complex of Sweden, more than 50 rock and ore samples were collected and measured for their seismic velocities. The samples are geographically from the northern and central parts of Sweden, ranging from metallic ore deposits, meta‐volcanic and meta‐intrusive rocks to deformed and metamorphosed rocks. First, ultrasonic measurements of P‐ and S‐wave velocities at both atmospheric and elevated pressures, using 0.5 MHz P‐ and S‐wave transducers were conducted. The ultrasonic measurements suggest that most of the measured velocities show positive correlation with the density of the samples with an exception of a massive sulphide ore sample that shows significant low P‐ and S‐wave velocities. The low P‐ and S‐wave velocities are attributed to the mineral texture of the sample and partly lower pyrite content in comparison with a similar type sample obtained from Norway, which shows significantly higher P‐ and S‐wave velocities. Later, an iron ore sample from the central part of Sweden was measured using a low‐frequency (0.1–50 Hz) apparatus to provide comparison with the ultrasonic velocity measurements. The low‐frequency measurements indicate that the iron ore sample has minimal dispersion and attenuation. The iron ore sample shows the highest acoustic impedance among our samples suggesting that these deposits are favourable targets for seismic methods. This is further demonstrated by a real seismic section acquired over an iron ore mine in the central part of Sweden. Finally, a laser‐interferometer device was used to analyse elastic anisotropy of five rock samples taken from a major deformation zone in order to provide insights into the nature of reflections observed from the deformation zone. Up to 10% velocity‐anisotropy is estimated and demonstrated to be present for the samples taken from the deformation zone using the laser‐interferometery measurements. However, the origin of the reflections from the major deformation zone is attributed to a combination of anisotropy and amphibolite lenses within the deformation zone.     

On correct application of one-step inversion of gravity data

One of the main problems on the numerical solution of integral equations is the resolution of input data. Among the integral equations used in geodesy we have the “onestep inversion” based on the first derivative of the Poisson integral, which transforms gravity values on the Earth’s surface to the gravity potential on the reference ellipsoid. In this study, it is shown that the required spatial resolution of the input gravity data on the Earth’s surface for correct one-step inversion depends on the height of the computational region, the fact that if overlooked can cause totally wrong results. Consequently the following two major questions are posed: (i) How could one know whether the spatial resolution of the input gravity data for correct one-step inversion is sufficient? (ii) What should be done if the spatial resolution is not sufficient? By studying the behaviour of the integral kernel, an algorithm is presented which enables an appropriate answer to the former question. In order to address the latter question, a method is proposed to modify the integral kernel which overcomes the adverse effect of insufficient spatial resolution of the input gravity data. Our answers, which possess the novelty of the study, are numerically verified by means of real and simulated gravity data. The numerical results approve the efficiency of the proposed method in solving the problem of insufficient spatial resolution of the input gravity data for correct one-step inversion.     

Spatial distribution of cadmium and lead in the sediments of the western Anzali wetlands on the coast of the Caspian Sea (Iran)

Spatial distribution patterns of total cadmium (Cd) and lead (Pb), their bioavailable fractions and total organic matter in sediment from Anzali wetlands are provided. Total sediment Pb was higher than Cd (34.95 versus 0.024 μg/g dry weight). The geoaccumulation index indicated that the sediment was “uncontaminated”, but some stations were categorized as “unpolluted” to “moderately polluted”. Less than 0.01 of Pb existed in exchangeable and carbonate fractions. The sum of exchangeable and carbonate-bound fractions of Cd was 42%, suggesting that Cd poses high risk to the aquatic ecosystems. Total Cd and Pb exhibited positive relationships with total organic matter. Considering spatial distribution maps of total and bioavailable fractions of metals suggested that high concentrations of metals does not necessarily indicate high bioavailable fraction. The methodologies we used in this study can be in more effective management of aquatic ecosystems, as well as ecological risk assessment of metals, and remediation programs.     

Application of first-arrival tomography to characterize a quick clay landslide site in Southwest Sweden

First-arrival traveltime tomography was applied to high-resolution seismic data acquired over a known quick-clay landslide scar near the Göta River in southwest Sweden in order to reveal the geometry and physical properties of clay-related normally consolidated sediments. Investigated area proved to be a challenging environment for tomographic imaging because of large P-wave velocity variations, ranging from 500 to 6000 m/s, and relatively steeply-dipping bedrock. Despite these challenges, P-wave velocity models were obtained down to ca. 150 m for two key 2D seismic profiles (each about 500-m long) intersecting over the landslide scar. The models portrait the sandwich-like structure of marine clays and coarse-grained consolidated sediments, but the estimated resolution (20 m) is too small to distinguish thin layers within this structure. Modelled velocity structures match well the results of reflection seismic processing and resistivity tomography available along the same profiles.     

A study of differentiation pattern and rare earth elements migration in geochemical and hydrogeochemical environments of Airekan and Cheshmeh Shotori areas (Central Iran)

The Airekan and Cheshmeh Shotori areas are located about 60 km northeast of Khour, in Isfahan province from Central Iran. Research on characteristics and rare earth elements (REE) pattern in hydrogeochemical environments of these areas suggests the same origin for the elements dissolved in groundwater in these areas. Investigation of migration pattern of REE in hydrogeochemical environments shows that the migration and transportation of REE has occurred through chloride complexes. REEs, leached by water/rock interaction from the Airekan granite, are transported by groundwater and then precipitated in the Cheshmeh Shotori area. Study of the Cheshmeh Shotori sediments shows the presence of a sequence of red oxidized and dark layers. Geochemical characteristics of these sediments reveal that their REE characteristics are mainly inherited from the Airekan granite. Changes in the REE pattern of these sediments with depth show that changes in oxidation and reduction process have not played a significant role in controlling their behavior. It is crucial to note that adsorption of REEs dissolved in water by hydrosilicate increases these elements in depth. The REE behavior shows water/rock interaction between the granitic rocks and groundwater as the main factor of solution, migration and precipitation of REEs in the Cheshmeh Shotori area.     

Integration of ASTER and landsat TM remote sensing data for chromite prospecting and lithological mapping in Neyriz ophiolite zone,south Iran

Chromite deposits in Iran are located in the ophiolite complexes, which have mostly podiform types and irregular in their settings. Exploration for podiform chromite deposits associated with ophiolite complexes has been a challenge for the prospectors due to tectonic disturbance and their distribution patterns. Most of Iranian ophiolitic zones are located in mountainous and inaccessible regions. Remote sensing approach could be applicable tool for choromite prospecting in Iranian ophiolitic zones with intensely rugged topography, where systematic sampling and conventional geological mapping are limited. In this study, Landsat Thematic Mapper (TM) and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite data were used for chromite prospecting and lithological mapping in the Neyriz ophiolitic zone in the south of Iran. Image transformation techniques, namely decorrelation stretch, band ratio and principal component analysis (PCA) were applied to Landsat TM and ASTER data sets for lithological mapping at regional scale. The RGB decorrelated image of Landsat TM spectral bands 7, 5, and 4, and the principal components PC1, PC2 and PC3 image of ASTER SWIR spectral bands efficiently showed the occurrence of major lithological units in the study area at regional scale. The band ratios of 5/3, 5/1, 7/5 applied on ASTER VNIR‐SWIR bands were very useful for discriminating most of rock units in the study area and delineation of the transition zone and mantle harzburgite in the Neyriz ophiolitic complex. Spectral Angle Mapper (SAM) technique was implemented to ASTER VNIR‐SWIR spectral bands for detecting minerals of rock units and especially delineation of the transition zone and mantle harzburgite as potential zones with high chromite mineralization in the Neyriz ophiolitic complex. The integration of information extracted from the image processing algorithms used in this study mapped most of lithological units of the Neyriz ophiolitic complex and identified potential areas of high chromite mineralization (transition zone and mantle harzburgite) for chromite prospecting targets in the future. Furthermore, image processing results were verified by comprehensive fieldwork and laboratory analysis in the study area. Accordingly, result of this investigation indicate that the integration of information extracted from the image processing algorithms using Landsat TM and ASTER data sets could be broadly applicable tool for chromite prospecting and lithological mapping in mountainous and inaccessible regions such Iranian ophiolitic zones.     

Flood frequency analysis based on simulated peak discharges

Flood frequency approaches vary from statistical methods, directly applied on the observed annual maximum flood series, to adopting rainfall–runoff simulation models that transform design rainfalls to flood discharges. Reliance on statistical flood frequency analysis depends on several factors such as the selected probability distribution function, estimation of the function parameters, possible outliers, and length of the observed flood series. Through adopting the simulation approach in this paper, watershed-average rainfalls of various occurrence probabilities were transformed into the corresponding peak discharges using a calibrated hydrological model. A Monte Carlo scheme was employed to consider the uncertainties involved in rainfall spatial patterns and antecedent soil moisture condition (AMC). For any given rainfall depth, realizations of rainfall spatial distribution and AMC conditions were entered as inputs to the model. Then, floods of different return periods were simulated by transforming rainfall to runoff. The approach was applied to Tangrah watershed in northeastern Iran. It was deduced that the spatial rainfall distribution and the AMCs exerted a varying influence on the peak discharge of different return periods. Comparing the results of the simulation approach with those of the statistical frequency analysis revealed that, for a given return period, flood quantiles based on the observed series were greater than the corresponding simulated discharges. It is also worthy to note that existence of outliers and the selection of the statistical distribution function has a major effect in increasing the differences between the results of the two approaches.     

Strong ground motion record selection for the reliable prediction of the mean seismic collapse capacity of a structure group

How to select a limited number of strong ground motion records (SGMRs) is an important challenge for the seismic collapse capacity assessment of structures. The collapse capacity is considered as the ground motion intensity measure corresponding to the drift‐related dynamic instability in the structural system. The goal of this paper is to select, from a general set of SGMRs, a small number of subsets such that each can be used for the reliable prediction of the mean collapse capacity of a particular group of structures, i.e. of single degree‐of‐freedom systems with a typical behaviour range. In order to achieve this goal, multivariate statistical analysis is first applied, to determine what degree of similarity exists between each selected small subset and the general set of SGMRs. Principal Component analysis is applied to identify the best way to group structures, resulting in a minimum number of SGMRs in a proposed subset. The structures were classified into six groups, and for each group a subset of eight SGMRs has been proposed. The methodology has been validated by analysing a first‐mode‐dominated three‐storey‐reinforced concrete structure by means of the proposed subsets, as well as the general set of SGMRs. The results of this analysis show that the mean seismic collapse capacity can be predicted by the proposed subsets with less dispersion than by the recently developed improved approach, which is based on scaling the response spectra of the records to match the conditional mean spectrum. Copyright © 2010 John Wiley & Sons, Ltd.