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.     

Data analysis of a revisited exploration geochemical dataset of quartz porphyrites from SW Jordan using GIS techniques

The spatial distribution of Cu and its pathfinder elements in quartz porphyrites along Wadi Araba was studied in order to gain an insight about the geologic framework of Cu mineralization in the area. A geochemical exploration data set of 242 samples was revisited. Multivariate statistical and geostatistical analysis were performed on the data set. The objectives of the study were to define Cu relationships with its possible pathfinder elements and to figure out whether geostatistical treatment of exploration data would help telescoping already existing Cu mineralization in Wadi Araba, SW Jordan. Factor analysis and hierarchical cluster analysis were performed on the data set. The results of which indicate that seven elements, namely, Ba, Cr, Mn, Ni, Sc, Sr and V can be considered as pathfinder elements for Cu. Significant correlations do exist between Cu and these elements. The encountered experimental semivariograms of the defined pathfinders were fit to exponential models with varying nugget effects with the exception of Mn which was fit to the spherical model. Ordinary kriging was found to be most appropriate to generate geochemical exploration maps. The applied methods succeeded in deriving valuable information and visualizing the already existing Cu mineralization in Wadi Abu Khusheiba area.     

Estimation of transit trip production using Factorial Kriging with External Drift: an aggregated data case study

Studies in transportation planning routinely use data in which location attributes are an important source of information. Thus, using spatial attributes in urban travel forecasting models seems reasonable. The main objective of this paper is to estimate transit trip production using Factorial Kriging with External Drift (FKED) through an aggregated data case study of Traffic Analysis Zones in São Paulo city, Brazil. The method consists of a sequential application of Principal Components Analysis (PCA) and Kriging with External Drift (KED). The traditional Linear Regression (LR) model was adopted with the aim of validating the proposed method. The results show that PCA summarizes and combines 23 socioeconomic variables using 4 components. The first component is introduced in KED, as secondary information, to estimate transit trip production by public transport in geographic coordinates where there is no prior knowledge of the values. Cross-validation for the FKED model presented high values of the correlation coefficient between estimated and observed values. Moreover, low error values were observed. The accuracy of the LR model was similar to FKED. However, the proposed method is able to map the transit trip production in several geographical coordinates of non-sampled values.     

Controlling factors on the morphology and spatial distribution of methane-related tubular concretions – Case study of an Early Eocene seep system

Up to 10 m in length and >1 m in diameter tubular, calcite-cemented sandstone concretions are hosted by the faulted Dikilitash unconsolidated sands and sandstones. These structures document shallow subsurface pathways of Early Eocene methane seepage in the Balkan Mountains foreland (NE Bulgaria). Their exceptional exposure allowed a unique study of the factors governing the morphology and spatial distribution of such fossilized fluid conduits. The large dimensions and subvertical, cylindrical shape of the most common tube type primarily reflects the buoyancy-driven, vertical path of an ascending gas-bearing fluid through permeable, mainly unconsolidated sandy host sediments. Tube morphology was also influenced by local stratigraphic anisotropies and might as well document differences in former seepage conditions. Mapping of >800 tubular concretions showed the NNW–NNE elongation and alignment of tube clusters and massive cemented sandstone structures. This suggests that Paleogene fault systems played a major role in directing the movement of fluids. However, within a single tube cluster, tubes are preferentially aligned, over distances up to 50 m along directions at an angle between 10° and 36° with respect to the inferred NNW–NNE, cluster parallel fault traces. In addition, cylindrical tubes of analogue dimensions are aligned over distances >100 m along N15° to N25°-oriented directions. It is hypothesized that this spatial geometry of tubular concretions reflects the complex geometry of deformations structures in fault damage zones along which fluids were preferentially channelled.