Open Access Article
Sivani
Baskaran†
and
Frank
Wania
*
Department of Physical and Environmental Sciences, Environmental Sciences and Department of Chemistry, University of Toronto Scarborough, 1265 Military Trail, Toronto, Ontario M1C 1A4, Canada. E-mail: frank.wania@utoronto.ca; Tel: +1-416-287-7225
First published on 7th June 2023
The octanol–air partition ratio (KOA), also referred to as the octanol–air partition coefficient, has a wide range of applications in environmental chemistry. In this review, we explore the historical context of using octanol as a surrogate for various types of organic matter. We examine in detail the single-parameter linear free energy relationships (spLFERs) that rely on the KOA to describe partitioning equilibria between the gas phase and vegetation, soil, particles, dust, surfaces, materials, and animal tissues. We further use poly-parameter linear free energy relationships (ppLFERs) to estimate how well octanol approximates the partitioning properties of these divergent phases. While the availability of ppLFERs for many environmentally and biologically relevant phases has rendered some of the spLFERs based on the KOA largely obsolete, the KOA still serves a useful purpose as a single parameter describing the tendency of a neutral organic chemical to partition from the gas phase into a wide variety of organic phases. As such, it is a well-defined, easy-to-comprehend and experimentally accessible descriptor of compound volatility from organic phases.
Environmental significanceThere is a long tradition of using the solvent octan-1-ol as a surrogate for a large variety of organic matter when seeking to predict the phase distribution of organic compounds. While not as ubiquitous and common as the octanol–water partitioning ratio KOW, the octanol–air partitioning ratio KOA is widely used when studying the partitioning of organic compounds between organic matter and the gas phase. With the increasing availability and acceptance of methods for predicting the partitioning of organic vapours from air into soil organic matter, atmospheric particles, biological tissues, and polymers, that do not depend on a surrogate organic solvent phase, it is valid to explore whether there is still a role to play for simple linear regressions involving the logarithm of the KOA. This review concludes that there is continued need for a single parameter quantifying organic compound volatility from organic phases and the KOA is well suited to fulfill this role. |
Chemical equilibrium partition ratios, K12, hereafter simply referred to as partition ratios, describe the relative chemical concentrations in two adjacent phases at equilibrium and are dependent on environmental conditions such as temperature. In many instances, partition ratios are referred to as partition coefficients; however, because this property is influenced by external conditions (i.e., temperature), IUPAC guidelines recommend the use of the term partition ratios.10 Because partition ratios of chemicals can span many orders of magnitude, they are often expressed on a base 10 logarithmic scale (i.e., log10
K12).
In this review we examine the role of the octanol–air partition ratio (KOA) and its relevance to our understanding of environmental chemical fate processes. In previous work we explored the different methods used to measure and estimate the KOA of an organic chemical.11 Here, we summarize all of the different phase equilibria that have been described with the help of the KOA including partitioning between the gas phase and vegetation, soil, particles, dust, surfaces, materials (e.g., polyurethane), and animal tissue. We further evaluate how good a surrogate octan-1-ol is for these phases and highlight the limitations of simple regressions between equilibrium partitioning ratios and the KOA. We also explore the use of the KOA in bioaccumulation and chemical risk assessments and as a coordinate in chemical space plots. Considering the availability of tools for estimating environmentally relevant phase equilibria involving the gas phase that do not rely on the KOA, we further address the question of why the KOA still serves a useful purpose.
; e.g., Berti et al.13), and Henry's law constant in octanol (koctH, Pa m3 mol−1; Roberts14). It has also been indirectly quantified by combining either a chemical's vapour pressure with the activity coefficients at infinite dilution in octanol (e.g., Hussam and Carr15) or the octanol–water (KOW) and the air–water partition ratios (KAW; e.g., Finizio et al.16).
K OA is the concentration of a chemical in octanol divided by the concentration of a chemical in air at equilibrium (eqn (1)).
![]() | (1) |
It was only in the early 1950s however that the use of solvent-water partition ratios was re-popularized by the work of Collander.23–25 In 1947, Collandar published his work “On Lipoid Solubility” where he compared the solvent-water partition ratios for a number of different solvents including octan-1-ol, referred to as octanol henceforth, olive oil and oleic acid.23 In particular, he noted that research has yet to identify a solvent that can dissolve a variety of solutes and represent the lipid phase well.23 Collander found that the partition ratios in one alcohol/water system can be used to calculate the partition ratios in another alcohol/water system.24 Octanol was favoured as it was, unlike ethyl ether, not very soluble in water and while hydrophilic substances and bases were more soluble in octanol, hydrophobic solutes were equally soluble in ethyl ether and octanol.24 Over the next two decades, different solvents, including various alcohols, were considered for the organic reference phase.26 In some of their earliest work, Hansch and colleagues found partition ratios had additive–constituent properties — which allows for fragment-based estimation techniques.27–29 Hansch et al. commonly used octanol as the reference organic phase in these works.27,28 Hansch later noted that multiple works showed good correlation between the octanol–water partition ratio and various organic bio-polymers.30 They reported on these developments and the applicability of the octanol–water phase ratio in drug design in some of the very first volumes of the Annual Reports in Medicinal Chemistry.30,31 The use of partition ratios, particularly the KOW, in drug design took off after the publication of various works by Hansch and their colleagues, including one of the earliest review papers on partitioning, by Leo et al. in 1971.29
It is important to note that initial measurements for partition ratios would have been extremely challenging in the absence of modern analytical instruments; early methods included the use of manometers12 and coulometry.27 The development of partition chromatography, including gas chromatography in the 1940s and 1950s,32 led to more measurements and increased interest in gas–liquid partition ratios. In 1959, in an effort to improve understanding of separation processes, Pierotti et al.33,34 measured the infinite dilution activity coefficient of various solutes in multiple solvents using gas–liquid chromatography.
The first directly reported KOA value, albeit referred to as an Ostwald coefficient, appears to be by Boyer and Bircher in 1960, while investigating the solubility of various gases in different alcohols using a modified Van Slyke-Neill blood gas apparatus.12,21 Another set of KOA data was published 13 years later in an effort to understand characteristics of solvents used in gas chromatography.35 The next reported measurement of KOA was published in 1978.36 This work built on work by Battino and co-workers exploring gas–liquid partitioning of gaseous compounds, driven in part by the desire to find a solvent that could act as a model for partitioning of compounds (including anesthetics) between the gas phase and cell membranes.37,38
By the 1980s, environmental chemists had adopted the use of KOW from the pharmaceutical industry for use in predicting the environmental fate and behaviour of chemicals (e.g., Smith et al.8). Paterson et al.39 were perhaps the first to identify the KOA as a useful descriptor for plant–air partitioning and noted the dearth of experimental KOA values. At the same time, octanol–air partitioning became increasingly of interest to the physical chemistry community in the mid 1980s as Pollack et al.,40 Carr (e.g., Hussam and Carr15), and Cabani (e.g., Berti et al.13) and colleagues explored the thermodynamics of gas–liquid equilibria, reporting Gibbs energies of solutions and infinite dilution activity coefficients for solutes in octanol.
Exploration of the KOA and other gas–liquid partition ratios to characterize the behaviour of inhaled anesthetics was renewed in the early 1990s (e.g., Taheri et al.41,42 and Liu et al.43). In the same decade, Harner and Carr began independently reporting KOA values for various solutes. Harner's focus was on the use of KOA to understand the environmental partitioning behaviour of chemicals between air and various types of organic matter.44 Dallas and Carr45 were interested in understanding the effect that the mutual solubility of octanol and water has on the KOW. A head space sampling system developed by Carr and colleagues made measuring vapour–liquid equilibria easier.15,45 In the following two decades interest in the KOA has not waned and many more measurements of the KOA have been reported.11
![]() | (2) |
is often expressed in units of m3 air µg−1 particles. The dimensionless QPA value can be obtained as:![]() | (3) |
If phase X is a biological phase, e.g., the foliage of a plant, the concentration ratio is sometimes referred to as a bioconcentration factor (BCF).
Generally, the CA in eqn (1) refers to the concentration in the gas phase. However, in some cases, the total concentration in air CTA, i.e., concentration in both the gas and particle phase (CP) has been measured and is used for the ratio calculation:
![]() | (4) |
| QXA = aKOAm | (5) |
log10 QXA = m log10 KOA + log10 a | (6) |
These types of equations are called single parameter linear free energy relationships (spLFERs). The implicit assumption underlying these relationships is that phase X, or a constituent of phase X, has solvation properties that resemble those of octanol. An example, where the entirety of phase X is assumed to have octanol-resembling solvation properties, are organic surface films. More commonly, phase X has different constituents that may contribute to the uptake of chemicals. In this case, the concentration ratio may be expressed as the sum of the concentration ratios in the bulk phase constituents weighted by the constituents' abundance in the bulk phase, e.g.:
| QXA = fAX + fWXKWA + afORXKOAm + fYXKYA + aSurface/XKSurface/A | (7) |
In many cases, the assumption is then made, explicitly or implicitly, that the term afORXKmOA is larger than all of the other terms in eqn (7), in which case we obtain:
| QXA = afORXKmOA | (8) |
log10 QXA = m log10 KOA + log10 fORX + log10 a | (9) |
Examples of such composite bulk phases are soil, atmospheric particles and biological materials.
In earlier papers, when the KOA had not been well established as a concept and when few experimental KOA values were available, KOW/KAW instead of KOA was used in relationships of type 5, 6, 8, and 9, i.e., the KOA was estimated as the ratio of the KOW and the equilibrium air water partition ratio KAW. It is believed that KOW/KAW in general has far higher uncertainty than KOA values, especially if predicted properties are being used. This has been observed in some correlations between empirical KOA and KOW/KAW values with QXA.46,47
![]() | (10) |
![]() | (11) |
![]() | (12) |
Incidentally, eqn (11) corresponds to a single parameter linear free energy relationship between KXA and PL. In other words, QXA can be, and often has been, related to PL in equations of the type:
| QXA = aPL−m | (13) |
log10 QXA = −m log10 PL + log10 a | (14) |
Equations involving KOA are believed to be preferable over those involving PL, because the ratio γO/γOR is less variable for different compounds than is γOR.49 Finizio et al.16 also noted that because the PL cannot be measured directly it can incur large errors for chemicals with high melting points that are solids at environmentally relevant temperatures. However, whether spLFERs involving PL or KOA are preferable may also depend on which of these two parameters can be obtained with higher reliability.50
Eqn (10) and (11) also illustrate a potential issue with spLFERs involving the KOA (or the PL). Both KOA and KXA are inversely proportional to PL. The PL varies over many orders of magnitude between different chemicals, whereas γO and γOR vary only over a small range. KXA and KOA (and also KXA and PL) are therefore bound to be highly correlated if the dataset includes chemicals of divergent volatility. This auto-correlation issue has been discussed by Paterson and Mackay51 and will be discussed further below.
The first is that QXA does not reflect equilibrium conditions. If a chemical's KXA is large, it can take a long time to reach equilibrium. If the chemical originates in the gas phase, the measured CX is lower than it should be at equilibrium and the QXA for such high KXA compounds accordingly is lower than KXA. This leads to slopes m that are lower than 1. If, however, the chemical originates in phase X, non-achievement of equilibrium leads to a measured CA that is lower than it should be at equilibrium. QXA for high KXA compounds then is higher than KXA and slopes m are more likely to exceed 1. An example of the latter would be a dust particle that includes polymer fragments containing the chemical of interest as an additive.
The KOA threshold for non-achievement of equilibrium depends on a number of factors, such as the size, viscosity, composition and homogeneity of the octanol-resembling constituent in phase X, the mass transfer kinetics between air and that phase and the time available for equilibration. However, often several studies indicate similar threshold values. For example, vegetation is likely to reach equilibrium with gas phase chemicals with a log10
KOA below 8.52,53 On the other hand, compounds with a log10
KOA in excess of 13 are unlikely to achieve equilibrium with atmospheric particles within the atmospheric lifetime of most particles.54
Disequilibrium may also be a result of temperature differences between air and the condensed phase X. Partitioning between gas and condensed phases are generally highly temperature dependent and the low heat capacity of air means that temperature can fluctuate far more quickly in air than in phase X. This, however, tends not to be an issue in laboratory experiments, where thermal equilibrium can typically be ensured. Some chemicals may also form non-exchangeable residues in the condensed phase, e.g., in soil organic matter, and thus are not available for partitioning but can be chemically extracted during analysis.
Another reason for slopes m being lower than 1 could be the use of total air concentrations CTA instead of gas phase concentrations CA (eqn (4)versuseqn (2)). In that case, QXA is too small (relative to KXA) for low volatility chemicals that partition to atmospheric particles.
The final reason for deviations of slopes m from 1 are differences in the solvation properties of octanol and the octanol-resembling phase. Goss and Schwarzenbach48 show that a slope of 1 is obtained if γO/γOR is constant. For substances within a compound class, slopes deviating from 1 can be obtained if γO/γOR varies proportionally with KOA.
| KVA = fA + fWKWA + fLKOA | (15) |
Subsequently, the KOA has been applied repeatedly to estimate chemical uptake in different plant species (e.g., Tremolda et al.,61 Dowdy and McKone,62 Keymeulen et al.63). Tolls and McLachlan64 reported a correlation of equation type 5 between log10
KVA measured for welsh ray grass and log10
KOA with R2 values of 0.89 to 0.90. KOA values that were consistently lower than the KVA were attributed to differences in the temperature at which the two variables had been measured.64 Brown et al.65 also found linear correlations between the log10 of the whole pine needle–air partition ratio (log10
KPine/A, R2 = 0.96) and the pine needle wax–air partition ratio (log10
KWax/A, R2 = 0.70) and the log10
KOA. The regression coefficients a and m in eqn (6) are thus specific to the plant species (e.g., Tolls and McLachlan,64 Kömp and McLachlan,66 McLachlan,53 Nizzetto et al.67). The intercept should further be dependent on whether the concentration in vegetation is expressed on a fresh or dry weight basis. Most spLFER equations are for 25 °C but exploration of the interspecies variability of KVA in relation to the KOA showed that KVA is more temperature dependent than the KOA.66 In a modelling effort, Taylor et al.68 used eqn (6) with m equal to 0.7 and log10
a equal to 0.15 for clover–air partitioning as a general model for plant–air partitioning whenever no plant specific model is available. We stress that most of these studies were limited to non-polar organic contaminants.
The KOA (or KOW and KAW combined) has also been related to the partitioning between specific parts of a plant and air, including leaf–air,55,65 bark–air,69,70 tree core–air71 and cuticle–air.63 The capacity to quantify these individual partition ratios in plants has advanced passive-sampling techniques in both indoor and outdoor environments for organic pollutants. For example, in indoor air, olive trees have been used to measure polycyclic aromatic hydrocarbons72 and in outdoor environments, bark, leaf, litter, and pine needles have been used as passive samplers for PBDEs,73 chlorobenzenes,74 and various volatile organic compounds.75
In addition to partitioning between air and individual plants and plant parts, the KOA can be used to describe bulk phase partitioning between the gas phase and entire forest canopies. Horstmann and McLachlan76 measured deposition of chemicals onto deciduous and coniferous canopies in Germany and found that the distribution of chemicals with log10
KOA < 10 approached equilibrium. Accordingly, they derived equations in the form of eqn (5) to relate pseudo canopy–air partition ratios KCanopy/A with the KOA, both for a deciduous canopy (m = 0.76, log10
a = 1.15) and a coniferous canopy (m = 0.69, log10
a = 1.58). Su et al.77 conducted a similar study in a deciduous forest in Canada and derived a similar correlation (m = 0.67, log10
a = 2.04). Many of these equations were derived from datasets comprising only a small subset of structurally related chemicals (e.g., PCBs) and thus these models likely work best when applied to chemicals of similar size and structure.
The relationships above assume that chemical equilibrium is established between the vegetation and the atmospheric gas phase. However, as mentioned in Section 2.4, this is generally only a defensible assumption for relatively volatile chemicals, whereas for less volatile chemicals the kinetics of uptake and plant growth play a more important role in controlling plant bioaccumulation and bioconcentration.78 McLachlan53 presented a conceptual framework, which illustrates how measurements of plant uptake can be used to delineate the KOA ranges within which different processes, including wet and dry deposition, are dominant. Volatile chemicals with low log10
KOA (<8.5) reach equilibrium and the KVA is directly correlated with the log10
KOA. For involatile chemicals with high log10
KOA (>11), particle bound deposition is important which is apparent in another correlation of the log10 of the particle–air partition ratio log10
KPA with log10
KOA.53 For chemicals with intermediate log10
KOA (8.5–11) uptake into plants from air is not a function of KOA, but rather a function of the kinetics of the transport from air to the plant storage compartment and the time of exposure. Later modelling work52 confirmed that non-volatile chemicals (log10
KOA > 8) do not approach equilibrium between plants and air and so the application of the KOA for estimating plant BCFs and KVA is more applicable to volatile chemicals. More volatile chemicals (log10
KOA < 5) also have the potential to volatilize from a plant's stem, trunk, and leaves.79
KOC) is highly correlated with the log10
KOW, with a slope of ∼1.82KOW/KAW and the organic carbon–air partition ratio KOC/A were found to be similarly related.81 Accordingly, Hippelein and McLachlan80 extended the Karickhoff equation to soil–air partitioning:![]() | (16) |
As more directly measured KOA became available, its role in soil–air partitioning became evident and the KOA was used directly to describe soil–air partitioning.80,83 For example, volatilization fluxes of pesticides from sludge-amended soil were found to correlate with KOA and the pesticide concentration in the sludge.84,85 He et al.86 found that when using directly measured KOA values, the predicted KSA value was ∼2.7 times greater than the experimental KSA, while KSA values predicted using experimental KOW/KAW values were in better agreement with experimental KSA values. This prompted He et al.86 to develop a modified Karickhoff equation:
![]() | (17) |
![]() | (18) |
Contrary to expectations, Hippelein and McLachlan88 observed that the absorption of organic chemicals in organic matter showed dependence on relative humidity and proposed a semi-empirical equation for predicting KSA:
![]() | (19) |
is the internal energy of soil–air phase change, which also varies with RH.88 Hippelein and McLachlan recommend an A value of 2 and a value of 0.0437 for B as they found B to be relatively constant across penta- and hexachlorobenzene and various PCBs.80 As both B and
are chemical specific, or at least homologue specific, using this semi-empirical equation for predicting KSA to describe general soil–air partitioning poses some challenges with data availability.89 Instead, Davie-Martin et al.89 suggest that compound group specific multiple linear regressions using KOA, RH, and temperature be used for estimating KSA because this approach requires fewer chemical specific properties (i.e., B and
) and performed much better than the general equation. However, there are limitations to this approach, as the fOC term is not included in the regression and KSA is assumed to be linearly correlated with RH.
While the relationships above all assume that partitioning of chemicals to soils from the air is solely governed by absorption on soil organic matter, this assumption is often not valid, especially in air-dry soil and soil with a low organic matter content. Under such circumstances, adsorption to mineral surfaces can become dominant.90,91 In particular, relative humidity plays a critical role in determining the relative importance of adsorption to the mineral–air interface. If only absorption processes (i.e., sorption into soil organic matter) is considered, estimates of the KSA are likely to be underestimated.91 In other words, whenever adsorption plays a role in soil–air partitioning, using octanol as a surrogate for soil will be insufficient.91,92
and log10
PL.95
Whereas in its original formulation, the Junge–Pankow relationship was based on the assumption that chemicals adsorb onto the surface of the particles, Pankow96 showed that the use of PL is also compatible with an absorption model of particle–air partitioning, which took the form of eqn (11) and related PL to
. If chemicals are absorbed into a particle's organic matter (OM) fraction, the KOA can serve as an alternative measure for estimating
(see Section 2.3 and Finizio et al.16 for details). Assuming an octanol density (ρO) of 820 kg m−3,16 a ratio MOγO/MOMγOM equal to 1 and that all OM present in the aerosol is available for partitioning,95 simplifies eqn (11) to an equation of type 6, where log10
a equals −11.91 and m is 1.
A strong relationship between
and KOA would suggest that interactions of neutral organic chemicals with the OM in the particles is similar to that in octanol, i.e., via absorption, rather than surface adsorption.16 Indeed, KOA and
values have been found to be correlated for a number of compound classes including PAHs, PCBs, organochlorine pesticides, and PCNs16,95 and the KOA is often considered to be a better predictor than PL for particle–air partitioning,95,97,98i.e., the KOA-based model performs better than the Junge–Pankow model.95,97,99 The KOA based model is believed to act as a more universal model for
, as structurally similar compounds tended to group together when regressing
against the KOA, but this did not occur when
was regressed against PL.95,97,98
While the regressions between KOA and
can be quite strong (R2 0.60 to 0.99), slopes of the regression are often less than 1 (e.g., Falconer and Harner,97 Finizio et al.,16 Harner and Bidleman,95 Radonic et al.100). Two explanations have been proposed to explain this: octanol and OM have different solvation properties or equilibrium has not been reached (see Section 2.4). It is likely that both explanations hold true. For example, the polarity of secondary aerosol particles can differ from octanol, such that the activity ratio of a chemical in octanol and organic matter deviates from 1, thus the slope of a KOA and
regression will also differ from 1.101 Yeo et al.102 suggested that the activity ratio assumption of 1 held true for multi-ortho PCBs but not for non/mono-ortho PCBs based on regressions of KOA and PL from Harner and Bidleman.95 Götz et al.103 noted that the KOA based model tended to overestimate the
and recommended using an MO/MOM value of 0.26 in eqn (11).
The presence of soot or elemental carbon can also alter the particle–air partitioning of chemicals because then adsorption to such sorbents can occur in addition to bulk phase absorption.104 In these instances a model based solely on KOA has a tendency to underestimate the
,104–106 and a second term to describe the adsorption to soot is needed:104
![]() | (20) |
is the soot–air partitioning ratio in units of L kg−1.104 The 1012 included in the denominator is for unit conversion. However, Helm and Bidleman107 found the opposite to be true, where the
model overpredicts the
whereas using KOA alone better agreed with field measurements. The differing conclusions from the various studies do not reject the suitability of KOA as parameter for predicting
, rather suggests that better characterization of the particle phase (e.g., the polarity) is needed to better understand chemical interactions between particle and gas phase.
Dust–air partitioning can also be defined by
which has units of m3 mg−1 and is often used in place of QDust/A (Weschler and Nazaroff;109 see Section 2.1). Whether the equilibrium partitioning of a chemical with dust is kinetically limited depends on both the KOA and size and viscosity of the particle, whereby chemicals with higher log10
KOA values take longer to equilibrate as do larger particles.108 Weschler and Nazaroff109 found a correlation between measured
and KOA values—however they note that for some points in the regression the
values are calculated using KOA. For chemicals with high log10KOA values, the
or QDust/A value is less than what would be expected which is attributed to non-equilibrium conditions (Weschler and Nazaroff,109 Zhang et al.;110 see Section 2.4).
Similarly, Shoeib et al.111,112 found that the KOA model for QDust/A tended to under-predict the concentration in the dust phase, which could be attributed to differences in how some compounds (e.g., MeFOSE and EtFOSE) interact with octanol and organic matter or non-equilibrium conditions.112 It is possible that the dust particles contain varying levels of contaminants, depending on whether they originate from the physical degradation of polymers containing the chemical as additive or whether the chemical is being taken up in the particles from the gas phase.112 Indeed, this was observed in house dust samples, where the bromine concentration for BDE 209 contaminated samples varied highly.113
Using Weschler and Nazaroff's109 equations for
(i.e., eqn (3) and (8)) and assuming values for fOM,Dust and ρDust, Li et al.114 present an equation for
in the form of eqn (6), where m equals 1 and log10
a is −13.0.
Chemicals with log10
KOA values between 8 and 11, had measured
values similar to predicted values.114 Chemicals with log10
KOA values less than 8 had predicted
values that were higher than the measured value, likely because air concentration levels were underestimated.114 Predicted
values lower than measured values for chemicals with log10KOA values greater than 11 were again attributed to disequilibrium conditions.114
Recently, Wei et al.115 combined eqn (4) and (6) parameterized for both KPA and KDust as functions of KOA, to estimate a dust–total air partition ratio KDust/TA:
log10 QDust/A = (m1 − m2)log10 KOA + log10 a1 − log10 a2 − log10 TSP | (21) |
QDust/A and log10
QPA against log10
KOA, respectively, and log10
a1 and log10
a2 are the intercepts of these regressions. If TSP × QPA is much less than 1, particle–air partitioning is negligible and the use of eqn (21) is redundant.115 Wei et al.115 observed good agreement between KOA and QDust/TA for phthalates, which were estimated to have TSP × QPA values less than one. Whereas PAHs which have a wider range in TSP × QPA values, showed a poorer correlation between KOA and QDust/TA.115
However, the organic film present on many surfaces is often treated as a bulk-phase, and therefore an absorption model using KOA has been applied to estimate film–air partition ratios QFilm/A (e.g., Harner and Bidleman,95 Weschler and Nazaroff108). Implicit in these approaches is either the assumption that surface adsorption is negligible or that adsorption to the film surface and absorption in the bulk film can be collectively described with one partitioning constant. The latter can be problematic, as the time scales for equilibration are often widely different for surface and bulk phase uptake. Also, adsorption scales with surface area and absorption with mass. In some instances, the term film–air partitioning has been used to describe partitioning between thin polymer films and air (e.g., ethylene vinyl acetate116); these instances are described in Section 2.5.6.
Q
Film/A values for both indoor117,118 and outdoor119 films were found to be correlated with, or similar to, KOA values. Subsequently, studies have used organic films as passive environmental samplers relating KOA and QFilm/A using eqn (8), where m and log10
a are both 1.120,121 While the measured QFilm/A correlated well with KOA, Wu et al.122 observed very shallow regression slopes particularly for films obtained from rural environments.
Bi et al.123 reported on the
of mirrors, plates, and windows in units of m and, assuming that organic films are present on all impervious surfaces, converted these values to dimensionless QFilm/A, by estimating a film thickness and percentage of organic material in the film. These QFilm/A values were in good agreement with KOA.
Li et al.114 presented a simplified relationship between KOA and a window film–air partition ratio. As with
, estimated
values for chemicals with log10
KOA values between 8 and 11 agreed well with experimental values, chemicals with log10
KOA < 8 were over-predicted, and those with log10
KOA values > 11 were under-predicted.114
The time needed for a chemical to reach equilibrium with a surface film is dependent on both the KOA, the thickness of the film (X; m), and the mass transfer coefficient from bulk air onto the film (MTCS in units of m h−1):108
![]() | (22) |
As discussed in Section 2.4, it can take several months for surface–air partitioning of chemicals with very high log10
KOA values to reach equilibrium. For chemicals with log10
KOA > 12, this can be upwards of 1 year.124
While modelling film growth by SVOCs deposition, Weschler and Nazaroff125 noted that the majority of chemicals in organic film have log10
KOA values between 10 and 13. Because chemicals with low KOA values equilibrate quickly between organic films and air, films can become enriched in chemicals with higher log10
KOA over time.125
The KOA has also been used to estimate a time-dependent uptake coefficient (ka) for gaseous compounds by liquid organic films:
![]() | (23) |
Ockenden and colleagues127 found that the KOA was linearly correlated with the sampling rate of PCBs into triolein containing semi-permeable membrane devices (SPMD) used as passive air samplers (PASs), thus the KOA could be used to estimate the amount of time required by non-polar substances to equilibrate with the PAS.
Multiple studies have shown correlations between KOA and the partitioning between air and polymer material used in other PASs, including polyurethane foam,128 low density polyethylene,129 ethylene vinyl acetate,130 and activated carbon and florisil.131 De Coensel132 noted that given the strong correlation between the KOA and KPDMS/A (polydimethyl siloxane–air partition ratio) for p-dichlorobenzene and naphthalene, KPDMS/A could even be used to estimate KOA. However such an approximation is likely to be limited to non-polar substances as the interactions of polar substances with PDMS and octanol is not the same.
The KOA has also been correlated with partitioning between various textile materials such as cotton and air, generally showing fairly strong correlations.133 Yu et al.134 observed that the correlation between the cotton–air partition ratio QCotton/A and KOA is specific to different homologous compound groups. However, Saini et al.135 found a low correlation between KOA and measured QCotton/A which was attributed to the irregularity and variance in the material itself, as polyester–air partition ratios KPolyester/A were more strongly correlated with KOA.
Won et al.136 also observed linear relationships between KOA and partition ratios between carpet, vinyl, and drywall material and air under different conditions. Partitioning between latex paint and air was also found to be correlated with KOA.137
Reppas-Chrysovitsinos et al.138 compiled published material–air partition ratios (KMA) and presented a general single parameter equation based on eqn (6), where m is equal to 1 and log10
a equal to −1.22. They also present KOA-based equations to estimate specific material–air partitioning ratios for polyethylene, polyurethane, PDMS, carbohydrates, polyoxymethylene (POM), and nylons.138
Arguably, adsorption is the most relevant sorption process to describe partitioning to non-porous and non-organic material and octanol is unlikely to be a relevant surrogate for surfaces and materials whenever adsorption is the dominant uptake process.
| KXA = KOA(fneutral lipids + 0.3fphosolipids) + KAW(fwater + 0.7fphospholipids) | (24) |
Hau and Connell correlated odour thresholds144 and nasal pungency thresholds145 of chemicals with KOW and KAW but not with KOA. This is likely because methods for measuring and estimating KOW and KAW were well-established and there were few measurements of KOA at the time. Later, Hau et al.146 showed biophase–air partitioning ratios to be related to KOA and as a result the minimum alveolar concentration (MAC), which describes the potency of anaesthetics, is also related to KOA. Slightly different correlations between KOA and MAC were observed for alkanes and alcohols, however if the MAC of chemicals with very high KOW values is corrected with the blood/brain partition ratio, then the regressions converge, and a single equation can describe the relationship between MAC and KOA:146
log10 MACcorr = −1.10 log10 KOA + 2.41 | (25) |
Nielsen et al.147 found the KOA to be a good proxy for sensory irritant receptor–air partitioning such that it can be used to estimate RD50, the concentration at which a chemical causes a 50% decrease in respiratory frequency in an air-breathing organism:
log10 RD50 = 0.8361 log10 KOA − 6.0879 | (26) |
Raines et al.148 similarly found a correlation between KOA and the concentration of aromatic inhaled drugs needed to inhibit the N-methyl-D-aspartate receptor by 50%.
More recently, the KOA has also been utilized in food science to estimate release behaviour of flavour compounds from an octanol–water emulsion.149,150 Measured blood–air and milk–air partition ratios151 and feces–air partition ratios152 were also shown to correlate with the KOA. KOA has also been used as a surrogate for lipid–air partitioning to predict respiratory elimination rates153 and human skin uptake108 in environmental chemical fate models. Kelly and Gobas154 noted that the estimated biomagnification factor of persistent neutral organic chemicals for wolves was correlated with KOA.
Bioaccumulation assessments based solely on octanol–water partition ratios (KOW) fail to take into account the role of respiratory elimination in air-breathing organisms.153 In fact, chemicals with log10
KOW values less than 5, and thus not considered bioaccumulative by various regulatory agencies (e.g., the Government of Canada155), can have the potential to biomagnify in air-breathing organisms. Modelling results from Gobas et al.153 indicate that persistent chemicals with log10
KOA greater than 5 and log10
KOW greater than 2 have the potential to bioaccumulate in terrestrial food chains. This was further supported by modelling work by Czub and McLachlan52 who found that in humans chemicals with log10
KOA between 6 and 10 and log10
KOW value between 2 and 9 had high environmental bioaccumulation potentials and Kelly et al.156 who noted that chemicals with log10
KOA values greater than 6 and log10
KOW values greater than 2 can biomagnify in air-breathing organisms. Of the chemicals listed on Canada's 2003 Domestic Substances List, 57% have an estimated log10
KOW value between 2 and 5 and an estimated KOA, value greater 5—and thus are potentially bioaccumulative.153 Providing evidence beyond modelling studies, Moses et al.157 found lipid-normalized concentrations of dieldrin, heptachlor epoxide and β-hexachlorocyclohexane to be 5–14 times higher in seals than in fish. These chemicals have log10
KOW values around 5 or lower and log10
KOA values greater than 8.
More recent modelling studies on compounds released indoors have also shown that inhalation exposure is the primary exposure route of persistent volatile organic compounds with log10
KOA values less than 6,158 while for chemicals with log10
KOA greater than 6, non-dietary ingestion (i.e., hand-to-mouth) is the primary exposure pathway, particularly in children.158,159
QXA (where X is an organic phase) and log10
KOA, it is not really possible to compare these relationships with each other, for example to assess whether the solvation properties of a particular phase more closely resemble those of octanol than those of another. This is because of divergent datasets used in those studies that include different types and numbers of chemicals and because of the unknown and variable uncertainty of the KOA and QXA values used in these relationships. Also, as mentioned before, it is not always assured that the QXA data refer to equilibrium conditions.
However, we can test how good a surrogate KOA is for other partitioning equilibria between an organic phase X and air by comparing KXA values that have either been obtained with a KOA-based spLFER prediction or by a poly-parameter linear free energy relationship (ppLFER) that has been calibrated directly for phase X. In some cases (i.e., plant cuticles, storage lipid, membrane lipid), we can directly compare the KOA with a ppLFER-predicted KXA. Furthermore, we sometimes have to rely on thermodynamic triangles to estimate a KXA from ppLFER-predicted KXW and KAW values (natural organic matter NOM, storage lipids, membrane lipids, tissue, POM). To avoid the uncertainty in KOA influencing the comparison, we can also predict the KOA at the appropriate temperature with ppLFERs.160,161 This approach is thus similar to the one used by Endo et al.,143 when probing the extent to which n-octanol is a good surrogate for various biomaterials. Table 1 details the partition ratios that are being compared. We used 1316 organic chemicals with complete sets of experimental solute descriptors, as supplied by the UFZ-LSER database.162 We note that this set of chemicals is limited to generally small and primarily non-polar substances. More complex multifunctional chemicals or chemicals with unusual substitutions, such as the polyflourinated alkyl substances, are not well-represented.
| Phase X | ppLFER predicted KXA | log10 KOA spLFER based prediction |
|---|---|---|
| NOM | Average log10 KNOM/A for 10 different organic matter systems including humic acids, fulvic acids, and natural organic matter, predicted with ppLFER by Niederer et al.167 |
log10 KOC/A predicted with spLFER by Borisover and Graber81 using KOA at 15 °C |
| OC | log10 KOC/A calculated as log10 KWA + log10 KOC, where KOC is predicted with ppLFER by Nguyen et al.168 |
log10 KOC/A predicted with spLFER by Borisover and Graber81 using KOA at 25 °C |
![]() |
||
| Plant cuticles | log10 KPlant/A predicted with ppLFER by Platts and Abraham163 |
log10 KOA at 25 °C |
![]() |
||
| Storage lipid | log10 KStorage-lipid/A calculated as log10 KWA + log10 KStorage-lipid/W where KStorage-lipid/W is predicted with ppLFER by Geisler et al.164 |
log10 KOA at 37 °C |
![]() |
||
| Membrane lipid | log10 KMembrane-lipid/A calculated as log10 KWA + log10 KMembrane-lipid/W where KMembrane-lipid/W is predicted with ppLFER by Endo et al.165 |
log10 KOA at 37 °C |
| Tissue | log10 KTissue/A calculated as log10 KTissue/W + log10 KWA where KTissue/W is predicted with a composition model by Endo et al.143 using ppLFERs for KStorage-lipid/W by Geisler et al.,164KMembrane-lipid/W by Endo et al.,165KAlbumin/W by Endo and Goss,169KProtein/W by Endo et al.170 |
log10 KTissue/A predicted using a fractional composition model with KOA to 37 °C; KTissue/A = (fStorage-lipid + fMembrane-lipid)KOA + fWKWA |
| Aerosols | Average of log10 KPA for aerosol sampled in Dübendorf fall & winter, Aspvretn, and Roost, predicted using ppLFERs by Arp et al.171 |
log10 KPA predicted with spLFER by Harner and Bidleman95 using KOA at 15 °C |
| POM | log10 KPOM/A calculated as log10 KPOM/W + log10 KWA where KPOM/W is predicted with ppLFER by Endo et al.172 |
log10 KPOM/A predicted with spLFER by Reppas-Chrysovitsinos et al.138 using KOA at 25 °C |
| PDMS | log10 KPDMS/A predicted with ppLFER by Sprunger et al.173 |
log10 KPDMS/A predicted with spLFER by Reppas-Chrysovitsinos et al.138 using KOA at 25 °C |
In general, there are very strong correlations (R2 > 0.89) between the KOA based spLFER and ppLFER estimated values of log10
KXA (Fig. 1). Such good correlations can at least in part be attributed to a statistical issue discussed by Paterson and Mackay,51 who assert that apparently significant correlations can be obtained by inadvertently correlating a quantity with itself. In the current case, log10
KXA and log10
KOA are almost inevitably correlated, because of the presence of the “solubility in air” in their denominator (eqn (10) and (11)). The vapour pressure, or “solubility in air”, of substances ranges over many orders of magnitude, whereas the solubilities of substances in organic phases vary little and thus contribute little to the variation in partition ratios between chemicals. For similar reasons the KOA is very highly correlated with the liquid saturation vapour pressure.50 In order to eliminate this auto-correlation, we also plot the difference between ppLFER and spLFER estimates of KXA against the log10
KOA used to generate the spLFER estimate (Fig. 2).
![]() | ||
Fig. 2 The difference between ppLFER and spLFER estimates for log10 KXA plotted against the log10 KOA used to derive the spLFER estimate. | ||
Thus, if the KOA is a good surrogate for a partitioning equilibrium, we expect to see a regression with a slope m of near unity in Fig. 1 (ref. 47) and a horizontal regression line in Fig. 2. By normalizing ppLFER estimated values of KXA by spLFER estimates of KXA we can better visualize the variation in KXA with KOA that is not easily observed in Fig. 1. For example, a negative slope in Fig. 2 indicates that as the KOA of a chemical increases, the KOA based estimate (i.e., the spLFER estimate) is greater than the ppLFER estimate. This suggests that when the slope is negative (<0), the activity ratio (γO/γX) is greater than 1 whereas a positive slope suggests that the activity ratio is less than 1. In Table 2 we present the slopes and intercepts of the regressions and coefficients of determination (R2) of both plots and the root mean squared error (RMSE) between the spLFER and ppLFER-based estimates.
| Phase X | Fig. 1 | Fig. 2 | RMSE | ||||
|---|---|---|---|---|---|---|---|
| Slope | Intercept | R 2 | Slope | Intercept | R 2 | ||
| NOM | 0.78 | −0.04 | 0.98 | −0.24 | 0.14 | 0.79 | 0.42 |
| OC | 0.83 | −0.22 | 0.89 | −0.18 | −0.08 | 0.25 | 1.02 |
| Plant cuticle | 1.07 | −0.19 | 0.99 | 0.07 | −0.19 | 0.28 | 0.39 |
| Storage lipids | 0.96 | −0.31 | 0.97 | −0.04 | −0.31 | 0.06 | 0.50 |
| Membrane lipids | 0.99 | −0.41 | 0.97 | −0.01 | −0.41 | 0.00 | 0.53 |
| Tissue | 1.12 | 0.97 | 0.99 | 0.12 | 0.52 | 0.95 | 0.26 |
| POM | 0.93 | −0.47 | 0.92 | −0.07 | −0.42 | 0.07 | 0.87 |
| PDMS | 0.84 | 0.58 | 0.97 | −0.16 | 0.72 | 0.57 | 0.45 |
| Aerosols | 0.93 | −0.50 | 0.99 | −0.07 | 0.37 | 0.31 | 0.33 |
Slopes near 1 in Fig. 1 and near 0 in Fig. 2 and RMSE values of 0.50 and 0.53 respectively, indicate that the KOA value is a good surrogate for the partitioning between storage and membrane lipids and the gas phase and a very good predictor for partitioning into tissue. The results for PDMS are similar to those of storage and membrane lipids. The KOA also appears to be a reasonable surrogate for the partitioning between aerosols95 and plant cuticles163 and the gas phase, with RMSE values of 0.33 and 0.39 respectively. In the case of partitioning to soil organic matter from the gas phase, the spLFERs based on the KOA have the highest RMSE and tend to give predictions that are biased low for compounds with high KOA values, indicated by the slopes in the plots in Fig. 2 being negative and the large scatter below the regression line. To a smaller extent that is also the case for the KOA-based spLFERs for POM and PDMS138 and the storage164 and membrane lipids.165 The opposite is the case for the KOA and KAW-based fractional composition model for KTissue/A where a positive slope in Fig. 2 suggest that relying on KOA to estimate partitioning into tissues leads to an overestimation for low volatility compounds.
Major progress in predicting environmental phase equilibria has been made in the decades since the KOA has first been introduced. In particular, poly-parameter linear free energy relationships (ppLFERs) have become a viable method for predicting partition ratios involving environmental phases, mostly, because (i) calibrations have been performed to determine system constants (the regression constants of the ppLFER equations) for most environmentally relevant phase equilibria or can be estimated,162,174 (ii) solute descriptors for numerous contaminants have been experimentally determined.175 When carefully calibrated system constants and experimental solute descriptors are combined, ppLFERs can achieve prediction accuracy on the order of 0.3
log10 units.176 Even for substances, for which no experimental solute descriptors exist, those descriptors can now be predicted with the help of QSPRs162,177 and still result in a prediction accuracy for environmental phase equilibria better than one order of magnitude.176KOA-based spLFERs would be hard pressed to match such accuracy, even if they were only applied to substances appropriate for a particular regression, i.e., similar to the substances used in the calibration of that spLFER. In other words, there is no longer the need for a parameter such as KOA to predict contaminant partitioning from the gas phase into organism lipids,164 soil organic matter,167 atmospheric particles,171 or sampling sorbents.178
The merit of ppLFER extends beyond the provision of estimates that tend to be more precise than estimates from KOA-based spLFERs. By describing the different intermolecular interactions involved in partitioning, ppLFERs also provide a mechanistic understanding of organic chemical solvation in various phases.179
The question therefore arises whether there is still a need for the KOA. The answer is two-fold. First there are still some partitioning equilibria between air and organic phases, which need to be predicted with the help of KOA-based spLFERs. Second, the utility of the KOA is not limited to the prediction of specific environmentally relevant phase equilibria.
While many environmentally relevant phase equilibria can now be predicted directly, the availability of tools to predict partition ratios at a variety of temperatures is much more limited. For example, while ppLFERs for predicting gas particle partitioning at 15 °C exist,171 atmospheric temperatures range over more than 60 K. The KOA for many contaminants has been determined as a function of temperature or can be predicted using a ppLFER for the internal energy of phase transfer from octanol to gas phase
.161,180 It is thus possible to estimate an environmentally relevant equilibrium partition ratio at any temperature by either (i) applying KOA-based spLFERs by inserting (empirical or predicted) KOA values at the relevant temperatures or by (ii) adjusting the environmentally relevant partition ratio to the desired temperature using van't Hoff's equation and the (empirical or predicted)
. Implicit in either approach is that the temperature dependence of the partitioning between the environmentally relevant phase and the gas phase is well described by the temperature dependence of octanol–air partitioning. Admittedly, one can obviously also adjust an equilibrium partitioning ratio derived from a ppLFER using such a
. Note, since the KOA and many other environmentally relevant phase equilibrium ratios are commonly defined using molar volumetric concentration in the gas phase (rather than a partial pressure), the
needs to be used in the van't Hoff equation.
can be obtained from the enthalpy of phase transfer
by deducting the product of the ideal gas constant and absolute temperature.
There is need to further explore the extent of agreement between the internal energy of phase transfer between various environmentally relevant phases X and the gas phase
and
. This is currently limited by the availability of measurements for KXA at different temperatures, but there are indications that
can deviate from
.66
• An equilibrium partition ratio involving a real environmental organic phase is not a good choice because environmental phases are of variable composition and therefore not sufficiently well-defined.
• While the saturation vapour pressure is the traditional parameter for expressing compound volatility, the solid state vapour pressure is not an appropriate parameter to describe the volatility of a substance from any organic phase other than its own pure solid phase, because of the strong influence of the lattice energy.
• Although this can be in principle corrected by using the saturation vapour pressure of the sub-cooled liquid, it has the drawback that this parameter cannot be determined directly by experiment. It also implies that the organic reference phase is a different one for every substance.
• While other organic solvents or polymers may be equally good (and possibly even better) surrogates of environmental organic phases relative to octanol, the knowledge of how good a surrogate octanol is, is much better than for any other organic solvent (see Section 3).
As such, the KOA can still play an important role as a measure of volatility from organic phases. It is an easy to comprehend, single chemical characteristic that can be used to judge a chemical's propensity to undergo a very wide range of environmental processes. Just by knowing a compound's KOA, it is possible to infer a tremendous amount of information on its environmental behaviour. This is illustrated in Fig. 3. For example, instead of requiring separate parameters to assess a chemical's equilibrium phase distribution in atmosphere and soils, its susceptibility to evaporate or be exhaled, and the importance of processes controlling its deposition to plants, this can all be deduced from the compound's KOA. At the same time, it can inform you about the suitability of a sampling sorbent for the compound and which human exposure pathways are likely to be relevant.
![]() | ||
Fig. 3 Relationship between KOA range and the environmental fate of common contaminants. Vertical lines are indicative of the log10 KOA, blue and brown bars indicate the relative dominant phase a chemical will be present in,52,53,181–183 red bars indicate the bioaccumulation potential of chemicals within a KOA range,52,158,181,184–186 and green bars indicate the KOA range of various compound groups. PCBs: polychorinated biphenyls; PAHs: polycyclic aromatic hydrocarbons; PBDEs: polybrominated diphenyl ethers; OCPs: organochlorine pesticides; cVMS: cyclic volatile methylsiloxanes; nPFAS: per- and poly-fluoroalkyl substances. | ||
Another example for the utility of such a parameter is that it can serve as a volatility threshold for classification, prioritization and possibly even regulation of chemicals (Section 2.5.7). For example, KOA is being used as a threshold for the potential to bioaccumulate in air-breathing species.187
Another advantage of KOA is the possibility to use it in analyses that require the display of model results as a function of the chemical partitioning space. By taking advantage of octan-1-ol being able to represent a wide variety of organic phases reasonably well, the KOA facilitates the display of model-derived contaminant behaviour in a chemical partitioning space defined by coordinates of hydrophobicity and volatility. For example, when Czub et al.188 plotted the model calculated combined Arctic Contamination and BioAccumulation Potential (AC-BAP) of organic chemicals as a function of KOA and KAW, octan-1-ol simultaneously serves to approximate the solvation properties of soil and sedimentary organic matter, forest canopies, organic aerosol particles and lipids in a variety of organisms (plants, food animals, humans). Similarly, in Li et al.'s158 analysis of human exposure pathway, the KOA is used simultaneously to estimate partition equilibria between a very large variety of phases and the gas phase. This includes soil and sediment organic matter, atmospheric particles, house dust, various indoor surfaces, plant and animal lipids and human skin. While it would be possible to perform the type of model calculations performed by Czub et al.188 and Li et al.158 using partitioning ratios that are not estimated from KOA, it would be much more difficult to perform the systematic investigation of how the model results depend on the partitioning characteristics of an organic chemical without relying on partitioning ratios involving octan-1-ol. It would clearly be impossible to display them in a two-dimensional partitioning space.
Footnote |
| † Current affiliation and address: Department of Environmental Engineering, Norwegian Geotechnical Institute, Sandakerveien 140, 0484, Oslo, Norway. |
| This journal is © The Royal Society of Chemistry 2023 |