Subscriber access provided by UNIV AUTONOMA DE COAHUILA UADEC
Environmental Modeling
Material and Site Specific Partition Coefficients for Beneficial Use Assessments Nawaf I. Blaisi, Kyle Clavier, Justin Roessler, Jaeshik Chung, Timothy G. Townsend, Souhail R Al-Abed, and Jean-Claude J Bonzongo Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.9b01756 • Publication Date (Web): 29 Jul 2019 Downloaded from pubs.acs.org on August 7, 2019
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 33
Environmental Science & Technology
1
Material and Site Specific Partition Coefficients for Beneficial Use Assessments
2
Nawaf I. Blaisia, Kyle A. Claviera, Justin G. Roesslera, Jaeshik Chunga,b, Timothy G.
3
Townsenda*, Souhail R. Al-Abedc , Jean-Claude J. Bonzongoa
4
a Department
5
PO Box 116450 Gainesville, FL 32611 – 6450, USA
6
7
b Center
for Water Resource Cycle, Korea Institute of Science and Technology, Seoul 136791, Republic of Korea
8
9
of Environmental Engineering Sciences, University of Florida,
c
National Risk Management Research Laboratory, U.S. Environmental Protection Agency,
10
26 West Martin Luther King Drive, Cincinnati, OH 45268, USA
11
Submitted to:
12
Environmental Science and Technology
13
14
15
16 17 18
* Corresponding
author: Phone: +1 352-392-0846, Fax: 352-392-3076, email
[email protected] ACS Paragon Plus Environment
Environmental Science & Technology
20
TOC ART
21
ACS Paragon Plus Environment
Page 2 of 33
Page 3 of 33
23
Environmental Science & Technology
Abstract
24
Partition coefficients (Kd) available in the literature are often used in fate and transport
25
modeling conducted as part of beneficial use risk assessments for industrial byproducts. Since
26
element partitioning depends on soil properties as well as characteristics of the byproduct
27
leachate, site-specific Kd may lead to more accurate risk assessment. In this study,
28
contamination risk to groundwater of beneficially reused byproducts was assessed using batch
29
leaching tests on waste to energy bottom ash (WTE BA) and coal combustion fly ash (CFA).
30
Leachates were equilibrated with 8 different soils to obtain the waste-soil-specific Kd,exp for the
31
metals of interest. The Kd,exp values were used as inputs in the Industrial Waste Management
32
Evaluation Model (IWEM) to demonstrate the degree to which Kd estimates affect risk
33
assessment outcome. Measured Kd,exp for the most part fell within the large range of Kd reported
34
in the literature, but IWEM results using default Kd for some types of soils resulted in
35
overestimated risk compared to those derived from Kd,exp. Modeled concentration at the
36
receptor location was much lower for some elements for those soils with high concentrations
37
of iron and aluminum.
38 39 40 41
Keywords: coal combustion fly ash, dilution attenuation factor, leachate, partitioning
42
coefficient, waste to energy bottom ash
ACS Paragon Plus Environment
Environmental Science & Technology
43
Introduction
44
Recycling of industrial byproducts through beneficial use applications such as aggregate
45
replacement, fill material, and soil amendment plays an important role in sustainable materials
46
management.
47
resource consumption and reduced land disposal requirements, must be weighed against the
48
potential risks posed by chemicals present in these byproducts [1, 2]. In addition to assessing
49
potential for direct human exposure, a reused byproduct’s environmental risk is assessed by
50
evaluating its potential to leach contaminants into groundwater or surface water [3–5]. This
51
type of risk assessment involves estimation of concentrations of contaminants of potential
52
concern (COPCs) using leaching tests and relies on fate and transport modeling to predict
53
COPC concentrations at target receptor location or compliance point [6,7]. Fate and transport
54
models for contamination of groundwater typically require inputs that quantify subsurface
55
environment characteristics such as: infiltration rate, aquifer and vadose zone thickness, aquifer
56
pH, soil partitioning coefficient (Kd), receptor distance, and COPC concentration, among others
57
[7–10]. Each of these inputs influence modeled COPC concentrations at receptor sites and most
58
have been the subject of considerable research [11–14].
Benefits of these materials management strategies, such as reduced natural
59
The solid to liquid partitioning coefficient (Kd) between a waste leachate and soil interface,
60
a measure of the ratio of a contaminant in the solid and liquid phase at equilibrium, is one input
61
parameter used in fate and transport modeling that is extensively studied in the literature
62
[6,15,16]. While most beneficial use evaluations employ leaching tests to determine waste-
63
specific COPC concentrations, determinations of site-specific Kd’s are less frequent. Most risk
64
assessment protocols utilize reference Kd values selected from the literature or Kd values
65
determined from programs such as MINTEQA2 [8]. Relying on reference Kd that are not
66
specific to site and waste characteristics when assessing beneficial reuse risk may result in less
ACS Paragon Plus Environment
Page 4 of 33
Page 5 of 33
Environmental Science & Technology
67
accurate contamination risk profiles, especially considering the broad range of Kd reported in
68
the literature and their dependence on material-specific properties [17]. Consider, for instance,
69
a beneficial use scenario in which the default Kd value for a given element used in a fate and
70
transport model is higher than the actual Kd of the soil on site. The model would predict greater
71
attenuation of that element in the soil and thus lower concentrations at the receptor site location;
72
this would result in an under-estimation of leaching risk associated with the beneficial use
73
application and potential harm to human health or the environment. Conversely, a low default
74
Kd value would overestimate leaching risk and might inappropriately disqualify a candidate
75
beneficial use material based on perceived risk. An ideal scenario involves using a Kd value
76
that is specific to a given application and will allow for more accurate quantification of risk. It
77
is highly unlikely that default partitioning coefficients account for all of the complex
78
interactions between waste and soil with leachate, such as organic matter, trace element
79
content, and soil pH. Thus the default Kd values may be inaccurate for specific scenarios and
80
lead to one of the above-mentioned scenarios.
81
Here we demonstrate a method to determine waste- and soil-specific Kd values (herein
82
referred to as Kd,exp) and use them in fate and transport models to provide a more representative
83
prediction of risks associated with beneficial reuse projects utilizing common industrial waste
84
byproducts.
85
Leaching concentrations measured from waste to energy (WTE) bottom ash (BA) and coal
86
combustion fly ash (CFA) were input into US EPA’s Industrial Waste Management Evaluation
87
Model (IWEM) along with waste- and soil- specific Kd,exp. The Kd,exp values were determined
88
by experimentally measuring partition coefficients for eight different soil sources contacted
89
with WTE BA and CFA leachate; the different soil sources help highlight the importance of
90
major soil properties (e.g. iron and aluminum content) on the outcome of beneficial use
91
determinations.
ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 33
92
Materials and Methods
93
Soil and Waste Sample Collection. Soil samples were collected at a depth of 5-25 cm below
94
the ground surface from eight different locations in Florida, US. Soil sampling locations were
95
selected to provide a wide range of soils which would differ in both chemical composition and
96
physical characteristics, so that the impact of these parameters on the Kd,exp could be evaluated.
97
Each soil sample was air dried at room temperature (~ 24±2 oC) for three days, passed through
98
a 2-mm sieve, and refrigerated at 4 oC until used in laboratory experiments. Two grab samples
99
of industrial waste byproducts were collected from facilities in Florida, US. CFA samples were
100
collected from a coal-fired power generation unit. The samples were collected from the
101
discharge of the facility’s fabric filters and did not contain the scrubber residues associated
102
with acid gas removal. Samples of WTE BA were collected from a mass burn WTE facility.
103
This facility subjects its bottom ash to ferrous (magnet) and non-ferrous (eddy current
104
separation) metals recovery following combustion. As commonly performed prior to beneficial
105
use applications, the material was aged for three months before collection.
106
Soil
Classification
and
Elemental
Content.
The
operationally
defined
total
107
environmentally available concentrations of metals in soil samples were determined after
108
subjecting six replicates of each soil sample to acid digestion following EPA Method 3050b.
109
This strong acid digestion dissolves almost all elements, except for those in forms that are
110
generally not considered “environmentally available,” such as those bound by silicate
111
structures [18]. The percentages of sand, silt and clay in each soil sample were determined
112
based on grain size distribution in accordance with the procedures outlined in ASTM D422
113
[19]. The pH of the soil samples was determined following EPA Method 9045D [20]. Soil
114
moisture content was measured by determining the initial and oven-dried mass in accordance
115
with ASTM D2216 [21]. The soil organic matter content was estimated through loss on ignition
116
(LOI) after four hours of heating at 550 °C as described by Santisteban et al. (2004) [22].
ACS Paragon Plus Environment
Page 7 of 33
Environmental Science & Technology
117
Batch Leaching and Soil-Liquid Partitioning Tests. Samples were first size-reduced to
118
pass a 4.75 mm sieve and leachates were generated using a modified EPA method 1316 test at
119
a set liquid to solid ratio (L/S) of 10 mL-reagent water/g-waste [23]. Each waste was leached
120
in eight replicates; a 200-g dry mass of waste was used to generate a sufficient quantity of
121
leachate for the soil sorption tests.
122
To determine the Kd,exp , the leachate generated from the 1316 test was contacted with each
123
of the eight soils (each soil sorption test was conducted in triplicate). The soil sorption testing
124
was performed following ASTM D4646, a procedure designed to measure a soils affinity for
125
select constituents in leachates. This test involves equilibrating a given mass of soil sample
126
with a waste solute, such as a laboratory extraction leachate, of known trace element
127
composition in a 1:20 soil-to-solution ratio. Adsorbed solute can be calculated by measuring
128
remaining solute concentration after the sorption test, adjusted for extract volume and soil mass
129
[24]. For the purposes of this study, soil sorption tests assume trace element contribution from
130
the clean soil is negligible. It is also assumed that sorption and desorption are reversible and
131
that the leachate and soil have equilibrated during the contact time prescribed by ASTM D4646.
132
Acid Digestion and Elemental Analysis. All leachate samples were digested prior to
133
analysis using an automated hot block digestion system following EPA Method 3010A [25].
134
The acid digestates were evaluated for their major and trace element content using Inductively
135
Coupled Plasma – Atomic Emission Spectrometry (ICP-AES), in accordance with EPA
136
Method 6010B [26]. Cl- and SO42- in leachates were measured using ion chromatography (IC)
137
following EPA Method 9056A [27].
138
Waste- and Soil-specific Partitioning Coefficient Determinations. Kd,exp was calculated
139
as the difference between constituent waste leachate concentration (Cinitial) from method 1316
140
and the equilibrium aqueous concentration following the soil sorption test (Cfinal) divided by
141
Cfinal, adjusted for soil mass (Msoil) and extract volume (Vsolution), assuming the linear region in
ACS Paragon Plus Environment
Environmental Science & Technology
142
the sorption isotherm. The Kd,exp’s reported in this study are all expressed in units of L/kg (see
143
equation 1).
144 145
Equation 1. Kd,exp calculation using data generated from L/S 10 1316 extractions and the soil
146
sorption test outlined in ASTM D4646. Units for Cinitial and Cfinal are mg/L, Vsolution is in L, and
147
Msoil is kg. Final units for Kd,exp are therefore L/kg.
148 149
𝐾
𝑑,exp =
(𝐶𝑖𝑛𝑖𝑡𝑖𝑎𝑙 ― 𝐶𝑓𝑖𝑛𝑎𝑙)(𝑉𝑠𝑜𝑙𝑢𝑡𝑖𝑜𝑛) (𝐶𝑓𝑖𝑛𝑎𝑙)(𝑀𝑠𝑜𝑖𝑙)
150 151
Example of Beneficial Use Assessment. To determine the scale at which COPCs were
152
present in the WTE BA and CFA, leachate concentrations were compared to Florida, US
153
Groundwater Cleanup Target Levels (GCTLs). Risk thresholds such as GCTL are commonly
154
used in beneficial use assessments; as a screening step, leach test results are compared to
155
GCTL, but a complete risk assessment requires using leach test results as an input to a fate and
156
transport model to assess whether GCTL are expected to be exceeded at receptor sites or
157
compliance points [28]. Critical components of this assessment technique include the Kd, waste
158
leachate concentrations, soil and hydrogeologic conditions, and the risk-based thresholds. To
159
illustrate the impact that the Kd,exp can have on fate and transport modeling outcomes, two
160
beneficial use scenarios were evaluated using IWEM. IWEM has been used as a modeling tool
161
in the US at the state and federal level for beneficial use assessments, and was recently
162
employed as a part of the decision-making efforts related to the regulatory status of coal
163
combustion residuals [8]. A detailed description of the program can be found in the Supporting
164
Information (SI). The two scenarios evaluated were the use of WTE BA as a roadway sub-base
165
(Table S2) and the use of CFA as a structural fill material (Table S1). Many of the input
ACS Paragon Plus Environment
Page 8 of 33
Page 9 of 33
Environmental Science & Technology
166
parameters (infiltration rate, fill volume, roadway thickness, subsurface hydrogeology) for
167
these models were set using values developed for two previous assessments conducted by
168
Benson and Edil, and the Electric Power Research Institute [29-30].
169
Input concentrations were set at the COPC concentrations measured in the Method 1316
170
leachates for both the CFA (structural fill) and WTE BA (road base). Each of the modeling
171
evaluations were performed at three different Kd: the default Kd for each element (from IWEM)
172
and values representing the lowest and the highest calculated Kd,exp for each element from the
173
soil sorption testing (on 8 different soils). The default Kd values in IWEM employ multiple sets
174
of sorption isotherms derived from the geochemical speciation model MINTEQA2. Additional
175
information on IWEM can be found in the IWEM Model 3.1 Technical Background Document
176
[8].
177
Results and Discussion
178
Soil Characterization. Soil characterization results of the 8 samples tested, including soil
179
pH, conductivity, percentage of organic matter, total metal (Al, Fe, and Mn) content, and soil
180
classification conducted in accordance with the methodology outlined by the US Department
181
of Agriculture (USDA) [31] are presented in Table 1. Total concentrations of Al, Fe and Mn
182
are especially pertinent because (hydr)oxide forms of those metals on the sorbent surface are
183
known to be efficient scavengers of some trace elements [32–36]. Iron concentrations in soil
184
samples ranged from 167 to 5,800 mg-Fe/kg-soil. Soils 2, 4, 5 and 7 had relatively high Fe
185
contents (>1,900 mg-Fe/kg-soil), while soils 1, 3, 6, and 8 had lower Fe contents (