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.

The tectonism, volcanism, and sedimentation along the East African Rift System (EARS) produced a series of rift basins with a rich paleoanthropological record, including a Late Miocene–present record of hominin evolution. To better understand the relationship between Earth system history and human evolution within the EARS, the Hominin Sites and Paleolakes Drilling Project (HSPDP) collected paleolake sediments near key paleoanthropological sites in Ethiopia and Kenya, compiling a multi-proxy, high-resolution geological and environmental record. As part of the HSPDP, I studied the detrital mineral record of the basins and evaluated tectonic and climatic controls on East African landscapes during the Plio-Pleistocene using samples from three of the drill sites, Chew Bahir: (CHB, ~620–present; Ethiopia), Northern Awash (NA, ~3.3–2.9 Ma; Ethiopia,), and West Turkana (WTK, ~1.9–1.4 Ma; Kenya). I employed laser ablation U/Pb and (U-Th)/He double dating (LADD) of detrital zircons, which yields paired U/Pb and (U-Th)/He dates, and (U-Th)/He dating of detrital apatites to evaluate sediment provenance and the cooling history of the source rocks. In addition, I used in situ 10Be cosmogenic radionuclide analyses to determine paleoerosion rates. Two chapters of this dissertation focus on results from the NA and WTK drill sites. Source units for the NA and WTK drill sites are largely Cenozoic volcanic rocks, and the detrital zircon record yields an extensive record of the timing of various phases of volcanism within the EARS. Exceptionally young zircon (U-Th)/He dates reflect partial resetting associated with late mafic volcanism and/or hydrothermal activity. Erosion rates are consistent and relatively low across the Plio-Pleistocene, despite significant tectonic and geomorphic shifts in the landscape. Two other chapters of this dissertation cover results from the CHB drill site. The Chew Bahir basin has significant exposures of Neoproterozoic and Early Paleozoic crystalline basement units, and the detrital zircon record yields one singular phase of volcanism in the EARS. The CHB erosion rates show an overall decreasing trend over time, consistent with an aridifying climate, and increased environmental variability after ~200 ka.

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.

The mountains of western North America are spectacular and diverse, from sheer walls of crumbling black limestone in the Canadian Rockies, to smooth glacially polished granite in the Wind River Range, to gargantuan ice-clad volcanoes in the Cascades. These great bastions of rock, snow, and ice, still very much wild and untamed, provide an incredible arena for adventure, exploration, and challenge. Over the past three years, I have devoted thousands of hours to exploring these vast wild places, climbing high peaks, steep cliffs, and frozen waterfalls. In doing so, I studied the rich geologic history of the mountains. This thesis project is a compilation of stories and images from those adventures, along with the stories of the mountains themselves: how the rocks were formed, thrust skyward, and sculpted over the ages into their present, glorious form. The photographic and detailed narrative of the geology and adventures is on a new website called Cloud Piercers, which currently features three geologically diverse mountain massifs: (1) Mount Rainier, an active volcano in the Cascade Range of Washington; (2) Mount Robson, the highest peak in the Canadian Rockies, within a terrain of folded Paleozoic sedimentary rocks; and (3) the Wind River Range of Wyoming, composed mostly of Archean metamorphic and granitic rocks. This website will be expanded in the future as the geologic studies and adventures continue.

Worldwide, rivers and streams make up dense, interconnected conveyor belts of sediment– removing carved away earth and transporting it downstream. The propensity of alluvial river beds to self-organize into complex trains of bedforms (i.e. ripples and dunes) suggests that the associated fluid and sediment dynamics over individual bedforms are an integral component of bedload transport (sediment rolled or bounced along the river bed) over larger scales. Generally speaking, asymmetric bedforms (such as alluvial ripples and dunes) migrate downstream via erosion on the stoss side of the bedform and deposition on the lee side of the bedform. Thus, the migration of bedforms is intrinsically linked to the downstream flux of bedload sediment. Accurate quantification of bedload transport is important for the management of waters, civil engineering, and river restoration efforts. Although important, accurate qualification of bedload transport is a difficult task that continues t elude researchers. This dissertation focuses on improving our understanding and quantification of bedload transport on the two spatial scales: the bedform scale and the reach (~100m) scale.

Despite a breadth of work investigating the spatiotemporal details of fluid dynamics over bedforms and bedload transport dynamics over flat beds, there remains a relative dearth of investigations into the spatiotemporal details of bedload transport over bedforms and on a sub-bedform scale. To address this, we conducted two sets of flume experiments focused on the two fundamental regions of flow associated with bedforms: flow separation/reattachment on the lee side of the bedform (Chapter 1; backward facing-step) and flow reacceleration up the stoss side of the next bedform (Chapter 2; two-dimensional bedform). Using Laser and Acoustic Doppler Velocimetry to record fluid turbulent events and manual particle tracking of high-speed imagery to record bedload transport dynamics, we identified the existence and importance of “permeable splat events” in the region proximal to flow reattachment.

These coupled turbulent and sediment transport events are integral to the spatiotemporal pattern of bedload transport over bedforms. Splat events are localized, high magnitude, intermittent flow features in which fluid impinges on the bed, infiltrates the top portion of bed, and then exfiltrates in all directions surrounding the point of impingement. This initiates bedload transport in a radial pattern. These turbulent structures are primarily associated with quadrant 1 and 4 turbulent structures (i.e. instantaneous fluid fluctuations in the streamwise direction that bring fluid down into the bed in the case of quadrant 1 events, or up away from the bed in the case of quadrant 4 events) and generate a distinct pattern of bedload transport compared to transport dynamics distal to flow reattachment. Distal to flow reattachment, bedload transport is characterized by relatively unidirectional transport. The dynamics of splat events, specifically their potential for inducing significant magnitudes of cross-stream transport, has important implications for the evolution of bedforms from simple, two dimensional features to complex, three-dimensional features.

New advancements in sonar technology have enabled more detailed quantification of bedload transport on the reach scale, a process paramount to the effective management of rivers with sand or gravel-dominated bed material. However, a practical and scalable field methodology for reliably estimating bedload remains elusive. A popular approach involves calculating transport from the geometry and celerity of migrating bedforms, extracted from time-series of bed elevation profiles (BEPs) acquired using echosounders. Using two sets of repeat multibeam sonar surveys from the Diamond Creek USGS gage station in Grand Canyon National Park with large spatio-temporal resolution and coverage, we compute bedload using three field techniques for acquiring BEPs: repeat multi-, single-, and multiple single-beam sonar. Significant differences in flux arise between repeat multibeam and single beam sonar. Mulitbeam and multiple single beam sonar systems can potentially yield comparable results, but the latter relies on knowledge of bedform geometries and flow that collectively inform optimal beam spacing and sampling rate. These results serve to guide design of optimal sampling, and for comparing transport estimates from different sonar configurations.

Located in the Sunbelt of the Southwestern United States, Phoenix Arizona finds itself in one of the hottest, driest places in the world. Thankfully, Phoenix has the Salt River, Gila River, Verde River, and a vast aquifer to meet the water demands of the municipal, industrial, and agricultural sectors. However, rampant groundwater pumping and over-allocation of these water supplies based on unprecedented, high flows of the Colorado River have created challenges for water managers to ensure adequate water supply for the future. Combined with the current 17-year drought and the warming and drying projections of climate change, the future of water availability in Phoenix will depend on the strength of water management laws, educating the public, developing a strong sense of community, and using development to manage population and support sustainability. As the prevalence of agriculture declines in and around Phoenix, a substantial amount of water is saved. Instead of storing this saved water, Phoenix is using it to support further development. Despite uncertainty regarding the abundant and continuous availability of Phoenix's water resources, development has hardly slowed and barely shifted directions to support sustainability. Phoenix was made to grow until it legally cannot expand anymore. In order to develop solutions, we must first understand the push for development in water-stressed Phoenix, Arizona.

There is a need to understand spatio-temporal variation of slip in active fault zones, both for the advancement of physics-based earthquake simulation and for improved probabilistic seismic hazard assessments. One challenge in the study of seismic hazards is producing a viable earthquake rupture forecast—a model that specifies the expected frequency and magnitude of events for a fault system. Time-independent earthquake forecasts can produce a mismatch among observed earthquake recurrence intervals, slip-per-event estimates, and implied slip rates. In this thesis, I developed an approach to refine several key geologic inputs to rupture forecasts by focusing on the San Andreas Fault in the Carrizo Plain, California. I use topographic forms, sub-surface excavations, and high-precision geochronology to understand the generation and preservation of slip markers at several spatial and temporal scales—from offset in a single earthquake to offset accumulated over thousands of years. This work results in a comparison of slip rate estimates in the Carrizo Plain for the last ~15 kyr that reduces ambiguity and enriches rupture forecast parameters. I analyzed a catalog of slip measurements and surveyed earth scientists with varying amounts of experience to validate high-resolution topography as a supplement to field-based active fault studies. The investigation revealed that (for both field and remote studies) epistemic uncertainties associated with measuring offset landforms can present greater limitations than the aleatoric limitations of the measurement process itself. I pursued the age and origin of small-scale fault-offset fluvial features at Van Matre Ranch, where topographic depressions were previously interpreted as single-event tectonic offsets. I provide new estimates of slip in the most recent earthquake, refine the centennial-scale fault slip rate, and formulate a new understanding of the formation of small-scale fault-offset fluvial channels from small catchments (<7,000 m2). At Phelan Creeks, I confirm the constancy of strain release for the ~15,000 years in the Carrizo Plain by reconstructing a multistage offset landform evolutionary history. I update and explicate a simplified model to interpret the geomorphic response of stream channels to strike-slip faulting. Lastly, I re-excavate and re-interpret paleoseismic catalogs along an intra-continental strike-slip fault (Altyn Tagh, China) to assess consistency of earthquake recurrence.

For this dissertation, three separate papers explore the study areas of the western Grand Canyon, the Grand Staircase (as related to Grand Canyon) and Desolation Canyon on the Green River in Utah.

In western Grand Canyon, I use comparative geomorphology between the Grand Canyon and the Grand Wash Cliffs (GWC). We propose the onset of erosion of the GWC is caused by slip on the Grand Wash Fault that formed between 18 and 12 million years ago. Hillslope angle and channel steepness are higher in Grand Canyon than along the Grand Wash Cliffs despite similar rock types, climate and base level fall magnitude. These experimental controls allow inference that the Grand Canyon is younger and eroding at a faster rate than the Grand Wash Cliffs.

The Grand Staircase is the headwaters of some of the streams that flow into Grand Canyon. A space-for-time substitution of erosion rates, supported by landscape simulations, implies that the Grand Canyon is the result of an increase in base level fall rate, with the older, slower base level fall rate preserved in the Grand Staircase. Our data and analyses also support a younger, ~6-million-year estimate of the age of Grand Canyon that is likely related to the integration of the Colorado River from the Colorado Plateau to the Basin and Range. Complicated cliff-band erosion and its effect on cosmogenic erosion rates are also explored, guiding interpretation of isotopic data in landscapes with stratigraphic variation in quartz and rock strength.

Several hypotheses for the erosion of Desolation Canyon are tested and refuted, leaving one plausible conclusion. I infer that the Uinta Basin north of Desolation Canyon is eroding slowly and that its form represents a slow, stable base level fall rate. Downstream of Desolation Canyon, the Colorado River is inferred to have established itself in the exhumed region of Canyonlands and to have incised to near modern depths prior to the integration of the Green River and the production of relief in Desolation Canyon. Analysis of incision and erosion rates in the region suggests integration is relatively recent.

The morphology of mountainous areas is strongly influenced by stream bed incision rates, but most studies of landscape evolution consider erosion at basin scales or larger. The research here attempts to understand the smaller-scale mechanics of erosion on exposed bedrock channels in the conceptual framework of an established saltation-abrasion model by Sklar and Dietrich [2004]. The recirculating flume used in this experiment allows independent control of bed slope, water discharge rate, sediment flux, and sediment grain size – all factors often bundled together in simple models of river incision and typically cross-correlated in natural settings. This study investigates the mechanics of erosion on exposed bedrock channels caused by abrasion of transported particles. Of particular interest are saltating particles, as well as sediment near the threshold between saltation and suspension - sediment vigorously transported but with significant interaction with the bed. The size of these erosive tools are varied over an order of magnitude in mean grain diameter, including a sand of D¬50 = 0.56 mm, and three gravel sizes of 3.39, 4.63, and 5.88 mm. Special consideration was taken to prevent any flow conditions that created a persistent alluvial cover. The erodible concrete substrate is fully exposed at all times during experiments reported here. Rates of erosion into the concrete substrate (a bedrock proxy) were measured by comparing topographic data before and after each experimental run, made possible by a precision laser mounted on a high speed computer-controlled cart. The experimental flume was able to produce flow discharge as high as 75 liters per second, sediment fluxes (of many varieties) up to 215 grams per second, and bed slopes up to 10%. I find a general positive correlation is found between erosion rate and bed slope, shear stress, grain size, and sediment flux.