Difference between revisions of "User:Jhurley/sandbox"

From Enviro Wiki
Jump to: navigation, search
(Groundwater Flow)
(Introduction to LEAF)
 
Line 1: Line 1:
Groundwater migrates from areas of higher [[wikipedia: Hydraulic head | hydraulic head]] toward lower hydraulic head, transporting dissolved solutes through the combined processes of [[wikipedia: Advection | advection]] and [[wikipedia: Dispersion | dispersion]].  Advection refers to the bulk movement of solutes carried by flowing groundwater. Dispersion refers to the spreading of the contaminant plume from highly concentrated areas to less concentrated areas.  In many groundwater transport models, solute transport is described by the advection-dispersion-reaction equation in which dispersion coefficients can be calculated as the sum of molecular diffusion, mechanical dispersion, and macrodispersion.   
+
==PFAS Leaching Characterization with the Leaching Environmental Assessment Framework (LEAF)==
 
+
[[Wikipedia: Firefighting_foam#Synthetic_foams | Aqueous film-forming foams (AFFFs)]] are a major source of [[Perfluoroalkyl and Polyfluoroalkyl Substances (PFAS) | per- and poly-fluoroalkyl substances (PFAS)]] impacts in soil and groundwater. Standardized tools are needed to rapidly assess the potential for retention, leaching, and transport of PFAS from the source zone to downgradient regions, so that this information can be applied towards critical facets of site management such as prioritizing PFAS-impacted sites for further investigation and remediation. Existing standard leaching methods were developed prior to concerns regarding PFAS. Therefore, studies are needed to ensure that leaching methods are compatible for use with PFAS and that resulting data are representative of the risk of PFAS leaching at impacted sites.   
 
<div style="float:right;margin:0 0 2em 2em;">__TOC__</div>
 
<div style="float:right;margin:0 0 2em 2em;">__TOC__</div>
  
 
'''Related Article(s):'''
 
'''Related Article(s):'''
  
*[[Dispersion and Diffusion]]
+
*[[Perfluoroalkyl and Polyfluoroalkyl Substances (PFAS)]]
*[[Sorption of Organic Contaminants]]
+
*[[PFAS Sources]]
*[[Plume Response Modeling]]
+
*[[PFAS Transport and Fate]]
  
'''CONTRIBUTOR(S):''' [[Dr. Charles Newell, P.E.|Dr. Charles Newell]] and  [[Dr. Robert Borden, P.E.|Dr. Robert Borden]]
+
'''Contributors:''' Dr. Jennifer L. Guelfo, Dr. David Kosson, Dr. Andy Garrabrants, Ms. Fangfei Liu, Mr. Darlington Yawson, Dr. Md. Isreq Real
  
'''Key Resource(s):'''
+
'''Key Resources:'''
 +
*Development of Leaching Tests for Materials Containing SVOCs and PFAS, EPA 600/R-23/382<ref name="GarrabrantsEtAl2024"/>
  
*[http://hydrogeologistswithoutborders.org/wordpress/1979-english/ Groundwater]<ref name="FandC1979">Freeze, A., and Cherry, J., 1979. Groundwater, Prentice-Hall, Englewood Cliffs, New Jersey, 604 pages. Free download from [http://hydrogeologistswithoutborders.org/wordpress/1979-english/ Hydrogeologists Without Borders].</ref>, Freeze and Cherry, 1979.
+
*[https://www.epa.gov/hw-sw846/leaching-environmental-assessment-framework-leaf-methods-and-guidance Leaching Environmental Assessment Framework (LEAF) Methods and Guidance] (EPA website)
*[https://gw-project.org/books/hydrogeologic-properties-of-earth-materials-and-principles-of-groundwater-flow/ Hydrogeologic Properties of Earth Materials and Principals of Groundwater Flow]<ref name="Woessner2020">Woessner, W.W., and Poeter, E.P., 2020. Properties of Earth Materials and Principals of Groundwater Flow, The Groundwater Project, Guelph, Ontario, 207 pages. Free download from [https://gw-project.org/books/hydrogeologic-properties-of-earth-materials-and-principles-of-groundwater-flow/ The Groundwater Project].</ref>, Woessner and Poeter, 2020.
 
  
==Groundwater Flow==
+
==Introduction to LEAF==
[[File:Newell-Article 1-Fig1r.JPG|thumbnail|left|400px|Figure 1. Hydraulic gradient (typically described in units of m/m or ft/ft) is the difference in hydraulic head from Point A to Point B (ΔH) divided by the distance between them (ΔL). In unconfined aquifers, the hydraulic gradient can also be described as the slope of the water table (Adapted from course notes developed by Dr. R.J. Mitchell, Western Washington University).]]
+
The [https://www.epa.gov/ U.S. Environmental Protection Agency (EPA)] [https://www.epa.gov/hw-sw846/leaching-environmental-assessment-framework-leaf-methods-and-guidance Leaching Environmental Assessment Framework (LEAF)] is a suite of standardized test methods for evaluating contaminant release from solids under environmentally relevant conditions (Table 1). The four leaching methods within LEAF were originally validated for inorganic constituents<ref>Garrabrants, A.C., Kosson, D.S., Stefanski, L., DeLapp, R., Seignette, P.F.A.B., van der Sloot, H.A., Kariher, P., Baldwin, M., 2012. Interlaboratory Validation of the Leaching Environmental Assessment Framework (LEAF) Method 1313 and Method 1316, EPA/600/R-12/623, U.S. Environmental Protection Agency, Air Pollution and Control Division. [[Media: EPA 600_R-12_623.pdf | Free Download EPA 600/R-12/623]]</ref><ref>Garrabrants, A.C., Kosson, D.S., DeLapp, R., Kariher, P., Seignette, P.F.A.B., van der Sloot, H.A., Stefanski, L., Baldwin, M., 2012. Interlaboratory Validation of the Leaching Environmental Assessment Framework (LEAF) Method 1314 and Method 1315, EPA/600/R-12/624, U.S. Environmental Protection Agency, Air Pollution and Control Division. [[Media: EPA 600_R-12_624.pdf | Free Download EPA 600/R-12/624]]</ref> and included as standard leaching methods EPA 1313 – 1316 within Update V of SW-846<ref>USEPA, 2026. Hazardous Waste Test Methods / SW-846. [https://www.epa.gov/hw-sw846 USEPA SW-846 website]</ref>. To address the need for standardized tests to evaluate PFAS leaching and mobility, LEAF methods have been optimized and demonstrated for use with PFAS (Methods 1313A-1316A)<ref name="GarrabrantsEtAl2024">Garrabrants, A.C., Liu, F., Warne, R., DeLapp, R., Brown, L., Rubin, Z., Yawson, D., Kosson, D.S., Guelfo, J.L., Real, M.I., van der Sloot, H.A., Touati, A., Thorneloe, S., 2024. Development of Leaching Tests for Materials Containing SVOCs and PFAS, EPA 600/R-23/382, USEPA, Washington, D.C. [[Media: EPA 600_R-23_382.PDF | Free Download EPA 600/R-23/382]]</ref>. This article will focus on Methods 1313A, 1314A, and 1316A. Demonstration of Method 1315A for compacted granular materials (including concrete and asphalt) is ongoing.
Groundwater flows from areas of higher [[wikipedia: Hydraulic head | hydraulic head]] (a measure of pressure and gravitational energy) toward areas of lower hydraulic head (Figure 1). The rate of change (slope) of the hydraulic head is known as the hydraulic gradient. If groundwater is flowing and contains dissolved contaminants it can transport the contaminants by advection from areas with high hydraulic head toward lower hydraulic head zones, or “downgradient”.
 
  
==Darcy's Law==
+
{| class="wikitable" style="float:right; margin-left:10px;"
{| class="wikitable" style="float:right; margin-left:10px;text-align:center;"
+
|+Table 1. EPA SW-846 methods that comprise the LEAF framework
|+Table 1. Representative values of total porosity (''n''), effective porosity (''n<sub>e</sub>''), and hydraulic conductivity (''K'') for different aquifer materials<ref name="D&S1997">Domenico, P.A. and Schwartz, F.W., 1997. Physical and Chemical Hydrogeology, 2nd Ed. John Wiley & Sons, 528 pgs. ISBN 978-0-471-59762-9.  Available from: [https://www.wiley.com/en-us/Physical+and+Chemical+Hydrogeology%2C+2nd+Edition-p-9780471597629 Wiley]</ref><ref>McWhorter, D.B. and Sunada, D.K., 1977. Ground-water hydrology and hydraulics. Water Resources Publications, LLC, Highlands Ranch, Colorado, 304 pgs. ISBN-13: 978-1-887201-61-2 Available from: [https://www.wrpllc.com/books/gwhh.html Water Resources Publications]</ref><ref name="FandC1979" />
 
|-
 
!Aquifer Material
 
!Total Porosity<br /><small>(dimensionless)</small>
 
!Effective Porosity<br /><small>(dimensionless)</small>
 
!Hydraulic Conductivity<br /><small>(meters/second)</small>
 
|-
 
| colspan="4" style="text-align: left; background-color:white;" |'''Unconsolidated'''
 
 
|-
 
|-
|Gravel||0.25 — 0.44||0.13 — 0.44||3×10<sup>-4</sup> — 3×10<sup>-2</sup>
+
!Method
 +
!Description
 
|-
 
|-
|Coarse Sand||0.31 — 0.46||0.18 — 0.43||9×10<sup>-7</sup> — 6×10<sup>-3</sup>
+
| 1313 || Liquid-solid partitioning as a function of '''''extract pH''''' using a parallel batch extraction (i.e., equilibrium) procedure (Figure 1)
 
|-
 
|-
|Medium Sand||—||0.16 — 0.46||9×10<sup>-7</sup> — 5×10<sup>-4</sup>
+
| 1314 || Liquid-solid partitioning as a function of '''''liquid-solid ratio (L/S)''''' for constituents in solid materials using an up-flow '''''percolation''''' column procedure (Figure 3)
 
|-
 
|-
|Fine Sand||0.25 — 0.53||0.01 — 0.46||2×10<sup>-7</sup> — 2×10<sup>-4</sup>
+
| 1315 || '''''Mass transfer rates''''' of constituents in monolithic or compacted granular materials using a semi-dynamic tank leaching procedure
 
|-
 
|-
|Silt, Loess||0.35 — 0.50||0.01 — 0.39||1×10<sup>-9</sup> — 2×10<sup>-5</sup>
+
| 1316 || Liquid-solid partitioning as a function of '''''L/S''''' using a parallel batch extraction (i.e., '''''equilibrium''''') procedure (Figure 2)
 
|-
 
|-
|Clay||0.40 — 0.70||0.01 — 0.18||1×10<sup>-11</sup> — 4.7×10<sup>-9</sup>
+
| colspan="2" style="background:white;" | Note: Text shown in '''''bold''''' indicates primary condition evaluated in each method.
|-
 
| colspan="4" style="text-align: left; background-color:white;" |'''Sedimentary and Crystalline Rocks'''
 
|-
 
|Karst and Reef Limestone||0.05 — 0.50||—||1×10<sup>-6</sup> — 2×10<sup>-2</sup>
 
|-
 
|Limestone, Dolomite||0.00 — 0.20||0.01 — 0.24||1×10<sup>-9</sup> — 6×10<sup>-6</sup>
 
|-
 
|Sandstone||0.05 — 0.30||0.10 — 0.30||3×10<sup>-10</sup> — 6×10<sup>-6</sup>
 
|-
 
|Siltstone||—||0.21 — 0.41||1×10<sup>-11</sup> — 1.4×10<sup>-8</sup>
 
|-
 
|Basalt||0.05 — 0.50||—||2×10<sup>-11</sup> — 2×10<sup>-2</sup>
 
|-
 
|Fractured Crystalline Rock||0.00 — 0.10||—||8×10<sup>-9</sup> — 3×10<sup>-4</sup>
 
|-
 
|Weathered Granite||0.34 — 0.57||—||3.3×10<sup>-6</sup> — 5.2×10<sup>-5</sup>
 
|-
 
|Unfractured Crystalline Rock||0.00 — 0.05||—||3×10<sup>-14</sup> — 2×10<sup>-10</sup>
 
 
|}
 
|}
In&nbsp;unconsolidated&nbsp;geologic settings (gravel, sand, silt, and clay) and highly fractured systems, the rate of groundwater movement can be expressed using [[wikipedia: Darcy's law | Darcy’s Law]]. This law is a fundamental mathematical relationship in the groundwater field and can be expressed this way:
 
  
[[File:Newell-Article 1-Equation 1rr.jpg|center|500px]]
+
==Method Development for PFAS==
 +
Complete details of the development of LEAF Methods 1313A, 1314A, and 1316A for use with PFAS are available in Garrabrants ''et al.'', 2024<ref name="GarrabrantsEtAl2024"/>. Representative method modifications include:
 +
* Materials of construction for experimental apparatus: containers used for leaching vessels (Methods 1313A, 1316A) and column construction materials (Method 1314A) evaluated for background PFAS and PFAS uptake.
 +
* Reagents and eluant composition: eluant composition was optimized to use 1 mM CaCl2 to reduce formation of colloidal matter; Method 1313A pH adjustment now conducted with nonoxidizing HCl.
 +
* Experimental conditions: Longer equilibration times (e.g., Method 1313, 1316) may be required due to slow desorption kinetics of certain PFAS from soil and organic matrices, implementation of settling to facilitate separation of solids from eluates.
 +
* Eluate processing (all methods): Use of centrifugation in lieu of eluate filtering, sonication of bottle prior to eluate subsampling.
  
::Where:
+
==Batch Test Demonstration Studies==
:::''Q'' = Flow rate (Volume of groundwater flow per time, such as m<sup>3</sup>/yr)
+
PFAS-specific adaptations were tested in batch test demonstration studies, which included triplicate implementation of Methods 1313A and 1316A in four AFFF-impacted site soils.
:::''A'' = Cross sectional area perpendicular to groundwater flow (length<sup>2</sup>, such as m<sup>2</sup>)
 
:::''V<sub>D</sub>'' = “Darcy Velocity”; describes groundwater flow as the volume of flow through a unit of cross-sectional area (units of length per time, such as ft/yr)
 
:::''K'' = Hydraulic Conductivity (sometimes called “permeability”) (length per time)
 
:::''ΔH'' = Difference in hydraulic head between two lateral points (length)
 
:::''ΔL'' = Distance between two lateral points (length)
 
  
[https://en.wikipedia.org/wiki/Hydraulic_conductivity Hydraulic conductivity] (Table 1 and Figure 2) is a measure of how easily groundwater flows through a porous medium, or alternatively, how much energy it takes to force water through a porous medium. For example, fine sand has smaller pores with more frictional resistance to flow, and therefore lower hydraulic conductivity compared to coarse sand, which has larger pores with less resistance to flow (Figure 2).  
+
[[File: GuelfoFig1.png | thumb | 500 px | Figure 1: Figure 1. a) Overview of LEAF Method 1313A and b) PFHxS leaching as a function of pH<ref name="GarrabrantsEtAl2024"/>. Definitions: lower limit of quantification (LLOQ) and method detection limit (MDL)]]
 +
'''Draft Method 1313A''' was used to evaluate pH-dependent leaching in PFAS-contaminated soils in parallel batch extractions where each set of batch reactors is prepared and equilibrated at different pH (Figure 1).  Short-chain PFAS (≤6 fluorinated carbons) generally showed little to no variation in leaching across the tested pH range of 2-13 (e.g., [[Wikipedia: Perfluorohexanesulfonic acid | PFHxS]], Figure 1), whereas long-chain PFAS exhibited increased leaching at higher pH<ref name="GarrabrantsEtAl2024"/>. This trend is consistent with previous findings showing that soil-water partitioning coefficients (''K<sub>d</sub>'') decrease as pH increases (e.g., Higgins and Luthy 2006)<ref name="HigginsLuthy2006">Higgins, C.P., Luthy, R.G., 2006. Sorption of Perfluorinated Surfactants on Sediments. Environmental Science and Technology, 40(23), pp. 7251–7256. [https://doi.org/10.1021/es061000n doi: 10.1021/es061000n]</ref>. The most pronounced pH effects were observed for perfluoroalkyl sulfonamides (FASAs) such as [[Wikipedia: Perfluorooctanesulfonamide | perfluorooctane sulfonamide (FOSA)]], which transition from neutral to anionic forms within the circumneutral pH range (~pH 6). The anionic form has a lower ''K<sub>d</sub>'' and results in higher leaching concentrations<ref name="GarrabrantsEtAl2024"/><ref name="NguyenEtAl2020">Nguyen, T.M.H., Bräunig, J., Thompson, K., Thompson, J., Kabiri, S., Navarro, D.A., Kookana, R.S., Grimison, C., Barnes, C.M., Higgins, C.P., McLaughlin, M.J., Mueller, J.F., 2020. Influences of Chemical Properties, Soil Properties, and Solution pH on Soil–Water Partitioning Coefficients of Per- and Polyfluoroalkyl Substances (PFASs). Environmental Science and Technology, 54(24), pp. 15883–15892. [https://doi.org/10.1021/acs.est.0c05705 doi: 10.1021/acs.est.0c05705]&nbsp; [[Media: NguyenEtAl2020.pdf | Open Access Article]]</ref>. For many site management scenarios where pH is circumneutral, variations in anionic PFAS leaching are expected to be small over the relevant pH range. In such cases, when testing time and costs are primary considerations, Method 1313A may be a lower priority relative to evaluating leaching as a function of L/S (Method 1316A, Method 1314A).  Different considerations may be needed where FASAs or PFAS with multiple, ionizable functional groups (i.e., [[Wikipedia: Zwitterion | zwitterions]]) are of concern.
  
[[File:AdvectionFig2.PNG|400px|thumbnail|left|Figure 2. Hydraulic conductivity of selected rocks<ref>Heath, R.C., 1983. Basic ground-water hydrology, U.S. Geological Survey Water-Supply Paper 2220, 86 pgs. [//www.enviro.wiki/images/c/c4/Heath-1983-Basic_groundwater_hydrology_water_supply_paper.pdf Report pdf]</ref>.]]
+
[[File: GuelfoFig2.png | thumb | 500 px | Figure 2: a) Overview of LEAF Method 1316A and b) PFHxS leaching as a function of L/S ratio evaluated in parallel batch leaching vessels<ref name="GarrabrantsEtAl2024"/>]]
Darcy’s Law was first described by Henry Darcy (1856)<ref>Brown, G.O., 2002. Henry Darcy and the making of a law. Water Resources Research, 38(7), p. 1106. [https://doi.org/10.1029/2001wr000727 DOI: 10.1029/2001WR000727] [//www.enviro.wiki/images/4/40/Darcy2002.pdf Report.pdf]</ref> in a report regarding a water supply system he designed for the city of Dijon, France. Based on his experiments, he concluded that the amount of water flowing through a closed tube of sand (dark grey box in Figure 3) depends on (a) the change in the hydraulic head between the inlet and outlet of the tube, and (b) the hydraulic conductivity of the sand in the tube. Groundwater flows rapidly in the case of higher pressure (ΔH) or more permeable materials such as gravel or coarse sand, but flows slowly when the pressure difference is lower or the material is less permeable, such as fine sand or silt.
+
'''Draft Method 1316A''' was used to evaluate L/S-dependent leaching of PFAS in impacted soils using parallel batch extractions where each set of batch reactors is prepared and equilibrated at a different L/S. (Figure 2). Methods 1314A and 1316A are similar in intent as they both evaluate leaching as a function of L/S; however, the experimental approach differs. Method 1314A uses a flow-through column configuration (Figure 3; discussed further below).  Method 1314A may better simulate field conditions, but Method 1316A is simpler and less costly to implement.  Trends in Method 1314A and 1316A are expected to be qualitatively similar but leaching concentrations are expected to exhibit differences. Despite this, leaching studies comparing Methods 1314A and 1316A for inorganics showed that cumulative release results were within one order of magnitude<ref>Lopez Meza, S., Garrabrants, A.C., van der Sloot, H., Kosson, D.S., 2008. Comparison of the Release of Constituents from Granular Materials under Batch and Column Testing. Waste Management, 28(10), pp. 1853–1867. [https://doi.org/10.1016/j.wasman.2007.11.009 doi: 10.1016/j.wasman.2007.11.009]</ref>.  
  
[[File:Newell-Article 1-Fig3..JPG|500px|thumbnail|right|Figure 3. Conceptual explanation of Darcy’s Law based on Darcy’s experiment (Adapted from course notes developed by Dr. R.J. Mitchell, Western Washington University).]]
+
Method 1316A and Method 1314A may also provide different insights into transport mechanisms. Because Method 1316A is performed using equilibrated batch reactors at varying L/S, results can be used to develop equilibrium desorption isotherms and calculate desorption coefficients (e.g., ''K<sub><small>d,desorption</small></sub>''). Studies have shown that ''K<sub><small>d,desorption</small></sub>'' values for PFAS may be greater than ''K<sub><small>d</small></sub>'', an effect often attributed to desorption hysteresis<ref>Schaefer, C.E., Nguyen, D., Christie, E., Shea, S., Higgins, C.P., Field, J., 2022. Desorption Isotherms for Poly- and Perfluoroalkyl Substances in Soil Collected from an Aqueous Film-Forming Foam Source Area. Journal of Environmental Engineering, 148(1), Article 04021074. [https://doi.org/10.1061/(ASCE)EE.1943-7870.0001952 doi: 10.1061/(ASCE)EE.1943-7870.0001952]</ref>. Consequently Method 1316A provides a straightforward method to estimate site-specific desorption parameters. Although sorption parameters can also be inferred from column (Method 1314A) data, interpretation is often complicated by nonequilibrium processes.  Conversely, the column data can be valuable for quantifying those additional mechanisms providing transport parameters that can describe rate-limited transport (e.g., fraction of non-equilibrium sorption sites and sorption rates) and other dynamic behavior.
Since&nbsp;Darcy’s&nbsp;time,&nbsp;Darcy’s Law has been extended to develop a useful variation of Darcy's formula that is used to to calculate the actual velocity that the groundwater is moving in units such as meters traveled per year. This quantity is called “interstitial velocity” or “seepage velocity” and is calculated by dividing the Darcy Velocity (flow per unit area) by the actual open pore area where groundwater is flowing, the “effective porosity”&nbsp;(Table 1):
 
  
[[File:Newell-Article 1-Equation 2r.jpg|400px]]
+
As anticipated, trends in PFAS leaching obtained during the Method 1316A and Method 1314A demonstrations were qualitatively similar. These are further discussed below.
  
:Where:
+
==Column Test Demonstration Studies==
::''V<sub>S</sub>'' = “interstitial velocity” or “seepage velocity” (units of length per time, such as m/sec)<br />
+
[[File: GuelfoFig3.png | thumb | 500 px | Figure 3. a) Overview of LEAF Method 1314A and b) PFHxS leaching as a function of ∑(L/S) evaluated in saturated up-flow column tests<ref name="GarrabrantsEtAl2024"/>. Note that acrylic here simply refers to the column material of construction used in this round of Method 1314 testing.]]
::''V<sub>D</sub>'' = “Darcy Velocity”; describes groundwater flow as the volume of flow per unit area per time (also units of length per time)<br />
+
PFAS-specific adaptations were tested in column test demonstration studies, which included triplicate implementation of Method 1314A in three AFFF-impacted site soils.
::''n<sub>e</sub>'' = Effective porosity - fraction of cross section available for groundwater flow (unitless)
 
  
Effective porosity is smaller than total porosity. The difference is that total porosity includes some dead-end pores that do not support groundwater. Typical values for total and effective porosity are&nbsp;shown&nbsp;in&nbsp;Table&nbsp;1.
+
'''Draft Method 1314A''' was implemented in saturated, up-flow columns to evaluate leaching of PFAS as a function of cumulative L/S (∑(L/S)); Figure 3; total volume of water that has passed through the column divided by the soil mass in the column). As noted, the intent of Method 1314a and 1316a is similar, and in both tests, similar qualitative results were observed. For example, short-chain PFAS exhibited high initial concentrations that decreased rapidly. However, in Method 1314a, these rapid drops in short-chain PFAS tended to occur by ∑(L/S) ≈  2 (e.g., Site 1 and 3 soils, Figure 3) whereas in some cases, such as for PFHxS, Method 1316A produced slightly flatter elution curves than Method 1314A (Figures 2b and 3b). Long-chain PFAS generally displayed flatter elution profiles than short-chain PFAS across both methods. These trends are consistent with chain length dependent sorption documented in the literature<ref name="HigginsLuthy2006"/><ref name="NguyenEtAl2020"/><ref>Guelfo, J.L., Higgins, C.P., 2013. Subsurface Transport Potential of Perfluoroalkyl Acids at Aqueous Film-Forming Foam (AFFF)-Impacted Sites. Environmental Science and Technology, 47(9), pp. 4164–4171. [https://doi.org/10.1021/es3048043 doi: 10.1021/es3048043]</ref>. Although column modeling is beyond the scope of this article, prior studies have shown that saturated transport can be influenced by rate-limited desorption, particularly for long-chain PFAS<ref>Doria-Manzur, A., Gray, E.P., Streets, S.S., Guelfo, J.L., 2025. Per- and Polyfluoroalkyl Substances (PFAS) Transport from Biosolids-Amended Soils: An Experimental and Numerical Approach. Water Research, 288(Part B), Article 124674. [https://doi.org/10.1016/j.watres.2025.124674 doi: 10.1016/j.watres.2025.124674]&nbsp; [[Media: Doria-ManzurEtAl2026.pdf | Open Access Article]]</ref><ref>Guelfo, J.L., Wunsch, A., McCray, J., Stults, J.F., Higgins, C.P., 2020. Subsurface Transport Potential of Perfluoroalkyl Acids (PFAAs): Column Experiments and Modeling. Journal of Contaminant Hydrology, 233, Article 103661. [https://doi.org/10.1016/j.jconhyd.2020.103661 doi: 10.1016/j.jconhyd.2020.103661]&nbsp; [[Media: GuelfoEtAl2020.pdf | Open Access Manuscript]]</ref>. As noted, data from Methods 1316A and 1314A can support estimation of transport parameters representing equilibrium and nonequilibrium behavior, respectively.
  
[[File:Newell-Article 1-Fig4.JPG|500px|thumbnail|left|Figure 4. Difference between Darcy Velocity (also called Specific Discharge) and Seepage Velocity (also called Interstitial Velocity).]]
+
==LEAF Screening Evaluations==
 +
[[File: GuelfoFig4.png | thumb | 500 px | Figure 4. Example screening assessment for perfluorooctane sulfonate (PFOS) using total content and data from Methods 1313A and 1314A.  Figure format adapted from Garrabrants ''et al''. 2021<ref>Garrabrants, A.C., Kosson, D.S., Brown, K.G., Fagnant, D.P., Helms, G., Thorneloe, S.A., 2021. Methodology for Scenario-Based Assessments and Demonstration of Treatment Effectiveness Using the Leaching Environmental Assessment Framework (LEAF). Journal of Hazardous Materials, 406, Article 124635. [https://doi.org/10.1016/j.jhazmat.2020.124635 doi: 10.1016/j.jhazmat.2020.124635]&nbsp; [[Media: GarrabrantsEtAl2021.pdf | Open Access Manuscript]]</ref>.]]
 +
LEAF provides a standardized, robust approach for evaluating PFAS release from impacted granular materials under a range of environmental conditions. The tests are complementary, capture a range of conditions, and vary in ease of implementation. This provides the flexibility for users to select the test or test combinations that best suit their project objectives, timeline, and budget. A common use of LEAF data is in screening level assessments.  These are stepwise assessments that establish increasingly refined maximum leaching concentrations, ''C<sub><small>leach,max</small></sub>'' (Figure 4), which can then be compared to regulatory limits such as maximum contaminant levels, when available. For example, a stepwise screening assessment might include:
 +
#Assume the maximum leaching concentration is represented by the total content (total initial mass of contaminant present) leaching into the first L/S.
 +
#Assume only the available content leaches into the first L/S where available content is the maximum mass released over pH 2-13 measured using Method 1313a. For many PFAS, total content is equal to available content meaning that all of the PFAS mass is available for leaching.
 +
#Assume the leaching concentration at natural pH (measured in Method 1313A at natural pH or Method 1316A at L/S of 10) is maximum leaching concentration adjusted to the first L/S.
 +
#Consider the maximum leaching concentration over the L/S range (Method 1314A or Method 1316A) and the upper estimate of leaching, or ''C<sub><small>leach,max</small></sub>'', is the concentration from either Step 3 or Step 4, whichever is greater.
  
==Darcy Velocity and Seepage Velocity==
+
Screening assessments may be sufficient to meet project goals, but when additional refinements of leaching estimates are needed, site-specific data (e.g., infiltration) can be combined with test data and computational approaches (e.g., fate and transport models) for more site-specific estimates of leaching. Example scenarios where LEAF may be used to evaluate PFAS-impacted solids include 1) estimating PFAS release from AFFF-impacted soils, 2) estimating PFAS release from biosolids-amended soils at land application sites, 3) providing transport parameters to model PFAS transport from the source zone to the saturated zone, and 4) evaluating PFAS release from treatment residuals such as soils or sediments treated by soil washing or thermal approaches.  
In&nbsp;groundwater&nbsp;calculations, Darcy Velocity and Seepage Velocity are used for different purposes. For any calculation where the actual flow rate in units of volume per time (such as liters per day or gallons per minute) is involved, then the original Darcy Equation should be used (calculate ''V<sub>D</sub>'', Equation 1) without using effective porosity. When calculating solute travel time however, the seepage velocity calculation (''V<sub>S</sub>'', Equation 2) must be used and an estimate of effective porosity is required. Figure 4 illustrates the differences between Darcy Velocity and&nbsp;Seepage&nbsp;Velocity.
 
  
==Mobile Porosity==
+
==Summary and Ongoing Research==
{| class="wikitable" style="float:right; margin-left:10px; text-align:center;"
+
The LEAF framework offers a reliable, replicable approach to evaluating PFAS release from solids. With recent adaptations for PFAS-specific considerations, LEAF methods provide valuable tools for regulators and practitioners in managing PFAS-contaminated materials and assessing long-term environmental risksHowever, there are key areas of ongoing research, including:
|+Table 2.  Mobile porosity estimates from 15 tracer tests<ref name="Payne2008">Payne, F.C., Quinnan, J.A. and Potter, S.T., 2008. Remediation Hydraulics. CRC Press. ISBN 9780849372490  Available from: [https://www.routledge.com/Remediation-Hydraulics/Payne-Quinnan-Potter/p/book/9780849372490 CRC Press]</ref>
+
*An interlab validation of Methods 1313A, 1314A, and 1316A in coordination with the EPA
|-
+
*Optimization and demonstration of Method 1315A for use with PFAS-impacted solids
!Aquifer Material
+
*Evaluation of an unsaturated Method 1314A protocol to assess the need for and ability of LEAF testing to capture air-water interfacial partitioning of PFAS
!Mobile Porosity<br /><small>(volume fraction)</small>
+
*Application of the total oxidizable precursor (TOP) assay for evaluating the maximum additional PFAS leaching that may occur as a result of polyfluoroalkyl precursor transformation
|-
+
*Comparison of LEAF testing data to previously collected field-scale enhanced flushing data collected from the same site
|Poorly sorted sand and gravel||0.085
 
|-
 
|Poorly sorted sand and gravel||0.04-0.07
 
|-
 
|Poorly sorted sand and gravel||0.09
 
|-
 
|Fractured sandstone||0.001-0.007
 
|-
 
|Alluvial formation||0.102
 
|-
 
|Glacial outwash||0.145
 
|-
 
|Weathered mudstone regolith||0.07-0.10
 
|-
 
|Alluvial formation||0.07
 
|-
 
|Alluvial formation||0.07
 
|-
 
|Silty sand||0.05
 
|-
 
|Fractured sandstone||0.0008-0.001
 
|-
 
|Alluvium, sand and gravel||0.017
 
|-
 
|Alluvium, poorly sorted sand and gravel||0.003-0.017
 
|-
 
|Alluvium, sand and gravel||0.11-0.18
 
|-
 
|Siltstone, sandstone, mudstone||0.01-0.05
 
|}
 
 
 
Payne&nbsp;et&nbsp;al.&nbsp;(2008)&nbsp;reported the results from multiple short-term tracer tests conducted to aid the design of amendment injection systems<ref name="Payne2008" />. In these tests, the dissolved solutes were observed to migrate more rapidly through the aquifer than could be explained with typically reported values of ''n<sub>e</sub>''. They concluded that the heterogeneity of unconsolidated formations results in a relatively small area of an aquifer cross section carrying most of the water, and therefore solutes migrate more rapidly than expected. Based on these results, they recommend that a quantity called “mobile porosity” should be used in place of ''n<sub>e</sub>'' in equation 2 for calculating solute transport velocities. Based on 15 different tracer tests, typical values of mobile porosity range from 0.02 to 0.10 (Table 2).   
 
 
 
A data mining analysis of 43 sites in California by Kulkarni et al. (2020) showed that on average 90% of the groundwater flow occurred in about 30% of cross sectional area perpendicular to groundwater flow.  These data provided “moderate (but not confirmatory) support for the&nbsp;mobile&nbsp;porosity&nbsp;concept.”<ref name="Kulkarni2020">Kulkarni, P.R., Godwin, W.R., Long, J.A., Newell, R.C., Newell, C.J., 2020. How much heterogeneity? Flow versus area from a big data perspective. Remediation 30(2), pp. 15-23. [https://doi.org/10.1002/rem.21639 DOI: 10.1002/rem.21639]  [//www.enviro.wiki/images/9/9b/Kulkarni2020.pdf Report.pdf]</ref>
 
 
 
==Advection-Dispersion-Reaction Equation==
 
The transport of dissolved solutes in groundwater is often modeled using the Advection-Dispersion-Reaction (ADR) equation. As shown below (Equation 3), the ADR equation describes the change in dissolved solute concentration (''C'') over time (''t'') where groundwater flow is oriented along the ''x'' direction.
 
 
 
{|
 
| || [[File:AdvectionEq3r.PNG|center|635px]]
 
|-
 
| Where: ||
 
|-
 
|
 
:''R''
 
| is the linear, equilibrium retardation factor (see [[Sorption of Organic Contaminants]]),
 
|-
 
|
 
:''D<sub>x</sub>, D<sub>y</sub>, and D<sub>z</sub>''&nbsp;&nbsp;
 
| are hydrodynamic dispersion coefficients in the ''x, y'' and ''z'' directions (L<sup>2</sup>/T),
 
|-
 
|
 
:''v''
 
| is the advective transport or seepage velocity in the ''x'' direction (L/T), and
 
|-
 
|
 
:''λ''
 
| is an effective first order decay rate due to combined biotic and abiotic processes (1/T).
 
|}
 
 
 
The term on the left side of the equation is the rate of mass change per unit volume.  On the right side are terms representing the solute flux due to dispersion in the ''x, y'', and ''z'' directions, the advective flux in the ''x'' direction, and the first order decay due to biotic and abiotic processes. Dispersion coefficients (''D<sub>x,y,z</sub>'') are commonly estimated using the following relationships (Equation 4):
 
 
 
{|
 
| || [[File:AdvectionEq4.PNG|center|360px]]
 
|-
 
| Where: ||
 
|-
 
|
 
:''D<sub>m</sub>''
 
| is the molecular diffusion coefficient (L<sup>2</sup>/T), and
 
|-
 
|
 
:''&alpha;<sub>L</sub>, &alpha;<sub>T</sub>'', and ''&alpha;<sub>V</sub>''&nbsp;&nbsp;
 
| are the longitudinal, transverse and vertical dispersivities (L), respectively.
 
|}
 
 
 
==ADR Applications==
 
[[File:AdvectionFig5.png | thumb | right | 350px | Figure 5. Curves of concentration versus distance (a) and concentration versus time (b) generated by solving the ADR equation for a continuous source of a non-reactive tracer with ''R'' = 1, λ = 0, ''v'' = 5 m/yr, and ''D<sub>x</sub>'' = 10 m<sup>2</sup>/yr.]]
 
The ADR equation can be solved to find the spatial and temporal distribution of solutes using a variety of analytical and numerical approaches.  The design tools [https://www.epa.gov/water-research/bioscreen-natural-attenuation-decision-support-system BIOSCREEN]<ref name="Newell1996">Newell, C.J., McLeod, R.K. and Gonzales, J.R., 1996. BIOSCREEN: Natural Attenuation Decision Support System - User's Manual, Version 1.3. US Environmental Protection Agency, EPA/600/R-96/087. [https://www.enviro.wiki/index.php?title=File:Newell-1996-Bioscreen_Natural_Attenuation_Decision_Support_System.pdf Report.pdf]  [https://www.epa.gov/water-research/bioscreen-natural-attenuation-decision-support-system BIOSCREEN website]</ref>, [https://www.epa.gov/water-research/biochlor-natural-attenuation-decision-support-system BIOCHLOR]<ref name="Aziz2000">Aziz, C.E., Newell, C.J., Gonzales, J.R., Haas, P.E., Clement, T.P. and Sun, Y., 2000. BIOCHLOR Natural Attenuation Decision Support System. User’s Manual - Version 1.0. US Environmental Protection Agency, EPA/600/R-00/008.  [https://www.enviro.wiki/index.php?title=File:Aziz-2000-BIOCHLOR-Natural_Attenuation_Dec_Support.pdf Report.pdf]  [https://www.epa.gov/water-research/biochlor-natural-attenuation-decision-support-system BIOCHLOR website]</ref>, and [https://www.epa.gov/water-research/remediation-evaluation-model-chlorinated-solvents-remchlor REMChlor]<ref name="Falta2007">Falta, R.W., Stacy, M.B., Ahsanuzzaman, A.N.M., Wang, M. and Earle, R.C., 2007. REMChlor Remediation Evaluation Model for Chlorinated Solvents - User’s Manual, Version 1.0. US Environmental Protection Agency. Center for Subsurface Modeling Support, Ada, OK.  [[Media:REMChlorUserManual.pdf | Report.pdf]]  [https://www.epa.gov/water-research/remediation-evaluation-model-chlorinated-solvents-remchlor REMChlor website]</ref> (see also [[REMChlor - MD]]) employ an analytical solution of the ADR equation.  [https://www.usgs.gov/software/mt3d-usgs-groundwater-solute-transport-simulator-modflow MT3DMS]<ref name="Zheng1999">Zheng, C. and Wang, P.P., 1999. MT3DMS: A Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems; Documentation and User’s Guide. Contract Report SERDP-99-1 U.S. Army Engineer Research and Development Center, Vicksburg, MS. [[Media:Mt3dmanual.pdf | Report.pdf]]  [https://www.usgs.gov/software/mt3d-usgs-groundwater-solute-transport-simulator-modflow MT3DMS website]</ref> uses a numerical method to solve the ADR equation using the head distribution generated by the groundwater flow model MODFLOW<ref name="McDonald1988">McDonald, M.G. and Harbaugh, A.W., 1988. A Modular Three-dimensional Finite-difference Ground-water Flow Model, Techniques of Water-Resources Investigations, Book 6, Modeling Techniques. U.S. Geological Survey, 586 pages. [https://doi.org/10.3133/twri06A1  DOI: 10.3133/twri06A1]  [[Media: McDonald1988.pdf | Report.pdf]]  Free MODFLOW download from: [https://www.usgs.gov/mission-areas/water-resources/science/modflow-and-related-programs?qt-science_center_objects=0#qt-science_center_objects USGS]</ref>.
 
 
 
Figures 5a and 5b were generated using a numerical solution of the ADR equation for a non-reactive tracer (''R'' = 1; λ = 0) with ''v'' = 5 m/yr and ''D<sub>x</sub>'' = 10 m<sup>2</sup>/yr.  Figure 5a shows the predicted change in concentration of the tracer, chloride, versus distance downgradient from the continuous contaminant source at different times (0, 1, 2, and 4 years).  Figure 5b shows the change in concentration versus time (commonly referred to as the breakthrough curve or BTC) at different downgradient distances (10, 20, 30 and 40 m).  At 2 years, the mid-point of the concentration versus distance curve (Figure 5a) is located 10 m downgradient (x = 5 m/yr * 2 yr).  At 20 m downgradient, the mid-point of the concentration versus time curves (Figure 5b) occurs at 4 years (t = 20 m / 5 m/yr).
 
  
 +
Additionally, there are key areas for consideration in future research:
 +
*Collection of paired field-laboratory data under ambient conditions to further validate the applicability of LEAF assessments for estimation of field-relevant PFAS leaching and mobility
 +
*Consideration of biotransformation in modeling and interpretation of LEAF data, as the state of the science regarding biotransformation of polyfluoroalkyl substances to terminal perfluoroalkyl acids advances
  
<br clear="left" />
+
==Other LEAF Resources==
 +
There are numerous  resources describing the development of the LEAF leaching methods for inorganics. They are not specific to PFAS, but are still valuable resources focused on LEAF implementation, applications, and management of LEAF data. They include:
 +
*[https://www.vanderbilt.edu/leaching/leach-xs-lite/ Leach XS Lite] - a tool for LEAF data management and visualization; free to download after registering for a free license key
 +
*[https://www.epa.gov/hw-sw846/how-guide-leaching-environmental-assessment-framework LEAF “How-To” Guide] - guidance on LEAF background, implementation, test result interpretation; includes case studies
 +
*[https://www.epa.gov/hw-sw846/leaching-environmental-assessment-framework-leaf-methods-and-guidance USEPA LEAF Methods and Guidance] homepage
 +
<br clear="right"/>
  
 
==References==
 
==References==
 
 
<references />
 
<references />
  
 
==See Also==
 
==See Also==
 
*[http://iwmi.dhigroup.com/solute_transport/advection.html International Water Management Institute Animations]
 
*[http://www2.nau.edu/~doetqp-p/courses/env303a/lec32/lec32.htm NAU Lecture Notes on Advective Transport]
 
*[https://www.youtube.com/watch?v=00btLB6u6DY MIT Open CourseWare Solute Transport: Advection with Dispersion Video]
 
*[https://www.youtube.com/watch?v=AtJyKiA1vcY Physical Groundwater Model Video]
 
*[https://www.coursera.org/learn/natural-attenuation-of-groundwater-contaminants/lecture/UzS8q/groundwater-flow-review Online Lecture Course - Groundwater Flow]
 

Latest revision as of 18:39, 13 August 2026

PFAS Leaching Characterization with the Leaching Environmental Assessment Framework (LEAF)

Aqueous film-forming foams (AFFFs) are a major source of per- and poly-fluoroalkyl substances (PFAS) impacts in soil and groundwater. Standardized tools are needed to rapidly assess the potential for retention, leaching, and transport of PFAS from the source zone to downgradient regions, so that this information can be applied towards critical facets of site management such as prioritizing PFAS-impacted sites for further investigation and remediation. Existing standard leaching methods were developed prior to concerns regarding PFAS. Therefore, studies are needed to ensure that leaching methods are compatible for use with PFAS and that resulting data are representative of the risk of PFAS leaching at impacted sites.

Related Article(s):

Contributors: Dr. Jennifer L. Guelfo, Dr. David Kosson, Dr. Andy Garrabrants, Ms. Fangfei Liu, Mr. Darlington Yawson, Dr. Md. Isreq Real

Key Resources:

  • Development of Leaching Tests for Materials Containing SVOCs and PFAS, EPA 600/R-23/382[1]

Introduction to LEAF

The U.S. Environmental Protection Agency (EPA) Leaching Environmental Assessment Framework (LEAF) is a suite of standardized test methods for evaluating contaminant release from solids under environmentally relevant conditions (Table 1). The four leaching methods within LEAF were originally validated for inorganic constituents[2][3] and included as standard leaching methods EPA 1313 – 1316 within Update V of SW-846[4]. To address the need for standardized tests to evaluate PFAS leaching and mobility, LEAF methods have been optimized and demonstrated for use with PFAS (Methods 1313A-1316A)[1]. This article will focus on Methods 1313A, 1314A, and 1316A. Demonstration of Method 1315A for compacted granular materials (including concrete and asphalt) is ongoing.

Table 1. EPA SW-846 methods that comprise the LEAF framework
Method Description
1313 Liquid-solid partitioning as a function of extract pH using a parallel batch extraction (i.e., equilibrium) procedure (Figure 1)
1314 Liquid-solid partitioning as a function of liquid-solid ratio (L/S) for constituents in solid materials using an up-flow percolation column procedure (Figure 3)
1315 Mass transfer rates of constituents in monolithic or compacted granular materials using a semi-dynamic tank leaching procedure
1316 Liquid-solid partitioning as a function of L/S using a parallel batch extraction (i.e., equilibrium) procedure (Figure 2)
Note: Text shown in bold indicates primary condition evaluated in each method.

Method Development for PFAS

Complete details of the development of LEAF Methods 1313A, 1314A, and 1316A for use with PFAS are available in Garrabrants et al., 2024[1]. Representative method modifications include:

  • Materials of construction for experimental apparatus: containers used for leaching vessels (Methods 1313A, 1316A) and column construction materials (Method 1314A) evaluated for background PFAS and PFAS uptake.
  • Reagents and eluant composition: eluant composition was optimized to use 1 mM CaCl2 to reduce formation of colloidal matter; Method 1313A pH adjustment now conducted with nonoxidizing HCl.
  • Experimental conditions: Longer equilibration times (e.g., Method 1313, 1316) may be required due to slow desorption kinetics of certain PFAS from soil and organic matrices, implementation of settling to facilitate separation of solids from eluates.
  • Eluate processing (all methods): Use of centrifugation in lieu of eluate filtering, sonication of bottle prior to eluate subsampling.

Batch Test Demonstration Studies

PFAS-specific adaptations were tested in batch test demonstration studies, which included triplicate implementation of Methods 1313A and 1316A in four AFFF-impacted site soils.

Figure 1: Figure 1. a) Overview of LEAF Method 1313A and b) PFHxS leaching as a function of pH[1]. Definitions: lower limit of quantification (LLOQ) and method detection limit (MDL)

Draft Method 1313A was used to evaluate pH-dependent leaching in PFAS-contaminated soils in parallel batch extractions where each set of batch reactors is prepared and equilibrated at different pH (Figure 1). Short-chain PFAS (≤6 fluorinated carbons) generally showed little to no variation in leaching across the tested pH range of 2-13 (e.g., PFHxS, Figure 1), whereas long-chain PFAS exhibited increased leaching at higher pH[1]. This trend is consistent with previous findings showing that soil-water partitioning coefficients (Kd) decrease as pH increases (e.g., Higgins and Luthy 2006)[5]. The most pronounced pH effects were observed for perfluoroalkyl sulfonamides (FASAs) such as perfluorooctane sulfonamide (FOSA), which transition from neutral to anionic forms within the circumneutral pH range (~pH 6). The anionic form has a lower Kd and results in higher leaching concentrations[1][6]. For many site management scenarios where pH is circumneutral, variations in anionic PFAS leaching are expected to be small over the relevant pH range. In such cases, when testing time and costs are primary considerations, Method 1313A may be a lower priority relative to evaluating leaching as a function of L/S (Method 1316A, Method 1314A). Different considerations may be needed where FASAs or PFAS with multiple, ionizable functional groups (i.e., zwitterions) are of concern.

Figure 2: a) Overview of LEAF Method 1316A and b) PFHxS leaching as a function of L/S ratio evaluated in parallel batch leaching vessels[1]

Draft Method 1316A was used to evaluate L/S-dependent leaching of PFAS in impacted soils using parallel batch extractions where each set of batch reactors is prepared and equilibrated at a different L/S. (Figure 2). Methods 1314A and 1316A are similar in intent as they both evaluate leaching as a function of L/S; however, the experimental approach differs. Method 1314A uses a flow-through column configuration (Figure 3; discussed further below). Method 1314A may better simulate field conditions, but Method 1316A is simpler and less costly to implement. Trends in Method 1314A and 1316A are expected to be qualitatively similar but leaching concentrations are expected to exhibit differences. Despite this, leaching studies comparing Methods 1314A and 1316A for inorganics showed that cumulative release results were within one order of magnitude[7].

Method 1316A and Method 1314A may also provide different insights into transport mechanisms. Because Method 1316A is performed using equilibrated batch reactors at varying L/S, results can be used to develop equilibrium desorption isotherms and calculate desorption coefficients (e.g., Kd,desorption). Studies have shown that Kd,desorption values for PFAS may be greater than Kd, an effect often attributed to desorption hysteresis[8]. Consequently Method 1316A provides a straightforward method to estimate site-specific desorption parameters. Although sorption parameters can also be inferred from column (Method 1314A) data, interpretation is often complicated by nonequilibrium processes. Conversely, the column data can be valuable for quantifying those additional mechanisms providing transport parameters that can describe rate-limited transport (e.g., fraction of non-equilibrium sorption sites and sorption rates) and other dynamic behavior.

As anticipated, trends in PFAS leaching obtained during the Method 1316A and Method 1314A demonstrations were qualitatively similar. These are further discussed below.

Column Test Demonstration Studies

Figure 3. a) Overview of LEAF Method 1314A and b) PFHxS leaching as a function of ∑(L/S) evaluated in saturated up-flow column tests[1]. Note that acrylic here simply refers to the column material of construction used in this round of Method 1314 testing.

PFAS-specific adaptations were tested in column test demonstration studies, which included triplicate implementation of Method 1314A in three AFFF-impacted site soils.

Draft Method 1314A was implemented in saturated, up-flow columns to evaluate leaching of PFAS as a function of cumulative L/S (∑(L/S)); Figure 3; total volume of water that has passed through the column divided by the soil mass in the column). As noted, the intent of Method 1314a and 1316a is similar, and in both tests, similar qualitative results were observed. For example, short-chain PFAS exhibited high initial concentrations that decreased rapidly. However, in Method 1314a, these rapid drops in short-chain PFAS tended to occur by ∑(L/S) ≈ 2 (e.g., Site 1 and 3 soils, Figure 3) whereas in some cases, such as for PFHxS, Method 1316A produced slightly flatter elution curves than Method 1314A (Figures 2b and 3b). Long-chain PFAS generally displayed flatter elution profiles than short-chain PFAS across both methods. These trends are consistent with chain length dependent sorption documented in the literature[5][6][9]. Although column modeling is beyond the scope of this article, prior studies have shown that saturated transport can be influenced by rate-limited desorption, particularly for long-chain PFAS[10][11]. As noted, data from Methods 1316A and 1314A can support estimation of transport parameters representing equilibrium and nonequilibrium behavior, respectively.

LEAF Screening Evaluations

Figure 4. Example screening assessment for perfluorooctane sulfonate (PFOS) using total content and data from Methods 1313A and 1314A. Figure format adapted from Garrabrants et al. 2021[12].

LEAF provides a standardized, robust approach for evaluating PFAS release from impacted granular materials under a range of environmental conditions. The tests are complementary, capture a range of conditions, and vary in ease of implementation. This provides the flexibility for users to select the test or test combinations that best suit their project objectives, timeline, and budget. A common use of LEAF data is in screening level assessments. These are stepwise assessments that establish increasingly refined maximum leaching concentrations, Cleach,max (Figure 4), which can then be compared to regulatory limits such as maximum contaminant levels, when available. For example, a stepwise screening assessment might include:

  1. Assume the maximum leaching concentration is represented by the total content (total initial mass of contaminant present) leaching into the first L/S.
  2. Assume only the available content leaches into the first L/S where available content is the maximum mass released over pH 2-13 measured using Method 1313a. For many PFAS, total content is equal to available content meaning that all of the PFAS mass is available for leaching.
  3. Assume the leaching concentration at natural pH (measured in Method 1313A at natural pH or Method 1316A at L/S of 10) is maximum leaching concentration adjusted to the first L/S.
  4. Consider the maximum leaching concentration over the L/S range (Method 1314A or Method 1316A) and the upper estimate of leaching, or Cleach,max, is the concentration from either Step 3 or Step 4, whichever is greater.

Screening assessments may be sufficient to meet project goals, but when additional refinements of leaching estimates are needed, site-specific data (e.g., infiltration) can be combined with test data and computational approaches (e.g., fate and transport models) for more site-specific estimates of leaching. Example scenarios where LEAF may be used to evaluate PFAS-impacted solids include 1) estimating PFAS release from AFFF-impacted soils, 2) estimating PFAS release from biosolids-amended soils at land application sites, 3) providing transport parameters to model PFAS transport from the source zone to the saturated zone, and 4) evaluating PFAS release from treatment residuals such as soils or sediments treated by soil washing or thermal approaches.

Summary and Ongoing Research

The LEAF framework offers a reliable, replicable approach to evaluating PFAS release from solids. With recent adaptations for PFAS-specific considerations, LEAF methods provide valuable tools for regulators and practitioners in managing PFAS-contaminated materials and assessing long-term environmental risks. However, there are key areas of ongoing research, including:

  • An interlab validation of Methods 1313A, 1314A, and 1316A in coordination with the EPA
  • Optimization and demonstration of Method 1315A for use with PFAS-impacted solids
  • Evaluation of an unsaturated Method 1314A protocol to assess the need for and ability of LEAF testing to capture air-water interfacial partitioning of PFAS
  • Application of the total oxidizable precursor (TOP) assay for evaluating the maximum additional PFAS leaching that may occur as a result of polyfluoroalkyl precursor transformation
  • Comparison of LEAF testing data to previously collected field-scale enhanced flushing data collected from the same site

Additionally, there are key areas for consideration in future research:

  • Collection of paired field-laboratory data under ambient conditions to further validate the applicability of LEAF assessments for estimation of field-relevant PFAS leaching and mobility
  • Consideration of biotransformation in modeling and interpretation of LEAF data, as the state of the science regarding biotransformation of polyfluoroalkyl substances to terminal perfluoroalkyl acids advances

Other LEAF Resources

There are numerous resources describing the development of the LEAF leaching methods for inorganics. They are not specific to PFAS, but are still valuable resources focused on LEAF implementation, applications, and management of LEAF data. They include:


References

  1. ^ 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 Garrabrants, A.C., Liu, F., Warne, R., DeLapp, R., Brown, L., Rubin, Z., Yawson, D., Kosson, D.S., Guelfo, J.L., Real, M.I., van der Sloot, H.A., Touati, A., Thorneloe, S., 2024. Development of Leaching Tests for Materials Containing SVOCs and PFAS, EPA 600/R-23/382, USEPA, Washington, D.C. Free Download EPA 600/R-23/382
  2. ^ Garrabrants, A.C., Kosson, D.S., Stefanski, L., DeLapp, R., Seignette, P.F.A.B., van der Sloot, H.A., Kariher, P., Baldwin, M., 2012. Interlaboratory Validation of the Leaching Environmental Assessment Framework (LEAF) Method 1313 and Method 1316, EPA/600/R-12/623, U.S. Environmental Protection Agency, Air Pollution and Control Division. Free Download EPA 600/R-12/623
  3. ^ Garrabrants, A.C., Kosson, D.S., DeLapp, R., Kariher, P., Seignette, P.F.A.B., van der Sloot, H.A., Stefanski, L., Baldwin, M., 2012. Interlaboratory Validation of the Leaching Environmental Assessment Framework (LEAF) Method 1314 and Method 1315, EPA/600/R-12/624, U.S. Environmental Protection Agency, Air Pollution and Control Division. Free Download EPA 600/R-12/624
  4. ^ USEPA, 2026. Hazardous Waste Test Methods / SW-846. USEPA SW-846 website
  5. ^ 5.0 5.1 Higgins, C.P., Luthy, R.G., 2006. Sorption of Perfluorinated Surfactants on Sediments. Environmental Science and Technology, 40(23), pp. 7251–7256. doi: 10.1021/es061000n
  6. ^ 6.0 6.1 Nguyen, T.M.H., Bräunig, J., Thompson, K., Thompson, J., Kabiri, S., Navarro, D.A., Kookana, R.S., Grimison, C., Barnes, C.M., Higgins, C.P., McLaughlin, M.J., Mueller, J.F., 2020. Influences of Chemical Properties, Soil Properties, and Solution pH on Soil–Water Partitioning Coefficients of Per- and Polyfluoroalkyl Substances (PFASs). Environmental Science and Technology, 54(24), pp. 15883–15892. doi: 10.1021/acs.est.0c05705  Open Access Article
  7. ^ Lopez Meza, S., Garrabrants, A.C., van der Sloot, H., Kosson, D.S., 2008. Comparison of the Release of Constituents from Granular Materials under Batch and Column Testing. Waste Management, 28(10), pp. 1853–1867. doi: 10.1016/j.wasman.2007.11.009
  8. ^ Schaefer, C.E., Nguyen, D., Christie, E., Shea, S., Higgins, C.P., Field, J., 2022. Desorption Isotherms for Poly- and Perfluoroalkyl Substances in Soil Collected from an Aqueous Film-Forming Foam Source Area. Journal of Environmental Engineering, 148(1), Article 04021074. doi: 10.1061/(ASCE)EE.1943-7870.0001952
  9. ^ Guelfo, J.L., Higgins, C.P., 2013. Subsurface Transport Potential of Perfluoroalkyl Acids at Aqueous Film-Forming Foam (AFFF)-Impacted Sites. Environmental Science and Technology, 47(9), pp. 4164–4171. doi: 10.1021/es3048043
  10. ^ Doria-Manzur, A., Gray, E.P., Streets, S.S., Guelfo, J.L., 2025. Per- and Polyfluoroalkyl Substances (PFAS) Transport from Biosolids-Amended Soils: An Experimental and Numerical Approach. Water Research, 288(Part B), Article 124674. doi: 10.1016/j.watres.2025.124674  Open Access Article
  11. ^ Guelfo, J.L., Wunsch, A., McCray, J., Stults, J.F., Higgins, C.P., 2020. Subsurface Transport Potential of Perfluoroalkyl Acids (PFAAs): Column Experiments and Modeling. Journal of Contaminant Hydrology, 233, Article 103661. doi: 10.1016/j.jconhyd.2020.103661  Open Access Manuscript
  12. ^ Garrabrants, A.C., Kosson, D.S., Brown, K.G., Fagnant, D.P., Helms, G., Thorneloe, S.A., 2021. Methodology for Scenario-Based Assessments and Demonstration of Treatment Effectiveness Using the Leaching Environmental Assessment Framework (LEAF). Journal of Hazardous Materials, 406, Article 124635. doi: 10.1016/j.jhazmat.2020.124635  Open Access Manuscript

See Also