Nitrogen Oxides in the North Atlantic Troposphere: Impacts of Boreal Wildfire and Anthropogenic Emissions

Nitrogen oxides play a crucial role in the budget of tropospheric ozone (O sub(3)) and the formation of the hydroxyl radical. Anthropogenic activities and boreal wildfires are large sources of emissions in the atmosphere. However, the influence of the transport of these emissions on nitrogen oxides and O sub(3) levels at hemispheric scales is not well understood, in particular due to a lack of nitrogen oxides measurements in remote regions. In order to address these deficiencies, measurements of NO, NO sub(2) and NO sub(y) (total reactive nitrogen oxides) were made in the lower free troposphere (FT) over the central North Atlantic region (Pico Mountain station, 38 degree N 28 degree W, 2.3 km asl) from July 2002 to August 2005. These measurements reveal a well-defined seasonal cycle of nitrogen oxides (NO sub(x) = NO+NO sub(2) and NO sub(y)) in the background central North Atlantic lower FT, with higher mixing ratios during the summertime. Observed NO sub(x) and NO sub(y) levels are consistent with long-range transport of emissions, but with significant removal en-route to the measurement site. Reactive nitrogen largely exists in the form of PAN and HNO sub(3) ( similar to 80-90% of NO sub(y)) all year round. A shift in the composition of NO sub(y) from dominance of PAN to dominance of HNO sub(3) occurs from winter-spring to summer-fall, as a result of changes in temperature and photochemistry over the region. Analysis of the long-range transport of boreal wildfire emissions on nitrogen oxides provides evidence of the very large-scale impacts of boreal wildfires on the tropospheric NO sub(x) and O sub(3) budgets. Boreal wildfire emissions are responsible for significant shifts in the nitrogen oxides distributions toward higher levels during the summer, with medians of NO sub(y) (117-175 pptv) and NO sub(x) (9-30 pptv) greater in the presence of boreal wildfire emissions. Extreme levels of NO sub(x) (up to 150 pptv) and NO sub(y) (up to 1100 pptv) observed in boreal wildfire plumes suggest that decomposition of PAN to NO sub(x) is a significant source of NO sub(x), and imply that O sub(3) formation occurs during transport. Ozone levels are also significantly enhanced in boreal wildfire plumes. However, a complex behavior of O sub(3) is observed in the plumes, which varies from significant to lower O sub(3) production to O sub(3) destruction. Long-range transport of anthropogenic emissions from North America also has a significant influence on the regional NO sub(x) and O sub(3) budgets. Transport of pollution from North America causes significant enhancements on nitrogen oxides year-round. Enhancements of CO, NO sub(y) and NO sub(x) indicate that, consistent with previous studies, more than 95% of the NO sub(x) emitted over the U.S. is removed before and during export out of the U.S. boundary layer. However, about 30% of the NO sub(x) emissions exported out of the U.S. boundary layer remain in the airmasses. Since the lifetime of NO sub(x) is shorter than the transport timescale, PAN decomposition and potentially photolysis of HNO sub(3) provide a supply of NO sub(x) over the central North Atlantic lower FT. Observed Delta O sub(3)/ Delta NO sub(y) and large NO sub(y) levels remaining in the North American plumes suggest potential O sub(3) formation well downwind from North America. Finally, a comparison of the nitrogen oxides measurements with results from the global chemical transport (GCT) model GEOS-Chem identifies differences between the observations and the model. GEOS-Chem reproduces the seasonal variation of nitrogen oxides over the central North Atlantic lower FT, but does not capture the magnitude of the cycles. Improvements in our understanding of nitrogen oxides chemistry in the remote FT and emission sources are necessary for the current GCT models to adequately estimate the impacts of emissions on tropospheric NO sub(x) and the resulting impacts on the O sub(3) budget.

INTEGRATED ASSESSMENT OF ANTHROPOGENIC, CLIMATE, AND POLICY INDUCED CHANGES OF PHOSPHORUS EXPORT IN THE UNITED STATES LAURENTIAN GREAT LAKES WATERSHEDS

Anthropogenic activities have increased phosphorus (P) loading in tributaries to the Laurentian Great Lakes resulting in eutrophication in small bays to most notably, Lake Erie. Changes to surface water quality from P loading have resulted in billions of dollars in damage and threaten the health of the world’s largest freshwater resource. To understand the factors affecting P delivery with projected increasing urban lands and biofuels expansion, two spatially explicit models were coupled. The coupled models predict that the majority of the basin will experience a significant increase in urban area P sources while the agriculture intensity and forest sources of P will decrease. Changes in P loading across the basin will be highly variable spatially. Additionally, the impacts of climate change on high precipitation events across the Great Lakes were examined. Using historical regression relationships on phosphorus concentrations, key Great Lakes tributaries were found to have future changes including decreasing total loads and increases to high-flow loading events. The urbanized Cuyahoga watersheds exhibits the most vulnerability to these climate-induced changes with increases in total loading and storm loading , while the forested Au Sable watershed exhibits greater resilience. Finally, the monitoring network currently in place for sampling the amount of phosphorus entering the U.S. Great Lakes was examined with a focus on the challenges to monitoring. Based on these interviews, the research identified three issues that policy makers interested in maintaining an effective phosphorus monitoring network in the Great Lakes should consider: first, that the policy objectives driving different monitoring programs vary, which results in different patterns of sampling design and frequency; second, that these differences complicate efforts to encourage collaboration; and third, that methods of funding sampling programs vary from agency to agency, further complicating efforts to generate sufficient long-term data to improve our understanding of phosphorus into the Great Lakes. The dissertation combines these three areas of research to present the potential future impacts of P loading in the Great Lakes as anthropogenic activities, climate and monitoring changes. These manuscripts report new experimental data for future sources, loading and climate impacts on phosphorus.

Invasion patterns of emerald ash borer and European earthworms in forested ecosystems

Invasive and exotic species present a serious threat to the health and sustainability of natural ecosystems. These species often benefit from anthropogenic activities that aid their introduction and dispersal. This dissertation focuses on invasion dynamics of the emerald ash borer, native to Asia, and European earthworms. These species have shown detrimental impacts in invaded forest ecosystems across the Great Lakes region, and continue to spread via human-assisted long distance dispersal and by natural modes of dispersal into interior forests from areas of introduction. Successful forest management requires that the impact and effect of invasive species be considered and incorporated into management plans. Understanding patterns and constraints of introduction, establishment, and spread will aid in this effort. To assist in efforts to locate introduction points of emerald ash borer, a multicriteria risk model was developed to predict the highest risk areas. Important parameters in the model were road proximity, land cover type, and campground proximity. The model correctly predicted 85% of known emerald ash borer invasion sites to be at high risk. The model’s predictions across northern Michigan can be used to focus and guide future monitoring efforts. Similar modeling efforts were applied to the prediction of European earthworm invasion in northern Michigan forests. Field sampling provided a means to improve upon modeling efforts for earthworms to create current and future predictions of earthworm invasion. Those sites with high soil pH and high basal area of earthworm preferred overstory species (such as basswood and maples) had the highest likelihood of European earthworm invasion. Expanding beyond Michigan into the Upper Great Lakes region, earthworm populations were sampled across six National Wildlife Refuges to identify potential correlates and deduce specific drivers and constraints of earthworm invasion. Earthworm communities across all refuges were influenced by patterns of anthropogenic activity both within refuges and in surrounding ecoregions of study. Forest composition, soil pH, soil organic matter, anthropogenic cover, and agriculture proximity also proved to be important drivers of earthworm abundance and community composition. While there are few management options to remove either emerald ash borer or European earthworms from forests after they have become well established, prevention and early detection are important and can be beneficial. An improved understanding the factors controlling the distribution and invasion patterns of exotic species across the landscape will aid efforts to determine their consequences and generate appropriate forest management solutions to sustain ecosystem health in the presence of these invaders.

Effects of elevated atmospheric CO2 and O3 on wood density, antomical proerties and decomposition of Northern Hardwoods

Anthropogenic activities continue to drive atmospheric CO2 and O3 concentrations to levels higher than during the pre-industrial era. Accumulating evidence indicates that both elevated CO2 and elevated O3 could modify the quantity and biochemistry of woody plant biomass. Anatomical properties of woody plants are largely influenced by the activity of the cambium and the growth characteristics of wood cells, which are in turn influenced by a range of environmental factors. Hence, alterations in the concentrations of atmospheric CO2 and / or O3 could also impact wood anatomical properties. Many fungi derive their metabolic resources for growth from plant litter, including woody tissue, and therefore modifications in the quantity, biochemistry and anatomical properties of woody plants in response to elevated CO2 and / or O3 could impact the community of wood-decaying fungi and rates of wood decomposition. Consequently carbon and nutrient cycling and productivity of terrestrial ecosystem could also be impacted. Alterations in wood structure and biochemistry of woody plants could also impact wood density and subsequently impact wood quality. This dissertation examined the long term effects of elevated CO2 and / or O3 on wood anatomical properties, wood density, wood-decaying fungi and wood decomposition of northern hardwood tree species at the Aspen Free-Air CO2 and O3 Enrichment (Aspen FACE) project, near Rhinelander, WI, USA. Anatomical properties of wood varied significantly with species and aspen genotypes and radial position within the stem. Elevated CO2 did not have significant effects on wood anatomical properties in trembling aspen, paper birch or sugar maple, except for marginally increasing (P < 0.1) the number of vessels per square millimeter. Elevated O3 marginally or significantly altered vessel lumen diameter, cell wall area and vessel lumen area proportions depending on species and radial position. In line with the modifications in the anatomical properties, elevated CO2 and O3, alone, significantly modified wood density but effects were species and / or genotype specific. However, the effects of elevated CO2 and O3, alone, on wood anatomical properties and density were ameliorated when in combination. Wood species had a much greater impact on the wood-decaying fungal community and initial wood decomposition rate than did growth or decomposition of wood in elevated CO2 and / or O3. Polyporales, Agaricales, and Russulales were the dominant orders of fungi isolated. Based on the current results, future higher levels of CO2 and O3 may have moderate effects on wood quality of northern hardwoods, but for utilization purposes these may not be considered significant. However, wood-decaying fungal community composition and decomposition of northern hardwoods may be altered via shifts in species and / or genotype composition under future higher levels of CO2 and O3.

Evaluation of flood risk in response to climate variability

Standard procedures for forecasting flood risk (Bulletin 17B) assume annual maximum flood (AMF) series are stationary, meaning the distribution of flood flows is not significantly affected by climatic trends/cycles, or anthropogenic activities within the watershed. Historical flood events are therefore considered representative of future flood occurrences, and the risk associated with a given flood magnitude is modeled as constant over time. However, in light of increasing evidence to the contrary, this assumption should be reconsidered, especially as the existence of nonstationarity in AMF series can have significant impacts on planning and management of water resources and relevant infrastructure. Research presented in this thesis quantifies the degree of nonstationarity evident in AMF series for unimpaired watersheds throughout the contiguous U.S., identifies meteorological, climatic, and anthropogenic causes of this nonstationarity, and proposes an extension of the Bulletin 17B methodology which yields forecasts of flood risk that reflect climatic influences on flood magnitude. To appropriately forecast flood risk, it is necessary to consider the driving causes of nonstationarity in AMF series. Herein, large-scale climate patterns—including El Niño-Southern Oscillation (ENSO), Pacific Decadal Oscillation (PDO), North Atlantic Oscillation (NAO), and Atlantic Multidecadal Oscillation (AMO)—are identified as influencing factors on flood magnitude at numerous stations across the U.S. Strong relationships between flood magnitude and associated precipitation series were also observed for the majority of sites analyzed in the Upper Midwest and Northeastern regions of the U.S. Although relationships between flood magnitude and associated temperature series are not apparent, results do indicate that temperature is highly correlated with the timing of flood peaks. Despite consideration of watersheds classified as unimpaired, analyses also suggest that identified change-points in AMF series are due to dam construction, and other types of regulation and diversion. Although not explored herein, trends in AMF series are also likely to be partially explained by changes in land use and land cover over time. Results obtained herein suggest that improved forecasts of flood risk may be obtained using a simple modification of the Bulletin 17B framework, wherein the mean and standard deviation of the log-transformed flows are modeled as functions of climate indices associated with oceanic-atmospheric patterns (e.g. AMO, ENSO, NAO, and PDO) with lead times between 3 and 9 months. Herein, one-year ahead forecasts of the mean and standard deviation, and subsequently flood risk, are obtained by applying site specific multivariate regression models, which reflect the phase and intensity of a given climate pattern, as well as possible impacts of coupling of the climate cycles. These forecasts of flood risk are compared with forecasts derived using the existing Bulletin 17B model; large differences in the one-year ahead forecasts are observed in some locations. The increased knowledge of the inherent structure of AMF series and an improved understanding of physical and/or climatic causes of nonstationarity gained from this research should serve as insight for the formulation of a physical-casual based statistical model, incorporating both climatic variations and human impacts, for flood risk over longer planning horizons (e.g., 10-, 50, 100-years) necessary for water resources design, planning, and management.

Assessing success of forest restoration efforts in degraded montane cloud forests in southern Mexico

Montane cloud forests are home to great biodiversity. However, non-sustainable anthropogenic activities have led to the loss of forest cover in southern Mexico. Increasing conservation, restoration and sustainable use of forest resources prevents the loss of cloud forests. In this study, success of forest restoration was evaluated in a degraded forest of Highlands Chiapas. The goal of this study was to assess the structure and composition of native tree species. We evaluated vegetation composition at three sites that had undergone enrichment plantings. Floristic composition and structure of the herbaceous, seedling, sapling, and overstory layers were measured. A total of sixty-six native tree species were recorded. Enrichment planting was found to have increased tree diversity. Moreover, 54% of the planted species were found in the understory, indicating that they were successfully recruiting. In conclusion, enrichment planting can aid in the conservation of forest cover in degraded areas.

Long-Range Pollution Transport: Trans-Atlantic Mechanisms and Lagrangian Modeling Methods

Over the past several decades, it has become apparent that anthropogenic activities have resulted in the large-scale enhancement of the levels of many trace gases throughout the troposphere. More recently, attention has been given to the transport pathway taken by these emissions as they are dispersed throughout the atmosphere. The transport pathway determines the physical characteristics of emissions plumes and therefore plays an important role in the chemical transformations that can occur downwind of source regions. For example, the production of ozone (O3) is strongly dependent upon the transport its precursors undergo. O3 can initially be formed within air masses while still over polluted source regions. These polluted air masses can experience continued O3 production or O3 destruction downwind, depending on the air mass's chemical and transport characteristics. At present, however, there are a number of uncertainties in the relationships between transport and O3 production in the North Atlantic lower free troposphere. The first phase of the study presented here used measurements made at the Pico Mountain observatory and model simulations to determine transport pathways for US emissions to the observatory. The Pico Mountain observatory was established in the summer of 2001 in order to address the need to understand the relationships between transport and O3 production. Measurements from the observatory were analyzed in conjunction with model simulations from the Lagrangian particle dispersion model (LPDM), FLEX-PART, in order to determine the transport pathway for events observed at the Pico Mountain observatory during July 2003. A total of 16 events were observed, 4 of which were analyzed in detail. The transport time for these 16 events varied from 4.5 to 7 days, while the transport altitudes over the ocean ranged from 2-8 km, but were typically less than 3 km. In three of the case studies, eastward advection and transport in a weak warm conveyor belt (WCB) airflow was responsible for the export of North American emissions into the FT, while transport in the FT was governed by easterly winds driven by the Azores/Bermuda High (ABH) and transient northerly lows. In the fourth case study, North American emissions were lofted to 6-8 km in a WCB before being entrained in the same cyclone's dry airstream and transported down to the observatory. The results of this study show that the lower marine FT may provide an important transport environment where O3 production may continue, in contrast to transport in the marine boundary layer, where O3 destruction is believed to dominate. The second phase of the study presented here focused on improving the analysis methods that are available with LPDMs. While LPDMs are popular and useful for the analysis of atmospheric trace gas measurements, identifying the transport pathway of emissions from their source to a receptor (the Pico Mountain observatory in our case) using the standard gridded model output, particularly during complex meteorological scenarios can be difficult can be difficult or impossible. The transport study in phase 1 was limited to only 1 month out of more than 3 years of available data and included only 4 case studies out of the 16 events specifically due to this confounding factor. The second phase of this study addressed this difficulty by presenting a method to clearly and easily identify the pathway taken by only those emissions that arrive at a receptor at a particular time, by combining the standard gridded output from forward (i.e., concentrations) and backward (i.e., residence time) LPDM simulations, greatly simplifying similar analyses. The ability of the method to successfully determine the source-to-receptor pathway, restoring this Lagrangian information that is lost when the data are gridded, is proven by comparing the pathway determined from this method with the particle trajectories from both the forward and backward models. A sample analysis is also presented, demonstrating that this method is more accurate and easier to use than existing methods using standard LPDM products. Finally, we discuss potential future work that would be possible by combining the backward LPDM simulation with gridded data from other sources (e.g., chemical transport models) to obtain a Lagrangian sampling of the air that will eventually arrive at a receptor.

FLOOD RISK ASSESSMENT UNDER HISTORICAL AND PREDICTED LAND USE CHANGE USING CONTINUOUS HYDROLOGIC MODELING

Current procedures for flood risk estimation assume flood distributions are stationary over time, meaning annual maximum flood (AMF) series are not affected by climatic variation, land use/land cover (LULC) change, or management practices. Thus, changes in LULC and climate are generally not accounted for in policy and design related to flood risk/control, and historical flood events are deemed representative of future flood risk. These assumptions need to be re-evaluated, however, as climate change and anthropogenic activities have been observed to have large impacts on flood risk in many areas. In particular, understanding the effects of LULC change is essential to the study and understanding of global environmental change and the consequent hydrologic responses. The research presented herein provides possible causation for observed nonstationarity in AMF series with respect to changes in LULC, as well as a means to assess the degree to which future LULC change will impact flood risk. Four watersheds in the Midwest, Northeastern, and Central United States were studied to determine flood risk associated with historical and future projected LULC change. Historical single framed aerial images dating back to the mid-1950s were used along with Geographic Information Systems (GIS) and remote sensing models (SPRING and ERDAS) to create historical land use maps. The Forecasting Scenarios of Future Land Use Change (FORE-SCE) model was applied to generate future LULC maps annually from 2006 to 2100 for the conterminous U.S. based on the four IPCC-SRES future emission scenario conditions. These land use maps were input into previously calibrated Soil and Water Assessment Tool (SWAT) models for two case study watersheds. In order to isolate effects of LULC change, the only variable parameter was the Runoff Curve Number associated with the land use layer. All simulations were run with daily climate data from 1978-1999, consistent with the 'base' model which employed the 1992 NLCD to represent 'current' conditions. Output daily maximum flows were converted to instantaneous AMF series and were subsequently modeled using a Log-Pearson Type 3 (LP3) distribution to evaluate flood risk. Analysis of the progression of LULC change over the historic period and associated SWAT outputs revealed that AMF magnitudes tend to increase over time in response to increasing degrees of urbanization. This is consistent with positive trends in the AMF series identified in previous studies, although there are difficulties identifying correlations between LULC change and identified change points due to large time gaps in the generated historical LULC maps, mainly caused by unavailability of sufficient quality historic aerial imagery. Similarly, increases in the mean and median AMF magnitude were observed in response to future LULC change projections, with the tails of the distributions remaining reasonably constant. FORE-SCE scenario A2 was found to have the most dramatic impact on AMF series, consistent with more extreme projections of population growth, demands for growing energy sources, agricultural land, and urban expansion, while AMF outputs based on scenario B2 showed little changes for the future as the focus is on environmental conservation and regional solutions to environmental issues.

EVALUATION OF NON-STATIONARITY IN ANNUAL MAXIMUM FLOOD SERIES OF MODERATELY IMPAIRED WATERSHEDS IN THE UPPER MIDWEST AND NORTHEASTERN UNITED STATES

United States federal agencies assess flood risk using Bulletin 17B procedures which assume annual maximum flood series are stationary. This represents a significant limitation of current flood frequency models as the flood distribution is thereby assumed to be unaffected by trends or periodicity of atmospheric/climatic variables and/or anthropogenic activities. The validity of this assumption is at the core of this thesis, which aims to improve understanding of the forms and potential causes of non-stationarity in flood series for moderately impaired watersheds in the Upper Midwest and Northeastern US. Prior studies investigated non-stationarity in flood series for unimpaired watersheds; however, as the majority of streams are located in areas of increasing human activity, relative and coupled impacts of natural and anthropogenic factors need to be considered such that non-stationary flood frequency models can be developed for flood risk forecasting over relevant planning horizons for large scale water resources planning and management.