|
|
||||||||
a Dep. of Soil Science, Madison, WI USA
b Dep. of Biol. System Eng., Univ. of Wisconsin, Madison, WI 53706-1299 USA
c Dep. of Agronomy, Purdue Univ., West Lafayette, IN USA
d Dep. of Agric. and Biol. Eng., Cornell Univ., Ithaca, NY 14850 USA
e Weed Science Lab., USDA-ARS, BARC-W, 10300 Baltimore Blvd., Beltsville, MD 20705-2350 USA
kung{at}calshp.cals.wisc.edu
| ABSTRACT |
|---|
|
|
|---|
4 h after irrigation started, yet benzoic acids, PFBA, o-TFMBA, and 2,6-DFBA, were detected in the tile drainage at 102 min, 42 min, and 18 min after their applications, respectively. Tracer mass recovery from tile drainage was Br- (7.04%), PFBA (13.9%), o-TFMBA, (18.7%), and 2,6-DFBA (19.7%) of applied mass. The faster arrival time and greater recovery of sequentially applied tracers confirmed that water movement and contaminant transport shifts toward increasingly larger pores of the preferential flow paths as soil becomes wet during a precipitation event. The breakthrough patterns of these tracers can be used to quantify the water flux distributions of preferential paths. Because
90% of the chemical leached from this precipitation event occurred during the first day, it was critical to intensively monitor contaminant transport during the first 24 h after a rainfall. A soil sampling protocol based on collecting soil cores at random locations once every several days is unsuitable for determining the deep leaching under field conditions.
Abbreviations: BTC, breakthrough curve HPLC, high performance liquid chromatography IC, ion chromatography o-TFMBA, o-trifluoromethylbenzoic acid PFBA, pentafluorobenzoic acid WT, water tracer 2,6-DFBA, 2,6-difluorobenzoic acid
| INTRODUCTION |
|---|
|
|
|---|
Currently, efforts to directly measure soil hydraulic properties have been gradually decreasing because field results show that, with limited budget and instruments, the soil hydraulic properties and their spatial and temporal variability are too complex to be grasped, and field-scale water flux distribution cannot be directly and adequately characterized. Furthermore, results from field experiments demonstrated that strongly adsorbing chemicals can move as fast and deep into a soil profile as more mobile chemicals. For example, using a tile drain as an alternative sampling protocol, Everts et al. (1989) observed that both non-adsorbing Br- and adsorbing rhodamine water tracer (WT) dye had fast breakthrough patterns within the first 4 h after the initiation of the irrigation. Kladivko et al. (1991, 1999) demonstrated that pesticides were most susceptible to leaching during rainstorms that occur shortly after chemical application in spring, and more adsorbing pesticides atrazine [2-chloro-4-(ethylamino)-6-(isopropylamino)-s- triazine] and alachlor [2-chloro-2',6'-diethyl-N-(methoxymethyl)-acetanilide] can be transported to shallow groundwater (between 0.81.0 m) as rapidly as more mobile chemicals such as carbofuran (2,3-dihydro-2,2-dimethyl-7-benzofuranol methylcarbamate) during this period. Kung et al. (2000) applied a pulse of chloride and rhodamine WT dye solution in a field experiment on a clay loam soil. They found that both tracers arrived at the tile after 13 min after application. Saturated soil hydraulic conductivity of the field, based on soil cores collected at different depths, ranged from 0.5 to 5.0 x 10-6 m s-1. According to the conventional conceptualization, neither the conservative nor the reactive tracers should have reached the tile drain within 2 d.
The term preferential flow paths was coined to describe the fast contaminant transport that bypassed the vast soil matrix (Kung, 1990). It is critical to comprehend preferential flow in order to develop conceptual and numerical models to predict contaminant transport in unsaturated soils. Ju and Kung (1997) conducted numerical experiments in sandy soils with layering structures to characterize the water dynamics of the funnel-type preferential flow pathways and their impact on contaminant transport in an unsaturated sandy soil profile. They justified their numerical approach by the following two reasons: (i) the impact of a funnel mechanism that is caused by layered structure on contaminant transport was well understood (Kung, 1993; Ju and Kung, 1993); and (ii) the multiple-layered soil structures in an unsaturated sandy profile that trigger funnel flow are permanent and can be characterized by the ground-penetrating radar technique (Kung and Lu, 1993; Casper and Kung, 1996). In contrast, in unsaturated soils with finer textures, field-scale preferential flow mechanisms are not expressly understood. The preferential flow paths have a dynamic formationdestruction cycle, and there are no instruments to accurately map or reveal these paths. The preferential flow paths are three-dimensional, occupy only a very small volumetric fraction of the soil matrix, and cannot be quantified by the conventional techniques designed to measure the hydraulic properties of the bulk soil matrix. Therefore, the impacts of preferential flow paths on contaminant transport in fine-textured soils still cannot be numerically examined.
Jury (1982) introduced the transfer function concept to bypass the necessity of accurately and directly measuring the soil hydraulic properties of unsaturated soils in order to predict contaminant transport. In his approach, the breakthrough curve (BTC) of a conservative tracer was used to reflect the convoluted impact of the water dynamics on solute transport under field conditions. Jury's concept was adopted in this study; i.e., BTCs of conservative tracers were used to quantify the impact of the water dynamics of preferential flow on solute transport. In Jury's approach, the breakthrough of a single pulse of a conservative tracer was used and it was implicitly assumed that the water flux distribution was a stationary property of a soil profile. The rainfall intensityduration from individual precipitation events were not considered as parameters. Instead, only the total net infiltration from all precipitation events during a period was used to drive the contaminant transport.
Experimental results from Kung et al. (2000) showed that as a soil profile became wetter during a precipitation event, the arrival time of both conservative and reactive tracers to a tile drain shortened. This indicated that the water flux distribution of preferential flow paths change during a precipitation event; i.e., some preferential flow paths were not hydraulically active when a soil profile was initially dry. This implied that the water flux distribution of preferential flow paths was not a stationary property of a soil profile but depended primarily on soil water content of the soil profile. As a result, each rainfall event can have a different impact on the contaminant transport through preferential paths. In this study, different conservative tracers were sequentially applied during a single precipitation event. As the soil surface became wetter, the breakthrough patterns of tracers applied at different times reflected the impact of the change of the pore spectrum of hydraulically active preferential flow paths on contaminant transport.
The success of using a chemical BTC to quantify the water dynamics of preferential flow paths hinged on whether a representative breakthrough pattern could be measured to accurately reflect the convoluted impact of water fluxes along preferential flow paths on contaminant transport. As demonstrated by Ghodrati and Jury (1990), Ju et al. (1997), and Kung et al. (2000), conventional sampling protocols (e.g., soil coring and soil lysimeter methods) are unsuitable to monitor the field-scale chemical leaching patterns through preferential flow paths. Richards and Steenhuis (1988), Kladivko et al. (1991, 1999), and Czapar et al. (1994) showed that tile drain facilities can be used as an alternative.
Multiple conservative tracers with similar properties must first be chosen in order to reflect the impact of preferential flow paths on contaminant transport as a soil profile becomes wetter. Currently, Br- and Cl- are regarded as the most suitable conservative tracers and are commonly used to indicate water movement and solute transport. Bowman and Gibbens (1992) and Jaynes (1994) explored the transport properties of several fluorobenzoic acids, comparing them with Br-. Among the tested fluorobenzoic acid isomers, three were recommended (because they had the least adsorption and degradation) as alternative conservative tracers: PFBA, 2,6-DFBA, and o-TFMBA. The breakthrough patterns of these alternative tracers were rigorously tested in laboratory soil column studies (Jaynes, 1994). Under field conditions, only their recoveries in soil matrix from soil cores were compared, and the deep leaching through preferential paths could not be assessed by the soil coring method.
In this study, the primary objective was to use a tile drain monitoring facility to quantify the water dynamics of the preferential flow pathways and their impact on contaminant transport in a fine-textured, unsaturated soil profile. A preliminary experiment was first conducted to determine whether the three fluorobenzoic acids mentioned above would have the same breakthrough patterns as those of Br- and Cl- under field conditions in a soil profile with preferential flow paths. Although not a conservative tracer, nitrate has recently been blamed as a major cause of oxygen-depleted water (hypoxic zone) in the Gulf of Mexico and the USA. U.S. Geological Survey data indicated that >70% of the total N delivered to the Gulf by the Mississippi River originates above the confluence of the Ohio and Mississippi Rivers (Alexander et al., 1995; USDA-NRCS, 1996). A secondary objective of the preliminary study was to determine whether Br-, Cl-, and these fluorobenzoic acids have the same breakthrough patterns as freshly applied nitrate in a soil profile with preferential flow paths. After examining the suitability of three alternative tracers, the objective of our main experiment was then to quantify the water dynamics of preferential flow paths and their impact on contaminant transport in an unsaturated silt loam soil.
| Materials and methods |
|---|
|
|
|---|
1% soil organic matter in the top 45-cm depth and is typical of the southern portion of very extensive Midwest Corn Belt soils, extending from southwestern Ohio, across southern Indiana and southern Illinois, through Missouri, into southern Iowa and eastern Kansas (Kladivko et al., 1991). The perched groundwater at this site is maintained at approximately 0.75- to 0.9-m depth by tile drainage.
The five conservative tracers Br-, Cl-, PFBA, 2,6-DFBA, o-TFMBA, as well as nitrate, were co-applied to determine whether they had similar breakthrough patterns. The total masses of tracers applied were Br- (2.0 kg), Cl- and nitrate (2.4 kg each), and the benzoic acids (0.45 kg each). Each chemical was first dissolved separately (pH of each benzoic acid solution was adjusted to 7 by slowly adding KOH solution), and then all were combined into a
115-L tank mix to be applied uniformly through a sprayer. In most published results from field experiments, chemicals were applied to an entire field. However, in this preliminary study, the six tracers were applied only to a 1.5-m by 24-m strip. The long side of the strip was parallel to the long side of the plot and offset 0.3 m from the center tile line. The down-slope end of the chemical strip was located 33 m from the lower end of the field. The main reason for applying chemicals only to a strip was to reduce cost, because the benzoic acids were too expensive to be used throughout the entire field. This preliminary study was only to determine whether the benzoic acids had similar BTCs to those of Br-, Cl-, and nitrate. Also, because the watertable near the tile line has the steepest gradient, applying chemicals to a narrow strip close to the tile line could minimize the time required for the lateral transport of chemicals within a saturated zone to the tile line, so that the travel time of contaminants through preferential flow paths within the unsaturated soil profile can be accurately estimated.
In order to bring the watertable up to the tile line, the site had been pre-irrigated with
15 mm water (6.2 mm h-1 intensity) one day before the experiment with a solid-set sprinkler system. Then on the morning of the experiment, a 1-h pre-irrigation added another 6.2 mm water. After waiting for an hour, tracers were applied. Shortly after the tracer application, irrigation with the same intensity was applied for 35 min to leach these chemicals into the soil profile. The chemical application strip was located in the center of the irrigated area (24 m x 60 m). Again, the long sides of the irrigated area and the plot were parallel. Because of the plot's configuration (border tiles 10 m offset from the center tile), a 24-m-wide irrigation width extended beyond the adjacent border tiles and ensured that leachate originated from the treated chemical strip would only enter the center tile being monitored. The irrigation water was from the local municipal water supply and had none of the applied chemicals, except for 2.5 mg L-1 Cl-.
After the chemical application and initial irrigation, there were a 15-mm rainfall (Day 1), a 6.2-mm irrigation (Day 3), 11.7-mm of precipitation (Day 5), and 5.8-mm of precipitation (Day 8). The tile flow was continuously measured by the tipping-bucket method. Water samples were collected from the tile flow to determine the mass flux of these chemicals. The frequency of water sample collection was determined by the flow rate of the tile drainage. For example, shortly after each irrigation or precipitation event, when tile flow was relatively high, a 100-mL sample was collected every 10 min and three such samples were composited into a final water sample. As tile flow started to decrease, 100-mL samples were collected at 20-min intervals and three samples were composited into a final water sample. About two days after each irrigation or precipitation event, a 50-mL sample was collected at 30-min intervals and six samples were composited.
Main Experiments of Sequentially Applied Conservative Tracers
Field experiments were conducted at another 20-m by 225-m no-till plot at the South East Purdue Agricultural Center in Butlerville, Indiana, in May of 1997, when corn was starting to emerge. The tile flow rate was
1 mL s-1 (caused by a 1.2-cm pre-irrigation event 2 d earlier) and there was no precipitation in the five-day weather forecast. Irrigation (3 mm h-1 intensity) was continuously applied for 10 h to simulate a prolonged mild spring shower by a solid-set sprinkler system. The sprinkler system was made of two lines that were 12 m apart. Within each line, the distance between two adjacent nozzles was 12 m. Nozzles between the two lines were staggered by 6 m. The low intensity was chosen to avoid any surface runoff so that mass recovery could be accurately calculated. The irrigated area (24 m x 60 m) was about 73 m from the lower end of the field, where the tile monitoring facility was installed.
The conservative tracers tested in the preliminary study were sequentially applied. Shortly before the irrigation started, Br- (1.93 kg in 20 L water) was uniformly sprayed on a 1.5-m by 24-m strip. The strip was positioned at the center of the irrigated area and the long side of the strip was parallel to the tile line with 0.3 m offset. At 2 h after irrigation commenced, PFBA (1.37 kg in 16 L water) was sprayed on the same strip. Then, at 4 and 6 h after the irrigation started, o-TFMBA (1.5 kg in 16 L water) and 2,6-DFBA (1.5 kg in 20 L water) were sprayed on the same strip, respectively. For each acid tracer solution, pH was adjusted to 7 by slowly adding KOH solution during mixing. Each tracer was applied through a sprayer within
5 min; irrigation was continuous when tracers were applied. This sequential tracer application scheme was intended to explore how water and tracers would move through the larger end of the pore spectrum of preferential flow paths as a soil profile became wetter.
During the first 12 h after the irrigation started, water samples were collected from the tile once every 6 min. From 12 h to 24 h, water samples were collected at 10-min intervals and sequentially composited into 30-min samples. From the second to the fourth days, samples were collected every 40 min, with three samples composited to give 2-h samples. This intense sampling frequency was designed to register the fast contaminant breakthrough observed by Kung et al. (2000). The samples were stored at 5°C until analyzed, except during a short period of shipping. Flow rate of the tile drain was continuously measured by using the tipping-bucket method.
Chemical Analyses
A 3-mL aliquot from each sample was filtered (0.2 µm) and then analyzed by high performance liquid chromatography (HPLC) for the benzoic acids and by ion chromatography (IC) for Br-, Cl-, and nitrate. The HPLC conditions were as follows: mobile-phase, acetonnitrile 15 mM KH2PO4 (titrated to pH 2.6 with phosphoric acid) (35:65, v/v); flow rate, 1.7 mL min-1; guard column, Spherisorb SAX (10 mm x 4.6 mm x 5-µm i.d.; Waters, Milford, MA); analytical column, Supelcosil SAX1 (25 cm x 4.6 mm x 5-µm i.d.; Sigma-Aldrich, St. Louis); injection volume, 50 µL; and UV detection, 205 nm. The detection limit for each benzoic acid was 50 µg L-1. The IC conditions were: mobile-phase, 1 mM C8H5KO4 (titrated to pH 2.6 with NaOH); guard column, Waters IC-Pak anion; analytical column, IC-Pak anion (5 cm x 4.6 mm x 5-µm i.d.); and injection volume, 50 µL. The detection limit for of Br-, Cl-, and nitrate were 0.25, 0.5, 0.5 mg L-1, respectively.
| Results and discussion |
|---|
|
|
|---|
|
Main Experiment: Sequentially Applied Tracers
Based on 55 rain gauges, the measured total water applied during the 10-h irrigation was 28.8 mm with the Christianson uniformity of 0.85. The measured breakthrough patterns of the four applied tracers and the tile hydrograph are shown in Fig. 2
. The tile hydrograph showed that tile flow started to increase at 2 h and 50 min after irrigation started. Bromide was applied shortly before the irrigation started, while the PFBA was applied 2 h after the irrigation started. According to the conventional conceptualization of Darcian flow and convectiondispersion transport, Br- should reach the shallow groundwater first (i.e., first in, first out); however, PFBA was detected in the tile drainage at 3 h and 42 min after the irrigation started (i.e., 102 min after its application), while Br- was not detected until 4 h after its application. On the other hand, the o-TFMBA was applied at 4 h after the irrigation and arrived at the tile 4 h 42 min after the irrigation started (i.e., only 42 min after its application). The fourth tracer, 2,6-DFBA, was applied at 6 h after the irrigation began and arrived at the tile at only 18 min after its application. Figure 2 shows that the mass flux of the 2,6-DFBA at the first detection is much higher than those of the other three tracers. The water samples were collected once every 6 min in the first 12 h. It was very possible that the 2,6-DFBA actually first appeared in the tile drain shortly after the water sample taken at 12 min after application.
|
A preferential flow path is defined as hydraulically active when water and water-borne contaminants can directly reach the shallow groundwater surface through the preferential flow pathway. Let us assume there is a preferential flow path labeled "A". The tracers applied earlier entered into the Pathway A and were later sucked into small matrix pores along wall of the Pathway A. On the other hand, let us assume that the tracers applied later entered into the identical preferential flow Pathway A and reached groundwater through the Pathway A. Although both tracers passed through certain identical portions of flow Pathway A, the tracers should be considered to have traveled through two different flow paths. Because Br- from the first tracer pulse could not directly reach the shallow groundwater surface through the Pathway A, the preferential flow Pathway A was not hydraulically active to transport Br-. Alternatively, it can be said that Br- had traveled through a preferential flow Pathway B, which was made of two parts. The first part was made of the upper part of preferential flow Pathway A, while the lower part was made of some matrix pores. Therefore, the Pathway B had a smaller overall equivalent pore size than that of Pathway A. Later, when the soil profile became wet enough, the preferential flow Pathway A became completely hydraulically active to the fourth tracer. In other words, it can be described that preferential flow paths with larger equivalent pore sizes became hydraulically active later, when a soil profile became wetter.
The mechanism in which applied water and tracers would enter pores of a wide size-spectrum was not considered as a parameter in the conventional conceptualization of Darcian flow and convectivedispersive transport. As a result, a tracer applied early was expected to arrive at the tile first. The fact that Br- had the slowest arrival time indicated that the larger end of the pore spectrum of the preferential flow paths were not hydraulically active when the soil surface was relatively dry. Therefore, it is critical to consider how and when contaminants enter preferential flow paths with different pore sizes. The overall results from four tracers indicated that, as a soil profile becomes wetter during precipitation, water movement and contaminant transport would continuously shift toward the larger end of the pore spectrum of the preferential flow paths.
The tracer recovery (percentage of total applied mass) from tile drainage during the 100-h period were Br- (7.04%), PFBA (13.9%), o-TFMBA, (18.7%), and 2,6-DFBA (19.7%). The mass leached increased among the sequentially applied tracers; i.e., the later a tracer was applied, the more mass was leached. The fact that Br- had the least mass leached and the 2,6-DFBA had the most was consistent with their arrival times; i.e., Br- arrived the latest and the 2,6-DFBA arrived the earliest after application. Because Br- was applied to a relatively dry soil surface before the irrigation started, it had the least chance of being transported through the preferential flow paths. As a result, most of the applied Br- entered into the smaller end of the pore spectrum, was transported more slowly, and still remained in the soil matrix after 100 h. The mass recovery patterns confirmed that, when a soil profile became wetter during a prolonged precipitation event, contaminant transport would be continuously shifted toward the larger end of the pore spectrum of the preferential flow paths. The primary reason that mass recovery of Br- from this study was higher than that observed in the preliminary study was because the rain duration was much shorter in the preliminary study. In a precipitation event with very short duration, matric potential gradient would induce a tracer to enter into and be transported through the soil matrix pores.
These high recoveries after such a mild irrigation with only 28.8 mm total water applied were surprising. Note that in the preliminary experiment where
39 mm water was applied in a 10-d period, only
2.5% of the applied mass of each tracer was leached out. This was because the rain duration was short during the preliminary experiment. As a result, the tracers applied in the preliminary experiment had a better chance to enter the smaller matrix pores and became less susceptible to deep leaching through preferential flow paths.
Figure 2 shows that vast majority of the tracer recovered from the tile drain during the 100-h sampling period leached out from the root zone during the first day after the irrigation started (i.e., Br-, 90%; PFBA, 87%; o-TFMBA, 89%; and 2,6-DFBA, 91%), which is also when most of the water flow occurred. According to the comparison between the BTC of adsorbing and non-adsorbing tracers, Kung et al. (2000) found that the tail of the adsorbing tracer BTC trailed much faster than those of the non-adsorbing tracers. In other words, an adsorbing chemical had more of its leaching occur during the first day than a non-adsorbing chemical.
Among the deterministic models which simulate fast solute leaching through soil structural pores, the two-region (or two-domain) mobileimmobile conceptualization (Skopp and Warrick, 1974) was first proposed to explain the fast leaching observed in aggregated soils under laboratory conditions. Gerke and van Genuchten (1993) proposed the dual porosity model, in which water and solutes moved through two distinct and interacting domains simultaneously. Similarly, Gwo et al. (1995) proposed a multiple-porosity model, where soil is partitioned into many interacting domains. The arrival times and mass recovery patterns from our sequentially applied tracers and those from Kung et al. (2000) showed that the spectrum of hydraulically active preferential flow paths expanded continuously during a precipitation event; i.e., preferential flow paths with larger pore sizes became hydraulically active and contributed to contaminant transport as a soil profile became wetter. Because 90% or more of the contaminants that leached from the root zone in the 100-h sampling period were leached during the first 24 h through preferential flow paths, the contribution of the so-called immobile, or matrix flow, or micropore domain to contaminant transport was almost irrelevant during this time frame. Whether those contaminants that entered into soil matrix pores could move back into preferential flow paths and reach the watertable during subsequent precipitation events was beyond the scope of the current study.
One can use the BTCs of the sequentially applied tracers to make inferences about the behavior of freshly applied nitrate fertilizers, since the short time scale of the vast leaching (i.e., 1 d) would not likely provide time for significant N transformations and uptake to occur. Our data suggested that, after the application of nitrate fertilizer in spring to a relatively dry soil surface, as much as 7% of nitrate could be leached out of a root zone by a single 10-h precipitation event as mild as 3 mm h-1. On the other hand, if the nitrate were applied to relatively wet soil surface, as much as 20% could be rapidly leached out from a root zone after a mild precipitation event. This strongly suggested that highly soluble agrichemicals such as nitrate should never be applied to a wet soil surface. The question then becomes how wet is "wet"? In most cases it is impractical for N fertilizer to be applied in the middle of a rainstorm or immediately following 1.8 cm rain, as was 2,6-DFBA, but how long a delay would be required before the chemical behavior matches that of the Br- in this experiment is unknown. The current experiment cannot answer that question. The data do clearly indicate that fertigation of nitrate on wet soils could lead to substantial losses of compounds through preferential flow.
The fast arrival time and the vast transport during the first day after precipitation strongly suggested that, in order to collect representative samples to assess total deep leaching, it was extremely critical to intensively monitor contaminant transport during the first day after a rainfall event. Soil sampling protocols based on the coring method of collecting soil samples at random locations once every several days is inadequate to determine contaminant leaching under field conditions. On the other hand, a field where the tiles are installed with even spacing and depth and border tiles are installed to partition the water input area would behave like a huge, enclosed and undisturbed lysimeter. Water and contaminants leached from such fields can be accurately and holistically assessed, regardless of how rapidly and preferentially the contaminants are being transported downward.
The mass flux of each tracer plotted against time since application is shown in Fig. 3 (for Br-, time coincides with the onset of irrigation). The results show that the BTCs shifted toward the left when a tracer was applied to wetter soil during a precipitation event. Jury (1982) proposed to derive the travel-time probability density function from the tracer breakthrough patterns. The physical factors that determine the travel time under field conditions encompass the size spectrum of hydraulically active pores, chemical retardation and transformation, uptake, and degradation. Results from this study showed that much of the mass that would be leached out from a root zone was leached through preferential flow paths within the first 24 h of a precipitation. Furthermore, Kung et al. (2000) showed that the initial transport of an adsorbing tracer was almost identical to that of a non-adsorbing tracer. These results together suggested that chemical retardation, transformation, uptake, and degradation had only secondary effects in determining the initial deep leaching of agrichemicals through preferential flow paths. Under this condition, the travel-time probability density function is dictated by the field-scale size spectrum of hydraulically active pores of the preferential flow paths.
|
| Conclusions |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
Received for publication February 2, 1999.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
H. H. Gerke, J. Dusek, T. Vogel, and J. M. Kohne Two-Dimensional Dual-Permeability Analyses of a Bromide Tracer Experiment on a Tile-Drained Field Vadose Zone J., August 23, 2007; 6(3): 651 - 667. [Abstract] [Full Text] [PDF] |
||||
![]() |
K.-J.S. Kung, E. J. Kladivko, C. S. Helling, T. J. Gish, T. S. Steenhuis, and D. B. Jaynes Quantifying the Pore Size Spectrum of Macropore-Type Preferential Pathways under Transient Flow Vadose Zone J., August 24, 2006; 5(3): 978 - 989. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. W. Stone and J. T. Wilson Preferential Flow Estimates to an Agricultural Tile Drain with Implications for Glyphosate Transport J. Environ. Qual., August 9, 2006; 35(5): 1825 - 1835. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Alletto, Y. Coquet, P. Vachier, and C. Labat Hydraulic Conductivity, Immobile Water Content, and Exchange Coefficient in Three Soil Profiles Soil Sci. Soc. Am. J., June 21, 2006; 70(4): 1272 - 1280. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Gaur, R. Horton, D. B. Jaynes, and J. L. Baker Measured and Predicted Solute Transport in a Tile Drained Field Soil Sci. Soc. Am. J., April 19, 2006; 70(3): 872 - 881. [Abstract] [Full Text] [PDF] |
||||
![]() |
K.-J. S. Kung, M. Hanke, C. S. Helling, E. J. Kladivko, T. J. Gish, T. S. Steenhuis, and D. B. Jaynes Quantifying Pore-Size Spectrum of Macropore-Type Preferential Pathways Soil Sci. Soc. Am. J., June 28, 2005; 69(4): 1196 - 1208. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Kohne and H. H. Gerke Spatial and Temporal Dynamics of Preferential Bromide Movement towards a Tile Drain Vadose Zone J., February 1, 2005; 4(1): 79 - 88. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. J. Gish, C. L. Walthall, C. S. T. Daughtry, and K.-J. S. Kung Using Soil Moisture and Spatial Yield Patterns to Identify Subsurface Flow Pathways J. Environ. Qual., January 1, 2005; 34(1): 274 - 286. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. W. Haws and P. S. C. Rao The Effect of Vertically Decreasing Macropore Fractions on Simulations of Non-Equilibrium Solute Transport Vadose Zone J., November 1, 2004; 3(4): 1300 - 1308. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. J. Kladivko, J. R. Frankenberger, D. B. Jaynes, D. W. Meek, B. J. Jenkinson, and N. R. Fausey Nitrate Leaching to Subsurface Drains as Affected by Drain Spacing and Changes in Crop Production System J. Environ. Qual., September 1, 2004; 33(5): 1803 - 1813. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. J. Gish, K.-J. S. Kung, D. C. Perry, J. Posner, G. Bubenzer, C. S. Helling, E. J. Kladivko, and T. S. Steenhuis Impact of Preferential Flow at Varying Irrigation Rates by Quantifying Mass Fluxes J. Environ. Qual., May 1, 2004; 33(3): 1033 - 1040. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. N. Pace, M. N. Pace, M. A. Mayes, P. M. Jardine, T. L. Mehlhorn, J. M. Zachara, and B. N. Bjornstad Quantifying the Effects of Small-Scale Heterogeneities on Flow and Transport in Undisturbed Cores from the Hanford Formation Vadose Zone J., November 1, 2003; 2(4): 664 - 676. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. J. Gish, W. P. Dulaney, K.-J. S. Kung, C. S. T. Daughtry, J. A. Doolittle, and P. T. Miller Evaluating Use of Ground-Penetrating Radar for Identifying Subsurface Flow Pathways Soil Sci. Soc. Am. J., September 1, 2002; 66(5): 1620 - 1629. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. B. Jaynes, S. I. Ahmed, K.-J. S. Kung, and R. S. Kanwar Temporal Dynamics of Preferential Flow to a Subsurface Drain Soil Sci. Soc. Am. J., September 1, 2001; 65(5): 1368 - 1376. [Abstract] [Full Text] [PDF] |
||||
![]() |
K.-J.S. Kung, T.S. Steenhuis, E.J. Kladivko, T.J. Gish, G. Bubenzer, and C.S. Helling Impact of Preferential Flow on the Transport of Adsorbing and Non-Adsorbing Tracers Soil Sci. Soc. Am. J., July 1, 2000; 64(4): 1290 - 1296. [Abstract] [Full Text] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||