Processing and mechanical properties of cast aluminum containing scandium, zirconium, and ytterbium

A series of aluminum alloys containing additions of scandium, zirconium, and ytterbium were cast to evaluate the effect of partial ytterbium substitution for scandium on tensile behavior. Due to the high price of scandium, a crucible-melt interaction study was performed to ensure no scandium was lost in graphite, alumina, magnesia, or zirconia crucibles after holding a liquid Al-Sc master alloy for 8 hours at 900 °C in an argon atmosphere. The alloys were subjected to an isochronal aging treatment and tested for conductivity and Vickers microhardness after each increment. For scandium-containing alloys, peak hardnesses of 520-790 MPa, and peak tensile stresses of 138-234 MPa were observed after aging from 150-350 °C for 3 hours in increments of 50 °C, and for alloys without scandium, peak hardnesses of 217-335 MPa and peak tensile stresses of 45-63 MPa were observed after a 3 hour, 150 °C aging treatment. The hardness and tensile strength of the ytterbium containing alloy was found to be lower than in the alloy with no ytterbium substitution.

Torque converter turbine noise and cavitation noise over varying speed ratio

These investigations will discuss the operational noise caused by automotive torque converters during speed ratio operation. Two specific cases of torque converter noise will be studied; cavitation, and a monotonic turbine induced noise. Cavitation occurs at or near stall, or zero turbine speed. The bubbles produced due to the extreme torques at low speed ratio operation, upon collapse, may cause a broadband noise that is unwanted by those who are occupying the vehicle as other portions of the vehicle drive train improve acoustically. Turbine induced noise, which occurs at high engine torque at around 0.5 speed ratio, is a narrow-band phenomenon that is audible to vehicle occupants currently. The solution to the turbine induced noise is known, however this study is to gain a better understanding of the mechanics behind this occurrence. The automated torque converter dynamometer test cell was utilized in these experiments to determine the effect of torque converter design parameters on the offset of cavitation and to employ the use a microwave telemetry system to directly measure pressures and structural motion on the turbine. Nearfield acoustics were used as a detection method for all phenomena while using a standardized speed ratio sweep test. Changes in filtered sound pressure levels enabled the ability to detect cavitation desinence. This, in turn, was utilized to determine the effects of various torque converter design parameters, including diameter, torus dimensions, and pump and stator blade designs on cavitation. The on turbine pressures and motion measured with the microwave telemetry were used to understand better the effects of a notched trailing edge turbine blade on the turbine induced noise.

Quantum modeling of bioconjugated nanomaterials

In recent years, the bio-conjugated nanostructured materials have emerged as a new class of materials for the bio-sensing and medical diagnostics applications. In spite of their multi-directional applications, interfacing nanomaterials with bio-molecules has been a challenge due to somewhat limited knowledge about the underlying physics and chemistry behind these interactions and also for the complexity of biomolecules. The main objective of this dissertation is to provide such a detailed knowledge on bioconjugated nanomaterials toward their applications in designing the next generation of sensing devices. Specifically, we investigate the changes in the electronic properties of a boron nitride nanotube (BNNT) due to the adsorption of different bio-molecules, ranging from neutral (DNA/RNA nucleobases) to polar (amino acid molecules). BNNT is a typical member of III-V compounds semiconductors with morphology similar to that of carbon nanotubes (CNTs) but with its own distinct properties. More specifically, the natural affinity of BNNTs toward living cells with no apparent toxicity instigates the applications of BNNTs in drug delivery and cell therapy. Our results predict that the adsorption of DNA/RNA nucleobases on BNNTs amounts to different degrees of modulation in the band gap of BNNTs, which can be exploited for distinguishing these nucleobases from each other. Interestingly, for the polar amino acid molecules, the nature of interaction appeared to vary ranging from Coulombic, van der Waals and covalent depending on the polarity of the individual molecules, each with a different binding strength and amount of charge transfer involved in the interaction. The strong binding of amino acid molecules on the BNNTs explains the observed protein wrapping onto BNNTs without any linkers, unlike carbon nanotubes (CNTs). Additionally, the widely varying binding energies corresponding to different amino acid molecules toward BNNTs indicate to the suitability of BNNTs for the biosensing applications, as compared to the metallic CNTs. The calculated I-V characteristics in these bioconjugated nanotubes predict notable changes in the conductivity of BNNTs due to the physisorption of DNA/RNA nucleobases. This is not the case with metallic CNTs whose transport properties remained unaltered in their conjugated systems with the nucleobases. Collectively, the bioconjugated BNNTs are found to be an excellent system for the next generation sensing devices.

Structural and surface property characterization of titanium dioxide nanotubes for orthopedic implants

This research focused on the to modification of the surface structure of titanium implants with nanostructured morphology of TiO2 nanotubes and studied the interaction of nanotubes with osteoblast cells to understand the parameters that affect the cell growth. The electrical, mechanical, and structural properties of TiO2 nanotubes were characterized to establish a better understanding on the properties of such nanoscale morphological structures. To achieve the objectives of this research work I transformed the titanium and its alloys, either in bulk sheet form, bulk machined form, or thin film deposited on another substrate into a surface of titania nanotubes using a low cost and environmentally friendly process. The process requires only a simple electrolyte, low cost electrode, and a DC power supply. With this simple approach of scalable nanofabrication, a typical result is nanotubes that are each approximately 100nm in diameter and have a wall thickness of about 20nm. By changing the fabrication parameters, independent nanotubes can be fabricated with open volume between them. Titanium in this form is termed onedimensional since electron transport is narrowly confined along the length of the nanotube. My Ph.D. accomplishments have successfully shown that osteoblast cells, the cells that are the precursors to bone, have a strong tendency to attach to the inside and outside of the titanium nanotubes onto which they are grown using their filopodia – cell’s foot used for locomotion – anchored to titanium nanotubes. In fact it was shown that the cell prefers to find many anchoring sites. These sites are critical for cell locomotion during the first several weeks of maturity and upon calcification as a strongly anchored bone cell. In addition I have shown that such a surface has a greater cell density than a smooth titanium surface. My work also developed a process that uses a focused and controllably rastered ion beam as a nano-scalpel to cut away sections of the osteoblast cells to probe the attachment beneath the main cell body. Ultimately the more rapid growth of osteoblasts, coupled with a stronger cell-surface interface, could provide cost reduction, shorter rehabilitation, and fewer follow-on surgeries due to implant loosening.

Tropospheric ozone and CO over North Atlantic for the past decade

We investigate how declines in US emissions of CO and O3 precursors have impacted the lower free troposphere over the North Atlantic. We use seasonal observations for O3 and CO from the PICO-NARE project for the period covering 2001 to 2010. Observations are used to verify model output generated by the GEOS-Chem 3-D global chemical transport model. Additional satellite data for CO from AIRS/Aqua and for O3 from TES/Aura were also used to provide additional comparisons; particularly for fall, winter, and spring when PICO-NARE coverage is sparse. We find GEOS-Chem captures the seasonal cycle for CO and O3 well compared to PICO-NARE data. For CO, GEOS-Chem is biased low, particularly in spring which is in agreement with findings from previous studies. GEOS-Chem is 24.7 +/- 5.2 ppbv (1-σ) low compared to PICO-NARE summer CO data while AIRS is 14.2 +/- 6.6 ppbv high. AIRS does not show nearly as much variation as seen with GEOS-Chem or the Pico data, and goes from being lower than PICO-NARE data in winter and spring, to higher in summer and fall. Both TES and GEOS-Chem match the seasonal ozone cycle well for all seasons when compared with observations. Model results for O3 show GEOS-Chem is 6.67 +/- 2.63 ppbv high compared to PICO-NARE summer measurements and TES was 3.91 +/- 4.2 ppbv higher. Pico data, model results, and AIRS all show declines in CO and O3 for the summer period from 2001 to 2010. Limited availability of TES data prevents us from using it in trend analysis. For summer CO Pico, GEOS-Chem, and AIRS results show declines of 1.32, 0.368, and 0.548 ppbv/year respectively. For summer O3, Pico and GEOS-Chem show declines of -0.726 and -0.583 ppbv/year respectively. In other seasons, both model and AIRS show declining CO, particularly in the fall. GEOS-Chem results show a fall decline of 0.798 ppbv/year and AIRS shows a decline of 0.8372 ppbv/year. Winter and spring CO declines are 0.393 and 0.307 for GEOS-Chem, and 0.455 and 0.566 for AIRS. GEOS-Chem shows declining O3 in other seasons as well; with fall being the season of greatest decrease and winter being the least. Model results for fall, winter, and spring are 0.856, 0.117, and 0.570 ppbv/year respectively. Given the availability of data we are most confident in summer results and thus find that summer CO and O3 have declined in lower free troposphere of the North Atlantic region of the Azores. Sensitivity studies for CO and O3 at Pico were conducted by turning off North American fossil fuel emissions in GEOS-Chem. Model results show that North America fossil fuel emissions contribute 8.57 ppbv CO and 4.03 ppbv O3 to Pico. The magnitude of modeled trends declines in all seasons without North American fossil fuel emissions except for summer CO. The increase in summer CO declines may be due to a decline of 5.24 ppbv/year trend in biomass burning emissions over the study period; this is higher than the 2.33 ppbv/year North American anthropogenic CO model decline. Winter O3 is the only season which goes from showing a negative trend to a positive trend.

Muffler characterization with implementation of the finite element method and experimental techniques

Determining how an exhaust system will perform acoustically before a prototype muffler is built can save the designer both a substantial amount of time and resources. In order to effectively use the simulation tools available it is important to understand what is the most effective tool for the intended purpose of analysis as well as how typical elements in an exhaust system affect muffler performance. An in-depth look at the available tools and their most beneficial uses are presented in this thesis. A full parametric study was conducted using the FEM method for typical muffler elements which was also correlated to experimental results. This thesis lays out the overall ground work on how to accurately predict sound pressure levels in the free field for an exhaust system with the engine properties included. The accuracy of the model is heavily dependent on the correct temperature profile of the model in addition to the accuracy of the source properties. These factors will be discussed in detail and methods for determining them will be presented. The secondary effects of mean flow, which affects both the acoustical wave propagation and the flow noise generation, will be discussed. Effective ways for predicting these secondary effects will be described. Experimental models will be tested on a flow rig that showcases these phenomena.

Balancing wildlife and bioenergy values on private forestlands in the eastern Upper Peninsula of Michigan

There is interest in developing a reliable, sustainable, domestic U.S. biofuels industry. A domestic biofuels industry has the potential to provide economic, environmental, and national security benefits on a local, regional, national, and global level. The Mascoma Corporation plans to develop a cellulosic ethanol facility in Michigan’s eastern Upper Peninsula. The primary feedstock of the plant site would be trees sourced within a 150 mile supply radius. In the eastern Upper Peninsula, this radius encompasses Alger, Chippewa, Delta, Luce, Mackinac, and Schoolcraft counties. In these six counties there are 1,320,500 acres of NIPF (non-industrial private forestlands). These acres account for 40% of the total timberland in these six counties. Thus it is likely that in order for the successful implementation of a cellulosic ethanol facility the support of local NIPF owners will be necessary. This thesis presents research on how eastern Upper Peninsula forest landowners think about and manage their land. It is based on 48 in-depth interviews with these landowners. The goal was to determine how landowner values and beliefs, on a variety of issues including wildlife management, land management, biofuels development, and climate change, are expressed through both their current management decisions, and possibly their future land management decisions. Some of the values articulated by the landowners in this study included biodiversity protection, conservation of healthy game populations, and the production of high-value timber. Understanding the values and beliefs of landowners in the eastern Upper Peninsula of Michigan is critical for successfully developing a sustainable regional woody bioenergy.

Development of a high performance parallel computing platform and its use in the study of nanostructures : clusters, sheets and tubes .

Small clusters of gallium oxide, technologically important high temperature ceramic, together with interaction of nucleic acid bases with graphene and small-diameter carbon nanotube are focus of first principles calculations in this work. A high performance parallel computing platform is also developed to perform these calculations at Michigan Tech. First principles calculations are based on density functional theory employing either local density or gradient-corrected approximation together with plane wave and gaussian basis sets. The bulk Ga2O3 is known to be a very good candidate for fabricating electronic devices that operate at high temperatures. To explore the properties of Ga2O3 at nonoscale, we have performed a systematic theoretical study on the small polyatomic gallium oxide clusters. The calculated results find that all lowest energy isomers of GamOn clusters are dominated by the Ga-O bonds over the metal-metal or the oxygen-oxygen bonds. Analysis of atomic charges suggest the clusters to be highly ionic similar to the case of bulk Ga2O3. In the study of sequential oxidation of these slusters starting from Ga2O, it is found that the most stable isomers display up to four different backbones of constituent atoms. Furthermore, the predicted configuration of the ground state of Ga2O is recently confirmed by the experimental result of Neumark's group. Guided by the results of calculations the study of gallium oxide clusters, performance related challenge of computational simulations, of producing high performance computers/platforms, has been addressed. Several engineering aspects were thoroughly studied during the design, development and implementation of the high performance parallel computing platform, rama, at Michigan Tech. In an attempt to stay true to the principles of Beowulf revolutioni, the rama cluster was extensively customized to make it easy to understand, and use - for administrators as well as end-users. Following the results of benchmark calculations and to keep up with the complexity of systems under study, rama has been expanded to a total of sixty four processors. Interest in the non-covalent intereaction of DNA with carbon nanotubes has steadily increased during past several years. This hybrid system, at the junction of the biological regime and the nanomaterials world, possesses features which make it very attractive for a wide range of applicatioins. Using the in-house computational power available, we have studied details of the interaction between nucleic acid bases with graphene sheet as well as high-curvature small-diameter carbon nanotube. The calculated trend in the binding energies strongly suggests that the polarizability of the base molecules determines the interaction strength of the nucleic acid bases with graphene. When comparing the results obtained here for physisorption on the small diameter nanotube considered with those from the study on graphene, it is observed that the interaction strength of nucleic acid bases is smaller for the tube. Thus, these results show that the effect of introducing curvature is to reduce the binding energy. The binding energies for the two extreme cases of negligible curvature (i.e. flat graphene sheet) and of very high curvature (i.e. small diameter nanotube) may be considered as upper and lower bounds. This finding represents an important step towards a better understanding of experimentally observed sequence-dependent interaction of DNA with Carbon nanotubes.

Ultraviolet digital imaging of volcanic plumes : implementation and application to magmatic processes at basaltic volcanoes

Magmatic volatiles play a crucial role in volcanism, from magma production at depth to generation of seismic phenomena to control of eruption style. Accordingly, many models of volcano dynamics rely heavily on behavior of such volatiles. Yet measurements of emission rates of volcanic gases have historically been limited, which has restricted model verification to processes on the order of days or longer. UV cameras are a recent advancement in the field of remote sensing of volcanic SO2 emissions. They offer enhanced temporal and spatial resolution over previous measurement techniques, but need development before they can be widely adopted and achieve the promise of integration with other geophysical datasets. Large datasets require a means by which to quickly and efficiently use imagery to calculate emission rates. We present a suite of programs designed to semi-automatically determine emission rates of SO2 from series of UV images. Extraction of high temporal resolution SO2 emission rates via this software facilitates comparison of gas data to geophysical data for the purposes of evaluating models of volcanic activity and has already proven useful at several volcanoes. Integrated UV camera and seismic measurements recorded in January 2009 at Fuego volcano, Guatemala, provide new insight into the system’s shallow conduit processes. High temporal resolution SO2 data reveal patterns of SO2 emission rate relative to explosions and seismic tremor that indicate tremor and degassing share a common source process. Progressive decreases in emission rate appear to represent inhibition of gas loss from magma as a result of rheological stiffening in the upper conduit. Measurements of emission rate from two closely-spaced vents, made possible by the high spatial resolution of the camera, help constrain this model. UV camera measurements at Kilauea volcano, Hawaii, in May of 2010 captured two occurrences of lava filling and draining within the summit vent. Accompanying high lava stands were diminished SO2 emission rates, decreased seismic and infrasonic tremor, minor deflation, and slowed lava lake surface velocity. Incorporation of UV camera data into the multi-parameter dataset gives credence to the likelihood of shallow gas accumulation as the cause of such events.

Teaching and aligning upper level mathematics classes to Michigan high school content expectations in an alternative education setting : approaches and issues

After teaching regular education secondary mathematics for seven years, I accepted a position in an alternative education high school. Over the next four years, the State of Michigan adopted new graduation requirements phasing in a mandate for all students to complete Geometry and Algebra 2 courses. Since many of my students were already struggling in Algebra 1, getting them through Geometry and Algebra 2 seemed like a daunting task. To better instruct my students, I wanted to know how other teachers in similar situations were addressing the new High School Content Expectations (HSCEs) in upper level mathematics. This study examines how thoroughly alternative education teachers in Michigan are addressing the HSCEs in their courses, what approaches they have found most effective, and what issues are preventing teachers and schools from successfully implementing the HSCEs. Twenty-six alternative high school educators completed an online survey that included a variety of questions regarding school characteristics, curriculum alignment, implementation approaches and issues. Follow-up phone interviews were conducted with four of these participants. The survey responses were used to categorize schools as successful, unsuccessful, and neutral schools in terms of meeting the HSCEs. Responses from schools in each category were compared to identify common approaches and issues among them and to identify significant differences between school groups. Data analysis showed that successful schools taught more of the HSCEs through a variety of instructional approaches, with an emphasis on varying the ways students learned the material. Individualized instruction was frequently mentioned by successful schools and was strikingly absent from unsuccessful school responses. The main obstacle to successful implementation of the HSCEs identified in the study was gaps in student knowledge. This caused pace of instruction to also be a significant issue. School representatives were fairly united against the belief that the Algebra 2 graduation requirement was appropriate for all alternative education students. Possible implications of these findings are discussed.

Digital processing and data compilation approach for using remotely sensed imagery to identify geological lineaments in hard-rock terrains : an application for groundwater explorations in Nicaragua

Sustainable yields from water wells in hard-rock aquifers are achieved when the well bore intersects fracture networks. Fracture networks are often not readily discernable at the surface. Lineament analysis using remotely sensed satellite imagery has been employed to identify surface expressions of fracturing, and a variety of image-analysis techniques have been successfully applied in “ideal” settings. An ideal setting for lineament detection is where the influences of human development, vegetation, and climatic situations are minimal and hydrogeological conditions and geologic structure are known. There is not yet a well-accepted protocol for mapping lineaments nor have different approaches been compared in non-ideal settings. A new approach for image-processing/synthesis was developed to identify successful satellite imagery types for lineament analysis in non-ideal terrain. Four satellite sensors (ASTER, Landsat7 ETM+, QuickBird, RADARSAT-1) and a digital elevation model were evaluated for lineament analysis in Boaco, Nicaragua, where the landscape is subject to varied vegetative cover, a plethora of anthropogenic features, and frequent cloud cover that limit the availability of optical satellite data. A variety of digital image processing techniques were employed and lineament interpretations were performed to obtain 12 complementary image products that were evaluated subjectively to identify lineaments. The 12 lineament interpretations were synthesized to create a raster image of lineament zone coincidence that shows the level of agreement among the 12 interpretations. A composite lineament interpretation was made using the coincidence raster to restrict lineament observations to areas where multiple interpretations (at least 4) agree. Nine of the 11 previously mapped faults were identified from the coincidence raster. An additional 26 lineaments were identified from the coincidence raster, and the locations of 10 were confirmed by field observation. Four manual pumping tests suggest that well productivity is higher for wells proximal to lineament features. Interpretations from RADARSAT-1 products were superior to interpretations from other sensor products, suggesting that quality lineament interpretation in this region requires anthropogenic features to be minimized and topographic expressions to be maximized. The approach developed in this study has the potential to improve siting wells in non-ideal regions.

Cladophora, mass transport, and the nearshore phosphorus shunt

Excessive Cladophora growth in the Great Lakes has led to beach fouling and the temporary closure of nuclear power plants and has been associated with avian botulism and the persistence of human pathogens. As the growth-limiting nutrient for Cladophora, phosphorus is the appropriate target for management efforts. Dreissenids (zebra and quagga mussels) have the ability to capture particulate phase phosphorus (otherwise unavailable to Cladophora) and release it in a soluble, available form. The significance of this potential nutrient source is, in part, influenced by the interplay between phosphorus flux from the mussel bed and turbulent mixing in establishing the phosphorus levels to which Cladophora is exposed. It is hypothesized that under quiescent conditions phosphorus will accumulate near the sediment-water interface, setting up vertical phosphorus gradients and favorable conditions for resource delivery to Cladophora. These gradients would be eliminated under conditions of wind mixing, reducing the significance of the dreissenid-mediated nutrient contribution. Soluble reactive phosphorus (SRP) levels were monitored over dreissenid beds (densities on the order of 350•m-2 and 3000∙m-2) at a site 8 m deep in Lake Michigan. Monitoring was based on the deployment of Modified Hesslein Samplers which collected samples for SRP analysis over a distance of 34 cm above the bottom in 2.5 cm intervals. Deployment intervals were established to capture a wind regime (calm, windy) that persisted for an interval consistent with the sampler equilibration time of 7 hours. Results indicate that increased mussel density leads to an increased concentration boundary layer; increased wind speed leads to entrainment of the concentration boundary layer; and increased duration of quiescent periods leads to an increased concentration boundary layer. This concentration boundary layer is of ecological significance and forms in the region inhabited by Cladophora

Erupted magma volume estimates at Santiaguito and Pacaya Volcanoes, Guatemala using digital elevation models

High resolution digital elevation models (DEMs) of Santiaguito and Pacaya volcanoes, Guatemala, were used to estimate volume changes and eruption rates between 1954 and 2001. The DEMs were generated from contour maps and aerial photography, which were analyzed in ArcGIS 9.0®. Because both volcanoes were growing substantially over the five decade period, they provide a good data set for exploring effective methodology for estimating volume changes. The analysis shows that the Santiaguito dome complex grew by 0.78 ± 0.07 km3 (0.52 ± 0.05 m3 s-1) over the 1954-2001 period with nearly all the growth occurring on the El Brujo (1958-75) and Caliente domes (1971-2001). Adding information from field data prior to 1954, the total volume extruded from Santiaguito since 1922 is estimated at 1.48 ± 0.19 km3. Santiaguito’s growth rate is lower than most other volcanic domes, but it has been sustained over a much longer period and has undergone a change toward more exogenous and progressively slower extrusion with time. At Santiaguito some of the material being added at the dome is subsequently transported downstream by block and ash flows, mudflows and floods, creating channel shifting and areas of aggradation and erosion. At Pacaya volcano a total volume of 0.21 ± 0.05 km3 was erupted between 1961 and 2001 for an average extrusion rate of 0.17 ± 0.04 m3 s-1. Both the Santiaguito and Pacaya eruption rate estimates reported here are minima, because they do not include estimates of materials which are transported downslope after eruption and data on ashfall which may result in significant volumes of material spread over broad areas. Regular analysis of high resolution DEMs using the methods outlined here, would help quantify the effects of fluvial changes to downstream populated areas, as well as assist in tracking hazards related to dome collapse and eruption.

Investigations into the degassing and eruption mechanisms of Nyamuragira volcano, Democratic Republic of the Congo (Africa)

One of two active volcanoes in the western branch of the East African Rift, Nyamuragira (1.408ºS, 29.20ºE; 3058 m) is located in the D.R. Congo. Nyamuragira emits large amounts of SO2 (up to ~1 Mt/day) and erupts low-silica, alkalic lavas, which achieve flow rates of up to ~20 km/hr. The source of the large SO2 emissions and pre-eruptive magma conditions were unknown prior to this study, and 1994-2010 lava volumes were only recently mapped via satellite imagery, mainly due to the region’s political instability. In this study, new olivine-hosted melt inclusion volatile (H2O, CO2, S, Cl, F) and major element data from five historic Nyamuragira eruptions (1912, 1938, 1948, 1986, 2006) are presented. Melt compositions derived from the 1986 and 2006 tephra samples best represent pre-eruptive volatile compositions because these samples contain naturally glassy inclusions that underwent less post-entrapment modification than crystallized inclusions. The total amount of SO2 released from the 1986 (0.04 Mt) and 2006 (0.06 Mt) eruptions are derived using the petrologic method, whereby S contents in melt inclusions are scaled to erupted lava volumes. These amounts are significantly less than satellite-based SO2 emissions for the same eruptions (1986 = ~1 Mt; 2006 = ~2 Mt). Potential explanations for this observation are: 1) accumulation of a vapor phase within the magmatic system that is only released during eruptions, and/or 2) syn-eruptive gas release from unerupted magma. Post-1994 Nyamuragira lava volumes were not available at the beginning of this study. These flows (along with others since 1967) are mapped with Landsat MSS, TM, and ETM+, Hyperion, and ALI satellite data and combined with published flow thicknesses to derive volumes. Satellite remote sensing data was also used to evaluate Nyamuragira SO2 emissions. These results show that the most recent Nyamuragira eruptions injected SO2 into the atmosphere between 15 km (2006 eruption) and 5 km (2010 eruption). This suggests that past effusive basaltic eruptions (e.g., Laki 1783) are capable of similar plume heights that reached the upper troposphere or tropopause, allowing SO2 and resultant aerosols to remain longer in the atmosphere, travel farther around the globe, and affect global climates.

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.