Thick Clay Research Articles

Analysing Double Shearing Mechanism in Fiberglass Rock Bolting Systems: a Comprehensive Analytical Model and Numerical Simulation Approach

An analytical method and numerical simulation were developed to investigate the shear performance of fiberglass rock bolts (20-tonne and 30-tonne) by conducting sixteen double-shearing tests with both clean and infilled shear interfaces. Following the preparation of the required samples, each test set-up was subjected to different ranges of pretension values. The infilled scenario involved 5 mm thick sandy clay infilled shear interfaces. The results of the double shearing tests unveiled that as pretension increased, so did the confining pressures at the shear interfaces for both clean and infilled joints. Also, an analytical model was developed utilising the Fourier transform, energy balance theory, and linear elastic theory. The result was an empirical relationship that could determine the double shear performance of fibreglass rock bolts in close agreement with the experimental data. Coefficients were incorporated to facilitate model calibration and tuning. Eventually, fast Lagrangian analysis of continua (FLAC) three-dimensional (3D) modelling was utilised to conduct numerical simulations of fibreglass rock bolts subjected to double shearing scenarios. The numerical model was calibrated against experimental data and then extended to conduct a sensitivity analysis on fibreglass rock bolts subjected to double shear test setup variations. Scenarios included rock bolt installation angles, shearing rates, and various host rock strengths. The results revealed that increasing the shear speed from the experimental test baseline yielded substantial displacement increases in the post-failure residual performance of the rock bolts. Changing the installation angle resulted in greater peak shear forces and extended residual zones. The least significant impacts were observed when changing the host rock UCS, suggesting neither rock bolt was drastically impacted by weak or strong host rocks.

Clay resource estimation using a geostatistical approach in Tam Binh district, Vinh Long province, Viet Nam

ABSTRACT Clay is an essential resource for improving the brick and ceramic industries. However, the demand for clay needs to grow. Estimating industrial mineral reserves requires an optimal and unbiased statistical tool. Geographical information systems and geostatistics are two promising tools that could be used to predict clay resources. The study objective was to apply these two tools to compare the performance of the methods for spatial clay resource prediction, considering clay thickness and spatial distribution. The demonstration site is in Tam Binh district, Vinh Long Province, Viet Nam. The results show that the Gaussian variation model at a 600 × 400 m distance is good for sampling and extrapolating clay thickness. This is based on 240 geolocated samples with different densities. Projected clay reserves match measurements, and the spatial statistical method is more accessible and valuable than the traditional method. The interpolated method’s total selected soil borehole number was optimized for clay resource estimation, and geostatistics can help manage and plan clay resources and save money. Future research should examine various spatial interpolation techniques and approaches in management, land use, and topography to improve the accuracy of calculations.

Effective Use of Catastrophe Multicriteria Decision Analysis in Delineating Groundwater Recharge Potential Zones

Multicriteria decision analysis based on Thom’s [1] Catastrophe Theory (CT) has been applied extensively in solving various social, physical and natural sciences problems. It has become a key tool for identifying the groundwater potentiality of an area with a Geographic Information System (GIS). This paper aimed to apply catastrophe theory effectively by utilizing the standardization formulae suggested by Li et al. [2] and the modified formulae provided by Hongwei Zhu [3]. Depending on the nature of a hydrogeological parameter, CT formulae were chosen to ensure that the least important features also get attention during standardization using Hongwei Zhu's formulae where applicable, which were not possible using Li’s formulae. The standardized values of the features of all the themes obtained using four formulae were sorted in ascending order to estimate the normalized values to ensure that no normalized value of a feature exceeds the others having very close but lower standardized values. The effective use of CT technique was employed in a GIS environment to delineate groundwater recharge potential zones of the middle-west part of Kushtia district, Bangladesh, by integrating influential recharge factors, such as land type, slope, drainage density, distance from surface water bodies, soil permeability, surface clay thickness, rainfall, topographic ruggedness index, topographic curvature index, topographic wetness index and topographic position index. The groundwater recharge potential zones were classified as very good (28.76%), good (32.17%), and moderately good (39.05%) for effective CT technique. But in case of the improper use of CT covering area were 17.24%, 54.05% and 28.71% respectively, and the respective most sensitive factors are rainfall and drainage density. Finally, recharge potential zones were validated using groundwater recharge data with a determinant coefficient of 0.92 and 0.84 for effective and improper use of CT technique respectively.

Bearing capacity of modified suction caissons in clay overlying sand under monotonic horizontal loading

The modified suction caisson (MSC) is an innovative foundation for fixing offshore wind turbines, consisting of an internal caisson (IC) and an external skirt (ES). In this study, the Plaxis3D is conducted on the bearing capacity of the MSC in clay overlying sand under monotonic horizontal loading. The parameters of the aspect ratio of the IC, the size of the ES, the clay thickness, and the soil strength are considered. Results show that the horizontal ultimate bearing capacity of the MSC decreases as the length ratio of the ES to IC decreases. Moreover, the width of the ES enhances the horizontal bearing capacity of the MSC higher than the length of the ES. With the increase of the clay thickness, the rotation center of the TSC and MSC exhibits two and three stage change patterns, respectively. When the clay and sand are the inhomogeneous soils, the horizontal ultimate bearing capacity of the MSC is enhanced by 8.96% compared to the homogeneous soils. Additionally, the passive and active earth pressures around the MSC are mainly generated at the front and rear side of the loading direction, respectively. This study will provide theoretical guidance for the design of the MSC.

Wykorzystanie metody radaru penetrującego grunt w badaniu gleby zanieczyszczonej węglowodorami w obszarze Navodari - Rumunia

Ground penetrating radar (GPR) is a very useful geophysical method for use in hydrogeologic and near-surface mapping studies. It can be used to study contaminants in groundwater, subsurface faulting, and underground cavities (natural or man-made), all of which pose potentially dangerous geological hazards. The GPR technique is similar in principle to seismic reflection and sonar techniques. The propagation of the radar signal depends on the frequency-dependent electrical properties of the ground.Electrical conductivity of the soil or rock materials along the propagation paths introduces significant absorptive losses which limit the depth of penetration into the earth formations and is primarily dependent upon the moisture content and mineralization present.Reflected signals are amplified, and transformed to the audio-frequency range, recorded, processed, and displayed. From the recorded display, subsurface features such as soil/ soil, soil/rock, and unsaturated/saturated interfaces can be identified. In addition, the presence of floating hydrocarbons on the water table, the geometry of contaminant plumes, and the location of buried cables, pipes, drums, and tanks can be detected. The GPR data are presented as a two-dimensional depth profile along a scanned traverse line in which the vertical axis is two-way travel time measured in nanoseconds. The location of hydrocarbon contamination in the ground using the GPR method is based mainly on information taken from reflected signals. In the cases investigated in Romania contaminated sites (Navodari area), such signals were very rarely recorded. A long time after spillage, contamination takes the form of plumes with different size and distribution, which depends on the geological and hydraulic properties of the ground. The survey discussed in this paper was carried out using the GPR system-Noggin with two antennas (250 and 500mHz) Data collected were processed using software(EKKO_Project™ GPR Data Analysis) to produce 2D radargram in time scale. The presence of contaminant plumes as well as the water table are observed in the GPR sections at depths approximately of 0.5 to 1.5 m. In the GPR section, the oil contaminated layer exhibits discontinuous, subparallel, and chaotic high amplitude reflection patterns. Promising results were also obtained in the GPR survey where three obvious reflection patterns representing the top sand-silt layer, oil-contaminated zone and, the underlying thick soft clay were detected in all 2D radargrams of the GPR traverse lines.

Numerical analysis of downward progressive landslides in long natural slopes with sensitive clay

Landslides occurring in sensitive clay often result in widespread destruction, posing a significant risk to human lives and property due to the substantial decrease in undrained shear strength during deformation. Assessing the consequences of these landslides is challenging and necessitates robust numerical methods to comprehensively investigate their failure mechanisms. While studies have extensively explored upward progressive landslides in sensitive clays, understanding downward progressive cases remains limited. In this study, we utilised the nodal integration-based particle finite element method (N-PFEM) with a nonlinear strain-softening model to analyse downward progressive landslides in sensitive clay on elongated slopes, induced by surcharge loads near the crest. We focused on elucidating the underlying failure mechanisms and evaluating the effects of different soil parameters and strain-softening characteristics. The simulation results revealed the typical pattern for downward landslides, which typically start with a localised failure in proximity to the surcharge loads, followed by a combination of different types of failure mechanisms, including single flow slides, translational progressive landslides, progressive flow slides, and spread failures. Additionally, inclined shear bands occur within spread failures, often adopting distinctive ploughing patterns characterised by triangular shapes. The sensitive clay thickness at the base, the clay strength gradient, the sensitivity, and the softening rate significantly influence the failure mechanisms and the extent of diffused displacement. Remarkably, some of these effects mirror those observed in upward progressive landslides, underscoring the interconnectedness of these phenomena. This study contributes valuable insights into the complex dynamics of sensitive clay landslides, shedding light on the intricate interplay of factors governing their behaviour and progression.

Effect of Vehicle Cyclic Loading on the Failure of Canal Embankment on Soft Clay Deposit

Road embankments along irrigation canals, constructed on soft Bangkok clay, have always been unstable. Numerous studies have shown that rapid drawdown of water level may be one of the main causes, while vehicle cyclic loading may also contribute to embankment failure. This study aims to investigate the impact of vehicle loading on the failure of embankments built on Bangkok soft clay. The behavior of soft Bangkok clay under vehicle load has been investigated by employing conventional and dynamic triaxial techniques, and finite element method (FEM). This study also examined the effects of soft clay thickness and cyclic loading with different magnitudes and frequencies. The laboratory testing results indicate that the threshold stress of the soft clay is estimated to be approximately three-fourths of the undrained shear strength of the soil. The reduction in effective stress in the soft clay is caused by varied frequencies and thicknesses of the clay. Based on the analysis results, it has been proven that the cyclic loads exerted by vehicles solely are insufficient to cause the embankment to collapse. Nevertheless, the repetitive loading of vehicles may result in a one-quarter decrease in the embankment’s factor of safety.

Prediction of the undrained base heave and post uplift in braced excavation of a typical subway station in thick soft clay

Construction of a typical subway station usually involves a narrow long excavation with a width of approximately 20 m. An undrained base heave is a common problem in deep excavation areas with soft clay, which may induce considerable post uplift and damage the safety of the retaining system. Post uplift involves complex interactions between soil, pile post foundation, and supported struts, and an accurate prediction of post uplift is challenging. This paper aimed to establish a practical method for predicting the base heave and post uplift in the construction of typical subway stations in areas with thick soft clay. Series soil–water coupled finite element method (FEM) analysis was conducted to investigate the characteristics of the base heave with different soil properties and wall stiffnesses. The results showed that the calculated base heave can be fitted by a normalized function, regardless of large variations in soil properties and wall design parameters. The normalized function was then applied to a soil–post–strut interaction analysis model to predict the post uplift at different excavation depths. The application of the proposed method to a real subway station project showed that the predicted results were in good agreement with the field observations.

Analysis of Potential Water Inflow Rates at an Underground Coal Mine Using a WOA-CNN-SVM Approach

The water yield of aquifers increases the risk of water inflow, threatens the safe production of coal mines, and even causes geological disasters and construction hazards. To predict water yield quickly and accurately, multiple composite factors are used to invert unit water inflow rates to judge water yield grade. Taking the typical representative of north China-type coal fields as an example, six factors are selected: aquifer thickness, the radius of influence, normalized drawdown, permeability coefficient, the core rate of drilling holes, and the proportion of clay thickness to the thickness of the lower group. The whale optimization algorithm (WOA)–convolutional neural network (CNN)–support vector machine (SVM) model is established with the unit water inflow rate as the forecast target, and different models are selected for comparison. The water yield zoning map is obtained by bringing the borehole data into the model for prediction. The findings indicate that the root mean square error and average absolute error of the composite predictive model models are 0.0318 and 0.0268, respectively, and the model outperforms alternative models. The predicted water yield zoning aligns well with the actual conditions, offering a novel paradigm for water yield assessment.

The Badenian/Sarmatian (Middle Miocene) boundary in the Central Paratethys (Kreminna, western Ukraine): Foraminiferal and palynological evidence

The Badenian/Sarmatian boundary in the Paratethys basin, that marks the transition from normal marine to restricted semi-marine conditions due to isolation of the basin from the world ocean at the onset of Sarmatian time, is still far from being fully understood. The Kreminna section is located at the northeastern margin of the Carpathian Foreland Basin (Central Paratethys) in the Medobory Hills region. The Miocene deposits that overlie here the Upper Cretaceous substratum comprise the >1 m thick upper Badenian marls and clays passing upwards into ~4 m thick Sarmatian marly limestones with intercalations of marls, clays and limestones, and >2 m thick limestones in the uppermost part of the exposure. Fifty-three species of benthic foraminifera and four species of planktonic foraminifera have been recorded. Six benthic foraminiferal assemblages are composed almost exclusively of calcareous forms; agglutinated taxa are practically lacking. Elphidium spp., miliolids, Lobatula lobatula and Ammonia spp. are the most common calcareous benthic foraminifera in the material studied. Planktonic foraminifera are represented only by species of Globigerina and occur rarely in the lowermost part of the section. A characteristic feature of palynofacies is the very low proportion of land-derived elements – sporomorphs and cuticles, which suggests a sedimentary setting without terrestrial influx, and taxonomical impoverishment of dinoflagellate cyst assemblages, which are either monospecific or consist mainly of two to three species: in a majority of samples, assemblages with Polysphaeridium zoharyi and P. subtile occur. Most δ18O values range from –1.5 to ~+0.5‰ VPDB and most δ13C values are between 0 and +2‰ VPDB. In general, the δ13C curve mirrors the δ18O changes in the section. The Badenian/Sarmatian boundary is placed at the level where the Cibicidoides ungerianus Assemblage is replaced by the Elphidium fichtelianum Assemblage. At the boundary, planktonic foraminifers and most abundant stenohaline benthic foraminifera disappeared. Bottom waters were well-oxygenated both in the latest Badenian and earliest Sarmatian in the Kreminna location.

Comprehensive Study on Deep Excavation and Deformation Control of Thick Soft Clay Layers in Urban Complex Environments

Comprehensive Study on Deep Excavation and Deformation Control of Thick Soft Clay Layers in Urban Complex Environments

Evaluation of Consolidation Degree from CPTu using Rahardjo (2016) Method – Case Study of Consolidating Soil in East Kalimantan

A geotechnical investigation comprised 13 CPTu and 23 dissipation testing were performed on a site in East Kalimantan due to landslide triggered by 15 to 30m thick soft clay underlying 75m high overburden embankment. The investigation covers an area of about 250m x 275m wide on the area next to the toe of the embankment. From geological map, the site situated at borders of Balikpapan Formation (Tmbp) and Kampung Baru Formation (Tpkb) with Alluvium (Qa) formation from Heliosen period nearby. Measurement from penetration showed high pore-water pressure higher than the hydrostatic pressure. Result from dissipation testing showed that the soft clay is still consolidating with residual excess pore pressure (uf) still exist. Rahardjo (2016) Method was developed using CPTu data as its basis especially soft clay data to determine overconsolidation ratio (OCR). However, it can also be used to determine the degree of consolidation. The method uses Pore Pressure Ratio (Bq) value, excess pore pressure divided by net cone resistance – Bq = Δu/(qt – σv0), which was obtained when performing CPTu test. The formula proposed was 1/(1.2Bq+0.1). The method stated that the value of Bq=0.75 equals to OCR=1 which is showing a normally consolidated clay. Bq value higher than 0.75 will show a degree of consolidation of a consolidating soil and for Bq value lower than 0.75 will show an OCR value of overconsolidated soil. Applying Rahardjo (2016) Method to the dissipation data performed showed an agreement on soft clay where the soil is still consolidating. However, when applied to over consolidated soil near surface this method will show a greater OCR value.

Tectonic Subsidence on the East China Coast Recorded by Magnetic Properties of Pliocene Red Clay in the Yangtze Delta

Thick red clay in northern China contains rich information about the uplifting of the Qinghai–Tibet Plateau, the drying process of Asian inland, the East Asian monsoon changes, and global cooling since 22 Ma. In comparison, the red clay widely distributed in southern China is generally much younger (<1 million years), thus limiting the paleoclimate and paleoenvironment reconstruction over a longer geological time. We conducted a comprehensive magnetic investigation on the Pliocene red clay of the core LQ11, located in the Yangtze Delta, to reveal its paleoclimate and paleoenvironment implications for the eastern China coast. Our results revealed that the Pliocene red clay in the Yangtze Delta has higher S-ratio and lower HIRM (Hard isothermal remanent magnetizations) values than Quaternary vermiculate red clay of hot–humid climate origin in southern China. This indicates a weaker transformation from maghemite to hematite during the process of pedogenesis. The lack of net-like white veins in the Pliocene red clay also indicates a relatively low intensity of pedogenesis. We believe that the Pliocene red clay, which is presently 250 m below the mean sea level, was formed in high-altitude topography before the Quaternary period, where paedogenic intensity was remarkably low. This finding shows rapid tectonic subsidence occurring on the eastern China coast since the late Pliocene and enriches the theoretical research on paleoenvironment reconstruction based on red clay.

Assessing the relationship between contributing factors and sinkhole occurrence in Mexico City

Sinkholes pose a significant hazard in Mexico City (CDMX), causing substantial economic damage. While the link between sinkhole formation and groundwater extraction has been studied, specific mechanisms vary by site. Our overall aim is to characterize the phenomenon of sinkholes in CDMX. To achieve this, we create a database with 13 influencing factors, including population density, well density, distance to faults, fractures, roads, streams, elevation, slope, clay thickness, lithology, subsidence rate, geotechnical zones, and soil texture. Sinkhole locations were obtained from CDMX’s Risk Atlas (2017–2019). We shaped a susceptibility map based on statistical regression methods derived from applying linear regression models. For the susceptibility map, results showed that 40% of variables are significantly correlated with sinkhole density. Despite the regression model explained 24% of sinkhole density variability, it helped choosing variables for the susceptibility map that correlate better (89.7%). Hence, we identified that the northeast CDMX was the most susceptible zone. Therefore, the compound assessment of environmental factors is useful for the evaluation of susceptibility maps to identify prone factors for the generation of sinkholes. This framework provides relevant information for better use of the territory throughout the development of public policies.

Investigation of the performance enhancement of building-integrated photovoltaic system using evaporative porous clay applied in different building's directions

The integration of photovoltaic (PV) modules with building's facades and roofs results in significant reduction in the PV module's efficiency owing to its temperature rising, as well as increasing in the building's indoor thermal gain. Thus, this study proposed cost-effective and eco-friendly passive cooling technique utilizing evaporative porous clay material to minimize the integrated PV module's temperature. An experimentally validated numerical study is being performed using different clay thickness and porosity to estimate the expected enhancement in the thermal and electrical performance of the building-integrated photovoltaic (BIPV) system. The experimental tests were carried out using prototype of small rooms integrated with PV panels installed with a water-saturated porous clay material to validate the presented transient numerical heat and mass transfer model. After that, the proposed system is being compared with another passive cooling mechanism using phase change material (PCM), and the results are then analyzed and compared. The findings showed a maximum reduction of 19.1 °C and 12.8 °C of peak PV temperature when the porous clay and PCM are utilized, while the building's side temperature reduced by around 11.2 °C and 10.9 °C, respectively. Further, porous clay achieved up to 9.8 % improvement in PV system's efficiency, with a maximum water consumption rate reached around 3.35 L/h.m2. Moreover, applying such a porous clay cooling approach resulted in a maximum reduction of 53.8 % in the building's annual thermal load. In the outcomes reported here, porous clay with the least thickness and highest porosity achieved more feasible cooling effects compared with PCM.

Integrated geophysical assessment of a municipal waste disposal site for its geological suitability in terms of the underlain material

Dumpsites are major sources of groundwater pollution as a result of leachates that drain out of the decomposed waste. If there are no underlain materials that could serve as a seal to stop the percolation of leachate, it finds its way to the groundwater. A properly designed landfill is expected to have a high leachate curtailment capacity to limit groundwater pollution. A suitable landfill is expected to have a specific thickness of clay which acts as a natural filter. This study aims to determine the subsurface material and the leachate curtailment of Oke–Diya dumpsites. Very low frequency-electromagnetic method was adopted as a reconnaissance survey, after which electrical resistivity and multichannel analysis of surface waves (MASW) methods were carried out. The resistivity values obtained were used to determine the lithological units of the study area while the MASW was employed to determine the seismic wave arrival times which was processed to obtain the shear wave velocities of the subsurface. The rigidity moduli were also obtained from the shear wave velocities, from which the lithological units of the subsurface were inferred. The integrated method appeared to be the ideal tool to characterize the dumpsite and adjudge the leachate curtailment capacity. The methods corroborated each other. Oke-Diya dumpsite, from the results, revealed the study area had low leachate curtailment capacity and should be evacuated.

Kapasitas Tiang Tunggal pada Tanah Lempung Berdasarkan Hasil Uji Model Fisik

Soft clay soil as a construction subgrade requires reinforcement from the pile to reduce settlement and increase bearing capacity. Alternative materials that can be used are the bamboo pile or wooden pile. The purpose of this study is to find out how much the bearing capacity of the pile is compared to concrete pile. The research was conducted using physical model tests in the laboratory. Piles with a diameter of 2 cm with different lengths are driven into the soft clay and then tested using load plates. Single pile capacity is determined from the relationship between the load and the settlement. The results showed that the capacities of bamboo and wooden piles were 1.193 and 1.086 times higher than the capacity of concrete pile, respectively. This can be caused because concrete pile has a greater weight than bamboo pile and wooden pile. The weight of the concrete increases the load, so the settlement is higher and the pile capacity is lower. Pile lengths that produce more effective capacities are obtained for piles with a minimum depth of 3/5 of the thickness of the soft clay.

Effectiveness of Jet-Grouted Slab on Excavation: A Case Study

Jet-grouted slab is a possible solution to reduce the wall deformation induced by the excavation, especially for cases in soft clay. By replacing a layer of a certain thickness of the soft clay inside the excavation fully by chemical or cement grout, it could limit the lateral displacement of the wall to remain the serviceability of the excavation. In this paper, a 16.7 m deep excavation using semi top-down construction method in Taipei was selected to examine the behaviour of the excavation. The jet-grouted slab was adopted in the excavation and then the effectiveness of jet-grouted slab was evaluated. Except the grouted slab, the pit was retained by a 1.2 m thick, 38 m deep reinforcement concrete diaphragm wall with 3-level of steel struts and slabs. A series of finite element analysis was conducted using soft-soil creep and hardening soil with small strain models to carry out after-event predictions and the results were verified by field observation data. It was found that this method could effectively reduce wall deformation up to 33%. Moreover, the enhancement of young’s modulus of the jet-grouted slab does not affect the wall deflection significantly. The jet-grouted slab also was simulated by adopting two different constitutive models namely linear elastic and Mohr-Coulomb (MC) to have a better understanding of the performances by the use of jet-grouted slab. It shows that both constitutive models give similar results.

Study on the Analysis of Pile Foundation Deformation and Control Methods during the Excavation of Deep and Thick Sludge Pits

This study aims to apply performance-based safety-assessment methods to the monitoring and numerical simulation of excavation engineering projects in order to comprehensively enhance engineering risk management and decision support. In this paper, a deep excavation project in Hefei with thick silty clay layers was studied. The analysis included the surface settlement, the deformation of support structures, the vertical and horizontal displacements of pile tops, axial forces in steel braces, settlement, and the horizontal displacement of a gravity retaining wall on the south side of the excavation using field-monitoring data. A refined three-dimensional finite element model was established to further analyze the distribution of uplift displacement at the bottom of the excavation, horizontal displacement, and bending moments of piles based on simulation results. The research findings indicate that phased excavation can reduce the spatial extent of disturbance to the surrounding soil caused by excavation. Additionally, the closer the location to the excavation and the thicker the underlying silty clay layer, the faster the rate of settlement change and the greater the surface settlement. The spatial structure formed by steel braces and pile foundations effectively reduced the horizontal displacement of the engineering piles. The study’s use of field monitoring and finite element simulation provided valuable insights into the deformation of support structures and the response of the surrounding soil to excavation, confirming the rationality and applicability of the support structure in this paper. The proposed method can serve as a reference for similar complex stratum excavation design and construction. The performance-based safety assessment is introduced, and the monitoring data, numerical simulation results, and performance targets are comprehensively analyzed to provide a reliable scientific basis for engineering decision making.

Evaluation of sedimentary aquifer in tehsil Zehri area using VES technique, Balochistan Pakistan

A geophysical resistivity survey was conducted in Zehri tehsil in the district of Khuzdar for the metering of subsurface hydrogeological properties of aquifers. The Schlumberger electrode arrangement was employed with the vertical electrical sounding (VES) method for subsurface exploration. Two profile lines were designated for the survey. The distance between the two profile lines was about 1500 ​m and each profile contained four VES stations. The results of profile one show the aquifers were encountered in fractured limestone of the Jamburo Formation of Paleocene age at an average depth of 127 ​m, with an average resistivity of 39.7 ​Ωm, and the average aquifer thickness was 54.025 ​m. The result of profile two shows the aquifers were encountered only in VES four in the sandstone lithology of the Ziddi Formation with a resistivity of 79.9 ​Ωm, at a depth of 120 ​m with, and a thickness of 40 ​m. There was no aquifer encountered in the remaining all VES. The main lithology of these VES is thick impervious compacted limestone, shale, and clay of Jurassic age from the Ziddi Formation. The current tube wells were also rapidly depleting in the study area. The main factor behind the decline of aquifers in the study area is an uncontrolled installation of tube wells and continuous pumping of tube wells for agricultural purposes badly puts stress on aquifers and directly deteriorates the quality and quantity of groundwater.