Geostatistical Seismic Inversion with Self-Updating of Local Probability Distributions

Mathematical Geosciences - Three-dimensional subsurface elastic models inverted from seismic reflection data are the basis of the geo-modeling workflow. These models are often used to predict the...     

Indicator-Based Geostatistical Models For Mapping Fish Survey Data

Marine research survey data on fish stocks often show a small proportion of very high-density values, as for many environmental data. This makes the estimation of second-order statistics, such as the variance and the variogram, non-robust. The high fish density values are generated by fish aggregative behaviour, which may vary greatly at small scale in time and space. The high values are thus imprecisely known, both in their spatial occurrence and order of magnitude. To map such data, three indicator-based geostatistical methods were considered, the top-cut model, min–max autocorrelation factors (MAF) of indicators, and multiple indicator kriging. In the top-cut and MAF approaches, the variable is decomposed into components and the most continuous ones (those corresponding to the low and medium values) are used to guide the mapping. The methods are proposed as alternatives to ordinary kriging when the variogram is difficult to estimate. The methods are detailed and applied on a spatial data set of anchovy densities derived from a typical fish stock acoustic survey performed in the Bay of Biscay, which show a few high-density values distributed in small spatial patches and also as solitary events. The model performances are analyzed by cross-validating the data and comparing the kriged maps. Results are compared to ordinary kriging as a base case. The top-cut model had the best cross-validation performance. The indicator-based models allowed mapping high-value areas with small spatial extent, in contrast to ordinary kriging. Practical guidelines for implementing the indicator-based methods are provided.     

Non-stationary Geostatistical Modeling Based on Distance Weighted Statistics and Distributions

A common assumption in geostatistics is that the underlying joint distribution of possible values of a geological attribute at different locations is stationary within a homogeneous domain. This joint distribution is commonly modeled as multi-Gaussian, with correlations defined by a stationary covariance function. This results in attribute maps that fail to reproduce local changes in the mean, in the variance and, particularly, in the spatial continuity. The proposed alternative is to build local distributions, variograms, and correlograms. These are inferred by weighting the samples depending on their distance to selected locations. The local distributions are locally transformed into Gaussian distributions embedding information on the local histogram. The distance weighted experimental variograms and correlograms are able to adapt to local changes in the direction and range of spatial continuity. The automatically fitted local variogram models and the local Gaussian transformation parameters are used in spatial estimation algorithms assuming local stationarity. The resulting maps are rich in nonstationary spatial features. The proposed process implies a higher computational effort than traditional stationary techniques, but if data availability allows for a reliable inference of the local distributions and statistics, a higher accuracy of estimates can be achieved.     

Multivariate Spatial Modeling for Geostatistical Data Using Convolved Covariance Functions

Soil pollution data collection typically studies multivariate measurements at sampling locations, e.g., lead, zinc, copper or cadmium levels. With increased collection of such multivariate geostatistical spatial data, there arises the need for flexible explanatory stochastic models. Here, we propose a general constructive approach for building suitable models based upon convolution of covariance functions. We begin with a general theorem which asserts that, under weak conditions, cross convolution of covariance functions provides a valid cross covariance function. We also obtain a result on dependence induced by such convolution. Since, in general, convolution does not provide closed-form integration, we discuss efficient computation. We then suggest introducing such specification through a Gaussian process to model multivariate spatial random effects within a hierarchical model. We note that modeling spatial random effects in this way is parsimonious relative to say, the linear model of coregionalization. Through a limited simulation, we informally demonstrate that performance for these two specifications appears to be indistinguishable, encouraging the parsimonious choice. Finally, we use the convolved covariance model to analyze a trivariate pollution dataset from California.     

Development and Evaluation of Geostatistical Methods for Non-Euclidean-Based Spatial Covariance Matrices

Mathematical Geosciences - Customary and routine practice of geostatistical modeling assumes that inter-point distances are a Euclidean metric (i.e., as the crow flies) when characterizing spatial...     

Geostatistical Simulation of Regionalized Pore-Size Distributions Using Min/Max Autocorrelation Factors

In many fields of the Earth Sciences, one is interested in the distribution of particle or void sizes within samples. Like many other geological attributes, size distributions exhibit spatial variability, and it is convenient to view them within a geostatistical framework, as regionalized functions or curves. Since they rarely conform to simple parametric models, size distributions are best characterized using their raw spectrum as determined experimentally in the form of a series of abundance measures corresponding to a series of discrete size classes. However, the number of classes may be large and the class abundances may be highly cross-correlated. In order to model the spatial variations of discretized size distributions using current geostatistical simulation methods, it is necessary to reduce the number of variables considered and to render them uncorrelated among one another. This is achieved using a principal components-based approach known as Min/Max Autocorrelation Factors (MAF). For a two-structure linear model of coregionalization, the approach has the attractive feature of producing orthogonal factors ranked in order of increasing spatial correlation. Factors consisting largely of noise and exhibiting pure nugget–effect correlation structures are isolated in the lower rankings, and these need not be simulated. The factors to be simulated are those capturing most of the spatial correlation in the data, and they are isolated in the highest rankings. Following a review of MAF theory, the approach is applied to the modeling of pore-size distributions in partially welded tuff. Results of the case study confirm the usefulness of the MAF approach for the simulation of large numbers of coregionalized variables.     

Towards Geostatistical Learning for the Geosciences: A Case Study in Improving the Spatial Awareness of Spectral Clustering

Mathematical Geosciences - The particularities of geosystems and geoscience data must be understood before any development or implementation of statistical learning algorithms. Without such...     

Correction to: Development and Evaluation of Geostatistical Methods for Non-Euclidean-Based Spatial Covariance Matrices

Mathematical Geosciences - The original version of this article unfortunately contained a mistake in equation 9.     

Geostatistical Rock Physics Inversion for Predicting the Spatial Distribution of Porosity and Saturation in the Critical Zone

Mathematical Geosciences - Understanding the subsurface structure and function in the near-surface groundwater system, including fluid flow, geomechanical, and weathering processes, requires...     

Quantitative Control of Migration: A Geostatistical Attempt

This paper is devoted to a geostatistical attempt at modeling migration errors when localizing a reflector in the ground. Starting with a probabilistic velocity model and choosing the simple geometrical optics background for the wave propagation in such media, we give the expression of the errors. This may be quantified provided the covariance of the velocity field is known. Variance of arrival times at constant offset is related to the covariance of the velocity field at hand. A practical application is given in the same paragraph. After that we give a typical schema for migration and uncertainty modeling: starting with seismic data, we make the weak seismic inversion. We then obtain the covariance of the velocity field that we use for simulating migration errors. The main issues of this methodology are discussed in the last paragraph.     

Optimized Sample Schemes for Geostatistical Surveys

Sample schemes used in geostatistical surveys must be suitable for both variogram estimation and kriging. Previously schemes have been optimized for one of these steps in isolation. Ordinary kriging generally requires the sampling locations to be evenly dispersed over the region. Variogram estimation requires a more irregular pattern of sampling locations since comparisons must be made between measurements separated by all lags up to and beyond the range of spatial correlation. Previous studies have not considered how to combine these optimized schemes into a single survey and how to decide what proportion of sampling effort should be devoted to variogram estimation and what proportion devoted to kriging An expression for the total error in a geostatistical survey accounting for uncertainty due to both ordinary kriging and variogram uncertainty is derived. In the same manner as the kriging variance, this expression is a function of the variogram but not of the sampled response data. If a particular variogram is assumed the total error in a geostatistical survey may be estimated prior to sampling. We can therefore design an optimal sample scheme for the combined processes of variogram estimation and ordinary kriging by minimizing this expression. The minimization is achieved by spatial simulated annealing. The resulting sample schemes ensure that the region is fairly evenly covered but include some close pairs to analyse the spatial correlation over short distances. The form of these optimal sample schemes is sensitive to the assumed variogram. Therefore a Bayesian approach is adopted where, rather than assuming a single variogram, we minimize the expected total error over a distribution of plausible variograms. This is computationally expensive so a strategy is suggested to reduce the number of computations required     

Correction: Geostatistical Rock Physics Inversion for Predicting the Spatial Distribution of Porosity and Saturation in the Critical Zone

Mathematical Geosciences -     

Measures of Parameter Uncertainty in Geostatistical Estimation and Geostatistical Optimal Design

Studies of site exploration, data assimilation, or geostatistical inversion measure parameter uncertainty in order to assess the optimality of a suggested scheme. This study reviews and discusses measures for parameter uncertainty in spatial estimation. Most measures originate from alphabetic criteria in optimal design and were transferred to geostatistical estimation. Further rather intuitive measures can be found in the geostatistical literature, and some new measures will be suggested in this study. It is shown how these measures relate to the optimality alphabet and to relative entropy. Issues of physical and statistical significance are addressed whenever they arise. Computational feasibility and efficient ways to evaluate the above measures are discussed in this paper, and an illustrative synthetic case study is provided. A major conclusion is that the mean estimation variance and the averaged conditional integral scale are a powerful duo for characterizing conditional parameter uncertainty, with direct correspondence to the well-understood optimality alphabet. This study is based on cokriging generalized to uncertain mean and trends because it is the most general representative of linear spatial estimation within the Bayesian framework. Generalization to kriging and quasi-linear schemes is straightforward. Options for application to non-Gaussian and non-linear problems are discussed.     

Geostatistical Data Integration Model for Contamination Assessment

Soil contamination assessments can be improved with new methods aimed at the accurate estimation of the volume and extension of contaminated soil to be remediated. Geostatistical models that use secondary information to characterize soil contamination are incorporated into a new integration model to provide accurate three-dimensional maps. The proposed integration model is based on a stochastic inversion approach and uses sequential indicator simulation. A two-dimensional reference image representing the areal extension of the contamination is combined with local measurements of contamination in the vertical direction, to render a three-dimensional contamination map. To demonstrate how well the integration model performs, the case study presented focuses on geophysical data and how it can be integrated with soil contamination measurements to improve the characterization of a contaminated site. The results show that the model reproduces successfully the reference image thus providing an accurate three-dimensional contamination map.     

Spatial Variability of Rock Depth in Bangalore Using Geostatistical, Neural Network and Support Vector Machine Models

Geospatial technology is increasing in demand for many applications in geosciences. Spatial variability of the bed/hard rock is vital for many applications in geotechnical and earthquake engineering problems such as design of deep foundations, site amplification, ground response studies, liquefaction, microzonation etc. In this paper, reduced level of rock at Bangalore, India is arrived from the 652 boreholes data in the area covering 220 km². In the context of prediction of reduced level of rock in the subsurface of Bangalore and to study the spatial variability of the rock depth, Geostatistical model based on Ordinary Kriging technique, Artificial Neural Network (ANN) and Support Vector Machine (SVM) models have been developed. In Ordinary Kriging, the knowledge of the semi-variogram of the reduced level of rock from 652 points in Bangalore is used to predict the reduced level of rock at any point in the subsurface of the Bangalore, where field measurements are not available. A new type of cross-validation analysis developed proves the robustness of the Ordinary Kriging model. ANN model based on multi layer perceptrons (MLPs) that are trained with Levenberg–Marquardt backpropagation algorithm has been adopted to train the model with 90% of the data available. The SVM is a novel type of learning machine based on statistical learning theory, uses regression technique by introducing loss function has been used to predict the reduced level of rock from a large set of data. In this study, a comparative study of three numerical models to predict reduced level of rock has been presented and discussed.     

Impacts of Hypoxia on Zooplankton Spatial Distributions in the Northern Gulf of Mexico

The northern Gulf of Mexico (NGOMEX) was surveyed to examine the broad-scale spatial patterns and inter-relationships between hypoxia (<2?mg?L^?1 dissolved oxygen) and zooplankton biovolume. We used an undulating towed body equipped with sensors for conductivity, temperature, depth, oxygen, fluorescence, and an optical plankton counter to sample water column structure, oxygen, and zooplankton at high spatial resolution (1?m??vertical; 0.25?C1?km??horizontal). We contrast the distribution of zooplankton during summer surveys with different freshwater input, stratification, and horizontal and vertical extent of bottom-water hypoxia. Bottom-water hypoxia did not appear to influence the total amount of zooplankton biomass present in the water column or the areal integration of zooplankton standing stock in the NGOMEX region surveyed. However, where there were hypoxic bottom waters, zooplankton shifted their vertical distribution to the upper water column during the day where they normally would reside in deeper and darker waters. When bottom waters were normoxic (>2?mg?L^?1 dissolved oxygen), the daytime median depth of the water column zooplankton was on average 7?m deeper than the median depth of zooplankton in water columns with hypoxic bottom waters. A reduction in larger zooplankton when there were hypoxic bottom waters suggests that if zooplankton cannot migrate to deeper, darker water under hypoxic conditions, they may be more susceptible to size-selective predation by visual predators. Thus, habitat compression in the northern Gulf of Mexico due to hypoxic bottom water may have implications for trophic transfer by increasing the contact between predators and prey.