Improving kriging of groundwater level data using nonlinear normalizing transformations—a field application

Abstract

The present research study investigates the application of nonlinear normalizing data transformations in conjunction with ordinary kriging (OK) for the accurate prediction of groundwater level spatial variability in a sparsely-gauged basin. We investigate three established normalizing methods, Gaussian anamorphosis, trans-Gaussian kriging and the Box-Cox method to improve the estimation accuracy. The first two are applied for the first time to groundwater level data. All three methods improve the mean absolute prediction error compared to the application of OK to the non-transformed data. In addition, a modified Box-Cox transformation is proposed and applied to normalize the hydraulic heads. The modified Box-Cox transformation in conjunction with OK is found to be the optimal spatial model based on leave-one-out cross-validation. The recently established Spartan semivariogram family provides the optimal model fit to the transformed data. Finally, we present maps of the groundwater level and the kriging variance based on the optimal spatial model.

Editor D. Koutsoyiannis; Associate editor A. Montanari

Citation Varouchakis, E.A., Hristopoulos, D.T., and Karatzas, G.P., 2012. Improving kriging of groundwater level data using nonlinear normalizing transformations—a field application. Hydrological Sciences Journal, 57 (7), 1404–1419.     

边界单元法在电法勘探中的应用

本文用加权剩余法导出点源二维地电断面的边界单元法积分方程。通过付氏反变换,可获得三维空间的电位。若干模型上正演计算结果与解析解完全一致。对实测的ρ_s曲线进行了地形攻正及选择法反演拟合,解释结果与实际地质断面符合得很好。     

大气水—地表水—土壤水—地下水相互转化关系的试验研究

本文重点研究了有作物条件下不同供水农田的“四水”相互转化,突出了有作物条件下降水入渗、潜水蒸发的差异性,土壤水在农业生态环境中的重要性,为水资源多目标开发,探索新路子。     

世界全图数字化数据处理的解析变换方法

鉴于世界全图数字化数据是各类专业信息系统广泛用于空间信息定位的载体和基础,本文详细讨论了1:1400万世界全图数字化数据处理解析变换方法的正解变换模型、统一投影坐标系模型和反解变换模型,并给出了算例。最后证明,本文提供的更解变换模型的计算精度完全能满足地图数据库要求。     

ELIMINATING COMPLEX EIGENVECTORS AND EIGENVALUES IN MULTIWAY ANALYSES USING THE DIRECT TRILINEAR DECOMPOSITION METHOD

The direct trilinear decomposition method(DTDM)is an algorithm for performing quantitative curveresolution of three-dimensional data that follow the so-called trilinear model,e.g.chromatography-spectroscopy or emission-excitation fluorescence.Under certain conditions complexeigenvalues and eigenvectors emerge when the generalized eigenproblem is solved in DTDM.Previouspublications never treated those cases.In this paper we show how similarity transformations can be usedto eliminate the imaginary part of the complex eigenvalues and eigenvectors,thereby increasing theusefulness of DTDM in practical applications.The similarity transformation technique was first used byour laboratory to solve the similar problem in the generalized rank annihilation method(GRAM).Because unique elution profiles and spectra can be derived by using data matrices from three or moresamples simultaneously,DTDM with similarity transformations is more efficient than GRAM in the casewhere there are many samples to be investigated.     

An examination of the use of the logratio transformation for the testing of endmember hypotheses

Given a compositional dataset in the absence of any prior information on any mixing process which may have formed it, a complete analysis of mixtures determines three distinct types of estimates in order. These are: (i) the estimate of the number of endmembers or fixed source compositions, of which all the sample compositions of the dataset must be approximate mixtures; (ii) the estimated compositions for each of these chosen number of endmembers; and (iii) the estimated contributions of each of these endmember estimates to each sample. Traditionally, the estimate for the number of endmembers has been assessed either by mapping or by inspection of the coefficients of determination between the observed and estimated variables. Mapping entails the plotting on a map of the region from which the samples were taken, either the contours of the contributions of each endmember to each sample, or some other portrayal of the distribution of endmember abundances. Because it requires the complete analysis, assessment by this method is too elaborate except for final confirmation and display. Alternatively, choosing a number of endmembers, which result in suitability high coefficients of determination for all or most variables, may account for elements which are not part of the conjectured mixing process or, worse, may result in the identification of endmembers which may never in fact have existed. Such an error is similar to overspecifying a multiple regression model. So, the obvious starting point from which to assess the validity, or otherwise choice of endmember numbers, is to examine the matrix of residuals. The differences between the logratio-transformed observed and estimated data form an array of residual logratios. A linear combination of these may be formed for each sample, which, under a random perturbation assumption, should follow a univariate normal distribution. Whether or not this scalar is normal can be readily tested. It can also be examined graphically for such desirable qualities as symmetry when the test for normality may be too severe. This procedure is employed to assess the decompositions of the U.S.G.S. Mid-Pacific data and the Nazca Plate Surface sediments.This paper was presented at the 18th Geochautauqua, Newark, Delaware, 13–14 October 1989.     

STUDY OF THE POSSIBILITIES OF USING AN AIR MASS TRANSFORMATION MODEL IN TAIYUAN

An AMT-model,consisting of a trajectory model and a one-dimensional boundary layer model,is tested fortrajectories arriving in Taiyuan to study the possibility of using it in Taiyuan.The sensitivity of the model tothe different processes was studied.Some parameters of the model were modified for the purpose of forecast-ing in specific mountainous terrain and dry climate conditions.Results of examples which we have workedout for Taiyuan circumstances for the periods of July(summer)1985 and January(winter)1986,show that the12h runs of the AMT-model are able to reproduce(on historical data)the sounding of Taiyuan.The AMT-modelcontributes fruitfully to short-range weather forecasts(12—36h ahead)during periods of severe air pollution andwhen cold waves occur.     

TWO－DIMENSIONAL MATHEMATICAL MODEL FOR STUDYING EROSION AND SEDIMENTATION OF A WANDERING STRETCH IN THE LOWER YELLOW RIVER

ＴＷＯ－ＤＩＭＥＮＳＩＯＮＡＬＭＡＴＨＥＭＡＴＩＣＡＬＭＯＤＥＬＦＯＲＳＴＵＤＹＩＮＧＥＲＯＳＩＯＮＡＮＤＳＥＤＩＭＥＮＴＡＴＩＯＮＯＦＡＷＡＮＤＥＲＩＮＧＳＴＲＥＴＣＨＩＮＴＨＥＬＯＷＥＲＹＥＬＬＯＷＲＩＶＥＲ￥ＺＨＡＮＧＳｈｉｑｉｒ（Ｓｅｎｉｏｒ...     

空间后交的一种新算法

根据矩阵理论，推导了空间后交的一种新算法．这种方法简单易行，无需线性化，其计算结果、精度与传统方法一致．     

1966年云南东川6.5级地震的热异常时空结构分析

本文分析了１９６６年云南东川６．５级地震前后浅层地温异常的时空结构及降水异常的水平结构。分析表明，震前数月浅层地温稳定维持正异常，发震前一个月激剧增温，发震当月达最大，震后降低；震前震中区为负地温距平区，随发震时刻临近，负距平区迅速演变为正距平区，而且强度和面积都不断增大，至发震时达最大，震后不断减小；热异常是从深层传向地表的；震中区震前为旱区，震后为涝区。本文还讨论了震源处热信息传到地表的条件和可能机制，认为在水热活动异常的断裂带上，至少对于浅源地震，对流传热是一种较为可能的传热机制。