Search published articles


General users only can access the published articles
Showing 40 results for Subject:

Behrooz Samadian, Ali Fakher,
Volume 13, Issue 1 (Vol. 13, No. 1 2019)
Abstract

Introduction
Geotechnical investigations merely through boring and engineering experiments are considered a difficult task as they are highly costly and time-consuming. The identification of large areas initially requires geological studies followed by the inclusion of geotechnical information. Finally, a geological and geotechnical classification is prepared for the entire area. This type of classifications is employed in strategic urban planning and quick selection of geotechnical variables in small-scale projects. The present research performed the steps involved in these investigations and classifications for the city of Sanandaj, Iran. Hence, the geological-geotechnical classification of the city of Sanandaj was presented by integrating the geological information of this city with the geotechnical data obtained from drilled boreholes as well as multiple wells at different locations in this city.
Materials and Methods
This study was conducted on the city of Sanandaj in six steps. The steps involved and their respective objectives are given in summary in Table 1.
Discussion
This study is applicable to those regions with insufficient information on their boreholes. The present study used only 211 boreholes, the distance
Table1. Steps involved in this study
Objective or result Title Step
Identifying the general geological characteristics General geological investigation of the considered region 1
Determining the rock units and soil layers as well as their outcrops and investigating their appearance Determining the appearance of the layers through field investigations 2
Determining the layer types and drawing the longitudinal and lateral profiles Identifying subsurface layers 3
Determining the characteristics of geological units and their origin of emergence Geological classification based on the steps involved in formation of units 4
a)Collecting the available information, b) controlling the available information, c) completing the information Determining the geotechnical attributes of geological units 5
a) Presenting geological-geotechnical classification, b) presenting geological identification criteria to determine the type of a given unit at the site of the project Presenting a geological-geotechnical classification for the considered region 6
between which was greater than 5 km in some areas of the Sanandaj city. Hence, although no sufficient information was available on many areas of Sanandaj, the proposed method in this study was able to identify the geotechnical attributes of all soil layers and rock units. This study emphasizes on geological and geotechnical classification and presents a step-by-step method to systematically relate geological and geotechnical studies. By integrating these classifications, geotechnical identification of extensive regions such as urban areas can be facilitated even if the number of boreholes is insufficient. Moreover, simple identification criteria can be extracted from this method, through which the engineering attributes of the layers at each point can be determined. This method can be used as an optimal and economical method for geotechnical identification of extensive areas.
Conclusion
The following summaries can be concluded from this study:
-The step-by-step procedure of integrating geological and geotechnical information was described, through which the geological-geotechnical classification for this city was obtained.
-The geological units identified for Sanandaj were shale, limestone, andesite, and Quaternary, which includes layers of alluvial clay, residual clay, and sand and gravel. The extent and distribution of each of the aforementioned units in Sanandaj were identified and plotted. Moreover, the physical and mechanical characteristics of each of the units as well as their geotechnical hazards were determined and presented.
-In this study, simple geotechnical criteria such as faults, altitude level, and distance from river were identified. These parameters were effective in identification of geological units in Sanandaj../files/site1/files/131/5Extended_Abstract.pdf
 
Reza Nassirzadeh Goorchi, Mehdi Amini, Hossein Memarian,
Volume 13, Issue 2 (Vol. 13, No. 2 2019)
Abstract

Introduction
One of the most sensitive and important issues in some civil engineering projects is slope design and application. The process of slope design always involve many uncertainties. Hence, it is impossible to accurately comment on its stability or instability. Most of the uncertainties in the slope stability analysis are related to the nature of materials, geometry, environmental conditions, model errors, and measuring errors as well. Therefore, the slope stability analysis with a deterministic approach which uses the concept of safety factor would often not result satisfactory. Consequently, the use of probabilistic methods is more advised. Accordingly, in recent years, the probability analysis has been used to slope stability analysis. In these analyses, the effective quantities of slope stability are considered as statistical distributions, and the reliability coefficient would then be a statistical distribution. Likewise, one of the approaches to simulate uncertainties in the probabilistic analysis is to use the variation coefficient. If the variation coefficient changes, the probability of failure will change accordingly. When the variation coefficient becomes a larger number, costly solutions are required to reduce the probability of failure. If the variation coefficient becomes low, the reliability will be increased and the required costs to reduce the probability of failure will be decreased. Therefore, determining the amount of variation coefficient in these analyses is very important. Furthermore, the correlation coefficient between the quantities is another effective parameter in computing the probability of failure.
Material and methods
In this research, the stability analysis of the slope facing the spillway of the Shiraz Kavar dam has been done in two probabilistic and deterministic methods. Since circular slip probability is more likely than other types of failure, in the analysis of the stability of this slope, the problem of circular failure is very important, and an appropriate equilibrium program should be used for circular failure analysis. Therefore, SLIDE software was used to slope stability analysis. For material behavior, the Hook-Brown failure criterion was applied. In order to determine the strength parameters of the criterion, Geological Strength Index (GSI), uniaxial compressive strength (UCS) and rock constant parameter mi were used. For crushed rock with a moderate quality of crushing, the GSI quality of the rock mass was about 23 to 38, which the average value of that for the rock mass of the overflow was assumed 35. Also, the uniaxial compressive strength of the rock was evaluated about 50 to 100 MPa with an average value of 75 MPa. In addition, the value of mi was 10, and due to mechanized drilling, the disturbance factor was considered to be 1. The amount of unit weight was assumed to be 22 kN/m3. The initial model used for deterministic and probabilistic analyses, is the Morgenstern-Price model. To conduct probabilistic analyses, Monte Carlo simulation was performed using random sampling method (RS-MC) and 200,000 sampling were used to converge the simulation results. To determine the coefficient of variation and the probability distribution of UCS, GSI and mi, the proposed values ​​of Hook (1998) were applied and for unit weight (γ) James Rodriguez and Sitar (2007) studies were used. Also, the minimum and maximum values ​​of UCS and GSI are determined based on the results of experiments, and Third Sigma rule was utilized for mi and γ quantities. Since the earthquake phenomenon is rarely of great intensity and the number of small earthquakes is higher, therefore the truncated exponential distribution function can be in good agreement with the results of the earthquake. Usually, the maximum magnitude of the earthquake acceleration coefficient is twice that of the average.
Results and discussion
In the presented paper evaluation denotes that the safety factor computed by probabilistic analysis is given as a distribution function. The function provides a clearer view of failure condition. However, a deterministic analysis only illustrates a certain value for the failure. In addition, the results of the probabilistic analysis show that it is possible to optimize the dip of the slope; such that it remains completely stable and the volume of earthwork is also minimized. Therefore, by using probabilistic analysis, the optimal dip of the slope was determined. In these circumstances, the amount of earthwork was decreased by 28,000 cubic meters. Also, the sensitivity analysis of the variation coefficient and correlation coefficient between parameters are analyzed. The results of the sensitivity analysis of the failure probability versus the variation coefficient of the quantities showed that the quantities of sensitivity factor for static conditions is greater than the corresponding pseudo-static, and the GSI amount is the highest, while the specific gravity has the least effect on the probability value. In addition, the analysis indicated that if the GSI coefficient of over 21% is selected, the probability of a static failure is higher than the permissible limit. Also, increasing the variation coefficient of quantities by as much as 50% exhibits that the probability of static failure is still below the permissible limit. Also, the correlation coefficient between UCS and GSI shows that the higher variation coefficient of the quantities is chosen, the more variations of failure probability compared to . In the case of pseudo-static conditions, variations in the failure probability are linear in relation to , while in static conditions, these changes are exponential for an increase of 50% in the variation coefficient. Also, to reduce the coefficient of variation by 50%, the probability of static failure for different values of  is approximately zero.
./files/site1/files/132/8Extended_Abstracts.pdf
Ali M Rajabi, Alireza Sajdeh,
Volume 13, Issue 4 (Vol. 13, No. 4 2019)
Abstract

Introduction
Concrete faced rockfill dams have been considered in recent years more than other types of dams due to their low dependency on the bed and the shape of the valley, as well as the simpler construction technology. In this regard, rockfill dams are a suitable substitute for embankment dams because of higher stability of the body and the availability of rock aggregates. On the other hand, because the permeability of rock aggregates is much higher than other materials, different methods are used to seal these types of dams. One of these methods is the use of non-impermeable concrete facing in the upstream of these dams. This particular type of gravel dams is called Concrete-Faced Rockfill Dams (CRFD). In this study, a contact element with a definition of elastic-plastic failure in the modeling process is proposed to simulate the surface separation and re-contact of the concrete face with the rockfill surface of the dam.
Method
In this paper, behavior of a concrete faced rockfill dam under earthquake loads is investigated. For this purpose, near-field earthquake records with focal depth lower than 15 km (for example Tabas earthquake 1978, M=7.4, and San Fernando earthquake 1970, M=6.6) are used. Moreover, to study the dam behavior under dynamic loads, interaction between concrete face and rockfill part of the dam is investigated and finally, some parameters including displacement, absorbed energy and base shear are evaluated. So, finite element method and Abaqus software is used for the study. Verification of the models is carried out using the results of previous researches by conducting modal analysis and determining natural vibration period. Then, the interaction between the concrete face and rockfill part as well as the effect of water level changes in stability of dam under dynamic load is investigated. Concrete behavior is simulated using concrete damaged plasticity. Therefore, concrete density, compressive strength and tensile strength and elasticity modulus are 2350 kg/m3, 25 MPa, 3 MPa and 29 GPa, respectively. Poisson’s ratio is assumed to be 0.2. Furthermore, 4-node shell elements are used to simulate concrete face and Drucker-Prager constitutive model is used to define rockfill material behavior.
The density and Poisson’s ratio for 2B, 3C and 3B layers are 2150 kg/m3 and 0.35, respectively. The shear modulus values for these layers are respectively 8.93, 2.89, and 3.85 GPa. In order to perform the simulation, the part of the dam structure beside the bed rock and the surrounding rock is considered as fixed bearing, and only the rockfill part and concrete face of the dam is simulated. Based on this assumption that the bed is rigid, there is no need to consider the dam foundation. This method is frequently used in literature review.
All the surfaces of the dam and bed rock are considered as fixed bearing to simulate the real condition where the dam is attached to bed rock and the surrounding rock. The interaction between dam layers is defined as tie. For defining the interaction between rockfill body and concrete face, tangential and normal contacts are defined using penalty method with friction coefficient equal to 0.5. In the next step, the model is meshed using 4-node shell elements for concrete face, 8-node brick and 4-node pyramid solid elements for rockfill body. Rayleigh damping is used to simulate the structure damping. The effective length of the dam reservoir has been determined by conducting several analyzes, so that the minimum required length for reservoir is reached in order to decrease the number of elements of the model.
Results and discussion
1. Interaction between concrete face and rockfill body
The results show that the increase of friction coefficient between concrete face and rockfill part from 0.5 to 0.7 has not affect the displacement of dam crown along the earthquake direction. However, when the concrete face is fixed to the rockfill part, significant changes are induced in dam crown displacement time history. In all cases, the deflection due to the dam weight is increased when the concrete face is attached to the rockfill body. The reason can be attributed to the tied interaction between these layers which results in similar deflection of concrete face with rockfill body and higher deflection of concrete dam crown. However, after the application of earthquake load, the displacement of the dam crown decreased in both analyses when tie interaction is defined between concrete face and rockfill body. In this study, due to the very high volume of analysis and its timeliness, it was not possible to examine the dam behavior in the free vibration regime, and therefore, it is not possible to assume the last displacement values at the end of analyses as the permanent displacement of dam. Figure 1 shows the relative displacement of the dam for the two selected earthquakes with a friction coefficient equal to 0.5 between the concrete face and the gravel body. According to Figure 1, the maximum displacement induced by the earthquake is related to Tabas and then, San Francisco earthquake. Furthermore, the high energy content of the Tabas record has been more effective in inducing greater displacement than the other record.
 
Figure 1. Lateral displacement of dam crown relative to the dam base for the selected earthquakes; Tabas and San Fernando.
The results also indicate that when the friction coefficient between concrete face and rockfill body is 0.5, the lowest damage occurs in the dam compared to that happens when friction coefficient is 0.7 or when the surfaces are tied. When the tied surfaces are used, the most damages takes place in concrete face, since all rockfill body displacement transmits to concrete face which results in much more concrete damages compared to the other interaction cases.
2. Effect of water level in reservoir on dam behavior
In this section, the effect of water level on seismic behavior of dam is investigated. For this purpose, the dam reservoir is analyzed in three cases including empty, half full and full (90% of dam height). Each study cases are examined under San Fernando and Tabas earthquakes. Figure 2 shows the relative displacement of dam crown in the three water level case for San Fernando and Tabas earthquakes.
 
Figure 2. Relative displacement of dam crown in three water level cases of empty, half and full for (a) Tabas and (b) San Fernando earthquakes
According to Figure 2, for both earthquakes, the dam crown displacement along the earthquake direction is significantly increased by increasing the water level, so that the maximum displacement in full case is 50% higher than empty case.
Conclusion
In this study, using the finite element method and simulation by Abaqus, the seismic behavior of concrete face rockfill dams has been investigated. For this purpose, the verification is firstly carried out using previous research results in literature. In the next step, nonlinear dynamical analysis is carried out, taking into account large displacements for the models under the earthquake record acceleration. The results illustrate that increasing the friction coefficient between the concrete face and the rockfill body from 0.5 to 0.7 has no significant effect on the displacement of the dam crown under earthquake load. Moreover, by using tie interaction between the concrete layer and the rockfill body, there is a substantial difference in the history of the relative displacement of the dam, and the displacement of the dam due to its weight has been increased. Furthermore, the results of this study exhibit that, with increasing the water level in dam reservoir, the deformation of the crown of the dam along the earthquake application direction has had a relatively significant increase, such that in the full state, the maximum displacement is increased by about 50% compared to that of the empty case. This is while the most damage of concrete is observed in the case when half height of dam in filled by water. Due to the more destructive power of near-field earthquakes and their impact nature, only near-fault earthquakes have been used in this research. Therefore, the results of this study are valid only for the behavior of dam under near-field earthquakes.
./files/site1/files/134/3.pdf
Mehdi Zare, Ali Moradzadeh, Abolghasem Kamkar-Rouhani, Faramarz Doulati-Ardejani,
Volume 14, Issue 3 (11-2020)
Abstract

Introduction
"Sulfide-carbonate" deposit is a term, which comprises a series of sulfide minerals such as Zn-Pb ore minerals, mainly considered as related to weathering of Zn-Pb sulfide concentrations and influence in sedimentary hosts (carbonate). There are more than 350 Zn-Pb deposits located in Iran, including world-class deposits such as Angouran, Mehdiabad and Irankouh. Due to the mining activity of these deposits, it creates a significant amount of mine waste that releases of these wastes in the environment causing severe problems. One of the main problems is the formation of Acid Mine Drainage (AMD). AMD is produced by oxidation of sulphide minerals, particularly pyrite (FeS2) in waste dump. Due to low pH and the ability to dissolve metals and other compounds, it can host a number of environmental problems. A phenomenon known as natural or alkaline mine drainage (NAMD) occurs at high pH values ​​when the neutralizing minerals are significantly present in the mine waste or when the oxidation of the sulfide minerals is poor. However, the metals and cationic species, such as Cu, Pb and Cd, are more soluble at low pH. In contrast, elements that form anionic species, such as Se, Cr, V, and Mo, tend to be more soluble at high pH and Ni, Zn, Co, As, and Sb, are soluble at near-neutral pH, and can potentially contaminate mine effluents, even without acidic conditions. Therefore Acid or Neutralization potential (AP&NP) of waste dump is significantly affects on the composition, transfers and fates of contaminations transmitted from waste dump. The aim of this study was to monitoring heavy metals concentrations and assessments of pollution potential of waste dumps in Anguran mine by static method and has been compared by mineralogical approach.
Material and methods
The Angouran Zn-Pb deposit is located in the 135 kilometers southwest part of Zanjan Province, NW Iran. This area belongs to the northwestern part of the Sanandaj-Sirjan Zone, a metamorphic belt related to the Zagros orogeny. Angouran mine is one of the most important carbonate hosted Zn-Pb deposits in Iran that mining activity has been created a significant amount of waste dump in around pit. To achieve the goals, the 47 samples taken from different surficial parts of the waste dump were analyzed by using the ICP-MS method to determine the concentration of elements and heavy metals. These elements and metals includes: Ca, Mg, S and As, Cd, Cr, Cu, Ni, Pb, Zn. The pollution index (PI) were modeled for heavy metal contamination risk zoning then modified Acid Base Accounting (ABA) static method was used to evaluate of acid and neutralization potential (AP&NP) of the waste dump samples and the results were modeled by Kriging method. At the end, mineralogical approach (Mg + Ca concentration) was used to determine the source of neutralization and to better interpret the static results.
Results and discussion
The results of contamination index showed that zinc, arsenic and cadmium had the highest average contamination index (18.89, 12.13 and 5.8, respectively) and the trend of total metal changes in the region as Zn> As> Cd> Pb > Ni> Cr> Cu was rated.
Datas measured in modified ABA method were modeled in 2D maps using the Kiriging method. Due to the low total sulfur content (less than 1%), all of the samples were Net Neutralization Potential (NNP) with a range of 49- 990 kg calcium carbonate per ton, and the study area was classified into three neutralization potential (NP), High (NP) and Very High (NP) levels. The mineralogical approach (Mg + Ca concentration) was used as a useful tool for better interpretation of modified ABA results and determines the neutralizing source. Mineralogical approach results indicate that calcite species are the main source of neutralization and have high correlation coefficient (R = 0.99) with the modified ABA method. In order to validate the results, the presence of mineral calcite was confirmed by XRD analysis on 4 samples.
Assessment of AP and NP of sulfide – carbonate waste dump in this research can be used as a basis model for other similar mines to control environmental problems and to identify the behavior and to transfer heavy metals in mine drainage in the future. Mineralogical approach results show that neutralizing potential and neutralizing source can be obtained without using expensive mineralogy analyses in this type of carbonate-sulfide deposit



 
Moslem Babaei, Ali Raeesi Estabragh, Jamal Abdollahi, Mohadeseh Amini, Gholamali Vakili,
Volume 14, Issue 3 (11-2020)
Abstract

Introduction
Expansive soils are a very common cause of extreme damages because they are susceptible to volume change due to a change in water content. Geotechnical problems associated with the expansive soils are well documented in different literature. As a result, a clear understanding of the behavior of such soils is required for the effective design of structures and infrastructures on these soils. The effects of hydrocarbon pollutants as a flooding fluid on the swelling potential of an expansive soil during wetting and drying cycles have not been considered in the previous researches. The aim of this research is to study the properties of an expansive soil with different flooding fluids, i.e. distilled water and solutions of glycerol with 10 and 20% through a number of cycles of wetting and drying tests under constant surcharge pressure.
Material and methods
The soil that was used in this work was a highly expansive clay soil (according to the classification by McKeen (1992)). It was prepared by mixing 20% bentonite and 80% kaolin. This soil was classified as a clay with high plasticity according to the Unified Soil Classification System (USCS). The optimum water content in the standard compaction test was 18.11% and the maximum dry unit weight was 16.27 kN/m3.
Distilled water and solutions of glycerol with concentrations of 10 and 20% were used for flooding the samples. To prepare the glycerol solutions, the required amount of glycerol was mixed with distilled water.
For making compacted samples for testing, the needed air-dried soil was weighed and the required water was added to it to reach the desired water content (4% below the optimum water content according to the compaction curve). The soil and water were mixed by hand and then was kept in a plastic bag for 24 hours to allow the uniform distribution of moisture in the soil. Samples were prepared by static compaction of the moist soil in a special mould.
A conventional oedometer was modified to allow the wetting and drying tests to be conducted under controlled surcharge pressure and temperature. During wetting and drying, the vertical deformation of the sample was measured by using a dial gauge. The variation of water content with void ratio during wetting and drying cycles was determined by using the information from the duplicated samples.
Results and discussion
Fig. 1 shows the variations of vertical deformation during wetting and drying cycles for samples that were flooded with distilled water and solutions of 10 and 20% glycerol. This figure illustrates that by increasing the number of cycles the amount of irreversible deformation is reduced until the equilibrium condition is achieved where the deformation due to wetting and drying is nearly the same. These results indicate that by increasing the concentration of glycerol the equilibrium condition with reversible deformation is reached in a fewer cycle of wetting and drying than the sample that was flooded with distilled water.

Figure 1. Wetting and drying cycles for different quality of flooding fluids
The results of void ratio versus water content at the equilibrium conditions for the samples flooded with distilled water and solutions of 10 and 20% glycerol (that were obtained from duplicated samples) are shown in Fig. 2. This figure displays that the paths of drying-wetting for different flooding fluids are nearly S-shaped curves. It is also seen in this figure that the order of the curves in this space is dependent on the percent of glycerol, the curves for the sample flooded with distilled water and 20% glycerol are located at the top and bottom of the space of void ratio against water content.

Figure 2. Water content-void ratio paths for different quality of flooding fluids
The change in the thickness of the diffuse double layer (DDL) affects on the swelling behavior of soil. The thickness of DDL is dependent on factors such as valency and concentration of cations, temperature, and dielectric constant. The value of dielectric constant for water is 80 and for solutions of 10 and 20% glycerol are 74.9 and 71.8, respectively. The magnitude of the attractive and repulsive forces between clay particles are inversely and directly depended on the value of the dielectric constant. The reduction in the value of the dielectric constant causes an increase in the attractive forces and leads to a reduction in the thickness of DDL. When the flooding fluid is a solution of glycerol, the initial chemical composition of pore fluid in the sample is changed. The chemical composition of pore fluid has different effects on the structure of clay soil such as changes in the thickness of DDL. When the flooding fluid is distilled water the pore fluid of samples has a dielectric constant of about 80. Therefore, the values of attractive and repulsive forces are not changed because of the same dielectric constant of flooding fluid and pore fluid. The results of tests on these samples (flooded with distilled water) show that by repeating the wetting and drying cycles the potential of swelling is reduced and after several cycles a reversible equilibrium condition is attained as depicted in Fig.1. When the pore fluid is the solution of glycerol, the attractive forces are increased due to the reduction of the dielectric constant of pore fluid and causes a reduction in the thickness of DDL. The shrinking of DDL is led to the formation of flocculated structure in the soil and results in pasting of particles together leading to the reduction potential of swelling. When the concentration solution of glycerol is increased the dielectric constant is decreased, the magnitude of attractive forces is increased and the degree of flocculation of the soil structure is increased that is yielded to a reduction of swelling potential.
Conclusion
Effect of different flooding fluids on the properties of an expansive soil during wetting and drying cycles were studied. The following conclusions can be drawn from the present research:
-After a number of wetting and drying cycles, the observed irreversible          deformation was diminished and equilibrium was achieved. The solution of glycerol causes more reduction in the potential of swelling than distilled water.
-The wetting and drying paths in the space of void ratio and water content are S-shaped curves. The variations in the void ratio of samples flooded with the solution of glycerol are smaller than distilled water../files/site1/files/142/babaei.pdf
 
 
Ali M. Rajabi, Shima Bakhshi Ardakani,
Volume 14, Issue 4 (12-2020)
Abstract

Introduction
Improving the geotechnical characteristics of soils including superficial or deep soils has always been a challenge to geotechnical engineers. Therefore, various physical and chemical methods are used to improve different types of soils. In general, any physical, chemical, biological or combination of methods are used to change the characteristics of natural soil mass in order to achieve engineering goals which is defined in the "soil stabilization." Among different types of additives for soil stabilization, the use of pozzolans has been investigated by researchers because of their chemical compatibility with the environment and the cementation products due to chemical reactions. Todays, a lot of researches has been done on the use of natural or artificial zeolites as pozzolanic materials for the production of cement mixtures. This material, as a pozzolan, increases the speed of the pozzolanic reactions and reduces the density of cement products. However, many studies have been done to investigate the effect of zeolite and sepiolite on the resistance of cement products such as concrete, but so far, the use of these additives has been less considered for soil improvement. On the other hand, because of the compatibility of zeolite and sepiolite with the environment and their unique physiochemical properties, it is necessary to pay attention to these additives in order to improve the soil. Therefore, in this research, the effect of zeolite and sepillot additives with different percentages at different treatment times have been investigated to determine the elasticity modulus and hydraulic conductivity with focus on soil microstructure behavior.
Materials and methods
1. The properties of the soils
In this research, two types of soil including clayey sand (with 20% clay) and sandy clay (with 51% of clay) were used. The studied soils were a mixture of clay and sand of Firoozkouh (a typical type of sand located in north of Iran). Some physiochemical properties of zeolite and sepiolite are presented in Table 1.
Table 1. Physiochemical properties of zeolite and sepiolite used in this study
L.O.I. Na2O K2O MgO CaO Fe2O3 Al2O3 SiO2   Chemical component
25.11 0.02 0.01 15.73 0.01 o.61 0.3 55.3   Sepiolite (%)S
11.94 0.13 - 0.87 2.45 1.26 13.54 69.74   Zeolite (%)
2. Experiments
The uniaxial compressive strength tests were performed at 0.1 mm/min according to ASTM D2166 standard. The stabilized soil samples were compacted at percentages of 0, 5, 10, 15, 20 and 25 in cylindrical molds (38mm × 76mm) in five layers to achieve the desired density. In order to investigate the effect of curing time, the samples were placed inside sealed containers and underwent the test at instantaneous, 7, 14, and 28 days and at the desired additive percentages. To investigate the effect of additives on the soil hydraulic conductivity, clayey sand soil with additives 5, 10, 15, 20, and 25% was prepared using dry mixing method. Then, the prepared mixture was poured from a specific height into the permeability mold with a height of 8.65 cm and diameter of 5 cm. In this way, the specific dry unit weight of all samples was obtained as 1.47 g/cm3, close to the minimum specific dry unit weight. In this research, concerning the considerable effect of fine-grained soils on hydraulic conductivity, falling head test was used to determine the permeability coefficient.
In order to the morphology of the clayey sand soil without additives and stabilized with additives 15% was examined through SEM test.
Discussion and results
1. Modulus of elasticity
In this study, after uniaxial tests in different percentages and ages, the stress-strain graphs were plotted and then the elasticity modulus was calculated. The results showed that, with increasing zeolite content, the modulus of elasticity has been increased and, with increasing curing time, except for a slight decrease, after 7 days, the modulus of elasticity increased. During the initial treatment (7 days), the hardness of the sandy clay soil decreased and then increased with increasing time. In general, hardness in both soils in the high percentages of zeolite is significantly is increased.
Also, the effect of sepiolite on the modulus of elasticity has been studied. The results indicate that with the increase in the percentage of additive and lengthening the curing time, the modulus of elasticity is increased. This increase in the stabilization of both sandy clay and clayey sand soil is almost the same. In addition, in the case of sepiolite modification, the elasticity of sandy clay and clayey sand is approximately equal to 5 times in comparison to the initial value of unstabilized soil. However, in zeolite, the modulus of elasticity in clayey sand soils is almost 2 times, and sandy clay is nearly 5 times higher.
2. Permeability
To investigate the effect of additives on the soil hydraulic conductivity, clayey sand soil with additives 5, 10, 15, 20, and 25% was prepared using dry mixing method. The samples were saturated in a short period and permeability test was carried out immediately. Permeability coefficient changes were mostly influenced by physical factors. Therefore, due to the fineness of both types of additives, the hydraulic conductivity decreases with increasing additive content. The amount of reduced hydraulic conductivity in sepiolite stabilization is greater than zeolite due to the structure of the sepiolite (fiber-shaped) compared to zeolite.
 
 
3. SEM imaging
In this study, attempts were made to examine the reasons behind the obtained results more carefully through SEM imaging.

c                                     b                              a
Figure 1. SEM image of non-stabilized clayey sand soil (a) soil stabilized with zeolite 15% (b) soil stabilized with sepiolite 15% (c) during the curing time of 28 days at magnifications 10000X
Figure 1a displays the SEM image of non-stabilized clayey sand soil. As can be seen in the figure, the soil structure is clear as layered and clay scales can be seen as laminated. Figure 1b demonstrates the SEM images of clayey sand soil stabilized with zeolite 15% during the curing time of 28 days. The sample has lost its layered structure in response to stabilization with zeolite during the curing time and changed into an integrated structure. This can be due to incidence of chemical reactions such as ion exchange and pozzolanic reactions in response to adding zeolite. Figure 1c demonstrates the SEM images of clayey sand soil stabilized with sepiolite 15% during the curing time of 28 days. As shown in the figure, the sepiolite has a fibrous-shaped structure that is longitudinally twisted. Also, with  curing time increase, complex structures have emerged that could be due to the occurrence of chemical reactions.
Conclusion
This study examined the effect of zeolite and sepiolite additives on strength parameter of clayey soils. Accordingly, uniaxial compressive strength test was performed on clayey sand and sandy clay soil at percentages of 0, 5, 10, 15, 20 and 25% of zeolite and sepiolite with instantaneous curing times of 7, 14 and 28 days. Further, permeability test was conducted at different percentages on stabilized clayey sand soil. Also, to investigate the effect of these materials on soil microstructure, SEM imaging was performed at 28 days. The results show that both additives increase the elastic modulus of clayey sand and sandy clay soils. Also, the results indicate a steady increase in the stiffness of the cured soil with sepiolite during processing time. However, reducing soil hardness can be seen in stabilizing with zeolite at lower rates and lower percentages. In permeability test, hydraulic conductivity decreases with increasing additive content. The rate of permeability reduction in sepiolite is higher than zeolite. SEM images show that chemical reactions create an integrated structure that ultimately increases uniaxial compressive strength and modulus of elasticity. Also, SEM imaging depicts physical changes along chemical reaction in soil stabilized with sepiolite. Ultimately, increasing soil strength resulting from additive alongside environmentally friendliness is recommended in superficial and deep improvement of soil../files/site1/files/144/Rajabi.pdf
 
Maziar Hosseini, Majid Taromi, Mahdi Saeidi, Vahid Soleimani, Mehdi Soltani Negar,
Volume 14, Issue 4 (12-2020)
Abstract

Introduction
Series A of coarse-grained alluvial deposits of Tehran are extended in eastern and north-eastern areas of Tehran. Analyzing and studying of these alluvial deposits from a geological point of view as well as their creation time and general characteristics such as the deposits’ mineral types, their source, and formation conditions, gives a better point of view to geotechnical engineers about exploring their characteristics as well as geotechnical aspects in underground structure design, excavations, and foundation design processes. On the other hand, in order to analyze stability, estimating the factor of safety and the seismic design of these structures, considering their location, which is in Tehran with a high seismic hazard area, the necessity of knowing the exact mechanical and dynamic properties of Tehran's alluvium is felt more than ever.
Material and methods
Due to the grain size of Tehran’s coarse-grained alluviums (series A) as well as high level of cementation of them, it is impossible (or maybe so difficult) to make undisturbed samples in order to do experiments. Such that it is excavated 23 boreholes with 30 to 140 meters depth as well as 17 test wells with 20 meters depth in an area which was extended in 10 kilometers in long which were located in Tehran’s No. 13 and No. 14 districts (as it can be seen in Figure 1). During the excavation of the entrance ramp and tunnel of eastern highway of Tehran, in-situ tests have been done in different sequences. Since it was important to investigate real behavior of these alluviums, different in-situ tests such as plate load test, in-situ shear test, pressuremeter test, and downhole test have been done as well as many laboratory and field tests. Furthermore, (1) X-Ray Diffraction (XRD) and (2) X-ray Fluorescence (XRF) as well as (3) Scanning Electron Microscopy (SEM) methods, have been used to explore the type of minerals and those used in cementation.
 
 
 
 

(ب)
 
 
 
Figure 1. a) Geological plan and the location of boreholes and test wells in the alignment of East Tehran Freeway
Results and discussion
Based on the results of XRD tests, it is quite clear that the largest weight percentages of tested samples are lime and silica.
Calcium and magnesium levels-as the high-power cations in flocculation process-in soil sample No. 1 (soil with high cementation level) are much more than soil sample No. 2 (soil with moderate cementation level).
This is the cause of high cementation level of soil sample No. 1 comparing with soils sample No. 2. A rapid increase in stress level can be seen in in-situ shear test results, in low shear displacements, up to reaching a maximum of τp (peak point) and afterwards reduction in shear stress with softening behavior.  
Cohesion and shear strength levels also increase by increasing the depth. According to the plate load tests results, an increase in soil modules changes can be seen in different depths by depth increasing.
Large tendencies to increase in volume and dilation can be seen in under shear load cemented soils, after applying a primary compression on them. A brittle behavior with the occurrence of a certain peak can be seen in cemented samples. The significant increase in strength is directly related to the severe dilation rate, which can be seen in cemented samples results.  The shear strength would be decreased, if this cement is broken during the particles’ displacements.
The results of downhole tests are shown in Figure 2. According to this figure, it has been explored that Vs,30 is about 600 m/s in moderate cemented soils while it is about 850 m/s in highly cemented soils.  Because of the homogeneity and uniformity of sedimentary deposits, shear wave velocity is increasing due to the higher density of the layers and high level of cementation in both of the soil types. However, this increase is not significant at depths above 25 meters.
Conclusion
Based on the results, cementation level of the eastern coarse-grain-alluvium of Tehran is moderate to high and minerals used in cementation of this type of soil are generally carbonated and especially calcite.
Investigating the level of cementation of soil as well as the results of chemical analysis and in-situ tests, it can be found that the strength and deformation parameters of the soil are directly related to the degree of its cementation.
Based on the obtained results, the deformation modulus increases by about 25%, the cohesion by about 55% and the shear wave velocity by about 30% with increasing the degree of cementation (Table 1).
Increases of these parameters are directly related to depth. However, the cementation level does not significantly affect the internal friction angle of the soil.
Table 1. Average results of in-situ shear tests
Deformation Modulus (MPa) Peak Friction Angle (deg.) Cohesion
(kPa)
USCS Depth
(m)
Sample
50-60 39 30-35 GW-GM 5 Moderately Cemented Soil
(M.C. Soil)
75-85 41 50-60 SP-SC 10
85-90 41 50-60 GW-GC 15
95-105 41 50-60 GW-GC 20
60-70 39 35-40 GW-GM 5 Highly Cemented Soil (H.C. Soil)
75-85 39 50-60 GW-GC 10
110-120 42 65-75 GW-GC 15
125-140 41 110-120 GC 20
 
 
, , ,
Volume 15, Issue 3 (12-2021)
Abstract

In this research work the effect of glycerol and acetone on the properties of a clay soil with low plasticity (CL) was studied through experimental tests. Contaminated soil was prepared by mixing 10, 15 and 20% pure glycerol or acetone with soil. A set of experimental tests including Atterberg limits, compaction, free swelling, Unconfined Compressive Strength (UCS), California Bearing Ratio (CBR), consolidation tests and Scanning Electron Microscopy (SEM) were conducted on uncontaminated and contaminated soil. The results showed that the effect of glycerol on contaminated soil was a reduction of Atterberg Limits, Free Swelling, optimum water content and increasing maximum dry unit weight but the effect of acetone was in the opposite trend of glycerol. These variations in the behaviors of soil were the function of percent of glycerol or acetone. In addition, the results indicated that both glycerol and acetone cause the reduction in the values of compressive strength and CBR number of soil but the reduction due to glycerol was more than the acetone. The results of consolidation test showed that the values of compression index (Cc) and swelling index (Cs) were independent of the percent of chemical fluids. The SEM results also reveal that these two contaminations cause some variations in the physical and mechanical properties of a clay soil becase of changing the structure of soil and the effect of these contaminations on the properties of soil was not the same.

Key words: hydrocarbon contamination, Atterberg limits, free swelling, Unconfined Compressive Strength (UCS), CBR and consolidation test.

./files/site1/files/%D8%A7%D9%81%D8%B3%D8%B1%DB%8C_%D8%B1%D8%A6%DB%8C%D8%B3%DB%8C(1).pdf
, , , ,
Volume 15, Issue 4 (12-2021)
Abstract
, , , , ,
Volume 15, Issue 4 (12-2021)
Abstract
Ali Kamali, Dr Mehdi Mokhberi, Dr Abbas Ghalandarzade,
Volume 16, Issue 2 (9-2022)
Abstract

 Marls are one of the problematic soils that undergo more erosion due to their deformability and sensitivity to humidity and weather conditions. The effects of these soils have been observed on subgrades of pavements, abutment of dams, foundation of high-rise structures, the interaction of soil-structure and etc. The strength and deformation of marls are more effective to moisture content. In addition, the dynamic parameters of theses soil are considerable as well as static geotechnical specifications. Soil dynamic parameters can be obtained from both laboratory experiments and field experiments. With the aim of understanding the dynamic behavior of marl soils, this study was carried out to evaluate the marls of northwest region of Shiraz City which is mixture of marls containing the expansive anhydride gypsum. For this purpose, the laboratory cyclic triaxial test and in situ downhole test has been performed. The results showed that for normal consolidated marls, with increasing the confined pressure from 400 to 600 kPa (increasing soil depth), the shear modulus increases from 50 to 200 kPa. In addition, the behavior of the damping ratio is relatively different for strains less than and greater than 1%, but in general, for strains higher than 1%, the damping ratio decreases from 0.21 to 0.18 with the increase confined pressure from 100 to 600 kPa.
 

Mr. Farhad Mollaei, Dr. Reza Mohebian, Dr. Ali Moradzadeh,
Volume 18, Issue 3 (12-2024)
Abstract

The brittlenessindex is one of the most important parameters in geomechanical analysis and modeling. Many methods have been proposed to estimate the brittleness index. One of the recently used methods is the  intelligent method. In this paper, firstly the aim is to introduce a new algorithm using deep learning algorithms to predict the brittleness index in one of the wells of the hydrocarbon field in southwest Iran. In this article, first, the effective features for the input of the algorithms were determined using Pearson's correlation coefficient, and then using (recurrent neural network + multi-layer perceptron neural network) (LSTM + MLP) and (convolutional neural network + recurrent neural network) (CNN+ LSTM) brittleness index was estimated and the mean error value (MSE) and coefficient of determination (R2) were calculated for the training and test data. For both training and test data, both algorithms have a coefficient of determination close to 1 and a very low error. Also, in order to ensure the results of the algorithms, a part of the data was set aside as blind data, and the error and coefficient of determination were calculated for this data, and the error was MSE CNN+LSTM =26.0425,  MSE LSTM+MLP =32.0751  and the coefficient of determination was R2 CNN+LSTM  =0.8064,  R2 LSTM+MLP  =0.7615 . The results show the effectiveness of the introduced deep learning algorithms as a new method in predicting the brittleness index, and comparing the two algorithms presented, the CNN+LSTM algorithm has higher accuracy and less error.

Dr Eisa Hajiradkouchak, Dr Behzad Rahnama, Dr Hasan Nasrollahzadeh, Mr Ali Shahbazi, Mr Reza Raeiji, Mr Kazem Babaei,
Volume 18, Issue 3 (12-2024)
Abstract

Many researchers believe that providing safe water, sanitary disposal and optimal management are the three axes of health, and in all these cases, while paying attention to the process of doing work, continuous control should also be done. This study was designed and implemented with the aim of seasonally investigating the physicochemical and microbial water quality of Qarasu River in Golestan province using the IRWQIsc index. 6 sampling stations were identified for Qarasu River and sampling was done once every month in four seasons of 1400. The measured parameters include pH, BOD, COD, dissolved oxygen (DO), electrical conductivity (EC), ammonium (NH4), nitrate (NO3), phosphate (PO4), total hardness (TH), turbidity and total suspended solids. It was a stool form. According to the measured parameters, Iran's surface water quality index IRWQISC was calculated. The results of the study based on the index showed that the quality of this index for all stations in all seasons was between 70.5 and 14.7 and according to the IRWQISC index, it was in five good categories (70-1.85), relatively good. (55-1/70), relatively bad (30-44-9), bad (15-29-9) and very bad (less than 15). The influencing parameters were total suspended solids, turbidity, nitrate, temperature and fecal coliform. It can be concluded that the amount of 70.5 with good quality is related to (Tuskestan village) in winter and the amount of 7.14 with very bad quality is related to (Pol Qara Tepe) in summer that the quality of the river water in The Gorgan to Aqqla road bridge station (Qorban Abad) is in bad condition in all seasons due to the entry of urban and industrial pollutants into this station, and Tuskestan village station has good and relatively good quality in most seasons because Tuskestan is in It is located in high altitudeand the entrance of clean running water  into thisarea is more and it is far from industrial and urban pollutants.

Soroush Mahdavian, Navid Rashidi, Ali Raeesi, Jamal Abdullahi,
Volume 19, Issue 1 (6-2025)
Abstract

Clay soils typically have low strength and a high swelling percentage. They are considered as problematic soils in Civil Engineering projects. This research study examined the effects of magnesium chloride (MgCl2) solution on the  clay soil  improvement  through conducting laboratory experiments. The experimental program included Atterberg limits, compaction, swelling, unconfined compression strength (UCS) and Scanning Electron Microscopy (SEM) tests. Available clay soil in the Lab was mixed with MgCl2 solution at weight percentages of 3%, 5%, 7% and 10%  Samples for the swelling and strength tests were made using thestatic compaction method. The moisture and dry unit weight of the prepared samples were the same as those of thecorresponding compaction curves. The strength test results showed that the final strengths of the samples with 3% MgCl₂ at 7-, 14-, and 28-day curing times were 1401, 2018, and 1848 kPa, respectively. The results also showed that a reduction in strength of the samples occurred with more than a 3% solution of MgCl₂. For samples with 10% MgCl2 solution, the strength decreased until 14 days of curing time, but increased thereafter. Additionally, the results indicated that the reduction in swelling percentage compared to natural soil was 4.95%, 3.98%, 2.8%, and 3.9% for samples with 3%, 5%, 7%, and 10% MgCl₂, respectively, showing that the reduction in swelling depends on the MgCl₂ percentage. Additionally, the SEM results showed that the improvement in the soil was due to chemical reactions between the soil and MgCl₂.

Maedeh Roshan Liarajdameh, Milad Davari Sarem,
Volume 19, Issue 2 (10-2025)
Abstract

Iran, due to its location between two active tectonic plates, has always been exposed to numerous earthquakes. The occurrence of more than 100 severe earthquakes in the past century indicates the country’s high level of vulnerability to this natural hazard. The aim of this research is to analyze the seismicity and assess the earthquake hazard in Shahid Rajaei Port, as the largest commercial port in Iran (located at the intersection of the North-South transit corridor), which will be a fundamental step in enhancing the resilience and sustainability of the vital infrastructures in this region. In this study, all seismic events occurring within a 200-kilometer radius of the site were used, along with the Knopoff and Ez-Frisk software. The statistical analysis of historical and instrumental earthquakes indicates a high level of seismicity in the region, characterized by moderate-magnitude earthquakes with short return periods, such that earthquakes with magnitudes between four and five on the Richter scale constitute a larger share. The probabilistic hazard assessment estimated the maximum horizontal and vertical accelerations as 0.385 and 0.290 (g), respectively. Additionally, the site response spectrum was prepared based on the accelerographs of the Tabas earthquake and the isoacceleration maps of the study area, generated at intervals of 1.0 degrees in both latitude and longitude directions. The results showed that the study area has a seismic hazard of 0.85 (g), which is higher than the standard values specified in Iran’s Code 2800, placing it within the very high relative hazard zone. Therefore, implementing risk-based approaches in infrastructure development helps optimize port design and reduce earthquake-related damages.
 
Majid Taromi, Abbas Eftekhari, Jafar Khademi Hamidi, Mahdi Saeidi,
Volume 19, Issue 3 (12-2025)
Abstract

The Sabzkuh Water Conveyance System was planned as an inter-basin water transfer project in the Sabzkuh region of the Chaharmahal and Bakhtiari province, a tributary of the Karun River. Initially, the project included a diversion dam, an open canal, and a tunnel. However, in the initial feasibility studies, the potential environmental risks to the Choghakhor lagoon, ecosystem and regional groundwater resources were not given sufficient consideration.Only due to land acquisition problems in the Chahartagh forest area, the tunnel route was changed and the channel was also removed. Consequently, the project was scaled down to comprise only a diversion dam and a 11-kilometre-long tunnel. This is while the geological and environmental risks of the project area were not fully investigated in the selection of the final route. After 484 meters of excavation, the Solaqan fault was encountered and a sudden inrush of water with a flow rate of 100 liters per second and the liquefaction of the excavation front caused the TBM to get stuck. Subsequently, the drainage of water caused the drying of the Aloqareh spring and the depletion of the Shirkoshteh spring, and negative environmental and social consequences for the region. To release the TBM, two side galleries were constructed and grouting and drainage were carried out. Ultimately, the tunnel excavation operation was stopped due to the problems that arose. The challenges that arose in this project indicate the need for a comprehensive study and planning in the selection of the route and TBM, as well as environmental, economic and social assessments at all stages and examining their interaction with each other.

Soroush Mahdavian, Ali Raeesi Estabragh, Shima Azadeh Ranjbar,
Volume 19, Issue 3 (12-2025)
Abstract

This research study investigated the impact of dimethyl phthalate (DMP) on the physical and mechanical properties of clay soil through experimental testing. Additionally, the impact of hydrated lime and magnesium oxide on improving the properties of clay soil was investigated. The contaminated soil was artificially produced in the laboratory. Natural and contaminated soils were mixed with the above agents at percentages of 5.0%, 10.0% and 15.0%, respectively. The experimental test programme for this study comprised: Atterberg limits, compaction, uniaxial compression tests (UCS) and scanning electron microscopy (SEM). Samples for the UCS tests, both with and without additives, were prepared using the static compaction method and tested at curing times of 7, 14 and 28 days. The results showed that, in general, the Atterberg limit, compaction parameter (maximum dry unit weight) and UCS values for the contaminated soil were lower than for the natural soil. The results also indicated a reduction in UCS values for a mixture of contaminated soil and 10% hydrated lime or magnesium oxide. A comparison of the final strength values of samples made from a mixture of contaminated soil and 10% hydrated lime or magnesium oxide showed that the strength was 34.4% and 63.8% lower than that of a mixture with 5% of these additives at the same curing time. The E50 values were calculated from the stress-strain curves of the different tests. Additionally, the SEM results showed that changes in the properties of the mixture prepared from these agents and natural or contaminated soil were due to a chemical reaction between the soil and the additives.

Reza Mohseni Afkham, Dr Mojtaba Bahaaddini, Dr Abbas Majdi,
Volume 19, Issue 3 (12-2025)
Abstract

Tensile strength is one of the most important mechanical properties of brittle materials and plays a decisive role in the stability of many civil and mining structures. The Brazilian test is the most common indirect method for determining tensile strength and is widely employed. In this test, it is generally assumed that a central tensile crack initiates and propagates along the loading axis. However, the actual fracture process in the Brazilian test remains a controversial issue, and using curved loading platens has been recommended to better concentrate tensile stresses at the center of specimen. This study investigated the influence of platen curvature on the estimated tensile strength and the fracture patterns. To this end, five types of platen with curvature ratios of 0, 0.50, 0.57, 0.67, and 0.80 were prepared. All tests were recorded using a high-speed camera to precisely capture the initiation and propagation of cracks. To minimize the effect of rock heterogeneity and obtain consistent results, synthetic specimens were used, and five samples were tested for each curvature ratio.The results indicated that increasing the platen curvature led to a higher estimated tensile strength. While the increase was negligible for curvature ratios up to 0.67, at the ratio of 0.80 the tensile strength was approximately 48% higher compared to 0.67. Analysis of fracture patterns revealed that at the curvature ratio of 0.80, the fracture mode shifted to an unstable and disturbed pattern, characterized by secondary shear cracks and the irregular propagation of the main crack.

Ms Haniye Yaghoubi, Dr. Reza Jahanshahi, Dr. Morteza Mozafari,
Volume 19, Issue 3 (12-2025)
Abstract

This study examines the hydrochemistry and contamination levels of groundwater resources in the urban area of Birjand in eastern Iran. Water quality was assessed and pollution sources were identified through sampling 22 wells, 12 qanats and 4 springs. The results showed that electrical conductivity varied from 300 to 8,000 µS/cm, while pH ranged from 7.23 to 8.71. According to the Piper diagram, the dominant hydrochemical facies were chloride, sulphate and bicarbonate types. In some of the samples, the nitrate concentration exceeded the permissible limit of 50 mg/L set by the World Health Organization, indicating the influence of urban wastewater and agricultural effluents. The ionic ratios reveal the influence of halite and gypsum dissolution processes, as well as ion exchange reactions, on the chemical composition of the water. A health risk assessment showed that, while most sources are within the safe range for adults, some wells and qanats pose a higher risk to infants and children. This study aims to provide a scientific framework for understanding the geochemical processes that control water quality, and to identify high-risk areas for the sustainable planning and management of groundwater resources in the Birjand plain.

Faeze Samadpoor, Morteza Mozafari, Majid Dashti Barmaki, Parisa Sharifi,
Volume 19, Issue 4 (12-2025)
Abstract

Groundwater plays a vital role in meeting the drinking and agricultural water needs of Kermanshah Plain. In order to protect the aquifer, it is important to evaluate its sustainability in the face of current and future demands and stresses. Groundwater sustainability indicators help ensure the sustainable management of these resources. This research aims to evaluate the sustainability of groundwater resources in the Kermanshah Plain using various indicators. To this end, AHP analysis was used to evaluate the sustainability indicator of this aquifer based on nine indicators in five quantitative, qualitative, environmental, social, and political sectors. First, the value of each indicator was calculated, and then its sustainability was evaluated using data transferred to GIS software and interpolation. Next, the weight and rank of each indicator and category were calculated to prepare an index-equivalent map. Then, using weighted overlap, the final sustainability map was obtained. Finally, the Receiver Operating Characteristic (ROC) curve was used to measure the accuracy of the results. The prepared sustainability map shows that indicators of groundwater storage changes and quality conditions are among the most important factors affecting the sustainability of the plain's groundwater resources. The results also show that the sustainability situation is weaker in the central areas and more favorable in the border areas (river headwaters) and southeast of the aquifer. To improve the sustainability of the region's groundwater resources, it is recommended that new water management policies be adopted with the participation of the people and based on scientific, principled solutions.


Page 2 from 2     

© 2026 CC BY-NC 4.0 | Journal of Engineering Geology

Designed & Developed by : Yektaweb