Transport of iodide in structured clay–loam soil under maize during irrigation experiments analyzed using HYDRUS model

Transport of radioactive iodide 131I− in a structured clay loam soil under maize in a final growing phase was monitored during five consecutive irrigation experiments under ponding. Each time, 27 mm of water were applied. The water of the second experiment was spiked with 200 MBq of 131I− tracer. Its activity was monitored as functions of depth and time with Geiger-Müller (G-M) detectors in 11 vertically installed access tubes. The aim of the study was to widen our current knowledge of water and solute transport in unsaturated soil under different agriculturally cultivated settings. It was supposed that the change in 131I− activity (or counting rate) is proportional to the change in soil water content. Rapid increase followed by a gradual decrease in 131I− activity occurred at all depths and was attributed to preferential flow. The iodide transport through structured soil profile was simulated by the HYDRUS 1D model. The model predicted relatively deep percolation of iodide within a short time, in a good agreement with the observed vertical iodide distribution in soil. We found that the top 30 cm of the soil profile is the most vulnerable layer in terms of water and solute movement, which is the same depth where the root structure of maize can extend.

The charge effect on the hindrance factors for diffusion and convection of a solute in pores II

The diffusion and convection of a solute suspended in a fluid across porous membranes are known to be reduced compared to those in a bulk solution, owing to the fluid mechanical interaction between the solute and the pore wall as well as steric restriction. If the solute and the pore wall are electrically charged, the electrostatic interaction between them could affect the hindrance to diffusion and convection. In this study, the transport of charged spherical solutes through charged circular cylindrical pores filled with an electrolyte solution containing small ions was studied numerically by using a fluid mechanical and electrostatic model. Based on a mean field theory, the electrostatic interaction energy between the solute and the pore wall was estimated from the Poisson-Boltzmann equation, and the charge effect on the solute transport was examined for the solute and pore wall of like charge. The results were compared with those obtained from the linearized form of the Poisson-Boltzmann equation, i.e.the Debye-Hückel equation. © 2012 The Japan Society of Fluid Mechanics and IOP Publishing Ltd.

Conditional osmotic stress in yeast:a system to study transport through aquaglyceroporins and osmostress signaling

The accumulation and transport of solutes are hallmarks of osmoadaptation. In this study we have employed the inability of the Saccharomyces cerevisiae gpd1Δ gpd2Δ mutant both to produce glycerol and to adapt to high osmolarity to study solute transport through aquaglyceroporins and the control of osmostress-induced signaling. High levels of different polyols, including glycerol, inhibited growth of the gpd1Δ gpd2Δ mutant. This growth inhibition was suppressed by expression of the hyperactive allele Fps1-AΔ of the osmogated yeast aquaglyceroporin, Fps1. The degree of suppression correlated with the relative rate of transport of the different polyols tested. Transport studies in secretory vesicles confirmed that Fps1-Δ1 transports polyols at increased rates compared with wild type Fps1. Importantly, wild type Fps1 and Fps1-Δ1 showed similarly low permeability for water. The growth defect on polyols in the gpd1Δ gpd2Δ mutant was also suppressed by expression of a heterologous aquaglyceroporin, rat AQP9. We surmised that this suppression was due to polyol influx, causing the cells to passively adapt to the stress. Indeed, when aquaglyceroporin-expressing gpd1Δ gpd2Δ mutants were treated with glycerol, xylitol, or sorbitol, the osmosensing HOG pathway was activated, and the period of activation correlated with the apparent rate of polyol uptake. This observation supports the notion that deactivation of the HOG pathway is closely coupled to osmotic adaptation. Taken together, our "conditional" osmotic stress system facilitates studies on aquaglyceroporin function and reveals features of the osmosensing and signaling system. © 2005 by The American Society for Biochemistry and Molecular Biology, Inc.

Macrosegregation of Impurities in Directionally Solidified Silicon

Directional solidification of molten metallurgical-grade Si was carried out in a vertical Bridgman furnace. The effects of changing the mold velocity from 5 to 110 mu m seconds(-1) on the macrosegregation of impurities during solidification were investigated. The macrostructures of the cylindrical Si ingots obtained in the experiments consist mostly of columnar grains parallel to the ingot axis. Because neither cells nor dendrites can be observed on ingot samples, the absence of precipitated particles and the fulfillment of the constitutional supercooling criterion suggest a planar solid-liquid interface for mold velocities a parts per thousand currency sign10 mu m seconds(-1). Concentration profiles of several impurities were measured along the ingots, showing that their bottom and middle are purer than the metallurgical Si from which they solidified. At the ingot top, however, impurities accumulated, indicating the typical normal macrosegregation. When the mold velocity decreases, the macrosegregation and ingot purity increase, changing abruptly for a velocity variation from 20 to 10 mu m seconds(-1). A mathematical model of solute transport during solidification shows that, for mold velocities a parts per thousand yen20 mu m seconds(-1), macrosegregation is caused mainly by diffusion in a stagnant liquid layer assumed at the solid-liquid interface, whereas for lower velocities, macrosegregation increases as a result of more intense convective solute transport.

An alternative approach for the determination of soil water mobility

A new laboratory method was proposed to establish an easily performed standard for the determination of mobile soil water close to real conditions during the infiltration and redistribution of water in a soil. It consisted of applying a water volume with a tracer ion on top of an undisturbed ring sample on a pressure plate under a known suction or pressure head. Afterwards, soil water mobility was determined by analyzing the tracer-ion concentration in the soil sample. Soil water mobility showed to be a function of the applied water volume. No relation between soil water mobility and applied pressure head could be established with data from the present experiment. A simple one- or two-parameter equation can be fitted to the experimental data to parameterize soil water mobility as a function of applied solute volume. Sandy soils showed higher mobility than loamy soils at low values of applied solute volumes, and both sandy and loamy soils showed an almost complete mobility at high applied solute volumes.

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.

Commentary: Using the convection-dispersion model and transit time density functions in the analysis of organ distribution kinetics

The convection-dispersion model and its extended form have been used to describe solute disposition in organs and to predict hepatic availabilities. A range of empirical transit-time density functions has also been used for a similar purpose. The use of the dispersion model with mixed boundary conditions and transit-time density functions has been queried recently by Hisaka and Sugiyanaa in this journal. We suggest that, consistent with soil science and chemical engineering literature, the mixed boundary conditions are appropriate providing concentrations are defined in terms of flux to ensure continuity at the boundaries and mass balance. It is suggested that the use of the inverse Gaussian or other functions as empirical transit-time densities is independent of any boundary condition consideration. The mixed boundary condition solutions of the convection-dispersion model are the easiest to use when linear kinetics applies. In contrast, the closed conditions are easier to apply in a numerical analysis of nonlinear disposition of solutes in organs. We therefore argue that the use of hepatic elimination models should be based on pragmatic considerations, giving emphasis to using the simplest or easiest solution that will give a sufficiently accurate prediction of hepatic pharmacokinetics for a particular application. (C) 2000 Wiley-Liss Inc. and the American Pharmaceutical Association J Pharm Sci 89:1579-1586, 2000.

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.

Development and modulation of magnetic responsive supported ionic liquid membranes

The work presented in this thesis aims at developing a new separation process based on the application of supported magnetic ionic liquid membranes, SMILMs, using magnetic ionic liquids, MILs. MILs have attracted growing interest due to their ability to change their physicochemical characteristics when exposed to variable magnetic field conditions. The magnetic responsive behavior of MILs is thus expected to contribute for the development of more efficient separation processes, such as supported liquid membranes, where MILs may be used as a selective carrier. Driven by the MILs behavior, these membranes are expected to switch reversibly their permeability and selectivity by in situ and non-invasive adjustment of the conditions (e.g. intensity, direction vector and uniformity) of an external applied magnetic field. The development of these magnetic responsive membrane processes were anticipated by studies, performed along the first stage of this PhD work, aiming at getting a deep knowledge on the influence of magnetic field on MILs properties. The influence of the magnetic field on the molecular dynamics and structural rearrangement of MILs ionic network was assessed through a 1H-NMR technique. Through the 1H-NMR relaxometry analysis it was possible to estimate the self-diffusion profiles of two different model MILs, [Aliquat][FeCl4] and [P66614][FeCl4]. A comparative analysis was established between the behavior of magnetic and non-magnetic ionic liquids, MILs and ILs, to facilitate the perception of the magnetic field impact on MILs properties. In contrast to ILs, MILs show a specific relaxation mechanism, characterized by the magnetic dependence of their self-diffusion coefficients. MILs self-diffusion coefficients increased in the presence of magnetic field whereas ILs self-diffusion was not affected. In order to understand the reasons underlying the magnetic dependence of MILs self-diffusion, studies were performed to investigate the influence of the magnetic field on MILs’ viscosity. It was observed that the MIL´s viscosity decreases with the increase of the magnetic field, explaining the increase of MILs self-diffusion according to the modified Stokes- Einstein equation. Different gas and liquid transport studies were therefore performed aiming to determine the influence of the magnetic behavior of MILs on solute transport through SMILMs. Gas permeation studies were performed using pure CO2 andN2 gas streams and air, using a series of phosphonium cation based MILs, containing different paramagnetic anions. Transport studies were conducted in the presence and absence of magnetic field at a maximum intensity of 1.5T. The results revealed that gas permeability increased in the presence of the magnetic field, however, without affecting the membrane selectivity. The increase of gas permeability through SMILMs was related to the decrease of the MILs viscosity under magnetic field conditions.(...)

Hydrogeophysical data integration at larger scales: Application of Bayesian sequential simulation for the characterization of heterogeneous alluvial aquifers

Knowledge of the spatial distribution of hydraulic conductivity (K) within an aquifer is critical for reliable predictions of solute transport and the development of effective groundwater management and/or remediation strategies. While core analyses and hydraulic logging can provide highly detailed information, such information is inherently localized around boreholes that tend to be sparsely distributed throughout the aquifer volume. Conversely, larger-scale hydraulic experiments like pumping and tracer tests provide relatively low-resolution estimates of K in the investigated subsurface region. As a result, traditional hydrogeological measurement techniques contain a gap in terms of spatial resolution and coverage, and they are often alone inadequate for characterizing heterogeneous aquifers. Geophysical methods have the potential to bridge this gap. The recent increased interest in the application of geophysical methods to hydrogeological problems is clearly evidenced by the formation and rapid growth of the domain of hydrogeophysics over the past decade (e.g., Rubin and Hubbard, 2005).

Regional-scale integration of multi-scale hydrological and geophysical data using a two-step Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale for the purpose of improving predictions of groundwater flow and solute transport. However, extending corresponding approaches to the regional scale still represents one of the major challenges in the domain of hydrogeophysics. To address this problem, we have developed a regional-scale data integration methodology based on a two-step Bayesian sequential simulation approach. Our objective is to generate high-resolution stochastic realizations of the regional-scale hydraulic conductivity field in the common case where there exist spatially exhaustive but poorly resolved measurements of a related geophysical parameter, as well as highly resolved but spatially sparse collocated measurements of this geophysical parameter and the hydraulic conductivity. To integrate this multi-scale, multi-parameter database, we first link the low- and high-resolution geophysical data via a stochastic downscaling procedure. This is followed by relating the downscaled geophysical data to the high-resolution hydraulic conductivity distribution. After outlining the general methodology of the approach, we demonstrate its application to a realistic synthetic example where we consider as data high-resolution measurements of the hydraulic and electrical conductivities at a small number of borehole locations, as well as spatially exhaustive, low-resolution estimates of the electrical conductivity obtained from surface-based electrical resistivity tomography. The different stochastic realizations of the hydraulic conductivity field obtained using our procedure are validated by comparing their solute transport behaviour with that of the underlying ?true? hydraulic conductivity field. We find that, even in the presence of strong subsurface heterogeneity, our proposed procedure allows for the generation of faithful representations of the regional-scale hydraulic conductivity structure and reliable predictions of solute transport over long, regional-scale distances.

Use of high-resolution geophysical data to characterize heterogeneous aquifers: influence of data integration method on hydrological predictions

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. The investigation of exactly how much benefit can be brought by geophysical data in terms of its effect on hydrological predictions, however, has received considerably less attention in the literature. Here, we examine the potential hydrological benefits brought by a recently introduced simulated annealing (SA) conditional stochastic simulation method designed for the assimilation of diverse hydrogeophysical data sets. We consider the specific case of integrating crosshole ground-penetrating radar (GPR) and borehole porosity log data to characterize the porosity distribution in saturated heterogeneous aquifers. In many cases, porosity is linked to hydraulic conductivity and thus to flow and transport behavior. To perform our evaluation, we first generate a number of synthetic porosity fields exhibiting varying degrees of spatial continuity and structural complexity. Next, we simulate the collection of crosshole GPR data between several boreholes in these fields, and the collection of porosity log data at the borehole locations. The inverted GPR data, together with the porosity logs, are then used to reconstruct the porosity field using the SA-based method, along with a number of other more elementary approaches. Assuming that the grid-cell-scale relationship between porosity and hydraulic conductivity is unique and known, the porosity realizations are then used in groundwater flow and contaminant transport simulations to assess the benefits and limitations of the different approaches.

Modeling water movement in horizontal columns using fractal theory

Fractal mathematics has been used to characterize water and solute transport in porous media and also to characterize and simulate porous media properties. The objective of this study was to evaluate the correlation between the soil infiltration parameters sorptivity (S) and time exponent (n) and the parameters dimension (D) and the Hurst exponent (H). For this purpose, ten horizontal columns with pure (either clay or loam) and heterogeneous porous media (clay and loam distributed in layers in the column) were simulated following the distribution of a deterministic Cantor Bar with fractal dimension H" 0.63. Horizontal water infiltration experiments were then simulated using Hydrus 2D software. The sorptivity (S) and time exponent (n) parameters of the Philip equation were estimated for each simulation, using the nonlinear regression procedure of the statistical software package SAS®. Sorptivity increased in the columns with the loam content, which was attributed to the relation of S with the capillary radius. The time exponent estimated by nonlinear regression was found to be less than the traditional value of 0.5. The fractal dimension estimated from the Hurst exponent was 17.5 % lower than the fractal dimension of the Cantor Bar used to generate the columns.