Detailed characterization of solute transport in a heterogeneous field soil

There is a necessity for improved physical understanding of solute transport processes in heterogeneous soil systems. In situ nondestructive techniques like time domain reflectometry (TDR) and fiber optic miniprobes (FOMPs) permit the collection of unique measurements of solute transport processes in soils for the purposes of model development and validation. This study examined the application of TDR and FOMPs to measure solute transport at various points laterally and at two depths in a heterogeneous clay-loam soil. A miscible displacement experiment was performed at a constant irrigation flux to examine the applicability of these probes to field soils. In their first application to a field soil, the FOMPs were successfully calibrated and performed well in measuring solute breakthrough curves. Two flow regimes were identified in the soil profile, the first where lateral spreading of the solute occurred in the surface horizon, followed by convergence into preferential flow pathways in the second transport zone. The measured transport response was heterogeneous with at least two identifiable vertical flow phases. It was demonstrated using transfer function modeling and data from a corresponding laboratory study that the FOMPs were measuring the slower phase, while the TDR probes captured a composite of the fast and slow phases. The combination of these two techniques may be a means to separate solute transport phases in heterogeneous media and relate laboratory column results to field studies.     

Field-scale preferential transport of water and chloride tracer by depression-focused recharge

A tracer study was initiated in November 1993 to investigate depression-focused recharge and to monitor solute movement through the vadose zone into the shallow ground water in southeastern North Dakota. Granular potassium chloride (KCl) was surface-applied to two areas overlying subsurface drains and to one area instrumented with soil solution samplers, ground water monitoring wells, time domain reflectometry (TDR) probes, and temperature probes. One of the subsurface drain tracer plots was located on level ground while the other two sites were in small topographic depressions. Formation of ground water mounds beneath the depressions indicated that these areas are recharge sites. The applied Cl- tracer was found to move rapidly to the shallow ground water under the depressional areas after infiltration of spring snowmelt in 1994. Excessive rainfall events were also responsible for focused recharge and the rapid transport of the applied Cl- tracer. Water flow through partially frozen soil at the bottom of the depressions during thaw enhanced preferential movement of the tracer.     

Solute transport modeling under cultivated sandy soils and transient water regime

Drainable lysimeters offer the possibility to integrate heterogeneous solute leaching conditions caused by row crops and transient water regime, and to conveniently measure water and solute fluxes at the drainage outlet. To compare solute leaching behavior in and around drainable lysimeters operating under a transient water regime in potato (Solanum tuberosum L.) fields, parameters of the convective lognormal transfer (CLT) function model were fitted using bromide (Br-) flux concentrations (Cf) measured in lysimeters and from Br- resident concentrations (Cr) measured in adjacent soil cores. Expected mean values Ez(I) obtained from Cr and Cf CLT parameters were equivalent and well correlated (R2 = 0.78). However, estimated median values mu of the CLT function were smaller when derived from Cr (1.05 to 1.28) compared with Cf (1.23 to 2.14). Most mu values were also smaller than previously reported values for a 30-cm reference depth, indicating that 50% of solute mass would leach more readily in these coarse sandy soils. Higher variance and dispersion of Cr compared with those of Cf could be related to a smaller sampling support (sample size/sampling area) in the case of Cr measured by soil coring, or to disruption of solute transport mechanisms in the repacked lysimeter. Retained Br- in the top soil layer after 12 to 17 cm of cumulative drainage was indicated by measured Cr. Neither CLT function simulated well residual topsoil Cr values, indicating that Br- plant cycling or preferential flow probably interfered even though tuber Br- uptake was relatively small.     

Differences in microbial biomass,organic carbon,and dye sorption between flow and no-flow regions of unsaturated soil

Transport models in which the liquid phase is partitioned between conducting and nonconducting regions allow the possibility that degradation and sorption are different in these regions. However, there is little information on biological or chemical differences between conducting and nonconducting regions of the soil matrix. Previous work by the authors on Br transport through unsaturated, intact soil cores of Dundee silty clay loam (fine-silty, mixed, active, thermic Typic Endoaqualf) indicated non-equilibrium conditions that could be well-described by a two-region model. Fitted parameters indicated little solute transfer between flow regions, suggesting that dye movement in unsaturated soil might delineate conducting and nonconducting regions of this soil. Steady-state, unsaturated flow was established in intact cores (10 by 30 cm) of the Dundee soil, then Br and erioglaucine dye were displaced through these cores. The soil cores were then sectioned into 5-cm segments and stained soil was separated from unstained soil. Microbial biomass C, organic C, and dye sorption K(D) (= g(sorbed) kg(-1)soil/g L(-1)) values for stained and unstained soil were determined. Stained soil had higher microbial biomass C but generally lower organic C and lower affinity for dye sorption than unstained soil from the same depth increment. Fraction of immobile water, dispersion, and mass transfer between conducting and nonconducting regions were consistent with previous results.     

Impact of preferential flow at varying irrigation rates by quantifying mass fluxes

Solute concentration and soluble dye studies inferring that preferential flow accelerates field-scale contaminant transport are common but flux measurements quantifying its impact are essentially nonexistent. A tile-drain facility was used to determine the influence of matrix and preferential flow processes on the flux of mobile tracers subjected to different irrigation regimes (4.4 and 0.89 mm h(-1)) in a silt loam soil. After tile outflow reached steady state either bromide (Br; 280 kg ha(-1)) or pentafluorobenzoic acid (PFBA; 121 kg ha(-1)) was applied through the irrigation system inside a shed (3.5 x 24 m). Bromide fluxes were monitored at an irrigation rate of 4.4 mm h(-1) while PFBA fluxes were monitored at an irrigation rate of 0.89 mm h(-1). At 4.4 mm h(-1) nearly one-third of the surface-applied Br was recovered in the tile line after only 124 mm of irrigation and was poorly fit by the one-dimensional convective-dispersive equation (CDE). On the other hand, the one-dimensional CDE fit the main PFBA breakthrough pattern almost perfectly, suggesting the PFBA transport was dominated by matrix flow. Furthermore, after 225 mm of water had been applied, less than 2% of the applied PFBA had been leached through the soil compared with more than 59% of the applied Br. This study demonstrates that the methodology of applying a narrow strip of chemical to a tile drain facility is appropriate for quantifying chemical fluxes at the small-field scale and also suggests that there may be a critical input flux whereby preferential flow is initiated.     

Preferential flow of bromide, bentazon, and imidacloprid in a Dutch clay soil

Leaching to ground water and tile drains are important parts of the environmental assessment of pesticides. The aims of the present study were to (i) assess the significance of preferential flow for pesticide leaching under realistic worst-case conditions for Dutch agriculture (soil profile with thick clay layer and high rainfall) and (ii) collect a high-quality data set that is suitable for testing pesticide leaching models. The movement of water, bromide, and the pesticides bentazon [3-isopropyl-1H-2, 1,3-benzothiadiazine-4(3H)-one-2,2-dioxide] and imidacloprid [1-[(6-chloro-3-pyridinyl)-methyl]-N-nitro-2-imidazolidinimine] was monitored in a clay soil for about 1 yr. The 1.2-ha field was located in the central part of the Netherlands (51 degrees 53' N, 5 degrees 43' E). The soil was a Eutric Fluvisol cropped with winter wheat (Triticum aestivum L.). Tile drains were present at a 0.8- to 0.9-m depth and the ground water level fluctuated between a 0.5- and 2-m depth. All chemicals were applied in spring. None of the soil concentration profiles showed bimodal concentration distributions. However, for each substance the highest concentration in drain water was found in the first drainage event after its application, which indicates preferential flow. This preferential flow is probably caused by permanent macropores that were present in the 0.3- to 1.0-m layer. At the time of the first drainage event, the drain water concentration of each substance was about an order of magnitude higher than its ground water concentration. Thus, the flux concentrations in drain water proved to be a more sensitive detector of preferential flow than the resident concentrations in the soil profile and the ground water.     

Sorption and transport of 17beta-estradiol and testosterone in undisturbed soil columns

Land-applied domestic animal wastes contain appreciable amounts of 17beta-estradiol (henceforth, estradiol) and testosterone. These sex hormones may be transported through soil to groundwater and streams, where they may adversely affect the environment. Previous column transport studies with these hormones used repacked soil and did not consider preferential flow. We, therefore, determined the sorption and transport characteristics of estradiol and testosterone in undisturbed soil columns (15-cm i.d. by 32-cm height). In the sorption experiment, isotherms for estradiol and testosterone were nonlinear with Freundlich exponents (n) less than one. Sorption of both hormones decreased with soil depth, and estradiol sorbed more strongly than testosterone. Average estradiol Freundlich sorption coefficients (K(f)) values were 36.9 microg(1 - n) mL(n) g(-1) for the 0- to 10-cm soil depth and 25.7 microg(1 - n) mL(n) g(-1) for the 20- to 30-cm soil depth. Average testosterone K(f) values were 26.7 microg(1 - n) mL(n) g(-1) for the 0- to 10-cm soil depth and 14.0 microg(1 - n) mL(n) g(-1) for the 20- to 30-cm soil depth. In the transport experiment, 27% of the estradiol and 42% of the testosterone leached through the soil columns. Approximately 50% of the remaining soil-bound hormones were sorbed in the top 10 cm of soil. In almost all instances, breakthrough concentrations of estradiol, testosterone, and a chloride tracer peaked simultaneously. Simultaneous breakthrough and HYDRUS-1D transport parameters indicated both chemical and physical nonequilibrium processes affected hormone transport. This suggests hormones placed on soil surfaces may contaminate groundwater under conditions of preferential flow.     

Copper and calcium transport through an unsaturated soil column

To determine the relative importance of the physical and chemical factors that influence the movement of heavy metals through soils, leaching experiments were carried out under conditions of constant molarity during unsaturated steady-state water flow through a Manawatu fine sandy loam (a Dystric Fluventic Eutochrept). The movement and exchange of copper was studied in a binary Cu-Ca system. The movement of the associated anions, namely chloride and sulfate, was also monitored. The measurements were compared with predictions from the convection-dispersion equation (CDE), linked with cation exchange theory. The agreement between the measured and predicted breakthrough of sulfate and copper was good. This indicates that copper retardation in the Manawatu soil is closely related to the cation exchange capacity, and that exchange between Ca and Cu is the main process of Cu retardation in the Manawatu soil. However, copper appeared slightly later in the effluent than predicted, indicating that non-exchange processes are also involved in copper transport. Measurements of suction cups could also be used to obtain the parameters for the CDE to describe sulfate movement through the soil. Time domain reflectometry (TDR) measurements of the bulk-soil electrical conductivity could be used to monitor the movement of both sulfate and copper. This indicates that TDR can also be used to monitor cation transport and exchange through the soil, provided the percolating solution causes a sufficient change in the electrical conductivity.     

Phosphorus exchangeability and leaching losses from two grassland soils

Although phosphate phosphorus (P) is strongly sorbed in many soils, it may be quickly transported through the soil by preferential flow. Under flood irrigation, preferential flow is especially pronounced and associated solute losses may be important. Phosphorus losses induced by flood irrigation were investigated in a lysimeter study. Detailed soil chemical analyses revealed that P was very mobile in the topsoil, but the higher P-fixing capacity of the subsoil appeared to restrict P mobility. Application of a dye tracer enabled preferential flow pathways to be identified. Soil sampling according to dye staining patterns revealed that exchangeable P was significantly greater in preferential flow areas as compared with the unstained soil matrix. This could be partly attributed to the accumulation of organic carbon and P, together with enhanced leaching of Al- and Fe-oxides in the preferential flow areas, which resulted in reduced P sorption. The irrigation water caused a rapid hydrologic response by displacement of resident water from the subsoil. Despite the occurrence of preferential flow, most of the outflowing water was resident soil water and very low in P. In these soils the occurrence of preferential flow per se is not sufficient to cause large P losses even if the topsoil is rich in P. It appears that the P was retained in lower parts of the soil profile characterized by a very high P-fixing capacity. This study demonstrates the risks associated with assessing potential P losses on the basis of P mobility in the topsoil alone.     

Transport and degradation of toluene and o-xylene in an unsaturated soil with dipping sedimentary layers

A lysimeter trench was established at the Gardermoen delta (50 km north of Oslo, Norway) to study the flow of water and transport and degradation of aromatic jet fuel components (toluene and o-xylene) in the undisturbed unsaturated zone. Site investigations with ground-penetrating radar revealed the presence of dipping sedimentary layers within the foreset unit. This study has shown that the foreset bed of the Gardermoen delta structure provided a preferential flow path for the transport of the solute plumes, but did not have dramatic effects on the degradation potential under the current conditions. The degradation potential for toluene and o-xylene in the unsaturated zone at Gardermoen was very high and almost all of the injected hydrocarbons were biodegraded before reaching the saturated zone. However, the horizontal displacement of the plume showed that knowledge about sedimentary structures in the unsaturated zone is important for a sufficient monitoring of contaminant transport and for remediation purposes. Carbon dioxide and O2 were measured in situ simultaneously with extraction of water samples, and indicated aerobic biodegradation of toluene and o-xylene. Overall, first-order degradation coefficients were calculated to be in the range of 0.19 to 0.21 d(-1) and 0.10 to 0.11 d(-1) for toluene and o-xylene, respectively.     

Enhanced biogeochemical cycling and subsequent reduction of hydraulic conductivity associated with soil-layer interfaces in the vadose zone

Biogeochemical dynamics in the vadose zone are poorly understood due to the transient nature of chemical and hydrologic conditions but are nonetheless critical to understanding chemical fate and transport. This study explored the effects of a soil layer on linked geochemical, hydrological, and microbiological processes. Three laboratory soil columns were constructed: a homogenized medium-grained sand, a homogenized organic-rich loam, and a sand-over-loam layered column. Upward and downward infiltration of water was evaluated during experiments to simulate rising water table and rainfall events, respectively. In situ collocated probes measured soil water content, matric potential, and Eh. Water samples collected from the same locations were analyzed for Br, Cl, NO, SO, NH, Fe, and total sulfide. Compared with homogeneous columns, the presence of a soil layer altered the biogeochemistry and water flow of the system considerably. Enhanced biogeochemical cycling was observed in the layered column over the texturally homogeneous soil columns. Enumerations of iron- and sulfate-reducing bacteria showed 1 to 2 orders of magnitude greater community numbers in the layered column. Mineral and soil aggregate composites were most abundant near the soil-layer interface, the presence of which likely contributed to an observed order-of-magnitude decrease in hydraulic conductivity. These findings show that quantifying coupled hydrologic-biogeochemical processes occurring at small-scale soil interfaces is critical to accurately describing and predicting chemical changes at the larger system scale. These findings also provide justification for considering soil layering in contaminant fate and transport models because of its potential to increase biodegradation or to slow the rate of transport of contaminants.     

Soil water repellency induced by long-term irrigation with treated sewage effluent

This study describes soil water repellency developed under prolonged irrigation with treated sewage effluent in a semiarid environment. Soil surface layer (0-5 cm) and soil profile (0-50 cm) transects were sampled at a high resolution at the close of the irrigation season and rainy winter season. Samples from 0- to 5-cm transects were subdivided into 1-cm slices to obtain fine scale resolution of repellency and organic matter distribution. Extreme to severe soil water repellency in the 0- to 5-cm soil surface layer persisted throughout the 2-yr study period in the effluent-irrigated Shamouti orange [Citrus sinensis (L.) Osbeck cv. Shamouti] orchard plot. Nearby Shamouti orange plots irrigated with tap water were either nonrepellent or only somewhat repellent. Repellency was very variable spatially and with depth, appearing in vertically oriented "repellency tongues." Temporal and spatial variability in repellency in the uppermost 5-cm soil surface layer was not related to seasonality, soil moisture content, or soil organic matter content. Nonuniform distribution of soil moisture and fingered flow were observed in the soil profile after both seasons, demonstrating that the repellent layer had a persistent effect on water flow in the soil profile. A lack of correlation between bulk density and volumetric water content in the soil profile demonstrates that the observed nonuniform spatial distribution of moisture results from preferential flow and not heterogeneity in soil properties. Soil water repellency can adversely affect agricultural production, cause contamination of underlying ground water resources, and result in excessive runoff and soil erosion.     

Simulating sulfadimidine transport in surface runoff and soil at the microplot and field scale

To prevent residues of veterinary medicinal products (VMPs) from contaminating surface waters and ground water, an environmental impact assessment is required before a new product is allowed on the market. Physically based simulation models are advocated for the calculation of predicted environmental concentrations at higher tiers of the assessment process. However, the validation status of potentially useful models is poor for VMP transport. The objective of this study was to evaluate the dual-permeability model MACRO for simulation of transport of sulfonamide antibiotics in surface runoff and soil. Special focus was on effects of solute application in liquid manure, which may alter the hydraulic properties at the soil surface. To this end we used data from a microplot runoff experiment and a field experiment, both conducted on the same clay loam soil prone to preferential flow. Results showed that the model could accurately simulate concentrations of sulfadimidine and the nonreactive tracer bromide in runoff and in soil from the microplot experiments. The use of posterior parameter distributions from calibrations using the microplot data resulted in poor simulations for the field data of total sulfadimidine losses. The poor results may be due to surface runoff being instantly transferred off the field in the model, whereas in reality re-infiltration may occur. The effects of the manure application were reflected in smaller total and micropore hydraulic conductivities compared with the application in aqueous solution. These effects could easily be accounted for in regulatory modeling.     

Atrazine and alachlor transport in claypan soils as influenced by differential antecedent soil water content

Increased attention to ground water contamination has encouraged an interest in mechanisms of solute transport through soils. Few studies have investigated the effect of the initial soil water content on the transport and degradation of herbicides for claypan soils. We investigated the effect of claypan soils at initial field capacity vs. permanent wilting level on atrazine and alachlor transport. The soil studied was Mexico silt loam (fine, smectitic, mesic Aeric Vertic Epiaqualf) with a subsoil clay content, primarily montmorillonite, of >40%. Strontium bromide, atrazine, and alachlor were applied to plots; half were at field capacity (Wet treatment), and half were near the permanent wilting point (Dry treatment). Soil cores were removed at selected depths and times, and cores were analyzed for bromide and herbicide concentrations. Bromide, atrazine, and alachlor were detected at the 0.90-m depth in dry plots within 15 d after experiment initiation. Bromide was detected 0.15 m deeper (P < 0.05) in the Dry compared with the Wet treatment at 1, 7, and 60 d after application and >0.30 m deeper (P < 0.01) in the Dry treatment at 15 and 30 d after application; similar treatment results were found for atrazine and alachlor, although on fewer dates with significant differences. The mobility order of the applied chemicals was bromide > atrazine > alachlor. The atrazine apparent half-life was significantly longer in the Dry plots compared with the Wet plots. The retardation factor determined from the relative velocity of each herbicide to that of bromide was higher for alachlor than for atrazine. This study identifies the impact that shrinkage cracks have for different moisture conditions on preferential transport of herbicides in claypan soils.     

PERCOLATE WATER AND BROMIDE MOVEMENT IN THE ROOT ZONE OF EFFLUENT IRRIGATION SITES1

ABSTRACT: A bromide tracer was used to evaluate percolate water and ion movement in the upper 1.2 m of soil at a proposed sewage effluent irrigation site located in the Missouri Ozarks. Two plots representing Doniphan silt loam and Crider silt loam soils were sprinkler irrigated with local ground water at a rate of 7.62 cm/week from June through August 1976. Soil water potential, percent soil moisture by volume, and background levels of bromide in soil water, ground water, and precipitation were measured at the study plots. Bromide exchange properties and saturated hydraulic conductivity of the soils were determined in the laboratory. During two selected time periods, irrigation water, was spiked with NaBr (5.0 mg/l Br). Bromide movement through the upper profile was quantified by soil water samples and post-sampling neutron activation analysis. Soil moisture was near saturatin in both soils when the Br tracer was applied. Bromide concentrations above background levels (0.023 mg/l Br, Doniphan silt loam and 0.016 mg/l Br, Crider silt loam) were detected within 2.60 hours at 0.9 m in the Doniphan soil and within 3.75 hours at that depth in the Crider soil. The rate of Br movement in the profile was greater in both soils than the measured saturated hydraulic conductivity, Bromide concentrations above background levels were present in soil water from the study plots for a minimum of 21 days after irrigation with the Br tracer.     

Micro-trench experiments on interflow and lateral pesticide transport in a sloped soil in northern Thailand

During recent decades, a change in land use in the mountainous regions of Northern Thailand has been accompanied by an increased input of agrochemicals. We identified lateral water flow and pesticide transport pathways and mechanisms in a Hapludult on a sloped litchi orchard in Northern Thailand. During two rainy seasons, two micro-trench experiments were performed at the plot scale (2 by 3 m). The first experiment was performed at the footslope of the orchard; the second was performed at a midslope position. Two salt tracers (bromide and chloride) and two pesticides {methomyl [S-methyl-N-(methylcarbamoyloxy)thioacetimidate] and chlorothalonil (2,4,5,6-Tetrachlor-1,3-benzdicarbonitril)} were applied in stripes parallel to the slope 150 and 300 cm away from the trench. At the trench, soil water was collected by wick samplers. Tensiometers and time-domain reflectometry probes were installed. At the end of the experiment, soil samples were taken and analyzed for residual concentrations of tracers and pesticides. Lateral subsurface flow of water occurred exclusively along preferential flow paths and was mainly observed at 0- to 30- and 60- to 90-cm depth. Lateral transport of pesticides was negligible, but both pesticides were found beneath the application area at 90 cm depth. Therefore, they may pose a groundwater contamination risk. The amount of wick flow and the location of interflow were mainly a function of rain amount and antecedent soil water suction. During dry periods, water flow was restricted to the topsoil. After heavy rain events and wet periods, interflow was mainly observed in the subsoil. The cumulative rain amount between samplings necessary to induce interflow was 20 mm. At the footslope, the interflow was seven times higher, and the network of water-bearing pores increased compared with the midslope position.     

Phosphorus leaching in relation to soil type and soil phosphorus content

Phosphorus losses from arable soils contribute to eutrophication of freshwater systems. In addition to losses through surface runoff, leaching has lately gained increased attention as an important P transport pathway. Increased P levels in arable soils have highlighted the necessity of establishing a relationship between actual P leaching and soil P levels. In this study, we measured leaching of total phosphorus (TP) and dissolved reactive phosphorus (DRP) during three years in undisturbed soil columns of five soils. The soils were collected at sites, established between 1957 and 1966, included in a long-term Swedish fertility experiment with four P fertilization levels at each site. Total P losses varied between 0.03 and 1.09 kg ha(-1) yr(-1), but no general correlation could be found between P concentrations and soil test P (Olsen P and phosphorus content in ammonium lactate extract [P-AL]) or P sorption indices (single-point phosphorus sorption index [PSI] and P sorption saturation) of the topsoil. Instead, water transport mechanism through the soil and subsoil properties seemed to be more important for P leaching than soil test P value in the topsoil. In one soil, where preferential flow was the dominant water transport pathway, water and P bypassed the high sorption capacity of the subsoil, resulting in high losses. On the other hand, P leaching from some soils was low in spite of high P applications due to high P sorption capacity in the subsoil. Therefore, site-specific factors may serve as indicators for P leaching losses, but a single, general indicator for all soil types was not found in this study.     

SOIL MOISTURE SENSORS FOR CONTINUOUS MONITORING1

ABSTRACT: Certain physical and chemical properties of soil vary with soil water content. The relationship between these properties and water content is complex and involves both the pore structure and constituents of the soil solution. One of the most economical techniques to quantify soil water content involves the measurement of electrical resistance of a dielectric medium that is in equilibrium with the soil water content. The objective of this research was to test the reliability and accuracy of fiberglass soil-moisture electrical resistance sensors (ERS) as compared to gravimetric sampling and Time Domain Reflectometry (TDR). The response of the ERS was compared to gravimetric measurements at eight locations on the USDA-ABS Walnut Gulch Experimental Watershed. The comparisons with TDR sensors were made at three additional locations on the same watershed. The high soil rock content (≥45 percent) at seven locations resulted in consistent overestimation of soil water content by the ERS method. Where rock content was less than 10 percent, estimation of soil water was within 5 percent of the gravimetric soil water content. New methodology to calibrate the ERS sensors for rocky soils will need to be developed before soil water content values can be determined with these sensors.     

Autonomous real-time adaptive management of soil salinity using a receding horizon control algorithm: a pilot-scale demonstration

Soil salinization is a potentially negative side effect of irrigation with reclaimed water. While optimization schemes have been applied to soil salinity control, these have typically failed to take advantage of real-time sensor feedback. This study incorporates current soil observation technologies into the optimal feedback-control scheme known as Receding Horizon Control (RHC) to enable successful autonomous control of soil salinization. RHC uses real-time sensor measurements, physically-based state prediction models, and optimization algorithms to drive field conditions to a desired environmental state by manipulating application rate or irrigation duration/frequency. A simulation model including the Richards equation coupled to energy and solute transport equations is employed as a state estimator. Vertical multi-sensor arrays installed in the soil provide initial conditions and continuous feedback to the control scheme. An optimization algorithm determines the optimal irrigation rate or frequency subject to imposed constraints protective of soil salinization. A small-scale field test demonstrates that the RHC scheme is capable of autonomously maintaining specified salt levels at a prescribed soil depth. This finding suggests that, given an adequately structured and trained simulation model, sensor networks, and optimization algorithms can be integrated using RHC to autonomously achieve water reuse and agricultural objectives while managing soil salinization.     

WAVE INDUCED FLOW AND TRANSPORT IN SEDIMENT BEDS1

ABSTRACT: The fate of contaminants in large water bodies is highly influenced by the transfer of flow and solutes across the water sediment interface. In this paper, an analytical model is presented where flow in both sediment bed and open channel is coupled at the interface through a boundary layer occupying the upper part of the sediment bed. The presence of this layer allows not only the capture of the inertia effects through a drag term in the generalized Darcy's equation, but also the specification of different soil parameters for the two porous zones. The flow is advective and driven by wave action along the water surface. The resulting system is solved for the pressure and flux in each sediment layer. The generated transport velocity fields are linked to a random walk simulation that is used to examine the trajectories of solute particles. Comparison of these trajectories against experimental tracer tests suggests a pattern very similar to the one attributed to the presence of surface mounds. The results clearly show the significance of the boundary layer and the drag term for soil with high permeability and the impact of both the thickness of the boundary layer and the length of the gravity wave relative to the depth of the water channel on the transport and exchange across the interface. The paper also examines the sensitivity of the mass exchange to the permeability of the two porous zones.