3.3: Pathways and processes determining chemical fate
- Page ID
- 294546
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)3.3. Pathways and processes determining chemical fate
Authors: Dik van de Meent, Michael Matthies
Reviewer: John Parsons
Leaning objectives:
You should be able to
- explain how emission of chemicals into the environment leads to exposure of ecosystems, populations, and organisms including man.
- understand how quantitative knowledge of process kinetics is of use in mathematical modeling of environmental fate
- understand the concept of characterization of fate processes in terms of rate constants and half-lives
- give examples of relevant fate processes and briefly describe them
Keywords: fate processes, degradation, transport, partitioning
Characterizing 'fate'
Chemicals can escape during all steps of their life cycle, e.g. manufacturing, processing, use, or disposal. Release of chemicals into the environment necessarily leads to exposure of ecosystems, populations, and organisms including man. Exposure assessment science seeks to analyze, characterize, understand and (quantitatively) describe the pathways and processes that link releases to exposure. Chemicals in the environment undergo various transport, transfer and degradation processes, which can be described and quantified in terms of loss rates, i.e. the rates at which chemicals are lost from the environmental compartment into which they are emitted or transferred from adjacent compartments. Exposure assessment science aims to capture the 'environmental fate' of chemicals in process descriptions that can be used in mass balance modeling, using mathematical expressions borrowed from thermodynamic laws and chemical reaction kinetics (Trapp and Matthies, 1998).
The 'fate' of a chemical in the environment can be viewed of as the net result of a suite of transport, transfer and degradation processes (see Section 3.4 on partitioning and partitioning constants, Section 3.6 on availability and bioavailability, Section 3.7 on degradation) that start to act on the chemical directly after its emission (see Section 3.2 on sources of emission) and during the subsequent environmental distribution. Environmental fate modeling (see Section 3.8 on multimedia mass balance modelling) builds on this knowledge by implementing the various degradation, transfer and transport processes derived in exposure assessment science in mathematical models that simulate 'fate of chemicals in the environment'.
First-order kinetics
In chemical reaction kinetics, the amount of chemical in a 'system' (for instance, a volume of surface water) is described by mass balance equations of the kind:
(eq. 1)
where is the rate of change (kg.s-1) of the mass (kg) of chemical in the system over time t (s), i is the input rate (kg.s-1) and k (s-1) is the reaction rate constant. Mathematically, this equation is a first-order differential equation in m, meaning that the loss rate of mass from the systemis proportional to the first power of m. Equation 1 is widely applied in description and characterization of environmental fate processes: environmental fate processes generally obey first-order kinetics, and can generally be characterized by a first-order reaction rate constant k1st:
(eq. 2)
Such loss rated equations can also be formulated in integral format, which is obtained by integration of equation (2) over time t with initial mass m0 = m(0):
(eq. 3)


As shown in Figure 1, first-order loss processes are expected to result in exponential decrease of mass from which concentration can be calculated by dividing m with the compartment volume. Using the value for t in equation 3, it follows directly that the value of
is inversely proportional to the first-order loss rate constant
:
(eq. 4)
which shows that half-life time constant, i.e. independent of the concentration of the chemical considered. This is the case for all environmental loss processes that obey first-order kinetics. First-order loss processes can therefore be sufficiently characterized by the time required for disappearance of 50% of the amount originally present.
The disappearance time DT50 is often used in environmental regulation but is only identically with the half-life if the loss process is of first order. Note that the silent assumption of constancy of half-life implies that the process considered is assumed to obey first-order kinetics.
Abiotic chemical reactions
Occurrence of true first-order reaction kinetics in chemistry is rare (see Section 3.7 on degradation). It occurs only when substances degrade spontaneously, without interaction with other chemicals. A good example is radio-active decay of elements, with a reaction rate proportional to the (first power of) the concentration (mass) of the decaying element, as in equation 3.
Most chemical reactions between two substances are of second order:
(eq. 5)
or, when a chemical reacts with itself:
(eq. 6)
because the reaction rate is proportional to the concentrations (masses) of both of the two reactants. It follows directly from equation 2. As the concentrations (masses) of both reactants decrease as a result of the reaction taking place, the reaction rate decreases during the reaction, more rapidly so at high initial concentrations. When second-order kinetics applies, half-life is not constant, but increases with ongoing reaction, when concentrations decrease. In principle, this is the case for most chemical reactions, in which the chemical considered is transformed into something else by reaction with a transforming chemical agent.
In the environment, the availability of second reactant (transforming agent) is usually in excess, so that its concentration remains nearly unaffected by the ongoing transformation reaction. This is the case, for oxidation (reaction with oxygen) and hydrolysis (reaction with water). In these cases, the rate of reaction decreases with the decreasing concentration of the first chemical only:
(eq. 7)
and reaction kinetics become practically first-order: so-called pseudo first-order kinetics. Pseudo first-order kinetics of chemical transformation processes is very common in the environment.
Biotic chemical reaction
Chemical reactions in the biosphere are often catalyzed by enzymes. This type of reaction is saturable and the kinetics can be described by the Michael-Menten kinetic model for single substrate reactions. At low concentrations, there is no effect of saturation of the enzyme and the reaction can be assumed to follow (pseudo) first order kinetics. At concentrations high enough to saturate the enzyme, the rate of reaction is independent of the concentrations (masses) of the reactants, thus constant in time during the reaction, and the reaction obeys zero-order kinetics. This is true for catalysis, where the reaction rate depends only on the availability of catalyst (usually the reactive surface area):
= constant. (eq. 8)
One could say that the rate is proportional to the zero-th power of the mass of reactant present. In case of zero-order kinetics, the half-life times are longer for greater initial concentrations of chemical.
An example of zero-order reaction kinetics is the transformation of alcohol (ethanol) in the liver. It has been worked out theoretically and experimentally that human livers remove alcohol from the blood at a constant rate, regardless the amount of alcohol consumed.
Microbial degradation
Microbial degradation (often referred to as biodegradation) is a special case of biotic transformation kinetics. Although this is an enzymatically catalysed process, the microbial transformation process can be viewed of as the result of the encounter of molecules of chemical with microbial cells, which should result in apparent second-order kinetics (first order with respect to the number of microbial cells present, and first order with respect to the mass of chemical present):
(eq. 9)
where mbio stands for the concentration (mass) of active bacteria present in natural surface water, and kdeg represents a pseudo-first-order degradation rate constant.
Advective and dispersive transport
Chemicals can be moved from one local point to another by wind and water currents. Advection means transport along the current axis whereas dispersion is the process of turbulent mixing in all directions. Advective processes are driven by external forces such as wind and water velocity, or gravity such as rain fall and leaching in soil. In most exposure models these processes are described in a simplified manner, e.g. the dispersive air plume model. An example of a first-order advective loss process is the outflow of a chemical from a lake:
(eq. 10)
where Q stands for the flow rate of lake water [m³/s] and V is the lake volume [m³]. Q/V is known as the renewal rate constant kadv of the transport medium, here water. More sophisticated hydrological, atmospheric, or soil leaching models consider detailed spatial and temporal resolution, which require much more data and higher mathematical computing effort (see sections 3.1.2 and 3.8.2).
Transfer and partitioning
Due to Fick's first law the rate of transfer through an interface between two media (e.g. water and air, or water and sediment) is proportional to the concentration difference of the chemical in the two media (see section 3.4 on partitioning, and Schwarzenbach et al., 2017 for further reading). As long as the concentration in one media is higher than in the other, the more molecules are likely to pass through the interface. Examples are volatilization of chemicals from water (to air) and gas absorption from air (to water or soil), adsorption from water (to sediments, suspended solids and biota) and desorption from sediments and other solid surfaces.
When two environmental media are in direct contact, (first-order) transfer can take place in two directions, in the case of water and air by volatilization and gas absorption: Each at a rate proportional to the concentration of chemical in the medium of origin and each with a (first-order) rate constant characteristic of the physical properties of the chemical and of the nature of the interface (area, roughness). This is known as physical intermedia partitioning (see section 3.4 on partitioning), usually represented by a chemical reaction formula:
(eq. 11)
where [M] stands for a (mass) concentration (unit mass per unit volume) and the double arrow represents forward and reverse transport. Intermedia partitioning proceeds spontaneously until the two media have come to thermodynamic equilibrium. In the state of equilibrium, forward and backward rates (here: volatilization from water to air and gas absorption from air to water) have become equal. At equilibrium, the total (Gibbs free) energy of the system has reached a minimum: the system has come to rest, so that
(eq. 12)
and the ratio of concentrations of the chemical in the two media has reached its (thermodynamic) equilibrium value, called equilibrium constant or partition coefficient (see section 3.4 on partitioning).
Challenge
Challenge to environmental chemists is to describe and characterize the various processes of chemical and microbial degradation and transformation, of intra-media transport and intermedia transfer rate constants and of equilibrium constants, in terms of (i) physical and chemical properties of the chemicals considered and (ii) of the properties of the environmental media.
References
Schwarzenbach, R.P., Gschwend, P.M., Imboden, D.M. (2017). Environmental Organic Chemistry, Third Edition, Wiley, ISBN 978-1-118-76723-8.
Trapp, S., Matthies, M. (1998). Chemodynamics and Environmental Modeling. An Introduction. Springer, Heidelberg, ISBN 3-540-63096-1.
Name and explain in your own words the essential environmentally relevant property of (pseudo) first-order kinetics.
Give examples of transformation or transport processes that obey zero-order, first-order, second-order and pseudo first-order kinetics.
Why is it useful to formulate environmental fate processes in terms of process rates and rate constants and equilibrium constants?