Numerical and experimental study of seepage in unconfined aquifers with a periodic boundary condition

The assessment of groundwater conditions within an unconfined aquifer with a periodic boundary condition is of interest in many hydrological and environmental problems. A two-dimensional numerical model for density dependent variably saturated groundwater flow, SUTRA (Voss, C.I., 1984. SUTRA: a finite element simulation model for saturated-unsaturated, fluid-density dependent ground-water flow with energy transport or chemically reactive single species solute transport. US Geological Survey, National Center, Reston, VA) is modified in order to be able to simulate the groundwater flow in unconfined aquifers affected by a periodic boundary condition. The basic flow equation is changed from pressure-form to mixed-form. The model is also adjusted to handle a seepage-face boundary condition. Experiments are conducted to provide data for the groundwater response to the periodic boundary condition for aquifers with both vertical and sloping faces. The performance of the numerical model is assessed using those data. The results of pressure- and mixed-form approximations are compared and the improvement achieved through the mixed-form of the equation is demonstrated. The ability of the numerical model to simulate the water table and seepage-face is tested by modelling some published experimental data. Finally the numerical model is successfully verified against present experimental results to confirm its ability to simulate complex boundary conditions like the periodic head and the seepage-face boundary condition on the sloping face. (C) 1999 Elsevier Science B.V. All rights reserved.

Providing evidence-based answers to clinical questions - A pilot information service for general practitioners

Objective: To pilot a clinical information service for general practitioners. Methods: A representative sample of 31 GPs was invited to submit clinical questions to a local academic department of general practice. Their views on the service and the usefulness of the information were obtained by telephone interview. Results: Over one month, nine GPs (29% of the sample, 45% of those stating an interest), submitted 20 enquiries comprising 45 discrete clinical questions. The median time to search for evidence, appraise it and write answers to each enquiry was 2.5 hours (range, 1.0-7.4 hours). The median interval between receipt of questions and dispatch of answers was 3 clays (range, 1-12 days). Conclusions: The GPs found the answers useful in clinical decision making; in four out of 20 cases patient management was altered.

Increasing the efficiency of solute leaching: impacts of flow interruption with drainage of the "preferential flow paths"

Most soils contain preferential flow paths that can impact on solute mobility. Solutes can move rapidly down the preferential flow paths with high pore-water velocities, but can be held in the less permeable region of the soil matrix with low pore-water velocities, thereby reducing the efficiency of leaching. In this study, we conducted leaching experiments with interruption of the flow and drainage of the main flow paths to assess the efficiency of this type of leaching. We compared our experimental results to a simple analytical model, which predicts the influence of the variations in concentration gradients within a single spherical aggregate (SSA) surrounded by preferential flow paths on leaching. We used large (length: 300 mm, diameter: 216 mm) undisturbed field soil cores from two contrasting soil types. To carry out intermittent leaching experiments, the field soil cores were first saturated with tracer solution (CaBr2), and background solution (CaCl2) was applied to mimic a leaching event. The cores were then drained at 25- to 30-cm suction to empty the main flow paths to mimic a dry period during which solutes could redistribute within the undrained region. We also conducted continuous leaching experiments to assess the impact of the dry periods on the efficiency of leaching. The flow interruptions with drainage enhanced leaching by 10-20% for our soils, which was consistent with the model's prediction, given an optimised equivalent aggregate radius for each soil. This parameter quantifies the time scales that characterise diffusion within the undrained region of the soil, and allows us to calculate the duration of the leaching events and interruption periods that would lead to more efficient leaching. Application of these methodologies will aid development of strategies for improving management of chemicals in soils, needed in managing salts in soils, in improving fertiliser efficiency, and in reclaiming contaminated soils. (C) 2000 Elsevier Science B.V. All rights reserved.

Petrographic and isotope constraints on the origin of authigenic carbonate minerals and the associated fluid evolution in Late Permian coal measures, Bowen Basin (Queensland), Australia

Authigenic carbonate minerals are ubiquitous throughout the Late Permian coal measures of the Bowen Basin, Queensland, Australia. In the northern Bowen Basin, carbonates include the following assemblages: siderite I (delta O-18(SMOW) = +11.4 to + 17%, delta C-13(PDB) = - 5.3 to + 120), Fe-Mg calcite-ankerite-siderite II mineral association (delta O-18(SMOW) = +7.2 to + 10.20, delta C-13(PDB) = 10.9 to - 1.80 for ankerite) and a later calcite (delta O-18(SMOW) = +5.9 to + 14.60, delta C-13(PDB) = -11.4 to + 4.40). In the southern Bowen Basin, the carbonate phase consists only of calcite (delta O-18(SMOW) = +12.5 to + 14.80, delta C-13(PDB) = -19.4 to + 0.80), where it occurs extensively throughout all stratigraphic levels. Siderite I occurs in mudrocks and sandstones and predates all other carbonate minerals. This carbonate phase is interpreted to have formed as an early diagenetic mineral from meteoric waters under cold climate and reducing conditions. Fe-Mg calcite-ankerite-siderite Il occur in sandstones as replacement of volcanic rock fragments. Clay minerals (illite-smectite, chlorite and kaolinite) postdate Ca-Fe-Mg carbonates, and precipitation of the later calcite is associated with clay mineral formation. The Ca-Fe-Mg carbonates and later calcite of the northern Bowen Basin are regarded as having formed as a result of hydrothermal activity during the latest Triassic extensional tectonic event which affected this part of the basin, rather than deep burial diagenesis during the Middle to Late Triassic as previously reported. This hypothesis is based on the timing relationships of the authigenic mineral phases and the low delta O-18 values of ankerite and calcite, together with radiometric dating of illitic clays and recently published regional geological evidence. Following the precipitation of the Ca-Fe-Mg carbonates from strongly O-18-depleted meteoric-hydrothermal fluids, continuing fluid circulation and water-rock interaction resulted in dissolution of these carbonate phases as well as labile fragments of volcaniclastic rocks. Subsequently, the later calcite and day minerals precipitated from relatively evolved (O-18-enriched) fluids. The nearly uniform delta O-18 values of the southern Bowen Basin calcite have been attributed to very low water/rock ratio in the system, where the fluid isotropic composition was buffered by the delta O-18 values of rocks. (C) 2000 Elsevier Science B.V. All rights reserved.

Sudden drawdown and drainage of a horizontal aquifer

Drainage of a saturated horizontal aquifer following a sudden drawdown is reanalyzed using the Boussinesq equation. The effect of the finite length of the aquifer is considered in detail. An analytical approximation based on a superposition principle yields a very good estimate of the outflow when compared to accurate numerical solutions. An illustration of the new analytical approach to analyze basin-scale field data is used to demonstrate possible field applications of the new solution.

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.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

Computer simulations of coupled problems in geological and geochemical systems

The finite element method is used to simulate coupled problems, which describe the related physical and chemical processes of ore body formation and mineralization, in geological and geochemical systems. The main purpose of this paper is to illustrate some simulation results for different types of modelling problems in pore-fluid saturated rock masses. The aims of the simulation results presented in this paper are: (1) getting a better understanding of the processes and mechanisms of ore body formation and mineralization in the upper crust of the Earth; (2) demonstrating the usefulness and applicability of the finite element method in dealing with a wide range of coupled problems in geological and geochemical systems; (3) qualitatively establishing a set of showcase problems, against which any numerical method and computer package can be reasonably validated. (C) 2002 Published by Elsevier Science B.V.

Antioxidant behaviour of thia fatty acids

Eight thia fatty acids and other sulfides have been studied as inhibitors of autoxidation of arachidonic acid. The inhibitors extend the lag phase of the oxidation, to varying degrees. A carboxyl group in the vicinity of the sulfur reduces the antioxidant activity, while unsaturated sulfides are more effective than their saturated analogues. The results are consistent with the sulfides acting to reduce fatty acid hydroperoxides, which otherwise accumulate during the early stages of reaction and propagate the free-radical oxidation process.

Visible and near-infrared luminescent Eu(3+) or Er(3+) doped laponite-derived xerogels and thick films: Structural and spectroscopic properties

Laponite-derived materials represent promising materials for optical applications. In this work, Eu(3+)- or Er(3+)-doped laponite xerogels and films were prepared from colloidal dispersion. Homogeneous, crack-free and transparent single layers were deposited on soda-lime substrates with a thickness of 10 mu m. Structural and spectroscopic properties were analyzed by thermal analyses, X-ray diffractometry, transmission electron microscopy, infrared spectroscopy, and luminescence spectroscopy. The addition of a rare earth ion to the laponite does not promote any changes in thermal stability or phase transition. Laponite clay was identified after annealing up to 500 degrees C, with a decrease in basal spacing when the annealing temperature is changed from 100 degrees C to 500 degrees C. Enstatite polymorphs and amorphous silicate phases were observed after heat treatment at 700 degrees C and 900 degrees C. Stationary and time-dependent luminescence spectra in the visible region for Eu(3+), and (5)D(0) lifetime are discussed in terms of thermal treatment and structural evolution. In the layered host, the Eu(3+) ions are distributed in many different local environments. However, Eu(3+) ions were found to occupy at least two symmetry sites, and the ions are preferentially incorporated into the crystalline enstatite for the materials annealed at 700 degrees C and 900 degrees C. A (5)D(0) lifetime of 1.3 ms and 3.1 ms was obtained for Eu(3+) ions in an amorphous silicate and crystalline MgSiO(3) local environment, respectively. Strong Er(3+) emission at the 1550 nm region was observed for the materials annealed at 900 degrees C, with a bandwidth of 44 nm. (C) 2008 Elsevier B.V. All rights reserved.

Interfacial interactions and structure of organic-inorganic nanohybrids

Understanding the interfacial interactions and structure is important to better design and application of organic-inorganic nanohybrids. This paper presents our recent molecular dynamic studies on organoclays and polymer nanocomposites, including the layering behavior of organoclays, structural and dynamic properties of dioctadecyldimethyl ammoniums in organoclays, and interfacial interactions and structure of polyurethane nanocomposites. The results demonstrate that the layering behaviors of organoclays are closely related to the chain length of quaternary alkyl ammoniums and cation exchangeable capacity of clays. In addition to typical layered structures such as monolayer, bilayer and pseudo-trilayer, a pseudo-quadrilayer structure was also observed in organoclays modified with dioctadecyldimethyl ammoniums (DODDMA). In such a structure, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer or even the next nearest layer. Moreover, the diffusion constants of nitrogen and methylene atoms increase with the temperature and methelene towards the tail groups. For polyurethane nanocomposite, the van der Waals interaction between apolar alkyl chains and soft segments of polyurethane predominates the interactions between organoclay and polyurethane. Different from most bulk polyurethane systems, there is no distinct phase-separated structure for the polyurethane.

MINERALOGY AND GENESIS OF SMECTITES IN AN ALKALINE-SALINE ENVIRONMENT OF PANTANAL WETLAND, BRAZIL

Smectite formation in alkaline-saline environments has been attributed to direct precipitation from solution and/or transformation from precursor minerals, but these mechanisms are not universally agreed upon in the literature. The objective of this work was to investigate the mineralogy of smectites in the soils surrounding a representative alkaline-saline lake of Nhecolandia, a sub-region of the Pantanal wetland, Brazil, and then to identify the mechanisms of their formation. Soils were sampled along a toposequence and analyzed by X-ray diffraction, transmission electron microscopy-energy dispersive X-ray analysis, and inductively coupled plasma-mass spectrometry. Water was collected along a transect involving the studied toposequence and equilibrium diagrams were calculated using the databases PHREEQC and AQUA. The fine-clay fraction is dominated by smectite, mica, and kaolinite. Smectites are concentrated at two places in the toposequence: an upper zone, which includes the soil horizons rarely reached by the lake-level variation; and a lower zone, which includes the surface horizon within the area of seasonal lake-level variation. Within the upper zone, the smectite is dioctahedral, rich in Al and Fe, and is classified as ferribeidellite. This phase is interstratified with mica and vermiculite and has an Fe content similar to that of the mica identified. These characteristics suggest that the ferribeidellite originates from transformation of micas and that vermiculite is an intermediate phase in this transformation. Within the lower zone, smectites are dominantly trioctahedral, Mg-rich, and are saponitic and stevensitic minerals. In addition, samples enriched in these minerals have much smaller rare-earth element (REE) contents than other soil samples. The water chemistry shows a geochemical control of Mg and saturation with respect to Mg-smectites in the more saline waters. The REE contents, water chemistry, and the presence of Mg-smectite where maximum evaporation is expected, suggest that saponitic and stevensitic minerals originate by chemical precipitation from the water column of the alkaline-saline lake.

A fructan: fructan fructosyltransferase activity from Lolium rigidum

Fructan:fructan fructosyltransferase (FFT) activity was purified about 300-fold from leaves of Lolium rigidura Gaudin by a combination of affinity chromatography, gel filtration, anion exchange and isoelectric focusing. The FFT activity was free of sucrose:sucrose fructosyltransferase and invertase activities. It had an apparent pI of 4.7 as determined by isoelectric focusing, and a molecular mass of about 50000 (gel filtration). The FFT activity utilized the trisaccharides 1-kestose and 6(G)-kestose as sole substrates, but was not able to use 6-kestose as sole substrate. The FFT activity was not saturated when assayed at concentrations of 1-kestose, 6(G)-kestose or (1,1)-kestotetraose of up to 400 mM The rate of reaction of the FFT activity was most rapid when assayed with 1-kestose and was less rapid when assayed with 6(G)-kestose, (1,1)-kestotetraose or (1,1,1)-kestopentaose. The FFT activity when assayed at a relatively high concentration of enzyme activity (approximately equivalent to about half the activity in crude extracts per gram fresh mass) did not synthesize fructan of degree of polymerization > 6, even during extended assays of up to 10 h. When assayed with a combination of 1-kestose and uniformly labelled [C-14]sucrose as substrates, the major reaction was the transfer of a fructosyl residue from 1-kestose to sucrose resulting in the re-synthesis of 1-kestose. Tetrasaccharide and 6(G)-kestose were also synthesized. When assayed with 6(G)-kestose and [C-14]sucrose as substrates, the major reaction of the FFT activity was the synthesis of tetrasaccharide. However, some synthesis of 1-kestose and re-synthesis of 6(G)-kestose also occurred. When 6, kestose was the sole substrate for the FFT activity, synthesis of tetrasaccharide was 2.7 to 3.4-fold slower than when 1-kestose was used as the sole substrate. Owing to differences in the fructan:sucrose fructosyltransferase activity of the FFT with each of the trisaccharides, net synthesis of tetrasaccharide by the FFT was altered significantly in the presence of sucrose. The magnitude of this effect depended on the concentration of the trisaccharides. In the presence of sucrose, 6(G)-kestose could be a substrate of equivalent importance to 1-kestose for synthesis of tetrasaccharide.

Beach water table fluctuations due to wave run-up: Capillarity effects

High-frequency beach water table fluctuations due to wave run-up and rundown have been observed in the field [Waddell, 1976]. Such fluctuations affect the infiltration/exfiltration process across the beach face and the interstitial oxygenation process in the beach ecosystem. Accurate representation of high-frequency water table fluctuations is of importance in the modeling of (1) the interaction between seawater and groundwater, more important, the effects on swash sediment transport and (2) the biological activities in the beach ecosystem. Capillarity effects provide a mechanism for high-frequency water table fluctuations. Previous modeling approaches adopted the assumption of saturated flow only and failed to predict the propagation of high-frequency fluctuations in the aquifer. In this paper we develop a modified kinematic boundary condition (kbc) for the water table which incorporates capillarity effects. The application of this kbc in a boundary element model enables the simulation of high-frequency water table fluctuations due to wave run-up. Numerical tests were carried out for a rectangular domain with small-amplitude oscillations; the behavior of water table responses was found to be similar to that predicted by an analytical solution based on the one-dimensional Boussinesq equation. The model was also applied to simulate the water table response to wave run-up on a doping beach. The results showed similar features of water table fluctuations observed in the field. In particular, these fluctuations are standing wave-like with the amplitude becoming increasingly damped inland. We conclude that the modified kbc presented here is a reasonable approximation of capillarity effects on beach water table fluctuations. However, further model validation is necessary before the model can confidently be used to simulate high-frequency water table fluctuations due to wave run-up.

Supersaturated solutions evaluated with an in vitro stratum corneum tape stripping technique

The flux of a compound across a membrane from any formulation, whether it contains penetration enhancers or not, is limited by its saturated solubility in the vehicle. Under such conditions the concentration of the permeant in the outer layers of the stratum corneum is also saturated. Consequently, when the permeation of a drug from a supersaturated solution leads to enhanced penetration, the concentration of the drug in the outer layers of the membrane is also supersaturated. Therefore, the stratum corneum may possess antinucleant properties which inhibit or retard the crystallisation process. In this study, the enhanced in vitro permeation of supersaturated solutions of piroxicam across human skin in diffusion cells was demonstrated. The amount of permeant in the stratum corneum was determined using a tape stripping technique. Supersaturated solutions up to four degrees of saturation were investigated which produced a linear relationship between the degree of saturation and the amount of piroxicam in the stratum corneum (R-2 = 0.970). Furthermore, the amount of piroxicam in the viable layers of the skin also increased with increasing degree of saturation. An analysis of the results suggested that enhanced penetration across human skin from supersaturated solutions of piroxicam may occur as a result of the antinucleating ability of the intercellular lipids of the stratum corneum. (C) 1997 Elsevier Science B.V.