CXTANNEAL: an improved program for estimating solute transport parameters

CXTANNEAL is a program for analysing contaminant transport in soils. The code, written in Fortran 77, is a modified version of CXTFIT, a commonly used package for estimating solute transport parameters in soils. The improvement of the present code is that it includes simulated annealing as the optimization technique for curve fitting. Tests with hypothetical data show that CXTANNEAL performs better than the original code in searching for optimal parameter estimates. To reduce the computational time, a parallel version of CXTANNEAL (CXTANNEAL_P) was also developed. (C) 1999 Elsevier Science Ltd. All rights reserved.

Effect of companion ion on nitrate displacement through transport parameters analysis

Nitrate is the main form of nitrogen associated with water contamination; the high mobility of this species in soil justifies the concern regarding nitrogen management in agricultural soils. Therefore, the objective of this research was to assess the effect of companion cation on nitrate displacement, by analyzing nitrate transport parameters through Breakthrough Curves (BTCs) and their settings made by numerical model (STANMOD). The experiment was carried out in the Soil and Water Quality Laboratory of the Department of Biosystems Engineering, "Luiz de Queiroz" College of Agriculture in Piracicaba (SP), Brazil. It was performed using saturated soil columns in steady-state flow condition, in which two different sources of inorganic nitrate Ca(NO3)2 and NH4NO3 were applied at a concentration of 50 mg L-1 NO3-. Each column was filled with either a Red-Yellow Oxisol (S1) or an Alfisol (S2). Results are indicative that the companion ion had no effect on nitrate displacement. However, nitrate transport was influenced by soil texture, particle aggregation, solution speed in soil and organic matter presence. Nitrate mobility was higher in the Alfisol (S2).

A dedicated setup for the measurement of the electron transport parameters in gases at large electric fields

Electron transport parameters are important in several areas ranging from particle detectors to plasma-assisted processing reactors. Nevertheless, especially at high fields strengths and for complex gases, relatively few data are published. A dedicated setup has been developed to measure the electron drift velocity and the first Townsend coefficient in parallel plate geometry. An RPC-like cell has been adopted to reach high field strengths without the risk of destructive sparks. The validation data obtained with pure Nitrogen will be presented and compared to a selection of the available literature and to calculations performed with Magboltz 2 version 8.6. The new data collected in pure Isobutane will then be discussed. This is the first time the electron drift velocity in pure Isobutane is measured well into the saturation region. Good agreement is found with expectations from Magboltz. (C) 2009 Elsevier B.V. All rights reserved.

Metal transport parameters of a gneiss saprolitic silty soil for liner design

Diffusion coefficients and retardation factors of two metal cations (Cd2+ and Pb2+) were measured for a compacted Brazilian saprolitic soil derived from gneiss, aiming to assess its geoenvironmental performance as a liner for waste disposal sites. This soil occurs extensively all over the country in very thick layers, but has not been used in liners because of its hydraulic conductivity, higher than 10(-9) m/s when compacted at optimum water content of standard Proctor energy, but which can be reduced by means of appropriate compaction techniques or additives. Batch, column, and diffusion tests were carried out with monospecies synthetic solutions at pH 1, 3, and 5.5. Measured diffusion coefficients varied between 0.5 and 4 X 10(-10) m(2)/s. Retardation factors show that cadmium, a very mobile cation, is not adsorbed at pH I but is significantly retained at pH 3 and pH 5.5, whereas lead is retained at all tested pH values though slightly at pH 1. Estimated retardation factors from batch tests were 1.3-2.3 times those resulting from column tests and at its highest when obtained by diffusion tests; whereas batch tests allow a more complete exposure of the soil grains to the solution, time-dependent nonspecific adsorption may take longer to occur. The importance of contact time was observed and should be considered in further investigations. Its significant retention of metals suggests a promising utilization of this soil as a bottom liner for wastes landfills.

Spectroelectrochemistry and investigation of charge transport mechanisms of iron poly(pyridyl) redox polymers

In this work, we describe the characterization of the complex [Fe(tpy-NH2)(2)](PF6)(2) (tpy-NH2 = bis[4`-(3-aminophenyl)-2, 2`:6`,2 ``-terpyridine]. The complex was oxidatively electropolymerized on glassy.-carbon electrodes in CH3CN/0.1 M tetraethylammonium perchlorate (TEAP) to generate polymer films that exhibit reversible oxidative electrochemical behavior in a wide potential range (0.0-1.6 V), as well as high conductivity and stability/durability. In situ spectrocyclic voltammetry of this modified electrode was carried out on a photodiode array spectrophotometer attached to a potentiostat, which provided UV-Vis absorption spectra of the redox species during the potential sweep. We determined charge transport parameters as a function of time and thickness of the modified electrode, and the results showed that poly-[[Fe(tpy-NH2)(2)](2+)](n) can be made to exhibit three regimes of charge transport behavior by manipulation of the film thickness and the experimental time-scale. Morphological characterization of the film was provided by atomic force microscopy. (C) 2008 Elsevier Ltd. All rights reserved.

Inferring transport characteristics in a fractured rock aquifer by combining single-hole ground-penetrating radar reflection monitoring and tracer test data

Investigations of solute transport in fractured rock aquifers often rely on tracer test data acquired at a limited number of observation points. Such data do not, by themselves, allow detailed assessments of the spreading of the injected tracer plume. To better understand the transport behavior in a granitic aquifer, we combine tracer test data with single-hole ground-penetrating radar (GPR) reflection monitoring data. Five successful tracer tests were performed under various experimental conditions between two boreholes 6 m apart. For each experiment, saline tracer was injected into a previously identified packed-off transmissive fracture while repeatedly acquiring single-hole GPR reflection profiles together with electrical conductivity logs in the pumping borehole. By analyzing depth-migrated GPR difference images together with tracer breakthrough curves and associated simplified flow and transport modeling, we estimate (1) the number, the connectivity, and the geometry of fractures that contribute to tracer transport, (2) the velocity and the mass of tracer that was carried along each flow path, and (3) the effective transport parameters of the identified flow paths. We find a qualitative agreement when comparing the time evolution of GPR reflectivity strengths at strategic locations in the formation with those arising from simulated transport. The discrepancies are on the same order as those between observed and simulated breakthrough curves at the outflow locations. The rather subtle and repeatable GPR signals provide useful and complementary information to tracer test data acquired at the outflow locations and may help us to characterize transport phenomena in fractured rock aquifers.

A Database of >20 keV Electron Green's Functions of Interplanetary Transport at 1 AU

We use interplanetary transport simulations to compute a database of electron Green's functions, i.e., differential intensities resulting at the spacecraft position from an impulsive injection of energetic (>20 keV) electrons close to the Sun, for a large number of values of two standard interplanetary transport parameters: the scattering mean free path and the solar wind speed. The nominal energy channels of the ACE, STEREO, and Wind spacecraft have been used in the interplanetary transport simulations to conceive a unique tool for the study of near-relativistic electron events observed at 1 AU. In this paper, we quantify the characteristic times of the Green's functions (onset and peak time, rise and decay phase duration) as a function of the interplanetary transport conditions. We use the database to calculate the FWHM of the pitch-angle distributions at different times of the event and under different scattering conditions. This allows us to provide a first quantitative result that can be compared with observations, and to assess the validity of the frequently used term beam-like pitch-angle distribution.

Temperature effect in potassium and nitrate ions in soil transport

When doing researches on solute dynamics in porous medium, the knowledge of medium characteristics and percolating liquids, as well as of external factors is very important. An important external factor is temperature and, in this sense, our purpose was determining potassium and nitrate transport parameters for different values of temperature, in miscible displacement experiments. Evaluated parameters were retardation factor (R), diffusion/dispersion coefficient (D) and dispersivity, at ambient temperature (25 up to 28 ºC), 40 ºC and 50 ºC. Salts used were potassium nitrate and potassium chlorate, prepared in a solution made up of 5 ppm nitrate and 2.000 ppm potassium, with Red-Yellow Latosol porous medium. Temperature exhibited a positive influence upon porous medium solution and upon dispersion coefficient.

On the observability of turbulent transport rates by Argo: supporting evidence from an inversion experiment

Although estimation of turbulent transport parameters using inverse methods is not new, there is little evaluation of the method in the literature. Here, it is shown that extended observation of the broad scale hydrography by Argo provides a path to improved estimates of regional turbulent transport rates. Results from a 20 year ocean state estimate produced with the ECCO v4 non-linear inverse modeling framework provide supporting evidence. Turbulent transport parameter maps are estimated under the constraints of fitting the extensive collection of Argo profiles collected through 2011. The adjusted parameters dramatically reduce misfits to in situ profiles as compared with earlier ECCO solutions. They also yield a clear reduction in the model drift away from observations over multi-century long simulations, both for assimilated variables (temperature and salinity) and independent variables (bio-geochemical tracers). Despite the minimal constraints imposed specifically on the estimated parameters, their geography is physically plausible and exhibits close connections with the upper ocean ocean stratification as observed by Argo. The estimated parameter adjustments furthermore have first order impacts on upper-ocean stratification and mixed layer depths over 20 years. These results identify the constraint of fitting Argo profiles as an effective observational basis for regional turbulent transport rates. Uncertainties and further improvements of the method are discussed.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

Mixing-cell boundary conditions and apparent mass balance errors for advective-dispersive solute transport

In many field or laboratory situations, well-mixed reservoirs like, for instance, injection or detection wells and gas distribution or sampling chambers define boundaries of transport domains. Exchange of solutes or gases across such boundaries can occur through advective or diffusive processes. First we analyzed situations, where the inlet region consists of a well-mixed reservoir, in a systematic way by interpreting them in terms of injection type. Second, we discussed the mass balance errors that seem to appear in case of resident injections. Mixing cells (MC) can be coupled mathematically in different ways to a domain where advective-dispersive transport occurs: by assuming a continuous solute flux at the interface (flux injection, MC-FI), or by assuming a continuous resident concentration (resident injection). In the latter case, the flux leaving the mixing cell can be defined in two ways: either as the value when the interface is approached from the mixing-cell side (MC-RT -), or as the value when it is approached from the column side (MC-RT +). Solutions of these injection types with constant or-in one case-distance-dependent transport parameters were compared to each other as well as to a solution of a two-layer system, where the first layer was characterized by a large dispersion coefficient. These solutions differ mainly at small Peclet numbers. For most real situations, the model for resident injection MC-RI + is considered to be relevant. This type of injection was modeled with a constant or with an exponentially varying dispersion coefficient within the porous medium. A constant dispersion coefficient will be appropriate for gases because of the Eulerian nature of the usually dominating gaseous diffusion coefficient, whereas the asymptotically growing dispersion coefficient will be more appropriate for solutes due to the Lagrangian nature of mechanical dispersion, which evolves only with the fluid flow. Assuming a continuous resident concentration at the interface between a mixing cell and a column, as in case of the MC-RI + model, entails a flux discontinuity. This flux discontinuity arises inherently from the definition of a mixing cell: the mixing process is included in the balance equation, but does not appear in the description of the flux through the mixing cell. There, only convection appears because of the homogeneous concentration within the mixing cell. Thus, the solute flux through a mixing cell in close contact with a transport domain is generally underestimated. This leads to (apparent) mass balance errors, which are often reported for similar situations and erroneously used to judge the validity of such models. Finally, the mixing cell model MC-RI + defines a universal basis regarding the type of solute injection at a boundary. Depending on the mixing cell parameters, it represents, in its limits, flux as well as resident injections. (C) 1998 Elsevier Science B.V. All rights reserved.

Multi-scale micro-structure generation strategy for up-scaling transport in clays

Clays and claystones are used as backfill and barrier materials in the design of waste repositories, because they act as hydraulic barriers and retain contaminants. Transport through such barriers occurs mainly by molecular diffusion. There is thus an interest to relate the diffusion properties of clays to their structural properties. In previous work, we have developed a concept for up-scaling pore-scale molecular diffusion coefficients using a grid-based model for the sample pore structure. Here we present an operational algorithm which can generate such model pore structures of polymineral materials. The obtained pore maps match the rock’s mineralogical components and its macroscopic properties such as porosity, grain and pore size distributions. Representative ensembles of grains in 2D or 3D are created by a lattice Monte Carlo (MC) method, which minimizes the interfacial energy of grains starting from an initial grain distribution. Pores are generated at grain boundaries and/or within grains. The method is general and allows to generate anisotropic structures with grains of approximately predetermined shapes, or with mixtures of different grain types. A specific focus of this study was on the simulation of clay-like materials. The generated clay pore maps were then used to derive upscaled effective diffusion coefficients for non-sorbing tracers using a homogenization technique. The large number of generated maps allowed to check the relations between micro-structural features of clays and their effective transport parameters, as is required to explain and extrapolate experimental diffusion results. As examples, we present a set of 2D and 3D simulations and investigated the effects of nanopores within particles (interlayer pores) and micropores between particles. Archie’s simple power law is followed in systems with only micropores. When nanopores are present, additional parameters are required; the data reveal that effective diffusion coefficients could be described by a sum of two power functions, related to the micro- and nanoporosity. We further used the model to investigate the relationships between particle orientation and effective transport properties of the sample.

An improved procedure for obtaining and maintaining well characterized partial water saturation states on concrete samples to be used for mass transport tests

A conditioning procedure is proposed allowing to install into the concrete specimens any selected value of water saturation degree with homogeneous moisture distribution. This is achieved within the least time and the minimum alteration of the concrete specimens. The protocol has the following steps: obtaining basic drying data at 50 °C (water absorption capacity and drying curves); unidirectional drying of the specimens at 50 °C until reaching the target saturation degree values; redistribution phase in closed containers at 50 °C (with measurement of the quasi-equilibrium relative humidities); storage into controlled environment chambers until and during mass transport tests, if necessary. A water transport model is used to derive transport parameters of the tested materials from the drying data, i.e., relative permeabilities and apparent water diffusion coefficients. The model also allows calculating moisture profiles during isothermal drying and redistribution phases, thus allowing optimization of the redistribution times for obtaining homogeneous moisture distributions.

Fate and Transport of Naproxen in a Sandy Aquifer Material: Saturated Column Studies and Model Evaluation

Naproxen-C14H14O3 is a nonsteroidal anti-inflammatory drug which has been found at detectable concentrations in wastewater, surface water, and groundwater. Naproxen is relatively hydrophilic and is in anionic form at pH between 6 and 8. In this study, column experiments were performed using an unconsolidated aquifer material from an area near Barcelona (Spain) to assess transport and reaction mechanisms of Naproxen in the aquifer matrix under different pore water fluxes. Results were evaluated using HYDRUS-1D, which was used to estimate transport parameters. Batch sorption isotherms for Naproxen conformed with the linear model with a sorption coefficient of 0.42 (cm3 g−1), suggesting a low sorption affinity. Naproxen breakthrough curves (BTCs) measured in soil columns under steady-state, saturated water flow conditions displayed similar behavior, with no apparent hysteresis in sorption or dependence of retardation (R, 3.85-4.24) on pore water velocities. Soil sorption did not show any significant decrease for increasing flow rates, as observed from Naproxen recovery in the effluent. Sorption parameters estimated by the model suggest that Naproxen has a low sorption affinity to aquifer matrix. Most sorption of Naproxen occurred on the instantaneous sorption sites, with the kinetic sorption sites representing only about 10 to 40% of total sorption.