Analyzing the effects of vertical acceleration and deceleration in multi-phase, debris-flow sediment-column experiments / by Christine Marie Williams.

Water-saturated debris flows are among some of the most destructive mass movements. Their complex nature presents a challenge for quantitative description and modeling. In order to improve understanding of the dynamics of these flows, it is important to seek a simplified dynamic system underlying their behavior. Models currently in use to describe the motion of debris flows employ depth-averaged equations of motion, typically assuming negligible effects from vertical acceleration. However, in many cases debris flows experience significant vertical acceleration as they move across irregular surfaces, and it has been proposed that friction associated with vertical forces and liquefaction merit inclusion in any comprehensive mechanical model. The intent of this work is to determine the effect of vertical acceleration through a series of laboratory experiments designed to simulate debris flows, testing a recent model for debris flows experimentally. In the experiments, a mass of water-saturated sediment is released suddenly from a holding container, and parameters including rate of collapse, pore-fluid pressure, and bed load are monitored. Experiments are simplified to axial geometry so that variables act solely in the vertical dimension. Steady state equations to infer motion of the moving sediment mass are not sufficient to model accurately the independent solid and fluid constituents in these experiments. The model developed in this work more accurately predicts the bed-normal stress of a saturated sediment mass in motion and illustrates the importance of acceleration and deceleration.

Vegetation and hydrologic influences on carbon and nitrogen in subsurface water of a forested riparian wetland

Vegetation communities affect carbon and nitrogen dynamics in the subsurface water of mineral wetlands through the quality of their litter, their uptake of nutrients, root exudation and their effects on redox potential. However, vegetation influence on subsurface nutrient dynamics is often overshadowed by the influences of hydrology, soils and geology on nutrient dynamics. The effects of vegetation communities on carbon and nitrogen dynamics are important to consider when managing land that may change vegetation type or quantity so that wetland ecosystem functions can be retained. This study was established to determine the magnitude of the influences and interaction of vegetation cover and hydrology, in the form of water table fluctuations, on carbon and nitrogen dynamics in a northern forested riparian wetland. Dissolved organic carbon (DOC), dissolved inorganic carbon (DIC), nitrate (NO3-) and ammonium (NH4+) concentrations were collected from a piezometer network in four different vegetation communities and were found to show complex responses to vegetation cover and water table fluctuations. Dissolved organic carbon, DIC, NO3- and NH4+ concentrations were influenced by forest vegetation cover. Both NO3- and NH4+ were also influenced by water table fluctuations. However, for DOC and NH4+ concentrations there appeared to be more complex interactions than were measured by this study. The results of canonical correspondence analysis (CCA) and analysis of variance (ANOVA) did not correspond in relationship to the significance of vegetation communities. Dissolved inorganic carbon was influenced by an interaction between vegetation cover and water table fluctuations. More hydrological information is needed to make stronger conclusions about the relationship between vegetation and hydrology in controlling carbon and nitrogen dynamics in a forested riparian wetland.

Carbon dioxide sequestration in cement kiln dust through miner carbonation

The feasibility of carbon sequestration in cement kiln dust (CKD) was investigated in a series of batch and column experiments conducted under ambient temperature and pressure conditions. The significance of this work is the demonstration that alkaline wastes, such as CKD, are highly reactive with carbon dioxide (CO2). In the presence of water, CKD can sequester greater than 80% of its theoretical capacity for carbon without any amendments or modifications to the waste. Other mineral carbonation technologies for carbon sequestration rely on the use of mined mineral feedstocks as the source of oxides. The mining, pre-processing and reaction conditions needed to create favorable carbonation kinetics all require significant additions of energy to the system. Therefore, their actual net reduction in CO2 is uncertain. Many suitable alkaline wastes are produced at sites that also generate significant quantities of CO2. While independently, the reduction in CO2 emissions from mineral carbonation in CKD is small (~13% of process related emissions), when this technology is applied to similar wastes of other industries, the collective net reduction in emissions may be significant. The technical investigations presented in this dissertation progress from proof of feasibility through examination of the extent of sequestration in core samples taken from an aged CKD waste pile, to more fundamental batch and microscopy studies which analyze the rates and mechanisms controlling mineral carbonation reactions in a variety of fresh CKD types. Finally, the scale of the system was increased to assess the sequestration efficiency under more pilot or field-scale conditions and to clarify the importance of particle-scale processes under more dynamic (flowing gas) conditions. A comprehensive set of material characterization methods, including thermal analysis, Xray diffraction, and X-ray fluorescence, were used to confirm extents of carbonation and to better elucidate those compositional factors controlling the reactions. The results of these studies show that the rate of carbonation in CKD is controlled by the extent of carbonation. With increased degrees of conversion, particle-scale processes such as intraparticle diffusion and CaCO3 micropore precipitation patterns begin to limit the rate and possibly the extent of the reactions. Rates may also be influenced by the nature of the oxides participating in the reaction, slowing when the free or unbound oxides are consumed and reaction conditions shift towards the consumption of less reactive Ca species. While microscale processes and composition affects appear to be important at later times, the overall degrees of carbonation observed in the wastes were significant (> 80%), a majority of which occurs within the first 2 days of reaction. Under the operational conditions applied in this study, the degree of carbonation in CKD achieved in column-scale systems was comparable to those observed under ideal batch conditions. In addition, the similarity in sequestration performance among several different CKD waste types indicates that, aside from available oxide content, no compositional factors significantly hinder the ability of the waste to sequester CO2.

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.

Simulating dioxane transport in a heterogeneous glacial aquifer system (Washtenaw County, Michigan) using publicly available models and data

The primary challenge in groundwater and contaminant transport modeling is obtaining the data needed for constructing, calibrating and testing the models. Large amounts of data are necessary for describing the hydrostratigraphy in areas with complex geology. Increasingly states are making spatial data available that can be used for input to groundwater flow models. The appropriateness of this data for large-scale flow systems has not been tested. This study focuses on modeling a plume of 1,4-dioxane in a heterogeneous aquifer system in Scio Township, Washtenaw County, Michigan. The analysis consisted of: (1) characterization of hydrogeology of the area and construction of a conceptual model based on publicly available spatial data, (2) development and calibration of a regional flow model for the site, (3) conversion of the regional model to a more highly resolved local model, (4) simulation of the dioxane plume, and (5) evaluation of the model's ability to simulate field data and estimation of the possible dioxane sources and subsequent migration until maximum concentrations are at or below the Michigan Department of Environmental Quality's residential cleanup standard for groundwater (85 ppb). MODFLOW-2000 and MT3D programs were utilized to simulate the groundwater flow and the development and movement of the 1, 4-dioxane plume, respectively. MODFLOW simulates transient groundwater flow in a quasi-3-dimensional sense, subject to a variety of boundary conditions that can simulate recharge, pumping, and surface-/groundwater interactions. MT3D simulates solute advection with groundwater flow (using the flow solution from MODFLOW), dispersion, source/sink mixing, and chemical reaction of contaminants. This modeling approach was successful at simulating the groundwater flows by calibrating recharge and hydraulic conductivities. The plume transport was adequately simulated using literature dispersivity and sorption coefficients, although the plume geometries were not well constrained.

Sonic log prediction in carbonates

This work is conducted to study the complications associated with the sonic log prediction in carbonate logs and to investigate the possible solutions to accurately predict the sonic logs in Traverse Limestone. Well logs from fifty different wells were analyzed to define the mineralogy of the Traverse Limestone by using conventional 4-mineral and 3-mineral identification approaches. We modified the conventional 3-mineral identification approach (that completely neglects the gamma ray response) to correct the shale effects on the basis of gamma ray log before employing the 3-mineral identification. This modification helped to get the meaningful insight of the data when a plot was made between DGA (dry grain density) and UMA (Photoelectric Volumetric Cross-section) with the characteristic ternary diagram of the quartz, calcite and dolomite. The results were then compared with the 4-mineral identification approach. Contour maps of the average mineral fractions present in the Traverse Limestone were prepared to see the basin wide mineralogy of Traverse Limestone. In the second part, sonic response of Traverse Limestone was predicted in fifty randomly distributed wells. We used the modified time average equation that accounts for the shale effects on the basis of gamma ray log, and used it to predict the sonic behavior from density porosity and average porosity. To account for the secondary porosity of dolomite, we subtracted the dolomitic fraction of clean porosity from the total porosity. The pseudo-sonic logs were then compared with the measured sonic logs on the root mean square (RMS) basis. Addition of dolomite correction in modified time average equation improved the results of sonic prediction from neutron porosity and average porosity. The results demonstrated that sonic logs could be predicted in carbonate rocks with a root mean square error of about 4μsec/ft. We also attempted the use of individual mineral components for sonic log prediction but the ambiguities in mineral fractions and in the sonic properties of the minerals limited the accuracy of the results.

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.

Paleomagnetic Study of the Portage Lake Volcanics exposed in the Quincy Mine

A detailed paleomagnetic and rock-magnetic investigation was conducted on thirty six basaltic flows of the ~1095 Ma Portage Lake Volcanics. The flows were sampled along the East Adit of the Quincy Mine (Hancock, MI). Thirty two flows yielded well-defined primary magnetization directions carried by magnetite. A secondary magnetization component carried by hematite was also found in twenty nine flows. After correction for serial correlation between the flows, nineteen independent mean directions were calculated. The corresponding paleomagnetic pole is located at 25.5 °N, 182.1 °W (A95 = 3.5°). The new pole overlaps with the pole from the ~1087 Ma Lake Shore Traps suggesting a standstill of the North American plate during that time period. The low angular dispersion of virtual geomagnetic poles (S = 7.9°) suggests that the flows were erupted within a short time period, or that the strength of geomagnetic secular variation was lower than that of the recent field.

Investigation of surface water/groundwater relationships in a rural watershed in Southern Honduras

Despite failed attempts at obtaining a potable water system, the village of El Caracol in Southern Honduras remains committed to improving access to water resources. To assist in this endeavor, an investigation of the hydrogeological characteristics of the local watershed was conducted. Daily precipitation was recorded to examine the relationship between precipitation and approximated river and spring discharges. A Thornthwaite Mather Water Balance Model was used to predict monthly discharges for comparison with observed values, and to infer the percentage of topographic watersheds contributing to the respective discharges. As aquifer porosity in this region is thought to be primarily secondary (i.e., fractures), field observed lineaments were compared with those interpreted from remote sensing imagery in an attempt to determine the usefulness of these interpretations in locating potential water sources for a future project.

Fossil bubble in porphyritic basaltic pyroclasts produced by small and large strombolian eruption at Pacaya, Guatemala

We investigate compositionally monotonous, but energetically diverse, tephra samples from Pacaya to see if fossil bubbles in pyroclasts could reflect eruptive style. Bubble size distributions (BSD) were determined for four ash to lapilli size tephra samples using an adapted version of stereology conversion by Sahagian and Proussevitch (1998). Eruptions range from very weak to very energetic. Hundreds of ESEM BSEs images were processed throughout ImageJ software for a robust and statistically correct data set of vesicles (minimum 700 bubbles per sample). Qualitative textural analysis and major element chemical compositions were also executed. There is higher vesicularity for explosive pyroclasts and an inverse correlation between bubble number density (NV) and explosivity.

TECHNOLOGICAL ADVANCEMENT IN REMOTE SENSING AND INNOVATIVE APPLICATION TO VOLCANOLOGY

The integration of remote monitoring techniques at different scales is of crucial importance for monitoring of volcanoes and assessment of the associated hazard. In this optic, technological advancement and collaboration between research groups also play a key role. Vhub is a community cyberinfrastructure platform designed for collaboration in volcanology research. Within the Vhub framework, this dissertation focuses on two research themes, both representing novel applications of remotely sensed data in volcanology: advancement in the acquisition of topographic data via active techniques and application of passive multi-spectral satellite data to monitoring of vegetated volcanoes. Measuring surface deformation is a critical issue in analogue modelling of Earth science phenomena. I present a novel application of the Microsoft Kinect sensor to measurement of vertical and horizontal displacements in analogue models. Specifically, I quantified vertical displacement in a scaled analogue model of Nisyros volcano, Greece, simulating magmatic deflation and inflation and related surface deformation, and included the horizontal component to reconstruct 3D models of pit crater formation. The detection of active faults around volcanoes is of importance for seismic and volcanic hazard assessment, but not a simple task to be achieved using analogue models. I present new evidence of neotectonic deformation along a north-south trending fault from the Mt Shasta debris avalanche deposit (DAD), northern California. The fault was identified on an airborne LiDAR campaign of part of the region interested by the DAD and then confirmed in the field. High resolution LiDAR can be utilized also for geomorphological assessment of DADs, and I describe a size-distance analysis to document geomorphological aspects of hummock in the Shasta DAD. Relating the remote observations of volcanic passive degassing to conditions and impacts on the ground provides an increased understanding of volcanic degassing and how satellite-based monitoring can be used to inform hazard management strategies in nearreal time. Combining a variety of satellite-based spectral time series I aim to perform the first space-based assessment of the impacts of sulfur dioxide emissions from Turrialba volcano, Costa Rica, on vegetation in the surrounding environment, and establish whether vegetation indices could be used more broadly to detect volcanic unrest.

PROBABILISTIC MODELING OF RAINFALL INDUCED LANDSLIDE HAZARD ASSESSMENT IN SAN JUAN LA LAGUNA, SOLOLÁ, GUATEMALA

The municipality of San Juan La Laguna, Guatemala is home to approximately 5,200 people and located on the western side of the Lake Atitlán caldera. Steep slopes surround all but the eastern side of San Juan. The Lake Atitlán watershed is susceptible to many natural hazards, but most predictable are the landslides that can occur annually with each rainy season, especially during high-intensity events. Hurricane Stan hit Guatemala in October 2005; the resulting flooding and landslides devastated the Atitlán region. Locations of landslide and non-landslide points were obtained from field observations and orthophotos taken following Hurricane Stan. This study used data from multiple attributes, at every landslide and non-landslide point, and applied different multivariate analyses to optimize a model for landslides prediction during high-intensity precipitation events like Hurricane Stan. The attributes considered in this study are: geology, geomorphology, distance to faults and streams, land use, slope, aspect, curvature, plan curvature, profile curvature and topographic wetness index. The attributes were pre-evaluated for their ability to predict landslides using four different attribute evaluators, all available in the open source data mining software Weka: filtered subset, information gain, gain ratio and chi-squared. Three multivariate algorithms (decision tree J48, logistic regression and BayesNet) were optimized for landslide prediction using different attributes. The following statistical parameters were used to evaluate model accuracy: precision, recall, F measure and area under the receiver operating characteristic (ROC) curve. The algorithm BayesNet yielded the most accurate model and was used to build a probability map of landslide initiation points. The probability map developed in this study was also compared to the results of a bivariate landslide susceptibility analysis conducted for the watershed, encompassing Lake Atitlán and San Juan. Landslides from Tropical Storm Agatha 2010 were used to independently validate this study’s multivariate model and the bivariate model. The ultimate aim of this study is to share the methodology and results with municipal contacts from the author's time as a U.S. Peace Corps volunteer, to facilitate more effective future landslide hazard planning and mitigation.

A STUDY OF THE IMPACTS OF FREEZE-THAW ON CLIFF RECESSION AT THE CALVERT CLIFFS IN CALVERT COUNTY, MARYLAND

The Calvert Cliffs, which form much of the western coastline of the Chesapeake Bay in Calvert County, Maryland, are actively eroding and destabilizing, resulting in a critical situation for many homes in close proximity to the slope's crest. Past studies have identified that where waves directly interact with the toe of the slope, wave action controls cliff recession; however, where waves do not regularly interact with the slope toe, the past work identified that freeze-thaw controls recession. This study investigated the validity of this second claim by analyzing the recession rate and freeze-thaw behavior of six study sites along the Calvert Cliffs that are not directly affected by waves. While waves do remove failed material from the toe, in these regions freeze-thaw is believed to be the dominant factor driving recession at the Calvert Cliffs. Past recession rates were calculated using historical aerial photographs and were analyzed together with a number of other variables selected to represent the freeze-thaw behavior of the Calvert Cliffs. The investigation studied sixteen independent variables and found that over 65% of recession at these study sites can be represented by the following five variables: (1) cliff face direction, (2 and 3) the percent of total cliff height composed of soil with freeze-thaw susceptibility F4 and F2, (4) the number of freeze-thaw cycles, and (5) the weighted shear strength. Future mitigation techniques at these sites should focus on addressing these variables and might include vegetation or addressing the presence of water along the face of the slope. Unmitigated, the Calvert Cliffs will continue to recede until a stable slope angle is reached and maintained.

INVESTIGATION OF STRATEGIES TO PROMOTE EFFECTIVE TEACHER PROFESSIONAL DEVELOPMENT EXPERIENCES IN EARTH SCIENCE

This dissertation serves as a call to geoscientists to share responsibility with K-12 educators for increasing Earth science literacy. When partnerships are created among K-12 educators and geoscientists, the synergy created can promote Earth science literacy in students, teachers, and the broader community. The research described here resulted in development of tools that can support effective professional development for teachers. One tool is used during the planning stages to structure a professional development program, another set of tools supports measurement of the effectiveness of a development program, and the third tool supports sustainability of professional development programs. The Michigan Teacher Excellence Program (MiTEP), a Math/Science Partnership project funded by the National Science Foundation, served as the test bed for developing and testing these tools. The first tool, the planning tool, is the Earth Science Literacy Principles (ESLP). The ESLP served as a planning tool for the two-week summer field courses as part of the MiTEP program. The ESLP, published in 2009, clearly describe what an Earth science literate person should know. The ESLP consists of nine big ideas and their supporting fundamental concepts. Using the ESLP for planning a professional development program assisted both instructors and teacher-participants focus on important concepts throughout the professional development activity. The measurement tools were developed to measure change in teachers’ Earth science content-area knowledge and perceptions related to teaching and learning that result from participating in a professional development program. The first measurement tool, the Earth System Concept Inventory (ESCI), directly measures content-area knowledge through a succession of multiple-choice questions that are aligned with the content of the professional development experience. The second measurement, an exit survey, collects qualitative data from teachers regarding their impression of the professional development. Both the ESCI and the exit survey were tested for validity and reliability. Lesson study is discussed here as a strategy for sustaining professional development in a school or a district after the end of a professional development activity. Lesson study, as described here, was offered as a formal course. Teachers engaged in lesson study worked collaboratively to design and test lessons that improve the teachers’ classroom practices. Data regarding the impact of the lesson study activity were acquired through surveys, written documents, and group interviews. The data are interpreted to indicate that the lesson study process improved teacher quality and classroom practices. In the case described here, the lesson study process was adopted by the teachers’ district and currently serves as part of the district’s work in Professional Learning Communities, resulting in ongoing professional development throughout the district.