Next Article in Journal
Spatio-Temporal Variability of Soil Properties and Nutrient Uptake for Sustainable Intensification of Rainfed Pigeon Pea (Cajanus cajana) in Semi-Arid Tropics of India
Next Article in Special Issue
Finite Element Analysis for the Mechanism of Stress Wave Propagation and Crack Extension Due to Blasting of a Frozen Rock Mass
Previous Article in Journal
Sense of Competence and Feelings of Stress of Higher Education Faculty in the Transition to Remote Teaching: What Can We Learn from COVID-19 Pandemic in the Long Run
Previous Article in Special Issue
Experimental Study on Shear Strength Parameters of Round Gravel Soils in Plateau Alluvial-Lacustrine Deposits and Its Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analyzing Geotechnical Characteristics of Soils in Erbil via GIS and ANNs

by
Zhvan Baqi Qader
1,*,
Zuheir Karabash
2 and
Ali Firat Cabalar
1
1
Civil Engineering Department, University of Gaziantep, Gaziantep 27310, Turkey
2
Department of Dams and Water Resources Engineering, College of Engineering, University of Mosul, Mosul 41002, Iraq
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(5), 4030; https://doi.org/10.3390/su15054030
Submission received: 12 January 2023 / Revised: 14 February 2023 / Accepted: 18 February 2023 / Published: 22 February 2023
(This article belongs to the Special Issue Advances in Rock Mechanics and Geotechnical Engineering)

Abstract

:
The aim of this study is to analyze and model the geotechnical characteristics of soils in Erbil city using Geographic Information Systems (GIS) and Artificial Neural Networks (ANNs). The study used GIS to analyze the geotechnical properties of soils by collecting data from 102 boreholes in three different depth levels (1.5 m–3.5 m, 3.5 m–6.5 m and 6.5 m–9.5 m) to visualize and analyze soil characteristics such as fines content, moisture content, soil plasticity, shear strength parameters, compressibility, Standard penetration test (SPT), and bearing capacity. The paper also establishes the prediction of SPT-N value and bearing capacity based on geotechnical properties of soils using ANN methods and made correlations between SPT values and shear strength parameters with the bearing capacity of the soil. The results analyzed via GIS indicated that the soil classification was silty clay with a small amount of sandy gravel (CL) in most of the study area. According to the SPT–N values, most of the soils in Erbil City ranged between 33 and 50; a higher SPT value generally indicates denser and stronger soil. The value of the shear strength parameter for the maximum friction angle of the soil layers was found to be 36°, and the predominant cohesion was approximately 100 kPa. The compression index of soils ranged between 0.11 to 0.31. The results showed that the ANN models were able to accurately predict the geotechnical parameters of the soil types in the study area. In addition, the use of GIS and ANN techniques allowed for a comprehensive analysis of the geotechnical characteristics of the soils in Erbil, providing valuable information for future construction and development projects.

1. Introduction

One of the most important steps before constructing infrastructure is the geotechnical site investigation. It provides information on the site suitability for design criteria and possible construction problems such as time and resources. There are many methods for site investigations and in-situ tests, including pressure meter test, dilatometer test, SPT, cone penetration test (CPT), and plate load test [1]. For the construction of multistory buildings, highways, bridges, and industrial facilities, a soil survey is required to determine the type of soil, consistency, index properties, relative density, groundwater level, shear strength parameters, (SPT) value and bearing capacity [2]. It is necessary to know the bearing capacity of the soil layers for design, the choice of the foundation type, and the foundation depth for any superstructure [2,3,4]. A geotechnical investigation provides valuable information on the physical and mechanical properties of the soil and rock at a site, which is necessary for safe and durable engineering structures. The information collected from a geotechnical investigation is used to make informed decisions about the design and construction of the foundation, and to identify potential hazards, such as soil liquefaction or instability of slopes, that could compromise the safety of the structure. Therefore, to obtain the geotechnical parameters required for the calculation of the soil bearing capacity and settlement, in situ testing is required in addition to the collecting of disturbed and undisturbed specimens at different depths. Thus, several geotechnical experiments are conducted on these specimens to determine various parameters that are typically used to design the foundations [5]. Researchers have studied the reliability of SPT to determine the bearing capacity of soil [6,7,8]. Currently, the SPT test is used to evaluate the bearing capacity to design foundations [9,10].
GISs and ANNs could be used together in analyzing the geotechnical characteristics, and to predict the shear strength, settlement, and bearing capacity of the soil from the index properties of soils. The GISs provide an analytical function that is time–consuming for developing model entry data at different spatial scales [10,11]. A GIS is an organization of data that people interact with to integrate, analyze, and visualize data, to identify relationships, patterns and trends, and to resolve complicated issues; GIS has been used by many researchers to analyze various data [12,13,14,15]. ArcGIS was designed to capture data, store, update, process and present data, and to conduct analyses [16]. GIS can help to recognize possible challenges to the completion of the project early in the design process, which can help to avoid time losses for a construction project. Therefore, a GIS is a modular instrument that can be used to support geotechnical site assessments. It has been used to guide land preparation and to integrate field data with existing data [17].
ANN is one of the prevalent algorithms among researchers nowadays, specifically in geotechnical issues. ANN holds three significant advantages: first, the counting speed is high. Second, it has a strong fault-tolerant capability. Third, it is proficient in dealing with problems with complicated problem-solving rules [18]. The technique of utilizing ANN could be a suggestion for predictions, especially in cases where theoretical modeling does not give foreseen outcomes [19]. ANN aims to model the behavior of the nervous system in the human brain. ANN is an adequate solution for solving complex and nonlinear data modeling. Ref. [20] presents the estimation of standard penetration test values on cohesive soil using an artificial neural network without data normalization. Some previous studies investigated the assessment of geotechnical properties and determination of shear strength parameters by unitized ANN [21,22,23]. In the geotechnical domain, the development situations generally have multiple variables, making them challenging to model employing conventional mathematics [24].
In this study, the test results of 102 boreholes were gathered, categorized, and analyzed and modeled using ArcGIS10.7 software. These data were used to construct models using ANNs to predict the SPT values and bearing capacity of soils. The data represent the area of Erbil City in Iraq, covering a depth of 9.5 m below the ground surface. Data were collected from the Andrea Engineering Test Laboratory and the construction laboratory in Erbil. Therefore, for a geotechnical engineer, this information can be used to classify areas into zones according to GIS results. The purpose of this study is to analyze and model the geotechnical parameters such as the fines content, moisture content, soil plasticity, shear strength parameters, compressibility parameters, SPT–N values and soil bearing capacity via GIS to create a group of maps in different layers. In addition, the prediction of SPT-N value and bearing capacity based on geotechnical parameters were modeled by ANN methods. Correlations between shear strength parameters, SPT values and bearing capacities of soils were made by Minitab 17 programming.
By combining the capabilities of GIS and ANNs, researchers can develop models to analyze geotechnical issues at different spatial scales, producing results that are more accurate and efficient compared to manual methods In summary, GIS and ANNs complement each other by providing an analytical function that is efficient for developing model data entry at different spatial scales for geotechnical issues, allowing for more accurate predictions and a better understanding of the relationships between soil properties and other factors. To the best of our knowledge, this is the first study to cover mapping and modelling all the geotechnical characteristics of soil in Erbil city.

2. Study Area

Erbil is located in the northwestern region of Iraq. Geologically, it is in the low–folded belt of northern Iraq in the structural trough with a NW–SE axial trend, and within the foothill zone, which is part of the stable shelf tectonic unit of Iraq. Erbil City has an area of approximately 250 km2 and GPS coordinates of 36°11′27.4″ N 44°00′33.7″ [25]. A location map of Erbil City is shown in Figure 1. From a geomorphological perspective, the area is flat with uncommon low–lying hills. In addition, Erbil City is stratigraphically covered by quaternary and Pleistocene deposits, which are dominated by clay, silt, and sand [26]. Erbil City is mainly covered by soils such as gravel and conglomerates with sand, clay, and silt. Conglomerates cover more than 80% of the study area [27].
In the last decade, Erbil city had extensive development in the construction of railways, internal and ring roads. Therefore, collecting data, building a database, producing GIS maps for soil properties, and developing a model for soil behaviors and the bearing capacity of the foundations would be very useful for site engineers to make immediate decisions regarding the selection of project positions. Finally, the geotechnical properties at different depths were collected from the study area (Erbil city), analyzed and tabulated for the 102 boreholes. The locations of the boreholes were selected to ensure a uniform distribution throughout the study area. The borehole locations are shown in Figure 2.

3. Methodology

The methodology of the study involved collecting data from the field in Erbil city center and analyzing it using both GIS and ANNs. The data was processed and analyzed using these tools to gain insights and conclusions about the study area. The flow chart represents in Figure 3 the methodology of this study as a tool to help readers understand the process used. A flow chart can show the different steps involved in integrating GIS, ANN, and lab analyses, making it easier for readers to follow and comprehend the study’s methodology.

3.1. Data Collection

In this study, the results of site investigations and numerous series of soil laboratory tests were collected from 102 borehole locations that cover the main part of the region. Soil investigation included drilling boreholes, taking both disturbed and undisturbed soil samples at 1–meter intervals from 1.0 m to 9.5 m depth, and various field tests. Laboratory work included a series of geotechnical tests to determine the soil index properties, sieve analysis, compressibility, settlement, and shear strength.
A sample of the data collected from the laboratory and field tests is presented in Table 1. This table includes some statistical information (e.g., min. values, max. values, average values, and standard deviation of the input and output data). All data of the study could find in Appendix A as a Table A1.

3.2. Geographical Information Systems (GISs)

The results of the soil investigation and field and laboratory tests were employed to create a digital database for the study region. A database of geotechnical properties was used to provide the input values for the mapping software. In this study, the data were analyzed and presented as maps using ArcGIS (10.7) software. Deterministic methods (inverse distance weighting) were used to create maps and models of spatial data, which rely on probability and uncertainty. Deterministic methods use a fixed set of rules or algorithms to create maps and models, unlike geostatistical methods used for analyzing and modeling spatial data [11]. These methods are commonly used to analyze patterns, relationships, and trends in large, complex datasets.
Some of the most commonly used deterministic methods in ArcGIS include:
  • Interpolation: This method is used to predict values at unsampled locations based on observed data. Interpolation methods in ArcGIS include inverse distance weighting, spline, and triangulated irregular network (TIN) interpolation.
  • Buffering: This method is used to create a polygon around a feature that represents a specified distance. Buffers are commonly used in spatial analysis to identify areas that are within a certain distance of a feature of interest.
  • Overlay: This method is used to combine two or more maps based on a set of rules or conditions. Overlays can be used to create a new map that shows the spatial relationships between features in the input maps.
  • Reclassification: This method is used to change the values of a raster or vector layer based on a set of rules or conditions. Reclassification is often used to simplify complex data or to create new data layers based on existing data.
  • Extraction: This method is used to select features from a map based on a set of conditions or rules. Extractions can be used to create new data layers that contain only the features that meet specific criteria [11].
The resulting digital maps illustrate the soil formation patterns, distribution, and geotechnical properties of soils at different depths. These maps simplify and help designers and site engineers make the right decisions in the construction of projects. The aim of drawing digital maps using the GIS method was to illustrate the bearing capacity of foundations at three different depths. In this study, the bearing capacity was estimated using two methods. The first Meyerhof method (1963) used shear strength parameters for (10 × 10 m). In the second method, standard penetration numbers were used for bearing capacity estimation [28].

3.3. Statistical Analysis

To make a correlation between the geotechnical properties, the data obtained from the soil investigation of all boreholes in the study area were correlated by MINITAB 17 software. For instance, correlations were made between the SPT values of soil strata with the shear strength parameters of soils and the ultimate bearing capacity. This method is used to find a correlation between the response (Y) and predictor (X) using regression analysis, which is an extensively used method for analyzing multifactor data.

3.4. Neural Network Model

The process for creating an artificial neural network is assumed by using the Matlab application. This study aims to make models by an artificial neural network. The network model developed was formed from data collection of geotechnical properties of soils in the study area. The ANN analysis result aims to predict SPT N-value and bearing capacity using the identical algorithm, the Back-propagation algorithm, and the same activation function. The network architecture was chosen using hidden layers and varying the number of neurons in the hidden layer. The relation number of neurons in the hidden layer is between 15 and 18 according to previous researches [29,30]. The network performance that has the smallest error and the correlation coefficient value that is proximate to 1 is most suitable for data predictions. Root Mean Squared Error (RMSE) is a commonly used evaluation metric in ANN models. RMSE measures the difference between the predicted and actual values, and it is expressed in the same units as the target variable. A low RMSE value indicates that the predictions are close to the actual values, while a high RMSE value displays that the predictions are far from the actual values. In ANN models, RMSE is used to evaluate the performance of the model and determine the quality of the predictions. A lower RMSE value shows a better fit between the predictions and the actual values and a more accurate model. In this research, the R2, RMSE, and MAE values of the estimated and actual target parameters are computed in the implementation evaluation of regression models. The R2, RMSE, and MAE represented mathematically as Equations (1)–(3):
R 2 = 1 I = 1 N ( y mea y pre ) 2 I = 1 N ( y mea y m ) 2
RMSE = I = 1 N ( y mea y pre ) 2 N
MAE = 1 N I = 1 N | y mea y pre |
where ymea, ypre, and ym represent the average of existing output, predicted output, and actual output, respectively. N denotes for all number of data. The degree of fitting is raised with R2 immediacy to 1. RMSE and MAE are utilized to assess the model’s prediction capability. For the RMSE and MAE, the prediction model will be more exact and its accuracy will be higher with a smallish value.
Figure 4 illustrates the structure of neural network models to predict SPT N value as output with two models: (a) using input as (LL%, PL%, PI%, WC, cohesion, ϕ, Fine content) and (b) using inputs as (LL%, PL%, PI%, WC, ϕ, Fine content). The structures of neural network models to predict ultimate bearing capacity are presented in Figure 5. The parameters were used as input in two models: (c) (LL%, PL%, PI%, WC, cohesion, ϕ, Fine content) and (d) (LL%, PL%, PI%, WC, ϕ, Fine content).

4. Results and Discussions

4.1. Modeling of Soil Properties Using GIS Maps

GIS maps of the soil properties of the study area were produced. These maps included the distribution of the soil fines content, natural water content, liquid limit, plastic limit, cohesion, angle of internal friction, compression index, rebound index, SPT–N values, and bearing capacity of foundations at three different levels (1.5 m–3.5 m, 3.5 m–6.5 m, and 6.5 m–9.5 m). The soil characteristics at different depths were interpolated for the survey area to show the distribution of these properties in a clear way. In general, the soil characteristics in all the maps were divided into six major legends, each of which was represented by a unique color.

4.1.1. Fines Content Model

Fines content in soils is one of the most significant parameters that affect soil behaviors such as shear strength, compressibility, plasticity, and indirectly, the bearing capacity of foundations [31]. The results presented in Figure 6 indicate that the majority of the study area in Erbil city center has a high fines content. According to the figure, the fines content in most of the area is greater than 81%, meaning that the proportion of soil particles that are smaller than 0.075 mm in diameter is high. This high fine content will likely have significant effects on the soil’s engineering properties, such as its shear strength, compressibility, and plasticity. The impact of soil properties on its ability to support loads, resist deformation, and transmit loads must be carefully considered in future construction and development projects in the area. The results from Figure 6 provide crucial information for engineers and planners in Erbil city center, emphasizing the need to consider the fines content of soil in decision-making for development projects. The compression level of soil becomes crucial when large particles are replaced by fine particles, and the impact of fines content is more pronounced when the soil is near saturation. This has been noted by other researchers in the field [31,32,33]. The relationship between natural water content and fines content is intertwined and cannot be evaluated separately. Sometimes, improving natural water content weakens the effect of fines content on soil shear strength, due to the sensitivity of fines to changes in natural water content. In dry conditions, fines do not significantly affect soil behavior due to the influence of suction [34,35].

4.1.2. Atterberg Limits

Atterberg limits can be used to characterize soil behavior and classification, including the swelling potential of expansive soils, the consistency and plasticity of the soil. Two essential index properties can be obtained from the values of the natural water content: liquid limit and plastic limit [36]. Soil water content affects its consistency. A high water content in clay makes it behave like a liquid due to the reduction of attraction between clay particles caused by the excess water between the particles [37]. The liquid limit was used to determine the consistency of the fine-grained soil. This measure of soil consistency is useful for estimating soil consolidation properties and calculating the acceptable bearing capacity and settlement of foundations [38].
The variations in the soil liquid limit throughout the study area at the three different depths are illustrated in Figure 7. The results of the analysis of liquid limit values in the study area show that a significant portion of the soils have a liquid limit range between 40% and 52%. This range indicates the presence of low plastic clay, which is a type of soil that has low resistance to deformation when subjected to stress and is prone to collapsing. The presence of low plastic clay in high percentages in the study area is an important factor to consider in construction and development projects, as it may impact the stability and integrity of structures built on these soils. Additionally, the observation of soils with high liquid limit values (greater than 53%) in relatively central regions is also important. These central zones are considered critical points in the study area, as soils with high liquid limit values are prone to deformation and instability. These critical points need to be carefully evaluated and addressed in any future construction and development projects in the area to ensure the stability and integrity of the structures built on these soils. Numerous researchers have investigated the relationship between the liquid limit of soils and swelling potential. Some types of clay minerals with a high cation exchange capacity (CEC) suffer from expansion and an increase in the volume of the available water [38,39,40,41].
Plasticity is one of the most important features of clay, and the crystallinity of clay minerals is the primary source of this plasticity [39,40,41]. Soil is plastic when the water content is below the liquid and plastic limits. The plastic range, which is the difference between the two limit values, is called the plasticity index [42]. The plastic limit provides geotechnical engineers with indirect information about the activity, toughness index, and optimum moisture content of soils. Figure 8 shows the plastic limit variation in the study area at different depths. The analysis of the figure reveals that half of the study area has a plastic limit of soils that ranges between 19% and 24%, while the other half has a plastic limit ranging between 25% and 30%. This indicates that the soils in the study area have different levels of plasticity, which is a measure of the soil’s ability to change shape without breaking under stress. The presence of soils with low plastic limits (Figure 8a) in small zones in the southeast direction is also an important observation. Soils with low plastic limits are less plastic and more brittle, making them more prone to cracking and failure under stress. These areas need to be carefully assessed and addressed in future construction projects in the area to ensure the stability and reliability of the structures built on these soils.

4.1.3. Natural Water Content Model

The natural water content can be considered a parameter that profoundly affects the geotechnical properties [43,44]. Figure 9 shows the variations in the natural water content in Erbil City. The figure presents the water content at depths of 1.5 m–3.5 m, 3.5 m–6.5 m, and 6.5 m–9.5 m. The results show that the natural water content in the study area is mostly found to range between 16% and 20% at the three levels of investigation depths. The natural water content of soils is an important factor to consider in construction and development projects, as it affects the soil’s stability, bearing capacity, and compressibility. Soils with high natural water content are more susceptible to instability, while soils with low natural water content are more prone to drying and cracking.
This is consistent with the findings of many researchers working in this region [43,44]. The water table, during the time of exploration, was in very high depths below the natural ground level (NGL), meaning the groundwater was relatively shallow. The water table fluctuates seasonally, with an increase during spring. The soil above the water table affects its strength and compressibility, as more moisture results in decreased strength and increased compressibility. Saturated soil below the water table creates settling issues as the consolidation process reduces the natural water content under stable load.

4.1.4. Shear Strength Parameters

The shear strength of soils is a crucial parameter in many foundation engineering designs. It refers to the ability of soil to resist forces that cause slipping or sliding along a plane within the soil. The shear strength of soils is important in determining the bearing capacity of shallow and deep foundations, slope stability, tunnels, and lateral pressure on structures [28]. The soil’s shear strength comes from its cohesive strength (c) and frictional strength, represented by the angle of internal friction (ϕ). Cohesive strength is due to the bonding force between soil grains and the binding material, while frictional strength arises from the friction, interlocking, and rolling of soil grains [44,45,46,47]. The strength of any soil decreases as the shear strain and expansion or contraction increase or decrease, respectively, with respect to the soil density due to applied loads. Shear strength parameters are widely utilized by different standard equations in the design of foundations, particularly in empirical equations. [48].
Figure 10 and Figure 11 show the variations in the cohesion and angle of internal friction, respectively. The combination of these two parameters produces the shear strength of the soil, its variation throughout the study area, and the soil depth. Most of the soils in the study area at shallow depths (1.5 m–3.5 m) had cohesion values between 76 kPa and 100 kPa. However, there were relatively small regions with cohesion values greater than 100 kPa and less than 50 kPa. The area covered with soil had cohesion values greater than 100 kPa, which increased with the depth. The results of the figures are in agreement with the distribution of the fines content in the study area. The fines content of soils in the majority of the study area was found to be greater than 80%, and this high fine content is likely contributing to the high cohesion values in the area.
Soils with high fines content produce higher cohesion values [49]. The angle of internal friction is a parameter of the soil shear strength and is employed in bearing capacity estimation, slope stability analyses, and estimation of soil lateral earth pressures. The soils in the study area were found to have an angle of internal friction between 2° and 6°, which was found to be similar at all three depths in this investigation. The east-south part of the study area had soils with an angle of internal friction ranging from 7° to 12°. These results indicate that the angle of friction values found in the study area are consistent with those found in similar studies in the region. The angle of internal friction is an important parameter in determining the shear strength of the soil, and these results suggest that the soils in the study area have moderate to low shear strength values. The results of the figures can be used to identify zones of high and low shear strength, which are important in determining the suitability of the soil for different types of structures. The information can also be used to design and construct structures that are appropriate for the soil conditions in the area [47,50].

4.1.5. Consolidation Parameter Model

The consolidation parameters (compression index and swelling index) of saturated clayey soil should be checked during the analysis and design of the foundations [51]. In this study, the distribution of the compression index and swelling index throughout the study area and their variation with the soil depth were investigated. Figure 12 shows that the compression index of soils in the study area, Erbil city center, has a range of values between 0.17 and 0.22, with a lower value found at greater depths. The compression index is an important factor in determining soil compression and consolidation. The results suggest that the shallow soil strata in the study area have high consolidation and settlement potential, while the settlement potential is expected to decrease with depth due to the variation in compression indices.
The rebound index is an important parameter in geotechnical engineering as it represents the soil’s tendency to expand or contract under changes in moisture content. A high rebound index value designates that the soil is more susceptible to expansion and contraction, whereas a low value indicates that the soil is more stable. The results shown in Figure 13 suggest that the soil in the study area has moderate to low rebound index values in the range 0.015–0.078, indicating that the soil is relatively stable and less likely to expand or contract under changes in moisture content. This indicated that there was no soil swelling or shrinkage potential in the study area. The minimum values showed a consistent trend across all levels in the study area (Figure 13a), with the highest swelling index values located north of Erbil at the first level and in the south at the second and third levels. Most of the study area had moderate parameter values of 0.015–0.078. The consolidation parameters were consistent with previous studies [47,50,52].

4.1.6. Standard Penetration Test Model

The SPT is one of the many standard in situ tests used to identify soil type, stratigraphy, and relative strength measures during site investigations [53,54]. The eastern part of Erbil had higher SPT values, which is attributed to the higher unit weight of soils and the presence of stiffer and stronger soil layers as shown in Figure 14. The SPT values of soils at a depth of 6.5 m–9.5 m was mostly between 35 and 60, and the range increased from the west to the east of the study area. The increase in the SPT values with increasing depth is due to the influence of several factors. Soil type affects the SPT values as different soil types have different characteristics such as density, porosity, and strength, which all impact the SPT results. Unit weight, or the weight per unit volume of soil, is also a factor, as a higher unit weight results in higher SPT values. Confining pressure, the pressure applied to the soil from the surrounding material, also increases with depth, leading to higher SPT values. Overall, the results of SPT tests provide a good indication of the soil strength and its variation throughout the study area. To confirm the relationship between the SPT values and shear strength parameters (cohesion and angle of internal friction), data obtained from the boreholes were correlated. Figure 15 and Figure 16 illustrate the correlation between the SPT values and shear strength parameters for the soils within the study area. The figures evidently show that there is a good correlation between the SPT values, cohesion, and angle of internal friction. Generally, most values of SPT–N distribution from the center to the west of Erbil City were found to be between 17 and 48, while in the northeast and southeast, the values were higher than 50. As mentioned in various studies [55,56,57], the results of SPT–N showed that the soil in Erbil city is medium to very dense.

4.1.7. Bearing Capacity

The design of a footing depends on its soil’s bearing capacity. Many methods for estimating soil bearing capacity exist, relying on factors such as soil shear strength, footing type, and SPT value [58]. In this study, the soil strata’s ultimate bearing capacity was estimated using two methods, one based on shear strength parameters using Meyerhof’s equation, and the other based on SPT-N values. The variation in the soil strata’s ultimate bearing capacity based on Meyerhof’s equation is shown in Figure 17.
Figure 18 shows the ultimate bearing capacity of the soil strata in the study area at three different depths that were estimated from the SPT values. Changes in the ultimate bearing capacity were observed throughout the study area. The ultimate bearing capacity of the soil strata in the study area is a significant factor in the design of shallow and deep foundations. The results from Figure 18 indicate that the majority of the study area at shallow depths had an ultimate bearing capacity be-tween 170 kPa and 940 kPa. However, there were some small areas with lower ultimate bearing capacities. The ultimate bearing capacity increased with increasing depth, which can be attributed to the increase in soil confining pressure and soil unit weight. These findings are useful in determining the design and load-bearing capacity of foundations in the study area.
Figure 19 presents the correlation between the ultimate bearing capacity from Meyerhof’s equation and the SPT values of the soils within the study area. This figure indicates a good correlation between the SPT values and the ultimate bearing capacity of the foundations. The correlation between the ultimate bearing capacity from the SPT values of the soil strata within the study area and soil cohesion (shear strength parameter). The ultimate bearing capacity from Meyerhof’s equation and the ultimate bearing capacity from the SPT values of the soil strata within the study area were correlated, and the results are presented in Figure 20. The estimation of ultimate bearing capacity from the Standard Penetration Test (SPT) values is a relatively simple method compared to other methods that require more experimental tests and complex equations established on the shear strength parameters of the soil. The estimation of ultimate bearing capacity from SPT values is based on empirical correlations and has been widely used in the field of geotechnical engineering. The advantage of this method is its simplicity; however, it may not accurately reflect the actual soil conditions, which can be affected by factors such as soil type, stratigraphy, and loading conditions [59].

4.2. Artificial Neural Network Models

ANN models are commonly used for regression and classification tasks, including prediction problems. ANN models consist of interconnected nodes, or artificial neurons, that process and transmit information. The nodes are organized into layers, and the connections between them can be adjusted during the training process to minimize the error between the predicted and actual outcomes. ANN models have been widely used in various fields due to their ability to handle complex relationships and make accurate predictions. The relationship between the predicted original SPT values and bearing capacity in the training data and test data was demonstrated in four models that have an acceptable R2 and the smallest error. After input and output data are gathered and structured, training and test sets were established. 70% of the data were used for training and 15% used for testing and 15% for validation of the total data of boreholes. Predictions for SPT-N value from two models were developed.

4.2.1. Validation of Interpolations Based on Semivariograms

The transformed data’s spatial autocorrelation is modeled using semivariograms/covariance modeling for SPT-N values at depth 1.5–3.5 m, depicting the similarity decrease between data points as their distance increases. Binned values (red dots) are generated by grouping semivariograms/covariance points using square cells, while average points (blue crosses) are generated by binning empirical points in angular sectors. Binned points show local variation, while average values show smooth variation. A stable type model (dark blue line) is fitted to the empirical variogram for the measured data points.
At h = 0, the semivariogram should be 0. However, at an infinitesimally small separation distance, the semivariogram often exhibits a nugget effect, which is some value greater than zero. In this case, the nugget effect exists which is zero. The range is the distance at which the model levels out (5750.634). Locations closer than the range are spatially auto correlated, while farther locations are not. The partial sill is the sill minus the nugget (462.1716). The lag size is the distance class size (674.7548) with 12 lags. Semivariogram values are shown in Figure 21 with higher values in orange/red and lower values in blue/green.
Kriging estimates unknown spatial values. The search neighborhood step involves selecting nearby points with significant influence on the prediction location, determined by spatial auto-correlation. The method eliminates irrelevant points and weights nearby points using a search neighborhood of adjacent points, radius, and number of sectors to estimate values at the unknown location. Accurate neighborhood identification and selection of nearby points are crucial for successful kriging. As shown in Figure 22, five neighboring points are selected and a circle with four sectors is selected. The points highlighted in the data give an indicator of the weights associated with each point, and these weights are used to estimate the value at the unknown location, which is at the center of the crosshair.

4.2.2. Prediction for SPT-N Value

In Figure 23 and Figure 24 the results of the ANN model (a) show a good agreement between the predicted and measured SPT-N values, with an R2 of 0.92 for the training data and 0.81 for the testing data. The model uses Atterberg limit values, water content, cohesion, and internal friction as input variables and predicts the SPT value as its output. On the other hand, to predict SPT N value as output with (LL%, PL%, PI%, WC, ϕ, Fine content) as input, the result showed R2 values of 0.90 and 0.8 for training and testing, respectively. As mentioned in previous studies [60], the predictions of the SPT values in model (a) and model (b) were conducted with more superficial R2 values that give a significant agreement for using ANN modeling in geotechnical engineering that helps the engineering to be utilized in the design of infrastructures.

4.2.3. Prediction of Ultimate Bearing Capacity

The prediction of the ultimate bearing capacity of soil using ANNs is a machine learning approach that involves training a neural network model with a dataset of soil properties and corresponding ultimate bearing capacity values. The trained model can then be used to predict the ultimate bearing capacity of new soil samples based on their properties. Advantages of ANN prediction include the handling of non-linear relationships between soil properties and ultimate bearing capacity, the ability to incorporate complex relationships between soil properties and ultimate bearing capacity, and the ability to handle large datasets with many input variables.
To determine the best prediction results for the ultimate bearing capacity, two models were used. The inputs and outputs of the models are listed in Table 2 along with the value of R2, which is a measure of the goodness of fit of the model to the data. R2 ranges from 0 to 1, with a higher value indicating a better fit between the model and the data. The R2 value provides an indication of how well the models are able to predict the ultimate bearing capacity based on the inputs.
The results of bearing capacity prediction using ANN modeling show similarities with the results of previous studies [60]. In Figure 25 and Figure 26 the relationship between the predicted ultimate bearing capacity and the original ultimate bearing capacity on the training and testing data are illustrated.
Model (c) has a high R2 value, indicating a strong correlation between the predicted and actual values, and a low variance in the residuals. The high R2 value for the training data (0.91) indicates that the model is able to fit the training data well, while the R2 value of 0.82 for the testing data suggests that the model has good generalization ability and can predict unseen data with a certain degree of accuracy. In Model (d), the decrease in R2 value is likely due to the absence of the cohesion value as an input. Cohesion is an important factor that affects the bearing capacity of soil, so its absence in the input could result in a decrease in the accuracy of the model. This highlights the importance of considering all relevant factors in the input variables of the model to improve its accuracy and predictability.
Soil cohesion is a measure of the shear strength of soil, which determines its ability to support loads. The results of research in geotechnical engineering have shown that soil cohesion is the most important factor in estimating soil bearing capacity. This is because it determines the resistance of soil to sliding or deformation under load, and is essential for ensuring the stability of structures built on the soil. Therefore, accurate determination of soil cohesion is crucial for safe and effective design of geotechnical structures.

4.2.4. Percentage Error of ANN Models

Figure 27 and Figure 28 demonstrate the error percentage lines for model (a) and (b) outcomes referring to the difference between the predicted values and the actual values. It is a measure of the accuracy of the model’s predictions. A lower error percentage indicates a more accurate model, while a higher error percentage means that the model is less accurate. The acceptable error percentage depends on the specific application and the acceptable level of error for that particular field. In some cases, a low error percentage, such as 1–2%, may be acceptable, while in others, a higher error percentage may be acceptable if the cost of a lower error is too high [61]. It is important to note that no model can be 100% accurate and some level of error is always present. The goal is to reduce the error to the lowest possible level while still making predictions that are useful and relevant.

4.2.5. Analysis of Models

Sensitivity analysis is an important aspect of model evaluation, as it helps to assess the impact of individual model parameters on the final results. The absence of a sensitivity analysis in the study you are referring to could limit the understanding of the model’s behavior and the confidence in the results obtained. Analysis of models is particularly important when using ANOVA, as the parameters of these models can have a significant impact on the spatial structure of the data. By conducting a sensitivity analysis, the researchers could determine the robustness of the model results. The sensitivity analysis to evaluate the effect of each model parameter based on the semivariograms shown in Table 3.
A p-value less than 0.05 is commonly used as a threshold for determining statistical significance in hypothesis testing. In this context, a p-value less than 0.05 for a particular model parameter means that the results suggest a statistically significant relationship between the parameter and the outcome variable being modeled.
In this study, a p-value less than 0.05 for a particular model parameter would indicate that the parameter has a significant impact on the model results. This information can be useful in understanding the underlying relationships in the data and in guiding further analysis or interpretation of the results. However, it is important to consider other factors such as the sample size, the quality of the data, and the overall fit of the model in evaluating the reliability and robustness of the results.

5. Conclusions

This study focuses on developing maps for soil geotechnical properties that are widely utilized by geotechnical engineers in foundation design capacity. The main conclusions are as follows. An artificial neural network (NN) model was established for estimating SPT-N value and ultimate bearing capacity:
  • GIS is an effective tool that can be used by engineers to analyze the preliminary exploration of geotechnical sites. Information from 102 boreholes, considering the main geotechnical properties, was collected, evaluated, and used as input data for GIS analysis.
  • This information suggests that a significant portion of Erbil city has soil with a high proportion of fine-grained materials, such as clay and silt. High fines content can impact the soil’s physical and engineering properties, such as its compressibility, permeability, and shear strength. The presence of high fines content can also increase the susceptibility of soil to swelling and shrinkage, which can lead to instability in structures built on or in the soil. The small zones in the southeast of the study area with lower fine contents may have different soil characteristics and may offer potential sites for structures that require more stable soil conditions. These findings are important for the design of infrastructure and buildings in the city.
  • Atterberg limits in most of Erbil City were found to be between 40% and 52%, and 19% and 30% for the liquid and plastic limits, respectively. This indicates the high presence of low–plasticity clay and clayey silt. The results of the analysis of liquid limit and plastic limit values in the study area provide important information for engineers and planners in Erbil city center. They highlight the presence of low plastic clay in high percentages in the study area, as well as the need to carefully evaluate critical points with high liquid limit and plastic limit in future construction and development projects.
  • Digital mapping of shear strength parameters showed that most soil strata at three different depths had an internal friction angle between 2° and 6°, and the cohesive strength ranged between 76 kPa and 130 kPa. The results of the cohesion values show that the soils in the study area at shallow depths have moderate to high cohesion values, and that the soils with high cohesion values tend to be located in areas with high fines content. However, the results of the angle of internal friction show that the soils in the study area have moderate to low shear strength values, with the soils in the east-south part of the area having slightly higher shear strength values. These findings are important in determining the suitability of the soil for different types of structures and in designing and constructing structures that are appropriate for the soil conditions in the area.
  • The soil in the study area mostly has a moderate compressibility and resilience, with a moderate to low amount of rebound. The compression index decreases with depth, suggesting that the soil becomes less compressible as one moves deeper into the ground. The rebound index indicates that the soil has a moderate to low ability to recover its original volume after being compressed. These findings provide valuable information for designing structures that are built on or into the soil in the study area.
  • SPT values in the study area indicate moderate soil strength in the shallow strata, with a range of 17 to 48. As the depth of the soil strata increases, the SPT values increase and become higher, covering large parts of the study area. This suggests that the soil becomes stronger with increasing depth. The SPT is a widely used in-situ test for measuring soil strength, and these results provide valuable information for the design of foundations and other structures that are supported by the soil. The higher SPT values at greater depths indicate improved soil strength characteristics and can influence the design of these structures in terms of load-bearing capacity and stability.
  • This conclusion suggests that the soil in Erbil City is not capable of supporting heavy loads without modification or special design measures. The ultimate bearing capacity is a measure of the maximum weight or load that a soil can support without failure. A value lower than 170 kPa indicates that the soil may not be suitable for supporting heavy structures, such as buildings and bridges, without additional treatment or specialized foundation design. Improving the soil, such as through compaction or stabilization, and utilizing special footing designs, such piles, can increase the soil’s bearing capacity and ensure the stability and safety of structures built on the soil.
  • At the preliminary design point, the completed digital geotechnical maps are vital. The designer could use the geotechnical parameters, consolidation characteristics and SPT as an effective visual display tool simply by using the digital values of these parameters for the proposed region, where the necessary decisions can be made.
  • The correlation between the SPT values and shear strength parameters for the soils within the study area demonstrated a strong relationship between them.
  • The results obtained from the models were compared with those measured from the field tests. It was found that predicted SPT-N values and Q-ultimate bearing capacity are quite close to the measured values. In order to check the prediction performance of the ANN model developed, several performance indices, such as R2, MAPE, and RMSE were also calculated. The ANN model has shown good prediction performance based on the performance indices. Thus, the developed ANN model can be used to predict SPT-N and Q-ultimate bearing capacity from the soil parameters and borehole coordinates. The ANN model’s implementation has also demonstrated that the neural network is a valuable tool to minimize the uncertainties encountered during geotechnical engineering projects. Therefore, using Artificial Neural Networks may provide new techniques and methodologies and minimize the potential inconsistency of correlations. ANN prediction is a useful tool for predicting the ultimate bearing capacity of soil, but it should be used in conjunction with other methods and validated with independent data to ensure accurate predictions.

Author Contributions

Conceptualization, K.Z., Z.B.Q. and A.F.C.; methodology, Z.B.Q.; software, Z.B.Q.; validation, Z.B.Q., A.F.C. and Z.K.; formal analysis, Z.B.Q.; investigation, Z.B.Q.; resources, Z.B.Q.; data curation, Z.B.Q.; writing—original draft preparation, Z.B.Q.; writing—review and editing, A.F.C., Z.B.Q. and Z.K.; visualization, Z.B.Q.; supervision, A.F.C. and Z.K.; project administration, A.F.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available from the conforming author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The data of all boreholes are provided demonstrated in Table A1.
Table A1. Full data of all boreholes of study area.
Table A1. Full data of all boreholes of study area.
BH
NO
DEPTH
m
XYLL%PL%PI%WC%c (kN/m2)Φ
(°)
Fine ContentSPT-N Value
(kN/m2)
11.5–3.5406,851.54,009,57046222428.351462.17
21.5–3.5405,972.94,009,284.24022181596394.154
31.5–3.5407,487.64,007,947.14725221899591.557
41.5–3.5406,121.14,008,204.740221815.096394.154
51.5–3.5407,443.24,008,702.750272320.8109592.941
61.5–3.5407,941.44,009,347.750282213.9108491.845
71.5–3.5408,217.94,008,512.2335201920.09649641
81.5–3.5408,997.44,008,917.839211816.3109589.445
91.5–3.5408,743.44,009,398.147252218.099591.557
101.5–3.5408,026.14,010,40645242116.696595.563
111.5–3.5409,617.84,011,321.646271920.299590.741
121.5–3.5408,588.14,011,780.245252019.097491.660
131.5–3.5409,762.14,012,805.448262219.794581.375
141.5–3.5411,184.24,011,787.249232619.6108496.333
151.5–3.5411,184.34,012,540.747252215.6105597.430
161.5–3.5411,247.74,013,522.948262218.458392.719
171.5–3.5409,969.94,010,643.858312720.591488.940
181.5–3.5409,662.94,009,97756302618.894469.142
191.5–3.5409,915.14,009,232.245242113.677395.929
201.5–3.5410,862.14,009,607.3544251918.369594.530
211.5–3.5411,994.64,009,216.942231919.270498.229
221.5–3.5411,166.14,010,438.950262417.995592.258
231.5–3.5412,126.94,010,308.351292919.895596.270
241.5–3.5411,550.54,011,441.244232114.895710080
251.5–3.5413,143.34,011,488.944251920.597493.541
261.5–3.5414,064.74,011,692.841231817.4110453.637
271.5–3.5415,317.34,010,885.450262415.9116568.523
281.5–3.5414,149.34,009,671.451272419.194487.319
291.5–3.5414,362.54,008,450.652242819.790486.417
301.5–3.5416,171.14,012,249.645242115.398679.843
311.5–3.5415,609.84,012,761.838211716.095570.645
321.5–3.5414,549.14,013,637382117 701297100
331.5–3.5417,922.54,010,46253272621.669392.625
341.5–3.5416,961.14,011,319.349232615.1112586.771
351.5–3.5416,542.14,010,233.543232017.891389.379
361.5–3.5416,909.44,008,75444242014.892466.285
371.5–3.5406,775.34,007,107.550262418.382492.118
381.5–3.5405,227.44,007,274.353272617.374396.425
391.5–3.5407,403.94,006,885.340221813.6105563.819
401.5–3.5407,715.14,007,488.543222119.190493.220
411.5–3.5408,445.44,007,069.445242115.858595.919
421.5–3.5408,736.74,005,101.944251917.660495.522
431.5–3.5408,248.54,006,694.743241919.969492.122
441.5–3.5407,872.44,006,657.2256253115.6116595.528
451.5–3.5409,311.94,005,78646242217.7112497.432
461.5–3.5408,775.594,006,631.2446222421.185495.819
471.5–3.5409,973.44,006,13353252820.987593.631
481.5–3.5409,048.644,007,069.3946242216.689496.229
491.5–3.5409,874.194,007,406.146242217.792594.137
501.5–3.5408,985.144,007,907.654262815.399494.139
511.5–3.5411,015.24,007,902.755243115.3112595.282
521.5–3.5418,521.84,007,163.457253219.6104492.977
531.5–3.5409,785.254,008,333.148222614.497495.522
541.5–3.5410,998.64,008,878.349262318.8105492.232
551.5–3.5410,109.14,006,834.4441192216.8101487.341
561.5–3.5410,419.84,006,00145212416.3123552.328
571.5–3.5412,049.024,007,215.9752242821.8109551.844
581.5–3.5411,2634,006,447.551272419.496494.233
591.5–3.5412,664.184,007,328.449252419.4124492.638
601.5–3.5411,9904,006,315.244251921.7136595.739
611.5–3.5413,164.94,008,77943232018.9119596.641
621.5–3.5412,846.744,006,559.8149262317.479493.539
631.5–3.5412,536.54,005,174.248272118.081493.535
641.5–3.5415,747.44,005,960.6948232516.896496.421
651.5–3.5414,814.84,007,156.940211917.997539.731
661.5–3.5415,071.94,004,866.947242319.095590.834
671.5–3.5417,971.234,001,685.0148272117.48883061
681.5–3.5418,5194,003,222.900017.785101837
691.5–3.5415,9264,002,87800019.08391450
701.5–3.5416,1674,004,572.300011.08081536
711.5–3.5419,2894,001,977.600018.611452540
721.5–3.5411,924.274,004,297.5957253216.4116594.772
731.5–3.5410,482.294,003,953.6348232513.4109591.474
741.5–3.5412,426.984,002,802.6945212412.1144495.623
751.5–3.5414,3654,003,40844212313.7135494.620
761.5–3.5411,975.734,001,202.9547212613.8118593.239
771.5–3.5414,022.94,002,38745222315.495494.744
781.5–3.5411,697.914,000,107.5741271411.955677.88
791.5–3.5410,2003,999,70635231212.647383.711
801.5–3.5413,372.734,001,528.3800013.79858.4100
811.5–3.5414,253.794,000,845.7647222517.586591.781
821.5–3.5415,598.54,000,25945222420.189393.583
831.5–3.5413,510.133,999,595.5747252220.581395.627
841.5–3.5412,160.933,999,051.8841221912.7123610045
851.5–3.5410,732.93,999,78540211910.815303347
861.5–3.5410,2533,998,4454323205.719354253
871.5–3.5413,3093,998,392.541212012.3108581.660
881.5–3.5410,136.164,001,490.7255262917.381396.222
891.5–3.5409,296.114,002,026.5146242222.196596.831
901.5–3.5408,932.64,003,851.950242614.890397.329
911.5–3.5408,515.584,002,310.9344301431.125290.118
921.5–3.5409,9394,003,009.644281629.229389.212
931.5–3.5407,821.054,002,866.5642212123.245491.619
941.5–3.5407,298.54,004,255.6249262316.738389.85
951.5–3.5405,9374,003,22443232013.390493.37
961.5–3.5407,774.754,005,380.138211716.998496.639
971.5–3.5407,844.334,005,870.1139221718.943456.719
981.5–3.5406,914.854,005,671.1545232214.5121694.623
991.5–3.5405,2424,004,431.746252115.6110498.125
1001.5–3.5405,684.544,005,717.4543212214.895494.341
1011.5–3.5408,025.123,999,907.546242215.7103496.254
1021.5–3.5406,400.44,001,12442261615.594496.175
1033.5–6.5406,851.54,009,57049232629.440483.511
1043.5–6.5405,972.94,009,284.248232529.1455879
1053.5–6.5407,487.64,007,947.150242619.2109595.172
1063.5–6.5406,121.14,008,204.746262019.485497.960
1073.5–6.5407,443.24,008,702.752262623.2112494.970
1083.5–6.5407,941.44,009,347.752292314.8103493.777
1093.5–6.5408,217.94,008,512.2340221820.5113510090
1103.5–6.5408,997.44,008,917.8362313 105598.473
1113.5–6.5408,743.44,009,398.144242019.3102595.186
1123.5–6.5408,026.14,010,40642231915.5102694.180
1133.5–6.5409,617.84,011,321.664333119.589480.597
1143.5–6.5408,588.14,011,780.248252320.193589.199
1153.5–6.5409,762.14,012,805.446252117.298572.5100
1163.5–6.5411,184.24,011,787.245242114.3104495.750
1173.5–6.5411,184.34,012,540.749282118.2132692.962
1183.5–6.5411,247.74,013,522.945271819.062410025
1193.5–6.5409,969.94,010,643.851272418.793594.391
1203.5–6.5409,662.94,009,97754292519.395594.289
1213.5–6.5409,915.14,009,232.247232417.587594.633
1223.5–6.5410,862.14,009,607.3546222420.782498.241
1233.5–6.5411,994.64,009,216.945212317.072597.535
1243.5–6.5411,166.14,010,438.954252918.691497.876
1253.5–6.5412,126.94,010,308.355262919.393492.179
1263.5–6.5411,550.54,011,441.246222414.8978100100
1273.5–6.5413,143.34,011,488.946252121.199495.275
1283.5–6.5414,064.74,011,692.860322821.2147586.650
1293.5–6.5415,317.34,010,885.450282220.298591.956
1303.5–6.5414,149.34,009,671.442231919.495541.6100
1313.5–6.5414,362.54,008,450.671363519.795559.743
1323.5–6.5416,171.14,012,249.638211716.895582.682
1333.5–6.5415,609.84,012,761.844251916.197685.883
1343.5–6.5414,549.14,013,63742231910.47511100100
1353.5–6.5417,922.54,010,46248262220.276596.137
1363.5–6.5416,961.14,011,319.3482325 114569.175
1373.5–6.5416,542.14,010,233.545261921.694594.987
1383.5–6.5416,909.44,008,75454322219.694489.483
1393.5–6.5406,775.34,007,107.553282519.378397.945
1403.5–6.5405,227.44,007,274.354282616.882498.632
1413.5–6.5407,403.94,006,885.338211717.398592.539
1423.5–6.5407,715.14,007,488.547232421.592596.155
1433.5–6.5408,445.44,007,069.444232119.381396.340
1443.5–6.5408,736.74,005,101.946232319.287493.838
1453.5–6.5408,248.54,006,694.752262622.393594.653
1463.5–6.5407,872.44,006,657.2257263118.7106595.661
1473.5–6.5409,311.94,005,78655262917.6103496.963
1483.5–6.5408,775.594,006,631.2447232515.896592.270
1493.5–6.5409,973.44,006,13343212220.5107487.969
1503.5–6.5409,048.644,007,069.3946232317.390494.650
1513.5–6.5409,874.194,007,406.145222317.59259473
1523.5–6.5408,985.144,007,907.656302618.2110593.372
1533.5–6.5411,015.24,007,902.763283521.810269491
1543.5–6.5418,521.84,007,163.456243216.6129596.186
1553.5–6.5409,785.254,008,333.151272418.2118696.440
1563.5–6.5410,998.64,008,878.349262318.799596.3100
1573.5–6.5410,109.14,006,834.4446192717.298590.330
1583.5–6.5410,419.84,006,00148202818.2126587.250
1593.5–6.5412,049.024,007,215.9754282617.8134696.463
1603.5–6.5411,2634,006,447.555243119.599594.279
1613.5–6.5412,664.184,007,328.446262017.2118496.622
1623.5–6.5411,9904,006,315.246252117.1121596.538
1633.5–6.5413,164.94,008,77947232416.3115669.136
1643.5–6.5412,846.744,006,559.8145242123.876494.745
1653.5–6.5412,536.54,005,174.247262118.989598.141
1663.5–6.5415,747.44,005,960.6945212420.4126598.4100
1673.5–6.5414,814.84,007,156.944222220.4120596.4100
1683.5–6.5415,071.94,004,866.946232318.9117590.8100
1693.5–6.5417,971.234,001,685.0145252017.49193976
1703.5–6.5418,5194,003,222.945242117.795116490
1713.5–6.5415,9264,002,87841231919.089113968
1723.5–6.5416,1674,004,572.340202418.07695599
1733.5–6.5419,2894,001,977.600018.6112578.892
1743.5–6.5411,924.274,004,297.5962263616.6125568.380
1753.5–6.5410,482.294,003,953.6349232616.6115617.669
1763.5–6.5412,426.984,002,802.6948252318.1140539.9100
1773.5–6.5414,3654,003,40839211817.592498.180
1783.5–6.5411,975.734,001,202.9561263516.9137597.272
1793.5–6.5414,022.94,002,38753262716.6119563.673
1803.5–6.5411,697.914,000,107.5744232120.398550.420
1813.5–6.5410,2003,999,70637251217.860694.237
1823.5–6.5413,372.734,001,528.3842212116.394687.529
1833.5–6.5414,253.794,000,845.7646222420.192476.591
1843.5–6.5408,932.64,000,25948232519.294494.7100
1853.5–6.5413,510.133,999,595.5748242417.485388.697
1863.5–6.5412,160.933,999,051.8839211814.13829100100
1873.5–6.5410,732.93,999,78538201810.883241100
1883.5–6.5410,2533,998,4454120218.7132944100
1893.5–6.5413,3093,998,392.546242213.494561.163
1903.5–6.5410,136.164,001,490.7258273118.493599.141
1913.5–6.5409,296.114,002,026.5150262421.0108493.539
1923.5–6.5409,7074,003,851.951282315.08739648
1933.5–6.5408,515.584,002,310.9350282228.238495.320
1943.5–6.5409,9394,003,009.641241725.440578.429
1953.5–6.5407,821.054,002,866.5645261924.751487.9100
1963.5–6.5407,298.54,004,255.6245232212.2102593.420
1973.5–6.5405,9374,003,22441231815.087389.826
1983.5–6.5407,774.754,005,380.148232519.9128691.481
1993.5–6.5407,844.334,005,870.1138261223.750466.731
2003.5–6.5406,914.854,005,671.1555262917.4112598.771
2013.5–6.5405,2424,004,431.748232517.2125593.870
2023.5–6.5405,684.544,005,717.4550272316.294398.145
2033.5–6.5408,025.123,999,907.548262217.781399.181
2043.5–6.5406,400.44,001,1244926232185498.8100
2056.5–9.5406,851.54,009,57047232426.941383.510
2066.5–9.5405,972.94,009,284.251242730.4513.593.812
2076.5–9.5407,487.64,007,947.156302619.1107493.583
2086.5–9.5406,121.14,008,204.754282616.4112595.780
2096.5–9.5407,443.24,008,702.749272213.61054.594.896
2106.5–9.5407,941.44,009,347.755262921.797497.582
2116.5–9.5408,217.94,008,512.234228018.31175100100
2126.5–9.5408,997.44,008,917.854262818.5110598.482
2136.5–9.5408,743.44,009,398.142231917.41086.597.395
2146.5–9.5408,026.14,010,40641281317.41105.591.782
2156.5–9.5409,617.84,011,321.646252118.41225.552.1100
2166.5–9.5408,588.14,011,780.242251719.11035.586.3100
2176.5–9.5409,762.14,012,805.447262121.499590.9100
2186.5–9.5411,184.24,011,787.247262121.91184.591.780
2196.5–9.5411,184.34,012,540.753252818.31034.59085
2206.5–9.5411,247.74,013,522.94728019.7674.510029
2216.5–9.5409,969.94,010,643.856302618.893490.785
2226.5–9.5409,662.94,009,97739211818.2954.582.490
2236.5–9.5409,915.14,009,232.244222216.291491.347
2246.5–9.5410,862.14,009,607.3555282717.485496.180
2256.5–9.5411,994.64,009,216.948232516.1874.598.578
2266.5–9.5411,166.14,010,438.953252821.7924.590.1100
2276.5–9.5412,126.94,010,308.356263012.991496.699
2286.5–9.5411,550.54,011,441.247212615.61088.5100100
2296.5–9.5413,143.34,011,488.941212021.7953.576.990
2306.5–9.5414,064.74,011,692.859283116.6130593.847
2316.5–9.5415,317.34,010,885.464333118.792464.274
2326.5–9.5414,149.34,009,671.456332318.192473.472
2336.5–9.5414,362.54,008,450.666293718.2964.561.655
2346.5–9.5416,171.14,012,249.632151716.1944.552.4100
2356.5–9.5415,609.84,012,761.836201618.2934.542.3100
2366.5–9.5414,549.14,013,63744212311.37814100100
2376.5–9.5417,922.54,010,46243251817.1623.594.940
2386.5–9.5416,961.14,011,319.346222417.21104.574.381
2396.5–9.5416,542.14,010,233.555312420.7116488.7100
2406.5–9.5416,909.44,008,75448242422.7110586.2100
2416.5–9.5406,775.34,007,107.552272517.885498.452
2426.5–9.5405,227.44,007,274.351262518.2854.594.354
2436.5–9.5407,403.94,006,885.345252020.5103494.445
2446.5–9.5407,715.14,007,488.549262320.1974.596.171
2456.5–9.5408,445.44,007,069.448262220.785493.182
2466.5–9.5408,736.74,005,101.954252916.6833.594.183
2476.5–9.5408,248.54,006,694.750272319.29759379
2486.5–9.5407,872.44,006,657.2251252619.4138493.360
2496.5–9.5409,311.94,005,78649262318.21304.597.570
2506.5–9.5408,775.594,006,631.2448232517.697594.969
2516.5–9.5409,973.44,006,13344242017.799592.275
2526.5–9.5409,048.644,007,069.3948242418.4974.590.490
2536.5–9.5409,874.194,007,406.146232315.898593.389
2546.5–9.5408,985.144,007,907.656253118.81084.596.791
2556.5–9.5411,015.24,007,902.745232214.5115593.299
2566.5–9.5418,521.84,007,163.448232513.51134.593.8100
2576.5–9.5409,785.254,008,333.146242218.8953.59443
2586.5–9.5410,998.64,008,878.349272218.61034.592.9100
2596.5–9.5410,109.14,006,834.4442281417.21085.592.159
2606.5–9.5410,419.84,006,00155233217.11204.591.272
2616.5–9.5412,049.024,007,215.9755243115.4130582.389
2626.5–9.5411,2634,006,447.555262918.4102492.781
2636.5–9.5412,664.184,007,328.448252316.41335.596.143
2646.5–9.5411,9904,006,315.249262316.5120589.771
2656.5–9.5413,164.94,008,77945271820.21235.58669
2666.5–9.5412,846.744,006,559.8145242119.6854.594.775
2676.5–9.5412,536.54,005,174.237181920.288498.084
2686.5–9.5415,747.44,005,960.6946222418.489491.2100
2696.5–9.5414,814.84,007,156.943222116.3128542.387
2706.5–9.5415,071.94,004,866.938211716.3122445100
2716.5–9.5417,971.234,001,685.0147252217.2979.542100
2726.5–9.5418,5194,003,222.937211616.8909.536100
2736.5–9.5415,9264,002,87844232117.0931042100
2746.5–9.5416,1674,004,572.343232020.180960100
2756.5–9.5419,2894,001,977.644242017.7109481.498
2766.5–9.5411,924.274,004,297.5956253114.4142596.3100
2776.5–9.5410,482.294,003,953.6300015.6131718.1100
2786.5–9.5412,426.984,002,802.6943232020.51364.593.1100
2796.5–9.5414,3654,003,40858322617.297593.7100
2806.5–9.5411,975.734,001,202.9551242718.11295.597.672
2816.5–9.5414,022.94,002,38781265515.3140595.179
2826.5–9.5411,697.914,000,107.5751282317.51314.594.240
2836.5–9.5410,2003,999,70648282216.8103561.445
2846.5–9.5413,372.734,001,528.3845212417.299693.742
2856.5–9.5414,253.794,000,845.7649262318.6964.582100
2866.5–9.5415,598.54,000,25941231821.089476.5100
2876.5–9.5413,510.133,999,595.5743232013.9873.526.693
2886.5–9.5412,160.933,999,051.8836171913.263210080
2896.5–9.5410,732.93,999,78543212217.0163145100
2906.5–9.5410,2533,998,4454222209.073041100
2916.5–9.5413,3093,998,392.559263315.7943.570.8100
2926.5–9.5410,136.164,001,490.7255262921.897595.979
2936.5–9.5409,296.114,002,026.5153252819.21154.595.242
2946.5–9.5408,932.64,003,851.950282217.594496.849
2956.5–9.5408,515.584,002,310.9345281725.9334.587.938
2966.5–9.5409,9394,003,009.643281526.4554.59333
2976.5–9.5407,821.054,002,866.5600019.7523.525.8100
2986.5–9.5407,298.54,004,255.6241231815.21194.594.740
2996.5–9.5405,9374,003,22446242219.3934.563.738
3006.5–9.5407,774.754,005,380.143212220.5130694.289
3016.5–9.5407,844.334,005,870.1143232019.8493.598.938
3026.5–9.5406,914.854,005,671.1548232518.7994.594.8100
3036.5–9.5405,2424,004,431.750262415.511559673
3046.5–9.5405,684.544,005,717.4548232518.799594.8100
3056.5–9.5408,025.123,999,907.545252022.8105492.691
3066.5–9.5406,400.44,001,12448252322.197569.7100

References

  1. Clayton, C.R.; Matthews, M.C.; Simons, N.E. Site Investigation; Granada: London, UK, 1982. [Google Scholar]
  2. Roy, S.; Bhalla, S.K. Role of geotechnical properties of soil on civil engineering structures. Resour. Environ. 2017, 7, 103–109. [Google Scholar]
  3. Skempton, A. Effective stress in soils, concrete and rocks. Sel. Pap. Soil Mech. 1984, 1032, 4–16. [Google Scholar]
  4. Lupiezowiec, M. The application of c-φ reduction method to estimate the bearing capacity of subsoil. ACEE Archit. Civ. Eng. Environ. 2013, 6, 35–43. [Google Scholar]
  5. Patel, A. Geotechnical Investigations and Improvement of Ground Conditions; Woodhead Publishing: Cambridge, UK, 2019. [Google Scholar]
  6. Bozbey, I.; Togrol, E. Correlation of standard penetration test and pressuremeter data: A case study from Istanbul, Turkey. Bull. Eng. Geol. Environ. 2010, 69, 505–515. [Google Scholar] [CrossRef]
  7. Yusof, N.Q.; Zabidi, H. Reliability of using standard penetration test (SPT) in predicting properties of soil. J. Phys. Conf. Ser. 2018, 1082, 012094. [Google Scholar] [CrossRef] [Green Version]
  8. Cabalar, A.F.; Akbulut, N. Evaluation of Actual and Estimated Hydraulic Conductivity of Sands with Different Gradation and Shape; Springer Plus: Berlin/Heidelberg, Germany, 2016; Volume 5, pp. 1–6. [Google Scholar]
  9. Player, R.S. Geographic information system (GIS) use in geotechnical engineering. In Proceedings of the GeoCongress 2006: Geotechnical Engineering in the Information Technology Age 2006, Atlanta, GA, USA, 26 February–1 March 2006; pp. 1–6. [Google Scholar]
  10. Ahmed, C.; Mohammed, A.; Saboonchi, A. ArcGIS mapping, characterisations and modelling the physical and mechanical properties of the Sulaimani City soils, Kurdistan Region, Iraq. Geomech. Geoeng. 2022, 17, 384–397. [Google Scholar] [CrossRef]
  11. Wan-Mohamad, W.N.; Abdul-Ghani, A.N. The use of geographic information system (GIS) for geotechnical data processing and presentation. Procedia Eng. 2011, 20, 397–406. [Google Scholar] [CrossRef] [Green Version]
  12. Antoniou, A.A.; Papadimitriou, A.G.; Tsiambaos, G. A geographical information system managing geotechnical data for Athens (Greece) and its use for automated seismic microzonation. Nat. Hazards 2008, 47, 369–395. [Google Scholar] [CrossRef]
  13. Cabalar, A.F.; Canbolat, A.; Akbulut, N.; Tercan, S.H.; Isik, H. Soil liquefaction potential in Kahramanmaras, Turkey. Geomat. Nat. Hazards Risk 2019, 10, 1822–1838. [Google Scholar] [CrossRef] [Green Version]
  14. Mancini, F.; Stecchi, F.; Gabbianelli, G. GIS-based assessment of risk due to salt mining activities at Tuzla (Bosnia and Herze-govina). Eng. Geol. 2009, 109, 170–182. [Google Scholar] [CrossRef]
  15. Naji, D.M.; Akin, M.K.; Cabalar, A.F. A comparative study on the VS30 and N30 based seismic site classification in Kahramanmaras, Turkey. Adv. Civ. Eng. 2020, 2020, 1–5. [Google Scholar] [CrossRef]
  16. ESRI Environmental Systems Research Institute. 2005. Available online: https://www.esri.com/ (accessed on 10 February 2022).
  17. Arnous, M.O. Geotechnical site investigations for possible urban extensions at Suez City, Egypt using GIS. Arab. J. Geosci. 2013, 6, 1349–1369. [Google Scholar] [CrossRef]
  18. Chao, Z.; Ma, G.; Zhang, Y.; Zhu, Y.; Hu, H. The application of artificial neural network in geotechnical engineering. In IOP Con-ference Series: Earth and Environmental Science; IOP Publishing: Bristol, UK, 2018; Volume 189, p. 022054. [Google Scholar]
  19. Sulewska, M.J. Applying artificial neural networks for analysis of geotechnical problems. Comput. Assist. Methods Eng. Sci. 2017, 18, 231–241. [Google Scholar]
  20. Nugroho, S.A.; Fernando, H.; Suryanita, R. Estimation of standard penetration test value on cohesive soil using artificial neural network without data normalization. Int. J. Artif. Intell. ISSN 2022, 2252, 8938. [Google Scholar] [CrossRef]
  21. Johora, F.T.; Hickey, C.J.; Yasarer, H. Predicting Geotechnical Parameters from Seismic Wave Velocity Using Artificial Neural Networks. Appl. Sci. 2022, 12, 12815. [Google Scholar] [CrossRef]
  22. Mermerdaş, K.; Arbili, M.M. Explicit formulation of drying and autogenous shrinkage of concretes with binary and ternary blends of silica fume and fly ash. Constr. Build. Mater. 2015, 94, 371–379. [Google Scholar] [CrossRef]
  23. Arbili, M.M.; Ghaffoori, F.K.; Mermerdaş, K. Statistical analysis of the performance of the soft computing based prediction model for shrinkage of concrete including mineral admixtures. ZANCO J. Pure Appl. Sci 2016, 28, 574–579. [Google Scholar]
  24. Ghoreishi, B.; Khaleghi Esfahani, M.; Alizadeh Lushabi, N.; Amini, O.; Aghamolaie, I.; Hashim, N.A.; Alizadeh, S.M. Assessment of geotechnical properties and determination of shear strength parameters. Geotech. Geol. Eng. 2021, 39, 461–478. [Google Scholar] [CrossRef]
  25. Hameed, H. Water Harvesting in Erbil Governorate, Kurdistan Region, Iraq: Detection of Suitable Sites Using Geographic in-Formation System and Remote Sensing. Master’s Thesis, Lund University, Lund, Sweden, 2013. [Google Scholar]
  26. Jassim, S.Z.; Goff, J.C. Geology of Iraq: Dolin; Prague and Moravian Museum: Brno, Czech Republic, 2006; p. 408. [Google Scholar]
  27. Al-Sanjari, M.N.; Al-Tamimi, M.A. Interpretation of water quality parameters for Tigris River within Mosul City by using principal components analysis. Tikrit J. Pure Sci. 2009, 14, 68–74. [Google Scholar]
  28. Das, B.M.; Sivakugan, N. Principles of Foundation Engineering; Cengage Learning: Belmont, CA, USA, 2018. [Google Scholar]
  29. Fernando, H.; Nugroho, S.A.; Suryanita, R.; Kikumoto, M. Prediction of SPT value based on CPT data and soil properties using ANN with and without normalization. Int. J. Artif. Intell. Res. 2021, 5, 123–131. [Google Scholar] [CrossRef]
  30. Chakraborty, A.; Goswami, D. Prediction of slope stability using multiple linear regression (MLR) and artificial neural network (ANN). Arab. J. Geosci. 2017, 10, 385. [Google Scholar] [CrossRef]
  31. Duong, T.V.; Tang, A.M.; Cui, Y.J.; Trinh, V.N.; Dupla, J.C.; Calon, N.; Canou, J.; Robinet, A. Effects of fines and water contents on the mechanical behavior of interlayer soil in ancient railway sub-structure. Soils Found. 2013, 53, 868–878. [Google Scholar] [CrossRef] [Green Version]
  32. Ebrahimi, A. Behavior or fouled ballast. In Railway Track and Structures; Simmons-Boardman Publishing Corporation: New York, NY, USA, 2011; Volume 107. [Google Scholar]
  33. Kim, D.; Sagong, M.; Lee, Y. Effects of fine aggregate content on the mechanical properties of the compacted decomposed gra-nitic soils. Constr. Build. Mater. 2005, 19, 189–196. [Google Scholar] [CrossRef]
  34. Inam, A.; Ishikawa, T.; Miura, S. Effect of principal stress axis rotation on cyclic plastic deformation characteristics of unsaturated base course material. Soils Found. 2012, 52, 465–480. [Google Scholar] [CrossRef] [Green Version]
  35. Wang, Q.; Tang, A.M.; Cui, Y.J.; Delage, P.; Barnichon, J.D.; Ye, W.M. The effects of technological voids on the hydro-mechanical behaviour of compacted bentonite–sand mixture. Soils Found. 2013, 53, 232–245. [Google Scholar] [CrossRef]
  36. Koester, J.P. The influence of test procedure on correlation of Atterberg limits with liquefaction in fine-grained soils. ASTM Geotech. Test. J. 1992, 15, 352–361. [Google Scholar]
  37. Mitchell, J.K.; Soga, K. Fundamentals of Soil Behavior; John Wiley & Sons: New York, NY, USA, 2005. [Google Scholar]
  38. Bogati, K. Ground Improvement by Jet Grouting Techniques. Bachelor’s Thesis, HAMK University of Applied Sciences Construction Engineering, Amberg, Germany, 2019. [Google Scholar]
  39. Carter, M.; Bentley, S.P. Correlations of Soil Properties; Pentech Press Publishers: London, UK, 1991. [Google Scholar]
  40. Likos, W.J.; Song, X.; Xiao, M.; Cerato, A.; Lu, N. Fundamental challenges in unsaturated soil mechanics. Geotech. Fundam. Addressing New World Chall. 2019, 209–236. [Google Scholar] [CrossRef]
  41. Spagnoli, G.; Feinendegen, M. Relationship between measured plastic limit and plastic limit estimated from undrained shear strength, water content ratio and liquidity index. Clay Miner. 2017, 52, 509–519. [Google Scholar] [CrossRef]
  42. Casagrande, A. Notes on the design of the liquid limit device. Geotechnique 1958, 8, 84–91. [Google Scholar] [CrossRef]
  43. Gardner, C.M.; Robinson, D.; Blyth, K.; Cooper, J.D. Soil water content. In Soil and Environmental Analysis; CRC Press: Boca Raton, FL, USA, 2000; pp. 13–76. [Google Scholar]
  44. Das, B.M.; Das, B.M. Advanced Soil Mechanics; Taylor & Francis: New York, NY, USA, 2008. [Google Scholar]
  45. Wu, W.; Mhaimeed, A.S.; Al-Shafie, W.M.; Ziadat, F.; Dhehibi, B.; Nangia, V.; De Pauw, E. Mapping soil salinity changes using re-mote sensing in Central Iraq. Geoderma Reg. 2014, 2, 21–31. [Google Scholar] [CrossRef]
  46. Khudhur, S.M.; Khudhur, N.S. Soil pollution assessment from industrial area of Erbil City. J. Zankoi Sulaimani 2015, 17, 225–238. [Google Scholar] [CrossRef]
  47. Aldefae, A.H.; Mohammed, J.; Saleem, H.D. Digital maps of mechanical geotechnical parameters using GIS. Cogent Eng. 2020, 7, 1779563. [Google Scholar] [CrossRef]
  48. Nam, S.; Gutierrez, M.; Diplas, P.; Petrie, J. Determination of the shear strength of unsaturated soils using the multistage direct shear test. Eng. Geol. 2011, 122, 272–280. [Google Scholar] [CrossRef]
  49. Bapeer, G.B. The study of infiltration rate and atterberg limits of soils in Koi Sanjaq City, Erbil Governorate, Kurdistan Region, North Iraq. Iraqi Bull. Geol. Min. 2011, 7, 41–55. [Google Scholar]
  50. Nareeman, B.J.; Fattah, M.Y. Effect of Soil Reinforcement on Shear Strength and Settlement of Cohesive Frictional Soil. Geomate J. 2012, 3, 308–313. [Google Scholar] [CrossRef]
  51. Gunduz, Z.; Arman, H. Possible relationships between compression and recompression indices of a low-plasticity clayey soil. Arab. J. Sci. Eng. 2007, 32, 179. [Google Scholar]
  52. Qurtas, S.S. Using groundwater levels and Specific Yield to Estimate the Recharge, South of Erbil, Kurdistan Region, Iraq. Acad. J. Nawroz Univ. 2018, 7, 191–196. [Google Scholar] [CrossRef] [Green Version]
  53. Bowles, J. Foundation Analysis and Design; McGraw-Hill Book Company: New York, NY, USA, 1982; 816p. [Google Scholar]
  54. Kumar, R.; Bhargava, K.; Choudhury, D. Estimation of engineering properties of soils from field SPT using random number generation. INAE Lett. 2016, 1, 77–84. [Google Scholar] [CrossRef] [Green Version]
  55. Hasan, A.M.; Mawlood, Y.; Ahmed, A.A.; Ibrahim, H. Correlation of Shear Wave Velocity with SPT-N for a Tower-Building Site at Erbil City. J. Duhok Univ. 2020, 23, 235–245. [Google Scholar] [CrossRef]
  56. Mahmoud, M.A. Reliability of using standard penetration test (SPT) in predicting properties of silty clay with sand soil. Int. J. Civ. Struct. Eng. 2013, 3, 545–556. [Google Scholar]
  57. Myslivec, A.; Kysela, Z. The Bearing Capacity of Building Foundations; Elsevier: Amsterdam, The Netherlands, 2014. [Google Scholar]
  58. Das, B.M.; Sivakugan, N. Fundamentals of Geotechnical Engineering; Cengage Learning: Belmont, CA, USA, 2016. [Google Scholar]
  59. Dauji, S.; Rafi, A. Spatial interpolation of SPT with artificial neural network. Eng. J. 2021, 25, 109–120. [Google Scholar] [CrossRef]
  60. Jasim, M.M.; Al-Khaddar, R.M.; Al-Rumaithi, A. Prediction of bearing capacity, angle of internal friction, cohesion, and plasticity index using ANN (case study of Baghdad, Iraq). Int. J. Civ. Eng. Technol. 2019, 10, 2670–2679. [Google Scholar]
  61. Ïpek, E.; McKee, S.A.; Caruana, R.; de Supinski, B.R.; Schulz, M. Efficiently exploring architectural design spaces via predictive modeling. ACM SIGOPS Oper. Syst. Rev. 2006, 40, 195–206. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of Erbil City, Iraq.
Figure 1. Location of Erbil City, Iraq.
Sustainability 15 04030 g001
Figure 2. Borehole locations in the study area.
Figure 2. Borehole locations in the study area.
Sustainability 15 04030 g002
Figure 3. Flow chart showing in a simple way the methodology adopted in the present study.
Figure 3. Flow chart showing in a simple way the methodology adopted in the present study.
Sustainability 15 04030 g003
Figure 4. Structure of neural network models (a) and (b) to predict SPT-N value.
Figure 4. Structure of neural network models (a) and (b) to predict SPT-N value.
Sustainability 15 04030 g004
Figure 5. Structure of neural network models (c) and (d) to predict Q-ultimate.
Figure 5. Structure of neural network models (c) and (d) to predict Q-ultimate.
Sustainability 15 04030 g005
Figure 6. Fines content at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Figure 6. Fines content at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Sustainability 15 04030 g006
Figure 7. Liquid limit of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Figure 7. Liquid limit of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Sustainability 15 04030 g007
Figure 8. Plastic limit of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Figure 8. Plastic limit of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Sustainability 15 04030 g008
Figure 9. Natural water content of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Figure 9. Natural water content of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Sustainability 15 04030 g009
Figure 10. Cohesion of soils (kPa) at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Figure 10. Cohesion of soils (kPa) at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Sustainability 15 04030 g010
Figure 11. Angle of internal friction of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Figure 11. Angle of internal friction of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Sustainability 15 04030 g011
Figure 12. Compression index, Cc, of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Figure 12. Compression index, Cc, of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Sustainability 15 04030 g012
Figure 13. Rebound index, Cr, of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Figure 13. Rebound index, Cr, of soils at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Sustainability 15 04030 g013
Figure 14. SPT-N values in site at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Figure 14. SPT-N values in site at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Sustainability 15 04030 g014
Figure 15. Correlation of SPT-N values with cohesion of soils in the study area.
Figure 15. Correlation of SPT-N values with cohesion of soils in the study area.
Sustainability 15 04030 g015
Figure 16. Correlation of SPT-N values with angle of internal friction (ϕ) of soils in the study.
Figure 16. Correlation of SPT-N values with angle of internal friction (ϕ) of soils in the study.
Sustainability 15 04030 g016
Figure 17. Ultimate bearing capacity according to Meyerhof’s equation at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Figure 17. Ultimate bearing capacity according to Meyerhof’s equation at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Sustainability 15 04030 g017
Figure 18. Ultimate bearing capacity from SPT-N value at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Figure 18. Ultimate bearing capacity from SPT-N value at depths (a) 1.5–3.5 m. (b) 3.5–6.5 m. (c) 6.5–9.5 m.
Sustainability 15 04030 g018
Figure 19. Correlation of ultimate bearing capacity (Meyerhof) with SPT-N values of soils in study area.
Figure 19. Correlation of ultimate bearing capacity (Meyerhof) with SPT-N values of soils in study area.
Sustainability 15 04030 g019
Figure 20. Correlation of ultimate bearing capacity (SPT-N) with soil cohesion of soils in study area.
Figure 20. Correlation of ultimate bearing capacity (SPT-N) with soil cohesion of soils in study area.
Sustainability 15 04030 g020
Figure 21. Semivariograms/Covariance Modeling.
Figure 21. Semivariograms/Covariance Modeling.
Sustainability 15 04030 g021
Figure 22. Searching Neighborhood Dialog Box of Ordinary Kriging.
Figure 22. Searching Neighborhood Dialog Box of Ordinary Kriging.
Sustainability 15 04030 g022
Figure 23. Model (a): relations between (LL%, PL%, PI%, WC, c, ϕ, Fine content) and SPT-N.
Figure 23. Model (a): relations between (LL%, PL%, PI%, WC, c, ϕ, Fine content) and SPT-N.
Sustainability 15 04030 g023
Figure 24. Model (b): relations between (LL%, PL%, PI%, WC, ϕ, Fine content) and SPT-N.
Figure 24. Model (b): relations between (LL%, PL%, PI%, WC, ϕ, Fine content) and SPT-N.
Sustainability 15 04030 g024
Figure 25. Model (c): relations between (LL%, PL%, PI%, WC, c, ϕ, Fine content) and Q ultimate.
Figure 25. Model (c): relations between (LL%, PL%, PI%, WC, c, ϕ, Fine content) and Q ultimate.
Sustainability 15 04030 g025
Figure 26. Model (d): relations between (LL%, PL%, PI%, WC, ϕ) and Q ultimate.
Figure 26. Model (d): relations between (LL%, PL%, PI%, WC, ϕ) and Q ultimate.
Sustainability 15 04030 g026
Figure 27. Error percentage lines of model (a) and (b) for prediction SPT-N value.
Figure 27. Error percentage lines of model (a) and (b) for prediction SPT-N value.
Sustainability 15 04030 g027
Figure 28. Error percentage lines of model (c) and (d) for prediction Q ultimate.
Figure 28. Error percentage lines of model (c) and (d) for prediction Q ultimate.
Sustainability 15 04030 g028
Table 1. Inputs and output of the present study.
Table 1. Inputs and output of the present study.
NoLL%PL%PI%WC%c kN/m2ϕFine ContentSPT-N Value
kN/m2
Q UL
kN/m2
ASTM D 4318ASTM D2216ASTM 3080ASTM D 6913ASTM D1586-
146222428.351462.17117
240221815.096394.154224
347252218.099591.557225
..........
..........
..........
30448232518.799594.8100292
30545252022.8105492.691296
30648252322.197569.7100288
Min0001229342574
Avarege48.0224.9022.8918.6992.674.2189.4864.69255.09
Max713637291366100100375
SD *7.513.675.592.8719.550.6211.6525.9360.98
SD * = Standard deviation.
Table 2. ANN models for Q-Ultimate prediction.
Table 2. ANN models for Q-Ultimate prediction.
NoModel No.InputOutputTrainingValidationTestingAdjust R2
1Model (c)LL%, PL%, PI%, WC, c, ϕ, Fine contentQ-Ultimate 91.583.88288.79
3Model (d)LL%, PL%, PI%, WC, ϕ, Fine content73.9734.886.9870.8
Table 3. Analysis of models and effect each parameter on the output of prediction.
Table 3. Analysis of models and effect each parameter on the output of prediction.
Model aModel bModel cModel d
ParametersF-Valuep-ValueParametersF-Valuep-ValueParametersF-Valuep-ValueParametersF-Valuep-Value
LL%1.630.204LL%0.030.862LL%0.170.677LL%2.020.157
PL%1.070.303PL%0.510.477PL%0.070.792PL%2.190.142
PI%0.460.497PI%0.270.607PI%1.010.316PI%0.240.627
WC2.520.115WC2.360.126WC1.120.291WC5.950.016
c157.850.000φ0.750.387c320.210.000φ0.030.873
φ10.290.002No.20046.730.000φ28.470.000No.20068.380.000
No.2001.810.180 No.2000.040.843
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Qader, Z.B.; Karabash, Z.; Cabalar, A.F. Analyzing Geotechnical Characteristics of Soils in Erbil via GIS and ANNs. Sustainability 2023, 15, 4030. https://doi.org/10.3390/su15054030

AMA Style

Qader ZB, Karabash Z, Cabalar AF. Analyzing Geotechnical Characteristics of Soils in Erbil via GIS and ANNs. Sustainability. 2023; 15(5):4030. https://doi.org/10.3390/su15054030

Chicago/Turabian Style

Qader, Zhvan Baqi, Zuheir Karabash, and Ali Firat Cabalar. 2023. "Analyzing Geotechnical Characteristics of Soils in Erbil via GIS and ANNs" Sustainability 15, no. 5: 4030. https://doi.org/10.3390/su15054030

APA Style

Qader, Z. B., Karabash, Z., & Cabalar, A. F. (2023). Analyzing Geotechnical Characteristics of Soils in Erbil via GIS and ANNs. Sustainability, 15(5), 4030. https://doi.org/10.3390/su15054030

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop