Geophysical and Geochemical Characterization of Solidwaste Dumpsite: A Case Study of Chowa Gujar, Peshawar (Part of Indus Basin)

: Open and non-engineered dumping is a typical method for solid waste disposal in most cities of Pakistan. This practice of waste dumping poses a serious threat to the surrounding ecosystem and human population due to the release and transport of decomposed organic matter, i.e., leachate from dumpsite into the groundwater. The present study was conducted over a non-engineered and open dumpsite (Chowa Gujar), located in the outskirts of the highly populated city of Peshawar by using integrated geophysical techniques such as electrical resistivity tomography (ERT) and ground-penetrating radar (GPR) and geochemical techniques. The main goal was to delineate the characterization and depth of buried waste, to map the subsurface extension of contaminant plumes towards agricultural land and groundwater table and the concentration of heavy metals (HMs) in dump and agricultural soil. Geophysical results showed that the thickness of buried waste (predominantly composed of domestic waste) was around 4 m and the leachate plumes have percolated to the adjacent agricultural land. A range of heavy metals (mg/kg) such as Cr (20.5–26.6), Cd (2.6–5.7), Pb (0.35–21.25), Ni (2.5–53.05), Cu (29–68.3), Zn (45.7–77), and Co (18.9–23.2) have been found in the agricultural land adjacent to the dumpsite. The ﬁndings demonstrated that combined use of ERT and GPR successfully characterize the buried waste and spread of pollutant plumes spatially and vertically from Chowa Gujar dumpsite. The anomalous geophysical signatures were conﬁrmed by geochemical characterization. The movement of leachate plumes towards agricultural land and groundwater table and the concentration of HMs in soil show that Chowa Gujar dumpsite is a potential source of contamination not only to the surrounding population but also to the agricultural land, surface (Bara River), and subsurface water bodies. In the study region, there is an urgency to take remediation and mitigation measures to reduce the level of pollution created by the dumpsite.


Introduction
Municipal solid waste (MSW) management is an internationally identified environmental issue in major populated areas, especially in developing countries, Therefore it has been highlighted in United Nations Sustainable development goals (SDG's); the main aim of these SDG's is to reduce such threats by 2030 [1]. The world's urban population is estimated to exceed 6 billion by 2045, specifically in the countries of the developing world [2]. Pakistan's most recent census, conducted in 2017, counted 208 million people, The specific objectives of the present study are (1) to determine geometrical dimensions and to characterize the buried waste, (2) to map the subsurface distribution of leachate plumes towards groundwater table and agricultural land, and (3) to investigate the possibility of contamination of the nearby agricultural soil due to seepages of leachate from this dumpsite. This study directly addresses some of the Sustainable Development Goals like Goal #3 (regarding good health) and Goal #6 (regarding clean water and sanitation). This study is also designed to explore the potential harmful impacts of solid waste The specific objectives of the present study are (1) to determine geometrical dimensions and to characterize the buried waste, (2) to map the subsurface distribution of leachate plumes towards groundwater table and agricultural land, and (3) to investigate the possibility of contamination of the nearby agricultural soil due to seepages of leachate from this dumpsite. This study directly addresses some of the Sustainable Development Goals like Goal #3 (regarding good health) and Goal #6 (regarding clean water and sanitation).
This study is also designed to explore the potential harmful impacts of solid waste open dumping site on the agricultural soils and potable water resources of a big city such as Peshawar by using the geophysical as well as geochemical techniques for environmental management and sustainability.

The Hydrogeological Setting of the Study Area
The study area (Chowa Gujar), located southeast of Peshawar (33 • 59 48.50" N, 71 • 37 13.91" E), is a part of the Peshawar Basin ( Figure 2) which is an intermountain basin formed around 2.3 million years ago [23]. Peshawar Basin has been divided into four broad geological units as (1) metamorphic and igneous rocks of northern mountains, (2) the Mesozoic rocks of south and southwestern part, (3) the Tertiary rocks of central and southeastern part, and (4) Upper Tertiary and Quaternary fill of intermountain basin [24].
Sustainability 2022, 14, x FOR PEER REVIEW 4 of 19 open dumping site on the agricultural soils and potable water resources of a big city such as Peshawar by using the geophysical as well as geochemical techniques for environmental management and sustainability.

The Hydrogeological Setting of the Study Area
The study area (Chowa Gujar), located southeast of Peshawar (33°59′48.50″ N, 71°37′13.91″ E), is a part of the Peshawar Basin ( Figure 2) which is an intermountain basin formed around 2.3 million years ago [23]. Peshawar Basin has been divided into four broad geological units as (1) metamorphic and igneous rocks of northern mountains, (2) the Mesozoic rocks of south and southwestern part, (3) the Tertiary rocks of central and southeastern part, and (4) Upper Tertiary and Quaternary fill of intermountain basin [24]. Peshawar Basin is surrounded by Atock-Chirat mountain ranges in the south, Swat Mountain ranges in the southeast, and Khyber mountain ranges in the northwest. A central part of the Peshawar Basin is covered with Quaternary sediments. Physiographical Quaternary deposits are classified as Flood Plain, Lacustrine, and Piedmont deposits [26]. Flood Plain deposits are probably deposited by Kabul and Indus Rivers. Piedmont deposits are observed along the foothills at the margins of Peshawar Basin. During the middle Pleistocene, Lacustrine deposits were formed, because of overflow blockage of the Indus River [27].
Hydro-geologically Peshawar Basin is classified into three types of aquifer based on main lithologies of the area such as Peshawar Piedmont, Peshawar Lacustrine, and Peshawar Flood Plain aquifers [28]. The aquifer with a groundwater level of 24 m is located down-gradient of the Chowa Gujar dumpsite. The lithological column of a borehole was drilled up to a depth of 27 m at a distance of about 1.2 km from the dumpsite, and shows Peshawar Basin is surrounded by Atock-Chirat mountain ranges in the south, Swat Mountain ranges in the southeast, and Khyber mountain ranges in the northwest. A central part of the Peshawar Basin is covered with Quaternary sediments. Physiographical Quaternary deposits are classified as Flood Plain, Lacustrine, and Piedmont deposits [26]. Flood Plain deposits are probably deposited by Kabul and Indus Rivers. Piedmont deposits are observed along the foothills at the margins of Peshawar Basin. During the middle Pleistocene, Lacustrine deposits were formed, because of overflow blockage of the Indus River [27].
Hydro-geologically Peshawar Basin is classified into three types of aquifer based on main lithologies of the area such as Peshawar Piedmont, Peshawar Lacustrine, and Peshawar Flood Plain aquifers [28]. The aquifer with a groundwater level of 24 m is located down-gradient of the Chowa Gujar dumpsite. The lithological column of a borehole was drilled up to a depth of 27 m at a distance of about 1.2 km from the dumpsite, and shows a 6 m thick top layer of clay followed by a 3 m strata of gravel boulder with a very high resistivity log value. At a depth of 9-15 m, thin lenses of clay and gravel are encountered alternatively, while at depth of 15-21 m sand with minor clay has been observed. At a depth of 21-27 m there is saturated strata gravel boulder with low resistivity log value as compared to dry strata, below which gravel sand strata is encountered (Figure 3).
Sustainability 2022, 14, x FOR PEER REVIEW 5 of 19 a 6 m thick top layer of clay followed by a 3 m strata of gravel boulder with a very high resistivity log value. At a depth of 9-15 m, thin lenses of clay and gravel are encountered alternatively, while at depth of 15-21 m sand with minor clay has been observed. At a depth of 21-27 m there is saturated strata gravel boulder with low resistivity log value as compared to dry strata, below which gravel sand strata is encountered (Figure 3). Groundwater flow is a valuable indicator for determining the flow patterns of potential contaminants in the groundwater system [29]. Therefore, the topography trend was considered an important factor during the designing of the geophysical survey. In the study area, the groundwater flow direction is from the north towards the east and from the southwest towards the eastern region [30], which is a replica of the topographic trend in the area (Figure 3).

Materials and Methods
Geophysical and geochemical techniques are useful tools for determining the geology and hydrology and geochemistry of the surface and subsurface horizons. Integration of different techniques provides a more synergistic interpretation of the sub-surface data, which frequently leads to a more accurate model of the subsurface complex structures and processes.

Electric Resistivity Tomography
ERT is a non-invasive geophysical technique that is employed at the ground surface by introducing an electric current into the soil, measuring the potential difference response, and then calculating the apparent electrical resistivity [31]. A high concentration of charged particles facilitates the electrolytic conduction process which enables the ERT to perform well.
The ERT technique has widely been used in landfill investigations that involve contaminated soils and aquifers due to percolated constituents of leachate. The conductive Groundwater flow is a valuable indicator for determining the flow patterns of potential contaminants in the groundwater system [29]. Therefore, the topography trend was considered an important factor during the designing of the geophysical survey. In the study area, the groundwater flow direction is from the north towards the east and from the southwest towards the eastern region [30], which is a replica of the topographic trend in the area (Figure 3).

Materials and Methods
Geophysical and geochemical techniques are useful tools for determining the geology and hydrology and geochemistry of the surface and subsurface horizons. Integration of different techniques provides a more synergistic interpretation of the sub-surface data, which frequently leads to a more accurate model of the subsurface complex structures and processes.

Electric Resistivity Tomography
ERT is a non-invasive geophysical technique that is employed at the ground surface by introducing an electric current into the soil, measuring the potential difference response, and then calculating the apparent electrical resistivity [31]. A high concentration of charged particles facilitates the electrolytic conduction process which enables the ERT to perform well. The ERT technique has widely been used in landfill investigations that involve contaminated soils and aquifers due to percolated constituents of leachate. The conductive nature of dumpsite material to the near-surface soil media and aquifers makes ERT an efficient geophysical tool in such types of studies [32].
The ERT measurements were conducted along with four profiles (ERT-1 to ERT-4) with 1 m and 2 m electrode spacing depending on the available space. The details of ERT profiles positions and acquisition parameters are given in Figure 1B and Table 1 respectively. Wenner-Schlumberger configuration was adopted by different researchers to estimate the depth of buried waste and map contaminant plumes, because of its high signal to noise ratio, sensitivity to both horizontal and vertical changes, and good penetration depth [10,21,33]. The RES2DINV software was used for data processing. The ERT data were inverted based on a least-squares algorithm to obtain a two-dimensional distribution model [34]. The resistivity raw data were mostly analyzed for five iterations because they reveal a more logical subsurface geological model with less root mean square (RMS) error. For a fair analogy, the color scale for all inversion models is the same.

Ground Penetrating Radar (GPR)
Ground Penetrating Radar emits high-frequency electromagnetic waves (1 MHz to 1000 MHz) [35]. GPR method has been applied to a variety of problems such as detection of voids in the subsurface and geological layer boundaries, location of installations, and archaeological and forensics studies [36,37]. The GPR applications to monitor landfill have been discussed by [38]. For better interpretation, it is usually used with the integration of electrical methods [39].
In the present study, GPR data were acquired in continuous mode MALA GPR system equipped with 50 MHz frequency un-shielded antenna. The acquisition parameters (time window, sampling frequency) were selected according to the type of antenna and site conditions. Time window was selected 30 ns, sampling frequency was set at 500 MHz (ten times the antenna frequency).
A REFLEXW program [40] was used to process the data. Different processing schemes were applied to improve the quality of results. For instance, to remove the effect of time misalignment and ringing effect, static correction and background filter were applied to the data. To intensify the amplitude linear and exponential gain functions were applied. In the last step subsurface velocity of 0.2 m/ns was acquired via the hyperbolae overlap method and used for time to depth conversion.

Sample Collection and Geochemical Analysis
Soil samples were collected randomly from the dumpsite and the nearby agricultural land. Five representative soil samples were collected from (0-20 cm) open and buried dump waste site while five samples were collected from adjacent agricultural land and five samples were collected from referenced area. These soil samples were collected according to the United States Environmental Protection Agency soil sampling protocols [41].
The organic matters from the soil samples were hand-picked and then sieved to 2-mm mesh. These were then pulverized in the vibrating disc mill. About 1 g of each powdered sample was placed in the Teflon beakers, 10 mL of hydrofluoric acid (HF) of analytical grade was added and heated on a hot plate in a fume hood. After heating for a while, 20 mL of aqua regia (3 HCL/1 HNO 3 ) of analytical grade was added and kept on heat. After complete dryness, then 20 mL HCL (2 N) was added and heated for a while and then the digested contents were filtered and diluted up to 50 mL. All the treated samples were then transferred in clean polythene bottles for further analysis. Final solutions were analyzed through an atomic absorption spectrometer (Perkin Elmer AA700) under standardized conditions.

Geo-Accumulation Index
The geogenic or anthropogenic accumulated pollution loads in soil are measured using the geo-accumulation index and have been calculated using the formula below.
where C n is a concentration of metals in the target area and B n is the concentration of metals in referenced area. The factor 1.5 is used to account for the possible variations in the referenced area values.

Enrichment Factor
The enrichment factor (EF) was determined for the level of elements enrichment in soil It can be calculated through the following formula [42]: where C n is the element concentration in the target area, C m is the lowest possible concentration in the target area of an element, C b is the element concentration in the referenced area, and C x is the lowest possible concentration of an element in the reference area.

Waste Characterization and Leachate Plumes Mapping Using ERT and GPR Imaging
Profile ERT-1 was acquired on the surface of the dumpsite to determine the extent of waste deposit and to detect sub-surface pathways of a plume. The inverted model ( Figure 4A) shows the waste-contained zone along with ERT-1 profile from 0 to 40 m characterized by 31 Ωm to 50 Ωm resistivity value within upper 4 m. Further 40 m to 60 m zone with resistivity value of 25 Ωm within~1 m is defined as natural ground. However, a homogenous resistivity zone having values about 50 to 80 Ωm is observed underneath the natural ground may be indicating old or less saturated domestic waste material marked by a black dotted line. Another waste-contained cell is observed in the N-E part of ERT-1 with a very high resistivity of about 80 to 125 Ωm extended a few meters below the ground surface. An important anomalous zone characterized by extremely low resistivity ≤8 Ωm is observed at depth of~6 m. This zone is interpreted as contaminant plumes, and it is extended up to 14 m in the sub-surface. The NE dipping anomaly and trend of elevation ( Figure 3) indicate that the plumes have their origin in the SW portion of ERT-1 and may probably migrate in the NE direction over time in the future. The zone between waste contained zone and leachate plumes is known as the intermittent zone, which is observed between the depth of 4 m to 6 m and is the lower limit of solid dump waste and a transition to a high conductive zone.  GPR section ( Figure 4B) reveals the continuous reflections on the top relating to the surface cover. The very near-surface of dumpsite was not imaged based on radar reflections except the boundaries of different zones (domestic waste, old buried or less saturated domestic waste and commercial waste) compared to ERT-1 which helped in distinguishing several features within this depth. There are substantial reflections within the top 4 m and laterally these strong reflections are prominent between 0 to 40 m and then 60 to 80 m. Furthermore, ref. [14] indicated that dumpsite produce discontinuous intense scattering regions on radargram. According to [43] strong diffraction hyperbola is produced because of the high reflectivity of buried waste and the high heterogeneity of buried waste disturb the soil media. The strong signals at different depths underneath those areas covered by domestic and commercial wastes are attributed to the targets (rubbles and metallic wastes) within a waste deposit. GPR section ( Figure 4B) reveals the continuous reflections on the top relating to the surface cover. The very near-surface of dumpsite was not imaged based on radar reflections except the boundaries of different zones (domestic waste, old buried or less saturated domestic waste and commercial waste) compared to ERT-1 which helped in distinguishing several features within this depth. There are substantial reflections within the top 4 m and laterally these strong reflections are prominent between 0 to 40 m and then 60 to 80 m. Furthermore, ref. [14] indicated that dumpsite produce discontinuous intense scattering regions on radargram. According to [43] strong diffraction hyperbola is produced because of the high reflectivity of buried waste and the high heterogeneity of buried waste disturb the soil media. The strong signals at different depths underneath those areas covered by domestic and commercial wastes are attributed to the targets (rubbles and metallic wastes) within a waste deposit.
Sub-surface section ERT-1 ( Figure 4A) and GPR-1 ( Figure 4B) have a favorable correlation in defining the thickness of the waste-containing zone as well as characterization of the waste. The studied dump thickness is about 4 m, below which resistivity value drop less than 8 Ωm and serious radar signal attenuation is observed, that is probably related to plumes of conductive leachate. The 0-40 m zone can be referred to as domestic waste with resistivity range 31 to 50 Ωm and some prominent GPR reflections are observed up to the depth of 4 m. No significant reflection is observed on the top section,~1 m of radargram between 40 to 60 m, which can be referred to as natural ground while below this depth GPR reflection is observed can be referred to as old buried domestic waste or less saturated domestic waste with resistivity range of 50 to 80 Ωm. Between 60 to 75 m zone concentration of GPR signal reflections is high as compared to 0 to 60 m zone of the radargram and high resistivity values 80 to 127 Ωm are also observed exactly in this zone, suggesting that this waste cell is filled with some commercial waste. Classification of waste based on resistivity variation and the maximum depth of each waste-contained zone is mentioned in Table 2. The subsurface cross-section ( Figure 4C) depicts the movement of the leachate plume toward the low gradient side of the study area.
Following the leachate accumulation point on the surface of low land adjacent to a dumpsite boundary profile, ERT-2 was acquired. The purpose of this profile was to see the leachate concentration on the very near surface. This section ( Figure 5A) is located close to the leachate accumulation point and the depth of leachate percolation is very shallow,~2 m as compared to ERT-1 ( Figure 4A). The near-surface resistivity of this section is <5 Ωm (white dotted line) lower than the ERT-1. Therefore, it can be concluded that the leachate mineralization is high at the shallow depth close to the leachate accumulation point as compared to the leachate of the waste-containing zone. Below the leachate saturated zone between the depth of 2-5 m, there are dry strata of gravel. Similar very low resistivity patterns have been reported for leachate saturated zones in a controlled environment [44,45].
GPR-2 ( Figure 5B) reveals that at the top of radargram signal attenuated significantly. High conductive nature of leachate plume is assumed to be responsible for severe attenuation. GPR signal dies out as it interacts with a highly conductive contaminated body, therefore no reflection is observed below the signal attenuation zone [46]. A similar pattern of subsurface contamination generated by dumpsites was also observed by [37,47,48]. In the current study, some diffraction signals are observed towards the agricultural landside, which could be due to uncontaminated land. GPR-2 ( Figure 5B) reveals that at the top of radargram signal attenuated significantly. High conductive nature of leachate plume is assumed to be responsible for severe attenuation. GPR signal dies out as it interacts with a highly conductive contaminated body, therefore no reflection is observed below the signal attenuation zone [46]. A similar pattern of subsurface contamination generated by dumpsites was also observed by [37,47,48]. In the current study, some diffraction signals are observed towards the agricultural landside, which could be due to uncontaminated land.
Subsurface cross-section based on ERT-2 and GPR-2 sections ( Figure 5C) revealed that the concentration of leachate plumes from the accumulation point is high towards the Bara Riverside as compared to the agricultural land. Subsurface cross-section based on ERT-2 and GPR-2 sections ( Figure 5C) revealed that the concentration of leachate plumes from the accumulation point is high towards the Bara Riverside as compared to the agricultural land.
Profile ERT-3 was acquired from leachate accumulation point towards agricultural land on the eastern side. The purpose of this profile was to map the spatial and vertical distribution of plumes and to delineate the potential leachate migration routes beyond the leachate accumulation point toward agricultural land. The physical distance between the leachate accumulation point and agricultural land is 80 m. ERT-3 section ( Figure 6A Profile GPR-3 was about 160 m long and the location of a profile is the same as the location of ERT-3. The signal attenuation is observed from 0 to 90 m on radargram while the diffraction signal is observed between the distance of 90 to 104 m, which could be related to the near-surface abnormalities or uncontaminated geological strata which can also be correlated with the ERT-3 section ( Figure 6A). Signal attenuation from 104 to 114 m along GPR-3 is observed which could be due to a leachate lens that is consistent with the very low conductive zone of the ERT-3 section.
Subsurface cross-section based on ERT-3 and GPR-3 ( Figure 6A,B) depicts leachate plumes that have penetrated up to 8 m with a lateral extension of about 30 m within the agricultural land and it has not yet reached the groundwater table.
Profile ERT-4 was acquired within agricultural land adjacent to open dump waste. Section ERT-4 ( Figure 7A) depicts that a low resistivity anomalous zone is observed Profile GPR-3 was about 160 m long and the location of a profile is the same as the location of ERT-3. The signal attenuation is observed from 0 to 90 m on radargram while the diffraction signal is observed between the distance of 90 to 104 m, which could be related to the near-surface abnormalities or uncontaminated geological strata which can also be correlated with the ERT-3 section ( Figure 6A). Signal attenuation from 104 to 114 m along GPR-3 is observed which could be due to a leachate lens that is consistent with the very low conductive zone of the ERT-3 section.
Subsurface cross-section based on ERT-3 and GPR-3 ( Figure 6A,B) depicts leachate plumes that have penetrated up to 8 m with a lateral extension of about 30 m within the agricultural land and it has not yet reached the groundwater table.
Profile ERT-4 was acquired within agricultural land adjacent to open dump waste. Section ERT-4 ( Figure 7A) depicts that a low resistivity anomalous zone is observed between 12 m and 55 m laterally with a vertical extension of more than 14 m. The low resistivity anomalous zone is interpreted as leachate concentrated area in the shallow sub-surface, which is the reason for the attenuation in the upper section of the radargram ( Figure 7B). Radar reflections are observed at the start and end of the radargram due to uncontaminated land. The source of contamination along these profiles could be the adjacent open dump waste. GPR and ERT results demonstrate that open dumpsite has resulted in polluting the agricultural land and may hit the shallow groundwater table in the future. To validate the near-surface saturation of leachate, a small portion of survey line ERT-4 and GPR-4 was dug out up to a depth of 0.5 cm for visual observation ( Figure 1G).
Based on resistivity sections and borehole data, a specific range of resistivity for different lithology and leachate plume in the study area is given in Table 3: To validate the near-surface saturation of leachate, a small portion of survey line ERT-4 and GPR-4 was dug out up to a depth of 0.5 cm for visual observation ( Figure 1G).
Based on resistivity sections and borehole data, a specific range of resistivity for different lithology and leachate plume in the study area is given in Table 3:  55.09 mg/kg, and 18.81 mg/kg, respectively. The trend of heavy metals in dump soil appears as Cu > Zn > Ni > Co > Cr > Pb > Cd and in agricultural soil as Zn > Cu > Ni > Cr > Co > Pb > Cd. In the case of land use type, the concentration of heavy metals (HMs) in dump soil is greater than that of agricultural soil ( Figure 8). However, the concentration of all the HMs in the target area was found higher than that of the reference area.   55.09 mg/kg, and 18.81 mg/kg, respectively. The trend of heavy metals in dump soil appears as Cu > Zn > Ni > Co > Cr > Pb > Cd and in agricultural soil as Zn > Cu > Ni > Cr > Co > Pb > Cd. In the case of land use type, the concentration of heavy metals (HMs) in dump soil is greater than that of agricultural soil (Figure 8). However, the concentration of all the HMs in the target area was found higher than that of the reference area.   The HMs concentrations were compared with the Soil Environmental Quality Standard (SEQS) such as Canadian Council of Ministers of Environment (CCME, 2007) [49], and Chinese Environmental Protection Agency (CEPA, 1995) [50], whereby Cd, Ni, and Cu concentration in the samples from dumpsite and agricultural site exceeded the quality standards of CCME, 2007 and CEPA, 1995 while the concentration of Cr, Co, Pb, and Zn in the soils of both the sites was found to be below the soil quality criteria of CCME, 2007 and CEPA, 1995. These results show that the soil has been contaminated with heavy metals as a result of leachate migration through unfenced and unlined waste dumpsite.
The geo accumulation index of HMs in the study area is given in Table 6. According to [51], geo accumulation index values are classified as G-0 (0 < Igeo ≤ 1) unpolluted, G-1 (1 < Igeo ≤ 2) as slightly polluted, G-2 (2 < Igeo ≤ 3) as moderately polluted, G-3 (3 < Igeo ≤ 4) as moderate to severely polluted, G-4 (4 < Igeo ≤ 5) as severely polluted, and G-5 (Igeo > 5) as severe to extremely polluted. The trend of mean Igeo index in the study area was Ni > Pb > Cd > Cu > Co > Zn > Cr. Heavy metals in the study area could be categorized as Cr, Zn, and Co in the unpolluted category; Cu in slightly polluted; Cd and Pb in moderate polluted; and Ni in moderate sever polluted category Table 6. The enrichment factor is given in the Table 6. Six categories were recognized on the basis of enrichment factor Table 5   On the basis of statistical analysis (Figure 9), such as geo accumulation index and enrichment factor, it has been noticed that the soil of the study area has been affected by soild waste. The solid waste is composed of mainly domestic waste, and a small portion is commercial as confirmed by geophysical results in Section 4.1. To confirm the source of commercial waste, various commercial activities and operational units such as pharmaceutical companies, steel industry, leather, cement, clothes industry, and food processing and beverages industry in the premises of the study area were observed during the reconnaissance survey. On the basis of statistical analysis (Figure 9), such as geo accumulation index and enrichment factor, it has been noticed that the soil of the study area has been affected by soild waste. The solid waste is composed of mainly domestic waste, and a small portion is commercial as confirmed by geophysical results in Section 4.1. To confirm the source of commercial waste, various commercial activities and operational units such as pharmaceutical companies, steel industry, leather, cement, clothes industry, and food processing and beverages industry in the premises of the study area were observed during the reconnaissance survey. The studied HMs concentrations in study area were compared with the published geochemical results similar type from four different countries such as China, Italy, Pakistan, and Nigeria (Table 8). It was found that the concentrations of Cr, Pb, and Zn in the studied samples are lower while Cu, Ni, and Cd are higher than the open dumpsite situated in North China [52]. However, HMs concentrations around the dumpsite in Italy were found to be higher for all elements except Cd as compared to that of the studied samples [53]. The high rainfall (900 mm) could provide facilitation in leachate generation at the dumpsite of Italy as compared to the study area, with an average rainfall 420 mm. Similarly, HMs concentrations of Mahmood Booti dumping site, Lahore, Pakistan were found lower than that of the studied samples [54]. Although the size of Mahmood Booti dumpsite is larger than the study area, the Mahmood Booti dumpsite was authorized as a dumpsite under the city district government Lahore (CDGL) and Turk company Oz Pak. The dumpsite was operated under proper vigilance as compare to Chowa Gujar dumpsite. The concentrations of Pb, Cd, Ni, and Zn in the dumpsites in Islamabad, Pakistan were found to be higher than that of the studied samples [55], because 96% of total domestic, commercial, and industrial waste generated was dumped at this site. High volume of waste as compared to the study site could be the reason for high contamination concentration. In Nigeria, HMs concentration at dumpsite were found to be higher for Cd, Cr, Ni, and Pb while lower for Zn as compared to that of the studied samples [56]. Electronic waste was dumped at this site, while the waste dumped at Chowa Gujar was heterogeneous as confirmed by the geophysical survey discussed in Section 4.1. Therefore, waste type could be the reason for high HMs concentration at dumpsite of Nigeria. HMs The studied HMs concentrations in study area were compared with the published geochemical results similar type from four different countries such as China, Italy, Pakistan, and Nigeria (Table 8). It was found that the concentrations of Cr, Pb, and Zn in the studied samples are lower while Cu, Ni, and Cd are higher than the open dumpsite situated in North China [52]. However, HMs concentrations around the dumpsite in Italy were found to be higher for all elements except Cd as compared to that of the studied samples [53]. The high rainfall (900 mm) could provide facilitation in leachate generation at the dumpsite of Italy as compared to the study area, with an average rainfall 420 mm. Similarly, HMs concentrations of Mahmood Booti dumping site, Lahore, Pakistan were found lower than that of the studied samples [54]. Although the size of Mahmood Booti dumpsite is larger than the study area, the Mahmood Booti dumpsite was authorized as a dumpsite under the city district government Lahore (CDGL) and Turk company Oz Pak. The dumpsite was operated under proper vigilance as compare to Chowa Gujar dumpsite. The concentrations of Pb, Cd, Ni, and Zn in the dumpsites in Islamabad, Pakistan were found to be higher than that of the studied samples [55], because 96% of total domestic, commercial, and industrial waste generated was dumped at this site. High volume of waste as compared to the study site could be the reason for high contamination concentration. In Nigeria, HMs concentration at dumpsite were found to be higher for Cd, Cr, Ni, and Pb while lower for Zn as compared to that of the studied samples [56]. Electronic waste was dumped at this site, while the waste dumped at Chowa Gujar was heterogeneous as confirmed by the geophysical survey discussed in Section 4.1. Therefore, waste type could be the reason for high HMs concentration at dumpsite of Nigeria. HMs concentration of Chifing, Inner Mongolia China dumpsite were found to be higher for Pb and Zn, and lower for Cd and Cu compared to the concentration of HMs in study area [57]. The concentrations of Pb and Zn were high because the survey was conducted around the non-ferrous mining waste dump.

Conclusions
The Chowa Gujar dumpsite in the City of Peshawar is the model study site. This study characterized the MSW dumpsite by employing geophysical and geochemical methods. ERT models showed different stratigraphic layers and geometry of the waste dump in the sub-surface. High-resolution ERT sections provided accurate results by distinguishing the waste type, waste pockets, the thickness of waste contained zone, and sub-surface distribution of the associated leachate plume in the vicinity of the dumpsite. ERT results showed leachate contaminant plumes that are extended up to a depth range of 8 m and 14 m from the south and east edges of the dumpsite. GPR did not permit the deeper structure at the studied site to delineate the accurate vertical extension of the leachate plume or mapping migration pathways due to the strong attenuation of electromagnetic waves in the highly conductive zone. However, GPR data helped in demarcating lateral boundaries of different waste components and they have provided a positive correlation with ERT. Similarly, radar results have confirmed the isolated patches of very high conductive leachate-filled areas along a few profiles. Despite several factors (heterogeneous nature of waste material, moisture variations, lithological characteristics, and distance from the Bara River), the geophysical properties (electrical resistivity values 31-80 Ωm) of the municipal waste were measured. The commercial waste was characterized by very high resistivity values (80-127 Ωm). The interpretation of the non-destructive geophysical data sets was made more accurate by comparing the results from GPR and ERT techniques. Geochemical data revealed that agricultural soil of study area is contaminated. The measured concentrations of HMs and geo-statistical results revealed that the soil of study area is contaminated. It is evident from the geophysical and geochemical assessment that Chowa Gujar dumpsite may pose an obvious threat to the agricultural soil and both surface and subsurface water in the study area, which could further deteriorate due to the continued dumping of a large volume of municipal and commercial wastes. Therefore, precautionary measures should be taken to avoid further contamination and keep the soil of the area environmentally sustainable for agricultural use. This study could be helpful to facilitate the policymakers and stake holders to develop the optimal strategies to restore the polluted dumpsite into a pollution free sustainable environment.