Material- and Site-Specific Partition Coefficients for Beneficial Use

Jul 29, 2019 - Partition coefficient (Kd) values available in the literature are often used in fate and transport modeling conducted as part of benefi...
0 downloads 0 Views 508KB Size
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 (