Prediction of rock alteration patterns: a potential tool in mineral exploration

Hydrothermal alteration of a quartz-K-feldspar rock is simulated numerically by coupling fluid flow and chemical reactions. Introduction of CO2 gas generates an acidic fluid and produces secondary quartz, muscovite and/or pyrophyllite at constant temperature and pressure of 300 degrees C and 200 MPa. The precipitation and/or dissolution of the secondary minerals is controlled by either mass-action relations or rate laws. In our simulations the mass of the primary elements are conserved and the mass-balance equations are solved sequentially using an implicit scheme in a finite-element code. The pore-fluid velocity is assumed to be constant. The change of rock volume due to the dissolution or precipitation of the minerals, which is directly related to their molar volume, is taken into account. Feedback into the rock porosity and the reaction rates is included in the model. The model produces zones of pyrophyllite quartz and muscovite due to the dissolution of K-feldspar. Our model simulates, in a simplified way, the acid-induced alteration assemblages observed in various guises in many significant mineral deposits. The particular aluminosilicate minerals produced in these experiments are associated with the gold deposits of the Witwatersrand Basin.

Exposure of rats to a high but not low dose of ethanol during early postnatal life increases the rate of loss of optic nerve axons and decreases the rate of myelination

Visual system abnormalities are commonly encountered in the fetal alcohol syndrome although the level of exposure at which they become manifest is uncertain. In this study we have examined the effects of either low (ETLD) or high dose (ETHD) ethanol, given between postnatal days 4-9, on the axons of the rat optic nerve. Rats were exposed to ethanol vapour in a special chamber for a period of 3 h per day during the treatment period. The blood alcohol concentration in the ETLD animals averaged similar to 171 mg/dl and in the ETHD animals similar to 430 mg/dl at the end of the treatment on any given day. Groups of 10 and 30-d-old mother-reared control (MRC), separation control (SC), ETLD and ETHD rats were anaesthetised with an intraperitoneal injection or ketamine and xylazine, and killed by intracardiac perfusion with phosphate-buffered glutaraldehyde. In the 10-d-old rat optic nerves there was a total of similar to 145000-165000 axons in MRC, SC and ETLD animals. About 4 % of these fibres were myelinated. The differences between these groups were not statistically significant. However, the 10-d-old ETHD animals had only about 75000 optic nerve axone (P < 0.05) of which about 2.8 % were myelinated. By 30 d of age there was a total of between 75000 90000 optic nerve axons, irrespective of the group examined. The proportion of axons which were myelinated at this age was still significantly lower (P < 0.001) in the ETHD animals (similar to 77 %) than in the other groups (about 98 %). It is concluded that the normal stages of development and maturation of the rat optic nerve axons, as assessed in this study, can be severely compromised by exposure to a relatively high (but not low) dose of ethanol between postnatal d 4 and 9.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Identifying rate-limiting nodes in large-scale cortical networks for visuospatial processing: an illustration using fMRI

With the advent of functional neuroimaging techniques, in particular functional magnetic resonance imaging (fMRI), we have gained greater insight into the neural correlates of visuospatial function. However, it may not always be easy to identify the cerebral regions most specifically associated with performance on a given task. One approach is to examine the quantitative relationships between regional activation and behavioral performance measures. In the present study, we investigated the functional neuroanatomy of two different visuospatial processing tasks, judgement of line orientation and mental rotation. Twenty-four normal participants were scanned with fMRI using blocked periodic designs for experimental task presentation. Accuracy and reaction time (RT) to each trial of both activation and baseline conditions in each experiment was recorded. Both experiments activated dorsal and ventral visual cortical areas as well as dorsolateral prefrontal cortex. More regionally specific associations with task performance were identified by estimating the association between (sinusoidal) power of functional response and mean RT to the activation condition; a permutation test based on spatial statistics was used for inference. There was significant behavioral-physiological association in right ventral extrastriate cortex for the line orientation task and in bilateral (predominantly right) superior parietal lobule for the mental rotation task. Comparable associations were not found between power of response and RT to the baseline conditions of the tasks. These data suggest that one region in a neurocognitive network may be most strongly associated with behavioral performance and this may be regarded as the computationally least efficient or rate-limiting node of the network.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

MHC class II expression is regulated in dendritic cells independently of invariant chain degradation

We have investigated the mechanisms that control MHC class II (MHC II) expression in immature and activated dendritic cells (DC) grown from spleen and bone marrow precursors. Degradation of the MHC II chaperone invariant chain (li), acquisition of peptide cargo by MHC II, and delivery of MHC II-peptide complexes to the cell surface proceeded similarly in both immature and activated DC. However, immature DC reendocytosed and then degraded the MHC II-peptide complexes much faster than the activated DC. MHC II expression in DC is therefore not controlled by the activity of the protease(s) that degrade Ii, but by the rate of endocytosis of peptide-loaded MHC II. Late after activation, DC downregulated MHC II synthesis both in vitro and in vivo.

Change in the rate of shell deposition and shell microstructure in response to shell borers in the abalone Haliotis rubra

Gastropod shells consist of two crystal types of calcium carbonate, an outer, prismatic calcite layer and an inner nacreous layer made of aragonite. In cross-section, the nacre of the nacreous layer appears to have a regular brick-like microstructure composed of thin laminae of aragonite crystals, separated by very thin sheets of protein (Lutz and Rhoads, 1980; Nakahara, 1983). In abalone (Genus, Haliotis) and other gastropods, thin layers of non-lamellar pigmented material occur within the nacre and have been termed alternatively, fine lines, growth rings or growth lines (Shepherd et al., 1995). It has been suggested that these pigmented layers are small, prismatic, calcite layers (Shepherd and Avalos-Borja, 1997; Zaremba et al., 1996) but investigations using a Raman laser in Haliotis rubra show that they contain aragonite rather than calcite (Hawkes et al, 1996). Day and Fleming (1992) suggest that the occurrence of pigmented layers is correlated with regular exogenous cues such as reproduction or temperature changes and indeed in some species, pigmented layers in the shell can be used to age abalone (review: Shepherd and Triantafillos, 1997). However, McShane and Smith (1992) suggest that pigmented layers can occur irregularly and therefore may be unreliable indicators of age.

Similitude applied to centrifugal scaling of unsaturated flow

Centrifuge experiments modeling single-phase flow in prototype porous media typically use the same porous medium and permeant. Then, well-known scaling laws are used to transfer the results to the prototype. More general scaling laws that relax these restrictions are presented. For permeants that are immiscible with an accompanying gas phase, model-prototype (i.e., centrifuge model experiment-target system) scaling is demonstrated. Scaling is shown to be feasible for Miller-similar (or geometrically similar) media. Scalings are presented for a more, general class, Lisle-similar media, based on the equivalence mapping of Richards' equation onto itself. Whereas model-prototype scaling of Miller-similar media can be realized easily for arbitrary boundary conditions, Lisle-similarity in a finite length medium generally, but not always, involves a mapping to a moving boundary problem. An exception occurs for redistribution in Lisle-similar porous media, which is shown to map to spatially fixed boundary conditions. Complete model-prototype scalings for this example are derived.

Mathematical models in percutaneous absorption

A number of mathematical models have been used to describe percutaneous absorption kinetics. In general, most of these models have used either diffusion-based or compartmental equations. The object of any mathematical model is to a) be able to represent the processes associated with absorption accurately, b) be able to describe/summarize experimental data with parametric equations or moments, and c) predict kinetics under varying conditions. However, in describing the processes involved, some developed models often suffer from being of too complex a form to be practically useful. In this chapter, we attempt to approach the issue of mathematical modeling in percutaneous absorption from four perspectives. These are to a) describe simple practical models, b) provide an overview of the more complex models, c) summarize some of the more important/useful models used to date, and d) examine sonic practical applications of the models. The range of processes involved in percutaneous absorption and considered in developing the mathematical models in this chapter is shown in Fig. 1. We initially address in vitro skin diffusion models and consider a) constant donor concentration and receptor conditions, b) the corresponding flux, donor, skin, and receptor amount-time profiles for solutions, and c) amount- and flux-time profiles when the donor phase is removed. More complex issues, such as finite-volume donor phase, finite-volume receptor phase, the presence of an efflux. rate constant at the membrane-receptor interphase, and two-layer diffusion, are then considered. We then look at specific models and issues concerned with a) release from topical products, b) use of compartmental models as alternatives to diffusion models, c) concentration-dependent absorption, d) modeling of skin metabolism, e) role of solute-skin-vehicle interactions, f) effects of vehicle loss, a) shunt transport, and h) in vivo diffusion, compartmental, physiological, and deconvolution models. We conclude by examining topics such as a) deep tissue penetration, b) pharmacodynamics, c) iontophoresis, d) sonophoresis, and e) pitfalls in modeling.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

BeachWin: modelling groundwater effects on swash sediment transport and beach profile changes

Field and laboratory observations have shown that a relatively low beach groundwater table enhances beach accretion. These observations have led to the beach dewatering technique (artificially lowering the beach water table) for combating beach erosion. Here we present a process-based numerical model that simulates the interacting wave motion on the beach. coastal groundwater flow, swash sediment transport and beach profile changes. Results of model simulations demonstrate that the model replicates accretionary effects of a low beach water table on beach profile changes and has the potential to become a tool for assessing the effectiveness of beach dewatering systems. (C) 2002 Elsevier Science Ltd. All rights reserved.

Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

Properties of the vertebrate skeletal muscle tubular system as a sealed compartment

Confocal imaging of impermeant fluorescent dyes trapped in the tubular (t-) system of skeletal muscle fibres of rat and cane toad was used to examine changes in the morphology of the t-system upon mechanical skinning, the time course of dye loss from the sealed t-systern in mechanically skinned fibres and the influence of rapid application and removal of glycerol on the morphology of the sealed t-system. In contrast to intact fibres, which have a t-systern open to the outside, the sealed t-systern of toad mechanically skinned fibres consistently displayed local swellings (vesicles). The occurrence of vesicles in the sealed t-system of rat-skinned fibres was infrequent. Application and removal of 200-400 mM glycerol to the sealed t-system did not produce any obvious changes in its morphology. The dyes fluo-3, fura-2 and Oregon green 488 were lost from the sealed t-system of toad fibres at different rates suggesting that the mechanism of organic anion transport across the tubular wall was not by indiscriminate bulk transport. The rate of fluo-3 and fura-2 loss from the sealed t-system of rat fibres was greater in rat than in toad fibres and could be explained by differences in surface area: volume ratio of the t-system in the two fibre types. Based on the results presented here and on other results from this laboratory, an explanation is given for the formation of numerous vesicles in toad-skinned fibres and lack of vesicle formation in rat-skinned fibres. This explanation can also help with better understanding the mechanism responsible for vacuole formation in intact fibres. (C) 2002 Elsevier Science Ltd. All rights reserved.