Chlorine activation and enhanced ozone depletion induced by wildfire aerosol.

Remarkable perturbations in the stratospheric abundances of chlorine species and ozone were observed over Southern Hemisphere mid-latitudes following the 2020 Australian wildfires1,2. These changes in atmospheric chemical composition suggest that wildfire aerosols affect stratospheric chlorine and ozone depletion chemistry. Here we propose that wildfire aerosol containing a mixture of oxidized organics and sulfate3-7 increases hydrochloric acid solubility8-11 and associated heterogeneous reaction rates, activating reactive chlorine species and enhancing ozone loss rates at relatively warm stratospheric temperatures. We test our hypothesis by comparing atmospheric observations to model simulations that include the proposed mechanism. Modelled changes in 2020 hydrochloric acid, chlorine nitrate and hypochlorous acid abundances are in good agreement with observations1,2. Our results indicate that wildfire aerosol chemistry, although not accounting for the record duration of the 2020 Antarctic ozone hole, does yield an increase in its area and a 3-5% depletion of southern mid-latitude total column ozone. These findings increase concern2,12,13 that more frequent and intense wildfires could delay ozone recovery in a warming world.

Global population exposure to landscape fire air pollution from 2000 to 2019.

Wildfires are thought to be increasing in severity and frequency as a result of climate change1-5. Air pollution from landscape fires can negatively affect human health4-6, but human exposure to landscape fire-sourced (LFS) air pollution has not been well characterized at the global scale7-23. Here, we estimate global daily LFS outdoor fine particulate matter (PM2.5) and surface ozone concentrations at 0.25° × 0.25° resolution during the period 2000-2019 with the help of machine learning and chemical transport models. We found that overall population-weighted average LFS PM2.5 and ozone concentrations were 2.5 µg m-3 (6.1% of all-source PM2.5) and 3.2 µg m-3 (3.6% of all-source ozone), respectively, in 2010-2019, with a slight increase for PM2.5, but not for ozone, compared with 2000-2009. Central Africa, Southeast Asia, South America and Siberia experienced the highest LFS PM2.5 and ozone concentrations. The concentrations of LFS PM2.5 and ozone were about four times higher in low-income countries than in high-income countries. During the period 2010-2019, 2.18 billion people were exposed to at least 1 day of substantial LFS air pollution per year, with each person in the world having, on average, 9.9 days of exposure per year. These two metrics increased by 6.8% and 2.1%, respectively, compared with 2000-2009. Overall, we find that the global population is increasingly exposed to LFS air pollution, with socioeconomic disparities.

Global fine-scale changes in ambient NO2 during COVID-19 lockdowns.

Nitrogen dioxide (NO2) is an important contributor to air pollution and can adversely affect human health1-9. A decrease in NO2 concentrations has been reported as a result of lockdown measures to reduce the spread of COVID-1910-20. Questions remain, however, regarding the relationship of satellite-derived atmospheric column NO2 data with health-relevant ambient ground-level concentrations, and the representativeness of limited ground-based monitoring data for global assessment. Here we derive spatially resolved, global ground-level NO2 concentrations from NO2 column densities observed by the TROPOMI satellite instrument at sufficiently fine resolution (approximately one kilometre) to allow assessment of individual cities during COVID-19 lockdowns in 2020 compared to 2019. We apply these estimates to quantify NO2 changes in more than 200 cities, including 65 cities without available ground monitoring, largely in lower-income regions. Mean country-level population-weighted NO2 concentrations are 29% ± 3% lower in countries with strict lockdown conditions than in those without. Relative to long-term trends, NO2 decreases during COVID-19 lockdowns exceed recent Ozone Monitoring Instrument (OMI)-derived year-to-year decreases from emission controls, comparable to 15 ± 4 years of reductions globally. Our case studies indicate that the sensitivity of NO2 to lockdowns varies by country and emissions sector, demonstrating the critical need for spatially resolved observational information provided by these satellite-derived surface concentration estimates.

Wetland emission and atmospheric sink changes explain methane growth in 2020.

Atmospheric methane growth reached an exceptionally high rate of 15.1 ± 0.4 parts per billion per year in 2020 despite a probable decrease in anthropogenic methane emissions during COVID-19 lockdowns1. Here we quantify changes in methane sources and in its atmospheric sink in 2020 compared with 2019. We find that, globally, total anthropogenic emissions decreased by 1.2 ± 0.1 teragrams of methane per year (Tg CH4 yr-1), fire emissions decreased by 6.5 ± 0.1 Tg CH4 yr-1 and wetland emissions increased by 6.0 ± 2.3 Tg CH4 yr-1. Tropospheric OH concentration decreased by 1.6 ± 0.2 per cent relative to 2019, mainly as a result of lower anthropogenic nitrogen oxide (NOx) emissions and associated lower free tropospheric ozone during pandemic lockdowns2. From atmospheric inversions, we also infer that global net emissions increased by 6.9 ± 2.1 Tg CH4 yr-1 in 2020 relative to 2019, and global methane removal from reaction with OH decreased by 7.5 ± 0.8 Tg CH4 yr-1. Therefore, we attribute the methane growth rate anomaly in 2020 relative to 2019 to lower OH sink (53 ± 10 per cent) and higher natural emissions (47 ± 16 per cent), mostly from wetlands. In line with previous findings3,4, our results imply that wetland methane emissions are sensitive to a warmer and wetter climate and could act as a positive feedback mechanism in the future. Our study also suggests that nitrogen oxide emission trends need to be taken into account when implementing the global anthropogenic methane emissions reduction pledge5.

A pause in Southern Hemisphere circulation trends due to the Montreal Protocol.

Observations show robust near-surface trends in Southern Hemisphere tropospheric circulation towards the end of the twentieth century, including a poleward shift in the mid-latitude jet1,2, a positive trend in the Southern Annular Mode1,3-6 and an expansion of the Hadley cell7,8. It has been established that these trends were driven by ozone depletion in the Antarctic stratosphere due to emissions of ozone-depleting substances9-11. Here we show that these widely reported circulation trends paused, or slightly reversed, around the year 2000. Using a pattern-based detection and attribution analysis of atmospheric zonal wind, we show that the pause in circulation trends is forced by human activities, and has not occurred owing only to internal or natural variability of the climate system. Furthermore, we demonstrate that stratospheric ozone recovery, resulting from the Montreal Protocol, is the key driver of the pause. Because pre-2000 circulation trends have affected precipitation12-14, and potentially ocean circulation and salinity15-17, we anticipate that a pause in these trends will have wider impacts on the Earth system. Signatures of the effects of the Montreal Protocol and the associated stratospheric ozone recovery might therefore manifest, or have already manifested, in other aspects of the Earth system.

Assessing the costs of ozone pollution in India for wheat producers, consumers, and government food welfare policies.

We assess wheat yield losses occurring due to ozone pollution in India and its economic burden on producers, consumers, and the government. Applying an ozone flux-based risk assessment, we show that ambient ozone levels caused a mean 14.18% reduction in wheat yields during 2008 to 2012. Furthermore, irrigated wheat was particularly sensitive to ozone-induced yield losses, indicating that ozone pollution could undermine climate-change adaptation efforts through irrigation expansion. Applying an economic model, we examine the effects of a counterfactual, "pollution-free" scenario on yield losses, wheat prices, consumer and producer welfare, and government costs. We explore three policy scenarios in which the government support farmers at observed levels of either procurement prices (fixed-price), procurement quantities (fixed-procurement), or procurement expenditure (fixed-expenditure). In pollution-free conditions, the fixed-price scenario absorbs the fall in prices, thus increasing producer welfare by USD 2.7 billion, but total welfare decreases by USD 0.24 billion as government costs increase (USD 2.9 billion). In the fixed-procurement and fixed-expenditure scenarios, ozone mitigation allows wheat prices to fall by 38.19 to 42.96%. The producers lose by USD 5.10 to 6.01 billion, but the gains to consumers and governments (USD 8.7 to 10.2 billion) outweigh these losses. These findings show that the government and consumers primarily bear the costs of ozone pollution. For pollution mitigation to optimally benefit wheat production and maximize social welfare, new approaches to support producers other than fixed-price grain procurement may be required. We also emphasize the need to consider air pollution in programs to improve agricultural resilience to climate change.

Plant biochemistry influences tropospheric ozone formation, destruction, deposition, and response.

Tropospheric ozone (O3) is among the most damaging air pollutant to plants. Plants alter the atmospheric O3 concentration in two distinct ways: (i) by the emission of volatile organic compounds (VOCs) that are precursors of O3; and (ii) by dry deposition, which includes diffusion of O3 into vegetation through stomata and destruction by nonstomatal pathways. Isoprene, monoterpenes, and higher terpenoids are emitted by plants in quantities that alter tropospheric O3. Deposition of O3 into vegetation is related to stomatal conductance, leaf structural traits, and the detoxification capacity of the apoplast. The biochemical fate of O3 once it enters leaves and reacts with aqueous surfaces is largely unknown, but new techniques for the tracking and identification of initial products have the potential to open the black box.

Isotopic constraint on the twentieth-century increase in tropospheric ozone.

Tropospheric ozone (O3) is a key component of air pollution and an important anthropogenic greenhouse gas1. During the twentieth century, the proliferation of the internal combustion engine, rapid industrialization and land-use change led to a global-scale increase in O3 concentrations2,3; however, the magnitude of this increase is uncertain. Atmospheric chemistry models typically predict4-7 an increase in the tropospheric O3 burden of between 25 and 50 per cent since 1900, whereas direct measurements made in the late nineteenth century indicate that surface O3 mixing ratios increased by up to 300 per cent8-10 over that time period. However, the accuracy and diagnostic power of these measurements remains controversial2. Here we use a record of the clumped-isotope composition of molecular oxygen (18O18O in O2) trapped in polar firn and ice from 1590 to 2016 AD, as well as atmospheric chemistry model simulations, to constrain changes in tropospheric O3 concentrations. We find that during the second half of the twentieth century, the proportion of 18O18O in O2 decreased by 0.03 ± 0.02 parts per thousand (95 per cent confidence interval) below its 1590-1958 AD mean, which implies that tropospheric O3 increased by less than 40 per cent during that time. These results corroborate model predictions of global-scale increases in surface pollution and vegetative stress caused by increasing anthropogenic emissions of O3 precursors4,5,11. We also estimate that the radiative forcing of tropospheric O3 since 1850 AD is probably less than +0.4 watts per square metre, consistent with results from recent climate modelling studies12.

The Antarctic ozone hole and the pattern effect on climate sensitivity.

Since about 1980, the tropical Pacific has been anomalously cold, while the broader tropics have warmed. This has caused anomalous weather in midlatitudes as well as a reduction in the apparent sensitivity of the climate associated with enhanced low-cloud abundance over the cooler waters of the eastern tropical Pacific. Recent modeling work has shown that cooler temperatures over the Southern Ocean around Antarctica can lead to cooler temperatures over the eastern tropical Pacific. Here we suggest that surface wind anomalies associated with the Antarctic ozone hole can cause cooler temperatures over the Southern Ocean that extend into the tropics. We use the short-term variability of the Southern Annular Mode of zonal wind variability to show an association between surface zonal wind variations over the Southern Ocean, cooling over the Southern Ocean, and cooling in the eastern tropical Pacific. This suggests that the cooling of the eastern tropical Pacific may be associated with the onset of the Antarctic ozone hole.

Estimate of OH trends over one decade in North American cities.

The hydroxyl radical (OH) is the most important oxidant on global and local scales in the troposphere. Urban OH controls the removal rate of primary pollutants and triggers the production of ozone. Interannual trends of OH in urban areas are not well documented or understood due to the short lifetime and high spatial heterogeneity of OH. We utilize machine learning with observational inputs emphasizing satellite remote sensing observations to predict surface OH in 49 North American cities from 2005 to 2014. We observe changes in the summertime OH over one decade, with wide variation among different cities. In 2014, compared to the summertime OH in 2005, 3 cities show a significant increase of OH, whereas, in 27 cities, OH decreases in 2014. The year-to-year variation of OH is mapped to the decline of the NO2 column. We conclude that these cities in this analysis are either in the NOx-limited regime or at the transition from a NOx suppressed regime to a NOx-limited regime. The result emphasizes that, in the future, controlling NOx emissions will be most effective in regulating the ozone pollution in these cities.

Long-term tropospheric ozone pollution disrupts plant-microbe-soil interactions in the agroecosystem.

Tropospheric ozone (O3 ) threatens agroecosystems, yet its long-term effects on intricate plant-microbe-soil interactions remain overlooked. This study employed two soybean genotypes of contrasting O3 -sensitivity grown in field plots exposed elevated O3 (eO3 ) and evaluated cause-effect relationships with their associated soil microbiomes and soil quality. Results revealed long-term eO3 effects on belowground soil microbiomes and soil health surpass damage visible on plants. Elevated O3 significantly disrupted belowground bacteria-fungi interactions, reduced fungal diversity, and altered fungal community assembly by impacting soybean physiological properties. Particularly, eO3 impacts on plant performance were significantly associated with arbuscular mycorrhizal fungi, undermining their contribution to plants, whereas eO3 increased fungal saprotroph proliferation, accelerating soil organic matter decomposition and soil carbon pool depletion. Free-living diazotrophs exhibited remarkable acclimation under eO3 , improving plant performance by enhancing nitrogen fixation. However, overarching detrimental consequences of eO3 negated this benefit. Overall, this study demonstrated long-term eO3 profoundly governed negative impacts on plant-soil-microbiota interactions, pointing to a potential crisis for agroecosystems. These findings highlight urgent needs to develop adaptive strategies to navigate future eO3 scenarios.

Multiphase Ozonolysis of Bisphenol A: Chemical Transformations on Surfaces in the Environment.

High global plastic production volumes have led to the widespread presence of bisphenol compounds in human living and working environments. The most common bisphenol, bisphenol A (BPA), despite being endocrine disruptive and estrogenic, is still not fully banned worldwide, leading to continued human exposure via particles in air, dust, and surfaces in both outdoor and indoor environments. While its abundance is well documented, few studies have addressed the chemical transformations of BPA, the properties of its reactive products, and their toxicity. Here, the first gas-surface multiphase ozonolysis experiment of BPA thin films, at a constant ozone mixing ratio of 100 ppb, was performed in a flow tube for periods up to 24 h. Three transformation products involving the addition of 1, 2, and 3 oxygen atoms to the molecule were identified by LC-ESI-HRMS analyses. Exposure of indoor air to thin BPA surface films and BPA-containing thermal paper over periods of days validated the flow tube experiments, demonstrating the rapid nature of this multiphase ozonolysis reaction at atmospherically relevant ozone levels. Multiple transformation pathways are proposed that are likely applicable to not only BPA but also emerging commercial bisphenol products.

How Does Personal Hygiene Influence Indoor Air Quality?

Humans are known to be a continuous and potent indoor source of volatile organic compounds (VOCs). However, little is known about how personal hygiene, in terms of showering frequency, can influence these emissions and their impact on indoor air chemistry involving ozone. In this study, we characterized the VOC composition of the air in a controlled climate chamber (22.5 m3 with an air change rate at 3.2 h-1) occupied by four male volunteers on successive days under ozone-free (∼0 ppb) and ozone-present (37-40 ppb) conditions. The volunteers either showered the evening prior to the experiments or skipped showering for 24 and 48 h. Reduced shower frequency increased human emissions of gas-phase carboxylic acids, possibly originating from skin bacteria. With ozone present, increasing the number of no-shower days enhanced ozone-skin surface reactions, yielding higher levels of oxidation products. Wearing the same clothing over several days reduced the level of compounds generated from clothing-ozone reactions. When skin lotion was applied, the yield of the skin ozonolysis products decreased, while other compounds increased due to ozone reactions with lotion ingredients. These findings help determine the degree to which personal hygiene choices affect the indoor air composition and indoor air exposures.

Squalene Depletion in Skin Following Human Exposure to Ozone under Controlled Chamber Conditions.

A major component of human skin oil is squalene, a highly unsaturated hydrocarbon that protects the skin from atmospheric oxidants. Skin oil, and thus squalene, is continuously replenished on the skin surface. Squalene is also quickly consumed through reactions with ozone and other oxidants. This study examined the extent of squalene depletion in the skin oils of the forearm of human volunteers after exposure to ozone in a climate chamber. Temperature, relative humidity (RH), skin coverage by clothing, and participants' age were varied in a controlled manner. Concentrations of squalene were determined in skin wipe samples collected before and after ozone exposure. Exposures to ozone resulted in statistically significant decreases in post-exposure squalene concentrations compared to pre-exposure squalene concentrations in the skin wipes when squalene concentrations were normalized by concentrations of co-occurring cholesterol but not by co-occurring pyroglutamic acid (PGA). The rate of squalene loss due to ozonolysis was lower than its replenishment on the skin surface. Within the ranges examined, temperature and RH did not significantly affect the difference between normalized squalene levels in post-samples versus pre-samples. Although not statistically significant, skin coverage and age of the volunteers (three young adults, three seniors, and three teenagers) did appear to impact squalene depletion on the skin surfaces.

Changing Responses of PM2.5 and Ozone to Source Emissions in the Yangtze River Delta Using the Adjoint Model.

China's industrial restructuring and pollution controls have altered the contributions of individual sources to varying air quality over the past decade. We used the GEOS-Chem adjoint model and investigated the changing sensitivities of PM2.5 and ozone (O3) to multiple species and sources from 2010 to 2020 in the central Yangtze River Delta (YRDC), the largest economic region in China. Controlling primary particles and SO2 from industrial and residential sectors dominated PM2.5 decline, and reducing CO from multiple sources and ≥C3 alkenes from vehicles restrained O3. The chemical regime of O3 formation became less VOC-limited, attributable to continuous NOX abatement for specific sources, including power plants, industrial combustion, cement production, and off-road traffic. Regional transport was found to be increasingly influential on PM2.5. To further improve air quality, management of agricultural activities to reduce NH3 is essential for alleviating PM2.5 pollution, while controlling aromatics, alkenes, and alkanes from industry and gasoline vehicles is effective for O3. Reducing the level of NOX from nearby industrial combustion and transportation is helpful for both species. Our findings reveal the complexity of coordinating control of PM2.5 and O3 pollution in a fast-developing region and support science-based policymaking for other regions with similar air pollution problems.

Evolution of Ozone Formation Sensitivity during a Persistent Regional Ozone Episode in Northeastern China and Its Implication for a Control Strategy.

In recent years, the magnitude and frequency of regional ozone (O3) episodes have increased in China. We combined ground-based measurements, observation-based model (OBM), and the Weather Research and Forecasting and Community Multiscale Air Quality (WRF-CMAQ) model to analyze a typical persistent O3 episode that occurred across 88 cities in northeastern China during June 19-30, 2021. The meteorological conditions, particularly the wind convergence centers, played crucial roles in the evolution of O3 pollution. Daily analysis of the O3 formation sensitivity showed that O3 formation was in the volatile organic compound (VOC)-limited or transitional regime at the onset of the pollution episode in 92% of the cities. Conversely, it tended to be or eventually became a NOx-limited regime as the episode progressed in the most polluted cities. Based on the emission-reduction scenario simulations, mitigation of the regional O3 pollution was found to be most effective through a phased control strategy, namely, reduction of a high ratio of VOCs to NOx at the onset of the pollution and lower ratio during evolution of the O3 episode. This study presents a new possibility for regional O3 pollution abatement in China based on a reasonable combination of OBM and the WRF-CMAQ model.

Satellite NO2 Retrieval Complicated by Aerosol Composition over Global Urban Agglomerations: Seasonal Variations and Long-Term Trends (2001-2018).

Tropospheric nitrogen dioxide (NO2) poses a serious threat to the environmental quality and public health. Satellite NO2 observations have been continuously used to monitor NO2 variations and improve model performances. However, the accuracy of satellite NO2 retrieval depends on the knowledge of aerosol optical properties, in particular for urban agglomerations accompanied by significant changes in aerosol characteristics. In this study, we investigate the impacts of aerosol composition on tropospheric NO2 retrieval for an 18 year global data set from Global Ozone Monitoring Experiment (GOME)-series satellite sensors. With a focus on cloud-free scenes dominated by the presence of aerosols, individual aerosol composition affects the uncertainties of tropospheric NO2 columns through impacts on the aerosol loading amount, relative vertical distribution of aerosol and NO2, aerosol absorption properties, and surface albedo determination. Among aerosol compositions, secondary inorganic aerosol mostly dominates the NO2 uncertainty by up to 43.5% in urban agglomerations, while organic aerosols contribute significantly to the NO2 uncertainty by -8.9 to 37.3% during biomass burning seasons. The possible contrary influences from different aerosol species highlight the importance and complexity of aerosol correction on tropospheric NO2 retrieval and indicate the need for a full picture of aerosol properties. This is of particular importance for interpreting seasonal variations or long-term trends of tropospheric NO2 columns as well as for mitigating ozone and fine particulate matter pollution.

Vertical Evolution of Ozone Formation Sensitivity Based on Synchronous Vertical Observations of Ozone and Proxies for Its Precursors: Implications for Ozone Pollution Prevention Strategies.

Photochemical ozone (O3) formation in the atmospheric boundary layer occurs at both the surface and elevated altitudes. Therefore, the O3 formation sensitivity is needed to be evaluated at different altitudes before formulating an effective O3 pollution prevention and control strategy. Herein, we explore the vertical evolution of O3 formation sensitivity via synchronous observations of the vertical profiles of O3 and proxies for its precursors, formaldehyde (HCHO) and nitrogen dioxide (NO2), using multi-axis differential optical absorption spectroscopy (MAX-DOAS) in urban areas of the Beijing-Tianjin-Hebei (BTH), Yangtze River Delta (YRD), and Pearl River Delta (PRD) regions in China. The sensitivity thresholds indicated by the HCHO/NO2 ratio (FNR) varied with altitude. The VOC-limited regime dominated at the ground level, whereas the contribution of the NOx-limited regime increased with altitude, particularly on heavily polluted days. The NOx-limited and transition regimes played more important roles throughout the entire boundary layer than at the surface. The feasibility of extreme NOx reduction to mitigate the extent of the O3 pollution was evaluated using the FNR-O3 curve. Based on the surface sensitivity, the critical NOx reduction percentage for the transition from a VOC-limited to a NOx-limited regime is 45-72%, which will decrease to 27-61% when vertical evolution is considered. With the combined effects of clean air action and carbon neutrality, O3 pollution in the YRD and PRD regions will transition to the NOx-limited regime before 2030 and be mitigated with further NOx reduction.

Temperature-Dependent Evaporative Anthropogenic VOC Emissions Significantly Exacerbate Regional Ozone Pollution.

The evaporative emissions of anthropogenic volatile organic compounds (AVOCs) are sensitive to ambient temperature. This sensitivity forms an air pollution-meteorology connection that has not been assessed on a regional scale. We parametrized the temperature dependence of evaporative AVOC fluxes in a regional air quality model and evaluated the impacts on surface ozone in the Beijing-Tianjin-Hebei (BTH) area of China during the summer of 2017. The temperature dependency of AVOC emissions drove an enhanced simulated ozone-temperature sensitivity of 1.0 to 1.8 µg m-3 K-1, comparable to the simulated ozone-temperature sensitivity driven by the temperature dependency of biogenic VOC emissions (1.7 to 2.4 µg m-3 K-1). Ozone enhancements driven by temperature-induced AVOC increases were localized to their point of emission and were relatively more important in urban areas than in rural regions. The inclusion of the temperature-dependent AVOC emissions in our model improved the simulated ozone-temperature sensitivities on days of ozone exceedance. Our results demonstrated the importance of temperature-dependent AVOC emissions on surface ozone pollution and its heretofore unrepresented role in air pollution-meteorology interactions.

Short-Term Increases in NO2 and O3 Concentrations during Pregnancy and Stillbirth Risk in the U.S.: A Time-Stratified Case-Crossover Study.

Associations between gaseous pollutant exposure and stillbirth have focused on exposures averaged over trimesters or gestation. We investigated the association between short-term increases in nitrogen dioxide (NO2) and ozone (O3) concentrations and stillbirth risk among a national sample of 116 788 Medicaid enrollees from 2000 to 2014. A time-stratified case-crossover design was used to estimate distributed (lag 0-lag 6) and cumulative lag effects, which were adjusted for PM2.5 concentration and temperature. Effect modification by race/ethnicity and proximity to hydraulic fracturing (fracking) wells was assessed. Short-term increases in the NO2 and O3 concentrations were not associated with stillbirth in the overall sample. Among American Indian individuals (n = 1694), a 10 ppb increase in NO2 concentrations was associated with increased stillbirth odds at lag 0 (5.66%, 95%CI: [0.57%, 11.01%], p = 0.03) and lag 1 (4.08%, 95%CI: [0.22%, 8.09%], p = 0.04) but not lag 0-6 (7.12%, 95%CI: [-9.83%, 27.27%], p = 0.43). Among participants living in zip codes within 15 km of active fracking wells (n = 9486), a 10 ppb increase in NO2 concentration was associated with increased stillbirth odds in single-day lags (2.42%, 95%CI: [0.37%, 4.52%], p = 0.02 for lag 0 and 1.83%, 95%CI: [0.25%, 3.43%], p = 0.03 for lag 1) but not the cumulative lag (lag 0-6) (4.62%, 95%CI: [-2.75%, 12.55%], p = 0.22). Odds ratios were close to the null in zip codes distant from fracking wells. Future studies should investigate the role of air pollutants emitted from fracking and potential racial disparities in the relationship between short-term increases in NO2 concentrations and stillbirth.