Understanding volcanic processes using UV camera measurements of sulfur dioxide and coincident infrasound and seismicity

The exsolution of volatiles from magma maintains an important control on volcanic eruption styles. The nucleation, growth, and connectivity of bubbles during magma ascent provide the driving force behind eruptions, and the rate, volume, and ease of gas exsolution can affect eruptive activity. Volcanic plumes are the observable consequence of this magmatic degassing, and remote sensing techniques allow us to quantify changes in gas exsolution. However, until recently the methods used to measure volcanic plumes did not have the capability of detecting rapid changes in degassing on the scale of standard geophysical observations. The advent of the UV camera now makes high sample rate gas measurements possible. This type of dataset can then be compared to other volcanic observations to provide an in depth picture of degassing mechanisms in the shallow conduit. The goals of this research are to develop a robust methodology for UV camera field measurements of volcanic plumes, and utilize this data in conjunction with seismoacoustic records to illuminate degassing processes. Field and laboratory experiments were conducted to determine the effects of imaging conditions, vignetting, exposure time, calibration technique, and filter usage on the UV camera sulfur dioxide measurements. Using the best practices determined from these studies, a field campaign was undertaken at Volcán de Pacaya, Guatemala. Coincident plume sulfur dioxide measurements, acoustic recordings, and seismic observations were collected and analyzed jointly. The results provide insight into the small explosive features, variations in degassing rate, and plumbing system of this complex volcanic system. This research provides useful information for determining volcanic hazard at Pacaya, and demonstrates the potential of the UV camera in multiparameter studies.

Effectiveness of a nondestructive evaluation technique for assessing standing timber quality

The research presented in this thesis was conducted to further the development of the stress wave method of nondestructively assessing the quality of wood in standing trees. The specific objective of this research was to examine, in the field, use of two stress wave nondestructive assessment techniques. The first technique examined utilizes a laboratory-built measurement system consisting of commercially available accelerometers and a digital storage oscilloscope. The second technique uses a commercially available tool that incorporates several technologies to determine speed of stress wave propagation in standing trees. Field measurements using both techniques were conducted on sixty red pine trees in south-central Wisconsin and 115 ponderosa pine trees in western Idaho. After in-situ measurements were taken, thirty tested red pine trees were felled and a 15-foot-long butt log was obtained from each tree, while all tested ponderosa pine trees were felled and an 8 1/2 -foot-long butt log was obtained, respectively. The butt logs were sent to the USDA Forest Products Laboratory and nondestructively tested using a resonance stress wave technique. Strong correlative relationships were observed between stress wave values obtained from both field measurement techniques. Excellent relationships were also observed between standing tree and log speed-of-sound values.

Recurrent voluminous sector collapses at Volcán Barú, Panama

Two volcanic debris avalanche deposits (VDADs), both attributed to sector collapse at Volcán Barú, Panama, have been identified after an investigation of deposits that covered more than a thousand square kilometers. The younger Barriles Deposit is constrained by two radiocarbon ages that are ~9 ka; the older Caisán Deposit is at or beyond the radiocarbon range, >43,500 ybp. The total runout length of the Caisán Deposit was ~50 km and it covers 1190 km2. The Barriles Deposit extended to about 45 km and covered an area of 966 km2, overlapping most of the Caisán. The VDADs are blanketed by pyroclastic deposits and contain a predominance of andesitic material likely representing volcanic dome rock which accumulated above the active vent at Barú before collapsing. Despite heavy vegetation in the field area, over 4000 individual hummocks were digitized from aerial photography. Statistical analysis of hummock locations and geometries depict flow patterns of highly- fragmented material reflecting the effects of underlying topography and also help to define the limit of Barriles’ shorter termination. Barriles and Caisán are primarily unconfined, subaerial volcanic deposits that are among the world’s most voluminous. Calculated through two different geospatial processes, thickness values from field measurements and inferences yield volumes >30 km23 for both deposits. VDADs of comparable scale come from Mount Shasta, USA; Socompa, Chile/Argentina; and Shiveluch, Russia. Currently, the modern edifice is 200-400m lower than the pre-collapse Barriles and Caisán summits and only 16-25% of the former edifice has been replaced since the last failure.

Mechanisms contributing to the deep chlorophyll maximum in Lake Superior

The seasonal appearance of a deep chlorophyll maximum (DCM) in Lake Superior is a striking phenomenon that is widely observed; however its mechanisms of formation and maintenance are not well understood. As this phenomenon may be the reflection of an ecological driver, or a driver itself, a lack of understanding its driving forces limits the ability to accurately predict and manage changes in this ecosystem. Key mechanisms generally associated with DCM dynamics (i.e. ecological, physiological and physical phenomena) are examined individually and in concert to establish their role. First the prevailing paradigm, “the DCM is a great place to live”, is analyzed through an integration of the results of laboratory experiments and field measurements. The analysis indicates that growth at this depth is severely restricted and thus not able to explain the full magnitude of this phenomenon. Additional contributing mechanisms like photoadaptation, settling and grazing are reviewed with a one-dimensional mathematical model of chlorophyll and particulate organic carbon. Settling has the strongest impact on the formation and maintenance of the DCM, transporting biomass to the metalimnion and resulting in the accumulation of algae, i.e. a peak in the particulate organic carbon profile. Subsequently, shade adaptation becomes manifest as a chlorophyll maximum deeper in the water column where light conditions particularly favor the process. Shade adaptation mediates the magnitude, shape and vertical position of the chlorophyll peak. Growth at DCM depth shows only a marginal contribution, while grazing has an adverse effect on the extent of the DCM. The observed separation of the carbon biomass and chlorophyll maximum should caution scientists to equate the DCM with a large nutrient pool that is available to higher trophic levels. The ecological significance of the DCM should not be separated from the underlying carbon dynamics. When evaluated in its entirety, the DCM becomes the projected image of a structure that remains elusive to measure but represents the foundation of all higher trophic levels. These results also offer guidance in examine ecosystem perturbations such as climate change. For example, warming would be expected to prolong the period of thermal stratification, extending the late summer period of suboptimal (phosphorus-limited) growth and attendant transport of phytoplankton to the metalimnion. This reduction in epilimnetic algal production would decrease the supply of algae to the metalimnion, possibly reducing the supply of prey to the grazer community. This work demonstrates the value of modeling to challenge and advance our understanding of ecosystem dynamics, steps vital to reliable testing of management alternatives.

APPLICATION OF AN IMPUTATION METHOD FOR GEOSPATIAL INVENTORY OF FOREST STRUCTURAL ATTRIBUTES ACROSS MULTIPLE SPATIAL SCALES IN THE LAKE STATES, U.S.A.

Credible spatial information characterizing the structure and site quality of forests is critical to sustainable forest management and planning, especially given the increasing demands and threats to forest products and services. Forest managers and planners are required to evaluate forest conditions over a broad range of scales, contingent on operational or reporting requirements. Traditionally, forest inventory estimates are generated via a design-based approach that involves generalizing sample plot measurements to characterize an unknown population across a larger area of interest. However, field plot measurements are costly and as a consequence spatial coverage is limited. Remote sensing technologies have shown remarkable success in augmenting limited sample plot data to generate stand- and landscape-level spatial predictions of forest inventory attributes. Further enhancement of forest inventory approaches that couple field measurements with cutting edge remotely sensed and geospatial datasets are essential to sustainable forest management. We evaluated a novel Random Forest based k Nearest Neighbors (RF-kNN) imputation approach to couple remote sensing and geospatial data with field inventory collected by different sampling methods to generate forest inventory information across large spatial extents. The forest inventory data collected by the FIA program of US Forest Service was integrated with optical remote sensing and other geospatial datasets to produce biomass distribution maps for a part of the Lake States and species-specific site index maps for the entire Lake State. Targeting small-area application of the state-of-art remote sensing, LiDAR (light detection and ranging) data was integrated with the field data collected by an inexpensive method, called variable plot sampling, in the Ford Forest of Michigan Tech to derive standing volume map in a cost-effective way. The outputs of the RF-kNN imputation were compared with independent validation datasets and extant map products based on different sampling and modeling strategies. The RF-kNN modeling approach was found to be very effective, especially for large-area estimation, and produced results statistically equivalent to the field observations or the estimates derived from secondary data sources. The models are useful to resource managers for operational and strategic purposes.

Exploration of Ground Penetrating Radar and Time Domain Reflectometry Methods for the Determination of Pavement Dielectric Constant

Traditionally, densities of newly built roadways are checked by direct sampling (cores) or by nuclear density gauge measurements. For roadway engineers, density of asphalt pavement surfaces is essential to determine pavement quality. Unfortunately, field measurements of density by direct sampling or by nuclear measurement are slow processes. Therefore, I have explored the use of rapidly-deployed ground penetrating radar (GPR) as an alternative means of determining pavement quality. The dielectric constant of pavement surface may be a substructure parameter that correlates with pavement density, and can be used as a proxy when density of asphalt is not known from nuclear or destructive methods. The dielectric constant of the asphalt can be determined using ground penetrating radar (GPR). In order to use GPR for evaluation of road surface quality, the relationship between dielectric constants of asphalt and their densities must be established. Field measurements of GPR were taken at four highway sites in Houghton and Keweenaw Counties, Michigan, where density values were also obtained using nuclear methods in the field. Laboratory studies involved asphalt samples taken from the field sites and samples created in the laboratory. These were tested in various ways, including, density, thickness, and time domain reflectometry (TDR). In the field, GPR data was acquired using a 1000 MHz air-launched unit and a ground-coupled unit at 200 and 500 MHz. The equipment used was owned and operated by the Michigan Department of Transportation (MDOT) and available for this study for a total of four days during summer 2005 and spring 2006. The analysis of the reflected waveforms included “routine” processing for velocity using commercial software and direct evaluation of reflection coefficients to determine a dielectric constant. The dielectric constants computed from velocities do not agree well with those obtained from reflection coefficients. Perhaps due to the limited range of asphalt types studied, no correlation between density and dielectric constant was evident. Laboratory measurements were taken with samples removed from the field and samples created for this study. Samples from the field were studied using TDR, in order to obtain dielectric constant directly, and these correlated well with the estimates made from reflection coefficients. Samples created in the laboratory were measured using 1000 MHz air-launched GPR, and 400 MHz ground-coupled GPR, each under both wet and dry conditions. On the basis of these observations, I conclude that dielectric constant of asphalt can be reliably measured from waveform amplitude analysis of GJPR data, based on the consistent agreement with that obtained in the laboratory using TDR. Because of the uniformity of asphalts studied here, any correlation between dielectric constant and density is not yet apparent.

Hall-effect thruster-cathode coupling : the effect of cathode position and magnetic field topology

Hall-effect thruster (HET) cathodes are responsible for the generation of the free electrons necessary to initiate and sustain the main plasma discharge and to neutralize the ion beam. The position of the cathode relative to the thruster strongly affects the efficiency of thrust generation. However, the mechanisms by which the position affects the efficiency are not well understood. This dissertation explores the effect of cathode position on HET efficiency. Magnetic field topology is shown to play an important role in the coupling between the cathode plasma and the main discharge plasma. The position of the cathode within the magnetic field affects the ion beam and the plasma properties of the near-field plume, which explains the changes in efficiency of the thruster. Several experiments were conducted which explored the changes of efficiency arising from changes in cathode coupling. In each experiment, the thrust, discharge current, and cathode coupling voltage were monitored while changes in the independent variables of cathode position, cathode mass flow and magnetic field topology were made. From the telemetry data, the efficiency of the HET thrust generation was calculated. Furthermore, several ion beam and plasma properties were measured including ion energy distribution, beam current density profile, near-field plasma potential, electron temperature, and electron density. The ion beam data show how the independent variables affected the quality of ion beam and therefore the efficiency of thrust generation. The measurements of near-field plasma properties partially explain how the changes in ion beam quality arise. The results of the experiments show that cathode position, mass flow, and field topology affect several aspects of the HET operation, especially beam divergence and voltage utilization efficiencies. Furthermore, the experiments show that magnetic field topology is important in the cathode coupling process. In particular, the magnetic field separatrix plays a critical role in impeding the coupling between cathode and HET. Suggested changes to HET thruster designs are provided including ways to improve the position of the separatrix to accommodate the cathode.

Effect of external electric field on hydrogen adsorption over activated carbon separated by dielectric materials

Energy crisis and worldwide environmental problem make hydrogen a prospective energy carrier. However, storage and transportation of hydrogen in large quantities at small volume is currently not practical. Lots of materials and devices have been developed for storage hydrogen, but to today none is able to meet the DOE targets. Activated carbon has been found to be a good hydrogen adsorbent due to its high surface area. However, the weak van der Waals force between hydrogen and the adsorbent has limited the adsorption capacity. Previous studies have found that enhanced adsorption can be obtained with applied electric field. Stronger interaction between the polarized hydrogen and the charged sorbents under high voltage is considered as the reason. This study was initiated to investigate if the adsorption can be further enhanced when the activated carbon particles are separated with a dielectric coating. Dielectric TiO2 nanoparticles were first utilized. Hydrogen adsorption measurements on the TiO2-coated carbon materials, with or without an external electric field, were made. The results showed that the adsorption capacity enhancement increased with the increasing amount of TiO2 nanoparticles with an applied electric field. Since the hydrogen adsorption capacity on TiO2 particles is very low and there is no hydrogen adsorption enhancement on TiO2 particles alone when electric field is applied, the effect of dielectric coating is demonstrated. Another set of experiments investigated the behavior of hydrogen adsorption over TiO2-coated activated carbon under various electric potentials. The results revealed that the hydrogen adsorption first increased and then decreased with the increase of electric field. The improved storage was due to a stronger interaction between charged carbon surface and polarized hydrogen molecule caused by field induced polarization of TiO2 coating. When the electric field was sufficient to cause considerable ionization of hydrogen, the decrease of hydrogen adsorption occurred. The current leak detected at 3000 V was a sign of ionization of hydrogen. Experiments were also carried out to examine the hydrogen adsorption performances over activated carbon separated by other dielectric materials, MgO, ZnO and BaTiO3, respectively. For the samples partitioned with MgO and ZnO, the measurements with and without an electric field indicated negligible differences. Electric field enhanced adsorption has been observed on the activated carbon separated with BaTiO3, a material with unusually high dielectric constant. Corresponding computational calculations using Density Functional Theory have been performed on hydrogen interaction with charged TiO2 molecule as well as TiO2 molecule, coronene and TiO2-doped coronene in the presence of an electric field. The simulated results were consistent with the observations from experiments, further confirming the proposed hypotheses.

Laser Thomson scattering measurements of electron temperature and density in a Hall-effect plasma

Hall-effect thrusters (HETs) are compact electric propulsion devices with high specific impulse used for a variety of space propulsion applications. HET technology is well developed but the electron properties in the discharge are not completely understood, mainly due to the difficulty involved in performing accurate measurements in the discharge. Measurements of electron temperature and density have been performed using electrostatic probes, but presence of the probes can significantly disrupt thruster operation, and thus alter the electron temperature and density. While fast-probe studies have expanded understanding of HET discharges, a non-invasive method of measuring the electron temperature and density in the plasma is highly desirable. An alternative to electrostatic probes is a non-perturbing laser diagnostic technique that measures Thomson scattering from the plasma. Thomson scattering is the process by which photons are elastically scattered from the free electrons in a plasma. Since the electrons have thermal energy their motion causes a Doppler shift in the scattered photons that is proportional to their velocity. Like electrostatic probes, laser Thomson scattering (LTS) can be used to determine the temperature and density of free electrons in the plasma. Since Thomson scattering measures the electron velocity distribution function directly no assumptions of the plasma conditions are required, allowing accurate measurements in anisotropic and non-Maxwellian plasmas. LTS requires a complicated measurement apparatus, but has the potential to provide accurate, non-perturbing measurements of electron temperature and density in HET discharges. In order to assess the feasibility of LTS diagnostics on HETs non-invasive measurements of electron temperature and density in the near-field plume of a Hall thruster were performed using a custom built laser Thomson scattering diagnostic. Laser measurements were processed using a maximum likelihood estimation method and results were compared to conventional electrostatic double probe measurements performed at the same thruster conditions. Electron temperature was found to range from approximately 1 – 40 eV and density ranged from approximately 1.0 x 1017 m-3 to 1.3 x 1018 m-3 over discharge voltages from 250 to 450 V and mass flow rates of 40 to 80 SCCM using xenon propellant.