Cancer Risk Estimates from AirToxScreen Data.
Air pollution cancer risk for a given carcinogen is estimated by multiplying the lifetime exposure concentration by the IUR value. AirToxScreen cancer risk estimates are based on modeled exposure concentrations and IUR values listed in USEPA Integrated Risk Information System (IRIS) database (
40). We use these same IUR values for estimating cancer risk from our measured concentrations (presented in the following section). AirToxScreen provides cancer risk estimates for 72 individual HAPs, including primary and secondary VOCs and particulate metals.
To contextualize our measurements, we first present an analysis of AirToxScreen risk estimates for the BR-NO corridor (outlined in green on map in
Fig. 1A), which consists of 530 total census tracts.
Fig. 1A shows a map of total estimated cancer risk from AirToxScreen for the BR-NO corridor at the census tract level. Nearly all (97%) of the 530 census tracts in the BR-NO corridor are above the 50th percentile nationwide for total cancer risk. Within the corridor, there are many areas of very-high risk relative to the rest of the country: 32% of corridor tracts are above the 95th percentile nationwide, and 7% of corridor tracts are above the 99th percentile. Two census tracts (0.4% of total in corridor) have estimated tract-level cancer risk above the maximum “acceptable” level of cancer risk advised in USEPA Superfund Site remediation of 100-in-one million (
41). Many tracts have levels of cancer risk just below this “acceptable” threshold as well: For example, 40 tracts (7.5% of total in corridor) have risk levels higher than 50-in-one million.
Despite AirToxScreen identifying the area as a clear hot spot for cancer risk, there is still a large degree of spatial variability in total cancer risk estimates at the census tract level. This variability spans an order of magnitude across the corridor, which is shown in
Fig. 1A. The tract with the highest cancer risk estimate was 248-in-one million, and the lowest was 21-in-one million. This variability is the result of substantial spatial variation in individual HAPs concentrations and associated risk.
Fig. 1B shows a map of census tracts colored according to which HAP has the single largest individual risk estimate for that tract. Four different HAPs are shown as having the largest individual contribution to total cancer risk for these tracts: formaldehyde, ethylene oxide, chloroprene, and benzidine.
For a large majority (95%) of tracts in the BR-NO corridor, formaldehyde has the highest individual contribution to total risk of all HAPs. Formaldehyde is linked to nasal cancers and myeloid leukemia (
31). Statewide, an even higher percentage (97%) of all tracts have formaldehyde as the largest contributor to risk (
SI Appendix, Fig. S2A). However, the magnitude of formaldehyde-attributed cancer risk is relatively uniform across the corridor; the full range spans only a factor of 1.59 (17- to 26-in-one million risk). Formaldehyde’s uniformity is explained by its formation largely from secondary photochemical and biogenic sources, not specifically industrial primary emissions (
42). Despite the homogeneity within the corridor, AirToxScreen estimates of formaldehyde cancer risk are more variable statewide (
SI Appendix, Fig. S2D), and the highest-estimated formaldehyde risk is outside of the industrial corridor entirely; indeed the maximum formaldehyde risk in the corridor only corresponds to the 37th percentile statewide.
Ethylene oxide is the largest individual contributor to total cancer risk for 14 total tracts in the corridor (3% of tracts). Ethylene oxide is thought to be a primary-only pollutant largely from industrial point sources and is linked to lymphoid and breast cancers (
29). The tracts with ethylene oxide as the largest contributor include parts of four different parishes (Ascension, Iberville, St. Charles, St. John the Baptist), with a maximum individual tract contribution of 70-in-one million, corresponding to an exposure concentration of roughly 8 ppt. The high-ethylene oxide tracts cover a large fraction of the area (
1,000 km
2) of the BR-NO corridor. There are additionally five tracts outside of the BR-NO corridor where ethylene oxide is the highest single contributor to total cancer risk, each of which are in the Lake Charles area where there are multiple known industrial emissions sources (
SI Appendix, Fig. S2A).
Chloroprene is the largest contributor to cancer risk for two tracts, both of which are in St. John the Baptist Parish, home to the only major production facility of chloroprene in the United States, which accounts for 96% of nationwide emissions according to TRI (year 2022). Chloroprene is a monomer used for the production of polychloroprene rubber (trade name e.g., “Neoprene”), a synthetic rubber material used in a variety of applications (e.g., protective and orthopedic garments, gaskets, and seals, etc.) due to its with high physical flexibility and resistance to water and heat (
30,
43). These tracts and the area around this facility are essentially the only place where cancer risk is attributed to chloroprene in the state, though the magnitude of risk in these tracts is very high. The maximum tract-level risk value for chloroprene is 187-in-one million, which is the highest tract-level value for any chemical in the BR-NO corridor in AirToxScreen estimates. Chloroprene is linked to a wide variety of cancers, including leukemia and those affecting the liver and lungs (
30).
Benzidine [IUR
6.7
10
−2g m
−3 (
44)] is the largest contributor to risk for one tract (tract ID: 22051027501), which is in Jefferson Parish. The benzidine-attributed risk in this tract is only slightly higher than that of formaldehyde (25- vs. 21-in-one million, respectively). This level of risk corresponds to very low concentration (
1 ppt), due to the very high IUR value for benzidine.
Fig. 1C shows that the spatial variation of chloroprene- and ethylene oxide-attributable cancer risk are high relative to formaldehyde, which is comparatively invariable across the corridor. We group these chemicals into categories in
Fig. 1C, where formaldehyde, chloroprene, and ethylene oxide each comprise their own category, while all particulate metals are listed in the “Metals” category and all other VOCs are in the “Other VOCs” category. We see that tracts with the highest estimated total cancer risk have varying, but substantial, contributions from both chloroprene and ethylene oxide, or from ethylene oxide alone.
Based on AirToxScreen estimates, formaldehyde, ethylene oxide, and chloroprene together account for the large majority of total cancer risk in the BR-NO corridor, especially for those tracts with the highest levels of total risk. For the 30 tracts with the highest risk estimates (
Fig. 1C), the relative contribution of the sum of formaldehyde, chloroprene, and ethylene oxide to the total ranged from 41 to 96%, with a median value of 84%. The tract with only 41% of the risk being attributed to these three chemicals is an outlier, as it has a much-larger risk attributed to Other VOCs compared to the other tracts. This is the tract discussed above (tract ID: 22051027501), where benzidine has the single highest risk estimate and also makes up the large majority of risk from the Other VOCs category.
Conversely, those tracts in the corridor with the lowest cancer risk estimates (e.g., the lowest 30 tracts,
Fig. 1C) have very little risk attributed to ethylene oxide relative to formaldehyde and essentially none attributed to chloroprene. Formaldehyde dominates the total risk estimate, contributing on average 71% to the total risk. Formaldehyde-attributed cancer risk across these low-risk tracts, however, is only 19% lower on average than it is for the high-risk tract. Again, this owes to relatively uniform regional background concentrations from photochemical production.
The median relative contribution to total cancer risk for all other VOCs across the corridor was 15% and for metals was 0.6%. All individual risk estimates for Other VOCs in AirToxScreen were below 5-in-one million level of risk, and all individual risk estimates for metals were below 1-in-one million level of risk. The relative contribution for these two categories is roughly similar between the high- and low-risk tracts, with the exception of a few tracts with high relative contributions from Other VOCs, as discussed above. While this paper does not contain any metals measurements, nor the comprehensive set of VOC HAPS measurements treated by AirToxScreen, we will show in the next section that AirToxScreen modeled concentrations can dramatically underestimate actual measured concentrations, which is an argument for targeted sampling even for species and/or areas not highlighted as high-risk in AirToxScreen.
Cancer Risk Estimates from Measurements and Comparison with AirToxScreen.
Using advanced instrumentation aboard a mobile laboratory during February 2023, we measured a large suite of HAPs in the heart of the BR-NO industrial corridor. Seventeen of these VOC HAPs are known carcinogens, and we use their measured concentrations to present the associated cancer risks here. In a previously published mobile monitoring study, we aggregated ethylene oxide measurements into representative concentrations for 500 m
500 m grid cells and at the census tract level (
39). We apply similar spatiotemporal aggregation methods here to the rest of the carcinogenic air pollutants that we measured, and use these representative concentrations as the basis for our cancer risk estimates, assuming that they represent lifetime exposure concentrations.
Table 1 contains the full list of measured compounds and their IUR values used for risk calculations. Details of our spatiotemporal aggregation methods are provided in
Materials and Methods.
Fig. 2 illustrates the magnitude of total and individual estimated cancer risk for the 15 census tracts where we made measurements.
Fig. 2A illustrates how each measured HAP category contributes to the cumulative risk of each tract, with the corresponding estimate from AirToxScreen for comparison. Note that our cancer risk estimates are derived from the 17 HAPs we measured, but we compare to all 72 HAPs with modeled cancer risk in AirToxScreen. We also show the average across all 15 tracts for both measured risk and AirToxScreen. These tracts span parts of four parishes (see
Inset map in
Fig. 2B).
Fig. 2C presents similar information but shows the number of census tracts with measured estimates of risk that are above different thresholds, both for risk from individual HAPs and total.
Like AirToxScreen, we find that formaldehyde, chloroprene, and ethylene oxide drive the large majority of the total cancer risk (range: 62.6 to 96.2%). However, for all tracts but one, our estimates were larger than AirToxScreen (range: 1 to 11.7; median: 5.1). The median tract-level cancer risk from our measurements was 310-in-one million, compared to 60-in-one million for AirToxScreen, and the maximum tract-level risk estimates were 560-in-one million and 250-in-one million, respectively.
The total cancer risk estimates shown in
Fig. 2A are mapped in
Fig. 3. We see that the tracts in the northwest corner of our sampling domain (parts of Iberville and Ascension parishes) had the highest total cancer risk. The middle of our sampling domain had total cancer risk values on the lower side of this distribution (parts of St. James Parish), and then the southeastern part of the domain had higher values as well (parts of St. James and St. John the Baptist parishes) though not as high as the northwestern-most tracts.
SI Appendix, Fig. S5 contains a map with both our tract-level total risk estimates and those for 500 m
500 m grid cells, which illustrates that there can be substantial within-tract variability of total risk as well, with some locations having substantially higher risk levels than the tract average. AirToxScreen estimates had a different spatial pattern compared to our tract-level estimates. For AirToxScreen, the highest-risk tracts were in the southeastern part of the domain (St. John the Baptist parish) and have elevated chloroprene and/or ethylene oxide concentrations relative to the other tracts, as discussed above.
Formaldehyde and Other VOCs contribute roughly equally to total cancer risk, with mean values of 23.9- and 39.4-in-one million, respectively. These estimates represent an average of 8.4% (range: 4.1 to 14.3%) and 10.5% (range: 3.8 to 37.4%) of the total cancer risk we estimate across all tracts. For all tracts, Other VOCs had larger measured risk estimates compared to AirToxScreen, and formaldehyde-associated risk between measurements and AirToxScreen was very similar.
There was little spatial variability in the risk attributed to formaldehyde from our measurements, which is similar to AirToxScreen. The relative SD (RSD), which is the SD divided by mean, in estimated cancer risk from measurements and AirToxScreen were 0.074 and 0.058, respectively. As discussed more fully in Materials and Methods, we applied a “seasonal correction factor” to our formaldehyde measurements in an effort to better represent an annual average, as formaldehyde concentrations exhibit strong seasonality (due to photochemical production) and our measurements were taken during winter. All AirToxScreen concentration estimates are annual averages.
Spatial variability was higher for Other VOCs compared to formaldehyde. The RSD in Other VOCs in AirToxScreen risk estimates was 0.2, while the RSD from measurements was considerably higher (1.273). We measured a few outlier tracts (e.g., “B,” “C,” and “D” in
Fig. 2B) that had high risk; RSD was 0.361 with these tracts removed from the calculation. There is a very high degree of variability within the Other VOCs category as to which specific chemicals drive risk across these tracts, reflecting the large number of different industrial emissions sources in the area. Risk in the Other VOCs category for these outlier tracts, for example, is driven by high levels of 1,2-dichloroethane and, to a lesser extent, chloroform. For example, the risk attributed to 1,2-dichloroethane in tract “B” corresponds to 174-in-one million. This example is important to highlight because the risk is well above the unacceptable threshold but is not highlighted in AirToxScreen at all as a HAP of concern in the area.
Chloroprene and ethylene oxide cancer risk estimates from our measurements are more spatially heterogenous than the Other VOCs category. The observed spatial variation for chloroprene and ethylene oxide is expected of primary pollutants with major industrial sources in the area, in contrast with VOCs dominated by secondary production (e.g., formaldehyde). The RSD in measured risk estimates for chloroprene and ethylene oxide are 1.317 and 0.505, respectively.
For chloroprene, measured concentrations were high in the vicinity of the major emissions source, and below the instrumental limit of detection elsewhere (
SI Appendix, Fig. S6A). The maximum tract-level cancer risk for chloroprene from our measurements was 68-in-one million, while for AirToxScreen it was 187-in-one million. We estimate that chloroprene contributes 7.9% on average to total cancer risk for those tracts in our sampling domain, but the relative contribution varied widely (range: 0.2 to 36.3%)
There are a mix of tracts where AirToxScreen cancer risk estimates for chloroprene are larger than our estimates (e.g., “L”, “N” from
Fig. 2A) and vice versa (e.g., “M”, “O”). However, the broad spatial pattern of chloroprene risk across the domain for the two estimates is similar, with a high-to-low gradient of values moving west away from the single major point source located in the easternmost part of the domain (see
SI Appendix, Fig. S6 for grid cell and tract-level maps of chloroprene concentration and risk). The magnitudes of the two risk estimates were relatively similar for each tract (
SI Appendix, Fig. S4), and the discrepancies can be explained by inaccuracy in underlying emissions estimates and/or potential spatiotemporal sampling bias.
Ethylene oxide, while spatially heterogenous like chloroprene and unlike formaldehyde, represents a much larger fraction of the total cancer risk we estimated compared to AirToxScreen. The maximum tract-level cancer risk for ethylene oxide from our measurements was 515-in-one million compared to 70-in-one million from AirToxScreen. There are multiple ethylene oxide hot spots that we identified with our mobile lab (e.g., in both the northwest and southeast portions of the domain), owing to the much-larger (compared to chloroprene) number of emissions sources in the vicinity. 12 of 14 TRI-listed ethylene oxide emitters in the state are in the BR-NO industrial corridor, and 10 of those are within 11 km of our sampling route.
AirToxScreen identifies ethylene oxide as an important, but not dominant, chemical driving cancer risk in the region. Our measurements, on the other hand, identify ethylene oxide as contributing a large majority (average: 73.2%, range: 39 to 91.7%) of the total risk in the area. Some of this discrepancy is due to an ethylene oxide background concentration of zero in AirToxScreen, while our measurements suggest an area background concentration that corresponds to roughly 200-in-one million cancer risk, which we discuss further in the next section. With the exception of one tract (“N,” which has a high contribution from chloroprene), the risk we estimate from ethylene oxide alone is larger than the total risk estimate from AirToxScreen for all other tracts.
The table in
Fig. 2C facilitates comparison of the spatial extent of impacts between different HAPs. Ethylene oxide is almost uniformly above 100-in-one million levels of risk, while formaldehyde is uniformly between 10- and 50-in-one million levels of risk. 1,2-dichloroethane is the only other VOC we measured which had risk above 100-in-one million anywhere, which it did for one tract. Additionally, three and four tracts for 1,2-dichloroethane and chloroform, respectively, were between 10- and 50-in-one million levels of risk, highlighting that other VOCs do contribute meaningfully to cumulative risk. A number of the chemicals we measured (the bottom seven entries in this table), had risk levels less than 1-in-one million across all tracts in the domain.