JWARPJournal of Water Resource and Protection1945-3094Scientific Research Publishing10.4236/jwarp.2014.64025JWARP-43934ArticlesEarth&Environmental Sciences Numerical Groundwater Flow Modeling of the Coastal Plain Sand Aquifer, Akwa Ibom State, SE Nigeria niekanEdet1*RamadanAbdelaziz2BroderMerkel2ChieduOkereke3ThereseNganje3Department of Geology, University of Calabar, Calabar, Nigeria1Department of Geology, University of Calabar, Calabar, Nigeria 2Departement for Geology, TU Bergakademie Freiberg, Freiberg, GermanyDepartement for Geology, TU Bergakademie Freiberg, Freiberg, Germany* E-mail:aniekanedet@yahoo.com, aniekanedet25@gmail.com(NE);2003201406041932011 December 20133 January 2014 2 February 2014© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

Akwa Ibom State is located in the Niger delta region of Nigeria. The coastal plain sand aquifer, which underlies the area, is the major source of potable water for domestic, agricultural, and industrial use. The increase in demand for water has led to an increased exploitation of the aquifer in the last few years. Hence it is essential to quantify the amount of exploitable water in the aquifer. To achieve this goal, a regional numerical groundwater flow model using MODFLOW was calibrated under steady-state conditions to determine the aquifer’s hydraulic conductivity (K) and the recharge characteristics of the area. The hydraulic conductivity and recharge values from the pilot points calibration ranged from 1 to 75 m/d and 5.89 × 10-6 m3/d to1.23 × 10-4 m3/d. The calculations of the model showed that the average recharge rate amounts 0.40l/s/km2 or 12.60 mm/year. The estimated recharge suggests a high groundwater potential for the area. This is the first time that such a study has been carried out in the area, and it will serve as an important basis for future groundwater management and simulation of transient groundwater flow modeling in Akwa Ibom State.

Akwa Ibom; Coastal Aquifer; Modeling; MODFLOW1. Introduction

Akwa Ibom State was founded 26 years ago. According to the census figures, the population was 2,409,613 inhabitants in 1991, while in 2006, 3,902,051 people lived in the area. This growing population alongside expansion of urban, agricultural and industrial activities has led to an increase in the demand for water in the area. The demand is partly achieved by the indiscriminate sinking of boreholes by individuals, communities, private organizations, government at local, state and federal levels and donor agencies including UNICEF, EU, ADB, the World Bank, etc. All these activities result in an increase in groundwater abstraction from the main hydrogeological unit, the Coastal Plain Sand aquifer.

What worsens the situation is the lack of a law enabling the control of the abstraction of groundwater in the area. Hence, there is a need for an effective management document to be in place for conserving and controlling the use of the available groundwater resource. This study therefore is the first attempt to numerically model groundwater flow under steady-state conditions in Akwa Ibom State. The main objective is to develop a regionalscale steady-state groundwater flow model for the area. Minor objectives include the determination of hydraulic conductivity and recharge characteristics using available hydrogeological information and applying inverse modeling technique based on the pilot point and zonation techniques.

2. Study Area2.1. Geographic Setting

Akwa Ibom State shown in Figure 1 [1] , is the one of the coastal states in the Niger delta region of Nigeria covering an area of about 7081 km2, extending between latitudes 4321 - 5331 north and longitudes 7251 - 8251 east. It is about 150 km long and its width varies between 60 km in the north and 100 km in the south with an average elevation of about 60 masl. The area is bounded on the south by the Atlantic Ocean, on the southwest by Imo River and southeast by Cross River. On the northwest, north and north east the area is bounded by Cretaceous to Tertiary sediments. Precipitation is high in the south (>3000 mm/a) and decreases northwards (<2000 mm/a) with an average of 2500 mm/a [2] . Air temperature (26˚C - 28˚C) is uniform throughout the year [2] [3] . During dry season, the water level in rivers and streams decreases suggesting that part of the river is recharging the aquifer besides precipitation. Pan evaporation rates vary between 4.6 and 5.0 mm/d [2] . The Akwa Ibom area lies within the region of Nigeria, which is characterized by equatorial climate and tropical rainforest [4] .

2.2. Geologic Setting

The geology of the Akwa Ibom area has been extensively studied by [5] -[9] . As presented in these publications, the area is underlain by late Cretaceous to Quaternary sediments. The lowermost sediment is lateritic sandstone with minor shale of the Maastrichtian Nsukka Formation. Throughout the Paleocene to Early Eocene, a thick shale intercalated with sandstone and limestone was deposited (Imo Shale Formation), while in the middle of the Eocene semi-consolidated sandstone, siltstone and minor shale were deposited (Bende Ameki Formation). Since the Miocene-Oligocene grit and sand with intercalations of clay and lignite beds were deposited constituting the Ogwashi-Asaba Formation. Overlying this formation is a thick sequence of gravel, sand, silt, and clay constituting the Benin Formation also known as the Coastal Plain Sands and alluvial ridges.

2.3. Hydrogeologic Setting

The area has three hydrostratigraphic units including the Alluvial Deposit, Coastal Plain Sand and Lower Sand Aquifers [10] [11] . The Coastal Plain Sand aquifer is the most extensive and is separated from the Lower Sand aquifer by the Imo Shale aquiclude. The aquifer is composed of sand, gravel, and sandstone with clay intercalations. The thickness of the aquifer varies between 35 m in the north and 200 m in the south. The aquifer is considered to be semi-confined or unconfined, especially in the southern parts of the area. The main groundwater flow direction is from north to south into the Atlantic Ocean. There are also variations in the northeast-southwest and northwest-east directions into the Imo River and the Cross River.

3. Groundwater Flow Model

Based on the knowledge of hydrogeological system, a steady state mathematical model of the coastal plain sand aquifer was developed. The area of study is 7081 km2, therefore the model has a regional character. The hydrogeological conditions were defined from more than 200 boreholes drill log. Borehole data were obtained from [1] .

ReferencesEmmanuel, C.C. (2013) Monitoring the Quality of Groundwater and Leachate Contamination near Waste Dump in Coastal Plain Sand Aqufer. M.Sc. Thesis, University of Calabar, Calabar.Tahal Consultants (1980) Rural Water Supply in Mainland (Cross River State) Master Plan. Unpublished.Petters, S.W. (1988) Akwa Ibom State: Soil and Landuse. Unpublished, 23-24.Iloeje, N.P. (1991) A New Geography of Nigeria. Longman Ltd., Benin.Reyment, R.A. (1965) Aspects of the Geology of Nigeria. University of Ibadan Press, Ibadan.Short K.C. and Stauble, A.J. ,et al. (1967)Outline of the Geology of Niger Delta American Association of Petroleum Geologist Bulletin 51, 661-779.Murat, R.C. (1970) Stratigraphy and Paleogeography of Cretaceous and Lower Tertiary in the Southern Nigeria. In: Dessauvagie, T.F.G. and Whiteman, A.J., Eds., African Geology, University of Ibadan Press, Ibadan, 251-266.Kogbe, C.A. (1989) Geology of Nigeria. Rock View (Nigeria) Limited, Jos, 325-334.Wright, J.B. (1989) Review of the Origin and Evolution of the Benue Trough in Nigeria. In: Kogbe, C.A., Eds., Geology of Nigeria, Rock View (Nigeria) Limited, Jos, 359-376.Edet, A.E. (1993) Groundwater Quality Assessment in Parts of Eastern Niger Delta. Environmental Geology, 22, 41-46. http://dx.doi.org/10.1007/BF00775283Esu E.O., Okereke C.S. and Edet, A.E. ,et al. (1999)A Regional Hydrostratigrahic Study of Akwa Ibom State Global Journal of Pure and Applied Sciences 5, 89-96.McDonald, M.G. and Harbaugh, A.W. (1988) A Modular Three-Dimensional Finite Difference Groundwater Flow Model. US GeolSurv Open File Rep 83-875, chap A1.Harbaugh, A.W. (2005) MODFLOW-2005. The US Geological Survey Modular Groundwater Model—The Groundwater Flow Process. US Geological Survey Techniques and Model Methods 6-A16.EMRL (2004) Groundwater Modeling System (GMS) version 5.0 Help File, Environmental Modeling Research Laboratory, Brigham Young University, Provo.Doherty, J. (1994) PEST-Model-Independent Parameter Estimation. Watermark Computing, Corinda.Doherty, J. (2003) Groundwater Model Calibration Using Pilot Points and Regularization. Groundwater, 41, 170-177.http://dx.doi.org/10.1111/j.1745-6584.2003.tb02580.xAbdelaziz, R. and Bakr, M.I. (2012) Inverse Modeling of Groundwater Flow of Delta Wadi El-Arish. Journal of Water Resource and Protection, 4, 432-438. http://dx.doi.org/10.4236/jwarp.2012.47050Ebraheem, A.W., Riad, S., Wycisk, P. and Sefelnasr, A.M. (2004) A Local Scale Groundwater Flow Model for Groundwater Resources Management in Dakhla Oasis, SW Egypt. Hydrogeology Journal, 12, 714-722.http://dx.doi.org/10.1007/s10040-004-0359-8Varni, R.M. and Usunoff, E.J. (1999) Simulation of Regional-Scale Groundwater Flow in the Azul River Basin, Buenos Aires Province, Argentina. Hydrogeology Journal, 7, 180-187.http://dx.doi.org/10.1007/s100400050190Yidana, S.M. and Chegbleh, L.P. (2013) The Hydraulic Conductivity Field and Groundwater Flow in the Unconfined Aquifer System of Keta Strip, Ghana. Journal of African Earth Sciences, 86, 45-52.http://dx.doi.org/10.1016/j.jafrearsci.2013.06.009Carrera, J., Mousavi, S., Usunoff, E., Sanchez-Vila, X. and Galarza, G. (1993) A Discussion on Validation of Hydrogeological Models. Reliability Engineering & System Safety, 42, 201-216.http://dx.doi.org/10.1016/0951-8320(93)90089-HSala, J.M., Kruse, E. and Aguglino, R. (1987) Investigaciónhidrológica de la Cuenca del Arroyo Azul, Provincia de Buenos Aires. CIC, informe 37, 235 p.