Lifetime and distribution of ozone and related pollutants in the eastern United States

Most major cities in the eastern United States have air quality deemed unhealthy by the EPA under a set of regulations known as the National Ambient Air Quality Standards (NAAQS). The worst air quality in Maryland is measured in Edgewood, MD, a small community located along the Chesapeake Bay and generally downwind of Baltimore during hot, summertime days. Direct measurements and numerical simulations were used to investigate how meteorology and chemistry conspire to create adverse levels of photochemical smog especially at this coastal location. Ozone (O3) and oxidized reactive nitrogen (NOy), a family of ozone precursors, were measured over the Chesapeake Bay during a ten day experiment in July 2011 to better understand the formation of ozone over the Bay and its impact on coastal communities such as Edgewood. Ozone over the Bay during the afternoon was 10% to 20% higher than the closest upwind ground sites. A combination of complex boundary layer dynamics, deposition rates, and unaccounted marine emissions play an integral role in the regional maximum of ozone over the Bay. The CAMx regional air quality model was assessed and enhanced through comparison with data from NASA’s 2011 DISCOVER-AQ field campaign. Comparisons show a model overestimate of NOy by +86.2% and a model underestimate of formaldehyde (HCHO) by –28.3%. I present a revised model framework that better captures these observations and the response of ozone to reductions of precursor emissions. Incremental controls on electricity generating stations will produce greater benefits for surface ozone while additional controls on mobile sources may yield less benefit because cars emit less pollution than expected. Model results also indicate that as ozone concentrations improve with decreasing anthropogenic emissions, the photochemical lifetime of tropospheric ozone increases. The lifetime of ozone lengthens because the two primary gas-phase sinks for odd oxygen (Ox ≈ NO2 + O3) – attack by hydroperoxyl radicals (HO2) on ozone and formation of nitrate – weaken with decreasing pollutant emissions. This unintended consequence of air quality regulation causes pollutants to persist longer in the atmosphere, and indicates that pollutant transport between states and countries will likely play a greater role in the future.

AN EXAMINATION OF HYDROXYL RADICAL: OUR CURRENT UNDERSTANDING OF THE OXIDATIVE CAPACITY OF THE TROPOSPHERE THROUGH EMPIRICAL, BOX, AND GLOBAL MODELING APPROACHES

Hydroxyl radical (OH) is the primary oxidant in the troposphere, initiating the removal of numerous atmospheric species including greenhouse gases, pollutants that are detrimental to human health, and ozone-depleting substances. Because of the complexity of OH chemistry, models vary widely in their OH chemistry schemes and resulting methane (CH4) lifetimes. The current state of knowledge concerning global OH abundances is often contradictory. This body of work encompasses three projects that investigate tropospheric OH from a modeling perspective, with the goal of improving the tropospheric community’s knowledge of the atmospheric lifetime of CH4. First, measurements taken during the airborne CONvective TRansport of Active Species in the Tropics (CONTRAST) field campaign are used to evaluate OH in global models. A box model constrained to measured variables is utilized to infer concentrations of OH along the flight track. Results are used to evaluate global model performance, suggest against the existence of a proposed “OH Hole” in the tropical Western Pacific, and investigate implications of high O3/low H2O filaments on chemical transport to the stratosphere. While methyl chloroform-based estimates of global mean OH suggest that models are overestimating OH, we report evidence that these models are actually underestimating OH in the tropical Western Pacific. The second project examines OH within global models to diagnose differences in CH4 lifetime. I developed an approach to quantify the roles of OH precursor field differences (O3, H2O, CO, NOx, etc.) using a neural network method. This technique enables us to approximate the change in CH4 lifetime resulting from variations in individual precursor fields. The dominant factors driving CH4 lifetime differences between models are O3, CO, and J(O3-O1D). My third project evaluates the effect of climate change on global fields of OH using an empirical model. Observations of H2O and O3 from satellite instruments are combined with a simulation of tropical expansion to derive changes in global mean OH over the past 25 years. We find that increasing H2O and increasing width of the tropics tend to increase global mean OH, countering the increasing CH4 sink and resulting in well-buffered global tropospheric OH concentrations.

Observed Teleconnections in Northern Winter: Subseasonal Evolution and Tropical Linkages

Teleconnections refer to the climate variability links between non-contiguous geographic regions, and tend to be associated with variability in both space and time of the climate’s semi-permanent circulation features. Teleconnections are well-developed in Northern winter, when they influence subseasonal-to-seasonal climate variability, notably, in surface temperature and precipitation. This work is comprised of four independent studies that improve understanding of tropical-extratropical teleconnections and their surface climate responses, subseasonal teleconnection evolution, and the utility of teleconnections in attribution of extreme climate events. After an introduction to teleconnection analysis as well as the major teleconnection patterns and associated climatic footprints manifest during Northern winter, the lagged impact of the Madden-Julian Oscillation (MJO) on subseasonal climate variability is presented. It is found that monitoring of MJO-related velocity potential anomalies is sufficient to predict MJO impacts. These impacts include, for example, the development of significant positive temperature anomalies over the eastern United States one to three weeks following an anomalous convective dipole with enhanced (suppressed) convection centered over the Indian Ocean (western Pacific). Subseasonal teleconnection evolution is assessed with respect to the Pacific-North America (PNA) pattern and the North Atlantic Oscillation (NAO). This evolution is analyzed both in the presence and absence of MJO-related circulation anomalies. It is found that removal of the MJO results only in small shifts in the centers of action of the NAO and PNA, and that in either case there is a small but significant lag in which the NAO leads a PNA pattern of opposite phase. Barotropic vorticity analysis suggests that this relationship may result in part from excitation of Rossby waves by the NAO in the Asian waveguide. An attempt is made to elegantly differentiate between the MJO extratropical response and patterns of variability more internal to the extratropics. Analysis of upper-level streamfunction anomalies is successful in this regard, and it is suggested that this is the preferred method for the real time monitoring of tropical-extratropical teleconnections. The extreme 2013-2014 North American winter is reconstructed using teleconnection analysis, and it is found that the North Pacific Oscillation-West Pacific (NPO/WP) pattern was the leading contributor to climate anomalies over much of North America. Such attribution is cautionary given the propensity to implicate the tropics for all midlatitude climate anomalies based on the El Niño-Southern Oscillation (ENSO) paradigm. A recent hypothesis of such tropical influence is presented and challenged.

Photochemistry and transport of tropospheric ozone and its precursors in urban and remote environments

Tropospheric ozone (O3) adversely affects human health, reduces crop yields, and contributes to climate forcing. To limit these effects, the processes controlling O3 abundance as well as that of its precursor molecules must be fully characterized. Here, I examine three facets of O3 production, both in heavily polluted and remote environments. First, using in situ observations from the DISCOVER-AQ field campaign in the Baltimore/Washington region, I evaluate the emissions of the O3 precursors CO and NOx (NOx = NO + NO2) in the National Emissions Inventory (NEI). I find that CO/NOx emissions ratios derived from observations are 21% higher than those predicted by the NEI. Comparisons to output from the CMAQ model suggest that CO in the NEI is accurate within 15 ± 11%, while NOx emissions are overestimated by 51-70%, likely due to errors in mobile sources. These results imply that ambient ozone concentrations will respond more efficiently to NOx controls than current models suggest. I then investigate the source of high O3 and low H2O structures in the Tropical Western Pacific (TWP). A combination of in situ observations, satellite data, and models show that the high O3 results from photochemical production in biomass burning plumes from fires in tropical Southeast Asia and Central Africa; the low relative humidity results from large-scale descent in the tropics. Because these structures have frequently been attributed to mid-latitude pollution, biomass burning in the tropics likely contributes more to the radiative forcing of climate than previously believed. Finally, I evaluate the processes controlling formaldehyde (HCHO) in the TWP. Convective transport of near surface HCHO leads to a 33% increase in upper tropospheric HCHO mixing ratios; convection also likely increases upper tropospheric CH3OOH to ~230 pptv, enough to maintain background HCHO at ~75 pptv. The long-range transport of polluted air, with NO four times the convectively controlled background, intensifies the conversion of HO2 to OH, increasing OH by a factor of 1.4. Comparisons between the global chemistry model CAM-Chem and observations show that consistent underestimates of HCHO by CAM-Chem throughout the troposphere result from underestimates in both NO and acetaldehyde.

NITROUS OXIDE EMISSIONS IN COVER CROP-BASED CORN PRODUCTION SYSTEMS

Nitrous oxide (N2O) is a potent greenhouse gas; the majority of N2O emissions are the result of agricultural management, particularly the application of N fertilizers to soils. The relationship of N2O emissions to varying sources of N (manures, mineral fertilizers, and cover crops) has not been well-evaluated. Here we discussed a novel methodology for estimating precipitation-induced pulses of N2O using flux measurements; results indicated that short-term intensive time-series sampling methods can adequately describe the magnitude of these pulses. We also evaluated the annual N2O emissions from corn-cover crop (Zea mays; cereal rye [Secale cereale], hairy vetch [Vicia villosa], or biculture) production systems when fertilized with multiple rates of subsurface banded poultry litter, as compared with tillage incorporation or mineral fertilizer. N2O emissions increased exponentially with total N rate; tillage decreased emissions following cover crops with legume components, while the effect of mineral fertilizer was mixed across cover crops.

THE IMPACT OF UPPER-LEVEL PROCESSES ON THE INTENSITY AND STRUCTURAL CHANGES OF HURRICANE SANDY (2012)

The first part of this study examines the relative roles of frontogenesis and tropopause undulation in determining the intensity and structural changes of Hurricane Sandy (2012) using a high-resolution cloud-resolving model. A 138-h simulation reproduces Sandy’s four distinct development stages: (i) rapid intensification, (ii) weakening, (iii) steady maximum surface wind but with large continued sea-level pressure (SLP) falls, and (iv) re-intensification. Results show typical correlations between intensity changes, sea-surface temperature and vertical wind shear during the first two stages. The large SLP falls during the last two stages are mostly caused by Sandy’s moving northward into lower-tropopause regions associated with an eastward-propagating midlatitude trough, where the associated lower-stratospheric warm air wraps into the storm and its surrounding areas. The steady maximum surface wind occurs because of the widespread SLP falls with weak pressure gradients lacking significant inward advection of absolute angular momentum (AAM). Meanwhile, there is a continuous frontogenesis in the outer region during the last three stages. Cyclonic inward advection of AAM along each frontal rainband accounts for the continued expansion of the tropical-storm-force wind and structural changes, while deep convection in the eyewall and merging of the final two survived frontal rainbands generate a spiraling jet in Sandy’s northwestern quadrant, leading to its re-intensification prior to landfall. The physical, kinematic and dynamic aspects of an upper-level outflow layer and its possible impact on the re-intensification of Sandy are examined in the second part of this study. Above the outflow layer isentropes are tilted downward with radius as a result of the development of deep convection and an approaching upper-level trough, causing weak subsidence. Its maximum outward radial velocity is located above the cloud top, so the outflow channel experiences cloud-induced long-wave cooling. Because Sandy has two distinct convective regions (an eyewall and a frontal rainband), it has multiple outflow layers, with the eyewall’s outflow layer located above that of the frontal rainband. During the re-intensification stage, the eyewall’s outflow layer interacts with a jet stream ahead of the upper-level trough axis. Because of the presence of inertial instability on the anticyclonic side of the jet stream and symmetric instability in the inner region of the outflow layer, Sandy’s secondary circulation intensifies. Its re-intensification ceases when these instabilities disappear. The relationship between the intensity of the secondary circulation and dynamic instabilities of the outflow layer suggests that the re-intensification occurs in response to these instabilities. Additionally, it is verified that the long-wave cooling in the outflow layer helps induce symmetric instability by reducing static stability.

OBSERVATIONS AND EMISSIONS OF ENERGY-ASSOCIATED OZONE PRECURSORS IN THE MID-ATLANTIC UNITED STATES

Surface ozone is formed in the presence of NOx (NO + NO2) and volatile organic compounds (VOCs) and is hazardous to human health. A better understanding of these precursors is needed for developing effective policies to improve air quality. To evaluate the year-to-year changes in source contributions to total VOCs, Positive Matrix Factorization (PMF) was used to perform source apportionment using available hourly observations from June through August at a Photochemical Assessment Monitoring Station (PAMS) in Essex, MD for each year from 2007-2015. Results suggest that while gasoline and vehicle exhaust emissions have fallen, the contribution of natural gas sources to total VOCs has risen. To investigate this increasing natural gas influence, ethane measurements from PAMS sites in Essex, MD and Washington, D.C. were examined. Following a period of decline, daytime ethane concentrations have increased significantly after 2009. This trend appears to be linked with the rapid shale gas production in upwind, neighboring states, especially Pennsylvania and West Virginia. Back-trajectory analyses similarly show that ethane concentrations at these monitors were significantly greater if air parcels had passed through counties containing a high density of unconventional natural gas wells. In addition to VOC emissions, the compressors and engines involved with hydraulic fracturing operations also emit NOx and particulate matter (PM). The Community Multi-scale Air Quality (CMAQ) Model was used to simulate air quality for the Eastern U.S. in 2020, including emissions from shale gas operations in the Appalachian Basin. Predicted concentrations of ozone and PM show the largest decreases when these natural gas resources are hypothetically used to convert coal-fired power plants, despite the increased emissions from hydraulic fracturing operations expanded into all possible shale regions in the Appalachian Basin. While not as clean as burning natural gas, emissions of NOx from coal-fired power plants can be reduced by utilizing post-combustion controls. However, even though capital investment has already been made, these controls are not always operated at optimal rates. CMAQ simulations for the Eastern U.S. in 2018 show ozone concentrations decrease by ~5 ppb when controls on coal-fired power plants limit NOx emissions to historically best rates.

The role of the Indian Ocean sector and sea surface salinity for prediction of the coupled Indo-Pacific system

The purpose of this dissertation is to evaluate the potential downstream influence of the Indian Ocean (IO) on El Niño/Southern Oscillation (ENSO) forecasts through the oceanic pathway of the Indonesian Throughflow (ITF), atmospheric teleconnections between the IO and Pacific, and assimilation of IO observations. Also the impact of sea surface salinity (SSS) in the Indo-Pacific region is assessed to try to address known problems with operational coupled model precipitation forecasts. The ITF normally drains warm fresh water from the Pacific reducing the mixed layer depths (MLD). A shallower MLD amplifies large-scale oceanic Kelvin/Rossby waves thus giving ~10% larger response and more realistic ENSO sea surface temperature (SST) variability compared to observed when the ITF is open. In order to isolate the impact of the IO sector atmospheric teleconnections to ENSO, experiments are contrasted that selectively couple/decouple the interannual forcing in the IO. The interannual variability of IO SST forcing is responsible for 3 month lagged widespread downwelling in the Pacific, assisted by off-equatorial curl, leading to warmer NINO3 SST anomaly and improved ENSO validation (significant from 3-9 months). Isolating the impact of observations in the IO sector using regional assimilation identifies large-scale warming in the IO that acts to intensify the easterlies of the Walker circulation and increases pervasive upwelling across the Pacific, cooling the eastern Pacific, and improving ENSO validation (r ~ 0.05, RMS~0.08C). Lastly, the positive impact of more accurate fresh water forcing is demonstrated to address inadequate precipitation forecasts in operational coupled models. Aquarius SSS assimilation improves the mixed layer density and enhances mixing, setting off upwelling that eventually cools the eastern Pacific after 6 months, counteracting the pervasive warming of most coupled models and significantly improving ENSO validation from 5-11 months. In summary, the ITF oceanic pathway, the atmospheric teleconnection, the impact of observations in the IO, and improved Indo-Pacific SSS are all responsible for ENSO forecast improvements, and so each aspect of this study contributes to a better overall understanding of ENSO. Therefore, the upstream influence of the IO should be thought of as integral to the functioning of ENSO phenomenon.

Multi-sensor Cloud and Aerosol Retrieval Simulator and Its Applications

Executing a cloud or aerosol physical properties retrieval algorithm from controlled synthetic data is an important step in retrieval algorithm development. Synthetic data can help answer questions about the sensitivity and performance of the algorithm or aid in determining how an existing retrieval algorithm may perform with a planned sensor. Synthetic data can also help in solving issues that may have surfaced in the retrieval results. Synthetic data become very important when other validation methods, such as field campaigns,are of limited scope. These tend to be of relatively short duration and often are costly. Ground stations have limited spatial coverage whilesynthetic data can cover large spatial and temporal scales and a wide variety of conditions at a low cost. In this work I develop an advanced cloud and aerosol retrieval simulator for the MODIS instrument, also known as Multi-sensor Cloud and Aerosol Retrieval Simulator (MCARS). In a close collaboration with the modeling community I have seamlessly combined the GEOS-5 global climate model with the DISORT radiative transfer code, widely used by the remote sensing community, with the observations from the MODIS instrument to create the simulator. With the MCARS simulator it was then possible to solve the long standing issue with the MODIS aerosol optical depth retrievals that had a low bias for smoke aerosols. MODIS aerosol retrieval did not account for effects of humidity on smoke aerosols. The MCARS simulator also revealed an issue that has not been recognized previously, namely,the value of fine mode fraction could create a linear dependence between retrieved aerosol optical depth and land surface reflectance. MCARS provided the ability to examine aerosol retrievals against “ground truth” for hundreds of thousands of simultaneous samples for an area covered by only three AERONET ground stations. Findings from MCARS are already being used to improve the performance of operational MODIS aerosol properties retrieval algorithms. The modeling community will use the MCARS data to create new parameterizations for aerosol properties as a function of properties of the atmospheric column and gain the ability to correct any assimilated retrieval data that may display similar dependencies in comparisons with ground measurements.