The proximity of Miocene Columbia River basalts to the "locally-erupted" coastal Miocene basalts in northwestern Oregon, and the compelling similarities between the two groups, suggest that the coastal basalts, rather than being locally erupted, may be the westward extension of plateau basalts derived from eastern Oregon and Washington. The local-origin hypothesis is based largely on the interpretation of coastal dikes and sills as representing vent areas; however, a complex mechanism, as yet unsatisfactorily defined, would be required to cause the eruption of virtually identical magmas simultaneously from source areas 500 km apart. This study, therefore, has investigated the coastal basalt intrusions both laterally and vertically. Geochemical and paleomagnetic analysis was used to determine the occurrence and distribution of basalt units; gravity surveys enabled an examination of the subsurface extensions of basalt intrusions in sedimentary rocks.

Uranium contamination of the subsurface remains a significant problem at the Department of Energy Hanford site. A series of column experiments were conducted on Hanford sediment saturated with simulated groundwater to study the effects of aqueous bicarbonate and microbial growth on the mobility of Uranium. Spectral induced polarization (SIP) measurements in the columns were conducted concurrently with pore water sampling in order to monitor changes occurring inside the sediment after the initiation of microbial growth induced by glucose injection. The microbial growth caused significant increases in the real component of the complex conductivity and is the result of ion release into the pore fluid. In addition, an increase in the imaginary conductivity was observed at low frequencies (Hz), which may be due to biotic processes. Due to the use of natural sediment, the SIP response is complex and difficult to understand. However, results across all columns with microbial growth are consistent. Pore water testing showed that microbial growth leads to sudden increases in uranium concentrations; however, microbes also eventually create reducing conditions in the sediment which transforms soluble U6+ to insoluble U4+. Bicarbonate leads to significant increases in uranium concentrations likely due to the formation of mobile uranyl carbonate complexes. For the purposes of field scale remediation, microbial growth in an oxic environment should be avoided. However, within reducing conditions present in the deep vadose zone and phreatic zone, microbial growth seems unlikely to significantly increase uranium mobility.

An interdisciplinary approach at the interface between geodesy and geophysics has recently resolved several Earth science problems at regional and global scales. I use the term “geophysical geodesy” to distinguish the technical and theoretical aspect of geodesy from geophysical applications of geodetic techniques. Using a wide range of Earth observation data, I study the spatio-temporal characteristics of Earth surface deformation in the United States associated with several geophysical processes, including natural and anthropogenic subsidence and uplift, regional relative sea-level rise, and continental hydrological loading. The theoretical portion of this dissertation applies loading theory and develops a new hybrid method to improve the estimate of hydrologically-induced vertical deformation at time scales from sub-annual to multi-annual. The application part of this dissertation benefits from GPS and other geodetic and geologic data sets to study and model Earth’s surface uplift due to CO2 injection at an oil reservoir in coastal Texas, and coastal subsidence and nuisance flooding along the Mississippi River Delta and eastern seaboard of the United States.

The subject of this dissertation is use of computer models as data analysis tools in several different geoscience settings, including integrated surface water/groundwater modeling, tephra fallout modeling, geophysical inversion, and hydrothermal groundwater modeling. The dissertation is organized into three chapters, which correspond to three individual publication manuscripts. In the first chapter, a linear framework is developed to identify and estimate the potential predictive consequences of using a simple computer model as a data analysis tool. The framework is applied to a complex integrated surface-water/groundwater numerical model with thousands of parameters. Several types of predictions are evaluated, including particle travel time and surface-water/groundwater exchange volume. The analysis suggests that model simplifications have the potential to corrupt many types of predictions. The implementation of the inversion, including how the objective function is formulated, what minimum of the objective function value is acceptable, and how expert knowledge is enforced on parameters, can greatly influence the manifestation of model simplification. Depending on the prediction, failure to specifically address each of these important issues during inversion is shown to degrade the reliability of some predictions. In some instances, inversion is shown to increase, rather than decrease, the uncertainty of a prediction, which defeats the purpose of using a model as a data analysis tool. In the second chapter, an efficient inversion and uncertainty quantification approach is applied to a computer model of volcanic tephra transport and deposition. The computer model simulates many physical processes related to tephra transport and fallout. The utility of the approach is demonstrated for two eruption events. In both cases, the importance of uncertainty quantification is highlighted by exposing the variability in the conditioning provided by the observations used for inversion. The worth of different types of tephra data to reduce parameter uncertainty is evaluated, as is the importance of different observation error models. The analyses reveal the importance using tephra granulometry data for inversion, which results in reduced uncertainty for most eruption parameters. In the third chapter, geophysical inversion is combined with hydrothermal modeling to evaluate the enthalpy of an undeveloped geothermal resource in a pull-apart basin located in southeastern Armenia. A high-dimensional gravity inversion is used to define the depth to the contact between the lower-density valley fill sediments and the higher-density surrounding host rock. The inverted basin depth distribution was used to define the hydrostratigraphy for the coupled groundwater-flow and heat-transport model that simulates the circulation of hydrothermal fluids in the system. Evaluation of several different geothermal system configurations indicates that the most likely system configuration is a low-enthalpy, liquid-dominated geothermal system.

The role of water and volatiles in the solar system is of critical interest in planetary science. Evidence for the past action of water or direct observation of water on a planetary body can indicate the potential to harbor life and is critical to human exploration of the solar system. We study two very different remote sensing techniques that address the issue of identifying water-related processes on the surface of other planetary bodies, and in particular, Mars. The first technique, combined thermal infrared and visible imaging, has been used extensively on Mars for determining the thermal inertia of surface materials. In the second part of this dissertation, we develop a technique that combines remote thermophysical and visible data sets with ground-based field investigations for the identification of sedimentary features at the surfaces of alluvial fans. Several methods for remotely identifying sedimentary features will be explored using thermal and visible images. We combine results from pre-existing ground-based studies with thermal images and ground-based field investigations to develop a robust technique to be used on a variety of alluvial fans. In the third part, we characterize the remote thermophysical and visible properties of specific classes of sedimentary features on alluvial fans using the technique developed in part two. The second remote sensing technique, neutron spectroscopy, has been used on many planetary spacecraft missions for the identification of hydrogen on planetary surfaces. The Dynamic Albedo of Neutrons (DAN) instrument on the upcoming Mars Science Laboratory rover mission represents a new type of neutron detector for planetary spacecraft, with the neutron detectors mounted to a rover on the Martian surface (as opposed to in orbit around the planetary body) and neutron counts that are binned by time, energy, and location (as opposed to just by energy and location). In chapter four, we model expected neutron energies and arrival times for geologic settings where water has altered the chemistry of the near surface using available geochemical data from the Mars Exploration Rovers (MER). Particle transport models are used to determine the sensitivity of neutron detection techniques to the variations in hydrogen abundance, hydrogen layering and chemical composition measured by MER.

Our understanding of tephra depositional processes is significantly improved by high-resolution ground-penetrating radar (GPR) data collected at Cerro Negro volcano, Nicaragua. The data reveal three depositional regimes: (1) a near-vent region on the cone itself, where 10 GPR radargrams collected on the western flank show quantifiable differences between facies formed from low energy normal Strombolian and higher energy violent Strombolian processes, indicating imaging of scoria cone deposits may be useful in distinguishing eruptive style in older cones where the proximal to distal tephra blanket has eroded away; (2) a proximal zone in which horizons identified in crosswind profiles collected at distances of 700 and 1,000 m from the vent exhibit Gaussian distributions with a high degree of statistical confidence, with tephra thickness decreasing exponentially downwind from the cone base (350 m) to ~ 1,200 m from the vent, and where particles fall from a height of less than ~2 km; and (3) a medial zone, in which particles fall from ~4 to 7 km and the deposit is thicker than expected based on thinning trends observed in the proximal zone of the deposit, indicating a transition from sedimentation dominated by fallout from plume margins to that dominated by fallout from the buoyant eruption cloud. Horizons identified in a crosswind profile at 1600 m from vent exhibit Gaussian distributions, again with high degrees of statistical confidence. True diffusion coefficients are calculated from Gaussian fits of crosswind profiles and do not show any statistical variation between zones (2) and (3). Data display thinning trends that agree with the morphology predicted by the advection-diffusion equation to a high degree of statistical confidence, validating the use of this class of models in tephra forecasting. One such model, the Tephra2 model, is reformulated for student use. A strategy is presented for utilizing this research-caliber model to introduce university undergraduates to key concepts in model literacy, encouraging students to develop a deeper understanding of the applicability and limitations of hazard models generally. For this purpose, the Tephra2 numerical model is implemented on the VHub.org website, a venture in cyberinfrastructure that brings together volcanological models and educational materials, and provides students with the ability to explore and execute sophisticated numerical models like Tephra2.

Abstract Sinkholes and sinkhole-related features in West-Central Florida (WCF) are commonly identified using geotechnical investigations such as standard penetration test (SPT) borings and geophysical methods such as ground penetrating radar (GPR) and electrical resistivity tomography (ERT). Geophysical investigation results can be used to locate drilling and field testing sites while geotechnical investigation can be used to ground truth geophysical results. Both methods can yield complementary information. Geotechnical investigations give important information about the type of soil, groundwater level and presence of low-density soils or voids at the test location, while geophysical investigations like GPR surveys have better spatial coverage and can resolve shallow stratigraphic indicators of subsidence. In GPR profiles collected at 103 residential sites in covered-karst terrain in WCF, sinkhole-related anomalies are identified using GPR and SPT methods. We analyze the degree to which the shallow features imaged in GPR correlate spatially with the N-values (blow counts) derived from SPTs at the 103 residential sites. GPR anomalies indicating sinkhole activity are defined as zones where subsurface layers show local downwarping, discontinuities, or sudden increases in amplitude or penetration of the GPR signal. "Low SPT values" indicating sinkhole activity are defined using an optimization code that searched for threshold SPT value showing optimum correlation between GPR and SPT for different optimal depth ranges. We also compared these criteria with other commonly used geotechnical criteria such as weight of rod and weight of hammer conditions. Geotechnical results were also used to filter the data based on site characteristics such as presence of shallow clay layers to study the effectiveness of GPR at different zones. Subsets of the dataset are further analyzed based on geotechnical results such as clay thickness, bedrock depth, groundwater conditions and other geological factors such as geomorphology, lithology, engineering soil type, soil thickness and prevalent sinkhole type. Results are used to examine (1) which SPT indicators show the strongest correlations with GPR anomalies, (2) the degree to which GPR surveys improve the placement of SPT borings, and (3) what these results indicate about the structure of sinkholes at these sites. For the entire data set, we find a statistically significant correlation between GPR anomalies and low SPT N-values with a confidence level of 90%. Logistic regression analysis shows that the strongest correlations are between GPR anomalies and SPT values measured in the depth range of 0-4.5 m. The probability of observing a GPR anomaly on a site will decrease by up to 84% as the minimum SPT value increases from 0 to 20 in the general study area. Boreholes drilled on GPR anomalies are statistically significantly more likely to show zones of anomalously low SPT values than boreholes drilled off GPR anomalies. We also find that the optimum SPT criteria result in better correlation with GPR than other simple commonly used geotechnical criteria such as weight of rod and weight of hammer. Better correlations were found when sites with poor GPR penetrations are filtered out from the dataset. The odds ratio showed similar result while the result varied with the depth range, statistics and threshold SPT value (low N- value with optimum correlation), with a maximum observed odds ratio of 3. Several statistical results suggest that raveling zones that connect voids to the surface may be inclined, so that shallow GPR anomalies are laterally offset from deeper zones of low N-values. Compared to the general study area, we found locally stronger correlation in some sub-regions. For example, the odds ratio found for tertiary hawthorn subgroup were 25 times higher than the odds ratio found for the general study area (WCF).

Concepción is the most active composite volcano in Nicaragua, and is located on Ometepe Island, within Lake Nicaragua. Moderate to small volcanic explosions with a volcanic explosivity index (VEI) of 1-2 have been characteristic of this volcano during the last four decades. Although its current activity is not violent, its volcanic deposits reveal stages of violent activity involving Plinian and sub-Plinian eruptions that deposited vast amounts of volcanic tephra in the Atlantic Ocean. These observations, together with the 31,000 people living on the island, make Concepción volcano an important target for volcanological research. My research focuses on the investigation of the stability of the volcano edifice of Concepción, using geophysical data such as gravity, geodetic global positioning system (GPS), sulphur dioxide (SO2) flux, real-time seismic amplitude (RSAM), and satellite remotely-sensed data. The integration of these data sets provides information about the short-term behavior of Concepción, and some insights into the volcano's long-term behavior. This study has provided, for the first time, information about the shallow dynamics of Concepción on time scales of days to weeks. I furnish evidence that this volcano is not gravitationally spreading in a continuous fashion as previously thought, that its bulk average density is comparable to that of a pile of gravel, that the volcano edifice is composed of two major distinctive lithologies, that the deformation field around the volcano is recoverable in a matter of days, and that the deformation source is located in the shallow crust. This source is also degassing through the relatively open magmatic conduit. There are, however, several remaining questions. Although the volcano is not spreading continuously there is the possibility that gravitational spreading may be taking place in a stick-slip fashion. This has important implications for slope stability of the volcano, and the associated hazards. The factors influencing the long term slope stability of the volcano are still not fully resolved, but internal volcanic processes and anthropogenic disturbances appear to be the major factors.

Telica Volcano, Nicaragua, is a persistently restless volcano with high rates of seismicity that can vary from less than ten events to over a thousand events per day. Low-frequency (LF) events dominate the seismic catalogue and seismicity rates at Telica show little clear correlation with periods of eruption. As such, traditional methods of forecasting of volcanic activity based on increases in seismicity and recognition of LF activity are not applicable. A single seismic station has been operating at Telica since 1993, and in 2010 we installed a broadband seismic and continuous GPS network (TESAND network) at Telica. In this study we investigate the seismic characteristics surrounding a nine-month period of phreatic to phreatomagmatic explosions in 1999, and also from the initial three-and-a-half year deployment of the TESAND network, including a three-month phreatic vulcanian eruptive period in 2011. We demonstrate that pertinent information can be obtained from analysis of single-station data, and while large seismic networks are preferable when possible, we note that for many volcanoes this is not possible. We find unusual patterns of seismicity before both eruptive periods; rather than a precursory increase in seismicity as is observed prior to many volcanic eruptions, we observe a decrease in seismicity many months prior to eruption. We developed a new program for cross-correlation of large seismic data catalogues and analysed multiplet activity surrounding both eruptive periods. We observed that the formation of new multiplets corresponds to periods of high event rates (during inter-eruptive periods) and high percentages of daily events that belong to a multiplet. We propose a model for the seismicity patterns observed at Telica, where changes in seismicity are related to a cyclic transition between open-system degassing and closed-system degassing. Periods of open-system degassing occur during non-eruptive episodes and are characterised by high event rates, a broad range of frequency content of events and high degrees of waveform correlation. A transition to closed-system degassing could be due to sealing of fluid pathways in the magmatic and/or hydrothermal system, or due to magma withdrawal. Periods of closed-system degassing are characterised by low event rates, higher frequency contents and low degrees of waveform correlation. Eruptive periods may then represent a transition from closed-system degassing to open-system degassing, however the system must also be capable of transitioning to open-system degassing without eruption. These observations have important implications for volcano monitoring and eruption forecasting at persistently restless volcanoes. Rather than a precursory increase in seismicity as is often observed prior to eruption at other volcanoes, our observations indicate that phreatic eruptions at Telica occur after a decrease in seismicity, a corresponding change in the frequency content of events, and a decrease in waveform correlation. These changes may represent a period of closed-system degassing that could culminate in phreatic eruptions. The inclusion of real-time analysis of variations in frequency content and multiplet activity provides critical information for volcano monitoring institutions.

This is a new method of determining magnetization strength, direction, and paleopole location from magnetic anomalies across edges of the equatorial band of magnetic sources on Mars. Different assumed locations of paleopoles result in different inducing field directions in the vicinity of an edge. Thus, with different paleopoles, the resulting magnetic fields from the edges of magnetic sources are different, and correlate differently with the observed fields. Best correlating observed and computed magnetic edge effect fields yield the potential paleopoles. The total gradient (TG) of the z-component magnetic field was used to identify the edges of magnetization boundaries. Three edge segments yielded meaningful paleopoles. They are: (15°S, 285°E)/(0°, 292.5°E) across the northeast trending edge east of Tharsis; (15°N, 195°E) across the southern boundary located in the Terra Cimmeria and Terra Sirenum area; and (15°S, 165°E)/(45°S, 180°E) across the northeast edge of Hellas Planitia. The effective magnetization in these regions (assuming 40km magnetic layer thickness) is 7.02, 42.13-94.79, and 2.63-3.51 A/m, respectively. Evidence from the overlap of regions of TG and chronostratigraphy suggests that the dynamo was active during the Noachian and may have been active in the early Hesperian.