Hydrogeochemical controls on uranium in aquifers of the Jacobsville Sandstone

The occurrence of elevated uranium (U) in sandstone aquifers was investigated in the Upper Peninsula of Michigan, focusing on aquifers of the Jacobsville Sandstone. The hydrogeochemical controls on groundwater U concentrations were characterized using a combination of water sampling and spectral gamma-ray logging of sandstone cliffs and residential water wells. 235U/238U isotope ratios were consistent with naturally occurring U. Approximately 25% of the 270 wells tested had U concentrations above the U.S. Environmental Protection Agency Maximum Contaminant Level (MCL) of 30 μg/L, with elevated U generally occurring in localized clusters. Water wells were logged to determine whether groundwater U anomalies could be explained by the heterogeneous distribution of U in the sandstone. Not all wells with relative U enrichment in the sandstone produced water with U above the MCL, indicating that the effect of U enrichment in the sandstone may be modified by other hydrogeochemical factors. Well water had high redox, indicating U is in its highly soluble (VI) valence. Equilibrium modeling indicated that aqueous U is complexed with carbonates. In general, wells with elevated U concentrations had low 235U/238U activity ratios. However, in some areas U concentrations and 235U/238U activity ratios were simultaneously high, possibly indicating differences in rock-water interactions. Limited groundwater age dating suggested that residence time may also help explain variations in well water U concentrations. Low levels of U enrichment (4 to 30 ppm) in the Jacobsville sandstone may make wells in its oxidized aquifers at risk for U concentrations above the MCL. On average, U concentrations were highest in heavy mineral and clay layers and rip up conglomerates. Uranium concentrations above 4 ppm also occurred in siltstones, sandstones and conglomerates. Uranium enrichment was likely controlled by deposition processes, sorption to clays, and groundwater flow, which was controlled by permeability variations in the sandstone. Low levels of U enrichment were found at deltaic, lacustrine and alluvial fan deposits and were not isolated to specific depositional environments.

Automated optimization of polymer extrusion dies using finite element analysis

An extrusion die is used to continuously produce parts with a constant cross section; such as sheets, pipes, tire components and more complex shapes such as window seals. The die is fed by a screw extruder when polymers are used. The extruder melts, mixes and pressures the material by the rotation of either a single or double screw. The polymer can then be continuously forced through the die producing a long part in the shape of the die outlet. The extruded section is then cut to the desired length. Generally, the primary target of a well designed die is to produce a uniform outlet velocity without excessively raising the pressure required to extrude the polymer through the die. Other properties such as temperature uniformity and residence time are also important but are not directly considered in this work. Designing dies for optimal outlet velocity variation using simple analytical equations are feasible for basic die geometries or simple channels. Due to the complexity of die geometry and of polymer material properties design of complex dies by analytical methods is difficult. For complex dies iterative methods must be used to optimize dies. An automated iterative method is desired for die optimization. To automate the design and optimization of an extrusion die two issues must be dealt with. The first is how to generate a new mesh for each iteration. In this work, this is approached by modifying a Parasolid file that describes a CAD part. This file is then used in a commercial meshing software. Skewing the initial mesh to produce a new geometry was also employed as a second option. The second issue is an optimization problem with the presence of noise stemming from variations in the mesh and cumulative truncation errors. In this work a simplex method and a modified trust region method were employed for automated optimization of die geometries. For the trust region a discreet derivative and a BFGS Hessian approximation were used. To deal with the noise in the function the trust region method was modified to automatically adjust the discreet derivative step size and the trust region based on changes in noise and function contour. Generally uniformity of velocity at exit of the extrusion die can be improved by increasing resistance across the die but this is limited by the pressure capabilities of the extruder. In optimization, a penalty factor that increases exponentially from the pressure limit is applied. This penalty can be applied in two different ways; the first only to the designs which exceed the pressure limit, the second to both designs above and below the pressure limit. Both of these methods were tested and compared in this work.

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.

SELECTED CHEMOAUTOTROPHIC PROCESSES IN WASTEWATER TREATMENT: LOW TEMPERATURE NITRIFICATION INHIBITION BY AN AZO DYE AND BIOFILTRATION FOR ODOR CONTROL OF HYDROGEN SULFIDE

Biochemical processes by chemoautotrophs such as nitrifiers and sulfide and iron oxidizers are used extensively in wastewater treatment. The research described in this dissertation involved the study of two selected biological processes utilized in wastewater treatment mediated by chemoautotrophic bacteria: nitrification (biological removal of ammonia and nitrogen) and hydrogen sulfide (H2S) removal from odorous air using biofiltration. A municipal wastewater treatment plant (WWTP) receiving industrial dyeing discharge containing the azo dye, acid black 1 (AB1) failed to meet discharge limits, especially during the winter. Dyeing discharge mixed with domestic sewage was fed to sequencing batch reactors at 22oC and 7oC. Complete nitrification failure occurred at 7oC with more rapid nitrification failure as the dye concentration increased; slight nitrification inhibition occurred at 22oC. Dye-bearing wastewater reduced chemical oxygen demand (COD) removal at 7oC and 22oC, increased i effluent total suspended solids (TSS) at 7oC, and reduced activated sludge quality at 7oC. Decreasing AB1 loading resulted in partial nitrification recovery. Eliminating the dye-bearing discharge to the full-scale WWTP led to improved performance bringing the WWTP into regulatory compliance. BiofilterTM, a dynamic model describing the biofiltration processes for hydrogen sulfide removal from odorous air emissions, was calibrated and validated using pilot- and full-scale biofilter data. In addition, the model predicted the trend of the measured data under field conditions of changing input concentration and low effluent concentrations. The model demonstrated that increasing gas residence time and temperature and decreasing influent concentration decreases effluent concentration. Model simulations also showed that longer residence times are required to treat loading spikes. BiofilterTM was also used in the preliminary design of a full-scale biofilter for the removal of H2S from odorous air. Model simulations illustrated that plots of effluent concentration as a function of residence time or bed area were useful to characterize and design biofilters. Also, decreasing temperature significantly increased the effluent concentration. Model simulations showed that at a given temperature, a biofilter cannot reduce H2S emissions below a minimum value, no matter how large the biofilter.

EXPLORATION OF THE MTSAT2 SATELLITE CAPABILITIES FOR REAL TIME DETECTION AND CHARACTERIZATION OF VOLCANIC EMISSIONS

In this report, we attempt to define the capabilities of the infrared satellite remote sensor, Multifunctional Transport Satellite-2 (MTSAT-2) (i.e. a geosynchronous instrument), in characterizing volcanic eruptive behavior in the highly active region of Indonesia. Sulfur dioxide data from NASA's Ozone Monitoring Instrument (OMI) (i.e. a polar orbiting instrument) are presented here for validation of the processes interpreted using the thermal infrared datasets. Data provided from two case studies are analyzed specifically for eruptive products producing large thermal anomalies (i.e. lava flows, lava domes, etc.), volcanic ash and SO2 clouds; three distinctly characteristic and abundant volcanic emissions. Two primary methods used for detection of heat signatures are used and compared in this report including, single-channel thermal radiance (4-µm) and the normalized thermal index (NTI) algorithm. For automated purposes, fixed thresholds must be determined for these methods. A base minimum detection limit (MDL) for single-channel thermal radiance of 2.30E+05 Wm- 2sr-1m-1 and -0.925 for NTI generate false alarm rates of 35.78% and 34.16%, respectively. A spatial comparison method, developed here specifically for use in Indonesia and used as a second parameter for detection, is implemented to address the high false alarm rate. For the single-channel thermal radiance method, the utilization of the spatial comparison method eliminated 100% of the false alarms while maintaining every true anomaly. The NTI algorithm showed similar results with only 2 false alarms remaining. No definitive difference is observed between the two thermal detection methods for automated use; however, the single-channel thermal radiance method coupled with the SO2 mass abundance data can be used to interpret volcanic processes including the identification of lava dome activity at Sinabung as well as the mechanism for the dome emplacement (i.e. endogenous or exogenous). Only one technique, the brightness temperature difference (BTD) method, is used for the detection of ash. Trends of ash area, water/ice area, and their respective concentrations yield interpretations of increased ice formation, aggregation, and sedimentation processes that only a high-temporal resolution instrument like the MTSAT-2 can analyze. A conceptual model of a secondary zone of aggregation occurring in the migrating Kelut ash cloud, which decreases the distal fine-ash component and hazards to flight paths, is presented in this report. Unfortunately, SO2 data was unable to definitively reinforce the concept of a secondary zone of aggregation due to the lack of a sufficient temporal resolution. However, a detailed study of the Kelut SO2 cloud is used to determine that there was no climatic impacts generated from this eruption due to the atmospheric residence times and e-folding rate of ~14 days for the SO2. This report applies the complementary assets offered by utilizing a high-temporal and a high-spatial resolution satellite, and it demonstrates that these two instruments can provide unparalleled observations of dynamic volcanic processes.