InterPore   >   Wikipore   >  Article

2021-04-30 00:52:29

Breakthrough curve

Application

At an aquifer scale, test materials, e.g., unreactive solute tracers, are injected into subsurface media through injection wells and the resultant plume is monitored through observation wells at the downstream of injection points. At a laboratory scale, materials are injected into a packed column or sand tank, typically using a peristaltic or syringe pump, from one side (for column test commonly from the bottom) and the effluent is collected from the other side using an automatic fraction collector. Depending on the injection strategy and the duration of monitoring, breakthrough curves can appear as continuous, single-pulse, or multi-pulse.

Concentration (shown on the vertical axis of the plot) is typically normalized to the initial concentration injected into media and the time (on the horizontal axis) is often nondimensionalized by presenting it as the number of pore volumes (PV) entered the system. A range of underlying mechanisms of material transport in porous media may be inferred from the shape of the breakthrough curve.

 

Inferring transport characteristics from the shape of breakthrough curves 

The shape of breakthrough curves may be used to deduce different underlying transport mechanisms for solutes and particles. This has especially become important for nanoparticle transport studies in recent years as there is a range of mechanisms that simultaneously affect the transport of these materials and having a knowledge of how these mechanisms change the shape of BTC may help to identify the underlying transport mechanisms This helps to select the right model for simulating the transport of materials in porous media (Babakhani et al., 2017). 

Compared to a standard conservative BTC which look like a rectangular, a round-shaped BTC which is broader at the bottom and narrower at the top generally represents a dispersive transport behaviour rather than advective. A BTC with tailing indicates detachment of retained particle or desorption of adsorbed solute during the flushing phase of single- or multi-pulse BTC. An ascending plateau of BTC signifies blocking of attachment or adsorption sites in porous media by one layer of particles or solute ions/molecules/macromolecules/polymers and thus facilitating the transport of subsequent materials entering the system. A descending plateau observed for colloidal particles and nanoparticles may point to the ripening mechanism which, in contrast to the blocking mechanism, results from further retention of subsequent particles following the deposition of initial particles. In other words, the retention of particles initially entered the system, provide more deposition sites for subsequent particles entering the system and thus the retention process is progressive with time in a BTC with an ascending plateau. 

A total reduction in the area under the BTC may indicate the irreversible retention or sorption of colloidal particles/nanoparticles or solutes, respectively. If the plateau of BTC is flat, the deposition mechanism may result from physicochemical deposition, i.e., deep interaction energy minima. This may also be a concomitant process with particle aggregation inside pores since the growth in aggregate sizes leads to deeper energy wells (Babakhani, 2019; Babakhani et al., 2018; Phenrat et al., 2009). A horizontal shift in BTC towards later times (retardation) results from equilibrium adsorption (which is mainly a solute transport behaviour) while a shift towards earlier times (accelerated BTC) may emanate from size exclusion (a particle transport behaviour) or charge exclusion (a solute transport behaviour) (Babakhani, 2019; Babakhani et al., 2017; Van Genuchten, 1981). 

Different shapes of a single-pulse breakthrough curve representing specific mechanisms of solute or particle transport in porous media. 

Number of pore volumes on the time axis of breakthrough curves

The horizontal axis of the breakthrough curve, which shows the time elapsed, is often presented as a dimensionless number called pore volume or the number of pore volumes (PV) [] determined as: 

 PV = Q.t/Vv,

where Q is the flow rate [L3.T-1], t is the time elapsed [T], and Vv is the total volume of voids in the system and can be determined by multiplying the total volume of the system by the porosity of the porous medium. For cylindrical column experiments, which aims to mimic one-dimensional transport of materials, Vv is: 

Vv=φ.π.rc2.Lc,

where φ is the porosity of the porous medium [], rc is the inner radius of the column [L], and Lc is the length of the column [L]. Using the continuity equation, PV can be determined as: 

PV = v.t/Lc

where v is the pore water velocity (also known as real water velocity, linear groundwater flow velocity, seepage water velocity, interstitial velocity, or advection velocity) [L.T-1], which is equal to porosity times the Darcy velocity (also known as the superficial velocity or approach velocity). 

 

Continuous, single-pulse, and multi-pulse injection 

If the duration of material injection is shorter than the monitoring period and the return towards the base-flow condition is observed in the effluent sample concentrations, the resultant BTC appears as a single- or multi-pulse curve, otherwise, it appears as a continuous curve. Single- or multi-pulse breakthrough curves give additional information about the desorption or detachment processes of materials based on the shape of descending limb of the curve that is overlooked in continuous injection (Mahmoudi et al., 2020). Multi-pulse breakthrough curves are mostly products of the step-wise changes in the condition of injected materials or background flow, e.g., changing the ionic strengths or pH, aimed at evaluating their variations in a single experiment (Torkzaban et al., 2010). If the flow condition (e.g., pore water velocity) is changed across different pulses, using a zero velocity should be avoided since it will result in no effluent samples during that period, and consequently zero pore volume on the horizontal axis of the breakthrough curve thereby ignoring the processes that may happen during the stop-flow period and afterwards. 

Multi-peak or multi-shoulder BTC may appear for single- or multi-pulse BTC, because of heterogeneity in the packed media, e.g., in a stratigraphic, heterogeneous packed tank, material plumes from different layers may arrive at the outlet at different times and thus appear as multiple peaks in BTC (Babakhani et al., 2018; Phenrat et al., 2010). 


  • The breakthrough curve is relevant to the transport of materials in various natural systems such as aquifer, sediments, and soils as well as engineered systems such as sand filters. It is of interest to a broad range of practitioners and researchers concerned with monitoring and delineation of the fate and transport of hazardous environmental pollutants, e.g. radionuclide and pathogens, and emerging contaminants e.g. engineered nanoparticles and nano-/micro-plastics which are released into the environment through waste streams. This is also of interest to those working on soil and groundwater remediation where agents, like engineered nanoparticles or bacteria, are intentionally injected into the subsurface environment and the plumes of both pollutants and the newly added remediation agents have to be tracked.

  • Experimentally, breakthrough curves are obtained using bench-scale packed chambers. There is also a broad range of modelling tools available to solve the advection-dispersion equation which can be used to describe various aspects of breakthrough curves. Commonly used tools to solve the advection-dispersion equation in one-, two-, and three-dimensional domains and to determine the model parameters by fitting the model to experimental data are MODFLOW (MT3D-USGS), PHREEQC, HYDRUS, CXTFIT, and PEST codes.  

  • Babakhani P, Bridge J, Doong R-a, Phenrat T. Continuum-based models and concepts for the transport of nanoparticles in saturated porous media: A state-of-the-science review. Advances in Colloid and Interface Science 2017; 246: 75-104.

    Phenrat T., Babakhani P., Bridge J., Doong R.A., and Lowry G. V., “Mechanistic, Mechanistic-Based Empirical, and Continuum-Based Concepts and Models for the Transport of Polyelectrolyte-Modified Nanoscale Zerovalent Iron (NZVI) in Saturated Porous Media”, In Nanoscale Zerovalent Iron Particles for Environmental Restoration, Phenrat T. and Lowry G. V. Springer: New York, NY 10013-1578, USA, 201-233, 2019.

  • Babakhani P. The impact of nanoparticle aggregation on their size exclusion during transport in porous media: One-and three-dimensional modelling investigations. Scientific Reports 2019; 9: 1-12.

    Babakhani P, Bridge J, Doong R-a, Phenrat T. Continuum-based models and concepts for the transport of nanoparticles in saturated porous media: A state-of-the-science review. Advances in Colloid and Interface Science 2017; 246: 75-104.

    Babakhani P, Fagerlund F, Shamsai A, Lowry GV, Phenrat T. Modified MODFLOW-based model for simulating the agglomeration and transport of polymer-modified Fe nanoparticles in saturated porous media. Environ Sci Pollut Res, 1-20, doi:10.1007/s11356-015-5193-0 2018.

    Mahmoudi D, Rezaei M, Ashjari J, Salehghamari E, Jazaei F, Babakhani P. Impacts of stratigraphic heterogeneity and release pathway on the transport of bacterial cells in porous media. Science of The Total Environment 2020; 729: 138804.

    Phenrat T, Cihan A, Kim H-J, Mital M, Illangasekare T, Lowry GV. Transport and deposition of polymer-modified Fe0 nanoparticles in 2-D heterogeneous porous media: Effects of particle concentration, Fe0 content, and coatings. Environmental Science & Technology 2010; 44: 9086-9093.

    Phenrat T, Kim H-J, Fagerlund F, Illangasekare T, Tilton RD, Lowry GV. Particle size distribution, concentration, and magnetic attraction affect transport of polymer-modified Fe0 nanoparticles in sand columns. Environmental science & technology 2009; 43: 5079-5085.

    Torkzaban S, Kim HN, Simunek J, Bradford SA. Hysteresis of colloid retention and release in saturated porous media during transients in solution chemistry. Environmental science & technology 2010; 44: 1662-1669.

    Van Genuchten MT. Non-equilibrium transport parameters from miscible displacement experiments.  1981.