Apparent Formation Factor for Leachate-Saturated Waste and Sediments: Examples from the USA and China

The formation factor relates bulk resistivity to pore fluid resistivity in porous materials. Understanding the formation factor is essential in using electrical and electromagnetic methods to monitor leachate accumulations and movements both within and around landfills. Specifically, the formation factor allows leachate resistivity, the degree of saturation, and, possibly, even the hydraulic conductivity of the waste to be estimated from non-invasive surface measurements. In this study, apparent formation factors are computed for three landfills with different types of waste as well as sediments contaminated by landfill leachate. Resistivity soundings at the closed Mallard North landfill in suburban Chicago (Illinois, USA) mapped leachate surfaces that were confirmed by monitoring wells. The resistivity of leachate-saturated waste from resistivity sounding inversions was then divided by the leachate resistivity values measured in-situ to compute apparent formation factors (Fa) ranging from 1.6 to 4.9. A global Fa of 3.0±1.9 was computed for the entire monitored portion of this landfill At a nearby mixed laboratory waste landfill, a 2D inverted resistivity section was used to compute an Fa of 2.9. Finally, a distinctly different Fa value of 10.6±2.8 was computed for leachate-saturated retorted oil and organic compounds. The Fa for aquifers containing contaminated groundwater fall in the same range as aquifers with normal groundwater, 1.7-3.9. However, models from inverted sounding curves over these contaminated areas exhibit unusually low resistivity layers, which may be diagnostic of contamination.     

Equivalent-linear seismic analyses of MSW landfills using DEEPSOIL

Responses of municipal solid waste (MSW) landfills during earthquakes are gaining worldwide importance due to the devastating nature of earthquake on landfills. Apart from the post-earthquake safety and serviceability issues which pose environmental and public health problems, other important concerns are related to the behavior of closed landfills during and after earthquake. In present study, one dimensional (1-D) equivalent-linear analysis was carried out to model the behaviour of MSW landfills subjected to seismic base accelerations using the DEEPSOIL software. Influence of foundation types, height and stiffness of MSW landfills and seismic base accelerations on the seismic responses in terms of surface accelerations, normalized stresses (i.e., shear stress/effective vertical stress) and spectral amplification are evaluated. The results showed that height and stiffness of landfills, type of foundation and amount of seismic base acceleration and period play important role in evaluating the seismic responses of MSW landfills. Assumption of constant unit weight and shear wave velocity for landfills underestimates maximum horizontal acceleration (MHA), normalized shear stresses and spectral amplification at the top of landfills. Landfill models with smaller heights (up to 40 m high) showed higher amplification ratio for low seismic base accelerations with mean period near to that of soil and landfill resonance for all sites. A complex behavior was observed at higher seismic base accelerations due to non-linear behavior of landfill materials.     

3D behaviour of contamination in landfill sites using 2D resistivity/IP imaging: case studies in Portugal

Aiming at defining a valid spatial contamination model, resistivity and induced polarization (IP) measurements were used to investigate contamination plumes in the vicinity of two municipal landfills (Ovar and ílhavo). Previous geophysical surveys and underground water samples confirmed the contamination. However 2D resistivity/IP surveys enabled in obtaining a more accurate spatial model. The Ovar survey consisted of two profiles with nine Wenner soundings each; the ílhavo survey was carried out along two individual lines using a Wenner standard pseudo-section. In both situations, negative IP values were found associated with positive IP values, which can be explained mainly by 2D or 3D geometric effects caused by the presence of the conductive plumes. The data were modelled using a 2D inversion program (RES2DINV) and the resulting resistivity and chargeability distributions were displayed as pseudo-sections. The resistivity and chargeability pseudo-sections define the contamination plumes and the sedimentary structure. These case studies illustrate the advantages of 2D resistivity/IP surveys for the mapping of shape and dimension of contamination associated with landfills.     

Application of electrical resistivity tomography to evaluate the variation in moisture content of waste during 2 months of degradation

Moisture content is considered to play an important role in the degradation of municipal solid waste (MSW). The present study was carried out to test the use of electrical resistivity tomography (ERT) to estimate the daily variations in the volumetric moisture content (VMC) during waste degradation. The degradation of MSW in a laboratory-scale reactor was monitored daily by ERT and time domain reflectometry (TDR) over a period of 55 days after the waste was landfilled. The bulk electrical resistivity, normalised to a standard temperature of 25 °C, was compared to the VMC obtained by TDR. The relationship between the bulk electrical resistivity and the VMC was not obvious in the upper unit, because there was a bias in the positions of the ERT and TDR measurements due to the settling of the waste. However, the bulk electrical resistivity depended on the VMC in the lower unit of the reactor. According to a power equation fitted in the lower unit, the 2-D distributions of the VMC are presented. Variations in the VMC reflected the processes of waste degradation and leachate transportation. Furthermore, the volume of water stored in the lower unit of the reactor was calculated and found to be consistent with that estimated from the gravimetric moisture content during waste sampling. The data showed that ERT could be used to estimate the variations in moisture content in the initial period following the landfilling of waste.     

Geophysical assessment of subsurface coastal sediments and their engineering implications

Surface geophysics and a priori information were employed to delineate the subsurface geology at Idi-oro in Abijo, Ibeju Lekki area of Lagos, Nigeria for foundation investigation purpose. Resistivity measurement was conducted using 1-D and 2-D resistivity probing techniques. The resistivity measurements were made with ABEM tetrameter model SAS 1000 system. The 1-D vertical electrical resistivity sounding data were obtained using the Schlumberger electrode array while the 2-D resistivity data were obtained using the dipole–dipole array. The interpreted results revealed three to five subsurface geological layers. This is made up of the top soil with resistivity values that vary from 132.4 to over 2,313.5 Ω?m and thickness values that range from 0.3 to 4.8 m, the fine sand with resistivity values that vary from 221.0 to 3,032.7 Ω?m and thickness values that range from 0.4 to 5.5 m, the medium sand with resistivity values that vary from 202.8 to 1,247.7 Ω?m and thickness values that range from 4.9 to 58.4 m. On the other hand, the clayey sand has the resistivity values that vary from 146.1 to 1,744.0 Ω?m and thickness values that vary from 2.2 to 26.3 m, while the coarse sand has resistivity values that vary from 238.3 to 14,313.9 Ω?m but with no thickness value because the current terminated in this layer. The resistivity data correlate well with borehole logs. On the whole, it is concluded that the investigated area has competent sand layer that can support medium to giant engineering structures with resistivity values that vary from 202 to 14,314 Ω?m and thickness values that vary from 0.8 to 58.4 m.     

Geochemical characterisation and geophysical mapping of Landfill leachates: the Marozzo canal case study (NE Italy)

An integrated approach including geochemical analyses and geophysical investigations, aimed at characterizing the chemical composition of leachate and at highlighting the extent and magnitude of the plume in terms of spatial distribution is presented. The reason for the specific induced polarisation behaviour over landfills was not well understood. In this work we propose an explanation for this behaviour, which is being used to detect the presence of possible leaks and to track and delimit the pollution plume for several years.As a case study the example of an abandoned municipal solid waste landfill, located in Comacchio County, NE Italy, is presented, where the operative conditions for any resistivity monitoring method were to be excluded because of the high salinity of the phreatic aquifer.     

Hydrogeophysical investigation and estimation of groundwater potentials of the Lower Palaeozoic to Precambrian crystalline basement rocks in Keffi area, north-central Nigeria, using resistivity methods

The use of resistivity sounding and two-dimensional (2-D) resistivity imaging was investigated with the aim of delineating and estimating the groundwater potential in Keffi area. Rock types identified are mainly gneisses and granites. Twenty-five resistivity soundings employing the Schlumberger electrode array were conducted across the area. Resistivity sounding data obtained were interpreted using partial curve matching approach and 1-D inversion algorithm, RESIST version 1.0. The 2-D resistivity imaging was also carried out along two traverses using dipole–dipole array, and the data obtained were subjected to finite element method modeling using DIPRO inversion algorithm to produce a two-dimensional subsurface geological model. Interpretation of results showed three to four geoelectrical layers. Layer thickness values were generally less than 2 m for collapsed zone, and ranged from 5 to 30 m for weathered bedrock (saprolite). Two major aquifer units, namely weathered bedrock (saprolite) aquifer and fractured bedrock (saprock) aquifer, have been delineated with the latter usually occurring beneath the former in most areas. Aquifer potentials in the area were estimated using simple schemes that involved the use of three geoelectrical parameters, namely: depth to fresh bedrock, weathered bedrock (saprolite) resistivity and fractured bedrock (saprock) resistivity. The assessment delineated the area into prospective high, medium and low groundwater potential zones.     

Geostatistical integration using 2-D electrical resistivity and 3-D gravity methods for detecting cavities in a karst area

Two-dimensional (2-D) electrical resistivity and three-dimensional (3-D) gravity explorations were undertaken to estimate the 3-D distribution of karst cavities at the area of Yongweol-ri in Korea where ground subsidence has occurred. Although the gravity method is a low-cost way of analyzing a 3-D continuous structure, its vertical resolving power is poor. In contrast, the electrical resistivity method can provide a 2-D and/or 3-D subsurface structure with a much higher lateral and vertical resolution than the gravity method. Accordingly, geostatistical methods and density information were used to enhance the 2-D resistivity structure revealed by the electrical resistivity method into a 3-D structure. The assumptions are, first, that each data set senses the same underlying geological structure in terms of different material properties and, secondly, that two different material properties are correlated locally or globally throughout the entire target area. As a result, the distribution of limestone cavities can be estimated under the assumption that they are mostly filled with groundwater and clayey soils and have abnormally low levels of resistivity and density. The estimated distribution corresponds with the grouting, borehole imaging and monitoring data. In this example, it can be seen that the integration analysis of 2-D electrical resistivity and 3-D gravity methods is a very powerful tool for 3-D subsurface imaging and that the method can provide enhanced imaging capabilities for 3-D cavities.     

Controlled Source Audio Magneto Telluric (CSAMT) studies for uranium exploration in Durgi area,Palnad sub-basin,Cuddapah basin,India

Cuddapah basin is known for hosting unconformity proximal uranium deposits viz., Lambapur, Peddagattu, Chitirial and Koppunuru along the northern margin of the basin. It is well known that these deposits are mostly associated with basement granitoids in Srisailam Sub-basin, and with cover sediments in Palnad subbasin where basement topography and fault/fracture system influence the fluid flow causing basement alteration and ore deposition. Geological setup, surface manifestation of uranium anomalies and association of the hydro-uranium anomalies near Durgi area in southern part of the Palnad sub-basin, have prompted detail investigation by geophysical methods to probe greater depths. Controlled Source Audio Magneto Telluric (CSAMT) survey conducted over five decades of frequency (0.1-9600 Hz) delineated the various lithounits of Kurnool and Nallamalai Groups along with their thicknesses as there exist an appreciable resistivity contrast. Interpretation of CSAMT sounding data are constrained by resistivity logs and litholog data obtained from the boreholes drilled within the basin indicated three to four layered structure. Sub-surface 2-D and 3-D geo-electrical models are simulated by stitching 1-D layered inverted resistivity earth models. Stitched 1-D inverted resistivity sections revealed the unconformity between the Kurnool Group and Nallamalai Group along with basement undulations. The faults/fractures delineated from the CSAMT data corroborated well with the results of gravity data acquired over the same area. Simulated 3-D voxel resistivity model helped in visualising the faults/fractures, their depth extent, thickness of the Banganapalle quartzite and basement configuration. Integrated interpretation of CSAMT, gravity and borehole data facilitated in delineating the unconformity and the structural features favourable for uranium mineralisation in deeper parts of the Palnad sub-basin.     

Resistivity imaging survey to delineate subsurface seepage of hydrocarbon contaminated water at Karbala Governorate,Iraq

This research is an attempt to accomplish a 3-D resistivity imaging survey, which was carried out near a water well contaminated with hydrocarbon materials in Karbala governorate. Two-dimensional resistivity imaging measurements were collected along four parallel profiles, using a Wenner array with electrode spacing of 1 m. The RES3DINV program was used to invert the apparent resistivity data. The results displayed a resistivity distribution of the subsurface in a three-dimensional volume. Thus, both the horizontal and vertical extents of the contaminated zone were displayed. This technique revealed a low resistivity zone at depth ranges from 3 to 6 m in the investigation area, but the seepage starts at depth ranges between 2 and 3 m and continues down depth (may be to the groundwater level). This low resistivity zone is the most likely location for a subsurface seepage of contaminated water. It is clear that the sufficient measurement points along 2-D lines in a small area can increase the 3-D imaging resolution, and nearly real 3-D imaging can be achieved, when the size of subsurface anomaly compared with the electrode spacing (a) of the Wenner array is taken into consideration.     

Stability Analyses of Municipal Solid Waste Landfills with Decomposition

A bioreactor landfill is operated to enhance refuse decomposition, gas production, and waste stabilization. Some of the potential advantages of bioreactor include rapid stabilization of waste, increased landfill gas generation, gain in landfill space, enhanced leachate treatment, and reduced post closure maintenance period. Due to the accelerated decomposition and settlement of solid waste, bioreactor landfills are gaining popularity as an alternative to the conventional Subtitle D landfills. However, the addition of leachate to accelerate the decomposition changes the physical and engineering characteristic of waste and therefore affects the geotechnical characteristics of waste mass. The changes in the physical and mechanical characteristics of solid waste with time and decomposition are expected to affect the shear strength of waste mass. The objective of this paper is to analyze the stability of solid waste slopes within the bioreactor landfills, as a function of time and decomposition. The finite element program PLAXIS is used for numerical modeling of bioreactor landfills. Stability analysis of bioreactor landfills was also performed using limit equilibrium program STABL. Finally the results from finite element program PLAXIS and limit equilibrium program STABL are compared. GSTABL predicted a factor of safety of more than 1 in all the cases analyzed, whereas PLAXIS predicted a factor of safety of less than 1 at advanced stages for a slope of 2:1. However, the interface failures between solid waste and landfill liners have not been considered in this paper.     

Geophysical-geochemical investigation of fire-prone landfills

This paper deals with the integration of electrical resistivity tomography and geochemical methods for studying four different fire-prone landfills. Landfill gas composition (CH₄, H₂S, O₂, CO, CO₂) and subsurface temperature were measured with the constant net 50 × 50 m from the depth 10–60 cm. 28 electrical resistivity tomography lines were surveyed, while Wenner and Sclumberger electrode arrays were employed for all measurements. At the studied sites the landfill gas and temperature measurements mapped gas and temperature anomalies over underground fire sources. 2D electrical resistivity tomography lines, performed over these anomalies, showed these fire sources as high-resistivity zones. The joint employment of the electrical imaging and geochemical survey seems to be a useful tool in carrying out diagnostic investigations at fire-prone landfills.     

Integrated use of hydrochemistry and resistivity methods in groundwater contamination caused by a recently closed solid waste site

 In order to investigate the groundwater contamination by solid waste disposal using both hydrochemical and geophysical methods, the Halkalı (I˙stanbul) solid waste disposal site which was closed in 1994 was investigated. The disposal site lies on a ridge between two valleys filled with alluvium. A total of six boreholes were drilled on two lines across the Menekşe valley adjacent to the Halkalı site. Groundwater samples collected from these boreholes were analyzed for various contaminant parameters. The results indicate that TDS and chloride concentrations decrease horizontally away from the waste site whereas they increase with depth. Electrical soundings carried out at 12 locations yielded high resistivity values at the upstream part of Menekşe valley while lower values were obtained from the locations near the leachate seepage points. Received: 11 November 1997 · Accepted: 23 February 1998     

The use of vertical electrical sounding resistivity method for the location of low salinity groundwater for irrigation in Chaj and Rachna Doabs

A geoelectrical resistivity survey using vertical electrical sounding (VES) was conducted at Chaj Doab (land between rivers Jhelum and Chenab, Pakistan) and Rachna Doab (land between rivers Chenab and Ravi, Pakistan), with the objective of investigating groundwater conditions. A total of 90 sites were selected with 43 sites in Chaj and 47 sites in Rachna Doabs. The resistivity meter (ABEM Terrameter SAS 4000, Sweden) was used to collect the VES data by employing a Schlumberger electrode configuration, with half current electrode spacings (AB/2) ranging from 2 to 180 m and the potential electrode (MN) from 1 to 40 m. The field data were interpreted using the Interpex IX1D computer software and the resistivity versus depth models for each location was estimated. The outputs of subsurface layers with resistivities and thickness presented in contour maps and 3-D views by using SURFER software were created. A total of 102 groundwater samples from nearby hydrowells at different depths were collected to develop a correlation between the aquifer resistivity of VES and the electrical conductivity (EC) of the groundwater and to confirm the resulted geophysical resistivity models. From the correlation developed, it was observed that the groundwater salinity in the aquifer may be considered low and so safe for irrigation if resistivity >45 Ω m, and marginally fit for irrigation having resistivity between 25 and 45 Ω m. The study area has resistivities from 3.9 to 2,222 Ω m at the top of the unsaturated layer, between 1.21 and 171 Ω m, in the shallow aquifers, and 0.14–152 Ω m in the deep aquifers of the study area. The results indicate that the quality of groundwater is better near the rivers and in the shallow layers compared to the deep layers.     

水平电偶源r方向电场值应用的理论研究

根据准静态极限条件下,水平电偶极源在均匀半空间表面产生的r方向电场分量E r表达式,分别在一维和三维条件下,对E r进行数值模拟并迭代计算出相应视电阻率,讨论E r的观测范围和应用优势。在一维条件下,通过对四种典型层状模型(两层D、G,三层H、K)的数值计算,并与相同条件下的E x视电阻率进行对比,结果表明:E r比E x受测量方位的影响更小;在三维条件下,选取均匀半空间存在异常体的模型,采用积分方程法对E r进行数值计算,与相同条件下的E x视电阻率进行对比,结果表明:当均匀半空间存在低阻异常体时,E r对低阻体y方向边界的反应比E x效果较好。     

Electrical imaging of the groundwater aquifer at Banting,Selangor, Malaysia

A geophysical study was carried out in the Banting area of Malaysia to delineate groundwater aquifer and marine clay layer of the alluvial Quaternary deposits of Beruas and Gula Formations. The Beruas Formation is formed by peat and clayey materials as well as silt and sands, whereas the Gula Formation consists of clay, silt, sand and gravels. Both Formations were deposited on top of the Carboniferous shale of the Kenny Hill Formation. A 2-D geoelectrical resistivity technique was used. Resistivity measurement was carried out using an ABEM SAS 4000 Terrameter. The 2-D resistivity data of subsurface material for each survey line was calculated through inverse modelling and then compared with borehole data. The resistivity images of all the subsurface material below the survey lines show similar pattern of continuous structure of layering or layers with some lenses with resistivity ranging from 0.1 to 50 Ωm. The upper layer shows resistivity values ranging from 0.1 to 10 Ωm, representing a clay horizon with a thickness up to 45 m. The second layer with depth varies from 45 to 70 m below surface and has resistivity values ranging from 10 to 30 Ωm. Borehole data indicate coarse sand with some gravels for this layer, which is also the groundwater aquifer in the study area. The lowermost layer at a depth of 70 m below ground level shows resistivity values ranging from 30–50 Ωm and can be correlated with metasedimentary rocks consisting of shale and metaquartzite.     

Magnetotelluric images of the crustal structure of Chyulu Hills volcanic field, Kenya

Electromagnetic experiments were conducted in 1995 as part of a multidisciplinary research project to investigate the deep structure of the Chyulu Hills volcanic chain on the eastern flank of the Kenya Rift in East Africa. Transient electromagnetic (TEM) and broadband (120–0.0001 Hz) magnetotelluric (MT) soundings were made at eight stations along a seismic survey line and the data were processed using standard techniques. The TEM data provided effective correction for static shifts in MT data. The MT data were inverted for the structure in the upper 20 km of the crust using a 2-D inversion scheme and a variety of starting models. The resulting 2-D models show interesting features but the wide spacing between the MT stations limited model resolution to a large extent. These models suggest that there are significant differences in the physical state of the crust between the northern and southern parts of the Chyulu Hills volcanic field. North of the Chyulu Hills, the resistivity structure consists of a 10–12-km-thick resistive (up to 4000 Ω m) upper crustal layer, ca. 10-km-thick mid-crustal layer of moderate resistivity (50 Ω m), and a conductive substratum. The resistive upper crustal unit is considerably thinner over the main ridge (where it is ca. 2 km thick) and further south (where it may be up to 5 km thick). Below this cover unit, steep zones of low resistivity (0.01–10 Ω m) occur underneath the main ridge and at its NW and SE margins (near survey positions 100 and 150–210 km on seismic line F of Novak et al. [Novak, O., Prodehl, C., Jacob, A.W.B., Okoth, W., 1997. Crustal structure of the southern flank of the Kenya Rift deduced from wide-angle P-wave data. In: Fuchs, K., Altherr, R., Muller, B., Prodehl, C. (Eds.), Structure and Dynamic Processes in the Lithosphere of the Afro-Arabian Rift System. Tectonophysics, vol. 278, 171–186]). These conductors appear to be best developed in upper crustal (1–8 km) and middle crustal (9–18 km) zones in the areas affected by volcanism. The low-resistivity anomalies are interpreted as possible magmatic features and may be related to the low-velocity zones recently detected at greater depth in the same geographic locations. The MT results, thus, provide a necessary upper crustal constraint on the anomalous zone in Chyulu Hills, and we suggest that MT is a logical compliment to seismics for the exploration of the deep crust in this volcanic-covered basement terrain. A detailed 3-D field study is recommended to gain a better understanding of the deep structure of the volcanic field.     

Dynamic Properties of Municipal Solid Waste in Bioreactor Landfills with Degradation

Bioreactor landfills are operated to enhance refuse decomposition, gas production, and waste stabilization. The major aspect of bioreactor landfill operation is the recirculation of collected leachate back through the refuse mass. Due to the presence of additional leachate and accelerated decomposition, the characteristics of Municipal Solid Waste (MSW) in bioreactor landfills are expected to change. About 50% of the continental United States comes under the designated seismic impact zone. The federal regulations have focused increase attention on seismic design of solid waste fills, and have mandated that the solid waste landfills located in the seismic impact zones should be designed to resist the earthquake. Accordingly, assessment of dynamic properties of landfills is one of the major geotechnical tasks in landfill engineering. In order to understand the changes in dynamic properties of bioreactor waste mass with time and decomposition, four small scale bioreactor landfills were simulated in laboratory and samples were prepared to represent each phase of decomposition. The state of decomposition was quantified by methane yield, pH, and volatile organic content (VOC). A number of Resonant Column (RC) tests were performed to evaluate the dynamic properties (stiffness and damping) of MSW. The test results indicated that the normalized shear modulus reduction and damping curves are significantly affected by the degree of decomposition. The shear modulus increased from 2.11 MPa in Phase I to 12.56 MPa in Phase IV. The increase was attributed to the breakdown of fibrous nature of solid waste particles as it degrades. Therefore, considering MSW properties to be uniform throughout the bioreactor landfill is not a reasonable assumption and the shear modulus reduction curves should be evaluated based on the degree of MSW decomposition, rather than the sample composition itself.     

Comparison study of methane emissions from landfills with different landfill covers

Atmospheric methane, a more effective heat-trapping gas than CO₂ that may affect climate change, has its greatest man-made source in the US from municipal solid waste (MSW) landfills. Consequently, the wise management of landfills can reduce these greenhouse gas emissions to the atmosphere. Methane from modern MSW landfills built with composite covers is frequently vented directly to the atmosphere. Biofiltration of landfill gas could oxidize CH₄ to CO₂ and water. Methane oxidation in old landfills with conventional soil covers can be effective in reducing the amount of CH₄ emitted. In this study, comparison of methane emissions from three different landfill covers was conducted. Methane emissions from old landfills constructed with conventional soil covers, modern MSW landfills constructed with composite covers, and modern MSW landfills constructed with composite covers plus biofilters were calculated using the calculated CH₄ oxidation rates. The results showed that an average of only 14% of the generated CH₄ was emitted from landfills with modern composite covers plus biofilters, and an average of 85% of the generated CH₄ was emitted from landfills with conventional covers when 100% of the generated CH₄ emissions to the atmosphere from landfills with modern composite covers was assumed. By comparing the CH₄ emission rates from three different landfill types, the use of a properly sized biofilter should be an effective technique to reduce CH₄ emissions from landfills across the USA and potentially in many other areas of the world.     

Monitoring of wastewater percolation in unsaturated sandy soil using geoelectrical measurements at Gabal el Asfar farm,northeast Cairo,Egypt

The use of wastewater for irrigation in sandy soil increases the pollution risk of the soil and may infiltrate to the shallow groundwater aquifer. In such environment, some important parameters need to be obtained for monitoring the wastewater in the unsaturated zone over the aquifer. These parameters include clay content, heterogeneities of the upper soils, depth to the aquifer and the variations of groundwater quality. In the present work, the efficiency of DC resistivity method in forms of 1-D and 2-D measurements was studied for wastewater monitoring in the Gabal el Asfar farm, northeast of Cairo, Egypt. Forty-one Schlumberger soundings (VES) were performed then followed by three pole-dipole 2-D profiles along some considered regions within the area. The resistivity measurements were integrated with the boreholes, hydrogeological and hydrochemical (surface and groundwater samples) information to draw a clear picture for the subsurface conditions. The obtained results were presented as cross sections and 3-D visualization to trace the clay intercalations within the unsaturated zone. In addition, a vulnerability map was created using the obtained results from 1-D Schlumberger survey and confirmed with the 2-D resistivity profiling. The obtained results have shown that the 2-D resistivity imaging technique is a powerful tool for mapping the small-scale variability within the unsaturated zone and the wastewater infiltration. However, limitations of resistivity techniques were observed in the area with limited resistivity contrast such as thin clay layers with brackish water background. Under that condition, the measured pattern of resistivity distributions depends on the applied electrode array, electrode spacing and using the available geological information during the inversion process.