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.

Investigation into the meteorological aspects of the variations in the level of Lake Victoria

The level of Lake Victoria has, since 1961, reached a height which caused serious flood damage. Already the financial implications are considerable for Kenya, Tanzania and Uganda. If further rises can be anticipated, expenditure on flood control measures to the tune of several million pounds sterling must be envisaged. If such rises should lead to uncontrolled discharge at the Owen Falls Dam site because of overshooting, downstream districts of Uganda and the Sudan may be seriously flooded. All this merits a thorough study, and any indication of the future behaviour of lake levels, even when associated with a low probability, must be taken into account. In these circumstances the Water Development Department of Kenya approached the East African Meteorological Department in November, 1964, on behalf of all parties concerned with the request to study the meteorological background of the Iake level variation, with a view to forecasting future behaviour.

A New Formulation of the Spectral Energy Budget of the Atmosphere, with Application to Two High-Resolution General Circulation Models

International audience

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.

La forme des récifs coralliens et le régime des vents alternants

This paper deals with the development of crescent-shaped reefs in the China Sea under control of northwest and southeast monsoon winds alternating in direction seasonally. The transport of sediment derived from windward side to leeward around reef ends has been studied. A description of the analogy of these large structures with individual coral colonies developing in similar manner under current action as micro atolls has given.

Seasonal changes of heat and moisture fluxes from the Caspian Sea surface and its effect on synoptic systems with heavy rainfalls in southern shores of the sea during 2005-2010

Caspian Sea with its unique characteristics is a significant source to supply required heat and moisture for passing weather systems over the north of Iran. Investigation of heat and moisture fluxes in the region and their effects on these systems that could lead to floods and major financial and human losses is essential in weather forecasting. Nowadays by improvement of numerical weather and climate prediction models and the increasing need to more accurate forecasting of heavy rainfall, the evaluation and verification of these models has been become much more important. In this study we have used the WRF model as a research-practical one with many valuable characteristics and flexibilities. In this research, the effects of heat and moisture fluxes of Caspian Sea on the synoptic and dynamical structure of 20 selective systems associated with heavy rainfall in the southern shores of Caspian Sea are investigated. These systems are selected based on the rainfall data gathered by three local stations named: Rasht, Babolsar and Gorgan in different seasons during a five-year period (2005-2010) with maximum amount of rainfall through the 24 hours of a day. In addition to synoptic analyses of these systems, the WRF model with and without surface flues was run using the two nested grids with the horizontal resolutions of 12 and 36 km. The results show that there are good consistencies between the predicted distribution of rainfall field, time of beginning and end of rainfall by the model and the observations. But the model underestimates the amounts of rainfall and the maximum difference with the observation is about 69%. Also, no significant changes in the results are seen when the domain and the resolution of computations are changed. The other noticeable point is that the systems are severely weakened by removing heat and moisture fluxes and thereby the amounts of large scale rainfall are decreased up to 77% and the convective rainfalls tend to zero.

Effects of an asymmetric friction on the nonlinear equilibration of a baroclinic system

Following some recent linear and nonlinear studies the authors examine, using numerical simulations of a classical two-layer model, the effect of an asymmetric friction on the nonlinear equilibrium of moderately unstable baroclinic systems, The results show that the presence of an asymmetric friction leads to a significant wave scale selection: ''long'' waves (in terms of their zonal wavelengths) emerge with a traditional asymmetric friction (with the upper layer less viscous than the lower layer), while only ''short'' waves dominate with a nontraditional asymmetric friction (with the lower layer less viscous than the upper layer). The role of the nonlinear interactions and. more precisely, the effects of an asymmetric friction on the wave-mean flow and wave-wave interactions; and their consequences on the wave scale selection are examined.

On the effect of Arabian Sea and Indian Ocean sea surface temperature anomaly on the monsoon rainfall of south-east of Iran

Southeast region of the country has hot and dry weather which causes to happen heavy rainfall in short time period of warm seasons and to occur river flooding. These precipitations are influenced by monsoon system of India ocean. In these thesis, It was tried to evaluate the relation between thermal anomaly of sea surface in India ocean and Arab sea which effects on southeast monsoon precipitations of Iran, For evaluation of this happening in southeast, data were collected from 7 synoptic observation stations of Bandar Abbas, Minab, Kerman , Bam, Chabahar, Iranshahr, Zahedan and 17 rain gauge stations during June to September of each year from 1980 to 2010. Rainy days were determine and then some information about synoptic circulation models, maps of average pressure of sea surface, geopotential height of 700hP surface, geopotential height of 500hP surface, temperature of 850 hPa surface, humidity of 700 hPa surface, vertical velocity of 700 hPa surface, vertical velocity of 500 hP and humidity of 2 meters height for 6 systems were extracted from NCEP/NCAR website for evaluation. By evaluation of these systems it was determined that the monsoon low pressure system tab brings needed humidity of these precipitations to this region from India ocean and Arab sea with a vast circulation. It is seen that warm air pool locates on Iran and cold air pool locates on west of India at 800 hPa surface. In a rainy day this warm air transfers to high latitudes and influences the temperature trough of southeast cold air pool of the country. In the middle surfaces of 700 and 500 hPa, the connection between low height system above India and low height system above the higher latitudes causes the low height system above India to be strength and developed. By evaluation of humidity at 2 meters height and 700 hPa surface we observe that humidity Increases in the southeast region. With penetrating of the low height system of India above the 700 and 500 hPa surfaces of southeast of Iran, the value of negative omega (Rising vertical velocity) is increased. In the second pace, it was shown the evaluation of how the correlation between sea surface temperature anomaly in India Ocean and Arab sea influences southeast monsoon precipitation of Iran. For this purpose the data of water surface temperature anomaly of Arab sea and India ocean, the data of precipitation anomaly of 7 synoptic stations , mentioned above, and correlation coefficient among the data of precipitation anomaly and water surface temperature anomaly of Arab Sea, east and west of India ocean were calculated. In conclusion it was shown that the maximum correlation coefficient of precipitation anomaly had belonged to India Ocean in June and no meaningful correlation was resulted in July among precipitation anomaly and sea surface temperature anomaly for three regions, which were evaluated.

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.