New Study: Airborne radioactivity increases downwind of frac’ing; Particles released could damage health of residents nearby. Marco Kaltofen: “This investigation backs up its big conclusions with big data.” Petros Koutrakis, Harvard TH Chan School of Public Health in Boston, led the study: “If you asked me to go and live downwind [of fracking sites], I would not go.”

Dangerous radiation levels from fracking; Dr David Lowry on the risk of the radioactive gas radon being pumped into homes Oct 15, 2020, Letters, The Guardian

I was very interested to read your report (13 October) on recent research by Harvard University scientists on radiation risks from fracking.

I have raised this concern in several letters to the Guardian over the past seven years. Indeed, seven years ago this month, Public Health England said in a review of potential risks that “there is … the potential for radon gas to be present in natural gas extracted from UK shale”.

Eight years ago, Dr Marvin Resnikoff, of Radioactive Waste Management Associates, estimated that radon levels from the Marcellus gas field in the eastern US were up to 70 times the average, and suggested that the radiation from some shale gas deposits was as much as 30 times as high as natural background levels.

Hence, there is undoubtedly a risk of radon gas being pumped into citizens’ homes as part of the shale gas stream. Unless the gas is stored for up to a month to allow the radon’s radioactivity to naturally reduce, this is potentially very dangerous.

Unconventional oil and gas development and ambient particle radioactivity by Longxiang LiAnnelise J. BlombergJohn D. SpenglerBrent A. CoullJoel D. Schwartz & Petros Koutrakis, Oct 13, 2020, Nature Communications volume 11, 5002 (2020)


Unconventional oil and natural gas development (UOGD) expanded extensively in the United States from the early 2000s. However, the influence of UOGD on the radioactivity of ambient particulate is not well understood. We collected the ambient particle radioactivity (PR) measurements of RadNet, a nationwide environmental radiation monitoring network. We obtained the information of over 1.5 million wells from the Enverus database. We investigated the association between the upwind UOGD well count and the downwind gross-beta radiation with adjustment for environmental factors governing the natural emission and transport of radioactivity. Our statistical analysis found that an additional 100 upwind UOGD wells within 20 km is associated with an increase of 0.024 mBq/m3 (95% confidence interval [CI], 0.020, 0.028 mBq/m3) in the gross-beta particle radiation downwind.

Living much more frac’d than that by Encana/Ovintiv in Rosebud Alberta Canada:

(Black dots show deeper frac’d wells and older conventional wells, most by Encana/Ovintiv)

Based on the published health analysis of PR, the widespread UOGD could induce adverse health effects to residents living close to UOGD by elevating PR.


The extraction of crude oil and natural gas from the low-permeability unconventional geological accumulating formation (known as unconventional oil and natural gas development [UOGD]) expanded extensively over the past decade. As of 2017, over 120,000 onshore UOGD wells had been drilled via a practice involving directional drilling combined with multistage high-volume hydraulic fracturing (fracking)1. Meanwhile, numerous controversies have arisen, partially due to the potential harmful impacts on the local environment2,3,4,5,6,7, and on the health of nearby residents8,9,10,11.

Naturally occurring radioactive material (NORM) is a common by-product in Oil and Gas (O&G) production industry. The concentration of Uranium-238 in sedimentary formation rich in organic matter, such as black shale, is significantly higher than the background level of in earth’s crust due to the natural attenuation process12,13. Before widespread UOGD, studies had detected above-background levels of Radium-226, a decay product of U-238, in the wastes of conventional oil and natural gas development (COGD)14,15. Regarding UOGD, enhanced levels of U-238 and Ra-226 have recently been detected in the produced water from unconventional hydrocarbon reservoirs16,17, in the drill cuttings from the lateral drilling within the unconventional formation18,19,20, in the impoundment sediments21, in the soil of brine spill accident scene22, and in the stream sediments near discharging sites23. Two studies in the Marcellus shale region found a positive association between UOGD activities and indoor levels of Radon-222, a gaseous decay product of Ra-22624,25.

However, the influence of UOGD on the radioactivity of ambient particles (referred to as particle radioactivity [PR]) is not well understood. The particle-bound progeny of Radon-222 (referred to as radon in this study) contribute to the majority of PR26,27. Radon firstly decays into a chain of short-lived particle-reactive progeny. These short-lived radionuclides quickly react with the water molecules and atmospheric gases passing by, form ultrafine clusters and finally attach to airborne particles28,29,30. The short-lived progeny on the ambient particles then decay into two long-lived progeny, Lead-210 and Polonium-210, which respectively account for most of the beta- and alpha-radiation emitted by the particulate26,31. UOGD could influence local PR level by increasing the emission rate of radon. There is an increasing interest in the health effects of PR because the particle-bound Lead-210 and Polonium-210 tend to be deposited on the bronchial epithelium and expose neighboring cells to high-energy alpha particles that induce the carcinogenesis process32,33. Short-term exposure to PR has been associated with adverse health outcomes, including a decrease in lung function34, an increase in blood pressure35, and increased levels in biomarkers of inflammation36,37.

In our study, we investigate the likely impact of UOGD on PR. Our statistical analysis demonstrates that upwind UOGD activities could significantly elevate the PR level in downwind communities. UOGD has a larger impact on PR, compared to COGD. The impact of UOGD on PR decreases gradually along with an increasing downwind distance. The results of our study contribute to the currently limited knowledge regarding the influence of UOGD on PR.


In this study, we analyzed 320,796 PR measurements carried out at 157 RadNet monitors across the continental United States from 2001 to 2017 (Fig. 1). We categorized these monitors into O&G RadNet monitors and other RadNet monitors, based on the existence of O&G extraction activities within 50 km. As summarized in Table 1, the national average PR level was 0.35 mBq/m3, with an interquartile range (IQR) from 0.22 mBq/m3 to 0.43 mBq/m3. O&G RadNet monitors had a higher average PR (0.39 mBq/m3; IQR: 0.26,0.47 mBq/m3), compared to the average PR of other O&G RadNet monitors (0.33 mBq/m3; IQR: 0.20,0.41 mBq/m3). Concerning PR emission-related environmental factors, O&G RadNet monitors had a higher ground surface U-238 level, and a higher percentage of air mass originated from the continent. For PR movement-dependent factors, O&G RadNet monitors had higher PM2.5 concentrations, higher planetary boundary layer height (PBLH), higher wind velocity, and lower relative humidity.

Fig. 1: The location of RadNet monitors and the UOGD wells (completed by 2017) in the continental U.S.

Table 1 Descriptive statistics of the ambient particle radioactivity and its environmental predictors.Full size table

After excluding wells without production records, there were 1,574,602 completed O&G wells by the end of 2017. Out of these O&G wells, 152,904 (9.7%) were UOGD wells, and 1,421,698 (91.3%) were COGD wells (Supplementary Fig. 1). Out of the UOGD wells, 4611 (3.0%) were within 20 km of RadNet monitors and 28,016 (18.3%) were within the 50 km buffer. UOGD expanded rapidly in all three subregions: Marcellus-Utica subregion, Permian-Haynesville subregion and Bakken-Niobrara subregion, during the study period (Supplementary Fig. 2, Supplementary Table 1). Fort Worth, Texas, had the highest average upwind UOGD count (mean 586, IQR: 504,661) within 20 km in 2017.

According to our linear mixed effect (LME) model, there was a statistically significant association between the downwind PR and the upwind UOGD activity. With adjustment for environmental factors regarding the natural emission and movement of PR, an additional 100 upwind UOGD wells within 20 km was associated with a 0.024 mBq/m3 increase in the level of PR (95% CI: 0.020, 0.028 mBq/m3) for a wind velocity of 1 m/s (Supplementary Table 2). Under the same wind condition, an additional 100 upwind COGD wells within 20 km was associated with a 0.004 mBq/m3 increase in PR (95% CI: 0.003, 0.004 mB1/m3). Regarding the magnitude of the impact, UOGD and COGD could elevate the PR level by up to 0.13 mBq/m3 and 0.029 mBq/m3, respectively. We also found significant negative interactions between wind velocity and the counts of both UOGD wells and COGD wells located upwind of RadNet monitors, suggesting lower influence when the wind is strong (Table 2).Table 2 The associations between PR and other environmental factors.Full size table

Furthermore, we found that PR is significantly associated with upwind UOGD at every buffer distance of our study (Fig. 2). The influence of an additional 100 UOGD wells decreased gradually as the buffer radius increased from 20 to 50 km. An additional 100 UOGD wells within 50 km was associated with a 0.002 mBq/m3 increase in PR (95% CI: 0.002, 0.003 mBq/m3). Meanwhile, the association between the upwind COGD well count and PR was not statistically significant when the buffer radius is >20 km.

Fig. 2: The association between upwind O&G production activities and downwind level of PR.

In the negative control analysis, we found PR level is also statistically associated with the number of UOGD wells within 20 km downwind of a RadNet monitor. However, the increase of PR associated with an additional 100 downwind UOGD wells within 20 km was 0.021 mBq/m3 (95% CI: 0.017, 0.024 mBq/m3), smaller than the impact of an additional 100 upwind UOGD well (0.024 mBq/m3, 95% CI: 0.020, 0.028 mBq/m3). We found that the per-well influence of upwind UOGD on downwind PR decreased gradually during our study period (Supplementary Note 2, Supplementary Table 3). After evaluating the sensitivity of our results, we argued that our results are not sensitive to the angle of the buffers (Supplementary Note 3, Supplementary Fig. 3). The upwind-downwind difference in the influence of UOGD wells is more significant when the angle of the circular sections is 60° compared to the quadrant buffer used in the primary analysis (Supplementary Table 4). We also found that our results are robust to omitting one of the RadNet monitors in the analysis (Supplementary Fig. 4).

We found that the PR level is statistically associated with the emission- and movement-dependent environmental covariates (Table 2). There were significant positive correlations between PR and the ground surface concentration of U-238, the proportion of continent-sourced air mass, the number of sunspots, PM2.5 concentration, and the inverse of PBLH. Meanwhile, PR is negatively associated with latitude, relative humidity, and soil moisture.

In our subregional analysis, we found significant heterogeneity among the three subregions regarding the influence of UOGD on PR (Table 3). In the Marcellus-Utica subregion, we found no evidence of a statistically significant association between upwind UOGD wells and downwind PR for any buffer distances investigated. Meanwhile, the association was significant for each buffer distance in the Permian-Haynesville subregion. In the Bakken-Niobrara subregion, the association was not statistically significant when the buffer radius is smaller than 30 km. However, when the buffer radius is 20 km, the impact of additional 100 UOGD wells in Marcellus-Utica subregion (0.180 mBq/m3, 95% CI: −0.031 mBq/m3, 0.390 mBq/m3), though not significant, was greater than the impact in the other two subregions (Table 3).Table 3 The associations between PR and upwind UOGD well count in three subregions of our study extent.Full size table


In this study, we analyzed the radioactivity of airborne particles collected at 157 RadNet monitors across the continental United States from 2001 to 2017 (Fig. 1). To characterize upwind UOGD activities, we used the position and production records of 152,904 UOGD wells and counted the daily number of upwind UOGD wells. Our results added to the limited literature by evaluating the influence of UOGD on the radioactivity of ambient particles.

These associations suggested the existence of some pathways by which UOGD activities could release NORM into the atmospheric environment. Likely mechanisms include the fugitive release of natural gas, which contains a higher-than-background level of radon at wellheads, compressor stations, pipelines, and other associated facilities38,39,40; the management, storage, discharge and disposal of flow-back and produced water which is rich in NORMs16,41,42,43; the accidental spill or beneficial use of produced water in nearby communities22; the handling, transport, management, and disposal of radioactive drill cuttings18,19. The results of our negative control analysis suggested the potential transport mechanism of PR independent of atmospheric movement. This association could be explained by the increased off-site radon emission, or by the nearly isostropic dispersion of radon released on site under low-wind condition. To distinguish the contributions of these surface activities, more continuous measurements of PR, especially for some specific radionuclides, are needed at a finer spatiotemporal resolution.

Our results showed a remarkable distinction between the impacts of UOGD and COGD on PR. UOGD-specific processes, such as hydraulic fracturing and directional drilling, could potentially explain the larger associated impacts. The high-volume hydraulic fracturing process produced large volumes of flow-back water and drilling mud, which are subsequently stored in the temporary reserve pit adjacent to the drilling site. Most UOGD production states allow the operator to close the reserve pit within up to one year after completing the drilling44. This practice potentially enables the NORMs in the produced water to decay into radon above the ground surface and release the radon into the ambient environment. The lateral drilling process produces large volumes of drill cuttings from the unconventional accumulating formation, whose levels of NORMs are higher than those produced during the vertical drilling stage. These drill cuttings are currently not considered hazardous wastes by U.S. EPA. The practice of beneficial use of drill cuttings and land treatment could potentially release radon into the ambient environment45.

Our subregional analysis demonstrates remarkable heterogeneity in the estimated influences of UOGD in the three subregions. Due to a lack of monitors with UOGD wells nearby, the Marcellus-Utica subregional model did not have enough power to detect statistically significant associations. For a buffer radius of 20 km, the estimated influence of an additional 100 upwind UOGD wells is 17.96 × 10−2 mBq/m3 (95% CI: −3.06, 38.99 × 10−2 mBq/m3) in the Marcellus-Utica subregion, approximately seven times the estimated effects of a nationwide model. The difference is likely caused by the relatively few UOGD wells near the RadNet monitors in the Marcellus-Utica subregion (Supplementary Fig. 2). In the Bakken-Niobrara subregion, only two RadNet monitors (Casper, WY, and Navajo Lake, NM) have active UOGD wells around when the buffer radius is smaller than 30 km. When we enlarged the radius to 30 km, two additional RadNet monitors (Denver, CO, and Grand Junction, CO) had UOGD within the buffer, enabling us to identify the significant association (Supplementary Fig. 3).

Our results show a monotonic declining impact of O&G wells on PR as the buffer radius increases from 20 to 50 km (Fig. 2). The attenuation of radon after being emanated can explain this pattern. The trend indicates a more significant influence on the PR level of communities close to intensive UOGD activities. Limited by the accuracy of RadNet monitor location information, we did not estimate the impacts on a buffer distance smaller than 20 km. To tentatively extrapolate our results to these neighborhoods, we modeled the estimated influences as a power function of the radiuses with a negative exponent (Supplementary Note 5, Supplementary Fig 5). Based on this tentative extrapolation, an additional 100 UOGD wells within 10 km would be associated with an increase of 0.14 mBq/m3. However, the result of this extrapolation should be interpreted cautiously. Monitors closer to UOGD wells are needed to validate this extrapolation.

One strength of our study is the nationwide monitor network of PR. The long-term measurement enables us to compare the current PR level with the baseline PR level in the absence of widespread UOGD. Furthermore, all filters were measured by NAREL using the same protocol, excluding the uncertainties induced by the heterogeneous devices operated by different labs. The other strength is the comprehensive database covering O&G activities. The Enverus database facilitated distinguishing the distinct impacts of UOGD and COGD on PR. In addition, we obtained diverse environmental covariates related to the natural emissions and transport of PR. Adjustment for these factors allowed us to draw conclusions explicitly related to the impacts of O&G development by explaining the natural variation of PR. The associations between PR and these environmental factors, as summarized in Table 2, are in agreement with the findings of previous studies26,46.

One limitation of this study is that we only associated PR with the existence of completed O&G wells. Other construction-dependent factors may also influence the emission rate. However, the O&G wells with a detailed construction record are rare in our database, making it difficult to know the duration of construction, thus limiting our ability to investigate the construction-dependent association. Another limitation of our study is the simplification of the particle transport process. Our circular-sectional buffers are designed based on the Gaussian Dispersion Model. This calculation assumed a steady-state meteorological condition, which is reasonable in this case due to the short downwind distances. However, this computation could be improved by introducing advanced atmospheric dispersion models. Finally, we used drilling type as a proxy to whether an O&G well is UOGD or not. The method inevitably missed some vertical wells uncommonly completed by high-volume hydraulic fracturing.

Our results indicate the significant influence of UOGD on PR, a previously overlooked property of PM2.5. Particulate-bound radon progeny continue releasing ionizing radiation after being inhaled and thus could induce systemic oxidative stress and inflammation, even at the levels observed in this study. Nyhan et al. (2018 and 2019) found that a 0.07 mBq/m3 increase in 28-day average gross-beta radiation is associated with a 2.95 mm-Hg increase in diastolic blood pressure, a 3.94 mm-Hg increase in systolic blood pressure, a 2.41% decrease in forced vital capacity, and a 2.41% decrease in forced expiratory volume in Normative Ageing Study (NAS) population34,35. Blomberg et al.37 reported that a 0.12 mBq/m3 increase in seven-day average gross-beta radiation is associated with a 4.9% increase in C-reactive protein, a 2.8% increase in intercellular adhesion molecule-1, and a 4.3% increase in vascular cell adhesion molecule-1 in the same study popultaion. Jointly with these associations, our results suggest that an increase in PR due to the extensive UOGD may cause adverse health outcomes in nearby communities by elevating PR level (Supplementary Note 6). Further studies, especially those based on PR measurements close to UOGD activities, are needed to validate this exposure pathway.

Our analysis demonstrates that upwind UOGD activities could significantly elevate the PR level in downwind communities. UOGD has a larger impact on PR, compared to COGD. Based on previously published health effect analysis of PR, it is possible that the widespread of UOGD could induce adverse health effects to residents in proximity by elevating the PR.


Ambient PR measurements

We obtained PR measurements carried out by the RadNet monitoring network, which is operated by the U.S. Environmental Protection Agency (EPA). This network measures the background environmental radiation levels in the air, precipitation, and drinking water under both routine and emergency conditions. During the study period from 2001 to 2017, 157 RadNet sites (Fig. 1) reported gross-beta measurements of various time lengths. Most RadNet monitors are located in metropolitan areas for better population coverage. At each site, total suspended particles (TSP) are collected using a high-volume sampler with a 4-inch diameter polyester fiber filter. Samplers are operated continuously for a 3- or 4-day integration. Filters are then sent to the National Analytical Radiation Environmental Laboratory (NAREL) for the measurement of gross-beta radiation47,48. To create quasi-daily values from the 3- or 4-day integrated samples, we assigned the same beta-radiation level to each day of the collection period.

Unconventional oil and gas development data

We obtained position and production information of O&G wells from Enverus (formerly, a third-party data vendor used by the U.S. Energy Information Administration (EIA) to prepare monthly fossil fuel production and marketing reports. The comprehensive data coverage of Enverus is achieved by compiling the permits, construction logs, and production records from state agencies. Details about this data source were presented in a previous study49. Our dataset includes information for 2,159,858 wells hydraulically fractured from 01/01/1949 to 12/31/2017. We used drilling type information as the primary indicator of whether a well is targeting an unconventional accumulation formation or not50. Specifically, we considered horizontally drilled wells as UOGD wells and vertically drilled wells as COGD wells. Directionally drilled wells and wells without drilling type information were classified into UOGD or COGD wells based on their proximity to nearby UOGD wells and other secondary information using a Random Forest model (Supplementary Note 7).

We used the number of completed wells to characterize the intensity of O&G production activity. We identified a well as completed when the operator received the well from the driller. If this information is unavailable, we used the first production date as a proxy. Considering the transport of airborne particles, we focused on the completed wells positioned upwind of the RadNet monitor. Specifically, we created a circular sectional buffer centered on the daily wind direction with an angle of 90° and a radius of 20 km (Fig. 3). We counted the numbers of UOGD and COGD wells, respectively, within the buffer on a daily basis to detect the different impacts of the two types of wells. We also created a series of circular sectional buffers at distances ranging from 25 to 50 km by 5 km intervals to evaluate the potential dependency of impact on the spatial scale. For security reasons, the exact location of RadNet monitors was not publicly accessible. We used the centroid of each RadNet city as a proxy to the sampling location (Fig. 3). Because of the potential spatial mismatch, we did not run any analysis on a spatial scale smaller than 20 km.

Fig. 3: Methods to calculate the number of UOGD wells positioned upwind of the RadNet monitor in two example days.

Predictors of PR

We collected data on environmental variables related to the emission of PR. To control for the emanation rate of radon from soil, we downloaded ground surface concentration of U-238 at a spatial resolution of 3 km from the United States Geological Survey51. PR is associated with the origin of air masses because the emanation rate of radon from the ocean is two orders lower than that from the continent27. To capture this pattern, we modeled four 72-hour back-trajectories (arrival time 06:00, 12:00, 18:00, and 24:00) of each RadNet monitoring sites using the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model52. The proportion of trajectory over the continent was used as a proxy to the origin of the air mass with 0% indicating a maritime air mass, while 100% indicating a continental air mass. Finally, we collected the monthly number of sunspots observed by the Royal Observatory of Belgium53, as an indicator of the strength of solar activity. This information was used to adjust for the contribution of beta-emitting cosmogenic radionuclides originating from the upper atmosphere.

We also obtained environmental factors influencing the transport of PR. Due to the scavenger effect of aerosol on the short-lived progeny of radon, PM2.5 concentration strongly influences the spatiotemporal distribution of radioactive isotopes in the atmosphere27. We downloaded daily PM2.5 concentrations measured at EPA air quality monitors located within 50 km of RadNet sites and calculated the daily average if more than one measurement was carried out. The behavior of PR is driven by multiple meteorological factors, including wind velocity, relative humidity, PBLH, temperature, and soil moisture26,27. We obtained these variables from the North America Regional Reanalysis (NARR) dataset, with a spatial resolution of 32 km54.

Statistical analysis

We applied LME models to investigate the association between PR and UOGD activities. Dependency between daily PR measurements from the same monitor was controlled for by including monitor-specific random intercepts. We controlled for relevant environmental factors as fixed effects. We also adjusted for long-term temporal trend and seasonality by including polynomial terms based on the calendar year and temperature. We applied the LME as follows:PRi,t=(c0+γi)+c1⋅Numi,t+c2⋅Ui+c3⋅origini,t+c4⋅sunt+c5⋅pmi,t+c6⋅pblh−1i,t+c7⋅rhumi,t+c8⋅soilmi,t+c9⋅veli,t+c10⋅lati+∑p=11,12cp⋅yearp−10t+∑p=13⋯15cp⋅tempp−12i,t+c16⋅sunspotst×lati+c17⋅soilmi,t×U+c18⋅Numi,t×veli,ti,PRi,t=(c0+γi)+c1⋅Numi,t+c2⋅Ui+c3⋅origini,t+c4⋅sunt+c5⋅pmi,t+c6⋅pblhi,t−1+c7⋅rhumi,t+c8⋅soilmi,t+c9⋅veli,t+c10⋅lati+∑p=11,12⁡cp⋅yeartp−10+∑p=13⋯15⁡cp⋅tempi,tp−12+c16⋅sunspotst×lati+c17⋅soilmi,t×U+c18⋅Numi,t×veli,ti,(1)

where PRi,t is PR level of site i on day t; Num represents the number of upwind wells within the circular sectional buffer; c2 to c4 are the coefficients for the emission-dependent variables, including U for the ground surface concentration of U-238, origin for the origin of the air mass and sun for the number of sunspots; c5 to c9 are the coefficients for transport-related environmental factors, including pm for the concentration of PM2.5, pblh−1 for the inverse of HPBL, rhum for relative humidity, soilm for the moisture of soil and vel for wind velocity; c10 is the coefficient for latitude-dependent spatial trend; c11 to c15 are the coefficients for temporal trends terms, represented by the polynomial terms of the calendar year (year) and temperature (temp); c16 to c18 are the coefficients for interactions terms between environmental factors.

In the primary analysis, we associated daily PR levels with the daily number of upwind UOGD wells within 20 km using LME. To investigate the magnitude of UOGD’s impact, we calculated the increase in PR associated with the 95% percentile of upwind UOGD well count. To investigate the dependency of effects on the transport distance, we counted the number of O&G wells within a series of circular sectional buffers at distances ranging from 25 to 50 km by 5 km intervals, and then estimated the effects for each buffer distances.

To investigate the influence of COGD wells, we associated daily PR with the upwind number of COGD wells within the same buffers. As a negative control, we counted the number of downwind UOGD wells within 20 km. We hypothesized that PR is associated with both upwind UOGD and downwind UOGD activities, but the influence should be smaller than that of upwind UOGD. Because of the pronounced changes in hydraulic fracturing practice in our study period50, we explored the variation in the influence by restricting our model to four different sub-periods. To identify the potential regional heterogeneity in the effects of UOGD, we restricted our analysis to three separate subregions named after the shale formations underneath: Marcellus-Utica subregions, Permian-Haynesville subregion, and Bakken-Niobrara subregion (Supplementary Fig. 2). For sensitivity analysis, we first re-evaluated the associations by calculating the upwind UOGD wells within circular sectional buffers using another two central angles (60° and 120°). We then performed a leave-one-out sensitivity analysis to assess whether our estimates were sensitive to the omission of any single RadNet site.

We used the LME methods implemented in lme4 package (version 1.1–21)55 in R (version 3.4.2)56 to fit the models. The significance test was based on confidence intervals instead of p-values. The analysis was conducted on the Cannon cluster, supported by the Research Computing Group at Harvard University, Faculty of Arts and Sciences.

Data availability

All source data for figures are available in the Supplementary Information. Explanatory data sources, including NARR meterological data (, ground surface uranium grid data (, monitor-based PM2.5 measurements (, and number of sun black spot (, are publicly accessible. The UOGD data from is subscription-only, thus cannot be shared publicly. If a token is provided, the raw data can be obtained from Enverus directly with the code we provided (55_Download_DI_Access.R). The aggregated data that support the findings of this study are available from the corresponding author upon reasonable request.

Code availability

All model codes are available at:


  1. 1.U.S. Energy Information Administration (EIA). The Distribution of U.S. Oil and Natural Gas Wells by Production Rate (2019).
  2. 2.Allen, D. T. Atmospheric emissions and air quality impacts from natural gas production and use. Annu. Rev. Chem. Biomol. Eng. 5, 55–75 (2014).CAS PubMed Article PubMed Central Google Scholar 
  3. 3.Cheadle, L. C. et al. Surface ozone in the Colorado northern Front Range and the influence of oil and gas development during FRAPPE/DISCOVER-AQ in summer 2014. Elem. Sci. Anth. 5, 61 (2017).Article Google Scholar 
  4. 4.Hill, E. & Ma, L. Shale gas development and drinking water quality. Am. Econ. Rev. Pap. Proc107, 522–525 (2017).
  5. 5.Olmstead, S. M., Muehlenbachs, L. A., Shih, J. S., Chu, Z. & Krupnick, A. J. Shale gas development impacts on surface water quality in Pennsylvania. Proc. Natl Acad. Sci. USA 110, 4962–4967 (2013).ADS CAS PubMed Article PubMed Central Google Scholar 
  6. 6.Blair, B. D., Brindley, S., Dinkeloo, E., McKenzie, L. M. & Adgate, J. L. Residential noise from nearby oil and gas well construction and drilling. J. Expo. Sci. Environ. Epidemiol. 28, 538–547 (2018).PubMed Article PubMed Central Google Scholar 
  7. 7.Franklin, M., Chau, K., Cushing, L. J. & Johnston, J. E. Characterizing flaring from unconventional oil and gas operations in south Texas using satellite observations. Environ. Sci. Technol. 53, 2220–2228 (2019).ADS CAS PubMed Article PubMed Central Google Scholar 
  8. 8.Casey, J. A. et al. Unconventional natural gas development and birth outcomes in Pennsylvania, USA. Epidemiology 27, 163–72 (2016).PubMed PubMed Central Article Google Scholar 
  9. 9.Rasmussen, S. G. et al. Association between unconventional natural gas development in the marcellus shale and asthma exacerbations. JAMA Intern. Med. 176, 1334 (2016).PubMed PubMed Central Article Google Scholar 
  10. 10.Koehler, K. et al. Exposure assessment using secondary data sources in unconventional natural gas development and health studies. Environ. Sci. Technol. 52, 6061–6069 (2018).ADS CAS PubMed PubMed Central Article Google Scholar 
  11. 11.McKenzie, L. M. et al. Relationships between indicators of cardiovascular disease and intensity of oil and natural gas activity in Northeastern Colorado. Environ. Res. 170, 56–64 (2019).CAS PubMed Article Google Scholar 
  12. 12.Commission, U. S. A. E. Geology and Geochemistry of Uranium in Marine Black Shales A Review Geology and Geochemistry of Uranium in Marine Black Shales A Review. (1961).
  13. 13.Cordeiro, C., Favas, P. J. C., Pratas, J., Sarkar, S. K. & Venkatachalam, P. Uranium accumulation in aquatic macrophytes in an uraniferous region: relevance to natural attenuation. Chemosphere 156, 76–87 (2016).ADS CAS PubMed Article Google Scholar 
  14. 14.Kolb, W. A. & Wojcik, M. Enhanced radioactivity due to natural oil and gas production and related radiological problems. Sci. Total Environ. 45, 77–84 (1985).ADS CAS PubMed Article Google Scholar 
  15. 15.Fisher, R. S. Geologic and geochemical controls on naturally occurring radioactive materials (NORM) in produced water from oil, gas, and geothermal operations. Environ. Geosci. 5, 139–150 (1998).MathSciNet Article Google Scholar 
  16. 16.Torres, L., Yadav, O. P. & Khan, E. Risk assessment of human exposure to Ra-226 in oil produced water from the Bakken Shale. Sci. Total Environ. 626, 867–874 (2018).ADS CAS PubMed Article Google Scholar 
  17. 17.Brown, V. J. Radionuclides in fracking wastewater: managing a toxic blend. Environ. Health Perspect. 122, A50–5 (2014).PubMed PubMed Central Google Scholar 
  18. 18.Eitrheim, E. S., May, D., Forbes, T. Z. & Nelson, A. W. Disequilibrium of naturally occurring radioactive materials (NORM) in drill cuttings from a horizontal drilling operation. Environ. Sci. Technol. Lett. 3, 425–429 (2016).CAS Article Google Scholar 
  19. 19.Pennsylvania Department of Environmental Protection. Technologically enhanced naturally occuring radioactive materials (TENORM) study report. (2016).
  20. 20.Zhang, T., Hammack, R. W. & Vidic, R. D. Fate of radium in marcellus shale flowback water impoundments and assessment of associated health risks. Environ. Sci. Technol. 49, 9347–9354 (2015).ADS CAS PubMed Article PubMed Central Google Scholar 
  21. 21.Rich, A. L. & Crosby, E. C. Analysis of reserve pit sludge from unconventional natural gas hydraulic fracturing and drilling operations for the presence of technologically enhanced naturally occurring radioactive material (TENORM). NEW Solut. 23, 117–135 (2013).PubMed Article PubMed Central Google Scholar 
  22. 22.Lauer, N. E., Harkness, J. S. & Vengosh, A. Brine spills associated with unconventional oil development in North Dakota. Environ. Sci. Technol. 50, 5389–5397 (2016).ADS CAS PubMed Article PubMed Central Google Scholar 
  23. 23.Lauer, N. E., Warner, N. R. & Vengosh, A. Sources of radium accumulation in stream sediments near disposal sites in Pennsylvania: implications for disposal of conventional oil and gas wastewater. Environ. Sci. Technol. 52, 955–962 (2018).ADS CAS PubMed Article PubMed Central Google Scholar 
  24. 24.Casey, J. A. et al. Predictors of indoor radon concentrations in Pennsylvania, 1989–2013. Environ. Health Perspect. 123, 1130–1137 (2015).CAS PubMed PubMed Central Article Google Scholar 
  25. 25.Xu, Y., Sajja, M. & Kumar, A. Impact of the hydraulic fracturing on indoor radon concentrations in Ohio: a multilevel modeling approach. Front. Public Heal. 7, 76 (2019).Article Google Scholar 
  26. 26.Hernández, F., Hernández-Armas, J., Catalán, A., Fernández-Aldecoa, J. C. & Karlsson, L. Gross alpha, gross beta activities and gamma emitting radionuclides composition of airborne particulate samples in an oceanic island. Atmos. Environ. 39, 4057–4066 (2005).ADS Article CAS Google Scholar 
  27. 27.Baskaran, M. Po-210 and Pb-210 as atmospheric tracers and global atmospheric Pb-210 fallout: a review. J. Environ. Radioact. 102, 500–513 (2011).CAS PubMed Article PubMed Central Google Scholar 
  28. 28.Porstendörfer, J. Properties and behaviour of radon and thoron and their decay products in the air. J. Aerosol Sci. 25, 219–263 (1994).ADS Article Google Scholar 
  29. 29.Gründel, M. & Porstendörfer, J. Differences between the activity size distributions of the different natural radionuclide aerosols in outdoor air. Atmos. Environ. 38, 3723–3728 (2004).ADS Article CAS Google Scholar 
  30. 30.Mohery, M., Abdallah, A. M., Al-Amoudi, Z. M. & Baz, S. S. Activity size distribution of some natural radionuclides. Radiat. Prot. Dosim. 158, 435–441 (2014).CAS Article Google Scholar 
  31. 31.Cabello, M., Dueñas, C., Liger, E., Gordo, E. & Cañete, S. Variables influencing the gross alpha and gross beta activities in airborne particulate samples in Málaga, Spain. J. Radioanal. Nucl. Chem. 315, 299–307 (2018).CAS Article Google Scholar 
  32. 32.Darby, S. et al. Radon in homes and risk of lung cancer: Collaborative analysis of individual data from 13 European case-control studies. Br. Med. J. 330, 223–226 (2005).CAS Article Google Scholar 
  33. 33.Duan, P. et al. Nonlinear dose-response relationship between radon exposure and the risk of lung cancer: evidence from a meta-analysis of published observational studies. Eur. J. Cancer Prev. 24, 267–277 (2015).CAS PubMed Article Google Scholar 
  34. 34.Nyhan, M. M. et al. Associations between ambient particle radioactivity and lung function. Environ. Int. 130, 104795 (2019).CAS PubMed PubMed Central Article Google Scholar 
  35. 35.Nyhan, M. M. et al. Associations between Ambient Particle Radioactivity and Blood Pressure: The NAS (Normative Aging Study). J. Am. Heart Assoc. 7, e008245 (2018).PubMed PubMed Central Article Google Scholar 
  36. 36.Li, W. et al. Recent exposure to particle radioactivity and biomarkers of oxidative stress and inflammation: the Framingham Heart Study. Environ. Int. 121, 1210–1216 (2018).CAS PubMed PubMed Central Article Google Scholar 
  37. 37.Blomberg, A. et al. The role of ambient particle radioactivity in inflammation and endothelial function in an elderly cohort. Epidemiology 31, 1 (2020).
  38. 38.Alvarez, R. A., Pacala, S. W., Winebrake, J. J., Chameides, W. L. & Hamburg, S. P. Greater focus needed on methane leakage from natural gas infrastructure. Proc. Natl Acad. Sci. USA 109, 6435–6440 (2012).ADS CAS PubMed Article Google Scholar 
  39. 39.Howarth, R. W., Santoro, R. & Ingraffea, A. Methane and the greenhouse-gas footprint of natural gas from shale formations. Climatic Change 106, 679–690 (2011).ADS CAS Article Google Scholar 
  40. 40.Rowan, E. L. Radon-222 Content of Natural Gas Samples from Upper and Middle Devonian Sandstone and Shale Reservoirs in Pennsylvania: Preliminary Data Open-File USGS National Produced Waters Geochemical Database View project. (2014).
  41. 41.Rowan, E. L., Engle, M. A., Kirby, C. S. & Kraemer, T. F. Radium Content of Oil-and Gas-Field Produced Waters in the Northern Appalachian Basin (USA): Summary and Discussion of Data (2011).
  42. 42.Brown, V. J. Radionuclides in Fracking Wastewater: Managing a ToxicBlend. Environ. Health Perspect. 122, A50 (2014).PubMed PubMed Central Google Scholar 
  43. 43.Hill, L. L., Czolowski, E. D., DiGiulio, D. & Shonkoff, S. B. C. Temporal and spatial trends of conventional and unconventional oil and gas waste management in Pennsylvania, 1991–2017. Sci. Total Environ. 674, 623–636 (2019).ADS CAS PubMed Article PubMed Central Google Scholar 
  44. 44.Ramirez, P. Earthen pits excavated adjacent to drilling rigs are commonly used for the disposal of drilling muds and fluids in natural gas or oil fieldsU.S. Fish and Wildlife Service (2009).
  45. 45.Association of State and Territorial Solid Waste Management Officials. Beneficial Use of Drill Cuttings, Produced Water and Fluids in the U.S. (2017).
  46. 46.Persson, B. R. R. & Holm, E. Polonium-210 and lead-210 in the terrestrial environment: a historical review. J. Environ. Radioact. 102, 420–429 (2011).CAS PubMed Article PubMed Central Google Scholar 
  47. 47.Fraass, R. In Nuclear Terrorism and National Preparedness (eds. Apikyan, S. & Diamond, D.) 117–123 (Springer, Dordrecht, 2015).
  48. 48.U.S. Environmental Protection Agency (EPA). Expansion and Upgrade of the RadNet Air Monitoring Network: Conceptual Plan and Implementation Process. Office of Radiation and Indoor Air. Vol. I, 2012. (2012).
  49. 49.Czolowski, E. D., Santoro, R. L., Srebotnjak, T. & Shonkoff, S. B. C. Toward consistent methodology to quantify populations in proximity to oil and gas development: a national spatial analysis and review. Environ. Health Perspect. 125, 086004 (2017).PubMed PubMed Central Article Google Scholar 
  50. 50.U.S. Enverionmental Preotection Agency (EPA). Hydraulic Fracturing For Oil And Gas: Impacts From The Hydraulic Fracturing Water Cycle On Drinking Water Resources In The United States (Final Report). (2016).
  51. 51.Duval, J. S., Carson, J. M., Holman, P. B. & Darnley, A. G. Terrestrial radioactivity and gamma-ray exposure in the United States and Canada: U.S. Geological Survey Open-File Report 2005-1413. (2005).
  52. 52.Stein, A. F. et al. NOAA’s HYSPLIT Atmospheric Transport and Dispersion Modeling System. Bull. Am. Meteorol. Soc. 96, 2059–2077 (2015).ADS Article Google Scholar 
  53. 53.Royal Observatory of Belgium. Sunspot data from the World Data Center SILSO. (2019).
  54. 54.Mesinger, F. et al. North American regional reanalysis. Bull. Am. Meteorol. Soc. 87, 343–360 (2006).ADS Article Google Scholar 
  55. 55.Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4J. Stat. Softw. 67, 1–48 (2015).Article Google Scholar 
  56. 56.R Core Team. R: A Language and Environment for Statistical Computing. (2017).

Download references


This publication is made possible by U.S. EPA grant RD-835872. Its contents are solely the responsibility of the grantee and do not necessarily represent the official view of the U.S. EPA. Further, U.S. EPA does not endorse the purchase of any commercial products or services mentioned in the publication. We thank Drew Michanowicz and Jonathan Buonocore for their contribution to the work.

Author information


  1. Department of Environmental Health, Harvard T.H Chan School of Public Health, Boston, MA, 02114, USALongxiang Li, Annelise J. Blomberg, John D. Spengler, Joel D. Schwartz & Petros Koutrakis
  2. Department of Biostatistics, Harvard T.H Chan School of Public Health, Boston, MA, 02114, USABrent A. Coull


P.K. initiated the study; L.L. synthesized data and performed research; L.L, B.A.C, and P.K. developed the model; and L.L., A.J.B., and P.K. wrote the manuscript. L.L., P.K., A.J.B., B.A.C., John D. Spengler, Joel D. Schwartz together interpreted the results.

Corresponding author

Correspondence to Petros Koutrakis.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Communications thanks Seth Shonkoff and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information

Peer Review File

Airborne radioactivity increases downwind of fracking, study finds, Particles released by drilling could damage the health of nearby residents, say scientists by Damian Carrington Environment editor, 13 Oct 2020, The Guardian

The radioactivity of airborne particles increases significantly downwind of fracking sites in the US, a study has found.

The Harvard scientists said this could damage the health of people living in nearby communities and that further research was needed to understand how to stop the release of the radioactive elements from under the ground.

The radioactivity rose by 40% compared with the background level in the most affected sites.

The increase will be higher for people living closer than 20km to the fracking sites, which was the closest distance that could be assessed with the available data.

The scientists used data collected from 157 radiation-monitoring stations across the US between 2001 and 2017. The stations were built during the cold war when nuclear war was a threat. They compared data with the position and production records of 120,000 fracking wells.

“Our results suggest that an increase in particle radioactivity due to the extensive [fracking development] may cause adverse health outcomes in nearby communities,” the team concluded.

Petros Koutrakis at the Harvard TH Chan School of Public Health in Boston, who led the study, said:

“If you asked me to go and live downwind [of fracking sites], I would not go. People should not go crazy, but I think it’s a significant risk that needs to be addressed.”

Previous work has shown that chemicals released during fracking could pose a health risk to children and the process has contaminated groundwater in some places. Fracking is also an issue in the forthcoming presidential election, particularly in swing states such as Pennsylvania. Donald Trump has falsely claimed Joe Biden will ban fracking but the Democratic presidential candidate is largely supportive of fracking and only backs a ban on federal lands and offshore.

The new research, published in the journal Nature Communications, examined the increases in the radioactivity of airborne particles when there were operational fracking wells within 20km upwind of a location. With 100 wells upwind the average rise in radioactivity was 7%, but some places had nearly 600 wells upwind.

The team took into account other factors, including weather and sunspot activity. Cosmic rays produced by sunspots increase the levels of particle radioactivity.

The researchers found that fracking resulted in a far bigger increase in particle radioactivity than conventional oil and gas operations.

This is because the initial source of the radioactivity is a uranium isotope in the rocks. Tapping a conventional oil and gas reservoir barely disturbs the rock. But in the shale formations targeted by frackers the oil and gas is trapped within the rock, which is blasted apart with high-pressure water and releases the uranium.

The uranium isotope decays to the gas radon, which itself decays to ultrafine radioactive particles containing polonium and lead. These are thought to become attached to particles already in the air and are then carried by the wind.

“The polonium isotopes are the ones which are very toxic,” said Koutrakis. The element was used as a poison to kill the former Russian spy Alexander Litvinenko in London in 2006. Koutrakis said previous studies have shown that increases in particle radioactivity of the scale seen in his work can have harmful effects on people.

The biggest increases were seen in Marcellus and Utica shale fields below Pennsylvania, Ohio and West Virginia, as these rocks contain more uranium than in other fracking regions, such as Texas and North and South Dakota.

This investigation backs up its big conclusions with big data, making the conclusions convincing,” said Marco Kaltofen, at the Worcester Polytechnic Institute in Massachusetts, who was not part of the research team.

Wastewater from fracking is also known to have elevated radioactivity, but Kaltofen said: “We shouldn’t ask how airborne radioactive risk from fracking compares to the waterborne risk. We should ask if it’s a good idea to add radioactive particulates to either air or water. This study suggests that the answer is no.”

Koutrakis said it is not known if the radon mostly leaks during initial drilling of the wells, as leaks of methane during production, or from the wastewater produced and stored at the sites. Finding this out could allow operations to be changed to stop the leaks, he said.

“Nothing is more important than the health and safety of our workers, the local environment and the communities where we live, operate and raise families, and we take each report of safety or health issues related to our operations very seriously,” Bullshit! If his yammering was true, frac’ers would not frac – anywhere – never mind in and under parks, communities, homes, schools, care facilities, hospitals, etc said Reid Porter, spokesperson for the American Petroleum Institute. “US natural gas and oil companies meet or exceed strict federal and state regulations and also undergo regular inspections to ensure that all materials are managed, stored, transported, and disposed of safely.” Pffffft. Standard boring oil & gas propaganda. Canadian companies and their enablers (regulators and politicians) spit out similar lies, Encana/now American Ovintiv included. The photo below shows how Encana rid itself of its toxic waste in my community – on foodland – with Alberta’s lying bullying law-violating regulators nowhere to be seen.

Refer also to:

Oct 15, 2020: Final EHP Webinar in Shale Gas Development and Cancer Series: Radioactivity in frac fluid and natural gas by David O. Carpenter, MD

2020 06 15: Cabot Oil and Gas gets frac’d! Pennsylvania AG Shapiro files 15 criminal charges, including 9 felonies, against Cabot for failing to fix gas wells that leaked methane over 9 square miles into 18 Dimock residential water supplies. “The Grand Jury presentments prove that Cabot took shortcuts that broke the law, and damaged our environment — harming our water supply and public health,”

2020 05 09: Gas stoves are a major source of indoor air pollution, exposing millions of people to “pollution in their homes that would be illegal outdoors.” No wonder industry fights to prevent regulation of stoves fueled with their foul products!

2018 01 17: UK, New Study: Fracking is one of the most harmful forms of energy production

2015 04 08: “It looks like fracking has unearthed an unbargained for and serious cancer risk in peoples’ homes.” John Hopkins study links radon levels in Pennsylvania homes to fracking: “These findings worry us”

2013 12 26: Yet Another Peer Reviewed Study Warning about the Health Hazards of Natural Gas: Home gas ranges produce toxic gases says Lawrence Berkeley Lab Study

2001 08 01: The Health Hazards of Natural Gas

This entry was posted in Global Frac News. Bookmark the permalink.