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.     

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

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

氧化镍、二氧化钛高压相变及电导率压力效应

采用金刚石对顶砧压腔静高压技术，对氧化镍、二氧化钛在０-３２ＧＰａ压力范围内的电导率进行了测量．结果表明，电导率－压力关系曲线上的一些特征点与同一物质的高压相交点具有一致性．     

Nonequilibrium thermodynamics of pressure solution

This paper is concerned with the thermodynamic theory of solution and precipitation processes in wet crustal rocks and with the mechanism of steady pressure-solution slip in contact zones, such as grain-to-grain contacts, fracture surfaces, and permeable gouge layers, that are infiltrated by a mobile aqueous solution phase. A local dissipation jump condition at the phase boundary is fundamental to identifying the thermodynamic force driving the solution and precipitation process and is used here in setting up linear phenomenological relations to model near-equilibrium phase transformation kinetics. The local thermodynamic equilibrium of a stressed pure solid in contact with its melt or solution phase is governed by Gibbs's relation, which is rederived here, in a manner emphasizing its independence of constitutive assumptions for the solid while neglecting surface tension and diffusion in the solid. Fluid-infiltrated contact zones, such as those formed by rough surfaces, cannot generally be in thermodynamic equilibrium, especially during an ongoing process of pressure-solution slip, and the existing equilibrium formulations are incorrect in overlooking dissipative processes tending to eliminate fluctuations in superficial free energies due to stress concentrations near asperities, defects, or impurities. Steady pressure-solution slip is likely to exhibit a nonlinear dependence of slip rate on shear stress and effective normal stress, due to a dependence of the contact-zone state on the latter. Given that this dependence is negligible within some range, linear relations for pressure-solution slip can be derived for the limiting cases of diffusion-controlled and interface-reaction-controlled rates. A criterion for rate control by one of these mechanisms is set by the magnitude of the dimensionless quantityk/2C _pD, wherek is the interfacial transfer coefficient, is the mean diffusion path length,C _p is the solubility at pressurep, andD is the mass diffusivity.