Automation of shear-wave splitting parameter determination of local earthquakes at Yellowstone : application as indicator of crustal stress and temporal variation

Shear-wave splitting can be a useful technique for determining crustal stress fields in volcanic settings and temporal variations associated with activity. Splitting parameters were determined for a subset of local earthquakes recorded from 2000-2010 at Yellowstone. Analysis was automated using an unsupervised cluster analysis technique to determine optimum splitting parameters from 270 analysis windows for each event. Six stations clearly exhibit preferential fast polarization values sub-orthogonal to the direction of minimum horizontal compression. Yellowstone deformation results in a local crustal stress field differing from the regional field dominated by NE-SW extension, and fast directions reflect this difference rotating around the caldera maintaining perpendicularity to the rim. One station exhibits temporal variations concordant with identified periods of caldera subsidence and uplift. From splitting measurements, we calculated a crustal anisotropy of ~17-23% and crack density ~0.12-0.17 possibly resulting from stress-aligned fluid filled microcracks in the upper crust and an active hydrothermal system.

Estimation of vertical groundwater fluxes into a streambed through continuous temperature profile monitoring and the relationship of groundwater fluxes to coaster brook trout spawning habitat

We hypothesized that the spatial distribution of groundwater inflows through river bottom sediments is a critical factor associated with the selection of coaster brook trout (a life history variant of Salvelinus fontinalis,) spawning sites. An 80-m reach of the Salmon Trout River, in the Huron Mountains of the upper peninsula of Michigan, was selected to test the hypothesis based on long-term documentation of coaster brook trout spawning at this site. Throughout this site, the river is relatively similar along its length with regard to stream channel and substrate features. A monitoring well system consisting of an array of 27 wells was installed to measure subsurface temperatures underneath the riverbed over a 13-month period. The monitoring well locations were separated into areas where spawning has and has not been observed. Over 200,000 total temperature measurements were collected from 5 depths within each of the 27 monitoring wells. Temperatures within the substrate at the spawning area were generally cooler and less variable than river temperatures. Substrate temperatures in the non-spawning area were generally warmer, more variable, and closely tracked temporal variations in river temperatures. Temperature data were inverted to obtain subsurface groundwater velocities using a numerical approximation of the heat transfer equation. Approximately 45,000 estimates of groundwater velocities were obtained. Estimated velocities in the spawning and non-spawning areas confirmed that groundwater velocities in the spawning area were primarily in the upward direction, and were generally greater in magnitude than velocities in the non-spawning area. In the non-spawning area there was a greater occurrence of velocities in the downward direction, and velocity estimates were generally lesser in magnitude than in the spawning area. Both the temperature and velocity results confirm the hypothesis that spawning sites correspond to areas of significant groundwater influx to the river bed.

EVALUATION OF THE EVOLVING STRESS FIELD OF THE YELLOWSTONE VOLCANIC PLATEAU FROM 1988 TO 2010 FROM EARTHQUAKE FIRST MOTIONS INVERSION

Within the Yellowstone National Park, Wyoming, the silicic Yellowstone volcanic field is one of the most active volcanic systems all over the world. Although the last rhyolite eruption occurred around 70,000 years ago, Yellowstone is still believed to be volcanically active, due to high hydrothermal and seismic activity. The earthquake data used in this study cover the period of time between 1988 and 2010. Earthquake relocations and a set of 369 well-constrained, double-couple, focal mechanism solutions were computed. Events were grouped according to location and time to investigate trends in faulting. The majority of the events has oblique, normal-faulting solutions. The overall direction of extension throughout the 0.64 Ma Yellowstone caldera looks nearly ENE, consistently with the direction of alignments of volcanic vents within the caldera, but detailed study revealed spatial and temporal variations. Stress-field solutions for different areas and time periods were calculated from earthquake focal mechanism inversion. A well-resolved rotation of σ3 was found, from NNE-SSW near the Hebgen Lake fault zone, to ENE-WSW near Norris Junction. In particular, the σ3 direction changed throughout the years in the Norris Junction area, from being ENE-WSW, as calculated in the study by Waite and Smith (2004), to NNE-SSW, while the other σ3 directions are mostly unchanged over time. The Yellowstone caldera was subject to periods of net uplift and subsidence over the past century, explained in previous studies as caused by expanding or contracting sills, at different depths. Based on the models used to explain these deformation periods, we investigated the relationship between variability in aseismic deformation and seismic activity and faulting styles. Focal mechanisms and P and T axes were divided into temporal and depth intervals, in order to identify spatial or temporal trends in deformation. The presence of “chocolate tablet” structures, with composite dilational faults, was identified in many stages of the deformation history both in the Norris Geyser Basin area and inside the caldera. Strike-slip component movement was found in a depth interval below a contracting sill, indicating the movement of magma towards the caldera.

Effects of Chemical Deicers on Amphibian Communites

Roads and highways present a unique challenge to wildlife as they exhibit substantial impacts on the surrounding ecosystem through the interruption of a number of ecological processes. With new roads added to the national highway system every year, an understanding of these impacts is required for effective mitigation of potential environmental impacts. A major contributor to these negative effects is the deposition of chemicals used in winter deicing activities to nearby surface waters. These chemicals often vary in composition and may affect freshwater species differently. The negative impacts of widespread deposition of sodium chloride (NaCl) have prompted a search for an `environmentally friendly' alternative. However, little research has investigated the potential environmental effects of widespread use of these alternatives. Herein, I detail the results of laboratory tests and field surveys designed to determine the impacts of road salt (NaCl) and other chemical deicers on amphibian communities in Michigan's Upper Peninsula. Using larval amphibians I demonstrate the lethal impacts of a suite of chemical deicers on this sensitive, freshwater species. Larval wood frogs (Lithobates sylvatica) were tolerant of short-term (96 hours) exposure to urea (CH4N2O), sodium chloride (NaCl), and magnesium chloride (MgCl2). However, these larvae were very sensitive to acetate products (C8H12CaMgO8, CH3COOK) and calcium chloride (CaCl2). These differences in tolerance suggest that certain deicers may be more harmful to amphibians than others. Secondly, I expanded this analysis to include an experiment designed to determine the sublethal effects of chronic exposure to environmentally realistic concentrations of NaCl on two unique amphibian species, L. sylvatica and green frogs (L. clamitans). L. sylvatica tend to breed in small, ephemeral wetlands and metamorphose within a single season. However, L. clamitans breed primarily in more permanent wetlands and often remain as tadpoles for one year or more. These species employ different life history strategies in this region which may influence their response to chronic NaCl exposure. Both species demonstrated potentially harmful effects on individual fitness. L. sylvatica larvae had a high incidence of edema suggesting the NaCl exposure was a significant physiologic stressor to these larvae. L. clamitans larvae reduced tail length during their exposure which may affect adult fitness of these individuals. In order to determine the risk local amphibians face when using these roadside pools, I conducted a survey of the spatial distribution of chloride in the three northernmost counties of Michigan. This area receives a relatively low amount of NaCl which is confined to state and federal highways. The chloride concentrations in this region were much lower than those in urban systems; however, amphibians breeding in the local area may encounter harmful chloride levels arising from temporal variations in hydroperiods. Spatial variation of chloride levels suggests the road-effect zone for amphibians may be as large as 1000 m from a salt-treated highway. Lastly, I performed an analysis of the use of specific conductance to predict chloride concentrations in natural surface water bodies. A number of studies have used this regression to predict chloride concentrations from measurements of specific conductance. This method is often chosen in the place of ion chromatography due to budget and time constraints. However, using a regression method to characterize this relationship does not result in accurate chloride ion concentration estimates.

IMPACTS OF CLIMATE CHANGE ON SOIL MICROORGANISMS IN NORTHERN HARDWOOD FORESTS

As global climate continues to change, it becomes more important to understand possible feedbacks from soils to the climate system. This dissertation focuses on soil microbial community responses to climate change factors in northern hardwood forests. Two soil warming experiments at Harvard Forest in Massachusetts, and a climate change manipulation experiment with both elevated temperature and increased moisture inputs in Michigan were sampled. The hyphal in-growth bag method was to understand how soil fungal biomass and respiration respond to climate change factors. Our results from phospholipid fatty acid (PLFA) analyses suggest that the hyphal in-growth bag method allows relatively pure samples of fungal hyphae to be partitioned from bacteria in the soil. The contribution of fungal hyphal respiration to soil respiration was examined in climate change manipulation experiments in Massachusetts and Michigan. The Harvard Forest soil warming experiments in Massachusetts are long-term studies with 8 and 18 years of +5 °C warming treatment. Hyphal respiration and biomass production tended to decrease with soil warming at Harvard Forest. This suggests that fungal hyphae adjust to higher temperatures by decreasing the amount of carbon respired and the amount of carbon stored in biomass. The Ford Forestry Center experiment in Michigan has a 2 x 2 fully factorial design with warming (+4-5 °C) and moisture addition (+30% average ambient growing season precipitation). This experiment was used to examine hyphal growth and respiration of arbuscular mycorrhizal fungi (AMF), soil enzymatic capacity, microbial biomass and microbial community structure in the soil over two years of experimental treatment. Results from the hyphal in-growth bag study indicate that AMF hyphal growth and respiration respond negatively to drought. Soil enzyme activities tend to be higher in heated versus unheated soils. There were significant temporal variations in enzyme activity and microbial biomass estimates. When microbial biomass was estimated using chloroform fumigation extractions there were no differences between experimental treatments and the control. When PLFA analyses were used to estimate microbial biomass we found that biomass responds negatively to higher temperatures and positively to moisture addition. This pattern was present for both bacteria and fungi. More information on the quality and composition of the organic matter and nutrients in soils from climate change manipulation experiments will allow us to gain a more thorough understanding of the mechanisms driving the patterns reported here. The information presented here will improve current soil carbon and nitrogen cycling models.