• The Journal of Engineering Geology
• /
• v.3 no.1
• /
• pp.39-49
• /
• 1993

Intrinsic random function of order k(IRF-k) was used to estimate groundwater levels of nonstationaav random functions. The accuracy of IRF-k was compared to that of ordraarv krigrng assuming that the data of groundwater levels compose a stafionarv random function. Cross validation and statistical errors show that IRF-k is superior to orcinarv '(riging for the estimation of water levels. IRF-k and ordinary kriging made different contour and 3-D surface maps. The maps of IRF-k are more accurate than those of ordinary kriging.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.352-361
• /
• 2003

We have evaluated the extent of saltwater intrusion from electrical resistivity distribution in a coastal aquifer system in the southeastern part of Busan, Korea. This aquifer system is divided into four layers according to the hydrogeologic characteristics and the horizontal extent of intruded saltwater is determined at each layer through the geostatistical interpretation of electrical resistivity data. In order to define the statistical structure of electrical resistivity data, variogram analysis is carried out to obtain best generalized covariance models. IRF-k (intrinsic random function of order k) kriging is performed with covariance models to produce the plane of spatial mean resistivities. The kriged estimates are evaluated by cross validation to show a good agreement with the true values and the statistics of cross validation represented low errors for the estimates. In the resistivity contour maps more than 5 m below the surface, we can see a dominant direction of saltwater intrusion beginning from the east side. The area of saltwater intrusion increases with depth. The northeast side has low resistivities less than 5 ohm-m due to the presence of saline water in the depth range of 20 m through 70 m. These results show that the application of geostatistical technique to electrical resistivity data is useful for assessing saltwater intrusion in a coastal aquifer system.

• The Journal of Engineering Geology
• /
• v.3 no.1
• /
• pp.51-61
• /
• 1993

Restricted maximum likelihood(RML) method was used to determine the parameters of generalized covariance, and universal krigig with RML was applied to estimate a groundwater level distribution of nonstationarv random function. Universal kriging with RML was compared to IRF-k with weighted least squares method for the comparison of their accuracies. Cross validation shows that two methods have nearly the same ability for the estimation of groundwater levels. Scattergram of estimates versus true values and contour maps of groundwater levels have nearly the same results. The reason why two methods produced the same results is thought to be the non-Gaussian distribution and the snaall number of sample data.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2001.09a
• /
• pp.92-96
• /
• 2001

Although the problem of seawater intrusion at the coastal aquifer was recognized before over one hundred years at the coastal aquifer, much groundwater keep on being salinitized by several reasons such as groundwater exhaustion, coastalline change, and human activities. The horizontal and vertical electrical soundings and geostatistical methods were used to define the local characteristics of saltwater intrusion and to estimate the saltwater interface in the southeastern area of the Pusan City. The 24 points of the Schlumberger vertical electrical soundings(VES) to loom depth and the 2 lines of dipole-dipole horizontal soundings are peformed. The resistivity data have lognormal distributions. The horizontal extents of saline water intrusion were estimated from the inversion of horizontal prospecting data. Lognormal ordinary kriging is used in A-A' resistivity profiles on May and July because the data have stationary models in semivariograms. Lognormal IRF-k kriging is used for the isopleth maps using vertical resistivity data. The 10 ohm-m resistivity line on the isopleth maps of 21m, 30m, 50m, and 70m depth using resisitivity data measured in July is sifted to the east, cpomparing that of the isopleth maps measured in May. The kriged vertical and horizontal resistivity isopleth maps suggested that the geostatistical methods can be used to define the variation of earth resistivity distribution at the saltwater interface.