Multiple Random Walk Simulation: A Fast Method to Map Grade Uncertainty with Large Datasets

Multiple Random Walk Simulation consists of a methodology adapted to run fast simulations if close-spaced data are abundant (e.g., short-term mining models). Combining kriging with the simulation of random walks attempts to approximate traditional simulation algorithm results but at a computationally faster way when there is a large amount of conditioning samples. This paper presents this new algorithm illustrating the situations where the method can be used properly. A synthetic study case is presented in order to illustrate the Multiple Random Walk Simulation and to analyze the speed and goodness of its results against the ones from using Turning Bands Simulation and Sequential Gaussian Simulation.     

Growth and mortality of juvenile winter flounder in two New England estuaries

A length-based model for calculating growth and mortality of juvenile winter flounder (Pseudopleuronectes americanus) populations has been developed. This model is based on work by Sullivan et al. (1990) and incorporates the von Bertalanffy growth equation, including stochasticity in growth, and a mortality rate that decreases exponentially with size. The length-based model was fit to observed size-frequency distributions, and model likelihood profiles were generated to produce 95% confidence intervals about parameter estimates. We analyzed size-frequency distributions of 3 to 15 cm juvenile winter flounder, collected with a 1-m beam trawl, at monthly intervals from June to October during 1993 and 1994. Growth rates were higher at a contaminated site, New Haven Harbor, than at a clean site, the Connecticut River estuary, however, the parameter estimates had overlapping 95% confidence intervals. Mortality rates were similar at the two sites.     

Salinisation de l’aquifère libre de la Chaouia côtière (Azemmour-Tnine Chtouka), Maroc

L’aquifère libre de la Chaouia côtière constitue un exemple des aquifères les plus exploités au Maroc. Ce travail est consacré à l’étude des processus de la salinisation des eaux souterraines par l’analyse physico-chimique de 39 puits répartis dans la zone. Deux types de faciès ont été révélés, l’un est chloruré-sodique dans la frange côtière ; l’autre faciès est de type chloruré bicarbonaté-calcique caractérisant les eaux exploitées dans les calcaires marneux du Crétacé. L’influence marine (aérosols et intrusion marine), la dissolution/précipitation de la roche aquifère et l’infiltration des eaux d’irrigation, sont parmi les causes de l’augmentation de la salinité des eaux souterraines, en plus de l’exploitation excessive de l’eau souterraine.     

Uncertainty analysis and probabilistic segmentation of electrical resistivity images: the 2D inverse problem

In this paper, we present the uncertainty analysis of the 2D electrical tomography inverse problem using model reduction and performing the sampling via an explorative member of the Particle Swarm Optimization family, called the Regressive‐Regressive Particle Swarm Optimization. The procedure begins with a local inversion to find a good resistivity model located in the nonlinear equivalence region of the set of plausible solutions. The dimension of this geophysical model is then reduced using spectral decomposition, and the uncertainty space is explored via Particle Swarm Optimization. Using this approach, we show that it is possible to sample the uncertainty space of the electrical tomography inverse problem. We illustrate this methodology with the application to a synthetic and a real dataset coming from a karstic geological set‐up. By computing the uncertainty of the inverse solution, it is possible to perform the segmentation of the resistivity images issued from inversion. This segmentation is based on the set of equivalent models that have been sampled, and makes it possible to answer geophysical questions in a probabilistic way, performing risk analysis.     

Climatic signals in lacustrine deposits of the Upper Yacoraite Formation,Western Argentina: Evidence from clay minerals,analcime, dolomite and fibrous calcite

Authigenic clays are an important control on reservoir quality in lacustrine carbonates but remain challenging to predict. Lacustrine depositional systems respond to climatic variations in rainfall, surface runoff and groundwater input, and evaporation, and result in rapid and frequent changes in lake volume; this is expressed through changing water depth and shoreline position. In the upper portion of the Early Palaeocene Yacoraite Formation of the Salta Basin in Argentina, extensive lacustrine deposits were deposited during the sag phase of rifting. Prior high-resolution stratigraphic studies have suggested that climatic factors control microbial carbonate sequences within a ‘balanced fill’ lake, with variation in the lake level having a major influence on facies association changes. This study characterizes the evolution of facies and mineralogy within the Yacoraite Formation, focusing on the distribution of clay minerals, making a link between the high, medium and low-frequency sequence stratigraphic cycles. The low-frequency transgressive hemicycle of the upper portion of the Yacoraite Formation is comprised of abundant siliciclastic facies, suggesting a wetter period. Microbialites occurring in this interval are coarse-grained and agglutinated. Detrital clay minerals such as illite and chlorite and associated siliciclastic sediments were input to the lake during high-frequency transgressive periods. During high-frequency regressive hemicycles, sedimentation was dominated by carbonate facies with Ca-rich dolomite and the authigenic clays are comprised of chlorite/smectite mixed-layers. By contrast, the low frequency regressive hemicycle records fine-grained agglutinated microbialite with horizons of fibrous calcite, more stoichiometric dolomite, barite and authigenic magnesian smectite. This indicates elevated ion concentrations in the lake under intense evaporation during an arid period. Understanding the conditions that are favourable for formation and preservation of authigenic clays within the lacustrine environment can improve understanding of reservoir quality in comparable economically important deposits.     

Separation of PP- and PS-wave reflected seismic data using two-dimensional finite offset common-reflection-surface traveltime approximation

Recently, the interest in PS-converted waves has increased for several applications, such as sub-basalt layer imaging, impedance estimates and amplitude-versus-offset analysis. In this study, we consider the problem of separation of PP- and PS-waves from pre-stacked multicomponent seismic data in two-dimensional isotropic medium. We aim to demonstrate that the finite-offset common-reflection-surface traveltime approximation is a good alternative for separating PP- and PS-converted waves in common-offset and common shot configurations by considering a two-dimensional isotropic medium. The five parameters of the finite-offset common-reflection-surface are firstly estimated through the inversion methodology called very fast simulated annealing, which estimates all parameters simultaneously. Next, the emergence angle, one of the inverted parameters, is used to build an analytical separation function of PP and PS reflection separation based on the wave polarization equations. Once the PP- and PS-converted waves were separated, the sections are stacked to increase the signal-to-noise ratio using the special curves derived from finite-offset common-reflection-surface approximation. We applied this methodology to a synthetic dataset from simple-layered to complex-structured media. The numerical results showed that the inverted parameters of the finite offset common-reflection-surface and the separation function yield good results for separating PP- and PS-converted waves in noisy common-offset and common shot gathers.     

Estimating free parameters in 2D inversion: example of gravity inversion in a rifted basement

We advance a principle directed to assigning numerical values to free parameters usually present in inversion methods. It may be formulated as: ‘Optimum estimates of free parameters in an inversion procedure must lead, in tests using synthetic data, to solutions whose geometrical expression reflects the main qualitative or semiquantitative geological characteristic of the study area.’ To this end, the interpreter should (i) specify a typical anomalous source geometry which incorporates the most relevant geological information for the study area, (ii) compute the corresponding gravity anomaly and (iii) invert the anomaly for the source geometry finding the numerical values of the free parameters that lead to a solution closest to the typical source. Application of the above methodology to synthetic and real data from the basement relief of a rift basin has asserted its efficacy.     

Toward the understanding of the treatment of textile industries’ effluents by clay: adsorption of anionic dye on kaolinite

Clays, particularly kaolinite, are promising adsorbents for the treatment of textile effluents, but there is a need of better understanding the mechanisms of adsorption, especially in the case of anionic dyes. Thus, the removal of RR120 anionic dye was investigated using Tunisian raw clay (TBK) composed of kaolinite and illite, and a standard kaolinite (KGa-2), and conducting batch experiments by varying different parameters (contact time, ionic strength, concentration, temperature). We investigated the clays’ surface charges by electrophoretic mobility measures and the dye-clay interactions during adsorption, by the streaming-induced potentials (SIP). The results showed that KGa-2 has higher adsorption capacity for RR120 dye than TBK clay, moreover enhanced by increasing the ionic strength and/or lowering the pH of the aqueous. The SIP results showed an increase of negative charges for both clays, reflecting the adsorption of the anionic dye on the positive charges of the amphoteric surfaces of the clays. The SIP magnitudes indicated a higher adsorption rate for KGa-2 in accordance with the kinetic study. The Sips model that described the best adsorption isotherms indicates lateral interactions of the dye molecules, stronger in the case of KGa-2 than TBK. Also, the dye molecules form a thinner layer on KGa-2 surfaces. In addition, the dye molecule’s structure was not altered, as verified by mass spectrometry. The adsorption process was feasible and spontaneous and favored at ambient temperature. Thus, kaolinite-rich clays are effective in the removal of anionic dyes in aqueous solution and potential good adsorbents in wastewater treatment.     

Ultraviolet spectroscopy of Asteroid (4) Vesta

We report a comprehensive review of the UV–visible spectrum and rotational lightcurve of Vesta combining new observations by Hubble Space Telescope and Swift Gamma-ray Burst Observatory with archival International Ultraviolet Explorer observations. The geometric albedos of Vesta from 220 nm to 953 nm are derived by carefully comparing these observations from various instruments at different times and observing geometries. Vesta has a rotationally averaged geometric albedo of 0.09 at 250 nm, 0.14 at 300 nm, 0.26 at 373 nm, 0.38 at 673 nm, and 0.30 at 950 nm. The linear spectral slope as measured between 240 and 320 nm in the ultraviolet displays a sharp minimum near a sub-Earth longitude of 20°, and maximum in the eastern hemisphere. This is consistent with the longitudinal distribution of the spectral slope in the visible wavelength. The photometric uncertainty in the ultraviolet is ∼20%, and in the visible wavelengths it is better than 10%. The amplitude of Vesta’s rotational lightcurves is ∼10% throughout the range of wavelengths we observed, but is smaller at 950 nm (∼6%) near the 1-μm band center. Contrary to earlier reports, we found no evidence for any difference between the phasing of the ultraviolet and visible/near-infrared lightcurves with respect to sub-Earth longitude. Vesta’s average spectrum between 220 and 950 nm can well be described by measured reflectance spectra of fine particle howardite-like materials of basaltic achondrite meteorites. Combining this with the in-phase behavior of the ultraviolet, visible, and near-infrared lightcurves, and the spectral slopes with respect to the rotational phase, we conclude that there is no global ultraviolet/visible reversal on Vesta. Consequently, this implies a lack of global space weathering on Vesta, as previously inferred from visible–near-infrared data.