Journal of Engineering Geology

---

Introduction
When loose sand is subjected to seismic shaking, it tends to volume reduction and settlement. The density of the under layers is revealed in the settlement of the ground surface that causes the destruction of the structures located on ground surface. In dry sand layer, settlement in severe shakings occurs under a constant and effective stress condition and very rapid stress. In this regard, the sand deposit settlement is completed before the end of an earthquake, but if the sandy soil layer is saturated and drainage is limited the condition is prepared of fixed volume situation and the major effect of the seismic shocks is generation of exceed pore water pressure. Therefore, the deposit settlement of saturated sand requires a longer time, varying from a few minutes to a few days, depending on the permeability and compressibility of the soil and the length of the drainage path. The main purpose of the present study is to evaluate liquefaction hazards along Tabriz Metro Line 2 with using Standard Penetration test (SPT) and shear wave velocity (Vs) methods. Also, the probable rate of settlement in the soil layers in study area and correlation with liquefaction potential index (LPI) according to both procedure have been determined and discussed in the following paragraphs.
Material and Methods
In order to evaluate the liquefaction potential of soils using two field methods, geotechnical information of 54 boreholes in Tabriz Metro Line 2 were collected. The types of soil and geotechnical properties can affect the liquefaction potential. In this study, the gravely sand, silty sand, silty and sandy soils were studied. Ground water level is one of the main parameters in in soil liquefaction potential evaluation of soils. Variation of water level in boreholes is 2 to 30 meters. The peak ground acceleration (PGA) is necessary for the analysis of boreholes to evaluate liquefaction potential of soils. PGA values were selected in each boreholes position according to the Iranian Code of Practice for Seismic Resistant Design of Buildings (Code-2800-ver.4) equal to 0.35g (for return period 475 years and design life 50 years). Liquefaction potential of soil layers based on SPT results with appliying Idriss and Boulanger (2010) method has been assessed. Andrus, Stokoe and Jung (2004) procedure was used in shear wave velocity (Vs) method (with assuming cementation and un cementation condition in soils). Liquefaction potential index (LPI) of soil layers was calculated for both field tests results. Then, probable rate of settlement due to liquefaction in saturate soil layers was determined. Tokimatsu and Seed (1978) method applied for SPT results, Yoshimine (1992), Yoshimine et al. (2006) and Yi (2009) procedures have been used in Vs test. Finally, correlation between rate of settlement and LPI results were determined.
Results and discussions
Outcomes of this study can be explained in below:
1. Results obtained from comparison of both methods in liquefaction potential evaluation have been showed, agreement between two methods have been happened rarely. Specially, with assuming cementation condition in soils, LPI obtained from Vs method is more than SPT. Although, different factors can be affected at uncertainties in SPT results such as type of drilling machine, energy efficiency and accuracy of test performing. Also, in shear wave velocity method, maximum velocity for occurring liquefaction in soil layers () related to fines content percentage. It is possible that boundary values in procedure not compatible with geotechnical properties in study area.
2. Evaluation of probable rate of settlements in soil layer in study area have been showed that settlement values obtained from Vs is more than SPT. This condition is compatible with LPI amounts.
Conclusions
In sum up, settlement due to liquefaction in saturate soil layers is one of the important phenomena in geotechnical earthquake engineering. Maximum rate of settlement in soil layers in study area is equal 0.45m based on SPT method and 0.9m according to Vs procedure which should be considered. Accordingly, serious damages can be inflicted to buildings, underground structures and life lines in study area.  Therefore, it is suggested in future researches with using empirical and numerical (or soft computing) methods based on field and experimental tests results a detailed assessment conducted and influence of various parameters on settlement of soil layers be determined and the items listed below should be considered:
- Cementation parameter (C) values of soils in shear wave velocity method maybe not compatible with geotechnical properties in study area. It should be evaluated exactly. 
- In this research, peak ground acceleration (PGA) value was selected based on code 2800-ver.5. As regard to Tabriz Metro Line 2 is beside to Tabriz North Fault, PGA value according to historical earthquake catalogue and seismic risk analysis should be evaluated and seismic hazard have to determine with accuracy. 

---

Introduction
The standard penetration test (SPT) is one of the most common tests in geotechnical investigations. The results of this test are known as a simple, inexpensive, and tangible criterion in geological and geotechnical engineering. Many computational methods and engineering judgments depend on the results of this test. In this research, estimation of physical and engineering properties of clay soils was carried out using statistical methods based on standard penetration test results. The scope of this case study is related to a variety of clayey soils in Tabriz (the northwest of Iran). The existing relationships were confirmed based on database of this study. After statistical analysis of the database, eight relationships including single and two-variable associations have been proposed to estimate the physical and engineering properties with better performance using nonlinear regression.
Material and Methods
Different types of clayey silt and marl layers spread in Tabriz were included for the purpose of this study. The geological age of these layers dates back to the Miocene and Pliocene era. This research was conducted in two sections of the field and analysis. Two machine boreholes were drilled, and, based on ASTM, a standard penetration test with other laboratory tests were performed on the soil specimens in order to determine the physical and plasticity properties. According to the results of this study and the existing data, a total of 107 series were prepared. Based on the soil properties, 11 variables were selected including the fine grain percentage (FGP), liquid limit (LL), plastic limit (PL), percentage of clay particles (C), plastic index (PI), consistency index (CI), activity (A), dry unit weight (γ_d), natural moisture content (w_n), initial void ratio (e_o), and effective vertical stress (σ'_v). The standard penetration tests were run for each meter in drilled boreholes. The results of this test were corrected according to NCEER method. The correlation between the variables and corrected standardized penetration test results (N₆₀) were studied by Spearman ranking coefficient. Verifications of the existing eight experimental relationships between  standard penetration and other soil properties, proposed by Kayabasi (2015) and Hoshmand et al., (2012), were checked out using the findings and data of the present study. The linear, exponential, logarithmic, and exponential regressions between each variable and N₆₀ were investigated using SPSS software, version 16. The best regression with the highest R² for each variable was selected. Eight new relationships were proposed. Performance of the suggested relationships was compared with the existing relationships.
Results and Discussion
The findings of the current study could be summarized as:
1. The clay soils of the studied area in Tabriz were classified into four categories including CH, MH, CL, and ML according to USCS classification. The range of changes in plastic index and liquid limits of the samples were 9.19 ~ 45% and 29 ~ 77%, respectively. The corrected standard penetration test results (N₆₀) changed from 9 to 28 showing that soil compression was low to high.
2. The highest positive and negative Spearman correlation coefficients were related to the consistency index (+0.772) and moisture content (-0.759), respectively.
3. The existing empirical relationships, based on the database of this study, were found to have better statistical coefficients in terms of consistency index, activity, moisture percentage, and fine grained percentage. In term of sample depth, the experimental relationship, showed the lowest statistical coefficient.
4. Four single-variable and two-variable relationships were proposed by nonlinear regression analysis. Using these relationships, clay soil properties including activity, moisture content, fine grain percentage, and consistency index were estimated based on N_60. In addition, two relations were proposed between sample depth (D) and vertical effective stress (σ'_v) with N₆₀. The statistical coefficients of the suggested relationships were better than the existing empirical relationships. The proposed relationship of estimating the consistency index with coefficient (R²) of 0.673 and regression line slope of about 1 had the best performance.
Conclusion
In general, the main objective of this study was to investigate the correlation between physical and plasticity properties of clay soils and N₆₀ on Tabriz clayey soils. Clay soils of the present study included various silty and marl layers. Sufficient correlation was observed between the physical and engineering properties of clay soils and N₆₀. The validation of the existing experimental relationships based on A, w_n, FGP, and D resulted in weak statistical coefficients (R² <5) employing the database of the current study. Six new experimental relationships were proposed to estimate A, FGP, w_n, and CI as well as two correlations of N₆₀ with effective stress and sample depth. Generally, the results have been revealed that the statistical coefficients of the proposed relationships were improved compared with the existing relationships. The most suitable relationship was the estimation of soil consistency index (R²~70) and root mean square error (RMSE=129). Finally, due to the novelty of this research topic, verification and development of the proposed relationships for the soils has been recommended in other areas.

---

Introduction
Pavement layers as a part of road structure play an important role and provide a flat and secure surface. Subgrade layer could act as a compacted embankment, natural or stabilized ground. Subgrade is a foundation of pavement layers, and it withstands all of loads due to vehicles that are transferred from upper layers (i.e., subbase, base and asphalt layers).Therefore, constructing pavements with bearing capability, high durability, quality, and maintenance in proper operating conditions is very important. However, suitable materials for constructing pavement layers are not available, and improvement techniques should be employed for them. Generally, different methods such as mechanical or chemical are available for improvement. Nowadays, geosynthetic materials such as geotextile and geogrid are used to optimize and enhance the bearing capacity of pavement layers. The present study is aimed to investigate the effects of geotextile applications on bearing capacity of clay-gravel mixtures in pavement layers.
Material and Methods
In this research, materials were prepared from Barandouz area. Clayey soil was mixed with gravel in 25, 50 and 75 percentages (by weight). Geotextile was woven and made of polypropylene (with commercial name Fibertex-F-32). Geotextile effects in mixture were evaluated in two conditions. Position number one indicates the arrangement of geotextile.  This means, at first, one geotextile layer was embedded in the middle of materials. Then, two and three geotextile layers in equal depths from each other were used in soil mixtures. Position number two shows the mixing pieces randomly. This means that geotextile pieces in 1×1 and 5×5 cm² were prepared and were randomly mixed with materials in 1, 2 and 3 percentages (by weight). For evaluating geotechnical behavior of improved clay-gravel mixtures, compaction and California bearing ratio test (CBR) (in dry and saturate conditions) based on ASTM were performed.        
It should be noted CBR test in dry and saturate conditions were carried out in three different compaction energies (i.e. 10, 25 and 56 blow count for per layer). Moreover, CBR was evaluated for piston penetration at 2.5 and 5 cm in the specimen.
Results and discussion
The findings of this study could be summarized as:
1. Results of compaction test showed that, in the unimproved position, with increasing gravel content in clay, maximum dry unit weight (γd_max) has been increased, while simultaneous optimum water content (w_opt) decreased.
In the improved position, in the first mode, when a geotextile layer was embedded in the middle of the specimens, γd_max reached to its upper value, whereas w_opt reached to its minimum value. On the other hand, with an increase in the number of geotextile layers in clay-gravel mixtures, dry density has been decreased, but optimum water content increased. Furthermore, in the second mode, when geotextile pieces with 1×1 and 5×5 cm² were randomly mixed in the specimens, the findings revealed that geotextile pieces with 1 cm² areas and 1% by weight in clay-gravel mixtures increases γd_max and reduces w_opt.
2. In dry and saturate conditions, California bearing ratio (CBR) test result displayed that in the unimproved condition, with an increase in gravel content in the clay, CBR value has been increased. In the improved situation, in the first mode, when a geotextile layer was embedded in the samples, CBR had a maximum value in all of the compaction energies even though it is reduced as the number of layers increased. In the second mode, when geotextile pieces in 1×1 cm dimensions with 1% (by weight) were randomly mixed with the specimens, CBR value reached at high.  In contrast, with increasing dimensions of pieces and percentages in the presence of geotextile in clay-gravel mixtures, CBR values declined.  Therefore, it can be concluded that, according to Code 234 (Iran Highway Asphalt Paving Code), the application of one geotextile sheet in the middle of materials or geotextile pieces in 1×1 cm dimensions with 1% (by weight) random mixing  is suitable for subbase and base layers in pavement design.
3. CBR test results in the saturate condition in clay-gravel mixtures illustrated that, in the non-reinforced condition, with an increase in clay content in specimens, swelling value keeps rising sharply. On the contrary, in the reinforced position with embedding a geotextile layer in the middle sector of samples or through adding geotextile pieces (1 cm²) with 1 % content  (by weight) to the specimens, the rate of swelling significantly decreased.   
Conclusion
To sum up, the main objective of the present study was to investigate the impact of geotextile applications on bearing capacity of clay-gravel mixtures in pavement layers. The findings demonstrated that when geotextile as a layer was embedded in the middle part of specimens or as pieces with 1×1 cm dimensions and 1% content (by weight) was randomly mixed with the mid materials, the bearing capacity of the reinforced specimens was enhanced.  In contrast, in the saturate condition, swelling potential significantly was reduced. It is noteworthy to mention that 1 cm² pieces of geotextile is more effective than the layers. This is due to the fact these pieces make aggregates closer to each other. Thereby, minimum void ratio (e_min) reaches its least value, the structure of grading improves, and the contacts between particles and geotextile pieces rise. As a suggestion for further research, it looks promising to evaluate the dynamic properties and the behavior of the improved materials with other geosyntheticses.
 

---

In this study, it is attempted to analyze sensitivity and reliability in order to evaluate the liquefaction potential in soil layers in Tabriz. 62 boreholes that had possible conditions for liquefaction were selected. Seismic mapping was simulated using finite fault method and then the effect of soil layers on PGA was estimated. In continue, the liquefaction potential index was estimated and the zoning map of liquefaction risk was presented. In final, through sensitivity and reliability analysis of the Monte Carlo method, the rate of density function against safety factor of the soil layers versus to liquefaction was determined.

---