An investigation of localised soil heterogeneities on solute transport using a multisegment percolation system

A multisegment percolation system (MSPS) consisting of 25 individual collection wells was constructed to study the effects of localised soil heterogeneities on the transport of solutes in the vadose zone. In particular, this paper discusses the transport of water and nutrients (NO3-, Cl-, PO43-) through structurally stable, free-draining agricultural soil from Victoria, Australia. A solution of nutrients was irrigated onto the surface of a large undisturbed soil core over a 12-h period. This was followed by a continuous irrigation of distilled water at a fate which did not cause pending for a further 18 days. During this time, the volume of leachate and the concentration of nutrients in the leachate of each well were measured. Very significant variation in drainage patterns across a small spatial scale was observed. Leaching of nitrate-nitrogen and chloride from the core occurred two days after initial application. However, less than 1% of the total applied phosphate-phosphorus leached from the soil during the 18-day experiment, indicating strong adsorption. Our experiments indicate considerable heterogeneity in water flow patterns and solute leaching on a small spatial scale. These results have significant ramifications for modelling solute transport and predicting nutrient loadings on a larger scale.

Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

Quantifying the influence of intra-aggregate concentration gradients on solute transport

The physical nonequilibrium of solute concentration resulting from preferential now of soil water has often led to models where the soil is partitioned into two regions: preferential flow paths, where solute transport occurs mainly by advection, and the remaining region, where significant solute transport occurs through diffusive exchange with the flow paths. These two-region models commonly ignore concentration gradients within the regions. Our objective was to develop a simple model to assess the influence of concentration gradients on solute transport and to compare model results with experiments conducted on structured materials. The model calculates the distribution of solutes in a single spherical aggregate surrounded by preferential now paths and subjected to alternating boundary conditions representing either an exchange of solutes between the two regions (a wet period) or no exchange but redistribution of solutes within the aggregate (a dry period). The key parameter in the model is the aggregate radius, which defines the diffusive time scales. We conducted intermittent leaching experiments on a column of packed porous spheres and on a large (300 mm long by 216 mm diameter) undisturbed field soil core to test the validity of the model and its application to field soils. Alternating wet and dry periods enhanced leaching by up to 20% for this soil, which was consistent with the model's prediction, given a fitted equivalent aggregate radius of 1.8 cm, If similar results are obtained for other soils, use of alternating wet and dry periods could improve management of solutes, for example in salinity control and in soil remediation.

Comparison of coherent and weakly incoherent transport models for the interlayer magnetoresistance of layered Fermi liquids

The interlayer magnetoresistance of layered metals in a tilted magnetic field is calculated for two distinct models for the interlayer transport. The first model involves coherent interlayer transport, and makes use of results of semiclassical or Bloch-Boltzmann transport theory. The second model involves weakly incoherent interlayer transport where the electron is scattered many times within a layer before tunneling into the next layer. The results are relevant to the interpretation of experiments on angular-dependent magnetoresistance oscillations (AMRO) in quasi-one- and quasi-two-dimensional organic metals. We find that the dependence of the magnetoresistance on the direction of the magnetic field is identical for both models except when the field is almost parallel to the layers. An important implication of this result is that a three-dimensional Fermi surface is not necessary for the observation of the Yamaji and Danner oscillations seen in quasi-two- and quasi-one-dimensional metals, respectively. A universal expression is given for the dependence of the resistance at AMRO maxima and minima on the magnetic field and scattering time (and thus the temperature). We point out three distinctive features of coherent interlayer transport: (i) a beat frequency in the magnetic oscillations of quasi-two-dimensional systems, (ii) a peak in the angular-dependent magnetoresistance when the field is sufficiently large and parallel to the layers, and (iii) a crossover from a linear to a quadratic field dependence for the magnetoresistance when the field is parallel to the layers. Properties (i) and (ii) are compared with published experimental data for a range of quasi-two-dimensional organic metals. [S0163-1829(99)02236-5].

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.

Transport properties of strongly correlated metals: A dynamical mean-field approach

The temperature dependence of the transport properties of the metallic phase of a frustrated Hubbard model on the hypercubic lattice at half-filling is calculated. Dynamical mean-held theory, which maps the Hubbard model onto a single impurity,Anderson model that is solved self-consistently, and becomes exact in the limit of large dimensionality, is used. As the temperature increases there is a smooth crossover from coherent Fermi liquid excitations at low temperatures to incoherent excitations at high temperatures. This crossover leads to a nonmonotonic temperature dependence for the resistance, thermopower, and Hall coefficient, unlike in conventional metals. The resistance smoothly increases from a quadratic temperature dependence at low temperatures to large values which can exceed the Mott-Ioffe-Regel value ha/e(2) (where a is a lattice constant) associated with mean free paths less than a lattice constant. Further signatures of the thermal destruction of quasiparticle excitations are a peak in the thermopower and the absence of a Drude peak in the optical conductivity. The results presented here are relevant to a wide range of strongly correlated metals, including transition metal oxides, strontium ruthenates, and organic metals.

Solute transport by surface runoff from low-angle slopes: theory and application

The removal of chemicals in solution by overland how from agricultural land has the potential to be a significant source of chemical loss where chemicals are applied to the soil surface, as in zero tillage and surface-mulched farming systems. Currently, we lack detailed understanding of the transfer mechanism between the soil solution and overland flow, particularly under field conditions. A model of solute transfer from soil solution to overland flow was developed. The model is based on the hypothesis that a solute is initially distributed uniformly throughout the soil pore space in a thin layer at the soil surface. A fundamental assumption of the model is that at the time runoff commences, any solute at the soil surface that could be transported into the soil with the infiltrating water will already have been convected away from the area of potential exchange. Solute remaining at the soil surface is therefore not subject to further infiltration and may be approximated as a layer of tracer on a plane impermeable surface. The model fitted experimental data very well in all but one trial. The model in its present form focuses on the exchange of solute between the soil solution and surface water after the commencement of runoff. Future model development requires the relationship between the mass transfer parameters of the model and the time to runoff: to be defined. This would enable the model to be used for extrapolation beyond the specific experimental results of this study. The close agreement between experimental results and model simulations shows that the simple transfer equation proposed in this study has promise for estimating solute loss to surface runoff. Copyright (C) 2000 John Wiley & Sons, Ltd.

Experimental examination of solute transport by surface runoff from low-angle slopes

The removal of chemicals in solution, by overland flow from agricultural land has the potential to be a significant source of chemical loss from zero-till and surface mulched farming systems. The objective of this study was to determine the magnitude of solute loss by surface runoff from agricultural systems. Previous experiments have enhanced the understanding of the exchange process, but the initial soil conditions together with the tracer application method in these experiments have meant that in some cases the results have limited applicability to field situations. In this study, two different sets of experiments were carried out to determine the magnitude of solute loss by surface runoff. These experiments entailed the surface application of bromide to (1) field scale plots 18 m long by 2 m wide and (2) repacked soil cores 236 mm in diameter; followed by the application of simulated rainfall in both cases. The most substantial finding of the field experiments was that the quantities of solute in surface runoff varied greatly with soil type and structure (0.07-14.9% of the applied bromide). Also, on some soils, large quantities of tracer were measured in the surface runoff even after several hours of infiltration. The experiments on soil cores showed that soil structure plays an important role in the quantity of chemical that may be transported in the surface runoff. These field results showed that, in certain systems, solute movement by overland flow is an important transport mechanism, which should be considered when budgeting for chemical loss. (C) 2000 Elsevier Science B.V. All rights reserved.

Identifying regional dust transport pathways: Application of kinematic trajectory modelling to a trans-Tasman case

PI kinematic trajectory model is used to investigate potential pathways of dust transport from Australia to New Zealand. Historically, these have been assumed to follow rather direct west-east trajectories spanning 2 to 3 days, often resulting in red snow events in the Southern Alps of New Zealand. However, results from the present study which examined the route taken by air parcels originating in southern Australia during dust storms on 24 and 25 May 1994, indicate that trans-Tasman dust transport trajectories are more diverse than previously thought, and display considerable variation during single events. These mon divergent pathways tie in more closely with aeolian dust sedimentation patterns identified by ocean coring in the Tasman Sea, and may account for the deposition of Australian dust on sub-Antarctic islands located well south of the Australian continent. Copyright 2000 John Wiley Sons, Ltd.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

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

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

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

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