Landfill leachate assessment by hydrological and geophysical data: case study NamSon, Hanoi, Vietnam

The main objective of this study was to assess the impact of the NamSon landfill on subsurface geological structure and hydrological environment by geophysical techniques and hydrochemical analysis of surface and groundwater. The electrical resistivity tomography (ERT), self-potential (SP) and very low frequency (VLF) methods were used for the investigation of geological structure near the landfill. Three profiles (900 m long in total) of the two-dimensional ERT, VLF density sections and 180 SP data points scattered throughout the study area near the disposal site constituted the basis of the data used in analysis. Additionally, surface water and groundwater samples were collected from six sites in the area for the chemical analysis. Interpretations of geophysical data show a low resistivity zone (< 15 Ωm), which appears to be a fully saturated zone with leachate from landfill. The results of the geophysical investigations are not always fully confirmed by the results of hydrochemical analysis. The quality of water in the vicinity of the landfill dramatically decreased over the year (2015–2016) and actions should be taken to inverse this negative trend.


Introduction
One of the major environmental problems caused by the municipal solid waste landfills is the leachate which is a consequence of precipitation, surface run-off, and infiltration or intrusion of groundwater percolating through a landfill, biochemical processes in wastes cells and the inherent water content of wastes themselves. Leachate is produced always when water or another liquid comes into contact with the wastes. The landfill leachate potentially results in contamination of the groundwater and surface water supplies and threatens human health when migrating from the landfill into neighbouring lands and waters. With the growth of population and development of the industry the problem of landfill leachate becomes more and more severe. Generally, the best method of controlling the pollution of environment by the landfill leachate is to treat leachate to remove the hazardous components before they enter the water systems. This is because, once the leachate enters the water bodies, they become very difficult and expensive to clean.
When analysing landfill leachate a number of complex factors should be taken into account, these are, inter alia, solid waste composition, age of the waste, operation of the landfill, hydrogeological conditions in vicinity of the landfill site, rate of the water movement through the waste, landfill temperature, moisture content, pH, landfill chemical and biological activities and seasonal weather variations [1][2][3][4]. The methods of characterisation and treatment of landfill leachate have been developing for last 40 years [5][6][7][8][9]. The most often used methods are based on two main catalogues of conventional and advanced technologies. The conventional technologies include physical, chemical and biological techniques of treatment. The advanced technologies are introduced to three integration groups of physico-chemical, multibiological and physico-chemical-biological processes [10][11][12][13]. Since it is difficult to obtain satisfactory quality using any of those methods by itself, a combination of physical, chemical and biological methods are employed for efficacious treatment of landfill leachate.
NamSon is the biggest landfill in Hanoi and receives more than 5000 tons of domestic waste per day. The landfill leachate is one of the most serious environmental problems in the vicinity of NamSon. During the last decade the landfill and its impact on environment have been the subject of the intensive studies and environmental projects. Hoang et al. [5] used the sequencing batch reactor technique for leachate treatment; Khanh [8], compared the landfill leachate treatment systems in Vietnam and abroad to choose the most suitable technique that would meet the expectations of Vietnamese standards for NamSon. Van et al. [11] published the results of the experimental investigations on using poly aluminium chlorite (PAC) and UV/O 3 that have been conducted for the determination of optimal pH value, reaction time and PAC concentration for the removal of chemical oxygen demand (COD) and colour for the NamSon landfill. Khai and Quynh Trang [12] analysed the recovery ability of N, P in leachate from the NamSon landfill for pre-treatment of ammonium. Nguyen et al. [14] studied the changes in the hydrogeochemical properties of groundwater during dry and rainy seasons in 2011 by analysing the major ions for Hanoi region. Recently, in 2012, a hydrogeological database was developed for Hanoi region, which was used Bui et al. [15,16] for characteristics of hydrogeological conditions in 240 boreholes covering the NamSon landfill area. Additionally, Thinh and Strokova [17] published the article presenting the land subsidence caused by groundwater exploitation in Hanoi by analysing complex data on geology, hydrology, soil properties, including SocSon district.
In this paper, to deeper understand the influence of the NamSon landfill activity on its closest vicinity, we propose the procedure of analysis of the leachate and assessment of contamination it creates within the context of variation of geological structure, hydrogeological and hydrological conditions nearby the landfill. The geophysical methods, which we apply, have proven to be an efficient and effective means for characterizing the subsurface geology and hydrology associated with landfills. A combination of various geophysical methods as well as geoelectric and electromagnetic techniques used in hydrology and hydrogeology was applied for the research in the north of the NamSon landfill. These techniques were applied in the lowest altitude terrain around the landfill to monitor the volatility of the soil and water environment that may be caused by the landfill operation. We also show the temporal deterioration of the water environment in the proximity of the landfill.
The geophysical characterization of landfills is now a routine practice. Even so, due to the inherent uniqueness of each site, it is necessary to select carefully the method and appropriately design the survey parameters to derive adequate benefit. In many cases, particularly when multiple objectives are to be met, it is best to use complimentary geophysical methods. A variety of techniques can be applied to define landfill properties, leachate circulation, or regional contamination from offsite seepage [18].
The article is built as follows: in the second chapter the study area is described in detail, next section deals with the geological features of the proximity of the NamSon landfill, the fourth chapter defines methodology of the field and laboratory research and describes the results of measurements and last two sections provide discussion of the results and present final conclusions.

Description of the study area
The NamSon landfill is located in two communes: NamSon and BacSon, SocSon district, Hanoi City. In 1999 (phase 1) the landfill was designed to occupy total land area of 83.5 ha with ten cells buried garbage and 2000 tons of garbage a day. The volume of wastes received steadily increases and currently exceeds 5000 tons/day. In 2014 the first part of the landfill has been filled up completely and now the NamSon phase 2 is being operating with the total area of 75 ha. Currently, solid wastes delivered to the NamSon landfill are not treated properly. While building, the bottom of the landfill was covered with a layer of geotextile which lays on soil and rocks beneath the broken stamping chapped. However, after a long time of use, then the geotextile layer has decomposed [8,12]. Therefore, due to the difference of pressures the water flows from the top of landfill (about 30 m) towards the floor and subterranean garbage and with the help of subtropical water seepage breaks through the bottom of the landfill aiming at the lower terrain. The fissured rock environment in the bottom of the landfill, with cracks and holes of size and shape different to define enables the water movement through the garbage to the environment. Burial is the most common treatment option of municipal solid waste in NamSon because this is a simple and low-cost technique. However, there are several hazardous components such as batteries, chemicals, heavy metals, etc., that, if not treated well, can be washed out by the water soaking through the landfilled wastes to the surface water, groundwater, causing serious environmental pollution and impact human life.
Since the year 2000 the NamSon landfill leachate has been the subject of the combination of froth flotation and biological treatment, but these processes are still not functioning effectively enough. According to the nearby LaiSon village inhabitants, the odours from the landfill are unbearable. In both the dry and rainy season the air the villagers inhale and water they use are highly contaminated due to high concentration of poisonous sediments in surface, groundwater and vapours in the proximity of the landfill. Before the landfill operation, no one from the villagers suffered from oncological diseases but nowadays, 40 people reported to have died in the past 3 years of cancer and about a dozen are diagnosed in terminal stage cancer.

Geological features of the study area
The topography of the case study terrain is fairly simple. The study area is generally hilly with a very soft transition from 10 to 20 m (see Fig. 1). The study area is located within the region of moist tropical monsoon heat. Average annual rainfall is 1670 mm and average air humidity is 84%. Contrary to the surface the geological and hydrogeological structure is quite complex (Fig. 2). The aquifers are of the holocene deposits (qh) and pleistocene deposits (qp) [16,19]. The water layer deposits of Triassic rift (t), exposed in the north of SocSon district, are characterized by slightly better water quality in the Quaternary sediments [18,20]. The silver-coloured soil groundwater samples sites ( + S1, + S3, + S5). c The map of elevation for NamSon landfill and the surroundings. Elevation in metres above sea level in some parts of the research area is the result of the ancient alluvial ferralitic products [21].
According to the Geology-Mineral map (scale 1:50,000) published in 2015 by the Hanoi Institute of Geology and Environment the structure of SocSon district reveals mainly in the upper Triassic system and the Jura conglomerate [22]. We can also distinguish two stratigraphic formations: the KhonLang (T2a kl) and the NaKhuat (T2 nk). Additionally, there are also two Quaternary stratigraphic units: The ThaiBinh formation (a albQ23, tb) and the Hanoi formation (a, dp, apQ12-3 hn) (Fig. 1).

Methodology and results of research
Our research was divided into two steps: the first one consists on using geophysical methods (supported by the 'classical' borehole measurements) to determine the geological structure, hydrogeological parameters and potential trajectory of leachate under and in the vicinity of the NamSon landfill; the second step requires the in-field and laboratory analysis of the basic physical, chemical and biological properties of surface and groundwater threatened by the landfill leachate. The investigation carried out in March 2015 was repeated in March 2016 and the differences were analysed.

Geophysical investigations
Electric and electromagnetic geophysical methods, such as ERT, SP, VLF methods have been widely used for investigating groundwater [6,24], geological [7,[25][26][27] and environmental features [13,[28][29][30][31]. Due to the complex morphology of the study area and the difficulty in access to the measurement place, the geoelectrical and electromagnetic methods were applied in this study for shallow geological mapping, particularly for the identification of fractured/saturated or conductive zones.

Self-potential (SP)
Recently the SP method has been adapted for hydrogeological and water engineering applications, with the use of highly sensitive sensors and the careful application of data correction processes. We adopted a simple SP survey which consists of a base electrode position and a roving electrode to determine potential differences along the profile line survey [32,33]. Self-potential method makes use of natural currents flowing in the ground that are generated by electro-chemical processes to locate shallow bodies of anomalous conductivity and water circulation. The measurement equipment is very simple and needs only electrodes, a wire and a precise mV meter. For our case study we used two non-polarizing electrodes constructed by the solution sursaturated in copper sulphate with copper electrode in the middle and we used mV Terrameter SAS 300C meter for collecting data [34]. Readings are taken with one electrode fixed at a base station and a second one, the mobile 'field' electrode moved along every sampled profile around the survey area. Reading stations are distributed at regular intervals by 5 m along linear profiles T1, T2 and T3 (Fig. 1b, c). Small potentials (a few mV) are produced by two electrolytic solutions of differing concentrations that are in direct contact and also by the flow of groundwater through the porous materials. The depth of investigation depends on the size of the mineralized body and the depth of the water table for a mineralization potential (generally shallow, less than 30 m). The depth of water table for study area is from 4 to 5 m in dry season [14,15]. The SP data interpretation is mainly qualitative (profile, map). The SP anomaly is positive (resurgence) where the hydrostatic pressure decreases (i.e., in the direction of the water flux) and SP anomaly is negative for the infiltration case. Other occurrences produce spontaneous potentials, which may be mapped to determine the information about the subsurface. Spontaneous potentials can be produced by differences in the mineralization of the ground, electro-chemical actions, geothermal activities and bioelectric generation of vegetation. Four different electrical potentials are recognized: 1. Electrokinetic (or streaming)-due to the flow of a fluid with certain electrical properties through a pipe or porous medium of different electrical properties [1,4]. 2. Liquid-junction (or diffusion)-caused by the displacement of ionic solutions of dissimilar concentrations. 3. Mineralization (or electrolytic contact)-produced at the contact of a conductor with another medium. 4. Shale-occurs when similar conductors have a solution of differing concentration [35].
The data collected in the NamSon landfill by the SP method were processed to eliminate noise exemplified by the low background voltages for recognition of different spontaneous-potential sources. After processing, the SP data are in the range of − 260 to + 30 mV for the three measured profiles. We aim to assess the relationship of the near-surface geological formation with environmental occurrences, and we use the processed SP data to map by contouring surficial voltages between base/reference electrode(s) and the mobile electrodes (Fig. 3a, b). On the map of equipotential based on SP data collected in March 2015 (Fig. 3a) one can distinguish three positive anomalies (marked in brown) and one negative anomaly (in dark blue). The values of the positive SP anomalies are within the range 10-25 mV and negative is equal to − 180 mV. After a year (Fig. 3b), however, the SP measurement values changed clearly by six negative SP anomalies with amplitude range from − 260 to − 180 to mV and three positive anomalies with amplitude of similar range to this indicated in Fig. 3a. These results may be biased by the groundwater flow which leads to the background voltage variation.

2D electrical resistivity tomography (ERT)
The electrical resistivity tomography (ERT) method belongs to the group of geoelectrical methods commonly used in geophysical engineering. It measures differences in the electrical resistivity of the ground [4,[36][37][38]. 2D measurements with the use of the ERT method are taken along profiles, in consecutive sequences [39]. The 2D model of the subsurface is more accurate when the resistivity changes in the vertical direction, as well as in the horizontal direction along the survey line [40]. We used SuperSting R1/IP equipment with 56 electrodes [41] for data collection along three profiles T1, T2 and T3 by Wenner array with 5-m electrode spacing (Fig. 1b, c). The measuring system used automatically broadens distance between electrodes which results in increase in the depth of the measurement. The electrodes are chosen from among the electrodes connected to the cable until there are no more programmed combinations left [42][43][44]. The NamSon measurement results are presented in 2D cross-sections showing electrical resistivity vs depth. Interpretation involves confronting given geoelectrical cross-section and its references with subsurface conditions using EarthImager 2D software [41]. The local picture has been ascertained from the electrical resistivity image profile T1(A), T2(A) and T3(A) which were analysed in March 2015 (Fig. 4).  environments, while it also provides an adequate resolution [35]. The RMS error for iterations of resistivity data varies by 12.78% for T1(A), 9.46% for T2(A) and 12.96% for T3(A) in the inverse model resistivity sections. The results suggest delineation of geological formations such as various subsurface lenses like clays and in situ saturated sand-clays up to the depth of several meters. The principal aquifers are formed in the unconsolidated alluvial formations, deposited under various sedimentary environments. The geoelectrical profiles (Fig. 4) show resistivity values that vary laterally with the depth; two levels can be distinguished in accordance with the obtained resistivity values.
After 1 year, also in March (dry season), we repeated the 2D electrical resistivity measurements by Wenner array with the same multi-electrode resistivity instrument and the measured resistivity profiles were interpreted with the use of the same software. Electrical resistivity values were obtained from three parallel lines (Fig. 5). The RMS error computed for iterations of resistivity data varies by 11.68% for T1(B), 9.11% for T2(B) and 9.57% for T3(B) in the inverse model resistivity sections. The results reveal the lenses at depths of about − 16 to − 20 m in the large are of the case study field with very low resistivity values less than 15 Ωm (dark blue) in section T1(B), T2(B) and T3(B). According to the 2D inverted resistivity sections, low values (< 20 Ωm) correspond to areas that may be occupied by contaminated groundwater from leachate. The high resistivity values (> 150 Ωm) in profile T1(B), T2(B) and T3(B) are associated with consolidated materials. Additionally, very high resistivity zone (200-290 Ωm) in profiles A and B are interpreted as consolidated formations, sandstone, siltstone, tuffite and clay shale (Fig. 2).

Very low frequency (VLF)
Very low frequency (VLF) is a geophysical method commonly used in engineering investigation and environmental studies [4,35,[45][46][47]. It is based on the induction of the secondary magnetic field Hs in the subsurface conductor of conductivity σ due to the effect of an artificially generated primary field Hp. This means that electric currents in the conducting body (e.g., a fracture) are generated when radio waves (EM field) pass through it, creating another magnetic field (Hs). If the source and receiver are placed near a more conductive zone, stronger eddy currents may be caused to circulate within this zone and an appreciable secondary magnetic field will thereby be created. In close neighborhood of the conductor the secondary or anomalous field may be compared in magnitude to the primary or regular field, in which case it can be detected by the receiver. The secondary field strength, Hs, is usually measured as a proportion of the primary field strength, Hp, at the receiver in percent. The inphase amplitude Hp is the real component, while the out of phase Hp is the imaginary component or quadrature component. The anomaly is defined as Hs/Hp. The ABEM WADI equipment used for VLF-EM measurements uses military transmitters as the source of primary electromagnetic waves Hp which is located several km away at the high powered military communication transmission stations. The transmitter's antenna transmits signals continuously at low radio frequency of 15-30 kHz.
The VLF measurements in NamSon were accomplished along profiles T1, T2 and T3 (Fig. 1) that ranged in length of 900 m by a spacing point of 10 m. The selected sites in the field were good criteria for VLF signal which could be acquired from the transmitters of 18-22 kHz frequency. We found it difficult to interpret three VLF curves, so we applied filtering techniques to enhance data and make tiltangle crossovers easier to identify by Karous-Hjelt filter [48]. The real and imaginary parts are presented in the plots by SECTOR software [49] and their interpretation can reveal different types of the VLF data that can be acquired. The real part always shows a positive peak above a conductor, while the imaginary part can show a positive or negative peaks, depending on the conditions of the overburden layer. An interpretation of in-phase (real part) anomaly of the VLF-EM data was carried out by the analysis of the transform of the in-phase component and the amplitude of its analytic signal. Further, the source and depth to the top of a subsurface conductive body was delineated. The amplitude of the analytical signal of the data was observed to mimic the conventional Fraser-filtered operation and was used to locate the exact location of the anomalous body. The inphase component of the transform yields an approximate depth of 20 m to the top of the conductor as shown in Fig. 6. The colour gradation indicates that the anomaly Hs/Hp in % is (the deeper the colour, the stronger the anomaly): red colour-positive Hs/Hp relationship; blue colour-negative Hs/Hp relationship; yellow/green colour-near zero Hs/Hp relationship. Negative Hs/Hp usually relates to the major resistivity layers, whereas positive anomaly indicates the major electricity conductivity layers (for instance, the presence of clay materials). Successive inflections of Hs/Hp indicate the fractured/saturated zones. The current density variation was assumed to indicate the presence of fractures, along the infiltration of leachate. After 1 year, we repeated the measurement of the 180 points by the [49], the results are presented in Fig. 6b (also Hs/Hp in %) but the picture shows the changes of VLF anomaly by far higher in the position and strength. The anomalous zone of the both times interpreted is located in the middle of the study area with direction from northwest to southeast (NW-SE). This direction indicates from upstream to downstream (direction of the canal LaiSon) and proves that the water infiltration through the landfill strongly influences the quality of surface water to groundwater to the depth of 20 m by fractured zone in the study area.

Hydrochemical and hydrophysical analysis
There are many sources of water in the landfill, for example, precipitation, groundwater infiltration, surface run-off, and the disposal of liquids or sludge. The amount of water entering the landfill depends on difficulty to evaluate factors such as the seasonal distribution of the water balance, topography, geological and hydrogeological properties of the site, soil type, land cover, human activity, etc. [50]. Water plays three major roles in land filling: as a medium essential for bacterial activity and for chemical reactions in organic matter, as a medium for movement of contaminants and as a solvent. Obviously, the landfill activities negatively influences the quality of groundwater and surface waters in the neighbourhood [10]. As water percolates through the landfill material, it carries the contaminants from the inorganic and organic waste components degraded by microbial activity to simpler easily soluble compounds. Obviously, percolating through the landfill, water deteriorates in quality, but also changes its colour and odour. The polluting potential of such water may be 10-100 times stronger than the one of the raw sewage. The main components of leachate are the ions of calcium, magnesium, potassium, iron, sodium, ammonium, bicarbonate, sulphate and chloride, trace metals such as manganese, zinc, copper, chromium, nickel, lead and cadmium, a wide variety of organic compounds which are usually measured as total organic carbon (TOC), chemical oxygen demand (COD) or biochemical oxygen demand (BOD). Of course, the individual compounds hazardous at very low concentrations should be considered, too, including pesticides, benzene, phenol, microbiological components [23,51,52].
In this research project, for analysis we sampled surface water in the LaiSon canal at three locations (N1, N3, N5) and groundwater in three boreholes drilled in the area of research (S1, S3, S5) (Fig. 1b) Tables 1 and 2, for surface water and groundwater, respectively.
The environmental indicators of water quality can be roughly divided into the physical indicators as temperature, electro-conductivity (EC), total dissolved solids (TDS) and total suspended solids (TSS), the chemical indicators as pH, biological oxygen demand (BOD), chemical oxygen demand (COD), dissolved oxygen (DO), nitrate and heavy metals and biological indicators as bacteria E. coli and coliform. Because direct measurements of water quality are expensive, we used on-site test kits and water samples analysed at laboratories by real-time monitoring.
The Table 1  As far as the groundwater is concerned (Table 2), the results of analysis show that in 2016 the values of most parameters reveal higher concentration than in 2015. As in case of the surface water, the groundwater is slightly alkaline with pH varying from 7.07 to 7.80. The concentration of the TDS in downstream (S3 and S5) is higher than upstream (S1). High values of the COD, Cadmium, Lead, Copper, and Coliform were also sensed in the samples. Furthermore, all the groundwater samples indicated Cyanide concentration exceeding permissible limits. High concentration of the TDS and EC was observed in all wells which mean that the groundwater in the landfill area is of poor quality.

Discussion
The main objectives of this paper are to present the integrated geophysical methodology of the landfill leachate investigation to assess the physical, chemical and biological temporal degradation of the near-surface geological structure and water bodies in the area of the NamSon landfill. We focused on the effects of the operation of the landfill imposed on the surface water environment, groundwater, sediment by analysis of the water in the formations near the ground, including the rocky land. To accomplish the abovementioned goals, we chose two points in time: March 2015 and March 2016 (dry season) to carry out the survey of water quality and observe the changes in environment around the landfill.

3
The methods chosen to survey the geological structure in the three profiles are based on the conductivity, resistivity, and dielectric constant of the geological environment.
When comparing the map of equipotential based on the self-potential (SP) data carried out on March 2015 (Fig. 3a) and the equivalent map for March 2016 (Fig. 3b) we can see very clearly the negative differences in the quantity, magnitude, as well as, overarching scope related to the location of the SP anomalies. The location of the anomaly peak SP in 2016 changed only slightly in comparison to the year 2015, but the amplitude of the anomalies changed significantly which points to the negative anomalies which appeared in large numbers nearby the landfill. This was confirmed by the results of the analysis of the quality of the surface-and groundwater.
The model of the electrical resistivity tomographic profiles T1(A), T2(A) and T3(A) in March 2015 (Fig. 4) and the profiles T1(B), T2(B) and T3(B) in March 2016 (Fig. 5) showed apparent difference in the resistivity values reaching as far as ca. 50 m deep below the ground surface. The investigation revealed that the geological structure is highly inhomogeneous vertically, but fairly homogeneous horizontally. As a result we can observe here good conditions for groundwater run-off, the formation of the lens containing water and different sediments.
The VLF-EM results mapped the shallow linear conductors that are probably the fractured/saturated zones of varying length in the area which is of significant hydrogeologic importance for groundwater bearing. Interpretations are based on the high amplitude signal, which suggests the presence of weathered or fractured zones. The analysis of the VLF data by in-phase component was used to delineate the source and depth to the top of a subsurface conductive body. The amplitude of the analytical signal of the data showed the location of the anomalous body. The anomaly observed in the VLF survey is also associated to the deviation of the measured signal from the regular level, and this deviation is a result of the response of the subsurface geological objects of interest. The integration of these results shows a main conductive anomalous zone in the northern part of the landfill site that up to 20-30 m deep. Because the measurements of geophysical survey carried out in dry season with the natural surface run-off lets us to draw a conclusion that the conductive anomalous body clearly confirms the presence of leachates. Both the VLF and SP data interpretation provided good response in fractured zones for water moving from upstream to downstream of LaiSon canal which is also reflected by the results of water analysis in 2015 and 2016. (See Fig. 6a, b the blue colour negative anomaly of VLF data create fractured zone by direction NW-SE and Fig. 3b from SP data by also negative anomaly).
In the sedimentary formation the sand, clay and a finegrained sand formation in the subsurface show minimum resistivity among all the aquifer materials. We observe also the high TDS, EC, Arsenic, Cadmium, Lead, Copper and Iron in S3, S5 and N3, N5 (Tables 1, 2). The higher concentrations may be the consequence of the infiltration of the landfill seepage water flowing to the aquifer and canal. The lower resistivity may be due to the flow of polluted water into the aquifer (see ERT version of three profiles for 2015 and 2016 to confirm of ERT interpretations).
The results of the hydrochemical and hydrophysical analysis showed that the values of pH for almost all samples ranged within 7.4-7.8, i.e., the water is slightly alkaline [compare 2], but still fits into the allowed safe range. On the other hand, the TDS values for surface and groundwater are very high. Although generally the TDS is not considered as a primary pollutant, it is used as an indication of aesthetic characteristics of domestic water and as an aggregated indicator of the presence of a number of chemical contaminants [10]. Primary sources of the TDS in waters are soil contaminants and polluted water from the landfill [3]. Bear in mind, however, that the concentration of some naturally occurring dissolved solids increases as a result of weathering and dissolution of rocks and soils which influences the results of analyses [53].
The analysis revealed that the concentration of the biological indicators, the E. coli and coliform bacteria, is very high especially in surface water. This phenomenon can be explained as follows: The canal always flows according to the slope in the elevation of the terrain, i.e., in this case from NW to SE. It flows through the study area which is not a residential area, but upstream there is a number of farms where the farmers keep livestock such as the pigs and cattle. As a result they discharge daily a certain quantity of livestock sewage to the water environment (even though they use water canal for drinking). They also deposit contaminations on the banks of canal that are flushed by the canal water or rain (in rainy season) which influences the high concentrations of Coliform and E. coli in samples of surface water. In the dry season, the concentration of feacal Coliform and E. coli downstream the canal is caused by the weak flow and accumulation.
Additionally, the eutrophication of the water in the vicinity of the landfill is one of the key environmental problems of the NamSon. Obviously, it results in the sudden growth of plants taking up the nutrients (nitrogen, phosphorus and others) from the water and soil [54]. The result of water analysis showed high values of P and N for surface water and groundwater (Tables 1, 2). However, as a result of the increased photosynthetic processes, more oxygen is produced as a byproduct, but in total the eutrophication is always a serious pollution problem. All in all, however, the values of total P and N were high in 2015, after 1 year they even increased.

Conclusions
An integrated, geophysical survey, physical and hydrochemical analysis methods were employed to assess the subsurface geologic formations, aquifer location, dynamics and intrusion of seepage water in the vicinity of the NamSon landfill. The case study area generally close to the ground surface is characterized by the non-homogenous formations which are largely responsible for the transfer of contaminated water into the underlying aquifer. Most of the ionic compositions present in the water effectively indicated the effect of the seepage through the landfill. In water we found particularly high concentration of the TDS, COD, EC, Arsenic, Cadmium, Lead, Copper and Iron. The values of the TDS and some heavy metals concentrations much exceed the standards due to in situ contamination of groundwater in buried soil and soft sediments.
The detected lowering of resistivity is due to the encroachment of very high conductivity of surface water drained by the aquifers and also infiltration of contaminated water from the landfill into the lens and fractured sediments. The equipotential map based on the SP data demonstrates the direction of the water.
The ERT surveys indicated that the subsurface geological structure created favourable conditions for the formation of the lens containing water as well as the tropical disruption by the blocks of sediments. Through such a structure, especially in dry season, the surface water from the LaiSon canal easily infiltrates to the aquifer, which is indicated by the temporal change of the resistivity values. In particular, the field research carried out in March 2016 reveals that the tropical high conductivity structure broadened over time, especially with regard to section T3(B) of the profile measurement located right next to the LaiSon canal. This can be explained by the fact that in the dry season the highly polluted water from the canal is always recharging the aquifer. Poor water quality in the canal is the direct effect of the seepage through the deposits from the landfill.
The results of our research clearly indicate the need of regular monitoring of groundwater and surface water quality in the vicinity of the Nam Son landfill. Also we propose to install the waste water treatment plant that would diminish the pollution of the seepage waters flowing directly to the canal LaiSon. The management of the landfill should also pay more attention to the environmental protection of groundwater resources from the domestic wells in the around landfill area.