This study focuses on mapping faults along the Creeping Section of the San Andreas Fault (CSAF) in California between San Juan Bautista (121.54°W 36.85°N) and Parkfield (120.41°W 35.87°N). I synthesize high-quality base data, including and lidar topography from B4, EarthScope, and USGS 3DEP, recent maps of decadal-scale along-fault shear strain, and aerial and satellite imagery. Using these data, I produced (covering 150 km at 1:10,000 scale) three geospatial map datasets with attributes: geomorphic indicators of faulting, surficial geology, and active fault traces.The CSAF's creeping movement, though likely not associated with large earthquakes, has the potential to cause damage to infrastructure. Accurate fault mapping facilitates fault displacement hazard assessment. This type of work is useful for California state regulations, particularly the Alquist-Priolo Act of 1972, providing insights for engineering site assessments and fault exclusion zones. I discern, categorize, and rank geomorphic indicators to support fault line placement. This approach contributes to the identification of surface expression of creeping faults where the surface has undergone alteration in response to displacement along the fault. I created a surficial geologic map spanning from San Juan Bautista to the southern extent of EarthScope lidar coverage (120.59°W 36.03°N). I categorized each fault as either a primary or secondary fault trace and further broke them into confidence levels based on interpretations of indicators along with structural geologic reasoning and topographic patterns. Accessible target areas containing initial low confidence mapping or interesting structures were visited in the field. Zones along the creeping section exhibit structures such as a pressure ridge found 25 km north of Parkfield, sigmoidal faults and sagponds observed near Paicines Ranch (121.29°W 36.68°N), en-echelon faults, horsetail splays and Riedel shear structures near Lewis Creek (120.87°W 36.29°N). Controls on the structural style along the CSAF are the results of geologic units through which the faults cut and fault zone width and trend.

Accurate fault maps are an important component in the assessment of hazard from fault displacement. Different mapping techniques, biases and ambiguous geomorphic evidence for faulting can drive even expert mappers to produce different fault maps. Another challenge is that future ruptures may not follow past ruptures, so available evidence in the landscape may not lead to accurate rupture prediction. The ultimate goal of my work is to develop a systematized approach for fault mapping so that resulting maps are more evidence-based and ultimately of higher quality I systematized the active fault mapping process and the documentation of evidence for potential fault rupture. I developed and taught a systematic mapping process based on geomorphic landforms evident in remote sensing datasets to undergraduate students, graduate students, and geologic professionals. My approach uses data acquired before historic ruptures to make and test “pre-rupture” fault traces based on the landscape morphology, geomorphology, and geology. The mappers used the Geomorphic Indicator Ranking system (GIR) to represent the geomorphic evidence for faulting such as scarps, triangular facets, offset features, beheaded drainages, and many more. I evaluated the approach in three ways: (1) To assess the geomorphology that best predicts future rupture, I compared the separation distance between the mapped geomorphologic features and the rupture. Scarps and lineaments performed best. (2) I compared the fault confidence chosen by the mapper versus that computed from GIR elements (i.e., mapped geomorphology) near the fault traces. Accurately characterizing fault confidence requires a balance between the mapper input and the calculated confidence rankings. (3) I conducted listening sessions with 21 participants to understand each participant’s approach to fault mapping to highlight best practices and challenges of geomorphic fault mapping. The terminology and mapping process vary by experience level. My approach works both as a teaching tool to introduce tectonic geomorphology and fault mapping to novice mappers, but also works in an industry setting to establish consistent documentation for fault maps. These higher quality fault maps have implications applications of fault mapping including easier dissemination of information, comparison between different fault maps, and hopefully more accurate fault locations for hazard mitigation.

This study presents an analysis of fault scarps, with a focus on implementing the Landlab computational toolkit to model fault scarp evolution and analyzing fault scarps under transport and production-limited conditions with linear and nonlinear diffusive transport laws. The aim of the study is to expand diffusion modeling of fault scarps from 1D to 2D by using Landlab toolkit. The study evaluated two fault scarps in western US (NE California): one representing an old fault scarp (Twin Butte) and the other representing a young fault scarp (Active Hat Creek Fault). High-resolution digital elevation models (DEMs) were used to generate 2D surfaces of the fault scarps, which were then converted to 1D profiles for morphological modeling and analysis. The accuracy of the models was evaluated using Root Mean Squared Error (RMSE), and the best-fit models were selected for further examination. The grid search of the non-linear diffusion model of the Twin Butte and Active Hat Creek fault scarps showed optimum values for transport constant (k) and scarp age (t) that aligned with the apparent ages of the rocks and associated fault scarps. For both fault scarps, the optimum k value was around 7.5 m2 /kyr, while the optimum t value was around 110 kyr for the Twin Butte scarp and around 26 kyr for the Active Hat Creek scarp. The results suggest that the geomorphic processes (influenced by climate and rock types) in both fault scarps are similar, despite the difference in age and location. Integrating tectonic displacement in the model helps to better capture the observed patterns of tectonic deformation. The expansion of the fault scarps diffusion model from 1D to 2D opens up a range of fascinating possibilities, as it enables us to model the lateral movement of particles that the 1D model typically overlooks. By incorporating this additional dimension, we can better understand the complex interplay between vertical and horizontal displacements, providing a more accurate representation of the geological processes at work. This advancement ultimately allows for a more comprehensive analysis of fault scarps and their development over time, enhancing our understanding of Earth's dynamic crustal movements.

Despite the rapid adoption of robotics and machine learning in industry, their application to scientific studies remains under-explored. Combining industry-driven advances with scientific exploration provides new perspectives and a greater understanding of the planet and its environmental processes. Focusing on rock detection, mapping, and dynamics analysis, I present technical approaches and scientific results of developing robotics and machine learning technologies for geomorphology and seismic hazard analysis. I demonstrate an interdisciplinary research direction to push the frontiers of both robotics and geosciences, with potential translational contributions to commercial applications for hazard monitoring and prospecting. To understand the effects of rocky fault scarp development on rock trait distributions, I present a data-processing pipeline that utilizes unpiloted aerial vehicles (UAVs) and deep learning to segment densely distributed rocks in several orders of magnitude. Quantification and correlation analysis of rock trait distributions demonstrate a statistical approach for geomorphology studies. Fragile geological features such as precariously balanced rocks (PBRs) provide upper-bound ground motion constraints for hazard analysis. I develop an offboard method and onboard method as complementary to each other for PBR searching and mapping. Using deep learning, the offboard method segments PBRs in point clouds reconstructed from UAV surveys. The onboard method equips a UAV with edge-computing devices and stereo cameras, enabling onboard machine learning for real-time PBR search, detection, and mapping during surveillance. The offboard method provides an efficient solution to find PBR candidates in existing point clouds, which is useful for field reconnaissance. The onboard method emphasizes mapping individual PBRs for their complete visible surface features, such as basal contacts with pedestals–critical geometry to analyze fragility. After PBRs are mapped, I investigate PBR dynamics by building a virtual shake robot (VSR) that simulates ground motions to test PBR overturning. The VSR demonstrates that ground motion directions and niches are important factors determining PBR fragility, which were rarely considered in previous studies. The VSR also enables PBR large-displacement studies by tracking a toppled-PBR trajectory, presenting novel methods of rockfall hazard zoning. I build a real mini shake robot providing a reverse method to validate simulation experiments in the VSR.

Rivers in steep mountainous landscapes control how, where, and when signals of base-level fall are transmitted to the surrounding topography. In doing so rivers play an important role in determining landscape evolution in response to external controls of tectonics and climate. However, tectonics and climate often covary and understanding how they influence landscape evolution remains a significant challenge. The Hawaiian Islands, where tectonics are minimized but climate signals are amplified, provide an opportunity to better understand how signals of climate are recorded by landscapes. Focusing on the Hawaiian Islands, I examine (1) how variability in rock mass properties and thresholds in sediment mobility determine where waterfalls form or stall along the Nāpali coast of Kauaʻi, (2) I then extend these findings to other volcanoes to test if observed physical limits in flood size, climate, and volcano gradient can determine where waterfalls form, and (3) I explore how thresholds in river incision below waterfalls limit information about the influence of climate on river incision rates. Findings from this analysis show that waterfalls form or stall where the maximum unit stream power is at or below a critical unit stream power for bedrock river incision. Climate appears to have little effect in determining where these conditions are met but where waterfalls stall or form does record information about discharge-area scaling for global maximum observed floods. Below waterfalls the maximum incision depth for rivers on the island of Kauaʻi (which formed ~ 4-5 million years ago) is approximately proportional to the inverse square root of mean annual rainfall. Though maximum river incision depths for some of the younger volcanoes do not exhibit the same dependency on mean annual rainfall rates they are comparable to the maximum incision depths observed on Kauaʻi even though they are a quarter to one-tenth the age of Kauaʻi. Importantly, these patterns of incision can be explained by thresholds in sediment mobility as recorded by river longitudinal profiles and indicate that the Hawaiian Islands are dominated by threshold conditions where signals of climate are recorded in the topography through controls on incision depth but not incision rates.

Mountain landscapes reflect competition between tectonic processes acting to build topography and erosive processes acting to wear it down. In temperate mountain landscapes, bedrock rivers are the primary erosional agent, setting both the pace of landscape evolution and form of the surrounding topography. Theory predicts that river steepness is sensitive to climatic, tectonic, and lithologic factors, which dictate the rates and mechanics of erosional processes. Thus, encoded into topography is an archive of information about forces driving landscape evolution. Decoding this archive, however, is fraught and climate presents a particularly challenging conundrum: despite decades of research describing theoretically how climate should affect topography, unambiguous natural examples from tectonically active landscapes where variations in climate demonstrably influence topography are elusive. In this dissertation, I first present a theoretical framework describing how the spatially varied nature of orographic rainfall patterns, which are ubiquitous features of mountain climates, complicate expectations about how climate should influence river steepness and erosion. I then apply some of these ideas to the northern-central Andes. By analyzing river profiles spanning more than 1500 km across Peru and Bolivia, I show that the regional orographic rainfall pattern this landscape experiences systematically influences fluvial erosional efficiency and thus topography. I also show how common simplifying assumptions built into conventional topographic analysis techniques may introduce biases that undermine detection of climatic signatures in landscapes where climate, tectonics, and lithology all covary – a common condition in mountain landscapes where these techniques are often used. I continue by coupling this analysis with published erosion rates and a new dataset of 25 cosmogenic 10Be catchment average erosion rates. Once the influence of climate is accounted for, functional relationships emerge among channel steepness, erosion rate, and lithology. I then use these functional relationships to produce a calibrated erosion rate map that spans over 300 km of the southern Peruvian Andes. These results demonstrate that accounting for the effects of climate significantly enhances the ability to decode channel steepness patterns. Along with this comes the potential to better understand rates and patterns of tectonic processes, and identify seismic hazards associated with tectonic activity using topography.

Rock traits (grain size, shape, orientation) are fundamental indicators of geologic processes including geomorphology and active tectonics. Fault zone evolution, fault slip rates, and earthquake timing are informed by examinations of discontinuities in the displacements of the Earth surface at fault scarps. Fault scarps indicate the structure of fault zones fans, relay ramps, and double faults, as well as the surface process response to the deformation and can thus indicate the activity of the fault zone and its potential hazard. “Rocky” fault scarps are unusual because they share characteristics of bedrock and alluvial fault scarps. The Volcanic Tablelands in Bishop, CA offer a natural laboratory with an array of rocky fault scarps. Machine learning mask-Region Convolutional Neural Network segments an orthophoto to identify individual particles along a specific rocky fault scarp. The resulting rock traits for thousands of particles along the scarp are used to develop conceptual models for rocky scarp geomorphology and evolution. In addition to rocky scarp classification, these tools may be useful in many sedimentary and volcanological applications for particle mapping and characterization.

Faults found in the arid to semi-arid Basin and Range Physiographic province of the southwestern US are given broad age definitions in terms of which features appear to be the oldest. Particularly in the northwestern corner of Arizona, detailed geomorphic studies on the tectonic history and timing of faulting are not widespread. At the base of the Virgin Mountains in northwestern Arizona is a fault scarp along the Piedmont Fault line. This normal fault crosses a series of alluvial fans that are filled with sediments of ambiguous ages. Previous studies that were done in this region find a broad, Miocene age for the exhumation and uplift of these surfaces, with some indications of Laramide faulting history. However, specific fault characteristics and a time constraint of the tectonic history of the Piedmont Fault scarp has yet to be established. Here, we aim to determine the age, fault-slip rate, seismic history, and potential hazard of the fault scarp near Scenic and Littlefield, Arizona through structure from motion (SfM) modeling, which is a form of photogrammetry using a drone. In addition, we distinguish the climatic and tectonic influences on the geomorphology observed along the scarp through analysis along the fault line. With data collected from a ~500 m section of the fault, we present results from a digital elevation model (DEM) and orthophotos derived through the SfM modelling. Based on field observations and morphologic dating, we determine that the Piedmont Fault experiences an approximately continuous fault-slip and an earthquake recurrence interval in the range of 7,000 years. The approximate age of the scarp is 16.0 ka ± 5 kyr. Therefore, we conclude that the earthquake hazard posed to nearby cities is minimal but not nonexistent. Future work includes further analysis of fault profiles due to uncertainty in the present one and Terrestrial Cosmogenic Nuclide (TCN) dating of samples taken from the tops of boulders in a residual debris flow sitting on faulted and unfaulted alluvia. Determining the ages for these boulder surfaces can hopefully further inform our knowledge of the tectonic activity present in the North Virgin Mountains.

Cinder cones are common volcanic structures that occur in fields, and on the flanks of shield volcanoes, stratovolcanoes, and calderas. Because they are common structures, they have a significant possibility of impacting humans and human environments. As such, there is a need to analyze cinder cones to get a better understanding of their eruptions and associated hazards. I will approach this analysis by focusing on volcanic bombs and ballistics, which are large clots of lava that are launched from the volcanic vent, follow ballistic trajectories, and can travel meters to a few kilometers from their source (e.g. Fagents and Wilson 1993; Waitt et al. 1995).
Tecolote Volcano in the Pinacate Volcanic Field in Mexico contains multiple vents within a horseshoe-shaped crater that have all produced various ejecta (Zawacki et al. 2019). The objectives of this research are to map ballistic distribution to understand the relationship between the source vent or vents and the bombs and ballistics that litter the region around Tecolote, and interpret the eruption conditions that ejected those bombs by using their distributions, morphologies, and fine-scale textures.
The findings of this work are that these bombs are apparently from the last stages of the eruption, succeeding the final lava flows. The interiors and exteriors of the bombs display different cooling rates which can are indicated by the fabric found within. Using this, certain characteristics of the bombs during eruption were extrapolated. The ‘cow pie’ bombs were determined to be the least viscous or contained a higher gas content at the time of eruption. Whereas the ribbon/rope bombs were determined to be the most viscous or contained a lesser gas content. Looking at the Southern Bomb Field site, it is dominated by large bombs that were during flight were molded into aerodynamic shapes. The Eastern Rim site is dominated by smaller bombs that appeared to be more liquid during the eruption. This difference in the two sites is a probable indication of at least two different eruptive events of different degrees of explosivity. Overall, aerodynamic bombs are more common and extend to greater distances from the presumed vent (up to 800 m), while very fluidal bombs are uncommon beyond 500 meters. Fluidal bombs (‘cow pie’, ‘ribbon’, ‘rope/spindle’) show a clear trend in decreasing size with distance from vent, whereas the size-distance trend is less dramatic for the aerodynamic bombs.

East African extensional basins have played a crucial role in revealing the evolution and characteristics of the early stages of continental rifting and for providing the geological context of hominin evolution and innovation. The numerous volcanic eruptions, rapid sedimentation and burial, and subsequent exposure through faulting and erosion, provide excellent conditions for the preservation of tectonic history, paleoenvironment data, and vertebrate fossils. The reconstruction of depositional environments and provision of geochronologic frameworks for hominin sites have been largely provided by geologic investigations in conjunction with paleontological studies, like the Ledi-Geraru Research Project (LGRP). High-resolution paleoclimate records that can be directly linked to hominin fossil outcrops have been developed by the Hominin Sites and Paleolakes Drilling Project (HSPDP) which collected sedimentary-paleolake cores at or near key hominin fossil sites.

Two chapters of this dissertation are a result of research associated with the HSPDP. For HSPDP, I establish a tephrostratigraphic framework for the drill cores from the Northern Awash (Afar, Ethiopia) and Baringo-Tugen Hills-Barsemoi (Kenya) HSPDP sites. I characterize and fingerprint tephra through glass shard and feldspar phenocryst geochemistry. From tephra geochemical analyses, I establish chronostratigraphic ties between the HSPDP cores’ high-resolution paleoclimate records to outcrop stratigraphy which are associated with hominin fossils sites.

Three chapters of this dissertation are a result of field work with the LGRP. I report new geological investigations (stratigraphic, tectonic, and volcanic) of two previously unmapped regions from the eastern Ledi-Geraru (ELG), Asboli and Markaytoli. Building upon this research I present interpretations from tephra analyses, detailed stratigraphic analyses, and geologic mapping, of the Pleistocene (~2.6 to < 2.45 Ma) basin history for the LGRP. My work with the LGRP helps to reconstruct a more complete Early Pleistocene depositional and geologic history of the lower Awash Valley.

Overall, this dissertation contributes to the reconstruction of hominin paleoenvironments and the geochronological framework of the Pliocene and Pleistocene faunal/hominin records. It further contributes to rift basin history in East Africa by elaborating the later structural and stratigraphic history of the lower Awash region.