Three Dimensional Modelling of Kharga Reservoir Water, New Valley- Egypt, Using Magnetotelluric Data

Egypt suffers a continuous decreasing in water resources, resulting from constructing many dams on the river Nile as in Sudan and Ethiopia. The river extends in Egypt for more than 1000 km and is a main source of water required for life and agricultural purposes. The Kharga Oasis lies between latitudes 24o00’ 26 o00’ N and longitudes 30o 00 31o 00 E. It represents a region of the new valley governorate which includes Siwa, Bahariya, Farafra, Dakhala and Baries Oases. It is considered as a depression (basin) and is characterized by surface springs, playa and sand dunes. The reservoir water extends for more than 125 km from north Kharga city to Baris Oasis in the south with a width of about 25 km, covering an area of 3125 km2. It lies more than 200 km far from the River Nile (Figure 1). As groundwater is the only source for agricultural areas, the evaluation of reservoir waters needs more information about its spatial extents and depths as well as about the sources of recharging.


Introduction
Egypt suffers a continuous decreasing in water resources, resulting from constructing many dams on the river Nile as in Sudan and Ethiopia. The river extends in Egypt for more than 1000 km and is a main source of water required for life and agricultural purposes. The Kharga Oasis lies between latitudes 24 o 00 ' -26 o 00 ' N and longitudes 30 º 00 -31º 00 E. It represents a region of the new valley governorate which includes Siwa, Bahariya, Farafra, Dakhala and Baries Oases. It is considered as a depression (basin) and is characterized by surface springs, playa and sand dunes. The reservoir water extends for more than 125 km from north Kharga city to Baris Oasis in the south with a width of about 25 km, covering an area of 3125 km 2 . It lies more than 200 km far from the River Nile ( Figure 1). As groundwater is the only source for agricultural areas, the evaluation of reservoir waters needs more information about its spatial extents and depths as well as about the sources of recharging.
It is well known that MT is sensitive for electrically conductive layers. In this study an attempt to evaluate the environmental and regional structures of the reservoir water at Kharga Oasis using 3D-MT and compare the results with previously constructed models based on geological studies supported by borehole logs.

Hydro-geology and Boreholes Information
The geology of Kharga Oasis has been reported by several authors [1,2]. It is subjected, as well, to comprehensive studies of hydro-geology, geochemistry and geophysics by [3][4][5][6]. According to these works, the Kharga reservoir is divided into three main basins as: Kharga, north Kharga and south Kharga as based on aerial photography and wells information. The structures at these three zones lead to conclude that; 1) A main N-S fault extension draw the longer axis of Kharga main

Abstract
Most of the Egyptian populations live along the two banks of the Nile River in order to access water for daily needs and other purposes. About 90% of Egyptian territories is a desert with little amount of rainfall. In these arid deserts, groundwater is the only possible source that would support cultivation and civil expansions. Several Oases including Kharga are distributed in the western desert of Egypt, in which the groundwater is the main source of irrigation and daily life. The New-Valley governorate is keen on reclaiming the western desert and developing new urban areas.
Previously available boreholes and geological information are used to preliminary evaluate the three dimensional (3D) subsurface structures including the reservoir water and its environment.
In this study we utilize magnetotelluric (MT) recordings to spot more light on the Nubian aquifer of the Kharga Oasis. A 3D-MT inversion is applied using "w3dinvmt" code. The resultant 3D-resistivity indicates a low resistive layer associated with a Quaternary aquifer and extends from the surface down to a depth more than 50 m. Low to moderate resistivity values are found to indicate a deep Cretaceous aquifer defined at depths from 250 m to 500 m. A resistive zone is found to exist between them that can be associated with solidified limestone and phosphate layers. Furthermore, a high resistivity value appears to belong to the basement complex of Precambrian rocks in the Oasis. A constructed 3D-model is well matching with major hydro-geological structure of the Kharga Oasis that has been inferred from previous works.   from 350-600 m 3 /day, which is suitable for irrigation and drinking purposes. The deeper wells have salinity range from 513 to 1444 ppm, while shallower wells and springs vary in salinity between 308 to 4078 ppm. The ground waters of deeper wells are costly expensive for drilling and energy consumption but recharge from the Nubian sandstone aquifers [1]. The saturated thicknesses of shallow and deep wells are different according to the water penetration rate [5].

Magnetotelluric Data
The MT method uses the Earth's natural occurring electromagnetic fields to investigate the electrical resistivity structure of the subsurface from shallow depths (tens of meters) to deep depths [8]. On the other hand, electrical resistivity of a lithology is dependent on the existence of fluids, clay minerals and porosity in the rocks. The MT data were collected using EDU-7e MT system [9] designed for broadband survey (0.001-1000 s), audio-magnetotelluric (AMT) and magnetotellric (MT). The MT system and acquisition in Kharga Oasis are illustrated in (Figure 4). The software for data processing is updated using (Procmt & Mapros) programs for processing [9]. The MT data is acquired during a field trip in two profiles in the northern part of the Kharga reservoir which consists of 13 stations after omitting noisy soundings ( Table 1). The station separations are different in distance (1-3 km), due to the hard topography and sand dunes. The time series are edited manually to remove sections contaminated by cultural noise. The MT data quality appeared then to be of a good quality ( Figure 5).
Preliminary images of electrical resistivity models ( Figure 6) using Bostick, Occam [10] and Rebocc codes [11] are estimated [12]. The results are compared with the boreholes information and subsurface structures. Both 1D-MT models (Bostick & Occam) show the variations of vertical resistivity with depth ( Figures 6 a and b), while the result of 2D-MT (Rebocc) is presented for several zones of distinctive resistivities     10 sec. The impedances are rotated to be consistent with the coordinate system which is oriented in cross-strike illustrating the orientation of (N-S) and (E-W).
The model grid dimensions are Nx=28, Ny=28 and NZ=23 (plus 7 air layers). The mesh is created with a vertical factor of 1.2 and with a top-layer thickness 100 m (Figure 7). The horizontal grid spacing is 1 km in the central part and increases in both side until total model domain spacing of 100 × 100 × 100 km. Data error is set for the inversion as 10 % root mean square of Zxy × Zyx.
The initial model is homogeneous half space where the resistivity is 100 ohm-m. The results of 3D inversion are obtained at 3 iterations. The normalized RMS misfit is 3.5 and took place after 120 hours. The 3D-resistivity horizontal slices of the reservoir model at different depths (50, 250, 500 and 1500 m) are shown in (Figure 8). Also, 3D-resistivity can describe in the vertical cross-sections slices with depths ( Figure 9).

Results & Conclusion The 3-D resistivity structures of the Kharga Oasis (Figures 8 & 9) are characterized as following
The conductive thin layer ( Figure 8 a and b) has resistivity (10-100 ohm-m) extends to more than 50 m and represent shallow aquifer of Quaternary sediments. It can be rechargeable from water of irrigation ( Figure 6c). In general, the resistivity increases with depth, except those of the near surface layers (low resistive) of Quaternary sediments and those with low to moderate resistivity of Cretaceous [12].

Three Dimensional modeling
The 3D-MT inversion modeling is applied using "wsinv3dmt" code provided by Siripunvarapurn [11]. The number of periods and responses should not exceed the parameter set in 'para.h' of the code. The number of periods are chosen 3 per decade, and 4 responses are supported to the off diagonal impedances (Zxy and Zyx) respectively. The 3D-MT forward and inverse modeling has been reported by many authors since 90's [13][14][15][16][17][18][19][20][21][22]. It has also been tested and developed with different data sets from land and marine environments [23][24][25][26][27].
The "wsinv3dmt" code is used not only for gridded data but also for data acquired along profiles [11,28]. They interpreted 2D-MT profile data using a 3D approach. The profiles can be put together or separately using a code for single processor machine. A 3D minimum structure inversion is based on a data-space variant of the Occam scheme [10,28,29].
Generally speaking, a 3D inversion is associated with long time consumption and requires a large accessory memory. All computations are performed in this study using a UNIX system with 8 GB RAMS. The "wsinv3dmt" code forms the full Jacobian of the model parameter data mapping. The 13 periods selected for the inversion are from 0.01 sec. to  which contains some salts and fertilizers. Also, it can be interpreted as a zone of clay, plya and shale. This kind of aquifer can be used in agriculture only. Also, mixed limestone and phosphate layer (Paleocene rocks) at depth 50-250 serve as barrier between fresh and salinity water. Also, it is play a role of the cap-rock of the reservoir water system in the region.
-There are conductive layers of deep aquifer (Figure 8 b and c) belongs to cretaceous rocks and has thick layer of sediments (750 m) which starts 250 m and reachs to basement rock. The salinity range (300-1000 ppm) which is suitable for irrigation and daily life. The lithology consists of siltstone, sandstone with shale. The source of recharge is Nubian Sandstone Aquifer (NSA) which is the largest water aquifer systems (5 million cubic feet) of groundwater located in the eastern portion of the Sahara Desert [5].
-The conductive zones at depths 1500-2000m (Figures 9a, b, c and d) have resistivity between (50-100 ohm-m), could be the source of surface natural springs and geothermal activity in the area.
-The high resistive layers (blue color) are associate with and basement complex of Pre-Cambrian rocks. The 3D resistivity model is reliable than previous 1D and 2D models and adds more detail information to the reservoir water of Kharga Oasis. The extension and change of resistivity with different depths are manifested.

Recommendation
It is recommended that to carry out more MT stations to cover the boundary of water reservoir and using different geophysical tools like seismic velocity analysis, gravity and magnetic methods.