Journal of Engineering Geology

---

Prediction of location of future earthquakes with event probability is useful in reduction of earthquake hazard. Determination of predicted locations has attracted more attention to design, seismic rehabilitation and reliability of structures in these sites. Many theories were proposed in the prediction of time of occurrence of earthquake. There is not a method for prediction time of future earthquakes. Many studies have been done in the prediction of magnitude of earthquakes, but there are not any investigations on prediction of earthquake hazard zonation. In this study, the locations that have probability of the event of future earthquake have been predicted by artificial neural networks in Qum and Semnan. Neural networks used in this study can extract to complicate properties of patterns by receipting the interval patterns. Furthermore, the map of earthquake hazard zonation has been drawn. Properties of occurred earthquake were collected since 1903. The most probable event of earthquake in Qum has been predicted 31.6% in center, and 28.9% in north of Semnan

---

Time history analyses as crucial means in many earthquake engineering applications are highly dependent to characteristics of the seismic excitation record so that the resulting responses may vary from case to case. Strong ground motion scaling is a known codified solution to reduce such a dependency and increase reliability of time history analyses. The well-known code practice may result in highly non-economic designs due to considerable error in the spectra scaled to match the target code spectrum. This problem is formulated here in an optimization framework with the scaling coefficients as the design variables. Harmony search as a recent meta-heuristic algorithm is utilized to solve the problem and is applied to the treated examples. Using a variety of target period ranges the scaling error is evaluated and studied after more unified via optimization. The effect of base structural period and interval variation on the scaling error is then studied in addition to considerable error decrease with respect to traditional code-based procedure. The results also show dependency of spectral matching error to the period-interval elongation/variation, the base-structural period and more error sensitivity for narrow-band resonance with the filtered records on softer soil types

---

Ground-Penetrating Radar (GPR) is a non-destructive and high-resolution geophysical method which uses high-frequency electromagnetic (EM) wave reflection off buried objects to detect them. In current research this method has been used to identify geometrical parameters of buried cylindrical targets such as tunnel structures. To achieve this aim, relationships between the geometrical parameters of cylindrical targets with the parameters of GPR hyperbolic response have been determined using two intelligent pattern recognition methods known as artificial neural network and template matching. To this goal GPR responses of synthetic cylindrical objects produced by 2D finite-difference method have been used as templates in the neural network and template matching algorithms. The structure of applied neural network has been designed based on extracting discriminant and unique features (eigenvalues and the norm of eigenvalues) from the GPR images and predicting all geometrical parameters of the targets, simultaneously. Also the template matching operation carried out using two diverse similarity approaches, spatial domain convolution and normalized cross correlation in 2D wave number domain. The results of the research show that both two employed intelligent methods can be applied for in situ, real-time, accurate and automatic interpretation of real GPR radargrams, however in general the neural network method has led to less error and better estimation than template matching to predict the geometrical parameters of the cylindrical tar

---

Extended Abstract
Introduction
Dimension stones market is considered as an important and profitable sector of mineral deposit business due to their share in national economic performance. There exist a number of technical reports highlighting a lack of rock quality control in the sequence of quarrying and dimension stones production procedures, which has lowered the production efficiency and consequently the profitability of this strategic mineral industry in Iran. The quality of dimension stones depends on several factors which fractures, joints, voids and fine beddings are the most important factors that down-grade the quality. Therefore, foremost the quality and desirability of the building stone must be precisely determined by sampling, compressive strength testing and preparing microscopic sections. All of the mentioned evaluation methods are destructive. Moreover, sampling and performing multiple tests on all parts of a quarry or on all quarried stone blocks, is not possible. Detection of fractures hidden into the dimension stone blocks is achievable using fast, low-cost, accurate and non-destructive ground-penetrating radar (GPR) method. GPR is a high-resolution geophysical method which uses electromagnetic waves with high-frequency in order to map structures and detect buried objects in subsurface without coring or any destruction of the medium.
 
 
Materials and methods
In current research, GPR method has been applied to evaluate the quality of quarried travertine blocks at Haji-Abad quarry complex in Mahallat district, Markazi province, before starting any processing operation. To achieve this goal, the 2-D GPR responses of synthetic models resembling cubic dimension stone blocks containing fine layering and discontinuities, were primarily simulated using a modified 2-D finite-difference forward modeling program in the frequency-domain coded in MATLAB. Among the variety of available numerical methods, the finite-difference time-domain (FDTD) method has paid more attention due to having the simple understanding of the concepts, flexibility, simulation and modeling of complex environments and the acceptability of its responses in the applied cases. In this research, the simulation has been implemented for both calcareous and dolomitic rocks (including travertine and marbles) and granites. In the study area, the GPR data acquisition was carried out using a GPR system equipped with shielded 250 MHz central frequency antenna, 0.5 m antenna distance and 2 cm sampling intervals by monostatic common-offset reflection profiling method. In order to process, analyze and interpretation of data, Ground Vision and Radexplorer software were employed. The most important pre-processing and processing operations applied to the data to provide the final sections, comprising time-zero correction, dewowing (removing very low frequency components from the data), DC shift removal, Butterworth filtering, running average, background removal and types of amplitude gain.
Results and Conclusion
The results of the forward modeling show that the GPR response of fine beddings interfaces and major discontinuities hidden in the volume of dimension stone blocks are clearly detectable. Interpretation of the actual radargrams taken from a real GPR case study (Haji-Abad quarry complex) after employing various B-scan pre-processing and filtering procedures, indicates that GPR method is highly capable to detect fine beddings and discontinuities in order to evaluate the quality of dimension stone before starting any quarrying process. Validation of the obtained results of the present research was carried out on one of the blocks with a predicted large oblique joint while the existence of the large joint was proven under the cutting saw in the stone processing plant. However, it should be noted that due to the existence of inherent heterogeneity encompassing fine beddings, in addition to noises from different sources and their associated multiple reflections in real radargrams, the response of shallow major discontinuities may mask the response of minor ones located underneath or deeper, so as a result may not be detectable with routing GPR radargrams../files/site1/files/121/Ahmadi_Abstract.pdf
 
Keywords: Dimension stone Blocks (cubes); dimension stones production; Ground Penetrating Radar (GPR); Forward modeling; Quality control; Haji-Abad mining complex in Mahallat
 

---

Extended Abstract
Paper pages (319-342)
Tunneling in soft grounds causes to changes in displacements and subsequently in-situ stresses around ground. These displacements may damage structural assets. Thus, estimating the magnitude and shape of settlement curve is necessary. There are several empirical and analytical methods for predicting settlement. For example, Peck’s empirical method is well known method for predicting settlement due to tunneling. Tunneling process is done by imposing volume loss in tunnel. Then, soil displacement is measured by using image processing technique and that data is fitted to Gaussian curve. By conducting tests in loose and dense sands, it is concluded that by increasing relative density of the soil, the magnitude of settlement decreases and the settlement trough width will be increased. Also soil volume loss is not the same as the tunnel volume loss.  
Introduction
Many researchers investigated settlement due to tunneling but there is a lack of research about the effect of relative density on settlement. Marshall et al. (2012) by conducting centrifuge tests in high density sandy soil, showed that settlement trough is affected by tunnel size, tunnel depth and tunnel volume loss. Zhou et al. (2014) by performing several tests in loose, medium and dense sand, examined the effect of relative density on settlement and showed that by decreasing the relative density the magnitude of settlement increases and settlement trough width will be decreased. In this paper by using 1g physical modeling (Figure 1) which is designed in Sahand University of Technology, the effect of relative density on settlement has been studied.
Material and methods
Simulation of tunnel volume loss is carried out by using two different diameter tubes as a shield and lining (Figure 2), while pulling out the larger tube volume loss is imposed. Also by changing tube diameter different volume losses have been applied. Measuring of soil displacements is achieved by image processing technique. For this purpose, different photos are taken from the whole process of the test by digital camera and by using Geo PIV, settlement of ground is determined. 
Results and discussions
Experiments were conducted in loose and dense silica sands and the measured data have been fitted to Gaussian curve. The result showed that Peck equation fitted well to surface and sub-surface settlement data. As shown in Figures 3 and 4, contour of displacement curve versus normalized tunnel depth and distance versus normalized tunnel diameter indicate that in dense sands most of the displacement occurred in the region which placed in distance of 1.25 times tunnel diameter and in loose sands in the region of 0.6 times tunnel diameter. Thus, settlement trough width in loose sands is narrower. Also by measuring soil volume loss in loose and dense sands at different levels (Figure5) it is concluded that in loose sands due to less dilation, more volume loss is transferred to higher levels.
Conclusion
The following main conclusion can be drawn:
1. Gaussian curve predicts well surface and subsurface transverse settlements but selection of its parameters requires more accuracy that may result in inaccurate prediction.
2. Settlement curve in loose sands is narrower than dense sands.
3. Displacement and soil volume loss in loose sand are more than dense sand../files/site1/files/121/Emami_Abstract.pdf
 

---

Drilling and cutting stones as types of the engineering operations have encountered a lot of extensive and determining applications in different technical and engineering aspects of the mining. Estimating the drillability and cutability of stones by using drilling equipment and diamond wire saw have important roles in estimating the expenses and also designing mines. In this article some samples of carbonate ornamental stones from different mines in Iran have been studied in order to estimate and predict the drilling and also cutability rate.
In order to evaluate the effect of the textural specifications on the rate of drilling and cutability, first a picture was provided from the thin microscopic surface of every stone sample and then the area, perimeter, diameter the longest diagonal and the shortest diagonal of the grains in the sections were determined and the other textural specifications were also determined through using mathematical relations and equations. After that the relationship between the abovementioned parameters with the drilling and cutting rate were determined by using univariate fitting. And finally to achieve more correlation coefficient multivariate fitting was applied for the data. Among the textural specifications affecting the drilling rate textural coefficient, the diameter of the grain, dequi, the ratio of the grain condition and the index of grain size homogeneity had a significant relationship with the drillability rate and also among those affecting the cutting rate, textural coefficient, the diameter of the grain, dequi, density, shape factor, index of interlocking, and the index of grain size homogeneity had significant relationships with the cutting rate and at the end the final equation to predicate the drillability and cutability was produced for these parameters../files/site1/files/121/AleeiAbstract.pdf

---

Introduction
     Ground-penetrating radar (GPR) method is a pretty new, non-destructive and high-resolution geophysical method that is widely used to identify the thickness of snow and ice layers and glaciers bed, because snow and ice are transparent for electromagnetic (EM) waves. Therefore, this method has been used to determine the thickness and basement topography of Alam-kooh glacier. In this research, only the GPR acquired data using unshielded antenna with central frequency of 25 MHz along one line in Alam-kuh glacier, Kelardasht- Mazandaran, have been processed and interpreted. The GPR data acquisition has been done by using common offset method, and transmitter-receiver separation of 6 meters. The final real radargram related to one of the surveyed GPR profiles in this region has been prepared after applying various processing operations containing signal saturation correction, amplitude gain, f-k migration filtering and static (topography) correction on the raw data. After applying processing sequences on the acquired data, the EM waves reflection off the interfaces of different layers including the reflections from the glacier basement have been detected, and by assigning a suitable EM wave velocity in the ice (0.16 m/ns), the thickness of 50 m for the ice layer laid under the survey line has been estimated. Also, in present research, forward and inverse modeling of GPR data have been performed to employ for snow, ice and glaciological investigations in the AlamKooh region of Mazandaran. To achieve this goal, GPR response of synthetic model corresponding to the real radargram was simulated first, by 2-D finite-difference time domain (FDTD) method. Afterward the inversion method by solving an optimization problem was employed to validate the interpretation of real GPR data.
Methodology and Approaches
     Based on the nature, physical and geometric properties of the subsurface target in the field data, their synthetic model have been built and their two-dimensional GPR responses forward modeling using ReflexW software and finite difference algorithms improved in the frequency domain, have been obtained. Also, it has been used an effective algorithm, coded in GUI environment of MATLAB programming software and as a result, a reliable and accurate inverse modeling has been carried out. In the present study, to simulate the behavior of the propagation of EM waves in GPR method, two-dimensional finite difference method has been used. The main advantage of this method is its comparative simplicity of the concept, high accuracy and simple implementation for complex and arbitrary models as well as easily adjusting the antenna when applied. In this study, acquisition of GPR field data and synthetic data modeling have been carried out in TM mode. The radargrams of the GPR data have been demonstrated using ReflexW software after performing necessary processing sequences.
Results and Conclusions
     The obtained results reveal that moraine materials covering the surface of the area are mainly fine-grained granite. The bed-rock or basement in the area is also granite. The polarity representing ice-bed rock is clearly seen on the GPR profiles. The topography of the glacier basement has successfully been detectable using just by GPR method. The electrical resistive nature of the glacier has caused large penetration depth of GPR pulses in this research. Furthermore, the results of the research for presented profiles on the basis of forward and inverse modeling output of GPR data in comparison with real GPR radargrams in the region validated the accuracy of GPR investigations in the area. Although with a quick glance, the error obtained by the inverse modeling for real GPR data seems unexpected and unacceptable, absolutely the high rate of error depends on many factors influencing on the real earth models containing various limitations existing in all forward modeling algorithms and software packages, impossibility of making forward modeling exactly according to the real models (due to the complex nature of the ground), taking into account the homogeneity and uniform host environment and targets in the modeling process unlike the diversity, the presence of different types of noises and so on. Therefore, making a controlled geophysical test site and trying performance of inverse modeling algorithm for field GPR data in this site, as well as determining the important physical parameters such as dielectric permittivity and electrical conductivity by experimental method through sampling from different depths for complex geological environments are suggested../files/site1/files/124/1ahmadpur%DA%86%DA%A9%DB%8C%D8%AF%D9%87.pdf

---

Introduction
Every year, numerous casualties and a large deal of financial losses are experienced due to earthquake events. The losses incurred by an earthquake vary depending on local site effect. Some well-known examples include the earthquake in Caracas 1967, Mexico city 1985, Kalamata 1986, Loma Prieta 1989, Roodbar 1990, Bam 2003, Jammu and Kashmir 2005, Sichuan 2008 and Haiti 2010. Therefore, in order to conquer drastic effects of an earthquake, one should evaluate urban districts in terms of the local site effect. Various methods are available for the evaluation of site effect. One of the most common methods includes ambient noise survey. Today, this approach is being used as fast, applicable, cost-effective method. Ambient seismic noise are feeble ground motions with displacement amplitudes of about 0.1–1 μm and that can be detected by seismograph with high magnification. Many investigations have been conducted to determine the nature of ambient noise. One of the possible sources of ambient noise can be human activity, such as traffic, industrial noises and nature activity, such as wind, ocean waves. The Babol city is one of the largest cities in the north of Iran (Mazandaran province). It lies on alluvium beds in the region presenting a high seismic potential. Therefore, comprehensive studies are necessary to introduce suitable solutions for minimizing earthquake damage and loss of life. For this reason, in Babol city, ambient noise survey has been performed at 60 stations and the obtained data were analyzed with Nakamura or H/V method (1989). The results were compared with local geological, geotechnical and seismic data to confirm their reliability for a seismo-stratigraphic.
Methodology and data collection
The analysis of ambient noise was initially proposed by Kanai and Takana (1961). Since then, many researchers have used ambient noise for site effect evaluation. As it is said before, one of the most popular techniques for estimation of site effects in the regions with low seismicity is ambient noise survey by Nakamura or H/V method (1989). Based on the literature review, the Nakamura method (1989) has been used in many places. Many theoretical and experimental studies show that, this method has the capability of estimation of fundamental frequency. Ambient noise survey was carried out at 60 stations in Babol city. Ambient noise was recorded using a velocity meter SARA. Two horizontal and a vertical components of ambient noise at each location are recorded for duration of 15 min with 100 samples per second. Because the environmental noise has an effect on ambient noise they are recoded between 10 p.m. to 6 a.m. The locations were determined by using GPS at the sites. The ambient noise survey in this study was made in compliance with the guidelines of SESAME (2004).
Results and discussion
The maximum and minimum values of fundamental frequency of the present research are 11.4 and 0.65 Hz, respectively. Also, the maximum and minimum values of amplitude of H/V peaks have been calculated as 3.71 and 1.19, respectively. The most significant point is that the fundamental frequency of the major part of Babol city is smaller than 1 Hz in agreement with the previous knowledge of the city geological setting. Another relevant point is the presence of some stations with very high (> 5 Hz) fundamental frequencies. In these cases, ambient noise recording and data analysis were repeated, but similar results were obtained. Considering the lack of sufficient geotechnical data (in some stations), the above phenomena need to be integrated with other methods. For this purpose, the data of electrical resistivity tomography (ERT) were used. The ERT showed that small regions of the north-west, west, and south areas have high resistance values probably related with the presence of hard deposits in the shallow subsoil. Babolrood river diversion in the west part and its return to the previous direction in the northern section is possibly due to the existence of these relatively hard deposits. By comparing these two tests, we observed that the ERT results correlate with the ambient data analysis. Therefore, we can conclude that the high-frequency peaks measured are reliable, but we need direct investigation to associate them to a specific shallow geological layer. To validate the results, fundamental frequencies obtained from ambient noise survey were compared with geotechnical data, numerical analysis and seismic data in the study area. A general review shows that the geotechnical data, equivalent-linear analysis results and seismic data have an acceptable conformity with the results of ambient noise survey.
Conclusion
The results show minimum and maximum fundamental frequencies 0.65 and 11.4 Hz, respectively. Assessment also reveals that the major parts of Babol city have the fundamental frequencies less than 1 Hz, which are in conformity with that of previous research. According to the results of seven cross sections, it can be concluded that fundamental frequency variations are in line with the geotechnical and geological data in the study area. It means that this method is the appropriate way to assess the local site effect in the Babol city. It is also observed that besides the soil layers, the soil stiffness and its shear wave velocity are effective factors in changing the fundamental frequency. Site frequencies were also estimated by preliminary 1-D site modeling using the equivalent-linear method. In general, a reasonable correspondence between the methods was obtained. Using seismic data, the HVSR of two strong ground motions have been calculated and the results have been compared with the nearest ambient noise recording station. Analyzing the spectral ratios demonstrates that the value of the fundamental frequency obtained by the H/V method (1.06 Hz) is very close to that of frequencies obtained by the seismic data (0.95 and 0.90 Hz)../files/site1/files/124/6rezaee%DA%86%DA%A9%DB%8C%D8%AF%D9%87.pdf

---

Introduction
Ground-penetrating radar (GPR) is a high-resolution geophysical method which uses electromagnetic waves with high-frequency in order to map structures and objects buried in subsurface without any destruction of the medium. In present research, choice of optimum parameters of real data acquisition for this method has been studied. The governed behavior on the GPR fields can be simulated by solving the Maxwell’s equations and the appropriate boundary conditions that form the basis of electromagnetic theory. Among the variety of available numerical methods, the finite-difference time-domain (FDTD) method has paid more attention due to having the simple understanding of the concepts, flexibility, simulation and modeling of complex environments and the acceptability of its responses in the applied cases. The purpose of this study is to identify what reasonable information can be obtained from field data under different environmental conditions and different survey parameters.
 
Materials and methods
To achieve the goal, first forward modeling of GPR data has been carried out for several synthetic models corresponding to common targets in subsurface installations, using 2-D finite-difference time-domain method by means of GPRMAX, ReflexW and Radexplorer softwares. The main purpose of the simulations is investigation of the effect of survey parameters such as spatial sampling intervals (trace interspacing) and temporal sampling frequency on the GPR response of targets with various physical and geometrical parameters. Also to select and design the most appropriate conditions and survey parameters for real GPR data, numerous field traverses were performed in Isfahan University of Technology campus over the pre-known buried cylindrical targets containing power cable, petro-gas pipe, water pipeline and waste water pipeline with diverse host media. In this operation due to having one monostatic GPR system equipped by shielded antenna with central frequency of 250 MHz, some of the survey parameters containing central frequency, antenna separation and antenna directivity are invariant. The most important investigated survey parameters are temporal sampling frequency, spatial sampling distance (trace intervals), time window and number of stacked traces.
 
Results and discussion
Regarding carried out investigations through field data acquisition, in only one case the GPR system failed to detect any understated targets which this mode is related to choice a sampling distance of 1 cm and a sampling frequency of 504 MHz. The sampling frequency of 504 MHz is just capable to detect the surface water pipeline (due to its low burial depth). Also only in three cases the GPR system is capable to detect all subsurface targets so that the first mode of the trace interval is 2 cm and the sampling frequency is 1954 MHz, whereas in the latter two, the trace interval is 1 cm and the sampling frequencies have been selected 1563 and 1954 MHz. At the end success or failure of the targets detection was investigated on the basis of selected survey parameters and the probability of successful target detection was determined depending on the temporal and spatial sampling frequency so that the maximum probability of target detection is regarding to temporal sampling frequency of 1954 MHz and trace interval of 1 cm. Regarding GPR field data acquisition, considering the relations between the central frequency of GPR measurement systems, the depth of penetration and resolution, the diversity of materials and various components of the host media of targets and their surface overburdens a range of dierse equipments with a variety of frequencies is needed, which all of them are not generally available.
 
Conclusion
As a general conclusion of this study, in order to reduce the risk in GPR data acquisition operation, optimal survey parameters are suggested as follows:
The sampling frequency should be about 7 to 8 times the central frequency of the employed system (should not be less than this value in order to avoid aliasing and on the other hand, due to reduction in the amount of data and thus the memory needed for storage and processing), trace interspacing equal to 1 cm (in order to detect all buried targets especially targets with small size), the number of stacked traces equal to 16 (to reduce the amount of computer memory required for processing and storing data) and time window according to the computational-empirical relation (1).
                                                                                                                                                                (1)
Where W is time window, D is the maximum depth and V is the minimum velocity.
The results of this research are not restricted to the investigated case, but in practice are applicable for cases with similar host environments, especially in urban areas (which application of non-destructive methods such as GPR is necessary)../files/site1/files/131/6Extended_Abstract(1).pdf

---

Introduction
This work intends to apply ASTER images to discriminate hydrothermal alteration zones in Kerman Cenozoiic Magmatic Belt (KCMB). Band ratio, principal component analysis, Crosta and color composite images are important methods to analyze satellite images. Previous researches showed that these techniques are not able to discriminate hydrothermal alteration zones and they usually detect vegetation covering as alteration zones. The reason is found in the spectral signature of vegetation and alteration minerals. It means that they present the same interaction when face with electromagnetic energy in different wavelengths. Hydroxyl-bearing minerals are the important products of hydrothermal alteration. Clays, which contain Al-OH- and Mg-OH-bearing minerals and hydroxides in alteration zones, are distinguished by absorption bands in the 2.1–2.4 µm range of ASTER data. Solving these problems is difficult when using standard image-processing techniques such as band rationing, principal component analysis, or spectral angle mapper. In recent years, several attempts were made to extract altered regions in the areas covered with vegetation. To overcome this problem, this research uses ASTER data by applying support vector machine (SVM) algorithmn. SVM is a new technique for data classification in remote sensing application. This paper aims to investigate the potential of SVM algorithm in mapping of hydrothermally altered areas. In many applications, SVM has been shown to provide higher performance than traditional learning machines and has been introduced as powerful tools for solving classification problems. The adopted dataset contains three ASTER scenes using SWIR and VNIR bands, covering the Meiduk porphyry copper deposit, Kader, Abdar and Iju occurrences located in Kerman Province, southeast Iran.
Material and methods
This work has been prepared on three ASTER level 1B scenes. Two scenes were acquired on 18^th April 2000 and another scenes on 15^th June 2007. These scenes were georeferenced by using an orthorectified ETM ⁺  image,  in  UTM projection and WGS-84 ellipsoid as a datum.  The first two data sets were corrected for Crosstalk. Atmospheric corrections were also performed by using Fast Line of Sight Atmospheric Analysis of Spectral Hypercubes (FLAASH). The data sets were then mosaicked.­­ Internal Average Relative Reflectance (IARR) correction was also applied. In this part, the training and test samples of the ASTER data are presented. The adopted image is a multispectral satellite image that contains 2204 training pixels which 516 pixels are related to arjillic zone, 1278 pixels are related to phyllic zone and 500 pixels are pertinent to propylitic zone (Fig. 1).

Fig. 1. Training pixels for learning SVM algorithm; Red pixels: arjillic; Green pixels: phyllic; Blue pixels: propylitic
Results and discussion
ASTER bands 4, 6, 7 and 8 were applied for determination of phyllic and arjilic zones and 9 bands of ASTER for propylitic alteration. In order to evaluate the developed algorirhm, confusion matrix was used and validation showed that discrimination of phylic and arjilic is not possible but propylitic zone could be identified by SVM. Also, the present research introduced a new error function, so called blind error, which is calculated using confusion matrix. Based on blind error, SVM did not classify 73.6 percent of the alteration pixels. But the remained pixels were classified with accuracy of 66.06%. Honarmand et al. (2011) and Mojedifar et al. (2013) studied the field samples of the present study area. Their studies showed that sericitization is the most widespread form of hydrothermal alteration at the Iju, Serenu, Chahfiroozeh, Meiduk, Parkam, Kader and Abdar porphyry copper deposits. Two types of phyllic alteration could be found in the study area including ferric-iron-rich and iron-oxide poor phyllic alteration. ASTER images were also analyzed by band rationing and principal component analysis (PCA) in order to compare their results with the SVM classified image. A comparison of the field data with altered areas mapped by PCA reveals errors in the classified map. Vegetation cover and sedimentary rocks are enhanced, which are erroneously identified as areas of alteration. The band ratio approach yields similar errors to those produced by the PCA method. These problems are less evident in the classified image obtained by SVM. The qualitative assessment of the accuracy of these methods indicates that SVM algorithm could be a reliable technique for alteration mapping, provided that the nature of the training areas is well known.
Conclusion
A comparison of the results obtained from traditional classification methods and support vector machine algorithm was performed in order to map hydrothermal alteration. Since the known occurrence of mineralization in the study area is consistent with the mapped distribution of hydrothermal alteration using SVM, this method is suggested to apply in exploring for hydrothermal alteration in other parts of the Iranian Cenozoic magmatic belt../files/site1/files/132/6Extended_Abstracts.pdf
 

---

Introduction
The nature of near-field earthquake records is very complicated and uncertain. Due to this complexity, the prediction of the real structural responses has become very difficult. Based on the analysis of the physical characteristics of near-field records, it is possible to use the simplified mathematical models. Near-field ground motions which are often associated with a progressive directional phenomenon due to their particular type of the causative fault, have much more destructive effects on the structures than the other quake tremors. The related research results show that under the influence of a strong near-field ground motion which contains forward directivity effects, the structural responses would be entered to a great nonlinear domain. On the other hand, due to the limited number of available near-field records, it is needed to prepare artificial acclerograms which can simulate the characteristics of the strong ground motions. Thus, it is possible to achieve a vast data base corresponding to wide range of powerful ground motions using mathematical wavelets. As a result, it provides a general overview of these types of artificial quake tremors and prepares an extended knowledge on the performance of structures in confronting these destructive movements.
Material and methods
The results obtained from the seismological studies on strong near-field records indicate that the most of these tremors contain large amounts of kinetic energy corresponding to the content of low frequency band. Additionally, by ignoring the high-frequency band the coherent velocity pulses can be detected with acceptable accuracy. In order to separate the high and low frequency bands, the empirical mode decomposition (EMD) method is used based on programming in MATLAB software. Various methods have been proposed for simulation of near-field records which most of them is based on using harmonic functions and the spectral assessment of the low frequency band of earthquake records.
In this regards, one of the best closed form evaluation has been performed by Mavroeidis and Papageorgiou (2003) which is to be formulated by making parametric changes to the so-called Gabor wavelet and replacing a simpler function instead of the Gaussian curve with a more efficient algebraic statement. Ghodrati Amiri et al. (2012) proposed another efficient formulation matched either of the benefits of Mavroeidis’s and Gabor wavelets. Both of the aforementioned models are based on the preparing of an efficient multi-statement parametric configuration of harmonic wavelets as noted above. In this study, in addition to calibrate the desired closed-formulations on the velocity pulses of the selected strong records, the accuracy of the notified simulation has also been investigated from the spectral and energy point of views.

Results and discussion
The band of high frequencies corresponding to the spectral content of strong near-field records can be ignored appropriately. This is because the major amount of the related kinetic energy is usually transmitted in the form of a low frequency pulse along with a number of high frequency spikes. Generally, these features are displayed over a relatively short time domain. In this study, the analytical attention to this subject is concentrated on the simulation of coherent multiple pulses via EMD method. The purpose of such simulation is to create a wide range of powerful and high-energy artificial motions. Moreover, due to the limited availability of natural near-field earthquake records, the proposed pulses can be used to evaluate the structural seismic performance.
Conclusion
Generally, strong near-field records contain a few consecutive pulses with different periods and spectral configurations. The essential effects of these pulses must not be ignored in conducting of nonlinear dynamic time history analyses. Obviously, the effects of these type of earthquake records on the seismic response of mid-to-high rise structures (with a large periodic range) will be significant. Furthermore, the probable occurrence of the resonant mode, may cause destructive effects on the seismic response of structural skeletons. The proposed pulses in this study were formulated through the EMD method as well as performing an analytical calibration process related to both bands of high and low frequencies. The spectral evaluations of the fitted mathematical closed-form pulses were accomplished for the selected earthquake records. The obtained results indicate a good analytical convergence and correlation with the physical parameters of the natural ground motions.


./files/site1/files/134/7.pdf

---

Introduction
 Analysis of time, location and magnitude of foreshocks and aftershocks has been one of the most important cases for experts in various scientific fields such as: seismology, structural engineering and crisis management, and other interrelated fields. Since this analysis and the result of studies on seismotectonic and cases of earthquakes help us identify the foreshocks and aftershocks with the goal of decreasing losses and nervious stress of the injured community in quake-stricken areas and skilled crisis management. The cause fault of earthquake plays the important role in foreshocks and aftershocks of the earthquake. So, study on fault behaviour is a suitable method for analyzing and routing the basic parameters of foreshocks and aftershocks. Also, foreshocks and aftershock are important parts of any earthquake in a seismic area. The analysis of the basic parameters of the foreshocks is one of the most practical researches for reducing the risk of earthquakes. The identification of behavioral pattern of foreshocks can help researchers detect the active fault conditions for the occurrence of earthquakes in different areas. The present study is concerned with the study of behavioral patterns earthquakes, foreshocks, and aftershock of Zohan earthquake. Experience of large aftershocks in different parts of the world indicates that, following earthquakes and depending on seismic-tectonic conditions, large aftershocks are likely to occur in the earthquake-effected zone, which will aggravate the damage caused by earthquakes (Omi et al., 2013). The main factor contributing to the worsening of damage caused by aftershocks is the performance of structures that are weakened but not destroyed by main earthquakes and are, thus, highly likely to be destroyed by large aftershocks (Saket and Fatemi Aghda, 2006).
Material and methods
The present paper makes use of data collected in a real earthquake and similar expriences in other earthquakes for presenting a practical pattern for predicting primary earthquake patterns, determining the location, magnitude, and time of aftershocks. The target of this case is decreasing the effects of earthquake. To this end, we used the results from studies on basic parameters of foreshocks and aftershocks of Zohan earthquake, and 2012 earthquake in South Khorasan province. The rationale for selecting the aforementioned studies is: location of event, the Zohan earthquake, had been identificated as an area with high risk for the occurrence of earthquakes, although there has been no wide-scale earthquake in this area in the last two decades. These conditions are important causes for more concentrated studies on this area because there is a high chance for wide-scale earthquakes striking this area.
Result and Discussion
In this part of research, we conduct a study on the location, magnitude and depth of foeshocks. Some of the world-wide research suggested that these data can help to predict the time of  mainshocks. Studies conducted on the variations of frequency in foreshocks can follow this goal.
In this paper, the available statistical data such as periodical variations of seismicity in the weeks leading up to the main shock can be used as a tool for estimating the approximate time of a future important earthquake. The weekly variations of seismicity before Zohan earthquake indicate a relative increase and then decrease within a 100 km radius around the epicenter of the main shock.
 
Table 1: Variations of frequency of foreshock based magnitude before Zohan earthquake

Week before main shock	Frequency of foreshock in the Radius of 100Km from main shock	Frequency of foreshock(with M>2.5) in the Radius of 100Km from main shock
6	0	0
5	1	1
4	1	0
3	2	0
2	5	3
1	2	0

 
Studies on numerous earthquakes in Iran and other regions in the world show that the distribution of aftershocks can be related to fault type or the direction of principal stress (Saket and Fatemi Aghda, 2006) and (King et al., 1994). Whereas maximum Coulomb stress change is related to maximum principal stress in earthquakes, the concentration of aftershocks can coincide with the direction of maximum principal stress (σ₁) of the causative fault in mainshock. Considering the direction of maximum principal stress and its adaptation to the scattering of aftershocks, the above hypothesis is confirmed.
Also studies on frequency changes and seismic quiescence of small aftershocks help us in predicting future aftershocks. The results the of presented research by Itawa (2008) on the World earthquake catalogue suggest that seismic quiescence theory is true for different regions of the world. Based on the results of the  study mentioned above, this case can be used as a tool for predicting large aftershocks in Zohan earthquake.

Fig 1. Adaptation of direction of maximum principal stress with scatering of the aftershocks of Zohan earthquake. a: direction of maximum principal stress (σ₁) of the causative fault in mainshock. b- scaterring of the aftershocks
Table 2: Seismic sequience versus magnitude of aftershocks

Row	Seismic Quiescence for aftershocks	Aftershock Magnitude	Data and Time of aftershocks
1	13	3.0	2012/12/05 17:21:03
2	36	3.4	2012/12/05 17:57:03
3	161	3.1	2012/12/05 20:38:09
4	3906	3.9	2012/12/08 13:44:19

In addation, frequency of aftershocks and certain time distance (seismic quiescence) between their can use precursors for detecting the time of large aftershocks. The relevant analysis in this study showed that methods such as: time series beside seismic quiescence can help in conducting a more accurate time forecast of large aftershocks.
Conclusion

• The results of this research suggest that we can identify some of the charactristics of the main shock by focusing on location, magnitude and depth of foeshocks.
• In Zohan earthquake, the direction of maximum principal stress is adpated to the scattering of aftershocks, and this case suggests that there is a specific relationship between them.
• The relevant analysis in this study showed that the methods such as: time series beside seismic quiescence can help conduct a more accurate time forecast of large aftershocks../files/site1/files/144/saket.pdf

---

In Yazd Darreh-Zereshk copper deposit geophysical data containing magnetic, resistivity and induced polarization have been surveyed and 25 boreholes have been drilled in the area. In the present research, inversion and processing of geophysical data as well as their qualitative and quantitative accordance with boreholes assay data have been carried out. To achieve the goal first, total magnetic intensity map after applying necessary filters and processing, was mapped to identify surface and deep expansion of anomalies on it. Drawing the anomaly profile of magnetic stations surveyed along 4 geoelectric profiles shows that most of the magnetic anomaly zones have high chargeability and low resistivity that indicates the qualitative compatibility of magnetic and geoelectric data, as a result increasing the probability of mineralization in the area. Afterward  on the basis of qualitative interpretation of geoelectrical sections, optimal locations of drilling on the each profile were proposed. Plotting mineral deposit cross-section along the geoelectrical profiles using the boreholes assay data, revealed that drilling of some boreholes located on the geophysical profiles haven’t been based on the results of geophysical operation, carried out without any right logic, purpose and design. In general, the qualitative accordance of the results of geoelectrical operation with the boreholes assay data showed a pretty good qualitative accordance. Also investigation of linear correlation coefficient value between inverted geophysical data and borehole assay in a specific same range after a same definite gridding and interpolation of their values, overall indicated a relatively good quantitative accordance (between 0.4 and 0.7)../files/site1/files/151/1.pdf

---

Chahar-Gonbad region of Kerman province is geologically located in the southern part of central Iran zone, dominantly in Uromieh-Dokhtar volcanic belt. In this region, many high potential prospects, specially Cu-Au mineralization, have been detected during large scale exploration and reconnaissance phases. In this paper, remote sensing and field geophysics were used for alteration mapping on the surface and ore body delineation on the subsurface, respectively. To this end, we used an ASTER satellite image and different maps were generated by spectral technics such as false color composites and spectral ratios. Results showed argillic (and phillic) alteration in Bab-Zangoeie area is surrounded by propylitic alteration, which could be a promising evidence for Cu mineralization. Integrating these results with previous exploration studies led to selecting target area selection for ground study and field geophysics. We used both induced polarization (IP) and resistivity (RS) methods as two powerful geoelectrical methods by a pole-dipole array along four profiles. After preprocessing analysis, forward and inverse models were constructed in 2D section and 3D overlay model of joint IP/RS anomalies were constructed. Based on the obtained results, the deposit in depth where we proposed drilling targets. Further drilling operation have proved the mineralization in our proposed targets../files/site1/files/152/%D9%86%D8%B8%D8%B1%DB%8C.pdf
 

---