Sapele, located in Delta State (Figure 1), Western Niger Delta, Nigeria, lies between longitude 5˚E - 5˚45'E and latitude 5˚30'N - 6˚N in geographic coordinate. The city hosts a flow station and an oil rig owned by one of the oil and gas companies. Figure 1 shows the map of Delta State while Figure 2 shows the picture of a fared gas taken during the survey.

Oil and gas production is usually accompanied by substantial discharge of wastewater in the form of brines. Constituents of brines include sodium, calcium, ammonia, boron, trace metals, and high total dissolved solids (TDS). Oil spillage is a result of leakages of hydrocarbon from the pipes, and to an extent, poor maintenance of oil pipelines and poor monitoring of pressure regimes of the fluids with respect to the strength of the pipe, equipment failure, operators error, corrosion, pigging operations, flow line replacement, flow station upgrade and the activities of vandals in search of crude and refined petroleum products [9] .

As a consequent of the environmental problems highlighted above, 2-dimensional electrical resistivity imaging was carried out in a residential area close to the flow station to see if contaminants have actually infiltrated into the groundwater resources in the study area.

The Niger Delta is a large curve shaped delta which is located in Southern Nigeria like some other deltaic environments in the world. It occupies an area lying between longitude 4˚E - 9˚E and latitude 4˚N - 6˚N. It is bounded in the west by the Calabar flank, in the north by the Anambra platform and in the south by the Atlantic Ocean under which it extends Figure 1 and Figure 3. Both marine and mixed continental depositional environment characterize the Niger Delta of Nigeria [10] . The Niger Delta covers an area of about 75,000 km² (28,957 mi²) in southern Nigeria, where the Niger Delta discharges its water into the Atlantic Ocean through a series of distributaries.

From the Eocene to the present, the Delta has prograded Southwest ward, forming depobelts that represent the most active portion of the Delta at each stage of its development [11] . These depobelts form one of the largest regressive deltas in the world with an area of some 300,000 km² [12] a sediment volume of 500,000 km³ [13] and a sediment thickness of over 10 km in the basin depocenter.

The Niger Delta consists of three main tertiary stratigraphic units overlain by Quaternary deposit [14] shown in Table 1. These three subsurface stratigraphic units in the Niger Delta are Benin, Agbada and Akata formations. The base is the Akata formation comprising mainly of marine shale and sand beds consisting of dark grey sandy, silty shale with plant remains at the top. It is over 4000 ft thick. The underlying Agbada formation is a sequence of sandstones and shales [15] . It consists of an upper predominantly sandy unit with minor shale intercalations and a lower shale unit which is thicker than the upper sandy unit. It is over 10,000 ft thick. The aquifer in the Benin formation is largely phreatic. These formations are overlain by various types of Quaternary deposits [15] . These areas are made up of top soil, red laterite, clay, fine sand, medium sand and coarse sand in form of

pebbles. The thickness is variable but generally exceeds 6000 ft [16] .

The deposits of the Freshwater Swamps and the Sombreiro-Warri Deltaic Plain are universally considered to be recent expressions of and a continuation of the Benin Formation. They result from the sediment laden discharges of the River Niger that is spread on the delta by its various tributaries. The sediment is generally an admixture of medium to coarse-grained sands, sandy clays, silts and clays that eventually settle in fluvial/tidal channel, tidal flat and mangrove swamp environments, a process that has been ongoing since the late Quaternary and is related to interglacial marine transgressions [17] -[19] . The Niger Delta is one of the most hydrocarbon-rich regions in the world. Exploration and exploitation of hydrocarbons has been going on in the region since 1956, when oil was discovered there [20] . The oil and gas production and a rapidly growing population have resulted in environmental degradation of the Delta [21] .

The study area consists of Fresh water swamp, Coastal Plain Sands, Mangrove swamps, and Sombreiro-Warri plains [22] . The water table in the study area is approximately 4 to 5 m beneath the surface, and the direction of flow is towards River Ethiope which drains into the Atlantic ocean through the Benin River (Figure 1).

2-dimensional model, where the resistivity changes in the vertical direction, as well as in the horizontal direction along the survey line was employed in the study area to image the subsurface. In this case, it is assumed that resistivity does not change in the direction that is perpendicular to the survey line.

Geophysical methods are frequently used to study subsoil contamination caused by industrial residues of different nature. The effectiveness of electrical methods for the characterization of oil contaminated subsoil has been reported by several researchers [23] -[26] .

Geoelectric method (Resistivity method) has been used by several researchers to study environmental problems and groundwater studies. Reference [27] used Geoelectric method in the evaluation of Olushosun landfill site Southwest Nigeria and its implication on Groundwater. Geophysical and hydrochemical assessment of groundwater pollution due to a dumpsite in Lagos State was also conducted by [28] .

References [9] [29] also found out that high hydrocarbon content of soils has been known to affect soil physiochemical properties, which in turn affect the agricultural potentials of such soils. Hence, this work attempts to find out if hydrocarbon has infiltrated into the groundwater in the study area.

Hydrocarbon pollution is very complex. Fresh spilled hydrocarbon gives different result from aged spills. The following results from some researchers give clue on what to expect from the investigation in the study area.

Reference [30] in their work on “Two-dimensional electrical imaging for detection of hydrocarbon contaminants” pointed out the usefulness and the pitfalls of electrical tomography in the characterization of underground leakage of hydrocarbons. Experimental evidence, obtained from a joint geochemical and geophysical investigation approach, indicated that subsoil which has been saturated with diesel oil for a long period (>20 years) exhibits an increased conductivity. It suggests that electrical tomography could be useful for monitoring the effects of induced biodegradation (bioremediation) through the repetition of the survey at different times, in order to observe any changes in the resistivity due to the increase of free ions resulting from hydrocarbon degradation.

Electrical resistivity models for oil contamination in some area showed high resistivity [31] and low resistivity [32] [33] . Recent oil pollution shows a high resistivity anomaly, while mature oil pollution produces a low resistivity anomaly [32] . Months after spill, a low resistivity anomaly is developed in the contaminated zone, with a strength that depends on the geological characteristics of the subsoil [32] [34] . The formation processes of such low resistivity anomaly are related to chemical reactions and to variations in the physical characteristics of the oil contaminated zone. The low resistivity anomaly is caused by an increase in the total dissolved solids (TDS), due to bacterial degradation of hydrocarbons in the lower part of the vadose zone. Reference [32] found that aged contamination appears as a low resistivity horizon slightly above groundwater table (GWT).

Traverse 1 reveals the resistivity image along the profile. Figure 7(a) and Figure 7(b) show the Gradient and Wenner array maps respectively.

Integration of the two 2-D ERT array results show that subsurface is characterized with soil material with resistivity ranging from 42 - 15,000 Ω∙m reflective of varying degree of conductivity associated with varying lithology and fluid type. Figure 5 and Figure 6 are the borehole data used for comparing the interpretation of the 2-D ERT section.

The ERT result shows that the subsurface is composed predominantly of sand material from the surface to a depth of about 60 m beneath the surface. The 2-D electrical resistivity imaging result (Figure 7(A)) show the

resistivity distribution over a lateral distance of 126 m from the surface to a depth of about 21 m beneath the surface. The 2-D section revealed an anomalously high resistivity (1500 - 15,000 Ω∙m) structure within a lateral distance of 2 - 32 m; 36 - 40 m; 40 - 52 m; 58 - 63 m; 76 - 86 m; and 114 - 118 m at a depth of 0.5 - 9 m; 0.1 - 1 m; 5 - 18 m; 7 - 11 m; 8 - 10 m; and 0 - 2 m respectively. These anomalously high resistivity structures are attributed to the presence of hydrocarbon within the subsurface which may be due to leakages from various pipelines within the study area or possible activity of the vandals.

Correlation with the borehole data (Figure 5 and Figure 6) shows that the subsurface is characterized predominantly with sand, and may compose of clayey sand/clay in some locations. The 2-D ERT result shows that the resistivity within the depth of investigation ranges from 45 - 15,000 Ohm-m. Since the formation in the study area is compose of sand with water table between 4 - 5 m beneath the surface, the high resistivity (≥1500 Ω∙m) is a possible indication that the aquifer within this area may have been polluted by hydrocarbon.

Figure 8(a) and Figure 8(b) show the Gradient and Wenner array 2D resistivity structures. The two resistivity structures also show that the subsurface is equally compose of varying degrees of resistivity as can be seen from the resistivity values (43 - 15,000 Ω∙m), revealing varying degree of conductivity associated with lithology and fluid type.

The 2-D section also revealed an anomalously high resistivity (2000 - 15,000 Ω∙m) structure within a varying lateral and vertical location distributed over the entire traverse. At a lateral distance of 2 - 58 m, possible hydrocarbon pollution were noticed at a depth of 0.5 - 6 m, while at a lateral distance of 60 - 126 m the depth of pollution vary from 0.5 - 21 m beneath the surface.

Integration of the two array maps and correlation with the borehole data (Figure 5 and Figure 6) show that the resistivity of the sand layer within the depth of investigation of the 2-D ERT ranges from 45 - 15,000 Ω∙m and may compose of fine to coarse sand/clay in some location.

Traverse three reveal the resistivity image along the survey line. Figure 9(a) and Figure 9(b) show the Gradient and Wenner array structures respectively.

Integration of the 2-D ERT array results show that the subsurface has soil materials with resistivity ranging from 42.5 - 15,000 Ω∙m indicative of changing degree of resistivity associated with varying lithology and fluid type. The 2-D section revealed an anomalously high resistivity (1500 - 15,000 Ω∙m) structure within a lateral distance 2 - 16 m within a depth of 0.1 - 6 m reflective of near surface hydrocarbon pollution. Similar occurrence was noticed within a lateral distance of 62 - 126 at a depth of 0.1 to 4 m; however within a lateral distance of 20 - 70 m and within 86 - 100 m the hydrocarbon pollution seems to have migrated to a depth of about 21 m beneath the surface. This may be a clear evidence of groundwater pollution by the hydrocarbon since the water table lies within 4 - 5 m within the study area.