Para a realização dessa pesquisa foi desenvolvido um estudo da espessura crustal nos Andes setentrionais e na parte norte dos Andes centrais (entre 10ºN e 9ºS). Foram usadas reflexões da Moho de seis sismos de profundidade média para achar a espessura crustal na área de estudo. As ondas pmP e smS (reflexões na Moho) são encontradas como precursoras das fases profundas pP e sS de sismos registrados em estações localizadas a distâncias telessísmicas. A metodologia utilizada exigia sismos com profundidade maiores que 100 km e magnitude maiores que 6. Entretanto, devido à baixa amplitude apresentada pelas ondas pmP e smS, foi necessário realizar um empilhamento de vários sismogramas para realçar as chegadas das ondas analisadas. A área estudada foi dividida em três partes: área norte entre 6ºN e 0º(na Colômbia), área central entre 2ºS e 5ºS (no Equador) e área sul entre 6ºS e 9ºS (na borda entre Peru e Brasil). Na área norte foram encontradas espessuras crustais entre 26 e 56 km em quatro pontos de reflexão de dois sismos; na área central foram encontradas espessuras da crosta entre 40 e 63 km para três pontos de reflexão de um sismo; finalmente, na região sul, foram encontradas espessuras crustais entre 35 e 40 km para 7 pontos de reflexão de 3 sismos. Foi obtido a relação Vp/Vs = 1.79 ± 0.16 para a região norte dos Andes, usando as diferenças de tempo sS-smS e pP-pmP. Os resultados obtidos vão aumentar o banco de dados de espessura crustal nos Andes setentrionais e na parte norte dos Andes centrais, como também possibilitaram uma melhora do modelo de espessura crustal desenvolvido por Assumpção et al. (2013) para América do Sul, permitindo assim a analise e comparação entre os dados obtidos nesse estudo e o modelo crustal de Assumpção et al. (2013). / In this work we studied the crustal thickness in northern Andes and northern central Andes (between 10ºN and 9ºS). We analyzed the reflections from the underside of the Moho for six intermediate earthquakes occurred in the study area to estimate the crustal thickness at the bounce point. The pmP and smS phases (reflections at the Moho) are identied as precursors to the depth phase pP and sS, respectively, when recorded at teleseismic distances. This method require events of magnitude > 6 and depth > 100 km. In order to better identify those two reflections, it was necessary to stacking traces from different stations. The study area was divided in three sections: northern area between 6ºN - 0º (Colombia), in this area we analyzed two earthquakes and the crustal thicknesses in four bounce points were found; central area between 2ºS- 5ºS (Ecuador), in this area we studied one earthquake and the crustal thicknesses in three bounce points were found; southern area between 6ºS - 9ºS (Peru-Brazil border), where we analyzed three earthquakes and the crustal thicknesses in seven bounce points were found. Crustal thickness found in the southern area varied from 35 to 40 km, from 40 to 63 km in the central area and from 26 to 56 km in the northern area. We estimated a Vp/Vs ratio of 1.79±0.16 for the northern Andes using the time differences sS-smS and pP-pmP relation. Our analysis complements the data base of Moho depths in the Northern Andes and also, it allows a comparison with the model of crustal thickness estimated by Assumpção et al. (2013) confirming the Moho depth on both studies.

A dissertation submitted to the Faculty of Science, University of the Witwatersrand in fulfilment of the requirements for the degree of Master of Science, School of Geosciences. Johannesburg, 2016. / The measurements of physical rock properties, seismic velocities in particular, associated with ore deposits and their host rocks are crucial in interpreting seismic data collected at the surface for mineral exploration purposes. The understanding of the seismic velocities and densities of rock units can help to improve the understanding of seismic reflections and thus lead to accurate interpretations of the subsurface geology and structures. This study aims to determine the basic acoustic properties and to better understand the nature of the seismic reflectivity of the world’s deepest gold-bearing reef, the Carbon Leader Reef (CLR). This was done by measuring the physical properties (ultrasonic velocities and bulk densities) as well as conducting mineralogical analyses on drill-core samples. Ultrasonic measurements of P- and S-wave velocities were determined at ambient and elevated stresses, up to 65 MPa. The results show that the quartzite samples overlying and underlying the CLR exhibit similar velocities (~ 5028 m/s-5480 m/s and ~ 4777 m/s-5211 m/s, respectively) and bulk densities (~ 2.68 g/cm3 and 2.66 g/cm3). This is due to similar mineralogy and chemical compositions observed within the units. However, the CLR has slightly higher velocity (~ 5070 m/s-5468 m/s) and bulk density (~ 2.78 g/cm3) than the surrounding quartzite units probably due to higher pyrite content in the reef, which increases the velocity. The hangingwall Green Bar shale exhibits higher velocity (5124 m/s-5914 m/s) and density values (~ 2.89 g/cm3-3.15 g/cm3) compared to all the quartzite units (including the CLR), as a result of its finer grain size and higher iron and magnesium content. In the data set it is found that seismic velocities are influence by silica, iron and pyrite content as well as the grain size of the samples, i.e., seismic velocities increase with (1) decreasing silica content, (2) increasing iron and pyrite content and (3) decreasing grain size. Reflection coefficients calculated using the seismic velocities and densities at the boundaries between the CLR and its hangingwall and footwall units range between ~0.02 and 0.05, which is below the suggested minimum of 0.06 required to produce a strong reflection between two lithological units. This suggests that reflection seismic methods might not be able to directly image the CLR as a prominent reflector, as observed from the seismic data. The influence of micro-cracks is observed in the unconfined uniaxial compressive stress tests where two regimes can be identified: (1) From 0 - 25 MPa the P-wave velocities increase with progressive loading, but at different rates in shale and quartzite rocks owing to the presence of micro-cracks and (2) above stresses of ~20 - 25 MPa, the velocity stress relationship becomes constant, possibly indicating total closure of micro-cracks. The second part of the study integrates 3D reflection seismic data, seismic attributes and information from borehole logs and underground mapping to better image and model important fault systems that might have a direct effect on mining in the West Wits Line goldfields. 3D seismic data have delineated first-, second- and third-order scale faults that crosscut key gold-bearing horizons by tens to hundreds of metres. Applying the modified seismic attribute has improved the imaging of the CLR by sharpening the seismic traces. Conventional interpretation of the seismic data shows that faults with throws greater than 25 m can be clearly seen. Faults with throws less than 25 m were identified through volumetric (edge enhancement and ant-tracking seismic attributes) and horizon-based (dip, dip-azimuth and edge detection seismic attributes) seismic attribute analysis. These attributes provided more accurate mapping of the depths, dip and strikes of the key seismic horizon (Roodepoort shale), yielding a better understanding of the relationship between fault activity, methane migration and relative chronology of tectonic events in the goldfield. The strato-structural model derived for the West Wits Line gold mines can be used to guide future mine planning and designs to (1) reduce the risks posed by mining activities and (2) improve the resource evaluation of the goldbearing reefs in the West Wits Line goldfields. / LG2017

The Pile Mound survey includes magnetometry paired with targeted ground-penetrating radar (GPR) and electromagnetic induction (EMI) surveys of the mound and testing of associated features over the ca. 6.5 ha site. The GPR survey discovered six rock features (five large rock features within the mound and one marking the outside of the mound). Knowledge of mounds in the Upper Cumberland Plateau (UCP) is lacking—the closest other studied sites are at the Corbin Site, Croley-Evans, Bell Site, and Beasley Mounds, approximately 75 – 100 km away. However, the most similar mound construction is found at Corbin and Cherokee sites, some 175 – 275 km away. In addition, the associated ceramic assemblage appears to reflect more similarity to the East Tennessee Valley rather than the Middle Cumberland region. These data provide a unique opportunity to better understand the Mississippian occupation in the UCP of Tennessee.

South Sudan, a new country formed in 2011, has been planning to develop its mineral sector by allocating exploration licenses to investors. This decision requires preliminary knowledge of geology and mineral occurrences, both of which are unavailable because the country has been engaged in a civil war for more than 50 years. Exploration of mineral resources in South Sudan has lagged behind its petroleum industry, except for artisanal gold mining, which is practiced intermittently by local communities. Freely available satellite gravity and remote-sensing data were used to map the basement architecture as well as zones of hydrothermal alteration in the Didinga Hills; both basement architecture and hydrothermal alteration are of prime importance in exploration and development of mineral resources in the study area. Qualitative interpretation of gravity data is consistent with the known geology of petroleum fields and the Precambrian basement complex. Remote-sensing data and techniques—optimal band combination, band ratioing, and principal component analysis—have been effective in extracting information related to lithology, hydrothermal alteration, and geologic structures. The resulting basic information and methods have identified additional prospective exploration areas where more detailed gravity, magnetic, electromagnetic, and seismic surveys should be carried out; this will assist decision makers in matters related to land use, mineral titles, and exploration of natural resources, and lead to prosperity for the new nation of South Sudan.

The Charleston Uplift (CU), a 30-km-long by 7-km-wide, N46°E-oriented subsurface geologic anomaly in the northern Mississippi embayment near Charleston, Missouri, exhibits up to 36 m of vertical relief across the Paleogene/Quaternary unconformity. Subsurface structural relief, along with the CU’s coincident boundary alignment with contemporary microseismicity and the New Madrid North Fault (NMNF), suggest a structural origin. Subsequent seismic soundings indicate vertical structural relief is present in Cretaceous and Paleozoic horizons, supporting the fault-controlled origin. The southern boundary (CU-s) had not been investigated, nor had any direct fault images been acquired. Integrated microgravity and seismic-reflection methods across the inferred CU-s establish the first image of this fault. Forward modeling indicated that the vertical variation of strata across the CU-s would induce a microgravity anomaly of 1.6 mGal. The observed microgravity anomaly survey across the southern boundary is 1.616 ± .004 mGal, and is consistent with the tectonic interpretation. A subsequently acquired seismic-reflection profile corroborates this interpretation. The imaged fault shows approximately 60, 35, and 35 meters of vertical down-to-the-south throw across the tops of Paleozoic, Cretaceous, and Tertiary horizons, respectively. This confirms the CU is not an erosional feature, but a structurally controlled extension of the NMNF.

The San Andreas Fault stretches for over 800 miles through California. Along the foothills of the San Bernardino Mountains, areas in close proximity to the San Andreas Fault Zone may be subject to site amplification of ground motion caused by seismic activity via wave propagation through the subsurface. These seismic hazards are being addressed via the Alquist-Priolo Earthquake Faulting Zone Act and the National Earthquake Hazards Reduction Program (NEHRP). Shear wave velocity of the subsurface has served as a proxy for ground motion amplification and is therefore a useful parameter to help analyze and reduce seismic hazards. Low shear wave velocities of the subsurface have been known to correlate with higher amplitude ground motion. This study focuses on refraction microtremor analysis (ReMi) of the subsurface in Northern San Bernardino; more specifically, areas encompassing California State University San Bernardino, in close proximity to the San Andreas Fault. The technique will resolve shear wave velocity values for the top 30 meters (Vs30) of the subsurface. This depth of investigation has proven to be an effective means in determining subsurface conditions. ReMi profiles were situated 0.25 to 2.0 miles away from the San Andreas Fault, and in some instances, strategically positioned next to housing developments and structures. Phase velocity dispersion curves were generated by processing ReMi seismic data and subsequently inverted to attain average shear wave velocity profiles with depth. The geologic units in the study area consist of very young wash deposits, young alluvial fan deposits and Pelonist schist deposits. These geologic units may be an indicator to how seismic waves behave in subsurface lithology. To highlight differences in Vs30 values across the project area, a microzonation map was constructed.

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.