Analysis of DNAPL infiltration in a medium with a low-permeable lens

In this paper we study the infiltration of DNAPL in a porous medium containing a single low-permeable lens. Our aim is to determine whether or not DNAPL infiltrates into the lens. A key role is played by the capillary pressure: DNAPL cannot infiltrate into the lens unless the capillary pressure exceeds the entry pressure of the lens. In the model this is reflected by an interface condition, the extended capillary pressure condition. To derive analytical approximations we first consider a steady-state DNAPL plume in a homogeneous medium. This results in an estimate of the DNAPL plume width as a function of depth, and an asymptotic solution for small saturations. Assuming that the extent of the lens is much larger than the width of the unperturbed DNAPL plume in the homogeneous medium, we derive an explicit criterion for DNAPL infiltration into the lens in terms of a critical inflow rate. A numerical algorithm is presented in which the extended capillary pressure condition is incorporated. The numerical and analytical results show good qualitative agreement. This revised version was published online in August 2006 with corrections to the Cover Date.     

A general correspondence principle for time-domain electromagnetic wave and diffusion fields

A general correspondence principle is presented that relates any time-domain electromagnetic diffusion field to an electromagnetic wavefield in a 'corresponding' configuration. The principle applies to arbitrarily inhomogeneous and anisotropic media and arbitrary transmitters and receivers. For the correspondence between the two types of electromagnetic fields to hold, the electric conductivity in the diffusive case and the permittivity in the wavefield case should have the same spatial variation, while the permeability distributions in space in the two cases are to be identical. Essential steps in the derivation of the correspondence principle are the use of the time Laplace transformation of causal signals, taken at real, positive values of the transform parameter, the Schouten-Van der Pol theorem in the theory of the Laplace transformation, and the reliance upon Lerch's theorem of the uniqueness of the interrelation between causal field quantities and their time-Laplace-transform representations at real, positive values of the transform parameter. Correspondence is then established between the tensorial Green's functions in the two cases, where the Green's functions are the point-receiver responses (either electric or magnetic field) to point-transmitter excitations (either electric- or magnetic-current source).
Through the correspondence principle, all transient electromagnetic wavefields (where losses are neglected) have as a counterpart a transient diffusive electromagnetic field (where the electric displacement current is neglected). The interrelation yields the tool to compare quantitatively the potentialities of the two types of fields in transient electromagnetic geophysical prospecting.
Finally, a general medium-parameter scaling law for time-domain electromagnetic wavefields is presented.     

Transport in a 2-D saturated porous medium: A new method for particle tracking

A new method for solving the transport equation based on the management of a large numbe of particles in a discretized 2-D domain is presented. The method uses numerical variables to represent the number of particles in a given mesh and is more complex than the 1-D problem. The first part of the paper focuses on the specific management of particles in a 2-D problem. The method also would be valid for three dimensions as long as the medium can be modeled similar to a layered system. As the particles are no longer tracked individually, the algorithm is fast and does not depend on the number of particles present. The numerical tests show that the method is nearly numerical dispersion free and permits accurate calculations even for simulations of low-concentration transport. Because each mesh is considered as a closed system between two successive time steps, it is easy to add adsorption phenomenon without any problem of numerical stability. The model is tested under conditions that are extremely demanding for its operating mode and gives a good fit to analytical solutions. The conditions in which it can be used to best advantage are discussed.     

Early Archean processes: evidence from the South African Kaapvaal craton and its greenstone belts (extended abstract)

Geologie en Mijnbouw -     

The ClBrI composition of 3.23 Ga modified seawater: implications for the geological evolution of ocean halide chemistry

Fluid inclusion leachates obtained from vug and vein quartz samples from an Archean (3.23 Ga) Fe-oxide hydrothermal deposit in the west-central part of the Barberton greenstone belt, South Africa, were analyzed by ion chromatography for chloride, bromide, and iodide. The deposit, known as the ironstone pods, formed by seafloor hydrothermal activity and fluid discharge. Quartz is dominated by type I liquid-vapor, aqueous inclusions with a bimodal salinity distribution (0–0.25 M_Cl⁻ and 0.9–1.8 M_Cl⁻). Bulk analytical salinities range from 0.45 to 0.99 M_Cl⁻ represent averages of type I inclusions. Bulk fluid inclusion bromide and iodide concentrations are 1.44–3.32 mM and 0.01–0.12 mM, respectively. For comparison, modern seawater has halogen contents of 590 mM chloride, 0.9 mM bromide, and 0.5 μM total iodine. In the fluids from the ironstone pods, bromide and iodide are enriched relative to chloride, when compared with modern seawater.Approximate Br⁻Cl⁻ and I⁻Cl⁻ ratios of 3.2 Ga Barberton seawater are 2.5 × 10⁻³ and 40 × 10⁻⁶, respectively. Dispersion to higher values was caused principally by reaction with organic sediments whose trends are similar to those seen for modern vent fluids at unsedimented and sedimented ridges, relative to modern seawater. These halide ratios are greater than those of modern seawater, suggesting a change in the halide ratios of seawater over geological time. The analytical data are consistent with a model in which marine organic sedimentation has fractionated bromine and iodine out of seawater relative to chloride, thereby causing the halide ratios of seawater to decrease from high early and mid-Archean values towards their present day values.     

Assessment of predictions of the response of the Hualien containment model during forced vibration tests

Blind predictions for the response of the 1/4-scale reinforced concrete Hualien (Taiwan) containment model during forced vibration tests are compared with the observed data. The predictions obtained by the CLASSI approach reflect the experimental conditions prior to and after backfill of the soil surrounding the embedded foundation. The experimental data show a strong and unexpected coupling between the response in the NS and EW directions which is not present in the results for the axisymmetric theoretical models. Also, significant differences can be seen between the experimental responses in the two orthogonal horizontal directions which minimize cross-axis coupling. Although these differences are not accounted for in the theoretical models, the discrepancies between predictions and observations are within the uncertainty of the structural and geotechnical data. The obtained differences between predictions and observations give an excellent measure of the prediction error that can be expected in this type of analysis from uncertainty in the data. A detailed assessment of the initial structural and geotechnical data based on extensive comparisons with the results of previous identification studies is also presented. Finally, comparisons between the observed response and calculations based on revised models for the structure and the soil show that current methods of analysis can account accurately for the observed response.     

Neuralizing target superresolution algorithms

Tatem et al. (2001) have designed a Hopfield network-based algorithm for superresolving discrete targets that are larger than the sample spacing of an image. The algorithm iteratively minimizes a criterion function that contains a sigmoidal activation term. We have altered their algorithm to bring it in line with Hopfield's original network by reducing the pseudotemperature of the sigmoid. We found that smaller values of the pseudotemperature lead to faster convergence to a solution and resulting solutions that are more accurate.     

Beryllium-10 in Australasian tektites: Constraints on the location of the source crater

By using accelerator mass spectrometry we have measured the ¹⁰Be concentrations of 86 Australasian tektites. Corrected to the time of tektite production ∼0.8 My ago, the ¹⁰Be concentrations (10⁶ atom/g) range from 59 for a layered tektite from Huai Sai, Thailand, to 280 for an australite from New South Wales, Australia. The average value is 143 ± 50. When tektites are sorted by country, their average measured ¹⁰Be concentrations increase slowly with increasing distance from Southeast Asia, the probable location of the tektite producing event, from 59 ± 9 for 6 layered tektites from Laos to 136 ± 20 for 20 splash-form tektites from Australia. The lowest ¹⁰Be concentrations for tektites fall on or within a contour centered off the shore of Vietnam, south of the Gulf of Tonkin (107°E; 17°N), but also encompassing two other locations in the area of northeastern Thailand previously proposed for the site of a single tektite-producing impact. The ¹⁰Be concentrations of layered tektites show only a weak anticorrelation (R ∼ −0.3) with the numbers of relict crystalline inclusions.Loosely consolidated, fine-grained terrestrial sediments or recently consolidated sedimentary rocks are the most likely precursor materials. Dilution of sediments with other kinds of rock raises problems in mixing and is not supported by petrographic data. Sedimentary columns that have the right range of ¹⁰Be concentrations occur off the coasts of places where sedimentation rates are high relative to those in the deep sea. A single impact into such a region, 15 to 300 m thick, could have propelled near-surface, high-¹⁰Be material farthest—to Australia—while keeping the deeper-lying, low-¹⁰Be layers closer to home. We do not rule out, however, other proposed mechanisms for tektite formation.