Association between spill-related exposure to fine particulate matter and peripheral motor and sensory nerve function among oil spill response and cleanup workers following the Deepwater Horizon oil spill

BACKGROUND

Burning/flaring of oil/gas during the Deepwater Horizon oil spill response and cleanup (OSRC) generated high concentrations of fine particulate matter (PM_2.5). Personnel working on the water during these activities may have inhaled combustion products. Neurologic effects of PM_2.5 have been reported previously but few studies have examined lasting effects following disaster exposures. The association of brief, high exposures and adverse effects on sensory and motor nerve function in the years following exposure have not been examined for OSRC workers.

OBJECTIVES

We assessed the relationship between exposure to burning/flaring-related PM_2.5 and measures of sensory and motor nerve function among OSRC workers.

METHODS

PM_2.5 concentrations were estimated from Gaussian plume dispersion models and linked to self-reported work histories. Quantitative measures of sensory and motor nerve function were obtained 4-6 years after the disaster during a clinical exam restricted to those living close to two clinics in Mobile, AL or New Orleans, LA (n = 3401). We obtained covariate data from a baseline enrollment survey and a home visit, both in 2011-2013. The analytic sample included 1186 participants.

RESULTS

We did not find strong evidence of associations between exposure to PM_2.5 and sensory or motor nerve function, although there was a suggestion of impairment based on single leg stance among individuals with high exposure to PM_2.5. Results were generally consistent whether we examined average or cumulative maximum exposures or removed individuals with the highest crude oil exposures to account for co-pollutant confounding. There was no evidence of exposure-response trends.

IMPACT STATEMENT

Remediating environmental disasters is essential for long-term human and environmental health. During the Deepwater Horizon oil spill disaster, burning and flaring of oil and gas were used to remove these pollutants from the environment, but led to potentially high fine particulate matter exposures for spill response workers working on the water. We investigate the potential adverse effects of these exposures on peripheral nerve function; understanding the potential health harm of remediation tactics is necessary to inform future clean up approaches and protect human health.

Fine particulate matter and incident coronary heart disease events up to 10 years of follow-up among Deepwater Horizon oil spill workers

BACKGROUND

During the 2010 Deepwater Horizon (DWH) disaster, in-situ burning and flaring were conducted to remove oil from the water. Workers near combustion sites were potentially exposed to burning-related fine particulate matter (PM2.5). Exposure to PM2.5 has been linked to increased risk of coronary heart disease (CHD), but no study has examined the relationship among oil spill workers.

OBJECTIVES

To investigate the association between estimated PM2.5 from burning/flaring of oil/gas and CHD risk among the DWH oil spill workers.

METHODS

We included workers who participated in response and cleanup activities on the water during the DWH disaster (N = 9091). PM2.5 exposures were estimated using a job-exposure matrix that linked modelled PM2.5 concentrations to detailed DWH spill work histories provided by participants. We ascertained CHD events as the first self-reported physician-diagnosed CHD or a fatal CHD event that occurred after each worker's last day of burning exposure. We estimated hazard ratios (HR) and 95% confidence intervals (95%CI) for the associations between categories of average or cumulative daily maximum PM2.5 exposure (versus a referent category of water workers not near controlled burning) and subsequent CHD. We assessed exposure-response trends by examining continuous exposure parameters in models.

RESULTS

We observed increased CHD hazard among workers with higher levels of average daily maximum exposure (low vs. referent: HR = 1.26, 95% CI: 0.93, 1.70; high vs. referent: HR = 2.11, 95% CI: 1.08, 4.12; per 10 μg/m3 increase: HR = 1.10, 95% CI: 1.02, 1.19). We also observed suggestively elevated HRs among workers with higher cumulative daily maximum exposure (low vs. referent: HR = 1.19, 95% CI: 0.68, 2.08; medium vs. referent: HR = 1.38, 95% CI: 0.88, 2.16; high vs. referent: HR = 1.44, 95% CI: 0.96, 2.14; per 100 μg/m3-d increase: HR = 1.03, 95% CI: 1.00, 1.05).

CONCLUSIONS

Among oil spill workers, exposure to PM2.5 from flaring/burning of oil/gas was associated with increased risk of CHD.

Associations between airborne crude oil chemicals and symptom-based asthma

RATIONALE

The 2010 Deepwater Horizon (DWH) oil spill response and cleanup (OSRC) workers were exposed to airborne total hydrocarbons (THC), benzene, toluene, ethylbenzene, o-, m-, and p-xylenes and n-hexane (BTEX-H) from crude oil and PM2.5 from burning/flaring oil and natural gas. Little is known about asthma risk among oil spill cleanup workers.

OBJECTIVES

We assessed the relationship between asthma and several oil spill-related exposures including job classes, THC, individual BTEX-H chemicals, the BTEX-H mixture, and PM2.5 using data from the Gulf Long-Term Follow-up (GuLF) Study, a prospective cohort of 24,937 cleanup workers and 7,671 nonworkers following the DWH disaster.

METHODS

Our analysis largely focused on the 19,018 workers without asthma before the spill who had complete exposure, outcome, and covariate information. We defined incident asthma 1-3 years following exposure using both self-reported wheeze and self-reported physician diagnosis of asthma. THC and BTEX-H were assigned to participants based on measurement data and work histories, while PM2.5 used modeled estimates. We used modified Poisson regression to estimate risk ratios (RR) and 95% confidence intervals (CIs) for associations between spill-related exposures and asthma and a quantile-based g-computation approach to explore the joint effect of the BTEX-H mixture on asthma risk.

RESULTS

OSRC workers had greater asthma risk than nonworkers (RR: 1.60, 95% CI: 1.38, 1.85). Higher estimated THC exposure levels were associated with increased risk in an exposure-dependent manner (linear trend test p < 0.0001). Asthma risk also increased with increasing exposure to individual BTEX-H chemicals and the chemical mixture: A simultaneous quartile increase in the BTEX-H mixture was associated with an increased asthma risk of 1.45 (95% CI: 1.35,1.55). With fewer cases, associations were less apparent for physician-diagnosed asthma alone.

CONCLUSIONS

THC and BTEX-H were associated with increased asthma risk defined using wheeze symptoms as well as a physician diagnosis.

Modeled Air Pollution from In Situ Burning and Flaring of Oil and Gas Released Following the Deepwater Horizon Disaster

The GuLF STUDY, initiated by the National Institute of Environmental Health Sciences, is investigating the health effects among workers involved in the oil spill response and clean-up (OSRC) after the Deepwater Horizon (DWH) explosion in April 2010 in the Gulf of Mexico. Clean-up included in situ burning of oil on the water surface and flaring of gas and oil captured near the seabed and brought to the surface. We estimated emissions of PM2.5 and related pollutants resulting from these activities, as well as from engines of vessels working on the OSRC. PM2.5 emissions ranged from 30 to 1.33e6 kg per day and were generally uniform over time for the flares but highly episodic for the in situ burns. Hourly emissions from each source on every burn/flare day were used as inputs to the AERMOD model to develop average and maximum concentrations for 1-, 12-, and 24-h time periods. The highest predicted 24-h average concentrations sometimes exceeded 5000 µg m-3 in the first 500 m downwind of flaring and reached 71 µg m-3 within a kilometer of some in situ burns. Beyond 40 km from the DWH site, plumes appeared to be well mixed, and the predicted 24-h average concentrations from the flares and in situ burns were similar, usually below 10 µg m-3. Structured averaging of model output gave potential PM2.5 exposure estimates for OSRC workers located in various areas across the Gulf. Workers located nearest the wellhead (hot zone/source workers) were estimated to have a potential maximum 12-h exposure of 97 µg m-3 over the 2-month flaring period. The potential maximum 12-h exposure for workers who participated in in situ burns was estimated at 10 µg m-3 over the ~3-month burn period. The results suggest that burning of oil and gas during the DWH clean-up may have resulted in PM2.5 concentrations substantially above the U.S. National Ambient Air Quality Standard for PM2.5 (24-h average = 35 µg m-3). These results are being used to investigate possible adverse health effects in the GuLF STUDY epidemiologic analysis of PM2.5 exposures.

Estimation of Aerosol Concentrations of Oil Dispersants COREXIT™ EC9527A and EC9500A during the Deepwater Horizon Oil Spill Response and Clean-up Operations

The April 2010 Deepwater Horizon drilling unit explosion at the Macondo oil well resulted in the release of approximately 779 million l of oil into the Gulf of Mexico. As part of the response effort to break up oil slicks on the water's surface, 6.81 million l of chemical dispersants COREXIT™ EC9500A and COREXIT™ EC9527A were applied by plane or vessel or injected near the seabed. The GuLF Long-term Follow-up Study is investigating possible adverse health effects of workers involved in the oil spill response and clean-up (OSRC). In this paper, we describe potential dispersant-related air concentrations generated from aerial spraying of dispersants to provide insight as to what concentrations OSRC workers may have been exposed under worst-case conditions. Personal exposure measurement data were not collected. Modeling, therefore, was conducted to estimate airborne concentrations of total aerosol to COREXIT™ EC9527A and EC9500A. Using the AgDISP model, we estimated air concentrations to dispersant total aerosols, defined as all components of the dispersant including active ingredients, surfactants, and water, resulting from aerial and vessel applications, as average 1-h and 2-min concentrations. For comparison, 1-h air concentrations associated with aerial spraying were estimated using another model, AERMOD. At 152 m horizontal to the flight path, average 1-h total aerosol concentrations associated with aerial applications were estimated to be as high as 49.3 µg m-3 (9527A) and 45.4 µg m-3 (9500A), and both decreased with increased distance from the flight line. The estimates for spraying 9500A from vessels indicated that total aerosol concentrations were potentially as high as 0.33 µg m-3 at 10 m from the nozzles. These results suggest that personal exposures to dispersant aerosols were negligible.

Assessing Exposures from the Deepwater Horizon Oil Spill Response and Clean-up

The GuLF Study is investigating adverse health effects from work on the response and clean-up after the Deepwater Horizon explosion and oil release. An essential and necessary component of that study was the exposure assessment. Bayesian statistical methods and over 135 000 measurements of total hydrocarbons (THC), benzene, ethylbenzene, toluene, xylene, and n-hexane (BTEX-H) were used to estimate inhalation exposures to these chemicals for >3400 exposure groups (EGs) formed from three exposure determinants: job/activity/task, location, and time period. Recognized deterministic models were used to estimate airborne exposures to particulate matter sized 2.5 µm or less (PM2.5) and dispersant aerosols and vapors. Dermal exposures were estimated for these same oil-related substances using a model modified especially for this study from a previously published model. Exposures to oil mist were assessed using professional judgment. Estimated daily THC arithmetic means (AMs) were in the low ppm range (<25 ppm), whereas BTEX-H exposures estimates were generally <1000 ppb. Potential 1-h PM2.5 air concentrations experienced by some workers may have been as high as 550 µg m-3. Dispersant aerosol air concentrations were very low (maximum predicted 1-h concentrations were generally <50 µg m-3), but vapor concentrations may have exceeded occupational exposure excursion guidelines for 2-butoxyethanol under certain circumstances. The daily AMs of dermal exposure estimates showed large contrasts among the study participants. The estimates are being used to evaluate exposure-response relationships in the GuLF Study.

Association of Deepwater Horizon Oil Spill Response and Cleanup Work With Risk of Developing Hypertension

IMPORTANCE

Exposure to hydrocarbons, fine particulate matter (PM2.5), and other chemicals from the April 20, 2010, Deepwater Horizon disaster may be associated with increased blood pressure and newly detected hypertension among oil spill response and cleanup workers.

OBJECTIVE

To determine whether participation in cleanup activities following the disaster was associated with increased risk of developing hypertension.

DESIGN, SETTING, AND PARTICIPANTS

This cohort study was conducted via telephone interviews and in-person home exams. Participants were 6846 adults who had worked on the oil spill cleanup (workers) and 1505 others who had completed required safety training but did not do cleanup work (nonworkers). Eligible participants did not have diagnosed hypertension at the time of the oil spill. Statistical analyses were performed from June 2018 to December 2021.

EXPOSURES

Engagement in cleanup activities following the Deepwater Horizon oil spill disaster, job classes, quintiles of cumulative total hydrocarbons exposure level, potential exposure to burning or flaring oil, and estimated PM2.5 were examined.

MAIN OUTCOMES AND MEASURES

Systolic and diastolic blood pressure measurements were collected during home exams from 2011 to 2013 using automated oscillometric monitors. Newly detected hypertension was defined as antihypertensive medication use or elevated blood pressure since the spill. Log binomial regression was used to calculate prevalence ratios (PR) and 95% CIs for associations between cleanup exposures and hypertension. Multivariable linear regression was used to estimate exposure effects on continuous blood pressure levels.

RESULTS

Of 8351 participants included in this study, 6484 (77.6%) were male, 517 (6.2%) were Hispanic, 2859 (34.2%) were non-Hispanic Black, and 4418 (52.9%) were non-Hispanic White; the mean (SD) age was 41.9 (12.5) years at enrollment. Among workers, the prevalence of newly detected hypertension was elevated in all quintiles (Q) of cumulative total hydrocarbons above the first quintile (PR for Q3, 1.29 [95% CI, 1.13-1.46], PR for Q4, 1.25 [95% CI, 1.10-1.43], and PR for Q5, 1.31 [95% CI, 1.15-1.50]). Both exposure to burning and/or flaring oil and gas (PR, 1.16 [95% CI, 1.02-1.33]) and PM2.5 from burning (PR, 1.26 [95% CI, 0.89-1.71]) for the highest exposure category were associated with increased risk of newly detected hypertension, as were several types of oil spill work including cleanup on water (PR, 1.34 [95% CI, 1.08-1.66]) and response work (PR, 1.51 [95% CI, 1.20-1.90]).

CONCLUSIONS AND RELEVANCE

Oil spill exposures were associated with newly detected hypertension after the Deepwater Horizon disaster. These findings suggest that blood pressure screening should be considered for workers with occupational hydrocarbon exposures.

Fine Particulate Matter and Lung Function among Burning-Exposed Deepwater Horizon Oil Spill Workers

BACKGROUND

During the 2010 Deepwater Horizon (DWH) disaster, controlled burning was conducted to remove oil from the water. Workers near combustion sites were potentially exposed to increased fine particulate matter [with aerodynamic diameter ≤2.5μm (PM2.5)] levels. Exposure to PM2.5 has been linked to decreased lung function, but to our knowledge, no study has examined exposure encountered in an oil spill cleanup.

OBJECTIVE

We investigated the association between estimated PM2.5 only from burning/flaring of oil/gas and lung function measured 1-3 y after the DWH disaster.

METHODS

We included workers who participated in response and cleanup activities on the water during the DWH disaster and had lung function measured at a subsequent home visit (n=2,316). PM2.5 concentrations were estimated using a Gaussian plume dispersion model and linked to work histories via a job-exposure matrix. We evaluated forced expiratory volume in 1 s (FEV1; milliliters), forced vital capacity (FVC; milliliters), and their ratio (FEV1/FVC; %) in relation to average and cumulative daily maximum exposures using multivariable linear regressions.

RESULTS

We observed significant exposure-response trends associating higher cumulative daily maximum PM2.5 exposure with lower FEV1 (p-trend=0.04) and FEV1/FVC (p-trend=0.01). In comparison with the referent group (workers not involved in or near the burning), those with higher cumulative exposures had lower FEV1 [-166.8mL, 95% confidence interval (CI): -337.3, 3.7] and FEV1/FVC (-1.7, 95% CI: -3.6, 0.2). We also saw nonsignificant reductions in FVC (high vs. referent: -120.9, 95% CI: -319.4, 77.6; p-trend=0.36). Similar associations were seen for average daily maximum PM2.5 exposure. Inverse associations were also observed in analyses stratified by smoking and time from exposure to spirometry and when we restricted to workers without prespill lung disease.

CONCLUSIONS

Among oil spill workers, exposure to PM2.5 specifically from controlled burning of oil/gas was associated with significantly lower FEV1 and FEV1/FVC when compared with workers not involved in burning. https://doi.org/10.1289/EHP8930.

Analysis of polycyclic aromatic hydrocarbons (PAHs) in air using passive sampling calibrated with active measurements

There are limited ambient air measurements of extended (beyond EPA Priority 16) lists of polycyclic aromatic hydrocarbons (PAHs). We measured air concentrations of 45 PAHs using passive and active air sampling at 15 sites in a central urban community and one rural site for two years. Passive sampling was conducted with cylindrical XAD-based samplers deployed to capture spatial variability. High volume active samplers with quartz fiber filters for particles and XAD-4 absorbent for gases were deployed at two urban sites and the rural site to calibrate the passive measurements directly. Estimated passive sampling rates (PSRs) were evaluated as functions of meteorological data, seasons, locations, study year, and compared with other studies. Possible particle collection by the passive samplers was evaluated using a variety of particle measurements (TSP, PM₁₀, PM_2.5 and ultrafines <100 nm). Total PAHs were statistically associated with ultrafine particle concentrations and to a lesser extent PM_2.5 and PM₁₀, but not TSP. PSRs were more variable when PAH mass loadings were lower and near method detection limits; this occurred more often at the rural site. The PSRs were not statistically associated with meteorological conditions in this study, but wind speed had the highest potential to impact PSR results. The resulting passive PAH measurements are reported with respect to proximity to major roadways and other known air emissions types. PSRs were quantifiable for some PAHs that were found predominantly in the particulate phase in active sampling. This information, together with particle fraction calculations from active sampling, were used to estimate the particulate PAH capture of the passive sampler. Summed PAH (∑PAH) passive concentrations were measured within the range of 10-265 ng/m3, with the highest concentrations from naphthalene and the lowest detected concentrations from anthracene. These results indicated a stronger seasonal signal within 200 m of a major roadway.

Asthma exacerbations and traffic: examining relationships using link-based traffic metrics and a comprehensive patient database

BACKGROUND

The Rochester Epidemiology Project (REP) is a unique community-based medical record data linkage system that provides individual patient address, diagnosis and visit information for all hospitalizations, as well as emergency department, urgent care and outpatient clinic visits for asthma. Proximity to traffic is known to be associated with asthma exacerbations and severity. Our null hypothesis was that there is no association between residential proximity to traffic and asthma exacerbations over eleven years of REP data.

METHODS

Spatial coordinates of the homes of 19,915 individuals diagnosed with asthma were extracted from the REP database. Three metrics of traffic exposure at residences were calculated from link-based traffic count data. We used exploratory statistics as well as logistic and Poisson regression to examine associations between three traffic metrics at the home address and asthma exacerbations.

RESULTS

Asthma exacerbations increased as traffic levels near the home increased. Proximity to traffic was a significant predictor of asthma exacerbations in logistic and Poisson regressions controlling for age, gender and block group poverty.

CONCLUSIONS

Over eleven years in a comprehensive county-wide data set of asthma patients, and after controlling for demographic effects, we found evidence that living in proximity to traffic increased the risk of asthma exacerbations.

Traffic, air pollution, minority and socio-economic status: addressing inequities in exposure and risk

Higher levels of nearby traffic increase exposure to air pollution and adversely affect health outcomes. Populations with lower socio-economic status (SES) are particularly vulnerable to stressors like air pollution. We investigated cumulative exposures and risks from traffic and from MNRiskS-modeled air pollution in multiple source categories across demographic groups. Exposures and risks, especially from on-road sources, were higher than the mean for minorities and low SES populations and lower than the mean for white and high SES populations. Owning multiple vehicles and driving alone were linked to lower household exposures and risks. Those not owning a vehicle and walking or using transit had higher household exposures and risks. These results confirm for our study location that populations on the lower end of the socio-economic spectrum and minorities are disproportionately exposed to traffic and air pollution and at higher risk for adverse health outcomes. A major source of disparities appears to be the transportation infrastructure. Those outside the urban core had lower risks but drove more, while those living nearer the urban core tended to drive less but had higher exposures and risks from on-road sources. We suggest policy considerations for addressing these inequities.

Quantifying traffic exposure

Living near traffic adversely affects health outcomes. Traffic exposure metrics include distance to high-traffic roads, traffic volume on nearby roads, traffic within buffer distances, measured pollutant concentrations, land-use regression estimates of pollution concentrations, and others. We used Geographic Information System software to explore a new approach using traffic count data and a kernel density calculation to generate a traffic density surface with a resolution of 50 m. The density value in each cell reflects all the traffic on all the roads within the distance specified in the kernel density algorithm. The effect of a given roadway on the raster cell value depends on the amount of traffic on the road segment, its distance from the raster cell, and the form of the algorithm. We used a Gaussian algorithm in which traffic influence became insignificant beyond 300 m. This metric integrates the deleterious effects of traffic rather than focusing on one pollutant. The density surface can be used to impute exposure at any point, and it can be used to quantify integrated exposure along a global positioning system route. The traffic density calculation compares favorably with other metrics for assessing traffic exposure and can be used in a variety of applications.

Validation of a novel air toxic risk model with air monitoring

Three modeling systems were used to estimate human health risks from air pollution: two versions of MNRiskS (for Minnesota Risk Screening), and the USEPA National Air Toxics Assessment (NATA). MNRiskS is a unique cumulative risk modeling system used to assess risks from multiple air toxics, sources, and pathways on a local to a state-wide scale. In addition, ambient outdoor air monitoring data were available for estimation of risks and comparison with the modeled estimates of air concentrations. Highest air concentrations and estimated risks were generally found in the Minneapolis-St. Paul metropolitan area and lowest risks in undeveloped rural areas. Emissions from mobile and area (nonpoint) sources created greater estimated risks than emissions from point sources. Highest cancer risks were via ingestion pathway exposures to dioxins and related compounds. Diesel particles, acrolein, and formaldehyde created the highest estimated inhalation health impacts. Model-estimated air concentrations were generally highest for NATA and lowest for the AERMOD version of MNRiskS. This validation study showed reasonable agreement between available measurements and model predictions, although results varied among pollutants, and predictions were often lower than measurements. The results increased confidence in identifying pollutants, pathways, geographic areas, sources, and receptors of potential concern, and thus provide a basis for informing pollution reduction strategies and focusing efforts on specific pollutants (diesel particles, acrolein, and formaldehyde), geographic areas (urban centers), and source categories (nonpoint sources). The results heighten concerns about risks from food chain exposures to dioxins and PAHs. Risk estimates were sensitive to variations in methodologies for treating emissions, dispersion, deposition, exposure, and toxicity.

Cumulative risk assessment and environmental equity in air permitting: interpretation, methods, community participation and implementation of a unique statute

In 2008, the statute authorizing the Minnesota Pollution Control Agency (MPCA) to issue air permits was amended to include a unique requirement to analyze and consider “cumulative levels and effects of past and current environmental pollution from all sources on the environment and residents of the geographic area within which the facility’s emissions are likely to be deposited.” Data describing the Statute Area suggest it is challenged by environmental and socioeconomic concerns, i.e., concerns which are often described by the phrase ‘environmental equity’. With input from diverse stakeholders, the MPCA developed a methodology for implementing a cumulative levels and effects analysis when issuing air permits in the designated geographic area. A Process Document was created defining explicit steps a project proposer must complete in the analysis. An accompanying Reference Document compiles all available environmental health data relevant to the Statute Area that could be identified. The final cumulative levels and effects methodology is organized by health endpoint and identifies hazard, exposure and health indices that require further evaluation. The resulting assessment is summarized and presented to decision makers for consideration in the regulatory permitting process. We present a description of the methodology followed by a case study summary of the first air permit processed through the “cumulative levels and effects analysis”.

Modeling community asbestos exposure near a vermiculite processing facility: Impact of human activities on cumulative exposure

Contaminated vermiculite ore from Libby, Montana was processed in northeast Minneapolis from 1936 to 1989 in a densely populated urban residential neighborhood, resulting in non-occupational exposure scenarios from plant stack and fugitive emissions as well as from activity-based scenarios associated with use of the waste rock in the surrounding community. The objective of this analysis was to estimate potential cumulative asbestos exposure for all non-occupationally exposed members of this community. Questionnaire data from a neighborhood-exposure assessment ascertained frequency of potential contact with vermiculite processing waste. Monte Carlo simulation was used to develop exposure estimates based on activity-based concentration estimates and contact durations for four scenarios: S1, moved asbestos-contaminated waste; S2, used waste at home, on lawn or garden; S3, installed/removed vermiculite insulation; S4, played in or around waste piles at the plant. The simulation outputs were combined with air-dispersion model results to provide total cumulative asbestos exposure estimates for the cohort. Fiber emissions from the plant were the largest source of exposure for the majority of the cohort, with geometric mean cumulative exposures of 0.02 fibers/cc × month. The addition of S1, S2 and S3 did not significantly increase total cumulative exposure above background exposure estimates obtained from dispersion modeling. Activity-based exposures were a substantial contributor to the upper end of the exposure distribution: 90th percentile S4 exposure estimates are ∼10 times higher than exposures from plant emissions. Pile playing is the strongest source of asbestos exposure in this cohort, with other activity scenarios contributing less than from plant emissions.

Multipathway screening factors for assessing risks from ingestion exposures to air pollutants

Regulatory agencies are frequently called upon to assess the potential for significant environmental impacts from air pollution emissions. These assessments often entail air dispersion modeling to estimate air concentrations that can be compared with standards or health benchmarks. Some air pollutants can also impact human health through pathways in media besides air. Risk assessment models are available that consider pollutant deposition, movement, uptake, and other processes on land and water and in biota, but they are typically effort-intensive. A screening-level assessment of potential multipathway effects would be useful. We developed multipathway screening factors (MPSFs) that can be applied to inhalation risk estimates to give screening estimates of risks via ingestion pathways. The MPSFs were generated using a generic multipathway risk assessment, consisting of air dispersion and deposition modeling followed by risk modeling for 42 persistent, bioaccumulative air pollutants. MPSFs are defined as the ratio of ingestion risks to inhalation risks. We report here the results of a sensitivity analysis that evaluates the effects on the MPSF ratio of varying inputs to the air dispersion and deposition modeling analysis. Model input parameters were systematically varied and multipathway risks recalculated. From the sensitivity analysis results, reasonable upper-bound values for the ratio of ingestion risks to inhalation risks for each pollutant were selected. The particle size distribution and the method of calculating particle deposition had the most disproportionate effect on inhalation versus ingestion risks and the greatest effect on MPSFs. Risk calculations are often done at the points of maximum air concentration and maximum deposition. In this study, the MPSFs were usually highest at the location of the maximum inhalation risk.

Community Exposure to Asbestos From a Vermiculite Exfoliation Plant in NE Minneapolis

Western Mineral Products/W. R. Grace operated a vermiculite plant in a mixed industrial/residential area of northeast Minneapolis from 1936 to 1989. The plant processed vermiculite ore contaminated with amphibole asbestos from a mine in Libby, MT. Air monitoring in the early 1970s found fiber concentrations in excess of 10 fibers per cubic centimeter of air (f/cc), indicating that worker exposure to asbestos was occasionally 100 times the current occupational standard. Residents of the surrounding community also had direct contact with vermiculite processing wastes (containing up to 10% amphibole asbestos) that were made freely available. Children played on waste piles and neighborhood residents hauled the wastes away for home use. In total, 259 contaminated residential properties have been found to date. Reported emission factors and plant process data were used as inputs to model airborne emissions from the plant over several operating scenarios using the U.S. Environmental Protection Agency (EPA) ISC-Prime model. Results estimate short-term air concentrations of asbestos fibers in residential areas nearest the plant may have at times exceeded current occupational standards. Exposure estimates for other pathways were derived primarily from assessments done in Libby by the U.S. EPA. The Northeast Minneapolis Community Vermiculite Investigation (NMCVI) was conducted by the Minnesota Department of Health to identify and characterize the exposures of a cohort of over 6000 people who live or lived in Northeast Minneapolis and may have been exposed to asbestos. This cohort is now being investigated in a respiratory health screening study conducted by the University of Minnesota and the Minnesota Department of Health.

Relationships between personal, indoor, and outdoor exposures to trace elements in PM(2.5)

Twenty-four hour average fine particle concentrations of 23 trace elements (TEs) were measured concurrently in (a) ambient air in three urban neighborhoods (Battle Creek-BCK; East St. Paul-ESP; and Phillips-PHI), (b) air inside residences of participants, and (c) personal air near the breathing zone of healthy, non-smoking adults. The outdoor (O), indoor (I), and personal (P) samples were collected in the Minneapolis/St. Paul metropolitan area over three seasons (Spring, Summer, Fall) using either the federal reference (O) or inertial impactor (I,P) inlets to collect PM(2.5). In addition to descriptive statistics, a hierarchical, mixed-effects statistical model was used to estimate the mutually adjusted effects of monitor location, community, and season on mean differences between monitoring locations while accounting for within-subject and within-monitoring period correlation. The relationships among P, I, and O concentrations varied across TEs. The O concentrations were usually higher than P or I for elements like Ca and Al that originate mainly from entrained crustal material, while P concentrations were often highest for other elements with non-crustal sources. Unadjusted mixed model results demonstrated that O monitors more frequently underestimated than overestimated P TE exposures for elements associated with non-crustal sources. This finding was true even though the O TE measurements were taken in the same neighborhoods as the P and I measurements. Further adjustment for community or season effects in the mixed models reduced the number of significant O-P and O-I differences compared to unadjusted models, but still indicated a tendency for underestimation of personal and indoor TE exposures by central site monitors, particularly in the PHI community. These results indicate that community and season are important covariates for developing long term TE exposure estimates, and that personal exposure to trace elements in PM(2.5) is likely to be underestimated by outdoor central site monitors.

Volatile organic compound emissions from dry mill fuel ethanol production

Ethanol fuel production is growing rapidly in the rural Midwest, and this growth presents potential environmental impacts. In 2002, the U.S. Environmental Protection Agency (EPA) and the Minnesota Pollution Control Agency (MPCA) entered into enforcement actions with 12 fuel ethanol plants in Minnesota. The enforcement actions uncovered underreported emissions and resulted in consent decrees that required pollution control equipment be installed. A key component of the consent decrees was a requirement to conduct emissions tests for volatile organic compounds (VOCs) with the goal of improving the characterization and control of emissions. The conventional VOC stack test method was thought to underquantify total VOC emissions from ethanol plants. A hybrid test method was also developed that involved quantification of individual VOC species. The resulting database of total and speciated VOC emissions from 10 fuel ethanol plants is relatively small, but it is the most extensive to date and has been used to develop and gauge compliance with permit limits and to estimate health risks in Minnesota. Emissions were highly variable among facilities and emissions units. In addition to the variability, the small number of samples and the presence of many values below detection limits complicate the analysis of the data. To account for these issues, a nested bootstrap procedure on the Kaplan-Meier method was used to calculate means and upper confidence limits. In general, the fermentation scrubbers and fluid bed coolers emitted the largest mass of VOC emissions. Across most facilities and emissions units ethanol was the pollutant emitted at the highest rate. Acetaldehyde, acetic acid, and ethyl acetate were also important emissions from some units. Emissions of total VOCs, ethanol, and some other species appeared to be a function of the beer feed rate, although the relationship was not reliable enough to develop a production rate-based emissions factor.

Estimating volatile organic compound concentrations in selected microenvironments using time-activity and personal exposure data

Repeated measures of personal exposure to 14 volatile organic compounds (VOC) were obtained over 3 seasons for 70 healthy, nonsmoking adults living in Minneapolis-St. Paul. Matched data were also available for participants' time-activity patterns, and measured VOC concentrations outdoors in the community and indoors in residences. A novel modeling approach employing hierarchical Bayesian techniques was used to estimate VOC concentrations (posterior mode) and variability (credible intervals) in five microenvironments: (1) indoors at home; (2) indoors at work/school; (3) indoors in other locations; (4) outdoors in any location; and (5) in transit. Estimated concentrations tended to be highest in "other" indoor microenvironments (e.g., grocery stores, restaurants, shopping malls), intermediate in the indoor work/school and residential microenvironments, and lowest in the outside and in-transit microenvironments. Model estimates for all 14 VOC were reasonable approximations of measured median concentrations in the indoor residential microenvironment. The largest predicted contributor to cumulative (2-day) personal exposure for all 14 VOC was the indoor residential environment. Model-based results suggest that indoors-at-work/school and indoors-at-other-location microenvironments were the second or third largest contributors for all VOC, while the outside-in-any-location and in-transit microenvironments appeared to contribute negligibly to cumulative personal exposure. Results from a mixed-effects model indicate that being in or near a garage increased personal exposure to o-xylene, m/p-xylene, benzene, ethylbenzene, and toluene, and leaving windows and doors at home open for 6 h or more decreased personal exposure to 13 of 14 VOC, all except trichloroethylene.

A field comparison of volatile organic compound measurements using passive organic vapor monitors and stainless steel canisters

Concurrent field measurements of 10 volatile organic compounds (VOCs) were made using passive diffusion-based organic vapor monitors (OVMs) and the U.S. Federal Reference Method, which comprises active monitoring with stainless steel canisters (CANs). Measurements were obtained throughout a range of weather conditions, repeatedly over the course of three seasons, and at three different locations in the Minneapolis/St. Paul metropolitan area. Ambient concentrations of most VOCs as measured by both methods were low compared to those of other large metropolitan areas. For some VOCs a considerable fraction of measurements was below the detection limit of one or both methods. The observed differences between the two methods were similar across measurement sites, seasons, and meteorological variables. A Bayesian analysis with uniform priors on the differences was applied, with accommodation of sometimes heavy censoring (nondetection) in either device. The resulting estimates of bias and standard deviation of the OVM relative to the CAN were computed by tertile of the canister-measured concentration. In general, OVM and CAN measurements were in the best agreement for benzene and other aromatic compounds with hydrocarbon additions (ethylbenzene, toluene, and xylenes). The two methods were not in such good agreement for styrene and halogenated compounds (carbon tetrachloride, p-dichlorobenzene, methylene chloride, and trichloroethylene). OVMs slightly overestimated benzene concentrations and carbon tetrachloride at low concentrations, but in all other cases where significant differences were found, OVMs underestimated relative to canisters. Our study indicates that the two methods are in agreement for some compounds, but not all. We provide data and interpretation on the relative performance of the two VOC measurement methods, which facilitates intercomparisons among studies.

Evaluating differences between measured personal exposures to volatile organic compounds and concentrations in outdoor and indoor air

Accurate estimation of human exposures to volatile organic compounds (VOCs) is a key element of strategies designed to protect public health from the adverse effects of hazardous air pollutants. The focus here is on examining the capability of three different exposure metrics (outdoor community concentrations, indoor residential concentrations, and a simple time-weighted model) to estimate observed personal exposures to 14 VOCs. The analysis is based on 2-day average concentrations of individual VOCs measured concurrently in outdoor (O) air in three urban neighborhoods, indoor (I) air in participant's residences, and personal (P) air near the breathing zone of 71 healthy, nonsmoking adults. A median of four matched P-I-O samples was collected for each study participant in Minneapolis/St. Paul over three seasons (spring, summer, and fall) in 1999 using charcoal-based passive air samplers (3M model 3500 organic vapor monitors). Results show a clear pattern for the 14 VOCs, with P > I > O concentrations. Intra-individual variability typically spanned at least an order of magnitude, and inter-individual variability spanned 2 or more orders of magnitude for each of the 14 VOCs. Although both O and I concentrations generally underestimated personal exposures, I concentrations provided a substantially better estimate of measured P concentrations. Mean squared error (MSE) as well as correlation measures were used to assess estimator performance at the subject-specific level, and hierarchical, mixed effects models were used to estimate the bias and variance components of MSE by tertile of personal exposure. Bias and variance both tended to increase in the upper third of the P exposure distribution for O versus P and I versus P. A simple time-weighted model incorporating measured concentrations in both outdoor community air and indoor residential air provided no improvement over I concentration alone for the estimation of P exposure.

Comparing air dispersion model predictions with measured concentrations of VOCs in urban communities

Air concentrations of nine volatile organic compounds were measured over 48-h periods at 23 locations in three communities in the Minneapolis-St. Paul metropolitan area. Concentrations at the same times and locations were modeled using a standard regulatory air dispersion model (ISCST3). The goal of the study was to evaluate model performance by comparing predictions with measurements using linear regression and estimates of bias. The modeling, done with mobile and area source emissions resolved to the census tract level and characterized as model area sources, represents an improvement over large-scale airtoxics modeling analyses done to date. Despite the resolved spatial scale, the model did not fully capture the spatial resolution in concentrations in an area with a sharp gradient in emissions. In a census tract with a major highway at one end of the tract (i.e., uneven distribution of emissions within the tract), model predictions atthe opposite end of the tract overestimated measured concentrations. This shortcoming was seen for pollutants emitted mainly by mobile sources (benzene, ethylbenzene, toluene, and xylenes). We suggest that major highways would be better characterized as line sources. The model also failed to fully capture the temporal variability in concentrations, which was expected since the emissions inventory comprised annual average values. Based on our evaluation metrics, model performance was best for pollutants emitted mainly from mobile sources and poorest for pollutants emitted mainlyfrom area sources. Important sources of error appeared to be the source characterization (especially location) and emissions quantification. We expect that enhancements in the emissions inventory would give the greatest improvement in results. As anticipated for a Gaussian plume model, performance was dramatically better when compared to measurements that were not matched in space or time. Despite the limitations of our analysis, we found thatthe regulatory air dispersion model was generally able to predict space and time matched 48-h average ambient concentrations of VOC species within a factor of 2 on average, results that meet regulatory model acceptance criteria.

Comparison of personal, indoor, and outdoor exposures to hazardous air pollutants in three urban communities

Two-day average concentrations of 15 individual volatile organic compounds (VOCs) were measured concurrently in (a) ambient air in three urban neighborhoods, (b) air inside residences of participants, and (c) personal air near the breathing zone of 71 healthy, nonsmoking adults. The outdoor (O), indoor (I), and personal (P) samples were collected in the Minneapolis/St. Paul metropolitan area over three seasons (spring, summer, and fall) in 1999 using charcoal-based passive air samplers (3M model 3500 organic vapor monitors). A hierarchical, mixed-effects statistical model was used to estimate the mutually adjusted effects of monitor location, community, and season while accounting for within-subject and within-time-index (monitoring period) correlation. Outdoor VOC concentrations were relatively low compared to many other urban areas, and only minor seasonal differences were observed. A consistent pattern of P > I > O was observed across both communities and seasons for 13 of 15 individual VOCs (exceptions were carbon tetrachloride and chloroform). Results indicate that ambient VOC measurements at central monitoring sites can seriously underestimate actual exposures for urban residents, even when the outdoor measurements are taken in their own neighborhoods.

Analysis of air toxics emission inventory: inhalation toxicity-based ranking

Air toxics emission inventories play an important role in air quality regulatory activities. Recently, Minnesota Pollution Control Agency (MPCA) staff compiled a comprehensive air toxics emission inventory for 1996. While acquiring data on the mass of emissions is a necessary first step, equally important is developing information on the potential toxicity of the emitted pollutants. To account for the toxicity of the pollutants in the emission inventory, inhalation health benchmarks for acute effects, chronic effects, and cancer were used to weight the mass of emissions. The 1996 Minnesota emissions inventory results were ranked by mass of emissions and by an index comprised of emissions divided by health benchmarks. The results show that six of eight pollutants ranked highest by toxicity were also the pollutants of concern indicated in environmental monitoring data and modeling data. Monitoring data and modeling results did not show high impacts of the other two pollutants that were identified by the toxicity-based emission ranking method. The biggest limitation in this method is the lack of health benchmark values for many pollutants. Despite uncertainties and limited information, this analysis provides useful information for further targeting pollutants and source categories for control.

An assessment of air toxics in Minnesota

We used monitoring and modeling to assess the concentrations of air toxics in the state of Minnesota. Model-predicted concentrations for 148 hazardous air pollutants were from the U.S. Environmental Protection Agency Cumulative Exposure Project (1990 data). Monitoring data consisted of samples of volatile organic compounds, carbonyls, and particulate matter [Less than and equal to] 10 microm in aerodynamic diameter collected at 25 sites throughout the state for varying periods of time (up to 8 years; 1991-1998). Ten pollutants exceeded health benchmark values at one or more sites by modeling, monitoring, or both (including acrolein, arsenic, benzene, 1,3-butadiene, carbon tetrachloride, chromium, chloroform, ethylene dibromide, formaldehyde, and nickel). Polycyclic organic matter also exceeded the benzo[a]pyrene health benchmark value assumed to represent this class of pollutants. The highest modeled and monitored concentrations of most pollutants were near the center of the Minneapolis-St. Paul metropolitan area; however, many smaller cities throughout the state also had elevated concentrations. Where direct comparisons were possible, monitored values often tended to exceed model estimates. Upper-bound excess lifetime inhalation cancer risks were estimated to range from 2.7 [times] 10(-5) to 140. 9 [times] 10(-5) (modeling) and 4.7 [times] 10(-5) to 11.0 [times] 10(-5) (using a smaller set of monitored carcinogens). Screening noncancer hazard indices summed over all end points ranged from 0.2 to 58.1 (modeling) and 0.6 to 2.0 (with a smaller set of monitored pollutants). For common sets of pollutants, the concentrations, cancer risks, and noncancer hazard indices were comparable between model-based estimates and monitored values. The inhalation cancer risk was apportioned to mobile sources (54%), area sources (22%), point sources (12%), and background (12%). This study provides evidence that air toxics are a public health concern in Minnesota.

Comparison of short-term variations (15-minute averages) in outdoor and indoor PM2.5 concentrations

Measurements of 15-min average PM2.5 concentrations were made with a real-time light-scattering instrument at both outdoor (central monitoring sites in three communities) and indoor (residential) locations over two seasons in the Minneapolis-St. Paul metropolitan area. These data are used to examine within-day variability of PM2.5 concentrations indoors and outdoors, as well as matched indoor-to-outdoor (I/O) ratios. Concurrent gravimetric measurements of 24-hr average PM2.5 concentrations were also obtained as a way to compare real-time measures with this more traditional metric. Results indicate that (1) within-day variability for both indoor and outdoor 15-min average PM2.5 concentrations was substantial and comparable in magnitude to day-to-day variability for 24-hr average concentrations; (2) some residences exhibited substantial variability in indoor aerosol characteristics from one day to the next; (3) peak values for indoor short-term (15-min) average PM2.5 concentrations routinely exceeded 24-hr average outdoor values by factors of 3-4; and (4) relatively strong correlations existed between indoor and outdoor PM2.5 concentrations for both 24-hr and 15-min averages.