Journal of Engineering Geology

---

Huge amount of scarp tires have made a big environmental problem that its reduction assessment is an important aim for researchers. Reusing of these materials in different branch of industrial is concentrated. Because of elastic behavior and absorbing the energy, waste tires are concerned as a famous damper. As a more effective manner to check the seismic forces is reducing the entered movements to the structure, it seems that using soil-scrap tire mixtures acts as seismic damper and these mixtures can reduce the earthquake forces acting to the structure. In this research this mixture is modeled as a damper under the foundation and the bed rock depth is analyzed in both time and frequency domains. Obtained results show that using soil-tire mixture reduces the resonance amplitude, energy and the acceleration reached to the ground surface. It should be mentioned that increasing the depth of the bed rock, the act of the mixture as a damper reduces. The results show that this mixture has a better operation in strong ground motions. Using these materials is cheaper compared with other base isolation systems, more available, not requiring high technology and also it reduces both earthquake forces and environmental problems.

---

Comprehensive laboratory tests were performed to assess fatigue behavior of Alvand monzogranite rock subjected to uniaxial cyclic loading. A series of static loading tests was done to obtain the required data for the fatigue tests. Three maximum load levels (85, 90, 95% uniaxial compressive strength (&sigmac)) at amplitudes 70% were used with 1Hz cyclic loading frequency. The results indicated that maximum stress level significantly influenced fatigue behavior of this rock. It was found that fatigue life decreases in a power function with increasing maximum stress level. Accumulative fatigue damage process shows three stages of behavior including crack initiation phase, uniform velocity phase and acceleration phase. Fatigue damage process were analyzed according to axial and lateral maximum and minimum strain, tangent and second modulus, toughness and hysteresis energy in both loading and unloading conditions. Among these parameters, lateral strain, axial strain and second modulus show the best three-stage fatigue damage behavior. Also, it should be noted that most of the cracks generated in parallel to loading direction and lateral strain are affected by more than axial strain.  

---

./files/site1/files/5Extended_Abstract.pdf Extended Abstract
 (Paper pages 91-114)
Introduction
Due to possibility of occurrence in various natural environments and the variety of natural and artificial factors that affect landslides, landslide has special importance in natural hazards. Depending on the landform, several factors can cause or accelerate the landslide. According to previous researches, Human activities, land morphology, geological setting, slope, aspect, climate conditions, proximity to some watershed features such as rivers and faults are the most important parameters. Landslides occur frequently each year and they can cause heavy losses which compensating some of them requires a lot of money and time.
Assessing landslide related hazards with only limited background information and data is a constant challenge for engineers, geologists, planners, landowners, developers, insurance companies, and government entities.
The landslide occurrence in terms of time and place are not easily predictable, for this reason, Landslide Hazard Zonation (LHZ) or Landslide Susceptibility Zonation (LSZ) maps are used to predict the happening of landslides. A landslide susceptibility map depicts areas likely to have landslides in the future by correlating some of the principal factors that contribute to landslides with the past distribution of slope failures. These maps are basic tools for land-use planning, especially in mountain areas. Landslide susceptibility mapping relies on a rather complex knowledge of slope movements and their controlling factors. The reliability of landslide susceptibility maps mostly depends on the amount and quality of available data, the working scale and the selection of the appropriate methodology of analysis and modeling.
Such maps are obtained by dividing of a region into near-homogeneous domains and weighting them according to the degree of possible hazard of a landslide. There are two ways to do landslide hazard zonation: (i) a qualitative approach that is based on expert knowledge of the target area and portrays susceptibility zoning in descriptive terms; and, (ii) a quantitative approach based on statistical algorithms. In the present study of landslide susceptibility zonation, bivariate statistical methods (information value, density area, LNRF, frequency ratio) were used. In bivariate statistical analysis, each factor map is combined with the landslide distribution map and weighting values based on landslide densities are calculated for each parameter class.
Materials and Methods
The best method for studying landslides, which has long been of interest to researchers, is hazard zonation. In this method due to the affecting factors in landslide occurrence, the study area is classified into areas with low to very high risk. Such zonation could be of great help in regional planning. Different methods have been developed for this purpose. In this research four bivariate statistical methods namely information value, density area, LNRF, and frequency ratio are used to investigate the hazard zonation in Poshtdarband region, Kermanshah province. The study began with the preparation of a landslide inventory map. The instability factors used in this study included geology, land use, normalized difference moisture index (NDMI), slope gradient, aspect, distance from faults, distance from surface water, distance from roads, profile curvature and plan curvature. Landslide area ratio was calculated in classes of effective factors maps and weighted by four bivariate statistical methods. In addition, landslide hazard zonation maps were obtained from algebraic sum of weighted maps with regard to breakpoints of frequency curve. Finally, by using density ratio (D_r) Index through all four methods hazard classes were compared and with the help of quality sum (Q_s) and precision (P) indexes these four methods were compared and evaluated.
Results and Discussion
If the landslide susceptibility analyses are performed effectively, they can help engineers, contractors, land use planners, etc. minimize landslide. In this study, bivariate statistical methods were applied to generate landslide susceptibility maps using the instability factors. The bivariate approach computes the frequency of landslides with respect to each input factor separately, and the final susceptibility map is a simple combination of all the factors irrespective of their relative significance in causing landslides in a particular region.
In table 1 subclasses of instability factors which had the highest value in different methods, are summarized.
The density ratio indexes (D_r), quality sum indices (Q_s) and precision indices (P) were used to compare the methods. By overlaying the landslide inventory map of the study area and landslide hazard zonation maps, quality sum (Q_s) and precision (P) indices introduce a suitable model for the studied region, and density ratio index (D_r) introduces division precision among the zones or hazard classes in each zonation model.
Table1. subclasses of instability factors in different methods which had the
highest value

factor methods	aspect	Slope	distance from surface water	land use	plan curvature	profile curvature	distance from fault	distance from the roads	NDMI
information value	N, NE	>40	>1000	forest	concave	concave	<500	>1000	-0.17_ -0.408
density area	N, NE	>40	>1000	forest	concave	concave	<500	>1000	-0.17_ -0.408
LNRF	SW, S	10-20	>1000	pasture	Convex	convex	<500	>1000	-0.17_ -0.408
frequency ratio	N, NE	>40	>1000	forest	concave	concave	<500	>1000	-0.17_ -0.408

The density ratio for information value method in the very high hazard class is accounted 1.700495. These values for density area, frequency ratio, and LNRF methods are, 3.407827, 3.402257, and 1.694628 respectively.
Method precision (P) values for information value, density area, frequency ratio, and LNRF methods are 0.160826, 0.241024, 0.240672 and 0.16942 respectively.
Conclusion

• Frequency ratio, density area and information value methods showed that forest land use, slope and slope shape factors have the highest impacts on a landslide occurrence.
• The LNRF method showed that geology factors, pasture land use and distance from surface water had the greatest role in landslide making.
• For frequency ratio, information value, and density area methods, the effective factors in landslide are the same, however through the LNRF method, the three factors which have the greatest impact on landslide happening, are generally different from the three other methods.
• The density ratio values show that density area and frequency ratio methods respectively have more accuracy and applicability within all used methods for separating hazard classes in the study area.
• The quality sum (Q_s) results indicate that although there are minor differences, the frequency ratio compared to the density area method was more accurate and more applicable for separating landslide hazard in the Poshtdarband region.
• The calculated results of P index indicated that among the used methods, the density area method with a nuance of the frequency ratio method is the most suitable method for the study area.

---

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
One of the main problems in the Golestan province watersheds is the high degree of erosion and soil degradation, so that the equilibrium between the soil process and the soil erosion is unbalanced, and the erosion rate increases from west to east. Among these, the gully erosion and piping have the highest role. Gully is a canal or stream with the headcut with active erosion, sharpened slope and steep walls that results from the destruction of surface flow (usually during or after the occurrence of precipitation), dissolution movements, and small mass movements. The extent of gully in the eastern parts of Golestan province has caused the land degradation of arable land and landscape and has increased the conservation cost and etc. Because of connecting upstream areas of the basin to the downstream areas, gully has particular importance, which provides the possibility of sediment and pollutant transport, road destruction and financial losses to agricultural lands. In order to prevent and control the development of gully processes from a small scale to large one, it is a versatile utility to identify and extract the areas prone to gully erosion.
Due to the high intensity of gully erosion and its increasing growth in the Gharnaveh watershed, the Garnaveh River has an unstable status and severe eroded gully, and in some areas it has a great depth and vertical lateral walls, as well. Therefore, in this research, the watershed of Garnaveh was selected to prepare the risk areas of gully erosion.
The aim of this research is to determine Gully Erosion Hazard zoning using Frequency Ratio and Gupta & Joshi methods (Gully Nominal Risk Factor-GNRF) in the Garnaveh watershed (Golestan province). Ultimately, the accuracy of the model has been evaluated using quality sum method and Kappa coefficient.
Material and methods
The study area is located in the northern part of Iran, Golestan province. The Garnaveh watershed with an area of about 78430 hectares lies between longitudes 370360 E and 414472 E, and latitudes of 4183819 N and 4155267 N (UTM Zone 40).
At first, gully erosion inventory map with the scale of 1:75,000 (dependent variable) for the Gharnaveh watershed has been prepared using multiple field surveys and satellite images. From total gullies, 70% have been selected randomly for building gully erosion hazard zoning model and the remaining ones (30%) have been used to validate the provided model.
In this research, seven data layers including slope percent, slope aspect, plan curvature, lithology formation, land use types, distance from rivers and distance from roads have been selected as gully erosion controlling factors (covariates/ independent variables) and then they have been digitized in ArcGIS software. The amount of Gully density of each factor class has been calculated from a combination of independent and dependent variables, and the rating of classes have done based on Frequency Ratio and Gully Nominal Risk Factor equations. Finally, the Gully erosion hazard zoning map has been drawn from the summation of weighting maps in ArcGIS. In this map, the value of each pixel is calculated by summing the weights of all the factors in that pixel. The pixel values are categorized based on the natural breaks classifier into very low, low, medium, high and very high hazard zones. Then, an accuracy of zoning map has been evaluated by quality sum method and Kappa coefficient.
Results and discussion
The result of affecting factors classification of the gullies shows that loess deposits formation, rangeland, areas with low distance from road and rivers, northwest aspect, low slope amplitude and concave slopes contain the most susceptibility to gullying. The results of frequency percent comparison of gullies in hazard classes show that from all gully zones in the validation step of the GNRF and frequency ratio models %74.52 and %78.11 of zones are located in the high and very high risk classes, respectively. The result of model validation using the quality sum method and a Kappa coefficient show that the frequency ratio model is a more appropriate model for gully erosion hazard zoning (with the quality sum and a Kappa coefficient of 3 and 0.89, respectively) than the GNRF model (having the quality sum and Kappa coefficient of 1.27 and 0.74, respectively).
Conclusion
In this research, the areas susceptible to gully erosion in the Gharnaveh watershed have been mapped with the frequency ratio and GNRF (for the first time) models. For this purpose, 7 affecting factors (independent variable) and 805 gully zones (dependent variable) were provided to measure the hazard maps of gully erosion. The following results are obtained from this study.
- The geology factors were identified as the most effective factors in the occurrence of gully erosion in the Gharnaveh watershed.
- Based on the gully erosion zoning hazard map of the Gharnaveh watershed, more than 70 percent of gullies are situated in the very high and high hazard classes.
- The produced gully erosion hazard map is useful for planners and engineers to reorganize the areas susceptible to gully erosion hazard, and offers appropriate methods for hazard reduction and management, as well../files/site1/files/131/4Extended_Abstract.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
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
One-dimensional site response analysis is widely performed to account for local site effects during an earthquake. Most of these approaches assume that dynamic soil properties are frequency independent. Laboratory test results as well as in-situ testing show that shear modulus and damping ratio are dependent on the frequency of loading. Although the amplification factor at ground surface with respect to frequency dependent dynamic properties of mixed alluvium materials under different near-fault motions with various velocity period is recognized, it is not well characterized and quantified.
Material and analysis methods
In this study, the tests results of samples which obtained from the drilling borehole (BH14) form Pardis city in Iran, are used. The soil is classified as clayey of high plasticity/clayey sand (CH/SC) and almost uniform and similar in the whole log profile.
Shear modulus and damping ratios versus shear strain curves (ASTM D3999) of CH/SC natural materials at effective confining pressures of 1, 2 and 5 kg/cm² with frequency of 0.5, 2, 5, and 10 Hz were used in one dimension response analyses using EERA Code.
Generally the damping ratio versus shear strain of the studied materials under low loading frequency (i.e. 0.5 Hz) almost falls in the range identified in literature. However, at higher loading frequencies (5 and 10 Hz) the damping ratios completely fall above the known upper bound trend. It is observed that, in general, the G and D values increase as loading frequency increases. Moreover, at certain strain G/G_max ratio decreases as loading frequency is increased.
Different dynamics behaviour curves were used in analyses, in isotropic consolidation conditions. In order to assess the amplification, acceleration spectra, acceleration spectra ratio, coefficient of B, at ground surface under eight well-known near-fault ground motions, 1728 one dimensional analyses were carried out with EERA code. The analyses have been performed for three base acceleration levels, namely, 0.1 g, 0.35 g and 1 g, using the simple time history scaling method. Field and laboratory test results of shear wave velocity were used in the analyses.
In this study, several well-known near-fault motion records are utilized for ground response analyses. Near-fault earthquakes records were selected from the strong motion database of the Pacific Earthquake Engineering Research Center (PEER) and Iran Strong Motion Network (ISMN) for specific reasons of location of the near faults sites.
In current building codes, the upper 30 m soil deposits overlying the higher impedance earth crust are regarded as the most relevant and significant in characterizing the seismic behavior of a site. Therefore, it is useful to accomplish investigations for obtaining their amplification and spectral acceleration for 30 m and even thicker (e.g. 60 m, for usual deep excavation in Iran), in order to have economical and safe designs and constructions.
Results and discussion
Figure 1 presents a comparison of normalized spectral acceleration (B factor) versus period for 30 m and 60 m thick profiles and Vs testing for frequencies dependent and independent analyses under input base acceleration of 0.35g for longitudinal component of used earthquakes. B factor of Iranian Standard 2800 and UBC97 also has been presented in the figure. The spectral acceleration at short period for frequency dependent analysis is higher than that of the frequency independent analysis. The  increases in frequency dependent analysis and higher thick profile (i.e. 60 m).
Conclusion
Results show that the effect of loading frequency has a considerable influence on the acceleration response at the ground surface. For both 30 m and 60 m soil columns, the increase of the loading frequency, decreases the amplification factor especially in the short period zone of the spectra. Based on the acceleration response spectra of near field strong motions derived for soils types of I and IV in this study, the period corresponding to  in the design spectrum of Iranian Standard 2800 should increase to 0.5 and 1.4, respectively. Therefore, selection of the appropriate G and D curves measured at frequency similar to those of the anticipated cyclic loading (e.g. seismic) has a paramount importance../files/site1/files/134/1.pdf

---

One of the effective parameters in the dynamic behavior of reinforced soil walls is the fundamental vibration frequency. In this paper, analytical expressions for the first three natural frequencies of a geosynthetic reinforced soil wall are obtained in the 3D domain, using plate vibration theory and the energy method. The interaction between reinforced soil and the wall is also considered by modeling the soil and the reinforcement as axial springs. The in-depth transverse vibration mode-shapes, which were impossible to analyze via 2D modeling, are also analyzed by employing plate vibration theory. Different behaviors of soil and reinforcements in tension and compression are also considered for the first time in a 3D analytical investigation to achieve a more realistic result. The effect of different parameters on the natural frequencies of geosynthetic reinforced soil walls are investigated, including the soil to reinforcement stiffness ratio, reinforcement to wall stiffness ratio, reinforcement length, backfill width and length to height ratio of the wall, using the proposed analytical expressions. Finally, the results obtained from the analytical expressions proposed are compared with results from the finite element software Abaqus and other researchers’ results, showing that the proposed method has high accuracy. The proposed method will be a beginning of the 3D analytical modeling of reinforced soil walls.
 

---

Introduction
Soil dynamic properties are used to evaluate the dynamic response of soils at different strain levels in geotechnical engineering. The shear modulus (G) and damping ratio (D) are among the most important dynamic properties of soils. In general, the factors affecting the dynamic behavior of soils are divided into two categories: first; soil type and characteristics such as water content, void ratio and soil plasticity and second; parameters of loads applied on the soil such as the number of loading cycles, loading frequency and loading waveform .Therefore, it is widely believed that the dynamic response of soils partially depends on the characteristics of the load. In this paper, 1-g shaking table tests were employed to investigate the effect of loading waveform and frequency content on dynamic properties of dry sands. The response obtained from soil samples during loading with different frequencies, input accelerations and waveforms were used to generate hysteresis loops of tested samples at different strain amplitudes. Then, hysteresis loops were used to determine the damping ratio and shear modulus at different strain levels. Finally, the effects of loading frequency and waveform on the changes of each parameter (G and D) were investigated.
Material and methods
A hydraulic shaking table with a single degree of freedom, designed and constructed at the Crisis Management Center of Urmia University, was used to conduct the experiments. Firoozkuh No. 161 sand was used in all the experiments. The Firoozkuh sand gradation curve is similar to that of Toyoura sand. In this study, accelerometers were used to measure the acceleration of the input to the sample as well as to record the acceleration caused by the input excitation at different depths of the soil sample. The displacement transducers (LVDT sensors) were also used to measure linear displacement. Each soil sample was constructed using dry Firoozkuh sand poured uniformly into the container from four equal heights of 150 mm to reach a total height of 600 mm. During the compaction process, the accelerometers A₁, A₂, and A₃ were placed at a depth of 150, 300 and 450 mm with respect to the bottom of container. Also, one accelerometer, A₀, was attached rigidly to the container base to measure base acceleration. A displacement transducer (L₁) was placed on the soil surface at a height of 600 mm from the floor of the container to measure the vertical displacement of the surface of the soil. In this study, 42 shaking table tests were performed to study the effect of loading frequency and waveform on dynamic properties of dry sand. The test samples were subjected to rectangular, sinusoidal and triangular loading at frequencies of 0.5 to 9 Hz and at input acceleration of 0.1g and 0.3 g.
Results and discussion
Given the importance of G-γ and D-γ curves in dynamic analyses, the changes in shear modulus with shear strain has been studied. The results show that the shear modulus increases as the frequency increases in all cases, and this increase is observed at lower frequencies and increases with increasing frequency. On the other hand, the shear modulus decreases with increasing shear strain. At a constant testing frequency, soil samples subjected to the rectangular waveform exhibited the largest shear modulus while the samples subjected to the triangular waveform showed the least shear modulus. The shear modulus of the samples under the sinusoidal waveform is barely more than the shear modulus of samples under triangular waveform. Moreover, by increasing the shear strain, the shear modulus values ​​of samples with different waveforms have become closer and the waveform effect is reduced. As for the effect of input acceleration on the shear modulus , increasing the input acceleration increases the shear strain and consequently, decreases the shear modulus in all states (the values ​​of shear modulus in various frequencies and the waveforms under the input acceleration of 0.1 g are larger than the shear modulus values ​​under the input acceleration of 0.3g). In the case of the damping ratio, the results show that, in all cases, damping ratio increases with shear strain. At low strain levels, the damping ratio values at various frequencies and waveforms are low and yet very close. At higher strain levels, the increase in frequency increases the damping ratio. This increase is more significant at higher frequencies. Also, the effect of waveform on the damping ratio is more apparent at larger shear strains, and at such shear strain levels, soil samples under rectangular loading exhibit the largest damping ratio. The damping ratio associated with the sinusoidal and triangular loading are also close to each other and it is a slightly larger for sinusoidal loading. On the other hand, the damping ratio increases with input acceleration. In addition, the effect of increased input acceleration on the increase in the damping ratio is more obvious at higher frequencies mainly due to the increase in shear strain.
Conclusion
In the present study, the effects of loading frequency and waveform on the dynamic properties of dry sand were investigated using shaking table tests. The following conclusions were drawn:

• The shear modulus increases with frequency. The trend is more obvious at larger frequencies. The effect of loading frequency on the damping ratio of the soil at low levels of strain is negligible, and at relatively large strain levels, damping ratio increases with loading frequency.
• Soil samples exhibit the highest shear modulus and damping ratio under rectangular loading. Therefore, in all the tested frequencies and input accelerations, the values of G and D for the rectangular waveforms are greater than those of the sinusoidal and triangular waveforms. The shear modulus and damping ratio for the sinusoidal waveforms are marginally greater than those of triangular waveforms, yet one can consider them practically similar.
• In all cases, the shear strain increased by increasing the amplitude of the input acceleration, and as a result, the shear modulus decreased and the damping ratio increased../files/site1/files/142/2.pdf