Solute transport as affected by surface land characteristics in gravely calcareous arid soils

Gravely calcareous soils cover approximately most of arid lands (in percent); however, the solute transport behavior in these soils remains a current issue. This research aimed at estimating and correlating the solute transport parameters in gravely calcareous soils as being affected by different land uses through the knowledge of the soil morphological, physical, and chemical properties. Four different land use sites were selected: irrigated trees and bare, range, and alluvial sediment lands. Solute transport parameters of soil pore water velocity (V), dispersion coefficient (D), and retardation factor (R) were estimated using bromide breakthrough curve tests for surface soil columns. In addition, field Brilliant Blue FCF dye tracing experiment was conducted to determine the maximum dimensional movements. Soil morphological analysis was able to explain the heterogeneity in the solute transport parameters. Conductive solute transport mechanism with V of 17.99 m/day was favored in a high continuous pore system observed under tree lands. Presence of high gravel and CaCO₃ contents under range lands increased pore system tortuosity and thus increased D magnitude up to 1,339.88 cm²/day. Existence of thin surface crusts at both bare soils and alluvial sediments had considerably restricted V down to 1.46 m/day. Dye staining technique aided the explanation of the existing variations by providing visual evidence on the preferential flow paths and patterns governing the solute transport mechanism at each site.     

Sulfate transport in a Coastal Plain confining unit, New Jersey, USA

 A transient 1-D, two-pathway non-equilibrium deterministic advective dispersion model was used to examine the distribution of chloride (43–100 mg/L) and sulfate (57–894 mg/L) concentrations in the 35-m-thick section of the Lower confining unit, Atlantic Coastal Plain, New Jersey, USA. The model was used to constrain hypotheses about how pore-water chemistry changed over time. Explanations of the solute concentrations were explored by inverse and direct methods given a few known constraints, including concentrations of pore-water constituents from 12 core samples, reported simulated flow rates, and estimated hydrogeologic properties. The hypothesis that is best supported by the model results is that the distribution of chloride and sulfate concentrations in the confining unit reflect the history of the aquifer system since it was filled with seawater at the last eustatic high, about 84×10³yr BP. The model simulates fresh-water flushing of the seawater-permeated silts at a steady upward pore-water flow velocity of 8.8×10^–6 m/d, with a dispersion coefficient of 9.2×10^–7 m²/d, a dimensionless partition expression for chloride, β_Cl=0.981, and a dimensionless exchange coefficient, ω_Cl=0.31×10^–2. Sulfate concentrations were simulated over the flow path using flow and dispersion values calculated for chloride transport plus a retardation term. Parameters for sulfate transport include retardation coefficient=4.51, β_SO4=0.994, and ω_SO4=0.31×10^–2. Sensitivity analysis indicates that the model is most sensitive to flow velocity, and that fresh-water flushing of the confining unit is best simulated by having seawater concentration levels at the inflow boundary of the confining unit exponentially decrease with a concentration half-life rate of 825 yr. Received, January 1997 / Revised, April 1998, October 1998, January 1999 / Accepted, January 1999     

The estimation of dispersion behavior in discrete fractured networks of andesite in Lan-Yu Island,Taiwan

This study was based on the discrete fracture model to investigate the influence of fracture parameters on the solute transport in the fractured rocks of andesite in Lan-Yu island, Taiwan. In the simulation cases, the centers of fractures, fracture lengths and apertures were assumed to have Poisson’s distribution, negative exponential distribution and lognormal distribution, respectively. With the above assumptions, constructing the discrete fracture model became practicable. Using the mass-balance equation with specified boundary conditions, the flow field in the rock was solved. Then particles were released under the flow field. Monte Carlo method was used assuming that the amount of particles was proportional to the flow rates to get the particle accumulated percentage breakthrough curve and to estimate the dispersion coefficient. On the basis of the discrete fracture model, it was possible to evaluate the property of dispersion behavior of andesite in Lan-Yu Island with flow and transport mechanism. Properties of the dispersion behavior such as the relation between distance and traveling-time (ln〈r ²〉 and ln 〈t〉), anisotropic behavior, and the overall dispersion coefficient in a fracture network were characterized: the slope value of ln〈r ²〉 and ln〈t〉 was 1.64 an indication of non-Fickian dispersion, the particles dispersion along the flow (D11) was bigger than that perpendicular to the flow (D22), and the dispersion coefficient by this study was 0.91 m comparing the value 1 m from Sauty’s method.     

A natural-gradient field tracer test for evaluation of pollutant-transport parameters in a porous-medium aquifer

This paper describes a natural-gradient field tracer test to characterise solute-transport properties in a sand and gravel aquifer in the Hebei Province, northern China. Some laboratory-scale column tests on aquifer material and a local-scale field borehole-dilution test have been conducted previously, but the field test reported herein represents the only large-scale tracer test in the aquifer, which is the sole water supply to the city of Shi Jiazhuang and which is threatened by urban pollution. The aim of the study was to quantify the transport behaviour of nonreactive pollutants in this aquifer. Little quantitative data are available concerning its solute-transport properties; thus, the results of the tracer test are significant and critical for understanding pollutant transport and fate. The in-situ tracer test was carried out in the aquifer using a slug injection of the geochemically conservative, radioactive iodine tracer ¹³¹I. The longitudinal (α _L ) and transverse (α _T ) hydrodynamic dispersivities for solute transport in the field are 1.72 and 0.0013 m, respectively. The ratio of longitudinal dispersivity α _L and the flow length at the field scale is 1:10. The ratio between α _L and α _T from the in-situ test (~1,300:1) demonstrates a dominant longitudinal dispersion in this fluvial sand and gravel aquifer. The tracer test further indicates a relatively short transit time for the aquifer (linear velocities ~13 m/d) under natural-gradient conditions. Electronic Publication     

Estimation of contaminant transport parameters for a tropical sand in a sand tank model

This research describes the goals, design and implementation of a quasi natural gradient, laboratory scale, sand tank (aquifer) model experiment. The model was used to study the transport of an inorganic tracer (Chloride) in groundwater, within a tropical aquifer (porous medium) material. Three-dimensional sand tank (1.8 m × 0.3 m × 0.8 m) experiments were conducted to investigate contaminant transport and natural attenuation within the sand tank. In all, 360 samples were collected during 24 sampling sessions, for the three days of the tracer experiments in the Sand Tank. The Owena sand is a poorly graded sand with 88.1 % sand and 11.9 % gravel. Geotechnical properties including; coefficient of uniformity C_u = 2.53, coefficient of gradation C_z = 0.181, hydraulic conductivity K = 5.76 × 10^?4 m/s, bulk density p = 1.9 Mg/m³, effective porosity n_e = 0.215 and median grain diameter D₅₀ = 0.55 mm, were determined. Other relevant hydraulic and solute transport parameters, such as dispersion coefficients and dispersivities were also established for the tropical soil.     

On matrix diffusion: formulations, solution methods and qualitative effects

 Matrix diffusion has become widely recognized as an important transport mechanism. Unfortunately, accounting for matrix diffusion complicates solute-transport simulations. This problem has led to simplified formulations, partly motivated by the solution method. As a result, some confusion has been generated about how to properly pose the problem. One of the objectives of this work is to find some unity among existing formulations and solution methods. In doing so, some asymptotic properties of matrix diffusion are derived. Specifically, early-time behavior (short tests) depends only on φ _m ² R _m D _m  / L _m ², whereas late-time behavior (long tracer tests) depends only on φ _m R _m , and not on matrix diffusion coefficient or block size and shape. The latter is always true for mean arrival time. These properties help in: (a) analyzing the qualitative behavior of matrix diffusion; (b) explaining one paradox of solute transport through fractured rocks (the apparent dependence of porosity on travel time); (c) discriminating between matrix diffusion and other problems (such as kinetic sorption or heterogeneity); and (d) describing identifiability problems and ways to overcome them. Received, October 1997 · Revised, November 1997 · Accepted, December 1997     

Interpretation of tracer tests performed in fractured rock of the Lange Bramke basin, Germany

 Two multitracer tests performed in one of the major cross-fault zones of the Lange Bramke basin (Harz Mountains, Germany) confirm the dominant role of the fault zone in groundwater flow and solute transport. Tracers having different coefficients of molecular diffusion (deuterium, bromide, uranine, and eosine) yielded breakthrough curves that can only be explained by a model that couples the advective–dispersive transport in the fractures with the molecular diffusion exchange in the matrix. For the scale of the tests (maximum distance of 225 m), an approximation was used in which the influence of adjacent fractures is neglected. That model yielded nearly the same rock and transport parameters for each tracer, which means that the single-fracture approximation is acceptable and that matrix diffusion plays an important role. The hydraulic conductivity of the fault zone obtained from the tracer tests is about 1.5×10^–2 m/s, whereas the regional hydraulic conductivity of the fractured rock mass is about 3×10^–7 m/s, as estimated from the tritium age and the matrix porosity of about 2%. These values show that the hydraulic conductivity along the fault is several orders of magnitude larger than that of the remaining fractured part of the aquifer, which confirms the dominant role of the fault zones as collectors of water and conductors of fast flow. Received, April 1997 Revised, January 1998, August 1998 Accepted, August 1998     

Determination of groundwater solute transport parameters in finite element modelling using tracer injection and withdrawal testing data

The groundwater tracer injection and withdrawal tests are often carried out for the determination of aquifer solute transport parameters. However, the parameter analyses encounter a great difficulty due to the radial flow nature and the variability of the temporal boundary conditions. An adaptive methodology for the determination of groundwater solute transport parameters using tracer injection and withdrawal test data had been developed and illustrated through an actual case. The methodology includes the treatment of the tracer boundary condition at the tracer injection well, the normalization of tracer concentration, the groundwater solute transport finite element modelling and the method of least squares to optimize the parameters. An application of this methodology was carried out in a field test in the South of Hanoi city. The tested aquifer is Pleistocene aquifer, which is a main aquifer and has been providing domestic water supply to the city since the French time. Effective porosity of 0.31, longitudinal dispersivity of 2.2 m, and hydrodynamic dispersion coefficients from D = 220 m²/d right outside the pumping well screen to D =15.8 m²/d right outside the tracer injection well screen have been obtained for the aquifer at the test site. The minimal sum of squares of the differences between the observed and model normalized tracer concentration is 0.00119, which is corresponding to the average absolute difference between observed and model normalized concentrations of 0.035 5 (while 1 is the worst and 0 is the best fit).     

The study of hydrodynamic behaviour of a complex karst system under low-flow conditions using natural and artificial tracers (the catchment of the Unica River,SW Slovenia)

In this study a multi-tracer test with fluorescent tracers was combined with time series analyses of natural tracers to characterize the dynamics of the solute transport through different recharge pathways and to study hydraulic behaviour of a binary karst system under low-flow conditions. Fluorescent tracer testing included the introduction of uranine, amidorhodamine G, or naphthionate at three injection points. Sampling and monitoring took place at two karst springs (Malenščica, Unica) and at two underground rivers (Pivka, Rak) recharging the Unica River at the Polje of Planina, SW Slovenia. Other monitored parameters included precipitation, spring or underground river discharge, water temperature, and electrical conductivity. Water samples were collected and analyzed for total organic carbon, Mg²⁺, SO₄ ²⁻, and NO₃ ⁻ in the laboratory. In the study area, results of the tracer test suggest that contaminant transport in karst may be retarded for several weeks during low-flow conditions followed by increases in contaminant concentrations after subsequent rainfall events. Based on interpretation of tracer concentration breakthrough curves, low apparent dominant flow velocities (i.e., between 5.8 and 22.8 m/h through the well developed karst conduits, and 3.6 m/h through the prevailing vadose zone with a dominant influence of a diffuse recharge) were obtained. Together with analyses of hydro-chemographs the artificial tracing identified different origins of water recharging the studied aquifer. During prolonged low-water conditions the Malenščica spring is mainly recharged from the karst aquifer and the Unica spring by the sinking Pivka River. After more intensive rainfall events allogenic recharge from Cerknica prevails in the Malenščica spring, while the Unica spring drains mainly the allogenic water from the Pivka Valley. These findings of alternating hydraulic connections and drainage areas due to respective hydrological conditions are important and should be considered when monitoring water quality, implementing groundwater protection measures, and optimizing future water exploitation.     

Experimental Study on Surface Reactions of Heavy Metal Ions With Quartz—Aqueous Ion Concentration Dependence

Adsorption of divalent metal ions, including Cu²⁺, Pb²⁺, Zn²⁺, Cd²⁺ and Ni²⁺, on quartz surface was measured as a function of metal ion concentration at 30°C under conditions of solution pH= 6. 5 and ion strength I = 0. 1mol/L. Results of the experimental measurements can be described very well by adsorption isotherm equations of Freudlich. The correlation coefficients (r) of adsorption isotherm lines are > 0. 96. Moreover, the experimental data were interpreted on the basis of surface complexation model. The experimental results showed that the monodentate-coordinated metal ion surface complex species (SOM⁺) are predominant over the bidentate-coordinated metal ion surface complex species [(SO)₂M] formed only by the ions Cu²⁺, Zn²⁺ and Ni²⁺. And the relevant apparent surface complexation constants are lgK_M = 2.2–3.3 in order of K_Cd≥K_Pb > K_Zn > K_Ni≥K_Cu, and lgβ_M = 5.9-6.8 in order of β_Ni > β_Zn > β_Cu. Therefore, the reactive ability of the ions onto mineral surface of quartz follows the order of Cd > Pb > Zn > Ni> Cu under the above-mentioned solution conditions. The apparent surface complexation constants, influenced by the surface potential, surface species and hydrolysis of metal ions, depend mainly on the Born solvation coefficient of the metal ions. This project was financially supported by the National Natural Science Foundation of China (No. 49572091).     

In situ push-pull tests for the determination of TCE degradation and permanganate consumption rates

A simple, single-well push-pull test was conducted at a TCE-contaminated site to estimate the site-specific TCE degradation and permanganate (MnO₄⁻) consumption rate. Known quantities of a conservative tracer (Br⁻) and permanganate were rapidly injected into a saturated aquifer then periodically sampled during extraction from the same well. Concentrations of Br⁻, TCE, and MnO₄⁻ were measured; breakthrough curves (BTCs) for all species of solute were determined. Data analysis of BTCs for bromide and TCE showed that the first-order rate constant of TCE degradation by MnO₄⁻ is 1.67 ± 0.152 h⁻¹. Further, the in situ MnO₄⁻ demand rate by TCE and aquifer materials is estimated to be 0.54 ± 0.371 h⁻¹. This study demonstrates that in situ push-pull tests are useful and economical tools for field investigations to determine contaminant reaction and oxidant consumption rates, which may then be used to optimize groundwater remediation designs.     

Ca–Mg diffusion in diopside: tracer and chemical inter-diffusion coefficients

We have experimentally determined the tracer diffusion coefficients (D*) of ⁴⁴Ca and ²⁶Mg in a natural diopside (~Di₉₆) as function of crystallographic direction and temperature in the range of 950–1,150 °C at 1 bar and f(O₂) corresponding to those of the WI buffer. The experimental data parallel to the a*, b, and c crystallographic directions show significant diffusion anisotropy in the a–c and b–c planes, with the fastest diffusion being parallel to the c axis. With the exception of logD*(²⁶Mg) parallel to the a* axis, the experimental data conform to the empirical diffusion “compensation relation”, converging to logD ~ −19.3 m²/s and T ~ 1,155 °C. Our data do not show any change of diffusion mechanism within the temperature range of the experiments. Assuming that D* varies roughly linearly as a function of angle with respect to the c axis in the a–c plane, at least within a limited domain of ~20° from the c-axis, our data do not suggest any significant difference between D*(//c) and D*(⊥(001)), the latter being the diffusion data required to model compositional zoning in the (001) augite exsolution lamellae in natural clinopyroxenes. Since the thermodynamic mixing property of Ca and Mg is highly nonideal, calculation of chemical diffusion coefficient of Ca and Mg must take into account the effect of thermodynamic factor (TF) on diffusion coefficient. We calculate the dependence of the TF and the chemical interdiffusion coefficient, D(Ca–Mg), on composition in the diopside–clinoenstatite mixture, using the available data on mixing property in this binary system. Our D*(Ca) values parallel to the c axis are about 1–1.5 log units larger than those Dimanov et al. (1996). Incorporating the effect of TF, the D(Ca–Mg) values calculated from our data at 1,100–1,200 °C is ~0.6–0.7 log unit greater than the experimental quasibinary D((Ca–Mg + Fe)) data of Fujino et al. (1990) at 1 bar, and ~0.6 log unit smaller than that of Brady and McCallister (1983) at 25 kb, 1,150 °C, if our data are normalized to 25 kb using activation volume (~4 and ~6 cm³/mol for Mg and Ca diffusion, respectively) calculated from theoretical considerations.     

Dualistic distribution coefficients of elements in the system mineral-hydrothermal solution. I. Gold accumulation in pyrite

The use of trace elements (TE) as geochemical indicators is complicated by the dualism of their distribution coefficients D due to the additional (i.e., above the concentrations of an isomorphic component) incorporation of elements at structural defects of various nature (including the surface of the crystal). A pressing problem in this situation is to determine the true D values that pertain to the structural component of an admixture D ^str and evaluate effects of other modes of TE occurrence. Only upon distinguishing D ^str in the bulk coefficient D ^bulk it is possible to evaluate the ore potential of fluid in terms of certain TE from the composition of a mineral containing the TE. Pyrite synthesized in solutions of variable pH at 450°C and 1 kbar (100 MPa) at fluid portions sampled in a trap is utilized to demonstrate the role of a surface nonautonomous phase (NP) in the incorporation of gold in this mineral. The distribution coefficient of gold between pyrite and hydrothermal solution is 0.14 for “pure” pyrite and 0.05 for As-bearing pyrite (containing 0.02–0.05 wt % As), and these coefficients for NP are 310 and 170, respectively. This increases the D ^bulk for evenly distributed (“invisible”) gold by factors of four and nine. In contrast to the results of earlier studies conducted at room temperature and pressure or parameters close to them, our data demonstrate that the accumulation of “invisible” Au in pyrite is controlled not only by reducing adsorption with the development of Au(0) particles and films but also by Au incorporation in NP developing in the surface layer of the crystal approximately 500 nm thick as chemically bound Au [most likely as Au(I)]. The possible reason for the high absorption capacity of NP is the defect (pyrrhotite-like) structure, which is not saturated with bonds of excess S and sulfoxi onions.     

An SF<Subscript>6</Subscript> Tracer Study of the Flow Dynamics in the Stockton Deep Water Ship Channel: Implications for Dissolved Oxygen Dynamics

A sulfur hexafluoride (SF₆) tracer release experiment was conducted in the Stockton Deep Water Ship Channel (DWSC) to quantify mixing and transport rates. SF₆ was injected in the San Joaquin River upstream of the DWSC and mapped for 8 days. From the temporal change in SF₆ distributions, the longitudinal dispersion coefficient (K _x ) was determined to be 32.7 ± 3.6 m² s⁻¹ and the net velocity was 1.75 ± 0.03 km day⁻¹. Based on the decrease in SF₆ inventory during the experiment, the pulsed residence time for waters in the DWSC was estimated at ∼17 days. Within the DWSC from Stockton downstream to Turner Cut, dissolved oxygen concentrations maintained a steady state value of 4 mg l⁻¹. These values are below water quality objectives for the time of year. The low flow rates observed in the DWSC and the inability of oxygen-rich waters from downstream to mix into the DWSC upstream of Turner Cut contribute to the low dissolved oxygen concentration.     

Experimental study of migration of potassium ion through a two-layer soil system

A barrier system based on the hydraulic trap design concept for a landfill was proposed. To study the field scenario in which a clay liner is underlain by a granular layer functioning as a secondary leachate drain layer, a laboratory advection–diffusion test was performed to investigate factors controlling the transport of contaminants in a two-layer soil system. The soils used for this study were Ariake clay and, the underlying layer, Shirasu soil from the Kyushu region of Japan. Potassium (K⁺) was selected as the target chemical species with an initial concentration of 905 mg L⁻¹. The effective diffusion coefficients (D _e) of K⁺ for Ariake clay and Shirasu soil were back-calculated using an available computer program, Pollute V 6.3. Values of D _e derived from this experiment are consistent with previously published ones. The Ariake clay has lower D _e than the Shirasu soil. The hypothesis that mechanical dispersion can be considered negligible is reasonable based on both the observation that the predicted values well fit the experimental data and the analyses of two dimensionless parameters. Parametric analyses show that transport of K⁺ through soils is controlled by advection–diffusion rather than diffusion only, whereas at low Darcy velocity (i.e., ≤10⁻⁹ m s⁻¹), transport of K⁺ will be controlled by diffusion. Applications of the test results and parametric analysis results in practical situations were reviewed.     

Point defects and cation tracer diffusion in (Ti<Subscript><Emphasis Type="Italic">x</Emphasis></Subscript> Fe<Subscript><Emphasis Type="Italic">1−x</Emphasis></Subscript>)<Subscript>3−δ</Subscript>O<Subscript>4</Subscript>. II. Cation tracer diffusion

 Cation tracer diffusion coefficients, D_Me ^*, for Me=Fe, Mn, Co and Ti, were measured using radioactive isotopes in the spinel solid solution (Ti _x Fe _1−x )_3−δO₄ as a function of the oxygen activity. Experiments were performed at different cationic compositions (x=0, 0.1, 0.2 and 0.3) at 1100, 1200, 1300 and 1400 °C. The oxygen activity dependence of all data for D_Me ^* at constant temperature and cationic composition can be described by equations of the type D_Me ^*=D^∘ _Me[V]. C_V·a _O2 ^2/3+D_Me[I] ^∘·a _O2 ^−2/3·D_Me[V] ^∘ and D_Me[I] ^∘ are constants and C_V is a factor of the order of unity which decreases with increasing δ. All log D_Me ^* vs. loga _O2 curves obtained for different values of x and for different temperatures go through a minimum due to a change in the type of point defects dominating the cation diffusion with oxygen activity. Cation vacancies prevail for the cation diffusion at high oxygen activities while cation interstitials become dominant at low oxygen activities. At constant values of x, D_Me[V] ^∘ decreases with increasing temperature while D_Me[I] ^∘ increases.     

Roughness length and drag coefficient at two MONTBLEX-90 tower stations

Using MONTBLEX-90 mean velocity data, roughness lengths and drag coefficients are estimated at Jodhpur and Kharagpur. At Jodhpur, since the surface is not uniform the roughness length is estimated separately in three different subsectors within the range of prevailing wind directions and averages to 1.23 cm in the sector between 200° and 230° which is relatively flat with no obstacles on the ground. At Kharagpur, where the terrain is more nearly homogeneous, the average value (for all prevailing wind directions) is 1.94 cm. The drag coefficient C_D at Jodhpur shows variation both with the roughness subsector and with wind speed, the average over all directions increasing rapidly as themean wind speed Ū₁₀ at 10m height drops according to the power lawC _D = 0.05 Ū ₁₀ ^t-1.09 in trie range 0.5 < Ū₁₀ < 7 m s⁻¹. At Kharagpur, the drag coefficient is smaller than at Jodhpur by nearly 50% for the same range of wind speeds (> 3 ms⁻¹).     

Experimental determination of three-dimensional dispersivities in homogeneous porous medium

Laboratory tracer experiments were conducted to investigate longitudinal dispersivity (α _x ) as well as the transversal (α _y ) and vertical (α _z ) dispersivities in homogeneous 3–5 mm sandy aquifer. The experiments were carried out in a channel 12-m long, 1.35-m wide and 0.60-m high which was built in the Hydraulics Laboratory of Civil Engineering Department in Dokuz Eylul University. NaCl was used as a tracer and conductivity values were measured at 220 measurement points. Mass Transport 3 Dimensional (Zheng and Wang in SERDP-99-1, US Army Engineer Research and Development Center, Vicksburg, MS, 1999; MT3DMS code) which is a three-dimensional solute transport simulation model incorporating finite differences solution option was used to solve the three-dimensional advective–dispersive transport equation. The estimated dispersivity values were modified until an acceptable compatibility between the observed and calculated concentrations at measurement points was reached. The best match was obtained for α _x  = 12 cm, α _y /α _x  = 0.2 and α _z /α _x  = 0.05. These values are compatible with those encountered in the literature.     

Nitrate vertical transport and simulation in soils of rocky desertification in Karst regions,Southwest China

The characteristics of nitrate vertical transport in soils collected from Libo and Puding in Guizhou Province were studied by simulating soil column in laboratory. The results were as follows: (1) Vertical transport velocity of nitrate decreased, and the breakthrough curves (BTCs) of nitrate were more dispersed, in each horizon from surface layer to bottom layer in every soil profile. As rocky desertification progressed, the BTCs experienced a gentle up and down trend, and tailing was more obvious. (2) An analytical solute transport model (CXTFIT 2.0) was used to estimate nitrate dispersion coefficient (D) and average pore water velocity (V) from the observed BTCs. The results showed that CXTFIT 2.0 model was suitable in fitting the nitrate transport in these soils. The dispersion coefficient was found to be a function of average pore water velocity. (3) The transport of nitrate was mainly affected by the soil structural coefficient. As soil structural coefficient decreased, nitrate outflow was retarded, and the peak concentration was reduced. Soil bulk density, organic matter, and clay also affected the vertical transport of nitrate. Low bulk density, clay content, and high organic matter content were each associated with faster nitrate transport.     

Oxygen diffusion in basalt and andesite melts: experimental results and discussion of chemical versus tracer diffusion

Chemical diffusion coefficients for oxygen in melts of Columbia River basalt (Ice Harbor Dam flow) and Mt. Hood andesite have been determined at 1 atm. The diffusion model is that of sorption or desorption of oxygen into a sphere of uniform initial concentration from a constant and semi-infinite atmosphere. The experimental design utilizes a thermogravimetric balance to monitor the rate of weight change arising from the response of the sample redox state to an imposed fO₂. Oxygen diffusion coefficients are approximately an order-ofmagnitude greater for basaltic melt than for andesitic melt. At 1260° C, the oxygen diffusion coefficients are: D=1.65×10^–6cm²/s and D=1.43×10^–7cm²/s for the basalt and andesite melts, respectively. The high oxygen diffusivity in basaltic melt correlates with a high ratio of nonbridging oxygen/tetrahedrally coordinated cations, low melt viscosity, and high contents of network-modifying cations. The dependence of the oxygen diffusion coefficient on temperature is: D=36.4exp(–51,600±3200/RT)cm²/s for the basalt and D=52.5exp(–60,060±4900/RT)cm²/s for the andesite (R in cal/deg-mol; T in Kelvin). Diffusion coefficients are independent of the direction of oxygen diffusion (equilibrium can be approached from extremely oxidizing or reducing conditions) and thus, melt redox state. Characteristic diffusion distances for oxygen at 1260° C vary from 10^-2 to 10² m over the time interval of 1 to 10⁶ years. A compensation diagram shows two distinct trends for oxygen chemical diffusion and oxygen tracer diffusion. These different linear relationships are interpreted as supporting distinct oxygen transport mechanisms. Because oxygen chemical diffusivities are generally greater than tracer diffusivities and their Arrhenius activation energies are less, transport mechanisms involving either molecular oxygen or vacancy diffusion are favored.