Hybrid Modeling Approach to Estimate Exposures of Hazardous Air

Oct 25, 2016 - *Phone: 919-541-4650; fax: 919-541- 4511; e-mail: ... Citation data is made available by participants in Crossref's Cited-by Linking se...
0 downloads 0 Views 1MB Size
Article

A hybrid modeling approach to estimate exposures of hazardous air pollutants (HAPs) for the National Air Toxics Assessment (NATA) Richard D. Scheffe, Madeleine Strum, James Thurman, Sharon B Phillips, Alison Eyth, Steve Fudge, Mark Morris, Ted Palma, and Richard Cook Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b04752 • Publication Date (Web): 25 Oct 2016 Downloaded from http://pubs.acs.org on October 26, 2016

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 free 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 accessible to all readers and 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.

Environmental Science & Technology 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 31

Environmental Science & Technology

1

A hybrid modeling approach to estimate exposures

2

of hazardous air pollutants (HAPs) for the National

3

Air Toxics Assessment (NATA)

4

Richard D. Scheffe1,*, Madeleine Strum1, Sharon B. Phillips1, James Thurman1 , Alison Eyth1,

5

Steve Fudge2, Mark Morris1, Ted Palma1, Richard Cook3

6

* Corresponding author email: [email protected]; Phone: 919-541-4650; Fax: 919-541-

7

4511

8

1

U.S. Environmental Protection Agency, Office of Air Quality Planning and Standards, Durham,

9

NC 27711, United States 3

10 11

3

EC/R Incorporated, Chapel Hill, NC 27514, United States

U.S. Environmental Protection Agency, Office of Transportation and Air Quality, Ann Arbor,

12

MI 48105, United States

13

KEYWORDS

14

Hazardous Air Pollutants

15

Hybrid Modeling

16

Air Toxics Risk

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 31

17

ABSTRACT

18

A hybrid air quality model has been developed and applied to estimate annual concentrations of

19

40 hazardous air pollutants (HAPs) across the continental United States (CONUS) to support the

20

2011 calendar year National Air Toxics Assessment (NATA).

21

transport model (CTM) with a Gaussian dispersion model, both reactive and non-reactive HAPs

22

are accommodated across local to regional spatial scales, through a multiplicative technique

23

designed to improve mass conservation relative to previous additive methods. The broad scope

24

of multiple pollutants capturing regional to local spatial scale patterns across a vast spatial

25

domain is

26

improved performance relative to the stand alone CTM and dispersion model. However, model

27

performance varies widely across pollutant categories and quantifiably definitive performance

28

assessments are hampered by a limited observation base and challenged by the multiple physical

29

and chemical attributes of HAPs. Formaldehyde and acetaldehyde are the dominant HAP

30

concentration and cancer risk drivers, characterized by strong regional signals associated with

31

naturally emitted carbonyl precursors enhanced in urban transport corridors with strong mobile

32

source sector emissions.

33

dominated source sectors creates largely similar concentration patterns across the majority of

34

HAPs. However, reactive carbonyls exhibit significantly less spatial variability relative to non-

35

reactive HAPs across the CONUS.

By combining a chemical

precedent setting within the air toxics community. The hybrid design exhibits

The multiple pollutant emission characteristics of combustion

36

ACS Paragon Plus Environment

2

Page 3 of 31

Environmental Science & Technology

37 38

INTRODUCTION

39

The 1990 Clean Air Act Amendments (CAA) list 187 hazardous air pollutants (HAPs) known

40

to cause or suspected of causing cancer as well as respiratory, neurological, reproductive and

41

other serious health effects. HAPs include a variety of volatile (e.g., formaldehyde, benzene,

42

1,3-butadiene) and semi-volatile (e.g., naphthalene, PAH congeners) organic compounds (VOCs

43

and SVOCs) and metals (e.g., arsenic, nickel, hexavalent chromium). These air toxics are

44

emitted by mobile sources (e.g., cars, trucks and construction equipment); large or major

45

stationary sources (e.g., chemical manufacturing, refineries and power plants); smaller, or area,

46

sources (e.g., gas stations and dry cleaners); and natural processes (biogenic VOC releases and

47

wildfires). Markedly different air quality management approaches are applied to air toxics

48

relative to the National Ambient Air Quality Standards (NAAQS) criteria pollutants (ozone,

49

particulate matter, nitrogen dioxide, sulfur dioxide, carbon monoxide and lead), despite similar

50

emission sources. The NAAQS set defined ambient air target concentration levels based on

51

human and environmental welfare risk assessments, and emission strategies are developed to

52

meet those targets. In contrast, air toxics emitted from stationary sources are regulated by

53

emissions targets through the National Emissions Standards for Hazardous Air Pollutants

54

(NESHAPS)1 authorized by section 112d of the CAA. Also, mobile source toxics are regulated

55

based on technological feasibility (under Section 202(l) of the CAA), and also reduced by

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 31

56

controls on hydrocarbons and particulate matter (for example, Tier 3 and 2007 heavy duty

57

vehicle standards), and the Mobile Sources Air Toxics (MSAT)2 rules. Consequently, regulatory

58

modeling to characterize air toxics concentrations is limited to individual source sector Risk and

59

Technology Reviews (RTRs) which attempt to quantify residual risk after implementation of the

60

NESHAPS Maximum Available Control Technology (MACT) rules.

61

National Air Toxics Assessment (NATA) complements RTR modeling by providing a highly

62

spatially resolved national view of air toxics concentrations and risk patterns of 178 HAPs, not

63

achievable through a limited observations network.

The non–regulatory

64 65

Here we describe the air quality modeling approach and results driving the 2011 calendar

66

year NATA which estimates HAPs inhalation risks associated with cancer incidences and

67

noncancer respiratory effects. NATA modeling provides census tract level risk values driven by

68

annual concentration estimates based on EPA’s triannually released national emissions inventory

69

(NEI). NATAs have been developed for calendar years 1996, 1999, 2002, and 2005 to provide a

70

screening of potential high risk areas across the U.S. that warrant further analysis.3

71 72

The NATA presents a comprehensive multiple-pollutant challenge to estimate HAPs with

73

a range of reactive and thermodynamic properties influencing exposures across local to regional

74

spatial scales.

75

previous urban scale studies that integrate a variety of modeling platforms and observational

76

data.4-10 Collectively, those urban scale applications expanded the range of exposures and

77

reduced exposure misclassification, relative to observation based exposure studies. This national

This national scale hybrid modeling application stands apart from several

ACS Paragon Plus Environment

4

Page 5 of 31

Environmental Science & Technology

78

scale application transforms that scope of exposures into a true multiple pollutant, multiple

79

spatial scale application over a vast geographic domain.

80 81

METHODS

82

Structure and overview. The 2011 NATA air quality model, referred to herein as the hybrid

83

method, utilizes the fine spatial scale and direct source attribution features of the AERMOD

84

dispersion model11-12 and the full treatment of chemistry and transport afforded by the

85

Community Multiscale Air Quality (CMAQ) model, version 5.02 with the Carbon Bond 05

86

(CB05) chemical mechanism13.

87

Eulerian air quality model designed to simulate the formation and fate of gaseous and particulate

88

species, including ozone, oxidant precursors, primary and secondary PM concentrations and

89

sulfur, nitrogen and mercury deposition over urban and regional spatial scales. In this

90

application, AERMOD treats all species as chemically non-reactive and excludes deposition.

91

AERMOD receptor locations are based on the centroids of populated census blocks, monitoring

92

site positions and five evenly distributed points within each 12 km horizontal CMAQ grid

93

(Figure S1) resulting in a minimum of 5 receptors to cells with over 10,000 receptors and 6.5

94

million receptors nationwide.

CMAQ is a comprehensive, three-dimensional grid-based

95 96

Equation 1 is used to calculate 2011 annual average air concentrations at receptor locations,

97

which are constrained by forcing the average of within grid cell receptor concentrations to equal

98

the CMAQ surface level grid value with AERMOD providing sub-grid scale spatial texture:

99

ACS Paragon Plus Environment

5

Environmental Science & Technology

100 101

C = AERMODREC × (CMAQPNFB /AERMODGRIDAVG) + CMAQSEC

Page 6 of 31

+

CMAQPBIOGENICS

CMAQPFIRES + (1)

102 103

Where;

104 105

C = concentration at a receptor,

106

AERMODREC = the AERMOD receptor concentration,

107

CMAQPNFB = CMAQ grid cell concentration contributed by primary emissions, excluding fires

108

and biogenics,

109

AERMODGRIDAVG = the average of all AERMOD results within a CMAQ grid, calculated

110

through surface interpolation of all AERMOD receptor locations across CMAQ grid cells to

111

eliminate concentration discontinuities,

112

CMAQSEC = the CMAQ grid cell contribution from atmospheric reactions, and

113

CMAQPFIRES and CMAQPBIOGENICS = CMAQ contribution from primary emissions of fires and

114

biogenics, respectively, which are not incorporated in AERMOD.

115 116

CMAQ separately tracks primary and secondary contributions, enabling the AERMOD

117

estimate at each receptor location to be normalized to the non-fire, non-biogenic CMAQ primary

118

contribution. By anchoring concentration averages to CMAQ, mass conservation is improved,

119

which can be violated in direct additive model combinations or when observations are directly

120

added to model results, a procedure used in previous NATA iterations. For example, the most

121

recent 2005 NATA was susceptible to duplicate counting as an added “background”

122

concentration based on ambient observations was added to AERMOD estimates. Diagnosing

ACS Paragon Plus Environment

6

Page 7 of 31

Environmental Science & Technology

123

model behavior based on paired model to measurement values can be compromised by the dual

124

use of observations and largely inconsistent model inputs driving AERMOD and CMAQ.

125

Therefore, the receptor average mass allocation constraint to CMAQ grid values imposed by

126

equation (1) is appropriate when combining results from vastly different model architectures as it

127

minimizes potential double counting. The emissions and meteorological data sets used to drive

128

CMAQ were processed further to generate AERMOD inputs consistent with CMAQ. This

129

hybrid approach, which builds on earlier area-specific applications to Philadelphia, PA6 and

130

Detroit, MI9, reflects an evolution of national scale HAPs modeling intended to optimize

131

characterization of non-reactive and reactive species across multiple spatial scales.

132 133

Meteorological data processing: The gridded meteorological data for 2011 at the 12 km

134

continental U.S. scale domain was derived from version 3.4 of the Weather Research and

135

Forecasting Model (WRF)14. The WRF model was initialized using the 12km North American

136

Model (12NAM) analysis product, initialized by all available observations15, provided by

137

National Climatic Data Center (NCDC). Where 12NAM data was unavailable, the 40km Eta

138

Data Assimilation System (EDAS) analysis (ds609.2) from the National Center for Atmospheric

139

Research (NCAR) was used. Landuse and land cover data were based on the 2006 National

140

Land Cover Database16. WRF meteorological outputs were processed using the Meteorology-

141

Chemistry Interface Processor (MCIP) package, version 4.1.3, to derive hourly specific inputs to

142

CMAQ: horizontal wind components (i.e., speed and direction), temperature, moisture, vertical

143

diffusion rates, and rainfall rates for each grid cell in each vertical layer17. CMAQ resolved the

144

vertical atmosphere with 25 layers, preserving greater resolution in the planetary boundary layer

145

(PBL). The meteorological inputs to AERMOD were based on the underlying prognostic data,

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 31

146

i.e. 12 km WRF data, that was used in CMAQ. The AERMOD meteorological inputs were

147

created using the Mesoscale Model Interface Program18 (MMIF). MMIF was used to extract

148

data at every 4th WRF grid cell, as well as over 700 WRF grid cells containing National Weather

149

Service (NWS) surface stations. Each AERMOD simulation of a source calculated the air

150

quality impacts of that source at receptors out to 50 km from the source. The closest WRF grid

151

cell that was extracted from MMIF was used for the respective AERMOD simulation. All

152

receptors received impacts from sources within 50 km and cumulative concentrations were

153

calculated by adding concentrations from all modeled sources impacting the receptor. Isakov et

154

al.6 demonstrated successful use of meteorological variables derived from prognostic modeling

155

to drive dispersion models, motivated primarily to address spatial gaps in meteorological

156

monitoring. Additionally, MMIF outputs have been shown to compare favorably well against

157

observed meteorological data when used in AERMOD19.

158

Emissions Processing:

The 2011 National Emissions Inventory (NEI) provided the root

159

emissions data for CMAQ and AERMOD20.

160

broad categories (emissions input resolution down to over 800 source classification codes is

161

retained – Text S1) with similar spatial and temporal delineation: major point sources, non-point

162

sources (excluding transportation, fires and biogenics), on-road mobile, and non-road mobile

163

(including locomotive, aviation and commercial marine vessels). Fires (combined wild,

164

prescribed and agricultural) and biogenic emissions are handled only through CMAQ. NEI data,

165

provided as specific point and aggregated county level annual estimates, are processed to hourly

166

values distributed over 12 km horizontal grids through the Sparse Matrix Operator Kernel

167

Emissions (SMOKE) modeling system21. Biogenic emissions were generated by version 3.60

168

the Biogenic Emissions Inventory System (BEIS)22. Hourly temporal allocations were developed

Emissions to AERMOD are grouped into four

ACS Paragon Plus Environment

8

Page 9 of 31

Environmental Science & Technology

169

for AERMOD, consistent with CMAQ, a departure from previous NATAs that were based on

170

annual average inputs.

171

typically were allocated spatially to population census tract resolution using a variety of

172

surrogates (e.g., land use classifications, population)19.

Non-point source, on-road and non-road emissions for AERMOD

173

Initial and Boundary Conditions (IC/BCs): The CMAQ lateral boundary and initial species

174

concentrations for benzene, formaldehyde and acetaldehyde were generated by a 2011 year GEOS-Chem

175

simulation23. Due to the scarcity of observations suitable for establishing BCs and the

176

extended calendar year simulation, zero value IC/BCs were used for the remaining hybrid

177

model HAPs.

178 179

Treatment of species. The hybrid model was applied to 40 (Table S1) high risk HAPs among

180

a total of 178 HAPs included in the 2011 NATA. The remaining 138 HAPs (Exhibit B-119) were

181

modeled just with AERMOD. Although this application focuses on HAPS, it reflects the second

182

major national scale U.S. EPA application of the multiple pollutant version of CMAQ, following

183

a national assessment of increased ethanol use associated with renewable fuels24. The

184

atmospheric chemistry treatments in chemical transport models (CTMs) such as CMAQ are

185

based on gas phase reaction processes optimized to characterize ozone, linked with a variety of

186

heterogeneous and thermodynamic processes to accommodate particulate matter (PM) formation.

187

Consequently, the inclusion of explicit HAP species in current chemical mechanisms is

188

predicated by its relative importance in ozone chemistry. Formaldehyde and acetaldehyde are

189

explicit species that generate significant amounts of peroxy radicals leading to enhanced ozone

190

production and secondary particulate matter formation, and also are high risk HAPs which

191

exemplify multiple pollutant linkages driven by atmospheric processes.

HAP species not

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 31

192

incorporated explicitly in the chemical mechanism are added as non-reactive tracers (e.g., several

193

halogenates and napthalene). Other important and reactive HAPs, including benzene, toluene,

194

xylene isomers, 1,3 butadiene and acrolein, were added to the CB05 chemical mechanism in

195

order to explicitly track their decay through reactions with time-varying radicals and oxidants.

196

Since the impact of these species on the chemistry and radical cycling in the mechanism is

197

already accounted for in CB05 lumped model species, this explicit representation consisted of

198

decay without modification of the radicals or oxidants, and in the case of 1,3-butadiene,

199

production of secondary acrolein. Future applications of the hybrid model will increase the

200

number of hybrid species as part of a continuous CMAQ development program, which will add

201

explicit treatment of napthalene, benzene, acrolein, 1,3 butadiene, xylenes, and toluene.

202

Additional non-reactive tracers to be added include carbon tetrachloride, ethylbenzene,

203

chloroprene and methyl chloride. AERMOD treats all HAPs as nonreactive and was applied to

204

the remaining HAPs not incorporated within CMAQ using previous methodology to develop the

205

complete 2011 NATA19.

206 207

Source Attribution: Estimates of the source contributions associated with primary emissions

208

were generated by the following ratio technique normalized to CMAQ concentrations for the

209

AERMOD source groups:

210 211

[CREC,J] = AERMODREC,J × [CMAQPNFB]/[AERMODGRIDAVG]

(2)

212 213

where [CREC,J] = the contribution to concentration at a receptor from source category J,

214

excluding secondary formation. This ratio approach provides an estimate of primary emission

ACS Paragon Plus Environment

10

Page 11 of 31

Environmental Science & Technology

215

contributions only. Primary emission contributions from biogenics and fires were processed

216

through CMAQ and all contributions from secondary formation processes were aggregated into

217

CMAQSEC, resulting in all within grid receptors receiving identical contributions from fires,

218

biogenics and secondary contributions.

219

Risk Calculations: Cancer and noncancer risks are estimated by relating annual exposure

220

concentration estimates to relevant health toxicity values.

Exposure concentrations were

221

generated by the Hazardous Air Pollutant Exposure Model19 (HAPEM), which modified hybrid

222

concentrations to account for movement of individuals across indoor, outdoor and commuting

223

environments. Differences between exposure and ambient concentrations were relatively minor

224

for most species. The toxicity values (Table S1) are quantitative expressions used to estimate the

225

likelihood of adverse health effects given an .estimated level and duration of exposure. Because

226

NATA is focused on long-term exposures, the toxicity values used are based on the results of

227

chronic dose-response studies when such data are available. Cancer risk is estimated from

228

inhalation exposure concentration using a cancer unit risk estimate (URE), which is the upper-

229

bound excess lifetime cancer risk estimated to result from continuous exposure to an agent at a

230

concentration of 1 µg/m3. Noncancer risk (hazard) for a pollutant is benchmarked to a Reference

231

Concentration (RfC), which is an estimate of a continuous inhalation exposure that is thought to

232

be without an appreciable risk of deleterious health effects over a lifetime. The resulting

233

noncancer metrics for each pollutant (the hazard quotient) are aggregated into an index if the

234

pollutants act by similar toxic modes of action or affect the same target organ. Only non-cancer

235

respiratory effects are reported here, which is the most significant for inhalation exposure to the

236

pollutants included in NATA.

ACS Paragon Plus Environment

11

Environmental Science & Technology

Page 12 of 31

237

Model Evaluation and Air Quality Observations: Annual average comparisons of paired

238

model estimates and observations for a subset of HAPs were developed using 2011 ambient

239

observations obtained from Phase 9 of EPA’s Air Toxics Monitoring Archive25, derived largely

240

from EPA’s Air Quality System (AQS). State and local agencies and Tribes (SLTs) are required

241

to submit data generated from EPA’s National Air Toxics Trends (NATTS) network and the

242

Urban Air Toxics Monitoring Program (UATMP) to the AQS. The UATMP provides centralized

243

laboratory analytical support for air toxics networks (Figure S2) 26-27, which are described in the

244

supporting materials (Text S2).

245 246

Results and Discussion

247

Model Evaluation

248 249

We focus on a core group of 8 HAPs (benzene, formaldehyde, 1,3 butadiene, acetaldehyde,

250

acrolein, nickel, arsenic and napthalene) among the 40 hybrid species, that aggregated represent over 90%

251

of the national HAPs associated cancer and noncancer risk and reflect the variety of chemical and

252

physical properties across the complete set of 178 HAPs. Paired average annual model to monitor

253

site comparisons are presented for the hybrid model along with stand-alone CMAQ and

254

AERMOD for several HAPs. Acrolein was excluded given the well-known and substantial

255

sampling artifact issues associated with acrolein formation inside canisters and stability of

256

calibration standards28. Inclusion of all three model results is intended to demonstrate the

257

merged attributes of the hybrid model.

258 259

Significant scatter occurs across all pollutants with hybrid model normalized mean error

260

(NME) ranging from 38 (acetaldehyde) to 62% (1,3 butadiene) for VOCs (Figure S3, Table S2).

ACS Paragon Plus Environment

12

Page 13 of 31

Environmental Science & Technology

261

Secondarily formed species, formaldehyde and acetaldehyde, have less NME, which is associated

262

with relatively smooth spatial gradients reducing spatial distribution errors. CMAQ alone and

263

the hybrid model exhibit similar NME (acetaldehyde 36%, 38%; formaldehyde 40%, 39%) and

264

NMB (acetaldehyde 17%, 21%; formaldehyde -37%, -33%) for secondarily formed VOCs. This

265

convergence of performance between CMAQ and the hybrid is expected for pollutants exhibiting

266

relatively low spatial variability. The low NMB of -6% for benzene might appear encouraging

267

from a risk perspective based on long term averaging. However, the low bias masks the

268

canceling of positive and negative errors. Because we have relatively high confidence in the

269

accuracy and precision of benzene observations relative to other HAPs, there is a sound

270

infrastructure for further diagnosing of model behavior to truly improve the performance of

271

relatively conservative species. Over 95% of benzene observations exceed minimum detection

272

limits (MDLs)27 and there is no evidence of substantial sampling artifacts. AERMOD under

273

estimates are expected for secondarily formed HAPs (e.g., -81 NMB for acetaldehyde) given the

274

exclusion of atmospheric chemistry. Error and bias statistics between the hybrid model and

275

CMAQ were considerably more aligned relative to either model compared to AERMOD. The

276

hybrid model exhibited improved NMB relative to CMAQ for non-secondarily formed VOCs -

277

benzene (CMAQ:hybrid -27%: -6%) and 1,3 butadiene (CMAQ:hybrid -49%:-20%). Benzene

278

appeared to exhibit over-predictions in the Northeast while no clear pattern of regional biases

279

were exhibited for carbonyls, an analysis hampered by the limited number of monitoring stations

280

(Figure S4). Very minor differences in performance statistics across the three modeling

281

approaches were exhibited for non-reactive species (napthalene, arsenic and nickel).

282

ACS Paragon Plus Environment

13

Environmental Science & Technology

283

Page 14 of 31

Formaldehyde under-predictions of nearly 30% are consistent with other modeling

284

studies29-31. Biogenically generated isoprene is the principal VOC precursor for formaldehyde.

285

BEIS generates less isoprene relative to the Model of Emissions of Gases and Aerosols from

286

Nature (MEGAN)32; however, MEGAN typically exhibits a strong positive bias for isoprene

287

leading to over estimates of formaldehyde30,33. Biogenic emissions characterization has been an

288

active research area for three decades, largely in support of ozone and particulate matter

289

modeling and, looking forward, perhaps reinvigorated by formaldehyde and acetaldehyde being

290

significant air toxics risk drivers.

291 292

Very high NMEs and NMBs for toluene and metals and order of magnitude and greater

293

underpredictions of halogenated VOCs (Table S2) are consistent with the generally poor

294

performance exhibited by non-secondarily formed species. Poor performance illustrates the

295

difficulty of characterizing HAPs that rely on a voluntary emissions reporting program lacking

296

the rigor associated with criteria pollutant emissions, compounded by species that push the

297

detection limits of current analysis methods with concentrations often below MDLs27.

298 299

This evaluation is broad in the scope of pollutant and geographic coverage relative to

300

previous NATAs and related air toxics model evaluations, which collectively lack any

301

standardized evaluation approach. For example, an analysis34 for the 1996 NATA in Baltimore,

302

MD demonstrated good comparisons based on the ratio of means, but their focus was on

303

demonstrating discrepancies between indoor and outdoor exposures. Lupo and Symanski35

304

reported ratios of model to monitor means ranging from 0.5 to 2.0 for 48% of paired

305

comparisons based on the earlier 1996 and 1999 NATAs throughout Texas urban areas, and

ACS Paragon Plus Environment

14

Page 15 of 31

Environmental Science & Technology

306

noted that concordance analysis suggested poorer model performance relative to a means

307

comparison. Logue et al.36 produced a rigorous area-specific evaluation of the 2005 NATA in

308

Pittsburgh, PA with a relatively rich observation base enabling comparison of nearly 50 HAPs.

309

Their findings illustrated a wide range of error across different HAPs, but concentrations for the

310

highest risk HAPs typically were within a factor of 2 of the observations. George et al.37

311

concluded that the NATA 2002 simulation exhibited good performance for benzene based

312

largely on a means evaluation for the Detroit metropolitan area using 2004-2007 Detroit

313

Exposure and Aerosol Research Study (DEARS) observations. Kimbrough et al.38 reported

314

NATA 2005 under predictions of formaldehyde, acetaldehyde and benzene, but not for 1,3

315

butadiene using 2008-2009 observations in Las Vegas. Vennam et al.39 applied CMAQ with 4

316

km resolution to an airport in Rhode Island for the 2005 calendar year and reported NME from

317

36-70% for eight VOC HAPs and napthalene, with modest error improvement (5-20%) relative

318

to 36 km grid resolution. Stroud et al.40 applied the Unified Regional Air-quality Modelling

319

System (AURAMS – nested 45 and 22.5 km grid cells) for 2006 over Canada for six VOCs and

320

reported NMB values of -15% and 26% for formaldehyde and acetaldehyde, respectively, and

321

generally better performance for carbonyls relative to primary VOC species - findings

322

directionally consistent with the hybrid model.

323 324

These model evaluation results demonstrate expected behavior of the hybrid model and

325

its two modeling components, reinforcing the original model design objectives to capture

326

concentration patterns of multiple HAPs across a variety of spatial scales. This diversity of

327

HAPs and spatial scales combined with a limited observation network challenges model

328

evaluation efforts. Carbonyls are compromised by a limited observations network with

ACS Paragon Plus Environment

15

Environmental Science & Technology

Page 16 of 31

329

sampling and analysis artifacts and evolving model process formulations that have not been

330

designed to explicitly address the intermediate reaction sets, of great significance to HAPs,

331

leading to ozone and secondary particulate matter formation. However, the limited spatial

332

variability of carbonyls, relative to non-reactive HAPs, results from atmospheric dynamics

333

somewhat analogous to those influencing ozone, perhaps setting a hypothetical reference frame

334

for model performance, contingent on improved observations and atmospheric process

335

characterizations. The relatively low bias in benzene comparisons may be a reasonable

336

performance goal for non-reactive HAPs accompanied by a reliable observation base.

337

Nevertheless, the high NME (50%) for benzene is a concern for such a high volume, ubiquitous

338

pollutant that benefits from a nationally consistent approach for generating transportation related

339

emissions. Arguably, the hybrid model has greater skill in capturing formaldehyde patterns as

340

much of the NME (39%) is associated with under-predictions. The poor performance of metals

341

and halogenated VOCs (Table S2) reflects the relative small investments in HAP emissions

342

inventories and ambient monitoring networks, limited by a lack of regulatory drivers. Unlike

343

ozone and PM precursors, HAP emissions reporting for State agencies is voluntary, resulting in

344

inconsistent methodologies. Relatedly, the robustness of HAP observations also is.

345 346

As HAP modeling has transitioned from a dispersion model centered framework to

347

greater reliance on CTMs, the evaluation methods also reflect those practices associated with

348

CTMs. There are no accepted benchmarks for HAP model evaluation, in contrast to ozone

349

model performance where error and bias for paired (monitoring site and hourly resolution)

350

comparisons typically are within 35 and 15%, respectively41-42. Such criteria are not realistic for

351

HAPs which are more analogous to ozone precursors, which also lack accepted performance

ACS Paragon Plus Environment

16

Page 17 of 31

Environmental Science & Technology

352

criteria. In contrast, a common assertion of “a factor of 2 agreement seems reasonable” often

353

accompanies air toxics modeling studies. In addition, aggregation to annual average estimates

354

associated with long term chronic health effects dampens errors that would arise from paired

355

temporal comparisons, but also compromises diagnostic power to improve model processes and

356

inputs.

357 358

Concentration Patterns of Selected HAPs

359 360

Pronounced national and regional scale patterns emerge for several high concentration

361

gases.

Acetaldehyde and formaldehyde, generated largely through secondary formation

362

processes, exhibit the highest annual average concentrations ranging from 0.5

363

aggregated across NOAA climate regions43 for cancer related HAPs (Table S3; Figures 1 and

364

S5).

365

concentrations of the remaining cancer benchmarked HAPs by one or more orders of magnitude.

366

Concentrations of toluene and xylene isomers, both lacking a cancer benchmark, are comparable

367

to carbonyl levels.

to 2 µg/m3

Aggregated benzene concentrations approach 1 µg/m3, exceeding the average

ACS Paragon Plus Environment

17

Environmental Science & Technology

Page 18 of 31

368 369 370 371

Figure 1. Annual average 2011 concentration estimates (µg/m3) of reactive HAPs generated by

372

the hybrid model for formaldehyde, acetaldehyde, 1,3 butadiene and acrolein, clockwise from

373

upper left panel, respectively.

374

All of these HAPs exhibit elevated levels in populated areas, indicative of similar source

375

signatures influenced strongly by onroad vehicles, as well as numerous combustion based

376

sources. However, the elevated southeast U.S. acetaldehyde and formaldehyde signals (Figures

377

1- 2) illustrate the combined influences of a rich source of biogenic precursors (isoprene),

378

atmospheric chemistry and transport. The regional nature of these HAPS governed largely by

379

chemistry and transport provides a smoothing of spatial gradients with relatively low spatial

ACS Paragon Plus Environment

18

Page 19 of 31

Environmental Science & Technology

380

coefficients of variation (Cv’s of 0.31 and 0.30 for formaldehyde and acetaldehyde), in contrast

381

to metals (Table S3) which have stronger local scale signals associated with point sources. Yu

382

and Stuart10 reported similar diffuse patterns of secondarily formed HAPs relative to primary

383

species using a hybrid modeling approach

384

estimates for the Tampa, FL area. Benzene is widely distributed across all regions and exhibits

385

less spatial variability compared to metals, although the influence of direct emissions results in

386

greater variability compared to secondarily formed species. 1,3 butadiene undergoes decay, but

387

is not formed secondarily, therefore exhibiting spatial variability (Cv = 0.81) similar to

388

conservative species benzene (Cv = 0.61) and naphthalene (Cv = 0.82).

389

reported similar relative patterns of spatial variability based on near roadway measurements of

390

benzene, 1,3 butadiene, acetaldehyde and formaldehyde.

that combined CMAQ and dispersion model

Kimbrough et al.38

391 392

Figure 2. Normalized concentrations, calculated as the ratio of the regional to national average,

393

of key HAPs by NOAA climate region43 to illustrate pollutant variability across regions. NOAA

394

has defined nine climatically consistent regions across the CONUS.

ACS Paragon Plus Environment

19

Environmental Science & Technology

Page 20 of 31

395

Acrolein is formed secondarily, but is handled simplistically in CMAQ v5.02 through 1,3

396

butadiene decay. Elevated acrolein in northern CA and a widespread region of acrolein in the

397

southeast (Figure 1) is associated with high agricultural fire emissions which are introduced into

398

the lowest CMAQ vertical layer. Emissions from wildfires and prescribed fires are allocated to

399

multiple vertical layers based on plume rise calculations, which increases dilution and transport,

400

perhaps explaining the enhanced effect of northern CA acrolein relative to larger. Because all

401

three carbonyls also are products of industrial, residential and transportation combustion sources,

402

their concentration patterns reflect an underlying structure based on anthropogenic source

403

distributions with regional enhancements associated with natural sources and fires.

404 405

Although metal HAPs exhibit greater spatial variability, their broad (Figures 1, 2 and S5)

406

and fine spatial scale (Figure 3) concentration patterns resemble non-reactive VOCs reflecting

407

the multiple pollutant attributes of ubiquitous combustion processes. Noted exceptions such as

408

elevated Ar and Ni associated with gold mining in Carlin, NV and manufacturing in southeast

409

Idaho (Figure S5), respectively, emerge from scattered metals extraction and processing

410

facilities.

411 412

Figure 3. Annual average concentration distributions of arsenic (left) and benzene(right) in the

413

New York city area illustrating similar distribution patterns (with different concentration ranges)

ACS Paragon Plus Environment

20

Page 21 of 31

Environmental Science & Technology

414

across these two HAPs, and the spatial resolution of the hybrid receptors at the census tract level

415

along with the CMAQ 12 km grid lines.

416 417 418

Risk and Source Attribution

419 420

Pollutant Risk. Formaldehyde, acetaldehyde and benzene represent 77% of total cancer

421

risk nationally, with over half the national risk associated with formaldehyde (Figure 4).

422

Formaldehyde dominates HAPs risk as it is the highest risk driver in over 99% of all U.S. census

423

tracts (Figure S6). Although not a hybrid model HAP, carbon tetrachloride (CCl4) is carried

424

along in the risk reports given the significant 8% contribution. CCl4 concentrations underlying

425

risk estimates are based on observations, justified by spatially homogeneous measurements,

426

centered about 0.6 µg/m3, reflecting a 30 year residence time and elimination of most CCl4

427

emission sources. Naphthalene and 1,3 butadiene also are significant national cancer risk

428

drivers.

ACS Paragon Plus Environment

21

Environmental Science & Technology

Page 22 of 31

429 430

Figure 4. Summary of relative cancer (left) and non-cancer (right) risk results with source

431

attribution (bottom). Note that all source categories reflect primary emission contributions only

432

and clearly understate the contributions from source categories, such as onroad and biogenetics,

433

contributing to secondarily formed HAPs.

434

fraction of the specified HAP to the aggregate risk of all HAPs. CCl4 (top left) is a non-hybrid

435

HAP that is delineated from the non-hybrid category because of its relatively high contribution to

436

risk. CCl4 will be treated as a hybrid HAP in future NATAs. “Other” represents collections of

437

hybrid HAPs which individually contribute less than 1% and 0.1 % to national cancer (top left)

438

and non-cancer risk (top right), respectively, and source groupings contributing less than 0.3%

439

to cancer risk (bottom left).

Pollutant risk percentages are calculated as the

440

ACS Paragon Plus Environment

22

Page 23 of 31

Environmental Science & Technology

441 442 443

Although formaldehyde dominates HAPs risk, other HAPs pose significant cancer risks

444

from a local scale perspective. Chloroprene, for example, emitted by a single facility in southeast

445

Louisiana generates the highest census tract risk (~800 in a million) in the nation. Of a total 178

446

HAPs considered in NATA, only 13 HAPs are high risk cancer drivers for any particular census

447

tract (Figure S6). Aside from benzene, formaldehyde and napthalene, these HAPs generally are

448

associated with localized point source facilities. In contrast, HAPs such as chloroprene and coke

449

oven emissions are modeled just with AERMOD as non-hybrid HAPs as they are limited to

450

specific point source facilities.

451 452

Noncancer risk also is driven by carbonyls with acrolein the dominant driver for 70% of

453

the risk across the U.S. (Figure 4). As noted above, acrolein will be treated explicitly in future

454

chemical mechanisms employed by the hybrid model. In addition, theses high non-cancer risk

455

results should promote research into and deployment of reliable measurement technology

456

enabling ground truthing of this important HAP. Diesel PM, based on direct particulate

457

emissions from diesel engines, accounts for nearly 5% of the noncancer related risk. Diesel

458

exhaust is a known carcinogen that currently is not included in EPA’s NATA cancer risk results.

459 460

Source attribution. The hybrid model provides limited capability to characterize

461

emissions sector contributions to concentration and risk as each AERMOD simulation is linked

462

to a specific source group and CMAQ separately tracks secondary production of carbonyls

463

(formaldehyde, acetaldehyde and acrolein) from primary contributions. The combination of high

ACS Paragon Plus Environment

23

Environmental Science & Technology

Page 24 of 31

464

concentrations, relatively high toxicity and carbonyl production through atmospheric

465

transformations leads to 47% of the cancer risk interpreted as “secondary” (Figure 4).

466

Anthropogenic contributors to secondary production are spread across many combustion based

467

sectors that emit VOC and inorganic (oxides of nitrogen and carbon monoxide) precursors for

468

carbonyl formation. Consequently, the 18% contribution of onroad sources noted in Figure 4 is

469

associated with primary emissions only, thus understating actual source attribution to sectors that

470

emit organic and inorganic carbonyl precursors.

471 472

Emissions strategy implications

473 474

Further source attribution delineation on secondary carbonyl production can be inferred

475

through instrumented CTM techniques such as direct decoupled method (DDM) and adjoint

476

methods, integrated source apportionment tracking, as well as through brute force sector zero-out

477

simulations; techniques that have been widely applied for ozone and particulate matter.

478

Dunker44 and Dunker et al.45 applied a path-integral method imbedded in the Comprehensive Air

479

Quality Model with Extensions (CAMx) to demonstrate that roughly 50 to 75% of formaldehyde

480

would naturally be derived in urban and rural locations, respectively44-45. The anthropogenic

481

increment in rural and urban locations was limited by NOx and VOC emissions, respectively,

482

partly reflecting the strong primary formaldehyde emissions in urban locations. Carlton and

483

Baker30 estimated that nearly 50% of biogenically derived secondary organic aerosol (SOA) is

484

controllable based on a series of CMAQ zero-out simulations, and reasoned that the strict

485

classification of a biogenic or anthropogenic contribution to SOA is insufficient. The same

486

reasoning is applicable to carbonyls in which the root material may be of biogenic (or fire based)

ACS Paragon Plus Environment

24

Page 25 of 31

Environmental Science & Technology

487

origin. Clearly there are analogies to ozone where “background” ozone has increased in relative

488

importance as the ozone NAAQS has tightened over time. However, two important differences

489

exist: (1) trans-continental transport of carbonyls, and their precursors, are probably insignificant

490

and (2) the long term chronic exposures associated with HAPs elevates the importance of a

491

stable pool of HAPs relative to acute based criteria pollutant exposures where peak

492

concentrations often drive risk.

493 494

For many years air toxics was viewed as a local scale or “hot spot” exposure issue.

495

However, this assessment demonstrates a dominant regional air toxics risk component. The

496

majority of HAPs are combustion derived products emitted collectively through various source

497

sectors, illustrated through the similar concentration patterns across pollutant groups (VOCs,

498

SVOCs and metals; Figures 1, 3 and S5). Consequently, HAPs fit well within a multiple

499

pollutant air management framework where the primary focus on ozone and fine particulate

500

matter abatement strategies likely impart directionally positive impacts on reducing HAPs

501

emissions. Nevertheless, the relevance of local scale exposures remains, as the occurrence of

502

anomalous high exposures to individual pollutants is augmented by the numerous HAPs of

503

concern.

504 505

The hybrid model is unique in its ability to characterize gradients of air pollutants with a

506

wide range of chemical and physical attributes across regional to local scales, reflecting progress

507

in developing a multiple pollutant air quality modeling platform. The findings reported here

508

indicating the importance of carbonyls as national risk drivers challenge perceptions that air

509

toxics is largely a local scale exposure phenomenon. Overall exposure characterization, and

ACS Paragon Plus Environment

25

Environmental Science & Technology

Page 26 of 31

510

consequent management, is not markedly different than that associated with strong regional

511

nature of ozone and fine particulate matter. Fortunately, the infrastructure developed to address

512

those NAAQS is well suited to handle the most dominant HAPs risk drivers.

513 514

Limitations and Recommendations

515 516

Limitations of the current hybrid approach clearly in addressing the multiple pollutant,

517

multiple spatial scale challenges are driven by basic data concerns. The scarcity of monitoring

518

sites, with observations challengings MDLs, combined with an evolving and voluntary emissions

519

system, challenges diagnosing the most efficient course for improvement.

520

capturing regional level exposures of high risk secondarily formed HAPs is encouraging.

521

Despite this success, added attention to characterizing carbonyls is warranted as they drive

522

national level cancer risk estimates.

523

modeling may reside in addressing local scale characterizations of less ubiquitous HAPs that

524

pose significant risks in local areas. The high spatial resolution afforded by AERMOD reveals

525

relatively minor differences in the spatial distribution of concentrations across pollutant groups

526

(Figures 1 and 3), suggesting potential efficiencies limiting the use of dispersion modeling to key

527

locations with unique source configurations. Because of these local scale concerns, a national

528

level modeling assessment at best can indicate areas requiring further analysis. Progress can be

529

catalyzed by (1) conducting focused urban and local scale field campaigns to evaluate emissions,

530

(2) requiring SLTs to report HAP emissions, (3) supporting satellite and ground based

531

technologies to improve the carbonyl observational data bases and (4) expanding the integration

The success in

Nevertheless, the most challenging aspect of HAPs

ACS Paragon Plus Environment

26

Page 27 of 31

Environmental Science & Technology

532

and coordination of HAPs emissions and modeling with criteria pollutants to optimize

533

infrastructure leveraging.

534

Acknowledgments.

535

The authors greatly appreciate the effort and insights from the anonymous reviewers that lead

536

to a greatly improved manuscript.

537

management and resource support, Karen Wesson for support in the early model development

538

stage, Chris Misenis for meteorological data processing and Halil Cakir for software support.

539 540

ASSOCIATED CONTENT

541

Supporting Information Available.

542

Figures S1- S6 and Tables S1-S3 provide additional information on model and risk results and

543

model evaluation statistics. Descriptions of monitoring methods are included. This material is

544

available free of charge via the Internet at http://pubs.acs.org.

545

AUTHOR INFORMATION

546

Corresponding Author

547

We also thank Tyler Fox and Marc Houyoux for

*E-mail: [email protected]; Phone: 919-541-4650; Fax: 919-541- 4511

548

Author Contributions

549

The manuscript was written through contributions of all authors. All authors have given approval

550

to the final version of the manuscript.

551

ACS Paragon Plus Environment

27

Environmental Science & Technology

Page 28 of 31

552

References

553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596

(1)

U.S. EPA. 2015a. National Emission Standards for Hazardous Air Pollutants, 40 CFR, Part 61. http://www.epa.gov/ttn/atw/mactfnlalph.html

(2)

U.S. EPA. 2007. Control of Hazardous Substances from Mobile Sources, 40 CFR, Vol. 7, No. 37.

(3)

U.S. EPA Website: http://www.epa.gov/national-air-toxics-assessment

(4)

Cook, R.; Isakov, V.; Touma, J.; Benjey, W.; Thurman, J.; Kinnee, E.; Ensley, D. Resolving Local Scale Emissions for Modeling Air Quality near Roadways; J. Air & Waste Manage. Assoc. 2008, 58, 451-461.

(5)

Dionisio, K.L., Isakov, V.; Baxter, L.K.; Sarnat, J.A.; Sarnat, S.E.; Burke, J.; Rosenbaum, A.; Graham, S.E.; Cook, R.; Mulholland J.; Ozkaynak, H. Development and evaluation of alternative approaches for exposure assessment of multiple air pollutants in Atlanta, Georgia, Journal of Exposure Science and Environmental Epidemiology. 2013, 23, 581–592.

(6)

Isakov, V., Irwin, J.S.; Ching, J. Using CMAQ for Exposure Modeling and Characterizing the Subgrid Variability for Exposure Estimates, J. Appl. Meteorol. Climatol. 2007, 46, 1354-1371.

(7)

Ozkaynak, H., Baxter, L.K.; Burke, J. Air pollution exposure prediction approaches used in air pollution epidemiology studies, Journal of Exposure Science and Environmental Epidemiology. 2013, 23, 566–572.

(8)

Stein, A.F.; Isakov, V.; Godowitch, J.; Draxler, R.R. A Hybrid Approach to Resolve Pollutant Concentrations in an Urban Area; Atmos. Environ. 2007, 41, 9410-9426.

(9)

Wesson, K., Fann, N.; Morris, M.; Fox, T.; Hubbell, B. A multi–pollutant, risk–based approach to air quality management: Case study for Detroit. Atmospheric Pollution Research. 2010, 1, 296‐304.

(10) Yu, H.; A.L. Stuart. Exposure and inequality for select urban air pollutants in the Tampa Bay area. Science of the Total Environment. 2016, 551–552, 474–483. (11) Cimorelli, A.J., Perry, S.G., Venkatram, A., Weil, J.C., Paine, R.J., Wilson, R.B., Lee, R.F., Peters, W.D., Brode, R.W. AERMOD: a dispersion model for industrial source applications. Part I: general model formulation and boundary layer characterization. Journal of Applied Meteorology, 2005, 44, 682–693. (12)

U.S. EPA. User's Guide for the AMS/EPA Regulatory Model – AERMOD. 2004, EPA454/B-03-001. U.S. Environmental Protection Agency, Research Triangle Park, NC 27711.

ACS Paragon Plus Environment

28

Page 29 of 31

597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641

Environmental Science & Technology

(13)

Byun D; Schere K. Review of the governing equations, computational algorithms, and other components of the models-3 community multiscale air quality (CMAQ) modeling system. Appl Mech Rev, 2006, 59:51–77

(14)

Skamarock, W.C., Klemp, J.B., Dudhia, J., Gill, D.O., Barker, D.M., Duda, M.G., Huang, X., Wang, W., Powers, J.G., 2008. A Description of the Advanced Research WRF Version 3

(15)

National Weather Service, Environmental Monitoring Center, http://www.emc.ncep.noaa.gov/index.php?branch=NAM

(16)

Multi-Resolution Land Characteriztics Consortium (MRLC), National Land Cover Data Base (NLCD), http://www.mrlc.gov/nlcd06_data.php

(17)

Otte T.L.; Pleim, J.E. The Meteorology-Chemistry Interface Processor (MCIP) for the CMAQ modeling system: updates through v3.4.1. Geoscientific Model Development 3, 2010, 243-256.

(18)

Brashers, B.; Emery, C. The Mesoscale Model Interface Program (MMIF), ENVIRON International Corporation, Air Sciences Group – Technical Report, 2015. 773 San Marin Drive Novato, CA 94998.

(19)

U.S. EPA, 2015a, Technical Support Document for EPA’s National-Scale Air Toxics Assessment, http://www3.epa.gov/airtoxics/natamain.

(20)

U.S. EPA, 2014a, Profile of the 2011 National missions Inventory, http://www.epa.gov/ttn/chief/net/lite_finalversion_ver10.pdf

(21)

Houyoux, M. R.; Vukovich, J. M.; Coats, C. J. Jr.; Wheeler, N. M.; Kasibhatla, P. S. Emission inventory development and processing for the Seasonal Model for Regional Air Quality (SMRAQ) project. J. Geophys. Res. 2000, 105, D7, pp 9079-9090.

(22)

Bash, J.O., Baker, K.R.; Beaver, M.R. Evaluation of improved land u se and canopy representation in BEIS v3.61 with biogenic VOC measurements in California, Geosci. Model Dev. Discuss. 2015, 8, 8117–8154.

(23)

Yantosca, B., Sulprizio, M.; Yannetti, M.; Lundgren, L.; Xu, J. GEOS-Chem v10-01 Online User’s Guide, Atmospheric Chemistry Modeling Group, School of Engineering and Applied Sciences, Harvard University, 2015, 29 Oxford Street, Cambridge, MA 02.

(24)

Cook, R., Phillips, S.; Houyoux, M.; Dolwick, P.; Mason, R.; Yanca, C.; Zawacki, M.; Davidson, K.; Michaels, H.; Harvey, C.; Somers, J.; Luecken, D. Air quality impacts of increased use of ethanol under the United States’ Energy Independence and Security Act, Atmospheric Environment. 2011, 45 7714-7724.

ACS Paragon Plus Environment

29

Environmental Science & Technology

642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687

Page 30 of 31

(25)

U.S. EPA website: http://www.epa.gov/ttn/amtic/toxdat.html#data

(26)

Scheffe R.D.; Brook J.R.; Demerjian, K.L. Air Quality Measurements. In: GM Hidy and JR Brook and KL Demerjian and LT Molina and WT Pennell and RD Scheffe (Eds), 2015, Technical Challenges of Multipollutant Air Quality Management, Springer, New York.

(27)

Strum, M.; Scheffe, R. National Review of Ambient Air Toxics Observations, Journal of the Air & Waste Management Association. 2016, 66 (2), 120-133; DOI: 10.1080/10962247.2015.1076538

(28)

U.S. EPA. 2010. Data Quality Evaluation Guidelines for Ambient Air Acrolein Measurements http://www.epa.gov/ttnamti1/files/ambient/airtox/20101217acroleindataqualityeval.pdf.

(29)

Luecken, D.J., Kutzell, W.T.; Strum, M.L.; Pouliot, G.A. Regional sources of atmospheric formaldehyde and acetaldehyde, and implications for atmospheric modeling, Atmospheric Environment. 2012, 47: 477-490.

(30)

Carlton, A.G.; Baker, K.R. Photochemical Modeling of the Ozark Isoprene Volcano: MEGAN, BEIS, and Their Impacts on Air Quality Predictions, | Environ. Sci. Technol. 2011, 45, 4438–4445; dx.doi.org/10.1021/es200050x.

(31)

Warneke, C.; deGouw, J.A.; DelNegro, L.; Brioude, J.; McKeen, S.; Stark, H.; Kustler, W.C.; Goldan, P.D.; Trainer, M.; Fehsenfeld, F.C.; Wiedinmyer, C.; Guenther, A.B.; Hansel, A.;.Wisthaler, A.; Atlas, E.; Holloway, J.S.; Ryerson, T.B.; Peischl, J.; Huey, L.G.; Case Hanks, A.T. Biogenic emission measurement and inventories determination of biogenic emissions in the eastern United States and Texas and comparison with biogenic emission inventories, J. Geophys. Res. 2010, 115, D00F18, oi:10.1029/2009JD012445.

(32)

Guenther, A.; Karl, T.; Harley, P.; Wiedinmyer , C.; Palmer, P.I.; Geron, C. Estimates of of global terrestrial isoprene emissions using MEGAN (model of emissions and gases from nature). Atmos. Chem. Phys. 2006, 6, 3181-3210.

(33)

Hogrefe, C., Isukapalli, S.S., Tang, X., Georgopoulos, P.G., He, S., Zalewsky, E.E., Hao, W., Ku, J-Y., Key, T., Sistla, T. Impact of Biogenic Emission Uncertainties on the Simulated Response of Ozone and Fine Particulate Matter to Anthropogenic Emission Reductions, Journal of the Air & Waste Management Association. 2011, 61, 92-108; DOI: 10.3155/1047-3289.61.1.92.

(34)

Payne-Sturges,D.C.; Burke, T.A.; Breysse, P.; 1 Diener-West, M.; Buckley, T.J. Personal Exposure Meets Risk Assessment: A Comparison of Measured and Modeled Exposures and Risks in an Urban Community. Environmental Health Perspectives. 2004, 112,589-598.

ACS Paragon Plus Environment

30

Page 31 of 31

688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730

Environmental Science & Technology

(35)

Lupo, P.J., Symanski, E. A Comparative Analysis of Modeled and Monitored Ambient Hazardous Air Pollutants in Texas: A Novel Approach Using Concordance Correlation, Journal of the Air & Waste Management Association. 2009, 59, 1278-1286; DOI: 10.3155/1047-3289.59.11.1278.

(36)

Logue, J.M., Small, M.J., Robinson, A.L. Evaluating the national air toxics assessment (NATA): Comparison of predicted and measured air toxics concentrations, risks, and sources in Pittsburgh, Pennsylvania, Atmospheric Environment. 2011, 45, 476-484.

(37)

George, B., Schultz, B.D., Palma, T., Vette, A.F., Whitaker, D.A., Williams, R.W. An evaluation of EPA’s National-Scale Air Toxics Assessment (NATA): Comparison with benzene measurements in Detroit, Michigan, Atmospheric Environment. 2011, 45(19), 3301-3308. Kimbrough, S; Palma, T.; Baldauf, R. Analysis of mobile source air toxics (MSATs) Near-road VOC and carbonyl concentrations. Journal of the Air & Waste Management Association. 2014, 64(3): 349-359.

(38)

(39)

Vennam, L.P., Vizuete, W.; Arunachalam, S. Evaluation of model-predicted Hazardous Air Pollutants (HAPs) near a mid-sized U.S. airport, Atmospheric Environment. 2015, 119, 100-117; dx.doi.org/10.1016/j.atmosenv.2015.08.015.

(40)

Stroud, C.A.; Zaganescu, C.; Chen, J.; McLinden, C.A.; Zhang.; Wang, D. Toxic volatile organic air pollutants across Canada: multi-year concentration trends, regional air quality modeling and source apportionment, J. Atmos. Chem. 2015; DOI 10.1007/s10874-0159319-z.

(41)

Russell, A.; Dennis, R. NARSTO critical review of photochemical models and modeling Atmospheric Environment . 2000, 34, 2283-2324.

(42)

Seigneur, C.; Dennis, R. Atmospheric Modeling. In: GM Hidy and JR Brook and KL Demerjian and LT Molina and WT Pennell and RD Scheffe (Eds) Technical Challenges of Multipollutant Air Quality Management, 2012, Springer, New York.

(43)

National Oceanic and Atmospheric Adminstration, National Centers for Environmental Information, https://www.ncdc.noaa.gov/monitoring-references/maps/us-climateregions.php

(44)

Dunker, A.M. Path-integral method for the source apportionment of photochemical pollutants, Geosci. Model Dev. 2015, 8, 1763–1773.

(45)

Dunker, A.M.; B. Koo.; G. Yarwood. Source Apportionment of the Anthropogenic Increment to Ozone, Formaldehyde, and Nitrogen Dioxide by the Path-Integral Method. Environ. Sci. Technol. 2015, 49, 6751−67.

ACS Paragon Plus Environment

31