Using time‐domain reflectometry to characterize shallow solute transport in an oak woodland hillslope in northern California,USA

The natural heterogeneity of water and solute movement in hillslope soils makes it difficult to accurately characterize the transport of surface‐applied pollutants without first gathering spatially distributed hydrological data. This study examined the application of time‐domain reflectometry (TDR) to measure solute transport in hillslopes. Three different plot designs were used to examine the transport of a conservative tracer in the first 50 cm of a moderately sloping soil. In the first plot, which was designed to examine spatial variability in vertical transport in a 1·2 m² plot, a single probe per meter was found to adequately characterize vertical solute travel times. In addition, a dye and excavation study in this plot revealed lateral preferential flow in small macropores and a transport pattern where solute is focused vertically into preferential flow pathways. The bypass flow delivers solute deeper in the soil, where lateral flow occurs. The second plot, designed to capture both vertical and lateral flow, provided additional evidence confirming the flow patterns identified in the excavation of the first plot. The third plot was designed to examine lateral flow and once again preferential flow of the tracer was observed. In one instance rapid solute transport in this plot was estimated to occur in as little as 3% of the available pore space. Finally, it was demonstrated that the soil anisotropy, although partially responsible for lateral subsurface transport, may also homogenize the transport response across the hillslope by decreasing vertical solute spreading. Copyright © 2002 John Wiley & Sons, Ltd.     

Characterizing heterogeneous soil water flow and solute transport using information measures

Heterogeneous water flow and solute transport in soils are an important phenomenon and difficult to be characterized. The objectives of this study were to investigate the heterogeneity of solute transport related to heterogeneous soil water flow using dye infiltration experiments, and to characterize heterogeneous water flow and solute transport in soils using the information theory. Field experiments of dye infiltration were performed in four plots. Various information measures were applied to characterize information content and complexity of water flow and solute transport in soils. Information contents and complexities of the maximum and apparent infiltration depths, and the mean and standard deviation of concentrations in the vertical direction of the plots were calculated. More heterogeneous processes of soil water flow and transport result in higher information/complexity values. The probability distributions of mean concentration were similar to those of the corresponding apparent infiltration depths for the plots, indicating that heterogeneity of dye concentrations was closely related to that of soil water flow. However, the range of information entropy and complexity of the water flow sequences was much narrower than that of the sequences of the concentrations. The results suggested that the transport processes were more heterogeneous than the water flow processes. Compared with the probability distributions of flow parameters, the information measures appeared to be a more versatile tool to describe flow and transport heterogeneities in soils.     

Temporal variability of solute transport under vadose zone conditions

The heterogeneity of the solute flux field in the horizontal plane at the field scale has been documented in several field studies. On the other hand, little information is available on the persistence of certain solute transport scenarios over consecutive infiltration cycles. This study was initiated to analyse the recurrence of solute leaching behaviour as estimated in two soil column tests emphasizing the preferential flow phenomenon. Twenty-four small-sized soil samples were subjected to two consecutive unsaturated steady-state flow leaching experiments with bromide as tracer. Observed breakthrough curves (BTCs) were analysed by the method of moments and by the advection–dispersion equation (ADE) to classify solute behaviour. Frequency distributions of the parameters indicating the solute velocity were heavily skewed or bimodal, reflecting the broad variability of the leaching scenarios, including some with pronounced preferential solute breakthrough. Exclusion of the preferential flow columns from our calculations revealed an average amount of 37% of immobile water. The large-scale BTCs derived from assembling the individual concentration courses of each run showed similar features, such as an early bromide breakthrough. However, two distinct apices, viz. one preferential and one matrix, were observed only in the first run, whereas the concentration decrease between the peaks was missing from the second run. A change in soil structure with continuous leaching was presumed to modify the interplay of the various flow domains, thereby altering the spreading of the BTCs. Correlation analysis between parameters of both tests suggests that preferential transport conditions are likely to occur at the same locations in the field over several infiltration cycles, whereas the ‘classical’ or expected matrix flow is time variant and therefore seems to be hardly predictable. © 1998 John Wiley & Sons, Ltd.     

Comparison of soil salinity and solute transport for different cultivated soil types in northeastern Egypt / Comparaison de la salinité du sol et du transport de solutés pour différents types de sols cultivés du nord-est de l'Egypte

Abstract

The soil salinity distribution and solute transport properties of three different soil types were investigated and compared within a project area in northeastern Egypt. For this purpose, dye tracer experiments and salinity sampling were carried out. The resulting salinity maps showed that the soil salinity in the cultivated western site of the project area is 8–10 times higher than that in the cultivated eastern site. However, the cultivated soil displayed significantly lower salinity with higher uniformity as compared to the uncultivated soil. The preferential flow phenomenon was less apparent in the cultivated soil. This is mainly due to tillage which disrupts the structure of the soil so that deep cracks are no longer connected to the soil surface. This reduces the risk for groundwater contamination through preferential flow. The study showed that careful and continuous monitoring of the salinity status is needed now and in the future.     

Effects of antecedent moisture and macroporosity on infiltration and water flow in frozen soil

Infiltration into frozen soil plays an important role in soil freeze–thaw and snowmelt-driven hydrological processes. To better understand the complex thermal energy and water transport mechanisms involved, the influence of antecedent moisture content and macroporosity on infiltration into frozen soil was investigated. Ponded infiltration experiments on frozen macroporous and non-macroporous soil columns revealed that dry macroporous soil produced infiltration rates reaching 10³ to 10⁴ mm day⁻¹, two to three orders of magnitude larger than dry non-macroporous soil. Results suggest that rapid infiltration and drainage were a result of preferential flow through initially air-filled macropores. Using recorded flow rates and measured macropore characteristics, calculations indicated that a combination of both saturated flow and unsaturated film flow likely occurred within macropores. Under wet conditions, regardless of the presence of macropores, infiltration was restricted by the slow thawing rate of pore ice, producing infiltration rates of 2.8 to 5.0 mm day⁻¹. Reduced preferential flow under wet conditions was attributed to a combination of soil swelling, due to smectite-rich clay (that reduced macropore volume), and pore ice blockage within macropores. In comparison, dry soil column experiments demonstrated that macropores provided conduits for water and thermal energy to bypass the frozen matrix during infiltration, reducing thaw rates compared with non-macroporous soils. Overall, results showed the dominant control of antecedent moisture content on the initiation, timing, and magnitude of infiltration and flow in frozen macroporous soils, as well as the important role of macropore connectivity. The study provides an important data set that can aid the development of hydrological models that consider the interacting effects of soil freeze–thaw and preferential flow on snowmelt partitioning in cold regions.     

Differences in preferential flow with antecedent moisture conditions and soil texture: Implications for subsurface P transport

Preferential flowpaths transport phosphorus (P) to agricultural tile drains. However, if and to what extent this may vary with soil texture, moisture conditions, and P placement is poorly understood. This study investigated (a) interactions between soil texture, antecedent moisture conditions, and the relative contributions of matrix and preferential flow and (b) associated P distributions through the soil profile when fertilizers were applied to the surface or subsurface. Brilliant blue dye was used to stain subsurface flowpaths in clay and silt loam plots during simulated rainfall events under wet and dry conditions. Fertilizer P was applied to the surface or via subsurface placement to plots of different soil texture and moisture condition. Photographs of dye stains were analysed to classify the flow patterns as matrix dominated or macropore dominated, and soils within plots were analysed for their water‐extractable P (WEP) content. Preferential flow occurred under all soil texture and moisture conditions. Dye penetrated deeper into clay soils via macropores and had lower interaction with the soil matrix, compared with silt loam soil. Moisture conditions influenced preferential flowpaths in clay, with dry clay having deeper infiltration (92 ± 7.6 cm) and less dye–matrix interaction than wet clay (77 ± 4.7 cm). Depth of staining did not differ between wet (56 ± 7.2 cm) and dry (50 ± 6.6 cm) silt loam, nor did dominant flowpaths. WEP distribution in the top 10 cm of the soil profile differed with fertilizer placement, but no differences in soil WEP were observed at depth. These results demonstrate that large rainfall events following drought conditions in clay soil may be prone to rapid P transport to tile drains due to increased preferential flow, whereas flow in silt loams is less affected by antecedent moisture. Subsurface placement of fertilizer may minimize the risk of subsurface P transport, particularily in clay.     

ABSENCE OF PREFERENTIAL FLOW IN THE PERCOLATING WATERS OF A CONIFEROUS FOREST SOIL

Evidence for the functioning of macropores and the presence of preferential flow in forest soils is equivocal. This is partly because many workers use only one diagnostic technique to indicate whether or not macropore flow occurs. In this paper three lines of evidence are used to suggest that preferential flow does not occur in the percolating waters of a coniferous forest soil under the range of hydrological conditions that prevail in the field. To simulate field conditions, realistic rainfall intensities were used in conservative solute transport experiments on four undisturbed soil columns. A method is described in which breakthrough data can be used to calculate the percentage of antecedent water displaced from a soil column during frontal-type breakthrough experiments. Calculations based on this method using the experimental data show that as little as five percent of the antecedent water was immobile. The simple form of the functional advection–dispersion equation, based on a single value for linear velocity and the dispersion coefficient was fitted to two of the breakthrough curves with reasonable accuracy, further suggesting that preferential flow did not occur in the experiments. Finally, soil moisture characteristic curves were determined for replicate soil samples from the forest soil. The operational water contents of the columns during the breakthrough experiments were compared with the soil moisture characteristics and it was found that pores exerting pressure heads greater than −0·5 kPa did not appear to contribute to flow through the columns, again suggesting an absence of preferential flow. © 1997 by John Wiley & Sons, Ltd.     

Simulation of field‐scale water flow and bromide transport in a cracked clay soil

Single domain models may seriously underestimate leaching of nutrients and pesticides to groundwater in clay soils with shrinkage cracks. Various two‐domain models have been developed, either empirical or physically based, which take into account the effects of cracks on water flow and solute transport. This paper presents a model concept that uses the clay shrinkage characteristics to derive crack volume and crack depth under transient field conditions. The concept has been developed to simulate field average behaviour of a field with cracks, rather than flow and transport at a small plot. Water flow and solute transport are described with basic physics, which allow process and scenario analysis. The model concept is part of the more general agrohydrological model SWAP, and is applied to a field experiment on a cracked clay soil, at which water flow and bromide transport were measured during 572 days. A single domain model was not able to mimic the field‐average water flow and solute transport. Incorporation of the crack concept considerably improved the simulation of water content and bromide leaching to the groundwater. Still deviations existed between the measured and simulated bromide concentration profiles. The model did not reproduce the observed bromide retardation in the top layer and the high bromide dispersion resulting from water infiltration at various soil depths. A sensitivity analysis showed that the amounts of bromide leached were especially sensitive to the saturated hydraulic conductivity of the top layer, the solute transfer from the soil matrix to crack water flow and the mean residence time of rapid drainage. The shrinkage characteristic and the soil hydraulic properties of the clay matrix showed a low sensitivity. Copyright © 2000 John Wiley & Sons, Ltd.     

Tracing solute infiltration using a combined method of dye tracer test and electrical resistivity tomography in an undisturbed forest soil profile

An accurate prediction of solute infiltration in a soil profile is important in the area of environmental science, groundwater and civil engineering. We examined the infiltration pattern and monitored the infiltration process using a combined method of dye tracer test and electrical resistivity tomography (ERT) in an undisturbed field soil (1 m × 1 m). A homogeneous matrix flow was observed in the surface soil (A horizon), but a preferential flow along macropores and residual rock structure was the dominant infiltration pattern in the subsurface soil. Saturated interflow along the slopping boundaries of A and C1 horizons and of an upper sandy layer and a lower thin clay layer in the C horizon was also observed. The result of ERT showed that matrix flow started first in A horizon and then the infiltration was followed by the preferential flows along the sloping interfaces and macropores. The ERT did not show as much detail as the dye‐stained image for the preferential flow. However, the area with the higher staining density where preferential flow was dominant showed a relatively lower electrical resistivity. The result of this study indicates that ERT can be applied for the monitoring of solute transportation in the vadose zone. Copyright © 2010 John Wiley & Sons, Ltd.     

Effect of anisotropy on solute transport in degraded fen peat soils

Peat soils are heterogeneous, anisotropic porous media. Compared to mineral soils, there is still limited understanding of physical and solute transport properties of fen peat soils. In this study, we aimed to explore the effect of soil anisotropy on solute transport in degraded fen peat. Undisturbed soil cores, taken in vertical and horizontal direction, were collected from one drained and one restored fen peatland both in a comparable state of soil degradation. Saturated hydraulic conductivity (K _s) and chemical properties of peat were determined for all soil cores. Miscible displacement experiments were conducted under saturated steady state conditions using potassium bromide as a conservative tracer. The results showed that (1) the K _s in vertical direction (K _sv) was significantly higher than that in horizontal direction (K_sh), indicating that K _s of degraded fen peat behaves anisotropically; (2) pronounced preferential flow occurred in vertical direction with a higher immobile water fraction and a higher pore water velocity; (3) the 5% arrival time (a proxy for the strength of preferential flow) was affected by soil anisotropy as well as study site. A strong correlation was found between 5% arrival time and dispersivity, K _s and mobile water fraction; (4) phosphate release was observed from drained peat only. The impact of soil heterogeneity on phosphate leaching was more pronounced than soil anisotropy. The soil core with the strongest preferential flow released the highest amount of phosphate. We conclude that soil anisotropy is crucial in peatland hydrology but additional research is required to fully understand anisotropy effects on solute transport.     

Microtopography,water storage and flow patterns in a fine‐textured soil under agricultural use

Agricultural use of soils implies tillage and often compaction and therefore influences processes on soil surface and affects infiltration of water into the subsoil. Although many studies on soil surface processes or flow patterns in soils exist, works relating both are rare in literature. We did two tracer experiments with Brilliant Blue FCF on a tilled and compacted plot and a non‐tilled one to investigate water storage on the soil surface during simulated rainfall and changes of soil microtopography, to analyse the associated flow patterns in the soil and to relate both to tillage and compaction. Our results show that storage was larger on the tilled and compacted plot than on the non‐tilled one. After tillage, transport processes above the plough pan were partly disconnected from those underneath because macropores were disrupted and buried by the tillage operation. However, preferential flow along cracks occurred on both plots and the macropores buried below the tillage pan still functioned as preferential flow paths. Therefore, we conclude that the studied soil is susceptible to deep vertical solute propagation at dry conditions when cracks are open, irrespective of tillage and compaction. Copyright © 2012 John Wiley & Sons, Ltd.     

Movement of Heavy Metals in Soil through Preferential Flow Paths under Different Rainfall Intensities

Preferential flow in soils deserves attention due to its potential role in accelerating the movement of contaminants to groundwater. This study investigates the movement of Cd, Cu, and Pb through preferential flow paths under different applied rainfall intensities. Artificial acid rain (pH of 4.1) containing CdCl₂, CuCl₂, and Pb(NO₃)₂ was applied to undisturbed soil and repacked sand columns at low and high intensities, and leachate metals and chloride concentrations were measured. Cd was found in the leachate at both low and high rates in all columns, while the increase in Cu concentrations in the leachate was detected only at the high rate of the undisturbed columns. Pb was retained in both columns. For undisturbed columns, the breakthrough curves of Cd and Cu were similar to those of Cl, showing early initial breakthrough by preferential flow and dependency on rainfall intensities. The Cd concentrations were detected in the leachate from repacked columns for high rate rainfall, implying that even homogeneous soil may not be perfectly able to retain metals and the initially wet condition is more harmful for subsurface contaminant transport. In conclusion, the study demonstrated that, despite its highly sorptive nature, the transport of some metals may be as fast as that of a tracer under preferential flow conditions, and the rainfall intensity is a significant factor for the degree of transport.     

Double funnelling in a mature coastal British Columbia forest: spatial patterns of stemflow after infiltration

Many studies have focused on the amount of stemflow in different forests and for different rainfall events, but few studies have focused on how stemflow intensity varies during events or the infiltration of stemflow into the soil. Stemflow may lead to higher water delivery rates at the base of the tree compared with throughfall over the same area and fast and deeper infiltration of this water along roots and other preferential flow pathways. In this study, stemflow amounts and intensities were measured and blue dye experiments were conducted in a mature coniferous forest in coastal British Columbia to examine double funnelling of stemflow. Stemflow accounted for only 1% of precipitation and increased linearly with event total precipitation. Funnelling ratios ranged from less than 1 to almost 20; smaller trees had larger funnelling ratios. Stemflow intensity generally was highest for periods with high‐intensity rainfall later in the event. The maximum stemflow intensities were higher than the maximum precipitation intensities. Dye tracer experiments showed that stemflow infiltrated primarily along roots and was found more frequently at depth than near the soil surface. Lateral flow of stemflow was observed above a dense clay layer for both the throughfall and stemflow experiments. Stemflow appeared to infiltrate deeper (122 cm) than throughfall (85 cm), but this difference was in part a result of site‐specific differences in maximum soil depth. However, the observed high stemflow intensities combined with preferential flow of stemflow may lead to enhanced subsurface stormflow. This suggests that even though stemflow is only a very minor component of the water balance, it may still significantly affect soil moisture, recharge, and runoff generation. Copyright © 2016 John Wiley & Sons, Ltd.     

The effect of alternating different water qualities on accumulation and leaching of solutes in a Mediterranean cracking soil

The relevance of bypass flow on water flow, solute or pesticide transport is becoming increasingly recognized. Recent investigations proved that soil salinization may be influenced by bypass flow, i.e. the rapid transport of water and solutes via macropores and/or shrinkage cracks to subsoil and groundwater. This paper explores the role of bypass flow in the process of accumulation and leaching of solutes, as well as of sodium, in a Mediterranean cracking soil irrigated with saline/sodic waters. The results of bypass flow experiments performed on undisturbed soil cores showed that leaching of solutes occurred in concomitance with bypass fluxes when a low salinity solution was alternated with a high salinity solution. Exchange of solutes between the incoming solution and the soil matrix occurred during the bypass flow events at the contact surfaces (cracks walls) between the solution and the soil matrix and where cracks terminated in the soil samples. Concomitant exchanges of sodium were indicated by measurements performed in the effluent solution during the bypass flow measurements. The amount of Sodium released from the soil during the bypass flow events, as well as that of the soluble salts leached from the soil, were found to depend on the degree of soil cracking. These results indicated that:

• 1 in management of irrigation in cracking soils, under the occurrence of bypass fluxes, alternating a low salinity/sodicity water with a high salinity/sodicity solution can be effective for preventing salinization and sodification:
• 2 greater efficiency of removal of sodium/soluble salts can be obtained if application of the leaching solution is performed when the soil is at a considerable degree of cracking.

Copyright © 2002 John Wiley & Sons, Ltd.     

A 3D ERT study of solute transport in a large experimental tank

A high resolution, cross-borehole, 3D electrical resistivity tomography (ERT) study of solute transport was conducted in a large experimental tank. ERT voxels comprising the time sequence of electrical images were converted into a 3D array of ERT estimated fluid conductivity breakthrough curves and compared with direct measurements of fluid conductivity breakthrough made in wells. The 3D ERT images of solute transport behaviour were also compared with predictions based on a 3D finite-element, coupled flow and transport model, accounting for gravity induced flow caused by concentration differences.The tank (dimensions 185×245×186 cm) was filled with medium sand, with a gravel channel and a fine sand layer installed. This heterogeneous system was designed to complicate solute transport behaviour relative to a homogeneous sand tank, and to thus provide a challenging but insightful analysis of the ability of 3D ERT to resolve transport phenomena. Four ERT arrays and 20 piezometers were installed during filling. A NaCl tracer (conductivity 1.34 S/m) was injected and intensively monitored with 3D ERT and direct sampling of fluid chemistry in piezometers.We converted the bulk conductivity estimate for 250 voxels in the ERT imaged volume into ERT estimated voxel fluid conductivity by assuming that matrix conduction in the tank is negligible. In general, the ERT voxel response is in reasonable agreement with the shape of fluid conductivity breakthrough observed in six wells in which direct measurements of fluid conductivity were made. However, discrepancies occur, particularly at early times, which we attribute to differences between the scale of the image voxels and the fluid conductivity measurement, measurement errors mapped into the electrical inversion and artificial image roughness resulting from the inversion.ERT images revealed the 3D tracer distribution at 15 times after tracer injection. The general pattern and timing of solute breakthrough observed with ERT agreed with that predicted from the flow/transport modelling. However, the ERT images indicate a vertical component of tracer transport and preferential flow paths in the medium sand. We attribute this to transient vertical gradients established during tracer injection, and heterogeneity caused by sorting of the sand resulting from the filling procedure. In this study, ERT provided a unique dataset of 250 voxel breakthrough curves in 1.04 m³. The use of 3D ERT to generate an array of densely sampled estimated fluid conductivity breakthrough curves is a potentially powerful tool for quantifying solute transport processes.     

Characterising deep vadose zone water movement and solute transport under typical irrigated cropland in the North China Plain

Understanding the dynamics and mechanisms of soil water movement and solute transport is essential for accurately estimating recharge rates and evaluating the impacts of agricultural activities on groundwater resources. In a thick vadose zone (0–15 m) under irrigated cropland in the piedmont region of the North China Plain, soil water content, matric potential, and solute concentrations were measured. Based on these data, the dynamics of soil water and solutes were analysed to investigate the mechanisms of soil water and solute transport. The study showed that the 0–15‐m vadose zone can be divided into three layers: an infiltration and evaporation layer (0–2 m), an unsteady infiltration layer (2–6 m), and a quasi‐steady infiltration layer (6–15 m). The chloride, nitrate, and sulphate concentrations all showed greater variations in the upper soil layer (0–1 m) compared to values in the deep vadose zone (below 2 m). The average concentrations of these three anions in the deep vadose zone varied insignificantly with depth and approached values of 125, 242, and 116 mg/L. The accumulated chloride, sulphate, and nitrate were 2,179 ± 113, 1,760 ± 383, and 4,074 ± 421 kg/ha, respectively. The soil water potential and solute concentrations indicated that uniform flow and preferential flow both occurred in the deep vadose zone, and uniform flow was the dominant mechanism of soil water movement in this study. The piston‐like flow velocity of solute transport was 1.14 m per year, and the average value of calculated leached nitrate nitrogen was 107 kg/ha?year below the root zone. The results can be used to better understand recharge processes and improve groundwater resources management.     

Leaching of Salmonella enterica in clay columns comparing two manure application methods

Transfer of zoonotic bacterial pathogens through intact soil columns was monitored in an outdoor lysimeter over 36 d. Manure spiked with Salmonella enterica serovar Senftenberg was applied to either the soil surface or injected 0.08 m into the soil to compare leaching associated with the two manure application methods. The highest concentrations of S. enterica (up to 60,000 S. enterica CFU/mL) were detected on Day 1 in the first drainage samples, with measurable but declining concentrations persisting for 10 to 36 d depending on replicate columns. The total recovery of leached S. enterica in drainage samples ranged from 0.08% to 13.8%. When comparing the two application methods, there was no statistically significant difference in the leaching concentration of S. enterica at each sampling time during the study period. In addition, comparison of enumerations by selective plating and real-time polymerase chain reaction yielded similar concentrations of S. enterica, indicating that mainly viable and culturable cells were leached from the columns. When the experiment was terminated, the fluorescent dye Acid Yellow was applied to four selected columns and the distribution of dye and size of active (dye-stained) pores were measured with a digital camera and visualization software. The profiles showed that the area covered by active pores ranged from 0.1% to 3.6%. The relatively small fraction of active pores in the soil profile was consistent with the evidence of rapid transport of S. enterica and chloride in the columns.     

Tracer and hydrometric study of preferential flow in large undisturbed soil cores from the Georgia Piedmont,USA

We studied the temporal patterns of tracer throughput in the outflow of large (30 cm diameter by 38 cm long) undisturbed cores from the Panola Mountain Research Watershed, Georgia. Tracer breakthrough was affected by soil structure and rainfall intensity. Two rainfall intensities (20 and 40 mm hr⁻¹) for separate Cl⁻ and Br⁻ amended solutions were applied to two cores (one extracted from a hillslope soil and one extracted from a residual clay soil on the ridge). For both low and high rainfall intensity experiments, preferential flow occurred in the clay core, but not in the hillslope core. The preferential flow is attributed to well‐developed interpedal macrochannels that are commonly found in structured clay soils, characteristic of the ridge site. However, each rainfall intensity exceeded the matrix infiltration capacity at the top of the hillslope core, but did not exceed the matrix infiltration capacity at the middle and bottom of the hillslope core and at all levels in the clay core. Localized zones of saturation created when rainfall intensity exceeds the matrix infiltration capacity may cause water and tracer to overflow from the matrix into macrochannels, where preferential flow occurs to depth in otherwise unsaturated soil. Copyright © 1999 John Wiley & Sons, Ltd.