Journal of Engineering Geology

---

Extended Abstract
Introduction
The presence of potentially toxic elements in the environment and especially in soil has been one of the greatest concerns due to their health implications. Potentially toxic elements from anthropogenic sources tend to be more mobile than those from lithogenic or pedogenic sources.  Generally, the distribution of potentially toxic elements is influenced by the nature of parent materials, climatic conditions, and their relative mobility depending on soil parameters, such as mineralogy, texture and class of soil. In the inhabited, and industrial areas, vicinity to the un-engineered landfills, excess accumulation of toxic elements in surface soils can directly threaten wellbeing of exposed inhabitants via ingestion, inhalation and dermal contact routes. A few studies conducted on risk assessment of potentially toxic elements in soils of Kermanshah province, west of Iran. Soil in the study area is susceptible to contamination by anthropogenic activities in the form of industrial wastewater, agricultural activities, solid waste, runoff, atmospheric deposition and especially un-engineered landfills. The presence of toxic elements in soil around of un-engineered landfills without proper consideration to the environmental protection measures, will certainly lead to a significant environmental hazard in Kermanshah province. Therefore, the main purposes of this study are to evaluate the contamination levels, health risk assessment, and source identification of As, Cd, Cr, Cu, Ni, Pb and Zn in the Gasre Shirin, Gilane Gharb, Paveh, Javanrood, Eslamshahr, Ravansar, Kermanshah and Sanghar un-engineered landfills.
 
Material and methods
     A total of 30 topsoil samples were collected (0-20 cm depth) from the eight un-engineered landfills of the Kermanshah province. In order to achieve a representative sample, composite samples were prepared by mixing the four subsamples taken at each corners of 2×2 m square cell because composite sampling yields homogenized samples for analyses. The subsamples were mixed and a final sample of 1 kg was taken by repeated coning and quartering. To determine background concentration of heavy metals, eight soil samples were collected from areas far from known sources of contamination (40-60 cm depth).
The collected samples were immediately stored in polyethylene bags and air-dried in the laboratory at room temperature. Then, samples passed through a 2mm stainless steel sieve. The <2mm fraction was ground in an agate mortar and pestle and passed through a 63 micron sieve. In order to determine the concentration of As, Cd, Cr, Cu, Ni, Pb and Zn complete dissolution of soil samples (approximately 1 g of each) was carried out using a mixture of HF, HNO₃, HClO₄ and H₂O₂ in a Teflon beaker on sand bath at atmospheric pressure. The concentrations of the selected elements were measured by an accredited commercial laboratory (Zar Azma Laboratory, Iran) using ICP-MS methods. Data quality was ensured through the use of internal duplicates, blanks, and HRM. The precision and accuracy of measurements are 95% and +/-5% respectively.
The assessment of soil contamination was carried out using geochemical indices including contamination factor (CF), modified degree of contamination (mC_d) and enrichment factor (EF). The methodology used for the health risk assessment was based on the guidelines and Exposure Factors Handbook of US Environmental Protection Agency. The average daily doses (ADDs) of heavy metals received through ingestion, inhalation, and dermal contact for both adults and children were calculated. In this study, hazard quotient (HQ), hazard index (HI) and carcinogenic risk (RI) methods were used to estimate non-carcinogenic and carcinogenic effects of heavy metals. The HQ was calculated by subdividing ADD of a heavy metal to its reference dose (RfD) for the same exposure pathway(s). If the ADD exceeds the RfD, HQ>1, it is likely that there will be adverse health effects, whereas if the ADD is less than the RfD, HQ<1, it is considered that there will be no adverse health effects. A hazard index (HI), the sum of HQs, which means the total risk of non- carcinogenic element via three exposure pathways for single element of <1 indicates no adverse health effects, while HI values >1 show possible adverse health effects. Carcinogenic risk is regarded as the probability of an individual developing any type of cancer in the whole life time due to exposure to carcinogenic hazards and was calculated for As and Cd as follows:
                                                                   ⁽¹⁾
The value of SF represents the probability of developing cancer per unit exposure level of mg/kg day. The acceptable risk range for carcinogens is set to 10^-6 to by the USEPA, so that RI values below 10^-6 do not require further action, while risks greater than 10^-4 are considered to be of concern and require additional action to reduce the exposure and resulting risk.
Results and discussion
The soil pH ranges from 7.01 to 8.06, with an average value of 7.51 suggesting neutral conditions. Organic carbon (OC) contents of soil samples ranged from 0.06% to 4.91% (average 1.59%). In this study, based on the USDA textural triangle the main soil textures are loamy, clay loam and sandy loam, respectively.
The average abundance order of selected elements content is: Zn>Ni>Pb>Cr>Cu>As>Cd. Comparison of mean concentration of the potentially toxic elements in the soil samples with mean worldwide values reveals higher Zn, Pb and Ni contents in this area.
The results of contamination factor indicate very high contamination for Cd, Cu, Pb and Zn. Modified Degree of Contamination (mCd) calculated based on background values proves very high degrees of contamination for selected trace elements in Gasre Shirin and Eslamshahr landfills soil samples The results of enrichment factor evaluation similarity to contamination factor indicate that Cd, Cr, Pb, Cu and Zn have more influence from anthropogenic sources. The maximum EF of Pb, Zn and Cd and Cu is 346.7,124 and 51.9 respectively, which means very high enrichment in Ghasre Shirin landfill soil samples.
Exposure doses of 7 heavy metals in soil samples of un-enggenerd landfills for children and adults were calculated. The total exposure HQs calculated based on adults from ingestion, dermal contact, and inhalation for Cd, Cu, Ni, Zn, As and Pb was less than 1(except Ghasreshirin landfill). The hazard quotient values based on the adult risk for Cr were greater than 1.0. The results show that HQ for Pb and As in children by dermal and ingestion pathway is exceeded 1.0 in soil samples of Paveh, Javanrood, Ravansar, Kermanshah and Sangher landfills and Ghasreshirin and Eslamshahr landfills, respectively.
Conclusion
The concentration, pollution level, potential sources and health risk of potentially toxic elements in eight landfills top soil of Kermanshah province were investigated in this study. The following conclusions were drawn from this research.
- Compared with the background values of As, Cd, Cr, Cu, Ni, Pb and Zn in soils of Kermanshah Province, landfills soil have elevated metal concentrations as a whole.
- According to high contamination level and health risk of some studied potentially toxic elements, and also due to the proximity of contamination sources to residential district of the study area, more attention should be paid to manage and reduce contamination.
- These results provide basic information of toxic elements pollution control and environment management in the area../files/site1/files/121/Shakerii_Abstract.pdf

---

Introduction
Estimation of Liquefaction is one of the main objectives in geotechnical engineering. For this purpose, several numerical and experimental methods have been proposed. An important stage to predict the liquefaction is the prediction of excess pore water pressure at a given point. In general, there are two important methods for soil dynamics analyses, fully coupled effective stress and uncoupled total stress analysis. The main purpose of this study is to evaluate the model capacity of the finite difference software, FLAC, based on effective stress analysis methods to predict the excess pore water pressure during seismic loading. A level ground centrifuge test conducted during the VELACS project on the Nevada sand with a density of 40%, was utilized to calibrate the numerical model. After the validation of the numerical model, a model was conducted to predict excess pore pressure and consequently the liquefaction for the site of Bandar Abbas Mosque.
Theoretical bases
A fully coupled u–P formulation, where pore pressures and displacements are computed simultaneously and interactively at each time step, is used in FLAC software. This feature is used to simulate the excess pore water pressure time histories during cyclic loading.
The finite difference based software, FLAC, used the Finn model that incorporates two equations correlating the volumetric strain induced by the cyclic shear strain and excess pore water pressure produced during cyclic loading. As mentioned above, the pore water pressure generation can be computed from two sets of equations: Martin et al. (1975) and the Byrne (1991) formulations in which the volumetric strain that was produced in any cycle of loading is depended on the shear strain that was formed during that cycle as well as the previously accumulated volumetric strain.
Modeling and Results
The VELACS model # 1 centrifuge test representing a level ground site constituted of the Nevada sand at 40% relative density has been numerically simulated in the current study to validate the numerical model. The centrifuge model contains a laminar box with slipping “rings” that allows differential horizontal displacements. This was simulated in the FLAC model by free-field boundary conditions which prevent reflection of the waves in the side walls. Figure 1 shows comparison of EPWP time histories ratio of numerical modeling and centrifuge test. Static analysis was carried out before dynamic analysis in order to find initial stress and strain state. At the next stage, the dynamic loads were applied at the base of the model and dynamic analysis was performed. The Bandar Abbas mosque project is located approximately 500 meters from the coast. In the project, due to the groundwater level and the existence of loose layers of silt, investigating the potential of liquefaction is necessary.
For numerical modeling the results of the general soil mechanics test on soil samples and standard penetration test performed on the site were used to calibrate the parameters and select the model constants.
Conclusion
The results of numerical modeling have been matched to experimental results of the centrifuge test using both Martin and Byrne formulations, except for the case of 5 m the numerical model has predicted lower excess pore water pressure values than the experimental values. This may be originated from the fundamental assumption of the Martin et al. (1975) EPWP theory, in which excess pore water pressure is directly related to the relevant volume changes. On the other hand, the Martin et al. (1975) model was adopted for one-dimensional measures of shear strain, while, in a 2D analysis under both horizontal and vertical shakings, there are three strain rate measures. FLAC uses some assumptions to solve this problem and it can affect the results.
The results of the numerical model showed liquefaction to a depth of about 5 meters that is almost compatible with the results from the lab, which has declared that the depth 2 to 5 m is liquefiable.
With careful selection of numerical model parameters one can generally use the simulation results to have a general sense on the pore water pressure generation and liquefaction prediction.
./files/site1/files/134/5.pdf

---

Introduction
The existence of problematic soils due to their geotechnical properties, such as low strength and stability, high compressibility, and swelling, is one of the most important issues and challenges that geotechnical and civil engineers are faced in urban environments, especially in metropolises. Various methods are used to stabilize and to improve the behavior of problematical soils such as compaction, consolidation, stone columns, jet grouting, biological procedures, and additive materials including nanomaterials. Because of their high specific surface, the use of nanoparticles is very effective to increase the shear and mechanical strength parameters of soil. Mashhad city is located on alluvial deposits of Mashhad Plain. A wide area of this city, especially the central and eastern areas where the Imam Reza holy shrine is located, has been built on weak and fine-grained deposits. Considering constructing high-rise buildings such as hotels and commercial complexes in these areas, as well as the need for restructuring the urban decay, the soil improvement will be inevitable. Given the significant application of these nanoparticles, the purpose of this study is to investigate the effects of nanoclay and nanosilica on each other and to find their optimal composition as a suitable alternative for traditional materials to improve the weak and problematic soils. This not only increases the bearing capacity and strength properties but also reduces the cost and time of project implementation.
Method and Materials
To achieve a hybrid with maximum strength and bearing capacity in executable projects, laboratory tests were performed on the soil picked up from the vicinity around Razavi holy shrine in Mashhad mixed with nanoclay and nanosilica. The type of soil is classified as CL-ML based on sieve and hydrometer tests. The nanoclay used in this research is the type of montmorillonite- K10, and the nanosilica is as a powdered shape with 99% purity.
At first, nanoclay and nanosilica were mixed independently with soil in six different weight ratios (0%, 0.1%, 0.5%, 1%, 2.5%, & 5%) and (0%, 0.1%, 0.25%, 0.5%, 0.75%, & 1%). Soil mechanical and strength properties, including compressive and shear strength, settlement, plasticity index, and swelling, were studied by standard laboratory tests on all specimens. After determining the optimum ratio of each nanoparticle, four hybrids consisting of nanosilica and nanoclay were made in four different combinations and then the effects of these four hybrids were investigated on the soil in the laboratory scale (Table 1).
Table 1. Composition of hybrids made with different percentages of nanomaterials

Nanomaterials composition	Hybrid Name
5% Nanoclay + 0.25% Nanosilica	5NC + 0.25NS
5% Nanoclay  1% Nanosilica	5NC + 1NS
2.5% Nanoclay + 0.25% Nanosilica	2.5NC + 0.25NS
2.5% Nanoclay + 1% Nanosilica	2.5NC + 1NS

Conclusion
The results of the Atterberg limit test on improved and pure soil indicate that the addition of nanoclay and nanosilica and the optimized ratios of these nanoparticles hybrid to increase the soil resistance parameters did not change the soil swelling index.
Evaluation of shear strength test results showed a significant synergistic effect of these nanoparticles on increasing the shear strength parameters. The nanoparticles hybrid of 2.5% nanosilica and 1% nanosilica increased the cohesion up to 106% and also hybrids of 5% nanosilica and 1% nanosilica increased the internal friction angle of soil up to 32%.
Examination of unconfined compressive strength tests presented a 134% increase in the compressive strength of the specimen improved with 2.5% nanoclay and a 620% increase in soil improved with 1% nanosilica. The optimum hybrid compositions of 5% nanoclay and 1% nanosilica increased significantly the compressive strength of the studied soil up to 785% and reduced the settlement of the soil by 60% compared to pure soil.

1. Laboratory studies of electron microscopy examination on ​​pure and improved soil samples with nanoparticle hybrid revealed the presence of these particles in pores of the improved soil. On the other hand, the high specific surface area of ​​the nanoparticles increased the interaction of the soil particles, and the effect of adding these nanoparticles on the refining process is observed in compressive strength increase.

As the nanoclay, nanosilica, and hybrid of nanoparticles are the results of soil processing, these particles are very effective to solve the environmental problems because of good compatibility with soil environments. In addition, low volumes of nanoclay, nanosilica, and hybrid in these nanoparticles are necessary to increase the compressive strength and decrease the settlement of soil. Therefore, using these nanoparticles at the project site reduces significantly the cost and execution time of the project.