Elsebaiya area has long attracted the attention of many geologists because of the presence of economic phosphorite deposits. The phosphate rocks are part of the Vast Upper Cretaceous marine phosphorite basin which is used as a source of phosphorus for fertilizers, and for manufacture of phosphoric acid and other chemicals .
The Egyptian phosphorite bearing rocks are geographically widely distributed as thin shallow marine deposits  .
The phosphate occurrences in Egypt, from North to South, may be subdivided into three East-West trending facies belts, Northern facies belt, Central facies belt and Southern facies belt . The phosphorite of the Central facies belt represents the most economic occurrences and was confined to the following localities, the Red Sea coast from Safaga to Quseir, the Nile Valley between Idfu and Qena and the Western Desert on Abu Tartour Plateau.
In this study, both ordinary kriging and cokriging methods were applied to investigate their precision on mapping the distribution of P2O5 % and thickness. So it may be possible to gain a clear picture of the distribution at the studied areas.
Statistical analysis starts with constructing of histograms and calculates their parameters such as: minimum, maximum, range, mean, standard deviation, variance, skewness, kurtosis, and coefficient of variations to illustrate the distribution of P2O5 % and thickness within the deposit while geostatistical analysis is used for modeling them in both sides of the River Nile in Elsebaiya area, Nile Valley, Egypt by using GS+ software.
Experimental variogram and cross variogram were constructed for two variables to characterize the spatial variability of the measured P2O5 % and thickness values depending on the available data. These data are X and Y coordinates, P2O5 % and ore thickness for each sample location.
Depending on the parameters of the variogram and cross variogram models of two variables, a comparison between ordinary kriging and cokriging techniques was used to present distribution model of P2O5 % which could be used in mine planning and production.
2. Description of the Studied Areas
Phosphate deposits are distributed in several locations in the Nile Valley. The studied areas are Um Tondoba mine in Elsebaiya East area and Western River Nile mine in Elsebaiya West area as shown in Figure 1. Um Tondoba mine is located between latitudes 25˚06' - 25˚15'N and longitudes 32˚44' - 33˚06'E while Western River Nile mine is located between latitudes 25˚05' - 25˚14'N and longitudes 31˚30' - 32˚10'E.
Ignoring the geology of the ore deposits is the most common source of errors when evaluating and producing them. Geostatistical evaluation of any ore body should be built on the basis of studying its geologic setting . It is possible to create a computer model from borehole data, but this model is meaningless if it is not compatible with the geologist’s understanding. The studied area is characterized by the presence of Upper Cretaceous - Paleogene succession which is
Figure 1. Location map of the studied areas.
considered as one of the most important geologic successions that have long been subjected to detailed geological studies either in Nile Valley or generally in Egypt.
In the studied areas, Duwi Formation consists of differed successive layers for the two areas as shown in Figure 2. The phosphate of this formation can be classified into two horizons where the bottom one represents the high economic phosphorite bed in Elsebaiya East area where it has largest thickness and high phosphorus content while in Elsebaiya West area, phosphate located as economic in two members (lower and upper from base to top)  .
A variogram is a graph illustrating the average variability between samples and the distance between them. The variogram can be described as the variation in values among samples some distance apart as a measure of their spatial correlation .
The characterization of phosphate ore deposits has been carried out through variograms.
The variogram function summaries the spatial continuity for all possible pairings of data for all lag distances h as:
Figure 2. Lithostratigraphic succession in the studied area.
where g(x) represents the value of any variable at location x and let g(x + h) represent the value of the same variable at h distance. The variogram is a function describing half of the expected squared differences between g(x) and g(x + h).
For cokriging analysis, the cross variogram should be determined by the following equation .
where u and v are the primary and covariate variables, respectively.
The variogram as a tool for structural analysis of the ore body reflects many parameters to represent characteristics of the studied regionalized variables, these parameters are: sill (C), nugget effect (C0) and range of influence (a).as shown in Figure 3.
3.2.1. Ordinary Kriging Method
Kriging is a geostatistical technique that is used to interpolate a surface from a scattered set of known points (X, Y plots) in which a continuous surface of values can be predicted between the known locations .
Ordinary kriging, which allows the mean of the measurements to vary spatially, was used in this study. The different parameters obtained from variograms
Figure 3. Variogram parameters.
were used to calculated sample property values based on weighting values assigned to sample values in nearby locations. The Ordinary kriging prediction Z(x0), is defined as:
With the constraint that the weights, Ordinary kriging, sum to 1 to ensure an unbiased prediction:
where Z(x0) is the estimated value at location x, Z(xa) is the known value at location xa, λa is the weight associated with the data.
3.2.2. Cokriging Method
Cokriging is an interpolation technique that allows one to better estimate map values if the distribution of a secondary variate is known. The secondary variate (also called a covariate or Z2) is sampled from the same locations as the primary variate Z, and also from a number of additional locations. If the primary variate is difficult or expensive to measure, then cokriging can greatly improve interpolation estimates without having to more intensely sample the primary variate . The kriging estimate is based not only on distance to nearby sample locations for Z and the variogram for Z, but also distance to nearby sample locations for Z2, the variogram for Z2, and the cross variogram for Z1 and Z2. This can provide a more robust estimate of Z at unsampled locations if Z1 and Z2 are sufficiently correlated.
4. Available Data
The present study was carried out at Phosphate Deposits on Both Sides of the River Nile in Elsebaiya Area, Egypt. The available data for the present study is based on sampling information derived from the location as shown in Figure 4 and Figure 5. The data used in this studied areas were obtained from Elnasr Mining Company. These data include technical report about mining works and site surveying map for the studied areas. Both ordinary kriging and cokriging techniques were used to achieve a goal of this study.
5. Conventional Statistical Analysis
The statistical analyses have been performed on two variables, P2O5 % and ore thickness, which gives the distribution of them and the standard parameters of statistics. Table 1 gives the summary statistics of the data sets for two variables at the studied areas. Figure 6 and Figure 7 shows the histograms of variables
Figure 4. Scatter plot of bore holes at east Sebaiya.
Figure 5. Scatter plot of bore holes at west Sebaiya.
Figure 6. Histograms of P2O5 % in the studied areas where, (a) Lower member of west sebaiya, (b) Upper member of west Sebaiya and (c) Um Tondoba mine in east Sebaiya.
Figure 7. Histograms of ore thickness in the studied areas where, (a) Lower member of west sebaiya, (b) upper member of west Sebaiya and (c) Um Tondoba mine in east Sebaiya.
Table 1. Statistical parameters of two variables in the studied areas.
that constructed for studied areas. The results show that the P2O5 % distribution is negatively skewed of all studied areas otherwise the distribution is positively skewed of the thickness at all studied areas. On the other hand, coefficient of variation values is the lowest in East Sebaiya for P2O5 % and thickness which shows a significant change in variables within the ore body. The higher range of P2O5 % occurs in the lower member of West Sebaiya, while the higher range for thickness occurs in the upper member.
6. Construction of Variograms
Depending on the available data for each variable, isotropic variograms and crossvariograms were constructed and fitted to a suitable model as shown in Figures 8-10. The variogram parameters for P2O5 % and thickness in the studied areas were summarized in Table.
As shown from the variogram models in Figures 8-10 and geostatistical parameters in Table 2 it is clear that P2O5 % distribution is characterized by the presence of nugget effect in the Lower Member of West Sebaiya and this is expected to the presence of intercalated within the deposits while in the Upper Member of West Sebaiya and Um Tondoba Mine in East Sebaiya, it is clear that high continuity in P2O5 % distribution while in thickness distribution the higher continuity occurs at the Western part of Sebaiya than the Eastern side. The cross variogram obtained strong spatial dependence of variables properties at all studied areas.
The screen effect ratio can be regarded as a criterion to classify the spatial dependence of phosphate properties. If the ratio is less than 25%, the variable has strong spatial dependence, between 25% and 75%, the variable has moderate spatial dependence and more than 75%, the variable shows weak spatial dependence  .
Figure 8. Variogram model for P2O5 % in the studied areas where, (a) lower member of west Sebaiya, (b) Upper member of west Sebaiya and (c) Um Tondoba mine in east Sebaiya.
Figure 9. Variogram model for ore thickness in the studied areas where, (a) lower member of west Sebaiya, (b) Upper member of west Sebaiya and (c) Um Tondoba mine in east Sebaiya.
Table 2. Variogram and cross variogram parameters of two variables in the studied areas.
Figure 10. Cross variogram model between assay and thickness in the studied areas where, (a) Lower member of west Sebaiya, (b) Upper member of west Sebaiya and (c) Um Tondoba mine in East Sebaiya.
7. Spatial Prediction
Both Ordinary Kriging and Cokriging techniques were used to interpolate the unknown locations by creating map analysis that showing the distribution of P2O5 % in the studied areas as shown in Figures 11-13. These maps are classified into classes, each class represents specific percentage of P2O5 % colured with specific colur, which helps in showing the distribution in different locations at the studied areas. The Ordinary Kriging maps indicate that the P2O5 % tend to be increased closed to middle part of the Lower Member of West Sebaiya where the phosphate medium grade (24% - 26%) is the dominant ratio and the high grade value of P2O5 % disappears.
While in the Upper Member of West Sebaiya the low grade value of P2O5 % appears and the phosphate medium grade (21% - 23%) is the dominant ratio which increases northward.
Also, it is clear from these models the majority of yellow and red colures which high grade value of P2O5 % appears (27% - 32%) at the Eastern part of Elsebaiya area which tend to be increased to the South-East direction. While the Cokriging technique indicates clear different in P2O5 % distribution along all studied areas where some parts of the studied areas turned to a lower concentration part of P2O5 % content and others to a higher concentration of P2O5 % content as shown in cokriging distribution maps. This information will assist mining engineers to take appropriate decisions during exploitation stage.
In this research, ordinary kriging and cokriging methods were used to spatial estimation of P2O5 % in phosphate deposits at Elsebaiya area.
Statistical analysis reflects the clear difference between the distributions and statistical parameters of phosphate within the studied areas. The appearance of high coefficients of variation values at the Upper Member of Western Sebaiya in P2O5 % and thickness is due to the high variability of them.
Figure 11. Kriging model showing P2O5 % distribution for lower member of west Sebaiya where, (a) Ordinary Kriging model, (b) Cokriging model.
Figure 12. Kriging model showing P2O5 % distribution for upper member of west Sebaiya where, (a) Ordinary Kriging model, (b) Cokriging model.
Figure 13. Kriging model showing P2O5 % distribution for um tondoba mine in east Sebaiya where, (a) Ordinary Kriging model, (b) Cokriging model.
Applying geostatistical techniques revealed different behavior of mineralization within the deposit as reflected by the different variogram parameters. The appearance of nugget effect on the variogram models for P2O5 % at the lower member of Western Sebaiya, referred to short scale distance variability to give an attention to vertical variation due to possibility of presence of other elements within the deposits, which could be modeled too.
Results obtained from cokriging method were compared with those from ordinary kriging method which illustrated the clear differences in distributions of P2O5 % in the deposits where it could be useful when planning the mine for production.
 Abbady, A.G.E., Uosif, M.A.M. and El-Taher, A. (2005) Natural Radioactivity and Dose Assessment for Phosphate Rocks from Wadi El-Mashash and El-Mahamid Mines, Egypt. Journal of Environmental Radioactivity, 84, 65-78.
 Zidan, I.H. (2013) Evaluation of Phosphorite and Uranium in Lower Phosphorite Member, Duwi Formation at Kummer Area, South Esna, West Nile Valley, Egypt. Journal of Sedimentological Society of Egypt, 21, 143-154.
 Ahmadi, S.H. and Sedghamiz, A. (2008) Application and Evaluation of Kriging and Cokriging Methods on Groundwater Depth Mapping. Environmental Monitoring and Assessment, 138, 357-368.
 Chang, Y.H., Scrimshaw, M.D., Emmerson, R.H.C. and Lester, J.N. (1998) Geostatistical Analysis of Sampling Uncertainty at the Tollesbury Managed Retreat Site in Blackwater Estuary, Essex, UK: Kriging and Cokriging Approach to Minimise Sampling Density. Science of the Total Environment, 221, 43-57.
 Chien, Y.L., Lee, D.Y., Guo, H.Y. and Houng, K.H. (1997) Geostatistical Analysis of Soil Properties of Mid-West Taiwan Soils. Soil Science, 162, 291-297.