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.     

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.     

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

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

Sediment and solute transport on soil slope under simultaneous influence of rainfall impact and scouring flow

Soil erosion and nutrient losses with surface runoff in the loess plateau in China cause severe soil quality degradation and water pollution. It is driven by both rainfall impact and runoff flow that usually take place simultaneously during a rainfall event. However, the interactive effect of these two processes on soil erosion has received limited attention. The objectives of this study were to better understand the mechanism of soil erosion, solute transport in runoff, and hydraulic characteristics of flow under the simultaneous influence of rainfall and shallow clear‐water flow scouring. Laboratory flume experiments with three rainfall intensities (0, 60, and 120 mm h⁻¹) and four scouring inflow rates (10, 20, 30, and 40 l min⁻¹) were conducted to evaluate their interactive effect on runoff. Results indicate that both rainfall intensity and scouring inflow rate play important roles on runoff formation, soil erosion, and solute transport in the surface runoff. A rainfall splash and water scouring interactive effect on the transport of sediment and solute in runoff were observed at the rainfall intensity of 60 mm h⁻¹ and scouring inflow rates of 20 l min⁻¹. Cumulative sediment mass loss (Ms) was found to be a linear function of cumulative runoff volume (Wr) for each treatment. Solute transport was also affected by both rainfall intensity and scouring inflow rate, and the decrease in bromide concentration in the runoff with time fitted to a power function well. Reynolds number (Re) was a key hydraulic parameter to determine erodability on loess slopes. The Darcy–Weisbach friction coefficients (f) decreased with the Reynolds numbers (Re), and the average soil and water loss rate (M_l) increased with the Reynolds numbers (Re) on loess slope for both scenarios with or without rainfall impact. Copyright © 2010 John Wiley & Sons, Ltd.     

The effect of non divergence-free velocity fields on field scale ground water solute transport

Solute plume subjected to field scale hydraulic conductivity heterogeneity shows a large dispersion/macrodispersion, which is the manifestation of existing fields scale heterogeneity on the solute plume. On the other hand, due to the scarcity of hydraulic conductivity measurements at field scale, hydraulic conductivity heterogeneity can only be defined statistically, which makes the hydraulic conductivity a random variable/function. Random hydraulic conductivity as a parameter in flow equation makes the pore flow velocity also random and the ground water solute transport equation is a stochastic differential equation now. In this study, the ensemble average of stochastic ground water solute transport equation is taken by the cumulant expansion method in order to upscale the laboratory scale transport equation to field scale by assuming pore flow velocity is a non stationary, non divergence-free and unsteady random function of space and time. Besides the stochastic explanation of macrodispersion and the velocity correction term obtained by Kavvas and Karakas (J Hydrol 179:321–351, 1996) before a new velocity correction term, which is a function of mean pore flow velocity divergence, is obtained in this study due to strict second order cumulant expansion (without omitting any term after the expansion) performed. The significance of the new velocity correction term is investigated on a one dimensional transport problem driven by a density dependent flow field.     

Calculating the macrodispersion coefficient of the ensemble averaged solute transport equation in the discrete domain

Recognizing the spatial heterogeneity of hydraulic parameters, many researchers have studied the solute transport by both groundwater and channel flow in a stochastic framework. One of the methodologies used to up‐scale the stochastic solute transport equation, from a point‐location scale to a grid scale, is the cumulant expansion method combined with the calculus for the time‐ordered exponential and the calculus for the Lie operator. When the point‐location scale transport equation is scaled up to the grid scale, using the cumulant expansion method, a new dispersion coefficient emerges in the dispersive term of the solute transport equation in addition to the molecular dispersion coefficient. This velocity driven dispersion is called ‘macrodispersion’. The macrodispersion coefficient is the integral function of the time‐ordered covariance of the random velocity field. The integral is calculated over a Lagrangian trajectory of the flow. The Lagrangian trajectory depends on the following: (i) the spatial origin of the particle; (ii) the time when the macrodispersion is calculated; and (iii) the mean velocity field along the trajectory itself. The Lagrangian trajectory is a recursive function of time because the location of the particle along the trajectory at a particular time depends on the location of the particle at the previous time. This recursive functional form of the Lagrangian trajectory makes the calculation of the macrodispersion coefficient difficult. Especially for the unsteady, spatially non‐stationary, non‐uniform flow field, the macrodispersion coefficient is a highly complex expression and, so far, calculated using numerical methods in the discrete domains. Here, an analytical method was introduced to calculate the macrodispersion coefficient in the discrete domain for the unsteady and steady, spatially non‐stationary flow cases accurately and efficiently. This study can fill the gap between the theory of the ensemble averaged solute transport model and its numerical implementations. Copyright © 2011 John Wiley & Sons, Ltd.     

Experimental and analytical studies of solute transport during run‐off over vegetated surfaces

Solute transport in overland flow is considered as one of the main contributors to water pollution. Although many models of pollutant transport mechanism from soil to run‐off water have been proposed, the characteristics of solute transport accompanying the water run‐off over vegetated surface have not been well studied. In this study, a series of laboratory experiments were conducted to study the solute transport over vegetated surfaces. Based on the experimental results, an idea of the “stationary water layer” in run‐off was proposed. Applying the complete mixing theory in the stationary water layer, an analytical solute transport model was developed with the assumption that the upper run‐off completely mixes with the underlying water in the stationary water layer for each site. The results show that the predictions made by the present model are in good agreement with the measured experimental data. For the vegetated surfaces, the depth of stationary water layer is related to the rainfall intensity, bed slope, and vegetation density. The analytical solution shows that the maximum solute transport occurs at the time of concentration. This study advances our understanding of the mechanisms of solute transport over vegetated areas.     

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.     

Contaminant transport in a small deformation aquitard affected by the delayed drainage phenomenon

This paper presents a formulation accounting for the effect of delayed drainage phenomenon (DDP) on the breakthrough of contaminant flux in an aquitard, by considering the movement of soil particles, porosity variation, hydraulic head variation, and transient flow during the consolidation. The water flow equation in an aquitard was based on the Terzaghi's consolidation theory, and the contaminant transport equation was derived on the basis of the mass balance law. Two cases were used to illustrate the effect of DDP on the contaminant transport in an aquitard of small deformation. It is found that the breakthrough time of contaminant in an aquitard is very long, which is mainly ascribed to the low permeability of aquitard and sorption of soil particles. It is also found that the increase of depletion, which is in general induced by the increase of thickness and specific storativity and the decrease of hydraulic conductivity, enhances the impact of DDP on the contaminant transport in an aquitard. A larger delay index (τ₀) of DDP gives a greater delay breakthrough time (DBT) of solute transport in an aquitard, which controls the difference of the breakthrough time of contaminant transport in aquitards with and without the occurrence of DDP. For the cases where advection plays a dominant role during the process of solute transport, τ₀ is almost linearly correlated with DBT, and the ratio of DBT over the breakthrough time without consideration of DDP also approximately shows a linear relationship with the ratio of specific storativity to porosity, given a fixed drawdown in the adjacent aquifer with the sorption being ignored.     

Impact of debris dams on hyporheic interaction along a semi‐arid stream

Hyporheic exchange increases the potential for solute retention in streams by slowing downstream transport and increasing solute contact with the substrate. Hyporheic exchange may be a major mechanism to remove nutrients in semi‐arid watersheds, where livestock have damaged stream riparian zones and contributed nutrients to stream channels. Debris dams, such as beaver dams and anthropogenic log dams, may increase hyporheic interactions by slowing stream water velocity, increasing flow complexity and diverting water to the subsurface. Here, we report the results of chloride tracer injection experiments done to evaluate hyporheic interaction along a 320 m reach of Red Canyon Creek, a second order stream in the semi‐arid Wind River Range of Wyoming. The study site is part of a rangeland watershed managed by The Nature Conservancy of Wyoming, and used as a hydrologic field site by the University of Missouri Branson Geologic Field Station. The creek reach we investigated has debris dams and tight meanders that hypothetically should enhance hyporheic interaction. Breakthrough curves of chloride measured during the field experiment were modelled with OTIS‐P, a one‐dimensional, surface‐water, solute‐transport model from which we extracted the storage exchange rate α and cross‐sectional area of the storage zone A_s for hyporheic exchange. Along gaining reaches of the stream reach, short‐term hyporheic interactions associated with debris dams were comparable to those associated with severe meanders. In contrast, along the non‐gaining reach, stream water was diverted to the subsurface by debris dams and captured by large‐scale near‐stream flow paths. Overall, hyporheic exchange rates along Red Canyon Creek during snowmelt recession equal or exceed exchange rates observed during baseflow at other streams. Copyright © 2005 John Wiley & Sons, Ltd.     

Bedload transport processes in a braided gravel-bed river model

A 1:50 scale hydraulic model was designed, based on Froude number similarity and using hydrological and sediment data from a small braided gravel-bed river (the North Branch of the Ashburton River, Canterbury, New Zealand). Eighteen experiments were conducted; seven using steady flows, and eleven using unsteady flows. The experiments were carried out in a 20 m × 3 m tilting flume equipped with a continuous sediment feed and an automated data acquisition and control system. In all experiments water at 30°C was used to reduce viscosity-related scale effects. Analyses of the experimental data revealed that bedload transport rates in braided channels are highly variable, with relative variability being inversely related to mean bedload transport rate. Variability was also found to be cyclic with short-term variations being caused by the migration of bedforms. Bedload transport was found to be more efficient under steady flow than under unsteady flow, and it was postulated that this is caused by a tendency for channel form to evolve towards a condition which maximizes bedload transport for the occurring flow. Average bedload transport rate was found to vary with channel form, although insufficient measurements were made to define a relationship.     

Modelling the process of soil detachment by rill flow on steep loessial hillslopes

Soil detachment by rill flow is a key process of rill erosion, modelling this process can help in understanding rill erosion mechanisms. However, many soil detachment models are established on conceptual assumptions rather than experimental data. The objectives of this study were to establish a model of soil detachment by rill flow based on flume experimental data and to quantitatively verify the model. We simulated the process of soil detachment by rill flow in flume experiments with a soil-feeding hopper using loessial soil on steep slopes. Seven flow discharges, six slopes and five sediment loads were combined. Soil detachment capacity, sediment transport capacity, and soil detachment rate by rill flow under different sediment loads were measured. The process of soil detachment by rill flow can be modelled by a dual power function based on soil detachment capacity and transport capacity deficit as variables. The established model exhibited high credibility (NSE=0.97; R²=0.97). The contributions of soil detachment capacity and transport capacity deficit to soil detachment rate by rill flow reached 60% and 36%, respectively. Soil detachment capacity exerted more influence on soil detachment rate than did transport capacity deficit. The performance of the WEPP rill erosion equation is also favourable (NSE=0.95; R²=0.97). The two power exponents in the model we established strengthen the role of soil detachment capacity in soil detachment rate and weaken that for transport capacity deficit. Soil detachment capacity and transport capacity deficit played important roles in the determination of soil detachment rate by rill flow. The results can be applied to implement the numerical modeling and prediction of rill erosion processes on steep loessial hillslopes. © 2019 John Wiley & Sons, Ltd.     

FLOW SEPARATION IN UNDISTURBED SOILS USING MULTIPLE ANIONIC TRACERS. PART 1. ANALYTICAL METHODS AND UNSTEADY RAINFALL AND RETURN FLOW EXPERIMENTS

Tracers provide one of the few ways of obtaining realistic information on the flows of water and solutes in undisturbed structured soils. Three fluorinated organic acids [pentafluorobenzoic acid, 2,6-diflurobenzoic acid and o-(trifluoromethyl) benzoic acid] and bromide were tested as anionic tracers in situations where the separation of different flow components is of interest. The fluorobenzoates were relatively conservative (approximately 90% for loam soil) although, in some instances, co-elution or complexation may lead to apparently non-conservative behaviour. Tracer mixtures applied to soils, were separated and quantified by high-performance liquid chromatography. The relative mobilities of the tracers were studied during unidirectional steady flow in large undisturbed soil columns brought back to the laboratory. The breakthrough curves showed highly dispersive behaviour, with very early breakthrough of solute and a long tail. Contributions to the column outflow of different simultaneous applications of rainfall and return flow, labelled with different tracers, were separated using multicomponent mixing equations. The results show the importance of preferential flow and relatively immobile storage in the transport process.     

Depth‐averaged,spatially distributed flow dynamic and solute transport modelling of a large‐scaled,subtropical constructed wetland

Constructed wetlands are being utilized worldwide to effectively reduce excess nutrients in agricultural runoff and wastewater. Despite their frequency, a multi‐dimensional, physically based, spatially distributed modelling approach has rarely been applied for flow and solute transport in treatment wetlands. This article presents a two‐dimensional hydrodynamic and solute transport modelling of a large‐scaled, subtropical, free water surface constructed wetland of about 8 km² in the Everglades of Florida, USA. In this study, MIKE 21 was adopted as the basic model framework. Field monitoring of the time series hydrological and chloride data, as well as spatially distributed data such as bathymetry and vegetation distribution, provided the necessary model input and testing data. Simulated water level profiles were in good agreement with the spatio‐temporal variations of measured ones. On average, the root‐mean‐square error of model calibration on annual water level fluctuations was 0·09 m. Manning's roughness coefficients for the dense emergent and submerged aquatic vegetation areas, which were estimated as a function of vegetation type, ranged from 0·67 to 1·0 and 0·12 to 0·15 s/m^1/3, respectively. The solute transport model calibration for four monitoring sites agreed well with the measured annual variations in chloride concentration with an average percent model error of about 15%. The longitudinal dispersivity was estimated to be about 2 m and was more than an order of magnitude higher than the transverse one. This study is expected to play the role of a stepping stone for future modelling efforts on the development and application of more advanced flow and transport models applicable to a variety of constructed wetland systems, as well as to the Everglades stormwater treatment areas in operation or in preparation. Copyright © 2010 John Wiley & Sons, Ltd.     

Watershed structural influences on the distributions of stream network water and solute travel times under baseflow conditions

Watershed structure influences the timing, magnitude, and spatial location of water and solute entry to stream networks. In turn, stream reach transport velocities and stream network geometry (travel distances) further influence the timing of export from watersheds. Here, we examine how watershed and stream network organization can affect travel times of water from delivery to the stream network to arrival at the watershed outlet. We analysed watershed structure and network geometry and quantified the relationship between stream discharge and solute velocity across six study watersheds (11.4 to 62.8 km²) located in the Sawtooth Mountains of central Idaho, USA. Based on these analyses, we developed stream network travel time functions for each watershed. We found that watershed structure, stream network geometry, and the variable magnitude of inputs across the network can have a pronounced affect on water travel distances and velocities within a stream network. Accordingly, a sample taken at the watershed outlet is composed of water and solutes sourced from across the watershed that experienced a range of travel times in the stream network. We suggest that understanding and quantifying stream network travel time distributions are valuable for deconvolving signals observed at watershed outlets into their spatial and temporal sources, and separating terrestrial and in‐channel hydrological, biogeochemical, and ecological influences on in‐stream observations. Copyright © 2016 John Wiley & Sons, Ltd.     

The role of eddies in solute transport and recovery in rock fractures: Implication for groundwater remediation

A better understanding of solute transport and retention mechanism in rock fractures has been challenging due to difficulty in their direct observations in microscale rough‐walled fractures. Six representative troughs in a rough‐walled fracture were selected for microscale observations of eddy formation with increasing flow velocity and its effect on spatiotemporal changes of solute concentration. This experimental study was enabled by a microscale visualization technique of micro particle image velocimetry. With increasing flow velocity (Re ≤ 2.86), no eddies were generated, and solutes along the main streamlines transported rapidly, whereas those near the wall moved slowly. A larger amount of solutes remained trapped at all troughs at Re = 2.86 than Re < 1. For Re = 8.57, weak eddies started to be developed at the troughs on the lee side, which little contributed to overall solute flushing in the fracture. Accordingly, a large of amount of water was needed for solute flushing. The flow condition of 1 < Re < 10, before a full development of eddies, was least favourable in terms of time and amount of remediation fluid required to reach a target concentration. After large eddies were fully developed at troughs on the lee side for Re = 17.13, solutes were substantially reduced by eddies with less amount of water. Fully developed eddies were found to enhance solute transport and recovery, as opposed to a general consensus that eddies trap and delay solutes. Direct inflow into troughs on the stoss side also made a great contribution to solute flushing out of the troughs. This study indicates that fully developed eddies or strong inflows at troughs are highly possible to form for Re > 10 and this flow range could be favourable for efficient remediation.     

Solute transport routing in a small stream

Abstract

The impact of pollution incidents on rivers and streams may be predicted using mathematical models of solute transport. Practical applications require an analytical or numerical solution to a governing solute mass balance equation together with appropriate values of relevant transport coefficients under the flow conditions of interest. This paper considers two such models, namely those proposed by Fischer and by Singh and Beck, and compares their performances using tracer data from a small stream in Edinburgh, UK. In calibrating the models, information on the magnitudes and the flow rate dependencies of the velocity and the dispersion coefficients was generated. The dispersion coefficient in the stream ranged between 0.1 and 0.9 m²/s for a flow rate range of 13–437 L/s. During calibration it was found that the Singh and Beck model fitted the tracer data a little better than the Fischer model in the majority of cases. In a validation exercise, however, both models gave similarly good predictions of solute transport at three different flow rates.     

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.