Multivariate analysis of heavy metal contents of sediments from Gumusler creek, Nigde, Turkey

In this study, heavy metal contents of samples from Gumusler creek in Turkey were studied and the metal contamination characteristics were investigated. In this respect, considering the pollutant sites in the area, systematic sediment samples were collected in a zone starting from the manifestation part of the water to the entrance of the Karasu creek in Gumusler town. Samples were taken from lower section of the river bed at 30 stations along Gumusler creek, 13 km in length and their heavy metal contents were analyzed with XRF Spectrometer. Correlation coefficients, element coefitic coefficient correlation, dendogram hierarchical cluster, model summary and Annova analysis statistical methods were applied to data. Strong positive correlations were determined for some elements which are believed to have possibly the same origin. In addition, mineralizations in the area are thought to cause variation in metal contents. Results of chemical analysis show that soil limit values and clark values were exceeded. The heavy metal accumulation in the creek is believed to be derived from non-operated Sb-Hg-W and Fe quarries.     

Spatial variability in petrophysical properties in Upper Maastrichtian chalk outcrops at Stevns Klint, Denmark

The spatial variability of porosity and permeability was determined for a section in the uppermost Maastrichtian chalk exposed in the Sigerslev quarry at Stevns Klint. The aims were threefold: (1) to quantify the spatial variability pattern and its link to geology by applying systematic variogram analysis procedure, (2) to detect potential cyclicity in the petrophysical properties in this pure and overtly homogeneous chalk, and (3) to compare the section with chalk reservoirs in the North Sea in order to asses to what extent this onshore sequence has a potential as an analogue of the offshore, more deeply buried chalk reservoirs. The interval is of comparable stratigraphic age to the uppermost reservoir zones in the Tor Formation of the North Sea chalk reservoirs. The variability and spatial correlation of porosity and permeability in both horizontal and vertical directions are interpreted in a geological context and show indications of small-scale heterogeneity at 15–25 cm scale, but the clear cyclic layering described from other chalk deposits is not recognised at this locality. The investigated outcrop is not a close analogue to North Sea reservoir chalk, but some aspects are common including basic material properties, porosity/permeability trends and the variability pattern. The outcrop has a potential as analogue for some of the onshore subsurface chalk successions at shallow burial depth that form important aquifers.     

Using Sequential Self-Calibration Method to Identify Conductivity Distribution: Conditioning on Tracer Test data

An iterative inverse method, the sequential self-calibration method, is developed for mapping spatial distribution of a hydraulic conductivity field by conditioning on nonreactive tracer breakthrough curves. A streamline-based, semi-analytical simulator is adopted to simulate solute transport in a heterogeneous aquifer. The simulation is used as the forward modeling step. In this study, the hydraulic conductivity is assumed to be a deterministic or random variable. Within the framework of the streamline-based simulator, the efficient semi-analytical method is used to calculate sensitivity coefficients of the solute concentration with respect to the hydraulic conductivity variation. The calculated sensitivities account for spatial correlations between the solute concentration and parameters. The performance of the inverse method is assessed by two synthetic tracer tests conducted in an aquifer with a distinct spatial pattern of heterogeneity. The study results indicate that the developed iterative inverse method is able to identify and reproduce the large-scale heterogeneity pattern of the aquifer given appropriate observation wells in these synthetic cases.     

Modèles géostatistiques de concentrations ou de débits le long des cours d'eau

Estimating concentrations or flow rates along a stream network requires specific models. Two classes of models, recently proposed in the literature, are generalized, to the intrinsic case in particular. We present a global construction by ‘streams’, i.e. on the whole set of paths between sources and outlet. Combining stationary or intrinsic one-dimensional random functions leads to stationary or intrinsic models on segments, with discontinuities at the forks. A construction from outlet to sources, leads to stationary or intrinsic models on each stream, without any discontinuity at the forks. The linear variogram is found as a particular case. The extension to the linear model of coregionalization is immediate, allowing a multivariate modelling of concentrations. To cite this article: C. de Fouquet, C. Bernard-Michel, C. R. Geoscience 338 (2006).     

Fracture and bedding plane control on groundwater flow in a chalk aquitard

Shallow groundwater in the northern Negev desert of Israel flows preferentially through a complex system of discontinuities. These discontinuities intersect what would otherwise be a massive, low-conductivity, high-porosity Eocene chalk. Vertical fractures and horizontal bedding planes were observed and mapped along approximately 1,200 m of scanline, 600 m of core and 30 two-dimensional trace planes. A bimodal distribution of size exists for the vertical fractures which occur as both single-layer fractures and multi-layer fractures. A bimodal distribution of log transmissivity was observed from slug tests conducted in packed-off, vertical intervals within the saturated zone. The different flow characteristics between the horizontal bedding planes and vertical-type fractures appear to be the cause of the bimodality. Two distinct conceptual models (discrete fracture network) were developed based on the fracture orientation, size, intensity and transmissivity statistics derived from field data. A correlation between fracture size and hydraulic aperture was established as the basis for calibrating the simulated model transmissivity to the field observations. This method of defining transmissivity statistically based on prior information is shown to be a reasonable and workable alternative to the usual conjecture approach towards defining transmissivity in a fractured-rock environment.     

Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data

This paper presents the extension of the self-calibrating method to the coupled inverse modelling of groundwater flow and mass transport. The method generates equally likely solutions to the inverse problem that display the variability as observed in the field and are not affected by a linearisation of the state equations. Conditioning to the state variables is measured by an objective function including, among others, the mismatch between the simulated and measured concentrations. Conditioning is achieved by minimising the objective function by gradient-based methods. The gradient contains the partial derivatives of the objective function with respect to: log conductivities, log storativities, prescribed heads at boundaries, retardation coefficients and mass sources. The derivatives of the objective function with respect to log conductivity are the most cumbersome and need the most CPU-time to be evaluated. For this reason, to compute this derivative only advective transport is considered. The gradient is calculated by the adjoint-state method. The method is demonstrated in a controlled, synthetic study, in which the worth of concentration data is analysed. It is shown that concentration data are essential to improve transport predictions and also help to improve aquifer characterisation and flow predictions, especially in the upstream part of the aquifer, even in the case that a considerable amount of other experimental data like conductivities and heads are available. Besides, conditioning to concentration data reduces the ensemble variances of estimated transmissivity, hydraulic head and concentration.     

Space-time mapping of soil salinity using probabilistic bayesian maximum entropy

The mapping of saline soils is the first task before any reclamation effort. Reclamation is based on the knowledge of soil salinity in space and how it evolves with time. Soil salinity is traditionally determined by soil sampling and laboratory analysis. Recently, it became possible to complement these hard data with soft secondary data made available using field sensors like electrode probes. In this study, we had two data sets. The first includes measurements of field salinity (EC_a) at 413 locations and 19 time instants. The second, which is a subset of the first (13 to 20 locations), contains, in addition to EC_a, salinity determined in the laboratory (EC_2.5). Based on a procedure of cross-validation, we compared the prediction performance in the space-time domain of 3 methods: kriging using either only hard data (HK) or hard and mid interval soft data (HMIK), and Bayesian maximum entropy (BME) using probabilistic soft data. We found that BME was less biased, more accurate and giving estimates, which were better correlated with the observed values than the two kriging techniques. In addition, BME allowed one to delineate with better detail saline from non-saline areas.     

基于地统计方法的气候要素空间插值研究

在回顾了地统计学的产生、发展及其基本原理的基础上,对于目前众多可以提供计算网格的气候要素的空间插值方法中,具体探讨了普通克里格法和协同克里格法。将甘肃省1961—1990年30年平均降水量和蒸发量作为区域化变量,根据不同的半变异函数理论模型,采用普通克里格法和双变量协同克里格法,通过对比分析得到：（1）不论是多年均降水量还是多年平均蒸发量在空间上都呈现明显的梯度变化,二者的空间变程都很大,而降水量变化幅度更大。降水量从东南部向西北部逐渐减少,蒸发量则相反,从东南向西北逐渐增加。（2）基于地统计的插值方法,根据半变异函数云图和试验方差最小的原理,选择合适的半变异函数理论模型进行变量的空间插值,能够较好地模拟区域化变量的空间连续分布格局,并取得较好的效果。对比普通克里格法和协同克里格法,后者增加了高度对降水量和蒸发量的影响,在空间分布上更为合理,插值的精度也要明显好于普通克里格法。（3）采用地统计方法虽然在总体上能够较好地反映气候要素的空间分布格局,但检验显示,两种方法空间插值的精度都还不是很高,插值的精度还有待进一步提高。     

青藏高原气温分布的空间插值方法比较

使用反距离平方、趋势面、Kriging插值、Cokriging插值和综合方法对青藏高原1961-1990年30年平均1月气温进行空间插值比较研究，其中后两种方法能够把影响青藏高原气温分布的关键因素——高程置于插值算法之中。反距离平方和趋势面插值的结果都与实际情况相差较远；普通Kriging插值能够反映出青藏高原气温分布的一定的空间结构，但结果依然不好；由于考虑了高度变量，Cokriging插值表现出一定程度的性能改进，但因为台站海拔高度偏低，插值结果依然不理想。前4种方法效果不好的原因主要是因为青藏高原上气象台站稀少且高原西北地区无气象台站。综合方法把气温分解为结构化分量和随机分量，使用直减率把气温订正到同一海拔高度后，再对它们做Kriging插值分析，其结果较为正确地反映了青藏高原气温空间分布的特征，误差远小于其它方法。研究结果表明：样本本身的空间分布是影响插值精度的重要因素，合理的采样设计是必要的前提；对于台站稀少的地区，必须把随机插值方法和确定性方法结合起来估计气候变量的空间分布。同时给出了青藏高原1月的气温空间分布状况。     

地统计学方法在土壤学中的应用

地统计学目前在土壤科学中得到广泛的应用和发展,成为认识土壤特征的一个重要工具.地统计学的空间变异函数和克里格插值等方法是土壤性状分析的主要手段,空间变异函数主要用于描述土壤理化性状空间变异特性,不同的插值方法可优化田间试验设计和田间采样方案,克里格插值方法则特别适用于未测量点土壤属性值的估测等.近年插值方法又被广泛应用于确定区域土壤环境容量和土壤质量标准,随机模拟则用于对土壤特性进行不确定性估计等.因而,地统计学方法对我国大量土壤学资料的整合与分析具有极大的应用前景.